Next Article in Journal
SPRi Biosensor for Simultaneous Determination of HIF-1α, Angiopoietin-2, and Interleukin-1β in Blood Plasma
Previous Article in Journal
IoT-Based Assessment of a Driver’s Stress Level
Previous Article in Special Issue
AI Optimization-Based Heterogeneous Approach for Green Next-Generation Communication Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ergodic Rate Analysis of Simultaneous Transmitting and Reflecting Reconfigurable Intelligent Surface-Assisted Rate-Splitting Multiple Access Systems Based on Discrete Phase Shifts

College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(17), 5480; https://doi.org/10.3390/s24175480
Submission received: 15 July 2024 / Revised: 19 August 2024 / Accepted: 20 August 2024 / Published: 23 August 2024
(This article belongs to the Special Issue Energy-Efficient Communication Networks and Systems: 2nd Edition)

Abstract

:
In this paper, we combine simultaneous transmitting and reflecting reconfigurable intelligent surface (STAR-RIS) with rate-splitting multiple access (RSMA) technology and investigate the ergodic rate performance of an STAR-assisted RSMA system. Considering the discrete phase shifts of the STAR-RIS in practice, the downlink performance of STAR-RIS-assisted RSMA with discrete phase shifts is compared to that with continuous phase shifts. Firstly, the cumulative distribution function of signal-to-interference-plus-noise ratio (SINR) of users is analyzed. Then, the total ergodic rate of the system and its approximate closed-form solution are, respectively, derived based on the cumulative distribution function of users. The simulation results validate the effectiveness of the theoretical analysis, showing good agreement between the derived theoretical ergodic rate and the corresponding simulations. Although the system performance with discrete phase shifts is inferior to that with continuous phase shifts due to quantization errors, the performance of the continuous phase shift system is well approximated when the quantization bit of the phase shift system reaches 3 in the simulations. Additionally, the impact of the number of STAR-RIS elements on the system’s performance is analyzed.

1. Introduction

The development of communication technology has always been a key factor driving progress in modern society. From early telegraphs and telephones to modern mobile communication and the internet, each technological leap has greatly changed people’s lifestyles and the way that information is disseminated. Since the advent of traditional wireless communication, people have always believed that the propagation medium is a physical entity of random behavior between transmitters and receivers [1]. The interaction of uncontrollable electromagnetic waves with objects between transmitters and receivers degrades the quality of received signals. In recent years, with the development of smart surface technology, the field of communication has again undergone revolutionary changes. Reconfigurable intelligent surfaces (RISs) are an emerging communication technology that dynamically adjust the propagation characteristics of signals by deploying a large number of passive reflecting elements on the wireless signal propagation path [2]. This optimization of the signal transmission path enhances communication quality, extends signal coverage, and reduces signal interference. The potential advantages of RIS technology lie in its flexibility and programmability.
In addition, notable potential features of the RIS include energy efficiency and cost effectiveness. Its reflecting elements only passively reflect the impinging signals, without requiring any transmit radio frequency (RF) chains, and thus can be operated with orders-of-magnitude-lower energy costs compared to traditional relays and amplifiers [3,4]. The RIS is able to optimize signal transmission paths by passively reflecting and actively adjusting the electromagnetic wave propagation environment, thereby reducing signal attenuation and interference in conventional systems. Moreover, the flexible deployment features of the RIS enable it to work in key locations within the network, to cover blind or weak signal areas. The RIS’s capabilities allow communication systems to achieve the same coverage and data rates at lower transmit power, directly reducing energy consumption [5].
With further research, scholars have proposed a simultaneous transmitting and reflecting reconfigurable intelligent surface (STAR-RIS) as a branch of RIS technology. This surface can switch between transmitting and reflecting, expanding the functionality of intelligent surfaces. STAR-RIS not only reflects signals to enhance signal strength but also transmits signals, allowing signals to pass through the surface for further propagation. This feature is particularly useful for improving signal penetration and coverage range, especially in densely built urban environments, where signals may attenuate due to obstacles like buildings [6]. By deploying an STAR-RIS, it is possible to effectively bypass obstacles and achieve continuous signal coverage.
Mobile communication technology has undergone several technological innovations since its inception, with multiple access technology being key to enabling multi-user communication. Orthogonal multiple access (OMA) was a widely used technology in early mobile communication systems; it ensured signal orthogonality by assigning unique time or frequency resources to each user to avoid interference. However, with an increasing number of users, spectrum scarcity becomes more pronounced. To address this issue, non-orthogonal multiple access (NOMA) has emerged. NOMA allows multiple users to share the same time-frequency resources, distinguishing users within a power domain to achieve multi-user access [7,8]. NOMA leverages the different channels between users and employs successive interference cancelation (SIC) at the receiver to decode interference completely, thereby reducing interference among multiple users. This technology enhances spectrum efficiency, enabling more users to communicate using the same resources. Compared to NOMA, which is suitable for scenarios with strong multi-user interference, rate-splitting multiple access (RSMA) has broader applicability. RSMA further optimizes the way multiple users share resources by splitting each user’s data streams into multiple sub-streams; independently coding, modulating, and allocating power to these sub-streams; and, after that, combining all the sub-streams into a composite signal for transmission [9]. This technology not only improves spectrum efficiency but also enhances system flexibility and fairness, enabling dynamic resource allocation based on users’ real-time requirements.
In [10,11,12], a detailed study on RSMA technology was provided. The principles and communication models of RSMA were introduced in reference [10], which compared and summarized the advantages and disadvantages of several multiple access technologies including OMA, NOMA, and RSMA in communication systems. It pointed out that the signal interference issue in 5G and 6G had the greatest impact on key technical indicators in communication systems. The relationship between RSMA, spatial division multiple access (SDMA), and NOMA was clarified in [11]. It proposed that RSMA acted as a bridge between SDMA and NOMA. SDMA treated interference completely as noise, whereas NOMA decoded interference entirely. Meanwhile, RSMA decoded part of the interference and treated part as noise. The combination of cooperative user relaying (CUR) with RSMA in wireless communication networks was explored in reference [12], which discussed the core challenges of CUR and RSMA and looked ahead to the prospects of their combination. Due to the promising future of RSMA and RIS, research on RIS and RSMA has attracted widespread attention. Reference [13] studied a wired and tethered unmanned aerial vehicle (UAV) communication system assisted by RIS and optimized the UAV position, RIS phase, and RSMA parameters to maximize the weighted sum-rate of users. In [14], a quantum machine learning (QML) method was proposed to maximize the energy efficiency of RIS-assisted RSMA communications. Compared to conventional optimization, QML achieved the same optimization performance with lower complexity. A hybrid airborne full-duplex (FD) relay system with a reconfigurable intelligent surface on the UAV was considered in [15]. Finite block length encoding RSMA was used to improve spectral efficiency and reduce latency. Reference [15] also employed the alternating optimization (AO) method to handle the non-convex problem of the weighted sum-rate maximization. The ergodic capacity in a multi-user RSMA-RIS communication system was investigated in [16], which compared RSMA-RIS with NOMA-RIS and showed that RSMA-RIS outperformed in both single-user ergodic capacity and multi-user average ergodic capacity.
Combining STAR-RIS with RSMA, for green communication and expanded communication coverage, has been studied in various papers. Closed-form solutions for user outage probability under Energy Splitting (ES) and Mode Switching (MS) STAR-RIS protocols were provided in [17]. Reference [17] also demonstrated the potential of RSMA and STAR-RIS technologies to enhance next-generation wireless communication networks. Reference [18] investigated spatially correlated Rayleigh channels in RSMA-assisted STAR-RIS communication, deriving new expressions for joint moments of spatially correlated Rayleigh channels. In [19], covert communication issues in an STAR-RIS-assisted RSMA system were explored. It highlighted the impact of STAR-RIS’s element reflectivity on covert performance relative to transmit power. When the transmit power was high enough, the detection error probability was a monotonically increasing function of the reflection coefficient. The potential of combining STAR-RIS with RSMA for assisting multi-device communication was underscored in [20], which achieved significant performance improvements through parameter adjustments. Although there exists a large number of studies on STAR-RIS and RSMA communication, they lack consideration of the discrete phase shifts of STAR-RIS, which is not realistic. Based on this, the main contributions of this paper are as follows.
(1)
We primarily consider a multi-user downlink RSMA communication system assisted by STAR-RIS. Firstly, we cluster the multiple users and ignore the interference between different user clusters. Secondly, considering that users are uniformly distributed within the range centered on the STAR-RIS, we quantify the phase shifts of the STAR-RIS. With this result, we further derive the theoretical values of the system ergodic rate under both the discrete and continuous phase shifts of the STAR-RIS for the performance evaluation.
(2)
Through Monte Carlo simulations, we verify the correctness of the theoretical derivations and demonstrate that the system’s ergodic rate with the discrete phase shifts of the STAR-RIS is inferior to that with the continuous phase shifts. However, increasing the number of quantization bits can make the performance of the discrete phase shifts approach that of the continuous phase shifts. Additionally, we show the impact of the number of STAR-RIS elements, the path loss exponent, and the different modes of STAR-RIS on the system performance.

2. System Model

2.1. Channel Model

In the STAR-RIS-assisted downlink RSMA system model, there is one wireless AP, K STAR-RIS units, and 2K users. As shown in Figure 1, each STAR-RIS consists of N elements, and both the AP and the users are equipped with only one antenna each. The 2K users are divided into K user clusters, each containing one STAR-RIS, with the transmitting user and reflecting user located on opposite sides of the STAR-RIS. It is assumed that the direct channels between the AP and the users are blocked by obstacles [21,22].
In this paper, it is assumed that there is a block between individual user clusters, so that the mutual influence between user clusters can be ignored. The channel between the AP and STAR-RIS is modeled as a Rayleigh channel, denoted by h k = ρ 0 / d A P , k α h ˜ k . In the statement, ρ 0 denotes the path loss at a reference unit distance, α is the path loss exponent, d A P , k refers to the distance between the AP and the k-th STAR-RIS, and the n-th element in h ˜ k N × 1 satisfies h ˜ k , n C N 0 , 1 . Due to the random distribution of the transmitting user T and the reflecting user R within a circular area, with the center point C, the STAR-RIS is located meters above point C. It is also assumed that the STAR-RIS can obtain the channel state information of the users. The STAR-RIS to the corresponding user is also modeled as a Rayleigh channel, denoted as g k , m = ρ 0 / d k , m 2 + L k 2 α / 2 g ˜ k , m , where d k , m denotes the distance from user m to the center point C, and g ˜ k , m , n is the n-th element in g ˜ k , m and satisfies g ˜ k , m , n C N 0 , 1 . m t , r denotes the type of user in the k-th user cluster, m = t denotes the transmitting user, and m = r is the reflecting user.
Since the user clusters are independent of each other, this paper specifically analyzes a single user cluster. It is assumed that the reflecting users are near users with better channel gains, and the transmitting users are far users. The positions of the users are modeled using a homogeneous Poisson point process [23]. The reflecting users are distributed within a circle of radius R 1 , with polar coordinate d k , r , ϑ k , r . The transmitting users are distributed within an annular region between radius R 1 and radius R 2 , where R 1 < R 2 , with polar coordinate d k , t , ϑ k , t . The probability density function for the distance from the users to the center point and the angle is as follows [23].
f d k , r x = x π x 2 π R 1 2 = 2 x R 1 2 0 < x < R 1 ,
f d k , t x = x π x 2 R 1 2 π R 2 2 R 1 2 = 2 x R 2 2 R 1 2 R 1 < x < R 2 ,
f ϑ k , l x = 1 π   - π 2 < x < π 2 .

2.2. STAR-RIS-Assisted RSMA Downlink Communication Model

Based on the working principle of the STAR-RIS, the signal is transmitted to the user T through the transmission module of the STAR-RIS and to the user R through the reflection module. The STAR-RIS works in the ES mode. Assuming that all elements on the STAR-RIS have the same amplitude coefficients, the transmission factor and reflection factor matrices for the k-th user cluster are, respectively, represented as
Θ k , t = β k , t d i a g e j θ k , t , 1 , e j θ k , t , 2 , , e j θ k , t , N ,
Θ k , r = β k , r d i a g e j θ k , r , 1 , e j θ k , r , 2 , , e j θ k , r , N ,
where β k , t and β k , r represent the transmission and reflection amplitude coefficients in the k-th user cluster, respectively, and satisfy β k , t + β k , r 1 , β k , t , β k , r 0 , 1 [22]. The phase shift is θ n k , t , θ n k , r 0 , 2 π and n 1 , 2 , N .
For the RSMA-assisted communication system, all the common messages are encoded into a single common stream which must be decoded by all users, and all the private messages are encoded into individual private streams which will be decoded by their respective receivers. By dividing the information into private and common parts, RSMA can make more efficient use of the available spectrum, especially if there is interference between users. RSMA can flexibly adjust the ratio of private and common signals according to channel conditions and user requirements, thus adapting to different network environments. Based on this, the transmitting signals are classified into private and common signals, and the transmitting signal at the AP can be represented as
x = a c p s s c + k = 1 K a k , t p p s s k , t p + k = 1 K a k , r p p s s k , r p ,
where p s is the total power of the signal transmitted by the AP. s c , s k , t p , and s k , r p are the common information and the private information of user T and user R in the k-th user cluster. The power allocation coefficients of the common and private information are denoted by a c and a k , m p , respectively, and satisfy a c + k = 1 K a k , t p + a k , r p = 1 .
This gives a receiving signal at user m in the k-th user cluster of
y k , m = g k , m H Θ k , m h k x + n k , m ,
where n k , m is the additive Gaussian white noise of user m, obeying n k , m C N 0 , σ n 2 .
According to the decoding order of RSMA, user m of the k-th user cluster will decode the common information first, and its signal-to-interference-plus-noise ratio (SINR) is
γ k , m c = a c p s g k , m H Θ k , m h k 2 p s k = 1 K a k , t p + a k , r p g k , m H Θ k , m h k 2 + σ n 2 .
After performing SIC on common information, the user detects their own private information. The remaining interference consists only of other users’ private information. The SINR for user m in the k-th user cluster is
γ k , m p = a k , p m p s g k , m H Θ k , m h k 2 p s i = 1 , i k K a i , m p + i = 1 K a i , m ¯ p g k , m H Θ k , m h k 2 + σ n 2 ,
in which m ¯ t , r . When m = t , m ¯ = r , and vice versa.

3. Ergodic Rate Analysis

3.1. Ergodic Rate under Discrete Phase Shifts

It is assumed that the STAR-RIS has known the instantaneous channel state information [20], which was obtainable through popular channel estimation techniques such as parallel factor decomposition. The parameters of STAR-RIS can be adjusted to maximize the SINR for achieving the highest ergodic rate. Considering that the transmitter has a single antenna, let the n-th element in g k , m be g k , m , n = g k , m , n e j ϕ k , m , n and the n-th element in h k be h k , n = h k , n e j φ k , n . The optimal phase θ k , m , n of the STAR-RIS is θ k , m , n * = ϕ k , m , n φ k , n . Due to the actual hardware limitation, θ k , m , n can only take a limited number of discrete values, represented by the set S = 0 , Δ θ , , 2 b 1 Δ θ , where Δ θ = 2 π / 2 b and b represents the number of quantization bits. The optimal discrete phase shift variable is then represented as [24]
θ ¯ k , m , n = Δ θ θ k , m , n * Δ θ + 1 2 .
From the above equation, it can be seen that the phase quantization error is θ k , m , n e = θ ¯ k , m , n θ k , m , n * , which obeys the uniform distribution of the interval on Δ θ / 2 , Δ θ / 2 , . According to the central limit theorem (CLT), when N 1 , X k , m = n = 1 N h ˜ k , n g ˜ k , n m e j θ k , m , n e / N obeys the complex Gaussian distribution, such that its real part U = X k , m and the imaginary part V = X k , m . The real and imaginary parts are independent of each other and obey U N μ , σ U 2 and V N 0 , σ V 2 , and there are the following relationships according to [25].
μ = φ 1 μ 0 ,
σ U 2 = 1 2 N 1 + φ 2 2 φ 1 2 μ 0 2 ,
σ V 2 = 1 2 N 1 φ 2 .
In Equations (11)–(13), φ p = E e j p θ k , m , n e is the characteristic function. Since the probability density function of θ k , m , n e is symmetric about the vertical axis, the calculation yields φ 1 = S a Δ θ / 2 and φ 2 = S a Δ θ , where S a x = sin x / x is the sampling function and μ 0 is the mean of n = 1 N h ˜ k , n g ˜ k , n m / N with value μ 0 = π / 4 .
Under conditions of N 1 and φ 1 > 0 , the distribution of X k , m = U 2 + V 2 can be approximated by the Nakagami-m distribution [25]. Denote Z k , m = X k , m 2 , then Z k , m obeys the gamma distribution and its cumulative distribution function is
F Z k , m z = γ λ , λ Ω z Γ λ ,
where Ω = μ 2 , λ = N φ 1 2 μ 0 2 / 1 + φ 2 2 φ 1 2 μ 0 2 / 2 and γ s , x is the lower incomplete gamma function [26]. The discrete phase shift matrix can be substituted into (8) and (9) and they can be simplified to obtain the SINR for the user to decode the common and private information as
γ k , m c = τ k , m a c Z k , m τ k , m 1 a c Z k , m + d k , m 2 + L k 2 α / 2 ,
γ k , m p = τ k , m a k , m p Z k , m τ k , m 1 a c a k , m p Z k , m + d k , m 2 + L k 2 α / 2 ,
where there is τ k , m = N 2 ρ 0 2 β k , m p s / σ n 2 d A P , k α . The ergodic rate of the common and private information of user m in the k-th user cluster can be calculated using the Shannon equation,
R ¯ k , m χ = E log 2 1 + γ k , m χ = 1 ln 2 0 1 F k , m χ x 1 + x d x ,
where χ c , p represents the rate of common and private messages sent by the corresponding user, respectively, and F k , m χ x denotes the cumulative distribution function of γ k , m χ . Substituting Equations (1)–(3) and (14) into (17), the ergodic rate of common and private messages of user m can be further derived as
R ¯ k , m χ = 1 ln 2 π / 2 π / 2 δ 1 δ 2 0 ϖ k , m χ 1 F Z k , m Τ k , m χ x Ψ y 1 + x f d k , m y f ϑ k , m z d x d y d z .
In the formula, there exist Τ k , m χ x = x / v k , m χ u k , m χ x and Ψ y = y 2 + L k 2 α / 2 . For the reflecting user, δ 1 = 0 and δ 2 = R 1 . The expressions for the common information rate are v k , r c = τ k , r a c , u k , r c = τ k , r 1 a c and ϖ k , r c = a c / 1 a c . The expressions for the private information rate are v k , r p = τ k , r a k , r p , u k , r p = τ k , r 1 a c a k , r p and ϖ k , r c = a k , r p / 1 a c a k , r p . For the transmitting user, δ 1 = R 1 and δ 2 = R 2 . The common information rate is calculated as v k , t c = τ k , t a c , u k , t c = τ k , t 1 a c and ϖ k , t c = a c / 1 a c , and the private information rate is calculated as v k , t p = τ k , t a k , t p , u k , t p = τ k , t 1 a c a k , t p and ϖ k , r c = a k , t p / 1 a c a k , t p . The calculation of the triple integral is
R ¯ k , m χ = 1 ln 2 δ 1 δ 2 0 ϖ k , m χ 1 F Z k , m Τ k , m χ x Ψ y 1 + x f d k , m y d x d y = 1 ln 2 δ 1 δ 2 0 ϖ k , m χ f d k , m y 1 + x d x d y Ι 1 1 ln 2 δ 1 δ 2 0 ϖ k , m χ F Z k , m Τ k , m χ x Ψ y 1 + x f d k , m y d x d y Ι 2 .
Calculating the results of terms Ι 1 and Ι 2 separately, it can be obtained that
Ι 1 = 1 δ 2 2 δ 1 2 ln 2 δ 1 δ 2 2 y 0 ϖ k , m χ 1 1 + x d x d y = log 2 1 + ϖ k , m χ ,
Ι 2 = ϖ k , m χ π 2 2 ln 2 δ 2 + δ 1 N p 2 Τ λ i = 1 N p j = 1 N p δ 2 δ 1 2 v j + δ 2 + δ 1 2 1 v i 2 1 v j 2 1 + ϖ k , m χ v i + 1 / 2 × γ λ , λ Ω Τ k , m χ ϖ k , m χ 1 + v i / 2 Ψ δ 2 δ 1 2 v j + δ 2 + δ 1 2 .
The specific calculation of Ι 2 is expanded. Firstly, the cumulative distribution function of Z k , m is substituted into the Ι 2 term of (19) with the expression
Ι 2 = 2 ln 2 δ 2 2 δ 1 2 Γ λ δ 1 δ 2 0 ϖ k , m χ y 1 + x γ λ , λ Ω Τ k , m χ x Ψ y d x d y .
Subsequently, a change of variables is performed for the inner integral and then the Chebyshev–Gauss quadrature formula is used for approximation, obtaining
0 ϖ k , m χ y 1 + x γ λ , λ Ω Τ k , m χ x Ψ y d x = ϖ k , m χ 2 1 1 y 1 + ϖ k , m χ x + 1 / 2 γ λ , λ Ω Τ k , m χ ϖ k , m χ x + 1 / 2 Ψ y d x = ϖ k , m χ π 2 N p i = 1 N p y 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 γ λ , λ Ω Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y ,
where v i = cos 2 i 1 2 N p π .
Thus, Ι 2 can be converted into a single integral with the expression
Ι 2 = ϖ k , m χ π ln 2 δ 2 2 δ 1 2 Γ λ N p i = 1 N p δ 1 δ 2 y 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 γ λ , λ Ω Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y Ι 3 d y .
Processing Ι 3 using a similar method of Chebyshev–Gauss quadrature, the expression can be rewritten as
Ι 3 = δ 1 δ 2 y 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 γ λ , λ Ω Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y d y = δ 2 δ 1 2 1 1 δ 2 δ 1 t / 2 + δ 2 + δ 1 / 2 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 × γ λ , λ Ω Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 t / 2 + δ 2 + δ 1 / 2 d t = δ 2 δ 1 π 2 N p j = 1 N p δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 1 v i 2 1 v j 2 1 + ϖ k , m χ v i + 1 / 2 × γ λ , λ Ω Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 ,
where v j = cos 2 j 1 2 N p π .
The conclusion of Equation (21) is obtained by collating Equations (24) and (25). Finally, substituting Equations (20) and (21) into Equation (19) yields an approximate closed expression for R ¯ k , m χ .

3.2. Ergodic Rate under Continuous Phase Shifts

When STAR-RIS has continuous phase shifts, the value of the optimal phase θ k , m , n is θ k , m , n * = ϕ k , m , n φ k , n . Since the channels of AP to STAR-RIS and STAR-RIS to users are independent of each other and both h ˜ k and g ˜ k , n m obey a complex Gaussian distribution with mean 0 and variance 1, h ˜ k , n and g ˜ k , n m obey a Rayleigh distribution. Let X k = n = 1 N h ˜ k , n g ˜ k , n m ; according to the CLT, it can be obtained that X k N N μ 0 , N σ 0 2 , where μ 0 = π 2 × π 2 = π 4 and σ 0 2 = 1 π 2 / 16 [27]. Let Y k , m = n = 1 N h ˜ k , n g ˜ k , n m 2 / N σ 0 2 ; then, Y k , m obeys the non-central chi square distribution of Y k , m χ 1 2 λ in which λ = N μ 0 2 / σ 0 2 . The cumulative distribution function is
F Y k , m = 1 Q 1 2 λ , x = 1 1 2 e r f c x + λ 2 + e r f c x λ 2 ,
As a result, the expressions (8) and (9) for decoding common and private information for user m in the k-th user cluster can be re-expressed as
γ k , m c = τ k , m a c Y k , m τ k , m 1 a c Y k , m + d k , m 2 + L k 2 α / 2 ,
γ k , m p = τ k , m a k , m p Y k , m τ k , m 1 a c a k , m p Y k , m + d k , m 2 + L k 2 α / 2 ,
where τ k , m = N σ 0 2 ρ 0 2 β k , m p s / σ n 2 d A P , k α . The ergodic rate for calculating the common and private information of user m in the k-th user cluster is given by
R ¯ k , m χ = E log 2 1 + γ k , m χ = 1 ln 2 0 1 F k , m χ x 1 + x d x .
Similar to the method of analyzing the ergodic rate of the system under discrete phase shifts, it can be further obtained that
R ¯ k , m χ = 1 ln 2 π / 2 π / 2 δ 1 δ 2 0 ϖ k , m χ 1 F Y k , m Τ k , m χ x Ψ y 1 + x f d k , m y f ϑ k , m z d x d y d z ,
where Τ k , m χ x = x / v k , m χ u k , m χ x and Ψ y = y 2 + L k 2 α / 2 . For the reflecting user, δ 1 = 0 and δ 2 = R 1 . The expressions for the common information rate are v k , r c = τ k , r a c , u k , r c = τ k , r 1 a c and ϖ k , r c = a c / 1 a c . The expressions for the private information rate are v k , r p = τ k , r a k , r p , u k , r p = τ k , r 1 a c a k , r p and ϖ k , r c = a k , r p / 1 a c a k , r p . For the transmitting user, δ 1 = R 1 and δ 2 = R 2 . The common information rate is calculated as v k , t c = τ k , t a c , u k , t c = τ k , t 1 a c and ϖ k , t c = a c / 1 a c , and the private information rate is calculated as v k , t p = τ k , t a k , t p , u k , t p = τ k , t 1 a c a k , t p , and ϖ k , r c = a k , t p / 1 a c a k , t p . The expression of the R ¯ k , m χ is
R ¯ k , m χ = 1 ln 2 δ 1 δ 2 0 ϖ k , m χ 1 F Y k , m Τ k , m χ x Ψ y 1 + x f d k , m y d x d y = 1 ln 2 δ 1 δ 2 0 ϖ k , m χ 1 F Y k , m Τ k , m χ x Ψ y 1 + x f d k , m y d x d y .
Substitute Equation (26) into (31) to obtain
R ¯ k , m χ = 1 ln 2 δ 2 2 δ 1 2 δ 1 δ 2 0 ϖ k , m χ y e r f c Τ k , m χ x Ψ y + λ 2 + e r f c Τ k , m χ x Ψ y λ 2 1 + x d x d y .
The inner integrals are permuted and processed using the Chebyshev–Gauss quadrature formula, which provides
ϖ k , m χ 2 1 1 y e r f c Τ k , m χ ϖ k , m χ x + 1 / 2 Ψ y + λ 2 + e r f c Τ k , m χ ϖ k , m χ x + 1 / 2 Ψ y λ 2 1 + ϖ k , m χ x + 1 / 2 d x d y = ϖ k , m χ π 2 N p i = 1 N p y 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y + λ 2 + e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y λ 2 .
Then, R ¯ k , m χ can be re-expressed as
R ¯ k , m χ = ϖ k , m χ π 2 ln 2 δ 2 2 δ 1 2 N p × i = 1 N p δ 1 δ 2 y 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y + λ 2 + e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y λ 2 d y Κ 1 .
After integrating Κ 1 by permutation, the term Κ 1 in the above Equation (34) can be re-written as
Κ 1 = δ 1 δ 2 y 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y + λ 2 + e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ y λ 2 d y = δ 2 δ 1 2 1 1 δ 2 δ 1 t / 2 + δ 2 + δ 1 / 2 1 v i 2 1 + ϖ k , m χ v i + 1 / 2 × e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 t / 2 + δ 2 + δ 1 / 2 + λ 2 + e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 t / 2 + δ 2 + δ 1 / 2 λ 2 d y .
Next, the Chebyshev–Gauss quadrature formula can be used to obtain the approximate closed-form solution as
Κ 1 = δ 2 δ 1 π 2 N p j = 1 N p δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 1 v i 2 1 v j 2 1 + ϖ k , m χ v i + 1 / 2 × e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 + λ 2 + e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 λ 2 ,
where v i = cos 2 i 1 2 N p π and v j = cos 2 j 1 2 N p π .
Substituting the expression for Κ 1 into Equation (34) leads to the approximate closed expression of R ¯ k , m χ as
R ¯ k , m χ = ϖ k , m χ π 2 4 ln 2 δ 2 + δ 1 N p 2 i = 1 N p j = 1 N p δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 1 v i 2 1 v j 2 1 + ϖ k , m χ v i + 1 / 2 × e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 + λ 2 + e r f c Τ k , m χ ϖ k , m χ v i + 1 / 2 Ψ δ 2 δ 1 v j / 2 + δ 2 + δ 1 / 2 λ 2 .
Thus, the sum ergodic rate under discrete and continuous phase shift systems is expressed as
R ¯ t o t = min R ¯ k , m c + k = 1 K R ¯ k , t p + R ¯ k , r p .

4. Simulation and Analysis

This paper analyzes the sum ergodic rate of an STAR-RIS-assisted RSMA downlink communication model under both discrete and continuous phase shift conditions. It also compares the impact of different factors on system performance and verifies the correctness of the ergodic rate derivation through Monte Carlo simulations. In this model, users are clustered at the transmitter, and their information is processed based on RSMA. Each user cluster is served by an STAR-RIS, and the mutual influence between different user clusters is not considered.
For simulation, six users and three user clusters are considered. It is assumed that the position of the AP with respect to each STAR-RIS is known and the distances between the AP and the STAR-RIS are d 1 = 30 m, d 2 = 40 m, and d 3 = 50 m. The same distribution of user positions is considered in each user cluster, which is R 1 = 10 m and R 2 = 20 m. The order of the Gaussian–Chebyshev polynomials N p is 50, and the number of Monte Carlo simulation times is set to 10 6 . The defaults for the rest of the simulation parameters are shown in Table 1.
Figure 2 shows the curves of system’s sum ergodic rate under different numbers of quantization bits and continuous phase shifts. In the simulation, the power ratio allocated to the common information part is set by default to a c = 0.6 , and the remaining transmit power is equally divided among the private information parts of all users. From the figure, it can be observed that the system ergodic rate at b = 1 is much lower than that at b = 2 and b = 3 . Additionally, it can be seen that when b = 3 , the system ergodic rate is already very close to the ergodic rate under continuous phase shifts. Therefore, we can conclude that when designing quantized phase shifts, the number of quantization bits only needs to reach 3 to achieve a rate performance close to that of continuous phase shifts.
Figure 3 and Figure 4 compare the ergodic rate under different numbers of STAR-RIS elements in the cases of continuous phase shifts and discrete phase shifts, respectively. As shown in Figure 3 and Figure 4, the theoretical curves of the system ergodic rate match well with the simulation curves. In both cases, as the number of elements increases, the total ergodic rate of the system increases accordingly. Thus, the more elements the STAR-RIS has, the greater the system gain it can bring. Additionally, since the impact of noise on the system can be ignored under high SNR, the system rate is a constant and is independent of the number of STAR-RIS elements under high SNR. Hence, the gain effect of STAR-RIS is more significant under the condition of low transmit power at the transmitter.
In Figure 5 and Figure 6, under the conditions of continuous phase shifts and discrete phase shifts of STAR-RIS, the value of the path loss exponent α is changed, with α taking values of 2.0, 2.2, and 2.4, respectively. It can be seen that the ergodic rate for the system increases as the pass loss exponent decreases. This occurs because, with the communication distance and reference path loss remaining unchanged, a smaller α results in higher channel gain of the Rayleigh channel, which in turn leads to better channel quality. Thus, the user’s receiving end in this channel has a higher SINR. From both Figure 5 and Figure 6, it can be seen that regardless of whether the phase shifts are discrete or continuous, minimizing the path loss exponent as much as possible will bring better performance gains to the system. Additionally, under low transmit power conditions, reducing the path loss exponent results in a more significant increase in the ergodic rate of the STAR-RIS continuous phase shift system.
Figure 7 compares the total system ergodic rate for both STAR-RIS ES and MS modes with continuous and discrete phase shifts. To ensure fairness, the MS transmission protocol has N / 2 elements for transmission and reflection, respectively. From the figure, it can be seen that the total ergodic rate in the ES mode is higher than that in the MS mode for both discrete phase shift and continuous phase shift cases. This is due to the fact that only subsets of the elements are selected for transmission and reflection, respectively, in the MS mode, and the amplitude coefficients for transmission and reflection are limited to binary values. As a result, MS is usually unable to achieve the same full-dimensional transmission and reflection beamforming gains as ES. It is evident that the ergodic sum rate in the MS mode with quantization bit b = 1 is lower than the ergodic sum rate with continuous phase shifts, which is consistent with the findings in Figure 2.

5. Conclusions

This paper mainly focuses on the analysis and comparison of the ergodic rate and performance of an STAR-RIS-assisted RSMA system under discrete and continuous phase shifts of STAR-RIS. Firstly, users are clustered at the transmitter, ensuring no interference between different user clusters. Additionally, it is considered that users are uniformly distributed within the range centered on the STAR-RIS. Secondly, the cumulative distribution function of the users’ SINR is derived based on the central limit theorem, and the approximate closed-form expression of the ergodic rate is obtained using Gaussian–Chebyshev integration. Monte Carlo simulations are conducted to show that the theoretical curves match the simulation curves, further validating the rationality and correctness of the theoretical derivation. During the simulations, by comparing the system performance under different phase shift quantization bits, it is concluded that to achieve performance close to the ideal continuous phase shift, the number of quantization bits for the phase shift should not be less than 3. Apart from this, it is found that increasing the number of elements and reducing the path loss exponent can improve the system performance to a certain extent.

Author Contributions

Conceptualization, T.L.; methodology, Y.Z.; software, Y.Z.; writing—original draft preparation, Y.Z.; writing—review and editing, T.L.; All authors have read and agreed to the published version of the manuscript.

Funding

This research and APC were funded by the Natural Science Foundation of China (No. 61601220).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

Explanations of meaning of parameters used in the paper.
ParametersExplanations
A × B complex   matrices   of   dimension   A × B
C N a , b complex Gaussian distribution with mean a and variance b
absolute value or modulus of a complex number or number of elements of a set
d i a g x diagonal matrix with diagonal elements as vectors x
E · average value of a variable
S a · sampling function
γ s , x lower incomplete gamma function
Γ · gamma function
Q M t , s Markum Q Function
e r f c · error complementarity function

References

  1. Gunasinghe, D.; Amarasuriya, G. Performance analysis of STAR-RIS for wireless communication. In Proceedings of the ICC 2022-IEEE International Conference on Communications, Seoul, Republic of Korea, 16–20 May 2022; pp. 3275–3280. [Google Scholar]
  2. 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]
  3. Wu, Q.; Zhang, S.; Zheng, B.; You, C.; Zhang, R. Intelligent Reflecting Surface-Aided Wireless Communications: A Tutorial. IEEE Trans. Commun. 2021, 69, 3313–3351. [Google Scholar] [CrossRef]
  4. Bjornson, E.; Ozdogan, O.; Larsson, E.G. Intelligent Reflecting Surface Versus Decode-and-Forward: How Large Surfaces are Needed to Beat Relaying? IEEE Wirel. Commun. Lett. 2019, 9, 244–248. [Google Scholar] [CrossRef]
  5. Tang, W.; Chen, M.Z.; Chen, X.; Dai, J.Y.; Han, Y.; Di Renzo, M.; Zeng, Y.; Jin, S.; Cheng, Q.; Cui, T.J. Wireless Communications With Reconfigurable Intelligent Surface: Path Loss Modeling and Experimental Measurement. IEEE Trans. Wirel. Commun. 2020, 20, 421–439. [Google Scholar] [CrossRef]
  6. Xu, J.; Liu, Y.; Mu, X.; Dobre, O.A. STAR-RISs: Simultaneous transmitting and reflecting reconfigurable intelligent surfaces. IEEE Commun. Lett. 2021, 25, 3134–3138. [Google Scholar] [CrossRef]
  7. Xiang, Z.; Tong, X.; Cai, Y. Secure transmission for NOMA systems with imperfect SIC. China Commun. 2020, 17, 67–78. [Google Scholar] [CrossRef]
  8. Lin, Z.; Lin, M.; Wang, J.B.; Cola, T.d.; 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]
  9. Lin, Z.; Lin, M.; Cola, T.d.; Wang, J.B.; Zhu, W.P.; Cheng, J. Supporting iot with rate-splitting multiple access in satellite and aerial-integrated networks. IEEE Internet Things J. 2021, 8, 11123–11134. [Google Scholar] [CrossRef]
  10. Mao, Y.; Dizdar, O.; Clerckx, B.; Schober, R.; Popovski, P.; Poor, H.V. Rate-splitting multiple access: Fundamentals, survey, and future research trends. IEEE Commun. Surv. Tutor. 2022, 24, 2073–2126. [Google Scholar]
  11. Mao, Y.; Clerckx, B.; Li, V.O. Rate-splitting multiple access for downlink communication systems: Bridging, generalizing, and outperforming SDMA and NOMA. EURASIP J. Wirel. Commun. Netw. 2018, 2018, 133. [Google Scholar] [CrossRef] [PubMed]
  12. Ho, C.M.; Nguyen, T.V.; Lee, C.; Nguyen, T.M.T.; Cho, S. Cooperative user relaying with rsma for 6g networks: Overview, research challenges and future trends. In Proceedings of the 2023 International Conference on Information Networking (ICOIN), Bangkok, Thailand, 11–14 January 2023; pp. 584–587. [Google Scholar]
  13. Rivera, M.; Jaafar, W.; Yanikomeroglu, H. Optimization of RIS-assisted RSMA-enabled tethered-UAV communications. In Proceedings of the 2023 IEEE Globecom Workshops (GC Wkshps), Kuala Lumpur, Malaysia, 4–8 December 2023; pp. 227–232. [Google Scholar]
  14. Narottama, B.; Aïssa, S. Quantum machine learning for performance optimization of RIS-assisted communications: Framework design and application to energy efficiency maximization of systems with RSMA. IEEE Trans. Wirel. Commun. 2024. early acces. [Google Scholar] [CrossRef]
  15. Singh, S.K.; Agrawal, K.; Singh, K.; Clerckx, B.; Li, C.-P. RSMA enhanced RIS-FD-UAV-aided short packet communications under imperfect SIC. In Proceedings of the 2022 IEEE Globecom Workshops (GC Wkshps), Rio de Janeiro, Brazil, 4–8 December 2022; pp. 1549–1554. [Google Scholar]
  16. Jolly, A.; Singh, K.; Biswas, S. RSMA for IRS aided 6G communication systems: Joint active and passive beamforming design. In Proceedings of the 2021 IEEE International Conference on Advanced Networks and Telecommunications Systems (ANTS), Hyderabad, India, 13–16 December 2021; pp. 7–12. [Google Scholar]
  17. Karim, F.; Singh, S.K.; Singh, K.; Prakriya, S.; Li, C.-P. Performance analysis for RSMA-empowered STAR-RIS-aided downlink communications. In Proceedings of the 2023 IEEE 34th Annual International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC), Toronto, ON, Canada, 5–8 September 2023; pp. 1–6. [Google Scholar]
  18. Dhok, S.; Sharma, P.K. Rate-splitting multiple access with STAR RIS over spatially-correlated channels. IEEE Trans. Commun. 2022, 70, 6410–6424. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Yang, L.; Li, X.; Guo, K.; Liu, H. Covert communications for star-ris-assisted industrial networks with a full duplex receiver and RSMA. IEEE Internet Things J. 2024, 11, 22483–22493. [Google Scholar] [CrossRef]
  20. Mohamed, Z.; Albaden, K.; Aïssa, S. Simultaneously transmitting and reflecting reconfigurable intelligent surface aided RSMA communications: Outage probability analysis. In Proceedings of the 2023 IEEE 34th Annual International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC), Toronto, ON, Canada, 5–8 September 2023; pp. 1–6. [Google Scholar]
  21. 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]
  22. Wu, C.; You, C.; Liu, Y.; Gu, X.; Cai, Y. Channel estimation for STAR-RIS-aided wireless communication. IEEE Commun. Lett. 2021, 26, 652–656. [Google Scholar] [CrossRef]
  23. Sobhi-Givi, S.; Shayesteh, M.G.; Kalbkhani, H. Energy-efficient power allocation and user selection for mmWave-NOMA transmission in M2M communications underlaying cellular heterogeneous networks. IEEE Trans. Veh. Technol. 2020, 69, 9866–9881. [Google Scholar] [CrossRef]
  24. Yang, S.; Zhang, J.; Xia, W.; Ren, Y.; Yin, H.; Zhu, H. On the discrete phase shifts design for distributed RIS-aided downlink MIMO-NOMA systems. In Proceedings of the 2022 IEEE Wireless Communications and Networking Conference (WCNC), Austin, TX, USA, 10–13 April 2022; pp. 363–368. [Google Scholar]
  25. Badiu, M.-A.; Coon, J.P. Communication through a large reflecting surface with phase errors. IEEE Wirel. Commun. Lett. 2019, 9, 184–188. [Google Scholar] [CrossRef]
  26. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  27. Anjana, B.; Jadhav, H.; Kumaravelu, V.B.; Soria, F.R.C.; Sayeed, M.S.; Murugadass, A. Smart reconfigurable intelligent surface with discrete phase shifter for next generation networks. In Proceedings of the 2022 International Conference on Wireless Communications Signal Processing and Networking (WiSPNET), Chennai, India, 24–26 March 2022; pp. 178–182. [Google Scholar]
Figure 1. System model of STAR-RIS.
Figure 1. System model of STAR-RIS.
Sensors 24 05480 g001
Figure 2. The system ergodic rate under different numbers of phase shift quantization bits.
Figure 2. The system ergodic rate under different numbers of phase shift quantization bits.
Sensors 24 05480 g002
Figure 3. The system ergodic rate of different element numbers under continuous phase shifts.
Figure 3. The system ergodic rate of different element numbers under continuous phase shifts.
Sensors 24 05480 g003
Figure 4. The system ergodic rate of different element numbers under discrete phase shifts.
Figure 4. The system ergodic rate of different element numbers under discrete phase shifts.
Sensors 24 05480 g004
Figure 5. The system ergodic rate of different path loss exponents under continuous phase shifts.
Figure 5. The system ergodic rate of different path loss exponents under continuous phase shifts.
Sensors 24 05480 g005
Figure 6. The system ergodic rate of different path loss exponents under discrete phase shifts.
Figure 6. The system ergodic rate of different path loss exponents under discrete phase shifts.
Sensors 24 05480 g006
Figure 7. The system ergodic rate of different transmission modes.
Figure 7. The system ergodic rate of different transmission modes.
Sensors 24 05480 g007
Table 1. Simulation parameters set based on discrete and continuous phase shifts of the STAR-RIS-assisted RSMA multi-user system.
Table 1. Simulation parameters set based on discrete and continuous phase shifts of the STAR-RIS-assisted RSMA multi-user system.
ParametersDefault ValueParametersDefault Value
number of STAR-RIS
elements
N = 64 noise power σ n = 80 dBm
path loss exponent α 0 = 2.2 path loss at a reference
unit distance
ρ 0 = 30 dB
number of
quantization bits
b = 1 amplitude coefficient
of STAR-RIS
β r = β t = 0.5
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

Liu, T.; Zhou, Y. Ergodic Rate Analysis of Simultaneous Transmitting and Reflecting Reconfigurable Intelligent Surface-Assisted Rate-Splitting Multiple Access Systems Based on Discrete Phase Shifts. Sensors 2024, 24, 5480. https://doi.org/10.3390/s24175480

AMA Style

Liu T, Zhou Y. Ergodic Rate Analysis of Simultaneous Transmitting and Reflecting Reconfigurable Intelligent Surface-Assisted Rate-Splitting Multiple Access Systems Based on Discrete Phase Shifts. Sensors. 2024; 24(17):5480. https://doi.org/10.3390/s24175480

Chicago/Turabian Style

Liu, Tao, and Yue Zhou. 2024. "Ergodic Rate Analysis of Simultaneous Transmitting and Reflecting Reconfigurable Intelligent Surface-Assisted Rate-Splitting Multiple Access Systems Based on Discrete Phase Shifts" Sensors 24, no. 17: 5480. https://doi.org/10.3390/s24175480

APA Style

Liu, T., & Zhou, Y. (2024). Ergodic Rate Analysis of Simultaneous Transmitting and Reflecting Reconfigurable Intelligent Surface-Assisted Rate-Splitting Multiple Access Systems Based on Discrete Phase Shifts. Sensors, 24(17), 5480. https://doi.org/10.3390/s24175480

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