deconvblind

I=imread（‘E:\test.jpg‘）;
I=im2double（I）;%转换原图格式uint至双精度浮点型
subplot（221）;
imshow（I）;
title（‘original‘）;
PSF=fspecial（‘gaussian‘710）;%构造高斯噪声点扩散函数，包含两个参数：模版尺寸与表示滤波器的标准差
V=.0001;
BlurredNoisy=imnoise（imfilter（IPSF）‘gaussian‘0V）;%添加均值为0，方差为V的高斯噪声至滤波掩模后的图上
subplot（222）;
imshow（BlurredNoisy）;
title（‘blurred‘）;
WEIGHT=zeros（size（I））;
WEIGHT（5:end-45:end-4）=1;
INITPSF