python scipy.signal_使用scipy.signal在Python中进行卷积和反卷积

deconvolve返回两个数组,商和余数.所以尝试：

f, r = signal.deconvolve(s, s_f)

docstring显示了deconvolve的使用示例.这是另一个(sig是scipy.signal,np是numpy)：

要去卷积的信号是z,滤波器系数在滤波器中：

In [9]: z

Out[9]:

array([ 0.5, 2.5, 6. , 9.5, 11. , 10. , 9.5, 11.5, 10.5,

5.5, 2.5, 1. ])

In [10]: filter = np.array([0.5, 1.0, 0.5])

应用deconvolve：

In [11]: q, r = sig.deconvolve(z, filter)

In [12]: q

Out[12]: array([ 1., 3., 5., 6., 5., 4., 6., 7., 1., 2.])

将过滤器应用于q以验证我们返回z：

In [13]: sig.convolve(q, filter)

Out[13]:

array([ 0.5, 2.5, 6. , 9.5, 11. , 10. , 9.5, 11.5, 10.5,

5.5, 2.5, 1. ])

通过构造,这是一个非常干净的例子.余数为零：

In [14]: r

Out[14]: array([ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])

当然,你不会总是得到这么好的结果.