请问关于 fftshift 的结果解释

yuyizhen2004

大虾可以运行下面代码：
a=1;
c=5;

[x,y]=meshgrid(1:128,1:128);
gauss=zeros(128,128);
b_x=64;
b_y=64;
gauss=a*exp(-((x-b_x).^2+(y-b_y).^2)/c^2);
P=fftshift(fft2(gauss));
Q=fftshift(fft2(fftshift(gauss)));
figure
subplot(121);imshow(P);title('fftshift(fft2(gauss))');
subplot(122);imshow(Q);title('fftshift(fft2(fftshift(gauss)))');
figure
subplot(121);mesh(P);title('fftshift(fft2(gauss))');
subplot(122);mesh(Q);title('fftshift(fft2(fftshift(gauss)))');

请问大虾这个结果如何解释？

[ 本帖最后由 eight 于 2008-5-12 18:59 编辑 ]

sigma665

请说清楚你的问题
什么样的结果
有无错误提示

xjzuo

试试如下代码：
%%%==============================%%%
clear all
a=1;
c=5;
[x,y]=meshgrid(1:128,1:128);
gauss=zeros(128,128);
b_x=64;
b_y=64;
gauss=a*exp(-((x-b_x).^2+(y-b_y).^2)/c^2);
P=fftshift(fft2(gauss));
Q=fftshift(fft2(fftshift(gauss)));
figure
subplot(121);imshow(P);title('fftshift(fft2(gauss))');
subplot(122);imshow(Q);title('fftshift(fft2(fftshift(gauss)))');
figure
subplot(121);mesh(real(P));title('fftshift(fft2(gauss))');
subplot(122);mesh(real(Q));title('fftshift(fft2(fftshift(gauss)))');
%%%============================================%%%

注：一般加abs较多．

[ 本帖最后由 xjzuo 于 2008-5-11 16:11 编辑 ]