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列表删除的三种方法代码分享
Oct 31 Python
浅谈python可视化包Bokeh
Feb 07 Python
Python实现将数据框数据写入mongodb及mysql数据库的方法
Apr 02 Python
对pandas的行列名更改与数据选择详解
Nov 12 Python
Python中查看变量的类型内存地址所占字节的大小
Jun 26 Python
Python的Lambda函数用法详解
Sep 03 Python
django框架中间件原理与用法详解
Dec 10 Python
Python:type、object、class与内置类型实例
Dec 25 Python
简单了解python调用其他脚本方法实例
Mar 26 Python
使用python-Jenkins批量创建及修改jobs操作
May 12 Python
Django 实现jwt认证的示例
Apr 30 Python
Pytorch 如何实现LSTM时间序列预测
May 17 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
PHP中的日期处理方法集锦
2007/01/02 PHP
php操作SVN版本服务器类代码
2011/11/27 PHP
PHP error_log()将错误信息写入一个文件(定义和用法)
2013/10/25 PHP
PHP处理Oracle的CLOB实例
2014/11/03 PHP
Codeigniter的dom类用法实例
2015/06/26 PHP
PHP下载远程图片的几种方法总结
2017/04/07 PHP
PHP如何根据文件头检测文件类型实例代码
2018/10/14 PHP
JavaScript的面向对象(一)
2006/11/09 Javascript
解决jquery submit()提交表单提示:f[s] is not a function
2013/01/23 Javascript
JavaScript动态插入script的基本思路及实现函数
2013/11/11 Javascript
JavaScript阻止事件冒泡示例分享
2014/12/28 Javascript
js实现背景图片感应鼠标变化的方法
2015/02/28 Javascript
jQuery仿Flash上下翻动的中英文导航菜单实例
2015/03/10 Javascript
JavaScript编写简单的计算器
2015/11/25 Javascript
jquery分页插件jquery.pagination.js使用方法解析
2016/04/01 Javascript
微信小程序 loading 详解及实例代码
2016/11/09 Javascript
使用mock.js随机数据和使用express输出json接口的实现方法
2018/01/07 Javascript
微信小程序登陆注册功能的实现代码
2019/12/10 Javascript
vue实现拖拽效果
2019/12/23 Javascript
微信浏览器下拉黑边解决方案 wScroollFix
2020/01/21 Javascript
以Flask为例讲解Python的框架的使用方法
2015/04/29 Python
Python3 循环语句(for、while、break、range等)
2017/11/20 Python
python使用代理ip访问网站的实例
2018/05/07 Python
将pip源更换到国内镜像的详细步骤
2019/04/07 Python
使用Python制作表情包实现换脸功能
2019/07/19 Python
matplotlib绘制多个子图(subplot)的方法
2019/12/03 Python
python通过nmap扫描在线设备并尝试AAA登录(实例代码)
2019/12/30 Python
Python阶乘求和的代码详解
2020/02/14 Python
pip安装tensorflow的坑的解决
2020/04/19 Python
美国汽车轮胎和轮毂销售网站:Tire Rack
2018/01/11 全球购物
世界闻名的衬衫制造商:Savile Row Company
2018/07/30 全球购物
英语感恩演讲稿
2014/01/14 职场文书
房地产公司见习自我鉴定
2014/04/28 职场文书
促销活动计划书
2014/05/02 职场文书
授权委托书怎么写
2014/09/25 职场文书
mysql查询结果实现多列拼接查询
2022/04/03 MySQL