存在I/Q不平衡的OFDM全双工双向译码转发中继系统及其性能分析

吴皓威; 赵俊波; 文格; 欧静兰

doi:10.11999/JEIT160545

摘要: 全双工技术可以使频谱利用率翻倍，是5 G系统的关键技术之一。采用直接变换结构的全双工系统中残余自干扰(Residual Self-Interference, RSI)和同相/正交(In-phase/Quadrature, I/Q)不平衡是限制系统性能的两大主要因素。该文针对存在I/Q不平衡的OFDM全双工双向中继系统，建立了译码转发中继模式下的全双工系统信号模型，分析了瑞利衰落信道下系统的中断性能，获得了系统中断概率的闭式表达式。仿真结果不仅验证了理论分析的正确性，还得到结论：随着I/Q不平衡程度和残余自干扰强度的降低，系统中断性能将得到改善；只有沿着最速下降路线降低I/Q不平衡或中继节点RSI，才能实现最优的性能提升；通过系统I/Q不平衡与RSI参数所在的坐标点和最速下降路线的相对位置关系，来确定改善全双工双向中继系统中断性能的最优措施。

关键词:

Abstract: The full-duplex transmission is one of the key technologies in the 5 G communication systems, due to the ability of improving spectrum efficiency. However, the performance of the full-duplex system, with the zero intermediate frequency structure, is badly impacted by the Residual Self-Interference (RSI) and In-phase/ Quadrature (I/Q) imbalance. In this paper, the OFDM full-duplex bidirectional relaying system under the RSI and I/Q Imbalance (IQI) is investigated, in a cooperative scenario where the Decode-and-Forward (DF) protocol is considered. The outage performance of the system and its closed-form expressions are derived under Rayleigh fading channels, and the influences of the IQI and RSI on system performances are analyzed, respectively. The simulation results verify the analysis, and the conclusions are given as follows. First, the outage performance improves as decreasing of the IQI and RSI. Second, the optimal improvement of outage performance is achieved by reducing the RSI and I/Q imbalance, according the route with the steepest descent method. Third, the best way for enhancing the outage performance is chosen, by the relative position between the steepest descent route and the current coordinate of IQI-RSI.

Key words:

1. 引言

星地融合网络(Integrated Satellite-Terrestrial Networks, ISTNs)具有覆盖范围广、传输数据通量高、可连接用户数多和抗毁性强等特性，是在遇到地震、台风等极端自然灾害之后重建通信链路的最有效途径^[1]。近年来，随着物联网(Internet-of-Things, IoT)、车联网(Internet-of-Vehicles, IoV)、远程医疗等其他技术的快速发展和急切需求，星地融合网络在学术界和工程领域都获得了巨大的关注^[2]。

服务质量(Quality of Service, QoS)和频谱效率是未来星地融合网络发展的两个关键因素。因此，为了提高频谱效率并增强用户的服务体验，非正交多址接入(Non-Orthogonal Multiple Access, NOMA)技术被引入到星地融合网络中^[3]。NOMA技术采用功率域多址方案，实现利用同一频率资源向多个传感器用户传输内容的目的，从而提高频谱效率，并通过功率分配因子提高传感器用户的服务质量，这有别于传统的正交多址接入( Orthogonal Multiple Access, OMA)技术^[4]。由于NOMA技术的优越性，许多研究分析了将NOMA技术应用于卫星通信(Satellite Communication, SatCom)后对于系统性能的提升。文献[5]研究了 NOMA辅助的星地融合网络的中断性能，并推导得到了系统的中断概率(Outage Probability, OP)的闭式解析表达式和渐近表达式。文献[6]建立了基于NOMA的卫星通信系统的功率分配优化模型，提高了传感器用户的服务质量。文献[7]将NOMA方案应用于基于内容分发的星地融合网络，并推导得到了中断概率和命中概率的解析表达式。

多天线技术同样是增强系统传输质量的关键技术。相比基于单天线的NOMA系统，多天线技术的应用可以有效提高系统的容量^[8]。文献[9]研究了多天线节点下的星地融合中继网络(Integrated Satellite-Terrestrial Relay Network, ISTRN)的性能，并验证了多天线技术对于系统性能提升的积极影响。文献[10]通过分析比较得出了多天线星地融合网络的中断概率更优。文献[11]讨论了多天线星地融合网络的可靠性和安全性，并通过数值结果验证了多天线技术的优越性。

在实际条件的影响下，基于NOMA的星地融合网络在传输和检测过程中会受到各种非完美因素的影响。在传输过程中，信道通常会经历严重的衰落，如雨、雾等天气的影响，因此，系统很难获得完美的信道状态信息(Channel State Information, CSI)^[12]。同时由于信道估计技术的限制，信道估计误差(Channel Estimation Errors, CEEs)不可避免地会出现在信道状态信息估计过程中^[13]。文献[14]在陆地移动卫星系统中考虑了不完美CSI，推导得到了系统的精确和渐近中断概率。文献[15]建立了非理想CSI下的卫星通信系统模型，并推导得到了系统的最优能量利用效率。此外，在卫星信号接收检测过程中，NOMA用户采用串行干扰消除(Successive Interference Cancellation, SIC)技术来获取叠加信号中每个传感器用户的目标信号^[16]。然而，由于接收器性能的限制，实际条件下很难实现完美的串行干扰消除^[17]。文献[18]分析了在非完美SIC下的基于NOMA的卫星通信系统性能，推导得出了每个传感器用户中断概率和渐近中断概率的解析表达式。然而，作者采用了一个独立的干扰因子来表示非完美SIC的影响，缺乏系统性的理论分析。此外，由于频谱资源的重复利用，同频干扰(Co-Channel Interference, CCI)同样会出现在信号传输过程中^[19]。文献[20]分别从遍历容量(Ergodic Capacity, EC)、中断概率、平均符号错误率和能量效率等方面评估了同频干扰对星地融合网络的影响。

针对以上问题，本文在非完美条件下，即信道估计误差，非完美SIC和同频干扰的影响下，对基于NOMA的星地融合网络的遍历容量性能进行了研究。本文的主要贡献如下：(1)建立了非完美条件下的星地融合NOMA网络的系统模型。由于接收机性能的限制，串行干扰消除是非完美的。由于非完美的信道状态信息，考虑了信道估计误差。同时，频率的重复利用导致了同频干扰。(2)通过遍历容量性能的理论推导，刻画了不同非完美条件参数对于星地融合NOMA网络的影响。(3)采用蒙特卡罗仿真进一步验证了遍历容量性能理论推导的正确性。

2. 系统模型和问题建模

如图1所示，本文考虑一个基于NOMA的多天线星地融合网络，其中地球同步卫星(Geosynchronous Earth Orbit, GEO)信号源 $S$ 采用NOMA技术通过直连链路与地面的传感器用户 ${U_i},i \in \left( {1,2, \cdots, p,q, \cdots ,\mathcal{N}} \right)$ 进行通信。地面传感器用户通过分组策略分为多个簇，不失一般性，本文针对一个卫星波束下的两个NOMA用户 ${U_i},i \in \left( {p,q} \right)$ 进行分析。两用户NOMA场景被3GPP(Third Generation Partnership Project)组织认可为可以增强传感器用户的频谱效率。因此，采用两用户NOMA方案，不仅可以简化系统模型并可以为未来研究多NOMA传感器用户场景提供基础^[5]。此外，假设两个地面传感器用户位于卫星多个波束中的同一个波束内。地面传感器用户为了增强天线增益配置了 $M$ 根天线。同时，由于地面设备频率复用，每一个传感器用户都会受到 ${I_c},c \in \left\{ {1,2, \cdots ,N} \right\}$ 个干扰的影响。

图 1 多天线星地融合NOMA网络

下载: 全尺寸图片幻灯片

2.1 信道模型

在星地融合网络中，采取阵列馈电反射(Array Fed Reflector, AFR)技术固定每个天线的辐射模式来减少控制系统的处理消耗。相比较直接辐射阵列(Direct Radiating Array, DRA)技术，阵列馈电反射技术可以获得更高的天线增益和能量效率^[6]。考虑到自由空间损耗、雨衰和天线增益，星地链路的信道分量表达式为

${{\boldsymbol{h}}_{{{\rm{SU}}_i}}} = \sqrt {{\mathcal{F}_i}{\mathcal{H}_i}} {{\boldsymbol{\zeta}} _i}^{ - \frac{1}{2}} \cdot {\boldsymbol{f}}_i^{\frac{1}{2}}{{\boldsymbol{g}}_{{{\rm{SU}}_i}}}$

(1)

其中， ${\mathcal{F}_i}$ 代表自由空间损耗，并可以表示为

${\mathcal{F}_i}{\text{ = }}{\left( {\frac{\mathcal{V}}{f}} \right)^2}\frac{1}{{d_{\text{G}}^2 + d_i^2}}$

(2)

其中， $\mathcal{V}{\text{ = }}{{\rm{c}}\mathord{\left/ {\vphantom {c {4{\text{π }}}}} \right. } {4{\pi}}}$ ， ${\rm{c}}$ 代表电磁波速度， $f$ 代表载波频率， ${d_{\text{G}}} \approx {\text{35786 km}}$ 代表GEO卫星高度， ${d_i}$ 代表传感器用户 ${U_i}$ 到卫星波束覆盖范围中心的距离。 ${\mathcal{H}_i}$ 代表地面传感器用户 ${U_i}$ 的天线增益，表达式为

${\mathcal{H}_i} \simeq \left\{ \begin{aligned} &{{\mathcal{H}_{\max }} - 2.5 \times {{10}^{ - 3}}{{\left( {\frac{{{d_{\rm{a}}}\theta }}{\lambda }} \right)}^2},}\;{0^\circ < \theta < {\theta _{\rm{a}}}} \\ &{2 + 15\lg \frac{{{d_{\text{a}}}}}{\lambda },}\;\qquad\qquad\qquad\;\, {{\theta _{\text{a}}} < \theta < {\theta _{\text{b}}}} \\ & {32 - 25\lg \theta ,}\qquad\qquad\qquad\;\;\; {{\theta _{\text{b}}} < \theta < 48^\circ } \\ & { - 10,}\qquad\qquad\qquad\qquad \quad\;\; \; {48^\circ < \theta < 180^\circ } \end{aligned} \right.$

(3)

其中， ${\mathcal{H}_{\max }}$ 代表用户天线的最大增益， ${d_{\text{a}}}$ 代表天线直径， $\lambda$ 代表信号波长， $\theta$ 代表离轴角度，其中 ${\theta _{\text{a}}} = \dfrac{{20\lambda }}{{{d_{\text{a}}}}}\;\sqrt {{\mathcal{H}_{\max }} - \left( {2 \;+\; 15\lg \dfrac{{{d_{\text{a}}}}}{\lambda }} \right)}$ , ${\theta _{\text{b}}} = 15.85 \times {\left( {{{{d_{\text{a}}}}}/{\lambda }} \right)^{ - 0.6}}$ 。此外， ${{\boldsymbol{\zeta}} _i}{\text{ = }}{\left[ {{\zeta _{i.1}},{\zeta _{i.2}}, \cdots ,{\zeta _{i.M}}} \right]^{\rm{T}}}$ 代表雨衰分量， ${\zeta _{i.m}} \sim {{\rm{CN}}}\left( {\mu ,\delta _\zeta ^2} \right)$ 服从对数随机分布， $m \in \left\{ M \right\}$ 。 ${{\boldsymbol{f}}_i}{\text{ = }}{\left[ {{{f}_{i.1}},{{f}_{i.2}}, \cdots ,{{f}_{i.M}}} \right]^{\text{T}}}$ 代表卫星天线增益向量，其中 ${{f}_{i.m}}$ 表达式为

${{f}_{i.m}} \simeq {{f}_{\max }}{\left(\frac{{{B_1}\left( {{r_{i,m}}} \right)}}{{2{r_{i,m}}}} + 36\frac{{{B_3}\left( {{r_{i,m}}} \right)}}{{r_{i,m}^3}}\right)^2}$

(4)

其中， ${{f}_{\max }}$ 代表卫星天线最大增益， ${B_1}\left( \cdot \right)$ 和 ${B_3}\left( \cdot \right)$ 分别表示次序1和3的第1类贝塞尔函数。 ${r_{i,m}} = {{\tau \sin {\varphi _{i,m}}} \mathord{\left/ {\vphantom {{\tau \sin {\varphi _{i,m}}} {\sin {\varphi _{3\;{\text{dB}}}}}}} \right. } {\sin {\varphi _{3\;{\text{dB}}}}}}$ , $\tau = 2.07123$ ，其中 ${\varphi _{i,m}}$ 代表第 $m$ 个天线波束中心与用户 ${U_i}$ 和卫星直连线路中心的夹角， ${\varphi _{3\;{\text{dB}}}}$ 代表3 dB角度。最后， ${{\boldsymbol{g}}_{{\text{S}}{{\text{U}}_i}}}$ 代表星地链路随机信道分量。目前存在许多的数学模型表示星地链路信道信息，例如Lutz模型、Markov模型以及Karasaw模型。本文采用阴影莱斯(Shadowed-Rician, SR)模型来描绘星地信道。SR分布的模型与实测数据十分吻合，同时计算复杂度相对较低。因此很多文献采用SR分布作为星地融合网络信道模型^[5]。根据SR分布， ${{\boldsymbol{g}}_{{\text{S}}{{\text{U}}_i}}}$ 的第 $m$ 个分量可以表示为

${g_{{\text{S}}{{\text{U}}_i},m}} = X\exp \left( {{\rm{j}}\varsigma } \right) + Y\exp \left( {{\rm{j}}\vartheta } \right)$

(5)

其中， $X$ 和 $Y$ 均是独立平稳随机过程，并分别代表直连链路和多径分量的振幅。直连链路和多径分量分别服从Nakagami- $m$ 分布和瑞利分布。 $\varsigma$ 代表直连链路的确定性分量， $\vartheta$ 代表平稳随机分量，并服从 $[0,2{\pi})$ 。

通过数学转换， ${g_{{\text{S}}{{\text{U}}_i},m}}$ 的平方振幅的概率密度函数(Probability Density Function, PDF)可以表示为

$\begin{split} {f_{{{\left\| {{g_{{\text{S}}{{\text{U}}_i},m}}} \right\|}^2}}}\left( x \right) =& {\alpha _{{\text{S}}{{\text{U}}_i},m}}{{\rm{e}}^{ - {\beta _{{\text{S}}{{\text{U}}_i},m}}x}}\\ & \cdot {}_1{F_1}\left( {{m_{{\text{S}}{{\text{U}}_i},m}};1;{\delta _{{\text{S}}{{\text{U}}_i},m}}x} \right) \end{split}$

(6)

其中， ${\alpha _{{\text{S}}{{\text{U}}_i},m}}{\text{ }} \triangleq {\text{ }}\dfrac{{{{\left( {\dfrac{{2{b_{{\text{S}}{{\text{U}}_i},m}}{m_{{\text{S}}{{\text{U}}_i},m}}}}{{2{b_{{\text{S}}{{\text{U}}_i},m}}{m_{{\text{S}}{{\text{U}}_i},m}} + {\varOmega _{{\text{S}}{{\text{U}}_i},m}}}}} \right)}^{{m_{{\text{S}}{{\text{U}}_i},m}}}}}}{{2{b_{{\text{S}}{{\text{U}}_i},m}}}}$ , ${\delta _{{\text{S}}{{\text{U}}_i},m}} \triangleq \dfrac{{{\varOmega _{{\text{S}}{{\text{U}}_i},m}}}}{{2{b_{{\text{S}}{{\text{U}}_i},m}}\left( {2{b_{{\text{S}}{{\text{U}}_i},m}}{m_{{\text{S}}{{\text{U}}_i},m}} + {\varOmega _{{\text{S}}{{\text{U}}_i},m}}} \right)}}$ , ${\beta _{{\text{S}}{{\text{U}}_i},m}} \triangleq \dfrac{1}{{2{b_{{\text{S}}{{\text{U}}_i},m}}}}$ , ${\varOmega _{{\text{S}}{{\text{U}}_i},m}},2{b_{{\text{S}}{{\text{U}}_i},m}}$ 和 ${m_{{\text{S}}{{\text{U}}_i},m}} \in (0,\infty )$ 分别代表直连分量的平均功率，多径分量的平均功率和衰落程度。

假设 ${m_{{\text{S}}{{\text{U}}_i},m}}$ 为正整数，通过代数变换^[21]，概率密度函数可以重新表示为

$\begin{split} {f_{{{\left\| {{g_{{\text{S}}{{\text{U}}_i},m}}} \right\|}^2}}}\left( x \right) =& {\alpha _{{\text{S}}{{\text{U}}_i},m}}\sum\limits_{k = 0}^{{m_{{\text{S}}{{\text{U}}_i},m}} - 1} {\zeta \left( k \right)} {x^k}\\ & \cdot{{\rm{e}}^{ - \left( {{\beta _{{\text{S}}{{\text{U}}_i},m}} - {\delta _{{\text{S}}{{\text{U}}_i},m}}} \right)x}} \end{split}$

(7)

其中， $\zeta \left( k \right) = {\left( { - 1} \right)^k}{\left( {1 - {m_{{\text{S}}{{\text{U}}_i},m}}} \right)_k}\delta _{{\text{S}}{{\text{U}}_i},m}^k/{\left( {k!} \right)^2}$ 。 ${\left( \cdot \right)_k}$ 代表Pochhammer符号。

2.2 信号模型

$S$ 采用叠加编码技术将叠加信号发送给两个地面传感器用户 ${U_i}$ ，信号表达式为

$x = \sqrt {{a_p}} {s_p} + \sqrt {{a_q}} {s_q}$

(8)

其中， ${s_i}$ 代表不同用户的目标信号，并满足 ${\rm{E}}[|{s_i}{|^2}] = 1$ 。 ${a_i}$ 是NOMA技术中的不同信号的功率分配系数，并满足 ${a_p} + {a_q} = 1$ 。本场景假设用户 ${U_p}$ 相比 ${U_q}$ 经历更严重的衰落，因此，更多的信号功率将分配给信道状态更差的用户，即 ${a_p} > {a_q}$ 。因此，用户接收到的信号可以表示为

$\begin{split} {y_i} =& {{\boldsymbol{h}}_{{\text{S}}{{\text{U}}_i}}}{\boldsymbol{w}}_1^{\rm{H}}\sqrt {{P_{\text{S}}}} (\sqrt {{a_p}} {s_p} + \sqrt {{a_q}} {s_q}) \\ & + \sum\limits_{I = 1}^N {{{\boldsymbol{g}}_{iI}}{\boldsymbol{w}}_1^{\rm{H}}\sqrt {{P_{\text{I}}}} {s_{iI}}} + {\boldsymbol{w}}_1^{\rm{H}}{{\boldsymbol{n}}_i} \end{split}$

(9)

其中， ${P_{\text{S}}}$ 代表卫星信号的传输功率， ${P_{\text{S}}} = \sigma P$ ， $P$ 代表卫星整体功率， $\sigma \in \left( {0,1} \right)$ 代表功率系数。 ${{\boldsymbol{h}}_{{\text{S}}{{\text{U}}_i}}}$ 代表卫星和用户间的信道分量，满足SR分布。 $\left\| {{{\boldsymbol{w}}_1}} \right\| = 1$ 代表波束成形分量。 ${s_{iI}}$ 代表干扰信号， ${P_{\text{I}}}$ 代表干扰信号功率。 ${{\boldsymbol{g}}_{iI}}$ 代表干扰和用户之间的信道分量，服从瑞利分布^[9]。 ${{\boldsymbol{n}}_i}$ 代表卫星和用户之间的加性高斯白噪声并服从 ${n_{i,m}}{\text{～}}\rm{C}\rm{N}\left( {0,\delta _{i,m}^2} \right)$ 。 $\delta _{i,m}^2 = {{\rm{K}}_{\rm{B}}}{T_i}$ , ${{\rm{K}}_{\rm{B}}} = 1.380649 \times {10^{ - 23}}\;{\rm{J}}/{\rm{K}}$ 代表Boltzmann常数， ${T_i}$ 代表用户噪声温度。

2.3 问题建模

考虑到实际系统条件，星地链路的电磁环境和气候环境极其复杂，导致系统无法获得完美的CSI。因此在信道估计过程中会产生信道估计误差。采用均方误差算法，信道估计模型可以表示为

${h_{{\text{S}}{{\text{U}}_i},m}} = {\tilde h_{{\text{S}}{{\text{U}}_i},m}} + {e_{{\text{S}}{{\text{U}}_i},m}}$

(10)

其中， ${h_{{\text{S}}{{\text{U}}_i},m}}$ 代表星地链路的SR分布的实际振幅， ${\tilde h_{{\text{S}}{{\text{U}}_i},m}}$ 代表星地链路的SR分布的检测振幅， ${e_{{\text{S}}{{\text{U}}_i},m}}$ 代表星地链路的信道估计误差。 ${\tilde h_{{\text{S}}{{\text{U}}_i},m}}$ 和 ${e_{{\text{S}}{{\text{U}}_i},m}}$ 相互正交。 ${e_{{\text{S}}{{\text{U}}_i},m}}$ 服从复杂高斯分布，即 ${e_{{\text{S}}{{\text{U}}_i},m}}{\text{～}} {\rm{CN}} \left( {0,{{ \mathcal{\overline V} }_{{e_{{\text{S}}{{\text{U}}_i},m}}}}} \right)$ 。采用导频符号来估计星地链路的信道状态信息，则 ${e_{{\text{S}}{{\text{U}}_i},m}}$ 的方差可以表示为

$\begin{split} {{ \mathcal{\overline V} }_{{e_{{\text{S}}{{\text{U}}_i},m}}}} &= {\rm{E}}\left\{ {{{\left| {{h_{{\text{S}}{{\text{U}}_i},m}}} \right|}^2}} \right\} - {\rm{E}}\left\{ {{{\left| {{{\tilde h}_{{\text{S}}{{\text{U}}_i},m}}} \right|}^2}} \right\} \\ & = \frac{1}{{{\mathcal{L}_{{\text{S}}{{\text{U}}_i},m}}{{ \mathcal{\overline V} }_{{\text{S}}{{\text{U}}_i},m}} + 1}} \end{split}$

(11)

其中， ${\mathcal{L}_{{\text{S}}{{\text{U}}_i},m}}$ 代表导频符号长度， ${ \mathcal{\overline V} _{{\text{S}}{{\text{U}}_i},m}}$ 代表导频符号在星地链路的平均信噪比。将 ${ \mathcal{\overline V} _{{\text{S}}{{\text{U}}_i},m}}$ 写作最小均方误差形式为

${ \mathcal{\overline V} _{{\text{S}}{{\text{U}}_i},m}} = {\rm{E}}\left\{ {{\mathcal{V}_{{\text{S}}{{\text{U}}_i},m}}} \right\} = {P_{\text{e}}}E\left\{ {{{\left| {{h_{{\text{S}}{{\text{U}}_i},m}}} \right|}^2}} \right\}/\delta _{i,m}^2$

(12)

其中， ${P_{\text{e}}}$ 代表导频符号功率， ${P_{\text{e}}} = \left( {1 - \sigma } \right)P$ 。本文采用 ${ \mathcal{\overline V} _{{e_{{\text{S}}{{\text{U}}_i},m}}}}$ 来表示系统的信道估计误差。

传感器用户通过最大比合并技术来增强星地融合网络的性能，其表达式为

${{\boldsymbol{w}}_1} = \frac{{{{\boldsymbol{h}}_{{\text{SU}}}}}}{{{{\left\| {{{\boldsymbol{h}}_{{\text{SU}}}}} \right\|}^2}}}$

(13)

根据上述分析，可以推导得到，在用户 ${U_q}$ 处检测信号 ${s_p}$ 的信干噪比为

${\gamma _{{\text{S}}{{\text{U}}_{q - p}}}} = \frac{{{a_p}\sigma {\lambda _{{\text{S}}{{\text{U}}_q}}}}}{{{a_q}\sigma {\lambda _{{\text{S}}{{\text{U}}_q}}} + \sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + {\lambda _{{l_q}}} + 1}}$

(14)

其中， ${\lambda _{{\text{S}}{{\text{U}}_q}}} = {\left| {{{\tilde {\boldsymbol{h}}}_{{\text{S}}{{\text{U}}_q}}}{\boldsymbol{w}}_1^{\rm{H}}} \right|^2}P/\delta _q^2 = {\bar \gamma _{{\text{S}}{{\text{U}}_q}}}{\left\| {{{\tilde {\boldsymbol{h}}}_{{\text{S}}{{\text{U}}_q}}}} \right\|^2}$ , ${\bar \gamma _j}$ , $j \in \left\{ {{\text{S}}{{\text{U}}_p},{\text{S}}{{\text{U}}_q}} \right\}$ 代表卫星到不同用户的平均信噪比。 ${\lambda _{{l_q}}} = \displaystyle\sum\nolimits_{c = 1}^N {{{\left| {{{\boldsymbol{g}}_{lq}}{\boldsymbol{w}}_1^{\rm{H}}} \right|}^2}{P_l}/\delta _q^2} = {\bar \gamma _q}\sum\nolimits_{c = 1}^N {{{\left\| {{{\boldsymbol{g}}_{lq}}} \right\|}^2}}$ , ${\bar \gamma _i}$ , $i \in \left\{ {p,q} \right\}$ 代表干扰到不同用户的平均信噪比。

在NOMA方案中，采用SIC技术来检测和接收信号，用户 ${U_q}$ 先检测信号 ${s_p}$ ，将其分离后再检测目标信号 ${s_q}$ 。但由于接收机性能的限制，假设非完美SIC在接收端产生，此时， ${U_q}$ 处检测信号 ${s_q}$ 的信干噪比为

${\gamma _{{\text{S}}{{\text{U}}_q}}} = \frac{{{a_q}\sigma {\lambda _{{\text{S}}{{\text{U}}_q}}}}}{{{a_p}\xi \sigma {\lambda _{{\text{S}}{{\text{U}}_q}}} + \sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + {\lambda _{{l_q}}} + 1}}$

(15)

其中， $\xi$ 代表非完美SIC系数，满足 $\xi \in \left( {0,1} \right)$ 。

最后，本文在 ${U_p}$ 处检测信号 ${s_p}$ ，将信号 ${s_q}$ 直接视作干扰，则信干噪比为

${\gamma _{{\text{S}}{{\text{U}}_p}}} = \frac{{{a_p}\sigma {\lambda _{{\text{S}}{{\text{U}}_p}}}}}{{{a_q}\sigma {\lambda _{{\text{S}}{{\text{U}}_p}}} + \sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + {\lambda _{{l_p}}} + 1}}$

(16)

其中， ${\lambda _{{\text{S}}{{\text{U}}_p}}} = {\left| {{{\tilde {\boldsymbol{h}}}_{{\text{S}}{{\text{U}}_p}}}{\boldsymbol{w}}_1^{\rm{H}}} \right|^2}P/\delta _p^2 = {\bar \gamma _{{\text{S}}{{\text{U}}_p}}}{\left\| {{{\tilde {\boldsymbol{h}}}_{{\text{S}}{{\text{U}}_p}}}} \right\|^2}$ , ${\lambda _{{l_p}}} = \displaystyle\sum\nolimits_{c = 1}^N {{{\left| {{{\boldsymbol{g}}_{lp}}{\boldsymbol{w}}_1^{\rm{H}}} \right|}^2}{P_l}/\delta _p^2} = {\bar \gamma _p}\sum\nolimits_{c = 1}^N {{{\left\| {{{\boldsymbol{g}}_{lp}}} \right\|}^2}}$ 。

值得说明的是文献[14,15,18,20]都只是分析了非完美条件的特例。本文采用导频信号来估计信道估计误差，这有别于文献[14]。在此基础上，考虑了服从瑞利分布的同频干扰。此外，在文献[18]中，非完美的SIC干扰被看作一个独立的参数，本文将非完美SIC与实际信号处理过程相关联。由此可知，本文所研究的非完美条件与之前所研究的文章有很大区别。

3. 性能分析

本节首先将会给出阴影莱斯信道和瑞利信道的信道统计特征的概率密度函数和累积分布函数(Cumulative Distribution Function, CDF)，在此基础上，将推导得到系统遍历容量的闭式解。

3.1 信道统计特征

星地链路服从阴影莱斯分布，在最大比合并技术的基础上， ${\gamma _j} = {\bar \gamma _j}{\left| {{h_j}} \right|^2}$ 的概率密度函数为

${f_{{\gamma _j}}}\left( x \right) = \frac{{\alpha _j^{{M}}}}{{\varGamma \left( {{M}} \right)}}\sum\limits_{k = 0}^{{{M}}\left( {{m_j} - 1} \right)} {\frac{{\zeta {{\left( k \right)}_{{M}}}}}{{{{\left( {{{\bar \gamma }_j}} \right)}^{k + {{M}}}}}}{x^{k + {{M}} - 1}}{{\rm{e}}^{ - {\varDelta _j}x}}}$

(17)

其中， ${\varDelta _j} = {{{\beta _j} - {\delta _j}} \mathord{\left/ {\vphantom {{{\beta _j} - {\delta _j}} {{{\bar \gamma }_j}}}} \right. } {{{\bar \gamma }_j}}}$ 。

经过数学变换^[21]， ${\gamma _j}$ 的累积分布函数为

$\begin{split} {F_{{\gamma _j}}}\left( x \right) =& 1 - \frac{{\alpha _j^{{M}}}}{{\varGamma \left( {{M}} \right)}}\\ & \cdot\sum\limits_{{k_j} = 0}^{{{M}}({m_j} - 1)} {\sum\limits_{t = 0}^{{k_j}} {\frac{{{{\left( {{{M}}\left( {1 - {m_j}} \right)} \right)}_{{k_1}}}{{\left( { - {\delta _j}} \right)}^{{k_j}}}}}{{{k_j}!{{\left( {{{\bar \gamma }_j}} \right)}^{{k_1} + {{M}}}}t!\varDelta _j^{{k_j} - t + {{M}}}}}{x^t}} {{\rm{e}}^{ - {\varDelta _j}x}}} \end{split}$

(18)

干扰到用户的地面链路服从瑞利分布，概率密度函数和累积分布函数可以分别表示为

${f_{{\gamma _i}}}\left( x \right) = \sum\limits_{s = 1}^{\varrho \left( \mathcal{A} \right)} {\sum\limits_{n = 1}^{{\tau _l}\left( \mathcal{A} \right)} {{\chi _{s,n}}\left( \mathcal{A} \right)\frac{{\mu _{\left\langle s \right\rangle }^{ - n}}}{{\left( {n - 1} \right)!}}} } {x^{n - 1}}{{\rm{e}}^{ - \frac{x}{{{\mu _{\left\langle s \right\rangle }}}}}}$

(19)

${F_{{\gamma _i}}}\left( x \right) = 1 - \sum\limits_{s = 1}^{\varrho \left( \mathcal{A} \right)} {\sum\limits_{n = 1}^{{\tau _l}\left( \mathcal{A} \right)} {\sum\limits_{k = 0}^{n - 1} {\frac{{{\chi _{s,n}}\left( \mathcal{A} \right)}}{{k!}}{{\left( {\frac{x}{{{\mu _{\left\langle s \right\rangle }}}}} \right)}^k}{{\rm{e}}^{ - \frac{x}{{{\mu _{\left\langle s \right\rangle }}}}}}} } }$

(20)

其中， $\mathcal{A} = {\rm{diag}}\left( {{\mu _1},{\mu _2}, \cdots ,{\mu _{{L}}}} \right)$ , $L \in \left\{ {N,M} \right\}$ , $\left\{ {{\mu _s}} \right\}_{s = 1}^{{L}}$ 代表CCI链路的平均信噪比， $\varrho \left( \mathcal{A} \right)$ 代表 $\mathcal{A}$ 中不同对角线元素的个数， ${\mu _{\left\langle 1 \right\rangle }} > {\mu _{\left\langle 2 \right\rangle }} > \cdots > {\mu _{\left\langle {\varrho \left( \mathcal{A} \right)} \right\rangle }}$ 代表不同对角线元素的增加序列， ${\tau _l}\left( \mathcal{A} \right)$ 代表 ${\mu _{\left\langle s \right\rangle }}$ 的多样性， ${\chi _{s,n}}\left( \mathcal{A} \right)$ 代表 $\mathcal{A}$ 中第 $\left( {s,n} \right)$ 个统计分量。

3.2 遍历容量分析

遍历容量是经常用来评估星地融合网络性能的重要指标，是指整个通信网络所有信道的容量上限。基于NOMA的多天线星地融合网络的EC定义为不同地面接收机的信干噪比的平均瞬时互信息之和，其表示为

$\begin{split} {{\rm{EC}}_{{\text{SU}}}} = & \left\{ {{\rm{E}}\left[ {{{\log }_2}\left( {1 + {\gamma _{{\text{S}}{{\text{U}}_{q - p}}}}} \right)} \right] + {\rm{E}}\left[ {{{\log }_2}\left( {1 + {\gamma _{{\text{S}}{{\text{U}}_q}}}} \right)} \right]} \right.\\ & \left. { + {\rm{E}}\left[ {{{\log }_2}\left( {1 + {\gamma _{{\text{S}}{{\text{U}}_p}}}} \right)} \right]} \right\} \\[-10pt] \end{split}$

(21)

将不同用户处的信干噪比式(14)、式(15)和式(16)代入EC表达式中，经过数学变换后，EC可以表示为

$\begin{split} & {{{\rm{EC}}_{{\text{SU}}}} } \\ & = \frac{1}{{\ln 2}}\left\{ {\rm{E}}\left[ {\ln \left( {\frac{\sigma }{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}{\lambda _{{\text{S}}{{\text{U}}_q}}} + \frac{{{\lambda _{{l_q}}}}}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}} + 1} \right)} \right] \right.\\ & \quad - {\rm{E}}\left[ {\ln \left( {\frac{{{a_q}\sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}{\lambda _{{\text{S}}{{\text{U}}_q}}} + \frac{{{\lambda _{{l_q}}}}}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}} + 1} \right)} \right] \\ & \quad + {\rm{E}}\left[ {\ln \left( {\frac{{\left( {{a_q} + {a_p}\xi } \right)\sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}{\lambda _{{\text{S}}{{\text{U}}_q}}} + \frac{{{\lambda _{{l_q}}}}}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}} + 1} \right)} \right] \\ & \quad - {\rm{E}}\left[ {\ln \left( {\frac{{{a_p}\xi \sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}{\lambda _{{\text{S}}{{\text{U}}_q}}} + \frac{{{\lambda _{{l_q}}}}}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}} + 1} \right)} \right] \\ & \quad + {\rm{E}}\left[ {\ln \left( {\frac{\sigma }{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + 1}}{\lambda _{{\text{S}}{{\text{U}}_p}}} + \frac{{{\lambda _{{l_p}}}}}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + 1}} + 1} \right)} \right] \\ & \quad \left. - {\rm{E}} \left[ {\ln \left( {\frac{{{a_q}\sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + 1}}{\lambda _{{\text{S}}{{\text{U}}_p}}} + \frac{{{\lambda _{{l_p}}}}}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + 1}} + 1} \right)} \right] \right\} \end{split}$

(22)

为了简化分析，采用变量替换， ${A_1} = \dfrac{\sigma }{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}$ , ${A_2} = \dfrac{{{a_q}\sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}$ , ${A_3} = \dfrac{{\left( {{a_q} + \xi {a_p}} \right)\sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}$ , ${A_4} = \dfrac{{\xi {a_p}\sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_q}}}}} + 1}}$ , ${A_5} = \dfrac{\sigma }{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + 1}}$ , ${A_6} = \dfrac{{{a_q}\sigma }}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + 1}}$ , ${B_1} = \dfrac{1}{{\sigma {{\overline V }_{{e_{{\text{SU}}q}}}} + 1}}$ , ${B_2} = \dfrac{1}{{\sigma {{\overline V }_{{e_{{\text{S}}{{\text{U}}_p}}}}} + 1}}$ 。此时，EC公式可以表示为

$\begin{split} {{\rm{EC}}_{{\text{SU}}}} = & \frac{1}{{\ln 2}}\left\{ {\rm{E}}\left[ {\ln \left( {{A_1}{\lambda _{{\text{S}}{{\text{U}}_q}}} + {B_1}{\lambda _{{l_q}}} + {\text{1}}} \right)} \right] \right.\\ & - {\rm{E}}\left[ {\ln \left( {{A_2}{\lambda _{{\text{S}}{{\text{U}}_q}}} + {B_1}{\lambda _{{l_q}}} + {\text{1}}} \right)} \right] \\ & + {\rm{E}}\left[ {\ln \left( {{A_3}{\lambda _{{\text{S}}{{\text{U}}_q}}} + {B_1}{\lambda _{{l_q}}} + {\text{1}}} \right)} \right] \\ & - {\rm{E}}\left[ {\ln \left( {{A_4}{\lambda _{{\text{S}}{{\text{U}}_q}}} + {B_1}{\lambda _{{l_q}}} + {\text{1}}} \right)} \right] \\ & + {\rm{E}}\left[ {\ln \left( {{A_5}{\lambda _{{\text{S}}{{\text{U}}_p}}} + {B_2}{\lambda _{{l_p}}} + {\text{1}}} \right)} \right] \\ & \left.- {\rm{E}}\left[ {\ln \left( {{A_6}{\lambda _{{\text{S}}{{\text{U}}_p}}} + {B_2}{\lambda _{{l_p}}} + {\text{1}}} \right)} \right] \right\} \end{split}$

(23)

令 ${\lambda _{{{\rm{SU}}_q}}} = {x_1}$ , ${\lambda _{{l_2}}} = {y_1}$ , ${A_1}{\lambda _{{\text{S}}{{\text{U}}_q}}} + {B_1}{\lambda _{{l_2}}} = {z_1}$ ，采用数学变换可得

${f_{{z_1}}}\left( z \right) = \int\limits_0^\infty {\frac{1}{{{A_1}{B_1}}}{f_Y}\left( {\frac{{z - u}}{{{B_1}}}} \right)} {f_X}\left( {\frac{u}{{{A_1}}}} \right){\rm{d}}u$

(24)

将SR分布和瑞利分布表达式(17)、式(19)代入，可得

$\begin{split} {f_{{z_1}}}\left( z \right) =& {{\alpha _{{\text{S}}{{\text{U}}_q}}^{{M}}}/ {\varGamma \left( {{M}} \right)}}\\ & \cdot\sum\limits_{k = 0}^{{{M}}\left( {{m_{{{\rm{SU}}_q}}} - 1} \right)} \frac{{{{\left( {{{M}}\left( {1 - {m_{{\text{S}}{{\text{U}}_q}}}} \right)} \right)}_k}{{\left( { - {\delta _{{\text{S}}{{\text{U}}_q}}}} \right)}^k}}}{{{{\left( {k!} \right)}^2}{{\left( {{{\bar \gamma }_{{\text{S}}{{\text{U}}_q}}}} \right)}^{k + 1}}}}\\ & \cdot\sum\limits_{s = 1}^{\varrho \left( A \right)} {\sum\limits_{n = 1}^{{\tau _l}\left( A \right)} {\sum\limits_{t = 0}^{n - 1} {\frac{{{{\left( { - 1} \right)}^t}{\chi _{s,n}}\left( A \right)\mu _{\left\langle s \right\rangle }^{ - n}}}{{\left( {n - 1} \right)!}}} } } \\ & \cdot \left( {\begin{array}{*{20}{c}} {n - 1} \\ t \end{array}} \right)\left( {k + t + {{M}} - 1} \right)!\frac{1}{{A_1^{k + {{M}}}B_1^n}}\\ & \cdot{\left( {\frac{{{\varDelta _{{\text{S}}{{\text{U}}_q}}}}}{{{A_1}}} - \frac{1}{{{B_1}{\mu _{\left\langle s \right\rangle }}}}} \right)^{ - \left( {k + t + {{M}}} \right)}}{z^{n - 1 - t}}{{\rm{e}}^{ - \frac{z}{{{B_1}{\mu _{\left\langle s \right\rangle }}}}}} \end{split}$

(25)

同理可以得到其余部分对应的PDF。

同时，本文将EC公式改写为

$\begin{split} {{\rm{EC}}_{{\text{SU}}}} =& \frac{1}{{\ln 2}}\left\{ \underbrace {\int\limits_0^\infty {\ln \left( {z + 1} \right){f_{{z_1}}}\left( z \right){\rm{d}}z} }_{{E_1}}\right. \\ & - \underbrace {\int\limits_0^\infty {\ln \left( {z + 1} \right){f_{{z_2}}}\left( z \right){\rm{d}}z} }_{{E_2}} \\ &+ \underbrace {\int\limits_0^\infty {\ln \left( {z + 1} \right){f_{{z_3}}}\left( z \right){\rm{d}}z} }_{{E_3}} \\ & - \underbrace {\int\limits_0^\infty {\ln \left( {z + 1} \right){f_{{z_4}}}\left( z \right){\rm{d}}z} }_{{E_4}} \\ &+ \underbrace {\int\limits_0^\infty {\ln \left( {z + 1} \right){f_{{z_5}}}\left( z \right){\rm{d}}z} }_{{E_5}} \\ & \left.- \underbrace {\int\limits_0^\infty {\ln \left( {z + 1} \right){f_{{z_6}}}\left( z \right){\rm{d}}z} }_{{E_6}} \right\} \end{split}$

(26)

采用文献[22]中的公式8.4.6.5，并利用Meijer-G函数^[21]，可以得到

$\ln \left( {1 + z} \right) = G_{22}^{12}\left( {z\left| {\begin{array}{*{20}{c}} {1}&1 \\ 1&0 \end{array}} \right.} \right)$

(27)

将式(27)和式(25)代入 ${E_1}$ 中，利用文献[22]中的公式2.24.3.1和8.2.2.14， ${E_1}$ 可以重新表示为

$\begin{split} {E_1} =& {{\alpha _{{\text{S}}{{\text{U}}_q}}^{{M}}}/ {\varGamma \left( {{M}} \right)}}\\ & \cdot\sum\limits_{k = 0}^{{{M}}\left( {{m_{{{\rm{SU}}_q}}} - 1} \right)} \frac{{{{\left( {{{M}}\left( {1 - {m_{{\text{S}}{{\text{U}}_q}}}} \right)} \right)}_k}{{\left( { - {\delta _{{\text{S}}{{\text{U}}_q}}}} \right)}^k}}}{{{{\left( {k!} \right)}^2}{{\left( {{{\bar \gamma }_{{\text{S}}{{\text{U}}_q}}}} \right)}^{k + 1}}}}\\ & \cdot\sum\limits_{s = 1}^{\varrho \left( A \right)} {\sum\limits_{n = 1}^{{\tau _l}\left( A \right)} {\frac{{{\chi _{s,n}}\left( A \right)\mu _{\left\langle s \right\rangle }^{ - t}}}{{\left( {n - 1} \right)!}}} } \sum\limits_{t = 0}^{n - 1} {\left( {\begin{array}{*{20}{c}} {n - 1} \\ t \end{array}} \right)} {\left( { - 1} \right)^t} \\ & \cdot \left( {k + t + {{M}} - 1} \right)!\frac{1}{{A_1^{k + {{M}}}B_1^t}}\\ & \cdot{\left( {\frac{{{\varDelta _{{\text{S}}{{\text{U}}_q}}}}}{{{A_1}}} - \frac{1}{{{B_1}{\mu _{\left\langle s \right\rangle }}}}} \right)^{ - \left( {k + t + {{M}}} \right)}}\\ & \cdot G_{32}^{13}\left( {{B_1}{\mu _{\left\langle s \right\rangle }}\left| {\begin{array}{*{20}{c}} {1 - n + t,1,1} \\ {1,0} \end{array}} \right.} \right) \\[-21pt] \end{split}$

(28)

同理可以得到其他分量的表达式。

在此基础上，经过简单的整理，最终EC的表达式可以推导为

$\begin{split} {{\rm{EC}}_{{\text{SU}}}} =& \frac{1}{{\ln 2}}\frac{{\alpha _{{{\rm{SU}}_q}}^{{M}}}}{{\varGamma \left( {{M}} \right)}}\\ & \cdot\sum\limits_{k = 0}^{{{M}}\left( {{m_{{\text{S}}{{\text{U}}_q}}} - 1} \right)} \frac{{{{\left( {{{M}}\left( {1 - {m_{{\text{S}}{{\text{U}}_q}}}} \right)} \right)}_k}{{\left( { - {\delta _{{\text{S}}{{\text{U}}_q}}}} \right)}^k}}}{{{{\left( {k!} \right)}^2}{{\left( {{{\bar \gamma }_{{\text{S}}{{\text{U}}_q}}}} \right)}^{k + 1}}}}\\ & \cdot\sum\limits_{s = 1}^{\varrho \left( A \right)} {\sum\limits_{n = 1}^{{\tau _l}\left( A \right)} {\sum\limits_{t = 0}^{n - 1} {\frac{{{{\left( { - 1} \right)}^t}{\chi _{s,n}}\left( A \right)}}{{\left( {n - 1} \right)!\mu _{\left\langle s \right\rangle }^t}}} } } \\ & \cdot G_{32}^{13}\left( {{B_1}{\mu _{\left\langle s \right\rangle }}\left| {\begin{array}{*{20}{c}} {1 - n + t,1,1} \\ {1,0} \end{array}} \right.} \right) \\ & \cdot \left( {\begin{array}{*{20}{c}} {n - 1} \\ t \end{array}} \right)\left( {k + t + {{M}} - 1} \right)!\\ & \cdot\left\{ \frac{1}{{A_1^{k + {{M}}}B_1^t}}{{\left( {\frac{{{\varDelta _{{\text{S}}{{\text{U}}_q}}}}}{{{A_1}}} - \frac{1}{{{B_1}{\mu _{\left\langle s \right\rangle }}}}} \right)}^{ - \left( {k + t + {{M}}} \right)}} \right.\\ & - \frac{1}{{A_2^{k + {{M}}}B_1^t}}{{\left( {\frac{{{\varDelta _{{\text{S}}{{\text{U}}_q}}}}}{{{A_2}}} - \frac{1}{{{B_1}{\mu _{\left\langle s \right\rangle }}}}} \right)}^{ - \left( {k + t + {{M}}} \right)}} \\ & + \frac{1}{{A_3^{k + {{M}}}B_1^t}}{{\left( {\frac{{{\varDelta _{{\text{S}}{{\text{U}}_q}}}}}{{{A_3}}} - \frac{1}{{{B_1}{\mu _{\left\langle s \right\rangle }}}}} \right)}^{ - \left( {k + t + {{M}}} \right)}} \\ & \left. - \frac{1}{{A_4^{k + {{M}}}B_1^t}}\left( \frac{{{\varDelta _{{\text{S}}{{\text{U}}_q}}}}}{{{A_4}}} - \frac{1}{{{B_1}{\mu _{\left\langle s \right\rangle }}}} \right)^{ - \left( {k + t + {{M}}} \right)} \right\} \\ &+ \frac{1}{{\ln 2}}\frac{{\alpha _{{\text{S}}{{\text{U}}_p}}^{{M}}}}{{\varGamma \left( {{M}} \right)}} \\ &\cdot \sum\limits_{k = 0}^{{{M}}\left( {{m_{{\text{S}}{{\text{U}}_p}}} - 1} \right)} {\frac{{{{\left( {{{M}}\left( {1 - {m_{{\text{S}}{{\text{U}}_p}}}} \right)} \right)}_k}{{\left( { - {\delta _{{\text{S}}{{\text{U}}_p}}}} \right)}^k}}}{{{{\left( {k!} \right)}^2}{{\left( {{{\bar \gamma }_{{\text{S}}{{\text{U}}_p}}}} \right)}^{k + 1}}}}}\\ & \cdot\sum\limits_{s = 1}^{\varrho \left( A \right)} {\sum\limits_{n = 1}^{{\tau _l}\left( A \right)} {\sum\limits_{t = 0}^{n - 1} {\frac{{{{\left( { - 1} \right)}^t}{\chi _{s,n}}\left( A \right)}}{{\left( {n - 1} \right)!\mu _{\left\langle s \right\rangle }^t}}} } } \\ & \cdot G_{32}^{13}\left( {{B_2}{\mu _{\left\langle s \right\rangle }}\left| {\begin{array}{*{20}{c}} {1 - n + t,1,1} \\ {1,0} \end{array}} \right.} \right) \\ & \cdot \left( {\begin{array}{*{20}{c}} {n - 1} \\ t \end{array}} \right)\left( {k + t} \right)!\\ & \cdot\left\{ \frac{1}{{A_5^{k + {{M}}}B_2^t}}{{\left( {\frac{{{\varDelta _{{\text{S}}{{\text{U}}_p}}}}}{{{A_5}}} - \frac{1}{{{B_2}{\mu _{\left\langle s \right\rangle }}}}} \right)}^{ - \left( {k + t + {{M}}} \right)}} \right.\\ & \left.- \frac{1}{{A_6^{k + {{M}}}B_2^t}}{{\left( {\frac{{{\varDelta _{{\text{S}}{{\text{U}}_p}}}}}{{{A_6}}} - \frac{1}{{{B_2}{\mu _{\left\langle s \right\rangle }}}}} \right)}^{ - \left( {k + t + {{M}}} \right)}} \right\} \end{split}$

(29)

4. 数据校验

采用蒙特卡罗仿真来验证非完美条件下多天线星地融合NOMA网络的性能。仿真中采用的系统参数如表1所示。卫星信道参数如表2所示。不失一般性，假设 $\delta _p^2 = \delta _q^2 = {\delta ^2}$ , ${\bar \gamma _{{\text{S}}{{\text{U}}_p}}} = {\bar \gamma _{{\text{S}}{{\text{U}}_q}}} = \bar \gamma$ , ${\bar \gamma _1} = {\bar \gamma _2} = {\bar \gamma _I}$ , ${k_{{\text{S}}{{\text{U}}_p}}} = {k_{{\text{S}}{{\text{U}}_q}}} = k$ , ${\mathcal{L}_{{\text{S}}{{\text{U}}_p}}} = {\mathcal{L}_{{\text{S}}{{\text{U}}_q}}} = \mathcal{L}$ 。

表 1 系统参数

参数	数值
卫星轨道	GEO
载波频率	18 GHz
载波带宽	50 MHz
3 dB角度	0.4°
最大波束增益	48 dB
接收天线增益	4 dB
噪声温度	300 K
雨衰	$\mu = - 3.125$ , $\delta _\zeta ^2 = 1.591$

下载: 导出CSV

| 显示表格

表 2 信道参数

衰落程度	${m_j}$	${b_j}$	${\varOmega _j}$
重度衰落(Frequent Heavy Shadowing, FHS)	1	0.063	0.0007
中度衰落(Average Shadowing, AS)	5	0.251	0.2790
轻度衰落(Infrequent Light Shadowing, ILS)	10	0.150	1.2900

下载: 导出CSV

| 显示表格

图2给出了遍历容量随功率分配因子的变化趋势。首先，可以发现仿真值与理论值相一致，证明了理论推导的正确性。其次，仿真结果随着功率分配因子 ${a_p}$ 的增大而升高。然而，从实际角度来看，当 ${a_p}$ 趋于1时，传感器用户 ${U_q}$ 的QoS和公平性要求无法满足。同时，由于功率分配的减少，用户 ${U_q}$ 的遍历容量也会损失。基于实际公平性考虑， ${a_p}$ 的值通常选取在 $\left( {0.7,0.8} \right)$ 之间，其余仿真均采用 ${a_p} = 0.75$ 作为仿真参数。最后，仿真结果显示遍历容量随着衰落程度的减轻而增强，这与实际相一致。

图 2 遍历容量随功率分配因子的变化

下载: 全尺寸图片幻灯片

图3分析了遍历容量随非完美SIC系数的变化趋势，其中 $\xi = 0$ 表示完美SIC。非完美SIC对于遍历容量的影响十分明显，随着 $\xi$ 的增大，系统的遍历容量明显变差。这是因为NOMA传感器用户在接收端采用SIC技术来接收信号，由于接收端的性能较差，即SIC技术非完美，这将对系统的传输能力造成很大的负面影响，从而减少系统的遍历容量。

图 3 遍历容量随非完美SIC的变化

下载: 全尺寸图片幻灯片

图4给出了遍历容量随信道估计误差的变化。由图4(a)可以看出，当系统分配更多的功率用于信号传输时，系统的遍历容量将会更高。由图4(b)可知，随着导频信号长度的增加，系统的信道估计误差将会减小，从而提升系统的遍历容量。此外，随着系统平均信噪比的增强，导频信号长度对于遍历容量的影响会减少。通过对不同信道估计误差参数的仿真可以得出，当系统具有较高的平均信噪比时，可以获得更好的CSI减轻估计误差对于系统性能的影响。

图 4 遍历容量随信道估计误差的变化

下载: 全尺寸图片幻灯片

图5反映了不同CCI情况对于遍历容量的影响，其中 ${\bar \gamma _I} = - \infty$ 代表不存在同频干扰。首先，存在同频干扰时，系统的遍历容量会有明显的减弱。此外，随着干扰功率的增大以及干扰数量的增多，系统的遍历容量会明显减小，这与现实情况相一致。

图 5 遍历容量随同频干扰的变化

下载: 全尺寸图片幻灯片

5. 结论

本文结合现实条件，包括信道估计误差、非完美串行干扰消除和同频干扰，建立了一个基于NOMA的星地融合网络。这些非完美条件将会影响基于NOMA的多天线星地融合网络的遍历容量性能。为了揭示每个非完美条件参数对于系统遍历容量的影响，本文推导得到了非完美条件下的遍历容量解析表达式。最后，通过蒙特卡罗仿真对理论推导的正确性进行了验证。仿真结果揭示了每个非完美条件参数对于星地融合NOMA网络遍历容量性能的影响，从侧面反映了具备完美系统场景对于实际网络的重要性，对工程实践具有重要的指导作用。

ZHANG X, CHENG W, and ZHANG H. Full-duplex transmission in PHY and MAC layers for 5 G mobile wireless networks[J]. IEEE Wireless Communications, 2015, 22(5): 112-121. doi: 10.1109/MWC.2015.7306545.

ZHANG Z, LONG K, VASILAKOS AV, et al. Full-duplex wireless communications: Challenges, solutions, and future research directions[J]. Proceedings of the IEEE, 2016, 104(7): 1369-1409. doi: 10.1109/JPROC.2015.2497203.

ZHANG Z, CHAI X, LONG K, et al. Full duplex techniques for 5 G networks: Self-interference cancellation, protocol design, and relay selection[J]. IEEE Communications Magazine, 2015, 53(5): 128-137. doi: 10.1109/MCOM.2015. 7105651.

HONG S, BRAND J, CHOI J, et al. Applications of self-interference cancellation in 5 G and beyond[J]. IEEE Communications Magazine, 2014, 52(2): 114-121. doi: 10. 1109/MCOM.2014.6736751.

SYRJALA V, VALKAMA M, ANTTILA L, et al. Analysis of oscillator phase-noise effects on self-interference cancellation in full-duplex OFDM radio transceivers[J]. IEEE Transactions on Wireless Communications, 2014, 13(6): 2977-2990. doi: 10.1109/TWC.2014.041014.131171.

KORPI D, RIIHONEN T, SYRJALA V, et al. Full-duplex transceiver system calculations: Analysis of ADC and linearity challenges[J]. IEEE Transactions on Wireless Communications, 2014, 13(7): 3821-3836. doi: 10.1109/ TWC.2014.2315213.

HORLIN F and BOURDOUX A. Digital Compensation for Analog Front-Ends[M]. England: Wiley, 2008: 71-95.

GUPTA A and JHA R K. A survey of 5 G network: Architecture and emerging technologies[J]. IEEE Access, 2015, 3: 1206-1232. doi: 10.1109/ACCESS.2015.2461602.

RATAJCZAK K, BAKOWSKI K, and WESOLOWSKI K. Two-way relaying for 5 G systems: Comparison of network coding and MIMO techniques[C]. IEEE Wireless Communications and Networking Conference (WCNC), Istanbul, Turkey, 2014: 376-381. doi: 10.1109/WCNC.2014. 6952037.

欧静兰, 吴皓威, 邹玉涛, 等. 过时信道状态下机会双向中继选择算法[J]. 北京邮电大学学报, 2014, 37(6): 44-48. doi: 10.13190/j.jbupt.2014.06.009.

OU Jinglan, WU Haowei, ZOU Yutao, et al. Opportunistic two-way relay selection scheme with outdated channel state information[J]. Journal of Beijing University of Posts and Telecommunications, 2014, 37(6): 44-48. doi: 10.13190/ j.jbupt.2014.06.009.

MOKHTAR M, GOMAA A, and AL-DHAHIR N. OFDM AF relaying under I/Q imbalance: Performance analysis and baseband compensation[J]. IEEE Transactions on Communications, 2013, 61(4): 1304-1313. doi: 10.1109/ TCOMM.2013.020813.120576.

MOKHTAR M, BOULOGEORGOS A, KARAGIANNIDIS G K, et al. OFDM opportunistic relaying under joint transmit/receive I/Q imbalance[J]. IEEE Transactions on Communications, 2014, 62(5): 1458-1468. doi: 10.1109/ TCOMM.2014.022314.130911.

MOKHTAR M, Al-DHAHIR N, and HAMILA R. I/Q imbalance and loop-back self-interference effects in full-duplex OFDM DF relays[C]. IEEE Wireless Communications and Networking Conference (WCNC), New Orleans, LA, USA, 2015: 81-86. doi: 10.1109/TCOMM.2014. 2325036.

MOKHTAR M, Al-DHAHIR N, and HAMILA R. OFDM full-duplex DF relaying under I/Q imbalance and loopback self-interference[J]. IEEE Transactions on Vehicular Technology, 2016, 65(8): 6737-6741. doi: 10.1109/TVT.2015. 2479257.

LI J, MATTHAIOU M and SVENSSON T. I/Q imbalance in two-way AF relaying[J]. IEEE Transactions on Communications, 2014, 62(7): 2271-2285. doi: 10.1109/ TCOMM.2014.2325036.

WANG Z, HUANG J, ZHOU S, et al. Iterative receiver processing for OFDM modulated physical-layer network coding in underwater acoustic channels[J]. IEEE Transactions on Communications, 2013, 61(2): 541-553. doi: 10.1109/TCOMM.2012.022513.120085.

期刊类型引用(3)

1.	安秀娟，王兆荣，张东坡. 非正交多址接入异构网络用户公平性问题研究. 现代电子技术. 2024(13): 61-65 . 百度学术
2.	谈苗苗，孔槐聪，黄庶沛，程铭，林敏. 基于上/下行NOMA的卫星通信系统遍历容量分析. 南京航空航天大学学报. 2023(05): 780-788 . 百度学术
3.	溪利亚，王舜燕. 无线网络中多址接入协议仿真测试. 计算机仿真. 2023(10): 249-253 . 百度学术

其他类型引用(1)

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

留言板