数值分析-矩阵特征值问题计算课件.pptx

1、1第八章第八章 矩阵特征值问题计算矩阵特征值问题计算 对对n 阶方阵阶方阵A求数求数 和非零向量和非零向量x ,使其满足使其满足Ax= x 这样的这样的 值称为矩阵值称为矩阵A的特征值，非零向量的特征值，非零向量 x 称为矩称为矩阵阵A的与特征值的与特征值 相对应的一个特征向量。相对应的一个特征向量。2定义定义1 设矩阵设矩阵A, B R n n，若有可逆阵若有可逆阵P，使使 则称则称A与与B相似相似。APPB1定理定理1 若矩阵若矩阵A, B R n n且相似且相似，则则（1）A与与B的特征值完全相同的特征值完全相同；（2）若若x是是B的特征向量的特征向量，则则Px便为便为A的特征向量的特征

2、向量。8.1 预备知识预备知识3定理定理2： 设设A R n n具有完全的特征向量系，即存在具有完全的特征向量系，即存在n个个线性无关线性无关 nDAPP 211其中其中 i为为A的特征值的特征值，P的各列为相应于的各列为相应于 i的特征向量的特征向量。 的特征向量构成的特征向量构成Rn的一组基底的一组基底，则经相似变换可化则经相似变换可化A为为对角阵，即有可逆阵对角阵，即有可逆阵P，使使4定理定理3 ：A R n n， 1, , n为为A的特征值的特征值，则则 niiniiiaAtr11)( （2）A的行列式值等于全体特征值之积，即的行列式值等于全体特征值之积，即nA 21)det( （1）

3、A的迹数等于特征值之和，即的迹数等于特征值之和，即5定理定理46定理定理5 设设A R n n为对称矩阵为对称矩阵，其特征值其特征值 1 2 n，则则 （1）对任意对任意A R n，x0，1),(),( xxxAxn),(),(min0 xxxAxxn（2）),(),(max01xxxAxx（3）7定理定理6 （Gerschgorin圆盘定理圆盘定理） 设设A R n n，则则niaaznijjijii, 2, 1,1 表示以表示以aii为中心为中心，以以 半径为的复平面上的半径为的复平面上的n个圆盘个圆盘。 nijjija1（2）如果矩阵如果矩阵A的的m个圆盘组成的并集个圆盘组成的并集S（连

4、通的连通的）与其余与其余（1）A的每一个特征值必属于下述某个圆盘之中的每一个特征值必属于下述某个圆盘之中，n m个圆盘不连接个圆盘不连接，则则S内恰包含内恰包含m个个A的特征值的特征值。 8910111213定理定理71415一一 幂法幂法1 基本思想基本思想01122nnva xa xa x 任取非零向量任取非零向量 v0 , 则可则可唯一表示为唯一表示为 设设n 阶矩阵阶矩阵A 的特征值的特征值 , 满足满足 ,且其对应有且其对应有n个线性无关的特个线性无关的特 征向量征向量 x1 , x2, , xn ，即，即), 2 , 1(nixAxiii), 2 , 1(nii |21n 8.2

5、幂法和反幂法幂法和反幂法16则则kknknnkknknnkknnkkkxaxaxaxaxaxaxaxaxaxaAvAv1111212211122211122110nnxaxaxav2211017其中其中nknnkkxaxa12122由假设条件由假设条件 , 211njj 所以当所以当k充分大时，有充分大时，有111xavkk1 11limkkkva x从而从而 0lim kk 所以所以18即为矩阵即为矩阵 A 的对应特征值的对应特征值 1 的近似特征向量。的近似特征向量。用用 (vk)i 表示表示 vk 的第的第 i 个分量个分量,则当则当k充分大时充分大时,有有 11ikikvv即为主特征值

6、的近似值。即为主特征值的近似值。111xavkkkkkkvxaAvv111111且且19nnRAn1设设有有主特征值主特征值满足满足个线性无关的特征向量，个线性无关的特征向量，n321则对任意非零初始向量则对任意非零初始向量，下面的式子成立，下面的式子成立定理定理1 11limkkkva x11)()(limikikkvv00v20 迭代公式迭代公式 (1)实质上是由矩阵实质上是由矩阵A的乘幂的乘幂 Ak与非零与非零向量向量 v0 相乘来构造向量序列相乘来构造向量序列 xk , 从而计算主特征从而计算主特征值及其对应的主特征向量值及其对应的主特征向量,故称这种方法为幂法。故称这种方法为幂法。2

7、1nnRAn1设设有有主特征值主特征值满足满足个线性无关的特征向量，个线性无关的特征向量，n321则对任意非零初始向量则对任意非零初始向量00v定理定理按照下述方法构造的向量序列按照下述方法构造的向量序列 ,ku kv), 2 , 1(),max(,0100kvuvAuvuvkkkkkkk)max(lim) 1 (11xxukk1lim)2(kk则有则有22 2. 幂法实用计算公式幂法实用计算公式), 2 , 1(),max(,0100kvuvAuvuvkkkkkkk)max(lim) 1 (11xxukk1lim)2(kk23例例1 求矩阵求矩阵 1634310232A的主特征值与其对应的特

8、征向量。的主特征值与其对应的特征向量。解取解取 v0=(0,0,1)T , 则则TTvuv25. 0 , 1, 5 . 01, 4,1 , 4 , 211111), 2 , 1(),max(,0100kvuvAuvuvkkkkkkk24直到直到k=8 时的计算结果见下表时的计算结果见下表0.5, 1, 0.750011.00005.5000, 11.0000, 8.250080.5, 1, 0.750011.00055.5002, 11.0005, 8.250170.5, 1, 0.750110.99745.4987, 10.9974, 8.249460.5, 1, 0.749411.0142

9、5.5075, 11.0142, 8.257650.5, 1, 0.753610.92235.4621, 10.9223, 8.230640.5, 1, 0.736011.44445.7222, 11.4444, 8.36130.5, 1, 0.861194.5, 9, 7.7520.5, 1, 0.25 4 2, 4, 1,1k从而从而kukvk25二、幂法的加速二、幂法的加速1、原点平移法、原点平移法 如果如果 是矩阵是矩阵 A 的特征值的特征值,则对任意的实数则对任意的实数p, 矩阵矩阵 A-pE 的特征值为的特征值为 -p,且且 A 与与 A-pE 的特征向量相同的特征向量相同.据此据

10、此, 如果要计算如果要计算 A 的主特征值的主特征值 1 , 只要选择合适的只要选择合适的数数 p,使使 1-p 为矩阵为矩阵 A-pE 的主特征值的主特征值,且且 1212max ppini那么那么,对矩阵对矩阵 A-pE 应用幂法求其主特征值应用幂法求其主特征值 1-p ,收敛收敛速度将会加快速度将会加快.这种通过求这种通过求 A-pE 的主特征值和特征的主特征值和特征向量向量,进而得到进而得到A的主特征值和特征向量的方法叫的主特征值和特征向量的方法叫原点原点平移法。平移法。26ppn 1且使且使min,max112 ppppn 显然显然,当当 2 - p = - ( n- p ),即即

11、P= ( 2+ n ) 2 时时,上式取上式取最小值最小值;如果希望计算如果希望计算 n , 类似的讨论可知应选取类似的讨论可知应选取 p= ( 1 + n-1)2 。则对任意实数则对任意实数p,矩阵矩阵 A-pE 的主特征值为的主特征值为 1-p或或 n-p , 要求要求 1 , 则选则选 p 使使27例例2 用原点平移加速法求例用原点平移加速法求例1中矩阵中矩阵A的主特征值与其的主特征值与其对应的特征向量。对应的特征向量。 解解 取取p=-2.5, 做平移变换做平移变换B=A-pE，则则 5 . 36345 . 510235 . 4B 1634310232A对对B应用幂法，仍应用幂法，仍取

12、取 x0=(0,0,1)T , 则则TTvuv875. 0 , 1, 5 . 01, 4,5 . 3 , 4 , 211111), 2 , 1(),max(,0100kvuvAuvuvkkkkkkk28迭代迭代5步的计算结果见下表步的计算结果见下表0.5, 1, 0.750013.50006.7500, 13.5000, 10.125050.5, 1, 0.750013.50076.7503, 13.5007, 10.125640.5, 1, 0.750713.51796.76, 13.5179, 10.140630.5, 1, 0.7545147, 14, 10.562520.5, 1, 0

13、.8754 2, 4, 3.51k可得到可得到B的主特征值的主特征值 1 13.5000 特征向量特征向量 v1 (0.5 ,1.0, 0.7500)T 因此因此,A的主特征值为的主特征值为 1 = 1 +p 11.0000, 特征向量仍为特征向量仍为v1 =(0.5,1,0.7500)T。kvkku293031设设A为为n阶实对称矩阵，称阶实对称矩阵，称 xxxAxxR, 为为向量向量 x 的瑞利商的瑞利商，其中，其中 ( x, x)= xT x 为内积。不为内积。不难证明，对实对称矩阵难证明，对实对称矩阵A，如果其特征值满足，如果其特征值满足nn 1212、瑞利商加速、瑞利商加速由幂法公式

14、生成的由幂法公式生成的 xk 的瑞利商满足的瑞利商满足 kkkkkkxxxAxxR2121, 由此可见，由此可见，R(xk) 比比 mk 更快的收敛于更快的收敛于 1 。32幂法的瑞利商加速迭代公式为幂法的瑞利商加速迭代公式为 kkkkkkkkkkmyxkxxxymAxy/, 2 , 1,1111其中其中A为为n阶实对称矩阵。阶实对称矩阵。 对给定的误差限对给定的误差限 ，当，当 | mk mk-1 | 时，取时，取.,11kkxvm 33三、反幂法三、反幂法 反幂法是用于求非奇异矩阵反幂法是用于求非奇异矩阵A的按模最小的特征值的按模最小的特征值和对应特征向量的方法和对应特征向量的方法. 而结

15、合原点平移法的反幂法而结合原点平移法的反幂法则可以求矩阵则可以求矩阵A的任何一个具有先验了解的特征值和的任何一个具有先验了解的特征值和对应的特征向量。对应的特征向量。设矩阵设矩阵A非奇异非奇异,其特征值其特征值 i (i=1,2,n) ,满足满足0121 nn 其相应的特征向量其相应的特征向量 x1 , x2, , xn 线性无关线性无关,则则 A-1 的特的特征值为征值为1/ i ，对应的特征向量仍为，对应的特征向量仍为 xi (i=1,2, ,n).iiiiiixxAxAx1134此时此时,A-1 的特征值满足的特征值满足11111 nn因此因此,对对 A-1 应用幂法应用幂法,可求出其主

16、特征值可求出其主特征值 1/ n 和特征向量和特征向量 xn uk , k从而求得从而求得A的按模最小特征值的按模最小特征值 n 1/和对应的特征向量和对应的特征向量 xn uk ， 这种方法称为这种方法称为反幂法反幂法。k35为了避免求为了避免求 A-1 ,可通过解线性方程组可通过解线性方程组A vk= uk-1 得到得到yk ，采用采用LU分解分解,即先对即先对 A 进行进行LU分解分解 A=LU , 此时反幂法此时反幂法的迭代公式为的迭代公式为 , 2 , 1/max,1kvuvvzUvzuLzkkkkkkkkkkk求出解求出解knknux ,1), 2 , 1(),max(,01100

17、kvuvuAvuvkkkkkkk36对给定的误差对给定的误差 ，当，当 | | j 时时, 故故 , 1 ji 0lim kkB设方阵设方阵X 的的QR 的分解式为的分解式为xxRQX 由由 )()(1UDBIRQAkkxxk )(1UDRRBRIQkxxkxx 49由由 知知, 对充分大的对充分大的 k , 非奇异非奇异, 它应有唯一的它应有唯一的 QR 分解式分解式 ,并且并且 0lim kkB1 xkxRBBIkkRQ,lim,limIRIQkkkk 于是于是)(1UDRRQQAkxkkxk 但上三角阵但上三角阵 的对角线元素不一定大于零的对角线元素不一定大于零.UDRRkxk为此为此,

18、引入对角阵引入对角阵)1, 1, 1( diagDk以便保证以便保证 的对交线元素都是正数的对交线元素都是正数)(UDRRDkxkk50从而得到从而得到 的的QR 分解式分解式kA1)(1UDRRDDQQAkxkkkkxk 由由 的的 QR 分解式的唯一性得到分解式的唯一性得到kA1 UDRRDRDQQQkxkkkkkxk 从而从而 kTkkQAQA1 kxxTxTkkDQAQQQD)( 由于由于 所以所以,11TxxxxQDRRQXDXA 1 xxxTxDRRAQQ51从而从而 kkxxTkkkDQDRRQDA)(11 其中其中 nxxDRRR *2110于是于是kkTkkkDQRQDA01

19、 因为因为 为上三角阵为上三角阵, 为对角阵为对角阵,且元素为且元素为1 或或-1, 所以所以0R,limIQkk 52 );( kaikii当当 );,(0 kjiakij且且当当 的极限不一定存在的极限不一定存在 )(ijakij 53例例 1 用用 QR 算法求矩阵算法求矩阵 特征值特征值. A的特征值的特征值为为-1,4,1+2i . 2100322023011525A54解解 令令 用施密特正交化过程将用施密特正交化过程将 分解分解为为,1AA 1A 8989. 03962. 0000740. 01761. 09813. 004192. 08804. 01887. 01961. 01

20、038. 01923. 00377. 09806. 0 111RQA 7822. 00002736. 35242. 2005288. 25852. 10381. 203922. 04912. 59612. 10992. 555将将 与与 逆序相乘逆序相乘, 求出求出1Q1R2A 7031. 03099. 0001294. 38525. 04770. 205361. 15171. 29401. 13977. 02390. 05922. 19508. 56517. 4用用 代替代替A 重复上面过程重复上面过程,计算计算11次得次得1A 0000. 1*0*1211. 03290. 10*5910.

21、 38789. 10*0000. 456由由 不难看出不难看出, 矩阵矩阵A 的一个特征值是的一个特征值是4 ,另一个特另一个特 征值是征值是-1 ,其它两个特征值是方程其它两个特征值是方程 12A的根的根.求得为求得为 .21i 01211. 03920. 15910. 38789. 1 57上上Hessenberg化化58596061626364用正交变换化对称矩阵为对称三对角阵用正交变换化对称矩阵为对称三对角阵65带原点位移的带原点位移的QR算法算法6667用单步用单步QR方法计算上方法计算上Hessenberg特征值特征值68697071727374757677787980818283

22、84858687隐式隐式QR算法算法888990919293定理定理6 设设A为为n阶实对称矩阵，则必存在正交矩阵阶实对称矩阵，则必存在正交矩阵P使得使得),(21nTdiagApp其中其中n21,为为A的的n个特征值。个特征值。证证设设A的互不相同的特征值为的互不相同的特征值为,s21它们它们的重数依次为的重数依次为).(,nrrrrrrs21s21 根据性质根据性质1和性质和性质3知，知，恰有恰有ir个线性无关的实特征向量，个线性无关的实特征向量，),(s21ii 对应于特征值对应于特征值94把它们标准正交化，把它们标准正交化，特征向量组，特征向量组，特征向量共有特征向量共有n个，个，并有

23、并有 ppApp11),(n21diag其中其中n21,为为A的的n个特征值。个特征值。ir个单位正交的个单位正交的就可得到就可得到nrrrs21 知，知，由由这样的这样的特征值的特征向量是正交的，特征值的特征向量是正交的，向量两两正交，以它们为列向量构成正交矩阵向量两两正交，以它们为列向量构成正交矩阵P，又由性质又由性质2知，知， A的属于不同的属于不同故这故这n个单位特征个单位特征95由定理由定理6可知，实对称矩阵的对角化问题，实质上可知，实对称矩阵的对角化问题，实质上是求正交矩阵是求正交矩阵P的问题，计算的问题，计算P的步骤如下：的步骤如下：（1）（2）求出齐次线性方求出齐次线性方程组程

24、组0 xEAi )(的基础解系，的基础解系，进行正交化和单位化，得到进行正交化和单位化，得到A对于对于i的一组的一组标准正交的特征向量，标准正交的特征向量，的个数恰好是的个数恰好是i作为作为A的特征值的重数；的特征值的重数；;,r21求出实对称矩阵求出实对称矩阵A的全部特征值的全部特征值,i对于各个不同的特征值对于各个不同的特征值对基础解系对基础解系这个向量组所含向量这个向量组所含向量96（3）量构成一组量构成一组nR的标准正交基的标准正交基;,n21ppp（4）则则P为正交矩阵且使得为正交矩阵且使得AppApp1T 为对角阵，对角线上的元素为为对角阵，对角线上的元素为相应特征向量的特征值。相

25、应特征向量的特征值。),(r21ii 的所有标准正交的特征向的所有标准正交的特征向将将),(n21pppp 取取97例例13设实对称矩阵设实对称矩阵, 310130004A求正交阵求正交阵., Appp1使使解解.)-)(4-(2)( 20864310130004EA2 得特征值得特征值42321 ,98997.2 对称对称QR方法方法对称矩阵的三对角化对称矩阵的三对角化100带原点位移的带原点位移的QR迭代迭代101102隐式隐式QR迭代迭代103104定义定义JacobiJacobi方法是用来计算实对称矩阵方法是用来计算实对称矩阵A 的全部的全部特征值及其相应特征向量的一种变换方法特征值及

26、其相应特征向量的一种变换方法. .Jacobi方法的基本思想是通过一系列的平面旋转矩方法的基本思想是通过一系列的平面旋转矩阵构成的正交变换将对称矩阵逐步化为对角阵阵构成的正交变换将对称矩阵逐步化为对角阵,从从而得到而得到A的全部特征值及其相应的特征向量的全部特征值及其相应的特征向量.定理定理 (1)(1)如果如果n 阶方阵阶方阵A满足满足 则称则称A为正交阵为正交阵. . ）TTAAIAA 1(即7.3 Jacobi7.3 Jacobi方法方法105(2) 设设A 是是n 阶实对称矩阵阶实对称矩阵, ,则则A 的特征值都是的特征值都是 实数实数, ,并且有互相正交的并且有互相正交的n 个特征向

27、量个特征向量. .(3) 相似矩阵具有相同的特征值相似矩阵具有相同的特征值. . (4) 设设A 是是n 阶实对称矩阵阶实对称矩阵, ,P 为为n 阶正交阵阶正交阵, ,则则 B=PAP 也是对称矩阵也是对称矩阵. .(5) n 阶正交矩阵的乘积是正交矩阵。阶正交矩阵的乘积是正交矩阵。(6)(6)设设A 是是n n 阶实对称矩阵阶实对称矩阵, ,则必有正交矩阵则必有正交矩阵P ， 使使 nTAPP 21106由由(6) 可知可知,对于任意的对于任意的 阶实对称矩阵阶实对称矩阵A , 只只要能求得一个要能求得一个 正交阵正交阵P ,使使 阵阵为为对对 角角 APPT则可得到则可得到A 的全部特征

28、值及其相应的特征向量的全部特征值及其相应的特征向量,这就是雅这就是雅克比方法的理论基础克比方法的理论基础.其中其中 的对角线元素是的对角线元素是A 的的n 个个特征值，正交特征值，正交矩阵矩阵P 的第的第i 列是列是A 的对应于特征值的对应于特征值 的特征的特征向量。向量。 i 下面我们详细介绍雅可比方法。首先引进下面我们详细介绍雅可比方法。首先引进 中中 的平面旋转变换的平面旋转变换.变换变换 jikyxyyxyyxkkiijiii,cossinsincos nR107 11cossinsincos1 ijP i j i j记为记为 ,其中其中yPxij 108X= , Tnxxx).(,2

29、,1Y = ,Tnyyy).(,2,1jixx, 称为称为 平面内的平面旋转矩阵平面内的平面旋转矩阵.容易得到如容易得到如 下下性质性质:ijpnRypijjixx, 则称则称 x= 为为 中中 平面内的一个旋转变换平面内的一个旋转变换, (2) 的主对角线元素中除第的主对角线元素中除第i 个与第个与第j 个元素个元素为为 外外,其它元素均为其它元素均为1;非对角元素中除第非对角元素中除第i 行第行第 j 列元素为列元素为 ,第第j 行第行第i 列元素为列元素为 外外,其它元素均为零其它元素均为零. ijp cos sin sin(1) 为正交矩阵为正交矩阵ijP109(3) 只改变只改变A的

30、第的第i行第行第j 行元素行元素, AP只改变只改变A 的的 第第i 列第列第j 列元素列元素,所以所以 只改变只改变A 的第的第i 行行, 第第j 行行,第第i 列列,第第j 列元素列元素. APTAPPT 设设A= 为为n 阶实对称矩阵阶实对称矩阵, 为一非对角线元素为一非对角线元素.令令 则则 为实对为实对 称矩阵称矩阵,且且 与与A 有同的特征值有同的特征值. 通过直接计算知通过直接计算知 3 nannij0 jiijaa nnijijTijaAPPA 111A1A110),(jik ),(jik ),(cossinsincos2cos2sin)(212sincossin2sinsin

31、cos) 1() 1() 1() 1() 1() 1() 1(22) 1(22) 1(jilkaaaaaaaaaaaaaaaaaaaaaaaklkljkjkjkjkjkikkiikijiijjjiijijjjiijjijjjiiii 111当取当取 满足关系式满足关系式 时时, 且且 )10(22tanjjiiijaaa , 011 jiijaa )11(),(2),(2212222121222121 jilkaaaaaaajikaaaaklklijjjiijjiijkikjkik由于在正交相似变换下由于在正交相似变换下, 矩阵元素的平方和不变矩阵元素的平方和不变,所以若所以若用用D(A) 表

32、示矩阵表示矩阵A的对角线元素的平方和的对角线元素的平方和,用用S(A) 表示表示A 的非对角元素平方和的非对角元素平方和,则由则由(11)式得式得112)12(2)()(2)()(2121 ijijaASASaADAD1A , , 的非对角元素平方和的非对角元素平方和A的非对角元素平方的非对角元素平方 且将事先选定的非对角元素消去了且将事先选定的非对角元素消去了 22ija1A这说明用这说明用 对对A 作正交相似变换化为作正交相似变换化为 后后, , 的的 对角线元素平方和比对角线元素平方和比A 的对角元素平方和增加了的对角元素平方和增加了 ijP1A 01 ija即即,22ija和减少了和减

33、少了113 因此因此,只要我们逐次地用这种变换只要我们逐次地用这种变换,就可以使得就可以使得 矩阵矩阵A 的非对角线元素平方和趋于零的非对角线元素平方和趋于零,也即使得也即使得矩阵矩阵A 逐步化为对角阵逐步化为对角阵. 这里需要说明一点这里需要说明一点:并不是对矩阵并不是对矩阵A的每一对非的每一对非对角线非零元素进行一次这样的变换就能得到对对角线非零元素进行一次这样的变换就能得到对角阵角阵.因为在用变换消去因为在用变换消去 的时候的时候,只有第只有第 i 行行,第第 j 行行,第第 i 列列,第第 j 列元素在变化列元素在变化,如果如果 或或 为零为零,经变换后又往往不是零了经变换后又往往不是

34、零了.ijaikajka114 雅克比方法就是逐步对矩阵雅克比方法就是逐步对矩阵A 进行正交相似变进行正交相似变换换,消去非对角线上的非零元素消去非对角线上的非零元素,直到将直到将A的非对角的非对角线元素化为接近与零为止线元素化为接近与零为止,从而求得从而求得A的去全特的去全特征征值值,把逐次的正交相似变换矩阵乘起来把逐次的正交相似变换矩阵乘起来,便是所要便是所要求的特征向量求的特征向量.Jacobi的计算步骤归纳如下的计算步骤归纳如下:第一步第一步 在矩阵在矩阵A 的非对角元素中选取一个非零的非对角元素中选取一个非零 元素元素 .一般来说一般来说,取绝对值最大的非对取绝对值最大的非对 角线元

35、素角线元素. ija115第二步第二步 由公式由公式 求出求出 , 从而得从而得 平面旋转矩阵平面旋转矩阵jjiiijaaa 2tan ijPP 1第三步第三步 , 的元素由公式的元素由公式(9)计算计算111APPAT 1A2A第四步第四步 以以 代替代替A ,重复第一重复第一,二二,三步求出三步求出 及及 继续重复这一过程继续重复这一过程,直到直到 的非对交的非对交 线元素全化为充分小线元素全化为充分小(即小于允许误差即小于允许误差) 时为止时为止. 1A2PmA第五步第五步 的对角化元素为的对角化元素为A的全部特征值的近的全部特征值的近 似值似值, 的第的第j 列为对应于特列为对应于特

36、征值征值 ( 为为 的对角线上第的对角线上第j 个元素个元素) 的特征向量的特征向量.mAmPPPP21 j j mA116例例 1 1 用雅克比方法求矩阵用雅克比方法求矩阵 的特征值与特征向量的特征值与特征向量 210121012 A解解 首先首先 i = 1 , j =2 , 由于由于 故故 取取 所以所以 22211 aa,4 1000212102121 121PP117 2212121302101 111APPAT再取再取i =1,=1,j =3,=3,由由 2212122tan 得得88808. 0cos,45969. 0sin 所以所以 88808. 0045969. 001045

37、969. 0088808. 0 2P118 36603. 262797. 0062797. 0332505. 0032505. 063398. 0 222APPAT继续做下去继续做下去,直到非对角线元素趋于零直到非对角线元素趋于零,进行九次变进行九次变换后换后,得得 41421. 300000. 000000. 000000. 000000. 200000. 000000. 000000. 058758. 0 9A 的对角线元素就是的对角线元素就是A 的特征值的特征值,即即 9A41421. 3,00000. 2,58758. 0321 119 36603. 262797. 0062797.

38、0332505. 0032505. 063398. 0 222APPAT继续做下去继续做下去,直到非对角线元素趋于零直到非对角线元素趋于零,进行九次变进行九次变换后换后,得得 41421. 300000. 000000. 000000. 000000. 200000. 000000. 000000. 058758. 0 9A 的对角线元素就是的对角线元素就是A 的特征值的特征值,即即 9A41421. 3,00000. 2,58758. 0321 120 用雅可比方法求得的结果都比较高用雅可比方法求得的结果都比较高,特别是求特别是求得的特征向量正交性很好得的特征向量正交性很好,所以雅可比方法是

39、求实所以雅可比方法是求实对称矩阵的全部特征值及其对应特征向量的一个对称矩阵的全部特征值及其对应特征向量的一个较好的方法较好的方法. 但由于上面介绍的雅可比方法每次迭代都选取但由于上面介绍的雅可比方法每次迭代都选取绝对值最大的非对角元素作为消去对象绝对值最大的非对角元素作为消去对象,花费很多机花费很多机器时间器时间.另外当矩阵是稀疏矩阵时另外当矩阵是稀疏矩阵时,进行正交相似变进行正交相似变换后并不能保证其稀疏的性质换后并不能保证其稀疏的性质,所以对阶数较高的矩所以对阶数较高的矩阵不一采用这种方法阵不一采用这种方法.121 目前常采用一种过关雅克比方法目前常采用一种过关雅克比方法.这种方法是这种方法是选取一个单调减小而趋于零的数列选取一个单调减小而趋于零的数列 作为极值作为极值,这些这些极称为极称为“关关”,对矩阵的非对角线元素规定一个顺对矩阵的非对角线元素规定一个顺序序(例如先行后列例如先行后列,自左至右的顺序自左至右的顺序). 首先对限值首先对限值 按规定的顺序逐个检查矩阵的非按规定的顺序逐个检查矩阵的非对角线元素对角线元素,碰到绝对值小于碰到绝对值小于 的元素就跳过去的元素就跳过去,否否则就作变换将其化为零则就作变换将其化为零.1a1a重复上述过程重复上述过程,直到所有的非对角元素的绝对值都直到所有的非对角元素的绝对值都小于小于 为止为止.1a