Python -高斯卷积

Question

我需要用以3934.8A为中心的特定参数的高斯函数来卷积下一条曲线。

​

​

我看到的问题是，我的曲线是一个离散的数组，而高斯函数是一个定义良好的连续函数。我怎么才能让它工作呢？

Answer 1

为此，您需要以与曲线相同的空间比例离散化高斯曲线，然后进行卷积。

具体地说，假设原始曲线具有沿x轴均匀分布的N点(其中N通常在50到10,000之间)。那么x轴上的点间距将是(physical range)/(digital range) = (3940-3930)/N，代码将如下所示：

dx = float(3940-3930)/N
gx = np.arange(-3*sigma, 3*sigma, dx)
gaussian = np.exp(-(x/sigma)**2/2)
result = np.convolve(original_curve, gaussian, mode="full")

在这里，这是一个以零为中心的高斯函数，不包括你引用的偏移量(对我来说，这只会增加混淆，因为卷积本质上是一个翻译操作，所以从已经翻译的东西开始是令人困惑的)。

我强烈建议将所有东西都放在真实的物理单元中，就像我上面所做的那样。然后很清楚，例如，高斯的宽度是多少，等等。