Can any body tell me how to implement KSVD in matlab?

Question

Yuva Kishore Reddy Nayannagari 2016-3-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/271507-can-any-body-tell-me-how-to-implement-ksvd-in-matlab

编辑： Royi Avital 2020-12-19

i am doing the project on "quaternion matrix analysis of color image using vector sparse model" so i need to implement ksvd so can please anynoe help me???

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Royi Avital 2016-3-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/271507-can-any-body-tell-me-how-to-implement-ksvd-in-matlab#answer_212331

编辑：Royi Avital 2020-12-19

在 MATLAB Online 中打开

Here is a function which implements the K-SVD algorithm.

function [ mD ] = DL_KSVD( mD, mX, paramCardinality )
      numberOfAtoms = size(mD, 2);
      for ii = 1:50
          %<! Update the Representations
          mA = Omp(mX, mD, paramCardinality);
          %<! Update the Dictionary
          for jj = 1:5
              mE = mX - (mD * mA);
              for kk = 1:numberOfAtoms
                  vP         = find(mA(kk, :));
                  mEP        = mE(:, vP) + mD(:, kk) * mA(kk, vP);
                  vA         = mD(:, kk)' * mEP;
                  mA(kk, vP) = vA;
                  mD(:, kk)  = mEP * vA' / (vA * vA');
              end
          end
          mD = bsxfun(@rdivide, mD, sqrt(sum(mD .^ 2, 1)));
          %%Progress (Debug):
  %         mA     = omp(mD' * mX, mD' * mD, cardinality);
  %         mX_hat = mD * mA;
  %         A = mean(abs(mX(:) - mX_hat(:)))
      end
end

The Omp() function mA = Omp(mX, mD, paramCardinality) solves mD * mX - mA_2 s.t. each column of mA has cardinality (L0 Pseudo Norm) less equal to paramCardinality. It uses the Orthogonal Matching Pursuit method.