第三章 第五题答案

答：A=imread('F:\论文图片\Lena256256','bmp');

B=rgb2gray(A);

B=double(B);

s=size(B);

M=s(1);N=s(2);

for u=0:M-1

    for v=0:N-1

        k=0;

        for x=0:M-1

            for y=0:N-1

                k=B(x+1,y+1)*exp(-j*2*pi*(u*x/M+v*y/N))+k;

            end

        end

        F(u+1,v+1)=k;

    end

end

subplot(1,3,1);imshow(A)

subplot(1,3,2);imshow(abs(F))