SPAMS 信号稀疏表示Matlab工具包windows版本

function tmp = displayPatches（D）

V=1;
[n K] = size（D）;
size（D）;
sizeEdge=sqrt（n/V）;
if floor（sizeEdge） ~= sizeEdge
   V=3;
   sizeEdge=sqrt（n/V）;
end

p=3.5;
        M=max（（D（:）））;
        m=min（（D（:）））;
        if （m >= 0）
           me=0;
           sig=sqrt（mean（（（D（:）））.^2））;
        else
           me=mean（D（:））;
           sig=sqrt（mean（（（D（:）-me））.^2））;
        end
        D=D-me;
        D=min（max（D-p*sig）p*sig）;
        M=max（（D（:）））;
        m=min（（D（:）））;
       D=（D-m）/（M-m）;
  %     D=1-D;

nBins=floor（sqrt（K））;

tmp = zeros（（sizeEdge+1）*nBins+1（sizeEdge+1）*nBins+1V）;
patch = zeros（sizeEdgesizeEdge1）;
mm=sizeEdge*sizeEdge;
for ii = 1:nBins
   for jj = 1:nBins
      patchCol=（D（1:n（ii-1）*nBins+jj））;
      patchCol=reshape（patchCol [sizeEdgesizeEdge V]）;

      M=max（（patchCol（:）））;
      m=min（（patchCol（:）））;
  %    patchCol=1.0*（patchCol-m）/（M-m）+0.0;


      tmp（（ii-1）*（sizeEdge+1）+2:ii*（sizeEdge+1）...
         （jj-1）*（sizeEdge+1）+2:jj*（sizeEdge+1）:）=patchCol;
   end
end

colormap（‘bone‘）;
imagesc（tmp）;