# Patent

## Title
Method and apparatus for two-user joint demodulation in a system having transmit diversity

## Abstract

Method and apparatus for two-user join demodulation in a system having transmit diversity. A receiver includes at least a detector unit and a channel estimation unit, and exploits transmit diversity employed by at least one of two users or transmitters on a channel. Symbol detection and channel estimation are performed in block or batch fashion, or recursively on a symbol-by-symbol basis. Joint detection makes use of transmit diversity by at least one of the transmitting users. Channel estimates can be updated using channel tracking. If the signals from the first user and the second user are asynchronous, detection can be accomplished in part through reference to a pulse-shape component. The receiver may be incorporated into a mobile communication terminal.

## Background
BACKGROUND OF THE INVENTION
The performance of a wireless communication system can be improved through the addition and exploitation of signal “diversity” at one or more points in the system. A signal is said to have diversity when it can be divided into multiple (usually two) different signals, each carrying the same information, but differing in some respect. The term “diversity” typically refers to spatial diversity that results when a signal is sampled at, or transmitted from two different points in space, as by two antennas. In such an arrangement, the signal is said to have a different “diversity branch” for each antenna.
Transmit diversity is an important method for improving reception in the situation where limited diversity is available at the receiver, for example, in the case of single-antenna handsets. If the handset has a priori knowledge of the transmit diversity scheme, the receiver in the handset can use the knowledge to demodulate a signal more effectively in the presence of noise. A space-time block code is used under the assumption of flat channels to achieve the same diversity benefit as would be obtained by maximal ratio combining of signals from two antennas at the receiver. This approach can be extended for dispersive channels.
FIG. 1 illustrates the approach. Base station 101 transmits from two antennas simultaneously. In FIG. 1, two symbols, a1 and a2, are transmitted from two antennas so that the received signal combines the symbols in a favorable manner. During the first transmission, a1 is transmitted from antenna 102 and a2 is transmitted from antenna 103. During the next transmission, −a2* is transmitted from antenna 102 and a1* is transmitted from antenna 103. Assuming that fading is flat for each signal, Nyquist signaling is used so that there is no inter-symbol interference (ISI) from the filtered and sampled signal, and the channel is constant over the two symbol transmissions. The two received, sampled signals, r1 and r2 can be described by the following vector/matrix form. In the following discussion, lowercase variables denote vector quantities and uppercase variables denote matrix quantities:




  
    [
    
      
        
          
            r
            1
          
        
      
      
        
          
            r
            2
          
        
      
    
    ]
  
  =
  
    
      
        [
        
          
            
              
                
                  a
                  1
                
              
            
            
              
                
                  -
                  
                    a
                    2
                    *
                  
                
              
            
          
          ⁢
          
            
              
                
                  a
                  2
                
              
            
            
              
                
                  a
                  1
                  *
                
              
            
          
        
        ]
      
      ⁡
      
        [
        
          
            
              
                c
                1
              
            
          
          
            
              
                c
                2
              
            
          
        
        ]
      
    
    +
    
      [
      
        
          
            
              ω
              1
            
          
        
        
          
            
              ω
              2
            
          
        
      
      ]
    
  




where c1 and c2 represent channel estimates between the two antennas and wireless terminal 104 of FIG. 1, and ω1 and ω2 represent extraneous noise. For some processing, it is useful to reformulate the channel estimates and symbol representations using r2* in place of r2, yielding:





  
    [
    
      
        
          
            r
            1
          
        
      
      
        
          
            r
            2
            *
          
        
      
    
    ]
  
  =
  
    
      
        [
        
          
            
              
                
                  c
                  1
                
              
            
            
              
                
                  c
                  2
                  *
                
              
            
          
          ⁢
          
            
              
                
                  c
                  2
                
              
            
            
              
                
                  -
                  
                    c
                    1
                    *
                  
                
              
            
          
        
        ]
      
      ⁡
      
        [
        
          
            
              
                a
                1
              
            
          
          
            
              
                a
                2
              
            
          
        
        ]
      
    
    +
    
      
        [
        
          
            
              
                ω
                1
              
            
          
          
            
              
                ω
                2
                *
              
            
          
        
        ]
      
      .
    
  




For the kth pair of transmitted symbols {a1(k); a2(k)}, and the kth received sample pair {r1(2k); r2(2k+1)}, the above equations are written as:

yk=Akck+vk,

rk=Ckak+ωk

respectively.

An example receiver used with the transmit diversity system shown above is illustrated in FIG. 2. The receiver includes a vector data generator, 201, a channel estimation unit, 202, and a detector unit, 203. The signal stream is represented by {a(k)}, a portion of which, sp, is used for channel estimation. The symbols ap represent symbols from a known pilot field or a known pilot carrier for use in obtaining the estimated channel coefficients. The received samples {r(k)} are used to first form the vector yp, which corresponds to the known pilot symbols. These received samples are used to estimate the channel coefficients ĉ. For example, the channel coefficients can be obtained via least-squares estimation since the symbol data are known. The estimate, ĉLS is given by:

ĉLS=(ApHAp)−1ApHyp.

This equation reduces to:

ĉLS=ApHyp,

in the case of flat fading and under the assumption that:

|a1|2=|a2|2=1.

Note, there may be more than one known symbol pair, and ĉLS can be easily modified to account for this, under the assumption that the channel is fixed for an estimation period. Once ĉLS is estimated, then the received sample pair rk corresponding to the unknown transmitted symbol pair ak is used to detect these symbols. In the case that no known pilot symbols exist, the above estimation approaches can be used iteratively to generate estimates of the channel and detected symbol values.

Increasingly, receiving systems must demodulate and extract information from two users (transmitters) occupying the same bandwidth or channel. In some cases, both signals contain desired information. More typically however, one user is transmitting a desired signal, and a second user represents co-channel interference. Co-channel interference has become a significant problem as the population density of people using wireless, cellular communication devices increases. Because of the traffic demand, it is necessary to re-use channels more frequently on a geographic basis. Joint demodulation schemes have been used successfully to demodulate signals from two users on the same channel. Such schemes do not make use of diversity and only work when the signal-to-noise ratio is fairly high. There are also co-channel interference canceling methods that make use of diversity at the receiver.
BRIEF SUMMARY OF THE INVENTION
The present invention provides methods and apparatus for joint demodulation of signals from two users or transmitters on the same channel. The method of the invention exploits transmit diversity employed by at least one of the users. In one embodiment of the invention, symbol detection and channel estimation are performed in block or batch fashion. The received signal corresponds to a plurality of transmitted symbol sequences or blocks, originating from both a first user and a second user. A stored symbol sequence corresponding to the first user, the first-user symbol sequence, as well as a stored symbol sequence corresponding to the second user, the second-user symbol sequence, are initialized, possibly through a semi-blind initialization procedure. Stored channel estimates are also initialized. In a block-by-block, iterative fashion, received first-user symbol sequences and second-user symbol sequences are jointly and iteratively detected using the channel estimates, and based on the use of transmit diversity by at least one of the users. Detected sequences are used to update stored sequences as all sequences are demodulated. Channel estimates may also be updated on a block-by-block basis, unless they are available from another source, such as an alternate pilot carrier in a code division multiple access (CDMA) system.
In another embodiment, symbols are updated, and channel estimates are updated, if necessary, recursively on a symbol-by-symbol basis. The received signal corresponds to a stream of a plurality of transmitted symbols from each user. A specific first-user symbol and a specific second-user symbol are detected at each of a plurality of symbol times, and stored symbols are updated. The detection process can use previously detected symbols as well as channel estimates. Again, joint detection makes use of transmit diversity by at least one of the transmitting users. Channel estimates can be updated at each current symbol time using channel tracking. If the signals from the first user and the second user are asynchronous, detection can be accomplished through reference to a pulse-shape component.
A receiving system which implements the invention includes a channel estimation unit operable to produce channel estimates for the first user and the second user, and a detector unit operatively connected to the channel estimation unit. The detector unit is operable to jointly detect a first-user symbol sequence and a second-user symbol sequence or first-user and second-user symbols at specific times. The receiving system may also include initialization logic. The detector unit may include a symbol update unit. In one embodiment, the symbol update unit and channel estimation unit are made up of interconnected reformulation operators and conjugation operators which can perform complex conjugate operations, as well as interconnected multipliers and adders. These components may be implemented in hardware, microcode, software, or a combination, and form the means to carry out the invention. In one embodiment they form at least a portion of the baseband logic of a processor controlled, wireless terminal.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a diagram illustrating the concept of transmit diversity through the use of two antennas.
FIG. 2 is a block diagram of a prior-art receiver that makes use of transmit diversity.
FIG. 3 is a flowchart illustrating the method of the invention according to one embodiment.
FIG. 4 is a flowchart illustrating the method of the invention according to one embodiment.
FIG. 5 is a flowchart illustrating a method of blind initialization that can be used in conjunction with the invention.
FIG. 6 is a flowchart illustrating the method of the invention according to another embodiment.
FIG. 7 is a flowchart that illustrates the operation of a detector unit according to an embodiment of the invention.
FIG. 8 is a signal processing flow diagram that illustrates a symbol update unit that may be used within the detector unit of FIG. 7.
FIG. 9 illustrates the operation of a channel estimation unit according to an embodiment of the invention in flowchart form.
FIG. 10 is a signal processing flow diagram that illustrates the operation of a channel estimation unit according to the invention.
FIG. 11 is a flowchart that illustrates the operation of a detector unit that handles inter-symbol interference (ISI) according to an embodiment of the invention.
FIG. 12 is a signal processing flow diagram that illustrates the operation of a symbol update unit that handles ISI according to an embodiment of the invention.
FIG. 13 is a flowchart that illustrates channel tracking according to an embodiment of the invention.
FIG. 14 is a signal processing flow diagram that further illustrates channel tracking according to an embodiment of the invention.
FIG. 15 is a block diagram of a mobile terminal that makes use of the present invention.



DETAILED DESCRIPTION OF THE INVENTION
The present invention is typically embodied in a receiving system for a communication network where information is transmitted and received as a stream of symbols modulated onto a carrier wave. One example would be a network employing time domain multiplexed access (TDMA). The embodiments disclosed are not restricted to any particular standard, however. The well-known Global System for Mobile (GSM) communication system is one example of a TDMA system which might employ the invention. The invention can also be used in a wireless data network employing code division multiplexed access (CDMA).
It should also be understood that not every feature of the receiving system described is necessary to implement the invention as claimed in any particular one of the appended claims. Various elements of receivers are described to fully enable the invention. It should also be understood that throughout this disclosure, where a process or method is shown or described, the steps of the method may be performed in any order or simultaneously, unless it is clear from the context that one step depends on another being performed first.
Some of the diagrams, which are used to illustrate the inventive concepts, are shown in software flowchart form, while others are shown as block diagrams and signal processing flow diagrams with various operators. The particular form of illustration is selected in each case to aid in understanding the invention. It is to be understood that the concepts discussed can be implemented in hardware or software. For example, a flowchart which is described as a detector unit or a detection method can be thought of as illustrating an algorithm performed by software, or the operation of a specific piece of hardware, such as an application specific integrated circuit (ASIC) that performs the illustrated task.
Finally the usage of terms such as “updating” and “initializing” typically refer to updating symbols or vectors stored or being stored in order to be operated on, decoded, or demodulated. The term user, as in first user and second user is colloquially being used to refer to a transmitter or sender of information, notwithstanding that most aspects of the invention are being implemented in a receiver or a receiver portion of a transceiver in a communication system.
The invention consists of the use of joint detection of at least two users to improve performance. In this approach, the objective is to minimize the quantity:

|r−{circumflex over (r)}(a,b)|2

under the same assumptions as the single user case, to obtain values for new estimates for users a and b using:





  a
  ,
  
    b
    =
    
      
        arg
        ⁢
        
          
        
        ⁢
        
          min
          
            a
            ,
            b
          
        
      
      |
      
        r
        -
        
          
            r
            ^
          
          ⁡
          
            (
            
              a
              ,
              b
            
            )
          
        
      
      ⁢
      
        |
        2
      
      .
    
  




An equivalent log-likelihood to be used for detection is:

γ=Re{aHγa+bHγb−aHHb}

where H=CHG and the above equation needs to be maximized.

Channel estimates can be obtained via least-squares estimation if the symbol data are known. In the case of flat fading, and assuming |a|2=1, they reduce to:

ĉLS=AHy

For the two-user scenario, with flat fading, the estimate for two user's channels, c and g, becomes:





  
    
      [
      
        
          
            
              c
              ^
            
          
        
        
          
            
              g
              ^
            
          
        
      
      ]
    
    LS
  
  =
  
    
      [
      
        
          
            
              A
              H
            
          
        
        
          
            
              B
              H
            
          
        
      
      ]
    
    ⁢
    y
  




and this equation shows each user interfering with the other user's channel estimate. An alternative estimation technique is the least-squares estimate with the interfering term removed from the received signal. Now, the estimates for each user's channel can be found using:

ĉLS=AH(y−Bg),

ĝLS=BH(y−Ac).

Since the interfering users channel estimate must be known in order to find the new estimate, values from prior iterations can be used to calculate new estimates.

One approach is to calculate the single-user estimates using the received signal with the interfering users term removed. Equations for these estimates are written as:




  
    
      a
      ^
    
    =
    
      
        1
        α
      
      ⁢
      
        
          C
          H
        
        ⁡
        
          (
          
            r
            -
            Gb
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      b
      ^
    
    =
    
      
        1
        β
      
      ⁢
      
        
          
            G
            H
          
          ⁡
          
            (
            
              r
              -
              Ca
            
            )
          
        
        .
      
    
  




An assumption is made that the channel estimates for both users are known and prior estimates of the symbol values are available from some prior iteration. The above procedure provides continuous value estimates for the symbols, which may be quantized (hard-detected) to the closest transmit symbol. The channel estimation and detection steps are iteratively computed until some convergence criterion is met or some maximum number of iterations is reached.

Another iterative approach can be used to compute the maximum a poste-riori estimate of the possible transmit symbols. This approach is described in U.S. Pat. No. 5,887,035, which is incorporated herein by reference.
One embodiment of the invention is based on a two-user situation in which the base stations are synchronized. This can easily be accomplished if the system is designed to allow for joint demodulation according to the invention. It is also the case if two users are transmitting from the same base station. In this case, the received signal can be written as:

rk=Ckak+Gkbk+ωk

where bk is the vector of symbols transmitted from a second base station and Gk is the corresponding channel response. Different signals are transmitted from the multiple transmit antennas and channel estimates will contain the effects of any transformation from the originating signal stream to the multiple antennas.

For this example, a form of joint detection is illustrated which assumes that both signals use transmit diversity, the medium response is flat faded, channel coefficients are not known, and there are no training symbols. To address this problem, an iterative solution to two suboptimal problems, detection and channel estimation is used. The cases where some of the above quantities are known require some modifications to the approach that can easily be accomplished by one of ordinary skill in the art.
Iteration is used to address the transmit diversity problem when multiple users (desired and interfering) are present. Iteration can be performed in a block (i.e. batch) mode, or in a recursive mode of operation. FIG. 3 is a flowchart illustrating one embodiment of the block mode. In this embodiment, a detector unit performs block detection at 301 for the symbols:

â={â1, . . . , âK}, and

{circumflex over (b)}={{circumflex over (b)}1, . . . , {circumflex over (b)}K},

given channel estimates:

ĉ={ĉ1 . . . ĉK}, and

ĝ={ĝ1, . . . ĝK}.

The terms ã and {tilde over (b)} are either initial values for the first iteration or prior detected values during subsequent iterations. In either case these can be referred to as “previous” values, simply meaning they are already stored and available. These values can be initially set at the initialization step, 300. Initially setting these values to zero results in single-user detection for the first pass of the detector unit. In this example, the channel estimates are represented as complete sequences; however, a small set of channel estimates may be used when the channel is not time-varying or only slowly time-varying.

At step 302 of FIG. 3, the values of ĉ and ĝ are copied to {tilde over (c)} and {tilde over (g)}, respectively for further processing in a channel estimation unit at step 303, where new channel estimates are computed. At step 305 a determination as to whether to perform another iteration is made. If another iteration is to be performed, the receiver copies the values â and {circumflex over (b)} to ã and {tilde over (b)}, respectively at step 304. The detector unit calculates the values of â and {circumflex over (b)} with ĉ and ĝ fixed. The values of â, {circumflex over (b)}, ã and {tilde over (b)} may be either hard or soft values.
FIG. 4 shows a similar embodiment, but initial values are now â, {circumflex over (b)}, {tilde over (c)} and {tilde over (g)}. Setting {tilde over (c)} and {tilde over (g)} to have zero initial values at initialization step 400 forces channel estimation at step 403 to perform single-user channel estimation for the first iteration. The rest of the processes are the same, with steps 401, 402, 404, and 405 being essentially the same processes that take place in steps 301, 302, 304, and 305 of FIG. 3. This embodiment is similar to the one shown in FIG. 3, but the iterative algorithm has different starting and termination points.
In the algorithm of FIG. 3, if the channel estimation unit 303 is skipped, then iteration only executes detection at step 301. This approach could be used when channel estimates are known or can be obtained by other means (such as the use of alternate pilot channels in a CDMA system). Similarly, in the second embodiment, if detection at step 401 is skipped during iteration, then only channel estimation at step 403 is executed. This might occur, for example, while detecting known pilot symbols that are being used for initial channel estimation.
The initialized values are assumed to be known in the above examples. However, it may be the case that the second user signal is an interferer and the values corresponding to the interferer are unknown. FIG. 5 illustrates a semi-blind initialization procedure that can be used in this case. Assuming that â is known and can be initialized at step 500, the single user channel estimate ĉ is computed at step 501. Next, some rough estimate of ĝ is generated at step 502. Then â is copied to ã at step 503. {tilde over (b)} is set to zero and the initial conditions are output at 504.
The rough estimate of ĝ discussed above can be obtained by training the desired user's channel estimate over known symbols corresponding to the desired user, while using a sequence estimator to compute the channel estimate for the interfering user. Such an approach is described in U.S. Pat. No. 6,304,618 which is incorporated herein by reference.
Rather than use a block estimation approach, an embodiment that uses recursive estimation can also be used. FIG. 6 shows a flowchart for the update at symbol time k. Steps 601 through 605 are similar to the corresponding steps 301–305 in FIG. 3, however, detection and channel estimation units only operate on the kth symbol or channel estimation term. At step 606 an update is made of the channel estimates from the kth to the (k+1)st symbol time. Known, conventional channel tracking techniques may be used to predict the channel estimates from one time to a future time. A recursive embodiment corresponding to the block approach shown in FIG. 4 can also be developed but is not shown here.
In the case of recursive detection where the signals are synchronous, the detection at each symbol time can be performed separately. There are multiple ways to do this, and one of ordinary skill in the art can select the appropriate way for a proposed system. For the kth symbol time, the problem is to find the values of âk and {tilde over (b)}k using the equation:




  
    
      a
      ^
    
    k
  
  ,
  
    
      
        b
        k
      
      ^
    
    =
    
      
        arg
        ⁢
        
          
        
        ⁢
        
          min
          
            ak
            ,
            bk
          
        
      
      |
      
        
          r
          k
        
        -
        
          
            
              r
              ^
            
            k
          
          ⁡
          
            (
            
              
                c
                k
              
              ,
              
                g
                k
              
            
            )
          
        
      
      ⁢
      
        |
        2
      
      .
    
  




An equivalent log likelihood to be used for detection is:

γ=Re{aHγa+bHγb−aHHkb},

where Hk=CkHGk. The object here is to detect the symbol from the first and second user at the current symbol time that maximizes the above equation. Although this can be done directly, each symbol can no longer be detected independently, and so the complexity of the detection procedure is significant.

A lower complexity approach is to calculate the single-user estimates using the received signal with the interfering user terms removed as previously discussed. In the recursive case the equations become:




  
    
      
        a
        ^
      
      k
    
    =
    
      
        1
        α
      
      ⁢
      
        
          C
          k
          H
        
        ⁡
        
          (
          
            
              r
              k
            
            -
            
              
                G
                k
              
              ⁢
              
                
                  b
                  ~
                
                k
              
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        b
        ^
      
      k
    
    =
    
      
        1
        β
      
      ⁢
      
        
          G
          k
          H
        
        ⁡
        
          (
          
            
              r
              k
            
            -
            
              
                C
                k
              
              ⁢
              
                
                  a
                  ~
                
                k
              
            
          
          )
        
      
    
  
  ,




where ãk and {tilde over (b)}k are previous estimates of ak and bk, respectively. A flowchart showing the method executed by a detector unit using this approach is shown in FIG. 7. After initialization at 700, an update unit updates of âk and {circumflex over (b)}k using channel estimates ĉk and ĝk, and prior detected values of ãk and {tilde over (b)}k at step 701. At step 702 a determination is made as to whether another iteration is needed. If not, final values of âk and {circumflex over (b)}k are output at 704. If so, current symbol values are copied to prior detected values at step 703.

FIG. 8 shows a signal processing flow/block diagram for a symbol update unit around which a detector unit according to one embodiment of the invention can be build. The values ĉk and ĝk are reformed into Ck and Gk, respectively, by reformulation operators 801 and 802. Conjugation operators 803 and 804 provide complex conjugate values. Logic such as that shown in FIG. 8 provides continuous, final estimates for the symbols. For symbol aligned signals, batch detection can be accomplished by calling a recursive detection unit once for each symbol time given the appropriate input values for each call.
Alternate detection approaches can be used. For example, rather than use previously detected values for ãk and {tilde over (b)}k, these terms could be enumerated. Then, the chosen values of âk and {circumflex over (b)}k will be determined by choosing hypotheses that correspond to the best metric. Another alternate detection approach is to compute the maximum a posteriori estimate of the possible transmit symbols. This is more complicated in that probabilities for each transmit symbol are estimated during each iteration. U.S. Pat. No. 5,887,035, previously discussed and incorporated by reference, provides useful background that can be used to develop alternate detection approaches.
The recursive channel estimation problem at time k is to find estimates of ĉk and ĝk from the equation:




  
    
      
        c
        ^
      
      k
    
    ⁢
    
      
        g
        ^
      
      k
    
  
  =
  
    
      arg
      ⁢
      
        
      
      ⁢
      
        min
        
          
            c
            k
          
          ⁢
          
            g
            k
          
        
      
    
    |
    
      
        y
        k
      
      -
      
        
          
            y
            ^
          
          k
        
        ⁡
        
          (
          
            
              a
              k
            
            ,
            
              b
              k
            
          
          )
        
      
    
    ⁢
    
      |
      2
    
    .
  




Since joint estimation was previously shown to result in the single-user estimator with interference terms, an approach similar to that employed in detection must be used. The channel update is performed according to the equations:

ĉk=AkH(yk−Bk{tilde over (g)}k),

ĝk=BkH(yk−Ak{tilde over (c)}k)

where {tilde over (c)}k and {tilde over (g)}k are previous estimates of ck and gk, respectively. Since the interfering user's channel estimate must be known in order to find the new estimate, values from prior iterations can be used to calculate these new estimates.

In the above approach, ĉk and ĝk depend only upon the single received sample vector yk, which may be subject to noise and other signal interference. Since there typically is correlation in time from one sample to the next, channel tracking can be used. Consider the least-means-square (LMS) update for the channel coefficient c1(k). In the case that the transmit diversity approach has two received samples at two different symbol times, kε{2j,2j+1}, the LMS updates are found using the following equations:

ĉ(2j+1)=ĉ(2j)+μa1(k)*e1(k),

ĉ(2j+2)=ĉ(2j+1)−μa2(k)e2(k).

A similar form exists for the other channel coefficients. Additional background on LMS tracking can be found in Haykin, Adaptive Filter Theory, third edition, Prentice Hall, 1996, which is incorporated herein by reference.

The term ĉk is the update given a new estimate of ĝk. Using exponential smoothing, this update is written as:

ĉk=(1−μ){tilde over (c)}k+μAkH(yk−Bk{tilde over (g)}k).

Letting:

ek=yk−Akĉk−Bkĝk,

gives the update:

ĉk=ĉk+μAkHek,

giving LMS-like updates for the vector channel estimates. When no previous estimate for ĉk is available, then estimates from previous time k−1 can be used to give:

ĉk=ĉk−1+μAkHek.

The update for ĝk is similar.

A flowchart for the above approach is shown if FIG. 9. All vectors are initialized at step 900. At step 901, the update as described immediately above is performed. Steps 902 and 903 provide iteration. Once the iteration is complete, channel tracking is performed at step 904 using the step size term as input. Final values are output at step 905.
FIG. 10 is a signal flow diagram for a channel estimation unit that is used in performing the method illustrated in FIG. 9. The values âk and {circumflex over (b)}k are reformed into Ak and Bk, respectively, by reformulation operators 1001 and 1002. Conjugation operators 1003 and 1004 provide complex conjugate values.
One approach for batch channel estimation is to assume that the channel is fixed so that joint least-squares (LS) estimation can be performed using multiple received samples. This approach can work well since there is no longer an undetermined system of equations that need to be solved. For example, if the received signals yk through yk+N are used to find the joint LS:




  
    
      
        [
        
          
            
              
                
                  c
                  ^
                
                k
              
            
          
          
            
              
                
                  g
                  ^
                
                k
              
            
          
        
        ]
      
      LS
    
    =
    
      
        
          (
          
            [
            
              
                
                  
                    A
                    k
                  
                
                
                  
                    B
                    k
                  
                
              
              
                
                  
                    A
                    
                      k
                      +
                      1
                    
                  
                
                
                  
                    B
                    
                      k
                      +
                      1
                    
                  
                
              
              
                
                  ⋮
                
                
                  ⋮
                
              
              
                
                  
                    A
                    
                      k
                      +
                      N
                    
                  
                
                
                  
                    B
                    
                      k
                      +
                      N
                    
                  
                
              
            
            ]
          
          )
        
        +
      
      ⁢
      
        y
        k
      
    
  
  ,




where ( )+ is the pseudo-inverse operation. Since the symbol terms may not be able to be pre-computed, a recursive approach may be desired. In this case, the recursive channel tracking approach described above can be used for each symbol.

An alternative approach is needed when the two transmitted user signals are flat faded, but have misaligned symbol timing from each other at the receiver, resulting in intersymbol interference (ISI) from one signal to the other. This ISI appears, even in flat fading, since the pulse-shape responses are aligned so that they no longer have the same zero crossings. With dispersion, sampling the received signal at a higher, but uniform, sampling rate helps make the resulting receiver more tolerant to ISI. A higher sampling rate is required in the case where the medium response is not flat faded and significant dispersion results.
For the scenario above, non-uniform sampling at a rate greater than the symbol rate, but taking into account transmit diversity is used. The approaches discussed previously are modified so that they can be used in this case. Sequence estimation approaches can be modified for use with multi-user transmit diversity, but the same method as previously described for recursive channel estimation is not used. Instead, batch estimation and a sub-optimal, recursive form for detecting symbols and computing channel estimates is employed.
Consider the two-user system described previously, now with symbol misalignment. The received signal is now sampled twice, each sampling time being optimal to one of the transmitted signals in the sense that no own-signal ISI is present for the corresponding signal. Let the sequences of samples for the first signal be represented as {rka} and {yka}. Similarly, the sequences of samples for the second signal are {rkb} and {ykb}. The signals ρa(j) and ρb(j) are introduced to represent the pulse-shape response for the first and second signals, respectively. The received signals are then written as:

rka=Ckak+Gkuka+ωka,

rkb=Ckuk+Gkb+ωkb,

and:

yka=Akck+ΓKbgka+vka,

ykb=ΓKack+Bkgk+vkb.

Terms listed below are given by:





  
    
      u
      k
      a
    
    =
    
      
        ∑
        j
      
      ⁢
      
        
          
            ρ
            a
          
          ⁡
          
            (
            j
            )
          
        
        ⁢
        
          a
          
            k
            +
            j
          
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      u
      k
      b
    
    =
    
      
        ∑
        j
      
      ⁢
      
        
          
            ρ
            b
          
          ⁡
          
            (
            j
            )
          
        
        ⁢
        
          b
          
            k
            +
            j
          
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      Γ
      k
      b
    
    =
    
      
        ∑
        j
      
      ⁢
      
        
          
            ρ
            a
          
          ⁡
          
            (
            j
            )
          
        
        ⁢
        
          A
          
            k
            +
            j
          
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      Γ
      k
      a
    
    =
    
      
        ∑
        j
      
      ⁢
      
        
          
            ρ
            b
          
          ⁡
          
            (
            j
            )
          
        
        ⁢
        
          B
          
            k
            +
            j
          
        
      
    
  
  ,




where the summations above are taken over significant non-zero components of the pulse shape response terms.

In batch mode, detection is performed in the same manner as in the synchronous case. For each symbol time k, the detection approach is used to update âk and {circumflex over (b)}k. The ordering of k does not necessarily have to be in natural order {1,2, . . . , K}. After detection, channel estimation is performed using the approach described below. Because of the nature of the channel fading process, updates in natural order are assumed for the recursive channel estimation approach.
For detection, the terms uka and ukb are composed of prior, detected symbol values from the kth as well as other symbols from times j≠k. A flowchart for this approach is shown in FIG. 11. The form of this flowchart is similar to previous flowcharts, with initialization taking place at step 1100, iteration being handled by steps 1102 and 1103, and final values being output at 1104. Now, however, step 1101 is a process that performs the symbol update based upon the following equations:




  
    
      
        a
        ^
      
      k
    
    =
    
      
        1
        α
      
      ⁢
      
        
          C
          H
        
        ⁡
        
          (
          
            
              r
              k
              a
            
            -
            
              
                G
                k
              
              ⁢
              
                u
                k
                b
              
            
          
          )
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      
        b
        ^
      
      k
    
    =
    
      
        1
        α
      
      ⁢
      
        
          
            G
            H
          
          ⁡
          
            (
            
              
                r
                k
                b
              
              -
              
                
                  C
                  k
                
                ⁢
                
                  u
                  k
                  a
                
              
            
            )
          
        
        .
      
    
  




The interfering terms are composed of symbols from prior iterations and the pulse shape component. FIG. 12 shows a signal processing flow diagram of a symbol update unit that can be used to perform step 1101 in the flowchart of FIG. 11. The symbol update unit includes reformulation operators 1201 and 1202 and conjugation operators 1203 and 1204 as before. In this case however, calculation units 1205 and 1206 calculate values for the u terms. The interconnections between various other components such as adders and multipliers are shown in the figure.

A recursive channel estimator uses a form similar to that of the synchronous case, and is described by the equations:

ĉk=ĉk−1+μaAH(yka−Γkbgk),

ĝk=ĝk−1+μbBH(ykb−Γkack).

FIG. 13 is a flowchart for the recursive channel estimation. Stored values are initialized at step 1300. At step 1301, the update the update according to the above equations is performed. Steps 1302 and 1302 provide iteration. Once the iteration is complete, channel tracking is performed at 1304 using the step size term as input. Final values are output at step 1305.

FIG. 14 is a signal processing flow diagram for the update unit that performs the update at step 1301 in FIG. 13. This unit includes operators similar to those previously discussed, including reformulation operators 1401 and 1402 and conjugation operators 1403 and 1404.
A recursive approach can be used with transmit diversity where there is ISI. For the update of the symbols the same detection unit described above is used, but with the following changes. In performing the update of âk using received data sample yka, the interference terms ukb are made up of hypothesized and/or feedback values of the terms {tilde over (b)}j. From these hypothesized values, different estimates of âk and {circumflex over (b)}k result, and the values that lead to the best metric are detected. Once each symbol is detected at update time k, the recursive approach above is used to perform channel tracking. The detected symbols are used for channel tracking, along with any undetected (hypothesized) symbols that can be obtained from the chosen hypotheses.
In the case that conventional transmission is used with systems that also contain transmit diversity, the invention can still be used. Assume one signal, say ak, is a traditional transmission, and the interferer uses transmit diversity. In this case, assuming a1 and a2 are transmitted at consecutive symbol times corresponding to the transmit diversity symbol time k, then term yk replaces Ak and ck as follows:




  
    y
    k
  
  =
  
    
      
        [
        
          
            
              
                
                  a
                  1
                
              
            
            
              
                
                  a
                  2
                
              
            
          
          ⁢
          
            
              
                0
              
            
            
              
                0
              
            
          
        
        ]
      
      ⁡
      
        [
        
          
            
              
                c
                1
              
            
          
          
            
              0
            
          
        
        ]
      
    
    +
    
      
        B
        k
      
      ⁢
      
        g
        k
      
    
    +
    
      
        v
        k
      
      .
    
  




Replacing rk, Ck, and ak results in:





  
    r
    k
  
  =
  
    
      
        [
        
          
            
              
                
                  c
                  1
                
              
            
            
              
                0
              
            
          
          ⁢
          
            
              
                0
              
            
            
              
                
                  c
                  1
                  *
                
              
            
          
        
        ]
      
      ⁡
      
        [
        
          
            
              
                a
                1
              
            
          
          
            
              
                a
                2
                *
              
            
          
        
        ]
      
    
    +
    
      
        G
        k
      
      ⁢
      
        b
        k
      
    
    +
    
      
        ω
        k
      
      .
    
  




The above forms can be used for detection and channel estimation according to the invention in a mixed scenario.

It is also important to understand that the invention can be applied where there are more than two users. In such cases, the technique is applied to two users at a time. Such would be the case, for example, with a multi-input, multi-output (MIMO) system with many users.
FIG. 15 is a block diagram of a mobile terminal that implements the invention. FIG. 15 illustrates a terminal with voice capability, such as a mobile telephone. This illustration is an example only, and the invention works equally well with mobile terminals that are dedicated to communicating with text or other forms of data. Some types of mobile terminals may not contain all of the components shown in FIG. 15. As shown in FIG. 15, the terminal includes radio block 1501, a baseband logic block, 1502, control logic block 1503 and an audio interface block, 1504. Within radio block 1501, the receive and transmit information is converted from and to the radio frequencies (RF) of the various carrier types, and filtering is applied, as is understood in the art. Radio block 1501 includes the preprocessor previously discussed. The terminal's antenna system, 1507, is connected to the radio block. In baseband logic block 1502, basic signal processing occurs, e.g., synchronization, channel coding, decoding and burst formatting. In this example, the baseband logic includes a channel estimation unit, 1512, and a detection unit, 1513, according to the invention. The baseband logic block can be implemented by one or more ASIC'S, or perhaps by a digital signal processor (DSP). Audio interface block 1504 handles voice as well as analog-to-digital (A/D) and D/A processing. It also receives input through microphone 1505, and produces output through speaker 1506. Control logic block 1503, coordinates the aforedescribed blocks and also plays an important role in controlling the human interface components (not shown) such as a key pad and liquid crystal display (LCD). The functions of the aforedescribed transceiving blocks are directed and controlled by one or more microprocessors or digital signal processors such as main processor 1508, shown for illustrative purposes. Program code, often in the form of microcode is stored in memory 1509 and controls the operation of the terminal through the processor or processors. The mobile terminal illustrated in FIG. 15 interfaces to a smart card identity module (SIM), 1511, through a smart card reader interface. The interconnection between the main processor, control logic, memory, and SIM is depicted schematically. The interface is often an internal bus.
A mobile terminal implementation of the invention does not have to be a traditional “cellular telephone” type of terminal, but may include a cellular radio-telephone with or without a multi-line display; a personal communications system (PCS) terminal that may combine a cellular radiotelephone with data processing, facsimile and data communications capabilities; a personal data assistant (PDA) that can include a radiotelephone, pager, Internet/intranet access, Web browser, organizer; and a conventional laptop and/or palmtop computer or other appliance that includes a radiotelephone transceiver. Mobile terminals are sometimes also referred to as “pervasive computing” devices.
Specific embodiments of an invention are described herein. One of ordinary skill in the networking and signal processing arts will quickly recognize that the invention has other applications in other environments. In fact, many embodiments and implementations are possible. In addition, the recitation “means for” is intended to evoke a means-plus-function reading of an element in a claim, whereas, any elements that do not specifically use the recitation “means for,” are not intended to be read as means-plus-function elements, even if they otherwise include the word “means.” The following claims are in no way intended to limit the scope of the invention to the specific embodiments described.

## Claims

1. A method of demodulating a received signal corresponding to a plurality of transmitted symbol sequences originating from both a first user and a second user, the method comprising:
initializing, using known pilot symbols for the first user and the second user, a first-user symbol sequence corresponding to the first user and a second-user symbol sequence corresponding to the second user and channel estimates far the first user and the second user wherein the first user and the second user can be asynchronous; and
jointly detecting the first-user symbol sequence and the second-user symbol sequence using, at each of a plurality of iterations previous first and second user symbol sequences and the channel estimates for the first user and the second user, and based on the use of transmit diversity by at least one of the first user and the second user and on a pulse shape component and symbols from prior iterations.

2. The method of claim 1 further comprising updating the channel estimates prior to each of the plurality of iterations.

3. The method of claim 1 wherein channel estimates are obtained using alternate pilot channels in a code division multiplexed access (CDMA) system.

4. A method of demodulating a received signal corresponding to a plurality of transmitted symbols originating from both a first user and a second user, the method comprising:
initializing a pulse-shape component, using known pilot symbols for the first user and the second user, a first-user symbol corresponding to the first user and a second-user symbol corresponding to the second user and channel estimates for the first user and the second user wherein the first user and the second user can be asynchronous;
jointly detecting a first-user symbol and a second-user symbol at a current symbol time using previous first and second user symbols and channel estimates for the first user and the second user and the pulse-shape component, and based on the use of transmit diversity by at least one of the first user and the second user; and
if additional symbol times remain, jointly detecting the first user symbol and the second user symbol at a next symbol time.

5. The method of claim 4 further comprising:
predicting the channel estimates for a next symbol time, at least in part through the use of channel tracking; and
updating the first-user symbol and the second-user symbol.

6. The method of claim 4 where the jointly detecting can be accomplished at least in part by estimating the plurality of symbols out of natural order.

7. Apparatus for demodulating a received signal corresponding to a plurality of transmitted symbol sequences originating from both a first user and a second user, the apparatus comprising:
means for initializing, using known pilot symbols for the first user and the second user, a first-user symbol sequence corresponding to the first user and a second-user symbol sequence corresponding to the second user wherein the first user and the second user can be asynchronous;
means for initializing channel estimates for the first user and the second user, using the known pilot symbols for the first user and the second user; and
means for iteratively and jointly detecting the first-user symbol sequence and the second-user symbol sequence using previous first and second user symbol sequences and the channel estimates for the first user and the second user, and based on the use of transmit diversity by at least one of the first user and the second user and on a pulse shape component and symbols from prior iterations.

8. A receiving system operable to demodulate a received signal corresponding to a plurality of transmitted symbol sequences originating from both a first user and a second user, the system comprising:
a channel estimation unit operable to produce channel estimates for the first user and the second user, using known pilot symbols for the first user and the second user, and
a detector unit operatively connected to the channel estimation unit, the detector unit operable to jointly detect a first-user symbol sequence and a second-user symbol sequence using previous first and second user symbol sequences, the known pilot symbols for the first user and the second user, and the channel estimates for the first user and the second user, wherein the first user and the second user can be asynchronous, and based on the use of transmit diversity by at least one of the first user and the second user and on a pulse shape component and symbols from prior iterations.

9. The system of claim 8 wherein the channel estimation unit updates the channel estimates prior to each of a plurality of iterations.

10. The system of claim 8 wherein the channel estimation unit obtains channel estimates using alternate pilot channels in a code division multiplexed access (CDMA) system.

11. A receiving system operable to demodulate a received signal corresponding to a plurality of transmitted symbols originating from both a first user and a second user, the system comprising:
a channel estimation unit operable to produce channel estimates for the first user and the second user, using known pilot symbols for the first user and the second user; and
a detector unit including a symbol update unit, the detector unit operatively connected to the channel estimation unit, the detector unit operable to jointly detect a first-user symbol and a second-user symbol at specific symbol times using previous first and second user symbols, the known pilot symbols for the first user and the second user, and the channel estimates for the first user and the second user, as well as a pulse-shape component and symbols from prior iteration, wherein the first user and the second user can be asynchronous, and based on the use of transmit diversity by at least one of the first user and the second user;
and wherein both the symbol update unit and channel estimation unit further comprise reformulation operators, conjugation operators connected to the reformulation operators, and multipliers and adders interconnected with the reformulation operators and conjugation operators.

12. The system of claim 11 wherein the channel estimation logic predicts channel estimates at least in part through channel tracking.

13. The receiving system of claim 11 wherein the detector unit is further operable to detect the first user symbol and the second user symbol out of natural order.

14. A mobile terminal comprising:
a processor system for controlling the operation of the mobile terminal;
a radio block operable to transmit signals, and also to receive signals corresponding to a plurality of transmitted symbols originating from both a first user and a second user;
baseband logic operatively connected to the radio block and the processor system, the baseband logic further comprising a channel estimation unit operable to produce channel estimates for the first user and the second user, using known pilot symbols for the first user and the second user, and a detector unit operatively connected to the channel estimation unit, the detector unit operable to jointly detect first-user symbols and second-user symbols based on the channel estimates for the first user and the second user, and using the known pilot symbols for the first user and the second user, wherein the first user and the second user can be asynchronous, and based on the use of transmit diversity by at least one of the first user and the second user and on a pulse shape component and symbols from prior iterations.

15. The mobile terminal of claim 14 wherein the channel estimation unit performs block estimation and the detector unit performs block detection of the first user symbols and the second user symbols.

16. The mobile terminal of claim 14 wherein the channel estimation unit performs recursive estimation using channel tracking and the detector unit performs recursive detection of the first user symbols and the second user symbols.

17. The mobile terminal of claim 16 wherein the detector unit further comprises a symbol update unit, and the symbol update unit and the channel estimation unit each further comprise:
reformulation operators;
conjugation operators connected to the reformulation operators; and
multipliers and adders interconnected with the reformulation operators and conjugation operators.

18. The mobile terminal of claim 14 wherein the channel estimation unit produces channel estimates by referring to alternate pilot channels in a code division multiplexed access (CDMA) system.

19. The mobile terminal of claim 14 wherein the detector unit further comprises a symbol update unit, and the symbol update unit and the channel estimation unit each further comprise:
reformulation operators;
conjugation operators connected to the reformulation operators; and
multipliers and adders interconnected with the reformulation operators and conjugation operators.

20. The mobile terminal of claim 14 wherein the detector unit is further operable to detect the first user symbol and the second user symbol out of natural order.