Question

跪求omp算法的matlab代码，用于稀疏编码的！

%A-稀疏系数矩阵
%D-字典/测量矩阵（已知）
%X-测量值矩阵（已知）
%K-稀疏度
%L-the maximal number of coefficient for representation
function A=OMP(D,X,L)
[n,P]=size(X);
[n,K]=size(D);
for k=1:1:P
a=[];
x=X(:,k);
residual=x;%残差
indx=zeros(L,1);%索引集
for j=1:1:L
proj=D'*residual;%D转置与residual相乘，得到与residual与D每一列的内积值
pos=find(abs(proj)==max(abs(proj)));%找到内积最大值的位置
pos=pos(1);%若最大值不止一个，取第一个
indx(j)=pos;%将这个位置存入索引集的第j个值
a=pinv(D(:,indx(1:j)))*x;%indx(1:j)表示第一列前j个元素
residual=x-D(:,indx(1:j))*a;
end
temp=zeros(K,1);
temp(indx)=a;
A(:,k)=temp;%只显示非零值及其位置
end
return

这里有一个代码，但是老是运行有错误，
Error using OMP (line 6)
Not enough input arguments.
请问是怎么回事啊？

Answer 1

你要从工作空间调用这个函数，OMP之后要输入参数的。。。

追问

你好，也就是说这个函数不能直接调用的是吗？函数本身是没有问题的对吗？

追答

对，你调用要加上参数 OMP(1,2,3)这种形式