Matlab程序Gauss列主元消去法.docxVIP

(完满版)Matlab程序Gauss列主元消去法 (完满版)Matlab程序Gauss列主元消去法 (完满版)Matlab程序Gauss列主元消去法 ．Gauss列主元消去法求解线性方程组 6 2 1 2 x1 4 2 5 0 2 x2 7 2 0 8 5 x3 , 1 1 3 2 7 x4 0 迭代法计算停止的条件为： max x(jk 1) x (j k ) 10 6． 1 j 3 采用 用用 Gauss列主元消去法时， Matlab计算程序为： clear clc A=[2 2 1 2;4 1 3 -1;-4 -2 0 1;2 3 2 3]; B=[1;2;1;0]; n=length(B); X=zeros(n,1); c=zeros(1,n); d1=0; for i=1:n-1 max=abs(A(i,i)); m=i; for j=i+1:n if maxabs(A(j,i)) max=abs(A(j,i)); m=j; end end if (m~=i) for k=i:n c(k)=A(i,k); A(i,k)=A(m,k); A(m,k)=c(k); end d1=B(i); B(i)=B(m); B(m)=d1; end for k=i+1:n for j=i+1:n A(k,j)=A(k,j)-A(i,j)*A(k,i)/A(i,i); end B