Abstract:
A satellite navigation receiver and associated methods are described that can provide improved integer ambiguity resolution and more accurate positioning information. A modified BIE process may be utilized to enable the receiver to perform the integer ambiguity resolution more optimally. The output of the modified BIE process may be time-domain smoothed to provide a solution which is smoother in ambiguity space, and therefore also provide a position solution that is smoother in time. Transitions between an ambiguity-determined solution to a float solution, when necessary, may be smoothed in time. A weighting scheme may dynamically blend the ambiguity-determined solution and the float solution to leverage the advantages of both solutions, such as faster pull-in, higher accuracy, and more stable and smooth performance.

Description:
CROSS-REFERENCE TO RELATED APPLICATION 
       [0001]    This application claims priority from U.S. Provisional Patent Application Ser. No. 62/310,297, which was filed on Mar. 18, 2016, which is incorporated herein by reference. 
     
    
     TECHNICAL FIELD 
       [0002]    This application generally relates to a satellite navigation receiver having improved ambiguity resolution. In particular, this application relates to a precise point positioning real-time kinematic (PPP-RTK) satellite navigation receiver and associated methods that perform integer ambiguity resolution that is stable, robust, and accurate with fast initialization. 
       BACKGROUND 
       [0003]    Global navigation satellite systems (GNSS) utilize satellites to enable a receiver to determine position, velocity, and time with very high accuracy and precision using signals transmitted from the satellites. Such GNSS include the Global Positioning System (GPS), GLONASS, and Galileo. The signals transmitted from the satellites include one or more carrier signals at separate known frequencies, such as a first carrier (L1), a second carrier (L2), and an additional third carrier (L5) in the GPS. A code, such as a pseudo-random (PN) noise code modulated with information, may modulate a carrier of the signal, and may be unique to each satellite. Because the satellites have known orbital positions with respect to time, the signals can be used to estimate the relative position between an antenna of a receiver and each satellite, based on the propagation time of one or more signals received from four or more of the satellites. In particular, the receiver can synchronize a local replica of the carrier and code transmitted in a signal to estimate the relative position. 
         [0004]    The most accurate GNSS systems are referred to as precise point positioning real-time kinematic (PPP-RTK) or global RTK. The algorithms used in PPP-RTK systems are a combination of the algorithms used in local RTK systems and PPP systems. Both local RTK systems and PPP systems can achieve high accuracy by determining carrier phase related ambiguities. In local RTK systems, a roving receiver receives real-time corrections from a nearby local reference station, such as through a radio link. Because the local reference station has a known precise location, it can help determine the precise location of the roving receiver. In PPP systems, a roving receiver receives corrections that are globally applicable, which eliminates the need for local reference stations. The corrections can include information regarding the position and clock error of satellites, so that the roving receiver can receive information regarding the precise location of the satellites to help determine the precise location of the receiver. PPP systems have a global network of reference stations that are used to develop the global corrections, which are then transmitted to the roving receiver. 
         [0005]    PPP-RTK systems involve integer ambiguity resolution at the global network of reference stations and at the roving receiver. PPP-RTK systems are often used in applications such as precision farming, military navigation, and marine offshore positioning, due to its improved navigation accuracy and simplified infrastructure (i.e., eliminating the need for local reference stations). However, current PPP-RTK systems do not typically have real-time integer ambiguity resolution that is simultaneously stable, robust, and accurate, and with fast initialization times. 
         [0006]    Accordingly, there is an opportunity for a satellite navigation receiver that addresses these concerns. More particularly, there is an opportunity for a satellite navigation receiver and associated methods that can provide improved integer ambiguity resolution and more accurate positioning information. 
       SUMMARY 
       [0007]    The systems and methods described herein may result in a mobile receiver with improved integer ambiguity resolution and more accurate positioning information. The real-time integer ambiguity resolution described herein may be simultaneously stable, robust, and accurate, and have fast initialization times. For example, a modified version of the best integer equivariant (BIE) process may enable the mobile receiver to perform the integer ambiguity resolution more optimally. The modified BIE process described below may compute sums of weights and weighted sums of candidate narrow lane integer ambiguities during a search of the candidate narrow lane integer ambiguities. This may eliminate the need to store a large number of candidates or choose artificial thresholds to control the number of candidates to explore or store. The modified BIE process may also utilize an adaptive weight scaling during the search of the candidate narrow lane integer ambiguities. This may mitigate potential numerical issues attributable to a possible large dynamic range in weight magnitudes of the candidates. The modified BIE process may further utilize dynamic thresholds to control when to terminate the search of the candidates. In this way, only those candidates whose weight is large enough to have a meaningful numerical impact are included. 
         [0008]    Other features described herein also enable the mobile receiver to perform the integer ambiguity resolution more optimally. For example, the output of the modified BIE process may also be time-domain smoothed to provide a solution which is smoother in ambiguity space, and therefore also provide a position solution that is smoother in time. As another example, transitions between an ambiguity-determined solution to a float solution, when necessary, may be smoothed in time. As a further example, a weighting scheme may dynamically blend the ambiguity-determined solution and the float solution to leverage the advantages of both solutions, such as faster pull-in, higher accuracy, and more stable and smooth performance. The weighting scheme may utilize particular figures-of-merit and other heuristics to perform the blending. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0009]      FIG. 1  is a block diagram of a satellite navigation receiver, in accordance with some embodiments. 
           [0010]      FIG. 2  is a system diagram illustrating a method of resolving integer ambiguities in a satellite navigation receiver, in accordance with some embodiments. 
           [0011]      FIG. 3  is a flowchart illustrating operations for determining estimated float narrow lane ambiguities in a satellite navigation receiver, in accordance with some embodiments. 
           [0012]      FIG. 4  is a flowchart illustrating operations for calculating an ambiguity determined position solution using a modified best integer equivariant process in a satellite navigation receiver, in accordance with some embodiments. 
           [0013]      FIG. 5  is a flowchart illustrating operations for computing corrections to a state vector and a covariance matrix of a float solution in a satellite navigation receiver, in accordance with some embodiments. 
           [0014]      FIG. 6  is a flowchart illustrating operations for generating a smoothed ambiguity determined position solution in a satellite navigation receiver, in accordance with some embodiments. 
           [0015]      FIG. 7  is a flowchart illustrating operations for updating smoothed ambiguity values in a satellite navigation receiver, in accordance with some embodiments. 
           [0016]      FIG. 8  is a flowchart illustrating operations for weighting an ambiguity determined position solution in a satellite navigation receiver, in accordance with some embodiments. 
       
    
    
     DETAILED DESCRIPTION 
       [0017]    The description that follows describes, illustrates and exemplifies one or more particular embodiments of the invention in accordance with its principles. This description is not provided to limit the invention to the embodiments described herein, but rather to explain and teach the principles of the invention in such a way to enable one of ordinary skill in the art to understand these principles and, with that understanding, be able to apply them to practice not only the embodiments described herein, but also other embodiments that may come to mind in accordance with these principles. The scope of the invention is intended to cover all such embodiments that may fall within the scope of the appended claims, either literally or under the doctrine of equivalents. 
         [0018]    It should be noted that in the description and drawings, like or substantially similar elements may be labeled with the same reference numerals. However, sometimes these elements may be labeled with differing numbers, such as, for example, in cases where such labeling facilitates a more clear description. Additionally, the drawings set forth herein are not necessarily drawn to scale, and in some instances proportions may have been exaggerated to more clearly depict certain features. Such labeling and drawing practices do not necessarily implicate an underlying substantive purpose. As stated above, the specification is intended to be taken as a whole and interpreted in accordance with the principles of the invention as taught herein and understood to one of ordinary skill in the art. 
         [0019]      FIG. 1  shows a satellite navigation receiver  11  capable of receiving signals transmitted by satellites  10  that include one or more carrier signals (e.g., a first carrier (L1), a second carrier (L2) and an additional third carrier (L5) of the Global Positioning System (GPS)) such that the receiver  11  can determine position, velocity, and time with very high accuracy and precision based on the received signals. The received signals may be transmitted from one or more satellites  10 , such as a GPS satellite, a Galileo-compatible satellite, or a Global Navigation Satellite System (GLONASS) satellite. The satellites  10  have approximately known orbital positions versus time that can be used to estimate the relative position between an antenna  17  of the receiver  11  and each satellite  10 , based on the propagation time of one or more received signals between four or more of the satellites  10  and the antenna  17  of the receiver  11 . 
         [0020]    In any of the above referenced drawings of this document, any arrow or line that connects any blocks, components, modules, multiplexers, memory, data storage, accumulators, data processors, electronic components, oscillators, signal generators, or other electronic or software modules may comprise one or more of the following items: a physical path of electrical signals, a physical path of an electromagnetic signal, a logical path for data, one or more data buses, a circuit board trace, a transmission line; a link, call, communication, or data message between software modules, programs, data, or components; or transmission or reception of data messages, software instructions, modules, subroutines or components. 
         [0021]    In embodiments, the receiver  11  described herein may comprise a computer-implemented system or method in which one or more data processors process, store, retrieve, and otherwise manipulate data via data buses and one or more data storage devices (e.g., accumulators or memory) as described in this document and the accompanying drawings. As used in this document, “configured to, adapted to, or arranged to” mean that the data processor or receiver  11  is programmed with suitable software instructions, software modules, executable code, data libraries, and/or requisite data to execute any referenced functions, mathematical operations, logical operations, calculations, determinations, processes, methods, algorithms, subroutines, or programs that are associated with one or more blocks set forth in  FIG. 1  and/or any other drawing in this disclosure. Alternately, separately from, or cumulatively with the above definition, “configured to, adapted to, or arranged to” can mean that the receiver  11  comprises one or more components described herein as software modules, equivalent electronic hardware modules, or both to execute any referenced functions, mathematical operations, calculations, determinations, processes, methods, algorithms, or subroutines. 
         [0022]    Precise point positioning (PPP) includes the use of precise satellite orbit and clock corrections provided wirelessly via correction data, rather than through normal satellite broadcast information (ephemeris and clock data) that is encoded on the received satellite signals, to determine a relative position or absolute position of a mobile receiver. PPP may use correction data that is applicable to a wide geographic area. Although the resulting positions can be accurate within a few centimeters using state-of-the-art algorithms, conventional precise point positioning can have a long convergence time of up to tens of minutes to stabilize and determine the float or integer ambiguity values necessary to achieve the purported (e.g., advertised) steady-state accuracy. Hence, such long convergence time is typically a limiting factor in the applicability of PPP. 
         [0023]    As shown in  FIG. 1 , the receiver  11  may include a receiver front-end module  40  coupled to an electronic data processing system  129 . Furthermore, a correction wireless device  44  (e.g., a receiver or transceiver) may provide correction data or differential correction data (e.g., PPP correction data) to enhance the accuracy of position estimates provided or estimated by the receiver  11 . 
         [0024]    In an embodiment, the receiver front-end module  40  includes a radio frequency (RF) front end  42  coupled to an analog-to-digital converter  46 . The receiver front-end module  40  or the RF front end  42  may receive a set of carrier signals from one or more satellite transmitters on satellites. The analog-to-digital converter  46  may convert the set of carrier signals into digital signals, such as digital baseband signals or digital intermediate frequency signals for processing by the electronic data processing system  129 . 
         [0025]    In an embodiment, the electronic data processing system  129  includes a baseband processing module  48  (e.g., baseband/intermediate frequency processing module) and a navigation positioning estimator  50 . For example, the baseband processing module  48  and the navigation positioning estimator  50  may be stored in a data storage device  155 . 
         [0026]    In an embodiment, the baseband processing module  48  may include a measurement module  161  that includes a carrier phase measurement module  151  and/or a code phase measurement module  153 . The carrier phase measurement module  151  may facilitate the measurement of the carrier phase of one or more carrier signals received by the receiver  11 . Similarly, the code phase measurement module  153  may facilitate the measurement of the code phase of one or more code signals that modulate the carrier signals received by the receiver  11 . 
         [0027]    The navigation positioning estimator  50  can use the carrier phase measurements and/or the code phase measurements to estimate the range between the receiver  11  and one or more satellites, or estimate the position (e.g., three dimensional coordinates) of the receiver  11  with respect to one or more satellites (e.g., four or more satellites). The ambiguities refer to the differences in the measurements, such as between-satellite single differences. The code phase measurements or carrier phase measurements can be converted from propagation times, between each satellite and the receiver  11  that is within reception range of the receiver, to distances by dividing the propagation time by the speed of light, for example. 
         [0028]    In the electronic data processing system  129 , the data storage device  155  may be coupled to a data bus  157 . An electronic data processor  159  may communicate with the data storage device  155  and the correction wireless device  44  via the data bus  157 . As used herein, the data processor  159  may include one or more of the following: an electronic data processor, a microprocessor, a microcontroller, an application specific integrated circuit (ASIC), digital signal processor (DSP), a programmable logic device, an arithmetic logic unit, or another electronic data processing device. The data storage device  155  may include electronic memory, registers, shift registers, volatile electronic memory, a magnetic storage device, an optical storage device, or any other device for storing data. 
         [0029]    In an embodiment, the navigation positioning estimator  50  includes a precise position estimator, such as a precise point position (PPP) estimator or a wide area differential global navigation satellite system (GNSS) position estimator. The navigation positioning estimator  50  may receive correction data from the correction wireless device  44 , which is a receiver or transceiver capable of communication with a wireless satellite communications device. 
         [0030]    The mobile receiver described herein assumes that two frequencies are available and used for navigation. However, it is contemplated that the described concepts may be extended to cover scenarios with more than two frequencies, and to be usable with any GNSS systems. 
         [0031]    The following description uses a notation system where individual terms may be related to a specific frequency, satellite, or receiver. The notation uses subscripts and superscripts to distinguish these elements and uses the location of the subscript or superscript to distinguish the elements, where frequency is designated by a numerical right subscript, receiver is designated by a left subscript, and satellite is designated by a left superscript. For example, the term  R   k X i  refers to a frequency I, a receiver R, and a satellite k. A right superscript retains the usual meaning of an exponent. However, not all subscripts and superscripts may be designated for each term. As such, when an element is not relevant to the context of a particular equation, the subscript and/or superscript may be dropped. 
         [0032]    In addition, the following description utilizes parameters and notation, including f i  as a frequency in hertz, λ i  as a wavelength of f i  in meters, P i  as a measured pseudorange in meters, Φ i  as a measured carrier phase in cycles, N as an integer number of ambiguity cycles,           as a non-integer (float) number of ambiguity cycles, and c as the speed of light in meters/second. Differences between satellite pairs are represented using ∇ with left superscripts indicate the satellites involved. For example, the term  i,j ∇x is equal to  i x− j x. 
         [0033]      FIG. 2  shows a system diagram of operations within the receiver  11  for resolving integer ambiguities. The operations are shown as functional blocks within the baseband processing module  48  and the navigation positioning estimator  50 . The text between the functional blocks generally denotes the output generated by a block. 
         [0034]    Baseband processing at block  301  may be performed by the baseband processing module  48 . The baseband processing may measure the pseudorange P i  and carrier phase Φ i  of one or more received satellite signals. The uncorrected pseudorange P i  and carrier phase Φ i  measurements, where the frequency i=1, 2, . . . , may be given as: 
         [0000]    
       
         
           
             
               
                 
                   
                     P 
                     i 
                   
                   = 
                   
                     D 
                     + 
                     
                       B 
                       i 
                         
                         
                         
                       k 
                     
                     + 
                     
                       B 
                       i 
                         
                         
                       R 
                         
                     
                     + 
                     
                       I 
                       
                         f 
                         i 
                         2 
                       
                     
                     + 
                     
                       ε 
                       
                         p 
                         , 
                         i 
                       
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       Φ 
                       i 
                     
                      
                     
                       λ 
                       i 
                     
                   
                   = 
                   
                     D 
                     + 
                     
                       b 
                       i 
                         
                         
                         
                       k 
                     
                     + 
                     
                       b 
                       i 
                         
                         
                       R 
                         
                     
                     - 
                     
                       I 
                       
                         f 
                         i 
                         2 
                       
                     
                     + 
                     
                       δ 
                        
                       
                           
                       
                        
                       
                         PW 
                         i 
                       
                     
                     + 
                     
                       
                         N 
                         i 
                       
                        
                       
                         λ 
                         i 
                       
                     
                     + 
                     
                       ε 
                       
                         Φ 
                         , 
                         i 
                       
                     
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
             
               
                 
                   
                     D 
                     = 
                     
                       ρ 
                       + 
                       
                         
                             
                           k 
                         
                          
                         τ 
                       
                       + 
                       
                         
                             
                           R 
                         
                          
                         τ 
                       
                       + 
                       T 
                       + 
                       
                         δ 
                          
                         
                             
                         
                          
                         PC 
                       
                       + 
                       
                         δ 
                          
                         
                             
                         
                          
                         T 
                       
                       + 
                       
                         δ 
                          
                         
                             
                         
                          
                         R 
                       
                       + 
                       
                         δ 
                          
                         
                             
                         
                          
                         S 
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where k=1, 2, . . . is the index of the satellite; ρ is the geometric range in meters; B is the pseudorange bias in meters; b is the carrier phase bias in meters; I is the ionospheric error in meters-Hz 2 ; ε p,i  is the pseudorange noise error in meters (including white noise, multipath, and remaining modeling errors); ε Φ,i  is the carrier measurement noise error in meters (including white noise, multipath, and remaining modeling errors); τ is the clock error in meters;  R τ is the receiver clock error and is specific to a given GNSS system; T is the tropospheric delay in meters; δPC is the phase center offset and variation in meters; δT is the error due to tidal forces and polar motion in meters; δR is the relativistic effect on satellite clock in meters; δS is the relativistic effect on signal propagation (the Shapiro delay) in meters; and δPW is the phase windup error in meters. 
         [0035]    An alternative to having a receiver clock error  R τ for each GNSS system is to estimate one clock for a designated primary constellation (e.g., GPS) and relative receiver clock offsets between the primary constellation and the other GNSS constellations. The tropospheric delay T is typically divided into a dry component T dry  and a wet component T wet . The dry component T dry  can be accurately modeled using an a priori troposphere model, such as GPT2 (Global Pressure and Temperature). The remaining wet component T wet  after removing an a priori wet model can be further estimated by one zenith bias with a mapping function bias and/or two additional horizontal gradient coefficients. 
         [0036]    In blocks  302  and  303 , the pseudorange P i  and the carrier phase Φ i  from the baseband processing block  301  may be processed with models and measurement combinations to eliminate and/or reduce a subset of the error terms in equations (1), (2), and (3). For the ranging codes and the navigation message to travel from a satellite  10  to the receiver, they must be modulated onto a carrier frequency. In the case of GPS, two frequencies are utilized: one at 1575.42 MHz (10.23 MHz×154) called L1; and a second at 1227.60 MHz (10.23 MHz×120) called L2. Both L1 and L2 are in the satellite L-band. 
         [0037]    The signals transmitted by GLONASS satellites  10  are derived from the fundamental frequencies of 1602 MHz for L1 and 1246 MHz for L2. Each GLONASS satellite  10  transmits on a different frequency using FDMA (frequency division multiple access) and according to a designated frequency channel number. The L1 center frequency for GLONASS is given by: 
         [0000]        k   f   1 =1602 MHz+ k   n× 0.5625 MHz,  (4)
 
         [0000]    where  k n is the frequency channel number of satellite k, and where  k nε{−7, −6, . . . , 6}. The L2 center frequency for GLONASS is given by: 
         [0000]        k   f   2 =1246 MHz+ k   n× 0.4375 MHz.  (5)
 
         [0038]    In PPP systems, a float solution is based on processing ionosphere-free (IF) combinations of both the pseudorange P i  and the carrier phase Φ i  on the two frequencies, as given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     P 
                     IF 
                   
                    
                   
                     = 
                     △ 
                   
                    
                   
                     
                       
                         
                           f 
                           1 
                           2 
                         
                         
                           ( 
                           
                             
                               f 
                               1 
                               2 
                             
                             - 
                             
                               f 
                               2 
                               2 
                             
                           
                           ) 
                         
                       
                        
                       
                         P 
                         1 
                       
                     
                     - 
                     
                       
                         
                           f 
                           2 
                           2 
                         
                         
                           ( 
                           
                             
                               f 
                               1 
                               2 
                             
                             - 
                             
                               f 
                               2 
                               2 
                             
                           
                           ) 
                         
                       
                        
                       
                         P 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
             
               
                 
                   
                     P 
                     IF 
                   
                   = 
                   
                     D 
                     + 
                     
                       B 
                       IF 
                         
                         
                         
                       k 
                     
                     + 
                     
                       B 
                       IF 
                         
                         
                       R 
                         
                     
                     + 
                     
                       ε 
                       
                         P 
                         , 
                         IF 
                       
                         
                         
                         
                       k 
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
             
               
                 
                   
                     Φ 
                     IF 
                   
                    
                   
                     = 
                     △ 
                   
                    
                   
                     
                       
                         
                           f 
                           1 
                           2 
                         
                         
                           ( 
                           
                             
                               f 
                               1 
                               2 
                             
                             - 
                             
                               f 
                               2 
                               2 
                             
                           
                           ) 
                         
                       
                        
                       
                         Φ 
                         1 
                       
                        
                       
                         λ 
                         1 
                       
                     
                     - 
                     
                       
                         
                           f 
                           2 
                           2 
                         
                         
                           ( 
                           
                             
                               f 
                               1 
                               2 
                             
                             - 
                             
                               f 
                               2 
                               2 
                             
                           
                           ) 
                         
                       
                        
                       
                         Φ 
                         2 
                       
                        
                       
                         λ 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       
                         Φ 
                         IF 
                       
                        
                       
                         λ 
                         IF 
                       
                     
                     = 
                     
                       D 
                       + 
                       
                         b 
                         IF 
                           
                           
                           
                         k 
                       
                       + 
                       
                         b 
                         IF 
                           
                           
                         R 
                           
                       
                       + 
                       
                         δ 
                          
                         
                             
                         
                          
                         
                           PW 
                           IF 
                         
                       
                       + 
                       
                         
                           N 
                           IF 
                         
                          
                         
                           λ 
                           IF 
                         
                       
                       + 
                       
                         ε 
                         
                           Φ 
                           , 
                           IF 
                         
                           
                           
                           
                         k 
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
           
         
       
     
         [0000]    where  R B IF  is the receiver ionosphere-free code bias, which is the ionosphere-free combination of the L1 receiver code bias and the L2 receiver code bias. There is a receiver ionosphere-free code bias per receiver and constellation for all visible CDMA satellites. 
         [0039]    For GLONASS satellites, an additional inter-channel code bias may need to be estimated, if the magnitude of the inter-channel code bias is significant. In this case, the ionosphere-free pseudorange measurement is given as: 
         [0000]        P   IF   =D+   k   B   IF + R   B   IF + R   k   C   GLN + k ε P,IF ,  (10)
 
         [0000]    where C GLN  is the GLONASS code bias in meters. 
         [0040]    In PPP systems, one goal is to have a coherent model for receiver clock and bias terms. The measurement compensation may include compensating the measurements using broadcast satellite ephemeris and clock, compensating the measurements for the deterministic terms (e.g., δPC, δT, δR, δS, and δPW), and compensating the measurements for the PPP corrections for the satellite orbit and clock. 
         [0041]    In blocks  302  and  303 , it is assumed there is a common receiver clock term for both the pseudorange P i  and the carrier phase Φ i . The receiver ionosphere-free code bias  R B IF  may be considered a nuisance parameter and be naturally absorbed into the receiver clock error  R τ. In addition, the PPP corrections for the satellite clocks inherently account for the satellite pseudorange bias terms  k B IF  (but not for the receiver-dependent GLONASS channel bias or the system bias between GPS and GLONASS). The receiver carrier phase bias  R b IF  may not be easily estimated separately in the float solution and is therefore absorbed into each of the resulting float ambiguity terms. The PPP corrections include additional terms for each satellite that enable compensation of each measurement for satellite carrier phase biases  k b IF , which are not constant over time. 
         [0042]    In block  304 , a recursive estimator (e.g., a Kalman filter) may compute a float solution and corresponding zero-difference ionospheric-free float ambiguity values. The float solution may consist of a state vector X and a covariance matrix P for terms such as position, clock bias, tropospheric delay, and floating ambiguity values. The position of the receiver may be updated at each interval (e.g., epoch) using the recursive estimator in block  304 , based on the compensated ionosphere-free measurements from block  303 . The float solution from block  304  does not itself involve ambiguity resolution. 
         [0043]    The float solution may be determined using simplified ionosphere-free measurement equations. For GPS, such equations are given as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       P 
                       ~ 
                     
                     IF 
                   
                   = 
                   
                     
                       
                         - 
                         
                           
                               
                             k 
                           
                            
                           
                             e 
                             ⇀ 
                           
                         
                       
                       · 
                       
                         x 
                         ⇀ 
                       
                     
                     + 
                     τ 
                     + 
                     
                       
                         M 
                          
                         
                           ( 
                           
                             
                                 
                               k 
                             
                              
                             E 
                           
                           ) 
                         
                       
                        
                       
                         T 
                         wet 
                       
                     
                     + 
                     
                       ε 
                       
                         
                           P 
                           ~ 
                         
                         IF 
                       
                         
                         
                         
                       k 
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       
                         Φ 
                         ~ 
                       
                       IF 
                     
                      
                     
                       λ 
                       IF 
                     
                   
                   = 
                   
                     
                       
                         - 
                         
                           
                               
                             k 
                           
                            
                           
                             e 
                             ⇀ 
                           
                         
                       
                       · 
                       
                         x 
                         ⇀ 
                       
                     
                     + 
                     τ 
                     + 
                     
                       M 
                        
                       
                         ( 
                         
                           
                               
                             k 
                           
                            
                           E 
                         
                         ) 
                       
                        
                       
                         T 
                         wet 
                       
                     
                     + 
                     
                       
                          
                         IF 
                           
                           
                           
                         k 
                       
                        
                       
                         λ 
                         NL 
                       
                     
                     + 
                     
                       
                         ε 
                         
                           
                             Φ 
                             ~ 
                           
                           IF 
                         
                           
                           
                           
                         k 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   12 
                   ) 
                 
               
             
           
         
       
     
         [0000]    For GLONASS, such equations are given as: 
         [0000]    
       
         
           
             
               
                 
                   
                       
                   
                    
                   
                     
                       
                         P 
                         ~ 
                       
                       IF 
                     
                     = 
                     
                       
                         
                           - 
                           
                             
                                 
                               k 
                             
                              
                             
                               e 
                               ⇀ 
                             
                           
                         
                         · 
                         
                           x 
                           ⇀ 
                         
                       
                       + 
                       τ 
                       + 
                       
                         Δτ 
                         GLN 
                       
                       + 
                       
                         C 
                         GLN 
                           
                           
                           
                         k 
                       
                       + 
                       
                         
                           M 
                            
                           
                             ( 
                             
                               
                                   
                                 k 
                               
                                
                               E 
                             
                             ) 
                           
                         
                          
                         
                           T 
                           wet 
                         
                       
                       + 
                       
                         ε 
                         
                           
                             P 
                             ~ 
                           
                           IF 
                         
                           
                           
                           
                         k 
                       
                     
                   
                 
               
               
                 
                   ( 
                   13 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       
                         Φ 
                         ~ 
                       
                       IF 
                     
                      
                     
                       λ 
                       IF 
                     
                   
                   = 
                   
                     
                       
                         - 
                         
                           
                               
                             k 
                           
                            
                           
                             e 
                             ⇀ 
                           
                         
                       
                       · 
                       
                         x 
                         ⇀ 
                       
                     
                     + 
                     τ 
                     + 
                     
                       Δτ 
                       GLN 
                     
                     + 
                     
                       M 
                        
                       
                         ( 
                         
                           
                               
                             k 
                           
                            
                           E 
                         
                         ) 
                       
                        
                       
                         T 
                         wet 
                       
                     
                     + 
                     
                       
                          
                         IF 
                           
                           
                           
                         k 
                       
                        
                       
                         λ 
                         NL 
                       
                     
                     + 
                     
                       
                         ε 
                         
                           
                             Φ 
                             ~ 
                           
                           IF 
                         
                           
                           
                           
                         k 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   14 
                   ) 
                 
               
             
           
         
       
     
         [0044]    In equations (11)-(14), {right arrow over (x)} is the receiver position,  k {right arrow over (x)} is the position of satellite k, and {right arrow over (e)} is the receiver to satellite line-of-sight vector, where {right arrow over (e)}={right arrow over (r)}/∥{right arrow over (r)}∥ and {right arrow over (r)}         k{right arrow over (x)}−{right arrow over (x)}. Also, τ is the receiver clock error (relative to GPS) in meters; T wet  is the residual zenith tropospheric wet delay; E is the elevation angle from the receiver to the satellite; M(•) is the elevation wet mapping function that maps zenith tropospheric delay to the line-of-sight;  k             IF  is the float ambiguity; λ NL  is the narrow lane wavelength and is defined as 
         [0000]    
       
         
           
             
               
                 λ 
                 NL 
               
                
               
                 = 
                 △ 
               
                
               
                 
                   c 
                   
                     ( 
                     
                       
                         f 
                         1 
                       
                       + 
                       
                         f 
                         2 
                       
                     
                     ) 
                   
                 
                 = 
                 
                   
                     λ 
                     IF 
                   
                    
                   
                     ( 
                     
                       
                         f 
                         1 
                       
                       - 
                       
                         f 
                         2 
                       
                     
                     ) 
                   
                 
               
             
             ; 
           
         
       
     
         [0000]    ΔT GLN  is a slowly varying term for the system bias between GPS and GLONASS;  k C GLN  is the GLONASS ionosphere-free code bias in meters;  k ε {tilde over (P)}     IF    is the ionosphere-free pseudorange measurement noise error in meters (including white noise, multipath, and remaining modeling errors); and  k ε {tilde over (Φ)}     IF    is the ionosphere-free carrier measurement noise error in meters (including white noise, multipath, and remaining modeling errors). The narrow lane wavelength λ NL  is utilized in equations (11)-(14) instead of an ionospheric-free wavelength because the ionospheric-free wavelength is relatively short, which causes difficulties in directly resolving integer ambiguities. Accordingly, the float ambiguity term satisfies: 
         [0000]        k             IF λ NL   =N   IF λ IF + R   b   IF   (15)
 
         [0045]    The state vector of Kalman filter may consist of the collection of elements: receiver position {right arrow over (x)}, receiver clock error τ, system bias Δτ GLN , ionosphere-free code bias  k C GLN , tropospheric delay T, and float ambiguity  k             IF . The total active states may include three states for the receiver velocity {dot over ({right arrow over (x)})} and be given by 6+1+1+N GLN +1+(N GPS +N GLN ). There may be a total of 9+N GPS +2N GLN  active states in the Kalman filter at any time, where N GPS  and N GLN  represent the number of GPS and GLONASS satellites in view of the receiver, respectively. The number of states may increase by two if tropospheric gradient terms are included. 
         [0046]    The Kalman filter may have time update and measurement update operations, as is known in the art. The process noise added to the states may include a small amount of fully correlated noise (e.g., 0.04 cycles 2  per second) because the receiver phase bias  R b IF  has been absorbed into the zero-difference floating ambiguity states. 
         [0047]    The state vector X and the covariance matrix P of the Kalman filter can be referred to as the float solution. The purpose of ambiguity determination as described herein is to result in corrections to the float solution (i.e., ΔX, ΔP). The corrected state vector may be given by X+ΔX with covariance P−ΔP. The corrected state vector may contain an improved estimate of the position, among other things, such as clock bias, tropospheric delay, and floating ambiguity values. 
         [0048]    Ambiguity determination may be performed using wide lane and narrow lane measurement combinations, due to the relatively small ionosphere-free wavelength λ IF , e.g., approximately 0.6 cm for GPS. The wide lane ambiguities may be determined in block  305  and the narrow lane ambiguities may be determined in block  306 . The wide lane wavelength may be given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     λ 
                     WL 
                   
                    
                   
                     = 
                     △ 
                   
                    
                   
                     c 
                     
                       ( 
                       
                         
                           f 
                           1 
                         
                         - 
                         
                           f 
                           2 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   16 
                   ) 
                 
               
             
           
         
       
     
         [0000]    and the narrow lane wavelength may be given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     λ 
                     NL 
                   
                    
                   
                     = 
                     △ 
                   
                    
                   
                     
                       c 
                       
                         ( 
                         
                           
                             f 
                             1 
                           
                           + 
                           
                             f 
                             2 
                           
                         
                         ) 
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   17 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Accordingly, the wide lane wavelength may be approximately 86.2 cm for GPS and 84.2 cm for GLONASS, and the narrow lane wavelength may be approximately 10.7 cm for GPS and 10.5 cm for GLONASS. 
         [0049]    The wide lane ambiguity N WL  is defined as: 
         [0000]        N   WL             N   1   −V   2   (18)
 
         [0000]    and may be determined first, and the narrow lane ambiguity N NL  may be determined based on the wide lane ambiguity N WL . In particular, because: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         λ 
                         IF 
                       
                        
                       
                         N 
                         IF 
                       
                     
                     = 
                     
                       
                         
                           λ 
                           NL 
                         
                          
                         
                           1 
                           
                             ( 
                             
                               
                                 f 
                                 1 
                               
                               - 
                               
                                 f 
                                 2 
                               
                             
                             ) 
                           
                         
                          
                         
                           N 
                           IF 
                         
                       
                       = 
                       
                         
                           λ 
                           NL 
                         
                          
                         
                           ( 
                           
                             
                               
                                 
                                   f 
                                   1 
                                 
                                 
                                   ( 
                                   
                                     
                                       f 
                                       1 
                                     
                                     - 
                                     
                                       f 
                                       2 
                                     
                                   
                                   ) 
                                 
                               
                                
                               
                                 N 
                                 1 
                               
                             
                             - 
                             
                               
                                 
                                   f 
                                   2 
                                 
                                 
                                   ( 
                                   
                                     
                                       f 
                                       1 
                                     
                                     - 
                                     
                                       f 
                                       2 
                                     
                                   
                                   ) 
                                 
                               
                                
                               
                                 N 
                                 2 
                               
                             
                           
                           ) 
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   19 
                   ) 
                 
               
             
           
         
       
     
         [0000]    the narrow lane ambiguity (i.e., any one of N 1 , N 2 , or N NL ) can be found using one of the following relationships: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       λ 
                       IF 
                     
                      
                     
                       N 
                       IF 
                     
                   
                   = 
                   
                     
                       
                         λ 
                         NL 
                       
                        
                       
                         N 
                         1 
                       
                     
                     - 
                     
                       
                         λ 
                         NL 
                       
                        
                       
                         
                           f 
                           2 
                         
                         
                           ( 
                           
                             
                               f 
                               1 
                             
                             - 
                             
                               f 
                               2 
                             
                           
                           ) 
                         
                       
                        
                       
                         N 
                         WL 
                       
                     
                   
                 
               
               
                 
                   ( 
                   20 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       λ 
                       IF 
                     
                      
                     
                       N 
                       IF 
                     
                   
                   = 
                   
                     
                       
                         λ 
                         NL 
                       
                        
                       
                         N 
                         2 
                       
                     
                     + 
                     
                       
                         λ 
                         NL 
                       
                        
                       
                         
                           f 
                           1 
                         
                         
                           ( 
                           
                             
                               f 
                               1 
                             
                             - 
                             
                               f 
                               2 
                             
                           
                           ) 
                         
                       
                        
                       
                         N 
                         WL 
                       
                     
                   
                 
               
               
                 
                   ( 
                   21 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       
                         λ 
                         IF 
                       
                        
                       
                         N 
                         IF 
                       
                     
                     = 
                     
                       
                         
                           λ 
                           NL 
                         
                          
                         
                           
                             N 
                             NL 
                           
                           2 
                         
                       
                       + 
                       
                         
                           λ 
                           WL 
                         
                          
                         
                           
                             N 
                             WL 
                           
                           2 
                         
                       
                     
                   
                   , 
                   
                     
 
                   
                    
                   
                     
                       N 
                       NL 
                     
                      
                     
                       = 
                       def 
                     
                      
                     
                       
                         N 
                         1 
                       
                       + 
                       
                         
                           N 
                           2 
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   22 
                   ) 
                 
               
             
           
         
       
     
         [0050]    Accordingly, once the wide lane ambiguity N WL  is determined, equations (20)-(22) may be utilized to find an expression to determine a narrow lane ambiguity. For example, using N 2  as the narrow lane ambiguity, N 2  can be found by rewriting equation (21) as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       λ 
                       NL 
                     
                      
                     
                       N 
                       2 
                     
                   
                   = 
                   
                     
                       
                         λ 
                         IF 
                       
                        
                       
                         N 
                         IF 
                       
                     
                     - 
                     
                       
                         λ 
                         NL 
                       
                        
                       
                         
                           f 
                           1 
                         
                         
                           ( 
                           
                             
                               f 
                               1 
                             
                             - 
                             
                               f 
                               2 
                             
                           
                           ) 
                         
                       
                        
                       
                         
                           N 
                           WL 
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   23 
                   ) 
                 
               
             
           
         
       
     
         [0000]    Because the narrow lane wavelengths λ NL  are much longer than the ionosphere-free wavelengths λ IF , ambiguity resolution is more easily performed. It should be noted that equations (20) and (22) could also be rewritten to find N 1  or N NL  as the narrow lane ambiguity. Without loss of generality, in the description that follows, N 2  is used as the narrow lane ambiguity. For simplicity and clarity, the subscript “NL” is used for the narrow lane ambiguity. 
         [0051]    Accordingly, at block  305 , the between-satellite single-difference wide lane ambiguities N WL  may be resolved using the Melbourne-Wubbena combination. The Melbourne-Wubbena combination is a geometry-free, ionospheric-free linear combination of phase and code measurements from a single receiver, and is given as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         λ 
                         WL 
                       
                        
                       MW 
                     
                      
                     
                       = 
                       Δ 
                     
                      
                     
                       
                         
                           λ 
                           WL 
                         
                          
                         
                           ( 
                           
                             
                               Φ 
                               1 
                             
                             - 
                             
                               Φ 
                               2 
                             
                           
                           ) 
                         
                       
                       - 
                       
                         ( 
                         
                           
                             
                               
                                 f 
                                 1 
                               
                               
                                 ( 
                                 
                                   
                                     f 
                                     1 
                                   
                                   + 
                                   
                                     f 
                                     2 
                                   
                                 
                                 ) 
                               
                             
                              
                             
                               P 
                               1 
                             
                           
                           + 
                           
                             
                               
                                 f 
                                 2 
                               
                               
                                 ( 
                                 
                                   
                                     f 
                                     1 
                                   
                                   + 
                                   
                                     f 
                                     2 
                                   
                                 
                                 ) 
                               
                             
                              
                             
                               P 
                               2 
                             
                           
                         
                         ) 
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   24 
                   ) 
                 
               
             
           
         
       
     
         [0000]    and can be written for GPS as: 
         [0000]      λ WL   MW=λ   WL   N   WL + k   B   WL + R   B   WL +ε WL ,  (25)
 
         [0000]    and written for GLONASS as: 
         [0000]      λ WL   MW=λ   WL   N   WL + k   B   WL + R   B   WL + k   IFB   WL +ε WL ,  (26)
 
         [0000]    where  k B WL  and  R B WL  are the satellite and receiver wide lane biases, respectively, that are a collection of the original biases with various scaling factors. The term  k IFB WL  represents the inter-frequency bias term that models the effect of GLONASS code channel biases on the wide lane measurement combination. The inter-frequency bias may vary from receiver to receiver, and may also vary in different installations (e.g., due to varying antenna and cabling setups). The magnitude of the inter-frequency bias is typically less than 0.1 cycles per frequency number difference. It can be assumed that the GLONASS code bias-related terms can be accurately modeled by a term that is linear in the GLONASS frequency number. As such, the inter-frequency bias  k IFB WL  may be approximately equal to K* k n, where  k nε{−7, 6, . . . , 6} and K is an unknown slowly varying coefficient for a given receiver. 
         [0052]    In block  305 , undifferenced Melbourne-Wubbena measurements may be used to estimate one wide lane ambiguity state per visible satellite. Typically, the wide lane satellite biases  k B WL  are broadcast in real-time within correction data and can be used for measurement compensation. The receiver wide lane biases  R B WL  may be lumped into the float wide lane ambiguity state            WL . Accordingly, the float wide lane ambiguity state            WL  is no longer an integer. However, the between-satellite single-differences  k,m ∇           WL , for GPS are integers and can be resolved in their single-difference form. For GLONASS, it is necessary to remove the inter-frequency bias contribution from the single-difference form in order to recover the integer. 
         [0053]    Because the float ambiguity states contain the wide lane ambiguity and the receiver bias, some amount of fully correlated process noise is typically applied in the dynamic update of the Kalman filter. The single differential wide-lane ambiguity and the variance-covariance can be derived based on the undifferenced float ambiguity states and variance-covariance in block  305  after a reference satellite for each constellation is chosen. A standard ambiguity resolution process can be applied for single-difference ambiguities  k,m ∇           WL . Techniques for solving this type of ambiguity resolution are known in the art. Ambiguity resolution validation may also be performed in block  305 . After ambiguity resolution validation is performed, a single-difference integer ambiguity constraint can be applied to wide lane float estimator. 
         [0054]    Narrow lane ambiguities may be determined based on the wide lane ambiguities determined in block  305  and the ionospheric-free float ambiguity values from block  304 . Ultimately, single-difference float narrow lane ambiguities may be computed that are used to fix or determine precise ambiguity values. The updated values may then be used to correct the state vector and covariance matrix of the float solution (i.e., compute ΔX, ΔP) in order to update the position of the receiver. 
         [0055]    In block  306 , initial estimates of the narrow lane ambiguities may be computed. The steps performed in block  306  are shown in the process  306  of  FIG. 3 . At step  402 , subsets of visible satellites may be selected as candidates for processing the narrow lane ambiguities and as a reference satellite for between-satellite single-differencing calculations. The subsets of the visible satellites may be selected based on measurement residuals from processing the respective ionosphere-free code and carrier measurements in the float Kalman filter (i.e., equations (13) and (14)), float ambiguity covariance for the respective float ambiguity values, PPP correction quality for the given satellites, and wide lane fixing status for the satellites (i.e., fixed or not fixed). For example, when the measurement residuals are high, this indicates possible problems with the associated satellites, and can result in not selecting those satellites and their estimates and/or measurements. 
         [0056]    At step  404 , it may be determined whether a condition exists for transitioning to the float solution from block  304  as the narrow lane ambiguity values. The conditions at step  404  may include whether the age of the PPP corrections is above a certain threshold (e.g., three minutes), whether there are not enough satellites available as candidates, and whether a suitable reference satellite is available. If such a condition exists at step  404 , then the process  306  continues to step  410  to transition to the float solution. 
         [0057]    Transitioning to the float solution at step  410  may be performed to ensure that the transition is relatively smooth and not too rapid. Because the state vector after applying the correction due to ambiguity determination is given by X+ΔX, the difference between the float solution and the solution after ambiguity determination is given by the offset ΔX and using the float solution is equivalent to setting the offset ΔX to zero. Therefore, an offset ΔX t  may be utilized to perform the transition, where the offset ΔX t  has been stored at a previous interval. Starting at the interval when the transition is begun, the subsequent position change may be limited at each interval by a term related to the offset ΔX t . In particular, the position change may be limited to not vary by more than a magnitude ε that is a small predetermined value. The transition may be made over N steps, where N is the rounded-off value of ∥ΔX t ∥/ε. Accordingly, at step  401 , for k=1 . . . N, the float solution can be transitioned to by the end of the transition period, where 
         [0000]    
       
         
           
             Δ 
              
             
                 
             
              
             
               X 
               
                 t 
                 + 
                 k 
               
             
              
             
               
                 N 
                 - 
                 k 
               
               N 
             
              
             Δ 
              
             
                 
             
              
             
               
                 X 
                 t 
               
               . 
             
           
         
       
     
         [0058]    Returning to step  404  in  FIG. 3 , if a transition condition does not exist, then the process  306  continues to step  406 . At step  406 , a reference satellite may be selected for each constellation. The selection of the reference satellite may be based on the float ambiguity covariance P and other heuristics, such as considering only the satellites for which wide lane ambiguity values have been determined and favoring the satellites which were used successfully for ambiguity determination in a previous interval. At step  408 , estimated float narrow lane ambiguities  k,m ∇           NL  may be determined that are between-satellite single-differences. For each constellation with m denoting the index of the selected reference satellite, the estimated float narrow lane ambiguities  k,m ∇           NL  may be given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                         
                       
                         k 
                         , 
                         m 
                       
                     
                      
                     
                       ∇ 
                       
                          
                         NL 
                       
                     
                   
                    
                   
                     = 
                     Δ 
                   
                    
                   
                     
                       
                           
                         
                           k 
                           , 
                           m 
                         
                       
                        
                       
                         ∇ 
                         
                            
                           IF 
                         
                       
                     
                     - 
                     
                       
                         
                           f 
                           1 
                         
                         
                           
                             f 
                             1 
                           
                           - 
                           
                             f 
                             2 
                           
                         
                       
                        
                       
                         
                           
                               
                             
                               k 
                               , 
                               m 
                             
                           
                            
                           
                             ∇ 
                             
                                
                               WL 
                             
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   27 
                   ) 
                 
               
             
           
         
       
     
         [0000]    The estimated float narrow lane ambiguities  k,m ∇           NL  may be considered as noisy measurements of the integers  k,m ∇N NL . 
         [0059]    Returning to  FIG. 2 , a modified BIE algorithm may be performed at block  307  to calculate the best estimate for narrow lane single-difference ambiguities  k,m ∇           NL,BIE , based on the estimated float narrow lane ambiguities  k,m ∇           NL  and the corresponding covariance matrix Q NL . Generally, these modified BIE single differences are non-integer. The modified BIE algorithm may be based on the LAMBDA technique to solve an integer least squares problem, which can optionally use a Z-transform and a reverse Z-transform. The model used by in block  307  may be written as: 
         [0000]        z=H   N   N+H   ξ ξ+η z ,  (28)
 
         [0000]    where z is the measurement vector, N is the integer ambiguity vector, ξ is the vector of real-valued parameters, η z , is measurement noise, Nε           q , ξε           p  with H N  and H ξ  being the corresponding design matrices and the noise assumed to be zero-mean normally distributed. The float solution after a least-squares adjustment may be given as: 
         [0000]    
       
         
           
             
               
                 
                   
                     ( 
                     
                       
                         
                           
                             N 
                             ^ 
                           
                         
                       
                       
                         
                           ξ 
                         
                       
                     
                     ) 
                   
                   , 
                   
                     
                       ( 
                       
                         
                           
                             
                               Q 
                               
                                 
                                   N 
                                   ^ 
                                 
                                  
                                 
                                   N 
                                   ^ 
                                 
                               
                             
                           
                           
                             
                               Q 
                               
                                 
                                   N 
                                   ^ 
                                 
                                  
                                 
                                   ξ 
                                   ^ 
                                 
                               
                             
                           
                         
                         
                           
                             
                               Q 
                               
                                 
                                   ξ 
                                   ^ 
                                 
                                  
                                 
                                   N 
                                   ^ 
                                 
                               
                             
                           
                           
                             
                               Q 
                               
                                 
                                   ξ 
                                   ^ 
                                 
                                  
                                 
                                   ξ 
                                   ^ 
                                 
                               
                             
                           
                         
                       
                       ) 
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   29 
                   ) 
                 
               
             
           
         
       
     
         [0060]    The BIE solution may be given by: 
         [0000]    
       
         
           
             
               
                 
                   
                       
                   
                    
                   
                     
                        
                       BIE 
                     
                     = 
                     
                       
                         
                           ∑ 
                           
                             z 
                             ∈ 
                             
                               
                                 ℤ 
                                  
                                 
                                     
                                 
                               
                               q 
                             
                           
                         
                          
                         
                             
                         
                          
                         
                           z 
                            
                           
                               
                           
                            
                           
                             w 
                              
                             
                               ( 
                               z 
                               ) 
                             
                           
                         
                       
                       
                         
                           ∑ 
                           
                             z 
                             ∈ 
                             
                               
                                 ℤ 
                                  
                                 
                                     
                                 
                               
                               q 
                             
                           
                         
                          
                         
                           w 
                            
                           
                             ( 
                             z 
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   30 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       χ 
                       z 
                       2 
                     
                      
                     
                       = 
                       def 
                     
                      
                     
                       
                          
                         
                           
                             N 
                             ^ 
                           
                           - 
                           z 
                         
                          
                       
                       
                         Q 
                         
                           
                             N 
                             ^ 
                           
                            
                           
                             N 
                             ^ 
                           
                         
                         
                           
                               
                           
                            
                           
                             - 
                             1 
                           
                         
                       
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   31 
                   ) 
                 
               
             
             
               
                 
                   
                     w 
                      
                     
                       ( 
                       z 
                       ) 
                     
                   
                   = 
                   
                     
                       exp 
                        
                       
                         ( 
                         
                           - 
                           
                             
                               
                                 ( 
                                 
                                   
                                     N 
                                     ^ 
                                   
                                   - 
                                   z 
                                 
                                 ) 
                               
                                
                               
                                 
                                   Q 
                                   
                                     
                                       N 
                                       ^ 
                                     
                                      
                                     
                                       N 
                                       ^ 
                                     
                                   
                                   
                                     
                                         
                                     
                                      
                                     
                                       - 
                                       1 
                                     
                                   
                                 
                                  
                                 
                                   ( 
                                   
                                     
                                       N 
                                       ^ 
                                     
                                     - 
                                     z 
                                   
                                   ) 
                                 
                               
                             
                             2 
                           
                         
                         ) 
                       
                     
                     = 
                     
                       
                         exp 
                         ( 
                         
                           - 
                           
                             
                               
                                  
                                 
                                   
                                     N 
                                     ^ 
                                   
                                   - 
                                   z 
                                 
                                  
                               
                               
                                 Q 
                                 
                                   
                                     N 
                                     ^ 
                                   
                                    
                                   
                                     N 
                                     ^ 
                                   
                                 
                                 
                                   
                                       
                                   
                                    
                                   
                                     - 
                                     1 
                                   
                                 
                               
                               2 
                             
                             2 
                           
                         
                         ) 
                       
                       = 
                       
                         exp 
                          
                         
                           ( 
                           
                             - 
                             
                               
                                 χ 
                                 z 
                                 2 
                               
                               2 
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   32 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       Δξ 
                       BIE 
                     
                      
                     
                       = 
                       def 
                     
                      
                     
                       
                         Q 
                         
                           
                             ξ 
                             ^ 
                           
                            
                           
                             N 
                             ^ 
                           
                         
                       
                        
                       
                         
                           Q 
                           
                             
                               N 
                               ^ 
                             
                              
                             
                               N 
                               ^ 
                             
                           
                           
                             
                                 
                             
                              
                             
                               - 
                               1 
                             
                           
                         
                          
                         
                           ( 
                           
                             
                               N 
                               ^ 
                             
                             - 
                             
                                
                               BIE 
                             
                           
                           ) 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   33 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       Δ 
                        
                       
                           
                       
                        
                       
                         Q 
                         BIE 
                       
                     
                      
                     
                       = 
                       def 
                     
                      
                     
                       
                         Q 
                         
                           
                             ξ 
                             ^ 
                           
                            
                           
                             N 
                             ^ 
                           
                         
                       
                        
                       
                         Q 
                         
                           
                             N 
                             ^ 
                           
                            
                           
                             N 
                             ^ 
                           
                         
                         
                           
                               
                           
                            
                           
                             - 
                             1 
                           
                         
                       
                        
                       
                         Q 
                         
                           
                             N 
                             ^ 
                           
                            
                           
                             ξ 
                             ^ 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   34 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       ξ 
                       BIE 
                     
                     = 
                     
                       
                         ξ 
                         ^ 
                       
                       + 
                       
                         Δξ 
                         BIE 
                       
                     
                   
                 
               
               
                 
                   ( 
                   35 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       Q 
                       
                         
                           ξ 
                           BIE 
                         
                          
                         
                           ξ 
                           BIE 
                         
                       
                     
                     = 
                     
                       
                         Q 
                         
                           
                             ξ 
                             ^ 
                           
                            
                           
                             ξ 
                             ^ 
                           
                         
                       
                       - 
                       
                         Δ 
                          
                         
                             
                         
                          
                         
                           Q 
                           BIE 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   36 
                   ) 
                 
               
             
             
               
                 
                   
                       
                   
                    
                   
                     
                       χ 
                       BIE 
                       2 
                     
                     = 
                     
                       
                         
                           
                             ∑ 
                             
                               z 
                               ∈ 
                               
                                 
                                   ℤ 
                                    
                                   
                                       
                                   
                                 
                                 q 
                               
                             
                           
                            
                           
                               
                           
                            
                           
                             
                               χ 
                               z 
                               2 
                             
                             · 
                             
                               w 
                                
                               
                                 ( 
                                 z 
                                 ) 
                               
                             
                           
                         
                         
                           
                             ∑ 
                             
                               z 
                               ∈ 
                               
                                 
                                   ℤ 
                                    
                                   
                                       
                                   
                                 
                                 q 
                               
                             
                           
                            
                           
                             w 
                              
                             
                               ( 
                               z 
                               ) 
                             
                           
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   37 
                   ) 
                 
               
             
           
         
       
     
         [0061]    The steps performed in block  307  are shown in the process  307  of  FIG. 4 . At step  504 , a tree search of the candidate narrow lane integer ambiguities  k,m ∇N NL  may be initialized, such as by computing a Z-transform, for example. The nodes in the search tree are integer vectors. During this search, at step  506 , the next candidate narrow lane integer ambiguity  k,m ∇N NL  may be visited in the tree and weighted sums Σz·w(z) of the candidate narrow lane integer ambiguities  k,m ∇N NL  and a sum of weights Σw(z) may be updated. By updating both the sum of weights Σw(z) and the weighted sum Σz·w(z) during the search, there is no need to store a large number of candidates or utilize thresholds to control the number of candidates. 
         [0062]    The weighting at step  506  may be adaptively scaled so that a large dynamic range in weight magnitudes of the candidates is avoided. For a given choice of χ 0   2 , scaled weights can be defined as 
         [0000]    
       
         
           
             
               
                 w 
                 0 
               
                
               
                 ( 
                 z 
                 ) 
               
             
              
             
               = 
               def 
             
              
             
               
                 exp 
                  
                 
                   ( 
                   
                     
                       - 
                       
                         
                           χ 
                           z 
                           2 
                         
                         2 
                       
                     
                     + 
                     
                       
                         χ 
                         0 
                         2 
                       
                       2 
                     
                   
                   ) 
                 
               
               . 
             
           
         
       
     
         [0000]    Presuming that for any χ 0   2 ε         , the ambiguity determined position solution 0.7V BIE  can be written as: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                        
                       BIE 
                     
                     = 
                     
                       
                         
                           ∑ 
                           
                               
                           
                            
                           
                             z 
                             · 
                             
                               w 
                                
                               
                                 ( 
                                 z 
                                 ) 
                               
                             
                           
                         
                         
                           ∑ 
                           
                             w 
                              
                             
                               ( 
                               z 
                               ) 
                             
                           
                         
                       
                       = 
                       
                         
                           
                             ∑ 
                             
                                 
                             
                              
                             
                               z 
                               · 
                               
                                 exp 
                                  
                                 
                                   ( 
                                   
                                     - 
                                     
                                       
                                         χ 
                                         z 
                                         2 
                                       
                                       2 
                                     
                                   
                                   ) 
                                 
                               
                             
                           
                           
                             ∑ 
                             
                                 
                             
                              
                             
                               exp 
                                
                               
                                 ( 
                                 
                                   - 
                                   
                                     
                                       χ 
                                       z 
                                       2 
                                     
                                     2 
                                   
                                 
                                 ) 
                               
                             
                           
                         
                         = 
                         
                           
                             
                               
                                 exp 
                                  
                                 
                                   ( 
                                   
                                     
                                       χ 
                                       0 
                                       2 
                                     
                                     2 
                                   
                                   ) 
                                 
                               
                                
                               
                                 ∑ 
                                 
                                     
                                 
                                  
                                 
                                   z 
                                   · 
                                   
                                     exp 
                                      
                                     
                                       ( 
                                       
                                         - 
                                         
                                           
                                             χ 
                                             z 
                                             2 
                                           
                                           2 
                                         
                                       
                                       ) 
                                     
                                   
                                 
                               
                             
                             
                               
                                 exp 
                                  
                                 
                                   ( 
                                   
                                     
                                       χ 
                                       0 
                                       2 
                                     
                                     2 
                                   
                                   ) 
                                 
                               
                                
                               
                                 ∑ 
                                 
                                     
                                 
                                  
                                 
                                   exp 
                                    
                                   
                                     ( 
                                     
                                       - 
                                       
                                         
                                           χ 
                                           z 
                                           2 
                                         
                                         2 
                                       
                                     
                                     ) 
                                   
                                 
                               
                             
                           
                           = 
                           
                             
                               ∑ 
                               
                                   
                               
                                
                               
                                 z 
                                 · 
                                 
                                   exp 
                                    
                                   
                                     ( 
                                     
                                       
                                         - 
                                         
                                           
                                             χ 
                                             z 
                                             2 
                                           
                                           2 
                                         
                                       
                                       + 
                                       
                                         
                                           χ 
                                           0 
                                           2 
                                         
                                         2 
                                       
                                     
                                     ) 
                                   
                                 
                               
                             
                             
                               ∑ 
                               
                                   
                               
                                
                               
                                 exp 
                                  
                                 
                                   ( 
                                   
                                     
                                       - 
                                       
                                         
                                           χ 
                                           z 
                                           2 
                                         
                                         2 
                                       
                                     
                                     + 
                                     
                                       
                                         χ 
                                         0 
                                         2 
                                       
                                       2 
                                     
                                   
                                   ) 
                                 
                               
                             
                           
                         
                       
                     
                   
                   , 
                 
               
               
                 
                   ( 
                   38 
                   ) 
                 
               
             
           
         
       
     
         [0000]    and therefore as: 
         [0000]    
       
         
           
             
               
                 
                   
                      
                     BIE 
                   
                   = 
                   
                     
                       
                         ∑ 
                         
                             
                         
                          
                         
                           z 
                           · 
                           
                             
                               w 
                               0 
                             
                              
                             
                               ( 
                               z 
                               ) 
                             
                           
                         
                       
                       
                         ∑ 
                         
                           
                             w 
                             0 
                           
                            
                           
                             ( 
                             z 
                             ) 
                           
                         
                       
                     
                     . 
                   
                 
               
               
                 
                   ( 
                   39 
                   ) 
                 
               
             
           
         
       
     
         [0000]    The choice of χ 0   2  may be dynamically changed during the search to be the minimum χ 0   2  of all the candidates that have been visited during the search at that point. For each candidate visited during the search, partial sums of the numerator and denominator terms of equation (29) can therefore be accumulated. 
         [0063]    In addition, because the sum of weights Σw(z) and the weighted sum Σz·w(z) are generated during the search, only the candidates with large enough weights to have a significant numerical impact may be included. For example, if z 1  represents the best solution (i.e., where X 1   2  is a minimum) found so far in a search with w(z 1 ) as its corresponding weight, and z c  is an integer vector candidate at a current node of the search with a corresponding weight w(z c ), then the candidate z c  can be considered not significant when w(z c )&lt;&lt;w(z 1 ). Candidates may be included during the search as long as: 
         [0000]        w ( z   c )≧ε w ( z   1 ),  (40)
 
         [0000]    where ε is a small threshold, such as 10 −6 .
 
This is equivalent to:
 
         [0000]    
       
         
           
             
               
                 
                   
                     χ 
                     c 
                     2 
                   
                   ≤ 
                   
                     
                       χ 
                       1 
                       2 
                     
                     + 
                     
                       ln 
                        
                       
                           
                       
                        
                       
                         
                           1 
                           
                             ɛ 
                             2 
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   41 
                   ) 
                 
               
             
           
         
       
     
         [0064]    During the search, at step  508 , determined ambiguity values  k,m ∇           NL,BIE  can be formed based on the weighted sums of the candidate narrow lane integer ambiguities and the sum of weights. It can also be determined at step  508  whether there are candidates remaining with weights greater than a predetermined threshold. If there are still candidates remaining at step  508 , then the process  307  may return to step  506  to continue the search and repeat step  506  on the next candidate. If there are no candidates remaining at step  508 , then the process  307  may continue to step  510  to finalize the search, such as by applying a reverse Z-transform, for example. The determined ambiguity values  k,m ∇           NL,BIE  may be utilized to form a constraint that can be applied to the float solution in order to calculate an ambiguity determined position solution (ΔX BIE , ΔP BIE ) at step  512 . 
         [0065]    The steps performed in step  512  are shown in the process  600  of  FIG. 5 . At step  602 , a float change vector ΔN may be calculated as the difference of the determined ambiguity values  k,m ∇           NL,BIE  (calculated at step  510 ) and the float ambiguity values  k,m ∇           NL . Accordingly, the float change vector ΔN may have elements consisting of  k,m ∇           NL,BIE − k,m ∇           NL . At step  604 , a design matrix H may be formed. The design matrix H may have the same row dimensions as the float change vector ΔN and a column dimension equal to the state size of the float state vector X. Accordingly, each row of the design matrix H may have a +1 coefficient at the state index of  k             IF  and a −1 coefficient at the state index of  m             IF , where m is the reference satellite for k. 
         [0066]    At step  606 , the Kalman gain K may be computed based on the design matrix H and the covariance matrix P of the Kalman filter, as given by: 
         [0000]        K=PH   T ( HPH   T ) −1 .  (42)
 
         [0000]    Correction terms may be formed at steps  608  and  610 . In particular, at step  608 , a state correction term ΔX BIE =KΔN, and at step  610  a covariance correction term ΔP BIE =KHP. The state correction term ΔX BIE  and the covariance correction term ΔP BIE  may form the ambiguity determined position solution. The state vector can be corrected by adding the state correction term ΔX BIE (i.e., X+ΔX BIE ), and the covariance matrix can be corrected by subtracting the covariance correction term ΔP BIE  (i.e., P−ΔP BIE ). 
         [0067]    Returning to  FIG. 2 , the determined ambiguity values  k,m ∇           NL,BIE  may be time-smoothed at block  308  to generate smoothed ambiguity values  k,m ∇           NL,SBIE  The steps performed in block  308  are shown in the process  308  of  FIG. 6 . At step  702 , a smoothed BIE estimator may be updated with the determined ambiguity values  k,m ∇           NL,BIE  to generate the smoothed ambiguity values  k,m ∇           NL,SBIE . The smoothed ambiguity values  k,m ∇           NL,SBIE  are generally smoothly varying in time, changing gradually from the estimated float narrow lane ambiguities  k,m ∇           NL  to the determined ambiguity values  k,m ∇           NL,BIE . 
         [0068]    The steps performed at step  702  are shown in the process  702  of  FIG. 7 . The following terms and notation are used in the discussed on the process  702 . In particular,  k             NL  represents the float narrow lane value for a satellite k in cycles;  k             NL,BIE  represents the ambiguity determined narrow lane value for a satellite k in cycles;  k             NL,SBIE  represents the smoothed ambiguity determined narrow lane value for a satellite k in cycles;  k              NL,BIAS  represents the narrow lane float bias value for a satellite k in cycles; W represents the window length for smoothing determined narrow lane values and is given as an integer number of intervals (e.g., epochs); and  k C represents a smoothing count for satellite k, is given as an integer number of intervals (e.g., epochs), is initialized to 0, and is incremented by one each interval. In addition, t may denote the current interval and t−1 may denote the previous interval. 
         [0069]    At step  802 , unavailable satellites may be removed from being used in the time-domain smoothing. An unavailable satellite may include satellites that the receiver can no longer receive signals from. At step  804 , it may be determined whether the reference satellite has changed from a previous interval. If the reference satellite has changed at step  804 , then the process  702  continues to step  814  to calculate time-domain smoothed ambiguity values for the new reference satellite and satellites other than the removed unavailable satellites. The old reference satellite may be denoted as m 0  and the new reference satellite may be denoted as m 1 . In addition, the difference in integer ambiguities between the new reference satellite m 1  and a given satellite k may be given by: 
         [0000]        k,m     1     ∇N   NL = k,m     0     ∇N   NL − m     1     m     0     ∇N   NL   (43)
 
         [0070]    At step  814 , if the new reference satellite was not used in a previous interval, then the time-domain smoothed ambiguity values for the old reference satellite m 0 , the new reference satellite m 1 , and a given satellite k may be calculated by: 
         [0000]        k             NL,SMIE ( t− 1)= k             NL,SBIE ( t− 1)− m     1     ,m     0   ∇           NL ( t ); k≠m   0   ,k≠m   1   (44)
 
         [0000]        k   C ( t− 1)=0  (45)
 
         [0000]        m     0               NL,SBIE ( t− 1)= m     1     ,m     0   ∇           NL ( t )  (46)
 
         [0000]        m     0     C ( t− 1)=0  (47)
 
         [0000]        m     1               NL,SBIE ( t− 1)=0  (48)
 
         [0000]        m     1     C ( t− 1)=0.  (49)
 
         [0071]    If the new reference satellite m 1  was used in a previous interval, then the time-domain smoothed ambiguity values for the old reference satellite m 0 , the new reference satellite m 1 , and a given satellite k may be adjusted by: 
         [0000]        k             NL,SBIE ( t− 1)= k             NL,SBIE ( t− 1)− m     1               NL,SBIE ( t− 1); k≠m   0   ,k≠m   1   (50)
 
         [0000]        k   C ( t− 1)=min( k   C ( t− 1), m     1     C ( t− 1))  (51)
 
         [0000]        m     0     C ( t− 1)=min( m     0     C ( t− 1), m     1     C ( t− 1))  (52)
 
         [0000]        m     1               NL,SBIE ( t− 1)=0  (53)
 
         [0000]        m     1     C ( t− 1) is unchanged.  (54)
 
         [0072]    Following step  814  or if the reference satellite has not changed at step  804 , then the process  702  continues to step  806 . At step  806 , it may be determined whether new satellites should be used in the time-domain smoothing. New satellites may include the satellites which were not used in utilized in a previous interval. If it is determined that no new satellites should be used at step  806 , then the process  702  continues to step  812 , as described below. However, if it is determined that new satellites should be used at step  806 , then the process  702  continues to step  808 . At step  808 , time-domain smoothed ambiguity values  k,m ∇           NL,SBIE  may be calculated for the new satellites k having a reference satellite m (where k≠m), as given by: 
         [0000]        k             NL,SBIE ( t− 1)= k,m ∇           NL,BIE ( t )  (55)
 
         [0000]        k   C ( t− 1)=0.  (56)
 
         [0073]    Following step  808 , the time-domain smoothed ambiguity values for the new satellites may be adjusted at step  810  by a bias between the estimated float narrow lane ambiguities  k,m ∇           NL , and time-domain smoothed ambiguity values  k,m ∇           NL,SBIE  from a previous interval. By adjusting with the bias, the initial bias may be minimized between the estimated float narrow lane ambiguities  k,m ∇           NL  and time-domain smoothed ambiguity values  k,m ∇           NL,SBIE  from the previous interval. This may be calculated by: 
         [0000]        k             NL,SBIE ( t− 1)= k             NL,SBIE ( t− 1)+ k,m     1   ∇Δ           NL,BIAS ( t− 1),  (57)
 
         [0000]      where  k,j ∇Δ           NL,BIAS ( t )           k Δ           NL,BIAS ( t )− j Δ           NL,BIAS ( t )  (58)
 
         [0000]    The term  k Δ           NL,BIAS (t) in equation (58) represents the respective undifferenced float ambiguity entries in the state vector increment ΔX SBIE  of the ambiguity determined position solution from the prior interval. 
         [0074]    At step  812 , the time-domain smoothed ambiguity values for all of the satellites may be updated based on the determined ambiguity values  k,m ∇           NL,BIE . Step  812  may be performed following step  810  or if no new satellites were determined to be needed at step  806 . In some embodiments, an exponential filter such as a recursive estimator may be utilized but in other embodiments, other techniques for smoothing may be utilized. The update of the time-domain smoothed ambiguity values  k,m ∇           NL,SBIE  may be performed according to the following for each satellite k: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                        
                       
                         NL 
                         , 
                         SBIE 
                       
                         
                         
                         
                       k 
                     
                      
                     
                       ( 
                       t 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                           W 
                           - 
                           1 
                         
                         W 
                       
                        
                       
                         ( 
                         
                           
                              
                             
                               NL 
                               , 
                               SBIE 
                             
                               
                               
                               
                             k 
                           
                            
                           
                             ( 
                             
                               t 
                               - 
                               1 
                             
                             ) 
                           
                         
                         ) 
                       
                     
                     + 
                     
                       
                         1 
                         W 
                       
                        
                       
                         ( 
                         
                           
                              
                             
                               NL 
                               , 
                               BIE 
                             
                               
                               
                               
                             k 
                           
                            
                           
                             ( 
                             t 
                             ) 
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   ( 
                   59 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       
                           
                         k 
                       
                        
                       C 
                     
                      
                     
                       ( 
                       t 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                             
                           k 
                         
                          
                         C 
                       
                        
                       
                         ( 
                         
                           t 
                           - 
                           1 
                         
                         ) 
                       
                     
                     + 
                     1 
                   
                 
               
               
                 
                   ( 
                   60 
                   ) 
                 
               
             
             
               
                 
                   
                     
                       σ 
                       
                         NL 
                         , 
                         SBIE 
                       
                       2 
                         
                         
                       k 
                     
                      
                     
                       ( 
                       t 
                       ) 
                     
                   
                   = 
                   
                     
                       1 
                       
                         
                           
                               
                             k 
                           
                            
                           C 
                         
                          
                         
                           ( 
                           t 
                           ) 
                         
                       
                     
                      
                     
                       
                         
                           ( 
                           
                             
                               
                                  
                                 
                                   NL 
                                   , 
                                   SBIE 
                                 
                                   
                                   
                                   
                                 k 
                               
                                
                               
                                 ( 
                                 t 
                                 ) 
                               
                             
                             - 
                             
                               
                                  
                                 
                                   NL 
                                   , 
                                   BIE 
                                 
                                   
                                   
                                   
                                 k 
                               
                                
                               
                                 ( 
                                 t 
                                 ) 
                               
                             
                           
                           ) 
                         
                         2 
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   61 
                   ) 
                 
               
             
           
         
       
     
         [0075]    Following step  702 , an ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) may be calculated at step  704  based on the smoothed ambiguity values  k,m ∇           NL,SBIE . Step  704  may include the steps described above with relation to the process  600  of  FIG. 5 , and may also correspond to block  309  of  FIG. 2 . One difference is at step  606  where the Kalman gain K is computed. In the case of using the smoothed ambiguity values  k,m ∇           NL,SBIE , the Kalman gain K includes uncertainty in the constraint, and is given by: 
         [0000]        K=PH   T ( HPH   T   +R ) −1   (62)
 
         [0000]    where the covariance matrix R=diag(var(           BIE )) is computed as a function of the weights used in the search performed at block  307  of  FIG. 2 . At step  706 , float ambiguity biases may be stored that are derived from the ambiguity determined position solution (ΔX SBIE , ΔP SBIE ) from step  704 . The float ambiguity biases are used as described above with relation to step  810  of the process  702  of  FIG. 7 . 
         [0076]    Returning to  FIG. 2 , following block  309 , a final position solution may be computed using block  350 . Block  350  may include several steps, as shown in the process  350  of  FIG. 8 . It should be noted that the weighting process described below may also be applied to the ambiguity determined position solution (ΔX BIE , ΔP BIE ) from block  307 , in some embodiments. At step  902  (embodied in block  310 ), the ambiguity determined position solution (ΔX SBIE , ΔP SBIE ) from block  309  may be blended with the float solution from block  304  to form a weighted smoothed ambiguity determined position solution (ΔX WSBIE , ΔP WSBIE ). In particular, a time-weighted weighting factor IT may be determined, where 0≦ w ≦1. The weighted smoothed ambiguity determined position solution (ΔX WSBIE , ΔP WSBIE ) may be computed from the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) as follows: 
         [0000]      Δ X   WSBIE   = w ΔX   SBIE   (63)
 
         [0000]      Δ P   WSBIE   = w     2   ΔP   SBIE   (64)
 
         [0077]    As can be seen from equations (53) and (54), when  w =0, the float solution will be utilized and when  w =1, the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) will he utilized. A larger weighting factor  w  is used when the confidence in the correctness of the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) increases. 
         [0078]    The weighting may be determined based on one or more factors. The factors may include a minimum variance solution, an acceptability factor α that indicates whether the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) is acceptable, a convergence indicator of the float solution, and/or a look-up table that is indexed by ranges of corresponding error variances and corresponding figures of merit of the float solution and the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ). 
         [0079]    The minimum variance solution factor may minimize an error variance g of a combination of the float solution and the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ). The error variances may be three dimensional position variances that are computed as a trace of the corresponding three dimensional position covariance matrix. The error variance of the float solution may be denoted as σ FLOAT   2  and the error variance of the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) may be denoted as σ SBIE   2 . The minimum variance weighting between the float solution and the ambiguity determined position solution (ΔX SBIE , ΔP SBIE ) may be given by: 
         [0000]      Δ X=gΔX   SBIE ,  (65)
 
         [0000]    where 
         [0000]    
       
         
           
             
               
                 
                   g 
                    
                   
                     = 
                     Δ 
                   
                    
                   
                       
                   
                    
                   
                     
                       
                         
                           0 
                           / 
                           
                             σ 
                             FLOAT 
                             2 
                           
                         
                         + 
                         
                           1 
                           / 
                           
                             σ 
                             SBIE 
                             2 
                           
                         
                       
                       
                         
                           1 
                           / 
                           
                             σ 
                             FLOAT 
                             2 
                           
                         
                         + 
                         
                           1 
                           / 
                           
                             σ 
                             SBIE 
                             2 
                           
                         
                       
                     
                     = 
                     
                       
                         
                           σ 
                           FLOAT 
                           2 
                         
                         
                           
                             σ 
                             SBIE 
                             2 
                           
                           + 
                           
                             σ 
                             FLOAT 
                             2 
                           
                         
                       
                       = 
                       
                         
                           1 
                           
                             
                               ( 
                               
                                 
                                   σ 
                                   SBIE 
                                   2 
                                 
                                 / 
                                 
                                   σ 
                                   FLOAT 
                                   2 
                                 
                               
                               ) 
                             
                             + 
                             1 
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   66 
                   ) 
                 
               
             
           
         
       
     
         [0000]    The error variance g will be limited according to 0≦g≦1, and will have a small value when the uncertainty of the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) is large, and conversely will have a large value when the uncertainty of the float solution is relatively large. In some embodiments, a time-filtered version g of the error variance  g  may be utilized due to noise. 
         [0080]    The acceptability factor α may indicate whether the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) is acceptable and be limited according to 0≦α≦1. The acceptability factor α may be smaller when the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) is unacceptable and may be larger when it is acceptable. The acceptability factor α may be dependent on a quality of the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ), such as based on the magnitude of the quadratic form of the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) and the ratio test. The ratio test can be defined as the ratio of the quadratic form of the second best solution to the best solution. When the ratio is large, it can indicate that the best solution is the correct solution. 
         [0081]    The quadratic form of the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) may be very large if there are problems with determining ambiguity values. In this case, the acceptability factor α may be 0. However, if the ratio test is large and the quadratic form is small, then the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) may be deemed more reliable and have an acceptability factor α of 1. In between these cases, the acceptability factor α may by found by a smooth function which decreases from 1 to 0 as the ratio decreases and the quadratic form increases. 
         [0082]    The convergence indicator factor may indicate whether the float solution has reached a steady state. In this case, the convergence indicator may be greater when the float solution has reached a steady state or is close to a steady state, and conversely may be less in other situations. 
         [0083]    In embodiments, the weighting factor  w  may be chosen based on a combination of these factors. For example, if the acceptability factor α&lt;1 or the float solution has converged and the ambiguity determined position solution (ΔX SBIE ,ΔP SBIE ) is determined to be not trustworthy, then the weighting factor  w  may be equal to a g . In other cases, the weighting factor  w  may be equal to the acceptability factor α. 
         [0084]    Returning to the process  350  of  FIG. 8 , it may be determined at step  904  if the weighted smoothed ambiguity determined position solution (ΔX WSBIE ,ΔP WSBIE ) has a poor quality. If the weighted smoothed ambiguity determined position solution (ΔX WSBIE ,ΔP WSBIE ) has a poor quality, then the process  350  may continue to step  912  to transition to the float solution. However, if the weighted smoothed ambiguity determined position solution (ΔX WSBIE ,ΔP WSBIE ) does not have a poor quality then the process  350  may continue to step  906  (also embodied by block  311 ). 
         [0085]    At step  906 , an estimated position jump in the weighted smoothed ambiguity determined position solution (ΔX WSBIE ,ΔP WSBIE ) may be calculated. The estimated position jump may be calculated by comparing the change in position since the last interval to the change according to an estimate using the integrated carrier phase that infers positions between the intervals using carrier phase time differences, as is known in the art. At step  908 , it may be determined whether the estimated position jump from step  906  is greater than a predetermined threshold. If the estimated position jump is greater than the threshold, then the process  350  may continue to step  912  to transition to the float solution. However, if the estimated position jump is not greater than the threshold, then the process  350  may continue to step  910 . 
         [0086]    At step  910 , the navigation output showing the position of the receiver may be updated at block  312  by adjusting the float solution from block  304  with the weighted smoothed ambiguity determined position solution (ΔX WSBIE ,ΔP WSBIE ). The state vector can be corrected by adding the state correction term ΔX WSBIE  (i.e., X+ΔX WSBIE ), and the covariance matrix can be corrected by subtracting the covariance correction term ΔP WSBIE  (i.e., P−ΔP WSBIE ). 
         [0087]    Transitioning to the float solution at step  912  may be performed to ensure that the transition is relatively smooth and not too rapid. At each interval, an offset ΔX WSBIE  is typically stored. The offset can be denoted as ΔX t  at a subsequent interval when the decision to transition to the float solution is made. Accordingly, at a time t, the transition may be started and ΔX t           ΔX WSBIE (t−1). Starting at the interval when the transition is begun, the position change may be limited at each interval due to the offset ΔX t . In particular, the position change may be limited to not vary by more than a magnitude E that is a small predetermined value. The transition may be made over N steps, where N is the rounded-off value of ∥ΔX t ∥/ε. Accordingly, at step  401 , for k=1 . . . N, the float solution can be transitioned to by the end of the transition period, where 
         [0000]    
       
         
           
             
               Δ 
                
               
                   
               
                
               
                 
                   X 
                   WSBIE 
                 
                  
                 
                   ( 
                   
                     t 
                     + 
                     k 
                   
                   ) 
                 
               
             
             = 
             
               
                 
                   N 
                   - 
                   K 
                 
                 N 
               
                
               Δ 
                
               
                   
               
                
               
                 
                   X 
                   t 
                 
                 . 
               
             
           
         
       
     
         [0088]    Any process descriptions or blocks in figures should be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps in the process, and alternate implementations are included within the scope of the embodiments of the invention in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those having ordinary skill in the art. 
         [0089]    This disclosure is intended to explain how to fashion and use various embodiments in accordance with the technology rather than to limit the true, intended, and fair scope and spirit thereof. The foregoing description is not intended to be exhaustive or to be limited to the precise forms disclosed. Modifications or variations are possible in light of the above teachings. The embodiment(s) were chosen and described to provide the best illustration of the principle of the described technology and its practical application, and to enable one of ordinary skill in the art to utilize the technology in various embodiments and with various modifications as are suited to the particular use contemplated. All such modifications and variations are within the scope of the embodiments as determined by the appended claims, as may be amended during the pendency of this application for patent, and all equivalents thereof, when interpreted in accordance with the breadth to which they are fairly, legally and equitably entitled.