# Line-Of-Sight Interference Statistics and User Scheduling with 3D Angular Separation

Authors: Yawei Yu; Jianhua Zhang; Lei Tian
Journal: International Journal of Antennas and Propagation
 (2018)

---

## Abstract

We focus on the interference statistical metrics of matched filter transmission over line-of-sight (LoS) downlink channels. The exact moments of the interference are derived, which can be used to make approximations for user rates and gain insights into the effects of system parameters such as angle spread (AS), antenna numbers, and array spacings. As the interference decreases with larger AS, we propose a user scheduling scheme by grouping users into different subspaces with corresponding orthogonal fraction of resource blocks, aiming to maximize the 3D angular separation among multiusers within each group; thus, a lower interference and higher sum rate would be expected. By performing this user scheduling scheme with 3D angular separation, a large sum rate gain is observed, which provides us a promising method in interference control and management.


---

## Body

## 1. Introduction

The growing importance of massive MIMO [1], small cell [2], and millimeter wave [3] communications is motivating an increased interest in line-of-sight (LoS) channels [4] and simple transmission schemes such as matched filter (MF) precoding [5]. However, interference will arise in such multiuser systems with MF transmission [6], which will degrade the system performance [7], thus drawing the researchers’ attention [8].Considering a single-cell massive MIMO system whereK single-antenna users are evenly distributed around the base station (BS) in the horizontal plane, the favorable propagation conditions of LoS channels, which are defined as mutual orthogonality among the vector-valued channels to the terminals with no mutual interference, are well investigated in [9]. As the mutual orthogonality among LoS channels is closely related to the antenna array at BS, the impact of different antenna parameters on asymptotic orthogonality, including the array geometries [10], antenna spacings [11], and antenna numbers [12], was studied. Using a uniform rectangular array (URA) at BS for the LoS MIMO system, results in [13] show that the ergodic rate is maximized for antenna spacing equal to integer multiples of one-half wavelength while the achievable sum rate with other antenna arrays can be found in [14]. In LoS channels, the phase of the received signal is determined by the relative location of the user (to BS); the work in [15] demonstrates that the user interference is partially determined by the LoS angle of arrival of signals. The deterministic dependence of interference on locations poses a significant challenge in modelling the interference distribution; the authors in [16] show that the main lobe distribution of LoS interference can be approximated by a beta mixture in the uplink LoS channel. As interference distribution can be used to compute signal-to-noise-plus-interference ratio (SINR) metrics such as outage probability and average throughput, prior work on LoS interference modeling is summarized in [17] and the interference model based on SINR is proposed and emphasized. The authors in [18] propose a statistical-eigenmode space-division multiple-access downlink transmission scheme based on the ergodic SINR and present the achievable sum rate of the multiuser system. This work also studies, for a different reason, the interference variables considered here.As the sum rate of the multiuser MIMO system with interference is closely related to the SINR at each user, optimizing the SINR is key to improve the system performance [19]. Framework on transceiver (precoder and decoder) designing for interference alignment (IA) is presented [20], and different IA schemes based on resource allocation in multiuser MIMO systems are summarized [21]. By grouping users with different mobilities to have different coherence times, the interference signal only occupies limited dimensions and the degrees of freedom for the LoS interference channel are improved due to this blind IA scheme [22]. Furthermore, IA via a time-indexed interference graph for LoS channels is studied in [23], a maximal independent set is presented via dynamic programing, and orthogonal interference subspaces are designed to combat the interference at each user. Considering the MIMO downlink channel in a multiuser system, an efficient suboptimal algorithm of user scheduling is given by assigning the coordinates of transmission space to different users, aiming to achieve the best performance in terms of the sum rate throughput [24]. By allocating various resource blocks (e.g., time, frequency, and power) to users in different groups, the multiuser interference would remain within acceptable limits [25].As we can summarize, all of the above work is closely related to the statistic characteristics of LoS interference; however, the derivation of a multiuser interference is difficult due to integrations of complicated components such as different elevation angle distributions and exponential antenna radiation patterns; thus, work on exact expressions of the interference moments is rare. Furthermore, conventional work all assumes that multiusers are evenly or randomly distributed around BS in the azimuth domain with no consideration of the elevation counterpart, which is not that realistic especially for 3D propagation scenarios with multiusers in high-rise buildings. For a 3D LoS propagation channel, the interference model is lacking and the effects of system parameters including the angle spread (AS), antenna spacing, and antenna numbers on the interference still remain unknown. Till now, there are various IA schemes such as resource allocations in the time, frequency, and power domains; the utilization of the degree of freedom in 3D space is rarely reported. The overall performance of multiuser MIMO systems is a complex joint multiobjective optimization problem since many variables and parameters have to be optimized, including the selection of users, the number and interval spacing of antennas, power allocation, and transmission technique. Based on the prior work, we make the following contributions:(1)
Considering the LoS downlink channel with a simple MF transmission scheme and multiusers are distributed in the 3D space, we derive the exact expressions of the interference moments and present the simulative interference results under different antenna numbers, antenna spacings, and AS values. As usual, the multiusers are assumed as evenly distributed around the BS in the azimuth domain and von Mises (VM) distribution in the elevation counterpart [26].
(2)
Comparative results of the interference against AS, antenna numbers, and antenna spacings demonstrate that the averaged interference decreases rapidly when the elevation angle spread (EAS) increases (or the antenna number decreases, antenna spacing increases), which enables us to reduce the interference once users are scheduled to maintain decent angular separations (or finite antenna activation, antenna configuration optimization). While the mean interference level falls with increasing AS, there remains large variability (see the results of interference variance), especially in the super dense area withEAS<10°, which poses a severe challenge to simple MF processing.
(3)
Except IA schemes achieved by user scheduling in time, frequency, and power domains as mentioned above, we present the user scheduling via 3D space utilization as well. On the criteria of maximizing the angular separations among users within each group, multiusers are divided into different groups and occupy corresponding orthogonal fraction of wireless resource blocks. By performing the user scheduling scheme via 3D angular separation, a large sum rate gain is obtained and will be further enhanced by larger EAS and antenna spacing.
(4)
We conclude that by performing the user scheduling scheme via 3D angular separation, a significant system sum rate gain can be achieved, which ensures the promising future of 3D space utilization in interference management.The remainder of this article is outlined as follows: in Section2, we give the system model and performance metrics. In Sections 3 and 4, we derive the interference moments and present the results of interference against AS and antenna spacing. In Section 5, the user scheduling scheme is put forwarded along with the simulative results of the system sum rate. Then, conclusions are drawn in Section 6.


## 2. System Model and Performance Metrics

Given the typical URA at BS, as shown in Figure1, there are M vertical elements spaced by d1 wavelength and N horizontal elements with d2 wavelength spacing. For simplicity, we assume K single-antenna users (UEs) being served in a pure LoS environment within a single cell and the LoS angle of UE k is denoted by (φk,θk). The channel impulse response (CIR) between UE k and the URA, Hk, is given by [27]
(1)Hk=βkHh,k⊗Hv,k,where
(2)Hv,k=1,e−juk,e−j2uk,…,e−jM−1ukT,Hh,k=1,e−jvk,e−j2vk,…,e−jN−1vk,uk=2πλd1cosθk,υk=2πλd2sinθkcosϕk.Figure 1
Multiuser LoS interference with URA at BS.βk is the link gain, and λ is the wavelength.We next rewrite theM×N matrix Hk as a vector hk in a dimension of 1×MN as follows:
(3)hk=βkak,(4)ak=Hh,k,e−jukHh,k,e−j2ukHh,k,…,e−jM−1ukHh,k.Thus, the received downlink signal at UEk is obtained
(5)yk=hkx+nk,where x is the transmitted signal and nk∼CN0,σz2 is the additive white Gaussian noise with mean value 0 and variance σz2.Assuming linear precoding at the BS,(6)x=Gd,where G=g1,g2,…,gK is a MN×K precoding matrix and d=d1,d2,…,dKT contains the transmitted data. If matched filter precoding is employed at the BS, then gk=hk+MNβk−1/2, and the SINR at UE k, assuming Edκ2=1, is given by
(7)SINRk=hkgk2∑j≠kKhkgj2+σz2=MNβkβk/MN∑j≠kKakaj+2+σz2.In (7), we observe that the intracell interference to UE k from the transmitted signal to UE j is βkIkj/MN, where
(8)Ikj=akaj+2.Thus, the total intracell interference to UEk is βkIk/MN, where
(9)Ik=∑j≠kKIkj.The desired signal power isMNβk; the SINR and the rate of UE k can be expressed as
(10)SINRk=MNβkβk/MNIk+σz2,Rk=log21+MNβkβk/MNIk+σz2.We note that the moments can be used to estimate the mean rate through the delta method [28], where the first- and second-order center moments of Ik are included:
(11)ERk≈log21+1Eδ+γ,where
(12)δ=1MNIkNM+σ2βk,(13)γ=Varδ2ln21Eδ2−1Eδ+12.Clearly, the performance is critically dependent on the distribution of the interferenceIk. Hence, we will derive the moments of Ik, that is, EIk and VarIk.The other important variable governing the performance is the distribution of the LoS angles,ϕk,θk,k=1,2,…,K. We assume the uniform distribution in azimuth domain with ϕk∼U−π,π and the von Mises distribution in elevation domain with θk∼VMκ,μ, where the density function is
(14)fθ=eκcosθ−μ2πI0κ,0≤θ≤π.Here, κ is a concentration parameter, μ is a location parameter, and I0· is the modified Bessel function of the first kind in the zeroth order.


## 3. Interference Moments

From (11)–(13), we observe that EIk and VarIk are required to compute the sum rate. Hence, we compute the first two moments of Ii based on the following lemma.Lemma 1.
Whenϕ∼U−π,π and θ∼VMκ,π/2 with fθ in (14), the approximation
(15)eκcosθ−π/2≈eκ−1sinθL+1,with the approximation accuracy validated in Figure 2, enables us to obtain
(16)Ee−jm−puke−jn−quk≜Fm−p,n−q,κ≈A+B+bJ0−2πdzm,pC+2b∑r=1RjrJr−2πdzm,pD,where
(17)A=aπJ0−2πdzm,pJ02πdyn,q,(18)B=2a∑r=1RjrJr−2πdzm,pπcosr2πJ−r/2πdyn,qJr/2πdyn,q,when L=2L0 is even in (15),
(19)C=∑l=0L0k0πcoslπJ−lπdyn,qJlπdyn,q,(20)D=∑l=0L0kl2πcosl1πJ−l1πdyn,qJl1πdyn,q+cosl2πJ−l2πdyn,qJl2πdyn,q.when L=2L1−1 is odd in (15),
(21)C=∑l=1L1k0πsinl−12πJ−l−1/2πdyn,qJl−1/2πdyn,q,(22)D=∑l=1L1kl2πsinl3πJ−l3πdyn,qJl3πdyn,q+sinl4πJ−l4πdyn,qJl4πdyn,q.Figure 2
Approximation accuracy ofeκsinθ/eκ≈eκ−1sinθ2L+1/eκ.The parameters mentioned in (16)–(22) are given below:
(i)
dyn,q and dzm,p are defined as
(23)dyn,q=n−qd2λ,dzm,p=m−pd1λ.
(ii)
a and b are defined as
(24)a=12πI0κ,b=ek−12πI0κ.
(iii)
l1, l2, l3, and l4 are defined as
(25)l1=2l+r2,l2=2l−r2,l3=2l−1+r2,l4=2l−1−r2.
(iv)
kl are coefficients of sinθL expansion in [29, eq. (1.320.1)]:when L=2L0 is even,
(26)sinθL=∑l=0L0klcos2lθ,(27)k0=122L0L02L0kl=−1l22L0−1L0−l2L0,L=1,2,…L0.when L=2L1−1 is odd,
(28)sinθL=∑l=1L1klsin2l−1θ,(29)kl=−1l−122L1−2L1−l2L1−1,L=1,2,…,L1.Proof 1.
See AppendixA.
Since the double integralEe−jm−puke−jn−qvk occurs throughout the moment derivations, we use the short notation, Fm−p,n−q,κ, defined in (16). From Lemma 1, the following results can be obtained.Theorem 1.
(30)EIk=K−1∑m,n,p,qFm−p,n−q,κ2,where ∑m,n,p,q denotes ∑m=0M−1∑n=0N−1∑p=0M−1∑q=0N−1.Proof 2.
See AppendixB.Theorem 2.
(31)EIk2=K−12∑p1q1m1n1∑p2q2m2n2Fm1+m2−p1−p2,n1+n2−q1−q2−q2,κ⋅F−m1−p1,−n1−q1,κ⋅F−m2−p2,−n2−q2,κ.Proof 3.
See AppendixC.
From Theorems1 and 2, the variance of Ik follows from
(32)VarIk=EIk2−EIk2.Note that, whenκ=0, we have the special case where elevation angle is uniformly distributed with θ∼U0,π.


## 4. Simulative Results of Interference and Analysis

### 4.1. Numerical and Simulation Results: Moments

The interference is closely related to the multiuser angle distribution in both elevation and azimuth domains. By following (16) in Lemma 1 and Theorem 1, we present the interference variation against elevation angle spread under various antenna spacings in Figure 3 when considering an even distribution of azimuth angle over [0,2π].Figure 3
Mean interferenceEIk against elevation angle spread: K=5.Specifically, results are shown for the rectangular antenna array with 100 (10 × 10) and 36 (6 × 6) elements at BS with antenna spacingd=0.25λ, d=0.5λ, and d=1λ (d1=d2=d) when K=5. For all cases, the simulative results of the approximate exact moments fit quite well with the true value; thus, the derivation of Lemma 1 based on the simple approximation in (15) provides us a general method to calculate the mean interference with high approximation accuracy. Obvious features of the results include a drop in the mean interference as EAS increases (κ decreases), as antenna spacing increases, and as antenna number at BS decreases. Furthermore, we note that for small EASs, the interference is massive and stays stable, which should account for the tiny variation of the term e−jm−puke−jn−qvk in (16) for small EAS (<10°) in the super dense area. When the EAS increases from 10° to 40°, defined as dense area, and the mean interference EIk decreases rapidly, then it reaches the bottom gradually with a long flat tail in the sparse area with large EASs (40°∼90°). Moreover, a value of MN2 makes the interference power equal to the signal power, and this is exceeded by many of the scenarios. Clearly, either low numbers of users or carefully separated users (via scheduling) or both are essential for the satisfactory performance.


### 4.2. Numerical and Simulation Results: Variance

In Figure4, we plot the spread of the interference by including levels which are one standard derivation above and below the mean, that is, EIk±stdIk, where stdIk=VarIk. The same parameters such as K=5 and a rectangular antenna array with 100 (10 × 10) elements at BS with d1=0.5λ and d2=0.5λ are used in the simulation.Figure 4
Interference variationEIk±stdIk against elevation angle spread: K=5 and a rectangular antenna array with 100 (10 × 10) elements at BS with d1=0.5λ and d2=0.5λ.While the mean interference level falls with increasing EAS, there remains large variability, especially in the super dense area with small EAS (<10°), where the one standard derivation level is approximately 2.5 times greater than the mean interference, which poses a severe challenge for the interference alignment in super dense area. With the EAS increases from 10° to 40°, the one standard derivation level falls in the same shape as the mean interference does, which indicates a stable interference variance in the dense area. At high EAS in the sparse area, the mean interference keeps a long flat tail and the one standard derivation level gets closer to the mean; the smaller spread of interference would contribute to a more stable system performance. To be noted, the uniform distribution of elevation angle is achieved when the EAS=90° (κ=0).


## 4.1. Numerical and Simulation Results: Moments

The interference is closely related to the multiuser angle distribution in both elevation and azimuth domains. By following (16) in Lemma 1 and Theorem 1, we present the interference variation against elevation angle spread under various antenna spacings in Figure 3 when considering an even distribution of azimuth angle over [0,2π].Figure 3
Mean interferenceEIk against elevation angle spread: K=5.Specifically, results are shown for the rectangular antenna array with 100 (10 × 10) and 36 (6 × 6) elements at BS with antenna spacingd=0.25λ, d=0.5λ, and d=1λ (d1=d2=d) when K=5. For all cases, the simulative results of the approximate exact moments fit quite well with the true value; thus, the derivation of Lemma 1 based on the simple approximation in (15) provides us a general method to calculate the mean interference with high approximation accuracy. Obvious features of the results include a drop in the mean interference as EAS increases (κ decreases), as antenna spacing increases, and as antenna number at BS decreases. Furthermore, we note that for small EASs, the interference is massive and stays stable, which should account for the tiny variation of the term e−jm−puke−jn−qvk in (16) for small EAS (<10°) in the super dense area. When the EAS increases from 10° to 40°, defined as dense area, and the mean interference EIk decreases rapidly, then it reaches the bottom gradually with a long flat tail in the sparse area with large EASs (40°∼90°). Moreover, a value of MN2 makes the interference power equal to the signal power, and this is exceeded by many of the scenarios. Clearly, either low numbers of users or carefully separated users (via scheduling) or both are essential for the satisfactory performance.


## 4.2. Numerical and Simulation Results: Variance

In Figure4, we plot the spread of the interference by including levels which are one standard derivation above and below the mean, that is, EIk±stdIk, where stdIk=VarIk. The same parameters such as K=5 and a rectangular antenna array with 100 (10 × 10) elements at BS with d1=0.5λ and d2=0.5λ are used in the simulation.Figure 4
Interference variationEIk±stdIk against elevation angle spread: K=5 and a rectangular antenna array with 100 (10 × 10) elements at BS with d1=0.5λ and d2=0.5λ.While the mean interference level falls with increasing EAS, there remains large variability, especially in the super dense area with small EAS (<10°), where the one standard derivation level is approximately 2.5 times greater than the mean interference, which poses a severe challenge for the interference alignment in super dense area. With the EAS increases from 10° to 40°, the one standard derivation level falls in the same shape as the mean interference does, which indicates a stable interference variance in the dense area. At high EAS in the sparse area, the mean interference keeps a long flat tail and the one standard derivation level gets closer to the mean; the smaller spread of interference would contribute to a more stable system performance. To be noted, the uniform distribution of elevation angle is achieved when the EAS=90° (κ=0).


## 5. User Scheduling with 3D Angular Separation

As indicated by the results in Figures3 and 4 that the interference will decrease rapidly as EAS increases, in this section, we focus on the user scheduling scheme via 3D angular separation and put forward the user scheduling algorithm to divide multiusers into different independent groups, aiming to maximize the mean angular separation within each group and achieve satisfactory system performance, that is, high sum rate.


### 5.1. User Scheduling Algorithm with 3D Angular Separation

Assume there areK LoS users in a single cell which are expected to be grouped into G independent groups based on the criteria of maximizing the angular separation of users within each group, which turn out to be a clustering problem. For each LoS user i, the location is characterized by (di,θi,φi), where di, θi, and φi denote the distance, the LoS elevation angle, and the LoS azimuth angle from user i to the BS, respectively. The azimuth angle of multiuser, φ, is assumed as a uniform distribution where ϕ∼U0,2π, while the elevation angle, θ, follows the typical VM distribution where θ∼VMμ,κ. Multiusers are located evenly between the reference distance for the antenna far-field, d0, and the cell radius, R. The path loss of the signal from the BS to user i is given by [30, eq (2.28)]
(33)Pi=PtiPad0dir,d>d0,where Pti is the signal transmission power to user i, Pa is a unitless constant which depends on the antenna characteristics and the average channel attenuation, and r denotes the path loss exponent. Performing the transformation from the spheric coordinate to the Cartesian coordinate, the location vectors of users i and j are given by
(34)xi=Pisinθicosϕi,Pisinθisinϕi,Picosθi,xj=Pjsinθjcosϕj,Pjsinθjsinϕj,Pjcosθj.Then, the angular separation between usersi and j, ∆ij, is calculated as
(35)△ij=arccosxi⋅xjT=xi⋅xjTxixj.Thus, we could obtain the angular separation matrix ∆ as follows:
(36)Δij=Δij,i=1,2,…,K;j=1,2,…,K,which will be used as the input of the user scheduling algorithm below:Algorithm 1: User scheduling with 3D angular separation.
InputΔ: angular separation matrix of all K LoS users
G: The intended user group number
OutputΛ1,Λ2,…,ΛG: user scheduling into G groups.
1:if G=1 then
2: No user scheduling (K users in the same group)
3:else
4:i1,i2=argminΔi1,i2    (minimum angular separation)
5:Λ1=i1,Λ2=i2 (user i1→group Λ1,i2→group Λ2)
6:IG=Λ1∪Λ2=i1,i2       (the set of grouped users)
7:Δi1,i2=inf                (update ∆)
8:for g=3:G do       (search the initial user for group 3 to G)
9:IL=1,2,…,K−IG         (users left to be grouped)
10:Si=sumΔi,IG        for each user i∈IL
11:il=argminSi,i∈IL (minimum sum angle separation)
12:Λl=il,IG=IG∪Λl      (new group with user il)
13:Δil,IG=inf                (update ∆)
14:end for                (obtain IG=i1,i2,…,iG)
15:for u=1:K−G do
16:i1,i2=argminΔi1,i2       (users to be grouped)
17:g1,g2=argmaxmeanΔΛg1,i1+ΔΛg2,i2
18:Λg1=Λg1∪i1,Λg2=Λg2∪i2  (user scheduling)
19:Δg1,i1=inf,Δg2,i2=inf          (update ∆)
20:IG=IG∪i1,i2             (update IG)
21:if lengthIG=K then
22:break          (ensure all K users are grouped)
23:end if
24:end for
25:end ifThe following are the parameters:Λk denotes the users selected into the kth group, IG denotes the users which were divided while IL denotes the users left to be grouped, Si denotes the sum of the angular separation between user i and all users in group G, and i1 and i2 denote the users to be grouped in the current step with the target group number g1 and g2, respectively.In Algorithm1, the users in G groups are initialized by steps 4∼14 where G users with the minimum angular separation among all K users are selected into G different groups, which needs a calculation of G×KK−1/2 complexity. Then, the left K−G users are allocated to the existing G groups in steps 15∼24 in an order of separating the two nearest users firstly (step 16). Note that, there are GG−1/2 cases to allocate users i1 and i2 to two different groups. For each case, the averaged angular separation of multiuser is calculated. By maximizing the averaged angular separation (step 17), users i1 and i2 are distributed to groups Λg1and Λg2, respectively; thus, the calculation complexity is K−GGG−1/2. Finally, by performing the iterative suboptimal user scheduling algorithm via 3D angular separation, these K users would be divided into expected G groups and the total calculation complexity is OG/2K2−GG−2/2K−G2G−1/2, which can be simplified as OG/2K2 as G≪K usually.In the following part of this section, the system sum rate will be calculated; then, the comparative results and analysis will be presented under different expected group numbers, that is,G=1,2,3,4,5. To be noted, when G=1, all users are in the same group and no user scheduling is performed. When G>1, as a reference, we present the results of the system sum rate for user random scheduling as well, which means K users are randomly allocated to G groups without selection.


### 5.2. System Sum Rate Calculation

Based on Algorithm1, users will be divided into different groups with the largest angular separation to avoid the potential interference. As users in different groups will be served in orthogonal resource blocks, the mutual interference between users in different groups is assumed 0; thus, the multiuser sum rate can be calculated as
(37)R¯=∑g=1G∑k∈Λglog21+MNβkβk/MNIkg+σz2,where Ikg denotes the interference to user k from other users in group Λg.We assume that the angles of all users are independent; the first momentEIk in (9) can be simplified as
(38)EIk=∑j≠kKEIkj=K−1EIkj.Thus, we would obtain
(39)EIkg=∑j≠k,j∈ΛgLgEIkj=Lg−1EIkj.


### 5.3. Simulative Results and Analysis

In the simulation, we assume the same antenna array, the URA in Figure1 with 100 elements (10 × 10) and half wavelength antenna spacing in both elevation and azimuth domains (d1=0.5 and d2=0.5), is used at BS. K=20 users are distributed evenly around the BS in the azimuth domain while the elevation angles obey a von Mises distribution in the LoS MIMO system. The comparative results of the system sum rate between user scheduling in Algorithm 1 and random scheduling without user selection are presented in Figure 5 under different EASs, that is, EAS=10°,30°,50°, which would provide insights on the effects of user scheduling via 3D angular separation on the system performance.Figure 5
Multiuser sum rate under different user group values. (a) Multiuser sum rate under different user group values:K=20,Φ∼U0,2π, and θ∼VM0,8.52 with EAS=10°. (b) Multiuser sum rate under different user group values: K=20,Φ∼U0,2π, and θ∼VM0,1.48 with EAS=30°. (c) Multiuser sum rate under different user group values: K=20,Φ∼U0,2π, and θ∼VM0,0.14 with EAS=50°.
(a)
(b)
(c)From the results in Figure5, an obvious feature is that a large sum rate gain is obtained by performing the user scheduling via 3D angular separation (solid red curves) as compared to the sum rate of no scheduling (solid red curve with G=1) and random scheduling (dotted dark curves). Furthermore, the more groups the multiusers are expected to be divided into, the higher the sum rate is observed, and this sum rate gain will be further enhanced. Taking the results in Figure 5(a) as an example, when the expected group number G increases from 1 to 2, 3, 4, and 5, the 0.5 CDF point sum rate gains harvested by user scheduling scheme are 3.9, 7.3, 10.3, and 13.3 bps/Hz, respectively, under EAS=10°. Once the multiusers are scheduling into different independent groups, the intergroup interference would be eliminated, which contributes to a lower total system interference and higher sum rate. Meanwhile, by performing the user scheduling algorithm with maximal angular separation among users within each group, the intragroup interference will be lowered; thus, the system sum rate is further improved.Increasing the EAS from 10° to 30° and 50° where users are better separated in the 3D space, the corresponding results of the system sum rate are illustrated in Figures5(b) and 5(c), respectively. As compared to the results in Figure 5(a), the system sum rate is improved by larger EAS as expected. A more specific result of the sum rate gain against EAS is presented in Table 1. As we can see, a large sum rate gain is obtained by user scheduling with 3D angular separation, and the gain ratio decreases as the EAS increases. For small EAS (10°), the sum rate gain can even be 3.5 times larger than that of no scheduling case, which highlights the importance of user scheduling for interference alignment in super dense area.Table 1
The sum rate gain and gain ratio (CDF=0.5) harvested by user scheduling with 3D angular separation when G=5 where Lg denotes the user number in group Λg while the moment of interference, EIkj, can be found in (B.3) in the Appendix B.
User scheduling (G=5) versus
No scheduling (G=1)
Random scheduling (G=5)
Sum rate gain (bps/Hz)/gain ratio
EAS=10°
14.0 bps/Hz/345%
4.6 bps/Hz/34%
EAS=30°
16.9 bps/Hz/114%
7.7 bps/Hz/32%
EAS=50°
13.4 bps/Hz/55%
6.1 bps/Hz/19%Although the above observations are for multiuser interference in a single cell, similar results of significant sum rate gain due to the proposed user scheduling algorithm can be expected for multiple cells in the future work. As we can conclude, in environments with considerable angular diversity, user scheduling based on angular separations will be an effective tool in interference alignment.


## 5.1. User Scheduling Algorithm with 3D Angular Separation

Assume there areK LoS users in a single cell which are expected to be grouped into G independent groups based on the criteria of maximizing the angular separation of users within each group, which turn out to be a clustering problem. For each LoS user i, the location is characterized by (di,θi,φi), where di, θi, and φi denote the distance, the LoS elevation angle, and the LoS azimuth angle from user i to the BS, respectively. The azimuth angle of multiuser, φ, is assumed as a uniform distribution where ϕ∼U0,2π, while the elevation angle, θ, follows the typical VM distribution where θ∼VMμ,κ. Multiusers are located evenly between the reference distance for the antenna far-field, d0, and the cell radius, R. The path loss of the signal from the BS to user i is given by [30, eq (2.28)]
(33)Pi=PtiPad0dir,d>d0,where Pti is the signal transmission power to user i, Pa is a unitless constant which depends on the antenna characteristics and the average channel attenuation, and r denotes the path loss exponent. Performing the transformation from the spheric coordinate to the Cartesian coordinate, the location vectors of users i and j are given by
(34)xi=Pisinθicosϕi,Pisinθisinϕi,Picosθi,xj=Pjsinθjcosϕj,Pjsinθjsinϕj,Pjcosθj.Then, the angular separation between usersi and j, ∆ij, is calculated as
(35)△ij=arccosxi⋅xjT=xi⋅xjTxixj.Thus, we could obtain the angular separation matrix ∆ as follows:
(36)Δij=Δij,i=1,2,…,K;j=1,2,…,K,which will be used as the input of the user scheduling algorithm below:Algorithm 1: User scheduling with 3D angular separation.
InputΔ: angular separation matrix of all K LoS users
G: The intended user group number
OutputΛ1,Λ2,…,ΛG: user scheduling into G groups.
1:if G=1 then
2: No user scheduling (K users in the same group)
3:else
4:i1,i2=argminΔi1,i2    (minimum angular separation)
5:Λ1=i1,Λ2=i2 (user i1→group Λ1,i2→group Λ2)
6:IG=Λ1∪Λ2=i1,i2       (the set of grouped users)
7:Δi1,i2=inf                (update ∆)
8:for g=3:G do       (search the initial user for group 3 to G)
9:IL=1,2,…,K−IG         (users left to be grouped)
10:Si=sumΔi,IG        for each user i∈IL
11:il=argminSi,i∈IL (minimum sum angle separation)
12:Λl=il,IG=IG∪Λl      (new group with user il)
13:Δil,IG=inf                (update ∆)
14:end for                (obtain IG=i1,i2,…,iG)
15:for u=1:K−G do
16:i1,i2=argminΔi1,i2       (users to be grouped)
17:g1,g2=argmaxmeanΔΛg1,i1+ΔΛg2,i2
18:Λg1=Λg1∪i1,Λg2=Λg2∪i2  (user scheduling)
19:Δg1,i1=inf,Δg2,i2=inf          (update ∆)
20:IG=IG∪i1,i2             (update IG)
21:if lengthIG=K then
22:break          (ensure all K users are grouped)
23:end if
24:end for
25:end ifThe following are the parameters:Λk denotes the users selected into the kth group, IG denotes the users which were divided while IL denotes the users left to be grouped, Si denotes the sum of the angular separation between user i and all users in group G, and i1 and i2 denote the users to be grouped in the current step with the target group number g1 and g2, respectively.In Algorithm1, the users in G groups are initialized by steps 4∼14 where G users with the minimum angular separation among all K users are selected into G different groups, which needs a calculation of G×KK−1/2 complexity. Then, the left K−G users are allocated to the existing G groups in steps 15∼24 in an order of separating the two nearest users firstly (step 16). Note that, there are GG−1/2 cases to allocate users i1 and i2 to two different groups. For each case, the averaged angular separation of multiuser is calculated. By maximizing the averaged angular separation (step 17), users i1 and i2 are distributed to groups Λg1and Λg2, respectively; thus, the calculation complexity is K−GGG−1/2. Finally, by performing the iterative suboptimal user scheduling algorithm via 3D angular separation, these K users would be divided into expected G groups and the total calculation complexity is OG/2K2−GG−2/2K−G2G−1/2, which can be simplified as OG/2K2 as G≪K usually.In the following part of this section, the system sum rate will be calculated; then, the comparative results and analysis will be presented under different expected group numbers, that is,G=1,2,3,4,5. To be noted, when G=1, all users are in the same group and no user scheduling is performed. When G>1, as a reference, we present the results of the system sum rate for user random scheduling as well, which means K users are randomly allocated to G groups without selection.


## 5.2. System Sum Rate Calculation

Based on Algorithm1, users will be divided into different groups with the largest angular separation to avoid the potential interference. As users in different groups will be served in orthogonal resource blocks, the mutual interference between users in different groups is assumed 0; thus, the multiuser sum rate can be calculated as
(37)R¯=∑g=1G∑k∈Λglog21+MNβkβk/MNIkg+σz2,where Ikg denotes the interference to user k from other users in group Λg.We assume that the angles of all users are independent; the first momentEIk in (9) can be simplified as
(38)EIk=∑j≠kKEIkj=K−1EIkj.Thus, we would obtain
(39)EIkg=∑j≠k,j∈ΛgLgEIkj=Lg−1EIkj.


## 5.3. Simulative Results and Analysis

In the simulation, we assume the same antenna array, the URA in Figure1 with 100 elements (10 × 10) and half wavelength antenna spacing in both elevation and azimuth domains (d1=0.5 and d2=0.5), is used at BS. K=20 users are distributed evenly around the BS in the azimuth domain while the elevation angles obey a von Mises distribution in the LoS MIMO system. The comparative results of the system sum rate between user scheduling in Algorithm 1 and random scheduling without user selection are presented in Figure 5 under different EASs, that is, EAS=10°,30°,50°, which would provide insights on the effects of user scheduling via 3D angular separation on the system performance.Figure 5
Multiuser sum rate under different user group values. (a) Multiuser sum rate under different user group values:K=20,Φ∼U0,2π, and θ∼VM0,8.52 with EAS=10°. (b) Multiuser sum rate under different user group values: K=20,Φ∼U0,2π, and θ∼VM0,1.48 with EAS=30°. (c) Multiuser sum rate under different user group values: K=20,Φ∼U0,2π, and θ∼VM0,0.14 with EAS=50°.
(a)
(b)
(c)From the results in Figure5, an obvious feature is that a large sum rate gain is obtained by performing the user scheduling via 3D angular separation (solid red curves) as compared to the sum rate of no scheduling (solid red curve with G=1) and random scheduling (dotted dark curves). Furthermore, the more groups the multiusers are expected to be divided into, the higher the sum rate is observed, and this sum rate gain will be further enhanced. Taking the results in Figure 5(a) as an example, when the expected group number G increases from 1 to 2, 3, 4, and 5, the 0.5 CDF point sum rate gains harvested by user scheduling scheme are 3.9, 7.3, 10.3, and 13.3 bps/Hz, respectively, under EAS=10°. Once the multiusers are scheduling into different independent groups, the intergroup interference would be eliminated, which contributes to a lower total system interference and higher sum rate. Meanwhile, by performing the user scheduling algorithm with maximal angular separation among users within each group, the intragroup interference will be lowered; thus, the system sum rate is further improved.Increasing the EAS from 10° to 30° and 50° where users are better separated in the 3D space, the corresponding results of the system sum rate are illustrated in Figures5(b) and 5(c), respectively. As compared to the results in Figure 5(a), the system sum rate is improved by larger EAS as expected. A more specific result of the sum rate gain against EAS is presented in Table 1. As we can see, a large sum rate gain is obtained by user scheduling with 3D angular separation, and the gain ratio decreases as the EAS increases. For small EAS (10°), the sum rate gain can even be 3.5 times larger than that of no scheduling case, which highlights the importance of user scheduling for interference alignment in super dense area.Table 1
The sum rate gain and gain ratio (CDF=0.5) harvested by user scheduling with 3D angular separation when G=5 where Lg denotes the user number in group Λg while the moment of interference, EIkj, can be found in (B.3) in the Appendix B.
User scheduling (G=5) versus
No scheduling (G=1)
Random scheduling (G=5)
Sum rate gain (bps/Hz)/gain ratio
EAS=10°
14.0 bps/Hz/345%
4.6 bps/Hz/34%
EAS=30°
16.9 bps/Hz/114%
7.7 bps/Hz/32%
EAS=50°
13.4 bps/Hz/55%
6.1 bps/Hz/19%Although the above observations are for multiuser interference in a single cell, similar results of significant sum rate gain due to the proposed user scheduling algorithm can be expected for multiple cells in the future work. As we can conclude, in environments with considerable angular diversity, user scheduling based on angular separations will be an effective tool in interference alignment.


## 6. Conclusion

In this article, we derive the exact moments of interference in 3D LoS downlink channel where the multiusers are assumed as a uniform distribution in the azimuth domain and von Mises distribution in the elevation counterpart. The results demonstrate that the mean interference decreases as the antenna number decreases, as the antenna spacing increases and the angular separation increases. For small EAS, a large interference is observed and the spread of interference is massive; thus, carefully separated users via scheduling are essential for the satisfactory performance. Based on the 3D angular separation, we put forward a user scheduling algorithm and group users into independent subspaces with corresponding orthogonal fraction of resource blocks. Simulative results show that the system performance is improved in terms of sum rate as compared to the no user scheduling and random scheduling cases, that is, by grouping 20 users into 5 subspaces, the average sum rate increases from 4.1 bps/Hz (no scheduling) and 13.5 bps/Hz (random scheduling) to 18.1 bps/Hz (user scheduling) atEAS=10°, where the sum rate gain ratio can be as high as 334% and 34%, respectively.Collectively, the LoS interference statistics provide insights into the effects of antenna number at BS, antenna spacing, and multiuser angular separation on the system performance. The large sum rate gain and gain ratio harvested by user scheduling highlight the importance of 3D angular separation for interference alignment, especially for scenarios with multiuser in super dense area. As we can conclude, the user scheduling algorithm based on 3D angular separation enables us a promising way in interference control and management.



---

*Source: 7218471-2018-07-08.xml*