Matlab实战：5分钟搞定OAM涡旋光束仿真（附完整源码）

从零到一：用Matlab亲手“点亮”你的第一束涡旋光

如果你正在光学实验室里，或者刚刚开始接触轨道角动量（OAM）这个听起来有点科幻的概念，那么“仿真”可能是你绕不开的第一步。毕竟，搭建一套能产生完美拉盖尔-高斯光束的光路，不仅需要精密的仪器，还需要反复调试的耐心。而Matlab，就像是我们手边的一个虚拟光学实验台，它让我们能在代码的世界里，自由地构建、操控和观察这些拥有奇特相位结构的光束。今天，我们不谈那些复杂的数学推导，也不去深究物理原理的每一个细节，我们就来干一件最实在的事：用Matlab，从无到有，亲手生成一束携带轨道角动量的涡旋光，并把它清晰地“画”出来。 这个过程，对于理解OAM光束的形态、探索不同模式叠加的效果，甚至是为后续的物理实验做预演，都至关重要。

1. 理解核心：涡旋光束到底是什么？

在动手写代码之前，我们得先搞清楚我们要“造”的是个什么东西。你可以把一束普通的高斯激光想象成一根光滑、均匀的光柱，它的波前（相位面）是一系列平行的平面。而涡旋光束，特别是拉盖尔-高斯（LG）光束，它的波前被“拧”成了螺旋状。想象一下DNA的双螺旋结构，或者一个旋转的楼梯，光场的相位沿着环绕光束中心的方位角方向连续变化。

这种螺旋相位结构带来了两个最直观的特征：

1. 中心暗核：在光束的中心点，所有相位的螺旋线汇聚于此，导致相位不确定（奇点），因此该点的光强为零。所以在光强分布图上，你会看到一个明亮的环，而中心是暗的。
2. 轨道角动量（OAM）：这是涡旋光束最核心的物理属性。不同于与光偏振相关的自旋角动量，OAM源于光场在空间上的螺旋相位分布。其“强度”由一个整数 拓扑荷数 l 来描述。l 可以是正数、负数或零（零就是普通高斯光束）。l 的绝对值越大，相位螺旋得越“紧”，理论上携带的OAM也越大。l 的正负则决定了螺旋的方向（左旋还是右旋）。

除了拓扑荷 l，描述一个LG光束还需要另一个整数 径向指数 p。p 决定了光束在径向（从中心向外）上有多少个光强环。p=0 时，只有一个亮环；p=1 时，有两个亮环，以此类推。

注意：我们常说的“OAM模式”，通常就是指一组特定的 (l, p) 值。在通信领域，不同的 l 值可以被用作独立的数据信道，这也就是OAM复用技术的基础。

所以，我们的Matlab仿真目标就很明确了：编写一个函数，输入我们想要的 l 和 p，以及一些光束的基本参数（如波长、束腰半径），它就能输出一个二维复数矩阵，这个矩阵精确代表了该LG光束在某个横截面上的复振幅分布（包含振幅和相位信息）。

2. 搭建你的Matlab光学仿真环境

工欲善其事，必先利其器。我们不需要任何特殊的工具箱，基础的Matlab环境就足够了。但为了代码清晰、易于管理和复用，我强烈建议你按照以下结构来组织你的项目文件夹。这是我踩过几次坑后总结出来的高效做法。

在你的工作目录下，创建如下几个文件：

你的项目文件夹/
├── main.m                 % 主脚本，用于控制整体流程和绘图
├── generateLGBeam.m       % 核心函数：生成单模LG光束
├── propagateBeam.m        % （可选）函数：模拟光束传播
├── plotBeamProfile.m      % 函数：绘制光强和相位图
└── utils/                 % 子文件夹，存放辅助函数
    ├── makeGrid.m         % 生成计算网格
    └── normalize.m        % 数据归一化

首先，我们从最基础的网格生成开始。所有的光束计算都需要在一个离散的二维网格上进行。utils/makeGrid.m 这个辅助函数会让后续的代码简洁很多。

% utils/makeGrid.m
function [X, Y, x, y] = makeGrid(N, L)
% 生成二维计算网格
% 输入：
%   N - 网格点数（单边），通常取256, 512, 1024等2的幂次
%   L - 物理尺寸（单边长度，单位自定，如米或毫米）
% 输出：
%   X, Y - 二维网格矩阵 (N x N)
%   x, y - 一维坐标向量 (1 x N)

    dx = L / N; % 网格间距
    x = linspace(-L/2, L/2, N); % 生成从 -L/2 到 L/2 的N个点
    y = x; % 通常x和y方向尺寸对称
    [X, Y] = meshgrid(x, y); % 生成二维网格
end

在主脚本 main.m 的开头，我们先定义光束的物理参数，并调用这个函数来创建我们的“计算画布”。

% main.m (部分)
clear; close all; clc;

%% 1. 定义基本参数
lambda = 632.8e-9;     % 波长，例如氦氖激光的632.8纳米
w0 = 1e-3;             % 束腰半径 (