matlab实现经典功率谱估计

clf fs=100 %采样点数 N=128 %数据点数 n=0:N-1;t=n/fs %数据对应的时间序列 f1=15;f2=40 x=sin(2*pi*f1*t)+2*sin(2*pi*f2*t) %信号 y=fft(x,N) mag=abs(y) f=n*fs/N %频率序列 subplot(221),plot(f,mag) %绘出随频率变化的振幅 xlabel('频率/Hz');ylabel('振幅');title('N=128');grid on % 振幅的大小与采样点数有关 subplot(222),plot(f(1:N/2),mag(1:N/2)) %绘出Nyquist频率之前随频率变化的振幅 xlabel('频率/Hz');ylabel('振幅') title('N=128');grid on %对信号采样数据为1024点的处理 fs=100 %采样点数 N=1024 %数据点数 n=0:N-1;t=n/fs %数据对应的时间序列 f1=15;f2=40 x=sin(2*pi*f1*t)+2*sin(2*pi*f2*t) %信号 y=fft(x,N) mag=abs(y) f=n*fs/N %频率序列 subplot(223),plot(f,mag) %绘出随频率变化的振幅 xlabel('频率/Hz');ylabel('振幅');title('N=1024');grid on subplot(224),plot(f(1:N/2),mag(1:N/2)) %绘出Nyquist频率之前随频率变化的振幅 xlabel('频率/Hz');ylabel('振幅') title('N=1024');grid on