Processing math: 100%
高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于上行非正交多址接入技术的星空地融合网络性能分析

袁祖霞 程铭 郭克锋

袁祖霞, 程铭, 郭克锋. 基于上行非正交多址接入技术的星空地融合网络性能分析[J]. 电子与信息学报, 2022, 44(8): 2666-2676. doi: 10.11999/JEIT220379
引用本文: 袁祖霞, 程铭, 郭克锋. 基于上行非正交多址接入技术的星空地融合网络性能分析[J]. 电子与信息学报, 2022, 44(8): 2666-2676. doi: 10.11999/JEIT220379
YUAN Zuxia, CHENG Ming, GUO Kefeng. Performance Analysis of Satellite-Aerial-Terrestrial Integrated Network Based on Uplink NOMA Technology[J]. Journal of Electronics & Information Technology, 2022, 44(8): 2666-2676. doi: 10.11999/JEIT220379
Citation: YUAN Zuxia, CHENG Ming, GUO Kefeng. Performance Analysis of Satellite-Aerial-Terrestrial Integrated Network Based on Uplink NOMA Technology[J]. Journal of Electronics & Information Technology, 2022, 44(8): 2666-2676. doi: 10.11999/JEIT220379

基于上行非正交多址接入技术的星空地融合网络性能分析

doi: 10.11999/JEIT220379
基金项目: 国家自然科学基金(62001517),江苏省研究生科研与实践创新计划项目(KYCX19_0896)
详细信息
    作者简介:

    袁祖霞:女,1983年生,博士生,研究方向为卫星通信、通信信号处理等

    程铭:男,1991年生,讲师,研究方向为毫米波通信、星天地融合网络、预编码技术等

    郭克锋:男,1990年生,讲师,研究方向为MIMO通信系统、卫星通信、多用户通信系统等

    通讯作者:

    袁祖霞 yuanzuxia@126.com

  • 中图分类号: TN929.1

Performance Analysis of Satellite-Aerial-Terrestrial Integrated Network Based on Uplink NOMA Technology

Funds: The National Natural Science Foundation of China (62001517), The Postgraduate Research Practice Innovation Program of Jiangsu Province (KYCX19_0896)
  • 摘要: 针对光电混合的星空地融合网络上行链路,该文研究了多天线波束成形技术和上行非正交多址接入(NOMA)技术相结合的系统遍历和速率性能。首先,在无人机采用多天线和上行NOMA技术条件下,为实现系统和速率最大化,提出了一种基于统计信道状态信息的波束成形方案。接着,假设卫星-无人机链路采用自由空间光链路且服从伽马-伽马衰落,无人机-地面用户链路采用射频链路且服从相关瑞利衰落,推导了系统和速率的闭合表达式。最后,通过数值仿真验证了理论分析的正确性。仿真结果表明,与正交多址接入(OMA)方案相比,所提方案提高了系统性能,并且与基准波束成形(BF)方案相比,所提方案具有更好的性能优势。
  • 第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))。

    图1所示,本文研究了基于上行NOMA技术的光电混合SATIN,其中,K个地面用户(D)与卫星(S)之间进行通信,由于遮挡效应,忽略卫星和地面用户之间的直达链路影响[6],因此利用无人机作为中继(R)协作地面用户与卫星进行通信,且无人机具有静态悬停特性。假设卫星和地面用户为单天线,无人机配置N根天线。所有地面用户分布于无人机的覆盖范围内,且信道相关性强和信道增益差异性大的用户被分为一组[10]K个用户被分配到M个NOMA组中,每组Kg个用户且采用NOMA技术接入。

    图 1  系统模型
    2.1.1   FSO链路

    考虑到频谱资源的短缺和用户容量的需求,假设卫星-无人机链路用FSO进行连接,因路径损耗、信道衰落和天线增益的影响,FSO链路信道表示为[12]

    hFSO,k=FSOgFSO,k (1)

    其中,FSO=GFt+GFrAtenFSAtenAtmLLensMe,其中,GFtGFr分别为发送和接收增益,AtenFSAtenAtmLLensMe分别为自由空间损耗、大气衰减、透镜损耗、系统边缘。因无人机-卫星的FSO链路传输受多普勒效应、大气湍流和指向误差等重要因素影响,文献[13]对FSO传输特性作了相关研究,但模型建立较复杂。根据文献[14]中FSO链路的建模,设信道增益gFSO,k服从伽马-伽马分布,该分布描述了弱到强湍流的闪烁统计信息,包含了大尺度和小尺度湍流涡旋未考虑指向误差[6],其概率密度函数(Probability Density Function, PDF)表达式为

    f|gFSO|(x)=αβτ2A0Γ(α)Γ(β)G3,01,3[αβA0x|τ2τ21,α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],其中,GRtGRr分别表示发送和接收增益,cfcdkαh分别表示光速、载波频率、无人机和地面用户之间的距离以及路径损耗因子。根据文献[16]的信道散射模型,信道增益矢量gk表示为

    gk=1LLl=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(N1)κdesinθd)]T,其中κ=2π/2π λλ为波数,λ为波长,de=λ/λ22为中继天线间隔,θd表示入射信号的到达角度(Angle-of-Arrival, AoA)。

    整个通信过程包含两个阶段,第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)=Kgk=1hm,kαm,kPuxm,k(t)+nR(t) (5)

    其中,hm,k为用户-无人机链路信道增益矢量,αm,k为第(m,k)用户发射功率系数,Pu为用户发射功率。多天线无人机对每一个NOMA组执行接收波束成形[6,10]wmCN×1为第m组接收波束成形权矢量。根据上行NOMA基本原理,无人机使用SIC技术,依次解码信道条件好和差的用户信号。用户k受到组内其他用户的干扰和非完美SIC产生的干扰。xm(t)xm(t)采用SIC技术的信号。因此,无人机处第(m,k)个用户的信号表示为[6,8]

    yRm,k(t)=wHmxm(t)=wHmhm,kαm,kPuxm,k(t)(m,k)+wHmk1i=1hm,iαm,iPuxm,i(t)+wHmKgj=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|2k1i=1αm,iPu|wHmhm,i|2+G2kη2Pr|hFSO,k|2Kgj=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=Pr2FSOσ2s, γ1m,k=αm,kPuσ2r|wHmhm,k|2=ˉγk|wHmgm,k|2, γ2m,k=k1i=1αm,iPu|wHmhm,i|2σ2r=k1i=1ˉγi|wHmgm,i|2, γ3m,k=Kgj=k+1ξm,jαm,jPu|wHmhm,j|2σ2r=Kgj=k+1ˉγj|wHmgm,j|2,信噪比ˉγi表示为

    ˉγi={αm,iˉγrk,ikξm,iαm,iˉγrk,i>k (9)

    其中,ˉγrk=Pu2kσ2r为用户-无人机链路平均信噪比,系数Ckm,k=1η2E[γ1m,k+γ2m,k+γ3m,k+1]。本系统的主要流程见图2系统框图。

    图 2  系统框图

    本节通过设计波束成形方案实现系统和速率最大化为目标的优化问题。由式(8)可得到第(m,k)个用户可达速率为Rm,k=1/122E[log2(1+γm,k)],因此,优化问题可表示为

     maxw1,w2,,wM Rsum12MMm=1Kgk=1E[log2(1+γm,k)] s.twm2F=1,m=1,2,,M (10)

    每一组用户的接收波束权向量wm(m = 1,2,,M)都是相互独立的,本文将优化问题式(10)分解成M个子问题。第m组用户和速率最大化可表示为

    maxwm RmKgk=112E[log2(1+γm,k)]s.twm2F=1 (11)

    与文献[17]类似,利用Jensen不等式,可得第(m,k)用户可达速率上界为

    12E[log2(1+γm,k)]12log2(1+E[γm,k]) (12)

    由式(8)可得

    E[γm,k]aE[γ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(Kgi=1,ikˉγikE[gm,igHm,i]+IM)wm=ˉγkwHmRm,kwmwHm(Kgi=1,ikˉγikRm,i+IM)wm (14)

    由式(11)和式(13)可知,难以直接求解第m组用户和速率最大化时的最优权向量wm。因此,采用次优化的波束成形方案,波束权向量使得NOMA用户组中信道增益条件最好的用户达到最大可达速率。根据式(14),利用广义瑞利商理论[19],可得次优方案权向量wm的表达式为

    wm=umax[(Kgi=1,ikˉγ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] = 1LLl=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,k12ln2E[ln(γFSO,kCkm,k(γ2m,k+γ3m,k+1)+1)]Rγ2m,k (16)

    定义U=γFSO,kCkm,k,式(16)中,Rγ1m,kRγ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)]=12ln20ess(1Mγ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)=1us0eusx(1Fγ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)0es0eusx(1Fγa,1m,k(x))dxI1dsI2du (22)

    其中,先推导Fγa,1m,k(x),再接着分别推导表达式I1, I2fU(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满足˜gkiCN(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,kQX˜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=Npi=1λi,pi|˜hm,kui,pi|2 (26)

    再利用Laplace变换,可得γm,k的PDF表达式为[16]

    fγm,k(x)=tipi=1vpiqi=1cpi,qiΓ(qi)(λi,piˉγm,k)qixqi1exp(xλi,piˉγm,k) (27)

    其中,t为不同的非零特征值数目,vpiλi,pi重复次数,满足tpi=1vpi等于非零特征值数量。系数cpi,qi

    cpi,qi=1(vpiqi)!λvpiqii,pivpiqisvpiqi[tin=1,npi11+sλi,n]|s=λ1i,pi

    由式(27)可得γm,k的MGF为

    Mγm,k(s)=0esxfγm,k(x)dx=tipi=1vpiqi=1cpi,qiΓ(qi)(λi,piˉγm,k)qi0xqi1exp[(s+1λi,piˉγm,k)x]dx=tipi=1vpiqi=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)=Kgi=1Mγm,i(s)=t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg=1vpKgqKg=1(Kgi=1cpi,qi(λi,piˉγm,i)qi)(Kgi=1(s+1λi,piˉγm,i)qi)Ψ (29)

    其中,Ψ的反Laplace变换表示为Kgi=1qic=1Ti,c(1λi,piˉγm,i)xqicexλi,piˉγm,iTi,c(s)=1Γ(qic+1)Γ(c)dc1dsc1(1/1Kgt=1,ti(s+1λt,ptˉγm,t)qtKgt=1,ti(s+1λt,ptˉγm,t)qt)

    因此,γam,k的PDF表达式为

    fγam,k(x)=t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg=1vpKgqKg=1(Kgi=1cpi,qi(λi,piˉγm,i)qi)Kgi=1qic=1Ti,c(1λi,piˉγm,i)xqicexλi,piˉγm,i (30)

    更进一步求得γam,k的CDF为

    Fγam,k(u)=1ufγa,1m,k(x)dx=1t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg=1vpKgqKg=1(Kgi=1cpi,qi(λi,piˉγm,i)qi)Kgi=1qic=1Ti,c(1λi,piˉγm,i)exp(uλi,piˉγm,i)qicd=0(qic)!(λi,piˉγm,i)qicd+1d!ud (31)

    因式(17)中γa,1m,k=γam,k+1,最终γa,1m,k的CDF表达式Fγa,1m,k(x)

    Fγa,1m,k(x)=1t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg=1vpKgqKg=1(Kgi=1cpi,qi(λi,piˉγm,i)qi)Kgi=1qic=1Ti,c(1λi,piˉγm,i)exp(1λi,piˉγm,i)qicd=0(qic)!(λi,piˉγm,i)qicd+1d!dw=0(dw)(1)dwxwexp(xλi,piˉγm,i) (32)

    式(22)中,根据式(32),接着推导积分表达式I1

    I1=0eusx(1Fγa,1m,k(x))dx=t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg=1vpKgqKg=1(Kgi=1cpi,qi(λi,piˉγi)qi)Kgi=1qic=1Ti,c(1λi,piˉγi)exp(1λi,piˉγi)qicd=0(qic)!(λi,piˉγi)qicd+1d!dw=0(dw)(1)dww!(us+1λi,piˉγi)w1 (33)

    将式(33)代入式(22)中,并根据公式

    (1+x)N1=1Γ(N+1)G1,11,1[x|N0] (34)

    和文献[11]中的式(7.813),最终可得积分I2的表达式为

    I2=t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg=1vpKgqKg=1(Kgi=1cpi,qi(λi,piˉγi)qi)Kgi=1qic=1Ti,c(1λi,piˉγi)exp(1λi,piˉγi) qicd=0(qic)!(λi,piˉγi)qicd+1d!dw=0(dw)(1)dww!(λi,piˉγi)w+1Γ(w+1)G1,22,1[λi,piˉγiu|0,w0] (35)

    为求解式(22)中Rγ1m,k,还需求解U的PDF fU(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αβτ2u122A0Γ(α)Γ(β)G3,01,3[αβA0Ckm,kˉγFSOu|τ2τ21,α1,β1] (36)

    其中,Ckm,k由式(8)和式(17)计算可得

    Ckm,k=1η2E[γam,k+1]=1η20xfγam,k(x)dx+1η2=1η2+1η2t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg = 1vpKgqKg = 1(Kgi=1cpi,qi(λi,piˉγi)qi)Kgi=1qic=1Ti,c(1λi,piˉγi)(qic+1)!(λi,piˉγi)qic+2 (37)

    由式(36)、式(37)和文献[20]中的式(07.34.21.0013.01),经数学计算得到Rγ1m,k

    Rγ1m,k=2α+β64πln2Ckm,kˉγFSOαβτ2A0Γ(α)Γ(β)Θ1 (38)

    其中

    Θ1=t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tKgpKg=1vpKgqKg=1(Kgi=1cpi,qi(λi,piˉγi)qi)Kgi=1qic=1Ti,c(1λi,piˉγi)exp(1λi,piˉγi)qicd=0(qic)!(λi,piˉγi)qicd12d!dw=0(dw)(1)dww!(λi,piˉγi)w+1Γ(w+1)G8,13,8[Ckm,k(αβ)216A20ˉγFSOλi,piˉγi|12,τ22,τ2+12]

    其中,τ212,τ22,α12,α2,β12,β2,12,w12

    同理,可得式(18)中Rγ2m,k表达式为

    Rγ2m,k=2α+β64πln2Ckm,kˉγFSOαβτ2A0Γ(α)Γ(β)Θ2 (39)

    其中

    Θ2=t1p1=1vp1q1=1t2p2 = 1vp2q2 = 1tk1pk1=1vpk1qk1tk+1pk+1=1vpk+1qk+1tKgpKg = 1vpKgqKg = 1(Kgi=1,ikcpi,qi(λi,piˉγi)qi)Kgi=1,ikqic=1Ti,c(1λi,piˉγi)exp(1λi,piˉγi)×qicd=0(qic)!(λi,piˉγi)qicd12d!dw=0(dw)(1)dww!(λ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πln2Ckm,kˉγFSOαβτ2A0Γ(α)Γ(β)(Θ1Θ2) (40)

    最终,可得系统遍历和速率表达式为

    Rsum=1MMm=1Kgk=1Rm,k (41)

    本节通过计算机仿真实验验证了系统遍历和速率性能分析的正确性与所提出的方案的优越性。仿真实验中,卫星-无人机和无人机-地面用户链路分别服从伽马-伽马衰落和相关瑞利衰落。无人机天线数N={8,32},所有用户K=6分布在无人机的覆盖范围内,且被分成M=2组。此外,假设FSO链路和RF链路的平均信噪比相同ˉγFSO=ˉγrk=ˉγ,噪声功率σ2r=σ2d=KbTBKb = 1.38×1023 J/KT = 300 K,噪声带宽B=20 MHz,其他系统参数如表1所示[6],蒙特卡罗仿真次数为106。为了表明本文所提方案的优势,仿真实验中分别与传统OMA策略和其他波束成形方案作了对比:

    表 1  系统参数
    参数数值参数数值
    卫星轨道GEOFSO链路路径损耗(dB)72
    卫星高度(×104 km)3.6指向损失A01
    无人机与用户距离dk(km)1.1FSO接收器接受孔径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角度θd3075
    下载: 导出CSV 
    | 显示表格

    (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增加有所提高。因此,可以通过增加无人机的天线数提高系统性能。

    图 4  NOMA/OMA策略下遍历和速率曲线

    图5对比了NOMA完美/非完美SIC条件下采用不同波束成形方案,系统遍历和速率随ˉγ的变化曲线,其中N = 8K=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种方案。

    图 5  NOMA完美/非完美SIC条件下采用不同BF方案遍历和速率曲线

    图6对比了不同波束成形方案下系统遍历和速率随ˉγ的变化曲线,其中N = 8K=6αm,j = 0.8。将NOMA完美SIC条件下BF方案与MRC/ZF, MRC/MRT两种BF方案作对比。由图可以发现,系统遍历和速率性能在所提BF方案下高于MRC/ZF、MRC/MRT的BF方案,且随ˉγ的增大性能显著提高,从而表明所提波束成形方案的优越性。

    图 6  不同波束成形方案下遍历和速率曲线

    本文针对RF/FSO混合星空地融合网络,研究了多天线波束成形技术和上行NOMA技术相结合的系统遍历与速率性能。首先,在无人机采用多天线和上行NOMA技术条件下,为实现系统遍历和速率最大化,提出了基于统计CSI的波束成形方案。接着,假设卫星-无人机链路服从伽马-伽马衰落,无人机-地面用户链路服从相关瑞利衰落,本文在考虑NOMA非完美SIC条件下,推导了系统和速率的闭合表达式。最后,数值仿真结果验证了理论分析的正确性。仿真结果表明,与OMA方案相比,所提方案提高了系统性能,并且与MRC/ZF, MRC/MRT波束成形方案相比,所提方案具有更好的性能优势。本文研究为NOMA技术在卫星通信中的应用提供了理论依据,并对实际通信系统的设计具有指导意义。此外,我们将在今后的工作中进一步地基于NOMA上下行混合链路技术的星空地融合网络作分析和研究。

  • 图  1  系统模型

    图  2  系统框图

    图  3  波束成形次优化方案的流程图

    图  4  NOMA/OMA策略下遍历和速率曲线

    图  5  NOMA完美/非完美SIC条件下采用不同BF方案遍历和速率曲线

    图  6  不同波束成形方案下遍历和速率曲线

    表  1  系统参数

    参数数值参数数值
    卫星轨道GEOFSO链路路径损耗(dB)72
    卫星高度(×104 km)3.6指向损失A01
    无人机与用户距离dk(km)1.1FSO接收器接受孔径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角度θd3075
    下载: 导出CSV
  • [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/JEIT200363

    XU 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)

  • 加载中
图(6) / 表(1)
计量
  • 文章访问数:  714
  • HTML全文浏览量:  243
  • PDF下载量:  86
  • 被引次数: 5
出版历程
  • 收稿日期:  2022-04-01
  • 修回日期:  2022-07-10
  • 录用日期:  2022-07-13
  • 网络出版日期:  2022-07-15
  • 刊出日期:  2022-08-17

目录

/

返回文章
返回