matlab光斑中心,Matlab学习手记——二维高斯曲面拟合法定位光斑中心
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学习手记——二维高斯曲面拟合法定位光斑中心相关推荐
- Matlab绘制三维曲面(以二维高斯函数为例)
原文地址为: Matlab绘制三维曲面(以二维高斯函数为例) 寒假学习了一下Python下的NumPy和pymatlab,感觉不是很容易上手.来学校之后,决定继续看完数字图像处理一书.还是想按照上学期 ...
- 二维高斯曲面拟合法求取光斑中心及算法的C++实现
(1)二维高斯去曲面拟合推导 一个二维高斯方程可以写成如下形式: 其中,G为高斯分布的幅值,,为x,y方向上的标准差,对式(1)两边取对数,并展开平方项,整理后为: 假如参与拟合的数据点有N个,则将这 ...
- 二维高斯曲面拟合法求取光斑中心
(1)二维高斯去曲面拟合推导 一个二维高斯方程可以写成如下形式: 其中,G为高斯分布的幅值,,为x,y方向上的标准差,对式(1)两边取对数,并展开平方项,整理后为: 假如参与拟合的数据点有N个,则将这 ...
- matlab画三维向量的曲面,Matlab绘制三维曲面(以二维高斯函数为例)
寒假学习了一下Python下的NumPy和pymatlab,感觉不是很容易上手.来学校之后,决定继续看完数字图像处理一书.还是想按照上学期的模式,边看边实现书中的算法.上学期看的时候,是用C语言实现的 ...
- MATLAB编程(4)——MATLAB绘制二维高斯函数的三维图
本篇博文记录使用MATLAB绘制二维高斯函数的三维图. 用到的MATLAB函数--mesh()(绘制三维线框图)和surf()(绘制三维表面图). MATLAB命令窗口输入>> doc 函 ...
- Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动
Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...
- 高斯 二维 matlab,MATLAB实例:二维高斯分布图
MATLAB实例:二维高斯分布图 问题描述:用MATLAB在一张图上同时画出二维高斯散点图与分布图. 数据来源:MATLAB高斯混合数据的生成 将二维数据保存为gauss_data.txt. 前期工作 ...
- matlab生成满足二维高斯(正…
原文地址:matlab生成满足二维高斯(正态)分布的随机数/作图程序作者:乐韵悠杨 产生满足二维高斯(正态)分布的随机数: mu=[0,2];%数学期望 sigma=[1 0;0,4];%协方差矩阵 ...
- 学习:二维码、QR码、J4L-QRCode、java
转自:http://baijinshan.iteye.com/blog/1004554 学习:二维码.QR码.J4L-QRCode.java Java企业应用网络应用Rails应用服务器 开源码 J ...
- 深度强化学习笔记(二)——Q-learning学习与二维寻路demo实现
深度强化学习笔记(二)--Q-learning学习与二维寻路demo实现 文章目录 深度强化学习笔记(二)--Q-learning学习与二维寻路demo实现 前言 理论 什么是Q-Learning 算 ...
最新文章
- 吴恩达深度学习笔记(67)-迁移学习(Transfer learning)
- 程序员经常去的 14 个顶级开发者社区(转)
- 一句话实现Mysql查询结果带行号
- SAP WebIDE里的JavaScript代码检查code check
- linux网络编程之inet_pton和inet_ntop函数
- 水晶报表 jar包版本过低_工具类学习-UReport报表设计器整合
- 七:构造方法与自定义构造方法
- 随便创建Servlet文件都出现,请求资源不可用的情况
- 设计模式入门-工厂模式
- 如何以管理员方式管理文件
- Android 学习思维导图
- strcmp函数的实现
- 当toB产品经理遇到数据分析
- Mac系统搭建C语言开发环境
- 身居乱世之中,重新审视“活法
- 35岁的大龄程序员都去哪里了?
- STIM300读取数据
- create Golang project with interface #1
- 链安团队漏洞分析连载(第一期)一一溢出漏洞
- 深度学习 -- SSD 算法流程详解