当前位置：首页 > news >正文

从基模到Ince高斯光束：MATLAB仿真中的模式演化与参数调控

news 2026/3/30 22:43:48

1. Ince高斯光束：光学世界里的"变形金刚"

第一次在实验室看到Ince高斯光束的仿真结果时，我盯着屏幕上那些不断变换的光斑图案看了足足十分钟。这种光束就像光学领域的"变形金刚"，只需调整几个关键参数，就能在厄密高斯光束和拉盖尔高斯光束之间自由切换。对于从事激光光学研究的朋友来说，掌握Ince高斯光束的仿真技术，就相当于获得了一把打开复杂光场调控大门的钥匙。

Ince高斯光束本质上是在椭圆坐标系下求解傍轴波动方程得到的特殊解。与常见的基模高斯光束相比，它多了三个关键控制参数：u阶数、v阶数和椭圆率e。这三个参数就像调节旋钮，u和v决定光束的空间模式结构，e则控制着光束从圆形到椭圆形的连续变化。当e趋近于0时，光束会展现出拉盖尔高斯光束的特性；当e趋近于无穷大时，又会转变为厄密高斯光束的模式。

在实际应用中，这种灵活的模式转换特性让Ince高斯光束在光镊、激光加工等领域大显身手。比如在光学微操纵中，不同形状的光束可以更精准地捕获不同形状的微粒；在激光材料加工中，椭圆光束可能更适合某些特殊形状的加工需求。而MATLAB作为强大的数值计算工具，为我们提供了一套完整的仿真实验平台，无需昂贵的光学设备就能探索这些奇妙的光学现象。

2. MATLAB仿真环境搭建与参数设置

2.1 基础参数配置：从波长到坐标系

开始仿真前，我们需要先搭建好MATLAB的计算环境。建议使用R2018b或更新版本，这些版本对矩阵运算和图形显示做了大量优化。在我的实际使用中，发现2018b版本在运行本文的椭圆坐标系转换时，速度比早期版本快了近30%。

首先定义基本光学参数，这段代码设置了交互式输入界面，方便反复调试：

parity = input('Ince高斯光束模式(奇1或偶0)：'); lambda = input('入射光波波长(mm)='); z = input('传播距离='); v = input('阶数v='); u = input('阶数u='); e = input('椭圆率e=');

这里有几个经验参数需要注意：

波长λ通常以毫米为单位输入，但实际计算中会自动转换为米制
传播距离z的取值建议在0-2倍瑞利长度之间，这个范围内光束变化最明显
椭圆率e的范围理论上为0到∞，但实际仿真时e=0.01～100已经能展现完整变化

接下来计算衍生参数，这部分直接关系到仿真的物理准确性：

k = 2*pi/lambda; % 波数 w0 = 1.5*lambda; % 束腰半径 ZR = k*w0^2/2; % 瑞利长度 wz = w0*sqrt(1+(z/ZR)^2); % z处束宽 Rz = z + ZR^2/z; % 波前曲率半径 Phi = atan(z/ZR); % 古伊相位

2.2 参数校验：避免常见的输入错误

在调试过程中，我发现大约40%的仿真失败都源于参数输入不当。Ince多项式对u和v的关系有严格要求，必须添加校验代码：

if parity==0 if (v<0)||(v>u) error('偶模式：0<=v<=u'); end else if (v<1)||(v>u) error('奇模式：1<=v<=u'); end end if (-1)^(v-u)~=1 error('u和v必须同奇偶'); end

特别提醒：当u和v都较大时（如u=15，v=13），计算量会急剧增加。在我的i7-11800H笔记本上，u超过20后仿真时间可能达到分钟级。建议初学者先从u=3～5的小参数开始实验。

3. 椭圆坐标系构建与Ince多项式实现

3.1 从笛卡尔坐标到椭圆坐标的魔法转换

Ince高斯光束的核心在于椭圆坐标系的使用。这个坐标系就像给光场戴上了一副"椭圆眼镜"，让我们能看到隐藏在常规直角坐标系中的特殊模式。转换函数如下：

function [xi,eta,x1,y1] = mesh_elliptic(f,L,N) [x1,y1] = meshgrid(linspace(-L,L,N)); xi = zeros(N); eta = zeros(N); en = acosh((x1(1:(N+1)/2,(N+1)/2:N)+1i*y1(1:(N+1)/2,(N+1)/2:N))/f); ee = real(en); nn = imag(en); nn = nn + (nn<0)*2*pi; xi(1:(N+1)/2,(N+1)/2:N) = ee; eta(1:(N+1)/2,(N+1)/2:N) = nn; % 其他象限处理（略） end

实际调用时，坐标尺寸L建议取4-5倍波长，取样点数N根据电脑性能选择。在我的测试中，N=1001能在精度和速度间取得较好平衡。值得注意的是，椭圆焦距f会随传播距离z变化：

f0 = sqrt(e/2)*w0; % z=0处椭圆焦距 f = f0*wz/w0; % z处椭圆焦距

3.2 Ince多项式：光束模式的"基因编码"

Ince多项式决定了光束的空间结构特性，分为偶多项式CInce和奇多项式SInce。以偶多项式为例，其实现核心是构建特征方程并求解：

function [IP,eta,coef,dIP] = CInce(u,v,q,z) if mod(u,2)==0 j = u/2; N = j+1; n = v/2+1; M = diag(q*(j+(1:N-1)),1) + diag([2*q*j,q*(j-(1:N-2))],-1)... + diag([0,4*((0:N-2)+1).^2]); [A,ets] = eig(M); ets = diag(ets); [ets,index] = sort(ets); A = A(:,index); % 归一化处理（略） r = 0:N-1; [R,X] = meshgrid(r,z); IP = cos(2*X.*R)*A(:,n); dIP = -2*R.*sin(2*X.*R)*A(:,n); else % u为奇数的情况（略） end end

调试技巧：当遇到"NaN"或异常结果时，建议先检查q值（即椭圆率e）是否过大。在我的经验中，e>100时可能需要调整归一化方式。另外，奇偶模式的选择会显著影响最终光强分布，建议对同一组(u,v,e)参数分别运行偶、奇模式，对比观察差异。

4. 光束合成与可视化：从公式到图像

4.1 电场表达式合成：拼出完整的光学拼图

有了前面的基础，现在可以组装完整的Ince高斯光束表达式了。这个过程就像用不同形状的积木搭建独特的光学结构：

E_Gauss = (w0./wz).*exp(-(r.^2./wz^2)).*exp(1i.*(k.*z+k.*r.^2./(2*Rz)-Phi)); if z==0 if parity==0 E_Ince = CInce(u,v,e,etha).*CInce(u,v,e,1i*xhi).*exp(-(r/w0).^2); else E_Ince = SInce(u,v,e,etha).*SInce(u,v,e,1i*xhi).*exp(-(r/w0).^2); end else if parity==0 E_Ince = (CInce(u,v,e,etha).*CInce(u,v,e,1i*xhi)).*E_Gauss.*exp((-u+1).*Phi); else E_Ince = (SInce(u,v,e,etha).*SInce(u,v,e,1i*xhi)).*E_Gauss.*exp((-u+1).*Phi); end end

归一化处理是保证不同模式间可比性的关键。这部分代码虽然看起来复杂，但其实是在计算一个与u、v、e相关的比例因子：

if parity == 0 if mod(u,2)==0 [C0,~,coef] = CInce(u,v,e,0); [Cp,~,~] = CInce(u,v,e,pi/2); Norm = (-1)^(v/2)*sqrt(2)*gamma(u/2+1)*coef(1)*sqrt(2/pi)/w0/C0/Cp; else % 奇数u归一化（略） end else % 奇模式归一化（略） end E_Ince = E_Ince*Norm;

4.2 三维可视化：让光场跃然屏上

MATLAB提供了丰富的可视化工具，这里我们使用三种方式展示Ince高斯光束：

I_Ince = E_Ince.*conj(E_Ince); normI_Ince = I_Ince/max(max(I_Ince)); figure() pcolor(x1,y1,normI_Ince); shading interp; colormap hot; colorbar; figure() pcolor(x1,y1,angle(E_Ince)); clim([-pi,pi]); % 相位范围固定 shading interp; colormap hot; colorbar; figure() surf(x1,y1,normI_Ince); shading interp; colormap hot; colorbar;

绘图技巧：