高斯牛顿相关,为了说明上一篇ICP NDT精度问题

伟大的大数学家泰勒

泰勒公式,是一个用函数在某点的信息描述其附近取值的公式。如果函数满足一定的条件,泰勒公式可以用函数在某一点的各阶导数值做系数构建一个多项式来近似表达这个函数。
它最重要的贡献就是将复杂问题简单化,将非线性问题转换为线性问题,这样我们对于线性问题的多种求解方式便可以大显身手,这是一种将未知问题转换为已知问题又不丢失太大精度的思想!

超定方程的求解

当问题的已知条件/约束/等式量达不到代求解问题维数的时候,自然无法求解,只有两把钥匙却要打开多把锁?
现在我们来把锁全部打开
首先,我们将问题记为Ax=b,其中x和b都是多维(n维)的,我们约束只有极少维度(m,m<n)
这个世界上有千千万万把钥匙,现在我们拿来多把可能会打开对应锁的钥匙x1
这个时候开锁的程度记为r,r=b-Ax1,我们只要找到最可能打开的几把钥匙就好了(极为相似的钥匙也可以打开同一把锁——配钥匙就是找近似解)
配钥匙的方法就是最小二乘,即r²的导数为0,此时构建了f(x)=(Ax1-b)²的问题,求导得A高斯牛顿相关,为了说明上一篇ICP NDT精度问题(Ax-b)=0,可得高斯牛顿相关,为了说明上一篇ICP NDT精度问题即可解得高斯牛顿相关,为了说明上一篇ICP NDT精度问题

更优解、更精确解的方式

牛顿法

首先泰勒展开,忽略二阶以上的项
高斯牛顿相关,为了说明上一篇ICP NDT精度问题
然后令导数等于0
高斯牛顿相关,为了说明上一篇ICP NDT精度问题
其中多维一阶导为称为雅克比矩阵记为J,多维二阶导记为H,得到高斯牛顿相关,为了说明上一篇ICP NDT精度问题

高斯牛顿法

优化问题高斯牛顿相关,为了说明上一篇ICP NDT精度问题

增量最优高斯牛顿相关,为了说明上一篇ICP NDT精度问题

一阶泰勒近似高斯牛顿相关,为了说明上一篇ICP NDT精度问题

导数为零高斯牛顿相关,为了说明上一篇ICP NDT精度问题

解得高斯牛顿相关,为了说明上一篇ICP NDT精度问题

高斯牛顿相关,为了说明上一篇ICP NDT精度问题则有高斯牛顿相关,为了说明上一篇ICP NDT精度问题

伟大的拉格朗日

引入一个新的参数 λ (即拉格朗日乘子),将约束条件函数与原函数联系到一起
将约束条件乘以λ加入原函数中,在约束条件不满足的情况下是不可能达到最优的,这种扰动可以使求解问题更快收敛。
题外话:另外对于有k个约束问题的最优化求解,拉格朗日乘子的加入可以增加方程数量以减少不定数量使方程更易求解,是一种过拟合更快求最优的思想

最常用的求精确解方式方法Levenberg-Marquardt

带约束的求解问题高斯牛顿相关,为了说明上一篇ICP NDT精度问题

构建拉格朗日乘子方程高斯牛顿相关,为了说明上一篇ICP NDT精度问题

同样的一阶近似然后对待求量x求导,令其为零
最终得到高斯牛顿相关,为了说明上一篇ICP NDT精度问题

图解一阶牛顿法

拟牛顿法是用Hessian矩阵的逆矩阵来代替Hessian矩阵
一阶牛顿法(一维实根f(x)=0):高斯牛顿相关,为了说明上一篇ICP NDT精度问题
以上的图为实根大多数还是非实根的,上述是最简单最简单的牛顿法原理,方便理解,各参数量对应的图形用色彩标记过了,正方向也标记了
一阶泰勒展开近似高斯牛顿相关,为了说明上一篇ICP NDT精度问题
求导令其为零即可得到高斯牛顿相关,为了说明上一篇ICP NDT精度问题
高斯牛顿相关,为了说明上一篇ICP NDT精度问题
二阶导记为H 海塞(Hessian)矩阵高斯牛顿相关,为了说明上一篇ICP NDT精度问题

最后一个问题拟牛顿法

通过对一阶泰勒近似的求导高斯牛顿相关,为了说明上一篇ICP NDT精度问题
得到高斯牛顿相关,为了说明上一篇ICP NDT精度问题,我们不要将其等于0
然后有高斯牛顿相关,为了说明上一篇ICP NDT精度问题
学术上我们将高斯牛顿相关,为了说明上一篇ICP NDT精度问题这个式子叫做拟牛顿条件
要保证Gk为正定矩阵,这是因为只有正定才能保证牛顿法的搜索方向是向下搜索的
然后我们将H矩阵替代为Gk就好了,这个替代的过程是大名鼎鼎的各种算法的由来

记号的改写

这里我们记高斯牛顿相关,为了说明上一篇ICP NDT精度问题
其中高斯牛顿相关,为了说明上一篇ICP NDT精度问题 就是一阶导

另外我们将所有的替代海塞矩阵的G记为B

DFP算法

它是Davidon、Fletcher、Powell三位牛人名字的首字母缩写
DFP替代高斯牛顿相关,为了说明上一篇ICP NDT精度问题

BFGS算法(目前认为最好的替代算法,也是最常用的)

BFGS替代高斯牛顿相关,为了说明上一篇ICP NDT精度问题

上一篇:ICP-MS特点及应用领域


下一篇:备案域名为什么会掉