# Patent

## Title
Time of arrival based positioning system

## Abstract

A TOA positioning system can be implemented that employs a calculated initial location of a wireless network device. For each of a plurality of reference wireless network devices, a distance between the wireless network device and the reference wireless network device is determined based, at least in part, on a round trip transit time between the wireless network device and the reference wireless network device. An initial location of the wireless network device can be calculated based, at least in part, on a location of each of the plurality of reference wireless network devices. A location of the wireless network device can be estimated based, at least in part, on the calculated initial location, the distance to each of the reference wireless network devices, and an initial distance calibration constant.

## Background
BACKGROUND
Embodiments of the inventive subject matter generally relate to the field of wireless communication and, more particularly, to a time of arrival (TOA) based positioning system.
A wireless communication device can use various position estimation techniques to determine an unknown location of the wireless communication device based on communicating with a plurality of reference wireless communication devices with known locations. For example, the wireless communication device can employ round trip time (RTT) based positioning techniques by determining the travel time of radio signals from the wireless communication device to the reference wireless communication devices. The wireless communication device can determine the distance to the reference wireless communication devices based on the determined travel time of the radio signals and can use time of arrival (TOA) positioning techniques to determine the unknown location.
SUMMARY
Various embodiments of a TOA based positioning mechanism are disclosed. In one embodiment, a wireless network device of a communication network determines a round trip transit time between the wireless network device and each of a plurality of reference wireless network devices of the communication network. For each of the plurality of reference wireless network devices, a distance between the wireless network device and the reference wireless network device is determined based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device. An initial location of the wireless network device is calculated based, at least in part, on a location of each of the plurality of reference wireless network devices. TOA computations are performed for estimating a location of the wireless network device based, at least in part, on the initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, and an initial distance calibration constant.



BRIEF DESCRIPTION OF THE DRAWINGS
The present embodiments may be better understood, and numerous objects, features, and advantages made apparent to those skilled in the art by referencing the accompanying drawings.
FIG. 1 is an example block diagram illustrating a mechanism for estimating the unknown location of a network device in a wireless communication network;
FIG. 2 is an example block diagram of one embodiment of the location calculation unit including a TOA positioning mechanism;
FIG. 3 is a flow diagram illustrating example operations of a TOA positioning mechanism;
FIG. 4 is an example block diagram of one embodiment of the location calculation unit including a hybrid TOA positioning mechanism;
FIG. 5 depicts a flow diagram illustrating example operations of a hybrid TOA positioning mechanism;
FIG. 6 is a continuation of FIG. 5 and also illustrates example operations of a hybrid TOA positioning mechanism; and
FIG. 7 is a block diagram of one embodiment of an electronic device including a mechanism for determining the location of the electronic device in a wireless communication network.



DESCRIPTION OF EMBODIMENT(S)
The description that follows includes exemplary systems, methods, techniques, instruction sequences, and computer program products that embody techniques of the present inventive subject matter. However, it is understood that the described embodiments may be practiced without these specific details. For instance, although examples refer to wireless local area network (WLAN) devices (e.g., 802.11n compatible devices) executing the location estimation techniques described herein, embodiments are not so limited. In other embodiments, various other devices and standards (e.g., WiMAX) can execute the location estimation techniques. In other instances, well-known instruction instances, protocols, structures, and techniques have not been shown in detail in order not to obfuscate the description.
A network device associated with an unknown location (“unknown network device”) can use a TOA based positioning algorithm to determine its location based on a plurality of reference network devices. Depending on the implementation, the TOA positioning algorithm may or may not require an accurately known distance calibration constant to estimate the location of the unknown network device. Conventional TOA positioning algorithms typically use, as inputs, an arbitrarily selected initial location of the unknown network device, an arbitrarily selected distance calibration constant, and the distance to and locations of each reference network device to iteratively calculate the location of the unknown network device. Typically, the initial location of the unknown network device is zero (e.g., where the X, Y, and Z coordinates are zero) or a randomly selected value. Likewise, the initial distance calibration constant is also typically zero or a randomly selected value. However, the conventional TOA positioning algorithms can be very sensitive to the initial location of the unknown network device. Furthermore, some conventional TOA positioning algorithms may be sensitive to the distance calibration constant. Employing a randomly selected value or a zero value of the initial location can render the conventional TOA positioning algorithms unsolvable (e.g., by making the value of one or more determinants approach zero or infinity) and can result in the conventional TOA positioning algorithms not converging to an appropriate solution (i.e., the location of the unknown network device). Also, not knowing or randomly selecting the distance calibration constant can result in poor performance of the conventional TOA positioning algorithms, can result in the conventional TOA positioning algorithms not converging to the appropriate location of the unknown network device, and/or can result in the conventional TOA positioning algorithms converging to the wrong solution.
In some implementations, a first TOA positioning system that employs a calculated (rather than an arbitrarily selected) initial location of the unknown network device can be implemented to minimize the sensitivity of some of the conventional TOA positioning algorithms to the initial location of the unknown network device. The first TOA positioning system is typically not sensitive to the distance calibration constant and may utilize an arbitrarily selected distance calibration constant for positioning calculations. The first TOA positioning system that can utilize the arbitrarily selected distance calibration constant is referred to herein as “TOA1 positioning system.” The initial location of the unknown network device can be calculated based, at least in part, on known locations of the reference network devices. The TOA1 positioning system can execute a single iteration (rather than multiple iterations) to calculate the location of the unknown network device and the distance calibration constant based, at least in part, on the calculated initial location of the unknown network device and the distance to each of the reference network devices. Calculating (rather than arbitrarily guessing) the initial location of the unknown network device can mitigate location sensitivity issues, ensure convergence of the TOA1 positioning system in determining the location of the unknown network device, and improve the performance and the positioning accuracy of the TOA1 positioning system. Moreover, by executing only a single iteration to estimate the location of the unknown network device, the TOA1 positioning system can minimize the resources and the time consumed to calculate the location of the unknown network device.
In some implementations, a hybrid TOA positioning system can improve the performance of estimating the location of the unknown network device. As part of the hybrid TOA positioning system, the TOA1 positioning system described above can execute a single iteration to estimate the location of the unknown network device (referred to herein as “intermediate location of the unknown network device”) and a distance calibration constant (referred to herein as “intermediate distance calibration constant”). The intermediate location of the unknown network device and the intermediate distance calibration constant (along with the location of the reference network devices and the distance to each of the reference network devices) can then be used to calculate a target (or estimated) distance calibration constant and a second intermediate location of the unknown network device. In addition, the hybrid TOA positioning system can also comprise a second TOA positioning system referred to herein as “TOA2 positioning system.” The second TOA positioning system is typically sensitive to the distance calibration constant and may not utilize an arbitrarily selected distance calibration constant for positioning calculations. Instead, the TOA2 positioning system can utilize the intermediate distance calibration constant calculated by the TOA1 positioning system for positioning calculations, as will be described below. The hybrid TOA positioning system can iteratively execute the TOA2 positioning system based, at least in part, on the target distance calibration constant and the second intermediate location to estimate the location (“estimated location”) of the unknown network device. Calculating (rather than guessing) the target distance calibration constant can minimize sensitivity to the distance calibration constant, minimize positioning error associated with determining the location of the unknown network device, increase positioning accuracy and improve performance gain, consequently improving the overall performance of the unknown network device.
FIG. 1 is an example block diagram illustrating a mechanism for estimating the unknown location of a network device in a wireless communication network 100. In one example, the wireless communication network 100 comprises a WLAN device 102 with an unknown location (“unknown WLAN device”) and four reference WLAN devices 104, 106, 108, and 112. The unknown WLAN device 102 comprises a location calculation unit 110. In some implementations, the location calculation unit 110 may be implemented in a communication unit of the unknown WLAN device 102 that implements protocols and functionality to enable WLAN communication in the wireless communication network 100. It is noted that although not depicted in FIG. 1, in some embodiments one or more of the reference WLAN devices 104, 106, 108, and 112 can also comprise a location calculation unit and corresponding functionality for determining their respective location. In some implementations, the unknown WLAN device 102 and the reference WLAN devices 104, 106, 108, and 112 can each be electronic devices with WLAN communication capabilities, such as a laptop computer, a tablet computer, a mobile phone, a smart appliance, a gaming console, an access point, or other suitable electronic devices. The unknown WLAN device 102 can execute operations to determine its location, as will be described below in stages A-F.
At stage A, the location calculation unit 110 determines the round trip transit time (RTT) between the unknown WLAN device 102 and each of the reference WLAN devices 104, 106, 108, and 112 in the wireless communication network 100. In one implementation, the location calculation unit 110 can transmit one or more control messages to the reference WLAN device 104 and can receive corresponding one or more response control messages (e.g., acknowledgement (ACK) messages) from the reference WLAN device 104. The location calculation unit 110 can also record the time instants at which the control messages were transmitted to and the corresponding response control messages were received from the reference WLAN device 104. The location calculation unit 110 can then calculate the RTT between the unknown WLAN device 102 and the reference WLAN device 104 (“RTT associated with the reference WLAN device”) as the elapsed time between transmitting the control messages and receiving the corresponding response control messages, as will further be described in blocks 304-308 of FIG. 3. Likewise, the location calculation unit 110 can also calculate the RTT associated with the reference WLAN device 106 and the RTT associated with the reference WLAN device 108. It is noted that in other implementations, other suitable techniques can be employed to determine the RTT associated with the reference WLAN devices 104, 106, 108, and 112.
At stage B, the location calculation unit 110 calculates the distance between the unknown WLAN device 102 and each of the reference WLAN devices 104, 106, 108, and 112 based on the corresponding RTT associated with the reference WLAN devices. Operations for calculating the distance between the unknown WLAN device 102 and each of the reference WLAN devices 104, 106, 108, and 112 will further be described in FIGS. 2 and 3.
At stage C, the location calculation unit 110 determines the location of each of the reference WLAN devices 104, 106, 108, and 112. In some implementations, the location calculation unit 110 can request and receive location coordinates associated with each of the reference WLAN devices 104, 106, 108, and 112. For example, the location calculation unit 110 can transmit a request for the location coordinates in the control messages transmitted at stage A and can receive the location coordinates in the response control messages. As another example, the location calculation unit 110 can transmit a location request message (distinct from the control messages transmitted at stage A) to each of the reference WLAN devices and can receive a corresponding location response message comprising the location coordinates associated with each of the reference WLAN devices. In another implementation, the location calculation unit 110 can query a centralized server (or can access a predetermined memory location) to determine the location coordinates associated with the reference WLAN devices 104, 106, 108, and 112. In another implementation, the reference WLAN devices 104, 106, 108, and 112 may broadcast their respective location coordinates at periodic intervals (e.g., in a beacon message or another suitable control message). The location calculation unit 110 can determine the location coordinates associated with the reference WLAN devices 104, 106, 108, and 112 based on receiving and analyzing the periodically received messages. It is noted that the location calculation unit 110 can determine two-dimensional (2-D) coordinates (e.g., X and Y coordinates), three-dimensional (3-D) coordinates (e.g., X, Y, and Z coordinates), latitudes and longitudes, spherical coordinates, and/or other suitable location indicators associated with the reference WLAN devices 104, 106, 108, and 112.
At stage D, the location calculation unit 110 calculates an initial location of the unknown WLAN device 102 based, at least in part, on the location of each of the reference WLAN devices 104, 106, 108, and 112. As one example, the location calculation unit 110 can calculate the initial location of the unknown WLAN device 102 as the average of the known locations of the reference WLAN devices 104, 106, 108, and 112, as will be described in FIGS. 2 and 3. The location calculation unit 110 can use the calculated initial location of the unknown WLAN device 102 as an input to a first time of arrival (TOA1) positioning algorithm to estimate the location of the unknown WLAN device 102, as will be described below in stage F and in FIGS. 2-3. In some implementations (e.g., in the hybrid TOA positioning system), the location calculation unit 110 can use the previously calculated location of the unknown WLAN device output (referred to herein as the “first intermediate location of the unknown WLAN device”) generated at the output of the TOA1 positioning algorithm to calculate a second intermediate location of the unknown WLAN device 102 based, at least in part, on the known locations of the reference WLAN devices 104, 106, 108, and 112, and the RTT associated with the reference WLAN devices 104, 106, 108, and 112, as will be described in FIGS. 4-6. In this implementation, the location calculation unit 110 can use the second intermediate location of the unknown WLAN device 102 as an input to a second time of arrival (TOA2) positioning algorithm to determine the estimated location of the unknown WLAN device 102, as will be described below in stage F and in FIGS. 4-7.
At stage E, the location calculation unit 110 determines a distance calibration constant. In some implementations, as will be described in accordance with FIG. 2, the initial distance calibration constant can be a predetermined value or a randomly selected value. The location calculation unit 110 can use the initial distance calibration constant as an input to the TOA1 positioning algorithm to estimate the location of the unknown WLAN device 102 and a distance calibration constant, as will be described below in stage F and in FIGS. 2-3. In other implementations, the location calculation unit 110 can use the distance calibration constant at the output of the TOA1 positioning algorithm (referred to herein as “intermediate distance calibration constant) to calculate a target (or estimated) distance calibration constant. In this implementation, the location calculation unit 110 can use the target distance calibration constant as an input to the TOA2 positioning algorithm to estimate the location of the unknown WLAN device 102, as will be described below in stage F and in FIGS. 4-7.
At stage F, the location calculation unit 110 estimates the location of the unknown WLAN device 102 based, at least in part, on the initial location of the unknown WLAN device 102, the distance calibration constant, and the distance to each of the reference WLAN devices 104, 106, 108, and 112. In some implementations, as will be described in FIGS. 2-3, calculating the location of the unknown WLAN device 102 can be a one-stage process. In this implementation, the location of the unknown WLAN device 102 can be calculated by executing a single iteration of a Taylor series-based TOA1 positioning algorithm using at least the calculated initial location of the unknown WLAN device 102, and the known locations of the reference WLAN devices as inputs. In some implementations, as will be described in FIGS. 5-7, calculating the location of the unknown WLAN device 102 can be a two-stage process. In this implementation, the location calculation unit 110 can first execute the Taylor series based TOA1 positioning algorithm (as described above) and determine a first intermediate location of the unknown WLAN device 102. Furthermore, the TOA1 positioning algorithm can also be used to determine an intermediate distance calibration constant. Next, the location calculation unit 110 can determine the second intermediate location of the unknown WLAN device 102 (as described in stage D) and the target distance calibration constant (as described above in stage E) based on the first intermediate location and the intermediate distance calibration constant. The location calculation unit 110 can then iteratively execute a Taylor series-based TOA2 positioning algorithm to estimate the location of the unknown WLAN device 102.
FIG. 2 is an example block diagram of one embodiment of the location calculation unit 110 including a TOA positioning mechanism. The location calculation unit 110 comprises a distance calculation unit 202, a TOA1 initial location calculation unit 204, and a TOA1 positioning unit 206. The distance calculation unit 202 and the TOA1 initial location calculation unit 204 are coupled with the TOA1 positioning unit 206.
The distance calculation unit 202 can calculate the distance between the unknown WLAN device 102 and each of the reference WLAN devices 104, 106, 108, and 112 based, at least in part, on the round-trip transit time (RTT) associated with each of the reference WLAN devices. As depicted in FIG. 2, the distance calculation unit 202 takes as an input parameter, the RTT 208 between the unknown WLAN device 102 and each of the reference WLAN devices 104, 106, 108, and 112. If i represents a counter for the each of the reference WLAN devices and N represents the number of reference WLAN devices in the communication network 100, then the input to the distance calculation unit 202 is RTTi (i.e., RTT1, RTT2, . . . RTTN) as depicted in FIG. 2. The output of the distance calculation unit 202 is the estimated distance 210 ({circumflex over (d)}i) between the unknown WLAN device 102 and the ith reference WLAN device. Thus, if RTTi represents the RTT associated with the ith reference WLAN device and c represents the speed of light, the distance {circumflex over (d)}i 210 can be calculated in accordance with Eq. 1a. It is noted that the actual distance (di) between the unknown WLAN device 102 and the ith reference WLAN device can be represented as the sum of a TOA1 distance calibration constant (dc—TOA1) and the estimated distance ({circumflex over (d)}i) between the unknown WLAN device 102 station and the ith reference WLAN device, as depicted in Eq. 1b.




  
    
      
        
          
            d
            ^
          
          i
        
        =
        
          c
          ×
          
            
              RTT
              i
            
            2
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          1
        
        ⁢
        
          
        
        ⁢
        a
      
    
  
  
    
      
        
          d
          i
        
        =
        
          
            
              d
              ^
            
            i
          
          +
          
            d
            
              c_TOA
              ⁢
              
                
              
              ⁢
              1
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          1
        
        ⁢
        
          
        
        ⁢
        b
      
    
  




The distance calibration constant (dc—TOA1) can represent the difference between the estimated distance ({circumflex over (d)}i) and the actual distance (di) and may be indicative of the internal processing time (or turnaround time) associated with the reference WLAN devices. For example, the distance calibration constant can account for the elapsed time between the reference WLAN device detecting a control message from the unknown WLAN device 102 and the reference WLAN device transmitting a response control message to the unknown WLAN device 102. The distance calibration constant may also depend on the type and configuration of the reference WLAN devices. In some implementations, the distance calibration constant may differ from one reference WLAN device to another, while in other implementations the distance calibration constant may not differ across the reference WLAN devices.
The TOA1 initial location calculation unit 204 can calculate the initial location of the unknown WLAN device 102 based, at least in part, on the location 212 of each of the reference WLAN devices 104, 106, 108, and 112. As depicted in FIG. 2, the known locations 212 of the reference WLAN devices 104, 106, 108, and 112 are provided as input parameters to the TOA1 initial location calculation unit 204. As described herein, i represents a counter for the reference WLAN devices and Xi represents the location 212 (e.g., also referred to as absolute location, actual location, or real location) of the ith reference WLAN device. In one implementation, the location of the ith reference WLAN device can be represented in the form, Xi={xis,s=0, 1, . . . , D−1}, where D represents the dimension of the coordinate system used to specify the location of the reference WLAN device (e.g., if the location of the reference WLAN devices is specified in terms of X, Y, and Z coordinates, then D=3), s represents a counter for the dimension of the coordinate system (i.e., the dimension index), and xis represents the value of the sth dimension of the location of the ith reference WLAN device. It is noted that the location of the reference WLAN devices can be expressed using other suitable coordinate system (e.g., Cartesian coordinates, spherical coordinates, geodesic coordinates, etc.). Furthermore, the location of the reference WLAN devices can be expressed in other suitable number of dimensions (e.g., 2-dimensional coordinates, 3-dimensional coordinates, etc.) As depicted in FIG. 2, the output of the TOA1 initial location calculation unit 204 is the initial location (YTOA1init) 214 of the unknown WLAN device 102. In one implementation, the initial location of the unknown WLAN device 102 can be represented as YTOA1init={yTOA1inits,s=0, 1, . . . , D−1}, where yTOA1inits represents the value of the sth dimension of the initial location of the unknown WLAN device 102. In some implementations, the initial location 214 of the unknown WLAN device 102 can be expressed as a function ƒ(.) of the known locations 212 of the N reference WLAN devices, as depicted in Eq. 2a. In some implementations, the initial location 214 of the unknown WLAN device 102 can be calculated as an average of the known locations 212 of the reference WLAN devices, as depicted in Eq. 2b. In another implementation, the initial location 214 of the unknown WLAN device 102 can be calculated by maximizing the determinant of a product of matrices generated during the execution of a Taylor series based TOA1 positioning algorithm, as depicted in Eq. 2c. Maximizing the expression det(BToA1TBTOA1) as depicted in Eq. 2c can ensure that the TOA1 positioning algorithm (described below) converges to the location of the unknown WLAN device 102. In Eq. 2c, BTOA1T represents the transpose of the BTOA1 matrix and the BTOA1 matrix will further be described below with reference to Exp. 11 of FIG. 2 and with reference to FIG. 3.




  
    
      
        
          
        
        ⁢
        
          
            Y
            
              TOA
              ⁢
              
                
              
              ⁢
              
                1
                init
              
            
          
          =
          
            f
            ⁡
            
              (
              
                
                  X
                  0
                
                ,
                
                  X
                  1
                
                ,
                
                  ...
                  ⁢
                  
                    
                  
                  ⁢
                  
                    X
                    
                      N
                      -
                      1
                    
                  
                
              
              )
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          2
        
        ⁢
        
          
        
        ⁢
        a
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            Y
            
              TOA
              ⁢
              
                
              
              ⁢
              
                1
                init
              
            
          
          =
          
            
              f
              ⁡
              
                (
                
                  
                    X
                    0
                  
                  ,
                  
                    X
                    1
                  
                  ,
                  
                    ...
                    ⁢
                    
                      
                    
                    ⁢
                    
                      X
                      
                        N
                        -
                        1
                      
                    
                  
                
                )
              
            
            =
            
              
                1
                N
              
              ⁢
              
                
                  ∑
                  
                    i
                    =
                    0
                  
                  
                    N
                    -
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                
                  X
                  i
                
              
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          2
        
        ⁢
        
          
        
        ⁢
        b
      
    
  
  
    
      
        
          Y
          
            TOA
            ⁢
            
              
            
            ⁢
            
              1
              init
            
          
        
        =
        
          
            f
            ⁡
            
              (
              
                
                  X
                  0
                
                ,
                
                  X
                  1
                
                ,
                
                  ...
                  ⁢
                  
                    
                  
                  ⁢
                  
                    X
                    
                      N
                      -
                      1
                    
                  
                
              
              )
            
          
          =
          
            
              max
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  
                    1
                    init
                  
                
              
            
            ⁢
            
              {
              
                det
                ⁡
                
                  (
                  
                    
                      B
                      
                        TOA
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                      T
                    
                    ⁢
                    
                      B
                      
                        TOA
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                    
                  
                  )
                
              
              }
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          2
        
        ⁢
        
          
        
        ⁢
        c
      
    
  




It is noted that in other implementations, the TOA1 initial location calculation unit 204 can employ other suitable techniques to calculate the initial location 214 of the unknown WLAN device 102. For example, the initial location 214 of the unknown WLAN device 102 can be calculated as a weighted average of the known locations 212 of the reference WLAN devices, where the weights are selected based on the confidence of the location of the reference WLAN devices and/or based the distance between the reference WLAN devices and the unknown WLAN device 102. As another example, the initial location 214 of the unknown WLAN device 102 can be calculated as a weighted (or non-weighted) combination of the known locations 212 of only a subset of the reference WLAN devices.
The TOA1 positioning unit 206 can execute a single iteration of a Taylor series based TOA1 positioning algorithm to calculate the location 220 of the unknown WLAN device 102 and to estimate a distance calibration constant 218 based, at least in part, on the initial location 214 of the unknown WLAN device 102, the estimated distances 210 to each of the reference WLAN devices, and the location 212 of each of the reference WLAN devices. As depicted in FIG. 2, the estimated distances 210 between the unknown WLAN device 102 and each of the reference WLAN devices (calculated by the distance calculation unit 202) are provided as one set of inputs to the TOA1 positioning unit 206. The initial location 214 of the unknown WLAN device 102 (YTOA1init) calculated by the TOA1 initial location calculation unit 204 is provided as another input to the TOA1 positioning unit 206. The known locations 212 (Xi) of the reference WLAN devices 104, 106, 108, and 112 are also provided as an input to the TOA1 positioning unit 206. Additionally, an initial distance calibration constant (dcinit) 216 can also be provided as an input to the TOA1 positioning unit 206. As described above, the initial distance calibration constant 216 may be indicative of internal processing time (or turnaround time) associated with the reference WLAN devices. In some implementations, the initial distance calibration constant 216 may be selected to be zero, while in other implementations the initial distance calibration constant 216 may be selected to be another suitable random value. In some implementations, the initial distance calibration constant 216 may be preconfigured, while in other implementations the initial distance calibration constant 216 may be dynamically selected (e.g., a random value may be dynamically selected for the initial distance calibration constant 216).
The TOA1 positioning unit 206 can then execute one iteration of a Taylor series based TOA1 positioning algorithm to solve the expression depicted by Eq. 3a and to estimate the location (YTOA1) 220 of the unknown WLAN device 102. In accordance with Eq. 3a, the TOA1 positioning unit 206 can attempt to identify the value of YTOA1 that minimizes the sum (across all the reference WLAN devices) of the estimation error between A) the measured distance to the ith reference WLAN device determined in Eq. 1 based on the RTT associated with the ith reference WLAN device and B) the estimated distance to the ith reference WLAN device based on the known location of the ith reference WLAN device. In Eq. 3a, ∥Xi−YTOA1∥ represents the estimated distance between the ith reference WLAN device and the unknown WLAN device 102. ∥Xi−YTOA1∥ can be calculated in accordance with Eq. 3b. The estimated distance between a reference WLAN device and the unknown WLAN device can be determined based on the actual location of the reference WLAN device and the estimated location (YTOA1) of the unknown WLAN device.




  
    
      
        
          Y
          
            TOA
            ⁢
            
              
            
            ⁢
            1
          
        
        =
        
          arg
          ⁢
          
            
              min
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
            
            ⁢
            
              {
              
                
                  ∑
                  
                    i
                    =
                    0
                  
                  
                    N
                    -
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                
                  (
                  
                    
                      
                        d
                        ^
                      
                      i
                    
                    +
                    
                      d
                      
                        c_TOA
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                    
                    -
                    
                      
                      
                        
                          X
                          i
                        
                        -
                        
                          Y
                          
                            TOA
                            ⁢
                            
                              
                            
                            ⁢
                            1
                          
                        
                      
                      
                    
                  
                  }
                
              
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          3
        
        ⁢
        
          
        
        ⁢
        a
      
    
  
  
    
      
        
          
          
            
              X
              i
            
            -
            
              Y
              
                TOA
                ⁢
                
                  
                
                ⁢
                1
              
            
          
          
        
        =
        
          
            [
            
              
                ∑
                
                  s
                  =
                  0
                
                
                  D
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                
                  (
                  
                    
                      x
                      i
                      s
                    
                    -
                    
                      y
                      
                        TOA
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                      s
                    
                  
                  )
                
                2
              
            
            ]
          
          
            1
            2
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          3
        
        ⁢
        
          
        
        ⁢
        b
      
    
  




In some implementations, by assuming that there is no estimation error {e.g., by assuming that ({circumflex over (d)}i+dc—TOA1−∥Xi−YTOA1∥)=0}, the TOA1 positioning unit 206 can simplify Eq. 3a and can instead solve Eq. 4 for all values of i to estimate the location of the unknown WLAN device 102.

{circumflex over (d)}i=∥Xi−YTOA1∥−dc—TOA1  Eq. 4

The TOA1 distance calibration constant (dc—TOA1) 218 can be expressed as the sum of the initial distance calibration constant (dcinit) 216 and a distance calibration error (Δdc), as depicted in Eq. 5a. In other words, the distance calibration error (Δdc) can be expressed as the difference (or error) between the estimated distance calibration constant (dc—TOA1) 218 and the initial distance calibration constant (dcinit) 216. Likewise, the location 220 of the unknown WLAN device 102 (YTOA1) can be expressed as the sum of the initial location 214 of the unknown WLAN device (YTOA1init) and a TOA1 location error (ΔYTOA1), as depicted in Eq. 5b. In other words, the TOA1 location error (ΔYTOA1) can be expressed as the difference (or error) between the estimated location 220 of the unknown WLAN device 102 (YTOA1) and the initial location 214 of the unknown WLAN device (YTOA1init). It is noted that in some implementations, the TOA1 location error can be represented as ΔYTOA1={ΔyTOA1s,s=0, 1, . . . , D−1}, where D represents the dimension of the coordinate system and ΔyTOA1s represents the value of the sth dimension of the TOA1 location error. The TOA1 positioning unit 206 can substitute Eq. 5a and Eq. 5b into Eq. 4 to yield Eq. 6.

dc—TOA1=dcinit+Δdc  Eq. 5a

YTOA1=YTOA1init+ΔYTOA1  Eq. 5b

{circumflex over (d)}i=∥Xi−YTOA1init−ΔYTOA1∥−dcinit−Δdc  Eq. 6

The TOA1 positioning unit 206 can then use Taylor series expansion and neglect the second and higher order terms (e.g., because the higher order terms have a zero or negligible value), to yield Eq. 7. The TOA1 positioning unit 206 can generate an equation in accordance with Eq. 7 for each of the reference WLAN devices 104, 106, 108, and 112 and can organize the plurality of generated equations in matrix form to yield Eq. 8.




  
    
      
        
          
            d
            ^
          
          i
        
        =
        
          
            
            
              
                X
                i
              
              -
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  
                    1
                    init
                  
                
              
            
            
          
          -
          
            d
            
              c
              init
            
          
          +
          
            
              ∑
              
                s
                =
                0
              
              
                D
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  -
                  
                    
                      
                        x
                        i
                        s
                      
                      -
                      
                        y
                        
                          TOA
                          ⁢
                          
                            
                          
                          ⁢
                          
                            1
                            init
                          
                        
                        s
                      
                    
                    
                      
                      
                        
                          X
                          i
                        
                        -
                        
                          Y
                          
                            TOA
                            ⁢
                            
                              
                            
                            ⁢
                            
                              1
                              init
                            
                          
                        
                      
                      
                    
                  
                
                )
              
              ⁢
              Δ
              ⁢
              
                
              
              ⁢
              
                y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
                s
              
            
          
          -
          
            Δ
            ⁢
            
              
            
            ⁢
            
              d
              c
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        7
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            A
            
              ToA
              ⁢
              
                
              
              ⁢
              1
            
          
          =
          
            
              
                B
                
                  ToA
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              ⁡
              
                [
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    Y
                    
                      TOA
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    c
                  
                
                ]
              
            
            T
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        8
      
    
  




In some implementations, the TOA1 location error matrix (ΔYTOA1) may be a 1×D row matrix that can represent the difference between the estimated location of the unknown WLAN device 102 and the initial TOA1 location of the unknown WLAN device 102 as depicted in Exp 9. The TOA1 distance error matrix (AToA1) can be an N×1 column matrix, where each element (i.e., each row) of the AToA1 matrix is represented by Exp. 10 and each row is associated with one of the reference WLAN devices. As depicted by Exp. 10, each element of the TOA1 distance error matrix represents a combination of A) the estimated distance between the unknown WLAN device and a reference WLAN device (e.g., based on RTT associated with the reference WLAN device as depicted in Eq. 1), B) the initial distance calibration constant 216, and C) the estimated distance between the unknown WLAN device and the reference WLAN device based on the known location of the reference WLAN device and the initial location 214 of the unknown WLAN device 102. The TOA coefficient matrix (BToA1) can be an N×(D+1) matrix, where N represents the number of reference WLAN devices and D represents the dimensions of the coordinate system that is used to specify the location of the reference WLAN devices. Each element of the TOA1 coefficient matrix can be represented in accordance with Exp. 11. The TOA1 coefficient matrix (BTOA1) can comprise the coefficients of the TOA1 location error matrix (ΔYToA1) and can represent the relationship between the TOA1 location error matrix (ΔYToA1) and the TOA1 distance error matrix (AToA1).





  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            Y
            
              TOA
              ⁢
              
                
              
              ⁢
              1
            
          
        
        =
        
          [
          
            
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    y
                    
                      TOA
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                    0
                  
                
              
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    y
                    
                      TOA
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                    1
                  
                
              
              
                ...
              
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    y
                    
                      TOA
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                    
                      D
                      -
                      2
                    
                  
                
              
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    y
                    
                      TOA
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                    
                      D
                      -
                      1
                    
                  
                
              
            
          
          ]
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        9
      
    
  
  
    
      
        
          
            i
            th
          
          ⁢
          
            
          
          ⁢
          element
          ⁢
          
            
          
          ⁢
          of
          ⁢
          
            
          
          ⁢
          
            A
            
              TOA
              ⁢
              
                
              
              ⁢
              1
            
          
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            
              d
              ^
            
            i
          
        
        +
        
          d
          
            c
            init
          
        
        -
        
          
          
            
              X
              i
            
            -
            
              Y
              
                TOA
                ⁢
                
                  
                
                ⁢
                
                  1
                  init
                
              
            
          
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        10
      
    
  
  
    
      
        
          
            
              (
              
                i
                ,
                s
              
              )
            
            th
          
          ⁢
          
            
          
          ⁢
          element
          ⁢
          
            
          
          ⁢
          of
          ⁢
          
            
          
          ⁢
          
            B
            
              ToA
              ⁢
              
                
              
              ⁢
              1
            
          
          ⁢
          
            :
          
        
        ⁢
        
          
        
        -
        
          
            
              x
              i
              s
            
            -
            
              y
              
                TOA
                ⁢
                
                  
                
                ⁢
                
                  1
                  init
                
              
              s
            
          
          
            
            
              
                X
                i
              
              -
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  
                    1
                    init
                  
                
              
            
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        11
      
    
  




In one implementation, the dimension (e.g., the value of s) can vary from column to column while the reference device under consideration (e.g., the value of i) can vary from row to row. The dimension can remain constant throughout a particular column while the reference device under consideration can remain constant throughout a particular row. The TOA1 positioning unit 206 can then calculate the TOA1 location error ΔYTOA1 and the distance calibration error Δdc in accordance with Eq. 12. It is noted that in some implementations, prior to evaluating Eq. 12, the TOA1 positioning unit 206 may first determine whether the condition det(BTOA1TBTOA1)≠0 is satisfied to ensure convergence of the Taylor series based TOA1 positioning algorithm.

[ΔYTOAΔdc]T=(BToA1TBToA1)−1BToA1TAToA1  Eq. 12

After the TOA1 positioning unit 206 calculates the value of TOA1 location error and the distance calibration error (in accordance with Eq. 12), the TOA1 positioning unit 206 can use Eq. 5a to calculate the TOA1 distance calibration constant (dc—TOA1) 218 by summing the distance calibration error (Δdc) with the initial distance calibration constant (dcinit). Likewise, the TOA1 positioning unit 206 can use Eq. 5b to calculate the estimated location 220 of the unknown WLAN device (YTOA1) by summing the TOA1 location error (ΔYTOA1) with the initial location 214 of the unknown WLAN device 102 (YTOA1init). In one implementation, the location 220 of the unknown WLAN device 102 at the output of the TOA1 positioning unit 206 can be represented as YTOA1={yTOA1s,s=0, 1, . . . , D−1}, where yTOA1s represents the value of the sth dimension of the location of the unknown WLAN device 102. It is noted that the TOA1 positioning unit 206 may execute the operations described herein with reference to Eq. 3a-Eq. 12 only once (i.e., only one iteration) to calculate the estimated location of the unknown WLAN device 102.
FIG. 3 depicts a flow diagram (“flow”) 300 illustrating example operations of a TOA positioning mechanism. The flow 300 begins at block 302.
At block 302, an unknown network device determines a plurality of reference network devices of a communication network based on which to estimate the location of the unknown network device. With reference to the example of FIG. 1, the location calculation unit 110 of the unknown WLAN device 102 can identify the reference WLAN devices 104, 106, 108, and 112 based on which to estimate the location of the unknown WLAN device 102. In some implementations, the location calculation unit 110 can identify the reference WLAN devices 104, 106, 108, and 112 from a plurality of WLAN devices in the wireless communication network 100 based on analyzing one or more performance measurements associated with the plurality of WLAN devices. For example, location calculation unit 110 can identify N WLAN devices in the wireless communication network 100 with the highest received signal strength indicator (RSSI), the lowest packet error rate (PER), etc. In some implementations, if the wireless communication network 100 comprises a plurality of access points, the access points may advertise (e.g., in a beacon message) their presence/location to other WLAN devices in the wireless communication network 100. The location calculation unit 110 can detect these beacon messages, identify the access points that transmitted the beacon messages, and select at least a subset of the access points as reference WLAN devices (e.g., based on analyzing the performance measurements associated with the beacon messages received from each of the access points). The flow continues at block 304.
At block 304, a control message is transmitted to each of the plurality of reference network devices. For example, the location calculation unit 110 can transmit a control message to the reference WLAN device 104. The control message can comprise one or more of an identifier associated with the unknown WLAN device 102, an identifier associated with the reference WLAN device 104, and a suitable payload (e.g., a predetermined combination of symbols, a NULL payload, etc.). The location calculation unit 110 can also record the time instant at which the control message was transmitted to the reference WLAN device 104. Likewise, with reference to FIG. 1, the location calculation unit 110 can also record the time instants at which control messages were transmitted to the reference WLAN devices 106, 108, and 112. The flow continues at block 306.
At block 306, a response control message is received from each of the plurality of reference network devices. For example, the location calculation unit 110 can receive the response control message from the reference WLAN device 104. The response control message can be a WLAN acknowledgment (ACK) message or any suitable message that indicate receipt of the control message (transmitted at block 304) at the reference WLAN device 104. The location calculation unit 110 can also record the time instant at which the response control message was received at the unknown WLAN device 102 from the WLAN device 104. Likewise, with reference to FIG. 1, the location calculation unit 110 can also record the time instants at which response control messages were received from the reference WLAN devices 106, 108, and 112. The flow continues at block 308.
At block 308, the round-trip transit time (RTT) associated with each of the reference network devices is determined. In one implementation, the location calculation unit 110 can determine the RTT associated with the reference WLAN device 104 based on the time instants at which the control message was transmitted to the reference WLAN device 104 (recorded at block 304) and the time instant at which the response control message was received from the reference network device 104 (recorded at block 306). The location calculation unit 110 can compute the RTT associated with the reference WLAN device 104 by subtracting the time instant at which the control message was transmitted from the time instant at which the response control message was received. It is noted that in other implementations, the location calculation unit 110 can employ other suitable techniques to determine the RTT associated with the reference WLAN device 104. Likewise, with reference to FIG. 1, the location calculation unit 110 can also determine the RTT associated with the reference WLAN devices 106, 108, and 112. The flow continues at block 310.
At block 310, for each of the plurality of reference network devices, the distance between the unknown network device and the reference network device is calculated based, at least in part, on the RTT associated with the reference network device. For example, the location calculation unit 110 (e.g., the distance calculation unit 202 of FIG. 2) can calculate the distance between the unknown WLAN device 102 and each of the reference network devices 104, 106, 108, and 112 based, at least in part, on the RTT associated with the corresponding reference WLAN device as described above with reference to Eq. 1a of FIG. 2. It is noted that in other implementations, the distance calculation unit 202 can use other suitable techniques to determine the distance between the unknown WLAN device 102 and each of the reference WLAN devices 104, 106, 108, and 112. The flow continues at block 312.
At block 312, the initial location of the unknown network device is calculated based, at least in part, on known locations of the reference network devices. For example, location calculation unit 110 (e.g., the TOA1 initial location calculation unit 204 of FIG. 2) can calculate the initial location 214 of the unknown WLAN device 102. As described above with reference to Eq. 2a, the initial location 214 of the unknown WLAN device 102 can be calculated as a function of the location 212 of the reference WLAN devices. In some examples, as described above in Eq. 2b, the initial location of the unknown WLAN device 102 can be calculated as an average of the locations of the reference WLAN devices. In other examples, as described above in Eq. 2c, the initial location of the unknown WLAN device 102 can be calculated as another suitable weighted (or non-weighted) combination of the locations of the reference WLAN devices. More specifically, as described above in FIG. 2, the location calculation unit 110 can determine the TOA1 location error matrix (ΔYTOA1) such that each element of the TOA1 location error matrix represents the difference between a coordinate of the estimated location of the unknown WLAN device (e.g., the estimated X-coordinate) and an initial location coordinate of the initial location of the unknown WLAN device 102 (e.g., the X-coordinate of the initial location of the unknown WLAN device 102). The location calculation unit 110 can calculate the TOA1 distance error matrix (AToA1) in accordance with Exp. 10. The location calculation unit 110 can then calculate the TOA1 coefficient matrix (BTOA1) in accordance with Exp 11. The location calculation unit 110 can calculate the initial location of the wireless network device based, at least in part, on maximizing the determinant of the product of the BTOA1 matrix and a transpose of the BTOA1 matrix, as depicted in Eq. 2b. The flow continues at block 314.
At block 314, a single iteration of a first time of arrival (TOA1) positioning algorithm is executed to estimated location of the unknown network device based, at least in part, on the initial location of the unknown network device and the distance between the unknown network device and each of the reference network devices. For example, the location calculation unit 110 (e.g., the TOA1 positioning unit 206 of FIG. 2) can execute a single iteration of a Taylor series based TOA1 positioning algorithm to calculate the location of the unknown WLAN device 102. In addition, the TOA1 positioning unit 206 can also receive (as an input) an initial distance calibration constant 216. As described above in FIG. 2, the initial distance calibration constant 216 can be zero, a predetermined value, or a randomly selected value. As described above in Eq. 3a-Eq. 12, the TOA1 positioning unit 206 can use the initial location 214 of the unknown WLAN device 102 (determined at block 312), the distance 210 to each of reference WLAN devices (determined at block 310), the initial distance calibration constant 216, and the known locations 212 of the reference WLAN devices as inputs to execute a single iteration of the Taylor series based TOA1 positioning algorithm and to calculate the location of the unknown WLAN device 102. Therefore, the TOA1 positioning unit 206 can also determine a TOA1 distance calibration constant 218 based, at least in part, on the initial distance calibration constant 216 (described above) by executing the single iteration of the Taylor series based TOA1 positioning algorithm. It is noted that in other embodiments, the TOA1 positioning unit 206 can employ other suitable techniques to estimate the location of the unknown WLAN device 102 based, at least in part, on the calculated initial location 214 of the unknown WLAN device 102 (e.g., as will be further described below). From block 314, the flow ends.
FIG. 4 is an example block diagram of one embodiment of the location calculation unit 110 including a hybrid TOA positioning mechanism. The location calculation unit 110 comprises the distance calculation unit 202, the TOA1 initial location calculation unit 204, and the TOA1 positioning unit 206 previously described in FIGS. 2-3. Additionally, the location calculation unit 110 also comprises a TOA2 initial value calculation unit 402 and a TOA2 positioning unit 404. As shown in the example of FIG. 4, the distance calculation unit 202 and the TOA1 initial location calculation unit 204 are coupled with the TOA1 positioning unit 206. The TOA1 positioning unit 206 and the distance calculation unit 202 are coupled with the TOA2 initial value calculation unit 402. The distance calculation unit 202 and the TOA2 initial value calculation unit 402 are coupled with the TOA2 positioning unit 404.
The distance calculation unit 202 can calculate the distance 210 between the unknown WLAN device 102 and each of the reference WLAN devices 104, 106, 108, and 112 based, at least in part, on the corresponding RTT 208 associated with the reference WLAN devices, as described above in accordance with Eq. 1a of FIG. 2. The TOA1 initial location calculation unit 204 can calculate the initial location 214 of the unknown WLAN device 102 based, at least in part, on the location 212 of each of the reference WLAN devices 104, 106, 108, and 112, as described above in accordance with Eq. 2a, Eq. 2b, and Eq. 2c of FIG. 2. The TOA1 positioning unit 206 can employ a single iteration of the Taylor series based TOA1 positioning algorithm to estimate the TOA1 location 220 of the unknown WLAN device 102 and the distance calibration constant 218 based, at least in part, on the initial location 212 of the unknown WLAN device 102, the initial distance calibration constant 216, and the distances 210, as described above in accordance with Eq. 3a-Eq. 12 of FIG. 2. In the hybrid TOA mechanism, the distance calibration constant 218 and the location 220 at the output of the TOA1 positioning unit 206 can be referred to as “intermediate distance calibration constant” and “first intermediate location of the unknown WLAN device” respectively.
The TOA2 initial value calculation unit 402 can calculate a distance calibration constant (“target distance calibration constant”) and an initial location (“second intermediate location of the unknown WLAN device”) that will be provided as inputs to the TOA2 positioning unit 404. The TOA2 initial value calculation unit 402 can calculate the target distance calibration constant and the second intermediate location of the unknown WLAN device based, at least in part, on the estimated distances ({circumflex over (d)}i) 210 to each of the reference WLAN devices 104, 106, 108, and 112 (determined by the distance calculation unit 202), the first intermediate location (YTOA1) 220 of the unknown WLAN device 102 and the intermediate distance calibration constant 218 (calculated by the TOA1 positioning unit 206). Additionally, the locations (Xi) 210 of the reference WLAN devices 104, 106, 108, and 112 are provided as input parameters to the TOA2 initial value calculation unit 402. As depicted in FIG. 4, the TOA2 initial value calculation unit 402 can estimate the target distance calibration constant (dc) 406 and the second intermediate location (YTOA2init) 408 of the unknown WLAN device 102. In some implementations, the target distance calibration constant 406 can be represented as a function g(.) of the intermediate distance calibration constant 218, the first intermediate location 220 of the unknown WLAN device, the known locations 212 of the reference WLAN devices, and the distances 210 to the reference WLAN devices, as depicted in Eq. 13a. In one example, the TOA2 initial value calculation unit 402 can calculate the value of the target distance calibration constant 406 value as the weighted sum of the difference in A) the estimated distance (∥Xi−YTOA1∥) between the unknown WLAN device and the reference WLAN devices based on known locations of the reference WLAN devices and the first intermediate location of the unknown WLAN device, and B) the measured distance between the unknown WLAN device and the reference WLAN device based on the RTT associated with the reference WLAN devices, as depicted by Eq. 13b. As another example, the TOA2 initial value calculation unit 402 can assign the intermediate distance calibration constant 220 as the target distance calibration constant 406, as depicted in accordance with Eq. 13c.




  
    
      
        
          d
          c
        
        =
        
          g
          ⁡
          
            (
            
              
                d
                
                  c
                  
                    TOA
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
              
              ,
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              ,
              
                X
                0
              
              ,
              
                X
                1
              
              ,
              
                ...
                ⁢
                
                  
                
                ⁢
                
                  X
                  
                    N
                    -
                    1
                  
                
              
              ,
              
                
                  d
                  ^
                
                0
              
              ,
              
                
                  d
                  ^
                
                1
              
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  d
                  ^
                
                
                  N
                  -
                  1
                
              
            
            )
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          13
        
        ⁢
        
          
        
        ⁢
        a
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            d
            c
          
          =
          
            
              ∑
              
                i
                =
                0
              
              
                N
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                β
                i
              
              ⁡
              
                (
                
                  
                    
                    
                      
                        X
                        i
                      
                      -
                      
                        Y
                        
                          TOA
                          ⁢
                          
                            
                          
                          ⁢
                          1
                        
                      
                    
                    
                  
                  -
                  
                    
                      d
                      ^
                    
                    i
                  
                
                )
              
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          13
        
        ⁢
        
          
        
        ⁢
        b
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            d
            c
          
          =
          
            d
            
              c_TOA
              ⁢
              
                
              
              ⁢
              1
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          13
        
        ⁢
        
          
        
        ⁢
        c
      
    
  




In Eq. 13b, βi represents the weighting factor associated with each of the reference WLAN devices and can be selected so that each of the weighting factors are greater than or equal to zero and the sum of all the weighting factors equals one. In one example, the weighting factor βi associated with each reference WLAN device can be calculated as the inverse of the number of reference WLAN devices, as depicted by Eq. 13d. In another example, the weighting factor βi associated with each reference WLAN device can be calculated in accordance with Eq. 13e.




  
    
      
        
          
            β
            i
          
          =
          
            1
            N
          
        
        ,
        
          i
          =
          0
        
        ,
        1
        ,
        ...
        ⁢
        
          
        
        ,
        
          N
          -
          1
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          13
        
        ⁢
        
          
        
        ⁢
        d
      
    
  
  
    
      
        
          
            β
            i
          
          =
          
            1
            
              
                
                
                  
                    X
                    i
                  
                  -
                  
                    Y
                    
                      TOA
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                  
                
                
              
              ⁢
              
                
                  ∑
                  
                    j
                    =
                    0
                  
                  
                    N
                    -
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                
                  1
                  
                    
                    
                      
                        X
                        j
                      
                      -
                      
                        Y
                        
                          TOA
                          ⁢
                          
                            
                          
                          ⁢
                          1
                        
                      
                    
                    
                  
                
              
            
          
        
        ,
        
          

        
        ⁢
        
          i
          -
          0
        
        ,
        1
        ,
        ...
        ⁢
        
          
        
        ,
        
          N
          -
          1
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          13
        
        ⁢
        
          
        
        ⁢
        e
      
    
  




With reference to Eq. 13e, the weighting factor can be determined based on distance measurements to (e.g., the confidence in location of) the reference WLAN devices. In other words, the distance (∥Xi−YTOA1∥) between the unknown WLAN device 102 and the ith reference WLAN device can be used to calculate the weighting factor and can be used as an indication of the confidence in the measurements associated with the ith reference WLAN device. For example, if the reference WLAN device 104 is close to the unknown WLAN device 102, then ∥Xi−YTOA1∥ has a small value, its inverse has a large value, and consequently, the reference WLAN device 104 is associated with a higher weighting factor (e.g., higher importance). As another example, if the reference WLAN device 108 is far away from the unknown WLAN device 102, then ∥Xi−YTOA1∥ has a large value, its inverse has a small value, and consequently, the reference WLAN device 108 is associated with a small weighting factor (e.g., low importance). The factor




  
    ∑
    
      j
      =
      0
    
    
      N
      -
      1
    
  
  ⁢
  
    
  
  ⁢
  
    1
    
      
      
        
          X
          j
        
        -
        
          Y
          
            TOA
            ⁢
            
              
            
            ⁢
            1
          
        
      
      
    
  




is a constant multiplicative factor that is the sum of the inverse of the distance between the unknown WLAN device 102 and each of the reference WLAN devices. It is noted that in other implementations, the TOA2 initial value calculation unit 402 can use other suitable techniques to calculate the weighting factor βi associated with each of the reference WLAN devices.

Additionally, the TOA2 initial value calculation unit 402 also determines the second intermediate location (YTOA2init) 408 of the unknown WLAN device 102 as a function h(.) of the first intermediate location (YTOA1) 220 of the unknown WLAN device 102 (calculated by the TOA1 positioning unit 206) and the locations 210 of the reference WLAN devices, as depicted in Eq. 14a. In some implementations, the second intermediate location 408 of the unknown WLAN device 102 can be calculated as an average of the locations 210 of the reference WLAN devices, as depicted in Eq. 14b. In another implementation, the first intermediate location (YTOA1) 220 of the unknown WLAN device 102 can be designated as the second intermediate location 408 of the unknown WLAN device 102, as depicted in Eq. 14c.




  
    
      
        
          T
          
            TOA
            ⁢
            
              
            
            ⁢
            
              2
              init
            
          
        
        =
        
          h
          ⁡
          
            (
            
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              ,
              
                X
                0
              
              ,
              
                X
                1
              
              ,
              
                ...
                ⁢
                
                  
                
                ⁢
                
                  X
                  
                    N
                    -
                    1
                  
                
              
            
            )
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          14
        
        ⁢
        
          
        
        ⁢
        a
      
    
  
  
    
      
        
          T
          
            TOA
            ⁢
            
              
            
            ⁢
            
              2
              init
            
          
        
        =
        
          
            h
            ⁡
            
              (
              
                
                  Y
                  
                    TOA
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
                ,
                
                  X
                  0
                
                ,
                
                  X
                  1
                
                ,
                
                  ...
                  ⁢
                  
                    
                  
                  ⁢
                  
                    X
                    
                      N
                      -
                      1
                    
                  
                
              
              )
            
          
          =
          
            
              1
              N
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  0
                
                
                  N
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                X
                i
              
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          14
        
        ⁢
        
          
        
        ⁢
        b
      
    
  
  
    
      
        
          T
          
            TOA
            ⁢
            
              
            
            ⁢
            
              2
              init
            
          
        
        =
        
          
            h
            ⁡
            
              (
              
                
                  Y
                  
                    TOA
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
                ,
                
                  X
                  0
                
                ,
                
                  X
                  1
                
                ,
                
                  ...
                  ⁢
                  
                    
                  
                  ⁢
                  
                    X
                    
                      N
                      -
                      1
                    
                  
                
              
              )
            
          
          =
          
            Y
            
              TOA
              ⁢
              
                
              
              ⁢
              1
            
          
        
      
    
    
      
        
          Eq
          .
          
            
          
          ⁢
          14
        
        ⁢
        
          
        
        ⁢
        c
      
    
  




It is noted that in some implementations, the second intermediate location of the unknown WLAN device 102 can be represented as YTOA2init={yTOA2inits,s=0, 1, . . . , D−1}, where yTOA2inits represents the value of the sth dimension of the second intermediate location of the unknown WLAN device 102. It is also noted that in some implementations, the TOA2 initial value calculation unit 402 can employ other suitable techniques to calculate the second intermediate location 408 of the unknown WLAN device 102. For example, the second intermediate location 408 of the unknown WLAN device 102 can be calculated as a weighted average of the locations of the reference WLAN devices where the weights are selected based on the confidence of the location of the reference WLAN devices and/or based the distance between the reference WLAN device and the unknown WLAN device 102. As another example, the second intermediate location 408 of the unknown WLAN device 102 can be calculated as a combination of the first intermediate location (YTOA1) 220 and the locations 212 of the reference WLAN devices.

The TOA2 positioning unit 404 can use a Taylor series based TOA2 positioning algorithm to iteratively calculate the location of the unknown WLAN device 102 based, at least in part, on the second intermediate location 408 of the unknown WLAN device 102, the target distance calibration constant 406, and the distances 210 to the reference WLAN devices 104, 106, 108, and 112. As depicted in FIG. 4, for each reference WLAN device i, the distances 210 calculated by the distance calculation unit 202 are provided as one set of inputs to the TOA2 positioning unit 404. The second intermediate location 408 of the unknown WLAN device 102 and the target distance calibration constant 406 calculated by the TOA2 initial value calculation unit 402 are also provided as inputs to the TOA2 positioning unit 404. Additionally, the known locations 212 of the reference WLAN devices 104, 106, 108, and 112 are also provided as inputs to the TOA2 positioning unit 404. The TOA2 positioning unit 404 can then iteratively execute the expression depicted by Eq. 15 to gradually approach the real location (YTOA2) 410 of the unknown WLAN device 102 after several iterations. In accordance with Eq. 15, the TOA2 positioning unit 404 can attempt to identify the value of YTOA2 that minimizes the sum of the estimation error between A) the measured distance to the ith reference WLAN device determined in Eq. 1 based on the RTT associated with the ith reference WLAN devices and B) the estimated distance to the ith reference WLAN device based on the known location of the ith reference WLAN device. In Eq. 15, ∥Xi−YTOA2∥ represents the distance between the location (Xi) of the ith reference WLAN device and the real location (i.e., the YTOA2 location to be determined) of the unknown WLAN device 102.




  
    
      
        
          Y
          
            TOA
            ⁢
            
              
            
            ⁢
            2
          
        
        =
        
          arg
          ⁢
          
            
              min
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
              
            
            ⁢
            
              {
              
                
                  ∑
                  
                    i
                    =
                    0
                  
                  
                    N
                    -
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                
                  (
                  
                    
                      
                        d
                        ^
                      
                      i
                    
                    +
                    
                      d
                      c
                    
                    -
                    
                      
                      
                        
                          X
                          i
                        
                        -
                        
                          Y
                          
                            TOA
                            ⁢
                            
                              
                            
                            ⁢
                            2
                          
                        
                      
                      
                    
                  
                  )
                
              
              }
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        15
      
    
  




In some implementations, by assuming that there is no estimation error {e.g., by assuming that ({circumflex over (d)}i+dc−∥Xi−YTOA2∥)=0}, the TOA2 positioning unit 404 can simplify Eq. 15 and can instead solve Eq. 16 for all values of i to estimate the location of the unknown WLAN device 102.

{circumflex over (d)}i+dc=∥Xi−YTOA∥  Eq. 16

The location 410 of the unknown WLAN device 102 (YTOA2) can be expressed as the sum of the second intermediate location 408 of the unknown WLAN device (YTOA2init) and a TOA2 location error (ΔYTOA2), as depicted in Eq. 17. It is noted that in some implementations, the TOA2 location error can be represented as ΔYTOA2={ΔyTOA2s,s=0, 1, . . . , D−1}, where ΔyTOA2s represents the value of the sth dimension of the TOA2 location error. The TOA2 positioning unit 404 can substitute Eq. 17 into Eq. 16 to yield Eq. 18.

YTOA2=YTOA2init+ΔYTOA2  Eq. 17

{circumflex over (d)}i+dc=∥Xi−YTOA2init−ΔYTOA2∥, i=0, 1, . . . , N−1  Eq. 18

The TOA2 positioning unit 404 can then use Taylor series expansion and can neglect the second and higher order terms to yield Eq. 19. The TOA2 positioning unit 404 can generate an equation in accordance with Eq. 19 for each of the reference WLAN devices 104, 106, 108, and 112, and can organize the plurality of generated equations in matrix form to yield Eq. 20.




  
    
      
        
          
            
              d
              ^
            
            i
          
          +
          
            d
            c
          
        
        =
        
          
            
            
              
                X
                i
              
              -
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  
                    2
                    init
                  
                
              
            
            
          
          +
          
            
              ∑
              
                s
                =
                0
              
              
                D
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  -
                  
                    
                      
                        x
                        i
                        s
                      
                      -
                      
                        y
                        
                          TOA
                          ⁢
                          
                            
                          
                          ⁢
                          
                            2
                            init
                          
                        
                        s
                      
                    
                    
                      
                      
                        
                          X
                          i
                        
                        -
                        
                          Y
                          
                            TOA
                            ⁢
                            
                              
                            
                            ⁢
                            
                              2
                              init
                            
                          
                        
                      
                      
                    
                  
                
                )
              
              ⁢
              Δ
              ⁢
              
                
              
              ⁢
              
                y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
                s
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        19
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            A
            
              ToA
              ⁢
              
                
              
              ⁢
              2
            
          
          =
          
            
              B
              
                ToA
                ⁢
                
                  
                
                ⁢
                2
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            
              Y
              
                TOA
                ⁢
                
                  
                
                ⁢
                2
              
              T
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        20
      
    
  




In some implementations, the matrix ΔYTOA2T may be a D×1 column matrix that represents the transpose of the TOA2 location error matrix (i.e., the difference between the estimated location of the unknown WLAN device 102 and the second intermediate location of the unknown WLAN device 102) associated with the unknown WLAN device 102. As described with reference to FIG. 2, the TOA2 distance error matrix (ATOA2) can be an N×1 column matrix, where each element (i.e., each row) of the AToA2 matrix is associated with one of the reference WLAN devices and is represented in accordance with Exp. 21a. The TOA2 coefficient matrix (BToA2) can be an N×D matrix, where N represents the number of reference WLAN devices and D represents the dimensions of the coordinate system that is used to specify the location of the reference WLAN devices. Each element of the N×D BTOA2 matrix can be represented by Exp 21b.





  
    
      
        
          
            i
            th
          
          ⁢
          
            
          
          ⁢
          element
          ⁢
          
            
          
          ⁢
          of
          ⁢
          
            
          
          ⁢
          
            A
            
              ToA
              ⁢
              
                
              
              ⁢
              2
            
          
          ⁢
          
            :
          
          ⁢
          
            
          
          ⁢
          
            
              d
              ^
            
            i
          
        
        +
        
          d
          c
        
        -
        
          
          
            
              X
              i
            
            -
            
              Y
              
                TOA
                ⁢
                
                  
                
                ⁢
                
                  2
                  init
                
              
            
          
          
        
      
    
    
      
        
          Exp
          .
          
            
          
          ⁢
          21
        
        ⁢
        
          
        
        ⁢
        a
      
    
  
  
    
      
        
          
            
              (
              
                i
                ,
                s
              
              )
            
            th
          
          ⁢
          
            
          
          ⁢
          element
          ⁢
          
            
          
          ⁢
          of
          ⁢
          
            
          
          ⁢
          
            B
            
              ToA
              ⁢
              
                
              
              ⁢
              2
            
          
          ⁢
          
            :
          
        
        ⁢
        
          
        
        -
        
          
            
              x
              i
              s
            
            -
            
              y
              
                TOA
                ⁢
                
                  
                
                ⁢
                
                  2
                  init
                
              
              s
            
          
          
            
            
              
                X
                i
              
              -
              
                Y
                
                  TOA
                  ⁢
                  
                    
                  
                  ⁢
                  
                    2
                    init
                  
                
              
            
            
          
        
      
    
    
      
        
          Exp
          .
          
            
          
          ⁢
          21
        
        ⁢
        
          
        
        ⁢
        b
      
    
  




As described above with reference to FIG. 2, in one implementation, the dimension (e.g., the value of s) can vary from column to column and can remain constant throughout a particular column. The reference device under consideration (e.g., the value of i) can vary from row to row and can remain constant throughout a particular row. The TOA2 positioning unit 404 can then calculate the TOA2 location error ΔYTOA2 in accordance with Eq. 22.

ΔYTOA2T=(BToA2TBToA2)−1BToA2TAToA2  Eq. 22

The TOA2 positioning unit 404 can then use Eq. 17 to iteratively estimate the location 410 of the unknown WLAN device (YTOA2) by summing the TOA2 location error (ΔYTOA2) determined in accordance with Eq. 22 with the second intermediate location 408 of the unknown WLAN device 102 (YTOA2init). In one implementation, the location 410 of the unknown WLAN device 102 at the output of the TOA2 positioning unit 404 can be represented as YTOA2={yTOA2s,s=0, 1, . . . , D−1}, where yTOA2s represents the value of the sth dimension of the location of the unknown WLAN device.
In the next iteration, the TOA2 positioning unit 404 can use the output of the TOA2 positioning unit 404 calculated at the previous iteration as the initial location of the unknown WLAN device 102 at the input of the TOA2 positioning unit 404 {i.e., YTOA2init (at next iteration)=YTOA2 (at previous iteration)} and use Eq. 17 and Eq. 22 to determine the TOA2 location error (ΔYTOA2T) at the next iteration and the estimated location of the unknown WLAN device 102 (YTOA2) until the TOA2 location error is below a threshold error value or until a threshold number of iterations have been executed.
FIG. 5 and FIG. 6 depict a flow diagram (“flow”) 500 illustrating example operations of a hybrid TOA positioning mechanism. The flow 500 begins at block 502 in FIG. 5.
At block 502, an unknown network device determines a plurality of reference network devices of a communication network based on which to estimate the location of the unknown network device. With reference to the example of FIG. 1, the location calculation unit 110 of the unknown WLAN device 102 can identify the reference WLAN devices 104, 106, 108, and 112 based on which to estimate the location of the unknown WLAN device 102. As described above in FIG. 3, the location calculation unit 110 can select the reference WLAN devices 104, 106, 108, and 112 from a plurality of WLAN devices as those that are associated with the best performance measurements (e.g., RSSI, PER, etc.). The flow continues at block 504.
At block 504, the round-trip transit time (RTT) associated with each of the plurality of reference network device is determined. In one implementation, as described above in blocks 304-308 of FIG. 3, the location calculation unit 110 can record the time instants at which control messages were transmitted to each of the reference WLAN devices and the time instants at which the corresponding response control messages were received at the unknown WLAN device 102. Accordingly, the location calculation unit 110 can determine the RTT associated with each of the reference WLAN devices 104, 106, 108, and 112. It is noted that in other implementations, the location calculation unit 110 can employ other suitable techniques to determine the RTT associated with the reference WLAN devices 104, 106, 108, and 112. The flow continues at block 506.
At block 506, for each of the plurality of reference network devices, the distance between the unknown network device and the reference network device is calculated based, at least in part, on the RTT associated with the reference network device. For example, the location calculation unit 110 (e.g., the distance calculation unit 202 of FIG. 4) can determine the distance between the unknown WLAN device 102 and each of the plurality of reference WLAN devices, as described above with reference to Eq. 1a of FIG. 2 and with reference to block 310 of FIG. 3. The flow continues at block 508.
At block 508, the initial location of the unknown network device is calculated based, at least in part, on known locations of the reference network devices. For example, the location calculation unit 110 (e.g., the TOA1 initial location calculation unit 204 of FIG. 4) can calculate the initial location of the unknown WLAN device 102, as described above with reference to Eqs. 2a-2c of FIG. 2 and block 312 of FIG. 3. The flow continues at block 510.
At block 510, an initial distance calibration constant associated with the reference network devices is determined. For example, the location calculation unit 110 can determine the initial distance calibration constant 216 associated with the reference network devices. As described above with reference to FIG. 2, the initial distance calibration constant 216 can be zero, a previously determined value, a randomly selected value, etc. The flow continues at block 512.
At block 512, a single iteration of a first time of arrival (TOA1) positioning algorithm is executed to estimate the first intermediate location of the unknown network device and the intermediate distance calibration constant based, at least in part, on the initial location of the unknown network device, the initial distance calibration constant, and the distance to each of reference network devices. The TOA1 positioning algorithm may be used to determine the intermediate value of the distance calibration constant that may be eventually provided to a TOA2 positioning algorithm, as will be described below. In some implementations, the location calculation unit 110 (e.g., the TOA1 positioning unit 206 of FIG. 2) can execute a single iteration of the Taylor series based TOA1 positioning algorithm to estimate the first intermediate location 220 of the unknown WLAN device 102 and the intermediate distance calibration constant 218, as described above in accordance with Eq. 3a to Eq. 12 of FIG. 2. It is noted that in other implementations, the TOA1 positioning unit 206 can employ other suitable techniques to estimate the first intermediate location 220 of the unknown WLAN device 102 and the intermediate distance calibration constant 218 based, at least in part, on the initial location 214 of the unknown WLAN device 102. The flow continues at block 514 in FIG. 6.
At block 514, a target distance calibration constant is calculated based, at least in part, on the intermediate distance calibration constant, the first intermediate location of the unknown network device, the locations of the reference network devices, and the distance to each of the reference network devices. For example, the location calculation unit 110 (e.g., the TOA2 initial value calculation unit 402) can calculate the target distance calibration constant 406, as described above with reference to Eq. 13a-Eq. 13e of FIG. 4. It is noted that the TOA2 initial value calculation unit 402 can use other suitable techniques to estimate the target distance calibration constant 406. The target distance calibration constant 406 can then be used as an input parameter to a second time of arrival (TOA2) positioning algorithm for iteratively estimating the location of the unknown WLAN device 102, as will be described in block 518. The flow continues at block 516.
At block 516, a second intermediate location of the unknown network device is determined based, at least in part, on the first intermediate location of the unknown network device and the locations of the reference network devices. For example, the location calculation unit 110 (e.g., the TOA2 initial value calculation unit 402) can calculate the second intermediate location 408 of the unknown WLAN device 102 as described above in accordance with Eqs 14a-14c of FIG. 4. It is noted that the TOA2 initial value calculation unit 402 can use other suitable techniques to determine the second intermediate location 408 of the unknown WLAN device 102. The second intermediate location 408 of the unknown WLAN device 102 can then be provided as an input parameter to the TOA2 positioning algorithm for iteratively estimating the location of the unknown WLAN device 102, as will be described in block 518. The flow continues at block 518.
At block 518, the TOA2 positioning algorithm is iteratively executed to calculate the estimated location of the unknown network device based, at least in part, on the second intermediate location of the unknown network device and the target distance calibration constant. For example, as described above with reference to Eq. 15-Eq. 22 of FIG. 4, the location calculation unit 110 (e.g., the TOA2 positioning unit 404) can execute multiple iterations of a Taylor-series based TOA2 positioning algorithm to calculate the location 410 of the unknown WLAN device 102 based, at least in part, on the second intermediate location of the unknown network device (determined at block 516), the target distance calibration constant (determined at block 514), and the locations of the reference network devices. It is noted that the TOA2 initial value calculation unit 402 can employ other suitable techniques to determine the estimated location 410 of the unknown WLAN device 102. From block 518, the flow ends.
It should be understood that FIGS. 1-6 are examples meant to aid in understanding embodiments and should not be used to limit embodiments or limit scope of the claims. Embodiments may perform additional operations, fewer operations, operations in a different order, operations in parallel, and some operations differently. It is noted that in some implementations, the location calculation unit 110 can determine (for each reference WLAN device) multiple RTT measurements and multiple distance measurements. For example, the location calculation unit 110 can transmit a predetermined number of control messages to the reference WLAN device 104 and can receive a corresponding number of response control messages. Accordingly, the location calculation unit 110 can calculate the predetermined number of RTT measurements associated with the reference WLAN device 104. In some implementations, the location calculation unit 110 can determine and use the average RTT value for subsequent operations. In other implementations, the location calculation unit 110 can determine a distance to the reference WLAN device 104, for each of the RTT measurements and can determine and use the average distance to the reference WLAN device 104 for subsequent operations.
Although the Figures refer to the unknown WLAN device 102 executing operations described herein to estimate the location of the unknown WLAN device 102, embodiments are not so limited. In other embodiments, functionality for determining the location of the unknown WLAN device 102 can be executed by one or more other suitable electronic devices. In some implementations, one or more of the reference WLAN devices, a centralized server, and/or another suitable electronic device can execute some/all of the operations for determining the location of the unknown WLAN device 102. For example, the unknown WLAN device 102 can determine the RTT associated with each of the reference WLAN devices 102, 106, 108, and 112 and can provide the RTT values to the centralized server (or another device to which the processing has been offloaded). The centralized server can then calculate the location of the unknown WLAN device 102 (as described by FIGS. 1-6) and can communicate the calculated location to the unknown WLAN device 102.
Embodiments may take the form of an entirely hardware embodiment, a software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining software and hardware aspects that may all generally be referred to herein as a “circuit,” “module” or “system.” Furthermore, embodiments of the inventive subject matter may take the form of a computer program product embodied in any tangible medium of expression having computer usable program code embodied in the medium. The described embodiments may be provided as a computer program product, or software, that may include a machine-readable medium having stored thereon instructions, which may be used to program a computer system (or other electronic device(s)) to perform a process according to embodiments, whether presently described or not, since every conceivable variation is not enumerated herein. A machine-readable medium includes any mechanism for storing or transmitting information in a form (e.g., software, processing application) readable by a machine (e.g., a computer). A machine-readable medium may be a machine-readable storage medium, or a machine-readable signal medium. A machine-readable storage medium may include, for example, but is not limited to, magnetic storage medium (e.g., floppy diskette); optical storage medium (e.g., CD-ROM); magneto-optical storage medium; read only memory (ROM); random access memory (RAM); erasable programmable memory (e.g., EPROM and EEPROM); flash memory; or other types of tangible medium suitable for storing electronic instructions. A machine-readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, an electrical, optical, acoustical, or other form of propagated signal (e.g., carrier waves, infrared signals, digital signals, etc.). Program code embodied on a machine-readable signal medium may be transmitted using any suitable medium, including, but not limited to, wireline, wireless, optical fiber cable, RF, or other communications medium.
Computer program code for carrying out operations of the embodiments may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C++ or the like and conventional procedural programming languages, such as the “C” programming language or similar programming languages. The program code may execute entirely on a user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN), a personal area network (PAN), or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider).
FIG. 7 is a block diagram of one embodiment of an electronic device 700 including a mechanism for determining the location of the electronic device in a wireless communication network. In some implementations, the electronic device 700 may be one of a laptop computer, a tablet computer, a netbook, a mobile phone, a smart appliance, a gaming console, or other electronic systems comprising wireless communication capabilities. The electronic device 700 includes a processor unit 702 (possibly including multiple processors, multiple cores, multiple nodes, and/or implementing multi-threading, etc.). The electronic device 700 includes a memory unit 706. The memory unit 706 may be system memory (e.g., one or more of cache, SRAM, DRAM, zero capacitor RAM, Twin Transistor RAM, eDRAM, EDO RAM, DDR RAM, EEPROM, NRAM, RRAM, SONOS, PRAM, etc.) or any one or more of the above already described possible realizations of machine-readable media. The electronic device 700 also includes a bus 710 (e.g., PCI, ISA, PCI-Express, HyperTransport®, InfiniBand®, NuBus, AHB, AXI, etc.), and network interfaces 704 that include at least one of a wireless network interface (e.g., a WLAN interface, a Bluetooth® interface, a WiMAX interface, a ZigBee® interface, a Wireless USB interface, etc.) and a wired network interface (e.g., a powerline communication interface, an Ethernet interface, etc.).
The electronic device 700 also includes a communication unit 708. The communication unit 708 comprises a location calculation unit 712. In some implementations, the location calculation unit 712 can execute a single iteration of a first TOA positioning algorithm to estimate the location of the electronic device 700 based, at least in part, on a calculated initial location of the electronic device 700, the distance between the electronic device 700 and each of a plurality of reference WLAN devices, and the location of each of the plurality of reference WLAN devices, as described with reference to FIGS. 1-3. In other implementations, the location calculation unit 712 can execute a hybrid TOA positioning algorithm to estimate the location of the electronic device 700, as described above with reference to FIGS. 1 and 4-6. Any one of these functionalities may be partially (or entirely) implemented in hardware and/or on the processor unit 702. For example, the functionality may be implemented with an application specific integrated circuit, in logic implemented in the processor unit 702, in a co-processor on a peripheral device or card, etc. Further, realizations may include fewer or additional components not illustrated in FIG. 7 (e.g., video cards, audio cards, additional network interfaces, peripheral devices, etc.). The processor unit 702, the memory unit 706, and the network interfaces 704 are coupled to the bus 710. Although illustrated as being coupled to the bus 710, the memory unit 706 may be coupled to the processor unit 702.
While the embodiments are described with reference to various implementations and exploitations, it will be understood that these embodiments are illustrative and that the scope of the inventive subject matter is not limited to them. In general, techniques for a TOA based positioning system or techniques for a hybrid TOA based positioning system as described herein may be implemented with facilities consistent with any hardware system or hardware systems. Many variations, modifications, additions, and improvements are possible.
Plural instances may be provided for components, operations, or structures described herein as a single instance. Finally, boundaries between various components, operations, and data stores are somewhat arbitrary, and particular operations are illustrated in the context of specific illustrative configurations. Other allocations of functionality are envisioned and may fall within the scope of the inventive subject matter. In general, structures and functionality presented as separate components in the exemplary configurations may be implemented as a combined structure or component. Similarly, structures and functionality presented as a single component may be implemented as separate components. These and other variations, modifications, additions, and improvements may fall within the scope of the inventive subject matter.

## Claims

1. A method comprising:
determining, at a wireless network device of a communication network, a round trip transit time between the wireless network device and each of a plurality of reference wireless network devices of the communication network;
for each of the plurality of reference wireless network devices, determining a distance between the wireless network device and the reference wireless network device based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device;
obtaining location coordinates associated with each of the plurality of reference wireless network devices; and
performing time of arrival (TOA) computations for estimating a location of the wireless network device by executing an iteration of a Taylor series based TOA positioning algorithm based, at least in part, on an initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, and an initial distance calibration constant, wherein the initial location is based, at least in part, on an average of the location coordinates associated with each of the plurality of reference wireless network devices.

2. The method of claim 1, wherein the initial distance calibration constant is one of a zero value, a predetermined value, and a dynamically selected random value.

3. The method of claim 1, wherein said performing the TOA computations for estimating the location of the wireless network device further comprises determining a distance calibration constant based, at least in part, on the initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, the initial distance calibration constant, and the location coordinates associated with each of the plurality of reference wireless network devices in response to said executing the single iteration of the Taylor series based TOA positioning algorithm.

4. The method of claim 1, further comprising:
calculating the initial location of the wireless network device as a non-weighted average of the location coordinates associated with each of the plurality of reference wireless network devices.

5. The method of claim 1, further comprising:
calculating the initial location of the wireless network device as a weighted average of the location coordinates of at least a subset of the plurality of reference wireless network devices.

6. The method of claim 5, wherein weights used in the weighted average are based upon a confidence metric associated with the location coordinates associated with the plurality of reference wireless network devices.

7. The method of claim 1, further comprising:
identifying the plurality of reference wireless network devices based, at least in part, on analyzing one or more performance measurements associated with at least each of the plurality of reference wireless network devices.

8. The method of claim 1, wherein said for each of the plurality of reference wireless network devices, said determining the distance between the wireless network device and the reference wireless network device based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device further comprises:
for each of the plurality of reference wireless network devices,
recording a first time instant at which a first control message is transmitted from the wireless network device to the reference wireless network device;
recording a second time instant at which a second control message is received at the wireless network device from the reference wireless network device in response to the first control message;
subtracting the first time instant from the second time instant to determine the round trip transit time between the wireless network device and the reference wireless network device; and
multiplying round trip transit time between the wireless network device and the reference wireless network device by a speed of light factor to yield the distance between the wireless network device and the reference wireless network device.

9. The method of claim 1, wherein the wireless network device and the plurality of reference wireless network devices comprise wireless local area network (WLAN) communication capabilities.

10. The method of claim 1, wherein said obtaining location coordinates associated with each of the plurality of reference wireless network devices comprises:
transmitting a request for the location coordinates to each of the plurality of reference wireless network devices; and
receiving the location coordinates from each of the plurality of reference wireless network devices.

11. The method of claim 10, wherein the request for the location coordinates is included in a control message that is also used to determine the round trip transmit time between the wireless network device and each of the plurality of reference wireless network devices.

12. The method of claim 1, wherein said obtaining location coordinates associated with each of the plurality of reference wireless network devices comprises:
transmitting a request for the location coordinates to a centralized server; and
receiving the location coordinates associated with each of the plurality of reference wireless network devices from the centralized server.

13. The method of claim 1, wherein said obtaining location coordinates associated with each of the plurality of reference wireless network devices comprises:
receiving a periodic broadcast message from each of the plurality of reference wireless network devices, the periodic broadcast message including the location coordinates from a corresponding one of the plurality of reference wireless network devices.

14. The method of claim 1, wherein the location coordinates associated with each of the plurality of reference wireless network devices comprises one of two-dimensional (2-D) coordinates, three-dimensional (3-D) coordinates, geographic latitude and longitude coordinates, or spherical coordinates.

15. A method comprising:
determining, at a wireless network device of a communication network, a round trip transit time between the wireless network device and each of a plurality of reference wireless network devices of the communication network;
for each of the plurality of reference wireless network devices, determining a distance between the wireless network device and the reference wireless network device based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device;
determining a location error matrix, wherein each element of the location error matrix represents a difference between an estimated location coordinate of the wireless network device and a corresponding initial location coordinate of an initial location of the wireless network device;
determining a distance error matrix, wherein each element of the distance error matrix represents a difference between a measured distance between the wireless network device and one of the plurality of reference wireless network devices and an estimated distance between the wireless network device and the one of the plurality of reference wireless network devices, wherein the measured distance is determined based, at least in part, on the round trip transit time between the wireless network device and the one of the plurality of reference wireless network devices, wherein the estimated distance is determined based, at least in part, on the location of the one of the plurality of reference wireless network devices and the initial location of the wireless network device;
determining a coefficient matrix as the product of the distance error matrix and an inverse of a transpose of the location error matrix;
calculating the initial location of the wireless network device based, at least in part, on maximizing the determinant of the product of the coefficient matrix and a transpose of the coefficient matrix; and
performing time of arrival (TOA) computations for estimating a location of the wireless network device based, at least in part, on the initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, and an initial distance calibration constant.

16. A wireless network device comprising:
a network interface; and
a location calculation unit coupled with the network interface, the location calculation unit operable to:
determine a round trip transit time between the wireless network device and each of a plurality of reference wireless network devices of a communication network;
for each of the plurality of reference wireless network devices, determine a distance between the wireless network device and the reference wireless network device based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device;
obtain location coordinates associated with each of the plurality of reference wireless network devices; and
perform time of arrival (TOA) computations to estimate a location of the wireless network device by executing an iteration of a Taylor series based TOA positioning algorithm based, at least in part, on an initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, and an initial distance calibration constant, wherein the initial location is based, at least in part, on an average of the location coordinates associated with each of the plurality of reference wireless network devices.

17. The wireless network device of claim 16, wherein the location calculation unit operable to perform the TOA computations to estimate the location of the wireless network device further comprises the location calculation unit operable to determine a distance calibration constant based, at least in part, on the initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, the initial distance calibration constant, and the location coordinates associated with each of the plurality of reference wireless network devices in response to the location calculation unit executing the single iteration of the Taylor series based TOA positioning algorithm.

18. The wireless network device of claim 11, wherein the location calculation unit is further operable to:
calculate the initial location of the wireless network device as a non-weighted average of the location coordinates associated with each of the plurality of reference wireless network devices, or
calculate the initial location of the wireless network device as a weighted average of the location of at least a subset of the plurality of reference wireless network devices.

19. A wireless network device comprising:
a network interface; and
a location calculation unit coupled with the network interface, the location calculation unit operable to:
determine a round trip transit time between the wireless network device and each of a plurality of reference wireless network devices of a communication network;
for each of the plurality of reference wireless network devices, determine a distance between the wireless network device and the reference wireless network device based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device;
determine a location error matrix, wherein each element of the location error matrix represents a difference between an estimated location coordinate of the wireless network device and a corresponding initial location coordinate of an initial location of the wireless network device;
determine a distance error matrix, wherein each element of the distance error matrix represents a difference between a measured distance between the wireless network device and one of the plurality of reference wireless network devices and an estimated distance between the wireless network device and the one of the plurality of reference wireless network devices, wherein the measured distance is determined based, at least in part, on the round trip transit time between the wireless network device and the one of the plurality of reference wireless network devices, wherein the estimated distance is determined based, at least in part, on the location of the one of the plurality of reference wireless network devices and the initial location of the wireless network device;
determine a coefficient matrix as the product of the distance error matrix and an inverse of a transpose of the location error matrix;
calculate the initial location of the wireless network device based, at least in part, on maximizing the determinant of the product of the coefficient matrix and a transpose of the coefficient matrix; and
perform time of arrival (TOA) computations to estimate a location of the wireless network device based, at least in part, on the initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, and an initial distance calibration constant.

20. One or more non-transitory machine-readable storage media having instructions stored therein, which when executed by one or more processors causes the one or more processors to perform operations that comprise:
determining, at a wireless network device of a communication network, a round trip transit time between the wireless network device and each of a plurality of reference wireless network devices of the communication network;
for each of the plurality of reference wireless network devices, determining a distance between the wireless network device and the reference wireless network device based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device;
obtaining location coordinates associated with each of the plurality of reference wireless network devices; and
performing time of arrival (TOA) computations for estimating a location of the wireless network device by executing an iteration of a Taylor series based TOA positioning algorithm based, at least in part, on an initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, and an initial distance calibration constant, wherein the initial location is based, at least in part, on an average of the location coordinates associated with each of the plurality of reference wireless network devices.

21. The non-transitory machine-readable storage media of claim 20, wherein said operation of performing the TOA computations for estimating the location of the wireless network device further comprises determining a distance calibration constant based, at least in part, on the initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, the initial distance calibration constant, and the location coordinates associated with each of the plurality of reference wireless network devices in response to said executing the single iteration of the Taylor series based TOA positioning algorithm.

22. The non-transitory machine-readable storage media of claim 20, wherein the instructions stored therein, which when executed by one or more processors causes the one or more processors to further perform operations that comprise:
calculating the initial location of the wireless network device as a non-weighted average of the location coordinates associated with each of the plurality of reference wireless network devices, or
calculating the initial location of the wireless network device as a weighted average of the location of at least a subset of the plurality of reference wireless network devices.

23. One or more machine-readable storage media having instructions stored therein, which when executed by one or more processors causes the one or more processors to perform operations that comprise:
determining, at a wireless network device of a communication network, a round trip transit time between the wireless network device and each of a plurality of reference wireless network devices of the communication network;
for each of the plurality of reference wireless network devices, determining a distance between the wireless network device and the reference wireless network device based, at least in part, on the round trip transit time between the wireless network device and the reference wireless network device;
determining a location error matrix, wherein each element of the location error matrix represents a difference between an estimated location coordinate of the wireless network device and a corresponding initial location coordinate of an initial location of the wireless network device;
determining a distance error matrix, wherein each element of the distance error matrix represents a difference between a measured distance between the wireless network device and one of the plurality of reference wireless network devices and an estimated distance between the wireless network device and the one of the plurality of reference wireless network devices, wherein the measured distance is determined based, at least in part, on the round trip transit time between the wireless network device and the one of the plurality of reference wireless network devices, wherein the estimated distance is determined based, at least in part, on the location of the one of the plurality of reference wireless network devices and the initial location of the wireless network device;
determining a coefficient matrix as the product of the distance error matrix and an inverse of a transpose of the location error matrix; and calculating the initial location of the wireless network device based, at least in part, on maximizing the determinant of the product of the coefficient matrix and a transpose of the coefficient matrix; and
performing time of arrival (TOA) computations for estimating a location of the wireless network device based, at least in part, on the initial location of the wireless network device, the distance between the wireless network device and each of the plurality of reference wireless network devices, and an initial distance calibration constant.