Method for precision dynamic differential positioning

• US 4,812,991 A
• Filed: 05/01/1986
• Issued: 03/14/1989
• Est. Priority Date: 05/01/1986
• Status: Expired due to Term

• Current Assignee: Mesc Electronic Systems, Inc., Hughes Electronics Corporation (Rtx Corporation)
• Original Assignee: Magnavox Government and Industrial Electronics Company
• Inventors: Hatch, Ronald R.
• Primary Examiner(s): Lall, Parshotam S.
• Assistant Examiner(s): Trans, V. N.
• Application Number:
  US06/858,206
• Time in Patent Office: 1,048 Days
• Field of Search: 364/444, 364/449, 364/451, 364/452, 364/458, 364/460, 364/571, 342/394, 342/356, 342/357, 342/429, 342/442, 342/450, 342/451
• US Class Current: 701/225
• CPC Class Codes: G01S 19/44 Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

Abstract

A method and apparatus for rapidly and accurately determining the position coordinates of a remote, movable receiver relative to a fixed reference receiver. The technique utilizes successive code measurements and carrier phase measurements of both the L₁ and L₂ carrier signals broadcast from four or more orbiting GPS satellites. Code measurements based on a weighted average of the individual L₁ and L₂ code measurements in each satellite/receiver link are adjusted in accordance with the corresponding carrier phase measurement for an L₁ -L₂ carrier difference signal and are further smoothed over time. This yields a rapid determination of the remote receiver'"'"'s position coordinates, with progressively increasing accuracy. After merely about two to three minutes of processing, a wide laning carrier phase procedure can be implemented to yield a position determination that is accurate to within about one centimeter.

Claims

1. A method for determining the coordinates of a remote receiver relative to a reference receiver having known coordinates, the method using the modulated L₁ and L₂ carrier signals broadcast from four or more GPS satellites, the method comprising steps of:

measuring the range to each of the four or more satellites from both the reference receiver and the remote receiver, at each of a succession of time points, the step of measuring including steps of

detecting the L₁ and L₂ codes, to produce a code measurement for each satellite/receiver pair, at each time point, the code measurement being based on a predetermined linear combination of the detected L₁ and L₂ codes, and

detecting the L₁ and L₂ carrier phases, to produce a carrier phase measurement for each satellite/receiver pair, at each time point, the carrier phase measurement being based on a predetermined linear combination of the detected L₁ and L₂ carrier phases,

wherein the predetermined linear combinations of the detected L₁ and L₂ codes and carrier phases are selected such that the ionosphere affects the resulting code and carrier phase measurements substantially equally;

smoothing the successive code measurements for each satellite/receiver pair in accordance with the corresponding carrier phase measurements for the same time point and in accordance with the code measurements and carrier phase measurements for previous time points, to produce a smoothed code measurement for each satellite/receiver pair, at each time point;

comparing the successive smoothed code measurements for the reference receiver with theoretical range values based on the known coordinates of the reference receiver and known orbits of the four or more satellites, to produce an error value for each smoothed code measurement for the reference receiver;

determining the errors in the internal clocks of each of the four or more satellites, based on the error values produced in the preceding step of comparing;

adjusting the successive smoothed code measurements for the remote receiver, to eliminate the effects of the satellite clock errors determined in the preceding step of determining and thereby produce a succession of corrected range measurements for each satellite; and

determining the particular coordinates of the remote receiver that minimize errors in the successive corrected range measurements.

2. A method as defined in claim 1, wherein the first step of detecting that is part of the step of measuring includes steps of:

detecting the individual L₁ and L₂ codes for each satellite/receiver pair; and

computing a weighted average of the detected individual L₁ and L₂, codes to produce the code measurement for each satellite/receiver pair, at each time point, each of the code measurements having less noise then the individual L₁ and L₂ codes on which it is based.

3. A method as defined in claim 1, wherein:

the second step of detecting that is part of the step of measuring includes steps of differencing the L₁ and L₂ carrier signals for each satellite/receiver pair, at each time point, and detecting the carrier phase of each resulting L₁ -L₂ difference signal, this difference being the carrier phase measurement; and

the step of smoothing includes steps of

computing an expected value for each code measurement, based on the corresponding smoothed code measurement for the previous time point and on the difference between the corresponding carrier phase measurements for the same time point and the previous time point, and

computing a weighted average of each code measurement and the corresponding expected code measurement, to produce the smoothed code measurement for each satellite/receiver pair, at each time point.

4. A method as defined in claim 3, wherein the final step of determining includes steps of:

providing theoretical range values from the remote receiver to the four or more satellites, at each time point, based on an estimate of the remote receiver'"'"'s position coordinates and the known satellite orbits;

differencing the theoretical range values, determined in the step of providing, and the corresponding corrected range measurements, produced in the step of adjusting, to produce a set of error values for each time point; and

adjusting the estimate of the remote receiver'"'"'s position coordinates at each time point, to minimize the corresponding set of error values.

5. A method as defined in claim 4, wherein:

the method steps are performed recursively; and

the estimate of the remote receiver'"'"'s position coordinates at each time point, used in the step of providing, is based on the adjusted estimate provided in the step of adjusting for the previous time point.

6. A method as defined in claim 3, wherein the final step of determining includes steps of:

differencing the L₁ -L₂ carrier phase measurements for the remote receiver, produced in the preceding step of measuring, and the corresponding corrected range measurements, produced in the step of adjusting, to produce a set of error values for each time point; and

determining the particular position coordinates for the remote receiver that minimize the set of error values at each time point.

7. A method as defined in claim 3, wherein the final step of determining includes steps of:

differencing the L₁ carrier phase measurements for the remote receiver, produced in the step of measuring, and the corresponding corrected range measurements, produced in the step of adjusting, to produce a set of error values for each time point; and

determining the particular position coordinates for the remote receive that minimize the set of error values at each time point.

8. A method as defined in claim 3, wherein the final step of determining includes steps of:

computing a weighted average of the L₁ and L₂ carrier phase measurements, produced in the step of measuring for the remote receiver, at each time point;

differencing the weighted averages, computed in the preceding step of computing, and the corresponding corrected range measurements, produced in the step of adjusting, to produce a set of error values for each time point; and

determining the particular position coordinates for the remote receiver that minimize the set of error values at each time point.

9. A method as defined in claim 1, wherein the method is performed recursively, in real time.

10. A method as defined in claim 1, wherein:

the step of comparing includes a step of differencing the successive smoothed code measurements for the reference receiver and the theoretical range values, based on the known coordinates of the reference site and known orbits of the four or more satellites, to produce the error value for each satellite, at each time point; and

the first step of determining includes a step of defining the satellite clock errors to be equal to the fractional portions of the corresponding error values produced in the preceeding step of differencing.