北航数值分析大作业第二次的.docxVIP

  1. 1、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
北航数值分析大作业第二次的 北航数值分析大作业第二次的 PAGE 北航数值分析大作业第二次的 数值分析 第二次大作业 姓名: 路子威 学号: S1 题目:使用带双步位移的QR分解法求矩阵的全部特征值,并对其中的每一个实特征值求相应的特征向量。已知: (i,j=1,2,……,10) 要求: 用幂法、反幂法和QR方法求矩阵的特征值时,要求迭代精度水平为e=10-12。 打印以下内容: 全部源程序; 矩阵A经过拟上三角化后所得的矩阵A(n-1); 对矩阵A(n-1)实行QR方法迭代结束后所得的矩阵; 矩阵A的全部特征值λi=(RI,Ii);(i=1,2,…,10) A的相应于实特征值的特征向量。 采用e型输出实型数,并且至少显示12位有效数字。 算法: 输入需要求解的矩阵A 对上述生成的矩阵进行拟上三角化 对拟上三角化后的矩阵进行QR分解(带双步位移)。 用函数characteristic()求解矩阵A(n-1)即A的所有特征值。 用函数characteristicvector()求解矩阵A(n-1)即A的所有特征向量。 输出A的所有特征值λ、A的所有实特征值对应的特征向量、拟上三角矩阵A(n-1)、及其Q、R和R*Q 说明: 为了减少求特征值和特征向量过程中的计算量,在对矩阵进行QR分解前先进行拟上三角化。拟上三角化之后矩阵中包含有大量的0元素,在进行QR迭代时可以减少大量的0元素乘法运算,节省计算时间。且上三角矩阵进行QR分解之后,通过RQ得到新的矩阵A仍然是上三角矩阵,因此可以在整个迭代过程中使用。而带双步位移的QR分解法可以加快收敛速度,使矩阵在更短的迭代次数下完成收敛。 源程序: #include #include #include #define e /*设置精度水平*/ #define N 10 /*设置矩阵阶数*/ #define L 2500 /*设置迭代次数*/ /*子程序*/ /*sgn函数* 0. 0 -1. 018099 -0. 0 1. 0. -0. 0. 0 832816 -0. 0 0 0 -0. 0 0 -0. 0 0 55235 0 0 0 生成的上三角矩阵R: 1. 0. 0 -0. 0 054472 0 2273214 0. 0039456 0 0 -0. 0 0 0 546304 0 0 0 56156e-017 0591709 0 0 0 054775e-028 0 -0. 0. 0 0 0 06929e-108 0 0 0 0 0 0 0 0 0 生成的正交矩阵Q: -1 365568e-273 0 0 99121e-208 0 0 10 0 0 0 -1 797898e-100 192 0 0 0 0

文档评论(0)

150****7735 + 关注
实名认证
文档贡献者

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

1亿VIP精品文档

相关文档