# Patent

## Title
Reduced complexity parametric covariance estimation for precoded MIMO transmissions

## Abstract

A model-based technique for estimating impairment covariance associated with a MIMO signal is disclosed. In an exemplary method, an impairment model is constructed for a received composite information signal comprising at least a first data stream transmitted from first and second antennas according to a first antenna weighting vector. The impairment model includes first and second model terms corresponding to the first and second antennas, respectively, but in several embodiments does not include a cross-antenna interference term. In another embodiment, an impairment model for a received MIMO signal is constructed by computing an impairment model term for each antenna and an additional term to account for preceding interference in a single-stream MIMO transmission scenario. The impairment terms are grouped so that only two associated scaling terms are unknown; values for the scaling terms are estimated by fitting the model to measured impairment covariance values.

## Background
FIELD OF THE INVENTION
The present invention relates generally to wireless communication systems, and in particular to a parametric system and method for covariance estimation in a wireless communication system employing precoded multiple-input multiple-output (MIMO) transmissions.
BACKGROUND
Spread-spectrum communication systems are well known in the art and widely deployed. A class of receivers well suited for use in spread-spectrum systems—such as those standardized in IS-95, IS-2000 (cdma2000), and the 3rd-Generation Partnership Project's (3GPP) Wideband Code-Division Multiple Access (W-CDMA) specifications—is the linear interference-whitening (LIW) receiver. LIW receivers suppress interference in addition to collecting signal energy for detection. One form of the LIW receiver is a transversal chip equalizer; another is a G-Rake receiver. The Rake receiver derives its name from its rake-like structure, wherein multiple receiver “fingers” are used to receive multiple signal images in a received multipath signal. By coherently combining the finger outputs in a weighted Rake combiner, the conventional Rake receiver can use multipath reception to improve the Signal to Interference-plus-Noise Ratio (SINR) of the received signal. A Generalized Rake (G-Rake) receiver improves interference suppression performance over a conventional Rake receiver using more sophisticated generation of the combining weights.
Recently, 2×2 Multiple-Input Multiple-Output (MIMO) technology has been standardized in Release 7 of the 3GPP specifications. The standardized scheme, referred to as Dual-Transmit Adaptive Arrays (D-TxAA), is similar to selective per-antenna rate control (S-PARC), except that adaptive unitary precoding is applied to each of the data streams, in this case to each of one or two High-Speed Downlink Shared Channel (HS-DSCH) data streams.
D-TxAA can be viewed as an extension of the previously standardized closed loop mode 1 (CL-1) transmit diversity scheme, in that the precoding vectors (which map a data stream to the multiple transmit antennas) used for each of the D-TxAA data streams are selected from the same codebook used for CL-1. In contrast to CL-1, however, D-TxAA includes two modes of operation—single-stream mode and dual-stream mode. In single-stream mode, one of the four possible precoding vectors from the CL-1 codebook is applied to a single data stream. In dual-stream mode, orthogonal pairs of precoding vectors (again selected from the CL-1 codebook) are applied to the two data streams. The use of precoding has a significant impact on the receiver, and in particular complicates the design of LIW receivers such as Rake receivers.
Earlier versions of the 3GPP W-CDMA specifications (i.e., prior to Release 7) define two transmit diversity modes, CL-1, and an open-loop mode known as STTD. U.S. patent application Ser. No. 10/800,167 (Pub. No. US 2005/0201447), titled “Method and Apparatus for Parameter Estimation in a Generalized Rake Receiver,” filed Mar. 12, 2004 by Cairns et al. (the “Cairns application”), assigned to the assignee of the present application and incorporated herein by reference in its entirety, discloses a solution for G-Rake receivers in a transmit diversity system. The solution describes a parametric approach to estimating an impairment covariance matrix used to form G-Rake combining weights. The parametric approach estimates the impairment covariance matrix as a sum of terms, including a separate term for each transmit antenna as well as a term corresponding to the sum of noise plus other-cell interference.
This solution works well for open-loop transmit diversity modes. In an open-loop mode, the impairments corresponding to each transmit antenna during a particular symbol period are uncorrelated, since different symbols are transmitted from the different antennas. In closed-loop mode, however, the mobile terminal specifies a phase offset, and the same symbol is transmitted by a primary antenna and simultaneously by a secondary antenna with the specified phase offset. In this case, the impairment due to each transmit antenna is highly correlated. This correlation may be exploited to improve interference suppression and receiver performance. U.S. patent application Ser. No. 11/751,109, titled “Receiver Parametric Covariance Estimation for Transmit Diversity,” filed May 21, 2007 by Jonsson et al. (the “Jonsson application”), assigned to the assignee of the present application and incorporated herein by reference in its entirety, discloses a parametric approach to estimating an impairment covariance matrix that accounts for the simultaneous transmission of the same symbols from a first and second antenna. In this approach the impairment covariance matrix for a system employing two transmit antennas is formulated as a sum of seven terms, including a term corresponding to each of the transmit antennas, a noise-plus-other-cell-interference term, plus four additional terms corresponding to the four possible precoding vectors in the CL-1 codebook. The terms are weighted by fitting parameters determined by fitting the parametrically modeled impairment covariance matrix to a measured impairment covariance matrix. An implicit assumption is that if one or more of the preceding vectors are not used by any user in the cell, then the corresponding fitting parameter will ideally be estimated as zero.
The CL-1 covariance estimation approach described in the Jonsson application applies to the transmission of only a single data stream, mapped according to a precoding vector to two (or more) antennas. In contrast, in D-TxAA, two data streams may be transmitted simultaneously, with both data streams sharing the same set of channelization codes. This creates additional self-interference, referred to as code-reuse interference, which affects the formulation of the impairment covariance. Code reuse is not accounted for in the formulation of the Jonsson application, since only one data stream is ever transmitted in CL-1.
In a co-pending patent application titled “Receiver Parametric Covariance Estimation for Precoded MIMO Transmission,” U.S. patent application Ser. No. 12/036,323 (the “Grant application”), the entire contents of which are incorporated by reference herein, a MIMO G-Rake receiver operating at the symbol level is disclosed that is based upon the most general G-Rake formulation for MIMO. For a 2×2 MIMO scenario, this receiver computes an impairment covariance matrix according to:

Ru=α1R11+α2R22+α12+R12++jα12−R12−βRn,  (1)

where R11 captures the interference due to a first transmit antenna, R22 captures the interference due to a second transmit antenna, R12+ and R12− represent cross-antenna interference, and Rn accounts for white noise passing through the receive filter. The weighting terms are given by:





  
    
      
        
          
            
              
                α
                1
              
              =
              
              ⁢
              
                
                  1
                  
                    
                      
                        γ
                        p
                      
                      ⁡
                      
                        (
                        1
                        )
                      
                    
                    ⁢
                    
                      N
                      p
                    
                  
                
                [
                
                  
                    
                      Γ
                      
                        D
                        /
                        P
                      
                    
                    ⁡
                    
                      (
                      
                        
                          
                            
                              γ
                              s
                            
                            ⁡
                            
                              (
                              1
                              )
                            
                          
                          ⁢
                          
                            
                              
                              
                                b
                                11
                              
                              
                            
                            2
                          
                        
                        +
                        
                          
                            
                              γ
                              s
                            
                            ⁡
                            
                              (
                              2
                              )
                            
                          
                          ⁢
                          
                            
                              
                              
                                b
                                12
                              
                              
                            
                            2
                          
                        
                      
                      )
                    
                  
                  +
                
              
            
          
        
        
          
            
              
              ⁢
              
                
                  
                    Γ
                    
                      V
                      /
                      P
                    
                  
                  ⁢
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      
                        K
                        v
                      
                    
                    ⁢
                    
                      
                        
                          γ
                          v
                        
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                      ⁢
                      
                        
                          
                          
                            v
                            
                              1
                              ⁢
                              i
                            
                          
                          
                        
                        2
                      
                    
                  
                
                +
                
                  
                    
                      γ
                      o
                    
                    ⁡
                    
                      (
                      1
                      )
                    
                  
                  ⁢
                  
                    Γ
                    
                      O
                      /
                      P
                    
                  
                
                +
                
                  
                    γ
                    p
                  
                  ⁡
                  
                    (
                    1
                    )
                  
                
              
              ]
            
          
        
        
          
            
              
                α
                2
              
              =
              
              ⁢
              
                
                  1
                  
                    
                      
                        γ
                        p
                      
                      ⁡
                      
                        (
                        2
                        )
                      
                    
                    ⁢
                    
                      N
                      p
                    
                  
                
                [
                
                  
                    
                      Γ
                      
                        D
                        /
                        P
                      
                    
                    ⁡
                    
                      (
                      
                        
                          
                            
                              γ
                              s
                            
                            ⁡
                            
                              (
                              1
                              )
                            
                          
                          ⁢
                          
                            
                              
                              
                                b
                                21
                              
                              
                            
                            2
                          
                        
                        +
                        
                          
                            
                              γ
                              s
                            
                            ⁡
                            
                              (
                              2
                              )
                            
                          
                          ⁢
                          
                            
                              
                              
                                b
                                22
                              
                              
                            
                            2
                          
                        
                      
                      )
                    
                  
                  +
                
              
            
          
        
        
          
            
              
              ⁢
              
                
                  
                    Γ
                    
                      V
                      /
                      P
                    
                  
                  ⁢
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      
                        K
                        v
                      
                    
                    ⁢
                    
                      
                        
                          γ
                          v
                        
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                      ⁢
                      
                        
                          
                          
                            v
                            
                              2
                              ⁢
                              i
                            
                          
                          
                        
                        2
                      
                    
                  
                
                +
                
                  
                    
                      γ
                      o
                    
                    ⁡
                    
                      (
                      2
                      )
                    
                  
                  ⁢
                  
                    Γ
                    
                      O
                      /
                      P
                    
                  
                
                +
                
                  
                    γ
                    p
                  
                  ⁡
                  
                    (
                    2
                    )
                  
                
              
              ]
            
          
        
        
          
            
              
                α
                12
                +
              
              =
              
              ⁢
              
                
                  1
                  
                    
                      N
                      p
                    
                    ⁢
                    
                      
                        
                          
                            γ
                            p
                          
                          ⁡
                          
                            (
                            1
                            )
                          
                        
                        ⁢
                        
                          
                            γ
                            p
                          
                          ⁡
                          
                            (
                            2
                            )
                          
                        
                      
                    
                  
                
                [
                
                  
                    Γ
                    
                      D
                      /
                      P
                    
                  
                  (
                  
                    
                      
                        
                          γ
                          s
                        
                        ⁡
                        
                          (
                          1
                          )
                        
                      
                      ⁢
                      
                        Re
                        [
                        
                          
                            b
                            11
                          
                          ⁢
                          
                            b
                            21
                            *
                          
                        
                        ]
                      
                    
                    +
                  
                
              
            
          
        
        
          
            
              
                
                  
                  ⁢
                  
                    
                      
                        γ
                        s
                      
                      ⁡
                      
                        (
                        2
                        )
                      
                    
                    ⁢
                    
                      Re
                      [
                      
                        
                          b
                          12
                        
                        ⁢
                        
                          b
                          22
                          *
                        
                      
                      ]
                    
                  
                  )
                
                +
                
                  
                    Γ
                    
                      V
                      /
                      P
                    
                  
                  ⁢
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      
                        K
                        v
                      
                    
                    ⁢
                    
                      
                        
                          γ
                          v
                        
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                      ⁢
                      
                        Re
                        ⁡
                        
                          [
                          
                            
                              v
                              
                                1
                                ⁢
                                i
                              
                            
                            ⁢
                            
                              v
                              
                                2
                                ⁢
                                i
                              
                              *
                            
                          
                          ]
                        
                      
                    
                  
                
              
              ]
            
          
        
        
          
            
              
                α
                12
                -
              
              =
              
              ⁢
              
                
                  1
                  
                    
                      N
                      p
                    
                    ⁢
                    
                      
                        
                          
                            γ
                            p
                          
                          ⁡
                          
                            (
                            1
                            )
                          
                        
                        ⁢
                        
                          
                            γ
                            p
                          
                          ⁡
                          
                            (
                            2
                            )
                          
                        
                      
                    
                  
                
                [
                
                  
                    Γ
                    
                      D
                      /
                      P
                    
                  
                  (
                  
                    
                      
                        
                          γ
                          s
                        
                        ⁡
                        
                          (
                          1
                          )
                        
                      
                      ⁢
                      
                        Im
                        [
                        
                          
                            b
                            11
                          
                          ⁢
                          
                            b
                            21
                            *
                          
                        
                        ]
                      
                    
                    +
                  
                
              
            
          
        
        
          
            
              
                
                  
                  ⁢
                  
                    
                      
                        γ
                        s
                      
                      ⁡
                      
                        (
                        2
                        )
                      
                    
                    ⁢
                    
                      Im
                      [
                      
                        
                          b
                          12
                        
                        ⁢
                        
                          b
                          22
                          *
                        
                      
                      ]
                    
                  
                  )
                
                +
                
                  
                    Γ
                    
                      V
                      /
                      P
                    
                  
                  ⁢
                  
                    
                      ∑
                      
                        i
                        =
                        1
                      
                      
                        K
                        v
                      
                    
                    ⁢
                    
                      
                        
                          γ
                          v
                        
                        ⁡
                        
                          (
                          i
                          )
                        
                      
                      ⁢
                      
                        Im
                        ⁡
                        
                          [
                          
                            
                              v
                              
                                1
                                ⁢
                                i
                              
                            
                            ⁢
                            
                              v
                              
                                2
                                ⁢
                                i
                              
                              *
                            
                          
                          ]
                        
                      
                    
                  
                
              
              ]
            
          
        
        
          
            
              β
              =
              
              ⁢
              
                N
                0
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




Here, Np is the pilot code spreading factor, γx (k) is the fraction of the total base station chip energy allocated to voice (x=v), data (x=s), overhead (x=o), or pilots (x=p) for antenna/stream k, Γx/P is the ratio of chip energies Ex/Ep, bij is the i,j-th element of pre-coding matrix B, and νij is the i-th element of the pre-coding vector v for the j-th voice user. Note that in WCDMA, the columns of B and the pre-coding vectors v are drawn from the codebook Φ={φ1,φ2,φ3,φ4}, where φk=[1 ej(2k-1)π/4 ]T for k=1, 2, 3, 4. The columns of B are chosen from orthogonal pairs of Φ in dual stream mode, whereas only one pre-coding vector is chosen for the first column of B in single stream mode, while the second column is set to the zero vector ([0 0]T).
The G-Rake receiver described in the Grant application utilizes the impairment covariance matrix and net channel estimates to compute combining weights. The combining weights for this receiver structure depend on whether one or two streams are being transmitted. For single stream mode, the combining weights wsingle are obtained by solving the following system of equations:

Ruwsingle=h(b),  (3)

where the notation h(b) indicates the “effective” net channel coefficients that depend on the pre-coding vector b. (b is the first column of B as described above for single stream mode.)

For dual stream mode, two sets of combining weights (w1dual,w2dual) must be computed. These weights may be obtained by solving the following systems of equations




  
    
      
        
          
            
              
                
                  
                    (
                    
                      
                        R
                        u
                      
                      +
                      
                        
                          
                            α
                            PC
                          
                          ⁡
                          
                            (
                            1
                            )
                          
                        
                        ⁢
                        
                          h
                          ⁡
                          
                            (
                            
                              b
                              2
                            
                            )
                          
                        
                        ⁢
                        
                          
                            h
                            H
                          
                          ⁡
                          
                            (
                            
                              b
                              2
                            
                            )
                          
                        
                      
                    
                    )
                  
                  ⁢
                  
                    w
                    1
                    dual
                  
                
                =
                
                  h
                  ⁡
                  
                    (
                    
                      b
                      1
                    
                    )
                  
                
              
            
          
          
            
              
                
                  
                    (
                    
                      
                        R
                        u
                      
                      +
                      
                        
                          
                            α
                            PC
                          
                          ⁡
                          
                            (
                            2
                            )
                          
                        
                        ⁢
                        
                          h
                          ⁡
                          
                            (
                            
                              b
                              1
                            
                            )
                          
                        
                        ⁢
                        
                          
                            h
                            H
                          
                          ⁡
                          
                            (
                            
                              b
                              1
                            
                            )
                          
                        
                      
                    
                    )
                  
                  ⁢
                  
                    w
                    2
                    dual
                  
                
                =
                
                  h
                  ⁡
                  
                    (
                    
                      b
                      2
                    
                    )
                  
                
              
            
          
        
        .
      
    
    
      
        (
        4
        )
      
    
  




Here, h(bn) is the effective net channel coefficient vector due to pre-coding for stream n, and αPC (n) is a per-code scaling factor that multiplies the outer product of the effective net coefficients to account for the code reuse interference (note: pre-coding vector b1 corresponds to the first column of matrix B while pre-coding vector b2 corresponds to the second column).

Symbol estimates for either single or dual stream mode are obtained by computing the inner product of the combining weights for the given stream with the despread traffic symbols.
As noted above, the MIMO G-Rake receiver formulation disclosed in the Grant application is the most general formulation. However, this solution is quite complex. Four fundamental matrix terms must be computed: R11, R22, R12, and Rn. (The matrix term R12 is used to compute R12+ and R12− for Equation (1)). Of these fundamental terms, R11, R22, and Rn are conjugate symmetric, so only slightly more than half of the matrix elements must be calculated. R12, on the other hand, is not conjugate symmetric, so all matrix elements must be computed. Those skilled in the art will appreciate that the calculations required to compute these matrix terms represent a considerable computation burden.
In addition, the formulation of Equation (1) requires that five parameters be estimated to form the impairment covariance matrix. This is compared to the estimation of only two parameters in a “baseline” non-MIMO G-Rake receiver. Even in these simpler receivers, post-estimation smoothing or other adjustment (e.g. clipping) of parameters has been found useful to obtain good overall receiver performance. Joint estimation of five parameters is likely to require similar (and likely even more complex) post-processing to yield good receiver performance.
SUMMARY
A reduced-complexity model-based technique for estimating impairment covariance associated with a MIMO signal is disclosed. In an exemplary method, an impairment model is constructed for a received composite information signal comprising at least a first data stream transmitted from first and second antennas according to a first antenna weighting vector. The impairment model includes first and second model terms corresponding to the first and second antennas, respectively, but in several embodiments omits cross-antenna interference terms. In some embodiments, the first and second model terms are each computed as functions of propagation channel estimates corresponding to the first and second antennas, respectively, and a third impairment model term is computed to account for noise and inter-cell interference. In a variation of this embodiment, the first, second, and third impairment model terms are grouped so that only two scaling parameters associated with the grouped first, second, and third impairment model terms are unknown, and values for the two unknown scaling parameters are estimated by fitting the impairment model to measured impairment covariance values.
In another embodiment, an impairment model for a received MIMO signal is constructed by computing an impairment model term for each antenna and an additional term to account for preceding interference in a single-stream MIMO transmission scenario. The impairment terms are grouped so that only two associated scaling terms are unknown; values for the scaling terms are estimated by fitting the model to measured impairment covariance values.
Another embodiment relates to a wireless communication receiver for use in a mobile terminal in a transmit diversity wireless communication system. The receiver includes a radio front-end circuit configured to provide a received signal of interest containing at least a first data stream transmitted simultaneously from a first and second antenna according to a first antenna weighting vector. The receiver circuit is configured to carry out one or more of the methods described herein for estimating impairment covariance.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a functional block diagram of a wireless communication system.
FIG. 2 is a functional block diagram of an exemplary symbol-level LIW receiver configured to process received signals transmitted from at least first and second transmit antennas.
FIG. 3 is a functional block diagram of an exemplary chip-level LIW receiver configured to process received signals transmitted from at least first and second transmit antennas.
FIG. 4 is a block diagram illustrating the composition of a transmitted signal.
FIG. 5 is a flow diagram of an exemplary method of estimating impairment covariance of a received signal of interest.
FIG. 6 is a flow diagram of an exemplary method for estimating channel quality, e.g. SINR.



DETAILED DESCRIPTION
FIG. 1 depicts an exemplary wireless communication system 100 employing closed-loop transmit diversity, such as CL-1, and/or multiple-input multiple-output (MIMO) transmissions, such as according to the D-TxAA specifications. Within a Radio Access Network (RAN) 102, a Radio Network Controller (RNC) 104 controls a plurality of base transceiver stations (BTS) 106, also known in the art as Node B's. Each Node B 106 provides radio communication services with subscriber mobile terminals 112 within a geographic area called a cell, which may be divided into sectors, as depicted in FIG. 1. The RNC 104 communicates with a Core Network (CN) 114, which in turn is connected to one or more external networks 116, such as the Public Switched Telephone Network (PSTN), the Internet, or the like.
Embodiments of the present invention are described herein with respect to WCDMA standards, including the specifications for D-TxAA, which is more fully described below. However, the invention is not so limited, and the inventive concepts disclosed and claimed herein may be advantageously applied to a wide array of transmit diversity systems.
Each base station 106 includes at least a primary transmit antenna 108 and a secondary transmit antenna 110 (either per-cell or per-sector, depending on the network 100 configuration), as shown in FIG. 2. The base station 106 may transmit an information signal, such as a precoded voice signal or a precoded High-Speed Downlink Packet Access (HSDPA) data signal using both antennas 108, 110. The signal transmitted on the secondary antenna 110 is weighted relative to the signal transmitted on the primary antenna 108, wherein the transmit weights may comprise phase offset only, or may more generally comprise a complex quantity having both phase and amplitude. The phase shift employed may be determined by feedback from the mobile terminal 112, thus forming a closed-loop transmit diversity system.
There are two main types of LIW receiver architectures. One uses symbol-level equalization, which is typically based on maximum likelihood estimation techniques. This type of receiver includes the G-Rake receiver 200, illustrated in block diagram form in FIG. 2. A radio processor 202 generates chip samples from a received signal, which includes the information signal transmitted from antennas 108 and 110 at base station 106. The chip samples are provided to a finger placement circuit 204, which determines the “finger delays,” usually including multipath delays, used to despread a received CDMA signal in a correlation unit 206. The finger delays are also provided to a weight computer 208 that computes combining weights which are used to combine the despread values in a combiner 210 to produce soft values, or estimates of the symbol values.
Another type of LIW receiver is a chip-level equalizer, which typically includes Minimum Mean-Square Error (MMSE) transversal chip equalization 300, as illustrated in block diagram form in FIG. 3. A radio processor 302 generates chip samples from a received signal. The chip samples are provided to a tap placement circuit 304, which determines the tap delays, related to multipath delays, for a Finite Impulse Response (FIR) filter 306. The selected tap delays are also provided to a weight calculator 308 that computes filter coefficients (or weights) for the FIR filter 306. The FIR filter 306 filters the chip samples to produce a signal that is despread by a correlator 310 to produce symbol estimates.
Both types of LIW receivers 200, 300 rely on an estimate of a covariance matrix. In the case of maximum likelihood G-Rake processing, the covariance matrix is an impairment covariance matrix. In the case of MMSE-based processing, a data covariance matrix, which is closely related to the impairment covariance matrix, is used. According to embodiments described herein, an impairment model is constructed, the impairment model generally including several model terms scaled by corresponding scaling parameters. As will be shown below, the model terms may be constructed so that each is independent of antenna weighting vectors (preceding vectors) used for transmitting the information signals. Effects of the antenna weighting vectors are captured by the scaling parameters. The scaling parameters may be calculated, in some embodiments, or may be jointly estimated (“fitted”) according to well-known techniques by fitting the impairment model to measured impairment covariance or data covariance matrices. For convenience, this process is described with respect to a G-Rake receiver and an impairment covariance matrix. However, the use of a data covariance matrix in the case of a chip equalizer is directly analogous.
In order to provide context for a detailed discussion of impairment models, additional background of the D-TxAA MIMO scheme standardized by 3GPP in Release-7 High-Speed Packet Access (HSPA) specifications is provided here. A high-level view of the scheme is shown in FIG. 4. In Release-6 (non-MIMO) HSPA, the theoretical maximum downlink data rate is 14.4 Mbps which is achieved using 15 codes, 16-QAM, and coding rate 1. With D-TxAA, this peak rate may doubled to 28.8 Mbps by transmitting a second, separately encoded data stream in parallel, when channel conditions warrant. Typically, this occurs at high signal-to-noise ratios (SNRs) and when the channel is full rank. In rank-deficient scenarios and/or lower SNRs, the second data stream may be switched off, as indicated in the figure, so that only a single HSPA data stream is transmitted.
D-TxAA uses a form of unitary precoding applied to the HS-DSCH data streams using the preceding weight matrix B. The weight vectors (columns of B) applied to each stream are drawn from the same codebook of four phase-only weights used for the closed-loop mode-1 (CL-1) transmit diversity option defined in Rel-99:




  
    
      
        
          
            u
            i
          
          =
          
            
              1
              
                2
              
            
            ⁡
            
              [
              
                
                  
                    1
                  
                
                
                  
                    
                      e
                      
                        j
                        ⁢
                        
                          π
                          4
                        
                        ⁢
                        
                          (
                          
                            
                              2
                              i
                            
                            -
                            1
                          
                          )
                        
                      
                    
                  
                
              
              ]
            
          
        
        ,
        
          i
          =
          
            1
            ,
            
              
            
            ⁢
            2
            ,
            
              
            
            ⁢
            3
            ,
            
              
            
            ⁢
            4.
          
        
      
    
    
      
        (
        5
        )
      
    
  




Those skilled in the art will appreciate that while phase-only weights are used in D-TxAA as currently standardized, MIMO transmission more generally may use amplitude weighting as well. Further, the codebook from which the weight vectors are drawn is by no means limited to four entries.

In any event, for the case of single-stream D-TxAA transmission, the weight vector is a single one of the four possible vectors in the CL-1 codebook. For the case of dual-stream transmission, the two used weight vectors are selected to be orthogonal; hence, B is unitary. With respect to the CL-1 codebook as defined in Equation (1), the orthogonal pairings are thus (1,3) and (2,4) and permutations thereof.
With this structure, there are fundamentally eight different transmitter modes corresponding to the different choices for B and the different modes—4 single-stream modes and 4 dual-stream modes. The single-stream modes are defined by

B ε{[u1 0],[u2 0],[u3 0],[u4 0]}  (6)

and the dual-stream modes by

Bε{[u1 u3],[u2 u4],[u3 u1],[u4 u2]}.  (7)

A selection of the “best” transmitter mode is made by the user equipment (UE) through maximizing some metric such as data rate, SINR, received power per-stream, or the like. The UE indicates its preferred preceding configuration through feedback, called preceding control information (PCI), transmitted to the Node B on the high-speed dedicated physical control channel (HS-DPCCH), as shown in FIG. 4. The PCI feedback 405 informs the scheduler in the Node B of the preferred number of data streams and the preferred preceding weight vector(s). The Node B signals the actually used preceding matrix B, as well as the actually transmitted number of streams, on the downlink high-speed shared control channel (HS-SCCH). This information is used to configure the receiver in the UE.
As can be seen in FIG. 4, a number of different signals form the composite transmitted signal on each antenna. These include the following: the one or two high-speed downlink shared channel (HS-DSCH) data streams 410 which are precoded, using matrix B, at block 420; a number of dedicated channel signals (e.g., voice and/or control), including voice signals 430, which may be configured using CL-1 transmit diversity using the matrix V at block 440; and a number of other overhead signals 450 (e.g., pilots, voice, control, etc.) that are not precoded. The preceding matrix V is given by

V=[v1v2 . . . vKv],  (8)

where Kv is the number of dedicated channels configured in CL-1 transmit diversity. Each precoding vector (column of V) is given by one of the 4 possible vectors in the CL-1 codebook.

It is important to note that when the transmitter is configured in dual-stream mode, the same set of channelization (spreading) codes is used for both streams. This creates self-interference that must be resolved by the receiver. In other words, when the receiver demodulates each stream, it must suppress interference from the other stream.
Various embodiments of the present invention exploit a new way of expressing the covariance term associated with a given precoded signal. In particular, the MIMO G-Rake receiver described in the Grant application may be simplified, reducing the number of terms required to estimate the impairment covariance matrix, and simplifying estimation of scaling parameters applied to the impairment covariance matrix terms to produce an impairment covariance matrix estimate. The presently disclosed techniques thus significantly reduce complexity compared to previous solutions.
If it is assumed that the overhead signals are uniformly distributed across antennas and that the voice signals are not pre-coded (i.e. vi=[1 0]T), it can be shown that Equation (1) may be simplified considerably. In fact, these assumptions have been shown to have very little impact on receiver performance under a range of conditions. Because the simplifications are mode-specific, the two D-TxAA modes are considered separately below.
First, 3GPP's WCDMA Release 7 specifications specify that, for dual stream MIMO, the data signal power is evenly split between the two streams, and that the pre-coding vectors b1 and b2 are orthogonal. Using these facts and the assumptions above, the impairment covariance matrix for two stream mode may be given by:




  
    
      
        
          R
          u
        
        =
        
          
            
              (
              
                
                  Γ
                  
                    O
                    /
                    P
                  
                
                +
                
                  Γ
                  
                    D
                    /
                    P
                  
                
              
              )
            
            ⁢
            
              (
              
                
                  
                    1
                    
                      2
                      ⁢
                      
                        N
                        p
                      
                      ⁢
                      
                        
                          γ
                          p
                        
                        ⁡
                        
                          (
                          1
                          )
                        
                      
                    
                  
                  ⁢
                  
                    R
                    11
                  
                
                +
                
                  
                    1
                    
                      2
                      ⁢
                      
                        N
                        p
                      
                      ⁢
                      
                        
                          γ
                          p
                        
                        ⁡
                        
                          (
                          2
                          )
                        
                      
                    
                  
                  ⁢
                  
                    R
                    22
                  
                
              
              )
            
          
          +
          
            
              1
              
                N
                p
              
            
            ⁢
            
              (
              
                
                  R
                  11
                
                +
                
                  R
                  22
                
              
              )
            
          
          +
          
            β
            ⁢
            
              
            
            ⁢
            
              
                R
                n
              
              .
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




Those skilled in the art will immediately appreciate that the formulation of Equation (9) omits any cross-antenna terms. Furthermore, the model terms R11, R22, and Rn are independent of the precoding vectors b1 and b2. The functional forms for the terms R11, R22, and Rn are given in the Cairns and Jonsson applications discussed above.

Although the impairment terms R11, R22, and Rn can be grouped in various ways, the grouping of Equation (9) is particularly advantageous, as Equation (9) can also be written as:

Ru=μ(k1R11+k2R22)+k3(R11+R22)+βRn,  (10)

where kn indicates a constant value that is known a priori by the receiver, and μ and β are unknown quantities that the receiver must estimate. Those skilled in the art will appreciate that three conjugate symmetric matrix terms must be computed, while only two parameters must be estimated. This represents a considerable computation savings over the more general formulation described in the Grant application. Those skilled in the art will also appreciate that performance advantages may also be anticipated, since only two parameters must be estimated.

For single-stream mode, the impairment covariance matrix is given by:




  
    
      
        
          R
          u
        
        =
        
          
            
              Γ
              
                O
                /
                P
              
            
            ⁡
            
              (
              
                
                  
                    1
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        N
                        p
                      
                      ⁢
                      
                        
                          γ
                          p
                        
                        ⁡
                        
                          (
                          1
                          )
                        
                      
                    
                  
                  ⁢
                  
                    R
                    11
                  
                
                +
                
                  
                    1
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        N
                        p
                      
                      ⁢
                      
                        
                          γ
                          p
                        
                        ⁡
                        
                          (
                          2
                          )
                        
                      
                    
                  
                  ⁢
                  
                    R
                    22
                  
                
              
              )
            
          
          +
          
            
              1
              
                N
                p
              
            
            ⁢
            
              (
              
                
                  R
                  11
                
                +
                
                  R
                  22
                
              
              )
            
          
          +
          
            
              Γ
              
                D
                /
                P
              
            
            ⁡
            
              (
              
                
                  1
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      N
                      p
                    
                  
                
                ⁢
                
                  
                    R
                    eff
                  
                  ⁡
                  
                    (
                    b
                    )
                  
                
              
              )
            
          
          +
          
            β
            ⁢
            
              
            
            ⁢
            
              R
              n
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




This can also be written as:

Ru=μ1(k1R11+k2R22)+k3(R11+R22)+μ2(k4Reff(b))+βRn,  (12)

where kn again indicates a constant value that is known a priori by the receiver, and μ1, μ2, and β are unknown quantities that the receiver must estimate. The term Reff (b) represents the interference due to pre-coding. This term is calculated in a similar fashion to R11 and R12, except that the medium coefficients used for the computation are “effective” medium coefficients given by





  
    
      
        
          
            g
            eff
          
          =
          
            
              ∑
              
                m
                =
                1
              
              N
            
            ⁢
            
              
                
                  
                    
                      γ
                      p
                    
                    ⁡
                    
                      (
                      1
                      )
                    
                  
                  
                    
                      γ
                      p
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
              
              ⁢
              
                b
                m
              
              ⁢
              
                g
                m
              
            
          
        
        ,
      
    
    
      
        (
        13
        )
      
    
  




where gm is a composite medium coefficient vector corresponding to the mth transmit antenna and bm is the mth element of the pre-coding vector b. Note that gm contains the stacked contributions of the medium coefficients corresponding to the propagation channel between the mth transmit antenna and the lth receive antenna (gm=[gm1T gm2T . . . gmLT]T).

Advantages of this approach and the grouping of terms above are that only four conjugate symmetric matrix terms (R11, R22, Reff(b), and Rn) must be computed for single-stream mode, and only three parameters must be estimated. This formulation also represents a considerable computational savings, and the overall receiver performance is better since only three parameters must be estimated. In addition, the single-stream formulation of Equation (12) can be simplified further by dropping the Reff(b) term, (making it identical to the dual-stream formulation), yielding additional complexity reduction while sacrificing little in terms of overall receiver performance.
FIG. 5 illustrates an exemplary method for estimating and applying impairment covariance according to one or more embodiments of the present invention. Method 500 can be viewed as a direct implementation of equations (9) and (11).
Method 500 begins at block 510, with the obtaining of MIMO information, which may include a MIMO mode (e.g., single-stream or dual-stream) and pre-coding vector(s) (b or B, depending on the mode), the allocation of pilot channel power across antennas (γp(1) and γp(2)), the allocation of data power across streams (γd(1) and γd(2) ), and the number of spreading codes used (K). Those skilled in the art will appreciate that Equations (9) and (11) are based on an assumption that data signal power is evenly split between the streams. Thus, in some embodiments, this even split (e.g., γd(1)=γd(1)=0.5) may simply be assumed. In others, a receiver may selectively employ method 500, based on the known or estimated values for γd(1) and γd(2).
At block 520, net and medium channel coefficients are obtained for each transmit antenna. Techniques for estimating these coefficients, including techniques based on pilot channel measurements, are well known in the art. In one embodiment, the measured net channel response hm is obtained from measurements of de-spread pilot symbols by:




  
    
      
        
          
            h
            m
          
          =
          
            
              1
              
                N
                P
              
            
            ⁢
            
              
                ∑
                
                  k
                  =
                  0
                
                
                  
                    N
                    P
                  
                  -
                  1
                
              
              ⁢
              
                
                  
                    x
                    
                      pilot
                      ,
                      m
                    
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
                ⁢
                
                  
                    s
                    m
                    *
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        14
        )
      
    
  




where Xpilot,m(k) is the vector of de-spread pilot symbols for the k-th symbol corresponding to the m-th transmit antenna in the current slot, sm*(k) is the complex conjugate of the known value for the k-th pilot symbol corresponding to the m-th transmit antenna, and NP is the number of symbols used in the measurement. Once hm has been derived from pilot symbol measurements, corresponding estimates for the medium channel response may be calculated, according to:

gm=(AHA)−1AHhm,  (15)

where A is an M×L (corresponding to M receiver fingers and L multipath delays; M≧L) conversion matrix that incorporates information characterizing the transmitter and receiver contributions to the composite, or net, channel response. The (i,j) element of A is thus given by:

αi,j=RTX/RX(di−τj).  (16)

Here, RTX/RX (λ) is the convolution of the transmit and receive filters evaluated at λ, di is the delay of the i-th receiver finger, and τj is the j-th path delay.

At block 530, the impairment model terms R11, R22, and Rn are computed, based on the medium channel response values gm. The functional form of R11 and R22 was given in the Cairns and Jonsson applications discussed above, and depends on the pulse shape autocorrelation function as well as the medium response channel estimates. Elements in these matrices correspond to pairs of fingers. For example, for finger f1 (associated with delay df1 and receive antenna l1) and for finger f2 (associated with delay df2 and receive antenna l2), the corresponding matrix element of Rmm (m is equal to 1 or 2) is given by




  
    
      
        
          
            r
            
              m
              ,
              m
            
          
          ⁡
          
            (
            
              
                f
                1
              
              ,
              
                f
                2
              
            
            )
          
        
        =
        
          
            ∑
            
              
                P
                1
              
              =
              0
            
            
              P
              -
              1
            
          
          ⁢
          
            
              ∑
              
                
                  P
                  2
                
                =
                0
              
              
                P
                -
                1
              
            
            ⁢
            
              
                
                  g
                  m
                
                ⁡
                
                  (
                  
                    
                      p
                      1
                    
                    ,
                    
                      l
                      1
                    
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    (
                    
                      
                        g
                        m
                      
                      ⁡
                      
                        (
                        
                          
                            p
                            2
                          
                          ,
                          
                            l
                            2
                          
                        
                        )
                      
                    
                    )
                  
                  *
                
                ·
                
                   
                  
                    
                  
                  ⁢
                  
                     
                    
                      
                        
                          ∑
                          
                            
                              k
                              =
                              
                                -
                                ∞
                              
                            
                            
                              k
                              ≠
                              0
                            
                          
                          ∞
                        
                        ⁢
                        
                          
                            x
                            ⁡
                            
                              (
                              
                                
                                  d
                                  
                                    f
                                    1
                                  
                                
                                -
                                
                                  τ
                                  
                                    p
                                    1
                                  
                                
                                -
                                
                                  kT
                                  c
                                
                              
                              )
                            
                          
                          ⁢
                          
                            
                              x
                              *
                            
                            ⁡
                            
                              (
                              
                                
                                  d
                                  
                                    f
                                    2
                                  
                                
                                -
                                
                                  τ
                                  
                                    p
                                    2
                                  
                                
                                -
                                
                                  kT
                                  c
                                
                              
                              )
                            
                          
                        
                      
                      ,
                    
                  
                
              
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




where P is the number of paths, gm (p,l) is the medium response channel coefficient that is associated with receive antenna l and path delay τp, x(τ) is the chip pulse shape autocorrelation function, and Tc is the chip period.

At block 540, impairment covariance or data covariance is measured, to provide an “instantaneous” estimate of the signal impairments. Those skilled in the art will appreciate that whether impairment covariance or data covariance is measured will depend on the receiver implementation. Typically, a G-Rake receiver is configured to perform impairment covariance measurements using pilot channel data. Measured data covariance can be obtained by summing outer products of several sample data vectors (vector elements corresponding to finger outputs), i.e.:




  
    
      
        
          
            
              R
              ~
            
            measured
          
          =
          
            
              ∑
              n
            
            ⁢
            
              
                y
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  y
                  H
                
                ⁡
                
                  (
                  n
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        18
        )
      
    
  




where n is a time index and y(n) is a vector of de-spread pilot values corresponding to time index n.

If the current transmission mode is dual-stream, as determined at block 550, then the three impairment model terms R11, R22, and Rn provide an adequate basis for modeling the signal impairments. In this case, then, the procedure continues immediately at block 560 with estimating values for only two scaling values, μ and β, in accordance with Equation (10). In some embodiments, the scaling values μ and β are “fitted” so that the fitted estimates {circumflex over (μ)} and {circumflex over (β)} satisfy the expression:

{tilde over (R)}measured−k3(R11+R22)={circumflex over (μ)}(k1R11+k2R22)+{circumflex over (β)}Rn.  (19)

(A corresponding fitting expression may be readily derived by those skilled in the art for fitting the impairment terms to measured data covariance, where the fitting equation includes an additional term corresponding to the outer product of the net channel response). As can easily be seen by comparing Equations (9) and (10),





  
    
      k
      1
    
    =
    
      1
      
        2
        ⁢
        
          
        
        ⁢
        
          N
          p
        
        ⁢
        
          
            γ
            p
          
          ⁡
          
            (
            1
            )
          
        
      
    
  
  ,
  
    
      k
      2
    
    =
    
      1
      
        2
        ⁢
        
          
        
        ⁢
        
          N
          p
        
        ⁢
        
          
            γ
            p
          
          ⁡
          
            (
            2
            )
          
        
      
    
  
  ,
  
    
      and
      ⁢
      
        
      
      ⁢
      
        k
        3
      
    
    =
    
      1
      
        N
        p
      
    
  
  ,




each of which can be easily calculated from known values. A least-squares or other approach may be used to fit the weighted sum of terms in Equation (19) to the measured covariance, where equations with complex quantities are treated as two real-value equations, as described in the Cairns application. Clipping of values at zero may be used, as well as other clipping procedures. The instantaneous values of the fitting parameters {circumflex over (μ)} and {circumflex over (β)} may additionally be smoothed, or filtered, to reduce the effects of transient channel conditions.

If, on the other hand, the current transmission is single-stream, then an additional impairment model term may be computed to account for interference caused by precoding. Thus, at block 570, “effective” medium coefficients geff are computed, from the medium channel coefficients obtained at block 520 and the preceding vector b currently in use, according to Equation (13). The effective medium coefficients are used at block 574 to compute an impairment model term Reff (b), which has the same form as R11 and R22, but is instead a function of geff rather than g1 or g2. Then, at block 578, values for three scaling terms, μ1, μ2, and β, are estimated, in accordance with Equation (12). Accordingly, in some embodiments, fitting parameter values {circumflex over (μ)}1, {circumflex over (μ)}2, and {circumflex over (β)} are jointly estimated so that:

{tilde over (R)}measured−k3(R11+R22)≈{circumflex over (μ)}1(k1R11+k2R22)+{circumflex over (μ)}2(k4Reff(b))+{circumflex over (β)}Rn.  (20)

As was the case with the dual-stream expression, k1, k2, and k3 are based on known quantities. In addition, k4 is known, i.e.





  
    
      k
      3
    
    =
    
      1
      
        2
        ⁢
        
          N
          p
        
      
    
  
  ,




but can in any event be accounted for in the estimation of {circumflex over (μ)}2. The values for {circumflex over (μ)}1, {circumflex over (μ)}2, and {circumflex over (β)} may be estimated using known techniques, e.g., the least-squares approach discussed above, and may optionally be smoothed or filtered based on successive estimated values.

Once the scaling parameters (two, in the case of dual-stream mode, or three, in the case of single-stream mode) are estimated, they may be applied to the impairment model terms to obtain a parametric estimate of impairment covariance matrix ({circumflex over (R)}u), according to Equation (10) or (12), as shown at block 580. At block 590, the estimated impairment covariance may be used to compute combining weights and determine symbol estimates.
Computing combining weights begins with calculating effective net channel estimates. For single stream, the effective net channel estimates ĥ(b) can be computed as:




  
    
      
        
          
            h
            ^
          
          ⁡
          
            (
            b
            )
          
        
        =
        
          
            ∑
            
              m
              =
              1
            
            N
          
          ⁢
          
            
              b
              m
            
            ⁢
            
              
                
                  h
                  ^
                
                m
              
              .
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




For dual stream:





  
    
      
        
          
            
              h
              ^
            
            ⁡
            
              (
              
                b
                k
              
              )
            
          
          =
          
            
              ∑
              
                m
                =
                1
              
              N
            
            ⁢
            
              
                b
                mk
              
              ⁢
              
                
                  h
                  ^
                
                m
              
            
          
        
        ,
        
          
            for
            ⁢
            
              
            
            ⁢
            k
          
          =
          
            1
            ,
            
              
            
            ⁢
            2.
          
        
      
    
    
      
        (
        22
        )
      
    
  




For single-stream mode, combining weights may be calculated by solving:

{circumflex over (R)}uwsingle=ĥ(b).  (23)

An iterative approach (Gauss-Seidel, conjugate gradient) or exact approach (LU, QR, Cholesky, etc.) may be used.

For dual-stream mode, per-code scaling terms are first computed via:




  
    
      
        
          
            
              
                α
                ^
              
              PC
            
            ⁡
            
              (
              n
              )
            
          
          =
          
            
              1
              K
            
            ⁢
            
              (
              
                
                  SF
                  data
                
                
                  SF
                  pilot
                
              
              )
            
            ⁢
            
              (
              
                
                  
                    γ
                    s
                  
                  ⁡
                  
                    (
                    n
                    )
                  
                
                
                  
                    γ
                    p
                  
                  ⁡
                  
                    (
                    1
                    )
                  
                
              
              )
            
            ⁢
            
              Γ
              
                D
                /
                P
              
            
          
        
        ,
        
          
            for
            ⁢
            
              
            
            ⁢
            n
          
          =
          
            1
            ,
            
              
            
            ⁢
            2.
          
        
      
    
    
      
        (
        24
        )
      
    
  




ΓD/P may be obtained by explicit signaling from the serving Node B, or by estimation. Given the per-code scaling, two sets of combining weights that account for cross-stream interference may be computed by solving:

{circumflex over (R)}stream1w1dual=ĥ(b1)

{circumflex over (R)}stream2w2dual=ĥ(b2),  (25)

where {circumflex over (R)}stream1={circumflex over (R)}u+{circumflex over (α)}PC(1)ĥ(b2)ĥH(b2) and {circumflex over (R)}stream2={circumflex over (R)}u+{circumflex over (α)}PC(2)ĥ(b1)ĥH(b1) represent the cross-stream interference at the first stream by the second stream, and at the second stream by the first stream, respectively.

The combining weight equations may be solved directly using iterative or exact approaches, or a two-step process can be employed, in which

{circumflex over (R)}uv1dual=ĥ(b1)

{circumflex over (R)}uv2dual=ĥ(b2)  (26)

are solved to obtain intermediate weights v1dual and v2dual. The final weights are computed according to:





  
    
      
        
          
            
              
                
                  w
                  1
                  dual
                
                =
                
                  
                    v
                    1
                    dual
                  
                  -
                  
                    
                      v
                      2
                      dual
                    
                    ⁡
                    
                      (
                      
                        
                          
                            
                              
                                α
                                ^
                              
                              PC
                            
                            ⁡
                            
                              (
                              2
                              )
                            
                          
                          ⁢
                          
                            
                              (
                              
                                v
                                2
                                dual
                              
                              )
                            
                            H
                          
                          ⁢
                          
                            
                              h
                              ^
                            
                            ⁡
                            
                              (
                              
                                b
                                1
                              
                              )
                            
                          
                        
                        
                          1
                          +
                          
                            
                              
                                
                                  α
                                  ^
                                
                                PC
                              
                              ⁡
                              
                                (
                                2
                                )
                              
                            
                            ⁢
                            
                              
                                (
                                
                                  v
                                  2
                                  dual
                                
                                )
                              
                              H
                            
                            ⁢
                            
                              
                                h
                                ^
                              
                              ⁡
                              
                                (
                                
                                  b
                                  2
                                
                                )
                              
                            
                          
                        
                      
                      )
                    
                  
                
              
            
          
          
            
              
                
                  w
                  2
                  dual
                
                =
                
                  
                    v
                    2
                    dual
                  
                  -
                  
                    
                      v
                      1
                      dual
                    
                    ⁡
                    
                      (
                      
                        
                          
                            
                              
                                α
                                ^
                              
                              PC
                            
                            ⁡
                            
                              (
                              1
                              )
                            
                          
                          ⁢
                          
                            
                              (
                              
                                v
                                1
                                dual
                              
                              )
                            
                            H
                          
                          ⁢
                          
                            
                              h
                              ^
                            
                            ⁡
                            
                              (
                              
                                b
                                2
                              
                              )
                            
                          
                        
                        
                          1
                          +
                          
                            
                              
                                
                                  α
                                  ^
                                
                                PC
                              
                              ⁡
                              
                                (
                                1
                                )
                              
                            
                            ⁢
                            
                              
                                (
                                
                                  v
                                  1
                                  dual
                                
                                )
                              
                              H
                            
                            ⁢
                            
                              
                                h
                                ^
                              
                              ⁡
                              
                                (
                                
                                  b
                                  1
                                
                                )
                              
                            
                          
                        
                      
                      )
                    
                  
                
              
            
          
        
        .
      
    
    
      
        (
        27
        )
      
    
  




The combining weights may then be used to detect and demodulate the transmitted symbols. The estimated covariance may also be used, for example, to estimate channel quality, i.e., to estimate or predict signal-to-noise-plus-interference or other channel quality metric. Periodically, such as at every WCDMA time slot, the process may be repeated, to dynamically adapt the covariance estimate to changing channel conditions and interference profiles.
The techniques described above may be adapted slightly to facilitate the estimation of SINR or other channel quality metrics for signal transmission scenarios other than a scenario currently employed, i.e. for a projected transmission scenario. Similarly, the techniques may be used to estimate SINR for projected transmission scenarios even when the mobile station is not scheduled for downlink transmissions at all. For example, Rstream1 and Rstream2 as estimated above correspond to a particular active transmission scenario (dual-stream) using a particular set of precoding vectors (b1 and b2). For Channel Quality Indicator (CQI) estimation purposes, the receiver may typically consider a number of transmitter modes that differ from the currently used mode. For instance, the current mode may be one of the four possible dual-stream modes shown in Equation (3); however, the receiver must also prepare a CQI report for one or more of the four possible single-stream modes shown in Equation (2). Such a requirement occurs in Release 7 of the UMTS standard, where the user equipment (UE) is required to periodically prepare single-stream CQI reports (referred to as Type B reports), regardless of what the UE's preferred mode is.
Those skilled in the art will appreciate that the impairment covariance model derived above for dual-stream mode (Equation (10)) is independent of the preceding vectors b1 and b2. As noted above, this same model may also be applied to single-stream modes, with some sacrifice in accuracy. As a result, an impairment covariance estimate obtained using the techniques described herein may be used to produce CQI estimates for any of the possible dual-stream or single-stream modes.
In one embodiment, an estimate for impairment covariance obtained, for example, according to:

{circumflex over (R)}u={circumflex over (μ)}(k1R11+k2R22)+k3(R11+R22)+{circumflex over (β)}Rn  (28)

may be used with effective net response coefficients ĥ(b) calculated for the target mode to estimate signal-to-noise-plus-interference (SINR) for the target mode. Thus, SINR for a target single-stream mode may be estimated according to:

SINR=(ĥ(b))H{circumflex over (R)}u−1ĥ(b)=wHĥ(b).  (29)

In this embodiment, the term Reff (b) may be ignored. For dual-stream mode:

SINRstream1=(ĥ(b1))H{circumflex over (R)}stream1−1ĥ(b1)=(w1dual)Hĥ(b1)

SINRstream2=(ĥ(b2))H{circumflex over (R)}stream2−1ĥ(b2)=(w2dual)Hĥ(b2),  (30)

where Rstream1 and Rstream2 are calculated as above according to the preceding matrix B for the projected mode.

If the quantity ΓD/P (the ratio of data energy to pilot energy across all antennas) is known or can be estimated, then another method may be used for determining CQI for projected transmission modes. If the impairment covariance is estimated for a single-stream mode, for example, then the estimate includes a term




  
    
      Γ
      
        D
        /
        P
      
    
    ⁡
    
      (
      
        
          1
          
            2
            ⁢
            
              N
              p
            
          
        
        ⁢
        
          
            R
            eff
          
          ⁡
          
            (
            
              b
              actual
            
            )
          
        
      
      )
    
  
  ,




where bactual represents the precoding vector for the currently scheduled mode. Because both ΓD/P and Np are known, this term can be calculated and removed (i.e., subtracted) from the overall impairment covariance estimate. Then, a new covariance estimate and associated SINR may be calculated for any projected transmission mode. For a single-stream mode, the projected impairment covariance is then:





  
    
      
        
          
            R
            
              u
              ⁢
              _
              ⁢
              projected
            
          
          =
          
            
              
                R
                ^
              
              u
            
            -
            
              
                Γ
                
                  D
                  /
                  P
                
              
              ⁡
              
                (
                
                  
                    1
                    
                      2
                      ⁢
                      
                        N
                        p
                      
                    
                  
                  ⁢
                  
                    
                      R
                      eff
                    
                    ⁡
                    
                      (
                      
                        b
                        actual
                      
                      )
                    
                  
                
                )
              
            
            +
            
              
                Γ
                
                  D
                  /
                  P
                
              
              ⁡
              
                (
                
                  
                    1
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        N
                        p
                      
                    
                  
                  ⁢
                  
                    
                      R
                      eff
                    
                    ⁡
                    
                      (
                      
                        b
                        projected
                      
                      )
                    
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        31
        )
      
    
  




where Ru is the estimated covariance for a currently scheduled single-stream mode employing precoding vector bactual, and bprojected is the precoding vector for a proposed, or projected single-stream mode. Those skilled in the art will immediately appreciate that SINR may be calculated for any desired mode using this general approach.

Accordingly, an exemplary method for estimating SINR for one or more projected transmission modes is illustrated in FIG. 6, and begins at block 610 with an estimated impairment covariance for a current mode. The estimated impairment covariance is adjusted at block 620 to remove the effects (if any) of a preceding vector corresponding to the current mode, and adjusted again at block 630 to add back in the effects of one or more preceding vectors corresponding to a projected transmission. The adjusted impairment covariance estimate is then used at block 640 to calculate a projected SINR for the projected transmission mode according to conventional techniques.
Although described herein in terms of own-cell interference, the present invention may be applied using additional model terms to model other-cell interference. For example, a single other-cell covariance term per dominant interfering base station may be added to equation (10) or (12), as described in the Cairns application. Additionally or alternatively, multiple covariance terms may be added to account for transmit diversity being used in other cells. Soft handoff may be handled, again as described in the Cairns application. Although described herein in terms of downlink reception, the present invention may be applied in the uplink as well.
Embodiments of the present invention thus provide a reduced-complexity method for MIMO interference suppression for both symbol-level (G-Rake) and chip-level (chip equalizer) LIW receivers, which are the two main architectures for advanced receivers in WCDMA systems. Those skilled in the art will appreciate that the particular design of a LIW receiver in accordance with the inventive techniques, and the associated nomenclature used in connection with such a receiver, may vary according to the network standard involved, but such variations are not germane to understanding or explaining the present invention. Moreover, it should be understood that the networks and radio devices illustrated and discussed herein are simplified; actual implementations likely will have additional entities omitted herein for clarity.
Nevertheless, an exemplary mobile terminal 112 includes one or both of the exemplary receiver circuits 200 or 300, illustrated in FIGS. 2 and 3, respectively. These receiver circuits may be implemented using various processing circuits, including A/D converters, filters, DSPs or other digital processors, memory, and the like. In at least one exemplary embodiment, mobile terminal 112 includes one or more DSPs and/or Application Specific Integrated Circuts (ASICS) or other programmable devices to implement receiver 112 including a G-Rake receiver as illustrated in FIG. 2. The processing circuits may be configured to include processing logic to carry out one or more of the methods described herein. It should thus be understood that at least a portion of the present invention's functionality may be embodied as stored computer instructions in the form of micro-code, firmware, software, etc.
More generally, the present invention can be implemented in hardware, software, or essentially any combination thereof, according to the needs of a particular design. Although the present invention has been described herein with respect to particular features, aspects and embodiments thereof, it will be apparent that numerous variations, modifications, and other embodiments are possible within the broad scope of the present invention, and accordingly, all variations, modifications and embodiments are to be regarded as being within the scope of the invention. The present embodiments are therefore to be construed in all aspects as illustrative and not restrictive and all changes coming within the meaning and equivalency range of the appended claims are intended to be embraced therein.

## Claims

1. A method of estimating impairment covariance associated with a first received composite information signal comprising at least a first data stream transmitted from first and second antennas according to a first antenna weighting vector, the method comprising:
constructing a first impairment model comprising at least first and second model terms corresponding to the first and second antennas respectively, wherein the model terms capture propagation channel effects and are independent of the first antenna weighting vector, and wherein the first impairment model omits cross-antenna interference terms; and
computing a first parametric estimate of the impairment covariance, using the first impairment model, by calculating the first and second model terms as a function of propagation channel estimates corresponding to the first and second antennas, respectively, measuring impairment covariance or data covariance associated with the received composite information signal, and estimating values for two or more unknown parameters of the first impairment model using the measured impairment covariance or data covariance and the computed first, second, and third model terms.

2. The method of claim 1, wherein constructing the first impairment model comprises computing the first model term as a function of first propagation channel estimates corresponding to the first antenna, the second model term as a function of second propagation channel estimates corresponding to the second antenna, and a third model term modeling noise and inter-cell interference.

3. The method of claim 2, wherein constructing the first impairment model further comprises:
grouping the computed first, second, and third model terms so that only two scaling parameters associated with the grouped first, second, and third model terms are unknown;
measuring impairment covariance or data covariance associated with the received composite information signal; and
estimating values for the two unknown scaling parameters using the measured impairment covariance or data covariance and the computed first, second, and third model terms.

4. The method of claim 1, further comprising calculating processing weights as a function of the first parametric estimate of the impairment covariance and propagation channel estimates corresponding to the first and second antennas.

5. The method of claim 4, wherein the processing weights comprise combining weights for use in a symbol-level equalizer.

6. The method of claim 4, wherein the processing weights comprise filter weights for use in a chip-level equalizer.

7. The method of claim 1, further comprising calculating a signal-to-interference-plus-noise ratio (SINR) estimate as a function of the first parametric estimate of the impairment covariance and first effective net channel response coefficients calculated from propagation channel estimates corresponding to each of the first and second antennas and the first antenna weighting vector.

8. The method of claim 7, further comprising calculating a projected SINR corresponding to a projected transmitted signal configuration, the projected SINR calculated as a function of second effective net channel response coefficients calculated from the propagation channel estimates and a second antenna weighting vector corresponding to the projected transmitted signal configuration.

9. The method of claim 1, wherein the first received composite information signal comprises a dual-stream MIMO signal and wherein the first parametric estimate of the impairment covariance is computed in a dual-stream mode, further comprising, for a second received composite information signal comprising a single-stream MIMO signal received in a single-stream mode:
computing a first impairment model term as a function of first propagation channel estimates corresponding only to the first antenna, a second impairment model term as a function of second propagation channel estimates corresponding only to the second antenna, and a third impairment model term as a function of the first and second propagation channel estimates;
grouping the computed first, second, and third impairment model terms so that only two scaling parameters associated with the grouped first, second, and third impairment model terms are unknown;
measuring impairment covariance or data covariance associated with the second received composite information signal;
estimating values for the two unknown scaling parameters using the measured impairment covariance or data covariance and the computed first, second, and third impairment model terms;
constructing a second impairment model based on the first, second, and third impairment model terms and the estimated values for the two unknown scaling parameters; and
computing a second parametric estimate of the impairment covariance for the second composite received signal using the second impairment model.

10. The method of claim 9, further comprising computing a fourth impairment model term corresponding to noise and other-cell interference and estimating a third scaling parameter corresponding to the fourth impairment model term jointly with estimating values for the two unknown scaling parameters, wherein the second impairment model is further based on the fourth impairment model term and the third scaling parameter.

11. The method of claim 9, further comprising calculating processing weights as a function of the second parametric estimate of the impairment covariance and propagation channel estimates corresponding to the first and second antennas.

12. The method of claim 9, further comprising calculating a signal-to-interference-plus-noise ratio (SINR) estimate as a function of the second parametric estimate of the impairment covariance and first effective net channel response coefficients calculated from propagation channel estimates corresponding to each of the first and second antennas and a single-stream antenna weighting vector mapping the single-stream MIMO signal to the first and second antennas.

13. The method of claim 12, further comprising calculating a projected SINR corresponding to a projected transmitted signal configuration by:
calculating a revised impairment covariance estimate based on the second parametric estimate of the impairment covariance, the single-stream antenna weighting vector, and a projected antenna weighting vector corresponding to the projected transmitted signal configuration; and
calculating the projected SINR as a function of the revised impairment covariance estimate.

14. A wireless communication receiver, comprising:
a radio front-end circuit configured to receive a first composite information signal, the first composite information signal comprising at least a first data stream transmitted from first and second transmit antennas according to a first antenna weighting vector; and
one or more processing circuits configured to:
construct a first impairment model for the impairment covariance comprising first and second model terms corresponding to the first and second antennas respectively, wherein the model terms capture propagation channel effects and are independent of the first antenna weighting vector, and wherein the impairment model does not include a cross-antenna interference term; and
compute a first parametric estimate of the impairment covariance, using the first impairment model, by calculating the first and second model terms as a function of propagation channel estimates corresponding to the first and second antennas, respectively, measuring impairment covariance or data covariance associated with the received composite information signal, and estimating values for two or more unknown parameters of the first impairment model using the measured impairment covariance or data covariance and the computed first, second, and third model terms.

15. The wireless communication receiver of claim 14, wherein the one or more processing circuits are configured to construct the first impairment model by computing the first model term as a function of first propagation channel estimates corresponding to the first antenna, the second model term as a function of second propagation channel estimates corresponding to the second antenna, and a third model term modeling noise and inter-cell interference.

16. The wireless communication receiver of claim 14, wherein the one or more processing circuits are further configured to construct the first impairment model by:
grouping the computed first, second, and third model terms so that only two scaling parameters associated with the grouped first, second, and third model terms are unknown;
measuring impairment covariance or data covariance associated with the first received composite information signal; and
estimating values for the two unknown scaling parameters using the measured impairment covariance or data covariance and the computed first, second, and third model terms.

17. The wireless communication receiver of claim 14, wherein the one or more processing circuits are further configured to calculate processing weights as a function of the first parametric estimate of the impairment covariance and propagation channel estimates corresponding to the first and second antennas.

18. The wireless communication receiver of claim 17, wherein the one or more processing circuits comprise a symbol-level equalizer, and wherein the processing weights comprise combining weights for use in the symbol-level equalizer.

19. The wireless communication receiver of claim 17, wherein the one or more processing circuits comprise a chip-level equalizer, and wherein the processing weights comprise filter weights for use in the chip-level equalizer.

20. The wireless communication receiver of claim 14, wherein the one or more processing circuits are further configured to calculate a signal-to-interference-plus-noise ratio (SINK) estimate as a function of the first parametric estimate of the impairment covariance.

21. The wireless communication receiver of claim 14, wherein the first received composite information signal comprises a dual-stream MIMO signal, wherein the one or more processing circuits are further configured to, for a second received composite information signal comprising a single-stream MIMO signal:
compute a first impairment model term as a function of first propagation channel estimates corresponding only to the first antenna, a second impairment model term as a function of second propagation channel estimates corresponding only to the second antenna, and a third impairment model term as a function of the first and second propagation channel estimates;
group the computed first, second, and third impairment model terms so that only two scaling parameters associated with the grouped first, second, and third impairment model terms are unknown;
measure impairment covariance or data covariance associated with the second received composite information signal;
estimate values for the two unknown scaling parameters using the measured impairment covariance or data covariance and the computed first, second, and third impairment model terms;
construct a second impairment model based on the first, second, and third impairment model terms and the estimated values for the two unknown scaling parameters; and
compute a second parametric estimate of the impairment covariance for the second composite received signal using the second impairment model.

22. The wireless communication receiver of claim 21, wherein the one or more processing circuits are further configured to compute a fourth impairment model term corresponding to noise and other-cell interference and to estimate a third scaling parameter corresponding to the fourth impairment model term jointly with estimating values for the two unknown scaling parameters, wherein the second impairment model is further based on the fourth impairment model term and the third scaling parameter.

23. The wireless communication receiver of claim 21, wherein the one or more processing circuits are further configured to calculate processing weights as a function of the second parametric estimate of the impairment covariance and propagation channel estimates corresponding to the first and second antennas.

24. The wireless communication receiver of claim 21, wherein the one or more processing circuits are further configured to calculate a signal-to-interference-plus-noise ratio (SINR) estimate as a function of the second parametric estimate of the impairment covariance and first effective net channel response coefficients calculated from propagation channel estimates corresponding to each of the first and second antennas and a single-stream antenna weighting vector mapping the single-stream MIMO signal to the first and second antennas.

25. The wireless communication receiver of claim 24, wherein the one or more processing circuits are further configured to calculate a projected SINR corresponding to a projected transmitted signal configuration by:
calculating a revised impairment covariance estimate based on the second parametric estimate of the impairment covariance, the single-stream antenna weighting vector, and a projected antenna weighting vector corresponding to the projected transmitted signal configuration; and
calculating the projected SINR as a function of the revised impairment covariance estimate.

26. A method of estimating impairment covariance associated with a first received composite information signal comprising at least a first data stream transmitted from first and second antennas according to an antenna weighting vector, the method comprising:
computing a first model term as a function of first propagation channel estimates corresponding to the first antenna, a second model term as a function of second propagation channel estimates corresponding to the second antenna, and a third model term modeling noise and inter-cell interference, wherein the model terms are independent of the antenna weighting vector;
measuring impairment covariance or data covariance associated with the received composite information signal; and
estimating values for two unknown scaling parameters for an impairment model, using the measured impairment covariance or data covariance and the computed first, second, and third model terms, wherein the impairment model comprises the first, second, and third model terms, grouped so that only two scaling parameters associated with the grouped first, second, and third model terms are unknown, and wherein the impairment model omits cross-antenna interference terms;
computing a first parametric estimate of the impairment covariance using the impairment model and the estimated values for the two unknown scaling parameters.