eig_qr.m

来自「使用QR分解方法计算矩阵特征值的matlab源码」· M 代码 · 共 17 行

function r=eig_qr(A)
%qr方法计算矩阵特征值
A=hessenberg(A);
%将矩阵转化为上Hessenberg矩阵
[q1,r1]=qr_decomposition(A);
%对上Hessenberg矩阵进行qr分解，直到得到的新矩阵的对角线元素趋近于其特征值
A1=r1*q1;
[q2,r2]=qr_decomposition(A1);
A2=r2*q2;
while norm(diag(A1-A2))>1e-5
%设定允许误差
[q1,r1]=qr_decomposition(A2);
A1=r1*q1;
[q2,r2]=qr_decomposition(A1);
A2=r2*q2;
end
r=diag(A2);