向量的导数

向量的导数

向量的导数取决于你对向量的上下文和所涉及的变量维度。常见情况下,我们主要讨论以下两种情况:

1. 标量函数对向量的导数

如果一个标量函数 f ( x ) f(\mathbf{x}) f(x) 是关于向量 x = [ x 1 , x 2 , … , x n ] ⊤ \mathbf{x} = [x_1, x_2, \dots, x_n]^\top x=[x1,x2,,xn] 的函数,其导数通常是梯度 ∇ f ( x ) \nabla f(\mathbf{x}) f(x),表示为一个向量:

∇ f ( x ) = [ ∂ f ∂ x 1 ∂ f ∂ x 2 ⋮ ∂ f ∂ x n ] . \nabla f(\mathbf{x}) = \begin{bmatrix} \frac{\partial f}{\partial x_1} \\ \frac{\partial f}{\partial x_2} \\ \vdots \\ \frac{\partial f}{\partial x_n} \end{bmatrix}. f(x)= x1fx2fxnf .

例如:若 f ( x ) = x 1 2 + x 2 2 f(\mathbf{x}) = x_1^2 + x_2^2 f(x)=x12+x22,则:

∇ f ( x ) = [ ∂ f ∂ x 1 ∂ f ∂ x 2 ] = [ 2 x 1 2 x 2 ] . \nabla f(\mathbf{x}) = \begin{bmatrix} \frac{\partial f}{\partial x_1} \\ \frac{\partial f}{\partial x_2} \end{bmatrix} = \begin{bmatrix} 2x_1 \\ 2x_2 \end{bmatrix}. f(x)=[x1fx2f]=[2x12x2].

2. 向量函数对向量的导数

如果一个向量函数 f ( x ) \mathbf{f}(\mathbf{x}) f(x) 是关于向量 x = [ x 1 , x 2 , … , x n ] ⊤ \mathbf{x} = [x_1, x_2, \dots, x_n]^\top x=[x1,x2,,xn] 的函数,导数是一个 雅可比矩阵(Jacobian matrix)。设:

f ( x ) = [ f 1 ( x ) f 2 ( x ) ⋮ f m ( x ) ] , \mathbf{f}(\mathbf{x}) = \begin{bmatrix} f_1(\mathbf{x}) \\ f_2(\mathbf{x}) \\ \vdots \\ f_m(\mathbf{x}) \end{bmatrix}, f(x)= f1(x)f2(x)fm(x) ,

则雅可比矩阵是一个 m × n m \times n m×n 的矩阵:

J = ∂ f ∂ x = [ ∂ f 1 ∂ x 1 ∂ f 1 ∂ x 2 … ∂ f 1 ∂ x n ∂ f 2 ∂ x 1 ∂ f 2 ∂ x 2 … ∂ f 2 ∂ x n ⋮ ⋮ ⋱ ⋮ ∂ f m ∂ x 1 ∂ f m ∂ x 2 … ∂ f m ∂ x n ] . J = \frac{\partial \mathbf{f}}{\partial \mathbf{x}} = \begin{bmatrix} \frac{\partial f_1}{\partial x_1} & \frac{\partial f_1}{\partial x_2} & \dots & \frac{\partial f_1}{\partial x_n} \\ \frac{\partial f_2}{\partial x_1} & \frac{\partial f_2}{\partial x_2} & \dots & \frac{\partial f_2}{\partial x_n} \\ \vdots & \vdots & \ddots & \vdots \\ \frac{\partial f_m}{\partial x_1} & \frac{\partial f_m}{\partial x_2} & \dots & \frac{\partial f_m}{\partial x_n} \end{bmatrix}. J=xf= x1f1x1f2x1fmx2f1x2f2x2fmxnf1xnf2xnfm .

例如:若 f ( x ) = [ x 1 2 x 1 x 2 e x 2 ] \mathbf{f}(\mathbf{x}) = \begin{bmatrix} x_1^2 \\ x_1x_2 \\ e^{x_2} \end{bmatrix} f(x)= x12x1x2ex2 ,则:

J = [ ∂ f 1 ∂ x 1 ∂ f 1 ∂ x 2 ∂ f 2 ∂ x 1 ∂ f 2 ∂ x 2 ∂ f 3 ∂ x 1 ∂ f 3 ∂ x 2 ] = [ 2 x 1 0 x 2 x 1 0 e x 2 ] . J = \begin{bmatrix} \frac{\partial f_1}{\partial x_1} & \frac{\partial f_1}{\partial x_2} \\ \frac{\partial f_2}{\partial x_1} & \frac{\partial f_2}{\partial x_2} \\ \frac{\partial f_3}{\partial x_1} & \frac{\partial f_3}{\partial x_2} \end{bmatrix} = \begin{bmatrix} 2x_1 & 0 \\ x_2 & x_1 \\ 0 & e^{x_2} \end{bmatrix}. J= x1f1x1f2x1f3x2f1x2f2x2f3 = 2x1x200x1ex2 .

3. 向量函数对标量的导数

如果 x ( t ) \mathbf{x}(t) x(t) 是标量 t t t 的函数(即 x \mathbf{x} x 是一个随时间变化的向量),其导数是一个向量:

d x d t = [ d x 1 d t d x 2 d t ⋮ d x n d t ] . \frac{d\mathbf{x}}{dt} = \begin{bmatrix} \frac{dx_1}{dt} \\ \frac{dx_2}{dt} \\ \vdots \\ \frac{dx_n}{dt} \end{bmatrix}. dtdx= dtdx1dtdx2dtdxn .

总结

  • 标量对向量的导数 → 梯度(向量)。
  • 向量对向量的导数 → 雅可比矩阵(矩阵)。
  • 向量对标量的导数 → 向量。
    具体形式取决于你分析的场景和维度关系!
### 计算离散矩阵表示的函数θ与f在方形区域边界上法向量导数的内积 为了计算离散矩阵表示的函数 θ 与 f 在 [0,1]×[0,1] 方形区域边界上的法向量导数的内积,可以采用数值微分的方法来近似偏导数,并利用这些导数来构建所需的表达式。 #### 数值梯度计算 对于给定的离散化后的二维网格中的函数 \( \theta(x,y) \) 和 \( f(x,y) \),首先需要通过 `gradient` 函数获得其沿 x 轴和 y 轴方向的一阶差商作为对偏导数的估计: ```matlab % 假设 theta 和 f 是大小相同的 n-by-n 矩阵 [nx, ny] = size(theta); dx = (1 / nx); % 定义空间步长 dy = dx; % 对于正方形域,假设 dy=dx % 使用 gradient() 来获取一阶导数 [dthetadx, dthetady] = gradient(theta, dx, dy); [dfdx, dfdy] = gradient(f, dx, dy); ``` #### 边界条件处理 接着,在四个边上分别考虑外法线的方向并据此调整相应的导数项。注意这里只关心位于矩形四条边处的数据点: - 上下两端 (\(x\) 不变): 法线垂直向上或向下; 因此可以根据位置选取合适的导数组合形成最终的结果。 #### 内积运算 最后一步就是沿着每一边执行乘积累加操作以完成整个积分过程: ```matlab integral_value = 0; for i = 1:nx integral_value += (-dthetady(i,1)*dfdx(i,1)); % Bottom edge end for j = 2:ny integral_value += (+dthetadx(nx,j)*dfdy(nx,j)); % Right edge end for i = nx:-1:1 integral_value += (+dthetady(i,ny)*dfdx(i,ny)); % Top edge end for j = ny-1:-1:1 integral_value += (-dthetadx(1,j)*dfdy(1,j)); % Left edge end disp(['The computed inner product is ', num2str(integral_value)]); ``` 上述代码实现了基于离散数据集的数值逼近方案用于求取所需物理量[^1]。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

xwhking

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

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

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

打赏作者

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

抵扣说明:

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

余额充值