[三] 测地线方程：历史、推导与现实应用

引言：最短路径的直觉

想象你在一颗巨大星球上行走：从赤道的一个点出发，走到另一个经度的点。如果你沿着纬线走，那只是最省力的直觉，却未必是最短的距离。真正的最短路径，是那条看起来“弯着走”的大圆弧。

在平坦世界里，最短路径就是直线。但在弯曲空间中，“最短”和“最直”变成了一个更深的几何问题：测地线。测地线方程是一条连接历史、数学与现实的主线，它告诉我们：自由运动的轨迹在曲率中如何被重新定义。

第一章：测地线到底是什么

测地线（geodesic）可以从两个角度理解：

1. 几何角度：曲面或流形上“最直”的曲线，即切向量沿自身平行移动。
2. 变分角度：使弧长泛函取极值的曲线。

设曲线由参数 $ \lambda$ 描述：

$$x^i = x^i(\lambda), \quad i=1,\dots,n$$

弧长为：

$$S = \int_{\lambda_1}^{\lambda_2} ds = \int_{\lambda_1}^{\lambda_2} \sqrt{g_{ij}(x) \dot{x}^i \dot{x}^j} , d\lambda$$

让 $S$ 取极值的曲线，就是测地线。

第二章：测地线方程的历史脉络

测地线的故事几乎和微积分一样古老。

2.1 17-18世纪：变分法的萌芽

• 1697，伯努利：研究凸曲面最短路径，提出几何条件。
• 1732，欧拉：给出隐式曲面的测地线方程。
• 1744，欧拉《Methodus inveniendi》：系统建立变分法。
• 1788，拉格朗日《解析力学》：发展欧拉-拉格朗日方程，为测地线提供通用框架。

2.2 19世纪：几何语言的形成

• 1854，黎曼：引入度量张量，奠定弯曲空间几何基础。
• 1869，克里斯托费尔：提出克里斯托费尔符号，描述坐标基的变化。
• 1896，里奇与列维-奇维塔：形成绝对微分学与协变导数。
• 1917，列维-奇维塔：以平行移动解释协变导数，测地线获得清晰几何意义。

2.3 20世纪：物理的舞台

• 1915，爱因斯坦：将测地线方程作为自由落体的运动定律。
• 由此，测地线不仅属于几何，也成为引力理论的核心。

第三章：测地线方程的完整推导

3.1 变分原理

我们从弧长泛函开始：

$$S = \int \sqrt{g_{ij} \dot{x}^i \dot{x}^j} , d\lambda$$

由于平方根带来计算困难，我们使用等价的作用量：

$$S’ = \frac{1}{2} \int g_{ij} \dot{x}^i \dot{x}^j , d\lambda$$

$S$ 与 $S’$ 具有相同的极值曲线，只要参数是仿射的。

3.2 欧拉-拉格朗日方程

令拉格朗日量：

$$L = \frac{1}{2} g_{ij}(x) \dot{x}^i \dot{x}^j$$

欧拉-拉格朗日方程是：

$$\frac{d}{d\lambda}\left(\frac{\partial L}{\partial \dot{x}^k}\right) - \frac{\partial L}{\partial x^k} = 0$$

计算：

$$\frac{\partial L}{\partial \dot{x}^k} = g_{kj} \dot{x}^j$$

$$\frac{\partial L}{\partial x^k} = \frac{1}{2} \partial_k g_{ij} \dot{x}^i \dot{x}^j$$

代入后得：

$$\frac{d}{d\lambda}(g_{kj} \dot{x}^j) - \frac{1}{2} \partial_k g_{ij} \dot{x}^i \dot{x}^j = 0$$

展开第一项：

$$g_{kj} \ddot{x}^j + (\partial_l g_{kj}) \dot{x}^l \dot{x}^j - \frac{1}{2} \partial_k g_{ij} \dot{x}^i \dot{x}^j = 0$$

整理指标并乘以 $g^{km}$，得到：

$$\ddot{x}^m + \Gamma^m_{ij} \dot{x}^i \dot{x}^j = 0$$

其中：

$$\Gamma^m_{ij} = \frac{1}{2} g^{mk}(\partial_i g_{kj} + \partial_j g_{ki} - \partial_k g_{ij})$$

这就是测地线方程。

3.3 几何表述

用协变导数描述更直观：

$$\nabla_{\dot{\gamma}} \dot{\gamma} = 0$$

测地线是切向量沿自身平行移动的曲线，这与“最直”完全一致。

第四章：一个可计算的例子——球面

在半径为 $R$ 的球面上，用 $(\theta, \phi)$ 表示位置：

$$ds^2 = R^2(d\theta^2 + \sin^2\theta , d\phi^2)$$

非零的克里斯托费尔符号为：

$$\Gamma^\theta_{\phi\phi} = -\sin\theta\cos\theta, \quad \Gamma^\phi_{\theta\phi} = \Gamma^\phi_{\phi\theta} = \cot\theta$$

因此测地线方程是：

$$\theta’’ - \sin\theta\cos\theta , \phi’^2 = 0$$

$$\phi’’ + 2\cot\theta , \theta’ \phi’ = 0$$

这组方程的解是球面上的大圆。换句话说，飞机与卫星的最佳航线，正是球面测地线。

第五章：现实世界的测地线

5.1 地理导航与测绘

地球表面最短路径是大圆航线，航空与航运路线规划直接依赖测地线。

5.2 广义相对论

自由落体物体满足：

$$\frac{d^2 x^\mu}{d\tau^2} + \Gamma^\mu_{\alpha\beta} \frac{dx^\alpha}{d\tau} \frac{dx^\beta}{d\tau} = 0$$

这说明引力不是力，而是时空的曲率。光的偏折、行星近日点进动，都来自测地线。

5.3 机器人与计算机图形学

• 机器人手臂的姿态空间是流形，最省力的姿态变化就是测地线。
• 在曲面网格上寻找最短路径时，数值方法往往离散测地线方程。

5.4 机器学习与数据分析

在流形学习中，测地线距离比欧氏距离更能反映真实的“数据曲率”。Isomap 与图嵌入方法正是通过近似测地线来发现低维结构。

结语：最直的曲线

测地线方程把“最短路径”提升为“最直运动”的几何定律。从欧拉的变分法到爱因斯坦的时空，从地球的航线到机器人的姿态规划，它是一条横跨三百年的数学曲线。

真正的直线，从来不只是一根尺子画出来的，它是一个空间对自身的回应。

附录：重要公式汇总

• 弧长泛函：$S = \int \sqrt{g_{ij} \dot{x}^i \dot{x}^j} , d\lambda$
• 等价作用量：$S’ = \frac{1}{2} \int g_{ij} \dot{x}^i \dot{x}^j , d\lambda$
• 测地线方程：$\ddot{x}^m + \Gamma^m_{ij} \dot{x}^i \dot{x}^j = 0$
• 克里斯托费尔符号：$\Gamma^m_{ij} = \frac{1}{2} g^{mk}(\partial_i g_{kj} + \partial_j g_{ki} - \partial_k g_{ij})$
• 几何形式：$\nabla_{\dot{\gamma}} \dot{\gamma} = 0$

系列导航

本文是广义相对论系列文章的第 [3] 篇。

本系列文章：