# Patent

## Title
Iterative detection and cancellation for wireless communication

## Abstract

Techniques for recovering a desired transmission in the presence of interfering transmissions are described. For iterative detection and cancellation, multiple groups of code channels are formed for a plurality of code channels for at least one sector. Processing is performed for the multiple groups of code channels in multiple iterations. For each iteration, data detection and signal cancellation are performed for the multiple groups of code channels in multiple stages, e.g., in a sequential order starting with the strongest group to the weakest group. Each stage of each iteration may perform data detection, signal reconstruction, and signal cancellation. Each stage of each iteration may also perform equalization, data detection, signal reconstruction, and signal cancellation.

## Background
BACKGROUND
I. Field
The present disclosure relates generally to communication, and more specifically to techniques for recovering transmission in wireless communication.
II. Background
A wireless multiple-access communication system can concurrently communicate with multiple wireless devices, e.g., cellular phones. Examples of such multiple-access systems include Code Division Multiple Access (CDMA) systems, Time Division Multiple Access (TDMA) systems, Frequency Division Multiple Access (FDMA) systems, and Orthogonal FDMA (OFDMA) systems.
A wireless multiple-access system typically includes many base stations that provide communication coverage for a large geographic area. Each base station may transmit data to one or more wireless devices located within its coverage area. A given wireless device may receive a desired transmission from a serving base station as well as interfering transmissions from nearby base stations. These interfering transmissions are intended for other wireless devices located within the coverage areas of the nearby base stations but act as interference to this given wireless device. The interference hinders the wireless device's ability to recover the desired transmission and has a large impact on performance.
There is therefore a need in the art for techniques to recover a desired transmission in the presence of interfering transmissions in a wireless communication system.
SUMMARY
Techniques for recovering a desired transmission in the presence of interfering transmissions are described herein. In an embodiment, multiple groups of code channels (which are also referred to as user groups) are formed for a plurality of code channels for at least one sector. Each group may include all code channels for one sector, a subset of all code channels for one sector, or multiple code channels for multiple sectors. Processing is performed for the multiple groups of code channels in multiple iterations. For each iteration, data detection and signal cancellation are performed for the multiple groups of code channels in multiple stages, e.g., in a sequential order starting with the first group having the strongest received powers and concluding with the last group having the weakest received powers.
In an embodiment, for each stage of each iteration, data detection is performed to obtain a detected signal for a group of code channels being processed by that stage. A signal for the group of code channels is reconstructed based on the detected signal. The reconstructed signal is canceled to obtain an output signal for the stage. In another embodiment, for each stage of each iteration, equalization is performed to obtain an equalized signal for a group of code channels being processed by that stage. Data detection is then performed on the equalized signal to obtain a detected signal for the group of code channels. A signal for the group of code channels is reconstructed based on the detected signal. The reconstructed signal is canceled to obtain an output signal for the stage.
Causal cancellation or both causal and non-causal cancellation may be performed prior to the iterative processing. Various details of the cancellation and the iterative processing are described below. Various aspects and embodiments of the disclosure are also described in further detail below.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a CDMA system.
FIG. 2 shows a block diagram of a base station and a wireless device.
FIG. 3 shows a block diagram of a CDMA modulator.
FIG. 4 shows a processor for iterative detection and cancellation.
FIG. 5 shows a stage that performs detection and cancellation.
FIG. 6 shows a stage that performs equalization, detection and cancellation.
FIG. 7 shows a processor for iterative detection and cancellation with causal and non-causal cancellation.
FIG. 8 shows a code channel gain estimation unit.
FIG. 9 shows a process for performing iterative detection and cancellation.
FIG. 10 shows a detection and cancellation process for one stage.
FIG. 11 shows an equalization, detection and cancellation process for one stage.



DETAILED DESCRIPTION
The techniques described herein may be used for various communication systems such as CDMA, TDMA, FDMA, OFDMA, and Single-Carrier FDMA (SC-FDMA) systems. A CDMA system may implement a radio technology such as cdma2000, Wideband-CDMA (W-CDMA), and so on. cdma2000 covers IS-2000, IS-856, and IS-95 standards. A TDMA system may implement a radio technology such as Global System for Mobile Communications (GSM). These various radio technologies and standards are known in the art. W-CDMA and GSM are described in documents from an organization named “3rd Generation Partnership Project” (3GPP). cdma2000 is described in documents from an organization named “3rd Generation Partnership Project 2” (3GPP2). 3GPP and 3GPP2 documents are publicly available. An OFDMA system utilizes OFDM to transmit symbols in the frequency domain on orthogonal frequency subcarriers. An SC-FDMA system transmits symbols in the time domain on orthogonal frequency subcarriers. For clarity, the techniques are described below for a CDMA system, which may be a cdma2000 system or a W-CDMA system.
FIG. 1 shows a CDMA system 100 with multiple base stations 110 and multiple wireless devices 120. A base station is generally a fixed station that communicates with the wireless devices and may also be called a Node B, an access point, or some other terminology. Each base station 110 provides communication coverage for a particular geographic area 102. The term “cell” can refer to a base station and/or its coverage area depending on the context in which the term is used. To improve system capacity, a base station coverage area may be partitioned into multiple smaller areas, e.g., three smaller areas 104a, 104b, and 104c. Each smaller area is served by a respective base transceiver subsystem (BTS). The term “sector” can refer to a BTS and/or its coverage area depending on the context in which the term is used. For a sectorized cell, the BTSs for all sectors of that cell are typically co-located within the base station for the cell. A system controller 130 couples to base stations 110 and provides coordination and control for these base stations.
The techniques described herein may be used for systems with sectorized cells as well as systems with un-sectorized cells. In the following description, the term “sector” can refer to (1) a BTS and/or its coverage area for a system with sectorized cells and (2) a base station and/or its coverage area for a system with un-sectorized cells. The terms “sector” and “base station” are used interchangeably herein.
Wireless devices 120 are typically dispersed throughout the system, and each wireless device may be stationary or mobile. A wireless device may also be called a mobile station, a user equipment, a terminal, a station, a subscriber unit, or some other terminology. A wireless device may be a cellular phone, a personal digital assistant (PDA), a wireless modem card, a handheld device, a laptop computer, and so on. A wireless device may communicate with zero, one, or multiple base stations on the forward and reverse links at any given moment. The forward link (or downlink) refers to the communication link from the base stations to the terminals, and the reverse link (or uplink) refers to the communication link from the terminals to the base stations. For simplicity, FIG. 1 shows only transmissions on the forward link.
FIG. 2 shows a block diagram of a base station 110 and a wireless device 120, which may be one of the base stations and one of the wireless devices shown in FIG. 1. For simplicity, FIG. 2 shows base station 110 having one transmit antenna and wireless device 120 having one receive antenna. In general, base station 110 and wireless device 120 may each be equipped with any number of antennas. For simplicity, FIG. 2 shows only the processing units for data transmission on the forward link.
At base station 110, a transmit (TX) data processor 210 receives traffic data for the wireless devices being served, processes (e.g., encodes, interleaves, and symbol maps) the traffic data to generate data symbols, and provides the data symbols to a CDMA modulator 220. As used herein, a data symbol is a modulation symbol for data, a pilot symbol is a modulation symbol for pilot, a modulation symbol is a complex value for a point in a signal constellation (e.g., for M-PSK or M-QAM), a symbol is generally a complex value, and pilot is data that is known a priori by both the base stations and the wireless devices. CDMA modulator 220 processes the data symbols and pilot symbols as described below and provides output chips. A transmitter (TMTR) 230 processes (e.g., converts to analog, amplifies, filters, and frequency upconverts) the output chip stream and generates a forward link signal, which is transmitted from an antenna 232.
At wireless device 120, an antenna 252 receives the forward link signals from base station 110 as well as other base stations and provides a received signal. A receiver (RCVR) 254 processes (e.g., filters, amplifies, frequency downconverts, and digitizes) the received signal and provides received samples to a processor 260. Processor 260 may perform iterative detection and cancellation, as described below. Antenna 252 may receive the forward link signal from base station 110 via one or more signal paths, and the received signal may include one or more signal instances (or multipaths) for base station 110. Rake receiver 270 may be used to process all multipaths of interest. Processor 260 or rake receiver 270 provides data symbol estimates, which are estimates of the data symbols sent by base station 110 to wireless device 120. A receive (RX) data processor 280 processes (e.g., symbol demaps, deinterleaves, and decodes) the data symbol estimates and provides decoded data. In general, the processing by processor 260/rake receiver 270 and RX data processor 280 is complementary to the processing by CDMA modulator 220 and TX data processor 210, respectively, at base station 110.
Controllers/processors 240 and 290 direct operation at base station 110 and wireless device 120, respectively. Memories 242 and 292 store data and program codes for base station 110 and wireless device 120, respectively.
In CDMA, multiple orthogonal code channels may be obtained with different orthogonal codes. The code channels may also be referred to as traffic channels, physical channels, data channels, and so on. For example, multiple orthogonal traffic channels are obtained with different Walsh codes in cdma2000, and multiple orthogonal physical channels are obtained with different orthogonal variable spreading factor (OVSF) codes in W-CDMA. The code channels may be used to send different types of data (e.g., traffic data, broadcast data, control data, pilot, and so on) and/or traffic data for different users. Data for the code channels is scaled, combined, and spectrally spread across the entire system bandwidth. The spectral spreading is performed with a spreading code, which is a pseudo-random number (PN) code in cdma2000 and a scrambling code in W-CDMA. In cdma2000, the channelization with Walsh codes is called “covering”, and the spectral spreading is called “spreading”. In W-CDMA, the channelization with OVSF codes is called “spreading”, and the spectral spreading is called “scrambling”. For clarity, cdma2000 terminology (e.g., covering, spreading, Walsh codes, and PN code) is used in the following description.
FIG. 3 shows a block diagram of CDMA modulator 220 within base station 110. For simplicity, the following description assumes that N code channels are available for each sector, and each code channel is assigned a different Walsh code of length N, where N may be equal to 16, 32, 64, 128, 256, or some other value. In general, orthogonal codes of different lengths may be used for the code channels, and N may correspond to the length of the longest orthogonal code. For simplicity, the following description assumes that the N code channels are for N users, and the terms “code channels” and “users” are used interchangeably. In actuality, some code channels are used for overhead, e.g., pilot, control data, broadcast data, etc.
CDMA modulator 220 includes N code channel processors 310a through 310n for the N code channels. Within each code channel processor 310, a multiplier 312 receives and scales the data or pilot symbols for code channel n of sector k with a gain of gk,n and provides scaled symbols. The gain gk,n may be set to zero if code channel n is not used by sector k. A Walsh cover unit 314 channelizes the scaled symbols with a Walsh code wn for code channel n. Unit 314 performs covering by repeating each scaled symbol to generate N replicated symbols and multiplying the N replicated symbols with the N chips of Walsh code wn to generate N data chips for that scaled symbol. A combiner 320 receives and adds the data chips for all N code channels. A PN spreader 322 multiplies the combined data chips with a PN code ck assigned to sector k and generates output chips.
The output chips for sector k in one symbol period may be expressed in matrix form as follows:

sk=CkWGkdk=Akdk,  Eq (1)

where dk is an N×1 vector of data symbols sent on the N code channels of sector k,

Gk is an N×N diagonal matrix of gains for the N code channels of sector k,
W is an N×N Hadamard matrix containing N Walsh codes in N columns,
Ck is an N×N diagonal matrix containing N chips of the PN code for sector k,
Ak is an N×N processing matrix for data vector dk, and
sk is an N×1 vector of output chips for sector k.
For clarity, vectors are denoted with bolded and underlined lower case text (e.g., d), and matrices are denoted with bolded and underlined upper case text (e.g., G). A diagonal matrix contains possible non-zero values along the diagonal and zeros elsewhere.
Vector dk contains N data symbols to be sent simultaneously on N code channels in one symbol period. Matrix Gk contains N gains for the N code channels along the diagonal and zeros elsewhere. The N gains determine the amount of transmit power used for the N code channels. Matrix W contains N Walsh codes for the N code channels in N columns. If the code channels have different Walsh code lengths, then N is equal to the longest Walsh code length for all code channels, and each shorter Walsh code is repeated in matrix W. Since the same Walsh matrix W is used for all sectors, subscript k is not used for W. Matrix Ck contains N PN chips along the diagonal and zeros elsewhere. These PN chips are from the PN code for sector k for one symbol period. Vector sk contains N output chips transmitted by sector k for all N code channels in one symbol period.
Matrix Ak represents all of the processing observed by data vector dk and may be expressed as:

Ak=CkWGk.  Eq (2)

The columns of Ak represent code channels/users, and the rows of Ak represent time.

Wireless device 120 receives the forward link signals from K sectors, which include the serving sector as well as interfering sectors. In general, K may be any value. The received signal for each sector k, without noise, may be expressed as:

xk=HkCkWGkdk=HkAkdk=Hksk,  Eq (3)

where Hk is an (N+Δ)×N channel response matrix for sector k,

xk is an (N+Δ)×1 vector of received samples for sector k, and
Δ is the delay spread of the wireless channel, in units of chips.
Matrix Hk contains channel gains for sector k along the diagonal and zeros elsewhere. Vector xk contains N+Δ received samples for sector k for one symbol period, in the absence of noise. For simplicity, the description herein is for the case in which dk covers one symbol period. dk may also cover multiple symbol periods (e.g., the previous, current, and next symbol periods) to account for intersymbol interference (ISI). In this case, the dimensions of the other matrices would increase correspondingly.
The received samples at wireless device 120 for all K sectors may be expressed as:




  
    
      
        
          
            y
            _
          
          =
          
            
              
                
                  ∑
                  
                    k
                    =
                    1
                  
                  K
                
                ⁢
                
                  
                    x
                    _
                  
                  k
                
              
              +
              
                n
                _
              
            
            =
            
              
                
                  
                    ∑
                    
                      k
                      =
                      1
                    
                    K
                  
                  ⁢
                  
                    
                      
                        H
                        _
                      
                      k
                    
                    ⁢
                    
                      
                        A
                        _
                      
                      k
                    
                    ⁢
                    
                      
                        d
                        _
                      
                      k
                    
                  
                
                +
                
                  n
                  _
                
              
              =
              
                
                  
                    ∑
                    
                      k
                      =
                      1
                    
                    K
                  
                  ⁢
                  
                    
                      
                        H
                        _
                      
                      k
                    
                    ⁢
                    
                      
                        s
                        _
                      
                      k
                    
                  
                
                +
                
                  n
                  _
                
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          4
          )
        
      
    
  




where y is an (N+Δ)×1 vector of received samples at wireless device 120, and

n is an (N+Δ)×1 vector of noise at wireless device 120.
For simplicity, the noise may be assumed to be additive white Gaussian noise (AWGN) with a zero mean vector and a covariance matrix of σn2I, where σn2 is the variance of the noise, and I is the identity matrix with ones along the diagonal and zeros elsewhere.
Wireless device 120 may perform iterative detection and cancellation to recover one or more transmissions from one or more sectors. For iterative detection and cancellation, the users in the K sectors (or the code channels for the K sectors) are processed iteratively multiple times to improve performance. For each iteration, the users in the K sectors are processed in multiple stages. Each stage performs detection and cancellation for one group of users and provides a canceled output to the next stage. The user groups may be sorted based on received powers. For each iteration, the first stage may process the strongest user group, the second stage may process the second strongest user group, and so on, and the last stage may process the weakest user group. Each stage may process a user group using various techniques, as described below.
The users in the K sectors may be arranged into M groups, where M may be any integer value. In general, each user group may include any number of users and any of the users in the K sectors. In an embodiment, each group contains a predetermined number of users, e.g., L users, where L may be any integer value. The user groups may be formed in various manners.
In an embodiment, which is referred to as sector-based processing, each group includes all users in one sector. In this embodiment, M user groups may be formed, with each user group containing L users in one sector, where M=K and L=N. The K sectors may be sorted from strongest to weakest. The first group may contain all users in the strongest sector, the second group may contain all users in the next strongest sector, and so on, and the last group may contain all users in the weakest sector.
In another embodiment, which is referred to as global processing, the received powers for all users in all sectors are determined. The users are then sorted from strongest to weakest and stored in a list. The first group may contain L strongest users in the list, the second group may contain the L next strongest users, and so on, and the last group may contain the L weakest users in the list. In this embodiment, a given group may contain users in the same sector or different sectors.
In yet another embodiment, which is referred to as local processing, each group contains a subset of the users in one sector. The users may be sorted based on their received powers. The first group may contain the L strongest users in the same sector, the second group may contain the L strongest remaining users in the same sector, and so on.
The received signal for user group m, without noise, may be expressed as:

xm=HmCmWmGmdm=HmAmdm=Tmdm,  Eq (5)

where dm is an L×1 data vector for the L users in group m,

Gm is an L×L gain matrix for the users in group m,
Wm is an N×L matrix of Walsh codes for the users in group m,
Cm is an N×N PN matrix for the users in group m
Hm is an (N+Δ)×N channel response matrix for the users in group m,
Tm is an (N+Δ)×L system matrix for data vector dm, and
xm is an (N+Δ)×1 vector of received samples for user group m.
Vector dm and matrices Gm and Wm contain the data symbols, the gains, and the Walsh codes, respectively, for the users in group m. These users may belong in the same sector or different sectors. Matrix Wm may contain duplicate columns if multiple users in group m are assigned code channels with the same Walsh code. Cm contains PN chips for all sectors transmitting to the users in group m. Hm contains complex channel gains for the users in group m. If the users in group m belong in one sector, then Cm and Hm contain PN chips and channel gains for one sector. If the users in group m belong in multiple sectors, then Cm and Hm are block diagonal matrices containing PN chips and channel gains for these multiple sectors, one diagonal channel gain matrix and one diagonal PN matrix for each sector. xm contains the received samples for all users in group m in the absence of noise.
The system matrix for user group m may be given as:

Tm=HmCmWmGm.  Eq (6)

System matrix Tm represents all of the processing as well as the channel response observed by data vector dm. The height of Tm is related to time (in number of chips), and the width of Tm is determined by the number of users in group m. A single system matrix T may be defined for all M·L users in all K sectors. However, the processing for this single large system matrix T would be computationally intensive.

The received samples at wireless device 120 for the M user groups may be expressed as:




  
    
      
        
          y
          _
        
        =
        
          
            
              
                ∑
                
                  m
                  =
                  1
                
                M
              
              ⁢
              
                
                  x
                  _
                
                m
              
            
            +
            
              n
              _
            
          
          =
          
            
              
                ∑
                
                  m
                  =
                  1
                
                M
              
              ⁢
              
                
                  
                    T
                    _
                  
                  m
                
                ⁢
                
                  
                    d
                    _
                  
                  m
                
              
            
            +
            
              
                n
                _
              
              .
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          7
          )
        
      
    
  




For the first iteration, each stage performs data detection based on its input signal and obtains a detected signal for the user group being processed by that stage. Each stage also reconstructs the signal for its user group and subtracts the reconstructed signal from its input signal to obtain an output signal for this stage, which is an input signal for the next stage. As an example, for the case in which M=3, the received signal may be expressed as:

y=T1d1+T2d2+T3d3+n.  Eq (8)

The received signal y is the input signal y1,0 for the first stage of the first iteration, or y1,0=y.

In the first iteration, the processing by the three stages for three user groups may be expressed as:

y1,1=y1,0−T1{circumflex over (d)}1,1=e1,1+T2d2+T3d3+n, 

y1,2=y1,1−T2{circumflex over (d)}1,2=e1,1+e1,2+T3d3+n, and

y1,3=y1,2−T3{circumflex over (d)}1,3=e1,1+e1,2+e1,3+n,  Eq (9)

where {circumflex over (d)}1,1 to {circumflex over (d)}1,3 are detected signals for user groups 1 to 3, respectively,

y1,1 to y1,3 are output signals from stages 1 to 3, respectively, and
e1,1 to e1,3 are residual errors for user groups 1 to 3, respectively.
Each line in equation set (9) represents the processing by one stage. Each stage m performs data detection to obtain {circumflex over (d)}1,m, then reconstructs the signal for user group m as Tm {circumflex over (d)}1,m, and then subtracts the reconstructed signal from the input signal y1,m−1 to obtain the output signal y1,m. The residual error e1,m is the difference between the actual signal and the reconstructed signal, or e1,m=Tm d1,m−Tm {circumflex over (d)}1,m.
A second iteration may then be performed for the three user groups. The output signal y1,3 from the last stage of the first iteration is used as the input signal y2,0 for the first stage of the second iteration, or y2,0=y1,3.
In the second iteration, the first stage may add the reconstructed signal for the first user group back to the input signal, as follows:

r2,1=y2,0+T1{circumflex over (d)}1,1=T1d1+e1,2+e1,3+n  Eq (10)

where r2,1 is a reconstituted signal containing the signal for the first user group as well as residual errors for the other user groups and noise. The first user group's signal in r2,1 will likely have higher quality than in y because r2,1 contains residual errors instead of actual signals for the other user groups.

The first stage performs data detection on the reconstituted signal r2,1 to obtain an improved detected signal for the first user group. The signal due to the first user group may be reconstructed and canceled from the reconstituted signal to obtain the output signal for the first stage, which is an input signal for the next stage.
For the second iteration, the processing by the three stages for the three user groups may be expressed as:

r2,1y2,0+T1{circumflex over (d)}1,1=T1d1+e1,2+e1,3+n and y2,1=r2,1−T1{circumflex over (d)}2,1=e2,1+e1,2+e1,3+n, 

r2,2=y2,1+T2{circumflex over (d)}1,2=e2,1+T2d2+e1,3+n and y2,2=r2,2−T2{circumflex over (d)}2,2=e2,1+e2,2+e1,3+n, 

r2,3=y2,2+T3{circumflex over (d)}1,3=e2,1+e2,2+T3d3+n and y2,3=r2,3−T3{circumflex over (d)}2,3=e2,1+e2,2+e2,3+n,  Eq (11)

where r2,1 to r2,3 are reconstituted signals for stages 1 to 3, respectively,

{circumflex over (d)}2,1 to {circumflex over (d)}2,3 are detected signals for user groups 1 to 3, respectively,
y2,1 to y2,3 are output signals from stages 1 to 3, respectively, and
e2,1 to e2,3 are residual errors for user groups 1 to 3, respectively.
Each line in equation set (11) represents the processing by one stage. Each stage m generates the reconstituted signal r2,m by adding back the reconstructed signal Tm {circumflex over (d)}1,m from the prior iteration, then performs data detection to obtain {circumflex over (d)}2,m, then reconstructs the signal for user group m as Tm {circumflex over (d)}2,m, and then subtracts the reconstructed signal from the reconstituted signal r2,m to obtain the output signal y2,m. The residual error is the difference between the actual signal and the newly reconstructed signal, or e2,m=Tmd1,m−Tm{circumflex over (d)}2,m. The detected signal {circumflex over (d)}2,m in the second iteration should be more accurate than the detected signal {circumflex over (d)}1,m in the first iteration. Hence, the residual error e2,m in the second iteration should be smaller than the residual error e1,m in the first iteration.
In general, any number of iterations may be performed for the M user groups. The output signal from the last stage of a given iteration i may be used as the input signal for the first stage of the next iteration i+1. For each stage in the second iteration and later, the reconstructed signal for a user group is added back to the input signal for that stage to obtain the reconstituted signal. For each stage in each iteration, the signal for a user group is detected, reconstructed, and subtracted from the input/reconstituted signal for that stage to obtain the output signal for the next stage. The signal quality generally improves after each iteration. Hence, the detected signals are generally more accurate after each iteration.
FIG. 4 shows a block diagram of a processor 260a for iterative detection and cancellation, which is an embodiment of processor 260 in FIG. 2. Processor 260a includes M processing stages 420a through 420M for M user groups. For sector-based processing, each stage performs processing for one sector. M may be equal to K if all sectors are processed or may be less than K if fewer sectors are processed, e.g., if weak sectors are skipped.
A multiplexer 410 is provided with the received signal y on one input and the output signal yi,M from the last stage 420M on the other input, where subscript i denotes the iteration number. Multiplexer 410 provides the received signal y for the first iteration and provides the output signal yi,M for each subsequent iteration.
The M stages 420a through 420M are coupled in cascaded, with the output of each stage being provided to the input of the next stage. Each stage 420 derives the reconstituted signal ri,m based on its input signal yi,m−1, performs data detection on the reconstituted signal to obtain the detected signal {circumflex over (d)}i,m for user group m being processed by that stage. Each stage 420 also reconstructs the signal for user group m and subtracts the reconstructed signal {circumflex over (x)}i,m from the reconstituted signal to obtain the output signal yi,m for the stage.
A memory 430 stores intermediate results from the M stages, e.g., the detected signals, reconstructed signals, etc. The intermediate results may be stored in the forms of chips, symbols, etc. Memory 430 provides the intermediate results to the M stages as needed.
Each stage may perform processing for its user group using various techniques such as a detection and cancellation technique and an equalization, detection and cancellation technique. These techniques are described below.
FIG. 5 shows an embodiment of a stage 420x that performs detection and cancellation. Stage 420x may be used for each of the M stages 420a through 420M in FIG. 4. In the following description, stage 420x is referred to as stage m and performs processing for user group m.
Within stage m, a summer 510 receives the input signal yi,m−1 for stage m and the reconstructed signal {circumflex over (x)}i−1,m for stage m from the prior iteration. Summer 510 sums the two signals and provides the reconstituted signal ri,m for stage m, which may be expressed as:




  
    
      
        
          
            
              
                r
                _
              
              
                i
                ,
                m
              
            
            =
            
              
                
                  y
                  _
                
                
                  i
                  ,
                  
                    m
                    -
                    1
                  
                
              
              +
              
                
                  
                    x
                    _
                  
                  ^
                
                
                  
                    i
                    -
                    1
                  
                  ,
                  m
                
              
            
          
          ,
          
            

          
          ⁢
          where
        
        ⁢
        
          

        
        ⁢
        
          
            
              
                x
                _
              
              ^
            
            
              
                i
                -
                1
              
              ,
              m
            
          
          =
          
            {
            
              
                
                  
                    0
                  
                  
                    
                      
                        for
                        ⁢
                        
                          
                        
                        ⁢
                        i
                      
                      =
                      1
                    
                  
                
                
                  
                    
                      
                        
                          T
                          _
                        
                        m
                      
                      ⁢
                      
                        
                          
                            d
                            _
                          
                          ^
                        
                        
                          
                            i
                            -
                            1
                          
                          ,
                          m
                        
                      
                    
                  
                  
                    
                      
                        for
                        ⁢
                        
                          
                        
                        ⁢
                        i
                      
                      >
                      1
                    
                  
                
              
              .
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          12
          )
        
      
    
  




For the first iteration, the input signal yi,m−1 contains the signal for user group m, {circumflex over (x)}i−1,m is equal to zero, and the reconstituted signal is equal to the input signal, or ri,m=yi,m−1. For each subsequent iteration, the input signal contains the residual error ei−1,m for user group m, and the reconstructed signal {circumflex over (x)}i−1,m for user group m from the prior iteration i−1 is added back to the input signal yi,m−1 to obtain the reconstituted signal ri,m, which then contains the signal for user group m.
A block 514 performs data detection on the reconstituted signal ri,m and provides the detected signal {circumflex over (d)}i,m for user group m. Data detection may be performed in various manners, as described below. A block 516 reconstructs the signal for user group m based on the detected signal and provides the reconstructed signal {circumflex over (x)}i,m for user group m, which may be expressed as:

{circumflex over (x)}i,m=Tm{circumflex over (d)}i,m.  Eq (13)

A summer 518 subtracts the reconstructed signal from the reconstituted signal and provides the output signal yi,m for stage m, which may be expressed as:

yi,m=ri,m−{circumflex over (x)}i,m.  Eq (14)

The reconstituted signal ri,m for stage m in iteration i may be expressed as:

ri,m=Tmdm+vi,m,  Eq (15)

where vi,m is the total noise and interference for stage m in iteration i.

The total noise and interference vi,m may be expressed as:




  
    
      
        
          
            
              v
              _
            
            
              1
              ,
              m
            
          
          =
          
            
              
                ∑
                
                  j
                  =
                  1
                
                
                  m
                  -
                  1
                
              
              ⁢
              
                
                  e
                  _
                
                
                  1
                  ,
                  j
                
              
            
            +
            
              
                ∑
                
                  j
                  =
                  
                    m
                    +
                    1
                  
                
                M
              
              ⁢
              
                
                  
                    T
                    _
                  
                  j
                
                ⁢
                
                  
                    d
                    _
                  
                  j
                
              
            
            +
            
              
                n
                _
              
              ⁢
              
                
              
              ⁢
              for
              ⁢
              
                
              
              ⁢
              the
              ⁢
              
                
              
              ⁢
              first
              ⁢
              
                
              
              ⁢
              iteration
            
          
        
        ,
        and
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          16
          )
        
      
    
  
  
    
      
        
          
            v
            _
          
          
            i
            ,
            m
          
        
        =
        
          
            
              
                ∑
                
                  j
                  =
                  1
                
                
                  m
                  -
                  1
                
              
              ⁢
              
                
                  e
                  _
                
                
                  i
                  ,
                  j
                
              
            
            +
            
              
                ∑
                
                  j
                  =
                  
                    m
                    +
                    1
                  
                
                M
              
              ⁢
              
                
                  e
                  _
                
                
                  
                    i
                    -
                    1
                  
                  ,
                  j
                
              
            
            +
            
              n
              ⁢
              
                
              
              ⁢
              for
              ⁢
              
                
              
              ⁢
              iteration
              ⁢
              
                
              
              ⁢
              i
            
          
          >
          1.
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          17
          )
        
      
    
  




For the first iteration, vi,m contains residual errors e1,j for user groups 1 through m−1 that have been canceled in prior stages as well as signals Tj dj for user groups m+1 through M to be processed in later stages. For each subsequent iteration, vi,m contains residual errors ei,j for user groups 1 through m−1 from the current iteration i as well as residual errors ei−1,j for user groups m+1 through M from the prior iteration i−1.
Data detection may be performed on the reconstituted signal using linear minimum mean square error (MMSE), least squares (LS), or some other data detection technique. Data detection may be performed for user group m based on the linear MMSE technique, as follows:

{circumflex over (d)}i,m=Rdd,mTmH(Tm Rdd,m TmH+Rvv,i,m)−1ri,m,  Eq (18)

where Rdd,m=E{dm dmH} is the covariance of data signal dm for user group m,

Rvv,i,m=E{vi,m vi,mH} is the covariance of total noise and interference vector vi,m,
E{ } denotes an expectation operation, and
“H” denotes a conjugate transpose.
The data symbols in dm may be assumed to be uncorrelated so that Rdd,m=I. The total noise and interference vi,m may be assumed to be AWGN so that Rvv,i,m=σv,i,m2I, where σv,i,m2 is the variance of the total noise and interference for user group m in iteration i. σv,i,m2 may be estimated as described below.
Equation (18) may then be expressed as:




  
    
      
        
          
            
              
                
                  
                    d
                    _
                  
                  ^
                
                
                  i
                  ,
                  m
                
              
              =
              
                
                  
                    (
                    
                      
                        
                          
                            T
                            _
                          
                          m
                          H
                        
                        ⁢
                        
                          
                            T
                            _
                          
                          m
                        
                      
                      +
                      
                        
                          σ
                          
                            v
                            ,
                            i
                            ,
                            m
                          
                          2
                        
                        ⁢
                        
                          I
                          _
                        
                      
                    
                    )
                  
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    T
                    _
                  
                  m
                  H
                
                ⁢
                
                  
                    r
                    _
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    Z
                    _
                  
                  
                    mmse
                    ,
                    i
                    ,
                    m
                  
                
                ⁢
                
                  
                    r
                    _
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          19
          )
        
      
    
  




where Zmmse,i,m=(TmHTm+σv,i,m2I)−1 TmH is an MMSE data detection filter for user group m in iteration i. Equation (19) is derived from equation (18) using the matrix inverse lemma and the assumptions for Rdd,m and Rvv,i,m described above.

Data detection may also be performed for user group m based on the least squares technique, as follows:




  
    
      
        
          
            
              
                
                  
                    d
                    _
                  
                  ^
                
                
                  i
                  ,
                  m
                
              
              =
              
                
                  
                    (
                    
                      
                        
                          T
                          _
                        
                        m
                        H
                      
                      ⁢
                      
                        
                          T
                          _
                        
                        m
                      
                    
                    )
                  
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    T
                    _
                  
                  m
                  H
                
                ⁢
                
                  
                    r
                    _
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    Z
                    _
                  
                  
                    ls
                    ,
                    i
                    ,
                    m
                  
                
                ⁢
                
                  
                    r
                    _
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          20
          )
        
      
    
  




where Zls,i,m=(TmH Tm)−1 TmH is a least squares data detection filter for user group m in iteration i.

FIG. 6 shows an embodiment of a stage 420y that performs equalization, detection and cancellation. Stage 420y may also be used for each of the M stages 420a through 420M in FIG. 4. In the following description, stage 420y is referred to as stage m and performs processing for user group m.
Within stage m, a summer 610 receives the input signal yi,m−1 and the reconstructed signal {circumflex over (x)}i−1,m, sums the two signals, and provides the reconstituted signal ri,m for stage m. A block 612 performs equalization on the reconstituted signal ri,m and provides an equalized signal ŝi,m for user group m. A block 614 performs data detection on the equalized signal and provides the detected signal {circumflex over (d)}i,m for user group m. A block 616 reconstructs the signal for user group m based on the detected signal and provides the reconstructed signal {circumflex over (x)}i,m for user group m. A summer 618 subtracts the reconstructed signal from the reconstituted signal and provides the output signal yi,m for stage m.
Equalization may be performed using linear MMSE, least squares, or some other equalization technique. A equalizer matrix may be derived based on the linear MMSE technique, as follows:

Mmmse,i,m=[HmHHm+σv,i,m2I]−1HmH,  Eq (21)

where Mmmse,i,m is an MMSE equalizer matrix for user group m in iteration i. Hm may be estimated based on the pilot received from the sector(s) for user group m. Mmmse,i,m is a function of σv,i,m2 and may be different for different iterations.

A equalizer matrix may also be derived based on the least squares technique, as follows:

Mls,m=[HmHHm]−1HmH,  Eq (22)

where Mls,m is a least squares equalizer matrix for user group m.

Equalization may be performed for user group m as follows:




  
    
      
        
          
            
              
                
                  
                    s
                    _
                  
                  ^
                
                
                  i
                  ,
                  m
                
              
              =
              
                
                  
                    M
                    _
                  
                  
                    i
                    ,
                    m
                  
                
                ⁢
                
                  
                    r
                    _
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    M
                    _
                  
                  
                    i
                    ,
                    m
                  
                
                ⁡
                
                  (
                  
                    
                      
                        
                          T
                          _
                        
                        m
                      
                      ⁢
                      
                        
                          d
                          _
                        
                        m
                      
                    
                    +
                    
                      
                        v
                        _
                      
                      
                        i
                        ,
                        m
                      
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
                
                  
                    
                      M
                      _
                    
                    
                      i
                      ,
                      m
                    
                  
                  ⁢
                  
                    
                      T
                      _
                    
                    m
                  
                  ⁢
                  
                    
                      d
                      _
                    
                    m
                  
                
                +
                
                  
                    V
                    _
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    
                      M
                      _
                    
                    
                      i
                      ,
                      m
                    
                  
                  ⁢
                  
                    
                      H
                      _
                    
                    m
                  
                  ⁢
                  
                    
                      A
                      _
                    
                    m
                  
                  ⁢
                  
                    
                      d
                      _
                    
                    m
                  
                
                +
                
                  
                    V
                    _
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          23
          )
        
      
    
  




where Mm is an equalizer matrix for user group m,

Vi,m is a matrix of total noise and interference for user group m, and
ŝi,m is an equalized signal for user group m in iteration i.
Equalizer matrix Mi,m may be derived using the linear MMSE or least squares technique. ŝi,m is an estimate of data signal sm for user group m.
Equation (23) may be approximated as follows:

ŝi,m≈Amdm+Vi,m.  Eq (24)

The total noise and interference Vi,m may be expressed as:

Vi,m=Mi,mvi,m.  Eq (25)

Data detection may be performed on the equalized signal using linear MMSE, least squares, or some other data detection technique. Data detection may be performed for user group m based on the linear MMSE technique, as follows:




  
    
      
        
          
            
              
                
                  
                    d
                    _
                  
                  ^
                
                
                  i
                  ,
                  m
                
              
              =
              
                
                  
                    (
                    
                      
                        
                          
                            A
                            _
                          
                          m
                          H
                        
                        ⁢
                        
                          
                            A
                            _
                          
                          m
                        
                      
                      +
                      
                        
                          R
                          _
                        
                        
                          vv
                          ,
                          i
                          ,
                          m
                        
                      
                    
                    )
                  
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    A
                    _
                  
                  m
                  H
                
                ⁢
                
                  
                    
                      s
                      _
                    
                    ^
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    (
                    
                      
                        
                          
                            G
                            _
                          
                          m
                          H
                        
                        ⁢
                        
                          
                            W
                            _
                          
                          H
                        
                        ⁢
                        
                          
                            C
                            _
                          
                          m
                          H
                        
                        ⁢
                        
                          
                            C
                            _
                          
                          m
                        
                        ⁢
                        
                          W
                          _
                        
                        ⁢
                        
                          
                            G
                            _
                          
                          s
                        
                      
                      +
                      
                        
                          R
                          _
                        
                        
                          vv
                          ,
                          i
                          ,
                          m
                        
                      
                    
                    )
                  
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    A
                    _
                  
                  m
                  H
                
                ⁢
                
                  
                    
                      s
                      _
                    
                    ^
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    (
                    
                      
                        N
                        ⁢
                        
                          
                            G
                            _
                          
                          m
                          H
                        
                        ⁢
                        
                          
                            G
                            _
                          
                          m
                        
                      
                      +
                      
                        
                          R
                          _
                        
                        
                          vv
                          ,
                          i
                          ,
                          m
                        
                      
                    
                    )
                  
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    A
                    _
                  
                  m
                  H
                
                ⁢
                
                  
                    
                      s
                      _
                    
                    ^
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          26
          )
        
      
    
  




where Rvv,i,m is a covariance matrix for the total noise and interference Vi,m. Gain matrix Gm and covariance matrix Rvv,i,m may be determined as described below.

The total noise and interference matrix Vi,m may be approximated as white. Data detection may then be performed as follows:




  
    
      
        
          
            
              
                
                  
                    d
                    _
                  
                  ^
                
                
                  i
                  ,
                  m
                
              
              =
              
                
                  
                    (
                    
                      
                        N
                        ⁢
                        
                          
                            G
                            _
                          
                          m
                          H
                        
                        ⁢
                        
                          
                            G
                            _
                          
                          m
                        
                      
                      +
                      
                        diag
                        ⁡
                        
                          (
                          
                            
                              R
                              _
                            
                            
                              vv
                              ,
                              i
                              ,
                              m
                            
                          
                          )
                        
                      
                    
                    )
                  
                  
                    -
                    1
                  
                
                ⁢
                
                  
                    A
                    _
                  
                  m
                  H
                
                ⁢
                
                  
                    
                      s
                      _
                    
                    ^
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
        
          
            
              =
              
                
                  
                    Z
                    _
                  
                  
                    mmse
                    ,
                    i
                    ,
                    m
                  
                
                ⁢
                
                  
                    
                      s
                      _
                    
                    ^
                  
                  
                    i
                    ,
                    m
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          27
          )
        
      
    
  




where Zmmse,i,m=(NGmHGm+diag(Rvv,i,m))−1 AmH. Zmmse,i,m is an MMSE data detection filter for user group m with the approximation of Vi,m as white. Since (NGmHGm+diag(Rvv,i,m)) is a diagonal matrix, the inverse operation becomes simple scalar operation, and the matrix inversion in equation (26) can be avoided.

Data detection may also be performed for user group m based on the least squares technique, as follows:

{circumflex over (d)}i,m=AmHŝi,m,  Eq (28)

where AmH is a least squares data detection filter for user group m.

In the embodiments shown in FIGS. 5 and 6, the reconstructed signal {circumflex over (x)}i,m is stored in memory 430 and provided to summer 510 or 610 in the next iteration. In another embodiment, the detected signal {circumflex over (d)}i,m is stored in memory 430 and used to generate the reconstructed signal {circumflex over (x)}i,m in the next iteration. Storing the detected signal {circumflex over (d)}i,m utilizes less memory. However, storing the reconstructed signal {circumflex over (x)}i,m avoids the regeneration of {circumflex over (x)}i,m from {circumflex over (d)}i,m.
The forward link signal transmitted by a given sector may reach the wireless device via one or more signal paths. These signal paths may include a line-of-sight path and/or reflected paths, which are created by reflections of radio waves in the environment. The multipath effect causes ISI, which is a phenomenon in which a data symbol in a received signal is smeared over time and acts as distortion to one or more other symbols in the received signal
Data for the previous, current, and next symbol periods may be processed concurrently in order to capture and remove ISI. The received signal may then be expressed as:

y=[TpTcTn][dpdcdn]T+n,  Eq (29)

where dp, dc and dn are data vectors for three symbol periods, and

Tp, Tc and Tn are system matrices for three symbol periods.
Data vector dc includes data symbols for all users in all K sectors for the current symbol period (“c”). System matrix Tc covers all users in all K sectors for the current symbol period. Similarly, dp and Tp cover all users for the prior symbol period (“p”), and dn and Tn cover all users for the next symbol period (“n”).
Equation (29) indicates that a three fold increase in the column space of the system matrix and a corresponding increase in computational complexity will result from considering the previous, current, and next symbol periods concurrently.
In an embodiment, which is referred to as causal cancellation, the detected symbols for the previous symbol period are subtracted from the received signal prior to the iterative processing. This embodiment mitigates computational increase while accounting for ISI due to the previous symbol. The detected symbols for symbol period t may be stored and used in the next symbol period t+1. In symbol period t+1, the stored detected symbols for the previous symbol period t may be subtracted from the input signal, as follows:

yc=y−Tp{circumflex over (d)}p≈[TcTn][dcdn]T+nc,  Eq (30)

where {circumflex over (d)}p is a vector of detected symbols for all users for the prior symbol period, and

yc is a received signal with causal cancellation.
In another embodiment, which is referred to as causal and non-causal cancellation, the detected symbols for the previous, current, and next symbol periods are subtracted from the received signal prior to the iterative processing. One round of processing may be performed, e.g., on the received signal yc with causal cancellation, to obtain detected symbols for all symbol periods of interest, which may be stored in memory. A received signal with causal and non-causal cancellation may be obtained as:

ycn=y−Tp{circumflex over ({circumflex over (d)}p−Tn{circumflex over (d)}n≈Tcdc+vi,m,  Eq (31)

where {circumflex over ({circumflex over (d)}p is a vector of detected symbols for all users for the prior symbol period, and

ycn is a received signal with causal and non-causal cancellation.
In the second round of processing, for the current symbol period, the detected symbols {circumflex over ({circumflex over (d)}p obtained in this round for the prior symbol period and the detected symbols {circumflex over (d)}c and {circumflex over (d)}n obtained in the first round for the current and next symbol period are subtracted from the received signal to obtain ycn. Iterative processing is then performed on ycn to obtain detected symbols {circumflex over ({circumflex over (d)}c for the current symbol period. {circumflex over ({circumflex over (d)}c is subtracted from the received signal in the next symbol period.
FIG. 7 shows a block diagram of a processor 260b for iterative detection and cancellation with causal and non-causal cancellation, which is another embodiment of processor 260 in FIG. 2. For the first round, a block 702 performs causal cancellation on the received signal y based on the reconstructed signal {circumflex over (x)}p=Tp {circumflex over (d)}p from a memory 730, e.g., as shown in equation (30), and provides the received signal yc with causal cancellation. For the second round, a block 704 performs causal and non-causal cancellation on the received signal based on the reconstructed signals {circumflex over ({circumflex over (x)}p=Tp {circumflex over ({circumflex over (d)}p and {circumflex over (x)}n=Tn {circumflex over (d)}n from memory 730, e.g., as shown in equation (31), and provides the received signal ycn with causal and non-causal cancellation.
A multiplexer 710 is provided with the received signal yc from block 702 on one input, the received signal ycn from block 704 on another input, and the output signal yi,M from the last stage 720M on a third input. Multiplexer 710 provides the received signal yc for the first iteration of the first round, the received signal ycn for the first iteration of the second round, and the output signal yi,M for each subsequent iteration. Processing stages 720a through 720M process the signal from multiplexer 710 as described above. For example, each processing stage 720 may be implemented as shown in FIG. 5 or 6. Memory 730 stores the detected signal and/or the reconstructed signal and provides these signals to block 702, block 704, and stages 720 as appropriate.
In each iteration, the M user groups may be processed sequentially from the strongest user group to the weakest user group. This may improve detection performance for each user group since the signals from stronger user groups (if any) have been canceled. The user groups may also be processed in other orders. In general, the sequential processing of the user groups may result in the signal quality progressively improving for each subsequently processed user group, since the interfering signals from earlier processed user groups have been removed.
The wireless device may desire to recover a signal from a single sector. In an embodiment, the user group containing the desired signal is processed in the last stage after canceling the signals from other user groups. In another embodiment, the M user groups are processed from strongest to weakest, as described above. If the desired user group is not the last user group that is processed, then reconstructed signal for the desired user group may be added back to the output signal from the last stage to obtain the reconstituted signal for the desired user group. The reconstituted signal may then be processed to detect the desired signal.
The wireless device may desire to recover signals from multiple sectors, e.g., for soft handoff. In an embodiment, the user groups containing the desired signals from these sectors are processed in the last few stages after canceling the signals from other user groups. In another embodiment, the M user groups are processed from strongest to weakest, as described above. For each desired user group, the reconstructed signal for that user group may be added back to the output signal from the last stage, and the reconstituted signal may be processed to recover the signal from that user group.
The successive processing of the M user groups, one user group at a time, may have certain advantages. First, the size of the vectors and matrices for one user group may be much smaller than the size of the vectors and matrices for all users in full multi-user detection (MUD). Second, the successive processing of the M user groups may result in the signal quality progressively improving for each subsequently processed user group. The iterative processing of the M user groups further improves detection performance.
In equation (19), the variance of the total noise and interference, σv,i,m2, for user group m in iteration i may be estimated as follows. The total received power at the wireless device, Ptotal—rx, may be expressed as:




  
    
      
        
          
            P
            total_rx
          
          =
          
            
              
                
                  ∑
                  
                    k
                    =
                    1
                  
                  K
                
                ⁢
                
                  (
                  
                    
                      ∑
                      i
                    
                    ⁢
                    
                      
                        h
                        
                          k
                          ,
                          i
                        
                        2
                      
                      ·
                      
                        
                          ∑
                          
                            n
                            =
                            1
                          
                          N
                        
                        ⁢
                        
                          g
                          
                            k
                            ,
                            n
                          
                          2
                        
                      
                    
                  
                  )
                
              
              +
              
                σ
                n
                2
              
            
            =
            
              
                
                  ∑
                  
                    k
                    =
                    1
                  
                  K
                
                ⁢
                
                  
                    ∑
                    i
                  
                  ⁢
                  
                    
                      ∑
                      
                        n
                        =
                        1
                      
                      N
                    
                    ⁢
                    
                      
                        h
                        
                          k
                          ,
                          i
                        
                        2
                      
                      ·
                      
                        g
                        
                          k
                          ,
                          n
                        
                        2
                      
                    
                  
                
              
              +
              
                σ
                n
                2
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          32
          )
        
      
    
  




where hk,i is a channel gain for signal path i of sector k. Each sector may be associated with any of number of signal paths at the wireless device, depending on the wireless environment.

For each stage m in the first iteration, the powers of the user groups processed by the prior stages have been removed (albeit not perfectly), and the power of the user group being processed by stage m is not part of the noise power. The total received power may be expressed as:

Ptotal—rx=Pprior—stages+Pcurrent—stage+Plater—stages,  Eq (33)

where Pprior—stages, Pcurrent—stage, and Plater—stages are the powers of the user groups processed by the prior, current, and later stages, respectively. Plater—stages may also include the powers of groups that are not processed. The powers may be different between the first iteration and each subsequent iteration.

Ideally, σv,1,m2=Plater—stages+σn2 since the powers of the prior stages have been removed. σv,1,m2 may then be expressed as:




  
    
      
        
          
            σ
            
              v
              ,
              1
              ,
              m
            
            2
          
          =
          
            
              
                ∑
                
                  k
                  ∈
                  
                    K
                    m
                  
                
                
                  
                
              
              ⁢
              
                
                  ∑
                  
                    n
                    ∈
                    
                      N
                      m
                    
                  
                
                ⁢
                
                  
                    h
                    k
                    2
                  
                  ·
                  
                    g
                    
                      k
                      ,
                      n
                    
                    2
                  
                
              
            
            +
            
              σ
              n
              2
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          34
          )
        
      
    
  




where Km and Nm include all user groups processed by later stages.

For each subsequent iteration, σv,i,m2 is determined by the residual errors ei,j for user groups 1 through m−1 from the current iteration i as well as the residual errors ei−1,j for user groups m+1 through M from the prior iteration i−1.
The thermal noise σn2 in equation (34) may be expressed as:




  
    
      
        
          σ
          n
          2
        
        =
        
          
            
              P
              total_rx
            
            -
            
              
                ∑
                
                  k
                  =
                  1
                
                K
              
              ⁢
              
                
                  ∑
                  
                    n
                    =
                    1
                  
                  N
                
                ⁢
                
                  
                    h
                    k
                    2
                  
                  ·
                  
                    g
                    
                      k
                      ,
                      n
                    
                    2
                  
                
              
            
          
          =
          
            
              E
              ⁢
              
                {
                
                  
                    y
                    _
                  
                  ⁢
                  
                    
                      y
                      _
                    
                    H
                  
                
                }
              
            
            -
            
              
                ∑
                
                  k
                  =
                  1
                
                K
              
              ⁢
              
                
                  ∑
                  
                    n
                    =
                    1
                  
                  N
                
                ⁢
                
                  
                    h
                    k
                    2
                  
                  ·
                  
                    
                      g
                      
                        k
                        ,
                        n
                      
                      2
                    
                    .
                  
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          35
          )
        
      
    
  




In equations (26) and (27), covariance matrix Rvv,i,m may be estimated as follows. Equation (23) may be rewritten as:

ŝi,m=Bi,m+Vi,m,  Eq (36)

where Bi,m=Mi,mHmAmdm. The covariances of ŝi,m, Bi,m and Vi,m may be expressed as:

Rss,i,m=Rbb,i,m+Rvv,i,m.  Eq (37)

The covariance of ŝi,m may be expressed as:

Rss,i,m=E{ŝi,mŝi,mH},  Eq (38)

Rss,i,m may be estimated by computing the outer product of ŝi,m and averaging over multiple symbol periods.

The covariance of Bi,m may be estimated as:

Rbb,i,m=E{Mi,mHmAmAmHHmHMi,mH}.  Eq (39)

Am, Hm and Mi,m may be estimated and used to derive Rbb,i,m.

The covariance of Vi,m may then be expressed as:

Rvv,i,m=E{Vi,mVi,mH}=Rss,i,m−Rbb,i,m.  Eq (40)

The gain of each user (or code channel) in each sector may be estimated and used to form the gain matrix Gm for each user group. To estimate the gains, received signal y may be processed as follows:




  
    
      
        
          
            
              
                
                  q
                  _
                
                k
              
              =
              
                
                  
                    W
                    _
                  
                  H
                
                ⁢
                
                  
                    C
                    _
                  
                  k
                  H
                
                ⁢
                
                  
                    H
                    _
                  
                  k
                  H
                
                ⁢
                
                  y
                  _
                
              
            
          
        
        
          
            
              =
              
                
                  
                    W
                    _
                  
                  H
                
                ⁢
                
                  
                    C
                    _
                  
                  k
                  H
                
                ⁢
                
                  
                    
                      H
                      _
                    
                    k
                    H
                  
                  (
                  
                    
                      
                        
                          H
                          _
                        
                        k
                      
                      ⁢
                      
                        
                          C
                          _
                        
                        k
                      
                      ⁢
                      
                        W
                        _
                      
                      ⁢
                      
                        
                          G
                          _
                        
                        k
                      
                      ⁢
                      
                        
                          d
                          _
                        
                        k
                      
                    
                    +
                    
                      
                        ∑
                        
                          
                            j
                            =
                            1
                          
                          ,
                          
                            j
                            ≠
                            k
                          
                        
                        K
                      
                      ⁢
                      
                        
                          
                            H
                            _
                          
                          j
                        
                        ⁢
                        
                          
                            C
                            _
                          
                          j
                        
                        ⁢
                        
                          W
                          _
                        
                        ⁢
                        
                          
                            G
                            _
                          
                          j
                        
                        ⁢
                        
                          
                            d
                            _
                          
                          j
                        
                      
                    
                    +
                    
                      n
                      _
                    
                  
                  )
                
              
            
          
        
        
          
            
              =
              
                
                  
                    
                      Ω
                      _
                    
                    k
                  
                  ⁢
                  
                    
                      G
                      _
                    
                    k
                  
                  ⁢
                  
                    
                      d
                      _
                    
                    k
                  
                
                +
                
                  
                    v
                    _
                  
                  k
                
              
            
          
        
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          41
          )
        
      
    
  




where Ωk=WHCkHHkHHkCkW,

vk is a vector of noise and interference for sector k, and
qk is a vector of decovered symbols for sector k.
The elements of qk may be expressed as:

E{|qk,n|2}=|αk,n|2gk,n2+σv,k2,  Eq (42)

where qk,n is the n-th element of qk,

αk,n and gk,n are the n-th diagonal elements of Ωk and Gk, respectively, and
σv,k2 is the variance of vk.
The square operation and expectation in equation (42) removes the data symbols in dk, which are assumed to be uncorrelated.
The noise and interference variance σv,k2 may be estimated by taking the difference of qk,pilot for the pilot channel for consecutive symbol periods, computing the squared magnitude of the difference, and filtering the squared magnitude to obtain the estimate of σv,k2. The estimated σv,k2 may then be subtracted from E{|qk,n|2} to obtain an estimate of |αk,n|2gk,n2, as follows:

Gk,n=E{|qk,n|2}−σv,k2=|αk,n|2gk,n2,  Eq (43)

where Gk,n is a scaled power gain for code channel n.

Scaled gains for the code channels may then be derived as follows:




  
    
      
        
          
            
              g
              
                k
                ,
                n
              
            
            
              g
              
                k
                ,
                pilot
              
            
          
          =
          
            
              
                G
                
                  k
                  ,
                  n
                
              
              
                G
                
                  k
                  ,
                  pilot
                
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          44
          )
        
      
    
  




where gk,pilot and Gk,pilot are the gain and the power gain, respectively, of the pilot channel for sector k.

In equation (44), the gains of the code channels are given relative to the gain gk,pilot of the pilot channel. This is a desired form since the channel response matrix Hk is also derived based on the pilot and includes the gain gk,pilot of the pilot channel, which would be canceled by the scaled gains from equation (44).
The processing shown in equations (43) to (44) may be performed for each sector. The gain matrix Gm for each user group may be formed based on the gains estimated for the users in that group.
FIG. 8 shows a block diagram of a code channel gain estimation unit 800. FIG. 8 shows the processing to estimate the gains of the N code channels for one sector k. Within gain estimation unit 800, a unit 810 performs channel matched filtering and multiplies the received samples in y with the complex conjugated channel gains in Hk for sector k. A multiplier 812 multiplies the output of unit 810 with the complex conjugated PN chips for sector k and provides despread samples. A unit 814 performs an N-point fast Hadamard transform (FHT) on N despread samples for each symbol period and provides N decovered symbols for N code channels, which are the N elements of qk. Unit 814 efficiently performs Walsh decovering for all N code channels.
A unit 820a computes the squared magnitude of the decovered symbol for each code channel. A filter 822a filters the output of unit 820a for each code channel. The output of filter 822a is an estimate of the expected value in equation (42).
The noise and interference variance is estimated based on the decovered symbols for the pilot channel. A unit 816 provides one symbol period of delay for each decovered symbol for the pilot channel. A summer 818 subtracts the delayed decovered symbol from the current decovered symbol and provides the difference. Since the pilot symbol is constant, taking the difference removes the pilot modulation while capturing the noise and interference, which are assumed to be random from symbol period to symbol period. A unit 820b computes the squared magnitude of the difference from summer 818 and further divides the result by two to account for the difference operation by summer 818. A filter 822b filters the output of unit 820b and provides the estimated noise and interference variance σv,k2.
A unit 824 subtracts the noise and interference variance from the output of filter 822a and provides the scaled power gain Gk,n for each code channel. A unit 826 determines a scaling factor 1/Gk,pilot based on the scaled power gain Gk,pilot for the pilot channel. A multiplier 828 multiplies the scaled power gain for each code channel with the scaling factor and provides Gk,n/Gk,pilot for each code channel. A unit 830 computes the square root of Gk,n/Gk,pilot for each code channel and provides the scaled gain gk,n/gk,pilot for that code channel.
The received power for each code channel may be determined based on the scaled gain gk,n for that code channel and the received power for sector k, as follows:




  
    
      
        
          
            P
            
              k
              ,
              n
            
          
          =
          
            
              
                (
                
                  
                    g
                    
                      k
                      ,
                      n
                    
                  
                  
                    g
                    
                      k
                      ,
                      pilot
                    
                  
                
                )
              
              2
            
            ⁢
            
              P
              
                k
                ,
                pilot
              
            
          
        
        ,
      
    
    
      
        Eq
        ⁢
        
          
        
        ⁢
        
          (
          45
          )
        
      
    
  




where Pk,pilot is the received pilot power for sector k, and

Pk,n is the received power for code channel n of sector k.
FIG. 9 shows an embodiment of a process 900 for performing iterative detection and cancellation. Multiple groups of code channels are formed for a plurality of code channels for at least one sector (block 910). Each group may include all code channels for one sector (for sector-based processing), a subset of all code channels for one sector (for local processing), or multiple code channels for multiple sectors (for global processing). Processing is performed for the multiple groups of code channels in multiple iterations (block 920). For each iteration, data detection and signal cancellation are performed for the multiple groups of code channels in multiple stages, e.g., in a sequential order starting with the first group having the strongest received powers and concluding with the last group having the weakest received powers (block 930).
FIG. 10 shows an embodiment of a detection and cancellation process 1000 performed by one stage in one iteration. Process 1000 may be used for block 930 in FIG. 9. Data detection is performed to obtain a detected signal for a group of code channels being processed by the stage (block 1012). A signal for the group of code channels is reconstructed based on the detected signal (block 1014). The reconstructed signal is canceled to obtain an output signal for the stage (block 1016).
For the data detection in block 1012, a system matrix may be determined for the group of code channels. The variance of noise and interference for the group of code channels may also be estimated. A data detection filter for the group of code channels may then be derived based on the system matrix and possibly the variance of the noise and interference, e.g., in accordance with the linear MMSE or least squares technique as shown in equation (19) or (20), respectively. Data detection is then performed with the data detection filter to obtain the detected signal for the group of code channels, e.g., as shown in equation (19) or (20).
FIG. 11 shows an embodiment of an equalization, detection and cancellation process 1100 performed by one stage in one iteration. Process 1100 may also be used for block 930 in FIG. 9. Equalization is performed to obtain an equalized signal for a group of code channels being processed by the stage (block 1110). Data detection is then performed on the equalized signal to obtain a detected signal for the group of code channels (block 1112). A signal for the group of code channels is reconstructed based on the detected signal (block 1114). The reconstructed signal is canceled to obtain an output signal for the stage (block 1116).
For the equalization in block 1110, an equalizer matrix may be derived for the group of code channels, e.g., in accordance with the linear MMSE or least squares technique as shown in equation (21) or (22), respectively. Equalization is then performed with the equalizer matrix to obtain the equalized signal.
For the data detection in block 1112, a gain matrix may be estimated for the group of code channels. The covariance of the noise and interference for the group of code channels may also be estimated. A data detection filter for the group of code channels may then be derived based on the gain matrix and possibly the covariance of the noise and interference, e.g., in accordance with the linear MMSE or least squares technique as shown in equation (27) or (28), respectively. Data detection is then performed with the data detection filter to obtain the detected signal for the group of code channels, e.g., as shown in equation (27) or (28).
For process 1000 in FIG. 10 and process 1100 in FIG. 11, for each subsequent iteration after the first iteration, a reconstructed signal obtained in a prior iteration for the group of code channels may be added back to an input signal for the stage to obtain a reconstituted signal. Data detection and signal cancellation may then be performed for the group of code channels based on the reconstituted signal.
Back in FIG. 9, for causal cancellation, a signal for the prior symbol period may be removed from the received signal. The processing in blocks 920 and 930 may then be performed after removing the signal for the prior symbol period. For causal and non-causal cancellation, a signal for the prior symbol period and a signal for the next symbol period may be removed from the received signal. The processing in blocks 920 and 930 may then be performed after removing the signals for the prior and next symbol periods.
The techniques described herein may be implemented by various means. For example, these techniques may be implemented in hardware, firmware, software, or a combination thereof. For a hardware implementation, the processing units used to perform the processing described herein 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, electronic devices, other electronic units designed to perform the functions described herein, or a combination thereof.
For a firmware and/or software implementation, the techniques may be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The firmware and/or software codes may be stored in a memory (e.g., memory 292 in FIG. 2) and executed by a processor (e.g., processor 290). The memory may be implemented within the processor or external to the processor.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the disclosure. 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 disclosure. Thus, the disclosure 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. An apparatus comprising:
at least one processor having multiple processing stages, configured to perform processing for multiple groups of code channels in multiple iterations and, for each iteration, to perform data detection and signal cancellation for each of the multiple groups of code channels in a respective one of said multiple processing stages;
wherein an output of each iteration for the multiple groups of code channels is input into a next following stage to be used to perform a next iteration.

2. The apparatus of claim 1, wherein for each stage of each iteration, the at least one processor is configured to perform data detection to obtain a detected signal for a group of code channels, to reconstruct a signal for the group of code channels based on the detected signal, and to cancel the reconstructed signal to obtain an output signal for the stage.

3. The apparatus of claim 2, wherein the at least one processor is configured to determine a system matrix for the group of code channels, to derive a data detection filter for the group of code channels based on the system matrix, and to perform data detection with the data detection filter to obtain the detected signal for the group of code channels.

4. The apparatus of claim 3, wherein the at least one processor is configured to estimate variance of noise and interference for the group of code channels and to derive the data detection filter further based on the variance of the noise and interference.

5. The apparatus of claim 3, wherein the at least one processor is configured to derive the data detection filter in accordance with linear minimum mean square error (MMSE) or least squares technique.

6. The apparatus of claim 1, wherein for each stage of each iteration, the at least one processor is configured to perform equalization to obtain an equalized signal for a group of code channels, to perform data detection on the equalized signal to obtain a detected signal for the group of code channels, to reconstruct a signal for the group of code channels based on the detected signal, and to cancel the reconstructed signal to obtain an output signal for the stage.

7. The apparatus of claim 6, wherein the at least one processor is configured to derive an equalizer matrix for the group of code channels and to perform equalization with the equalizer matrix to obtain the equalized signal.

8. The apparatus of claim 7, the at least one processor is configured to derive the equalizer matrix in accordance with linear minimum mean square error (MMSE) or least squares technique.

9. The apparatus of claim 6, wherein the at least one processor is configured to estimate a gain matrix for the group of code channels, to derive a data detection filter for the group of code channels based on the gain matrix, and to perform data detection on the equalized signal with the data detection filter to obtain the detected signal for the group of code channels.

10. The apparatus of claim 9, wherein the at least one processor is configured to estimate covariance of noise and interference for the group of code channels and to derive the data detection filter further based on the covariance of the noise and interference.

11. The apparatus of claim 1, wherein for each stage of each iteration after a first iteration, the at least one processor is configured to add a first signal obtained in a prior iteration for a group of code channels to an input signal to obtain a second signal, and to perform data detection and signal cancellation for the group of code channels based on the second signal.

12. The apparatus of claim 1, wherein the at least one processor is configured to remove a signal for a prior symbol period and to perform processing for the multiple groups of code channels in multiple iterations after removing the signal for the prior symbol period.

13. The apparatus of claim 1, wherein the at least one processor is configured to remove a signal for a prior symbol period and a signal for a next symbol period and to perform processing for the multiple groups of code channels in multiple iterations after removing the signals for the prior and next symbol periods.

14. The apparatus of claim 1, wherein for each of the multiple groups of code channels, the at least one processor is configured to obtain decovered symbols for the code channels in the group, to determine an averaged squared magnitude of the decovered symbols for each code channel in the group, and to determine a gain for each code channel in the group based on the averaged squared magnitude of the decovered symbols for the code channel.

15. The apparatus of claim 14, wherein for each of the multiple groups of code channels, the at least one processor is configured to estimate variance of noise and interference based on decovered symbols for a pilot channel, and to determine the gain for each code channel further based on the variance of the noise and interference.

16. The apparatus of claim 1, wherein the at least one processor is configured to determine received powers of a plurality of code channels for at least one sector, and to form the multiple groups of code channels based on the received powers.

17. The apparatus of claim 16, wherein the at least one processor is configured to perform data detection and signal cancellation for the multiple groups of code channels in a sequential order for each iteration, starting with a first group having strongest received powers and concluding with a last group having weakest received powers.

18. The apparatus of claim 1, wherein each group includes all code channels for one sector.

19. The apparatus of claim 1, wherein each group includes a subset of all code channels for one sector.

20. The apparatus of claim 1, wherein each group includes multiple code channels for multiple sectors.

21. A method comprising:
performing processing for multiple groups of code channels in multiple iterations; and performing data detection and signal cancellation for each of the multiple groups of code channels in a respective one of multiple stages for each iteration;
wherein an output of each iteration for the multiple groups of code channels is input into a next following stage to be used to perform a next iteration.

22. The method of claim 21, wherein the performing data detection and signal cancellation comprises, for each stage of each iteration,
performing data detection to obtain a detected signal for a group of code channels,
reconstructing a signal for the group of code channels based on the detected signal, and
canceling the reconstructed signal to obtain an output signal for the stage.

23. The method of claim 21, wherein the performing data detection and signal cancellation comprises, for each stage of each iteration,
performing equalization to obtain an equalized signal for a group of code channels,
performing data detection on the equalized signal to obtain a detected signal for the group of code channels,
reconstructing a signal for the group of code channels based on the detected signal, and
canceling the reconstructed signal to obtain an output signal for the stage.

24. The method of claim 21, wherein the performing data detection and signal cancellation comprises, for each stage of each iteration after a first iteration,
adding a first signal obtained in a prior iteration for a group of code channels to an input signal to obtain a second signal, and
performing data detection and signal cancellation for the group of code channels based on the second signal.

25. The method of claim 21, further comprising:
determining received powers of a plurality of code channels for at least one sector; and
forming the multiple groups of code channels based on the received powers, and
wherein the performing data detection and signal cancellation comprises performing data detection and signal cancellation for the multiple groups of code channels in a sequential order for each iteration, starting with a first group having strongest received powers and concluding with a last group having weakest received powers.

26. An apparatus comprising:
means for performing processing for multiple groups of code channels in multiple iterations; and
means for performing data detection and signal cancellation for each of the multiple groups of code channels in a respective one of multiple stages for each iteration;
wherein an output of each iteration for the multiple groups of code channels is input into a next following stage to be used to perform a next iteration.

27. The apparatus of claim 26, wherein the means for performing data detection and signal cancellation comprises, for each stage of each iteration,
means for performing data detection to obtain a detected signal for a group of code channels,
means for reconstructing a signal for the group of code channels based on the detected signal, and
means for canceling the reconstructed signal to obtain an output signal for the stage.

28. The apparatus of claim 26, wherein the means for performing data detection and signal cancellation comprises, for each stage of each iteration,
means for performing equalization to obtain an equalized signal for a group of code channels,
means for performing data detection on the equalized signal to obtain a detected signal for the group of code channels,
means for reconstructing a signal for the group of code channels based on the detected signal, and
means for canceling the reconstructed signal to obtain an output signal for the stage.

29. The apparatus of claim 26, wherein the means for performing data detection and signal cancellation comprises, for each stage of each iteration after a first iteration,
means for adding a first signal obtained in a prior iteration for a group of code channels to an input signal to obtain a second signal, and
means for performing data detection and signal cancellation for the group of code channels based on the second signal.

30. The apparatus of claim 26, further comprising:
means for determining received powers of a plurality of code channels for at least one sector; and
means for forming the multiple groups of code channels based on the received powers, and
wherein the means for performing data detection and signal cancellation comprises means for performing data detection and signal cancellation for the multiple groups of code channels in a sequential order for each iteration, starting with a first group having strongest received powers and concluding with a last group having weakest received powers.

31. A processor-readable memory including processor-executable instructions stored therein for performing a method comprising:
performing processing for multiple groups of code channels in multiple iterations; and performing data detection and signal cancellation for each of the multiple groups of code channels in a respective one of multiple stages for each iteration;
wherein an output of each iteration for the multiple groups of code channels is input into a next following stage to be used to perform a next iteration.