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,但是对六阶实对称阵而言,现象不明显;对六阶复对称真而言,没有什么规律可言.是什么原因?急,不尽感激!
已针对矩阵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,假设他俩是差不多的,就是对应元素相同,现在要对其分析画图等,即使验证这两个矩阵是接近相同的,通常有哪几种表述方法,比较直观呢?
再问: 是的,前面已完成了QR分解,就剩求特征值和特征向量,并完成迭代次数对特征值影响的分析了,是要用matlab的语言编写的,能在我编写的程序上略加指点吗?因为比较急,可能没时间看参考书了。 另外问一下:对A采用了两种QR基本算法得到的两个正交阵Q1和Q2,假设他俩是差不多的,就是对应元素相同,现在要对其分析画图等,即使验证这两个矩阵是接近相同的,通常有哪几种表述方法,比较直观呢?
matlab中,已知矩阵A,已完成对A的QR分解,下一步求A的特征值和特征向量,程序怎么编写?
怎样用matlab程序求矩阵的最大特征值和特征向量
用MATLAB编写矩阵A的LU分解程序.
用matlab求解矩阵a的最大特征值及其特征向量,并将其特征向量归一化,对矩阵a进行
在MATLAB中求矩阵特征值和特征向量的代码
matlab中如何求矩阵的特征值和特征向量
用matlab求解矩阵a的最大特征值及其特征向量,并将其特征向量归一化,对矩阵a进行一致性检验.速求
求矩阵A=(1100)的特征值和特征向量,并判断是否可对角化
已知三阶矩阵A的特征值,以及对应的三个特征向量,求矩阵A.
已知矩阵A的一个特征值为λ,求矩阵E+A的一个特征向量
如何用C语言编写求对称矩阵的特征值和特征向量的程序
已知3阶矩阵A的3个特征值和对应的特征向量,如何求矩阵A?