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

Wireless Performance Evaluation of Building Layouts: Closed-Form Computation
of Figures of Merit

Jiliang Zhang, Senior Member IEEE, Andrés Alayón Glazunov, Senior Member IEEE, and Jie Zhang, Senior Member IEEE Jiliang Zhang is with the Department of Electronic and Electrical Engineering, The University of Sheffield, Sheffield, UK. Andrés Alayón Glazunov is with Department of Electrical Engineering, University of Twente, Enschede, Overijssel, Netherlands, and also with Department of Electrical Engineering, Chalmers University of Technology, Gothenburg, Sweden. Jie Zhang is with the Department of Electronic and Electrical Engineering, the University of Sheffield, Sheffield, UK, and also with Ranplan Wireless Network Design Ltd., Cambridge, UK. The research is funded by EUROSTARS Project BuildWise (11088).
Abstract

This paper presents a part of our ground-breaking work on evaluation of buildings in terms of wireless friendliness in the building-design stage. The main goal is to devise construction practices that provide for a good performance of wireless networks deployed in buildings. In this paper, the interference gain (IG) and power gain (PG) are defined as two figures of merit (FoM) of the wireless performance of buildings. The FoMs bridge the gap between building design and wireless communications industries. An approach to derive exact closed-form equations for these FoMs is proposed for the first time. The derived analytic expressions facilitate straightforward and more computationally efficient numerical evaluation of the proposed FoMs as compared to Monte Carlo simulations for well-known indoor propagation models. It is shown that the derived closed-form expression can be readily employed to evaluate the impact of building properties, such as the sizes and the aspect ratios (ARs) of rooms, on the wireless performance. The proposed approach sheds light to architects on evaluation and design of wireless-friendly building layouts.

Index Terms:
Smart buildings, building design, wireless performance evaluation, 5G networks.

I Introduction

It has been predicted that mobile traffic will increase by up to 10001000 times in the next decade, with over 80%80\% of mobile traffic taking place indoors [1]. 6G systems will meet the required mobile traffic through the combined effects of densification of the wireless network through small cells and the large-scale deployment of massive multiple-input multiple-output (MIMO) array antennas [2]. The large number of small cells are anticipated to enhance coverage and spectral efficiency [3, 4, 5, 6, 7, 8, 9], whereas in massive MIMO, array antennas are scaled up to increase reliability and spectral efficiency [11, 12, 13, 10]. On the other hand, wireless communication plays an important role in the realization of the “smart building” paradigm, or more broadly, in the “smart city” vision [14]. Different types of buildings have an intrinsic wireless performance which is independent of how densely small cells are deployed, or the massive MIMO array antenna deployed [15, 16]. The present work is the first step towards providing guidance, to both civil engineers and architects, on incorporating the wireless performance of buildings at the design stage.

Building design and construction has mainly focused on improving safety, visual and thermal comfort, and indoor air quality [17]. In the last decade, energy efficiency became an important metric in building design tools [18]. The operational energy required by a building can be evaluated by taking into account the lighting, the heating, the ventilation and the cooling as well as the provision of hot water [19]. At present, mobile operators are already among the top energy consumers [20], and an increasing trend of power consumption of wireless networks can be observed in different generations of wireless systems [21, 22]. Therefore, the energy saving is worthwhile to be considered in the design of green buildings. However, the indoor power consumption of wireless networks has not been considered.

One important factor impacting the power consumption in small-cell networks is the building-blockage. In [23, 24, 25], the LOS probability is analytically and empirically modelled. In [26], an analytic expression for the coverage rate is provided by approximating the line-of-sight (LOS) region as a fixed LOS ball. In [27], four random wall generation methods to place blockage objects in an indoor scenario were proposed and analyzed, where walls are either distributed randomly, semi-deterministically or heuristically. However, previous evaluations of the blockage effect proposed that buildings form a random process of shapes. The LOS probability was there assumed to be irrelevant for different directions of the link. This model seems to be promising in providing insights on the configuration of a wireless network in a homogeneous indoor scenario, where walls are assumed to be uniformly distributed in a building. However, in specific building designs, whether a wireless link is blocked by structures not only depends on the distance between the transmitter and the receiver, but also on their relative position and direction in the floor plan of the building. That is, the blockage condition of the wireless signal depends on the considered direction. Therefore, the homogeneous assumption is not suitable for evaluating specific building designs.

Rather than random blockages, ray-tracing/ray-launching techniques are employed to evaluate the indoor performance of wireless networks (i.e., for users located indoors) for specific building designs in [28, 29, 30, 31, 32, 33, 34]. However, results attained through this approach are only valid for a specific network deployment, and are not general enough to evaluate the wireless performance of buildings, where all possible positions of transmitters and receivers in the building under design are needed to be evaluated. In this case, ray-tracing/ray-launching techniques become computationally expensive, and therefore it is not practical to evaluate the wireless performance of a building layout in a reasonable time.

Moreover, wireless networks are traditionally designed and deployed when the buildings have been built. Before the networks are deployed, wireless communications engineers have to carry out ray-tracing/ray-launching based simulations for specific deployments of wireless networks. In other words, wireless communications engineers adapt a network to suit its deployment scenarios without modifying any property of the building layout. However, for an existing building with a specific intrinsic wireless performance determined by the building layout, the indoor wireless networks may not achieve an acceptable performance regardless how densely small cells or MIMO array antenna elements are deployed. The wireless performance of building layout has to be taken into account carefully in the building design stage, which means that building designers, both civil engineers and architects, requires guidance on how to predict and evaluate the wireless performance of a building under design.

Taking wireless performance of buildings into account will become indispensable for the design of future smart buildings in which appliances and many other devices and “things” will all be connected wirelessly. To address the challenges posed by the building evaluation in terms of wireless performance, ground-breaking works are needed. We are engaged in the collaborative project Build-Wise [35]. Supported by this project, figures of merits (FoMs) reflecting the intrinsic wireless performance of buildings have been defined considering the impact of blockage, bounce reflections, and waveguide effects introduced by indoor building structures employing the open space scenario as a benchmark for wireless network performance evaluation [15, 16]. The proposed FoMs bridge the gap between building design and wireless communications, laying the foundation to an improved wireless service indoors by optimising building layouts. However, the computation of the FoMs is based on Monte Carlo simulations, which is computationally expensive since an infinite number of transmit elements is assumed in the FoM definition. Therefore, an analytic framework is critically needed for fast computation of the FoMs.

The main contribution of this paper is a mathematical framework to evaluate the wireless performance of buildings. The proposed approach facilitates straightforward and more computationally efficient numerical evaluation of the proposed FoMs. Firstly, two FoMs, i.e., the interference gain (IG) and power gain (PG) are defined to capture the impact of the building layout on the power of interference signals and intended signals, respectively. Secondly, an approach to derive exact FoMs in closed form is proposed for polygonal-shaped rooms. By using this approach, the FoMs can be quickly and exactly computed given building layout tiled by polygonal rooms. Thirdly, the analytic approach is validated through Monte Carlo simulations for typical indoor path gain models. This proposed analytic framework provides a direct and practical quantification of the wireless performance of a building. Finally, the impact of the building layout, e.g., their sizes and the aspect ratios (AR) of rectangular rooms, on wireless performance FoMs is analyzed based on the proposed analytic approach.

The remainder of this paper is organized as follows. Section II summarizes the notation and main assumptions made in this paper. In Section III, two FoMs, i.e., the IG, and the PG are introduced to measure the wireless performance of a building. In Section IV, an approach to compute the closed-form expressions of the proposed FoMs is provided. In Section V, the approach is validated through Monte Carlo simulations. In Section VI, the impact of building layout on FoMs is analyzed. Conclusions are provided in Section VII.

II Notation and assumption

II-A Notation

We use the following notations throughout this paper. 1) Polylogarithm Lis(z)k=1zkks\mathrm{Li}_{s}(z)\triangleq\sum\limits_{k=1}^{\infty}{\frac{z^{k}}{k^{s}}}. 2) Beta function B(x,y)01tx1(1t)y1dt=Γ(x)Γ(y)Γ(x+y)B(x,y)\triangleq\int_{0}^{1}t^{x-1}(1-t)^{y-1}\mathrm{d}t=\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}. 3) Gamma function Γ(z)0xz1exdx=1zn=1(1+1n)z1+zn\Gamma(z)\triangleq\int_{0}^{\infty}x^{z-1}e^{-x}\mathrm{d}x=\frac{1}{z}\prod\limits_{n=1}^{\infty}\frac{\left(1+\frac{1}{n}\right)^{z}}{1+\frac{z}{n}}. 4) Hypergeometric function F12(a,b;c;z)1B(b,cb)01xb1(1x)cb1(1zx)adx=n=0(a)n(b)n(c)nznn!{}_{2}F_{1}(a,b;c;z)\triangleq\frac{1}{B(b,c-b)}\int_{0}^{1}x^{b-1}(1-x)^{c-b-1}(1-zx)^{-a}\mathrm{d}x=\sum\limits_{n=0}^{\infty}\frac{(a)_{n}(b)_{n}}{(c)_{n}}\frac{z^{n}}{n!}, where (x)nk=0n1(xk)(x)_{n}\triangleq\prod\limits_{k=0}^{n-1}(x-k). 5) “\land”, “\lor”, and “¬\lnot” denote the logical conjunction, the logical disjunction, and the logical complement, respectively. “\top” and “\bot” denote the logical tautology and logical contradiction, respectively. 6) Im[z]\mathrm{Im}[z] denotes the imaginary part of a complex number zz. 7) In the open space scenario, IOI_{\mathrm{O}} and POP_{\mathrm{O}} denote the power of interference signals and the power of intended signals, respectively. In the LOS scenario, ILI_{\mathrm{L}} and PLP_{\mathrm{L}} denote the power of interference signals and the power of intended signals, respectively. In the non-LOS (NLOS) scenario, INI_{\mathrm{N}} and PNP_{\mathrm{N}} denote the power of interference signals and the power of intended signals, respectively. 8) γO\gamma_{\mathrm{O}} and γB\gamma_{\mathrm{B}} denote the signal-to-interference-plus-noise-ratio (SINR) in the open space and the indoor scenario, respectively. 9) GOG_{\mathrm{O}}, GLG_{\mathrm{L}}, and GNG_{\mathrm{N}} denote path gain models in the open space, the LOS, and the NLOS scenarios, respectively. 10) PTP_{\mathrm{T}} [Wm2][\mathrm{W}\mathrm{m}^{-2}] denotes the transmit power density, i.e., transmit power per unit area. 11) PthP_{\mathrm{th}} [Wm2][\mathrm{W}\mathrm{m}^{-2}] denotes the threshold power level determined by the sensitivity of the receiver also defined per unit area. 12) ROR_{\mathrm{O}}, RLR_{\mathrm{L}} and RNR_{\mathrm{N}} denote the coverage distances under open space, LOS and NLOS scenarios, respectively. 13) gIg_{\mathrm{I}} and gPg_{\mathrm{P}} denote the IG and the PG of buildings, respectively. 14) δ(ξ)\delta(\xi) is defined as a function of inequality ξ\xi. If ξ\xi holds, δ(ξ)=1\delta(\xi)=1. Otherwise, δ(ξ)=0\delta(\xi)=0.

II-B Assumptions

We use the following assumptions throughout this paper.

  1. 1.

    The FoMs proposed in this paper is designed to evaluate a building layout in terms of wireless performance. Instead of any specifically defined radio network, the idea of the approach is to evaluate the maximum achievable performance gain of indoor wireless networks in a building. In this paper, PTP_{\mathrm{T}} [W/m2][\mathrm{W}/\mathrm{m}^{2}] denotes the transmit power density, i.e., transmit power per unit area. For a small cell network with an intensity of Λ[BSs/m2]\Lambda\ \mathrm{[BSs/m^{2}]}, the power of each small cell base station is PT/Λ[W/BS]P_{\mathrm{T}}/\Lambda\ \mathrm{[W/BS]}. To figure out the maximum achievable performance gain, we need to use the best case of indoor wireless network deployment even though it might be not realistic in existing networks. It is predictable that the average SINR for a downlink reference user equipment (UE) with a fixed sensitivity increases with the intensity of transmit elements with a fixed total transmit power. Therefore, an infinite number of transmit elements is assumed in our derivations, i.e., small cells [2] or equivalent antenna elements [39] are assumed to be uniformly distributed in both the indoor and the outdoor environment.

  2. 2.

    Multi-slope path gain models [36] for the open space, the LOS and the NLOS scenarios are employed in this paper. Path gain is defined in the general form as

    G(R)PRPT,\displaystyle G(R)\triangleq\frac{P_{\mathrm{R}}}{P_{\mathrm{T}}}, (1)

    where PRP_{\mathrm{R}} and PTP_{\mathrm{T}} denote the received and the transmit power densities, respectively, and RR denotes the horizontal distance from the transmitter to the receiver.

    The path gain in the open space scenario is given by [37, Eqs. (4,5)]

    GO(R)=min{1,(λ4π)2R2,(hThR)2R4},\displaystyle G_{\mathrm{O}}(R)=\min\left\{1,\left(\frac{\lambda}{4\pi}\right)^{2}R^{-2},\left(h_{\mathrm{T}}h_{\mathrm{R}}\right)^{2}{R^{-4}}\right\}, (2)

    where hTh_{\mathrm{T}} and hRh_{\mathrm{R}} are the heights of the transmit and the receive antenna, respectively; fcf_{\mathrm{c}} is the center frequency, λcfc\lambda\triangleq\frac{c}{f_{\mathrm{c}}} denotes the corresponding wavelength, c=3×108[m/s]c=3\times 10^{8}\mathrm{[m/s]} is the speed of light.

    Path gain models for the LOS and NLOS indoor scenarios can then be computed as

    GL(R)={GO(R),R1m,(λ4π)2RnL,R>1m.\displaystyle G_{\mathrm{L}}(R)=\left\{\begin{array}[]{l l}G_{\mathrm{O}}(R),&R\leq 1~{}\mathrm{m},\\ \left(\frac{\lambda}{4\pi}\right)^{2}R^{-n_{\mathrm{L}}},&R>1~{}\mathrm{m}.\end{array}\right. (5)
    GN(R)={GO(R),R1m,(λ4π)2RnN,R>1m,\displaystyle G_{\mathrm{N}}(R)=\left\{\begin{array}[]{l l}G_{\mathrm{O}}(R),&R\leq 1~{}\mathrm{m},\\ \left(\frac{\lambda}{4\pi}\right)^{2}R^{-n_{\mathrm{N}}},&R>1~{}\mathrm{m},\end{array}\right. (8)

    where the numerical values of nLn_{\mathrm{L}} and nNn_{\mathrm{N}} are determined by the type of environment including building materials, shape and size of rooms, etc. Following the energy conservation law, path gains GO/L/NG_{\mathrm{O/L/N}} have to be less than one.

  3. 3.

    In this paper, we specialize our results to nL=1.73n_{\mathrm{L}}=1.73, and nN=3.19n_{\mathrm{N}}=3.19 following the 3GPP indoor channel model working at 0.5-100 GHz [38, Table 7.4.1-1]. Since the 3GPP model is defined for R>1mR>1~{}\mathrm{m}, we assume that GL(R)=GN(R)=GO(R)G_{\mathrm{L}}(R)=G_{\mathrm{N}}(R)=G_{\mathrm{O}}(R), at R1mR\leq 1~{}\mathrm{m}.

  4. 4.

    We constrain ourselves to Rs>1R_{s}>1 m\rm m, s{O,L,N}s\in\left\{\mathrm{O},\mathrm{L},\mathrm{N}\right\}, and therefore Pth<PT(λ4π)2P_{\mathrm{th}}<P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}.

  5. 5.

    In this paper, the detectable power is defined as the receive power from transmit elements satisfying the condition

    PTGs(R)>Pth,s{O,L,N},\displaystyle P_{\mathrm{T}}G_{s}(R)>P_{\mathrm{th}},s\in\{\mathrm{O,L,N}\}, (9)

    where PthP_{\mathrm{th}} [Wm2][\mathrm{W}\mathrm{m}^{-2}] is a threshold power level determined by the sensitivity of the UE receiver.

  6. 6.

    ROR_{\mathrm{O}}, RLR_{\mathrm{L}} and RNR_{\mathrm{N}} denote the coverage distances under open space, LOS and NLOS scenarios, respectively and satisfy the identities PTGO(RO)=PthP_{\mathrm{T}}G_{\mathrm{O}}(R_{\mathrm{O}})=P_{\mathrm{th}}, PTGL(RL)=PthP_{\mathrm{T}}G_{\mathrm{L}}(R_{\mathrm{L}})=P_{\mathrm{th}}, and PTGN(RN)=PthP_{\mathrm{T}}G_{\mathrm{N}}(R_{\mathrm{N}})=P_{\mathrm{th}}. Hence, it can be shown after some straightforward algebraic manipulations that

    RO={PTPthλ4π,PTPth<(4π)4(hThR)2λ4,(PTPth)14hThR,PTPth(4π)4(hThR)2λ4,\displaystyle R_{\mathrm{O}}=\left\{\begin{array}[]{l l}\sqrt{\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}}\frac{\lambda}{4\pi},&\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}<\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}},\\ \left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)^{\frac{1}{4}}\sqrt{h_{\mathrm{T}}h_{\mathrm{R}}},&\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\geq\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}},\end{array}\right. (12)
    RL=(PTPth)1nL(λ4π)2nL,\displaystyle R_{\mathrm{L}}=\left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)^{\frac{1}{n_{\mathrm{L}}}}\left(\frac{\lambda}{4\pi}\right)^{\frac{2}{n_{\mathrm{L}}}}, (13)
    RN=(PTPth)1nN(λ4π)2nN.\displaystyle R_{\mathrm{N}}=\left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)^{\frac{1}{n_{\mathrm{N}}}}\left(\frac{\lambda}{4\pi}\right)^{\frac{2}{n_{\mathrm{N}}}}. (14)
  7. 7.

    In the best case, the reference UE can make use of all the detectable power, which is constrained by the sensitivity of its receiver. The receive power is in turn assumed to be the result of maximum ratio transmissions (MRT) in distributed massive MIMO networks or coordinated multipoint transmission (CoMP) in ultra-dense small-cell networks [40, 41], where all network elements transmit the same signal for the reference UE continuously at full power. To simplify the analysis, we further assume that all transmit elements use the same frequency band.

  8. 8.

    In cooperative ultra-dense small-cell/massive MIMO systems, channel responses are smoothed by the extremely large spatial diversity as a result of the favorable action of the law of large numbers. In essence, small-scale fading is negligible [42], therefore, only large-scale fading is considered in this paper.

With all the assumptions, we first observe the impact of the room structure on a reference UE located in the center of a circular room with a radius of RWR_{\mathrm{W}}. As shown in Fig. 1, the SINR of the reference UE varies significantly for different values of RWR_{\mathrm{W}}. Thus, it is necessary to consider this impact on the building design stage. In practice, rooms in a building are not circular in general, and the evaluation of the building is therefore much more complicated. In the following Section III, we will design the FoMs to evaluate the wireless friendliness of building layouts.

Refer to caption
Figure 1: Impact of the radii RWR_{\mathrm{W}} of a circular room on the SINR of the reference UE at the center of it. σ=93\sigma=-93 dB, PT=30P_{\mathrm{T}}=-30 dBW/m2, Pth=100P_{\mathrm{th}}=-100 dBW/m2, and hT=hR=1.2h_{\mathrm{T}}=h_{\mathrm{R}}=1.2 m.

III FoM definition

In order to meaningfully evaluate the wireless performance of buildings, the impact of buildings on wireless networks has to be well understood. To help in the analysis we employ the open space scenario with path gain model (2) as a benchmark or reference. Then, the impact of buildings on the wireless networks can be evaluated by comparing the SINR of the signals received indoors with that in the open space scenario.

The expression of the SINR in the open space is given as follows. Using a polar-coordinate system (R,θ)(R,\theta) by choosing the UE as the reference point and an arbitrary angle as the reference 0-angle in the open space, the power of the intended signal is given by the integral of the receive power from all transmit elements with a distance that is less than ROR_{\mathrm{O}}, i.e.,

PO=02π0ROPTGO(R)RdRdθ.\displaystyle P_{\mathrm{O}}=\int_{0}^{2\pi}\int_{0}^{R_{\mathrm{O}}}P_{\mathrm{T}}G_{\mathrm{O}}(R)R\mathrm{d}R\mathrm{d}\theta. (15)

Whereas, the power of the interference signal is given by

IO=02πRO+PTGO(R)RdRdθ.\displaystyle I_{\mathrm{O}}=\int_{0}^{2\pi}\int_{R_{\mathrm{O}}}^{+\infty}P_{\mathrm{T}}G_{\mathrm{O}}(R)R\mathrm{d}R\mathrm{d}\theta. (16)

Thus, the SINR in an open space scenario is given by

γO=POIO+σ2.\displaystyle\gamma_{\mathrm{O}}=\frac{P_{\mathrm{O}}}{I_{\mathrm{O}}+\sigma^{2}}. (17)

Then, the SINR of a UE in the considered building is given as follows. For an arbitrary transmit element located at (R,θ)(R,\theta), the transmission from the element to the UE could be either LOS or NLOS. If the transmission is LOS, the received signal is considered as an intended signal if R<RLR<R_{\mathrm{L}}. Otherwise, if RRLR\geq R_{\mathrm{L}}, the received signal is considered as an interference signal. Similarly, if the transmission is NLOS, the intended signal satisfies R<RNR<R_{\mathrm{N}}, and the interference signal satisfies RRNR\geq R_{\mathrm{N}}. For a given (R,θ)(R,\theta), we use a Boolean variable 𝒱(R,θ)\mathcal{V}(R,\theta) to denote that “the transmission from the transmit element at (R,θ)(R,\theta) to the UE is LOS”. Accordingly, ¬𝒱(R,θ)\lnot\mathcal{V}(R,\theta) denotes that “the transmission from the transmit element at (R,θ)(R,\theta) to the UE is NLOS”. If a signal transmitted from (R,θ)(R,\theta) is an intended signal, we have [𝒱(R,θ)(R<RL)][¬𝒱(R,θ)(R<RN)][\mathcal{V}(R,\theta)\land(R<R_{L})]\lor[\lnot\mathcal{V}(R,\theta)\land(R<R_{N})]. Then, the power of intended signals in the considered building is given by

PB=𝒱(R,θ)(R<RL)PTGL(R)RdRdθPL+¬𝒱(R,θ)(R<RN)PTGN(R)RdRdθPN,\displaystyle P_{\mathrm{B}}=\underbrace{\int_{\mathcal{V}(R,\theta)\land(R<R_{L})}P_{\mathrm{T}}G_{\mathrm{L}}(R)R\mathrm{d}R\mathrm{d}\theta}_{P_{\mathrm{L}}}+\underbrace{\int_{\lnot\mathcal{V}(R,\theta)\land(R<R_{N})}P_{\mathrm{T}}G_{\mathrm{N}}(R)R\mathrm{d}R\mathrm{d}\theta}_{P_{\mathrm{N}}}, (18)

where PLP_{\mathrm{L}} and PNP_{\mathrm{N}}, respectively, denote receive power of intended signals from LOS and NLOS small cells/antenna elements, and can be computed separately since 𝒱(R,θ)(R<RL)\mathcal{V}(R,\theta)\land(R<R_{L}) and ¬𝒱(R,θ)(R<RN)\lnot\mathcal{V}(R,\theta)\land(R<R_{N}) are exclusive. Similarly, the power of the interference signal is given by

IB=𝒱(R,θ)(RRL)PTGL(R)RdRdθIL+¬𝒱(R,θ)(RRN)PTGN(R)RdRdθIN,\displaystyle I_{\mathrm{B}}=\underbrace{\int_{\mathcal{V}(R,\theta)\land(R\geq R_{L})}P_{\mathrm{T}}G_{\mathrm{L}}(R)R\mathrm{d}R\mathrm{d}\theta}_{I_{\mathrm{L}}}+\underbrace{\int_{\lnot\mathcal{V}(R,\theta)\land(R\geq R_{N})}P_{\mathrm{T}}G_{\mathrm{N}}(R)R\mathrm{d}R\mathrm{d}\theta}_{I_{\mathrm{N}}}, (19)

where ILI_{\mathrm{L}} and INI_{\mathrm{N}}, respectively, denote the LOS and the NLOS received powers of interference signals. Therefore, we have the expression of the indoor SINR as

γB=PBIB+σ2=PL+PNIL+IN+σ2.\displaystyle\gamma_{\mathrm{B}}=\frac{P_{\mathrm{B}}}{I_{\mathrm{B}}+\sigma^{2}}=\frac{P_{\mathrm{L}}+P_{\mathrm{N}}}{I_{\mathrm{L}}+I_{\mathrm{N}}+\sigma^{2}}. (20)

Combining (17) and (20), we obtain

γB=IO+σ2IL+IN+σ2PL+PNPOγO=gIgPγO,\displaystyle\gamma_{\mathrm{B}}=\frac{{I_{\mathrm{O}}+\sigma^{2}}}{I_{\mathrm{L}}+I_{\mathrm{N}}+\sigma^{2}}\frac{P_{\mathrm{L}}+P_{\mathrm{N}}}{P_{\mathrm{O}}}\gamma_{\mathrm{O}}=g_{\mathrm{I}}g_{\mathrm{P}}\gamma_{\mathrm{O}}, (21)

and we define two FoMs of a building’s wireless performance as [15, 16]

gIIO+σ2IL+IN+σ2,\displaystyle g_{\mathrm{I}}\triangleq\frac{{I_{\mathrm{O}}+\sigma^{2}}}{I_{\mathrm{L}}+I_{\mathrm{N}}+\sigma^{2}}, (22)
gPPL+PNPO,\displaystyle g_{\mathrm{P}}\triangleq\frac{P_{\mathrm{L}}+P_{\mathrm{N}}}{P_{\mathrm{O}}}, (23)

where gIg_{\mathrm{I}} and gPg_{\mathrm{P}} denote the IG and the PG, respectively. The variables IOI_{\mathrm{O}} and POP_{\mathrm{O}} denote the power of interference signals and the power of intended signals in the open space scenario, respectively. In the LOS scenario, ILI_{\mathrm{L}} and PLP_{\mathrm{L}} denote the power of interference signals and the power of intended signals, respectively. In the NLOS scenario, INI_{\mathrm{N}} and PNP_{\mathrm{N}} denote the power of interference signals and the power of intended signals from the transmit elements in the floor where the UE located, respectively. σ2\sigma^{2} denotes the power of thermal noise.

Remark 1.

The FoMs gIg_{\mathrm{I}} and gPg_{\mathrm{P}} are intrinsic wireless performance characteristics of buildings. They represent the effective change in the SINR when coverage is provided within the building relative to the SINR in the open space scenario. On one hand, gIg_{\mathrm{I}} captures the impact of the building on the power of interference signals due to blockage, while, on the other hand, gPg_{\mathrm{P}} captures the impact of the building on the power of intended signals.

The proposed FoMs and their analytic computation approach, will provide civil engineers and architects guidance to evaluate the wireless friendliness of a building layout under their design. If FoMs are lower than expected, an acceptable wireless performance will not be achievable. For communication engineers, the FoM is a performance gain upper bound provided by the buildings under design. To achieve the upper bound, Assumptions 1) and 7) stated in Section II need to be approached.

In our model we consider an infinite number of small cells/antenna elements. Therefore, analytic expressions are comparatively more computationally efficient than performing Monte Carlo simulations to approximately compute the FoMs. Also, analytic results can provide accurate results and lead to a more straightforward comparison and optimization of the building layout. Therefore, we propose next an approach to derive gIg_{\mathrm{I}} and gPg_{\mathrm{P}} in closed form for evaluating the wireless performance of buildings.

IV Analytic computation of FoM

Next we provide analytic equations for intended and interference signals starting with the open space scenario (IOI_{\mathrm{O}} and POP_{\mathrm{O}}) and later for the LOS and NLOS scenarios too (i.e., ILI_{\mathrm{L}}, PLP_{\mathrm{L}}, INI_{\mathrm{N}}, and PNP_{\mathrm{N}}, respectively).

Theorem 1.

In the open space scenario, the power of interference signals IOI_{\mathrm{O}} and the power of intended signals POP_{\mathrm{O}} are computed by (26) and (29), respectively.

IO={PTλ28π{12+ln(16Pthπ2hThRPTλ2)},PTPth<(4π)4(hThR)2λ4,πPthPThThR,PTPth(4π)4(hThR)2λ4.\displaystyle I_{\mathrm{O}}=\left\{\begin{array}[]{l l}\frac{P_{\mathrm{T}}\lambda^{2}}{8\pi}\left\{\frac{1}{2}+\ln\left(\frac{16\sqrt{P_{\mathrm{th}}}\pi^{2}h_{\mathrm{T}}h_{\mathrm{R}}}{\sqrt{P_{\mathrm{T}}}\lambda^{2}}\right)\right\},&\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}<\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}},\\ \pi\sqrt{P_{\mathrm{th}}P_{\mathrm{T}}}h_{\mathrm{T}}h_{\mathrm{R}},&\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\geq\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}}.\end{array}\right. (26)
PO={PTλ216π[1+ln(PTPth)],PTPth<(4π)4(hThR)2λ4,PTλ216π+PTλ28πln(16π2hThRλ2)+λ4PT29π3(hThR)2πPth2,PTPth(4π)4(hThR)2λ4.\displaystyle P_{\mathrm{O}}=\left\{\begin{array}[]{l l}\frac{P_{\mathrm{T}}\lambda^{2}}{16\pi}\left[1+\ln\left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)\right],&\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}<\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}},\\ \frac{P_{\mathrm{T}}\lambda^{2}}{16\pi}+\frac{P_{\mathrm{T}}\lambda^{2}}{8\pi}\ln\left(\frac{16\pi^{2}h_{\mathrm{T}}h_{\mathrm{R}}}{\lambda^{2}}\right)+\frac{\lambda^{4}{P_{\mathrm{T}}}}{2^{9}\pi^{3}(h_{\mathrm{T}}h_{\mathrm{R}})^{2}}-\frac{\pi P_{\mathrm{th}}}{2},&\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\geq\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}}.\end{array}\right. (29)
Proof:

See Appendix A. ∎

The parameters of the equations have been defined above in Section II and are omitted here for the sake of compactness.

Refer to caption
Figure 2: An example of the toy model and its ΩPL\Omega_{\mathrm{PL}}, ΩPN\Omega_{\mathrm{PN}}, ΩIL\Omega_{\mathrm{IL}}, and ΩIN\Omega_{\mathrm{IN}}.
Refer to caption
Figure 3: Examples of dividing a room into toy models.

The computation of ILI_{\mathrm{L}}, PLP_{\mathrm{L}}, INI_{\mathrm{N}}, and PNP_{\mathrm{N}} is not straightforward since the indoor scenarios are not homogeneous. Therefore, in order to facilitate their derivation, we propose a toy model (TM) as shown in Fig. 3. In the TM, we consider a circular sector with an infinite radius. A reference (or probing) UE is located at the vertex of the sector. A wall of finite length is modeled as a straight line segment across the sector with endpoints lying on the radii of the circular sector. The distance from the UE to the wall is denoted as D0D_{0}. A polar-coordinate system (R,θ)(R,\theta) is introduced by choosing the UE as the reference point. The perpendicular direction to the wall is assumed to be the reference 0-angle. θl\theta_{\mathrm{l}} and θr\theta_{\mathrm{r}} denote angles of the left and right end points of the wall, respectively. With these definitions, π2<θl<θr<π2-\frac{\pi}{2}<\theta_{\mathrm{l}}<\theta_{\mathrm{r}}<\frac{\pi}{2}.

It is worthwhile to note that the proposed TM of a wall can be used to represent any polygonal-shaped room. Indeed, any room can be modelled by adding several TMs, i.e., one for each wall that can be seen by the reference UE in the room. Some examples are given in Fig. 3. For a specific direction, whether an area is LOS or NLOS is determined by the closest wall.

For a room consisting of NTMN_{\mathrm{TM}} TMs, we then have

IL=nTM=1NTMIL,nTM,\displaystyle I_{\mathrm{L}}=\sum\limits_{n_{\mathrm{TM}}=1}^{N_{\mathrm{TM}}}I_{\mathrm{L},n_{\mathrm{TM}}}, (30)
PL=nTM=1NTMPL,nTM,\displaystyle P_{\mathrm{L}}=\sum\limits_{n_{\mathrm{TM}}=1}^{N_{\mathrm{TM}}}P_{\mathrm{L},n_{\mathrm{TM}}}, (31)
IN=nTM=1NTMIN,nTM,\displaystyle I_{\mathrm{N}}=\sum\limits_{n_{\mathrm{TM}}=1}^{N_{\mathrm{TM}}}I_{\mathrm{N},n_{\mathrm{TM}}}, (32)
PN=nTM=1NTMPN,nTM.\displaystyle P_{\mathrm{N}}=\sum\limits_{n_{\mathrm{TM}}=1}^{N_{\mathrm{TM}}}P_{\mathrm{N},n_{\mathrm{TM}}}. (33)

For known closed-form expressions of IL,nTMI_{\mathrm{L},n_{\mathrm{TM}}}, PL,nTMP_{\mathrm{L},n_{\mathrm{TM}}}, IN,nTMI_{\mathrm{N},n_{\mathrm{TM}}}, and PN,nTMP_{\mathrm{N},n_{\mathrm{TM}}}, it is straightforward to compute ILI_{\mathrm{L}}, PLP_{\mathrm{L}}, INI_{\mathrm{N}}, and PNP_{\mathrm{N}} by (30-33). For an arbitrary nTMn_{\mathrm{TM}} (i.e., TM of a wall contribution to the computed signals) we have

PL,nTM=ΩPLPTGL(R)dΩ,\displaystyle P_{\mathrm{L},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{PL}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}\mathrm{d}\Omega, (34)
IL,nTM=ΩILPTGL(R)dΩ,\displaystyle I_{\mathrm{L},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{IL}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}\mathrm{d}\Omega, (35)
PN,nTM=ΩPNPTGN(R)dΩ,\displaystyle P_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{PN}}}{P_{\mathrm{T}}}{G_{\mathrm{N}}(R)}\mathrm{d}\Omega, (36)
IN,nTM=ΩINPTGN(R)dΩ,\displaystyle I_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{IN}}}{P_{\mathrm{T}}}{G_{\mathrm{N}}(R)}\mathrm{d}\Omega, (37)

where ΩPL\Omega_{\mathrm{PL}}, ΩIL\Omega_{\mathrm{IL}}, ΩPN\Omega_{\mathrm{PN}}, and ΩIN\Omega_{\mathrm{IN}} denote the areas contributing to the LOS intended signals, the LOS interference signals, the NLOS intended signals, and the NLOS interference signals at the reference UE, respectively. Expressed in polar-coordinates, the straight line equation of the wall is given by R=D0cosθR=\frac{D_{0}}{\cos{\theta}}. For any infinitesimal area at position (R,θ)(R,\theta), the area is a LOS area if R<D0cosθ.R<\frac{D_{0}}{\cos{\theta}}. Else, the area is NLOS if RD0cosθ.R\geq\frac{D_{0}}{\cos{\theta}}. Therefore, ΩIL\Omega_{\mathrm{IL}}, ΩPL\Omega_{\mathrm{PL}}, ΩIN\Omega_{\mathrm{IN}}, and ΩPN\Omega_{\mathrm{PN}} are, respectively, defined as

ΩPL{(R,θ)|(R<min{RL,D0cosθ})(π2<θl<θ<θr<π2)},\displaystyle\Omega_{\mathrm{PL}}\triangleq\left\{\left(R,\theta\right)\left|\left(R<\min\left\{R_{\mathrm{L}},\frac{D_{0}}{\cos{\theta}}\right\}\right)\land\left(-\frac{\pi}{2}<\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}<\frac{\pi}{2}\right)\right.\right\}, (38)
ΩIL{(R,θ)|(RLR<D0cosθ)(π2<θl<θ<θr<π2)},\displaystyle\Omega_{\mathrm{IL}}\triangleq\left\{\left(R,\theta\right)\left|\left(R_{\mathrm{L}}\leq R<\frac{D_{0}}{\cos{\theta}}\right)\land\left(-\frac{\pi}{2}<\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}<\frac{\pi}{2}\right)\right.\right\}, (39)
ΩPN{(R,θ)|(D0cosθR<RN)(π2<θl<θ<θr<π2)},\displaystyle\Omega_{\mathrm{PN}}\triangleq\left\{\left(R,\theta\right)\left|\left(\frac{D_{0}}{\cos{\theta}}\leq R<R_{\mathrm{N}}\right)\land\left(-\frac{\pi}{2}<\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}<\frac{\pi}{2}\right)\right.\right\}, (40)
ΩIN{(R,θ)|(Rmax{RN,D0cosθ})(π2<θl<θ<θr<π2)}.\displaystyle\Omega_{\mathrm{IN}}\triangleq\left\{\left(R,\theta\right)\left|\left(R\geq\max\left\{R_{\mathrm{N}},\frac{D_{0}}{\cos{\theta}}\right\}\right)\land\left(-\frac{\pi}{2}<\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}<\frac{\pi}{2}\right)\right.\right\}. (41)

In a TM with a given set of (PTP_{\mathrm{T}}, PthP_{\mathrm{th}}, θr\theta_{\mathrm{r}}, θl\theta_{\mathrm{l}}, D0D_{0}, λ\lambda), the integrals PL,nTMP_{\mathrm{L},n_{\mathrm{TM}}}, IL,nTMI_{\mathrm{L},n_{\mathrm{TM}}}, PN,nTMP_{\mathrm{N},n_{\mathrm{TM}}} and IN,nTMI_{\mathrm{N},n_{\mathrm{TM}}} are computed by following Theorems 2-5.

Theorem 2.

For an arbitrary nTMn_{\mathrm{TM}}, PL,nTMP_{\mathrm{L},n_{\mathrm{TM}}} is computed by TABLE I in closed form, where PL,1P_{\mathrm{L},1}, PL,2P_{\mathrm{L},2}, and PL,3P_{\mathrm{L},3} are given by the following expressions

PL,1=PT(λ4π)2(θrθl)[12ln(λ4π)]+PT(λ4π)2Z0(θl,θr,1,RL,nL1),\displaystyle P_{\mathrm{L},1}=P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}\left(\theta_{\mathrm{r}}-\theta_{\mathrm{l}}\right)\left[\frac{1}{2}-\ln\left(\frac{\lambda}{4\pi}\right)\right]+P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{0}(\theta_{\mathrm{l}},\theta_{\mathrm{r}},1,R_{\mathrm{L}},n_{\mathrm{L}}-1), (42)
PL,2=PT(λ4π)2Z1(max{θl,θL1},min{θr,θL1},RL,D0,nL1),\displaystyle P_{\mathrm{L},2}=P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\max\{\theta_{\mathrm{l}},-\theta_{\mathrm{L1}}\},\min\{\theta_{\mathrm{r}},\theta_{\mathrm{L1}}\},R_{\mathrm{L}},D_{0},n_{\mathrm{L}}-1\right), (43)
PL,3=PT(λ4π)2Z1(max{θl,θL2},min{θr,θL2},1,D0,1)\displaystyle P_{\mathrm{L},3}=P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\max\{\theta_{\mathrm{l}},-\theta_{\mathrm{L2}}\},\min\{\theta_{\mathrm{r}},\theta_{\mathrm{L2}}\},1,D_{0},1\right)
PT(λ4π)2Z1(max{θl,θL2},min{θr,θL2},1,D0,nL1),\displaystyle\hskip 32.52127pt-P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\max\{\theta_{\mathrm{l}},-\theta_{\mathrm{L2}}\},\min\{\theta_{\mathrm{r}},\theta_{\mathrm{L2}}\},1,D_{0},n_{\mathrm{L}}-1\right), (44)
θL1=arccos(D0RL),\displaystyle\theta_{\mathrm{L1}}=\arccos\left(\frac{D_{0}}{R_{\mathrm{L}}}\right), (45)
θL2=arccos(D0),\displaystyle\theta_{\mathrm{L2}}=\arccos\left(D_{0}\right), (46)
Z0(z1,z2,z3,z4,z5)={(z2z1)ln(z4z3),z5=1,(z2z1)(z41z5z31z5)1z5,z51,\displaystyle Z_{0}(z_{1},z_{2},z_{3},z_{4},z_{5})=\left\{\begin{array}[]{ll}(z_{2}-z_{1})\ln\left(\frac{z_{4}}{z_{3}}\right),&z_{5}=1,\\ \frac{(z_{2}-z_{1})(z_{4}^{1-z_{5}}-z_{3}^{1-z_{5}})}{1-z_{5}},&z_{5}\neq 1,\end{array}\right. (49)
Z1(z1,z2,z3,z4,z5)={z42[tan(z2)tan(z1)]2+(z1z2)z322,z5=1,z4ln(tan(z2)+sec(z2)tan(z1)+sec(z1))+(z1z2)z3,z5=0,(z2z1)ln(2z4z3)+Im[Li2(e2jz2)Li2(e2jz1)]2,z5=1,{B(12,z52)[sgn(z2)sgn(z1)]z41z52(1z5)+(z1z2)z31z51z5+z41z5(1z5)z5×[sgn(z1)cosz5(z1)2F1(z52,12;z5+22,cos2(z1))sgn(z2)cosz5(z2)2F1(z52,12;z5+22,cos2(z2))],else,\displaystyle Z_{1}(z_{1},z_{2},z_{3},z_{4},z_{5})=\left\{\begin{array}[]{ll}\frac{z_{4}^{2}[\tan(z_{2})-\tan(z_{1})]}{2}+\frac{(z_{1}-z_{2})z_{3}^{2}}{2},&z_{5}=-1,\\ z_{4}\ln\left(\frac{\tan(z_{2})+\sec(z_{2})}{\tan(z_{1})+\sec(z_{1})}\right)+(z_{1}-z_{2})z_{3},&z_{5}=0,\\ (z_{2}-z_{1})\ln\left(\frac{2z_{4}}{z_{3}}\right)+\frac{\mathrm{Im}\left[\mathrm{Li}_{2}(-e^{2jz_{2}})-\mathrm{Li}_{2}(-e^{2jz_{1}})\right]}{2},&z_{5}=1,\\ \left\{\begin{array}[]{l}\frac{B\left(\frac{1}{2},\frac{z_{5}}{2}\right)\left[\mathrm{sgn}(z_{2})-\mathrm{sgn}(z_{1})\right]z_{4}^{1-z_{5}}}{2(1-z_{5})}\\ +\frac{(z_{1}-z_{2})z_{3}^{1-z_{5}}}{1-z_{5}}+\frac{z_{4}^{1-z_{5}}}{(1-z_{5})z_{5}}\times\\ \left[\mathrm{sgn}(z_{1})\cos^{z_{5}}(z_{1})\ _{2}F_{1}\left(\frac{z_{5}}{2},\frac{1}{2};\frac{z_{5}+2}{2},\cos^{2}(z_{1})\right)\right.\\ \left.-\mathrm{sgn}(z_{2})\cos^{z_{5}}(z_{2})\ _{2}F_{1}\left(\frac{z_{5}}{2},\frac{1}{2};\frac{z_{5}+2}{2},\cos^{2}(z_{2})\right)\right],\end{array}\right.&\mathrm{else},\end{array}\right. (58)

where π2<z1<z2<π2-\frac{\pi}{2}<z_{1}<z_{2}<\frac{\pi}{2}, z3>0z_{3}>0, z4>0z_{4}>0, and z5>1z_{5}>-1.

Proof:

See Appendix B. ∎

TABLE I: Computation of PL,nTMP_{\mathrm{L},n_{\mathrm{TM}}}, where PL,1P_{\mathrm{L},1}, PL,2P_{\mathrm{L},2}, and PL,3P_{\mathrm{L},3} are respectively computed by (42), (43), and (2).
Case D0D_{0} θ\theta PL,nTMP_{\mathrm{L},n_{\mathrm{TM}}}
1 1D0<RL1\leq D_{0}<R_{\mathrm{L}} (θl<θL1)(θr>θL1)(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L1}}) PL,1+PL,2P_{\mathrm{L},1}+P_{\mathrm{L},2}
2 λ4π<D0<1\frac{\lambda}{4\pi}<D_{0}<1 (θl<θL2)(θr>θL2)(\theta_{\mathrm{l}}<\theta_{\mathrm{L2}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L2}}) PL,1+PL,2+PL,3P_{\mathrm{L},1}+P_{\mathrm{L},2}+P_{\mathrm{L},3}
3 λ4π<D0<1\frac{\lambda}{4\pi}<D_{0}<1 (θl>θL2)(θr<θL2)[(θl<θL1)(θr>θL1)](\theta_{\mathrm{l}}>\theta_{\mathrm{L2}})\lor(\theta_{\mathrm{r}}<-\theta_{\mathrm{L2}})\land\left[(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L1}})\right] PL,1+PL,2P_{\mathrm{L},1}+P_{\mathrm{L},2}
4 else PL,1P_{\mathrm{L},1}
Theorem 3.

IL,nTMI_{\mathrm{L},n_{\mathrm{TM}}} is computed by TABLE II in closed form, where

IL,1=PT(λ4π)2Z1(θl,θr,RL,D0,nL1),\displaystyle I_{\mathrm{L},1}=P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\theta_{\mathrm{l}},\theta_{\mathrm{r}},R_{\mathrm{L}},D_{0},n_{\mathrm{L}}-1\right), (59)
IL,2=PT(λ4π)2Z1(θl,θL1,RL,D0,nL1),\displaystyle I_{\mathrm{L},2}=P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\theta_{\mathrm{l}},-\theta_{\mathrm{L1}},R_{\mathrm{L}},D_{0},n_{\mathrm{L}}-1\right), (60)
IL,3=PT(λ4π)2Z1(θL1,θr,RL,D0,nL1),\displaystyle I_{\mathrm{L},3}=P_{T}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\theta_{\mathrm{L1}},\theta_{\mathrm{r}},R_{\mathrm{L}},D_{0},n_{\mathrm{L}}-1\right), (61)

where θL1\theta_{\mathrm{L1}} and Z1Z_{1} are computed by (45) and (107), respectively.

Proof:

See Appendix C. ∎

TABLE II: Computation of IL,nTMI_{\mathrm{L},n_{\mathrm{TM}}}, where IL,1I_{\mathrm{L},1}, IL,2I_{\mathrm{L},2}, and IL,3I_{\mathrm{L},3} are respectively computed by (59), (60), and (61).
Case D0D_{0} θ\theta IL,nTMI_{\mathrm{L},n_{\mathrm{TM}}}
1 RLD0R_{\mathrm{L}}\leq D_{0} θl<θr\theta_{\mathrm{l}}<\theta_{\mathrm{r}} IL,1I_{\mathrm{L},1}
2 RL>D0R_{\mathrm{L}}>D_{0} (θr<θL1)(θl>θL1)\left(\theta_{\mathrm{r}}<-\theta_{\mathrm{L1}}\right)\lor\left(\theta_{\mathrm{l}}>\theta_{\mathrm{L1}}\right) IL,1I_{\mathrm{L},1}
3 RL>D0R_{\mathrm{L}}>D_{0} (θl<θL1)(θL1<θr)\left(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}}\right)\land\left(-\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right) IL,2δ(θl<θL1)+IL,3δ(θl>θL1)\begin{array}[]{l}I_{\mathrm{L},2}\delta(\theta_{\mathrm{l}}<-\theta_{\mathrm{L1}})\\ +I_{\mathrm{L},3}\delta(\theta_{\mathrm{l}}>\theta_{\mathrm{L1}})\end{array}
4 else 0
Theorem 4.

PN,nTMP_{\mathrm{N},n_{\mathrm{TM}}} is computed by TABLE III in closed form, where PN,1P_{\mathrm{N},1} and PN,2P_{\mathrm{N},2} are respectively computed by

PN,1=PT(λ4π)2Z1(max{θl,θN1},min{θr,θN1},RN,D0,nN1),\displaystyle P_{\mathrm{N},1}=-P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\max\{\theta_{\mathrm{l}},-\theta_{\mathrm{N1}}\},\min\{\theta_{\mathrm{r}},\theta_{\mathrm{N1}}\},R_{\mathrm{N}},D_{0},n_{\mathrm{N}}-1\right), (62)
PN,2=PT(λ4π)2Z1(max{θl,θN2},min{θr,θN2},1,D0,1)\displaystyle P_{\mathrm{N},2}=-P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\max\{\theta_{\mathrm{l}},-\theta_{\mathrm{N2}}\},\min\{\theta_{\mathrm{r}},\theta_{\mathrm{N2}}\},1,D_{0},1\right)
+PT(λ4π)2Z1(max{θl,θN2},min{θr,θN2},1,D0,nN1),\displaystyle\hskip 32.52127pt+P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\max\{\theta_{\mathrm{l}},-\theta_{\mathrm{N2}}\},\min\{\theta_{\mathrm{r}},\theta_{\mathrm{N2}}\},1,D_{0},n_{\mathrm{N}}-1\right), (63)
θN1=arccos(D0RN),\displaystyle\theta_{\mathrm{N1}}=\arccos\left(\frac{D_{0}}{R_{\mathrm{N}}}\right), (64)
θN2=arccos(D0),\displaystyle\theta_{\mathrm{N2}}=\arccos\left(D_{0}\right), (65)

where Z0Z_{0} and Z1Z_{1} are computed by (87) and (107), respectively.

Proof:

See Appendix D. ∎

TABLE III: Computation of PN,nTMP_{\mathrm{N},n_{\mathrm{TM}}}, where PN,1P_{\mathrm{N},1} and PN,2P_{\mathrm{N},2} are computed by (62) and (4), respectively.
Case D0D_{0} θ\theta PN,nTMP_{\mathrm{N},n_{\mathrm{TM}}}
1 1D0<RN1\leq D_{0}<R_{\mathrm{N}} (θl<θN1)(θr>θN1)(\theta_{\mathrm{l}}<\theta_{\mathrm{N1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{N1}}) PN,1P_{\mathrm{N},1}
2 λ4π<D0<1\frac{\lambda}{4\pi}<D_{0}<1 (θl<θN2)(θr>θN2)(\theta_{\mathrm{l}}<\theta_{\mathrm{N2}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{N2}}) PN,1+PN,2P_{\mathrm{N},1}+P_{\mathrm{N},2}
3 λ4π<D0<1\frac{\lambda}{4\pi}<D_{0}<1 (θl>θN2)(θr<θN2)[(θl<θN1)(θr>θN1)](\theta_{\mathrm{l}}>\theta_{\mathrm{N2}})\lor(\theta_{\mathrm{r}}<-\theta_{\mathrm{N2}})\land\left[(\theta_{\mathrm{l}}<\theta_{\mathrm{N1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{N1}})\right] PN,1P_{\mathrm{N},1}
4 else 0
Theorem 5.

IN,nTMI_{\mathrm{N},n_{\mathrm{TM}}} is computed by TABLE IV in closed form, where

IN,1=PT(λ4π)2Z1(θl,θr,+,D0,nN1),\displaystyle I_{\mathrm{N},1}=-P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\theta_{\mathrm{l}},\theta_{\mathrm{r}},+\infty,D_{0},n_{\mathrm{N}}-1\right), (66)
IN,2=PT(λ4π)2Z1(max{θl,θN1},min{θr,θN1},RN,D0,nN1),\displaystyle I_{\mathrm{N},2}=P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}Z_{1}\left(\max\{\theta_{\mathrm{l}},-\theta_{\mathrm{N1}}\},\min\{\theta_{\mathrm{r}},\theta_{\mathrm{N1}}\},R_{\mathrm{N}},D_{0},n_{\mathrm{N}}-1\right), (67)

where θN1\theta_{\mathrm{N1}}, and Z1Z_{1} are computed by (64) and (107), respectively.

Proof:

See Appendix E. ∎

TABLE IV: Computation of IN,nTMI_{\mathrm{N},n_{\mathrm{TM}}}, where IN,1I_{\mathrm{N},1} and IN,2I_{\mathrm{N},2} are computed by (66), and (67), respectively.
Case D0D_{0} θ\theta IN,nTMI_{\mathrm{N},n_{\mathrm{TM}}}
1 RND0R_{\mathrm{N}}\leq D_{0} θl<θr\theta_{\mathrm{l}}<\theta_{\mathrm{r}} IN,1I_{\mathrm{N},1}
2 RN>D0R_{\mathrm{N}}>D_{0} (θr<θN1)(θl>θN1)\left(\theta_{\mathrm{r}}<-\theta_{\mathrm{N1}}\right)\lor\left(\theta_{\mathrm{l}}>\theta_{\mathrm{N1}}\right) IN,1I_{\mathrm{N},1}
3 RN>D0R_{\mathrm{N}}>D_{0} (θl<θN1)(θN1<θr)\left(\theta_{\mathrm{l}}<\theta_{\mathrm{N1}}\right)\land\left(-\theta_{\mathrm{N1}}<\theta_{\mathrm{r}}\right) IN,1+IN,2I_{\mathrm{N},1}+I_{\mathrm{N},2}
Refer to caption
Figure 4: The procedure of computing gIg_{\mathrm{I}} and gPg_{\mathrm{P}}.
Remark 2.

The computation of the gIg_{\mathrm{I}} and gPg_{\mathrm{P}} follows the flow chart shown in Fig. 4. The input parameters are the building layout, the reference UE position, nLn_{\mathrm{L}}, nNn_{\mathrm{N}}, fcf_{c}, PTP_{\mathrm{T}} and PthP_{\mathrm{th}}. First, Theorems 2, 3, 4, and 5 are used to compute signals corresponding to the open space, LOS and NLOS scenarios. Second, the obtained (30-33) and (26-29) are substituted into (22-23) to obtain the FoM.

V Comparison theoretical results with Monte-Carlo simulations

In this section, we compare the results obtained by the derived closed-form expressions and the Monte Carlo simulations. The analytic results used in the following Figs. 5-11 are summarized in TABLE V. In Sections V-VI, we use hT=hR=1.2h_{\mathrm{T}}=h_{\mathrm{R}}=1.2 m without further specification.

TABLE V: Analytic results used in Figs. 5-11
Figures associating analytical results.
Fig. 5 (26) and (29)
Figs. 7-7 Theorems 2-5
Figs. 11-11 Remark 2

In Fig. 5, we compare the results obtained by the derived closed-form expressions of POP_{\mathrm{O}} and IOI_{\mathrm{O}}, i.e., (26) and (29), with Monte Carlo simulations. The numerical results show an excellent agreement between the analytic, i.e., closed-form, results and the Monte Carlo simulations.

Refer to caption
(a) fc=1f_{c}=1 GHz
Refer to caption
(b) fc=28f_{c}=28 GHz
Figure 5: Comparison of analytic POP_{\mathrm{O}} and IOI_{\mathrm{O}}, computed by Theorem 1, with Monte-Carlo simulations. PT=30dBWm2P_{\mathrm{T}}=-30\ \mathrm{dBWm^{-2}}.
Refer to caption
(a) θr=0.4\theta_{\mathrm{r}}=-0.4
Refer to caption
(b) θr=0\theta_{\mathrm{r}}=0
Refer to caption
(c) θr=0.4\theta_{\mathrm{r}}=0.4
Refer to caption
(d) θr=1\theta_{\mathrm{r}}=1
Figure 6: Comparison of analytic PLP_{\mathrm{L}}, PNP_{\mathrm{N}}, ILI_{\mathrm{L}} and INI_{\mathrm{N}}, computed by Theorems 2-5, with Monte-Carlo simulations. fc=1GHzf_{c}=1\mathrm{GHz} PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=75P_{\mathrm{th}}=-75 dBWm2\mathrm{dBWm^{-2}}, and θl=1\theta_{\mathrm{l}}=-1. Markers are generated from Monte Carlo simulations. Solid lines are analytic results.
Refer to caption
(a) θr=0.4\theta_{\mathrm{r}}=-0.4
Refer to caption
(b) θr=0\theta_{\mathrm{r}}=0
Refer to caption
(c) θr=0.4\theta_{\mathrm{r}}=0.4
Refer to caption
(d) θr=1\theta_{\mathrm{r}}=1
Figure 7: Comparison of analytic PLP_{\mathrm{L}}, PNP_{\mathrm{N}}, ILI_{\mathrm{L}} and INI_{\mathrm{N}}, computed by Theorems 2-5, with Monte-Carlo simulations. fc=28GHzf_{c}=28\mathrm{GHz} PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}, and θl=1\theta_{\mathrm{l}}=-1 rad. Markers are generated from Monte Carlo simulations. Solid lines are analytic results.

Let’s now consider two generic wireless systems operating at the 11 GHz and the 2828 GHz bands, respectively. The parameters PLP_{\mathrm{L}}, ILI_{\mathrm{L}}, PNP_{\mathrm{N}}, and INI_{\mathrm{N}} are validated in Figs. 7-7, respectively, where analytic PLP_{\mathrm{L}}, ILI_{\mathrm{L}}, PNP_{\mathrm{N}}, and INI_{\mathrm{N}} are provided in TABLEs I-IV. The analytic results match Monte Carlo simulations very well. With an increasing D0D_{0}, more transmitted power is considered as LOS power, and therefore PLP_{\mathrm{L}} and ILI_{\mathrm{L}} increase, while PNP_{\mathrm{N}} and INI_{\mathrm{N}} decrease.

Following Remark 2, the analytic gIg_{\mathrm{I}} and gPg_{\mathrm{P}} of a rectangular room of size 5m×10m5\ \mathrm{m}\times 10\ \mathrm{m} are validated in Figs. 11 and 11 for the 1 GHz and the 28 GHz bands, respectively. Also, gIg_{\mathrm{I}} and gPg_{\mathrm{P}} of an L-shaped corner room are validated in Figs. 11 and 11. Floor plans of the considered rooms are given in Fig. 3. Numerical results show that the analytic gIg_{\mathrm{I}} and gPg_{\mathrm{P}} agree well with the Monte Carlo simulations.

Refer to caption
(a) gIg_{\mathrm{I}}
Refer to caption
(b) gPg_{\mathrm{P}}
Figure 8: Comparison of analytic gains gIg_{\mathrm{I}} and gPg_{\mathrm{P}} with Monte-Carlo simulations at 1 GHz in a 5m×10m5\mathrm{m}\times 10\mathrm{m} rectangular room. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=90P_{\mathrm{th}}=90 dBWm2\mathrm{dBWm^{-2}}. Markers are generated from Monte Carlo simulations. Solid lines are computed by Remark 2.
Refer to caption
(a) gIg_{\mathrm{I}}
Refer to caption
(b) gPg_{\mathrm{P}}
Figure 9: Comparison of analytic gains gIg_{\mathrm{I}} and gPg_{\mathrm{P}} with Monte-Carlo simulations at 28 GHz in a 5m×10m5\mathrm{m}\times 10\mathrm{m} rectangular room. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}. Markers are generated from Monte Carlo simulations. Solid lines are computed by Remark 2.
Refer to caption
(a) gIg_{\mathrm{I}}
Refer to caption
(b) gPg_{\mathrm{P}}
Figure 10: Comparison of analytic gains gIg_{\mathrm{I}} and gPg_{\mathrm{P}} with Monte-Carlo simulations at 1 GHz in the L-shaped corner in Fig. 3. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=90P_{\mathrm{th}}=-90 dBWm2\mathrm{dBWm^{-2}}. Markers are generated from Monte Carlo simulations. Solid lines are computed by Remark 2.
Refer to caption
(a) gIg_{\mathrm{I}}
Refer to caption
(b) gPg_{\mathrm{P}}
Figure 11: Comparison of analytic gains gIg_{\mathrm{I}} and gPg_{\mathrm{P}} with Monte-Carlo simulations at 28 GHz in the L-shaped corner in Fig. 3. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}. Markers are generated from Monte Carlo simulations. Solid lines are computed by Remark 2.

In the following sections, we use only the analytic results to propose some examples of building wireless performance analysis since the computation of gIg_{\mathrm{I}} and gPg_{\mathrm{P}} by Monte Carlo simulation is computationally expensive.

VI Simulations

In this section, we give examples of the evaluation of building wireless performance of two reference floor plans: a single rectangular room layout and an office building layout with rectangular rooms.

VI-A Performance evaluation for a single room

The impact of the area of a rectangular room on gIg_{\mathrm{I}} and gPg_{\mathrm{P}} is analyzed. It is needed to emphasize that gIg_{\mathrm{I}} and gPg_{\mathrm{P}} have different values at different positions in the room. To facilitate the comparison of rooms in terms of wireless performance, for each considered room, we compute the average gIg_{\mathrm{I}} and the average gPg_{\mathrm{P}} for that room.

Refer to caption
(a) gIg_{\mathrm{I}}
Refer to caption
(b) gPg_{\mathrm{P}}
Figure 12: Impact of room size on gIg_{\mathrm{I}} and gPg_{\mathrm{P}} at 1 GHz in a rectangular room. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=75P_{\mathrm{th}}=-75 dBWm2\mathrm{dBWm^{-2}}.
Refer to caption
(a) gIg_{\mathrm{I}}
Refer to caption
(b) gPg_{\mathrm{P}}
Figure 13: Impact of room size on gIg_{\mathrm{I}} and gPg_{\mathrm{P}} at 28 GHz in a rectangular room. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}.
Refer to caption
(a) gIg_{\mathrm{I}}
Refer to caption
(b) gPg_{\mathrm{P}}
Figure 14: Impact of room size on gIg_{\mathrm{I}} and gPg_{\mathrm{P}} at 1 GHz in a rectangular room. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}.

Figs. 14-14 show the average gIg_{\mathrm{I}} and the average gPg_{\mathrm{P}} as a function of the indoor coverage area for two different AR, i.e., AR{1,2}AR\in\{1,2\} at 1 GHz and 28 GHz. While covering a small room in the open space scenario, the power of the interference signals becomes larger in proportion to the intended signals due to the strong attenuation, and the limited receive sensitivity. Therefore, the significant increase of INI_{\mathrm{N}} takes place resulting in a decrease of gIg_{\mathrm{I}}. Whereas, in a larger room, the interference is more likely to be LOS, then a larger ILI_{\mathrm{L}} significantly reduces gIg_{\mathrm{I}}. Therefore, it is observed that gIg_{\mathrm{I}} can be optimized subject to the size of the room by balancing INI_{\mathrm{N}} and ILI_{\mathrm{L}} in an appropriate way. For different values of AR, the optimum size of a room and the maximized gIg_{\mathrm{I}} are different. With an increasing AR, the maximum gIg_{\mathrm{I}} decreases, and therefore a better wireless performance is achieved with a lower AR.

In Fig. 14, shown results correspond to RO=4.2R_{\mathrm{O}}=4.2 m, RL=5.3R_{\mathrm{L}}=5.3 m, and RN=2.5R_{\mathrm{N}}=2.5 m due to (12)-(14). In Fig. 14, we have RO=2.7R_{\mathrm{O}}=2.7 m, RL=3.2R_{\mathrm{L}}=3.2 m, and RN=1.9R_{\mathrm{N}}=1.9 m. Since the coverage area is in general less than the size of the room, the intended signal is less likely to be blocked. Thus, the value of the average gPg_{\mathrm{P}} is close to one. Therefore, with a small coverage area, e.g., a relative small PTPth\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}, the impact of gPg_{\mathrm{P}} is relatively slight. On the other hand, for a large PTPth\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}, e.g., PT=30dBWm2P_{\mathrm{T}}=-30\ \mathrm{dBWm^{-2}}, Pth=100dBWm2P_{\mathrm{th}}=-100\ \mathrm{dBWm^{-2}} at 1 GHz, the average gIg_{\mathrm{I}} and the average gPg_{\mathrm{P}} are shown in Fig. 14. Following (12)-(14), we have RO=75R_{\mathrm{O}}=75 m, RL=148R_{\mathrm{L}}=148 m, and RN=15R_{\mathrm{N}}=15 m. Even though for a room with an area of 50 m2 and an AR=2=2, the maximum length of LOS link is 525\sqrt{2} m, which is less than RNR_{\mathrm{N}}. Thus, all the interference power is NLOS, i.e., IL=0I_{\mathrm{L}}=0. Also, the INI_{\mathrm{N}} becomes independent of the size and the AR of the room because the NLOS interference is not impacted by any wall. Therefore, if the maximum possible distance between two points of walls around the polygon room is less than RNR_{\mathrm{N}}, INI_{\mathrm{N}} does not change with the size and the AR of the room and gIg_{\mathrm{I}} is a constant.

However, when the size of a room is less than RNR_{\mathrm{N}}, the blockage of the room attenuates the power of intended signals significantly. Reducing the size of the room, more intended signals are blocked and thus the gPg_{\mathrm{P}} decreases. Therefore, under this condition, a decreasing size reduces the wireless performance of a room.

Refer to caption
(a) gIg_{\mathrm{I}} at 1 GHz
Refer to caption
(b) gIg_{\mathrm{I}} at 28 GHz
Refer to caption
(c) gPg_{\mathrm{P}} at 1 GHz
Refer to caption
(d) gPg_{\mathrm{P}} at 28 GHz
Figure 15: gIg_{\mathrm{I}} and gPg_{\mathrm{P}} in a typical office scenario. PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}.

VI-B An example of building evaluation

In this subsection, the wireless performance evaluation for the A1 scenario of WINNER II channel model [43, Fig. 2.1] is taken as an example. In the computations, we use PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}}, and Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}. For each room or corridor, gIg_{\mathrm{I}} and gPg_{\mathrm{P}} are computed by evoking Remark 2. The obtained gIg_{\mathrm{I}} and gPg_{\mathrm{P}} are illustrated in Fig. 15. The average gI=12.03g_{\mathrm{I}}=12.03 and the average gP=0.74g_{\mathrm{P}}=0.74, respectively. Whereas for fc=28f_{\mathrm{c}}=28 GHz, i.e., at mm-Wave band, the average gI=7.08g_{\mathrm{I}}=7.08 and the average gP=1.01g_{\mathrm{P}}=1.01, respectively. For PT=30P_{\mathrm{T}}=-30 dBWm2\mathrm{dBWm^{-2}} and Pth=100P_{\mathrm{th}}=-100 dBWm2\mathrm{dBWm^{-2}}, the WINNER II A1 at 1 GHz outperforms that at 28 GHz in terms of gIg_{\mathrm{I}} because it blocks almost all interference signals. Whereas in terms of gPg_{\mathrm{P}}, the WINNER II A1 building has better wireless performance at 28 GHz. Therefore, for a building wireless performance evaluation in the future, the tradeoff between gIg_{\mathrm{I}}, gPg_{\mathrm{P}}, and fcf_{\mathrm{c}} has to be taken into account carefully.

VII Conclusions

The intrinsic wireless performance characteristics of buildings have been formally defined as the interference gain and the power gain figures of merit. An analytical approach has been proposed to compute these figures of merit reflecting the wireless performance of buildings with polygonal-shaped rooms. The analytic approach has been validated through Monte Carlo simulations for typical indoor environments with well-established empirical path gain models showing a good agreement between them. The proposed analytic approach allows to quantify the impact of building properties, e.g., the aspect ratio and the area of rooms, on both the interference gain and the power gain. Numerical examples are provided based on different propagation scenarios at different frequencies, i.e., 1 GHz and 28 GHz. Based on the presented analysis it is shown that wireless building design need to consider the trade-off between the interference gain, the power gain and the frequency band of the wireless system deployed in a specific building. This finding has tremendous implications on the way new buildings shall be designed taking into consideration wireless performance already at the design phase. Indeed, it is expected that architects and civil engineers will now be able to carefully take the interference gain and the power gain into account before designing or retrofitting a building to provide adequate performance margins for future indoor wireless networks. The proposed approach bridges the building design and wireless communications industries. It has the potential to open up a wide range of innovative solutions offering wireless friendly buildings, leading to improved wireless services indoors where most of the communications traffic takes place.

Appendices

VII-A Proof of Theorem 1

VII-A1 Closed-form IOI_{\mathrm{O}} (26)

IOI_{\mathrm{O}} is defined by [15, 16]

IO=02πRO+PTGO(R)RdRdθ.\displaystyle I_{\mathrm{O}}=\int_{0}^{2\pi}\int_{R_{\mathrm{O}}}^{+\infty}P_{\mathrm{T}}G_{\mathrm{O}}(R)R\mathrm{d}R\mathrm{d}\theta. (68)

Substituting (12) into (68), considering the case PTPth<(4π)4(hThR)2λ4\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}<\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}}, we have

IO=02πλ4πPTPth+PTGO(R)RdRdθ.\displaystyle I_{\mathrm{O}}=\int_{0}^{2\pi}\int_{\frac{\lambda}{4\pi}\sqrt{\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}}}^{+\infty}{P_{\mathrm{T}}}{G_{\mathrm{O}}(R)}R\mathrm{d}R\mathrm{d}\theta. (69)

Substituting (2) into (69), we obtain the closed-form IOI_{\mathrm{O}} as

IO=PTλ28π{12+ln(16Pthπ2hThRPTλ2)}.\displaystyle I_{\mathrm{O}}=\frac{P_{\mathrm{T}}\lambda^{2}}{8\pi}\left\{\frac{1}{2}+\ln\left(\frac{16\sqrt{P_{\mathrm{th}}}\pi^{2}h_{\mathrm{T}}h_{\mathrm{R}}}{\sqrt{P_{\mathrm{T}}}\lambda^{2}}\right)\right\}. (70)

For PTPth(4π)4(hThR)2λ4\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\geq\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}}, we have

IO=02π(PTPth)14hThR+PTGO(R)RdRdθ.\displaystyle I_{\mathrm{O}}=\int_{0}^{2\pi}\int_{\left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)^{\frac{1}{4}}\sqrt{h_{\mathrm{T}}h_{\mathrm{R}}}}^{+\infty}{P_{\mathrm{T}}}{G_{\mathrm{O}}(R)}R\mathrm{d}R\mathrm{d}\theta. (71)

Substituting (2) into (71) and following some algebraic manipulations, we obtain

IO=πPthPThThR.\displaystyle I_{\mathrm{O}}=\pi\sqrt{P_{\mathrm{th}}P_{\mathrm{T}}}h_{\mathrm{T}}h_{\mathrm{R}}. (72)

Combining (70), and (72), we obtain (26).

VII-A2 Closed-form POP_{\mathrm{O}} (29)

POP_{\mathrm{O}} is defined by [15]

PO=02π0ROPTGO(R)RdRdθ.\displaystyle P_{\mathrm{O}}=\int_{0}^{2\pi}\int_{0}^{R_{\mathrm{O}}}P_{\mathrm{T}}G_{\mathrm{O}}(R)R\mathrm{d}R\mathrm{d}\theta. (73)

Substituting (12) into (73), for PTPth<(4π)4(hThR)2λ4\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}<\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}}, we have

PO=02π0λ4πPTPthPTGO(R)RdRdθ.\displaystyle P_{\mathrm{O}}=\int_{0}^{2\pi}\int_{0}^{\frac{\lambda}{4\pi}\sqrt{\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}}}{P_{\mathrm{T}}}{G_{\mathrm{O}}(R)}R\mathrm{d}R\mathrm{d}\theta. (74)

Substituting (2) into (74), POP_{\mathrm{O}} is then computed by

PO=2π0λ4πPTRdR+2πλ4πλ4πPTPthPT(λ4π)2R1dR.\displaystyle P_{\mathrm{O}}=2\pi\int_{0}^{\frac{\lambda}{4\pi}}P_{\mathrm{T}}R\mathrm{d}R+2\pi\int_{\frac{\lambda}{4\pi}}^{\frac{\lambda}{4\pi}\sqrt{\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}}}P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}R^{-1}\mathrm{d}R. (75)

Following some algebraic manipulations, we obtain (76).

PO=PTλ216π[1+ln(PTPth)].\displaystyle P_{\mathrm{O}}=\frac{P_{\mathrm{T}}\lambda^{2}}{16\pi}\left[1+\ln\left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)\right]. (76)

For PTPth(4π)4(hThR)2λ4\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\geq\frac{(4\pi)^{4}\left({h_{\mathrm{T}}h_{\mathrm{R}}}\right)^{2}}{\lambda^{4}}, we have

PO=02π0(PTPth)14hThRPTGO(R)RdRdθ.\displaystyle P_{\mathrm{O}}=\int_{0}^{2\pi}\int_{0}^{\left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)^{\frac{1}{4}}\sqrt{h_{\mathrm{T}}h_{\mathrm{R}}}}{P_{\mathrm{T}}}{G_{\mathrm{O}}(R)}R\mathrm{d}R\mathrm{d}\theta. (77)

Substituting (2) into (77), we obtain

PO=2π0λ4πPTRdR+2πλ4π4πhThRλPT(λ4π)2R1dR\displaystyle P_{\mathrm{O}}=2\pi\int_{0}^{\frac{\lambda}{4\pi}}P_{\mathrm{T}}R\mathrm{d}R+2\pi\int_{\frac{\lambda}{4\pi}}^{\frac{4\pi h_{\mathrm{T}}h_{\mathrm{R}}}{\lambda}}P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}R^{-1}\mathrm{d}R (78)
+2π4πhThRλ(PTPth)14hThRPT(hThR)2R3dR.\displaystyle\hskip 25.29494pt+2\pi\int_{\frac{4\pi h_{\mathrm{T}}h_{\mathrm{R}}}{\lambda}}^{\left(\frac{P_{\mathrm{T}}}{P_{\mathrm{th}}}\right)^{\frac{1}{4}}\sqrt{h_{\mathrm{T}}h_{\mathrm{R}}}}{P_{\mathrm{T}}}{(h_{\mathrm{T}}h_{\mathrm{R}})^{2}R^{-3}}\mathrm{d}R. (79)

Following some algebraic manipulations, we obtain

PO=PTλ216π+PTλ28πln(16π2hThRλ2)πPth2+λ4PT29π3(hThR)2.\displaystyle P_{\mathrm{O}}=\frac{P_{\mathrm{T}}\lambda^{2}}{16\pi}+\frac{P_{\mathrm{T}}\lambda^{2}}{8\pi}\ln\left(\frac{16\pi^{2}h_{\mathrm{T}}h_{\mathrm{R}}}{\lambda^{2}}\right)-\frac{\pi P_{\mathrm{th}}}{2}+\frac{\lambda^{4}{P_{\mathrm{T}}}}{2^{9}\pi^{3}(h_{\mathrm{T}}h_{\mathrm{R}})^{2}}. (80)

Combining (76), and (80), we obtain (29).

Refer to caption
Figure 16: Computation of PLP_{\mathrm{L}} for 1D0<RL1\leq D_{0}<R_{\mathrm{L}}.
Refer to caption
Figure 17: Computation of PLP_{\mathrm{L}} for λ4πD0<1\frac{\lambda}{4\pi}\leq D_{0}<1.

VII-B Proof of Theorem 2

Here, we derive PL,nTMP_{\mathrm{L},n_{\mathrm{TM}}} in a closed form. Substituting (38) into (34), we have

PL,nTM=θlθr0min{RL,D0cosθ}PTGL(R)RdRdθ.\displaystyle P_{\mathrm{L},n_{\mathrm{TM}}}=\int_{\theta_{\mathrm{l}}}^{\theta_{\mathrm{r}}}\int_{0}^{\min\left\{R_{\mathrm{L}},\frac{D_{0}}{\cos\theta}\right\}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}R\mathrm{d}R\mathrm{d}\theta. (81)

VII-B1 D0RLD_{0}\geq R_{\mathrm{L}}

For D0>RLD_{0}>R_{\mathrm{L}}, we have min{RL,D0cosθ}=RL\min\left\{R_{\mathrm{L}},\frac{D_{0}}{\cos\theta}\right\}=R_{\mathrm{L}}, and therefore

PL,nTM=θlθr0RLPTGL(R)RdRdθ.\displaystyle P_{\mathrm{L},n_{\mathrm{TM}}}=\int_{\theta_{\mathrm{l}}}^{\theta_{\mathrm{r}}}\int_{0}^{R_{\mathrm{L}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}R\mathrm{d}R\mathrm{d}\theta. (82)

Following some algebraic manipulations, using (2), and the following Lemma 1, we obtain

PL,nTM=PL,1,\displaystyle P_{\mathrm{L},n_{\mathrm{TM}}}=P_{\mathrm{L,1}}, (83)

where PL,1P_{\mathrm{L,1}} is given in (42).

Lemma 1.

The function Z0Z_{0} is defined as

Z0(z1,z2,z3,z4,z5)z1z2z3z4Rz5dRdθ,\displaystyle Z_{0}(z_{1},z_{2},z_{3},z_{4},z_{5})\triangleq\int_{z_{1}}^{z_{2}}\int_{z_{3}}^{z_{4}}R^{-z_{5}}\mathrm{d}R\mathrm{d}\theta, (84)

where π2<z1<z2<π2-\frac{\pi}{2}<z_{1}<z_{2}<\frac{\pi}{2}, z3>0z_{3}>0, z4>0z_{4}>0, and z5>1z_{5}>-1. The closed-form expression is given by

Z0(z1,z2,z3,z4,z5){(z2z1)ln(z4z3),z5=1,(z2z1)(z41z5z31z5)1z5,z51.\displaystyle Z_{0}(z_{1},z_{2},z_{3},z_{4},z_{5})\equiv\left\{\begin{array}[]{ll}(z_{2}-z_{1})\ln\left(\frac{z_{4}}{z_{3}}\right),&z_{5}=1,\\ \frac{(z_{2}-z_{1})(z_{4}^{1-z_{5}}-z_{3}^{1-z_{5}})}{1-z_{5}},&z_{5}\neq 1.\end{array}\right. (87)

The proof of this Lemma is straightforward and therefore omitted here.

VII-B2 1D0<RL1\leq D_{0}<R_{\mathrm{L}}

For 1D0<RL1\leq D_{0}<R_{\mathrm{L}}, an example of the toy model is shown in Fig. 17. Under this condition,

PL,nTM=ΩPL1PTGL(R)dΩΩPL2PT(λ4π)2RnLdΩ,\displaystyle P_{\mathrm{L},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{PL1}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{PL2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{L}}}}\mathrm{d}\Omega, (88)

where ΩPL1\Omega_{\mathrm{PL1}} and ΩPL2\Omega_{\mathrm{PL2}} are given by

ΩPL1{(R,θ)|(R<RL)(π2<θl<θ<θr<π2)},\displaystyle\Omega_{\mathrm{PL1}}\triangleq\left\{\left(R,\theta\right)\left|\left(R<R_{\mathrm{L}}\right)\land\left(-\frac{\pi}{2}<\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}<\frac{\pi}{2}\right)\right.\right\}, (89)
ΩPL2{(R,θ)|(D0cosθR<RL)(max{θL1,θl}<θ<min{θL1,θr})},\displaystyle\Omega_{\mathrm{PL2}}\triangleq\left\{\left(R,\theta\right)\left|\left(\frac{D_{0}}{\cos\theta}\leq R<R_{\mathrm{L}}\right)\land\left(\max\{-\theta_{\mathrm{L1}},\theta_{\mathrm{l}}\}<\theta<\min\{\theta_{\mathrm{L1}},\theta_{\mathrm{r}}\}\right)\right.\right\}, (90)

where θL1\theta_{\mathrm{L1}} is defined in (45). Next, following some algebraic manipulations, we obtain

ΩPL1PTGL(R)dΩ=PL,1,\displaystyle\int_{\Omega_{\mathrm{PL1}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}\mathrm{d}\Omega=P_{\mathrm{L,1}}, (91)

and

ΩPL2PTGL(R)dΩ={max{θL1,θl}min{θL1,θr}D0cosθRLPT(λ4π)2RnL1dRdθ,(θl<θL1)(θr>θL1),0,else,\displaystyle\int_{\Omega_{\mathrm{PL2}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}\mathrm{d}\Omega=\left\{\begin{array}[]{ll}\int_{\max\{-\theta_{\mathrm{L1}},\theta_{\mathrm{l}}\}}^{\min\{\theta_{\mathrm{L1}},\theta_{\mathrm{r}}\}}\int_{\frac{D_{0}}{\cos\theta}}^{R_{\mathrm{L}}}\frac{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{n_{\mathrm{L}}-1}}\mathrm{d}R\mathrm{d}\theta,&(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L1}}),\\ 0,&\textrm{else},\end{array}\right. (94)
={PL,2,(θl<θL1)(θr>θL1),0,else,\displaystyle\hskip 93.95122pt=\left\{\begin{array}[]{ll}-P_{\mathrm{L,2}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L1}}),\\ 0,&\textrm{else},\end{array}\right. (97)

where PL,2P_{\mathrm{L,2}} is given in (43), and Z1Z_{1} is computed by the following Lemma 2.

Lemma 2.

The function Z1Z_{1} is defined as follows

Z1(z1,z2,z3,z4,z5)z1z2z3z4cos(θ)Rz5dRdθ,\displaystyle Z_{1}(z_{1},z_{2},z_{3},z_{4},z_{5})\triangleq\int_{z_{1}}^{z_{2}}\int_{z_{3}}^{\frac{z_{4}}{\cos(\theta)}}R^{-z_{5}}\mathrm{d}R\mathrm{d}\theta, (98)

where π2<z1<z2<π2-\frac{\pi}{2}<z_{1}<z_{2}<\frac{\pi}{2}, z3>0z_{3}>0, z4>0z_{4}>0, and z51z_{5}\geq-1. The closed-form expression is given by

Z1(z1,z2,z3,z4,z5){z42[tan(z2)tan(z1)]2+(z1z2)z322,z5=1,z4ln(tan(z2)+sec(z2)tan(z1)+sec(z1))+(z1z2)z3,z5=0,(z2z1)ln(2z4z3)+Im[Li2(e2jz2)Li2(e2jz1)]2,z5=1,{B(12,z52)[sgn(z2)sgn(z1)]z41z52(1z5)+(z1z2)z31z51z5+z41z5(1z5)z5×[sgn(z1)cosz5(z1)2F1(z52,12;z5+22,cos2(z1))sgn(z2)cosz5(z2)2F1(z52,12;z5+22,cos2(z2))],else.\displaystyle Z_{1}(z_{1},z_{2},z_{3},z_{4},z_{5})\equiv\left\{\begin{array}[]{ll}\frac{z_{4}^{2}[\tan(z_{2})-\tan(z_{1})]}{2}+\frac{(z_{1}-z_{2})z_{3}^{2}}{2},&z_{5}=-1,\\ z_{4}\ln\left(\frac{\tan(z_{2})+\sec(z_{2})}{\tan(z_{1})+\sec(z_{1})}\right)+(z_{1}-z_{2})z_{3},&z_{5}=0,\\ (z_{2}-z_{1})\ln\left(\frac{2z_{4}}{z_{3}}\right)+\frac{\mathrm{Im}\left[\mathrm{Li}_{2}(-e^{2jz_{2}})-\mathrm{Li}_{2}(-e^{2jz_{1}})\right]}{2},&z_{5}=1,\\ \left\{\begin{array}[]{l}\frac{B\left(\frac{1}{2},\frac{z_{5}}{2}\right)\left[\mathrm{sgn}(z_{2})-\mathrm{sgn}(z_{1})\right]z_{4}^{1-z_{5}}}{2(1-z_{5})}\\ +\frac{(z_{1}-z_{2})z_{3}^{1-z_{5}}}{1-z_{5}}+\frac{z_{4}^{1-z_{5}}}{(1-z_{5})z_{5}}\times\\ \left[\mathrm{sgn}(z_{1})\cos^{z_{5}}(z_{1})\ _{2}F_{1}\left(\frac{z_{5}}{2},\frac{1}{2};\frac{z_{5}+2}{2},\cos^{2}(z_{1})\right)\right.\\ \left.-\mathrm{sgn}(z_{2})\cos^{z_{5}}(z_{2})\ _{2}F_{1}\left(\frac{z_{5}}{2},\frac{1}{2};\frac{z_{5}+2}{2},\cos^{2}(z_{2})\right)\right],\end{array}\right.&\mathrm{else}.\end{array}\right. (107)

The Lemma is straightforward to verify and the proof of this Lemma is omitted here.

Substituting (91) and (94) into (88), we obtain

PL,nTM={PL,1+PL,2,(θl<θL1)(θr>θL1),PL,1,else.\displaystyle P_{\mathrm{L},n_{\mathrm{TM}}}=\left\{\begin{array}[]{ll}P_{\mathrm{L,1}}+P_{\mathrm{L,2}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L1}}),\\ P_{\mathrm{L,1}},&\textrm{else}.\end{array}\right. (110)

VII-B3 λ4πD0<1\frac{\lambda}{4\pi}\leq D_{0}<1

An example of the toy model is shown in Fig. 17, for λ4πD0<1\frac{\lambda}{4\pi}\leq D_{0}<1. Under this condition, PL,nTMP_{\mathrm{L},n_{\mathrm{TM}}} is given by

PL,nTM=ΩPL1PTGL(R)dΩΩPL2PT(λ4π)2RnLdΩ+ΩPL3PT(λ4π)2RnLdΩΩPL3PT(λ4π)2R2dΩ,\displaystyle\begin{array}[]{l}P_{\mathrm{L},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{PL1}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{PL2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{L}}}}\mathrm{d}\Omega\\ \hskip 36.135pt+\int_{\Omega_{\mathrm{PL3}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{L}}}}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{PL3}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}R^{-2}\mathrm{d}\Omega,\end{array} (113)

where the first term ΩPL1PTGL(R)dΩ\int_{\Omega_{\mathrm{PL1}}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}\mathrm{d}\Omega and the second term ΩPL2PT(λ4π)2RnLdΩ\int_{\Omega_{\mathrm{PL2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{L}}}}\mathrm{d}\Omega are computed following steps of computing (91) and (94), respectively. For the third and the forth terms, ΩPL3\Omega_{\mathrm{PL3}} is defined as

ΩPL3{(R,θ)|(D0cosθR<1)(max{θL2,θl}<θ<min{θL2,θr})},\displaystyle\Omega_{\mathrm{PL3}}\triangleq\left\{\left(R,\theta\right)\left|\left(\frac{D_{0}}{\cos\theta}\leq R<1\right)\land\left(\max\{-\theta_{\mathrm{L2}},\theta_{\mathrm{l}}\}<\theta<\min\{\theta_{\mathrm{L2}},\theta_{\mathrm{r}}\}\right)\right.\right\}, (114)

where θL2\theta_{\mathrm{L2}} is defined in (46). Further substituting (114) into

ΩPL3PT(λ4π)2RnLdΩΩPL3PT(λ4π)2R2dΩ,\displaystyle\int_{\Omega_{\mathrm{PL3}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{L}}}}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{PL3}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}R^{-2}}\mathrm{d}\Omega, (115)

and following some straightforward mathematical manipulations, we have

ΩPL3PT(λ4π)2RnLdΩΩPL3PT(λ4π)2R2dΩ={PL,3,(θl<θL2)(θr>θL2),0,else,\displaystyle\int_{\Omega_{\mathrm{PL3}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{L}}}}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{PL3}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}R^{-2}\mathrm{d}\Omega=\left\{\begin{array}[]{ll}P_{\mathrm{L,3}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{L2}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L2}}),\\ 0,&\textrm{else},\end{array}\right. (118)

where PL,3P_{\mathrm{L,3}} is given in (2). Therefore, substituting (91), (94) and (118) into (113), we have

PL,nTM={PL,1+PL,2+PL,3,(θl<θL2)(θr>θL2),PL,1+PL,2,(θl<θL1)(θr>θL1)[(θl>θL2)(θr<θL2)],PL,1,else.\displaystyle P_{\mathrm{L},n_{\mathrm{TM}}}=\left\{\begin{array}[]{ll}P_{\mathrm{L,1}}+P_{\mathrm{L,2}}+P_{\mathrm{L,3}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{L2}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L2}}),\\ P_{\mathrm{L,1}}+P_{\mathrm{L,2}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L1}})\land\left[(\theta_{\mathrm{l}}>\theta_{\mathrm{L2}})\lor(\theta_{\mathrm{r}}<-\theta_{\mathrm{L2}})\right],\\ P_{\mathrm{L,1}},&\textrm{else}.\end{array}\right. (122)

Combining (83), (110) and (122), we obtain TABLE I and therefore Theorem 2 is proved.

VII-C Proof of Theorem 3

Here, we derive IL,nTMI_{\mathrm{L},n_{\mathrm{TM}}} in closed form. Since RL>1R_{\mathrm{L}}>1 is assumed, for R>RLR>R_{\mathrm{L}}, we have GL(R)=(λ4π)2RnLG_{\mathrm{L}}(R)=\left(\frac{\lambda}{4\pi}\right)^{2}R^{-n_{\mathrm{L}}}. Substituting (39) into (35), we have

IL,nTM=(θl<θ<θr)(RL<D0cosθ)RLD0cosθPT(λ4π)2R1nLdRdθ.\displaystyle I_{\mathrm{L},n_{\mathrm{TM}}}=\int_{(\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}})\land(R_{\mathrm{L}}<\frac{D_{0}}{\cos\theta})}\int_{R_{\mathrm{L}}}^{\frac{D_{0}}{\cos\theta}}{P_{\mathrm{T}}}\left(\frac{\lambda}{4\pi}\right)^{2}R^{1-n_{\mathrm{L}}}\mathrm{d}R\mathrm{d}\theta. (123)

VII-C1 RLD0R_{\mathrm{L}}\leq{D_{0}}

For RL<D0R_{\mathrm{L}}<{D_{0}}, (RL<D0cosθ)=(R_{\mathrm{L}}<\frac{D_{0}}{\cos\theta})=\top, and

IL,nTM=θlθrRLD0cosθPTGL(R)RdRdθ=IL,1,\displaystyle I_{\mathrm{L},n_{\mathrm{TM}}}=\int_{\theta_{\mathrm{l}}}^{\theta_{\mathrm{r}}}\int_{R_{\mathrm{L}}}^{\frac{D_{0}}{\cos\theta}}{P_{\mathrm{T}}}{G_{\mathrm{L}}(R)}R\mathrm{d}R\mathrm{d}\theta=I_{\mathrm{L},1}, (124)

where IL,1I_{\mathrm{L},1} is computed by (59).

VII-C2 RL>D0R_{\mathrm{L}}>{D_{0}}

For RL>D0R_{\mathrm{L}}>{D_{0}}, we can show (125) by using straightforward mathematical manipulations.

(θl<θ<θr)(RL<D0cosθ)=(θl<θ<θr)(|θ|>|θL1|)\displaystyle\left(\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}\right)\land\left(R_{\mathrm{L}}<\frac{D_{0}}{\cos\theta}\right)=\left(\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}\right)\land\left(|\theta|>|\theta_{\mathrm{L1}}|\right) (125)
=[(θl<θL1<θr)(θl<θ<θL1)]\displaystyle\hskip 151.76744pt=\left[\left(\theta_{\mathrm{l}}<-\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right)\land\left(\theta_{\mathrm{l}}<\theta<-\theta_{\mathrm{L1}}\right)\right]
[(θl<θL1<θr)(θL1<θ<θr)]\displaystyle\hskip 151.76744pt\lor\left[\left(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right)\land\left(\theta_{\mathrm{L1}}<\theta<\theta_{\mathrm{r}}\right)\right]
{[(θr<θL1)(θL1<θl)](θl<θ<θr)}.\displaystyle\hskip 151.76744pt\lor\left\{\left[\left(\theta_{\mathrm{r}}<-\theta_{\mathrm{L1}}\right)\lor\left(\theta_{\mathrm{L1}}<\theta_{\mathrm{l}}\right)\right]\land\left(\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}\right)\right\}.

Substituting (125) into (123), we obtain the closed-form expression of IL,nTMI_{\mathrm{L},n_{\mathrm{TM}}} for RL>D0R_{\mathrm{L}}>{D_{0}} as

IL,nTM={IL,1,(θr<θL1)(θL1<θl),IL,2,(θl<θL1<θr)(θl<θL1<θr),IL,3,(θl<θL1<θr)(θl<θL1<θr),IL,2+IL,3,(θl<θL1<θr)(θl<θL1<θr).\displaystyle I_{\mathrm{L},n_{\mathrm{TM}}}=\left\{\begin{array}[]{l l}I_{\mathrm{L},1},&\left(\theta_{\mathrm{r}}<-\theta_{\mathrm{L1}}\right)\lor\left(\theta_{\mathrm{L1}}<\theta_{\mathrm{l}}\right),\\ I_{\mathrm{L},2},&\left(\theta_{\mathrm{l}}<-\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right)\lor\left(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right),\\ I_{\mathrm{L},3},&\left(\theta_{\mathrm{l}}<-\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right)\lor\left(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right),\\ I_{\mathrm{L},2}+I_{\mathrm{L},3},&\left(\theta_{\mathrm{l}}<-\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right)\lor\left(\theta_{\mathrm{l}}<\theta_{\mathrm{L1}}<\theta_{\mathrm{r}}\right).\end{array}\right. (130)

where IL,2I_{\mathrm{L},2} and IL,3I_{\mathrm{L},3} are respectively computed by (60) and (61). Combining (124) and (130), we obtain TABLE II, and therefore Theorem 3 is proved.

VII-D Proof of Theorem 4

We now derive PN,nTMP_{\mathrm{N},n_{\mathrm{TM}}} in a closed form. Substituting (40) into (36), we have

PN,nTM=(θl<θ<θr)(RN>D0cosθ)D0cosθRNPTGN(R)RdRdθ.\displaystyle P_{\mathrm{N},n_{\mathrm{TM}}}=\int_{(\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}})\land(R_{\mathrm{N}}>\frac{D_{0}}{\cos\theta})}\int_{\frac{D_{0}}{\cos\theta}}^{R_{\mathrm{N}}}{P_{\mathrm{T}}}{G_{\mathrm{N}}(R)}R\mathrm{d}R\mathrm{d}\theta. (131)

VII-D1 D0RND_{0}\geq R_{\mathrm{N}}

For D0>RND_{0}>R_{\mathrm{N}}, (RN>D0cosθ)=(R_{\mathrm{N}}>\frac{D_{0}}{\cos\theta})=\bot, and therefore

PN,nTM=0.\displaystyle P_{\mathrm{N},n_{\mathrm{TM}}}=0. (132)

VII-D2 1D0<RN1\leq D_{0}<R_{\mathrm{N}}

For 1D0<RN1\leq D_{0}<R_{\mathrm{N}}

PN,nTM=ΩPN1PT(λ4π)2RnNdΩ,\displaystyle P_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{PN1}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega, (133)

where θL1\theta_{\mathrm{L1}} is defined in (45), and ΩPN1\Omega_{\mathrm{PN1}} is defined as

ΩPN1{(R,θ)|(D0cosθR<RN)(max{θN1,θl}<θ<min{θN1,θr})}.\displaystyle\Omega_{\mathrm{PN1}}\triangleq\left\{\left(R,\theta\right)\left|\left(\frac{D_{0}}{\cos\theta}\leq R<R_{\mathrm{N}}\right)\land\left(\max\{-\theta_{\mathrm{N1}},\theta_{\mathrm{l}}\}<\theta<\min\{\theta_{\mathrm{N1}},\theta_{\mathrm{r}}\}\right)\right.\right\}. (134)

Following some algebraic manipulations, we obtain

PN,nTM=ΩPN1PT(λ4π)2RnNdΩ={PN,1,(θl<θN1)(θr>θN1),0,else,\displaystyle P_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{PN1}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega=\left\{\begin{array}[]{ll}P_{\mathrm{N,1}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{N1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{N1}}),\\ 0,&\textrm{else},\end{array}\right. (137)

where Z1Z_{1} is defined in (98) and computed by (107). PN,1P_{\mathrm{N,1}} is computed by (62).

VII-D3 λ4πD0<1\frac{\lambda}{4\pi}\leq D_{0}<1

For λ4πD0<1\frac{\lambda}{4\pi}\leq D_{0}<1, an example of the toy model is shown in Fig. 19. Under this condition,

PN,nTM=ΩPN1PT(λ4π)2RnNdΩΩPN2PT(λ4π)2RnNdΩ+ΩPN2PT(λ4π)2R2dΩ,\displaystyle P_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{PN1}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{PN2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega+\int_{\Omega_{\mathrm{PN2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}R^{-2}\mathrm{d}\Omega, (138)

where θN2\theta_{\mathrm{N2}} is defined in (65), ΩPN2PT(λ4π)2RnNdΩ\int_{\Omega_{\mathrm{PN2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega is computed by (137), and ΩPN2\Omega_{\mathrm{PN2}} is defined as

ΩPN2{(R,θ)|(D0cosθR<1)(max{θN2,θl}<θ<min{θN2,θr})}.\displaystyle\Omega_{\mathrm{PN2}}\triangleq\left\{\left(R,\theta\right)\left|\left(\frac{D_{0}}{\cos\theta}\leq R<1\right)\land\left(\max\{-\theta_{\mathrm{N2}},\theta_{\mathrm{l}}\}<\theta<\min\{\theta_{\mathrm{N2}},\theta_{\mathrm{r}}\}\right)\right.\right\}. (139)
Refer to caption
Figure 18: Computation of PNP_{\mathrm{N}} for λ4πD0<1\frac{\lambda}{4\pi}\leq D_{0}<1.
Refer to caption
Figure 19: Computation of INI_{\mathrm{N}} for 1D0<RN1\leq D_{0}<R_{\mathrm{N}}.

Following some straightforward computations by substituting (139) into the second and the third terms of (138), we have

ΩPN2PT(λ4π)2R2dΩΩPN2PT(λ4π)2RnNdΩ={PN,2,(θl<θL2)(θr>θL2),0,else,\displaystyle\int_{\Omega_{\mathrm{PN2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-2}}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{PN2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega=\left\{\begin{array}[]{ll}P_{\mathrm{N,2}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{L2}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L2}}),\\ 0,&\textrm{else},\end{array}\right. (142)

where PN,2P_{\mathrm{N,2}} is computed by (4). Therefore, for λ4πD0<1\frac{\lambda}{4\pi}\leq D_{0}<1, PN,nTMP_{\mathrm{N},n_{\mathrm{TM}}} is given by

PN,nTM={PN,1+PN,2,(θl<θN2)(θr>θN2),PN,1,(θl<θN1)(θr>θN1)[(θl>θN2)(θr<θN2)],0,else.\displaystyle P_{\mathrm{N},n_{\mathrm{TM}}}=\left\{\begin{array}[]{ll}P_{\mathrm{N,1}}+P_{\mathrm{N,2}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{N2}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{N2}}),\\ P_{\mathrm{N,1}},&(\theta_{\mathrm{l}}<\theta_{\mathrm{N1}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{N1}})\land\left[(\theta_{\mathrm{l}}>\theta_{\mathrm{N2}})\lor(\theta_{\mathrm{r}}<-\theta_{\mathrm{N2}})\right],\\ 0,&\textrm{else}.\end{array}\right. (146)

Combining (132), (137) and (146), we obtain TABLE III and therefore Theorem 4 is proved.

VII-E Proof of Theorem 5

Here, we derive IN,nTMI_{\mathrm{N},n_{\mathrm{TM}}} in a closed form. Substituting (41) into (37), we have

IN,nTM=θlθrmax{RL,D0cosθ}+PTGN(R)RdRdθ.\displaystyle I_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\theta_{\mathrm{l}}}^{\theta_{\mathrm{r}}}\int_{\max\left\{R_{\mathrm{L}},\frac{D_{0}}{\cos\theta}\right\}}^{+\infty}{P_{\mathrm{T}}}{G_{\mathrm{N}}(R)}R\mathrm{d}R\mathrm{d}\theta. (147)

VII-E1 RN<D0R_{\mathrm{N}}<{D_{0}}

For RN<D0R_{\mathrm{N}}<{D_{0}}, we have (max{RL,D0cosθ}=D0cosθ)=\left(\max\left\{R_{\mathrm{L}},\frac{D_{0}}{\cos\theta}\right\}=\frac{D_{0}}{\cos\theta}\right)=\top. Therefore, we obtain

IN,nTM=θlθrD0cosθ+PT(λ4π)2R1nNdRdθ=IN,1,\displaystyle I_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\theta_{\mathrm{l}}}^{\theta_{\mathrm{r}}}\int_{\frac{D_{0}}{\cos\theta}}^{+\infty}{P_{\mathrm{T}}}\left(\frac{\lambda}{4\pi}\right)^{2}R^{1-n_{\mathrm{N}}}\mathrm{d}R\mathrm{d}\theta=I_{\mathrm{N},1}, (148)

where IN,1I_{\mathrm{N},1} is computed by (66).

VII-E2 RN>D0R_{\mathrm{N}}>{D_{0}}

For RN>D0R_{\mathrm{N}}>{D_{0}}, an example of the toy model is shown in Fig. 19. Under this condition,

IN,nTM=ΩIN1PT(λ4π)2RnNdΩΩIN2PT(λ4π)2RnNdΩ,\displaystyle I_{\mathrm{N},n_{\mathrm{TM}}}=\int_{\Omega_{\mathrm{IN1}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega-\int_{\Omega_{\mathrm{IN2}}}{P_{\mathrm{T}}\left(\frac{\lambda}{4\pi}\right)^{2}}{R^{-n_{\mathrm{N}}}}\mathrm{d}\Omega, (149)

where ΩIN1\Omega_{\mathrm{IN1}} and ΩIN2\Omega_{\mathrm{IN2}} are defined as

ΩIN1{(R,θ)|(RD0cosθ)(π2<θl<θ<θr<π2)},\displaystyle\Omega_{\mathrm{IN1}}\triangleq\left\{\left(R,\theta\right)\left|\left(R\geq\frac{D_{0}}{\cos\theta}\right)\land\left(-\frac{\pi}{2}<\theta_{\mathrm{l}}<\theta<\theta_{\mathrm{r}}<\frac{\pi}{2}\right)\right.\right\}, (150)
ΩIN2{(R,θ)|(D0cosθR<RN)(max{θN1,θl}<θ<min{θN1,θr})}.\displaystyle\Omega_{\mathrm{IN2}}\triangleq\left\{\left(R,\theta\right)\left|\left(\frac{D_{0}}{\cos\theta}\leq R<R_{\mathrm{N}}\right)\land\left(\max\left\{-\theta_{\mathrm{N1}},\theta_{\mathrm{l}}\}<\theta<\min\{\theta_{\mathrm{N1}},\theta_{\mathrm{r}}\right\}\right)\right.\right\}. (151)

Following some algebraic manipulations, we have

IN,nTM={IN,1+IN,2,(θl<θL2)(θr>θL2),IN,1,(θl>θL2)(θr<θL2),\displaystyle I_{\mathrm{N},n_{\mathrm{TM}}}=\left\{\begin{array}[]{ll}I_{\mathrm{N},1}+I_{\mathrm{N},2},&(\theta_{\mathrm{l}}<\theta_{\mathrm{L2}})\land(\theta_{\mathrm{r}}>-\theta_{\mathrm{L2}}),\\ I_{\mathrm{N},1},&(\theta_{\mathrm{l}}>\theta_{\mathrm{L2}})\lor(\theta_{\mathrm{r}}<-\theta_{\mathrm{L2}}),\end{array}\right. (154)

where IN,2I_{\mathrm{N},2} is computed by (67). Combining (148) and (154), we obtain TABLE IV, and therefore Theorem 5 is proved.

References

  • [1] Cisco, “Cisco vision: 5G-thriving indoors,” [Online], https://www.cisco.com/c/dam/en/us/solutions/collateral/service-provider/ultra-services-platform/5g-ran-indoor.pdf.
  • [2] M. Giordani, et al., “Toward 6G networks: Use cases and technologies” IEEE Commun. Mag., vol. 58, no. 3, pp. 55-61, 2020.
  • [3] X. Ge, et al., “5G ultra-dense cellular networks,” IEEE Wireless Commun., vol. 23, no. 1, pp. 72-79, 2016.
  • [4] L. Wang, “Wireless powered dense cellular networks: how many small cells do we need?” IEEE J. Sel. Areas Commun., vol. 35, no. 9, pp. 2010-2024, 2017.
  • [5] J. Chen, et al., “Coverage and handoff analysis of 5G fractal small cell networks,” IEEE Trans. Wireless Commun., vol 18, no. 2, pp. 1263-1276, 2019.
  • [6] X. Ge, et al., “Small cell networks with fractal coverage characteristics,” IEEE Trans. Commun., vol. 66, no. 11, pp. 5457-5469, 2018.
  • [7] J. Zhang, et al., “Femtocells-technologies and deployment,” Wiley, 2010.
  • [8] D. Lopez, et al., “OFDMA femtocells: a roadmap on interference avoidance,” IEEE Commun. Mag., vol. 47, no 9, pp. 41-48, 2009.
  • [9] A. H. Jafari, et al., “Performance analysis of dense small cell networks with practical antenna heights under Rician fading,” IEEE Access, vol. 6, pp. 9960-9974, 2017.
  • [10] X. Ge, et al., “Joint optimization of computation and communication power in multi-user massive MIMO systems,” IEEE Trans Wireless Commun., Vol. 17, No. 6, pp. 4051-4063, 2018.
  • [11] J. Hoydis, et al., “Massive MIMO in the UL/DL of cellular networks: How many antennas do we need?” IEEE J. Sel. Area Commun., vol. 31, no. 2, pp. 160-171, 2013.
  • [12] E. G. Larsson, et al., “Massive MIMO for next generation wireless systems,” IEEE Commun. Mag. vol. 52, no. 2, pp. 186-195, 2014.
  • [13] L. Lu, et al., “An overview of massive MIMO: Benefits and challenges,” IEEE J. Sel. Topics Sig. Proc., vol. 8, no. 5, pp. 742-758, 2014.
  • [14] T. Nam and T. A. Pardo, “Conceptualizing smart city with dimensions of technology, people, and institutions,” in Proc. Annual International Conference on Digital Government Research, pp. 282-291, 2011.
  • [15] Ranplan Wireless Network Design,“Method for evaluating the building wireless performance,” Patent, 2019.
  • [16] J. Zhang, et al., “Wireless energy efficiency evaluation for buildings under design based on analysis of interference gain,” IEEE Trans. Veh. Tech., vol. 69, no. 6, pp. 6310-6324, 2020.
  • [17] D. Yagaanbuant, et al., “Some design issues of high-rise buildings,” in Proc. IFOST, 2007.
  • [18] L. Pérez-Lombard, et al., “A review on buildings energy consumption information,” Energy and buildings, vol. 40, no. 3, pp. 394-398, 2008.
  • [19] CIBSE, “TM54: evaluating operational energy performance of buildings at the design stage,” CIBSE, 2013.
  • [20] D. Feng, et al., “A survey of energy-efficient wireless communications,” IEEE Commun. S. & Tutor., vol. 15, no. 1, pp. 167-178, 2013.
  • [21] C. L. I, et al., “Energy-efficient 5G for a greener future,” Nat. Electron., vol. 3, pp. 182-184, 2020.
  • [22] A. Abrol and R. K. Jha, “Power optimization in 5G Networks: a step towards green communication,” IEEE Access, vol. 4, pp. 1355-1374, 2016.
  • [23] H. Zheng, et al., “Exact line-of-sight probability for channel modelling in typical indoor environments,” IEEE Ant. Wireless Prop. Lett., vol. 17, no. 7, pp. 1359-1362, 2018.
  • [24] W. Yang, et al., “Line-of-sight probability for channel modeling in 3D indoor environments,” IEEE Ant. Wireless Prop. Lett., vol. 19, no. 7, pp. 1182-1186, 2020.
  • [25] W. Yang, et al., “Machine learning based indoor line-of-sight probability prediction,” in Proc. IEEE ISAP, 2019.
  • [26] T. Bai, et al., “Coverage and rate analysis for millimeter-wave cellular networks,” IEEE Trans. Wireless Commun., vol. 4, no. 2, pp. 1100-1114, 2015.
  • [27] M. K. Mu¨\rm\ddot{u}ller, et al., “Analyzing wireless indoor communications by blockage models,” IEEE Access, vol. 5, pp. 2172-2186, 2017.
  • [28] K. Haneda, et al., “Spatial coexistence of millimeter-wave distributed indoor channels,” in Proc. IEEE VTC-Spring, 2015.
  • [29] L. Azpilicueta, et al., “A hybrid ray launching-diffusion equation approach for propagation prediction in complex indoor environments,” IEEE Ant. Wireless Prop. Lett., vol. 16, pp. 214-217, 2017.
  • [30] P. H. Tseng, et al.,“Ray-tracing-assisted fingerprinting based on channel impulse response measurement for indoor positioning,” IEEE Trans. Ins. Mea., vol. 66, no. 5, pp. 1032-1045, 2017.
  • [31] F. Casino, et al., “Optimized wireless channel characterization in large complex environments by hybrid ray launching-collaborative filtering approach,” IEEE Ant. Wireless Prop. Lett., vol. 16, pp. 780-783, 2017.
  • [32] Z. Lai, et al., “Modelling the mmWave channel based on intelligent ray launching model,” in Proc. Eucap, 2015.
  • [33] W. Yang, et al., “Indoor measurement based verification of ray launching algorithm at the Ka-band,” in Proc. URSI GASS, 2020.
  • [34] Z. Lai, et al., “Intelligent ray launching algorithm for indoor scenarios,” Radioengineering, vol. 20, no. 2, 2011.
  • [35] EUROSTARS Project 11088 BuildWise: The development of a software tool to design buildings with tailored wireless performance (Ranked the 16th among 331 applications), [Online]. https://www.eurostars-eureka.eu/project/id/11088, 2017.
  • [36] X. Zhang and J. Andrews. “Downlink cellular network analysis with multi-slope path loss models,” IEEE Trans. Commun., vol. 63, no. 5, pp. 1881-1894, 2015.
  • [37] C. Sommer and F. Dressler, “Using the right two-ray model? a measurement based evaluation of PHY models in VANETs,” in Proc. ACM MobiCom, 2011.
  • [38] 3GPP, “Study on channel model for frequencies from 0.5 to 100 GHz,” 3GPP TR 38.901 version 14.0.0 Release 14, 2017.
  • [39] U. Madhow, et al.,“Distributed massive MIMO: algorithms, architectures and concept systems,” in Proc. IEEE ITA, 2014.
  • [40] S. Fu, et al.,“Distributed transmission scheduling and power allocation in CoMP,” IEEE Sys. J., vol. 12, no. 4, pp. 3096-3107, 2018.
  • [41] E. Pateromichelakis et al., “On the evolution of multi-cell scheduling in 3GPP LTE/LTE-A,” IEEE Commun. Surveys Tuts., vol. 15, no. 2, pp. 701-717, 2013.
  • [42] J. G. Andrews, et al., “What will 5G be?” IEEE J. Sel. Area Commun., vol. 32, no. 6, pp. 1065-1082, 2014.
  • [43] D. S. Baum, et al., “IST-WINNER D5.4, final report on link and system level channel models,” https://www.researchgate.net/publication/229031750_IST-2003-507581_WINNER_D5_4_v_14_Final_Report_on_Link_Level_and_System_Level_Channel_Models, 2005.