特征值数值解法

  • 格式:doc
  • 大小:57.00 KB
  • 文档页数:2

第八章 矩阵特征值问题的数值解法
矩阵特征值总是有广泛的应用背景. 例如在科学技术领域中,动力系统和结构系统中的
振动问题、电力系统的静态稳定分析上、工程设计中的某些临界值的确定等,都归结为矩阵

特征值问题. 本章介绍n阶实矩阵nnRA的特征值与特征向量的求解方法,即求参数λ和
相应的非零向量x,使Ax=λx,即(A-λI)x=0,并称λ为A的特征值,x为相应于λ的特征
向量.而

0)(xIA

有非零解的充分必要条件是
,0)(det)(111nnnnaaaIA

其中
),,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( 111111101kxx
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. 由于
,1111111011ikiniikkkkxxzAAzz




当k充分大时, 11111xzkk,于是可知