SciPy - 适合高斯包络

时间:2017-08-26 12:19:10

标签: python numpy scipy signal-processing gaussian

我有一个数值计算的光脉冲的时间形状,存储在NumPy数组中。我想找到它的高斯包络函数的参数(特别是FWHM),但首先我必须适应它。
是否有内置的SciPy函数或函数集用于信封拟合,或者我必须自己编写此代码? (由于准确性不是很重要,我可能只需要采用局部最大值并在其上拟合高斯...)

光脉冲看起来像这样:

The light pulse in question

1 个答案:

答案 0 :(得分:1)

如果你只想找到最大值,我建议使用Scipy的argrelextrema

但在我看来,你会想要使用Scipy的Hilbert transform 为了找到信封的分析表格。然后,您可以为包含高斯的信封生成一个数组。

对于高斯拟合,有一个很好的答案here