列主元高斯消去法

一、题目：

解法：1、程序A=[2,1,-3,-1;3,1,0,7;-1,2,4,-2;1,0,-1,5]b=[1;2;-1;5]n=4 ;A1=[A,b];for i = 1:n-1[XX,r] = max(abs(A1(i:n,i)));r = r + i - 1;if r>itemp=A1(i,:);A1(i,:)=A1(r,:);A1(r,:)=temp;endif A1(i,i)==0 endfor p = i+1:nA1(p,:)=A1(p,:)-A1(i,:)*A1(p,i)/A1(i,i);endendA = A1(:,1:n); b = A1(:,n+1);x(n) = b(n)/A(n,n);for i = n-1:-1:1x(i)=b(i);for p=n:-1:i+1x(i) = x(i)-A(i,p)*x(p);endx(i)=x(i)/A(i,i);endx

运行结果：X=-3.38  3.71  -1.95294117647059  1.27

2、程序a=input('请输入线性方程组的系数矩阵a=');b=input('请输入b=');[m,n]=size(a);if m~=n    error('此矩阵非方阵，运行截止')endfor t=1:n-1    i=t;h=a(t,t);    for s=t+1:n        if abs(h)

运行结果：经3次消元后的增广矩阵为    3.0000    1.0000         0    7.0000    2.0000         0    2.3333    4.0000    0.3333   -0.3333         0   -0.0000   -3.5714   -5.7143   -0.2857         0    0.0000         0    3.4000    4.3200 线性方程组的解向量是x =    -3.3059    3.0235   -1.9529    1.2706

个人写作，仅供参考

如若喜欢，望查阅本人下篇