Abstract:
A simplified de-correlation method in CDMA multi-user detection comprises: a. receive wireless symbols S; b. obtain a channel correlation matrix R, take a part from R to get a partial correlation matrix R P ; c. do inversion operation to the partial correlation matrix R P , then obtain matrix V (m) ; d. recover original data symbols D from received symbols S by V (m)  that the location of original data symbols D corresponds to.

Description:
FIELD OF THE INVENTION  
       [0001]     The invention relates to a de-correlation method and relative device in CDMA, especially relates to a simplified de-correlation method and relevant device in TD-SCDMA multi-user detection.  
       BACKGROUND OF THE INVENTION  
       [0002]     Under the collective influence of channel multi-path dispersion, Inter-Symbol Interference and Multiple Access Interference, it becomes a key of that the multi-user detection is practical or not, that how to depress the complexity of the multi-user detection cancellation algorithm, namely the de-correlation method, to acceptable degree. A received symbol vector can be expressed as the sum of a noise vector and the product that a transmitted symbol vector D left multiply by a channel correlation matrix R. Among the algorithm of demodulating the received symbol vector and getting the estimation {circumflex over (D)} of the transmitted symbol, the most complex step is getting the inverse matrix R −1  of the channel correlation matrix R. For example, usually existing technology adopts the following method: in downlink spreading factor is SF=16. Assuming there are K users in one time slot and 1≦K≦16, in each of the data part, there are K*N data symbols altogether, each user has N=22 data symbols. Let the n th  transmitted data symbols of all the K users be 
 
D (n) =(d 1   (n) ,d 2   (n) ,d 3   (n) , . . . d k   (n) ), n=1, . . . ,K  (1) 
 
         [0003]     The symbol sequence of each user consists of N elements with intervals Ts. The elements are taken from a complex alphabet (1,j,−1,−j). Each data symbols of user K is multiplied by the user specific signature sequence. The specific signature sequence is expressed as: 
 
C (k) =(c 1   (k) , c 2   (k) , . . . , c Q   (k) ) T , k=1, . . . ,K  (2) 
 
         [0004]     Here [·] T  means vector/matrix transposition. The user specific signature sequence consists of Q chips at chip interval T c that is equal to Ts/Q. Each of the K channels is characterized by its discrete impulse response. The discrete impulse response of the k th  channel is expressed as: 
 
H (k) =(h 1   (k) , h 2   (k) , . . . , h w   (k) ) T , k=1, . . . ,K  (3) 
 
         [0005]     The discrete impulse response consists of W samples at chip rate interval Tc. The channel impulse response is assumed to be unchanged during one time slot.  
         [0006]     It is very obviously that ISI arises when W&gt;1, at one time, MAI arises due to channel distortion and non-orthogonal symbol C (k) . The combined channel impulse response is defined by the convolution: 
 
 B   (k) =( b   1   (k)   ,b   2   (k)   , . . . ,b   Q+W−1   (k) ) T   =C   (k)   *H   (k)   , k= 1 , . . . ,K   (4) 
 
         [0007]     The received sequence e is obtained at the chip rate. It is a sum of K sequences, each of length (N·Q+W−1), that are assumed to be arrive synchronously, perturbed by an noise sequence with same length: 
 
 n =(n 1 ,n 2   , . . . , n   N·Q+W−1 ) T   (5) 
 
         [0008]     The received sequence can be written as 
 
 e =( e   1   , e   2   , . . . e   N·Q+W−1 ) T   =A·D+n   (6) 
 
         [0009]     With transmitted data vector 
 
D=(D (1) , D (2) , . . . , D (N) ) T   ={d   1   (1)   ,d   2   (1)   , . . . d   K   (1)   , . . . d   1   (n)   ,d   2   (n)   , . . . d   1   (N)   ,d   2   (N)   , . . . ,d   K   (N) } T =( d   1   , d   2   , . . . , d   KN ) T   (7) 
 
         [0010]     Where d 1   (n) , d 2   (n) , . . . d k   (n)  is the n th  symbol of all the K users, and 
 
 d   j   def   =d   k   (n)   ; j=k+K ·( n− 1), k=1 , . . . K, n= 1  . . . N   (8) 
 
         [0011]     With the matrix  
                   A   =     (     a   ij     )       ;     i   =     1   ⁢     …   ⁡     (       N   ·   Q     +   W   -   1     )             ,     j   =     1   ⁢           ⁢   …   ⁢           ⁢     K   ·   N           ⁢     
     ⁢       a         Q   ⁡     (     n   -   1     )       +   1     ,     k   +     K   ⁡     (     n   -   1     )             =     {             B   l     (   k   )       ,           when   ⁢     {               n   =     1   ⁢           ⁢   …   ⁢           ⁢   N       ,     k   =     1   ⁢           ⁢   …   ⁢           ⁢   K         ⁢                       1   =     1   ⁢           ⁢   …   ⁢           ⁢     (     Q   +   W   -   1     )                             0   ,         other                     (   9   )             
 
         [0012]     The received sequence e has to be processed to obtain a decision on the transmitted data symbol D under the assumption that the user-specific signature sequences C (k)  and the channel cross-correlation H (k) , k=1 . . . K are known at receivers.  
         [0013]     At the same time zero-forcing method is used to eliminate MAI and ISI. It is based on minimizing 
 
∥A{circumflex over (D)}−e∥ 2   (10) 
 
         [0014]     Here {circumflex over (D)} is the estimation of transmitted symbols of all the K users. So 
 
 {circumflex over (D)} =( A   H   A ) −1   A   H   e=D +( A   H   A ) −1   A   H   n   (11) 
 
         [0015]     Here [·] H  means Hermit transposition and A H  e is the output of match filter, (A H A) −1 A H n is noise term. Above {circumflex over (D)} contains desired output D and noise, without MAT and ISI. Its covariance matrix δ 2 (A H A) −1  gives the correlation of the noise term. Generally, the variance of the noise term is more than the noise term which is obtained because of using match filter method. The SNR per-symbol at the output of the de-correlation is equal to 
 
γ( k, n )=/δ 2 [( A   H   A ) −1 ] i,j   j=n+N ·( k− 1),  k= 1  . . . K, n= 1  . . . N   (12) 
 
 Let 
 
R=A H A  (13) 
 
         [0016]     Then (Eq.11) becomes 
 
 {circumflex over (D)} =( R ) −1 ( RD+A   H   n )= D +( R ) −1   A   H   n   (14) 
 
         [0017]     Because the difficulty exists in the matrix inversion operation of R, the calculation burden in TD-SCDMA system is bigger.  
       SUMMARY OF THE INVENTION  
       [0018]     The invention wants to resolve the technical problem that provides a simplified de-correlation method to greatly reduce the complexity degree in multi-user detection calculation.  
         [0019]     For resolving above problem, the invention adopts following technical project:  
         [0020]     Provided de-correlation contains following steps:  
         [0021]     a: Receive wireless symbols S;  
         [0022]     b: Obtain a channel correlation matrix R, take one part from R and get a partial correlation matrix R P ;  
         [0023]     c: Do inversion operation of the partial correlation matrix R P  to obtain matrix V (m) ;  
         [0024]     d: Recover transmitted data symbols D respectively from received symbols S according to V (m)  that the corresponding location of the transmitted data symbols D corresponds to.  
         [0025]     The advice according to above method includes  
         [0026]     relative calculation advice to define the partial correlation matrix R P ;  
         [0027]     matrix inversion device and drawing out device to define new matrix V (m) ;  
         [0028]     matrix-vector multiplication to multiply received wireless symbols S by said matrix V (m) , which are connect each with other in turn.  
         [0029]     Based on the character of channel correlation matrix, the invention has de-modulated only according to the inverse matrix of the diagonal partial matrix R p  of the correlation matrix R so as to reduce greatly the calculation in obtaining inverse matrix and reduce the loss of de-modulate performance. Assuming p=2, the most delay-time is 4 chips, multi-path is 4, when total users are 12, the performance is same as the un-simplified. Because the calculation complexity of the matrix with size M*M is directly proportional to the M 3 , so when P=22, N=22, the simplified inversion has only (5/22) 3  times calculation complexity of the original un-simplified matrix inversion. This is quite a great reduction. Equally, in the symbol recover calculation, the complexity is also reduced to (5/22) of the original operation. 
     
    
     BRIEF DESCRIPTION OF THE DRAWING FIGURES  
       [0030]      FIG. 1 . is adopted advice according to the simplified de-correlation method of the invention. 
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0031]     Original transmitted data symbol can be expressed as  
         [0000]     D=[ . . . d 1   (n−p) , d 2   (n−p) , . . . , d k   (n−p) , . . . , d 1   (n) , d 2   (n)  . . . , d k   (n)  . . . , d 1   (n+p) , d 2   (n+p) , . . . ,d k (n+p) . . . ] 
         [0032]     Above d 1   (n−p) , d 2   (n−p) , . . . , d k   (n−p)  means (n−p) th  symbol of all K users. Above d 1   (n) , d 2   (n)  . . . , d k   (n)  means n th  symbol of all K users. Above d 1   (n+p) , d 2   (n+p) , . . . , d k   (n+p)  means (n+p) th  symbol of all K users.  
         [0033]     Received symbol S can be expressed as  
         [0034]     S=[ . . . ŝ 1   (n−p) , ŝ 2   (n−p) , . . . , ŝ k   n−p) , . . . , ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ k   (n)  . . . , ŝ 1   (n+p) , ŝ 2   (n+p) , . . . , ŝ k   (n+p)  . . . ]. Above ŝ 1   (n−p) , ŝ 2   (n−p) , . . . , ŝ k   (n−p)  means (n−p) th  symbol of all K users. Above ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ k   (n)  means n th  symbol of all K users. Above ŝ 1   (n+p) , ŝ 2   (n+p) , . . . , ŝ k   (n+p)  means (n+P) th  symbol of all K users. Let us define ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ k   (n)  as S (n) . Since the multi-path spread in TD-SCDMA is assumed to be less than 12 chips, and the spreading factor for downlink is 16, so there are no ISI from symbols more than 2 symbols away. Also the MAI and ISI part from any particular other user symbol is always much less than desired user symbol amplitude. To recover D (n) =d 1   (n) , d 2   (n)  . . . , d k   (n) , influence from symbols in S (t)  (t&lt;n−p, or t&gt;n+p) can be neglected. It means D (n)  can be recovered from received symbols  
         [0035]     S p   (n) =[ŝ 1   (n−p) , ŝ 2   (n−p) , . . . , ŝ k   (n−p) , . . . , ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ k   (n)  . . . , ŝ 1   (n+p) , ŝ 2   (n+p) , . . . , ŝ k   (n+p) ]. Above ŝ 1   (n−p) , ŝ 2   (n−p) , . . . , ŝ k   (n−p)  means (n−p) th  symbol of all K users. Above ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ k   (n)  means n th  symbol of all K users. Above ŝ 1   (n+p) , ŝ 2   (n+p) , . . . , ŝ k   (n+p)  means (n+p) th  symbol of all K users. Assume P is large enough.  
         [0036]     We here define a partial correlation matrix 
 
 R   p   ={r   i,j   }, i,j= 1 . . . (2 P+ 1) K  
 
         [0037]     Here r i,j  is the element of R in the i th  row and j th  column.  
         [0038]     Define a new matrix V (m) ={v i,j   (m) }
 
 v   i,j   (m) =( R   −1   p ) i+(m−1)K,j , i=1 . . . K, j=1 . . . (2P+1)K, m=1 . . . 2P+1 
 
         [0039]     Here R −1   p  means matrix inversion. When P+1≦n≦N−P, D (n)  can be recovered as 
 
 {circumflex over (D)}   (n)   =V   (P+1)   S   p   (n)  
 
         [0040]     {circumflex over (D)} (n)  is the estimation of transmitted symbol. When n≦P, i.e. the interference symbol before D (n)  is less than P, D (n)  can be recovered as 
 
 {circumflex over (D)}   (n)   =V   (n)   S   p   (n) . 
 
         [0041]     Here S p   (n)  is defined as 
 
 S   p   (n)   =[ŝ   1   (1) , ŝ 2   (1) , . . . , ŝ K   (1) , . . . , ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ K   (n)  . . . , ŝ 1   (2P+1) , ŝ 2   (2P+1) , . . . ,ŝ K   (2P+1) ]. 
 
         [0042]     Above ŝ 1   (1) , ŝ 2   (1) , . . . , ŝ k   (1)  is first symbol of all K users. Above ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ k   (n)  is n th  symbol of all K users. Above ŝ 1   (2P+1) , ŝ 2   (2P+1) , . . . , ŝ k   (2P+1)  is 2P+1 th  symbol of all K users.  
         [0043]     When n≧N+1−P, and the interference symbols after D (n)  is less P, D (n)  can be recovered as 
 
 {circumflex over (D)}   (n)   =V   (2P+1+n−N)   S   p   (n) . 
 
         [0044]     Here S p   (n)  is defined as  
         [0045]     S p   (n) =[ŝ 1   (N−2P) , ŝ 2   (N−2P) , . . . ,ŝ k   (N−2P) , . . . , ŝ 1   (n) ,ŝ 2   (n)  . . . ,ŝ k   (n)  . . . ,ŝ 1   (N),ŝ   2   (N) , . . . ,ŝ k   (N) ].  
         [0046]     Above ŝ 1   (N−2P) , ŝ 2   (N−2P) , . . . , ŝ K   (N−2P)  is N−2P th  symbol of all K users. Above ŝ 1   (n) , ŝ 2   (n)  . . . , ŝ K   (n)  is n th  symbol of all K users. Above ŝ 1   (N) , ŝ 2   (N) , . . . , ŝ K   (N) ], is N th  symbol of all K users.  
         [0047]     In addition, the simplified de-correlation method can be named as partial de-correlation whose meaning is that use de-correlation method to eliminate the influence of the symbols which are earlier or later than the current on current symbols only according to P group of symbols which are earlier or later than current symbols.  
         [0048]      FIG. 1  is the structure sketch map of the adopted advice according to the simplified de-correlation method.  
         [0049]     It includes K match filters  1 , K buffer storage  2  which are correspondingly connected with K match filters, relative calculation device  3  to define partial correlation matrix R p , matrix inversion device  4  and drawing out device  5  which tag after the relative calculation device to define new matrix V (m) . The output of the mentioned buffer storage  2  and the output of the drawing out device  3  are connected with the input of the matrix-vector multiplication  6 .  
         [0050]     The mentioned baseband symbols of K users respectively are inputted to K match filters  1  and K buffer storage  2  to get received symbol S p   n .  1  and  2  are connected with each other correspondingly The mentioned user signature sequence C (k) =(c 1   (k) ,c 2   (k) , . . . C Q   (k) ) T , channel impulse response H (1) ,H (2) , . . . ,H (k) , channel specific multiplication factor is respectively inputted to the relative calculation device  3  to obtain partial correlation matrix, then obtain the inverse matrix V of matrix R p  through matrix inversion device  4  and drawing out device  5 , and draw out partial matrix V (m)  from it. Through matrix vector multiplication  6  symbol estimation value {circumflex over (D)} (n)  can be obtained from mentioned received symbol S p   (n)  and V (m) .  
         [0051]     The present invention is not restricted to above advice and method, so all the technical changes known by ordinary technician of the present fields should fall into the protective scope of the present invention.