?? eig_qr.m
字號:
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);
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -