求高手指正用matlab求解AX=b的预处理公轭梯度法

tangle1984

求解的问题是Ax=b，其中A为28*792，b为28*1，解x为792*1的矩阵。
防真结果出现的问题是图象质量结果不理想，而且时间要得比没预处理的时间更长！
高手指教啊，请发到tanger1984@163.com或者就在论坛回帖指教。十万火急！

1. load F:\biyesheji\data2\zxL\da_ss01%Initialization
   
2. load F:\biyesheji\data2\zxL\da_cnij1.txt
   
3. J=ss01';
   
4. A=J'*J;
   
5. I=eye(792,792);
   
6. [U,S,V]=svd(A+0.5*I);
   
7. [m,n]=size(S);
   
8. M=zeros(n,m);
   
9. for i=1:n
   
10. if S(i,i)>0.01
    
11. M(i,i)=1/S(i,i);
    
12. end
    
13. end %%%%%%%%%%%%%%得到预处理矩阵M
    
14. max_it=40;
    
15. b=J'*da_cnij1;
    
16. x=zeros(length(J),1);
    
17. r=b;
    
18. z=M*b;
    
19. p=z;
    
20. error=norm(r);
    
21. resvec(1)=error;
    
22. tic
    
23. % Begin iteration
    
24. for iter=1:max_it
    
25. ap=A*p;
    
26. normm=r'*z;
    
27. alaf=normm/(ap'*p) ;
    
28. x=x+alaf*p; %Update solution
    
29. r=r-alaf*ap;%compute residual
    
30. z=M*r;
    
31. error=norm(z); %check convergence
    
32. resvec=[resvec;error];
    
33. beta=(r'*z)/normm;
    
34. p=z+beta*p;
    
35. end
    
36. toc
    
37. save F:\biyesheji\data2\CL\tangle.txt x -ascii

复制代码

tangle1984

沉了,郁闷,有没有人有这个好的程序
给我发下啊