微分几何 - 曲率与曲线坐标系

背景: 坐标不是向量的分量!

笛卡尔坐标系太成功了，它成功到让我们产生了一个错误的等价关系：

“坐标 = 向量的分量”

其实这个等价只在笛卡尔坐标系和极少数全局线性坐标系里成立。在一般的坐标系（比如极坐标）里，坐标只是“门牌号”，不是“尺子上的刻度”。

我们来把这堵墙拆了，然后你会感到一种无比清爽的认知升级。

先破：为什么你会觉得 $(r, \theta)$ 应该是“$r$ 倍的什么 + $\theta$ 倍的什么”？

因为笛卡尔的 $(x, y)$ 就是这样：

• 有一个全局均匀的基 $\{\hat{\mathbf{i}}, \hat{\mathbf{j}}\}$。
• 任何位置向量 $\vec{p} = x\hat{\mathbf{i}} + y\hat{\mathbf{j}}$。
• 坐标 $(x, y)$ 就是这两个基向量的系数。

所以你的大脑被训练成：“如果给我一个坐标对 $(a, b)$，那就一定存在两个固定的基向量 $\mathbf{e}_1, \mathbf{e}_2$，使得我的向量是 $a\mathbf{e}_1 + b\mathbf{e}_2$。”

在极坐标里，这个想法完全错在这里：

• 确实存在一个 $r$，可以乘在一个基向量上。
• 但是 $\theta$ 根本不是一个系数，它是一个角度，量纲就不是长度，你没法说“$\theta$ 倍的某个长度基向量”。如果硬写“$\theta \mathbf{e}_2$”，那岂不是把角度和长度相加？这违反量纲原理，没有意义。

再立：那极坐标的 $(r, \theta)$ 到底是什么？它的“基”又是什么？

我们分两层理解：

第一层：坐标是“地址”，不是“分量”

想象你住在一个环形城市，每栋房子的地址是**（距离中心广场的距离，从正东起算的角度）**。比如你家住 $(500\text{m}, 30^\circ)$。

• 这个地址本身没有在描述一个向量，它只是在告诉你“怎么走过去”：先往东偏30°方向走500米。
• 你不会说：“我的位置向量是 500 × 某个向量 + 30 × 另一个向量”，因为30°不是一个长度。

所以极坐标的 $(r, \theta)$ 只是参数化平面上的点的一种方式，就像地图上的经纬度。

第二层：在每一个点上，我们安装“尺子”（局部基向量）

为了在某个点上做向量的计算（比如速度、加速度），我们需要在那个点处安装一组互相垂直的、长度为1的标架：

• $\hat{\mathbf{u}}_r$ ：从原点指向该点（径向向外）
• $\hat{\mathbf{u}}_\theta$：垂直于径向，逆时针方向

关键：这组基向量是“贴在点上的”，不同点的 $\hat{\mathbf{u}}_r, \hat{\mathbf{u}}_\theta$ 方向不同！

现在，从原点指向该点的位置向量，它恰好完全沿着该点的径向方向，所以它在这个该点的局部标架下，分量是 $(r, 0)$。

用公式写：

$$\vec{p} = r \cdot \hat{\mathbf{u}}_r(\theta) + 0 \cdot \hat{\mathbf{u}}_\theta(\theta)$$

你看，这里出现了 $r$ 倍的东西，但 $\theta$ 并没有作为系数出现。$\theta$ 的作用是决定了 $\hat{\mathbf{u}}_r$ 和 $\hat{\mathbf{u}}_\theta$ 的朝向。

那么“$r$ 倍的什么 + $\theta$ 倍的什么”这种想法到底错在哪？用一个比喻彻底揭穿

比喻：地球上的经纬度

北京的经纬度近似是 $(39^\circ N, 116^\circ E)$。

• 你能说“北京的位置向量 = 39° × (某个指向北的基) + 116° × (某个指向东的基)”吗？
• 不行！因为39°是角度，不是长度，你没法拿它去乘一个长度基向量。

那你怎么用经纬度描述北京的位置？

• 从地心出发，先找到一个依赖于这些角度的方向，然后沿这个方向走地球半径 $R$。
• 位置向量 = $R \cdot \hat{\mathbf{u}}_r(39^\circ N, 116^\circ E)$。
• 这和你极坐标里的 $r \cdot \hat{\mathbf{u}}_r(\theta)$ 完全一样。

所以，$(r, \theta)$ 和 $(纬度, 经度)$ 一样，是“曲线坐标”，它们的作用是确定基向量的朝向，而不是直接当系数。

一句让你永远不再迷糊的口诀

请默念三遍：

“坐标是门牌，基向量是门牌旁边的方向牌。位置向量就是门牌号中的距离，沿着方向牌指的方向走那么远。”

$(r, \theta)$ 里：

• $\theta$ 告诉你“方向牌”指向哪；
• $r$ 告诉你沿着方向牌走多远；
• 另一个方向牌（切向）在这次定位中没用到，所以它的系数是0。

你会发现，这样一想，笛卡尔坐标不过是“所有门牌旁边的方向牌都指向同一个方向”的特例而已。

这个认知墙一旦倒了，你就能真正理解为什么极坐标下速度是 $(\dot r, r\dot\theta)$，加速度是 $(\ddot r - r\dot\theta^2, r\ddot\theta + 2\dot r \dot\theta)$，以及为什么再往后学，连“位置向量”这个概念本身在弯曲空间（比如球面）上就失效了——因为球面上根本没有全局的位置向量，只有局部的“门牌”和“方向牌”。

Manifold上某点的切空间

1. 定义

一个球面或任何弯曲的空间（数学上叫流形（Manifold））。在它上面某一点 $p$，你可以做一个“切平面”，这个切平面就是点 $p$ 的切空间（Tangent Space），记作 $T_pM$。切空间里的元素叫切向量（Tangent Vector）。

2. 方向导数

在不同的数学领域，对“方向导数”这个名词有两种常见的约定：

定义 A（微积分课本常用）

方向导数 只对单位向量定义.
即规定 $\lVert\vec{v}\rVert=1$，此时 $\frac{\partial f}{\partial \vec{v}}(p)$ 称为函数在 $p$ 点沿方向 $\vec{v}$ 的方向导数。

按照这个定义，$(3,4)$ 不是单位向量（长度是5），就不能直接带入。
正确的做法是先归一化：方向是 $\vec{u}=(\frac35,\frac45)$，方向导数为

$$\frac{\partial f}{\partial \vec{u}}(p) = 2\cdot\frac35 + 4\cdot\frac45 = \frac{6}{5} + \frac{16}{5} = \frac{22}{5} = 4.4.$$

定义 B（微分几何、流形上常用）

方向导数 对任意切向量定义，不需要归一化。
即给定向量 $\vec{v}$，直接公式 $\frac{\partial f}{\partial \vec{v}}(p)=\nabla f(p)\cdot\vec{v}$，允许 $\vec{v}$ 有任意长度。

这时候，“方向导数”这个名称其实带有误导性。它更准确的名字应该是“沿着切向量的导数”或者“切向量作用在函数上的结果”。

3. 切向量的三种等价身份

在微分几何或多元微积分中，一个“切向量”可以等价地展现为：

1. 几何箭头（仅限欧氏空间或嵌入子流形）
   比如 $\vec{v}=(3,4)$，它就画在平面上，是一个有向线段。

2. 导子（derivation）
   在抽象流形上，切向量被 定义 为一个满足莱布尼兹律的导子，例如
   $v = 3\frac{\partial}{\partial x}\Big|_{p} + 4\frac{\partial}{\partial y}\Big|_{p}.$
   这个导子作用在函数 $f$ 上得到方向导数：$v(f) = 3\frac{\partial f}{\partial x} + 4\frac{\partial f}{\partial y}$。

3. 坐标数组
   一旦选定了基底（比如坐标基底 $\{\frac{\partial}{\partial x},\frac{\partial}{\partial y}\}$），每个导子就唯一对应一个坐标数组 $(3,4)$，并且这种对应是线性同构的。

所以，当说“切向量 $(3,4)$”时，我实际上是用坐标数组 代表 了导子 $3\partial_x + 4\partial_y$。这在数学物理中是标准简化.

4. 通过Metric度规赋予内积, 晋升为内积空间

黎曼度规（Riemannian metric）：

在光滑流形 $M$ 上，一个黎曼度规 $g$ 是在每一点 $p\in M$ 的切空间 $T_pM$ 上赋予一个内积 $g_p$，

具体来说，$g_p$ 是一个二元运算：
$g_p : T_pM \times T_pM \longrightarrow \mathbb{R}$
满足：

• 双线性（bilinear）：对两个输入都是线性的；
• 对称（symmetric）：$g_p(v,w) = g_p(w,v)$；
• 正定（positive-definite）：$g_p(v,v) \ge 0$，且等号当且仅当 $v=0$ 时成立。

这样，每一切空间 $(T_pM, g_p)$ 就变成了一个内积空间（inner product space）。

注：有时也会用 $\langle v, w \rangle_p$ 或 $\langle v, w\rangle$ 来表示 $g_p(v,w)$。

在局部坐标下如何计算内积？

设 $(U, \varphi = (x^1,\dots,x^n))$ 是 $p$ 附近的一个坐标卡，那么切空间 $T_pM$ 有一组坐标基底：
$\left\{ \frac{\partial}{\partial x^1}\Big|_p,\ \dots,\ \frac{\partial}{\partial x^n}\Big|_p \right\}.$

度规 $g_p$ 作用在这些基底上的值称为度规的分量：
$g_{ij}(p) := g_p\!\left( \frac{\partial}{\partial x^i},\ \frac{\partial}{\partial x^j} \right).$

• 由对称性，$g_{ij} = g_{ji}$；

• 由正定性，矩阵 $[g_{ij}]$ 在所有点都是正定的。对于任意两个切向量（用分量表示）
  $v = \sum_{i=1}^n v^i \frac{\partial}{\partial x^i}, \qquad w = \sum_{j=1}^n w^j \frac{\partial}{\partial x^j},$

• 它们的内积由双线性展开得到：
  $\langle v, w \rangle = \sum_{i=1}^n \sum_{j=1}^n g_{ij}(p)\, v^i\, w^j.$

（通常采用爱因斯坦求和约定：$\langle v, w \rangle = g_{ij} v^i w^j$，自动对重复指标求和。）

5. 弧长微元 ds 自然诱导出一个度规Metric和一个内积

上节只是讲到度规的通用形式, 下面看具体的一种实现:

内积的坐标计算公式为：

$$\langle v, v \rangle = \sum_{i=1}^n \sum_{j=1}^n g_{ij}\, v^i v^j.$$

展开这个双和就是：

$$\langle v, v \rangle = g_{11} (v^1)^2 + g_{12} v^1 v^2 + g_{21} v^2 v^1 + g_{22} (v^2)^2.$$

因为 $g_{12}=g_{21}$，所以通常写成

$$\langle v, v \rangle = g_{11} (v^1)^2 + 2g_{12} v^1 v^2 + g_{22} (v^2)^2.$$

现在只要知道 $g_{ij}$ 的具体数值，代入分量就算出来了。

R^2 平面直角坐标系

• 切向量 $(dx\frac{\partial}{\partial x},\ dy\frac{\partial}{\partial y})$
• 弧长微元 $(ds)^2=g_{11}(dx)^2 + 2g_{12}dxdy + g_{22}(dy)^2 = dx^2 + dy^2$

自动诱导了一个度规:

$[g_{ij}] = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} \quad \text{即} \quad g_{xx}=1,\ g_{yy}=1,\ g_{xy}=g_{yx}=0.$

平面极坐标系

• 切向量 $(dr\frac{\partial}{\partial r},\ d_\theta\frac{\partial}{\partial \theta})$
• 弧长微元 $(ds)^2=g_{11}(d_r)^2 + 2g_{12}d_r d_\theta + g_{22}(d_\theta)^2 = d_r^2 + r^2 d_\theta^2$

自动诱导了一个度规:

$[g_{ij}] = \begin{pmatrix} 1 & 0 \\ 0 & r^2 \end{pmatrix} \quad \text{即} \quad g_{rr} = 1,\ g_{\theta\theta} = r^2, \ g_{r\theta} = g_{\theta r} = 0.$

此时，我们就说：通过 $(ds)^2$ 的表达式，自然地给出了度规张量（内积矩阵）在极坐标下的分量。这就是“弧长平方自然诱导度规张量”的原始含义——它并不神秘，仅仅是换坐标后，长度公式变成了一个关于坐标微分的二次型，其系数就是内积矩阵。

6. 物理基底 (正交归一基底)

在自然基底下, 以平面极坐标为例, 某个切向量 $(v^r\partial{r}, v^\theta \partial{\theta})$, 等价表示为$(v^r, v^\theta)$ 或 $v^r \partial{r} + v^\theta \partial{\theta}$

从度规$[g_{ij}] = \begin{pmatrix} 1 & 0 \\ 0 & r^2 \end{pmatrix}$ 看出自然基虽然垂直, 但是$e_{\theta}$不是归一化, 可以进行归一化形成一个新基 $e_{\theta} = \frac{1}{r} \partial{\theta}$ , 配合 $e_r = \partial{r}$ 构成一组正交归一的新基底, 叫做物理基底

此时:切向量 $(v^r e_r, r v^\theta e_{\theta})$, 等价表示为$(v^r, r v^\theta)$ 或 $v^r e_r + r v^\theta e_{\theta}$

物理对应

如果 $v^r, v^\theta$ 分别赋予物理上的径向微元和角度微元量纲 $d_r, d_{\theta}$, 那么 $d_r \mathbf{e_r} + r d_{\theta} \mathbf{e_{\theta}}$ 就是表示 $\mathbf{ds}$

7. 自然基底对时间的导数

曲线坐标系通常不是全局线性不变的, 对于流行上某个位置 P, 基底自然不同, 因此基底也随着时间变化. 自然基底对时间的导数完全由度规张量唯一决定 (通过克氏符来描述), 也就是说, 物理基底同样完全由自然基底的度规张量唯一决定, 而根据链式法则, 自然(或物理)基底对时间的求导可以先对坐标求导然后由坐标再对时间求导来得到. 对于某个曲线坐标系, 通常有一个速查表格, 举例如下:

圆柱坐标中：

• 自然基：$\partial_r,\ \partial_\theta,\ \partial_z$
• 物理基：$\mathbf{e}_r=\partial_r,\ \mathbf{e}_\theta=\frac1r\partial_\theta,\ \mathbf{e}_z=\partial_z$

	$\partial/\partial r$	$\partial/\partial \theta$	$\partial/\partial z$
$\mathbf{e}_r$	0	$\mathbf{e}_\theta$	0
$\mathbf{e}_\theta$	0	$-\mathbf{e}_r$	0
$\mathbf{e}_z$	0	0	0

8. 位置,速度,和加速度向量

• 速度向量和加速度向量天然地生长在切空间里.
• 位置向量在平坦空间里，也可以被“请进”切空间，所以看起来和速度、加速度用了同一套基底；但在弯曲空间里，这个“请”的动作就再也做不到了。

  注意: 区分坐标系的弯曲性和空间的弯曲性. 如: 柱坐标, 极坐标等等可以是在平直空间里的弯曲坐标系.

• 通常我们研究 $\mathbb{R^3}$ 中的物理问题,位置向量可以被视作一个切向量

二维极坐标下的物理量

位置: $\mathbf{r} = r\mathbf{e}_r$
$(r, 0)$ 为在 P 点处的局部物理曲线坐标下的位置向量

速度: $\mathbf{\dot r} = \dot r \mathbf{e}_r + r \frac{d\mathbf{e}_r}{d\theta} \frac{d\theta}{dt} = \dot r \mathbf{e}_r + r\dot \theta \mathbf{e}_\theta$
$\dot r$ 为径向速度, $r \dot \theta$ 为切向速度

加速度: $\mathbf{\ddot r} = (\ddot{r} - r\dot{\theta}^2)\mathbf{e}_r + (r\ddot{\theta} + 2\dot{r}\dot{\theta})\mathbf{e}_\theta$
$\ddot r$ 为径向上的单纯由于径向速度大小变化产生的加速度
$-r \dot \theta^2$ 为径向上单纯由于角速度w导致的向心加速度, 方向指向中心, 非惯性系中的离心力（centrifugal force）
$r\ddot{\theta}$ 为角速度大小变化导致的切向加速度
$2\dot{r}\dot{\theta}$ 径向速度 + 角速度并存共同导致的横向加速度, 非惯性系中的科里奥利力（Coriolis force）

9. 嵌入在三维欧氏空间的一维流形

• 三维欧氏空间 $\mathbb{R}^3$。这是一个平直的三维流形，带有全局笛卡尔坐标 $(x,y,z)$。
• 一条嵌入在 $\mathbb{R}^3$ 中的光滑曲线 $C$。曲线本身就是一个一维流形。
• 曲线上的“门牌号”：
  选用弧长参数 $s$ 作为这一维流形的曲线坐标。
  在 $C$ 上的每个点 $p$，都贴着一个唯一的门牌 $s$。

对于曲线上的一个点 $p$，我们要区分两个完全不同的切空间：

切空间	维度	坐标基底	作用
曲线自己的切空间 $T_p C$	1 维	$\frac{\partial}{\partial s}$	只负责沿曲线方向的切向量
外围空间 $\mathbb{R}^3$ 的切空间 $T_p \mathbb{R}^3$	3 维	$\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z}$	包含所有从 $p$ 点出发的（三维）切向量

显然，$T_p C$ 是 $T_p \mathbb{R}^3$ 的一个一维线性子空间。
既然 $T_p \mathbb{R}^3$ 是三维的，可以在它里面挑三个正交归一的基向量，使得其中一个恰好沿曲线切向，另两个垂直于曲线。

曲线自己的“内蕴度规”

曲线作为一维流形，如果我们用弧长 $s$ 作为坐标，自然基底 $\frac{\partial}{\partial s}$ 满足：

$g\Big(\frac{\partial}{\partial s}, \frac{\partial}{\partial s}\Big) = 1 \quad \text{（因为 }ds^2 = 1\cdot ds^2\text{）}$

所以自然基底本身就是单位长度的, 定义:

$\boxed{\mathbf{T}(s) := \frac{\partial}{\partial s}}$

在一维曲线（比如一根弯弯曲曲的铁丝）上, 空间是弯曲的, 位置向量这个概念再也无法被"请进来"，但可以定义速度向量和加速度向量:

曲线上点 P 完全由一个参数 s=s(t) 来唯一确定:

$\mathbf{v}_{\text{内蕴}} = \frac{ds}{dt}\,\frac{\partial}{\partial s} = \dot{s}\,\mathbf{T}$

$\mathbf{a}_{\text{内蕴}} = \nabla_{\mathbf{v}}\mathbf{v} = \frac{d}{dt}(\dot{s})\,\mathbf{T} + \dot{s}\,\underbrace{\nabla_{\mathbf{v}}\mathbf{T}}_{=0} = \ddot{s}\,\mathbf{T}$

难点: 涉及到弯曲流形上的协变导数.

在高维外空间里“填充”出完整标架

我们需要再找两个与 $\mathbf{T}$ 正交且互相正交的单位向量 $\mathbf{N}, \mathbf{B}$，它们都在 $T_p \mathbb{R}^3$ 里，但不属于 $T_p C$。 方法: 标架向量对弧长 s 求导:

因为 $\mathbf{T}$ 是单位向量，$\frac{d\mathbf{T}}{ds}$ 必定垂直于 $\mathbf{T}$。
定义曲率 $\kappa(s) = \left\|\frac{d\mathbf{T}}{ds}\right\|$，并令主法向量：

$$\mathbf{N} := \frac{1}{\kappa}\frac{d\mathbf{T}}{ds}$$

同理可以继续求 $\frac{d\mathbf{N}}{ds}$，它必定垂直于 $\mathbf{N}$，可以分解成 $\mathbf{T}$ 和另一个垂直于 $\mathbf{T},\mathbf{N}$ 的方向 $\mathbf{B}$ 的组合。定义挠率 $\tau$ 使得：
$\frac{d\mathbf{B}}{ds} = -\tau \mathbf{N}$

这整套步骤与“求物理基底对坐标的偏导数”逻辑完全一致。得到的 Frenet–Serret 公式：

$$\frac{d}{ds} \begin{pmatrix} \mathbf{T} \\ \mathbf{N} \\ \mathbf{B} \end{pmatrix} = \begin{pmatrix} 0 & \kappa & 0\\ -\kappa & 0 & \tau\\ 0 & -\tau & 0 \end{pmatrix} \begin{pmatrix} \mathbf{T} \\ \mathbf{N} \\ \mathbf{B} \end{pmatrix}$$

本质上就是一维参数 $s$ 下的活动标架导数表，和您的圆柱坐标物理基底导数表：

$$\frac{\partial}{\partial \theta} \begin{pmatrix} \mathbf{e}_r \\ \mathbf{e}_\theta \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix} \begin{pmatrix} \mathbf{e}_r \\ \mathbf{e}_\theta \end{pmatrix}$$

区别仅在于：

• 极坐标是二维流形上的二维标架，参数有两个 $(r,\theta)$，您需要偏导数表。
• Frenet 是一维流形上的三维标架，参数只有一个 $s$，只需要全导数表。

	$d/ds$
$\mathbf{T}$	$\kappa \mathbf{N}$
$\mathbf{N}$	$-\kappa \mathbf{T} + \tau \mathbf{N}$
$\mathbf{B}$	$-\tau \mathbf{B}$

一维流形曲线上的位置, 速度, 和加速度

位置向量 — 假设为已知

设曲线用弧长 $s$ 参数化：

$$\mathbf{r} = \mathbf{r}(s).$$

如果我们改用时间 $t$ 作为参数，则位置向量写为

$$\boxed{\mathbf{r}(t) = \mathbf{r}(s(t))}$$

速度定义为位置对时间的导数:

$$\mathbf{v}(t) = \frac{d\mathbf{r}}{dt}$$

根据链式法则：

$$\boxed{\mathbf{v}(t) = \frac{d\mathbf{r}}{ds}\,\frac{ds}{dt}=v\mathbf{T}}$$

加速度定义:

$$\mathbf{a}(t) = \frac{d\mathbf{v}}{dt} = \frac{d}{dt}\bigl(v\,\mathbf{T}\bigr)$$

我们用乘积法则：

$$\mathbf{a} = \frac{dv}{dt}\,\mathbf{T} + v\,\frac{d\mathbf{T}}{dt}.$$

用链式法则：

$$\frac{d\mathbf{T}}{dt} = \frac{d\mathbf{T}}{ds}\,\frac{ds}{dt} = \kappa \mathbf{N}v$$

于是加速度成为：

$$\boxed{\mathbf{a} = \frac{dv}{dt}\,\mathbf{T} \;+\; \kappa\,v^{2}\,\mathbf{N}}.$$

解读加速度分量

• 切向分量 (tangential component)：$\displaystyle a_T = \frac{dv}{dt}$
  这是速率改变引起的加速度，方向沿 $\mathbf{T}$。
  若速率不变（匀速运动），此项为 $0$。

• 法向分量 (normal component)：$\displaystyle a_N = \kappa\,v^{2}$
  这是方向改变引起的加速度，方向沿主法向量 $\mathbf{N}$（指向弯曲内侧）。
  对于平面曲线，这就是“向心加速度”，其大小可以写成 $\frac{v^{2}}{R}$，其中 $R = \frac{1}{\kappa}$ 是曲率半径 (radius of curvature)。

注意：加速度向量完全没有沿着副法向量 $\mathbf{B}$ 的分量！
这是因为速度只沿 $\mathbf{T}$，而 $\mathbf{T}$ 的变化只产生在 $\mathbf{N}$ 方向，所以加速度始终在 $\mathbf{T}$ 和 $\mathbf{N}$ 张成的密切平面 (osculating plane) 内。