deconvolve(signal, divisor) function of
> Deconvolves `divisor` out of `signal`. Returns the quotient and remainder such that `signal = convolve(divisor, quotient) + remainder`
The core codes of
num = atleast_1d(signal) den = atleast_1d(divisor) N = len(num) D = len(den) input = ones(N - D + 1, float) input[1:] = 0 quot = lfilter(num, den, input)
It's quite different from any deconvolve methods I ever know, such as the frequency spectrum division.
Can anyone explain those behind the codes?