光学系统的矩阵表示法简介
在几何光学的传统分析中，追踪光线在复杂系统中的传播路径通常依赖于基于几何关系的代数计算。随着光学系统日益复杂，包含多个透镜、反射镜和间隔元件，传统的代数追踪法显得繁琐且易错。为了系统化、结构化地描述光线在光学系统中的变换过程，矩阵表示法（又称矩阵光学）被引入到几何光学领域。该方法通过线性代数工具，将光线的传播与折射行为抽象为矩阵乘法运算，为光学系统的总览设计与分析提供了极其高效且规范的数学框架。
矩阵光学的核心思想是将光线视为线性系统中的状态向量，而光学元件或传播空间则被视为对该向量进行变换的线性算子。在近轴近似的条件下，即光线与光轴的夹角足够小，使得正弦值可用弧度值线性替代时，光线的行为完全遵循线性规律。
在子午面（包含光轴的平面）内，一条光线在系统某一点的状态可以由两个参数唯一确定：
- 光线在该点距离光轴的高度，记为 $y$。
- 光线与光轴的夹角（或光线与光轴夹角的正切值），记为 $\theta$。
这两个参数构成了光线的状态列向量 $\begin{bmatrix} y \ \theta \end{bmatrix}$。当光线穿过一个光学元件或在空间中传播一段距离时，其状态发生线性改变。这种改变可以通过一个 $2 \times 2$ 的变换矩阵 $M$ 来描述。设入射光线状态为 $\begin{bmatrix} y_1 \ \theta_1 \end{bmatrix}$，出射光线状态为 $\begin{bmatrix} y_2 \ \theta_2 \end{bmatrix}$，则两者的关系可表示为：
$$
\begin{bmatrix} y_2 \ \theta_2 \end{bmatrix} = M \begin{bmatrix} y_1 \ \theta_1 \end{bmatrix} = \begin{bmatrix} A & B \ C & D \end{bmatrix} \begin{bmatrix} y_1 \ \theta_1 \end{bmatrix}
$$
其中，矩阵 $M = \begin{bmatrix} A & B \ C & D \end{bmatrix}$ 被称为该光学元件的射线传递矩阵（Ray Transfer Matrix）或 ABCD 矩阵。矩阵的四个元素具有明确的物理量纲和物理意义：$A$ 和 $D$ 是无量纲的，分别表示高度的线性放大率和角放大率；$B$ 具有长度量纲，表示光线高度对入射角的依赖；$C$ 具有长度倒数量纲，表示出射角对入射高度的依赖，与元件的光焦度直接相关。
光学系统变换的横向对比
矩阵表示法的最大优势在于其能够以统一的数学形式横向对比不同物理机制对光线的作用。在几何光学中，光线传播的基本规律主要包含直线传播、折射与反射。在矩阵框架下，这些规律被高度抽象为不同特征的 ABCD 矩阵。
自由空间传播：光线在均匀介质中沿直线传播。经过长度为 $d$ 的传播，高度增加 $d \cdot \theta$，而角度保持不变。其变换矩阵为：
$$
M_{space} = \begin{bmatrix} 1 & d \ 0 & 1 \end{bmatrix}
$$
该矩阵的特征是下三角元素 $C=0$，表明自由空间不具备改变光线角度（即聚焦或发散）的能力。折射界面：光线穿过曲率半径为 $R$ 的球面折射界面。根据近轴折射定律，高度在折射瞬间不发生改变，而角度的改变与界面的光焦度相关。其变换矩阵为：
$$
M_{refract} = \begin{bmatrix} 1 & 0 \ -\frac{n_2-n_1}{n_2 R} & \frac{n_1}{n_2} \end{bmatrix}
$$
其中 $n_1$ 和 $n_2$ 分别为入射侧和出射侧的折射率。该矩阵的特征是上三角元素 $B=0$，表明折射发生在同一空间位置上。反射界面：光线在球面反射镜上发生反射。在矩阵光学中，反射被视为折射的一种特例，即令 $n_2 = -n_1$。其变换矩阵为：
$$
M_{reflect} = \begin{bmatrix} 1 & 0 \ -\frac{2}{R} & 1 \end{bmatrix}
$$
通过这种统一的矩阵形式，设计者可以在同一计算流程中无缝交替处理包含透镜（折射）和反射镜的复杂折反系统。
复杂系统的级联与矩阵乘法
实际的光学仪器通常由多个元件和间隔组成。矩阵表示法处理级联系统时遵循线性代数中的矩阵乘法法则：系统总传递矩阵等于各个元件传递矩阵按逆序相乘。
假设一个光学系统由元件 1、元件 2 和元件 3 依次组成，其各自的传递矩阵为 $M_1$、$M_2$、$M_3$，光线首先经过元件 1，最后经过元件 3。则整个系统的总矩阵 $M_{total}$ 为：
$$
M_{total} = M_3 \cdot M_2 \cdot M_1
$$
这种级联特性使得极其复杂的光学系统可以被压缩为一个单一的 $2 \times 2$ 矩阵。一旦获得了系统的总矩阵，即可立即判断系统的宏观成像属性，而无需逐步追踪中间光线的轨迹。例如，若系统总矩阵的 $C$ 元素为零，则该系统为望远系统；若 $B$ 元素为零，则表明输入面与输出面构成了一对物像共轭面。
应用全景与系统级评估
在光学系统的全景设计与评估中，矩阵表示法不仅用于光线追踪，更是推导系统级通用原理的基石。以下是其在系统级应用的几个典型方面：
- 系统等效焦距计算：对于由多个透镜组成的复合系统，通过计算总矩阵的 $C$ 元素，可以直接求出系统的等效焦距和后焦距，从而评估系统的聚光能力。
- 高斯光束变换：虽然矩阵表示法起源于几何光学，但其数学结构同样适用于波动光学中的近轴高斯光束传播。利用 ABCD 定律，可以快速计算激光光束在复杂谐振腔或传输系统中的束腰位置和光斑尺寸。
- 光学谐振腔稳定性判据：在激光技术中，谐振腔由两个反射镜和中间的自由空间组成。将一个往返周期的变换矩阵求出后，根据矩阵的迹（$A+D$）的绝对值是否小于 2，即可作为谐振腔稳定性的普适判据。这是矩阵光学在系统级评估中最经典的应用之一。
综上所述，光学系统的矩阵表示法通过高度抽象的线性代数语言，将几何光学中光的传播、折射与反射规律统一在 ABCD 矩阵框架下。它避免了繁杂的中间几何计算，使设计者能够站在系统总览的高度，快速评估和优化复杂光学仪器的宏观特性，是现代光学工程不可或缺的基础工具。