rqrtz.m

特征值计算的几种算法

function l = rqrtz(A,M)

A = hess(A);
for(i=1:M)
    N = size(A);
    n = N(1,1);
    u = A(n,n);
    [q,r]=qr(A-u*eye(n,n));
    A = r*q+u*eye(n,n);
    l = diag(A);
end