幂法 反幂法 matlab,幂法反幂法求解矩阵最大最小特征值及其对应的特征向量

数值计算解矩阵的按模最大最小特征值及对应的特征向量

一.幂法

1. 幂法简介：

当矩阵A 满足一定条件时，在工程中可用幂法计算其主特征值(按模最大)及其特征向量。矩阵A 需要满足的条件为： (1) 的特征值为A i n λλλλ,0||...||||21

≥≥≥>

(2) 存在n 个线性无关的特征向量，设为n x x x ,...,,21

1.1计算过程：

i n

i i i u x

x αα,1

)

0()0(∑==，有对任意向量不全为0，则有

1

11111221

12111

1

1

11

1

011)()(...u u a u a u λu λαu αA x A Ax x k n n k n k k n

i i

k i i n

i i i k )(k (k))(k αλλλλλα++++=+=+++≈?

?

????+++======∑∑Λ 可见，当||1

2

λλ越小时，收敛越快；且当k 充分大时，有1)1111)11111λαλαλ=???

???==+++(k )(k k

(k k )(k x x u x u x ，对应的特征向量即是)(k x 1+。

2 算法实现

.

,, 3,,1 , ).5()

5(,,,,||).4();max(,).3()

(max(;0,1).2(,).1()()

()

(停机否则输出失败信息转置若转否则输出若计算最大迭代次数，误差限，初始向量输入矩阵βλβεβλβλε←+←<

←←k k N k y x Ay x x abs x y k N x A k k k

3 matlab 程序代码