CN105680936B - 一种非直视紫外光通信单次散射过程路径损耗计算方法 - Google Patents

Abstract

本发明公开了一种非直视紫外光通信单次散射过程路径损耗计算方法，具体按照以下步骤实施：步骤1：设定发射端和接收端的参数，定义区域V′，公共散射体求出V′的范围；步骤2：将区域V′分为若干个微元V"，在球坐标系中，得到微元的坐标；步骤3：遍历区域V′分割成的若干个微元V"，判断微元V"是否在公共散射体V内，如果在公共散射体V内计算被微元V"散射后被Rx接收到的能量，将被微元V"散射后被Rx接收到的能量相加，计算得到Rx接收到的总能量，求出路径损耗。本发明一种非直视紫外光通信单次散射过程路径损耗计算方法，仿真时间短，能够很容易的计算出通信系统的路径损耗，与MC方法的仿真结果拟合度高。

Description

一种非直视紫外光通信单次散射过程路径损耗计算方法

技术领域

本发明属于紫外光通信技术领域，具体涉及一种非直视紫外光通信单次散射过程路径损耗计算方法。

背景技术

紫外光通信就是把紫外光作为信息传输的载体，将需要传输的信息加载到紫外光上，以实现信息的发送和接收。在无线激光通信中，采用的是直视的通信方式，即发射端与接收端必须对准。非直视通信是指发射端发射出的紫外光子在大气中传输时，由于紫外光波长很短，被大气粒子散射而能够绕过障碍物到达接收端的过程。只要发射端发散角能够与接收端接收视场角在空中形成公共散射体，那么非直视(non-line-of-sight，NLOS)通信就变得可能。当直视通信无法实现时，NLOS通信则能保证发射端和接收端之间通信的畅通。

关于NLOS紫外光大气传输的理论模型现在采用的主要方法有单次散射近似法和蒙特卡罗(Monte Carlo，MC)方法。经典的单次散射模型通过三重积分求通信系统的路径损耗(path loss，PL)，但是，在仿真过程中复杂的三重积分很难实现。如果用简化公式来代替三重积分，要求公共散射体是闭合的，发散角和接收视场角比较小。MC方法是一种以概率统计理论为基础的计算方法，可以通过MC方法来对单次散射传输过程进行仿真，但是存在计算量大，耗时久的问题。

发明内容

本发明的目的是提供一种非直视紫外光通信单次散射过程路径损耗计算方法，解决了现有的经典单次散射模型仿真困难和MC算法仿真时间过长的问题。

本发明所采用的技术方案是，一种非直视紫外光通信单次散射过程路径损耗计算方法，具体按照以下步骤实施：

步骤1：设定发射端和接收端的参数，定义区域V′，公共散射体求出V′的范围；

步骤2：将区域V′分为若干个微元V"，在球坐标系中，得到微元的坐标；

步骤3：遍历区域V′分割成的若干个微元V"，判断微元V"是否在公共散射体V内，如果在公共散射体V内计算被微元V"散射后被Rx接收到的能量，将被微元V"散射后被Rx接收到的能量相加，计算得到Rx接收到的总能量，求出路径损耗。

本发明的特点还在于，

步骤1具体为：

在紫外光NLOS通信中，C_T和C_R分别表示发射端Tx的发散角圆锥和接收端Rx的视场角FOV圆锥，点T为C_T的顶点，点R为C_R的顶点，点H为C_T的任意横切面的圆心，点L为C_R的任意横切面的圆心，即TH为C_T的中心轴线，RL为C_R的中心轴线，发射端发散半角为发射仰角为TH和其在XOY平面投影的夹角θ_t，坐标原点O与点T重合，接收端视场半角为接收仰角为RL和其在XOY平面投影的夹角θ_r，TH在XOY平面的投影与X轴的夹角α_t是C_T的偏轴角，RL在XOY平面的投影与X轴的夹角α_r是C_R的偏轴角，d是点T到点R的直线距离，(θ_t,α_t)和(θ_r,α_r)确定了C_T和C_R的方向；

发散角圆锥和FOV圆锥的公共区域为公共散射体V，即V＝C_T∩C_R，点S为公共散射体V内的散射点，发射端发射一个光子，光子在点S被散射，在以坐标原点为参考点，由方位角、顶角和距离构成球面坐标系下，OS与Z轴的夹角θ为光子入射的顶角，OS在XOY平面的投影与X轴的夹角α为方位角，从发射端T到点S的距离为r，光子的入射方向和指向接收端的散射方向构成的夹角为散射角β_s，ζ为散射点S和点R的连线与RL构成的夹角，r'是从点S到接收端R的距离；

光子在传输的过程中，遇到V内的体积微元δV发生散射，到达Rx的能量dE_r为：

其中是Rx的立体角，E_t是发射光束能量，A_r是接收探测面面积，P(cosβ_s)是散射相函数，k_e是消光系数，吸收系数k_a和散射系数k_s之和构成了通信过程中大气的消光系数k_e，即k_e＝k_s+k_a；

散射相函数P(cosβ_s)是瑞利散射相函数P^R(cosβ_s)和米氏散射相函数P^M(cosβ_s)的加权和：

其中，k_s ^R是瑞利散射的散射系数，k_s ^M是米氏散射的散射系数，k_s＝k_s ^R+k_s ^M；瑞利散射相函数P^R(cosβ_s)为：

米氏散射相函数P^M(cosβ_s)为：

其中，γ、g和f是模型参数；

在球坐标系中，体积微元为δV＝r²sinθδθδαδr，

则，被公共散射体V散射后被Rx接收到的总能量E_r为：

定义区域V′，范围为[θ_min,θ_max]、[α_min,α_max]、[r_min,r_max]，使

θ_min和θ_max为：

过点H做线GH平行于XOY平面，G点为线GH与C_T的交点，G点在XOY平面的投影为G'，在线TH上选取任意一点E，E点在XOY平面的投影为E'，应用三角定理，则，α_min和α_max为：

对于共面的情况，公共散射体的体积是最大的，TP是发散角圆锥的中心轴线，RQ是FOV圆锥的中心轴线，P点在XOY平面的投影为P'，Q点在XOY平面的投影为Q'，点K、M、U、W为C_T和C_R的四个交点，∠PTP′＝θ_t,∠QRQ′＝θ_r,∠UTP和∠PTM等于∠KRQ和∠QRM等于所以，对于△UTR，由正弦函数可得，UT为：

同样的，

另

r_min＝min[UT,MT],r_max＝max[WT,KT](10)

对于无界的公共散射体V，r_max→∞，

至此，区域V′的范围可得。

步骤2具体为：

设置3个整数N_r、N_α、N_θ，区域[r_min,r_max]、[α_min,α_max]和[θ_min,θ_max]分别被分为N_r、N_α、N_θ等份，因此区域V′被分为了N_rN_αN_θ个微元V"，另N_r、N_α、N_θ均等于N，那么区域V′就被分为了N³个微元V"，假设3个整数i、j、k，且1≤i,j,k≤N，i为区域[r_min,r_max]被分的第i份，j为区域[α_min,α_max]被分的第j份，k为区域[θ_min,θ_max]被分的第k份，[i,j,k]代表一个微元V"，A、B、C、D为区域V′内的一微元V"的四个顶点，点S是V"的中心，设置3个变量r_u，α_u，θ_u，BC的长度r_u＝(r_max-r_min)/N；∠ATB＝(θ_max-θ_min)/N＝θ_u；点B和点D在XOY平面的投影分别为B'和D'，∠B′TD′＝(α_max-α_min)/N＝α_u，在球坐标系中，点S、A、B、C、D的坐标可由如下公式得到：

[r_S,α_S,q_S]＝[r_min+r_u(i-1/2),α_min+α_u(j-1/2),q_min+q_u(k-1/2)]

[r_A,α_A,q_A]＝[r_min+r_u(i-1),α_min+α_uj,q_min+q_u(k-1)]

[r_B,α_B,q_B]＝[r_min+r_u(i-1),α_min+α_uj,q_min+q_uk] (11)

[r_C,α_C,q_C]＝[r_min+r_ui,α_min+α_uj,q_min+q_uk]

[r_D,α_D,q_D]＝[r_min+r_u(i-1),α_min+α_u(j-1),q_min+q_uk]

在笛卡尔坐标系中，点S的坐标可由如下公式得到：

[x_S,yS,z_S]＝[r_S sin(q_S)cos(α_S),r_S sin(q_S)sin(α_S),r_S cos(q_S)] (12)。

步骤3具体为：

已知点T的坐标为(0,0,0)，点R的坐标为(d,0,0)，定义[T→E]和[R→F]是点T到点E和点R到点F的单元方向矢量，定义[T→S]和[R→S]是点T到点S和点R到点S的矢量，公式如下所示：

[T→E]＝[cosq_tcosα_t,cosq_tsinα_t,sinq_t]

[R→F]＝[cosq_rcosα_r,cosq_rsinα_r,sinq_r] (13)

[T→S]＝[x_S,y_S,z_S]

[R→S]＝[x_S-d,y_S,z_S]

则

如果为真，那么中心点S在区域V内；

因为微元V"非常小，可以将微元近似认为就是点S，则散射角β_s近似为：

其中[S→R]＝[d-x_S,-yS,-z_S]，

由公式(2)可得被微元V"散射后被Rx接收到的能量E_V"r为：

其中，ζ＝∠SRF,r'＝|R→S|；

遍历区域V′分割成的N³个微元V"，如果为真，由公式(16)可以计算得出E_V"r，因此，Rx接收到的总能量E_r为：

路径损耗PL可为：

PL＝10lg(E_t/E_r) (18)。

本发明的有益效果是：本发明一种非直视紫外光通信单次散射过程路径损耗计算方法，实现了NLOS紫外光通信单次散射过程路径损耗的计算，不需要根据发射端和接收端的几何参数，把三重积分分成多种情况，对模型中的几何参数无限制。路径损耗计算方法中遍历微元法中积分的上限和下限为定值，被积函数被简化，将分区域三重积分运算转化为乘法运算的累加和。本发明一种非直视紫外光通信单次散射过程路径损耗计算方法，仿真时间短，能够很容易的计算出通信系统的路径损耗，与MC方法的仿真结果拟合度高。

附图说明

图1是本发明路径损耗计算方法中NLOS非共面传输模型图；

图2是本发明路径损耗计算方法中NLOS共面传输模型图；

图3是紫外光非直视非共面传输模型的仿真图；

图4是本发明实施例中紫外光单次散射遍历微元法路径损耗仿真结果对比图。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。

本发明一种非直视紫外光通信单次散射过程路径损耗计算方法，具体按照以下步骤实施：

步骤1：设定发射端和接收端的参数，定义区域V′，公共散射体求出V′的范围。

在紫外光NLOS(non-line-of-sight，非直视)通信中，NLOS非共面单次散射模型图如图1所示。C_T和C_R分别表示发射端(Tx)的发散角圆锥和接收端(Rx)的视场角(field ofview，FOV)圆锥。点T为C_T的顶点，点R为C_R的顶点，点H为C_T的任意横切面的圆心，点L为C_R的任意横切面的圆心，即TH为C_T的中心轴线，RL为C_R的中心轴线。发射端发散半角为发射仰角为TH和其在XOY平面投影的夹角θ_t，坐标原点O与点T重合。接收端视场半角为接收仰角为RL和其在XOY平面投影的夹角θ_r。TH在XOY平面的投影与X轴的夹角α_t是C_T的偏轴角，RL在XOY平面的投影与X轴的夹角α_r是C_R的偏轴角。d是点T到点R的直线距离。(θ_t,α_t)和(θ_r,α_r)确定了C_T和C_R的方向。

定义发散角圆锥和FOV圆锥的公共区域为公共散射体V，也就是V＝C_T∩C_R。定义点S为公共散射体V内的散射点。发射端发射一个光子，光子在点S被散射。在以坐标原点为参考点，由方位角、顶角和距离构成球面坐标系下，定义关于这一段射线的各个参数：OS与Z轴的夹角θ为光子入射的顶角，OS在XOY平面的投影与X轴的夹角α为方位角，从发射端T到点S的距离定义为r。定义光子的入射方向和指向接收端的散射方向构成的夹角为散射角β_s。ζ为散射点S和点R的连线和RL构成的夹角。r'是从点S到接收端R的距离。

光子在传输的过程中，遇到V内的体积微元δV发生散射，到达Rx的能量可由如下公式得到

其中是Rx的立体角，E_t是发射光束能量，A_r是接收探测面面积，P(cosβ_s)是散射相函数，k_e是消光系数，吸收系数k_a和散射系数k_s之和构成了通信过程中大气的消光系数k_e，也就是k_e＝k_s+k_a。

散射相函数P(cosβ_s)是瑞利散射(分子)相函数P^R(cosβ_s)和米氏散射(气溶胶)相函数P^M(cosβ_s)的加权和：

其中，k_s ^R是瑞利散射的散射系数，k_s ^M是米氏散射的散射系数，k_s＝k_s ^R+k_s ^M。瑞利散射和米氏散射的相函数分别为：

其中，γ、g和f是模型参数。波长为266nm时，(γ,g,f)为(0.017,0.72,0.5)，(k_a,k_s ^R,k_s ^M)为(0.74,0.24,0.25)×10^-3m^-1。

在球坐标系中，体积微元为δV＝r²sinθδθδαδr。因此，被公共散射体V散射后被Rx接收到的总能量E_r为：

定义区域V′，范围为[θ_min,θ_max]、[α_min,α_max]、[r_min,r_max]，让从图1可知θ_min和θ_max为

接下来，分析α_min和α_max。如图1所示，过点H做线GH平行于XOY平面，G点为线GH与C_T的交点，G点在XOY平面的投影为G'，在线TH上选取任意一点E，E点在XOY平面的投影为E'，应用三角定理，所以，我们可以得到α_min和α_max，公式如下所示：

最后，我们给出r_min和r_max的设置。对于共面的情况，如图2所示，此时公共散射体的体积是最大的，TP是发散角圆锥的中心轴线，RQ是FOV圆锥的中心轴线，P点在XOY平面的投影为P'，Q点在XOY平面的投影为Q'。点K、M、U、W为C_T和C_R的四个交点。∠PTP′＝θ_t,∠QRQ′＝θ_r,∠UTP和∠PTM等于∠KRQ和∠QRM等于所以， 对于△UTR，由正弦函数可得，UT为

同样的，

让

r_min＝min[UT,MT],r_max＝max[WT,KT] (10)

对于无界的公共散射体V，r_max→∞。

至此，区域V′的范围可得。

步骤2：将区域V′分为若干个微元V"，在球坐标系中，得到微元的坐标

设置3个整数N_r、N_α、N_θ，区域[r_min,r_max]、[α_min,α_max]和[θ_min,θ_max]分别被分为N_r、N_α、N_θ等份，因此区域V′被分为了N_rN_αN_θ个微元V"，如果让N_r、N_α、N_θ都等于N，那么区域V′就被分为了N³个微元V"，假设3个整数i、j、k，且1≤i,j,k≤N，i代表的是区域[r_min,r_max]被分的第i份，j代表的是区域[α_min,α_max]被分的第j份，k代表的是区域[θ_min,θ_max]被分的第k份，[i,j,k]代表一个微元V"。如图1所示，A、B、C、D为区域V′内的一微元V"的四个顶点，点S是V"的中心。设置3个变量r_u，α_u，θ_u，BC的长度r_u＝(r_max-r_min)/N；∠ATB＝(θ_max-θ_min)/N＝θ_u；点B和点D在XOY平面的投影分别为B'和D'，∠B′TD′＝(α_max-α_min)/N＝α_u。在球坐标系中，点S、A、B、C、D的坐标可由如下公式得到：

[r_S,α_S,q_S]＝[r_min+r_u(i-1/2),α_min+α_u(j-1/2),q_min+q_u(k-1/2)]

[r_A,α_A,q_A]＝[r_min+r_u(i-1),α_min+α_uj,q_min+q_u(k-1)]

[r_B,α_B,q_B]＝[r_min+r_u(i-1),α_min+α_uj,q_min+q_uk] (11)

[r_C,α_C,q_C]＝[r_min+r_ui,α_min+α_uj,q_min+q_uk]

[r_D,α_D,q_D]＝[r_min+r_u(i-1),α_min+α_u(j-1),q_min+q_uk]

在笛卡尔坐标系中，点S的坐标可由如下公式得到：

[x_S,y_S,z_S]＝[r_Ssin(q_S)cos(α_S),r_Ssin(q_S)sin(α_S),r_Scos(q_S)] (12)

步骤3：遍历区域V′分割成的若干个微元V"，判断微元V"是否在公共散射体V内，如果在公共散射体V内计算被微元V"散射后被Rx接收到的能量，将被微元V"散射后被Rx接收到的能量相加，计算得到Rx接收到的总能量，求出路径损耗。

已知点T的坐标为(0,0,0)，点R的坐标为(d,0,0)，定义[T→E]和[R→F]是点T到点E和点R到点F的单元方向矢量，定义[T→S]和[R→S]是点T到点S和点R到点S的矢量，公式如下所示：

[T→E]＝[cosq_tcosα_t,cosq_tsinα_t,sinq_t]

[R→F]＝[cosq_rcosα_r,cosq_rsinα_r,sinq_r] (13)

[T→S]＝[x_S,y_S,z_S]

[R→S]＝[x_S-d,y_S,z_S]

则

如果为真，那么中心点S在区域V内。

图3为紫外光非直视非共面传输模型的仿真图，图中*表示的是微元的中心点，判断微元是否在区域V内，如果不在区域V内，则不在图中显示；如果在区域V内，则在图中显示。

因为微元V"非常小，可以将微元近似认为就是点S，则散射角β_s近似为：

其中[S→R]＝[d-x_S,-y_S,-z_S]，

由公式(2)可得被微元V"散射后被Rx接收到的能量E_V"r为：

其中，ζ＝∠SRF,r'＝|R→S|。

遍历区域V′分割成的N³个微元V"，如果为真，由公式(16)可以计算得出E_V"r，所以将该方法定义为遍历微元法。因此，Rx接收到的总能量E_r可由如下公式计算得到

路径损耗可由如下公式计算得到

PL＝10lg(E_t/E_r) (18)

实施例

为了验证遍历微元法的性能，我们将该算法与MC单次散射方法进行了相关的路径损耗的仿真和比较。

设置Rx的探测面半径为1.5×10^-2m，d＝100m，θ_t＝30°，θ_r＝60°，α_t＝0°，α_r＝{0°，10°，20°，30°，40°}。仿真结果对比图如图4所示，左边的纵坐标是路径损耗(PL)，右边的纵坐标是公共散射体体积|V|。图中的四条曲线分别表示的是利用蒙特卡罗MC多次散射方法得到的路径损耗的变化曲线、利用MC单次散射方法得到的路径损耗变化曲线、利用遍历微元法得到的路径损耗变化曲线以及利用遍历微元法得到的公共散射体体积变化曲线。注意到随着α_r的增加，PL增加，|V|逐渐减小，由遍历微元法得到的PL曲线与MC单次散射方法得到的PL曲线匹配度很高，并且略大于MC多次散射方法得到的PL曲线。

在MC仿真中，设置发射端发射的光子数为10⁶，通过MC方法计算PL所需的平均时间为107.514秒。基于之前对遍历微元法的分析，将N设置为60。用同一台电脑，设置相同的参数，采用遍历微元法计算PL仅需11.078秒。结果显示遍历微元法能够节省89.70％的时间。

通过验证可以看出遍历微元法仿真结果与MC方法仿真结果匹配度非常高，并且仿真时间更短。

Claims (3)

1.一种非直视紫外光通信单次散射过程路径损耗计算方法，其特征在于，具体按照以下步骤实施：

步骤1：设定发射端和接收端的参数，定义区域V′，公共散射体求出V′的范围；

所述步骤1具体为：

在紫外光NLOS通信中，C_T和C_R分别表示发射端Tx的发散角圆锥和接收端Rx的视场角FOV圆锥，点T为C_T的顶点，点R为C_R的顶点，点H为C_T的任意横切面的圆心，点L为C_R的任意横切面的圆心，即TH为C_T的中心轴线，RL为C_R的中心轴线，发射端发散半角为发射仰角为TH和其在XOY平面投影的夹角θ_t，坐标原点O与点T重合，接收端视场半角为接收仰角为RL和其在XOY平面投影的夹角θ_r，TH在XOY平面的投影与X轴的夹角α_t是C_T的偏轴角，RL在XOY平面的投影与X轴的夹角α_r是C_R的偏轴角，d是点T到点R的直线距离，(θ_t,α_t)和(θ_r,α_r)确定了C_T和C_R的方向；

发散角圆锥和FOV圆锥的公共区域为公共散射体V，即V＝C_T∩C_R，点S为公共散射体V内的散射点，发射端发射一个光子，光子在点S被散射，在以坐标原点为参考点，由方位角、顶角和距离构成球面坐标系下，OS与Z轴的夹角θ为光子入射的顶角，OS在XOY平面的投影与X轴的夹角α为方位角，从发射端T到点S的距离为r，光子的入射方向和指向接收端的散射方向构成的夹角为散射角β_s，ζ为散射点S和点R的连线与RL构成的夹角，r'是从点S到接收端R的距离；

光子在传输的过程中，遇到V内的体积微元δV发生散射，到达Rx的能量δE_r为：

其中是Rx的立体角，E_t是发射光束能量，A_r是接收探测面面积，P(cosβ_s)是散射相函数，k_e是消光系数，吸收系数k_a和散射系数k_s之和构成了通信过程中大气的消光系数k_e，即k_e＝k_s+k_a；

散射相函数P(cosβ_s)是瑞利散射相函数P^R(cosβ_s)和米氏散射相函数P^M(cosβ_s)的加权和：

其中，k_s ^R是瑞利散射的散射系数，k_s ^M是米氏散射的散射系数，k_s＝k_s ^R+k_s ^M；瑞利散射相函数P^R(cosβ_s)为：

米氏散射相函数P^M(cosβ_s)为：

其中，γ、g和f是模型参数；

在球坐标系中，体积微元为δV＝r²sinθδθδαδr，

则，被公共散射体V散射后被Rx接收到的总能量E_r为：

定义区域V′，范围为[θ_min,θ_max]、[α_min,α_max]、[r_min,r_max]，使

θ_min和θ_max为：

过点H做线GH平行于XOY平面，G点为线GH与C_T的交点，G点在XOY平面的投影为G'，在线TH上选取任意一点E，E点在XOY平面的投影为E'，应用三角定理，则，α_min和α_max为：

对于共面的情况，公共散射体的体积是最大的，TP是发散角圆锥的中心轴线，RQ是FOV圆锥的中心轴线，P点在XOY平面的投影为P'，Q点在XOY平面的投影为Q'，点K、M、U、W为C_T和C_R的四个交点，∠PTP′＝θ_t,∠QRQ′＝θ_r,∠UTP和∠PTM等于∠KRQ和∠QRM等于所以，对于△UTR，由正弦函数可得，UT为：

同样的，

令

r_min＝min[UT,MT],r_max＝max[WT,KT] (10)

对于无界的公共散射体V，r_max→∞，

至此，区域V′的范围可得；

步骤2：将区域V′分为若干个微元V"，在球坐标系中，得到微元的坐标；

步骤3：遍历区域V′分割成的若干个微元V"，判断微元V"是否在公共散射体V内，如果在公共散射体V内计算被微元V"散射后被Rx接收到的能量，将被微元V"散射后被Rx接收到的能量相加，计算得到Rx接收到的总能量，求出路径损耗。

2.根据权利要求1所述的一种非直视紫外光通信单次散射过程路径损耗计算方法，其特征在于，所述步骤2具体为：

设置3个整数N_r、N_α、N_θ，区域[r_min,r_max]、[α_min,α_max]和[θ_min,θ_max]分别被分为N_r、N_α、N_θ等份，因此区域V′被分为了N_rN_αN_θ个微元V"，令N_r、N_α、N_θ均等于N，那么区域V′就被分为了N³个微元V"，假设3个整数i、j、k，且1≤i,j,k≤N，i为区域[r_min,r_max]被分的第i份，j为区域[α_min,α_max]被分的第j份，k为区域[θ_min,θ_max]被分的第k份，[i,j,k]代表一个微元V"，A、B、C、D为区域V′内的一微元V"的四个顶点，点S是V"的中心，设置3个变量r_u，α_u，θ_u，BC的长度r_u＝(r_max-r_min)/N；∠ATB＝(θ_max-θ_min)/N＝θ_u；点B和点D在XOY平面的投影分别为B'和D'，∠B′TD′＝(α_max-α_min)/N＝α_u，在球坐标系中，点S、A、B、C、D的坐标可由如下公式得到：

在笛卡尔坐标系中，点S的坐标可由如下公式得到：

[x_S,y_S,z_S]＝[r_Ssin(θ_S)cos(α_S),r_Ssin(θ_S)sin(α_S),r_Scos(θ_S)] (12)。

3.根据权利要求2所述的一种非直视紫外光通信单次散射过程路径损耗计算方法，其特征在于，所述步骤3具体为：

已知点T的坐标为(0,0,0)，点R的坐标为(d,0,0)，定义[T→E]和[R→F]是点T到点E和点R到点F的单元方向矢量，定义[T→S]和[R→S]是点T到点S和点R到点S的矢量，公式如下所示：

则

如果为真，那么中心点S在区域V内；

因为微元V"非常小，可以将微元近似认为就是点S，则散射角β_s近似为：

其中[S→R]＝[d-x_S,-y_S,-z_S]，

由公式(2)可得被微元V"散射后被Rx接收到的能量E_V"r为：

其中，ζ＝∠SRF,r'＝|R→S|；

遍历区域V′分割成的N³个微元V"，如果为真，由公式(16)可以计算得出E_V"r，因此，Rx接收到的总能量E_r为：

路径损耗PL可为：

PL＝10lg(E_t/E_r) (18)。