《数值分析ppt第8章矩阵特征值问题计算》由会员分享，可在线阅读，更多相关《数值分析ppt第8章矩阵特征值问题计算（84页珍藏版）》请在金锄头文库上搜索。
1、第8章 矩阵特征问题的计算,8.1 引言 8.2 幂法及反幂法 8.3 豪斯霍尔德方法 8.4 QR方法,8.1 引 言,工程技术中有多种振动问题，如桥梁或建筑物的振动，机械零件、飞机机翼的振动，及一些稳定性分析和相关分析在数学上都可转化为求矩阵特征值与特征向量的问题.,下面先复习一些矩阵的特征值和特征向量的基础知识.,定义1 已知n阶矩阵A=(aij)，则,称为A的特征多项式.,一般有n个根(实的或复的，复根按重数计算)称为A的特征值. 用(A)表示A的所有特征值的集合.,A的特征方程, 设为A的特征值，相应的齐次方程组,注：当A为实矩阵时， ()=0为实系数n次代数方程，其复根是共轭成对出
2、现.,的非零解x称为矩阵A的对应于的特征向量.,例1 求A的特征值及特征向量，其中,解 矩阵A的特征方程为,求得矩阵A的特征值为：,对应于各特征值矩阵A的特征向量分别为：,定理1 设为ARnn的特征值, 且Ax=x (x0)，则有, -p为A-pI的特征值，即(A-pI)x=(-p)x ；, c为的cA特征值(c0为常数)；,下面叙述有关特征值的一些结论：, k为Ak的特征值，即Akx=kx ；, 设A为非奇异矩阵，那么0 , 且-1为A-1的特征值，即A-1x=-1x .,定理2 设i(i=1,2,n)为n阶矩阵A=(aij)的特征值，则有, 称为A的迹；,定理3 设ARnn，则有,定理4
3、设A 为分块上三角矩阵，即,其中每个对角块Aii均为方阵，则,定理5 设A与B为相似矩阵（即存在非奇异矩阵P使B=P-1AP），则,定理5说明，一个矩阵A经过相似变换，其特征值不变.,一个亏损矩阵是一个没有足够特征向量的矩阵，亏损矩阵在理论上和计算上都存在困难., A与B有相同的特征值；, 如果y是B的特征向量，则Py是A的特征向量.,定义2 如果实矩阵A有一个重数为k的特征值, 且对应于的A的线性无关的特征向量个数 k，则A称为亏损矩阵.,定理6 ARnn可对角化，即存在非奇异矩阵P使,的充分必要条件是A具有n个线性无关的特征向量., 如果ARnn有 m个 (mn) 不同的特征值1,2,m，
4、则对应的特征向量 x1,x2, xm 线性无关.,定理7(对称矩阵的正交约化) 设ARnn为对称矩阵，则, 存在一个正交矩阵P使的,且1,2,n为A的特征值，而P(u1,u2,un) 列向量uj为A的对应于j 的单位特征向量., A的特征值均为实数；, A有n个线性无关的特征向量；,定义3 设n阶矩阵A=(aij)，令,下面讨论矩阵特征值界的估计., ；, 集合 称为复平面上以aii为圆心，以ri为半径的n阶矩阵A的n个Gerschgorin圆盘.,定理8 (Gerschgorin圆盘定理),特别地，如果A的一个圆盘Di是与其它圆盘分离(即孤立圆盘)，则Di中精确地包含A的一个特征值., 设n
5、阶矩阵A(aij)，则A的每一个特征值必属于下面某个圆盘之中, 如果A有m个圆盘组成一个连通的并集S，且S与余下n-m个圆盘是分离的，则S内恰包含A的m个特征值.,或者说 A的特征值都在n个圆盘的并集中.,证明 只就给出证明. 设为A的特征值，即,Ax=x，其中x=(x1,x2, xn)T0.,或,记 ，考虑Ax=x的第k个方程，即,于是,即,这说明，A的每一个特征值必位于A的一个圆盘中，并且相应的特征值一定位于第k个圆盘中(其中k是对应特征向量x绝对值最大的分量的下标).,利用相似矩阵性质，有时可以获得A的特征值进一步的估计，即适当选取非奇异对角阵,并做相似变换 .适当选取 可使某些圆盘半径
6、及连通性发生变化.,例2 估计矩阵A的特征值范围，其中,解 矩阵A的3个圆盘为,由定理8，可知A的3个特征值位于3个圆盘的并集中，由于D1是孤立圆盘，所以D1内恰好包含A的一个特征值1(为实特征值)，即,A的其它两个特征值2, 3包含在D2, D3的并集中.,现在取对角阵,做相似变换,矩阵A1的3个圆盘为,显然，3个圆盘都是孤立圆盘，所以，每一个圆盘都包含A的一个特征值(为实特征值)，且有估计,当A为实矩阵，如果限制用正交相似变换，由于A有复的特征值， A不能用正交相似变换约化为上三角阵. 用正交相似变换能约化到什么程度呢？,定理9 (Schur定理) 设ARnn，则存在酉矩阵U使,其中rii
7、(i=1,2,n)为 A的特征值.,下面给出理论上有关通过酉相似变换及正交变换可以约化一般矩阵A到什么程度的问题.,其中Rii(i=1,2,m)为一阶或二阶方阵，且每个一阶Rii是A的实特征值，每个二阶对角块Rii的两个特征值是 A的两个共轭复特征值.,定理10 (实Schur分解) 设ARnn，则存在正交矩阵Q使,定义4 设ARnn为对称矩阵，对于任一非零向量x，称,我们转向实Schur型的实际计算.,为对应于向量x的瑞利(Rayleigh)商.,定理11 设ARnn为对称矩阵(其特征值次序记为12n)，则,1. (对任何非零xRn)；,2. ；,3. .,证明 只证1，关于2, 3自己作练
8、习.,由于A为实对称矩阵，可将 1,2,n 对应的特征向量 x1,x2,xn 正交规范化，则有 (xi, xj)=ij，设x0为Rn中任一向量，则有,于是,从而1成立. 结论1说明瑞利商必位于n和1之间.,关于计算矩阵A的特征值问题，当n2,3时，我们还可按行列式展开的办法求()=0的根. 但当n较大时，如果按展开行列式的办法，首先求出()的系数，再求()的根，工作量就非常大，用这种办法求矩阵的特征值是不切实际的，由此需要研究求A的特征值及特征向量的数值解法.,本章将介绍一些计算机上常用的两类方法，一类是幂法及反幂法（迭代法），另一类是正交相似变换的方法（变换法）.,幂法与反幂法都是求实矩阵的
9、特征值和特征向量的向量迭代法，所不同的是幂法是计算矩阵的主特征值(矩阵按模最大的特征值称为主特征值，其模就是该矩阵的谱半径)和相应特征向量的一种向量迭代法，而反幂法则是计算非奇异(可逆)矩阵按模最小的特征值和相应特征向量的一种向量迭代法. 下面分别介绍幂法与反幂法.,8.2 幂法及反幂法,现讨论求1及x1的方法.,设实矩阵A=(aij)有一个完全的特征向量组，即A有n个线性无关的特征向量，设矩阵A的特征值为1,2,n, 相应的特征向量为x1,x2,xn. 已知A的主特征值1是实根，且满足条件,8.2.1 幂法（又称乘幂法）,幂法的基本思想是: 任取非零的初始向量v0 , 由矩阵A构造一向量序列
10、vk,称为迭代向量，由假设，v0可唯一表示为,于是,其中,由假设 故 从而,为1的特征向量.,所以当k充分大时，有,即为矩阵A的对应特征值1 的一个近似特征向量.,用(vk)i 表示vk的第i个分量，则当k充分大时，有,即为A的主特征值1的近似值.,由于,这种由已知非零向量v0及矩阵A的乘幂Ak构造向量序列vk以计算A的主特征值1(2.7)及相应特征向量(2.5)的方法就称为幂法.,迭代公式实质上是由矩阵A的乘幂 Ak与非零向量v0相乘来构造向量序列vk=Akv0，从而计算主特征值1及其对应的特征向量，这就是幂法的思想.,的收敛速度由比值,来确定，r越小收敛越快，但当r1时收敛可能很慢.,定理
11、12 设ARnn有n个线性无关的特征向量，主特征值1满足条件,|1|2|n|，,则对任何非零向量v0(a10)，幂法的算式成立.,又设A有n个线性无关的特征向量，1对应的r个线性无关的特征向量为x1,x2,xr，则由(2.2)式有,如果A的主特征值为实的重根, 即1=2=r, 且,|r|r+1|n|，,为A的特征向量，这说明当A的主特征值是实的重根时，定理5的结论还是正确的.,应用幂法计算A的主特征值1及其对应的特征向量时，如果|1|1(或|1|1)，迭代向量 vk的各个不等于零的分量将随 k 而趋向于无穷(或趋向于零) ，这样在计算机实现时就可能 “溢出”. 为克服这个缺点，就需要将迭代向量
12、加以规范化.,设有一向量v0，将其规范化得向量为,其中max(v)表示v的绝对值最大的分量. 即如果有,则max(v)=vq，且q为所有绝对值最大的分量中的最小下标.,在定理12的条件下幂法可这样进行：任取一初始向量v00(a10)，构造规范化向量序列为,由(2.3)式,收敛速度由比值r=|2/1|确定. 总结上述结论，有,同理，可得到,定理13 设ARnn有n个线性无关的特征向量，主特征值1满足|1|2|n|，则对任意非零初始向量v0=u0(a10)，有幂法计算公式为,例1 用幂法计算矩阵,的主特征值与其对应的特征向量.,解 取 v0=u0=(0,0,1)T , 则,直到k=8 时的计算结果
13、见下表,从而,见书p303-例3.,8.2.2 幂法的加速方法,1、原点平移法,由前面讨论知道，应用幂法计算A的主特征值的收敛速度主要由比值 r=|2/1|来决定，但当r 接近于1时，收敛可能很慢. 这时，一个补救办法是采用加速收敛的方法.,其中p为参数，设A的特征值为i，则对矩阵B的特征值为i-p ，而且A, B的特征向量相同.,引进矩阵 B=A-pI .,如果要计算A的主特征值1, 只要选择合适的数p，使1-p为矩阵B=A-pI 的主特征值，且,那么，对矩阵B=A-pI应用幂法求其主特征值1-p, 收敛速度将会加快. 这种通过求B=A-pI的主特征值和特征向量，而得到A的主特征值和特征向量
14、的方法叫原点平移法. 对于A的特征值的某种分布，它是十分有效的.,例4 设AR44有特征值,比值r=|2/1|0.9. 做变换,B=A-12I (p=12),则B的特征值为,应用幂法计算B的主特征值1的收敛速度的比值为,虽然常常能够选择有利的p值, 使幂法得到加速, 但设计一个自动选择适当参数p的过程是困难的.,下面考虑当A的特征值是实数时，怎样选择p使采用幂法计算1得到加速.,且使收敛速度的比值,设A的特征值都是实数，且满足,则对实数p，使矩阵A-pI的主特征值为1-p或n-p时，当我们计算1及x1时，首先应选取p使,显然，当2-p=-(n-p )时，即 P=(2+n)/2P* 时为最小值，
15、这时收敛速度的比值为,当希望计算n时，应选取 p=(1+n-1)/2P* 使得应用幂法计算n得到加速.,当A的特征值都是实数，满足,且2, n能初步估计出来，我们就能确定P*的近似值.,例2 用原点平移加速法求例1中矩阵A的主特征值与其对应的特征向量.,对B应用幂法，仍取 v0=(0,0,1)T , 则,解 取p=-2.5, 做平移变换B=A-pI，则,迭代5步的计算结果见下表,可得到B的主特征值为 113.5000, 主特征向量为 v1 (0.5 ,1.0, 0.7500)T , 因此，A的主特征值为 1 = 1 +p 11.0000, 主特征向量仍为 x1 =(0.5,1,0.7500)T .,原点位移的加速方法，是一个矩阵变换方法. 这种变换容易计算，又不破坏矩阵A的稀疏性，但p的
