A Novel 3D Beamforming Scheme for LTE-Advanced System Yu-Shin Cheng 1, Chih-Hsuan Chen 2 Wireless Communications Lab, Chunghwa Telecom Co, Ltd No 99, Dianyan Rd, Yangmei City, Taoyuan County 32601, Taiwan E-mail: 1 cyshin@chtcomtw, 2 cchsuan@chtcomtw Abstract Beamforming is a well-known signal processing technique to increase the received signal strength to a chosen direction Recently, the three dimensional (3D) beamforming technique has gained a growing interest due to its potential to enable various strategies like user specific elevation beamforming and vertical sectorization Compared with conventional horizontal beamforming, 3D beamforming exploits the channel s degrees of freedom in the elevation direction with the active antenna system (AAS) Currently, 3GPP is on the study phase of this advanced MIMO technique and is working on the 3D channel model specification In this paper, we propose a new 3D beamforming algorithm which combines conventional horizontal beamforming and elevation beamforming Simulations are used to evaluate our proposed beamforming algorithm in urban macro environment with different inter-site distance (ISD) Index Terms LTE-Advanced, Beamforming, 3D channel I INTRODUCTION Multi-input multi-output (MIMO) is one of the most important technologies in wireless communication Long-Term Evolution (LTE) specified by 3GPP adopts this technology to improve spectrum efficiency and throughput LTE has defined several transmission modes (TM) for different MIMO transmission schemes For example, TM 3 is open-loop codebookbased precoding scheme, and TM 4 is closed-loop codebookbased precoding scheme Both of these transmission modes use PMI to choose precoding weights from codebook to apply on antenna elements Because of the limitation of codebook size, beam pattern generated by limited precoding weights cannot be adjusted precisely to optimize the performance Thanks to the non-codebook-based precoding scheme introduced in LTE- Advanced, the precoding weights could form more accurate beam pattern direct to the destination This non-codebookbased precoding technology is also known as beamforming, and the precoding weights are also called beamforming weights There are some methods to find beamforming weights [2], such as eigenstructure method which needs full channel state information to set the eigenvector of signal subspace as beamforming weights, and null-steering beamforming which needs only direction of arrival (DOA) to design beamforming weights Since beamforming weights depend on channel state information, the channel model plays an important role in the design of beamforming weights, ie, an accurate channel model aids the design of precise beamforming weights Due to the appearance of more and more high rise buildings, mobile users are sometimes located at the same horizontal angle but different elevation angles Tranditional channel models such as 3GPP spatial channel model (SCM) in [1] only describe channel in 2D (only horizontal angle) Recently, 3GPP introduces a new 3D channel model in TR 36873 [8] for LTE-Advanced system However, conventional beamforming algorithm can only be used on the 2D channel model A new 3D beamforming algorithm which considers both elevation angle and horizontal angle has been proposed The new 3D beamforming algorithm uses active antenna system (AAS) technology to adjust not only antenna weights of horizontal antenna ports but each antenna elements in the vertical direction In this way, the channel s elevation degree of freedom can be exploited The transmit power can be more concentrated on the target UEs and the interferences to other cells can be reduced In this paper, we propose a new 3D beamforming algorithm which jointly considers the antenna elements in horizontal and vertical directions The performance of our proposed algorithm is verified in different scenarios The rest of the paper is organized as follows In section 2, we describe the system model of our simulation environment Details of the 3D beamforming algorithm are shown in section 3 Section 4 shows the simulation results of the proposed algorithm in different scenarios and the comparision with conventional 2D beamforming Section 5 is the conclusion II SYSTEM MODEL There are N sites in the cellular network Each site is sectorized into three sectors, and each sector contains K users The total number of sectors in the network is M(=3N) The cellular network model is depicted in Fig 1 For the kth UE served by mth sector, the received signals can be written as : y (k) m = H (k) m w (k) m x(k) m + M K n=1,n m i=1 H (k) n w (i) n x(i) n + n (k) m (1) where x n (i) C S 1 is the transmit signals of the nth sector to the ith UE, w (i) n is the precoding vector of the nth sector to the ith UE Notice that the second term of (1) is the interference from other cells and n m (k) C S 1 is the additive white Gaussian noise In this paper, the precoding vector will be designed as beamforming weights H (k) m C T S is the channel matrix from the mth sector to the kth user, S is the number of antenna ports on BS and T is the number of antenna ports on UE The 2D antenna array model shown in Fig 2 is defined in 3GPP TR 37840 [9] Copyright IEICE - Asia-Pacific Network Operation and Management Symposium (APNOMS) 2014
ports can be written as: H (k) m = m,1,1 m,1,2 m,2,1 h(k) m,s,1 h(k) m,s,2 m,1,t h(k) m,s,t (2) where h(k) m,s,t = V i=1 ŵ s,i (3) m,s,t is the channel coefficient from the sth antenna port of the mth sector to the tth antenna port of the kth UE (we assume that each antenna port of UE has only one antenna element) ŵ s,i is the weight of the ith antenna element on the sth antenna port which affects the antenna pattern in the vertical direction Notice that in traditional passive antenna, these weights ŵ s,i are fixed, ie, they cannot be adjusted (k) h according to channel state information is the channel coefficient from the ith antenna element on the sth antenna port of the mth BS to the tth antenna port of the kth UE Fig 1 Each site is sectorized into three sectors and the arrows shows the main beam direction of each sector 1, 2 and 3 are the orders of three sectors of the cell Fig 2 Two dimensional array antenna model defined in 3GPP TR 37840 Since each antenna port contains V vertical antenna elements, the channel coefficient H (k) m composed of S antenna III THREE DIMENSIONAL BEAMFORMING ALGORITHM A Overview of Conventional Beamforming Algorithms Beamforming is a signal processing technique which applies the beamforming weights to adjust the phase and the amplitude of signals to form the beam pattern toward the desired direction The beamforming weights are applied on each antenna elements as shown in Fig 2 w s,i represents the beamforming weight of the ith element on the sth antenna port There are some methods to design beamforming weights such as eigenstruct method and null steering method [3] [4] Eigenstruct method is to consider the eigenvector with maximum eigenvalue as the beamforming weights to concentrate power to the desired direction On the other hand, null steering method is to find the steering vector which can null the interference from other sectors The concept of eigenstruct method is shown below At first we defined our channel model as (2), and we rearrange all channel coefficients in each column to form a channel vector, ie, H s =[ h s,1, h s,2,, h s,t ] T represents channel coefficient vector of sth port Therefore, we can rewrite (2) as: H = [ ] H1 H2 HS (4) To apply Eigenstruct beamforming method, we perform sigular value decomposition (SVD) on channel matrix H to find the eigenvector After performing SVD, the channel matrix H can be expressed as: H = U Σ V H (5) where U is a T T unitary matrix, Σ is a T S rectangular diagonal matrix, and V H (the conjugate transpose of V) is a S S unitary matrix T is the number of receiver antenna and S is the number of transmit antenna Then We
take the first right singular vector V 1 (the first column of matrix V) as our beamforming weights Thus in (1), w m (k) =V 1 Since the singular vectors of V are mutually orthogonal, ie, V 1 (V s ) H = 0, s 1 The transmit signal power can be concentrated on the largest singular value of matrix Σ Thus the UE received signal strength can be maximized However, conventional beamforming methods only consider two dimensional (2D) channel model Fig 3 shows the beam pattern comparison between 2D beamforming and 3D beamforming Fig 4 Compared with 3D beamforming case (right side), UE 1 will surfer from more leakage of power in fixed θ etilt case (left side) Fig 3 The comparision of beam pattern in 2D beamforming case and in 3D beamforming case We can observe in Fig 3 that the fixed vertical beam pattern in 2D beamforming will cause the decrease of SINR performance of UE 1 and UE 3 Therefore, 3D beamforming has been proposed to improve the performance of beamforming For tradtional passive antenna systems, the BS is configured with a fixed electrical downtilt angle θ etilt However, fixed θ etilt cannot satisfy all UEs with different elevation angles Furthemore, in nowadays, UEs may sometimes be served in the same building but in different floors With AAS, 3D beamforming weights can be designed to replace the beamforming weights of the fixed electrical downtilt angle, so that we can improve the performance of the scenarios as shown in Fig 4 We introduce two methods to find the beamforming weights for UEs in 3D channel model B 3D Beamforming Algorithm based on Direction of Arrival Since direction of arrival (DOA) information is much easier to estimate than full channel state information, DOA-based beamforming has been wildly discussed [5] The concept of DOA-based beamforming is that BS estimates the reference signals of the desired UE to find the DOA information of UE and utilize DOA information to design the corresponding beamforming weights The DOA based beamforming weights are the steering vectors contained DOA information and are defined as: w = 1 N 1 exp( j 2π λ (2 1)d cos( φ DOA )) exp( j 2π λ (N 1)d cos( φ DOA )) where N is the number of BS antennas, d is the BS antenna elements spacing and φ DOA is the estimation DOA of the UE The performance of received signals can be improved because steering vectors which contain DOA information will form beam pattern to the desired direction φ DOA In 3D channel model, the zenith angles between BS and UEs are different so that we need to estimate one more dimension of DOA to form beam pattern more precisely There are some methods [6] to estimate two dimensional DOA, ie, the estimated zenith angles θ DOA and the estimated azimuth φ DOA Due to the increase of DOA dimension, the beamforming weights should also be two dimensions We define horizontal beamforming weights as w H and vertical beamforming weights as w V The vertical beamforming weights are formed in the same way as (6) except that φ DOA is replaced by estimated zenith angle θ DOA and N is replaced by the number of vertical antenna elements N V And the horizontal beamforming weights are designed as : w H = 1 NH (6) 1 exp(j 2π λ (2 1)d sin( θ DOA ) sin( φ DOA )) exp(j 2π λ (N H 1)d sin( θ DOA ) sin( φ DOA )) where N H is the number of horizontal antenna elements and θ DOA is the estimated zenith degree Noticed that the horizontal beamforming weights take both azimuth and zenith DOA information into account After we obtain both w H and w V, we replace the fixed downtilt beamforming weights ŵ s,i in (3) with DOA based vertical beamforming weights w V, and DOA based horizontal beamforming weights w H is substituted for the precoding vector w (k) m in (1) With these (7)
two beamforming vectors w H and w V, we can still achieve beamforming gain in the case of 3D channel model C Proposed 3D Beamforming Algorithm based on Eigen structure method For conventional passive antenna system, the weight ŵ s,i in (3) cannot be adjusted dynamically With AAS, we now can adjust the weight according to each UE s condition In our proposed 3D beamforming algorithm, we replace fixed downtilt weight ŵ by vertical beamforming weights w v to dynamically adjust beam pattern according to every user s height and then also apply horizontatal beamforming to form beam more precisely The procedure of our proposed beamforming algorithm is as follows : Step 1: To obtain vertical beamforming weights of the sth port w V,s, we rearrange all antenna element channel coefficient into antenna port channel matrix P(k) m,s (8) and then replace H in (4) with P (k) m,s weights w V,s (P (k) m,s) H = (P (k) m,s) H w V,s = = 1,m,s,1 1,m,s,2 to find vertical beamforming 2,m,s,1 V,m,s,1 1,m,s,T V,m,s,T Step 2: We regard the vertical beamforming weights w V,s as dynamic eletrical downtilt angle of the antenna port Therefore, we combine the antenna element channel coefficient (k) of the sth antenna port into h m,s,t which represents the sth antenna port channel coefficient from the mth BS to the tth port of the kth UE, and the combination can be written as: 1,m,s,1 2,m,s,1 V,m,s,1 w V,s,1 1,m,s,2 w V,s,2 1,m,s,T w V,s,V V,m,s,T m,s,1 m,s,2 m,s,t (9) Step 3: After applying vertical beamfotrming, we have to obtain horizontal beamforming weights We rearrange the (k) channel coefficients h m,s,t which is obtained from previous steps to get the channel coefficient vector of sth port H(k) m,s = [ (k) (k) m,s,1, h m,s,2,, h m,s,t ] and then we rearrange all channel (k) vectors H m,s from all ports of BS into the complete channel matrix H (k) m (10) Then we can perform SVD on H (k) m to get the horizontal beamforming weights w H in the same way as the eigenstruct beamforming method introduced in section IIIA) (8) H (k) m = H (k) m,1 H (k) m,2 H (k) m,s T = m,1,1 m,1,2 m,1,t m,2,1 h(k) m,s,1 m,2,2 IV SIMULATION RESULTS h(k) m,s,2 m,2,t h(k) m,s,t (10) In this section, we show the performance of different beamforming algorithms in terms of post receiver signals to interference and noise ratio (SINR) The detailed system simulation parameters are defined in 3GPP TR 36873 as shown in Table 1 We assume that the sites can perfectly estimate the channel state information and DOA information of all UEs TABLE I SIMULATION CONFIGURATION PARAMETERS Parameters Values Environment 3D-UMa[8] Cellular layout 19 sites(57 sectors) Inter-site distance(isd) 500 m / 1732 m Carrier frequency 2 GHz Number of PDSCH RBs 50 Bandwidth 10MHz Channel model defined in 3GPP TR 36873 Shadowing std 7dB TX antenna port 4 TX elevation elements 4 / 10 per antenna port RX antenna port 2 RX elevation elements 1 per antenna port Antenna configuration cross-polarized Antenna spacing 05λ Maximum antenna gain 8 dbi BS TX power 46 dbm BS height 25 m UE height 3(n-1)+15, n U(1,N fl ),in meters where N fl U(4,8) UE distribution Uniform in cell UE number 10 UEs per sector UE speed 3 km/h Downtilt angle 12 (500 m) / 35 (1732 m) There are 19 hexagonal sites which is sectorized into 3 sectors with 10 UEs uniformly distributed in each sector and UEs choose which sectors to attach by geographical distance based wrapping Fig 5 shows the SINR result of ISD 500 (m) case with different number of vertical elements We compare the performance of 3 different beamforming algorithms, ie, 3D dynamic beamforming designed by Yan Li [7] (Li BF), DOA-based beamforming algorithms (DOA BF) and the proposed algorithm introduced in section II (proposed- BF) with no beamforming case (no BF) As can be seen in Fig 5, the SINR performance of our proposed beamforming is much better than DOA BF and Li BF We also show the simulation result of diffetent antenna models in Fig 6 where the vertical elements of a port are 10 instead of 4 Our propose beamforming algorithm also outperforms other beamforming
algorithms in 10 elements case The numerical result of ISD 500 (m) case is shown in Table 2 The average SINR improve more in 10 element case This is because more vertical antenna elements will make beam pattern more precise Fig 7 CDF of post receiver SINR of ISD 1732 (m) case with 4 antenna Fig 5 CDF of post receiver SINR of ISD 500 (m) case with 4 antenna Fig 8 CDF of post receiver SINR of ISD 1732 (m) case with 10 antenna TABLE III SIMULATION RESULT OF SINR IN THE SCENARIO OF ISD 1732(M) Fig 6 CDF of post receiver SINR of ISD 500 (m) case with 10 antenna TABLE II SIMULATION RESULT OF SINR IN THE SCENARIO OF ISD 500(M) no BF Li BF[7] 4 element(mean) 28317 db 38182 db 4 element (cell edge) -56637 db -45049 db 10 element(mean) 60435 db 108825 db 10 element(cell edge) -40020 db -25722 db DOA BF proposed BF 4 element(mean) 69036 db 170083 db 4 element (cell edge) -37003 db 92625 db 10 element(mean) 162070 db 266644 db 10 element(cell edge) -02175 db 218808 db Fig 7, Fig 8, and Table 3 show the performance of ISD 1732 (m) case Compared with ISD 500 (m) case, the SINR performance obtain less gain This is because when the radius become larger, the downtilt angle configuration become less no BF Li BF[7] 4 element(mean) 27673 db 34576 db 4 element (cell edge) -47611 db -44544 db 10 element(mean) 33459 db 42521 db 10 element(cell edge) -42878 db -42346 db DOA BF proposed BF 4 element(mean) 64643 db 169272 db 4 element (cell edge) -39186 db 83580 db 10 element(mean) 66245 db 199923 db 10 element(cell edge) -37566 db 72770 db apparent In spite of this, our proposed algorithm still obtain obvious gain in these two scenarios V CONCLUSION In this paper, we introduce a new 3D beamforming algorithm for 3D urban macro scenario in LTE-A network With our proposed beamforming algorithm, we form the beam pattern not only in azimuth angle but also in zenith angle so that the performance in 3D urban macro scenario can be improved Simulation results show that our proposed algorithm
obtain 20 db gain in mean SINR and 25 db gain in cell edge SINR in ISD 500 (m) case and nearly 14 db gain in mean SINR and 15 db gain in cell edge SINR in ISD 1732 (m) case We also show that increasing the number of vertical antenna elements can enhance the beamforming gain REFERENCES [1] Salo, Jari, et al MATLAB implementation of the 3GPP spatial channel model (3GPP TR 25996), on-line, 2005 Jan [2] Godara, Lal Chand Application of antenna arrays to mobile communications II Beam-forming and direction-of-arrival considerations, Proceedings of the IEEE, vol8, no8, 1195-1245, 1997 [3] Paulraj, Arogyaswami, and Thomas Kailath Eigenstructure methods for direction of arrival estimation in the presence of unknown noise fields, IEEE Transactions on Acoustics, Speech and Signal Processing, vol 34, no 1, 13-20, 1986 [4] Mouhamadou, Moctar, Patrick Vaudon, and Mohammed Rammal Smart antenna array patterns synthesis: Null steering and multi-user beamforming by phase control, Progress In Electromagnetics Research,vol 60, 95-106, 2006 [5] Krishnaveni, V, and T Kesavamurthy Beamforming for Direction-of- Arrival (DOA) Estimation-A Survey, International Journal of Computer Applications vol 61, no 11, 4-11, 2013 [6] Shahbazpanahi, Shahram, et al Robust adaptive beamforming for general-rank signal models, IEEE Transactions on Signal Processing, vol 51, no 9, 2257-2269, 2003 [7] Li, Yan, et al Dynamic Beamforming for Three-Dimensional MIMO Technique in LTE-Advanced Networks, International Journal of Antennas and Propagation, 2013 [8] 3GPP TR 36873 V200 3rd Generation Partnership Project,Technical Specification Group Radio Access Network,Study on 3D channel model for LTE (Release 12) [9] 3GPP TR 37840 V1210 Technical Specification Group Radio Access Network, Study of Radio Frequency (RF) and Electromagnetic Compatibility (EMC) requirements for Active Antenna Array System (AAS) base station (Release 12)