特征值数值解法
- 格式:doc
- 大小:57.00 KB
- 文档页数:2
第八章 矩阵特征值问题的数值解法
矩阵特征值总是有广泛的应用背景. 例如在科学技术领域中,动力系统和结构系统中的
振动问题、电力系统的静态稳定分析上、工程设计中的某些临界值的确定等,都归结为矩阵
特征值问题. 本章介绍n阶实矩阵nnRA的特征值与特征向量的求解方法,即求参数λ和
相应的非零向量x,使Ax=λx,即(A-λI)x=0,并称λ为A的特征值,x为相应于λ的特征
向量.而
0)(xIA
有非零解的充分必要条件是
,0)(det)(111nnnnaaaIA
其中
),,2,1(nia
i
为常数. 由于上面方程是λ的n次多项式,因此它有n个根(实根或复
根). 除特殊情况外(如n=2,3或A为上(下)三角矩阵,一般不直接求解,原因是这样的算法往
往不稳定.在计算上常用的方法是乘幂法与反幂法(迭代法)和相似变换方法(变换法). 本章
只介绍求矩阵特征值与特征向量的这两种基本方法.
第一节 乘幂法及反幂法
一、乘幂法
设矩阵nnRA的n个特征值满足
0n321
(1.1)
且有相应的n个线性无关的特征向量.,,21nxxx乘幂法是计算矩阵按模最大特征值及相应
特征向量的迭代法,其基本思想是对任给的非零向量,0nRz用矩阵A连续左乘,构造迭
代过程,具体过程是:
由假设知niiixz110),0(用A左乘两边得
niiiiniiixAxAzz11
01
.
再用A左乘上式,得
niiixzAAzz1
2102
12
.
一直这样做下去,一般地有
).,2,1( 111111101kxx
xzAAzz
n
i
i
k
i
i
k
n
i
ikikkk
我们只讨论21的情况,对其他情况的讨论可根据参考文献[2]参阅有关资料.由(1.2)知
,lim111xzkkk
(1.3)
于是对充分大的k有
.111xzkk
(1.4)
(1.3)表明序列kkz1越来直接近A的相应于1的特征向理(11,0x是A的相应于1的特
征向量的近似向量,其收敛速度取决于比值12.
下面我们来计算1. 由于
,1111111011ikiniikkkkxxzAAzz
当k充分大时, 11111xzkk,于是可知