六线性方程组超松弛迭代方法.docVIP

下载本文档

1
0
约3.11千字
约 4页
2017-06-12 发布于重庆
举报
版权申诉

六线性方程组超松弛迭代方法.doc

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

六线性方程组超松弛迭代方法

1、用谱半径判别超松弛迭代法产生的迭代序列的脸敛散性。 function H=ddpbj(A,om) D=diag(diag(A));U=-triu(A,1); L=-tril(A,-1);iD=inv(D-om*L); B2=iD*(om*U+(1-om)*D); H=eig(B2);mH=norm(H,inf); if mH=1 disp(谱半径不小于1，序列发散，谱半径和特征值如下：) else disp(谱半径小于1，序列收敛，谱半径和特征值如下：) end mH A=[5 1 -1 -2;2 8 1 3;1 -2 -4 -1;-1 3 2 7]; H=ddpbj(A,1.15) 谱半径小于1，序列收敛，谱半径和特征值如下： mH = 0.1596 H = 0.1049 + 0.1203i 0.1049 - 0.1203i -0.1295 + 0.0556i -0.1295 - 0.0556i H=ddpbj(A,5) 谱半径不小于1，序列发散，谱半径和特征值如下： mH = 14.1082 H = -14.1082 -2.5107 0.5996 + 2.6206i 0.5996 - 2.6206i 2、超松弛迭代法的程序。 function X=sordd(A,b,X0,P,om,wucha,max1) D=diag(diag(A));U=-triu(A,1);L=-tril(A,-1);dD=det(D); if dD==0 disp(因为对角矩阵D奇异，所以此方程组无解。) else disp(因为对角矩阵D非奇异，所以此方程组有解。) iD=inv(D-om*L);B2=iD*((1-om)*D+om*U);f2=om*iD*b;jX=A\b;X=X0; [n m]=size(A); for k=1:max1 X1=B2*X+f2;djwcX=norm(X1-X,P); xdwcX=djwcX/(norm(X,P)+eps); if (djwcXwucha)|(xdwcXwucha) return else k,X1,k=k+1;X=X1; end end if (djwcXwucha)|(xdwcXwucha) disp(收敛，A的分解矩阵D，U,L和方程组的精确解jX和近似解X如下;) else disp(没有达到给定的精度，迭代次数超过最大次数，方程组的精确解jX和近似解X如下;) X=X;jX=jX end end X=X;D,U,L,jX=jX A=[5 1 -1 -2;2 8 1 3;1 -2 -4 -1;-1 3 2 7]; b=[4;2;6;-3]; X0=[0;0;0;0]; X=sordd(A,b,X0,inf,1.15,0.00001,10) 因为对角矩阵D非奇异，所以此方程组有解。 k = 1 ans = 0.9200 0.0230 -1.4737 0.1312 k = 2 ans = 0.4981 0.2961 -1.5687 -0.0612 k = 3 ans = 0.3882 0.3834 -1.5809 -0.0894 k = 4 ans = 0.3689 0.3898 -1.5802 -0.0917 k = 5 ans = 0.3694 0.3896 -1.5794 -0.0915 k = 6 ans = 0.3697 0.3893 -1.5793 -0.0913 k = 7 ans = 0.3698 0.3892 -1.5793 -0.0913 X = 0.3698 0.3892 -1.5793 -0.0913 X=sordd(A,b,X0,inf,5,0.00001,10) 因为对角矩阵D非奇异，所以此方程组有解。 k = 1 ans = 4.0000 -3.7500 6.8750 -1.0714 k = 2 ans = -3.5179 18.3594 -83.9565 80.2264 k =