GS相位恢复代码

根据Phase retrieval algorithms: a comparison J. R. Fienup写的GS算法，使用软件Matlab

clear;
clc;
im=imread（‘peppers.png‘）;
img=double（im）;
f=fft2（img）;
f=fftshift（f）;
m=abs（f）;
xn=ones（256256）;
for n=1:200
    G=fft2（xn）;
    G=fftshift（G）;
   
phase_G=angle（G）;
G=m.*cos（phase_G）+m.*sin（phase_G）.*（-1i）;
g=ifft2（G）;

phase_g=angle（g）;
xn=img.*cos（phase_g）+img.*sin（phase_g）.*（-1i）;
end
x=uint8（real（xn））;
imshow（x）
p=snr（imx）;

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

     文件      37559  2006-04-17 11:20  GS\Cameraman.png

     文件        382  2016-05-10 15:41  GS\main.asv

     文件        381  2016-05-10 16:59  GS\main.m

     文件      40181  2005-11-14 08:14  GS\peppers.png

     文件        182  2016-05-10 15:36  GS\snr.asv

     文件        109  2016-05-10 16:06  GS\snr.m

     文件         64  2016-05-10 16:09  GS\Untitled3.m

     目录          0  2016-05-17 15:23  GS

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

                78858                    8