测试《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

更正后的效果

相关推荐
M78佐菲21 小时前
c语言学习笔记:排序与查找方法整理
linux·c语言·笔记·学习·算法
leisoo80971 天前
ig50数据落盘ClickHousevsTimescaleDBvsDuckDB实测对比 IG50免费开源股票数据API接口
开发语言·jvm·数据库·python·json
卷无止境1 天前
FastAPI 的 Metadata 到底是什么,又牵动了哪些核心概念
后端·python
卷无止境1 天前
FastAPI 调试实战,从断点到生产环境的排错心法
后端·python
yaoxin5211231 天前
503. Java 反射 - 编写 ServiceFactory 类
java·开发语言·python
AI备案指南-满满1 天前
人工智能拟人化互动服务安全自评估报告的评估要点有哪些?
人工智能·算法·安全·机器人·大模型备案·算法备案
kyrie_sakura1 天前
python学习笔记3 -- 流程控制语句结构
笔记·python·学习
程序员小八7771 天前
Java 快速转 Go
java·python·golang
土司大王1 天前
LeetCode hot100——两两交换链表中的节点
算法·leetcode·职场和发展