Abstract
Due to the shortage of spectrum resources, a nonorthogonal multiple access- (NOMA-) based overlay cognitive integrated satellite-terrestrial relay network is established in this paper, where cognitive radio (CR) and NOMA scheme are utilized to enhance spectral efficiency. In this system, a secondary terrestrial transmitter assists the transmission of primary signal and transmits its own signal simultaneously. To maximize the achievable rate of a secondary network on the premise of ensuring the quality of service of primary network, the decoding order and power allocation are jointly optimized. Particularly, we first derive the optimal power allocation factor for a given decoding order. After searching all different decoding orders, the globally optimal joint design of the decoding order and power allocation is obtained. To gain deep insights, the exact and asymptotic expressions for outage probability of both primary and secondary networks are derived. Numerical results demonstrate the accuracy of our derivation and indicate the impacts of key parameters on the proposed system.
1. Introduction
Although the standard of the next-generation wireless network has not been unified, the introduction of satellite communication (SatCom) into the existing wireless communication network has become the consensus of experts and engineers [1, 2], which has the capabilities of all-weather operation, long transmission distance, and seamless coverage. In addition, the integrated satellite-terrestrial relay networks (ISTRNs) are becoming the promising architecture of the next-generation communication network and the Internet of things (IoT), which can offer comprehensive and convenient services [3, 4]. With the rapid development of ISTRN, the establishment of large-scale networks and the exponential growth of the number for terminals will cause serious shortage of spectrum resources. To solve the above problem, the application of cognitive radio (CR) and nonorthogonal multiple access (NOMA) into ISTRN have been regarded as effective approaches to improve spectrum utilization and deal with large-scale users [5, 6]. In general, CR enables unauthorized users to reasonably utilize the authorized spectrum [7, 8]. On the other hand, the NOMA technique allows different users to multiplex the same time/code/frequency [9, 10].
1.1. Related Works
Over recent years, the exploration of ISTRN has been investigated in various research areas. Based on fractional programming, a power allocation optimization scheme for downlink land mobile satellite (LMS) was proposed in [11], and the optimal average energy efficiency (EE) was derived. The authors of [12] investigated the performance of the LMS system, which considered the inherent hardware impairments (HIs) and jamming signal. In [13], the authors discussed the physical layer security of ISTRN, in which the secrecy outage probability (SOP) and average security capacity were derived under multiple illegal eavesdroppers and a threshold-based user scheduling scheme. A variety of methods to enhance the physical layer security and reliability of ISTRN were discussed in [14], and a kind of effective achievable rate and a new beamforming (BF) scheme were developed. The authors of [15] carried out joint BF and optimization for reconfigurable intelligent surface- (RIS-) aided ISTRN.
To deal with the problem of spectrum shortage caused by increasingly scarce spectrum resources, many existing works have considered the combination of CR and ISTRN. The authors of [16] investigated the cognitive ISTRN with multiple secondary networks, and the outage probability (OP) of the primary network was minimized by the proposed secondary network selection scheme. In [17], the authors proposed a secure BF scheme to deal with the situation of multiple eavesdroppers in rate-splitting multiple access-based cognitive ISTRN. In the case of full channel state information (CSI) and statistical CSI, the authors proposed two BF schemes to maximize the instantaneous rate of terrestrial users in cognitive ISTRN by using hybrid zero forcing and partial zero forcing methods in [18]. The authors of [19] utilized terrestrial base stations and cooperative terminals to enhance the secrecy performance of the cognitive ISTRN system. On this basis, a joint artificial noise produce and cooperative interference BF scheme was proposed.
For the sake of solving the service request of ever-increasing large-scale users, NOMA is considered a significant scheme in ISTRN [20]. Joint influence of channel estimation error and HIs on the security performance of cognitive integrated satellite-terrestrial relay networks (CISTRN) was investigated in [21]. The authors of [22] investigated the security performance of ISTRN with the NOMA scheme, where the HIs and multiple relays were considered and the SOP of the system was derived. In [23], the performance of ISTRN which considered cache and NOMA was studied, the OP and hit probability of the considered system were derived based on stochastic geometry. The performance of NOMA-based ISTRN with content delivery was discussed in [24], where the transmission delay can be reduced.
1.2. Motivations
The system performance of the NOMA-based system can be enhanced by designing the decoding order or power allocation scheme [25–27]. The authors of [25] minimized the overall OP under the cases of different decoding orders. In [26], to maximize the achievable rate of the secondary network with energy harvesting, the power allocation and time portion were optimized. In [27], the fair throughput for the NOMA-based system and orthogonal multiple access- (OMA-) based system was maximized. However, the above works cannot be simply utilized to optimize the performance of a priority-based NOMA system. In the traditional NOMA system, power is always allocated based on channel conditions. However, the priority-based NOMA system should guarantee the quality of service (QoS) of users with high priority. In fact, the QoS of high-priority users can be guaranteed and the system performance can be improved by optimizing the power allocation. On this basis, the optimization of the decoding order can make the system achieve optimal performance. A few existing works have investigated the joint optimization of the decoding order and power allocation [28–30], but the QoS of high-priority users was not strictly considered. Moreover, the authors of [31] considered the QoS of high-priority users, but the method was carried out in terrestrial networks, and there is no similar work in ISTRN.
Motivated by the above discussion, a NOMA-based overlay CISTRN is considered. As far as we know, no prior work has discussed the NOMA-based overlay CISTRN with joint optimization of the decoding order and power allocation.
1.3. Our Contributions
The main contributions of this paper can be summarized as follows: (1)Firstly, an overlay CISTRN with the NOMA scheme is established, in which the secondary network assists the transmission of the primary network in exchange for the opportunity to utilize the authorized spectrum of the primary network to transmit its own signals simultaneously. Besides, the satellite signal cannot reach the primary user and secondary user due to some severe shadowing(2)Secondly, to improve the achievable performance of the secondary network, the decoding order and power allocation are optimized with the aim of maximizing the achievable rate of the secondary network under the predefined constraints of the primary satellite network. Specifically, the closed-form expressions for the optimal decoding order and power allocation factors are derived, which is of great significance for the study of NOMA-based CISTRN(3)Finally, basing on the joint optimal design of the decoding order and power allocation, the exact OPs of both the primary network and secondary network are obtained. To get further insights, the asymptotic expressions and diversity order (DO) of the two networks are also derived. In addition, the numerical results validate our theoretical analysis. On this foundation, the impacts of many key parameters on the system performance are revealed
1.4. Paper Organization
The structure of the rest of this paper is as follows. Section 2 provides the system model and problem formulation. The joint decoding order and power allocation design are given in Section 3. The performance of the proposed system is discussed in Section 5. Section 4 gives the numerical results and reveals the impacts of key parameters on system performance. Finally, the full text summary is given in Section 6.
Notations. denotes the Gaussian distribution with mean and variance . is the expectation operator, and is the gamma function (Eq. 8.339.1 in [32]). The probability density function (PDF) and the cumulative distribution function (CDF) of variable are, respectively, represented by and . is the Pochhammer symbol (p. xliii in [32]). stands for the second kind modified Bessel function (Eq. 8.432.6 in [32]).
2. System Model and Problem Formulation
As shown in Figure 1, a NOMA-based overlay CISTRN is considered in our paper, which consists of a primary satellite source () and a secondary terrestrial transmitter () with their corresponding receivers ( and ). It is assumed that there is no direct satellite (DS) link between and due to severe shadowing effects, such as obstacles. Therefore, the secondary network is asked for assisting the transmission from to as a relay [6]. As a reward, utilizes the authorized spectrum of to communicate with . The NOMA scheme is adopted in to transmit the expected signals of and simultaneously, and the two-user pair is considered in our paper to balance the complexity and effectiveness. Each node works in half-duplex mode and applies an omnidirectional antenna to reduce the complexity of the system [16], and the decode-and-forward (DF) protocol is employed by . Specifically, the channel coefficients of - link, - link, and - link are denoted by , , and , respectively. Moreover, the satellite-terrestrial channel experiences Shadowed-Rician (SR) fading, and the terrestrial channels undergo Nakagami- fading. Without loss of generality, it is assumed that the channel coefficients are variable but remain constant in a time slot (furthermore, the perfect hardware impairments and decoding errors are assumed in our proposed system to simplify the analysis. Besides, the perfect CSI can be obtained by feedback and training from users). In addition, all receivers sustain additive white Gaussian noise (AWGN) with .

2.1. Problem Formulation
There are two phases in the entire transmission process, which is shown in Figure 2.

During the first phase, transmits its signal to , and the signal received by is represented by where is the transmit power of , denotes the satellite message with , and is AWGN. From (1), the achievable rate of at can be expressed as where and is the signal-to-noise ratio (SNR) of the satellite. The target rate of the primary network is denoted as . To decode the primary signal successfully, the following constraint must be satisfied, namely, .
In the second phase, integrated primary and its own signal with the superposition coding technique (SCT), and the signal received by and can be expressed by where , is the transmit power of , denotes the expected signal of with , is the power allocation factor with , and represents AWGN.
At the receiving ends, there are two decoding orders. The first is that decodes the desired signal directly, and performs successive interference cancellation (SIC). The second is that detects its expected signal, while adopts SIC. Besides, the achievable rates of at and at under the conditions of two decoding orders can be given as follows, respectively.
2.1.1. The First Decoding Order
In this order, decodes while is regarded as the serial interference. Then, the achievable rate of at can be expressed as where , is the SNR of the secondary transmitter, and the tag “1” indicates the first decoding order.
Recalling to the properties of the NOMA scheme, SIC is adopted at ; thus, is decoded by with as the cochannel interference (CCI); then, the decoded message is removed from . Finally, decodes from the remaining signal. Hence, the achievable rates of and at can be, respectively, represented as where .
2.1.2. The Second Decoding Order
In the second decoding order, SIC is applied at , and is treated as the serial interference at both and firstly. Therefore, the achievable rates of at and are expressed as, respectively, where the tag “2” indicates the second decoding order.
If can be decoded by successfully, the achievable rate of at can be represented as
2.2. Channel Model
In this section, the statistical characteristics of SR fading and Nakagami- fading are provided. Firstly, the PDF of is given by [33] where , , and . denotes the Nakagami- parameter, and and are the average power of the multipath component and that of the line-of-sight (LOS) component, respectively.
When is an integer, with the help of Eq. 07.20.03.009.01 and Eq. 07.02.03.0014.01 in [34], we can get
By taking (10) into (9), it can be rewritten as
Owing to , the PDF of can expressed as where and .
With the utilization of Eq. 3.351.2 in [32], we can get CDF of as
In addition, the PDF and CDF of are expressed as, respectively [35], where , is the fading severity parameter, and denotes the average fading power, in which is the distance from to or and denotes the path loss exponent.
3. Joint Decoding Order and Power Allocation Design
In this section, by jointly designing the decoding order and power allocation factor, we maximize the achievable rate of the secondary network under the predefined constraints of the primary network. Firstly, we find the optimal power allocation factor under the two decoding orders. Then, the optimal allocation factors are utilized to obtain the corresponding achievable rates. By comparing the calculated achievable rates under different decoding orders, the desired decoding order corresponding to the optimal achievable rate can be obtained.
3.1. The First Decoding Order
Under the first decoding order, the optimization problem can be formulated as follows.
Problem 1. where (17) guarantees that the primary signal can be decoded correctly and SIC can be performed successfully. Besides, (17) can be reexpressed as
where .
From (18), we can obtain with . It is worth noting that if the constraints of (17) and (18) cannot be satisfied, all power is allocated to transmit the message of the primary network.
It can be observed that is a decreasing function of ; thus, the optimal power allocation factor can be expressed as
3.2. The Second Decoding Order
Under the second decoding order, the optimization problem can be formulated as follows.
Problem 2. where is the target rate of the secondary network, (22) and (23) ensure that the primary network can decode the expected signal successfully, and we can get
where . Moreover, to ensure the feasible region of this problem, not an empty set, should be guaranteed. Hence, we can get .
From (21), we can obtain that is a strict decreasing function of . Thus, the optimal power allocation factor is expressed as
Only when can the system work under the second order. Thus, we just need to compare the optimal achievable rates under the first and second order with . After comparison, we find that when , . Therefore, the joint optimal design of the decoding order and power allocation is given by where represents the decoding order.
According to the above joint optimization, we can get the achievable rates of and at and as, respectively, where and .
Remark. The above derivations show that the user with low priority should be decoded firstly when the channel gain of the high-priority user is larger than that of the low-priority user, as well as the specific value. In great majority of existing works on the NOMA system, the high-priority users should be decoded firstly, which is diverse with our observation.
4. Performance Analysis
In this section, the exact and asymptotic expressions of OP for both the primary network and secondary network are obtained based on the derivations of achievable rates, where OP is defined as the probability of the achievable rate lower than the target rate. Besides, we derive the DOs to provide more insights.
4.1. Primary Network
The primary network will be interrupted if the achievable rate of at or is lower than ; thus, the OP can be expressed as
Then, by taking (13) and (15) into (30), we can get
To derive the asymptotic OP, we first give the asymptotic CDF of and in the high SNR regime as [35]
Finally, we substitute (32) and (33) into (30); the asymptotic expression for OP of the primary system can be expressed as
According to (34), the DO of the primary network can be expressed as .
4.2. Secondary Network
In the secondary network, OP means the probability of achievable rates of at lower than or that of at lower than , which can be given by
We can easily get
From (29), can be reexpressed as
To simplify the calculation, some mathematical processing is conducted for (37) as follows [31].
By merging (38a) and (38b), (38c) and (38d), and (38e) and (38f), we can get
Similar to the derivation of , with the utilization of (15), we can obtain
Then, the derivation of is carried on, which is given by
Let ; (46) can be rewritten as
By substituting (14) and (15) into , with the help of Eq. 1.111 and Eq. 3.471.9 in [32], the expression of can be derived as
Let ; can be rewritten as
By utilizing the Gaussian-Chebyshev quadrature [33], we can get where and is the Gaussian-Chebyshev approximation.
By taking similar derivations of and , and are given by, respectively,
Therefore, the expression of OP for the secondary network is given by
Similar to the asymptotic OP of the primary network, we can easily obtain the asymptotic OP of the secondary network as
In addition, the DO of the secondary network is .
5. Numerical Results
In this section, numerical results are provided to prove the correctness of our theoretical derivation. The simulation tool is MATLAB 2019a. Besides, the impacts of key parameters are revealed. Particularly, we set , , and [3]. Moreover, the system and channel parameters of SR fading are provided in Table 1 [12, 36].
Figure 3 illustrates that the achievable rates versus different target rates of the primary system with setting , , and . It can be observed that our theoretical analysis is consistent with the Monte Carlo (MC) simulations, which demonstrates the accuracy of the analytical results. We can find that the achievable rate of the primary network increases with and reaches a fixed value, namely, . Noting that when achieves , begins to increase. This is due to the fact that if the primary network cannot correctly decode its expected signal, a secondary transmitter allocates all power to the primary network. Besides, decreases with stricter target rates. In a high SNR regime, it can be clearly found that the achievable rate of the secondary network with our method is higher than that with the method in [6], which proves that our method is superior to the fixed power allocation method. Moreover, in a low SNR regime, the achievable rate of the secondary network is lower than that in [6]; it is because our method should ensure the QoS of the primary network.

Figure 4 depicts that OP versus different SR fading parameters with , , and . It can be observed that the primary network outperforms the secondary network. This is because we give higher priority to the primary network in the optimal design for ensuring its QoS. Moreover, the curves in Figure 3 display that lighter SR channel fading can lead to better system performance.

Figure 5 shows that OP versus different Nakagami- fading parameters with setting and . We can also observe that the performance of the primary network is superior to that of the secondary network. Besides, the larger Nakagami- fading parameters mean lighter channel fading, which can improve the system performance.

Figure 6 plots the OP versus different distances from to with and . It can be clearly seen that whether it is the primary network or the secondary network, the increase in the distance between the secondary transmitter and the primary user improves the OP of the system. It is not difficult to understand that longer distance degrades the communication quality.

However, in Figure 7, which illustrates the OP versus different distances from to with and , the distance between the secondary transmitter and the secondary user only affects the OP of the secondary network. It indicates that the OP of the primary network is not related to , while that of the secondary network will rise with the increase in both and .

Figure 8 depicts the OP versus different target rates with setting and . We can find that the larger target rate of the primary network will increase the OP of both the primary network and the secondary network. It can be explained by the fact that a larger target rate means higher communication requirements, which is more difficult to achieve. Besides, the stricter target rate of the secondary network only reduces the system performance of the secondary network. Like , has no effect on the performance of the primary network.

6. Conclusions
In this paper, we established an overlay CISTRN with the NOMA scheme, where the secondary network multiplexed spectrum resources by accessing the spectrum of the primary network. In exchange, the secondary network helped the information transmission of the primary network as a relay. To improve the system performance of the secondary network, we jointly optimized the decoding order and power allocation, where the QoS of the primary network has been ensured. In addition, we derived the exact OP of both the primary network and secondary network. To get further insights, the asymptotic OPs and DOs were provided. From the theoretical results, it can be observed that when the channel coefficient of the high-priority user is higher than that of the low-priority user as well as a fixed value, the system will decode the message of the low-priority user first, which will make the system have better performance. This conclusion is diverse from great majority of existing works on a NOMA-based system. Finally, our theoretical results were validated by simulation results, and the impacts of key system parameters on optimal design as well as system performance were revealed.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding the publication of this article.
Acknowledgments
This work was supported by the National Science Foundation of China under Grants 61901502 and 62001517, the National Postdoctoral Program for Innovative Talents under Grant BX20200101, and the Research Project of National University of Defense Technology (NUDT) under Grants ZK18-02-11 and 18-QNCXJ-029 and in part by the Research Project of Space Engineering University under Grants 2020XXAQ01 and 2019XXAQ05 and in part by the Science and Technology Innovation Cultivation Fund of Space Engineering University. Thanks are due to all the authors for their contribution to this paper.