# Patent

## Title
Efficient computation of spatial filter matrices for steering transmit diversity in a MIMO communication system

## Abstract

Techniques for efficiently computing spatial filter matrices are described. The channel response matrices for a MIMO channel may be highly correlated if the channel is relatively static over a range of transmission spans. In this case, an initial spatial filter matrix may be derived based on one channel response matrix, and a spatial filter matrix for each transmission span may be computed based on the initial spatial filter matrix and a steering matrix used for that transmission span. The channel response matrices may be partially correlated if the MIMO channel is not static but does not change abruptly. In this case, a spatial filter matrix may be derived for one transmission span l and used to derive an initial spatial filter matrix for another transmission span m. A spatial filter matrix for transmission span m may be computed based on the initial spatial filter matrix, e.g., using an iterative procedure.

## Background
BACKGROUND
I. Field
The present invention relates generally to communication, and more specifically to spatial processing for data transmission in a multiple-input multiple-output (MIMO) communication system.
II. Background
A MIMO system employs multiple (NT) transmit antennas at a transmitting entity and multiple (NR) receive antennas at a receiving entity for data transmission. A MIMO channel formed by the NT transmit antennas and NR receive antennas may be decomposed into NS spatial channels, where NS≦min {NT, NR}. The NS spatial channels may be used to transmit data in parallel to achieve higher throughput and/or redundantly to achieve greater reliability.
Each spatial channel may experience various deleterious channel conditions such as, e.g., fading, multipath, and interference effects. The NS spatial channels may also experience different channel conditions and may achieve different signal-to-noise-and-interference ratios (SNRs). The SNR of each spatial channel determines its transmission capacity, which is typically quantified by a particular data rate that may be reliably transmitted on the spatial channel. For a time variant wireless channel, the channel conditions change over time and the SNR of each spatial channel also changes over time.
To improve performance, the MIMO system may utilize some form of feedback whereby the receiving entity evaluates the spatial channels and provides feedback information indicating the channel condition or the transmission capacity of each spatial channel. The transmitting entity may then adjust the data transmission on each spatial channel based on the feedback information. However, this feedback information may not be available for various reasons. For example, the system may not support feedback transmission from the receiving entity, or the wireless channel may change more rapidly than the rate at which the receiving entity can estimate the wireless channel and/or send back the feedback information. In any case, if the transmitting entity does not know the channel condition, then it may need to transmit data at a low rate so that the data transmission can be reliably decoded by the receiving entity even with the worst-case channel condition. The performance of such a system would be dictated by the expected worst-case channel condition, which is highly undesirable.
To improve performance (e.g., when feedback information is not available), the transmitting entity may perform spatial processing such that the data transmission does not observe the worst-case channel condition for an extended period of time, as described below. A higher data rate may then be used for the data transmission. However, this spatial processing represents additional complexity for both the transmitting and receiving entities.
There is therefore a need in the art for techniques to efficiently perform spatial processing to improve performance in a MIMO system.
SUMMARY
Techniques for efficiently computing spatial filter matrices used for spatial processing by a receiving entity are described herein. A transmitting entity may transmit data via a MIMO channel using either full channel state information (“full-CSI”) or “partial-CSI” transmission, as described below. The transmitting entity may also utilize steering transmit diversity (STD) for improved performance. With STD, the transmitting entity performs spatial processing with different steering matrices so that the data transmission observes an ensemble of effective channels and is not stuck on a “bad” channel realization for an extended period of time. The receiving entity performs the complementary receiver spatial processing for either full-CSI or partial-CSI transmission and for steering transmit diversity. The spatial filter matrices used for receiver spatial processing may be efficiently computed if the MIMO channel is relatively static or does not change abruptly.
If the MIMO channel is relatively static over a range of transmission spans (e.g., a range of symbol periods or frequency subbands), then the channel response matrices for the MIMO channel over these transmission spans may be highly correlated. In this case, an initial spatial filter matrix may be derived based on a channel response matrix and a selected receiver processing technique, as described below. A spatial filter matrix for each transmission span within the static range may then be computed based on the initial spatial filter matrix and the steering matrix used for that transmission span.
If the MIMO channel is not static but does not change abruptly, then the channel response matrices for different transmission spans may be partially correlated. In this case, a spatial filter matrix Mx(l) may be derived for a given transmission span l and used to derive an initial spatial filter matrix for another transmission span m. A spatial filter matrix Mx(m) for transmission span m may then be computed based on the initial spatial filter matrix, e.g., using an iterative procedure. The same processing may be repeated over a range of transmission spans of interest, so that each newly derived spatial filter matrix may be used to compute another spatial filter matrix for another transmission span.
The steering matrices may be defined such that the computation of the spatial filter matrices can be simplified. Various aspects and embodiments of the invention are described in further detail below.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a transmitting entity and a receiving entity in a MIMO system;
FIG. 2 shows a model for data transmission with steering transmit diversity;
FIGS. 3A and 3B show data transmission in a single-carrier MIMO system and a multi-carrier MIMO system, respectively;
FIGS. 4 and 5 show processes to compute spatial filter matrices for fully and partially correlated channel response matrices, respectively;
FIG. 6 shows a block diagram of an access point and a user terminal; and
FIG. 7 shows a block diagram of a processor for spatial filter matrix computation.



DETAILED DESCRIPTION
The word “exemplary” is used herein to mean “serving as an example, instance, or illustration.” Any embodiment described herein as “exemplary” is not necessarily to be construed as preferred or advantageous over other embodiments.
FIG. 1 shows a simple block diagram of a transmitting entity 110 and a receiving entity 150 in a MIMO system 100. At transmitting entity 110, a transmit (TX) spatial processor 120 performs spatial processing on data symbols (denoted by a vector s(m)) to generate transmit symbols (denoted by a vector x(m)). As used herein, a “data symbol” is a modulation symbol for data, a “pilot symbol” is a modulation symbol for pilot (which is data that is known a priori by both the transmitting and receiving entities), a “transmit symbol” is a symbol to be sent from a transmit antenna, a “received symbol” is a symbol obtained from a receive antenna, and a modulation symbol is a complex value for a point in a signal constellation used for a modulation scheme (e.g., M-PSK, M-QAM, and so on). The spatial processing is performed based on steering matrices V(m) and possibly other matrices. The transmit symbols are further conditioned by a transmitter unit (TMTR) 122 to generate NT modulated signals, which are transmitted from NT transmit antennas 124 and via a MIMO channel.
At receiving entity 150, the transmitted modulated signals are received by NR receive antennas 152, and the NR received signals are conditioned by a receiver unit (RCVR) 154 to obtain received symbols (denoted by a vector r(m)). A receive (RX) spatial processor 160 then performs receiver spatial processing (or spatial matched filtering) on the received symbols with spatial filter matrices (denoted by Mx(m)) to obtain “detected” data symbols (denoted by a vector {circumflex over (s)}(m)). The detected data symbols are estimates of the data symbols sent by transmitting entity 110. The spatial processing at the transmitting and receiving entities are described below.
The spatial filter matrix computation techniques described herein may be used for a single-carrier MIMO system as well as a multi-carrier MIMO system. Multiple carriers may be obtained with orthogonal frequency division multiplexing (OFDM), discrete multi tone (DMT), some other multi-carrier modulation techniques, or some other construct. OFDM effectively partitions the overall system bandwidth into multiple (NF) orthogonal subbands, which are also referred to as tones, subcarriers, bins, and frequency channels. With OFDM, each subband is associated with a respective subcarrier that may be modulated with data.
In MIMO system 100, the MIMO channel formed by the NT transmit antennas at transmitting entity 110 and the NR receive antennas at receiving entity 150 may be characterized by an NR×NT channel response matrix H(m), which may be given as:




  
    
      
        
          
            
              H
              _
            
            ⁡
            
              (
              m
              )
            
          
          =
          
            [
            
              
                
                  
                    
                      h
                      
                        1
                        ,
                        1
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                
                  
                    
                      h
                      
                        1
                        ,
                        2
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                
                  ⋯
                
                
                  
                    
                      h
                      
                        1
                        ,
                        
                          N
                          T
                        
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
              
              
                
                  
                    
                      h
                      
                        2
                        ,
                        1
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                
                  
                    
                      h
                      
                        2
                        ,
                        2
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                
                  ⋯
                
                
                  
                    
                      h
                      
                        2
                        ,
                        
                          N
                          T
                        
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
              
              
                
                  ⋮
                
                
                  ⋮
                
                
                  ⋰
                
                
                  ⋮
                
              
              
                
                  
                    
                      h
                      
                        
                          N
                          R
                        
                        ,
                        1
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                
                  
                    
                      h
                      
                        
                          N
                          R
                        
                        ,
                        2
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                
                  ⋯
                
                
                  
                    
                      h
                      
                        
                          N
                          R
                        
                        ,
                        
                          N
                          T
                        
                      
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
              
            
            ]
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          1
          )
        
      
    
  




where entry hi,j(m), for i=1 . . . NR and j=1 . . . NT, denotes the coupling or complex channel gain between transmit antenna j and receive antenna i for transmission span m. A transmission span may cover time and/or frequency dimensions. For example, in a single-carrier MIMO system, a transmission span may correspond to one symbol period, which is the time interval to transmit one data symbol. In a multi-carrier MIMO system, a transmission span may correspond to one subband in one symbol period. A transmission span may also cover multiple symbol periods and/or multiple subbands. For simplicity, the MIMO channel is assumed to be full rank with NS=NT≦NR.

The MIMO system may support data transmission using one or more operating modes such as, for example, a “calibrated” mode and an “uncalibrated” mode. The calibrated mode may employ full-CSI transmission whereby data is transmitted on orthogonal spatial channels (or “eigenmodes”) of the MIMO channel. The uncalibrated mode may employ partial-CSI transmission whereby data is transmitted on spatial channels of the MIMO channel, e.g., from individual transmit antennas.
The MIMO system may also employ steering transmit diversity (STD) to improve performance. With STD, the transmitting entity performs spatial processing with steering matrices so that a data transmission observes an ensemble of effective channels and is not stuck on a single bad channel realization for an extended period of time. Consequently, performance is not dictated by the worst-case channel condition.
1. Calibrated Mode—Full-CSI Transmission
For full-CSI transmission, eigenvalue decomposition may be performed on a correlation matrix of H(m) to obtain NS eigenmodes of H(m), as follows:

R(m)=HH(m)·H(m)=E(m)·Λ(m)·EH(m),  Eq (2)

where

    
    
        R(m) is an NT×NT correlation matrix of H(m);
        E(m) is an NT×NT unitary matrix whose columns are eigenvectors of R(m);
        Λ(m) is an NT×NT diagonal matrix of eigenvalues of R(m); and “H” denotes a conjugate transpose.

A unitary matrix U is characterized by the property UH·U=I, where I is the identity matrix. The columns of a unitary matrix are orthogonal to one another, and each column has unit power. The matrix E(m) may be used for spatial processing by the transmitting entity to transmit data on the NS eigenmodes of H(m). The eigenmodes may be viewed as orthogonal spatial channels obtained through decomposition. The diagonal entries of Λ(m) are eigenvalues of R(m), which represent the power gains for the NS eigenmodes. Singular value decomposition may also be performed to obtain matrices of left and right eigenvectors, which may be used for full-CSI transmission.

    
    


The transmitting entity performs spatial processing for full-CSI transmission with steering transmit diversity, as follows:

xf(m)=E(m)·V(m)·s(m),  Eq (3)

where s(m) is an NT×1 vector with up to NS data symbols to be sent in transmission span m;

    
    
        V(m) is an NT×NT steering matrix for transmission span m;
        E(m) is the matrix of eigenvectors for transmission span m; and
        xf(m) is an NT×1 vector with NT transmit symbols to be sent from the NT transmit antennas in transmission span m.

As shown in equation (3), each data symbol in s(m) is effectively spatially spread with a respective column of V(m). If NS<NT, then NS data symbols in s(m) are spatially spread with an NS×NS matrix V(m) to obtain NS“spread” symbols. Each spread symbol includes a component of each of the NS data symbols. The NS spread symbols from the spatial spreading are then sent on the NS eigenmodes of H(m). Each steering matrix V(m) is a unitary matrix and may be generated as described below.

    
    


The receiving entity obtains received symbols from the NR receive antennas, which may be expressed as:

rf(m)=H(m)·xf(m)+n(m)=H(m)·E(m)·V(m)·s(m)+n(m)=Hf—eff(m)·s(m)+n(m)  Eq (4)

where

    
    
        rf(m) is an NR×1 vector with NR received symbols obtained via the NR receive antennas in transmission span m;
        n(m) is a noise vector for transmission span m; and
        Hf—eff(m) is an NR×NT“effective” MIMO channel response matrix observed by the data vector s(m) for full-CSI transmission with steering transmit diversity, which is:

Hf—eff(m)=H(m)·E(m)·V(m).  Eq (5)

For simplicity, the noise is assumed to be additive white Gaussian noise (AWGN) with a zero mean vector and a covariance matrix of φnn=σ2·I, where σ2 is the variance of the noise and I is the identity matrix.

    
    


The receiving entity can recover the data symbols in s(m) using various receiver processing techniques. The techniques applicable for full-CSI transmission include a full-CSI technique and a minimum mean square error (MMSE) technique.
For the full-CSI technique, the receiving entity may derive a spatial filter matrix Mfsci(m) as follows:

Mfcsi(m)=VH(m)·Λ−1(m)·EH(m)·HH(m)  Eq (6)

The receiving entity may perform receiver spatial processing using Mfcsi(m), as follows:





  
    
      
        
          
            
              
                
                  s
                  _
                
                ^
              
              fcsi
            
            ⁡
            
              (
              m
              )
            
          
          =
          
            
              
                
                  M
                  _
                
                fcsi
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  r
                  _
                
                f
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                
                  V
                  _
                
                H
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  Λ
                  _
                
                
                  -
                  1
                
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  E
                  _
                
                H
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  H
                  _
                
                H
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              

            
            ⁢
            
              
            
            ⁢
            
              [
              
                
                  
                    
                      H
                      _
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ·
                  
                    
                      E
                      _
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ·
                  
                    
                      V
                      _
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ·
                  
                    
                      s
                      _
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                +
                
                  
                    n
                    _
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              ]
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                s
                _
              
              ⁡
              
                (
                m
                )
              
            
            +
            
              
                
                  n
                  _
                
                f
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          7
          )
        
      
    
  




where

    
    
        {circumflex over (s)}fcsi(m) is an NT×1 vector with NS detected data symbols; and
        nf(m) is the post-detection noise after the receiver spatial processing.
    
    


For the MMSE technique, the receiving entity may derive a spatial filter matrix Mf—mmse(m) as follows:

Mf—mmse(m)=[HfeffH(m)·Hf—eff(m)+σ2·I]−1·Hf—effH(m).  Eq (8)

The spatial filter matrix Mf—mmse(m) minimizes the mean square error between the symbol estimates from the spatial filter and the data symbols in s(m).

The receiving entity may perform MMSE spatial processing, as follows:




  
    
      
        
          
            
              
                
                  s
                  _
                
                ^
              
              f_mmse
            
            ⁡
            
              (
              m
              )
            
          
          =
          
            
              
                
                  D
                  _
                
                f_mmse
                
                  -
                  1
                
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  M
                  _
                
                f_mmse
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  r
                  _
                
                f
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                
                  D
                  _
                
                f_mmse
                
                  -
                  1
                
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  M
                  _
                
                f_mmse
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              

            
            ⁢
            
              
            
            ⁢
            
              [
              
                
                  
                    
                      
                        H
                        _
                      
                      f_eff
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ·
                  
                    
                      s
                      _
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                +
                
                  
                    n
                    _
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              ]
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                
                  
                    D
                    _
                  
                  f_mmse
                  
                    -
                    1
                  
                
                ⁡
                
                  (
                  m
                  )
                
              
              ·
              
                
                  
                    M
                    _
                  
                  f_mmse
                
                ⁡
                
                  (
                  m
                  )
                
              
              ·
              
                

              
              ⁢
              
                
              
              ⁢
              
                
                  
                    H
                    _
                  
                  f_eff
                
                ⁡
                
                  (
                  m
                  )
                
              
              ·
              
                
                  s
                  _
                
                ⁡
                
                  (
                  m
                  )
                
              
            
            +
            
              
                
                  n
                  _
                
                f_mmse
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          9
          )
        
      
    
  




where

    
    
        Df—mmse(M) is a diagonal matrix containing the diagonal elements of

Mf—mmse(m)·Hf—eff(m), or Df—mmse(m)=diag[Mf—mmse(m)·Hf—eff(m)]; and

        nf—mmse(m) is the MMSE filtered noise.

The symbol estimates from the spatial filter Mf—mmse(m) are unnormalized estimates of the data symbols. The multiplication with the scaling matrix Df—mmse(m) provides normalized estimates of the data symbols.

    
    


Full-CSI transmission attempts to send data on the eigenmodes of H(m). However, a full-CSI data transmission may not be completely orthogonal due to, for example, an imperfect estimate of H(m), error in the eigenvalue decomposition, finite arithmetic precision, and so on. The MMSE technique can account for (or “clean up”) any loss of orthogonality in the full-CSI data transmission.
Table 1 summarizes the spatial processing at the transmitting and receiving entities for full-CSI transmission with steering transmit diversity.









TABLE 1





Entity
Calibrated Mode - Full-CSI Transmission

















Transmitter
xf (m) = E (m) · V (m) · s (m)
Spatial




Processing



Hf—eff (m) = H (m) · E (m) · V (m)
Effective




Channel


Receiver
Mfcsi (m) = VH (m) · Λ−1 (m) · EH (m) ·
Spatial


full-CSI
HH (m)
Filter Matrix



ŝfcsi (m) = Mfcsi (m) · rf (m)
Spatial




Processing


Receiver
Mf—mmse (m) = [Hf—effH (m) ·
Spatial


MMSE
Hf—eff (m) + σ2 · I]−1 · Hf—effH (m)
Filter Matrix



Df—mmse (m) = diag [Mf—mmse (m) ·



Hf—eff (m)]



ŝf—mmse (m) = Df—mmse−1 (m) ·
Spatial



Mf—mmse (m) · rf (m)
Processing









2. Uncalibrated Mode—Partial-CSI Transmission

For partial-CSI transmission with steering transmit diversity, the transmitting entity performs spatial processing as follows:

xP(m)=V(m)·s(m),  Eq (10)

where xp(m) is the transmit data vector for transmission span m. As shown in equation (10), each data symbol in s(m) is spatially spread with a respective column of V(m). The NT spread symbols resulting from the multiplication with V(m) are then transmitted from the NT transmit antennas.

The receiving entity obtains received symbols, which may be expressed as:




  
    
      
        
          
            
              
                r
                _
              
              p
            
            ⁡
            
              (
              m
              )
            
          
          =
          
            
              
                
                  
                    H
                    _
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                ·
                
                  
                    
                      x
                      _
                    
                    p
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              +
              
                
                  n
                  _
                
                ⁡
                
                  (
                  m
                  )
                
              
            
            ⁢
            
              

            
            ⁢
            
              
            
            =
            
              
                
                  
                    H
                    _
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                ·
                
                  
                    V
                    _
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                ·
                
                  
                    s
                    _
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              +
              
                
                  n
                  _
                
                ⁡
                
                  (
                  m
                  )
                
              
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                
                  
                    H
                    _
                  
                  p_eff
                
                ⁡
                
                  (
                  m
                  )
                
              
              ·
              
                
                  s
                  _
                
                ⁡
                
                  (
                  m
                  )
                
              
            
            +
            
              
                n
                _
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          11
          )
        
      
    
  




where

    
    
        rp(m) is the received symbol vector for transmission span m; and
        Hp—eff(m) is an NR×NT effective MIMO channel response matrix observed by s(m) for partial-CSI transmission with steering transmit diversity, which is:

Hp—eff(m)=H(m)·V(m).  Eq (12)

    
    


The receiving entity can recover the data symbols in s(m) using various receiver processing techniques. The techniques applicable for partial-CSI transmission include a channel correlation matrix inversion (CCMI) technique (which is also commonly called a zero-forcing technique), the MMSE technique, and a successive interference cancellation (SIC) technique.
For the CCMI technique, the receiving entity may derive a spatial filter matrix Mccmi(m), as follows:

Mccmi(m)=[Hp—effH(m)·Hp—eff(m)]−1·Hp—effH(m)=Rp—eff−1(m)·Hp—effH(m).  Eq (13)

The receiving entity may perform CCMI spatial processing, as follows:





  
    
      
        
          
            
              
                
                  s
                  _
                
                ^
              
              ccmi
            
            ⁡
            
              (
              m
              )
            
          
          =
          
            
              
                
                  M
                  _
                
                ccmi
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  r
                  _
                
                p
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                
                  R
                  _
                
                p_eff
                
                  -
                  1
                
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  H
                  _
                
                p_eff
                H
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              

            
            ⁢
            
              
            
            ⁢
            
              [
              
                
                  
                    
                      
                        H
                        _
                      
                      p_eff
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ·
                  
                    
                      s
                      _
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                
                +
                
                  
                    n
                    _
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
              
              ]
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                s
                _
              
              ⁡
              
                (
                m
                )
              
            
            +
            
              
                
                  n
                  _
                
                ccmi
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          14
          )
        
      
    
  




where nccmi(m) is the CCMI filtered noise. Due to the structure of Rp—eff(m), the CCMI technique may amplify the noise.

For the MMSE technique, the receiving entity may derive a spatial filter matrix Mp—mmse(m), as follows:

Mp—mmse(m)=[Hp—effH(m)·Hp—eff(m)+σ2·I]−1Hp—effH(m).  Eq (15)

Equation (15) for the partial-CSI transmission has the same form as equation (8) for the full-CSI transmission. However, Hp—eff(m) (instead of Hf—eff(m)) is used in equation (15) for partial-CSI transmission.

The receiving entity may perform MMSE spatial processing, as follows:




  
    
      
        
          
            
              
                
                  s
                  _
                
                ^
              
              p_mmse
            
            ⁡
            
              (
              m
              )
            
          
          =
          
            
              
                
                  D
                  _
                
                p_mmse
                
                  -
                  1
                
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  M
                  _
                
                p_mmse
              
              ⁡
              
                (
                m
                )
              
            
            ·
            
              
                
                  r
                  _
                
                p
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
        
          

        
        ⁢
        
          
        
        ⁢
        
          =
          
            
              
                
                  
                    D
                    _
                  
                  p_mmse
                  
                    -
                    1
                  
                
                ⁡
                
                  (
                  m
                  )
                
              
              ·
              
                
                  
                    M
                    _
                  
                  p_mmse
                
                ⁡
                
                  (
                  m
                  )
                
              
              ·
              
                

              
              ⁢
              
                
              
              ⁢
              
                
                  
                    H
                    _
                  
                  p_eff
                
                ⁡
                
                  (
                  m
                  )
                
              
              ·
              
                
                  s
                  _
                
                ⁡
                
                  (
                  m
                  )
                
              
            
            +
            
              
                
                  n
                  _
                
                p_mmse
              
              ⁡
              
                (
                m
                )
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          16
          )
        
      
    
  




where Dp—mmse(m)=diag [Mp—mmse(m)·Hp—eff(m)] and np—mmse(m) is the MMSE filtered noise for partial-CSI transmission.

For the SIC technique, the receiving entity recovers the data symbols in s(m) in successive stages. For clarity, the following description assumes that each element of s(m) and each element of rp(m) corresponds to one data symbol stream. The receiving entity processes the NR received symbol streams in rp(m) in NS successive stages to recover the NS data symbol streams in s(m). Typically, the SIC processing is such that one packet is recovered for one stream, and then another packet is recovered for another stream, and so on. For simplicity, the following description assumes NS=NT.
For each stage l, where l=1 . . . NS, the receiving entity performs receiver spatial processing on NR input symbol streams rpl(m) for that stage. The input symbol streams for the first stage (l=1) are the received symbol streams, or rpl(m)=rp(m). The input symbol streams for each subsequent stage (l=2 . . . NS) are modified symbol streams from a preceding stage. The receiver spatial processing for stage l is based on a spatial filter matrix Mxl(m), which may be derived based on a reduced effective channel response matrix Hp—effl(m) and further in accordance with the CCMI, MMSE, or some other technique. Hpeffl(m) contains NS−l+1 columns in Hp—eff(m) corresponding to NS−l+1 data symbol streams not yet recovered in stage l. The receiving entity obtains one detected data symbol stream {ŝl′} for stage l and further processes (e.g., demodulates, deinterleaves, and decodes) this stream to obtain a corresponding decoded data stream {{circumflex over (d)}l}.
The receiving entity next estimates the interference that data symbol stream {sl} causes to the other data symbol streams not yet recovered. To estimate the interference, the receiving entity processes (e.g., re-encodes, interleaves, and symbol maps) the decoded data stream {{circumflex over (d)}l} in the same manner performed by the transmitting entity for this stream and obtains a stream of “remodulated” symbols {{hacek over (s)}l}, which is an estimate of the data symbol stream {sl} just recovered. The receiving entity then performs spatial processing on the remodulated symbol stream with steering matrices V(m) and further multiplies the result with channel response matrices H(m) to obtain NR interference components il(m) caused by stream {sl}. The receiving entity then subtracts the NR interference components il(m) from the NR input symbol streams rpl(m) for the current stage l to obtain NR input symbol streams rpl+1(m) for the next stage, or rpl−1(m)=rpl(m)−il(m). The input symbol streams rpl+1(m) represent the streams that the receiving entity would have received if the data symbol stream {sl} had not been transmitted, assuming that the interference cancellation was effectively performed. The receiving entity then repeats the same processing on the NR input symbol streams rpl+1(m) to recover another data stream. However, the effective channel response matrix Hp—effl+1(m) for the subsequent stage l+1 is reduced by one column corresponding to the data symbol stream {sl} recovered in stage l.
For the SIC technique, the SNR of each data symbol stream is dependent on (1) the receiver processing technique (e.g., CCMI or MMSE) used for each stage, (2) the specific stage in which the data symbol stream is recovered, and (3) the amount of interference due to the data symbol streams recovered in later stages. In general, the SNR progressively improves for data symbol streams recovered in later stages because the interference from data symbol streams recovered in prior stages is canceled. This may then allow higher rates to be used for data symbol streams recovered in later stages.
Table 2 summarizes the spatial processing at the transmitting and receiving entities for partial-CSI transmission with steering transmit diversity. For simplicity, the SIC technique is not shown in Table 2.









TABLE 2






Uncalibrated Mode - Partial-CSI



Entity
Transmission

















Transmitter
xp (m) = V (m) · s (m)
Spatial




Processing



Hp—eff (m) = H (m) · V (m)
Effective




Channel


Receiver
Mccmi (m) = [Hp—effH (m) · Hp—eff (m)]−1 ·
Spatial


CCMI
Hp—effH (m)
Filter Matrix



ŝccmi (m) = Mccmi (m) · rp (m)
Spatial




Processing


Receiver
Mp—mmse (m) = [Hp—effH (m) ·
Spatial


MMSE
Hp—eff (m) + σ2 · I]−1 · Hp—effH (m)
Filter Matrix



Dp—mmse (m) = diag [Mp—mmse (m) ·



Hp—eff (m)]



ŝp—mmse (m) = Dp—mmse−1 (m) ·
Spatial



Mp—mmse (m) · rp (m)
Processing









FIG. 2 shows a model for data transmission with steering transmit diversity. Transmitting entity 110 performs spatial processing (or spatial spreading) for steering transmit diversity (block 220) and spatial processing for either full-CSI or partial-CSI transmission (block 230). Receiving entity 150 performs receiver spatial processing for full-CSI or partial-CSI transmission (block 260) and receiver spatial processing (or spatial despreading) for steering transmit diversity (block 270). As shown in FIG. 2, the transmitting entity performs spatial spreading for steering transmit diversity prior to the spatial processing (if any) for full-CSI and partial-CSI transmission. The receiving entity may perform the complementary receiver spatial processing for full-CSI or partial-CSI transmission followed by spatial despreading for steering transmit diversity.
3. Spatial Filter Matrix Computation
With steering transmit diversity, different steering matrices V(m) may be used for different transmission spans to randomize the effective MIMO channel observed by a data transmission. This may then improve performance since the data transmission does not observe a “bad” MIMO channel realization for an extended period of time. The transmission spans may correspond to symbol periods for a single-carrier MIMO system or subbands for a multi-carrier MIMO system.
FIG. 3A shows a partial-CSI transmission with steering transmit diversity for a single-carrier MIMO system. For this system, the transmission span index m may be equal to a symbol period index n (or m=n). One vector s(n) of data symbols may be transmitted in each symbol period n and spatially spread with a steering matrix V(n) selected for that symbol period. Each data symbol vector s(n) observes an effective MIMO channel response of Hp—eff(n)=H(n)·V(n) and is recovered using a spatial filter matrix Mx(n).
FIG. 3B shows a partial-CSI transmission with steering transmit diversity in a multi-carrier MIMO system. For this system, the transmission span index m may be equal to a subband index k (or m=k). For each symbol period, one vector s(k) of data symbols may be transmitted in each subband k and spatially spread with a steering matrix V(k) selected for that subband. Each data symbol vector s(k) observes an effective MIMO channel response of Hp—eff(k)=H(k)·V(k) and is recovered using a spatial filter matrix Mx(k). The vector s(k) and the matrices V(k), H(k), and Mx(k) are also a function of symbol period n, but this is not shown for simplicity.
As shown in FIGS. 3A and 3B, if different steering matrices are used for different transmission spans, then the spatial filter matrices used by the receiving entity are a function of the transmission span index m. This is true even if the channel response matrix H(m) is fixed or constant over a range of transmission spans. For example, in a multi-carrier MIMO system, H(k) may be fixed across a set of subbands for a flat fading MIMO channel with a flat frequency response. As another example, in a single-carrier MIMO system, H(n) may be fixed over a given time interval for a MIMO channel with no temporal fading. This time interval may correspond to all or a portion of the time duration used to transmit a block of data symbols that is coded and decoded as a block.
A degree of correlation typically exists between the channel response matrices for adjacent transmission spans, e.g., between H(m) and H(m±1). This correlation may be exploited to simplify the computation for the spatial filter matrices at the receiving entity. The computation is described below for two cases—full-correlation and partial-correlation.
A. Full Correlation
With full-correlation, the channel response matrix for the MIMO channel is fixed over a range of transmission span indices of interest, e.g., for m=1 . . . M, where M may be any integer value greater than one. Thus, H(1)=H(2)= . . . =H(M)=H.
For the full-CSI technique, the spatial filter matrix Mfcsi(m) with fully correlated channel response matrices may be expressed as:

Mfcsi(m)=VH(m)·Λ−1·EH·HH.  Eq (17)

The spatial filter matrix Mfcsi(m) may then be computed as:

Mfcsi(m)=VH(m)·Mfcsibase, for m=1 . . . M,  Eq (18)

where Mfcsi—base=Λ−1·EH·HH is a base spatial filter matrix, which is the spatial filter matrix for the full-CSI technique without steering transmit diversity. The base spatial filter matrix Mfcsi—base is not a function of transmission span m because the channel response matrix H is fixed. Equation (18) indicates that the spatial filter matrix Mfcsi(m) for each transmission span m may be obtained by pre-multiplying the base spatial filter matrix Mfcsi—base with the steering matrix VH(m) used for that transmission span.

Alternatively, the spatial filter matrix Mfcsi(m) may be computed as:

Mfcsi(m)=W1(m)·Mfcsi(1), for m=2 . . . M,  Eq (19)

where Mfcsi(1)=VH(1)·Λ−1·EH·HH and W1(m)=VH(m)·V(1), Equation (19) indicates that the spatial filter matrix Mfcsi(m) for each transmission span m may be obtained by pre-multiplying the spatial filter matrix Mccmi(1) for transmission span 1 with the matrix W1(m). The matrices W1(m), for m=2 . . . M, are unitary matrices, each of which is obtained by multiplying two unitary steering matrices V(m) and V(1). The matrices W1(m) may be pre-computed and stored in a memory.

For the MMSE technique for full-CSI transmission, the spatial filter matrix Mf—mmse(m) with fully correlated channel response matrices may be expressed as:




  
    
      
        
          
        
        ⁢
        
          
            
              
                
                  
                    
                      M
                      _
                    
                    f_mmse
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                =
                
                ⁢
                
                  
                    
                      [
                      
                        
                          
                            
                              H
                              _
                            
                            f_eff
                            H
                          
                          ⁢
                          
                            
                              (
                              m
                              )
                            
                            ·
                            
                              
                                
                                  H
                                  _
                                
                                f_eff
                              
                              ⁡
                              
                                (
                                m
                                )
                              
                            
                          
                        
                        +
                        
                          
                            σ
                            2
                          
                          ·
                          
                            I
                            _
                          
                        
                      
                      ]
                    
                    
                      -
                      1
                    
                  
                  ·
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      
                        H
                        _
                      
                      f_eff
                      H
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ,
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  
                    
                      [
                      
                        
                          
                            
                              
                                V
                                _
                              
                              H
                            
                            ⁡
                            
                              (
                              m
                              )
                            
                          
                          ·
                          
                            
                              E
                              _
                            
                            H
                          
                          ·
                          
                            
                              H
                              _
                            
                            H
                          
                          ·
                          
                            H
                            _
                          
                          ·
                          
                            E
                            _
                          
                          ·
                          
                            
                              V
                              _
                            
                            ⁡
                            
                              (
                              m
                              )
                            
                          
                        
                        +
                        
                          
                            σ
                            2
                          
                          ·
                          
                            I
                            _
                          
                        
                      
                      ]
                    
                    
                      -
                      1
                    
                  
                  ·
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      
                        
                          V
                          _
                        
                        H
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                    ·
                    
                      
                        E
                        _
                      
                      H
                    
                    ·
                    
                      
                        H
                        _
                      
                      H
                    
                  
                  ,
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  
                    
                      
                        V
                        _
                      
                      H
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ·
                  
                    
                      [
                      
                        
                          
                            
                              E
                              _
                            
                            H
                          
                          ·
                          
                            
                              H
                              _
                            
                            H
                          
                          ·
                          
                            H
                            _
                          
                          ·
                          
                            E
                            _
                          
                        
                        +
                        
                          
                            σ
                            2
                          
                          ·
                          
                            I
                            _
                          
                        
                      
                      ]
                    
                    
                      -
                      1
                    
                  
                  ·
                  
                    
                      E
                      _
                    
                    H
                  
                  ·
                  
                    
                      
                        H
                        _
                      
                      H
                    
                    .
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          20
          )
        
      
    
  




Equation (20) is derived using the properties: (A·B)−1=B−1·A−1 and V·VH=I. The term within bracket in the second equality in equation (20) may be expressed as:




  
    
    ⁢
    
      [
      
        
          
            
              V
              _
            
            H
          
          ·
          
            
              E
              _
            
            H
          
          ·
          
            
              H
              _
            
            H
          
          ·
          
            H
            _
          
          ·
          
            E
            _
          
          ·
          
            V
            _
          
        
        +
        
          
            σ
            2
          
          ·
          
            I
            _
          
        
      
      ]
    
  




  
    
      
        
          
            =
            
            ⁢
            
              [
              
                
                  
                    
                      V
                      _
                    
                    H
                  
                  ⁡
                  
                    (
                    
                      
                        
                          
                            E
                            _
                          
                          H
                        
                        ·
                        
                          
                            H
                            _
                          
                          H
                        
                        ·
                        
                          H
                          _
                        
                        ·
                        
                          E
                          _
                        
                      
                      +
                      
                        
                          σ
                          2
                        
                        ·
                        
                          V
                          _
                        
                        ·
                        
                          I
                          _
                        
                        ·
                        
                          
                            V
                            _
                          
                          H
                        
                      
                    
                    )
                  
                
                ·
                
                  V
                  _
                
              
              ]
            
          
          ,
        
      
    
    
      
        
          
            =
            
            ⁢
            
              [
              
                
                  
                    
                      
                        V
                        _
                      
                      H
                    
                    (
                    
                      
                        
                          
                            E
                            _
                          
                          H
                        
                        ·
                        
                          
                            H
                            _
                          
                          H
                        
                        ·
                        
                          H
                          _
                        
                        ·
                        
                          E
                          _
                        
                      
                      +
                      
                        
                          σ
                          2
                        
                        ·
                        
                          
                            I
                            )
                          
                          _
                        
                      
                    
                    ]
                  
                  
                    -
                    1
                  
                
                ·
                
                  V
                  _
                
              
              ]
            
          
          ,
        
      
    
  



where “(m)” has been omitted for clarity. The inverse of the term in the second equality above may then be expressed as:

[VH(EH·HH·H·E+σ2·I)·V]−1=VH[(EH·HH·H·E+σ2·I)−1·V],

where VH=V−.

The spatial filter matrix Mf—mmse(m) may be computed as:

Mf—mmse(m)=VH(m)·Mf—mmse—base, for m=1 . . . M,  Eq (21)

where Mf—mmse—base=[EH·HH·H·E+σ2·I]−1·EH·HH. Similar to the full-CSI technique, the spatial filter matrix Mf—mmse(m) for transmission span m may be obtained by pre-multiplying the base spatial filter matrix Mf—mmse—base with the steering matrix VH(m). The spatial filter matrix Mf—mmse(m) may also be computed as:

Mf—mmse(m)=W1(m)·Mf—mmse(1) for m=2 . . . M,  Eq (22)

where Mf—mmse(1)=VH(1)·[EH·HH·H·E+σ2·I]−1·EH·HH.

For the CCMI technique, the spatial filter matrix Mccmi(m) with fully correlated channel response matrices may be expressed as:




  
    
      
        
          
        
        ⁢
        
          
            
              
                
                  
                    
                      
                        M
                        _
                      
                      ccmi
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  =
                  
                  ⁢
                  
                    
                      
                        [
                        
                          
                            
                              H
                              _
                            
                            p_eff
                            H
                          
                          ⁢
                          
                            
                              (
                              m
                              )
                            
                            ·
                            
                              
                                
                                  H
                                  _
                                
                                p_eff
                              
                              ⁡
                              
                                (
                                m
                                )
                              
                            
                          
                        
                        ]
                      
                      
                        -
                        1
                      
                    
                    ·
                    
                      
                        
                          H
                          _
                        
                        p_eff
                        H
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                  
                
                ,
              
            
          
          
            
              
                
                  =
                  
                  ⁢
                  
                    
                      
                        [
                        
                          
                            
                              V
                              _
                            
                            H
                          
                          ⁢
                          
                            
                              (
                              m
                              )
                            
                            ·
                            
                              
                                H
                                _
                              
                              H
                            
                            ·
                            
                              H
                              _
                            
                            ·
                            
                              
                                V
                                _
                              
                              ⁡
                              
                                (
                                m
                                )
                              
                            
                          
                        
                        ]
                      
                      
                        -
                        1
                      
                    
                    ·
                    
                      
                        
                          V
                          _
                        
                        H
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                    ·
                    
                      
                        H
                        _
                      
                      H
                    
                  
                
                ,
              
            
          
          
            
              
                
                  =
                  
                  ⁢
                  
                    
                      
                        [
                        
                          
                            
                              
                                V
                                _
                              
                              H
                            
                            ⁡
                            
                              (
                              m
                              )
                            
                          
                          ·
                          
                            R
                            _
                          
                          ·
                          
                            
                              V
                              _
                            
                            ⁡
                            
                              (
                              m
                              )
                            
                          
                        
                        ]
                      
                      
                        -
                        1
                      
                    
                    ·
                    
                      
                        
                          V
                          _
                        
                        H
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                    ·
                    
                      
                        H
                        _
                      
                      H
                    
                  
                
                ,
              
            
          
          
            
              
                
                  =
                  
                  ⁢
                  
                    
                      
                        
                          V
                          _
                        
                        
                          -
                          1
                        
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                    ·
                    
                      
                        R
                        _
                      
                      
                        -
                        1
                      
                    
                    ·
                    
                      
                        [
                        
                          
                            
                              V
                              _
                            
                            H
                          
                          ⁡
                          
                            (
                            m
                            )
                          
                        
                        ]
                      
                      
                        -
                        1
                      
                    
                    ·
                    
                      
                        
                          V
                          _
                        
                        H
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                    ·
                    
                      
                        H
                        _
                      
                      H
                    
                  
                
                ,
              
            
          
          
            
              
                
                  =
                  
                  ⁢
                  
                    
                      
                        
                          V
                          _
                        
                        H
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                    ·
                    
                      
                        R
                        _
                      
                      
                        -
                        1
                      
                    
                    ·
                    
                      
                        H
                        _
                      
                      H
                    
                  
                
                ,
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          23
          )
        
      
    
  




where [VH(m)]−1=V(m) because V(m) is a unitary matrix.

The spatial filter matrix Mccmi(m) may thus be computed as:

Mccmi(m)=VH(m)·Mccmi—base, for m=1 . . . M,  Eq (24)

where Mccmi—base=R−1·H. The spatial filter matrix Mccmi(m) may also be computed as:

Mccmi(m)=W1(m)·Mccmi(1), for m=2 . . . M,  Eq (25)

where Mccmi(1)=VH(1)·R−1·HH.

For the MMSE technique for partial-CSI transmission, the spatial filter matrix Mp—mmse(m) with fully correlated channel response matrices may be expressed as:




  
    
      
        
          
        
        ⁢
        
          
            
              
                
                  
                    
                      M
                      _
                    
                    p_mmse
                  
                  ⁡
                  
                    (
                    m
                    )
                  
                
                =
                
                ⁢
                
                  
                    
                      [
                      
                        
                          
                            
                              H
                              _
                            
                            p_eff
                            H
                          
                          ⁢
                          
                            
                              (
                              m
                              )
                            
                            ·
                            
                              
                                
                                  H
                                  _
                                
                                p_eff
                              
                              ⁡
                              
                                (
                                m
                                )
                              
                            
                          
                        
                        +
                        
                          
                            σ
                            2
                          
                          ·
                          
                            I
                            _
                          
                        
                      
                      ]
                    
                    
                      -
                      1
                    
                  
                  ·
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      
                        H
                        _
                      
                      p_eff
                      H
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ,
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  
                    
                      [
                      
                        
                          
                            
                              
                                V
                                _
                              
                              H
                            
                            ⁡
                            
                              (
                              m
                              )
                            
                          
                          ·
                          
                            
                              H
                              _
                            
                            H
                          
                          ·
                          
                            H
                            _
                          
                          ·
                          
                            
                              V
                              _
                            
                            ⁡
                            
                              (
                              m
                              )
                            
                          
                        
                        +
                        
                          
                            σ
                            2
                          
                          ·
                          
                            I
                            _
                          
                        
                      
                      ]
                    
                    
                      -
                      1
                    
                  
                  ·
                
              
            
          
          
            
              
                
                ⁢
                
                  
                    
                      
                        
                          V
                          _
                        
                        H
                      
                      ⁡
                      
                        (
                        m
                        )
                      
                    
                    ·
                    
                      
                        H
                        _
                      
                      H
                    
                  
                  ,
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  
                    
                      
                        V
                        _
                      
                      H
                    
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ·
                  
                    
                      [
                      
                        
                          
                            
                              H
                              _
                            
                            H
                          
                          ·
                          
                            H
                            _
                          
                        
                        +
                        
                          
                            σ
                            2
                          
                          ·
                          
                            I
                            _
                          
                        
                      
                      ]
                    
                    
                      -
                      1
                    
                  
                  ·
                  
                    
                      
                        H
                        _
                      
                      H
                    
                    .
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          26
          )
        
      
    
  




Equation (26) may be derived in similar manner as equation (20) above.

The spatial filter matrix Mpmmse(m) may be computed as:

Mp—mmse(m)=VH(m)·Mpmmsebase, for m=1 . . . M,  Eq (27)

where Mpmmse—base=[HH·H+σ2·I]−1·HH. The spatial filter matrix Mp—mmse(m) may also be computed as:

Mp—mmse(m)=W1(m)·Mpmmse(1), for m=2 . . . M,  Eq (28)

where Mp—mmse(1)=VH(1)·[HH·H+H=σ2·I]−1·HH.

Table 3 summarizes the computation for the spatial filter matrices for full-CSI and partial-CSI transmissions with fully correlated channel response matrices over transmission spans m=1 . . . M .







TABLE 3







Spatial Filter Matrices with Full Correlation









Mode
Spatial Filter Matrix
Technique





Full-CSI
Mfcsi—base = Λ−1 · EH · HH, and
Full-CSI



Mfcsi (m) = VH (m) · Mfcsi—base



Mf—mmse—base = [EH · HH · H · E +
MMSE



σ2 · I]−1 · EH · HH, and



Mf—mmse—base (m) = VH (m) · Mf—mmse—base


Partial-CSI
Mccmi—base = R−1 · HH, and
CCMI



Mccmi (m) = VH (m) · Mccmi—base



Mp—mmse—base = [HH · H + σ2 · I]−1 · HH, and
MMSE



Mp—mmse (m) = VH (m) · Mp—mmse—base









In general, the spatial filter matrix for transmission span m may be computed as Mx(m)=VH(m)·Mx—base, where the subscript “x” denotes the receiver processing technique and may be “fcsi”, “f_mmse”, “ccmi”, or “p_mmse”. The base spatial filter matrix Mx—bass may be computed as if steering transmit diversity was not used.
FIG. 4 shows a flow diagram of a process 400 to compute spatial filter matrices with fully correlated channel response matrices over transmission spans m=1 . . . M . An initial spatial filter matrix Mxinit is first computed (block 412). This initial spatial filter matrix may be the base spatial filter matrix Mx—base that is derived based on (1) the channel response matrix H and (2) the receiver processing technique selected for use (e.g., full-CSI, MMSE for full-CSI, CCMI, or MMSE for partial-CSI). Alternatively, the initial spatial filter matrix may be the spatial filter matrix Mx(1) for transmission span m=1, which may be derived based on H and V(1).
The transmission span index m is then set to 1 if Mx—init=Mx—base (as shown in FIG. 4) or set to 2 if Mx—init=Mx(1) (block 414). The spatial filter matrix Mx(m) for transmission span m is then computed based on the initial spatial filter matrix Mx—init and the steering matrix V(m) used for transmission span m (block 416). In particular, Mx(m) may be computed based on either Mx—base and V(m) or Mx(1) and W1(m), as described above. A determination is then made whether m<M (block 420). If the answer is ‘yes, then the index m is incremented (block 422), and the process returns to block 416 to compute the spatial filter matrix for another transmission span. Otherwise, if m=M in block 420, then the spatial filter matrices Mx(1) through Mx(M) are used for receiver spatial processing of received symbol vectors rx(1) through rx(M), respectively (block 424). Although not shown in FIG. 4 for simplicity, each spatial filter matrix may be used for receiver spatial processing as soon as both the spatial filter matrix Mx(m) is generated and the received symbol vector rx(m) are obtained.
For full-CSI transmission, the spatial processing at the transmitting entity may also be simplified as: xf(m)=E·V(m)·s(m). A matrix E·V(m) may be computed for each transmission span m based on the steering matrix V(m) for that transmission span and the matrix E, which is not a function of transmission span for the full correlation case.
B. Partial Correlation
With partial-correlation, the channel response matrices for the MIMO channel are less than fully correlated across a range of transmission span indices of interest. In this case, a spatial filter matrix computed for a transmission span l may be used to facilitate the computation of a spatial filter matrix for another transmission span m.
In an embodiment, a base spatial filter matrix Mx—base(l) for transmission span l is obtained from a spatial filter matrix Mx(l) computed for transmission span l by removing the steering matrix V(l) used for transmission span l, as follows:

Mx—base(l)=V(l)·Mx(l).  Eq (29)

The base spatial filter matrix Mx—base(l) is then used to derive a base spatial filter matrix Mx—base(m) for transmission span m (e.g., m=l±1). Mx—base(m) may be computed, e.g., using an iterative procedure or algorithm that iteratively performs a set of computations on Mx—base(l) to obtain a final solution for Mx—base(m). Iterative procedures for computing an MMSE solution (e.g., adaptive MMSE algorithms, gradient algorithm, lattice algorithms, and so on) are known in the art and not described herein. The spatial filter matrix Mx(m) for transmission span m may be computed as:

Mx(m)=VH(m)·Mx—base(m).  Eq (30)

The processing order for this embodiment may thus be given as: Mx(l)→Mx—base(l)Mx—base(m)→Mx(m), where “→” denotes a direct computation and “” denotes possible iterative computation. The base spatial filter matrices Mx—base(l) and Mx—base(m) do not contain steering matrices, whereas the spatial filter matrices Mx(l) and Mx(m) contain steering matrices V(l) and V(m) used for transmission spans l and m, respectively.

In another embodiment, the spatial filter matrix Mx(m) for transmission span m is computed using an iterative procedure that iteratively performs a set of computations on an initial guess {tilde over (M)}x(m). The initial guess may be derived from the spatial filter matrix Mx(l) derived for transmission span l, as follows:

{tilde over (M)}x(m)=Wl(m)·Mx(l),  Eq (31)

where Wl(m)=VH(m)·V(l). The processing order for this embodiment may be given as: Mx(l)→{tilde over (M)}x(m)Mx(m). The spatial filter matrices {tilde over (M)}x(m) and Mx(m) both contain the steering matrix V(m) used for transmission span m.

For the above embodiments, Mx—base(l) and {tilde over (M)}x(m) may be viewed as the initial spatial filter matrices used to derive the spatial filter matrix Mx(m) for a new transmission span m. In general, the amount of correlation between Mx(l) and Mx(m) is dependent on the amount of correlation between Mx—base(l) and Mx—base(m), which is dependent on the amount of correlation between H(l) and H(m) for transmission spans l and m. A higher degree of correlation may result in faster convergence to the final solution for Mx(l).
FIG. 5 shows a flow diagram of a process 500 to compute spatial filter matrices with partially correlated channel response matrices for transmission spans m=1 . . . M . The indices for the current and next transmission spans are initialized as l=1 and m=2 (block 512). A spatial filter matrix Mx(l) is computed for transmission span f in accordance with the receiver processing technique selected for use (block 514). An initial spatial filter matrix Mx—init for transmission span m is then computed based on the spatial filter matrix Mx(l) and the proper steering matrix/matrices V(l) and V(m), e.g., as shown in equation (29) or (31) (block 516). The spatial filter matrix Mx(m) for transmission span m is then computed based on the initial spatial filter matrix Mx—init, X e.g., using an iterative procedure (block 518).
A determination is then made whether m<M (block 520). If the answer is ‘yes’, then the indices l and m are updated, e.g., as l=m and m=m+1 (block 522). The process then returns to block 516 to compute a spatial filter matrix for another transmission span. Otherwise, if all spatial filter matrices have been computed, as determined in block 520, then the spatial filter matrices Mx(1) through Mx(M) are used for receiver spatial processing of received symbol vectors rx(1) through rx(M), respectively (block 524).
For simplicity, FIG. 5 shows the computation of M spatial filter matrices for M consecutive transmission spans m=1 . . . M. The transmission spans do not need to be contiguous. In general, a spatial filter matrix derived for one transmission span l is used to obtain an initial guess of a spatial filter matrix for another transmission span m, where l and m may be any index values.
4. Steering Matrices
A set of steering matrices (or transmit matrices) may be generated and used for steering transmit diversity. These steering matrices may be denoted as {V}, or V(i) for i=1 . . . L, where L may be any integer greater than one. Each steering matrix V(i) should be a unitary matrix. This condition ensures that the NT data symbols transmitted simultaneously using V(i) have the same power and are orthogonal to one another after the spatial spreading with V(i).
The set of L steering matrices may be generated in various manners. For example, the L steering matrices may be generated based on a unitary base matrix and a set of scalars. The base matrix may be used as one of the L steering matrices. The other L−1 steering matrices may be generated by multiplying the rows of the base matrix with different combinations of scalars. Each scalar may be any real or complex value. The scalars are selected to have unit magnitude so that steering matrices generated with these scalars are unitary matrices.
The base matrix may be a Walsh matrix. A 2×2 Walsh matrix W2×2 and a larger size Walsh matrix W2N×2N may be expressed as:




  
    
      
        
          
            W
            _
          
          
            2
            ×
            2
          
        
        =
        
          
            
              [
              
                
                  
                    1
                  
                  
                    
                      
                        
                      
                      ⁢
                      1
                    
                  
                
                
                  
                    1
                  
                  
                    
                      -
                      1
                    
                  
                
              
              ]
            
            ⁢
            
              
            
            ⁢
            and
            ⁢
            
              
            
            ⁢
            
              
                W
                _
              
              
                2
                ⁢
                N
                ×
                2
                ⁢
                N
              
            
          
          =
          
            
              [
              
                
                  
                    
                      
                        W
                        _
                      
                      
                        N
                        ×
                        N
                      
                    
                  
                  
                    
                      
                        
                      
                      ⁢
                      
                        
                          W
                          _
                        
                        
                          N
                          ×
                          N
                        
                      
                    
                  
                
                
                  
                    
                      
                        W
                        _
                      
                      
                        N
                        ×
                        N
                      
                    
                  
                  
                    
                      -
                      
                        
                          W
                          _
                        
                        
                          N
                          ×
                          N
                        
                      
                    
                  
                
              
              ]
            
            .
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          32
          )
        
      
    
  




Walsh matrices have dimensions that are powers of two (e.g., 2, 4, 8, and so on).

The base matrix may also be a Fourier matrix. For an N×N Fourier matrix DN×N, the element dn,m in the n-th row and m-th column of DN×N may be expressed as:




  
    
      
        
          
            d
            
              n
              ,
              m
            
          
          =
          
            e
            
              
                -
                j2π
              
              ⁢
              
                
                  
                    (
                    
                      n
                      -
                      1
                    
                    )
                  
                  ⁢
                  
                    (
                    
                      m
                      -
                      1
                    
                    )
                  
                
                N
              
            
          
        
        ,
        
          
            for
            ⁢
            
              
            
            ⁢
            n
          
          =
          
            
              
                {
                
                  1
                  ⁢
                  
                    
                  
                  ⁢
                  ...
                  ⁢
                  
                    
                  
                  ⁢
                  N
                
                }
              
              ⁢
              
                
              
              ⁢
              and
              ⁢
              
                
              
              ⁢
              m
            
            =
            
              
                {
                
                  1
                  ⁢
                  
                    
                  
                  ⁢
                  ...
                  ⁢
                  
                    
                  
                  ⁢
                  N
                
                }
              
              .
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          33
          )
        
      
    
  




Fourier matrices of any square dimension (e.g., 2, 3, 4, 5, and so on) may be formed. Other matrices may also be used as the base matrix.

For an N×N base matrix, each of rows 2 through N of the base matrix may be independently multiplied with one of K different possible scalars. KN−1 different steering matrices may be obtained from KN−1 different permutations of the K scalars for N−1 rows. For example, each of rows 2 through N may be independently multiplied with a scalar of +1, −1, +j, or −j. For N=4 and K=4, 64 different steering matrices may be generated from a 4×4 base matrix with four different scalars. In general, each row of the base matrix may be multiplied with any scalar having the form ejθ, where θ may be any phase value. Each element of a scalar-multiplied N×N base matrix is further scaled by 1/√{square root over (N)} to obtain an N×N steering matrix having unit power for each column.
Steering matrices derived based on a Walsh matrix (or a 4×4 Fourier matrix) have certain desirable properties. If the rows of the Walsh matrix are multiplied with scalars of ±1 and ±j, then each element of a resultant steering matrix is +1, −1, +j, or −j. In this case, the multiplication of an element (or “weight”) of a spatial filter matrix with an element of the steering matrix may be performed with just bit manipulation. If the elements of the L steering matrices belong in a set composed of {+1, −1, +j, −j}, then the computation to derive the spatial filter matrices for the full correlation case can be greatly simplified.
5. MIMO System
FIG. 6 shows a block diagram of an access point 610 and a user terminal 650 in a MIMO system 600. Access point 610 is equipped with Nap antennas that may be used for data transmission and reception, and user terminal 650 is equipped with Nut antennas, where Nap>1 and Nut>1.
On the downlink, at access point 610, a TX data processor 620 receives and processes (encodes, interleaves, and symbol maps) traffic/packet data and control/overhead data and provides data symbols. A TX spatial processor 630 performs spatial processing on the data symbols with steering matrices V(m) and possibly eigenvector matrices E(m) for the downlink, e.g., as shown in Tables 1 and 2. TX spatial processor 630 also multiplexes in pilot symbols, as appropriate, and provides Nap streams of transmit symbols to Nap transmitter units 632a through 632ap. Each transmitter unit 632 receives and processes a respective transmit symbol stream and provides a corresponding downlink modulated signal. Nap downlink modulated signals from transmitter units 632a through 632ap are transmitted from Nap antennas 634a through 634ap, respectively.
At user terminal 650, Nut antennas 652a through 652ut receive the transmitted downlink modulated signals, and each antenna provides a received signal to a respective receiver unit 654. Each receiver unit 654 performs processing complementary to that performed by receiver unit 632 and provides received symbols. An RX spatial processor 660 performs receiver spatial processing on the received symbols from all Nut receiver units 654a through 654ut, e.g., as shown in Tables 1 and 2, and provides detected data symbols. An RX data processor 670 processes (e.g., symbol demaps, deinterleaves, and decodes) the detected data symbols and provides decoded data for the downlink.
The processing for the uplink may be the same or different from the processing for the downlink. Traffic and control data is processed (e.g., encoded, interleaved, and symbol mapped) by a TX data processor 688, spatially processed by a TX spatial processor 690 with steering matrices V(m) and possibly eigenvector matrices E(m) for the uplink, and multiplexed with pilot symbols to generate Nut transmit symbol streams. Nut transmitter units 654a through 654ut condition the Nut transmit symbol streams to generate Nut uplink modulated signals, which are transmitted via Nut antennas 652a through 652ut. 
At access point 610, the uplink modulated signals are received by Nap antennas 634a through 634ap and processed by Nap receiver units 632a through 632ap to obtain received symbols for the uplink. An RX spatial processor 644 performs receiver spatial processing on the received symbols and provides detected data symbols, which are further processed by an RX data processor 646 to obtain decoded data for the uplink.
Processors 638 and 678 perform channel estimation and spatial filter matrix computation for the access point and user terminal, respectively. Controllers 640 and 680 control the operation of various processing units at the access point and user terminal, respectively. Memory units 642 and 682 store data and program codes used by controllers 630 and 680, respectively.
FIG. 7 shows an embodiment of processor 678, which performs channel estimation and spatial filter matrix computation for user terminal 650. A channel estimator 712 obtains received pilot symbols and derives a channel response matrix for each transmission span in which received pilot symbols are available. A filter 714 may perform time-domain filtering of the channel response matrices for the current and prior transmission spans to obtain a higher quality channel response matrix H(m). A unit 716 then computes an initial spatial filter matrix Mx—init.
For fully correlated H(m), the initial spatial filter matrix Mx—init may be (1) a base spatial filter matrix Mx—base computed based on H(m) and the selected receiver processing technique or (2) a spatial filter matrix Mx(1) for transmission span 1 computed based on H(1), V(1), and the selected receiver processing technique. For partially correlated H(m), the initial spatial filter matrix Mx—init may be an initial guess Mx—base(l) or {tilde over (M)}x(m) that is obtained based on a spatial filter matrix Mx(l) computed for another transmission span t. A unit 718 computes the spatial filter matrix Mx(m) for transmission span m based on the initial spatial filter matrix Mx—init and the steering matrix V(m) used for that transmission span. For partially correlated H(m), unit 718 may implement an iterative procedure to compute for Mx(m) based on the initial spatial filter matrix, which is an initial guess of Mx(m).
Processor 638 performs channel estimation and spatial filter matrix computation for access point 610 and may be implemented in similar manner as processor 678.
The spatial filter matrix computation techniques described herein may be implemented by various means. For example, these techniques may be implemented in hardware, software, or a combination thereof. For a hardware implementation, the processing units for spatial filter matrix computation may be implemented within one or more application specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), processors, controllers, micro-controllers, microprocessors, other electronic units designed to perform the functions described herein, or a combination thereof.
For a software implementation, the spatial filter matrix computation may be performed with modules (e.g., procedures, functions, and so on). The software codes may be stored in memory units (e.g., memory units 642 and 682 in FIG. 6) and executed by processors (e.g., controllers 640 and 680 in FIG. 6). The memory unit may be implemented within the processor or external to the processor, in which case it can be communicatively coupled to the processor via various means as is known in the art.
Headings are included herein for reference and to aid in locating certain sections. These headings are not intended to limit the scope of the concepts described therein under, and these concepts may have applicability in other sections throughout the entire specification.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

## Claims

1. A method of deriving spatial filter matrices in a wireless multiple-input multiple-output (MIMO) communication system, comprising:
determining an initial spatial filter matrix; and
deriving a plurality of spatial filter matrices for a plurality of transmission spans based on the initial spatial filter matrix and a plurality of steering matrices used for the plurality of transmission spans.

2. The method of claim 1, wherein the initial spatial filter matrix is determined based on a channel response matrix for a MIMO channel.

3. The method of claim 2, wherein the spatial filter matrix for each of the plurality of transmission spans is derived based on the initial channel response matrix and a steering matrix used for the transmission span.

4. The method of claim 2, wherein the initial spatial filter matrix is further determined based on a steering matrix for one of the plurality of transmission spans.

5. The method of claim 3, wherein the spatial filter matrix for each of the plurality of transmission spans is derived based on the initial channel response matrix, the steering matrix used to determine the initial spatial filter matrix, and a steering matrix used for the transmission span.

6. The method of claim 1, wherein data is transmitted on orthogonal spatial channels of a MIMO channel, and wherein the initial channel response matrix is determined in accordance with a full channel state information (full-CSI) technique.

7. The method of claim 1, wherein data is transmitted on orthogonal spatial channels of a MIMO channel, and wherein the initial channel response matrix is determined in accordance with a minimum mean square error (MMSE) technique.

8. The method of claim 1, wherein data is transmitted on spatial channels of a MIMO channel, and wherein the initial channel response matrix is determined in accordance with a channel correlation matrix inversion (CCMI) technique.

9. The method of claim 1, wherein data is transmitted on spatial channels of a MIMO channel, and wherein the initial channel response matrix is determined in accordance with a minimum mean square error (MMSE) technique.

10. The method of claim 1, wherein the plurality of steering matrices are used by a transmitting entity to spatially process data to achieve transmit diversity.

11. The method of claim 1, wherein elements of the plurality of steering matrices are members of a set comprised of +1, −1, +j, and −j, where j is a square root of −1.

12. The method of claim 1, wherein the plurality of transmission spans correspond to a plurality of symbol periods.

13. The method of claim 1, wherein the plurality of transmission spans correspond to a plurality of frequency subbands.

14. The method of claim 1, further comprising:
performing spatial processing on symbols received for the plurality of transmission spans with the plurality of spatial filter matrices.

15. An apparatus in a wireless multiple-input multiple-output (MIMO) communication system, comprising:
a processor operative to determine an initial spatial filter matrix and to derive a plurality of spatial filter matrices for a plurality of transmission spans based on the initial spatial filter matrix and a plurality of steering matrices used for the plurality of transmission spans; and
a memory operative to store the plurality of steering matrices.

16. The apparatus of claim 15, wherein the initial spatial filter matrix is determined based on a channel response matrix for a MIMO channel, and wherein the spatial filter matrix for each of the plurality of transmission spans is derived based on the initial channel response matrix and a steering matrix used for the transmission span.

17. The apparatus of claim 15, wherein the initial channel response matrix is determined in accordance with a full channel state information (full-CSI) technique, a minimum mean square error (MMSE) technique, or a channel correlation matrix inversion (CCMI) technique.

18. The apparatus of claim 15, wherein elements of the plurality of steering matrices are members of a set comprised of +1, −1, +j, and −j, where j is a square root of −1.

19. The apparatus of claim 15, further comprising:
a spatial processor operative to perform spatial processing on symbols received for the plurality of transmission spans with the plurality of spatial filter matrices.

20. An apparatus in a wireless multiple-input multiple-output (MIMO) communication system, comprising:
means for determining an initial spatial filter matrix; and
means for deriving a plurality of spatial filter matrices for a plurality of transmission spans based on the initial spatial filter matrix and a plurality of steering matrices used for the plurality of transmission spans.

21. The apparatus of claim 20, wherein the initial spatial filter matrix is determined based on a channel response matrix for a MIMO channel, and wherein the spatial filter matrix for each of the plurality of transmission spans is derived based on the initial channel response matrix and a steering matrix used for the transmission span.

22. The apparatus of claim 20, wherein the initial channel response matrix is determined in accordance with a full channel state information (full-CSI) technique, a minimum mean square error (MMSE) technique, or a channel correlation matrix inversion (CCMI) technique.

23. The apparatus of claim 20, wherein elements of the plurality of steering matrices are members of a set comprised of +1, −1, +j, and −j, where j is a square root of −1.

24. The apparatus of claim 20, further comprising:
means for performing spatial processing on symbols received for the plurality of transmission spans with the plurality of spatial filter matrices.

25. A method of deriving spatial filter matrices in a wireless multiple-input multiple-output (MIMO) communication system, comprising:
deriving a first spatial filter matrix for a first transmission span;
determining a first initial spatial filter matrix for a second transmission span based on the first spatial filter matrix; and
deriving a second spatial filter matrix for the second transmission span based on the first initial spatial filter matrix.

26. The method of claim 25, wherein the first spatial filter matrix is derived based on a channel response matrix obtained for a MIMO channel in the first transmission span and further in accordance with a receiver spatial processing technique.

27. The method of claim 25, wherein the determining the first initial spatial filter matrix comprises
processing the first spatial filter matrix to remove a first steering matrix used for the first transmission span, and wherein the first initial spatial filter matrix is equal to the first spatial filter matrix with the first steering matrix removed.

28. The method of claim 25, wherein the determining the first initial spatial filter matrix comprises:
processing the first spatial filter matrix to remove a first steering matrix used for the first transmission span and to include a second steering matrix used for the second transmission span, and wherein the first initial spatial filter matrix is equal to the first spatial filter matrix with the first steering matrix removed and the second steering matrix included.

29. The method of claim 25, wherein the second spatial filter matrix is derived using an iterative procedure that iteratively performs a set of computations on the first initial spatial filter matrix to obtain a final solution for the second spatial filter matrix.

30. The method of claim 25, further comprising:
determining a second initial spatial filter matrix for a third transmission span based on the second spatial filter matrix; and
deriving a third spatial filter matrix for the third transmission span based on the second initial spatial filter matrix.

31. The method of claim 25, wherein the first and second transmission spans correspond to two different symbol periods.

32. The method of claim 25, wherein the first and second transmission spans correspond to two different frequency subbands.

33. An apparatus in a wireless multiple-input multiple-output (MIMO) communication system, comprising:
a processor operative to derive a first spatial filter matrix for a first transmission span, determine a first initial spatial filter matrix for a second transmission span based on the first spatial filter matrix, and derive a second spatial filter matrix for the second transmission span based on the first initial spatial filter matrix.

34. The apparatus of claim 33, wherein the processor is operative to process the first spatial filter matrix to remove a first steering matrix used for the first transmission span, and wherein the first initial spatial filter matrix is equal to the first spatial filter matrix with the first steering matrix removed.

35. The apparatus of claim 33, wherein the processor is further operative to determine a second initial spatial filter matrix for a third transmission span based on the second spatial filter matrix, and to derive a third spatial filter matrix for the third transmission span based on the second initial spatial filter matrix.

36. An apparatus in a wireless multiple-input multiple-output (MIMO) communication system, comprising:
means for deriving a first spatial filter matrix for a first transmission span;
means for determining a first initial spatial filter matrix for a second transmission span based on the first spatial filter matrix; and
means for deriving a second spatial filter matrix for the second transmission span based on the first initial spatial filter matrix.

37. The apparatus of claim 36, wherein the means for determining the first initial spatial filter matrix comprises
means for processing the first spatial filter matrix to remove a first steering matrix used for the first transmission span, and wherein the first initial spatial filter matrix is equal to the first spatial filter matrix with the first steering matrix removed.

38. The apparatus of claim 36, further comprising:
means for determining a second initial spatial filter matrix for a third transmission span based on the second spatial filter matrix; and
means for deriving a third spatial filter matrix for the third transmission span based on the second initial spatial filter matrix.