4. 自由空间光束传输建模
在激光光学系统设计中，我们经常需要回答一个核心问题：一束已知参数的高斯光束，经过若干自由空间距离、透镜、反射镜或介质界面后，在目标位置的光斑半径、波前曲率半径和束腰位置分别是多少？如果每次都用菲涅耳衍射积分逐点计算，工程上既繁琐又难以获得解析洞察。因此，本节围绕一条从波动光学到矩阵光学的建模主线展开：近轴波动方程给出光束传播的基本微分方程，菲涅耳衍射提供标量积分描述，ABCD定律与传输矩阵则把复杂光学系统压缩为2×2矩阵，使高斯光束参数可以像几何光线一样被高效追踪。
一、近轴波动方程与高斯光束解
标量波动方程在缓变振幅近似下可化为近轴波动方程。设光场沿z轴传播，复振幅为U(x,y,z)，则近轴近似为：
∂²U/∂x² + ∂²U/∂y² - j2k ∂U/∂z = 0，其中 k = 2π/λ该方程最基本的高斯光束解可写为：
U(r,z) = (w₀/w(z)) · exp(-r²/w²(z)) · exp[-j(kz + kr²/(2R(z)) - ζ(z))]其中w(z)为光斑半径，R(z)为波前曲率半径，ζ(z)为古伊相移。三个参数随z的变化为：
w(z) = w₀ · sqrt(1 + (z/z_R)²)
R(z) = z · [1 + (z_R/z)²]
z_R = πw₀²/λ这里z_R称为瑞利长度，是判断光束处于准直区还是发散区的关键尺度。当|z|≤z_R时，光斑半径变化不超过束腰的√2倍，工程上常把这一段视为准直工作区。例如λ=1064 nm、w₀=0.5 mm时，z_R≈0.738 m；若目标工作距离为3 m，则z/z_R≈4.07，光斑将明显扩展，必须重新评估聚焦或准直方案。
二、菲涅耳衍射：从积分到工程判据
菲涅耳衍射积分描述光场从孔径面到观察面的传播：
U(x,y,z) = (e^{jkz}/(jλz)) · ∬ U₀(x₀,y₀) · exp{jk[(x-x₀)²+(y-y₀)²]/(2z)} dx₀dy₀工程上更常用菲涅耳数N_F来判断衍射区：N_F = a²/(λz)，a为孔径半径。N_F≫1时可用几何光学近似；N_F≈1时进入菲涅耳衍射区；N_F≪1时进入夫琅禾费远场区。对于高斯光束，由于其场分布本身是近轴方程的本征解，自由空间传播不会改变高斯形状，只改变w、R和ζ，这正是ABCD定律能够精确描述高斯光束的前提。
三、ABCD定律与传输矩阵
近轴光学系统中，光线矢量(r, θ)经过光学元件后线性变换：
[r₂] [A B] [r₁]
[θ₂] = [C D] [θ₁]对高斯光束，复光束参数q(z)定义为：
1/q(z) = 1/R(z) - jλ/(πw²(z))ABCD定律给出：
q₂ = (A q₁ + B)/(C q₁ + D)自由空间距离L的矩阵为[[1, L],[0, 1]]；薄透镜焦距f的矩阵为[[1, 0],[-1/f, 1]]。多个元件级联时，总矩阵按从输出面到输入面的顺序依次左乘。例如自由空间L₁、薄透镜f、自由空间L₂的系统矩阵为：
M = [[1, L₂],[0, 1]] · [[1, 0],[-1/f, 1]] · [[1, L₁],[0, 1]]
= [[1-L₂/f, L₁+L₂-L₁L₂/f],[-1/f, 1-L₁/f]]得到M后，将入射q₁代入q₂=(Aq₁+B)/(Cq₁+D)，再分离实部与虚部，即可求出输出面的R₂和w₂。若需要求束腰位置，可令1/q的实部为零，反解出z。
四、透镜变换与工程应用要点
薄透镜对高斯光束的作用可概括为：透镜不改变光斑半径w，但会改变波前曲率半径。若入射到透镜处的波前曲率半径为R₁，则出射曲率半径R₂满足1/R₂=1/R₁-1/f。当R₁→∞（入射为准直光）时，R₂=-f，光束被聚焦。聚焦光斑半径的经典公式为：
w_f ≈ λf/(πw_L)其中w_L为透镜处的入射光斑半径。该式说明：要获得更小的聚焦光斑，应增大入射光束在透镜上的光斑尺寸，或缩短焦距。但w_L不能超过透镜有效孔径，否则会引入截断衍射，实际聚焦光斑反而变大。工程上通常取w_L≈0.8×透镜半径作为折中。
五、建模流程与结果判读
实际建模时，建议按以下顺序执行：第一，明确入射高斯光束的w₀、z₀和λ，写出q₁；第二，按光路顺序列出每个元件的ABCD矩阵，注意矩阵乘法顺序；第三，计算总矩阵M；第四，用q₂=(Aq₁+B)/(Cq₁+D)求输出q₂；第五，从1/q₂分离出w₂和R₂；第六，若关注束腰，令Re(1/q)=0反解位置。结果判读时重点看三个量：w₂是否满足目标光斑要求，R₂是否接近无穷大（准直）或等于目标曲率，以及z_R是否覆盖工作距离。
下面用一张信号链与参数关系图，把“入射q参数→传输矩阵→输出q参数→可观测光束参数”的建模链路和关键公式集中呈现，便于在工程计算时快速对照。
该图把建模链路分为输入、矩阵、输出三段，并在下方集中列出瑞利长度、聚焦光斑、菲涅耳数和透镜变换四个关键判据。实际使用时，先由λ和w₀算出z_R，再根据光路写出总矩阵M，代入q₁得到q₂，最后从1/q₂分离出w₂和R₂。若w₂与目标不符，可调整透镜焦距f或透镜位置L₁；若焦深不足，则需要减小w_f或改用更长焦距配合更大w_L。整个过程中，ABCD定律把复杂的衍射积分转化为矩阵代数，是自由空间光束传输建模最实用的工程工具。