Performance Analysis of Satellite-Aerial-Terrestrial Integrated Network Based on Uplink NOMA Technology
-
摘要: 针对光电混合的星空地融合网络上行链路,该文研究了多天线波束成形技术和上行非正交多址接入(NOMA)技术相结合的系统遍历和速率性能。首先,在无人机采用多天线和上行NOMA技术条件下,为实现系统和速率最大化,提出了一种基于统计信道状态信息的波束成形方案。接着,假设卫星-无人机链路采用自由空间光链路且服从伽马-伽马衰落,无人机-地面用户链路采用射频链路且服从相关瑞利衰落,推导了系统和速率的闭合表达式。最后,通过数值仿真验证了理论分析的正确性。仿真结果表明,与正交多址接入(OMA)方案相比,所提方案提高了系统性能,并且与基准波束成形(BF)方案相比,所提方案具有更好的性能优势。Abstract: The system Ergodic Sum Rate (ESR) performance of a mixed Free Space Optical/ Radio Frequency (FSO/RF) Satellite-Aerial-Terrestrial Integrated Network (SATIN) combining multi-antenna beamforming scheme and uplink Non-Orthogonal Multiple Access (NOMA) technology is investigated in this paper. Firstly, a beamforming scheme using statistical channel state information is proposed to maximize the ESR of the considered system, where multi-antenna UAV is deployed with uplink NOMA technology. Secondly, under the assumption that the satellite-UAV link with FSO undergoes Gamma-Gamma fading while the UAV-terrestrial links adopted RF experience correlated Rayleigh fading, the closed-form expression for ESR of the considered SATIN is derived. Finally, numerical results are provided to validate the theoretical analysis. The simulation results show that, compared with the Orthogonal Multiple Access (OMA) scheme, the proposed scheme improves the system ESR, and the proposed scheme has the performance superiority over the other benchmark beamforming scheme.
-
1. 引言
第6代移动通信(6G)网络旨在扩展通信系统的覆盖范围和实现万物互联,而卫星通信具有覆盖范围广、通信容量大、能实现偏远地区网络接入等优点,是实现这一愿景的关键技术之一[1]。然而,卫星与地面用户建立互联网连接时,受到遮蔽效应、路径损耗、衰落等影响,链路难以保持视距传输。为提高卫星通信的可靠性,可利用中继协助卫星和地面用户之间进行通信[2]。但考虑实际情况,在人口稀少地区建立大量地面中继成本高、回报低,而无人机(Unmanned Aerial Vehicle, UAV)成本低、效益高、灵活性高,可作为空中中继。因此,将无人机辅助卫星通信的星空地融合(Satellite-Aerial-Terrestrial Integrated Network, SATIN)网络逐渐成为研究热点[3]。例如,文献[3]分析了采用放大转发协议的无人机中继辅助SATIN系统能效性能。
随着互联网快速发展,现有的射频(Radio Frequency, RF)通信系统已经难以满足高速率业务发展要求。自由空间光(Free-Space Optical, FSO)通信相比于RF通信有很多优点,如未分配的频谱、容量大、功耗低、抗干扰能力强、安全性高,光电混合的网络得到了大量研究[4,5]。将RF/FSO混合链路应用于星空地融合网络也做了初步的探索性研究[6],如文献[6]研究了多用户星空地融合网络遍历容量性能。在RF/FSO混合星空地融合网络中,地面用户通过RF链路接入中继,而FSO被用于连接中继和卫星,该方案尤其适合为偏远地区地面用户提供高速接入。
非正交多址接入(Non-Orthogonal Multiple Access, NOMA)技术作为提高用户接入数量和资源利用率的关键技术之一,在卫星通信领域引起了学者的广泛关注[7]。例如,文献[7]卫星通信系统采用NOMA技术,在满足发射功率和用户服务质量的约束条件下,实现了遍历容量的最大化。需要指出的是,文献[7]主要针对下行NOMA卫星通信系统,而上行NOMA技术在卫星通信网络中的研究成果极少。同时,现有文献中对NOMA技术的研究工作通常建立在完美串行干扰消除(Successive Interference Canceler, SIC)条件下[7,8],而非完美SIC更符合卫星通信实际应用场景,只有极少数文献对其做了相关研究。如文献[9]针对单跳下行NOMA卫星通信系统,分析了非完美SIC条件下系统中断概率和遍历容量性能。
此外,星空地融合网络中采用多天线波束成形(BeamForming, BF)技术有利于提高系统的频谱效率与系统容量[3]。将BF技术与NOMA技术相结合可以同时利用两者的优势[10]。文献[10]针对多波束卫星服务多个NOMA用户组场景,在满足功率约束条件下,通过设计波束成形方案解决了系统遍历容量最大化的问题。文献[10]采用BF技术提高了NOMA系统性能,但其BF方案是基于准确信道状态信息(Channel State Information, CSI)条件,然而实际信号传输时易受到干扰且衰减严重,同时准确CSI需实时反馈,反馈开销较大。相比于准确CSI,统计CSI更易获得。
总体来看,基于RF/FSO混合链路的星空地融合网络的研究对提高系统速率具有重要意义,然而,目前学术界对其开展的研究工作较少。此外,现有研究表明NOMA技术能提高卫星通信资源利用率和系统性能,但目前未有文献研究上行NOMA技术在星空地融合网络中应用。本文针对RF/FSO混合星空地融合网络,研究了多天线波束成形技术和NOMA技术相结合的系统遍历和速率性能。具体而言,系统采用了放大转发(Amplify-and-Forward, AF)协议的UAV作为空中中继,以辅助地面多用户与卫星通信。首先,在无人机采用多天线和上行NOMA技术条件下,为实现系统遍历和速率最大化,与文献[10]不同,本文提出了一种基于统计CSI的波束成形方案。接着,假设卫星-无人机链路采用自由空间光链路且服从伽马-伽马衰落,无人机-地面用户链路采用射频链路且服从相关瑞利衰落,与文献[8]完美SIC条件不同,本文在考虑NOMA非完美SIC条件下,推导了系统遍历和速率的闭合表达式。最后,数值仿真结果验证了理论分析的正确性,同时仿真表明所提方案具有更好的性能。
本文符号说明:
(⋅)H 表示共轭转置,‖⋅‖ 表示欧几里得范数,IN 表示N×N 单位矩阵,CN(μ,σ2) 表示均值为μ ,方差为σ2 的复高斯分布,E[⋅] 表示数学期望,Gm,np,q[⋅|⋅] 表示Meijer-G函数(文献[11], Eq.(9.301))。2. 系统模型
如图1所示,本文研究了基于上行NOMA技术的光电混合SATIN,其中,
K 个地面用户(D)与卫星(S)之间进行通信,由于遮挡效应,忽略卫星和地面用户之间的直达链路影响[6],因此利用无人机作为中继(R)协作地面用户与卫星进行通信,且无人机具有静态悬停特性。假设卫星和地面用户为单天线,无人机配置N 根天线。所有地面用户分布于无人机的覆盖范围内,且信道相关性强和信道增益差异性大的用户被分为一组[10],K 个用户被分配到M 个NOMA组中,每组Kg 个用户且采用NOMA技术接入。2.1 信道模型
2.1.1 FSO链路
考虑到频谱资源的短缺和用户容量的需求,假设卫星-无人机链路用FSO进行连接,因路径损耗、信道衰落和天线增益的影响,FSO链路信道表示为[12]
hFSO,k=ℓFSOgFSO,k (1) 其中,
ℓFSO=GFt+GFr−AtenFS−AtenAtm−LLens−Me ,其中,GFt 和GFr 分别为发送和接收增益,AtenFS ,AtenAtm ,LLens ,Me 分别为自由空间损耗、大气衰减、透镜损耗、系统边缘。因无人机-卫星的FSO链路传输受多普勒效应、大气湍流和指向误差等重要因素影响,文献[13]对FSO传输特性作了相关研究,但模型建立较复杂。根据文献[14]中FSO链路的建模,设信道增益gFSO,k 服从伽马-伽马分布,该分布描述了弱到强湍流的闪烁统计信息,包含了大尺度和小尺度湍流涡旋未考虑指向误差[6],其概率密度函数(Probability Density Function, PDF)表达式为f|gFSO|(x)=αβτ2A0Γ(α)Γ(β)G3,01,3[αβA0x|τ2τ2−1,α−1,β−1] (2) 其中,
A0 为指向损失,τ 为等效波束半径与指向误差偏移量之间的比率,α ,β 分别为与大气环境有关的大尺度和小尺度单元的有效数量,受到高度的影响[14]。2.1.2 RF链路
假设无人机-地面用户的链路为RF链路,考虑天线阵列的相关性,设信道模型服从相关瑞利分布[15],信道模型表示为
hk=ℓkgk (3) 其中,
ℓk 为自由空间路径损耗,表示为ℓk=GRt+GRr+10αhlgdk+20lg(4π/4π cc)+20lgfc [6],其中,GRt ,GRr 分别表示发送和接收增益,c ,fc ,dk ,αh 分别表示光速、载波频率、无人机和地面用户之间的距离以及路径损耗因子。根据文献[16]的信道散射模型,信道增益矢量gk 表示为gk=1√LL∑l=1ρlb(θd) (4) 在散射模型中,设有
L 路信号分布在角度θd∈[ˉθd−Δθd/Δθd22,ˉθd+Δθd/Δθd22] 范围内,角度扩展为Δθd ,因此,gk 信道矢量表示为L 路信号的叠加,ρl 为衰落系数,阵列导向矢量b(θd) 表示为b(θd)=[1,exp(jκdesinθd),⋯,exp(j(N−1)κdesinθd)]T ,其中κ=2π/2π λλ 为波数,λ 为波长,de=λ/λ22 为中继天线间隔,θd 表示入射信号的到达角度(Angle-of-Arrival, AoA)。2.2 信号模型
整个通信过程包含两个阶段,第1阶段,
M 个NOMA用户组采用时分多址接入(Time Division Multiple Access, TDMA)技术在不同时隙发送信号给无人机。在第m 个NOMA组中,用户信道增益条件按照大小排序(‖hm,1‖<‖hm,2‖<⋯<‖hm,Kg‖ ),第m 组第k 个用户,即第(m,k) 个用户Dm,k ,以最大或受限的发射功率αm,kPu 发送信号xm,k(t) 给无人机。无人机接收信号表示为xm(t)=Kg∑k=1hm,k√αm,kPuxm,k(t)+nR(t) (5) 其中,
hm,k 为用户-无人机链路信道增益矢量,αm,k 为第(m,k) 用户发射功率系数,Pu 为用户发射功率。多天线无人机对每一个NOMA组执行接收波束成形[6,10],wm∈CN×1 为第m 组接收波束成形权矢量。根据上行NOMA基本原理,无人机使用SIC技术,依次解码信道条件好和差的用户信号。用户k 受到组内其他用户的干扰和非完美SIC产生的干扰。x′m(t) 为x′m(t) 采用SIC技术的信号。因此,无人机处第(m,k) 个用户的信号表示为[6,8]yRm,k(t)=wHmx′m(t)=wHmhm,k√αm,kPuxm,k(t)⏟第(m,k)个用户信号+wHmk−1∑i=1hm,i√αm,iPuxm,i(t)⏟组内用户干扰+wHmKg∑j=k+1hm,j√ξm,jαm,jPuxm,j(t)⏟非完美SIC产生的干扰+wHmnR(t) (6) 其中,
ξm,j 为非完美SIC系数,nR(t)~ CN(0N,σ2rIN) 为无人机处的高斯噪声,Kg 为第m 组总用户数。第2阶段,无人机将接收到的RF信号经放大转发(AF),然后采用副载波强度调制技术转为光信号[6],再发送至卫星。卫星接收到第
(m,k) 用户的信号为ym,k(t)=Gkη√PrhFSO,kyRm,k(t)+ns(t) (7) 其中,
hFSO,k 为无人机-卫星链路信道增益因子,Pr 为无人机的发射功率,η 为电光转换系数,ns(t)~ CN(0,σ2s) 为卫星处噪声,放大增益因子Gk [3]表示为Gk=√1/1E[∑Kgi=1αm,iPu|wHmhm,i|2+σ2r]E[∑Kgi=1αm,iPu|wHmhm,i|2+σ2r] ,卫星接收到第(m,k) 用户信号的输出信干噪比(Signal-to-Interference-plus-Noise-Ratio, SINR)为γm,k=G2kη2αm,kPuPr|hFSO,k|2|wHmhm,k|2G2kη2Pr|hFSO,k|2k−1∑i=1αm,iPu|wHmhm,i|2+G2kη2Pr|hFSO,k|2Kg∑j=k+1ξm,jαm,jPu|wHmhm,j|2+G2kη2Pr|hFSO,k|2σ2r+σ2d=γFSO,kγ1m,kγFSO,k(γ2m,k+γ3m,k+1)+Ckm,k (8) 其中,
γFSO,k=Prσ2s|hFSO,k|2 = ˉγFSO|gFSO,k|2 ,ˉγFSO=Prℓ2FSOσ2s ,γ1m,k=αm,kPuσ2r|wHmhm,k|2=ˉγk|wHmgm,k|2 ,γ2m,k=k−1∑i=1αm,iPu|wHmhm,i|2σ2r=k−1∑i=1ˉγi|wHmgm,i|2 ,γ3m,k=Kg∑j=k+1ξm,jαm,jPu|wHmhm,j|2σ2r=Kg∑j=k+1ˉγj|wHmgm,j|2 ,信噪比ˉγi 表示为ˉγi={αm,iˉγrk,i≤kξm,iαm,iˉγrk,i>k (9) 其中,
ˉγrk=Puℓ2kσ2r 为用户-无人机链路平均信噪比,系数Ckm,k=1η2E[γ1m,k+γ2m,k+γ3m,k+1] 。本系统的主要流程见图2系统框图。2.3 波束成形
本节通过设计波束成形方案实现系统和速率最大化为目标的优化问题。由式(8)可得到第
(m,k) 个用户可达速率为Rm,k=1/122E[log2(1+γm,k)] ,因此,优化问题可表示为maxw1,w2,⋯,wM Rsum≜12MM∑m=1Kg∑k=1E[log2(1+γm,k)] s.t. ‖wm‖2F=1,m=1,2,⋯,M (10) 每一组用户的接收波束权向量
wm(m = 1,2,⋯,M) 都是相互独立的,本文将优化问题式(10)分解成M 个子问题。第m 组用户和速率最大化可表示为maxwm Rm≜Kg∑k=112E[log2(1+γm,k)]s.t. ‖wm‖2F=1 (11) 与文献[17]类似,利用Jensen不等式,可得第
(m,k) 用户可达速率上界为12E[log2(1+γm,k)]≤12log2(1+E[γm,k]) (12) 由式(8)可得
E[γm,k]a≈E[γFSO,k]E[γ1m,k]E[γFSO,k]E[γ2m,k+γ3m,k+1]+1η2E[γ1m,k]+1η2E[γ2m,k+γ3m,k+1]=E[γFSO,k](E[γ1m,k]E[γ2m,k+γ3m,k+1])E[γFSO,k]+1η2E[γ1m,k]E[γ2m,k+γ3m,k+1]+1η2 (13) 其中,步骤a利用了Mullen不等式[18],由式(13)可知,
E[γm,k] 随着E[γFSO,k] 和E[γ1m,k]E[γ2m,k+γ3m,k+1] 单调递增,比例系数与波束权向量wm 有关,可表示为E[γ1m,k]E[γ2m,k+γ3m,k+1]=ˉγkwHmE[gm,kgHm,k]wmwHm(Kg∑i=1,i≠kˉγikE[gm,igHm,i]+IM)wm=ˉγkwHmRm,kwmwHm(Kg∑i=1,i≠kˉγikRm,i+IM)wm (14) 由式(11)和式(13)可知,难以直接求解第
m 组用户和速率最大化时的最优权向量wm 。因此,采用次优化的波束成形方案,波束权向量使得NOMA用户组中信道增益条件最好的用户达到最大可达速率。根据式(14),利用广义瑞利商理论[19],可得次优方案权向量wm 的表达式为wm=umax[(Kg∑i=1,i≠kˉγikRm,i+IM)−1Rm,k] (15) 其中,
umax{A} 表示矩阵A 的主特征向量,Rm,k=E[gm,kgHm,k] 为gm,k 的相关矩阵(Channel Correlation Matrices, CCMs),由式(4)可得其计算公式为Rm,k=E[gm,kgHm,k] = 1L∑Ll=1E[|ρl|2]b(θd)bH(θd) 。因此,Rm,k 可由信号AOA信息θd 获得。对于实际通信网络,由于信道状态信息变化快且存在反馈延迟,无人机很难获得用户的完美CSI。本文提出了基于统计CSI的波束成形方案,该统计CSI可从角度信息AOA中获得,并且更容易实现,从而避免了大量计算gm,kgHm,k 的平均值。波束成形次优化方案流程图见图3。3. 遍历和速率性能分析
基于NOMA技术和所提波束成形方案,本小节将分析系统的遍历和速率(Ergodic Sum Rate, ESR)性能。根据式(8)中的信噪比
γm,k 表达式,可得第(m,k) 用户的可达速率为Rm,k=12E[log2(1+γm,k)]=12ln2E[ln(1 + γFSO,kγ1m,kγFSO,k(γ2m,k+γ3m,k+1)+Ckm,k)]=12ln2E[ln(γFSO,kCkm,k(γ1m,k + γ2m,k+γ3m,k+1)+1)]⏟Rγ1m,k−12ln2E[ln(γFSO,kCkm,k(γ2m,k+γ3m,k+1)+1)]⏟Rγ2m,k (16) 定义
U=γFSO,kCkm,k ,式(16)中,Rγ1m,k 和Rγ2m,k 分别表示为Rγ1m,k=12ln2E[ln(Uγa,1m,k+1)]=12ln2E[ln(Uγ1m,k+1)] (17) 和
Rγ2m,k = 12ln2E[ln(Uγb,1m,k+1)] = 12ln2E[ln(Uγ2m,k+1)] (18) 其中,
γa,1m,k=γam,k+1 ,γam,k=γ1m,k + γ2m,k+γ3m,k=∑Kgi=1ˉγk|wHkhm,i|2 ,γb,1m,k=γbm,k+1=γ2m,k+γ3m,k+1 。式(17)中,
U ,γ1m,k 都是变量,Rγ1m,k 不能直接计算得到。因此,可以先计算条件U=u 下的Rγ1m,k ,再求U 的期望,计算表达式为Rγ1m,k=∞∫0Rγ1m,k|ufU(u)du (19) 其中,
Rγ1m,k|u 表示为Rγ1m,k|u=12ln2E[ln(1+γ1m,k|u)]=12ln2∞∫0e−ss(1−Mγ1m,k|u(s|u))ds (20) 其中,
Mγ1m,k|u(s|u) 由MGF函数(Moment Generating Function, MGF)计算可得Mγ1m,k|u(s|u)=Mγa,1m,k|u(us|u)=1−us∞∫0e−usx(1−Fγa,1m,k(x))dx (21) 其中,
Fγa,1m,k(x) 为γa,1m,k 的累积分布函数(Cumulative Distribution Function, CDF)。将式(20)、式(21)代入式(19)可得
Rγ1m,k Rγ1m,k=12ln2ufU(u)⋅∞∫0e−s∞∫0e−usx(1−Fγa,1m,k(x))dx⏟I1ds⏟I2du (22) 其中,先推导
Fγa,1m,k(x) ,再接着分别推导表达式I1 ,I2 和fU(u) 。Fγa,1m,k(x) 由信道输出信噪比γm,k 的PDF,根据MGF方法获得。基于Kronecker模型,式(4)中信道增益矢量
gk 也可表示为gm,k=Φ1/122m,k˜gm,k (23) 其中,
Φm,k=E[gHm,kgm,k] 为信道相关矩阵,˜gm,k=[˜gk1,˜gk2,⋯,˜gkN]T 满足˜gki~CN(0,1) 。用户-无人机链路的信道矢量
gm,k 服从相关瑞利分布,输出信噪比γm,k 表示为γm,k=ˉγk|wHmgm,k|2=ˉγkwHm(Φ1/122m,k˜gm,k)(Φ1/122m,k˜gm,k)Hwm=ˉγk˜gHm,kΦH/H22m,kwmwHmΦ1/122m,k⏟QX˜gm,k (24) 其中,
ˉγk 为信道平均信噪比,QX 为非负定矩阵。根据谱定理[16],QX 可表示为QX=[ui,1,ui,2,⋯,ui,N]diag(λi,1,λi,2,⋯,λi,N)⋅[ui,1,ui,2,⋯,ui,N]H (25) 其中,
λi,pi(pi=1,2,⋯,N) 为按照递减顺序排列的特征值,ui,pi 为对应的特征向量。因此,γm,k 可表示为γm,k=N∑pi=1λi,pi|˜hm,kui,pi|2 (26) 再利用Laplace变换,可得
γm,k 的PDF表达式为[16]fγm,k(x)=ti∑pi=1vpi∑qi=1cpi,qiΓ(qi)(λi,piˉγm,k)qixqi−1⋅exp(−xλi,piˉγm,k) (27) 其中,
t 为不同的非零特征值数目,vpi 为λi,pi 重复次数,满足∑tpi=1vpi 等于非零特征值数量。系数cpi,qi 为cpi,qi=1(vpi−qi)!λvpi−qii,pi∂vpi−qi∂svpi−qi⋅[ti∏n=1,n≠pi11+sλi,n]|s=−λ−1i,pi 由式(27)可得
γm,k 的MGF为Mγm,k(s)=∞∫0e−sxfγm,k(x)dx=ti∑pi=1vpi∑qi=1cpi,qiΓ(qi)(λi,piˉγm,k)qi∞∫0xqi−1⋅exp[−(s+1λi,piˉγm,k)x]dx=ti∑pi=1vpi∑qi=1cpi,qi(λi,piˉγm,k)qi(s+1λi,piˉγm,k)−qi (28) NOMA用户组每组包含
Kg 个用户,输出信噪比为γam,k=∑Kgi=1γm,i 。因此,γam,k 的MGF可以表示为Mγam,k(s)=Kg∏i=1Mγm,i(s)=t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg=1vpKg∑qKg=1(Kg∏i=1cpi,qi(λi,piˉγm,i)qi)(Kg∏i=1(s+1λi,piˉγm,i)−qi)⏟Ψ (29) 其中,
Ψ 的反Laplace变换表示为∑Kgi=1∑qic=1Ti,c(−1λi,piˉγm,i)xqi−ce−xλi,piˉγm,i ,Ti,c(s)=1Γ(qi−c+1)Γ(c)dc−1dsc−1(1/1Kg∏t=1,t≠i(s+1λt,ptˉγm,t)qtKg∏t=1,t≠i(s+1λt,ptˉγm,t)qt) 因此,
γam,k 的PDF表达式为fγam,k(x)=t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg=1vpKg∑qKg=1(Kg∏i=1cpi,qi(λi,piˉγm,i)qi)⋅Kg∑i=1qi∑c=1Ti,c(−1λi,piˉγm,i)xqi−ce−xλi,piˉγm,i (30) 更进一步求得
γam,k 的CDF为Fγam,k(u)=1−∞∫ufγa,1m,k(x)dx=1−t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg=1vpKg∑qKg=1(Kg∏i=1cpi,qi(λi,piˉγm,i)qi)⋅Kg∑i=1qi∑c=1Ti,c(−1λi,piˉγm,i)⋅exp(−uλi,piˉγm,i)⋅qi−c∑d=0(qi−c)!(λi,piˉγm,i)qi−c−d+1d!ud (31) 因式(17)中
γa,1m,k=γam,k+1 ,最终γa,1m,k 的CDF表达式Fγa,1m,k(x) 为Fγa,1m,k(x)=1−t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg=1vpKg∑qKg=1(Kg∏i=1cpi,qi(λi,piˉγm,i)qi)⋅Kg∑i=1qi∑c=1Ti,c(−1λi,piˉγm,i)⋅exp(1λi,piˉγm,i)⋅qi−c∑d=0(qi−c)!(λi,piˉγm,i)qi−c−d+1d!⋅d∑w=0(dw)(−1)d−wxwexp(−xλi,piˉγm,i) (32) 式(22)中,根据式(32),接着推导积分表达式
I1 为I1=∞∫0e−usx(1−Fγa,1m,k(x))dx=t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg=1vpKg∑qKg=1⋅(Kg∏i=1cpi,qi(λi,piˉγi)qi)Kg∑i=1qi∑c=1Ti,c(−1λi,piˉγi)⋅exp(1λi,piˉγi)qi−c∑d=0(qi−c)!(λi,piˉγi)qi−c−d+1d!⋅d∑w=0(dw)(−1)d−ww!(us+1λi,piˉγi)−w−1 (33) 将式(33)代入式(22)中,并根据公式
(1+x)−N−1=1Γ(N+1)G1,11,1[x|−N0] (34) 和文献[11]中的式(7.813),最终可得积分
I2 的表达式为I2=t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg=1vpKg∑qKg=1⋅(Kg∏i=1cpi,qi(λi,piˉγi)qi)Kg∑i=1qi∑c=1Ti,c(−1λi,piˉγi)⋅exp(1λi,piˉγi)⋅ qi−c∑d=0(qi−c)!(λi,piˉγi)qi−c−d+1d!⋅d∑w=0(dw)(−1)d−ww!(λi,piˉγi)w+1Γ(w+1)⋅G1,22,1[λi,piˉγiu|0,−w0] (35) 为求解式(22)中
Rγ1m,k ,还需求解U 的PDFfU(u) 。因U=γFSO,kCkm,k=ˉγFSO|gFSO,k|2Ckm,k ,且已知gFSO,k 的PDF表达式为式(2),根据概率论知识,计算fU(u) 得fU(u)=(√Ckm,kˉγFSOu)′fgFSO,k(√Ckm,kˉγFSOu)=√Ckm,kˉγFSOαβτ2u−122A0Γ(α)Γ(β)⋅G3,01,3[αβA0√Ckm,kˉγFSOu|τ2τ2−1,α−1,β−1] (36) 其中,
Ckm,k 由式(8)和式(17)计算可得Ckm,k=1η2E[γam,k+1]=1η2∞∫0xfγam,k(x)dx+1η2=1η2+1η2t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg = 1vpKg∑qKg = 1⋅(Kg∏i=1cpi,qi(λi,piˉγi)qi)Kg∑i=1qi∑c=1Ti,c(−1λi,piˉγi)⋅(qi−c+1)!(λi,piˉγi)qi−c+2 (37) 由式(36)、式(37)和文献[20]中的式(07.34.21.0013.01),经数学计算得到
Rγ1m,k 为Rγ1m,k=2α+β64πln2√Ckm,kˉγFSOαβτ2A0Γ(α)Γ(β)Θ1 (38) 其中
Θ1=t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tKg∑pKg=1vpKg∑qKg=1⋅(Kg∏i=1cpi,qi(λi,piˉγi)qi)Kg∑i=1qi∑c=1Ti,c(−1λi,piˉγi)⋅exp(1λi,piˉγi)⋅qi−c∑d=0(qi−c)!(λi,piˉγi)qi−c−d−12d!⋅d∑w=0(dw)(−1)d−ww!(λi,piˉγi)w+1Γ(w+1)⋅G8,13,8[Ckm,k(αβ)216A20ˉγFSOλi,piˉγi|−12,τ22,τ2+12ℑ] 其中,
ℑ≜τ2−12,τ22,α−12,α2,β−12,β2,−12,w−12 。同理,可得式(18)中
Rγ2m,k 表达式为Rγ2m,k=2α+β64πln2√Ckm,kˉγFSOαβτ2A0Γ(α)Γ(β)Θ2 (39) 其中
Θ2=t1∑p1=1vp1∑q1=1t2∑p2 = 1vp2∑q2 = 1⋯tk−1∑pk−1=1vpk−1∑qk−1tk+1∑pk+1=1vpk+1∑qk+1⋯tKg∑pKg = 1vpKg∑qKg = 1(Kg∏i=1,i≠kcpi,qi(λi,piˉγi)qi)⋅Kg∑i=1,i≠kqi∑c=1Ti,c(−1λi,piˉγi)exp(1λi,piˉγi)×qi−c∑d=0(qi−c)!(λi,piˉγi)qi−c−d−12d!⋅d∑w=0(dw)(−1)d−ww!(λi,piˉγi)w+1Γ(w+1)⋅G8,13,8[Ckm,k(αβ)216A20ˉγFSOλi,piˉγi|−12,τ22,τ2+12ℑ] 进一步地,将式(38)、式(39)代入式(16),得第
(m,k) 用户的可达速率Rm,k 为Rm,k=2α+β64πln2√Ckm,kˉγFSOαβτ2A0Γ(α)Γ(β)(Θ1−Θ2) (40) 最终,可得系统遍历和速率表达式为
Rsum=1MM∑m=1Kg∑k=1Rm,k (41) 4. 仿真结果与分析
本节通过计算机仿真实验验证了系统遍历和速率性能分析的正确性与所提出的方案的优越性。仿真实验中,卫星-无人机和无人机-地面用户链路分别服从伽马-伽马衰落和相关瑞利衰落。无人机天线数
N={8,32} ,所有用户K=6 分布在无人机的覆盖范围内,且被分成M=2 组。此外,假设FSO链路和RF链路的平均信噪比相同ˉγFSO=ˉγrk=ˉγ ,噪声功率σ2r=σ2d=KbTB ,Kb = 1.38×10−23 J/K ,T = 300 K ,噪声带宽B=20 MHz ,其他系统参数如表1所示[6],蒙特卡罗仿真次数为106 。为了表明本文所提方案的优势,仿真实验中分别与传统OMA策略和其他波束成形方案作了对比:表 1 系统参数参数 数值 参数 数值 卫星轨道 GEO FSO链路路径损耗(dB) 72 卫星高度(×104 km) 3.6 指向损失A0 1 无人机与用户距离dk(km) 1.1 FSO接收器接受孔径Dr(m) 0.3[13] 半径r0(m) 500 电光转换系数η 1 RF链路载波频率fc(GHz) 2 等效波束半径与指向误差偏移量之间的比率τ 6 RF链路路径损耗(dB) 114~115 与大气环境有关的大尺度单元的有效数量α 2.902 FSO链路波长(nm) 1550[13] 与大气环境有关的小尺度单元的有效数量β 2.51 光接收带宽B0(GHz) 30[12] AOA角度θd 30∘∼75∘ (1) 所提NOMA策略与OMA策略对比:将NOMA策略和OMA策略两种情况下的性能作比较。在星地融合网络中,地面用户采用正交多址接入(Orthogonal Multiple Access, OMA)[9],仿真图中表示为“OMA 策略”。本文中NOMA策略仿真了非完美SIC/完美SIC的情况,图中表示为“所提NOMA策略-非完美SIC/完美SIC”。
(2) 波束成形方案:在卫星-中继-用户的星地融合网络中,将所提出的BF方案与其他接收/发送BF方案相比,如最大比合并/发送(Maximal Ratio Combining/Transmission, MRC/ MRT)[21]和最大比合并/迫零(MRC/Zero-Forcing, MRC/ZF)[22] 两种BF方案,图中分别表示为“MRC/MRT BF方案”和“MRC/ZF BF方案”。
图4所示为NOMA和OMA不同策略下系统遍历和速率随平均信噪比
ˉγ 增大而显著提高,设天线数N={8,32} ,K=6 ,αm,j = 0.8 。图中理论推导和仿真结果一致,证明理论推导式(41)的正确性。由图可知,NOMA策略下完美SIC的系统遍历和速率性能优于OMA策略,这是因为采用NOMA技术可以同时服务多个用户,具有更好的频谱效率。同时,在NOMA和OMA策略下,可以发现天线数配置不同时,系统遍历和速率都随着天线数N 增加有所提高。因此,可以通过增加无人机的天线数提高系统性能。图5对比了NOMA完美/非完美SIC条件下采用不同波束成形方案,系统遍历和速率随
ˉγ 的变化曲线,其中N = 8 ,K=6 ,αm,j = 0.8 ,ξm,j = {0,0.04} 。由图可以看出,随着影响因子ξm,j 增大,非完美SIC遍历和速率比完美SIC差,当ˉγ 在低信噪比区域时,其性能与完美SIC相差不大,但是在高信噪比区域随ˉγ 增大与完美SIC差距越来越大,且达到饱和状态。图中还比较了同一NOMA组中采用两种波束成形方案的性能,由用户信道增益条件好坏决定,第1种为权向量wm,1 满足信道增益条件最好的用户达到最大可达速率,第2种为权向量wm,2 满足信道增益条件最差的用户达到最大可达速率,波束成形权向量如式(15)。通过对比这两种方案可以看出,完美SIC情况下系统遍历和速率几乎接近,而非完美SIC情况下,第1种方案性能高于第2种方案。因此,本文采用了第1种方案。图6对比了不同波束成形方案下系统遍历和速率随
ˉγ 的变化曲线,其中N = 8 ,K=6 ,αm,j = 0.8 。将NOMA完美SIC条件下BF方案与MRC/ZF, MRC/MRT两种BF方案作对比。由图可以发现,系统遍历和速率性能在所提BF方案下高于MRC/ZF、MRC/MRT的BF方案,且随ˉγ 的增大性能显著提高,从而表明所提波束成形方案的优越性。5. 结束语
本文针对RF/FSO混合星空地融合网络,研究了多天线波束成形技术和上行NOMA技术相结合的系统遍历与速率性能。首先,在无人机采用多天线和上行NOMA技术条件下,为实现系统遍历和速率最大化,提出了基于统计CSI的波束成形方案。接着,假设卫星-无人机链路服从伽马-伽马衰落,无人机-地面用户链路服从相关瑞利衰落,本文在考虑NOMA非完美SIC条件下,推导了系统和速率的闭合表达式。最后,数值仿真结果验证了理论分析的正确性。仿真结果表明,与OMA方案相比,所提方案提高了系统性能,并且与MRC/ZF, MRC/MRT波束成形方案相比,所提方案具有更好的性能优势。本文研究为NOMA技术在卫星通信中的应用提供了理论依据,并对实际通信系统的设计具有指导意义。此外,我们将在今后的工作中进一步地基于NOMA上下行混合链路技术的星空地融合网络作分析和研究。
-
表 1 系统参数
参数 数值 参数 数值 卫星轨道 GEO FSO链路路径损耗(dB) 72 卫星高度(×104 km) 3.6 指向损失A0 1 无人机与用户距离dk(km) 1.1 FSO接收器接受孔径Dr(m) 0.3[13] 半径r0(m) 500 电光转换系数η 1 RF链路载波频率fc(GHz) 2 等效波束半径与指向误差偏移量之间的比率τ 6 RF链路路径损耗(dB) 114~115 与大气环境有关的大尺度单元的有效数量α 2.902 FSO链路波长(nm) 1550[13] 与大气环境有关的小尺度单元的有效数量β 2.51 光接收带宽B0(GHz) 30[12] AOA角度θd 30∘∼75∘ -
[1] KODHELI O, LAGUNAS E, MATURO N, et al. Satellite communications in the new space era: A survey and future challenges[J]. IEEE Communications Surveys & Tutorials, 2021, 23(1): 70–109. doi: 10.1109/COMST.2020.3028247 [2] 徐常志, 靳一, 李立, 等. 面向6G的星地融合无线传输技术[J]. 电子与信息学报, 2021, 43(1): 28–36. doi: 10.11999/JEIT200363XU Changzhi, JIN Yi, LI Li, et al. Wireless transmission technology of satellite-terrestrial integration for 6G mobile communication[J]. Journal of Electronics &Information Technology, 2021, 43(1): 28–36. doi: 10.11999/JEIT200363 [3] HUANG Qingquan, LIN Min, WANG Junbo, et al. Energy efficient beamforming schemes for satellite-aerial-terrestrial networks[J]. IEEE Transactions on Communications, 2020, 68(6): 3863–3875. doi: 10.1109/TCOMM.2020.2978044 [4] KAMGA G N, AÏSSA S, RASETHUNTSA T R, et al. Mixed RF/FSO communications with outdated-CSI-based relay selection under double generalized Gamma turbulence, generalized pointing errors, and Nakagami-m fading[J]. IEEE Transactions on Wireless Communications, 2021, 20(5): 2761–2775. doi: 10.1109/TWC.2019.2954866 [5] XU Guanjun and SONG Zhaohui. Performance analysis for mixed κ-μ fading and M-distribution dual-hop radio frequency/free space optical communication systems[J]. IEEE Transactions on Wireless Communications, 2021, 20(3): 1517–1528. doi: 10.1109/TWC.2020.3034104 [6] HUANG Qingquan, LIN Min, ZHU Weiping, et al. Uplink massive access in mixed RF/FSO satellite-aerial-terrestrial networks[J]. IEEE Transactions on Communications, 2021, 69(4): 2413–2426. doi: 10.1109/TCOMM.2021.3049364 [7] WANG Lina, WU Yanan, ZHANG Haijun, et al. Resource allocation for NOMA based space-terrestrial satellite networks[J]. IEEE Transactions on Wireless Communications, 2021, 20(2): 1065–1075. doi: 10.1109/TWC.2020.3030704 [8] JAMALI M V and MAHDAVIFAR H. Uplink non-orthogonal multiple access over mixed RF-FSO systems[J]. IEEE Transactions on Wireless Communications, 2020, 19(5): 3558–3574. doi: 10.1109/TWC.2020.2974947 [9] ZHAO Jinlong, YUE Xinwei, KANG Shaoli, et al. Joint effects of imperfect CSI and SIC on NOMA based satellite-terrestrial systems[J]. IEEE Access, 2021, 9: 12545–12554. doi: 10.1109/ACCESS.2021.3051306 [10] LIN Zhi, LIN Min, WANG Junbo, et al. Joint beamforming and power allocation for satellite-terrestrial integrated networks with non-orthogonal multiple access[J]. IEEE Journal of Selected Topics in Signal Processing, 2019, 13(3): 657–670. doi: 10.1109/JSTSP.2019.2899731 [11] GRADSHTEYN I S and RYZHIK I M. Table of Integrals, Series, and Products[M]. 7th ed. San Diego, USA: Academic Press, 2007. [12] ANTONINI M, BETTI S, CARROZZO V, et al. Feasibility analysis of a HAP-LEO optical link for data relay purposes[C]. 2006 IEEE Aerospace Conference, Big Sky, USA, 2006: 1–7. [13] LI Mi, HONG Yifeng, ZENG Cheng, et al. Investigation on the UAV-to-satellite optical communication systems[J]. IEEE Journal on Selected Areas in Communications, 2018, 36(9): 2128–2138. doi: 10.1109/JSAC.2018.2864419 [14] TRINH P V, DANG N T, and PHAM A T. All-optical relaying FSO systems using EDFA combined with optical hard-limiter over atmospheric turbulence channels[J]. Journal of Lightwave Technology, 2015, 33(19): 4132–4144. doi: 10.1109/JLT.2015.2466432 [15] ERTEL R B, CARDIERI P, SOWERBY K W, et al. Overview of spatial channel models for antenna array communication systems[J]. IEEE Personal Communications, 1998, 5(1): 10–22. doi: 10.1109/98.656151 [16] LIN Min, YANG Lvxi, ZHU Weiping, et al. An open-loop adaptive space-time transmit scheme for correlated fading channels[J]. IEEE Journal of Selected Topics in Signal Processing, 2008, 2(2): 147–158. doi: 10.1109/JSTSP.2008.922482 [17] YUAN Jide, MATTHAIOU M, JIN Shi, et al. Tightness of Jensen’s bounds and applications to MIMO communications[J]. IEEE Transactions on Communications, 2017, 65(2): 579–593. doi: 10.1109/TCOMM.2016.2623945 [18] MULLEN K. The Teacher's corner: A note on the ratio of two independent random variables[J]. The American Statistician, 1967, 21(3): 30–31. doi: 10.1080/00031305.1967.10479818 [19] GOLUB G H and VAN LOAN C F. Matrix Computations[M]. 2nd ed. Baltimore: The Johns Hopkins University Press, 1989. [20] WOLFRAM I. Mathematica edition: Version 12.0 [OL]. https://www.wolfram.com, 2020. [21] ARTI M K and BHATNAGAR M R. Beamforming and combining in hybrid satellite-terrestrial cooperative systems[J]. IEEE Communications Letters, 2014, 18(3): 483–486. doi: 10.1109/LCOMM.2014.012214.132738 [22] HUANG Qingquan, LIN Min, AN Kang, et al. Secrecy performance of hybrid satellite-terrestrial relay networks in the presence of multiple eavesdroppers[J]. IET Communications, 2018, 12(1): 26–34. doi: 10.1049/iet-com.2017.0948 期刊类型引用(1)
1. 王雪,孟姝宇,钱志鸿. 面向6G全域融合的智能接入关键技术综述. 电子与信息学报. 2024(05): 1613-1631 . 本站查看
其他类型引用(4)
-