matlab 简单二维函数中心谱

在一个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,[]);