tmp610.m

来自「信号分析与处理的一些源代码」· M 代码 · 共 21 行

%MATLAB PROGRAM 6-10
clf
%Estimate PSD by Multitaper Method
%=================================
N=1024;
Nfft=256;
Fs=1000;
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);
subplot(211)
plot(f,10*log10(Pxx1))
xlabel('Frequency (Hz)');
ylabel('Power Spectrum (dB)');
title('Multitaper Method(MTM) nw=4')
grid

subplot(212)
pmtm(xn,4,Nfft,Fs,[],0.5);