- 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章 矩阵特征值与特征向量的计算
--------学习小结
本章学习体会
通过本章的学习,我知道了求矩阵的特征值和特征向量的问题是代数计算的重要课题,在这一章,我了解到了直接计算矩阵的特征值和特征向量的MATLAB程序、间接计算矩阵的特征值和特征向量的幂法、反幂法、Jacobi方法、QR方法及MATLAB计算程序。我了解到自己对数值分析及MATLAB的掌握还很肤浅,了解到了自己的不足,同时意识到自己知识点薄弱的地方,还有对知识的理解有偏差。有的知识点理解的不透彻,自己可以动手做题,但编程实现还需要一定的编程语言知识以及数学知识和机器语言之间的转换。
四种方法各有其特点和适用范围。幂法主要用于计算矩阵按模最大的特征值及其相应的特征向量;反幂法主要用于计算矩阵按模最小的特征值及其相应的特征向量;Jacobi方法用于求实对称矩阵的全部特征值和特征向量的方法;QR方法则适用于计算一般实矩阵的全部特征值,尤其适用于计算中小型实矩阵的全部特征值。归结起来,这四种方法亦有其共同点,那就是都是用了迭代的方法来求矩阵的特征值和特征向量。此外,用MATLAB自带的解法求解特征值和特征向量也非常快速,而且不用编辑函数建立m文件。本章对于矩阵的特征值和特征向量的算法提出了新的思路,如幂法和反幂法、Jacobi、QR方法等。本章的小结主要从方法的思想,以及一些定理展开。 各种方法的运用范围1、幂法:主要用于计算矩阵按模最大的特征值和其相应的特征向量; 2、反幂法:主要计算矩阵按模最小的特征值以及其相应的特征向量; 3、Jacobi方法:用于求实对称矩阵的全部特征值和特征向量的方法;4、QR方法:适用于计算一般实矩阵的全部特征值,尤其适用于计算中小型实矩阵的全部特征值。
一、乘幂法
1、基本思想
2、一般算法
1)任意给定初始向量
2)对于k=1,2,...
3)如果 ,则
3、三种迭代公式
(1)使用范数
(2)使用范数
(3)表示的绝对值最大的分量。
二、反幂法(逆迭代)
对用乘幂法计算的按模最大的特征值与相应的
特征向量称为反幂法。
三、带原点位移的反幂法
依据:
设且
则为矩阵按模最小的特征值。
反幂法的特点
幂法和反幂法也有一定的局限性,由于幂法和反幂法的迭代是否收敛依赖于特征值的分布情况,因此实际使用时很不方便,特别是不适合于自动计算。只在矩阵阶数非常高,无法利用其他更有效的算法时,才用幂法计算按模最大的特征值和相应的特征向量,而用反幂法计算按模最小的特征值和相应的特征向量。
3.2 Jacobi方法
一、Jacobi方法的基本思想
迭代公式:
二、经典Jacobi方法的计算步骤
1、在的非对角元找按模最大的元素
2、求正交矩阵Pk使
3、控制迭代终止的条件
4、计算正交阵
三、平面旋转变换
1、初等旋转阵 (Givens矩阵)
2、初等旋转阵的性质
(1)左乘向量
(2)与矩阵相乘
左乘: 右乘: 左右乘:
四、经典Jacobi方法
实用Jacobi方法
1、按行循环消元
2、变容限循环消元法
3.3 QR方法
矩阵的QR分解
QR方法是求一般矩阵的全部特征值和特征向量的一种迭代法。
Q ——正交矩阵 R——上三角矩阵
1、Householder矩阵(镜面反射阵)
H为对称正交矩阵
设有非零向量 和单位向量 是实数并且 ,用反幂法求矩阵
解:对A-2.93I作三角分解得
按算法迭代 ,与准确值 , 与准确值 比较,残差 。
原创力文档


文档评论(0)