I = imread('cameraman.tif'); 
PSF = fspecial('gaussian',60,10);
%PSF點擴展函數
J = edgetaper(I,PSF);
%對圖像邊緣模糊化處理
subplot(1,2,1);imshow(I,[]);title('原始圖像');
subplot(1,2,2);imshow(J,[]);title('邊緣模糊化');