US 6,169,958 AGrant

Ionospheric correction for single frequency GPS receivers using three satellites

Issue Date:2001-01-02

•9 Claims

•1 Drawing Sheets

Abstract

Systems and methods are disclosed for ionospheric correction in a system employing a single GPS frequency receiver for determining the geographic location of an object on the earth's surface. The receiver receives signals transmitted at the GPS L₁frequency from at least first, second and third GPS satellites, the first, second and third satellites having respective orbital positions relative to the receiver such that the total electron count (TEC) contribution to ionoshperic interference to signals transmitted by the respective satellites and received by the receiver is approximately the same. Respective measured distances of the three satellites to the receiver are determined based on the actual signal transmission times. True distances of the respective satellites are then calculated based on the assumption that the TEC contribution to the interference from each satellite is approximately the same.

Metadata

Assignee

VSIS, Inc.

Inventors

Vijay Nagasamy
Mohammad Usman
James Sun

Application Information

Application Number:US 09/289,552

Filing Date:1999-04-09

Priority Date:1999-04-09

Art Unit:7

Classifications

IPC:

G06F16500

Field of Search:

701207701213701214701215701216701221701224701225701226701300342357013423570234235703342357063423570834235709

Patent Drawings (1 sheets)

Description

Field of the Invention

[0002] The present invention pertains generally to the field of locating systems employing GPS satellites, and more specifically to methods and apparatus for ionospheric correction in a locating system employing a single frequency GPS receiver.

Background of the Invention

[0003] Satellite-based global positioning systems are well known. For example, U.S. Pat. No. 5,210,540 (the “'540 patent”), issued to Masumota discloses a global positioning system for locating a mobile object, such as an automobile, in a global geometrical region. As described therein, the system includes a Global Positioning System (GPS) receiver for receiving radio waves from a plurality of satellites and outputting either two or three dimensional position data indicative of the present position of the mobile object. For greater accuracy, an altimeter is also employed to detect the mobile object's relative altitude. U.S. Pat. No. 5,430,654 (the “'654 patent”), issued to Kyrtsos et al discloses apparatus and methods for determining the position of a vehicle at, or near the surface of the earth using a satellite-based navigation system, wherein “precise” position estimates are achieved by reducing the effective receiver noise. The respective '540 and '654 patents are incorporated herein by reference for all they disclose and teach.

[0004] In particular, the Global Positioning System is a satellite-based navigation system that was designed and paid for by the U.S. Department of Defense. The GPS consists of twenty-four satellites, which orbit the earth at known coordinates. The particulars of the GPS are described in section 3.4.2 of “Vehicle Location and Navigation Systems,” by Yilin Zhao, Artech House, Inc., 1997, which is fully incorporated herein by reference. As noted therein, the observation of at least four GPS satellites simultaneously will permit determination of three-dimensional coordinates of a receiver located on the earth's surface, as well as the time offset between the receiver and the respective satellites.

[0005] One problem encountered in global positioning systems is ionospheric interference. The ionosphere is a dispersive medium, which lies between seventy and one thousand kilometers above the earth's surface. The ionosphere effects a certain, frequency dependent propagation delay on signals transmitted from the respective GPS satellites. The ionosphere also effects GPS signal tracking by the receiver. Notably, the ionospheric delay of a transmitted GPS signal can cause an error of up to ten meters when calculating the exact geographic position on the earth's surface of the receiver.

[0006] As demonstrated below, delay from ionospheric interference can be almost completely corrected for by using multiple frequency observations, i.e., by transmitting and receiving signals at two different, known GPS frequencies, L₁and L₂, from a respective satellite. However, for security reasons, most GPS receivers do not receive the L₂frequency. Instead, these single (i.e., L₁) frequency receivers can employ a model to estimate and correct for transmission delay due to ionospheric interference.

[0007] For example, the Global Positioning System, Interface Control Document, ICD-GPS-200, Revision C, Initial Release, Oct. 10, 1993, provides a method for ionospheric correction based on an “approximate atmospheric” model, which is dependent on a “total electron content” (TEC) model. In accordance with this model, and using only the L₁signal transmission frequency, it can be shown that the error ΔS₁in the true satellite to receiver/user distance ρ is:

\begin{matrix} Δ S_{1} = F_{pp} 40.3 \frac{TEC}{L_{1}^{2}} & (1) \end{matrix}

[0008] where F_pp, is an “obliquity factor”:

\begin{matrix} F_{pp} = {([1 - {(\frac{R_{e} \cos φ}{R_{e} + H_{r}})}^{2}])}^{- \frac{1}{2}} & (2) \end{matrix}

[0009] where R_eis the radius of the earth, H_ris the height of the maximum electron density in the ionosphere from the earth's surface, and φ is the angle between the respective satellite and a plane tangent to the earth's surface at the receiver's position.

[0010] Notably, the true TEC value of the ionosphere is very difficult to model and is highly sensitive to variables, such as time of day, solar activity and relative incident angle of the satellite with respect to the sunlight trajectory (if any) toward the receiver location, etc. In particular, the TEC nominal value varies widely, within a range of between 10¹⁶to 10¹⁹. As a result, the above ionospheric correction model has been shown to adequately correct for no more than 50% of the ionospheric transmission delay.

[0011] As noted above, a dual frequency receiver can virtually eliminate ionospheric errors by computing the pseudo-range distance of the respective satellite on both the L₁and L₂frequencies. For purposes of illustration, a short derivation of such a dual frequency correction methodology is as follows:

[0012] Let ΔS₁and ΔS₂represent the error in the pseudo-range distances computed at frequencies L₁and L₂, respectively. Then:

\begin{matrix} Δ S_{1} = - F_{pp} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L_{1}} - ρ & (3) \end{matrix}

[0013] and

\begin{matrix} Δ S_{2} = - F_{pp} 40.3 \frac{TEC}{L_{2}^{2}} = λ_{L2} - ρ & (4) \end{matrix}

[0014] where λ_L1and λ_L2are the respective pseudo-range distances computed at frequencies L₁and L₂, respectively, and ρ is the true satellite to receiver distance. Dividing equation (3) by equation (4), results in:

\begin{matrix} \frac{Δ S_{1}}{Δ S_{2}} = \frac{L_{2}^{2}}{L_{1}^{2}} & (5) \end{matrix}

[0015] Subtracting equation (4) from equation (3), gives:

[0016] ΔS₁−ΔS₂=λ_L1−λ_L2 (6)

[0017] Substituting equation (5) into equation (6), and after some minor algebraic manipulation, provides:

\begin{matrix} Δ S_{1} = (\frac{L_{2}^{2}}{L_{1}^{2} - L_{1}^{2}}) (λ_{L_{2}} - λ_{L_{1}}) & (7) \end{matrix}

[0018] Importantly, all quantities in the above expression (7) are either known by the receiver, or can be measured, with the TEC value totally canceled out of the equation.

[0019] Of course, in a GPS-based locating system having only a single (i.e., L₁) frequency receiver, the above-described ionospheric correction model based on both the L₁and L₂transmission frequencies can not be employed.

Summary of the Invention

[0020] The present invention is directed to satellite-based systems and methods for correcting for ionospheric interference in a single-frequency receiver system for determining the geographic location of an object on the earth's surface without requiring ionospheric TEC modeling.

[0021] In accordance with one aspect of the invention, a method using satellite signal transmission for determining the geographic location of a receiver on the earth's surface, includes:

[0022] receiving a first signal transmitted at a known frequency from a first satellite having a known orbital position;

[0023] receiving a second signal transmitted at the same frequency as the first signal from a second satellite having a known orbital position;

[0024] receiving a third signal transmitted at the same frequency as the first signal from a third satellite having a known orbital position;

[0025] calculating measured distances λ¹, λ²and λ³of the respective first, second and third satellites from the receiver based at least in part on the transmission time of the third signal; and

[0026] calculating actual distances ρ¹, ρ²and ρ³of the respective first, second and third satellites from the receiver based on the measured distances λ¹, λ²and λ³, according to the relationships

\frac{F_{pp}^{1}}{F_{pp}^{2}} (λ^{2} - ρ^{2}) = λ^{1} - ρ^{1}, \frac{F_{pp}^{1}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{1} - ρ^{1}, and

\frac{F_{pp}^{2}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{2} - ρ^{2},

[0027] where F_pp¹, F_pp², and F_pp³are obliquity factors for the respective first, second and third satellites based on the respective angles φ₁, φ₂and φ₃they form with a plane tangent to the earth's surface at the geographic location of the receiver, with

F_{pp}^{1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

[0028] where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

[0029] As will be apparent to those skilled in the art, other and further aspects and advantages of the present invention will appear hereinafter.

Brief Description of the Drawings

[0030] Preferred embodiments of the present invention are illustrated by way of example, and not by way of limitation, in the accompanying drawing, in which:

[0031] FIG. 1 is a diagrammatic illustration of a GPS-based locating system for locating the geographic position of a single frequency GPS receiver on the earth's surface.

Detailed Description of the Preferred Embodiments

[0032] An exemplary GPS-based locating system is illustrated in FIG. 1. For ease in illustration of the present inventions disclosed and described herein, only three satellites and a single receiver are shown in the locating system. It is assumed that the timing offset between the respective satellites and receiver is either already known or otherwise calculated by use of an additional satellite (not shown).

[0033] In particular, a GPS receiver 50 is at a particular location on the surface 22 of the earth 24, e.g., within an object to be located. A first GPS satellite 56 is located above the earth's atmosphere at an angle φ₁between the satellite 56 and a plane 52 tangent to the earth's surface 22 at the geographic location of the GPS receiver 50. The actual distance between the first satellite 56 and the receiver 50 is designated as ρ^GPS1. Likewise, a second GPS satellite 57 is located above the earth's atmosphere at an angle φ₂formed between satellite 57 and plane 52, with the actual distance between the second satellite 57 and the receiver 50 designated as ρ^GPS2. A third GPS satellite 58 is located above the earth's atmosphere at an angle φ₃formed between satellite 58 and plane 52, with the actual distance between the third satellite 58 and the receiver 50 designated as ρ^GPS3. Although the GPS satellites 56, 57 and 58 broadcast at both the L_{1 and L}₂frequencies, in accordance with the inventions disclosed herein, it is assumed that only the L₁frequency is received by the receiver 50. Notably, the respective L₁signals travel through the earth's ionosphere 32 before they reach the GPS receiver 50.

[0034] More particularly, the ionospheric TEC is mainly a function of the time of day, solar activities, and relative position of the satellite relative to the respective sunlight trajectory (if any) towards the receiver 50. Thus, at a given instant while satellites 56, 57 and 58 are in orbit over the receiver 50, the ionospheric TEC value for each satellite is approximately the same. Based on signal transmission measurements taken from each satellite 56, 57 and 58 at the same time of day and, thus, with the same solar activity, the only variable is the particular relative location of each satellite with respect to the receiver 50.

[0035] In a preferred embodiment, the respective satellites 56, 57 and 58 are selected for making a location calculation by the receiver 50 based on their relative proximity in space to each other, thereby increasing the likelihood that the respective nominal TEC values for each satellite can be assumed to be approximately the same.

[0036] Accordingly, from above equation (1), the respective errors due to ionospheric interference in the true satellite to receiver distances ρ^GPS1, ρ^GPS2and ρ^GPS3are as follows:

\begin{matrix} Δ S_{1} = - F_{pp}^{GPS1} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L1}^{GPS1} - ρ^{GPS1} & (8) \end{matrix}

\begin{matrix} Δ S_{2} = - F_{pp}^{GPS2} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L_{1}}^{GPS2} - ρ^{GPS2} & (9) \end{matrix}

[0037] and

\begin{matrix} Δ S_{3} = - F_{pp}^{GPS3} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L_{1}}^{GPS3} - ρ^{GPS3} & (10) \end{matrix}

[0038] where λ_L1^GPS1, λ_L1^GPS2, λ₁^GPS3are the measured pseudo-range distances between the respective satellites 56, 57 and 58 and the receiver 50, and where F_pp^GPS1, F_pp^GPS2, and F_pp^GPS3are the obliquity factors for satellites 56, 57 and 58 based on the respective angles φ₁, φ₂and φ₃they form with the tangent plane 52. Following from equation (2):

\begin{matrix} F_{pp}^{GPS1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}} & (11) \end{matrix}

\begin{matrix} F_{pp}^{GPS2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}} & (12) \end{matrix}

[0039] and

\begin{matrix} F_{pp}^{GPS3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}} & (13) \end{matrix}

[0040] where R_eis the radius of the earth, H_ris the height of the maximum electron density in the ionosphere from the earth's surface, which are known to the receiver 50.

[0041] Considering the relationship between the first two satellites 56 and 57, and dividing equation (8) by equation (9):

\begin{matrix} \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS2}} = \frac{λ_{L1}^{GPS1} - ρ^{GPS1}}{λ_{L1}^{GPS2} - ρ^{GPS2}} & (14) \end{matrix}

[0042] Rearranging:

\begin{matrix} \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS2}} (λ_{L_{1}}^{GPS2} - ρ^{GPS2}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1} & (15) \end{matrix}

[0043] By the same approach, it follows that:

\begin{matrix} \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1} & (16) \end{matrix}

[0044] and

\begin{matrix} \frac{F_{pp}^{GPS2}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS2} - ρ^{GPS2} & (17) \end{matrix}

[0045] Thus, the three equations (15), (16) and (17) contain only three unknowns ρ^GPS1, ρ^GPS2and ρ^GPS3, whereby the true geographic location on the earth's surface 22 of the receiver 50 may be readily determined from the respective measured pseudo range distances of the three satellites 56, 57 and 58 at the receiver 50. Notably, the solution will not depend on the respective TEC values.

[0046] Alternatively, correction for ionospheric interference may be done by measuring signals from just two satellites, e.g., 56 and 57.

[0047] Because the exact orbit position coordinates of the satellites 56 and 57 (and, thus φ₁, and φ₂) are known to the receiver 50, the actual distance ρ^GPS1-2between the two satellites can be calculated as follows:

\begin{matrix} ρ^{GPS1 - 2} = \sqrt{{(x^{GPS1} - x^{GPS2})}^{2} + {(y^{GPS1} - y^{GPS2})}^{2} + {(z^{GPS1} - z^{GPS2})}^{2}} & (19) \end{matrix}

[0048] where (x^GPS1, y^GPS1, z^GPS1) are the coordinates in space of satellite 56, and (x^GPS2, y^GPS2, z^GPS2) are the coordinates in space of satellite 57.

[0049] Referring again to FIG. 2, and considering the triangle formed by ρ^GPS1-2, ρ^GPS1and ρ^GPS2:

[0050] (ρ^GPS1-2)²=(ρ^GPS1)²+(ρ^GPS2)²−2·ρ^GPS1·ρ^GPS2·cos α (20)

[0051] where α is the angle between ρ^GPS1and ρ^GPS2.

[0052] Letting β be the “measured angle” between λ_L1^GPS1and λ_L1^GPS2, since λ_L1^GPS1−ρ^GPS1<<λ_L1^GPS1and λ_L1^GPS2−ρ^GPS2<<λ_L1^GPS2, it is reasonable to assume cosα is approximately equal to cosβ. Therefore, equation (20) can be rewritten as

[0053] (ρ^GPS1-2)²=(ρ^GPS1)²+(ρ^GPS2)²−2·ρ^GPS1·ρ^GPS2·cos β (21)

[0054] Thus, the two equations, (21) and (15) contain only two unknowns, ρ^GPS1and ρ^GPS2, whereby the true geographic location on the earth's surface 22 of the receiver 50 may be readily determined from the respective measured pseudo range distances λ_L1^GPS1and λ_L1^GPS2of satellites 56 and 57, and, thus, the measured angle Ø formed by satellites 56 and 57 with respect to the receiver 50. Again, the solution will not depend on the respective TEC values.

[0055] As will be apparent to those skilled in the art, the above methods for ionospheric correction may also be applied to a GLONASS, or combined GPS/GLONASS receiver, so long as the respective transmission timing, frequency and orbit coordinates are known by the receiver for each respective satellite.

[0056] While preferred systems and methods for correcting for ionospheric interference in a single frequency GPS receiver system have been shown and described, as would be apparent to those skilled in the art, many modifications and applications are possible without departing from the inventive concepts herein. Thus, the scope of the disclosed invention is not to be restricted except in accordance with the appended claims.

Claims

What is claimed is:

1. A method using satellite signal transmission for determining the geographic location of a receiver on the earth's surface, comprising:

where F_pp¹, F_pp², and F_pp³are obliquity factors for the respective first, second and third satellites based on the respective angles φ₁, φ₂and φ₃they form with a plane tangent to the earth's surface at the geographic location of the receiver, with

F_{pp}^{1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

receiving a first signal transmitted at a known frequency from a first satellite having a known orbital position;

receiving a second signal transmitted at the same frequency as the first signal from a second satellite having a known orbital position;

receiving a third signal transmitted at the same frequency as the first signal from a third satellite having a known orbital position;

calculating measured distances λ¹, λ²and λ³of the respective first, second and third satellites from the receiver based at least in part on the transmission time of the respective first, second and third signals; and

calculating actual distances ρ¹, ρ²and ρ³of the respective first, second and third satellites from the receiver based on the measured distances λ¹, λ²and λ³, according to the relationships

\frac{F_{pp}^{1}}{F_{pp}^{2}} (λ^{2} - ρ^{2}) = λ^{1} - ρ^{1}, \frac{F_{pp}^{1}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{1} - ρ^{1}, and

\frac{F_{pp}^{2}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{2} - ρ^{2},

2. The method of claim 1, wherein the first second and third satellites are global positioning system (GPS) satellites.

3. The method of claim 2, wherein the first second and third signals are transmitted at the GPS L₁signal frequency.

4. The method of claim 1, wherein the first, second and third satellites have respective orbital positions relative to the receiver such that the total electron count (TEC) contribution to ionoshperic interference to the transmission of the first, second and third signals is approximately the same.

5. A system for determining the geographic location of objects on the earth's surface, comprising:

F_{pp}^{1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

a receiver configured to receive signals transmitted at a selected frequency from first, second and third satellites having known orbital positions, the receiver further configured to

calculate measured distances λ¹, λ²and λ³of the respective first, second and third satellites from the receiver based at least in part on the transmission time of the received signals, and

calculate actual distances ρ¹, ρ²and ρ³of the respective first, second and third satellites from the receiver based on the measured distances λ¹, λ²and λ³, according to the relationships

\frac{F_{pp}^{1}}{F_{pp}^{2}} (λ^{2} - ρ^{2}) = λ^{1} - ρ^{1}, \frac{F_{pp}^{1}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{1} - ρ^{1}, and

\frac{F_{pp}^{2}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{2} - ρ^{2},

6. The system of claim 5, wherein the first second and third satellites are global positioning system (GPS) satellites.

7. The system of claim 6, wherein the first second and third signals are transmitted at the GPS L₁signal frequency.

8. The system of claim 5, wherein the first, second and third satellites have respective orbital positions relative to the receiver such that the total electron count (TEC) contribution to ionoshperic interference to the transmission of the first, second and third signals is approximately the same.

9. A system for determining the geographic location of objects on the earth's surface, comprising:

where F_pp^GPS1, F_pp^GPS2, and F_pp^GPS3are obliquity factors for the respective first, second and third satellites based on the respective angles φ₁, φ₂and φ₃they form with a plane tangent to the earth's surface at the geographic location of the receiver, with

F_{pp}^{GPS1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{GPS2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{GPS3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

a single frequency global positioning system (GPS) receiver configured to receive signals transmitted at the GPS L₁frequency from first, second and third GPS satellites, the respective receiver and satellites having synchronized clocks, the first, second and third satellites having respective orbital positions relative to the receiver such that the total electron count (TEC) contribution to ionoshperic interference to signals transmitted by the respective satellites and received by the receiver is approximately the same, the receiver further configured to calculate measured distances λ¹, λ²and λ³of the respective first, second and third satellites from the receiver based at least in part on the transmission time of the received signals, and

calculate actual distances ρ^GPS1, ρ^GPS2and ρ^GPS3of the respective first, second and third satellites from the receiver based on the measured distances λ_L1^GPS1, λ_L1^GPS2and λ_L1^GPS3according to the relationships

\frac{F_{pp}^{GPS1}}{F_{pp}^{GPS2}} (λ_{L_{1}}^{GPS2} - ρ^{GPS2}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1}, \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1}, and

\frac{F_{pp}^{GPS2}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS2} - ρ^{GPS2},

Patent Citations (7)

Patent	Date	Inventor	Cited By
US5210540	1993-05-01	Masumoto	Applicant
US5430654	1995-07-01	Kyrtsos et al.	Applicant
US5506587	1996-04-01	Lans	Examiner
US5825326	1998-10-01	Semler et al.	Examiner
US5831545	1998-11-01	Murray et al.	Examiner
US5999878	1999-12-01	Hanson et al.	Examiner
US6032108	2000-02-01	Seiple et al.	Examiner

Non-Patent Literature (7)

Klobuchar, John A., Air Force Geophysics Laboratory, Ionospheric Effects on GPS, GSP World, Apr. 1991, pp. 48-51.Applicant
Wu, B., et al.,Mayflower Communications Co., Inc. and T.R. Jenkins, Wright Laboratory, WPAFB, Ionospheric Error Compensation for GPS Receivers Using Real-Time Ionospheric Model, pp. 575-583.Applicant
Beser, Jacques and Balendra, Anushia, 3S Navigation, Enhanced Ionospheric Delay Compensation Using GLONASS, pp. 1609-1617.Applicant
Lin, Lao-Sheng, School of Geomatic Engineering, The University of New South Wales, “A Novel Approach to Improving the Accuracy of Real-Time Ionospheric Delay Estimation Using GPS”, pp. 169-178.Applicant
Robert E. Daniell, Jr., et al., Computational Physics, Inc., “A New, Improved Ionospheric Correction Algorithm for Single Frequency GPS Receivers,” pp. 635-640.Applicant
Zhao, Yilin, “Vehicle Location and Navigation Systems,” Section, 3.4.2 Global Positioning System, pp. 63-75, ISBN 0-89006-861-5, (alk.paper); TE228.3Z45, 1997, Library of Congress Card No. 97-4200; Artech House, Inc., Norwood, MA 02062,.Applicant
Christie, Jock R.I., et al., Department of Aeronautics and Astronautics, Stanford University, The Effects of the Ionosphere and C/A Frequency on GPS Signal Shape: Considerations for GNSS-2, The Effects of the Onosphere and C/A Frequency on GPS Signal Shape: Considerations for GNSS-2, pp. 647-653.Applicant

Abstract

Description

Field of the Invention

Background of the Invention

\begin{matrix} Δ S_{1} = F_{pp} 40.3 \frac{TEC}{L_{1}^{2}} & (1) \end{matrix}

[0008] where F_pp, is an “obliquity factor”:

\begin{matrix} F_{pp} = {([1 - {(\frac{R_{e} \cos φ}{R_{e} + H_{r}})}^{2}])}^{- \frac{1}{2}} & (2) \end{matrix}

[0012] Let ΔS₁and ΔS₂represent the error in the pseudo-range distances computed at frequencies L₁and L₂, respectively. Then:

\begin{matrix} Δ S_{1} = - F_{pp} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L_{1}} - ρ & (3) \end{matrix}

[0013] and

\begin{matrix} Δ S_{2} = - F_{pp} 40.3 \frac{TEC}{L_{2}^{2}} = λ_{L2} - ρ & (4) \end{matrix}

\begin{matrix} \frac{Δ S_{1}}{Δ S_{2}} = \frac{L_{2}^{2}}{L_{1}^{2}} & (5) \end{matrix}

[0015] Subtracting equation (4) from equation (3), gives:

[0016] ΔS₁−ΔS₂=λ_L1−λ_L2 (6)

[0017] Substituting equation (5) into equation (6), and after some minor algebraic manipulation, provides:

\begin{matrix} Δ S_{1} = (\frac{L_{2}^{2}}{L_{1}^{2} - L_{1}^{2}}) (λ_{L_{2}} - λ_{L_{1}}) & (7) \end{matrix}

[0018] Importantly, all quantities in the above expression (7) are either known by the receiver, or can be measured, with the TEC value totally canceled out of the equation.

Summary of the Invention

[0021] In accordance with one aspect of the invention, a method using satellite signal transmission for determining the geographic location of a receiver on the earth's surface, includes:

[0022] receiving a first signal transmitted at a known frequency from a first satellite having a known orbital position;

[0023] receiving a second signal transmitted at the same frequency as the first signal from a second satellite having a known orbital position;

[0024] receiving a third signal transmitted at the same frequency as the first signal from a third satellite having a known orbital position;

[0025] calculating measured distances λ¹, λ²and λ³of the respective first, second and third satellites from the receiver based at least in part on the transmission time of the third signal; and

\frac{F_{pp}^{1}}{F_{pp}^{2}} (λ^{2} - ρ^{2}) = λ^{1} - ρ^{1}, \frac{F_{pp}^{1}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{1} - ρ^{1}, and

\frac{F_{pp}^{2}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{2} - ρ^{2},

F_{pp}^{1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

[0028] where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

[0029] As will be apparent to those skilled in the art, other and further aspects and advantages of the present invention will appear hereinafter.

Brief Description of the Drawings

[0030] Preferred embodiments of the present invention are illustrated by way of example, and not by way of limitation, in the accompanying drawing, in which:

[0031] FIG. 1 is a diagrammatic illustration of a GPS-based locating system for locating the geographic position of a single frequency GPS receiver on the earth's surface.

Detailed Description of the Preferred Embodiments

[0036] Accordingly, from above equation (1), the respective errors due to ionospheric interference in the true satellite to receiver distances ρ^GPS1, ρ^GPS2and ρ^GPS3are as follows:

\begin{matrix} Δ S_{1} = - F_{pp}^{GPS1} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L1}^{GPS1} - ρ^{GPS1} & (8) \end{matrix}

\begin{matrix} Δ S_{2} = - F_{pp}^{GPS2} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L_{1}}^{GPS2} - ρ^{GPS2} & (9) \end{matrix}

[0037] and

\begin{matrix} Δ S_{3} = - F_{pp}^{GPS3} 40.3 \frac{TEC}{L_{1}^{2}} = λ_{L_{1}}^{GPS3} - ρ^{GPS3} & (10) \end{matrix}

\begin{matrix} F_{pp}^{GPS1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}} & (11) \end{matrix}

\begin{matrix} F_{pp}^{GPS2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}} & (12) \end{matrix}

[0039] and

\begin{matrix} F_{pp}^{GPS3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}} & (13) \end{matrix}

[0040] where R_eis the radius of the earth, H_ris the height of the maximum electron density in the ionosphere from the earth's surface, which are known to the receiver 50.

[0041] Considering the relationship between the first two satellites 56 and 57, and dividing equation (8) by equation (9):

\begin{matrix} \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS2}} = \frac{λ_{L1}^{GPS1} - ρ^{GPS1}}{λ_{L1}^{GPS2} - ρ^{GPS2}} & (14) \end{matrix}

[0042] Rearranging:

\begin{matrix} \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS2}} (λ_{L_{1}}^{GPS2} - ρ^{GPS2}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1} & (15) \end{matrix}

[0043] By the same approach, it follows that:

\begin{matrix} \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1} & (16) \end{matrix}

[0044] and

\begin{matrix} \frac{F_{pp}^{GPS2}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS2} - ρ^{GPS2} & (17) \end{matrix}

[0046] Alternatively, correction for ionospheric interference may be done by measuring signals from just two satellites, e.g., 56 and 57.

\begin{matrix} ρ^{GPS1 - 2} = \sqrt{{(x^{GPS1} - x^{GPS2})}^{2} + {(y^{GPS1} - y^{GPS2})}^{2} + {(z^{GPS1} - z^{GPS2})}^{2}} & (19) \end{matrix}

[0048] where (x^GPS1, y^GPS1, z^GPS1) are the coordinates in space of satellite 56, and (x^GPS2, y^GPS2, z^GPS2) are the coordinates in space of satellite 57.

[0049] Referring again to FIG. 2, and considering the triangle formed by ρ^GPS1-2, ρ^GPS1and ρ^GPS2:

[0050] (ρ^GPS1-2)²=(ρ^GPS1)²+(ρ^GPS2)²−2·ρ^GPS1·ρ^GPS2·cos α (20)

[0051] where α is the angle between ρ^GPS1and ρ^GPS2.

[0053] (ρ^GPS1-2)²=(ρ^GPS1)²+(ρ^GPS2)²−2·ρ^GPS1·ρ^GPS2·cos β (21)

Claims

What is claimed is:

1. A method using satellite signal transmission for determining the geographic location of a receiver on the earth's surface, comprising:

F_{pp}^{1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

receiving a first signal transmitted at a known frequency from a first satellite having a known orbital position;

receiving a second signal transmitted at the same frequency as the first signal from a second satellite having a known orbital position;

receiving a third signal transmitted at the same frequency as the first signal from a third satellite having a known orbital position;

\frac{F_{pp}^{1}}{F_{pp}^{2}} (λ^{2} - ρ^{2}) = λ^{1} - ρ^{1}, \frac{F_{pp}^{1}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{1} - ρ^{1}, and

\frac{F_{pp}^{2}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{2} - ρ^{2},

2. The method of claim 1, wherein the first second and third satellites are global positioning system (GPS) satellites.

3. The method of claim 2, wherein the first second and third signals are transmitted at the GPS L₁signal frequency.

5. A system for determining the geographic location of objects on the earth's surface, comprising:

F_{pp}^{1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

a receiver configured to receive signals transmitted at a selected frequency from first, second and third satellites having known orbital positions, the receiver further configured to

calculate measured distances λ¹, λ²and λ³of the respective first, second and third satellites from the receiver based at least in part on the transmission time of the received signals, and

calculate actual distances ρ¹, ρ²and ρ³of the respective first, second and third satellites from the receiver based on the measured distances λ¹, λ²and λ³, according to the relationships

\frac{F_{pp}^{1}}{F_{pp}^{2}} (λ^{2} - ρ^{2}) = λ^{1} - ρ^{1}, \frac{F_{pp}^{1}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{1} - ρ^{1}, and

\frac{F_{pp}^{2}}{F_{pp}^{3}} (λ^{3} - ρ^{3}) = λ^{2} - ρ^{2},

6. The system of claim 5, wherein the first second and third satellites are global positioning system (GPS) satellites.

7. The system of claim 6, wherein the first second and third signals are transmitted at the GPS L₁signal frequency.

9. A system for determining the geographic location of objects on the earth's surface, comprising:

F_{pp}^{GPS1} = - {[1 - {(\frac{R_{e} \cos φ_{1}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, F_{pp}^{GPS2} = - {[1 - {(\frac{R_{e} \cos φ_{2}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}}, and

F_{pp}^{GPS3} = {[1 - {(\frac{R_{e} \cos φ_{3}}{R_{e} + H_{r}})}^{2}]}^{- \frac{1}{2}},

where R_eis the radius of the earth and H_ris the height of maximum electron density in the ionosphere surrounding the earth's surface.

\frac{F_{pp}^{GPS1}}{F_{pp}^{GPS2}} (λ_{L_{1}}^{GPS2} - ρ^{GPS2}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1}, \frac{F_{pp}^{GPS1}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS1} - ρ^{GPS1}, and

\frac{F_{pp}^{GPS2}}{F_{pp}^{GPS3}} (λ_{L_{1}}^{GPS3} - ρ^{GPS3}) = λ_{L_{1}}^{GPS2} - ρ^{GPS2},

Patent

Date

Inventor

Cited By

US5210540

1993-05-01

Masumoto

Applicant

US5430654

1995-07-01

Kyrtsos et al.

Applicant

US5506587

1996-04-01

Lans

Examiner

US5825326

1998-10-01

Semler et al.

Examiner

US5831545

1998-11-01

Murray et al.

Examiner

US5999878

1999-12-01

Hanson et al.

Examiner

US6032108

2000-02-01

Seiple et al.

Examiner

Non-Patent Literature (7)

Klobuchar, John A., Air Force Geophysics Laboratory, Ionospheric Effects on GPS, GSP World, Apr. 1991, pp. 48-51.Applicant

Wu, B., et al.,Mayflower Communications Co., Inc. and T.R. Jenkins, Wright Laboratory, WPAFB, Ionospheric Error Compensation for GPS Receivers Using Real-Time Ionospheric Model, pp. 575-583.Applicant

Beser, Jacques and Balendra, Anushia, 3S Navigation, Enhanced Ionospheric Delay Compensation Using GLONASS, pp. 1609-1617.Applicant

Lin, Lao-Sheng, School of Geomatic Engineering, The University of New South Wales, “A Novel Approach to Improving the Accuracy of Real-Time Ionospheric Delay Estimation Using GPS”, pp. 169-178.Applicant

Robert E. Daniell, Jr., et al., Computational Physics, Inc., “A New, Improved Ionospheric Correction Algorithm for Single Frequency GPS Receivers,” pp. 635-640.Applicant

Zhao, Yilin, “Vehicle Location and Navigation Systems,” Section, 3.4.2 Global Positioning System, pp. 63-75, ISBN 0-89006-861-5, (alk.paper); TE228.3Z45, 1997, Library of Congress Card No. 97-4200; Artech House, Inc., Norwood, MA 02062,.Applicant

Christie, Jock R.I., et al., Department of Aeronautics and Astronautics, Stanford University, The Effects of the Ionosphere and C/A Frequency on GPS Signal Shape: Considerations for GNSS-2, The Effects of the Onosphere and C/A Frequency on GPS Signal Shape: Considerations for GNSS-2, pp. 647-653.Applicant