I=imread('moon.bmp'); 
[m n]=size(I);
[x,y]=meshgrid(-n/2:n/2-1,-m/2:m/2-1);
c1=-4*(pi)*(x.^2+y.^2);
figure, imshow(c1,[])

F=fftshift(fft2(I));
cF1=F.*c1; IcF1=ifft2(cF1);

IF=mat2gray(abs(IcF1));
Isharp=mat2gray(I)+IF;

subplot(1,4,1);imshow(I,[]), title('?????? e????a fe??a????');
subplot(1,4,2);imshow(mat2gray(log(1+abs(F)))), title('??F t?? e????a?');
subplot(1,4,3);imshow(mat2gray(log(1+abs(cF1)))), title('?ap?as?a?? f??t??');
subplot(1,4,4);imshow(Isharp), title(' ????s? e????a? µe Laplacian');
