CN113395101B - 空天地融合网络中最小化传输时延的用户接入方法及装置 - Google Patents

空天地融合网络中最小化传输时延的用户接入方法及装置 Download PDF

Info

Publication number
CN113395101B
CN113395101B CN202110615508.2A CN202110615508A CN113395101B CN 113395101 B CN113395101 B CN 113395101B CN 202110615508 A CN202110615508 A CN 202110615508A CN 113395101 B CN113395101 B CN 113395101B
Authority
CN
China
Prior art keywords
user
satellite
time slot
signal
representing
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
CN202110615508.2A
Other languages
English (en)
Other versions
CN113395101A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202110615508.2A priority Critical patent/CN113395101B/zh
Publication of CN113395101A publication Critical patent/CN113395101A/zh
Application granted granted Critical
Publication of CN113395101B publication Critical patent/CN113395101B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/18502—Airborne stations
    • H04B7/18506—Communications with or from aircraft, i.e. aeronautical mobile service
    • H04B7/18508—Communications with or from aircraft, i.e. aeronautical mobile service with satellite system used as relay, i.e. aeronautical mobile satellite service
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/12—Arrangements for detecting or preventing errors in the information received by using return channel
    • H04L1/16—Arrangements for detecting or preventing errors in the information received by using return channel in which the return channel carries supervisory signals, e.g. repetition request signals
    • H04L1/18—Automatic repetition systems, e.g. Van Duuren systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W28/00—Network traffic management; Network resource management
    • H04W28/02—Traffic management, e.g. flow control or congestion control
    • H04W28/08—Load balancing or load distribution
    • H04W28/09—Management thereof
    • H04W28/0958—Management thereof based on metrics or performance parameters
    • H04W28/0967—Quality of Service [QoS] parameters
    • H04W28/0975—Quality of Service [QoS] parameters for reducing delays
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W74/00—Wireless channel access, e.g. scheduled or random access
    • H04W74/02—Hybrid access techniques
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W76/00—Connection management
    • H04W76/10—Connection setup

Abstract

一种空天地融合网络中最小化传输时延的用户接入方法和装置，涉及空天地融合网络技术领域。目前，采用空天地融合网络辅助边远地区的数据传输过程汇总，因为链路状态始终在发生动态变化，导致通信链路的可靠性变差。本发明所述的方法通过分别得到飞机和卫星接收端的信号接收形式并计算飞机端的接收信噪比，再得到地面用户‑飞机和卫星通信链路的中断概率，然后计算用户‑飞机和卫星的数据包平均发送次数，得到数据传输时延，在用户面对多种空中平台时，规划以最小化传输时延为目标的用户接入方法，将优化问题分解为一些独立的子问题，并根据本方法提出的启发式用户接入算法得到用户接入结果。该方法可以广泛应用在空天地融合网络领域中。

Description

空天地融合网络中最小化传输时延的用户接入方法及装置
技术领域
本发明涉及空天地融合网络技术领域，特别涉及一种空天地融合网络中最小化传输时延的用户接入方法。
背景技术
当前移动通信已经发展到第五代，通信系统的有效性和可靠性已经发展到前所未有的新阶段。出于对商业利益的考虑，移动网络运营商倾向于在人口密集的区建立传统蜂窝基站。当前，地面宽带无线信号的覆盖率仅有10％，边远地区的数据传输仍然较为困难，距离实现通信无缝覆盖的目标仍有很远的距离。
空天地融合网络由于引入卫星、飞机等空中平台进行辅助通信，极大的扩展了传统地面蜂窝网络的通信覆盖范围，使得边远地区的通信成为可能。与蜂窝移动通信不同，卫星、飞机等空中平台始终处于高速移动状态；因此，用户-平台的链路状态始终在发生动态变化，导致通信链路的可靠性变差，数据无法成功传输的概率也相应增加。目前，无论是针对飞机通信的研究，还是在卫星通信场景下考虑地面用户-空中平台通信链路中断情况的研究都十分有限。此外，在考虑数据重传的情况下，用户如何选择空中平台才能实现传输时延最小化的目标也成为了亟待解决的问题之一。
发明内容
为了解决现有技术中存在的由于空天地融合网络中地面用户-空中平台通信链路可能会发生中断的问题，本申请提供了空天地融合网络中最小化传输时延的用户接入方法及装置。
一种空天地融合网络中最小化传输时延的用户接入方法，包括如下步骤：
步骤一：根据t时刻空天地融合网络的网络拓扑结构、发射功率、衰落损耗、收发端天线增益等参数，得到飞机接收端的信号接收形式并计算飞机端的接收信噪比；同时得到卫星接收端的信号接收形式并计算卫星端的接收信噪比；
步骤二：根据信道条件、接收端信噪比门限值，利用莱斯信道系数的概率密度函数计算得到地面用户-飞机通信链路的中断概率和地面用户-卫星通信链路的中断概率；其中，地面用户-飞机通信链路简称G2A，地面用户-卫星通信链路简称G2S；
步骤三：引入ARQ机制，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，得到数据传输时延；同时根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，得到数据传输时延；
步骤四：用户面对多种空中平台时，规划以最小化传输时延为目标的用户接入方法，将优化问题分解为一些独立的子问题，并根据启发式用户接入算法得到用户接入结果。
上述步骤一中所述飞机和卫星端的接收信噪比为：
若时隙t地面用户m在飞机c的覆盖范围内，且用户通过G2A链路向飞机发射信号，则飞机端的信噪比的计算公式为：
Figure GDA0003685218520000021
其中，SNRm,c(t)表示时隙t内用户m向飞机c发射信号时、飞机c处的信噪比，
Figure GDA0003685218520000022
表示用户m向飞机c的发射功率，
Figure GDA0003685218520000023
表示用户m的发射天线增益，
Figure GDA0003685218520000024
表示飞机c的接收天线增益，
Figure GDA0003685218520000025
表示时隙t内用户m到飞机c的自由传输损耗，其中λc表示向飞机c发射的信号波长，dm,c(t)表示时隙t内用户m到飞机c的距离，La表示附加的损耗因子，hm,c(t)表示时隙t内用户m向飞机c传输时的快衰落系数，N0表示噪声功率；
若时隙t地面用户m通过G2S链路向卫星发射信号，则卫星端的信噪比的计算公式为：
Figure GDA0003685218520000026
其中，SNRm,s(t)表示时隙t内用户m向卫星s发射信号时、卫星s处的信噪比，
Figure GDA0003685218520000027
表示用户m向卫星s的发射功率
Figure GDA0003685218520000028
表示卫星s的接收天线增益，
Figure GDA0003685218520000029
表示用户m到卫星s的自由传输损耗，其中λs表示向卫星s的发射的信号波长，dm,s(t)表示时隙t内用户m到卫星s的距离，
Figure GDA00036852185200000210
表示时隙t内用户m向飞机c传输时的快衰落系数。
步骤二中所述的G2A的中断概率的计算方法为：当所述接收端信噪比大于设定门限值γth时，认为数据在该条通信链路能够成功传输；因此，G2A链路中断概率
Figure GDA00036852185200000211
为：
Figure GDA00036852185200000212
其中，
Figure GDA00036852185200000213
表示时隙t内用户m向飞机c发射信号时链路的中断概率，Pr表示概率，γth表示设定的信噪比门限值。λc表示向飞机a的发射的信号波长。
由于|hm,c(t)|2服从莱斯分布，因此，|hm,c(t)|2的概率密度函数为：
Figure GDA0003685218520000031
其中，f|h|(h)表示|h|的概率密度函数，h表示概率密度函数的自变量，|h|表示快衰落因子h的模，K表示莱斯分布中的莱斯因子，I0(·)表示第一类0阶修正贝塞尔函数,e为自然常数；
G2A链路中断概率为：
Figure GDA0003685218520000032
其中，dh表示对自变量h的微分，Kc表示G2A链路中快衰落因子的模服从的莱斯分布因子，
Figure GDA0003685218520000033
为马库姆Q函数，其中a表示常数，b表示马库姆Q函数设定的积分下限。
步骤二中所述的G2S的中断概率的计算方法为：
Figure GDA0003685218520000034
其中，Ks表示G2S链路中快衰落因子的模服从的莱斯分布因子；式中其它各参数的定义与前文相同。
步骤三中在G2A链路中利用ARQ机制进行数据传输时，平均传输次数的计算方式为：
Figure GDA0003685218520000035
其中，
Figure GDA0003685218520000036
表示时隙t内用户m向飞机c发送数据的平均传输次数，Pr表示概率，N表示设置的ARQ机制传输次数上限，n表示某一次的传输序号，
Figure GDA0003685218520000037
表示时隙t内第n次传输成功，
Figure GDA0003685218520000038
表示时隙t内第n次传输失败。
所以，用户利用G2A链路进行数据传输的平均传输次数为：
Figure GDA0003685218520000041
时隙t内G2A链路的数据传输时延为：
Figure GDA0003685218520000042
其中，D1,m,c(t)为时隙t内用户m向飞机c发送数据的传输时延，
Figure GDA0003685218520000043
为时隙t用户m发送的数据包个数，R为每个数据包所包含的比特数，B为通信带宽；
时隙t内G2S链路的数据传输时延的计算方法为：
Figure GDA0003685218520000044
其中，D2,m,s(t)为时隙t内用户m向卫星s发送数据的传输时延，
Figure GDA0003685218520000045
表示时隙t内用户m向卫星s发送数据的平均传输次数；式中其它各参数的定义与前文相同。
步骤四中所述优化问题规划如下：
P1-1为最小化系统内所有用户总传输时延的优化问题；由于每个时隙的优化方法相对独立，将P1-1分解为T个独立的子问题，令第t个时刻的子问题为P1-2。
P1-1:
Figure GDA0003685218520000046
上式为P1-1的目标函数，表示最小化整个研究时间段内所有用户发送数据的传输时延。其中，a1,m,c(t)表示在时隙t内用户m是否接入飞机m的判决变量，a2,m,s(t)表示在时隙t内用户m是否接入卫星s的判决变量；
Figure GDA0003685218520000047
上式为P1-1的限制条件1，表示每个时刻每个用户只能与通信范围内的一架飞机或一个卫星相连。其中，ωm,c(t)表示时隙t内用户m是否在飞机c覆盖范围内的二进制判决变量，M表示系统内用户总数量，m表示用户序号；
Figure GDA0003685218520000048
上式为P1-1的限制条件2，限制了每个时隙内每架飞机可接入的用户数量上限，其中，
Figure GDA0003685218520000049
表示每架飞机每个时刻所能接入的最大用户数量，C表示系统内飞机总数量，c表示飞机序号；
Figure GDA0003685218520000051
上式为P1-1的限制条件3，限制了每个时隙内每颗卫星可接入的用户数量上限，其中，
Figure GDA0003685218520000052
表示每颗卫星每个时刻所能接入的最大用户数量，S表示为系统内卫星总数量，s表示某颗卫星序号；
α1,m,c(t)∈{0,1},1≤t≤T,1≤m≤M,1≤c≤C
上式为P1-1的限制条件4，表示在时隙t内用户m是否接入飞机m的判决变量α1,m,c(t)是一个二进制变量，其中，T表示整个时间段内的总时隙个数，t表示某一时隙序号；
α2,m,s(t)∈{0,1},1≤t≤T,1≤m≤M,1≤s≤S
上式为P1-1的限制条件5，表示在时隙t内用户m是否接入卫星s的判决变量α2,m,s(t)是一个二进制变量。
步骤四中所述启发式用户接入算法为：
首先对优化变量进行凸松弛，即：0≤α1,m,c(t)≤1和0≤α2,m,s(t)≤1。讨论时隙t内用户m在飞机c的通信范围内的情况，即ωm,c(t)＝1；利用拉格朗日对偶函数，将P1-2转化为：
Figure GDA0003685218520000053
其中，μ1,c(t)表示P1-1的限制条件2所对应的拉格朗日乘子，μ2,s(t)表示P1-1的限制条件3所对应的拉格朗日乘子，
Figure GDA0003685218520000054
和
Figure GDA0003685218520000055
为相应的拉格朗日乘子最优解。
一种空天地融合网络中最小化传输时延的用户接入装置，包括：
信噪比计算装置，所述信噪比计算装置用于根据t时刻空天地融合网络的网络拓扑结构、发射功率、衰落损耗、收发端天线增益等参数，得到飞机接收端的信号接收形式并计算飞机端的接收信噪比；同时该装置还用于得到卫星接收端的信号接收形式并计算卫星端的接收信噪比；
中断概率计算装置，所述中断概率计算装置用于根据信道条件、接收端信噪比门限值，利用莱斯信道系数的概率密度函数计算得到地面用户-飞机通信链路的中断概率和地面用户-卫星通信链路的中断概率；其中，地面用户-飞机通信链路简称G2A，地面用户-卫星通信链路简称G2S；
数据传输时延装置，所述数据传输时延装置用于在引入ARQ机制后，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，得到数据传输时延；同时该装置还用于根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，得到数据传输时延；
结果确定装置，所述结果确定装置用于在用户面对多种空中平台时，规划以最小化传输时延为目标的用户接入方法，将优化问题分解为一些独立的子问题，并根据启发式用户接入算法得到用户接入结果。
所述信噪比计算装置还包括：飞机端信噪比计算模块和卫星端信噪比计算模块，其中，飞机端信噪比计算模块用于得到飞机接收端的信号接收形式并计算飞机端的接收信噪比，具体采用下述公式计算获得：
若时隙t地面用户m在飞机c的覆盖范围内，且用户通过G2A链路向飞机发射信号，则飞机端的接收信噪比的计算公式为：
Figure GDA0003685218520000061
其中，SNRm,c(t)表示时隙t内用户m向飞机c发射信号时飞机c处的信噪比，
Figure GDA0003685218520000062
表示用户m向飞机c的发射功率，
Figure GDA0003685218520000063
表示用户m的发射天线增益，
Figure GDA0003685218520000064
表示飞机c的接收天线增益，
Figure GDA0003685218520000065
表示用户m到飞机c的自由传输损耗，其中λc表示向飞机c发射的信号波长，dm,c(t)表示时隙t内用户m到飞机c的距离。La表示附加的损耗因子，hm,c(t)表示用户m向飞机c传输时的快衰落系数，N0表示噪声功率；
所述的卫星端信噪比模块用于得到卫星端的信号接收形式并计算卫星端的接收信噪比，具体采用下述公式计算获得：
若时隙t地面用户m通过G2S链路向卫星发射信号，则卫星端的接收信噪比的计算公式为：
Figure GDA0003685218520000066
其中，SNRm,s(t)表示时隙t内用户m向卫星s发射信号时卫星s处的信噪比，
Figure GDA0003685218520000067
表示卫星s的接收天线增益，
Figure GDA0003685218520000071
表示用户m到卫星s的自由传输损耗，其中λs表示向卫星s的发射的信号波长，dm,s(t)表示时隙t内用户m到卫星s的距离。式中其它各参数的定义与前文相同。
所述中断概率计算装置还包括：G2A的中断概率的计算模块和G2S的中断概率的计算模块，其中G2A的中断概率的计算模块用于计算得到地面用户-飞机通信链路的中断概率，其中地面用户-飞机通信链路简称G2A，所述的G2A的中断概率的计算模块计算G2A中断概率具体采用下述公式计算获得：
当所述接收端信噪比大于设定门限值γth时，认为数据在该条通信链路能够成功传输；因此，G2A链路中断概率
Figure GDA0003685218520000072
为：
Figure GDA0003685218520000073
其中，
Figure GDA0003685218520000074
表示时隙t内用户m向飞机c发射信号时链路的中断概率，Pr表示概率，γth表示设定的信噪比门限值。
由于|hm,c(t)|2服从莱斯分布，因此，|hm,c(t)|2的概率密度函数为：
Figure GDA0003685218520000075
其中，f|h|(h)表示|h|的概率密度函数，h表示概率密度函数的自变量，|h|表示快衰落因子h的模，K表示莱斯分布中的莱斯因子，I0(·)表示第一类0阶修正贝塞尔函数，e为自然常数。
G2A链路中断概率为：
Figure GDA0003685218520000076
其中，dh表示对自变量h的微分，Kc表示G2A链路中快衰落因子的模服从的莱斯分布因子，
Figure GDA0003685218520000077
为马库姆Q函数，其中a表示常数，b表示马库姆Q函数设定的积分下限。
所述的G2S的中断概率的计算模块用于计算得到地面用户-卫星通信链路的中断概率，其中，地面用户-卫星通信链路简称G2S，所述的G2S的中断概率的计算模块计算G2S中断概率具体采用下述公式计算获得：
Figure GDA0003685218520000081
其中，Ks表示G2S链路中快衰落因子的模服从的莱斯分布因子。式中其它各参数的定义与前文相同。
所述数据传输时延装置还包括：G2A平均传输次数模块和G2S平均传输次数模块，其中G2A平均传输次数模块用于在引入ARQ机制后，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，得到数据传输时延，所述的G2A平均传输次数模块计算传输次数具体采用下述公式计算获得：
Figure GDA0003685218520000082
其中，
Figure GDA0003685218520000083
表示时隙t内用户m向飞机c发送数据的平均传输次数，Pr表示概率，N表示设置的ARQ机制传输次数上限，n表示某一次的传输序号，
Figure GDA0003685218520000084
表示时隙t内第n次传输成功，
Figure GDA0003685218520000085
表示时隙t内第n次传输失败。
所以，用户利用G2A链路进行数据传输的平均传输次数为：
Figure GDA0003685218520000086
时隙t内G2A链路的数据传输时延为：
Figure GDA0003685218520000087
其中，D1,m,c(t)为时隙t内用户m向飞机c发送数据的传输时延，
Figure GDA0003685218520000088
为时隙t用户m发送的数据包个数，R为每个数据包所包含的比特数，B为通信带宽；
所述的G2S平均传输次数模块用于在引入ARQ机制后，根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，得到数据传输时延，所述的G2S平均传输次数模块计算时隙t内G2S链路的数据传输时延具体采用下述公式计算获得：
Figure GDA0003685218520000091
其中，D2,m,s(t)为时隙t内用户m向卫星s发送数据的传输时延，
Figure GDA0003685218520000092
表示时隙t内用户m向卫星s发送数据的平均传输次数；
所述的结果确定装置还包括：优化问题模块和得到结果模块，其中优化问题模块用于将优化问题分解为一些独立的子问题，所述的优化问题模块优化问题的规划如下：
P1-1为最小化系统内所有用户总传输时延的优化问题。由于每个时隙的优化方法相对独立，将P1-1分解为T个独立的子问题，令第t个时刻的子问题为P1-2。
P1-1:
Figure GDA0003685218520000093
上式为P1-1的目标函数，表示最小化整个研究时间段内所有用户发送数据的传输时延。其中，a1,m,c(t)表示在时隙t内用户m是否接入飞机m的判决变量，a2,m,s(t)表示在时隙t内用户m是否接入卫星s的判决变量；
Figure GDA0003685218520000094
上式为P1-1的限制条件1，表示每个时刻每个用户只能与通信范围内的一架飞机或一个卫星相连。其中，ωm,c(t)表示时隙t内用户m是否在飞机c覆盖范围内的二进制判决变量，M表示系统内用户总数量，m表示某个用户序号；
Figure GDA0003685218520000095
上式为P1-1的限制条件2，限制了每个时隙内每架飞机可接入的用户数量上限，其中，
Figure GDA0003685218520000096
表示每架飞机每个时刻所能接入的最大用户数量，C表示系统内飞机总数量，c表示某架飞机序号；
Figure GDA0003685218520000097
上式为P1-1的限制条件3，限制了每个时隙内每颗卫星可接入的用户数量上限，其中，
Figure GDA0003685218520000098
表示每颗卫星每个时刻所能接入的最大用户数量，S表示为系统内卫星总数量，s表示某颗卫星序号；
α1,m,c(t)∈{0,1},1≤t≤T,1≤m≤M,1≤c≤C
上式为P1-1的限制条件4，表示在时隙t内用户m是否接入飞机m的判决变量α1,m,c(t)是一个二进制变量，其中，T表示整个时间段内的总时隙个数，t表示某一时隙序号；
α2,m,s(t)∈{0,1},1≤t≤T,1≤m≤M,1≤s≤S
上式为P1-1的限制条件5，表示在时隙t内用户m是否接入卫星s的判决变量α2,m,s(t)是一个二进制变量。
所述的得到结果模块用于根据启发式用户接入算法得到用户接入结果，所述的得到结果模块中的启发式用户接入算法为：
首先对优化变量进行凸松弛，即：0≤α1,m,c(t)≤1和0≤α2,m,s(t)≤1。讨论时隙t内用户m在飞机c的通信范围内的情况，即ωm,c(t)＝1。利用拉格朗日对偶函数，将P1-2转化为：
Figure GDA0003685218520000101
其中，μ1,c(t)表示P1-1的限制条件2所对应的拉格朗日乘子，μ2,s(t)表示P1-1的限制条件3所对应的拉格朗日乘子，
Figure GDA0003685218520000102
和
Figure GDA0003685218520000103
为相应的拉格朗日乘子最优解。
附图3为地面用户信息的上行传输过程，地面用户8通过飞机11将数据发送至飞机网关5，再由飞机网关5将数据发送至数据中心9进行处理，或者地面用户8通过卫星10将数据发送至卫星网关4，再由卫星网关4将数据发送至数据中心9进行处理。
本发明的优点：
本发明所述方法在性能接近最优的同时，具有较低的算法复杂度，能够为空天地融合网络提供一种可行的用户接入方法和装置,能够为空天地融合网络提供一种可行的用户接入方案。
本发明重点研究空天地融合网络中5G信号覆盖区域外的边远地区信息传输问题。由于这类区域的地面用户无法与传统地面网络基站建立连接，为了将这些数据传递到数据中心进行处理，需要借助该区域上空的卫星、飞机进行数据传输。在地面用户-空中平台通信链路发生中断时，引入反馈重传机制，提出以最小化传输时延为目标的用户接入方案，并根据提出的启发式用户接入算法得到用户接入结果。
本发明所述的方法及装置适用于边远地区信息传输。
附图说明
图1为本发明所述的空天地融合网络中最小化传输时延的用户接入方法流程图；
图2为现有空天地融合网络架构示意图；
图3为现有地面用户信息的上行传输过程；
图4为空天地融合网络中引入飞机前后平均时延对比图；
图5为空天地融合网络中引入飞机前后吞吐量对比图；
图6为采用本发明所述的方法与现有同类技术实现通信的平均时延对比图；
图7为不同门限值下是否引入ARQ的吞吐量对比图；
附图标记说明：低轨卫星1、民航飞机2、地面用户节点3、卫星网关4、飞机网关5、G2A通信电路6、G2S通信电路7、地面用户8、数据中心9、卫星10、飞机11。
具体实施方式
实施方式一：参照图1说明本实施方式。本实施方式所述的一种空天地融合网络中最小化传输时延的用户接入方法，包括如下步骤：
步骤一：根据t时刻空天地融合网络的网络拓扑结构、发射功率、衰落损耗、收发端天线增益等参数，得到飞机接收端的信号接收形式并计算飞机端的接收信噪比；同时得到卫星接收端的信号接收形式并计算卫星端的接收信噪比；
步骤二：根据信道条件、接收端信噪比门限值，利用莱斯信道系数的概率密度函数计算得到地面用户-飞机通信链路的中断概率和地面用户-卫星通信链路的中断概率；其中，地面用户-飞机通信链路简称G2A，地面用户-卫星通信链路简称G2S；
步骤三：引入ARQ机制，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，并得到数据传输时延；同时根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，并得到数据传输时延；
步骤四：用户面对多种空中平台时，规划以最小化传输时延为目标的用户接入方法，将优化问题分解为一些独立的子问题，并根据本方法提出的启发式用户接入算法得到用户接入结果。
上述空天地融合网络的网络拓扑结构参见图2所示，该网络架构有两种通信链路，分别为G2A通信链路6和G2S通信电路7，其中G2A通信链路6是指任意地面用户节点3的数据通过民航飞机2和飞机网关5实现数据传输的链路；G2S通信链路7是指任意地面用户节点3的数据通过低轨卫星1和卫星网关4实现数据传输的链路。
其中，地面用户信息的上行传输过程参见图3所述，地面用户信息的上行传输过程分别为通过飞机上行传输和通过卫星上行传输，其中通过飞机上行传输的过程为：地面用户8通过飞机11将数据发送至飞机网关5，再由飞机网关5将数据发送至数据中心9进行处理；通过卫星上行传输的过程为：地面用户8通过卫星10将数据发送至卫星网关4，再由卫星网关4将数据发送至数据中心9进行处理。
上述步骤一中所述飞机和卫星端的接收信噪比为：
若时隙t地面用户m在飞机c的覆盖范围内，且用户通过G2A链路向飞机发射信号，则飞机端的信噪比的计算公式为：
Figure GDA0003685218520000121
其中，SNRm,c(t)表示时隙t内用户m向飞机c发射信号时飞机c处的信噪比，
Figure GDA0003685218520000122
表示用户m向飞机c的发射功率，
Figure GDA0003685218520000123
表示用户m的发射天线增益，
Figure GDA0003685218520000124
表示飞机c的接收天线增益，
Figure GDA0003685218520000125
表示用户m到飞机c的自由传输损耗，其中λc表示向飞机c发射的信号波长，dm,c(t)表示时隙t内用户m到飞机c的距离，La表示附加的损耗因子，hm,c(t)表示用户m向飞机c传输时的快衰落系数，N0表示噪声功率；
若时隙t地面用户m通过G2S链路向卫星发射信号，则卫星端的信噪比的计算公式为：
Figure GDA0003685218520000126
其中，SNRm,s(t)表示时隙t内用户m向卫星s发射信号时卫星s处的信噪比，
Figure GDA0003685218520000127
表示卫星s的接收天线增益，
Figure GDA0003685218520000128
表示用户m到卫星s的自由传输损耗，其中λs表示向卫星s的发射的信号波长，dm,s(t)表示时隙t内用户m到卫星s的距离；式中其它各参数的定义与前文相同。
步骤二中所述的G2A的中断概率的计算方法为：当所述接收端信噪比大于设定门限值γth时，认为数据在该条通信链路能够成功传输；因此，G2A链路中断概率
Figure GDA0003685218520000129
为：
Figure GDA00036852185200001210
其中，
Figure GDA00036852185200001211
表示时隙t内用户m向飞机c发射信号时链路的中断概率，Pr表示概率，γth表示设定的信噪比门限值。
由于|hm,c(t)|2服从莱斯分布，因此，|hm,c(t)|2的概率密度函数为：
Figure GDA0003685218520000131
其中，f|h|(h)表示|h|的概率密度函数，h表示概率密度函数的自变量，|h|表示快衰落因子h的模，K表示莱斯分布中的莱斯因子，I0(·)表示第一类0阶修正贝塞尔函数，e为自然常数；
G2A链路中断概率为：
Figure GDA0003685218520000132
其中，dh表示对自变量h的微分，Kc表示G2A链路中快衰落因子的模服从的莱斯分布因子，
Figure GDA0003685218520000133
为马库姆Q函数，其中a表示常数，b表示马库姆Q函数设定的积分下限。
步骤二中所述的G2S的中断概率的计算方法为：
Figure GDA0003685218520000134
其中，Ks表示G2S链路中快衰落因子的模服从的莱斯分布因子；式中其它各参数的定义与前文相同。
步骤三中在G2A链路中利用ARQ机制进行数据传输时，平均传输次数的计算方式为：
Figure GDA0003685218520000135
其中，
Figure GDA0003685218520000136
表示时隙t内用户m向飞机c发送数据的平均传输次数，Pr表示概率，N表示设置的ARQ机制传输次数上限，n表示某一次的传输序号，
Figure GDA0003685218520000137
表示时隙t内第n次传输成功，
Figure GDA0003685218520000138
表示时隙t内第n次传输失败。
所以，用户利用G2A链路进行数据传输的平均传输次数为：
Figure GDA0003685218520000141
时隙t内G2A链路的数据传输时延为：
Figure GDA0003685218520000142
其中，D1,m,c(t)为时隙t内用户m向飞机c发送数据的传输时延，
Figure GDA0003685218520000143
为时隙t用户m发送的数据包个数，R为每个数据包所包含的比特数，B为通信带宽；
同时时隙t内G2S链路的数据传输时延的计算方法为：
Figure GDA0003685218520000144
其中，D2,m,s(t)为时隙t内用户m向卫星s发送数据的传输时延，
Figure GDA0003685218520000145
表示时隙t内用户m向卫星s发送数据的平均传输次数；式中其它各参数的定义与前文相同。
步骤四中所述优化问题规划如下：
P1-1为最小化系统内所有用户总传输时延的优化问题；由于每个时隙的优化方法相对独立，将P1-1分解为T个独立的子问题，令第t个时刻的子问题为P1-2。
P1-1:
Figure GDA0003685218520000146
上式为P1-1的目标函数，表示最小化整个研究时间段内所有用户发送数据的传输时延。其中，a1,m,c(t)表示在时隙t内用户m是否接入飞机m的判决变量，a2,m,s(t)表示在时隙t内用户m是否接入卫星s的判决变量；
Figure GDA0003685218520000147
上式为P1-1的限制条件1，表示每个时刻每个用户只能与通信范围内的一架飞机或一个卫星相连。其中，ωm,c(t)表示时隙t内用户m是否在飞机c覆盖范围内的二进制判决变量，M表示系统内用户总数量，m表示某个用户序号；
Figure GDA0003685218520000148
上式为P1-1的限制条件2，限制了每个时隙内每架飞机可接入的用户数量上限，其中，
Figure GDA0003685218520000149
表示每架飞机每个时刻所能接入的最大用户数量，C表示系统内飞机总数量，c表示某架飞机序号；
Figure GDA0003685218520000151
上式为P1-1的限制条件3，限制了每个时隙内每颗卫星可接入的用户数量上限，其中，
Figure GDA0003685218520000152
表示每颗卫星每个时刻所能接入的最大用户数量，S表示为系统内卫星总数量，s表示某颗卫星序号；
α1,m,c(t)∈{0,1},1≤t≤T,1≤m≤M,1≤c≤C
上式为P1-1的限制条件4，表示在时隙t内用户m是否接入飞机m的判决变量α1,m,c(t)是一个二进制变量，其中，T表示整个时间段内的总时隙个数，t表示某一时隙序号；
α2,m,s(t)∈{0,1},1≤t≤T,1≤m≤M,1≤s≤S
上式为P1-1的限制条件5，表示在时隙t内用户m是否接入卫星s的判决变量α2,m,s(t)是一个二进制变量。
步骤四中所述启发式用户接入算法为：
首先对优化变量进行凸松弛，即：0≤α1,m,c(t)≤1和0≤α2,m,s(t)≤1。讨论时隙t内用户m在飞机c的通信范围内的情况，即ωm,c(t)＝1；利用拉格朗日对偶函数，将P1-2转化为：
Figure GDA0003685218520000153
其中，μ1,c(t)表示P1-1的限制条件2所对应的拉格朗日乘子，μ2,s(t)表示P1-1的限制条件3所对应的拉格朗日乘子，
Figure GDA0003685218520000154
和
Figure GDA0003685218520000155
为相应的拉格朗日乘子最优解。
本实施方式提出的启发式用户接入算法在性能接近最优的同时，具有较低的算法复杂度，能够为空天地融合网络提供一种可行的用户接入方法。
实施方式二：参照附图1说明本实施方式，本实施方式提供一种空天地融合网络中最小化传输时延的用户接入装置，包括：
用于根据t时刻空天地融合网络的网络拓扑结构、发射功率、衰落损耗、收发端天线增益等参数，得到飞机接收端的信号接收形式并计算飞机端的接收信噪比的装置；同时该装置还能得到卫星接收端的信号接收形式并计算卫星端的接收信噪比；
用于根据信道条件、接收端信噪比门限值，利用莱斯信道系数的概率密度函数计算得到地面用户-飞机通信链路的中断概率和地面用户-卫星通信链路的中断概率的装置；其中，地面用户-飞机通信链路简称G2A，地面用户-卫星通信链路简称G2S；
用于在引入ARQ机制后，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，得到数据传输时延的装置；同时该装置还能根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，得到数据传输时延；
用于在用户面对多种空中平台时，规划以最小化传输时延为目标的用户接入方法，将优化问题分解为一些独立的子问题，并根据本方法提出的启发式用户接入算法得到用户接入结果的装置。
所述飞机和卫星端的接收信噪比为：
若时隙t地面用户m在飞机c的覆盖范围内，且用户通过G2A链路向飞机发射信号，则飞机端的信噪比的计算公式为：
Figure GDA0003685218520000161
其中，SNRm,c(t)表示时隙t内用户m向飞机c发射信号时飞机c处的信噪比，
Figure GDA0003685218520000162
表示用户m向飞机c的发射功率，
Figure GDA0003685218520000163
表示用户m的发射天线增益，
Figure GDA0003685218520000164
表示飞机c的接收天线增益，
Figure GDA0003685218520000165
表示用户m到飞机c的自由传输损耗，其中λc表示向飞机c发射的信号波长，dm,c(t)表示时隙t内用户m到飞机c的距离。La表示附加的损耗因子，hm,c(t)表示用户m向飞机c传输时的快衰落系数，N0表示噪声功率；
若时隙t地面用户m通过G2S链路向卫星发射信号，则卫星端的信噪比的计算公式为：
Figure GDA0003685218520000166
其中，SNRm,s(t)表示时隙t内用户m向卫星s发射信号时卫星s处的信噪比，
Figure GDA0003685218520000167
表示卫星s的接收天线增益，
Figure GDA0003685218520000168
表示用户m到卫星s的自由传输损耗，其中λs表示向卫星s的发射的信号波长，dm,s(t)表示时隙t内用户m到卫星s的距离。式中其它各参数的定义与前文相同。
所述的G2A的中断概率的计算方法为：
当所述接收端信噪比大于设定门限值γth时，认为数据在该条通信链路能够成功传输；因此，G2A链路中断概率
Figure GDA0003685218520000171
为：
Figure GDA0003685218520000172
其中，
Figure GDA0003685218520000173
表示时隙t内用户m向飞机c发射信号时链路的中断概率，Pr表示概率，γth表示设定的信噪比门限值。
由于|hm,c(t)|2服从莱斯分布，因此，|hm,c(t)|2的概率密度函数为：
Figure GDA0003685218520000174
其中，f|h|(h)表示|h|的概率密度函数，h表示概率密度函数的自变量，|h|表示快衰落因子h的模，K表示莱斯分布中的莱斯因子，I0(·)表示第一类0阶修正贝塞尔函数，e为自然常数。
G2A链路中断概率为：
Figure GDA0003685218520000175
其中，dh表示对自变量h的微分，Kc表示G2A链路中快衰落因子的模服从的莱斯分布因子，
Figure GDA0003685218520000176
为马库姆Q函数，其中a表示常数，b表示马库姆Q函数设定的积分下限。
步骤二中所述的G2S的中断概率的计算方法为：
Figure GDA0003685218520000177
其中，Ks表示G2S链路中快衰落因子的模服从的莱斯分布因子。式中其它各参数的定义与前文相同。
所述的在G2A链路中利用ARQ机制进行数据传输时，平均传输次数的计算方式为：
Figure GDA0003685218520000181
其中，
Figure GDA0003685218520000182
表示时隙t内用户m向飞机c发送数据的平均传输次数，Pr表示概率，N表示设置的ARQ机制传输次数上限，n表示某一次的传输序号，
Figure GDA0003685218520000183
表示时隙t内第n次传输成功，
Figure GDA0003685218520000184
表示时隙t内第n次传输失败。
所以，用户利用G2A链路进行数据传输的平均传输次数为：
Figure GDA0003685218520000185
时隙t内G2A链路的数据传输时延为：
Figure GDA0003685218520000186
其中，D1,m,c(t)为时隙t内用户m向飞机c发送数据的传输时延，
Figure GDA0003685218520000187
为时隙t用户m发送的数据包个数，R为每个数据包所包含的比特数，B为通信带宽；
同时时隙t内G2S链路的数据传输时延的计算方法为：
Figure GDA0003685218520000188
其中，D2,m,s(t)为时隙t内用户m向卫星s发送数据的传输时延，
Figure GDA0003685218520000189
表示时隙t内用户m向卫星s发送数据的平均传输次数；
所述的优化问题规划如下：
P1-1为最小化系统内所有用户总传输时延的优化问题。由于每个时隙的优化方法相对独立，将P1-1分解为T个独立的子问题，令第t个时刻的子问题为P1-2。
P1-1:
Figure GDA00036852185200001810
上式为P1-1的目标函数，表示最小化整个研究时间段内所有用户发送数据的传输时延。其中，a1,m,c(t)表示在时隙t内用户m是否接入飞机m的判决变量，a2,m,s(t)表示在时隙t内用户m是否接入卫星s的判决变量；
Figure GDA0003685218520000191
上式为P1-1的限制条件1，表示每个时刻每个用户只能与通信范围内的一架飞机或一个卫星相连。其中，ωm,c(t)表示时隙t内用户m是否在飞机c覆盖范围内的二进制判决变量，M表示系统内用户总数量，m表示某个用户序号；
Figure GDA0003685218520000192
上式为P1-1的限制条件2，限制了每个时隙内每架飞机可接入的用户数量上限，其中，
Figure GDA0003685218520000193
表示每架飞机每个时刻所能接入的最大用户数量，C表示系统内飞机总数量，c表示某架飞机序号；
Figure GDA0003685218520000194
上式为P1-1的限制条件3，限制了每个时隙内每颗卫星可接入的用户数量上限，其中，
Figure GDA0003685218520000195
表示每颗卫星每个时刻所能接入的最大用户数量，S表示为系统内卫星总数量，s表示某颗卫星序号；
α1,m,c(t)∈{0,1},1≤t≤T,1≤m≤M,1≤c≤C
上式为P1-1的限制条件4，表示在时隙t内用户m是否接入飞机m的判决变量α1,m,c(t)是一个二进制变量，其中，T表示整个时间段内的总时隙个数，t表示某一时隙序号；
α2,m,s(t)∈{0,1},1≤t≤T,1≤m≤M,1≤s≤S
上式为P1-1的限制条件5，表示在时隙t内用户m是否接入卫星s的判决变量α2,m,s(t)是一个二进制变量。
所述的启发式用户接入算法为：
首先对优化变量进行凸松弛，即：0≤α1,m,c(t)≤1和0≤α2,m,s(t)≤1。讨论时隙t内用户m在飞机c的通信范围内的情况，即ωm,c(t)＝1。利用拉格朗日对偶函数，将P1-2转化为：
Figure GDA0003685218520000201
其中，μ1,c(t)表示P1-1的限制条件2所对应的拉格朗日乘子，μ2,s(t)表示P1-1的限制条件3所对应的拉格朗日乘子，
Figure GDA0003685218520000202
和
Figure GDA0003685218520000203
为相应的拉格朗日乘子最优解。
以下为本实施方式提出的启发式算法的可行性证明：
命题1：定义数据集合
Figure GDA0003685218520000204
集合中最小项所对应元素的序号为βm(t)＝arg min{Gm(t)}，则优化问题的最优解为
Figure GDA0003685218520000205
证明：
对于给定时刻和给定用户，最多有一个α1,m,c(t)或α2,m,s(t)为1。由于集合Gm(t)的每一项都是正值，根据线性函数的性质，βm(t)对应那一项的优化变量应为1，其它项为0。
证毕。
当ωm,c(t)＝0时，优化问题的求解方法与ωm,c(t)＝1时类似。
当ωm,c(t)＝1时，命题一提出了针对优化问题P1-1的一种启发式用户接入求解算法，形式较为简单。为了清晰，故给出证明，解释为什么最优解形式是命题一的形式。
当ωm,c(t)＝0时，表示用户m不在飞机c的覆盖范围内的情况，优化问题的求解方法与ωm,c(t)＝1时类似。具体地，命题1中的数据集合修改为
Figure GDA0003685218520000206
优化问题的最优解为
Figure GDA0003685218520000207
本实施方式提出的空天地融合网络中最小化传输时延的用户接入装置在性能接近最优的同时，具有较低的算法复杂度，能够为空天地融合网络提供一种可行的用户接入方法。
实施方式三：本实施方式是采用仿真的方式将本发明所述的空天地融合网络中最小化传输时延的用户接入方法与现有同类方法进行比较，进而验证本发明与现有同类技术相比较的优点。
表1所示为在仿真时用到的参数设定。
表1仿真参数设置
Figure GDA0003685218520000208
Figure GDA0003685218520000211
在上述仿真参数的情况下的仿真结果为：
附图4为空天地融合网络中引入飞机前后平均时延对比图；为了简便，在仿真时只考虑飞机和卫星同时处在地面用户节点通信范围内的情况。可以发现，当网络中只存在卫星时，平均时延最大；只存在飞机时，平均时延最小。当网络中同时存在飞机和卫星时，如果用户数量小于门限值，网络平均时延与只存在飞机时的情况相同；如果超过该门限值，平均时延上升，时延值介于前两种场景之间。这说明了当飞机和卫星同时能够为地面用户提供服务时，用户优先选择接入飞机。当飞机无法满足用户通信需求时，用户才会选择接入卫星。
本实施方式提供了一种空天地融合网络中最小化传输时延的用户接入方法，通过该方法引入飞机后，在用户数量少于150个时，平均延时降低了125ms左右，在用户数量高于150个，少于300个时，平均延时降低了100ms左右。
附图5为空天地融合网络中引入飞机前后吞吐量对比图；可以发现，三种场景下，当用户数量小于各自门限值时，吞吐量几乎随着用户数量的增加呈线性增长；当超过门限值后，吞吐量趋于稳定。可以发现，由于引入了飞机作为空中平台，网络的饱和点上升，能够为更多用户提供服务，吞吐量远高于具有单一类型空中平台场景的吞吐量因此，在引入飞机作为空中平台之后，网络的饱和点上升，在用户数量达到门限值之后吞吐量能够趋于稳定，为更多客户提供服务。
附图6为采用本发明所述的方法与现有同类技术实现通信的平均时延对比图；即：本发明所述的方法、现有intlinprog算法发以及随机接入算法。其中，intlinprog算法利用了现有优化问题工具包，能够求得最优解，它对应的是性能上界。可以发现，启发式算法和intlinprog算法下的平均时延远小于随机接入算法。此外，表2对比了网络中有8架飞机和4颗卫星时，本发明所述方法和intlinprog算法的CPU运行时间。可以发现，本发明所述方法在性能接近上界的情况下，还具有较低的计算复杂度。
表2当C＝8,S＝4时CPU运行时间对比
用户数量 本发明所述方法 intlinprog算法
M＝100 3.4442s 3.6613s
M＝200 6.1808s 6.4362s
M＝300 8.9783s 9.4611s
根据图6和表2，本发明所述的方法和现有intlinprog算法下的平均时延远小于随机接入算法，并且本发明所述的方法在性能接近上界的情况下，还具有较低的计算复杂度。
附图7为不同门限值下是否引入ARQ的吞吐量对比图；由于飞机在提供地面通信的过程中还要受到地面导航站、机场等信号的干扰，因此，设定飞机端的接收信噪比门限值高于卫星端的门限值。可以发现，对于给定的信噪比门限值，由于采用ARQ机制，增加了用户-空中平台的通信链路成功建立的概率，因此吞吐量更高。此外，当信噪比门限值更高时，有ARQ机制要比没有ARQ机制的性能更好。这是因为当门限值更高时，用户-空中平台的通信链路更容易发生中断。

Claims (10)

1.一种空天地融合网络中最小化传输时延的用户接入方法，其特征在于，所述方法包括如下步骤：
步骤一：根据t时刻空天地融合网络的网络拓扑结构、发射功率、衰落损耗、收发端天线增益参数，得到飞机接收端的信号接收形式并计算飞机端的接收信噪比；同时得到卫星接收端的信号接收形式并计算卫星端的接收信噪比；
步骤二：根据信道条件、接收端信噪比门限值，利用莱斯信道系数的概率密度函数计算得到地面用户-飞机通信链路的中断概率和地面用户-卫星通信链路的中断概率；其中，地面用户-飞机通信链路简称G2A，地面用户-卫星通信链路简称G2S；
步骤三：引入ARQ机制，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，得到数据传输时延；同时根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，得到数据传输时延；
步骤四：用户面对多种空中平台时，规划以最小化传输时延为目标的用户接入方法，将优化问题分解为一些独立的子问题，并根据启发式用户接入算法得到用户接入结果。
2.根据权利要求1所述的一种空天地融合网络中最小化传输时延的用户接入方法，其特征为，步骤一中所述飞机和卫星端的接收信噪比采用下述方法获得：
若时隙t地面用户m在飞机c的覆盖范围内，且用户通过G2A链路向飞机发射信号，则飞机端的信噪比的计算公式为：
Figure FDA0003697282060000011
其中，SNRm,c(t)表示时隙t内用户m向飞机c发射信号时、飞机c处的信噪比，
Figure FDA0003697282060000012
表示用户m向飞机c的发射功率，
Figure FDA0003697282060000013
表示用户m的发射天线增益，
Figure FDA0003697282060000014
表示飞机c的接收天线增益，
Figure FDA0003697282060000015
表示时隙t内用户m到飞机c的自由传输损耗，其中λc表示向飞机c发射的信号波长，dm,c(t)表示时隙t内用户m到飞机c的距离，La表示附加的损耗因子，hm,c(t)表示时隙t内用户m向飞机c传输时的快衰落系数，N0表示噪声功率；
若时隙t地面用户m通过G2S链路向卫星发射信号，则卫星端的信噪比的计算公式为：
Figure FDA0003697282060000021
其中，SNRm,s(t)表示时隙t内用户m向卫星s发射信号时、卫星s处的信噪比，
Figure FDA0003697282060000022
表示用户m向卫星s的发射功率
Figure FDA0003697282060000023
表示卫星s的接收天线增益，
Figure FDA0003697282060000024
表示用户m到卫星s的自由传输损耗，其中λs表示向卫星s的发射的信号波长，dm,s(t)表示时隙t内用户m到卫星s的距离，
Figure FDA0003697282060000025
表示时隙t内用户m向飞机c传输时的快衰落系数。
3.根据权利要求2所述的一种空天地融合网络中最小化传输时延的用户接入方法，其特征为，步骤二中所述的G2A的中断概率的计算方法为：当所述接收端信噪比大于设定门限值γth时，认为数据在该条通信链路能够成功传输；因此，G2A链路中断概率
Figure FDA0003697282060000026
为：
Figure FDA0003697282060000027
其中，
Figure FDA0003697282060000028
表示时隙t内用户m向飞机c发射信号时链路的中断概率，Pr表示概率，γth表示设定的信噪比门限值，λc表示向卫星s的发射的信号波长；
由于|hm,c(t)|2服从莱斯分布，因此，|hm,c(t)|2的概率密度函数为：
Figure FDA0003697282060000029
其中，f|h|(h)表示|h|的概率密度函数，h表示概率密度函数的自变量，|h|表示快衰落因子h的模，K表示莱斯分布中的莱斯因子，I0(·)表示第一类0阶修正贝塞尔函数，e为自然常数；
G2A链路中断概率为：
Figure FDA00036972820600000210
其中，dh表示对自变量h的微分，Kc表示G2A链路中快衰落因子的模服从的莱斯分布因子，
Figure FDA00036972820600000211
为马库姆Q函数，其中a表示常数，b表示马库姆Q函数设定的积分下限；
步骤二中所述的G2S的中断概率的计算方法为：
Figure FDA0003697282060000031
其中，Ks表示G2S链路中快衰落因子的模服从的莱斯分布因子。
4.根据权利要求3所述的一种空天地融合网络中最小化传输时延的用户接入方法，其特征为，步骤三中所述的在G2A链路中利用ARQ机制进行数据传输时，平均传输次数的计算方式为：
Figure FDA0003697282060000032
其中，
Figure FDA0003697282060000033
表示时隙t内用户m向飞机c发送数据的平均传输次数，Pr表示概率，N表示设置的ARQ机制传输次数上限，n表示某一次的传输序号，
Figure FDA0003697282060000034
表示时隙t内第n次传输成功，
Figure FDA0003697282060000035
表示时隙t内第n次传输失败；
所以，用户利用G2A链路进行数据传输的平均传输次数为：
Figure FDA0003697282060000036
时隙t内G2A链路的数据传输时延为：
Figure FDA0003697282060000037
其中，D1,m,c(t)为时隙t内用户m向飞机c发送数据的传输时延，
Figure FDA0003697282060000038
为时隙t用户m发送的数据包个数，R为每个数据包所包含的比特数，B为通信带宽；
时隙t内G2S链路的数据传输时延的计算方法为：
Figure FDA0003697282060000039
其中，D2,m,s(t)为时隙t内用户m向卫星s发送数据的传输时延，
Figure FDA00036972820600000310
表示时隙t内用户m向卫星s发送数据的平均传输次数。
5.根据权利要求4所述的一种空天地融合网络中最小化传输时延的用户接入方法，其特征为，步骤四中所述的优化问题规划如下：
P1-1为最小化系统内所有用户总传输时延的优化问题；由于每个时隙的优化方法相对独立，将P1-1分解为T个独立的子问题，令第t个时刻的子问题为P1-2，
Figure FDA0003697282060000041
上式为P1-1的目标函数，表示最小化整个研究时间段内所有用户发送数据的传输时延，其中，a1,m,c(t)表示在时隙t内用户m是否接入飞机m的判决变量，a2,m,s(t)表示在时隙t内用户m是否接入卫星s的判决变量；
Figure FDA0003697282060000042
上式为P1-1的限制条件1，表示每个时刻每个用户只能与通信范围内的一架飞机或一个卫星相连，其中，ωm,c(t)表示时隙t内用户m是否在飞机c覆盖范围内的二进制判决变量，M表示系统内用户总数量，m表示用户序号；
Figure FDA0003697282060000043
上式为P1-1的限制条件2，限制了每个时隙内每架飞机可接入的用户数量上限，其中，
Figure FDA0003697282060000044
表示每架飞机每个时刻所能接入的最大用户数量，C表示系统内飞机总数量，c表示飞机序号；
Figure FDA0003697282060000045
上式为P1-1的限制条件3，限制了每个时隙内每颗卫星可接入的用户数量上限，其中，
Figure FDA0003697282060000046
表示每颗卫星每个时刻所能接入的最大用户数量，S表示为系统内卫星总数量，s表示某颗卫星序号；
α1,m,c(t)∈{0,1},1≤t≤T,1≤m≤M,1≤c≤C
上式为P1-1的限制条件4，表示在时隙t内用户m是否接入飞机m的判决变量α1,m,c(t)是一个二进制变量，其中，T表示整个时间段内的总时隙个数，t表示某一时隙序号；
α2,m,s(t)∈{0,1},1≤t≤T,1≤m≤M,1≤s≤S
上式为P1-1的限制条件5，表示在时隙t内用户m是否接入卫星s的判决变量α2,m,s(t)是一个二进制变量；
步骤四中所述的启发式用户接入算法为：
首先对优化变量进行凸松弛，即：0≤α1,m,c(t)≤1和0≤α2,m,s(t)≤1，讨论时隙t内用户m在飞机c的通信范围内的情况，即ωm,c(t)＝1；利用拉格朗日对偶函数，将P1-2转化为：
Figure FDA0003697282060000051
其中，μ1,c(t)表示P1-1的限制条件2所对应的拉格朗日乘子，μ2,s(t)表示P1-1的限制条件3所对应的拉格朗日乘子，
Figure FDA0003697282060000052
和
Figure FDA0003697282060000053
为相应的拉格朗日乘子最优解。
6.一种空天地融合网络中最小化传输时延的用户接入装置，其特征在于，所述装置包括：
信噪比计算装置，所述信噪比计算装置用于根据t时刻空天地融合网络的网络拓扑结构、发射功率、衰落损耗、收发端天线增益参数，得到飞机接收端的信号接收形式并计算飞机端的接收信噪比；同时该装置还用于得到卫星接收端的信号接收形式并计算卫星端的接收信噪比；
中断概率计算装置，所述中断概率计算装置用于根据信道条件、接收端信噪比门限值，利用莱斯信道系数的概率密度函数计算得到地面用户-飞机通信链路的中断概率和地面用户-卫星通信链路的中断概率；其中，地面用户-飞机通信链路简称G2A，地面用户-卫星通信链路简称G2S；
数据传输时延装置，所述数据传输时延装置用于在引入ARQ机制后，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，得到数据传输时延；同时该装置还用于根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，得到数据传输时延；
结果确定装置，所述结果确定装置用于在用户面对多种空中平台时，规划以最小化传输时延为目标的用户接入方法，将优化问题分解为一些独立的子问题，并根据启发式用户接入算法得到用户接入结果。
7.根据权利要求6所述的一种空天地融合网络中最小化传输时延的用户接入装置，其特征在于，所述信噪比计算装置还包括：飞机端信噪比计算模块和卫星端信噪比计算模块，其中，飞机端信噪比计算模块用于得到飞机接收端的信号接收形式并计算飞机端的接收信噪比，具体采用下述公式计算获得：
若时隙t地面用户m在飞机c的覆盖范围内，且用户通过G2A链路向飞机发射信号，则飞机端的接收信噪比的计算公式为：
Figure FDA0003697282060000054
其中，SNRm,c(t)表示时隙t内用户m向飞机c发射信号时飞机c处的信噪比，
Figure FDA0003697282060000061
表示用户m向飞机c的发射功率，
Figure FDA0003697282060000062
表示用户m的发射天线增益，
Figure FDA0003697282060000063
表示飞机c的接收天线增益，
Figure FDA0003697282060000064
表示用户m到飞机c的自由传输损耗，其中λc表示向飞机c发射的信号波长，dm,c(t)表示时隙t内用户m到飞机c的距离，La表示附加的损耗因子，hm,c(t)表示用户m向飞机c传输时的快衰落系数，N0表示噪声功率；
所述的卫星端信噪比模块用于得到卫星端的信号接收形式并计算卫星端的接收信噪比，具体采用下述公式计算获得：
若时隙t地面用户m通过G2S链路向卫星发射信号，则卫星端的接收信噪比的计算公式为：
Figure FDA0003697282060000065
其中，SNRm,s(t)表示时隙t内用户m向卫星s发射信号时卫星s处的信噪比，
Figure FDA0003697282060000066
表示卫星s的接收天线增益，
Figure FDA0003697282060000067
表示用户m到卫星s的自由传输损耗，其中λs表示向卫星s的发射的信号波长，dm,s(t)表示时隙t内用户m到卫星s的距离，式中其它各参数的定义与前文相同。
8.根据权利要求7所述的一种空天地融合网络中最小化传输时延的用户接入装置，其特征在于，所述中断概率计算装置还包括：G2A的中断概率的计算模块和G2S的中断概率的计算模块，其中G2A的中断概率的计算模块用于计算得到地面用户-飞机通信链路的中断概率，其中地面用户-飞机通信链路简称G2A，所述的G2A的中断概率的计算模块计算G2A中断概率具体采用下述公式计算获得：
当所述接收端信噪比大于设定门限值γth时，认为数据在该条通信链路能够成功传输；因此，G2A链路中断概率
Figure FDA0003697282060000068
为：
Figure FDA0003697282060000069
其中，
Figure FDA00036972820600000610
表示时隙t内用户m向飞机c发射信号时链路的中断概率，Pr表示概率，γth表示设定的信噪比门限值；
由于|h|服从莱斯分布，因此，|h|的概率密度函数为：
Figure FDA0003697282060000071
其中，f|h|(h)表示|h|的概率密度函数，h表示概率密度函数的自变量，|h|表示快衰落因子h的模，K表示莱斯分布中的莱斯因子，I0(·)表示第一类0阶修正贝塞尔函数；
G2A链路中断概率为：
Figure FDA0003697282060000072
其中，dh表示对自变量h的微分，Kc表示G2A链路中快衰落因子的模服从的莱斯分布因子，
Figure FDA0003697282060000073
为马库姆Q函数，其中a表示常数，b表示马库姆Q函数设定的积分下限；
所述的G2S的中断概率的计算模块用于计算得到地面用户-卫星通信链路的中断概率，其中，地面用户-卫星通信链路简称G2S，所述的G2S的中断概率的计算模块计算G2S中断概率具体采用下述公式计算获得：
Figure FDA0003697282060000074
其中，Ks表示G2S链路中快衰落因子的模服从的莱斯分布因子。
9.根据权利要求8所述的一种空天地融合网络中最小化传输时延的用户接入装置，其特征在于，所述数据传输时延装置还包括：G2A平均传输次数模块和G2S平均传输次数模块，其中G2A平均传输次数模块用于在引入ARQ机制后，根据数据包传输次数上限计算用户-飞机的数据包平均发送次数，得到数据传输时延，所述的G2A平均传输次数模块计算传输次数具体采用下述公式计算获得：
Figure FDA0003697282060000075
其中，
Figure FDA0003697282060000081
表示时隙t内用户m向飞机c发送数据的平均传输次数，Pr表示概率，N表示设置的ARQ机制传输次数上限，n表示某一次的传输序号，
Figure FDA0003697282060000082
表示时隙t内第n次传输成功，
Figure FDA0003697282060000083
表示时隙t内第n次传输失败；
所以，用户利用G2A链路进行数据传输的平均传输次数为：
Figure FDA0003697282060000084
时隙t内G2A链路的数据传输时延为：
Figure FDA0003697282060000085
其中，D1,m,c(t)为时隙t内用户m向飞机c发送数据的传输时延，
Figure FDA0003697282060000086
为时隙t用户m发送的数据包个数，R为每个数据包所包含的比特数，B为通信带宽；
所述的G2S平均传输次数模块用于在引入ARQ机制后，根据数据包传输次数上限计算用户-卫星的数据包平均发送次数，得到数据传输时延，所述的G2S平均传输次数模块计算时隙t内G2S链路的数据传输时延具体采用下述公式计算获得：
Figure FDA0003697282060000087
其中，D2,m,s(t)为时隙t内用户m向卫星s发送数据的传输时延，
Figure FDA0003697282060000088
表示时隙t内用户m向卫星s发送数据的平均传输次数。
10.根据权利要求9所述的一种空天地融合网络中最小化传输时延的用户接入装置，其特征在于，所述的结果确定装置还包括：优化问题模块和得到结果模块，其中优化问题模块用于将优化问题分解为一些独立的子问题，所述的优化问题模块优化问题的规划如下：
P1-1为最小化系统内所有用户总传输时延的优化问题，由于每个时隙的优化方法相对独立，将P1-1分解为T个独立的子问题，令第t个时刻的子问题为P1-2；
Figure FDA0003697282060000089
上式为P1-1的目标函数，表示最小化整个研究时间段内所有用户发送数据的传输时延；其中，a1,m,c(t)表示在时隙t内用户m是否接入飞机m的判决变量，a2,m,s(t)表示在时隙t内用户m是否接入卫星s的判决变量；
Figure FDA0003697282060000091
上式为P1-1的限制条件1，表示每个时刻每个用户只能与通信范围内的一架飞机或一个卫星相连；其中，ωm,c(t)表示时隙t内用户m是否在飞机c覆盖范围内的二进制判决变量，M表示系统内用户总数量，m表示某个用户序号；
Figure FDA0003697282060000092
上式为P1-1的限制条件2，限制了每个时隙内每架飞机可接入的用户数量上限，其中，
Figure FDA0003697282060000093
表示每架飞机每个时刻所能接入的最大用户数量，C表示系统内飞机总数量，c表示某架飞机序号；
Figure FDA0003697282060000094
上式为P1-1的限制条件3，限制了每个时隙内每颗卫星可接入的用户数量上限，其中，
Figure FDA0003697282060000095
表示每颗卫星每个时刻所能接入的最大用户数量，S表示为系统内卫星总数量，s表示某颗卫星序号；
α1,m,c(t)∈{0,1},1≤t≤T,1≤m≤M,1≤c≤C
上式为P1-1的限制条件4，表示在时隙t内用户m是否接入飞机m的判决变量α1,m,c(t)是一个二进制变量，其中，T表示整个时间段内的总时隙个数，t表示某一时隙序号；
α2,m,s(t)∈{0,1},1≤t≤T,1≤m≤M,1≤s≤S
上式为P1-1的限制条件5，表示在时隙t内用户m是否接入卫星s的判决变量α2,m,s(t)是一个二进制变量；
所述的得到结果模块用于根据启发式用户接入算法得到用户接入结果，所述的得到结果模块中的启发式用户接入算法为：
首先对优化变量进行凸松弛，即：0≤α1,m,c(t)≤1和0≤α2,m,s(t)≤1，讨论时隙t内用户m在飞机c的通信范围内的情况，即ωm,c(t)＝1，利用拉格朗日对偶函数，将P1-2转化为：
Figure FDA0003697282060000101
其中，μ1,c(t)表示P1-1的限制条件2所对应的拉格朗日乘子，μ2,s(t)表示P1-1的限制条件3所对应的拉格朗日乘子，
Figure FDA0003697282060000102
和
Figure FDA0003697282060000103
为相应的拉格朗日乘子最优解。
CN202110615508.2A 2021-06-02 2021-06-02 空天地融合网络中最小化传输时延的用户接入方法及装置 Active CN113395101B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110615508.2A CN113395101B (zh) 2021-06-02 2021-06-02 空天地融合网络中最小化传输时延的用户接入方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110615508.2A CN113395101B (zh) 2021-06-02 2021-06-02 空天地融合网络中最小化传输时延的用户接入方法及装置

Publications (2)

Publication Number Publication Date
CN113395101A CN113395101A (zh) 2021-09-14
CN113395101B true CN113395101B (zh) 2022-07-22

Family

ID=77619982

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110615508.2A Active CN113395101B (zh) 2021-06-02 2021-06-02 空天地融合网络中最小化传输时延的用户接入方法及装置

Country Status (1)

Country Link
CN (1) CN113395101B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114826940B (zh) * 2022-04-25 2023-04-25 北京理工大学 一种天地融合网络的效能评估方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105939184A (zh) * 2016-03-04 2016-09-14 哈尔滨工业大学深圳研究生院 基于无迹卡尔曼滤波的空天DTN网络bundle传输时延估计算法
CN110300059A (zh) * 2019-07-19 2019-10-01 中国科学技术大学 最小化信息在轨传输时延的方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102711262B (zh) * 2012-07-02 2014-10-08 北京交通大学 一种异构网络下行链路资源分配方法
EP3879881A1 (en) * 2014-11-26 2021-09-15 IDAC Holdings, Inc. Beam switching in wireless systems
CN105704824B (zh) * 2016-01-18 2018-11-30 中国科学院计算技术研究所 一种无线网络的多维资源分配的方法
CN107294593B (zh) * 2017-06-21 2020-05-01 哈尔滨工业大学深圳研究生院 基于geo骨干中继的深空下行链路多跳传输方法及系统
US10405215B2 (en) * 2017-06-27 2019-09-03 Veniam, Inc. Self-organized fleets of autonomous vehicles to optimize future mobility and city services
US10419106B1 (en) * 2017-09-29 2019-09-17 Amazon Technologies, Inc. Satellite-based content delivery network (CDN) in an extraterrestrial environment
US10455475B1 (en) * 2018-05-29 2019-10-22 Hughes Network Systems, Llc Inter-layer communications in wireless networks including a high latency connection
CN109309524B (zh) * 2018-09-25 2021-06-08 中国人民解放军空军工程大学 一种基于ieee802.11的低轨道卫星至同步轨道数据中继卫星的接入方法
CN111526090B (zh) * 2020-07-06 2020-12-25 北京大学深圳研究生院 一种天地一体化信息网络的路由方法
CN112040510B (zh) * 2020-09-04 2023-10-20 中国电子科技集团公司第二十研究所 基于QoS约束和负载均衡的卫星网络接入选择方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105939184A (zh) * 2016-03-04 2016-09-14 哈尔滨工业大学深圳研究生院 基于无迹卡尔曼滤波的空天DTN网络bundle传输时延估计算法
CN110300059A (zh) * 2019-07-19 2019-10-01 中国科学技术大学 最小化信息在轨传输时延的方法

Also Published As

Publication number Publication date
CN113395101A (zh) 2021-09-14

Similar Documents

Publication Publication Date Title
JP5571245B2 (ja) 航空機および船でのインターネット・サービスのためのハイブリッド衛星メッシュ・ネットワーク・システム
US11509391B2 (en) Adaptive multiple access scheme in integrated access and backhaul networks
She et al. UAV-assisted uplink transmission for ultra-reliable and low-latency communications
Huang et al. Multi-hop D2D communications with network coding: From a performance perspective
Shi et al. Joint gateway selection and resource allocation for cross-tier communication in space-air-ground integrated IoT networks
Lagunas et al. Carrier allocation for hybrid satellite-terrestrial backhaul networks
CN105072676A (zh) 基于tdma协议的航空自组网功率控制方法
CN116112060A (zh) 基于缓冲中继的空天地通信系统资源分配方法及装置
Hassan et al. Joint FSO fronthaul and millimeter-wave access link optimization in cloud small cell networks: A statistical-QoS aware approach
Guo et al. Throughput maximization in a UAV-enabled two-way relaying system with multi-pair users
CN113395101B (zh) 空天地融合网络中最小化传输时延的用户接入方法及装置
Kakar et al. UAV-aided multi-way communications
Karavolos et al. HST-NNC: A novel hybrid satellite-terrestrial communication with NOMA and network coding systems
Zhang et al. Statistical QoS provisioning for UAV-enabled emergency communication networks
Wang et al. Performance analysis of cooperative aerial base station-assisted networks with non-orthogonal multiple access
Alagha Satellite air interface evolutions in the 5G and IoT era
Bliss et al. Uplink NOMA for heterogeneous NTNs with LEO satellites and high-altitude platform relays
US11799539B2 (en) Relaying in a wireless communication network
Hu et al. Optimal resource allocation in ground wireless networks supporting unmanned aerial vehicle transmissions
CN102891736B (zh) 一种应用于空天信息网络的异构链路帧长自适应方法
CN116390117A (zh) 无人机中继上行通信的资源分配与轨迹优化方法及装置
Joo et al. Dynamic cross-layer transmission control for station-assisted satellite networks
CN114531188A (zh) 基于非正交多址的星地融合网络缓存方法
FI20205794A1 (en) Point-to-multipoint transmission using a multi-antenna system
CN116455456B (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
GR01 Patent grant
GR01 Patent grant