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

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

更正后的效果

相关推荐
小蒜学长4 小时前
基于Python的微博热搜话题可视化分析系统的设计与实现(代码+数据库+LW)
后端·python·django·数据可视化·情感分析·微博热搜
梦帮科技4 小时前
从提示词到可复现音乐工作流:RNOISE 2.0 的 Prompt Engineering 架构
人工智能·python·mysql·ci/cd·架构·node.js·numpy
飞Link4 小时前
【Numpy】 第三部分:矩阵运算与线性代数(计算篇)
python·numpy
电商API_180079052474 小时前
电商平台数据分析实战_帖子
开发语言·爬虫·python·数据采集·京东
AI 小老六4 小时前
Agent 记忆系统难在取舍
人工智能·算法·架构·agent·memory·harness
大熊背4 小时前
《Color constancy by characterization of illumination chromaticity》之色度色域最大化算法(一)
数码相机·算法·白平衡·色度色域
不会就选b4 小时前
算法日常・每日刷题--<贪心>13
算法
带多刺的玫瑰5 小时前
Leecode#35刷题之搜索插入位置
java·python·算法
QQ14220784495 小时前
IIS + wfastcgi 部署 Flask 应用实战:从连续 500 到稳定 200 的完整踩坑实录
python
WiKiLeaks_successor5 小时前
Scipy库里的众数函数不严谨,我把它重构了。
python·scipy