详解利用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浅拷贝与深拷贝用法实例
May 09 Python
在Python的Django框架中更新数据库数据的方法
Jul 17 Python
实例解析Python设计模式编程之桥接模式的运用
Mar 02 Python
python3读取excel文件只提取某些行某些列的值方法
Jul 10 Python
python中多个装饰器的执行顺序详解
Oct 08 Python
Django中Middleware中的函数详解
Jul 18 Python
Form表单及django的form表单的补充
Jul 25 Python
详解Python self 参数
Aug 30 Python
Python下利用BeautifulSoup解析HTML的实现
Jan 17 Python
Python如何使用ElementTree解析xml
Oct 12 Python
python 制作一个gui界面的翻译工具
May 14 Python
Python Pandas pandas.read_sql_query函数实例用法分析
Jun 21 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中的session永不过期的解决思路及实现方法分享
2011/04/20 PHP
php入门学习知识点八 PHP中for循环基本应用之九九乘法口绝表
2011/07/14 PHP
smarty模板中使用get、post、request、cookies、session变量的方法
2014/04/24 PHP
php微信开发之自定义菜单完整流程
2016/10/08 PHP
PHP耦合设计模式实例分析
2018/08/08 PHP
js left,right,mid函数
2008/06/10 Javascript
分享一个用Mootools写的鼠标滑过进度条改变进度值的实现代码
2011/12/12 Javascript
jquery获取特定name所有选中的checkbox,支持IE9标准模式
2013/03/18 Javascript
jQuery实现瀑布流的取巧做法分享
2015/01/12 Javascript
Jquery幻灯片特效代码分享--鼠标点击按钮时切换(1)
2015/08/15 Javascript
JS实现可关闭的对联广告效果代码
2015/09/14 Javascript
Bootstrap图片轮播组件使用实例解析
2016/06/30 Javascript
javascript 判断是否是微信浏览器的方法
2016/10/09 Javascript
jquery精度计算代码 jquery指定精确小数位
2017/02/06 Javascript
微信小程序之MaterialDesign--input组件详解
2017/02/15 Javascript
Bootstrap3 多个模态对话框无法显示的解决方案
2017/02/23 Javascript
基于JavaScript实现微信抢红包功能
2017/07/20 Javascript
简单实现js鼠标跟随效果
2020/08/02 Javascript
vue-router权限控制(简单方式)
2018/10/29 Javascript
微信小程序实用代码段(收藏版)
2019/12/17 Javascript
python实现保存网页到本地示例
2014/03/16 Python
Python中pygal绘制雷达图代码分享
2017/12/07 Python
python 多线程中子线程和主线程相互通信方法
2018/11/09 Python
django 2.2和mysql使用的常见问题
2019/07/18 Python
python Dijkstra算法实现最短路径问题的方法
2019/09/19 Python
浅谈Python3多线程之间的执行顺序问题
2020/05/02 Python
python中如何进行连乘计算
2020/05/28 Python
我们在web应用开发过程中经常遇到输出某种编码的字符,如iso8859-1等,如何输出一个某种编码的字符串?
2014/03/30 面试题
团委书记的竞聘演讲稿
2014/04/24 职场文书
离婚协议书包括哪些内容
2014/10/16 职场文书
运动员入场词
2015/07/18 职场文书
生日宴会家属答谢词
2015/09/29 职场文书
机械生产实习心得体会
2016/01/22 职场文书
送给火锅店的创意营销方案!
2019/07/08 职场文书
话题作文之生命的旋律
2019/12/17 职场文书
MySQL pt-slave-restart工具的使用简介
2021/04/07 MySQL