qrgivens.m

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

function [Q,R,c,s]= qrgivens(H)
[m,n]=size(H);
for k=1:n-1
    [c(k),s(k)]=givcos(H(k,k),H(k+1,k));
    H=garow(H,c(k),s(k),k,k+1,k,n);
end
R=H; Q=prodgiv(c,s,n);
return