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);%使用约束最小二乘方滤波复原图像