untitled2.m

来自「电机电流信号的文本载入」· M 代码 · 共 22 行

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函数进行功率谱分析
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));                          %指定采用的窗函数及窗长
[Pxx,f]=periodogram(x,window,nfft,Fs);             %periodogram函数直接法的功率谱估计
plot_Pxx=10*log10(Pxx(index+1));
subplot(2,1,2);
plot(k,plot_Pxx);                                   %绘制功率谱的分析波形
title('periodogram函数计算的功率谱')