zhuhanshu.m

在对矩阵进行拟上三角化的基础上利用带双步位移的QR分解法求解矩阵的特征值

clc
clear
format long
for i=1:10
    for j=1:10
        if i~=j
            A(i,j)=sin(0.5*i+0.2*j);
        else
            A(i,j)=1.5*cos(i+1.2*j);
        end
    end
end
[iter,D]=qralg(A);
disp('特征值是：')
D'