qr方法计算中小型矩阵的全部特征值

《qr方法计算中小型矩阵的全部特征值》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《计算方法》课程设计报告学生姓名:学号：学院:班级:题目:QR方法计算中小型矩阵的全部特征值指导教师：职称：教授讲师实验师2015年12月31日II目录目录……………………………………………………………………………I一、选题背景…………………………………………………………………11.1QR方法………………………………………………………………11.2矩阵的特征值………………………………………………………1二、算法设计…………………………………………………………………12.1QR方法的理论………………………………………………………12.2

2、基本QR方法………………………………………………………22.3HouseholderQR分解……………………………………………22.4带原点位移的QR方法……………………………………………3三、程序设计及功能说明……………………………………………………43.1主要程序以及主要功能……………………………………………43.1.1矩阵约化为上海森伯格矩阵………………………………43.1.1QR方法求实矩阵全部特征值………………………………4四、结果分析…………..……………………………………………………6五、总结及心得体会…..………………

3、……………………………………19参考文献…………….…..……………………………………………………20源程序………………..…..……………………………………………………21II一、选题背景1.1QR方法矩阵是高等数学中的常用工具，在很多方面都有重要运用，而矩阵特征值问题在许多领域的研究中有重要的地位，矩阵特征值的一些基本计算方法，研究不同种类矩阵的计算方法和最优计算方法.其中求解矩阵的普通方法包括传统的求法以及初等变换求矩阵的特征值方法；其他的一些优化方法包括幂法、反幂法、Jacobi方法、QR方法.在实际的求解矩阵特征值的问题，根

4、据矩阵的不同特点，选择最快速的方法求解，从而达到最优化解决实际问题。QR分解法是目前求一般矩阵全部特征值的最有效并广泛应用的方法，一般矩阵先经过正交相似变化成为Hessenberg矩阵，然后再应用QR方法求特征值和特征向量。它是将矩阵分解成一个正规正交矩阵Q与上三角形矩阵R，所以称为QR分解法，与此正规正交矩阵的通用符号Q有关。目前QR方法主要用来计算：(1)Hessenberg矩阵的全部特征值问题；(2)计算对称三对角矩阵的全部特征值问题。1.2矩阵的特征值关于计算矩阵A的特征值问题，当n=2,3时，我们还可按行列式展开的办法求，我

5、们还可按行列式展开的办法求det(λ)=0的根，但当n较大时，如果按展开行列式的办法，首先求出det(λ)的系数，再求det(λ)的根，工作量就非常大，用这种办法求矩阵用这种办法求矩阵的特征值是不切实际的的特征值是不切实际的，由此需要研究求A的特征值及特征向量的数值解法。二、算法设计2.1QR方法的理论QR方法的理论：对任意一个非奇异矩阵（可逆矩阵）A，可以把它分解成一个正交阵Q和一个上三角阵R的乘积，称为对矩阵A的QR分解，即A=QR。如果规定R的对角元取正实数，这种分解是唯一的。若A是奇异的，则A有零特征值。任取一个不等于A的特征

6、值的实数μ，则A-μI是非奇异的。只要求出A-μI的特征值和特征向量就容易求出矩阵A的特征值和特征向量，所以假设A是非奇异的，不失一般性。212.2基本QR方法设A=A1，对A1作QR分解，得A1=Q1R1，，交换该乘积的次序，得，由于Q1正交矩阵，A1到A2的变换为正交相似变换，于是A1和A2就有相同的特征值。一般的令A1=A，对k=1，2，3，…..A(k)=QkRk（QR分解）A(k-1)=RkQk(迭代定义)这样，可得到一个迭代序列{Ak}，这就是QR方法的基本过程。2.3HouseholderQR分解从QR算法的构造过程可以

7、看到算法的主要计算量出现在QR分解上，如果直接对矩阵A用QR方法求全部特征值，那麽涉及的计算量是很大的，因此应该先对A作预处理。应用中常先对做正交相似变换将其化为上Hessenberg矩阵H，然后再对H采用QR方法，可以大大减少计算量，这里Hessenberg矩阵也称为拟三角矩阵，它的非零元素比三角矩阵多了一条次对角线，其形式为：上Hessenberg矩阵下Hessenberg矩阵一般矩阵相似约化到Hessenberg矩阵的方法。定理:任取非零向量X=(X1,X2,.....,Xn)T∈Rn，可以选择一个Householder矩阵P，

8、使Px=－σｅ1式中e1=(1,0,......,0)T是Rn的单位向量，证：作u=x+σｅ1，用u做一个Householder矩阵P=I－β-1uuT,因为而证毕21定理中β=σ(σ+x1)，为避免σ+x1出现两个相似