信息处理 - 使用 FFT 进行振动分析 - 吾爱随笔录

我有 6 个列向量被加速度计用来分析结构的振动。每列是一个值向量，总时间为秒。我做了一些事情，但我不知道如何解决它。 $(\mathbf A_x,\mathbf A_y, \mathbf A_z, \mathbf G_x, \mathbf G_y, \mathbf G_z)$ $64998\times 1$ $664$

Fs=97.885;
L=64998;
NFFT=2^nextpow2(64998);

AXdft = fft(Ax);
AXdft=AXdft(1:NFFT/2+1,:);

freq = linspace(0,Fs,NFFT/2+1);

plot(freq,2*abs(AXdft));
xlabel('Hz');ylabel('Magnitude');

figure;
plot(freq,20*log10(abs(AXdft)));
xlabel('Hz');ylabel('dB');
title('Plot in dB');

顺便说一句，我将通过添加每一列hold on

得到一个数字，但缺少它的大小。频率没有出现在图中，我想我必须做 PSD，有人可以帮忙吗？