用matlab求信号的频谱

Question

现在有一信号，E=cos((1.9*10e+14)*t)*(1900*cos(3*10e+10)*t+pi/2)
要求这个信号的频谱，并作图，各路高手给点意见，要给出matlab的代码，先谢过大家了。
虽然对高手来说很简单，可是我很急啊，大家帮帮忙吧！
E=cos((1.9*10e+14)*t)*(1900*cos((3*10e+10)*t+pi/2))
上面漏了个括号，看这个

Answer 1

fs=1500; %自己设置采样频率
N=4000; %自己设置采样点数

t = (0:N-1)/fs; %间隔
NFFT = 2^nextpow2(N);%转化为2的基数倍
f= fs/2*linspace(0,1,NFFT/2); %求出FFT转化频率

E=cos((1.9e14).*t).*(1900*cos(3e10).*t+pi/2); %函数

E_change=fft(E,NFFT)/N; %进行FFT变换

plot(f,2*abs(E_change(1:NFFT/2)),'b');
%画出频谱特性图
grid on;

上面是比较正规的fft变换，有的文章中没有求FFT转化频率（这样就少了一些步骤），即没有转化为2的基数倍，虽然可以画图也不会出错，但是如果详细阅读matlab自带的help帮助的话，就可以发现它们还是有差别的。本程序自己可以设置采样频率和采样点数。

由于matlab进行FFT变换后画出的图是个对称图，这点你可参阅其fft的帮助，根据帮助文献进行了(1:NFFT/2)处理，画出一半即可。

如果还想画出相频特性图，可以代码后面加上：

ph_E=180*angle(E*N)/pi;
plot(f,ph_E(1:NFFT/2));
grid on;
%画出相频特性图

还有，你的函数由于不能直接复制到matlab里面运行，在matlab里
1.9e14表示 1.9*10^14