矩阵特征值计算.pptxVIP

第8章矩阵特征问题的计算8.1引言8.2幂法及反幂法8.3豪斯霍尔德方法8.4QR方法

工程技术中有多种振动问题，如桥梁或建筑物的振动，机械零件、飞机机翼的振动，及一些稳定性分析和相关分析在数学上都可转化为求矩阵特征值与特征向量的问题.下面先复习一些矩阵的特征值和特征向量的基础知识.8.1引言

定义1⑴已知n阶矩阵A=(aij)，则称为A的特征多项式.01一般有n个根(实的或复的，复根按重数计算)称为A的特征值.用λ(A)表示A的所有特征值的集合.02A的特征方程03

设λ为A的特征值，相应的齐次方程组注：当A为实矩阵时，?(λ)=0为实系数n次代数方程，其复根是共轭成对出现.的非零解x称为矩阵A的对应于λ的特征向量.例1求A的特征值及特征向量，其中

求得矩阵A的特征值为：对应于各特征值矩阵A的特征向量分别为：解矩阵A的特征方程为

定理1设λ为A∈Rn×n的特征值,且Ax=λx(x?0)，则有01cλ为的cA特征值(c≠0为常数)；03设A为非奇异矩阵，那么λ≠0,且λ-1为A-1的特征值，即A-1x=λ-1x.05λ-p为A-pI的特征值，即(A-pI)x=(λ-p)x；02λk为Ak的特征值，即Akx=λkx；04下面叙述有关特征值的一些结论：

定理2设λi(i=1,2,?,n)为n阶矩阵A=(aij)的特征值，则有⑴称为A的迹；⑵定理3设A∈Rn×n，则有定理4设A为分块上三角矩阵，即其中每个对角块Aii均为方阵，则

1定理5设A与B为相似矩阵（即存在非奇异矩阵P使B=P-1AP），则2定理5说明，一个矩阵A经过相似变换，其特征值不变.3一个亏损矩阵是一个没有足够特征向量的矩阵，亏损矩阵在理论上和计算上都存在困难.4A与B有相同的特征值；5如果y是B的特征向量，则Py是A的特征向量.6定义2如果实矩阵A有一个重数为k的特征值λ,且对应于λ的A的线性无关的特征向量个数k，则A称为亏损矩阵.

定理6⑴A∈Rn×n可对角化，即存在非奇异矩阵P使的充分必要条件是A具有n个线性无关的特征向量.如果A∈Rn×n有m个(m≤n)不同的特征值λ1,λ2,?,λm，则对应的特征向量x1,x2,?,xm线性无关.

A有n个线性无关的特征向量；A的特征值均为实数；且λ1,λ2,?,λn为A的特征值，而P＝(u1,u2,?,un)列向量uj为A的对应于λj的单位特征向量.存在一个正交矩阵P使的定理7(对称矩阵的正交约化)设A∈Rn×n为对称矩阵，则

下面讨论矩阵特征值界的估计.定义3设n阶矩阵A=(aij)，令⑴；⑵集合称为复平面上以aii为圆心，以ri为半径的n阶矩阵A的n个Gerschgorin圆盘.

定理8(Gerschgorin圆盘定理)特别地，如果A的一个圆盘Di是与其它圆盘分离(即孤立圆盘)，则Di中精确地包含A的一个特征值.设n阶矩阵A＝(aij)，则A的每一个特征值必属于下面某个圆盘之中如果A有m个圆盘组成一个连通的并集S，且S与余下n-m个圆盘是分离的，则S内恰包含A的m个特征值.或者说A的特征值都在n个圆盘的并集中.

证明只就⑴给出证明.设λ为A的特征值，即Ax=λx，其中x=(x1,x2,?,xn)T?0.或记，考虑Ax=λx的第k个方程，即于是即

这说明，A的每一个特征值必位于A的一个圆盘中，并且相应的特征值λ一定位于第k个圆盘中(其中k是对应特征向量x绝对值最大的分量的下标).利用相似矩阵性质，有时可以获得A的特征值进一步的估计，即适当选取非奇异对角阵并做相似变换.适当选取可使某些圆盘半径及连通性发生变化.

1例2估计矩阵A的特征值范围，其中2解矩阵A的3个圆盘为3由定理8，可知A的3个特征值位于3个圆盘的并集中，由于D1是孤立圆盘，所以D1内恰好包含A的一个特征值λ1(为实特征值)，即4A的其它两个特征值λ2,λ3包含在D