This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Configuring Intelligent Reflecting Surface with Performance Guarantees: Optimal Beamforming

Yaowen Zhang, , Kaiming Shen, ,
Shuyi Ren, , Xin Li, Xin Chen, Zhi-Quan Luo
Manuscript received . This work was supported in part by the National Natural Science Foundation of China (NSFC) under Grant 62001411 and in part by the Huawei Technologies. This article was presented in part at the IEEE Global Communications Conference (GLOBECOM), Madrid, Spain, 2021. Y. Zhang, K. Shen, and S. Ren are with the School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, 518172, China (e-mail: [email protected]; [email protected]; [email protected]). X. Li and X. Chen are with Huawei Technologies (e-mail: [email protected]; [email protected]). Z.-Q. Luo is both with the School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, 518172, China and with Shenzhen Research Institute of Big Data, China (e-mail: [email protected]).
Abstract

This work proposes linear time strategies to optimally configure the phase shifts for the reflective elements of an intelligent reflecting surface (IRS). Specifically, we show that the binary phase beamforming can be optimally solved in linear time to maximize the received signal-to-noise ratio (SNR). For the general KK-ary phase beamforming, we develop a linear time approximation algorithm that guarantees performance within a constant fraction (1+cos(π/K))/2(1+\cos(\pi/K))/2 of the global optimum, e.g., it can attain over 85%85\% of the optimal performance for the quadrature beamforming with K=4K=4. According to the numerical results, the proposed approximation algorithm for discrete IRS beamforming outperforms the existing algorithms significantly in boosting the received SNR.

Index Terms:
Intelligent reflecting surface (IRS), linear time algorithm for discrete beamforming, global optimum, approximation ratio.

I Introduction

Intelligent reflecting surface (IRS) is an emerging 6G technology that uses a large array of low-cost electromagnetic “mirrors”, i.e., reflective elements, to orient the impinging radio waves toward the intended receiver, thereby boosting the spectral efficiency, energy efficiency, and reliability of wireless transmission [1, 2]. From a practical standpoint, the choice for the phase shift induced by each reflective element is normally restricted to a set of discrete values because of hardware constraints. Although IRS has been studied extensively in the literature to date, it remains a challenging task to optimally coordinate discrete phase shifts across reflective elements.

This work proposes efficient strategies for discrete beamforming for IRS in order to maximize the received signal-to-noise ratio (SNR) boost, the computational complexity of which is only linear in the number of reflective elements. For the binary phase beamforming with each phase shift confined to the set {0,π}\{0,\pi\}, we show that the global optimal solution can be computed in linear time; for the general KK-ary phase beamforming case with each phase shift confined to the set {ω,2ω,,Kω}\{\omega,2\omega,\ldots,K\omega\} where ω=2π/K\omega=2\pi/K, we propose a linear time approximation algorithm that guarantees an SNR boost within a constant fraction (1+cos(π/K))/2(1+\cos(\pi/K))/2 of the global optimum, e.g., it can attain over 85%85\% of the optimal SNR boost for the quadrature beamforming with K=4K=4.

Unlike the existing multi-antenna devices, IRS is a passive device that does not perform active signal transmitting. This passive trait of IRS has spurred considerable research interests over the past few years in adapting the traditional methods of active beamforming, covering a variety of system models ranging from point-to-point communication [3, 4, 5, 6] to downlink cellular network [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19], uplink cellular network [20, 21], interference channel [22], and wiretap channel [23]. While the prior studies in this area mostly assume that the phase shift can be chosen arbitrarily for every reflective element, there is a growing trend toward practical beamforming that restricts the choice for phase shift to a set of discrete values [24, 25, 26, 27, 28, 29, 30, 31, 32, 33]. Two common objectives of IRS beamforing as adopted in the literature are to maximize the throughput [8, 3, 11, 10, 20, 18, 19, 15, 13, 23] and to minimize the power consumption under the quality-of-service (QoS) constraints [7, 12, 21, 17, 16], other works accounting for the generalized degree-of-freedom (GDoF) [22], energy efficiency [9], and outage probability [34].

Optimization is a key aspect of the IRS beamforming. In the realm of continuous IRS beamforming, two standard optimization tools—semidefinite relaxation (SDR) [35] and fractional programming (FP) [36, 37]—have been brought to the fore to address the nonvexity of the IRS beamforming problem. For instance, [7, 38, 14, 12, 20, 15, 19] suggest using SDR because of the quadratic programming form of the IRS beamforming problem, while [11, 23, 17, 34, 21, 5, 18, 23, 13] invoke FP to deal with the fractional function of the signal-to-interference-plus-noise ratio (SINR). Moreover, [22] develops a unified Riemannian conjugate gradient algorithm to enable interference alignment in an IRS-assisted system, [10, 16] cope with the nonconvex beamforming problem via successive convex approximation, [6] relies on the minorization-maximization (MM) to address the hardware impairments in the IRS beamforming problem, and [3, 19] suggest a decentralized beamforming scheme based on alternating direction method of multipliers (ADMM).

In comparison to the above continuous beamforming algorithms for IRS, the existing discrete beamforming algorithms lag somewhat in sophistication. Few attempts have been made for the discrete case to date. To achieve the global optimum, [32] uses the exhaustive search and [33] uses the branch-and-bound algorithm, neither of which is computationally scalable. The remaining works [24, 25, 26, 27, 28, 29, 30] typically relax the discrete constraint and then round the relaxed continuous solution to the closest point in the discrete constraint set, namely the closest point projection algorithm. However, this heuristic method can lead to arbitrarily bad performance as shown in this paper. We shall give a linear time beamforming algorithm with provable performance guarantees. The main contributions of the paper are two-fold:

I-1 Binary phase beamforming

We consider maximizing the SNR boost when each phase shift is selected from {0,π}\{0,\pi\}. Although the corresponding problem appears intractable, we show that the global optimum can be achieved in linear time. In addition, we show that the common heuristic of projecting the relaxed continuous solution to the closest point in the discrete constraint set [24, 25, 26, 27, 28, 29, 30, 31] can result in arbitrarily bad performance.

I-2 General KK-ary phase beamforming

We then assume that each phase shift is selected from a fixed set of KK discrete values. The proposed linear time algorithm has provable performance of reaching the global optimum to within a constant fraction 1/2+cos(π/K)/21/2+\cos(\pi/K)/2, thus guaranteeing a tighter approximation ratio than the closest point projection algorithm [24, 25, 26, 27, 28, 29, 30, 31] and the SDR algorithm [7, 38, 14, 12, 20, 15, 19].

The above results are of theoretical significance because they shatter the common beliefs that the discrete IRS beamforming problem must be NP-hard and that the closest point projection is the best strategy for coordinating phase shifts in practice. Furthermore, it is worth mentioning that the proposed algorithm in this work can be extended to blind beamforming without channel estimation. This topic is pursued in the companion paper [39] to the present.

Throughout the paper, we use the boldface lower-case letter to denote a vector, the bold upper-case letter a matrix, and the calligraphy upper-case letter a set. For a matrix 𝑨\bm{A}, 𝑨𝖳\bm{A}^{\mathsf{T}} refers to the transpose, 𝑨𝖧\bm{A}^{\mathsf{H}} the conjugate transpose, and 𝑨1\bm{A}^{-1} the inverse. For a vector 𝒂\bm{a}, 𝒂\|\bm{a}\| refers to the Euclidean norm, 𝒂𝖳\bm{a}^{\mathsf{T}} the transpose, and 𝒂𝖧\bm{a}^{\mathsf{H}} the conjugate transpose. The cardinality of a set 𝒬\mathcal{Q} is denoted as |𝒬||\mathcal{Q}|. The set of real numbers is denoted as \mathbb{R}. The set of complex numbers is denoted as \mathbb{C}. For a complex number uu\in\mathbb{C}, Re{u}\mathrm{Re}\{u\}, Im{u}\mathrm{Im}\{u\}, and Arg(u)\mathrm{Arg}(u) refer to the real part, the imaginary part, and the principal argument of uu, respectively. The rest of the paper is organized as follows. Section II describes the system model and problem formulation. Section III discusses binary phase beamforming. Section IV discusses KK-ary phase beamforming. Section V presents the numerical results. Section VI concludes the paper.

II System Model

Consider a pair of transmitter and receiver, along with an IRS that facilitates the data transmission between them. The IRS consists of NN passive reflective elements. Let hnh_{n}\in\mathbb{C}, n=1,,Nn=1,\ldots,N, be the cascaded channel from the transmitter to the receiver that is induced by the nnth reflective element; let h0h_{0}\in\mathbb{C} be the superposition of all those channels from the transmitter to the receiver that are not related to the IRS, namely the background channel. Each channel can be rewritten in an exponential form as

hn=βnejαn,n=0,,N,h_{n}=\beta_{n}e^{j\alpha_{n}},\;n=0,\ldots,N, (1)

with the magnitude βn(0,1)\beta_{n}\in(0,1) and the phase αn[0,2π)\alpha_{n}\in[0,2\pi).

Denote the IRS beamformer as 𝜽=(θ1,,θN)\bm{\theta}=(\theta_{1},\ldots,\theta_{N}), where each θn[0,2π)\theta_{n}\in[0,2\pi) refers to the phase shift of the nnth reflective element. The choice of each θn\theta_{n} is restricted to the discrete set

ΦK={ω,2ω,,Kω}\Phi_{K}=\big{\{}\omega,2\omega,\ldots,K\omega\big{\}} (2)

with the distance parameter

ω=2πK.\omega=\frac{2\pi}{K}. (3)

Let XX\in\mathbb{C} be the transmit signal with the mean power PP, i.e., 𝔼[|X|2]=P\mathbb{E}[|X|^{2}]=P. The received signal YY\in\mathbb{C} is given by

Y\displaystyle Y =(h0+n=1Nhnejθn)X+Z,\displaystyle=\Bigg{(}h_{0}+\sum^{N}_{n=1}h_{n}e^{j\theta_{n}}\Bigg{)}X+Z, (4)

where an i.i.d. random variable Z𝒞𝒩(0,σ2)Z\sim\mathcal{CN}(0,\sigma^{2}) models the additive thermal noise. The received SNR can be computed as

𝖲𝖭𝖱\displaystyle\mathsf{SNR} =𝔼[|YZ|2]𝔼[|Z|2]\displaystyle=\frac{\mathbb{E}[|Y-Z|^{2}]}{\mathbb{E}[|Z|^{2}]} (5a)
=P|β0ejα0+n=1Nβnej(αn+θn)|2σ2.\displaystyle=\frac{P\Big{|}\beta_{0}e^{j\alpha_{0}}+\sum^{N}_{n=1}\beta_{n}e^{j(\alpha_{n}+\theta_{n})}\Big{|}^{2}}{\sigma^{2}}. (5b)

The baseline SNR without IRS is

𝖲𝖭𝖱0=Pβ02σ2.\mathsf{SNR}_{0}=\frac{P\beta^{2}_{0}}{\sigma^{2}}. (6)

The SNR boost is defined as

f(𝜽)\displaystyle f(\bm{\theta}) =𝖲𝖭𝖱𝖲𝖭𝖱0\displaystyle=\frac{\mathsf{SNR}}{\mathsf{SNR}_{0}} (7a)
=1β02|β0ejα0+n=1Nβnej(αn+θn)|2.\displaystyle=\frac{1}{\beta^{2}_{0}}\left|\beta_{0}e^{j\alpha_{0}}+\sum^{N}_{n=1}\beta_{n}e^{j(\alpha_{n}+\theta_{n})}\right|^{2}. (7b)

We seek the optimal 𝜽\bm{\theta} to maximize the SNR boost:

maximize𝜽\displaystyle\underset{\bm{\theta}}{\text{maximize}} f(𝜽)\displaystyle\quad f(\bm{\theta}) (8a)
subject to θnΦK,n=1,,N.\displaystyle\quad\theta_{n}\in\Phi_{K},\;\forall n=1,\ldots,N. (8b)

The above problem is difficult to tackle directly because of the discrete constraint on 𝜽\bm{\theta}.

III Binary Phase Beamforming with K=2K=2

As KK\rightarrow\infty, the discrete beamforming problem in (8) reduces to the continuous, in which case the SNR boost is maximized when every hnh_{n} is perfectly aligned with h0h_{0}, so the optimal relaxed solution θn\theta^{\infty}_{n} equals the phase difference between h0h_{0} and hnh_{n}:

θn=α0αn.\theta^{\infty}_{n}=\alpha_{0}-\alpha_{n}. (9)

It is then tempting to believe that the problem becomes harder when KK decreases. But it turns out that the problem can be efficiently solved when K=2K=2, as stated below.

Theorem 1

The binary phase beamforming problem in (8) with K=2K=2 can be optimally solved in O(N)O(N) time.

Before proceeding to the proof of the above theorem, we first cast (8) to a quadratic programming form. By letting

𝒙=(1,x1,,xN)𝖳withxn=ejθn,\bm{x}=(1,x_{1},\ldots,x_{N})^{\mathsf{T}}\;\;\text{with}\;\;x_{n}=e^{j\theta_{n}}, (10)

we rewrite the problem in (8) as

maximize𝒙\displaystyle\underset{\bm{x}}{\text{maximize}} 𝒙𝖧𝑪𝒙\displaystyle\quad\bm{x}^{\mathsf{H}}\bm{C}\bm{x} (11a)
subject to xn{ejω,,ejKω},n,\displaystyle\quad x_{n}\in\{e^{j\omega},\ldots,e^{jK\omega}\},\;\forall n, (11b)

where the matrix 𝑪=[Cmn](N+1)×(N+1)\bm{C}=[C_{mn}]_{(N+1)\times(N+1)} is given by

Cmn=hm𝖧hn,m,n=0,,N.C_{mn}=h^{\mathsf{H}}_{m}h_{n},\;\forall m,n=0,\ldots,N. (12)

It can be shown that rank(𝑪)2\mathrm{rank}(\bm{C})\leq 2. We deal with the rank-one case and the rank-two case separately in the following.

III-1 rank(𝑪)=1\mathrm{rank}(\bm{C})=1

In this case, each hnh_{n} is a multiple of h0h_{0}, i.e., the quotient

rn=hnh0.r_{n}=\frac{h_{n}}{h_{0}}\in\mathbb{R}. (13)

When K=2K=2, by substituting (13) back into (11), we obtain

maximize𝒙\displaystyle\underset{\bm{x}}{\text{maximize}} (1+r1x1++rNxN)2\displaystyle\quad(1+r_{1}x_{1}+\ldots+r_{N}x_{N})^{2} (14a)
subject to xn{1,1},n.\displaystyle\quad x_{n}\in\{-1,1\},\;\forall n. (14b)

It can be readily seen that the optimal 𝒙\bm{x} would let every term rnxnr_{n}x_{n} be positive, i.e.,

xn=sgn(rn),n=1,,N,x^{\star}_{n}=\mathrm{sgn}(r_{n}),\;n=1,\ldots,N, (15)

where sgn(z)\mathrm{sgn}(z) equals 11 if z0z\geq 0 and equals 1-1 otherwise. We then recover the optimal phase shift as θn=arccos(xn)\theta_{n}^{\star}=\arccos(x_{n}^{\star}). The above method runs in O(N)O(N) time.

III-2 rank(𝑪)=2\mathrm{rank}(\bm{C})=2

Rewrite each hnh_{n} as a vector 𝒗n2\bm{v}_{n}\in\mathbb{R}^{2}:

𝒗n=[Re{hn}Im{hn}],n=0,,N,\bm{v}_{n}=\left[\begin{matrix}\mathrm{Re}\{h_{n}\}\\ \mathrm{Im}\{h_{n}\}\end{matrix}\right],\;\;n=0,\ldots,N, (16)

and write

𝑽=[𝒗0,,𝒗N].\bm{V}=\big{[}\bm{v}_{0},\ldots,\bm{v}_{N}\big{]}. (17)

With an auxiliary variable 𝒚2\bm{y}\in\mathbb{R}^{2}, we transform the objective function in (11) as follows:

𝒙𝖧𝑪𝒙\displaystyle\bm{x}^{\mathsf{H}}\bm{C}\bm{x} =𝑽𝒙2\displaystyle=\|\bm{V}\bm{x}\|^{2} (18a)
=max𝒚=1𝑽𝒙2𝒚2\displaystyle=\underset{\|\bm{y}\|=1}{\max}\;\|\bm{V}\bm{x}\|^{2}\cdot\|\bm{y}\|^{2} (18b)
=max𝒚=1((𝑽𝒙)𝖳𝒚)2\displaystyle=\underset{\|\bm{y}\|=1}{\max}\;\big{(}(\bm{V}\bm{x})^{\mathsf{T}}\bm{y}\big{)}^{2} (18c)
=max𝒚=1(𝒗0𝖳𝒚+x1𝒗1𝖳𝒚++xN𝒗N𝖳𝒚)2,\displaystyle=\underset{\|\bm{y}\|=1}{\max}\;\big{(}\bm{v}_{0}^{\mathsf{T}}\bm{y}+x_{1}\bm{v}^{\mathsf{T}}_{1}\bm{y}+\ldots+x_{N}\bm{v}^{\mathsf{T}}_{N}\bm{y}\big{)}^{2}, (18d)

where (18c) follows by the Cauchy-Schwarz inequality. Inspection of the new objective function in (18d) shows that the optimal choice of xn{1,1}x_{n}\in\{-1,1\} given 𝒚\bm{y} is to let xn𝒗n𝖳𝒚x_{n}\bm{v}^{\mathsf{T}}_{n}\bm{y} have the same sign as 𝒗0𝖳𝒚\bm{v}^{\mathsf{T}}_{0}\bm{y}, i.e.,

xn=sgn(𝒗0𝖳𝒚𝒗n𝖳𝒚),n=1,,N,x^{\star}_{n}=\mathrm{sgn}(\bm{v}^{\mathsf{T}}_{0}\bm{y}\bm{v}^{\mathsf{T}}_{n}\bm{y}),\;n=1,\ldots,N, (19)

so the joint optimization of (𝒙,𝒚)(\bm{x},\bm{y}) hinges on how to optimize 𝒚\bm{y}.

We treat each 𝒗n\bm{v}_{n}, n=0,1,,Nn=0,1,\ldots,N, as a normal vector and consider its tangent line as shown in Fig. 1. Since multiple 𝒗n\bm{v}_{n}’s may be associated with a common tangent line, the total number of tangent lines MN+1M\leq N+1. Denote an arbitrary tangent line as 1\ell_{1}. Starting from 1\ell_{1}, denote the remaining M1M-1 tangent lines as 2,,M\ell_{2},\ldots,\ell_{M} in counterclockwise order. Let 𝒱m\mathcal{V}_{m} be the set of 𝒗n\bm{v}_{n}’s associated with m\ell_{m}, m=1,,Mm=1,\ldots,M. As illustrated in Fig. 1, these tangent lines partition the unit circle into a total of 2M2M circular segments. Choose a circular segment that first intercepts 1\ell_{1} then intercepts 2\ell_{2} counterclockwise; denote it as 1\mathcal{L}_{1}. The remaining circular segments are denoted as 2,,M\mathcal{L}_{2},\ldots,\mathcal{L}_{M} in counterclockwise order.

Refer to caption

Figure 1: The rank 2 case of binary phase beamforming with N=2N=2. The normal vectors 𝒗0,𝒗1,𝒗2\bm{v}_{0},\bm{v}_{1},\bm{v}_{2} correspond to a tangent line each. The tangent lines, denoted as the dashed lines, partition the unit circle into 6 circular segments.

Notice that on each circular segment xnx^{\star}_{n} in (19) is fixed, so (𝒗0𝖳𝒚+x1𝒗1𝖳𝒚++xN𝒗N𝖳𝒚)\big{(}\bm{v}_{0}^{\mathsf{T}}\bm{y}+x^{\star}_{1}\bm{v}^{\mathsf{T}}_{1}\bm{y}+\ldots+x^{\star}_{N}\bm{v}^{\mathsf{T}}_{N}\bm{y}\big{)} is a piecewise linear function of 𝒚\bm{y}. When 𝒚\bm{y} is restricted to a particular segment m\mathcal{L}_{m}, the optimization problem of 𝒚\bm{y} in (18d) can be rewritten as

maximize𝒚2,𝒚=1\displaystyle\underset{\bm{y}\in\mathbb{R}^{2},\,\|\bm{y}\|=1}{\text{maximize}} 𝒘m𝖳𝒚\displaystyle\quad\bm{w}^{\mathsf{T}}_{m}\bm{y} (20)

with the coefficient vector

𝒘m=n=0Nsgn(𝒗n𝖳𝒚)𝒗n,\bm{w}_{m}=\sum^{N}_{n=0}\mathrm{sgn}(\bm{v}^{\mathsf{T}}_{n}\bm{y}^{\prime})\bm{v}_{n}, (21)

where 𝒚2\bm{y}^{\prime}\in\mathbb{R}^{2} is an arbitrary vector located in the interior of m\mathcal{L}_{m}. A naive idea is to compute every 𝒘m\bm{w}_{m} by using (21), requiring O(N)O(N) time; this would result in O(N2)O(N^{2}) time in total for all m=1,,Mm=1,\ldots,M. But we show that the time can be reduced to O(N)O(N) via iterative updating. When updating xnx^{\star}_{n} from one circular segment m\mathcal{L}_{m} to the next circular segment m+1\mathcal{L}_{m+1}, we just need to change the values of those xnx^{\star}_{n} variables related to the tangent line m+1\ell_{m+1} that separates the two circular segments, i.e., {xn:n𝒱m+1}\{x^{\star}_{n}:\forall n\in\mathcal{V}_{m+1}\}, so 𝒘m+1\bm{w}_{m+1} can be computed iteratively based on 𝒘m\bm{w}_{m} as

𝒘m+1=𝒘mn𝒱m+12sgn(𝒗n𝖳𝒚)𝒗n,\bm{w}_{m+1}=\bm{w}_{m}-\sum_{n\in\mathcal{V}_{m+1}}2\cdot\mathrm{sgn}(\bm{v}^{\mathsf{T}}_{n}\bm{y}^{\prime})\bm{v}_{n}, (22)

with an arbitrary 𝒚2\bm{y}^{\prime}\in\mathbb{R}^{2} located in the interior of m\mathcal{L}_{m}. Thus, while 𝒘1\bm{w}_{1} is still computed as in (21) in O(N)O(N) time, the other coefficient vectors 𝒘2,,𝒘M\bm{w}_{2},\ldots,\bm{w}_{M} can be obtained sequentially according to (22) in O(m=2M|𝒱m|)=O(N)O\big{(}\sum^{M}_{m=2}|\mathcal{V}_{m}|\big{)}=O(N) time in total.

Moreover, the optimal 𝒚\bm{y} in problem (20) for each mm can be obtained in O(1)O(1) time as

𝒚m=argmin𝒚2,𝒚=1𝒚𝒘m𝒘m,m=1,,M,\bm{y}^{\star}_{m}=\arg\min_{\bm{y}\in\mathbb{R}^{2},\,\|\bm{y}\|=1}\left\|\bm{y}-\frac{\bm{w}_{m}}{\|\bm{w}_{m}\|}\right\|,\;m=1,\ldots,M, (23)

i.e., 𝒚m=𝒘m/𝒘m\bm{y}^{\star}_{m}=\bm{w}_{m}/\|\bm{w}_{m}\| if 𝒘m/𝒘mm\bm{w}_{m}/\|\bm{w}_{m}\|\in\mathcal{L}_{m}; otherwise 𝒚m\bm{y}^{\star}_{m} is located at one of the endpoints of the circular segment m\mathcal{L}_{m}. As a result, it requires O(N)O(N) time overall to solve a sequence of problems in (20) for m=1,,Mm=1,\ldots,M.

Finally, the actual optimal 𝒚\bm{y}^{\star} is given by

𝒚=max{𝒚1,,𝒚M}.\bm{y}^{\star}=\max\{\bm{y}^{\star}_{1},\ldots,\bm{y}^{\star}_{M}\}. (24)

We then compute the corresponding xnx^{\star}_{n} as in (19) and recover the optimal phase shift by θn=arccos(xn)\theta^{\star}_{n}=\arccos(x^{\star}_{n}). The entire procedure requires O(N)O(N) time in total. The proof of Theorem 1 is then completed. We summarize the above steps in Algorithm 1.

We remark that the binary phase beamforming case is fairly special in that the variable xnx_{n} in (10), either 1-1 or 11, is real-valued. When K>2K>2, xnx_{n} is complex-valued in general and thus the above linear time algorithm no longer applies.

In the existing literature [24, 25, 26, 27, 28, 29, 30], a common way of discrete beamforming is to round the relaxed continuous solution θn\theta^{\infty}_{n} in (9) to the closest point in the discrete set ΦK\Phi_{K}, i.e.,

θnCPP=argminθnΦK|θnθn|,{\theta}^{\text{CPP}}_{n}=\arg\min_{\theta_{n}\in\Phi_{K}}\big{|}\theta_{n}-\theta^{\infty}_{n}\big{|}, (25)

referred to as the Closest Point Projection (CPP) algorithm. The following example however shows that the above algorithm can result in arbitrarily bad performance when K=2K=2.

Example 1

Consider the binary phase beamforming with K=2K=2. Assume that all the reflected channels have equal magnitudes, and that half of the reflected channels have the phase αn=α0ϵ+π/2\alpha_{n}=\alpha_{0}-\epsilon+\pi/2 while the other half have the phase αn=α0+ϵπ/2\alpha_{n}=\alpha_{0}+\epsilon-\pi/2. As ϵ0\epsilon\rightarrow 0, the closest point projection algorithm leads to f(𝛉CPP)1f(\bm{\theta}^{\text{CPP}})\rightarrow 1, i.e., IRS does not bring any improvements.

Algorithm 1 Proposed Binary Phase Beamforming Method
1:input: 𝒉0,𝒉1,,𝒉N\bm{h}_{0},\bm{h}_{1},\ldots,\bm{h}_{N}
2:for n=0,1,,Nn=0,1,\ldots,N do
3:     let 𝒗n=[Re{hn},Re{hn}]𝖳\bm{v}_{n}=\big{[}\mathrm{Re}\{h_{n}\},\mathrm{Re}\{h_{n}\}\big{]}^{\mathsf{T}}
4:     treat 𝒗n\bm{v}_{n} as the normal vector and draw its tangent line
5:end for
6:obtain MNM\leq N distinct tangent lines 1,,M\ell_{1},\ldots,\ell_{M}; let 𝒱m\mathcal{V}_{m} be the set of 𝒗n\bm{v}_{n}’s associated with m\ell_{m}, m=0,1,,Mm=0,1,\ldots,M
7:the unit circle is partitioned into MM circular segments 1,,M\mathcal{L}_{1},\ldots,\mathcal{L}_{M} as shown in Fig. 1
8:compute 𝒘1\bm{w}_{1} as in (21)
9:compute 𝒚1\bm{y}^{\star}_{1} as in (23)
10:for m=2,,Mm=2,\ldots,M do
11:     compute 𝒘m\bm{w}_{m} as in (22)
12:     compute 𝒚m\bm{y}^{\star}_{m} as in (23)
13:end for
14:𝒚=max{𝒚1,,𝒚M}\bm{y}^{\star}=\max\{\bm{y}^{\star}_{1},\ldots,\bm{y}^{\star}_{M}\}
15:for n=1,,Nn=1,\ldots,N do
16:     xn=sgn(𝒗0𝖳𝒚𝒗n𝖳𝒚)x^{\star}_{n}=\mathrm{sgn}(\bm{v}^{\mathsf{T}}_{0}\bm{y}^{\star}\bm{v}^{\mathsf{T}}_{n}\bm{y}^{\star})
17:     θn=arccos(xn)\theta^{\star}_{n}=\arccos(x^{\star}_{n})
18:end for
19:output: 𝜽=(θ1,,θN)\bm{\theta}^{\star}=(\theta^{\star}_{1},\ldots,\theta^{\star}_{N})

Refer to caption

Figure 2: The four sectors 𝒮1\mathcal{S}_{1} to 𝒮4\mathcal{S}_{4}. In the proof of Theorem 2, we first rotate μ1\mu_{1} clockwise by an angle of ω/2\omega/2, then combine it with μ3\mu_{3} to obtain μ13\mu_{13}.

IV General KK-Ary Phase Beamforming

This section considers the general KK-ary phase beamforming in (8). We begin with a sectorization scheme in order to reduce the search space of the optimal beamformer. This new idea then yields a linear time approximation algorithm that guarantees performance within a constant fraction of the global optimum. It is further shown that the proposed algorithm has a tighter approximation ratio than the existing algorithms.

IV-A Optimal KK-Ary Phase Beamforming

The following sectorization scheme is the building block of our approach to the KK-ary phase beamforming. Consider four sectors around h0h_{0} on the complex plane as illustrated in Fig. 2:

𝒮i={u:α0+(2i)ω2Arg(u)α0+(3i)ω2},i=1,2,3,4.\mathcal{S}_{i}=\bigg{\{}u\in\mathbb{C}:\alpha_{0}+\frac{(2-i)\omega}{2}\leq\mathrm{Arg}(u)\leq\alpha_{0}+\frac{(3-i)\omega}{2}\bigg{\}},\\ \forall i=1,2,3,4. (26)

The overall channel superposition under the optimal beamformer 𝜽\bm{\theta}^{\star} is denoted as

g=h0+n=1Nhnejθn.g^{\star}=h_{0}+\sum^{N}_{n=1}h_{n}e^{j\theta_{n}^{\star}}. (27)

We remark that 𝜽\bm{\theta}^{\star} must satisfy

|Arg(h0)Arg(n=1Nhnejθn)|ω2;\left|\mathrm{Arg}(h_{0})-\mathrm{Arg}\left(\sum^{N}_{n=1}h_{n}e^{j\theta^{\star}_{n}}\right)\right|\leq\frac{\omega}{2}; (28)

otherwise we could find an integer kk such that

|Arg(h0)Arg(n=1Nhnej(θn+kω))|ω2,\left|\mathrm{Arg}(h_{0})-\mathrm{Arg}\left(\sum^{N}_{n=1}h_{n}e^{j(\theta^{\star}_{n}+k\omega)}\right)\right|\leq\frac{\omega}{2}, (29)

and thus further increase the SNR boost. The bound in (28) implies that g(𝒮2𝒮3)g^{\star}\in(\mathcal{S}_{2}\cup\mathcal{S}_{3}), so the possible values of 𝜽\bm{\theta}^{\star} are limited, as stated in the following proposition.

Proposition 1

For the KK-ary phase beamforming problem in (8), the optimal solution 𝛉\bm{\theta}^{\star} is contained in either 𝒢(𝒮1𝒮2𝒮3)\mathcal{G}(\mathcal{S}_{1}\cup\mathcal{S}_{2}\cup\mathcal{S}_{3}) or 𝒢(𝒮2𝒮3𝒮4)\mathcal{G}(\mathcal{S}_{2}\cup\mathcal{S}_{3}\cup\mathcal{S}_{4}), where the function 𝒢()\mathcal{G}(\cdot) is given by

𝒢(𝒳)={𝜽ΦKN|hnejθn𝒳,n=1,,N}\mathcal{G}(\mathcal{X})=\left\{\bm{\theta}\in\Phi^{N}_{K}\,\big{|}\,h_{n}e^{j\theta_{n}}\in\mathcal{X},\;\forall n=1,\ldots,N\right\} (30)

with the subset 𝒳\mathcal{X}\subseteq\mathbb{C}.

Proof:

If gg^{\star} is already known, then it is optimal to rotate every hnh_{n} to the closest possible position to gg^{\star}. As a result, every hnejθnh_{n}e^{j\theta^{\star}_{n}} must be in 𝒮1𝒮2𝒮3\mathcal{S}_{1}\cup\mathcal{S}_{2}\cup\mathcal{S}_{3} if g𝒮2g^{\star}\in\mathcal{S}_{2}, and in 𝒮2𝒮3𝒮4\mathcal{S}_{2}\cup\mathcal{S}_{3}\cup\mathcal{S}_{4} if g𝒮3g^{\star}\in\mathcal{S}_{3}. ∎

In light of the above proposition, it suffices to search for 𝜽\bm{\theta}^{\star} in 𝒢(𝒮1𝒮2𝒮3)𝒢(𝒮2𝒮3𝒮4)\mathcal{G}(\mathcal{S}_{1}\cup\mathcal{S}_{2}\cup\mathcal{S}_{3})\cup\mathcal{G}(\mathcal{S}_{2}\cup\mathcal{S}_{3}\cup\mathcal{S}_{4}), i.e.,

𝜽=argmax𝜽𝒟f(𝜽),{\bm{\theta}}^{\star}=\arg\max_{\bm{\theta}\in\mathcal{D}}\,f(\bm{\theta}), (31)

where

𝒟=𝒢(𝒮1𝒮2𝒮3)𝒢(𝒮2𝒮3𝒮4).\mathcal{D}=\mathcal{G}(\mathcal{S}_{1}\cup\mathcal{S}_{2}\cup\mathcal{S}_{3})\cup\mathcal{G}(\mathcal{S}_{2}\cup\mathcal{S}_{3}\cup\mathcal{S}_{4}). (32)

Note that the search space size |𝒟|=2N|\mathcal{D}|=2^{N} is much smaller than the full solution space size |ΦKN|=KN|\Phi^{N}_{K}|=K^{N} for large KK. But the running time is still exponential. We aim at a more efficient algorithm in the next subsection.

IV-B Proposed Approximation Algorithm

Algorithm 2 Proposed KK-Ary Phase Beamforming Method
1:input: 𝒉0,𝒉1,,𝒉N\bm{h}_{0},\bm{h}_{1},\ldots,\bm{h}_{N}
2:𝒟=\mathcal{D}^{\prime}=\emptyset
3:compute the sectors 𝒮1,𝒮2,𝒮3,𝒮4\mathcal{S}_{1},\mathcal{S}_{2},\mathcal{S}_{3},\mathcal{S}_{4} as in (26)
4:for i=1,2,3i=1,2,3 do
5:     find 𝜽\bm{\theta} that lets every hnejθn𝒮i𝒮i+1h_{n}e^{j\theta_{n}}\in\mathcal{S}_{i}\cup\mathcal{S}_{i+1}
6:     𝒟𝒟{𝜽}\mathcal{D}^{\prime}\rightarrow\mathcal{D}^{\prime}\cup\{\bm{\theta}\}
7:end for
8:𝜽APX=argmax𝜽𝒟f(𝜽){\bm{\theta}^{\text{APX}}}=\arg\max_{\bm{\theta}\in{\mathcal{D}^{\prime}}}\,f(\bm{\theta})
9:output: 𝜽APX\bm{\theta}^{\text{APX}}

The premise behind the algorithm in (31) and (32) is that the optimally phase-shifted channels {h1ejθ1,,hNejθN}\{h_{1}e^{j\theta_{1}^{\star}},\ldots,h_{N}e^{j\theta_{N}^{\star}}\} can spread up to three consecutive sectors, i.e., 𝒮1𝒮2𝒮3\mathcal{S}_{1}\cup\mathcal{S}_{2}\cup\mathcal{S}_{3} or 𝒮2𝒮3𝒮4\mathcal{S}_{2}\cup\mathcal{S}_{3}\cup\mathcal{S}_{4}. It guarantees the global optimality but results in the search space 𝒟\mathcal{D} in (32) being exponentially large.

To achieve a trade-off between complexity and performance, we propose letting the phase-shifted channels {h1ejθ1,,hNejθN}\{h_{1}e^{j\theta_{1}},\ldots,h_{N}e^{j\theta_{N}}\} spread at most two consecutive sectors, i.e., 𝒮1𝒮2\mathcal{S}_{1}\cup\mathcal{S}_{2} or 𝒮2𝒮3\mathcal{S}_{2}\cup\mathcal{S}_{3} or 𝒮3𝒮4\mathcal{S}_{3}\cup\mathcal{S}_{4}. We now decide the IRS beamformer as

𝜽APX=argmax𝜽𝒟f(𝜽){\bm{\theta}^{\text{APX}}}=\arg\max_{\bm{\theta}\in{\mathcal{D}^{\prime}}}\,f(\bm{\theta}) (33)

with a new search space

𝒟=𝒢(𝒮1𝒮2)𝒢(𝒮2𝒮3)𝒢(𝒮3𝒮4).{\mathcal{D}^{\prime}}=\mathcal{G}(\mathcal{S}_{1}\cup\mathcal{S}_{2})\cup\mathcal{G}(\mathcal{S}_{2}\cup\mathcal{S}_{3})\cup\mathcal{G}(\mathcal{S}_{3}\cup\mathcal{S}_{4}). (34)

Algorithm 2 summarizes the details.

Notice that |𝒟|3|\mathcal{D}^{\prime}|\leq 3, so the search in (33) requires O(1)O(1) time. Moreover, the new search space 𝒟\mathcal{D}^{\prime} can be obtained in O(N)O(N) time. Thus, the total running time is O(N)O(N). Most importantly, the resulting SNR boost f(𝜽APX)f(\bm{\theta}^{\text{APX}}) is within a constant fraction of the highest possible boost.

Theorem 2

For the KK-ary phase beamforming problem in (8), the IRS beamformer 𝛉APX\bm{\theta}^{\text{APX}} in (33) by the approximation algorithm and the optimal IRS beamformer 𝛉\bm{\theta}^{\star} satisfy

1+cos(π/K)2f(𝜽)f(𝜽APX)f(𝜽).\frac{1+\cos(\pi/K)}{2}f(\bm{\theta}^{\star})\leq f(\bm{\theta}^{\text{APX}})\leq f(\bm{\theta}^{\star}). (35)
Proof:

Without loss of generality, assume that the optimal channel superposition gg^{\star} in (27) is located in 𝒮2\mathcal{S}_{2}. We partition gg^{\star} into three components:

g=μ1+μ2+μ3,g^{\star}=\mu_{1}+\mu_{2}+\mu_{3}, (36)

where μi\mu_{i}\in\mathbb{C} represents the superposition of those channels in {h0,h1ejθ1,,hNejθN}\{h_{0},h_{1}e^{j\theta^{\star}_{1}},\ldots,h_{N}e^{j\theta^{\star}_{N}}\} that are located in 𝒮i\mathcal{S}_{i}. Defining an auxiliary variable (see Fig. 2 for illustration)

μ13=μ1ejω+μ3,\mu_{13}=\mu_{1}e^{-j\omega}+\mu_{3}, (37)

we can bound the value of f(𝜽)f(\bm{\theta}^{\star}) from above as

f(𝜽)\displaystyle f(\bm{\theta}^{\star}) =1β02|μ1+μ2+μ3|2\displaystyle=\frac{1}{\beta^{2}_{0}}\big{|}\mu_{1}+\mu_{2}+\mu_{3}\big{|}^{2}
1β02(|μ1+μ3|+|μ2|)2\displaystyle\leq\frac{1}{\beta^{2}_{0}}\big{(}\big{|}\mu_{1}+\mu_{3}\big{|}+\big{|}\mu_{2}\big{|}\big{)}^{2}
1β02(|μ13|+|μ2|)2,\displaystyle\leq\frac{1}{\beta^{2}_{0}}\big{(}\big{|}\mu_{13}\big{|}+\big{|}\mu_{2}\big{|}\big{)}^{2}, (38)

where the last inequality follows since the rotation ejωe^{-j\omega} of μ1\mu_{1} in (37) reduces the angle between μ1\mu_{1} and μ3\mu_{3}.

Because 𝜽APX\bm{\theta}^{\text{APX}} is located in either 𝒢(𝒮1𝒮2)\mathcal{G}(\mathcal{S}_{1}\cup\mathcal{S}_{2}) or 𝒢(𝒮2𝒮3)\mathcal{G}(\mathcal{S}_{2}\cup\mathcal{S}_{3}) when g𝒮2g^{\star}\in\mathcal{S}_{2}, we can lower bound the value of f(𝜽APX)f(\bm{\theta}^{\text{APX}}) as

f(𝜽APX)\displaystyle f(\bm{\theta}^{\text{APX}})
=1β02max{|μ1ejω+μ2+μ3|2,|μ1+μ2+μ3ejω|2}\displaystyle=\frac{1}{\beta^{2}_{0}}\cdot\max\big{\{}\big{|}\mu_{1}e^{-j\omega}+\mu_{2}+\mu_{3}\big{|}^{2},\big{|}\mu_{1}+\mu_{2}+\mu_{3}e^{j\omega}\big{|}^{2}\big{\}}
=1β02max{|μ13+μ2|2,|μ13ejω+μ2|2}\displaystyle=\frac{1}{\beta^{2}_{0}}\cdot\max\big{\{}\big{|}\mu_{13}+\mu_{2}\big{|}^{2},\big{|}\mu_{13}e^{j\omega}+\mu_{2}\big{|}^{2}\big{\}}
1β02min|μ2|=|μ2|max{|μ13+μ2|2,|μ13ejω+μ2|2}\displaystyle\geq\frac{1}{\beta^{2}_{0}}\cdot\min_{|\mu^{\prime}_{2}|=|\mu_{2}|}\max\big{\{}\big{|}\mu_{13}+\mu^{\prime}_{2}\big{|}^{2},\big{|}\mu_{13}e^{j\omega}+\mu^{\prime}_{2}\big{|}^{2}\big{\}}
=1β02(|μ2|2+|μ13|2+2|μ2μ13|cos(ω/2)),\displaystyle=\frac{1}{\beta^{2}_{0}}\cdot\left(|\mu_{2}|^{2}+|\mu_{13}|^{2}+2\big{|}\mu_{2}\mu_{13}\big{|}\cos(\omega/2)\right), (39)

where the last equality follows by fact that μ2=(|μ2|/|μ13|)μ13ejω/2\mu^{\prime}_{2}=(|\mu_{2}|/|\mu_{13}|)\mu_{13}e^{j\omega/2} is the solution to the above min-max problem. With λ=|μ13|/|μ2|\lambda=|\mu_{13}|/|\mu_{2}|, we establish the following upper bound by combining (IV-B) and (IV-B):

f(𝜽)f(𝜽APX)\displaystyle\frac{f({\bm{\theta}}^{\star})}{f(\bm{\theta}^{\text{APX}})} (λ+1)2λ2+1+2λcos(ω/2)()21+cos(ω/2),\displaystyle\leq\frac{(\lambda+1)^{2}}{\lambda^{2}+1+2\lambda\cos(\omega/2)}\overset{(*)}{\leq}\frac{2}{1+\cos(\omega/2)}, (40)

with equality in ()(*) if and only if λ=1\lambda=1. Finally, plugging ω=2π/K\omega=2\pi/K in (40) completes the proof. ∎

Remark 1

The proposed algorithm can be extended to blind beamforming without any channel information as discussed in [39]. The work in [39] further proves that a polynomial number of random samples suffice to guarantee the optimal quadratic SNR boost.

IV-C Other KK-Ary Phase Beamforming Algorithms

The proposed approximation algorithm is now compared with the existing methods for the KK-ary phase beamforming. We begin with the closest point projection algorithm.

Proposition 2

For the KK-ary phase beamforming problem in (8), the IRS beamformer 𝛉CPP\bm{\theta}^{\text{CPP}} in (25) by the closest point projection algorithm and the optimal beamformer 𝛉\bm{\theta}^{\star} satisfy

cos2(π/K)f(𝜽)f(𝜽CPP)f(𝜽).\cos^{2}(\pi/K)f(\bm{\theta}^{\star})\leq f(\bm{\theta}^{\text{CPP}})\leq f(\bm{\theta}^{\star}). (41)
Proof:

Clearly, we can upper bound the optimal SNR boost as f(𝜽)(1/β02)(n=0Nβn)2f(\bm{\theta}^{\star})\leq(1/\beta^{2}_{0})\big{(}\sum^{N}_{n=0}\beta_{n}\big{)}^{2} by assuming that all the channels can be aligned exactly. We also have

f(𝜽CPP)\displaystyle f(\bm{\theta}^{\text{CPP}}) =1β02|β0ejα0+n=1Nβnej(θnCPP+αn)|2\displaystyle=\frac{1}{\beta^{2}_{0}}\cdot\left|\beta_{0}e^{j\alpha_{0}}+\sum^{N}_{n=1}\beta_{n}e^{j\left(\theta^{\text{CPP}}_{n}+\alpha_{n}\right)}\right|^{2} (42a)
=1β02|β0+n=1Nβnej(θnCPPθn)|2\displaystyle=\frac{1}{\beta^{2}_{0}}\cdot\left|\beta_{0}+\sum^{N}_{n=1}\beta_{n}e^{j\left(\theta^{\text{CPP}}_{n}-\theta^{\infty}_{n}\right)}\right|^{2} (42b)
1β02|β0+n=1Nβncos(θnCPPθn)|2\displaystyle\geq\frac{1}{\beta^{2}_{0}}\cdot\left|\beta_{0}+\sum^{N}_{n=1}\beta_{n}\cos\left(\theta^{\text{CPP}}_{n}-\theta^{\infty}_{n}\right)\right|^{2} (42c)
1β02|β0+n=1Nβncosω2|2\displaystyle\geq\frac{1}{\beta^{2}_{0}}\cdot\Bigg{|}\beta_{0}+\sum^{N}_{n=1}\beta_{n}\cos\frac{\omega}{2}\Bigg{|}^{2} (42d)
cos2(ω/2)1β02|n=0Nβn|2\displaystyle\geq\cos^{2}(\omega/2)\frac{1}{\beta^{2}_{0}}\cdot\Bigg{|}\sum^{N}_{n=0}\beta_{n}\Bigg{|}^{2} (42e)
cos2(ω/2)f(𝜽),\displaystyle\geq\cos^{2}(\omega/2)f(\bm{\theta}^{\star}), (42f)

where (42d) follows since |θnCPPθn|ω/2|\theta^{\text{CPP}}_{n}-\theta^{\infty}_{n}|\leq\omega/2 by the closest point projection. The proof is then completed. ∎

Furthermore, we remark that the lower bound in (41) can be tight as shown in Fig. 3.

Refer to caption

Figure 3: The worst-case scenario of the closest point projection algorithm. If μ1=0\mu_{1}=0 and |μ2|=|μ3|>0|\mu_{2}|=|\mu_{3}|>0, then f(𝜽CPP)=(4/β02)|μ2|2cos2(π/K)f(\bm{\theta}^{\text{CPP}})=(4/\beta^{2}_{0})|\mu_{2}|^{2}\cos^{2}(\pi/K) and f(𝜽)=(4/β02)|μ2|2f(\bm{\theta}^{\star})=(4/\beta^{2}_{0})|\mu_{2}|^{2} as γ=β0=0\gamma=\beta_{0}=0, so f(𝜽CPP)=cos2(π/K)f(𝜽)f(\bm{\theta}^{\text{CPP}})=\cos^{2}(\pi/K)f(\bm{\theta}^{\star}).

While the previous works [7, 38, 14, 12, 20, 15, 19] have used SDR for continuous beamforming, we show that SDR works for discrete beamforming as well. Toward this end, we first rewrite the problem as a complex KK-ary quadratic problem in (11), and then the standard SDR method applies. The next proposition is a direct result from the existing studies on SDR [40, 35].

Proposition 3

The KK-ary phase beamforming problem in (8) can be rewritten as a complex KK-ary quadratic problem in (11); the new problem can be approximately solved by the standard SDR method [40, 35]. The resulting IRS beamformer 𝛉SDR\bm{\theta}^{\text{SDR}} and the optimal IRS beamformer 𝛉\bm{\theta}^{\star} satisfy

(Ksin(π/K))24πf(𝜽)f(𝜽SDR)f(𝜽).\frac{(K\sin(\pi/K))^{2}}{4\pi}f(\bm{\theta}^{\star})\leq f(\bm{\theta}^{\text{SDR}})\leq f(\bm{\theta}^{\star}). (43)

We conclude this section by comparing the approximation ratios of the various algorithms in Fig. 4. It shows that the guaranteed performance of the proposed approximation algorithm is much closer to the global optimum.

Refer to caption

Figure 4: Algorithm 2 (APX) vs. closest point projection (CPP) algorithm vs. semidefine relaxation (SDR) algorithm in terms of the approximation ratio.

V Simulation Results

Refer to caption

Figure 5: The cumulative distribution of the SNR boost when K=2K=2.

Refer to caption

Figure 6: The cumulative distribution of the SNR boost when K=4K=4.

Refer to caption

Figure 7: The number of reflective elements NN vs. the 1st percentile SNR boost.

Refer to caption

Figure 8: The background-to-reflection ratio β0/βn\beta_{0}/\beta_{n} vs. the average SNR boost.

Refer to caption

Figure 9: The channel ratio β0/βn\beta_{0}/\beta_{n} vs. the 1st percentile SNR boost.

In this section we validate the performance of the proposed algorithms in simulations. The channel model follows the previous works [7, 29, 41]. The background channel h0h_{0} is given by

h0=10𝖯𝖫0/20ζ0,h_{0}=10^{-\mathsf{PL}_{0}/20}\cdot\zeta_{0}, (44)

where the transmitter-to-receiver pathloss 𝖯𝖫0\mathsf{PL}_{0} (in dB) is computed as 𝖯𝖫0=32.6+36.7log10(d0)\mathsf{PL}_{0}=32.6+36.7\log_{10}(d_{0}), with d0d_{0} in meters denoting the distance between the transmitter and the receiver, while the Rayleigh fading component ζ0\zeta_{0} is drawn from the Gaussian distribution 𝒞𝒩(0,1)\mathcal{CN}(0,1). The reflected channel hnh_{n} is given by

hn=10(𝖯𝖫1+𝖯𝖫2)/20ζn,n=1,,N,h_{n}=10^{-(\mathsf{PL}_{1}+\mathsf{PL}_{2})/20}\cdot\zeta_{n},n=1,\ldots,N, (45)

where 𝖯𝖫1\mathsf{PL}_{1} and 𝖯𝖫2\mathsf{PL}_{2} are both based on the pathloss model 𝖯𝖫0=30+22log10(d)\mathsf{PL}_{0}=30+22\log_{10}(d), with dd in meters respectively denoting the transmitter-to-IRS distance and the IRS-to-receiver distance, while the Rayleigh fading component ζn\zeta_{n} is drawn from the Gaussian distribution 𝒞𝒩(0,1)\mathcal{CN}(0,1) independently across n=1,,Nn=1,\ldots,N. We use the following parameters unless otherwise stated. The transmit power level PP equals 3030 dBm. The background noise power level σ2\sigma^{2} equals 90-90 dBm. The locations of the transmitter, IRS, and receiver are denoted by the 3-dimensional coordinate vectors (50,200,20)(50,-200,20), (2,1,0)(-2,-1,0), and (0,0,0)(0,0,0) in meters, respectively. The number of reflective elements NN equals 200200. The channels are estimated via an ON-OFF strategy in [42].

Fig. 5 shows the cumulative distribution of the SNR boost for the binary phase beamforming case with K=2K=2. When N=100N=100, it can be seen that the proposed approximation algorithm “APX” outperforms the closest point projection algorithm “CPP”, especially in the low SNR boost regime. For instance, the 5th percentile SNR boost by APX is more than 2 dB higher than that by CPP. This performance gap remains when the number of reflective elements raises to 200. Moreover, we consider the quadrature beamforming case with K=4K=4 in Fig. 6. Observe that APX and CPP now have similar performance. Thus, APX is more suited for the binary phase beamforming for IRS.

We further look into the low SNR boost regime by comparing the 1st percentile SNR boosts (i.e., the 1% lowest SNR boost across the random channel realizations) in Fig. 7. Different values of NN are considered here. The figure shows that APX outperforms CPP significantly in the low SNR boost regime when K=2K=2. For instance, when N=200N=200, APX improves upon CPP by about 5 dB. Observe that the gap between APX and CPP decreases with NN. When KK is increased to 44, the advantage of APX over CPP becomes marginal.

Fig. 8 shows how the reflected channel strength impacts the SNR boost. We put the IRS sequentially at the following 6 positions: (1,1,0)(-1,-1,0), (2,1,0)(-2,-1,0), (2.5,1,0)(-2.5,-1,0), (3,1,0)(-3,-1,0), (3.5,1,0)(-3.5,-1,0), (4,1,0)(-4,-1,0), so the reflected channel magnitude βn\beta_{n} varies among these IRS locations, while the background channel magnitude β0\beta_{0} is fixed. We evaluate the channel ratio β0/βn\beta_{0}/\beta_{n} for each IRS location. As shown in Fig. 8, the gap between APX and CPP in terms of the average SNR boost increases with β0/βn\beta_{0}/\beta_{n}, in both the binary phase beamforming case and the quadrature beamforming case. Furthermore, Fig. 9 shows the low SNR boost performance with respect to the different values of β0/βn\beta_{0}/\beta_{n}. The figure also suggests that the gap between APX and CPP becomes larger when β0/βn\beta_{0}/\beta_{n} is increased. Summarizing the above observations, we conclude that APX is more suited for the scenario where the reflected channels are relatively weak in comparison to the background channel.

Refer to caption

Figure 10: The cumulative distributions of the SNR boost under the different noise powers when K=2K=2.

The noise power σ2\sigma^{2} is fixed at 90-90 dBm in all the above simulations. We now raise σ2\sigma^{2} to 50-50 dBm and plot the resulting cumulative distributions of SNR boosts in Fig. 10. Due to more severe noise, the channel estimation becomes less accurate, and thus can spoil the IRS beamforming algorithms. Observe from Fig. 10 that APX is less affected by the increased noise as compared to CPP. Actually, APX is far superior to CPP when σ2=50\sigma^{2}=-50 dBm, whereas the two algorithms are close when σ2=90\sigma^{2}=-90 dBm. The tolerance to channel estimation error enables APX to perform much better than CPP in the presence of strong interference and background noise.

VI Conclusion

Despite the practical discrete phase value constraints for the IRS beamforming, this work shows that the optimal SNR boost can be reached by the proposed algorithm with a linear time complexity in the number of reflective elements, if the phase value is binary. For the general KK-ary phase beamforming, we propose a linear time algorithm with a provable approximation accuracy, whereas the existing closest point projection algorithm can result in arbitrarily bad performance. Simulation results demonstrate the superior performance of the proposed approximation algorithm over the existing closest point projection algorithm in enhancing the received SNR, especially when channel estimation is error-prone because of interference and noise.

References

  • [1] Q. Wu and R. Zhang, “Towards smart and reconfigurable environment: Intelligent reflecting surface aided wireless network,” IEEE Commun. Mag., vol. 58, no. 1, pp. 106–112, Jan. 2020.
  • [2] C. Huang, A. Zappone, G. C. Alexandropoulos, M. Debbah, and C. Yuen, “Reconfigurable intelligent surface for energy efficiency in wireless communication,” IEEE Trans. Wireless Commun., vol. 18, no. 8, pp. 4157–4170, Aug. 2019.
  • [3] B. Ning, Z. Chen, W. Chen, and J. Fang, “Beamforming optimization for intelligent reflecting surface assisted MIMO: A sum-path-gain maximization approach,” IEEE Wireless Commun. Lett., vol. 9, no. 7, pp. 1105–1109, Jul. 2020.
  • [4] Y. Han, S. Zhang, L. Duan, and R. Zhang, “Cooperative double-RIS aided communication: Beamforming design and power scaling,” IEEE Wireless Commun. Lett., vol. 9, no. 8, pp. 1206–1210, Aug. 2020.
  • [5] Z. Zhang, L. Dai, X. Chen, C. Liu, F. Yang, R. Schober, and H. V. Poor, “Active RIS vs. passive RIS: Which will prevail in 6G?” 2021, [Online]. Available: https://arxiv.org/abs/2103.15154.
  • [6] H. Shen, W. Xu, S. Gong, C. Zhao, and D. W. K. Ng, “Beamforming optimization for IRS-aided communications with transceiver hardware impairments,” IEEE Trans. Commun., vol. 69, no. 2, pp. 1214–1227, Feb. 2021.
  • [7] Q. Wu and R. Zhang, “Intelligent reflecting surface enhanced wireless network via joint active and passive beamforming,” IEEE Trans. Wireless Commun., vol. 18, no. 11, pp. 5394–5409, Nov. 2019.
  • [8] H. Guo, Y.-C. Liang, J. Chen, and E. G. Larsson, “Weighted sum-rate maximization for intelligent reflecting surface enhanced wireless networks,” in IEEE Global Commun. Conf. (GLOBECOM), Dec. 2019.
  • [9] F. Fang, Y. Xu, Q.-V. Pham, and Z. Ding, “Energy-efficient design of IRS-NOMA networks,” IEEE Trans. Veh. Technol., vol. 69, no. 11, pp. 14 088–14 092, Nov. 2020.
  • [10] X. Mu, Y. Liu, L. Guo, J. Lin, and N. Al-Dhahir, “Exploiting intelligent reflecting surfaces in NOMA networks: Joint beamforming optimization,” IEEE Trans. Wireless Commun., vol. 19, no. 10, pp. 6884–2020, Oct. 2020.
  • [11] Y. Cao, T. Lv, and W. Ni, “Intelligent reflecting surface aided multi-user mmWave communications for coverage enhancement,” in IEEE Ann. Int. Symp. Personal Indoor Mobile Radio Commun., Aug. 2021.
  • [12] G. Zhou, C. Pan, H. Ren, K. Wang, M. D. Renzo, and A. Nallanathan, “Robust beamforming design for intelligent reflecting surface aided MISO communication systems,” IEEE Wireless Commun. Lett., vol. 9, no. 10, pp. 1658–1662, Oct. 2020.
  • [13] J. Hu, Y.-C. Liang, and Y. Pei, “Reconfigurable intelligent surface enhanced multi-user MISO symbiotic radio system,” IEEE Trans. Commun., vol. 69, no. 4, pp. 2359–2371, Apr. 2021.
  • [14] B. Zheng, C. You, and R. Zhang, “Double-IRS assisted multi-user MIMO: Cooperative passive beamforming design,” IEEE Trans. Wireless Commun., vol. 20, no. 7, pp. 4513–4526, Jul. 2021.
  • [15] H. Xie, J. Xu, and Y.-F. Liu, “Max-min fairness in IRS-aided multi-cell MISO systems with joint transmit and reflective beamforming,” IEEE Trans. Wireless Commun., vol. 20, no. 2, pp. 1379–1393, Feb. 2021.
  • [16] M.-M. Zhao, A. Liu, Y. Wan, and R. Zhang, “Two-timescale beamforming optimization for intelligent reflecting surface aided multiuser communication with QoS constraints,” IEEE Trans. Signal Process., 2021.
  • [17] J. Zhu, Y. Huang, J. Wang, K. Navaie, and Z. Ding, “Power efficient IRS-assisted NOMA,” IEEE Trans. Commun., vol. 69, no. 2, pp. 900–913, Feb. 2021.
  • [18] Z. Zhang and L. Dai, “A joint precoding framework for wideband reconfigurable intelligent surface-aided cell-free network,” IEEE Trans. Signal Process., vol. 69, pp. 4085–4101, Jun. 2021.
  • [19] S. Huang, Y. Ye, M. Xiao, H. V. Poor, and M. Skoglund, “Decentralized beamforming design for intelligent reflecting surface-enhanced cell-free networks,” IEEE Wireless Commun. Lett., vol. 10, no. 3, pp. 673–677, Mar. 2021.
  • [20] M. Zeng, X. Li, G. Li, W. Hao, and O. A. Dobre, “Sum rate maximization for IRS-assisted uplink NOMA,” IEEE Commun. Lett., vol. 25, no. 1, pp. 234–238, Jan. 2021.
  • [21] Y. Cao, T. Lv, Z. Lin, and W. Ni, “Delay-constrained joint power control, user detection and passive beamforming in intelligent reflecting surface-assisted uplink mmWave system,” IEEE Trans. Cognitive Commun. Netw., vol. 7, no. 2, pp. 482–495, Jun. 2021.
  • [22] M. Fu, Y. Zhou, and Y. Shi, “Reconfigurable intelligent surface for interference alignment in MIMO device-to-device networks,” in IEEE Int. Conf. Commun. (ICC) Workshops, Jun. 2021.
  • [23] K. Feng, X. Li, Y. Han, S. Jin, and Y. Chen, “Physical layer security enhancement exploiting intelligent reflecting surface,” IEEE Commun. Lett., vol. 25, no. 3, pp. 734–738, Mar. 2021.
  • [24] L. You, J. Xiong, D. W. K. Ng, C. Yuen, W. Wang, and X. Gao, “Energy efficiency and spectral efficiency tradeoff in RIS-aided multiuser MIMO uplink transmission,” IEEE Trans. Signal Process., vol. 69, pp. 1407–1421, Mar. 2021.
  • [25] W. Zhang, J. Xu, W. Xu, D. W. K. Ng, and H. Sun, “Cascaded channel estimation for IRS-assisted mmWave multi-antenna with quantized beamforming,” IEEE Commun. Lett., vol. 25, no. 2, pp. 593–597, Feb. 2021.
  • [26] C. You, B. Zheng, and R. Zhang, “Channel estimation and passive beamforming for intelligent reflecting surface: Discrete phase shift and progressive refinement,” IEEE J. Sel. Areas Commun., vol. 38, no. 11, pp. 2604–2620, Nov. 2020.
  • [27] X. Yu, D. Xu, and R. Schober, “Optimal beamforming for MISO communications via intelligent reflecting surfaces,” in IEEE Int. Workshop Sig. Process. Advances Wireless Commun. (SPAWC), May 2020.
  • [28] H. Wang, N. Shlezinger, Y. C. Eldar, S. Jin, M. F. Imani, I. Yoo, and D. R. Smith, “Dynamic metasurface antennas for MIMO-OFDM receivers with bit-limited ADCs,” IEEE Trans. Commun., vol. 69, no. 4, pp. 2643–2659, Apr. 2021.
  • [29] S. Abeywickrama, R. Zhang, Q. Wu, and C. Yuen, “Intelligent reflecting surface: Practical phase shift model and beamforming optimization,” IEEE Trans. Commun., vol. 68, no. 9, pp. 5849–5863, Sep. 2020.
  • [30] Q. Wu and R. Zhang, “Beamforming optimization for wireless networks aided by intelligent reflecting surface with discrete phase shifts,” IEEE Trans. Commun., vol. 68, no. 3, pp. 1838–1851, Mar. 2020.
  • [31] V. Arun and H. Balakrishnan, “RFocus: Beamforming using thousands of passive antennas,” in USENIX Symp. Netw. Sys. Design Implementation (NSDI), Feb. 2020, pp. 1047–1061.
  • [32] W. Cai, M. Li, and Q. Liu, “Practical modeling and beamforming for intelligent reflecting surface aided wideband systems,” IEEE Commun. Lett., vol. 24, no. 7, pp. 1568–1571, Jul. 2020.
  • [33] B. Di, H. Zhang, L. Song, Y. Li, Z. Han, and H. V. Poor, “Hybrid beamforming for reconfigurable intelligent surface based multi-user communications: Achievable rates with limited discrete phase shifts,” IEEE J. Sel. Areas Commun., vol. 38, no. 8, pp. 1809–1822, Aug. 2020.
  • [34] T. Shafique, H. Tabassum, and E. Hossain, “Optimization of wireless relaying with flexible UAV-borne reflecting surfaces,” IEEE Trans. Commun., vol. 69, no. 1, pp. 309–325, Jan. 2021.
  • [35] Z.-Q. Luo, W.-K. Ma, A. M. So, Y. Ye, and S. Zhang, “Semidefinite relaxation of quadratic optimization problems,” IEEE Signal Process. Mag., vol. 27, no. 3, pp. 20–34, May 2010.
  • [36] K. Shen and W. Yu, “Fractional programming for communication systems—Part I: Power control and beamforming,” IEEE Trans. Signal Process., vol. 66, no. 10, pp. 2616–2630, May 2018.
  • [37] ——, “Fractional programming for communication systems—Part II: Uplink scheduling via matching,” IEEE Trans. Signal Process., vol. 66, no. 10, pp. 2631–2644, May 2018.
  • [38] D. Mishra and H. Johansson, “Channel estimation and low-complexity beamforming design for passive intelligent surface assisted MISO wireless energy transfer,” in IEEE Int. Conf. Acoust. Speech Signal Process. (ICASSP), May 2019.
  • [39] S. Ren, K. Shen, Y. Zhang, X. Li, X. Chen, and Z.-Q. Luo, “Configuring intelligent reflecting surface with performance guarantees: Blind beamforming,” 2021, [Online]. Available https://kaimingshen.github.io/doc/BlindIRS.pdf.
  • [40] A. M.-C. So, J. Zhang, and Y. Ye, “On approximating complex quadratic optimization problems via semidefinite programming relaxations,” Math. Program., vol. 110, no. 1, pp. 93–110, 2007.
  • [41] T. Jiang, H. V. Cheng, and W. Yu, “Learning to reflect and to beamform for intelligent reflecting surface with implicit channel estimation,” IEEE J. Sel. Areas Commun., vol. 39, no. 6, pp. 1913–1945, Jul. 2021.
  • [42] T. L. Jensen and E. de Carvalho, “An optimal channel estimation scheme for intelligent reflecting surfaces based on a minimum variance unbiased estimator,” in IEEE Int. Conf. Acoust. Speech Signal Process. (ICASSP), May 2020.