请教在matlab中计算fft的问题

刺猬精灵

ts=0.000000008;
t=0:ts:0.000001;
x1=cos(2*pi*5000000*t);
x2=cos(2*pi*20000000*t);
x3=x1+x2;

N1=length(x1);
ssf1=(-N1/2:N1/2-1)/(ts*N1);
fx1=fft(x1(1:N1));
subplot(311)
plot(ssf1,abs(fx1))

N2=length(x2);
ssf2=(-N2/2:N2/2-1)/(ts*N2);
fx2=fft(x2(1:N2));
subplot(311)
plot(ssf2,abs(fx2))


N3=length(x3);
ssf3=(-N3/2:N3/2-1)/(ts*N3);
fx3=fft(x3(1:N3));
subplot(311)
plot(ssf3,abs(fx3))

信号的频主应该为5M，和20M，怎么绘制出来不对。请问一下哪出错了？

频普图

南宫守护

plot(ssf3,fftshift(abs(fx3)))

刺猬精灵

回复 2# 南宫守护

非常感谢。对了。不过再作fftshift变换，是什么原理呢？

wangzeyue

这个，完全看不懂啊