数值分析实验作业,gauss消去法的数值稳定性分析.pdf

下载文档 降价啦

33
0
约1.25万字
约 9页
2020-09-01 发布于境外
举报
版权申诉
保障服务

数值分析实验作业,gauss消去法的数值稳定性分析.pdf

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

实验 3.1 Gauss 消去法的数值稳定性试验实验目的 : ( k ) 观察和理解 Gauss消元过程中出现小主元 ( 即 a kk 很小 ) 时引起的方程组解的数值不稳定性。实验内容：求解方程组 Ax b ，其中 - 15 0.3 10 59. 14 3 1 59 .17 5.291 - 6.130 - 1 2 46.78 （1）A1 ， b1 ； 11.2 9 5 2 1 1 2 1 1 2 10 7 0 1 8 3 2.099999999999 6 2 5.900000000001 A2 b2 5 1 5 1 5 0 1 0 2 1 （2） , . 实验要求：（1）计算矩阵的条件数，判断系数矩阵是良态的还是病态的。 4 x ,x R （2）用 Gauss列主元消去法求得 L 和 U及解向量 1 2 . ~ ~ ~ ~ 4 x ,x R （3）用不选主元的 Gauss消去法求得 L 和 U 及解向量 1 2 . （4 ）观察小主元并分析其对计算结果的影响 . 程序如下 : 计算矩阵条件数及 Gauss列主元消去法 : format longeng A1=[0.3e-15 59.14 3 1;5.291 -6.130 -1 2;11.2 9 5 2;1 2 1 1]; b1=[59.17;46.78;1;2]; n=4; k2=cond(A1) %k2 为矩阵的条件数； for k=1:n-1 a=max(abs(A1(k:n,k))); [p,k]=find(A1==a); B=A1(k,:);c=b1(k); A1(k,:)=A1(p,:);b1(k)=b1(p); A1(p,:)=B;b1(p)=c; if A1(k,k)~=0 A1(k+1:n,k)=A1(k+1:n,k)/A1(k,k); A1(k+1:n,k+1:n)=A1(k+1:n,k+1:n)-A1(k+1:n,k)*A1(k,k+1:n); else break end end L1=tril(A1,0); for i=1:n L1(i,i)=1; end L=L1 U=triu(A1,0) for