在一个512*512像素尺寸的黑色背景上叠加一个20*40像素尺寸的白色矩形
此图像在进行二维傅里叶变换之前乘上了(-1)^(x+y),从而使得频率谱关于中心对称
img=zeros(512,512);
img(246:266,236:276)=1;
sft=zeros(512,512);
for i=1:512
for j=1:512
sft(i,j)=(-1)^(i+j);
end
end
Y=fft2(img.*sft);
imgo=log(abs(Y));
figure
imshow(imgo,[]);