n=0:0.1:200;%设定信号时间长度为0到200秒,采样间隔0.1,则采样频率为10HZ,点数2001
y=sin(2*pi*0.2*n)+sin(2*0.213*n);
Y=fft(y);%FFT
Pyy=Y.*conj(Y)/2000;%信号功率谱
f=10*(0:1000)/2000;%计算横轴频率值
figure(1)
subplot(2,1,1),plot(n,y),title('信号'),xlabel('时间(S)')
subplot(2,1,2),plot(f,Pyy(1:1001)),title('信号功率谱'),xlabel('频率(Hz)')