基于Wringter Flow的相位恢复算法

基于Wirtinger流的相位恢复问题。简单来说就是在对光场相位进行恢复时，无须对光场函数进行凸集定义，而是用光谱的方法获得一个最优的初始化值，然后应用新的迭代规则调整该初始值从而获得理想结果。

% Implementation of the Wirtinger Flow （WF） algorithm presented in the paper 
% “Phase Retrieval via Wirtinger Flow: Theory and Algorithms“ 
% by E. J. Candes X. Li and M. Soltanolkotabi

% The input data are coded diffraction patterns about a random complex
% valued 1D signal. 

clear all
close all
clc
%% Make signal
n = 128;
x = randn（n1） + 1i*randn（n1）;

%% Make masks and linear sampling operators

L = 6;                   % Number of masks  

% Sample phases: each symbol in alphabet {1 -1 i  -i} has equal prob. 
Masks = randsrc（nL[1i -1i 1 -1]）;

% Sample magnitudes and make masks 
temp = rand（size（Masks））;
Masks = Masks .* （ （temp <= 0.2）*sqrt（3） + （temp > 0.2）/sqrt（2） ）;

% Make linear operators; A is forward map and At its scaled adjoint （At（Y）*nume

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

    .......      2320  2016-07-16 15:41  WF_1D_CDP.m

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

                 2320                    1