img?=?double（imread（‘standard_lena.bmp‘））;
n=512;
uo=zeros（nn）;
uo=img（1:n1:n）;
init?=?2055615866;
randn（‘seed‘?init）;
sigma?=?1;
figure;imagesc（uo）;colormap（gray）;
un?=sigma?*?randn（size（uo））;
u?=?uo?+25*un;
SNR1（uou）
figure;imagesc（u）;colormap（gray）;
u1=fft2（u）;
a=1.5;
t=4^（-a）;
j=sqrt（-1）;
g1=zeros（nn）;
g2=zeros（nn）;
g3=zeros（nn）;
g4=zeros（nn）;
for?i=1:n
????p（i）=（1-exp（-2*3.1415926535*j*（i-1）/n））^a*exp（2*3.1415926535*j*（i-1）/n）^（a/2）;
????q（i）=conj（p（i））;
end
for?k=1:3300
????k
????ty=0;
????for?i=1:n
????????g1（:i）=p（i）*u1（:i）;
????????g2（i:）=p（i）*u1（i:）;
????end
????g1=real（ifft2（g1））;
????g2=real（ifft2（g2））;
????g3=g1./abs（1+g1.^2+g2.^2）;
????g4=g2./abs（1+g1.^2+g2.^2）;
????g3=fft2（g3）;
????g4=fft2（g4）;
????for?i=1:n
????????g1（:i）=q（i）*g3（:i）;
????????g2（i:）=q（i）*g4（i:）;
????end
????u1=u1-（g1+g2）*t;
????ut=ifft2（u1）;
????SNR1（uoreal（ut））
end
ut=ifft2（u1）;
figure;imagesc（real（ut））;colormap（gray）;
SNR1（uout）
SNR1（uereal（ut））

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-11-22?22:56??frac\
?????文件??????263222??2009-06-11?13:13??frac\BARB.BMP
?????文件?????????249??2009-06-11?13:13??frac\SNR1.m
?????文件????????7168??2011-07-18?12:45??frac\Thumbs.db
?????文件????????1012??2009-11-27?15:29??frac\dqzb.m
?????文件????????1114??2009-06-11?13:13??frac\fdqz.m
?????文件?????????164??2009-06-11?13:13??frac\fractiond2.m
?????文件?????????172??2009-06-11?13:13??frac\fractiondc2.m
?????文件??????263222??2009-06-11?13:13??frac\standard_lena.bmp