从空间域到频域的代码的问题

读入一张图片，输出移中后的频域的图片。以下是我根据离散的傅里叶变换编的代码，一运行MATLAB就死了，请指出错误或给个正确的代码。
function DFT
p=imread('图片');
q=p;
pi=3.14;
e=2.72;
j=sqrtm(-1);
[m,n]=size(p);
for i=1:m
    for j=1:n
        for k=1:m
            for l=1:n
                q(i,j)=(-1)^(k+l)*p(k,l)*e^(-j*2*pi*(double(i)*k/m+double(j)*l/n))/(m*n);    %(-1)^(k+l)移中
            end
        end
    end
end
imshow(q);

[ 本帖最后由 碳水化合物 于 2011-5-19 18:26 编辑 ]