I=imread（‘E:\test.jpg‘）;
I=im2double（I）;%轉(zhuǎn)換原圖格式uint至雙精度浮點(diǎn)型
subplot（221）;
imshow（I）;
title（‘original‘）;
PSF=fspecial（‘gaussian‘710）;%構(gòu)造高斯噪聲點(diǎn)擴(kuò)散函數(shù)，包含兩個參數(shù)：模版尺寸與表示濾波器的標(biāo)準(zhǔn)差
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