Matlab模拟OAM螺旋谱展开，包括光束本征态、各类湍流、衍射、干涉下的OAM谱分布

Matlab模拟OAM螺旋谱展开，包括光束本征态、各类湍流、衍射、干涉下的OAM谱分布

江湖传闻有个叫OAM（轨道角动量）的光学黑科技，能用来搞量子通信和超分辨率成像。这玩意儿最骚的操作就是螺旋相位分布，今天咱们用Matlab整点活，看看不同场景下OAM频谱怎么蹦迪。

先整点基础款——拉盖尔-高斯光束生成。这货的相位像披萨旋转，每圈相位转2πl次：

lambda = 632e-9; % 红光 w0 = 1e-3; % 束腰 l = 3; % 拓扑荷数 x = linspace(-5e-3,5e-3,512); [X,Y] = meshgrid(x); [phi,r] = cart2pol(X,Y); E = (r/w0).^abs(l) .* exp(-r.^2/w0^2) .* exp(1i*l*phi); % 振幅+相位 imagesc(angle(E)); % 显示相位涡旋

这段代码里cart2pol把直角坐标转极坐标是关键，l参数控制相位旋转速度。跑出来的相位图像蚊香盘似的，l越大螺旋越密。

接下来搞事情——大气湍流模拟。用随机相位屏整活：

Cn2 = 1e-14; % 湍流强度 L0 = 10; % 外尺度(m) N = 512; delta = 10e-3/N; [kx,ky] = meshgrid((-N/2:N/2-1)*delta); k_sq = kx.^2 + ky.^2; phase_psd = 0.033*(2*pi)^2*Cn2*k_sq.^(-11/6).*exp(-k_sq/(k0^2)); % 功率谱密度 turb_phase = real(ifft2(fft2(randn(N)).*sqrt(phase_psd))); % 生成随机相位 E_turb = E .* exp(1i*2*pi*turb_phase); % 叠加湍流

这里的骚操作是用高斯随机场经过功率谱滤波生成湍流相位。k_sq.^(-11/6)对应Kolmogorov湍流模型，像给光波加了动态马赛克。

衍射效应更带劲，直接上角谱法：

z = 1000; % 传输距离(m) k = 2*pi/lambda; H = exp(1i*k*z*sqrt(1-(lambda*kx).^2 - (lambda*ky).^2)); % 传递函数 E_far = fftshift(ifft2(fft2(ifftshift(E)) .* H)); % 角谱衍射

fftshift和ifftshift这对好基友处理频域位移，H函数实现近场到远场变换。传输后的光斑会像被熊孩子撕开的棉花糖。

重头戏是OAM谱分析，用方位角傅里叶变换：

[~,r_idx] = max(max(abs(E),[],1)); % 找最大环 E_ring = E(:,r_idx); spectrum = fft(E_ring);

这个操作相当于把光环展开成直线再做FFT，l=3时会在第四个频点出现尖峰（Matlab下标从1开始）。就像把蚊香拉直了数圈数。

最后整点高阶玩法——多模式干涉：

E1 = gen_OAM(3); % 生成l=3模式 E2 = gen_OAM(-2); % 生成l=-2模式 interferogram = abs(E1 + E2).^2; % 干涉图样 spectrum = fftshift(fft2(interferogram)); % 频谱分析

这时候频谱会出现3-2=1和3+2=5的边带，像两个不同转速的齿轮咬合产生的振动模式。搞通信的可以拿这个做模式分复用，量子佬们能用来搞纠缠态制备。

实战中发现个玄学现象：强湍流下OAM谱会像摇骰子一样弥散，这时候需要自适应光学或者深度学习来矫正。有个邪道玩法是用卷积神经网络直接从光斑图预测拓扑荷数，准确率能到90%以上，不过那是另外的价钱了。