Low‑complexity cooperative active and passive beamforming multi‑RIS‑assisted communication networks

.

reflecting elements that can be adjusted to manipulate the phase shifts and amplitudes of electromagnetic waves.By collaboratively modifying the reflecting coefficients, RIS enhances the performance of non-line-of-sight links and overall system performance.RIS is cost-effective and power-efficient compared to conventional transmitters or relays.Numerous studies have proposed models and analyzed performance metrics to improve wireless network coverage areas by deploying RIS and its reflective elements.In the beginning, the primary focus of the research was on exploring the potential improvements in wireless communications achieved by utilizing a single RIS.These studies investigated various aspects such as energy efficiency [4][5][6], weighted sum rate [7][8][9], network coverage [10][11][12].The key outcomes of these studies underscore the substantial RIS deploying impact and leveraging its reflective elements to enhance network coverage.
Tons of papers have studied and designed RIS-assisted wireless networks because of the great impact of utilizing RIS capabilities in enhancing faded channels and blocking problems.The research begins by investigating joint active and passive beamforming (JAPBF) between BS and single RIS and then upgraded to utilizing cooperation of passive beamforming between multi-RISs in different topologies, i.e., parallel, cascaded, and hybrid ways [13][14][15][16][17][18][19].

Motivation and contribution
Unfortunately, assisting users with all deployed RISs is not an efficient way, because the signals reflected by the RISs are severely attenuated as shown in Fig. 1.Because of the mobility of users, in many cases, some RISs may be far from than user.By this motivation, it is critical to merge a new role as RIS selection (RIS-S) in multi-RIS assisting.RIS-S strategy is investigated through two schemes: a single selected RIS and a set of RISs to serve the user.RIS-S based on a single selected surface, authors focused on selecting the main and best-reflected path, [20][21][22][23][24]. Unlike RIS-S based on a set of Fig. 1 The conventional multi-hop RIS selected surfaces, authors focused on choosing the best-reflected paths [25,26].Figures 1 and 2 show the conventional multi-hop RIS assisted without any selection and multi-hop RIS-assisted network after investigating the RIS-S strategy, respectively.
The authors in [20] proposed two multi-RIS-aided schemes called exhaustive RISaided (ERA) and opportunistic RIS aided (ORA).It highlighted the advantages and drawbacks in the cases of exploiting all surfaces in assisting as ERA scheme or one scheduled RIS is exploited in assisting as ORA scheme.ORA scheme achieved low complexity of backhaul and fronthaul load.By investigating the parallel topology of multi-RIS assisted, [21] exploited jointly the RIS-S to active and passive beamforming (APBF) optimization by selecting the optimal combination between RISs and users.By assuming a RIS-S matrix includes all possible combinations set and by alternating manner RIS-S matrix is updated by the optimal combinations based on the APBF.Unlike in [21], cooperative optimization between the APBF and a RIS-S strategy is presented in [22], where the graph neural network (GNN) investigated to help in reaching the optimal RISs to construct the optimal topology of the wireless network.In a different way, authors in [23] optimized the outage probability and average symbol error probability by selecting the best RIS.The selection criteria of RIS-S were to select the RIS which has the best signal-to-noise ratio (SNR).Then, the authors extended this work to UAV assisted in case of multi-RISs distributed, by studying the optimal RIS placement to enhance the outage probability and the received signal.Hence, the selection strategy is based on the received signal power at UAV, not the destination, [24].In [25], authors differ from the previous work by assuming a set of RISs is selected as the best surfaces for assisting users, instead of exploiting all surfaces.Motivated by minimizing the overall consumed power in the network, a dynamic controlling of involving surface in assisting i.e., RIS has binary status on-off, the energy efficiency of the network is maximized in [26].
Unlike [21] and [22], a RIS-S strategy proposed to be applied as a pre-stage before APBF optimization.This approach reduces the complexity of the optimization problem Fig. 2 The proposed RIS-S strategy in multi-hop RIS by selecting the proper RISs for each user.Moreover, multi-hop RIS investigated to enhance guiding the signal to the endpoint by mitigating the scattered path losses, because the signals reflected by the far RISs are severely attenuated.
Aforementioned, [21][22][23][24][25], authors utilized different optimization techniques to relax and decouple the problem into solved subproblems.In this study, the objective function involves a ratio of two functions.However, there are several other types of optimization problems that are commonly encountered to solve it such as convex optimization, stochastic optimization, and dynamic programming.Due to the non-convexity of the problem, fractional programming investigated to reformulate the problem into a convex and decouple the problem into subproblems that could be analyzed.Then quadratic transform applied to reshape the objective function in a way that is solved.
In summary, our contributions are: We propose a RIS-S algorithm that specifies the best channels between the RISs and users based on a significant SNR.This selection role is a pre-stage before the APBF optimization, unlike previous works which selected the proper RISs.
1.In contrast to the previous work, all users are associated with a single RIS.We maximize the SINR, which by default leads to maximization of the sum rate by comprising the channel model and active beamforming (ABF) of BS, but it seems to be nonconvex and multivariate.For the purpose of resolving this problem, by employing closed-form fractional programming methods, we are able to effectively decouple the original problem into three independent sub-problems.In a similar way [18], by introducing new variables, we are able to split the sum-of-logarithms-of-ratio problem into two terms, each of which may then be simplified independently using the quadratic transform.2. Due to performing the RIS-S as a pre-stage before the APBF optimization, the complexity of the APBF optimization design is reduced.
The following sections will organize the remainder of the paper: In Sect.2, both the model of the system and the RIS selection scheme are introduced and summarized in the proposed algorithm.In Sect.3, we present the suggested problem formulation of the JAPBF technique based on RIS-S for multi-hop RIS-assisted multi-user wireless communication networks.The results of the simulation are presented in Sect. 4. In Sect.5, we summarize our study.

System and channel model
In our proposed mmWave wireless communication network with multiple RISs, as shown in Fig. 2, establishing direct connections between the BS and users faces significant obstacles.Our approach is like a multi-arm network, where N distributed RISs assist in transmitting data from the BS, equipped with M antennas, to K users, each with a single antenna.As mentioned earlier, we anchor the first RIS hop, denoted as RIS l,1 for any l-th virtual route where (l = 1, 2, ..., L), very close to the BS to ensure a line-of-sight (LoS) channel directly above the BS.Also, each virtual route represents Q of RISs chain which equals N /L .The remaining RISs, denoted as RIS l,n (i.e., RIS n = 2, 3, ..., Q ), are distributed in a sequential manner to create multi-hop routes in a multi-arm configuration.Each RIS is composed of N e reflecting elements, with each element connected to a controller responsible for controlling the phase shifts.
In order to tackle the issue of multi-hop RIS routes and the presence of arbitrary reflection paths.We let G l ∈ C N e ×M , H l,n ∈ C N e ×N e and g l,n,k ∈ C 1×N e indicate the channels from the BS to RIS l,1 , and RIS l,n to RIS l,n+1 and from RIS L,Q to user k, respectively.The interpretation of the reflection coefficient matrix of RIS l,n can be described as l,n ∈ C N e ×N e , for more specificity, it can be further represented as l,n = diag(φ l,n ) .φ l,n = [φ l,n,1 , φ l,n,2 , …, φ l,n,N e ] T is the reflection vector of the (l, n)-th RIS, where φ l,n,n e = β l,n,n e e jθ l,n,ne , n e = 1, 2, …, N e .In addition, θ l,n,n e ∈ [0, 2π ] and β l,n,n e ∈ [0, 1] rep- resents the alteration in phase and amplitude caused by the n e -th reflecting element of RIS l,n on the incoming signals.We assume that the elements of the RIS are fully reflective, which implies that β = 1.
where β G l represents the free space path loss; F G l is the Rician factor; the matrices G l,LoS and G l,NLoS are the line-of-sight (LoS) and non-LoS (NLoS) components of the channel G l , respectively.By the same, we can express the channel g l,Q,k from RIS L,Q to user k as where β g l,Q,k represents the free space path loss; F g l,Q,k is the Rician factor; the matri- ces g l,Q,k LoS and g l,Q,k NLoS are the LoS and NLoS components of the channel g l,Q,k , respectively.
Similarly, the remaining channels can be established using the same procedure.Once they are properly installed, high-quality LoS channels can be formed between them.Consequently, it is possible to assume that the channel state information (CSI) for the link between the BS and the first RIS, as well as the inter-RIS links, is already known.In RIS-assisted systems, there are two methods to obtain CSI for the RIS-user channel, depending on whether the reflective elements of the RIS have sensing capabilities.As RISs are usually positioned in areas conducive to establishing the LoS connections with users, such as the facades of tall buildings, the CSI can be acquired using location-based channel estimation algorithms [27][28][29], even if the RIS is only purely passive.

RIS selection strategy algorithm
As aforementioned, RIS-S role adds ability to choose the optimal RISs for association with the UE, considering the channel conditions.Several studies have proposed channel estimation methods in the presence of RIS by analyzing the cascaded channel from the BS to RIS and from RIS to user equipment (UE), [29][30][31].Efficient algorithms were proposed to estimate the RIS-related channels.For instance, in [30], a compressive sensing (1) (2) and deep learning approach was proposed by deploying a few active antenna elements on the RIS.We assume that perfect CSI is available.Moreover, we assume that all users have the possibility to associate with all RISs simultaneously [18].Due to the random scattering of signals, we assume the far surfaces can receive the signals through multi-hop RIS to mitigate the losses.Hence, the distributed RISs can simultaneously reflect the beam toward other surfaces and UEs by playing two roles routing and accessing.So, it could be divided into two groups: (1) Access RIS (ARIS): It is the RIS in the RIS chain.It reflects the beam to user equipment UE. (2) Routing RIS (RRIS): Its role is to focus the electromagnetic beams toward the next RIS hop.So, a chain of dynamically aligned RISs plays the dedicated role of routing.However, it minimizes interference efficiently.
As known to the reader, the interaction channel between two RISs is affected by some parameters such as the orientation angle of the surface, the incident angle of the beam, and the surface area.Also, the channel between the RIS and the user is affected by the arrival and departure of angles besides RIS and UE locations.Without loss of generality, all channels are assumed as Rician modeled, and we took into consideration the arrival and departure angles to include the orientation angles of the surface implicitly, while the routing ability between the surfaces is introduced and studied in [11,19] to construct a chain of RISs due to guiding signals into poor coverage areas.On the basis, to simply represent the selection matrix between RISs and user k is s k = s k,1 , . . ., s k,N ∈ Z 1×N consisting of binary variables, and s k,n is defined as For clarity, the selection decision role depends on the SNR of the cascaded channel between the BS to RIS and the RIS to UE, respectively.By this way, we ensure that the selected RIS satisfies high received power and at the same time is closer to UE.However, the channel of the nearest surface for user k may be a good path, this surface is far away from BS, and the received power will be attenuated.
Thus, specifying the γ T value is more critical to decide.We have two options to deal with this point.The threshold value could be a fixed or dynamic changed value.Option 1, γ T is a fixed value, it is not fair decision to select the more effective surfaces due to user mobility and channel condition variations.This means, if the fixed value is high and the surrounding channels are noisy, the selected surfaces are very low and may not be the best selection.So, practically the dynamic changed value is more convenient.On the basis, this study assumed that γ T is adaptively specified, by estimating the mean value between the highest and lowest values of possible γ n,k for user k .On the basis, the γ T and the number of selected surfaces are changed corresponding to each user k.
Initially, we assume that all surfaces have a possibility of reflection to user k , so the selection matrix s k could be one's entries.SNR is estimated for each RIS to user k .Then, due to updating s k , we define γ T as a SNR threshold.s k,n could be updated by "1" for all (3) 1, user k is assisted by RIS l,n .0, otherwise.
values higher than γ T .In this way, we ensure that all the effective surfaces are selected, and the un-useful surfaces are suppressed.The above steps are summarized in the RIS-S algorithm below.
Algorithm 1 RIS selection

Problem formulation
In this section, we start by denoting the downlink signal from the BS to user k as x k with E |x k | 2 = 1 , k = 1, 2, . . .K .Subsequently, we can express the K transmitted signals as where w k ∈ C M×1 indicates the active beamforming vector for user k .By modeling the channel, which incorporates all potential paths and routes originating from BS to user k can be expressed as matrix h k where the vector A ∈ C (1×Q) represents all transmitted beams from BS to RIS l,1 and the vector B ∈ C (1×Q) represents the paths and routes between RIS l,n−1 to RIS l,n and RIS l,Q to user k Then, the updated channel model, after selecting the effective paths and routes and excluding weak paths by investigating the selection matrix as (4) based on, the k-th received signal corresponding to the updated channel form h H s,k as where n k ∼ C N 0,σ 2 0 represents the additive white Gaussian noise (AWGN) at the k-th user.Additionally, the k-th user considers the interference caused by other users' signals.Therefore, we can express the Signal-to-Interference-plus-Noise Ratio (SINR) at user k as: with K k=1 �w k � 2 ≤ P T represents the power constraint, where P T is the maximum trans- mitted power at the BS.Maximizing the SINR relies on the BS leveraging spatial diversity of power to direct signals for user k , i.e., �w k � 2 ≤ P T .Consequently, to maximize the SINR there is γ max can be achieved, and as a result, we have 0 represent the transmit beamforming matrix.We aim to maximize the total user rate during downlink transmission by jointly designing the transmit beamforming matrix W at the BS and the phase shift matrices at φ l,n each RIS l,n , while ensuring compliance with the transmit power at the BS and the phase shift constraints of the RISs.In this way, the optimization problem can be represented as (P1) is formulated to address the ABF at BS and the PBF for RIS l,n , W , and φ 1,1 , . . ., φ L,Q , respectively.The objective is to minimize inter-user interference on the receiving end and maximize the sum rate for multiple users.However, due to the coupling of the objective variables in the joint APBF optimization problem, finding the optimal solution becomes challenging.In the following section, we will adopt a similar approach as in [18] to alleviate this challenge and obtain a sub-optimal solution for (P1).

Optimal active and passive beamforming
As previously mentioned, our proposed JAPBF aims to maximize the sum rate of multiple users while minimizing inter-user interference.However, we face challenges due to the interdependence between the variables W and φ l,n , as well as the presence of a sum-of-log- of-ratio term, K k=1 log(1 + γ k ) .To address these challenges efficiently, we adopt a similar (8) approach as in [18], where we solve the optimized problem by introducing auxiliary variables.Additionally, we explore closed-form fractional programming techniques [32] to transform problem (P1) into a form with significantly reduced complexity.By employing the Lagrange dual transform, we can transform the problem (P1) while ensuring the equivalence of ̺ k = γ k as The vector ̺ = [̺ 1 , ̺ 2 , . . ..̺K ] T represents the auxiliary variable ̺ k for each user k .In problem (P2), by introducing the auxiliary variable ϱ and considering the QT technique, we can separate the summation-of-ratio log side and the summation of log as mentioned in Eq. ( 12).This decomposition greatly facilitates the attainment of the optimal solution.Before analyzing the JAPBF for the two subproblems, we focus on the auxiliary variable ϱ in (P2), assuming the variables W and φ 1,1 , . . ., φ L,Q are fixed.Consequently, we can optimize the problem Based on this proposition, we proceed to find the optimal solution by introducing and reformulating the two separate subproblems, denoted as W and φ 1,1 , . . ., φ L,Q .
(1) The analytical solution for the optimal active beamforming matrix W : Assuming that ϱ is held constant, our attention is directed toward the third term in Eq. ( 12), which contains the variable to be optimized.Based on this premise, we can reformulate the corresponding optimization problem as follows where As mentioned in (14), problem (P2.2) can be classified as a multiple-ratio fractional programming problem.To address this, we need to reformulate the problem as a biconvex optimization problem by introducing a new objective function that leverages the QT in [32], as where χ is introduced as a new auxiliary variable, where χ = [χ 1 , χ 2 , . . .., χ K ] T repre- sents the auxiliary variable vector after transforming the optimization problem into a biconvex problem.Consequently, we can solve the biconvex problem by considering W and χ as variables and addressing the following problem: finding the solution to (P2.2) while subject to the constraint of W .To deal with and analyze the problem in Eq. ( 15), we can investigate proposition 2 in [18], to represent the optimal beamforming vector w k .
(2) The analytical solution for the optimal passive beamforming matrix φ 1,1 , . . ., φ L,Q : The problem (P2) is transformed into a relaxed form, denoted as (P2.1), by introducing the auxiliary variable ̺ .Unlike, consider the active beamforming matrix W to be constant and concentrate on iteratively optimizing the reflection vector of each RIS.Our goal is to enhance the received signal at the user's end by utilizing each RIS as both a relay and a scatter point, thereby increasing diversity.However, this approach poses challenges when it comes to optimizing the RIS phase shift matrix.Therefore, we treat the RIS phase shift matrix as an objective variable in the channel user k .Given that the reformulated channel is represented as where φ s l,n is the updated row vector after using the selected matrix S k that contains the reflecting vectors of φ H l,n by Q−n + 1 which indicates the times of appearance of reflect- ing vector in case of individually RIS l,n optimization, i.e., φ l,n ∈ C 1×(Q−n+1)N e .Further- more, the matrix α l,n represents the portion of the channel that is associated with the optimized reflecting vector of the RIS l,n , where α l,n ∈ C (Q−n+1)N e ×M .On the other hand, the matrix β l,n represents the residual portion of the channel that is not directly linked to the optimized reflecting variable, denoted as β l,n ∈ C 1×M .Moving forward, we can focus on the part of the channel for user k , which includes the optimized reflecting vec- tor of the RIS l,n representing the route and path as follows: (15) where Once we have reformulated h H s,k for the passive beamforming scenario, the objective function in (13) can be represented as a function of φ l,n as where δ 1 , δ 2 , …, δ K are the newly introduced complex auxiliary variables and δ = [δ 1 , δ 2 , . . ., δ K ] T is the auxiliary variable vector.Subsequently, as demonstrated in [32], optimizing f 3.1 (φ l,n , δ) with respect to φ l,n is equivalent to solving the following problem with respect to and δ.By investigating proposition 3 in [18], we can analyze the problem in Eq. ( 22), to achieve the optimal expression φ l,n,n e of the reflecting vec- tor φ l,n is updated.In this way, we can ensure the convergence of the objective function, by updating the phase shift matrix of each RIS l,n iteratively.From the complexity point of view, however, we compensate for the network performance by using more routes, the proposed scheme reduces the complexity passive beamforming between deployed RISs.In this way, the convergence of the objective function is achieved by iteratively updating the phase shift of each RIS l,n .(18)

Results and discussion
To demonstrate the effectiveness of the proposed strategy and using γ T adaptively changed, we consider two different scenarios for user allocation with fixed RISs location, as shown in Figs. 3 and 4.This helps us understand how users impact the RIS-S scheme.The distribution of multiple RISs resembles the arms of an octopus, ensuring that the signal is injected toward the user's side.In the mmWave frequency band, we assume that there are N e equals 49 reflected elements per surface.Using the Friis transmission law [33], we calculate the channels G l and g l,k , with a Rician factor F of 6, and the power of noise is set to − 85 dBm.Table 1 provides a summary of the simulation parameters.
As shown in Fig. 3, case 1 illustrates a wide distance distribution of users.Based on the RIS selection, it is noted that the RIS 4 and RIS 6 suppressed.That is due to long distances between the surface to BS and UE, respectively.So, it means that these surfaces are not effective to involve in assisting users or enhance the communication network.On the contrary, the optimization complexity increases without noted utility.On the other hand, case 2 in Fig. 4 illustrates a short-distance distribution and the suppressed surfaces are RIS 5 and RIS 6 for the same reason.
Simulation results were simulated by MATLAB 2021b on Dell desktop 10 th g with RAM 16G and processor 3.75 GHz.We conducted an averaging across 300 independent realizations of wireless channels.For evaluating the proposed technique, we compared proposed strategy with one scheme employing RIS-S based on GNN [22] and two schemes: without RIS-S1 [10] and without RIS-S2 [18] which signify cooperation among all RISs with parallel and cascade topologies, respectively.
As shown in Figs. 5 and 6 the sum rate is measured in bit per channel use (bpcu) as a unit.Initially, the sum rate convergence of the RIS-S scheme appears superior to the benchmark schemes in two scenarios.However, while the proposed RIS-S scheme ultimately achieves better convergence compared to the GNN-based scheme, the GNN-based approach slightly outperforms at the sixth iteration.This variation arises from the fact that the GNN-based scheme conducts selection internally during the optimization process, unlike the proposed technique, which performs it as a preliminary stage.Ultimately, the RIS-S strategy yields the optimal sum rate value, even with changes in user locations.The selection matrix shifts are evident in Figs. 5 and 6.The Rician factor F 6 The power of noise σ 2 0 − 85 dBm The frequency of carrier Additionally, two schemes without selection are executed in a simplified design from a path perspective.Consequently, each user is served by all selected RISs, while weak reflection paths may not be necessary for optimization.As noted in Figs. 5 and 6 the convergence of all schemes stabilized after iteration 4 except GNN-based at iteration 6.
Figures 7 and 8 illustrate the influence of the number of elements concerning the selection strategy.It is widely acknowledged that augmenting the elements per surface improves the communication link.However, this improvement is constrained by factors such as surface location and orientation angle.Therefore, the most effective approach to enhance the communication link is to prioritize the selection process before expanding the surface area without notable benefits.As depicted in Figs.7 and  8, the proposed RIS-S technique surpasses all other schemes.Moreover, it is noteworthy that the proposed technique achieves a higher sum rate compared to GNN-based methods, contrary to previous findings in Figs. 4 and 5.This discrepancy stems from the fact that GNN-based approaches select the optimal surface among deployed RISs, whereas the proposed technique selects a group or multiple surfaces from candidate RISs.Fig. 8 Sum rate based on number of elements (N e ) for the proposed RIS-S scheme versus benchmark schemes (case 2) Fig. 9 The consumed time needed for optimal rate (case 1) Figures 9 and 10 depict the time required for each scheme to converge and attain the optimal active weights and phase shifts of RISs in two scenarios.The RIS-S scheme exhibits the shortest convergence time, attributed to the optimization process being conducted between S surfaces rather than N RISs, as seen in benchmark schemes.Conversely, the GNN-based approach, despite performing the selection role, consumes more time to achieve optimal rates due to the involvement of selection within the optimization process.Furthermore, the two schemes without selection require extended time to reach optimal rates at each iteration, as JAPBF is established among all RISs.Ultimately, the significance of RIS-S as a preliminary stage before JAPBF optimization is performed, emphasizing its pivotal role.
We estimate the complexity of the proposed system as introduced in [18].In Table 2, we note that two parameters that indicate the difference in comparison of complexity are the number of selected RISs ( S ) and the number of convergence itera- tions (I).Parallel and cascaded topology utilization idea is based on using all deployed RISs ( N ) without checking if the surface is effective or not.Because the main object of parallel topology is achieving more spatial diversity by utilizing each surface to reflect directly to the user, also the cascaded topology constructs a chain of cascaded surfaces to guide the signal to the user.On the other hand, the RIS-S strategy facilitates choosing the more effective surfaces ( S ) for all deployed surfaces ( N ).At the most case, if all surfaces ( N ) are effective and this case is rarely, selection strat- egy works as the parallel topology.In addition, the convergence iteration number is affected by reaching the optimal maximum sum rate of the system, in order to compute and estimate the computational complexity of the JAPBF algorithm.Likewise the work proposed in [18] which represented the computational complexity as the same method for showing the complexity reduction.We can divide the estimation of complexity based on the active BF and the passive BF.To begin, the complexity values are O(K ) and O M 3 K , respectively, as a result of updating the auxiliary's variable χ op k and calculating the complexity of W for the BS.As a result, the complexity of active BF is O M 3 K + K .Second, by calculating the auxiliary variable, the complexity equals O(K ) for all users K .Also, by iteratively updating the phase shift of each element of RIS l,n , the complexity equals O(N e ) .In the same way, the complexity of passive BF is estimated to be O(N e + K ) .When N RISs are used as multi-hop points, the complexity equals O( N n=1 N e + NK ) .Finally, we can estimate the overall computational complexity for JAPBF as O (M 3 + 2 + N )K + N n=1 N e I c where I c is the number of iterations needed to converge.However, our proposed model achieved the same computational complexity as [18], but the JAPBF algorithm converges and reaches the optimal solution faster.

Conclusion
We present a novel technique to joint active and passive beamforming (JAPBF) that differs from conventional methods.Instead of utilizing all deployed RISs, we introduce a new strategy for selecting RISs.This allows us to investigate the impact of spatial diversity gain by serving users with multiple RISs, while keeping the complexity low.By examining RIS-S before applying JAPBF, we can maximize the overall 8 achievable sum rate for users.The problem formulation is non-convex, so we employ fractional programming and quadratic transform techniques to relax and break down the problem into simpler subproblems.Through simulation results, we demonstrate that the JAPBF approach based on our proposed RIS-S strategy outperforms JAPBF approaches that rely on all RIS cooperation schemes.Additionally, our design significantly reduces computational complexity.
The following ideas are suggested for future research: • RIS selection in case of more obstructed channel as Nakagami-m distribution.
• Studying the performance of the proposed RIS selection from the energy efficiency and spectral efficiency points to verify how the impact of utilized RIS reduction on the power consumption.

Table 1
Simulation parameters