System and method for performing time difference of arrival location without requiring a common time base or clock calibration

[0001] This application is a continuation of U.S. application Ser. No. 11/367,144, filed on Mar. 3, 2006 now U.S. Pat. No. 7,477,191.

[0002] The present invention relates generally to a system and method for performing Time Difference Of Arrival (TDOA) and more specifically to a system and method for performing TDOA without requiring a common time base or clock calibration.

[0003] Time Difference Of Arrival (TDOA) refers to determining the position of a wireless transceiver (e.g. radio) by transmitting a signal from the wireless transceiver and then estimating the signal arrival time at multiple reference receivers, for example access points ‘APs’ in a Wireless Local Area Network (WLAN) system with known positions. The TDOA algorithm then compares the relative arrival time of the signal at the various APs to determine the radio's location. If the various APs do not share a common time base, clock calibration is required to determine the relative difference between the time bases of the APs. Performing clock calibration in practice is problematic for several reasons. First, it consumes time that could be used for data transmission. Second, the clock frequencies of the various APs will vary with time and therefore require periodic recalibration. Third, any inaccuracy in the clock frequency calibration will result in location estimation error.

[0004] In accordance with an aspect of the present invention, there is disclosed herein a TDOA implementation that eliminates the need for clock calibration and therefore also eliminates the problems associated with performing clock calibration. The implementation includes a method and an apparatus configured to implement the method. The method relies on a packet transmitted from a reference wireless device with a known propagation delay, which serves as a common reference point for all TDOA estimates.

[0005] In accordance with an aspect of the present invention, there is described herein a method for determining a time difference of arrival from a signal sent by a third wireless device that is received by a first wireless device and a second wireless device. The method assumes that a signal was also received by the second wireless device from the first wireless device and that the propagation delay between the first wireless device and the second wireless device is known. The method comprises determining a first difference, the first difference being the difference between when a packet transmitted by a first wireless device was received by a second wireless device. The method further comprises determining a second difference, the second difference being the difference between when a packet transmitted by a first wireless device was received and when a packet from a third wireless device was received by the second wireless device. The method also comprises determining the relative frequency difference between the first wireless device and the second wireless device. The method still further comprises determining a product, the product being the product of the relative frequency difference and the second difference. The method still further comprises determining a time difference of arrival from the difference between the first difference and the product. Aspects of the present invention further comprise an apparatus for implementing methodologies of the present invention.

[0006] Still other objects of the present invention will become readily apparent to those skilled in this art from the following description wherein there is shown and described a preferred embodiment of this invention, simply by way of illustration of at least one of the best modes best suited to carry out the invention. As it will be realized, the invention is capable of other different embodiments and its several details are capable of modifications in various obvious aspects all without departing from the invention. Accordingly, the drawing and descriptions will be regarded as illustrative in nature and not as restrictive.

[0014] Throughout this description, the preferred embodiment and examples shown should be considered as exemplars, rather than limitations, of the present invention. An aspect of the present invention is for a method for performing Time Difference Of Arrival (TDOA) that eliminates the need for a common time base or clock calibration. The method relies on receiving a packet transmitted from a reference wireless receiver with a known propagation delay, which serves as a common reference point for TDOA estimates.

[0015] Referring to FIG. 1, there is illustrated a network 100 for implementing TDOA. The network comprises three access points (APs) AP1 102, AP2 104 and AP3 106 that can be employed to triangulate the position of client 118. AP1 102, AP2 104 and AP3 106 are coupled by a distribution network 108. Distribution network 108 may suitably comprise a wired segment, wireless segment or a combination of wired and wireless segments. Distribution network 108 may suitably comprise other network devices including but not limited to additional access points, an authentication authorization and accounting (AAA) server, a wireless domain server (WDS), or other computing devices (not shown). Coupled to network 108 is TDOA Logic 116. “Logic”, as used herein, includes but is not limited to hardware, firmware, software and/or combinations of each to perform a function(s) or an action(s), and/or to cause a function or action from another component. For example, based on a desired application or need, logic may include a software controlled microprocessor, discrete logic such as an application specific integrated circuit (ASIC), a programmable/programmed logic device, memory device containing instructions, or the like, or combinational logic embodied in hardware. Logic may also be fully embodied as software. TDOA Logic 116 is configured to compute TDOA based on signals received by at least two of AP1 102, AP2 106 and AP3 106. Although TDOA logic 116 is illustrated as being a standalone network component, it is also contemplated in alternative embodiments that TDOA logic can be co-located with any network component, for example AP2 104 and/or AP3 106.

[0016] Client 118 is illustrated as being in two-way communication with AP1 102 via path 110. AP2 106 and AP3 106 also can receive signals from client 118 as illustrated by 114 and 115 respectively. Furthermore, AP2 106 can receive signals from AP1 102 as illustrated by path 112. In the example used herein, the TDOA method is described using two APs, e.g. AP1 102 and AP2 104. The use of two APs for describing the method is for ease of illustration and those skilled in the art should readily appreciate that aspects of the present invention are suitable for any number of additional APs or wireless devices capable of receiving a wireless signal from a reference AP (e.g. AP1 102) and client 118.

[0017] An example of a method in accordance with an aspect of the present invention will be now be described with reference to FIG. 1. For example, AP1 102 transmits a packet to the client to be located (i.e. client 118). Client 118 transmits a packet in response. It is assumed that AP2 104 and AP3 106 (or any other AP used to calculate TDOA for client 118) knows the propagation delay between itself and AP1 102, has some time base that it uses to report TDOA, knows the time of arrival (TOA) of the packet from AP1 102 and knows the TOA of the packet from client 118.

[0018] In prior art systems, unless the APs (AP1 102, AP2 104 and AP3 106) share a common time base or calibration information is available to “align” the time bases of the APs, TOA estimates cannot be directly compared between APs. However, in accordance with an aspect of the present invention, because AP2 104 (and/or AP3) knows the propagation delay between itself and AP1 102 and knows the time of arrival of the packet from AP1 102, each AP can determine the time of departure (TOD) for the packet sent by AP1 102, which provides a common reference point. By subtracting the TOD of the packet from AP1 102 from the TOA of the packet sent from client 118, a TOA estimate that is relative to the same reference point can be obtained and thus TOA estimates can be directly compared between APs.

[0019] The propagation delay between each of the APs can either be estimated (e.g. the location of the AP is known) or can be predetermined by measuring the distance between APs or can be calculated dynamically, see e.g. U.S. Patent Application 2004/0203877 to Golden, hereby incorporated by reference, for a method of calculating propagation delay. It should also be appreciated that it is not necessary that AP1 102 transmits a packet prior to the packet sent by client 118. For example, radio tags transmit periodically rather than in response to an AP, AP1 102 could transmit immediately after it observes a client packet. In a preferred embodiment, each AP used by the TDOA algorithm is able to receive signals from both AP1 102 and client 118.

[0020] Referring to FIG. 2, with continued reference to FIG. 1, there is illustrated a timing diagram indicating how the TDOA for client 118 is computed between AP1 102 and AP2 104. In FIG. 2, t_{AP1 TX} is the start of packet time for the packet transmitted from AP1 102, t_{C TX} is the start of packet time for the packet transmitted from the client, t_{AP1 TX+DAP1:AP2} is the time when the packet from AP1 102 arrives at AP2 104, t_{C TX+DAP1:C} is the time when the packet from the client arrives at AP1 102, t_{C TX+DAP2:C} is the time when the packet from the client arrives at AP2 104, D_AP1:AP2 is the propagation delay between AP1 102 and AP2 104, D_AP1:C is the propagation delay between AP1 102 and the client 118, and D_AP2:C is the propagation delay between AP2 104 and the client 118. An aspect of the present invention is to determine the TDOA (D_AP1:C−D_AP2:C) between AP1 102 and AP2 104.

[0021] If both APs (AP1 102 and AP2 104) share a common time base, determining the TDOA is as simple as subtracting t_{C TX+DAP2:C} from t_{C TX+DAP1:C}. In a typical WLAN system, however, each AP has its own time base derived from its own reference clock. A simplified schematic of an AP 300 is shown in FIG. 3. AP 300 is suitable adapted to be employed as AP1 102 (FIG. 1), AP 104 (FIG. 1) and AP 106 (FIG. 1).

[0022] AP 300 sends and receives wireless signals via antenna 302. Switch 304 switches to couple either the transmit circuit or the receive circuit to antenna 302. As illustrated, transmit circuit comprises transmit buffer 316 and a digital to analog converter (DAC) 318. Transmit circuit further comprises any other circuitry (not shown) desired for amplification, encoding, modulation and frequency conversion, which are not essential for explaining aspects the present invention. The receive circuit comprises an analog to digital converter (ADC) 306, a clock (CLK) 308 coupled to ADC 308 and a Time of Arrival (TOA) Algorithm 310 coupled to ADC 306 and CLK 308 for determining a TOA for a packet received via antenna 302. TOA algorithm outputs a signal 314 indicative of a timing estimate of the TOA. In addition, receiver circuit further comprises a receive (RX) buffer 312 and any other circuitry (not shown) desired for amplification, decoding, demodulating and frequency conversion, which are not essential for explaining aspects of the present invention.

[0023] In operation, when a signal is received via antenna 302, switch 304 routes the packet to ADC 306, which converts the signal from an analog signal to a digital signal. The digital signal is processed by TOA algorithm 310 employing CLK 308 to determine a TOA for the signal.

[0024] Referring to FIG. 3, with continued reference to FIGS. 2 and 1, an AP 300 has a time base related to “absolute” time by:

[0025] $\begin{array}{cc}{t}^{\mathrm{CLK}}=\left(t-t_{\mathrm{CLK}\mathrm{ON}}\right)\left(1+\frac{f_{\mathrm{CLK}}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)& \mathrm{Equation}1\end{array}$

[0026] where t^CLK denotes time as perceived by AP 300, t denotes absolute time, t_{CLK ON} denotes the absolute time at which the clock and timer of the AP 300 were powered on and thus the timer started tracking time, f_CLK is the clock frequency, and f_nom denotes the nominal (or expected) clock frequency.

[0027] From AP1 102, we have time estimates of when the packet from AP1 102 was transmitted (t_{AP1 TX}^CLK1) and when the client packet was received (t_CTX+DAP1:C^CLK1) by AP2 104. Using Equation 1,

[0028] $\begin{array}{cc}{t}_{\mathrm{AP}1\mathrm{TX}}^{\mathrm{CLK}1}=\left(t_{\mathrm{AP}1\mathrm{TX}}-t_{\mathrm{CLK}1\mathrm{ON}}\right)\left(1+\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)\mathrm{and}& \mathrm{Equation}2\\ t_{\mathrm{CTX}+D_{\mathrm{AP}1:C}}^{\mathrm{CLK}1}=\left(t_{\mathrm{CTX}}+D_{\mathrm{AP}1:C}-t_{\mathrm{CLK}1\mathrm{ON}}\right)\left(1+\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right).& \mathrm{Equation}3\end{array}$

[0029] Subtracting Equation 2 from Equation 3 gives

[0030] $\begin{array}{cc}\left(t_{\mathrm{CTX}}+D_{\mathrm{AP}1:C}-t_{\mathrm{AP}1\mathrm{TX}}\right)\left(1+\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right).& \mathrm{Equation}4\end{array}$

[0031] From AP2 104, we have time estimates of when the packet from AP1 102 was received (t_{AP1 TX+DAP1:AP2}^CLK2) and when the client 118 packet was received (t_{C TX+DAP2:C}^CLK2). Using Equation 1,

[0032] $\begin{array}{cc}{t}_{\mathrm{AP}1\mathrm{TX}+D_{\mathrm{AP}1:\mathrm{AP}2}}^{\mathrm{CLK}2}=\left(t_{\mathrm{AP}1\mathrm{TX}}+D_{\mathrm{AP}1:\mathrm{AP}2}-t_{\mathrm{CLK}2\mathrm{ON}}\right)\left(1+\frac{f_{\mathrm{CLK}2}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)\mathrm{and}& \mathrm{Equation}5\\ t_{C\mathrm{TX}+D_{\mathrm{AP}2:C}}^{\mathrm{CLK}2}=\left(t_{\mathrm{CTX}}+D_{\mathrm{AP}2:C}-t_{\mathrm{CLK}2\mathrm{ON}}\right)\left(1+\frac{f_{\mathrm{CLK}2}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right).& \mathrm{Equation}6\end{array}$

[0033] Subtracting Equation 5 from Equation 6 gives

[0034] $\begin{array}{cc}\left(t_{\mathrm{CTX}}+D_{\mathrm{AP}2:C}-t_{\mathrm{AP}1\mathrm{TX}}-D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(1+\frac{f_{\mathrm{CLK}2}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right).& \mathrm{Equation}7\end{array}$

[0035] Subtracting Equation 7 from Equation 4 gives

[0036] $\begin{array}{cc}{D}_{\mathrm{AP}1:C}-D_{\mathrm{AP}2:C}+D_{\mathrm{AP}1:\mathrm{AP}2}+\left(t_{\mathrm{CTX}}-t_{\mathrm{AP}1\mathrm{TX}}\right)\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{CLK}2}}{f_{\mathrm{nom}}}\right)+D_{\mathrm{AP}1:C}\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)+\left(D_{\mathrm{AP}1:\mathrm{AP}2}-D_{\mathrm{AP}2:C}\right)\left(\frac{f_{\mathrm{CLK}2}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right).& \mathrm{Equation}8\end{array}$

[0037] which consists of the desired result (D_AP1:C−D_AP2:C) plus several error terms given by

[0038] $\begin{array}{cc}\mathrm{Error}=D_{\mathrm{AP}1:\mathrm{AP}2}+\left(t_{\mathrm{CTX}}-t_{\mathrm{AP}1\mathrm{TX}}\right)\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{CLK}2}}{f_{\mathrm{nom}}}\right)+D_{\mathrm{AP}1:C}\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)+\left(D_{\mathrm{AP}1:\mathrm{AP}2}-D_{\mathrm{AP}2:C}\right)\left(\frac{f_{\mathrm{CLK}2}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)& \mathrm{Equation}9\end{array}$

[0039] The first error term (D_AP1:AP2) can be determined in one of two ways. The first method is to estimate it from the known AP (102, 104) positions according to

[0040] $\begin{array}{cc}{D}_{\mathrm{AP}1:\mathrm{AP}2}=\frac{\sqrt{{\left(X_{\mathrm{AP}1}-X_{\mathrm{AP}2}\right)}^2+{\left(Y_{\mathrm{AP}1}-Y_{\mathrm{AP}2}\right)}^2+{\left(Z_{\mathrm{AP}1}-Z_{\mathrm{AP}2}\right)}^2}}{c}& \mathrm{Equation}10\end{array}$
where (X_AP1, Y_AP1, Z_AP1) are the Cartesian coordinates of AP1 102, (X_AP2, Y_AP2, Z_AP2) are the Cartesian coordinates of AP2 104, and c is the speed of light in free space. One disadvantage of this approach is that the propagation delay could actually be larger than that predicted by Equation 10 due to multipath. A more accurate method is described in US Patent Application Publication No. 2004/0203877, which measures the actual propagation delay. Subtracting D_AP1:AP2 from Equation 9 leaves us with an overall error of

[0041] $\begin{array}{cc}\mathrm{Error}=\left(t_{\mathrm{CTX}}-t_{\mathrm{AP}1\mathrm{TX}}\right)\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{CLK}2}}{f_{\mathrm{nom}}}\right)+D_{\mathrm{AP}1:C}\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)+\left(D_{\mathrm{AP}1:\mathrm{AP}2}-D_{\mathrm{AP}2:C}\right)\left(\frac{f_{\mathrm{CLK}2}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)& \mathrm{Equation}11\end{array}$

[0042] For a typical WLAN system where the maximum distance between radios is on the order of 1000 feet and the clock errors are confined to no more than 20 parts per million (PPM), the error term associated with

[0043] $\begin{array}{cc}{D}_{\mathrm{AP}1:C}\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)+\left(D_{\mathrm{AP}1:\mathrm{AP}2}-D_{\mathrm{AP}2:C}\right)\left(\frac{f_{\mathrm{CLK}2}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)& \mathrm{Equation}12\end{array}$
should be less than 0.1 nS. The error term associated with

[0044] $\begin{array}{cc}\left(t_{\mathrm{CTX}}-t_{\mathrm{AP}1\mathrm{TX}}\right)\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{CLK}2}}{f_{\mathrm{nom}}}\right),& \mathrm{Equation}13\end{array}$
however, would be 4 nS for t_{C TX}−t_{AP1 TX}=100 μS (typical for DSSS/CCK operation) if the two clocks have frequency differences of 40 PPM (f_{CLK1 Err}=+20 PPM and f_{CLK2 Err}=−20 PPM).

[0045] To avoid the error term given by Equation 13, the relative frequency difference between AP1 102 and AP2 104 can be computed. This is typically done as part of the demodulation process when a packet is received. Thus, AP2 104 computes the relative frequency difference between f_CLK1 and f_CLK2 when it receives the packet from AP1 102. Assume the relative frequency estimate is given by

[0046] $\begin{array}{cc}\frac{f_{\mathrm{CLK}1}}{f_{\mathrm{CLK}2}}+\frac{f_{\mathrm{Err}}}{f_{\mathrm{nom}}}& \mathrm{Equation}14\end{array}$
where f_Err denotes the frequency estimation error term. Multiplying Equation 7 and Equation 14 gives

[0047] $\begin{array}{cc}\left(t_{\mathrm{CTX}}{D}_{\mathrm{AP}2:C}-t_{\mathrm{AP}1\mathrm{TX}}-D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(1+\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}+\frac{f_{\mathrm{Err}}{f}_{\mathrm{CLK}2}}{f_{\mathrm{nom}}^2}\right).& \mathrm{Equation}15\end{array}$

[0048] Subtracting Equation 15 from Equation 4 gives

[0049] $\begin{array}{cc}{D}_{\mathrm{AP}1:C}-D_{\mathrm{AP}2:C}+D_{\mathrm{AP}1:\mathrm{AP}2}+\left(D_{\mathrm{AP}1:C}-D_{\mathrm{AP}2:C}+D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)-\left(t_{\mathrm{CTX}}+D_{\mathrm{AP}2:C}-t_{\mathrm{AP}1\mathrm{TX}}-D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(\frac{f_{\mathrm{Err}}{f}_{\mathrm{CLK}2}}{f_{\mathrm{nom}}^2}\right).& \mathrm{Equation}16\end{array}$

[0050] The error term is now given by

[0051] $\begin{array}{cc}\mathrm{Error}=D_{\mathrm{AP}1:\mathrm{AP}2}+\left(D_{\mathrm{AP}1:C}-D_{\mathrm{AP}2:C}+D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)-\left(t_{\mathrm{CTX}}+D_{\mathrm{AP}2:C}-t_{\mathrm{AP}1\mathrm{TX}}-D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(\frac{f_{\mathrm{Err}}{f}_{\mathrm{CLK}2}}{f_{\mathrm{nom}}^2}\right).& \mathrm{Equation}17\end{array}$

[0052] As previously discussed, D_AP1:AP2 can either be estimated or measured. The error term given by

[0053] $\begin{array}{cc}\left(D_{\mathrm{AP}1:C}-D_{\mathrm{AP}2:C}+D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(\frac{f_{\mathrm{CLK}1}-f_{\mathrm{nom}}}{f_{\mathrm{nom}}}\right)& \mathrm{Equation}18\end{array}$
should be less than 0.1 nS similar to the error term of Equation 12. The frequency estimation error is typically limited to several kilohertz and the error term associated with

[0054] $\begin{array}{cc}\left(t_{\mathrm{CTX}}+D_{\mathrm{AP}2:C}-t_{\mathrm{AP}1\mathrm{TX}}-D_{\mathrm{AP}1:\mathrm{AP}2}\right)\left(\frac{f_{\mathrm{Err}}{f}_{\mathrm{CLK}2}}{f_{\mathrm{nom}}^2}\right)& \mathrm{Equation}19\end{array}$
is several orders of magnitude less than the error term associated with Equation 13.

[0055] FIG. 4 is a block diagram that illustrates a computer system 400 upon which an embodiment of the invention may be implemented. Computer system 400 includes a bus 402 or other communication mechanism for communicating information and a processor 404 coupled with bus 402 for processing information. Computer system 400 also includes a main memory 406, such as random access memory (RAM) or other dynamic storage device coupled to bus 402 for storing information and instructions to be executed by processor 404. Main memory 406 also may be used for storing a temporary variable or other intermediate information during execution of instructions to be executed by processor 404. Computer system 400 further includes a read only memory (ROM) 408 or other static storage device coupled to bus 402 for storing static information and instructions for processor 404. A storage device 410, such as a magnetic disk or optical disk, is provided and coupled to bus 402 for storing information and instructions.

[0056] The invention is related to the use of computer system 400 for performing TDOA location without requiring a common time base or clock calibration. According to one embodiment of the invention, performing TDOA location without requiring a common time base or clock calibration is provided by computer system 400 in response to processor 404 executing one or more sequences of one or more instructions contained in main memory 406. Such instructions may be read into main memory 406 from another computer-readable medium, such as storage device 410. Execution of the sequence of instructions contained in main memory 406 causes processor 404 to perform the process steps described herein. One or more processors in a multi-processing arrangement may also be employed to execute the sequences of instructions contained in main memory 406. In alternative embodiments, hard-wired circuitry may be used in place of or in combination with software instructions to implement the invention. Thus, embodiments of the invention are not limited to any specific combination of hardware circuitry and software.

[0057] The term “computer-readable medium” as used herein refers to any medium that participates in providing instructions to processor 404 for execution. Such a medium may take many forms, including but not limited to non-volatile media, volatile media, and transmission media. Non-volatile media include for example optical or magnetic disks, such as storage device 410. Volatile media include dynamic memory such as main memory 406. Transmission media include coaxial cables, copper wire and fiber optics, including the wires that comprise bus 402. Transmission media can also take the form of acoustic or light waves such as those generated during radio frequency (RF) and infrared (IR) data communications. Common forms of computer-readable media include for example floppy disk, a flexible disk, hard disk, magnetic cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, an EPROM, a FLASHPROM, any other memory chip or cartridge, a carrier wave as described hereinafter, or any other medium from which a computer can read.

[0058] Various forms of computer-readable media may be involved in carrying one or more sequences of one or more instructions to processor 404 for execution. For example, the instructions may initially be borne on a magnetic disk of a remote computer. The remote computer can load the instructions into its dynamic memory and send the instructions over a telephone line using a modem. A modem local to computer system 400 can receive the data on the telephone line and use an infrared transmitter to convert the data to an infrared signal. An infrared detector coupled to bus 402 can receive the data carried in the infrared signal and place the data on bus 402. Bus 402 carries the data to main memory 406 from which processor 404 retrieves and executes the instructions. The instructions received by main memory 406 may optionally be stored on storage device 410 either before or after execution by processor 404.

[0059] Computer system 400 also includes a communication interface 418 coupled to bus 402. Communication interface 418 provides a two-way data communication coupling to a network link 420 that is connected to a local network 422. For example, communication interface 418 may be an integrated services digital network (ISDN) card or a modem to provide a data communication connection to a corresponding type of telephone line. As another example, communication interface 418 may be a local area network (LAN) card to provide a data communication connection to a compatible LAN. Wireless links may also be implemented. In any such implementation, communication interface 418 sends and receives electrical, electromagnetic, or optical signals that carry digital data streams representing various types of information.

[0060] Network link 420 typically provides data communication through one or more networks to other data devices. For example, network link 420 may provide a link to other access points (e.g. AP1 102, AP2 104 and/or AP3 106). As another example, computer system 400 can be used to implement TDOA Logic 116 in FIG. 1. Computer system 400 can obtain data from AP1 102, AP2, 104 and/or AP3 106 to perform TDOA calculations for client 118 through communication interface 418.

[0061] In view of the foregoing structural and functional features described above, methodologies in accordance with various aspects of the present invention will be better appreciated with reference to FIGS. 5 and 6. While, for purposes of simplicity of explanation, the methodologies of FIGS. 5 and 6 are shown and described as executing serially, it is to be understood and appreciated that the present invention is not limited by the illustrated order, as some aspects could, in accordance with the present invention, occur in different orders and/or concurrently with other aspects from that shown and described herein. Moreover, not all illustrated features may be required to implement a methodology in accordance with an aspect the present invention. Embodiments of the present invention are suitably adapted to implement the methodology in hardware, software, or a combination thereof.

[0062] Referring to FIG. 5, there is illustrated a methodology 500 in accordance with an aspect of the present invention. The description for methodology 500 will also refer to network 100 of FIG. 1 and timing diagram 200 of FIG. 2 for ease of illustration. It should be readily appreciated by one skilled in the art that aspects of methodology 500 are adaptable to any suitable wireless network and can be used to determine TDOA between signals received by more than two wireless devices. The description given herein refers to two devices for receiving a signal from client 118 for ease of illustration and should not be construed as limited to two devices.

[0063] Referring to FIG. 2, AP1 transmits a signal (AP1 TX) at t_{AP1 TX}. Client transmits a signal (Client TX) at t_{C TX}. AP1 receives the cleint signal at T_CTX+DAP1:C, where DAP1:C is the propagation delay from the client to AP1. AP2 receives the signal AP1 TX at t_{AP1 TX+DAP1:AP2}, where DAP1:AP2 is the known propagation delay between AP1 and AP2. AP2 receives the client signal Client TX at t_{C TX+DAP2:C}, where DAP2:C is the propagation delay between the client and AP2.

[0064] At 502, a first difference is computed. The first difference is the difference of signals perceived by AP1 102 (FIG. 2). The difference is between when AP1 102 receives the client packet and when the packet was transmitted by AP 1 102 (t_{C TX+DAP1:C}^CLK1−t_{AP1 TX}^CLK1 FIG. 2).

[0065] At 504, a second difference is computed. The second difference is derived from the time estimate from when the packet was received from AP1 102 (t_{AP1 TX+DAP1:AP2}^CLK2 FIG. 2) which is subtracted from the time when the client packet was received (t_CTX+DAP2:C^CLK2 FIG. 2).

[0066] At 506, the relative frequency difference between the clocks of AP1 102 and AP2 104 are computed. AP2 104 computes the relative frequency difference between f_CLK1 and f_CLK2 when it receives the packet from AP1 102, which can be estimated by Equation 14 described herein supra:

[0067] $\begin{array}{cc}\frac{f_{\mathrm{CLK}1}}{f_{\mathrm{CLK}2}}+\frac{f_{\mathrm{Err}}}{f_{\mathrm{nom}}}.& \end{array}$

[0068] At 508, a product is computed between the relative frequency difference (determined at 506) and the signals received by AP2 104 (determined at 504). At 510, the TDOA is calculated by the difference of signals for AP1 102 (determined at 502) is subtracted from the product determined at 508. It should be noted that method 500 does not depend on the order that the signals were received (e.g. whether the AP1 102 transmits first or client 118 transmits first). Methodology 500 assumes the propagation delay between AP1 102 and AP2 104 is known. The propagation delay can be a predetermined delay (e.g. measured previously) or can be dynamically calculated by an exchange of signals between AP1 102 and AP2 104. The estimated error for methodology 500 is determined using Equation 17 described herein supra.

[0069] Referring now to FIG. 6 there is illustrated a methodology 600 in accordance with an aspect of the present invention. As with methodology 500, methodology 600 assumes the propagation delay between AP1 102 and AP2 104 is known. Methodology 600 performs calculations based on two signals, wherein the propagation delay for one of the signals is known. A signal is transmitted by AP1 at TXAP1 and received by AP2 at RXAP2AP1. Another signal is transmitted by the client and received by AP1 at RXAP1C and by AP2 at RXAP2C.

[0070] When AP2 receives a signal from AP1, it can determine a transmit time relative to the clock for AP2. This is illustrated in step 602 where TX AP2 is calculated by subtracting the known propagation delay from RXAP2AP1.TX AP2 is used to calibrate the clock differences between AP1 and AP2. As shown at 604, CLKDIFFAP1AP2 (the time adjustment factor) is given by TXAP2−TXAP1. At 606, the TDOA is computed after adjusting RXAP2C. The TDOA algorithm uses RXAP1C and RXAP2C−CLKDIFFAP1AP2. In an alternate embodiment, the TDOA algorithm at 602 can use RXAP1C adjusted by the time difference factor and RXAP2C.

[0071] As with methodology 500 in FIG. 5, the order in which the signals from AP1 and the client are received is unimportant. The propagation delay between AP1, AP2 can be a predetermined fixed delay (e.g. a measured delay) or can be determined dynamically, such as by exchanging signals between AP1, AP2 to determine the propagation delay.

[0072] What has been described above includes exemplary implementations of the present invention. It is, of course, not possible to describe every conceivable combination of components or methodologies for purposes of describing the present invention, but one of ordinary skill in the art will recognize that many further combinations and permutations of the present invention are possible. Accordingly, the present invention is intended to embrace all such alterations, modifications and variations that fall within the spirit and scope of the appended claims interpreted in accordance with the breadth to which they are fairly, legally and equitably entitled.