Python SciPy signal.deconvolve用法及代码示例

本文简要介绍 python 语言中 scipy.signal.deconvolve 的用法。

用法:

scipy.signal.deconvolve(signal, divisor)#

使用逆滤波从signal 中解卷积divisor。

返回商和余数使得signal = convolve(divisor, quotient) + remainder

参数 ：：

signal： (N,) 数组

信号数据，通常是记录的信号

divisor： (N,) 数组

除数数据，通常是应用于原始信号的脉冲响应或滤波器

返回 ：：

quotient： ndarray

商，通常是恢复的原始信号

remainder： ndarray

余

例子：

对已过滤的信号进行反卷积：

>>> from scipy import signal
>>> original = [0, 1, 0, 0, 1, 1, 0, 0]
>>> impulse_response = [2, 1]
>>> recorded = signal.convolve(impulse_response, original)
>>> recorded
array([0, 2, 1, 0, 2, 3, 1, 0, 0])
>>> recovered, remainder = signal.deconvolve(recorded, impulse_response)
>>> recovered
array([ 0.,  1.,  0.,  0.,  1.,  1.,  0.,  0.])