arxiv: v1 [cs.it] 13 Jun 2014

Transcription

1 SUBMITTED TO IEEE SIGNAL PROCESSING LETTERS, OCTOBER 8, Piot Signa Design for Massive MIMO Systems: A Received Signa-To-Noise-Ratio-Based Approach Jungho So, Donggun Kim, Yuni Lee, Student Members, IEEE, and Youngchu Sung, Senior Member, IEEE arxiv: v1 [cs.it] 13 Jun 2014 Abstract In this paper, the piot signa design for massive MIMO systems to maximize the training-based received signa-tonoise ratio SNR) is considered under two channe modes: bock Gauss-Markov and bock independent and identicay distributed ) channe modes. First, it is shown that under the bock Gauss-Markov channe mode, the optima piot design probem reduces to a semi-definite programming SDP) probem, which can be soved numericay by a standard convex optimization too. Second, under the bock channe mode, an optima soution is obtained in cosed form. Numerica resuts show that the proposed method yieds noticeaby better performance than other existing piot design methods in terms of received SNR. Index Terms Channe estimation, piot design, Gauss-Markov mode, Kaman fiter, massive MIMO I. INTRODUCTION Efficient channe estimation is a crucia probem for massive mutipe-input mutipe-output MIMO) systems [1] and there is active research going on in this area [1] [4]. Whie much research is conducted on time-division dupexing TDD) massive MIMO systems [1] [4], recenty some researchers considered the probem of efficient channe estimation and piot signa design for more chaenging frequency-division dupexing FDD) massive MIMO systems in which the number of channe parameters to estimate may be much arger than the resource aocated to training. To quicky acquire a reasonabe channe estimate with imited training resources, the authors in [5] [7] expoited the channe s spatia and tempora correation under the framework of Kaman fitering with the state-space channe mode. In particuar, the authors in [5], [6] considered the piot signa design under the state-space i.e., Gauss- Markov) channe mode to minimize the channe estimation error, and showed that the channe can be estimated efficienty by propery designing the piot signa and expoiting the channe statistics. However, minimizing the channe estimation error is not the utimate metric of data communication. Hence, in this paper, we consider the optima piot signa design under the framework of the state-space channe mode to maximize the received SNR for data transmission, which is sometimes a fina goa of data communication. Corresponding author The authors are with Dept. of Eectrica Engineering, KAIST, Daejeon , South Korea. E-mai:{jhso, dg.kim, yuniee, and ycsung}@kaist.ac.kr. This research was supported by Basic Science Research Program through the Nationa Research Foundation of Korea NRF) funded by the Ministry of Education 2013R1A1A2A ). In the mutipe-input singe-output MISO case, the training-based capacity is a monotone increasing function of the training-based received SNR [8]. A training approach based on received SNR was considered in the context of feedback in [7], [9]. The difference of this paper from [7], [9] is that we here obtained an optima piot signa under the state-space channe mode based on the training-based received SNR defined in [8], which is different from the SNR definition used in [7]. Notation: We wi make use of standard notationa conventions. Vectors and matrices are written in bodface with matrices in capitas. A vectors are coumn vectors. For a matrix A, A T, A H, A 1, TrA), ranka), λ i A), and Ai, j) indicate the transpose, conjugate transpose, inverse, trace, rank, i-th argest eigenvaue, and i, j)-th eement of A, respectivey. LA) denotes the inear subspace spanned by the coumns of A, and L A) is the orthogona compement of LA). For a random vector x, E{x} denotes the expectation of x, and x CNµ, Σ) means that x is circuary-symmetric compex Gaussian-distributed with mean µ and covariance matrix Σ. I and O denote an identity matrix and an a-zero matrix, respectivey. II. SYSTEM MODEL AND BACKGROUND In this paper, we consider the same massive MISO system as that considered in [5], [7], [10]. The transmitter has N t transmit antennas, the receiver has a singe receive antenna N t 1), and each transmit-receive antenna pair has fat fading. Under this mode the received signa y i at symbo time i is given by y i = s H i hi) +n i, i = 1,2,..., 1) where s i is the N t 1 transmit signa vector at symbo time i, h i) is the N t 1 channe vector at symbo time i, and n i is the additive Gaussian noise at symbo time i from n i CN0,σ 2 ) with the noise variance σ 2. For the channe mode, we assume the stationary bock Gauss-Markov vector process [5], [7]. That is, the channe vector is constant over one bock and changes to a different state at the next bock according to the foowing mode: h +1 = ah + 1 a 2 b, h 0 CN0,R h ), = 0,1,...2) where h is the channe vector for the -th bock, a [0,1] is the tempora fading coefficient, and b CN0,R h ) is the innovation vector at the -th bock independent of {h 0,...,h }. We assume that one bock consists of T symbos: The first T t symbos are used for training and the foowing T d = T T t symbos are used for unknown data transmission. Thus, we have h i) = h for i = T + m, m = 1,2,,T. It is easy to verify the assumed timewise stationarity, i.e., R h = E{h 0 h H 0 } = E{h 1 h H 1 } =, for the considered channe parameter setup. R h captures the spatia correation of the channe and depends on the antenna geometry and the scattering environment [13]. We assume that a and R h are known to the system. Pease see [5] We assume that stationarity hods at east ocay [11], [12]. That is, the channe statistics vary much sowy than channe s fast fading.

2 2 SUBMITTED TO IEEE SIGNAL PROCESSING LETTERS, OCTOBER 8, 2018 regarding this assumption.) Let R h = UΛU H be the eigendecomposition of R h, where U is a N t R c matrix composed of orthonorma coumns and the R c R c matrix Λ contains a the non-zero eigenvaues of R h. Since a {h, = 0,1,...} are contained in the same subspace LU), we can mode the -th bock channe as h = Ug because of the assumed stationarity. Then, the channe dynamic 2) can be rewritten in terms of g as g +1 = ag + 1 a 2 e, g 0 CN0,Λ), = 0,1, 3) with e CN0,Λ). This random vector process is again a stationary process with E{g g H } = Λ for a ). By stacking the symbo-wise received signa in 1) corresponding to the training period of each bock, we have y = S H h +n, 4) where y = [y T+1,y T+2,,y T+Tt ] T, S = [s T+1,,s T+Tt ], and n = [n T+1,n T+2,,n T+Tt ] T. The tota power aocated to the training period of each bock is given by TrS S H ) ρt t, which means that each piot symbo has power ρ on average. Since h LU), there is no oss in setting s i = U s i because the signa power aocated to L U) wi simpy be ost without affecting the received signa y i. Hence, we have S = [U s T+1,,U s T+Tt ] = U S, 5) where S is a R c T t matrix and we assume R c T t, i.e., the number of symbos contained in one channe coherence time is smaer than the channe rank as in typica massive MIMO systems. Then, the measurement mode 4) is rewritten as y = U S ) H Ug )+n = S H g +n, 6) and the power constraint on S is given by Tr S SH ) = TrS S H ) ρt t. Thus, the origina state-space mode 2) and 4) is equivaent to the new mode 3) and 6) under the known stationary subspace condition h = Ug. Under the statespace mode 3) and 6), the optima minimum mean-squareerror MMSE) channe estimation is given by Kaman fitering [14]. That is, the MMSE estimate ĝ and its estimation error covariance matrix P are updated as foows [14]: K = P 1 S σ 2 I+ S H P 1 S, ĝ = ĝ 1 +K y S H ĝ 1 ), P = I K SH )P 1, ĝ 1 = aĝ 1 1, P 1 = a 2 P a 2 )Λ, 7) where ĝ := E{g y 0,y 1,,y }, P := E{g ĝ )g ĝ ) H }, ĝ 0 1 = 0, and P 0 1 = Λ. III. PROBLEM FORMULATION In this section, we consider the piot design probem to maximize the received SNR for the data transmission period under the assumption that T and T t are given and the transmit beamforming is used for the considered MISO channe during the data transmission period, i.e., s i = w i d i = U w i d i, i = T+T t +m, m = 1,,T d, 8) wherew i and d i are the transmit beamforming vector and data symbo for symbo time i. Here, we assume E{d i } = 0 and E{ d i 2 } = σ 2 d. From here on, we set σ2 = 1 for simpicity. Again due to h LU), we can set w i = U w i without any performance oss. From now on, we use i) instead of i for i = T+m, m = 1,,T. First, foowing the framework in [8], we derive the received SNR during the data transmission period. The true channe at symbo time i) is expressed as h i) = ĥi) i) + h i), 9) where i) is the bock number corresponding to symbo time i), ĥi) i) := Uĝ i) i) with ĝ i) i) obtained from 7) is the MMSE estimate for h i) = Ug i) ) this is true because TrE{g ĝ )g ĝ ) H }) = TrE{g ĝ )g ĝ ) H }U H U) = TrE{h ĥ )h ĥ ) H })), and h i) is the channe estimation error. Substituting 8) and 9) into 1), we have y i) = d i) w H i) ĥi) i) + h i) )+n i), = d i) w H i ĝi) i) +d i) w H i g i) +n i) ). 10) The key point in [8] is that in the right-hand side RHS) of 10), the term w i)ĝi) i) H is known to the receiver and the terms w i) H g i) and n i) are unknown. Hence, the trainingbased received SNR is defined as [5], [8] ) w i) H ĝ i) i) ĝi) i) H w i) SNR i) = Pi) i) +γ 1 I ), 11) w i) w H i) where γ := σd 2/σ2, since P i) i) = E{ g i) gi) H }. The optima beamforming vector that maximizes SNR i) is given by soving a generaized eigenvaue probem. In genera, a cosed-form soution to a generaized eigenvaue probem is not avaiabe. However, since the rank of ĝ i) i) ĝi) i) H in the numerator of the RHS of 11) is one, one can easiy sove the probem in this case, and the optima beamforming vector w i) and the corresponding optima SNR i) are given by w i) = P i) i) +γ 1 I ĝi) i), 12) SNR i) = ĝh i) i) Pi) i) +γ 1 I ĝi) i). 13) Note that the optima received SNR is the same for a data symbosi = T+T t +m,m = 1,,T d of each bock. Hence, we sha use the notation SNR for SNR i). Aso, note from 13) that the optima SNR is a function of symbo SNR γ, the error covariance matrix P i) i) and the channe estimate ĝ i) i). Hence, simpy minimizing the trace of P i) i) may not be optima to maximize the received SNR due to the term ĝ i) i). Using the fact that both P i) i) and ĝ i) i) are functions of the piot signa S, as seen in 7), we can express the optima SNR as a function of S, given by SNR = ĝ 1 +K y S H ĝ 1 )) H I K SH )P 1 +γ 1 I ĝ 1 +K y S H ĝ 1 )). 14) Our goa is to design the sequence{ S, = 0,1,2, } of piot matrices to maximize SNR. However, SNR is a function of

3 SO, KIM, LEE AND SUNG 3 a previous piot signa matrices via P 1 and ĝ 1, and the design probem is a compicated joint probem. Thus, as in [5], [10], we adopt the greedy sequentia approach and the design probem is expicity formuated as foows. Probem 1: Given the channe statistics information, a and R h, and a previous piot matrices{ S 0, S 1,, S 1 }, design S such that max S E{SNR } subject to Tr S SH ) ρt t. 15) Here, the expectation in 15) is to average out the randomness in the random vector y. IV. THE PROPOSED DESIGN METHOD To sove Probem 1, we begin with the foowing proposition. Proposition 1: The piot design probem 15) is equivaent to the foowing optimization probem: min Tr A B + S ) SH S 16) subject to Tr S SH ) ρt t, where A = γĝ 1 ĝ 1 H +γp 1+I andb = γi+p 1 1. Note that A and B are not functions of the design variabe S. Proof: From 13) the average received SNR,E{SNR }, with the optima beamforming vector wi) can be expressed as [ P E{SNR } = Tr +γ 1 I ) ] 1 E{ĝ ĝ H }. 17) Since ĝ is a Gaussian random vector with mean ĝ 1 and covariance matrix Q given by Q =P 1 S I+ S H P 1 S SH P 1 =P 1 S SH P S SH, 18) where the second equaity hods by the matrix inversion emma, E{ĝ ĝ H } is given by E{ĝ ĝ H } = ĝ 1ĝ H 1 +P 1 S SH P S SH. 19) The error covariance matrix P is expressed as P = P 1 K S P 1 = P 1 P 1 S SH P S i SH i. 20) Substituting 19) and 20) to 17), we have [ Tr P 1 P 1 S SH P S SH +γ 1 I ĝ 1 ĝ 1 H +P 1 S SH P S SH )] [ =Tr P 1 +γ 1 I ) P S SH ) P 1 S SH ) ] ĝ 1 ĝ 1 H P 1 1 +ĝ 1ĝ 1 H +P 1) S SH [ =Tr γ γi+p S { SH ĝ 1 ĝ 1 H P 1 1 +ĝ 1 ĝ H 1 +P 1)γI+P S SH ) } ] ĝ 1 ĝ 1 H +P 1)γI+P 1 1 ) [ ] [ =γtr ĝ 1 ĝ 1 H +P 1 1 γtr γĝ 1 ĝ 1 H +γp 1 +I)γI+P S SH ]. 21) Here, we used TrAB) = TrBA) and AB = B 1 A 1. Since the first term of the RHS of 21) is independent of S and the second term of the RHS of 21) is TrA B + S SH ) with A and B defined in the proposition, the probem 15) is equivaent to the probem 16). Note that the probem 16) is not a convex optimization probem. To tacke the probem 16), we use the semi-definite reaxation SDR) technique [15]. First, introducing a new variabe X := S SH, we change the optimization probem 16) as min Tr A B +X ) X subject to TrX ) ρt t, 22) X 0, rankx ) T t. Then, dropping the rank constraint in the probem 22), we change the probem to the foowing optimization probem: min Tr A B +X ) X subject to TrX ) ρt t, 23) X 0. Since A and B are positive-definite matrices, the probem 23) is a convex optimization probem and can be soved by a standard convex optimization sover. To obtain a soution matrix S of size R c T t from the soution X of 23), we use a randomization technique. That is, we generate T t random vectors according to the distribution CN0,X ). After the generation of these random vectors, we stack the vectors to make a R c T t matrix S. Since A and B can be obtained by the standard Kaman recursion, ony soving the probem 23) and appying the randomization technique are additionay necessary to design the received-snr-optimized piot sequence. A. The Bock I.I.D. Channe Case The bock channe case [13] is a specia case of the mode 2) or 3) with a = 0. Under this mode, the Kaman recursion 7) is sti vaid athough the recursion does not propagate, i.e., ĝ 1 = 0 and P 1 = Λ for every. Hence, Proposition 1 is vaid under the bock channe mode. In this case, P 1 = Λ is a diagona matrix and thus, the matrices A and B in Proposition 1 are diagona. In this case, the optimization probem 16) can be soved efficienty without soving 22) based on the foowing proposition. Proposition 2: There exists an optima soution to the probem 16) in the form of S = ΠD, where Π is a R c R c permutation matrix and D is a R c T t diagona matrix in the form of D = δ O 0 0 δ Tt T, δ i 0 i, 24)

4 4 SUBMITTED TO IEEE SIGNAL PROCESSING LETTERS, OCTOBER 8, 2018 when A and B are diagona matrices. Proof: The proof is simiar to that of [13, Theorem 3]. Since A is a positive definite matrix, the objective function of the probem 16) can be rewritten as Tr A 1 2 B A H 2 +A 1 2 S SH A H 2 ), 25) where A = A 1/2 A H/2. Let C := A 1 2 B A H 2 + S SH A H 2, λc ) := [λ 1 C ),,λ Rc C )] T and A 1 2 dc ) := [C 1,1),,C R c,r c )] T. Then, the objective function 25) can be rewritten as fλc )) := N 1 i=1 λ ic ), since the trace of a matrix is the sum of its eigenvaues. It is shown in [13, Theorem 3] that fλc )) is ower bounded by fdc )), i.e. fλc )) fdc )), based on the Schur convexity of f ). This ower bound can be achieved when C is a diagona matrix. To make C a diagona matrix, X = S SH shoud be a diagona matrix, since A and B are diagona matrices. Therefore, the minimum vaue of the objective function can be achieved when S SH is a diagona matrix. By decomposing the R c R c diagona matrix X = S SH of rank ess than or equa tot t, we have a soution to 16) in the form of S = ΠD. The ocations of the nonzero eements of X determine Π.) Using Proposition 2, the Lagrange mutipier technique and the fact that A = γb γi +I, we obtain the optima diagona eements {x i } of X = S SH given by ) B i,i) x i = max B i,i)+ νb i,i) γ), 0 26) ) 1 = max γ λ i R h ) + γλi R h )+1, 0. 27) ν Since the object function in 16) can be rewritten as Rc A i,i) i=1 B i,i)+x i and the term A i,i) B i,i)+x i is a monotone increasing function of B i,i), the indices with the smaest T t B i,i) vaues shoud be seected for possiby non-zero T t x i s. Let this index set be denoted by I. Then, the Lagrange mutipier ν is obtained to satisfy the power constraint i I x i = ρt t by the bisection method. The proposed index seection here corresponds to seecting the T t dominant eigendirections of R h since B = γi + P 1 1 = γi + Λ 1. Interestingy, this index seection method coincides with the resut in [13] minimizing the channe estimation MSE. The channe estimation MSE minimizing probem is equivaent to 16) with redefined A := I and B := Λ 1.) In both received SNR maximization and channe estimation MSE minimization, the T t dominant channe eigen-directions shoud be used for piot patterns, but the power aocation is a bit different. Remark 1: By Proposition 2, in MISO systems with the bock channe mode, a received-snr-optima piot signa is given by S = UΠD. Hence, there is no need to mix mutipe channe eigen-directions at a symbo time to improve the performance. At each symbo time, it is sufficient to use one coumn of U. On the other hand, in the bock-correated channe case a 0), the optima soution X to 22) is not diagona in genera and thus, mixing mutipe channe eigendirections at a symbo time can improve the received SNR performance. V. NUMERICAL RESULT In this section, we provide some numerica resuts to evauate our piot design method. We set 2 GHz carrier frequency, T s = 100µs symbo duration, bock size T = 10 with three training symbos per bock T t = 3), and the pedestrian mobie speed v = 3km/h a = ). The tempora fading coefficient a is given by a = J 0 2πf d T s T) by Jakes mode [16], where f d is the maximum dopper frequency and J 0 is the 0-th order Besse function.) For the channe spatia correation matrix R h, we consider the exponentia correation mode given by R h i,j) = r 2 i j with r = 0.9. NMSE Received SNR db) Orthogona Random Noh et a[5] Noh et a[6] Proposed PerfectCSI BockIndex ) BockIndex ) Fig. 1. NMSE and received SNR versus bock index : N = 16, T t = 3, ρ/σ 2 = 10dB, γ = 10dB, and v = 3km/h Fig. 1 shows the performance of the proposed piot design, when ρ/σ 2 = γ = 10dB and N t = 16. The normaized MSE NMSE) is defined as g ĝ 2 g. The resut is averaged 2 over 100 random reaizations of the channe process with ength 40 bocks. For comparison, we consider orthogona and random beam patterns for N t = 16. In addition, we consider the piot design agorithms minimizing the channe estimation MSE in [5], [6]. It is seen that the proposed method noticeaby outperforms other methods in terms of received SNR and especiay yieds quick convergence at the eary stage of channe earning, athough its MSE performance is worse than the methods in [5], [6]. Athough the resut is not shown here due to space imitation, it is observed in the bock channe case that the proposed piot design method in Section IV-A yieds sighty better performance than the method in [13] in terms of received SNR. VI. CONCLUSION In this paper, we have considered the piot signa design for massive MIMO systems to maximize the received SNR under the bock Gauss-Markov and bock channe modes. We have shown that the proposed design method yieds noticeaby better performance in terms of received SNR than channe estimation MSE-based methods. Furthermore, we have shown that using the T t dominant eigen-vectors of the channe covariance matrix without mixing as the piot signa provides an optima soution even for received SNR maximization under the bock channe mode. The extension to the MIMO case is eft as future work.

5 SO, KIM, LEE AND SUNG 5 REFERENCES [1] F. Rusek, D. Persson, B. K. Lau, E. G. Larsson, O. Edfors, F. Tufvesson, and T. L. Marzetta, Scaing up MIMO: Opportunities and chaenges with very arge arrays, IEEE Signa Process. Mag., vo. 30, pp , Jan., [2] T. L. Marzetta, Noncooperative ceuar wireess with unimited numbers of base station antennas, IEEE Trans. Wireess Commun., vo. 9, pp , Nov [3] J. Hoydis, S. ten Brink, and M. Debbah, Massive MIMO in the UL/DL of ceuar networks: How many antennas do we need?, IEEE J. Se. Areas Commun., vo. 31, pp , Feb [4] C. Shepard, H. Yu, N. Anand, L. E. Li, T. L. Marzetta, R. Yang, and L. Zhong, Argos: Practica many-antenna base stations, in Proc. MobiCom, Istanbu, Turkey), Aug [5] S. Noh, M. D. Zotowski, Y. Sung, and D. J. Love, Piot beam pattern design for channe estimation in massive MIMO systems, to appear in IEEE J. Se. Topics Signa Process., Oct avaiabe at IEEExpore). [6] S. Noh, M. D. Zotowski, Y. Sung, and D. J. Love, Tranining signa design for channe estimation in massive MIMO systems, in Proc. ICASSP to appear), Forence, Itay), May [7] J. Choi, D. J. Love, and P. Bidigare, Downink training techniques for FDD massive MIMO systems: Open-oop and cosed-oop training with memory, to appear in IEEE J. Se. Topics Signa Process., Oct avaiabe at IEEExpore). [8] B. Hassibi and B. M. Hochwad, How much training is needed in mutipe-antenna wireess inks?, IEEE Trans. Inf. Theory, vo. 49, pp , Apr [9] D. J. Love, J. Choi, and P. Bidigare, A cosed-oop training approach for massive MIMO beamforming systems, in Proc. IEEE CISS, Johns Hopkins Univ., Maryand), Mar [10] S. Noh, M. D. Zotowski, Y. Sung, and D. J. Love, Optima piot beam pattern design for massive MIMO systems, in Proc. Asiomar, Pacific Grove, CA), Nov [11] S. Stein, Fading channe issues in system engineering, IEEE J. Se. Areas Commun., vo. 5, pp , Feb [12] A. F. Moish, Wireess Communications. New York, NY: Wiey, [13] J. H. Kotecha, and A. M. Sayeed, Transmit signa design for optima estimation of correated MIMO channes, IEEE Trans. Signa Process.,, vo. 52, pp , Feb [14] T. Kaiath, A. H. Sayed, and B. Hassibi, Linear Estimation. Upper Sadde River, New Jersey: Prentice-Ha, [15] Z. Luo, W. Ma, A. M. So, Y. Ye, and S. Zhang, Semidefinite reaxation of quadratic optimization probems, IEEE Signa Process. Mag., vo. 27, pp , May [16] W. C. Jakes, Microwave Mobie Communication. New York, NY: Wiey, 1974.