边缘化中FEJ图例的理解

在这里插入图片描述

        如图所示,在解释为什么需要FEJ(First Estimation Jacobian)时,通常会将这个图拿出来说事。但是,很多时候只是一笔带过,这个图看的云里雾里的,不是非常明白(可能是我理解力的问题),所以,今天想专门写个笔记记录下对这张图的理解。

         非线性优化中,\Delta x取到最优的点就是:J^{T}J\Delta x =-J^{T}b,它是从F(\Delta x)的泰勒展开式的导数推导出来的,即(J\Delta x - b)^{T}(J\Delta x -b) =(\Delta x - J^{+}b)^{T}J^{T}J(\Delta x -J^{+}b),J^{+}J的广义逆。从这个式子可以看到,假如把\Delta x的分布看成正态分布的话,J^{T}J就是该正态分布的\sigma,即协方差,它决定了\Delta x的分布是一条直线还是一个椭圆。

        协方差代表了分布在主方向和非主方向上的方差,具体原理可以参考 CodingLabs - PCA的数学原理,如果协方差的通过特征值分解能变成

\begin{pmatrix} a & 0\\ 0 & b \end{pmatrix}

,a,b是该协方差的特征值,如果ab都不为0,它的分布就是个椭圆。否则有一个为0,它的分布就是直线。

        对于E_{1}^{'},在(0.5, 1.4)处展开,它的

J^{T}J=\begin{pmatrix} \frac{\partial E}{\partial x} \\ \frac{\partial E}{\partial y} \end{pmatrix} * \begin{pmatrix} \frac{\partial E}{\partial x} & \frac{\partial E}{\partial y} \end{pmatrix} = \begin{pmatrix} -0.3\\-0.84 \end{pmatrix} * \begin{pmatrix} -0.3,-0.84 \end{pmatrix}=\begin{pmatrix} 0.09&0.252\\0.252 &0.7056\end{pmatrix}

,求它的特征值:

\begin{vmatrix} 0.09-\lambda & 0.252\\0.252 &0.7056- \lambda \end{vmatrix}=(0.09-\lambda)(0.7056-\lambda)-0.252*0.252=\lambda(\lambda-0.7956)

即它的特征值是0和0.7956,那么它的分布是条直线,即图2所示的样子。

        同理,对于E_{2}^{'},在(1.2, 0.5)处展开,它的J^{T}J=\begin{pmatrix} -0.24\\-0.4 \end{pmatrix} * \begin{pmatrix} -0.24 & -0.4 \end{pmatrix}=\begin{pmatrix} 0.0576&0.096\\0.096 &0.16\end{pmatrix}

求得它的特征值是:0和0.7956,它的分布是一条直线,即图3所示的样子。

实际上不论在哪个地方进行线性化,分布都是直线。

        对于E^{'},它的J^{T}J就是对应E_{1}^{'}E_{2}^{'}J^{T}J的和,即\begin{pmatrix} 0.09&0.252\\0.252 &0.7056\end{pmatrix}+\begin{pmatrix} 0.0576&0.096\\0.096 &0.16\end{pmatrix}=\begin{pmatrix} 0.1476&0.348\\0.348 &0.8656\end{pmatrix},它的特征值:

\begin{vmatrix} 0.1476-\lambda & 0.348\\0.348 &0.8656- \lambda \end{vmatrix}=(0.1476-\lambda)(0.8656-\lambda)-0.348*0.348=\lambda^{2}-1.0132\lambda-0.07664262912

太复杂了,不求了。

总之,它有两个不为0的特征值,分布是一个椭圆,即图4所示的样子

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

神气爱哥

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值