# Patent

## Title
Method for estimating time of arrival of received signals for ultra wide band impulse radio systems

## Abstract

A method estimates a time of arrival of a signal received in a wireless communication system. An energy in a frame of a received signal is measured to determine a block in the frame, the block representing a coarse time of arrival of the received signal. Multiple time-delayed versions of a template signal are combined with the block of the received signal to select a particular template signal. The particular template signal identifies a particular chip representing a fine time of arrival of the received signal.

## Background
FIELD OF THE INVENTION
The present invention relates generally to radio communication systems, and more particularly to determining a time of arrival of a received signal in a wireless communications network.
BACKGROUND OF THE INVENTION
To estimate a distance between a transmitter and a receiver in a wireless communications network, the transmitter can send a signal to the receiver at a time t1. The receiver, as soon as possible, returns a reply signal to the transmitter. The transmitter measures the time of arrival (TOA) of the reply signal at time t2. An estimate of the distance between the transmitter and the receiver is the time for the signal to make the round trip divided by two and multiplying by the speed of light, i.e.:




  D
  =
  
    
      
        
        
          
            t
            1
          
          -
          
            t
            2
          
        
        
      
      2
    
    ⁢
    
      c
      .
    
  




Accurate time resolution of ultra wideband (UWB) signals facilitates very precise positioning capabilities based on signal TOA measurements. Although a theoretical lower bound for TOA estimation can be achieved by using maximum likelihood methods, those prior art methods are impractical due to a need for extremely high sampling rates and a large number of multipath components of the signal.
Another method is correlation-based. That method serially searches possible signal delays of a signal received via a first signal path and takes a very long time to estimate the TOA of the received signal.
Moreover, the signal received from the first path does not always have a strongest correlation output, which can result in an inaccurate TOA estimate by the prior art correlation-based methods.
Therefore, there is a need for a time of arrival estimation method that overcomes the problems of the prior art.
SUMMARY OF THE INVENTION
The invention provides a method for estimating a time of arrival of a signal received in a wireless communication system. An energy in a frame of a received signal is measured to determine a block in the frame, the block representing a coarse time of arrival of the received signal. Time-delayed versions of a template signal are combined with the block of the received signal to select a particular template signal. The particular template signal identifies a particular chip representing a fine time of arrival of the received signal.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram of frame-based time intervals;
FIG. 2 is a flow diagram of the method for estimating a time of arrival according to the invention;
FIG. 3 is a block diagram of determining a coarse time of arrival of a received signal according to the invention; and
FIG. 4 is a flow diagram for determining a fine time of arrival of a received signal according to the invention.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
System Structure and Method Operation
Our invention provides a method for estimating a time of arrival (TOA) of a signal at a radio transceiver in a wireless communications network. For the purpose of this description, the transceiver estimates the TOA for a received signal. However, it should be understood that the transceiver can transmit and receive.
As shown in FIG. 1, wireless impulse radio transceivers mark time in terms of frames 101, blocks 102, and chips 103. Frames are longer than blocks, which are longer than chips. Each frame includes multiple blocks. Each block includes multiple chips.
A binary phase shift keying (BPSK) time hopping (TH) impulse radio (IR) transmitted signal can be represented by:




  
    
      
        
          
            
              s
              tx
            
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              E
            
            ⁢
            
              
                ∑
                
                  j
                  =
                  
                    -
                    ∞
                  
                
                ∞
              
              ⁢
              
                
              
              ⁢
              
                
                  a
                  j
                
                ⁢
                
                  b
                  
                    ⌊
                    
                      j
                      /
                      
                        N
                        f
                      
                    
                    ⌋
                  
                
                ⁢
                
                  
                    w
                    tx
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        jT
                        f
                      
                      -
                      
                        
                          c
                          j
                        
                        ⁢
                        
                          T
                          c
                        
                      
                    
                    )
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        1
        )
      
    
  




where wtx(t) is a transmitted ultra wideband (UWB) pulse having a duration Tc, E is a transmitted pulse energy, Tf is a frame time, Nf is a number of pulses representing one information symbol, Ts=NfTf is a symbol interval, and b└j/Nf┘ε{+1,−1} is a binary information symbol. In order to smooth a power spectrum of the transmitted signal and to allow the channel to be shared by multiple users without causing collisions, a time-hopping (TH) sequence cjε{0, 1, . . . , Nc−1} is assigned to each transmitter, where Nc is the number of chips per frame interval, that is, Nc=Tf=Tc.
Additionally, random polarity codes, aj's, can be included. The polarity codes are binary random variables taking values ±1 with equal probability. The values are known at the receiver.
Consider the following channel model




  
    
      
        
          
            r
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                ∑
                
                  l
                  =
                  1
                
                L
              
              ⁢
              
                
              
              ⁢
              
                
                  E
                
                ⁢
                
                  α
                  l
                
                ⁢
                
                  
                    s
                    rx
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        
                          (
                          
                            l
                            -
                            1
                          
                          )
                        
                        ⁢
                        
                          T
                          c
                        
                      
                      -
                      
                        τ
                        TOA
                      
                    
                    )
                  
                
              
            
            +
            
              n
              ⁡
              
                (
                t
                )
              
            
          
        
        ,
      
    
    
      
        (
        2
        )
      
    
  




where al is a channel coefficient for an lth path, L is a number of multipath components, and τTOA is a TOA of the received signal. From equations (1) and (2), and considering effects of antennas, the received signal can be represented by:




  
    
      
        
          
            r
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                ∑
                
                  l
                  =
                  1
                
                L
              
              ⁢
              
                
              
              ⁢
              
                
                  E
                
                ⁢
                
                  α
                  l
                
                ⁢
                
                  
                    s
                    rx
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        
                          (
                          
                            l
                            -
                            1
                          
                          )
                        
                        ⁢
                        
                          T
                          c
                        
                      
                      -
                      
                        τ
                        TOA
                      
                    
                    )
                  
                
              
            
            +
            
              n
              ⁡
              
                (
                t
                )
              
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  




where srx(t) is given by




  
    
      
        
          
            
              s
              rx
            
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              ∑
              
                j
                =
                
                  -
                  ∞
                
              
              ∞
            
            ⁢
            
              
            
            ⁢
            
              
                a
                j
              
              ⁢
              
                b
                
                  ⌊
                  
                    j
                    /
                    
                      N
                      f
                    
                  
                  ⌋
                
              
              ⁢
              
                
                  w
                  rx
                
                ⁡
                
                  (
                  
                    t
                    -
                    
                      jT
                      f
                    
                    -
                    
                      
                        c
                        j
                      
                      ⁢
                      
                        T
                        c
                      
                    
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        4
        )
      
    
  




with wrx(t) denoting the received UWB pulse with unit energy. Assuming a data aided TOA estimation method using a training sequence, we consider a training sequence of bj=1∀j.
In this case, equation (4) can be expressed as




  
    
      
        
          
            s
            rx
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              j
              =
              
                -
                ∞
              
            
            ∞
          
          ⁢
          
            
              a
              j
            
            ⁢
            
              
                w
                rx
              
              ⁡
              
                (
                
                  t
                  -
                  
                    jT
                    f
                  
                  -
                  
                    
                      c
                      j
                    
                    ⁢
                    
                      T
                      c
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




For simplicity, we assume that the signal always arrives in one frame duration, i.e., τTOA<Tf, and there is no inter-frame interference (IFI); that is, Tf≧(L+cmax)Tc or, equivalently, Nc≧L+cmax, where cmax is a maximum value of the TH sequence. Note that the assumption of τTOA<Tf does not restrict the invention. In fact, it is enough to have τTOA<Ts for the invention to work when the frame is sufficiently large and predetermined TH codes are used.
Moreover, even if τTOA<Ts, an initial energy detection can be used to determine the arrival time within a symbol uncertainty.
Two Step TOA Estimation Method
One of the most challenging tasks in time of arrival estimation is to obtain a reliable TOA estimate in an acceptable time interval under the constraint of sampling rate. In order to have a low power and low complexity receiver, we use symbol-rate sampling in our preferred embodiment.
The invention provides a method for estimating a TOA that can perform TOA estimation from symbol-rate samples in less time than prior art methods, and at chip-level resolution.
As shown in FIG. 2, the invention estimates 200 a TOA of a received signal 201 at a particular chip 202 of a particular block 301 in a frame.
A first step according to the invention determines 300 a block representing a coarse TOA 301 of the received signal 201 based on a signal measurement of energy of the received signal.
As shown in FIG. 3, the signal 201 received during a frame 101 with Nb blocks 102 has a highest measured energy during a third block 301. Therefore, the coarse time of arrival is determined to be during the third block 301.
In a second step, a fine time of arrival 202 of the received signal is estimated by applying a change detection method 400, which combines multiple, time-delayed versions of a template signal to the received signal to identify a fine time of arrival 202. The template signals are transmitted signals corresponding to the received signal. In the preferred embodiment, the time delays are one-chip intervals.
FIG. 4 illustrates the change detection method 400 for determining the fine TOA 202 according to the invention.
As shown in FIG. 4, after the third block 301 is determined, the received signal is combined 410 with each of multiple, time-delayed versions 406 of the template signal 405 according to the block 301. The combining 410 produces a set of combined values 411, one for each combination of a time-delayed template with the received signal. The template signal associated with a combined value that matches a correlation value 415 is selected 420 to identify the chip that represents the fine TOA 202.
We express the TOA (τTOA) in equation (3) as follows:

τTOA=kTc=kbTb+kcTc  (6)

where kε[0,Nc−1] is the TOA in terms of the chip interval Tc, Tb is the block interval including B chips (Tb=BTc), and kbε[0,Nc/B−1] and kcε[0,B−1] are integers that determine, respectively, in which block and chip the signal arriving via the first signal path.
For simplicity, the TOA is assumed to be an integer multiple of the chip duration Tc. In a practical application, sub-chip resolution can be obtained by employing a delay-lock-loop (DLL) after the TOA estimation with chip-level uncertainty.
Coarse TOA Estimation
As shown in FIG. 3, the coarse arrival time of the signal in the frame interval is determined 300, i.e., which block in the frame represents the arrival time of the received signal, e.g. a UWB pulse. Without loss of generality, we assume that the frame time Tf is an integer multiple Nb of block size Tb, i.e., Tf=NbTb. In order to have reliable decision variables in this step, energies from N1 different frames of the incoming signal for each block can be combined. Hence, the decision variables are expressed as




  
    
      
        
          
            Y
            i
          
          =
          
            
              ∑
              
                j
                =
                0
              
              
                
                  N
                  1
                
                -
                1
              
            
            ⁢
            
              Y
              
                i
                ,
                j
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            for
            ⁢
            
              
            
            ⁢
            i
          
          =
          0
        
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            N
            b
          
          -
          1
        
        ,
        
          
        
        ⁢
        where
      
    
    
      
        (
        7
        )
      
    
  
  
    
      
        
          Y
          
            i
            ,
            j
          
        
        =
        
          
            ∫
            
              
                jT
                f
              
              +
              
                iT
                b
              
              +
              
                
                  c
                  j
                
                ⁢
                
                  T
                  c
                
              
            
            
              
                jT
                f
              
              +
              
                
                  (
                  
                    i
                    +
                    1
                  
                  )
                
                ⁢
                
                  T
                  b
                
              
              +
              
                
                  c
                  j
                
                ⁢
                
                  T
                  c
                
              
            
          
          ⁢
          
            
              
                
                
                  r
                  ⁡
                  
                    (
                    t
                    )
                  
                
                
              
              2
            
            ⁢
            
              
            
            ⁢
            
              
                d
                t
              
              .
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




Then, kb in (6) is estimated as




  
    
      
        
          
        
        ⁢
        
          
            
              k
              ^
            
            b
          
          =
          
            arg
            ⁢
            
              
                
                  
                
                ⁢
                max
              
              
                0
                ≤
                i
                ≤
                
                  
                    N
                    b
                  
                  -
                  1
                
              
            
            ⁢
            
              
                Y
                i
              
              .
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




In other words, we select the block with the largest signal energy.
The parameters of this step that can be optimized are the block size Tb(Nb) and the number of frames N1, from which energy is measured.
Fine TOA Estimation from Low Rate Correlation Outputs
After determining the coarse arrival time, the second step estimates the fine TOA kc according to equation (6). Ideally, chips kcε[0,B−1] need to be searched for fine TOA estimation, which corresponds to searching kε[{circumflex over (k)}bB,({circumflex over (k)}b+1)B−1], with chips {circumflex over (k)}b determined from equation (9).
However, in some cases, the first signal path can be in one of the blocks preceding the block having the highest energy level due to multipath effects. Therefore, instead of searching a single block, multiple blocks kε[{circumflex over (k)}bB−M1,({circumflex over (k)}b+1)B−1], with M1≧0, can be searched for the fine TOA in order to increase the probability of detecting of a first path of the received signal. In other words, in addition to the block with the largest signal energy, we perform an additional search over M1 chips by applying templates with relatively long time delays.
For notational simplicity, let U={ns, ns+1, . . . , ne} denote an uncertainty region, i.e., possible TOA of the first path of the received signal, where ns={circumflex over (k)}bB−M1 and ne=({circumflex over (k)}b+1)B−1 are the start and end points of the uncertainty region in the frame.
In order to estimate the fine TOA, i.e., a TOA at a chip-level resolution, we consider combinations of the received signal with time-shifted versions 406 of the template signal 405. For delay iTc, we obtain the following output:




  
    
      
        
          
            z
            i
          
          =
          
            
              ∫
              
                iT
                c
              
              
                
                  iT
                  c
                
                +
                
                  
                    N
                    2
                  
                  ⁢
                  
                    T
                    f
                  
                
              
            
            ⁢
            
              
                r
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                
                  s
                  temp
                
                ⁡
                
                  (
                  
                    t
                    -
                    
                      iT
                      c
                    
                  
                  )
                
              
              ⁢
              
                
              
              ⁢
              
                d
                t
              
            
          
        
        ,
      
    
    
      
        (
        10
        )
      
    
  




where N2 is the number of frames over which the combination output is obtained, and stemp(t) is the template signal given by




  
    
      
        
          
            s
            temp
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              j
              =
              0
            
            
              
                N
                2
              
              -
              1
            
          
          ⁢
          
            
              a
              j
            
            ⁢
            
              
                w
                ⁡
                
                  (
                  
                    t
                    -
                    
                      jT
                      f
                    
                    -
                    
                      
                        c
                        j
                      
                      ⁢
                      
                        T
                        c
                      
                    
                  
                  )
                
              
              .
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




From the combination outputs for different delays, the object is to determine the chip in which the first signal path has arrived. According to the block interval Tb and considering the multipath components in the received signal, which is typical for indoor UWB systems, we can assume that the block starts with a number of chips with noise-only components and the remaining chips start with signal plus noise components. Assuming that the statistics of the signal paths do not change in the uncertainty region, we can express different hypotheses approximately as follows:

H0:zi=ηi,i=ns, . . . , nf,

Hk:zi=ηi,i=ns, . . . , k−1,

zi=N2√{square root over (E)}αi−k+1+ηi,i=k, . . . ,nf,  (12)

for kεU, where ηn's denote the i.i.d. output noise distributed as N(0, αn2), and σn2=N2N0/2, α1, . . . , αnf−k+1 are independent channel coefficients, assuming nf−ns+1≦L, and nf=ne+M2, with M2 being the number of combination outputs that are considered out of the uncertainty region in order to have reliable estimates of the unknown parameters of α.
Due to very high resolution of UWB signals, it is appropriate to model the channel coefficients approximately as




  
    
      
        
          
            α
            1
          
          =
          
            
              d
              1
            
            ⁢
            
              
              
                α
                1
              
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            α
            l
          
          =
          
            {
            
              
                
                  
                    
                      
                        
                          d
                          l
                        
                        ⁢
                        
                          
                          
                            α
                            l
                          
                          
                        
                      
                      ,
                    
                  
                  
                    p
                  
                
                
                  
                    
                      0
                      ,
                    
                  
                  
                    
                      1
                      -
                      p
                    
                  
                
              
              ,
              
                l
                =
                2
              
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  n
                  f
                
                -
                
                  n
                  s
                
                +
                1
              
              ,
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




where p is a probability that a channel tap arrives in a given chip, dl is the phase of αl, which is ±1 with equal probability, and |αl| is the amplitude of αl, which is modeled as a Nakagami-m distributed random variable with parameter Ω; that is,




  
    
      
        
          
            p
            ⁡
            
              (
              α
              )
            
          
          =
          
            
              2
              
                Γ
                ⁡
                
                  (
                  m
                  )
                
              
            
            ⁢
            
              
                (
                
                  m
                  Ω
                
                )
              
              m
            
            ⁢
            
              α
              
                
                  2
                  ⁢
                  m
                
                -
                1
              
            
            ⁢
            
              e
              
                -
                
                  
                    mα
                    2
                  
                  Ω
                
              
            
          
        
        ,
      
    
    
      
        (
        14
        )
      
    
  




for α≧0, m≧0.5 and Ω≧0, where Γ is the Gamma function.
According to equation (12), the TOA estimation problem can be considered as a change detection problem. Let θ denote the unknown parameters of the distribution of α; that is, θ=[p m α]. Then, the log-likelihood ratio (LLR) is determined as:




  
    
      
        
          
            
              S
              k
              
                n
                f
              
            
            ⁡
            
              (
              θ
              )
            
          
          =
          
            
              ∑
              
                i
                =
                k
              
              
                n
                f
              
            
            ⁢
            
              log
              ⁢
              
                
                  
                    p
                    θ
                  
                  ⁡
                  
                    (
                    
                      
                        z
                        i
                      
                      ⁢
                      
                        ❘
                      
                      ⁢
                      
                        H
                        k
                      
                    
                    )
                  
                
                
                  p
                  ⁡
                  
                    (
                    
                      
                        z
                        i
                      
                      ⁢
                      
                        ❘
                      
                      ⁢
                      
                        H
                        0
                      
                    
                    )
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        15
        )
      
    
  




where pθ(zi|Hk) denotes the probability distribution function (p.d.f) of the correlation output under hypothesis Hk, with unknown parameters given by θ, and p(zi|H0) denotes the p.d.f. of the correlation output under hypothesis H0. Because θ is unknown, the maximum likelihood (ML) estimate can be obtained first for a given hypothesis Hk and then that estimate can be used in the LLR expression. In other words, the generalized LLR approach can be taken, where the TOA estimate is expressed as




  
    
      
        
          
            k
            ^
          
          =
          
            arg
            ⁢
            
              max
              
                k
                ∈
                U
              
            
            ⁢
            
              
                S
                k
                
                  n
                  f
                
              
              (
              
                
                  
                    θ
                    ^
                  
                  ML
                
                ⁡
                
                  (
                  k
                  )
                
              
              )
            
          
        
        ,
        
          

        
        ⁢
        where
      
    
    
      
        (
        16
        )
      
    
  
  
    
      
        
          
            
              θ
              ^
            
            ML
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          arg
          ⁢
          
            
          
          ⁢
          
            sup
            θ
          
          ⁢
          
            
              
                S
                k
                
                  n
                  f
                
              
              ⁡
              
                (
                θ
                )
              
            
            .
          
        
      
    
    
      
        (
        17
        )
      
    
  




However, the ML estimate is usually complicated. Therefore, simpler estimators, such as a method of moments (MM) estimator can be used to obtain the parameters. The nth moment of a random variable X having the Nakagami-m distribution with parameter Ω is given by




  
    
      
        
          E
          ⁢
          
            {
            
              X
              n
            
            }
          
        
        =
        
          
            
              Γ
              ⁡
              
                (
                
                  m
                  +
                  
                    n
                    /
                    2
                  
                
                )
              
            
            
              Γ
              ⁡
              
                (
                m
                )
              
            
          
          ⁢
          
            
              
                (
                
                  Ω
                  m
                
                )
              
              
                n
                /
                2
              
            
            .
          
        
      
    
    
      
        (
        18
        )
      
    
  




Then, from the correlator outputs {zi}inf=k+1, the MM estimates for the unknown parameters can be determined by:




  
    
      
        
          
            p
            MM
          
          =
          
            
              
                γ
                1
              
              ⁢
              
                γ
                2
              
            
            
              
                2
                ⁢
                
                  
                
                ⁢
                
                  γ
                  2
                  2
                
              
              -
              
                γ
                3
              
            
          
        
        ,
        
          
            m
            MM
          
          =
          
            
              
                2
                ⁢
                
                  
                
                ⁢
                
                  γ
                  2
                  2
                
              
              -
              
                γ
                3
              
            
            
              
                γ
                3
              
              -
              
                γ
                2
                2
              
            
          
        
        ,
        
          
            Ω
            MM
          
          =
          
            
              
                2
                ⁢
                
                  γ
                  2
                  2
                
              
              -
              
                γ
                3
              
            
            
              γ
              2
            
          
        
        ,
        
          

        
        ⁢
        where
      
    
    
      
        (
        19
        )
      
    
  
  
    
      
        
          
            γ
            1
          
          ⁢
          
            =
            △
          
          ⁢
          
            
              1
              
                EN
                2
                2
              
            
            ⁢
            
              (
              
                
                  μ
                  2
                
                -
                
                  σ
                  n
                  2
                
              
              )
            
          
        
        ,
        
          

        
        ⁢
        
          
            γ
            2
          
          ⁢
          
            =
            △
          
          ⁢
          
            
              1
              
                
                  E
                  2
                
                ⁢
                
                  N
                  2
                  4
                
              
            
            ⁢
            
              (
              
                
                  
                    
                      μ
                      4
                    
                    -
                    
                      3
                      ⁢
                      
                        
                      
                      ⁢
                      
                        σ
                        n
                        2
                      
                    
                  
                  
                    γ
                    1
                  
                
                -
                
                  6
                  ⁢
                  
                    EN
                    2
                    2
                  
                  ⁢
                  
                    σ
                    n
                    2
                  
                
              
              )
            
          
        
        ,
        
          

        
        ⁢
        
          
            γ
            3
          
          ⁢
          
            =
            △
          
          ⁢
          
            
              1
              
                
                  E
                  3
                
                ⁢
                
                  N
                  2
                  6
                
              
            
            ⁢
            
              (
              
                
                  
                    
                      μ
                      6
                    
                    -
                    
                      15
                      ⁢
                      
                        
                      
                      ⁢
                      
                        σ
                        n
                        6
                      
                    
                  
                  
                    γ
                    1
                  
                
                -
                
                  15
                  ⁢
                  
                    
                  
                  ⁢
                  
                    E
                    2
                  
                  ⁢
                  
                    N
                    2
                    4
                  
                  ⁢
                  
                    γ
                    2
                  
                  ⁢
                  
                    σ
                    n
                    2
                  
                
                -
                
                  45
                  ⁢
                  
                    EN
                    2
                    2
                  
                  ⁢
                  
                    σ
                    n
                    4
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        20
        )
      
    
  




with μj denoting the jth sample moment given by




  
    
      
        
          μ
          j
        
        =
        
          
            1
            
              
                n
                f
              
              -
              k
            
          
          ⁢
          
            
              ∑
              
                i
                =
                
                  k
                  +
                  1
                
              
              
                n
                f
              
            
            ⁢
            
              
                z
                i
                j
              
              .
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




Then, the chip having the first signal path can be obtained as




  
    
      
        
          
            k
            ^
          
          =
          
            arg
            ⁢
            
              
            
            ⁢
            
              
                max
                
                  k
                  ∈
                  U
                
              
              ⁢
              
                
                  S
                  k
                  
                    n
                    f
                  
                
                (
                
                  
                    
                      θ
                      ^
                    
                    MM
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        22
        )
      
    
  




where θMM(k)=[pMM mMMΩMM] is the MM estimate for the unknown parameters.
Let p1(z) and p2(z), respectively, denote the distributions of η and N2√{square root over (E)}d|α|+η. Then, the generalized LLR for the kth hypothesis is given by




  
    
      
        
          
            
              S
              k
              
                n
                f
              
            
            ⁡
            
              (
              
                θ
                ^
              
              )
            
          
          =
          
            
              log
              ⁢
              
                
                  
                    p
                    2
                  
                  ⁡
                  
                    (
                    
                      z
                      k
                    
                    )
                  
                
                
                  
                    p
                    1
                  
                  ⁡
                  
                    (
                    
                      z
                      k
                    
                    )
                  
                
              
            
            +
            
              
                ∑
                
                  i
                  =
                  
                    k
                    +
                    1
                  
                
                
                  n
                  f
                
              
              ⁢
              
                log
                ⁢
                
                  
                    
                      
                        pp
                        2
                      
                      ⁡
                      
                        (
                        
                          z
                          i
                        
                        )
                      
                    
                    +
                    
                      
                        (
                        
                          1
                          -
                          p
                        
                        )
                      
                      ⁢
                      
                        
                          p
                          1
                        
                        ⁡
                        
                          (
                          
                            z
                            i
                          
                          )
                        
                      
                    
                  
                  
                    
                      p
                      1
                    
                    ⁡
                    
                      (
                      
                        z
                        i
                      
                      )
                    
                  
                
              
            
          
        
        ,
        
          

        
        ⁢
        where
      
    
    
      
        (
        23
        )
      
    
  
  
    
      
        
          
            
              p
              1
            
            ⁡
            
              (
              z
              )
            
          
          =
          
            
              1
              
                
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    π
                  
                
                ⁢
                
                  
                
                ⁢
                
                  σ
                  n
                
              
            
            ⁢
            
              e
              
                -
                
                  
                    z
                    2
                  
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      σ
                      n
                      2
                    
                  
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        24
        )
      
    
  
  
    
      
        
          
            
              p
              2
            
            ⁡
            
              (
              z
              )
            
          
          =
          
            
              
                v
                1
              
              
                
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    π
                  
                
                ⁢
                
                  
                
                ⁢
                
                  σ
                  n
                
              
            
            ⁢
            
              e
              
                -
                
                  
                    z
                    2
                  
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      σ
                      n
                      2
                    
                  
                
              
            
            ⁢
            
              Φ
              (
              
                m
                ,
                
                  
                    1
                    2
                  
                  ;
                  
                    
                      z
                      2
                    
                    
                      v
                      2
                    
                  
                
              
              )
            
          
        
        ,
        
          

        
        ⁢
        with
      
    
    
      
        (
        25
        )
      
    
  
  
    
      
        
          
            v
            1
          
          ⁢
          
            =
            △
          
          ⁢
          
            
              
                2
                ⁢
                
                  π
                
                ⁢
                
                  Γ
                  ⁡
                  
                    (
                    
                      2
                      ⁢
                      m
                    
                    )
                  
                
              
              
                
                  Γ
                  ⁡
                  
                    (
                    m
                    )
                  
                
                ⁢
                
                  Γ
                  ⁡
                  
                    (
                    
                      m
                      +
                      0.5
                    
                    )
                  
                
              
            
            ⁢
            
              
                (
                
                  4
                  +
                  
                    
                      2
                      ⁢
                      
                        EN
                        2
                        2
                      
                      ⁢
                      Ω
                    
                    
                      m
                      ⁢
                      
                        
                      
                      ⁢
                      
                        σ
                        n
                        2
                      
                    
                  
                
                )
              
              
                -
                m
              
            
          
        
        ,
        
          

        
        ⁢
        
          
            v
            2
          
          ⁢
          
            =
            △
          
          ⁢
          
            2
            ⁢
            
              
            
            ⁢
            
              
                σ
                n
                2
              
              (
              
                1
                +
                
                  2
                  ⁢
                  m
                  ⁢
                  
                    
                      σ
                      n
                      2
                    
                    
                      
                        EN
                        2
                        2
                      
                      ⁢
                      Ω
                    
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        26
        )
      
    
  




and Φ denoting a confluent hypergeometric function given by [7]




  
    
      
        
          Φ
          ⁡
          
            (
            
              
                β
                1
              
              ,
              
                
                  β
                  2
                
                ;
                x
              
            
            )
          
        
        =
        
          1
          +
          
            
              
                β
                1
              
              
                β
                2
              
            
            ⁢
            
              x
              
                1
                !
              
            
          
          +
          
            
              
                
                  β
                  1
                
                ⁡
                
                  (
                  
                    
                      β
                      1
                    
                    +
                    1
                  
                  )
                
              
              
                
                  β
                  2
                
                ⁡
                
                  (
                  
                    
                      β
                      2
                    
                    +
                    1
                  
                  )
                
              
            
            ⁢
            
              
                x
                2
              
              
                2
                !
              
            
          
          +
          
            
              
                
                  
                    β
                    1
                  
                  ⁡
                  
                    (
                    
                      
                        β
                        1
                      
                      +
                      1
                    
                    )
                  
                
                ⁢
                
                  (
                  
                    
                      β
                      1
                    
                    +
                    2
                  
                  )
                
              
              
                
                  
                    β
                    2
                  
                  ⁡
                  
                    (
                    
                      
                        β
                        2
                      
                      +
                      1
                    
                    )
                  
                
                ⁢
                
                  (
                  
                    
                      β
                      2
                    
                    +
                    2
                  
                  )
                
              
            
            ⁢
            
              
                x
                3
              
              
                3
                !
              
            
          
          +
          ⋯
        
      
    
    
      
        (
        27
        )
      
    
  




Note that the p.d.f. of N2√{square root over (E)}d|α|+η,p2(z), is obtained from equations (14), (24) and the fact that d is ±1 with equal probability.
The TOA estimation rule can be expressed as




  
    
      
        
          k
          ^
        
        =
        
          arg
          ⁢
          
            
          
          ⁢
          
            
              max
              
                k
                ∈
                u
              
            
            ⁢
            
              
                {
                
                  
                    log
                    ⁡
                    
                      [
                      
                        
                          v
                          1
                        
                        ⁢
                        
                          Φ
                          ⁡
                          
                            (
                            
                              m
                              ,
                              
                                0.5
                                ;
                                
                                  
                                    z
                                    k
                                    2
                                  
                                  
                                    v
                                    2
                                  
                                
                              
                            
                            )
                          
                        
                      
                      ]
                    
                  
                  +
                  
                    
                      ∑
                      
                        n
                        =
                        
                          k
                          +
                          1
                        
                      
                      
                        n
                        f
                      
                    
                    ⁢
                    
                      log
                      ⁡
                      
                        [
                        
                          
                            
                              pv
                              1
                            
                            ⁢
                            
                              Φ
                              ⁡
                              
                                (
                                
                                  m
                                  ,
                                  
                                    0.5
                                    ;
                                    
                                      
                                        z
                                        k
                                        2
                                      
                                      
                                        v
                                        2
                                      
                                    
                                  
                                
                                )
                              
                            
                          
                          +
                          1
                          -
                          p
                        
                        ]
                      
                    
                  
                
                }
              
              .
            
          
        
      
    
    
      
        (
        28
        )
      
    
  




Note that equation (12) assumes that the block always starts with noise-only components, followed by the arriving signal. However, in practice, there can be cases where the first step determines a block of all noise components. By combining large number of frames; that is, by choosing a large N1 in equation (7), the probability of selecting a block 301 having only noise can be reduced. However, a large N1 also increases the estimation time. Hence, there is a trade-off between the estimation error and the estimation time.
In order to prevent erroneous TOA estimation when a noise-only block is selected, a one-sided test can be applied using the known distribution of the noise outputs. Because the noise outputs have a Gaussian distribution, the test compares the average energy of the outputs after the estimated change instant to a threshold.
In other words, if




  
    
      
        1
        
          
            n
            f
          
          -
          
            k
            ^
          
          +
          1
        
      
      ⁢
      
        
          ∑
          
            i
            =
            
              k
              ^
            
          
          
            n
            f
          
        
        ⁢
        
          z
          i
          2
        
      
    
    <
    
      δ
      1
    
  
  ,




then the block is considered a noise-only block, and the two-step process is repeated.

Another improvement of the invention can be achieved by checking whether the block 301 includes signal from all paths; that is, determining if the fine TOA is actually prior to block 301. Again, by following a one-sided test approach, we can check the average energy of the correlation outputs before the estimated TOA against a threshold and detect an all-signal block if the threshold is exceeded.
However, for very small values of the TOA estimate {circumflex over (k)}, there can be a significant probability that the signal from the first path arrives before the current observation region because the distribution of the correlation output after the first path includes both the noise distribution and the signal plus noise distribution with some probabilities as expressed by equation (13).
Hence, the test can fail even though the block is an all-signal block. Therefore, some additional correlation outputs before {circumflex over (k)} can be employed as well when calculating the average power before the TOA estimate. In other words, if




  
    
      1
      ⁢
      
        
      
      ⁢
      
        1
        
          
            k
            ^
          
          -
          
            n
            s
          
          +
          
            M
            3
          
        
      
      ⁢
      
        
          ∑
          
            i
            =
            
              
                n
                s
              
              -
              
                M
                3
              
            
          
          
            
              k
              ^
            
            -
            1
          
        
        ⁢
        
          z
          i
          2
        
      
    
    >
    
      δ
      2
    
  
  ,




the block is considered as an all-signal block, where M3≧0 additional outputs are used depending on {circumflex over (k)}.

When it is determined that the block includes all signal outputs, the TOA is expected to be in one of the previous blocks. Therefore, the uncertainty region, i.e., the observation block, is shifted backwards, and the change detection method is repeated.
Although the invention has been described by way of examples of preferred embodiments, it is to be understood that various other adaptations and modifications may be made within the spirit and scope of the invention. Therefore, it is the object of the appended claims to cover all such variations and modifications as come within the true spirit and scope of the invention.

## Claims

1. A method for estimating a time of arrival (TOA) of an ultra wide band (UWB) signal received in a wireless communication system, comprising:
measuring an energy in a frame of the UWB received signal to determine a block in the frame, the block representing a coarse time of arrival of the UWB received signal; and
combining a plurality of time-delayed versions of a template signal with the block of the UWB received signal to select a particular template signal, the particular template signal identifying a particular chip representing a fine time of arrival of the UWB received signal, and in which the template signal is a transmitted signal corresponding to the UWB received signal, in which the fine TOA of the UWB received signal, τTOA, is represented by is
τTOA=kTc=kbTb+kcTc, and where kε[0, Nc−1] is the TOA in terms of the pulse duration Tc, Nc is the number of chips per frame interval, Tb is the block interval including B pulses (Tb=BTc), and kbε[0, NcB−1], and kcε[0,B−1] are integers that determine, respectively, in which block and pulse the UWB signal is received via a first signal path.

2. The method of claim 1, in which the time delays correspond to an interval of one-chip.

3. The method of claim 1, in which the combining produces a set of combined values, the set having one combined value for each combination of the template signal with the UWB received signal.

4. The method of claim 3, further comprising: selecting the particular template signal associated with the combined value that matches a correlation value to identify the chip that represents the fine time of arrival of the UWB received signal.

5. A method for estimating a time of arrival (TOA) of an ultra wide band (UWB) signal received in a wireless communication system, comprising:
receiving the ultra wide band (UWB) signal, in which the UWB signal includes a plurality of frames, and each frame includes a plurality of blocks, and each block includes a plurality of UWB pulses;
measuring an energy of the pulses in each block to determine an energy of each block in the frame;
selecting the block with a highest energy as a coarse time of arrival (TOA) of the received UWB signal;
combining a plurality of time-delayed versions of a template signal with the selected block of the received (UWB) signal to select a particular template signal, the particular template signal identifying a particular pulse representing a fine TOA of the received UWB signal, in which the UWB signal is represented by:





  
    
      
        s
        tx
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        E
      
      ⁢
      
        
          ∑
          
            j
            =
            
              -
              ∞
            
          
          ∞
        
        ⁢
        
          
            a
            j
          
          ⁢
          
            b
            
              ⌊
              
                j
                ⁢
                
                  /
                
                ⁢
                
                  N
                  f
                
              
              ⌋
            
          
          ⁢
          
            
              w
              tx
            
            ⁡
            
              (
              
                t
                -
                
                  jT
                  f
                
                -
                
                  
                    c
                    j
                  
                  ⁢
                  
                    T
                    c
                  
                
              
              )
            
          
        
      
    
  
  ,




where wtx(t) is a transmitted pulse having a duration Tc, E is a transmitted pulse energy, Tfis a frame time, Nf is a number of pulses representing one information symbol, Ts=NfTf is a symbol interval, and b└j/Nf┘ε{+1, −1} is a binary information symbol, and the received (UWB) signal is




  
    
      r
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        
          ∑
          
            l
            =
            1
          
          L
        
        ⁢
        
          
            E
          
          ⁢
          
            α
            l
          
          ⁢
          
            
              s
              rx
            
            ⁡
            
              (
              
                t
                -
                
                  
                    (
                    
                      l
                      -
                      1
                    
                    )
                  
                  ⁢
                  
                    T
                    c
                  
                
                -
                
                  τ
                  TOA
                
              
              )
            
          
        
      
      +
      
        n
        ⁡
        
          (
          t
          )
        
      
    
  
  ,




where αl is a channel coefficient for an lth path, L is a number of multipath components, τTOA is the fine TOA of the received (UWB) signal, and n(t) is noise, and where




  
    
      
        
          s
          rx
        
        ⁡
        
          (
          t
          )
        
      
      ⁢
      
        
      
      ⁢
      is
      ⁢
      
        
      
      ⁢
      
        
          s
          rx
        
        ⁡
        
          (
          t
          )
        
      
    
    =
    
      
        ∑
        
          j
          =
          
            -
            ∞
          
        
        ∞
      
      ⁢
      
        
          a
          j
        
        ⁢
        
          b
          
            ⌊
            
              j
              ⁢
              
                /
              
              ⁢
              
                N
                f
              
            
            ⌋
          
        
        ⁢
        
          
            w
            rx
          
          ⁡
          
            (
            
              t
              -
              
                jT
                f
              
              -
              
                
                  c
                  j
                
                ⁢
                
                  T
                  c
                
              
            
            )
          
        
      
    
  
  ,



with wrx(t) denoting a received UWB pulse with unit energy and a training sequence is bj=1∀j.

6. The method of claim 5, in which the time delays correspond to an interval of one-chip.

7. The method of claim 5, in which the fine TOA is an integer multiple of the pulse duration Tc.

8. The method of claim 7, further comprising:
obtaining a sub-pulse TOA resolution by employing a delay-lock-loop (DLL) after estimating the fine TOA.