DOI QR코드

DOI QR Code

RIS Selection and Energy Efficiency Optimization for Irregular Distributed RIS-assisted Communication Systems

  • Xu Fangmin (Key Laboratory of Data Storage and Transmission Technology of Zhejiang Province, Institute of Communications Engineering, Hangzhou Dianzi University) ;
  • Fu Jinzhao (Key Laboratory of Data Storage and Transmission Technology of Zhejiang Province, Institute of Communications Engineering, Hangzhou Dianzi University) ;
  • Cao HaiYan (Key Laboratory of Data Storage and Transmission Technology of Zhejiang Province, Institute of Communications Engineering, Hangzhou Dianzi University) ;
  • Hu ZhiRui (Key Laboratory of Data Storage and Transmission Technology of Zhejiang Province, Institute of Communications Engineering, Hangzhou Dianzi University)
  • Received : 2023.04.11
  • Accepted : 2023.06.25
  • Published : 2023.07.31

Abstract

In order to improve spectral efficiency and reduce power consumption for reconfigurable intelligent surface (RIS) assisted wireless communication systems, a joint design considering irregular RIS topology, RIS on-off switch, power allocation and phase adjustment is investigated in this paper. Firstly, a multi-dimensional variable joint optimization problem is established under multiple constraints, such as the minimum data requirement and power constraints, with the goal of maximizing the system energy efficiency. However, the proposed optimization problem is hard to be resolved due to its property of nonlinear nonconvex integer programming. Then, to tackle this issue, the problem is decomposed into four sub-problems: topology design, phase shift adjustment, power allocation and switch selection. In terms of topology design, Tabu search algorithm is introduced to select the components that play the main role. For RIS switch selection, greedy algorithm is used to turn off the RISs that play the secondary role. Finally, an iterative optimization algorithm with high data-rate and low power consumption is proposed. The simulation results show that the performance of the irregular RIS aided system with topology design and RIS selection is better than that of the fixed topology and the fix number of RISs. In addition, the proposed joint optimization algorithm can effectively improve the data rate and energy efficiency by changing the propagation environment.

Keywords

1. Introduction

With the rapid growth of mobile users and new wireless applications, communication connections are becoming almost ubiquitous. This has led to higher demands on spectral efficiency and connectivity. Driven by these requirements, researchers have begun to develop the sixth generation (6G) wireless communication technology. Compared with 5G, 6G must greatly improve the communication network performance in terms of coverage, capacity, delay, and reliability, which results in extremely challenging energy efficiency and cost design [1-2].

In recent years, reconfigurable intelligent surface (RIS) technology has been considered as a potential candidate for 6G due to its ability to achieve programmable propagation environments with high energy efficiency and low hardware cost. Specifically, RIS is a reflective surface with integrated electronic circuits, which consists of a number of components controlled by programmable varactor diodes. Through intelligent control, RIS has the advantages of changing the propagation environment of electromagnetic waves, enriching the channel scattering conditions, solving signal propagation in the case of obstacles, and improving the safety of electromagnetic wave propagation [3-5]. The existing literature has studied RIS from the aspects of precoding, power allocation, phase shift adjustment and so on. The results show that RIS performs well in improving the spectral efficiency and coverage probability of communication systems. However, for the system model where multiple irregular RISs are deployed around the user, how to plan the RIS deployment, and how to select the appropriate RIS to provide services to the user and optimize the system energy efficiency have not been studied. Therefore, this paper focuses on the energy efficiency optimization problem of this novel system.

1.1 Related works

In order to further improve the performance of RIS-assisted communication systems, a large number of studies have been conducted in the existing literature on adjusting the RIS phase shift, transmission power, and beamforming [6-16]. Literature [6] combines the transmit power of the base station (BS) and phase adjustment to optimize the available rate of RIS-assisted communication system. In [7], the asymptotic performance of large RIS-based antenna array system is studied, the channel hardening effect is analyzed, and the performance bound is derived. [8] proposes a resource allocation algorithm taking into account the passive beamforming and modulation schemes. In [9], authors study the RIS-aided multi-user communication system, and propose corresponding resource allocation algorithms for the single-user scenario and the multi-user scenario respectively. Numerical results show that RIS improves system performance significantly. In [11], authors propose a coverage maximization algorithm and give the optimal RIS layout in closed form by analyzing the relationship between RIS layout and the coverage area of RIS. [15] analyzes the relationship between the performance and the number of RISs in a distributed RIS environment. It shows that although the increase of the number of RISs will improve the system capacity, the beamforming complexity, channel estimation overhead and channel feedback overhead will also be increased correspondingly, so the number of RISs should not be too large. However, the RIS components studied in these literatures are regularly deployed. [16] proposes a novel RIS topology of irregular array in which a given number of RIS elements are irregularly arranged on the surface of RIS. The choice of RIS element arrangement will increase the RIS degree of freedom and spatial diversity, thus improving the system capacity.

1.2 Motivation and contributions

Inspired by the advantages that the irregular configuration of phased array has, this paper studies the RIS on-off selection and energy efficiency optimization problem under the irregular RIS topology, which has not been studied in other literature so far.

Specifically, this paper first establishes a model of multiple irregular RIS-assisted downlink communication systems, and proposes an optimization problem with the objective of maximizing the energy efficiency (EE), which combines multiple-variable optimization such as RIS topology design, RIS switch selection, and power allocation. Since it is difficult to obtain the optimal solution of the joint optimization problem, we then decompose the optimization problem into the following sub-problems: topological structure optimization of RIS, RIS phase shift matrix, power allocation and RIS switch state selection sub-problems, so as to obtain the sub-optimal solution of system energy efficiency. Finally, the superiority of the proposed algorithm is verified by simulation analysis. The main work of this paper is summarized as follows.

1) We propose a distributed irregular RIS-assisted communication system, which can adapt to network requirements by designing RIS element topology structure and dynamically opening or closing RISs. In order to maximize the energy efficiency of the system, we jointly optimize RIS element topology, the phase shift of RIS, transmit power and RIS switching state vector. Under the constraint of user minimum rate, transmit power and RIS phase shift unit modulus, we propose an optimization problem aiming at energy efficiency maximization.

2) To deal with the proposed non-convex optimization problem, we decompose the problem into four sub-problems of low complexity. Then, the optimal solution is obtained by solving the four joint sub-problems iteratively.

3) We apply a modified Tabu search algorithm to obtain the solution of the topology design sub-problem. The phase shift is optimized by SCA, and power allocation in derived in a close form. For RIS switch optimization sub-problem, the greedy algorithm is used to find the optimal solution.

4) To evaluate the performance, we prove the convergence of the proposed algorithm. In addition, the complexity is analyzed theoretically and EE is evaluated by simulations.

2. System Model and Question Description for RIS-assisted Communication Systems

2.1 System Model

This paper investigates a novel downlink system model of multiple irregular RIS-assisted communication systems, as shown in Fig. 1. In the system, there is one base station (BS) with 𝑀 antennas and L irregular RIS evenly distributed 𝑑0 meters away from a single antenna user, among which the 𝑁𝑁𝑙𝑙 reflection elements of the l th irregular RIS are sparsely distributed on the grid points of RIS. The small yellow rectangle of RISs in the figure represents the reflection element that works, and if all elements work, it’s a common regular RIS. Here we focus on how to select RIS when there are many RIS around the user and how to design the topology of the reflecting element on RISs to optimize system energy efficiency.

E1KOBZ_2023_v17n7_1823_f0001.png 이미지

Fig. 1. The downlink of Irregular RIS-aided communication systems (the on-off status of RIS is controlled by BS)

As shown in Fig. 1, the received signals of the user equipment (UE) consists of the signal from direct link and the signals from reflected link. Assuming that the BS has the channel status information of all channels [17-18], the signal received at UE can be expressed as:

y = (ℎ𝐻𝑑 + ∑𝐿𝑙=1𝑥𝑙𝐻𝑙𝑍𝑙Θ𝑙𝐺𝑙)𝑤s + σ2       (1)

where 𝑤 ∈ ℂ𝑀 represents the transmitted precoding vector from BS, σ2 is the variance of AWGN, 𝑠 ∈ ℂ is the transmitted signal, ℎ𝑙 ∈ ℂ𝑁𝑆 is the channel vector between the irregular RIS l and the UE, 𝐺𝑙 ∈ ℂ𝑁𝑆×𝑀 is a 𝑁𝑆 × 𝑀 channel matrix between BS and the l th RIS, ℎ𝑑 ∈ ℂ𝑀 is the channel from BS to UE, Θ𝑙 = diag(𝜃𝑙,1, … 𝜃𝑙,𝑛 … , 𝜃𝑙,𝑁𝑆) is the phase shift matrix at the lth RIS where θ𝑙,𝑛 = 𝑒𝑙,𝑛 denotes the reflection coefficient of the nth elements of RIS l, 𝜑𝑙,𝑛 is the phase shift. Γ = {1, … , 𝐿} and θ𝑙 = [𝜃𝑙,1, … 𝜃𝑙,𝑛 … , 𝜃𝑙,𝑁𝑆]𝑇 are the index set of RISs and phase shifts of RIS l respectively. ℂ denotes the set of complex numbers. In order to reduce power consumption, some RISs are not always selected as a serving RIS for the UE. Hence, we introduce the variable 𝑥𝑙 ∈ {0,1} here to indicate whether the corresponding RIS is selected. That is to say, 𝑥𝑙=1 indicates that the l th RIS is selected, otherwise, it is not selected. 𝑍𝑙 = diag(𝐳𝑙) represents the topology selection matrix of the lth irregular RIS with 𝐳𝑙 = [𝑧𝑙,1, … 𝑧𝑙,𝑛 … , 𝑧𝑙,𝑁𝑆]𝑇. If the RIS element is deployed at the ith grid point of the l irregular RIS, then 𝑧𝑙,𝑖 = 1, otherwise, 𝑧𝑙,𝑖 = 0.

Assuming that the BS adopts the maximum ratio transmission, the data rate and energy efficiency can be expressed respectively as:

\(\begin{aligned}R_{t}=B \log _{2}\left(1+\frac{p_{\text {down }}\left|h_{d}^{H}+\sum_{l=1}^{L} x_{l} h_{l}^{H} Z_{l} \operatorname{diag}\left(\theta_{l}\right) G_{l}\right|^{2}}{\sigma^{2}}\right)\end{aligned}\)       (2)

and

\(\begin{aligned}\mathrm{EE}=\frac{R_{t}}{\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} x_{l} N_{l} P_{R}}\end{aligned}\)       (3)

where B is bandwidth, 𝑝down the transmission power of the BS. ρ = 𝜏−1𝑝down + 𝑃𝑈 + 𝑃𝐵S + ∑𝐿𝑙=1𝑥𝑙𝑁𝑙𝑃𝑅 is the power consumption, 𝜏−1 is transmitted power amplifier efficiency of BS. 𝑃𝐵S, 𝑃𝑈, and 𝑃𝑅 are the hardware static power consumption at BS, UE and RIS, respectively.

2.2 Question formulation

Based on the above analysis, we focus our attention on how to maximize the energy efficiency of the system by considering the optimal design of the RIS location, the topology structure of the reflective components, and the phase shift. Hence, the optimization problem is expressed as follows

\(\begin{aligned}\max _{\boldsymbol{\theta}, p_{\text {down }}, x, Z} E E=\frac{R_{t}}{\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} x_{l} N_{l} P_{R}}\end{aligned}\)       (4)

s.t.

𝐶1 : 𝐵log2(1 + 𝑝down𝛾) ≥ 𝑅min,

𝐶2: 𝑝down ≤ 𝑃max,

𝐶3: 𝜑𝑙,𝑛 ∈ [0,2𝜋]; 𝑙 = 1, … , 𝐿; 𝑛 = 1, … , 𝑁𝑆

𝐶4: 𝑥𝑙 ∈ {0,1}; 𝑙 = 1, … , 𝐿,

𝐶5: 𝑧𝑙,𝑛(𝑧𝑙,𝑛 − 1) = 0; 𝑙 = 1, … , 𝐿; 𝑛 = 1, … , 𝑁𝑆

𝐶6 : ∑𝑛 𝑧𝑙,𝑛 = 1; 𝑙 = 1, … , 𝐿;

where 𝛉 = [θ1, … θ𝑙 … , θ𝐿], \(\begin{aligned}\gamma=\frac{\left|h_{d}^{H}+\sum_{l=1}^{L} x_{l} h_{l}^{H} Z_{l} \operatorname{diag}\left(\theta_{l}\right) G_{l}\right|^{2}}{\sigma^{2}},\end{aligned}\), 𝐶1 ensures the UE’s minimum data rate, 𝐶2 is the transmission power constraint of the BS, 𝐶4 indicates whether RIS l is selected, 𝐶5 and 𝐶6 are sparse deployment constraints of irregular RISs where N diagonal elements of the matrix Z are assigned 1. Because the constraints 𝐶3, 𝐶4, 𝐶5 and 𝐶6 are non-convex, the optimization problem (4) is hard to solve directly.

3. Energy Efficiency Optimization Strategy

To tackle the above-mentioned issue, we decompose the optimization problem (4) into 4 sub-problems. Specifically, we first optimize topology matrix Z when the phase shift, power, and on-off variable of RISs are fixed. Then, we employ Taylor expansion to optimize of phase shift and power allocation with other variables fixed. Furthermore, we select APs by greedy algorithm to maximize the EE. Finally, we iteratively update the solutions until convergence.

3.1 RIS reflection element topology optimization based on the modified Tabu search

Tabu search, as a low-complexity strategy for combinatorial optimization problems, has received a lot of attention in numerous areas, such as scheduling, spatial planning, and graph coloring. The Tabu search algorithm starts with an initial solution vector and finds a local neighborhood around it. It selects the best vector among neighboring solutions even if the current solution is the worst. Therefore, it is prone to deviate from the local optimal value. Motivated by the success applications of Tabu search, a modified Tabu search based optimization for the topology of the RIS is proposed. Different from the traditional Tabu search, in the proposed algorithm, RIS elements are not arranged in the two adjacent grid points and the moving distance is adjusted dynamically according to the number of the system parameters, such as reflection elements of RIS, the number of iterations, etc. Before introducing the proposed algorithm, we introduce some definitions.

Definition of moving distance: Define p as the number of 0 and 1 interchanges per time, which will be dynamically adjusted. The larger p is, the faster the exchange is, and the smaller p is, the slower the exchange is. Therefore, if a smaller search range is needed, p is initialized to a smaller number and then adjusted according to the system parameters. Conversely, p is initialized to the larger number.

Definition of Neighborhood: The neighborhood of a given topology vector z𝑙 is defined as a set 𝑁(z𝑙). If the cardinality of 𝑁(z𝑙) is Q, it has Q neighboring vectors. Each neighboring vector has 2 p different value from the former one.

Tabu list: Tabu list records the previous solutions and avoids getting stuck in a local optimum. Any solution worse than the elements in the Tabu list will not be adopted. In addition, if the Tabu list is full, the first vector in the list will be removed.

Next, we propose the steps of topology planning based on the modified Tabu list as follows:

Step 1: We initialize the topology of each RIS l as z(0)𝑙 with randomly select N discontinuous ones, and set the remaining 𝑁𝑆 − 𝑁 elements to 0. Then, the energy efficiency is calculated. To provide better diversity performance, we add all cases with N consecutive ones directly to the Tabu list.

Step 2 starts the iteration.

Step 3: We generate a neighbor list with Q neighbors by exchanging p zeros and p ones to obtain the topology of RIS denoted by z(𝑖)𝑙,1, … , z(𝑖)𝑙,𝑄. Then, calculate the energy efficiency of all candidates according to equation (3).

Step 4: The candidate with the largest energy efficiency is then selected, i.e., \(\begin{aligned}\mathrm{z}_{l, *}^{(i)}=\underset{q}{\operatorname{argmax}}\left\{E E\left(\mathrm{z}_{l, q}^{(i)}\right)\right\}\end{aligned}\), which is saved as a new topology for the next iteration. That is to say, in the neighbor list, we select the topology with the highest energy efficiency as z(𝑖+1)𝑙,1, and add z(𝑖)𝑙,∗ to the tabu list to avoid circular search. If the list is full, we delete the earliest Tabu solutions in the current list.

Step 5: The optimization value and Tabu list will be updated accordingly. After reaching the iteration threshold, we can obtain the history optimal topology of RIS.

The aforementioned topology optimization of RIS is briefly described in Algorithm 1.

Algorithm 1 Topology optimization of RIS

Input: M, N, 𝑁𝑆, L, the number of neighbors 𝑄, moving distance p, the number of iterations T, the length of Tabu list.

Output: Topology 𝑍0

1) Initialize the topology of each RIS z(0)𝑙, Tabu list (the cases with N consecutive ‘one’ are added to the initial Tabu list), randomly select N ones that are not adjacent, 𝑁𝑆 − 𝑁 zeors, the number of neighbors 𝑄, and moving distance p.

2) for i=1:T

3) Adjust p adaptively according to i, exchange p zeros and p ones that are not adjacent to change the topology of RIS, and generate a neighbor list z(𝑖)𝑙,1, … , z(𝑖)𝑙,𝑄.

4) In the neighbor list, select the topology with the highest energy efficiency as z(𝑖+1)𝑙,1.

5) Add \(\begin{aligned}\mathrm{z}_{l, *}^{(i)}=\underset{q}{\operatorname{argmax}}\left\{E E\left(\mathrm{z}_{l, q}^{(i)}\right)\right\}\end{aligned}\) to the Tabu list and update the Tabu list.

6) end for

7) Output: The history optimal topology of RIS l.

3.2 Optimization for phase shift and transmit power

By fixing RIS topology and RIS switching variable x, we jointly optimize subproblem of RIS phase coefficient and transmitting power to obtain the suboptimal solution of the objective function. Since the denominator of optimization problem (4) has nothing to do with phase optimization, the phase optimization problem can be transformed into the following problem:

\(\begin{aligned}\max _{\boldsymbol{\theta}}\left|h_{d}^{H}+\sum_{l=1}^{L} x_{l} h_{l}^{H}\left(\theta_{l}\right)^{T} U_{l}\right|^{2}\end{aligned}\)       (5)

𝑠.𝑡.𝜑𝑙,𝑛 ∈ [0,2𝜋]; 𝑙 = 1, … , 𝐿; 𝑛 = 1, … , 𝑁S

where (𝜃𝑙)𝑇𝑈𝑙 = ℎ𝐻𝑙𝑍𝑙diag(𝜃𝑙)𝐺𝑙, 𝑈𝑙 = diag(ℎ𝐻𝑙𝑍𝑙)𝐺𝑙.

In order to make RIS selection variable x, topology variables 𝑍𝑙 and the corresponding phase coefficient separation, we denote 𝑼 = [𝑥1𝑈1, … , 𝑥𝐿𝑈𝐿] ∈ ℂ𝑄×𝑀 and 𝒗 = [𝜃1*, … , 𝜃𝐿*] ∈ ℂ𝑄. Then, problem (5) reduces to the optimization problem of variable 𝑼 and variable 𝒗:

\(\begin{aligned}\max _{v}\left|h_{d}^{H}+\boldsymbol{U}^{H} \boldsymbol{v}\right|^{2}\end{aligned}\)       (6)

s.t.,

|𝑣𝑞| = 1; 𝑞 ∈ {1, … ,𝑄}

where 𝑄 = 𝐿𝑁𝑆.

In order to solve problem (6), the first-order Taylor expansion of its objective function is carried out,

|ℎ𝐻𝑑 + 𝑼𝐻𝒗|2 ≈ 2𝑅((ℎ𝐻𝑑 + 𝑼𝐻𝒗𝒏−𝟏)𝐻𝑼𝐻𝒗) + |ℎ𝐻𝑑 + 𝑼𝐻𝒗𝒏−𝟏|2 − 2𝑅((ℎ𝐻𝑑 + 𝑼𝐻𝒗𝒏−𝟏)𝐻𝑼𝐻𝒗𝒏−𝟏)       (7)

Then, by relaxing the constraints, problem (6) is transformed into a solvable convex problem according to (7) as follows:

\(\begin{aligned}\max _{\boldsymbol{v}} 2 R\left(\left(h_{d}^{H}+\boldsymbol{U}^{H} \boldsymbol{v}^{\boldsymbol{n}-\mathbf{1}}\right)^{H} \boldsymbol{U}^{H} \boldsymbol{v}\right)\end{aligned}\)       (8)

s.t.

|𝑣𝑞| ≤ 1; 𝑞 ∈ {1, … ,𝑄}

Therefore, if the following conditions are satisfied:

𝐼m([(ℎ𝐻𝑑 + 𝑼𝐻𝒗𝒏−𝟏)𝐻𝑼𝐻]𝑞)[𝒗]𝑞 = 0; 𝑞 ∈ {1, … , 𝑄}

is a real number, and [𝒗]𝑞=1, then 𝒗 has an optimal solution:

𝒗 = 𝑒𝑗∡𝑼𝐻(ℎ𝐻𝑑+𝑼𝐻𝒗𝒏−𝟏)       (9)

Here, [∙]𝑞 is the q th element of [∙]. The target value obtained by the above algorithm is monotone non-decreasing, and the iteratively obtained value 𝒗𝑛 converges to the point satisfying the Karush-Kuhn-Tucker (KKT) optimal condition of the original non-convex problem (6).

In order to optimize the power, we substitute the optimized phase (9) into formula (4) and further simplify the optimization problem as follows:

\(\begin{aligned}\max _{p_{1}} E E=\frac{B \log _{2}\left(1+p_{\text {down }} \gamma\right)}{\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} x_{l} N_{l} P_{R}}\end{aligned}\)       (10)

s.t.

𝐶1:𝐵log2(1 + 𝑝down𝛾) ≥ 𝑅min

𝐶2: 𝑝down ≤ 𝑃max

Since the objective function of Equation (10) is upper convex, the point where its derivative equals 0 is the power point that maximizes the objective function. We take the derivative of the objective function of problem (10) with respect to pdown as

\(\begin{aligned}\frac{\partial \frac{B \log _{2}\left(1+p_{\text {down }} \gamma\right)}{\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} x_{l} N_{l} P_{R}}}{\partial p_{\text {down }}}=\frac{B\left(\gamma\left(\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} x_{l} N_{l} P_{R}\right)-\tau^{-1} \cdot\left(1+p_{\text {down }} \gamma\right) \ln \left(1+p_{\text {down }} \gamma\right)\right)}{\ln 2 \cdot\left(1+p_{\text {down }} \gamma\right)\left(\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} x_{l} N_{l} P_{R}\right)^{2}}\end{aligned}\)       (11)

The optimal power can be obtained by solving the following formula

𝛾(𝜏−1𝑝down + 𝑃𝑈 + 𝑃BS + ∑𝐿𝑙=1𝑥𝑙𝑁𝑙𝑃𝑅) − 𝜏−1 ∙ (1 + 𝑝down𝛾)ln(1 + 𝑝down𝛾)=0       (12)

Let 𝑡 = 1 + 𝑝down𝛾, and from (12), we have

𝑓(𝑡) = 𝑡(lnt − 1) = 𝛾𝜏(𝑃𝑈 + 𝑃BS + ∑𝐿𝑙=1𝑥𝑙𝑁𝑙𝑃𝑅)       (13)

From (13), we have

𝑝down = [𝑓−1(𝛾𝜏(𝑃𝑈 + 𝑃BS + ∑𝐿𝑙=1𝑥𝑙𝑁𝑙𝑃𝑅)) − 1]/𝛾       (14)

where 𝑓−1() is the inverse function of 𝑓(𝑡).

3.3 RIS selection and switching optimization

After optimizing the topology structure, power and phase shift, the optimization objective function (4) is simplified to

\(\begin{aligned}\max _{x} E E=\frac{B \log _{2}\left(1+\frac{p_{\text {down }}\left|h_{d}^{H}+\sum_{l=1}^{L} x_{l} h_{l}^{H} z_{l} \operatorname{diag}\left(\theta_{l}\right) G_{l}\right|^{2}}{\sigma^{2}}\right)}{\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} x_{l} N_{l} P_{R}}\end{aligned}\)       (15)

However, (15) is a nonlinear integer problem and can’t be solved directly. Therefore, we adopt greedy algorithm (shown in Algorithm 2) to calculate the corresponding objective function by trying to deactivate each RIS. If the disable of one RIS helps to improve the energy efficiency of the system, it will be closed; otherwise, it will not be closed. The RIS selection and switching optimization steps are described as follows:

Sept 1: We first initialize the switch vector of RISs to 𝑥 = (𝑥1, 𝑥2, … , 𝑥𝐿)𝑇 = 1𝐿, denote the set of activated RIS as Γ′, and calculate the corresponding energy efficiency as

\(\begin{aligned}E E_{0}=\frac{B \log _{2}\left(1+p_{\text {down }} \gamma\right)}{\tau^{-1} p_{\text {down }}+P_{U}+P_{B S}+\sum_{l=1}^{L} N_{l} P_{R}}\end{aligned}\)

Step 2: We calculate the energy efficiency 𝐸𝐸𝑙 by deactivating RIS 𝑙 ∈ Γ′ according to (15)

𝐸𝐸𝑙 = EE|𝑥𝑙=0

Step 3: If switching off RIS 𝑙*makes the system most energy efficient, turn it off. Then, update initial energy efficiency 𝐸𝐸0 and the set Γ′ for the next iteration.

Step 4: If switching off RIS does not improve EE, break and output the optimized EE and optimized vector 𝑥.

Algorithm 2: RIS selection and switching optimization based on greedy algorithm.

Input: 𝒁0, 𝛉, 𝑝down

1) Initialize the switch vector of RISs, i.e. 𝒙 = 𝟏𝐿, denote the initial energy efficiencyas 𝐸𝐸0; Denote the set of active RIS as Γ′.

2) while the cardinality of Γ′ >1 do;

3) for each active RIS 𝑙 ∈ Γ′;

4) Deactivate RIS l, i.e. 𝑥𝑙 = 0, and calculate the energy efficiency 𝐸𝐸𝑙 by (15);

5) If 𝐸𝐸𝑙 − 𝐸𝐸0 > 𝜀 > 0 , and 𝑙* = argmax{𝐸𝐸𝑙}, let 𝑥𝑙* = 0. Update theinitial energy efficiency 𝐸𝐸0 and the set of active RIS Γ′.

6) else 𝑥𝑙* = 1 and break.

7) end for;

8) end while

9) Output: The optimized EE and optimized RIS switch vector x.

3.4 The framework of the joint optimization algorithm

According to the previous topology optimization, phase and power optimization, and RIS switch optimization, we obtain the overall framework of the algorithm in this paper. The specific optimization scheme is given in Algorithm 3. In this framework, the variables are alternately updated until the convergence of the objective function is achieved. In addition, since the objective function in (4) is increasing at each iteration and the target value is bounded, this ensures the convergence of algorithm 3.

Algorithm 3: The Proposed Joint optimization Framework

Input: Initialize RIS phase element 𝛉𝟎, transmit power 𝑝0down, RIS topology 𝒁0 and RIS switching vector 𝑥0; Set m=1;

Step 1: fix the phase 𝛉𝑚−1, transmit power 𝑝𝑚−1down, RIS switching vector 𝑥𝑚−1, and optimized RIS topology of the (m-1) th iteration, and optimize the matrix of the mth topology 𝒁𝑚.

Step 2: fix the RIS topology 𝒁𝑚, transmit power 𝑝𝑚−1down, and RIS on-off vector 𝑥𝑚−1, and optimize the phase element of RIS 𝛉𝑚.

Step 3:fix the RIS topology 𝒁𝑚, RIS on-off vector 𝑥𝑚−1, 𝛉𝑚, and optimize the transmit power 𝑝𝑚down under the maximum transmit power constraint ;

Step 4:fix the RIS topology 𝒁𝑚, the optimized transmit power 𝑝𝑚down, 𝛉𝑚, and optimize 𝑥𝑚 (which indicates whether a RIS is enabled or not);

m=m+1;

Repeat the above steps until EE converges;

Output: 𝒁opt, 𝑝optdown, 𝛉opt, 𝑥opt and the optimized EE.

3.5 Complexity analysis

The complexity of Algorithm 3 mainly depends on the updates of the three variables 𝑍, θ, and 𝑥. Specifically, according to Algorithm 1, the complexity of updating the topology 𝑍 is 𝑂(𝑇1𝑄L). According to (9), the complexity of solving θ is 𝑂(𝑇2𝐴M), where A is the total number of elements of the L RISs. According to Algorithm 2, the complexity of solving 𝑥 mainly depends on calculating the objective function in (15) and the cardinality of Γ. Hence, the complexity of this part is O(𝐿(𝐿 − 1)𝐴M), where 𝑂(𝐴M) represents the complexity of calculating the EE in (15), and O(𝐿(𝐿 − 1)) is the iteration number of Algorithm 2. 𝑇1 and 𝑇2 are the number of iterations for convergence of topology and phase shift optimization, respectively. We ignore the complexity of power allocation since it is a closed-form solution. In summary, the complexity of our proposed algorithm is 𝑂(𝑇1𝑄L + 𝑇2𝐴M + 𝐿(𝐿 − 1)𝐴M). By contrast, the optimization of topology structure by exhaustive method has a higher complexity, i.e., \(\begin{aligned}L C_{N_{S}}^{N}=\frac{L N_{S} !}{N !\left(N_{S-N}\right) !}\end{aligned}\). The complexity of optimization 𝑥 by exhaustive method is ∑𝐿𝑙=1𝐶𝑙𝐿.

4. Simulation and Analysis

4.1 Parameters Setup

In this section, we consider a 200-meter square area with an 8-antenna BS in the center and one user randomly distributed in the area, as shown in Fig. 2. The performance of the proposed algorithm will be analyzed and compared with the following two algorithms: a centralized regular RIS (CRIS) scheme of RIS: the joint optimization of phase and power with one centralized regular RIS located at a distance of 𝑑0 meters from the user; Distributed regular RIS scheme (DRIS) scheme: the joint optimization of AP selection, phase and power with L regular RISs evenly distributed 𝑑0 meters away from the user. The proposed algorithm: the joint optimization proposed in Section 3 with L irregular distributed RISs deployed in the system.

E1KOBZ_2023_v17n7_1823_f0002.png 이미지

Fig. 2. Simulation scenario with one user randomly distributed in the 200-meter square area

The channel gain is determined by both large-scale fading and small-scale fading [19]. The Rayleigh fading channel model is used to account for the small-scale fading, i.e. ℎ𝑙, ℎ𝑑, 𝐺𝑙~𝐶N(0,1). The large-scale fading model is based on the distance between the transmitter and the receiver which is calculated by \(\begin{aligned}\beta=\frac{\gamma_{0}}{d^{a}}\end{aligned}\), where 𝛾0 is the reference pathloss at one-meter distance. 𝑎 is the path loss exponent. We set \(\begin{aligned}\beta=\frac{10^{-3.53}}{d^{3.76}}\end{aligned}\) for the following simulations.

In our simulation, the irregular RIS has 𝑁𝑆 = 20 points over its surface and N elements. In the distributed RIS case, each regular RIS has N elements, and in the centralized RIS case, each regular RIS has LN elements. Transmit power amplifier efficiency 𝜏−1 = 1/8, Bandwidth B=1MHz. 𝑃𝑅 = 10dBm, 𝑃𝐾 = 10𝑑𝐵m. The maximum transmit power of the BS is 𝑃max = 50dBm.

4.2 Performance Analysis

Fig. 3 describes the relationship between the average EE and the maximum transmission power. It shows that the EE of the proposed algorithm is superior to other algorithms. Compared with centralized RIS schemes, the proposed algorithm and distributed algorithm with RIS selection improve the energy efficiency by about 9.2% and 4.7% respectively when 𝑃max = 30dBm. Furthermore, as observed, the EE of the system increases first and then becomes stable with the increase of transmission power. This is because the data rate cannot keep up with the transmission power growth rate.

E1KOBZ_2023_v17n7_1823_f0003.png 이미지

Fig. 3. EE v.s. the maximal transmission power (𝑑0 = 20𝑚, 𝑁 = 4, 𝐿 = 8)

Fig. 4 describes the CDF of EE with different schemes. As observed, the probability that the user is in low energy efficiency is much less than the other two cases.

E1KOBZ_2023_v17n7_1823_f0004.png 이미지

Fig. 4. the CDF of EE (𝑑0 = 20𝑚, 𝑁 = 4, 𝐿 = 8, 𝑃max = 30)

Fig. 5 shows the relationship between the sum-rate versus the maximum BS transmission power. It can be seen that as the maximum transmit power increases, the sum-rate increases. Compared with the CRIS scheme, DRIS improves sum-rate effectively. This is due to the advantages of distributed deployment. Deploying multiple RIS in a DRIS network increases the received signal path for both user-RIS links and RIS-BS links, which provides greater diversity gain and thus increases data rate. However, in the case of CRIS, there is only one RIS, which results in a single path direction for user-RIS and RIS-BS, providing a less effective path diversity. Fig. 5 also shows that the sum-rate of the algorithm combined with irregular RIS is better than the case with regular RIS.

E1KOBZ_2023_v17n7_1823_f0005.png 이미지

Fig. 5. Sum-rate v.s. the maximal transmission power (𝑑0 = 20𝑚, 𝑁 = 4, 𝐿 = 8)

Fig. 6 gives the relationship between the user's minimum rate requirement and the EE. As observed, with the increase of the minimum rate requirements, the EE first stabilizes and then decreases. This is because, when the minimum rate threshold is relatively small, the system can easily achieve the rate requirement without additional increase in transmit power. However, when the threshold is large, the transmitted power needs to be increased to meet the minimum rate requirement, and the contribution brought by the rate increase is not enough to offset the increase in power consumption, so the energy efficiency is reduced. In addition, and the proposed algorithm is significantly better than other algorithms. This comes from the fact that the power consumption of the two distributed RIS is much lower than the transmission power of the centralized RIS, and the achievable rate is also lower than that of the centralized RIS. Therefore, the EE of the two distributed RIS is obviously better than that of the centralized RIS. Meanwhile, the proposed algorithm optimizes RIS topology structure and effectively improves EE.

E1KOBZ_2023_v17n7_1823_f0006.png 이미지

Fig. 6. the EE v.s. the minimal demand of date-rate (𝑑0 = 20𝑚, 𝑁 = 4, 𝐿 = 8)

Fig. 7 provides the relationship between the EE and the number of reflecting elements of RIS (N) in the 8-RIS scenario. It shows that as N increases, the EE of the proposed algorithm first increases and then decreases. This is because when N is relatively small, the contribution from the rapid growth of the system rate compensates for the increase in power consumption as N increases, which leads to the increase in EE. However, when N is relatively large, the increase in the system rate cannot catch up with the increase in power consumption, which leads to a decrease in EE. Meanwhile, the system complexity and deployment cost are also increasing with the increase of N. Therefore, the number of RIS elements should not be blindly increased in practice. Note that N here represents the number of active RIS elements, not the number of optional RIS primitives. Otherwise, EE will increase first and then becomes stable.

E1KOBZ_2023_v17n7_1823_f0007.png 이미지

Fig. 7. the number of reflecting elements of RIS v.s. EE (𝑑0 = 20𝑚)

Fig. 8 gives the relationship between EE and the user-RIS distance with different maximal power. It shows that energy efficiency decreases first and then increases as distance d0 increases. This is because the greater the user-RIS distance, the lower the data rate. Later, as the distance further increases, the RIS-BS link enhancement compensates for the loss of EE, so the EE increases.

E1KOBZ_2023_v17n7_1823_f0008.png 이미지

Fig. 8. EE v.s. the user-RIS distance d0

Fig. 9 depicts the effect of the number of RISs (L) on EE. As observed, EE increases with the increase of L. This is because more RISs provide more freedom of selection. Furthermore, when L continues to increase, Algorithm 2 will help select an appropriate number of RISs from L RISs to provide service for the user. Therefore, EE finally tends to be flat here. Compared with Fig. 7 and Fig. 9, EE shows different trends. This is because the N in Fig. 7 represents the number of RISs reflection elements participating in the service, and the higher the N, the higher the power consumption, which leads to the decrease of EE. However, L in Fig. 9 represents the number of the optional RISs. RIS will not be selected if there are already enough activated RISs, so EE will not decrease after reaching the maximum value.

E1KOBZ_2023_v17n7_1823_f0009.png 이미지

Fig. 9. EE v.s. the number of RISs L

Fig. 10 (a) compares the convergence performance of topology designs of different algorithms. It shows that the modified Tabu search has a larger initial value. This is because, we initialize to avoid starting from the possible case of low diversity gain. As observed, the modified Tabu search based topology design has better convergence than the traditional Tabu search. It can be seen that after a certain number of iterations, the modified algorithm changes the step size to accelerate the convergence. Fig. 10 (b) gives the convergence performance of the alternate iteration in Algorithm 3. It shows that the proposed algorithm converges after about 4 iterations.

E1KOBZ_2023_v17n7_1823_f0010.png 이미지

Fig. 10. EE v.s. the number of iterations (𝑁𝑆 = 12, 𝑁 = 4)

In Fig. 11, we activate 4 of the 12 RISs elements by exhaustive search, with a total of 495 possibilities, of which the maximum is 2.94 Mbits/Joule. By comparing Fig. 10 (a), it can be concluded that the algorithm proposed in this paper achieves near-optimal performance with low complexity.

E1KOBZ_2023_v17n7_1823_f0011.png 이미지

Fig. 11. the maximum value of exhaustive search (𝑁𝑆 = 12, 𝑁 = 4)

5. Conclusion

In this paper, for multiple irregular RIS-assisted communication systems, we investigate the joint optimization problem with the goal of maximizing energy efficiency under the minimum user data-rate constraint. However, the joint optimization problem in terms of 4-dimensional variables considering RIS topology structure, RIS switch, phase shift and power control is a nonlinear non-convex optimization problem and is very hard to be solved. To this end, we first decompose the hybrid optimization problem into four sub-problems. Then, the four sub-problems are iteratively optimized until convergence. The simulation results show that the proposed irregular RIS assisted communication system improves the system energy efficiency and sum-rate effectively.

6. Future research

In the communication scenario where multiple RISs are deployed around the user and RIS elements are irregularly distributed, this paper establishes the system model and proposes the optimization framework of energy efficiency, which is beneficial for the application of RIS in the sixth generation (6G). However, there are still some shortcomings in the paper, which can be further studied from the following aspects:

1) The research in this paper is based on the assumption of perfect CSI, but in the practical communication systems, there are always errors in channel estimation. Therefore, the optimization of imperfect CSI will be a part of our work in the future.

2) In this paper, we focus on the optimization for the single cell scenario, where the interference between the cells is ignored. Optimization research in multi-cell scenario will also be the focus of future research.

3) Recently, with the development of metasurfaces, the concept of simultaneous transmit and reflect RIS (STAR-RIS) has been proposed which has ability of providing full-space coverage in wireless environments. In our future work, EE optimization for irregular STAR-RIS will also be considered.

References

  1. Saad W, Bennis M, Chen M, "A Vision of 6G Wireless Systems: Applications, Trends, Technologies, and Open Research Problems," IEEE Network, vol. 34, no. 3, pp. 134-142, May/June 2020. https://doi.org/10.1109/MNET.001.1900287
  2. Mao B, Tang F, Kawamoto Y, et al, "AI Models for Green Communications Towards 6G," IEEE Communications Surveys & Tutorials, vol. 24, no. 1, pp. 210-247, Firstquarter 2022. https://doi.org/10.1109/COMST.2021.3130901
  3. J. Rao et al., "A Novel Reconfigurable Intelligent Surface for Wide-Angle Passive Beamforming," IEEE Transactions on Microwave Theory and Techniques, vol. 70, no. 12, pp. 5427-5439, Dec. 2022. https://doi.org/10.1109/TMTT.2022.3195224
  4. Renzo M D, Zappone A, Debbah M, et al, "Smart radio environments empowered by reconfigurable intelligent surfaces: How it works, state of research, and the road ahead," IEEE Journal on Selected Areas in Communications, vol. 38, no. 11, pp. 2450-2525, Nov. 2020. https://doi.org/10.1109/JSAC.2020.3007211
  5. Shen H, Xu W, Gong S, et al, "Secrecy rate maximization for intelligent reflecting surface assisted multiantenna communications," IEEE Communications Letters, vol. 23, no. 9, pp. 1488-1492, Sept. 2019. https://doi.org/10.1109/LCOMM.2019.2924214
  6. Huang C, Zappone A, Debbah M, et al, "Achievable Rate Maximization by Passive Intelligent Mirrors," in Proc. of 2018 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Calgary, AB, Canada, pp. 3714-3718, 2018.
  7. Jung M, Saad W, Jang Y, et al, "Performance Analysis of Large Intelligent Surfaces (LISs): Asymptotic Data Rate and Channel Hardening Effects," IEEE Transactions on Wireless Communications, vol. 19, no. 3, pp. 2052-2065, March 2020.
  8. Jung M, Saad W, Debbah M, et al, "On the Optimality of Reconfigurable Intelligent Surfaces (RISs): Passive Beamforming, Modulation, and Resource Allocation," IEEE Transactions on Wireless Communications, vol. 20, no. 7, pp. 4347-4363, July 2021.
  9. Papazafeiropoulos A, Abdullah Z, Kourtessis P, et al, "Coverage Probability of STAR-RIS-Assisted Massive MIMO Systems With Correlation and Phase Errors," IEEE Wireless Communications Letters, vol. 11, no. 8, pp. 1738-1742, Aug. 2022. https://doi.org/10.1109/LWC.2022.3179653
  10. Huang C, Zappone A, Alexandropoulos G C, et al, "Reconfigurable Intelligent Surfaces for Energy Efficiency in Wireless Communication," IEEE Transactions on Wireless Communications, vol. 18, no. 8, pp. 4157-4170, Aug. 2019. https://doi.org/10.1109/TWC.2019.2922609
  11. Zeng S, Zhang H, Di B, et al, "Reconfigurable Intelligent Surface (RIS) Assisted Wireless Coverage Extension: RIS Orientation and Location Optimization," IEEE Communications Letters, vol. 25, no. 1, pp. 269-273, Jan. 2021. https://doi.org/10.1109/LCOMM.2020.3025345
  12. Zhao J, Chen M, Chen M, et al, "Energy Efficient Full-Duplex Communication Systems with Reconfigurable Intelligent Surface," in Proc. of 2020 IEEE 92nd Vehicular Technology Conference (VTC2020-Fall), Victoria, BC, Canada, pp. 1-5, 2020.
  13. Basar E, Di Renzo M, De Rosny J, et al, "Wireless Communications Through Reconfigurable Intelligent Surfaces," IEEE Access, vol. 7, pp. 116753-116773, 2019. https://doi.org/10.1109/ACCESS.2019.2935192
  14. R. Su, L. Dai, J. Tan, M. Hao, and R. MacKenzie, "Capacity enhancement for irregular reconfigurable intelligent surface-aided wireless communications," in Proc. IEEE Global Commun. Conf. (IEEE GLOBECOM'20), Taipei, pp. 1-6, Dec. 2020.
  15. Yang Z, Chen M, Saad W, et al, "Energy-Efficient Wireless Communications With Distributed Reconfigurable Intelligent Surfaces," IEEE Transactions on Wireless Communications, vol. 21, no. 1, pp. 665-679, Jan. 2022. https://doi.org/10.1109/TWC.2021.3098632
  16. R. Su, L. Dai, J. Tan, M. Hao and R. MacKenzie, "Capacity Enhancement for Reconfigurable Intelligent Surface-Aided Wireless Network: from Regular Array to Irregular Array," IEEE Transactions on Vehicular Technology, vol. 72, no. 5, pp. 6392-6403, May 2023. https://doi.org/10.1109/TVT.2023.3236179
  17. A. Taha, M. Alrabeiah and A. Alkhateeb, "Enabling Large Intelligent Surfaces with Compressive Sensing and Deep Learning," IEEE Access, vol. 9, pp. 44304-44321, 2021. https://doi.org/10.1109/ACCESS.2021.3064073
  18. B. Zheng and R. Zhang, "Intelligent Reflecting Surface-Enhanced OFDM: Channel Estimation and Reflection Optimization," IEEE Wireless Communications Letters, vol. 9, no. 4, pp. 518-522, April 2020. https://doi.org/10.1109/LWC.2019.2961357
  19. O. Ozdogan, E. Bjornson and E. G. Larsson, "Intelligent Reflecting Surfaces: Physics, Propagation, and Pathloss Modeling," IEEE Wireless Communications Letters, vol. 9, no. 5, pp. 581-585, May 2020. https://doi.org/10.1109/LWC.2019.2960779