如何在python中将x轴加速器信号转换为频谱图？

Question

我的数据集（患者编号、时间/毫秒、x、y、z、标签）

1,15,70,39,-970,0
1,31,70,39,-970,0
1,46,60,49,-960,0
1,62,60,49,-960,0
1,78,50,39,-960,0
1,93,50,39,-960,0
.
.
.

我试图在预处理阶段使用谱图作为 x 轴信号，然后将其用作机器学习模型的输入数据，而不是使用原始的原始 x 轴数据

这是我尝试做的事情

import matplotlib.pyplot as plt
import numpy as np

dt = 0.0005
t = np.arange(0.0, 20.0, dt)

data = np.loadtxt("trainingdataset.txt", delimiter=",")
x = data[:]

NFFT = 1024       # the length of the windowing segments
Fs = int(1.0/dt)  # the sampling frequency

ax1 = plt.subplot(211)
plt.plot(x)
plt.subplot(212, sharex=ax1)
Pxx, freqs, bins, im = plt.specgram(x, NFFT=NFFT, Fs=Fs, noverlap=900)
plt.show()

它给我以下错误

Warning (from warnings module):
  File "C:\Users\hadeer.elziaat\AppData\Local\Programs\Python\Python36\lib\site-packages\matplotlib\axes\_axes.py", line 7221
    Z = 10. * np.log10(spec)
RuntimeWarning: divide by zero encountered in log10

Answer 1

如果x 是您的信号并且您可以假设您的采样率是time/millisecond 的平均值，那么您可能可以使用librosa 库来计算使用librosa.feature.melspectrogram 的梅尔谱图，还有其他用于计算信号相关特征的实用程序。