%checkerboard产生测试板图像,第一个参数为每个正方形一边的像素数,第二个为参数行数,第三个为列数 f=checkerboard(8);%产生一个一边为8个正方形的测试板 subplot(2,2,1),imshow(f),title('原图像') PSF=fspecial('motion',7,45)%用fspecial处理图像的运动模糊,PSF为空间滤波器 gb=imfilter(f,PSF,'circular');%减少边界效应 subplot(2,2,2),imshow(gb);title('运动模糊图像') noise=imnoise(zeros(size(f)),'gaussian',0,0.001);%高斯噪声 subplot(2,2,3),imshow(noise);title('运动模糊图像') g=gb+noise;%添加运动模糊的高斯噪声 subplot(2,2,4),imshow(pixeldup(g,8),[]);title('运动模糊+高斯噪声图像')%大图像运算过程过慢,选用小图像来节省时间 1234567891011121314