跳到正文
格致开物MATHWIKI

测地线

测地线(geodesic)是曲面或更一般的黎曼流形上,速度向量沿曲线自身保持平行的参数曲线。对欧氏空间中的光滑曲面,这意味着曲线的空间加速度没有切向分量,只有法向分量。它可以在三维空间中弯曲,例如圆柱上的等速螺线;“直”指的是曲面内部的几何,而不是空间中的直尺形状。

测地线与测地距离有关,但两者不是同一个对象。足够短的测地线段具有最短路性质;延长后的整段可能已不是两端之间的最短路径。参数也属于定义的一部分:沿同一路线忽快忽慢地运动,一般不再满足测地线方程。

圆柱上的路线与加速度

先看半径为一的无限长圆柱面 X(u,z)=(cos⁡u,sin⁡u,z),u,z∈ℝ. 这里 u 是绕轴转过的角,z 是高度。参数 u 相差 2π 时对应同一位置,因而这个式子是周期参数表示,不是整个圆柱的一张一一对应坐标图;限制角度区间后才得到局部坐标。

取曲线 γ(t)=(cos⁡(at+b),sin⁡(at+b),ct+d). 其中 a,b,c,d 为常数。求两次导数: γ′(t)=(−asin⁡(at+b),acos⁡(at+b),c), γ″(t)=−a2(cos⁡(at+b),sin⁡(at+b),0). 圆柱的外单位法向是 n=(cos⁡u,sin⁡u,0)。因此 γ″=−a2n,加速度完全垂直于曲面;速度大小则恒为 a2+c2。当 a,c 都非零时得到螺线,a=0 时是母线,c=0 时是水平圆。这三种曲线在所给等速参数下都是测地线。一个水平圆能是圆柱的测地线,并不意味着它是包含该圆的水平平面的测地线:所用曲面不同,允许的切向方向也不同。

单位圆柱上的四分之一圈螺线,由A到B升高二,选定点的加速度沿圆柱内法向
曲线从 A=(1,0,0) 到 B=(0,1,2)。金色箭头表示选定点的加速度方向;箭头长度作示意,不表示加速度大小。

这个例子也说明速度保持不变还不够。在平面上匀速绕圆运动时,向心加速度仍在该平面的切空间内,所以圆不是平面的测地线。

协变导数与测地线方程

曲面上不同点的切平面通常不同,不能把不同切平面中的速度分量直接相减而认为已经得到曲面内的变化。对欧氏空间中的曲面,先在空间中求导,再投影回当前切平面,得到沿曲线的协变导数。若 V(t) 是沿曲线的切向量,则 DVdt=projTγ(t)MV′(t). 在一般光滑黎曼流形上,使用由度量确定的 Levi-Civita 联络。它给出同样的“沿曲线比较切向量”的求导规则;该规则与内积相容,且无挠。本文中的测地线均采用这一联络。[1]

定义。 设 (M,g) 为光滑黎曼流形,I⊆ℝ 为区间。二次连续可微曲线 γ:I→M 满足 Dγ′dt=∇γ′γ′=0 时,称为测地线。常值曲线也包含在定义内;非恒定测地线的速度不会降到零。

用局部坐标 x1,…,xm 计算时,令 ∂i=∂/∂xi 表示坐标方向的切向量,写 ∇∂i∂j=∑k=1mΓijk∂k. 系数 Γijk 称为 Christoffel 系数。它们记录基向量本身怎样随位置变化,并不是额外施加的力。若 gij=g(∂i,∂j),而 (gij) 为矩阵 (gij) 的逆,则 Γijk=12∑ℓ=1mgkℓ(∂igjℓ+∂jgiℓ−∂ℓgij). 下标 i,j 的对称性对应无挠性。这个表达式来自联络的度量相容性和无挠性;此处用它把已有度量转成可计算的微分方程。

写 γ′(t)=∑jx˙j(t)∂j。对系数和基向量分别求导,得到 Dγ′dt=∑k(x¨k+∑i,jΓijk(x(t))x˙ix˙j)∂k. 因此测地线方程为 x¨k+∑i,j=1mΓijk(x(t))x˙ix˙j=0,k=1,…,m. 它是二阶常微分方程组。光滑系数保证给定初始位置 p 与初始速度 v∈TpM 后,存在唯一的局部解,并可延拓成唯一的最大解;这项结论使用常微分方程的局部存在唯一性定理,不保证最大解对所有实数时间都有定义。

恒速与仿射参数

度量相容性允许像普通内积那样求导: ddtg(γ′,γ′)=2g(Dγ′dt,γ′)=0. 所以测地线恒速;若初始速度非零,就始终非零。这里证明了“测地线必恒速”,没有证明反命题。

再令 β(s)=γ(φ(s)),其中 φ 二次连续可微且 φ′(s)≠0。链式法则给出 Dβ′ds=(φ′)2Dγ′dt+φ″γ′=φ″γ′. 对非恒定测地线,最后一式为零当且仅当 φ″=0,即 φ(s)=As+B、A≠0。因此保留测地线方程的正则换参恰为仿射换参。按弧长行进是其中一种常用选择。

同一条平面直线上等参数间隔采样,匀速点间距相等,平方换参点间距逐渐增大
上排为 γ(t)=(t,0),下排为 β(s)=(s²,0),取相同的四分之一参数间隔。下排几何轨迹仍是直线,但 β=(2,0),不满足仿射参数的测地线方程;正则性比较取 0<s<1。

在欧氏平面的直角坐标中,度量为常数矩阵,所有 Christoffel 系数都为零,方程退化为 x¨=y¨=0,所以解是 (x,y)=(a1t+b1,a2t+b2)。这正是带恒定速度的直线。若换成极坐标,度量变成 dr2+r2dθ2,虽然平面没有因此弯曲,方程却写成 r¨−rθ˙2=0,θ¨+2rr˙θ˙=0(r>0). 这里 Γθθr=−r、Γrθθ=Γθrθ=1/r。非零系数可以只来自坐标基的变化,不能单凭它们判断空间是否具有曲率。

圆柱展开与绕行次数

回到半径为 R>0 的无限长圆柱 X(u,z)=(Rcos⁡u,Rsin⁡u,z). 求导得到 Xu=(−Rsin⁡u,Rcos⁡u,0)、Xz=(0,0,1),内积为 guu=R2,guz=0,gzz=1. 因此长度元为 ds2=R2du2+dz2。取展开坐标 s=Ru 后成为 ds2+dz2,在不拉伸圆柱的局部展开中,长度与平面相同;测地线对应展开平面的等速直线。这也由所有系数恒定、测地线方程 u¨=z¨=0 直接推出。

给定 A=X(0,0)、B=X(α,h),其中可取 −π≤α≤π,则 B 在展开平面有无穷多个代表点 Bk=(R(α+2πk),h),k∈ℤ. 整数 k 区分绕轴的路线类别。沿任何从 A 出发的连续圆柱路径,都可连续选择角度:在每个足够小的局部片上选一支角度,跨片时加减 2π 使数值连续。由于起点角度固定为零,终点角度只能是 α+2πk。对分段光滑路径,这样得到的展开路径也分段光滑,且长度保持不变。

展开后连接 (0,0) 与 Bk 的任意路径 c 满足 L(c)=∫‖c′(t)‖dt≥‖∫c′(t)dt‖=R2(α+2πk)2+h2. 等速线段达到这个下界,卷回去得到 γk(t)=X((α+2πk)t,ht),0≤t≤1. 故圆柱上的内蕴距离为 d(A,B)=mink∈ℤR2(α+2πk)2+h2. 这给出了该圆柱情形的完整最短路证明:既比较了同一绕行类别内的路径,也比较了所有可能的绕行类别。它没有把“某一条曲线满足测地线方程”直接当作全局最短的证明。

展开平面中起点A连向同一个圆柱终点的三个代表点B负一、B零、B一,横坐标分别负三π除二、π除二、五π除二,高度都为二
单位圆柱上同一对端点对应三条展开直线。实线的绕行量最小;两条虚线仍对应测地线,但整段较长。

可复算的端点例子

取 R=1、α=π/2、h=2,即前图的 A=(1,0,0)、B=(0,1,2)。三条路径的长度分别为 L0=π2/4+4≈2.543109, L−1=9π2/4+4≈5.119239,L1=25π2/4+4≈8.104630. 三者都满足测地线方程,只有 k=0 给出全局最短路径。这里坐标中的圆柱半径与高度采用同一长度单位。

变式。 把终点改为 B=(−1,0,2),其角度取 α=π。此时 k=0 与 k=−1 的角度变化分别为 π 和 −π,两条不同的最短螺线长度都为 π2+4≈3.724192。其他整数给出绝对值更大的角度变化,长度严格更大。这是最短路不唯一的具体例子。

局部最短、延拓与适用范围

在光滑黎曼流形内部,每个点都有足够小的正常邻域,使从该点出发的适当短测地线段实现该邻域内两端之间的距离。一般证明使用指数映射和 Gauss 引理,此处不由圆柱的特殊展开推广得出。[2] “足够短”依赖所在位置和路线。圆柱上绕轴多圈的例子已经表明,局部最短不能直接延伸成任意长段的全局最短。

测地线能延拓多久还依赖空间是否完整。例如在带欧氏度量的开圆盘中,从原点沿径向匀速行进,有限时间就到达未包含在空间内的边界,因而不能继续作为该空间中的曲线。若把边界也加入,允许沿边界行走的最短路问题又需要单独说明约束,不能直接套用无边界流形内部的方程。

与平均曲率比较时也应区分对象:测地线描述曲面内部的曲线路径;平均曲率描述曲面本身在环境空间中的弯曲。圆柱与平面局部等距,测地线可通过展开对应,而它们的平均曲率并不相同。

参考资料

  1. ↑ Ross Ogilvie,Riemannian Geometry,Universität Mannheim,2024 年春季,第 4 章,§4.1 定义 4.1、定理 4.5、练习 4.7 及嵌入曲面的切向投影;另见第 3 章 §3.6 的 Levi-Civita 联络与 Christoffel 系数。访问日期:2026-10-25。
  2. ↑ Ross Ogilvie,前引讲义第 4 章 §4.4“Exponential Maps”,引理 4.31 与定理 4.34。