求助，关于图像傅里叶变换加窗的问题

[tjucruiser]

求图像频谱，图像大小256*256，分别加上hann. barret，welch窗等，不知道这个窗怎么做

用的是matlab编程

我写了代码，但是感觉不对
%生成窗
w=hann(256);
[M,N]=meshgrid(w,w);
Win_han=M.*N;

%图像加窗IMG.*Win_han
f=fft2(IMG.*Win_han,256,256);
F=fftshift(f);
NPS=F.*conj(F)/(256*256);

请高手指点！

[beyond071]

f = fft2(IMG, 256, 256);%图像二维傅立叶变换
F = fftshift(f);%中心变换
Q = F.*Win_han;%滤波
R = ifftshift(Q);%中心逆变换
U = ifft2(R);%二维傅立叶逆变换
imshow(real(U));%显示实部

[beyond071]

通过上面可以验证下LZ加的窗是否正确，求频谱就在二维傅立叶逆变换之前求就OK了