空地毫米波通信链路传播路径损耗计算方法与流程

技术领域：
：本发明涉及一种空地毫米波通信链路传播路径损耗计算方法，特别针对不同地面场景、高度、频段和距离的空地链路传播路径损耗的计算，属于无线信息传输领域。
背景技术：
：：随着毫米波应用越来越广泛，对毫米波通信系统的可靠性和有效性要求趋于严格。为实现频谱资源的高效利用，需要根据毫米波信道特性采取一系列抗衰落和干扰的措施。因此，加强对毫米波信道的研究使得我们能够更加深刻地理解其传播损耗特性。随着通信、微电子、材料等领域技术的兴起和发展，无人机(unmannedaerialvehicle,uav)技术发展迅速，得到全球各个国家的关注。基于uav平台的空地毫米波通信技术，得到了学术界和工业界的广泛重视。比如，uav毫米波空中基站被认为是下一代移动通信系统的典型应用。空地毫米波通信由于传播地理环境导致的反射和绕射等影响，地面接收机接收到的信号通常包含直射路径、地面反射分量及周围建筑物绕射分量。深入研究空地信道尤其是传播损耗特性反映了毫米波信号在特定通信场景中的大尺度衰落，是空地毫米波通信网络规划和系统设计的重要参考依据，对空中和地面节点的通信系统设计和优化具有重要意义。技术实现要素：：本发明提供了一种基于场景分类的空地毫米波通信链路传播路径损耗计算方法，该方法无需提供精确的传播场景数字地图，而是通过前期典型场景大数据分析结果获取场景、高度、频率和距离依赖的模型参数。本发明所采用的技术方案有：一种空地毫米波通信链路传播路径损耗计算方法，包括如下步骤：第一步，用户输入空地毫米波通信链路的参数，包括地面场景类型参数、空中和地面节点位置的坐标参数、以及通信频率参数，利用如下方法进行传播路径平均损耗计算lsum＝plos·llos+pref·lref+pdiff·ldiff+pns·lns(17)其中，plos、pref、pdiff和pns分别为该场景下传播过程出现视距路径、无视距路径但存在反射路径、无视距和反射路径但存在绕射路径和接收不到信号情况的统计概率，llos、lref、ldiff和lns分别为该场景下视距、反射、绕射和无传播路径情况下对应的信号损耗统计值；第二步，针对分类场景下空地毫米波通信链路出现视距、反射、绕射和无传播路径概率的具体方法如下：1)存在视距传播路径的概率计算方法如下，其中，d2d是地面节点和空中节点的二维距离，通过地面和空中节点的坐标参数计算得到，其中，ha是空中节点相对地面节点的高度，参数a1_los，a2_los，b1_los，b2_los，b3_los通过查视距传播路径概率的参数取值表获得；2)不存在视距传播路径，但是存在反射传播路径的概率计算方法如下，其中，其中，参数a1_ref，a2_ref，a3_ref，b1_ref，b2_ref，b3_ref，c1_ref，c2_ref，c3_ref通过查反射路径概率的参数取值表得到；3)不存在视距和反射传播路径，但是存在绕射传播路径的概率计算方法如下，其中，其中，参数a1_diff，a2_diff，a3_diff，b1_diff，b2_diff，b3_diff，c1_diff，c2_diff，c3_diff通过查绕射传播路径概率的参数取值表得到；4)无传播路径的概率计算方法如下，pns＝1-plos-pref-pdiff(29)其中，plos、pref和pdiff可通过前三步得到；第三步，针对视距情况下的空地毫米波通信路径损耗的计算方法如下，其中，fc为通信频率，d3d是地面节点和空中节点的三维距离，为零均值高斯随机变量，标准差单位为db；第四步，针对反射和绕射情况下的空地毫米波通信路径损耗计算方法如下，其中，n为传播损耗指数，为零均值高斯随机变量，单位为db。进一步地，第三步中：产生方法如下，首先生成一个均值为0，方差为1的标准正态分布χ，令其中，σlos的值可通过查视距情况σlos参数的取值表得到。进一步地，第四步中：1)传播损耗指数具体计算方法如下n＝j·hak(32)其中，j和k是与环境有关的参数，本发明不同场景的j和k取值如传播损耗模型参数的取值表所示；2)的产生方法如下，首先生成一个均值为0，方差为1的标准正态分布χ，令其中，σref/diff的取值可通过查反射和绕射情况σref/diff参数的取值表得到。本发明具有如下有益效果：(1).本发明提出的毫米波空地传播路径损耗计算方法，无需用户提供详细的传播场景数据地图，计算方法简单，计算结果的统计偏差较小；(2).本发明提出的毫米波空地传播路径损耗计算方法，考虑了空中节点的高度参数，比传统方法更适用于空地传播场景。附图说明：图1为本发明空地毫米波通信传播场景示意图。具体实施方式：下面结合附图对本发明作进一步的说明。本发明空地毫米波通信链路传播路径损耗计算方法，包括如下步骤：第一步，用户输入空地毫米波通信链路的参数，包括地面场景类型参数、空中和地面节点位置的坐标参数、以及通信频率参数，本发明利用如下方法进行传播路径平均损耗计算lsum＝plos·llos+pref·lref+pdiff·ldiff+pns·lns(33)其中，plos、pref、pdiff和pns分别为该场景下传播过程出现视距路径(简称视距情况)、无视距路径但存在反射路径(简称反射情况)、无视距和反射路径但存在绕射路径(简称绕射情况)和接收不到信号情况(简称无传播路径情况)的统计概率，llos、lref、ldiff和lns分别为该场景下视距、反射、绕射和无传播路径情况下对应的信号损耗统计值，其中无传播路径情况下损耗值默认设置为有效信号功率分贝值减去接收机灵敏度；第二步，本发明针对分类场景下空地毫米波通信链路出现视距、反射、绕射和无传播路径概率的具体方法如下：1)存在视距传播路径的概率计算方法如下，其中，d2d是地面节点和空中节点的二维距离，可以通过地面和空中节点的坐标参数计算得到，其中，ha是空中节点相对地面节点的高度，参数a1_los，a2_los，b1_los，b2_los，b3_los可通过查表1获得；表1视距传播路径概率的参数取值2)不存在视距传播路径，但是存在反射传播路径的概率计算方法如下，其中，其中，参数a1_ref，a2_ref，a3_ref，b1_ref，b2_ref，b3_ref，c1_ref，c2_ref，c3_ref可通过查表2得到；表2反射路径概率的参数取值场景郊区市区密集城市高层城市a1_ref0.26510.14920.16950.2079a2_ref8.7315.9440.6483.66a3_ref17.8329.431.6668.5b1_ref770.6809.8742.7773.2b2_ref37.149.0680.49174.2b3_ref41.8644.2187.29201.8c1_ref346296.2304.5316.8c2_ref22.7329.3638.2794.93c3_ref24.3626.1251.5694.263)不存在视距和反射传播路径，但是存在绕射传播路径的概率计算方法如下，其中，其中，参数a1_diff，a2_diff，a3_diff，b1_diff，b2_diff，b3_diff，c1_diff，c2_diff，c3_diff可通过查表3得到；表3绕射传播路径概率的参数取值场景郊区市区密集城市高层城市a1_diff0.7050.5680.33740.4009a2_diff7.4719.3848.8688.88a3_diff20.6728.1629.76104.2b1_diff882.2888.9781.9756.7b2_diff42.3664.8180.71214.1b3_diff90.1787.8380.11336.2c1_diff467.4353.7335.8325.2c2_diff25.013357.3160.12c3_diff64.324042.05135.24)无传播路径的概率计算方法如下，pns＝1-plos-pref-pdiff(45)其中，plos、pref和pdiff可通过前三步得到；第三步，本发明针对视距情况下的空地毫米波通信路径损耗的计算方法如下，其中，fc为通信频率，d3d是地面节点和空中节点的三维距离，为零均值高斯随机变量，标准差单位为db；进一步地，产生方法如下，首先生成一个均值为0，方差为1的标准正态分布χ，令其中，σlos的值可通过查表4得到；表4视距情况σlos参数的取值场景σlos郊区2.24市区1.44密集城市1.91高层城市2.18第四步，本发明针对反射和绕射情况下的空地毫米波通信路径损耗计算方法如下，其中，n为传播损耗指数，为零均值高斯随机变量，单位为db；1)传播损耗指数具体计算方法如下n＝j·hak(48)其中，j和k是与环境有关的参数，本发明不同场景的j和k取值如表5所示；表5传播损耗模型参数的取值2)的产生方法如下，首先生成一个均值为0，方差为1的标准正态分布χ，令其中，σref/diff的取值可通过查表6得到；表6反射和绕射情况σref/diff参数的取值下面以不同场景固定位置的空地节点为例，对技术方案进行清楚、完整的描述。第一步，用户输入空地毫米波通信链路的参数，包括地面场景类型，空中和地面节点位置的坐标参数，以及通信频率参数；本案例假设地面节点位置lg＝[0,0,2]，空中节点位置la＝[300,300,32]，通信频率28ghz，接收机灵敏度减去信号功率为255db；第二步，计算不同场景空地毫米波通信链路出现视距、反射和绕射情况的概率如下：1)结合地面和空中节点坐标，通信频率等参数，再通过查表1得到不同场景的a1_los，a2_los，b1_los，b2_los，b3_los，进而利用公式(34)计算出现视距情况的概率，在郊区场景、市区场景、密集城市场景和高层城市场景下分别为0.917、0.807、0.563、0.330；2)结合地面和空中节点坐标，通信频率等参数，再通过查表2得到不同场景的a1_ref，a2_ref，a3_ref，b1_ref，b2_ref，b3_ref，c1_ref，c2_ref，c3_ref，进而利用公式(37)计算出现反射情况的概率，在郊区场景、市区场景、密集城市场景和高层城市场景下分别为0.022、0.059、0.133、0.108；3)结合地面和空中节点坐标，通信频率等参数，再通过查表3得到不同场景的a1_diff，a2_diff，a3_diff，b1_diff，b2_diff，b3_diff，c1_diff，c2_diff，c3_diff，进而利用公式(41)计算出现绕射情况的概率，在郊区场景、市区场景、密集城市场景和高层城市场景下分别为0.061、0.134、0.184、0.240；4)结合plos、pref、pdiff和公式(45)计算可得到无传播路径的概率，在郊区场景、市区场景、密集城市场景和高层城市场景下分别为0.000、0.000、0.120、0.322；第三步，结合地面和空中节点坐标，通信频率等参数，通过查表4得到σlos，进而利用公式(46)计算视距情况下的空地毫米波通信路径损耗，在郊区场景、市区场景、密集城市场景和高层城市场景下分别为113.918、113.918、113.918、113.918db；第四步，计算反射和绕射情况下的空地毫米波通信路径损耗，具体包括：1)结合地面和空中节点的坐标、公式(48)和表5，计算获得不同场景下反射和绕射的传播损耗指数，反射的传播损耗指数在郊区场景、市区场景、密集城市场景和高层城市场景下分别为2.313、2.367、2.361、2.383，绕射的传播损耗指数在郊区场景、市区场景、密集城市场景和高层城市场景下分别为3.106、3.909、4.362、4.210；2)本案例中在郊区场景、市区场景、密集城市场景和高层城市场景下分别为2.13、2.92、2.77、3.01db，在郊区场景、市区场景、密集城市场景和高层城市场景下分别为4.55、4.37、5.98、5.25db；3)将地面和空中节点的坐标，传播损耗指数和零均值高斯随机变量代入公式(47)，可得到不同场景反射和绕射情况下的路径损耗，反射情况下的路径损耗在郊区场景、市区场景、密集城市场景和高层城市场景下分别为124.284、126.485、126.166、126.988db，绕射情况下的路径损耗在郊区场景、市区场景、密集城市场景和高层城市场景下分别为147.550、168.472、181.996、177.270db；第五步，将前述步骤获得的路径概率和平均损耗值代入公式(33)，最终获得本案例在郊区场景、市区场景、密集城市场景和高层城市场景下的毫米波传播路径损耗值分别为116.203、121.959、145.016、176.022db。以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明原理的前提下还可以作出若干改进，这些改进也应视为本发明的保护范围。当前第1页12