3. 几何光学与矩阵光学基础：光线传播矩阵（ABCD矩阵）、透镜与反射镜的矩阵表示、谐振腔的往返矩阵构建
在激光谐振腔的分析与设计中，几何光学与矩阵光学提供了最直观且高效的数学工具。当谐振腔的横向尺寸远大于工作波长，且我们主要关注光线在腔内的宏观传播轨迹（而非衍射极限下的场分布）时，可以将激光束近似为几何光线，用近轴光线传播的ABCD矩阵来描述其通过自由空间、透镜、反射镜等光学元件的变换规律。本章将系统建立光线传播矩阵（ABCD矩阵）的数学框架，详细推导透镜与反射镜的矩阵表示，并重点讲解如何构建谐振腔的往返矩阵，为后续分析腔的稳定性、模式分布及高斯光束传播奠定坚实的数学基础。
上图展示了本章的核心知识框架：从光线矢量与传播矩阵的基本定义出发，建立自由空间、透镜、反射镜等基本元件的矩阵表示，最终通过矩阵连乘构建谐振腔的往返矩阵，并利用其矩阵元判断腔的稳定性。下面我们将逐一深入展开。
3.1 光线传播矩阵（ABCD矩阵）的核心原理
在近轴光学中，一条光线在某一横截面上的状态可以用两个参数唯一确定：光线离光轴的垂直距离 \( r \) 和光线与光轴的夹角 \( \theta \)（近轴近似下取 \(\sin\theta \approx \theta\)）。这两个参数构成一个二维列向量 \(\begin{pmatrix} r \\ \theta \end{pmatrix}\)，称为光线矢量。当光线通过一个光学系统（如一段自由空间、一个透镜或一面反射镜）后，出射光线矢量与入射光线矢量之间满足线性变换关系：
\[ \begin{pmatrix} r_2 \\ \theta_2 \end{pmatrix} = \begin{pmatrix} A & B \\ C & D \end{pmatrix} \begin{pmatrix} r_1 \\ \theta_1 \end{pmatrix} \]
其中 \(2 \times 2\) 矩阵 \(\begin{pmatrix} A & B \\ C & D \end{pmatrix}\) 即为该光学系统的光线传播矩阵，通常称为ABCD矩阵。矩阵元具有明确的物理意义：A 表示入射高度对出射高度的影响（角放大率的倒数），B 表示入射角对出射高度的贡献（等效光程），C 表示入射高度对出射角的贡献（光焦度），D 表示入射角对出射角的贡献（角放大率）。对于无源无损耗光学系统，矩阵行列式满足 \(AD - BC = 1\)，这是能量守恒在几何光学中的体现。
ABCD矩阵的最大优势在于其级联性：若光线依次通过矩阵为 \(M_1, M_2, \ldots, M_n\) 的多个元件，则整个系统的等效矩阵为 \(M = M_n \cdots M_2 \cdot M_1\)。注意矩阵乘法的顺序必须与光线传播的先后顺序一致，即后通过的元件矩阵写在左边。这一性质使得复杂光学系统的分析简化为矩阵连乘运算，非常适合计算机编程实现。
3.2 透镜与反射镜的矩阵表示
在谐振腔中，最常见的两个光学元件是薄透镜和球面反射镜。它们的ABCD矩阵形式简洁且物理意义清晰。
（1）自由空间的传播矩阵
光线在折射率为 \(n\) 的均匀介质中传播距离 \(L\) 时，高度增加 \(L \cdot \theta\)，而角度保持不变。因此自由空间的传播矩阵为：
\[ T(L) = \begin{pmatrix} 1 & L \\ 0 & 1 \end{pmatrix} \]
若介质折射率为 \(n\)，则需将 \(L\) 替换为光学距离 \(L/n\)。在谐振腔中，通常默认腔内介质为空气（\(n \approx 1\)），直接使用几何距离。
（2）薄透镜的矩阵
对于焦距为 \(f\) 的薄透镜，光线通过时高度不变，但角度发生偏转，偏转量与离轴高度成正比：\(\theta_2 = \theta_1 - r_1/f\)。因此薄透镜矩阵为：
\[ L(f) = \begin{pmatrix} 1 & 0 \\ -1/f & 1 \end{pmatrix} \]
焦距 \(f\) 为正表示会聚透镜，为负表示发散透镜。在谐振腔中，有时用透镜等效替代热透镜效应或腔内其他聚焦元件。
（3）球面反射镜的矩阵
球面反射镜对光线的变换与薄透镜类似，但需注意反射引起的光线方向反转。对于曲率半径为 \(R\) 的球面反射镜（凹面镜 \(R>0\)，凸面镜 \(R<0\)），在近轴条件下，其反射矩阵为：
\[ M_{mirror} = \begin{pmatrix} 1 & 0 \\ -2/R & 1 \end{pmatrix} \]
这里 \(2/R\) 是反射镜的光焦度。对比薄透镜矩阵可以发现，球面反射镜等效于一个焦距 \(f = R/2\) 的薄透镜。这一等效关系在分析谐振腔时非常有用，可以将反射镜“展开”为透镜来处理。
（4）平面镜的矩阵
平面镜（\(R \to \infty\)）的矩阵为单位矩阵 \(\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}\)，即对光线不产生任何变换，仅改变传播方向。在构建往返矩阵时，平面镜通常只作为参考面，不引入额外的矩阵元。
| 光学元件 | 符号 | ABCD矩阵 | 物理参数 |
|---|---|---|---|
| 自由空间传播 | \(T(L)\) | \(\begin{pmatrix} 1 & L \\ 0 & 1 \end{pmatrix}\) | \(L\): 传播距离 |
| 薄透镜 | \(L(f)\) | \(\begin{pmatrix} 1 & 0 \\ -1/f & 1 \end{pmatrix}\) | \(f\): 焦距 |
| 球面反射镜 | \(M(R)\) | \(\begin{pmatrix} 1 & 0 \\ -2/R & 1 \end{pmatrix}\) | \(R\): 曲率半径 |
| 平面镜 | \(M(\infty)\) | \(\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}\) | 单位矩阵 |
3.3 谐振腔的往返矩阵构建
谐振腔的往返矩阵描述了光线在腔内完成一次完整往返（从某一参考面出发，经过所有元件，回到同一参考面且方向相同）后的变换关系。构建往返矩阵是分析腔稳定性和模式分布的关键步骤。
构建步骤：
- 选择参考面：通常选在某一反射镜表面（如M₁），作为往返的起点和终点。
- 确定往返路径：从参考面出发，依次经过腔内所有元件，到达另一反射镜（M₂），反射后沿原路返回，再次经过所有元件，回到参考面。注意反射镜的矩阵只计算一次（在反射点处）。
- 按顺序连乘矩阵：将往返路径上所有元件的矩阵按光线传播的先后顺序从左到右排列并相乘。例如，对于两镜腔（M₁和M₂，腔长L），从M₁表面出发，先经过自由空间L到达M₂，经M₂反射，再经过自由空间L回到M₁，最后经M₁反射（但M₁的反射矩阵通常不包含在往返矩阵中，因为参考面选在M₁表面，往返矩阵描述的是从M₁表面出发并回到M₁表面的变换，M₁的反射效应通过边界条件处理）。因此往返矩阵为：
\[ M = T(L) \cdot M_{mirror2} \cdot T(L) = \begin{pmatrix} 1 & L \\ 0 & 1 \end{pmatrix} \begin{pmatrix} 1 & 0 \\ -2/R_2 & 1 \end{pmatrix} \begin{pmatrix} 1 & L \\ 0 & 1 \end{pmatrix} \]
计算得到：
\[ M = \begin{pmatrix} 1 - \frac{2L}{R_2} & 2L - \frac{2L^2}{R_2} \\ -\frac{2}{R_2} & 1 - \frac{2L}{R_2} \end{pmatrix} \]
若参考面选在M₂表面，则往返矩阵为 \(M' = T(L) \cdot M_{mirror1} \cdot T(L)\)，形式类似。两种选择得到的矩阵虽然不同，但矩阵的迹（\(A+D\)）相同，因此稳定性判断结果一致。
稳定性条件：谐振腔的稳定性要求光线在腔内往返任意多次后不逸出腔外，即往返矩阵的幂次 \(M^n\) 保持有界。数学上等价于矩阵的迹满足：
\[ \left| \frac{A+D}{2} \right| \le 1 \]
定义谐振腔的 \(g\) 参数：\(g_1 = 1 - \frac{L}{R_1}\)，\(g_2 = 1 - \frac{L}{R_2}\)，则稳定性条件可改写为熟悉的 \(0 \le g_1 g_2 \le 1\)。这一判据在腔设计中极为常用。
3.4 案例拆解：典型两镜腔的往返矩阵与稳定性分析
我们以一个典型的平凹腔为例，演示往返矩阵的构建和稳定性判断的完整流程。
案例背景：某固体激光器采用平凹腔设计，输出镜为平面镜（M₁，\(R_1 = \infty\)），后腔镜为凹面镜（M₂，\(R_2 = 500 \text{ mm}\)），腔长 \(L = 300 \text{ mm}\)。需要验证该腔是否稳定，并计算往返矩阵的具体形式。
分析过程：
第一步，确定参考面。选择在平面镜M₁表面作为参考面。
第二步，写出往返路径上的元件序列：从M₁出发 → 自由空间传播L → M₂反射 → 自由空间传播L → 回到M₁。由于M₁是平面镜，其矩阵为单位阵，不改变光线状态，因此往返矩阵为：
\[ M = T(L) \cdot M_{mirror2} \cdot T(L) \]
代入 \(L = 300 \text{ mm}\)，\(R_2 = 500 \text{ mm}\)：
\[ M = \begin{pmatrix} 1 & 300 \\ 0 & 1 \end{pmatrix} \begin{pmatrix} 1 & 0 \\ -2/500 & 1 \end{pmatrix} \begin{pmatrix} 1 & 300 \\ 0 & 1 \end{pmatrix} \]
先计算后两个矩阵的乘积：
\[ \begin{pmatrix} 1 & 0 \\ -0.004 & 1 \end{pmatrix} \begin{pmatrix} 1 & 300 \\ 0 & 1 \end{pmatrix} = \begin{pmatrix} 1 & 300 \\ -0.004 & -0.2 \end{pmatrix} \]
再与第一个矩阵相乘：
\[ M = \begin{pmatrix} 1 & 300 \\ 0 & 1 \end{pmatrix} \begin{pmatrix} 1 & 300 \\ -0.004 & -0.2 \end{pmatrix} = \begin{pmatrix} 1 - 1.2 & 300 - 60 \\ -0.004 & -0.2 \end{pmatrix} = \begin{pmatrix} -0.2 & 240 \\ -0.004 & -0.2 \end{pmatrix} \]
第三步，稳定性判断。计算矩阵的迹：\(A + D = -0.2 + (-0.2) = -0.4\)，则 \((A+D)/2 = -0.2\)，其值 \(0.2 \le 1\)，满足稳定性条件。同时用 \(g\) 参数验证：\(g_1 = 1 - L/R_1 = 1 - 0 = 1\)，\(g_2 = 1 - L/R_2 = 1 - 300/500 = 0.4\)，\(g_1 g_2 = 0.4\)，在0到1之间，腔稳定。
第四步，结果分析。往返矩阵的B元为240 mm，表示光线在往返一次后，出射高度对入射角的敏感程度。C元为-0.004 mm⁻¹，表示光焦度。这些参数可用于后续计算腔的本征模式（高斯光束的q参数变换）。
复盘与扩展：若将腔长增加到 \(L = 600 \text{ mm}\)，则 \(g_2 = 1 - 600/500 = -0.2\)，\(g_1 g_2 = -0.2\)，不满足稳定性条件，腔变为非稳腔。这说明同一对反射镜，腔长变化可能导致稳定性突变。在实际激光器设计中，必须通过往返矩阵或g参数预先验证工作点位于稳定区内。
3.5 常见误区与总结复盘
在学习和应用ABCD矩阵分析谐振腔时，以下误区需要特别警惕：
- 矩阵顺序错误：光线先经过的元件矩阵写在右边，后经过的写在左边。很多初学者习惯从左到右写，导致结果完全错误。
- 反射镜符号混淆：凹面镜的曲率半径 \(R\) 在矩阵中取正值，凸面镜取负值。但若反射镜的凹面朝向腔外，则符号需反转。务必结合光路图确认。
- 忽略参考面选择的影响：往返矩阵的具体形式依赖于参考面的位置，但稳定性判据（迹的一半）与参考面无关。若计算出的稳定性结果随参考面变化，说明矩阵构建有误。
- 将往返矩阵与单程矩阵混淆：单程矩阵只描述光线从一镜到另一镜的传播，往返矩阵则包含完整的往返过程。在计算腔的损耗或模式时，必须使用往返矩阵。
- 近轴近似的适用范围：ABCD矩阵仅适用于近轴光线（\(\theta\) 很小）。对于大角度光线或大离轴距离，几何光学本身失效，需采用波动光学或光线追迹软件。
总结本章，我们建立了光线传播矩阵的完整理论框架：从光线矢量定义出发，给出了自由空间、薄透镜、球面反射镜的ABCD矩阵，并详细演示了谐振腔往返矩阵的构建流程和稳定性判据。这些内容不仅是后续分析腔模式、高斯光束传输的基础，也是激光器设计中最常用的工程工具。熟练掌握矩阵连乘和g参数分析，是每一位激光工程师的基本功。