Question

我正在尝试在Python中设置一个过滤器，我陷入了一个基本的问题。

这是我必须实现的功能：

y[n] = (1 − λ) · x[n] + λ · y[n − 1]

我正在读取m文件并应用过滤器。

以下是代码：

import sys
import numpy as np
import matplotlib.pyplot as plt

y = np.loadtxt('acs712_192us.m')
size = len(y)

x = np.arange(0, size)

out = []

lamb = 0.9 

for i in range(0, len(y)):
    out.append(((1-lamb) * y[i]) + (lamb * out[i - 1]))

plt.plot(x, y)
plt.plot(x, out)
plt.show()

运行此程序时，出现以下错误：

文件“main_LI.py”，第15行，在out.append（（（1-lamb）* y [i]）+（lamb * out [i-1]））IndexError：列表索引超出范围< / p>

我知道我无法访问[-1]位置（在第一个循环中），这可能是问题吗？

任何人都有任何提示吗？

谢谢你们！

Answer 1

首先，您可以在Python中使用negative indexes。问题是您提供的索引不存在。

out最初是一个空列表，所以问题是你试图访问第一次迭代中不存在的索引：

import sys
import numpy as np
import matplotlib.pyplot as plt

y = np.loadtxt('acs712_192us.m')
size = len(y)

x = np.arange(0, size)

out = []

lamb = 0.9 

for i in range(len(y)):
    if i == 0:
       out.append(value) # value is what you want out[0] to be
    else:
        out.append(((1-lamb) * y[i]) + (lamb * out[i - 1]))

plt.plot(x, y)
plt.plot(x, out)
plt.show()

Answer 2

您可以使用列表推导来轻松完成此任务。

l = 0.9
x = np.arange(0, 10)
filtered = [x[i] if i == 0 else (1-l)*x[i] + l*x[i-1] for i in range(0,len(x))]

>>> print(filtered)
[0, 0.09999999999999998, 1.1, 2.1, 3.1, 4.1, 5.1, 6.1000000000000005, 7.1, 8.1]

这是你期望的输出吗？

在python中访问数组的位置[i-1] - 过滤

2 个答案: