def test_lfilter_gen(self):
x = np.random.normal(0, 1, 1000)
hb = np.array([0, 0, 1, 0], dtype=np.float)
f = signal.lfilter_gen(hb, 1)
y = [f.send(v) for v in x]
self.assertArrayEqual(np.append([0, 0], x[:-2]), y)
hb, ha = sp.iirfilter(4, 0.01, btype='lowpass')
y1 = sp.lfilter(hb, ha, x)
f = signal.lfilter_gen(hb, ha)
y2 = [f.send(v) for v in x]
self.assertArrayEqual(y1, y2, precision=6)
评论列表
文章目录