我有一个数值计算的光脉冲的时间形状,存储在NumPy数组中。我想找到它的高斯包络函数的参数(特别是FWHM),但首先我必须适应它。
是否有内置的SciPy函数或函数集用于信封拟合,或者我必须自己编写此代码? (由于准确性不是很重要,我可能只需要采用局部最大值并在其上拟合高斯...)
光脉冲看起来像这样:
答案 0 :(得分:1)
如果你只想找到最大值,我建议使用Scipy的argrelextrema。
但在我看来,你会想要使用Scipy的Hilbert transform 为了找到信封的分析表格。然后,您可以为包含高斯的信封生成一个数组。
对于高斯拟合,有一个很好的答案here。