数字信号处理(一)
时频域信号转换与时域滤波
现在有这样一个项目需求:在传感器端只能获取到时域信号相关振动数据(速度、加速度等),要针对振动信号进行故障分析,而要实现它,第一步就是要实现时频域信号间的互相转换,并叠加时域滤波。
FFT快速傅里叶变换
了解FFT相关基础可以参考以下知乎大佬的文章:
要实现时频域信号间的高效转换,就要用到FFT快速傅里叶变换,通过python的第三方scipy库中的fft函数就能快速实现需求。
scipy中有两个模块scipy.fft和scipy.fftpack都可以实现FFT,这两个模块非常相似,但scipy.fft模块较新,应优先于scipy.fftpack使用。
import numpy as np
import matplotlib.pyplot as plt
from scipy import fft, signal
时域图像展示
利用with open() as f方法导入时域信号数据后,可以利用以下代码展示时域图像:
def timeWave(SAMPLE_RATE, DURATION, Slist)