如何检查加速度列中的峰值？答案

【问题标题】：How to check peak in the acceleration column?如何检查加速度列中的峰值？
【发布时间】：2019-10-10 06:57:03
【问题描述】：

我有一个数据框，其中包含“AccelerationG”列作为时间戳。我想查看设备的振动/加速度并检查它是否超过峰值？

AccelerationG
0.95
0.93
1.12
1.12
0.95
0.93
1.12
0.95
1.12
1.12
0.93
0.93
1.12
1.12
0.95
5.42
10.66
14.39



How can I approach this?

【问题讨论】：

df['AccelerationG'].gt(threshold).sum() 你可以根据自己的需要改变阈值，它给出了振动交叉的发生次数
df[df['AccelerationG']>threshold] 为您提供超过阈值的行。
我们可以把它画成图表吗？我的意思是我们可以做快速傅立叶变换，检查峰值和其他东西吗？
你试图绘制什么？

标签： python pandas fft accelerometer vibration

【解决方案1】：

dat=[0.95,0.93,1.12,1.12,0.95,0.93,1.12,0.95,1.12,1.12,
0.93,0.93,1.12,1.12,0.95,5.42,10.66,14.39]

import numpy as np
import matplotlib.pyplot as p
%matplotlib inline

p.figure(figsize=(10,10))
p.subplot(221)
p.plot(dat, '.-')
p.title('all data')

p.subplot(222)
p.plot(dat[:-3],'.-')
p.title('data truncated')

p.subplot(223)
mn=np.mean(dat[:-3])       # DC 
p.plot(dat[:-3]- mn,'.-')  # subtract DC
p.title('DC removed')

p.subplot(224)
p.psd(dat[:-3]-mn,12,1/0.01);  
p.title('power spectral density')

如果数据是在 100 Hz 时采集的，超过 140 毫秒，那么在 35 Hz 处会有一个频率峰值。

【讨论】：