H_QR.m

function [Q,R]=H_QR(A)
% household 的QR方法,A=QR，matlab源程序已经具备，大家可以按照这个参考一下计算原理
[m,n]=size(A);Q=eye(n);
%Ak=zeros(n0,n0);
   for j=1:n-1 
      sigma=sign(A(j,j))*norm(A(j:n,j),2);
      u=A(j:n,j)+sigma*eye(n-j+1,1);
      beta=sigma*(sigma+A(j,j));
      I=eye(n-j+1);
      Hj=I-1/beta*(u*u');
      Hk=blkdiag(eye(j-1),Hj);
      Q=Q*Hk;
      A=Hk*A;
   end 
      R=A;