python滤波器 lfilter和filtfilt 对比与应用 基于Butterworth Filter

Code

# initial filter param  
def init_filter(self, order, cutoff, fs):
	wn = 2.0*cutoff/fs
	self._b, self._a = sig.butter(order, wn, 'lowpass', analog = False)
	self._dss = sig.lfilter_zi(self._b, self._a)

def lfilter_filter(self, f_x, f_y, f_z):

    f = [f_x, f_y, f_z]
    res = [0,0,0]

    for i in range (0,3):
        res[i], self._dss[i] = sig.lfilter(self._b, self._a, [f[i]], zi=self._dss[i])

    return res[0][0], res[1][0], res[2][0]    

Plot

init_filter(order=3, cutoff=4, fs=100)
在这里插入图片描述
init_filter(order=3, cutoff=6, fs=100)
在这里插入图片描述

Ref

一些文字:在scipy.signal中应用过滤器:使用lfilter还是filtfilt?
一些代码:scipy中 filtfilt和lfilter 区别

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值