数学实验报告利用MALTAB计算线性代数方程组的数值解法..docxVIP

下载本文档

40
0
约7.43千字
约 14页
2017-01-13 发布于重庆
举报
版权申诉

数学实验报告利用MALTAB计算线性代数方程组的数值解法..docx

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

数学实验报告利用MALTAB计算线性代数方程组的数值解法.

实验三线性代数方程组的数值解法一、迭代法求解方程组㈠问题描述给定方程组的矩阵A，通过迭代法求解方程组。1、选取不同的初始向量和不同的右端项向量，给定误差要求，用两种迭代法计算；2、去顶右端项向量和初始向量，将A的主对角线元素成倍增长若干次，非主对角线元素不变，用雅克比迭代法计算。㈡方法与公式1、雅克比迭代法2、高斯-赛德尔迭代法㈢结果与分析1、不同初始向量、不同右端项向量、不同精度要求(1)初始向量定为zeros(n,1);①b=zeros(n,1)迭代次数为0，直接得到结果。②b = ones(n,1)精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1013172023262933高斯-赛德尔79111314161820③b = 1:n精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1013162023262932高斯-赛德尔79111315171820④b = n:1精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1013162023262932高斯-赛德尔68101214161819(2)初始向量定为ones(n,1)①b=zeros(n,1)精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比10411041104110411041104110411041高斯-赛德实验知，此时雅可比迭代法速度非常慢。事实上，迭代100次时，所得结果约为10^-32，已经可以认为是0，但是由于没有达到精度要求，故不算收敛。②b=ones(n,1)精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比912161922252932高斯-赛德尔79111314161820③b = 1:n精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1013161923262932高斯-赛德尔68101314161819④b = n:1精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1013161923262932高斯-赛德尔79111315171820(3)初始向量定为1:n①b=zeros(n,1)精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比10451045104510451045104510451045高斯-赛德此时又出现了雅可比迭代法收敛速度极慢的情况。②b=ones(n,1)精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1417202427303336高斯-赛德尔79111314161820③b = 1:n精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比912151922252832高斯-赛德尔68101214161819④b = n:1精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1013161922252832高斯-赛德尔79111315171820(4)初始向量定为n:1①b=zeros(n,1)精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比10451045104510451045104510451045高斯-赛德此时又出现了雅可比迭代法收敛速度极慢的情况。②b=ones(n,1)精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1417202427303336高斯-赛德尔79111314161820③b = 1:n精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比1013161922252832高斯-赛德尔68101214161819④b = n:1精度10^-310^-410^-510^-610^-710^-810^-910^-10雅可比912151922252832高斯-赛德尔79111315171820 (5)简要小结a.在个别情况下雅可比迭代法收敛速度极慢，但事实上没有达到收敛时其计算结果已经可以接受；b.要求的精度越高，迭代的次数越多，迭代的次数与所要求的精度的对数值近似呈线性，也就是说两者近似呈指数关系；b.高斯-赛德尔迭代法有着比雅可比更好的迭代特性；2、更改A的主对角线元素(1) b =20:1;初值 1:20对角线倍数12345678910迭代次数19119