第7章 求矩阵特征值的数值方法

⎤
x
j
⎥ ⎥⎦
→ a1x1 + a2 x2 (k → ∞)
max(a1x1 + a2 x2 )
在 a1 + a2 ≠ 0 时，这是一个主特征向量。
18
② 当 λ1 = −λ2 时， λ 2 是 A2 的主特征值，故
λ1,2
=
± lim k →∞
max( A2zk ) ．
计算时，迭代可能发散，但是
A2 zk = Ayk+1 = mk+1 Azk+1 = mk+1mk+2 zk+2
⎜ ⎝
⎡
n
max ⎢⎢⎣a1x1
+ a2 x2
+
aj
j =3
λj λ1
⎞k −1 ⎟ ⎠
x
⎛ ⎜ ⎝
λj λ1
⎞k −1 ⎟ ⎠
⎤
j
⎥ ⎥⎦
xj
⎤ ⎥ ⎥⎦
⎫ ⎪ ⎪⎪ ⎬ ⎪ ⎪ ⎭⎪
，
∑∑ ( ) =
λ1
⎡
m
ax
⎢ ⎢⎣
a1
x1
+
a2
x
2
+
n j=2
a
j
⎛ ⎜ ⎝
λj λ1
⎞k ⎟ ⎠
⎤
x
j
⎥ ⎥⎦
+
n
aj
j=2
⎛ ⎜ ⎝
λj λ1
⎞k −1 ⎟ ⎠
⎤
x
j
⎥ ⎥⎦
→
λ1
k → ∞ （主特征值），
15
∑∑ ( ) zk
=
a1 x1
+
n j=2
aj
⎛ ⎜ ⎝
λj λ1
⎞k ⎟ ⎠
⎡ max ⎢⎢⎣a1x1
+
n
aj
j=2
⎛ ⎜ ⎝
λj λ1
xj
⎞k ⎟ ⎠
xj
⎤ ⎥ ⎥⎦
→
x1 max( x1)
k → ∞, a1 ≠ 0
⎡
m
ax
⎢ ⎢⎣
a1
x1
+ a2 x2
+
n
aj
j=2
⎛ ⎜ ⎝
λj λ1
⎞ k −1 ⎟ ⎠
⎤
x
j
⎥ ⎥⎦
→
λ1
k→∞
17
∑∑ zk
=
a1 x1
+ a2 x2
+
n
aj
j =3
⎛ ⎜ ⎝
λj λ1
⎞k ⎟ ⎠
⎡ max ⎢⎢⎣a1 x1
+ a2 x2
+
n
aj
j=3
⎛ ⎜ ⎝
λj λ1
xj
⎞k ⎟ ⎠
也是属于 λ 的特征向量。
证明： A(a1x1 + a2 x2 ) = a1Ax1 + a2 Ax2 = λ (a1x1 + a2 x2 ) ，
所以． a1x1 + a2 x2 是属于 λ 的特征向量。
8
2．幂法
幂法是计算实矩阵按模最大的特征值及其对应特征 向量的一种迭代方法，主要适用于中小型矩阵和大型稀疏 矩阵。
k = 1, 2," ，
26
迭代计算得
k
0
yk
mk
1
zk
1
1
1 0.875 -0.25 0.0417
2 2.0417 -1.298 0.3512
3 2.4071 -1.679 0.4834
4 2.4725 -1.748 0.5076
5 2.4825 -1.758 0.5113
0.875 2.0417 2.4071 2.4725 2.4825
则称 A 为正定矩阵。
2
定理 1． 设 A = (aij ) ∈ R n×n ，则下列条件等价 （1） A 为正定矩阵； （2） A 的所有特征值都是正数； （3） A 的各阶顺序主子式均大于零。
3
定义 3． 设 A = (aij ) ∈ C n×n ，且 AT A = I ，则称 A 为正交矩阵。
5
定理 3． λ 为方阵 A 的特征值，则 λ 是方程 λI − A =0
的根。
注．显然，若 λ 为方阵 A 的特征值，则存在向量 x ≠ 0 ，使得：
Ax = λ x ，即齐次线性方程组 (λ I − A) x = 0 有非零解，所以系
数行列式 λ I − A = 0 。属于 λ 的特征向量是 (λ I − A) x = 0 的
第七章 求矩阵特征值的数值方法
1
1、预备知识
定义 1 设 A = (aij )n×n , 如果 AT = A ，则称 A 为对称矩阵。 定义 2 设 A = (aij ) ∈ R n×n 是对称矩阵，且对 ∀x ∈ Rn , x ≠ 0 ，
都有
n
∑ xT Ax = aij xi x j > 0 ， i, j=1
1 -0.286 0.0476
1 -0.636 0.172
1 -0.697 0.2008
1 -0.707 0.2053
1 -0.708 0.2059
27
所以，特征值
λ3
=
1 2.4825
=
0.4028
，
相应特征向量 x = (1, −0.708, −0.2059)T 。
28
4、原点位移技术
幂法的收敛速度取决于收敛因子即次大特征值与 最大特征值之比，当收敛因子越小，收敛越快。由于 一般情况下的收敛因子不是很小，所以幂法的收敛速 度是很慢的。如果经若干次幂法迭代后得到特征值的 一个粗糙近似值a和对应的粗糙特征向量x，那么取初 始位移为a,以x为初始向量，作动态原点位移的反幂 法，就可以加速收敛性。
非零解．
6
定理 4． λ 为方阵 A 的特征值， P (t ) 是一个多项式，则 P (λ ) 是 P ( A) 的一个特征值。 证明：设 P (t ) = amt m + am−1t m−1 + ...... + a1t + a0 ，x 是属于 λ 的
特征向量。则
( ) P A = am Am + am−1Am−1 + ...... + a1A + a0 I ， P ( A) x = am Am x + am−1Am−1x + ...... + a1Ax + a0 x
为了避免求 A−1 ，通常利用解线性方程组方法。其迭代格式
( ) 为，取初始向量 z0 z0 ∞ = 1 ， 解线性方程组得yk
⎧⎪⎨mAyk k==mzakx−1( yk ) k = 1, 2,"
⎪⎩zk = yk / mk
24
当 λ1 ≥ λ2 ≥ ... ≥ λn−1 > λn 时，有
mk
→
+ a2 x2
+
n
aj
j=3
⎛ ⎜ ⎝
λj λ1
⎞k −1 ⎟ ⎠
xj
⎤ ⎥ ⎥⎦
⎡ max ⎢⎢⎣a1x1
+ a2 x2
+
n
aj
j=3
⎛ ⎜ ⎝
λj λ1
⎞k −1 ⎟ ⎠
xj
⎤ ⎥ ⎥⎦
可得：
Azk
+
λ1zk
→
2λ1a1 x1
max (a1x1 + a2 x2 )
(k
→
∞)
Azk
− λ1zk
→
−2λ1a2 x2
max (a1x1 + a2 x2 )
(k
→
∞)
（属于 λ1 和 λ2 = −λ1 的特征向量）
21
例
1．
用幂法求矩阵
A
=
⎜⎛ ⎜
2 3
4 9
6 15
⎟⎞ ⎟
的按模最大的特征
⎜⎝ 4 16 36⎟⎠
值和对应的特征向量。
解：取初始向量 z0 = (1,1,1)T ，由
⎧⎪⎨mykk
= Azk−1 = max(
4
定理 2．关于正交矩阵有如下结论
（1）单位矩阵是正交矩阵。
（2）若 A 为正交矩阵，则 A−1 = AT 。 （3）若 A 为正交矩阵，则 AT 也是正交矩阵。 （4）若 A 为正交矩阵，则 A = 1 或 A = −1 。 （5）若 A, B 为同阶正交矩阵，则 AB 与 BA 也是正交矩阵。
= amλ m x + am−1λ m−1x + ...... + a1λ x + a0 x = P (λ ) x
即， P (λ ) 是 P ( A) 的一个特征值
7
定理 5． λ 为方阵 A 的特征值，若 x1, x2 都是属于 λ 的特
征向量，则
( ) a1x1 + a2 x2 a1 + a2 ≠ 0
又设 z0 = a1x1 + a2 x2 +" + an xn ，
Ax j = λj x j ( j = 1, 2,", n) ，
∑ 因此有： zk
=
λ1k
⎡ ⎢a1 x1 ⎣
+
j
n =2
a
j
(
λj λ1
)k
⎤ xj ⎥
⎦
。
所以, zk 渐近于特征方向,但是 zk 可能无限增长或收敛于零。
10
⎜⎛ 2 0 1⎟⎞ 例如，矩阵 A = ⎜ 1 1 0⎟ ，取 z0 = (1, 0,1)T ，利用上述方法