t1=1:0.001:2*pi;
w=sin(2*pi*30*t1)+2*sin(2*pi*50*t1);
W=fftshift(fft(w,256));%%傅里叶变换，并将频谱中心搬移到原点，变换点数是256点
fs=linspace(-1000/2,1000/2,256);%%产生频率值，1000是采样间隔，来自于0.001的倒数
plot(fs,abs(W))%%频谱图，总共四个峰值，-50,50，-30,30处。
