qrsolv.m

用matlab写的一些数值算法

function  x = qrsolv(H,R, b)
% Solve an overdetermined system Ax = b by means of Housholder 
% transformation.  H,R  is the output from qrfac(A)

% Version 11.12.2003.  INCBOX

[m,n] = size(H);
for  k = 1 : n
  ii = k:m;  b(ii) = apphouse(H(ii,k),b(ii));
end
x = R \ b(1:n);