测试《A Simple Algorithm for Fitting a Gaussian Function》拟合

原创已于 2026-01-02 18:39:35 修改 · 213 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#算法 #python

于 2025-12-27 22:26:08 首次发布

https://github.com/JohannesMeyersGit/1D-Gaussian-Fitting/blob/main/Itterativ_1D_Gaussian_Fit.py

源码每次迭代采样不同子区间，error(拟合的均值-实际均值) 先减低后增，改成样本点不变

error 曲线看上去正常,看上去 A 的值离实际很大，不知道是不是和论文的结论一致？

补充:github上代码可能是作者公式写错了

def get_gauss_params(a, b, c):
    """
    Calculate gaussian coefficients from parabola coefficients.
    See Eqs. (5-7)
    """
    mu = -b / (2 * c)
    sigma = np.sqrt((-1 / (2 * c)))
    #公式写错了
    #A = np.exp(((a - b ** 2) / 4 * c))
    A = np.exp(a - (b ** 2) / (4 * c))
    return A, mu, sigmax

更正后的效果