Next Article in Journal
A Smart Visual Sensing Concept Involving Deep Learning for a Robust Optical Character Recognition under Hard Real-World Conditions
Next Article in Special Issue
FDSS-Based DFT-s-OFDM for 6G Wireless Sensing
Previous Article in Journal
Modeling and Control of a Spherical Robot in the CoppeliaSim Simulator
Previous Article in Special Issue
Guaranteeing QoS for NOMA-Enabled URLLC Based on κμ Shadowed Fading Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Resource Allocation in Multi-Carrier Multiplexed NOMA Cooperative System

1
Beijing National Research Center for Information Science and Technology, Tsinghua University, Beijing 100084, China
2
Department of Electronic Engineering, Tsinghua University, Beijing 100084, China
3
State Key Laboratory of Wireless Mobile Communications, China Academy of Telecommunication Technology, Beijing 100084, China
4
School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
5
The National Key Laboratory of Science and Technology on Communications, University of Electronic Science and Technology of China, Chengdu 611731, China
6
School of Information and Communication Engineering, Beijing University of Posts and Telecommunications, Beijing 100876, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(16), 6023; https://doi.org/10.3390/s22166023
Submission received: 29 June 2022 / Revised: 28 July 2022 / Accepted: 10 August 2022 / Published: 12 August 2022

Abstract

:
Non-orthogonal multiple access (NOMA) cooperative communication technology can combine the advantages of NOMA and cooperative communication, providing high spectrum efficiency and increasing user coverage for next-generation wireless systems. However, the research on NOMA cooperative communication technology is still in a preliminary stage and has mainly concentrated on the scenario of fewer users. This paper focuses on a user-centered NOMA collaboration system in an ultra-dense network, and it constructs a resource allocation optimization problem to meet the demands of each user. Then, this paper decomposes the optimization problem into two subproblems; one is the grouping match among multiple relays and users, and the other is jointly allocating power and subcarrier resources. Accordingly, a dynamic packet matching algorithm based on Gale–Shapley and an iterative algorithm based on the difference of convex functions programing are proposed. Compared with existing schemes, the proposed algorithms can improve system throughput while ensuring the quality of service of users.

1. Introduction

It is predictable that with the expansion in the Internet of Things (IoT) and the development of communication, a large number of wireless connections and huge data traffic will pose challenges for next-generation wireless systems [1]. The demand for spectrum efficiency and network capacity has grown rapidly [2]; since orthogonal resources are limited, traditional orthogonal multiple access (OMA) has difficulty meeting the multiple UEs demands [3,4]. Recently, due to the superior spectral efficiency, non-orthogonal multiple access (NOMA) has attracted tremendous attention in industry and academia [5]. Compared with traditional OAM, NOMA can reuse non-orthogonal superposition and then assume large-scale connectivity [6,7].
Previous research on single-carrier NOMA technology is relatively complete, and performance evaluations at the link and system levels have proven that NOMA has a better transmission rate and lower error rate than OMA systems [8]. For multi-carrier NOMA technology, the research is not sufficient and has mainly focused on the typical sparse code multiple access (SCMA) and pattern division multiple access (PDMA) technology. Research on SCMA technology mainly focuses on codebook design, channel transmission rate assessment and receiver design [9,10,11,12]. Research on PDMA technology has primarily focused on pattern design, interruption probability analysis and receiver design [13,14,15].
Cooperative communication is an important method for combating channel fading and path loss and for reducing the shadow effect [16,17]. It was first introduced in [18] and became an important component in the previous communication network. With the distributed transmission of relay at the transmitting site and the combination of signals at the receiving site, cooperative communication technology can obtain a cooperative diversity gain similar to the spatial diversity gain of the multiple input multiple output (MIMO) [19]. Because of the above feature, cooperative communication is widely used in IoT scenarios where the size of massive UEs is small and difficult to implement multiple antennas [20,21].
NOMA cooperative communication technology was introduced in [22]. It can combine the advantages of NOMA and cooperative communication to meet the demands of next-generation communication. However, there have been few studies addressing the combination of the two technologies, especially on NOMA multi-carrier cooperative communication technology.
Multi-carrier multiplexed NOMA cooperative technology can effectively obtain multiplexing diversity gain and improve system performance [23]. Luo et al. [24] studied the resource allocation optimization of an SCMA cooperative system, and the optimization was carried out with the weighted sum of power, codebook and subcarrier pairing as the alternative objective. Han et al. [25] studied the method of self-interference mitigation for an SCMA cooperative system with a large-scale transceiver antenna. Tang et al. [26] proposed an uplink PDMA collaboration system with a half-duplex decoder and relay, and they analyzed the interrupt performance. Tang et al. [27] analyzed the outage probability of a downlink PDMA collaboration system with decoding and forwarding of a half-duplex relay. However, in [26,27], the derived closed-form expression of interrupt probability was simplified to only the scenario with three users and two carriers while not considering the return link discontinuity. The more general interrupt probability expression and the full-duplex scenario still need to be studied and analyzed. In addition, Sun et al. [28] studied the multi-carrier NOMA cooperative system, which adopts a full-duplex base station for upstream and downstream simultaneous transmission, but the resource allocation optimization algorithm also assumes that the stack number of upstream or downstream users on each subcarrier does not exceed 2.
There are still several problems in the previous research studies. First, the theoretical boundaries of capacity and outage probability in a general scenario are not clear. Second, scheduling and resource allocation under multi-user and multi-relay with multi-carrier are extended to larger dimensions. Furthermore, effective multi-dimensional constellations, such as those in SCMA and PDMA, are more difficult to design and optimize than single-dimensional power segmentation. Finally, the widespread deployment of multi-carrier NOMA in conjunction with existing orthogonal frequency division multiple access lacks viable applications and validation.
In this paper, we studied multi-carrier multiplexed NOMA cooperative technology in a super-dense network, and the main contributions of this paper are listed below.
  • We designed a user-centered multi-carrier multiplexed NOMA cooperative system that can fully combine the advantages of multi-carrier NOMA and cooperative communication technology to meet the abundant UE demands.
  • We constructed a problem to optimize throughput while ensuring multiple users’ demands and decomposed it into two subproblems. Then, we proposed the corresponding dynamic grouping matching algorithm and iterative algorithm based on the difference of convex functions programing (DCP) to solve them.
  • Simulations were used to verify the effectiveness of the proposed NOMA cooperative network framework and the corresponding algorithms. Compared with two existing schemes, combining the dynamic grouping matching algorithm with an iterative algorithm improved system throughput while ensuring user quality of service (QoS).
The remainder of the paper is organized as follows. Section 2 describes the system model, including the signaling model and the throughput model. In Section 3, the problem is formulated as an optimization problem. In Section 4, we propose the resource allocation algorithm of the cooperative network. Section 5 presents the simulation results, which prove the effectiveness of the proposed algorithm. Finally, the conclusions are drawn in Section 6.

2. System Model

In this section, we describe the downlink NOMA-based cooperative network setting [29] consisting of a base station (BS), M relays, and N UEs, as shown in Figure 1. Each node is equipped with a transmit antenna and a receive antenna. The system frequency band is divided into K subcarriers. The signals of different UEs or different packets of a UE can be superposed in one subcarrier to transmit simultaneously. In addition, all the relays are connected to the BS in the backhaul stage. We assume that the UEs and relays follow two independent Poisson point processes (PPPs) with the densities of λ u and λ r , respectively. The notions about the system model in this section are listed in Table 1.
The BS can support the backhaul stage for relays and provide access services for UEs. In this paper, the UEs are assumed to be served dynamically via dense relays. The BS and the relays share the same frequency band, and the relays work in the time division duplex (TDD) mode. The signals passed into the backhaul stage and forward stage do not affect each other. At the backhaul stage, the signal is transmitted from the BS to the relays, and then, the relays decode the information of users and transmit it to corresponding users during the forward stage. Therefore, the downlink transmission can be divided into two processes, backhaul transmission and forward transmission, as illustrated in Figure 1.

2.1. Signaling Model

In the forward stage, we take the nth UE as an example to illustrate the signals in the downlink cooperative network, n { 1 , 2 , , N } . The inter-group interference can be avoided by properly grouping relays and allocating subcarriers [30]. We use x m n k to denote the signal, transmitted from the mth relay to the nth UE on the kth subcarrier. The transit power of x m n k is p m n k . Meanwhile, h m n k and ϖ m n are the small-scale and large-scale fading channel coefficients from the mth relay to the nth UE on the kth subcarrier, respectively. The channels of the forward and backhaul stages are independent Rayleigh fading channels, and the path loss exponent is α . ϖ m n = d m n α , where d m n denotes the distance from the mth relay to the nth UE. Then, the signal received on the kth subcarrier of the nth UE can be written as
y n A , k = p m n k ϖ m n h m n k x m n k + l m M p l n k ϖ m n h l n k x l n k + z n k ,
where z n k is the AWGN at the receiver of the nth UE on the kth subcarrier with mean zero and variance σ 2 .
In the backhaul stage, we take the mth relay ( m { 1 , 2 , , M } ) as an example to analyze the backhaul stage signals in the considered downlink cooperative network. Suppose that the signal transmitted from the BS to the mth relay on the kth subcarrier is x m k , and the power of the signal is q m k . Additionally, h m k and ϖ m denote the small-scale and large-scale channel coefficients from the BS to the mth relay on the kth subcarrier. Then, the received signal on the kth subcarrier of the mth relay can be written as
y m B , k = q m k ϖ m h m k x m k + h m k l m M q l k ϖ l x l k interference of other relays + z m k ,
where z m k denotes the additive white Gaussian noise (AWGN) with mean zero and variance σ 2 at the receiver of the mth relay on the kth subcarrier.

2.2. Throughput Model

Since the transmission process is divided into two stages, the throughput is analyzed separately at the two stages.
In the forward stage, since the successive interference cancellation (SIC) technique is applied at the receivers to decode the signals from different relays, we assume that the channel coefficients meet | H 1 n k | | H 2 n k | | H M n k | . H m n k = ϖ m h m n k represents the channel coefficient between the nth UE and the mth relay on the kth subcarrier. Then, the decoding order is consistent with the relay indexes. Thus, the received signal-to-interference-plus-noise ratio (SINR) of the nth UE served by the mth relay on the kth subcarrier is written as
S I N R m n A , k = ϖ m n p m n k | h m n k | 2 σ 2 + l = m + 1 M | h l n k | 2 ϖ l n p l n k .
The corresponding throughput of the nth UE is given by
R ˜ n A = k = 1 K m = 1 M c m n k β k log 2 ( 1 + S I N R m n A , k ) .
In the backhaul stage, without loss of generality, we assume that the channel coefficients yield | H 1 k | | H 2 k | | H M k | . Here, H m k = ϖ m h m k represents the channel coefficient between the BS and the mth relay on the kth subcarrier. Then, the decoding is carried out in the reverse order of the relay indexes. The SINR of the mth relay on the kth subcarrier is given by
S I N R m B , k = ϖ m q m k | h m k | 2 σ 2 + ϖ m | h m k | 2 l = 1 m 1 q l k .
Correspondingly, the backhaul throughput for the nth UE can be given by
R ˜ n B = k = 1 K m = 1 M c m n k ( 1 β k ) log 2 ( 1 + S I N R m B , k ) ,
where c m n k indicates whether the nth UE is served by the mth relay on the kth subcarrier or not, and c m n k { 0 , 1 } . If c m n k =1, the nth UE is served by the mth relay on the kth subcarrier. β k ( 0 β k 1 ) denotes the proportion of a time slot occupied by the forward stage on the kth subcarrier, and 1 β k is the proportion of the time slot used for backhaul.
Since the signal needs to be transmitted to the relays first and then forwarded to the UEs, the system throughput is given by
R ˜ = n = 1 N R ˜ n A , R ˜ n A R ˜ n B .

3. Problem Formulation

In this section, we maximize the system throughput under QoS constraints. Because the SIC technique is applied at the receiver, the complexity of the receiver grows with the number of superposed signals on a subcarrier. We assume that each UE can be served by up to Q relays on a subcarrier to harness the complexity of the receiver at the UEs. The constraints are written as
C 1 ( a ) : c m n k { 0 , 1 } , m , n , k ( b ) : m = 1 M c m n k Q , n , k .
Apart from the overall power constraints of the system, the power allocation in the NOMA system needs to satisfy the threshold for SIC decoding at the receiver (cf., the OMA system). Therefore, constraints C2 and C3 are given by
C 2 ( a ) : k = 1 K p m n k P m , max , m ( b ) : m = 1 M k = 1 K q m k P max BS ,
and
C 3 ( a ) : ϖ m n p m n k | h m n k | 2 q = m + 1 M ϖ q n | h q n k | 2 c q n k p q n k p thr , m , n , k ( b ) : ϖ m | h m k | 2 ( q m k l = 1 m 1 q l k ) p thr , m , k ,
where P max BS and P m , max denote the maximum available powers of the BS and the mth relay, respectively, and p thr is the decoding power threshold for the SIC receiver.
In terms of the QoS of each UE, we consider
C 4 : R ^ n target R ˜ n A R ˜ n B , n ,
where R ^ n target denotes the target data rate of the nth UE. Additionally, the time slot assignment coefficient between the access stage and backhaul stage on an arbitrary kth subcarrier needs to meet
C 5 : 0 β k 1 , k .
Therefore, the optimization problem can be formulated as
maximize c , fi , p , q R ˜ ( c , fi , p , q ) = n = 1 N R ˜ n A s . t . C 1 , C 2 , C 3 , C 4 , C 5 ,
where q R + K M × 1 and p R + N K M × 1 collect the power q m k allocated on the BS and the power p m n k allocated on the relays, respectively. c Z N K M × 1 and fi R K × 1 collect the variables c m n k and β k , respectively.

4. Resource Allocation Algorithms

Problem (13) is a mixed integer non-linear programming problem. It is challenging to derive a global-optimal solution [31]. In this paper, a low-complexity suboptimal solution is developed in the presence of multiple relays and UEs. Problem (13) is divided into two subproblems. First, we apply a dynamic group matching algorithm to map each UE with relays. Then, an iterative algorithm is proposed based on the D.C. programming to achieve a suboptimal solution for the joint power and subcarrier allocation.

4.1. Dynamic Group Matching for UEs and Relays

The grouping process of relays and UEs is a matching process between each UE and a set of relays serving the UE. To maximize the system throughput, we apply a deferred-acceptance strategy from the Gale–Shapley algorithm to balance the two-side matching priority of the UEs and relays. Let Φ ( m , n ) represent the matched pair of the mth relay and the nth UE, and let Φ denote the set of matched pairs. | Φ ( m , n ) | = 1 denotes that the nth UE is matched with the mth relay; otherwise, | Φ ( m , n ) | = 0 . We define an evaluation model of the pair between the nth UE and the mth relay as
R m n = log 2 ( 1 + S I N R m n A ) ,
where
S I N R m n A = ϖ m n | h m n k | 2 σ 2 + l = m + 1 M ϖ l n | h l n k | 2 .
With a two-sided competitive selection of the UEs and relays, each node has its matching priority list to match with others. We denote the matching priority sets of UEs and relays as
{ M P _ U E } = { M P _ U E 1 , , M P _ U E n , , M P _ U E N } ,
{ M P _ R E } = { M P _ R E 1 , , M P _ R E m , , M P _ R E M } ,
where M P _ U E n is the matching priority list that the nth UE matches with its nearby relays; similarly, M P _ R E m is the matching priority list of the nearby UEs that the mth relay can match with. They can be further represented as
M P _ U E n = { M P _ U E n ( 1 ) , , M P _ U E n ( m n ) , M P _ U E n ( M n ) } ,
M P _ R E m = { M P _ R E m ( 1 ) , , M P _ R E m ( n m ) , M P _ R E m ( N m ) } ,
where M n and N m are the number of relays near the nth UE and the number of UEs near the mth relay, respectively; M P _ U E n ( m n ) denotes the relay whose matching priority of the nth UE is m n , and M P _ R E m ( n m ) denotes the UE whose matching priority of the mth relay is n m . If M P _ U E n ( m n ) > M P _ U E n ( l n ) , it signifies that the matching priority of the nth UE with the m n th relay is higher than the matching priority of the nth UE with the l n th relay. We also define the relay with the highest matching priority of the nth UE as M P _ U E n highest . Correspondingly, we define the UE with the highest matching priority of the mth relay as M P _ R E m highest . In this paper, to maximize system throughput, we have
M P _ U E n highest = arg m max m M n h m n ,
and
M P _ R E m highest = arg n max R m n .
The reason for our choice of the throughput R m n and small-scale channel coefficient h m n as the priority judgment criteria of relays and UEs is that they are our optimization function or one of the parameters of the optimization function, and the results screened by these criteria are more conducive to the maximization of throughput.
With the above illustration, the dynamic grouping matching algorithm between UEs and relays can be described as follows. First, we initialize the matching priority according to the available CSI. Then, we divide the grouping process into two matching processes. The first process is to guarantee that each UE can be served by a relay, and the second process is to group the relays for each UE.
In the first process, each UE requests matching the relay that prioritizes the UE over the other UEs. Then, each relay that has received the matching request from the UEs matches the UE which prioritizes the relay over the other relays, and then, it rejects the other UEs. This process is repeated until all UEs are served by at least one relay.
In the second process, each UE requests matching the unmatched relay that has the highest priority to the UE. Subsequently, the relay that has received matching requests from UEs selects the UE according to its matching priority if the number of relays in a group is below Q. When the number of relays in a group is Q, we determine whether the UE sending this matching request is more effective for improving the throughput than the other UEs in the group.
If this is the case, then we update the matched pair; otherwise, we reject the matching request. This process is repeated until all the relays are grouped or no UEs request matching with any relays. The details of the dynamic grouping matching algorithm are provided in Algorithm 1.
Algorithm 1 Dynamic Group Matching Algorithm.
1:
Initialization: Initialize the matched pairing set of the UEs and relays Φ = , the unmatched set of the UEs and relays U _ U E , U _ R E , and initialize the matching priority sets of the UEs and relays { M P _ U E } , { M P _ R E } through the Equations (20) and (21)
2:
while  U _ U E  do
3:   
Each UE in U _ U E requests to match its highest matching priority relay from U _ R E according to the matching priority set { M P _ U E }
4:   
for relay m = 1, 2,…, M do
5:     
Each relay matches the UE with the highest priority according to the matching priority set of relays { M P _ R E } and rejects the other UEs
6:     
The rejected UEs remove the mth relay from its matching priority set { M P _ U E }
7:     
Add the matched paring Φ ( m , n ) to the set Φ and remove the mth relay and nth UE from U _ R E and U _ U E , respectively
8:   
end for
9:
end while
10:
while { M P _ U E } or U _ R E  do
11:   
Each UE requests to match its highest matching priority relay from U _ R E according to the updated set { M P _ U E }
12:   
for relay m = 1, 2,…, M do
13:     
The mth relay makes the following judgment for its highest matching priority UE according to its matching priority set { M P _ R E }
14:     
if  m = 1 M | Φ ( m , n ) | < Q  then
15:        
Add the matched paring Φ ( m , n ) to the set Φ and remove the mth relay and nth UE from U _ R E and U _ U E , respectively
16:     
else
17:        
The relay matches with the nth UE when there exists a relay that satisfies ψ m n > ψ l n and | Φ ( l , n ) | = 1 ; then, it updates Φ and removes the lth relay into U _ R E
18:        
Otherwise, reject the matching request of the nth UE and remove the mth relay from its matching priority set { M P _ U E }
19:
     end if
20:
   end for
21:
end while

4.2. Joint Power and Subcarrier Allocation Algorithm for the Cooperative Network

Given the matching outcome described in Section IV, we propose the joint power and subcarrier allocation algorithm based on the D.C. programming to optimize the system throughput of the cooperative network. We denote the assignment c m n k as c m n k = | Φ ( m , n ) | b n k , where | Φ ( m , n ) | { 0 , 1 } denotes the matched pair of the nth UE and the mth relay, and b n k { 0 , 1 } denotes whether the nth UE is served on the kth subcarrier or not. Problem (13) can be rewritten as
maximize b , β , p , q R ˜ ( b , β , p , q ) s . t . C 1 , C 2 , C 3 , C 4 , C 5 ,
where b Z N K × 1 collects variables b n k , n , k .
We combine the mixed integer constraint C 1 with constraint C 5 , as given by
C 1 ( a ) : 0 b n k β k = u n k 1 , n , k ( b ) : 0 b n k ( 1 β k ) = v n k 1 , n , k ( c ) : u n k + v n k { 0 , 1 } , n , k ,
The matching between the UEs and relays in constraint C1 is obtained by Algorithm 1. Only b n k remains to be solved in constraint C1. The integer constraint C 1 c is equivalent to the following expression:
C 6 ( a ) : n = 1 N k = 1 K ( ( u n k + v n k ) ( u n k + v n k ) 2 ) 0 , ( b ) : 0 u n k + v n k 1 , n , k .
Now, the optimization with the integer constraints is transformed to a continuous-value problem. We define u R N K × 1 , and v R N K × 1 to collect the variables u n k and v n k , respectively. Problem (22) can be reformulated as:
minimize u , v , p , q R ˜ ( u , p ) s . t . C 1 a , C 1 b , C 2 , C 3 , C 4 , C 6 .
According to the theorem of monotone optimization [28], the equivalent problem of (25) can be formed as:
minimize u , v , p , q R ˜ ( u , p ) + η { n = 1 N k = 1 K ( ( u n k + v n k ) ( u n k + v n k ) 2 ) } , s . t . C 1 a , C 1 b , C 2 , C 3 , C 4 , C 6 b .
where η is a sufficiently large penalty factor if u n k + v n k is neither 0 nor 1, and η 1 .
Then, we transform the decoding threshold constraint C3 into a maximum interference [32] constraint C3’ by
C 3 ( a ) : l = m + 1 M ϖ l n | h l n k | 2 c l n k p l n k ξ a , m , n , k ( b ) : ϖ m | h m k | 2 l = 1 m 1 q l k ξ l , m , k .
The new constraint C3′ is a convex set. However, the problem is still a non-convex problem, since neither the objective function nor constraint C4 is convex. Nevertheless, the following equivalent form always holds,
log 2 ( 1 + S I N R m n A , k ) = log 2 ( σ 2 + l = m M ϖ l n | h l n k | 2 p l n k ) log 2 ( σ 2 + l = m + 1 M ϖ l n | h l n k | 2 p l n k ) f m n A , k ( p ) g m n A , k ( p ) .
Therefore, we derive that
R ˜ n A = k = 1 K m = 1 M u n k ( f m n A , k ( p ) g m n A , k ( p ) ) F n A ( u , p ) G n A ( u , p ) .
Similarly, we have
R ˜ n B = k = 1 K m = 1 M v n k ( f m B , k ( q ) g m B , k ( q ) ) F n B ( v , q ) G n B ( v , q ) ,
where
f m B , k ( q ) = log 2 ( σ 2 + ϖ m | h m k | 2 l = 1 m q l k ) , and
g m B , k ( q ) = log 2 ( σ 2 + ϖ m | h m k | 2 l = 1 m 1 q l k ) .
Then, the non-convex constraint C4 can be rewritten as
C 4 ( a ) : F n A ( u , p ) + G n B ( v , q ) ( G n A ( u , p ) + F n B ( v , q ) ) 0 ( b ) : R ^ n , target + G n A ( u , p ) F n A ( u , p ) 0 .
Constraint C4′ is the difference of two convex Functions (31)–(33). Additionally, we have
R ˜ ( u , p ) = n = 1 N F n A ( u , p ) n = 1 N G n A ( u , p ) F A ( u , p ) G A ( u , p ) .
Therefore, we can rewrite (26) as
minimize u , v , p , q F A ( u , p ) G A ( u , p ) + η ( H ( u , v ) M ( u , v ) ) s . t . C 1 a , C 1 b , C 2 , C 3 , C 4 , C 6 b ,
where H ( u , v ) = n = 1 N k = 1 K ( u n k + v n k ) , M ( u , v ) = n = 1 N k = 1 K ( u n k + v n k ) 2 .
Note that F n A ( u , p ) , G n A ( u , p ) , F n B ( v , q ) , G n B ( v , q ) , H ( u , v ) , and M ( u , v ) are convex functions. Therefore, problem (35) is a D.C. program. We can implement successive convex approximation to obtain a suboptimal solution of the problem [33,34]. Given the differentiability of the convex functions F n A ( u , p ) , G n A ( u , p ) , F n B ( v , q ) , and M ( u , v ) , for any feasible point u ( τ ) , v ( τ ) , p ( τ ) , and q ( τ ) , we have
F n A ( u , p ) F n A ( u ( τ ) , p ( τ ) ) + u F n A ( u ( τ ) , p ( τ ) ) ( u u ( τ ) ) + p F n A ( u ( τ ) , p ( τ ) ) ( p p ( τ ) ) u ( τ ) , p ( τ ) affine F n A ( u , p ) ,
G n A ( u , p ) G n A ( u ( τ ) , p ( τ ) ) + u G n A ( u ( τ ) , p ( τ ) ) ( u u ( τ ) ) + p G n A ( u ( τ ) , p ( τ ) ) ( p p ( τ ) ) u ( τ ) , p ( τ ) affine G n A ( u , p ) ,
F n B ( v , q ) F n B ( v ( τ ) , q ( τ ) ) + v F n B ( v ( τ ) , q ( τ ) ) ( v v ( τ ) ) + q F n B ( v ( τ ) , q ( τ ) ) ( q q ( τ ) ) v ( τ ) , q ( τ ) affine F n B ( v , q ) ,
and
M ( u , v ) M ( u ( τ ) , v ( τ ) ) + u M ( u ( τ ) , v ( τ ) ) ( u u ( τ ) ) + v M ( u ( τ ) , v ( τ ) ) ( v v ( τ ) ) u ( τ ) , v ( τ ) affine M ( u , v ) .
In (36)–(39), u ( τ ) , p ( τ ) affine F n A ( u , p ) , u ( τ ) , p ( τ ) affine G n A ( u , p ) , v ( τ ) , q ( τ ) affine F n B ( v , q ) , and u ( τ ) , v ( τ ) affine M ( u , v ) are affine functions of F n A ( u , p ) , G n A ( u , p ) , F n B ( v , q ) , and M ( u , v ) , respectively. The gradients in the affine functions can be given by
u F n A ( u ( τ ) , p ( τ ) ) ( u u ( τ ) ) = k = 1 K m = 1 M f m n A , k ( p ( τ ) ) ( u n k u n k ( τ ) ) ,
p F n A ( u ( τ ) , p ( τ ) ) ( p p ( τ ) ) = k = 1 K m = 1 M u n k ( τ ) q = m M ϖ q n | h q n k | 2 ( p q n k p q n k ( τ ) ) ( σ 2 + q = m M ϖ q n | h q n k | 2 p q n k ( τ ) ) ln 2 ,
u G n A ( u ( τ ) , p ( τ ) ) ( u u ( τ ) ) = k = 1 K m = 1 M g m n A , k ( p ( τ ) ) ( u n k u n k ( τ ) ) ,
p G n A ( u ( τ ) , p ( τ ) ) ( p p ( τ ) ) = k = 1 K m = 1 M u n k ( τ ) q = m + 1 M ϖ q n | h q n k | 2 ( p q n k p q n k ( τ ) ) ( σ 2 + q = m + 1 M ϖ q n | h q n k | 2 p q n k ( τ ) ) ln 2 ,
v F n B ( v ( τ ) , q ( τ ) ) ( v v ( τ ) ) = k = 1 K m = 1 M f m B , k ( q ( τ ) ) ( v n k v n k ( τ ) ) ,
q F n B ( v ( τ ) , q ( τ ) ) ( q q ( τ ) ) = k = 1 K m = 1 M v n k ( τ ) ϖ m | h m k | 2 q = 1 m ( p q k p q k ( τ ) ) ( σ 2 + ϖ m | h m k | 2 q = 1 m p q k ( τ ) ) ln 2 ,
u M ( u ( τ ) , v ( τ ) ) ( u u ( τ ) ) = n = 1 N k = 1 K 2 ( u n k ( τ ) + v n k ( τ ) ) ( u n k u n k ( τ ) ) ,
and
v M ( u ( τ ) , v ( τ ) ) ( v v ( τ ) ) = n = 1 N k = 1 K 2 ( u n k ( τ ) + v n k ( τ ) ) ( v n k v n k ( τ ) ) .
For a given feasible point u ( τ ) , v ( τ ) , p ( τ ) , and q ( τ ) , we can achieve the upper bound of (35) by solving the following convex optimization problem:
minimize u , v , p , q F A ( u , p ) u ( τ ) , p ( τ ) affine G A ( u , p ) + η ( H ( u , v ) u ( τ ) , v ( τ ) affine M ( u , v ) ) s . t . C 1 a , C 1 b , C 2 , C 3 , C 6 b C 4 ( a ) : F n A ( u , p ) + G n B ( v , q ) u ( τ ) , p ( τ ) affine G n A ( u , p ) v ( τ ) , q ( τ ) affine F n B ( v , q ) 0 ( b ) : R ^ n , target + G n A ( u , p ) u ( τ ) , p ( τ ) affine F n A ( u , p ) 0 ,
where u ( τ ) , p ( τ ) affine G A ( u , p ) = n = 1 N u ( τ ) , p ( τ ) affine G n A ( u , p ) .
Generally, the convex problem in (48) can be readily settled by standard convex program solvers, and it can be solved by standard convex programming solvers such as CVX [33]. We propose a successive convex approximation to tighten the upper bound solution in (48) by an iterative algorithm, i.e., Algorithm 2. It can generate a sequence of feasible solutions continuously and achieve a locally optimal solution in polynomial time [34].
Algorithm 2 Iterative Algorithm for Resource Allocation.
1:
Initialization: Initialize the maximum number of iterations T max and set iteration index τ = 1 ; and initialize the current feasible point u ( 1 ) , v ( 1 ) , p ( 1 ) , q ( 1 ) and a penalty factor η 1
2:
repeat
3:   
Set the variables u , v , p , and q to be solved by the standard convex program solvers
4:   
Evaluate the convex functions F n A ( u , p ) , G n A ( u , p ) , G n B ( v , q ) , and H ( u , v )
5:   
According to the current point u ( τ ) , v ( τ ) , p ( τ ) , and q ( τ ) , evaluate the affine functions u ( τ ) , p ( τ ) affine F n A ( u , p ) , u ( τ ) , p ( τ ) affine G n A ( u , p ) , v ( τ ) , q ( τ ) affine F n B ( v , q ) , and u ( τ ) , v ( τ ) affine M ( u , v )
6:   
Substitute them into (48) to solve the convex problem for getting the upper-bound optimal point u opt , v opt , p opt , and q opt
7:   
Update iteration index τ = τ + 1
8:   
Update the next iteration point u ( τ ) = u opt , v ( τ ) = v opt , p ( τ ) = p opt , q ( τ ) = q opt
9:
until ( u ( τ 1 ) , v ( τ 1 ) , p ( τ 1 ) , q ( τ 1 ) ) = ( u ( τ ) , v ( τ ) , p ( τ ) , q ( τ ) ) or τ = T max
10:
Output the suboptimal solution u * = u ( τ ) , v * = v ( τ ) , p * = p ( τ ) , q * = q ( τ )

4.3. Computational Complexity Analysis

The computational complexity of an exhaustive search in the grouping matching algorithm is O ( 2 M N ) . The exhaustive search scheme is user-centric, as it divides each UE into a group, and each relay can either belong to the group of the UE or not. Thus, the solution to all groupings is 2 M N , and the computational complexity of the exhaustive search is O ( 2 M N ) . The computational complexity of Algorithm 1 is O ( N M 2 ) . Specifically, N M steps are needed, while each UE matches with a relay for grouping in the proposed grouping algorithm, and the steps for the grouping process are less than M · N M steps. Therefore, the total computational complexity of the proposed grouping algorithm is O ( N M 2 ) . The computational complexity of the D.C. programming is O ( T max M ) ; as T max is no more than Q M , the computational complexity of the D.C. programming is O ( M ) . Thus, the computational complexity is O ( N M 3 ) .

4.4. Convergence of Algorithm 1

We divide the algorithm into two processes, and the first process guarantees UE communications. The system performance is slightly degraded to satisfy QoS. The second process of the algorithm is convergent, and the proof is as follows.
Proof. 
When UE n a matches with relays m 2 , m 3 , …, and m Q (the descending order of priority) in the second process, assuming that there exists a relay m q matching UE n b , the priority of m q for n a is higher than that of m Q . n a is higher than n b in the priority list of relay m q simultaneously. □
When UE n a matches relay m Q not m q , there are two situations: UE n a sends a request to m q , and relay m q is rejected, which indicates that the priority of n b is higher than n a in the priority list of m q . n a does not send the request to m q . We can conclude that the priority of m Q is higher than m q in the priority list of n a . The two situations of the hypothesis cannot exist simultaneously, and thus, the hypothesis is not true. There are no better matched pairs, so the matched pair obtained is stable.

5. Simulation Results and Analysis

In this section, we evaluate the proposed framework and algorithms in terms of the system throughput through simulations. To make a fair comparison, we try to use the same system configuration in OMA, Co-OMA and the traditional NOMA system with the proposed scheme. We deploy the BS in the middle of a 1000 m × 1000 m area. The UEs and relays are modeled as independent PPP with density λ u and λ r . UEs are generated in the area randomly, as illustrated in the top of Figure 2. Other simulation parameters are summarized in Table 2. To show the grouping directly, we provide the schematic diagram when the maximum stack semaphore Q in our proposed scheme is set as 3. The matched pairs of relays and UEs chosen by our proposed algorithms are shown at the bottom of Figure 2. The codes are developed on MATLAB using the CVX toolbox and are executed on a 64-bit operating system with 16 GB RAM and Intel CORE i7, 3.4 GHz.
In Figure 3, we compared the system throughput of the proposed NOMA-based cooperative network scheme, the OMA-based cooperative network (Co-OMA) scheme, the NOMA, and the OMA scheme under different densities of UEs. The density of relays was set to 300. According to Figure 3, we can see that the proposed scheme achieved the highest system throughput, and it exhibited a 50% gain when the density of UEs exceeded 200 compared to the Co-OMA scheme. This is because with the increasing density of UEs, the spectrum resources is limited, and the proposed scheme shows the advantage of using non-orthogonal resources.
In addition, NOMA and OMA solutions without collaborative communication technology can provide high overall throughput when the density of user nodes becomes too high. This is because to ensure the QoS of some weak channel UEs, a large number of system resources are sacrificed to the specific UEs, which leads to the slow decline of the total throughput of the system and eventually tends to be stable. The comparison between NOMA and cooperative OMA in [35] shows the relationship between backhaul capacity and micro-area access number. When the number of UEs served by the system exceeds the threshold value, the system performance will decline OMA with sufficient system resources (low user node density). However, the user node density increases, and the cooperative OMA can improve the throughput by using the channel gain of the backhaul link and the multiplexing gain of a large number of relay nodes, thus exceeding the throughput performance of a non-collaborative NOMA scenario.
To verify the effectiveness of our proposed algorithms, we compared it with the following benchmarks in the NOMA-based cooperative network: (1) opportunistic with DCP-based allocation (ORG-DCPA), (2) dynamic relay grouping with fixed transmit power allocation (DRG-FTPA), and (3) opportunistic relay grouping with fixed transmit power allocation (ORG-FTPA). The four algorithms can be seen as the combination of DRG, ORG, DCPA and FTPA. The computational complexity of DRG is O ( N · M 2 ) as shown in the analyses above; the worst case of ORG is exhaustive search, so the computational complexity is O ( N · 2 M ) . The computational complexity of DCPA and FTPA are O ( M ) and O ( 1 ) , respectively. Thus, the complexity of DRG-DCPA, ORG-DCPA, DRG-FTPA, and ORG-FTPA can be present as O ( N · M 3 ) , O ( M · N · 2 M ) , O ( N · M 2 ) , and O ( N · 2 M ) . It needs to be emphasized that although DRG-FTPA has the lowest complexity, due to the unchangeable natuer of the transmit power, the DRG-FTPA shows the worst performance during the simulation. Then, we showed the simulation results in terms of the maximum power of relays, the density of relays and the density of UEs.
Figure 4 shows the system throughput diagram of several algorithm schemes under different relay transmitting powers, where the densities of UEs and relays are 200 and 500, respectively. The proposed DRG-DCPA resource allocation algorithm can obtain the highest system throughput, which is followed by the ORG-DCPA algorithm and ORG-FTPA algorithm. In particular, compared with the system throughput of other resource allocation algorithms, the DRG-DCPA algorithm can at least double the system throughput when the maximum power is more than 15 dBm.
Figure 5 shows the system throughput comparison of the above resource allocation algorithm schemes under different relay densities where the user node density is set to 50. When the density of relays is large enough, the system throughput of the DRG-DCPA resource allocation algorithm improved by more than 60%, 100% and 200%. When the density of intermediate relays increases, the system throughput of the DRG-DCPA algorithm improves faster, indicating that the algorithm has a stronger ability to utilize relay resources and thus can obtain more multiplexing gain.
Figure 6 shows the system throughput comparison of the above resource allocation algorithm schemes for different densities of UE. The density of the relay is set to 300. When the relay node density is sufficient, the throughput of the DRG-DCPA algorithm will increase with increasing user density, and when the user node density is 200, the throughput of the DRG-DCPA algorithm will increase by more than 60%, 120% and 210% compared with the other three algorithms. In addition, the use of a fixed percentage of the power allocation algorithm under the condition of excessive user node density decreases. This is because the above allocation algorithm may spend too much power to weak channel users to ensure the QoS. This causes the system to not effectively utilize resources; when the user node is too saturated, the system throughput will deduce.
In Figure 7, we present the effect in the proposed collaborative system on throughput with different maximum numbers of signals superposed per subcarrier. When the maximum relay power is 15 dBm, the proposed system can improve the throughput more than twice under QoS constraints. In addition, we use Q to denote the maximum number of signals superposed per subcarrier; when Q is larger, the proposed scheme can achieve a higher throughput. In particular, the throughput increases rapidly when the intermediate maximum power is higher. One possible reason is that the number of overstacked signals at low SINR will make the interference in the SIC decoder too large to meet the constraint C3; thus, the gain of the throughput is reduced.
It should be noted that the influence of the direct path is not considered in the reference, and the channel condition is poor. The simulation comparison result is the lower bound of the system and algorithm performance, but it can also effectively demonstrate the effectiveness of the proposed collaborative system and optimization algorithm. To confirm the above conjecture, we give a further simulation. In this simulation, we consider the direst path scenario. In addition, if the impact of the direct path on the system is considered, there is generally a direct path between the BS and the relay, but there is no direct path between the relay and the user. At this point, the channel between the base station and the relay is better, which is generally regarded as the Rician fading channel, and the improvement of the channel condition can help relax the constraint C3. The simulation shows that in the low SINR, Q = 3 with direct path performances better than Q = 4, but as the power increased, the Q = 4 without a direct solution plays a non-orthogonal higher resource utilization. The simulation result can further improve our conjecture.

6. Conclusions

In this paper, the NOMA-based multi-user and multi-relay cooperative network has been studied. To maximize the system throughput, we have designed the resource allocation algorithm as a mixed integer non-linear programming problem. To improve its tractability, we have divided the problem between (1) dynamic group matching of relays and UEs and (2) DCP-based joint allocation of power and subcarriers. Simulation results have confirmed that higher system throughput can be achieved through the proposed algorithm. Compared with Co-OMA, OMA, and NOMA, the proposed algorithm had the highest throughput. The proposed algorithm can also increase the system throughput substantially when the maximum power of the relays is high. The superiority of the proposed algorithm was substantiated by comparing it with different algorithms under various user density and relay density configurations. Simulation results confirmed that the proposed algorithms can be appropriately applied to IoT scenarios with massively small UEs.

Author Contributions

Conceptualization, J.Z.; methodology, J.Z.; software, J.S. and J.M.; validation, J.Z., J.S. and Y.S.; writing—original draft preparation, J.S. and Y.S.; writing—review and editing, J.Z., J.S. and Y.S.; supervision, J.Z., J.M., T.L. and S.Z.; project administration, J.Z., T.L. and S.Z.; funding acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 62001264), the Natural Science Foundation of Beijing (No. L192025), and the National Key R & D Program (No. 2018YFB1801102).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
NOMANon-orthogonal multiple access
IoTInternet of things
OMAOrthogonal multiple access
SCMASparse code multiple access
PDMAPattern division multiple access
MIMOMultiple input multiple output
DCPDifference of convex functions programming
QoSQuality of service
BSBase station
UEUser equipment
PPPsPoisson point processes
TDDTime division duplex
AWGNAdditive white Gaussian noise
SICSuccessive interference cancellation
SINRSignal-to-interference-plus-noise ratio
Co-OMAOMA-based cooperative network
ORG-DCPAOpportunistic with DCP-based allocation
DRG-FTPADynamic relay grouping with fixed transmit power allocation
ORG-FTPAOpportunistic relay grouping with fixed transmit power allocation

References

  1. Chen, S.; Liang, Y.C.; Sun, S.; Kang, S.; Cheng, W.; Peng, M. Vision, requirements, and technology trend of 6G: How to tackle the challenges of system coverage, capacity, user data-rate and movement speed. IEEE Wirel. Commun. 2020, 27, 218–228. [Google Scholar] [CrossRef]
  2. Di, B.; Zhang, H.; Song, L.; Li, Y.; Han, Z.; Poor, H.V. Hybrid beamforming for reconfigurable intelligent surface based multi-user communications: Achievable rates with limited discrete phase shifts. IEEE J. Sel. Areas Commun. 2020, 38, 1809–1822. [Google Scholar] [CrossRef]
  3. Shahab, M.B.; Abbas, R.; Shirvanimoghaddam, M.; Johnson, S.J. Grant-free non-orthogonal multiple access for IoT: A survey. IEEE Commun. Surv. Tutor. 2020, 22, 1805–1838. [Google Scholar] [CrossRef]
  4. Liaqat, M.; Noordin, K.A.; Abdul Latef, T.; Dimyati, K. Power-domain non orthogonal multiple access (PD-NOMA) in cooperative networks: An overview. Wirel. Netw. 2020, 26, 181–203. [Google Scholar] [CrossRef]
  5. Ye, N.; Li, X.; Yu, H.; Zhao, L.; Liu, W.; Hou, X. DeepNOMA: A Unified Framework for NOMA Using Deep Multi-Task Learning. IEEE Trans. Wirel. Commun. 2020, 19, 2208–2225. [Google Scholar] [CrossRef]
  6. Dai, L.; Wang, B.; Ding, Z.; Wang, Z.; Chen, S.; Hanzo, L. A survey of non-orthogonal multiple access for 5G. IEEE Commun. Surv. Tutor. 2018, 20, 2294–2323. [Google Scholar] [CrossRef]
  7. Zuo, J.; Liu, Y.; Al-Dhahir, N. Reconfigurable intelligent surface assisted cooperative non-orthogonal multiple access systems. IEEE Trans. Commun. 2021, 69, 6750–6764. [Google Scholar] [CrossRef]
  8. Ding, Z.; Lei, X.; Karagiannidis, G.K.; Schober, R.; Yuan, J.; Bhargava, V.K. A survey on non-orthogonal multiple access for 5G networks: Research challenges and future trends. IEEE J. Sel. Areas Commun. 2017, 35, 2181–2195. [Google Scholar] [CrossRef]
  9. Zhang, X.; Zhang, D.; Yang, L.; Han, G.; Chen, H.H.; Zhang, D. Scma codebook design based on uniquely decomposable constellation groups. IEEE Trans. Wirel. Commun. 2021, 20, 4828–4842. [Google Scholar] [CrossRef]
  10. da Silva, B.F.; Silva, D.; Uchôa-Filho, B.F.; Le Ruyet, D. A multistage method for SCMA codebook design based on MDS codes. IEEE Wirel. Commun. Lett. 2019, 8, 1524–1527. [Google Scholar] [CrossRef]
  11. Han, S.; Huang, Y.; Meng, W.; Li, C.; Xu, N.; Chen, D. Optimal power allocation for SCMA downlink systems based on maximum capacity. IEEE Trans. Commun. 2018, 67, 1480–1489. [Google Scholar] [CrossRef]
  12. Xiang, L.; Liu, Y.; Xu, C.; Maunder, R.G.; Yang, L.L.; Hanzo, L. Iterative receiver design for polar-coded SCMA systems. IEEE Trans. Commun. 2021, 69, 4235–4246. [Google Scholar] [CrossRef]
  13. Dai, X.; Zhang, Z.; Bai, B.; Chen, S.; Sun, S. Pattern division multiple access: A new multiple access technology for 5G. IEEE Wirel. Commun. 2018, 25, 54–60. [Google Scholar] [CrossRef]
  14. Tang, W.; Kang, S.; Ren, B.; Yue, X. Uplink grant-free pattern division multiple access (GF-PDMA) for 5G radio access. China Commun. 2018, 15, 153–163. [Google Scholar] [CrossRef]
  15. Zhang, P.; Yang, X.; Chen, J.; Huang, Y. A survey of testing for 5G: Solutions, opportunities, and challenges. China Commun. 2019, 16, 69–85. [Google Scholar]
  16. Shah, A.S.; Islam, M.S. A survey on cooperative communication in wireless networks. Int. J. Intell. Syst. Appl. 2014, 6, 66. [Google Scholar] [CrossRef]
  17. Lv, L.; Chen, J.; Ni, Q.; Ding, Z.; Jiang, H. Cognitive non-orthogonal multiple access with cooperative relaying: A new wireless frontier for 5G spectrum sharing. IEEE Commun. Mag. 2018, 56, 188–195. [Google Scholar] [CrossRef]
  18. Nosratinia, A.; Hunter, T.E.; Hedayat, A. Cooperative communication in wireless networks. IEEE Commun. Mag. 2004, 42, 74–80. [Google Scholar] [CrossRef]
  19. Wang, M.; Gao, F.; Jin, S.; Lin, H. An overview of enhanced massive MIMO with array signal processing techniques. IEEE J. Sel. Top. Signal Process. 2019, 13, 886–901. [Google Scholar] [CrossRef]
  20. Pan, J.; Ye, N.; Yu, H.; Hong, T.; Alrubaye, S.; Mumtaz, S.; Chih-Lin, I. AI-Driven Blind Signature Classification for IoT Connectivity: A Deep Learning Approach. IEEE Trans. Wirel. Commun. 2022, 21, 6033–6047. [Google Scholar] [CrossRef]
  21. Ye, N.; Yu, J.; Wang, A.; Chen, S.; Liang, Y.C.; Sun, S.; Kang, S.; Cheng, W.; Peng, M. Help from space: Grant-free massive access for satellite-based IoT in the 6G era. Digit. Commun. Netw. 2022, 8, 215–224. [Google Scholar] [CrossRef]
  22. Ding, Z.; Peng, M.; Poor, H.V. Cooperative non-orthogonal multiple access in 5G systems. IEEE Commun. Lett. 2015, 19, 1462–1465. [Google Scholar] [CrossRef]
  23. Zeng, M.; Hao, W.; Dobre, O.A.; Ding, Z. Cooperative NOMA: State of the art, key techniques, and open challenges. IEEE Netw. 2020, 34, 205–211. [Google Scholar] [CrossRef]
  24. Luo, L.; Li, L.; Su, X. Optimization of resource allocation in relay assisted multi-user SCMA uplink network. In Proceedings of the 2017 International Conference on Computing, Networking and Communications (ICNC), Santa Clara, CA, USA, 26–29 January 2017; pp. 282–286. [Google Scholar]
  25. Han, S.; Zhang, J.; Guo, C.; Liu, N. Full-duplex MIMO relay system design based on SCMA. In Proceedings of the 2017 IEEE International Conference on Communications (ICC), Paris, France, 21–25 July 2017; pp. 1–5. [Google Scholar]
  26. Tang, W.; Kang, S.; Ren, B. Performance analysis of cooperative pattern division multiple access (Co-PDMA) in uplink network. IEEE Access 2017, 5, 3860–3868. [Google Scholar] [CrossRef]
  27. Tang, W.; Kang, S.; Fu, X.; Yue, X.; Zhang, X. On the performance of PDMA with decode-and-forward relaying in downlink network. IEEE Access 2018, 6, 20113–20124. [Google Scholar] [CrossRef]
  28. Sun, Y.; Ng, D.W.K.; Ding, Z.; Schober, R. Optimal joint power and subcarrier allocation for full-duplex multicarrier non-orthogonal multiple access systems. IEEE Trans. Commun. 2017, 65, 1077–1091. [Google Scholar] [CrossRef]
  29. Fu, Y.; Salaün, L.; Sung, C.W.; Chen, C.S. Subcarrier and power allocation for the downlink of multicarrier NOMA systems. IEEE Trans. Veh. Technol. 2018, 67, 11833–11847. [Google Scholar] [CrossRef]
  30. Li, X.; Wang, Q.; Liu, M.; Li, J.; Peng, H.; Piran, M.J.; Li, L. Cooperative wireless-powered NOMA relaying for B5G IoT networks with hardware impairments and channel estimation errors. IEEE Internet Things J. 2020, 8, 5453–5467. [Google Scholar] [CrossRef]
  31. Burer, S.; Letchford, A.N. Non-convex mixed-integer nonlinear programming: A survey. Surv. Oper. Res. Manag. Sci. 2012, 17, 97–106. [Google Scholar] [CrossRef]
  32. Liu, Y.; Li, X.; Yu, F.R.; Ji, H.; Zhang, H.; Leung, V.C. Grouping and cooperating among access points in user-centric ultra-dense networks with non-orthogonal multiple access. IEEE J. Sel. Areas Commun. 2017, 35, 2295–2311. [Google Scholar] [CrossRef]
  33. Grant, M.; Boyd, S. CVX: Matlab Software for Disciplined Convex Programming, Version 2.1. 2014. Available online: http://cvxr.com/cvx (accessed on 20 June 2020).
  34. Razaviyayn, M. Successive Convex Approximation: Analysis and Applications. Ph.D. Thesis, University of Minnesota, Minneapolis, MN, USA, 2014. [Google Scholar]
  35. Ju, J.; Duan, W.; Sun, Q.; Gao, S.; Zhang, G. Performance analysis for cooperative NOMA with opportunistic relay selection. IEEE Access 2019, 7, 131488–131500. [Google Scholar] [CrossRef]
Figure 1. An illustration of the two-stage transmission based on NOMA.
Figure 1. An illustration of the two-stage transmission based on NOMA.
Sensors 22 06023 g001
Figure 2. The deployment of BS, UEs, and relays in a 1000 m × 1000 m area. The cell boundaries form a Voronoi tessellation.
Figure 2. The deployment of BS, UEs, and relays in a 1000 m × 1000 m area. The cell boundaries form a Voronoi tessellation.
Sensors 22 06023 g002
Figure 3. The system throughput versus the density of UEs for different schemes.
Figure 3. The system throughput versus the density of UEs for different schemes.
Sensors 22 06023 g003
Figure 4. System throughput versus the maximum power of the relays with different resource allocation algorithms.
Figure 4. System throughput versus the maximum power of the relays with different resource allocation algorithms.
Sensors 22 06023 g004
Figure 5. The system throughput versus the density of relays for using different optimization algorithms.
Figure 5. The system throughput versus the density of relays for using different optimization algorithms.
Sensors 22 06023 g005
Figure 6. The system throughput versus the density of UEs for using different optimization algorithms.
Figure 6. The system throughput versus the density of UEs for using different optimization algorithms.
Sensors 22 06023 g006
Figure 7. The system throughput versus the maximum power of the relays with different maximum numbers of superposed signals.
Figure 7. The system throughput versus the maximum power of the relays with different maximum numbers of superposed signals.
Sensors 22 06023 g007
Table 1. Notations of symbols.
Table 1. Notations of symbols.
SymbolsDescription
λ u The density of UE
λ r The density of relay
α Path loss exponent
β k The proportion of a time slot occupied by the access stage on the kth subcarrier
x m n k Signal transmitted from the mth relay to the nth UE on the kth subcarrier
p m n k Transit power of the signal x m n k
h m n k Small-scale fading channel coefficients of x m n k
ϖ m n Large-scale fading channel coefficients of x m n k
y n A , k Signal received on the kth subcarrier of the nth UE
x m k Signal transmitted from the BS to the mth relay on the kth subcarrier
q m k The power of x m k
h m k Small-scale fading channel coefficients of x m k
ϖ m Large-scale fading channel coefficients of x m k
y m B , k Signal received on the kth subcarrier of the mth relay
S I N R m n A , k The received SINR of the nth UE served by the mth relay on the kth subcarrier
R ˜ n A The forward throughput of the nth UE
S I N R m B , k The signal-to-interference-plus-noise ratio (SINR) of the mth relay on the kth subcarrier
R ˜ n B The backhaul throughput for the nth UE
R ˜ The whole system throughput
Table 2. Parameters setting.
Table 2. Parameters setting.
ParametersValue
The system bandwidth3 MHz
Carrier frequency2 GHz
The subcarrier bandwidth15 KHz
Path loss exponent α 3.5
Noise power spectral density−174 dBm/Hz
The density of relays λ r 500 per/km2
The density of UEs λ u 200 per/km2
The maximum power of BS46 dBm
The maximum power of relays0:3:30 dBm
The target data rate of UEs2 bps/Hz
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zeng, J.; Sun, J.; Song, Y.; Mei, J.; Lv, T.; Zhou, S. Resource Allocation in Multi-Carrier Multiplexed NOMA Cooperative System. Sensors 2022, 22, 6023. https://doi.org/10.3390/s22166023

AMA Style

Zeng J, Sun J, Song Y, Mei J, Lv T, Zhou S. Resource Allocation in Multi-Carrier Multiplexed NOMA Cooperative System. Sensors. 2022; 22(16):6023. https://doi.org/10.3390/s22166023

Chicago/Turabian Style

Zeng, Jie, Jiaying Sun, Yuxin Song, Jiajia Mei, Tiejun Lv, and Shidong Zhou. 2022. "Resource Allocation in Multi-Carrier Multiplexed NOMA Cooperative System" Sensors 22, no. 16: 6023. https://doi.org/10.3390/s22166023

APA Style

Zeng, J., Sun, J., Song, Y., Mei, J., Lv, T., & Zhou, S. (2022). Resource Allocation in Multi-Carrier Multiplexed NOMA Cooperative System. Sensors, 22(16), 6023. https://doi.org/10.3390/s22166023

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop