MOPKSVD等算法的图像稀疏表示代码

图像稀疏表示代码，实现输入图像的稀疏表示，有MOP,KSVD等算法。

 function [Ax ] = KSVD（ycodebook_sizeerrGoal）
%==============================
%input parameter
%  y - input signal
%  codebook_size - count of atoms
%output parameter
%  A - dictionary
%  x - coefficent
%==============================
%SIZ = size（y2）;
if（size（y2） < codebook_size）
   disp（‘codebook_size is too large or training samples is too small‘）;
   return;
end
% initialization
[rowscols]=size（y）;
r=randperm（cols）;   %随机排列
A=y（:r（1:codebook_size））;
A=A./repmat（sqrt（sum（A.^21））rows1）;
ksvd_iter=10;
% main loop
for k=1:ksvd_iter
    % sparse coding
    if nargin==2
        x=OMP（Ay5.0/6*rows）;
      
    elseif nargin==3
        x=OMPerr（AyerrGoal）;
    end
    for m=1:codebook_size
        mindex=find（x（m:））;
        if ~isempty（mindex）
            mx=x（:mindex）;
            mx（m:）=0;
            my=A*mx;
            resy=y（:mindex）;
            mE=resy-my;
            [usv]=svds（mE1）;
            A（:m）=u;
            x（mmindex）=s*v‘;
        end
    end
end


end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1057  2014-08-25 16:16  S_R\KSVD.m

     文件       1074  2014-08-25 11:13  S_R\MOD.m

     文件        536  2014-08-21 14:05  S_R\M_RGB_SPA.m

     文件       1662  2014-08-25 16:20  S_R\M_SPA.m

     文件        710  2014-06-22 15:44  S_R\OMP.m

     文件       1123  2014-06-22 15:40  S_R\OMPerr.m

     目录          0  2014-08-25 17:58  S_R

----------- ---------  ---------- -----  ----

                 6162                    7