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