clear;clc;clf;
I=imread（‘war01.bmp‘）;
subplot（221）;imshow（I）;title（‘原圖像‘）;
[mn]=size（I）;
F=fftshift（fft2（I））;
k=0.0025;
for?u=1:m?
????for??v=1:n
????????H（uv）=exp（（-k）*（（（u-m/2）^2+（v-n/2）^2）^（5/6）））;
????end
end
G=F.*H;
I0=real（ifft2（fftshift（G）））;
I1=imnoise（uint8（I0）‘gaussian‘00.001）;
subplot（222）;?
imshow（uint8（I1））;title（‘模糊退化且添加高斯噪聲