imfilter

PSF=fspecial('motion',40,75)%%%運動模糊算子，有兩個參數，表示攝像物體逆時針方向以

%%%%%%theta角度運動了len個像素，len的默認值為9，theta的默認值為

MF=imfilter(I,PSF,'circular'); %%%圖像空域濾波

noise=imnoise(zeros(size(I)),'gaussian',0,0.001);%%生成各種噪聲圖像

MFN=imadd(MF,im2uint8(noise));%%運動模糊圖像'

NSR=sum(noise(:).^2)/sum(MFN(:).^2);

imshow(deconvwnr(MFN,PSF,NSR),[])%逆濾波復原'

[reg1 LAGRA]=deconvreg(MFN,PSF,NP/3.0);%使用約束最小二乘方濾波復原圖像