matlab高斯法实验题答案.docVIP

下载本文档

1
0
约4.07千字
约 8页
2017-08-24 发布于河南
举报
版权申诉

matlab高斯法实验题答案.doc

1、本文档共8页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

高斯-若当全主元消去法解线性方程组我们知道对于一个方程组而言，如果做相应的变换使系数矩阵变为单位阵的话那么就能得到它的解了，高斯-若当法就是这样的，但是为了减少误差，使用了全选主元法,下面是程序。保存为gauss.m在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 indexR(i)=irow; indexC(i)=icol; if(a(icol,icol)==0) disp(error); end pivinv=1/a(icol,icol); for l=1:n a(icol,l)=a(icol,l)*pivinv; end for l=1:m b(icol,l)=b(icol,l)*pivinv; end for ll=1:n if(ll~=icol) dum=a(ll,icol); a(ll,icol)=0.0; for l=1:n a(ll,l)=a(ll,l)-a(icol,l)*dum; end for l=1:m b(ll,l)=b(ll,l)-b(icol,l)*dum; end end end end for l=n:-1:1 %已经进行过交换，现在交换回来 if(indexR(l)~=indexC(l)) for k=1:n temp=a(k,indexR(l)); a(k,indexR(l))=a(k,indexC(l)); a(k,indexC(l))=a(k,indexR(l)); end end end ans=b; 下面解这个简单的方程组 a=[1 0 0 0;0 2 0 0 ;0 0 3 0 ;0 0 0 4] a = 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 b=[10;40;60;80] b = 10 40 60 80 gauss(a,4,b) ans = 10 20 20 20 昨天才回答过这个问题..你可以再搜搜的 Gauss消去法的分析。其包括两个过程：消去过程：把方程组系数矩阵A化为同解的上三角矩阵；回代过程：按相反的顺序，从xn至x1逐个求解上三角方程组。 %高斯消去法的MATLAB