所以我想知道如何使卷积和解卷积正常工作。我的代码是:import numpy as np
from scipy import interpolate
from scipy import signal
import matplotlib.pyplot as plt
a = np.arange(0, 50, 1)
b = np.array([1,2,1])
aconvolved = signal.convolve(np.sin(a), b)
adeconvolved = signal.deconvolve(aconvolved, b)
plt.figure()
plt.plot(a, np.sin(a), 'g')
plt.plot(a, aconvolved[0], 'r')
plt.plot(a, adeconvolved[0], 'b')
plt.draw()
当我把“a”卷积成一个点(比如b=[2])而不是数组时,它工作得很好。但是当我试图在一个以上的点或者一个数组上做这个,它就是不起作用。我收到的错误包括:
^{pr2}$
如果我改为使用:a = np.arange(1, 50, 1)
b = np.array([1,2,1])
即使在这种情况下,我知道它们都是从1开始的。如果我做同样的事,但是b = np.array([1,2,1])
我得到一个错误:"ValueError: BUG: filter coefficient a[0] == 0 not supported yet"
有谁知道我做错了什么吗?在