MATLAB仿真：涡旋光束经大气湍流后的光强及相位分布基于蒙特卡洛相位屏法

在自由空间涡旋光通信中，大气作为传输信道，由于湍流效应的影响，涡旋光束相位产生畸变，导致在接收端接收到的光束信号质量下降，无法准确的识别光信息。在研究湍流传输问题时，需要建立符合实际的大气湍流环境模型，对于湍流的建模，数值模拟是一种非常重要并且有效的方法。当光束在大气湍流中传输时，可以认为大气湍流是类似于一种屏的形式对传输的光束进行影响，这种方法就是相位屏仿真方法。目前生成相位屏的方法中，使用最多的是泽尼克多项式法和功率谱反演法。本文的原理部分参考文章《大气湍流中涡旋光束轨道角动量的高精度识别_史晨寅》

大气湍流原理：

本文产生大气湍流相位屏的方式为功率谱反演法，该方法以快速傅里叶变换为基础。首先利用湍流折射率功率谱对随机复高斯矩阵进行滤波，然后进行快速傅里叶变换来得到湍流的相位扰动。相位功率谱Fq(k)和各向异性折射率功率Fn(k)谱的关系为：

式中，k=2p/l为波数，Dz为相位屏之间的间隔。通过选择不同的折射率模型，如柯尔莫哥洛夫、冯卡曼和改进的冯卡曼模型等可以得到对应的相位功率谱，分别为：

这篇文章中我们利用改进的冯卡曼模型。其中km=5.92/l0，k0=2p/L0，l0为内尺度，L0为外尺度，通过选择km和k0的值是的量纲分析预测的小尺度和大尺度行为相匹配。r0为大气相干直径也称弗里德参数，可以表示为：

Cn2为大气结构常数，可以表示大气湍流的强弱。我们在这里利用蒙特卡洛相位屏法进行模拟。通过在传播路径上设置多个湍流屏模拟光束经过大气湍流的过程，模拟过程以及产生的湍流屏如图1所示。

图1：多个湍流屏模拟光束经过大气湍流的过程以及湍流相位屏

仿真结果展示：

我们设置入射光束为拉盖尔高斯光束：

固定参数：L0=1.5; l0=1.0e-3; z=1000; w=0.05m; lambda=632.8e-9;湍流屏个数设置10个。如图2展示了拓扑荷为2的拉盖尔高斯涡旋光束在不同的大气结构常数下的光强和相位分布：可以看到随着大气结构常数增加，大气湍流强度越强，光束受到的影响越大，相位畸变也越大。

图2：拓扑荷为2的拉盖尔高斯涡旋光束在不同的大气结构常数下的光强和相位分布

图3展示了拓扑荷为3、4、5的涡旋光束在不同大气结构常数下的光强分布，图中可以看出随着拓扑荷增大在强湍流情况下，大的拓扑荷仍然可以保持环状结构。

图3：拓扑荷为3、4、5的涡旋光束在不同大气结构常数下的光强分布

MATLAB仿真代码：

喜欢的话给小编点个关注点个赞呀

有任何问题可以在评论区留言

更多完整代码进T店博士生牛马赚点小钱

clc;clear;
%% 坐标
N=512;
delta1=0.5/N;
[x,y]=meshgrid((-N/2:1:N/2-1).*delta1);
r1 = sqrt(x.^2+y.^2);
phi= angle(x+1i*y);
%% 光束参数
lambda=632.8e-9;%m
k = 2 * pi / lambda;
w=0.05;
l=3;
E1= (sqrt(r1.^2)./w).^abs(l).*exp(-(r1.^2)./(w.^2)).*exp(1i.*l.*phi);
%%  传输
Cn2=1.0e-13;
L0=1.5;
l0=1.0e-3;
z=1000;
n=10;
dz=z/n;
r0=(0.423 * k^2 * Cn2 *dz) .^(-3/5);
更多完整代码见t店博士生牛马赚点小钱