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

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.


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 where K 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 signalto-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 spacedivision 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 with EAS < 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 Section 2, 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.

System Model and Performance Metrics
Given the typical URA at BS, as shown in Figure 1, there are M vertical elements spaced by d 1 wavelength and N horizontal elements with d 2 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 where β k is the link gain, and λ is the wavelength.We next rewrite the M × N matrix H k as a vector h k in a dimension of 1 × MN as follows: Thus, the received downlink signal at UE k is obtained where x is the transmitted signal and n k ∼ CN 0, σ 2 z is the additive white Gaussian noise with mean value 0 and variance σ 2 z .Assuming linear precoding at the BS, where , and the SINR at UE k, assuming E d κ 2 = 1, is given by In ( 7), we observe that the intracell interference to UE k from the transmitted signal to UE j is β k I kj / MN , where Thus, the total intracell interference to UE k is β k I k / MN , where The desired signal power is MNβ k ; the SINR and the rate of UE k can be expressed as We note that the moments can be used to estimate the mean rate through the delta method [28], where the firstand second-order center moments of I k are included: where 2  13 Clearly, the performance is critically dependent on the distribution of the interference I k .Hence, we will derive the moments of I k , that is, E I k and Var I k .
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 Here, κ is a concentration parameter, μ is a location parameter, and I 0 • is the modified Bessel function of the first kind in the zeroth order.

Interference Moments
From ( 11)-( 13), we observe that E I k and Var I k are required to compute the sum rate.Hence, we compute the first two moments of I i based on the following lemma.14), the approximation with the approximation accuracy validated in Figure 2, enables us to obtain

16
where 3 International Journal of Antennas and Propagation πd n,q y J r/2 πd n,q y , 18 when L = 2L 0 is even in (15), The parameters mentioned in ( 16)-( 22) are given below: (i) d n,q y and d m,p z are defined as

23
(ii) a and b are defined as , and l 4 are defined as (1.320.1)]: when L = 2L 0 is even, Since the double integral E e −j m−p u k e −j n−q v k occurs throughout the moment derivations, we use the short notation, F m − p, n − q, κ , defined in (16).From Lemma 1, the following results can be obtained. where Theorem 2.
From Theorems 1 and 2, the variance of I k follows from Note that, when κ = 0, we have the special case where elevation angle is uniformly distributed with θ ∼ U 0, π .

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π].Specifically, results are shown for the rectangular antenna array with 100 (10 × 10) and 36 (6 × 6) elements at BS with antenna spacing d = 0 25λ, d = 0 5λ, and d = 1λ (d 1 = d 2 = 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 −j m−p u k e −j n−q v k 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 E I k 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 MN 2 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.

Numerical and Simulation
Results: Variance.In Figure 4, we plot the spread of the interference by including levels which are one standard derivation above and below the mean, that is, E I k ± std I k , where std I k = Var I k .The same parameters such as K = 5 and a rectangular antenna array with 100 (10 × 10) elements at BS with d 1 = 0 5λ and d 2 = 0 5λ are used in the simulation.
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 5 International Journal of Antennas and Propagation 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).

User Scheduling with 3D Angular Separation
As indicated by the results in Figures 3 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.

User Scheduling Algorithm with 3D Angular Separation.
Assume there are K 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 (d i , θ i , φ i ), where d i , θ 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 ϕ ∼ U 0, 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, d 0 , and the cell radius, R. The path loss of the signal from the BS to user i is given by [30, eq (2.28)] where P ti is the signal transmission power to user i, P a 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 x i = P i sin θ i cos ϕ i , P i sin θ i sin ϕ i , P i cos θ i , x j = P j sin θ j cos ϕ j , P j sin θ j sin ϕ j , P j cos θ j 34 Then, the angular separation between users i and j, Δ ij , is calculated as Thus, we could obtain the angular separation matrix Δ as follows: which will be used as the input of the user scheduling algorithm below: The following are the parameters: Λ k denotes the users selected into the kth group, I G denotes the users which were divided while I L denotes the users left to be grouped, S i denotes the sum of the angular separation between user i and all users in group G, and i 1 and i 2 denote the users to be grouped in the current step with the target group number g 1 and g 2 , respectively.
In Algorithm 1, 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 × K K − 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 G G − 1 /2 cases to allocate users i 1 and i 2 to two different groups.For each case, the averaged angular separation of multiuser is calculated.By maximizing the averaged angular separation (step 17), users i 1 and i 2 are distributed to groups Λ g 1 and Λ g 2 , respectively; thus, the calculation complexity is K − G G G − 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 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.

System Sum Rate Calculation.
Based on Algorithm 1, 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 where I g k denotes the interference to user k from other users in group Λ g .

6
International Journal of Antennas and Propagation We assume that the angles of all users are independent; the first moment E I k in ( 9) can be simplified as Thus, we would obtain 5.3.Simulative Results and Analysis.In the simulation, we assume the same antenna array, the URA in Figure 1 with 100 elements (10 × 10) and half wavelength antenna spacing in both elevation and azimuth domains (d 1 = 0 5 and d 2 = 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.
From the results in Figure 5, 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 Figures 5(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 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 (minimum angular separation) 5: (the set of grouped users) 7: Δ i 1 ,i 2 = inf (update Δ) 8: for g = 3 G do (search the initial user for group 3 to G) 9: (users left to be grouped) 10: (users to be grouped) 17: break (ensure all K users are grouped) 23: end if 24: end for 25: end if Algorithm 1: User scheduling with 3D angular separation.7 International Journal of Antennas and Propagation 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.
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.

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) at E AS = 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.

Appendix A. Proof of Lemma 1
Given the uniform distribution in azimuth domain with ϕ ∼ U 0, 2π , we have For simplicity, we ignore the subscript k in the derivation.Then, substituting f θ in (14) with θ ∼ VM κ, π/2 and the approximation (15) into (A.1),we would obtain where a and b are defined in (24).Further, expanding and substituting (A.3) into (A.2),we would obtain (A.4).Then, the integral calculation of (A.4) based on lemmas in [29, eq. ( 6.6818, 6.6819)] would lead to (A.5), which can be further simplified as (A.6) with the desired parameters A and B in ( 17) and ( 18), respectively.Thus, the remaining work in the proof of Lemma 1 will be the derivation of the exact expressions of C and D in ( 19)-( 22) based on the preliminary definitions below: 9 International Journal of Antennas and Propagation C = π 0 sin θ L J 0 2πd n,q y sin θ dθ, A 7 D = π 0 sin θ L J 0 2πd n,q y sin θ cos rθ dθ A 8 Case 1.When L = 2L 0 is even, substituting (26) into (A.7) and (A.8) will lead to the result in (A.9) and (A.10), and then performing the integral calculation based on [29, eq. ( 6.6818, 6.6819)] will lead to the desired result in (19) and (20).Case 2. When L=2 L 1 − 1 is odd, substituting (28) into (A.7) and (A.8) will lead to the result in (A.11) and (A.12), respectively, and then performing the integral calculation based on [29, eq. ( 6.6818, 6.6819)] will lead to the desired result in (21) and (22).

B. Proof of Theorem 1
As we can see from (8), we need to obtain the term a k a + j based on the definition in (4).e −j m−p u k −u j e −j n−q v k −v j = E 〠 m,n,p,q e −j m−p u k e −j n−q v k e j m−p u j e j n−q v j , B 2 where Σ m,n,p,q denotes Σ M−1 m=0 Σ N−1 m=0 Σ M−1 p=0 Σ N−1 q=0 As users k and j are assumed to be independent, then (B.2) can be further simplified as E I kj = 〠 m,n,p,q E e −j m−p u k e −j n−q v k E e j m−p u j e j n−q v j = 〠 m,n,p,q E e −j m−p u k e −j n−q v k 2 B 3 Substituting E I kj of (B.3) into (38), we can obtain Using Lemma 1 for the expectation expressions gives the desired result of (30).

C. Proof of Theorem 2
Assuming the angle distributions of different users are independent, then from (9), we have Define d n 1 n 2 ,q 1 q 2 y and d m 1 m 2 ,p 1 p 2 y as Similar as (B.3), we can have the following result after several simplifications u k e −j2πd n 1 n 2 ,q 1 q 2 y v k E e j2πd m 1 , p1 z u j e j2πd n 1, q 1 y v j E e j2πd m 2 ,p 2 z u i e j2πd n 2 ,q 2 y v i

C 3
Using Lemma 1 for the expectation expressions gives the desired result of (31).

Figure 1 :
Figure 1: Multiuser LoS interference with URA at BS.

ae 1
k a + j = H h,k, H h,k e −ju k , … , H h,k e −j M−1 u k ⋅ H * h,j , H * h , j e ju j , … , H * h,k e j M−1 u j T = H h,k H * h,j 1 + e −j u k −u j + ⋯ + e −j M−1 u k −u j = 〠 N−1 n=0 e −jn v k −v j 〠 M−1 m=0 e −jm u k −u j = 〠 −jm u k −u j e −jn v k −v j BThus, substituting (B.1) into (8) will giveE I kj = E 〠 u k −u j e −jn v k −v j 〠 k −u j e jq v k −v j

E
I kj I ki = K − 1 2 E I kj I ki C 1

2
[27]rnational Journal of Antennas and Propagationk is denoted by (φ k , θ k ).The channel impulse response (CIR) between UE k and the URA, H k , is given by[27]

Table 1 :
The sum rate gain and gain ratio (CDF = 0 5) harvested by user scheduling with 3D angular separation when G = 5 where L g denotes the user number in group Λ g while the moment of interference, E I kj , can be found in (B.3) in the Appendix B.