# Patent

## Title
Signal detection by a receiver in a multiple antenna time-dispersive system

## Abstract

In a MIMO system the bit error rate floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) equalizer for all of the respective transmit antenna—receive antenna pairings that are possible in the MIMO system. The resulting joint equalization compensates not only for the impact of the channel on the transmit antenna—receive antenna pairings but also for the interference of the other transmit antennas on any given receive antenna. The joint equalization outperforms simply replicating the prior art minimum mean square error (MMSE) equalizer for each transmit antenna—receive antenna pairings.

## Background
TECHNICAL FIELD
This invention relates to the art of wireless communications, and more particularly, to wireless communication systems using multiple antennas at the transmitter and multiple antennas at the receiver, so called multiple-input multiple-output (MIMO) systems.
BACKGROUND OF THE INVENTION
It has previously been assumed for MIMO systems that the time dispersion between one or more of the transmit antennas and one or more of the receive antennas was negligible, i.e., the various paths were nonresolvable, as often occurs in low bandwidth systems, because the pulse width is longer than the channel time dispersion. However, it has been recognized that under typical urban (TU) conditions, i.e., the conditions of the so-called “TU” model, that the time dispersion between one or more of the transmit antennas and one or more of the receive antennas is nonnegligible. Such a non negligible time dispersion causes the various paths to interfere with each other, resulting in a bit error rate floor, i.e., minimum, and so the resulting bit error rate is unacceptable.
I. Ghauri and D. Slock have shown, in “Linear Receivers for the DS-CDMA Downlink Exploiting Orthogonality of Spreading Codes”, 32nd Asilomar Conference, Nov. 1-4, 1998 pp. 650-4, that a minimum mean square error (MMSE) equalizer operating on received code division multiple access (CDMA) chips can be employed to compensate for time dispersion in a single transmit, single receive antenna system, thus reducing the bit error rate floor and improving performance.
SUMMARY OF THE INVENTION
We have recognized, in accordance with the principles of the invention, that in a MIMO system the bit error rate floor caused by time dispersion can be reduced by employing a joint equalizer for all of the respective transmit antenna—receive antenna pairings that are possible in the MIMO system. Advantageously, the resulting joint equalization compensates not only for the impact of the channel on the transmit antenna—receive antenna pairings but also for the interference of the other transmit antennas on any given receive antenna. In a particular embodiment of the invention, the joint equalizer is a joint minimum mean square error (MMSE) equalizer, and in such an embodiment the joint equalization outperforms simply replicating the prior art minimum mean square error (MMSE) equalizer for each transmit antenna—receive antenna pairing.
In one embodiment of the invention, which is especially useful for CDMA, after the equalization is completed the resulting chip streams, one for each transmit antenna, are despread in the conventional manner and then the resulting depread symbols may be further processed in the conventional manner. Alternatively, instead of further processing the despread symbols in the conventional manner, the despread symbols may be processed, in accordance with an aspect of the invention, so as to have their soft bits computed through the use of a posteriori probability (APP) metric. Prior to computing the soft bits, the despread symbols may be spatially whitened using a spatial whitening filter.
In another embodiment of the invention, which is also especially useful for CDMA, the equalizer is iteratively computed so that a symbol from one transmit antenna is determined during each iteration. Initially the received samples are stored in a memory. After a symbol for an antenna is determined, the received samples for each of the remaining antennas are then recomputed by subtracting out the determined symbol from the samples as they existed prior to determining the symbol. Once all the symbols for all of the transmit antennas are determined for a symbol period are determined the operation begins anew with the samples corresponding to the next symbol period.



BRIEF DESCRIPTION OF THE DRAWING
In the drawing:
FIG. 1 shows an embodiment a multiple-input multiple-output (MIMO) wireless system in which the bit error rate floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) equalizer for all of the respective transmit antenna—receive antenna pairings that are possible in the MIMO system, in accordance with the principles of the invention;
FIG. 2 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 1;
FIG. 3 shows in more detail the process by which the weights employed by the joint equalizer of FIG. 1 are determined;
FIG. 4 shows another exemplary process, in flow chart form, for the overall operation of system of FIG. 1;
FIG. 5 shows another embodiment a multiple-input multiple-output (MIMO) wireless system in which the bit error rate floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) equalizer for all of the respective transmit antenna—receive antenna pairings that are possible in the MIMO system, in accordance with the principles of the invention;
FIG. 6 shows a more detailed version of a buffer-subtractor of FIG. 5;
FIG. 7 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 5;
FIG. 8 shows a particular embodiment of the joint equalizer of FIG. 1, in which the equalization is performed in the discrete frequency domain, in accordance with an aspect of the invention; and
FIG. 9 shows a particular embodiment of the joint equalizer of FIG. 1, in which the equalization is calculated in the discrete frequency domain and applied in the time domain.



DETAILED DESCRIPTION
The following merely illustrates the principles of the invention. It will thus be appreciated that those skilled in the art will be able to devise various arrangements which, although not explicitly described or shown herein, embody the principles of the invention and are included within its spirit and scope. Furthermore, all examples and conditional language recited herein are principally intended expressly to be only for pedagogical purposes to aid the reader in understanding the principles of the invention and the concepts contributed by the inventor(s) to furthering the art, and are to be construed as being without limitation to such specifically recited examples and conditions. Moreover, all statements herein reciting principles, aspects, and embodiments of the invention, as well as specific examples thereof, are intended to encompass both structural and functional equivalents thereof. Additionally, it is intended that such equivalents include both currently known equivalents as well as equivalents developed in the future, i.e., any elements developed that perform the same function, regardless of structure.
Thus, for example, it will be appreciated by those skilled in the art that any block diagrams herein represent conceptual views of illustrative circuitry embodying the principles of the invention. Similarly, it will be appreciated that any flow charts, flow diagrams, state transition diagrams, pseudocode, and the like represent various processes which may be substantially represented in computer readable medium and so executed by a computer or processor, whether or not such computer or processor is explicitly shown.
The functions of the various elements shown in the FIGS., including any functional blocks labeled as “processors”, may be provided through the use of dedicated hardware as well as hardware capable of executing software in association with appropriate software. When provided by a processor, the functions may be provided by a single dedicated processor, by a single shared processor, or by a plurality of individual processors, some of which may be shared. Moreover, explicit use of the term “processor” or “controller” should not be construed to refer exclusively to hardware capable of executing software, and may implicitly include, without limitation, digital signal processor (DSP) hardware, network processor, application specific integrated circuit (ASIC), field programmable gate array (FPGA), read-only memory (ROM) for storing software, random access memory (RAM), and non-volatile storage. Other hardware, conventional and/or custom, may also be included. Similarly, any switches shown in the FIGS. are conceptual only. Their function may be carried out through the operation of program logic, through dedicated logic, through the interaction of program control and dedicated logic, or even manually, the particular technique being selectable by the implementor as more specifically understood from the context.
In the claims hereof any element expressed as a means for performing a specified function is intended to encompass any way of performing that function including, for example, a) a combination of circuit elements which performs that function or b) software in any form, including, therefore, firmware, microcode or the like, combined with appropriate circuitry for executing that software to perform the function. The invention as defined by such claims resides in the fact that the functionalities provided by the various recited means are combined and brought together in the manner which the claims call for. Applicant thus regards any means which can provide those functionalities as equivalent as those shown herein.
Software modules, or simply modules which are implied to be software, may be represented herein as any combination of flowchart elements or other elements indicating performance of process steps and/or textual description. Such modules may be executed by hardware which is expressly or implicitly shown.
Unless otherwise explicitly specified herein, the drawings are not drawn to scale.
In the description, identically numbered components within different ones of the FIGS. refer to the same components.
FIG. 1 shows an exemplary embodiment of a multiple-input multiple-output (MIMO) wireless system in which the bit error rate floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) equalizer for all of the respective transmit antenna—receive antenna pairings that are possible in the MIMO system, in accordance with the principles of the invention. Shown in FIG. 1 are a) transmitter 101; b) transmit antennas 103; including transmit antennas 103-1 through 103-M; c) receive antennas 105; including receive antennas 105-1 through 105-N; d) receiver front-end 107; e) joint equalizer 109; f) optional despreaders 111; g) soft bit mapper 112, which may include optional spatial whitening filter 113; and optional a posteriori probability (APP) metric processor 115.
Transmitter 101 is a MIMO transmitter, e.g., one in which an original data stream is divided into substreams and each resulting substream is transmitted as a modulated radio signal via an individual one of transmit antennas 103. The transmitted signals pass to the receiver over a time dispersive channel such that signals from each transmit antenna 103 reach each of receive antennas 105.
Receive antennas 105 convert the radio signals impinging upon them into electrical signals, which are supplied to receiver front-end 107. Receiver front-end 107 operates conventionally to produce a plurality of streams of binary numbers, each stream representing samples of the radio signals received at an associated one of antennas 105. Typically receiver front-end 107 performs radio frequency downconversion, filtering, sampling, and analog-to-digital conversion. The resulting samples are provided to joint equalizer 109.
Joint equalizer 109 compensates for the effects of the transmit signals from each of antennas 103 having passed through the channel as well as the interference that results from transmitting via multiple antennas simultaneously. The output of joint equalizer 109 is M, i.e., the number of transmit antennas, streams of corrected symbols, or in the case of CDMA, streams of corrected chips which when properly combined form symbols. Operation of joint equalizer 109 will be described more fully hereinbelow. If CDMA is employed, the output of joint equalizer 109 is supplied to optional despreaders 111, which produces symbols from the stream of chips supplied by joint equalizer 109.
The symbols produced may then be further processed in the conventional manner for a MIMO system, e.g., soft bits may be developed by soft bit mapper 112 for use in a decoder, e.g., the well known “Turbo decoder”. Alternatively, the symbols may be supplied within soft bit mapper 112 to optional spatial whitening filter 113, which makes the noise equal on each branch. Note that the whitening is performed only in the space domain. If whitening is performed in the time domain some temporal dispersion will be introduced into the signal. The symbols, or whitened symbols if optional spatial whitening filter 113 is employed, may further be supplied to optional a posteriori probability (APP) metric processor 115 within soft bit mapper 112, in accordance with an aspect of the invention. APP metric processor 115 performs a particular type of mapping from symbols to soft bits. Operation of APP metric processor 115 is described more fully hereinbelow.
FIG. 2 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 1. Prior to performing the process of FIG. 2, initial values of parameters M, N, L, P, E, and d must be determined. What each of these parameters represents is listed in Table 1. Additionally, prior to performing the process of FIG. 2 it is necessary to determine noise covariance Rpp, which contains samples of the autocorrelation of the chip pulse shape r(t), when CDMA is employed, or the symbol pulse shape autocorrelation, when CDMA is not employed.







TABLE 1





Parameter Definitions


















M
number of transmit antennas



N
number of receive antennas



L
length of the channel impulse response (chips)



P
over-sampling factor



E
length of the equalizer (chips)



d
equalizer delay (chips)



σn2
power of interference + noise (RF bandwidth)



σx2
power of downlink signal (RF bandwidth)



G
Number of CDMA chips per symbol










The process of FIG. 2 is executed periodically, with a periodicity that preferably does not exceed the coherence time, Tco, which is the time duration for which the channel properties are substantially constant. The process is entered in step 201, in which L*P discrete channel estimations ho to hLP-1 are developed. The channel estimate may be obtained using any conventional technique, e.g., by using correlators tuned to the pilot channel. Each of discrete channel estimations is taken with a time spacing of the chip duration divided by P, for CDMA, or symbol duration divided by P for TDMA. Also, in step 201, the background noise plus interference power σn2 and the power of the downlink from the base station to the terminal σx2 are determined in the conventional manner. In step 203, the weights employed by joint equalizer 109 are determined, as will be described fully further hereinbelow.
Next, in step 205, a set of P samples from each antenna is obtained. Thereafter, in step 207, the determined weights of joint equalizer 109 are applied to the samples from each antenna by joint equalizer 109. The samples are then despread by despreader 111, if CDMA was employed, in optional step 209. In optional step 211, conventional soft mapping of the symbols to soft bits is performed, and the resulting soft bits are supplied as an output for use by a decoder. Thereafter, conditional branch point 213 tests to determine whether one coherence time has elapsed since the previous execution of step 201. If the test result in step 213 is NO, indicating that the channel is believed to still remain substantially the same as when it was last estimated, control passes to step 205, and the process continues as described above. If the test result in step 213 is YES, indicating that sufficient time has passed such that the channel may have changed enough so as not to be considered substantially the same as when it was last estimated, control passes back to step 201 and the process continues as described above.
FIG. 3 shows in more detail the process of step 203 by which the weights employed by joint equalizer 109. Note that the process requires the use of several matrices, and the dimensions of the various matrices are given in Table 2.










TABLE 2







Matrix
Dimension









Rp
EP × EP



Rpp
NEP × NEP



Γ(hn,m)
EP × (E + L − 1)



Γ(H)
NEP × M(E + L − 1)



Γ(Hm)
NEP × M(E + L − 1)



ed
(E + L − 1) × 1



A
M × M(E + L − 1)



am
1 × M(E + L − 1)



W
M × NEP



wm
1 × NEP



Q
M × M










In step 301, the channel estimates ho to hLP-1 for each transmit and receive pair obtained in step 201 are arranged in a respective matrix hn,m as shown in equation 1. In step 303, matrix convolution operator Γ(hn,m) is then formed for each respective matrix hn,m, as shown in equation 2. Thereafter, in step 305, MIMO convolution operator Γ(H) is then formed from the various matrix convolution operators as shown in equation 3.




  
    
      
        
          h
          
            n
            ,
            m
          
        
        =
        
          [
          
            
              
                
                  h
                  0
                
              
              
                
                  h
                  P
                
              
              
                ⋯
              
              
                
                  h
                  
                    
                      (
                      
                        L
                        -
                        1
                      
                      )
                    
                    ⁢
                    P
                  
                
              
            
            
              
                ⋮
              
              
                ⋮
              
              
                
                  
                
              
              
                ⋮
              
            
            
              
                
                  h
                  
                    P
                    -
                    1
                  
                
              
              
                
                  h
                  
                    
                      2
                      ⁢
                      P
                    
                    -
                    1
                  
                
              
              
                ⋯
              
              
                
                  h
                  
                    LP
                    -
                    1
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          1
          )
        
      
    
  
  
    
      
        
          Γ
          ⁡
          
            (
            
              h
              
                n
                ,
                m
              
            
            )
          
        
        =
        
          [
          
            
              
                
                  h
                  
                    n
                    ,
                    m
                  
                
              
              
                
                  0
                  p
                
              
              
                ⋯
              
              
                
                  0
                  p
                
              
            
            
              
                
                  0
                  p
                
              
              
                
                  
                
              
              
                
                  
                
              
              
                
                  
                
              
            
            
              
                ⋮
              
              
                
                  
                
              
              
                ⋰
              
              
                ⋮
              
            
            
              
                
                  0
                  p
                
              
              
                ⋯
              
              
                
                  0
                  p
                
              
              
                
                  h
                  m
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          2
          )
        
      
    
  
  
    
      
        
          Γ
          ⁡
          
            (
            H
            )
          
        
        =
        
          [
          
            
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        1
                        ,
                        1
                      
                    
                    )
                  
                
              
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        1
                        ,
                        2
                      
                    
                    )
                  
                
              
              
                ⋯
              
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        1
                        ,
                        M
                      
                    
                    )
                  
                
              
            
            
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        2
                        ,
                        1
                      
                    
                    )
                  
                
              
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        2
                        ,
                        2
                      
                    
                    )
                  
                
              
              
                
                  
                
              
              
                ⋮
              
            
            
              
                ⋮
              
              
                
                  
                
              
              
                ⋰
              
              
                
                  
                
              
            
            
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        N
                        ,
                        1
                      
                    
                    )
                  
                
              
              
                ⋯
              
              
                
                  
                
              
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        N
                        ,
                        M
                      
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          3
          )
        
      
    
  




In step 307, delay vector ed is formed as shown in equation 4. Delay vector ed is a one dimensional vector with E+L−1 elements which are all zero except for the single value at the E+L−1−d location, which has a value of 1. A typical value for d, which is selectable by the implementor, is such that the location which has a value of 1 is at the center of the vector. The purpose of delay vector ed, is to impose the overall equalizer delay d onto the equalizer. Thereafter, in step 309, delay matrix A is computed from equation 5, in which IM is an identity matrix of size M×M.

ed=[0 . . . 100]  equation (4)

A=IM⊕ed  equation (5)

Finally, in step 311, equalizer weight matrix, W, is computed in accordance with equation 6, in which XH means the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X.




  
    
      
        W
        =
        
          A
          ⁢
          
            
          
          ⁢
          
            
              Γ
              ⁡
              
                (
                H
                )
              
            
            H
          
          ⁢
          
            
              (
              
                
                  
                    
                      Γ
                      ⁡
                      
                        (
                        H
                        )
                      
                    
                    H
                  
                  ⁢
                  
                    Γ
                    ⁡
                    
                      (
                      H
                      )
                    
                  
                
                +
                
                  
                    
                      σ
                      n
                      2
                    
                    
                      σ
                      x
                      2
                    
                  
                  ⁢
                  
                    R
                    pp
                  
                
              
              )
            
            
              -
              1
            
          
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          6
          )
        
      
    
  




In one embodiment of the invention, execution of step 205 is such that the P samples r of antenna n are initially arranged as a vector c shown in equation 7, where k is the current received chip time index if CDMA is employed, or the symbol time index if CDMA is not employed. E consecutive in time vectors c for antenna n are then arranged as shown in equation 8, and the E consecutive in time vectors c for all of the N antennas are further arranged as shown in equation 9, forming a vector of received samples at time k.




  
    
      
        
          
            c
            n
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          [
          
            
              
                
                  
                    r
                    n
                  
                  ⁡
                  
                    (
                    
                      kT
                      c
                    
                    )
                  
                
              
            
            
              
                
                  
                    r
                    n
                  
                  ⁡
                  
                    (
                    
                      
                        kT
                        c
                      
                      +
                      
                        T
                        s
                      
                    
                    )
                  
                
              
            
            
              
                ⋮
              
            
            
              
                
                  
                    r
                    n
                  
                  ⁡
                  
                    (
                    
                      
                        kT
                        c
                      
                      +
                      
                        
                          (
                          
                            P
                            -
                            1
                          
                          )
                        
                        ⁢
                        
                          T
                          s
                        
                      
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          7
          )
        
      
    
  
  
    
      
        
          
            r
            n
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          [
          
            
              
                
                  
                    c
                    n
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
            
              
                
                  
                    c
                    n
                  
                  ⁡
                  
                    (
                    
                      k
                      -
                      1
                    
                    )
                  
                
              
            
            
              
                ⋮
              
            
            
              
                
                  
                    c
                    n
                  
                  ⁡
                  
                    (
                    
                      k
                      -
                      E
                      +
                      1
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          8
          )
        
      
    
  
  
    
      
        
          r
          ⁡
          
            (
            k
            )
          
        
        =
        
          [
          
            
              
                
                  
                    r
                    1
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
            
              
                
                  
                    r
                    2
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
            
              
                ⋮
              
            
            
              
                
                  
                    r
                    N
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          9
          )
        
      
    
  




The application of the determined weights by joint equalizer 109 in step 207 may be performed as shown in equation 10, where y(k) is a resulting vector of size M×1 which contains the equalized chips if CDMA is employed, or symbols if CDMA is not employed.

y(k)=Wr(k)  equation (10)

FIG. 4 shows another exemplary process, in flow chart form, for the overall operation of system of FIG. 1. The version of the process shown in FIG. 4 is similar to that shown in FIG. 2, but it is for embodiments of the invention that include optional spatial whitening filter 113 (FIG. 1) and optional a posteriori probability (APP) metric processor 115. Unless otherwise noted, all variables and parameters employed in the process of FIG. 4 are the same as described for FIG. 2. As with the process of FIG. 2, prior to performing the process of FIG. 4, initial values of parameters M, N, L, P, E, and d must be determined. Additionally, prior to performing the process of FIG. 4 it is necessary to determine Rpp.
The process of FIG. 4 is executed periodically, with a periodicity that preferably does not exceed the coherence time, Tco, which is the time duration for which the channel properties are substantially constant. The process is entered in step 401, in which L*P discrete channel estimations ho to hLP-1 are developed. Also in step 401, the background noise plus interference power σn2 and the power of the downlink from the base station to the terminal σx2 are determined in the conventional manner. In step 403, the weights W employed by joint equalizer 109 are determined, as described hereinabove. Additionally, in step 403, the effective channel matrix Heff and, optionally, spatial whitening filter Q are determined, in accordance with an aspect of the invention. More specifically, Heff is determined as shown in equation 11, and Q is determined as shown in equation 12.

Heff(m,n)=enTWΓ(H)emm=1:M, n=1:N  equation (11)

Q=(WWH)−1/2  equation (12)

Next, in step 405 a set of P samples from each antenna is obtained. Thereafter, in step 407, the determined weights of joint equalizer 109 are applied to the samples from each antenna by joint equalizer 109. The samples are then despread by despreader 111, if CDMA was employed, in optional step 409. In step 411, whitening filter Q is applied to the despread, equalizer outputs of step 409. In step 413, APP softbits are computed, in accordance with an aspect of the invention, using equations 13 and 14. The softbits are the output of the process of FIG. 4, and they may be made available to a decoder, e.g., a turbo decoder.




  
    
      
        
          
            
              
                
                  L
                  D
                
                ⁡
                
                  (
                  
                    
                      x
                      k
                    
                    ❘
                    y
                  
                  )
                
              
              =
              
              ⁢
              
                
                  log
                  e
                
                ⁡
                
                  (
                  
                    
                      Pr
                      ⁡
                      
                        (
                        
                          
                            x
                            k
                          
                          =
                          
                            1
                            ❘
                            y
                          
                        
                        )
                      
                    
                    
                      Pr
                      ⁡
                      
                        (
                        
                          
                            x
                            k
                          
                          =
                          
                            
                              -
                              1
                            
                            ❘
                            y
                          
                        
                        )
                      
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  ln
                  ⁡
                  
                    (
                    
                      
                        
                          ∑
                          
                            
                              x
                              ^
                            
                            ,
                            
                              
                                x
                                k
                              
                              =
                              1
                            
                          
                        
                        ⁢
                        
                          
                            p
                            ⁡
                            
                              (
                              
                                
                                  y
                                  ❘
                                  X
                                
                                =
                                
                                  x
                                  ^
                                
                              
                              )
                            
                          
                          ⁢
                          
                            
                              ∏
                              
                                j
                                ≠
                                k
                              
                            
                            ⁢
                            
                              p
                              ⁡
                              
                                (
                                
                                  X
                                  =
                                  
                                    x
                                    j
                                  
                                
                                )
                              
                            
                          
                        
                      
                      
                        
                          ∑
                          
                            
                              x
                              ^
                            
                            ,
                            
                              
                                x
                                k
                              
                              =
                              
                                -
                                1
                              
                            
                          
                        
                        ⁢
                        
                          
                            p
                            ⁡
                            
                              (
                              
                                
                                  y
                                  ❘
                                  X
                                
                                =
                                
                                  x
                                  ^
                                
                              
                              )
                            
                          
                          ⁢
                          
                            
                              ∏
                              
                                j
                                ≠
                                k
                              
                            
                            ⁢
                            
                              p
                              ⁡
                              
                                (
                                
                                  X
                                  =
                                  
                                    x
                                    j
                                  
                                
                                )
                              
                            
                          
                        
                      
                    
                    )
                  
                
                +
              
            
          
        
        
          
            
              
              ⁢
              
                ln
                ⁡
                
                  (
                  
                    
                      p
                      ⁡
                      
                        (
                        
                          
                            x
                            k
                          
                          =
                          1
                        
                        )
                      
                    
                    
                      p
                      ⁡
                      
                        (
                        
                          
                            x
                            k
                          
                          =
                          
                            -
                            1
                          
                        
                        )
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        13
      
    
  
  
    
      
        
          p
          ⁡
          
            (
            
              
                y
                ❘
                X
              
              =
              
                x
                k
              
            
            )
          
        
        =
        
          
            
              (
              
                1
                
                  
                    
                      2
                      ⁢
                      π
                    
                  
                  ⁢
                  σ
                
              
              )
            
            
              N
              r
            
          
          ⁢
          
            exp
            ⁡
            
              (
              
                -
                
                  
                    
                      
                      
                        y
                        -
                        
                          
                            H
                            eff
                          
                          ⁢
                          
                            s
                            ^
                          
                        
                      
                      
                    
                    2
                  
                  
                    2
                    ⁢
                    
                      σ
                      2
                    
                  
                
              
              )
            
          
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        14
      
    
  




Thereafter, conditional branch point 415 tests to determine whether one coherence time has elapsed since the previous execution of step 401. If the test result in step 415 is NO, indicating that the channel is believed to still remain substantially the same as when it was last estimated, control passes to step 405, and the process continues as described above. If the test result in step 415 is YES, indicating that sufficient time has passed such that the channel may have changed enough so as not to be considered substantially the same as when it was last estimated, control passes back to step 401 and the process continues as described above.
FIG. 5 shows another exemplary embodiment of a multiple-input multiple-output (MIMO) wireless system in which the bit error rate floor caused by time dispersion is reduced by employing a joint minimum mean square error (MMSE) equalizer for all of the respective transmit antenna—receive antenna pairings that are possible in the MIMO system, in accordance with the principles of the invention. Shown in FIG. 5 are a) transmitter 501; b) transmit antennas 503; including transmit antennas 503-1 through 503-M; c) receive antennas 505; including receive antennas 505-1 through 505-N; d) receiver front-end processor 507; e) buffer-subtractor 521; f) joint equalizer 523; g) optional despreader 525; h) soft bit mapper 527; i) space-time regenerator 529; j) order controller 531 and k) switch 533.
Transmitter 501 is a MIMO transmitter, e.g., one in which an original data stream is divided into substreams and each resulting substream is transmitted as a modulated radio signal via an individual one of transmit antennas 503. The transmitted signals pass to the receiver over a time dispersive channel such that signals from each transmit antenna 503 reach each of receive antennas 505.
Receive antennas 505 convert the radio signals impinging upon them into electrical signals, which are supplied to receiver front-end processor 507. Receiver front-end processor 507 operates conventionally to produce a stream of binary numbers representing samples of the radio signals received at antennas 505. Typically receiver front-end processor 507 performs radio frequency downconversion, filtering, sampling, and analog-to-digital conversion. The resulting samples are provided to buffer-substractor 521.
Buffer-subtractor 521 is shown in more detail in FIG. 6. Buffer subtractor 521 includes buffer 601, memory 603, and subtractor 605. Buffer 601 stores a time consecutive set of samples from each of antennas 505 as those samples become available from front-end processor 507. Buffer 601 supplies the samples it receives from front-end processor 507 to memory 603 when an entire set is stored, i.e., the contents of buffer 601 are quickly dumped to memory 603. Once the set of samples are stored in memory 603 they may be independently accessed and sent to joint equalizer 523 (FIG. 5). Additionally, space time samples from a second input of buffer-subtractor 521 are supplied to subtractor 605. Subtractor 605 is cable of forming the difference between a specified location in memory 603 and a second input to buffer-subtractor 521. The resulting difference is stored in the specified location in memory 603.
Joint equalizer 523 performs M passes through memory 603 with a different equalizer weight wm in each pass. Each weight is chosen to emphasize transmit antenna m and to suppress transmit antennas m+1 through M.
The output of joint equalizer 523 is supplied to conventional soft bit mapper 527, via optional conventional despreader 525 if CDMA is employed. The soft bits developed by soft bit mapper 527 are then supplied as an output, such as may be used by a decoder, e.g., the well known “Turbo decoder”.
The same output that is supplied to conventional soft bit mapper 527 may also be supplied to space-time regenerator 529 via switch 533. Space-time regenerator 529 forms a set of time consecutive samples for each receive antenna assuming the soft symbol is correct. In other words, assuming a particular soft symbol had been the actual symbol transmitted by a particular transmit antenna, space-time regenerator 529 creates the corresponding effect that such a symbol would have caused on each of receive antennas 505 given the channel characteristics. Operation of space-time regenerator 529 will be explained more fully hereinbelow.
Order controller 531 determines, based on channel estimates, the signal from which transmit antenna will be processed at any particular time, as will be explained more fully hereinbelow.
FIG. 7 shows an exemplary process, in flow chart form, for the overall operation of system of FIG. 5 in which switch 533 is connected between despreader 525 and space-time regenerator 529. Prior to performing the process of FIG. 7, initial values of parameters M, N, L, P, E, and d must be determined. What each of these parameters represents is listed in Table 1. Additionally, prior to performing the process of FIG. 7 it is necessary to determine Rpp, which is samples of the autocorrelation of the chip pulse shape r(t), when CDMA is employed, or the symbol pulse shape autocorrelation, when CDMA is not employed.
The process of FIG. 7 is executed periodically, with a periodicity that preferably does not exceed the coherence time, Tco. The process is entered in step 701, in which L*P discrete channel estimations ho to hLP-1 are developed. The channel estimate may be obtained using any conventional technique, e.g., by using correlators tuned to the pilot channel. Each discrete channel estimation is taken with a time spacing of the chip duration divided by P, for CDMA, or symbol duration divided by P for TDMA. Also in step 701, the background noise plus interference power σn2 and the power of the downlink from the base station to the terminal σx2 are determined in the conventional manner.
In step 703, order controller 531 determines, according to equations 15 and 16, the order in which the signals from the various transmit antennas will be processed, with the signal from a respective transmit antenna being processed for each execution of joint equalizer 523. Sort is a function that rearranges the elements of vector P so that they run from largest to smallest and order is a list of all the antenna transmit antenna numbers as they should be processed by joint equalizer 523. It is preferable to process the so-called “strong” signals first. However, the particular characteristic, or set of characteristics, which are used to define the “strength” of a signal is at the discretion of the implementer. In the particular embodiment shown herein, estimated signal powers are employed as the strength. Further note that although herein the order for all the antennas is determined simultaneously, those of ordinary skill in the art will readily recognize that it is possible to successively determine which antenna to process.

P=∥diag(WΓ(H))∥2  equation (15)

[P′,order]=sort(P)  equation (16)

In step 705, the equalizer weights for the particular antenna currently being processed, m, as specified by the order, is determined according to equations 17 and 18, in which delay vector am is the mth row of delay matrix A which was described hereinabove.




  
    
      
        
          w
          m
        
        =
        
          
            a
            m
          
          ⁢
          
            
              Γ
              ⁡
              
                (
                
                  H
                  m
                
                )
              
            
            H
          
          ⁢
          
            
              (
              
                
                  
                    
                      Γ
                      ⁡
                      
                        (
                        
                          H
                          m
                        
                        )
                      
                    
                    H
                  
                  ⁢
                  
                    Γ
                    ⁡
                    
                      (
                      
                        H
                        m
                      
                      )
                    
                  
                
                +
                
                  
                    
                      σ
                      n
                      2
                    
                    
                      σ
                      x
                      2
                    
                  
                  ⁢
                  
                    R
                    pp
                  
                
              
              )
            
            
              -
              1
            
          
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          17
          )
        
      
    
  
  
    
      
        
          Γ
          ⁡
          
            (
            
              H
              m
            
            )
          
        
        =
        
          [
          
            
              
                0
              
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        1
                        ,
                        2
                      
                    
                    )
                  
                
              
              
                0
              
              
                
                  ⋯Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        1
                        ,
                        M
                      
                    
                    )
                  
                
              
            
            
              
                0
              
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        2
                        ,
                        2
                      
                    
                    )
                  
                
              
              
                0
              
              
                ⋮
              
            
            
              
                ⋮
              
              
                ⋮
              
              
                ⋮
              
              
                ⋮
              
            
            
              
                0
              
              
                
                  Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        N
                        ,
                        2
                      
                    
                    )
                  
                
              
              
                0
              
              
                
                  ⋯Γ
                  ⁡
                  
                    (
                    
                      h
                      
                        N
                        ,
                        M
                      
                    
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          18
          )
        
      
    
  




Note that for each iteration m of equation 18, columns corresponding to order(1) through order(m−1) are set to the block zero matrix 0. Note that order is vector which contains a listing of the M antenna numbers in the order in which they will be processed. Order may be the result of the well known function sort of MatLab®. Doing so accounts for the fact that signals from transmit antennas 1 through m−1 have already been subtracted from the signal remaining to be processed for this set of samples.
Next, in step 707 a set of samples that span at least the duration of a data symbol is obtained from each receive antenna. Thereafter, in step 709, counter variable m is initialized to 1. In step 711, the weights determined for joint equalizer 523 in step 705 are applied to the samples from each receive antenna by joint equalizer 523, in accordance with equation 19.

y(k)=WmTr(k)  equation (19)

The samples are then despread by despreader 525, if CDMA was employed, in optional step 713. In optional step 715, conventional soft mapping of the symbols to soft bits is performed, and the resulting soft bits are supplied as an output for use by a decoder.
In step 717, samples are produced by space-time regenerator 529 according to equations 20 and 21.




  
    
      
        
          
            
              x
              ^
            
            m
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          
            
              
                d
                ^
              
              m
            
            ⁡
            
              (
              k
              )
            
          
          ⁡
          
            [
            
              
                
                  
                    s
                    ⁡
                    
                      (
                      
                        k
                        ⁢
                        
                          
                        
                        ⁢
                        G
                      
                      )
                    
                  
                
              
              
                
                  
                    s
                    ⁡
                    
                      (
                      
                        
                          k
                          ⁢
                          
                            
                          
                          ⁢
                          G
                        
                        +
                        1
                      
                      )
                    
                  
                
              
              
                
                  ⋮
                
              
              
                
                  
                    s
                    ⁡
                    
                      (
                      
                        
                          k
                          ⁢
                          
                            
                          
                          ⁢
                          G
                        
                        +
                        G
                        -
                        1
                      
                      )
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          20
          )
        
      
    
  
  
    
      
        
          
            y
            m
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          
            [
            
              
                
                  
                    Γ
                    ⁡
                    
                      (
                      
                        h
                        
                          1
                          ,
                          m
                        
                      
                      )
                    
                  
                
              
              
                
                  
                    Γ
                    ⁡
                    
                      (
                      
                        h
                        
                          2
                          ,
                          m
                        
                      
                      )
                    
                  
                
              
              
                
                  ⋮
                
              
              
                
                  
                    Γ
                    ⁡
                    
                      (
                      
                        h
                        
                          N
                          ,
                          m
                        
                      
                      )
                    
                  
                
              
            
            ]
          
          ⁢
          
            
              
                x
                ^
              
              m
            
            ⁡
            
              (
              k
              )
            
          
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          21
          )
        
      
    
  




In step 719 the output of space-time regenerator 529 is subtracted from the contents of memory 603 (FIG. 6), as shown by equation 22.

r(k)=r(k)−ym(k)  equation (22)

Thereafter, conditional branch point 721 tests to determine if m is equal to M. If the test result in step 721 is NO, indicating that not all of the transmit antennas have yet had their signal contribution processed, m is incremented in step 723. Thereafter, control passes back to step 711 and the process continues as described above. If the test result in step 712 is YES, indicating that all of the transmit antennas have had their signal contribution processed, control passes to conditional branch point 725, which tests to determine whether one coherence time has elapsed since the previous execution of step 701. If the test result in step 725 is NO, indicating that the channel is believed to still remain substantially the same as when it was last estimated, control passes to step 707, and the process continues as described above. If the test result in step 725 is YES, indicating that sufficient time has passed such that the channel may have changed enough so as not to be considered substantially the same as when it was last estimated, control passes back to step 701 and the process continues as described above.
In an alternative configuration of FIG. 5, switch 533 is connected to the output of an error correction decoder which is either directly or indirectly connected to soft bit mapping 527. Doing so may improve performance.
FIG. 8 shows a particular embodiment of joint equalizer 109, in which the equalization is performed in the discrete frequency domain, in accordance with an aspect of the invention. Shown in FIG. 8 making up joint equalizer 109 are a) fast Fourier transform (FFT) processors 801, including EFT processors 801-1 through 801-N, where N is the number of receive antennas; b) channel estimators 803, including channel estimators 803-1 through 803-N; c) fast Fourier transform (FFT) processors 805, including FFT processors 805-1 through 805-N*M, where M is the number of transmit antennas; d) MMSE detection per frequency bin processor 807; and e) inverse fast Fourier transform (IFFT) processors 809, including IFFT processors 809-1 through 809-M.
Each of FFT processors 801 receives from front-end 107 a signal of time domain digital samples corresponding to a respective one of receive antennas 105-N and performs the FFT algorithm on a set of consecutive samples to convert the time domain samples to samples in the discrete frequency domain, rn(ω), where ω is a particular discrete frequency. The number of samples F is at the discretion of the implementer based on a tradeoff between the performance and complexity of FFT processors 801. Typically the number of samples is a power of 2, e.g., 128, although the more samples the more accurate the equalization process will be. The possible values for ω are determined as




  
    ω
    =
    
      
        2
        ⁢
        π
        ⁢
        
          
        
        ⁢
        n
      
      
        
          T
          s
        
        ⁢
        F
      
    
  
  ,




where n ranges from 0 to F−1. The resulting discrete frequency samples for each receive antenna are supplied to MMSE detection per frequency bin processor 807.

Each of channel estimators 803 also receives from front-end 107 a signal of time domain digital samples corresponding to a respective one of receive antennas 105-N and performs a channel estimate for the channel between its respective receive antenna and each of the M transmit antennas, thereby producing M channel estimates. Each channel estimate is a series of complex numbers that defines the impulse response of the channel. Fast Fourier transform (FFT) processors 805 each converts a respective channel estimate into the discrete frequency domain representation thereof, and supplies the resulting discrete frequency domain representation of the channel estimates hn,m(ω) to MMSE detection per frequency bin processor 807.
MMSE detection per frequency bin processor 807 performs the equalization in the frequency domain by computing

z(ω)=(H(ω)HH(ω)+σ2I)−1H(ω)Hr(ω)  equation (23)

where




  
    
      
        
          H
          ⁡
          
            (
            ω
            )
          
        
        =
        
          [
          
            
              
                
                  
                    h
                    
                      1
                      ,
                      1
                    
                  
                  ⁡
                  
                    (
                    ω
                    )
                  
                
              
              
                ⋯
              
              
                
                  
                    h
                    
                      1
                      ,
                      M
                    
                  
                  ⁡
                  
                    (
                    ω
                    )
                  
                
              
            
            
              
                ⋮
              
              
                ⋰
              
              
                ⋮
              
            
            
              
                
                  
                    h
                    
                      N
                      ,
                      1
                    
                  
                  ⁡
                  
                    (
                    ω
                    )
                  
                
              
              
                ⋯
              
              
                
                  
                    h
                    
                      N
                      ,
                      M
                    
                  
                  ⁡
                  
                    (
                    ω
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          24
          )
        
      
    
  
  
    
      
        
          r
          ⁡
          
            (
            ω
            )
          
        
        =
        
          [
          
            
              
                
                  
                    r
                    1
                  
                  ⁡
                  
                    (
                    ω
                    )
                  
                
              
            
            
              
                
                  
                    r
                    2
                  
                  ⁡
                  
                    (
                    ω
                    )
                  
                
              
            
            
              
                ⋮
              
            
            
              
                
                  
                    r
                    N
                  
                  ⁡
                  
                    (
                    ω
                    )
                  
                
              
            
          
          ]
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          25
          )
        
      
    
  
  
    
      
        
          
            σ
            2
          
          =
          
            
              σ
              n
              2
            
            
              σ
              x
              2
            
          
        
        ,
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          26
          )
        
      
    
  




σn2 is the background noise plus interference power,
σx2 is the power of the downlink signal from the base station to the terminal, and
I is the identity matrix.
Each of the resulting M components of resulting vector z(ω) are then inverse frequency transformed from the discrete frequency domain into the time domain by inverse fast Fourier transform (IFFT) processors 809. The time domain equalized outputs are then supplied as the output of joint equalizer 109.
FIG. 9 shows a particular embodiment of joint equalizer 109 in which the equalization is calculated in the discrete frequency domain and applied in the time domain, in accordance with an aspect of the invention. Shown in FIG. 9 making up joint equalizer 109 are a) matrix finite impulse response (FIR) filter 901; b) channel estimators 903, including channel estimators 903-1 through 903-N; c) fast Fourier transform (FFT) processors 905, including FFT processors 905-1 through 905-N*M, where M is the number of transmit antennas; d) MMSE tap weight calculator 907; and e) inverse fast Fourier transform (IFFT) processors 909, including IFFT processors 909-1 through 909-N*M
Matrix finite impulse response (FIR) filter 901 continuously receives from front-end 107 a signal of time domain digital samples corresponding to a respective one of receive antennas 105-N. The number of taps of matrix FIR filter 901 is at the discretion of the implementer based on a tradeoff between performance and complexity. Typically the number of samples is a power of 2, e.g., 128.
Each of channel estimators 903 also receives from front-end 107 a signal of time domain digital samples corresponding to a respective one of receive antennas 105-N and performs a channel estimate for the channel between its respective receive antenna and each of the M transmit antennas, thereby producing M channel estimates. Each channel estimate is a series of complex numbers that defines the impulse response of the channel. Fast Fourier transform (FFT) processors 905 each converts a respective channel estimate into the discrete frequency domain, and supplies the resulting discrete frequency domain representation of the channel estimates hn,m(ω) to MMSE tap weight calculator 907. The number of samples employed by FFT processors 905 for each conversion should be the same as the number of taps in matrix FIR filter 901.
MMSE tap weight calculator 907 develops frequency domain representations of the weights necessary to perform the equalization in the time domain by computing

S(ω)=(H(ω)HH(ω)+σ2I)−1H(ω)H  equation (27)

where H(ω), σ2 are defined as explained hereinabove in connection with FIG. 8.

Each of the resulting M components of resulting vector S(ω) are grouped by frequency into frequency vectors, and are then the frequency vectors are inverse frequency transformed from the discrete frequency domain to become filter weights in the time domain by inverse fast Fourier transform (IFFT) processors 909. The weights are then supplied to matrix FIR filter 901 which utilizes them to perform equalization in the time domain on the signals received from front-end 107 as shown by equation 28,




  
    
      
        
          y
          ⁡
          
            (
            k
            )
          
        
        =
        
          
            ∑
            
              j
              =
              0
            
            
              F
              -
              1
            
          
          ⁢
          
            
              S
              j
            
            ⁢
            
              r
              ⁡
              
                (
                
                  k
                  -
                  j
                
                )
              
            
          
        
      
    
    
      
        equation
        ⁢
        
          
        
        ⁢
        
          (
          28
          )
        
      
    
  




where y (k) is the vector output at time k, y having M components—one for each transmit antenna—, Sj is the M×N filter matrix for delay j, which is the inverse Fourier transform of S(ω), r(k)is the vector input signal which is received by matrix FIR filter 901 as defined in equation (9), and F is the number of samples taken for each FFT.
Given the foregoing, those of ordinary skill in the art will readily recognize that other equalizer algorithms which approximate the operation and performance of MMSE, such as least mean square (LMS), recursive least squares (RLS), or minimum intersymbol interference (ISI) subject to an anchor condition, can be employed in a joint manner, e.g., in a space manner, in the implementation of joint equalizer 109.
Those of ordinary skill in the art will readily recognize that the techniques of the instant invention may be employed in systems in which the various transmit antennas are transmitting at different data rates, e.g., using different encoding rates and/or transmit constellations, such as quaternary phase-shift keying (QPSK) or 16-ary quadrature amplitude modulation (16-QAM). In such a situation, if the embodiment of FIG. 5 is employed, order controller 531 is not required because the antennas must always be processed in an a priority determined order.

## Claims

1. A method for compensating for time dispersion in a receiver of a wireless system that has a plurality of transmit antennas and a plurality of a receive antennas, the method comprising the steps of:
receiving samples for each receive antenna;
determining a joint equalizer solution using channel information for at least one pairing of at least one of said transmit antennas and said receive antennas and said received samples of at least two of said receive antennas; and
applying said determined joint equalizer solution to said received samples from at least one of said receive antennas to develop equalized samples;
wherein said step of determining a joint equalizer solution is performed at least partly in the discrete frequency domain.

2. The invention as defined in claim 1 wherein said joint equalizer solution is a joint minimum mean square error (MMSE) solution.

3. The invention as defined in claim 1 further comprising the step of estimating a channel for said at least one pairing of at least one of said transmit antennas and said receive antennas.

4. The invention as defined in claim 1 said step of applying said determined joint equalizer solution is performed in the frequency domain.

5. A method for compensating for time dispersion in a receiver of a wireless system that has a plurality of transmit antennas and a plurality of receive antennas, the method comprising the steps of:
receiving samples for each receive antenna;
determining a joint equalizer solution using channel information for at least one pairing of at least one of said transmit antennas and said receive antennas and said received samples of at least two of said receive antennas;
applying said determined joint equalizer solution to said received samples from at least one of said receive antennas to develop equalized samples; and
despreading said equalized samples.

6. The invention as defined in claim 1 wherein at least two of said transmit antennas transmit at different rates.

7. The invention as defined in claim 1 wherein at least two of said transmit antennas transmit using different transmit constellations.

8. The invention as defined in claim 1 further comprising the step of performing soft mapping using a version of said equalized samples.

9. The invention as defined in claim 8 wherein said version of said equalized samples are despread samples.

10. The invention as defined in claim 8 wherein said step of performing soft mapping further comprises the step of spatial whitening said version of said equalized samples.

11. The invention as defined in claim 8 herein said step of performing soft mapping further comprises the step of performing a posteriori probability (APP) metric processing an said version of said equalized samples.

12. The invention as defined in claim 1 wherein said determining step is performed multiple times, once for each one of said transmit antennas.

13. The invention as defined in claim 1 wherein said determining and applying steps are iterated multiple times over a symbol period, and iteration for each one of said transmit antennas, and said method further comprises, for each iteration, the steps of:
generating a representation of signals that would have arrived had a particular symbol for a currently being processed transmit antenna had been transmitted; and
subtracting said representation from said samples received for each receive antenna.

14. The invention as defined in claim 1 wherein said joint equalizer solution is one from the group consisting of: a joint least mean square (LMS) solution, a joint recursive least squares (RLS) solution, or a joint minimum intersymbol interference (ISI) subject to an anchor condition solution.

15. A receiver for use in a multiple-input multiple-output (MIMO) system in which a plurality of signal detectors receive signals transmitted by a plurality of signal sources, said receiver comprising:
a joint equalizer that develops a joint equalizer solution using channel information for at least one pairing of at least one of said signal sources and said signal detectors and received samples of at least two of said signal detectors and supplies as an output a signal that includes at least said joint equalizer solution applied to a signal received by at least one of said signal detectors;
a soft bit mapper for developing soft bits from said joint equalizer output; and
a despreader interposed between said joint equalizer and said soft bit mapper.

16. The invention as defined in claim 15 wherein said joint equalizer solution is a joint minimum mean square error (MMSE) equalizer solution.

17. The invention as defined in claim 15 wherein said soft bit mapper further comprises an a posteriori probability (APP) metric processor.

18. The invention as defined in claim 15 wherein said soft bit mapper further comprises a spatial whitening unit.

19. The invention as defined in claim 15 where at least two of said transmit sources transmits signals at different rates.

20. The invention as defined in claim 15 wherein at least two of said transmit sources transmits signals using different transmit constellations.

21. The invention as defined in claim 15 further comprising:
a space time regenerator coupled to said joint equalizer; and
a buffer-subtractor coupled between said signal detectors and said joint equalizer and between said space time regenerator and said joint equalizer.

22. The invention as defined in claim 21 further comprising a front end processor coupled to said buffer-subtractor.

23. The invention as defined in claim 1 further comprising
an error correction decoder coupled to said soft bit mapper;
a space time regenerator coupled to said error correction decoder; and
a buffer-subtractor coupled between said signal detectors and said joint equalizer and between said space time regenerator and said joint equalizer.

24. The invention as defined in claim 23 further comprising a front end processor coupled to said buffer-subtractor.

25. The invention as defined in claim 15 further comprising an order controller for determining an order in which signals from said signal detectors will be processed by said joint equalizer.

26. A receiver for use in a multiple-input multiple-output (MIMO) system in which a plurality of signal detectors receive signals transmitted by a plurality of signal sources, said receiver comprising:
a joint equalizer that develops a joint equalizer solution using channel information for at least one pairing of said at least one of said signal sources and said signal detectors and received samples of at least two of said signal detectors and supplies as an output a signal that includes at least said equalizer solution applied to a signal received by at least one of said signal detectors; and
a soft bit mapper for developing soft bits from said joint equalizer output;
wherein said joint equalizer solution is a joint minimum mean square error (MMSE) equalizer solution and wherein said join equalizer further comprises:
a first plurality of fast Fourier transform processors, each of said fast Fourier transform processors being coupled to receive a respective input corresponding to a signal received by one of said signal detectors and supplying as an output a discrete frequency domain representation thereof;
a plurality of channel estimation units each of which is coupled to receive a is respective input corresponding to a signal received by one of said signal detectors which develops a channel estimate for each channel between each respective signal source and each respective signal detector;
a second plurality of fast Fourier transform processors, each of said second plurality of fast Fourier transform processors coupled to receive a respective input corresponding to a channel estimate for a respective one of said channels between said signal sources and said signal detectors and supplying as an output a discrete frequency domain representation thereof;
an MMSE detection per frequency bin processor coupled to receive as inputs said outputs from said first plurality of fast Fourier transform processors and from said second plurality of fast Fourier transform processors to produce a discrete frequency domain representation of an application of said joint minimum mean square error (MMSE) equalizer solution to said signals received by each of said signal detectors; and
a plurality of inverse fast Fourier transform processors which convert said discrete frequency domain representation of an application of said joint minimum mean square error (MMSE) equalizer solution to the time domain.

27. The invention as defined in claim 26 wherein M is the number of signal sources and N is the number of signal detectors, and
wherein said MMSE detection per frequency bin processor performs the equalization in the frequency domain by computing


z(ω)=(H(ω)HH(ω)+σ2I)−1H(ω)Hr(ω)


where




  
    
      H
      ⁢
      
        (
        ω
        )
      
    
    =
    
      [
      
        
          
            
              
                h
                
                  1
                  ,
                  1
                
              
              ⁡
              
                (
                ω
                )
              
            
          
          
            ⋯
          
          
            
              
                h
                
                  1
                  ,
                  M
                
              
              ⁡
              
                (
                ω
                )
              
            
          
        
        
          
            ⋮
          
          
            ⋰
          
          
            ⋮
          
        
        
          
            
              
                h
                
                  N
                  ,
                  1
                
              
              ⁡
              
                (
                ω
                )
              
            
          
          
            ⋯
          
          
            
              
                h
                
                  N
                  ,
                  M
                
              
              ⁡
              
                (
                ω
                )
              
            
          
        
      
      ]
    
  




  
    
      r
      ⁢
      
        (
        ω
        )
      
    
    =
    
      [
      
        
          
            
              
                r
                1
              
              ⁡
              
                (
                ω
                )
              
            
          
        
        
          
            
              
                r
                2
              
              ⁡
              
                (
                ω
                )
              
            
          
        
        
          
            ⋮
          
        
        
          
            
              
                r
                N
              
              ⁡
              
                (
                ω
                )
              
            
          
        
      
      ]
    
  




  
    
      σ
      2
    
    =
    
      
        σ
        n
        2
      
      
        σ
        x
        2
      
    
  


σ2 is the background noise plus interference power,
σx2 is the sum of the power received by all said signal detectors from all of said signal sources,
each r(ω) is said output of a one of said first plurality of first Fourier transform processors,
each h(ω) is said output of a one of said second plurality of first Fourier transform processors,
I is the identity matrix,
XH means the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X, and
z(ω) is the equalized output vector in the discrete frequency domain for frequency ω.

28. A receiver for use in a multiple-input multiple-output (MIMO) system in which a plurality of signal detectors receive signals transmitted by a plurality of signal sources, said receiver comprising:
a joint a equalizer that develops a joint equalizer solution using channel information for at least one pairing of said at least one of said signal sources and said signal detectors and received samples of at least two of said signal detectors and supplies as an output a signal that includes at least said equalizer solution applied to a signal received by at least one of said signal detectors; and
a soft bit mapper for developing soft bits from said joint equalizer output;
wherein said joint equalizer solution is a joint minimum mean square error (MMSE) equalizer solution and wherein said joint equalizer further comprises:
a plurality of channel estimation units each of which coupled to receive a respective input corresponding to a signal received by one of said signal detectors which develops a channel estimate for each channel between each respective signal source and each respective signal detector;
a plurality of fast Fourier transform processors, each of said plurality of fast Fourier transform processors being coupled to receive a respective input corresponding to a channel estimate for a respective one of said channels between said signal sources and said signal detectors and supplying as an output a discrete frequency domain representation thereof;
an MMSE tap weight calculator coupled to receive as inputs said outputs from said plurality of fast Fourier transform processors to produce a discrete frequency domain representation of a joint minimum mean square error (MMSE) equalizer solution to said signals received by each of said signal detectors;
a plurality of inverse fast Fourier transform processors which convert said discrete frequency domain representation of a joint minimum mean square error (MMSE) equalizer solution to matrices of filter coefficients in the time domain; and
a matrix finite impulse response (FIR) fitter coupled to apply said matrices of filter coefficients in the time domain to said signals received by said signal detectors.

29. The invention as defined in claim 28 wherein M is the number of signal sources and N is the number of signal detectors, and wherein said MMSE equalizer solution is developed by computing:


S(ω)=(H(ω)HH(ω)+σ2I)−1H(ω)H

where




  
    
      H
      ⁢
      
        (
        ω
        )
      
    
    =
    
      [
      
        
          
            
              
                h
                
                  1
                  ,
                  1
                
              
              ⁡
              
                (
                ω
                )
              
            
          
          
            ⋯
          
          
            
              
                h
                
                  1
                  ,
                  M
                
              
              ⁡
              
                (
                ω
                )
              
            
          
        
        
          
            ⋮
          
          
            ⋰
          
          
            ⋮
          
        
        
          
            
              
                h
                
                  N
                  ,
                  1
                
              
              ⁡
              
                (
                ω
                )
              
            
          
          
            ⋯
          
          
            
              
                h
                
                  N
                  ,
                  M
                
              
              ⁡
              
                (
                ω
                )
              
            
          
        
      
      ]
    
  




  
    
      σ
      2
    
    =
    
      
        σ
        n
        2
      
      
        σ
        x
        2
      
    
  


σn2 is the background noise plus interference power,
σx2 is the sum of the power received by all said signal detectors from all of said signal sources,
each h(ω) is said output of a one of said plurality of fast Fourier transform processors,
I is the identity matrix,
XH means the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X, and
S(ω) is a frequency domain matrix representing the equalization for frequency ω.

30. The invention as defined in claim 29 wherein said matrix FIR filter applies said matrices of filter coefficients in the time domain to said signals received by said signal detectors by computing




  
    y
    ⁢
    
      (
      k
      )
    
  
  =
  
    
      ∑
      
        j
        =
        0
      
      
        F
        -
        1
      
    
    ⁢
    
      
        S
        j
      
      ⁢
      
        r
        ⁢
        
          (
          
            k
            -
            j
          
          )
        
      
    
  



where y(k) is a vector output at time k, y having one component for each of said signal sources,
Sj is a M×N filter matrix for delay j, which is the inverse Fourier transform of S(ω),
r(k) is a vector of received samples at time k, and
F is the number of samples taken for each fast Fourier transform by each of said plurality of fast Fourier transform processors.

31. The invention as defined in claim 15 wherein said MIMO system is a wireless system, said signal sources are transmit antennas and said detectors are antennas of said receiver.

32. The invention as defined in claim 15 wherein said joint equalizer develops said joint equalizer solution as a function of estimates of the channels between each of said signal sources and said signal detectors.

33. The invention as defined in claim 15 wherein said joint equalizer develops and applies said joint equalizer solution in a time domain.

34. The invention as defined in claim 15 wherein said joint equalizer develops and applies said joint equalizer solution in a frequency domain.

35. The invention as defined in claim 16 wherein said joint equalizer develops and applies said joint equalizer solution in a discrete frequency domain and applies said joint minimum mean square error (MMSE) equalizer solution in a time domain.

36. The invention as defined in claim 16 wherein said joint equalizer develops said joint minimum mean square error (MMSE) equalizer solution by computing




  W
  =
  
    A
    ⁢
    
      
    
    ⁢
    
      
        (
        
          Γ
          ⁢
          
            (
            H
            )
          
        
        )
      
      H
    
    ⁢
    
      
        (
        
          
            
              
                (
                
                  Γ
                  ⁢
                  
                    (
                    H
                    )
                  
                
                )
              
              H
            
            ⁢
            
              Γ
              ⁢
              
                (
                H
                )
              
            
          
          +
          
            
              
                σ
                n
                2
              
              
                σ
                x
                2
              
            
            ⁢
            
              R
              pp
            
          
        
        )
      
      
        -
        1
      
    
  



where
Γ(H) is a MIMO convolution operator,
XH means the Hermitian transpose of X, which is the complex conjugate transpose of the vector or matrix X,
A is a delay matrix,
σn2 is the background noise plus interference power,
σx2 is the sum of the power received by all said signal detectors from all of said signal sources,
noise covariance Rpp,
and said joint equalizer applies said joint minimum mean square error (MMSE) equalizer solution by computing


y(k)=Wr(k)


where r(k) is a vector of received samples at time k, and
is y(k) is the vector output at time k representing the application of said joint minimum mean square error (MMSE) equalizer to said vector of received samples at time k.

37. The invention as defined in claim 15 wherein said joint equalizer solution is one from the group consisting of: a joint least mean square (LMS) solution, a joint recursive least squares (RLS) solution, or a joint minimum intersymbol interference (ISI) subject to an anchor condition solution.

38. A receiver for use in a multiple-input multiple-output (MIMO) system in which a plurality of receive antennas receive signals transmitted by a plurality of transmit antennas, said receiver comprising:
means for (i) developing a joint equalizer solution using channel information for at least one pairing of at least one of said transmit antennas and said receive antennas and said received samples of at least two of said receive antennas, said joint equalizer solution being developed at least partly in a frequency domain, and (ii) supplying as an output a signal that includes at least said joint equalizer solution applied to signal received by at least one of said receive antennas; and
means for developing soft bits from said output signal.