2-10 例:用幂法求矩阵A=02-1的按模 0-12 最大的特征值和相应的特征向量 取x)=(0,0,1),E≤10-3 解:y(0)=x0)=(0,0,1) (0,-1,2),=2, =(O,-0.5,1) (0.5,-2,2.5),O=2.5
(0) 3 (0) (0) (1) (0) (1) (1) (2) (1) 2 1 0 0 2 1 0 1 2 (0,0,1) , 10 . (0,0,1) , (0, 1, 2) , 2, (0, 0.5,1) , (0.5, 2, 2.5) , T T T T T A x y x x Ay x y x Ay − − = − − = = = = = − = = = − = = − = 例:用幂法求矩阵 的按模 最大的特征值和相应的特征向量。 取 解: 2.5
y=(27650948,-2.9981848,2.9999024) 0=2.9990924 y8)=(0,9219772,-0996973,1) y=(28436517,-29993946,29996973) 由29996973-29990924=00006049<10 故λ1≈2996973.相应特征向量为 l(2.8436517,-29993946,2.9996973) 事实上,的特征值=3,2=2,4=1 与对应的特征向量为(1,-1,1) 此例中比值为
(8) (7) (8) (9) (8) 3 1 (2.7650948, 2.9981848,2.9990924) 2.9990924 (0.9219772, 0.9996973,1) (2.8436517, 2.9993946,2.9996973). 2.9996973 2.9990924 0.0006049 10 . 2.9996973. x Ay y x Ay − = = − = = − = = − − = 由 故 相应特 1 2 3 1 2 1 (2.8436517, 2.9993946,2.9996973) 3, 2, 1 1 -1,1 2 . 3 T u A − = = = = 征向量为 。 事实上, 的特征值 , 与 对应的特征向量为(, )。 此例中比值为
两种特殊情况 前面假定|λ|>λ如果按模最大的特征值有多个 甲4|=2|=…=|2m}>12m≥…≥|2 幂法是否有效? (1)是m重根,即1=22=…=n,矩阵A仍有n个 线性无关的特征向量。此时有 (+y) x[1l1+…+cn2l m+1 k+ k+1 +1m+1 ·· 显然,只要a12…,an不全为零,当k充分大时,就有 x(6+1)≈A4+(x1l1+…+anln)
两种特殊情况 1 2 1 2 1 1 1 2 ( 1) 1 1 1 1 1 1 1 1 1 1 1 . 1 , [ ( ) ( ) ] m m n m k k m m m n k k m m n n m A n x u u u u + + + + + + + + = = = = = = = + + + + + 前面假定 如果按模最大的特征值有多个, 即 幂法是否有效? ( ) 是 重根,即 矩阵 仍有 个 线性无关的特征向量。此时有 显然 1 ( 1) 1 1 1 1 , , ( ) m k k m m k x u u + + + + ,只要 不全为零,当 充分大时,就有
a1+…+ann.是矩阵A相应于的特征向量 Aa, u=Ma, Aa,u2=na,u2,., Aamum=n,mamu →A(14+a21l2+…+ann)=41al1+2a2l2+…+nann 1(a1 +CL2+…+C,l m m 因α+…+αnu也是矩阵A相应于λ的特征向量,故有 m x(6)为相应的特征向量,即对这种情况幂法仍然有效
1 1 1 1 1 1 1 1 2 2 2 2 2 1 1 2 12 1 1 1 2 2 2 1 1 1 2 2 1 1 1 , , , , ( ) ( ) m m m m m m m m m m m m m m m m u u A A u u A u u A u u A u u u u u u u u u u u A + + = = = + + + = + + + = + + + + + 也是矩阵 相应于 的特征向量 因 也是矩阵 相应于 的特征向量,故有 ( 1) 1 2 ( ) ( 1) k i m k i k x x x + + = = = 为相应的特征向量,即对这种情况幂法仍然有效
(2)=-12,4>2,且矩阵A有n个线性无关的特征向量。 4+=2H[au4+(-1+a22+ C22++ 由上式可知,x是个摆动序列,当k充分大时,有 x(k-1) (ax1-a2l2)x26)≈2(al4a2l2) (k+2) (k+2)/(k 又由x4+1)≈4[ax4+(-1)a2l2 (k) n[c1+(-1)a22 +1)+21x()≈2+a1 41)-1x()22+(-1)+a2u2 故在这种情况下,仍可按幂法产生向量序列
1 2 1 3 ( 1) 1 1 1 1 3 1 1 1 2 2 3 3 1 1 ( 1) (2 1) 2 1 (2 ) 2 1 1 1 2 2 1 1 1 2 2 ( 2) 2 ( 1 1,2 ( ) 2 , , [ ( 1) ( ) ( ) ] ( ) ( ) k k k k k n n n k k k k k k k i k i i A n x u u u u x k x u u x u u x x x + + + + + + − − + = − = + − + + + − + ( ) 且矩阵 有 个线性无关的特征向量。 由上式可知, 是个摆动序列,当 充分大时,有 2) ( ) ( 1) 1 1 1 1 1 2 2 ( ) 1 1 1 2 2 ( 1) ( ) 1 1 1 1 1 ( 1) ( ) 1 1 1 1 2 2 / [ ( 1) ] [ ( 1) ] 2 2 ( 1) k i k k k k k k k k k k k k k x x u u x u u x x u x x u + + + + + + + + + + − + − + − − 又由 故在这种情况下,仍可按幂法产生向量序列