基于Robust-L1 的相位恢复

压缩感知领域的，基于Robust-L1 的相位恢复算法。

clear all;
clc;
%% 产生随机稀疏图片
x0 =zeros（1616）;
x0（37）=255;
x0（62）=255;
x0（410）=255;
x0（1012）=255;
x0（137）=255;
imshow（x0[]）;
x0 =double（x0）;
%% 选omige
N =256;
M =N/2;
ind =randperm（N）;
omige1 =ind（1:N-M）;
omige1 =omige1‘;
omige2 =ind（N-M+1:N）;
omige2 =omige2‘;
%% 随机取幅值
fx =fft2（x0）; 
b1 =abs（fx）;
b2 =reshape（b1[N 1]）;
b2（omige2） =0;
b =reshape（b2[1616]）;
%% 更新Z
alafa =0.1;
beta =0.75;
y =0;
lamda =8;
Iter =500;

z0 =randn（1616）;
fz =fft2（z0）;
z0 =b*（fz/abs（fz））;
for k =1:Ite

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1068  2015-05-22 16:43  random.m