CN111555822A - 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 - Google Patents

一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 Download PDF

Info

Publication number
CN111555822A
CN111555822A CN202010350908.0A CN202010350908A CN111555822A CN 111555822 A CN111555822 A CN 111555822A CN 202010350908 A CN202010350908 A CN 202010350908A CN 111555822 A CN111555822 A CN 111555822A
Authority
CN
China
Prior art keywords
photon
photons
phase screen
phase
scattering
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
Application number
CN202010350908.0A
Other languages
English (en)
Other versions
CN111555822B (zh
Inventor
张建磊
杨祎
寇琳琳
贺锋涛
段作梁
陆蓉
王烨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Growerm Photoelectric Technology Co ltd
Original Assignee
Xian University of Posts and Telecommunications
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian University of Posts and Telecommunications filed Critical Xian University of Posts and Telecommunications
Priority to CN202010350908.0A priority Critical patent/CN111555822B/zh
Publication of CN111555822A publication Critical patent/CN111555822A/zh
Application granted granted Critical
Publication of CN111555822B publication Critical patent/CN111555822B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/70Photonic quantum communication
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B13/00Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
    • H04B13/02Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/0082Monitoring; Testing using service channels; using auxiliary channels
    • H04B17/0087Monitoring; Testing using service channels; using auxiliary channels using auxiliary channels or channel simulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

一种基于相位屏的水下无线光传输蒙特卡洛仿真方法,它是由下述步骤组成:A、仿真系统初始化;B、对高斯光源进行采样,生成光子的初始坐标与方向;C、选定海洋功率谱,根据谱反演方法,生成仿真所用相位屏;D、计算光子随机步长;E、根据光子坐标,方向与随机步长,判断光子是否与相位屏相交,是否达到接收面;F、若光子不与相位屏相交,未到达接收面,则更新光子的权重、位置、方向;G、若光子与相位屏相交,未到达接收面,则计算光子与相位屏相交的坐标,更新光子坐标、方向;E、循环执行步骤D~G,直到光子与接收面相交;更加全面的研究光信号在海水中传输的多径效应,可推广应用到水下无线光通信领域。

Description

一种基于相位屏的水下无线光传输蒙特卡洛仿真方法
技术领域
本发明属于水下无线光通信领域,尤其涉及一种基于相位屏的水下无线光传输蒙特卡洛仿真方法。
背景技术
水下通信是探索海洋的关键技术之一。由于电磁波在水下衰减严重,声学通信是目前水下应用最广泛的技术。水下声学通信虽然传输距离远,但受到传输带宽限制,通信速率限制在Kbps量级。研究发现,蓝绿波段可见光可作为可见光通信窗口。近年来,蓝绿光水下光通信由于速率高、方向性好受到研究者广泛关注。
海水介质成分复杂,含有叶绿素、黄色物质及各种矿物质。另外,由于水下盐度场、温度场及水流密度场的存在导致海洋湍流现象。复杂的海洋环境将对蓝绿信号光传输造成严重降质。水中杂质对传输光束具有散射作用,这是水下光脉冲展宽及码间串扰的主要原因。利用散射相函数可以模拟水中粒子的散射特性。进一步借助蒙特卡洛方法可以进行水下无线光信道模拟,可以得到大量光子在水下传输时由于散射引起的时间展宽效应以及光斑空间扩展现象。仅由粒子散射效应导致的时间展宽效应以及光斑扩展现象是时间不变的。湍流效应不会衰减光信号,但会造成海水折射率的不均匀性,改变光传播方向,造成接收端的光强闪烁,是接收机误码率特性的主要因素之一。湍流效应随着传输距离的增加,变得更加剧烈。尤其对于长距离的水下蓝绿光通信,湍流效应不可忽略。
真实海水环境对传输光信号的影响是光吸收散射、湍流效应的综合作用。但是对于水下无线光信道建模,水下粒子散射效应和湍流效应在信道仿真过程难以统一。水下湍流效应无法融入到模拟吸收散射的蒙特卡洛仿真框架。
解决上述技术问题采用的技术方案是由下述步骤组成:
A、仿真系统初始化;
以Z轴正方向为光传输方向,Z=0为光源所在平面,设置传输距离为10~200m,波长λ为400nm~600nm;光子数N为1000~10000000;高斯光束的束宽ω0为0.5~50mm,Δz为相位屏在Z轴上的间隔,取值范围0.5~10m,相位屏的大小为0.5~3m×0.5~3m;方位角
Figure BDA0002471802030000021
为光子散射方向在XOY面上的投影与X轴正半轴的夹角;散射角θ为光子传播方向与X轴正半轴之间夹角;
设定水体系数a为0.1~3/m,散射系数b为0.01~3/m,衰减系数c=a+b,海水信道的单次散射率Tsca为散射系数b与衰减系数c之比,光子初始权重为1;
B.对高斯光源进行采样,确定光子的初始坐标和初始方向,采样方法如下:
(B1)生成取值范围为(0,1)的随机数r1和r2
(B2)确定初始散射角
Figure BDA0002471802030000022
(B3)确定方位角
Figure BDA0002471802030000023
(B4)确定径向距离
Figure BDA0002471802030000024
(B5)确定光子初始坐标(x0,y0,z0)为
Figure BDA0002471802030000025
(B6)确定光子初始方向(μxyz)为
Figure BDA0002471802030000026
(B7)重复B1~B6,即可生成所有光子的初始坐标与初始方向;
C、生成相位屏,具体步骤如下:
(C1)选定海洋功率谱φ(κxy):
Figure BDA0002471802030000027
其中,κx和κy分别是空间频率在x方向和y方向上的分量,ε为单位体积海水的能动耗散功率;χT为均方温度耗散率,ω为温度导致的海洋湍流与盐度导致的海洋湍流的比值,η为Kolmogorov尺度;其中,
Figure BDA0002471802030000031
AT=1.863×10-2,As=1.9×10-4,ATS=9.41×10-3
(C2)确定相位屏谱
Figure BDA0002471802030000032
Figure BDA0002471802030000033
其中,k=2π/λ为光束的波数,Δz为相位屏在Z轴上的间隔;
(C3)利用
Figure BDA0002471802030000034
对高斯随机复矩阵h(κxy)进行滤波和傅里叶变换,即可得到相位屏
Figure BDA0002471802030000035
Figure BDA0002471802030000036
其中,C为控制相位屏方差的常数因子,
Figure BDA0002471802030000037
其中,Δκx、Δκy为取样间隔;
(C4)重复C1~C3,即可生成全部相位屏;
D、计算光子的随机步长d:
Figure BDA0002471802030000038
其中,r3是范围为(0,1)的随机数;
E、根据光子坐标(x,y,z)、光子方向(μxyz)与随机步长d,判断光子是否与相位屏相交,是否达到接收面;
F、若光子不与相位屏相交,未到达接收面,则
(F1)更新光子的权重Wi+1
Wi+1=Wi·Tsca
其中,Wi+1为散射i+1次后光子的权重,Wi为前i次光子权重损失之和,Tsca为海水信道的单次散射率;
(F2)根据HG散射相函数计算散射角θ与方位角
Figure BDA0002471802030000039
Figure BDA0002471802030000041
其中,g为非对称因子,取g=0.924,r4和r5是取值范围为(0,1)的随机数;
(F3)更新光子的方向
Figure BDA0002471802030000042
Figure BDA0002471802030000043
Figure BDA0002471802030000044
其中,
Figure BDA0002471802030000045
(F4)更新光子的坐标(xi+1,yi+1,zi+1):
Figure BDA0002471802030000046
G、若光子与相位屏相交,未到达接收面,则
(G1)计算光子与相位屏相交的坐标(xp,yp,xp)为:
Figure BDA0002471802030000047
其中,(xi,yi,zi)和(μxiyizi)分别为穿过相位屏前光子的坐标与方向,zphase是相位屏在Z轴上的位置;
(G2)若光子在穿过相位屏时,传输方向发生改变,根据广义折射定律,更新光子的方向(mxap,myap,mzap):
Figure BDA0002471802030000051
(G3)光子通过相位屏,下一次散射发生前光子的坐标(xi+1,yi+1,zi+1)为:
Figure BDA0002471802030000052
(G4)根据HG散射相函数计算散射角θ与方位角
Figure BDA0002471802030000053
Figure BDA0002471802030000054
其中,g为非对称因子,取g=0.924,r4和r5是取值范围为(0,1)的随机数;
(G5)更新光子的方向
Figure BDA0002471802030000055
Figure BDA0002471802030000056
Figure BDA0002471802030000061
其中,
Figure BDA0002471802030000062
E、循环执行步骤D~G,直到光子与接收面相交,记录光子与接收面相交的位置坐标,记录光子的权重。
本发明上述步骤A中高斯光束的束宽为ω0=20mm,波长λ=532nm;光子数为100000;在Z轴上以Δz=2.5m为间隔,设置7张相位屏,位置分别为z=2.5m,z=5m,z=7.5m,z=10m,z=12.5m,z=15m,z=17.5m;相位屏大小为1m×1m。
本发明上述的步骤(C1)中
χT=10-7K2/s,ε=10-4m2/s3,η=5mm,ω=-0.3。
与本发明相比,现有技术仅能对单纯的水下散射效应仿真方法或者水下湍流效应仿真方法,针对散射效应和湍流效应所采用的仿真方法不同,故无法对两种效应进行同时仿真模拟,本发明能够同时对水下散射效应与湍流效应进行仿真,模拟产生的接收光斑更加真实,能够为水下无线光通信系统信道分析提供支撑。
附图说明
图1为本发明的原理图,图中所展示为光子在水下迁移时,与水中粒子进行碰撞产生吸收散射,以及通过相位屏时传输方向改变的示意图。其中,1、高斯光源;2、相位屏;3、接收平面。
图2为本发明的整体流程图。
具体实施方式
下面结合附图和实施例对本发明做进一步详细说明,但本发明不限于这些实施例。
实施例1
在图1、2中,本发明一种基于相位屏的水下无线光传输蒙特卡洛仿真方法,是由下述步骤组成:
A、仿真系统初始化;
以Z轴正方向为光传输方向,Z=0为光源所在平面,设置传输距离为20m,高斯光束的束宽为ω0=20mm,波长λ=532nm;光子数为100000;在Z轴上以Δz=2.5m为间隔,设置7张相位屏,位置分别为z=2.5m,z=5m,z=7.5m,z=10m,z=12.5m,z=15m,z=17.5m;相位屏大小为1m×1m。
方位角
Figure BDA0002471802030000071
为光子散射方向在XOY面上的投影与X轴正半轴的夹角;散射角θ为光子传播方向与X轴正半轴之间夹角;
本实施例中,水体系数为a=0.069/m,散射系数b=0.08/m,衰减系数c=a+b=0.149/m,海水信道的单次散射率Tsca为散射系数b与衰减系数c之比,光子初始权重为1;
B.对高斯光源进行采样,确定光子的初始坐标和初始方向,采样方法如下:
(B1)生成取值范围为(0,1)的随机数r1和r2
(B2)确定初始散射角
Figure BDA0002471802030000072
(B3)确定方位角
Figure BDA0002471802030000073
(B4)确定径向距离
Figure BDA0002471802030000074
(B5)确定光子初始坐标(x0,y0,z0)为
Figure BDA0002471802030000075
(B6)确定光子初始方向(μxyz)为
Figure BDA0002471802030000076
(B7)重复B1~B6,即可生成所有光子的初始坐标与初始方向;
C、生成相位屏,具体步骤如下:
(C1)选定海洋功率谱φ(κxy):
Figure BDA0002471802030000081
其中,κx和κy分别是空间频率在x方向和y方向上的分量,ε为单位体积海水的能动耗散功率,取值范围为10-10~10-1m2/s3;χT为均方温度耗散率,取值范围为10-10~10-4K2/s,ω为温度导致的海洋湍流与盐度导致的海洋湍流的比值,取值范围为(-5<ω<0),η为Kolmogorov尺度,取值范围6×10-5~0.01m;
Figure BDA0002471802030000082
AT=1.863×10-2,As=1.9×10-4,ATS=9.41×10-3;本实施例中,χT=10-7K2/s,ε=10-4m2/s3,η=5mm,ω=-0.3。
(C2)确定相位屏谱
Figure BDA0002471802030000083
Figure BDA0002471802030000084
其中,k=2π/λ为光束的波数,Δz为相位屏在Z轴上的间隔;
(C3)利用
Figure BDA0002471802030000085
对高斯随机复矩阵h(κxy)进行滤波和傅里叶变换,即可得到相位屏:
Figure BDA0002471802030000086
其中,C为控制相位屏方差的常数因子,
Figure BDA0002471802030000087
对上述相位屏表达式进行离散化,在空间域内,x=mΔx,y=nΔy,Δx=Δy=1mm为采样间隔,m、n为整数;在波数域内,κx=m′Δκx,κy=n′Δκy,m′、n′为整数,Δκx、Δκy为取样间隔,Δκx=2π/(NΔx),Δκy=2π/(NΔy),N为1000。
Figure BDA0002471802030000088
可得相位屏表达式
Figure BDA0002471802030000089
Figure BDA00024718020300000810
(C4)重复C1~C3,即可生成全部相位屏;
D、计算光子的随机步长d:
Figure BDA0002471802030000091
其中,r3是范围为(0,1)的随机数;
E、根据光子坐标(x,y,z)、光子方向(μxyz)与随机步长d,判断光子是否与相位屏相交,是否达到接收面;
F、若光子不与相位屏相交,未到达接收面,则
(F1)更新光子的权重Wi+1
Wi+1=Wi·Tsca
其中,Wi+1为散射i+1次后光子的权重,Wi为前i次光子权重损失之和,Tsca为海水信道的单次散射率;
(F2)根据HG散射相函数计算散射角θ与方位角
Figure BDA0002471802030000092
Figure BDA0002471802030000093
其中,g为非对称因子,取g=0.924,r4和r5是取值范围为(0,1)的随机数;
(F3)更新光子的方向
Figure BDA0002471802030000094
Figure BDA0002471802030000095
Figure BDA0002471802030000096
其中,
Figure BDA0002471802030000101
(F4)更新光子的坐标(xi+1,yi+1,zi+1):
Figure BDA0002471802030000102
G、若光子与相位屏相交,未到达接收面,则
(G1)计算光子与相位屏相交的坐标(xp,yp,xp)为:
Figure BDA0002471802030000103
其中,(xi,yi,zi)和(μxiyizi)分别为穿过相位屏前光子的坐标与方向,zphase是相位屏在Z轴上的位置;
(G2)若光子在穿过相位屏时,传输方向发生改变,根据广义折射定律,更新光子的方向(mxap,myap,mzap):
Figure BDA0002471802030000104
(G3)光子通过相位屏,下一次散射发生前光子的坐标(xi+1,yi+1,zi+1)为:
Figure BDA0002471802030000111
(G4)根据HG散射相函数计算散射角θ与方位角
Figure BDA0002471802030000112
Figure BDA0002471802030000113
其中,g为非对称因子,取g=0.924,r4和r5是取值范围为(0,1)的随机数;
(G5)更新光子的方向
Figure BDA0002471802030000114
Figure BDA0002471802030000115
Figure BDA0002471802030000116
其中,
Figure BDA0002471802030000117
E、循环执行步骤D~G,直到光子与接收面相交,记录光子与接收面相交的位置坐标,记录光子的权重。
实施例2
本实施例的步骤A中,以Z轴正方向为光传输方向,Z=0为光源所在平面,设置传输距离为10m,高斯光束的束宽为ω0=0.5mm,波长λ=450nm;光子数为1000;在Z轴上以Δz=0.5m为间隔,设置5张相位屏,位置分别为z=0.5m,z=1m,z=1.5m,z=2m,z=2.5m;相位屏大小为0.5m×0.5m。
其余各步骤与实施例1完全相同。
实施例3
本实施例的步骤A中,以Z轴正方向为光传输方向,Z=0为光源所在平面,设置传输距离为200m,高斯光束的束宽为ω0=50mm,波长λ=589nm;光子数为10000000;在Z轴上以Δz=10m为间隔,设置20张相位屏,位置分别为z=10m,z=120m,z=30m,z=40m,z=50m,z=60m,z=70m,z=80m,z=90m,z=100m,z=110m,z=120m,z=130m,z=140m,z=150m,z=160m,z=170m,z=180m,z=190m,z=200m;相位屏大小为3m×3m。
其余各步骤与实施例1完全相同。

Claims (3)

1.一种基于相位屏的水下无线光传输蒙特卡洛仿真方法,其特征在于它是由下述步骤组成:
A、仿真系统初始化;
以Z轴正方向为光传输方向,Z=0为光源所在平面,设置传输距离为10~200m,波长λ为400nm~600nm;光子数N为1000~10000000;高斯光束的束宽ω0为0.5~50mm,Δz为相位屏在Z轴上的间隔,取值范围0.5~10m,相位屏的大小为0.5~3m×0.5~3m;方位角
Figure FDA0002471802020000011
为光子散射方向在XOY面上的投影与X轴正半轴的夹角;散射角θ为光子传播方向与X轴正半轴之间夹角;
设定水体系数a为0.1~3/m,散射系数b为0.01~3/m,衰减系数c=a+b,海水信道的单次散射率Tsca为散射系数b与衰减系数c之比,光子初始权重为1;
B.对高斯光源进行采样,确定光子的初始坐标和初始方向,采样方法如下:
(B1)生成取值范围为(0,1)的随机数r1和r2
(B2)确定初始散射角
Figure FDA0002471802020000012
(B3)确定方位角
Figure FDA0002471802020000013
(B4)确定径向距离
Figure FDA0002471802020000014
(B5)确定光子初始坐标(x0,y0,z0)为
Figure FDA0002471802020000015
(B6)确定光子初始方向(μxyz)为
Figure FDA0002471802020000016
(B7)重复B1~B6,即可生成所有光子的初始坐标与初始方向;
C、生成相位屏,具体步骤如下:
(C1)选定海洋功率谱φ(κxy):
Figure FDA0002471802020000021
其中,κx和κy分别是空间频率在x方向和y方向上的分量,ε为单位体积海水的能动耗散功率;χT为均方温度耗散率,ω为温度导致的海洋湍流与盐度导致的海洋湍流的比值,η为Kolmogorov尺度;其中,
Figure FDA0002471802020000022
AT=1.863×10-2,As=1.9×10-4,ATS=9.41×10-3
(C2)确定相位屏谱
Figure FDA0002471802020000023
Figure FDA0002471802020000024
其中,k=2π/λ为光束的波数,Δz为相位屏在Z轴上的间隔;
(C3)利用
Figure FDA0002471802020000025
对高斯随机复矩阵h(κxy)进行滤波和傅里叶变换,即可得到相位屏
Figure FDA0002471802020000026
Figure FDA0002471802020000027
其中,C为控制相位屏方差的常数因子,
Figure FDA0002471802020000028
其中,Δκx、Δκy为取样间隔;
(C4)重复C1~C3,即可生成全部相位屏;
D、计算光子的随机步长d:
Figure FDA0002471802020000029
其中,r3是范围为(0,1)的随机数;
E、根据光子坐标(x,y,z)、光子方向(μxyz)与随机步长d,判断光子是否与相位屏相交,是否达到接收面;
F、若光子不与相位屏相交,未到达接收面,则
(F1)更新光子的权重Wi+1
Wi+1=Wi·Tsca
其中,Wi+1为散射i+1次后光子的权重,Wi为前i次光子权重损失之和,Tsca为海水信道的单次散射率;
(F2)根据HG散射相函数计算散射角θ与方位角
Figure FDA0002471802020000031
Figure FDA0002471802020000032
其中,g为非对称因子,取g=0.924,r4和r5是取值范围为(0,1)的随机数;
(F3)更新光子的方向
Figure FDA0002471802020000033
Figure FDA0002471802020000034
Figure FDA0002471802020000035
其中,
Figure FDA0002471802020000036
(F4)更新光子的坐标(xi+1,yi+1,zi+1):
Figure FDA0002471802020000037
G、若光子与相位屏相交,未到达接收面,则
(G1)计算光子与相位屏相交的坐标(xp,yp,xp)为:
Figure FDA0002471802020000041
其中,(xi,yi,zi)和(μxiyizi)分别为穿过相位屏前光子的坐标与方向,zphase是相位屏在Z轴上的位置;
(G2)若光子在穿过相位屏时,传输方向发生改变,根据广义折射定律,更新光子的方向(mxap,myap,mzap):
Figure FDA0002471802020000042
(G3)光子通过相位屏,下一次散射发生前光子的坐标(xi+1,yi+1,zi+1)为:
Figure FDA0002471802020000043
(G4)根据HG散射相函数计算散射角θ与方位角
Figure FDA0002471802020000044
Figure FDA0002471802020000045
其中,g为非对称因子,取g=0.924,r4和r5是取值范围为(0,1)的随机数;
(G5)更新光子的方向
Figure FDA0002471802020000051
Figure FDA0002471802020000052
Figure FDA0002471802020000053
其中,
Figure FDA0002471802020000054
E、循环执行步骤D~G,直到光子与接收面相交,记录光子与接收面相交的位置坐标,记录光子的权重。
2.根据权利要求1所述的一种基于相位屏的水下无线光传输蒙特卡洛仿真方法,其特征在于:上述步骤A中高斯光束的束宽为ω0=20mm,波长λ=532nm;光子数为100000;在Z轴上以Δz=2.5m为间隔,设置7张相位屏,位置分别为z=2.5m,z=5m,z=7.5m,z=10m,z=12.5m,z=15m,z=17.5m;相位屏大小为1m×1m。
3.根据权利要求1所述的一种基于相位屏的水下无线光传输蒙特卡洛仿真方法,其特征在于:上述的步骤(C1)中
χT=10-7K2/s,ε=10-4m2/s3,η=5mm,ω=-0.3。
CN202010350908.0A 2020-04-28 2020-04-28 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 Active CN111555822B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010350908.0A CN111555822B (zh) 2020-04-28 2020-04-28 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010350908.0A CN111555822B (zh) 2020-04-28 2020-04-28 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法

Publications (2)

Publication Number Publication Date
CN111555822A true CN111555822A (zh) 2020-08-18
CN111555822B CN111555822B (zh) 2021-08-20

Family

ID=72008242

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010350908.0A Active CN111555822B (zh) 2020-04-28 2020-04-28 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法

Country Status (1)

Country Link
CN (1) CN111555822B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112583480A (zh) * 2020-12-24 2021-03-30 微感科技(南通)有限公司 一种高精确度无线光通信链路的实现方法
CN112636821A (zh) * 2020-12-18 2021-04-09 南京先进激光技术研究院 水下无线光通信光信道仿真方法及仿真系统
CN113078957A (zh) * 2021-03-16 2021-07-06 河海大学 一种基于空间坐标变换技术的uwoc系统信道特征mcns方法
CN113128034A (zh) * 2021-04-07 2021-07-16 西安邮电大学 一种基于gpu的水下无线光信道并行仿真方法
CN113297643A (zh) * 2021-04-14 2021-08-24 西安理工大学 一种水下蓝绿激光通信中光斑面积的计算方法
CN114578376A (zh) * 2022-05-05 2022-06-03 中国科学院西安光学精密机械研究所 一种基于海洋湍流的单光子成像仿真方法
CN118094964A (zh) * 2024-04-26 2024-05-28 山东省科学院海洋仪器仪表研究所 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7069068B1 (en) * 1999-03-26 2006-06-27 Oestergaard Leif Method for determining haemodynamic indices by use of tomographic data
US20160334793A1 (en) * 2015-04-09 2016-11-17 University Of New Hampshire POSE DETECTION AND CONTROL OF UNMANNED UNDERWATER VEHICLES (UUVs) UTILIZING AN OPTICAL DETECTOR ARRAY
CN107018105A (zh) * 2017-04-18 2017-08-04 东莞信大融合创新研究院 水下可见光通信方法及系统
CN108023652A (zh) * 2017-10-27 2018-05-11 西安邮电大学 一种应用于海水信道的激光传输特性的模拟方法
US20190007137A1 (en) * 2017-04-27 2019-01-03 Nec Laboratories America, Inc. Oam based physical layer security using hybrid free-space optical-terahertz technology
US20190129312A1 (en) * 2017-10-26 2019-05-02 Government Of The United States, As Represented By The Secretary Of The Air Force Generation of Vector Partially Coherent Optical Sources Using Phase-Only Spatial Light Modulators
CN109995427A (zh) * 2019-03-25 2019-07-09 西安电子科技大学 一种水下上行激光通信的蒙特卡洛仿真方法
US20190266712A1 (en) * 2018-02-24 2019-08-29 United States Of America As Represented By The Administrator Of The Nasa System and method for imaging underwater environments using fluid lensing
CN111027231A (zh) * 2019-12-29 2020-04-17 杭州科洛码光电科技有限公司 水下阵列相机成像方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7069068B1 (en) * 1999-03-26 2006-06-27 Oestergaard Leif Method for determining haemodynamic indices by use of tomographic data
US20160334793A1 (en) * 2015-04-09 2016-11-17 University Of New Hampshire POSE DETECTION AND CONTROL OF UNMANNED UNDERWATER VEHICLES (UUVs) UTILIZING AN OPTICAL DETECTOR ARRAY
CN107018105A (zh) * 2017-04-18 2017-08-04 东莞信大融合创新研究院 水下可见光通信方法及系统
US20190007137A1 (en) * 2017-04-27 2019-01-03 Nec Laboratories America, Inc. Oam based physical layer security using hybrid free-space optical-terahertz technology
US20190129312A1 (en) * 2017-10-26 2019-05-02 Government Of The United States, As Represented By The Secretary Of The Air Force Generation of Vector Partially Coherent Optical Sources Using Phase-Only Spatial Light Modulators
CN108023652A (zh) * 2017-10-27 2018-05-11 西安邮电大学 一种应用于海水信道的激光传输特性的模拟方法
US20190266712A1 (en) * 2018-02-24 2019-08-29 United States Of America As Represented By The Administrator Of The Nasa System and method for imaging underwater environments using fluid lensing
CN109995427A (zh) * 2019-03-25 2019-07-09 西安电子科技大学 一种水下上行激光通信的蒙特卡洛仿真方法
CN111027231A (zh) * 2019-12-29 2020-04-17 杭州科洛码光电科技有限公司 水下阵列相机成像方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
FAN WANG ET AL: "Light transmission characteristic analyses of a", 《APPLIED OPTICS》 *
JIANLEI ZHANG ET AL: "Monte-Carlo-based optical wireless underwater channel modeling with", 《OPTICS COMMUNICATIONS》 *
RENZHI YUAN ET AL: "Monte-Carlo Integration Models for Multiple", 《IEEE TRANSACTIONS ON COMMUNICATIONS》 *
ZAHRA VALI ET AL: "Modeling turbulence in underwater wireless", 《JOURNAL OF THE OPTICAL SOCIETY OF AMERICA A》 *
杨天星等人: "海洋湍流随机相位屏模型", 《光学学报》 *
杨曈等人: "基于蒙特卡罗法的水下激光光幕探测性能研究", 《应用光学》 *
牛超君等人: "相位屏法模拟高斯阵列光束海洋湍流传输特性", 《光学学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112636821A (zh) * 2020-12-18 2021-04-09 南京先进激光技术研究院 水下无线光通信光信道仿真方法及仿真系统
CN112583480A (zh) * 2020-12-24 2021-03-30 微感科技(南通)有限公司 一种高精确度无线光通信链路的实现方法
CN113078957A (zh) * 2021-03-16 2021-07-06 河海大学 一种基于空间坐标变换技术的uwoc系统信道特征mcns方法
CN113078957B (zh) * 2021-03-16 2022-06-21 河海大学 一种基于空间坐标变换技术的uwoc系统信道特征mcns方法
CN113128034A (zh) * 2021-04-07 2021-07-16 西安邮电大学 一种基于gpu的水下无线光信道并行仿真方法
CN113297643A (zh) * 2021-04-14 2021-08-24 西安理工大学 一种水下蓝绿激光通信中光斑面积的计算方法
CN113297643B (zh) * 2021-04-14 2023-02-03 西安理工大学 一种水下蓝绿激光通信中光斑面积的计算方法
CN114578376A (zh) * 2022-05-05 2022-06-03 中国科学院西安光学精密机械研究所 一种基于海洋湍流的单光子成像仿真方法
CN118094964A (zh) * 2024-04-26 2024-05-28 山东省科学院海洋仪器仪表研究所 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法
CN118094964B (zh) * 2024-04-26 2024-07-23 山东省科学院海洋仪器仪表研究所 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法

Also Published As

Publication number Publication date
CN111555822B (zh) 2021-08-20

Similar Documents

Publication Publication Date Title
CN111555822B (zh) 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法
CN109327271B (zh) 基于部分相干光的水下光通信系统平均信道容量确定方法
Liu et al. Special Smith–Purcell radiation from an open resonator array
CN114884599B (zh) 蓝绿涡旋光气海跨介质准直入射单程传输模型的构建方法
Fu et al. Effect of aperture averaging on mean bit error rate for UWOC system over moderate to strong oceanic turbulence
Angara et al. Underwater wireless optical communication system channel modelling with oceanic bubbles and water constituents under different wind conditions
Lou et al. Basic principles of underwater acoustic communication
CN106644104A (zh) 一种基于谱反演法的离散雨滴介质的相位屏建模方法
CN114578376B (zh) 一种基于海洋湍流的单光子成像仿真方法
CN116054955A (zh) 一种基于相干检测的涡旋光模式识别装置、方法及程序
Ma et al. Research on channel modeling method of underwater optical communication
CN114204985A (zh) 无线紫外光非直视通信中光子检测概率快速估算方法
CN114884571A (zh) 一种幂指数型Airy光束轨道角动量模式在等离子体湍流中传输的计算方法
Aref et al. Design and simulation of a new model for shallow water multipath acoustic channel in the Persian Gulf
Cheng et al. Influence of the complex refractive index of underwater suspended particles on the transmission characteristics of polarized light
CN118094964B (zh) 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法
CN114745046B (zh) 一种分析从随机波动海面出射激光光束指向偏差的方法
CN116663377A (zh) 基于谱反演法的离散海洋悬浮粒子介质的相位屏建模方法
Yuksel et al. A geometrical optics approach for modeling atmospheric turbulence
CN116073926A (zh) 一种包含气泡散射的水下无线光通信信道建模方法及系统
CN114791992B (zh) 一种基于射线理论的深海目标前向散射声场计算方法
CN117938286A (zh) 蓝绿激光大气-海洋跨介质传输模型的构建方法
Ferlic et al. Propagation of Laguerre-Gaussian beams through underwater optical turbulence
Li et al. Pulse response of underwater optical wireless channel based on TTHG model
CN116388867A (zh) 部分相干反常涡旋光束在各向异性湍流大气中传输的仿真方法

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230625

Address after: 710000 Room 023, F2003, 20th Floor, Building 4-A, Xixian Financial Port, Fengdong New City Energy Jinmao District, Xixian New District, Xi'an City, Shaanxi Province

Patentee after: Shaanxi Growerm Photoelectric Technology Co.,Ltd.

Address before: 710121 South Campus of Xi'an University of Posts and telecommunications, Weiguo Road, Chang'an District, Xi'an City, Shaanxi Province

Patentee before: XI'AN University OF POSTS & TELECOMMUNICATIONS