加速度时域信号如何转化为加速度功率谱密度?

zhuhuanuaa

加速度时域信号如何转化为加速度功率谱密度?是用PSD函数吗？大家看看下面这样对吗？
%绘制加速度
figure(1)
t1=ScopeData21.time;
s1=ScopeData21.signals.values;
tt=0:0.005:10;
plot(t1,s1,'b:');
%%%绘制车身加速度PSD图
figure(11)
xn1=s1;
window=boxcar(length(xn1));%加窗
nfft=1024;  %%%这个数值不知道怎么定的？
fs=200;%采样频率
[pxx,f]=periodogram(xn1,window,nfft,fs);
title('车身加速度功率谱');
ylabel('Acceleration PSD(m^2/s^4)/HZ');
xlabel('Frequency (Hz)');
plot(f,pxx,'b--','linewidth',2)

kyu16866

这难道用的就是传说中的周期图法？那个值一般取2的几次幂。代表取样总数，f=fs/nfft;表示频率分辨率