python - 如何从 fft 函数中获取频率轴?
问题描述
所以,我可能遗漏了一些明显的东西,但是我搜索了很多教程和文档,似乎找不到一个直接的答案。如何找到在 Python 中执行 fft 的函数的频率轴(特别是 scipy 库中的 fft)?
我正在尝试获取原始 EMG 信号,对其执行带通滤波器,然后执行 fft 以查看剩余的频率分量。但是,我不确定如何找到准确的 x 组件列表。我目前正在处理的特定信号以 1000 Hz 采样,有 5378 个样本。
它只是创建一个从 0 开始到 fft 数据长度的线性 x 吗?我看到很多人创建一个从 0 到采样点乘以采样间距的 linspace。但是在这种情况下我的样本间距是多少?它只是样本/采样率吗?还是完全是别的东西?
解决方案
这是一个例子。
首先创建一个预先确定采样间隔的正弦波。我们将合并两个频率为 20 和 40 的正弦波。请记住,如果时间间隔很大,高频可能会混叠。
#Import the necessary packages
from scipy import fftpack
import matplotlib.pyplot as plt
import numpy as np
# sampling freq in herts 20Hz, and 40Hz
freq_sampling1 = 10
freq_sampling2 = 20
amplitude1 = 2 # amplitude of first sine wave
amplitude2 = 4 # amplitude of second sine wave
time = np.linspace(0, 6, 500, endpoint=True) # time range with total samples of 500 from 0 to 6 with time interval equals 6/500
y = amplitude1*np.sin(2*np.pi*freq_sampling1*time) + amplitude2*np.sin(2*np.pi*freq_sampling2*time)
plt.figure(figsize=(10, 4))
plt.plot(time,y, 'k', lw=0.8)
plt.xlim(0,6)
plt.show()
注意图中两个正弦波是叠加的。一个带频率。10和幅度2和另一个与频率。20和幅度4。
# apply fft function
yf = fftpack.fft(y, time.size)
amp = np.abs(yf) # get amplitude spectrum
freq = np.linspace(0.0, 1.0/(2.0*(6/500)), time.size//2) # get freq axis
# plot the amp spectrum
plt.figure(figsize=(10,6))
plt.plot(freq, (2/amp.size)*amp[0:amp.size//2])
plt.show()
请注意,在幅度谱中,两个频率被恢复,而在其他频率处幅度为零。幅度值也分别为 2 和 4。
您可以改为使用tom10fftpack.fftfreq
建议的频率轴
因此,代码更改为
yf = fftpack.fft(y, time.size)
amp = np.abs(yf) # get amplitude spectrum
freq = fftpack.fftfreq(time.size, 6/500)
plt.figure(figsize=(10,6))
plt.plot(freq[0:freq.size//2], (2/amp.size)*amp[0:amp.size//2])
plt.show()
我们只绘制幅度谱的正部分[0:amp.size//2]
推荐阅读
- javascript - document.title 元素在刷新时消失
- spring - 如何在 Angular 8 中调用 Rest API Post 方法。出现 CORS 错误
- javascript - 如何解决 CSS 没有反映在 Air-Table 的 iframe 上的问题?
- amazon-web-services - 使用 AWS Amplify 将 React Native 应用程序发布到云
- curl - 使用 curl 将 post 请求发送到 mlflow api 到多条记录
- nuxt.js - Nuxt.Js 在 Azure 上部署没有 node_modules 文件夹
- java - 如何在不将活动切换到android中的默认呼叫应用程序的情况下拨打电话?
- azure - Azure - 无法获得“网络配置文件”
- android - CameraX 预览未显示在 PreviewView 中
- distance - 测量两个概率分布之间的重叠