%求四阶线性方程组的MATLAB程序clearAb=0.001 2 1 5 1; 3 -4 0.1 -2 2; 2 -1 2 0.01 3; 1.1 6 2.3 9 4;%增广矩阵num=1 2 3 4;%未知量x的对应序号for i=1:3A=abs(Ab(i:4,i:4);%系数矩阵取绝对值r,c=find(A=max(A(:);r=r+i-1;%最大值对应行号c=c+i-1;%最大值对应列号q=Ab(r,:),Ab(r,:)=Ab(i,:),Ab(i,:)=q;%行变换w=Ab(:,c),Ab(:,c)=Ab(:,i),Ab(:,i)=w;%列变换n=num(i),num(i)=num(c),num(c)=n;%列变换引起未知量x次序变化 for j=i:3 Ab(j+1,:)=-Ab(j+1,i)*Ab(i,:)/Ab(i,i)+Ab(j+1,:);%消去过程 endend %最后得到系数矩阵为上三角矩阵%回代算法求解上三角形方程组x(4)=Ab(4,5)/
