I=im2double(imread('miss01.jpg'));
subplot(2,2,1),imshow(I);title('原始图像');
PSF=fspecial('motion',31,11);
blurred=imfilter(I,PSF,'conv','circular');
subplot(2,2,2),imshow(blurred);title('模糊仿真');
noise_mean=0;
noise_var=0.0001;
blurred_noisy=imnoise(blurred,'gaussian',noise_mean,noise_var);
subplot(2,2,3),imshow(blurred_noisy);title('模糊与噪声仿真');
estimated_nsr=noise_var/var(I(:));
wnr=deconvwnr(blurred_noisy,PSF,estimated_nsr);
subplot(2,2,4),imshow(wnr);title('维纳滤波复原图像');
I=im2double(imread('miss01.jpg'));
PSF=fspecial('gaussian',7,10);
V=.01;
BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V);
subplot(1,2,1),imshow(BlurredNoisy),title('模糊与噪声');
NOISEPOWER=V*prod(size(I));
J=deconvreg(BlurredNoisy,PSF,NOISEPOWER);
subplot(1,2,2);imshow(J),title('约束最小二乘方复原图像');
I=im2double(imread('lena01.jpg'));
subplot(2,3,1),imshow(I);title('原始图像');
PSF=fspecial('motion',3,1);
blurred=imfilter(I,PSF,'conv','circular');
subplot(2,3,2),imshow(blurred);title('模糊仿真');
noise_mean=0;
noise_var=0.01;
blurred_noisy=imnoise(blurred,'gaussian',noise_mean,noise_var);
subplot(2,3,3),imshow(blurred_noisy);title('模糊与噪声仿真');
NOISEPOWER=noise_var*prod(size(I));
J=deconvreg(blurred_noisy,PSF,NOISEPOWER);
subplot(2,3,4);imshow(J),title('约束最小二乘方复原图像');
estimated_nsr=noise_var/var(I(:));
wnr=deconvwnr(blurred_noisy,PSF,estimated_nsr);
subplot(2,3,5),imshow(wnr);title('维纳滤波复原图像');