如果您无法下载资料，请参考说明：
1、部分资料下载需要金币，请确保您的账户上有足够的金币
2、已购买过的文档，再次下载不重复扣费
3、资料包下载后请先用软件解压，在使用对应软件打开
functionx=Gauss_lzy(A,b)%求解方程组的列主元Gauss消元法%A为方程组系数矩阵，b为方程组右端向量，x为方程组的解[n,m]=size(A);nb=length(b);ifn~=merror('%系数矩阵必须是方的！');endifm~=nberror('%b的维数与方程的行数不匹配！');endfork=1:n-1%选主元a_max=0;fori=k:nifabs(A(i,k))>a_maxa_max=abs(A(i,k));r=i;endendifa_max<1e-15error('%系数矩阵奇异，无法求解方程组');end%交换两行ifr>kforj=k:nz=A(k,j);A(k,j)=A(r,j);A(r,j)=z;endz=b(k);b(k)=b(r);b(r)=z;end%消元过程fori=k+1:nm=A(i,k)/A(k,k);forj=k+1:nA(i,j)=A(i,j)-m*A(k,j);endb(i)=b(i)-m*b(k);endend%回代过程ifabs(A(n,n))<1e-15error('%系数矩阵奇异，无法求解方程组！');endx=zeros(size(b));fork=n:-1:1forj=k+1:nb(k)=b(k)-A(k,j)*x(j);endx(k)=b(k)/A(k,k);end
立即下载
