I?=?imread（‘peppers.png‘）;%讀取圖像
I?=?I（10+[1:256]222+[1:256]:）;%選定圖像的范圍
figure（1）;??imshow（I）;%顯示原圖像
title（‘Original?Image‘）;
LEN?=?31;?THETA?=?11;
PSF?=?fspecial（‘motion‘LENTHETA）;%產(chǎn)生運動模糊函數(shù)
blurred?=?imfilter（IPSF‘circular‘‘conv‘）;%產(chǎn)生運動模糊圖像
figure（2）;?subplot（221）?
imshow（blurred）;%顯示運動模糊圖像
title（‘Blurred‘）;
wnr1?=?deconvwnr（blurredPSF）;%使用維納濾波器進行濾波
subplot（222）;?imshow（wnr1）;%顯示維納濾波器后濾波的圖像
title（‘Restored?True?PSF‘）;
wnr2?=?deconvwnr（blurred...
???fspecial（‘motion‘2*LENTHETA））;%使用兩倍長度的點擴散函數(shù)進行維納濾波
subplot（223）;?imshow（wnr2）;?%顯示恢復的圖像
title（‘Restored?“Long“?PSF‘）;
wnr3?=?deconvwnr（blurred...
???fspecial（‘motion‘LEN2*THETA））;%使用兩倍角度的點擴散函數(shù)進行維納濾波
subplot（224）;?imshow（wnr3）;?%顯示恢復的圖像
title（‘Restored?Steep‘）;
noise?=?0.1*randn（size（I））;%生成隨機噪聲
blurredNoisy