文献综述：关于方程axd=λx的求解

《文献综述：关于方程axd=λx的求解》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、沪J、学本科生毕业论文（设计）文献综述姓名与学号3100102197张蹇煜指导教师程晓良年级与专业2010信息与计算科学所在学院数学系1、求解特征值和特征向量的方法1.1非对称特征值问题的计算方法(1)岸法。岸法是就是那种一个矩阵的模最大特征值和对应的特征向量的一种迭代方法。假定A是可对角化的矩阵。那么可以得到一个式子uk=Ak*uO/Xk,其中uk是A的特征向量空间中的某个向量，那么当k充分大时,uk就是特征值Ak的对应的一个特征向量。可以由此得到一个迭代格式:yk=Auk-1.pk=<(k)j,<(k)j是yk的模最大分量,uk=yk/pk,其中uO是任意给定的初始向量，通常

2、要求

3、

4、u0

5、

6、-二1。如果们想要接着求接下来的特征向量和特征值，那么就需要利用已经求得的特征值把原矩阵降阶。最简单实用的收缩技巧是利用正交变换。假设Axl=Xlxl并假设酉矩阵P使得Pxl=del这里el=(l,0,・・・,0)T。将上式带入上上式并整理可得PAP*el=Xlel即PAP*的右下为一个n-1的矩阵，它的特征值就是A除了入1以外的特征值。那么对PAP*的右下矩阵继续做幕法即可。而得到P的变换可以使用复的Householder变换来实现。(2)反幕法。反幕法又称作反迭代法，就是应用幕法作用在A-1上，来求A的模最小特征值和对应的特征向量。因此，其基本迭代格式为Ayk

7、二zk-1,

8、Jk=Q,是yk的模最大分量,zk=yk/pk。反幕法主要用来求特征向量，是在用某种方法求得A的某个特征值Xi的近似值pi之后,应用反幕法于A-pil上。也就是说，在实际计算中，常用的是带位移的反幕法。设p是给定的位移。带原点位移p的反岸法的迭代格式如下：(A-pI)vk=zk-1,zk=vk/

9、

10、vk

11、

12、2,从上述的迭代格式上可以看出，反幕法每迭代一次就需要解一个线性方程组，这比幕法的运算量大得多。但是，由于方程组的系数矩阵不随着k变化，所以，可以先把系数矩阵做LU分解,然后就只需要解两个三角矩阵方程组就可以了。(3)QR方法。QR方法是自电子计算机问世以来矩阵计

13、算的重大进展之一,也是目前计算一般矩阵的全部特征值和特征向量的最有效方法之一。QR方法是利用正交相似变换将一个给定的矩阵逐步约化为上三角矩阵或拟上三角矩阵的一种迭代方法，其基本收敛速度是二次的，当然原矩阵实对称时,可以达到三次收敛。QR算法的基本迭代格式如下：Am-1=QmRmzAm=RmQmf其中Qm为酉矩阵,Rm为上三角矩阵。QR方法有很多加速和优化的方法,例如双重步位移的QR迭代方法。考虑如下的迭代格式(带位移的QR方法)：H1二QOTAQO,Hk-

14、JkI=QkRk,Hk+l二RkQk+pkl,这样可以在迭代的过程中会有复数的产生。如果一开始就给出两个p值，就可以使用双重

15、步位移的QR方法。可以提供的代码：functionh=DoubleStepQR(H,n)%双重歩位移QR迭代法%输入数据为Hessenberg矩阵H以及矩阵H的阶数neps=lE-10;m=n-1;s=H(m,m)+H(n,n)；t=H(m,m)*H(n,n)-H(m,n)*H(n,m);x=H(1,1)*H(1,1)+H(1,2)*H⑵1)・s*H(1,1)+t;y二H(2,1)*(H(1,1)+H(2,2)・s);z二H(2,1)*H(3,2);fork=0:n-3[v,beta]=house([x,y,z]');if(k>l)elseendi=eye(length(v))-

16、beta*(v*v,);H(k+1:k+3,q:n)=i*H(k+1:k+3,q:n);if(k+4

17、l))