Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Cost-Driven Method for Deep-Learning-Based Hardware Trojan Detection
Previous Article in Journal
ECG-Based Stress Detection and Productivity Factors Monitoring: The Real-Time Production Factory System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Design of a Simultaneous Reflection and Transmission RIS in Mode-Switching Mode to Assist NOMA Systems

The College of Information, Mechanical and Electrical Engineering, Shanghai Normal University, Shanghai 200234, China
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(12), 5504; https://doi.org/10.3390/s23125504
Submission received: 30 April 2023 / Revised: 23 May 2023 / Accepted: 8 June 2023 / Published: 11 June 2023
(This article belongs to the Section Communications)

Abstract

:
Simultaneous transmitting and reflecting reconfigurable intelligent surfaces (STAR-RISs) can reflect signals and transmissive signals simultaneously and can extend the coverage of signals. A conventional RIS mainly focuses on the case where the signal source and the target are on the same side. In this paper, a STAR-RIS-assisted non-orthogonal multiple access (NOMA) downlink communication system is considered to maximize the achievable rate for users by jointly optimizing the power-allocation coefficients, active beamforming, and STAR-RIS beamforming under the mode-switching (MS) protocol. The critical information of the channel is first extracted using the Uniform Manifold Approximation and Projection (UMAP) method. Based on the key extracted channel features, STAR-RIS elements and users are clustered individually using the fuzzy C-mean clustering (FCM) method. The alternating optimization method decomposes the original optimization problem into three sub-optimization problems. Finally, the sub-problems are converted to unconstrained optimization methods using penalty functions for the solution. Simulation results show that when the number of elements of RIS is 60, the achievable rate of the STAR-RIS–NOMA system is about 18% higher than that of the RIS–NOMA system.

1. Introduction

With the rapid development of modern mobile communication technologies, the capacity of communication networks will grow significantly in the next decade, and the high complexity of networks, the high cost of hardware, and the increasing consumption of resources have become the key issues facing wireless communication in the future [1]. Reconfigurable intelligent surfaces (RISs) are low-cost, programmable, easy to deploy, and low-energy, which makes RISs stand out in competition with other technologies. RISs and their variants [2,3] have become promising technologies for sixth-generation (6G) wireless networks. Simultaneous transmitting and reflecting (STAR)-RIS is a new technique in RISs; whereas a conventional RIS focuses on the case where the signal source and the target are located on the same side of the RIS [4,5,6], a STAR-RIS is capable of simultaneous transmission and reflection, and the target can be distributed on both sides of the signal source, allowing for full spatial coverage and enabling better-aided system communication [7,8,9].
Due to the significant capacity growth of future communication networks, non-orthogonal multiple access (NOMA) [10,11,12,13,14] has also received considerable attention as a promising technology to enhance spectrum efficiency and support large-scale connectivity. NOMA transmits the information of multiple users on the same sub-channel. Interference information is introduced, in contrast to orthogonal multiple access (OMA) systems, where each user has a separate resource in the time, frequency, or code domain [15]. NOMA works through sharing time-frequency resources by separating users into other domains [16]. Power-domain-based NOMA maximizes performance by assigning different power levels to users. The receiver-side detection process of power-domain NOMA highly depends on successive interference cancellation (SIC) [17]. The SIC receiver performs interference cancellation in a specific order according to different users’ signal power levels and distinguishes between other users [18]. The main advantages of NOMA are: (1) NOMA’s resource allocation method is based on a non-orthogonal, and this approach can improve spectrum efficiency. (2) NOMA uses a scheduling-free access scheme, which can better meet latency and reliability requirements. (3) NOMA can be more easily combined with other coding schemes. Combining RIS with NOMA can provide services to more users and result in higher network performance [19,20].
In the study of assisted communication, the design of the beamforming is a major concern for researchers. Ref. [21] considered joint beamforming for RIS-assisted hybrid satellite–terrestrial relay networks to minimize the total transmit power of satellites and base stations (BSs) while ensuring user rate requirements. In [22], a confidential and energy-efficient hybrid beamforming scheme for a satellite-ground integrated network was studied, and a hybrid beamformer for the base station and a digital beamformer for the satellite was jointly designed. The problem of maximizing spectral efficiency was considered in [23] for RIS-assisted single-user communication systems. In [24], the authors applied active RISs to collaboratively enhance secure transmission from the base station to the mobile user to maximize the achievable secrecy rate subject to transmit power constraints and interference thresholds. The authors in [25] used RISs to augment artificial noise to assist multiple-input multiple-output (MIMO) secure communication systems. In [26], a Compressed Sensing (CS)-based channel estimation method was proposed by exploiting the sparse representation of cascaded channels. In [27], a globally optimal solution for a single-user system’s transmit-power-minimization problem was obtained using a branch-and-bound approach. The authors of [28,29] used semidefinite relaxation and streamlined optimization, respectively, to maximize the received signal-to-noise ratio of the multiple-input single-output (MISO) communication systems for RISs. Ref. [30] solved the RIS-assisted NOMA downlink-reachability and rate-maximization problem using semi-positive definite planning (SDP). Ref. [31] derived uplink-outage probability expressions for blind RIS–NOMA using RIS as an innovative reflector and RIS as an access point (AP). The results showed that adding RIS elements enhanced the total capacity and interruption performance. Some researchers have also conducted in-depth studies on STAR-RIS-assisted communication. In [32], the authors studied STAR-RIS-assisted dual-user communication systems and proposed the problem of maximizing the total coverage of NOMA and OMA systems. In [33], the power-minimization problem for STAR-RIS-assisted unicast and multicast systems was studied for different operating protocols of STAR-RIS, and the active-passive beamforming was jointly optimized. In [34], the authors proposed the index modulation (IM)-assisted STAR-RIS–NOMA system to improve spectral efficiency. In [35], the authors considered the STAR-RIS-assisted full-duplex (FD) wireless communication system. They derived analytical closed-form expressions for the outage probability and achievable throughput for uplink and downlink communication. Ref. [36] studied the STAR-RIS-assisted NOMA downlink communication system, maximizing the signal-to-noise ratio (SNR) to design the reflected and transmitted phase-shift matrices. In [37], the power of the STAR-RIS-enabled uplink NOMA system was minimized by a joint optimization scheme. The authors of [38] investigated the energy-saving potential of STAR-RIS in a MIMO-enabled NOMA system. However, the studies mentioned above on STAR-RIS are all processed directly on the high-dimensional channel, which is challenging to solve.
Based on the above background, integrating STAR-RIS with the NOMA system is of great practical importance. The main contributions of this paper are summarized as follows.
In the STAR-RIS–NOMA system, we jointly optimize the user’s power-allocation coefficient, active beamforming at the BS, and passive beamforming at the STAR-RIS to maximize the achievable sum rate. Since the number of antennas at the BS and the number of STAR-RIS elements are large, and thus the associated channel dimension is high, we use Uniform Streaming Approximation and Projection (UMAP) to extract the key features of the channel and reduce the complexity of the solution. The elements of the STAR-RIS are clustered to improve the efficiency of the whole communication system.
We decouple the original problem into three sub-problems for solving a given decoding order. Regarding the optimization problem of the power-allocation coefficients, we use SCA with a penalty function to solve it. After adopting the SDR relaxation constraint, the block coordinate descent (BCD) method is used to solve the active beamforming at the BS and the passive beamforming at the STAR-RIS.
Simulation results show: (1) combining the STAR-RIS with NOMA can significantly improve the system’s performance, and the system performance increases with the number of STAR-RIS elements; (2) the performance of the STAR-RIS–NOMA system is better than that of the reflection/transmission RIS–NOMA-only system.
The rest of the paper is organized as follows: Section 2 presents the STAR-RIS–NOMA system’s communication model and problem representation. Section 3 performs feature extraction of the channel and clusters the users and STAR-RIS elements. A new algorithm for alternating optimization is proposed in Section 4 to solve the joint optimization problem. Numerical results are given in Section 5. Section 6 concludes the whole paper.
Notations: Scalars, vectors and matrices are represented in lowercase, bold lowercase and uppercase letters, respectively. N r i s × 1 is a complex vector of dimension N r i s . d i a g u represents a diagonal matrix whose diagonal elements are the corresponding elements of the vector u . u H denotes the conjugate transpose of the vector u . T r U and R a n k U denote the trace and rank of matrix U , respectively. C N 0 , σ 2 denotes a complex Gaussian distribution with mean zero and variance σ 2 . U 0 denotes U as a semi-positive definite matrix.

2. System Model and Problem Formulation

2.1. System Mode

As shown in Figure 1, consider the downlink transmission of the STAR-RIS–NOMA system, where the base station (BS) communicates with M single-antenna users, the set of users is M 1 , , M , and the direct link between the users and the BS is blocked. Assume that the BS is equipped with N t transmit antennas, while the STAR-RIS is equipped with N r i s elements and the set of elements is N r i s 1 , , N r i s . All perfect CSIs are available at the BS. H N r i s × N t is the channel from the BS to the STAR-RIS, and g m N r i s × 1 is the channel from the STAR-RIS to the m - th user. We use fuzzy C-mean (FCM) clustering to group M users into L clusters. The set of users in the l - th cluster is denoted by L l . Thus, the number of users in the l - th cluster can be denoted by L l , where l = 1 L L l = M . The detailed process on user clustering will be discussed in Section 3.
This paper considers the user-achievable sum-rate maximization problems under the mode-switching (MS) protocol. Under the MS protocol, STAR-RIS elements are divided into two groups: one working in reflection mode (R mode) and the other in transmission mode (T mode). Let u t = β 1 t e j ϕ 1 t , , β N r i s t e j ϕ N r i s t H and u r = β 1 r e j ϕ 1 r , , β N r i s r e j ϕ N r i s r H , where β n r i s t , β n r i s r 0 , 1 , β n r i s t + β n r i s r = 1 and ϕ n r i s t , ϕ n r i s r [ 0 , 2 π ) ,   n r i s N r i s . The corresponding STAR-RIS transmission and reflection coefficient matrices are Φ t = d i a g u t and Φ r = d i a g u r ., respectively.
In the STAR-RIS–NOMA system, v l denotes the active beamforming vector of cluster l , and s l , c denotes the signal sent by the BS to the c - th user of the l - th cluster. Therefore, the received signal of the c - th user of the l - th user cluster is given by
y l , c = g l , c H Φ l , c H v l η l , c s l , c + i c , i L l g l , c H Φ l , c H v l η l , i s l , i + j l , j L i L j g l , c H Φ l , c H v j η j , i s j , i + n l , c ,
where η l , c denotes the power-allocation coefficient of the c - th user of the l - th user -cluster, n l , c represents the additive Gaussian white noise, and n l , c C N 0 , σ 2 , σ 2 is the noise power. Due to the presence of multiple users, each user uses the SIC sequentially to eliminate interference. This paper considers the user-achievable sum-rate maximization problem in the case of fixed decoding order. For simplicity and without loss of generality, we assume that the equivalent channel gain order in the l - th user cluster is
g l , 1 H Φ l , 1 H v l 2 g l , 2 H Φ l , 2 H v l 2 g l , L l H Φ l , L l H v l 2 .
Thus, the achievable rate for the c - th user of the l - th user cluster can be given as
R l , c = log 2 1 + g l , c H Φ l , c H v l 2 η l , c g l , c H Φ l , c H v l 2 i < c η l , i + j l , j L i L j g l , c H Φ l , c H v j 2 η l , i + σ 2 .

2.2. Problem Formulation

To improve the overall data rate, the joint power-allocation coefficients, active beamforming, and transmission and reflection beamforming optimization problems at the RIS are proposed to maximize the achievable rate for users. The problem formulations are
max η l , c , v l , u o l L c L l R l , c
s . t .   R l , c R l , c min ,
l L v l 2 P max ,
ϕ n r i s o [ 0 , 2 π ) ,   n r i s N r i s ,
β n r i s o 0 , 1 , β n r i s t + β n r i s r = 1 ,
where o { t , r } . Constraint (4b) guarantees a minimum quality of service (QoS) requirement for each user, and constraint (4c) indicates a total transmit power budget of P max . Since the objective function is nonconvex, the constraint (4b)–(4d) is also nonconvex. Therefore, problem (4) is a nonconvex optimization problem, which is more difficult to solve directly. Compared with reflection-only RIS, the STAR-RIS has a more complex optimization problem because of the presence of both reflection and transmission modes.

3. Channel Feature Extraction and Clustering of Users and STAR-RIS Elements

The number of transmitting antennas arranged at the BS is usually high. The number of STAR-RIS elements is significant, resulting in a high dimensionality of the relevant channels, which is challenging to optimize directly. The Uniform Manifold Approximation and Projection (UMAP) [39] can identify the critical structures in the high-dimensional channel space and preserve these structures in the low-dimensional embedding. Therefore, we first process the channel from the BS to the STAR-RIS and the STAR-RIS to the user using UMAP to extract the critical information and then perform clustering for STAR-RIS elements and system users, respectively.

3.1. Channel Feature Extraction Using UMAP

According to the three assumptions of UMAP about the data, a manifold with a fuzzy topology can be modeled. The embedding is found by searching the low-dimensional projection of the data with the closest equivalent fuzzy topology. The process of channel feature extraction using UMAP is shown in Figure 2.
The STAR-RIS-user channel matrix is G = [ g 1 , , g M ] . To extract the features of the STAR-RIS-user channel, the first step is to establish the relationship between the high-dimensional samples. In the STAR-RIS-user channel matrix G , the probability-density function between each point can be given as
p i | j ( G ) = e d ( g i g j ) ρ i ( G ) θ i ( G ) ,
where d denotes the Euclidean distance between sample points and ρ i ( G ) is the distance of the nearest point to point i . To ensure the symmetry of the distance, the joint probability is established as
p i j ( G ) = p i | j ( G ) + p j | i ( G ) p i | j ( G ) p j | i ( G ) .
Fix the R ( G ) nearest neighbors of the point i . Solve for θ i ( G ) according to R ( G ) = 2 i p i j ( G ) . The relationship between the low-dimensional points can be established by Equation (7).
q i j ( G ) = ( 1 + a ( G ) ( g ˜ i g ˜ j ) b ( G ) ) 1 ,
where a ( G ) , b ( G ) is the hyperparameter and g ˜ i is the downscaled channel. In order to describe the relationship between the high-dimensional sample points and the reduced-dimensional sample points, the cross-entropy-loss function is constructed as
C E ( g i j , g ˜ i j ) = i j p i j ( G ) ( g i j ) log ( p i j ( G ) ( g i j ) q i j ( G ) ( g ˜ i j ) ) + ( 1 p i j ( G ) ( g i j ) ) log ( 1 p i j ( G ) ( g i j ) 1 q i j ( G ) ( g ˜ i j ) ) ,
where g i j denotes the distance between the high-dimensional sample points and g ˜ i j denotes the distance between the low-dimensional sample points. Therefore, the low-dimensional space g ˜ i , corresponding to g i in the high-dimensional space, can be obtained by solving the following problem:
min   C E ( g i j , g ˜ i j )
We can use (10) to solve problem (9) using a stochastic gradient descent (SGD) update to obtain the reduced channel matrix G ˜ = g ˜ 1 , , g ˜ M .
g ˜ i j t + 1 = g ˜ i j t α 1 g ˜ q i j ( G ) ( g ˜ i j ) ,
where α 1 is the step size.
The channel from the BS to the STAR-RIS is H = h 1 , , h N r i s H . First, it is necessary to establish the relationship between the high-dimensional samples. The probability-density function between each point in the BS to the STAR-RIS channel matrix H can be given as
p i | j ( H ) = e d ( h i h j ) ρ i ( H ) θ i ( H ) ,
where d denotes the Euclidean distance between sample points, and ρ i ( H ) is fixed as the distance of the nearest point to point i . To ensure the symmetry of the distance, the joint probability is established:
p i j ( H ) = p i | j ( H ) + p j | i ( H ) p i | j ( H ) p j | i ( H ) .
Fix the R ( H ) nearest neighbors of the point i . Solve for θ i ( H ) according to R ( H ) = 2 i p i j ( H ) . The relationship between the low-dimensional points can be established by (13).
q i j ( H ) = ( 1 + a ( H ) ( h ˜ i h ˜ j ) b ( H ) ) 1 ,
where a ( H ) , b ( H ) is the hyperparameter and h ˜ i is the downscaled channel. To describe the relationship between the high-dimensional sample points and the downscaled low-dimensional sample points, a cross-entropy-loss function is constructed:
C E ( h i j , h ˜ i j ) = i j p i j ( H ) ( h i j ) log ( p i j ( H ) ( h i j ) q i j ( H ) ( h ˜ i j ) ) + ( 1 p i j ( H ) ( h i j ) ) log ( 1 p i j ( H ) ( h i j ) 1 q i j ( H ) ( h ˜ i j ) ) ,
where h i j denotes the distance between the high-dimensional sample points and h ˜ i j denotes the distance between the low-dimensional sample points. Therefore, h i in the high-dimensional space corresponds to the low-dimensional space h ˜ i , which can be obtained by solving the following problem:
min   C E ( h i j , h ˜ i j )
Problem (15) is solved by SGD using (16) to obtain the channel matrix H ˜ = h ˜ 1 , , h ˜ N r i s Τ of the reduced BS to the STAR-RIS.
h ˜ i j t + 1 = h ˜ i j t α 2 h ˜ q i j ( H ) ( h ˜ i j ) ,
where α 2 is the step size.

3.2. User Clustering

After extracting the STAR-RIS-user channel features using UMAP and then using FCM-clustering to subdivide the user clusters, it is possible to classify them more accurately with less complexity. The process of clustering users by FCM is shown in Figure 3. Using FCM, we divide the users into L clusters and the sub-cluster set into L = 1 , , L . The users are thus classified by the optimization problem (17).
min   J F C M ( u s e r ) Z , K = j = 1 L i = 1 M z j i χ ( u s e r ) g ˜ i k j 2
s . t .   j = 1 L z j i = 1 ,
where χ ( u s e r ) is an affiliation factor showing the importance of whether a sample belongs to a class or not, Z is the affiliation matrix z j i χ ( u s e r ) indicating the affiliation of sample i to class j , and K = k 1 , , k L is the cluster center. The constraint (17b) is combined into the objective function using the Lagrange multiplier method:
J F C M ( u s e r ) Z , K , ϑ m ( u s e r ) = j = 1 L i = 1 M z j i χ ( u s e r ) g ˜ i k j 2 + m = 1 M ϑ m ( u s e r ) j = 1 L z j m 1 ,
where ϑ m ( u s e r ) is the Lagrange multiplier. Taking the derivative of (18) with respect to z j i and k j , respectively, yields (19) and (20).
J z j i = χ ( u s e r ) g ˜ i k j 2 z j i χ ( u s e r ) 1 + ϑ i ( u s e r ) = 0 .
J k j = i = 1 M 2 z j i χ ( u s e r ) g ˜ i k j = 0 .
Solving Equations (19) and (20), we can obtain
z j i = r = 1 L g ˜ i k j 2 χ ( u s e r ) 1 g ˜ i k r 2 χ ( u s e r ) 1 1 .
k j = i = 1 M z j i χ ( u s e r ) g ˜ i i = 1 M z j i χ ( u s e r ) .
Equations (21) and (22) are iterated and, finally, the users are divided into L clusters. The set of users in the l - th user cluster is L l , and the number of users is L l .

3.3. STAR-RIS Element Clustering

Let X = G ˜ H H ˜ H = x 1 , , x N r i s ; then, also using FCM, divide the elements of the STAR-RIS into γ clusters and the set of sub-clusters into γ = 1 , , γ . The process of FCM-clustering of STAR-RIS elements is shown in Figure 4, and the classification results of STAR-RIS elements are obtained by the optimization problem (23):
min   J F C M ( r i s ) A , D = j = 1 γ i = 1 N r i s a j i χ ( r i s ) x ˜ i d j 2
s . t .   j = 1 γ a j i = 1 ,
where χ ( r i s ) is an affiliation factor showing the importance of whether a sample belongs to a class or not, A is the affiliation matrix, a j i χ ( r i s ) indicates the affiliation of sample i belonging to class j , and D = d 1 , , d L r i s is the cluster center. The constraint (23b) is combined into the objective function using the Lagrange multiplier method:
J F C M ( r i s ) A , D , ϑ m ( r i s ) = j = 1 γ i = 1 N r i s a j i χ ( r i s ) x ˜ i d j 2 + n = 1 N r i s ϑ n ( r i s ) j = 1 γ a j n 1 ,
where ϑ n ( r i s ) is the Lagrange multiplier. Equation (24) can be obtained by deriving Equations (25) and (26) for a j i and d j respectively.
J a j i = χ ( r i s ) x ˜ i d j 2 z j i χ ( r i s ) 1 + ϑ i ( r i s ) = 0 .
J d j = i = 1 N r i s 2 a j i χ ( r i s ) x ˜ i d j = 0 .
Solving Equations (25) and (26), we obtain
a j i = r = 1 γ x ˜ i d j 2 χ ( r i s ) 1 x ˜ i d r 2 χ ( r i s ) 1 1 .
d j = i = 1 N r i s a j i χ ( r i s ) x ˜ i i = 1 N r i s a j i χ ( r i s ) .
Equations (27) and (28) are iterated and, finally, the N r i s elements of the STAR-RIS are divided into γ clusters. The set of elements in the r - th RIS component cluster is γ r , and the number of elements is γ r .
After the channel feature extraction and the clustering operation, the received signal of the c - th user of the l - th user cluster transmitted/reflected through the r - th STAR-RIS element cluster is
y r , l , c = g ˜ r , l , c H Φ r , l , c H ˜ r v r , l η r , l , c s r , l , c + i c , i L l g ˜ r , l , c H Φ r , l , c H ˜ r v r , l η r , l , i s r , l , i + j l , j L i L j g ˜ r , l , c H Φ r , l , c H ˜ r v r , j η r , j , i s r , j , i + n r , l , c .
We assume that the power-allocation coefficients satisfy c L l η r , l , c = 1 . Thus, the achievable rate of the c - th user of the l - th user cluster transmitted/reflected through the r - th STAR-RIS cluster elements can be given as Equation (30).
R r , l , c = log 2 1 + g ˜ r , l , c H Φ r , l , c H ˜ r v r , l 2 η r , l , c g ˜ r , l , c H Φ r , l , c H ˜ r v r , l 2 i < c η r , l , i + j l , l L g ˜ r , l , c H Φ r , l , c H ˜ r v r , j 2 + σ 2 .
After clustering the STAR-RIS elements and users separately, the STAR-RIS element clusters correspond to the user clusters, in which case β n r i s t and β n r i s r are determined. Therefore, to maximize the achievable sum rate of M users, problem (4) is reformulated as
max η r , l , c , v r , l , u o l L c L l R r , l , c
s . t .   R r , l , c R r , l , c min ,
l L v r , l 2 P max ,
ϕ n r i s t , ϕ n r i s r [ 0 , 2 π ) ,   n r i s N r i s ,
c L l η r , l , c = 1 ,
where o { t , r } . Let U r . l . c = u r , l , c u r , l , c H , b r , l , c = d i a g g ˜ r , l , c H H ˜ r , B r . l . c = b r , l , c b r , l , c H and V r , l = v r , l v r , l H , then, the achievable rate is
R r , l , c = log 2 1 + η r , l , c T r U r , l , c B r . l . c V r , l T r U r , l , c B r . l . c V r , l i < c η r , l , i + j l , j L T r U r , l , c B r . l . c V r , j + σ 2 .
Problem (31) can be reformulated as
min η r , l , c , V r , l , U r , l , c l L c L l f 1 η r , l , c , U r , l , c , V r , l + l L c L l f 2 η r , l , c , U r , l , c , V r , l
s . t .   f 3 η r , l , c , U r , l , c , V r , l 0 ,
l L T r V r , l P max ,
R a n k V r , l = 1 ,
R a n k U r , l , c = 1 ,
V r , l 0 ,
U r , l , c 0 ,
( 31 e ) ,
where f 1 η r , l , c , U r , l , c , V r , l is as in (34), f 2 η r , l , c , U r , l , c , V r , l is as in Equation (35), and f 3 η r , l , c , U r , l , c , V r , l is as in Equation (36). Problem (33) is still nonconvex, and we use SDR to relax the rank-one constraint and then use a first-order Taylor expansion to approximate the nonconvex part of the objective function with the constraint present.
f 1 η r , l , c , U r , l , c , V r , l = log 2 T r U r , l , c B r . l . c V r , l i < c η r , l , i + j l , j L T r U r , l , c B r . l . c V r , j + σ 2 .
f 2 η r , l , c , U r , l , c , V r , l = log 2 T r U r , l , c B r . l . c V r , l i c η r , l , i + j l , j L T r U r , l , c B r . l . c V r , j + σ 2 .
f 3 η r , l , c , U r , l , c , V r , l = T r U r , l , c B r . l . c V r , l i < c η r , l , i + j l , j L T r U r , l , c B r . l . c V r , j + σ 2 2 R r , l , c min 1 . η r , l , c T r U r , l , c B r . l . c V r , l

4. Solution of the Problem

In this section, we propose a new alternating optimization iterative algorithm to solve the joint optimization problem of power-allocation coefficients, active beamforming, and STAR-RIS passive beamforming.

4.1. Power-Allocation Coefficient Optimization

When fixing the STAR-RIS coefficients and the BS beamforming vectors, problem (33) can be reformulated as
min η r , l , c l L c L l f 1 η r , l , c + l L c L l f 2 η r , l , c
s . t .   f 3 η r , l , c 0 ,
( 31 e ) .
At this point, the constraint is convex, and the first-order Taylor expansion of f 1 η r , l , c is
f 1 η r , l , c f 1 η r , l , c τ + T r η f 1 η r , l , c τ H η r , l , c η r , l , c τ f 1 η r , l , c , η r , l , c τ ,
where η r , l , c τ denotes the value of η r , l , c at the τ - th iteration and, according to the SCA principle, Equation (37a) can be replaced by
min η r , l , c l L c L l f 1 η r , l , c , η r , l , c τ + l L c L l f 2 η r , l , c
Using the penalty-function method, problem (37) can be transformed into an unconstrained optimization problem:
min η r , l , c l L c L l f 1 η r , l , c , η r , l , c τ + l L c L l f 2 η r , l , c + ζ 1 c L l η r , l , c 2 1 + max 0 , f 3 η r , l , c 2
where ζ 1 is the penalty factor. Thus, the gradient-descent method can be used to solve for η r , l , c by using Equation (41).
η r , l , c ( τ + 1 ) = η r , l , c ( τ ) α 3 l L c L l η f 1 τ η r , l , c , η r , l , c τ 1 + α 3 l L c L l η f 2 τ η r , l , c + α 3 ζ 1 2 c L l η r , l , c τ + η τ max 0 , f 3 η r , l , c 2 ,
where α 3 is the step size.

4.2. Active Beamforming Optimization

Fixing the STAR-RIS coefficients and power-allocation coefficients, problem (33) can be formulated as
min V r , l l L c L l f 1 V r , l + l L c L l f 2 V r , l
s . t .   f 3 V r , l 0 ,
( 33 c ) , ( 33 d ) , ( 33 f ) .
Using the SDR relaxation rank-one constraint (33d), the first-order Taylor expansion of f 1 V r , l is
f 1 V r , l f 1 V r , l τ + T r V f 1 V r , l τ H V r , l V r , l τ   f 1 V r , l , V r , l τ ,
where V r , l τ denotes the value of V r , l at the τ - th iteration. Thus, the objective function (42a) can be converted to
min V r , l l L c L l f 1 V r , l , V r , l τ + l L c L l f 2 V r , l
After removing the constraints (42b) and (33c) using the penalty function, problem (42) can be converted to
min V r , l l L c L l f 1 V r , l , V r , l τ + l L c L l f 2 V r , l + ζ 2 max 0 , f 3 V r , l , l L T r V r , l P max 2
s . t .   ( 33 f ) ,
Problem (45) is a semidefinite programming problem (SDP), which is solved by using the BCD [40] method. First, problem (45) is converted into an unconstrained optimization problem:
min V r , l l L c L l f 1 V r , l , V r , l τ + l L c L l f 2 V r , l + ζ 2 max 0 , f 3 V r , l , l L T r V r , l P max 2 ϖ 1 log det V r , l
where ϖ 1 > 0 is the barrier factor. The key to applying the coordinate-descent method to the problem of (46) is the following block determinant:
det V r , l = det S det y x T S 1 x ,
where V r , l = S x x T y , and therefore, problem (46) can be transformed by solving
min x l L c L l f 1 x + l L c L l f 2 x + ζ 2 max 0 , f 3 x , l L C l T x P max 2 ϖ 1 log y x T S 1 x
where f 1 x is (49), f 2 x is (50), and f 3 x is Equation (51).
f 1 x = log 2 w l , c T x i < c η r , l , i + j l , j L w l , j T x + σ 2 .
f 2 x = log 2 w l , c T x i c η r , l , i + j l , j L w l , j T x + σ 2 .
f 3 x = w l , c T x i < c η r , l , i + j l , j L w l , j T x + σ 2 R r , l , c M I N η r , l , c w l , c T x ,
where R r , l , c M I N = 2 R r , l , c min 1 .
The gradient-descent method is used to solve for x according to Equation (52).
x τ + 1 = x τ α 4 l L c L l x f 1 τ x , x τ 1 + l L c L l x f 2 τ x + α 4 ζ 2 x τ max 0 , f 3 x , l L w l T x P max 2 + α 4 2 ϖ 1 y x T S 1 x S 1 x ,
where α 4 is the step size.

4.3. STAR-RIS Coefficient Optimization

Fixing the power-allocation coefficients and BS beamforming vectors, problem (33) can be formulated as
min U r , l , c l L c L l f 1 U r , l , c + l L c L l f 2 U r , l , c
s . t .   f 3 U r , l , c 0 ,
( 33 e ) , ( 33 g ) .
Regarding the relaxation of constraint (33e) using SDR, the first-order Taylor expansion of f 1 U r , l , c is
f 1 U r , l , c f 1 U r , l , c τ + T r U f 1 U r , l , c τ H U r , l , c U r , l , c τ f 1 U r , l , c , U r , l , c τ ,
where U r , l , c τ denotes the value of U r , l , c at the τ - th iteration. The objective function (53a) can be re-expressed as
min U r , l , c l L c L l f 1 U r , l , c , U r , l , c τ + l L c L l f 2 U r , l , c
Convert problem (53) to problem (56) using the penalty function.
min U r , l , c l L c L l f 1 U r , l , c , U r , l , c τ + l L c L l f 2 U r , l , c + ζ 3 max 0 , f 3 U r , l , c 2
s . t .   ( 33 g ) ,
where ζ 3 is the penalty factor. Problem (56) is an SDP problem, and the solution using the BCD method can convert (56) into an unconstrained problem.
min U r , l , c l L c L l f 1 U r , l , c , U r , l , c τ + l L c L l f 2 U r , l , c + ζ 3 max 0 , f 3 U r , l , c 2 ϖ 2 log det U r , l , c
where ϖ 2 > 0 is the barrier factor. The key to applying the coordinate-descent method to the problem of (57) is the following block determinant:
det U r , l , c = det O det y ^ x ^ T O 1 x ^ ,
where U r , l , c = O x ^ x ^ T y ^ ; therefore, problem (57) can be transformed by solving
min x ^ l L c L l f ^ 1 x ^ + l L c L l f ^ 2 x ^ + ζ 3 max 0 , f ^ 3 x ^ 2 ϖ 2 log y ^ x ^ T O 1 x ^
where f ^ 1 x ^ is shown in Equation (60), f ^ 2 x ^ is shown in (61), and f ^ 3 x ^ is shown in Equation (62).
f ^ 1 x ^ = log 2 w ^ l , c T x ^ i < c η r , l , i + j l , j L w ^ l , j T x ^ + σ 2 .
f ^ 2 x ^ = log 2 w ^ l , c T x ^ i c η r , l , i + j l , j L w ^ l , j T x ^ + σ 2 .
f ^ 3 x ^ = w ^ l , c T x ^ i < c η r , l , i + j l , j L w ^ l , j T x ^ + σ 2 R r , l , c M I N η r , l , c w ^ l , c T x ^ .
The gradient-descent method is used to solve for x ^ according to Equation (63).
x ^ τ + 1 = x ^ τ α 5 l L c L l x ^ f ^ 1 τ x ^ , x ^ τ 1 + l L c L l x ^ f ^ 2 τ x ^   + α 5 ζ 3 x ^ τ max 0 , f ^ 3 x ^ 2 + 2 α 5 ϖ 2 y ^ x ^ T O 1 x ^ O 1 x ^ ,
where α 5 is the step size factor.
In Algorithm 1, the complexity of computing the power allocation is O 1 O τ 1 max M 3.5 log 2 1 ε 1 , where τ 1 max is the number of iterations of Algorithm 1, and ε 1 denotes the accuracy of Algorithm 1. In Algorithm 2, the complexity of the active beam formation is calculated as O 2 O τ 2 max N t log N t , where τ 2 max is the number of iterations of the algorithm. In Algorithm 3, the complexity of calculating the STAR-RIS passive beam formation is O 3 O τ 3 max γ r log γ r , where τ 3 max is the number of iterations of the algorithm. Thus, the total complexity of the algorithm is O τ 4 max O 1 + O 2 + O 3 , where τ 4 max denotes the number of outer iterations.
Algorithm 1: Proposed Penalty-Based Iterative Algorithm for Solving Problem (37)
1: Initialize feasible point η r , l , c ( 0 ) ,   and   penalty   factor   ζ 1 .
2: repeat:
3:  Set the iteration index τ = 0 .
4:  repeat:
5:   Use   f 1 η r , l , c , η r , l , c τ   to   approximate   f 1 η r , l , c .
6:   Update   η r , l , c ( τ + 1 ) by solving problem (40) using (41).
7:   τ = τ + 1 .
8:    until the fraction of objective function values decreases below a predefined threshold or the maximum number of internal iterations is reached.
9:  Update   η r , l , c ( 0 ) .
10:   Update   penalty   factor   ζ 1 .
11: until the constraint violation is below a predefined threshold.
Algorithm 2: Proposed Penalty-Based Iterative Algorithm for Solving Problem (42)
1: Initialize the feasible point a V r , l ( 0 ) ,   penalty   factor   ζ 2 .
2: repeat:
3:  Set   the   iteration   index   τ = 0 .
4:  repeat:
5:   For   a   given   V r , l ( τ ) ,   set   the   barrier   factor   ϖ 1 , and use BCD to convert the solution of problem (46) to problem (48).
6:   Update   x τ + 1 by (52).
7:   Update   V r , l ( τ + 1 ) with the obtained optimal solution.
8:   τ = τ + 1 .
9:  until the fraction of objective function values decreases below a predefined threshold or the maximum number of internal iterations is reached.
10 :   Update   V r , l ( 0 ) .
12: until the constraint violation is below a predefined threshold.
Algorithm 3: Proposed Penalty-Based Iterative Algorithm for Solving Problem (53)
1: Initialize feasible point U r , l , c ( 0 ) ,   and   penalty   factor   ζ 3 .
2: repeat:
3:  Set   the   iteration   index   τ = 0 .
4:  repeat:
5:   For   a   given   U r , l , c ( τ ) ,   set   the   barrier   factor   ϖ 2 and convert problem (57) to problem (59).
6:   Update   x ^ τ + 1 using (63).
7:   Update   U r , l , c ( τ + 1 ) with the obtained optimal solution.
8:   τ = τ + 1 .
9:  until the fraction of objective function values decreases below a predefined threshold or the maximum number of internal iterations is reached.
10 :   Update   U r , l , c ( 0 ) .
12: until the constraint violation is below a predefined threshold.

5. Numerical Results

In this section, numerical simulations are performed to evaluate the performance of the proposed algorithm. The specific simulation scenario is shown in Figure 5. Without loss of generality, it is assumed that the BS and the STAR-RIS are located at (0, 0, 50) and (0, 40, 50), respectively. The STAR-RIS–NOMA system has M = 12 users. The users are divided into L = 4 clusters with cluster center coordinates of (0, 20, 0), (0, 30, 15), (0, 50, 0), and (0, 60, −15), respectively, and the users are randomly distributed in clusters corresponding to a radius of 6. It is assumed that STAR-RIS is equipped with N r i s = N R I S h N r i s v elements, where N R I S h = 5 denotes the number of elements along the horizontal direction, N R I S v = 5 represents the number of elements along the vertical path, and STAR-RIS elements are divided into γ = 4 clusters. The narrowband quasi-static fading channel from the BS to the STAR-RIS and from the STAR-RIS to the user is modeled as a Rician-fading channel as follows:
H = ς 0 ξ B R ω B R B B R B B R + 1 H L o s + 1 B B R + 1 H N L o s ,
g l , c = ς 0 ξ R U , l , c ω R U B R U B R U + 1 g l , c L o s + 1 B R U + 1 g l , c N L o s ,
where ξ B R and ξ R U , l , c denote the distance between the BS and the STAR-RIS and the distance between the STAR-RIS and the user, respectively; ς o denotes the path loss at a reference distance of 1 m; ω B R and ω R U denote the corresponding path loss index; B B R and B R U denote the Rician factors; and H L o s and g l , c L o s are the determined line-of-sight (LoS) components. H N L o s and g l , c N L o s are the random non-line-of-sight (NLoS) components modeled using Rician fading. In order not to lose generality, assume that users have the same QoS requirements and set R l , c min = 0.1   bits s Hz . The system parameters used are shown in Table 1.
Figure 6 investigates the convergence of the power-allocation algorithm. With P max = 35 dBm , N r i s = 25 , and N t = 5 , the user-achievable sum rate of the algorithm converges in about 11 iterations as the number of iterations increases.
Figure 7 shows the relationship between the user-achievable sum rate and the number of RIS elements. A total of four cases are considered. The first case is the use of the STAR-RIS in the NOMA system, where the number of elements of STAR-RIS is N r i s = N . In the second case, two conventional RISs are used in the system: one for reflection and the other for transmission, with N elements per RIS. The third case also uses conventional RISs, but each RIS has N 2 elements. The fourth type is the NOMA system without RIS assistance. The results show that the performance of the STAR-RIS is higher than that of the conventional RIS for the same number of RIS elements. The higher version of the dual RIS than the STAR-RIS is because the dual RIS has N reflected and N transmitted elements at the same moment, while the transmitted/reflected elements of the STAR-RIS add up to a total of N .
Figure 8 shows the user-achievable sum rate versus the BS transmit antenna in the case of conventional RIS-assisted NOMA and STAR-RIS-assisted NOMA communications. The experimental results show that the use of STAR-RIS-assisted communication is superior to a conventional RIS when the number of BS antennas is the same, given the same number of RIS elements. It can be observed that using RIS in NOMA can significantly improve the system’s performance. Furthermore, as the number of antennas at the BS increases, the performance of the RIS–NOMA system also increases.
Figure 9 shows the user-achievable sum rate versus the total transmitted power Pmax during conventional RIS and STAR-RIS-assisted NOMA communications. The user-achievable sum rate increases with the total transmission power. RIS-assisted communications significantly improve user achievability and rates over NOMA systems without RIS assistance. It was observed that the performance difference between the system using two conventional RISs and the STAR-RIS system was relatively small.
Figure 10 considers the relationship between the user-achievable sum rate and the minimum user rate for three cases: a STAR-RIS system with N elements, a RIS system with 2 N elements, and a RIS system with N elements. The results show that the user-achievable sum rate decreases as the minimum user rate increases. In all three cases, the RIS system with 2 N components has the highest user-achievable sum rate, and the STAR-RIS system with N elements performs better than the RIS system with N elements.
Figure 11 shows how the proposed method in this paper compares with four other methods or schemes. The first scheme is where no RIS is configured in the system. The second scheme is where each cell of the STAR-RIS is randomly configured with a phase. The third scheme uses a pairwise-rising method [41] to determine the phase shift of the STAR-RIS. The fourth scheme is to use the Riemann conjugate-gradient method [42] to solve the phase-shift matrix of the STAR-RIS. As can be seen from the results in Figure 11, when the number of RIS elements is 60, the user-achievable sum rate of the proposed method in this paper is 13.1 bits/s/Hz. The achievable sum rate of the random phase-shift method is only 10.8 bits/s/Hz. The achievable sum rates using the pairwise-ascent method and the Riemann conjugate-gradient method are 12.8 bits/s/Hz and 12 bits/s/Hz, respectively.
Figure 12 shows the variation of achievable users and rates with the number of antennas at the BS. When the number of BS transmitting antennas is 10, the user-achievable sum rate of the random phase shift is 7.3 bits/s/Hz. The user-achievable sum rates of the pairwise-ascent method and Riemann conjugate-gradient method are 10.8 bits/s/Hz and 11.5 bits/s/Hz, respectively. The achievable sum rate of the proposed method is 12.6 bits/s/Hz. Among the five schemes, the worst performance is achieved without RIS assistance, followed by the random RIS phase shift. The Riemann conjugate-gradient method is closer to the case of the pairwise-ascent method, and the proposed method in this paper has the best performance. Therefore, this shows the effectiveness of the proposed algorithm.

6. Conclusions

In this paper, we study the downlink transmission system of the STAR-RIS–NOMA. A new alternate optimization method is proposed for the MS protocol to jointly optimize the power-allocation coefficient of the user, active beamforming at the BS, and passive beamforming of the STAR-RIS so that the achievable rate of the user is maximized. The channel is first feature-extracted, and the user and STAR-RIS elements are clustered according to the extracted features to reduce the solution difficulty. The SDR technique and penalty function are used to solve the nonconvex factors present in the problem. The effectiveness of the proposed algorithm can be seen in the simulation results.
This paper considers the STAR-RIS in an ideal state, where the transmission and reflection phase-shift coefficients can be adjusted independently. However, in real applications, coupled transmission and reflection phase-shift coefficients are more common, and the related research problems are more complex. Therefore, the design of coupled STAR-RIS phase-shift coefficients is a crucial research direction for the future.

Author Contributions

Conceptualization, X.Z. and H.W.; methodology, H.W.; software, H.W.; validation, X.Z. and J.C.; formal analysis, X.Z. and H.W.; investigation, H.W.; resources, J.C.; data curation, J.C.; writing—original draft preparation, H.W.; writing—review and editing, H.W. and X.Z.; visualization, J.C.; supervision, X.Z.; project administration, X.Z. and J.C.; funding acquisition, X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Shanghai Capacity Building Projects in Local Institutions under Grants 19070502900 and Science and Technology Commission of Shanghai Municipality under Grants 22142201900.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tariq, F.; Khandaker, M.R.A.; Wong, K.-K.; Imran, M.A.; Bennis, M.; Debbah, M. A Speculative Study on 6G. IEEE Wirel. Commun. 2020, 27, 118–125. [Google Scholar] [CrossRef]
  2. di Renzo, M.; Zappone, A.; Debbah, M.; Alouini, M.S.; Yuen, C.; de Rosny, J.; Tretyakov, S. Smart Radio Environments Empowered by Reconfigurable Intelligent Surfaces: How It Works, State of Research, and The Road Ahead. IEEE J. Sel. Areas Commun. 2020, 38, 2450–2525. [Google Scholar] [CrossRef]
  3. Huang, C.; Hu, S.; Alexandropoulos, G.C.; Zappone, A.; Yuen, C.; Zhang, R.; di Renzo, M.; Debbah, M. Holographic MIMO Surfaces for 6G Wireless Networks: Opportunities, Challenges, and Trends. IEEE Wirel. Commun. 2020, 27, 118–125. [Google Scholar] [CrossRef]
  4. Huang, C.; Zappone, A.; Alexandropoulos, G.C.; Debbah, M.; Yuen, C. Reconfigurable Intelligent Surfaces for Energy Efficiency in Wireless Communication. IEEE Trans. Wirel. Commun. 2019, 18, 4157–4170. [Google Scholar] [CrossRef] [Green Version]
  5. Akyildiz, I.F.; Han, C.; Nie, S. Combating the Distance Problem in the Millimeter Wave and Terahertz Frequency Bands. IEEE Commun. Mag. 2018, 56, 102–108. [Google Scholar] [CrossRef] [Green Version]
  6. Huang, C.; Mo, R.; Yuen, C. Reconfigurable Intelligent Surface Assisted Multiuser MISO Systems Exploiting Deep Reinforcement Learning. IEEE J. Sel. Areas Commun. 2020, 38, 1839–1850. [Google Scholar] [CrossRef]
  7. Liu, Y.; Mu, X.; Xu, J.; Schober, R.; Hao, Y.; Poor, H.V.; Hanzo, L. STAR: Simultaneous Transmission and Reflection for 360° Coverage by Intelligent Surfaces. IEEE Wirel. Commun. 2021, 28, 102–109. [Google Scholar] [CrossRef]
  8. Xu, J.; Liu, Y.; Mu, X.; Zhou, J.T.; Song, L.; Poor, H.V.; Hanzo, L. Simultaneously Transmitting and Reflecting Intelligent Omni-Surfaces: Modeling and Implementation. IEEE Veh. Technol. Mag. 2022, 17, 46–54. [Google Scholar] [CrossRef]
  9. Jian, M.; Alexandropoulos, G.C.; Basar, E.; Huang, C.; Liu, R.; Liu, Y.; Yuen, C. Reconfigurable intelligent surfaces for wireless communications: Overview of hardware designs, channel models, and estimation techniques. Intell. Converg. Netw. 2022, 3, 1–32. [Google Scholar] [CrossRef]
  10. Ding, Z.; Liu, Y.; Choi, J.; Sun, Q.; Elkashlan, M.; Chih-Lin, I.; Poor, H.V. Application of Non-Orthogonal Multiple Access in LTE and 5G Networks. IEEE Commun. Mag. 2017, 55, 185–191. [Google Scholar] [CrossRef] [Green Version]
  11. Zheng, B.; Wu, Q.; Zhang, R. Intelligent Reflecting Surface-Assisted Multiple Access with User Pairing: NOMA or OMA? IEEE Commun. Lett. 2020, 24, 753–757. [Google Scholar] [CrossRef] [Green Version]
  12. Mu, X.; Liu, Y.; Guo, L.; Lin, J.; Al-Dhahir, N. Exploiting Intelligent Reflecting Surfaces in Multi-Antenna Aided NOMA Systems. arXiv 2019, arXiv:1910.13636. [Google Scholar]
  13. Liu, F.; Mähönen, P.; Petrova, M. Proportional fairness-based power allocation and user set selection for downlink NOMA systems. In Proceedings of the 2016 IEEE International Conference on Communications (ICC), Kuala Lumpur, Malaysia, 22–27 May 2016; pp. 1–6. [Google Scholar] [CrossRef]
  14. Lin, Z.; Lin, M.; Wang, J.-B.; de Cola, T.; Wang, J. Joint Beamforming and Power Allocation for Satellite-Terrestrial Integrated Networks with Non-Orthogonal Multiple Access. IEEE J. Sel. Top. Signal Process. 2019, 13, 657–670. [Google Scholar] [CrossRef] [Green Version]
  15. Ding, Z.; Yang, Z.; Fan, P.; Poor, H.V. On the Performance of Non-Orthogonal Multiple Access in 5G Systems with Randomly Deployed Users. IEEE Signal Process. Lett. 2014, 21, 1501–1505. [Google Scholar] [CrossRef] [Green Version]
  16. Balasubramanya, N.M.; Gupta, A.; Sellathurai, M. Combining Code-Domain and Power-Domain NOMA for Supporting Higher Number of Users. In Proceedings of the 2018 IEEE Global Communications Conference (GLOBECOM), Abu Dhabi, United Arab Emirates, 9–13 December 2018; pp. 1–6. [Google Scholar] [CrossRef]
  17. Ali, S.; Hossain, E.; Kim, D.I. Non-Orthogonal Multiple Access (NOMA) for Downlink Multiuser MIMO Systems: User Clustering, Beamforming, and Power Allocation. IEEE Access 2017, 5, 565–577. [Google Scholar] [CrossRef]
  18. Ding, Z.; Schober, R.; Poor, H.V. Unveiling the Importance of SIC in NOMA Systems—Part 1: State of the Art and Recent Findings. IEEE Commun. Lett. 2020, 24, 2373–2377. [Google Scholar] [CrossRef]
  19. Liu, Y.; Liu, X.; Mu, X.; Hou, T.; Xu, J.; di Renzo, M.; Al-Dhahir, N. Reconfigurable Intelligent Surfaces: Principles and Opportunities. IEEE Commun. Surv. Tutor. 2021, 23, 1546–1577. [Google Scholar] [CrossRef]
  20. de Sena, A.S.; Carrillo, D.; Fang, F.; Nardelli, P.H.J.; da Costa, D.B.; Dias, U.S.; Ding, Z.; Papadias, C.B.; Saad, W. What Role Do Intelligent Reflecting Surfaces Play in Multi-Antenna Non-Orthogonal Multiple Access? IEEE Wirel. Commun. 2020, 27, 24–31. [Google Scholar] [CrossRef]
  21. Lin, Z.; Niu, H.; An, K.; Wang, Y.; Zheng, G.; Chatzinotas, S.; Hu, Y. Refracting RIS-Aided Hybrid Satellite-Terrestrial Relay Networks: Joint Beamforming Design and Optimization. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 3717–3724. [Google Scholar] [CrossRef]
  22. Lin, Z.; Lin, M.; Champagne, B.; Zhu, W.-P.; Al-Dhahir, N. Secrecy-Energy Efficient Hybrid Beamforming for Satellite-Terrestrial Integrated Networks. IEEE Trans. Commun. 2021, 69, 6345–6360. [Google Scholar] [CrossRef]
  23. Jung, M.; Saad, W.; Jang, Y.; Kong, G.; Choi, S. Performance Analysis of Large Intelligent Surfaces (LISs): Asymptotic Data Rate and Channel Hardening Effects. IEEE Trans. Wirel. Commun. 2020, 19, 2052–2065. [Google Scholar] [CrossRef] [Green Version]
  24. Niu, H.; Lin, Z.; An, K.; Liang, X.; Hu, Y.; Li, D.; Zheng, G. Active RIS-Assisted Secure Transmission for Cognitive Satellite Terrestrial Networks. IEEE Trans. Veh. Technol. 2023, 72, 2609–2614. [Google Scholar] [CrossRef]
  25. Hong, S.; Pan, C.; Ren, H.; Wang, K.; Nallanathan, A. Artificial-Noise-Aided Secure MIMO Wireless Communications via Intelligent Reflecting Surface. IEEE Trans. Commun. 2020, 68, 7851–7866. [Google Scholar] [CrossRef]
  26. Wang, P.; Fang, J.; Duan, H.; Li, H. Compressed Channel Estimation for Intelligent Reflecting Surface-Assisted Millimeter Wave Systems. IEEE Signal Process. Lett. 2020, 27, 905–909. [Google Scholar] [CrossRef]
  27. Yu, X.; Xu, D.; Schober, R. Optimal Beamforming for MISO Communications via Intelligent Reflecting Surfaces. In Proceedings of the 2020 IEEE 21st International Workshop on Signal Processing Advances in Wireless Communications (SPAWC), Atlanta, GA, USA, 26–29 May 2020; pp. 1–5. [Google Scholar] [CrossRef]
  28. Wu, Q.; Zhang, R. Intelligent Reflecting Surface Enhanced Wireless Network via Joint Active and Passive Beamforming. IEEE Trans. Wirel. Commun. 2019, 18, 5394–5409. [Google Scholar] [CrossRef] [Green Version]
  29. Yu, X.; Xu, D.; Schober, R. MISO Wireless Communication Systems via Intelligent Reflecting Surfaces: (Invited Paper). In Proceedings of the 2019 IEEE/CIC International Conference on Communications in China (ICCC), Changchun, China, 11–13 August 2019; pp. 735–740. [Google Scholar] [CrossRef]
  30. Yang, G.; Xu, X.; Liang, Y.-C. Intelligent Reflecting Surface Assisted Non-Orthogonal Multiple Access. In Proceedings of the 2020 IEEE Wireless Communications and Networking Conference (WCNC), Seoul, Republic of Korea, 25–28 May 2020; pp. 1–6. [Google Scholar] [CrossRef]
  31. Kumaravelu, V.B.; Imoize, A.L.; Soria, F.R.C.; Velmurugan, P.G.S.; Thiruvengadam, S.J.; Murugadass, A.; Gudla, V.V. Outage Probability Analysis and Transmit Power Optimization for Blind-Reconfigurable Intelligent Surface-Assisted Non-Orthogonal Multiple Access Uplink. Sustainability 2022, 14, 13188. [Google Scholar] [CrossRef]
  32. Wu, C.; Liu, Y.; Mu, X.; Gu, X.; Dobre, O.A. Coverage Characterization of STAR-RIS Networks: NOMA and OMA. IEEE Commun. Lett. 2021, 25, 3036–3040. [Google Scholar] [CrossRef]
  33. Mu, X.; Liu, Y.; Guo, L.; Lin, J.; Schober, R. Simultaneously Transmitting and Reflecting (STAR) RIS Aided Wireless Communications. IEEE Trans. Wirel. Commun. 2022, 21, 3083–3098. [Google Scholar] [CrossRef]
  34. Zhu, J.; Gao, P.; Chen, G.; Xiao, P.; Quddus, A. Index Modulation for STAR-RIS Assisted NOMA System. IEEE Commun. Lett. 2023, 27, 716–720. [Google Scholar] [CrossRef]
  35. Karim, F.; Singh, S.K.; Singh, K.; Flanagan, M.F. STAR-RIS aided Full Duplex Communication System: Performance Analysis. In Proceedings of the GLOBECOM 2022—2022 IEEE Global Communications Conference, Rio de Janeiro, Brazil, 4–8 December 2022; pp. 3114–3119. [Google Scholar] [CrossRef]
  36. Chen, J.; Yu, X. Ergodic Rate Analysis and Phase Design of STAR-RIS Aided NOMA With Statistical CSI. IEEE Commun. Lett. 2022, 26, 2889–2893. [Google Scholar] [CrossRef]
  37. Ma, H.; Wang, H.; Li, H.; Feng, Y. Transmit Power Minimization for STAR-RIS-Empowered Uplink NOMA System. IEEE Wirel. Commun. Lett. 2022, 11, 2430–2434. [Google Scholar] [CrossRef]
  38. Fang, F.; Wu, B.; Fu, S.; Ding, Z.; Wang, X. Energy-Efficient Design of STAR-RIS Aided MIMO-NOMA Networks. IEEE Trans. Commun. 2023, 71, 498–511. [Google Scholar] [CrossRef]
  39. McInnes, L.; Healy, J. UMAP: Uniform Manifold Approximation and Projection for Dimension Reduction. arXiv 2018, arXiv:1802.03426. [Google Scholar]
  40. Waldspurger, I.; d’Aspremont, A.; Mallat, S. Phase recovery, MaxCut and complex semidefinite programming. Math. Program. 2012, 149, 47–81. [Google Scholar] [CrossRef] [Green Version]
  41. Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers; Now Foundations and Trends: Norwell, MA, USA, 2011. [Google Scholar]
  42. Huang, W. Optimization Algorithms on Riemannian Manifolds with Applications. Ph.D. Thesis, Florida State University, Tallahassee, FL, USA, 2013. [Google Scholar]
Figure 1. STAR-RIS-assisted NOMA downlink communication system.
Figure 1. STAR-RIS-assisted NOMA downlink communication system.
Sensors 23 05504 g001
Figure 2. UMAP dimensionality reduction process.
Figure 2. UMAP dimensionality reduction process.
Sensors 23 05504 g002
Figure 3. User-clustering process.
Figure 3. User-clustering process.
Sensors 23 05504 g003
Figure 4. STAR-RIS-element clustering process.
Figure 4. STAR-RIS-element clustering process.
Sensors 23 05504 g004
Figure 5. STAR-RIS simulation scenario.
Figure 5. STAR-RIS simulation scenario.
Sensors 23 05504 g005
Figure 6. Convergence of the power-allocation algorithm.
Figure 6. Convergence of the power-allocation algorithm.
Sensors 23 05504 g006
Figure 7. Achievable sum rate relative to the number of RIS elements for P max = 35   dBm and N t = 5 .
Figure 7. Achievable sum rate relative to the number of RIS elements for P max = 35   dBm and N t = 5 .
Sensors 23 05504 g007
Figure 8. Achievable sum rate relative to the number of BS elements for P max = 35   dBm and N r i s = 25 .
Figure 8. Achievable sum rate relative to the number of BS elements for P max = 35   dBm and N r i s = 25 .
Sensors 23 05504 g008
Figure 9. Achievable sum rate relative to total transmitting power for N r i s = 25 and N t = 5 .
Figure 9. Achievable sum rate relative to total transmitting power for N r i s = 25 and N t = 5 .
Sensors 23 05504 g009
Figure 10. Achievable sum rate versus minimum user rate for P max = 35   dBm , N r i s = 25 and N t = 5 .
Figure 10. Achievable sum rate versus minimum user rate for P max = 35   dBm , N r i s = 25 and N t = 5 .
Sensors 23 05504 g010
Figure 11. Achievable sum rate relative to the number of RIS elements under different methods.
Figure 11. Achievable sum rate relative to the number of RIS elements under different methods.
Sensors 23 05504 g011
Figure 12. Achievable sum rate relative to the number of BS antennas under different methods.
Figure 12. Achievable sum rate relative to the number of BS antennas under different methods.
Sensors 23 05504 g012
Table 1. System Parameters.
Table 1. System Parameters.
ParameterParameter ExplanationValue
B B R Rician factor for BS to STAR-RIS channels3 dB
B R U Rician factor of STAR-RIS to user channel3 dB
ς o Path loss at a reference distance of 1 m−30 dB
ω B R Path loss index for BS to STAR-RIS channels2.2
ω R U Path Loss Index for Path Loss Index2.2
ζ 1 , ζ 2 , ζ 3 , ζ 4 Initialize penalty factor10−4
τ max Maximum number of inner-layer iterations25
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhou, X.; Wang, H.; Chen, J. Joint Design of a Simultaneous Reflection and Transmission RIS in Mode-Switching Mode to Assist NOMA Systems. Sensors 2023, 23, 5504. https://doi.org/10.3390/s23125504

AMA Style

Zhou X, Wang H, Chen J. Joint Design of a Simultaneous Reflection and Transmission RIS in Mode-Switching Mode to Assist NOMA Systems. Sensors. 2023; 23(12):5504. https://doi.org/10.3390/s23125504

Chicago/Turabian Style

Zhou, Xiaoping, Hanqi Wang, and Jiajia Chen. 2023. "Joint Design of a Simultaneous Reflection and Transmission RIS in Mode-Switching Mode to Assist NOMA Systems" Sensors 23, no. 12: 5504. https://doi.org/10.3390/s23125504

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