最大似然法实现图像盲复原，可以随意换图像

% 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[]）;