曲线求长与不可求长曲线（海岸线居然算不出长度？）

曲线求长与不可求长曲线

小象实战讲义 · 人工智能数学基础

定积分不仅是计算面积和体积的利器，它还能帮助我们解决一个直观而深刻的问题：如何精确度量一条曲线的长度？本节将系统学习光滑曲线弧长的计算公式，从参数方程到直角坐标再到极坐标，并揭示一个反直觉的数学事实：并非所有看起来“有长度”的曲线都能被有限度量，由此我们将触及“分形”这一现代数学的迷人领域。

💡 核心导读

• 光滑曲线与弧长定义：理解曲线可求长的前提——光滑性（一阶导函数连续），并掌握用折线逼近曲线求弧长的基本思想。
• 三大弧长公式：掌握参数方程、直角坐标方程和极坐标方程下的弧长（弧微分）计算公式及其推导。
• 不可求长曲线：认识科赫雪花曲线这类具有“自相似性”的分形，理解其周长无限而面积有限的悖论，了解“分数维”概念。
• 典型计算：学习旋轮线、椭圆等经典曲线的弧长计算，了解“椭圆积分”这一非初等函数表达。

光滑曲线的弧长

1. 问题与定义

我们希望计算一条平面曲线 $L$ 的长度。如果这条曲线可以由参数方程描述： $$\begin{cases} x = x(t) \ y = y(t) \end{cases} \quad \alpha \le t \le \beta$$ 并且函数 $x(t)$ 和 $y(t)$ 在 $[\alpha, \beta]$ 上具有连续的一阶导函数（即 $x(t), y(t) \in C^1$），那么我们称这样的曲线为光滑曲线。光滑曲线的切线方向是连续变化的，这是其长度可以被精确定义和计算的关键前提。

2. 弧长公式的推导（参数方程情形）

计算思路源于定积分的基本思想：分割、近似、求和、取极限。

1. 分割：将参数区间 $[\alpha, \beta]$ 用分点 $t_0=\alpha < t_1 < \dots < t_n=\beta$ 任意分割成 $n$ 个小区间。
2. 近似：对应地，曲线上也得到一系列点 $P_i(x(t_i), y(t_i))$。用连接相邻点的弦（直线段）来近似代替对应的曲线弧。第 $i$ 段弦的长度为： $$|P_{i-1}P_i| = \sqrt{[x(t_i) - x(t_{i-1})]^2 + [y(t_i) - y(t_{i-1})]^2}$$
3. 求和与取极限：所有弦长的总和 $S_n = \sum_{i=1}^{n} |P_{i-1}P_i|$ 是曲线弧长的一个近似。当分割越来越细（即 $\max \Delta t_i \to 0$）时，如果 $S_n$ 的极限存在且有限，则该极限值就定义为曲线的弧长 $S$。

为了计算这个极限，利用微分进行近似。因为 $x(t), y(t)$ 可导，当 $\Delta t_i = t_i - t_{i-1}$ 很小时，有： $$\Delta x_i = x(t_i) - x(t_{i-1}) \approx x’(t_{i-1}) \Delta t_i, \quad \Delta y_i = y(t_i) - y(t_{i-1}) \approx y’(t_{i-1}) \Delta t_i$$ 代入弦长公式： $$|P_{i-1}P_i| \approx \sqrt{[x’(t_{i-1}) \Delta t_i]^2 + [y’(t_{i-1}) \Delta t_i]^2} = \sqrt{[x’(t_{i-1})]^2 + [y’(t_{i-1})]^2} \ \Delta t_i$$ 于是，弧长 $S$ 就是如下黎曼和的极限： $$S = \lim_{\max \Delta t_i \to 0} \sum_{i=1}^{n} \sqrt{[x’(t_{i-1})]^2 + [y’(t_{i-1})]^2} \ \Delta t_i$$ 根据定积分的定义，我们得到参数方程下的弧长公式：

$$\boxed{S = \int_{\alpha}^{\beta} \sqrt{[x’(t)]^2 + [y’(t)]^2} \ dt}$$

其中，被积表达式 $\sqrt{[x’(t)]^2 + [y’(t)]^2} \ dt$ 称为弧微分，记作 $ds$。因此弧长公式也常写作 $S = \int_L ds$。

3. 特殊情形：直角坐标方程

若曲线由直角坐标方程 $y = f(x) \ (a \le x \le b)$ 给出，可将其视为参数方程的特例：$x = t, \ y = f(t)$。代入上述公式，得到直角坐标下的弧长公式：

$$\boxed{S = \int_{a}^{b} \sqrt{1 + [f’(x)]^2} \ dx}$$

4. 特殊情形：极坐标方程

若曲线由极坐标方程 $r = r(\theta) \ (\alpha \le \theta \le \beta)$ 给出。利用极坐标与直角坐标的转换关系： $$x = r(\theta)\cos\theta, \quad y = r(\theta)\sin\theta$$ 将 $\theta$ 视为参数，计算 $x’(\theta)$ 和 $y’(\theta)$： $$\begin{aligned} x’(\theta) &= r’(\theta)\cos\theta - r(\theta)\sin\theta \ y’(\theta) &= r’(\theta)\sin\theta + r(\theta)\cos\theta \end{aligned}$$ 则 $$[x’(\theta)]^2 + [y’(\theta)]^2 = [r’(\theta)]^2 + [r(\theta)]^2$$ 代入参数方程弧长公式，得到极坐标下的弧长公式：

$$\boxed{S = \int_{\alpha}^{\beta} \sqrt{[r(\theta)]^2 + [r’(\theta)]^2} \ d\theta}$$

import numpy as np
import sympy as sp

# 示例1：使用sympy验证直角坐标下圆弧弧长公式
# 计算上半圆 y = sqrt(R^2 - x^2) 在 [-R, R] 上的弧长
x, R = sp.symbols('x R', positive=True)
y = sp.sqrt(R**2 - x**2)
y_prime = sp.diff(y, x)
arc_length_expr = sp.sqrt(1 + y_prime**2)
# 尝试积分（应得到 πR）
S = sp.integrate(arc_length_expr, (x, -R, R))
print("上半圆弧长公式推导结果：", S.simplify())

# 示例2：数值计算一段旋轮线的弧长
# 旋轮线参数方程: x = r*(t - sin(t)), y = r*(1 - cos(t)), t in [0, 2π]
def cycloid_arc_length(r, n=10000):
    t = np.linspace(0, 2*np.pi, n)
    x_prime = r * (1 - np.cos(t))  # dx/dt
    y_prime = r * np.sin(t)        # dy/dt
    ds = np.sqrt(x_prime**2 + y_prime**2)
    dt = 2*np.pi / (n-1)
    S_numerical = np.sum(ds) * dt  # 矩形法近似积分
    return S_numerical

r = 2.0
S_num = cycloid_arc_length(r)
S_exact = 8 * r  # 旋轮线一拱弧长的精确公式
print(f"\n旋轮线（r={r}）一拱弧长：")
print(f"  数值积分结果：{S_num:.6f}")
print(f"  理论精确值：{S_exact:.6f}")
print(f"  绝对误差：{abs(S_num - S_exact):.6f}")

反直觉的发现：不可求长曲线

1. 问题的提出

根据黎曼积分的定义，曲线可求长的本质是：用内接折线去逼近曲线时，折线总长度序列的极限存在且有限。然而，并非所有我们直觉上认为“有长度”的曲线都满足这个条件。

2. 科赫雪花曲线

科赫曲线是一个经典的反例，它通过一个简单的递归规则生成：

1. 从一条长度为 $L_0$ 的线段开始。
2. 将线段三等分，去掉中间一段，并用一个等边三角形的两条边替代（凸起）。
3. 对新生图形的每一条边，重复步骤2。

经过无限次迭代后得到的极限曲线，就是科赫曲线。将其首尾相连，即得到科赫雪花。

关键性质：

• 自相似性：图形的局部与整体具有相似的形态。
• 周长无限：设初始正三角形边长为 $a$，周长为 $L_0 = 3a$。每迭代一次，每条边长度变为原来的 $4/3$ 倍。因此，第 $n$ 级图形的周长 $L_n = 3a \times (4/3)^n$。当 $n \to \infty$ 时，$L_n \to \infty$。
• 面积有限：可以证明，科赫雪花所围成的面积是一个有限值 $S = \frac{2\sqrt{3}}{5}a^2$。
• 处处不可导：曲线极其曲折，在任何一点都没有切线。
• 分数维：其分形维数 $D = \frac{\log 4}{\log 3} \approx 1.262$，大于1（曲线）而小于2（面）。

3. 海岸线问题与分形

科赫曲线揭示的“有限面积，无限周长”现象，并非纯粹的数学游戏。它直接关联到著名的“海岸线测量悖论”：不同精度下测量的海岸线长度会显著不同，测量尺度越小，测得的长度越长。数学家曼德博指出，海岸线具有统计意义上的自相似性，其长度依赖于测量尺度，在理论上，当尺度无限小时，长度趋于无穷。因此，用传统的“长度”来描述海岸线是不完备的，需要引入“分形维数”这一概念。

分形理论现已广泛应用于自然科学、计算机图形学、数据压缩等多个领域。科赫曲线是分形的一个标准模型。

典型曲线弧长计算示例

1. 旋轮线一拱的长度

旋轮线的参数方程为： $$\begin{cases} x = r(\theta - \sin\theta) \ y = r(1 - \cos\theta) \end{cases} \quad 0 \le \theta \le 2\pi$$ 其中 $r$ 是生成圆的半径。计算其一拱的弧长： $$\begin{aligned} S &= \int_{0}^{2\pi} \sqrt{[x’(\theta)]^2 + [y’(\theta)]^2} \ d\theta \ &= \int_{0}^{2\pi} \sqrt{[r(1-\cos\theta)]^2 + [r\sin\theta]^2} \ d\theta \ &= r \int_{0}^{2\pi} \sqrt{2 - 2\cos\theta} \ d\theta \ &= r \int_{0}^{2\pi} \sqrt{4\sin^2(\theta/2)} \ d\theta \quad (\text{利用半角公式}) \ &= 2r \int_{0}^{2\pi} |\sin(\theta/2)| \ d\theta \ &= 2r \int_{0}^{2\pi} \sin(\theta/2) \ d\theta \quad (\text{在 } [0, 2\pi] \text{上 } \sin(\theta/2) \ge 0) \ &= 2r \left[ -2\cos(\theta/2) \right]_{0}^{2\pi} \ &= 8r \end{aligned}$$ 因此，半径为 $r$ 的圆所生成的旋轮线，其一拱的弧长为 $8r$。

2. 椭圆的周长

椭圆的标准方程为 $\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1 \ (a \ge b > 0)$。利用对称性，只需计算第一象限的弧长再乘以4。 由方程可得第一象限函数表达式：$y = \frac{b}{a} \sqrt{a^2 - x^2}$。其弧微分： $$ds = \sqrt{1 + (y’)^2} dx = \sqrt{1 + \frac{b^2 x^2}{a^2(a^2 - x^2)}} dx = \sqrt{\frac{a^4 - (a^2 - b^2)x^2}{a^2(a^2 - x^2)}} dx$$ 令 $c^2 = a^2 - b^2$（$c$ 为半焦距），则椭圆周长 $S$ 为： $$S = 4 \int_{0}^{a} \sqrt{\frac{a^4 - c^2 x^2}{a^2(a^2 - x^2)}} dx$$ 这个积分无法用初等函数（有理函数、指数、对数、三角函数等）表示。历史上，这类积分被专门研究，称为椭圆积分。因此，椭圆周长没有像圆周长 $2\pi r$ 那样的简洁初等表达式，其精确值通常用级数或特殊函数表示。一个常用的近似公式是 $S \approx \pi [3(a+b) - \sqrt{(3a+b)(a+3b)}]$。

椭圆的面积则容易计算：$A = \pi ab$。

📝 动手练一练

1. 计算圆弧长：利用直角坐标下的弧长公式，验证半径为 $R$ 的上半圆 $y = \sqrt{R^2 - x^2} \ (-R \le x \le R)$ 的弧长为 $\pi R$。

   参考答案： 由 $y’ = \frac{-x}{\sqrt{R^2 - x^2}}$，则 $ds = \sqrt{1 + \frac{x^2}{R^2 - x^2}} dx = \frac{R}{\sqrt{R^2 - x^2}} dx$。 弧长 $S = \int_{-R}^{R} \frac{R}{\sqrt{R^2 - x^2}} dx = R \left[ \arcsin(\frac{x}{R}) \right]_{-R}^{R} = R(\frac{\pi}{2} - (-\frac{\pi}{2})) = \pi R$。

2. 计算心形线长度：求心形线 $r = a(1 + \cos\theta) \ (a>0)$ 的完整长度。

   参考答案： 心形线关于极轴对称。弧微分 $ds = \sqrt{r^2 + (r’)^2} d\theta = \sqrt{a^2(1+\cos\theta)^2 + a^2\sin^2\theta} \ d\theta = a\sqrt{2+2\cos\theta} \ d\theta = 2a|\cos(\theta/2)| d\theta$。 完整长度 $S = 2 \int_{0}^{\pi} 2a\cos(\theta/2) d\theta = 8a \left[ \sin(\theta/2) \right]_{0}^{\pi} = 8a$。

本章小结

本节我们深入探讨了定积分在几何上的另一核心应用——曲线求长。

要点回顾：

• 可求长条件：光滑曲线（$C^1$ 参数方程）的弧长可以通过积分精确定义和计算。
• 三大公式：参数方程 $S=\int \sqrt{x’^2+y’^2}dt$；直角坐标 $S=\int \sqrt{1+y’^2}dx$；极坐标 $S=\int \sqrt{r^2+r’^2}d\theta$。核心思想是弧微分 $ds$ 的积分。
• 认知突破：科赫雪花曲线揭示了存在“有限面积、无限周长”的图形，引出了分形与分数维的概念，打破了我们对维度和测量的传统认知。
• 计算实践：掌握了旋轮线等规则曲线弧长的计算，并了解到椭圆周长等积分无法用初等函数表达，属于椭圆积分范畴。

行动清单：

1. 推导一遍：亲手推导一遍从折线逼近到积分公式的完整过程，巩固对定积分定义的理解。
2. 编程验证：使用 Python（如 SymPy, NumPy）对讲义中的弧长公式进行符号推导和数值验证，加深对公式的记忆和应用能力。
3. 拓展阅读：搜索“分形”、“曼德博集”、“海岸线悖论”等关键词，了解这一数学分支的更多奇妙现象和应用。

— 小象教研组