线性调频信号的FFT谱

laughing

　　线性调频信号x(t)=exp(-0.3t)cos(pi*t^2+2*pi*t)
　　其Fourier幅值谱图为：


　　小弟的弱弱问题是：如何用matlab画出上述x(t)的Fourier幅值谱图?
　　小弟的程序是：
　　t=0:10/1000:10;
　　x=exp(-0.3*t).*cos(pi*t.^2+2*pi*t);
　　Xf=fft(x);
　　plot(Xf,x)
　　作出的图形是：



　　明显不对。请问错在哪儿?怎么改?

[ 本帖最后由 zhangnan3509 于 2007-7-25 16:52 编辑 ]

eight

原帖由 laughing 于 2007-7-25 15:31 发表
线性调频信号x(t)=exp(-0.3t)cos(pi*t^2+2*pi*t)
其Fourier幅值谱图为：
28253

小弟的弱弱问题是：如何用matlab画出上述x(t)的Fourier幅值谱图？
小弟的程序是：
t=0:10/1000:10;x=exp(-0.3*t).*cos(pi* ...

至少应该是：plot(real(Xf))

[ 本帖最后由 eight 于 2007-7-25 16:09 编辑 ]

zhangnan3509

咱们版里有很多FFT的例子，找找对照一下吧:handshake

sagan

t=0:10/1000:10;
x=exp(-0.3*t).*cos(pi*t.^2+2*pi*t);
Xf=fft(x);
X=sqrt(Xf.*conj(Xf));
y=X(1:100)
plot(y);

VibrationMaster

plot(Xf,x)%-->以Xf为横坐标(而且还是复数),x为纵坐标画图,这种图画出来的意思是什么

plot(abs(Xf));%幅值谱
plot(real(Xf));%实部
...........

saintlei

没取幅度

kerrylau2000

t=0:10/1000:10;
x=exp(-0.3*t).*cos(pi*t.^2+2*pi*t);
Xf=fft(x);
y=abs(Xf);
plot(y(1:500));%幅值谱

%you can try it.