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