Question

我正在尝试对数据使用低通滤波器。我的样本数据是每分钟的数据，我想用1440（每天）的时间常数对其进行过滤。它只是编写一个for循环，用于计算过滤并将数据附加到y（已过滤数据）列表中。

以下是我编写的代码：

OD_deltaT=df["OD_deltaT"].values.tolist()
#a=1/(timeconstant)+1
a=1/1441
u, y=[0], []
for i, x in enumerate(OD_deltaT,1):
    u.append(x)
    y[1]=u[1]
    y[i]=(1-a)*y[i-1]+(a)*u[i]

但我收到此错误：

list index out of range.

提前谢谢！

Answer 1

您需要附加数据，就像对u所做的那样。您还可以以更简洁的方式编写代码：

u = [0]
y = [0]
for x in OD_deltaT:
    u.append(x)
    y.append((1-a)*y[-1])+a*x)

当然，scipy过滤器可以更快地完成;）

低通滤波器-Python

1 个答案: