数值实验3_3病态的线性方程组求解.docVIP

下载本文档

30
0
约2.31千字
约 2页
2017-10-10 发布于重庆
举报
版权申诉

数值实验3_3病态的线性方程组求解.doc

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

数值实验3_3 病态的线性方程组求解自动化系李琳琳 2004211068 选择问题的维数为6，分别用Gauss消去法、J迭代方法、GS迭代方法和SOR迭代方法求解并与问题的解比较，所得结果见下表 n=6 Gauss法 J迭代 GS迭代 SOR迭代实际解X X(1) 1.0e+005*0.0021 1.0e+308* -0.7702 0.9999 1.0000 1 X(2) 1.0e+005*-0.0647 0e+308*-1.6621 1.0009 1.0003 1 X(3) 1.0e+005*0.4598 -Inf 0.9981 1.0010 1 X(4) 1.0e+005*-1.2366 -Inf 0.9974 0.9922 1 X(5) 1.0e+005 *1.3981 -Inf 1.0089 1.0126 1 X(6) 1.0e+005*-0.5610 -Inf 0.9947 0.9939 1 逐步增大问题的维数，仍用上述方法求解，结果如下： n=8 Gauss法 J迭代 GS迭代 SOR迭代实际解X X(1) 1.0e+008*-0.0002 1.0e+308 * 1.1898 1.0001 1.0002 1 X(2) 1.0e+008 *0.0094 Inf 0.9974 0.9964 1 X(3) 1.0e+008*-0.1220 Inf 1.0136 1.0179 1 X(4) 1.0e+008 * 0.6586 Inf 0.9794 0.9726 1 X(5) 1.0e+008*-1.7698 Inf 0.9982 1.0031 1 X(6) 1.0e+008 *2.5001 Inf 1.0141 1.0130 1 X(7) 1.0e+008*-1.7767 Inf 1.0101 1.0091 1 X(8) 1.0e+008 * 0.5007 Inf 0.9870 0.9877 1 N=10 Gauss法 J迭代 GS迭代 SOR迭代实际解X X(1) 1.0e+011*0.0000 1.0e+308*-0.5027 1.0001 1.0001 1 X(2) 1.0e+011*-0.0012 1.0e+308*-1.1763 0.9982 0.9989 1 X(3) 1.0e+011 *0.0254 1.0e+308 -1.6465 1.0056 1.0022 1 X(4) 1.0e+011*-0.2342 -Inf 0.9993 1.0042 1 X(5) 1.0e+011 *1.1282 -Inf 0.9916 0.9903 1 X(6) 1.0e+011 *1.1282 -Inf 0.9968 0.9965 1 X(7) 1.0e+011 * 5.1657 -Inf 1.0052 1.0042 1 X(8) 1.0e+011*-5.0212 -Inf 1.0087 1.0080 1 X(9) 1.0e+011 *2.6492 -Inf 1.0039 1.0039 1 X(10) 1.0e+011*-0.5851 -Inf 0.9904 0.9916 1 N=15: Gauss法 J迭代 GS迭代 SOR迭代实际解X X(1) 1.0e+015*-0.0000 1.0e+308 * 1.6060 0.9999 0.9998 1 X(2) 1.0e+015* 0.0000 Inf 1.0035 1.0036 1 X(3) 1.0e+015* 0.0000 Inf 0.9821 0.9830 1 X(4) 1.0e+015*-0.0006 Inf 1.0236 1.0223 1 X(5) 1.0e+015* 0.0089 Inf 1.0074 1.0035 1 X(6) 1.0e+015*-0.0656 Inf 0.9902 0.9937 1 X(7) 1.0e+015* 0.2681 Inf 0.9856 0.9893 1 X(8) 1.0e+015*-0.6188 Inf 0.9905 0.9924 1 X(9) 1.0e+015* 0.6456 Inf 0.9990 0.9984 1 X(10) 1.0e+015* 0.4318 Inf 1.0068 1.0045 1 X(11) 1.0e+015*-2.3874 Inf 1.0113 1.0083 1 X(12) 1.0e+015* 3.5025 Inf 1.0113 1.0089 1 X(13) 1.0e+015*-2.7037 Inf 1.0066 1.0056 1 X(14) 1.0e+015* 1.1122 Inf 0.