当前位置: 首页 > 专利查询>清华大学专利>正文

一种无线光非视距通信链路中存在反射面的信道特性分析方法技术

技术编号:29460736 阅读:29 留言:0更新日期:2021-07-27 17:30
本发明专利技术公开了一种信道特性分析方法,该方法针对在无线光非视距通信链路中存在反射面时的信道特性,以蒙特卡洛方法为基础,通过生成大量光子进行模拟光线传输的特性分析,从而获得无线光非视距通信链路中存在反射面时的信道特性,所述特性包括通信链路路径损耗和任意时间段内的通信链路信道响应结果。本发明专利技术对估计无线光非视距通信链路的性能具有十分重要的意义。

A channel characteristic analysis method with reflector in wireless optical NLOS communication link

【技术实现步骤摘要】
一种无线光非视距通信链路中存在反射面的信道特性分析方法
技术介绍
传统无线光非视距通信链路特性的分析模型中仅考虑了光子与大气中气溶胶粒子和分子发射生的散射与吸收作用。然而,在实际应用中,由于空间环境存在各种建筑墙面、树冠面、以及地面等,光子无法避免与这些表面发生相互作用。找寻一种可以考虑这些相互作用,进行更准确科学的信道特性分析方法具有重要的现实意义。专利技术目的本专利技术的目的即在于克服现有技术中所存在的困难,针对在无线光非视距通信链路中存在反射面时的信道特性提出了一种分析方法,基于蒙特卡洛统计算法,通过发射大量光子,来模拟通信链路传输场景,并可以获得链路路径损耗和信道脉冲响应结果,从而可以考虑到光子与相关物体表面的相互作用,对估计无线光非视距通信链路的性能具有十分重要的意义。
技术实现思路
本专利技术提供了一种无线光非视距通信链路中存在反射面的信道特性分析方法,在所述无线光非视距通信链路模型中,接收机Rx位于坐标系原点,发射机Tx位于任意坐标(xt,yt,zt)处;同时假设Tx所发出的光束Beam与Rx的视场FOV均为一个圆锥,则Tx所发出的光束轴线指向的单位方向向量如式(1)所示:μT=[sinθTcosφT,sinθTsinφT,cosθT](1),Rx的视场轴线指向的单位方向向量如式(2)所示:μR=[sinθRcosφR,sinθRsinφR,cosθR](2),θT和θR表示光束轴线与视场轴线的倾斜角,角度方向为从z轴正方向到轴线;φT和φR表示光束轴线与视场轴线的方位角,角度方向为从x轴正方向逆时针到轴线在xy平面中的投影;βT和βR分别表示Tx光束的束散角和Rx视场的视场角;在空间中假设存在任意平面P,其法向量为nP=(xp,yp,zp),与z轴的交点Pz=(0,0,h),则平面的方程被写为P:xpx+ypy+zpz=zph;通过生成大量光子进行模拟光线传输的特性分析,从而获得无线光非视距通信链路中存在反射面时的信道特性包括通信链路路径损耗和任意时间段内的通信链路信道响应结果。所述无线光非视距通信链路中存在反射面时的信道特性分析方法包括以下步骤:步骤1、设置算法初始化的参数,包括系统模型参数:Tx的位置坐标(xt,yt,zt),θT,θR,φT,φR,βT,βR,平面P的法向量nP=(xp,yp,zp),以及其与z轴的交点Pz=(0,0,h);总共产生光子的个数Np;光子发生CIE的次数Nc;初始化记录光子每次发生CIE被接收机探测到的概率矩阵Pr的元素为0,该矩阵的尺寸为Np行,Nc列;初始化记录光子每次发生CIE被接收机探测到经过的路径矩阵Lr的元素为0,该矩阵的尺寸为Np行,Nc列;光子最低存活能量ξ;光子初始能量ξ0;光子的消光系数ke,散射系数ks,平面P的反射系数ρ;光源空间分布函数fLED(Ω),其中Ω为函数变量;光子经过散射事件后的空间分布函数fsca(Ω);光子经过平面反射事件后的空间分布函数fref(Ω);接收机的有效探测面积SR;步骤2、进入第一层循环,遍历所有Np个光子,初始循环变量为kp=1;步骤3、初始化本次循环光子的变量,包括:vr:用于记录光子的位置的变量,此时设置vr=(xt,yt,zt),初始化为发射机Tx的坐标,表示光子正在发射机出等待发送;μ:用于记录光子出射的方向的方向向量,初始化为μT,即发射机Tx的轴线的方向向量;d:光子经过的距离,初始化为0;ξ0=1,即归一化光子出射的能量;定义标志变量δ,用于指示记录光子改变方向的原因,δ为0表示从发射机出射,δ为1表示经过散射事件,δ为2表示经过反射事件,此处初始化δ=0;步骤4、进入第二层循环,遍历本次循环光子设定的发生CIE的次数,初始循环变量为kc=1;步骤5、根据δ的值,随机生成本次光子的出射向量μ,即:如果δ=0,则根据分布函数fLED(Ω)生成随机的出射向量μ;如果δ=1,则根据分布函数fsca(Ω)生成随机的出射向量μ;如果δ=2,则根据分布函数fref(Ω)生成随机的出射向量μ;步骤6、根据光子的消光系数ke随机生成本次光子的运动距离:dt;步骤7、计算新的光子的位置,即vr_new=vr+dμ;步骤8、通过光子前后两次位置,并结合平面P的方程xpx+ypy+zpz=zph,判断光子从vr运动到vr_new时是否穿过平面:如果穿过平面:计算该交点,并令vr_new等于该交点的坐标,设置δ=2;如果未穿过平面:设置δ=1;步骤9、更新变量vr=vr_new,d=d+dt,Lr(kp,kc)=d+||vr||,其中Lr(kp,kc)表示矩阵Lr的第kp行,第kc列的值;||·||表示求向量的模长;步骤10、计算光子经过CIE后,到达接收机的条件概率:如果δ=1,如果δ=2,其中,·表示向量的点乘,arccos表示反余弦函数;步骤11、如果表示光子在接收机视场内,则记录此时光子被成功接收的概率;如果δ=1,则其中Pr(kp,kc)表示矩阵Pr的第kp行,第kc列的值;如果δ=2,则其中Pr(kp,kc)表示矩阵Pr的第kp行,第kc列的值;步骤12、更新光子剩余能量:如果δ=1,如果δ=2,ξ0=ξ0(1-Prc)ρ;步骤13、进行kc=kc+1的操作,如果kc>Nc,则结束第二层循环,进入步骤14,否则返回步骤5;步骤14、进行kp=kp+1的操作,如果kp>Np,则结束第一层循环,进入步骤15,否则返回步骤3;步骤15、计算路径损耗PL,如下所示:路径损耗路径损耗PL的单位为dB,即分贝;步骤16、计算从光子从发射机出射时到记录结束时间tz内的信道脉冲响应,时间区间步长ts,包括以下过程:初始化信道脉冲响应向量CIR为1行,列的0向量,其中表示向上取整;计算光子传输时间矩阵其中c为光速;计算时间矩阵中元素对应脉冲响应向量CIR中的下标值矩阵,即通过循环遍历整个所述下标值矩阵,计算信道脉冲响应向量CIR,包括以下子步骤:子步骤S1:进入第一层循环,遍历所有Np个光子,初始循环变量为kp=1;子步骤S2:进入第二层循环,遍历发生CIE的次数Nc,初始循环变量为kc=1;子步骤S3:CIR(PphotoTimeIndex(kp,kc))=CIR(PphotoTimeIndex(kp,kc))+Pr(kp,kc);子步骤S4:进行kc=kc+1的操作,如果kc>Nc,则结束第二层循环,进入子步骤S5,否则进入子步骤S3;子步骤S4:进行kp=kp+1的操作,如果kp>Np,则结束第一层循环,否则进入子步骤S2;子步骤S5:对CIR进行归一化,即附图说明图1是无线光非视距通信链路中存在反射面本文档来自技高网...

【技术保护点】
1.一种无线光非视距通信链路中存在反射面的信道特性分析方法,其特征在于,在所述无线光非视距通信链路模型中,接收机Rx位于坐标系原点,发射机Tx位于任意坐标(x

【技术特征摘要】
1.一种无线光非视距通信链路中存在反射面的信道特性分析方法,其特征在于,在所述无线光非视距通信链路模型中,接收机Rx位于坐标系原点,发射机Tx位于任意坐标(xt,yt,zt)处;同时假设Tx所发出的光束Beam与Rx的视场FOV均为一个圆锥,则Tx所发出的光束轴线指向的单位方向向量如式(1)所示:
μT=[sinθTcosφT,sinθTsinφT,cosθT](1),
Rx的视场轴线指向的单位方向向量如式(2)所示:
μR=[sinθRcosφR,sinθRsinφR,cosθR](2),
θT和θR表示光束轴线与视场轴线的倾斜角,角度方向为从z轴正方向到轴线;φT和φR表示光束轴线与视场轴线的方位角,角度方向为从x轴正方向逆时针到轴线在xy平面中的投影;βT和βR分别表示Tx光束的束散角和Rx视场的视场角;
在空间中假设存在任意平面P,其法向量为nP=(xp,yp,zp),与z轴的交点Pz=(0,0,h),则平面的方程被写为P:xpx+ypy+zpz=zph;
通过生成大量光子进行模拟光线传输的特性分析,从而获得无线光非视距通信链路中存在反射面时的信道特性包括通信链路路径损耗和任意时间段内的通信链路信道响应结果,具体包括以下步骤:
步骤1、设置算法初始化的参数,包括系统模型参数:Tx的位置坐标(xt,yt,zt),θT,θR,φT,φR,βT,βR,平面P的法向量nP=(xp,yp,zp),以及其与z轴的交点Pz=(0,0,h);总共产生光子的个数Np;光子发生CIE的次数Nc;初始化记录光子每次发生CIE被接收机探测到的概率矩阵Pr的元素为0,该矩阵的尺寸为Np行,Nc列;初始化记录光子每次发生CIE被接收机探测到经过的路径矩阵Lr的元素为0,该矩阵的尺寸为Np行,Nc列;光子最低存活能量ξ;光子初始能量ξ0;光子的消光系数ke,散射系数ks,平面P的反射系数ρ;光源空间分布函数fLED(Ω),其中Ω为函数变量;光子经过散射事件后的空间分布函数fsca(Ω);光子经过平面反射事件后的空间分布函数fref(Ω);接收机的有效探测面积SR;
步骤2、进入第一层循环,遍历所有Np个光子,初始循环变量为kp=1;
步骤3、初始化本次循环光子的变量,包括:
vr:用于记录光子的位置的变量,此时设置vr=(xt,yt,zt),初始化为发射机Tx的坐标,表示光子正在发射机出等待发送;
μ:用于记录光子出射的方向的方向向量,初始化为μT,即发射机Tx的轴线的方向向量;
d:光子经过的距离,初始化为0;
ξ0=1,即归一化光子出射的能量;
定义标志变量δ,用于指示记录光子改变方向的原因,δ为0表示从发射机出射,δ为1表示经过散射事件,δ为2表示经过反射事件,此处初始化δ=0;
步骤4、进入第二层循环,遍历本次循环光子设定的发生CIE的次数,初始循环变量为kc=1;
步骤5、根据δ的值,随机生成本次光子的出射向量μ,即:
如果δ=0,则根据分布函数fLED(Ω)生成随机的出射向量μ;
如果δ=1,则根据分布函数...

【专利技术属性】
技术研发人员:曹天张洪明宋健
申请(专利权)人:清华大学
类型:发明
国别省市:北京;11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1