%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