简短摘要:如何快速计算两个阵列的有限卷积?我试图得到两个函数f(x),g(x)的有限卷积。为了实现这一点,我获取了函数的离散示例,并将它们转换为长度为steps的数组。steps)]garray = [g(x) for x in xarray]
然后,我尝试用scipy.signal.convolve函数计算卷积这里,Riemann表示一个简单的Riemann和,trapezoidal是一个改进的Riemann算法,它使用梯形规则,s
为了更好地理解离散有限卷积是如何工作的(阅读教育目的),我编写了卷积函数的全python实现。其想法是提供与numpy.convolve相同的输出,包括mode选项。我把代码推广到n维卷积,而不仅仅是一维卷积。这意味着它更接近于scipy.signal.convolve的行为。def py_nd_convolve(s, k, mode='full'):
# All python implementation of n-dimensional