原创力文档- 1、本文档共13页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
北航数值分析大作业一(纯原创,高分版)
数值分析上机实习作业一
题目：
设有501×501 的实对称矩阵A ，
1
⋱ ⋱ ⋱
= ⋱ ⋱ ⋱
⋱ ⋱ ⋱
[ 501 ]
0.1
( ) ( ) ( )
其中 = 1.64 − 0.024 sin 0.2 − 0.64 , = 1,2, ⋯ ,501 , = 0.16, = −0.64。矩
阵A 的特征值为， ( 1,2, ⋯ ,501, 并且有
| | | |
1 ≤ 2 ≤ ⋯ ≤ 501, =
1≤≤501
1. 求λ , 和λ 的值。
1 501
2. 求A 的与数 = 1 + 501−1最接近的特征值 ( = 1,2, ⋯ ,39).
40
3. 求A 的（普范数）条件数cond(A)和行列式detA 。
说明
1. 在所用的算法中，凡是要给出精度水平的ε，都取 = 10−12 。
2. 选择算法的时候应使矩阵A 的所有零元素都不存储。
3. 打印以下内容：
(1)算法设计方案和思路。
(2)全部源程序。
( ) ( )
(3)特征值 , , , = 1,2, ⋯ ,39 以及 , det 的值(采用e 型输出实型数，
1 501 2
并且至少显示12 位有效数字) 。
(4)发现的现象与问题。
算法设计思路和方案
因为题目只要求求部分特征值，及最大的特征值和最小的特征值以及距离某些给定实数
最近的特征值，而矩阵 A 是实对称矩阵，所有特征值均为实数，故我们可以用幂法和反幂
法结合适当的平移量求解。其中求解最大的特征值和最小的特征值，采用幂法并结合适当的
平移量求解，而求解λ ，和距离μ 最近的特征值采用带平移的反幂法求解。
s k
算法分析
一、幂法求解,
| |
我们知道幂法适用于求解矩阵 A 的绝对值最大的特征值，设矩阵A 的特征值为λ1 ≥
| | | | | |
文档评论（0）