I = checkerboard(8);
% 創建一個棋盤圖像
PSF = fspecial('gaussian',7,10);
%創建點擴散函數PSF
V = .01;
BlurredNoisy = imnoise(imfilter(I,PSF),'gaussian',0,V);
%對模糊圖像添加高斯噪聲
NOISEPOWER = V*prod(size(I));
%噪聲功率
[J LAGRA] = deconvreg(BlurredNoisy,PSF,NOISEPOWER);
subplot(221); imshow(BlurredNoisy);
title('A = Blurred and Noisy');
subplot(222); imshow(J);
title('[J LAGRA] = deconvreg(A,PSF,NP)');
subplot(223); imshow(deconvreg(BlurredNoisy,PSF,[],LAGRA/10));
title('deconvreg(A,PSF,[],0.1*LAGRA)');
subplot(224); imshow(deconvreg(BlurredNoisy,PSF,[],LAGRA*10));
title('deconvreg(A,PSF,[],10*LAGRA)');