mp 算法 matlab代码,MP 算法 matlab实现

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

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

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

function a=mp(y,phi,K,iware)

[~,N]=size(phi);

a = zeros(N,1);

for i = 1:N %原子归一化

phi(:,i) = phi(:,i)/norm(phi(:,i));

end

v =y; %初始残差

numIts = 0;

iware=2*iware;

while numIts

numIts = numIts + 1;

g=phi'*v; %计算内积

[~,pos]=max(abs(g)); %选取最大的作为最匹配原子

a(pos)=a(pos)+g(pos); %记录位置

v = v -g(pos)*phi(:,pos); %更新残差

end

a=a/iware;

end