l1 魔术错误 - matlab

l1 magic error - matlab

我正在尝试使用比信号维度(压缩感知)更少的观测值来重建图像。我正在尝试 运行 以下代码 -

A = imread('cameraman.png');

x_i = 37;
y_i = 95;
s = 35;

A = A([x_i:x_i+s],[y_i:y_i+s]);
x = double(A(:));
figure(1),imshow(A)
xlabel('original')
n=length(x);
m=floor(n/3);
Phi=randn(m,n);   %Measurment Matrix
Psi=dftmtx(n);   %sensing Matrix( or can be dct(eye(n)) )
y=Phi*x;  %compressed signal
Theta=Phi*Psi;
%Initial Guess:  y=Theta*s => s=Theta\y
s2=Theta\y;
%Solution
s1=l1qc_logbarrier(s2, Theta,[], y, 1e-1, 1e-1);
%Reconstruction
x1=Psi*s1;
figure,imshow(reshape(x1,size(A)), [0 256] ),xlabel('OMP')

但是在 运行 编译代码时出现以下错误。

Error using linsolve
Matrix must be positive definite.

Error in l1qc_newton (line 92)
    [dx,hcond] = linsolve(H11p, w1p, opts);

Error in l1qc_logbarrier (line 104)
  [xp, up, ntiter] = l1qc_newton(x, u, A, At, b, epsilon, tau, newtontol,
  newtonmaxiter, cgtol, cgmaxiter);

Error in cs_image2 (line 23)
s1=l1qc_logbarrier(s2, Theta,[], y, 1e-1, 1e-1);

上面这段代码恢复了感知矩阵为dct矩阵而不是dft矩阵时的解。有人可以指出错误所在吗?这是l1-magic固有的问题吗?使用不同的求解器对我有用吗?

注意: l1qc_logbarrier 是库 l1 magic 中的一个函数。 http://users.ece.gatech.edu/justin/l1magic/index.html

l1qc_logbarrier 解决 -

错误似乎是因为使用 dftmtx 输出的矩阵,其元素是 complex 数字,而 l1-magic 不适用于这种情况。我试图解决的问题是 l1 最小化 -

这里 is complex. So l1-magic fails. To get over it, I made the following adjustments. Let , where R is the real and C the imaginary part of and similarly, let 。这给出了,

我们知道b是实数,因此,

合并 2 个约束,我们得到,

这给了我们一个 real 形式的约束 可以使用 l1-magic 来解决。