hessqrtz.m

来自「关于矩阵特征值问题的的源代码」· M 代码 · 共 8 行

function l = hessqrtz(A,M)

A = hess(A);
for(i=1:M)
    [q,r]=qr(A);
    A = r*q;
    l = diag(A);
end