【光学】利用Matlab仿真拉盖尔-高斯光束的轨道角动量特性及其应用探索

1. 从“甜甜圈”光斑说起：什么是拉盖尔-高斯光束？

如果你玩过激光笔，或者看过一些科幻电影里的激光武器，你印象中的光束可能就是一个明亮的光点。但在光学研究的世界里，光束的形状可以千变万化，其中有一种光束，它打出来的光斑不是实心的圆点，而是一个亮环，中心反而是暗的，看起来就像一个发光的“甜甜圈”。这种光束，就是拉盖尔-高斯光束，我们通常简称它为LG光束。

我第一次在实验室里用空间光调制器生成出这种“甜甜圈”光斑时，感觉特别神奇。这不仅仅是形状上的不同，它背后隐藏着一个更酷的特性：轨道角动量。你可以把它想象成光在“旋转”着前进。普通的光束，比如我们常见的激光，它的角动量主要来自光波的偏振方向（自旋角动量），就像地球的自转。而LG光束携带的轨道角动量，则更像是地球绕着太阳的公转，是整个光波前在空间上绕着传播轴做螺旋运动。

这种螺旋运动，体现在它的相位分布上。如果你能“看到”光的相位，你会发现LG光束的相位图就像一个螺旋的楼梯，或者一个开瓶器的螺纹，我们称之为螺旋相位。这个螺旋的“圈数”，在物理上用拓扑荷数 l 来表示，它决定了光束轨道角动量的大小。l可以是0，±1，±2……当l=0时，它就是普通的高斯光束，没有甜甜圈，也没有螺旋相位。当l不为零时，中心就会出现相位奇点，光强为零，从而形成了那个标志性的中空暗斑。

那么，我们为什么要研究这种看起来有点“奇怪”的光束呢？因为它太有用了。这种独特的相位结构和携带的轨道角动量，让LG光束在光通信、光学微操控、超分辨成像等领域大放异彩。比如在光通信里，不同拓扑荷数的LG光束可以当作独立的信道来同时传输数据，理论上能极大提升通信容量，这就是轨道角动量复用技术。在光学操控里，用这种光束照射微小的颗粒，可以产生特殊的力学效应，像光镊一样“抓住”并旋转粒子。所以，无论是做理论研究还是工程应用，掌握LG光束的仿真和生成都是非常关键的一步。

对于刚接触这个领域的朋友，或者想用Matlab做光学仿真的同学，理解并可视化LG光束的光强和相位，就是叩开这扇大门的第一个台阶。接下来，我就带你一步步用Matlab把它“造”出来，并看看我们能用它做些什么有意思的探索。

2. 手把手用Matlab构建你的第一个LG光束

理论说再多，不如动手跑一遍代码来得实在。用Matlab仿真LG光束的好处是，我们不需要昂贵的实验设备，在电脑上就能直观地看到各种参数下的光束形态，这对于理解和设计光学系统至关重要。下面，我就把我常用的仿真步骤拆开揉碎了讲给你听。

2.1 仿真环境与核心思路搭建

首先，我们得在Matlab里创建一个用于计算的“舞台”，也就是我们的仿真平面。这个平面是一个二维的网格，覆盖我们关心的空间区域。通常，我们会使用meshgrid函数来生成这个网格。

%% 1. 定义仿真区域和网格
N = 512; % 网格点数，决定了图像分辨率，512或1024都是常用值
delta = 0.003; % 网格间距（米），决定了仿真区域的实际物理大小
% 生成从 -N/2*delta 到 (N/2-1)*delta 的坐标轴
x = (-N/2 : N/2-1) * delta;
y = (-N/2 : N/2-1) * delta;
% 创建二维网格坐标
[X, Y] = meshgrid(x, y);

这里有个小经验，N取2的整数次幂（如256, 512, 1024）有时在后续做傅里叶变换时会比较高效。delta的大小需要根据你关心的实际尺度来定，比如模拟实验室中毫米尺度的光斑，delta可以设为几十到几百微米。

有了直角坐标(X, Y)，我们通常还需要极坐标(ρ, θ)，因为LG光束的数学表达式在极坐标下更简洁美观。用Matlab自带的cart2pol函数可以轻松转换：

%% 2. 转换为极坐标
[Theta, Rho] = cart2pol(X, Y); % Theta是方位角，Rho是径向距离

现在，舞台搭好了。LG光束的复振幅表达式是它的灵魂，其核心包含三部分：高斯包络、拉盖尔多项式、以及最重要的螺旋相位项。

%% 3. 定义LG光束参数并计算复振幅
w0 = 0.3; % 束腰半径（米），控制光束的粗细
lambda = 632.8e-9; % 光波长，例如常用的氦氖激光波长632.8纳米
k = 2 * pi / lambda; % 波数
z = 0; % 观测面距离光束腰部的距离，这里先看束腰处的面
l = 4; % 拓扑荷数，决定相位螺旋的“圈数”，也是轨道角动量的阶数
p = 0; % 径向指数，通常从0开始，决定径向亮环的数目

% 计算在距离z处的光束半径w(z)和曲率半径R(z)
wz = w0 * sqrt(