Guass列主元消去法.doc

Guass列主元消去法 （已验证过，保证可用正确） function y=Guass(A,b) n=length(A); u=zeros(1,n); for k=1:n-1 p=find(abs(A(k:n,k))==max(abs(A(k:n,k))));% p(1)=p(1)+k-1; u(k)=p(1); dd=A(p(1),1:n); % 中间变量 A(p(1),1:n)=A(k,1:n); A(k,1:n)=dd; if A(k,k)~=0 A(k+1:n,k)=A(k+1:n,k)/A(k,k); A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n); else stop % 矩阵奇异 end end %-------------------- U=zeros(n); for i=1:n for j=i:n U(i,j)=A(i,j); end end L=A-U; for i=1:n L(i,i)=1; end P=eye(n); for i=1:n-1 dd=P(i,:); % 中间变量 P(i,:)=P(u(i),:); P(u(i),:)=dd; end [n,m]=size(L); b1=P*b; for j=1:n-1 b1(j)=b1(j)/L(j,j); b1(j+1:n)=b1(j+1:n)-b1(j)*L(j+1:n,j); end b1(n)=b1(n)/L(n,n); [m,n]=size(U); y=b1; for j=n:-1:2 y(j)=y(j)/U(j,j); y(1:j-1)=y(1:j-1)-y(j)*U(1:j-1,j); end y(1)=y(1)/U(1,1); 示例：求9阶的Ax=b的解Xs为用Guass列主元消去法得到的近似解，可以看出在这种情况下，用Guass列主元消去法得到解与精确解的误差10-14左右，因为本例才用的是随机数，所以数据和误差会不同 在matlab命令窗口中输入 A=rand(9) A = 0.844392156527205 0.462449159242329 0.721758033391102 0.269061586686018 0.644442781431336 0.607303940685635 0.832916819075216 0.939829470344921 0.993704624120852 0.344462411301042 0.424349039815375 0.473485992965320 0.765500016621438 0.647617630172684 0.450137696965896 0.256440992229147 0.645551874972524 0.218676632399634 0.780519652731358 0.460916366028964 0.152721200438232 0.188661976791491 0.679016754093202 0.458725493648868 0.613460736812875 0.479463224948888 0.105798273250228 0.675332065747000 0.770159728608609 0.341124607049109 0.287498173066131 0.635786710514084 0.661944751905652 0.582249164527227 0.639316961040108 0.109697464523194 0.006715314318477 0.322471807186779 0.607389213768347 0.091113463686535 0.945174113109401 0.770285514803660 0.540739337124410 0.544716110526763 0.063591370975106 0.602170487581795 0.784739294760742 0.191745255461798 0.576209380663007 0.208934922426023 0.