arXiv:2004.13318v4 [cs.IT] 15 May 2021 EETASCIN NWRLS COMMUNICATIONS WIRELESS ON TRANSACTIONS IEEE nelgn eetn ufc:Sse oeigand Modeling System Surface: Reflecting Intelligent sgnrlysbpia ntrso hogptprui cost unit per throughput of while terms in costs, suboptimal individual densit generally IRS/BS their zero is (i.e., given IRS without cost network conventional deployment total a an exists ratio there that unveiled is it Moreover, h EEb edn eus [email protected] to b request must a purposes other sending any by for IEEE material the this use to permission niern,Ntoa nvriyo igpr,Singapore Singapore, 3610 of China [email protected]). University Electri University, National of Department Xiamen Information the Engineering, with is Education), Marine Zhang R. of [email protected]). and (Ministry Communication ogy Acoustic Underwater Nation R-26 and the R-261-518-005-720 and Grant Resea Research 20720190008), under Fundamental (No. Singapore the Universities Scienc 2019J05002), Central Natural (No. the the Province 61771017), No. Fujian and of 61801408 (No. China of eut nonly in results rwho ieesdt aebe otnal prigen- communi- spurring wireless innovative continually and been new inventing have in data thusiasm wireless of growth iepsiewrls ewr,sailtruhu,stoch analysis. throughput, performance spatial ometry, network, wireless tive/passive yrdntoktruhu gis h ovninlnetwor conventional which the IRS, against enhanci demonstr in without throughput and IRSs network analysis distributed deploying hybrid our of validate effectiveness numeric to Extensive the presented large. becomes are elements effects results IRS hardening of channel spatially number exhibit under the network which the channels, interf in and related need IRSs signal analysis distributed the performance by our both reflected of only, the distributions BSs power characterizing active the derive with on networks performance works wireless on key of prior based pertinent therein to BSs/IRSs other locati Compared random deployed as and the fading hybrid well of channel both new as over thro averaged a downlink metrics spatial paper achievable the its (BSs this in characterize stations and in base future. solutio IRSs, active the study cost-effective passive both in comprising we a network growth problem, be wireless capacity this can sustainable networks address their wireless achieve in to deplo IRSs large-scale the whether of problem syst open wireless an passiv remains IRS-aided still link-le various low-cost the for in optimization massive advances formance recent tuning Despite elements. via fav a reflecting constructing channels spectral by networks, the communication wireless improve of substantially efficiency energy to paradigm promising .Lui ihteSho fIfrais n e aoaoyo Laboratory Key and Informatics, of School the with is Lyu J. oyih 22 EE esnlueo hsmtra spermi is material this of use Personal IEEE. ©2021 Copyright hswr a upre npr yteNtoa aua Scien Natural National the by part in supported was work This yrdAtv/asv ieesNtokAddby Aided Network Wireless Active/Passive Hybrid h rlfrto fmbl plctosadexplosive and applications mobile of proliferation The Terms Index Abstract htmxmzstehbi ewr hogptsbetto subject throughput network hybrid the maximizes that Itlietrflcigsrae(R)i e and new a is (IRS) surface reflecting —Intelligent agnlyicesdinterference increased marginally Itlietrflcigsrae yrdac- hybrid surface, reflecting —Intelligent infiatybot h inlpower signal the boosts significantly .I I. NTRODUCTION inbnLyu, Jiangbin efrac Analysis Performance pia R/Sdensity IRS/BS optimal tcatcgeometry stochastic ebr IEEE Member, ntenetwork. the in a n Computer and cal lUiest of University al bandfrom obtained e . 3-000-E86-112. 153(email: 117583 td However, tted. c ud for Funds rch eFoundation ce Foundation e 5(e-mail: 05 si ge- astic e per- vel Technol- ratio) y ughput m,it ems, gthe ng erence orable yment and ) when cor- to s ons . but the ate nd To al k n e f . n u Zhang, Rui and , civ t utial aaiygot ntefuture. t the solution in cost-effective growth a capacity be sustainable its would achieve network nodes/components active wireless more the and in exis more mmWave the adding whether by fre- at doubtful paradigm is operating radio it when Therefore, costly bands. especially requiring frequency chains, antennas (RF) active mounti quency more as well with as base network more the them of in deployment relays the and/or (BSs) to stations due consumptio cost, energy hardware more higher increasingly and incurred also Although net they wireless SE, [1]. network the enhanced significantly communication ultra-dense technologies (mmWave) these include wave an decade (MIMO), millimeter multiple-output last technologies multiple-input massive wireless the (UDN), prominent in some proposed others, Among deployment/operati affordable cost. an spectral/en at higher yet (SE/EE) achieve efficiency would that technologies cation o-rhgnlmlil ces(OA 1][7,physica [14], [15]–[17], [13] (NOMA) access system multiple (OFDM) as non-orthogonal orthogona multiplexing such IRS-aided division setups, [6]–[12], frequency different design under beamforming and/or passive aspects s various wireless IRS-aided from investigating in recently and interest flexibility of high with as networks well wireless as in cost scalable deployed compatibility. and densely low be a to can at need IRS the such, standar operating without As and systems infrastructure seamless current WiFi and their or modify ceilings), cellular and object into walls existing integrated the (e.g., from environment removed the or in to attached easily Moreove be [2]. can nee (SIC) the cancellation without self-interference (FD) costly full-duplex in operat- efficiently t yet spectral consumption leading energy ing and thus cost reflection, hardware wave lower transmission/recept much passive signal leverages for simply chain but RF active any does IRS require relaying/beamforming, active conventional Co the cancellation. to beamform and/or nulling passive collaboratively interference (3D) spatial thus ing, three-dimensional as them, des such various and/or by achieve functions to amplitude reflected propagation the signal signal elements the alter reconfiguring the reflecting to of of dynamically passive tuned consisting phase low-cost be surface of can planar that number a uncontrollabl massive is and a random IRS be Specifically, to [2]–[5]. deemed traditionally environm is propagation new that wireless the on the based reconfiguring solution of promising concept a as emerged recently has otcl hscalne nelgn eetn ufc (I surface reflecting intelligent challenge, this tackle To h peln datgso R aeatatdagetdeal great a attracted have IRS of advantages appealing The elw IEEE Fellow, mpared ystems s[2]. ds ,IRS r, work ergy onal of d ting RS) ired not ion ent ng ly o o d n e s 1 - l l IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS 2 layer security [18]–[21], wireless information and power trans- fading effect. To our best knowledge, the modeling of a general fer [22]–[24], and so on. The above works on IRS-aided wire- multi-cell hybrid wireless network aided by randomly located less systems mainly aim to optimize the system performance IRSs and the characterization of the distribution of the users’ at the link level with one or more IRSs at fixed locations, achievable signal-to-interference-plus-noise ratio (SINR) as which show that the IRS-aided system can achieve significant well as the spatial throughput of the hybrid network have not energy efficiency [12] and/or spectral efficiency improvement been investigated yet in the literature. over the traditional system without IRS, with optimized IRS Motivated by the above, in this paper, we model a hybrid reflection coefficients. In [25], the authors investigate a multi- active/passive wireless network under the general multi-cell user system aided by multiple intelligent surfaces (equivalent setup and derive the distributions of the signal power, inter- to a single large IRS) co-located at a random location in ference power, and thereby the users’ achievable SINR in the the network. However, since IRS typically serves users in its network, with the ultimate goal of characterizing the spatial proximity, distributed IRSs should be deployed in the network throughput of the network, defined as the achievable rate per to serve distant groups of users and thereby boost the network user equipment (UE) averaged over both the wireless channel throughput. Motivated by this, in our prior work [26], the fading and the random BS/IRS locations. Thus, this work is a spatial throughput of a single-cell multi-user system aided by substantial extension of our prior work [26] under the single- distributed IRSs located at random locations is characterized, cell setup to the more general multi-cell setup. We focus on which is compared favorably with the conventional system the downlink communication from the BSs to the UEs while aided by distributed relays but with significantly reduced active the proposed analytical framework can be similarly extended antennas, under their respectively optimized deployment. to the uplink communication, which is left for our future Furthermore, for large-scale deployment of IRSs in future work. Compared to other prior works on characterizing the wireless systems, one critical issue is the modeling, design performance of wireless networks with active BSs only (see, and performance characterization of the IRS-aided multi-cell e.g., [33] and references therein), our analysis needs to derive hybrid wireless network comprising both distributed active the power distributions of both the signal and interference BSs and passive IRSs subjected to the inter-cell interference. reflected by distributed IRSs in the network under spatially There have been some recent works (e.g., [27]–[32]) along correlated channels, which exhibit channel hardening effects this line. Considering a finite number of co-channel/interfering when the number of IRS elements becomes large. Our main BSs, joint active/passive beamforming design with a cell-edge contributions are summarized as follows. IRS is investigated for the users’ weighted-sum-rate maximiza- tion [27] or minimum-rate maximization [28], respectively. In • First, we model the random BS/IRS locations by indepen- [29], the authors consider the quasi-static phase-shift design dent homogeneous Poisson point processes (HPPPs) and of one IRS in the presence of one interfering BS, based propose a practical UE-to-IRS association rule when they on the statistical channel state information (CSI) assuming are in close proximity. Then, for a typical BS-IRS-UE link given BS/IRS locations. In [30], the sum rate of multiple with their given locations, we derive its channel power transmit-receive (Tx-Rx) pairs aided by multiple distributed distribution in terms of the number of reflecting elements IRSs at given locations is maximized. However, the above per IRS, denoted by N, based on which the mean channel works only consider a given number of BSs and IRSs at fixed power is shown to scale with N in the order of O(N 2) locations, but do not investigate the impact of their spatial and O(N) for the cases with reflect beamforming by the random locations on the performance of large-scale hybrid associated IRS and random scattering by non-associated active/passive wireless networks. IRSs, respectively. Furthermore, we define the network Due to practical space constraints and heterogeneous/dense coverage probability and spatial throughput in terms of BS deployment, modern cellular networks typically exhibit an key system parameters including the BS/IRS densities and increasing degree of spatial irregularity, for which the con- network loading factor. ventional grid-based BS deployment models become no more • Next, we propose an analytical framework for the IRS- suitable. As field trials are costly and system-level simulations aided hybrid network based on stochastic geometry, and are time-consuming, stochastic geometry has been extensively address its new challenges. In particular, for a typical UE applied as a tractable analytical tool to model the spatial 0 with randomly distributed IRSs nearby, its nearest IRS 0 distribution of heterogeneously/densely/irregularly deployed (and hence the IRS 0-UE 0 distance d0) typically has the wireless nodes, which provides meaningful performance lower dominant impact on the mean signal or interference power bounds and scaling laws for practical wireless networks [33] compared to other (farther) IRSs, under practical values of [34]. Based on stochastic geometry, the authors in [31] model the IRS density. Moreover, with reflect beamforming by the IRSs by boolean line segments in a large-scale network the associated IRS, the signal link exhibits channel hard- and derive the probability that a given IRS is capable of ening when N becomes large, while the extent of channel providing an indirect path for a given Tx-Rx pair (i.e., the hardening varies with the distance d0, rendering it difficult reflection probability). The authors in [32] further exploit to characterize the signal power distribution and thus the the deployment of IRSs for providing indirect line-of-sight SINR distribution. To overcome this difficulty, we propose (LoS) paths for blocked links, thus improving the coverage to approximate the conditional signal power distribution probability in a large-scale network. However, these two works by the Gamma distribution, whose shape parameter kS do not consider the inter-cell interference and the small-scale specifies the extent of channel hardening conditioned on IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS 3
d0. The conditional SINR distribution is then obtained in BS 0 terms of the interference power Laplace transform and its signal derivatives up to integer-order kS. Moreover, we propose an interpolation method for non-integer kS, and apply the IRS j reflect normal approximation of the signal power in the case with scatteringg beamforming large kS in order to reduce the computational complexity. IRS 0 These new analytical methods jointly yield accurate and UE 0 efficient characterization of the network SINR distribution and hence its spatial throughput. • Finally, extensive numerical results are provided to validate our analytical results. It is found that increasing IRS density interference in a hybrid wireless network can significantly enhance the signal power but with only marginally increased interfer- BS m ence, thus greatly improving its throughput as compared Fig. 1: IRS-aided multi-cell wireless network in the downlink. to the traditional wireless network with active BSs only, especially when the BS density, network loading factor or N is large. Moreover, it is unveiled that there exists an division multiple access (OFDMA) scheme and assume that optimal IRS/BS density ratio ζ∗ for maximizing the spatial the transmission bandwidth and each time slot are equally throughput of the new hybrid network under a given total divided into orthogonal resource blocks (RBs), each randomly deployment cost, where ζ∗ is shown to increase with the assigned by a BS to one of its served UEs, over which the BS/IRS cost ratio and the network loading factor, while the channel is assumed to be frequency-flat and constant, while the conventional network without IRS (i.e., zero IRS/BS density channels may vary over different frequency bands or different ratio), or the hybrid network with excessively large ζ (where time slots. We assume that the network has a homogeneous the BS density is too low to provide enough signal power for traffic load, where all BSs have a common loading factor p effective IRS passive beamforming), is generally suboptimal (0 < p ≤ 1), i.e., each of the RBs is active with probability in terms of throughput per unit cost. Furthermore, it is p independently. Consider one typical RB used by the typical shown that the maximum spatial throughput of the hybrid UE 0 located at the origin, which is associated with its nearest network with the optimal IRS/BS density ratio ζ∗ grows BS 0 with distance l0, as shown in Fig. 1. As a result, the BSs ′ almost linearly with the total cost, thus providing a new that transmit on the same RB form a thinned HPPP ΛB with ′ , and cost-effective approach to achieve sustainable capacity density λB pλB. growth for future wireless networks. We consider that distributed IRSs are deployed to assist the BS-UE communications in the network. Assume that all The rest of this paper is organized as follows. The new IRSs are of the same height equal to m,2 while the model of the proposed hybrid wireless network is presented HI IRSs’ horizontal locations are modeled by a 2D HPPP in Section II. The distributions/mean values of the signal ΛI (independent of ) on the ground plane with given density and interference powers are then characterized in Section III ΛB IRSs/m2. Denote the set of IRS horizontal locations as and Section IV, respectively. Next, the SINR distribution and λI , w R2 , where w is the 2D coordinate the network spatial throughput are obtained in Section V. W { j ∈ |j ∈ ΛI} j of an IRS . Denote as the horizontal distance Numerical results are provided in Section VI. Finally, we j ∈ ΛI dj between UE 0 and IRS , as shown in Fig. 1. Since the IRS conclude the paper in Section VII. j typically provides signal enhancement via reflect beamforming in a local region [26], we consider the practical scenario II. SYSTEM MODEL where UE 0 is associated with its nearest IRS 0 for dedicated In this paper, we consider an IRS-aided multi-cell wire- reflect beamforming, if IRS 0 is within a certain threshold 3 less network shown in Fig. 1, and focus on the downlink distance D1, i.e., d0 ≤ D1. On the other hand, if there communication from the BSs to UEs. Assume that the BSs is no IRS within D1, then UE 0 is served by BS 0 only 4 are of the same height equal to HB meters (m), while the without any associated IRS. Moreover, for the purpose of BSs’ horizontal locations are modeled by a 2-dimensional exposition, we assume that each IRS is always on and reflects (2D) HPPP ΛB on the ground plane with given density 2 2 λB BSs/m . Assume that the UE locations follow another Our analytical framework is applicable to any given height of BSs/IRSs, independent HPPP on the ground plane,1 such that we can which can also be extended to account for their random heights by employing the 3-dimensional point process where the height of each BS/IRS is randomly focus on one typical UE 0 in this HPPP to analyze the average set within a certain range. 3 UE performance without changing the location distribution The threshold D1 is practically set such that each IRS serves a finite of other UEs according to the Slivnyak’s theorem [33]. To number of UEs in its neighborhood only. 4To implement this, each IRS controller can sense the nearby UEs and facilitate our analysis, we consider the orthogonal frequency decide to associate with them or not based on their signal strengths, then send the associated UEs’ identification back to the BS (via a separate control 1The HPPP assumption implies the uniformly random UE distribution, link) for RB allocations. Note that such IRS-UE associations change only which is more accurate for UEs of homogeneous distribution and/or higher when the UEs move in/out from the coverage of each IRS, which usually mobility, but in general can serve as a good baseline to evaluate the network happens not so frequently for typical IRS-aided scenarios (e.g., hotspot with performance with heterogeneous user distribution/mobility in practice. local users). IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS 4 the received wave at all time, regardless of whether there is any framework by treating it as equivalent random perturbation in UE associated with it.5 As a result, UE 0 receives the reflected the locations of the BSs/IRSs (see Sections III.G and VI.A in signal/interference by all IRSs (including its associated IRS 0 [35] and the references therein), which is ignored in this work if any). To model the randomly reflected signal/interference for simplicity. The channel power gain from BS m to UE 0 by non-serving IRSs accurately while maintaining analytical is thus given by Approximation 1 tractability, we adopt the that only the IRSs 2 , 2 2 −α/2 |hd,m| gd,mξd,m = β(lm + HB) ξd,m, (2) within a sufficiently large threshold distance D2 (D2 > D1) , from UE 0, denoted by the set J {j ∈ ΛI|dj ≤ D2}, will where gd,m is the average channel power gain, ξd,m accounts contribute the signal/interference to it. Finally, to maximize for channel fading, lm denotes the BS-UE horizontal distance, the passive beamforming gain of the IRS to each served UE, 4πfc −2 and β = ( c ) denotes the average channel power gain at we assume that its served UEs are assigned in orthogonal- a reference distance of 1 m based on the free-space path-loss time RBs, i.e., time division multiple access (TDMA) or time model, with fc denoting the carrier frequency, and c denoting sharing is adopted for the UEs served by the same IRS.6 the speed of light. Similarly, the channel power gains from BS m to the n-th element of IRS j, and from the latter to UE 0 A. Channel Model are given by
Assume for simplicity that the BSs and UEs are each (j) 2 (j) (j) 2 2 −α/2 (j) |h | , g ξ = β r + (HB − HI) ξ , equipped with a single antenna, while each IRS has N i,m,n i,m i,m,n m,j i,m,n (3) reflecting elements. The baseband equivalent channels from and −α/2 BS m to IRS j, from IRS j to UE 0, and from BS m to UE 0 |h(j)|2 , g(j)ξ(j) = β d2 + H2 ξ(j), (4) are denoted by h(j) , [h(j) , ··· ,h(j) ]T ∈ CN×1, h(j) , r,n r r,n j I r,n i,m i,m,1 i,m,N r (j) (j) T CN×1 C 7 where g(j) and g(j) denote the average channel power gains, [hr,1 , ··· ,hr,N ] ∈ , and hd,m ∈ , respectively, i,m r C T (j) (j) where denotes the set of complex numbers and [·] denotes ξi,m,n and ξr,n account for channel fading, rm,j and dj denote (j) , (j) (j) 9 the matrix transpose. Let φ [φ1 , ··· , φN ] and further the BS-IRS and IRS-UE horizontal distances, respectively. (j) iφ(j) iφ(j) denote Φ , diag{[e 1 , ··· ,e N ]} (with i denoting the imaginary unit) as the phase-shifting matrix of IRS j, (j) B. BS-IRS-UE Channel Power Statistics where φn ∈ [0, 2π) is the phase shift by element n of the IRS on the incident signal,8 and diag{x} denotes a diagonal In this subsection, we derive the BS-IRS-UE cascaded matrix with each diagonal element being the corresponding channel power statistics, which is new for the IRS-aided hybrid element in x. Each element of the IRS receives the superposed network and essential to our subsequent performance analysis multi-path signals from the BS, and scatters the combined (j) (j) for it. Assume that the channels hd,m, hi,m,n and hr,m,n, signal with adjustable phase as if from a single point source. m ∈ ΛB, j ∈ ΛI, n = 1, ··· ,N are independent. For the Therefore, the cascaded BS-IRS-UE channel can be modeled cascaded BS m-IRS j-UE 0 link in (1), the channel reflected as a concatenation of three components, namely, BS-IRS link, through each element n is given by IRS reflecting with phase shifts, and IRS-UE link, given by (j) (j) i (j) , (j) φn [6] hir,m,n hi,m,nhr,ne i (j) ∠ (j) ∠ (j) N (j) (j) φn + hi,m,n+ hr,n (j) (j) (j) i (j) = |hi,m,n||hr,n|e , (5) , h T Φ(j)h(j) (j) φn hir,m [ i,m] r = hi,m,nhr,ne ,m ∈ ΛB, nX=1 where the channel amplitude |h(j) | , |h(j) ||h(j)| resem- (1) ir,m,n i,m,n r,n (j) (j) (j) i (j) ∠ (j) , , φn bles a doubly-faded RV while the channel phase h where hir,m,n hi,m,nhr,ne denotes the BS m-IRS j-UE ir,m,n (j) ∠ (j) ∠ (j) 0 channel reflected by element n, n =1, ··· ,N. φn + hi,m,n + hr,n is adjustable via controlling the phase (j) For the BS-IRS, IRS-UE and BS-UE links, we assume a shift φn exerted by IRS j. simplified fading channel model without shadowing, which dist. For simplicity, assume that ξd,m = ξ ∼ Exp(1) is an consists of distance-dependent path loss with path-loss expo- exponential random variable (RV) with unit mean accounting nent α ≥ 2 and an additional random term ξ accounting for for the small-scale Rayleigh fading. Therefore, the ampli- small-scale fading. Note that the shadowing effect can also tude |hd | follows the Rayleigh distribution [36] with scale be incorporated into the stochastic geometry based analytical ,m parameter gd,m/2, denoted by R gd,m/2 , while hd,m 5The results in this paper can be extended to the general case where each follows thep circularly symmetric complex p Gaussian (CSCG) IRS is independently on or off with a certain probability. Nevertheless, we distribution [36] with mean zero and covariance gd,m, denoted consider that all IRSs are on to characterize the worst-case interference. by . We also assume Rayleigh faded channel for 6 CN (0,gd,m) It is shown in [17] that for IRS-aided multiple access, the TDMA scheme (j) (j) dist. is in general superior over the FDMA scheme due to the hardware limitation of the BS-IRS and IRS-UE links, i.e., ξi,m,n, ξr,n = ξ ∼ Exp(1), IRS passive reflection, which can be made time-selective, but not frequency- in order to investigate the worse-case propagation condition selective [2]. for the IRS and hence characterize the achievable performance 7The subscripts “i”, “r” and “d” represent the BS-to-IRS channel, IRS- reflected channel (i.e., the IRS-to-UE channel), and direct BS-to-UE channel, respectively. 9For the purpose of exposition, we consider far-field propagation for all 8 In this paper, we assume (maximum) unit amplitude for each reflection links, and accordingly assume HB ≥ 1 m and HI ≥ 1 m, which also avoid coefficient to maximize the IRS beamforming gain to its served UE [6]. unbounded power gain when the horizontal distance lm or dj becomes zero. IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS 5 lower bound for the IRS-aided hybrid network.10 Therefore, The average BS 0-IRS 0-UE 0 signal power is the second (j) (j) (j) (j) (0) moment of |hir | and thus given by we have |hi,m,n| ∼ R gi,m/2 and |hr,n| ∼ R gr /2 . ,0 q q As a result, the channel amplitude (j) in (5) is a (0) , E (0) 2 E (0) 2 (0) |hir,m,n| gir,0 {|hir,0| } = {|hir,0|} + var{|hir,0|} double-Rayleigh RV. Note that for a double-Rayleigh dis- 2 2 π π (0) tributed RV Y = X X with independent X ∼ R(δ ) and = N 2 + 1 − N g g(0), (12) 1 2 1 1 16 16 i,0 r X2 ∼ R(δ2), its mean and variance are respectively given by [37] which is proportional to the average channel power product E , (0) (0) , π2 2 {Y } πδ1δ2/2, (6) gi,0 gr , with the beamforming gain coefficient Gbf 16 N + π2 2 , 2 2 2 1 − 16 N that grows with N in the order of O(N ). var{Y } 4δ1δ2 (1 − π /16). (7) On the other hand, for any IRS j that does not provide Therefore, the mean and variance of the channel amplitude reflect beamforming for UE 0 (including IRS 0 if the dis- (j) |hir,m,n| are respectively given by tance d0 > D1), it scatters the incoming signal from BS m without passive beamforming, thus resulting in uniformly E (j) , π (j) (j) (j) hir gi gr , (8) random channel phase ∠h due to the uniformly random ,m,n 4 q ,m ir,m,n phases ∠h(j) and ∠h(j). In this case, we have the following ( j) , 2 (j) (j) i,m,n r,n var hir,m,n (1 − π /16)gi,mgr . (9) proposition.
In the case where IRS 0 provides reflect beamforming ser- Proposition 2. With random scattering by IRS j, the BS m- vice for the desired signal from BS 0 to UE 0 (i.e., d0 ≤ D1), IRS j-UE 0 channel in (1) for practically large N can be ∠ (0) (0) we assume that the cascaded channel phase hi,0,nhr,n for approximated by the CSCG distribution, i.e., each reflected path n = 1, ··· ,N can be obtained via IRS- N customized channel estimation methods (please refer to [13] approx. h(j) = h(j) ∼ CN 0,Ng(j) g(j) . (13) [14] for more details). As a result, IRS 0 can then adjust ir,m ir,m,n i,m r nX=1 the phase shift φ(0) such that the N reflected paths of the desired signal are of the same phase at UE 0’s receiver by (j) Proof: The channel hir,m,n reflected by each element n (0) ∠ (0) (0) 11 setting φn = − hi,0,nhr,n ,n = 1, ··· ,N. Therefore, has zero mean and independent in-phase and quadrature- the amplitude of the BS 0-IRS 0-UE 0 channel is given by 1 (j) (j) phase components each with variance 2 gi,mgr , respec- N tively, with the detailed derivations given in Appendix A. (0) h(0) T h(0) (0) (j) |hir,0| = | i,0 | | r | = hir,0,n . (10) Moreover, since the channels hir,m,n, n = 1, ··· ,N are nX=1 i.i.d., based on Approximation 2, the independent in-phase
Furthermore, we adopt the Approximation 2 that by the and quadrature-phase components of the combined channel (j) N (j) central limit theorem (CLT), given N i.i.d. RVs X1, ··· ,XN hir,m = n=1 hir,m,n can be each approximated by an inde- each with mean and variance 2, the sum N P 1 (j) (j) µ ω Y = n=1 Xn pendent normal distribution N (0, 2 Ngi,mgr ) for practically can be approximated by the normal distribution N (Nµ,NωP 2) large N. As a result, the combined BS m-IRS j-UE 0 channel for sufficiently large N. As a result, we have can be approximated by the CSCG distribution given by (13). Proposition 2 is thus proved. Proposition 1. With reflect beamforming by IRS 0, the BS 0- IRS 0-UE 0 channel amplitude in (10) for practically large Therefore, the average BS m-IRS j-UE 0 channel power is N can be approximated by the normal/Gaussian distribution, given by (j) , E (j) 2 (j) (j) i.e., gir,m {|hir,m| } = Ngi,mgr , (14) which is proportional to the average channel power product (0) approx. E (0) (0) |hir,0| ∼ N N hir,0,n ,N var hir,0,n (j) (j) , gi gr , with the scattering gain coefficient Gsc N that ,m 2 grows linearly with N. π (0) (0) π (0) (0) = N N g gr ,N 1 − g g . (11) 4 i,0 16 i,0 r q Proof: Based on Approximation 2 with large N,12 the BS 0-IRS 0-UE 0 channel amplitude in (10) is the sum of N i.i.d. C. SINR, Coverage Probability, and Spatial Throughput double-Rayleigh RVs h(0) , n =1, ··· ,N, each with mean ir,0,n Denote the downlink transmit power on each RB as P0. and variance given by (8) and (9), respectively. Therefore, Then, the overall signal power (normalized by ) from BS 0 P0 Proposition 1 is proved. to UE 0 is given by 2 10 The proposed analytical method in this paper can be extended to other S , h + h(j) , (15) fading channel models such as Rician fading. d,0 ir,0 11 jX∈J For the ease of practical implementation, we consider reflect beamforming for enhancing the desired signal power only, instead of nulling any co-channel which accounts for the direct path and reflected paths via all interference. 12According to our simulations, when N > 25, this approximation is IRSs j ∈ J . Further denote the total received interference already quite accurate for our considered setup. power (normalized by P0 as well) from all co-channel BSs IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS 6
′ m ∈ ΛB \{0} by and IV, respectively, and then deriving the SINR distribution 2 and hence the spatial throughput in Section V. In particular, I , I = h + h(j) . (16) for the signal (or interference) power characterization, we first m d,m ir,m m∈XΛ′ \{0} m∈XΛ′ \{0} jX∈J derive its conditional distribution conditioned on the distances B B l0 and d0, based on which we are then able to obtain the SINR The received SINR at UE 0 is thus given by distribution. Moreover, we also characterize the unconditional S mean signal/interference power averaged over the channel γ , , (17) I + W fading and random BS/IRS locations, in order to reveal the impact of IRS on them, which helps illustrating their respective where W , σ2/P , and the receiver noise is assumed to be 0 effects to the SINR distribution. additive white Gaussian noise (AWGN) with power σ2. The corresponding achievable rate in bits/second/Hz III. SIGNAL POWER DISTRIBUTION (bps/Hz) is given by In this section, we characterize the conditional signal power , distribution as well as the (unconditional) mean signal power. R log2(1 + γ). (18) Since in our considered IRS-aided hybrid network, the exact Note that the signal power S, interference power I, and signal power distribution entails a more complicated form (as thus SINR γ and achievable rate R are all RVs depending will be shown in Section III-A) as compared to the conven- on the random channel fading as well as random BS/IRS tional case without IRS, the well-known analytical method locations. An outage event occurs when the rate R is lower proposed in [33] for deriving the SINR distribution directly than a minimum required target R¯. The coverage probability cannot be applied in our context. Therefore, we propose of the typical UE 0 is then defined as the average non- to approximate the conditional signal power distribution by outage probability over the random channel fading and random the Gamma distribution (which belongs to the exponential BS/IRS locations, i.e., distribution family [36]) based on its first and second moments conditioned on l0 and d0, under three different cases based on Pcov , P{R ≥ R¯} = P{γ ≥ γ¯} , 1 − Fγ (¯γ), (19) d0, i.e., with IRS reflect beamforming (d0 ≤ D1), with IRS , R¯ where γ¯ 2 − 1 denotes the corresponding minimum re- scattering only (D1 < d0 ≤ D2), and without any nearby quired SINR, and Fγ (·) is the cumulative distribution function IRS (d0 >D2). Such Gamma approximation incorporates the (cdf)13 of γ. Accordingly, we can define the spatial throughput BS-IRS-UE channel power statistics derived in Section II-B of the network in bps/Hz/m2 as as well as the impact of randomly distributed IRSs, which is tailored to the new IRS-aided hybrid network and thus differs ν , P Rλ¯ ′ = P Rpλ¯ . (20) cov B cov B from that for the conventional network without IRS. Moreover, In order to obtain the spatial throughput ν, we need to for the case with IRS reflect beamforming, we characterize the characterize the cdf of the SINR γ that depends on the impact of d0 on the mean signal power and the d0-dependent distributions of the signal power S and interference power channel hardening effect to draw useful insights. Finally, the I. Note that the instantaneous S and I are independent due mean signal power is obtained by integrating the conditional to their independent small-scale fading, while their large-scale mean signal power over the distributions of l0 and d0. statistics averaged over fading are dependent in general due Note that the first and second moments of the signal power to the common BS and IRS locations. Specifically, under the S depend on the average channel power gains gd,0 (BS 0-UE 0 (j) (j) distance-based association rule, the distribution of the BS 0- link), gi,0 (BS 0-IRS j link) and gr (IRS j-UE 0 link), which UE 0 link distance l0 affects not only the mean signal power, further depend on the corresponding horizontal link distance 2 but also the mean interference power since the interfering BSs l0, r0,j and dj that are related by the cosine law, i.e., r0,j = 2 2 are located at distances more than l0 from UE 0. Moreover, l0 + dj − 2l0dj cos ϕ, where ϕ is the BS 0-UE 0-IRS j angle for the IRS-aided downlink communication, both the signal projected on the ground plane. The exact expressions for the and interference are reflected by the same set (J ) of IRSs moments of S conditioned on l0 and d0, though not expressible within D2 of UE 0, and hence the path-losses of the BS-IRS in closed-form, can be obtained by numerical integrals over and/or IRS-UE links are random but correlated in general. ϕ and/or dj (as illustrated in Appendix B). Nevertheless, for Such correlation introduced by randomly distributed IRSs near simplicity in this section, we apply the Approximation 3 that (j) UE 0 imposes new difficulty to the system-level performance r0,j ≈ l0 and hence gi,0 ≈ gd,0, j ∈ J , in order to obtain analysis, which is challenging to deal with. In addition, the closed-form approximations for the moments of S|l0,d0 . Such IRS 0-UE 0 distance d0 determines whether there is reflect approximation is reasonable since the IRSs considered here are beamforming provided by IRS 0 and thus the IRS 0-UE 0 in the local region of UE 0 (and hence with small dj ), which channel power gain, which also has a significant impact on is also verified by our numerical examples later in Section VI the system performance. where the analytical results are shown to match well with the To tackle the above challenges, we decompose the perfor- Monte Carlo (MC) simulation results. mance analysis into three parts, by first characterizing the sig- nal power and interference power distributions in Sections III A. The Case with IRS Reflect Beamforming 1) Gamma Approximation with Moment Matching: First, 13 The cdf of an RV X is defined as FX (x) , P{X