Matlab学习手记——二维高斯曲面拟合法定位光斑中心

Matlab学习手记——二维高斯曲面拟合法定位光斑中心

终端探测器探测到的光斑灰度分布函数可近似看做高斯分布,因此可以通过二维高斯函数进行拟合,模型表示为:

效果图:

为了方便计算,做一步变换,两边取对数,得到:

展开并进一步变形为:

求解问题变为二元二次多项式拟合问题。

此问题一般基于最小二乘来求解,这里提供两种思路,第一种思路是通用方法,第二种是针对此类问题的特殊解法,更为简单。

1)基于梯度

最小二乘求解,使残差平方和最小:

根据最小值条件,对每个变量求偏导并令其为零,得到方程组:

进一步变换:

得到线性方程组:

求解此线性方程组得到高斯函数的各项参数:

2)二次多项式拟合的简单方法

可以参考另外两篇博客:

1. Matlab学习手记——二次多项式曲面拟合

2. Matlab学习手记——基于最小二乘的多项式曲面拟合

确定待拟合的x、y和z,写出系数矩阵即可求解。

xfit2 = xfit .* xfit;

yfit2 = yfit .* yfit;

A = [xfit2, yfit2, xfit, yfit, ones(num, 1)];

result = A \ zfit;

Matlab学习手记——二维高斯曲面拟合法定位光斑中心相关教程

matlab光斑中心,Matlab学习手记——二维高斯曲面拟合法定位光斑中心相关推荐

  1. Matlab绘制三维曲面(以二维高斯函数为例)

    原文地址为: Matlab绘制三维曲面(以二维高斯函数为例) 寒假学习了一下Python下的NumPy和pymatlab,感觉不是很容易上手.来学校之后,决定继续看完数字图像处理一书.还是想按照上学期 ...

  2. 二维高斯曲面拟合法求取光斑中心及算法的C++实现

    (1)二维高斯去曲面拟合推导 一个二维高斯方程可以写成如下形式: 其中,G为高斯分布的幅值,,为x,y方向上的标准差,对式(1)两边取对数,并展开平方项,整理后为: 假如参与拟合的数据点有N个,则将这 ...

  3. 二维高斯曲面拟合法求取光斑中心

    (1)二维高斯去曲面拟合推导 一个二维高斯方程可以写成如下形式: 其中,G为高斯分布的幅值,,为x,y方向上的标准差,对式(1)两边取对数,并展开平方项,整理后为: 假如参与拟合的数据点有N个,则将这 ...

  4. matlab画三维向量的曲面,Matlab绘制三维曲面(以二维高斯函数为例)

    寒假学习了一下Python下的NumPy和pymatlab,感觉不是很容易上手.来学校之后,决定继续看完数字图像处理一书.还是想按照上学期的模式,边看边实现书中的算法.上学期看的时候,是用C语言实现的 ...

  5. MATLAB编程(4)——MATLAB绘制二维高斯函数的三维图

    本篇博文记录使用MATLAB绘制二维高斯函数的三维图. 用到的MATLAB函数--mesh()(绘制三维线框图)和surf()(绘制三维表面图). MATLAB命令窗口输入>> doc 函 ...

  6. Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动

    Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...

  7. 高斯 二维 matlab,MATLAB实例:二维高斯分布图

    MATLAB实例:二维高斯分布图 问题描述:用MATLAB在一张图上同时画出二维高斯散点图与分布图. 数据来源:MATLAB高斯混合数据的生成 将二维数据保存为gauss_data.txt. 前期工作 ...

  8. matlab生成满足二维高斯(正…

    原文地址:matlab生成满足二维高斯(正态)分布的随机数/作图程序作者:乐韵悠杨 产生满足二维高斯(正态)分布的随机数: mu=[0,2];%数学期望 sigma=[1 0;0,4];%协方差矩阵 ...

  9. 学习:二维码、QR码、J4L-QRCode、java

    转自:http://baijinshan.iteye.com/blog/1004554 学习:二维码.QR码.J4L-QRCode.java Java企业应用网络应用Rails应用服务器  开源码 J ...

  10. 深度强化学习笔记(二)——Q-learning学习与二维寻路demo实现

    深度强化学习笔记(二)--Q-learning学习与二维寻路demo实现 文章目录 深度强化学习笔记(二)--Q-learning学习与二维寻路demo实现 前言 理论 什么是Q-Learning 算 ...

最新文章

  1. 吴恩达深度学习笔记(67)-迁移学习(Transfer learning)
  2. 程序员经常去的 14 个顶级开发者社区(转)
  3. 一句话实现Mysql查询结果带行号
  4. SAP WebIDE里的JavaScript代码检查code check
  5. linux网络编程之inet_pton和inet_ntop函数
  6. 水晶报表 jar包版本过低_工具类学习-UReport报表设计器整合
  7. 七:构造方法与自定义构造方法
  8. 随便创建Servlet文件都出现,请求资源不可用的情况
  9. 设计模式入门-工厂模式
  10. 如何以管理员方式管理文件
  11. Android 学习思维导图
  12. strcmp函数的实现
  13. 当toB产品经理遇到数据分析
  14. Mac系统搭建C语言开发环境
  15. 身居乱世之中,重新审视“活法
  16. 35岁的大龄程序员都去哪里了?
  17. STIM300读取数据
  18. create Golang project with interface #1
  19. 链安团队漏洞分析连载(第一期)一一溢出漏洞
  20. 深度学习 -- SSD 算法流程详解

热门文章

  1. dmp文件导入 Oracle 步骤
  2. windows10自带我的文档等路径修改
  3. Google谷歌新手SEO优化教程篇【1】
  4. Java工程师的职业规划,上个月成功拿到阿里P7offer
  5. 启动Kylin出现的问题
  6. 2021-2027全球与中国筒式过滤器外壳市场现状及未来发展趋势
  7. android 强制关闭键盘,Android关闭输入软键盘无效的问题
  8. 锁定计算机后怎么解锁,电脑锁住了怎么解锁
  9. RuoYi(若依) 微服务分离版 启动及常见问题总结
  10. 计算机组成原理区分正负数实验,计算机组成原理实验报告(4个).doc