qr_decomposition.m

来自「使用QR分解方法计算矩阵特征值的matlab源码」· M 代码 · 共 24 行

function [q,r]=qr_decomposition(A)
%利用Givens变换实现对上Heissenberg矩阵的qr分解
[m,m]=size(A);
for k=1:m-1
%反复使用Givens变换
[c,s]=givens(A(k:k+1,k));
G1=[c,s;-s,c];
G2=eye(m);
for  i=k:k+1
   for  j=k:k+1
       G2(i,j)=G1(i-k+1,j-k+1);
   end
end
%Givens变换矩阵
if k==1
G=G2;
else
G=G2*G;
end
A(k:k+1,k:m)=G1*A(k:k+1,k:m);
end
q=G';
r=A;
end