clf;
clc,clear,close all;I=imread('peppers.png'); % 读入图像
I1=rgb2gray(I); % 图像灰度转换
figure;imshow(I);title('Original Image');
figure;imshow(I1);title('Original Grayimage');
PSF=fspecial('gaussian',7,10); % 设置高斯低通滤波器
Blurred=imfilter(I1,PSF,'symmetric','conv'); % 图像卷积计算,symmetric图像大小通过镜像反射其边界来扩展
figure;imshow(Blurred); title('Blurred Image');
UNDERPSF = ones(size(PSF)-4); % 欠尺寸PSF,比原PSF小4
[J1,P1] = deconvblind(Blurred,UNDERPSF); % J1表示图形,P1表示PSF模型
figure;imshow(J1);title('Deblurring with Undersized PSF'); % 欠尺寸PSF,图像复原,比原PSF小4
OVERPSF = padarray(UNDERPSF,[4 4],'replicate','both'); % 过尺寸PSF,比原PSF大4
[J2,P2] = deconvblind(Blurred,OVERPSF); % padarray填充图像或填充数组
figure;imshow(J2);title('Deblu