Fs=512;                                 %采樣頻率
t=0:1/Fs:1.99805;
load dtdld03.txt;                       %添加信號文本文件
x=dtdld03(1:1024,4);                    %設置輸入信號x
figure
subplot(2,1,1);
plot(t,x);                              %繪制濾波前的信號波形
grid;                 
title('濾波前')
%periodogram函數(shù)進行功率譜分析
nfft=512;                                           %設定FFT算法的長度
n1=0:1/Fs:1;
Nn=length(n1);
index=0:round(nfft/2-1);
k=index*Nn/nfft;
window=boxcar(length(x));                          %指定采用的窗函數(shù)及窗長
[Pxx,f]=periodogram(x,window,nfft,Fs);             %periodogram函數(shù)直接法的功率譜估計
plot_Pxx=10*log10(Pxx(index+1));
subplot(2,1,2);
plot(k,plot_Pxx);                                   %繪制功率譜的分析波形
title('periodogram函數(shù)計算的功率譜')