clc;
clear;
close?all;
%%?原圖像
RGB=imread（‘Daxing?Airport.jpg‘）;
image_o=rgb2gray（RGB）;
subplot（221）;imshow（image_o）;
title（‘原圖像‘）;

%%?圖像退化+噪聲
f=im2double（image_o）;%圖像轉換為雙精度，像素值縮放到?[0?1]
F=fft2（f）;%快速傅里葉變換
F=fftshift（F）;%中心化，將低頻分量移至頻譜中心
[MN]=size（F）;
[uv]=meshgrid（1:M1:N）;%生成二維坐標系
h=exp（-0.0025*?（?（u-M/2）.^2+（v-N/2）.^2）.^（5/6）?）;%大氣湍流模型
H=h‘;
F=F.*H;
X=ifftshift（F）;%逆中心化
x=ifft2（X）;%逆快速傅里葉變換
x=uint8（abs（x）*256）;
image_t=imnoise（x‘gaussian‘00.005）;%添加高斯噪聲
subplot（222）;imshow（image_t）;
title（‘帶噪退化圖像‘）;

%%?維納濾波
I=im2double（image_t）;
f_I=fft2（I）;
f_I=fftshift（f_I）;
fH_I=f_I;

D=abs（H）;
D=D.^2;
[M1N1]=size（fH_I）;
K=0.05;
threshold=60;
????%對一定半徑范圍內進行濾波
????for?i=1:M1
????????for?j=1:N1
?????????????if?sqrt（（i-M1/2）.^2+（j-N1/2）.^2）?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????129923??2019-12-19?14:58??Daxing?Airport.jpg

?????文件???????1676??2020-06-25?17:43??homework.m

?????文件????????499??2019-12-19?22:09??psnr_mse.m

-----------?---------??----------?-----??----

???????????????132098????????????????????3