作业帮 > 数学 > 作业

matlab中,已知矩阵A,已完成对A的QR分解,下一步求A的特征值和特征向量,程序怎么编写?

来源:学生作业帮 编辑:拍题作业网作业帮 分类:数学作业 时间:2024/04/29 18:32:16
matlab中,已知矩阵A,已完成对A的QR分解,下一步求A的特征值和特征向量,程序怎么编写?
已针对矩阵A用自己编写函数[q,r]=qrhs(A)完成了QR分解,要求A的特征值和特征向量该怎么求呢?只会求其特征值,特征向量不知道从哪写起.我的算法:
function I = tzzhs(A,M,j)%基于householder法求特征值的函数
%六阶矩阵:A 求矩阵特征值:I
for i=1:M %多次迭代
[q,r]=qrhs(A);
A = r*q;
I = diag(A);
end
如何在上述程序的基础上求出特征向量?求程序
2 对求的特征向量验证正确性
3 我要分析迭代次数M对特征值的影响,取M=5,10,50,100,但是对六阶实对称阵而言,现象不明显;对六阶复对称真而言,没有什么规律可言.是什么原因?急,不尽感激!
楼主的问题是自己写程序完成矩阵的QR分解,既然是迭代实现QR分解,就与矩阵论中说的计算特征值和特征向量的方法有些区别了.大体的步骤应该是首先将矩阵化成双对角矩阵,然后追赶计算特征值和特征向量,程序代码可以参考 徐士良编的 常用数值算法 c语言描述
再问: 是的,前面已完成了QR分解,就剩求特征值和特征向量,并完成迭代次数对特征值影响的分析了,是要用matlab的语言编写的,能在我编写的程序上略加指点吗?因为比较急,可能没时间看参考书了。 另外问一下:对A采用了两种QR基本算法得到的两个正交阵Q1和Q2,假设他俩是差不多的,就是对应元素相同,现在要对其分析画图等,即使验证这两个矩阵是接近相同的,通常有哪几种表述方法,比较直观呢?