Python基于scipy实现信号滤波功能


Posted in Python onMay 08, 2019

Python基于scipy实现信号滤波功能Python基于scipy实现信号滤波功能

1.背景介绍

在深度学习中,有时会使用Matlab进行滤波处理,再将处理过的数据送入神经网络中。这样是一般的处理方法,但是处理起来却有些繁琐,并且有时系统难以运行Matlab。Python作为一种十分强大的语言,是支持信号滤波滤波处理的。

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

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

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

2.实战演练

首先我们使用到了scipy模块,可以通过下述命令进行安装:(我使用的Python==3.6)

pip install scipy

1).低通滤波

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

from scipy import signal
b, a = signal.butter(8, 0.8, 'lowpass')  #配置滤波器 8 表示滤波器的阶数
filtedData = signal.filtfilt(b, a, data) #data为要过滤的信号

2).高通滤波

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

from scipy import signal
b, a = signal.butter(8, 0.2, 'highpass')  #配置滤波器 8 表示滤波器的阶数
filtedData = signal.filtfilt(b, a, data) #data为要过滤的信号

3).带通滤波

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

from scipy import signal
b, a = signal.butter(8, [0.2,0.8], 'bandpass')  #配置滤波器 8 表示滤波器的阶数
filtedData = signal.filtfilt(b, a, data) #data为要过滤的信号

4).带阻滤波

这里假设采样频率为1000hz,信号本身最大的频率为500hz,要滤除100hz以上,400hz以下频率成分,即截至频率为100,400hz,则wn1=2*100/1000=0.2,Wn1=0.2; wn2=2*400/1000=0.8,Wn2=0.8。Wn=[0.02,0.8],和带通相似,但是带通是保留中间,而带阻是去除。

from scipy import signal
b, a = signal.butter(8, [0.2,0.8], 'bandstop')  #配置滤波器 8 表示滤波器的阶数
filtedData = signal.filtfilt(b, a, data) #data为要过滤的信号

3.函数介绍

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'

总结

以上所述是小编给大家介绍的Python基于scipy实现信号滤波功能,希望对大家有所帮助,如果大家有任何疑问欢迎给我留言,小编会及时回复大家的!

Python 相关文章推荐
用Python编写分析Python程序性能的工具的教程
Apr 01 Python
python类中super()和__init__()的区别
Oct 18 Python
用python找出那些被“标记”的照片
Apr 20 Python
老生常谈Python startswith()函数与endswith函数
Sep 08 Python
对python中的six.moves模块的下载函数urlretrieve详解
Dec 19 Python
基于Python实现用户管理系统
Feb 26 Python
Python实现的统计文章单词次数功能示例
Jul 08 Python
对django views中 request, response的常用操作详解
Jul 17 Python
numpy:找到指定元素的索引示例
Nov 26 Python
JupyterNotebook 输出窗口的显示效果调整方法
Apr 13 Python
python操作toml文件的示例代码
Nov 27 Python
python 利用matplotlib在3D空间中绘制平面的案例
Feb 06 Python
python实现抽奖小程序
Apr 15 #Python
Python常见数据类型转换操作示例
May 08 #Python
Python数据类型之Number数字操作实例详解
May 08 #Python
利用PyCharm Profile分析异步爬虫效率详解
May 08 #Python
Python数据类型之String字符串实例详解
May 08 #Python
Python数据类型之List列表实例详解
May 08 #Python
Python3使用TCP编写一个简易的文件下载器功能
May 08 #Python
You might like
SONY ICF-SW55的电路分析
2021/03/02 无线电
php学习之简单计算器实现代码
2011/06/09 PHP
PHP动态创建Web站点的方法
2011/08/14 PHP
从PHP $_SERVER相关参数判断是否支持Rewrite模块
2013/09/26 PHP
利用switch语句进行多选一判断的实例代码
2016/11/14 PHP
thinkPHP多表查询及分页功能实现方法示例
2017/07/03 PHP
laravel框架使用极光推送消息操作示例
2020/02/15 PHP
一些常用的Javascript函数
2006/12/22 Javascript
使用JavaScript 编写简单计算器
2014/11/24 Javascript
JavaScript驾驭网页-获取网页元素
2016/03/24 Javascript
Bootstrap项目实战之首页内容介绍(全)
2016/04/25 Javascript
javascript中利用柯里化函数实现bind方法【推荐】
2016/04/29 Javascript
基于HTML+CSS+JS实现增加删除修改tab导航特效代码
2016/08/05 Javascript
基于Phantomjs生成PDF的实现方法
2016/11/07 Javascript
前端编码规范(3)JavaScript 开发规范
2017/01/21 Javascript
javascript编程实现栈的方法详解【经典数据结构】
2017/04/11 Javascript
详解Angular5 服务端渲染实战
2018/01/04 Javascript
node.js 基于cheerio的爬虫工具的实现(需要登录权限的爬虫工具)
2019/04/10 Javascript
layui的数据表格+springmvc实现搜索功能的例子
2019/09/28 Javascript
[02:57]DOTA2亚洲邀请赛小组赛第四日 赛事回顾
2015/02/02 DOTA
[04:09]2018年度DOTA2社区贡献奖-完美盛典
2018/12/16 DOTA
Python实现随机生成有效手机号码及身份证功能示例
2017/06/05 Python
Java实现的执行python脚本工具类示例【使用jython.jar】
2018/03/29 Python
Python决策树之基于信息增益的特征选择示例
2018/06/25 Python
利用html5 canvas破解简单验证码及getImageData接口应用
2013/01/25 HTML / CSS
日本最佳原创设计品牌:Felissimo(芬理希梦)
2019/03/19 全球购物
软件工程师面试题
2012/06/25 面试题
4s店机修工岗位职责
2013/12/20 职场文书
企业文明单位申报材料
2014/05/16 职场文书
初中学校对照检查材料
2014/08/19 职场文书
房地产营销活动策划方案
2014/09/15 职场文书
大学生入党自传2015
2015/06/26 职场文书
升学宴来宾致辞
2015/07/27 职场文书
2015秋季田径运动会广播稿
2015/08/19 职场文书
在Java中Collection的一些常用方法总结
2021/06/13 Java/Android
HTML5中的DOCUMENT.VISIBILITYSTATE属性详解
2023/05/07 HTML / CSS