雅可比矩阵

Jacobian矩阵和Hessian矩阵

发表于 2012 年 8 月 8 日

1. Jacobian

在向量分析中, 雅可比矩阵是一阶偏导数以一定方式排列成的矩阵, 其行列式称为雅可比行列式. 还有, 在代数几何中, 代数曲线的雅可比量表示雅可比簇:伴随该曲线的一个代数群, 曲线可以嵌入其中. 它们全部都以数学家卡尔·雅可比(Carl Jacob, 1804年10月4日-1851年2月18日)命名;英文雅可比量”Jacobian”可以发音为[ja ?ko bi ?n]或者[?? ?ko bi ?n].

雅可比矩阵

雅可比矩阵的重要性在于它体现了一个可微方程与给出点的最优线性逼近. 因此, 雅可比矩阵类似于多元函数的导数.

假设\(F\): \({R_n} \to {R_m}\)是一个从欧式n维空间转换到欧式m维空间的函数. 这个函数由m个实函数组成: y1(x1,…,xn), …, ym(x1,…,xn). 这些函数的偏导数(如果存在)可以组成一个m行n列的矩阵, 这就是所谓的雅可比矩阵:

\begin{bmatrix} \frac{\partial y_1}{\partial x_1} & \cdots & \frac{\partial y_1}{\partial x_n} \\ \vdots & \ddots & \vdots \\ \frac{\partial y_m}{\partial x_1} & \cdots & \frac{\partial y_m}{\partial x_n}  \end{bmatrix}

此矩阵表示为: \({J_F}({x_1}, \ldots ,{x_n})\), 或者\(\frac{{\partial ({y_1}, \ldots ,{y_m})}}{{\partial ({x_1}, \ldots ,{x_n})}}\).

这个矩阵的第i行是由梯度函数的转置yi(i=1,…,m)表示的.

如果\({\bf{p}}\)是\({R_n}\)中的一点, \(F\)在\({\bf{p}}\)点可微分, 那么在这一点的导数由\({J_F}({\bf{p}})\)给出(这是求该点导数最简便的方法). 在此情况下, 由\(F({\bf{p}})\)描述的线性算子即接近点\({\bf{p}}\)的\(F\)的最优线性逼近, \({\bf{x}}\)逼近于\({\bf{p}}\):

\(F({\bf{x}}) \approx F({\bf{p}}) + {J_F}({\bf{p}}) \cdot ({\bf{x}} – {\bf{p}})\)

雅可比行列式

如果m = n, 那么\(F\)是从n维空间到n维空间的函数, 且它的雅可比矩阵是一个方块矩阵. 于是我们可以取它的行列式, 称为雅可比行列式.

在某个给定点的雅可比行列式提供了 在接近该点时的表现的重要信息. 例如, 如果连续可微函数\(F\)在\({\bf{p}}\)点的雅可比行列式不是零, 那么它在该点附近具有反函数. 这称为反函数定理. 更进一步, 如果\({\bf{p}}\)点的雅可比行列式是正数, 则\(F\)在\({\bf{p}}\)点的取向不变;如果是负数, 则\(F\)的取向相反. 而从雅可比行列式的绝对值, 就可以知道函数\(F\)在\({\bf{p}}\)点的缩放因子;这就是为什么它出现在换元积分法中.

对于取向问题可以这么理解, 例如一个物体在平面上匀速运动, 如果施加一个正方向的力\(F\), 即取向相同, 则加速运动, 类比于速度的导数加速度为正;如果施加一个反方向的力\(F\), 即取向相反, 则减速运动, 类比于速度的导数加速度为负.

时间: 2024-12-07 09:20:58

雅可比矩阵的相关文章

vins-mono:雅可比矩阵的推导

imu残差公式 我们计算雅可比采用加扰动形式 以下是位置增量偏差对旋转qbk的雅可比推导 对应的雅可比矩阵的代码 jacobian_pose_i.block<3, 3>(O_P, O_R) = Utility::skewSymmetric(Qi.inverse() * (0.5 * G * sum_dt * sum_dt + Pj - Pi - Vi * sum_dt)); 原文地址:https://www.cnblogs.com/feifanrensheng/p/10723536.html

双视图下Bundle Adjustment重投影误差的雅可比矩阵

(一)建立问题模型 以第一帧的相机坐标系为世界坐标系,世界坐标系中存在n个空间点P,其在第二帧的投影点为p,现在希望求出相机的位姿R,t a a 原文地址:https://www.cnblogs.com/npulzb/p/11643498.html

深入理解图优化与g2o:图优化篇

前言 本节我们将深入介绍视觉slam中的主流优化方法——图优化(graph-based optimization).下一节中,介绍一下非常流行的图优化库:g2o. 关于g2o,我13年写过一个文档,然而随着自己理解的加深,越发感觉不满意.本着对读者更负责任的精神,本文给大家重新讲一遍图优化和g2o.除了这篇文档,读者还可以找到一篇关于图优化的博客: http://blog.csdn.net/heyijia0327 那篇文章有作者介绍的一个简单案例,而本文则更注重对图优化和g2o的理解与评注. 本

Machine Learning 学习笔记 (2) —— 另一种寻找极值的方法:牛顿法

[总目录]http://www.cnblogs.com/tbcaaa8/p/4415055.html 1. 用牛顿法解方程 牛顿法是一种求解方程的迭代算法,也可以用于方程组的求解.其思想是利用方程(尤其是非线性方程)的线性部分,对原方程进行近似.不失一般性,考虑方程f(x)=0.对f(x)在x=t处进行泰勒展开,可得f(x)=f(t)+f'(t)(x-t)+... 取线性部分代替f(x),带入方程f(x)=0,可得f(t)+f'(t)(x-t)=0 ,进而解出x=t-f(t)/f'(t).将方程

g2o求解BA 第10章

1.g2o_bal_class.h1.1 projection.hg2o还是用图模型和边,顶点就是相机和路标,边就是观测,就是像素坐标.只不过这里的相机是由旋转(3个参数,轴角形式,就是theta*nx,theta*ny,theta*ny),位移(3个参数),f,k1,k2.就是之前BA模型的实现.但是这里归一化平面坐标取得是负值,而且最后没有加cx,cy.具体实现在projection.h中的CamProjectionWithDistortion函数实现的,而其中世界坐标转成相机坐标是由rot

相机标定:关于用Levenberg-Marquardt算法在相机标定中应用

LM算法在相机标定的应用共有三处. (1)单目标定或双目标定中,在内参固定的情况下,计算最佳外参.OpenCV中对应的函数为findExtrinsicCameraParams2. (2)单目标定中,在内外参都不固定的情况下,计算最佳内外参.OpenCV中对应的函数为calibrateCamera2. (3)双目标定中,在左右相机的内外参及左右相机的位姿都不固定的情况下,计算最佳的左右相机的内外参及最佳的左右相机的位姿矩阵.OpenCV中对应的函数为stereoCalibrate. 本文文阅读前提

直接法

前言 直接法是视觉里程计另一主要分支,它与特征点法有很大不同.随着SVO.LSD-SLAM等直接法SLAM方案的流行,直接法本身也得到越来越多的关注.特征点法与直接究竟谁更好一些,是近年视觉里程计研究领域一个非常有趣的问题.本讲,我们将介绍直接法的原理,并利用g2o实现基于直接法的视觉里程计. 1. 直接法的引出 尽管特征点法在视觉里程计中占据主流地位,研究者们认识它至少有以下几个缺点: 关键点的提取与描述子的计算非常耗时.实践当中,SIFT目前在CPU上是无法实时计算的,而ORB也需要近20毫

Matlab:非线性热传导(抛物方程)问题

函数文件1:real_fun.m 1 function f=real_fun(x0,t0) 2 %精确解 3 f=4*x0*(1-x0)*sin(t0); 函数文件2:F.m 1 function f=F(N,u,U,t,h1,h2) 2 %非线性方程组 3 %h1是x的步长,h2是t的步长 4 %u表示迭代节点,上一时刻的数值解 5 %h表示时间节点上的步长 6 %N表示空间节点的步数 7 a0=0.5*t^4*h2*N^2; 8 f(1,1)=a0*(U(2)^2-2*U(1)^2)+h2*

视觉SLAM中的数学基础 第三篇 李群与李代数

视觉SLAM中的数学基础 第三篇 李群与李代数 前言 在SLAM中,除了表达3D旋转与位移之外,我们还要对它们进行估计,因为SLAM整个过程就是在不断地估计机器人的位姿与地图.为了做这件事,需要对变换矩阵进行插值.求导.迭代等操作.例如,在经典ICP问题中,给定了两组3D点,我们要计算它们之间的变换矩阵.假设第一组的3D点为$\mathbf{P}=\{ \mathbf{p}_i | i = [1,2, \ldots, N] \}$,第二组3D点为$\mathbf{Q}=\{ \mathbf{q}