Householder变换的matlab实现

happy

1. function [v,b]=Householder(x)
   
2. if size(x,1)==1
   
3.     x=x';
   
4. end
   
5. n=length(x);
   
6. ninf=max(abs(x));
   
7. if ninf==0
   
8.     v=x;
   
9.     b=0;
   
10. return;
    
11. end
    
12. x=x/ninf;
    
13. n2=dot(x(2:n),x(2:n));
    
14. v=x;
    
15. v(1)=1;
    
16. if n2==0
    
17.     b=0;
    
18. else
    
19.     a=sqrt(x(1)*x(1)+n2);
    
20. if x(1)<=0
    
21.     v(1)=x(1)-a;
    
22. else
    
23.     v(1)=-n2/(x(1)+a);
    
24. end
    
25. b=2*v(1)*v(1)/(n2+v(1)*v(1));
    
26. v=v/v(1);
    
27. end
    
28. return

复制代码

happy

1. function [A,d]=HouseholderQR(A)
   
2. [m,n]=size(A);
   
3. I=diag(ones(m,1));
   
4. d=zeros(n,1);
   
5. for j=1:n
   
6.     [v,b]=Householder(A(j:m,j));
   
7.     A(j:m,j:n)=(I(j:m,j:m)-b*v*v')*A(j:m,j:n);
   
8.     d(j)=b;
   
9.     if j<m
   
10.         A(j+1:m,j)=v(2:m-j+1);
    
11.     end;
    
12. end;

复制代码

花如月

这么好的帖子，顶一下。顺便问以下2个函数的功能是一样的吧:loveliness:

chengtianmmz

感谢资源共享，又学到了东西:@)