( )
§1 幂法和逆幂法(求部分eigenvalues and eigenveltors)
一. 幂法 记λ (A ) = {λ 1, λ2,L, λn }满足 : λi ui = Aui , i = 1,2, L , n.
设 λ1 ≥ λ2 ≥ L ≥ λn (1), 且u1 , u 2 , L , u n 线性无关.任取初始向量 x0 ≠ 0. 由xk +1 = Au xk , k = 0,1, L , (2 )得{xk }.由x0 ≠ 0知∃不全为0的数α1 , α 2 , L , α n , 使 x0 = α1 u1 + α 2 u 2 + L + α n u n (3)
Τ
9 0 用平移原点的方法,p = − 1, B = A − pI = A + I = 18 1 18 1
k
y = xk k xk x = Ay k k +1
∞
, k = 0,1, L
列表计算 :
0
1
1 .0 0 0 2 .0 0 0 2 .0 0 0
只要 α 1 ≠ 0, 则当k → ∞ 时
而xk +1 = Axk
∴ Axk
λ1xk 即λ1为近似的特征值.
xk 就是近似的特征向量.
结论:由xk+1=Axk ( k=0,1, ) 得{xk} , 若序列中相邻二向量的对应分量 L uu r 成比例,即xk+1 λxk ,则λ, k为A的近似特征值和相应的特征向量 x
T
(矩阵实对称时,如何加速收敛)
(
)
x∗ Ax 称 ∗ 为A的关于的瑞利商. xx
Q AT = A
ur uu r ∴ λ ( A ) ∈ R, 且A有n个两两相交的特征向量 : u1 , u2 ,L , uu ur uu uu r r r urT uu 0, i ≠ j r 并设un .u1 , u2 ,L un已被标准化, 即 ui , u j = 1, i = j