Systems and methods for real time kinematic satellite positioning

1. A method for Real Time Kinematic satellite positioning comprising:

at a mobile receiver, receiving a first navigation satellite carrier signal from a first navigation satellite, receiving a second navigation satellite carrier signal from a second navigation satellite, receiving a third navigation satellite carrier signal from a third navigation satellite, and receiving a fourth navigation satellite carrier signal from a fourth navigation satellite;

at the mobile receiver, receiving first, second, third, and fourth phase correction signals from a reference station;

at the mobile receiver, calculating a set of integer phase ambiguities from double-differenced measurements of pseudo-range and phase; wherein calculating the set of integer phase ambiguities comprises generating a first set of integer phase ambiguity hypotheses using a Kalman filter and performing hypothesis testing on the first set of integer phase ambiguity hypotheses; wherein generating the first set of integer phase ambiguity hypotheses comprises generating the first set of integer phase ambiguity hypotheses using means and covariances generated by a Bierman-Thornton filter and at least one of a LAMBDA and an MLAMBDA algorithm; wherein generating the first set of integer phase ambiguity hypotheses comprises generating the first set of integer phase ambiguity hypotheses from a measurement equation that relates phase change and pseudo-range to integer ambiguity without inclusion of a baseline vector; wherein generating the first set of integer phase ambiguity hypotheses comprises generating the first set of integer phase ambiguity hypotheses without a dynamic transition model; and

at the mobile receiver, calculating a position of the mobile receiver from the set of integer phase ambiguities and the double-differenced measurements of pseudo-range and phase.

2. The method of claim 1 , wherein receiving a navigation satellite carrier signal comprises receiving a first GPS signal on an L 1 frequency.

3. The method of claim 2 , wherein receiving a navigation satellite carrier signal further comprises receiving a second GPS signal on an L 2 frequency and combining the first and second GPS signals to create a beat signal.

4. The method of claim 1 , further comprising receiving a set of additional phase correction signals from a second reference station; wherein calculating the position of the mobile receiver comprises calculating the position of the mobile receiver at least in part on the set of additional phase correction signals.

5. The method of claim 1 , wherein receiving the first, second, third, and fourth phase correction signals comprises receiving the first, second, third, and fourth phase correction signals from a UHF radio.

6. The method of claim 1 , wherein generating the first set of integer phase ambiguity hypotheses further comprises performing a decorrelating reparameterization of a hypothesis search space.

7. The method of claim 6 , wherein performing hypothesis testing comprises performing hypothesis testing using a Bayesian update algorithm.

8. The method of claim 7 , wherein performing hypothesis testing comprises performing hypothesis testing comprises removing a hypothesis from further testing based on a pseudo-likelihood of the hypothesis passing a threshold value.

9. The method of claim 8 , wherein performing hypothesis testing comprises ceasing hypothesis testing when a ratio of probability of a most likely hypothesis to probability of a second-most likely hypothesis passes a threshold value.

10. The method of claim 9 , further comprising generating a second set of integer phase ambiguity hypotheses after hypotheses from the first set of integer phase ambiguity hypotheses have been removed from further testing.

11. The method of claim 10 , wherein generating the second set of integer phase ambiguity hypotheses comprises generating the second set of integer phase ambiguity hypotheses in response to a hypothesis search space becoming smaller than a threshold value.

12. The method of claim 11 , wherein calculating the set of integer phase ambiguities comprises tracking relative probabilities in logarithmic space.

13. A method for Real Time Kinematic satellite positioning comprising:

at a mobile receiver, receiving a first navigation satellite carrier signal from a first navigation satellite, receiving a second navigation satellite carrier signal from a second navigation satellite, receiving a third navigation satellite carrier signal from a third navigation satellite, and receiving a fourth navigation satellite carrier signal from a fourth navigation satellite;

at the mobile receiver, transmitting carrier signal data corresponding to the first, second, and third navigation satellites from the mobile receiver to a remote computer;

at the remote computer, receiving first, second, third, and fourth phase correction signals from a reference station;

at the remote computer, calculating a set of integer phase ambiguities from double-differenced measurements of pseudo-range and phase; wherein calculating a set of integer phase ambiguities comprises generating a first set of integer phase ambiguity hypotheses using means and covariances generated by a Bierman-Thornton filter and at least one of a LAMBDA and an MLAMBDA algorithm; and performing hypothesis testing on the first set of integer phase ambiguity hypotheses; wherein generating the first set of integer phase ambiguity hypotheses comprises generating the first set of integer phase ambiguity hypotheses from a measurement equation that relates phase change and pseudo-range to integer ambiguity without inclusion of a baseline vector; wherein generating the first set of integer phase ambiguity hypotheses comprises generating the first set of integer phase ambiguity hypotheses without a dynamic transition model; and

at the remote computer, calculating a position of the mobile receiver from the set of integer phase ambiguities and the double-differenced measurements of pseudo-range and phase.

14. The method of claim 13 , wherein generating the first set of integer phase ambiguity hypotheses further comprises performing a decorrelating reparameterization of a hypothesis search space.

15. The method of claim 14 , wherein performing hypothesis testing comprises performing hypothesis testing using a Bayesian update algorithm.