Gram-Schmidt算法的matlab程序

野菊

Gram-Schmidt算法方面的一些matlab程序

pasuka

matlab的PDE工具自带的代码，仅供参考

1. function  [y,h]=pdegrmsc(r,v)
   
2. %PDEGRMSC Gram Schmidt phase of Arnoldi algorithm.
   
3. 
   
4. %       A Ruhe 8-08-94.
   
5. %       Copyright 1994-2001 The MathWorks, Inc.
   
6. 
   
7. maxort=4;
   
8. qmax=0.7;
   
9. [n j]=size(v);
   
10. rnn=norm(r);
    
11. hc=zeros(j,1);
    
12. if j>0,
    
13.   for iort=1:maxort,
    
14.     hv=v'*r;
    
15.     hc=hc+hv;
    
16.     r=r-v*hv;
    
17.     rnorm=rnn;
    
18.     rnn=norm(r);
    
19.     if rnn > qmax*rnorm, break, end;
    
20.   end;
    
21. end;
    
22. h=[hc;rnn];
    
23. y=r/rnn;
    

复制代码

Zoyie

谢谢楼主分享~~~~~~

怪咖先生

感谢大家分享