4、ABCD矩阵法：光线传输矩阵、谐振腔稳定性分析、本征模式求解

好，咱们今天聊点硬核的——ABCD矩阵法。

说实话，我刚入行那会儿，看到一堆2×2矩阵就头疼。心想，这玩意儿不就是把几何光学那套搬过来吗？后来真做项目了才发现，嗯，真香。尤其是分析谐振腔的时候，没有矩阵法，你根本没法系统性地判断腔体稳不稳、模式长啥样。

4.1 光线传输矩阵：从几何光学到矩阵语言

先说说基本概念。ABCD矩阵，也叫光线传输矩阵，描述的是光线在光学系统中的传播行为。

一条光线，我们用两个参数描述：

• r：光线离光轴的距离
• θ：光线与光轴的夹角（小角度近似下，θ ≈ tanθ）

光线通过一个光学元件后，输出参数和输入参数的关系可以写成：

[ r2 ]   =   [ A  B ] [ r1 ]
[ θ2 ]       [ C  D ] [ θ1 ]

说白了，就是用一个2×2矩阵，把光线从“入口”映射到“出口”。

我个人习惯把常见元件的矩阵背下来，因为项目里翻书太慢了。这里列几个最常用的：

光学元件	ABCD矩阵	说明
自由空间传播（距离L）	[1 L; 0 1]	光线走一段直线
薄透镜（焦距f）	[1 0; -1/f 1]	只改变角度，不改变位置
球面反射镜（曲率半径R）	[1 0; -2/R 1]	和透镜类似，但注意符号
介质界面（折射率n1→n2）	[1 0; 0 n1/n2]	平面界面，只改变角度

4.2 谐振腔稳定性分析：一个矩阵搞定一切

谐振腔稳不稳定，说白了就是光线在腔内来回反射会不会跑出去。用ABCD矩阵分析，这事儿变得特别简单。

对于一个由两个球面镜（曲率半径R1、R2，间距L）组成的谐振腔，光线往返一周的传输矩阵可以算出来。我直接给结论：

M = [ A  B ]
    [ C  D ]

其中：
A = 2g1g2 - 1
B = 2g2L
C = 2g1(g2-1)/L
D = 2g1g2 - 1

g1 = 1 - L/R1
g2 = 1 - L/R2

稳定条件是什么？很简单：

为什么会这样？因为矩阵的特征值决定了光线轨迹的发散程度。如果特征值的模大于1，光线每走一圈就离轴更远，最终跑出腔外——这就是不稳定腔。

我记得有一次做固体激光器设计，客户要求腔长特别短。我按经验选了个平平腔（两个平面镜），结果一算g1·g2=1，刚好在边界上。实际装调时稍微有点对准误差，激光就出不来。后来换成平凹腔，g1·g2=0.5，稳得很。

4.3 本征模式求解：从矩阵到光斑大小

稳定腔不光要稳，还得知道腔内光斑长什么样。这就是本征模式求解干的事。

对于基模（TEM00），光斑半径ω和波前曲率半径R可以通过ABCD矩阵的元素求出来。这里有个关键公式：

q参数定义：
1/q = 1/R - iλ/(πω²)

经过ABCD矩阵变换：
q2 = (A·q1 + B) / (C·q1 + D)

在谐振腔中，自洽条件要求光线往返一周后q参数不变：

q = (A·q + B) / (C·q + D)

解这个方程，得到：

1/q = (D - A) / (2B) ± i·√[1 - ((A+D)/2)²] / B

实部对应波前曲率，虚部对应光斑大小。你想想看，只要知道腔的ABCD矩阵，光斑半径直接算出来，多方便。

我实际项目中常用这个公式来估算腔内的光斑尺寸，尤其是设计端面泵浦固体激光器时，泵浦光斑和振荡光斑要匹配，不然效率上不去。

4.4 知识体系总览

下面这张图是我自己整理的ABCD矩阵法知识框架，帮你理清思路：

嗯，这张图把ABCD矩阵法的三个核心模块串起来了。从光线传输矩阵出发，往左走是稳定性分析，往右走是本征模式求解，最终都落到实际腔型设计上。

4.5 实战小例子

最后给个具体算例。假设一个平凹腔：平面镜R1=∞，凹面镜R2=500mm，腔长L=300mm。

g1 = 1 - 300/∞ = 1
g2 = 1 - 300/500 = 0.4
g1·g2 = 0.4 → 在0~1之间，稳定！

往返矩阵：
A = 2×1×0.4 - 1 = -0.2
B = 2×0.4×300 = 240
C = 2×1×(0.4-1)/300 = -0.004
D = -0.2

(A+D)/2 = -0.2 → 绝对值0.2 ≤ 1，确认稳定

然后算q参数，就能得到凹面镜上的光斑半径。具体数值我就不展开了，留给你自己练练手。

好了，ABCD矩阵法就聊到这儿。这个方法贯穿了整个激光器光路设计，从稳定性到模式计算，一套矩阵全搞定。你想想看，是不是比追迹几百条光线要优雅得多？

---