samp9_12.m

最大熵法与Welch法估计的功率谱的比较。 多信号分类法估计的功率谱。

%Samp9_12
clf;Fs=1000;   %采样频率
N=1024;Nfft=256;n=0:N-1;t=n/Fs;   %数据长度、分段长度和时间序列
window=hanning(256);   %采用窗口
randn('state',0);      %设置产生随机数的初始状态
xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);   %原始信号
[Pxx1,f]=pmem(xn,14,Nfft,Fs);   %采用最大熵法，采用滤波器阶数14，估计功率谱
subplot(2,1,1),plot(f,10*log10(Pxx1));   %绘制功率谱
xlabel('频率/Hz');ylabel('功率谱/dB');
title('最大熵法(MTM) Order=14');grid on
%采用Welch方法估计功率谱
noverlap=128;    %重叠数据
dflag='none';
subplot(2,1,2)
psd(xn,Nfft,Fs,window,noverlap,dflag);  %采用Welch方法估计功率谱
xlabel('频率/Hz');ylabel('功率谱/dB')
title('Welch方法估计的功率谱');grid on