qrtz.m

来自「<精通MATLAB科学计算>一书的例子代码」· M 代码 · 共 8 行

function l = qrtz(A,M)

for(i=1:M)                 %M为迭代步数
    [q,r]=qr(A);
    A = r*q;
    l = diag(A);
end