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

更正后的效果

相关推荐
飞Link14 分钟前
什么是 PX4?一文了解 PX4 在空地协同中的作用
算法·无人机
一去不复返的通信er7 小时前
力扣-图论-岛屿数量
算法·leetcode·图论
wangchunyu1149 小时前
JAVA CompletableFuture异步编程实战
java·开发语言·python
程序 代码狂人10 小时前
切片表与分区表的区别(数据开发面试版)
java·大数据·开发语言·数据仓库·算法
HDU-Dade10 小时前
【校招笔试真题】简要介绍 Self-Attention 的前沿算法
人工智能·算法
2501_9445215910 小时前
Immer 跨端三方库仓颉适配全流程
python·eclipse·django
信誓旦旦的程序猿10 小时前
【Python量化实战宝典 #26】止损与回撤控制手册:移动止损、最大回撤(Python 可跑)
python·1024程序员节·股票数据api·股票数据api接口·股票api数据接口·股票量化数据api·股票量化数据接口
2501_9445215910 小时前
Deepmerge 跨端三方库仓颉适配全流程
java·python·eclipse·github
Patrick在香港11 小时前
同一行 strftime,两台机器跑出两种日期
python·datetime·时区·locale·日期处理
2601_9679873711 小时前
用Python+OpenCV+YOLO写个游戏脚本:从自动瞄准到自动采集的保姆级教程
python·opencv·yolo·目标检测·游戏自动化