CN115276859A - 一种考虑高空平台物理结构的空地信道建模方法 - Google Patents
一种考虑高空平台物理结构的空地信道建模方法 Download PDFInfo
- Publication number
- CN115276859A CN115276859A CN202210527355.0A CN202210527355A CN115276859A CN 115276859 A CN115276859 A CN 115276859A CN 202210527355 A CN202210527355 A CN 202210527355A CN 115276859 A CN115276859 A CN 115276859A
- Authority
- CN
- China
- Prior art keywords
- line
- altitude platform
- sight
- receiving end
- ground
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000010363 phase shift Effects 0.000 claims abstract description 60
- 238000004891 communication Methods 0.000 claims abstract description 9
- 239000013598 vector Substances 0.000 claims description 107
- 238000004364 calculation method Methods 0.000 claims description 36
- 230000008859 change Effects 0.000 claims description 18
- 230000014509 gene expression Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000004083 survival effect Effects 0.000 claims description 6
- 230000008033 biological extinction Effects 0.000 claims description 4
- 235000007164 Oryza sativa Nutrition 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 235000009566 rice Nutrition 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims 1
- 230000033001 locomotion Effects 0.000 description 7
- 238000005562 fading Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L43/00—Arrangements for monitoring or testing data switching networks
- H04L43/08—Monitoring or testing based on specific metrics, e.g. QoS, energy consumption or environmental parameters
- H04L43/0852—Delays
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Environmental & Geological Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种考虑高空平台物理结构的空地信道建模方法,包括:针对考虑物理结构的高空平台对地通信场景,对第q根发射天线和第p根接收天线之间的信道进行建模,得到信道模型;实时计算信道模型中高空平台发射端和地面接收端非视距路径的有效数量;实时计算信道模型中视距路径和第n条非视距路径的时延;实时计算信道模型中高空平台发射端的离开角的俯仰角和方向角,及地面接收端的到达角的俯仰角和方向角;实时计算信道模型中视距路径的总相移和非视距路径的总相移。本发明中考虑到HAP与其天线阵列的空间关系,用天线位置代替传统发射端与接收端位置计算演进,支持时变的信道模型参数,并能够保证产生的信号角度的连续性。
Description
技术领域
本发明属于无线信息传输领域,具体涉及一种考虑高空平台物理结构的空地信道建模方法。
背景技术
高空平台(High Altitude Platform,HAP)以其易维护、高覆盖率等优点被广泛应用通信领域。实际场景中,高空平台自身会对信号传播造成遮挡,导致散射效应。此外,高空平台自身存在三维的姿态旋转,对高空平台的信道特性有显著影响。因此,为了更好地设计、优化和评估基于高空平台的无线通信系统,需要构建考虑高空平台物理结构影响的通信信道模型。
目前,非平稳信道建模方法大多考虑了发射端与接收端的移动,但是这些方法都将高空平台视作质点进行建模,忽略了高空平台自身物理结构给信道带来的影响。此外,大量空地通信实测数据表明,高空平台对地的信道是复杂多变的,而大部分建模方法只考虑了信道参数的时变性,没有考虑参数的连续动态演进。
有鉴于此,实有必要提供一种考虑高空平台物理结构的空地信道建模方法,该方法能够综合考虑高空平台的物理结构影响、三维姿态以及三维任意运动等因素,并对信道参数进行了连续动态演进,能够有效提高信道模型的真实性和准确性。
发明内容
针对于上述现有技术的不足,本发明的目的在于提供一种考虑高空平台物理结构的空地信道建模方法,以解决现有技术中将高空平台发射端视为质点导致的信道参数计算不精确及演进不连续的问题。
为达到上述目的,本发明采用的技术方案如下:
本发明的一种考虑高空平台物理结构的空地信道建模方法,步骤如下:
1)针对考虑物理结构的高空平台对地通信场景,对第q根发射天线和第p根接收天线之间的信道进行建模,得到信道模型;
2)实时计算信道模型中高空平台发射端和地面接收端非视距路径的有效数量;
3)实时计算信道模型中视距路径和第n条非视距路径的时延;
4)实时计算信道模型中高空平台发射端的离开角的俯仰角和方向角,及地面接收端的到达角的俯仰角和方向角;
5)实时计算信道模型中视距路径的总相移和非视距路径的总相移。
进一步地,所述步骤1)中信道模型的表达式为:
式中,hqp(t,τ)表示第q根发射天线和第p根接收天线间的信道冲激响应,t表示时间,τ表示时延,K表示莱斯因子(为视距路径功率与非视距路径功率的比值,参考值为8),e表示自然对数的底数,j表示虚数单位,表示视距路径的总相移,δ(·)表示单位冲激函数,表示视距路径的时延,NTx(t)和NRx(t)分别表示高空平台发射端和地面接收端的非视距路径数量,M表示每条非视距路径的散射支路数,表示第n条非视距路径的总相位,表示第n条非视距路径的时延。
其中,信号无遮挡地在发射端与接收端之间进行直线传播的路径称为视距路径,经过反射、散射及衍射后到达的信号的路径称为非视距路径,传播过程中导致信号传播方向发生改变的障碍物统称为散射体,如高空平台自身的物理结构,地面环境中的建筑物、植被等。
进一步地,所述步骤2)具体包括:
21)分别假设初始时刻高空平台发射端和地面接收端在t=t0时刻非视距路径的有效数量为及其中,分别表示高空平台发射端路径新生概率和消亡概率,参考值分别为0.9、0.3;和分别表示地面接收端路径新生概率和消亡概率,参考值分别为0.8、0.1;
其中,Δrp和Δrq分别为高空平台发射端与地面接收端的天线阵列间距,δd为由准平稳段确定的时间系数,和分别为时域和天线阵列下的场景相关系数;高空平台发射端和地面接收端经过Δt时刻非视距路径的有效数量为:
其中,E{·}表示平均值算符。
进一步地,所述步骤3)具体包括:
31)假设高空平台初始位置LTx(t0)、高空平台发射端第q根发射天线的初始位置为重构高空平台的物理模型,并在高空平台下方布置多个接收天线,利用射线跟踪方法,得到初始时刻高空平台发射端附近的散射体的位置
32)实时计算高空平台发射端和地面接收端的位置矢量,方法如下:
321)假设初始时刻高空平台发射端位置矢量为LTx(t0)和地面接收端位置矢量为LRx(t0);
322)迭代计算t+Δt时刻高空平台发射端和地面接收端的位置矢量,计算方法为:
式中,vTx(t)和vRx(t)分别表示高空平台发射端和地面接收端的速度矢量;
323)令t=t+Δt,返回步骤322),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
33)实时计算第q根发射天线和第p根接收天线的位置矢量,方法如下:
332)迭代计算t+Δt时刻第q根发射天线和第p根接收天线的位置矢量,计算方法为:
式中,RTx(t)表示高空平台发射端的旋转矢量,RRx(t)表示地面接收端的旋转矢量,RP(t)表示高空平台姿态矩阵,分别记为:
333)令t=t+Δt,返回步骤332),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
34)利用高空平台发射端和地面接收端的拓扑关系,实时计算视距路径的时延,方法如下:
341)计算初始时刻第q根发射天线和第p根接收天线之间的距离,即初始时刻视距路径的传播距离,计算方法如下:
其中,||·||代表一个向量的欧氏距离;
342)利用步骤33)中计算得到的第q根发射天线和第p根接收天线的位置矢量,迭代计算t+Δt时刻第q根发射天线和第p根接收天线之间的距离,即t+Δt时刻视距路径的传播距离为:
则视距路径t+Δt时刻的时延为:
式中,c表示光速;
343)令t=t+Δt,返回步骤342),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
35)实时计算第n条非视距路径的时延,方法如下:
351)初始时刻高空平台发射端和地面接收端附近第一个散射体的位置矢量为和由于高空平台发射端附近的散射体位置相对于高空平台所配备的天线位置是固定的,故第q根发射天线指向高空平台发射端附近第一个散射体的向量在任意时刻是固定的,记为:
初始时刻第p根发射天线指向地面接收端附近第一个散射体的向量记为:
初始时刻高空平台发射端附近第一个散射体指向地面接收端附近第一个散射体的向量记为:
则初始时刻第n条非视距信号的总传播距离记为:
其中,为接收端附近散射体的速度矢量;使用式(7)-(8)和式(17)-(18)中得到的位置矢量,计算经过Δt时刻的第p根发射天线指向地面接收端附近第一个散射体的向量ap,n,m(t+Δt),以及高空平台发射端附近第一个散射体指向地面接收端附近第一个散射体的向量cn,m(t+Δt),计算方法为:
然后利用步骤351)中的几何拓扑关系,实时计算第n条非视距路径的传播总距离为:
则每一条非视距路径的时延为:
式中,M为非视距路径的散射支路数;
353)令t=t+Δt,重复步骤352),直至高空平台发射端和地面接收端停止运动与姿态变化时结束。
进一步地,所述步骤4)具体包括:
41)实时计算高空平台发射端视距信号的离开角的方向角和俯仰角,及地面接收端视距信号的到达角和俯仰角,方法如下:
利用步骤34)中通过式(12)迭代计算得到的第q根发射天线和第p根接收天线之间的距离dLoS(t),实时计算高空平台发射端视距信号的离开角的方向角和俯仰角地面接收端视距信号的到达角的方向角和俯仰角计算方法为:
其中,arctan2{·}为反正切函数,ex、ey和ez分别为直角坐标系x轴、y轴和z轴的单位向量;
42)实时计算高空平台发射端非视距信号的离开角的方向角和俯仰角方法如下:由于第q根发射天线指向高空平台发射端附近第一个散射体的向量bq,n,m(t)在任意时刻是固定的,故任意时刻高空平台发射端非视距信号的离开角也是固定的,计算方法如下:
43)实时计算地面接收端非视距信号的到达角的方向角和俯仰角,方法如下:利用步骤35)中用式(19)迭代计算得到的向量ap,n,m(t),即第p根发射天线指向地面接收端附近第一个散射体的向量,实时计算地面接收端非视距信号的到达的方向角和俯仰角计算方法为:
进一步地,所述步骤5)具体包括:
将步骤4)中计算得到的高空平台发射端的离开角的方向角和俯仰角以及地面接收端的到达角的方向角和俯仰角带入式(32)和式(33),实时计算视距路径中离开信号的方向矢量和到达信号的方向矢量非视距路径中离开信号方向矢量和到达信号的方向矢量计算方法为:
其中,k=2πf0/c为波数,f0为用户输入的载波频率;
53)实时计算视距路径的空间相移和非视距路径的多普勒相移方法为:将步骤33)中通过式(7)和式(8)迭代计算得到的天线位置矢量和高空平台发射端和地面接收端的旋转矢量RTx(t)和RRx(t)、高空平台姿态矩阵RP(t),以及式(32)和式(33)计算得到的各信号的方向矢量,带入下式(36)和式(37);
54)计算视距路径的总相移和非视距路径的总相移计算方法为:将步骤51)计算得到的视距路径的随机初始相移和非视距路径的随机初始相移步骤52)计算得到的视距路径的多普勒相移和非视距路径的多普勒相移以及步骤53)计算得到的视距路径的空间相移和非视距路径的多普勒相移分别带入下式(38)和式(39);
本发明的有益效果:
本发明中将HAP的建模由质点转为实际物理模型的建模,该方法综合考虑了HAP的物理结构影响、三维姿态以及三维任意运动等因素,能够有效提高信道模型的真实性和准确性。
本发明中考虑到HAP与其天线阵列的空间关系,用天线位置代替传统发射端与接收端位置计算演进,支持时变的信道模型参数,并能够保证产生的信号角度的连续性。
附图说明
图1为本发明采用的建模方法的流程图。
图2为本发明中考虑高空平台物理结构情况下U2V通信场景示意图。
图3为本发明采用射线跟踪方法得到的散射体的位置示意图。
图4a为本发明产生的信号的离开角和到达角的方向角示意图。
图4b为本发明产生的信号的离开角和到达角的俯仰角示意图。
图5为本发明产生的信道衰落系数示意图。
具体实施方式
为了便于本领域技术人员的理解,下面结合实施例与附图对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。
参照图1-图3所示,本发明的一种考虑高空平台物理结构的空地信道建模方法,示例中具体如下:
1)用户输入包含高空平台的型号、材质、初始位置、速度、三维运动轨迹、三维姿态、运动方向的俯仰角和方位角,地面接收端的初始位置、速度大小、三维运动轨迹、运动方向的俯仰角和方位角,高空平台发射端及地面接收端的天线阵列型号大小、天线的初始位置、信道状态更新间隔、采样率、载波频率的参数信息;用户输入参数如表1所示,高空平台发射端与地面接收端轨迹与姿态变化如表2所示;
表1
表2
2)针对考虑物理结构的高空平台对地通信场景,对第q根发射天线和第p根接收天线之间的信道进行建模,得到信道模型,其表达式为:
式中,hqp(t,τ)表示第q根发射天线和第p根接收天线间的信道冲激响应,t表示时间,τ表示时延,K表示莱斯因子,为视距路径功率与非视距路径功率的比值,参考值为8,e表示自然对数的底数,j表示虚数单位,表示视距路径的总相移,δ(·)表示单位冲激函数,表示视距路径的时延,NTx(t)和NRx(t)分别表示高空平台发射端和地面接收端的非视距路径数量,M表示每条非视距路径的散射支路数,表示第n条非视距路径的总相位,表示第n条非视距路径的时延。
其中,信号无遮挡地在发射端与接收端之间进行直线传播的路径称为视距路径,经过反射、散射及衍射后到达的信号的路径称为非视距路径,传播过程中导致信号传播方向发生改变的障碍物统称为散射体,如高空平台自身的物理结构,地面环境中的建筑物、植被等。
3)实时计算信道模型中高空平台发射端和地面接收端非视距路径的有效数量;
31)分别假设初始时刻高空平台发射端和地面接收端在t=t0时刻非视距路径的有效数量为及其中,分别表示高空平台发射端路径新生概率和消亡概率,参考值分别为0.9、0.3;和分别表示地面接收端路径新生概率和消亡概率,参考值分别为0.8、0.1;
其中,Δrp和Δrq分别为高空平台发射端与地面接收端的天线阵列间距,δd为由准平稳段确定的时间系数,和分别为时域和天线阵列下的场景相关系数;高空平台发射端和地面接收端经过Δt时刻非视距路径的有效数量为:
其中,E{·}表示平均值算符。
4)实时计算信道模型中视距路径和第n条非视距路径的时延;
41)假设高空平台初始位置LTx(t0)、高空平台发射端第q根发射天线的初始位置为重构高空平台的物理模型,并在高空平台下方布置多个接收天线,(由于高空平台与其所配备的天线的相对位置是不变的,因此高空平台发射端附近的散射体位置相对于高空平台所配备的天线的位置是也是不变的)利用射线跟踪方法,得到初始时刻高空平台发射端附近的散射体的位置
42)实时计算高空平台发射端和地面接收端的位置矢量,方法如下:
421)假设初始时刻高空平台发射端位置矢量为LTx(t0)和地面接收端位置矢量为LRx(t0);
422)迭代计算t+Δt时刻高空平台发射端和地面接收端的位置矢量,计算方法为:
式中,vTx(t)和vRx(t)分别表示高空平台发射端和地面接收端的速度矢量;
423)令t=t+Δt,返回步骤422),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
43)实时计算第q根发射天线和第p根接收天线的位置矢量,方法如下:
432)迭代计算t+Δt时刻第q根发射天线和第p根接收天线的位置矢量,计算方法为:
式中,RTx(t)表示高空平台发射端的旋转矢量,RRx(t)表示地面接收端的旋转矢量,RP(t)表示高空平台姿态矩阵,分别记为:
433)令t=t+Δt,返回步骤432),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
44)利用高空平台发射端和地面接收端的拓扑关系,实时计算视距路径的时延,方法如下:
441)计算初始时刻第q根发射天线和第p根接收天线之间的距离,即初始时刻视距路径的传播距离,计算方法如下:
其中,||·||代表一个向量的欧氏距离;
442)利用步骤43)中计算得到的第q根发射天线和第p根接收天线的位置矢量,迭代计算t+Δt时刻第q根发射天线和第p根接收天线之间的距离,即t+Δt时刻视距路径的传播距离为:
则视距路径t+Δt时刻的时延为:
式中,c表示光速;
443)令t=t+Δt,返回步骤442),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
45)实时计算第n条非视距路径的时延,方法如下:
451)初始时刻高空平台发射端和地面接收端附近第一个散射体的位置矢量为和由于高空平台发射端附近的散射体位置相对于高空平台所配备的天线位置是固定的,故第q根发射天线指向高空平台发射端附近第一个散射体的向量在任意时刻是固定的,记为:
初始时刻第p根发射天线指向地面接收端附近第一个散射体的向量记为:
初始时刻高空平台发射端附近第一个散射体指向地面接收端附近第一个散射体的向量记为:
则初始时刻第n条非视距信号的总传播距离记为:
其中,为接收端附近散射体的速度矢量;使用式(7)-(8)和式(17)-(18)中得到的位置矢量,计算经过Δt时刻的第p根发射天线指向地面接收端附近第一个散射体的向量ap,n,m(t+Δt),以及高空平台发射端附近第一个散射体指向地面接收端附近第一个散射体的向量cn,m(t+Δt),计算方法为:
然后利用步骤451)中的几何拓扑关系,实时计算第n条非视距路径的传播总距离为:
则每一条非视距路径的时延为:
式中,M为非视距路径的散射支路数;
453)令t=t+Δt,重复步骤452),直至高空平台发射端和地面接收端停止运动与姿态变化时结束。
5)实时计算信道模型中高空平台发射端的离开角的俯仰角和方向角,及地面接收端的到达角的俯仰角和方向角;
51)实时计算高空平台发射端视距信号的离开角的方向角和俯仰角,及地面接收端视距信号的到达角和俯仰角,方法如下:
利用步骤44)中通过式(12)迭代计算得到的第q根发射天线和第p根接收天线之间的距离dLoS(t),实时计算高空平台发射端视距信号的离开角的方向角和俯仰角地面接收端视距信号的到达角的方向角和俯仰角计算方法为:
其中,arctan2{·}为反正切函数,ex、ey和ez分别为直角坐标系x轴、y轴和z轴的单位向量;
52)实时计算高空平台发射端非视距信号的离开角的方向角和俯仰角方法如下:由于第q根发射天线指向高空平台发射端附近第一个散射体的向量bq,n,m(t)在任意时刻是固定的,故任意时刻高空平台发射端非视距信号的离开角也是固定的,计算方法如下:
53)实时计算地面接收端非视距信号的到达角的方向角和俯仰角,方法如下:利用步骤45)中用式(19)迭代计算得到的向量ap,n,m(t),即第p根发射天线指向地面接收端附近第一个散射体的向量,实时计算地面接收端非视距信号的到达的方向角和俯仰角计算方法为:
将步骤5)中计算得到的高空平台发射端的离开角的方向角和俯仰角以及地面接收端的到达角的方向角和俯仰角带入式(32)和式(33),实时计算视距路径中离开信号的方向矢量和到达信号的方向矢量非视距路径中离开信号方向矢量和到达信号的方向矢量计算方法为:
其中,k=2πf0/c为波数,f0为用户输入的载波频率;
63)实时计算视距路径的空间相移和非视距路径的多普勒相移方法为:将步骤43)中通过式(7)和式(8)迭代计算得到的天线位置矢量和高空平台发射端和地面接收端的旋转矢量RTx(t)和RRx(t)、高空平台姿态矩阵RP(t),以及式(32)和式(33)计算得到的各信号的方向矢量,带入下式(36)和式(37);
64)计算视距路径的总相移和非视距路径的总相移计算方法为:将步骤61)计算得到的视距路径的随机初始相移和非视距路径的随机初始相移步骤62)计算得到的视距路径的多普勒相移和非视距路径的多普勒相移以及步骤63)计算得到的视距路径的空间相移和非视距路径的多普勒相移分别带入下式(38)和式(39);
7)将步骤2)中计算得到的高空平台发射端和地面接收端非视距路径的有效数量NTx(t)和NRx(t)、步骤3)中计算得到的视距路径和第n条非视距路径的时延τLoS(t)和以及步骤5)中得到的视距路径的总相移和非视距路径的总相移输入信道模型,实时计算输出非平稳信道衰落。
本实例得到的效果可以通过仿真实验中所获得的具体数据进行进一步说明。本实例的效果分2部分展示,高空平台物理结构影响通过图4a、图4b说明,三维姿态的影响通过图5说明。从图4a、图4b、图5中可观察到:1)图4a给出了信号离开角和到达角的方向角,在方向角的维度中,视距路径与非视距路径的到达角都是时变的,由模型生成的随机地面散射环境决定,视距路径的离开角也是时变的,由高空平台的姿态、位置变化决定,而非视距路径的离开角则是固定的,由高空平台的自身物理结构决定;图4b给出了信号离开角和到达角的俯仰角,在俯仰角的维度中,视距路径的离开角与到达角几乎不发生变化,这是由于高空平台距地面的距离较远导致,非视距路径的离开角是固定的,由高空平台的自身物理结构决定,仅有非视距路径的到达角是时变的,由模型生成的随机地面散射环境决定;2)图5给出了本发明产生的信道衰落系数,可以看到,本发明提出模型产生的信道衰落系数随着姿态角度的变化而变化,而未考虑姿态的信道衰落系数则无法体现这一特性。此外,本发明提出的模型产生的信道衰落系数相较于未考虑姿态的情况起伏变化更加频繁,这是由于姿态的变化将导致信道中信号传播复杂性的增加。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (6)
1.一种考虑高空平台物理结构的空地信道建模方法,其特征在于,步骤如下:
1)针对考虑物理结构的高空平台对地通信场景,对第q根发射天线和第p根接收天线之间的信道进行建模,得到信道模型;
2)实时计算信道模型中高空平台发射端和地面接收端非视距路径的有效数量;
3)实时计算信道模型中视距路径和第n条非视距路径的时延;
4)实时计算信道模型中高空平台发射端的离开角的俯仰角和方向角,及地面接收端的到达角的俯仰角和方向角;
5)实时计算信道模型中视距路径的总相移和非视距路径的总相移。
4.根据权利要求3所述的考虑高空平台物理结构的空地信道建模方法,其特征在于,所述步骤3)具体包括:
31)假设高空平台初始位置LTx(t0)、高空平台发射端第q根发射天线的初始位置为重构高空平台的物理模型,并在高空平台下方布置多个接收天线,利用射线跟踪方法,得到初始时刻高空平台发射端附近的散射体的位置
32)实时计算高空平台发射端和地面接收端的位置矢量,方法如下:
321)假设初始时刻高空平台发射端位置矢量为LTx(t0)和地面接收端位置矢量为LRx(t0);
322)迭代计算t+Δt时刻高空平台发射端和地面接收端的位置矢量,计算方法为:
式中,vTx(t)和vRx(t)分别表示高空平台发射端和地面接收端的速度矢量;
323)令t=t+Δt,返回步骤322),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
33)实时计算第q根发射天线和第p根接收天线的位置矢量,方法如下:
332)迭代计算t+Δt时刻第q根发射天线和第p根接收天线的位置矢量,计算方法为:
式中,RTx(t)表示高空平台发射端的旋转矢量,RRx(t)表示地面接收端的旋转矢量,RP(t)表示高空平台姿态矩阵,分别记为:
333)令t=t+Δt,返回步骤332),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
34)利用高空平台发射端和地面接收端的拓扑关系,实时计算视距路径的时延,方法如下:
341)计算初始时刻第q根发射天线和第p根接收天线之间的距离,即初始时刻视距路径的传播距离,计算方法如下:
式中,||·||代表一个向量的欧氏距离;
342)利用步骤33)中计算得到的第q根发射天线和第p根接收天线的位置矢量,迭代计算t+Δt时刻第q根发射天线和第p根接收天线之间的距离,即t+Δt时刻视距路径的传播距离为:
则视距路径t+Δt时刻的时延为:
式中,c表示光速;
343)令t=t+Δt,返回步骤342),直至高空平台发射端和地面接收端停止运动与姿态变化时结束;
35)实时计算第n条非视距路径的时延,方法如下:
351)初始时刻高空平台发射端和地面接收端附近第一个散射体的位置矢量为和由于高空平台发射端附近的散射体位置相对于高空平台所配备的天线位置是固定的,故第q根发射天线指向高空平台发射端附近第一个散射体的向量在任意时刻是固定的,记为:
初始时刻第p根发射天线指向地面接收端附近第一个散射体的向量记为:
初始时刻高空平台发射端附近第一个散射体指向地面接收端附近第一个散射体的向量记为:
则初始时刻第n条非视距信号的总传播距离记为:
式中,为接收端附近散射体的速度矢量;使用式(7)-(8)和式(17)-(18)中得到的位置矢量,计算经过Δt时刻的第p根发射天线指向地面接收端附近第一个散射体的向量ap,n,m(t+Δt),以及高空平台发射端附近第一个散射体指向地面接收端附近第一个散射体的向量cn,m(t+Δt),计算方法为:
然后利用步骤351)中的几何拓扑关系,实时计算第n条非视距路径的传播总距离为:
则每一条非视距路径的时延为:
式中,M为非视距路径的散射支路数;
353)令t=t+Δt,重复步骤352),直至高空平台发射端和地面接收端停止运动与姿态变化时结束。
5.根据权利要求4所述的考虑高空平台物理结构的空地信道建模方法,其特征在于,所述步骤4)具体包括:
41)实时计算高空平台发射端视距信号的离开角的方向角和俯仰角,及地面接收端视距信号的到达角和俯仰角,方法如下:
利用步骤34)中通过式(12)迭代计算得到的第q根发射天线和第p根接收天线之间的距离dLoS(t),实时计算高空平台发射端视距信号的离开角的方向角和俯仰角地面接收端视距信号的到达角的方向角和俯仰角计算方法为:
其中,arctan2{·}为反正切函数,ex、ey和ez分别为直角坐标系x轴、y轴和z轴的单位向量;
43)实时计算地面接收端非视距信号的到达角的方向角和俯仰角,方法如下:利用步骤35)中用式(19)迭代计算得到的向量ap,n,m(t),即第p根发射天线指向地面接收端附近第一个散射体的向量,实时计算地面接收端非视距信号的到达的方向角和俯仰角计算方法为:
6.根据权利要求5所述的考虑高空平台物理结构的空地信道建模方法,其特征在于,所述步骤5)具体包括:
将步骤4)中计算得到的高空平台发射端的离开角的方向角和俯仰角以及地面接收端的到达角的方向角和俯仰角带入式(32)和式(33),实时计算视距路径中离开信号的方向矢量和到达信号的方向矢量非视距路径中离开信号方向矢量和到达信号的方向矢量计算方法为:
式中,k=2πf0/c为波数,f0为用户输入的载波频率;
53)实时计算视距路径的空间相移和非视距路径的多普勒相移方法为:将步骤33)中通过式(7)和式(8)迭代计算得到的天线位置矢量和高空平台发射端和地面接收端的旋转矢量RTx(t)和RRx(t)、高空平台姿态矩阵RP(t),以及式(32)和式(33)计算得到的各信号的方向矢量,带入下式(36)和式(37);
54)计算视距路径的总相移和非视距路径的总相移计算方法为:将步骤51)计算得到的视距路径的随机初始相移和非视距路径的随机初始相移步骤52)计算得到的视距路径的多普勒相移和非视距路径的多普勒相移以及步骤53)计算得到的视距路径的空间相移和非视距路径的多普勒相移分别带入下式(38)和式(39);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210527355.0A CN115276859A (zh) | 2022-05-16 | 2022-05-16 | 一种考虑高空平台物理结构的空地信道建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210527355.0A CN115276859A (zh) | 2022-05-16 | 2022-05-16 | 一种考虑高空平台物理结构的空地信道建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115276859A true CN115276859A (zh) | 2022-11-01 |
Family
ID=83759300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210527355.0A Pending CN115276859A (zh) | 2022-05-16 | 2022-05-16 | 一种考虑高空平台物理结构的空地信道建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115276859A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108418645A (zh) * | 2018-01-26 | 2018-08-17 | 南京航空航天大学 | 一种非平稳移动通信信道建模及参数平滑演进方法 |
CN108768472A (zh) * | 2018-03-19 | 2018-11-06 | 西安电子科技大学 | 一种临近空间高空平台空间-极化分集mimo信道建模方法 |
US20220006674A1 (en) * | 2018-10-26 | 2022-01-06 | Fondation B-Com | Method for estimating a wireless communication channel, device for estimating a wireless communication channel and computer program therefor |
-
2022
- 2022-05-16 CN CN202210527355.0A patent/CN115276859A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108418645A (zh) * | 2018-01-26 | 2018-08-17 | 南京航空航天大学 | 一种非平稳移动通信信道建模及参数平滑演进方法 |
CN108768472A (zh) * | 2018-03-19 | 2018-11-06 | 西安电子科技大学 | 一种临近空间高空平台空间-极化分集mimo信道建模方法 |
US20220006674A1 (en) * | 2018-10-26 | 2022-01-06 | Fondation B-Com | Method for estimating a wireless communication channel, device for estimating a wireless communication channel and computer program therefor |
Non-Patent Citations (1)
Title |
---|
邹劲柏;谢浩;艾渤;郭爱煌;: "高速移动环境大规模MIMO信道建模与性能分析", 铁道学报, no. 04, 15 April 2018 (2018-04-15) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109412673B (zh) | 面向无人机通信的几何随机信道模型的实时仿真方法 | |
CN109471080A (zh) | 基于simulink的高速平台雷达回波信号模拟系统 | |
CN112070894B (zh) | 真实环境导航多径实时仿真方法、装置、介质及电子设备 | |
CN107479038B (zh) | 一种高精度雷达目标回波实时模拟方法 | |
CN106291489B (zh) | 适用于多种发射信号波形的合成孔径雷达回波仿真方法 | |
CN108418645B (zh) | 一种非平稳移动通信信道建模及参数平滑演进方法 | |
CN110850376B (zh) | 箔条干扰模拟方法及箔条干扰模拟器 | |
CN105445718B (zh) | 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法 | |
CN105978647B (zh) | 一种三维双移动mimo传播信道的建模及仿真方法 | |
CN106054187B (zh) | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 | |
CN108459308A (zh) | 一种基于时变rcs数据的回波模拟方法及装置 | |
CN111507008B (zh) | 一种基于参数化的通用雷达模型建模方法 | |
CN103941250B (zh) | 一种被动声纳阵列信号仿真中的海洋信道仿真方法 | |
CN206058609U (zh) | 一种雷达电磁环境模拟系统 | |
CN108872971A (zh) | 一种基于运动单阵列的目标定位方法与装置 | |
CN114124266A (zh) | 一种基于irs辅助无人机与无人船通信的信道建模方法 | |
CN107092017A (zh) | 一种基于近场mimo成像的rcs外推方法 | |
CN108169746A (zh) | 线性调频脉冲半主动雷达导引头信号处理方法 | |
CN111181659B (zh) | 基于5g的无线信道仿真器的多普勒效应实现方法 | |
Deng et al. | A novel PE/FDTD hybrid model for predicting echo signals of radar targets in large-scale complex environments | |
CN114039684A (zh) | 一种基于几何的宽带uav mimo信道建模方法 | |
CN105548959B (zh) | 一种基于稀疏重建的多传感器多目标的定位方法 | |
CN112422214B (zh) | 一种适用于航空信道的通信效果演示验证系统 | |
CN103869297B (zh) | 高超声速平台载雷达的杂波生成系统及其方法 | |
CN115276859A (zh) | 一种考虑高空平台物理结构的空地信道建模方法 |
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 |