雅可比(Jacobi)矩阵的若干应用
2009-01-07 20:17阅读:
雅可比矩阵
以
m个
n元函数
u=
ui(
x1,
x2, ...,
xn)的偏导数

(其中,i=1,..,m;
j=1,...,n)为元素的矩阵:
该方程组是一个关于微分的线性方程组,中其系数矩阵即是上面的雅可比矩阵, 记为
J,
因此可写成矩阵形式:
应用一:机器人学
下图是
一个两自由度的平面机械手

容易求得未端点的位置为:

对其微分, 可得:

写成矩阵形式:

简写为
dx=
Jdθ,
其中的J称之为机械手的雅可比矩阵,反映了关节微小位移dθ与手部(手爪)微小运动dx之间的关系.
它可看成是从关节空间到操作空间运动速度的传动比,同时也可用来表示两空间之间力的传递关系.
假设关节速度为

,未端点速度为

,对
dx=
Jdθ两端同除以dt,
得

因此,
机械手的雅可比矩阵定义为它的操作空间速度与关节空间速度的线性变换,
v为机械手未端在操作空间的广义速度,
或称操作速度,
w为关节速度.
以上可以扩展到三维空间多自由度的情形.
应用二,非线性最小二乘
如果定义
ri(x)=yi-f(ti,x)
(其中i=1,...,m)为残留函数。
则问题为函数:
F(x)=r(x)
r(x)
的最小化问题。
这里的
ri(x),即是前面的ui。只不过(ti,
yi)作为不同的采样点,产生不同的ui,但它们形式是一样。采样点个数即ui的个数。
F(x)在xk
处的泰勒展开为:
这里的
为梯度向量,Ak为F(x)的Hessian矩阵(二阶导矩阵)。
根据牛顿法有:
由于

,其中

为
ri(x)的Hessian矩阵。
因此,步长
可通过解关于r(xk)、r(xk)的Jaccobi矩阵和Hessian矩阵的线性方程求解:
如果把上式中
的
一项去掉,就得到高斯-牛顿法;如果用uk