Abstract:
An adaptive interference suppression receiver for a DS/CDMA system includes: input signal generators outputting a part of received signals of receiver antennas that has a coding block&#39;s transfer code; adaptive filter blocks filtering a complex signal output from the input signal generators with a tab weight updated based on a constraint MMSE criterion and estimating phase and amplitude components of a specific user channel using the filtering output signal to generate a channel-estimated signal; a signal restorer combining the channel-estimated signal and the filtering output signal and restoring an original signal; a signal selector for selecting known training data signal or the restored signal; and error generators comparing the channel-estimated signal and the filtering output signal and calculating an error between the two signals.

Description:
BACKGROUND OF THE INVENTION  
         [0001]    (a) Field of the Invention  
           [0002]    The present invention relates to a receiver for a MIMO (Multiple-Input Multiple-Output) antenna CDMA (Code Division Multiple Access) system. More specifically, the present invention relates to an adaptive MMSE (Minimum Mean Square Error) interference suppression receiving apparatus for a space-time block coded DS (Direct Sequence)/CDMA communication system.  
           [0003]    (b) Description of the Related Art  
           [0004]    In general, the MIMO antenna system is a communication system that transmits space-time block coded data via multiple transmitter antennas after data spreading with a different code by users, and receives data from users via multiple receiver antennas.  
           [0005]    The conventional technologies concerning the structure of a receiver for the space-time block coded MIMO antenna system are disclosed in the following documents: Ref. 1: S. M. Alamounti, A simple transmitter diversity scheme for wireless communications, IEEE J. Select. Atras Commun., vol. 16, pp.1451-1458, October 1998; Ref. 2: V. Tarokh, H. Jafarkhani, and A. R. Calderbank, Space-time block codes from orthogonal designs, IEEE Trans. Inform. Theory, vol. 45, pp.1456-1467, July 1999; and Ref. 3: A. F. Naguib, N. Seshadri, and A. R. Calderbank, Increasing Data Rate Over Wireless Channels, IEEE Signal Processing Magazine, vol.17, No.3, pp.76-92, May 2000.  
           [0006]    Ref. 1 and Ref. 2 disclose a space-time block coded TDMA system using two transmitter antennas, an optimal block coding transmission method for a system using at least two transmitter antennas, and a maximum likelihood receiver for the same.  
           [0007]    The maximum likelihood receiver proposed in the two references can be applied for space-time block coded DS/CDMA systems and can provides an antenna diversity gain proportional to the product of the number of transmitter antennas and the number of receiver antennas, thereby guaranteeing an excellent data detection performance relative to general Rake receivers.  
           [0008]    But the maximum likelihood receiver, which is devised on the assumption that the interference signal and the background noise other than the desired user&#39;s transfer signals are AWGN (Additive White Gaussian Noise), has a deterioration of performance in a MAI (Multiple Access Interference) environment.  
           [0009]    On the other hand, the receiver of Ref. 3 is designed to enhance a reception performance by eliminating interference signals of other users under a circumstance in which multiple users are present in the TDMA system after space-time block coding. But the receiver has a loss of receiver antenna diversity in the order of as much as the number of interference signals to be removed, causing a deterioration of reception performance.  
           [0010]    The above reception methods for the space-time block coded TDMA communication system can be expanded to the reception method for the space-time block coded DS/CDMA communication system with a little modification. However, the maximum likelihood receiver has a deterioration of performance due to the MAI that inevitably occurs in the CDMA environment when multiple users share a single channel, and the interference suppression receiver of Ref. 3 has a deterioration of data detection performance caused by a loss of diversity gain.  
         SUMMARY OF THE INVENTION  
         [0011]    It is an object of the present invention to solve the problems with the prior art and to provide an adaptive MMSE interference suppression receiver that has an excellent data detection performance even in the multipath fading channel environment, which is an actual mobile communication channel environment, in the space-time block coded DS/CDMA system.  
           [0012]    In one aspect of the present invention, there is provided an adaptive interference suppression receiver for a DS/CDMA communication system, which receives via a multipath fading channel a space-time block coded signal spread with another spreading code by users and despreads the received signal, the adaptive interference suppression receiver including: at least two input signal generators coupled to at least two receiver antennas for signal reception, respectively, the input signal generators outputting a part of the received signals of the respective receiver antennas that has a transfer code of a corresponding coding block; adaptive filter blocks coupled to the at least two input signal generators, respectively, the adaptive filter blocks filtering a complex signal output from the input signal generators with a tab weight updated based on at least one constraint MMSE (Minimum Mean Square Error) criterion and estimating phase and amplitude components of a specific user channel using the resultant filtering output signal to generate a channel-estimated signal, the number of the adaptive filter blocks being the same as the number of multipaths per one input signal generator; a signal restorer for combining the channel-estimated signal and the filtering output signal received from the adaptive filter blocks for all multipath components and restoring an original signal that is intended to be transferred by the specific user; a signal selector for selecting either one of a known training data signal or the restored signal of the signal restorer; and error generators for comparing the channel-estimated signal and the filtering output signal received from the adaptive filter blocks with the signal from the signal selector and calculating an error between the two signals for update of a tap weight in the adaptive filter blocks, the number of the error generators being the same as the number of the adaptive filter blocks.  
           [0013]    The adaptive filter blocks include: an adaptive filter for filtering the output signal of the input signal generators according to the tap weight; a simple mean channel estimator for performing simple mean channel estimation using the filtering output signal of the adaptive filter and the output signal of the signal selector; and a multiple constraint criterion filter coefficient updater for regulating the tap weight of the adaptive filter based on at least two constraint MMSE criteria using the error signal of the error generators and the output signal of the adaptive filter.  
           [0014]    Here, the at least two constraint MMSE criteria substantially constrain the inner product of the tap weight w ml1 ,w ml2  of the adaptive filter block for the l-th multipath component of the m-th receiver antenna and a spreading code vector c l  of the corresponding multipath component to “1” and the inner product of the tap weight w ml1 ,w ml2  and a spreading code vector c i (i≠l) of another multipath component to “0”, thereby minimizing the error calculated by the error generators, the at least two constraint MMSE criteria being defined as the following equations:  
           minimize     w     m                 l                 1              E        [                    h   ^       1      m                 l            (   n   )              b   1          (   n   )         +           h   ^       2                 m                 l            (   n   )              b   2          (   n   )         -       w     m                 l                 1     H            r     m                 1            (   n   )                2     ]                     subject                 to                   C   H          w     m                 l                 1         =     g   l                 minimize     w     m                 l                 2              E        [                    h   ^       2      m                 l            (   n   )              b   1   *          (   n   )         -           h   ^       1                 m                 l            (   n   )              b   2   *          (   n   )         -       w     m                 l                 2     H            r     m                 2            (   n   )                2     ]                     subject                 to                   C   H          w     m                 l                 2         =     g   l                           
 
           [0015]    wherein E is a mean value; ĥ 1ml (n) is a channel estimated for the l-th multipath component between the first transmitter antenna and the m-th receiver antenna; ĥ 2ml (n) is a channel estimated for the l-th multipath component between the second transmitter antenna and the m-th receiver antenna; b 1 (n),b 2 (n) is the output signal of the signal selector for a corresponding coding block; r m1 (n) is a vector corresponding to an interval of receiving odd transfer symbols of the n-th coding block; r m2 (n) is a vector corresponding to an interval of receiving even transfer symbols of the n-th coding block; g 1  is a vector having a size of L×1 where the l-th term is “1” and the others are all “0”, the vector being a constraint criterion added to extract only properly synchronized signals by the respective multipaths from a frequency-selective fading channel and to remove the effect of the desired user&#39;s other multipaths; C=[c 1  c 2  . . . c L ]; and the subscript H represents a Hermitian operation.  
           [0016]    The adaptive filter blocks update the tap weight w ml (n) that combines errors for the respective symbols according to the following equation and satisfies at least one constraint criterion, based on that the tap weight w ml (n)w ml1 (n),w ml2 (n) satisfying the MMSE criterion has a same solution:  
             w   ml ( n+ 1)= w   ml ( n )+μ P   C   ⊥   [e   ml,1 *( n ) r   ml,1 ( n )+ e   ml,2 *( n ) r   ml,2 ( n )] 
           [0017]    wherein P C   ⊥ =I−C(C H C) −1 C H  represents an orthogonal projection matrix that projects all the vectors on a space orthogonal to C;  
             e   m/1   =ĥ   1ml ( n ) b   1 ( n )+ ĥ   2ml ( n ) b   2 ( n )− w   ml   H   r   m1 ( n )  
             e   m/2   =ĥ   2ml ( n ) b   1 *( n )− ĥ   1ml ( n ) b   2 *( n )− w   ml   H   r   m2 ( n )  
           [0018]    wherein μ is a step size that determines how fast the tap weight is changed; and the subscript * represents a complex conjugate operation.  
           [0019]    The simple mean channel estimator estimates, for a predetermined number of received pilot symbols, a channel of the l-th multipath component between each transmitter antenna and the m-th receiver antenna by multiplying the filtering output signal of the adaptive filter by a complex conjugate of the output signal of the signal selector and then averaging the multiplication results according to the following equation:  
                   h   ^       1      m                 l            (   n   )       =       1     2        N   P                ∑     i   =   1       N   P                           w     m                 l     H          (     n   -     i                 Q       )            (           d   1   *          (     n   -     i                 Q       )              r     m                 1            (     n   -     i                 Q       )         -         d   2   *          (     n   -     i                 Q       )              r     m                 2            (     n   -     i                 Q       )           )                             h   ^       2      m                 l            (   n   )       =       1     2        N   P                ∑     i   =   1       N   P                           w     m                 l     H          (     n   -     i                 Q       )            (           d   2   *          (     n   -     i                 Q       )              r     m                 1            (     n   -     i                 Q       )         +         d   1   *          (     n   -     i                 Q       )              r     m                 2            (     n   -     i                 Q       )           )                                       
 
           [0020]    wherein N P  is the number of consecutive pilot symbols used for channel estimation; and Q is a pilot symbol insertion period.  
           [0021]    The adaptive filter blocks include: an adaptive filter for filtering the output signal of the input signal generators according to the tap weight; a maximum likelihood channel estimator for performing maximum likelihood channel estimation using the filtering output signal of the adaptive filter and the output signal of the signal selector; and a single constraint criterion filter coefficient updater for regulating the tap weight of the adaptive filter based on a single constraint MMSE criterion using the error signal of the error generators and the output signal of the adaptive filter.  
           [0022]    Here, the single constraint MMSE criterion substantially constrains the inner product of the tap weight w ml1 ,w ml2  of the adaptive filter block for the l-th multipath component of the m-th receiver antenna and a spreading code vector c l  of the corresponding multipath component to “1”, thereby minimizing the error calculated by the error generators, the constraint MMSE criterion being defined as the following equations:  
           minimize     w     m                 l                 1              E        [                  h     1      m                 l            (   n   )              b   1          (   n   )         +         h     2                 m                 l            (   n   )              b   2          (   n   )         -       w     m                 l                 1     H            r     m                 1            (   n   )                2     ]                     subject                 to                   w     m                 l                 1     H          c   l       =   1               minimize     w     m                 l                 2              E        [                  h     2      m                 l            (   n   )              b   l   *          (   n   )         -         h     1                 m                 l            (   n   )              b   2   *          (   n   )         -       w     m                 l                 2     H            r     m                 2            (   n   )                2     ]                     subject                 to                   w     m                 l                 2     H          c   l       =   1                         
 
           [0023]    wherein E is a mean value; h 1ml (n) is a channel estimated for the l-th multipath component between the first transmitter antenna and the m-th receiver antenna; h 2ml (n) is a channel estimated for the l-th multipath component between the second transmitter antenna and the m-th receiver antenna; b 1 (n),b 2 (n) is the output signal of the signal selector for a corresponding coding block; r m1 (n) is a vector corresponding to an interval of receiving odd transfer symbols of the n-th coding block; r m2 (n) is a vector corresponding to an interval of receiving even transfer symbols of the n-th coding block; and the subscript H represents a Hermitian operation.  
           [0024]    The adaptive filter blocks update the tap weight w ml (n) that combines errors for the respective symbols according to the following equation and satisfies the single constraint criterion, based on that the tap weight w ml (n)=w ml1 (n),w ml2 (n) satisfying the MMSE criterion has a same solution:  
             w   ml ( n+ 1)= w   ml ( n )+μ P   c     l     ⊥   [e   ml,1 *( n ) r   ml,1 ( n )+ e   ml,2 *( n ) r   ml,2 ( n )] 
           [0025]    wherein P c     l     ⊥ =I−c l c l   H /(c l   H c l ) represents an orthogonal projection matrix that projects all the vectors on a space orthogonal to C l ; and  
             e   ml1   =ĥ   1ml ( n ) b   1 ( n )+ ĥ   2ml ( n ) b ( n )− w   ml   H   r   m1 ( n )  
             e   ml2   =ĥ   2ml ( n ) b   1 *( n )− ĥ   1ml ( n ) b   2 *( n )− w   ml   H   r   m2 ( n )  
           [0026]    wherein μ is a step size that determines how fast the tap weight is changed; and the subscript * represents a complex conjugate operation.  
           [0027]    The maximum likelihood channel estimator estimates, for a predetermined number of received pilot symbols, a channel of the multipath component between each transmitter antenna and the m-th receiver antenna by multiplying the filtering output signal of the adaptive filter by a complex conjugate of the output signal of the signal selector and then averaging the multiplication results according to the following equation:  
         [             h   ^       1      m                 1               h   ^       2      m                 1                   h   ^       1      m                 2               h   ^       2      m                 2               ⋮       ⋮               h   ^       1      m                 L               h   ^       2      m                 L             ]     =       1     2        N   P                ∑     i   =   1       N   p                             [         1             w     m                 1     H          (     n   -     i                 Q       )            c   2           ⋯             w     m                 1     H          (     n   -     i                 Q       )            c   L                     w   m2   H          (     n   -     i                 Q       )            c   1           1       ⋯             w     m                 2     H          (     n   -     i                 Q       )            c   L               ⋮       ⋮       ⋰       ⋮                 w     m                 L     H          (     n   -     i                 Q       )            c   1                 w     m                 L     H          (     n   -     i                 Q       )            c   2           ⋯       1         ]       -   1       ·     W   m   H          A                               
 
           [0028]    wherein  
           W m =[w m1  w m2  . . . w mL ],  
             A=[d   1 *( n−iQ ) r   m1 ( n−iQ )− d   2 *( n−iQ ) r   m2 ( n−iQ )  d   2 *( n−iQ ) r   m1 ( n−iQ )+ d   1 *( n−iQ ) r   m2 ( n−iQ )] 
           [0029]    The signal restorer includes: a coupler for receiving the output signal y ml,1 (n),y ml,2 (n) of the adaptive filter and the resultant signal of channel estimation ĥ 1ml (n),ĥ 2ml (n) and combining them for all the multipath components to generate a combined signal {overscore (d)} 1 (n),{overscore (d)} 2 (n); and a data determiner for restoring an original signal ĥ 1 (n),ĥ 2 (n) that is intended to be transferred by the specific user, using the combined signal {overscore (d)} 1 (n),{overscore (d)} 2 (n) of the coupler according to the following equation:  
                 d   ^     1     =     arg            min             d   ^                         ∑     m   =   1     M                       ∑     l   =   1     L                     [           y       m                 l     ,   1            (   n   )                h   ^       1                 m                 l     *          (   n   )         +         y       m                 l     ,   2     *          (   n   )                h   ^       2                 m                 l            (   n   )           ]         -     d   ^                              d   ^     2     =     arg            min             d   ^                         ∑     m   =   1     M                       ∑     l   =   1     L                     [           y       m                 l     ,   1            (   n   )                h   ^       2                 m                 l     *          (   n   )         -         y       m                 l     ,   2     *          (   n   )                h   ^       1                 m                 l            (   n   )           ]         -     d   _                                          
 
           [0030]    wherein {circumflex over (d)} is a symbol on a constellation corresponding to the combined signal {overscore (d)} 1 (n),{overscore (d)} 2 (n). 
       
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0031]    The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate an embodiment of the invention, and, together with the description, serve to explain the principles of the invention;  
         [0032]    [0032]FIG. 1 illustrates an example of a space-time block coded DS/CDMA system using two transmitter antennas in accordance with an embodiment of the present invention;  
         [0033]    [0033]FIG. 2 is a detailed block diagram of an adaptive interference suppression receiver for the space-time block coded DS/CDMA MIMO antenna system shown in FIG. 1;  
         [0034]    [0034]FIG. 3 is a diagram of the adaptive filter block shown in FIG. 2;  
         [0035]    [0035]FIG. 4 illustrates another example of the adaptive filter block shown in FIG. 2;  
         [0036]    [0036]FIG. 5 is a detailed block diagram of the coupler shown in FIG. 2;  
         [0037]    [0037]FIG. 6 is a diagram of the energy calculating block shown in FIG. 2;  
         [0038]    [0038]FIG. 7 is a detailed block diagram of the data determiner shown in FIG. 2;  
         [0039]    [0039]FIG. 8 illustrates the error generating block shown in FIG. 2;  
         [0040]    [0040]FIG. 9 is a graph showing the result of a performance simulation of the adaptive interference suppression receiver for a space-time block coded DS/CDMA system according to an embodiment of the present invention; and  
         [0041]    [0041]FIG. 10 is a graph showing the result of a performance simulation of the adaptive interference suppression receiver for a space-time block coded DS/CDMA system according to an embodiment of the present invention. 
     
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS  
       [0042]    In the following detailed description, only the preferred embodiment of the invention has been shown and described, simply by way of illustration of the best mode contemplated by the inventor(s) of carrying out the invention. As will be realized, the invention is capable of modification in various obvious respects, all without departing from the invention. Accordingly, the drawings and description are to be regarded as illustrative in nature, and not restrictive.  
         [0043]    [0043]FIG. 1 shows an example of a space/time block coded DS/CDMA system using two transmitter antennas.  
         [0044]    The space/time block coded DS/CDMA MIMO antenna system shown in FIG. 1 performs space/time block coding  10  and spreading  20  of data to be transferred from its transmitter, sends the space/time block coded and spread data via two transmitter antennas  30 , and receives transfer data {circumflex over (d)} 1 (n),{circumflex over (d)} 2 (n) using M (where M is an integer being at least 2) receiver antennas  40  at a receiver  50  and restores them to reception data.  
         [0045]    In this system, the fading channel at time t is given by the following equation:  
               h     (   t   )       =       ∑     i   =   1     2                       ∑     m   =   1     M                       ∑     l   =   1     L                       h       (   t   )     ,     i                 m                 l              δ     (     t   -     τ   l       )                       Equation                 1                               
 
         [0046]    where h (t),iml  is a channel coefficient of the l-th path received on the m-th receiver antenna from the i-th transmitter antenna  40 , and τ l  is a transmission delay of the l-th path and is expediently given by τ l =(l−1)T c  (where T c  is the chip period).  
         [0047]    Assuming that the channel is not changed in the interval of transmitting the respective coding blocks, the channel coefficient is given by h (t),iml =h iml (n) when 2nT=t&lt;(2n+1)T where T is the symbol period.  
         [0048]    [0048]FIG. 2 is a detailed block diagram of the adaptive interference suppression receiver 50 for the space-time block coded DS/CDMA MIMO antenna system shown in FIG. 1.  
         [0049]    The adaptive interference suppression receiver  50  comprises, as shown in FIG. 2, chip-matched filters  210  and  220 ; input signal generating buffers  310  and  320 ; adaptive filter blocks  410 ,  420 ,  430 , and  440 ; error generating blocks  510 ,  520 ,  530 , and  540 ; a coupler  600 ; a data determiner  700 ; and a reference signal selector  800 .  
         [0050]    First, the chip-matched filters  210  and  220  perform matched-filtering of the signals received from the M antennas  40  on the basis of the chip period.  
         [0051]    The input signal generating blocks  310  and  320  generate vectors corresponding to the signal blocks output from the chip-matched filters  210  and  220  and output the generated vectors according to the symbol period.  
         [0052]    The adaptive filter blocks  410 ,  420 ,  430 , and  440  filter the output signals of the input signal generating buffers  310  and  320  according to a tap weight that is regulated with a specific cycle.  
         [0053]    The coupler  600  combines the output signals of the adaptive filter blocks  410 ,  420 ,  430 , and  440  for all the multipath components.  
         [0054]    The data determiner  700  restores the output signal of the coupler  600  to the original signal that is intended to be transferred by a specific user.  
         [0055]    The reference signal selector  800  selectively outputs either the output signal of the data determiner  700  or a known training data signal.  
         [0056]    The error generating blocks  510 ,  520 ,  530 , and  540  use the output signals of the reference signal selector  800  and the adaptive filter blocks  410 ,  420 ,  430 , and  440  and generate an error signal for regulating the tap weight to be used in the adaptive filter blocks  410 ,  420 ,  430 , and  440 .  
         [0057]    Here, the number of the adaptive filter blocks  410 ,  420 ,  430 , and  440  are L×M, where L is the number of the multipath components and M is the number of the receiver antennas  40 .  
         [0058]    Let r m (j) be the output of the chip-matched filters  210  and  220  for the received signals of the m-th receiver antenna  40 , and these signals are fed into the input signal generating buffers  310  and  320 . The input signal generating buffers  310  and  320  generate a vector r m1 (n) corresponding to an interval of receiving odd transfer symbols of the n-th coding block and a vector r m2 (n) corresponding to an interval of receiving even transfer symbols of the n-th coding block from the output signals of the chip-matched filters  210  and  220 , and transmitting the generated vectors to the corresponding adaptive filter blocks  410 ,  420 ,  430 , and  440  according to the symbol period. Let c=[c l  . . . c L ] T  be the normalized spreading code for the desired user&#39;s l-th path and c l =[0 l−1   T c T 0 L−1   T ] T  (where O l  is the I         1 0 vector) be the spreading code for the l-th path, then the received signal matrix r m (n) of the m-th antenna for the n-th coding block can be expressed by the following equation:  
                       r   m          (   n   )       =     [         r     m                 1            (   n   )              r     m                 2            (   n   )         ]                 =     [             r   m          (     2      n                 N     )               r   m          (       2      n                 N     +   N     )               ⋮       ⋮               r   m          (       2      n                 N     +   N   +   L   -   2     )               r   m          (       2      n                 N     +     2      N     +   L   -   2     )             ]                 =           [           c   1         ⋯         c   L           ]            ≡   c                [             h     1      m                 1            (   n   )               h     2      m                 1            (   n   )               ⋮       ⋮               h     Lm                 1            (   n   )               h     Lm                 1            (   n   )             ]            ≡       H   m          (   n   )                    [             d   1          (   n   )             -       d   2   *          (   n   )                     d   2          (   n   )               d   1   *          (   n   )             ]            ≡     D        (   n   )               +       [             u     m                 1            (   n   )               u   m2          (   n   )             ]            ≡       u   m          (   n   )                           =       c                     H   m          (   n   )            D        (   n   )         +       u   m          (   n   )                       Equation                 2                               
 
         [0059]    where u m (n) is a matrix having a size of (N+M−1)         2 that represents a signal impaired by MAI (Multiple-Access Interference), ISI (Inter-Symbol Interference), and AWGN (Additive White Gaussian Noise). Expediently, it is assumed that the earlier multipath component is processed with the previous adaptive filter block (i.e., having the smaller index value).  
         [0060]    The maximum likelihood receiver for the space-time block coded CDMA MIMO antenna system in the frequency-selective fading channel environment can be induced from the receivers of Ref. 1 and Ref. 2 as follows. The matched filter output vector {tilde over (y)} m (n) having a size of 2L         1 for r m (n) can be expressed by the following equation:  
                         y   ~     m          (   n   )       =                [             c   1   H            r     m                 1            (   n   )                     (       c   1   H            r     m                 2            (   n   )         )     *             ⋮               c   L   H            r     m                 1            (   n   )                     (       c   L   H            r     m                 2            (   n   )         )     *           ]            [             h     1      m                 1            (   n   )               h     2      m                 1            (   n   )                   h     2      m                 1     *          (   n   )             -       h     1      m                 1     *          (   n   )                 ⋮       ⋮               h     l                 m                 L            (   n   )               h     2      m                 L            (   n   )                   h     2      m                 L     *          (   n   )             -       h     1      m                 L     *          (   n   )               ]            ≡         H   ~     m          (   n   )                    [             d   1          (   n   )                   d   2          (   n   )             ]            ≡     d        (   n   )               +         u   ~     m          (   n   )                     =                    H   ~     m          (   n   )            d        (   n   )         +         u   ~     m          (   n   )                       Equation                 3                               
 
         [0061]    where ũ m (n) is a signal impaired by AWGN, MAI, and ISI. Working with the normal distribution of ũ m (n), the maximum likelihood receiver can be induced in the same manner as in the TDMA system as expressed by the following equation:  
                   d   ^       M                 L            (   n   )       =     arg            min             d   ^                  ∑     m   =   1     M                       y   ~     m          (   n   )       -           H   ~     m          (   n   )       ·     d   ^              2                 Equation                 4                               
 
         [0062]    To summarize the above cost function, it can be divided into two cost functional equations related to d 1 (n),d 2 (n), respectively, as shown in the following equation, and used for data detection:  
                   d   ^       1   ,     M                 L              (   n   )       =     arg            min             d   ^                         ∑     m   =   1     M                       ∑     l   =   1     L          [         c   1   H            r     m                 1            (   n   )              h     1      m                 l     *          (   n   )         +         (       c   1   H            r     m                 1            (   n   )         )     *            h     2      m                 l            (   n   )           ]         -     d   ^                      Equation                 5                               
 
         [0063]    In case of data despreading with spreading codes as in the maximum likelihood receiver, the despread signals are impaired as ũ m (n) by MAI and inter-path interference in addition to AWGN, causing a performance deterioration.  
         [0064]    The embodiment of the present invention involves despreading using adaptive filter coefficients optimal in the aspect of the CMMSE (Constrained Minimum Mean Square Error) instead of using the spreading codes to effectively eliminate all the interference signals, thus increasing the reception performance.  
         [0065]    Let w ml1  be an adaptive filter coefficient for odd received signals for the l-th multipath component of the m-th antenna and w ml2  be an adaptive filter coefficient for even received signals, then in the first embodiment, the tap weights of the respective adaptive filters are determined so as to minimize the CMMSE cost function of the following equation 6.  
         [0066]    Namely, the inner product of the coefficient w ml1 ,w ml2  of the adaptive filter for the l-th multipath component of the m-th receiver antenna and the spreading code vector c l  of the corresponding multipath component is substantially constrained to “1”, and the inner product of the filter coefficient w ml1 ,w ml2  and the spreading code vector c i (i≠l) of another multipath component is substantially constrained to “0”, thereby minimizing the CMMSE cost product.  
                   minimize     w     m                 l                 1              E        [                    h   ^       1      m                 l            (   n   )              b   1          (   n   )         +           h   ^       2                 m                 l            (   n   )              b   2          (   n   )         -       w     m                 l                 1     H            r     m                 1            (   n   )                2     ]                     subject                 to                   C   H          w     m                 l                 1         =     g   l            
              minimize     w     m                 l                 2              E        [                    h   ^       2      m                 l            (   n   )              b   1   *          (   n   )         -           h   ^       1                 m                 l            (   n   )              b   2   *          (   n   )         -       w     m                 l                 2     H            r     m                 2            (   n   )                2     ]                     subject                 to                   C   H          w     m                 l                 2         =     g   l               Equation                 6                               
 
         [0067]    where E is a mean value; ĥ 1ml (n) is a channel estimated for the l-th multipath component between the first transmitter antenna and the m-th receiver antenna; ĥ 2ml (n) is a channel estimated for the l-th multipath component between the second transmitter antenna and the m-th receiver antenna; b 1 (n),b 2 (n) is an output signal of the reference signal selector  800  for the corresponding coding block; g l  is a vector having a size of L×1 where the l-th term is 1 and the others are all 0, the vector being a constraint criterion added to extract only properly synchronized signals by the respective multipaths from a frequency-selective fading channel and to remove the effect of the desired user&#39;s other multipaths; C=[c 1  c 2  . . . c L ]; and the subscript H represents the Hermitian operation.  
         [0068]    This eliminates the inter-path interference caused by the user&#39;s other multipaths as well as MAI and determines an unbiased channel estimate even by a simple mean channel estimation method using the adaptive filter output, thereby increasing reception performance.  
         [0069]    The optimal filter coefficients w ml1    0  and w ml1    0  that minimize a given constraint cost function can be calculated by the Lagrangian multiplier method according to the following equation:  
           w   ml1   0   =R   r     m1     −1   C ( C   H   R   r     m1     −1   C ) −1   g   l    Equation 7  
           w   ml2   0   =R   r     m2     −1   C ( C   H   R   r     m2     −1   C ) −1   g   l    
         [0070]    where R r     m1   =E[r m1 (n)r ml   H (n)] and R r     m2   =E[r m2 (n)r m2   H (n)] can be rewritten to have the same value as in the following equation, giving the same optimal solution:  
               R     r     m                 1         =       R     r     m                 2         =       R     u   m       +       ∑     l   =   1     L                     E        ⌊         (                h     1      m                 l            (   n   )            2     +              h     2      m                 l            (   n   )            2       )     ·     c   l            c   l   H       ⌋                     Equation                 8                               
 
         [0071]    As the two adaptive filters have the same optimal solution as in the equation 8, it is more effective to average two error values and use the mean value in updating the filter coefficient rather than to independently adapt the filters.  
         [0072]    By using the adaptive filter coefficient w ml (n) for the l-th multipath component of the m-th receiver antenna in filtering for both odd and even symbols, the error values between the adaptive filter outputs of odd and even received signals and the reference signal are properly combined to update the filter coefficient according to the following equation:  
           w   ml ( n+ 1)= w   ml ( n )+μ P   C   ⊥   [e   ml,1 *( n ) r   ml,1 ( n )+ e   ml,2 *( n ) r   ml,2 ( n )]  Equation 9  
         [0073]    where P C   ⊥ =I−C(C H C) −1 C H  represents an orthogonal projection matrix that projects all the vectors on a space orthogonal to C;  
           e   m/1   =ĥ   1ml ( n ) b   1 ( n )+{circumflex over ( h )} 2ml ( n ) b   2 ( n )− w   ml   H   r   m1 ( n )  
           e   m/2   =ĥ   2ml ( n ) b   1 *( n )−{circumflex over ( h )} 1ml ( n ) b   2 *( n )− w   ml   H   r   m2 ( n )  
         [0074]    where μ is a step size that determines how fast the tap weight of the adaptive filter is changed; and the subscript * represents a complex conjugate operation.  
         [0075]    The error generating blocks  510 ,  520 ,  530 , and  540  shown in FIG. 2 combine the outputs y ml ,ĥ ml  of the adaptive filter blocks  410 ,  420 ,  430 , and  440  and the output b(n) of the reference signal selector  800  to generate two error signals e ml1 ,e ml2 .  
         [0076]    [0076]FIG. 3 illustrates the adaptive filter block shown in FIG. 2, in which an adaptive filter block for the l-th path of the m-th receiver antenna is shown.  
         [0077]    The adaptive filter block comprises, as shown in FIG. 3, an adaptive filter  411 , a simple mean channel estimator  412 , and a multiple constraint criterion filter coefficient updater  413 .  
         [0078]    First, the adaptive filter  411  receives the output signals r m1 (n),r m2 (n) of the input signal generating buffers  310 ,  320 ,  330 , and  340  to generate adaptive filter signals y ml1 ,y ml2  for the received signals of the m-th antenna. Let y m (n) be a vector composed of those signals, then it can be expressed by the following equation 10.  
         [0079]    The filter output comprising the vector y m (n) is acquired by adapting a tap weight to the input signal and can be defined as w m1   H r m1 (n)≡y ml,1 ,w m1   H r m2 (n)≡y ml,2 (i.e., y ml ≡└y ml,1 y ml,2 ┘).  
                       y   m          (   n   )       =              [             w     m                 1     H            r     m                 1            (   n   )                     (       w     m                 1     H            r     m                 2            (   n   )         )     *             ⋮               w     m                 1     H            r     m                 1            (   n   )                     (       w     m                 L     H            r     m                 2            (   n   )         )     *           ]     =       [             y       m                 1     ,   1            (   n   )                   y       m                 1     ,   2     *          (   n   )               ⋮               y       m                 L     ,   1            (   n   )                   y       m                 L     ,   2     *          (   n   )             ]     =           [             h     1      m                 1            (   n   )               h     2      m                 1            (   n   )                   h     2      m                 1     *          (   n   )             -       h     1      m                 1     *          (   n   )                 ⋮       ⋮               h     1                 m                 L            (   n   )               h     2      m                 L            (   n   )                   h     2      m                 L     *          (   n   )             -       h     1      m                 L     *          (   n   )               ]            ≡         H   ~     m          (   n   )                    [             d   1          (   n   )                   d   2          (   n   )             ]            ≡     d        (   n   )               +         u   _     m          (   n   )                         =                    H   ~     m          (   n   )            d        (   n   )         +         u   _     m          (   n   )                       Equation                 10                               
 
         [0080]    where {overscore (u)} m (n) is a noise signal that is obtained by removing the received signals of MAI and inter-path interference, and it is presumed to be AWGN.  
         [0081]    From the above signal model, the maximum likelihood receiver that employs an interference suppression receiver can be readily induced according to the following equation:  
                   d   ^       CMMSE   +   ML            (   n   )       =     arg            min             d   ^                  ∑     m   =   1     M                     y   m          (   n   )       -           H   ^     m          (   n   )       ·     d   ^              2                 Equation                 11                               
 
         [0082]    The cost function in the equation 11 can be divided into two cost functions related to d 1 (n) and d 2 (n), respectively, as expressed by the following equation 12, and used for data detection:  
                         d   ^       1   ,     CMMSE   +   ML              (   n   )       =            arg            min             d   ^                         ∑     m   =   1     M                       ∑     l   =   1     L          [         w     m                 l     H            r     m                 1            (   n   )                h   ^       1      m                 l     *          (   n   )         +         (       w     m                 l     H            r     m                 1            (   n   )         )     *              h   ^       2      m                 l            (   n   )           ]         -     d   ^                          =            arg            min             d   ^                         ∑     m   =   1     M                       ∑     l   =   1     L          [           y       m                 l     ,   1            (   n   )                h   ^       1      m                 l     *          (   n   )         +           y       m                 l     ,   2     *          (   n   )       *              h   ^       2      m                 l            (   n   )           ]         -     d   ^                                d   ^       2   ,     CMMSE   +   ML              (   n   )       =            arg            min             d   ^                         ∑     m   =   1     M                       ∑     l   =   1     L          [         w     m                 l     H            r     m                 1            (   n   )                h   ^       2      m                 l     *          (   n   )         -         (       w     m                 l     H            r     m                 1            (   n   )         )     *              h   ^       1      m                 l            (   n   )           ]         -     d   ^                          =            arg            min             d   ^                         ∑     m   =   1     M                       ∑     l   =   1     L          [           y       m                 l     ,   1            (   n   )                h   ^       2      m                 l     *          (   n   )         -           y       m                 l     ,   2     *          (   n   )       *              h   ^       1      m                 l            (   n   )           ]         -     d   ^                            Equation                 12                               
 
         [0083]    On the other hand, the simple mean channel estimator  412  estimates channel parameters using pilot symbols received from the transmitter.  
         [0084]    In the first embodiment, the following equation 13 is applied to determine channel estimates ĥ 1ml (n) and ĥ 2ml (n) for h 1ml (n) and h 2ml (n). To reduce the effect of noises, preliminary channel estimates obtained from multiple pilot symbols are averaged to give a final channel estimate.  
                         h   ^       1      m                 l            (   n   )       =       1     2        N   P                ∑     i   =   1       N   P                           w     m                 l     H          (     n   -     i                 Q       )            (           d   1   *          (     n   -     i                 Q       )              r     m                 1            (     n   -     i                 Q       )         -         d   2   *          (     n   -     i                 Q       )              r     m                 2            (     n   -     i                 Q       )           )                             h   ^       2      m                 l            (   n   )       =       1     2        N   P                ∑     i   =   1       N   P                           w     m                 l     H          (     n   -     i                 Q       )            (           d   2   *          (     n   -     i                 Q       )              r     m                 1            (     n   -     i                 Q       )         +         d   1   *          (     n   -     i                 Q       )              r     m                 2            (     n   -     i                 Q       )           )                         Equation                 13                               
 
         [0085]    where N P  is the number of consecutive pilot symbols used for channel estimation; and Q is a pilot symbol insertion period.  
         [0086]    The above channel estimate is a maximum likelihood (ML) channel estimate that is unbiased and acquired on the assumption that the channel is estimated with multiple adaptive filters having multiple constraint criteria.  
         [0087]    The multiple constraint criterion filter coefficient updater  413  uses the outputs e ml (n) of the error generating blocks  510 ,  520 ,  530 , and  540  and the outputs y ml,1 (n),y ml,2 (n) of the adaptive filters  411  to update the filter coefficient w ml (n) of the adaptive filter  411  according to the equation 9.  
         [0088]    [0088]FIG. 5 is a detailed block diagram of the coupler  600  shown in FIG. 2.  
         [0089]    The coupler  600  is for calculating the parts  
         ∑     m   =   1     M                       ∑     l   =   1     L            [           y       m                 l     ,   1            (   n   )                h   ^       1      m                 l     *          (   n   )         +           y       m                 l     ,   2     *          (   n   )       *              h   ^       2      m                 l            (   n   )           ]                   and                 ∑     m   =   1     M                       ∑     l   =   1     L            [           y       m                 l     ,   1            (   n   )                h   ^       2      m                 l     *          (   n   )         -           y       m                 l     ,   2     *          (   n   )       *              h   ^       1      m                 l            (   n   )           ]                   of                             
 
         [0090]    the equation 12 using the output y ml (n) of the adaptive filter blocks  410 ,  420 ,  430 , and  440  and the channel estimation result ĥ ml (n), and it comprises, as shown in FIG. 5, L×M energy calculating blocks  610 ,  620 ,  630 , and  640  and two summators  650  and  660 .  
         [0091]    [0091]FIG. 6 is a diagram of the energy calculating blocks  610 ,  620 ,  630 , and  640  shown in FIG. 5, in which an energy calculating block for the l-th path of the m-th receiver antenna is illustrated.  
         [0092]    The energy calculating blocks  610 ,  620 ,  630 , and  640  comprise, as shown in FIG. 6, three complex conjugate calculators  611 ,  612 , and  613 , four multipliers  614 ,  615 ,  616 , and  617 , and two summators  618  and  619 .  
         [0093]    First, the complex conjugate calculators  611 ,  612 , and  613  calculate complex conjugates for the output y ml,2 (n) of the adaptive filter  411  and the output ĥ 2m,l (n)ĥ 1m,l (n) of the simple mean channel estimator  412 .  
         [0094]    The multiplier  614  multiplies the outputs of the complex conjugate calculators  611  and  612  to output y ml,2 *(n)·ĥ 2m,l (n); the multiplier  615  multiplies the output of the adaptive filter  411  by the output of the complex conjugate calculator  612  to output y ml,1 (n)·ĥ 2m,l *(n); the multiplier  616  multiplies the output of the adaptive filter  411  by the output of the conjugate calculator  613  to output y ml,1 (n)·ĥ 1m,l *(n); and the multiplier  617  multiplies the outputs of the complex conjugate calculators  611  and  613  to output y ml,2 *(n)ĥ 1m,l *(n).  
         [0095]    The summator  618  sums the outputs of the multipliers  614  and  616  to output y ml,1 *(n)·ĥ 1m,l *(n)+y ml,2 *(n)·ĥ 2m,l *(n) to the summator  650 ; and the summator  619  sums the outputs of the multipliers  615  and  617  to output y ml,1 (n)·ĥ 2m,l *(n)+y ml,2 *(n)·ĥ 1m,l *(n) to the summator  660 .  
         [0096]    Finally, the summator  650  sums the outputs of the energy calculating blocks  610 ,  620 ,  630 , and  640  for L multipath components and M receiver antennas  40  to output  
           ∑     m   =   1     M                       ∑     l   =   1     L            [           y       m                 l     ,   1            (   n   )                h   ^       1      m                 l     *          (   n   )         +           y       m                 l     ,   2     *          (   n   )       *              h   ^       2      m                 l            (   n   )           ]                   as                       d   _     1          (   n   )             ;                         
 
         [0097]    and the summator  660  sums the outputs of the energy calculating blocks  610 ,  620 ,  630 , and  640  for L multipath components and M receiver antennas  40  to output  
         ∑     m   =   1     M                       ∑     l   =   1     L            [           y       m                 l     ,   1            (   n   )                h   ^       2      m                 l     *          (   n   )         -           y       m                 l     ,   2     *          (   n   )       *              h   ^       1      m                 l            (   n   )           ]                   as                         d   _     2          (   n   )       .                               
 
         [0098]    [0098]FIG. 7 is a detailed block diagram of the data determiner  700  shown in FIG. 2.  
         [0099]    The data determiner  700  calculates a distance between the output {overscore (d)} 1 (n),{overscore (d)} 2 (n) of the coupler  600  and codes on the respective constellation to determine data with a code of the minimum distance, and determines the remaining part of the equation 12. The data determiner  700  comprises, as shown in FIG. 7, a real part calculator  710 , an imaginary part calculator  720 , a distance calculator  730 , and a code selector  740 .  
         [0100]    First, the real part calculator  710  extracts the real part of the complex signal {overscore (d)} 1 (n),{overscore (d)} 2 (n) fed into the data determiner  700  from the coupler  600 , and the imaginary part calculator  720  extracts the imaginary part of the complex signal.  
         [0101]    The distance calculator  730  receives the output signals of the real part calculator  710  and the imaginary part calculator  720  and calculates the distances between the output signals and all the codes on the modified constellation.  
         [0102]    The code selector  740  determines a code having the least deviation from the distance calculated by the distance calculator  730  and outputs the complex value {circumflex over (d)} 1 (n),{circumflex over (d)} 2 (n) of the corresponding code.  
         [0103]    Generally, a training process for adapting filter coefficients is necessary for the purpose of converging the adaptive receiver on the optimal solution. The reference signal selector  800  shown in FIG. 2 selects known data d(n) irrespective of the output data of the data determiner  700  and outputs them as b(n) in the training interval, because the receiver  50  has the information of the transmitted code value in the training interval or the pilot transmission interval. If it is not in the training interval, the reference signal selector  800  selects the output {circumflex over (d)}(n) of the data determiner  700  and generates the selected signal as b(n). In the embodiment of the present invention, additional training data necessary to general adaptive receivers are not needed because pilot symbols regularly inserted are used instead of the training data.  
         [0104]    [0104]FIG. 8 is a diagram of the error generating blocks  510 ,  520 ,  530 , and  540  shown in FIG. 2, in which an error generating block for the l-th path of the m-th receiver antenna is illustrated.  
         [0105]    The error generating blocks  510 ,  520 ,  530 , and  540  comprise, as shown in FIG. 8, two complex conjugate calculators  511  and  512 , four multipliers  513 ,  514 ,  515 , and  516 , and two summators  517  and  518 .  
         [0106]    First, the complex conjugate calculators  511  and  512  calculate complex conjugates for the output {circumflex over (b)} 1 (n),{circumflex over (b)} 2 (n) of the reference signal selector  800 .  
         [0107]    The multipliers  513  and  514  multiply the output of the simple mean channel estimator  412  by the output of the reference signal selector  800  to output ĥ 1m,l ·{circumflex over (b)} 1 (n) and ĥ 2m,l ·{circumflex over (b)} 2 (n), respectively; and the multipliers  515  and  516  multiply the outputs of the complex conjugate calculators  511  and  512  by the output of the simple mean channel estimator  412  to output ĥ 1m,l ·{circumflex over (b)} 2 *(n) and ĥ 2m,l ·{circumflex over (b)} 1 *(n), respectively.  
         [0108]    The summator  517  sums the outputs of the multipliers  513  and  514  and subtracts the output y ml,1 (n) of the adaptive filter  411  from the sum to output ĥ 1m,l ·{circumflex over (b)} 1 (n)+ĥ 2m,l ·{circumflex over (b)} 2 (n)−y ml,1 (n) as an error signal e ml,1 (n) to the multiple constraint criterion filter coefficient updater  413 ; and the summator  518  subtracts the output of the multiplier  515  and the output y ml,2 (n) of the adaptive filter  411  from the output of the multiplier  516  to output ĥ 2m,l ·{circumflex over (b)} 1 *(n)−ĥ 1m,l ·{circumflex over (b)} 2 *(n)−y ml,2 (n) as an error signal e ml,2 (n) to the multiple constraint criterion filter coefficient updater  413 .  
         [0109]    [0109]FIG. 4 is a diagram of the adaptive filter block according to a second embodiment of the present invention, in which an adaptive filter block for the l-th path of the m-th receiver antenna is illustrated. The adaptive filter block comprises an adaptive filter  411 , a maximum likelihood channel estimator  414  and a single constraint criterion filter coefficient updater  415 . The same reference numeral is assigned to the same component having the same function as described in the first embodiment with reference to FIG. 3. Namely, the function of the adaptive filter shown in FIG. 4 is the same as that of the adaptive filter shown in FIG. 3.  
         [0110]    As illustrated in FIG. 4, the second embodiment is the same as the first embodiment, excepting the adaptive filter coefficient updating method and the channel estimation method.  
         [0111]    Let w ml1  be an adaptive filter coefficient for odd received signals for the l-th multipath component of the m-th antenna and w ml2 be an adaptive filter coefficient for even received signals, then the tap weights of the respective adaptive filters are determined so as to minimize the CMMSE cost function of the following equation:  
                   minimize     w     m                 l                 1              E        [                  h     1      m                 l            (   n   )              b   1          (   n   )         +         h     2                 m                 l            (   n   )              b   2          (   n   )         -       w     m                 l                 1     H            r     m                 1            (   n   )                2     ]                     subject                 to                   w     m                 l                 1     H          c   l       =   1          
              minimize     w     m                 l                 2              E        [                  h     2      m                 l            (   n   )              b   l   *          (   n   )         -         h     1                 m                 l            (   n   )              b   2   *          (   n   )         -       w     m                 l                 2     H            r     m                 2            (   n   )                2     ]                     subject                 to                   w     m                 l                 2     H          c   l       =   1             Equation                 14                               
 
         [0112]    where w ml1    0  and w ml2   0  are optimal filter coefficients that minimize a given constrained cost function. w ml1   0  and w ml2    0  are calculated by the Lagrangian multiplier method according to the following equation:  
           w   ml1   0   =R   r     m1     −1   c   l   /c   l   H   R   r     m1     −1   c   l    Equation 15  
         
       w 
       ml2 
       0 
       =R 
       r 
       
         m2 
       
       −1 
       c 
       l 
       /c 
       l 
       H 
       R 
       r 
       
         m2 
       
       −1 
       c 
       l  
     
         [0113]    where R r     m1   =E[r m1 (n)r m1   H (n)] and R r     m2   =E[r m2 (n)r m2   H (n)] can be written to have the same value as the following equation, giving the same optimal solution:  
               R     r   m1       =       R     r   m2       =       R     u   m1       +       ∑     l   =   1     L          E        ⌊         (                h     1      ml            (   n   )            2     +              h     2      ml            (   n   )            2       )     ·     c   l            c   l   H       ⌋                     Equation                 16                               
 
         [0114]    As the two adaptive filters have the same optimal solution as in the equation 16, it is more effective to average two error values and use the mean value in updating the filter coefficient rather than to independently adapt the filters.  
         [0115]    By using the adaptive filter coefficient w ml (n) for the l-th multipath component of the m-th receiver antenna for the l-th multipath component of the m-th receiver antenna in filtering for both odd and even symbols, the error values between the outputs of the adaptive filters from the odd and even received signals and the reference signal are properly combined to update the filter coefficient as the following equation 17. This updating process of the filter coefficient is performed with the single constraint criterion filter coefficient updater  415 .  
           w   ml ( n+ 1)= w   ml ( n )+μ P   c     l     ⊥   [e   ml,1 *( n ) r   ml,1 ( n )+ e   ml,2 *( n ) r   ml,2 ( n )]  Equation 17  
         [0116]    where P c     l     ⊥ =I−c l c l   H /(c l   H c l ) represents an orthogonal projection matrix that projects all the vectors on a space orthogonal to the c l  vector; and e ml1  and e ml2  are as defined in the first embodiment.  
         [0117]    When adapting the adaptive filter  411  according to the equation 17, the inner product of the coefficient of the adaptive filter  411  for the corresponding path and the desired user&#39;s other multipath components is not always zero. When using the simple mean channel estimation method according to the equation 13 used in the first embodiment, the channel estimate is biased as in the following equation 18, causing a deterioration of channel estimation performance and hence a performance deterioration.  
                     E        [         h   ^       1      ml            (   n   )       ]       =       1     2        N   p                ∑     i   =   1       N   p              ∑     l   =   1     L              h     1      mi            (   n   )              w   mi   H          (     n   -   iQ     )            c   l                           E        [         h   ^       2      ml            (   n   )       ]       =       1     2        N   p                ∑     i   =   1       N   p              ∑     l   =   1     L              h     2      mi            (   n   )              w   mi   H          (     n   -   iQ     )            c   l                           Equation                 18                               
 
         [0118]    Accordingly, in the second embodiment, the maximum likelihood channel estimation method according to the following equation 19 is used to obtain unbiased channel estimates, and the channel estimation is performed with the maximum likelihood channel estimator  414 .  
               [             h   ^       1      ml               h   ^       2      ml                   h   ^       1      m2               h   ^       2      m2               ⋮       ⋮               h   ^       1      mL               h   ^       2      mL             ]     =       1     2        N   p                ∑     i   =   1       N   p                  [         1             w   ml   H          (     n   -   iQ     )            c   2           ⋯             w   ml   H          (     n   -   iQ     )            c   L                     w   m2   H          (     n   -   iQ     )            c   1           1       ⋯             w   m2   H          (     n   -   iQ     )            c   L               ⋮       ⋮       ⋰       ⋮                 w   mL   H          (     n   -   iQ     )            c   1                 w   mL   H          (     n   -   iQ     )            c   2           ⋯       1         ]       -   1       ·     W   m   H          A                 Equation                 19                               
 
         [0119]    where  
         W m =[w m1  w m2  . . . w mL ],  
           A=[d   1 *( n−iQ ) r   m1 ( n−iQ )− d   2 *( n−iQ ) r   m2 ( n−iQ )  d   2 *( n−iQ ) r   m1 ( n−iQ )+ d   1 *( n−iQ ) r   m2 ( n−iQ )] 
         [0120]    The adaptive interference suppression receiver for the space-time block coded DS/CDMA system of Ref. 2 using more than two transmitter antennas can be realized with an adequate modification of the first and second embodiments of the present invention. The input signal generating buffers  310  and  320  transmit signals corresponding to the respective symbols of the whole coding block according to the symbol period separately by the signal receiving intervals, and the coupler  600  is properly modified to generate a combined signal suitable for the respective coding method. The individual channel estimators  412  and  414  properly combine the adaptive filter outputs by as many as the number of transmitter antennas and reference signals to perform channel estimation in the same manner as described in the first and second embodiments.  
         [0121]    [0121]FIG. 9 is a graph showing the result of a simulation of the adaptive interference suppression receiver  50  for the space-time block coded DS/CDMA system according to the embodiment of the present invention, in which the result of a comparison (in regard to bit error rate) with the maximum likelihood receiver of Ref. 2 is shown.  
         [0122]    For the simulation, a Rayleigh fading channel with one (refer to (a) of FIG. 9) or three (refer to (b) of FIG. 9) multipaths is used. The simulation involves fifty users having the same transmission power and a 3GPP short scrambling code having a period of 256 chips. The simulation is performed at a fading change rate of 10 −3  that can be expressed by the multiplication f D T of Doppler frequency f D  by the code period Time T, while switching the number of receiver antennas between one and two.  
         [0123]    In (a) and (b) of FIG. 9, A is the performance result of the maximum likelihood receiver (with one receiver antenna) of Ref. 2; B is the performance result of the interference suppression receiver (with one receiver antenna) of the present invention; C is the performance result of the maximum likelihood receiver (with two receiver antennas) of Ref. 2; and D is the performance result of the interference suppression receiver (with two receiver antennas) of the present invention.  
         [0124]    It can be seen from the result of the simulation that the present invention is much superior in data detection performance to the conventional maximum likelihood receiver of Ref. 2, irrespective of the number of multipaths.  
         [0125]    [0125]FIG. 10 is a graph showing the result of a simulation of the adaptive interference suppression receiver  50  for the space-time block coded DS/CDMA system according to the embodiment of the present invention, in which the result of a comparison (in regard to bit error rate) with a receiver expanded from the maximum likelihood receiver of Ref. 1 and the receiver of Ref. 3 in the CDMA environment is shown.  
         [0126]    For the simulation, a single-path fading channel having a fading change rate of 10 −3  that can be expressed by the multiplication f D T of Doppler frequency f D  by the code period Time T is used. In this channel environment, the simulation is performed with two receiver antennas, one desired user and one interfering user having a 100-fold power. The spreading code used in the simulation is a gold code having a period of 63.  
         [0127]    In FIG. 10, E is the performance result of the maximum likelihood receiver of Ref. 1; F is the performance result of a receiver expanded from the receiver of Ref. 3 in the CDMA environment; and G is the performance result of the interference suppression receiver (with two receiver antennas) of the present invention.  
         [0128]    It can be seen from the result of the simulation that the present invention is much superior in data detection performance to the receiver expanded from the maximum likelihood receiver of Ref. 1 and the receiver of Ref. 3 in the CDMA environment.  
         [0129]    While this invention has been described in connection with what is presently considered to be the most practical and preferred embodiment, it is to be understood that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.  
         [0130]    As described above, the present invention, if applied to a CDMA receiver, has a higher data detection performance than the conventional maximum likelihood receiver in the multipath fading channel environment that is an actual mobile communication environment, providing high-speed and high-quality services and increasing the service area of one base station.  
         [0131]    Also, the present invention effectively eliminates MAI without a loss of the reception diversity gain.  
         [0132]    Also, the present invention lowers the transmission power of the terminal, increases the life time of the battery in the terminal, and removes MAI to eliminate a stern power control.  
         [0133]    Furthermore, the present invention only uses information concerning the desired user signals (i.e., spreading code and synchronization information) without a need for signals of all the users, and can be readily applied to the terminal as well as the base station.