?clear?all;
I=imread（‘peppers.png‘）;
I=I（125+[1:256]1:256:）;????%圖像剪切
subplot（231）;imshow（I）;
xlabel（‘（a）原始圖像‘）;
%模擬運動模糊和噪聲
LEN=11;??????????????????????%設(shè)置長度
THETA=5;????????????????????%設(shè)置角度
PSF=fspecial（‘motion‘LENTHETA）;??%生成濾波器
blurred=imfilter（IPSF‘conv‘）;???%圖像卷積運算
subplot（232）;imshow（blurred）;
xlabel（‘（b）模糊圖像‘）;
%添加高斯噪聲
V=0.02;
BlurredNoise=imnoise（blurred‘gaussian‘0V）;???%添加高斯白噪聲
subplot（233）;imshow（BlurredNoise）;
xlabel（‘（c）帶高斯噪聲模糊圖像‘）
%使帶高斯噪聲模糊圖像復(fù)原
NP=V*prod（size（I））;????????%噪聲功率
[J1LAGRA]=deconvreg（BlurredNoisePSFNP）;?%圖像復(fù)原
subplot（234）;imshow（J1）;
xlabel（‘（d）使用真實NP復(fù)原‘）
J2=deconvreg（BlurredNoisePSFNP*1.3）;???%采用放大1.3倍的NP
subplot（235）;imshow（J2）;
xlabel（‘（e）采用放大1.3倍的NP復(fù)原‘）
%使用小于真實NP值的復(fù)原
J3=deconvreg（BlurredNoisePSFN