维纳滤波实现散焦模糊恢复

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（‘模糊退化且添加高斯噪声