二维离散傅里叶变换DFT基函数可视化

N x N图像的第( k,l  )个傅里叶基函数：x方向的周期： k/N；y方向的周期：l/N

% 用图像绘制：
figure(1);
M = 16;%图片的宽度
N = 16;%图片的高度
A = zeros(M,N);

for k = 0:M/2
    for l = 0:N/2
            for x = 1:M
                for y = 1:N
                  subplot(M/2+1,N/2+1,(k*(N/2+1)+l+1)); %画子图
                  A(x,y) = cos(2*pi*k*x/M + 2*pi*l*y/N) + sin(2*pi*k*x/M + 2*pi*l*y/N);%求第（k,l）个傅里叶基函数在（x,y)处的值。为可视化虚部也看做实数
                end
            end
            fprintf('%d',k)
            fprintf('%d\n',l)
            colormap('Cool');%选择色调，这里为Cool
            imagesc(A); % 由A(x,y)的值确定颜色
            %set(gca,'xaxislocation','top');  % 把x轴换到上方
            set(gca,'ydir','reverse');       % 把y轴刻度倒置
            title( {['B(', num2str(k), ',',num2str(l),')']})
    end
end