CN113297643A - 一种水下蓝绿激光通信中光斑面积的计算方法 - Google Patents
一种水下蓝绿激光通信中光斑面积的计算方法 Download PDFInfo
- Publication number
- CN113297643A CN113297643A CN202110401992.9A CN202110401992A CN113297643A CN 113297643 A CN113297643 A CN 113297643A CN 202110401992 A CN202110401992 A CN 202110401992A CN 113297643 A CN113297643 A CN 113297643A
- Authority
- CN
- China
- Prior art keywords
- photon
- photons
- formula
- scattering
- area
- 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
- 238000004891 communication Methods 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000005540 biological transmission Effects 0.000 claims abstract description 41
- 239000002245 particle Substances 0.000 claims abstract description 32
- 239000013535 sea water Substances 0.000 claims abstract description 24
- 238000010521 absorption reaction Methods 0.000 claims abstract description 19
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 239000000126 substance Substances 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- 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
-
- 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/80—Optical aspects relating to the use of optical transmission for specific applications, not provided for in groups H04B10/03 - H04B10/70, e.g. optical power feeding or optical transmission through water
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Geometry (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种水下蓝绿激光通信中光斑面积的计算方法,具体为:空间直角坐标系下,确定光源发射光子传输方向的余弦函数、光子的随机步长和光子第一次遇到粒子发生吸收或散射时的坐标;再得到光子遇到粒子发生吸收和散射的概率,重复以上步骤,直到光子经过n次散射到达接收面,光子在接收面上的坐标为光子发生第n+1次吸收或散射时的坐标,根据光子在接收面上的坐标分布,找出最外层光子,用线连接得到一个闭合曲线,闭合曲线围出的面积为接收面接收到的光斑的面积;利用紫外光在大气中的散射传输计算蓝绿激光在海水中传输,获得蓝绿激光在水下传输到达接收面上的光斑面积,从而为水下蓝绿激光通信中,接收端的透镜尺寸设计提供理论依据。
Description
技术领域
本发明属于水下无线激光通信技术领域,具体涉及一种水下蓝绿激光通信中光斑面积的计算方法。
背景技术
与陆上通信相比,电磁波在水中传播会发生严重的衰减,通信距离受到严重限制。水声通信虽然传输距离远,但是会产生严重的时延,且传输速率很慢。光通信弥补了电磁波在水中传输的高强度衰减以及声波在水中传输的高时延等缺点,是水下通信中一种新兴的通信技术。水下光通信有水下激光通信和水下可见光通信两种方式,使用的光波段均为450-550nm的蓝绿光。
光在海水中传输会受到海水中水分子、各种溶解物质和悬浮粒子的散射影响,改变光子的传输路径,使得光斑大小发生变化。以激光为例,蓝色激光器发出的蓝色激光在海水中传输,若在理想状态下,接收面接收到的光斑大小与初始激光器发出的光斑大小相一致。但在实际情况下并不是这样,接收面接收到的激光光斑大小与激光在海水中的传输方式(激光平行于海面传输、垂直于海面传输和斜程传输)和传输距离有关。相同传输距离下,不同的传输方式,光斑的大小不同;相同传输方式下,不同的传输距离,光斑的大小也不同。
发明内容
本发明的目的是提供一种水下蓝绿激光通信中光斑面积的计算方法,根据计算得到的激光光斑面积,为设计水下蓝绿激光通信中接收端的透镜尺寸提供理论依据。
本发明所采用的技术方案是,一种水下蓝绿激光通信中光斑面积的计算方法,具体按照以下步骤实施:
步骤1,在空间直角坐标系O-xyz下,光源的坐标为(x,y,z),确定光源发射光子传输方向的余弦函数(ux,uy,uz);
步骤2,确定光子的随机步长和光子第一次遇到粒子发生吸收或散射时的坐标;
步骤3:确定光子遇到粒子发生吸收和散射的概率,光子遇到粒子发生吸收的概率为P,则散射的概率为1-P;
步骤4,设激光器发出的光子有N个,光子第一次遇到粒子时,被吸收的光子数为NP,被散射的光子数为N(1-P);根据Henyey-Greenstein相位函数确定光子被粒子散射时,散射角θ的余弦函数;
步骤5,根据光子发生散射时散射角的余弦函数,确定光子新的传输方向的方向余弦;
步骤6,重复步骤2-5,直到光子经过n次散射到达接收面,得到光子第n次与粒子相遇时的坐标(xn,yn,zn),光子的数目N(1-P)n-1,以及光子发生第n次吸收或散射后,光子的数目N(1-P)n,光子传输方向的方向余弦uxn,uyn,uzn;
步骤7,光子在接收面上的坐标为光子发生第n+1次吸收或散射时的坐标(xn+1,yn+1,zn+1),光子在接收面上的数目为N(1-P)n,根据N(1-P)n个光子在接收面上的坐标分布,找出最外层光子,并用线连接,得到一个闭合曲线,闭合曲线围出的面积为接收面接收到的光斑的面积,最后用MATLAB中的polyarea函数计算出光斑面积。
本发明的特点还在于,
步骤1中,光子传输方向的余弦函数(ux,uy,uz),如式(1)、式(2)及式(3)所示;
ux=sin θ1 cosψ (1);
uy=sin θ1 sinψ (2);
uz=cos θ1 (3);
其中,ψ为光子的方位角,θ1为光子的偏转角,θ1在[0,φ/2]内取值。
步骤2中,具体为:由Bouguer-Lambert定律,结合蒙特卡罗方法,对随机步长进行抽样,如式(4)所示;
Δs=-lnξ(s)/bp (4);
其中,ξ(s)是在[0,1]上均匀分布的随机数,bp是海水的前向散射系数;ξ(s)的值由MATLAB中的rand函数随机生成;
确定光子传输的随机步长后,根据光子的传输方向,计算光子第一次遇到粒子发生吸收或散射时的坐标(x1,y1,z1),如式(5)、式(6)及式(7)所示;
x1=x+uxΔs (5);
y1=y+uyΔs (6);
z1=z+uzΔs (7)。
步骤3中,光子遇到粒子发生吸收的概率为P,如式(8)所示;
其中,a是海水的吸收系数,bp是海水的前向散射系数。
步骤4中,散射角θ的余弦函数如式(9)所示;
其中,Bp是后向散射比,Bp的计算公式如式(11)所示;
其中,bbp是海水的后向散射系数。
步骤5中,光子新的传输方向的方向余弦,式(12)、(13)及式(14)所示;
本发明的有益效果是,利用紫外光在大气中的散射传输计算蓝绿激光在海水中传输,获得蓝绿激光在水下传输到达接收面上的光斑面积,从而为水下蓝绿激光通信中,接收端的透镜尺寸设计提供理论依据。
附图说明
图1是蓝色激光在海水中传输受海水中水分子、各种溶解物质和悬浮粒子的吸收和散射影响的传输示意图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种水下蓝绿激光通信中光斑面积的计算方法,根据光子在接收面上的坐标分布,找出最外层光子,并用线连接,得到一个闭合曲线,闭合曲线围出的面积就是接收面接收到的光斑面积,最后,用MATLAB中的polyarea函数计算出光斑面积;
图1是蓝色激光平行于海面在海水中传输时,受海水中水分子、各种溶解物质和悬浮粒子的吸收和散射影响的传输示意图。图中,灰色粒子表示海水中水分子、各种溶解物质和悬浮粒子对激光的吸收,黑色粒子表示海水中水分子、各种溶解物质和悬浮粒子对激光的散射。
具体按照以下步骤实施:
在空间直角坐标系O-xyz下,光源的坐标为(x,y,z),确定光源发射光子传输方向的余弦函数(ux,uy,uz),如式(1)、式(2)及式(3)所示;
ux=sin θ1 cosψ (1);
uy=sin θ1 sinψ (2);
uz=cos θ1 (3);
其中,ψ为光子的方位角,θ1为光子的偏转角,θ1在[0,φ/2]内取值;
步骤2,确定光子的随机步长和光子第一次遇到粒子发生吸收或散射时的坐标;
具体为:由Bouguer-Lambert定律,结合蒙特卡罗方法,对随机步长进行抽样,如式(4)所示;
Δs=-lnξ(s)/bp (4);
其中,ξ(s)是在[0,1]上均匀分布的随机数,bp是海水的前向散射系数;ξ(s)的值由MATLAB中的rand函数随机生成;
确定光子传输的随机步长后,根据光源坐标(x,y,z)和光子的传输方向,计算光子第一次遇到粒子发生吸收或散射时的坐标(x1,y1,z1),如式(5)、式(6)及式(7)所示;
x1=x+uxΔs (5);
y1=y+uyΔs (6);
z1=z+uzΔs (7);
步骤3:确定光子遇到粒子发生吸收和散射的概率,光子遇到粒子发生吸收的概率为P,如式(8)所示,则散射的概率为1-P;
其中,a是海水的吸收系数;
步骤4,设激光器发出的光子有N个,光子第一次遇到粒子时,被吸收的光子数为NP,被散射的光子数为N(1-P);根据Henyey-Greenstein相位函数确定光子被粒子散射时,散射角θ的余弦函数,如式(9)所示;
其中,Bp是后向散射比,Bp的计算公式如式(11)所示;
其中,bbp是海水的后向散射系数,bp是海水的前向散射系数。
步骤5:根据光子发生散射时散射角的余弦函数,确定光子新的传输方向的方向余弦,如式(12)、(13)及式(14)所示;
步骤6,重复步骤2-5,直到光子经过n次散射到达接收面,则可得到光子第n次与粒子相遇时的坐标(xn,yn,zn),光子的数目N(1-P)n-1,以及光子发生第n次散射后,光子的数目N(1-P)n,光子传输方向的方向余弦uxn,uyn,uzn;
步骤7,光子在接收面上的坐标为光子发生第n+1次吸收或散射时的坐标(xn+1,yn+1,zn+1),光子在接收面上的数目为N(1-P)n。根据N(1-P)n个光子在接收面上的坐标分布,找出最外层光子,并用线连接,得到一个闭合曲线,闭合曲线围出的面积就是接收面接收到的光斑的面积,最后用MATLAB中的polyarea函数计算出光斑面积。
Claims (6)
1.一种水下蓝绿激光通信中光斑面积的计算方法,其特征在于,具体按照以下步骤实施:
步骤1,在空间直角坐标系O-xyz下,光源的坐标为(x,y,z),确定光源发射光子传输方向的余弦函数(ux,uy,uz);
步骤2,确定光子的随机步长和光子第一次遇到粒子发生吸收或散射时的坐标;
步骤3:确定光子遇到粒子发生吸收和散射的概率,光子遇到粒子发生吸收的概率为P,则散射的概率为1-P;
步骤4,设激光器发出的光子有N个,光子第一次遇到粒子时,被吸收的光子数为NP,被散射的光子数为N(1-P);根据Henyey-Greenstein相位函数确定光子被粒子散射时,散射角θ的余弦函数;
步骤5:根据光子发生散射时散射角的余弦函数,确定光子新的传输方向的方向余弦;
步骤6,重复步骤2-5,直到光子经过n次散射到达接收面,得到光子第n次与粒子相遇时的坐标(xn,yn,zn),光子的数目N(1-P)n-1,以及光子发生第n次吸收或散射后,光子的数目N(1-P)n,光子传输方向的方向余弦uxn,uyn,uzn;
步骤7,光子在接收面上的坐标为光子发生第n+1次吸收或散射时的坐标(xn+1,yn+1,zn+1),光子在接收面上的数目为N(1-P)n,根据N(1-P)n个光子在接收面上的坐标分布,找出最外层光子,并用线连接,得到一个闭合曲线,闭合曲线围出的面积为接收面接收到的光斑的面积,最后用MATLAB中的polyarea函数计算出光斑面积。
2.根据权利要求1所述的一种水下蓝绿激光通信中光斑面积的计算方法,其特征在于,所述步骤1中,光子传输方向的余弦函数(ux,uy,uz),如式(1)、式(2)及式(3)所示;
ux=sinθ1cosψ (1);
uy=sinθ1sinψ (2);
uz=cosθ1 (3);
其中,ψ为光子的方位角,θ1为光子的偏转角,θ1在[0,φ/2]内取值。
3.根据权利要求2所述的一种水下蓝绿激光通信中光斑面积的计算方法,其特征在于,所述步骤2中,具体为:由Bouguer-Lambert定律,结合蒙特卡罗方法,对随机步长进行抽样,如式(4)所示;
Δs=-lnξ(s)/bp (4);
其中,ξ(s)是在[0,1]上均匀分布的随机数,bp是海水的前向散射系数;ξ(s)的值由MATLAB中的rand函数随机生成;
确定光子传输的随机步长后,根据光子的传输方向,计算光子第一次遇到粒子发生吸收或散射时的坐标(x1,y1,z1),如式(5)、式(6)及式(7)所示;
x1=x+uxΔs (5);
y1=y+uyΔs (6);
z1=z+uzΔs (7)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110401992.9A CN113297643B (zh) | 2021-04-14 | 2021-04-14 | 一种水下蓝绿激光通信中光斑面积的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110401992.9A CN113297643B (zh) | 2021-04-14 | 2021-04-14 | 一种水下蓝绿激光通信中光斑面积的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113297643A true CN113297643A (zh) | 2021-08-24 |
CN113297643B CN113297643B (zh) | 2023-02-03 |
Family
ID=77319743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110401992.9A Active CN113297643B (zh) | 2021-04-14 | 2021-04-14 | 一种水下蓝绿激光通信中光斑面积的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113297643B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2533538C1 (ru) * | 2013-08-19 | 2014-11-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тверской государственный университет" | Способ раздельного определения вероятностей поглощения и рассеяния фотонов на единицу пути в твердых оптических материалах |
US20150345999A1 (en) * | 2014-05-28 | 2015-12-03 | Pruftechnik Dieter Busch Ag | Method for determining a closed trajectory by means of a laser and a laser light sensor and apparatus for determining a closed trajectory curve |
CN109995427A (zh) * | 2019-03-25 | 2019-07-09 | 西安电子科技大学 | 一种水下上行激光通信的蒙特卡洛仿真方法 |
CN110661568A (zh) * | 2019-10-28 | 2020-01-07 | 桂林电子科技大学 | 水下到达激光信号3dB强度光斑半径的计算方法 |
CN111555822A (zh) * | 2020-04-28 | 2020-08-18 | 西安邮电大学 | 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 |
CN112235044A (zh) * | 2020-09-21 | 2021-01-15 | 西安理工大学 | 一种水下激光通信系统的信道建模方法 |
CN112539698A (zh) * | 2020-11-09 | 2021-03-23 | 北京工业大学 | 一种激光光束作用材料内部在线跟踪与实时反馈的方法 |
-
2021
- 2021-04-14 CN CN202110401992.9A patent/CN113297643B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2533538C1 (ru) * | 2013-08-19 | 2014-11-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тверской государственный университет" | Способ раздельного определения вероятностей поглощения и рассеяния фотонов на единицу пути в твердых оптических материалах |
US20150345999A1 (en) * | 2014-05-28 | 2015-12-03 | Pruftechnik Dieter Busch Ag | Method for determining a closed trajectory by means of a laser and a laser light sensor and apparatus for determining a closed trajectory curve |
CN109995427A (zh) * | 2019-03-25 | 2019-07-09 | 西安电子科技大学 | 一种水下上行激光通信的蒙特卡洛仿真方法 |
CN110661568A (zh) * | 2019-10-28 | 2020-01-07 | 桂林电子科技大学 | 水下到达激光信号3dB强度光斑半径的计算方法 |
CN111555822A (zh) * | 2020-04-28 | 2020-08-18 | 西安邮电大学 | 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 |
CN112235044A (zh) * | 2020-09-21 | 2021-01-15 | 西安理工大学 | 一种水下激光通信系统的信道建模方法 |
CN112539698A (zh) * | 2020-11-09 | 2021-03-23 | 北京工业大学 | 一种激光光束作用材料内部在线跟踪与实时反馈的方法 |
Non-Patent Citations (3)
Title |
---|
ELDAR MUSA: "Line laser-based break sensor that detects light spots on yarns", 《OPTICS AND LASERS IN ENGINEERING》 * |
刘力: "跨介质蓝绿激光上行通信系统研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
詹恩奇等: "星载激光脉冲水下光斑分布的研究", 《武汉理工大学学报(交通科学与工程版)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113297643B (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guo et al. | Channel-parameter estimation for satellite-to-submarine continuous-variable quantum key distribution | |
US20190271059A1 (en) | Laser Shock Peening Apparatus for Surface of Workpiece, and Laser Shock Peening Method | |
CN111555822B (zh) | 一种基于相位屏的水下无线光传输蒙特卡洛仿真方法 | |
CN109327271B (zh) | 基于部分相干光的水下光通信系统平均信道容量确定方法 | |
CN108023652B (zh) | 一种应用于海水信道的激光传输特性的模拟方法 | |
CN112235044B (zh) | 一种水下激光通信系统的信道建模方法 | |
Choudhary et al. | Pathloss analysis of NLOS underwater wireless optical communication channel | |
JP2012531824A (ja) | 地表上及び地表下の波面管理 | |
Khalighi et al. | Underwater visible light communications, channel modeling and system design | |
CN113297643B (zh) | 一种水下蓝绿激光通信中光斑面积的计算方法 | |
Zhang et al. | Angle of arrival analysis for underwater wireless optical links | |
ZA202310738B (en) | Target object cross-media imaging simulation method based on monte carlo | |
Fang et al. | Investigation of non-line-of-sight underwater optical wireless communications with wavy surface | |
CN113078957B (zh) | 一种基于空间坐标变换技术的uwoc系统信道特征mcns方法 | |
Guerra et al. | Effect of absorbing and non‐diffracting particles in UWOC links | |
Zhou et al. | Single scattering impulse response modeling of underwater wireless optical channels | |
CN111025272B (zh) | 具备隧道效应抑制能力的平面声学基阵超宽覆盖波束发射方法 | |
Jones et al. | Tailoring underwater laser acoustic pulses | |
Yokoyama et al. | Examination for the application of Laguerre-Gaussian beams for underwater optical wireless communication | |
Wang et al. | Performance research for quantum key distribution based on real air-water channel | |
CN114791992A (zh) | 一种基于射线理论的深海目标前向散射声场计算方法 | |
CN118094964B (zh) | 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 | |
Ma et al. | Research on channel modeling method of underwater optical communication | |
CN113740873B (zh) | 一种基于高斯卷积的海洋激光雷达快速仿真方法 | |
Schildkraut et al. | Investigating the Role of Spatiotemporal Optical Beam Profiles in Mixed Layer Oceanic Communication Channels |
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 |