%?2D最大似然法（空域），卷積函數：g=imfilter?（f?h?‘conv‘?‘replicate‘?）。

clear
%?f=imread（‘cell.tif‘）;?????%?讀入圖像f
f=imread（‘Lena256.tif‘）;?????%?讀入圖像f
f=im2double（f）;???????????%?從8?位型轉為雙精度以便后面的運算
figure（1）subplot（131）imshow（f）;??
R=3;
h=fspecial（‘disk‘R）;?????%?創建離散點擴散函數h
%?h=fspecial（‘gaussian‘51）;????%?創建離散點擴散函數h
g=imfilter?（f?h?‘conv‘?‘replicate‘?）;?%?原圖像f與點擴散函數h卷積，得散焦模糊圖像g。加邊界處理，所以無黑邊
figure（1）subplot（132）imshow（g[]）;???