CN113364544B - 一种非对称毫米波亚毫米波无线信道仿真方法 - Google Patents

一种非对称毫米波亚毫米波无线信道仿真方法 Download PDF

Info

Publication number
CN113364544B
CN113364544B CN202110635458.4A CN202110635458A CN113364544B CN 113364544 B CN113364544 B CN 113364544B CN 202110635458 A CN202110635458 A CN 202110635458A CN 113364544 B CN113364544 B CN 113364544B
Authority
CN
China
Prior art keywords
scatterer
hop
mobile station
station side
channel
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.)
Active
Application number
CN202110635458.4A
Other languages
English (en)
Other versions
CN113364544A (zh
Inventor
张焱
杜川
张万成
何遵文
赵雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202110635458.4A priority Critical patent/CN113364544B/zh
Publication of CN113364544A publication Critical patent/CN113364544A/zh
Application granted granted Critical
Publication of CN113364544B publication Critical patent/CN113364544B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种非对称毫米波亚毫米波室外信道仿真方法，涉及无线通信领域。本发明充分考虑了信道分布环境呈现出非对称性的特点，设计了一种毫米波亚毫米波信道仿真方法，该方法能够生成非对称条件下的毫米波亚毫米波无线信道冲激响应和参数。通过结合毫米波亚毫米波具有空气分子吸收和粗糙表面漫散射的传播特点，将信道冲激响应分为视距径、单跳和双跳，从而准确地实现上下行信道的参数刻画与信道模拟。并且，非对称上下行信道的相关性得到了充分利用，这有效地降低了信道响应和信道参数生成的复杂度，从而增强了其在系统仿真验证应用中的实用性。本发明可应用于非对称毫米波亚毫米波通信系统的设计和部署中，并提供理论与模型依据。

Description

一种非对称毫米波亚毫米波无线信道仿真方法
技术领域
本发明涉及无线通信领域，具体涉及一种非对称的毫米波和亚毫米波频段的无线信道仿真方法。
背景技术
随着信息社会的不断发展，未来的后5G和6G移动通信系统将需要在数据速率、时延、覆盖能力等多个方面达到更高的要求。当前频率较低的频段资源几乎已经被各种通信业务占用而变得拥挤不堪，难以挖掘出可加以利用的频谱资源。而毫米波和亚毫米波频段具有超大的、连续可利用的带宽，能满足通信系统上百Gbps的数据传输速率和百微秒级的时延需求。且毫米波和亚毫米波波长短，使得在基站侧或用户侧都可以在有限的空间内部署大规模的多输入多输出(Multiple-Input-Multiple-Output，MIMO)天线阵列。当前的大规模MIMO移动通信中，绝大多数基站侧和移动台侧的多波束天线阵列采用对称结构进行设计。在对称结构中，同一侧的发射和接收通道数量相等，与之对应的接收和发送波束宽度以及增益也相同。但对称的收发阵列结构存在计算复杂度大，硬件成本高，能量功率消耗大等缺点，从而限制了高频通信的大范围部署和应用。
在非对称结构的阵列中，同一侧(基站侧和移动台侧)用于发送和接收的阵列规模不相等。采用较大规模的高增益波束发送阵列和较小规模的低增益波束接收阵列，在获得更高的系统容量和链路增益的同时，收发端也能更容易实现波束的动态跟踪和管理。在收发阵列不对称的情况下，即使是时分双工(Time-Division Duplex，TDD)系统，也不再满足上下行信道的互易性。但本质上，信道冲激响应由传播环境、收发波束宽度等因素决定，相似的传播环境使得上下行信道间仍存在一定相关性，不能简单地用两个独立的信道模型加以生成。如何在建模与仿真过程中，充分刻画这种相关性，并据此降低信道仿真过程的复杂度，目前仍存在较大难度。针非对称结构的对毫米波和亚毫米波频段的无线信道构建模型，形成流程化的建模方法，兼顾模型准确性和效率，对未来毫米波和亚毫米波移动通信系统的设计和部署具有十分重要的意义。
发明内容
本发明公开的一种非对称毫米波亚毫米波无线信道仿真方法，该方法解决了现有无线信道仿真技术尚未充分考虑传播环境非对称性的问题；解决了非对称毫米波亚毫米波上下行无线信道准确、高效的联合生成问题；解决了非对称的上下行信道之间不存在互易性，对称信道的仿真方法不能直接用于非对称信道的问题；能够应用于无线通信领域，为非对称毫米波亚毫米波通信系统设计、部署、优化提供理论与模型依据。
本发明的技术方案如下：
一种非对称毫米波亚毫米波无线信道仿真方法，包括以下步骤：
步骤一、建立散射体模型，并通过筛选建立非对称有效散射体模型；
非对称散射体模型，其散射体在球面上的分布可建模为服从von Mises-Fisher随机分布。von Mises-Fisher是一种应用于球上的方向随机分布，三维单位向量x分布在立体球面上，其von Mises-Fisher概率密度函数表达式为：
f3(x；μ,ζ)＝C3(ζ)exp(ζμTx) (1)
上式中，向量μ表示平均方向，满足||μ||＝1。ζ为聚集因子，且ζ≥0。C3(ζ)表示归一化常量，由下式给出：
Figure BDA0003103918810000021
在散射体球状分布的情况下，随机变量x由散射体的方位角θ和俯仰角φ构成：x＝[cosφcosθ,cosφsinθ,sinφ]。此外，向量μ＝[cosφ0cosθ0,cosφ0sinθ0,sinφ0]，这里的θ0，φ0分别为散射体的平均方位角和俯仰角。基站侧和移动台侧的散射体球状分布概率密度函数为：
Figure BDA0003103918810000022
有效散射体的生成过程包括以下步骤：
1.选择传输侧u(u为基站侧T或移动台侧R，即u∈{T,R})，根据散射体球状分布概率密度函数，该侧共计Nu个散射体中，生成的第i个(其中i∈{1,2,…,Nu})散射体记为
Figure BDA0003103918810000023
2.判断传输侧u的散射体
Figure BDA0003103918810000024
是否在链路方向
Figure BDA0003103918810000025
(
Figure BDA0003103918810000026
为上行链路UL或下行链路DL，即
Figure BDA0003103918810000027
的有效散射体内，即判断散射体
Figure BDA0003103918810000028
所对应的水平方位角
Figure BDA0003103918810000029
和垂直俯仰角
Figure BDA00031039188100000210
是否满足：
Figure BDA00031039188100000211
其中
Figure BDA00031039188100000212
和
Figure BDA00031039188100000213
分别表示波束中心方向的水平方位角和垂直俯仰角，
Figure BDA00031039188100000214
和
Figure BDA00031039188100000215
分别表示在水平和垂直方向上的波束宽度。
3.如果散射体
Figure BDA0003103918810000031
在有效范围内，则记入对应的链路方向
Figure BDA0003103918810000032
与传输侧u的有效散射体集合
Figure BDA0003103918810000033
否则不计入；
4.更改散射体i、链路方向
Figure BDA0003103918810000034
与传输侧u，重复步骤1至3，直到所有散射体都被遍历，并得到对应的有效散射体集合
Figure BDA0003103918810000035
步骤二、建立三维波束模型，将上行和下行的定向性波束分别建模为三维空间中的椎体形状；
三维波束模型为水平方向上角度宽度为α，竖直方向上角度宽度为β的椎体。波束的中心方向用水平方位角θB和垂直俯仰角φB描述，因此可以将该波束表示为B(θB,φB,α,β)，其天线方向图可记为F(θB,φB)。使用该三维波束模型，可以得到天线增益
Figure BDA0003103918810000036
步骤三、根据非对称有效散射体模型和三维波束模型，建立非对称毫米波亚毫米波信道的信道冲激响应。
非对称毫米波亚毫米波信道的信道冲激响应，可由矩阵
Figure BDA0003103918810000037
进行描述。其中
Figure BDA0003103918810000038
中的p(p＝1,2,…,P)表示基站侧共P个天线阵元中的第p个天线阵元，其位置向量表示为
Figure BDA0003103918810000039
q(q＝1,2,…,Q)表示移动台侧共Q个天线阵元中的第q个天线阵元，其位置向量表示为
Figure BDA00031039188100000310
所有的几何关系可以统一在以O为原点的xyz直角坐标系中表示，因此满足位置向量
Figure BDA00031039188100000311
其中OT和OR分别表示基站侧和移动台侧的散射体球心，移动台侧散射体分布的球体半径记为RT，移动台侧散射体分布的球半径记为RR，移动台与基站之间的水平距离记为D。
由于移动台的移动，角度和位置向量均具有时变性，因此散射体位置向量可以表示为：
Figure BDA00031039188100000312
其中散射体
Figure BDA00031039188100000313
表示基站侧的第i个散射体，该散射体相对于基站侧散射体球心OT的位置向量记为
Figure BDA00031039188100000314
散射体
Figure BDA00031039188100000315
表示移动台侧的第j个散射体，该散射体相对于移动台侧散射体球心OR的位置向量记为
Figure BDA0003103918810000041
VMS表示移动台的移动速度矢量，t为瞬时时刻，则视距径路径向量可以表示为：
Figure BDA0003103918810000042
因此，从移动台侧散射体球心OR到基站侧散射体
Figure BDA0003103918810000043
以及从基站侧散射体球心OT到移动台侧散射体
Figure BDA0003103918810000044
的单跳路径向量表示为：
Figure BDA0003103918810000045
从基站侧散射体
Figure BDA0003103918810000046
到移动台侧散射体
Figure BDA0003103918810000047
的双跳路径向量
Figure BDA0003103918810000048
表示为：
Figure BDA0003103918810000049
考虑到毫米波亚毫米波在非视距径条件下，每发生一次反射，都会产生数十dB的漫散射损耗，因此在模型中对三次及更高次的反射忽略不计。
Figure BDA00031039188100000410
由视距径，单跳，双跳分量叠加组成：
Figure BDA00031039188100000411
上式中K为莱斯因子，表示视距径的功率占总功率的比例。将所有的几何关系都统一在以O为原点的xyz直角坐标系中表示。SBu(u∈{T,R})表示两种情况下的单跳，u＝T时单跳的散射体位于基站侧，u＝R时散射体位于移动台侧。
Figure BDA00031039188100000412
分别表示
Figure BDA00031039188100000413
方向上的从基站侧的第p个天线阵元到移动台侧第q个天线阵元的视距径、单跳和双跳分量的信道冲激响应函数。
信道冲激响应函数的视距径、单跳和双跳分量中，均涉及到收发天线增益
Figure BDA00031039188100000414
以及由于自由空间扩散和空气分子吸收效应电磁波信号所经历的衰减A(d)，可分别由下式计算得到：
Figure BDA00031039188100000415
Figure BDA00031039188100000416
Figure BDA00031039188100000417
为基站侧第p个阵元的天线方向图，
Figure BDA00031039188100000418
为移动台侧第q个阵元的天线方向图。d表示传播距离，Kf表示空气吸收系数，主要由空气的成分组成决定，且随着信号频率f的变化而变化，c表示光速。
对于视距径分量，
Figure BDA0003103918810000051
可由以下步骤得出：
1、收发天线增益
Figure BDA0003103918810000052
可由(10)式计算得到，将d＝WLoS(t)代入(11)式可得到视距径下信号所经历的自由空间扩散和空气分子吸收带来的总衰减ALoS(t)。
2、计算视距径下移动台相对于基站的多普勒频移：
Figure BDA0003103918810000053
其中
Figure BDA0003103918810000054
为视距径向量WLoS(t)的单位向量；
3、计算视距径下移动台相对于基站的时延：
Figure BDA0003103918810000055
其中WLoS(t)为视距径向量WLoS(t)的模；
4、计算
Figure BDA0003103918810000056
Figure BDA0003103918810000057
其中
Figure BDA0003103918810000058
为波数，ΦLoS为视距径随机相位，在(0,2π]范围内均匀分布。
对于单跳路径分量，
Figure BDA0003103918810000059
可由以下步骤得出：
1、根据(10)式计算天线增益
Figure BDA00031039188100000510
将
Figure BDA00031039188100000511
代入(11)式可得到有效散射体
Figure BDA00031039188100000512
单跳所经历的自由空间扩散和空气分子吸收的总衰减
Figure BDA00031039188100000513
其中标量
Figure BDA00031039188100000514
和
Figure BDA00031039188100000515
分别为向量
Figure BDA00031039188100000516
和
Figure BDA00031039188100000517
所对应的模长；
2、计算单跳下移动台相对于基站的多普勒频移：
Figure BDA00031039188100000518
其中
Figure BDA00031039188100000519
为单跳下径向量
Figure BDA00031039188100000520
的单位向量；
3、计算单跳下移动台相对于基站的时延：
Figure BDA0003103918810000061
4、计算单跳漫散射功率损耗χSB；
单跳及双跳漫散射功率损耗χ，可采用如下步骤计算。其中χSB代表单跳漫散射功率损耗，χDB代表双跳漫散射功率损耗。
计算单跳漫散射功率损耗χSB。其中假设入射波以θin入射到平面上，以入射点为原点，散射平面为xoy平面，建立三维直角坐标系，z轴为法线方向。入射波所在的平面称为入射面，入射角记为θin，垂直方向出射角记为θout，水平方向出射角
Figure BDA0003103918810000062
假设这三个角均在一定角度范围内均匀分布，即：θin～U(0,π/2)，θout～U(0,π/2)和
Figure BDA0003103918810000063
考虑信道中存在N∈种散射体，第∈(∈＝1,…,N∈)种散射体占总散射体的概率为P∈(满足
Figure BDA0003103918810000064
)。其对应的粗糙表面均方根高度记为ρh,∈，折射率记为nt,∈。将这N∈种散射体的粗糙表面均方根高度和折射率的期望值，作为所有散射体粗糙表面的平均特征，即
Figure BDA0003103918810000065
根据θin，θout和
Figure BDA0003103918810000066
之间的关系，χSB分为反射和散射两种情况进行计算。
①对于反射的情况(θin＝θout且
Figure BDA0003103918810000067
)，使用下式计算。其中瑞利粗糙因子记为ρr，菲涅尔反射系数记为Fr。
Figure BDA0003103918810000068
其中，瑞利粗糙因子ρr和菲涅尔反射系数Fr使用下式计算。
Figure BDA0003103918810000069
Figure BDA00031039188100000610
②对于散射的情况(θin≠θout或
Figure BDA00031039188100000611
)，使用下式计算。其中散射损耗因子记为ρs，菲涅尔反射系数记为Fr。
Figure BDA00031039188100000612
5、计算
Figure BDA0003103918810000071
Figure BDA0003103918810000072
Figure BDA0003103918810000073
为
Figure BDA0003103918810000074
单跳的随机相位，在(0,2π]范围内均匀分布；
Figure BDA0003103918810000075
表示
Figure BDA0003103918810000076
方向上u侧的有效散射体总数；
Figure BDA0003103918810000077
分别为单跳路径向量
Figure BDA0003103918810000078
对应的单位向量，后两者分别对应于从u侧的p天线到
Figure BDA0003103918810000079
以及从
Figure BDA00031039188100000710
到u侧的q天线的路径向量。
对于双跳路径分量，
Figure BDA00031039188100000711
可由以下步骤得出：
1、天线增益
Figure BDA00031039188100000712
可由(10)式计算得到，将
Figure BDA00031039188100000713
代入(11)式可得到有效散射体
Figure BDA00031039188100000714
和
Figure BDA00031039188100000715
双跳对应的空气分子吸收效应导致的衰减
Figure BDA00031039188100000716
其中
Figure BDA00031039188100000717
为从
Figure BDA00031039188100000718
到
Figure BDA00031039188100000719
的径向量
Figure BDA00031039188100000720
的模长；
2、计算双跳下接收端的多普勒频移：
Figure BDA00031039188100000721
3、计算双跳下移动台相对于基站的时延；
Figure BDA00031039188100000722
4、双跳漫散射功率损耗χDB可近似为：
χDB＝(χSB)2 (24)
5、计算
Figure BDA00031039188100000723
Figure BDA00031039188100000724
Figure BDA00031039188100000725
为随机相位，在(0,2π]范围内均匀分布；
Figure BDA00031039188100000726
表示基站侧和移动台侧的有效散射体总数；
Figure BDA00031039188100000727
分别为
Figure BDA00031039188100000728
对应的单位向量，后两者分别对应于从u侧的p天线到
Figure BDA0003103918810000081
以及从
Figure BDA0003103918810000082
到u侧的q天线的路径向量。
有益效果：
1、本发明充分考虑了信道分布环境呈现出非对称性的特点，设计了一种毫米波亚毫米波信道仿真方法，该方法能够生成非对称条件下的毫米波亚毫米波无线信道冲激响应和参数。
2、本发明结合毫米波亚毫米波具有空气分子吸收和粗糙表面漫散射的传播特点，分为视距径、单跳和双跳对上下行信道中每条多径的参数进行刻画，从而能够更为准确地实现对实际信道的模拟。
3、本发明利用非对称条件下上下行信道的相关性，有效降低了生成信道参数和信道冲激响应的复杂度，从而增强了其在系统仿真验证应用中的实用性。
附图说明
图1为本发明一种非对称毫米波亚毫米波室外信道仿真方法的流程图。
图2为非对称毫米波亚毫米波室外信道模型几何关系图。
图3为天线波束和对应的三维椎体模型示意图。
图4为在默认波束宽度组合条件下非对称有效散射体的三维空间分布，其中(a)为基站侧，(b)为移动台侧。
图5为非对称毫米波亚毫米波室外信道模型仿真下行功率角度谱三维空间分布，其中(a)为基站侧的分布情况，(b)为移动台侧的分布情况。
图6为非对称毫米波亚毫米波室外信道模型仿真上行功率角度谱三维空间分布，其中(a)为基站侧的分布情况，(b)为移动台侧的分布情况。
图7为非对称毫米波亚毫米波室外信道模型仿真上下行信道功率时延谱分布，其中(a)为下行信道分布情况，(b)为上行信道分布情况。
具体实施方式
下面结合附图并举实施例，对本发明的一种非对称毫米波亚毫米波室外信道仿真方法进行详细描述。
本发明提供的一种非对称毫米波亚毫米波室外信道仿真方法，可应用于室外场景下非对称的毫米波亚毫米波频段无线信道的建模、设计和优化中，对非对称移动通信系统在室外的大规模快速部署提供信道模型支撑。
本发明提供的非对称毫米波亚毫米波室外信道仿真方法，其总的建模流程如图1所示。
首先，建立三维波束模型，将上行和下行的定向性波束分别建模为三维空间中的椎体形状。天线波束和对应的三维椎体模型示意图如图3所示。仿真选取的默认波束宽度组合如表1所示。
表1三维波束模型的默认波束宽度组合
Figure BDA0003103918810000091
本案例中在有效散射区域EA内的收发天线波束增益F(α,β)视为是相等的。因此对于水平角度宽度和垂直角度宽度分别为α,β的波束，其天线增益可由下式计算得出。
Figure BDA0003103918810000092
然后，建立散射体模型，并筛选建立非对称有效散射体模型。生成三维空间中散射体的von Mises-Fisher分布在基站侧的参数设置为：ζ＝4,θ0＝0,φ0＝-π/4，移动台侧的参数设置为：ζ＝2,θ0＝π,φ0＝π/8。基站侧散射体分布的球体半径RT＝40m，球上总的散射体数量设为NT＝100；移动台侧球半径为RR＝10m，散射体总数NR＝200。
在默认波束宽度组合条件下根据散射体球状分布概率密度函数f3(x；μ,ζ)，可以得到非对称有效散射体三维空间分布，如图4所示。经过统计，上下行有效散射体个数分别为：
Figure BDA0003103918810000093
可以看出，上下行信道波束宽度的不同，会导致上下行信道中有效散射体的数量存在显著差异，使得上下行信道不再具有互易性。
最后，根据三维波束模型和非对称有效散射体模型，建立非对称毫米波亚毫米波信道的各个参数以及信道冲激响应。非对称毫米波亚毫米波室外信道模型几何关系图如图2所示。信道模型仿真中各个场景参数的具体设置为，基站与移动台之间的水平初始距离D＝70m，基站的高度HT＝50m，移动台的高度HR＝1.5m。仿真频率为f＝350GHz，自由空间传播环境，分子吸收系数Kf＝0.00161m-1，莱斯因子取K＝30dB。假设信道中主要考虑建筑外墙，树木，水泥地面三种散射体类型，每种散射体的占比P∈分别设为{0.6,0.2,0.2}，均方根表面高度ρh分别设为{0.05,0.08,0.25}mm，折射率nt,∈分别为{2.5,1.3,2.1}。
基站侧的天线阵列配置为Px×Pz＝12×12，移动台侧为Qx×Qz＝3×3，阵元间距设为d＝λ/2。对于第(px,pz)个天线阵元，其位置向量为：
Figure BDA0003103918810000101
散射体
Figure BDA0003103918810000102
的位置向量为
Figure BDA0003103918810000103
其对应的单位向量为
Figure BDA0003103918810000104
从p天线到
Figure BDA0003103918810000105
的波出发方向视为与
Figure BDA0003103918810000106
相同，因此：
Figure BDA0003103918810000107
分别计算收发双方各天线阵元间非对称双向链路的信道参数，包括视距径、单跳径和双跳径的天线增益
Figure BDA0003103918810000108
衰减A(d)、多普勒频移fd、径时延τ以及漫散射功率损耗χ等。根据各径的参数结果计算出视距径、单跳和双跳的信道冲激响应
Figure BDA0003103918810000109
和
Figure BDA00031039188100001010
进一步根据(9)式得到收发双方各天线阵元间非对称双向链路的信道冲激响应
Figure BDA00031039188100001011
并由
Figure BDA00031039188100001012
计算得到整个非对称信道的信道冲激响应矩阵
Figure BDA00031039188100001013
由于上下行信道环境仍然存在一定的相关性，利用同一组信道参数支持同时生成上下行的信道参数和信道冲激响应函数，从而对信道做进一步的分析，可以提高建模效率，减小计算量。
通过同时生成上下行非对称信道的角度功率谱分布和功率时延谱，可以验证该信道仿真方法的可行性。功率角度谱的计算方式为：
Figure BDA00031039188100001014
其中
Figure BDA00031039188100001015
表示
Figure BDA00031039188100001016
方向的信道中径的总数，按照在基站侧和移动台侧发生双跳的散射体簇之间采用一一匹配的规则，
Figure BDA00031039188100001017
aLoS,aSB,aDB分别为信道冲激响应中视距径分量、单跳路径分量和多跳路径分量对应的幅度信息。θ,φ分别为每条径对应的水平方位角和垂直方位角。
功率时延谱的计算方式为：
Figure BDA00031039188100001018
其中，τLoS为LoS径的时延，
Figure BDA00031039188100001019
为单跳
Figure BDA00031039188100001020
对应的时延，τj为第j对双跳对应的时延。
如图4所示，由于上下行信道环境包含了重叠的有效散射体，因此两者间仍存在一定的关联性，因此可以通过同一组信道参数同时生成上下行信道的信道冲激响应，并得到上下行功率角度谱的三维空间分布。在采用默认波束宽度组合的条件下，得到了与之对应的下行(如图5所示)和上行(如图6所示)功率角度谱三维空间分布图，其中图5(a)和图6(a)为基站侧，图5(b)和图6(b)为移动台侧。此外，同一次生成的下行(如图7(a)所示)和上行(如图7(b)所示)的信道功率时延谱分布。从图5、图6和图7中可看出，在基站侧的单跳比移动台侧具有更大的空间和时延分布范围，且上行信道高于下行信道。随着接收端波束宽度变宽，时延扩展也越大。这与在实际非对称信道中，基站周围的簇距离基站通常较远，即使在定向性波束条件下，各条多径的时延仍具有一定差异是相吻合的。
在本案例中，利用非对称的上下行信道之间散射环境的关联性，可以在同一组信道环境参数的条件下生成上下行信道的各个参数及信道冲激响应。相比于独立地对上行和下行信道进行分析，本发明有效降低了信道仿真分析的复杂度，对实际的非对称毫米波亚毫米波频段的室外无线通信信道仿真和设计具有一定的实用性。
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种非对称毫米波亚毫米波无线信道仿真方法，其特征在于，包括如下步骤：
步骤一、建立散射体模型，并通过筛选建立非对称有效散射体模型；
非对称散射体模型，其散射体在球面上的分布可建模为服从von Mises-Fisher随机分布；von Mises-Fisher是一种应用于球上的方向随机分布，三维单位向量x分布在立体球面上，其von Mises-Fisher概率密度函数表达式为：
f3(x；μ,ζ)＝C3(ζ)exp(ζμTx) (1)
上式中，向量μ表示平均方向，满足||μ||＝1；ζ为聚集因子，且ζ≥0；C3(ζ)表示归一化常量，由下式给出：
Figure FDA0003818346140000011
在散射体球状分布的情况下，随机变量x由散射体的方位角θ和俯仰角φ构成：x＝[cosφcosθ,cosφsinθ,sinφ]；此外，向量μ＝[cosφ0cosθ0,cosφ0sinθ0,sinφ0]，这里的θ0，φ0分别为散射体的平均方位角和俯仰角；基站侧和移动台侧的散射体球状分布概率密度函数为：
Figure FDA0003818346140000012
有效散射体的生成过程包括以下步骤：
1)选择传输侧u(u为基站侧T或移动台侧R，即u∈{T,R})，根据散射体球状分布概率密度函数，该侧共计Nu个散射体中，生成的第i个(其中i∈{1,2,…,Nu})散射体记为
Figure FDA0003818346140000013
2)判断传输侧u的散射体
Figure FDA0003818346140000014
是否在链路方向
Figure FDA0003818346140000015
(
Figure FDA0003818346140000016
为上行链路UL或下行链路DL，即
Figure FDA0003818346140000017
)的有效散射体内，即判断散射体
Figure FDA0003818346140000018
所对应的水平方位角
Figure FDA0003818346140000019
和垂直俯仰角
Figure FDA00038183461400000110
是否满足：
Figure FDA00038183461400000111
其中
Figure FDA00038183461400000112
和
Figure FDA00038183461400000113
分别表示波束中心方向的水平方位角和垂直俯仰角，
Figure FDA00038183461400000114
和
Figure FDA00038183461400000115
分别表示在水平和垂直方向上的波束宽度；
3)如果散射体
Figure FDA00038183461400000116
在有效范围内，则记入对应的链路方向
Figure FDA00038183461400000117
与传输侧u的有效散射体集合
Figure FDA0003818346140000021
否则不计入；
4)更改散射体i、链路方向
Figure FDA0003818346140000022
与传输侧u，重复步骤1至3，直到所有散射体都被遍历，并得到对应的有效散射体集合
Figure FDA0003818346140000023
步骤二、建立三维波束模型，将上行和下行的定向性波束分别建模为三维空间中的椎体形状；
步骤三、根据非对称有效散射体模型和三维波束模型，建立非对称毫米波亚毫米波信道的信道冲激响应。
2.如权利要求1所述，一种非对称毫米波亚毫米波无线信道仿真方法，其特征在于，步骤二的实现方法为：
三维波束模型为水平方向上角度宽度为α，竖直方向上角度宽度为β的椎体；波束的中心方向用水平方位角θB和垂直俯仰角φB描述，因此可以将该波束表示为B(θB,φB,α,β)，其天线方向图可记为F(θB,φB)；使用该三维波束模型，可以得到天线增益
Figure FDA00038183461400000210
3.如权利要求1所述，一种非对称毫米波亚毫米波无线信道仿真方法，其特征在于，步骤三的实现方法为：
非对称毫米波亚毫米波信道的信道冲激响应，可由矩阵
Figure FDA0003818346140000024
进行描述；其中
Figure FDA0003818346140000025
中的p(p＝1,2,…,P)表示基站侧共P个天线阵元中的第p个天线阵元，其位置向量表示为
Figure FDA0003818346140000026
q(q＝1,2,…,Q)表示移动台侧共Q个天线阵元中的第q个天线阵元，其位置向量表示为
Figure FDA0003818346140000027
所有的几何关系可以统一在以O为原点的xyz直角坐标系中表示，因此满足位置向量
Figure FDA0003818346140000028
其中OT和OR分别表示基站侧和移动台侧的散射体球心，移动台侧散射体分布的球体半径记为RT，移动台侧散射体分布的球半径记为RR，移动台与基站之间的水平距离记为D；
由于移动台的移动，角度和位置向量均具有时变性，因此散射体位置向量可以表示为：
Figure FDA0003818346140000029
其中散射体
Figure FDA0003818346140000031
表示基站侧的第i个散射体，该散射体相对于基站侧散射体球心OT的位置向量记为
Figure FDA0003818346140000032
散射体
Figure FDA0003818346140000033
表示移动台侧的第j个散射体，该散射体相对于移动台侧散射体球心OR的位置向量记为
Figure FDA0003818346140000034
VMS表示移动台的移动速度矢量，t为瞬时时刻，则视距径路径向量可以表示为：
Figure FDA0003818346140000035
因此，从移动台侧散射体球心OR到基站侧散射体
Figure FDA0003818346140000036
以及从基站侧散射体球心OT到移动台侧散射体
Figure FDA0003818346140000037
的单跳路径向量表示为：
Figure FDA0003818346140000038
从基站侧散射体
Figure FDA0003818346140000039
到移动台侧散射体
Figure FDA00038183461400000310
的双跳路径向量
Figure FDA00038183461400000311
表示为：
Figure FDA00038183461400000312
考虑到毫米波亚毫米波在非视距径条件下，每发生一次反射，都会产生数十dB的漫散射损耗，因此在模型中对三次及更高次的反射忽略不计；
Figure FDA00038183461400000313
由视距径，单跳，双跳分量叠加组成：
Figure FDA00038183461400000314
上式中K为莱斯因子，表示视距径的功率占总功率的比例；将所有的几何关系都统一在以O为原点的xyz直角坐标系中表示；SBu(u∈{T,R})表示两种情况下的单跳，u＝T时单跳的散射体位于基站侧，u＝R时散射体位于移动台侧；
Figure FDA00038183461400000315
分别表示
Figure FDA00038183461400000316
方向上的从基站侧的第p个天线阵元到移动台侧第q个天线阵元的视距径、单跳和双跳分量的信道冲激响应函数；
信道冲激响应函数的视距径、单跳和双跳分量中，均涉及到收发天线增益
Figure FDA00038183461400000317
以及由于自由空间扩散和空气分子吸收效应电磁波信号所经历的衰减A(d)，可分别由下式计算得到：
Figure FDA00038183461400000318
Figure FDA0003818346140000041
Figure FDA0003818346140000042
为基站侧第p个阵元的天线方向图，
Figure FDA0003818346140000043
为移动台侧第q个阵元的天线方向图；d表示传播距离，Kf表示空气吸收系数，主要由空气的成分组成决定，且随着信号频率f的变化而变化，c表示光速；
对于视距径分量，
Figure FDA0003818346140000044
可由以下步骤得出：
1)收发天线增益
Figure FDA0003818346140000045
可由(10)式计算得到，将d＝WLoS(t)代入(11)式可得到视距径下信号所经历的自由空间扩散和空气分子吸收带来的总衰减ALoS(t)；
2)计算视距径下移动台相对于基站的多普勒频移：
Figure FDA0003818346140000046
其中
Figure FDA0003818346140000047
为视距径向量WLoS(t)的单位向量；
3)计算视距径下移动台相对于基站的时延：
Figure FDA0003818346140000048
其中WLoS(t)为视距径向量WLoS(t)的模；
4)计算
Figure FDA0003818346140000049
Figure FDA00038183461400000410
其中
Figure FDA00038183461400000411
为波数，ΦLoS为视距径随机相位，在(0,2π]范围内均匀分布；
对于单跳路径分量，
Figure FDA00038183461400000412
可由以下步骤得出：
1)根据(10)式计算天线增益
Figure FDA00038183461400000413
将
Figure FDA00038183461400000414
代入(11)式可得到有效散射体
Figure FDA00038183461400000415
单跳所经历的自由空间扩散和空气分子吸收的总衰减
Figure FDA00038183461400000416
其中标量
Figure FDA00038183461400000417
和
Figure FDA00038183461400000418
分别为向量
Figure FDA00038183461400000419
和
Figure FDA00038183461400000420
所对应的模长；
2)计算单跳下移动台相对于基站的多普勒频移：
Figure FDA00038183461400000421
其中
Figure FDA0003818346140000051
为单跳下径向量
Figure FDA0003818346140000052
的单位向量；
3)计算单跳下移动台相对于基站的时延：
Figure FDA0003818346140000053
4)计算单跳漫散射功率损耗χSB；
单跳及双跳漫散射功率损耗χ，可采用如下步骤计算；其中χSB代表单跳漫散射功率损耗，χDB代表双跳漫散射功率损耗；
计算单跳漫散射功率损耗χSB；其中假设入射波以θin入射到平面上，以入射点为原点，散射平面为xoy平面，建立三维直角坐标系，z轴为法线方向；入射波所在的平面称为入射面，入射角记为θin，垂直方向出射角记为θout，水平方向出射角
Figure FDA0003818346140000054
假设这三个角均在一定角度范围内均匀分布，即：θin～U(0,π/2)，θout～U(0,π/2)和
Figure FDA0003818346140000055
考虑信道中存在N∈种散射体，第∈(∈＝1,…,N∈)种散射体占总散射体的概率为P∈(满足
Figure FDA0003818346140000056
)；其对应的粗糙表面均方根高度记为ρh,∈，折射率记为nt,∈；将这N∈种散射体的粗糙表面均方根高度和折射率的期望值，作为所有散射体粗糙表面的平均特征，即
Figure FDA0003818346140000057
根据θin，θout和
Figure FDA00038183461400000512
之间的关系，χSB分为反射和散射两种情况进行计算；
①对于反射的情况(θin＝θout且
Figure FDA00038183461400000513
)，使用下式计算；其中瑞利粗糙因子记为ρr，菲涅尔反射系数记为Fr；
Figure FDA0003818346140000058
其中，瑞利粗糙因子ρr和菲涅尔反射系数Fr使用下式计算；
Figure FDA0003818346140000059
Figure FDA00038183461400000510
②对于散射的情况(θin≠θout或
Figure FDA00038183461400000511
)，使用下式计算；其中散射损耗因子记为ρs，菲涅尔反射系数记为Fr；
Figure FDA0003818346140000061
5)计算
Figure FDA0003818346140000062
Figure FDA0003818346140000063
Figure FDA0003818346140000064
为
Figure FDA0003818346140000065
单跳的随机相位，在(0,2π]范围内均匀分布；
Figure FDA0003818346140000066
表示
Figure FDA0003818346140000067
方向上u侧的有效散射体总数；
Figure FDA0003818346140000068
分别为单跳路径向量
Figure FDA0003818346140000069
对应的单位向量，后两者分别对应于从u侧的p天线到
Figure FDA00038183461400000610
以及从
Figure FDA00038183461400000611
到u侧的q天线的路径向量；
对于双跳路径分量，
Figure FDA00038183461400000612
可由以下步骤得出：
1)天线增益
Figure FDA00038183461400000613
可由(10)式计算得到，将
Figure FDA00038183461400000614
代入(11)式可得到有效散射体
Figure FDA00038183461400000615
和
Figure FDA00038183461400000616
双跳对应的空气分子吸收效应导致的衰减
Figure FDA00038183461400000617
其中
Figure FDA00038183461400000618
为从
Figure FDA00038183461400000619
到
Figure FDA00038183461400000620
的径向量
Figure FDA00038183461400000621
的模长；
2)计算双跳下接收端的多普勒频移：
Figure FDA00038183461400000622
3)计算双跳下移动台相对于基站的时延；
Figure FDA00038183461400000623
4)双跳漫散射功率损耗χDB可近似为：
χDB＝(χSB)2 (24)
5)计算
Figure FDA00038183461400000624
Figure FDA00038183461400000625
Figure FDA0003818346140000071
为随机相位，在(0,2π]范围内均匀分布；
Figure FDA0003818346140000072
表示基站侧和移动台侧的有效散射体总数；
Figure FDA0003818346140000073
分别为
Figure FDA0003818346140000074
对应的单位向量，后两者分别对应于从u侧的p天线到
Figure FDA0003818346140000075
以及从
Figure FDA0003818346140000076
到u侧的q天线的路径向量。
CN202110635458.4A 2021-06-07 2021-06-07 一种非对称毫米波亚毫米波无线信道仿真方法 Active CN113364544B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110635458.4A CN113364544B (zh) 2021-06-07 2021-06-07 一种非对称毫米波亚毫米波无线信道仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110635458.4A CN113364544B (zh) 2021-06-07 2021-06-07 一种非对称毫米波亚毫米波无线信道仿真方法

Publications (2)

Publication Number Publication Date
CN113364544A CN113364544A (zh) 2021-09-07
CN113364544B true CN113364544B (zh) 2022-11-04

Family

ID=77533106

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110635458.4A Active CN113364544B (zh) 2021-06-07 2021-06-07 一种非对称毫米波亚毫米波无线信道仿真方法

Country Status (1)

Country Link
CN (1) CN113364544B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113890583B (zh) * 2021-11-12 2024-06-14 北京理工大学 一种毫米波亚毫米波频段非对称信道参数预测方法
CN114172529B (zh) * 2021-12-13 2023-03-21 华南理工大学 一种适用于毫米波移动通信系统的收发非对称装置和方法
CN114640414B (zh) * 2022-02-15 2022-12-02 北京交通大学 一种基于角度联合概率的非对称波束信道多径演化方法
CN114726464B (zh) * 2022-03-07 2024-01-30 东南大学 一种上下行非对称信道模型参数生成方法
CN115225179B (zh) * 2022-07-14 2024-02-09 重庆邮电大学 一种应用面向高机动平台的短波宽带移动信道模型的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108365903B (zh) * 2018-01-29 2021-02-02 哈尔滨工程大学 一种基于随机散射簇的三维Massive MIMO信道建模方法
CN108736995A (zh) * 2018-06-11 2018-11-02 北京科技大学 一种毫米波无线信道建模方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
5G毫米波QuaDRiGa平台信道仿真与验证研究;李树等;《电波科学学报》;20170415(第02期);全文 *
基于OFDM的毫米波通信多径信道模型分析;陈婷等;《火控雷达技术》;20200625(第02期);全文 *

Also Published As

Publication number Publication date
CN113364544A (zh) 2021-09-07

Similar Documents

Publication Publication Date Title
CN113364544B (zh) 一种非对称毫米波亚毫米波无线信道仿真方法
Ju et al. Millimeter wave and sub-terahertz spatial statistical channel model for an indoor office building
KR102204783B1 (ko) 딥러닝 기반의 빔포밍 통신 시스템 및 방법
Hemadeh et al. Millimeter-wave communications: Physical channel models, design considerations, antenna constructions, and link-budget
Wu et al. Interference and coverage analysis for terahertz networks with indoor blockage effects and line-of-sight access point association
Zhang et al. 5G millimeter-wave antenna array: Design and challenges
Zhang et al. Measurement-based 5G millimeter-wave propagation characterization in vegetated suburban macrocell environments
Raghavan et al. Millimeter wave channel measurements and implications for PHY layer design
Xu et al. Measurements and models for 38-GHz point-to-multipoint radiowave propagation
Ertel et al. Overview of spatial channel models for antenna array communication systems
JP5139462B2 (ja) セルラー無線ネットワーク内におけるマルチアンテナの配備によって得られる容量および有効範囲の利得を近似し最適化する方法
Nguyen et al. Evaluation of empirical ray-tracing model for an urban outdoor scenario at 73 GHz E-band
Anjinappa et al. Base station and passive reflectors placement for urban mmWave networks
US6173191B1 (en) Localization of shaped directional transmitting and transmitting/receiving antenna array
Cheng et al. Modeling and simulation for UAV air-to-ground mmWave channels
Mukherjee et al. Analytical calculation of Rician K-factor for indoor wireless channel models
Nawaz et al. Effect of directional antenna on the Doppler spectrum in 3-D mobile radio propagation environment
Huang et al. A novel 3D GBSM for mmWave MIMO channels
Anjinappa et al. Millimeter-wave V2X channels: Propagation statistics, beamforming, and blockage
Khan et al. A generalized model for the spatial characteristics of the cellular mobile channel
Olenko et al. Analytically derived uplink/downlink TOA and 2-D-DOA distributions with scatterers in a 3-D hemispheroid surrounding the mobile
Hasan et al. Performance estimation of massive MIMO drop-based propagation channel model for mmWave communication
WO2023169596A1 (zh) 一种上下行非对称信道模型参数生成方法
Guey et al. Modeling and evaluation of beam tracking in mobile millimeter wave communication
Tsoulos et al. Adaptive antennas for microcellular and mixed cell environments with DS-CDMA

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