MP 算法 matlab实现

同学委托写个MP算法，找了一堆资料都没找到好用的，全是OMP什么的，自己记下来以后自用

        MP算法的介绍我上篇博客介绍的很清楚

        MP算法的关键就在于要归一化，matlab第一次用，搞了好久才搞定

[plain] view plain copy print ?

1. function a=mp(y,phi,K,iware)  
2. [~,N]=size(phi);  
3. a = zeros(N,1);  
4. for i = 1:N  %原子归一化  
5.     phi(:,i) = phi(:,i)/norm(phi(:,i));  
6. end  
7. v =y; %初始残差  
8. numIts = 0;  
9. iware=2*iware;  
10. while numIts<iware; %迭代计数  
11.       numIts = numIts + 1;                                 
12.       g=phi'*v; %计算内积  
13.       [~,pos]=max(abs(g)); %选取最大的作为最匹配原子  
14.       a(pos)=a(pos)+g(pos);  %记录位置  
15.       v = v -g(pos)*phi(:,pos); %更新残差  
16.         
17. end   
18. a=a/iware;   
19. end