In this paper, a group-blind multiuser receiver is investigated for multicarrier code division multiple access (MC-CDMA) systems. In uplink transmissions, the receiver has the knowledge of the spreading sequences of the users within the cell, while the spreading sequences of the users from other cells are unknown. By imposing linear vector constraints on the multiuser detector, group-blind detection can efficiently suppress the known interferers, while mitigating the unknown interferers based on the minimum mean square error (MMSE) principle. A subspace-based channel estimation algorithm is proposed to blindly estimate the channel information. An adaptive implementation of the group-blind multiuser detection is also developed by using signal and noise subspace tracking algorithms. Numerical results show that the performance of the group-blind detector is much better than the subspace-based blind linear MMSE detector.published_or_final_versio

Chan, SC

Cheng, H

HKU Scholars Hub

Title A group-blind multiuser receiver for MC-CDMA systemsAuthor(s) Cheng, H; Chan, SCCitation Proceedings - Ieee International Symposium On Circuits AndSystems, 2004, v. 3, p. III597-III600Issued Date 2004URL http://hdl.handle.net/10722/46432Rights Creative Commons: Attribution 3.0 Hong Kong LicenseA GROUP-BLIND MULTIUSER RECEIVER FOR MC-CDMA SYSTEMS H. Cheng and S. C. Chan Dept. of Electrical & Electronic Engineering The University of Hong Kong, Pokfulam Road, Hong Kong ABSTRACTIn this paper, a group-blind multiuser receiver is investigated for multicarrier code division multiple access (MC-CDMA) systems. In uplink transmissions, the receiver has the knowledge of the spreading sequences of the users within the cell, while the spreading sequences of the users from other cells are unknown. By imposing linear vector constraints on the multiuser detector, group-blind detection can efficiently suppress the known interferers, while mitigating the unknown interferers based on the minimum mean square error (MMSE) principle. A subspace-based channel estimation algorithm is proposed to blindly estimate the channel information. An adaptive implementation of the group-blind multiuser detection is also developed by using signal and noise subspace tracking algorithms. Numerical results show that the performance of the group-blind detector is much better than the subspace-based blind linear MMSE detector. 1. INTRODUCTION MC-CDMA was recently proposed as an efficient multicarrier transmission scheme for supporting multiple access commu-nications [1]. It combines CDMA and orthogonal frequency division multiplexing (OFDM) techniques. It has received considerable attention because of its advantages in frequency diversity, multipath fading resilience and etc [2]. Similar to conventional CDMA systems, MC-CDMA is an interference-limited system. It is rather sensitive to multiple access interference (MAI), which arises from the other users signals transmitted using the same frequency band. Recently, a subspace-based blind MMSE detection has been proposed in [3] for MAI suppression in CDMA systems. By estimating signal subspace from the received signals, the linear MMSE detector can be blindly obtained with the known spreading sequence of the desired user. This method is especially attractive for the CDMA downlink transmissions, where the receiver only knows its own spreading sequence. However, in uplinks, the base station receiver has the knowledge of the users within the cell, but not that of the users from other cells. In [4], subspace-based group-blind multiuser detection is proposed for CDMA uplinks, which zero-forces the interferers within the cell, while suppresses the unknown interferers from other cells with the MMSE principle. By imposing linear vector constraints, the group-blind multiuser detection considerably outperforms the subspace-based blind method in CDMA uplinks.  In [5] [6], subspace-based blind detection and channel estimation techniques have been proposed for MC-CDMA downlinks. In this paper, the group-blind receiver in MC-CDMA uplinks is studied. To investigate the performance of the group-blind and blind detection in practical applications, the adaptive implementation of the two detectors is developed based on the signal and noise subspace trackers. A large number of subspace trackers with various complexities and performance have been proposed [7]. By computational complexity, subspace tracking algorithms can be classified in )( 2rMO , )( 2MrO  and )(MrOtechniques, where M  is the input vector dimension and r  is the number of eigenvectors to be tracked. In particular, the bi-iteration square-root SVD subspace tracking algorithm (Bi-SVD) [8] with a complexity of  )( 2sMrO  ( sr is the rank of signal subspace), and the projection approximation subspace tracking with deflation (PASTd) algorithm [9] with a complexity of )( sMrO  are used for tracking principle eigenvectors and eigenvalues; while the square-root QR inverse iteration method [10] with a complexity )( 2nMrO  ( nr is the rank of noise subspace, and Mrr ns   ) is adopted for minor subspace tracking. In the numerical simulations, we found that the convergence of the PASTd algorithm is quite slow when the rank of signal subspace is large and input signal-to-noise ratio (SNR) is high. Whereas, the Bi-SVD signal subspace tracker presents performance close to the solution with exact SVD.  The rest of the paper is organized as follows: Section 2 describes the principle of MC-CDMA systems. Section 3 is devoted to the group-blind multiuser receiver, while the Bi-SVD signal subspace tracking algorithm is presented in Section 4. Section 5 gives the numerical examples, and finally, conclusion is drawn in Section 6.2. SYSTEM MODEL We consider a MC-CDMA system for uplink transmission from the terminal users to the base station. Fig. 1 shows the structure of an MC-CDMA transmitter. The original data stream of the kth  user is first converted into P  parallel data sequences. Each serial/ parallel (S/P) converted output spreads with a spreading sequence of length M . The data chips after spreading are S/P converted into M  parallel subcarriers. The total PMN   chips are modulated by the IDFT and converted back into serial data. The cyclic prefix is inserted between symbols to avoid inter-symbol interference caused by multipath fading. Let T  be the symbol duration of the original data stream, and sT  be the one at the subcarrier, then it has TPTs  . The frequency separation between the successive subcarriers is sTf /1 ' . The complex,,,;,((( ,6&$6equivalent transmitted signal of user k  can be written as tfpPmjPpMmsimkipkiikk eiTtcbPts'  f ¦ ¦ ¦ )(21010)(,)(, )()(S\ , (1)where kP  is the chip energy, )(,ipkb  is the transmitted thi  data bit in the thp  data stream, )(,imkc  is the thm  chip of the given spreading sequence, and )(t\  represents the rectangular pulse waveform. Note that there are P  data streams for the kth  user, and the data bits on different streams use the same spreading sequence. In order to achieve the maximum frequency diversity, the subcarrier separation for )(,npkb  is sTPf /'  ' . The kthuser’s signal )(tsk  propagates through a frequency-selective Rayleigh fading channel with L  received paths. The channel is represented by a time-variant impulse response  ¦   Lllklkk tttg1,, )()()( WGD , (2)where k  is the user index, )(, tlkD is the lth  path gain which is independent zero mean, complex Gaussian random process, and lk ,W  is the propagation delay for the lth  path, respectively.  In MC-CDMA systems, because of the large symbol duration, each subcarrier approximately experiences flat fading. Assume that there are K  simultaneous users in the system, all using the same selection of P  and M . At the receiver, the signal is sampled at a rate sTN / , the samples corresponding to the cyclic prefix are discarded, and a DFT of size N  is performed. For simplicity, we will consider a synchronous model, where all the users arrive at the same time. The received signal TpMPpp iyiyi )](,),([)( )1(  y after DFT operation with respect to the data symbol )(,ipkb  can be given as )()()(10)(,)(,10)(,)(,ibPibPiKk kipkipkkKkipkipkkkpndnhCy ¦  ¦    , (3)where (for convenience, the subscript i  is omitted) pkb ,  is the data symbol corresponding to the pth  data stream; kP  is the chip energy; pk ,h  denotes channel frequency response, whose element is 1,,0),(1)(210,,  ¦   MmeNh NpPmljLllkpPmk SD ;]}1[,],1[],0[{  Mcccdiag kkkk C  is the code matrix; pk ,dTpMPkMkpkkhchc ],,[ )1(,1,0   denotes the effective signa-ture; and TpMPpPp ],,,[ )1(  KKK n  represents complex additive white Gaussian with zero mean and variance 2vV .3. GROUP-BLIND MULTIUSER DETECTION Assume that the data stream for each user is independent, i.e., 0)}()({ *  tbtbE ji , and the ambient noise is white. Hence, referring to the linear model in Eq. (3), the data correlation  matrix yR  of the received signal can be given by MMnKkkH PnnE u ¦  ǿddyyR Hkky210)}()({ V , (4)Following the eigenvalue decomposition, yR  is given as ,[HnnnHsssHnHsnsnsHUȁUUȁUUUȁ00ȁ]UUUUȁR y »»¼º««¬ª»¼º«¬ª  (5)where ),,( 21 Kdiag OOO  sȁ  contains K  largest eigenvalues of yR , which corresponds to the K  active users,  KMuCUscontains K  signal eigenvectors which span the signal subspace, and )( KMMnuCU  contains KM   noise eigenvectors which span the noise subspace.  Consider an uplink transmission, where the base station receiver has the knowledge of the users within the cell, while the spreading sequences of users from other cells are unknown. In group-blind multiuser detection, it is assumed that the given spreading sequences of the first Kˆ  users,  KK dˆ , are known to the receiver. Denote ],,[1ˆ0  KddD   and T]0,0,1[0  f ,where 1ˆ0u Kf . Without loss of generality, the subscript pwill be dropped from the model in Eq. (3), and user 1 is the desired user. Then the group-blind linear MMSE detector for user 1 can be formulated as the following multiple constrained optimization problem: THHwtsnnbPEJ 0200 ..},|)()({|min)( fDwyww   . (6)By the method of Lagrange multipliers, the solution of Eq. (6) can be written as [4] 0111 )( fDUȁUDDUȁUw   HsssHHsss . (7)If only the spreading sequence of desired user is known to the receiver, it has 0dD   and 10  f . Consequently, the subspace-based blind MMSE detector is formulated as  1..},|)()({|min)( 0200   dwywwHHwtsnnbPEJ . (8)Then the linear detector w  can be solved as [3] 010101dUȁUdUȁUdw HsssHsssH . (9)kc1TPTS  12P)2cos( 1tf PS)2cos( 1)1( tf MP Skc2kMc)2cos( 1tfS¦s (t)IDFTother usersS/P)(TDataFig. 1. Transmitter structure for MC-CDMA.,,,A blind subspace-based channel estimation algorithm is developed as follows. Consider a channel with L  multipaths, and kg  is the channel impulse response of the kth  user. It has kmk gFh  , where LMmuCF  is a DFT matrix. We then have kmkk gFCd ˆ , where kC  is the code matrix, and kgˆ  denotes the estimation of kg . The orthogonality between the noise and signal subspace implies 0  kmkHnkHn gFCUdU . Then the channel vector kg  can be estimated by least-squares approach: kkmkHnnHkHmHkLkk gFCUUCFggȍCg  	 }{minargˆ1u .(10)It is seen that kg  is the eigenvector corresponding to the smallest eigenvalue of matrix kȍ . If ],,,,[ 1ˆ210 Km dddFC has full column rank, the channel estimation is unique and exact up to a complex constant. The residual scalar ambiguity inherent in blind algorithms can be resolved by transmitting short pilot symbols or using differential encoding/decoding.   4. BI-SVD SIGNAL SUBSPACE TRACKING We adopted Bi-SVD and PASTd subspace tracking algorithms to recursively estimate signal subspace. It is found that the convergence of the PASTd is rather slow. While the Bi-SVD subspace tracker presents much faster convergence, and its performance is very close to the solution using exact SVD. In this section, Bi-SVD algorithm is presented for signal subspace tracking. Let a Mn u )1( data matrix )]()0([)( nnH yyY  denote received signals at the first n  data samples. Define a )1()1( u nn  matrix }1,,{)( 1  nndiagn OOȁ , then the exponentially windowed data sample covariance matrix yR  is  ),()()]()([)]()([)()()(0nnnnnniinHHHniinXXYȁYȁyyR y  ¦  O(11)where»»¼º««¬ª )1()()(nnnHXyXO. Let»»¼º««¬ª 0ǿQ srA )0( . Consider the following bi-iteration applied on the data matrix )(nX : [8] .:)()()(),1()()(,)1(:)()()(),1()()(ionfactorizatQRrMnnnnnnionfactorizatQRrnnnnnnnsAABHsBBAu  u  RQAQXARQBQXB(12)It can be shown that )(nAQ will converge toward the srdominant right eigenvectors of )(nX , which is identical to the signal subspace basis matrix of the covariance matrix yR .However, the vertical dimension of the data matrix )(nX keeps increasing during the updating process. In order to eliminate the growing matrices in Eq. (12), a key step in the bi-iteration subspace tracker is to approximate the data matrix )(nX with a low-rank approximation )(~nX , i.e., )1()()()(~  nnnn HABB QRQX . (13)This approximation leaves the auxiliary recursion matrix )(nBin Eq. (12) unaltered: )1()(~)1()()(   nnnnn AA QXQXB .A square-root bi-iteration subspace tracking algorithm is summarized in Table 1. More details about the bi-iteration subspace tracking algorithm can be referred to [8]. TABLE 1 A SQUARE-ROOT BI-ITERATION SUBSPACE TRACKER)).()()()(()();()()()1()(),()(:*)(*)()(,)()()()()()1()()(),(:)()()(),(:)()()(,)()()()(),1()1()(),(/)()(,)()(),()1()()(),()1()(:;)1()1(;)1(:2nnnndiagnnnnnnnandnextractnnnnnennnnnnnsolvetoonsubstitutibacknnnnsolvetoonsubstitutibacknnnnnnnnnnnennnnennnnnnnUpdatennnInitializeHBHAABsHAAAHAHAAHRyHRRAAARBHRHRHBHRHBBAByyAHArABrA ssRșșRȁfyșQQfșfșGhhhHRG0RHHRHhhRhHhG0RQRHyyyhQyyyQhǿșR0ǿQ  »»¼º««¬ª »»¼º««¬ª  »¼º«¬ª  »»¼º««¬ª »¼º«¬ª       »»¼º««¬ª O5. NUMERICAL RESULTS Numerical simulations have been implemented to illustrate the performance of the group-blind and blind multiuser detector in MC-CDMA systems. We consider a synchronous MC-CDMA system in uplink channels with random spreading sequence. The input SNR is defined as the average transmitted bit energy to noise ratio 2/ vkP V . The multiple access interferers (MAI’s) are 10-dB, i.e., 10/ 0  PPk . The scalar ambiguity of the blind channel estimator is compensated with a complex constant, which makes the first element of the estimation p,0hˆ  and that of the true channel vector p,0h  identical. First, the steady-state performance of the detectors is investigated. In the case of static channel, we set 32 M  and 3 P . There are 10 K  users, of where 7ˆ  K  are known users. A multipath channel with six paths is considered. The channel coefficients are randomly generated according to a complex Gaussian distribution. The simulation results are averaged over 200 Monte-Carlo trials. Fig. 2(a) and (b) illustrate the steady-state signal-to-interference-plus ratio (SINR) of the subspace-based group-blind, blind MMSE receivers and the training MMSE detector versus data samples and the input SNR. Assume that the transmitting symbols among the users are zero mean and i.i.d. sequences. Consequently, we have    ¦    112000 )(KkinsnHkkkH PP RRǿddddR y V . (14),,,The output SINR is defined as )/()( wRwwRw inHsHSINR  .It shows that the group-blind detection substantially outperforms the blind receiver in terms of output SINR. In the simulations, the initial estimation of the signal subspace required for the PASTd is obtained by applying an SVD to the first 50 data vectors. It can be seen that the convergence of the PASTd is rather slow even with this initialization. Whereas, the Bi-SVD presents much faster convergence, and its performance is very close to that of the detectors using an exact SVD.    Next, a slowly time-varying channel with 6 L  Rayleigh fading paths in indoor environments is considered, where the multipath intensity profile decays exponentially. The maximum Doppler frequency ( Df ) is 52 Hz corresponding to a terminal speed smv /3  with a carrier frequency of 5.2 GHz. The subcarrier spacing ( f' ) is 312.5 KHz, and the time duration at the subcarrier is sTs P2.3 . Hence, the normalized fading rate is 410664.1 u sDTf . Due to the slow convergence of the PASTd, the Bi-SVD is used for signal subspace tracking. We average simulation results over 500 Monte-Carlo trials. Fig. 3 depicts the BER comparison for the subspace-based group-blind and blind receivers. To save simulation time, we set 16 M  and 8 P . There are 7 K  users, of which 5ˆ  K  are known users. Again, it is seen that the group-blind detector performs much better than the blind linear detector. Moreover, the performance of the blind multiuser detectors with signal and noise subspace tracking is very close to that of the receivers using an exact SVD.  6. CONCLUSIONS In this paper, a group-blind detector for MC-CDMA systems is studied. By imposing linear vector constraints, the group-blind detector zero-forces the interferers within the cell, and suppresses the unknown interferes from other cells by MMSE criteria.  Hence, it considerably outperforms the subspace-based blind detector. Signal and noise subspace tracking algorithms are investigated to recursively updating the coefficients of the group-blind and blind detectors. Numerical results show that the group-blind detector performs much better than the blind detector, and the performance of detectors using the Bi-SVD is very close to that of detectors using an exact SVD.  10 15 20 25 3010-410-310-210-1SNR(dB)BERSub(Bi-iteration)Sub(SVD)GSub(Bi-iteration)GSub(SVD)Fig. 3. BER comparison for the group-blind and blind detectors. 16 M , 8 P , 7 K / 5ˆ  K , 96.0 O , and 1.0 G .REFERENCES [1] N. Yee, J. P. Linnartz, and G. Fettweis, “Multicarrier CDMA in indoor wireless radio networks”, PIMRC, pp. 109-133, vol. 1, 1993.[2] S. Hara and R. Prasad, “Overview of multicarrier CDMA”, IEEE Commun. Magazine, vol. 35, pp. 126-133, Dec. 1997.[3] X. D. Wang and H. Vincent, “Blind multiuser detection: A subspace approach”, IEEE Trans. Inform. Theory, vol. 44, pp. 677-690, March 1998.[4] X. D. Wang and A. Host-Madsen, “Group-blind multiuser detection for uplink CDMA”, IEEE J. Select. Areas Commun., vol. 17, pp. 1971-1984, Nov. 1999. [5] J. Namgoong, T.F. Wong and J.S. Lehnert, “Subspace MMSE receiver for multicarrier CDMA”, IEEE WCNC '99, pp. 90-94, New Orleans, LA, 1999. [6] W.  Sun and H. B. Li, ““Blind channel identification for multi-carrier CDMA systems with transmit diversity”, ICC'02, pp. 727-731, NY, 2002.” [7] P. Comon and G. H. Golub, “Tracking a few extreme singular values and vectors in signal processing”, Proc. IEEE, vol. 78, pp. 1327-1343, Aug. 1990. [8] P. Strobach, “Bi-iteration SVD subspace tracking algorithms”, IEEE Trans. Signal Processing, vol. 45, pp. 1222-1240, May. 1997. [9] B. Yang, “Projection approximation subspace tracking”, IEEE Trans. Signal Processing, vol. 44, pp. 95-107, Jan. 1995. [10] P. Strobach, “Square-root QR inverse iteration for tracking the minor subsapce”, IEEE Trans. Signal Processing, vol. 48, pp. 2994-2999, Nov. 2000. 500 1000 1500 2000 250046810121416182022Number of samplesAverage SINR(dB)Training MMSEGroupSub(SVD)GroupSub(Bi-SVD)GroupSub(PASTd)Sub(SVD)Sub(Bi-SVD)Sub(PASTd)(SNR = 22 dB/ 10 users) 10 15 20 25 3051015202530SNR(dB)Average SINR(dB)Training MMSEGroupSub(SVD)GroupSub(Bi-SVD)GroupSub(PASTd)Sub(SVD)Sub(Bi-SVD)Sub(PASTd)Fig. 2. Output SINR for the group-blind, blind detectors and training MMSE detector versus: (a) data samples. (b) input SNR. 32 M , 3 P , 10 K / 7ˆ  K , 995.0 O , and 1.0 G .,,,

A group-blind multiuser receiver for MC-CDMA systems

http://hub.hku.hk/bitstream/10722/46432/1/90064.pdf

A group-blind multiuser receiver for MC-CDMA systems

Abstract

Similar works

Full text

Available Versions

HKU Scholars Hub