数学实验第3次作业_线性方程组的数值解法..docxVIP

数学实验第3次作业_线性方程组的数值解法..docx

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

线性方程组的数值解法一 实验目的1 学会用MATLAB软件求解线性代数方程组,对迭代法的收敛性和解的稳定性做初步分析;2 通过实例学习线性代数方程组解决简化的实际问题。二 实验内容1 已知方程组,其中,定义为:试通过迭代法求解此方程组,认识迭代法收敛的含义以及迭代初值和方程组系数矩阵性质对收敛速度的影响。实验要求:选取不同的初始向量和不同的方程组右端项向量,给定迭代误差要求,用雅可比迭代法和高斯-赛德尔迭代法计算,观测得到的迭代向量序列是否均收敛?若收敛,记录迭代次数,分析计算结果并得出你的结论。取定右端项向量和初始向量,将的主对角线元素成倍增长若干次,非对角线元素不变,每次用雅可比迭代法计算,要求迭代误差满足,比较收敛速度,分析现象并得出你的结论。初步解决:首先建立利用雅克比迭代法和高斯-赛德尔迭代法计算的函数M文件。雅克比迭代法:高斯-赛德尔迭代法:关于迭代法的收敛性:由原式化简可得,也就是说,研究此处迭代法是否收敛实际上就是研究有原矩阵经变换之后所得矩阵的性质。而由于该题第一问中,保持不变,所以对于收敛性只需研究一次。在命令栏中输入以下命令:构造题目中的矩阵:提取对角矩阵、上三角矩阵、下三角矩阵:对于两种迭代方法分别求:分别计算两个矩阵的三种范数:分别是矩阵、的1-范数、2-范数和范数。由书上定理,矩阵的谱半径不超过任何一个范数,即,而由图中可以看出,两个矩阵的六个范数没有一个大于1,所以两个矩阵的谱半径一定都小于1,所以此时两种迭代方法均收敛。利用两个函数研究原方程组:对于相同的初始值,选取不同的右端项向量。令,分别取、、、,分别观察两种迭代方法最后的迭代次数。(迭代误差拟定为0.001)1313131210998相关截图如下:从表中的数据可以看出,当初始值不变,右端项向量变化时,使用雅克比迭代法和高斯-赛德尔迭代法都可以得到最后的结果,但是明显高斯-赛德尔迭代法比雅克比迭代法的收敛速度要快。从表中还可以看出,的变化是会影响迭代次数的,因为变化以后方程组的解发生了变化,使得计算得过程也发生了变化,进而影响了迭代次数。对于相同的右端项向量,选取不同的初始值。令,分别取、、,分别观察两种迭代方法最后的迭代次数。(迭代误差拟定为0.001)13131413109109相关截图如下:从表中可以看出,初始值的变化会影响迭代次数,但是当迭代误差不大的时候,初始值的变化对迭代次数的变化的影响也不会非常大。但是通过实验可以发现,将最后一组数据的迭代误差改为0.0001是,用雅克比迭代法的迭代次数变为16,而高斯-赛德尔迭代法的迭代次数变为11。但是总体而言,高斯-赛德尔迭代法的收敛速度还是比雅克比迭代法的收敛速度快。下面研究第二问。将初始值定为,右端向量定为。改变的的主对角元素,分别变为原来的2倍、3倍、5倍、10倍,观察迭代次数变化。(用雅克比迭代法,迭代误差定为)主对角元素69153010865相关截图如下:从表格中的数据可以看出,当改变主对角元素的值时,对迭代次数的影响很大,而且在此时给出的条件下,当主对角元素的值增大时,迭代次数同步的减小。2 对于5.1.2节假设,第一年有50颗植物,且第50年后有600颗植物。试分别用追赶法、稀疏系数矩阵和满矩阵求解;若的误差,估计对结果的影响。初步解决:用三种不同的方法解决问题。首先对数据做简单处理,由题目中,第一年有50颗植物,且第50年后有600颗植物,这些条件可知,追赶法:在命令栏中输入以下命令:输出植物数量的每年数量:161.5030285736184264.0030285736184367.0781800022994470.2783314309803573.6322404310953677.1461570026443780.8277690241990884.6850768743313988.72646532554121092.96071916925781197.397042435534812102.04507839399813106.91493051577414112.01718443547415117.36293096126316122.96379018303717128.83193673110018134.98012624025219141.42172307680720148.17072938881921155.24181554266022162.65035201210123170.41244278923424178.54496038983925187.06558252930026195.99283054879227205.34610967525728215.14575120269729225.41305668646030236.17034424659531247.44099

文档评论(0)

stzs + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档