CN110661568B - 水下到达激光信号3dB强度光斑半径的计算方法 - Google Patents
水下到达激光信号3dB强度光斑半径的计算方法 Download PDFInfo
- Publication number
- CN110661568B CN110661568B CN201911031117.5A CN201911031117A CN110661568B CN 110661568 B CN110661568 B CN 110661568B CN 201911031117 A CN201911031117 A CN 201911031117A CN 110661568 B CN110661568 B CN 110661568B
- Authority
- CN
- China
- Prior art keywords
- photon
- scattering
- receiving plane
- arrival
- power value
- 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Optical Communication System (AREA)
Abstract
本发明公开一种水下到达激光信号3dB强度光斑半径的计算方法,通过建立水下激光信号传输过程的MC仿真平台,在给定信道条件下,统计接收平面上到达光子的坐标信息和功率值信息,便以此为依据进行光子二维分布情况统计,由此得到到达激光信号的光斑。本发明能够为海洋激光通信系统中光学接收镜头(含接收镜头阵列)的尺寸设计及其最优位置分布策略提供数据和技术参考。
Description
技术领域
本发明涉及水下激光通信技术领域,具体涉及一种水下到达激光信号3dB强度光斑半径的计算方法。
背景技术
激光在散射海水信道的传输过程,可以视为大量光子在其中的传播过程。在水下激光信道研究领域中,由于海水复杂而随机的时空变特性,海水中水分子、浮游植物和有色溶解有机物等的吸收效应以及水分子、浮游植物和悬浮颗粒物等散射效应,激光信号在海水中传播时,会遭受严重的衰减和偏离原有传播路径。
由于海水散射效应的影响,到达激光信号会在接收平面上弥散成比其自然发散大很多的信号光斑,但海水信道下到达激光信号的主要光功率却分布在3dB强度光斑以内,此时若将全部信号光斑分布全部纳入统计范围,计算复杂度必然骤增。同时,实际光学接收镜头的尺寸也远小于到达激光信号光斑的尺寸。因此在给定合理尺寸接收平面的条件下,计算该3dB强度光斑的半径,既可以表征信道和光学系统参数对其产生的影响,又可以为海洋激光通信系统光学接收镜头的尺寸设计及其最优位置分布策略提供数据和技术参考。
虽然目前国内许多研究团队对水下激光信道的特性开展了相关应用研究。但现有研究资料中,仍然没有全面表征出到达激光信号光斑的二维分布特性的相关研究。
发明内容
本发明针对光子到达接收平面的二维空间分布,提供一种水下到达激光信号3dB强度光斑半径的计算方法。
为解决上述问题,本发明是通过以下技术方案实现的:
水下到达激光信号3dB强度光斑半径的计算方法,其具体包括步骤如下:
步骤1、将激光通信系统发收天线的准直光轴视为z轴,以实际过程中光子到达平面即接收平面xOy建立xyz坐标系,并在给定信道条件下,基于MC仿真方法对海水中激光信号的传输过程进行仿真,以MC仿真过程中各个光子到达平面即统计平面x′O′y′建立x′y′z坐标系;
步骤2、对于每个光子,在MC仿真过程中,记录该光子在光源处出的发射位置并记录该光子在运动过程中的每一次散射的到达坐标(x′i,y′i,z′i)及其对应的单位方向向量当该光子的每一次散射的散射步长li在光轴上投影的矢量之和即散射步长矢量和Lz大于等于传输距离Z0时,表明该光子到达统计平面x′O′y′;
步骤3、基于步骤2所记录的仿真结果,计算每个光子到达接收平面xOy的到达坐标(xM,yM,Z0),其中:
式中,为光子的发射位置,为光子第i次散射的单位方向向量,为光子到达统计平面x′O′y′的单位方向向量,li为第i次散射的散射步长,Lz为散射步长矢量和,Z0传输距离,i=1,2,…,M,M为光子到达接收平面时所经历的散射次数;
步骤4、计算每个光子到达接收平面xOy的功率值ω,其中:
ω=ω0uM
式中,ω0为光子发射时所携带的功率值,u=Ks(λ)/Katt(λ)为反照率,Ks(λ)为散射系数,Katt(λ)为衰减系数;M为光子到达接收平面时所经历的散射次数;
步骤5、将接收平面xOy的几何中心与激光通信系统发收天线的准直光轴重合,并将接收平面xOy分割成N×N个正方形的面积微元,其中N为能被2整除的正整数;
步骤6、基于MC仿真过程所得到的每个光子到达接收平面xOy的到达坐标(xM,yM,Z0)和功率值ω,统计每个面积微元内所有光子的功率值RInty(l,m),并将所有面积微元内所有光子的功率值RInty(l,m)进行累加后得到激光信号到达接收平面上的总功率值Rtotal;其中l,m=-N/2,1-N/2,…,N/2-1,N/2;
步骤7、以接收平面xOy的几何中心为原点,以当前迭代半径r3dB(k)作圆,并统计圆内所有光子的功率值之和RInty(k);其中
r3dB(k)=(2RPD/N)×k
式中,2RPD为接收平面xOy的边长,k为迭代次数,k=1,2,…,且k≤N/2;
步骤8、求解满足10lg[RInty(k)/Rtotal]≥-3的最小k值,则当前迭代半径r3dB(k)即为所求的水下到达激光信号3dB强度光斑半径。
与现有技术相比,本发明通过建立水下激光信号传输过程的MC仿真平台,在给定信道条件下,统计接收平面上到达光子的坐标信息和功率值信息,便以此为依据进行光子二维分布情况统计,由此得到到达激光信号的光斑。本发明能够为海洋激光通信系统中光学接收镜头(含接收镜头阵列)的尺寸设计及其最优位置分布策略提供数据和技术参考。
附图说明
图1为光子到达坐标示意图。
图2为3dB强度光斑半径示意图示意图。
图3为纯净海水3dB强度光斑半径与传输距离曲线图。
图4为清洁海水3dB强度光斑半径与传输距离曲线图。
图5为沿岸海水3dB强度光斑半径与传输距离曲线图。
图6为港口海水3dB强度光斑半径与传输距离曲线图。
图7为纯净海水3dB强度光斑半径与接收平面半孔径曲线图。
图8为清洁海水3dB强度光斑半径与接收平面半孔径曲线图。
图9为海岸海水3dB强度光斑半径与接收平面半孔径曲线图。
图10为港口海水3dB强度光斑半径与接收平面半孔径曲线图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实例,并参照附图,对本发明进一步详细说明。
一种水下到达激光信号3dB强度光斑半径的计算方法,其具体如下步骤:
步骤1、如附图1所示,将激光通信系统发收天线的准直光轴视为z轴,以实际过程中光子到达平面即接收平面xOy建立xyz坐标系,并在给定信道条件下,基于MC仿真方法对海水中激光信号的传输过程进行仿真,以MC仿真过程中各个光子到达平面即统计平面x′O′y′建立x′y′z坐标系。
MC仿真方法是研究散射媒介中光波辐射传输方程的通用方法,也常用来追踪光子的传播轨迹。相对于实验方法,它能灵活调整激光信号的各种参数和需要统计大量的光子或光子包,从而可以全面准确的揭示水下激光信道的散射特征。因此,采用MC仿真方法对海水中激光信号的传输过程进行仿真。
步骤2、对于每个光子,在MC仿真过程中,记录该光子在光源处出的发射位置并记录该光子在运动过程中的每一次散射的到达坐标(x′i,y′i,z′i)及其对应的单位方向向量当该光子的每一次散射的散射步长li在光轴上投影的矢量之和即散射步长矢量和Lz大于等于传输距离Z0时,表明该光子到达统计平面x′O′y′,此时光子的第M次散射完成,到达坐标为(x′M,y′M,z′M),单位方向向量为
步骤3、基于步骤2所记录的仿真结果,计算每个光子到达接收平面xOy的到达坐标(xM,yM,Z0),其中:
式中,为光子的发射位置,为光子第i次散射的单位方向向量,为光子到达统计平面x′O′y′的单位方向向量,li为第i次散射的散射步长,Lz为散射步长矢量和,Z0传输距离,i=1,2,…,M,M为光子到达接收平面时所经历的散射次数。
步骤4、计算每个光子到达接收平面xOy的功率值ω,其中:
ω=ω0uM
式中,ω0为光子发射时所携带的功率值,取值为1,u=Ks(λ)/Katt(λ)为反照率,Ks(λ)为散射系数,Katt(λ)为衰减系数。
步骤5、设定接收平面为一个面积有限的正方形,将接收平面xOy的几何中心与激光通信系统发收天线的准直光轴重合,并将接收平面xOy分割成N×N个正方形的面积微元;其中N为一个取值足够大且能被2整除的正整数。
步骤6、以平面的几何中心为原点建立xyz三维坐标系,其中,z轴为光信号的光轴,平面x、y轴方向的半边长为RPD,接收视场角为ΨR。当光子到达功率值ω大于或等于预先设定的阈值Th,到达xOy平面的坐标(xM,yM,Z0)在接收平面尺寸以内且其到达角度小于或等于接收视场半角时,判定光子能够被接收,具体的判断条件如下:
基于MC仿真得到每个光子在接收平面xOy上的到达坐标(xM,yM,Z0)和功率值ω,统计到达每个面积微元内所有光子的功率值和,则第(l,m)个面积微元所覆盖的功率值可以记为RInty(l,m),其中l,m=-N/2,1-N/2,…,N/2-1,N/2,据此,激光信号到达接收平面上的总功率值为:
步骤7、令迭代次数k=1,2,…,且k≤N/2;如附图2所示,以接收平面xOy的几何中心为原点,以当前迭代半径r3dB(k)作圆,并统计圆内所有光子的功率值之和RInty(k);其中
r3dB(k)=(2RPD/N)×k
式中,2RPD为接收平面xOy的边长。
步骤8、求解满足10lg[RInty(k)/Rtotal]≥-3的最小k值,即从1开始不断增加迭代次数k的取值,一但10lg[RInty(k)/Rtotal]≥-3满足时,则当前迭代半径r3dB(k)即为所求的水下到达激光信号3dB强度光斑半径。
为了更好地说明上述方法,表一、表二和表三分别给出不同类型的信道参数和用于MC仿真的系统参数。
表一 信道参数配置表
表二 MC仿真参数配置表(用于研究3dB光斑半径与传输距离的关系)
表三 MC仿真参数配置表(用于研究3dB光斑半径与接收平面尺寸的关系)
根据上述表一和表二的仿真参数,通过表四至表七分别给出了纯净、清洁、沿岸和港口海水信道的若干传输节点处的3dB强度光斑半径大小,对应的曲线关系图由附图3~6给出。
表四 纯净海水信道若干传输节点处3dB强度光斑半径
表五 清洁海水信道若干传输节点处3dB强度光斑半径
表六 沿岸海水信道若干传输节点处3dB强度光斑半径
表七 港口海水信道若干传输节点处3dB强度光斑半径
根据上述表一和表三的仿真参数,通过表八至表十一分别给出了纯净、清洁、沿岸和港口海水信道的若干接收平面半孔径取值下的3dB强度光斑半径大小,对应的曲线关系图由附图7~10给出。
表八 纯净海水信道若干接收平面半孔径节点处的3dB强度光斑半径
表九 清洁海水信道若干接收平面半孔径节点处的3dB强度光斑半径
表十 沿岸海水信道若干接收平面半孔径节点处的3dB强度光斑半径
表十一 港口海水信道若干接收平面半孔径节点处的3dB强度光斑半径
通过表三至表十一及附图3~10的数据和曲线图说明,接收平面到达激光信号的3dB强度光斑半径的大小与信道类型、信号发散角、接收视场角、传输距离和接收平面统计孔径密切相关。故本发明说明了上述各个因素对3dB强度光斑半径的影响情况。
本发明通过建立水下激光信号传输过程的MC仿真平台,在给定信道条件下,统计大小有限的接收平面上到达光子的坐标信息,输出接收平面上光子到达ωi的二维分布情况,提出了一种3dB强度光斑半径计算方法,通过计算给定通信距离和接收平面尺寸条件下该而为分布的3dB强度光斑半径,给出了几种典型信道类型的3dB强度光斑半径与传输距离、接收平面尺寸的曲线关系,为海洋激光通信系统中光学接收镜头的尺寸设计及其最优位置分布策略提供数据和技术参考。
需要说明的是,尽管以上本发明所述的实施例是说明性的,但这并非是对本发明的限制,因此本发明并不局限于上述具体实施方式中。在不脱离本发明原理的情况下,凡是本领域技术人员在本发明的启示下获得的其它实施方式,均视为在本发明的保护之内。
Claims (2)
1.水下到达激光信号3dB强度光斑半径的计算方法,其特征是,其具体包括步骤如下:
步骤1、将激光通信系统发收天线的准直光轴视为z轴,以实际过程中光子到达平面即接收平面xOy建立xyz坐标系,并在给定信道条件下,基于MC仿真方法对海水中激光信号的传输过程进行仿真,以MC仿真过程中各个光子到达平面即统计平面x′O′y′建立x′y′z坐标系;
步骤2、对于每个光子,在MC仿真过程中,记录该光子在光源处的发射位置并记录该光子在运动过程中的每一次散射的到达坐标(xi′,yi′,zi′)及其对应的单位方向向量当该光子的每一次散射的散射步长li在光轴上投影的矢量之和即散射步长矢量和Lz大于等于传输距离Z0时,表明该光子到达统计平面x′O′y′;
步骤3、基于步骤2所记录的仿真结果,计算每个光子到达接收平面xOy的到达坐标(xM,yM,Z0),其中:
式中,为光子的发射位置,为光子第i次散射的单位方向向量,为光子到达统计平面x′O′y′的单位方向向量,li为第i次散射的散射步长,Lz为散射步长矢量和,Z0传输距离,i=1,2,…,M,M为光子到达接收平面时所经历的散射次数;
步骤4、计算每个光子到达接收平面xOy的功率值ω,其中:
ω=ω0uM
式中,ω0为光子发射时所携带的功率值,u=Ks(λ)/Katt(λ)为反照率,Ks(λ)为散射系数,Katt(λ)为衰减系数;M为光子到达接收平面时所经历的散射次数;
步骤5、将接收平面xOy的几何中心与激光通信系统发收天线的准直光轴重合,并将接收平面xOy分割成N×N个正方形的面积微元,其中N为能被2整除的正整数;
步骤6、基于MC仿真过程所得到的每个光子到达接收平面xOy的到达坐标(xM,yM,Z0)和功率值ω,统计每个面积微元内所有光子的功率值RInty(l,m),并将所有面积微元内所有光子的功率值RInty(l,m)进行累加后得到激光信号到达接收平面上的总功率值Rtotal:
其中,RInty(l,m)为第(l,m)个面积微元所覆盖的功率值,l,m=-N/2,1 -N/2,…,N/2-1,N/2;
步骤7、以接收平面xOy的几何中心为原点,以当前迭代半径r3dB(k)作圆,并统计圆内所有光子的功率值之和RInty(k);其中
r3dB(k)=(2RPD/N)×k
式中,2RPD为接收平面xOy的边长,k为迭代次数,k=1,2,…,且k≤N/2;
步骤8、求解满足10lg[RInty(k)/Rtotal]≥-3的最小k值,则当前迭代半径r3dB(k)即为所求的水下到达激光信号3dB强度光斑半径。
2.根据权利要求1所述的水下到达激光信号3dB强度光斑半径的计算方法,其特征是,光子发射时所携带的功率值ω0为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911031117.5A CN110661568B (zh) | 2019-10-28 | 2019-10-28 | 水下到达激光信号3dB强度光斑半径的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911031117.5A CN110661568B (zh) | 2019-10-28 | 2019-10-28 | 水下到达激光信号3dB强度光斑半径的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110661568A CN110661568A (zh) | 2020-01-07 |
CN110661568B true CN110661568B (zh) | 2022-04-15 |
Family
ID=69041885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911031117.5A Active CN110661568B (zh) | 2019-10-28 | 2019-10-28 | 水下到达激光信号3dB强度光斑半径的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110661568B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113297643B (zh) * | 2021-04-14 | 2023-02-03 | 西安理工大学 | 一种水下蓝绿激光通信中光斑面积的计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102594440A (zh) * | 2012-02-22 | 2012-07-18 | 大连大学 | 一种光子传输性能仿真方法 |
CN107018105A (zh) * | 2017-04-18 | 2017-08-04 | 东莞信大融合创新研究院 | 水下可见光通信方法及系统 |
CN108023652A (zh) * | 2017-10-27 | 2018-05-11 | 西安邮电大学 | 一种应用于海水信道的激光传输特性的模拟方法 |
CN109995427A (zh) * | 2019-03-25 | 2019-07-09 | 西安电子科技大学 | 一种水下上行激光通信的蒙特卡洛仿真方法 |
-
2019
- 2019-10-28 CN CN201911031117.5A patent/CN110661568B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102594440A (zh) * | 2012-02-22 | 2012-07-18 | 大连大学 | 一种光子传输性能仿真方法 |
CN107018105A (zh) * | 2017-04-18 | 2017-08-04 | 东莞信大融合创新研究院 | 水下可见光通信方法及系统 |
CN108023652A (zh) * | 2017-10-27 | 2018-05-11 | 西安邮电大学 | 一种应用于海水信道的激光传输特性的模拟方法 |
CN109995427A (zh) * | 2019-03-25 | 2019-07-09 | 西安电子科技大学 | 一种水下上行激光通信的蒙特卡洛仿真方法 |
Non-Patent Citations (4)
Title |
---|
UWOC OR underwater wireless optical communicaiton;Gabriel, C (Gabriel, Chadi);《JOURNAL OF OPTICAL COMMUNICATIONS AND NETWORKING》;20130101;第5卷(第1期);全文 * |
Zhang Yingluo ; Wang Yingmin ; Huang Aiping.Analysis of Underwater Laser Transmission Characteristics Under Monte Carlo Simulation.《2018 OCEANS - MTS/IEEE Kobe Techno-Oceans (OTO)》.2018, * |
海水水质对激光接收功率的影响分析;李晨等;《光电技术应用》;20170215(第01期);全文 * |
王菲 ; 赵彦.海水激光通信信道的特性仿真与分析.《通信电源技术》.2019, * |
Also Published As
Publication number | Publication date |
---|---|
CN110661568A (zh) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ding et al. | Modeling of non-line-of-sight ultraviolet scattering channels for communication | |
Gabriel et al. | Channel modeling for underwater optical communication | |
Choudhary et al. | Pathloss analysis of NLOS underwater wireless optical communication channel | |
Jagadeesh et al. | Characterization of channel impulse responses for nlos underwater wireless optical communications | |
CN110661568B (zh) | 水下到达激光信号3dB强度光斑半径的计算方法 | |
Elamassie et al. | Channel modeling and performance characterization of underwater visible light communications | |
CN111555822A (zh) | 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 | |
Vali et al. | Use of Gaussian beam divergence to compensate for misalignment of underwater wireless optical communication links | |
Ali | Characteristics of optical channel for underwater optical wireless communication system | |
Lin et al. | Dynamic optical wireless communication channel characterization through air-water interface | |
Dong et al. | Monte-Carlo-based impulse response modeling for underwater wireless optical communication | |
Qadar et al. | Monte Carlo–based channel estimation and performance evaluation for UWOC links under geometric losses | |
Geldard et al. | A study of spatial and temporal dispersion in turbulent underwater optical wireless channel | |
CN109302232A (zh) | 基于平行多跳架构的水下无线光通信系统误码率计算方法 | |
Mahmoud et al. | Statistical studies using goodness-of-fit techniques with dynamic underwater visible light communication channel modeling | |
Zhang et al. | A uniform spatial channel model for underwater wireless optical communication links | |
Qadar et al. | Monte Carlo based estimation and performance evaluation of temporal channel behavior of UWOC under multiple scattering | |
Zhang et al. | Angle of arrival analysis for underwater wireless optical links | |
CN116800364A (zh) | 基于射线追踪法的luna月球基地通信信道方法 | |
CN111723478A (zh) | 清澈海洋水质uwoc系统信道冲激响应拟合函数求解方法及系统 | |
Ma et al. | Modeling for underwater optical channel in relatively turbid waters based on the Monte Carlo method | |
Zhang et al. | On integrated stochastic channel model for underwater optical wireless communications | |
Wang et al. | Simulation of Photon Distribution in Underwater Laser Communication | |
Kumar et al. | Study of under-water Sonar System for change in propagation speed, depth of water, bottom loss and estimating optimal PDFs | |
Ao et al. | Study on Laser Propagation Characteristic of Downlink Channel in Non-uniform Seawater |
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 |