samp9_11.m

MATLAB7.x数字信号处理 光盘内容

%Samp9_11
clf;Fs=1000;   %采样频率
N=1024;Nfft=256;n=0:N-1;t=n/Fs;  %数据长度、分段数据长度，时间序列
randn('state',0);  %设置产生随机数的初始状态
xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N); %原始信号
[Pxx1,f]=pmtm(xn,4,Nfft,Fs); %用多窗口法(NW=4)估计功率谱
subplot(2,1,1),plot(f,10*log10(Pxx1));  %绘制功率谱
xlabel('频率/Hz');ylabel('功率谱/dB');
title('多窗口法(MTM) nw=4');grid on
[Pxx,f]=pmtm(xn,2,Nfft,Fs); %用多窗口法(NW=2)估计功率谱
subplot(2,1,2),plot(f,10*log10(Pxx));  %绘制功率谱
xlabel('频率/Hz');ylabel('功率谱/dB');
title('多窗口法(MTM) nw=2');grid on