详解利用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中查找excel某一列的重复数据 剔除之后打印
Feb 10 Python
Python开发之快速搭建自动回复微信公众号功能
Apr 22 Python
浅谈Python类里的__init__方法函数,Python类的构造函数
Dec 10 Python
python实现简单神经网络算法
Mar 10 Python
Python实现的文轩网爬虫完整示例
May 16 Python
win8.1安装Python 2.7版环境图文详解
Jul 01 Python
Django缓存系统实现过程解析
Aug 02 Python
python 有效的括号的实现代码示例
Nov 11 Python
Python获取对象属性的几种方式小结
Mar 12 Python
Python Flask上下文管理机制实例解析
Mar 16 Python
python安装后的目录在哪里
Jun 21 Python
python 获取计算机的网卡信息
Feb 18 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字符串处理的10个简单方法
2010/06/30 PHP
PHP详解ASCII码对照表与字符转换
2011/12/05 PHP
PHP的几个常用数字判断函数代码
2012/04/24 PHP
PHP编程中的常见漏洞和代码实例
2014/08/06 PHP
php中最简单的字符串匹配算法
2014/12/16 PHP
php微信高级接口调用方法(自定义菜单接口、客服接口、二维码)
2016/11/28 PHP
PHP编程快速实现数组去重的方法详解
2017/07/22 PHP
jquery 经典动画菜单效果代码
2010/01/26 Javascript
一个简单的js渐显(fadeIn)渐隐(fadeOut)类
2010/06/19 Javascript
js获取select标签的值且兼容IE与firefox
2013/12/30 Javascript
jQuery中:first选择器用法实例
2014/12/30 Javascript
jquery+ajax实现跨域请求的方法
2015/01/20 Javascript
限制上传文件大小和格式的jQuery插件实例
2015/01/24 Javascript
angularJS中router的使用指南
2015/02/09 Javascript
jQuery实现类似标签风格的导航菜单效果代码
2015/08/25 Javascript
JavaScript检测原始值、引用值、属性
2016/06/20 Javascript
ES6中Symbol类型用法实例详解
2017/04/06 Javascript
vuejs父子组件之间数据交互详解
2017/08/09 Javascript
带你了解session和cookie作用原理区别和用法
2017/08/14 Javascript
vue 全选与反选的实现方法(无Bug 新手看过来)
2018/02/09 Javascript
Vue2.0中集成UEditor富文本编辑器的方法
2018/03/03 Javascript
jQuery中$原理实例分析
2018/08/13 jQuery
jQuery easyui datagird编辑行删除行功能的实现代码
2018/09/20 jQuery
VUE table表格动态添加一列数据,新增的这些数据不可以编辑(v-model绑定的数据不能实时更新)
2020/04/03 Javascript
python队列通信:rabbitMQ的使用(实例讲解)
2017/12/22 Python
python中的常量和变量代码详解
2018/07/25 Python
python 对给定可迭代集合统计出现频率,并排序的方法
2018/10/18 Python
Python实现平行坐标图的绘制(plotly)方式
2019/11/22 Python
美国首屈一指的礼品篮供应商:GiftTree
2018/01/06 全球购物
Holland & Barrett爱尔兰:英国领先的健康零售商
2019/03/31 全球购物
Rowdy Gentleman服装和配饰:美好时光
2019/09/24 全球购物
给老婆的搞笑检讨书
2014/01/12 职场文书
求职信需要的五点内容
2014/02/01 职场文书
元旦文艺汇演主持词
2014/03/26 职场文书
2015年办公室人员工作总结
2015/05/15 职场文书
Python django中如何使用restful框架
2021/06/23 Python