# Patent

## Title
Common coordinate-quartz loop for reducing the impact of shock and vibration on global navigation satellite system measurements

## Abstract

A navigation receiver operably coupled to an antenna can determine location by receiving and processing radiofrequency navigation signals from global navigation satellites.
The navigation receiver includes a quartz crystal oscillator that serves as a reference frequency source for a local oscillator signal that is mixed with the radiofrequency navigation signals to generate intermediate signals at intermediate frequencies lower than the radiofrequencies.
In particular applications, the quartz crystal oscillator and the antenna are subjected to vibration and shock.
A shift in the frequency or phase of the intermediate signals can result, and performance of phase-lock loops can degrade.
Performance is improved with a phase-lock loop that processes a combination of individual channel control signals and common control signals.
The common control signals are generated by a common coordinate-quartz loop discriminator that processes signals generated from the entire group of navigation signals received by the navigation receiver.

## Background
BACKGROUND OF THE INVENTION
The present invention relates generally to global navigation satellite systems, and more particularly to method and apparatus for reducing the impact of shock and vibration on global navigation satellite system measurements.
Global navigation satellite systems (GNSSs) can determine locations with high accuracy. Currently deployed global navigation satellite systems are the United States Global Positioning System (GPS) and the Russian GLONASS. Other global navigation satellite systems, such as the European GALILEO system, are under development. In a GNSS, a navigation receiver receives and processes radio signals transmitted by satellites located within a line-of-sight of the receiver. The satellite signals comprise carrier signals modulated by pseudo-random binary codes. The receiver measures the time delays of the received signals relative to a local reference clock or oscillator. Code measurements enable the receiver to determine the pseudo-ranges between the receiver and the satellites. The pseudo-ranges differ from the actual ranges (distances) between the receiver and the satellites due to various error sources and due to variations in the time scales of the satellites and the receiver. If signals are received from a sufficiently large number of satellites, then the measured pseudo-ranges can be processed to determine the code coordinates and coordinate time scales at the receiver. This operational mode is referred to as a stand-alone mode, since the measurements are determined by a single receiver. A stand-alone system typically provides meter-level accuracy.
To improve the accuracy, precision, stability, and reliability of measurements, differential navigation (DN) systems have been developed. In a DN system, the position of a user is determined relative to a base station, also referred to as a base. The base is typically fixed, and the coordinates of the base are precisely known; for example, by surveying. The base contains a navigation receiver that receives satellite signals and that can determine the coordinates of the base by GNSS measurements.
The user, whose position is to be determined, can be stationary or mobile; in a DN system, the user is often referred to as a rover. The rover also contains a navigation receiver that receives satellite signals. Signal measurements processed at the base are transmitted to the rover via a communications link. To accommodate a mobile rover, the communications link is often a wireless link. The rover processes the measurements received from the base, along with measurements taken with its own receiver, to improve the accuracy of determining its position. Accuracy is improved in the differential navigation mode because errors incurred by the receiver at the rover and by the receiver at the base are highly correlated. Since the coordinates of the base are accurately known, measurements from the base can be used to compensate for the errors at the rover. A differential global positioning system (DGPS) computes locations based on pseudo-ranges only.
The location determination accuracy of a differential navigation system can be further improved by supplementing the code pseudo-range measurements with measurements of the phases of the satellite carrier signals. If the carrier phases of the signals transmitted by the same satellite are measured by both the navigation receiver in the base and the navigation receiver in the rover, processing the two sets of carrier phase measurements can yield a location determination accuracy to within several percent of the carrier's wavelength. A differential navigation system that computes locations based on real-time carrier signals, in addition to the code pseudo-ranges, is often referred to as a real-time kinematic (RTK) system. Processing carrier phase measurements to determine coordinates includes the step of ambiguity resolution; that is, determining the integer number of cycles in the carrier signal received by the navigation receiver from an individual satellite.
Global navigation satellite systems have been used for control of construction vehicles, such as bulldozers: an antenna and a receiver are mounted on the vehicle. During operation of the construction vehicle, the antenna and the receiver can be subjected to shock and vibration, which can disturb the GNSS measurements and decrease the accuracy of the coordinate measurements. Method and apparatus for reducing the impact of shock and vibration on GNSS measurements would be advantageous.
BRIEF SUMMARY OF THE INVENTION
Disclosed are a method and an apparatus for improving the performance of a navigation receiver operably coupled to an antenna, in which the navigation receiver and the antenna are subjected to vibration and shock.
The navigation receiver receives a group of navigation signals. Each individual navigation signal in the group of navigation signals corresponds to an individual global navigation satellite in a group of global navigation satellites, and each individual navigation signal in the group of navigation signals has a corresponding radiofrequency. A local oscillator generates a local oscillator signal having a local oscillator frequency; the local oscillator signal is based on the quartz crystal oscillator signal, and the local oscillator frequency is based on the quartz crystal oscillator frequency. A group of intermediate signals based on the group of navigation signals and based on the local oscillator signal is generated. Each individual intermediate signal in the group of intermediate signals corresponds to an individual navigation signal in the group of navigation signals, and each individual intermediate signal in the group of intermediate signals has a corresponding intermediate frequency based on a difference between the corresponding radiofrequency and the local oscillator frequency.
A group of sequences of digital signals is generated. Each individual sequence of digital signals in the group of sequences of digital signals corresponds to a digital representation of an individual intermediate signal in the group of intermediate signals. For each individual sequence of digital signals in the group of sequences of digital signals, a phase of the individual sequence of digital signals is tracked with an individual channel phase-lock loop. Tracking is performed by the following steps.
An individual channel correlator, operating at a first clock frequency, generates a sequence of in-phase correlation signals and a sequence of quadrature-phase correlation signals; the sequence of in-phase correlation signals and the sequence of quadrature-phase correlation signals are based on the sequence of individual digital signals and based on a sequence of individual channel reference signals received from an individual channel numerically-controlled oscillator. A common coordinate-quartz loop discriminator, operating at a second clock frequency, generates a sequence of common control signals. An individual channel accumulator, operating at a third clock frequency, generates a sequence of accumulated in-phase correlation signals and a sequence of accumulated quadrature-phase correlation signals; the sequence of accumulated in-phase correlation signals is based on the sequence of in-phase correlation signals, and the sequence of accumulated quadrature-phase correlation signals is based on the sequence of quadrature-phase correlation signals.
An individual channel discriminator, operating at the third clock frequency, generates a sequence of phase-lock loop tracking error signals based on the sequence of accumulated in-phase correlation signals and based on the sequence of accumulated quadrature-phase correlation signals. An individual channel loop filter, operating at the third clock frequency, generates a sequence of individual channel control signals. The individual channel numerically-controlled oscillator, operating at a fourth clock frequency, generates the sequence of individual channel reference signals; the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual channel control signals. The individual channel numerically-controlled oscillator transmits the sequence of individual channel reference signals to the individual channel correlator.
The sequence of common control signals can be (a) a sequence of common phase-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the group of sequences of digital signals; (b) a sequence of common frequency-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the group of sequences of digital signals; or (c) a sequence of common phase-control signals and a sequence of common frequency-control signals, in which the sequence of common phase-control signals and the sequence of common frequency-control signals are based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the group of sequences of digital signals.
The sequence of individual channel control signals can be (a) a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals; (b) a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, in which the sequence of individual channel frequency-control signals and the sequence of individual channel phase-control signals are based on the sequence of phase-lock loop tracking error signals; (c) a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals; or (d) a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, in which the sequence of individual channel frequency-control signals is based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals and in which the sequence of individual channel phase-control signals is based on the sequence of phase-lock loop tracking error signals.
These and other advantages of the invention will be apparent to those of ordinary skill in the art by reference to the following detailed description and the accompanying drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows an example of a global navigation satellite system;
FIG. 2 shows a schematic functional block diagram of a navigation receiver;
FIG. 3A and FIG. 3B show schematic functional block diagrams of a digital channel processor;
FIG. 4A shows a schematic functional block diagram of a first embodiment of a discriminator for a common coordinate-quartz loop;
FIG. 4B shows a schematic functional block diagram of a second embodiment of a discriminator for a common coordinate-quartz loop;
FIG. 5A shows a schematic functional block diagram of a first embodiment of a loop filter for a phase-lock loop;
FIG. 5B shows a schematic functional block diagram of a second embodiment of a loop filter for a phase-lock loop;
FIG. 6 shows a schematic functional block diagram of a loop filter for an embodiment of a common coordinate-quartz loop;
FIG. 7 shows a schematic functional block diagram of an analog radiofrequency processing unit;
FIG. 8 shows a functional block diagram of a numerically-controlled oscillator; and
FIG. 9 shows plots of frequency shifts as a function of vibration frequency.



DETAILED DESCRIPTION
FIG. 1 shows a high-level schematic of a global navigation satellite system (GNSS). The GNSS is used to determine, for example, the three-dimensional coordinates of the rover 102 travelling along the ground surface 130; in this example, the rover is a construction vehicle—in particular, a bulldozer. A constellation of global navigation satellites 120 transmits navigation signals. Shown are five representative navigation satellites, denoted navigation satellite 120A, navigation satellite 120B, navigation satellite 120C, navigation satellite 120D, and navigation satellite 120E.
Refer to the rover 102. In this example, the antenna 110 and the navigation receiver 116 are mounted on the blade 106 via the mast 112; therefore, the antenna 110 and the navigation receiver 116 can be used to determine the three-dimensional coordinates of the blade 106. Navigation signals received by the antenna 110 are processed by the navigation receiver 116, which transmits measurements via a signal cable (not shown) or via a wireless link to a control unit 114, typically mounted inside the cab 104. In the example shown, navigation signal 121A, navigation signal 121B, navigation signal 121C, navigation signal 121D, and navigation signal 121E are received by the antenna 110 from navigation satellite 120A—navigation satellite 120E, respectively.
Both the antenna and the navigation receiver can be subjected to shock and vibration. The engine powering the rover 102 (typically a heavy-duty diesel engine for a bulldozer) can be a source of vibration. Shock and vibration can arise from movement of the rover 102 over rough terrain, and earthmoving operations with the blade 106 can subject the antenna 110 and the navigation receiver 116 to severe shock and vibration. Shock and vibration can degrade GNSS measurements; details are discussed below. Embodiments of the invention described herein reduce the impact of shock and vibration on GNSS measurements.
Each satellite in a global navigation satellite system can transmit signals on one or more frequency bands (for example, on the L1 and L2 frequency bands). A signal transmitted by a specific satellite on a specific frequency band corresponds to a specific satellite channel. Satellite channels are processed in a navigation receiver as described below. To simplify the discussion, processing of signals for a single frequency band is described; however, embodiments of the invention can be implemented for multiple frequency bands.
FIG. 2 shows a high-level schematic functional block diagram of the navigation receiver 116. The received signal 201 represents the total signal (also referred to as the combined signal, aggregate signal, or composite signal) received by the antenna 110. For the example shown in FIG. 1, the received signal 201 includes the navigation signal 121A—navigation signal 121E. The received signal 201 is first inputted into the analog radiofrequency (RF) processing unit 202, details of which are described below. The output analog signal 203 is inputted into the analog-digital converter (ADC) 204, which digitizes the analog signal 203. The output digital signal 205 is then inputted into the digital channel processing unit 206, which processes navigation data such as phase measurements, code measurements, and information symbols used for solving navigation tasks. The output digital signal 207 is inputted into the user data processing unit 208, which processes user data such as coordinates and velocity of the rover. The output digital signal 213 is inputted, for example, into the control system 114 (FIG. 1). The control system 114, which can also receive signals from other sensors mounted on the dozer body or the blade 106, can be used to control the elevation and orientation of the blade 106. The user data processing unit 208 also calculates the output 209 which is fed back to the digital channel processing unit 206. The output 209 is described below.
Analog RF processing units of varying complexity are well-known in the art. FIG. 7 shows a high-level functional schematic of a simplified example of the analog RF processing unit 202. The received signal 201 (FIG. 2) is inputted into the low-noise amplifier (LNA) 702. The amplified output signal 701 is inputted into the RF bandpass filter 704. The filtered RF signal 703, with frequency ƒRF, is inputted into the RF mixer 706. The quartz crystal oscillator (QCO) 708 generates the QCO reference signal 705 with frequency ƒq. The QCO reference signal 705 is inputted into the local oscillator 710, which generates the local oscillator (LO) signal 707 with frequency ƒLO. The value and stability of the frequency ƒLO depends on the value and stability of the frequency ƒq. If the frequency ƒq, varies, then there is a correlated variation in the frequency ƒLO. The local oscillator signal 707 is inputted into the RF mixer 706.
The output signal 709 from the RF mixer 706 comprises signals at the two heterodyne frequencies (ƒRF+ƒLO) and (ƒRF−ƒLO). The output signal 709 is inputted into the intermediate frequency (IF) bandpass filter 712, which rejects the signal at the upconverted frequency (ƒRF+ƒLO) and passes the signal at the downconverted frequency (ƒRF−ƒLO). The output signal 203 is the downconverted analog signal with an intermediate frequency (IF) ƒIF=(ƒRF−ƒLO) that is lower than the RF frequency of the received signal 201. The output signal 203 is transmitted to the analog-digital converter 204 (FIG. 2) for further processing. The QCO 708 also transmits the QCO reference signal 211 to the digital channel processing unit 206 and the analog-digital converter 204 (FIG. 2). The QCO reference signal 211 serves as a precision timing (clock) signal.
In some GNSSs, such as GPS, within a single frequency band (such as L1), all the satellites transmit on the same carrier frequency (the value of ƒRF is the same for all carrier signals). The navigation signal from a specific satellite can be identified by a unique pseudo-random number code. In other GNSSs, such as GLONASS, within a single frequency band (such as L1), each satellite transmits on a different carrier frequency (the value of ƒRF is different for each satellite). The navigation signal from a specific satellite can be identified by a unique carrier frequency. To simplify the discussion, in the description below, all the carrier signals have the same carrier frequency. With proper transformation of variables, however, the embodiments described below can be adapted for a GNSS in which each carrier signal has a different carrier frequency.
A set of functional steps can be grouped into user-specified functional blocks and functional sub-blocks. Herein, a module refers to a user-specified set of hardware, firmware, and software that implements (executes) a user-specified set of functional steps. The same set of functional steps can be implemented by different user-specified modules; that is, different combinations of hardware, firmware, and software.
Processing steps that require high-speed execution are typically performed by hard-wired circuits; for example, application-specific integrated circuits (ASICs) and field-programmable gate arrays (FPGAs) and firmware.
In processing each satellite channel, two tracking systems track changes in parameters of the received signal. The first tracking system, referred to as the delay-lock loop (DLL), tracks changes in the pseudo-random number (PRN) modulation code delay of the received signal. The second tracking system, referred to as the phase-lock loop (PLL), tracks changes in the carrier phase of the received signal.
In a tracking system, the input signal is multiplied by different reference signals generated in the navigation receiver, and the multiplication products are accumulated. Devices that perform multiplication are referred to as multipliers. Devices that perform several combinations of multiplications are referred to as combined multipliers. Devices that perform accumulation are referred to as accumulators or counters. Devices that perform multiplication followed by accumulation are referred to as correlators. The output of a correlator is the cross-correlation function of the input signal and the reference signal; for simplicity, the output is referred to as the correlation signal.
A specific satellite channel is typically processed via three parallel processing paths. The first path calculates the in-phase correlation signal I. This signal is obtained when the reference signal in the correlator is in-phase relative to the incoming signal carrier. The reference code in the reference signal is a replica of the PRN code modulating the incoming signal. The output of this correlator (correlator-I) is the in-phase correlation signal I.
The phase of the reference signal can differ from that of the carrier phase of the incoming signal; similarly, the reference code delay can differ from the modulation code delay of the incoming signal. If the phase shift between the phase of the reference carrier and the phase of the incoming carrier is denoted as φ, and the time shift between the code delay time of the reference code and the code delay time of the incoming signal is denoted as τ, then the in-phase correlation signal I can be approximated by:

I=k·US·μ·Ro(τ)·cos(φ)+Iin,  (E1)

where:

    
    
        Ro(τ) is the normalized cross-correlation function of the incoming PRN code after filtering in the analog RF processing unit and a reference code which is a locally generated copy of the PRN code modulating the incoming signal;
        US is the amplitude of the incoming signal;
        μ=±1 is a binary information symbol modulating the incoming signal;
        k is a proportionality coefficient; and
        Iin is the interference at the output of correlator-I caused by additive interference at the receiver input.

[To simplify the notation, the satellite channel index (see below) is omitted in (E1)-(E5).]

    
    


The second correlator (correlator-Q) calculates the quadrature correlation signal Q. This signal is obtained when the second reference carrier is shifted by τ/2 from the first reference carrier. The reference code is the same as the one in the first correlator (correlator-I). The quadrature correlation signal Q formed at the output of correlator-Q can be approximated by:

Q=k·US·μ·Ro(τ)·sin(φ)+Qin,  (E2)

where Qin is the interference at the output of correlator-Q caused by additive interference at the receiver input.

The third correlator (correlator-dI) calculates the code correlation signal dI that is used to generate an error signal to control the system tracking modulation code delays (DLL). To obtain this signal, the first reference carrier, in-phase with the input carrier, and a reference code that consists of a sequence of short strobe pulses (strobe sequence) are used. Strobe pulses correspond to the changes of sign of input PRN-code elements (chips). Note that the sign of a strobe pulse matches the sign of a chip which follows the strobe pulse. The code correlation signal dI can be approximated by:

dI=k·US·μ·ΔRo(τ)·cos(φ)+dIin,  (E3)

where ΔRo(τ) is the cross-correlation function of the incoming PRN code after filtering in the analog RF processing unit and the reference code generated in the form of a sequence of short strobe pulses, and dIin is the interference at the output of correlator-dI caused by the additive interference at the receiver input.

During receiver operation, values in expressions (E1)-(E3) change over time, and the corresponding correlation signals change respectively. Correlation signals I, Q, and dI are used for coordinating (co-working) the PLL and DLL tracking systems.
Expressions (E1)-(E3) hold under the condition that PLL errors φ and DLL errors τ remain practically the same over the time (duration) of accumulating in correlators. In tracking mode, this condition is true to an adequate accuracy. In addition, in tracking mode, values φ and τ are small, and Ro(τ)·cos(φ)≅1. From (E1), it follows that the in-phase correlation signal I produces a binary information symbol sequence μ=±1. These symbols convey messages from the navigation satellites to the navigation receiver. The messages include satellite data on satellite coordinates, expected parameters of signal propagation, and other information that is used for coordinate determination.
The PLL and DLL digital tracking systems are implemented by closed circuits comprising four main components: correlator, discriminator, loop filter, and numerically-controlled oscillator (NCO). Depending on the purpose and type of the tracking system (PLL or DLL), these tracking systems use the I, Q, or dI correlation signals in different combinations. The correlators include combined multipliers that multiply the incoming signal by corresponding reference signals and accumulators that accumulate corresponding multiplication products over a specified time interval. Combined multipliers multiply the incoming signal by an in-phase reference carrier and reference code, by a quadrature-phase reference carrier and reference code, and by an in-phase reference carrier and reference strobe.
A typical PLL discriminator forms an error signal Zd according to the algorithm:

Zd=arctan(Q/I).  (E4)

From (E1), (E2), and (E4), error signal Zd is proportional to error φ within ±π/2, and then it periodically repeats. The mean value of Zd(φ) forms the PLL discriminator characteristic. The error signal is inputted to the loop filter, which in turn outputs a control signal and closes the PLL circuit, thereby controlling the frequency and phase of the PLL reference signal carrier NCO (also referred to as the carrier NCO). The carrier NCO produces reference signals for correlator-I and correlator-Q.

The DLL circuit operates in a similar way and includes similar components. The DLL discriminator forms tracking error signal Vd using the algorithm:

Vd=dI/I.  (E5)

The DLL reference signal code and strobe NCO (also referred to as the code and strobe NCO) generates two reference signals (reference oscillations) for the correlator-dI: one oscillation (code reference signal or code reference oscillation) in the form of a modulation code replica (or code sequence) to generate the in-phase correlation signal I, and the other (strobe reference signal or strobe reference oscillation) in the form of a strobe sequence which is sent to correlator-dI.

Discriminators operating according to expressions (E4) and (E5) can be modified if the in-phase correlation signal I is subjected to a preliminary smoothing process to reduce the interference component Iin. Smoothing the in-phase correlation signal I increases the operational stability of the tracking loops at low signal-to-noise ratios.
The PLL tracking system tries to make the tracking error φ equal to zero. As a result, one of the reference signals generated by the carrier NCO will be in-phase with the incoming signal, and its phase can be used to estimate the carrier phase of the incoming signal. This estimate is referred to as a phase measurement of the receiver.
The DLL tracking system tries to reduce the tracking error τ to zero. As a result, reference signals generated by the code and strobe NCO will correspond to the modulation code of the incoming signal. Delays of these reference signals can be used to estimate modulation code time delay of the incoming signal. This estimate is referred to as a code measurement of the receiver.
Phase and code measurements obtained from signals from a number of satellites are used to determine, for example, coordinates and speed of the rover. Due to external effects on tracking systems, the tracking errors fluctuate close to steady balance points in which error signals are equal to zero. At normal conditions in tracking modes, these errors are negligible. One source of tracking errors is the inertia of the tracking system. External disturbances (such as radio interference) are another source of tracking errors. Both types of errors depend on the inertia of the tracking system. In the tracking mode, the inertia can be estimated by the bandwidth of the equivalent filter, operation of which is equivalent to the corresponding tracking circuit (loop). The bandwidth of this filter characterizes the inertia of the tracking system and its filtering parameters. The bandwidth should be chosen based on a trade off between dynamic and interference errors.
Normal errors, such as caused by low-level noise, are accommodated during normal operation of the tracking systems; normal errors determine the accuracy of the measurements. Abnormal errors disrupt the normal operation of the tracking systems and can lead to corruption or loss of measurements. Abnormal errors can result from strong radio interference and strong multipath signals. Abnormal errors can also result from dynamic disturbances, such as shock and vibration, on the antenna and receiver. A particular type of abnormal phase error is a sudden change in the phase of the PLL carrier frequency (cycle slip). For a cycle slip of a short duration, the PLL transfers to a new stable point, and tracking continues; however, an abnormal error corresponding to an integer number of half cycles remains in the measurement of full phases. For a cycle slip of a long duration, or for a series of short-duration cycle slips, PLL tracking can be lost, and the tracking system needs to re-enter a search mode.
Embodiments of the invention disclosed herein reduce the impact of shock and vibration on measurements. Shock and vibration affect both the antenna and the quartz crystal oscillator in the receiver; however, mechanical forces have significantly different effects on the operation of the antenna and the quartz crystal oscillator.
Mechanical forces applied to the antenna cause displacements of the antenna and, consequently, phase shifts of the carrier signals received by the antenna. For a carrier signal transmitted from a specific satellite, the phase shift is proportional to the projection of the displacement of the antenna onto the radius vector between the rover and the specific satellite; that is, the phase shift depends on the specific satellite channel. As described in more detail below, the amplitude of the phase shift is proportional to the amplitude of the antenna oscillations (vibrations) and is not dependent on the frequency of the antenna oscillations. Derivatives (with respect to time) of the antenna oscillations depend on both the amplitude and frequency of the antenna oscillations. The derivatives of the antenna oscillations whose order is equal to or greater than the order of astaticism of the tracking system cause dynamic tracking errors and, in some cases, can result in cycle slips and interruptions in tracking for one or more satellite channels. But even if there is no tracking loss, dynamic errors are different for different satellite channels and can affect estimates of the rover coordinates.
As discussed above, the navigation receiver includes a quartz crystal oscillator (QCO) 708 (FIG. 2) that serves as both a reference frequency source and a reference time source (clock). Mechanical forces applied to the quartz crystal oscillator result in both a local oscillator frequency shift and a receiver clock offset, which are inter-related. Slight mechanical oscillations of the quartz crystal oscillator cause errors in the measurement of time delays; however, the resulting errors in the measurement of the coordinates are not significant, since receiver clock offsets are corrected in solving the navigation task.
The local oscillator frequency shifts are tracked by the PLL. Large amplitude mechanical oscillations, however, can disrupt the tracking operation. As rule, the PLL is used as a guiding tracking loop, and the DLL is used as a guided tracking loop. Consequently, the PLL normally has second order and generates frequency guiding indications (target designations) for the DLL. Thus the DLL is free from dynamic loads, and its order can be reduced up to the first order. The PLL has more stable tracking characteristics compared to the DLL; but if there is a loss of tracking in the PLL, there is a loss of tracking in the DLL as well.
As mentioned above, mechanical forces applied to the quartz crystal oscillator cause a shift of the local oscillator frequency ƒLO and, consequently, a shift of the intermediate frequencies of the satellite channels and their phases; the resulting phase errors are common to all the satellite channels. As discussed above, the signal of the intermediate frequency is generated as a difference of the input signal and the local oscillator signal. In the intermediate signal, therefore, there are both phase offsets caused by vibration of the antenna and phase offsets caused by vibration of the quartz crystal oscillator. To simplify the analysis, phase offsets of the intermediate signal can be treated as equivalent phase offsets of the input signal.
Consider the effects of a harmonic mechanical vibration when the antenna is mounted on the receiver case; the harmonic mechanical vibration is assumed to be identical for both the antenna and the quartz crystal oscillator. Consider the component of the mechanical vibration projected onto the radius vector between the rover and a specific satellite. The displacement of the receiver due to harmonic vibration can be represented by

Av=A0 sin(2πƒvt)=A0 sin(ωvt),  (E6)

where Av is the displacement (of the projected component), A0 is the amplitude of vibration (of the projected component), ƒv is the frequency of vibration, ωv=2πƒv is the angular frequency of vibration, and t is time. The speed of vibration is given by the first derivative with respect to time:

A′v=ωvA0 cos(ωvt);  (E7)

and the acceleration of vibration is given by the second derivative with respect to time:

A′′V=−ωv2A0 sin(ωvt).  (E8)

To simplify the calculation, ignore the frequency change in the input signal due to satellite movement. The equivalent perturbed frequency of the input signal, when both the antenna and the quartz crystal oscillator are subjected to the same vibration, can then be expressed as

ƒin=ƒ0+Δƒa+Δƒq,  (E9)

where:

    
    
        ƒin is the equivalent perturbed frequency of the input signal;
        ƒ0 is the nominal (unperturbed) frequency of the input signal (for the GPS C/A channel, for example, ƒ0=1575.42 MHz);
        Δƒa is the frequency change of the input signal due to the effects of vibration on the antenna; and
        Δƒq is the frequency change of the input signal due to the effects of vibration on the quartz crystal oscillator.
    
    


The frequency change of the input signal due to the effects of vibration on the antenna is proportional to the speed of vibration:




  
    
      
        
          
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  f
                  a
                
              
              =
              
              ⁢
              
                
                  A
                  v
                  ′
                
                λ
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  A
                  v
                  ′
                
                ⁡
                
                  (
                  
                    
                      f
                      0
                    
                    c
                  
                  )
                
              
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  
                    (
                    
                      1
                      c
                    
                    )
                  
                  ⁢
                  
                    f
                    0
                  
                  ⁢
                  
                    ω
                    v
                  
                  ⁢
                  
                    A
                    0
                  
                  ⁢
                  
                    cos
                    ⁡
                    
                      (
                      
                        
                          ω
                          v
                        
                        ⁢
                        t
                      
                      )
                    
                  
                
              
              ,
            
          
        
      
    
    
      
        (
        E10
        )
      
    
  




where λ is the (unperturbed) wavelength of the input signal, and C is the speed of light.

The frequency change of the input signal due to the effects of vibration on the quartz crystal oscillator is proportional to the acceleration of vibration:




  
    
      
        
          
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  f
                  q
                
              
              =
              
              ⁢
              
                
                  α
                  q
                
                ⁢
                
                  f
                  0
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  k
                  q
                
                ⁢
                
                  A
                  v
                  ′′
                
                ⁢
                
                  f
                  0
                
              
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  
                    -
                    
                      f
                      0
                    
                  
                  ⁢
                  
                    k
                    q
                  
                  ⁢
                  
                    A
                    0
                  
                  ⁢
                  
                    ω
                    v
                    2
                  
                  ⁢
                  
                    sin
                    ⁡
                    
                      (
                      
                        
                          ω
                          v
                        
                        ⁢
                        t
                      
                      )
                    
                  
                
              
              ,
            
          
        
      
    
    
      
        (
        E11
        )
      
    
  




where





  
    α
    q
  
  =
  
    
      δ
      ⁢
      
        
      
      ⁢
      
        f
        q
      
    
    
      f
      
        0
        ⁢
        q
      
    
  




is the shift of the frequency of the quartz crystal oscillator relative to the nominal value of the frequency of the quartz crystal oscillator;

    
    
        ƒ0q is the nominal value of the frequency of the quartz crystal oscillator;
        δƒq is the shift of the frequency of the quartz crystal oscillator from ƒ0q; and
        kq is a coefficient of proportionality that characterizes the sensitivity of the quartz crystal oscillator to acceleration effects.

The coefficient kq depends on various factors, such as the direction of acceleration relative to the crystal axes of the quartz crystal, the specific type of quartz crystal, and the means by which the quartz crystal is mounted in the receiver case. For quartz crystal oscillators typically used in receivers, kq≈−0.102×10−9.

    
    


Let dƒa≡ΔƒaMAX represent the maximum value of Δƒa, and dƒq≡ΔƒqMAX represent the maximum value of Δƒq; then




  
    
      
        
          
            
              
                
                  df
                  a
                
                ≡
                
                ⁢
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    f
                    aMAX
                  
                
              
            
          
          
            
              
                =
                
                ⁢
                
                  
                    (
                    
                      1
                      c
                    
                    )
                  
                  ⁢
                  
                    f
                    0
                  
                  ⁢
                  
                    A
                    0
                  
                  ⁢
                  
                    ω
                    v
                  
                
              
            
          
          
            
              
                
                  =
                  
                  ⁢
                  
                    
                      (
                      
                        
                          2
                          ⁢
                          π
                        
                        c
                      
                      )
                    
                    ⁢
                    
                      f
                      0
                    
                    ⁢
                    
                      A
                      0
                    
                    ⁢
                    
                      f
                      v
                    
                  
                
                ,
              
            
          
        
        ⁢
        
          

        
        ⁢
        and
      
    
    
      
        (
        E12
        )
      
    
  
  
    
      
        
          
            
              
                df
                q
              
              ≡
              
              ⁢
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  f
                  qMAX
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  f
                  0
                
                ⁢
                
                  k
                  q
                
                ⁢
                
                  A
                  0
                
                ⁢
                
                  ω
                  v
                  2
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                4
                ⁢
                
                  π
                  2
                
                ⁢
                
                  f
                  0
                
                ⁢
                
                  k
                  q
                
                ⁢
                
                  A
                  0
                
                ⁢
                
                  
                    f
                    v
                    2
                  
                  .
                
              
            
          
        
      
    
    
      
        (
        E13
        )
      
    
  




Refer to FIG. 9. Plot 902 shows a representative plot of dƒa (Hz) as a function of ƒv (Hz) according to (E12); similarly, plot 904 shows a representative plot of dƒq (Hz) as a function of ƒv according to (E13). For these plots, an amplitude of A0=0.2 m was used. From (E12), dƒa varies linearly as a function of ƒv; from (E13), dƒq varies quadratically as a function of ƒv. For ƒv< ̃5 Hz, dƒa>dƒq; for ƒv> ̃5 Hz, dƒq grows increasingly larger than dƒa.

If the antenna is moving at a constant velocity, the carrier frequency is offset proportionally to the component the velocity vector projected onto the radius vector between the rover and a specific satellite; that is, the carrier frequency offset depends on the specific satellite channel. Shock and vibration cause acceleration of the antenna. The effect of shock and vibration on GNSS measurements depends on the projection of the acceleration vector onto the radius vector between the antenna and a specific satellite; that is, the effect of shock and vibration also depends on the specific satellite channel.
The frequency of the quartz crystal oscillator is offset proportionally to the acceleration of the quartz crystal oscillator; if the quartz crystal oscillator is moving at a constant velocity (zero acceleration), then the frequency offset of the quartz crystal oscillator is zero. A quartz crystal is anisotropic, and the frequency offset depends on the acceleration vector relative to the crystal axes. Typically, the frequency offset is the most sensitive along one crystal axis (the main crystal axis), and the frequency offset is proportional to the component of the acceleration vector projected onto the main crystal axis. The frequency offset of the quartz crystal oscillator is therefore the same for all satellite channels.
FIG. 3A and FIG. 3B show a schematic functional block diagram of the phase-lock loops, according to an embodiment of the invention, in the digital channel processing unit 206 (FIG. 2). The input digital signal 205 (FIG. 2) is processed by a set of individual digital channel processors; as discussed above, each specific satellite channel corresponds to navigation signals transmitted from a specific navigation satellite on a specific frequency band. In the example shown in FIG. 1, there are five satellite channels, corresponding to navigation satellite 120A—navigation satellite 120E, respectively. In general, on a single frequency band, there are N satellite channels, where N is an integer referring to the number of satellites in view from the rover. A specific satellite channel is indexed as satellite channel j, where j is an integer j=1, . . . , N.
The structure of three representative digital channel processors is shown in FIG. 3A. Operations for all digital channel processors are similar; detailed operation of the j-th digital channel processor is described in reference to FIG. 3B.
First refer to FIG. 3A. The first digital channel processor includes the first channel correlator 302-1, the first reset accumulator 304-1 (ACC1), the first channel discriminator 306-1, the first loop filter 308-1 (LF1), and the first numerically-controlled oscillator 310-1 (NCO1). Similarly, the N-th digital channel processor includes the N-th channel correlator 302-N, the N-th reset accumulator 304-N (ACCN), the N-th channel discriminator 306-N, the N-th loop filter 308-N (LFN), and the N-th numerically-controlled oscillator 310-N (NCON).
FIG. 3B shows the structure of the j-th digital channel processor. The j-th digital channel processor includes the j-th channel correlator 302-j, the j-th channel reset accumulator 304-j (ACCj), the j-th channel discriminator 306-j, the j-th channel loop filter 308-j (LFj), and the j-th channel numerically-controlled oscillator 310-j (NCOj). The digital signal 205 is inputted into the correlator 302-j, which processes the digital signals corresponding to the j-th satellite. The correlator 302-j outputs the in-phase correlation signal 301 Ij [according to (E1)] and the quadrature-phase correlation signal 303 Qj [according to (E2)]. The correlator 302-j outputs the signal Ij and the signal Qj at a high clock frequency (for example, 1 kHz). The values Ij and Qj are inputted into the reset accumulator 304-j and the common coordinate-quartz loop (CCQL) discriminator 320. More details of the CCQL discriminator are described below.
The reset accumulator 304-j performs an accumulate-and-reset operation, also referred to as an accumulate-and-dump operation. A reset accumulator initializes its value to zero, receives input values of a variable, sums the input values over an accumulation interval ΔT, and outputs the sum. The reset accumulator then repeats the cycle of operation for the next set of input values. The reset accumulator outputs the signal 305 I*j; and the signal 307 Q*j(defined below) at a much lower frequency than the clock frequency for Ij and Qj; for example, a frequency of 200 Hz. The reset accumulator 304-j receives the in-phase correlation signal and outputs the accumulated in-phase correlation signal I*j=ΣIj (information symbols are removed prior to accumulation). Similarly, the reset accumulator 304-j receives the quadrature-phase correlation signal Qj and outputs the accumulated quadrature-phase correlation signal Q*j=ΣQj (information symbols are removed prior to accumulation).
The values of I*j and Q*j are inputted into the discriminator 306-j, which produces the PLL tracking error signal 311,




  
    Zdr
    j
  
  =
  
    
      arctan
      ⁡
      
        (
        
          
            Q
            j
            *
          
          
            I
            j
            *
          
        
        )
      
    
    .
  




The PLL tracking error signal Zdrj is inputted into the loop filter 308-j. The loop filter 308-j determines the PLL order, mostly with one or two accumulators assigning the second or third order, respectively. In some embodiments, the common phase-control signal 335 Zφcj, described in more detail below, is also inputted from the CCQL discriminator 320 into the loop filter 308-j. 

In the embodiment shown in FIG. 3B, the loop filter 308-j outputs the individual channel frequency-control signal 321 Zθrj and the individual channel phase-control signal 323 Zφrj. The CCQL discriminator 320 outputs the common frequency-control signal 331 Zƒcj and the common phase-control signal 333 Zφcj. These four control signals are used to control operational parameters of the NCO 310-j (NCOj), as described in more detail below in reference to FIG. 8. Other embodiments, as described below, use fewer than four control signals. Details of the control signals are also described below.
FIG. 4A shows a schematic functional block diagram for a first embodiment of the CCQL discriminator 320 [this embodiment is referenced as the CCQL discriminator 320(A)]. A vector of initial channel phase tracking errors (mismatches) is first calculated from the in-phase correlation signals and the quadrature-phase correlation signals; this vector is denoted

Zdc=(Zdc1, . . . ,Zdcj, . . . ZdcN).  (E14)

The vector components are calculated schematically as shown in FIG. 4A. For the first channel, the in-phase correlation signal I1 and the quadrature-phase correlation signal Q1 are inputted into the operator 402-1 to calculate (CALC) Zdc1; the output is the value Zdc1. For the j-th channel, the in-phase correlation signal Ij and the quadrature-phase correlation signal Qj are inputted into the operator 402-j to calculate (CALC) Zdcj; the output is the value Zdcj. For the N-th channel, the in-phase correlation signal IN and the quadrature-phase correlation signal QN are inputted into the operator 402-N to calculate (CALC) ZdcN; the output is the value ZdcN.

In an embodiment, the vector component Zdcj is calculated as




  
    
      
        
          Zdc
          j
        
        =
        
          
            arctan
            ⁡
            
              (
              
                
                  Q
                  j
                
                
                  I
                  j
                
              
              )
            
          
          .
        
      
    
    
      
        (
        E15
        )
      
    
  




In another embodiment, the vector component Zdcj is calculated as

Zdcj=Ij*Qj.  (E16)

The vector Zdc 401 and the matrix H 209A (see below) are inputted into the operator 404, in which a four-dimensional vector Zs 411 is calculated using the least-mean-squares (LSM) method. The vector Zs=(Zx,Zy,Zz,Zq) is a mean-weighted tracking error vector, where the three components (Zx,Zy,Zz) correspond to the projections of the tracking error vector onto the (x,y,z) Cartesian axes; and Zq is the time component of the tracking error vector, Zq=cΔt, where is the speed of light, and Δt is the tracking error of the PLL for the time coordinate t.
The matrix H 209A is calculated in the user data processing unit 208 and fed back to the digital channel processing unit 206 as part of the output 209 (FIG. 2). The matrix H is the matrix of directional cosines supplemented by a unit column:




  
    
      
        H
        =
        
          
            
            
              
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      α
                      1
                    
                  
                
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      β
                      1
                    
                  
                
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      γ
                      1
                    
                  
                
                
                  1
                
              
              
                
                  ...
                
                
                  ...
                
                
                  ...
                
                
                  ...
                
              
              
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      α
                      j
                    
                  
                
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      β
                      j
                    
                  
                
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      γ
                      j
                    
                  
                
                
                  1
                
              
              
                
                  ...
                
                
                  ...
                
                
                  ...
                
                
                  ...
                
              
              
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      α
                      N
                    
                  
                
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      β
                      N
                    
                  
                
                
                  
                    cos
                    ⁢
                    
                      
                    
                    ⁢
                    
                      γ
                      N
                    
                  
                
                
                  1
                
              
            
            
          
          .
        
      
    
    
      
        (
        E17
        )
      
    
  




The angles (αj,βj,γj) refer to the angles between the radius vector from the rover to the j-th satellite and the three corresponding axes of the (x,y,z) Cartesian coordinate system.
The operator 404 first calculates a matrix W, an N×N matrix of weight estimates for measurements. The matrix W depends on how Zdcj above is calculated. For




  
    
      
        
          
            Zdc
            j
          
          =
          
            arctan
            ⁡
            
              (
              
                
                  Q
                  j
                
                
                  I
                  j
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        E15
        )
      
    
  




the matrix W is given by:





  
    
      
        
          W
          ij
        
        =
        
          {
          
            
              
                
                  
                    I
                    ~
                  
                  j
                  2
                
              
              
                
                  
                    i
                    =
                    j
                  
                  ;
                
              
            
            
              
                0
              
              
                
                  
                    i
                    ≠
                    j
                  
                  ;
                
              
            
          
        
      
    
    
      
        (
        E18
        )
      
    
  




where Ĩj is the smoothed in-phase component of the j-th channel. For Zdcj=Ij*Qj (E16), the matrix W is given by:





  
    
      
        
          W
          ij
        
        =
        
          {
          
            
              
                
                  I
                  j
                  2
                
              
              
                
                  i
                  =
                  
                    j
                    .
                  
                
              
            
            
              
                0
              
              
                
                  i
                  ≠
                  
                    j
                    .
                  
                
              
            
          
        
      
    
    
      
        (
        E19
        )
      
    
  




From the matrix W and the matrix H, the matrix G is calculated, where

G=(HT*W*H)−1*HT*W.  (E20)

For an arbitrary matrix M, MT refers to the transpose of matrix M, and M−1 refers to the inverse of matrix M. The vector Zs is then calculated as

Zs=G*Zdc.  (E21)

The vector Zs=(Zx,Zy,Zy,Zq) is then transformed by an operator F=(Fx,Fy,Fy,Fq). The transformation performed by the operator Fs depends on the CCQL order and on the CCQL control mode. There are three CCQL control modes: frequency control (FC), phase control (PC), and frequency and phase control (FPC). These control modes are described in more detail below. For frequency control, the operator Fs transforms the vector Zs into the vector Zs1, where Zs1=(Zx1,Zy1,Zz1,Zq1). For phase control, the operator Fs transforms the vector Zs into the vector Zs2, where Zs2=(Zx2,Zy2,Zz2,Zq2). For frequency and phase control, the operator Fs transforms the vector Zs into the pair of vectors Zs1 and Zs2.
The operator Fs acts as a loop filter whose output is a smoothed function of the input. The vector Zs1=(Zx1,Zy1,Zz1,Zq1) is the smoothed tracking error vector. The vector Zs2=(Zx2,Zy2,Zz2,Zq2) is the smoothed tracking error rate vector, where the tracking error rate is proportional to the time derivative of the tracking error.
The individual loops and the common quartz loop solve different tasks; therefore, their orders can be different. Individual loops track slow phase changes caused, for example, by satellite acceleration and smooth accelerations of the rover (for example, during smooth start of the vehicle); a third-order PLL can be used in the individual loops. If these accelerations are negligible and can be neglected, then a second-order PLL can be used in individual loops. If there are intense dynamic disturbances, then depending on their nature, a common quartz loop of a specific order can be used.
FIG. 4A schematically shows the transformation of Zs 411 by the operator Fs 406, which generates both Zs1 413 and Zs2 415. The vector component Zx is inputted into the operator Fx 406-x; the outputs are Zx1 and Zx2. The vector component Zy is inputted into the operator Fy 406-y; the outputs are Zy1 and Zy2. The vector component Zz is inputted into the operator Fz 406-z; the outputs are Zz1 and Zz2. The vector component Zq is inputted into the operator Fq 406-q; the outputs are Zq1 and Zq2. The output vector components are inputted into the group operator 420, which yields the vector Zs1=(Zx1,Zy1,Zz1,Zq1) 413 and the vector Zs2=(Zx2,Zy2,Zz2,Zq2) 415. In general, the operators Fx,Fy,Fz, and Fq are independent. In some embodiments, all of the operators are the same. Embodiments of loop filters for implementing the operator Fs are described below.
The vector Zs1 413, the vector Zs2 415, and the matrix H 209A are inputted into the operator 408, which calculates the common phase-control vector Zφc and the common frequency-control vector Zƒc. The common phase-control vector and the common frequency-control vector are calculated by projecting the vectors Zs1 and Zs2, respectively, onto the radius vectors between the rover and the satellites:

Zφc=H*Zs1;  (E22)

and

Zƒc=H*Zs2,  (E23)

where

Zφc=(Zφc1, . . . ,Zφcj, . . . ,ZφcN),  (E24)

and

Zƒc=(Zƒc1, . . . ,Zƒcj, . . . ,ZƒcN).  (E25)

Representative components are referenced for the j-th channel: the common frequency-control vector Zƒcj 331(A) and the common phase-control vector Zφcj 333(A) are inputted into the NCO 310-j [here the common frequency-control vector 331(A) is a first embodiment of the common frequency-control vector 331 in FIG. 3B, and the common phase-control vector 333(A) is a first embodiment of the common phase-control vector 333 in FIG. 3B]. In this embodiment, there is no common phase-control vector 335 inputted into the loop filter 308-j. 

As discussed above, in the embodiment shown in FIG. 4A, the operator 406 generates both the vector Zs1 and the vector Zs2, and the operator 408 calculates both the common phase-control vector Zφc and the common frequency-control vector Zƒc. In another embodiment, the operator 406 generates only the vector Zs1, and the operator 408 calculates only the common phase-control vector Zφc. In a further embodiment, the operator 406 generates only the vector Zs2, and the operator 408 calculates only the common frequency-control vector Zƒc. FIG. 6, described below, shows an embodiment of the operator 406 that, through the operation of two switches, can generate all three combinations of outputs: only the vector Zs1, only the vector Zs2, or both the vector Zs1 and the vector Zs2.
FIG. 4B shows a schematic functional block diagram for a second embodiment of the CCQL discriminator 320 [this embodiment is referenced as the CCQL discriminator 320(B)]. Compared to the CCQL discriminator 320(A) shown in FIG. 4A, the operator Fs 416 is simpler than the operator Fs 406, and the operator 418 is simpler than the operator 408. The CCQL discriminator 320(B) outputs only common phase-control signals (no common frequency-control signals are outputted). The operator Fs 416 is implemented as a loop filter. The vector Zs 411 is inputted into the operator Fs 416, which multiplies Zs by a coefficient K1,c to generate the output vector Zs1 431. The vector Zs1 431 and the matrix H 209A are inputted into the operator 418, which calculates the common phase-control vector Zφc according to (E22) above.
Each component of Zφc is inputted in parallel to two on/off switches. In FIG. 4B, three representative components are shown: Zφc1 is inputted to switch 410-1 and switch 412-1; Zφcj is inputted to switch 410-j and switch 412-j; and ZφcN is inputted to switch 410-N and switch 412-N. Let the switch group PC refer to the group of switches (410-1, . . . , 410-j, . . . , 410-N); and let the switch group FC refer to the group of switches (412-1, . . . , 412-j, . . . , 412-N). As discussed in more detail below, there are three control modes [phase control (PC), frequency control (FC), and frequency and phase control (FPC) mode] for the numerically-controlled oscillator 310-j (FIG. 3B). Depending on the desired control mode for the numerically-controlled oscillator, the switch group PC and the switch group FC are opened or closed. When the PC mode is implemented, the switch group PC is closed, and the switch group FC is opened. When the FC mode is implemented, the switch group FC is closed, and the switch group PC is opened. When the FPC mode is implemented, the switch group PC and the switch group FC are both closed. Representative components are referenced for the j-th channel: when the switch 410-j is closed, the common phase-control vector Zφcj 333(B) is inputted into the NCO 310-j [here the common phase-control vector 333(B) is a second embodiment of the common phase-control vector 333 in FIG. 3B]. In this embodiment, when the switch 412-j is closed, the common phase-control vector 335(B) is inputted into the loop filter 308-j. 
FIG. 5A shows a schematic functional block diagram for a first embodiment of the loop filter 308-j for the individual j-th satellite channel (FIG. 3B). This loop filter is referenced as the loop filter 308-j (A). The solid lines represent the functional blocks to implement a second-order PLL. The dotted lines represent additional functional blocks (in addition to those for the second-order loop filter) to implement a third-order PLL. A second-order PLL can be used when the rover moves along an approximately straight-line trajectory. A third-order PLL can be used when the rover moves along an approximately parabolic trajectory (assuming the time derivative of the rover acceleration is sufficiently small). The second-order loop filter is described first, followed by a description of the third-order loop filter. The loop filter 308-j (A) operates according to a clock with an update period of Tcr; the update frequency is ƒcr=1/Tcr. Each cycle of operation of the loop filter 308-j (A) occurs at a set of discrete time instants tcr; the period between consecutive time instants tcr, is the update period Tcr tcr,i+1−tcr,i, where i is an arbitrary integer index.
The tracking error signal 311 Zdrj outputted from the j-th channel discriminator 306-j (FIG. 3B) is inputted into the operator 502, which multiplies Zdrj by the coefficient K1,r. The output signal 503 (K1,rZdrj) is inputted into the switch 504, which has two output positions: switch output 504A for frequency control (FC) of the individual phase-lock loop and switch output 504B for frequency and phase control (FPC) of the individual phase-lock loop. Older generation navigation receivers typically use only frequency control; more modern navigation receivers typically use frequency and phase control, which, in some instances, yield smaller errors than frequency control alone. A single module can be implemented for both FC and FPC modes; the switch 504 allows selection of the desired mode. In other embodiments, a module is implemented for only FC mode or only FPC mode.
For a second-order loop filter with FC operation, the switch 504 is switched to the switch output 504A, and the output signal 503 is inputted into the sum operator 506. The tracking error signal 311 Zdrj is also inputted into the operator 520, which multiplies Zdrj by the coefficient K2,r. The output signal 521 (K2,rZdrj) is inputted into the accumulator 522, which accumulates values of the output signal 521. Note: In FIG. 5A, FIG. 5B, and FIG. 6B, the accumulators represented by the symbol Σ are not reset accumulators. The output signal 523 (accumulated values of output signal 521) is inputted into the sum operator 506. The output signal 505 (sum of output signal 503 and output signal 523) is inputted into the operator 508, which multiplies the output signal 505 by the value 1/Tcr.
In the embodiment shown, the output signal 507 is inputted into the sum operator 510. A predicted Doppler shift 209B ƒd SV j for the j-th satellite is also inputted into the sum operator 510. Refer to FIG. 2. The predicted Doppler shift 209B is calculated in the user data processing unit 208 (which calculates the radial velocity of a satellite using ephemeris data and solutions of the navigation task) and is included in the output 209 fed back from the user data processing unit 208 to the digital channel processing unit 206. The predicted Doppler shift 209B unloads the accumulator 522 (reduces the values stored in the accumulator and increases the accuracy of representation of these values). The output signal 321(A) Zƒrj is outputted to the NCO 310-j [here the output signal 321(A) is a first embodiment of the output signal 321 in FIG. 3B]. Calculation of the predicted Doppler shift 209B, however, increases the computational load and is not used in some embodiments; then the output signal 321(A) Zƒrj is the same as the output signal 507.
For a second-order loop filter with FPC operation, the switch 504 is switched to the switch output 504B, and the output signal 503 is outputted to the NCO 310-j as the output signal 323(A) Zφrj [here the output signal 323(A) is a first embodiment of the output signal 323 in FIG. 3B]. The output signal 523 is inputted into the sum operator 506. Since the switch 504 is switched to the switch output 504B, the output signal 505 is the same as the output signal 523. The output signal 505 is then processed as described above in reference to FC operation.
For a third-order loop filter with FC operation, the switch 504 is switched to the switch output 504A, and the output signal 503 is inputted into the sum operator 506. The tracking error signal 311 Zdrj is also inputted into the operator 520, which multiplies Zdrj by the coefficient K2,r. The output signal 521 (K2,rZdrj) is inputted into the accumulator 522. The tracking error signal 311 Zdrj is further inputted into the operator 530, which multiplies Zdrj by the coefficient K3,r. The output signal 531 (K3,rZdrj) is inputted into the accumulator 532, which accumulates the output signal 531. The output signal 533 (accumulated values of output signal 531) is inputted into the accumulator 522. The output signal 523 (accumulated values of output signal 521 and output signal 533) is inputted into the sum operator 506. The output signal 505 is then processed as described above in reference to the second-order loop filter with FC operation.
For a third-order loop filter with FPC operation, the switch 504 is switched to the switch output 504B, and the output signal 503 is outputted to the NCO 310-j as the output signal 323(A) Zφrj. The output signal 523 is inputted into the sum operator 506. Since the switch 504 is switched to the switch output 504B, the output signal 505 is the same as the output signal 523. The output signal 505 is then processed as described above in reference to the second-order loop filter with FPC operation.
FIG. 5B shows a schematic functional block diagram for a second embodiment of the loop filter 308-j for the individual j-th satellite channel (FIG. 3B). This loop filter is referenced as the loop filter 308-j (B). The solid lines represent the functional blocks to implement a second-order PLL. The dotted lines represent additional functional blocks (in addition to those for the second-order loop filter) to implement a third-order PLL. The second-order loop filter is described first; a description of the third-order loop filter then follows. The loop filter 308-j (B) operates according to a clock with an update period of Tcr; the update frequency is ƒcr=1/Tcr. Each cycle of operation of the loop filter 308-j (B) occurs at a set of discrete time instants tcr; the period between consecutive time instants tcr is the update period Tcr=(tcr,i+1−tcr,i, where i is an arbitrary integer index.
The tracking error signal 311 Zdrj (FIG. 3B) is inputted into the operator 502, which multiplies Zdrj by the coefficient K1,r. The output signal 503 (K1,rZdrj) is inputted into the switch 504, which has two output positions: switch output 504A for frequency control (FC) of the individual phase-lock loop and switch output 504B for frequency and phase control (FPC) of the individual phase-lock loop.
For a second-order loop filter with FC control, the switch 504 is switched to the switch output 504A, and the output signal 503 is inputted into the sum operator 506. The tracking error signal 311 Zdrj is also inputted into the operator 520, which multiplies Zdrj by the coefficient K2,r. The output signal 521 (K2,rZdrj) is inputted into the accumulator 522. The common phase-control signal 335(B) Zφcj from the CCQL discriminator 320 (FIG. 3B and FIG. 4B) is inputted into the operator 542, which multiplies Zφcj by the coefficient K2,c/K1,c. The output signal 543 [(K2,c/K1,c)Zφcj] is inputted into the accumulator 522.
The output signal 523 (accumulated values of output signal 521 and output signal 543) is inputted into the sum operator 506. The output signal 505 (sum of output signal 503 and output signal 523) is then processed as described above in reference to the loop filter 308-j (A). The output signal 321(B) Zƒrj is outputted to the NCO 310-j [here the output signal 321(B) is a second embodiment of the output signal 321 in FIG. 3B].
For a second-order loop filter with FPC operation, the switch 504 is switched to the switch output 504B, and the output signal 503 is outputted to the NCO 310-j as the output signal 323(B) Zφrj [here the output signal 323(B) is a second embodiment of the output signal 323 in FIG. 3B]. The output signal 523 is inputted into the sum operator 506. Since the switch 504 is switched to the switch output 504B, the output signal 505 is the same as the output signal 523. The output signal 505 is then processed as described above in reference to the loop filter 308-j (A). The output signal 321(B) Zƒrj is outputted to the NCO 310-j. 
For a third-order loop filter with FC operation, the switch 504 is switched to the switch output 504A, and the output signal 503 is inputted into the sum operator 506. The tracking error signal 311 Zdrj is also inputted into the operator 520, which multiplies Zdrj by the coefficient K2,r. The output signal 521 (K2,rZdrj) is inputted into the accumulator 522. The common phase-control signal 335(B) Zφcj from the CCQL discriminator 320 (FIG. 3B and FIG. 4B) is inputted into the operator 542, which multiplies Zφcj by the coefficient K2,c/K1,c. The output signal 543 [(K2,c/K1,c)Zφcj] is inputted into the accumulator 522.
The tracking error signal 311 Zdrj is also inputted into the operator 530, which multiplies Zdrj by the coefficient K3,r. The output signal 531 (K3,rZdrj) is inputted into the accumulator 532. The common phase-control signal 335(B) Zφcj is also inputted into the operator 552, which multiplies Zφcj by the coefficient K3,c/K1,c. The output signal 553 [(K3,c/K1,c)Zφcj] is inputted into the accumulator 532. The accumulator 532 accumulates values of output signal 531 and output signal 553.
The output signal 533 (accumulated values of output signal 531 and output signal 553) is inputted into accumulator 522. The accumulator 522 accumulates values of output signal 521, output signal 543, and output signal 533. The output signal 523 (accumulated values of output signal 521, output signal 543, and output signal 533) is inputted into the sum operator 506. The output signal 505 (sum of output signal 503 and output signal 523) is then processed as described above in reference to the loop filter 308-j (A). The output signal 321(B) is outputted to the NCO 310-j (FIG. 3B).
For a third-order PLL filter with FPC operation, the switch 504 is switched to the switch output 504B, and the output signal 503 is outputted to the NCO 310-j as the output signal 323(B) Zφrj. The output signal 523 is inputted into the sum operator 506. Since the switch 504 is switched to the switch output 504B, the output signal 505 is the same as the output signal 523. The output signal 505 is then processed as described above in reference to the loop filter 308-j (A). The output signal 321(B) is outputted to the NCO 310-j (FIG. 3B).
FIG. 6 shows a schematic functional block diagram for an embodiment of the Fs operator 406 in the CCQL discriminator 320 (FIG. 4A). The loop filter 406 operates according to a clock with an update period of Tcc; the update frequency is ƒcc=1/Tcc. Each cycle of operation of the loop filter 406 occurs at a set of discrete time instants tcc; the period between consecutive time instants tcc is the update period Tcc=tcc,i+1−tcc,i, where i is an arbitrary integer index.
The input to the loop filter 406 is the vector Zs=(Zx,Zy,Zz,Zq) 411, which is received from the operator 404 (FIG. 4A). The outputs of the loop filter 406 are determined by the position of the switch 640 and by the position of the switch 604 (FIG. 6). First consider the configuration in which the switch 640 is closed. When the switch 604 is switched to the switch output 604A [frequency control (FC)], Fs transforms the vector Zs=(Zx,Zy, Zy,Zq) into the single vector Zs2=(Zx2,Zy2,Zz2,Zq2). When the switch 604 is switched to the switch output 604B (FPC), Fs transforms the vector Zs=(Zx,Zy,Zy,Zq) into two vectors: the vector Zs1=(Zx1,Zy1,Zz1,Zq1) and the vector Zs2=(Zx2, Zy2, Zz2, Zq2). Now consider the configuration in which the switch 640 is open. In this instance, the switch 604 is switched only to the switch output 604(B) [phase control (PC)]; Fs transforms the vector Zs=(Zx,Zy,Zy,Zq) into the single vector Zs1=(Zx1,Zy1,Zz1,Zq1). Details of the signal processing are described below.
First consider the configuration in which the switch 640 is closed. The solid lines represent the functional blocks to implement a second-order CCQL. The dotted lines represent additional functional blocks (in addition to those for the second-order CCQL) to implement a third-order CCQL. The second-order loop filter is described first; a description of the third-order loop filter then follows.
The input signal 411 Zs is inputted into the operator 602, which multiplies Zs by the coefficient K1,sc. The output signal 603 (K1,scZs) is inputted into the switch 604, which has two output positions: switch output 604A for frequency control (FC) of the NCOs and switch output 604B for frequency and phase control (FPC) of the NCOs. For a second-order loop filter with FC operation, the switch 604 is switched to the switch output 604A, and the output signal 603 is inputted into the sum operator 606. The input signal 411 Zs is also inputted into the operator 620, which multiplies Zs by the coefficient K2,sc. The output signal 621 (K2,scZs) is inputted into the accumulator 622, which accumulates the output signal 621. The output signal 623 (accumulated values of the output signal 621) is inputted into the sum operator 606. The output signal 605 (sum of output signal 603 and output signal 623) is inputted into the operator 608, which multiplies the output signal 605 by the value 1/Tcc. The output signal 415 Zs2 is outputted to the operator 408 (FIG. 4A).
For a second-order loop filter with FPC operation, the switch 604 is switched to the switch output 604B, and the output signal 603 is outputted as the output signal 413 Zs1 to the operator 408 (FIG. 4A). The output signal 623 is inputted into the sum operator 606. Since the switch 604 is switched to the switch output 604B, the output signal 605 is the same as the output signal 623. The output signal 605 is then processed as described above in reference to FC operation. The output signal 415 Zs2 is outputted to the operator 408 (FIG. 4A).
For a third-order loop filter with FC operation, the switch 604 is switched to the switch output 604A, and the output signal 603 is inputted into the sum operator 606. The input signal 411 Zs is also inputted into the operator 620, which multiplies Zs by the coefficient K2,sc. The output signal 621 (K2,scZs) is inputted into the accumulator 622. The input signal 411 Zs is further inputted into the operator 630, which multiplies Zs by the coefficient K3,sc. The output signal 631 (K3,scZs) is inputted into the accumulator 632. The output signal 633 (accumulated values of output signal 631) is inputted into the accumulator 622.
The output signal 623 (accumulated values of output signal 621 and output signal 633) is inputted into the sum operator 606. The output signal 605 (sum of output signal 603 and output signal 623) is inputted into the operator 608, which multiplies the output signal 605 by the value 1/Tcc. The output signal 415 Zs2 is outputted to the operator 408 (FIG. 4A).
For a third-order loop filter with FPC operation, the switch 604 is switched to the switch output 604B, and the output signal 413 Zs1 is outputted to the operator 408 (FIG. 4A). The output signal 623 is inputted into the sum operator 606. Since the switch 604 is switched to the switch output 604B, the output signal 605 is the same as the output signal 623. The output signal 605 is then processed as described above in reference to the third-order loop filter with FC operation. The output signal 415 Zs2 is outputted to the operator 408 (FIG. 4A).
Now consider the configuration in which the switch 640 is open and the switch 604 is switched to the switch output 604B. The phase control (PC) mode of a first-order CCQL is implemented. The input signal 411 Zs is inputted into the operator 602, which multiplies Zs by the coefficient K1,sc, and the output signal 603 (K1,scZs) is outputted as the output signal 413 Zs1 to the operator 408 (FIG. 4A).
FIG. 8 shows a schematic functional block diagram of a representative numerically-controlled oscillator (NCO). Shown is the NCO 310-j (NCOj) for the j-th satellite channel (FIG. 3B). The NCO 310-j operates according to a clock with an update period of Tsampl; the update (sampling) frequency is ƒsampl=1/Tsampl. The clock generator 870 receives the QCO reference signal 211 from the quartz crystal oscillator 708 (FIG. 7) and generates the clock signal 341 with a clock frequency of ƒsampl. Each cycle of operation of the NCO 310-j occurs at a set of discrete time instants ts; the period between consecutive time instants ts is the update period Tsampl=ts,i+1−ts,i, where is an arbitrary integer index. The NCO 310-j generates the output reference signal 351 (RefSigj), which is inputted to the channel correlator 302-j (FIG. 3B). The RefSigj has the components cos(ΦNCOj) and sin(ΦNCOj).
The RefSigj has the following mathematical representation. Assume A(t) represents a harmonic oscillation in which the amplitude has been normalized to one:

A(t)=sin(2πƒt),  (E26)

    
    
        where
        t is the time,
        ƒ is the frequency, and
        Φ=2πƒt is the full phase of the harmonic oscillation.

The harmonic oscillation is a function with period of 2π and can be written as

A(t)=sin(2πƒt)=sin(2πƒt−2π·n),  (E27)

where n is an arbitrary integer number. Let

    
    






  
    n
    =
    
      E
      ⁡
      
        [
        
          
            2
            ⁢
            π
            ⁢
            
              
            
            ⁢
            ft
          
          
            2
            ⁢
            π
          
        
        ]
      
    
  
  ,




where E[x] is the integer part of x. Then





  
    
      
        
          
            
              
                Φ
                NCO
              
              =
              
              ⁢
              
                (
                
                  
                    2
                    ⁢
                    
                      π
                      ⁢
                      ft
                    
                  
                  -
                  
                    2
                    ⁢
                    
                      π
                      ·
                      
                        E
                        ⁡
                        
                          [
                          
                            
                              2
                              ⁢
                              
                                π
                                ⁢
                                ft
                              
                            
                            
                              2
                              ⁢
                              π
                            
                          
                          ]
                        
                      
                    
                  
                
                )
              
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  FP
                  ⁡
                  
                    [
                    
                      
                        2
                        ⁢
                        
                          π
                          ⁢
                          ft
                        
                      
                      
                        2
                        ⁢
                        π
                      
                    
                    ]
                  
                
              
              ,
            
          
        
      
    
    
      
        (
        E28
        )
      
    
  




where FP[x] is the fractional part of x. FP[Φ] is referred to as the fractional part of the full phase Φ. For the j-th satellite channel, ΦNCOj is the fractional part of the full phase Φj. The RefSigj then has the components cos(ΦNCOj) and sin(ΦNCOj).

The two principal components of the NCO 310-j are the phase accumulator 802 and the phase-to-amplitude converter 830. The phase accumulator 802 includes the sum operators 860 and the register 812. Various input values, described below, are summed in the sum operators 860. The output 827 of the sum operators 860 is recorded in the register 812 at the time instants ts. Assume that time instants ts (and tcr and tcc below) are referenced to a common system clock. The output 829 of the register 812 at the time instant ts is the value ΦNCOj, which is inputted into the phase-to-amplitude converter 830, which then generates the value of the in-phase reference signal 833 cos(ΦNCOj) and the value of the quadrature-phase reference signal 831 sin(ΦNCOj), which are components of the RefSigj 351. To decrease computational load, the phase-to-amplitude converter 830 can include a lookup table, stored in read-only memory (ROM), of values of cos(ΦNCOj) and sin(ΦNCOj) as functions of ΦNCOj.
In the embodiment shown in FIG. 8, four input values are accumulated in the phase accumulator 802. In other embodiments, fewer than four input values are accumulated; for example, the CCQL discriminator 320(B) shown in FIG. 6 generates the common phase-control signal Zφcj but not the common frequency-control signal Zƒcj. The phase accumulator 802 receives the following inputs:

    
    
        The individual channel frequency-control signal 321 Zƒrj. As described above in reference to FIG. 3B, FIG. 5A, and FIG. 5B, Zƒrj is outputted from the individual loop filter 308-j. Zƒrj is updated at time instants tcr. The period between consecutive time instants tcr is the update period of Tcr; the update frequency is ƒcr=1/Tcr. Zƒrj represents a frequency in units of Hz (number of cycles per time period of one sec). Zƒrj is first inputted into the scale converter 840; the output 841 represents a frequency in units of number of cycles per time period of Tsampl.
        The individual channel phase-control signal 323 Zφrj. As described above in reference to FIG. 3B, FIG. 5A, and FIG. 5B, Zφrj is outputted from the individual loop filter 308-j. Zφrj is updated at time instants tcr. Zφrj represents a phase shift in units of number of cycles.
        The common frequency-control signal 331 Zƒcj. As described above in reference to FIG. 3B and FIG. 4A, Zƒcj is outputted from the CCQL discriminator 320. Zƒcj is updated at time instants tcc. The period between consecutive time instants tcc is the update period of Tcc; the update frequency is ƒcc=1/Tcc Hz). Zƒcj represents a frequency in units of Hz (number of cycles per time period of one sec). Zƒcj is first inputted into the scale converter 850; the output 851 represents a frequency in units of number of cycles per time period of Tsampl.
        The common phase-control signal 333 Zφcj. As described above in reference to FIG. 3B, FIG. 4A, and FIG. 4B, Zφcj is outputted from the CCQL discriminator 320. Zφcj is updated at time instants tcc. Zφcj represents a phase shift in units of number of cycles.
    
    


The value of the (re-scaled) individual channel frequency-control signal 841 Zƒrj is inputted into the sum operator 804 at every time instant tcr; the most recent value of Zƒrj is stored in the sum operator 804 until it is updated. The value of the (re-scaled) common frequency-control signal 851 Zƒcj is inputted into the sum operator 804 at every time instant tcc; the most recent value of Zƒcj is stored in the sum operator 804 until it is updated.
The value of the individual channel phase-control signal 323 Zφrj is inputted into the sum operator 808 at every time instant tcr. The value of Zφrj is not stored but is reset to zero until it is updated. The value of the common channel phase-control signal 333 Zφcj is inputted into the sum operator 810 at every time instant tcc. The value of Zφcj is not stored but is reset to zero until it is updated.
The operation of the phase accumulator 802 proceeds as follows. To simplify the discussion, assume that the phase accumulator has already operated for a time t, where t>Tcr and t>Tcc. At a particular time instant ts,i, where i is an arbitrary integer index, the most recent value of the (re-scaled) individual channel frequency-control signal 841 Zƒrj and the most recent value of the (re-scaled) common frequency-control signal 851 Zƒcj are inputted into the sum operator 804. The sum 821 is inputted into the sum operator 806; the most recent value of the sum 821 is stored in the sum operator 806 and changes only when a new value of Zƒrj or Zƒcj is received. The value of ΦNCOj outputted from the register 812 at the previous time instant ts,i−1 is also inputted into the sum operator 806. The sum 823 is inputted into sum operator 808.
If ts,i coincides with a time instant tcr, then the current value of the individual channel phase-control signal 323 Zφrj is inputted into the sum operator 808. If ts,i does not coincide with a time instant tcr, then a value of zero is inputted into the sum operator 808. The sum 825 is inputted into the sum operator 810. If ts,i coincides with a time instant tcc, then the current value of the common channel phase-control signal 333 Zφcj is inputted into the sum operator 810. If ts,i does not coincide with a time instant tcc, then a value of zero is inputted into the sum operator 810. The sum 827 is then recorded in the register 812 to generate the phase.
The register 812, however, has a maximum value of 2π (1 cycle=2π); that is, the register 812 performs addition modulo 2π. The value stored in the register at ts,i then does not correspond to the (cumulative) full phase but to the fractional part of the full phase: the output 829 of the register 812 is then ΦNCOj, which is then converted to the RefSigj 351, as described above.
The frequency shift of the quartz crystal oscillator caused by external effects (shock and vibration) has a stronger impact on the local oscillator frequency ƒLO than on clock frequencies, such as ƒsampl. For example, typical values of various operating frequencies are the following:

    
    
        Quartz crystal oscillator frequency, ƒq=2×107 Hz.
        Received RF frequency, ƒRF=1.5×109 Hz.
        Local oscillator frequency, ƒLO=1.49×109 Hz.
        Intermediate frequency, ƒIF=1×107 Hz.
        Sampling frequency, ƒsampl=5×107 Hz.

Assume, due to external effects, such as heavy shock, that the shift in the quartz crystal oscillator frequency is 2 Hz. The resulting shift in the local oscillator frequency and the resulting shift in the intermediate frequency is 150 Hz, which is much wider than the PLL tracking bandwidth; consequently, tracking is lost. The resulting shift in the sampling frequency is 5 Hz and does not affect the tracking process.

    
    


Examples of mean values for signal generation frequencies (clock frequencies) are listed in Table 1 below.










TABLE 1







Signals
Signal Generation Frequency









Ij, Qj
about 1000 Hz



Ij*, Qj*
about 200 Hz



Zdrj
about 200 Hz



Zfrj, Zφrj
about 200 Hz



fsampl
from about 40 MHz




to about 60 MHz



Zfcj, Zφcj
from about 700 Hz




to about 1000 Hz










The foregoing Detailed Description is to be understood as being in every respect illustrative and exemplary, but not restrictive, and the scope of the invention disclosed herein is not to be determined from the Detailed Description, but rather from the claims as interpreted according to the full breadth permitted by the patent laws. It is to be understood that the embodiments shown and described herein are only illustrative of the principles of the present invention and that various modifications may be implemented by those skilled in the art without departing from the scope and spirit of the invention. Those skilled in the art could implement various other feature combinations without departing from the scope and spirit of the invention.

## Claims

1. A method for improving the performance of a navigation receiver operably coupled to an antenna, wherein the navigation receiver and the antenna are subjected to vibration and shock, and wherein the navigation receiver comprises a quartz crystal oscillator generating a quartz crystal oscillator signal having a quartz crystal oscillator frequency, the method comprising the steps of:
receiving a plurality of navigation signals, wherein each individual navigation signal in the plurality of navigation signals corresponds to an individual global navigation satellite in a plurality of global navigation satellites, and wherein each individual navigation signal in the plurality of navigation signals has a corresponding radiofrequency;
generating a local oscillator signal having a local oscillator frequency, wherein the local oscillator signal is based on the quartz crystal oscillator signal, and wherein the local oscillator frequency is based on the quartz crystal oscillator frequency;
generating a plurality of intermediate signals based on the plurality of navigation signals and based on the local oscillator signal, wherein each individual intermediate signal in the plurality of intermediate signals corresponds to an individual navigation signal in the plurality of navigation signals, and wherein each individual intermediate signal in the plurality of intermediate signals has a corresponding intermediate frequency based on a difference between the corresponding radiofrequency and the local oscillator frequency;
generating a plurality of sequences of digital signals, wherein each individual sequence of digital signals in the plurality of sequences of digital signals corresponds to a digital representation of an individual intermediate signal in the plurality of intermediate signals; and
for each individual sequence of digital signals in the plurality of sequences of digital signals, tracking a phase of the individual sequence of digital signals with an individual channel phase-lock loop, wherein the step of tracking the phase of the individual sequence of digital signals with an individual channel phase-lock loop comprises the steps of:
generating, with an individual channel correlator operating at a first clock frequency, a sequence of in-phase correlation signals and a sequence of quadrature-phase correlation signals, wherein the sequence of in-phase correlation signals and the sequence of quadrature-phase correlation signals are based on the sequence of individual digital signals and based on a sequence of individual channel reference signals received from an individual channel numerically-controlled oscillator;
generating, with a common coordinate-quartz loop discriminator operating at a second clock frequency, a sequence of common control signals;
generating, with an individual channel accumulator operating at a third clock frequency, a sequence of accumulated in-phase correlation signals and a sequence of accumulated quadrature-phase correlation signals, wherein the sequence of accumulated in-phase correlation signals is based on the sequence of in-phase correlation signals, and wherein the sequence of accumulated quadrature-phase correlation signals is based on the sequence of quadrature-phase correlation signals;
generating, with an individual channel discriminator operating at the third clock frequency, a sequence of phase-lock loop tracking error signals based on the sequence of accumulated in-phase correlation signals and based on the sequence of accumulated quadrature-phase correlation signals;
generating, with an individual channel loop filter operating at the third clock frequency, a sequence of individual channel control signals; and
generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at a fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual channel control signals; and
transmitting the sequence of individual channel reference signals to the individual channel correlator.

2. The method of claim 1, wherein:
the corresponding radiofrequency for each individual navigation signal in the plurality of navigation signals is the same radiofrequency; and
the corresponding intermediate frequency for each individual intermediate signal in the plurality of intermediate signals is the same intermediate frequency.

3. The method of claim 1, wherein:
the corresponding radiofrequency for each individual navigation signal in the plurality of navigation signals is a different individual radiofrequency; and
the corresponding intermediate frequency for each individual intermediate signal in the plurality of intermediate signals is a different individual intermediate frequency.

4. The method of claim 1, wherein the sequence of common control signals is selected from the group consisting of:
a sequence of common phase-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;
a sequence of common frequency-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals; and
a sequence of common phase-control signals and a sequence of common frequency-control signals, wherein the sequence of common phase-control signals and the sequence of common frequency-control signals are based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals; and
wherein the sequence of individual channel control signals is selected from the group consisting of:
a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals;
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals and the sequence of individual channel phase-control signals are based on the sequence of phase-lock loop tracking error signals;
a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals; and
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals is based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals, and wherein the sequence of individual channel phase-control signals is based on the sequence of phase-lock loop tracking error signals.

5. The method of claim 4, wherein the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at a fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:
the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at the fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on a sequence of control signals selected from the group consisting of:
the sequence of individual channel frequency-control signals and the sequence of common phase-control signals;
the sequence of individual channel frequency-control signals and the sequence of common frequency-control signals;
the sequence of individual channel frequency-control signals, the sequence of common phase-control signals, and the sequence of common frequency-control signals;
the sequence of individual channel frequency-control signals, the sequence of individual channel phase-control signals, and the sequence of common phase-control signals;
the sequence of individual channel frequency-control signals, the sequence of individual channel phase-control signals, and the sequence of common frequency-control signals; and
the sequence of individual channel frequency-control signals, the sequence of individual channel phase-control signals, the sequence of common phase-control signals, and the sequence of common frequency-control signals.

6. The method of claim 1,
wherein the sequence of common control signals comprises:
a sequence of common phase-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;

wherein the sequence of individual channel control signals comprises:
a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals; and

wherein the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at a fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:
the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at the fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common phase-control signals and based on the sequence of individual channel frequency-control signals.

7. The method of claim 1,
wherein the sequence of common control signals comprises:
a sequence of common phase-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;

wherein the sequence of individual channel control signals comprises:
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals is based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals, and wherein the sequence of individual channel phase-control signals is based on the sequence of phase-lock loop tracking error signals; and

wherein the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at a fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:
the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at the fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common phase-control signals, based on the sequence of individual channel frequency-control signals, and based on the sequence of individual channel phase-control signals.

8. The method of claim 1, wherein:
wherein the sequence of common control signals comprises:
a sequence of common phase-control signals and a sequence of common frequency-control signals, wherein the sequence of common phase-control signals and the sequence of common frequency-control signals are based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;

wherein the sequence of individual channel control signals comprises:
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals is based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase control signals, and wherein the sequence of individual channel phase-control signals is based on the sequence of phase-lock loop tracking error signals; and

wherein the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at a fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:
the step of generating, with the individual channel numerically-controlled oscillator, wherein the individual channel numerically-controlled oscillator operates at the fourth clock frequency, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common phase-control signals, based on the sequence of common frequency-control signals, based on the sequence of individual channel phase-control signals, and based on the sequence of individual channel frequency-control signals.

9. A navigation receiver operably coupled to an antenna, the navigation receiver comprising:
a low-noise amplifier configured to receive a plurality of navigation signals, wherein each individual navigation signal in the plurality of navigation signals corresponds to an individual global navigation satellite in a plurality of global navigation satellites, and wherein each individual navigation signal in the plurality of navigation signals has a corresponding radiofrequency;
a quartz crystal oscillator configured to generate a quartz crystal oscillator signal having a quartz crystal oscillator frequency;
a local oscillator configured to generate a local oscillator signal having a local oscillator frequency, wherein the local oscillator signal is based on the quartz crystal oscillator signal, and wherein the local oscillator frequency is based on the quartz crystal oscillator frequency;
a radiofrequency mixer configured to generate a plurality of intermediate signals based on the plurality of navigation signals and based on the local oscillator signal, wherein each individual intermediate signal in the plurality of intermediate signals corresponds to an individual navigation signal in the plurality of navigation signals, and wherein each individual intermediate signal in the plurality of intermediate signals has a corresponding intermediate frequency based on a difference between the corresponding radiofrequency and the local oscillator frequency;
an analog-to-digital converter configured to generate a plurality of sequences of digital signals, wherein each individual sequence of digital signals in the plurality of sequences of digital signals corresponds to a digital representation of an individual intermediate signal in the plurality of intermediate signals; and
a digital channel processing unit comprising:
a plurality of individual channel correlators configured to operate at a first clock frequency;
a common coordinate-quartz loop discriminator configured to operate at a second clock frequency;
a plurality of individual channel accumulators configured to operate at a third clock frequency;
a plurality of individual channel discriminators configured to operate at the third clock frequency;
a plurality of individual channel loop filters configured to operate at the third clock frequency; and
a plurality of individual channel numerically-controlled oscillators configured to operate at a fourth clock frequency;

wherein the digital processing unit is configured to receive, from the analog-to-digital converter, the plurality of sequences of digital signals; and
wherein the digital processing unit is configured to, for each individual sequence of digital signals in the plurality of sequences of digital signals, perform the steps of:
generating, with an individual channel correlator in the plurality of individual channel correlators, a sequence of in-phase correlation signals and a sequence of quadrature-phase correlation signals, wherein the sequence of in-phase correlation signals and the sequence of quadrature-phase correlation signals are based on the sequence of individual digital signals and based on a sequence of individual channel reference signals received from an individual channel numerically-controlled oscillator in the plurality of individual channel numerically-controlled oscillators;
generating, with the common coordinate-quartz loop discriminator, a sequence of common control signals;
generating, with an individual channel accumulator in the plurality of individual channel accumulators, a sequence of accumulated in-phase correlation signals and a sequence of accumulated quadrature-phase correlation signals, wherein the sequence of accumulated in-phase correlation signals is based on the sequence of in-phase correlation signals, and wherein the sequence of accumulated quadrature-phase correlation signals is based on the sequence of quadrature-phase correlation signals;
generating, with an individual channel discriminator in the plurality of individual channel discriminators, a sequence of phase-lock loop tracking error signals based on the sequence of accumulated in-phase correlation signals and based on the sequence of accumulated quadrature-phase correlation signals;
generating, with an individual channel loop filter in the plurality of individual channel loop filters, a sequence of individual channel control signals;
generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual channel control signals; and
transmitting the sequence of individual channel reference signals from the individual channel numerically-controlled oscillator to the individual channel correlator.

10. The navigation receiver of claim 9, wherein:
the corresponding radiofrequency for each individual navigation signal in the plurality of navigation signals is the same radiofrequency; and
the corresponding intermediate frequency for each individual intermediate signal in the plurality of intermediate signals is the same intermediate frequency.

11. The navigation receiver of claim 9, wherein:
the corresponding radiofrequency for each individual navigation signal in the plurality of navigation signals is a different individual radiofrequency; and
the corresponding intermediate frequency for each individual intermediate signal in the plurality of intermediate signals is a different individual intermediate frequency.

12. The navigation receiver of claim 9,
wherein the sequence of common control signals is selected from the group consisting of:
a sequence of common phase-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;
a sequence of common frequency-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals; and
a sequence of common phase-control signals and a sequence of common frequency-control signals, wherein the sequence of common phase-control signals and the sequence of common frequency-control signals are based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals; and

wherein the sequence of individual channel control signals is selected from the group consisting of:
a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals;
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals and the sequence of individual channel phase-control signals are based on the sequence of phase-lock loop tracking error signals;
a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals; and
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals are based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals, and wherein the sequence of individual channel phase-control signals are based on the sequence of phase-lock loop tracking error signals.

13. The navigation receiver of claim 12, wherein the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:
the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on a sequence of control signals selected from the group consisting of:
the sequence of individual channel frequency-control signals and the sequence of common phase-control signals;
the sequence of individual channel frequency-control signals and the sequence of common frequency-control signals;
the sequence of individual channel frequency-control signals, the sequence of common phase-control signals, and the sequence of common frequency-control signals;
the sequence of individual channel frequency-control signals, the sequence of individual channel phase-control signals, and the sequence of common phase-control signals;
the sequence of individual channel frequency-control signals, the sequence of individual channel phase-control signals, and the sequence of common frequency-control signals; and
the sequence of individual channel frequency-control signals, the sequence of individual channel phase-control signals, the sequence of common phase-control signals, and the sequence of common frequency-control signals.

14. The navigation receiver of claim 9,
wherein the sequence of common control signals comprises:
a sequence of common phase-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;

wherein the sequence of individual channel control signals comprises:
a sequence of individual channel frequency-control signals based on the sequence of phase-lock loop tracking error signals; and

wherein the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:
the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common phase-control signals and based on the sequence of individual channel frequency-control signals.

15. The navigation receiver of claim 9,
wherein the sequence of common control signals comprises:
a sequence of common phase-control signals based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;

wherein the sequence of individual channel control signals comprises:
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals are based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase-control signals, and wherein the sequence of individual channel phase-control signals are based on the sequence of phase-lock loop tracking error signals; and

wherein the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:
the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common phase-control signals, based on the sequence of individual channel frequency-control signals, and based on the sequence of individual channel phase-control signals.

16. The navigation receiver of claim 9,
wherein the sequence of common control signals comprises:
a sequence of common phase-control signals and a sequence of common frequency-control signals, wherein the sequence of common phase-control signals and the sequence of common frequency-control signals are based on the sequences of in-phase correlation signals and the sequences of quadrature-phase correlation signals generated for all of the sequences of digital signals in the plurality of sequences of digital signals;

wherein the sequence of individual channel control signals comprises:
a sequence of individual channel frequency-control signals and a sequence of individual channel phase-control signals, wherein the sequence of individual channel frequency-control signals are based on the sequence of phase-lock loop tracking error signals and based on the sequence of common phase control signals, and wherein the sequence of individual channel phase-control signals are based on the sequence of phase-lock loop tracking error signals; and

wherein the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common control signals and based on the sequence of individual control signals comprises:

the step of generating, with the individual channel numerically-controlled oscillator, the sequence of individual channel reference signals, wherein the sequence of individual channel reference signals is based on the sequence of common phase-control signals, based on the sequence of common frequency-control signals, based on the sequence of individual channel phase-control signals, and based on the sequence of individual channel frequency-control signals.