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