matlab高斯法实验题答案.docxVIP

matlab高斯法实验题答案 matlab高斯法实验题答案 PAGE / NUMPAGES matlab高斯法实验题答案 高斯 - 若当全主元消去法解线性方程组 我们知道关于一个方程组而言， 假如做相应的变换使系数矩阵变为单位阵的话那么就能获得它的解了，高斯 - 若当法就是这样的，可是为了减少偏差，使用了全选主元法 , 下边是程序。保留为 在 matlab 中运转 function ans=gauss(a,n,b) indexC=[]; % 用于纪录主元 indexR=[]; ipiv=[]; m=1; for j=1:n ipiv(j)=0; end for i=1:n big=0; for j=1:n if(ipiv(j)~= 1) for k=1:n if abs(a(j,k))big % big=abs(a(j,k)); irow=j; icol=k; end  找寻最大的主元 end end end ipiv(icol)=ipiv(icol)+1; if(irow~=icol)% 需要调动 for l=1:n temp=a(irow,l); a(irow,l)=a(icol,l); a(icol,l)=temp; end for l=1:m temp=b(irow,l); b(irow,l)=b(icol,l); b(icol,l)=temp; end end inde