function []=blurt(image)
c=image(:,:,1);
T=size(c);
fx=-T(1)/2+1:1:T(1)/2;
fy=-T(2)/2+1:1:T(2)/2;
% fx=1:1:T(1);
% fy=1:1:T(2);
a=1/100;b=1/100;
MTF=sinc(a.*fx)'*sinc(b.*fy);
f_image=(fftshift(fft2(c,T(1),T(2))));%.*MTF;
%f_image=fft2(c,T(1),T(2));
c0=f_image.*MTF;
c1=ifft2(f_image);
figure,imshow(c0,[]);
figure,imshow(f_image,[])
figure,imshow(MTF,[]);
figure,imshow(c1,[]);
figure,imshow(c,[])
% image = double(imread('lena.bmp'));
% blurt(image);
