1.   引 言

卫星通信凭借其覆盖范围广、通信距离远、不受地理条件限制等众多优点，已经被广泛应用于偏远地区通信以及导航定位、抗震抢险等领域^[1-3]，并将成为下一代无线通信的关键技术之一。然而，卫星与地面用户之间存在大时延和大路径损耗以及遮蔽效应导致的视距传输受阻等问题，使得卫星系统的用户体验有时无法得到保证。在这种情况下，基于中继转发的星地协作传输技术被认为是提升卫星通信服务质量的有效手段之一^[4]。在大多研究的星地协作传输网络中，通常采用地面中继将接收到的卫星信号转发给地面用户。例如，文献[5]研究了单用户场景下地面中继采用放大转发(Amplify and Forward, AF)协议的星地协作传输网络的误码性能与中断性能。进一步，针对地面多用户的星地协作传输网络，文献[6]在采用最优用户选择方案的情况下，推导得到用户的中断概率闭合表达式。需要指出的是，虽然使用地面中继可以建立卫星与地面用户之间的高质量通信链路，但对人口稀少的偏远地区而言，建造地面中继站存在高成本、低回报等问题，因此需要探索其他更加实用的解决方案^[7]。

跟地面中继相比，无人机由于其机动性好、通信组网方式灵活等优势，作为空中中继协助卫星与地面用户通信^[8]，可以实现增强接收信号功率、提高系统容量等目的，并且得到了学术界和工业界的高度重视。例如，文献[9]研究了基于无人机中继转发的星地协作网络中多用户传输场景下系统的中断性能；文献[10]研究了无人机中继采用AF协议和多用户调度方案下的星地协作网络性能；文献[11]针对无人机辅助卫星通信系统，分析了系统的中断性能。然而，考虑到卫星服务用户数量越来越多，需要进一步提高系统资源利用率以及服务用户通信质量。现有的正交多址 (Orthogonal Multiple Access, OMA)技术已经无法满足上述需求。近年来，非正交多址(Non-Orthogonal Multiple Access, NOMA)技术以其可大大提高系统频谱资源利用率和用户公平性等独特优势，已经成为极具发展前景的新型多址技术^[12]。在这种情况下，已经有学者研究如何将NOMA技术应用于卫星通信系统。例如，文献[13,14]研究了两用户场景下基于NOMA的星地协作传输系统性能；文献[15,16]针对基于NOMA的卫星通信系统中多用户传输场景，分析了系统性能；文献[17]针对基于NOMA的无人机辅助卫星通信系统，分析了无人机中继采用AF协议下的系统中断性能。

总的来看，上述文献对星地协作传输网络进行了深入的研究，验证了中继协作技术能够显著提升卫星通信系统存在遮蔽效应下的性能，但是它们主要存在以下问题：一是大多数文献，例如文献[5]仅研究了单用户场景下星地协作传输网络的性能；二是虽然也有相关文献研究多用户场景，但是大多数文献，例如文献[6,10]都是在假设准确信道状态信息(Channel State Information, CSI)已知的情况下，研究了采用用户调度或空分多址(Space Division Multiple Access, SDMA)方案下星地协作传输网络的性能；三是将NOMA技术应用于卫星通信系统的大多数文献中，例如文献[13-16]都是采用地面中继将接收到的卫星信号转发给地面用户；四是虽然也有相关文献提出基于NOMA的无人机辅助卫星通信系统，但是大多数文献，例如文献[17]仅考虑无人机配置单天线且地面用户簇内仅包含远近两个用户的通信场景，没有充分利用空间资源和考虑卫星服务用户数量越来越多的实际情况。在这种情况下，本文针对无人机辅助的卫星通信系统下行链路，研究基于SDMA和协作NOMA相结合的多用户传输系统。首先，配置多根天线的无人机作为中继站辅助卫星通信，采用NOMA技术服务多个地面用户，并得到用户的输出信干噪比(Signal-to-Interference-plus-Noise Ratio, SINR)表达式。然后，建立以平均SINR最大化为准则的优化问题，并在无人机仅已知信道角度信息的情况下，得到无人机对卫星信号的接收波束成形权矢量。进一步，为了降低算法复杂度，提出基于角度信息的迫零波束成形方案，得到无人机对多个地面用户的发射波束成形权矢量。其次，在卫星-无人机链路服从相关阴影莱斯分布，而无人机-地面用户链路服从Nakagami-m分布的条件下，分别推导出系统的中断概率闭合表达式和近似表达式。最后，仿真结果验证了所提方案的优越性和理论分析的正确性。

2.   系统模型

如图1所示，本文研究无人机辅助的卫星通信系统下行链路，其中静止轨道卫星S通过无人机中继R转发信号，与地面用户D进行通信。假设卫星S采用点波束技术，无人机中继R配置N元的均匀直线阵(Uniform Linear Array, ULA)，地面用户D配置单天线。为了实现多用户同时通信，将无人机覆盖范围内的地面用户中信道相关性强和信道增益差异大的用户划分为一簇^[18]，从而将地面用户划分为L簇，且簇内用户采用NOMA技术提高频谱利用率。跟针对单用户场景以及多用户场景下采用用户调度或SDMA方案的星地协作传输网络的文献相比^[6,10]，本文的研究更具一般性。为了便于理解，本文将首先介绍信道模型和信号模型，然后对所提的波束成形方案进行描述。

图 1  系统模型图

下载: 全尺寸图片 幻灯片

2.1   信道模型

2.1.1   S-R信道模型

在卫星通信中，为了克服远距离传输导致的大路径损耗，通常采用高增益的点波束技术。在这种情况下，将卫星的波束增益、电波传播过程中的自由空间损耗和小尺度衰落的影响加以考虑之后，卫星到无人机链路，即S-R链路信道矢量可以表示为

其中，${\zeta _{{{\rm{SR}}} }}$为信道系数，计算公式为${\zeta _{{\rm{SR}}}} = {\lambda \sqrt {{G_{\rm{S}}}} }/ {\left( {4\pi {d_{\rm{R}}}} \right)}$，其中$\lambda$为波长，${d_{\rm{R}}}$为卫星到无人机的距离，${G_{{\rm{S}}} }$为卫星波束增益，通常建模为^[19]

其中，$G_{{\rm{S}}} ^{{\rm{max}}}$为卫星波束的最大增益，${J}_{1}(\cdot)$和${J}_{3}(\cdot)$分别表示1阶和3阶的第1类贝塞尔函数，且${u_{\rm{s}}} = 2.07123 {{\sin {\theta _{{\rm{S}}} }} \mathord{\left/ {\vphantom {{2.07123\sin {\theta _{{\rm{S}}} }} {\sin {\theta _{{S} ,3dB}}}}} \right. } {\sin {\theta _{{{\rm{S}}} ,3{\rm{\,dB}}}}}}$，${\theta _{{\rm{S}}} }$是无人机偏离卫星波束最大增益方向的角度，${\theta _{{{\rm{S}}} ,3\,{\text{dB}}}}$是卫星波束的3 dB波束角度。

此外，式(1)中${{\boldsymbol{g}}_{{\rm{SR}}}}$表示S-R链路的小尺度衰落，通常服从相关阴影莱斯分布^[19]，表示为

其中，${{\boldsymbol{\bar g}}_{{\rm{SR}}}}$为直达径(Line-of-Sight, LoS)分量，表示为

其中，${Z_{{{\rm{SR}}} }}$是一个服从Nakagami-m分布的随机变量，衰落参数为${m_{{\rm{R}}} }$，平均功率为${\varOmega _{\rm{R}}}$。此外，由于无人机配置了ULA，阵列导向矢量${\boldsymbol{a}}\left( {{\theta _{\rm{r}}}} \right)$由式(5)给出

其中，${\theta _{{\rm{r}}} }$为卫星信号的到达角，${d_e}$表示ULA的阵元间隔。此外，式(3)中${\boldsymbol{\tilde \varPhi }}_{\rm{R}}^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}{{\boldsymbol{\tilde g}}_{{\rm{SR}}}}$表示多径导致的散射分量，${{\boldsymbol{\tilde \varPhi }}_{\rm{R}}}$表示空间相关矩阵(Spatial Correlation Matrix, SCM)，${{\boldsymbol{\tilde g}}_{{\rm{SR}}}} \in {\mathbb{C}^{N \times 1}}$表示独立同分布的复高斯随机矢量${{\boldsymbol{\tilde g}}_{{\rm{SR}}}} \sim \mathbb{C}\mathbb{N}\left( {0,2b{I_N}} \right)$，$2b$为散射分量的平均功率。

2.1.2   R-D信道模型

跟前面类似，考虑自由空间路径损耗以及小尺度衰落的影响，无人机与地面第$l$簇中第$n$个用户${D_{l,n}}$之间的信道矢量可以表示为

其中，$L'$为自由空间路径损耗，其计算公式为$L' = {\lambda \mathord{\left/ {\vphantom {\lambda {4\pi {d_{l,n}}}}} \right. } {4\pi {d_{l,n}}}}$，${d_{l,n}}$为无人机到地面用户${D_{l,n}}$的距离。此外，${{\boldsymbol{g}}_{l,n}}$为小尺度衰落，可以表示为^[10]

其中，信道系数${\rho _{l,n}}$是服从衰落参数为${m_{l,n}}$，平均功率为${\varOmega _{l,n}}$的Nakagami-m分布的随机变量。此外，${\boldsymbol{a}}\left( {{\theta _{l,n}}} \right) \in {\mathbb{C}^{N \times 1}}$表示阵列的导向矢量，其中${\theta _{l,n}}$为无人机到地面用户${D_{l,n}}$的离开角。

2.2   信号模型

如图1所示，卫星首先将叠加信号${x_{{\rm{SR}}}} = \displaystyle\sum\nolimits_{l = 1}^L {{x_l}}$发送给无人机，无人机对接收信号进行波束成形后，其输出信号可以表示为

其中，${P_{\rm{S}}}$为卫星发射功率，${x_{l,n}}$和${x_{j,k}}$分别为用户${D_{l,n}}$与用户${D_{j.k}}$的信号，${\alpha _{l,n}}$和${\alpha _{j,k}}$分别为用户${D_{l,n}}$与用户${D_{j.k}}$的功率分配系数，满足$\displaystyle\sum\nolimits_{n = 1}^M {{\alpha _{l,n}}} = 1$, $\displaystyle\sum\nolimits_{k = 1}^M {{\alpha _{j,k}}} = 1$。此外，${{\boldsymbol{w}}_{{{\rm{SR}}} }} \in {\mathbb{C}^{N \times 1}}$是接收波束成形权矢量，${{\boldsymbol{n}}_{{\rm{SR}}}} \sim \mathbb{C}\mathbb{N}\left( {0,\sigma _{{\rm{SR}}}^2{I_{\rm{N}}}} \right)$为加性高斯白噪声(Additive White Gaussian Noise, AWGN)。

与译码转发协议(Decode and Forward, DF)相比，AF协议是一种更加容易实现的中继协议^[20]，更加适用于本文基于无人机辅助的卫星通信场景。因此，无人机中继采用AF协议，通过固定增益因子$G{\text{ = }}\sqrt {{{{P_{{\rm{R}}} }} \mathord{\left/ {\vphantom {{{P_{R} }} {\left( {{{\rm{E}}} \left[ {{P_S}{{\left| {{\boldsymbol{w}}_{SR}^H{{\boldsymbol{h}}_{{SR} }}} \right|}^2}} \right] + \sigma _{SR}^2} \right)}}} \right. } {\left( {{{\rm{E}}} \left[ {{P_{\rm{S}}}{{\left| {{\boldsymbol{w}}_{{\rm{SR}}}^{\rm{H}}{{\boldsymbol{h}}_{{{\rm{SR}}} }}} \right|}^2}} \right] + \sigma _{{\rm{SR}}}^2} \right)}}}$放大其接收信号后转发给地面用户，用户${D_{l,n}}$的接收信号表示为

其中，${{\boldsymbol{w}}_l} \in {\mathbb{C}^{N \times 1}}$是无人机对第$l$簇的发射波束成形权矢量，${P_{\rm{R}}}$是无人机中继的发射功率，${n_{l,n}}$是均值为0，方差为$\sigma _{l,n}^2$的AWGN。

为了进一步消除簇内用户干扰，采用串行干扰消除(Successive Interference Cancellation, SIC)技术消除干扰。假设第$l$簇中用户的信道增益大小满足${\left\| {{{\boldsymbol{h}}_{l,1}}} \right\|^2} \lt {\left\| {{{\boldsymbol{h}}_{l,2}}} \right\|^2} \lt \cdots \lt {\left\| {{{\boldsymbol{h}}_{l,M}}} \right\|^2}$。针对无人机辅助的卫星通信系统下行链路，NOMA技术为信道质量差的用户分配较高的发送功率，而为信道质量好的用户分配较低的发送功率，由此用户${D_{l,n}}$主要受到的干扰来自其他簇以及第$l$簇中信道质量差的用户，则第$l$簇中用户的功率分配系数大小满足${\alpha _{l,1}} \gt {\alpha _{l,2}} \gt \cdots \gt {\alpha _{l,M}}$。因此接收端通过SIC技术消除簇内用户干扰后，${D_{l,n}}$的输出SINR为

其中，${\bar \gamma _{{\rm{SR}}}}{\text{ = }}{{{P_{\rm{S}}}\zeta _{{\rm{SR}}}^2} \mathord{\left/ {\vphantom {{{P_S}\zeta _{SR}^2} {\sigma _{SR}^2}}} \right. } {\sigma _{{\rm{SR}}}^2}}$, ${\bar \gamma _{l,n}}{\text{ = }}{{{P_{{\rm{R}}} }{{L'}^{{\text{ }}2}}} \mathord{\left/ {\vphantom {{{P_{R} }{{L'}^{{\text{ }}2}}} {\sigma _{l,n}^2}}} \right. } {\sigma _{l,n}^2}}$, ${\varLambda _{{\rm{SR}}}} = {\text{E}}\left[ {{{\bar \gamma }_{{{\rm{SR}}} }}{{\left| {{\boldsymbol{w}}_{{\rm{SR}}}^{\rm{H}}{{\boldsymbol{g}}_{{\rm{SR}}}}} \right|}^2}} \right] + 1$。

2.3   无人机波束成形方案

在无线通信中，考虑到信道的随机性，一般通过使用户输出平均SINR最大来达到系统传输性能最优的目的。因此根据式(10)，本文建立以用户输出平均SINR最大化为准则的优化问题，在数学上可以表示为^[19]

考虑到式(11)中${\text{E}}\left[ {{\gamma _{l,n}}} \right]$是很难计算得到的，因此将其近似表示为^[19]

于是，优化问题式(11)可以重新表示为

其中，${{\boldsymbol{\varPhi }}_{{\rm{SR}}}} = {{\rm E}_{{{\boldsymbol{g}}_{{\rm{SR}}}}}}\left[ {{{\boldsymbol{g}}_{{\rm{SR}}}}{\boldsymbol{g}}_{{\rm{SR}}}^{\rm{H}}} \right]$, ${{\boldsymbol{\varPhi }}_{l,n}} = {{\rm E}_{{{\boldsymbol{g}}_{l,n}}}}\left[ {{{\boldsymbol{g}}_{l,n}}{\boldsymbol{g}}_{l,n}^{{\rm{H}}} } \right]$分别为${{\boldsymbol{g}}_{{\rm{SR}}}}$和${{\boldsymbol{g}}_{l,n}}$的信道相关矩阵(Channel Correlation Matrices, CCM)，CCM通常由信道的多次采样统计得到。但是考虑到无人机能量受限的特性，本文假设仅已知S-R和R-D链路的信道角度信息${\theta _{r} }$和${\theta _{l,n}}$，相比于假设准确CSI已知带来的高反馈负载，仅获取信道角度信息可降低无人机的能耗。进一步，根据式(3)和式(7)，可以计算得到${{\boldsymbol{\varPhi }}_{{\rm{SR}}}}$和${{\boldsymbol{\varPhi }}_{l,n}}$，分别表示为

考虑到优化问题式(13)的目标函数随${\boldsymbol{w}}_{{\rm{SR}}}^{\rm{H}}{{\boldsymbol{\varPhi }}_{{\rm{SR}}}}{{\boldsymbol{w}}_{{\rm{SR}}}}$单调递增，且${{\boldsymbol{w}}_{{\rm{SR}}}}$和${{\boldsymbol{w}}_l}$相互独立。因此，我们首先利用阵列信号处理技术来获得归一化的无人机接收波束成形权矢量${{\boldsymbol{w}}_{{\rm{SR}}}}$。然后，提出基于信道角度信息的迫零波束成形方案，得到无人机发射波束成形权矢量${{\boldsymbol{w}}_l}$。

首先，建立优化问题式(16)获得归一化的无人机接收波束成形权矢量${{\boldsymbol{w}}_{{{\rm{SR}}} }}$

利用式(14)，对优化问题式(16)求解可以得到无人机接收波束成形权矢量${\boldsymbol{w}}_{{\rm{SR}}}^*$为

然后，将式(17)代入优化问题式(13)，可以将优化问题式(13)简化为

其中，$Q = {\boldsymbol{a}}^{{\rm{H}}}{\left( {{\theta _r}} \right) }{{\boldsymbol{\varPhi }}_{{\rm{SR}}}}{\boldsymbol{a}}\left( {{\theta _r}} \right)$。

由于优化问题式(18)的目标函数中${{\boldsymbol{w}}_l}$和${{\boldsymbol{w}}_j}$相互耦合，该问题数学上难以直接求解。于是，为了降低计算复杂度，本文提出一种次优方案。通过基于信道角度信息的迫零波束成形方案消除簇间干扰以实现无人机多波束全频率复用，提高系统频谱效率。因此优化问题式(18)可以重新表示为

利用式(15)并根据矩阵零空间理论，对优化问题式(19)进行求解可以得到${\boldsymbol{w}}_l^*$为

其中，${{\boldsymbol{H}}_l}{\text{ = }}\left[ {{{{\boldsymbol{\hat H}}}_1},{{{\boldsymbol{\hat H}}}_2}, \cdots ,{{{\boldsymbol{\hat H}}}_{l - 1}},{{{\boldsymbol{\hat H}}}_{l + 1}}, \cdots ,{{{\boldsymbol{\hat H}}}_L}} \right]$, ${{\boldsymbol{\hat H}}_1} = \left. {\left[ {{\boldsymbol{a}} \left( {{\theta _{l,1}}} \right),} \right.{\boldsymbol{a}}\left( {{\theta _{l,2}}} \right), \cdots ,{\boldsymbol{a}}\left( {{\theta _{l,M}}} \right)} \right]$, ${{\boldsymbol{\varPi }}_l} = {I_N} - {{\boldsymbol{H}}_l}{\left( {{\boldsymbol{H}}_l^{{\rm{H}}} {{\boldsymbol{H}}_l}} \right)^{ - 1}} {\boldsymbol{H}}_l^{{\rm{H}}}$。

于是，将式(17)和式(20)代入式(10)后，用户${D_{l,n}}$的输出SINR可以简化为

其中，${\varLambda '_{{\rm{SR}}}} = {\left\| {{{\boldsymbol{\pi }}_l}{\boldsymbol{a}}\left( {{\theta _{l,n}}} \right)} \right\|^2}\left\{ {{\rm E}\left[ {{{\bar \gamma }_{{{\rm{SR}}} }}{{\left| {{\boldsymbol{a}}{{\left( {{\theta _r}} \right)}^{\rm{H}}}{{\boldsymbol{g}}_{{{\rm{SR}}} }}} \right|}^2}} \right] + N} \right\}$。

需要指出的是，跟大多数采用基于准确CSI的波束成形方案，例如文献[10]不同的是，本文利用信道角度信息进行波束成形设计，避免了信道估计和反馈等需要额外消耗无线资源的过程，从而更加适合于无人机通信场景。接下来，将进一步对系统的中断性能进行分析。

3.   系统的中断性能分析

中断概率是衡量无线通信服务质量(Quality of Service, QoS)的一项重要指标，定义为信号输出SINR低于某一特定门限值的概率。因此，根据式(21)用户${D_{l,n}}$的中断概率表示为

其中，$u = {{N{\gamma _{{\rm{th}},n}}}/ {{{\bar \gamma }_{{\rm{SR}}}}\left( {{\alpha _{l,n}} - {\gamma _{{\rm{th}},n}}\displaystyle\sum\nolimits_{i = n + 1}^M {{\alpha _{l,i}}} } \right)}}$, $A = {\bar \gamma _{{\rm{SR}}}}{\bar \gamma _{l,n}}\left( {{\alpha _{l,n}} - {\gamma _{{\rm{th}},n}}\displaystyle\sum\nolimits_{i = n + 1}^M {{\alpha _{l,i}}} } \right)$, $X = \left| {\boldsymbol{a}}^{\rm{H}}{{\left( {{\theta _r}} \right)}} {{\boldsymbol{g}}_{{\rm{SR}}}} \right|^2$, ${Y_{l,n}} = {\left| {{\rho _{l,n}}} \right|^2}$, ${\varLambda ''_{{{\rm{SR}}} }} = {\rm E}\left[ {{{\bar \gamma }_{{\rm{SR}}}}{{\left| {{\boldsymbol{a}}{{\left( {{\theta _r}} \right)}^{{\rm{H}}} }{{\boldsymbol{g}}_{{\rm{SR}}}}} \right|}^2}} \right] + N$。

为了进一步得到${I_1}$的闭合表达式，首先需要得到随机变量$X$的概率密度函数(Probability Density Function, PDF)，并进一步得到累积分布函数(Cumulative Distribution Function, CDF)。根据文献[19]，随机变量$X$的PDF表达式可以表示为

其中，${b_{\rm{R}}} = b{\left\| {{\boldsymbol{a}}^{{\rm{H}}}{{\left( {{\theta _r}} \right)} }{\boldsymbol{\tilde \varPhi }}_{{\rm{R}}} ^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} \right\|^2}$, ${}_1{F_1}\left( {{m_{\rm{R}}},1,{a_2}x} \right) = \exp \left( {{a_2}x} \right)\;\displaystyle\sum\nolimits_{{t_1} = 0}^{{m_{{\rm{R}}} } - 1} {\left[ {{{{{\left( {1 - {m_{{\rm{R}}} }} \right)}_{{t_1}}}\;{{\left( { - {a_2}x} \right)}^{{t_1}}}} \mathord{\left/ {\vphantom {{{{\left( {1 - {m_{{\rm{R}}} }} \right)}_{{t_1}}}{{\left( { - {a_2}x} \right)}^{{t_1}}}} {\left( {{t_1}!{{\left( 1 \right)}_{{t_1}}}} \right)}}} \right. } {\left( {{t_1}!{{\left( 1 \right)}_{{t_1}}}} \right)}}} \right]}$，${m_{\rm{R}}}$为整数，${a_1}$和${a_2}$分别表示为

进一步，可以得到随机变量$X$的CDF表达式

其中，$\delta = {1 /{\left( {2{b_{\rm{R}}}} \right) - {a_2}}}$。

另一方面，由文献[21]中式(3.351.3)和式(23)，可将式(22)中${\varLambda ''_{{{\rm{SR}}} }}$进一步计算得到

接下来，为了计算${I_2}$的闭合表达式，需要推导式(22)中${F_{{Y_{l,n}}}}\left( y \right)$为随机变量${Y_{l,n}}$的CDF。根据顺序统计量原理，${F_{{Y_{l,n}}}}\left( y \right)$可以表示为

其中，${b_n} = {{M!} \mathord{\left/ {\vphantom {{M!} {\left( {\left( {n - 1} \right)!\left( {M - n} \right)!} \right)}}} \right. } {\left( {\left( {n - 1} \right)!\left( {M - n} \right)!} \right)}}$。由于${\left| {{\rho _{l,n}}} \right|^2}$服从伽马分布，因此${F_{{{\left| {{\rho _{l,n}}} \right|}^2}}}\left( y \right)$表示为

将式(29)代入式(28)中，${F_{{Y_{l,n}}}}\left( y \right)$可以进一步表示为

将式(30)和式(23)(代入式(22)中，${I_2}$可以表示为

其中，${y_1} = {{{\gamma _{{\rm{th}},n}}{{\varLambda ''}_{{\rm{SR}}}}{Z_a}} / {\left( {Ax - N{\gamma _{{{\rm{th}}} ,n}}{{\bar \gamma }_{l,n}}} \right)}}$, ${Z_a} = {{{{{\left\| {{{\boldsymbol{\varPi }}_l}{\boldsymbol{a}}\left( {{\theta _{l,n}}} \right)} \right\|}^2}} / {\left| {{\boldsymbol{a}}^{\rm{H}}{{\left( {{\theta _{l,n}}} \right)}}{{\boldsymbol{\varPi }}_l}{\boldsymbol{a}}\left( {{\theta _{l,n}}} \right)} \right|}}^2}$。另外，式(31)中求和项$\displaystyle\sum\nolimits_{l = 0}^{{m_{l,n}} - 1} {{{{{\left( {{{{m_{l,n}}{y_1}} / {{\varOmega _{l,n}}}}} \right)}^l}}/ {l!}}}$为一个$\left( {{m_{l,n}} - 1} \right)$次的多项式。由文献[21]中式(0.314)，可以简化为

其中，$z = {{{m_{l,n}}{y_1}}/{{\varOmega _{l,n}}}}$。而$c_l^j$由式(33)—式(35)递归计算得

将式(32)代入式(31)中，${I_2}$进一步表示为

令式(36)中$B = Ax - N{\gamma _{{\rm{th}},n}}{\bar \gamma _{l,n}}$，并根据文献[21]中式(3.471.9)，${I_2}$的闭合式可以表示为

将式(23)、式(26)和式(37)代入式(22)，用户${D_{l,n}}$的中断概率闭合表达式最终表示为

4.   系统的渐进性能分析

为了进一步揭示关键参数与系统性能之间的影响，在得到中断概率闭合表达式的基础上，进一步讨论高信噪比下中断概率的渐进表达式。用户${D_{l,n}}$中断概率的近似表达式表示为

对式(23)进行级数展开，可以得到

进一步得到$X$的CDF渐进表达式，即表示为

由式(29)，$F_{{{\left| {{\rho _{l,n}}} \right|}^2}}^\infty \left( y \right)$可以表示为

将式(42)代入式(28)中，可以进一步得到$F_{{Y_{l,n}}}^\infty \left( y \right)$表示为

根据式(40)和式(43)，$I_2^\infty$进一步表示为

其中，${\varLambda ''_{{{\rm{SR}}} }} \approx N + \dfrac{{{{\bar \gamma }_{{{\rm{SR}}} }}{a_1}}}{{{\delta ^2}}}$，积分部分$\displaystyle\int_u^\infty \left( Ax - N{\gamma _{{\rm{th}},n}}{{\bar \gamma }_{l,n}} \right)^{ - {m_{l,n}} - n}{\rm{d}}x$可以进一步表示为

将式(41)、式(44)和式(45)代入式(39)中，得到中断概率渐进表达式为

假设${\bar \gamma _{{\rm{SR}}}} = {\bar \gamma _{l,n}} = \bar \gamma$。当$\bar \gamma \to \infty$时，中断概率的渐进表达式为^[20]

其中，${G_{\rm{a}}}$和${G_{\rm{d}}}$分别表示系统的阵列增益和分集度。则式(46)可以改写为

其中

于是，系统的分集度和阵列增益分别为

从以上可以看出，卫星到无人机中继以及无人机中继到地面用户的链路信道参数只影响系统的阵列增益，而对系统的分集度不产生影响。

5.   仿真验证

本节通过计算机仿真验证理论分析的正确性，同时定量分析了系统参数对卫星通信系统性能的影响。此外，为了验证本文所提传输方案的优越性，还跟传统的OMA方案和文献[17]中的方案进行对比。其中OMA方案是一个用户独占一个时间/频率资源块，文献[17]采用无人机配置单天线且簇内仅包含远近两个用户的下行NOMA传输方案。仿真中考虑系统包含6个用户，将所有地面用户分为两个簇，每个簇包含3个用户。假设S-R链路经历中度阴影衰落，信道参数为$\left\{ {{m_{\rm{R}}},b,{\varOmega _{\rm{R}}}} \right\} = \left\{ 5,0.126, 0.835 \right\}$，而R-D链路服从Nakagami-m分布，衰落参数为$\left\{ {{m_{l,1}},{m_{l,2}},{m_{l,3}}} \right\} = \left\{ {2,4,6} \right\}$, $\left\{ {{\varOmega _{l,1}},{\varOmega _{l,2}},{\varOmega _{l,3}}} \right\} = \left\{ 2,2, 2 \right\}$。此外，假设${P_{\rm{S}}} = {P_{\rm{R}}} = P$，噪声功率$\sigma _{{\rm{SR}}}^2 = \sigma _{l,n}^2 = \kappa {\rm{BT}}$, $\kappa = 1.38 \times {10^{ - 23}}\;{{\rm{J}}/{\rm{K}}}$为Boltzmann常数。其他主要的仿真参数设置见表1。

表 1  系统的参数设置

参数	数值
卫星波束最大增益$G_{\rm{S}}^{{\rm{max}}}{\text{ } }\left( { {\text{dB} } } \right)$	53
3 dB角度${\theta _{ { {\rm{S} } } ,3\;{\text{dB} } } }$(o)	0.4
带宽$B{\text{ }}\left( {{\text{MHz}}} \right)$	5
无人机与地面用户间距离${d_{l,n}}{\text{ }}\left( {\text{m}} \right)$	500
噪声温度$T{\text{ }}\left( {\text{K}} \right)$	300

下载: 导出CSV 

| 显示表格

图2给出了第2个簇中3个用户的中断概率随发射功率$P$的变化情况。仿真中假设无人机配置天线数为$N = 8$，无人机到用户的离开角为$\left\{ {{\theta _{2,1}},{\theta _{2,2}},{\theta _{2,3}}} \right\} = \left\{ {{{45}^\circ},{{60}^\circ},{{75}^\circ}} \right\}{\text{ }}$，功率分配系数为$\left\{ {{\alpha _{2,1}},{\alpha _{2,2}},{\alpha _{2,3}}} \right\} = \left\{ {0.6,0.3,0.1} \right\}$，目标速率为$\left\{ {{R_{{\rm{th1}}}},{R_{{\rm{th2}}}},{R_{{\rm{th3}}}}} \right\} = \left\{ 0.65, 1.2, 2.4 \right\}{{ {\rm{bit}}/({\rm{s}}\cdot{\rm{Hz}})}}$。从图2可以看出Monte Carlo仿真结果与系统的理论值相吻合，证明了系统理论分析的正确性。

图 2  不同用户中断概率随发射功率P的变化情况

下载: 全尺寸图片 幻灯片

图3给出了不同发射功率$P$的情况下，用户${D_{2,3}}$的中断概率随目标速率$R$的变化情况，仿真中$P$分别取$0{\text{ dBW}}$和$3{\text{ dBW}}$。从中可以看出，NOMA方案下的用户${D_{2,3}}$中断性能明显优于OMA方案。另外，在$P{\text{ = }}3{\text{ dBW}}$条件下系统的中断性能优于$P{\text{ = }}0{\text{ dBW}}$，这与图2中用户中断概率随发射功率$P$的增大而减小的结论是相符合的。

图 3  不同发射功率P下用户${D_{2,3}}$中断概率随目标速率的变化情况

下载: 全尺寸图片 幻灯片

图4为改变无人机配置天线数的条件下，用户${D_{2,3}}$中断概率随发射功率$P$的变化曲线，并与传统的OMA方案和文献[17]的方案进行比较。仿真中无人机配置天线数分别为$N = 8$和$N = 16$。从中可以看出，无人机配置天线数目越多，用户${D_{2,3}}$中断性能越好。这是因为多天线技术可以提供阵列增益，从而提高地面用户的接收信号强度，提升系统性能。另外，仿真结果表明，在考虑簇内包含3个用户的情况下，用户${D_{2,1}}$中断性能与文献[17]中簇内远用户的中断性能相接近。由此可以体现出本文工作的优越性。

图 4  无人机配置不同天线数下用户${D_{2,3}}$中断概率随发射功率P的变化情况

下载: 全尺寸图片 幻灯片

图5为用户${D_{2,1}}$中断概率随功率分配系数的变化趋势图，可以进一步分析功率分配系数对用户中断性能的影响。仿真中为保证功率分配系数满足${\alpha _{2,1}} \gt {\alpha _{2,2}} \gt {\alpha _{2,3}}$，设定${\alpha _{2,3}} = 0.1$，则${\alpha _{2,1}}$取值范围为0.51～0.89, ${\alpha _{2,2}} = 0.9 - {\alpha _{2,1}}$。从图中可以看出，随着${\alpha _{2,1}}$的增加，用户${D_{2,1}}$的中断概率减小。由此说明，提高用户功率分配系数，可以改善用户的中断性能。这是因为随着${\alpha _{2,1}}$的增加而增大了式(21)表示的用户${D_{2,1}}$的SINR，从而使其中断概率减小。

图 5  用户${D_{2,1}}$中断概率随功率分配系数的变化情况

下载: 全尺寸图片 幻灯片

6.   结束语

针对无人机辅助的卫星通信系统下行链路，分析了基于SDMA和协作NOMA相结合的多用户传输系统的中断性能。首先，配置多根天线的无人机作为中继站辅助卫星通信，采用NOMA技术服务多个地面用户，得到地面用户的输出SINR表达式。然后，基于用户平均SINR最大化准则，提出了利用信道角度信息的波束成形方案。接着，进一步推导得到系统的中断概率闭合表达式和高信噪比下系统的中断概率近似表达式。最后，计算机仿真验证了本文所提方案的优越性以及理论分析的正确性，并且定量分析了相关参数对用户性能的影响，为进一步探索NOMA技术在无人机辅助卫星通信系统中的应用提供了有益的参考。