如何运用matlab实现快速傅立叶变化FFT？

咋画出x2(t)=sin(2pi15t) + sin(2pi40t)这个信号的幅度谱和相位谱啊。我自己运行出来的跟手画的咋不一样