计算方法(李桂成著),实验报告,矩阵特征值问题计算.docVIP

计算方法(李桂成著),实验报告,矩阵特征值问题计算 特征值 矩阵 计算方法 实验 计算 如何计算矩阵特征值 excel计算矩阵特征值 篇一：矩阵特征值与特征向量的计算上机报告 矩阵特征值和特征向量的计算 （一）问题 设： 621A= 231 111 利用反幂法求这一矩阵在2附近的特征值和对应的特征向量。 （二）解决问题算法 题设中已经要求使用反幂法求矩阵的特征值和特征向量。|A|=9≠0，故A-1存在，可以使用反幂法求解。下面回顾反幂法： 反幂法，就是应用幂法于A-1上求A的模的最小特征值和对应的特征向量。因此其基本的迭代格式： Ayk=zk?1 μk=yk模的最大分量 yk zk= k 若A的特征值为： |λn |λn?1 ≤?≤|λ1| 1 则{zk}收敛于A对应的λn的一个特征向量，{μk}则收敛于λ?n，其收敛速度依赖于|λn/λn?1|的大小。 i，应用反幂法与A-λ iI上来求其对应的某个特征向在求得A的某个特征值λi的近似值λ 量的： iI）yk=zk?1 （A?λ zk= yk k i是λi很好的近似时，只需迭代一次就可得到很好的近似特征向量。 当λ 就本题而言，先计算f(y)=|A- yI |，确定其零点的大致分布： 可见A的三个特征值互不相等且均为实数，分别在0.5，2，7附近。 所以只使用一次反幂法只能求出在0.5附近的特征值和特征向量，为了求出题设要求的在2附近的特征值和特征向量，必须要再使用收缩技巧。 一个简单实用的收缩技巧是正交变换法：假设Ax1=λ1x1，先计算一个Householder变换H使得： Hx1=αe1 因为本题的特征值均为实数，由实数域的封闭性，每一个特征值对应的特征向量分量也均为实数。然后可以得到HAH有如下形式： λ HAH′= 1 ? A1 由此A1是n-1阶矩阵，包含A的其余n-1个特征值，再用反幂法于A1，就可以得到A的另一个特征值，即为要求的特征值，然后使用反幂法求出要求的特征向量即可。 也可以使用带位移的反幂法先求出2附近的特征向量，在求出特征值。 （三）使用的软件 IDL （四）数值结果 所求出的在2附近的特征值：2.13307 该特征值对应的三维特征向量： (?0.497425,?0.819589,0.284327)T （五）数值结果分析 1 由上表可知，A的模最小的特征值为对应的特征向量为(?0.04615，?0.37492，1.00000)T 因为|λn/λn?1|≈3.7，所以上面使用反幂法收敛速度还是比较快的，迭代9次以后得到的近似值就在精度要求的范围内与实际值相等。 2 由上表可知，题设中要求的A的的特征值为迭代11次以后得到的近似值就在精度要求的范围内与实际值相等。 因为 λn?1 ≈3.4|λ n?2 λ λn n?1 |，所以在相同精度要求下第二次迭代需要的次数比第一次多，也说 λn n?1 明了反幂法的收敛速度与 λ 是成正相关的。 3、使用反幂法求特征值为2附近的特征向量： 当位移量选为2当位移量选为2.13307时： 迭代一次就可以得到很好的特征向量。 4、直接用带位移的反幂法求解： 当位移量选为2时： 得到特征向量u=从而Au=(?1.06104,?1.74824,0.606491)T 由Au=λu，可得λ=2.13307 可见直接使用带位移的反幂法，可以非常简便地求出要求的特征值与特征向量。 篇二：矩阵特征值问题的计算 本科毕业论 论文题目： 矩阵特征值问题的计算 学生姓名： 项田鹏 学号： 200500820328 专业： 信息与计算科学 指导教师： 尹 哲 学 院： 数学科学学院 2009 年 05月27 日 文 毕业论文（设计）内容介绍 目 录 中文摘要…………………………………………………………………3 英文摘要…………………………………………………………………3 一．引言…………………………………………………………………4 二．幂法…………………………………………………………………4 （一）幂法的算法及收敛定理…………………………………………4 （二）举例………………………………………………………………6 （三）对幂法的分析……………………………………………………7 二．反幂法………………………………………………………………7 （一）反幂法的算法……………………………………………………7 （二）举例………………………………………………………………8 （三）对反幂法的分析…………………………………………………9 三．雅可比法……………