function [x1] = OMP(K,A,y)
%UNTITLED 此处显示有关此函数的摘要
% 此处显示详细说明
%Initialize
x=0;
r=y;%残差向量
I=[];%索引集合
B=[];%储存更新之后的A列向量
k=0;%循环索引
while(k<K)
in_product=A'*r;
[val pos]=max(abs(in_product));%返回内积绝对值中,最大元素的值及其位置
I=[I;pos];%更新索引
B=[B A(:,pos)];%采样矩阵中的列
%size(B)
x=pinv(B)*y;%最小二乘法求目标信号
r=y-B*x;%更新残差向量
k=k+1;
end
[M N]=size(A);
x1=zeros(N,1);
x1(I)=x;
end
OMP正交匹配追踪算法
于 2020-10-30 18:16:11 首次发布