Calibration of positioning systems
Summary by NHIP
Mobile Station Calibration
The method calibrates mobile station locations by estimating bias errors without added hardware. It derives position approximations using an n-dimensional bias vector calculated from derivatives of measurement functions and refines positions via a specific covariance equation.
Claim Score by NHIP
Abstract
Mobile station positioning methods are calibrated using an error bias estimation for refining the determination of the mobile station location which does not require any added hardware to the telecommunications system. The position of the mobile station is calculated assuming no bias errors and a first order approximation of the mobile station position is derived as a function of the bias error. The bias error is then estimated and used to refine the previously calculated mobile station position.

Term
Term ended
Expired 23 August 2022, 4.1 years ago.
- Priority
- Filed
- Granted
- Expired
- Today
16 claims: 2 independent, 14 dependent
- 1A method of calibrating a mobile station location system in a cellular telecommunications network, the method comprising the steps of:determining a position (z k0 ) of a mobile station in the telecommunications network assuming no bias error;deriving, concurrently, an approximation of the position (z k ) of a k th mobile station as a function of an n-dimensional vector (b) of unknown biases according to a calibration bias estimation equation b = b 0 - { ∑ k = 1 N H k T P k - 1 H k - ∑ k = 1 N ( H k T P k - 1 G k ( G k P k - 1 G k ) - 1 G k T P k - 1 H k ) } - 1 ( ∑ k = 1 N H k T P k - 1 ( Y k - r ( z k0 , b 0 ) ) ) wherein b is an updated estimate of b 0 , b 0 is the apriori value of a bias vector, Y k is a generic m-dimensional measurement related to the k th mobile station, z k is the position of the k th mobile station represented by the column vector (xy), H k is a derivative of r(z,b) with respect to bias and G k is a derivative of the function of r(z,b) with respect to position;and refining the position (z k0 ) of the mobile station, using the updated estimate of b, according to the equation Y k =r ( z k ,b )+ e k where k=1, . . . ,N, and e k is an additive noise term with covariance P k .
- 9Broadest claimClaim Score 12, narrow(NHIP)A system for calibrating a mobile station location system in a cellular telecommunications network, the system comprising:means for determining a position (z k0 ) of a mobile station in the telecommunications network assuming no bias error;means for deriving, concurrently, an approximation of the position (z k ) of a k th mobile station as a function of an n-dimensional vector (b) of unknown biases according to a calibration bias estimation equation b = b 0 - { ∑ k = 1 N H k T P k - 1 H k - ∑ k = 1 N ( H k T P k - 1 G k ( G k P k - 1 G k ) - 1 G k T P k - 1 H k ) } - 1 ( ∑ k = 1 N H k T P k - 1 ( Y k - r ( z k0 , b 0 ) ) ) wherein b is an updated estimate of b 0 ,b 0 is the apriori value of a bias vector, Y k is a generic m-dimensional measurement related to the k th mobile station, z k is the position of the k th mobile station represented by the column vector (xy), H k is a derivative of r(z,b) with respect to bias and G k is a derivative of the function of r(z,b) with respect to position;and means for refining the position (z k0 ) of the mobile station using the updated estimate of b, according to the equation Y k =r(z k ,b )+ e k where k=1, . . . ,N, and e k is an additive noise term with covariance P k .
Independent claims2
34 paragraphs in 5 sections, as filed
0001This application claims priority under 35 U.S.C. §§ 119 and/or 365 to U.S. Provisional Application Ser. No. 60/201,424 filed in the United States of America on May 3, 2000; the entire content of which is hereby incorporated by reference.
FIELD OF INVENTION
0002The present invention relates to the positioning of mobile stations within a cellular telecommunications network, and more particularly, to the calibration of positioning systems within a cellular telecommunications network.
BACKGROUND
0003Legal as well as commercial demands are driving a great interest in developing methods for positioning mobile stations (MS's). While mobile positioning may be achieved through the use of external systems, like Global Positioning Satellite (GPS) systems, mobile positioning can also be achieved by means internal to the cellular system. Several cellular positioning methods have been proposed, each of which relies upon measurements such as time of arrival (TOA), round trip delay, or angle of arrival of communication signals. Cellular positioning methods can be divided into uplink or downlink methods, i.e., whether the base station (BS) performs measurements on mobile station originating signals (uplink methods) or the mobile station performs measurements on base station originating signals (downlink methods). One example of a proposed downlink method is the Enhanced-Observed Time Difference (E-OTD) method. The E-OTD method is based on measuring the time-of-arrival (TOA) of bursts transmitted from a BS on its broadcast control channel (BCCH) carrier.
0004All location or positioning methods require knowledge of a number of parameters, some time-varying and others of more or less constant nature. These parameters can be TOA measurements, angle of arrival measurements, signal strength measurement or other parameters known to be used in mobile station positioning. For example, the E-OTD method requires that the location of the base stations and the real time difference (RTD) between transmissions from pairs of base stations are precisely known. Whenever these conditions are not fulfilled the location accuracy can be degraded.
0005As discussed above, the Enhanced-Observed Time Difference is one example of a downlink positioning method. The E-OTD method is based on three quantities: the observed time difference (OTD); the real time difference (RTD); and the geometric time difference (GTD). These three quantities are related by the basic E-OTD equation: <br /><i>OTD=RTD+GTD</i> (1)<br /> where OTD is the time difference between the reception of signals from two base stations, RTD is the timing or synchronization difference between two base stations, and GTD is the difference of propagation path length between two base stations. <figref idref="DRAWINGS">FIG. 1</figref> depicts these quantities which are used in the E-OTD method.
0006Referring to <figref idref="DRAWINGS">FIG. 1</figref>, if bursts from base station BS<b>1</b> (e.g. serving base station) are received by a mobile station MS<b>1</b> at time t<b>1</b> and bursts from a neighbor base station BS<b>2</b> are received by the mobile station at time t<b>2</b>, then the observed time difference is calculated as <br /><i>OTD=t</i><b>2</b>−<i>t</i><b>1</b>. (2)<br /> Likewise, if BS<b>1</b> transmits its bursts at a time t<b>3</b> and BS<b>2</b> transmits its bursts at time t<b>4</b>, then the real time difference is calculated as <br /><i>RTD=t</i><b>4</b>−<i>t</i><b>3</b>. (3)<br /> Finally, if d<b>1</b> is the propagation path length from BS<b>1</b> to MS<b>1</b>, and d<b>2</b> is the propagation path length between BS<b>2</b> and MS<b>1</b>, then the geometric time difference is calculated as <br /><i>GTD=</i>(<i>d</i><b>2</b>−<i>d</i><b>1</b>)/<i>c,</i> (4)<br /> where d<b>1</b> and d<b>2</b> are the distance of the propagation path from BS<b>1</b> and BS<b>2</b>, respectively, and c is the speed of the radio waves, usually taken as the speed of light in vacuum.
0007As illustrated in equation 5 below, if the OTD and RTD is known, the GTD can be determined in order to derive the location (i.e., position) information. <br /><i>GTD=OTD−RTD</i> (5)<br /> Since the GTD is based on the difference in distance between a mobile station and two base stations, a known GTD defines a hyperbola upon which the Mobile station can be located. The position of the mobile station is given by the intersection of the hyperbolas defined by the GTDs, as shown in FIG. <b>2</b>. Accordingly, at least two GTD values are required in order to locate the mobile station.
0008Once a mobile station has measured the OTD values, in order to obtain the desired GTD values, the RTDs must be known. In addition, in order to obtain the location of the mobile station from the GTDs, the location of the base stations must be known. There are several methods for obtaining the RTD and the base station location.
0009One method of obtaining the RTD values is to measure the real time differences using Location Measurement Units (LMUs). In its simplest form, an E-OTD LMU consists of an E-OTD capable mobile station placed at a fixed known location in the network. The E-OTD LMU measures the OTDs between pairs of base stations, and assuming that the locations of the LMUs and the base stations are known (i.e., the GTDs are known), RTD values can be calculated using the following equation: <br /><i>RTD=OTD−GTD</i> (6)<br /> Depending on the network and the location of the LMUs (i.e., co-located at a base station site or at a location between base stations), the ratio between base stations and LMUs in the network is estimated to be 3:1 (i.e., for every 3 base stations in the network one LMU is required). Drawbacks of the LMU approach are that it requires new hardware to be deployed in the network, and depending on the stability of the base station clocks, the LMU's may have to report RTDs quite often, which increases the network load.
0010An alternative method of determining the RTD values is to synchronize the network. This is the case, for example, in IS95 where GPS is used as a timing source. Also for GSM and TDMA, synchronization is being discussed and is indeed a requirement for deployment of a technology called compact EDGE. In a synchronized network, the base stations (or LMUs) don't need to report the RTDs since they are assumed to be constant. However, the synchronization is usually obtained in the digital baseband domain and not the air interface. It will be recognized that synchronization in the air interface is important for accurately determining location as a number of additional delay sources are present in the air interface which must be accounted for, including transmit filters, combiners and transmission lines. One solution for obtaining air interface synchronization is to use an LMU which measures the transmission on the air interface from its serving cell and relates the base station time to absolute time. The alternative of synchronizing the network partially overcomes the drawbacks associated with the LMU method, however; the synchronization needs to take the air interface into account. The required level of synchronization on the air is 5 μs for IS95, which is enough for communication, but far too lax for location services (5 μs corresponds to a traveled distance of 1500 m). Factory calibration of all delay sources, including transmit filters, combiners and transmission lines is one possible solution to accounting for the air interface, but this may be prohibitive for cost reasons.
0011The base station and LMU locations are normally not known precisely by the operator today, or are not accurate to the level required by location systems. This is due to the fact that for communication purposes exact locations are not that critical. To deploy location systems, therefore, it is necessary to measure precisely the BS and LMU locations, possibly by maintenance personnel equipped with differential GPS receivers.
0012The problems with base station location calibration is that maintenance personnel need to visit each site and be equipped with differential GPS receivers. This may be costly and furthermore it may not be possible to obtain GPS coverage everywhere, e.g., inside buildings. In some regions in the world, differential GPS might not even be available. In these regions, the operator will have to rely on the accuracy provided by the standard service mode of GPS, which gives location accuracy of around 100 m, 95% of the time.
0013Accordingly, the invention provides a method for calibration of a positioning system within a telecommunications network which overcomes the drawbacks associated with known calibration methods. The method of the invention has the advantage that no new hardware is required, as it is entirely based on measurements which are made by multiple mobile stations.
SUMMARY OF THE INVENTION
0014As a solution to the above described problems, a method is disclosed that provides an estimation of the bias errors in the parameters that are relevant for mobile positioning.
0015According to embodiments of the invention, a method and system for determining mobile positioning is provided which accounts for bias errors in the relevant positioning parameters. This is achieved performing a plurality of measurements associated with a plurality of mobile stations and creating calibration parameters based on estimated mobile positions and the plurality of measurements. The estimated position of the mobiles is then updated using the estimated calibration parameters.
0016According to another embodiment of the invention, a method for estimating bias errors in relevant positioning parameters is provided.
BRIEF DESCRIPTION OF THE DRAWINGS
0017These and other objects, features and advantages of the present invention will be readily apparent to one skilled in the art from the following written description, read in conjunction with the drawings, in which:
0018<figref idref="DRAWINGS">FIG. 1</figref> illustrates the quantities which define OTD, RTD and GTD;
0019<figref idref="DRAWINGS">FIG. 2</figref> illustrates the principles of the E-OTD location method;
0020<figref idref="DRAWINGS">FIG. 3</figref> illustrates the steps of the bias estimation procedure;
0021<figref idref="DRAWINGS">FIGS. 4A and 4B</figref> illustrate resulting location accuracy for different noise variances, using timing calibration in accordance with the present invention; and
0022<figref idref="DRAWINGS">FIG. 5</figref> illustrates an exemplary embodiment of the present invention.
DETAILED DESCRIPTION OF THE INVENTION
0023In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular techniques, steps, circuits, circuit components etc. in order to provide a thorough understanding of the present invention. However, it will be apparent to one of ordinary skill in the art that the present invention may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known methods, devices, and circuits are omitted so as not to obscure the description of the present invention with unnecessary detail.
0024These and other aspects of the invention will now be described in greater detail in connection with a number of exemplary embodiments. To facilitate an understanding of the invention, many aspects of the invention are described in terms of sequences of actions to be performed by elements of a communication system. It will be recognized that in each of the embodiments, the various actions could be performed by specialized circuits, by program instructions being executed by one or more processors, or by a combination of both. Moreover, the invention can additionally be considered to be embodied entirely within any form of computer readable storage medium having stored therein an appropriate set of computer instructions that would cause a processor to carry out the techniques described herein. Thus, the various aspects of the invention may be embodied in many different forms, and all such forms are contemplated to be within the scope of the invention.
0025The calculation of bias error estimates according to an embodiment of the invention is formulated and solved in a generic framework, such that the procedure is applicable to various positioning methods. In the generic case, the measurement equation for determining the location of a mobile station can be defined as <br /><i>Y</i><sub>k</sub><i>=r</i>(<i>z</i><sub>k</sub><i>,b</i>)+<i>e</i><sub>k</sub><i>k=</i>1<i>, . . . , N</i> (7)<br /> where Y<sub>k </sub>is a generic m-dimensional measurement related to the kth mobile station, b is an n-dimensional vector of unknown biases, z<sub>k </sub>is the position of the kth mobile station represented by the column vector (x<sub>k </sub>y<sub>k</sub>)<sup>T</sup>, and e<sub>k </sub>is an additive noise term with covariance P<sub>k</sub>. Equation (7) can be written in compact form as <br /><i>Y=f</i>(<i>z,b</i>)+<i>e</i> (8)<br /> with <br /><i>Y=[Y</i><sub>1</sub><sup>T</sup><i>, . . . , Y</i><sub>N</sub><sup>1</sup>]<sup>1</sup><i>, f</i>(<i>z,b</i>)=[<i>r</i>(<i>z</i><sub>1</sub><i>,b</i>)<sup>T</sup><i>, . . . ,r</i>(<i>z</i><sub>N</sub><i>,b</i>)<sup>T</sup>]<sup>T</sup><i>, e=[e</i><sub>1</sub><sup>T</sup><i>, . . . ,e</i><sub>N</sub><sup>T</sup>]<sup>T</sup>.<br /> The goal is to estimate both the location z<sub>k </sub>and the bias vector b. This is achieved by first ensuring that z<sub>k </sub>and b can be uniquely estimated.
0026Generally it is known that the number of equations must be larger than the number of unknowns. This is a necessary requirement, but as shown later, not sufficient. Accordingly, the present invention requires that mN>2N+n−1. For uniqueness, it is also required that <br /><i>f</i>(<i>z</i><sub>1</sub><i>,b</i><sub>1</sub>)=<i>f</i>(<i>z</i><sub>2</sub><i>,b</i><sub>2</sub>)=><i>z</i><sub>1</sub><i>−z</i><sub>2</sub><i>,b</i><sub>1</sub><i>−b</i><sub>2</sub>. (9)<br /> As an example, consider the case of range measurements with uncertain locations. To simplify notation, the coordinates are represented as complex numbers, i.e., z=x+iy. In this case, r becomes <br /><i>r</i>(<i>z</i><sub>k</sub><i>,b</i><sub>s</sub>)=<i>abs</i>(<i>z</i><sub>k</sub><i>−b</i><sub>s</sub>) (10)<br /> where abs( ) represents the elementwise absolute values of the quantity in brackets and b<sub>s </sub>denotes the unknown sensor locations. Now b<sub>s2 </sub>is selected such that <br /><i>b</i><sub>s2</sub>=exp(<i>i</i>Φ)(<i>b</i><sub>s</sub><i>+c</i>) (11)<br /> where c is a complex arbitrary constant and Φ is an arbitrary rotation angle. This process amounts to a translation and rotation of the sensor positions. Then <maths id="MATH-US-00001" num="00001"><math overflow="scroll"><mtable><mtr><mtd><mtable><mtr><mtd><mrow><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><msub><mi>b</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mi>abs</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>-</mo><msub><mi>b</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mrow><mi>abs</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>-</mo><mrow><msup><mrow><mi>exp</mi><mo></mo><mrow><mo>(</mo><mrow><mrow><mo>-</mo><mi>ⅈ</mi></mrow><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>Φ</mi></mrow><mo>)</mo></mrow></mrow><mo>*</mo></msup><mo></mo><msub><mi>b</mi><mi>s2</mi></msub></mrow><mo>+</mo><mi>c</mi></mrow><mo>)</mo></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mrow><mi>abs</mi><mo>(</mo><mrow><mrow><mi>exp</mi><mo></mo><mrow><mo>(</mo><mrow><mrow><mo>-</mo><mi>ⅈ</mi></mrow><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>Φ</mi></mrow><mo>)</mo></mrow></mrow><mo></mo><mrow><mo>(</mo><mrow><mrow><mrow><mi>exp</mi><mo></mo><mrow><mo>(</mo><mrow><mi>ⅈ</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>Φ</mi></mrow><mo>)</mo></mrow></mrow><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>+</mo><mi>c</mi></mrow><mo>)</mo></mrow></mrow><mo>-</mo><msub><mi>b</mi><mi>s2</mi></msub></mrow><mo>)</mo></mrow></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><mrow><mrow><mi>exp</mi><mo></mo><mrow><mo>(</mo><mrow><mi>ⅈ</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>Φ</mi></mrow><mo>)</mo></mrow></mrow><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>+</mo><mi>c</mi></mrow><mo>)</mo></mrow></mrow><mo>,</mo><msub><mi>b</mi><mi>s2</mi></msub></mrow><mo>)</mo></mrow></mrow></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Hence the range measurement corresponding to the mobile station at z<sub>k </sub>and sensors at b<sub>s </sub>could equally well have been originated from the mobile station located at exp(iΦ)(z<sub>k</sub>+c) and sensors located at b<sub>s2</sub>. It can easily be verified that this holds for all measurements, i.e., f(z,b<sub>s1</sub>)=f(exp(iΦ)(z+c),b<sub>s2</sub>). Accordingly, the full bias vector can not be estimated. However, if, for example, one sensor location is accurately known and the orientation relative to one other sensor is known, then c and Φ can be selected such that b<sub>s2 </sub>contains only n−3 unknown parameters.
0027Now turning to the estimation of z<sub>1</sub>, . . . z<sub>N</sub>, b, this can be achieved in a least squares sense by minimizing the loss function defined in equation (13) using a Newton-Raphson approach. <maths id="MATH-US-00002" num="00002"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mi>V</mi><mo></mo><mrow><mo>(</mo><mrow><mi>z</mi><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow><mo>=</mo><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><mo></mo><msup><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>k</mi></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow><mi>T</mi></msup><mo></mo><mrow><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>k</mi></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> With θ=(z<sub>1</sub><sup>T</sup>, . . . , z<sub>N</sub><sup>T</sup>,b<sup>T</sup>)<sup>T</sup>, one iteration step looks like <br />θ=θ<sub>0</sub><i>−V″</i><sup>−1</sup>(θ<sub>0</sub>)<i>V</i>′(θ<sub>0</sub>) (14)<br /> In addition, with the first derivative of the loss function with respect to position defined as <br /><i>dV/dz</i><sub>k</sub><i>=−G</i><sub>k</sub><sup>T</sup><i>P</i><sub>k</sub><sup>−1</sup>(<i>Y</i><sub>k</sub><i>−r</i>(<i>z</i><sub>k</sub><i>,b</i>)) (15)<br /> where G<sub>k </sub>is the derivative of the function r(z,b) with respect to position, defined as <br /><i>G</i><sub>k</sub><i>=dr</i>(<i>z,b</i>)/<i>dz</i><sup>T</sup> (16)<br /> and, with the first derivative of the loss function with respect to the bias vector defined as <maths id="MATH-US-00003" num="00003"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><mo>ⅆ</mo><mi>V</mi></mrow><mo>/</mo><mrow><mo>ⅆ</mo><mi>b</mi></mrow></mrow><mo>=</mo><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><mrow><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>k</mi></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> where H<sub>k </sub>is the derivative of the function r(z,b) with respect to bias, defined as <br /><i>H</i><sub>k</sub><i>=dr</i>(<i>z</i><sub>k</sub><i>,b</i>)/<i>db</i><sup>T</sup> (18)<br /> The derivative vector of the loss function of equation (13) can therefore be written as shown in equation (19). <maths id="MATH-US-00004" num="00004"><math overflow="scroll"><mtable><mtr><mtd><mrow><msup><mi>r</mi><mi>′</mi></msup><mo>=</mo><mrow><mrow><mo>(</mo><mi>Θ</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>[</mo><mrow><mrow><msup><mrow><mo>(</mo><mrow><msubsup><mi>G</mi><mn>1</mn><mi>T</mi></msubsup><mo></mo><mrow><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mn>1</mn></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mn>1</mn></msub><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow><mi>T</mi></msup><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>…</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><msup><mrow><mo>(</mo><mrow><msubsup><mi>G</mi><mi>N</mi><mi>T</mi></msubsup><mo></mo><mrow><msubsup><mi>P</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>N</mi></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>N</mi></msub><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow><mi>T</mi></msup></mrow><mo>,</mo><mrow><mo>(</mo><mrow><munder><mo>∑</mo><mo>=</mo></munder><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><mrow><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>k</mi></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> The second derivatives are <br /><i>dV/dz</i><sub>k</sub><i>dz</i><sub>k</sub><sup>T</sup><i>=G</i><sub>k</sub><sup>T</sup><i>P</i><sub>k</sub><sup>−1</sup><i>G</i><sub>k</sub> (20)<br /><i>dV/dz</i><sub>k</sub><i>dz</i><sub>j</sub><sup>T</sup>=0 (21)
0028<br /><i>dV/dz</i><sub>k</sub><i>db</i><sup>T</sup><i>=G</i><sub>k</sub><sup>T</sup><i>P</i><sub>k</sub><sup>−1</sup><i>H</i><sub>k</sub> (22) <maths id="MATH-US-00005" num="00005"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><mrow><mo>ⅆ</mo><mi>V</mi></mrow><mo>/</mo><mrow><mo>ⅆ</mo><mi>b</mi></mrow></mrow><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mrow><mo>ⅆ</mo><msup><mi>b</mi><mi>T</mi></msup></mrow></mrow><mo>=</mo><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mi>k</mi></msub></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>23</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths> <i>dV/dbdz</i><sub>k</sub><sup>T</sup><i>=H</i><sub>k</sub><sup>T</sup><i>P</i><sub>k</sub><sup>−1</sup><i>G</i><sub>k</sub> (<b>24</b>) <br /> The second derivative matrix can thus be written as <maths id="MATH-US-00006" num="00006"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><msup><mi>V</mi><mi>″</mi></msup><mo></mo><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mo></mo><mtable><mtr><mtd><mrow><msubsup><mi>G</mi><mn>1</mn><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mn>1</mn></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><msubsup><mi>G</mi><mn>1</mn><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><msubsup><mi>G</mi><mi>N</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mi>N</mi></msub></mrow></mtd><mtd><mrow><msubsup><mi>G</mi><mi>N</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mi>N</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>H</mi><mn>1</mn><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mn>1</mn></msub></mrow></mtd><mtd><mi>⋯</mi></mtd><mtd><mrow><msubsup><mi>H</mi><mi>N</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mi>N</mi></msub></mrow></mtd><mtd><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mi>k</mi></msub></mrow></mrow></mtd></mtr></mtable><mo></mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>25</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Iteration Steps: <br /> The second derivative matrix can be partitioned as <maths id="MATH-US-00007" num="00007"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><msup><mi>V</mi><mi>″</mi></msup><mo>=</mo><mrow><mo></mo><mtable><mtr><mtd><mi>A</mi></mtd><mtd><mi>B</mi></mtd></mtr><mtr><mtd><mi>D</mi></mtd><mtd><mi>C</mi></mtd></mtr></mtable><mo></mo></mrow></mrow><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mi>where</mi><mo></mo><mstyle><mtext></mtext></mstyle><mo></mo><mi>A</mi><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><msubsup><mi>G</mi><mn>1</mn><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mn>1</mn></msub></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><msubsup><mi>G</mi><mi>N</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mi>N</mi></msub></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo>,</mo><mstyle><mtext> </mtext></mstyle><mo></mo><mrow><mi>B</mi><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><msubsup><mi>G</mi><mn>1</mn><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mi>⋯</mi></mtd></mtr><mtr><mtd><mrow><msubsup><mi>G</mi><mi>N</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mi>N</mi></msub></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow><mo>,</mo><mstyle><mtext></mtext></mstyle><mo></mo><mrow><mi>D</mi><mo>=</mo><mrow><mo>⌈</mo><mtable><mtr><mtd><mrow><msubsup><mi>H</mi><mn>1</mn><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mn>1</mn></msub></mrow></mtd><mtd><mi>⋯</mi></mtd><mtd><mrow><msubsup><mi>H</mi><mi>N</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>N</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mi>N</mi></msub></mrow></mtd></mtr></mtable><mo>⌉</mo></mrow></mrow><mo>,</mo><mstyle><mtext> </mtext></mstyle><mo></mo><mrow><mrow><mi>and</mi><mo></mo><mstyle><mtext> </mtext></mstyle><mo></mo><mi>C</mi></mrow><mo>=</mo><mrow><mrow><mo></mo><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mi>k</mi></msub></mrow></mrow><mo></mo></mrow><mo>.</mo></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>26</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> By using results on inverses of the partitioned matrices, it can be shown that <maths id="MATH-US-00008" num="00008"><math overflow="scroll"><mtable><mtr><mtd><mrow><msup><mrow><mo>⌊</mo><mtable><mtr><mtd><mi>A</mi></mtd><mtd><mi>B</mi></mtd></mtr><mtr><mtd><mi>D</mi></mtd><mtd><mi>C</mi></mtd></mtr></mtable><mo>⌋</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>=</mo><mrow><mrow><mo>⌊</mo><mtable><mtr><mtd><mi>U</mi></mtd><mtd><mi>V</mi></mtd></mtr><mtr><mtd><mi>W</mi></mtd><mtd><mi>X</mi></mtd></mtr></mtable><mo>⌋</mo></mrow><mo>.</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>27</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Similarly, the first derivative matrix can be partitioned into <br />V′=[E<sup>T</sup>F<sup>T</sup>]<sup>T</sup> (28)<br /> with <br /><i>E</i>=[(<i>G</i><sub>1</sub><sup>T</sup><i>P</i><sub>1</sub><sup>−1</sup>(<i>Y</i><sub>1</sub><i>−r</i>(<i>z</i><sub>1</sub><i>,b</i>))<sup>T</sup>, . . . , (<i>G</i><sub>N</sub><sup>T</sup><i>P</i><sub>N</sub><sup>−1</sup>(<i>Y</i><sub>N</sub><i>−r</i>(<i>z</i><sub>N</sub><i>,b</i>))<sup>T</sup>]<sup>T</sup> (29)<br /> and <maths id="MATH-US-00009" num="00009"><math overflow="scroll"><mtable><mtr><mtd><mrow><mi>F</mi><mo>=</mo><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><mrow><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>k</mi></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k</mi></msub><mo>,</mo><mi>b</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>30</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> This means that the update for b will be as shown in equation (31). <br /><i>b=b</i><sub>0</sub><i>−WE−XF</i> (31)<br /> As shown in <figref idref="DRAWINGS">FIG. 3</figref>, the procedure for estimating the bias error begins as step <b>301</b> by determining z<sub>k </sub>assuming no bias, i.e, b=0. This is achieved by selecting z<sub>k </sub>to minimize V(z<sub>k</sub>,b<sub>0</sub>), where b<sub>0 </sub>is the apriori value of the bias vector. Call these values z<sub>k0</sub>. In this case the minimizing z<sub>k0 </sub>satisfies V′(z<sub>k0</sub>,b<sub>0</sub>)=0. Then, at step <b>303</b>, an approximation of z<sub>k </sub>as a function of b is derived as discussed above. But since the vector E consists of V′(z<sub>ko</sub>,b<sub>0</sub>), k=1, . . . , N it follows that E=0, hence equation (31) reduces to that shown in equation (32). <br /><i>b=b</i><sub>0</sub><i>−XF</i> (32)<br /> Then at step <b>305</b>, the approximation derived in step <b>303</b> is solved for b. It is well known in the literature that an explicit expression for X can be defined as shown in equation (33), provided the inverses exist. <br /><i>X</i>=(<i>C−DA</i><sup>−1</sup><i>B</i>)<sup>−1</sup> (33)<br /> The final expression for equation (32) then becomes <maths id="MATH-US-00010" num="00010"><math overflow="scroll"><mtable><mtr><mtd><mtable><mtr><mtd><mrow><mi>b</mi><mo>=</mo><mi /><mo></mo><mrow><msub><mi>b</mi><mn>0</mn></msub><mo>-</mo><mrow><mo>{</mo><mrow><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mi>k</mi></msub></mrow></mrow><mo>-</mo></mrow></mrow></mrow></mrow></mtd></mtr><mtr><mtd><msup><mrow><mi /><mo></mo><mrow><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><mo>(</mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msup><mrow><msub><mi>G</mi><mi>k</mi></msub><mo></mo><mrow><mo>(</mo><mrow><msub><mi>G</mi><mi>k</mi></msub><mo></mo><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>G</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo></mo><msubsup><mi>G</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><msub><mi>H</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>}</mo></mrow></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mtd></mtr><mtr><mtd><mrow><mi /><mo></mo><mrow><mo>(</mo><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo></mo><mrow><msubsup><mi>H</mi><mi>k</mi><mi>T</mi></msubsup><mo></mo><mrow><msubsup><mi>P</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo></mo><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>k</mi></msub><mo>-</mo><mrow><mi>r</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>z</mi><mi>k0</mi></msub><mo>,</mo><msub><mi>b</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mo>(</mo><mn>34</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Finally, in step <b>307</b>, z<sub>k </sub>is determined using the updated estimate of b derived in step <b>305</b>.
0029This procedure for estimating bias errors can be utilized in known or proposed positioning methods, for example, a E-OTD positioning method with timing calibration error. In this E-OTD positioning method the mobile station performs pseudorange measurements so that the basic measurement equation is <br />ρ<sub>k</sub>=ρ(<i>z</i><sub>k</sub>)+<i>c+d</i><sub>k</sub><i>+v</i><sub>k</sub> (35)<br /> with <br />ρ(<i>z</i>)=√{square root over ((<i>x−x</i><sub>s</sub>)<sup>2</sup>+(<i>y−y</i><sub>s</sub>)<sup>2</sup>)}{square root over ((<i>x−x</i><sub>s</sub>)<sup>2</sup>+(<i>y−y</i><sub>s</sub>)<sup>2</sup>)} (36)<br /> where c is the unknown timing calibration bias, d is related to the unknown transmit time, and v is a noise term with covariance N, and x<sub>s </sub>and y<sub>s </sub>of equation (36) are the base station coordinate vectors. The goal is to estimate c and z (and d) simultaneously. The unknown scalar d is usually eliminated by defining the differences of the pseudorange measurements, T, as <maths id="MATH-US-00011" num="00011"><math overflow="scroll"><mtable><mtr><mtd><mrow><mi>T</mi><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd></mtr><mtr><mtd><mi>⋯</mi></mtd><mtd><mi>⋯</mi></mtd><mtd><mi>⋯</mi></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>37</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> with the remaining parameters defined as b=Tc, r(z,b)=Tρ(z)+b, c=Tv, cov(e)=P=TNT<sup>T</sup>, then the resulting equation is exactly (7). Note that a consequence of the reformulation above is that b and not c will be possible to estimate. It is straightforward to show that in this special case, G of equation (16), becomes <br /> G=TJ (38) <br /> where J is the derivative of ρ, defined as <maths id="MATH-US-00012" num="00012"><math overflow="scroll"><mtable><mtr><mtd><mrow><mi>J</mi><mo>=</mo><msup><mrow><mo>[</mo><mtable><mtr><mtd><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub></mrow><msub><mi>ρ</mi><mn>1</mn></msub></mfrac></mtd><mtd><mi>⋯</mi></mtd><mtd><mfrac><mrow><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>n</mi></msub></mrow><msub><mi>ρ</mi><mi>n</mi></msub></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub></mrow><msub><mi>ρ</mi><mn>1</mn></msub></mfrac></mtd><mtd><mi>⋯</mi></mtd><mtd><mfrac><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>n</mi></msub></mrow><msub><mi>ρ</mi><mi>n</mi></msub></mfrac></mtd></mtr></mtable><mo>]</mo></mrow><mi>T</mi></msup></mrow></mtd><mtd><mrow><mo>(</mo><mn>39</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Furthermore, H of equation (18), in this case, the derivative of ρ with respect to b becomes simply the identity matrix I. <br />H=I (40)
0030To illustrate the above procedure, a small network consisting of 5 base stations was simulated. The base stations were arranged in a circle, and 300 mobile stations were randomly placed in the vicinity of the base stations. Pseudoranges ρ were estimated by computing the MS-BS distances, adding a random calibration bias which was the same for all mobile stations and an individual transmission time bias term. White noise was added to the pseudoranges ρ. The locations were estimated assuming zero bias and in the second step the calibration bias was estimated. The mobile station locations were then recalculated using the estimated bias terms. The resulting location accuracy is illustrated in <figref idref="DRAWINGS">FIGS. 4A and 4B</figref> for different values of the noise variance and the magnitude of the bias errors. The results without timing calibration (curves <b>401</b> and <b>404</b>), with timing bias elimination (curves <b>402</b> and <b>405</b>) and, for comparison, from a bias free run (curves <b>403</b> and <b>406</b>) are shown. As illustrated by curves <b>402</b> an <b>405</b>, the bias is almost completely eliminated in both cases using the method in accordance with the invention.
0031This procedure can also be applied to E-OTD methods where there is uncertainty in the BTS locations. In this E-OTD positioning method (i.e., with BTS location uncertainty), the measurement equation is identical to the previous case except <br /><i>r</i>(<i>z</i><sub>k</sub><i>,b</i>)=<i>Tρ</i>(<i>z</i><sub>k</sub><i>,b</i>) (41)<br /> where <maths id="MATH-US-00013" num="00013"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mi>ρ</mi><mo></mo><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><mrow><mi>x</mi><mo>-</mo><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>s</mi></msub><mo>+</mo><msub><mi>b</mi><mi>x</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><mrow><mo>(</mo><mrow><msub><mi>y</mi><mi>s</mi></msub><mo>+</mo><msub><mi>b</mi><mrow><mi>y</mi><mo></mo><mstyle><mtext> </mtext></mstyle></mrow></msub></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mtd><mtd><mrow><mo>(</mo><mn>42</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> and b=[b<sub>x</sub>b<sub>y</sub>]<sup>T</sup>. As discussed above, it is not possible to estimate z<sub>k </sub>and b uniquely. To restrict the search space it is assumed that the location of the first two base stations are well known. In this case, the matrix G is also identical to equation (38), however; the matrix H is defined as <br />H=TK (43)<br /> where K is the derivative of ρ(z) with respect to b<sub>x </sub>and b<sub>y</sub>, defined as <maths id="MATH-US-00014" num="00014"><math overflow="scroll"><mtable><mtr><mtd><mrow><mi>K</mi><mo>=</mo><mrow><mo>⌊</mo><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><mrow><mo>-</mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub></mrow><mo>)</mo></mrow></mrow><msub><mi>ρ</mi><mn>3</mn></msub></mfrac></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd><mtd><mfrac><mrow><mo>-</mo><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>3</mn></msub></mrow><mo>)</mo></mrow></mrow><msub><mi>ρ</mi><mn>3</mn></msub></mfrac></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>⋯</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mfrac><mrow><mo>-</mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>-</mo><msub><mi>x</mi><mi>n</mi></msub></mrow><mo>)</mo></mrow></mrow><msub><mi>ρ</mi><mi>n</mi></msub></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mfrac><mrow><mo>-</mo><mrow><mo>(</mo><mrow><mi>y</mi><mo>-</mo><msub><mi>y</mi><mi>n</mi></msub></mrow><mo>)</mo></mrow></mrow><msub><mi>ρ</mi><mi>n</mi></msub></mfrac></mtd></mtr></mtable><mo>⌋</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>44</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths>
0032<figref idref="DRAWINGS">FIG. 5</figref> illustrates the principle of operation of the invention in connection with timing calibration bias. With reference to <figref idref="DRAWINGS">FIG. 5</figref>, the Location Services (LCS) application <b>503</b> sends a request for a location to a mobile location center, which is illustrated in <figref idref="DRAWINGS">FIG. 5</figref> as position calculation function (PCF) <b>507</b>. The PCF <b>507</b> upon receiving the request sends an E-OTD measurement request message to a mobile station, illustrated in <figref idref="DRAWINGS">FIG. 5</figref> as E-OTD MS#k <b>509</b>. The Mobile station performs the necessary E-OTD measurements (e.g., pseudoranges) and sends the results back to PCF <b>507</b> in an E-OTD Measurement Response message. The PCF <b>507</b> computes the position of the Mobile station using, for example, equation (5). The resulting estimate is sent to the location application requesting the location. In parallel, the location estimate, along with the previous estimate of the calibration bias and mobile station reported pseudoranges, are sent with a calibration data message to the Timing calibration function <b>501</b>. The Timing Calibration Function <b>501</b> collects received location estimates from PCF <b>507</b> until sufficient data is available to perform calibration bias estimation according to equation (34). The resulting estimate is forwarded to PCF <b>507</b> in a Timing Calibration report. Periodically, or depending on the synchronization architecture, PCF <b>507</b> also receives RTD reports from Location Measurement Units (LMUs) or Base Stations (BS) <b>505</b>. The PCF <b>507</b> uses initially the RTD reports in the position calculation, and later includes the updated calibration reports from the Timing calibration function <b>501</b>. It will be recognized that <figref idref="DRAWINGS">FIG. 5</figref> represents functional elements used in accordance with exemplary embodiments of the present invention and that these functional elements may be contained in one node in a radio communication system, or may be distributed across several nodes.
0033The invention has been described with reference to particular embodiments. However, it will be readily apparent to those skilled in the art that it is possible to embody the invention in specific forms other than those of the preferred embodiments described above. This may be done without departing from the spirit of the invention.
0034Thus, the preferred embodiment is merely illustrative and should not be considered restrictive in any way. The scope of the invention is given by the appended claims, rather than the preceding description, and all variations and equivalents which fall within the range of the claims are intended to be embraced therein.
Contents5
22 sheets
Sheet 1 Sheet 2 Sheet 3 Sheet 4 Sheet 5 Sheet 6 Sheet 7 Sheet 8 Sheet 9 Sheet 10 Sheet 11 Sheet 12 Sheet 13 Sheet 14 Sheet 15 Sheet 16 Sheet 17 Sheet 18 Sheet 19 Sheet 20 Sheet 21 Sheet 22
Every citation, both ways
| Document | Relation | Office | Cited during |
|---|---|---|---|
| US2003125045A1 | Cited by | United States of America | Pre-grant |
| US10750309B2 | Cited by | United States of America | Applicant |
| US2004082344A1 | Cited by | United States of America | Pre-grant |
| US9942705B1 | Cited by | United States of America | Applicant |
| US10952180B2 | Cited by | United States of America | Applicant |
| US9854402B1 | Cited by | United States of America | Applicant |
| US9702721B2 | Cited by | United States of America | Applicant |
| US8532567B2 | Cited by | United States of America | Applicant |
| US10299071B2 | Cited by | United States of America | Applicant |
| US9893917B2 | Cited by | United States of America | Applicant |
| US2011244891A1 | Cited by | United States of America | Pre-grant |
| US7646337B2 | Cited by | United States of America | Applicant |
| US8437251B2 | Cited by | United States of America | Applicant |
| US8514771B2 | Cited by | United States of America | Applicant |
| US9338767B2 | Cited by | United States of America | Applicant |
| US9544860B2 | Cited by | United States of America | Applicant |
| US10750310B2 | Cited by | United States of America | Applicant |
| US9891055B2 | Cited by | United States of America | Applicant |
| US10165059B2 | Cited by | United States of America | Applicant |
| US7528772B2 | Cited by | United States of America | Search report |
| US8503938B2 | Cited by | United States of America | Applicant |
| US9615204B1 | Cited by | United States of America | Applicant |
| US9654921B1 | Cited by | United States of America | Applicant |
| US11356799B2 | Cited by | United States of America | Applicant |
| US10805866B2 | Cited by | United States of America | Applicant |
| US8811348B2 | Cited by | United States of America | Applicant |
| US11221221B2 | Cited by | United States of America | Applicant |
| US10841739B2 | Cited by | United States of America | Applicant |
| US9451491B2 | Cited by | United States of America | Applicant |
| US10368199B2 | Cited by | United States of America | Applicant |
| US10856099B2 | Cited by | United States of America | Applicant |
| US9702709B2 | Cited by | United States of America | Applicant |
| US8787943B2 | Cited by | United States of America | Search report |
| US9661519B2 | Cited by | United States of America | Applicant |
| US8514692B2 | Cited by | United States of America | Applicant |
| US8600297B2 | Cited by | United States of America | Applicant |
| US2010093377A1 | Cited by | United States of America | Pre-grant |
| US9883360B1 | Cited by | United States of America | Applicant |
| US8478228B2 | Cited by | United States of America | Applicant |
| US10149092B1 | Cited by | United States of America | Applicant |
| US2005046613A1 | Cited by | United States of America | Pre-grant |
| US11778415B2 | Cited by | United States of America | Applicant |
| US10750311B2 | Cited by | United States of America | Applicant |
| US11419092B2 | Cited by | United States of America | Applicant |
| US9979776B2 | Cited by | United States of America | Applicant |
| US10341809B2 | Cited by | United States of America | Applicant |
| US9967704B1 | Cited by | United States of America | Applicant |
| US7605747B1 | Cited by | United States of America | Search report |
| US9736618B1 | Cited by | United States of America | Applicant |
| US9148795B2 | Cited by | United States of America | Applicant |
| US7626534B1 | Cited by | United States of America | Search report |
| US10341808B2 | Cited by | United States of America | Applicant |
| US10159006B2 | Cited by | United States of America | Applicant |
| US11665665B2 | Cited by | United States of America | Applicant |
| US9125093B2 | Cited by | United States of America | Applicant |
| US10313826B2 | Cited by | United States of America | Applicant |
| US10645693B2 | Cited by | United States of America | Applicant |
| US9854394B1 | Cited by | United States of America | Applicant |
| US10791414B2 | Cited by | United States of America | Applicant |
| US9749790B1 | Cited by | United States of America | Applicant |
| US10200811B1 | Cited by | United States of America | Applicant |
| US10508921B2 | Cited by | United States of America | Applicant |
| US2010283683A1 | Cited by | United States of America | Pre-grant |
| US10412703B2 | Cited by | United States of America | Applicant |
| US10064158B2 | Cited by | United States of America | Applicant |
| US8587479B2 | Cited by | United States of America | Applicant |
| US10959120B2 | Cited by | United States of America | Applicant |
| US7280831B2 | Cited by | United States of America | Search report |
| RU2699831C1 | Cited by | Russian Federation | Search report |
| US9137771B2 | Cited by | United States of America | Applicant |
| US9955298B1 | Cited by | United States of America | Applicant |
| US4494119A | Cites | United States of America | Search report |
| US5974329A | Cites | United States of America | Search report |
| US6014102A | Cites | United States of America | Search report |
| US6026304A | Cites | United States of America | Search report |
| US6047192A | Cites | United States of America | Search report |
| US6108555A | Cites | United States of America | Search report |
| US6128501A | Cites | United States of America | Search report |
| US6154657A | Cites | United States of America | Search report |
| US6246884B1 | Cites | United States of America | Search report |
| US6249253B1 | Cites | United States of America | Search report |
| US6334059B1 | Cites | United States of America | Search report |
| US6381463B1 | Cites | United States of America | Search report |
| US6424837B1 | Cites | United States of America | Search report |
| US6489923B1 | Cites | United States of America | Search report |
| WO9801768A1 | Cites | World Intellectual Property Organization (WIPO) | Applicant |
| WO9955018A1 | Cites | World Intellectual Property Organization (WIPO) | Applicant |
6 priority claims, no other members on record
Priority claims6
| Document | Office | Kind | Date |
|---|---|---|---|
| 20142400 | United States of America | P | |
| 20142400 | United States of America | P | |
| 84735701 | United States of America | A | |
| 60201424 | – | – | – |
| US20000201424P | – | – | – |
| US20010847357 | – | – | – |
54 transactions on the USPTO file
Allowed after 2 non-final rejections, 1 final rejection and 1 RCE.
- Non-final rejections
- 2
- Final rejections
- 1
- RCEs
- 1
- Appeals
- 0
Over time
Point at a mark for the transactionTransactions
| Event | |
|---|---|
| Expire Patent | |
| Post Issue Communication - Certificate of Correction | |
| Recordation of Patent Grant Mailed | |
| Patent Issue Date Used in PTA CalculationAllowed | |
| Issue Notification MailedAllowed | |
| Receipt into Pubs | |
| Dispatch to FDC | |
| Dispatch to FDC | |
| Application Is Considered Ready for Issue | |
| Receipt into Pubs | |
| Issue Fee Payment Verified | |
| Issue Fee Payment Received | |
| Workflow - File Sent to Contractor | |
| Mail Notice of AllowanceAllowed | |
| Notice of Allowance Data Verification CompletedAllowed | |
| Case Docketed to Examiner in GAU | |
| Date Forwarded to Examiner | |
| Date Forwarded to Examiner | |
| Disposal for a RCE / CPA / R129 | |
| Workflow - Request for RCE - Finish | |
| Workflow - Request for RCE - Finish | |
| Workflow - Request for RCE - Begin | |
| Miscellaneous Incoming Letter | |
| Request for Continued Examination (RCE) | |
| Mail Final Rejection (PTOL - 326)Final rejection | |
| Final RejectionFinal rejection | |
| IFW TSS Processing by Tech Center Complete | |
| Case Docketed to Examiner in GAU | |
| Date Forwarded to Examiner | |
| Response after Non-Final Action | |
| Workflow incoming amendment IFW | |
| Mail Non-Final RejectionNon-final rejection | |
| Non-Final RejectionNon-final rejection | |
| Date Forwarded to Examiner | |
| Response after Non-Final Action | |
| Request for Extension of Time - Granted | |
| Mail Non-Final RejectionNon-final rejection | |
| Non-Final RejectionNon-final rejection | |
| Case Docketed to Examiner in GAU | |
| Case Docketed to Examiner in GAU | |
| Case Docketed to Examiner in GAU | |
| Case Docketed to Examiner in GAU | |
| Application Dispatched from OIPE | |
| Application Is Now Complete | |
| Application Is Now Complete | |
| Oath or Declaration Filed (Including Supplemental) | |
| Reference capture on IDS | |
| Information Disclosure Statement (IDS) Filed | |
| Information Disclosure Statement (IDS) Filed | |
| New or Additional Drawing Filed | |
| Notice Mailed--Application Incomplete--Filing Date Assigned | |
| Correspondence Address Change | |
| IFW Scan & PACR Auto Security Review | |
| Initial Exam Team nn |
8 legal events, as the office reported them to INPADOC
Over the term
Point at a mark for the eventEvents
| Event | Code | |
|---|---|---|
| Lapsed due to failure to pay maintenance feeLapsedFP | FP | |
| Lapse for failure to pay maintenance feesLapsedPATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.)LAPS | LAPS | |
| Information on status: patent discontinuationPATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362STCH | STCH | |
| Maintenance fee reminder mailedREMI | REMI | |
| Fee paymentFPAY | FPAY | |
| Fee paymentFPAY | FPAY | |
| Certificate of correctionCC | CC | |
| AssignmentAS | AS |
Numbers
- Publication
- 06957072
- Publication, DOCDB
- 6957072
- Publication, EPODOC
- US6957072
- Application
- 9847357
- Application, DOCDB
- 84735701
- Application, EPODOC
- US20010847357
Titles
- English
- Calibration of positioning systems
Patent term adjustment
- A delay
- +508 daysthe office missed an examination deadline
- Applicant delay
- −31 days
- Net adjustment
- 477 days
Classification
- CPC, 4
- H04W64/00
- G01S5/021
- G01S5/0242
- G01S5/0244
- IPC, 4
- G01S5 02
- G01S19 09
- G01S19 46
- H04W64 00
- USPC, 4
- 455456100
- 342450000
- 342464000
- 455456500