hessqr.m

来自「various matlab programs to slove various」· M 代码 · 共 11 行

function [T,Q,R]=hessqr(A,niter)
n=max(size(A));
[T,Qhess]=houshess(A);
for j=1:niter
    [Q,R,c,s]= qrgivens(T);  
    T=R;
    for k=1:n-1,
        T=gacol(T,c(k),s(k),1,k+1,k,k+1);
    end
end
return