原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
线性方程组迭代法2
线性方程组的迭代法上机试验1 问题描述2 问题分析迭代法部分上机题有三个要求：1、SOR迭代法求解。2、共轭梯度法求解。3、迭代法和直接方法结果比较。3 求解及结果3.1 SOR迭代法求解及结果为了方便后面跟直接解法比较，针对n=2、5、8、10、15等5种情况进行数值实验，取1.25结果如下：n=2时，迭代次数k=30；nx分量SOR法求解绝对误差(1.0e-15 *)211.0000 0.444121.0000 -0.6661n=5时，迭代次数k= 55464;nx分量SOR法求解绝对误差511.0000 021.0002 0.000230.9993 -0.000741.0011 0.001150.9995 -0.0005n=8时，迭代次数k= 113389；nx分量SOR法求解绝对误差811.0000 0.0000 20.9998 -0.0002 31.0006 0.0006 41.0011 0.0011 50.9954 -0.0046 61.0021 0.0021 71.0039 0.0039 80.9971 -0.0029 n=10时，迭代次数k= 58951；nx分量SOR法求解绝对误差1011.0000 0.000021.0006 0.0000 30.9955 0.0006 41.0118 -0.0045 50.9915 0.0118 60.9950 -0.0085 71.0013 -0.0050 81.0059 0.0013 91.0039 0.0059 100.9945 0.0039 n=15时，迭代次数k= 132512；nx分量SOR法求解绝对误差1511.0000 0.0000 20.9999 -0.0001 31.0017 0.0017 40.9929 -0.0071 51.0078 0.0078 61.0025 0.0025 70.9980 -0.0020 80.9957 -0.0043 90.9966 -0.0034 100.9994 -0.0006 111.0022 0.0022 121.0039 0.0039 131.0036 0.0036 141.0006 0.0006 150.9950 -0.0050 3.2 CG迭代法求解及结果为了方便后面跟直接解法比较，针对n=2、5、8、10、15等5种情况进行数值实验，允许误差esp取值1.0e-6，结果如下：n=2时，迭代次数k=2：nx分量CG法求解绝对误差（1.0e-14 *）211.0000 0.288721.0000 0.1554n=5时，迭代次数k=4：nx分量CG法求解绝对误差511.0000 021.0005 0.000530.9976 -0.002441.0036 0.003650.9983 -0.0017n=8时，迭代次数k=6：nx分量CG法求解绝对误差 811.0000 020.9997 -0.000331.0014 0.001440.9982 -0.001850.9994 -0.000661.0014 0.001471.0015 0.001580.9985 -0.0015n=10时，迭代次数k=6：nx分量CG法求解绝对误差10110.0000 20.99920.0006 31.0031-0.0045 40.99780.0118 50.9973-0.0085 60.9998-0.0050 71.00220.0013 81.00290.0059 91.0010.0039 100.9966-0.0055n=15时，迭代次数k=5：nx分量CG法求解绝对误差1511.00020.000220.9967-0.003331.00830.008341050.9945-0.005560.9942-0.005870.9969-0.003181.00060.000691.00370.0037101.00550.0055111.00580.0058121.00430.0043131.00130.0013140.9968-0.0032150.9909-0.00913.2 迭代法与直接法比较1、n=8时几种方法的比较：n=8时，直接解法的误差控制在1e-7之内；而迭代法的误差达到1e-3，误差更大。且SOR法计算时间较直接解法长，其迭代次数为K=113389。2、n=15时几种方法的比较：n=15时，直接解法出现了较大的误差，超过了300%，其结算结果几乎不可用了，而SOR法和CG法都能得到比较好的结果。误差控制在1e-3之内，而SOR法和CG法求解的精度相差无几，但是CG法收敛的更快，只需要小于10次迭代就可以收敛，而SOR法则需要计算10万次以上才能计算出满足要求的解。如图1所示：3、直接法，在当n=15以后，即方程高度病态以后，就求解不出精
文档评论（0）