详解利用Python scipy.signal.filtfilt() 实现信号滤波


Posted in Python onJune 05, 2019

本文将以实战的形式基于scipy模块使用Python实现简单滤波处理,包括内容有1.低通滤波,2.高通滤波,3.带通滤波,4.带阻滤波器。具体的含义大家可以查阅大学课程,信号与系统。简单的理解就是低通滤波指的是去除高于某一阈值频率的信号;高通滤波去除低于某一频率的信号;带通滤波指的是类似低通高通的结合保留中间频率信号;带阻滤波也是低通高通的结合只是过滤掉的是中间部分。上面所说的内容会在实战部分加以介绍,可以对比理解一下。

如何实现的呢?我的理解,是通过时域转换为频域,在频域信号中去除相应频域信号,最后在逆转换还原为时域型号。具体的内容还是要查阅大学课程,信号与系统。自己学的很一般就不班门弄斧了。

有什么作用呢?My Opinions,可以消除一些干扰信号,以低通滤波为例,例如我们如果只是统计脉搏信号波形,应该在1Hz左右,却发现波形信号上有很多噪音,这些噪音都是成百上千Hz的,这些对于脉搏信号波形就属于无用的噪音,我们就可以通过低通滤波器将超出某一阈值的信号过滤掉,此时得到的波形就会比较平滑了。

在使用Python进行信号处理过程中,利用 scipy.signal.filtfilt()可以快速帮助实现信号的滤波。

1.函数的介绍

(1).滤波函数

scipy.signal.filtfilt(b, a, x, axis=-1, padtype='odd', padlen=None, method='pad', irlen=None)

输入参数:

b: 滤波器的分子系数向量

a: 滤波器的分母系数向量

x: 要过滤的数据数组。(array型)

axis: 指定要过滤的数据数组x的轴

padtype: 必须是“奇数”、“偶数”、“常数”或“无”。这决定了用于过滤器应用的填充信号的扩展类型。{‘odd', ‘even', ‘constant', None}

padlen:在应用滤波器之前在轴两端延伸X的元素数目。此值必须小于要滤波元素个数- 1。(int型或None)

method:确定处理信号边缘的方法。当method为“pad”时,填充信号;填充类型padtype和padlen决定,irlen被忽略。当method为“gust”时,使用古斯塔夫森方法,而忽略padtype和padlen。{“pad” ,“gust”}

irlen:当method为“gust”时,irlen指定滤波器的脉冲响应的长度。如果irlen是None,则脉冲响应的任何部分都被忽略。对于长信号,指定irlen可以显著改善滤波器的性能。(int型或None)

输出参数:

y:滤波后的数据数组

(2).滤波器构造函数(仅介绍Butterworth滤波器)

scipy.signal.butter(N, Wn, btype='low', analog=False, output='ba')

输入参数:

N:滤波器的阶数

Wn:归一化截止频率。计算公式Wn=2*截止频率/采样频率。(注意:根据采样定理,采样频率要大于两倍的信号本身最大的频率,才能还原信号。截止频率一定小于信号本身最大的频率,所以Wn一定在0和1之间)。当构造带通滤波器或者带阻滤波器时,Wn为长度为2的列表。

btype : 滤波器类型{‘lowpass', ‘highpass', ‘bandpass', ‘bandstop'},

output : 输出类型{‘ba', ‘zpk', ‘sos'},

输出参数:

b,a: IIR滤波器的分子(b)和分母(a)多项式系数向量。output='ba'

z,p,k: IIR滤波器传递函数的零点、极点和系统增益. output= 'zpk'

sos: IIR滤波器的二阶截面表示。output= 'sos'

2.函数的使用

信号滤波中最常用的无非低通滤波、高通滤波和带通滤波。下面简单介绍这三种滤波的使用过程:

(1).高通滤波

#这里假设采样频率为1000hz,信号本身最大的频率为500hz,要滤除10hz以下频率成分,即截至频率为10hz,则wn=2*10/1000=0.02

from scipy import signal

b, a = signal.butter(8, 0.02, 'highpass')
filtedData = signal.filtfilt(b, a, data)#data为要过滤的信号

(2).低通滤波

#这里假设采样频率为1000hz,信号本身最大的频率为500hz,要滤除10hz以上频率成分,即截至频率为10hz,则wn=2*10/1000=0.02

from scipy import signal

b, a = signal.butter(8, 0.02, 'lowpass') 
filtedData = signal.filtfilt(b, a, data)    #data为要过滤的信号

(3).带通滤波

#这里假设采样频率为1000hz,信号本身最大的频率为500hz,要滤除10hz以下和400hz以上频率成分,即截至频率为10hz和400hz,则wn1=2*10/1000=0.02,wn2=2*400/1000=0.8。Wn=[0.02,0.8]

from scipy import signal

b, a = signal.butter(8, [0.02,0.8], 'bandpass')
filtedData = signal.filtfilt(b, a, data)  #data为要过滤的信号

 参考:

1.https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.signal.filtfilt.html

2.https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.signal.butter.html#scipy.signal.butter

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python图书管理系统
Apr 05 Python
对python中raw_input()和input()的用法详解
Apr 22 Python
windows下安装Python的XlsxWriter模块方法
May 03 Python
numpy返回array中元素的index方法
Jun 27 Python
python 图片去噪的方法示例
Jul 09 Python
python 使用pdfminer3k 读取PDF文档的例子
Aug 27 Python
PyQT5 emit 和 connect的用法详解
Dec 13 Python
解决pytorch-yolov3 train 报错的问题
Feb 18 Python
30行Python代码实现高分辨率图像导航的方法
May 22 Python
Python使用tkinter实现摇骰子小游戏功能的代码
Jul 02 Python
python读取excel进行遍历/xlrd模块操作
Jul 12 Python
如何将json数据转换为python数据
Sep 04 Python
分享8个非常流行的 Python 可视化工具包
Jun 05 #Python
Django组件cookie与session的具体使用
Jun 05 #Python
使用python3调用wxpy模块监控linux日志并定时发送消息给群组或好友
Jun 05 #Python
python3下载抖音视频的完整代码
Jun 05 #Python
把JSON数据格式转换为Python的类对象方法详解(两种方法)
Jun 04 #Python
Django集成搜索引擎Elasticserach的方法示例
Jun 04 #Python
python添加菜单图文讲解
Jun 04 #Python
You might like
PHP自定义多进制的方法
2016/11/03 PHP
解决php写入数据库乱码的问题
2019/09/17 PHP
实用的JS正则表达式(手机号码/IP正则/邮编正则/电话等)
2013/01/11 Javascript
js左右弹性滚动对联广告代码分享
2014/02/19 Javascript
js单词形式的运算符
2014/05/06 Javascript
PHP中CURL的几个经典应用实例
2015/01/23 Javascript
JavaScript设置body高度为浏览器高度的方法
2015/02/09 Javascript
基于javascript实现简单计算器功能
2016/01/03 Javascript
JavaScript兼容浏览器FF/IE技巧
2016/08/14 Javascript
Vue+Webpack完美整合富文本编辑器TinyMce的方法
2018/11/30 Javascript
详解vue移动端项目代码拆分记录
2019/03/15 Javascript
vue语法自动转typescript(解放双手)
2019/09/18 Javascript
javascrpt密码强度校验函数详解
2020/03/18 Javascript
[01:19:46]EG vs Secret 2019国际邀请赛淘汰赛 胜者组 BO3 第二场 8.21.mp4
2020/07/19 DOTA
用python实现的可以拷贝或剪切一个文件列表中的所有文件
2009/04/30 Python
python使用在线API查询IP对应的地理位置信息实例
2014/06/01 Python
Python入门篇之字符串
2014/10/17 Python
python读写ini配置文件方法实例分析
2015/06/30 Python
python 列表降维的实例讲解
2018/06/28 Python
python中的tcp示例详解
2018/12/09 Python
Python3 ffmpeg视频转换工具使用方法解析
2020/08/10 Python
Python使用struct处理二进制(pack和unpack用法)
2020/11/12 Python
html5的canvas实现3d雪花飘舞效果
2013/12/27 HTML / CSS
德国电子产品购物网站:TechInTheBasket德国
2018/12/07 全球购物
会计助理的岗位职责
2013/11/29 职场文书
小学毕业家长寄语
2014/01/19 职场文书
大学生旷课检讨书
2014/01/22 职场文书
开办大学饮食联盟创业计划书
2014/01/29 职场文书
《我要的是葫芦》教学反思
2014/02/23 职场文书
2014年大堂经理工作总结
2014/11/21 职场文书
南湾猴岛导游词
2015/02/09 职场文书
学校德育工作总结2015
2015/05/11 职场文书
道歉的话怎么说
2015/05/12 职场文书
python的变量和简单数字类型详解
2021/09/15 Python
JVM的类加载器和双亲委派模式你了解吗
2022/03/13 Java/Android
PostgreSQL出现死锁该如何解决
2022/05/30 PostgreSQL