4. 费马原理与光程：费马原理的表述、光程概念、由费马原理推导反射和折射定律

几何光学中，反射定律与折射定律是描述光传播行为的两条基本经验规则。然而，这两条定律背后隐藏着一个更为深刻且统一的物理原理——费马原理（Fermat's Principle）。费马原理以“光程”为核心概念，将光的传播路径问题转化为极值问题，不仅能够统一推导反射与折射定律，还为后续理解透镜成像、像差分析乃至光学系统设计提供了坚实的理论基础。本节将从费马原理的表述出发，系统讲解光程的物理意义，并严格推导反射与折射定律，帮助学习者建立从“现象规律”到“原理统一”的认知跃迁。

4.1 费马原理的表述与物理内涵

费马原理由法国数学家皮埃尔·德·费马于1657年提出，其最经典的表述为：光在任意两点之间传播时，总是沿着光程为极值（通常为极小值）的路径传播。更精确的数学表述为：光线的实际路径使光程的一阶变分为零，即 δL = 0，其中 L 为光程。这里的“极值”不仅包括极小值，也包括极大值和恒定值情形，因此更严谨的说法是“光程取驻定值”。

理解费马原理需要把握三个层次：

• 物理本质：光程 L = n·s，其中 n 为介质折射率，s 为几何路程。由于光速 v = c/n，光程 L = n·s = c·(s/v) = c·t，即光程正比于光传播所需的时间。因此，费马原理等价于“光沿所需时间最短（或极值）的路径传播”，这也被称为“最短时间原理”。
• 变分思想：费马原理将几何光学问题转化为变分问题——在所有可能的路径中，寻找使光程泛函取极值的那一条。这与力学中的最小作用量原理在数学结构上高度一致，体现了物理学中“自然选择最优”的哲学思想。
• 适用边界：费马原理适用于均匀介质、分层介质以及折射率连续变化的非均匀介质。对于各向同性介质，光程与路径无关；对于各向异性介质（如晶体），需使用更一般的光程定义。

4.2 光程的数学定义与计算方法

光程（Optical Path Length, OPL）是费马原理的核心物理量。在折射率为 n 的均匀介质中，光程定义为几何路程 s 与折射率 n 的乘积：L = n·s。若光线穿越多种介质，则总光程为各段光程之和：

L = Σ nᵢ · sᵢ   （分段均匀介质）

若介质折射率连续变化，即 n = n(x, y, z)，则光程需用线积分表示：

L = ∫ n(x, y, z) ds   （连续非均匀介质）

其中 ds 为路径上的微元弧长。在实际工程计算中，光程常被用来衡量光线经过不同路径后的相位延迟。由于相位变化 Δφ = (2π/λ₀)·L，其中 λ₀ 为真空中的波长，因此光程差（OPD）直接决定了干涉、衍射现象的条纹分布。

4.3 由费马原理推导反射定律

反射定律的推导是费马原理最直观的应用。考虑一个平面反射镜，设入射点 A 的坐标为 (x₁, y₁)，反射点 B 的坐标为 (x₂, y₂)，反射面为 x 轴（y = 0），光线在反射面上的入射点 P 的坐标为 (x, 0)。由于光线在均匀介质（空气，n ≈ 1）中传播，光程 L 等于几何路程之和：

L(x) = |AP| + |PB| = √[(x - x₁)² + y₁²] + √[(x₂ - x)² + y₂²]

根据费马原理，实际光线的路径使 L(x) 取极值，即 dL/dx = 0。对 L(x) 求导：

dL/dx = (x - x₁) / √[(x - x₁)² + y₁²] - (x₂ - x) / √[(x₂ - x)² + y₂²] = 0

令 dL/dx = 0，得到：

(x - x₁) / |AP| = (x₂ - x) / |PB|

观察几何关系可知，(x - x₁)/|AP| 恰好等于入射角 θᵢ 的正弦值 sinθᵢ，而 (x₂ - x)/|PB| 等于反射角 θᵣ 的正弦值 sinθᵣ。因此：

sinθᵢ = sinθᵣ  ⟹  θᵢ = θᵣ

同时，由于 A、P、B 三点共面且入射光线与反射光线位于法线两侧，反射光线、入射光线和法线必然位于同一平面内。由此，反射定律的两条内容（反射角等于入射角；三线共面）均从费马原理中自然导出。

4.4 由费马原理推导折射定律（斯涅尔定律）

折射定律的推导比反射情形稍复杂，因为光线穿越两种不同介质。设介质 1 的折射率为 n₁，介质 2 的折射率为 n₂，分界面为 x 轴。入射点 A 在介质 1 中，坐标为 (x₁, y₁)（y₁ > 0），出射点 B 在介质 2 中，坐标为 (x₂, y₂)（y₂ < 0）。光线在分界面上的入射点 P 的坐标为 (x, 0)。

此时光程为两段介质中的光程之和：

L(x) = n₁·|AP| + n₂·|PB| = n₁·√[(x - x₁)² + y₁²] + n₂·√[(x₂ - x)² + y₂²]

对 L(x) 求导并令其为零：

dL/dx = n₁·(x - x₁)/|AP| - n₂·(x₂ - x)/|PB| = 0

整理得：

n₁·(x - x₁)/|AP| = n₂·(x₂ - x)/|PB|

同样地，(x - x₁)/|AP| = sinθ₁（入射角的正弦），(x₂ - x)/|PB| = sinθ₂（折射角的正弦）。代入得：

n₁·sinθ₁ = n₂·sinθ₂

这就是著名的斯涅尔定律（Snell's Law），即折射定律。它表明：入射角与折射角的正弦之比等于两种介质折射率的反比。当 n₁ > n₂ 时，θ₂ > θ₁，光线远离法线；当入射角增大到临界角 θ_c = arcsin(n₂/n₁) 时，折射角达到 90°，继续增大入射角则发生全反射。

4.5 案例拆解：利用费马原理分析光纤中的光线传播

为了更好地理解费马原理的实际应用，我们以阶跃折射率多模光纤为例进行案例拆解。光纤纤芯折射率 n₁ = 1.46，包层折射率 n₂ = 1.44，纤芯半径 a = 25 μm，工作波长 λ₀ = 850 nm。

背景与目标：设计一个光纤耦合系统，需要确定光线在光纤中传播的最大入射角（数值孔径 NA），并分析不同入射角下光线的光程差异对模式色散的影响。

分析过程：

1. 全反射条件：根据折射定律，光线在纤芯-包层界面发生全反射的临界角 θ_c = arcsin(n₂/n₁) = arcsin(1.44/1.46) ≈ 80.6°（相对于法线）。因此，光线与光纤轴线的最大夹角（即最大入射角）为 α_max = 90° - θ_c ≈ 9.4°。
2. 数值孔径：NA = n₀·sinα_max = √(n₁² - n₂²) = √(1.46² - 1.44²) ≈ 0.24。这意味着光纤只能接收入射角小于约 14°（空气中）的光线。
3. 光程分析：沿轴线传播的光线（入射角为 0°）光程最短，L_min = n₁·L（L 为光纤长度）。以最大入射角传播的光线，其路径为锯齿形，实际几何路程更长，光程也更大。设光纤长度为 1 m，则最大光程差 ΔL ≈ n₁·L·(1/cosα_max - 1) ≈ 1.46 × 1 × (1/cos9.4° - 1) ≈ 0.0197 m。
4. 模式色散估算：光程差导致不同模式的光到达终点的时间差 Δt = ΔL/c ≈ 0.0197 / (3×10⁸) ≈ 65.7 ps。对于 850 nm 波长的多模光纤，这一色散限制了传输带宽。

操作步骤与结果：在实际耦合系统设计中，我们根据 NA = 0.24 选择合适的光纤耦合透镜，确保聚焦光锥的数值孔径小于光纤 NA。同时，通过光程差分析，我们预估了 1 km 光纤的脉冲展宽约为 65.7 ns，这决定了系统的最大传输速率约为 15 Mb/s（按脉冲展宽不超过码元周期的 1/4 估算）。

复盘与经验：费马原理在此案例中不仅帮助我们理解了全反射条件，更重要的是提供了光程差计算的数学工具。实际工程中，光程差分析是评估多模光纤带宽、设计波分复用系统以及优化光纤耦合效率的基础。若改用单模光纤（纤芯直径仅 8 μm），则模式色散几乎为零，但耦合难度显著增加，这正是设计中的权衡点。

4.6 常见误区与总结复盘

在学习费马原理及其应用时，学习者常陷入以下几个误区，需要特别警惕：

本节总结：费马原理以光程为核心概念，通过变分极值思想统一了反射定律与折射定律。光程 L = n·s 将几何路程与介质属性耦合在一起，其物理本质是光传播时间的度量。由费马原理推导反射和折射定律的过程，展示了从“原理”到“定律”的演绎逻辑，也为后续学习透镜成像、像差理论以及光学系统设计奠定了方法论基础。在实际工程中，光程差分析是光学设计、干涉测量和光纤通信中不可或缺的工具。掌握费马原理，意味着从“记住公式”跃升到“理解原理”，这是光学工程师从入门到精通的关键一步。