一共三个function文件,分别如下
function [H,Q,c,s,B]=givensqr(H)
[n,m]=size(H);
Q=eye(n);
for k=1:m
a=H(k,k);b=H(k+1,k);
[c,s]=givens(a,b);
H=gmatmat(H,c,s,k);
M=eye(n);
M(k,k)=c;
M(k,k+1)=-s;
M(k+1,k)=s;
M(k+1,k+1)=c;
Q=QM’;
end
end
%---------------------------------
function B=gmatmat(A,c,s,k)
G(1,1)=c;
G(1,2)=-s;
G(2,1)=s;
G(2,2)=c;
A([k,k+1],:)=GA([k,k+1]😅;
B=A;
end
%---------------------------------
function [c,s]=givens(a,b)
r=sqrt(a2+b2);
c=a/r;
s=-b/r;
end