CN118094964A - 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 - Google Patents
基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 Download PDFInfo
- Publication number
- CN118094964A CN118094964A CN202410508362.5A CN202410508362A CN118094964A CN 118094964 A CN118094964 A CN 118094964A CN 202410508362 A CN202410508362 A CN 202410508362A CN 118094964 A CN118094964 A CN 118094964A
- Authority
- CN
- China
- Prior art keywords
- photon
- scattering
- phase screen
- photons
- updating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 94
- 230000005540 biological transmission Effects 0.000 title claims abstract description 40
- 239000013598 vector Substances 0.000 title claims abstract description 32
- 238000000342 Monte Carlo simulation Methods 0.000 title claims abstract description 23
- 238000001228 spectrum Methods 0.000 claims abstract description 15
- 239000013535 sea water Substances 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 16
- 230000010287 polarization Effects 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 description 12
- 230000003287 optical effect Effects 0.000 description 10
- 230000000694 effects Effects 0.000 description 7
- 238000010521 absorption reaction Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Optical Communication System (AREA)
Abstract
本发明公开了基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，包括生成光子的初始坐标与方向；确定入射光的Stokes参量和散射方向；更新Stokes参量和散射方向；选定海洋功率谱，根据功率谱反演方法和高阶次谐波补偿法，联合生成仿真所用湍流相位屏；计算光子的随机步长；根据光子坐标、光子方向余弦与随机步长，判断光子是否与相位屏相交、是否达到接收面；若光子未与相位屏相交、未到达接收面，则更新碰撞后光子的坐标、方向、权重及Stokes矢量；若光子与相位屏相交，未到达接收面，则计算光子与相位屏相交的坐标，更新光子坐标、方向及stokes矢量以及光子权重。本发明更加全面的模拟偏振光在海水中的传输过程。
Description
技术领域
本发明涉及水下激光偏振探测领域，尤其是基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法。
背景技术
光学探测是人类对水下资源探测的主要手段之一，20 世纪 60 年代海水中发现蓝绿透射窗口以及光波相对声波方向性好、定位准确和准直度高的优势，为水下光学探测提供了理论基础。蓝绿激光具有穿透能力强、探测精度高等优点，被广泛应用于水深探测、海底地形地貌绘制及目标探测等领域，蓝绿激光海洋探测技术已成为一项极具发展潜力的热点技术。
海水对光信号的传播产生的复杂影响，包括由悬浮粒子、浮游生物和溶解盐引起的光信号吸收和散射效应，以及海水盐度、温度和密度变化导致的湍流效应，这些因素共同作用，导致光传输路径扭曲和光的偏振态随机变化，进而影响水下光通信和偏振光探测技术的效率和可靠性。吸收会导致光信号的能量衰减，而散射则会使光线偏离原始传播路径，引起到达接收器的光强减弱。由于散射路径的多样性，还会引入时延，导致信号失真。另外，光与水中的微粒相互作用会导致光的偏振态随机变化，在水下激光偏振探测系统中，偏振信息是重要的信号属性之一。如果偏振态在传播过程中发生变化，那么在接收端，原始的偏振信息可能部分或完全丢失，从而降低信号的可识别性和信息传输的准确性。湍流会影响光的传播方向，导致光路径产生扭曲，从而影响到达接收器的光信号的质量和强度。在极端情况下，湍流可以导致光束的显著偏移，极大降低接收信号的可靠性。
真实海水环境对传输光信号的影响是光吸收散射、湍流效应的综合效果。目前的方法是基于广义折射定律将湍流效应融入到模拟吸收散射的蒙特卡洛仿真框架下进行建模研究。但现有的这种复合传输模型，对于关注水下光传输偏振特性的研究工作并不适用。
发明内容
为了克服现有技术中存在的上述问题，本发明提出基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法。
本发明解决其技术问题所采用的技术方案是：基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，包括如下步骤：
步骤1，建立三维坐标系，设定Z轴正方向为光传输方向，Z＝0为光源所在平面，进行模拟参数设置，确定系统初始化；
步骤2，根据高斯分布的特性来采样光子的初始位置和方向；
步骤3，确定入射光的stokes参量，舍选法确定散射角θ和方位角；
步骤4，更新Stokes参量和散射方向；
步骤5，选定海洋功率谱，根据功率谱反演方法和高阶次谐波补偿法，联合生成海洋湍流相位屏;
步骤6，计算光子的随机步长；
步骤7，根据光子坐标、光子方向余弦与随机步长，判断光子是否与相位屏相交、是否达到接收面；
步骤8，若光子未与相位屏相交、未到达接收面，则更新碰撞后光子的坐标、方向、权重及Stokes矢量；
步骤9，若光子与相位屏相交，未到达接收面，则计算光子与相位屏相交的坐标，更新光子坐标、方向及stokes矢量以及光子权重；
步骤10，循环执行步骤6-步骤9，直到光子与接收面相交。
上述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，所述步骤2具体包括：
步骤2.1，生成在(0,1)上均匀分布的随机数和；
步骤2.2，确定初始散射角；
其中，λ为波长，为初始束宽；
步骤2.3，确定初始方位角；
步骤2.4，确定距离；
步骤2.5，确定光子初始坐标为：
；
步骤2.6，确定光子初始方向为：
；
步骤2.7，步骤2.1-步骤2.6，即可生成所有光子的初始坐标和初始方向。
上述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，所述步骤3具体包括：设定为入射波总光强，为入射波的x和y两个方向光强的相差值；为± 45°线偏振光分量；表示偏振光中圆偏振的比例，入射波的Stokes矢量可以表示为，其联合概率密度函数可表示为：
；
其中，θ为散射角，是光子传播方向与X轴正半轴之间夹角；β为方位角，是光子散射方向在XOY面上的投影与X轴正半轴的夹角；参数和为Mueller矩阵的两个元素，当散射体为各向同性球体时，散射矩阵可以化简为：
；
参数、、、可以由散射振幅和得到：
；
其中，散射振幅和的大小与光子半径和光子复折射率有关；
选定舍选法选择出下一次散射的方位角和散射角θ，对联合概率密度函数中的两个随机变量，先均匀抽样，分别为第t时刻抽样的散射角、方位角和联合概率密度函数，其对应的取值范围分别为(0,π)、(0,2π)和(0,1)；若，则选取成功，进入步骤4，否则重新选取。
上述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，所述步骤4具体包括：
步骤4.1，从入射子午面旋转到散射面,引入旋转矩阵，将子午面的Stokes参数旋转到散射面入射,旋转后的Stokes参量由入射Stokes参量S与旋转矩阵相乘获得，其中；
步骤4.2，将Stokes参量在散射面进行两点转动参考面不变，此时的转换矩阵为Mueller矩阵；同时更新光子的方向：
当时，
；
其中，为散射前光子的方向；
当时，
；
步骤4.3，由散射面回到出射的新的子午面，此时的旋转矩阵，旋转角的表达式：
；
最终散射后的Stokes参量表示为：；
式中，表示散射前的Stokes参量；由于Stokes参量是一个统计结果，最后总的偏振光的状态为：
；
其中，为光子数。
上述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，所述步骤5具体包括：
步骤5.1，选定海洋湍流功率谱模型：
；
其中，k为空间波数；表示海水单位质量湍流动能耗散率；是Kolmogorov内尺度；；为均方温度耗散率；表示海洋湍流功率谱中温度和盐度起伏引起折射率变化贡献的比率；
步骤5.2，利用湍流折射率谱和复高斯随机数矩阵产生相空间复随机场，进行逆傅里叶变换获得二维相位的空间分布：
；
式中，和分别代表x方向和y方向的离散空间频率，，，其中和分别代表x方向和y方向的格点尺寸；表示高斯随机分布且均值为零的复随机矩阵；为相位屏在z轴上的间隔，k为空间波数；
步骤5.3，设定谐波次数，选用低频次谐波补偿法对相位屏的进行低频补偿，生成新的相位屏，表达式为：
；
式中，为相位屏低频补偿项，p为次谐波级数，且次谐波频率间隔；
步骤5.4，重复步骤5.1-步骤5.3，即可生成全部相位屏。
上述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，所述步骤6具体计算公式为：
光子的随机步长；
其中，为在(0,1)之间均匀分布的伪随机数；c为衰减系数；
若光子未碰撞之前的位置为和传播方向z上的角余弦为，光子位置的更新方法为：
。
上述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，所述步骤8具体包括：利用舍选法选择出散射的方位角和散射角θ，更新光子的方向：
当时，
；
当时，
；
更新Stokes参量：
；
，，；
；
更新光子坐标：
；
更新光子权重：
；
式中，为更新之前的光子权重，bb为散射系数，c为衰减系数。
上述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，所述步骤9具体包括：计算光子与相位屏相交的坐标为：
；
其中，和分别为穿过光子相位屏前的坐标和方向，为相位屏在z轴方向上的位置；
光子穿过相位屏时，其传输方向会发生变化，根据广义斯涅尔定律，新的光子方向表示为：
；
式中，为光波数；为相位屏的空间相位函数；
光子穿过相位屏后，下次散射发生前的光子坐标表示为：
；
更新光子的方向：
当时，
；
当时，
；
同时更新Stokes参量：；
，，
，；
更新光子权重：
；
其中，为更新之前的光子权重，bb为散射系数，c为衰减系数。
本发明的有益效果是，本发明能够同步考虑水体的衰减和湍流效应对偏振光传输的影响，更加全面的模拟偏振光在海水中的传输过程；
本发明模拟过程中采用高阶次谐波补偿法对基于FFT法生成的湍流相位屏进行低频补偿，重新生成的湍流相位屏更加真实可靠；
本发明通过对光子偏振状态的追踪，可以获取接收端光子的偏振状态，能够更加直观、全面、系统地分析水体的光学传输特性对水下偏振探测系统的影响。可为水下偏振探测系统的设计优化提供可视化参考依据。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1为本发明原理示意图；
图2为本发明整体流程图。
具体实施方式
为使本领域技术人员更好的理解本发明的技术方案，下面结合附图和具体实施方式对本发明作详细说明。
参见图1所示，本发明一种基于湍流相位屏的水下偏振光传输的矢量蒙特卡洛模拟方法，包括以下步骤：
步骤1，模拟系统初始化；以Z轴正方向为光传输方向，发射光源所在平面为Z＝0，设置光传输距离为18m，波长λ为550nm；光子数N为108；初始束宽为30mm，在Z轴方向上的相位屏间隔Δz=2m，位置分别为z＝2m ,z＝4m ,z＝6m ,z＝8m ,z＝10m ,z＝12m ,z＝14m, z＝16m；相位屏尺寸为1.1m× 1.1m；方位角为光子散射方向在XOY面上的投影与X轴正半轴的夹角；散射角θ为光子传播方向与X轴正半轴之间夹角；设定海水的吸收系数，散射系数bb=0.26m-1，衰减系数c=0.46m-1；光子的初始权重为W0=1。
步骤2，根据高斯分布的特性来采样光子的初始位置和方向，方法如下：
所述步骤2具体包括：
步骤2.1，生成在(0,1)上均匀分布的随机数和；
步骤2.2，确定初始散射角；
其中，λ为波长，为初始束宽；
步骤2.3，确定初始方位角；
步骤2.4，确定距离；
步骤2.5，确定光子初始坐标为：
；
步骤2.6，确定光子初始方向为：
；
步骤2.7，步骤2.1-步骤2.6，即可生成所有光子的初始坐标和初始方向。
步骤3，确定散射角θ和方位角；设定为入射波总光强，为入射波的x和y两个方向光强的相差值；为± 45°线偏振光分量；表示偏振光中圆偏振的比例，入射波的Stokes矢量可以表示为，其联合概率密度函数可表示为：
；
由上式可得偏振光由于是非轴对称的，线偏振光的联合概率密度函数与散射角θ和方位角都有关联，圆偏振光因为其轴对称性，联合概率密度函数仅与θ有关。式中，参数参数和为Mueller矩阵的两个元素，当散射体为各向同性球体时，散射矩阵可以化简为：
；
参数、、、可以由散射振幅和得到：
；
其中，散射振幅和的大小与光子半径和光子复折射率有关；
选定舍选法选择出下一次散射的方位角和散射角θ，对联合概率密度函数中的两个随机变量，先均匀抽样，分别为第t时刻抽样的散射角、方位角和联合概率密度函数，其对应的取值范围分别为(0,π)、(0,2π)和(0,1)；若，则选取成功，进入步骤4，否则重新选取。
步骤4，更新Stokes参量和散射角；方法如下：
步骤4.1，从入射子午面旋转到散射面,引入旋转矩阵，将子午面的Stokes参数旋转到散射面入射,旋转后的Stokes参量由入射Stokes参量S与旋转矩阵相乘获得，其中；
步骤4.2，将Stokes参量在散射面进行两点转动参考面不变，此时的转换矩阵为Mueller矩阵；同时更新光子的方向：
当时，
；
其中，为散射前光子的方向；
当时，
；
步骤4.3，由散射面回到出射的新的子午面，此时的旋转矩阵，旋转角的表达式：
；
最终散射后的Stokes参量表示为：；
式中，表示散射前的Stokes参量；由于Stokes参量是一个统计结果，最后总的偏振光的状态为：
；
其中，为光子数。
步骤5，生成海洋湍流相位屏；方法如下：
步骤5.1，选定海洋湍流功率谱模型：
；
其中，k为空间波数；是Kolmogorov内尺度；；表示海水单位质量湍流动能耗散率，取值范围为；为均方温度耗散率，取值范围为；表示海洋湍流功率谱中温度和盐度起伏引起折射率变化贡献的比率，其取值范围为[-5,0]；本实例中,,,。
步骤5.2，利用湍流折射率谱和复高斯随机数矩阵产生相空间复随机场，进行逆傅里叶变换获得二维相位的空间分布：
；
式中，和分别代表x方向和y方向的离散空间频率，，，其中和分别代表x方向和y方向的格点尺寸；表示高斯随机分布且均值为零的复随机矩阵；为相位屏在z轴上的间隔，k为空间波数；
步骤5.3，设定谐波次数，选用低频次谐波补偿法对相位屏的进行低频补偿，生成新的相位屏，表达式为：
；
式中，为相位屏低频补偿项，p为次谐波级数，且次谐波频率间隔；
步骤5.4，重复步骤5.1-步骤5.3，即可生成全部相位屏。
步骤6的具体实施方式如下：
计算光子的随机步长光子的随机步长；
其中，为在(0,1)之间均匀分布的伪随机数；c为衰减系数；
若光子未碰撞之前的位置为和传播方向z上的角余弦为，光子位置的更新方法为：
。
步骤7，根据光子坐标、光子方向余弦与随机步长L，判断光子是否与相位屏相交、是否达到接收面。
步骤8，若光子未与相位屏相交、未到达接收面，则：
利用舍选法选择出散射的方位角和散射角θ，更新光子的方向：
当时，
；
当时，
；
更新Stokes参量：
；
，，；
；
更新光子坐标：
；
更新光子权重：
；
式中，为更新之前的光子权重，bb为散射系数，c为衰减系数。
步骤9，若光子与相位屏相交，未到达接收面，则
步骤9.1，计算光子与相位屏相交的坐标为：
；
其中，和分别为穿过光子相位屏前的坐标和方向，为相位屏在z轴方向上的位置；
光子穿过相位屏时，其传输方向会发生变化，根据广义斯涅尔定律，新的光子方向表示为：
；
式中，为光波数；为相位屏的空间相位函数；
光子穿过相位屏后，下次散射发生前的光子坐标表示为：
；
更新光子的方向：
当时，
；
当时，
；
同时更新Stokes参量：；
，，
，；
更新光子权重：
；
其中，为更新之前的光子权重，bb为散射系数，c为衰减系数。
步骤10中，循环执行上述过程步骤6~步骤9，直至光子与接收面相交，确认并记录光子与接收面相交的位置坐标、stokes矢量以及光子权重。
实施例
本实施例中的步骤1中，光传输方向为Z轴正方向，光源所在平面为Z＝0平面，设置光传输距离为100m，波长λ为480nm；光子数N为107；初始束宽为40mm，在Z轴方向上的相位屏间隔Δz=10m，位置分别为z＝10m ,z＝20m ,z＝30m ,z＝40m ,z＝50m ,z＝60m ,z＝70m, z＝80m，z＝90m，z＝100m；相位屏尺寸为2m× 2m；设定海水的吸收系数，散射系数bb=0.28m-1，衰减系数c=0.51m-1；光子的初始权重为W0=1。
本实施例中的步骤5中，；；；；。
其余各步骤与实施例1完全相同。
以上实施例仅为本发明的示例性实施例，不用于限制本发明。本领域技术人员可以在本发明的实质和保护范围内，对本发明做出各种修改或等同替换，这种修改或等同替换也应视为落在本发明的保护范围内。
Claims (8)
1.基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，包括如下步骤：
步骤1，建立三维坐标系，设定Z轴正方向为光传输方向，Z＝0为光源所在平面，进行模拟参数设置，确定系统初始化；
步骤2，根据高斯分布的特性来采样光子的初始位置和方向；
步骤3，确定入射光的stokes参量，舍选法确定散射角θ和方位角β；
步骤4，更新Stokes参量和散射方向；
步骤5，选定海洋功率谱，根据功率谱反演方法和高阶次谐波补偿法，联合生成海洋湍流相位屏;
步骤6，计算光子的随机步长；
步骤7，根据光子坐标、光子方向余弦与随机步长，判断光子是否与相位屏相交、是否达到接收面；
步骤8，若光子未与相位屏相交、未到达接收面，则更新碰撞后光子的坐标、方向、权重及Stokes矢量；
步骤9，若光子与相位屏相交，未到达接收面，则计算光子与相位屏相交的坐标，更新光子坐标、方向及stokes矢量以及光子权重；
步骤10，循环执行步骤6-步骤9，直到光子与接收面相交。
2.根据权利要求1所述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，所述步骤2具体包括：
步骤2.1，生成在(0,1)上均匀分布的随机数和；
步骤2.2，确定初始散射角；
其中，λ为波长，为初始束宽；
步骤2.3，确定初始方位角；
步骤2.4，确定距离；
步骤2.5，确定光子初始坐标为：
；
步骤2.6，确定光子初始方向为：
；
步骤2.7，步骤2.1-步骤2.6，即可生成所有光子的初始坐标和初始方向。
3.根据权利要求1所述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，所述步骤3具体包括：设定为入射波总光强，为入射波的x和y两个方向光强的相差值；为± 45°线偏振光分量；表示偏振光中圆偏振的比例，入射波的Stokes矢量可以表示为，其联合概率密度函数可表示为：
；
其中，θ为散射角，是光子传播方向与X轴正半轴之间夹角；β为方位角，是光子散射方向在XOY面上的投影与X轴正半轴的夹角；参数和为Mueller矩阵的两个元素，当散射体为各向同性球体时，散射矩阵可以化简为：
；
参数、、、可以由散射振幅和得到：
；
其中，散射振幅和的大小与光子半径和光子复折射率有关；
选定舍选法选择出下一次散射的方位角和散射角θ，对联合概率密度函数中的两个随机变量，先均匀抽样，分别为第t时刻抽样的散射角、方位角和联合概率密度函数，其对应的取值范围分别为(0,π)、(0,2π)和(0,1)；若，则选取成功，进入步骤4，否则重新选取。
4.根据权利要求1所述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，所述步骤4具体包括：
步骤4.1，从入射子午面旋转到散射面,引入旋转矩阵，将子午面的Stokes参数旋转到散射面入射,旋转后的Stokes参量由入射Stokes参量S与旋转矩阵相乘获得，其中；
步骤4.2，将Stokes参量在散射面进行两点转动参考面不变，此时的转换矩阵为Mueller矩阵；同时更新光子的方向：
当时，
；
其中，为散射前光子的方向；
当时，
；
步骤4.3，由散射面回到出射的新的子午面，此时的旋转矩阵，旋转角的表达式：
；
最终散射后的Stokes参量表示为：；
式中，表示散射前的Stokes参量；由于Stokes参量是一个统计结果，最后总的偏振光的状态为：
；
其中，为光子数。
5.根据权利要求1所述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，所述步骤5具体包括：
步骤5.1，选定海洋湍流功率谱模型：
；
其中，k为空间波数；表示海水单位质量湍流动能耗散率；是Kolmogorov内尺度；；为均方温度耗散率；表示海洋湍流功率谱中温度和盐度起伏引起折射率变化贡献的比率；
步骤5.2，利用湍流折射率谱和复高斯随机数矩阵产生相空间复随机场，进行逆傅里叶变换获得二维相位的空间分布：
；
式中，和分别代表x方向和y方向的离散空间频率，，，其中和分别代表x方向和y方向的格点尺寸；表示高斯随机分布且均值为零的复随机矩阵；为相位屏在z轴上的间隔，k为空间波数；
步骤5.3，设定谐波次数，选用低频次谐波补偿法对相位屏的进行低频补偿，生成新的相位屏，表达式为：
；
式中，为相位屏低频补偿项，p为次谐波级数，且次谐波频率间隔；
步骤5.4，重复步骤5.1-步骤5.3，即可生成全部相位屏。
6.根据权利要求1所述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，所述步骤6具体计算公式为：
光子的随机步长；
其中，为在(0,1)之间均匀分布的伪随机数；c为衰减系数；
若光子未碰撞之前的位置为和传播方向z上的角余弦为，光子位置的更新方法为：
。
7.根据权利要求1所述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，所述步骤8具体包括：利用舍选法选择出散射的方位角和散射角θ，更新光子的方向：
当时，
；
当时，
；
更新Stokes参量：
；
，，；
其中；
更新光子坐标：
；
更新光子权重：
；
式中，为更新之前的光子权重，bb为散射系数，c为衰减系数。
8.根据权利要求1所述的基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法，其特征在于，所述步骤9具体包括：计算光子与相位屏相交的坐标为：
；
其中，和分别为穿过光子相位屏前的坐标和方向，为相位屏在z轴方向上的位置；
光子穿过相位屏时，其传输方向会发生变化，根据广义斯涅尔定律，新的光子方向表示为：
；
式中，为光波数；为相位屏的空间相位函数；
光子穿过相位屏后，下次散射发生前的光子坐标表示为：
；
更新光子的方向：
当时，
；
当时，
；
同时更新Stokes参量：；
，，
，；
更新光子权重：
；
其中，为更新之前的光子权重，bb为散射系数，c为衰减系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410508362.5A CN118094964B (zh) | 2024-04-26 | 2024-04-26 | 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410508362.5A CN118094964B (zh) | 2024-04-26 | 2024-04-26 | 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118094964A true CN118094964A (zh) | 2024-05-28 |
| CN118094964B CN118094964B (zh) | 2024-07-23 |
Family
ID=91155098
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410508362.5A Active CN118094964B (zh) | 2024-04-26 | 2024-04-26 | 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118094964B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6958812B1 (en) * | 2003-08-19 | 2005-10-25 | The United States Of America As Represented By The Secretary Of The Army | Systems and methods for analyzing particle systems of surface facets using polarized scattered light |
| CN108763610A (zh) * | 2018-03-28 | 2018-11-06 | 华北电力大学（保定） | 一种基于谱反演法的大气湍流相位屏仿真方法 |
| CN111555822A (zh) * | 2020-04-28 | 2020-08-18 | 西安邮电大学 | 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 |
| AU2020102396A4 (en) * | 2020-09-23 | 2020-11-05 | Changchun University Of Science And Technology | Machine learning-based method for analyzing characteristics of laser beam propagation through turbulent atmosphere |
| CN113868847A (zh) * | 2021-09-14 | 2021-12-31 | 西安理工大学 | 一种偏振光在火星沙尘暴中传输的蒙特卡洛计算方法 |
| CN114741903A (zh) * | 2022-05-18 | 2022-07-12 | 西安理工大学 | N×n型采样方式的次谐波海洋湍流相位屏模拟方法 |
| CN117411577A (zh) * | 2023-10-16 | 2024-01-16 | 西安电子科技大学 | 一种水下无线光通信垂直湍流信道特性分析模型的构建方法、系统、设备及介质 |
-
2024
- 2024-04-26 CN CN202410508362.5A patent/CN118094964B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6958812B1 (en) * | 2003-08-19 | 2005-10-25 | The United States Of America As Represented By The Secretary Of The Army | Systems and methods for analyzing particle systems of surface facets using polarized scattered light |
| CN108763610A (zh) * | 2018-03-28 | 2018-11-06 | 华北电力大学（保定） | 一种基于谱反演法的大气湍流相位屏仿真方法 |
| CN111555822A (zh) * | 2020-04-28 | 2020-08-18 | 西安邮电大学 | 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 |
| AU2020102396A4 (en) * | 2020-09-23 | 2020-11-05 | Changchun University Of Science And Technology | Machine learning-based method for analyzing characteristics of laser beam propagation through turbulent atmosphere |
| CN113868847A (zh) * | 2021-09-14 | 2021-12-31 | 西安理工大学 | 一种偏振光在火星沙尘暴中传输的蒙特卡洛计算方法 |
| CN114741903A (zh) * | 2022-05-18 | 2022-07-12 | 西安理工大学 | N×n型采样方式的次谐波海洋湍流相位屏模拟方法 |
| CN117411577A (zh) * | 2023-10-16 | 2024-01-16 | 西安电子科技大学 | 一种水下无线光通信垂直湍流信道特性分析模型的构建方法、系统、设备及介质 |
Non-Patent Citations (2)
| Title |
|---|
| 张建强;翟焱望;付时尧;高春清;: "径向偏振矢量光束在大气湍流下的传输分析", 光学学报, no. 11, 10 June 2020 (2020-06-10) * |
| 李玉杰等: "非Kolmogorov大气湍流随机相位屏模拟", 红外与激光工程, 25 December 2016 (2016-12-25), pages 1 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118094964B (zh) | 2024-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Applying the parabolic equation to tropospheric groundwave propagation: A review of recent achievements and significant milestones | |
| Zaporozhets et al. | Bistatic RCS calculations with the vector parabolic equation method | |
| CN110399680B (zh) | 一种浅海弹性结构辐射声场计算方法 | |
| CN108519602A (zh) | 海底冷泉快速成像系统及其阵列信号处理方法 | |
| CN114779170A (zh) | 一种浅海近场声源定位方法 | |
| CN110046374B (zh) | 一种基于高斯型声束的声场计算方法 | |
| CN108508446A (zh) | 基于海底冷泉成像数据的扇形变换方法 | |
| Huang et al. | A Gaussian finite-element method for description of sound diffraction | |
| Nakamura et al. | Sound pressure fields focused using biconcave acoustic lens for normal incidence | |
| CN110208764A (zh) | 基于电磁散射计算的动态海洋大型场景回波仿真方法 | |
| CN117970291A (zh) | 基于自聚焦光场的大视场低慢小目标激光探测方法及装置 | |
| CN118094964B (zh) | 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 | |
| CN111679248A (zh) | 一种基于海底水平l型阵列的目标方位和距离联合稀疏重构定位方法 | |
| He et al. | Kirchhoff approximations for the forward-scattering target strength of underwater objects | |
| Lamonaca et al. | Non-reflective hard source method for multiple physically extended sources and scattering bodies | |
| CN116125401A (zh) | 基于电磁散射计算的动态海面回波仿真方法 | |
| Kiehbadroudinezhad et al. | The spatial correlation of a multiple-input multiple-output and channel model using huygens-fresnel principle for underwater acoustic | |
| CN107356924A (zh) | 针对凹腔类结构的sar图像仿真方法 | |
| Weinberg | CASS roots | |
| Zou et al. | Passive estimation method for motion parameters of underwater near-field moving target | |
| Xu et al. | Multi-beam Line Measurement Problem Based on Geometric Modeling and Genetic Algorithm | |
| Hu et al. | A method of sound propagation loss calculation based on Gaussian beams | |
| Wang | A Deep-Water Bottom Reverberation Model Based on Ray Theory | |
| CN114791992B (zh) | 一种基于射线理论的深海目标前向散射声场计算方法 | |
| Wu et al. | Beam-forming and imaging using acoustic lenses: Some simulation and experimental results |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |