GS GS算法

GS算法，图像经过多次傅里叶变换，最后提取图像的相位信息，加载空间光调制器，观察图像。

clear
clc
A=imread（‘Eason.jpg‘）;
B= imresize（A[512 512]）;
C = im2double（rgb2gray（B））;
C1=abs（C）;
for p=1:100;
    C=fftshift（C）; 
    D=ifft2（C）;       
    D=fftshift（D）;
    D2=angle（D）;
    D=exp（i.*D2）;    %给D振幅赋值1，相位不变
    D=fftshift（D）;
    C=fft2（D）;       %D傅里叶变换回到C
    C=fftshift（C）;
    C2=angle（C）;
    C=C1.*exp（i.*C2）;   %C振幅不变，仍然赋值初始振幅C1
    p=p+1;
end
    C=fftshift（C）; 
    D=ifft2（C）;       
    D=fftshift（D）;
    D=angle（D）;
    imwrite（D‘杨民强.bmp‘‘bmp‘）;
figureimshow（D[]）    

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

     文件      44807  2012-09-12 21:13  GS\Eason.jpg

     文件      30739  2012-09-12 21:14  GS\Eason512.jpg

     文件      29100  2012-09-12 21:17  GS\Eason512gray.jpg

     文件        561  2012-09-18 16:16  GS\GS.m

     文件     263222  2012-09-18 16:16  GS\杨民强.bmp

     目录          0  2012-09-12 23:52  GS

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

               368429                    6