曲率与测地线

1. 曲率的定义

1.1 曲线的曲率

• 曲线的曲率定义：
  曲率描述了曲线的弯曲程度。在二维平面中，曲线的曲率可以定义为曲线在某一点的变化率，即曲线的法向量相对于该点的变化速度。

  • 数学公式：
    对于一条参数化曲线 $\gamma (t)=(x(t),y(t))$，曲率 $\kappa (t)$ 定义为：
    $$\kappa (t)=\frac{|{\gamma }^{\prime \prime }(t)|}{(1+({\gamma }^{\prime }(t){)}^2{)}^{3/2}}$$
    其中，(\gamma’(t)) 是曲线的一阶导数，(\gamma’'(t)) 是二阶导数。

  • 例子：
    对于单位圆 $\gamma (t)=(\cos (t),\sin (t))$，计算曲率。由公式可得，对于单位圆，曲率为常数 $\kappa =1$。

1.2 曲面的高斯曲率

• 高斯曲率的定义：
  高斯曲率 $K$ 描述了曲面在某一点的弯曲程度。它是由该点的主曲率的乘积给出的。

  • 数学公式：
    对于一个光滑的曲面，假设 $k_1$ 和 $k_2$ 是该点处的主曲率（曲面在该点的最大和最小弯曲度），则高斯曲率定义为：
    $$K=k_1\cdot k_2$$

  • 例子：
    对于球面，假设球的半径为 (R)，球面的主曲率 (k_1 = k_2 = \frac{1}{R})，因此球面上的高斯曲率 $K=\frac{1}{R^2}$。

1.3 局部与全局曲率的关系

• 局部曲率与全局曲率：
  局部曲率描述了曲线或曲面在某一点的弯曲程度，而全局曲率则考虑了曲线或曲面的整体形状。在微分几何中，全局性质如总曲率是通过积分局部曲率得到的。

2. 测地线

2.1 测地线的定义

• 测地线简介：
  测地线是流形上最短路径的自然推广。它是在流形上局部直线的类比，类似于平面中的直线。测地线是流形上弯曲的路径，满足“沿路径的加速度为零”的条件。

  • 数学公式：
    测地线方程是一个二阶常微分方程，通常写成：
    $$\frac{d^2{x}^i}{d{t}^2}+{\Gamma }_{jk}^i\frac{d{x}^j}{dt}\frac{d{x}^k}{dt}=0$$
    其中 ${\Gamma }_{jk}^i$ 是流形的克里斯托费尔符号，$x^i$ 是流形上的坐标，$t$ 是参数。

2.2 测地线的几何意义

• 几何意义：
  测地线的几何意义是，在流形上，沿着测地线移动时，物体不会偏离其原来的方向，即沿路径的加速度为零。测地线是最短路径，但在曲面上，测地线的长度不仅取决于路径本身，还与曲面的形状和弯曲程度有关。

2.3 练习案例

• 例子：
  计算球面上的大圆路径（即球面上的测地线）。设球面半径为 (R)，球面上的两点 A 和 B 在大圆上的距离可以通过测地线方程进行求解。对于球面，测地线是大圆。

3. 极小曲面与最短路径

3.1 极小曲面

• 极小曲面简介：
  极小曲面是曲面的一类，它在其每一点的高斯曲率为零，通常是通过局部最小化曲面面积来定义。常见的例子有双曲面和悬链线。

  • 数学公式：
    极小曲面满足以下方程：
    $$H=0$$
    其中 $H$ 是曲面的平均曲率。

3.2 最短路径

• 最短路径的定义：
  最短路径是两点之间曲面上的路径，其路径长度最小。测地线就是最短路径的数学模型。

  • 练习案例：
    计算球面上一条给定路径的长度，验证其是否为测地线（最短路径）。路径的长度可以通过积分公式计算：
    $$L={\int }_{\gamma }\sqrt{g_{ij}\frac{d{x}^i}{dt}\frac{d{x}^j}{dt}}dt$$
    其中 $g_{ij}$ 是流形的度量张量，$\gamma$ 是路径。

4. 课堂活动与案例

活动1：计算二维流形（如球面）的高斯曲率

• 任务： 学生通过计算球面上的高斯曲率来理解曲率的局部性质。

  • 步骤：

    1. 给定球面半径 (R)，计算其主曲率。
    2. 计算球面上的高斯曲率 $K$。
    3. 分析球面上不同点的曲率变化。

  • 答案： 对于半径为 (R) 的球面，主曲率 (k_1 = k_2 = \frac{1}{R})，因此高斯曲率为 $K=\frac{1}{R^2}$。

活动2：测地线的性质与物理学中的应用

• 任务： 学生讨论测地线的性质，并通过实例展示其在物理学中的应用，特别是测地线在引力场中的作用。

  • 步骤：
    1. 讨论在一个弯曲的时空中如何利用测地线描述物体的自由下落路径。
    2. 通过计算简单的引力场中的测地线来说明其几何意义。

Python代码实现示例

计算球面上的高斯曲率

import numpy as np
import matplotlib.pyplot as plt

# 球面半径
R = 1

# 球面上的高斯曲率 K = 1 / R^2
K = 1 / R**2
print(f"球面上的高斯曲率: {K}")

# 可视化球面上的一部分
theta = np.linspace(0, np.pi, 100)
phi = np.linspace(0, 2*np.pi, 100)
Theta, Phi = np.meshgrid(theta, phi)
X = R * np.sin(Theta) * np.cos(Phi)
Y = R * np.sin(Theta) * np.sin(Phi)
Z = R * np.cos(Theta)

# 绘制球面
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(X, Y, Z, rstride=5, cstride=5, color='b', alpha=0.6)
ax.set_title("球面")
plt.show()

计算球面上的测地线

from scipy.integrate import odeint

# 球面上的测地线方程
def geodesic_eqn(y, t, R):
    x, y_val, z = y
    dxdt = -np.sin(x) * np.cos(y_val)  # x方向的速度
    dydt = np.cos(x) * np.sin(y_val)  # y方向的速度
    dzdt = np.cos(x) * np.cos(y_val)  # z方向的速度
    return [dxdt, dydt, dzdt]

# 初始条件
R = 1
y0 = [0, 0, R]

# 时间范围
t = np.linspace(0, np.pi, 100)

# 解测地线方程
sol = odeint(geodesic_eqn, y0, t, args=(R,))

# 可视化测地线
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot(sol[:, 0], sol[:, 1], sol[:, 2], color='r')
ax.set_title("球面上的测地线")
plt.show()

总结

这节课涵盖了曲率与测地线的基本概念，详细解释了曲率的定义、测地线方程、极小曲面与最短路径的关系。每个部分都提供了数学公式、实例计算以及Python代码，帮助学生深入理解曲率和测地线的几何性质及其应用。