实验二：Gauss列主元消去法.doc

下载文档 降价啦

4
0
约3.62千字
约 3页
2017-08-06 发布于浙江
举报
版权申诉
保障服务

实验二：Gauss列主元消去法.doc

1、本文档共3页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

实验二：Gauss列主元消去法

实验二：Gauss列主元消去法程序1：Gauss列主元消去法 A=input(请输入线性方程组的增广矩阵A=); n=length(A)-1; x=zeros(n,1); aa=zeros(n,1); for j=1:n for i=1:(n+1) AA(j,i)=abs(A(j,i)); end end for k=1:(n-1) for i=k:n aa(i-(k-1))=AA(i,k); end for i=k:n if AA(i,k)==max(aa) break end end if AA(i,k)==0 break fprintf(方程组系数矩阵奇异\n); else for j=k:(n+1) jh=A(i,j); A(i,j)=A(k,j); A(k,j)=jh; end end fenzi=A(k,k); for j=k:(n+1) A(k,j)=A(k,j)/fenzi; end for p=(k+1):n jj=A(p,k); for j=k:(n+1) A(p,j)=A(p,j)-jj*A(k,j); end end end if k==(n-1) x(n)=A(n,(n+1))/A(n,n); for i=(n-1):(-1):1 he=0; for j=(i+1):n he=he+A(i,j)*x(j); end x(i)=A(i,(n+1))-he; end end x 用Gauss列主元消去法解方程组： 1．请输入线性方程组的增广矩阵A=[1e-008,2,3,1;-1,3.172,4.623,2;-2,1.072,5.643,3] x = -0.4653 -0.0700 0.3800 2．请输入线性方程组的增广矩阵A=[4,-2,4,10;-2,17,10,3;-4,10,9,-7]; x = 2.9464 0.6071 -0.1429 3．请输入线性方程组的增广矩阵A=[0.3e-020,1,0.7;1,1,0.9] x = 0.2000 0.7000 程序2：不选主元的高斯消去法 A=input(请输入线性方程组的增广矩阵A=); n=length(A)-1; x=zeros(n,1); for k=1:(n-1) if A(k,k)==0 break fprintf(方程组不能用普通的高斯消去法解\n); else fenzi=A(k,k); for j=k:(n+1) A(k,j)=A(k,j)/fenzi; end for p=(k+1):n jj=A(p,k); for j=k:(n+1) A(p,j)=A(p,j)-jj*A(k,j); end end x(n)=A(n,(n+1))/A(n,n); for i=(n-1):(-1):1 he=0; for j=(i+1):n he=he+A(i,j)*x(j); end x(i)=A(i,(n+1))-he; end end end x 用不选主元的Gauss消去法解方程组： 1．请输入线性方程组的增广矩阵A=[4,-2,4,10;-2,17,10,3;-4,10,9,-7]; x = 2.9464 0.6071 -0.1429 2．请输入线性方程组的增广矩阵A=[1e-008,2,3,1;-1,3.172,4.623,2;-2,1.072,5.643,3]; x = -0.4653 -0.0700 0.3800 3．请输入线性方程组的增广