相关疑难解决方法(0)

我试图用python过滤嘈杂的心率信号.因为心率不应该是每分钟约220次心跳,我想过滤掉220bpm以上的所有噪音.我将220 /分钟转换为3.66666666赫兹,然后将赫兹转换为rad/s得到23.0383461 rad/sec.

采集数据的芯片采样频率为30Hz,因此我将其转换为rad/s,得到188.495559 rad/s.

在线查找了一些东西之后,我发现了一些带通滤波器的功能,我想把它变成低通.这是带通代码的链接,所以我将其转换为:

from scipy.signal import butter, lfilter
from scipy.signal import freqs

def butter_lowpass(cutOff, fs, order=5):
    nyq = 0.5 * fs
    normalCutoff = cutOff / nyq
    b, a = butter(order, normalCutoff, btype='low', analog = True)
    return b, a

def butter_lowpass_filter(data, cutOff, fs, order=4):
    b, a = butter_lowpass(cutOff, fs, order=order)
    y = lfilter(b, a, data)
    return y

cutOff = 23.1 #cutoff frequency in rad/s
fs = 188.495559 #sampling frequency in rad/s
order = 20 #order of …

我正在尝试使用带有128点Hamming窗口的带通滤波器,其截止频率为0.7-4Hz(python).我从图像中获取样本信号.(1个样本= 1个图像).fps经常变化.

怎么能在python中完成？我读到这个:http://mpastell.com/2010/01/18/fir-with-scipy/但我发现firwin相当令人困惑.如何使用此变量fps完成此操作？

我正在阅读一篇论文，试图重现论文的结果。在本文中，他们对原始数据使用低通切比雪夫 I 型滤波器。他们给出了这些参数。

采样频率= 32Hz，Fcut=0.25Hz，Apass=0.001dB，Astop=-100dB，Fstop=2Hz，滤波器阶数=5。我找到了一些材料帮助我理解这些参数

但是当我看一下 scipy.signal.cheby1 时。该函数所需的参数不同。

cheby1(N, rp, Wn, btype='low', analog=False, output='ba')

这里N：滤波器的阶数；btype：滤波器的类型，在我的例子中，它是“低通”；Analog=False，因为数据是采样的，所以是数字的；输出：指定输出的类型。但我不确定rp、Wn。

在文档中，它说：

rp ：浮点 通带中允许低于单位增益的最大纹波。以分贝为单位指定，为正数。

Wn : array_like 给出临界频率的标量或长度为 2 的序列。对于 I 型滤波器，这是过渡带中增益首次降至 -rp 以下的点。对于数字滤波器，Wn 从 0 到 1 归一化，其中 1 是奈奎斯特频率，pi 弧度/样本。（因此，Wn 的单位是半周期/样本。）对于模拟滤波器，Wn 是角频率（例如rad/s）。

根据这个问题： How To apply a filter to a signal in python

我知道如何使用过滤器。但我不知道如何创建一个具有与上述相同参数的过滤器。我不知道如何转换这些参数并将它们提供给Scipy中的函数。