gpt4 book ai didi

python - FFT - 峰峰值、峰值、RMS

转载 作者:行者123 更新时间:2023-12-01 09:03:47 27 4
gpt4 key购买 nike

我的工作涉及振动,我试图从 FFT 幅度中获取以下信息:

  • 峰值到峰值
  • 峰值
  • RMS

我正在对一个简单的正弦波函数执行 FFT,考虑汉宁窗。请注意,正弦波函数的“全振幅”为 5,运行 FFT 下面的代码会得到 2.5 的振幅结果。因此,在这种情况下,我从 FFT 中获得峰值。峰峰值和 RMS 怎么样?

P.-S. - 我对带宽频率的 RMS 不感兴趣(即 parsevall 定理)。我对每个峰值的 RMS 感兴趣,这通常在振动软件中看到。

import numpy as np
import matplotlib.pyplot as plt


f_s = 100.0 # Hz sampling frequency
f = 1.0 # Hz
time = np.arange(0.0, 10.0, 1/f_s)
x = 5 * np.sin(2*np.pi*f*time)
N = len(time)
T = 1/f_s

# apply hann window and take the FFT
win = np.hanning(len(x))
FFT = np.fft.fft(win * x)
n = len(FFT)
yf = np.linspace(0.0,1.0/(2.0*T),N//2)

plt.figure(1)
plt.plot(yf,2.0/N * np.abs(FFT[0:N//2]))
plt.grid()

plt.figure(2)
plt.plot(time,x)
plt.xlabel('time')
plt.ylabel('Amplitude')
plt.grid()
plt.show()

最佳答案

您在频域中得到的峰值为 2.5,因为这是加窗信号的平均幅度,并且您没有补偿窗口权重。使用以下方法对频域结果进行归一化以考虑窗口后:

plt.plot(yf,2.0/win.sum() * np.abs(FFT[0:N//2]))

您应该得到 5 的振幅,就像在时域中一样。请注意,如果输入信号频率是 f_s/N 的精确倍数(在您的情况下为 0.1Hz),并且输入信号是纯音的基本假设,则此操作有效或由频率充分分离的音调组成是有效的。

峰峰值只是幅度的两倍,因此在您的示例中为 10。

对于 RMS 值,您可能对相应时域正弦音调分量的 RMS 值感兴趣(假设输入信号确实由其频率在频率上充分分离的正弦分量组成)。幅度为 A 的时域正弦曲线的 RMS 为 A/sqrt(2),因此您只需除以 sqrt(2) > 从幅度值中获取相应的等效 RMS 值,因此示例中为 5/sqrt(2) ~ 3.53

关于python - FFT - 峰峰值、峰值、RMS,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/52243955/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com