P=imread('shuzi.jpg');
imshow(P)
P1=im2bw(P);
P2 = ~P1;
s=size(P2)
I= imread('2.png');
I1=im2bw(I);
I2= ~I1;
C=fftshift(ifft2(fft2(P2).*conj(fft2(I2))));
figure;
imshow(C);
thresh=170;
figure;
imshow(C>thresh);
F=conj(C).*C;
figure;
mesh(F);
本程序摘自《MATLAB仿真及其在光学课程中的应用》一书第六章