householder.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 16 行

function A=householder(A)
[m,n]=size(A);
for k=1:n
x=A(k:m,k);
y=x'*x;
y=sqrt(y);
v(k:m,k)=x;
if x(1)>=0
   v(k,k)=y+x(1);
else v(k,k)=x(1)-y;
end
r=v(k:m,k)'*v(k:m,k);
r=sqrt(r);
v(k:m,k)=v(k:m,k)/r;
A(k:m,k:n)=A(k:m,k:n)-2*v(k:m,k)*(v(k:m,k)'*A(k:m,k:n));
end