1.一种基于规则几何的UAV‑V非平稳信道建模方法,其特征在于,包括以下步骤:步骤1,针对城市街道场景,仿真建立无人机到车非平稳信道的三维半椭球、圆柱、单球复合几何模型,使用圆柱模型表征无人机及移动发射端Tx侧散射体分布,用单球模型表征移动车辆及移动接收端Rx侧散射体分布,用半椭球模型模拟城市街道环境,并且使用底面同心圆模型对地面一次反射分布进行了模拟;
步骤2,根据移动到移动信道特征,建立信道模型,并计算脉冲信道响应,包括视距LoS分量与非视距NLoS分量;
步骤3,基于几何的随机模型,推导几何模型空时互相关函数和多普勒功率谱密度的表达式,建立基于规则几何的UAV‑V非平稳信道模型;
步骤2包括:对于多输入多输出MIMO信道,计算脉冲信道响应,即移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之间的信道关系,p=1,2,…,MT,q==1,2,…,MR;其中MT表示移动发射端Tx天线总数,MR表示移动接收端Rx天线总数;
所述脉冲信道响应包括视距LoS分量和非视距NLoS分量,其中非视距NLoS分量根据信号从移动发射端Tx到移动接收端Rx经历的过程分为单弹SB分量,双弹DB分量和地面一次反射SBG分量,具体为经过移动发射端Tx圆柱模型上的散射体的SB1射线分量,经过移动接收端Rx单球模型上的散射体的SB2射线分量,经过半椭球体模型上散射体的SB3射线分量,以及分别经过移动发射端Tx圆柱模型上的散射体和移动接收端Rx单球模型上的散射体的DB1射线分量,分别经过移动发射端Tx圆柱模型上的散射体和半椭球体模型上的散射体的DB2射线分量,分别经过移动接收端Rx单球模型上的散射体和半椭球体模型上的散射体的DB2射线分量,以及经过地面反射的SBG射线分量;
多输入多输出MIMO信道的第p个移动发射端Tx天线单元与第q个移动接收端Rx天线单元之间的复杂信道冲激响应hpq(t)表示为:其中, 为t时刻视距LoS分量, 为t时刻单弹分量, 为t时刻双弹分量,是t时刻地面一次反射分量;
步骤2还包括:所述视距LoS分量 的具体表达式为:
其中K指莱斯因子,λ指载波波长,εpq(t)指第p个移动发射端Tx天线单元到第q个移动接收端Rx天线单元之间的距离,fTm,fRm分别为移动发射端经过静态散射体导致的最大多普勒频率和移动接收端Rx的最大多普勒频率, 分别为移动发射端视距径的方位角和俯仰角, 分别为接收端视距径的方位角和俯仰角,△t代表时间延迟;γT表示移动发射端Tx的移动方位角,ξ表示移动发射端Tx的移动俯仰角,γR表示移动接收端Rx的移动方位角,e是自然常数,j是虚数单位;
参数ΔT是第p根天线元件和UAV天线阵列中心中间的距离,参数ΔR是第q根GS端天线阵列中心之间的距离,对于均匀性阵列,定义为:其中δT、δR分别为移动发射端Tx的天线阵列间距和移动接收端Rx的天线阵列间距;
移动发射端Tx和移动接收端Rx采用均匀线性天线阵列且分别配有MT和MR根子天线,OT和OR分别表示移动发射端Tx天线阵列的中心和移动接收端Rx天线阵列的中心,设定圆柱模型上分布着N1个动态散射体,第n1个散射体表示Sn1,单球模型上分布着N2个动态散射体,第n2个散射体表示Sn2,其中n1=1,2,…,N1;n2=1,2,…,N2;
半椭球模型上分布着N3个静止散射体,第n3个散射体表示Sn3,以及表示地面反射的圆形模型上分布着N4个动态散射体,第n4个散射体表示Sn4,其中n3=1,2,…,N3;n4=1,2,…,N4;
分别由 和 表示经过散射体Sni射线的离开方位角AAoD和离开俯仰角EAoD,分别用 和 表示相应的到达方位角AAoA和到达俯仰角EAoA;i=1,2,3,4;
对于单弹分量 表达式为:
其中 指各单弹分量的功率,εpq,ni(t)指经过散射体Sni的单弹路径的路径总和,ni(t)表示t时刻第ni个散射体, 表示散射体个数趋于无穷大,fTm,n,fRm,n分别为移动发射端Tx经过动态散射体导致的最大多普勒频率和移动接收端Rx的最大多普勒频率,所述路径总和包括三条:经过移动发射端Tx圆柱模型上的散射体的SB1射线,经过移动接收端Rx单球模型上的散射体的SB2射线,经过半椭球体模型上散射体的SB3射线;
分别为移动发射端到散射点的方位角和俯仰角, 分别为移动接收端
到散射点的方位角和俯仰角;
对于双弹分量,共有三条双弹路径,分别是:分别经过移动发射端Tx圆柱模型上的散射体和移动接收端Rx单球模型上的散射体的DB1射线 分别经过移动发射端Tx圆柱模型上的散射体和半椭球体模型上的散射体的DB2射线 分别经过移动接收端Rx单球模型上的散射体和半椭球体模型上的散射体的DB2射线 各自的表达式为:其中 表示DB1射线分量的能量相关参数, 表示DB2分量的能量相关参数, 表示DB3分量的能量相关参数, 表示散射体Sn1与散射体Sn2之间的距离, 表示经过散射体Sn1射线的离开方位角, 表示经过散射体Sn1射线的离开俯仰角, 表示经过散射体Sn2射线的到达方位角, 表示经过散射体Sn2射线的到达俯仰角;
表示散射体Sn1与散射体Sn3之间的距离, 表示经过散射体Sn1射线的离开方位角, 表示经过散射体Sn3射线的到达方位角, 表示经过散射体Sn3射线的到达俯仰角; 表示散射体Sn1与散射体Sn3之间的距离, 表示经过散射体Sn3射线的离开方位角, 经过散射体Sn3射线的离开俯仰角;
对于地面一次反射分量 表达式为:
其中ηSBG表示表示SBG射线分量的能量相关参数,εpq,n4(t)表示经过散射体Sn4的路径总和, 表示经过散射体Sn4射线的离开方位角, 经过散射体Sn4射线的离开俯仰角, 表示经过散射体Sn4射线的到达方位角, 表示经过散射体Sn4射线的到达俯仰角;
步骤2还包括:由于移动散射体的存在,相应的最大多普勒频移也会带来改变,fTm,n,fRm,n计算公式为:fTm,n=|vTcosγT‑vSR|/λ,
fRm,n=|vRcosγR‑vSR|/λ,
其中vT表示移动发射端Tx的移动速度,γT表示移动发射端Tx的移动方位角,vSR表示散射体的移动速度;
vR表示移动接收端Rx的移动速度,γR表示移动接收端Rx的移动方位角;
步骤2还包括:方位角AoD和俯仰角AoA以及地面一次反射模型圆的半径R的离散变量形式转换为连续变量,同时设定方位角AoD、俯仰角AoA和地面一次反射模型圆的半径R是独立的;使用冯米塞斯分布来描述经过散射体Sni射线的离开方位角 和到达方位角 的分布,定义为:其中 表示经过散射体Sni射线的离开方位角, 表示经过散射体Sni射线的离开方位角 的概率密度函数PDF, 表示离开方位角 的均值,ki是控制方位角相对于平均值分布浓度的实值参数,ki≥0,ki越小,散射体分布越趋近于均匀分布,ki越大,则越趋近于标准正态分布,I0(ki)表示0阶第一类修正贝塞尔函数;
其中 表示经过散射体Sni射线的到达方位角, 表示经过散射体Sni射线的到达方位角 的概率密度函数PDF, 表示到达方位角 的均值;
对于俯仰角的概率密度函数PDF,使用车对车V2V场景下常用的余弦分布来描述俯仰角的分布情况,具体表达式为:其中 表示经过散射体Sni射线的离开俯仰角, 表示经过散射体Sni射线的离开俯仰角 的分布函数, 表示离开俯仰角 的均值, 表示离开俯仰角 偏离均值的最大范围;
其中 表示经过散射体Sni射线的到达俯仰角, 表示经过散射体Sni射线的到达俯仰角 的分布函数, 表示到达俯仰角 的均值, 表示到达俯仰角 偏离均值的最大范围;
表征地面反射经过的散射体分布位置的半径R的概率密度函数f(R)为:
2
f(R)=2R/(RR) ,
其中RR表示围绕移动接收端Rx的球体模型的半径。
2.根据权利要求1所述的方法,其特征在于,步骤3包括:空时互相关函数表示为:其中, 表示移动发射端Tx的第p根天线与移动接收端Rx的
第q根天线之间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道的空时互相关函数, 表示空时互相关函数的视距LoS分量,表示空时互相关函数的单弹分量, 表示空时
互相关函数的双弹分量, 表示空时互相关函数的地面一次反射分
量,δT,δR分别表示移动发射端Tx的天线阵列间距和移动接收端Rx的天线阵列间距;
通过设置Δt=0得到空间互相关函数CCF,表示为:
其中 表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线
之间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道的空间互相关函数, 表 示空间互相关函数的视距LoS 分量 ,表示空间互相关函数的单弹分量, 表示空间互
相关函数的双弹分量, 表示空间互相关函数的地面一次反射分量;
通过设置δT=δR=0,即在移动发射端Tx和移动接收端Rx使用相同的天线单元,得到时间自相关函数ACF来衡量信道的时间相关性,表示为:其中 表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之
间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道的时间自相关函数, 表示时间自相关函数的视距LoS分量, 表示时间自相关函数的单弹分量, 表示时间自相关函数的双弹分量,
表示时间自相关函数的地面一次反射分量;
其中,对于视距LoS分量,具体表达式为:
其中 表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之
间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道的时间自相关函数视距LoS分量,εp′q′(t)表示移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的距离,εpq(t)表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之间的距离;
对于单弹信道中经过移动发射端Tx圆柱模型上的散射体的SB1射线分量,经过移动接收端Rx圆柱模型上的散射体的SB2射线分量,经过半椭球体模型上散射体的SB3射线分量,具体表达式为:其中 表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之
间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道经过移动发射端Tx圆柱模型上的散射体的时间自相关函数SB1射线分量,ηSB1表示SB1射线分量的能量相关参数, 表示经过散射体Sn1射线的离开俯仰角的最小值, 表示经过散射体Sn1射线的离开俯仰角的最大值, 表示移动发射端Tx的第p根天线与散射体Sn1的距离,表示移动发射端Tx的第p’根天线与散射体Sn1的距离, 表示移动接收端Rx的第q根天线与散射体Sn1的距离, 表示移动接收端Rx的第q’根天线与散射体Sn1的距离,表示移动发射端Tx经过散射体Sn1导致的最大多普勒频率, 表示移动接收端Rx经过散射体Sn1导致的最大多普勒频率, 表示经过散射体Sn1射线的离开方位角, 表示经过散射体Sn1射线的离开俯仰角, 表示经过散射体Sn1射线的到达方位角, 表示经过散射体Sn1射线的到达俯仰角, 表示 的概率密度函数, 表示 的概率密度函数;
对于单弹信道中经过移动接收端Rx单球模型上的散射体的SB2射线分量,具体表达式为:其中 表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之
间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道经过移动接收端Rx单球模型上的散射体的时间自相关函数SB2射线分量,ηSB2表示SB2射线分量的能量相关参数, 表示经过散射体Sn2射线的到达俯仰角的最小值, 表示经过散射体Sn2射线的到达俯仰角的最大值, 表示移动发射端Tx的第p根天线与散射体Sn2的距离,表示移动发射端Tx的第p’根天线与散射体Sn2的距离, 表示移动接收端Rx的第q根天线与散射体Sn2的距离, 表示移动接收端Rx的第q’根天线与散射体Sn2的距离,表示移动发射端Tx经过散射体Sn2导致的最大多普勒频率, 表示移动接收端Rx经过散射体Sn2导致的最大多普勒频率, 表示经过散射体Sn2射线的离开方位角, 表示经过散射体Sn2射线的离开俯仰角, 表示经过散射体Sn2射线的到达方位角, 表示经过散射体Sn2射线的到达俯仰角, 表示 的概率密度函数, 表示 的概率密度函数;
对于经过半椭球体模型上散射体的SB3射线分量,具体表达式为:
其中 表示为移动发射端Tx的第p根天线与移动接收端Rx的第q根天线
之间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道经过半椭球体模型上散射体的时间自相关函数SB3射线分量,ηSB3表示SB3射线分量的能量相关参数, 表示经过散射体Sn3射线的到达俯仰角的最小值, 表示经过散射体Sn3射线的到达俯仰角的最大值, 表示移动发射端Tx的第p根天线与散射体Sn3的距离, 表示移动发射端Tx的第p’根天线与散射体Sn3的距离, 表示移动接收端Rx的第q根天线与散射体Sn3的距离, 表示移动接收端Rx的第q’根天线与散射体Sn3的距离, 表示移动发射端Tx经过散射体Sn3导致的最大多普勒频率, 表示移动接收端Rx经过散射体Sn3导致的最大多普勒频率, 表示经过散射体Sn3射线的离开方位角, 表示经过散射体Sn3射线的离开俯仰角, 表示经过散射体Sn3射线的到达方位角, 表示经过散射体Sn3射线的到达俯仰角, 表示 的概率密度函数, 表示 的概率密度函数;
对于经过地面反射的SBG射线分量,具体表达式为:
其中, 表示为移动发射端Tx的第p根天线与移动接收端Rx的第q根天
线之间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道经过地面反射的时间自相关函数SBG射线分量,ηSB4表示SBG射线分量的能量相关参数,表示移动发射端Tx的第p根天线与散射体Sn4的距离, 表示移动发射端Tx的第p’根天线与散射体Sn4的距离, 表示移动接收端Rx的第q根天线与散射体Sn4的距离, 表示移动接收端Rx的第q’根天线与散射体Sn4的距离, 表示移动发射端Tx经过散射体Sn4导致的最大多普勒频率, 表示移动接收端Rx经过散射体Sn4导致的最大多普勒频率,表示经过散射体Sn4射线的离开方位角, 表示经过散射体Sn4射线的离开俯仰角, 表示经过散射体Sn4射线的到达方位角,H0表示移动发射端Tx和移动接收端Rx之间的高度差,表示 的概率密度函数;
对于双弹信道中经过移动发射端Tx圆柱模型上的散射体和经过移动接收端Rx圆柱模型上的散射体的DB1射线分量,具体表达式为:其中, 表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线
之间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道经过移动发射端Tx圆柱模型上的散射体和经过移动接收端Rx圆柱模型上的散射体的时间自相关函数DB1射线分量, 表示经过散射体Sn1射线的离开俯仰角的最小值, 表示经过散射体Sn1射线的离开俯仰角的最大值, 表示经过散射体Sn2射线的到达俯仰角的最小值,表示经过散射体Sn2射线的到达俯仰角的最大值, 表示天线经过散射体Sn1与散射体Sn2之间的距离;
对于双弹信道中经过移动发射端Tx圆柱模型上的散射体和经过半椭球体模型上散射体的DB2射线分量,具体表达式为:其中, 表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线
之间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道经过移动发射端Tx圆柱模型上的散射体和经过半椭球体模型上散射体的时间自相关函数DB2射线分量, 表示经过散射体Sn1射线的离开俯仰角的最小值, 表示经过散射体Sn1射线的离开俯仰角的最大值, 表示经过散射体Sn3射线的到达俯仰角的最小值, 表示经过散射体Sn3射线的到达俯仰角的最大值, 表示天线经过散射体Sn1与散射体Sn3之间的距离;
对于双弹信道中经过移动接收端Rx圆柱模型上的散射体和经过半椭球体模型上散射体的DB3射线分量,具体表达式为:其中, 表示为移动发射端Tx的第p根天线与移动接收端Rx的第q根天
线之间的信道和移动发射端Tx的第p’根天线与移动接收端Rx的第q’根天线之间的信道经过移动接收端Rx圆柱模型上的散射体和经过半椭球体模型上散射体的时间自相关函数DB3射线分量, 表示经过散射体Sn3射线的到达俯仰角的最小值, 表示经过散射体Sn3射线的到达俯仰角的最大值, 表示经过散射体Sn2射线的到达俯仰角的最小值, 表示经过散射体Sn2射线的到达俯仰角的最大值, 表示天线经过散射体Sn2与散射体Sn3之间的距离;
能量参数需满足如下公式:
ηSB1+ηSB2+ηSB3+ηDB1+ηDB2+ηDB3+ηSBG=1。
3.根据权利要求2所述的方法,其特征在于,步骤3还包括:多普勒功率谱密度表达式为:其中,S(f,t)表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之间的信道的多普勒功率谱密度,Rp,q(Δt,t)表示移动发射端Tx的第p根天线与移动接收端Rx的第q根天线之间的信道的时间自相关函数。
4.一种电子设备,其特征在于,包括处理器和存储器,所述存储器存储有程序代码,当所述程序代码被所述处理器执行时,使得所述处理器执行如权利要求1至3中任一项所述的方法的步骤。
5.一种存储介质,其特征在于,存储有计算机程序或指令,当所述计算机程序或指令在计算机上运行时,执行如权利要求1至3中任一项所述的方法的步骤。