I = imread('cameraman.tif');           %读入matlab自带图像cameraman.tif，I

 subplot(2,2,1);imshow(I);              %显示图像I,位置在2行2列的第一个位置

 title('Original Image');               %设置图像I标题为"Original Image"

 axis on;                               %显示图像I像素坐标

 H = fspecial('motion',20,45);          %按角度45°移动20像素的运动滤波器

 MotionBlur = imfilter(I,H,'replicate');%图像大小通过复制外边界的值来扩展

 subplot(2,2,2); imshow(MotionBlur);    %I为输入图像，H为滤波掩模

 title('Motion Blurred Image');         %因为滑动窗口在边缘的时候不好处理，

 axis on;                             %所以输出的图像就直接用原来的边界值赋值

 H = fspecial('disk',10);               %有（10*2+1）个边的圆形均值滤波器

 blurred = imfilter(I,H,'replicate');

 subplot(2,2,3);imshow(blurred); 

 title('Blurred Image');

 axis on;

 H = fspecial('unsharp');               %根据alpha决定拉氏算子创建的模板

 sharpened = imfilter(I,H,'replicate');

 subplot(2,2,4);imshow(sharpened);

 title('Sharpened Image');

 axis on;