#### Abstract

A unified nonorthogonal multiple access (NOMA) framework for satellite-terrestrial network is investigated over Shadowed-Rician fading channels, which can be applied to both power domain NOMA (PD-NOMA) and code domain NOMA (CD-NOMA). Under the realistic assumption of hardware impairments (HIs) and imperfect successive interference cancellation (SIC), the achievable outage probability of the paired users is analyzed first, and the closed-form expressions of exact and asymptotic outage probability are derived. Next, the ergodic sum rate is investigated in different levels of HIs and residual interference. In addition, the energy efficiency in the delay-tolerant transmission mode is studied to characterize the performance of the satellite-terrestrial network. Finally, the Monte Carlo simulation results demonstrate the accuracy of theoretical derivation and showcase the effects of a wide range of satellite channel parameters and the realistic assumption on outage probability, ergodic rate, and energy efficiency.

#### 1. Introduction

The satellite-terrestrial network has been envisioned as one of the research focuses to provide global service to users in future wireless communication networks [1]. A number of recently published studies have demonstrated its superiority over traditional satellite communications. Nonorthogonal multiple access (NOMA) is an effective approach to provide massive connectivity and is considered a new paradigm for the multiple access technology [2–4]. Therefore, introducing NOMA technology into the satellite-terrestrial network can improve spectrum efficiency and resource utilization [5, 6].

The application of NOMA technology to satellite communications is a promising way to further expand the integration of space and earth. The NOMA-based hybrid satellite-terrestrial relay system was investigated, in which the outage probability, asymptotic outage probability, and achievable rate were derived in [7]. In order to analyze the multilayer heterogeneous satellite network, the authors in [8] proposed a NOMA-enabled two-layer geostationary earth orbit/low earth orbit (LEO) satellite network, where the ergodic capacity of the proposed system was investigated in the case of coexistence of frequencies. In [9], a dedicated half-duplex (HD) terrestrial decode-and-forward (DF) relay was used to assist a user with weak channel condition to complete communications with the satellite. Moreover, in [10], the authors extended the work of [9] to investigate an amplify-and-forward multiantenna relay scenario and performed the performance analyses for fixed gain and variable gain relaying protocols, in which the exact and asymptotic outage probability expressions were obtained. In order to guarantee the transmission delay and transmission rate for satellite industrial Internet of things, the authors in [11] proposed the satellite-terrestrial integrated industrial Internet of things and optimized the power distribution ratio of terrestrial users, thus achieving the purpose of reducing transmission cost. Considering the various delay quality of service (QoS) requirements for satellite Internet of things (S-IoT) applications, the authors in [12] investigated the delay-limited performance of NOMA-based downlink S-IoT networks by analyzing the effective capacity and proved the effects of different power distribution schemes on the proposed networks. In [13], the authors investigated the NOMA-assisted integrated satellite terrestrial networks with the partial relay selection scheme under imperfect successive interference cancellation (SIC), which can improve the spectral efficiency and expand the transmission coverage area. In [14], uplink and downlink NOMA were considered in intersatellite communication networks, where the optimum control methods were employed to delay suppression, reduce bit error rate, and increase throughput. The authors in [15] adopted an unmanned aerial vehicle using the DF protocol to assist satellite communicated with the terrestrial NOMA users; specifically, the outage probability and system throughput in delay-limited transmission mode were deduced to validate the superiority of the proposed system. In order to maximize the system capacity and energy efficiency, the NOMA architecture-enabled space-terrestrial satellite network was investigated, in which the satellite and terrestrial base stations provided services in a collaborative manner for users, and the NOMA scheme was only implemented for terrestrial networks [16]. More recently, a great deal of works have focused on the performance of hybrid satellite-terrestrial networks in cognitive radio (CR) scenarios [17–20]. In [17], the NOMA-based cognitive hybrid satellite-terrestrial overlay networks were proposed, in which the outage probability performance was investigated. One of the assumptions was that the power distribution scheme of NOMA users was determined by instantaneous channel conditions. The partial secondary network selection schemes were analyzed for NOMA-based overlay cognitive integrated satellite-terrestrial relay networks in [18], which proved the superiority of the proposed scheme. In [19], the authors derived the closed-form expression of outage probability for the primary satellite network, and the secondary transmitter was regarded as a relay to assist the primary transmission. Furthermore, the authors of [20] investigated the performance of the secondary network in terms of outage probability, sum throughput, and ergodic rate in which the nearby NOMA user operated in full-duplex mode and adopted DF relay protocol to improve the performance of the remote NOMA user.

Unfortunately, for some practical reasons such as quantization error, phase noise, in-phase/quadrature imbalance, and limitations on transceiver performance, the transmission nodes are vulnerable to hardware impairments (HIs) [21, 22]. However, the HIs from different sources have a certain randomness, resulting in these impairments which cannot be eliminated satisfactorily [23]. In [24–26], the authors proposed a realistic scenario in which the effects of HIs and imperfect channel state information (CSI) on the performance of simultaneous wireless information and power transfer NOMA relaying network were considered. The authors in [27] studied the outage performance of the NOMA-based integrated satellite-terrestrial relay network in the presence of HIs. In [28], the authors derived a closed-from expression for the outage probability in the NOMA-based satellite communication network considering HIs and channel estimation errors. The authors of [29] investigated the impacts of HIs on the secrecy performance of the NOMA-based integrated satellite multiple-terrestrial relay network.

However, existing NOMA schemes can be classified into two categories: power domain NOMA (PD-NOMA) and code domain power NOMA (CD-NOMA) according to the characteristics of resource mapping and spread spectrum mode. This division method is not conducive to the unification of the NOMA standard. Therefore, in order to facilitate the unification of NOMA technical solution, NOMA is divided into single-carrier NOMA and multicarrier NOMA according to the number of resource element (RE) or subcarrier occupied by users. CD-NOMA is a typical multicarrier NOMA, which uses a sparse matrix to map the data flow of multiple users to multiple REs or subcarriers. PD-NOMA can be regarded as a single carrier NOMA to realize user reuse by power adjustment. Therefore, from this perspective, PD-NOMA can be considered a special case of CD-NOMA, and CD-NOMA can be regarded as a special extension of PD-NOMA. While the aforementioned works have provided us with a better understanding of the application of PD-NOMA technology in satellite-terrestrial networks, the analysis for the unified NOMA framework-based satellite-terrestrial networks is far from well understood, which inspires the construction of this paper. In addition, from the perspective of practical application, it is of great significance to investigate the effects of HIs and imperfect SIC on satellite-terrestrial networks. In this paper, we propose a unified NOMA framework for satellite-terrestrial networks and analyze the performance of the proposed system with HIs and imperfect SIC under Shadowed-Rician fading channels. The essential contributions of this paper are summarized as follows: (1)We propose a unified NOMA technology framework for satellite-terrestrial networks, which integrates two design patterns, PD-NOMA and CD-NOMA. A sparse matrix is applied to map the data stream of users to REs, and the unified NOMA framework can be transformed into PD-NOMA scheme () and CD-NOMA scheme () according to the different numbers of RE occupied by users(2)Considering two imperfections, HIs and imperfect SIC, we investigate the outage performance of the proposed system and derive the closed-form expressions of exact outage probability for the paired NOMA users. More specifically, the asymptotic outage probability is provided to get more insights into the high SNR region. Additionally, the impacts of a wide range of satellite parameters on network performance are studied(3)We study the ergodic rate of the unified NOMA framework for the satellite-terrestrial network, and the effects of HIs and residual interference on ergodic rate are discussed. In addition, the energy efficiency in the delay-tolerant transmission mode is investigated to characterize the performance of the proposed system(4)Numerical results validate the theoretical analysis and the advantages of the satellite-terrestrial network adopting the NOMA scheme. The HIs and residual interference have significant negative effects on the outage performance, ergodic sum rate, and energy efficiency. Moreover, a pair of terrestrial users is able to achieve better system performance as the number of REs increases

The remainder of this article is organized as follows. In Section 2, the unified NOMA framework for satellite-terrestrial network is presented and the instantaneous SINR is derived. Section 3 provides the exact and asymptotic outage performance for a pair of terrestrial NOMA users. Then, in Section 4, the ergodic sum rate of satellite-terrestrial NOMA networks is analyzed. Section 5 obtains the numerical Monte Carlo results to verify the theoretical analysis. This paper is summarized in Section 6, and the mathematical proof is shown in Appendix.

The main mathematical symbols used in this paper are shown as follows: and denote the cumulative distribution function (CDF) and the probability density function (PDF) of a random variable , respectively; is the expectation operator. denotes the Euclidean two norms of a vector; is the dot product of two matrix elements.

#### 2. System Model

Consider a unified downlink NOMA scenario for the satellite-terrestrial network, where the satellite () transmits the superimposed signals to the terrestrial users directly or with the aid of the DF relay (), as depicted in Figure 1. In particular, maps the data flows of users into time/frequence REs by using a sparse spread spectrum matrix , where denotes the RE indication assigned to and is the -th column of the matrix . It is noteworthy that, depending on the number of RE occupied by users, the unified NOMA scheme can be converted into the PD-NOMA scheme () and CD-NOMA scheme (). The channel vector pertaining to the satellite link () is denoted by , where represents the channel response between the -th RE of and , and each channel response is independent and follows the Shadowed-Rician fading distribution, whereas the channel vector for the terrestrial link () is represented by , where represents the channel response between the -th RE of and with . Meanwhile, all the receiving nodes are disturbed by additive white Gaussian noise (AWGN) with zero mean and average power of . We assume that and terrestrial nodes are equipped with a single antenna and terrestrial users are divided into orthogonal paired user groups, where two prepaired users are ordered according to their respective channel gain [19, 30]. To facilitate analysis, it is assumed that and are randomly selected (user selection can further improve system performance, which inspires us to investigate more effective user selection strategies in future work) to execute the NOMA scheme [31], and the satellite link condition of is worse than that of due to the influence of shadow fading. In addition, due to the propagation characteristics and large transmission delay of the satellite channel, it is difficult to obtain perfect CSI for the satellite channel, which has been studied in [6, 32]. However, the main aim of this paper is to investigate the performance of the unified NOMA framework for satellite-terrestrial networks with HIs and imperfect SIC. Therefore, perfect CSI is assumed in this paper, and the results of this paper will be the basis for future research on imperfect CSI.

Two phases are used throughout the transmission process. In the first phase, sends the superposed information to and terrestrial users, where the destination user’s data flow is spread over a column of . In addition, is the transmit power at , , denotes the message transmitted to with unit power, and . and are the power allocation coefficients for and that satisfy and to guarantee the users’ fairness. , is the gain of the satellite antenna, denotes the beam gain of node , is the antenna gain of node , and denotes the included angle between and beam center in relation to the satellite [32]. , denotes the 3 dB angle, and is the first class-modified Bessel function with order . Moreover, is the free space loss between the satellite and node , where , , and denote the speed of light, operating frequency, and the distance between the satellite and node , respectively. Accordingly, the received signals at and terrestrial users can be given by where denotes the AWGN vector for the link between the satellite and node and is an identity matrix of size . is the HIs between and node , where is the level of HIs [29, 32], , and denotes the HI indication. Since maximal ratio combiner (MRC) is adopted at over REs, we have that , where represents the receive weight vector. Assuming and applying the SIC principle, the instantaneous received signal to interference and noise ratio (SINR) of decoding can be given by where . By applying the SIC scheme, begins decoding after decoding . However, taking into consideration the influence of imperfect SIC, the instantaneous received SINR of decoding can be represented as where represents the level of residual interference from the previous SIC process and (it is worth noting that we use the imperfect SIC model in [33, 34], assuming that is a constant. However, obeys the Gaussian distribution in some complicated cases [35, 36], which is beyond the scope of our research, but we will consider it in future work).

Similar to (2) and (3), we can express the instantaneous received SINRs via the - link and - link at respective nodes as where represents the level of residual interference from and .

In the second phase, transmits the superimposed signal to and with power and power allocation coefficient , where , , denotes the distances from to , and is the pathloss exponent. By spreading the terrestrial user’s data stream over a column of the matrix , the received signal at is represented as where , the distortion noise from to is denoted as , and is the AWGN at .

Assume that and ; the instantaneous SINR of decoding its own message can be given by

Following the SIC scheme, decodes firstly with the received SINR as

Subsequently, by subtracting the interfering term associated with from , the received SINR of decoding its own message can be given by where represents the level of residual interference from the preceding SIC cycle and .

#### 3. Channel Statistical Properties

Assuming that the channel link from satellite to the -th RE of terrestrial node is modeled as Shadowed-Rician fading distributions, the PDF and CDF of are given by [32]

with , , and , where , , and denote the average power of the multipath component, the shape parameter of Nakagami-, and the average power of the line-of-sight component, respectively. and denotes the Pochhammer symbol ([37], p. xliii). Additionally, the terrestrial link from relay to the -th RE of node undergoes identically distributed (i.i.d.) Rayleigh fading distribution.

Lemma 1. *Let be sum of squared identically distributed Shadowed-Rician random variables; the PDF of is given by
where and denotes the beta function ([37], Eq. (8.384.1)).*

*Proof. *See Appendix A.

Assume that all of the terrestrial links follow the Rayleigh fading distribution and , , and have the same column weights; we derive . With the aid of Eq. (7) of [38], the CDF and PDF of are given by
respectively.

#### 4. Outage Probability Analysis

##### 4.1. Exact Outage Probability

In this section, we derive the closed-form expressions for the exact outage probabilities. The outage event is defined as the instantaneous SINR is inferior to the predefined threshold , where and is the predefined rate for . Assume that the terrestrial users process the data received from and nodes by selection combining. Thus, the outage probability of can be formulated as

Theorem 2. *The expression for the outage probability of for is given by
**where , , , , , , , , , and . (17) is obtained in the conditions of , , , and .*

*Proof. *See Appendix B.

Corollary 3. *The outage probability expression for with is given by
*

*Proof. *By substituting (2)–(5), (9), and (10) into (16) and then by some algebraic application, we can obtain the outage probability expression for by means of (12). The proof is completed.

Next, we focus attention on calculating the outage probability of . falls into outage if it cannot decode its own signal at both phases or cannot decode ’s signal in the first phase. Hence, the outage probability of can be formulated as

Theorem 4. *The expression for the outage probability of for is given by
**where , , and . The above variables need to satisfy the conditions that and .*

*Proof. *See Appendix C.

Corollary 5. *The outage probability expression for with is given by
*

##### 4.2. Asymptotic Outage Probability

In order to reduce computational complexity and better understand the effects of key system indicators and channel parameters on the performance of the proposed system at a high SNR regime, the closed-form expressions of asymptotic outage probabilities for and are derived. As approaches infinity, the first summation term in equation (13) has a greater effect on the system performance than the other terms, which are approximately 0. Thus, we only retain the first summation term in equation (13) and employ the Maclaurin series expansion. The approximate expression of PDF for can be expressed as where is a higher-order infinitesimal of . Based on (22), the approximate expression CDF of for at the high SNR regime can be expressed as

In addition, we can easily acquire the approximate expression CDF of for at the high SNR regime which can be expressed as [32]

In the following, the approximate expression CDF of will be derived. By applying the Taylor series expansion, can be expanded as . Thereby, in (25) can be replaced by . By substituting into (25), when , the approximate expression CDF of for can be expressed as

For the special case of , when , with the approximation of can be expressed as where denotes the channel mean power.

By substituting (23)–(26) into (B.1) and (C.1), respectively, we can get the asymptotic outage probability of and at the high SNR regime as follows: respectively, where .

According to the conclusions of (27) and (28), the diversity orders of and in CD-NOMA () and PD-NOMA () are both /three.

#### 5. Ergodic Sum Rate and Energy Efficiency

##### 5.1. Ergodic Rate

The ergodic sum rate is an important indicator of system performance, which depends on the conditions of the channel links and the decoding process at all the terrestrial nodes. Hence, in this section, the ergodic sum rate of the proposed system is investigated.

The achievable rate of and can be given by respectively.

Thus, the ergodic rate of can be written as

Since it is difficult to obtain the exact expression of the ergodic sum rate, we derive the approximate ergodic sum rate.

Theorem 6. *The approximate expression for the ergodic rate of for is given by
**where and .*

*Proof. *See Appendix D.

Corollary 7. *The approximate expression for the ergodic rate of for is given by
**where and .*

In what follows, the ergodic rate of can be written as

By substituting (2), (4), (6), (8), and (9) into (34), after some mathematical operations, the approximate expression for the ergodic rate of for can be formulated as

In (35), and can be expressed as where .

Due to the proof process being similar to (32), we will not elaborate on it.

For the special case with , the approximate expression for the ergodic rate of is calculated as and in (38) can be expressed as respectively, where .

##### 5.2. Energy Efficiency

Energy efficiency is another important index to evaluate the system performance. According to the definition in [39], energy efficiency can be expressed as where is the fixed energy consumed by the transceiver and denotes the efficiency of the power amplifier. is the ergodic sum rate, which can be obtained from (32), (33), (35), and (38).

#### 6. Numerical Results

In this section, we provide numerical results to verify the analysis result and further evaluate the performance of the satellite-terrestrial networks. Assume that the type of satellite is low-orbit earth satellite, satellite links obey the Shadowed-Rician fading distribution, and channel parameters are shown in Table 1 [29, 40]. In this simulations, we set that the power allocation coefficients for and are and , the distances from to terrestrial users are and , and the pathloss exponent of the terrestrial link is . In addition, let the terrestrial channel mean power . The target rates are bit/s/Hz and bit/s/Hz [32], respectively. The Monte Carlo simulation parameters are depicted in Table 2 [5].

Figure 2 illustrates the effects of the transmit SNR on the outage probability with different numbers of REs, where the satellite links undergo FHS. In this example, we set and . For comparison purposes, we choose the PD-NOMA case () as the benchmark. It can be seen that the calculated result curves perfectly match well with the Monte Carlo simulation curves, and the asymptotic outage probabilities are approximated with the calculated results in the medium-high SNR region. The outage probabilities of and decrease as the value of transmit SNR increases. Additionally, the outage performance improves significantly when the number of RE increases from to and . The reason for this phenomenon is that CD-NOMA () has a higher diversity order than PD-NOMA, which is consistent with the analysis of (27) and (28). Thus, this important observation also indicates that the diversity gains of the paired NOMA users are closely related to the number of REs.

Figure 3 shows the impacts of the transmit SNR and HIs on the outage probabilities of paired users. In this simulation, we set , , and . We can observe that the exact analytical and asymptotic results are well aligned to the Monte Carlo results. From the outage probability versus SNR curves, it can also be seen that the outage performances for and increase as the transmit SNR increases. In particular, different levels of HIs are set to be from 0.01 to 0.03. As can be observed from the figure, the outage performances of paired NOMA users are becoming worse with the value of increasing from 0 to 0.03. Another important observation is that the outage probability of is more sensitive to HIs than that of .

Figure 4 depicts the outage probabilities of paired NOMA users versus the imperfect SIC and transmit SNR in three different satellite shadowing conditions. For comparison, the outage probabilities of paired users experiencing FHS fading are selected as the benchmark. It is observed that the outage performance for both NOMA users improves as the degree of shadow fading experienced by satellite links decreases. The reason for this phenomenon is that the paired users can achieve a better outage performance when they undergo a better shadowing scenario. Additionally, the effect of residual interference from the SIC process on the outage probability is focused, where satellite links are subject to FHS. In this simulation, we set and . Apparently, the outage performance of the will be significantly degraded if the level of residual interference increases by 0.03 from zero. Therefore, it is necessary to develop a reasonable SIC scheme in a practical scenario.

Figure 5 reveals the outage probability of the considered system against the transmit SNR and power allocation factor under the FHS scenario. We assume that the system performance is analyzed under ideal conditions () and set the transmit SNR at and to be 30 dB. As can be seen from Figure 5, the outage probability decreases for and increases for as the value of increases. However, the outage performance of starts to decline once goes beyond 0.4. This phenomenon indicates that the increase in , the corresponding power allocation factor , will decrease and will affect the decoding in the SIC process when . Further, there exists a value that enables NOMA systems to yield better outage performance. Consequently, to improve the outage performance and fairness of paired users, can be set between 0.17 and 0.33 when dB.

Figure 6 illustrates the ergodic sum rate versus the transmit SNR with different numbers of REs (i.e., , , and ). In this simulation, we set and , and all satellite links undergo FHS fading. A special case, (PD-NOMA), is selected as a benchmark to evaluate the ergodic rate of the proposed system. The curves of ergodic rate for and can be plotted according to (32), (33), (35), and (38). As can be observed from the figure, the approximate analysis results are close to the Monte Carlo simulation results. Apparently, as the number of REs increases, the ergodic rates of the paired NOMA users gradually increase, and the slope of the curve in the case of is steeper than that in the case of . This observation means that CD-NOMA () is capable of outperforming PD-NOMA on the ergodic sum rate. From different perspectives, the spectral efficiency of the proposed system is negligible in the low SNR regime, regardless of the value of .

In Figure 7, the ergodic sum rates of paired NOMA users versus the transmit SNR are presented in three different levels of HIs. In this example, we set , , and . As illustrated, the approximate analysis results are basically consistent with the numerical simulation. In order to facilitate comparison, the curves of ergodic sum rate under ideal conditions () is taken as a benchmark. As the level of HIs increases from 0 to 0.03, the ergodic rates of and are getting worse and exacerbating. In addition, it can be also observed that the ergodic rate of is more sensitive to HIs than that of , which is corresponding to the analysis results in Figure 3. Therefore, the design of the effective compensation algorithm and calibration method is of great significance to alleviate the influences of HIs on the performance of NOMA-based satellite-terrestrial networks in practical application.

Figure 8 illustrates the ergodic rate versus the transmit SNR with different levels of residual interference from the SIC cycle. In this simulation, we set the number of REs and the level of HIs . Then, the cases of residual interference level , , and are considered for the proposed system. As can be observed, the ergodic sum rate increases as the transmit SNR increases. Moreover, the ergodic rate of decreases significantly when the level of residual interference increases from 0 to 0.01, 0.02, and 0.03. Hence, it is necessary to design a reasonable multiuser reception algorithm to improve the performance of NOMA-based satellite-terrestrial networks in practical scenarios.

To illustrate the effects of different levels of HIs and different levels of residual interference and power distribution factors on the energy efficiency, Figure 9 depicts the energy efficiency versus the transmit SNR in delay-tolerant transmission mode, where and . The power amplifier and fixed energy loss are set as and W, respectively [6, 39]. As can be seen from Figure 9(a), when increases from 0 to 0.03, the energy efficiency of the proposed system decreases significantly, and the transmit SNR that enables the considered system to obtain the maximum energy efficiency moves to the right. This happens because the ergodic rate achieved by the terrestrial users decreases under the influence of HIs, which is consistent with the analysis results in Figure 7. In Figure 9(b), the energy efficiency under perfect SIC condition is obviously better than that under imperfect SIC condition, which also corresponds to the analysis in Figure 8. In addition, it is observed that the energy efficiency of the considered system is improved with the increase in the power distribution factor for . Although the energy efficiency can be improved to a certain extent by increasing the power allocation factor of , it will greatly sacrifice the performance of , thus losing the fairness of users in the weak channel. This phenomenon indicates that optimizing the power distribution factor can further enhance the performance of the paired NOMA users. Moreover, as illustrated in Figures 9(a) and 9(b), energy efficiency increases with the number of REs.

**(a)**

**(b)**

#### 7. Conclusion

This paper proposed a unified NOMA framework for the satellite-terrestrial network under Shadowed-Rician fading channels. Considering the realistic scenario with HIs and imperfect SIC, the performance of the considered system has been characterized in terms of outage probability, ergodic sum rate, and energy efficiency. More specifically, the closed-form expressions of exact outage probability for the paired users were derived. In order to gain more insights, the asymptotic outage probability was provided in the high SNR region. The analytical results showed that a pair of NOMA users was able to achieve better outage performance as the number of REs increased, and the CD-NOMA was capable of outperforming the PD-NOMA on the outage performance. Furthermore, the approximate expression of ergodic sum rate was obtained due to the integral being difficult to solve. We can observe that the HIs and residual interference had significant negative effects on ergodic sum rate, and the ergodic rate of is more sensitive to HIs than that of . The energy efficiency of the considered system with imperfect SIC and HIs has been discussed in delay-tolerant transmission mode. Finally, the Monte Carlo simulation was presented to verify the theoretical derivation. The unified NOMA framework applied to the satellite-terrestrial network provides a novel design idea, which is a relatively promising research topic.

#### Appendix

#### A. Proof of Lemma 1

For mathematical manageability, we assume that , , and have the same column weights. That means that , , and undergo the same distribution. Let ; can be represented as where the symbol represents the statistical convolution between independent distributions. And then we can use mathematical induction to prove the PDF of . When , we can get , and the corresponding PDF of is given by

Substituting (11) into (A.1) and utilizing Eq. (8.191.1) of [37], we can obtain

When , the corresponding PDF of is given by

Similarly, by plugging (A.3) and (11) into (A.4), can be obtained as follows:

Assuming that and , the PDF of is derived, as shown in (13).

The proof is completed.

#### B. Proof of Theorem 2

By assuming and and substituting (2)–(5), (9), and (10) into (16), the outage probability of can be expressed as

By plugging (13) into , while utilizing Eq. (3.351.2) of [37], can be evaluated as

Similar to the calculation of , the expression of can be expressed as follows:

To proceed forward, we focus on analyzing . With the aid of (25), can be evaluated as

By plugging (B.2)–(B.4) into (B.1), the closed-form expression of can be derived.

The proof is completed.

#### C. Proof of Theorem 4

By applying the assumptions in Appendix B and plugging (2), (6), and (8) into (19), the outage probability of can be expressed as

Using steps similar to Appendix B of Theorem 2, with the help of (13) and Eq. (3.351.2) of [37], and are evaluated as respectively. After some algebraic manipulations, can be expressed as

Plugging (C.2)–(C.4) into (C.1), we can derive (20) and complete the proof.

#### D. Proof of Theorem 6

By substituting (3), (5), and (11) into (31), the expression for the ergodic rate of for can be expressed as

Following the inequality ([21], Eq. (35)) , (D.1) can be approximated as

As a further development, we assume that and