求大神帮助

[对月临风]

I=imread('SATURN.bmp');
imgpinpu(I);


function imgpinpu(I)
imshow(I);
fftI=fft2(I);
sfftI=fftshift(fftI);
rr=real(sfftI);
ii=imag(sfftI);
A=sqrt(rr.^2+ii.^2);
A=(A-min(min(A)))/(max(max(A))-min(min(A)))*255;
figure;
imshow(A);






function imgpinpu(I) 这个函数怎么编写 谢谢了 有关傅里叶变换的