用matlab求信号的频谱

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