在向量分析中,雅可比矩阵(也称作Jacobi矩阵,英语:Jacobian matrix)是函数的一阶偏导数以一定方式排列成的矩阵。
当其为方形矩阵时,其行列式称为雅可比行列式(Jacobi determinant)。要注意的是,在英文中雅可比矩阵跟雅可比行列式都可称作Jacobian。[1]
其重要性在于,如果函数 f : ℝn → ℝm 在点 x 可微的话,在点 x 的雅可比矩阵即为该函数在该点的最佳线性逼近,也代表雅可比矩阵是单变数实数函数的微分在向量值多变数函数的推广,在这种情况下,雅可比矩阵也被称作函数 f 在点 x 的微分或者导数。
在代数几何中,代数曲线的雅可比行列式'表示雅可比簇:伴随该曲线的一个代数群,曲线可以嵌入其中。
它们全部都以普鲁士数学家卡尔·雅可比命名。
雅可比矩阵
假设某函数从 f : ℝn → ℝm, 从 x ∈ ℝn 映射到 向量 f(x) ∈ ℝm, 其雅可比矩阵是一 m×n 的矩阵,换句话讲也就是从 ℝn 到 ℝm 的线性映射,其重要意义在于它表现了一个多变数向量函数的最佳线性逼近。因此,雅可比矩阵类似于单变数函数的导数。
此函数 f 的雅可比矩阵 J 为 m×n 的矩阵,一般由以下方式定义:
![{\displaystyle \mathbf {J} ={\begin{bmatrix}{\dfrac {\partial \mathbf {f} }{\partial x_{1}}}&\cdots &{\dfrac {\partial \mathbf {f} }{\partial x_{n}}}\end{bmatrix}}={\begin{bmatrix}{\dfrac {\partial f_{1}}{\partial x_{1}}}&\cdots &{\dfrac {\partial f_{1}}{\partial x_{n}}}\\\vdots &\ddots &\vdots \\{\dfrac {\partial f_{m}}{\partial x_{1}}}&\cdots &{\dfrac {\partial f_{m}}{\partial x_{n}}}\end{bmatrix}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/74e93aa903c2695e45770030453eb77224104ee4)
矩阵的分量可表示成:
![{\displaystyle \mathbf {J} _{ij}={\frac {\partial f_{i}}{\partial x_{j}}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/2e44836c4c7fef55a3ef28bbc0fddfa3770b65dc)
雅可比矩阵的其他常用符号还有:
、
、
或者 ![{\displaystyle {\frac {\partial (f_{1},\ldots ,f_{m})}{\partial (x_{1},\ldots ,x_{n})}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/d4549450a0cc1be5355876ff3911861fb3b1a4ba)
此矩阵的第
行是由函数
的梯度函数所表示的,
。
如果
是
中的一点,
在
点可微分,根据数学分析,
是在这点的导数。在此情况下,
这个线性映射即
在点
附近的最优线性逼近,也就是说当
足够靠近点
时,我们有
![{\displaystyle f(x)\approx f(p)+\mathbf {J} _{\mathbf {f} }(p)\cdot (x-p)}](//wikimedia.org/api/rest_v1/media/math/render/svg/c72fa49ba51fb457d6530fcc53f523ac470c2802)
讲更详细点也就是:
![{\displaystyle \mathbf {f} (\mathbf {x} )=\mathbf {f} (\mathbf {p} )+\mathbf {J} _{\mathbf {f} }(\mathbf {p} )(\mathbf {x} -\mathbf {p} )+o(\|\mathbf {x} -\mathbf {p} \|)}](//wikimedia.org/api/rest_v1/media/math/render/svg/a267bfe872cecd8691135b1cf44d22e2d0f2cd30)
其中,o 代表小o符号,‖x − p‖ 为 x 与 p 之间的距离。
例子
例一
由球坐标系到直角坐标系的转化由 F: ℝ+ × [0, π) × [0, 2π) → ℝ3 函数给出,其分量为:
![{\displaystyle {\begin{aligned}x&=r\sin \theta \cos \varphi ;\\y&=r\sin \theta \sin \varphi ;\\z&=r\cos \theta \end{aligned}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/251af2f9e253d45778b79ea234e261246f0e57b0)
此坐标变换的雅可比矩阵是
![{\displaystyle \mathbf {J} _{\mathbf {F} }(r,\theta ,\varphi )={\begin{bmatrix}{\dfrac {\partial x}{\partial r}}&{\dfrac {\partial x}{\partial \theta }}&{\dfrac {\partial x}{\partial \varphi }}\\[1em]{\dfrac {\partial y}{\partial r}}&{\dfrac {\partial y}{\partial \theta }}&{\dfrac {\partial y}{\partial \varphi }}\\[1em]{\dfrac {\partial z}{\partial r}}&{\dfrac {\partial z}{\partial \theta }}&{\dfrac {\partial z}{\partial \varphi }}\end{bmatrix}}={\begin{bmatrix}\sin \theta \cos \varphi &r\cos \theta \cos \varphi &-r\sin \theta \sin \varphi \\\sin \theta \sin \varphi &r\cos \theta \sin \varphi &r\sin \theta \cos \varphi \\\cos \theta &-r\sin \theta &0\end{bmatrix}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/db70522712b56ee256247603ef074d20313e5ba1)
其雅可比行列式为 r2 sin θ。以体积元变换为例,由于 dV = dx dy dz,如果做变数变换,则其体积元(Volume element,dV),会变成:dV = r2 sin θ dr dθ dφ。
例二
F : ℝ3 → ℝ4,其各分量为
![{\displaystyle y_{1}=x_{1}\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/d5dd9c8e322ef5d8d0b9779e2ab5c20f7bd9cf4d)
![{\displaystyle y_{2}=5x_{3}\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/9d498e10016f6b52729a2ea3cee4dadc70a839a2)
![{\displaystyle y_{3}=4x_{2}^{2}-2x_{3}\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/b6301612fb299c9efc6227e8cc1fe3d4952cb845)
![{\displaystyle y_{4}=x_{3}\sin x_{1}\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/5cb3f4e3af6e6d0e0196fc81be2bebdac56828e4)
其雅可比矩阵为:
![{\displaystyle J_{F}(x_{1},x_{2},x_{3})={\begin{bmatrix}{\frac {\partial y_{1}}{\partial x_{1}}}&{\frac {\partial y_{1}}{\partial x_{2}}}&{\frac {\partial y_{1}}{\partial x_{3}}}\\[3pt]{\frac {\partial y_{2}}{\partial x_{1}}}&{\frac {\partial y_{2}}{\partial x_{2}}}&{\frac {\partial y_{2}}{\partial x_{3}}}\\[3pt]{\frac {\partial y_{3}}{\partial x_{1}}}&{\frac {\partial y_{3}}{\partial x_{2}}}&{\frac {\partial y_{3}}{\partial x_{3}}}\\[3pt]{\frac {\partial y_{4}}{\partial x_{1}}}&{\frac {\partial y_{4}}{\partial x_{2}}}&{\frac {\partial y_{4}}{\partial x_{3}}}\\\end{bmatrix}}={\begin{bmatrix}1&0&0\\0&0&5\\0&8x_{2}&-2\\x_{3}\cos x_{1}&0&\sin x_{1}\end{bmatrix}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/ebf91f57907e1656c70decff4f02a0acf8fbdacf)
此例子说明雅可比矩阵不一定为方阵。
在动力系统中
考虑形为
的动力系统,
。如果
,那么
是一个临界点。系统接近临界点时的行为跟
的特征值相关。
雅可比行列式
如果 m = n,那么 F 是从 ℝn 映射到 ℝn 的函数,且它的雅可比矩阵是一个方阵。于是我们可以取它的行列式,称为雅可比行列式。
在某个给定点的雅可比行列式提供了 F 在接近该点时的表现的重要资讯。例如,如果连续可微函数 F 在 p 点的Jacobi行列式不等于零,那么它在该点附近有 F 的反函数。这称为反函数定理。更进一步,如果 p 点的Jacobi行列式是正数,则 F 在 p 点保持定向(preserves orientation);如果是负数,则 F 逆转定向(reverses orientation)。而从Jacobi行列式的绝对值,就可以知道函数 F 在 p 点附近是放大或缩小体积;这就是它出现在换元积分法中的原因。
例子一
设有函数 F : ℝ3 → ℝ3,其分量为:
![{\displaystyle y_{1}=5x_{2}\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/64c8c296bd6d1f2f70cd7b3eca245787094e8ae4)
![{\displaystyle y_{2}=4x_{1}^{2}-2\sin(x_{2}x_{3})\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/85b38bf8221947f50508fda77d37f1b7fb015d14)
![{\displaystyle y_{3}=x_{2}x_{3}\,}](//wikimedia.org/api/rest_v1/media/math/render/svg/85b37c5c6c9cacc0857ca25e1778113631b97c00)
则它的Jacobi行列式为:
![{\displaystyle {\begin{vmatrix}0&5&0\\8x_{1}&-2x_{3}\cos(x_{2}x_{3})&-2x_{2}\cos(x_{2}x_{3})\\0&x_{3}&x_{2}\end{vmatrix}}=-8x_{1}\cdot {\begin{vmatrix}5&0\\x_{3}&x_{2}\end{vmatrix}}=-40x_{1}x_{2}}](//wikimedia.org/api/rest_v1/media/math/render/svg/d853ffdff65f9e7e189d418aca5450b9beb0fa14)
从中我们可以看到,当 x1 和 x2 同号时,F 逆转定向;该函数处处具有反函数,除了在 x1 = 0 或 x2 = 0 的点。
例子二
这是一个与巴塞尔问题
较为相似的级数
的求解方法,首先可以转化为二重积分(在这里 D1 指 x 与 y 皆为从 0 到 1 的正方形区域):
![{\displaystyle \sum _{n=0}^{\infty }{\frac {1}{(2n+1)^{2}}}=\iint \limits _{D_{1}}\sum _{n=1}^{\infty }(xy)^{2n}\mathrm {d} x\mathrm {d} y=\iint \limits _{D_{1}}{\frac {\mathrm {d} x\mathrm {d} y}{1-x^{2}y^{2}}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/d8e81ad8bfc54ea67f977535883996559d944432)
此时定义映射 F : ℝ2 → ℝ2,满足:
![{\displaystyle {\begin{cases}u=\arctan \left(x{\sqrt {\dfrac {1-y^{2}}{1-x^{2}}}}\right)\\v=\arctan \left(y{\sqrt {\dfrac {1-x^{2}}{1-y^{2}}}}\right)\end{cases}}\iff {\begin{cases}x={\dfrac {\sin u}{\cos v}}\\y={\dfrac {\sin v}{\cos u}}\end{cases}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/74648b359311a31d66916b87dd484967701115ae)
于是有相应的雅可比行列式:
![{\displaystyle {\begin{vmatrix}{\dfrac {\partial x}{\partial u}}&{\dfrac {\partial x}{\partial v}}\\{\dfrac {\partial y}{\partial u}}&{\dfrac {\partial y}{\partial v}}\end{vmatrix}}={\begin{vmatrix}{\dfrac {\cos u}{\cos v}}&{\dfrac {\sin u\sin v}{\cos ^{2}v}}\\{\dfrac {\sin u\sin v}{\cos ^{2}u}}&{\dfrac {\cos v}{\cos u}}\end{vmatrix}}=1-{\frac {\sin ^{2}u\sin ^{2}v}{\cos ^{2}u\cos ^{2}v}}=1-x^{2}y^{2}}](//wikimedia.org/api/rest_v1/media/math/render/svg/d9a74e95ae1caf1ae7ce9ffcdd11ab50795946af)
因此
,并且将正方形 D1 映射成 u>0、v>0、u+v<π/2 的等腰直角三角形,记为 D2,得到:
![{\displaystyle \iint \limits _{D_{1}}{\frac {\mathrm {d} x\mathrm {d} y}{1-x^{2}y^{2}}}=\iint \limits _{D_{2}}\mathrm {d} u\mathrm {d} v=\int _{0}^{\frac {\pi }{2}}\left(\int _{0}^{{\frac {\pi }{2}}-v}\mathrm {d} u\right)\mathrm {d} v={\frac {\pi ^{2}}{8}}}](//wikimedia.org/api/rest_v1/media/math/render/svg/dfab3a91cfa61801e8750601131d3e546aaa6f59)
逆矩阵
根据反函数定理,一个可逆函数(存在反函数的函数)的雅可比矩阵的逆矩阵即为该函数的反函数的雅可比矩阵。即,若函数
在点
的雅可比矩阵是连续且可逆的,则
在点
的某一邻域内也是可逆的,且有
![{\displaystyle J_{F^{-1}}\circ f=J_{F}^{-1}}](//wikimedia.org/api/rest_v1/media/math/render/svg/9c06c7494056b97e0940ee96b6af7c168b7a7917)
成立。相反,倘若雅可比行列式在某一个点不为零,那么该函数在这个点的某一邻域内可逆(存在反函数)。
一个多项式函数的可逆性与未经证明的雅可比猜想有关。其断言,如果函数的雅可比行列式为一个非零实数(相当于其不存在复零点),则该函数可逆且其反函数也为一个多项式。
参看
参考资料
W., Weisstein, Eric. Jacobian. mathworld.wolfram.com. [2 May 2018]. (原始内容存档于3 November 2017).
外部链接