雅可比矩阵

雅可比矩阵[编辑]

维基百科,自由的百科全书

向量分析中,雅可比矩阵一阶偏导数以一定方式排列成的矩阵,其行列式称为雅可比行列式

还有,在代数几何中,代数曲线雅可比量表示雅可比簇:伴随该曲线的一个代数群,曲线可以嵌入其中。

它们全部都以数学家卡尔·雅可比命名;英文雅可比量"Jacobian"可以发音为[ja ˈko bi ən]或者[ʤə ˈko bi ən]

雅可比矩阵[编辑]

雅可比矩阵的重要性在于它体现了一个可微方程与给出点的最优线性逼近。因此,雅可比矩阵类似于多元函数的导数。 假设F:Rn→Rm 是一个从n维欧氏空间转换到m维欧氏空间的函数。这个函数由m个实函数组成: y1(x1, ..., xn), ..., ym(x1, ...,xn)。这些函数的偏导数(如果存在)可以组成一个m行n列的矩阵(m by 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)表示的

如果p是Rn中的一点,Fp点可微分,那么在这一点的导数由JF(p)给出(这是求该点导数最简便的方法)。在此情况下,由JF(p)描述的线性算子即接近点pF的最优线性逼近,x逼近于p

F(\mathbf{x}) \approx F(\mathbf{p}) + J_F(\mathbf{p})\cdot (\mathbf{x}-\mathbf{p})

例子[编辑]

球坐标系到直角坐标系的转化由F函数给出:R × [0,π] × [0,2π] → R3

 x_1 = r \sin\theta \cos\phi \,
 x_2 = r \sin\theta \sin\phi \,
 x_3 = r \cos\theta \,

此坐标变换的雅可比矩阵是

J_F(r,\theta,\phi) =\begin{bmatrix}\frac{\partial x_1}{\partial r} & \frac{\partial x_1}{\partial \theta} & \frac{\partial x_1}{\partial \phi} \\[3pt]\frac{\partial x_2}{\partial r} & \frac{\partial x_2}{\partial \theta} & \frac{\partial x_2}{\partial \phi} \\[3pt]\frac{\partial x_3}{\partial r} & \frac{\partial x_3}{\partial \theta} & \frac{\partial x_3}{\partial \phi} \\\end{bmatrix}=\begin{bmatrix} \sin\theta \cos\phi & r \cos\theta \cos\phi  & -r \sin\theta \sin\phi \\\sin\theta \sin\phi &  r \cos\theta \sin\phi & r \sin\theta \cos\phi \\ \cos\theta            & -r \sin\theta            &  0                    \end{bmatrix}.

R4的f函数:

 y_1 = x_1 \,
 y_2 = 5x_3 \,
 y_3 = 4x_2^2 - 2x_3 \,
 y_4 = x_3 \sin(x_1) \,

其雅可比矩阵为:

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}.

此例子说明雅可比矩阵不一定为方矩阵。

在动力系统中[编辑]

考虑形为x' = F(x)的动力系统F : Rn → Rn。如果F(x0) = 0,那么x0是一个驻点。系统接近驻点时的表现通常可以从JF(x0)的特征值来决定。

雅可比行列式[编辑]

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

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

例子[编辑]

设有函数F : R3 → R3,其分量为:

 y_1 = 5x_2 \,
 y_2 = 4x_1^2 - 2 \sin (x_2x_3) \,
 y_3 = x_2 x_3 \,

则它的雅可比行列式为:

\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.

从中我们可以看到,当x1x2同号时,F的取向相反;该函数处处具有反函数,除了在x1 = 0和x2 = 0时以外。

  • 1
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
雅可比矩阵是一个矩阵,其中每个元素是一个函数的偏导数。在神经网络推理中,输出通常是一个向量,例如类别概率。当我们想要了解每个输入变量对输出的影响时,可以使用雅可比矩阵。 在Python中,我们可以使用不同的方法来计算雅可比矩阵。一个常用的方法是使用自动微分库,例如NumPy或PyTorch。这些库提供了方便的函数来计算函数的导数。 例如,在PyTorch中,我们可以使用backward函数计算雅可比矩阵。可以通过调用该函数并传递一个单位矩阵作为参数来计算所有输入变量的导数。这将返回一个大小为(输出维度,输入维度)的雅可比矩阵。 另一个方法是使用符号计算库,例如SymPy。这允许我们在符号级别上定义函数,并使用diff函数来计算每个变量的导数。这样可以得到一个表示雅可比矩阵的符号表达式。 总之,Python提供了多种方法来计算雅可比矩阵,具体取决于您使用的库和所需的功能。您可以根据自己的需求选择合适的方法。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [在Python中计算神经网络的雅可比矩阵](https://blog.csdn.net/weixin_39971138/article/details/110837157)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值