This disclosure document relates to a multi-constellation, multi-frequency GNSS system for interference mitigation for a satellite navigation receiver.
The electromagnetic spectrum is limited for wireless communications. As wireless communications are engineered to support greater data transmission throughput for end users, the potential for interference to satellite navigation receivers tend to increase. Interference may be caused by various technical factors, such as inadequate frequency spacing or spatial separation between wireless transmitters, intermodulation distortion between wireless signals, receiver desensitization, or deviation from entirely orthogonal encoding of spread-spectrum signals, outdated radio or microwave frequency propagation modeling of government regulators, among other things. Accordingly, there is need to ameliorate interference through an interference mitigation system.
In accordance with one embodiment, a receiver system with interference mitigation, the receiver system comprises an antenna for receiving a radio frequency signal. A primary downconverter is configured to convert the received electromagnetic signal (e.g., radio frequency signal) to an analog intermediate frequency signal. An analog-to-digital converter is configured to convert the received radio frequency signal to a digital intermediate signal. A wideband interference mitigation module is coupled to an output of the primary downconverter to process the digital intermediate frequency signal. A selective filtering module is associated with a secondary downconverter that comprises a digital harmonic-resistant translator. The selective filtering module comprises: (a) a low-pass filter that is configured as an anti-aliasing digital filter consistent with a target receive bandwidth to suppress aliasing associated with the analog-to-digital conversion, and (b) narrow band rejection filter to filter the digital baseband signal to reduce or to mitigate electromagnetic interference, where the narrow band rejection filter is configured for adaptive control responsive to detection by the wideband interference mitigation module of certain interference in the received radio frequency signal.
As used in this document, adapted to, arranged to or configured to means that one or more data processors, logic devices, digital electronic circuits, delay lines, or electronic devices are programmed with software instructions to be executed, or are provided with equivalent circuitry, to perform a task, calculation, estimation, communication, or other function set forth in this document.
An electronic data processor means a microcontroller, microprocessor, an arithmetic logic unit, a Boolean logic circuit, a digital signal processor (DSP), a programmable gate array, an application specific integrated circuit (ASIC), or another electronic data processor for executing software instructions, logic, code or modules that are storable in any data storage device.
As used in this document, a radio frequency signal comprises any electromagnetic signal or wireless communication signal in the millimeter frequency bands, microwave frequency bands, ultra-high-frequency bands, or other frequency bands that are used for wireless communications of data, voice, telemetry, navigation signals, and the like.
The receiver system 100 represents an illustrative example of one possible reception environment for a radio receiver such as a global navigation satellite system (GNSS) receiver. The satellite 101 (e.g., satellite vehicle) transmits the satellite signal 102 on multiple frequencies, such that the collective set of signals may be referred to as a composite signal. For example, in
Precise point positioning (PPP) includes the use of precise satellite orbit and clock corrections provided wirelessly via correction data, rather than through normal satellite broadcast information (ephemeris and clock data) that is encoded on the received satellite signals, to determine a relative position or absolute position of a mobile receiver. PPP may use correction data that is applicable to a wide geographic area. Although the resulting positions can be accurate within a few centimeters using state-of-the-art algorithms, conventional precise point positioning can have a long convergence time of up to tens of minutes to stabilize and determine the float or integer ambiguity values necessary to achieve the purported (e.g., advertised) steady-state accuracy. Hence, such long convergence time is typically a limiting factor in the applicability of PPP.
In accordance with one embodiment,
Here, a dual band system is described as an example with a first analog signal path within the first analog module 111 corresponding to the first radio frequency signal (e.g., a GNSS channel or set of GNSS channels) and a second analog signal path of a second analog module 131 associated with second radio frequency signal (e.g., a GNSS channel or set of GNSS channels), although in other configurations multiple parallel signal paths for corresponding different frequency bands may be used. For example, a dual band system includes a low-band and a high band, where the low-band has a lower frequency range than the high band does. For the Global Positioning System (GPS), the transmitted L1 frequency signal of a satellite 101 may comprise the high band; the transmitted L2 frequency signal may comprise the low band signal. Further, the L1 carrier is at 1,575.42 MHz, which is modulated with the P(Y) code (pseudo random noise code) and M code that occupies a target reception bandwidth on each side of the carrier. Meanwhile, the L2 carrier is at 1,227.6 MHz and modulated with the C/A (coarse acquisition) code, P(Y) code (pseudo random noise code) and M code that occupies a target reception bandwidth on each side of the carrier. The splitter 107 (e.g., diplexer) splits the composite signal into the first signal path (e.g., upper signal path or high-band path) and the second signal path (e.g., lower path or low-band path).
In one embodiment, the signal splitter 107 or hybrid may split the received signal into two received radio frequency signals for processing by a first analog module 111 and a second analog module 131. Each analog module (111, 131) converts the received radio frequency (RF) or microwave frequency signal to an intermediate frequency (IF) signal or a quasi-baseband signal, which means a baseband signal that is equal to or greater than the zero frequency or direct current (DC) signal in the frequency domain. For example, the intermediate frequency (IF) signal or a quasi-baseband signal may fall within a certain frequency range that excludes zero frequency or direct current DC signal, such as a certain frequency range that has a lower limit that is equal to or greater than approximately 1 Megahertz (MHz).
The first analog module 111 may comprise an optional first pre-amplifier 141 or a low-noise amplifier (LNA) for amplifying the received signal. Similarly, the second analog module 131 may comprise an optional second pre-amplifier 151 or a low-noise amplifier (LNA) for amplifying the received signal. To simplify the receiver analog filtering design, the front-end of a typical modern GNSS receiver uses a wideband front-end design to receive multiple GNSS signals using two/three wideband filters (not shown), where each band targets a target bandwidth (e.g., 140-300 MHz).
In the first analog module 111, a first downconverter 142 is configured to convert the (amplified) first radio frequency signal to an analog intermediate frequency signal. For example, the first analog module 111 comprises a first downconverter 142, such as the combination of a mixer 156 and a local oscillator 154, which moves or shifts the high band (L1, G1, B1, or similar frequency associated with a GNSS) radio frequency (RF) into the analog intermediate frequency (IF). The first downconverter 142 is coupled to the first analog-to-digital (ADC) converter 112. For example, the first downconverter 142 may be coupled to the first ADC 112 via an optional first automatic gain control 143 and an optional first anti-aliasing filter 158, such as an analog low pass filter. Further, the optional first anti-aliasing filter 158 may be configured as an analog low pass filter or analog bandpass filter (BPF) that rejects the image band of the mixer output from the first downconverter 142.
In the first analog module 111, an optional first automatic gain control (AGC) 143 is coupled to the first ADC 112 and the first downconverter 142. The optional first AGC 143 is illustrated in dashed lines to indicate that it is optional and can be deleted in certain embodiments. For example, in one configuration of the first signal path, an optional first automatic gain control (AGC) 143 is coupled to the first ADC 112, the first downconverter 142, and the first pre-amplifier 141. The optional first automatic gain control (AGC) 143 may control the gain (e.g., root mean square (RMS) amplitude) of an input signal to the corresponding first analog-to-digital converter (ADC) 112 to be constant or within a target range (e.g., despite fluctuations in ambient radio frequency noise and the interference signal 103). The first AGC 143 receives gain-related feedback (as indicated by the dashed arrow) from the first ADC 112 to adjust the gain setting of first downconverter 142 (and/or the first pre-amplifier 141).
In the second analog signal path 157, a second downconverter 152 is configured to convert the (amplified) second radio frequency signal to an analog intermediate frequency signal. For example, the second analog module 131 comprises a second downconverter 152, such as the combination of a mixer 166 and a local oscillator 164, which moves the low band (L2, or similar frequency associated with a GNSS) radio frequency (RF) into the intermediate frequency (IF). In one configuration, the second downconverter 152 may be coupled to the second ADC converter 132 via an optional second automatic gain control 153 and an optional second anti-aliasing filter 168, such as an analog low pass filter. Further, the optional second anti-aliasing filter 168 may be configured as an analog low pass filter or analog bandpass filter (BPF) that rejects the image band of the mixer output from the second downconverter 152.
As used in this document, each of the first downconverter 142 and the second downconverter 152 is individually referred to as a primary downconverter; the first downconverter 142 and second downconverter 152 are collectively referred to as primary downconverters. Meanwhile, as used in this document, each of the first harm resistant frequency translator 201 (
In the second analog signal path 157, an optional second automatic gain control (AGC) 153 is coupled to the second analog-to-digital converter (ADC) 132 and the second downconverter 152. As indicated by the dashed lines in
Each analog-to-digital converter (ADC) (112, 132) may be coupled to its optional corresponding automatic gain control (AGC) (143, 153) that provides variable gain amplification. In turn, each optional AGC is coupled to its corresponding downconverter (142, 152). In one embodiment, the automatic gain control AGC provides a feedback signal to the downconverter (142, 152) or intermediate frequency (IF) filter (e.g., analog IF filter) that is associated with the downconverter. The downconverter (142, 152) or its analog IF filter adapts the signal voltage (pea-to-peak) within the first ADC 112 to be commensurate with its operational range.
In one embodiment, each ADC (112, 132) samples the analog received signal (from the corresponding downconverter (142, 152)) using a predefined sampling rate, which, per Nyquist theorem, should be equal to or greater than two times the bandwidth (e.g., target reception bandwidth) for the practical sampling design. In an alternate embodiment, each ADC (112, 132) samples the analog signal using a predefined sampling rate, which, per Nyquist theorem, should be greater than the one times bandwidth for complex signal sampling and two times bandwidth for the real signal sampling. The bandwidth of the ADC determines maximum tolerable interference at a given quantization loss. The resulting digital sequence or filter input of digital signal (113, 133) (e.g., intermediate frequency (IF) signal) reconstructs the received signal, such as the first signal (e.g., high-band RF signal) and the second signal (e.g., low-band RF signal) to the digital intermediate frequency signal (or alternately a quasi-digital baseband signal or digital baseband signal) with a corresponding baseband bandwidth or range.
Quantization refers to the sampling of an analog signal, such as a received satellite signal, by an analog-to-digital converter to produce a number of finite or discrete digital values (e.g., quantum level) per unit time. Interference or aliasing from electromagnetic interference can cause quantization distortion or quantization noise, where quantization distortion can represent a sampling error during the analog to digital conversion and where quantization noise can represent a difference between the digital values (e.g., quantum level) and the actual amplitude, phase or both of sampled analog signal (e.g., for analog signals of weak amplitude or in the presence of multipath or electromagnetic interference).
Aliasing sometimes refers to a false signal, misleading signal, or processing artifact of an ADC (112, 132) or RF-sampling ADC that uses a sampling rate that is too low relative to the frequency of the received signal, such as received satellite signal at a carrier frequency or fundamental frequency. The sampling rate should be selected in accordance with Nyquist theorem and/or Shannon's sampling theorem to avoid aliasing (e.g., in the absence of electromagnetic interference signals). For example, Shannon's sampling theorem provides that the analog signal must be sampled at a rate or frequency that is at least twice the frequency of the highest components contained within received signal, such as the carrier frequency and the full modulation bandwidth (e.g., at least half power bandwidth or minimum percent (e.g., 95 percent) of transmitted energy or mean transmitted power for spread spectrum signals transmitted between a lower frequency limit and an upper frequency limit) with respect to the carrier frequency.
In general, an alias signal is a false digital signal that represents a subharmonic or harmonic component of the true digital signal. In addition to improper selection of the sampling frequency, aliasing can result from an interference signal that is processed by an analog to digital converter, particularly in the absence of proper interference mitigation filtering.
The ADC (112, 132) inherently produces aliasing based on the following factors: the received frequency of the received satellite signal (e.g., with a corresponding fundamental frequency or carrier frequency), one or more harmonics of frequency of analog received satellite signal (e.g., second harmonic representing two multiplied by the fundamental frequency and the third harmonic representing three multiplied by the fundamental frequency), the sampling frequency (e.g., greater than two times the received satellite frequency or fundamental frequency) of the analog to digital converter, the reference frequency of the mixer, the modulation and bandwidth (e.g., half-power bandwidth) of the received satellite signal, and filter attenuation of filters or decimating filters that filter the analog received satellite signal and its harmonics. Prior to the ADC (112, 132), an optional analog anti-aliasing filter or low pass filter (158, 168) may be tuned or configured to attenuate the alias signals, such as one or more harmonics of analog received satellite signal, or the reference frequency of the mixer.
In terms of the AGC feedback control from its respective ADC (112, 132), the AGC feedback control can be done either in analog or digital domain. For example, an envelope detector is typically used for AGC and variable gain control if analog control is used. Because of advances in digital processing theory and practice, the digital processing for the AGC feedback control, can be based on a statistical processes, such as digital analysis of a histogram of the sample digital stream or filter inputs of digital signal (113, 133) at the output of the corresponding analog-to-digital converter, ADC (112, 132) to generate a feedback signal to control the AGC (143, 153), such as the first AGC 143 associated with the corresponding first analog signal module 111 and the second AGC 153 associated with the second analog signal path 157. Each AGC is coupled to the downconverter (142, 152), which in practice may comprise the downconverter and IF filter module with inherent gain/amplification adjustments.
In an alternate embodiments, the analog-to-digital converter (ADC) (e.g., first ADC, second ADC, or both) may comprise a radio frequency ADC or RF-sampling ADC that replaces some components of the analog module, such as the local oscillator (e.g. numerically controller oscillator), mixer, any analog intermediate frequency amplifier, and analog band pass filter. Further, in certain alternate embodiments, the analog ADC may comprise an RF ADC followed by an integral digital-downconverter (DDC) with a digital local oscillator (NCO), digital mixer and digital baseband filter. For example, the DDC comprises a tuning NCO that provides a reference signal to one or more digital mixers, which can be filtered by a digital low pass filter, a digital bandpass filter, or both, such as one or more finite-impulse response filters (e.g., for low pass, bandpass, each filter having decimation rate by n), amplified (optional) and down-sampled).
WBI Mitigation Module
A first analog-to-digital converter 112 is configured to convert the analog intermediate frequency signal to a digital intermediate frequency signal, or to convert an analog baseband signal to a digital quasi-baseband signal. A first wide band interference mitigation module 117 is coupled to a digital output of the first analog-to-digital filter 112. The first WBI mitigation module 117 is configured to mitigate or attenuate WBI in the digital baseband signal or quasi-baseband signal. The first WBI mitigation module 117 detects the short period ADC saturation resulting from the WBI and processes or blanks those saturation period to mitigate its impact. Blanking refers to a process of disabling a signal, data stream or channel for one or more sampling intervals, such as any of the following: (a) retaining (e.g., storing and holding) a previous value of the disabled signal for a corresponding sampling interval; (b) discarding or rejecting a sampled value for a sampling interval; (c) discarding and averaging a sampled value over discarded and adjacent sampling time intervals, and/or (d) assigning a null value or predefined logic value to digital signal for a discarded sampling interval.
A second analog-to-digital converter 132 is configured to convert the analog intermediate frequency signal to a digital intermediate frequency signal, or to convert an analog baseband signal to a digital quasi-baseband signal. A second wide band interference mitigation module 137 is coupled to a digital output of the second analog-to-digital filter 132. The second WBI mitigation module 137 is configured to mitigate or attenuate WBI in the digital baseband signal or quasi-baseband signal. The second WBI mitigation module 137 detects the short period ADC saturation resulting from the WBI and processes or blanks those saturation period to mitigate its impact.
The first selective filtering module 114 extracts a targeted component, such as selected band, an upper band, a high band or a sub-band, from the inputted digital intermediate frequency signal, or digital (upper band spectrum) signal 113. For example, the first selective filtering module 114 may comprise a bandpass filter for a band component or sub-band component of the received satellite GNSS signal in its derivative form of the digital intermediate frequency signal or quasi-base band signal. The first selective filtering module 114 outputs a resultant digital baseband signal, which in one configuration comprises a GNSS signal at a band of interest (e.g. L1 or G1 or B1), the in-GNSS-band NBI, and the noise.
In one embodiment, as illustrated in
The first harmonic-resistant frequency translator 201 or secondary downconverter represents a digital downconverter, which can operate in the digital frequency domain, to process, mix or shift an inputted digital intermediate frequency signal (e.g., quasi baseband signal) to a digital baseband signal. The output of the first harm resistant frequency translator is inputted to a first digital low pass filter such as a band low pass filter or bandpass low pass filter that passes the upper band or high band of the received satellite signal (e.g., received GNSS signal). The first digital low pass filter or bandpass filter (BPF) is configured to reject the image band (e.g., to reduce or attenuate aliased signal components) of the output or mixer output component from the secondary downconverter.
The second selective filtering module 134 extracts a targeted component, such as selected band, an low band, a low band or a sub-band, from the inputted digital intermediate frequency signal, from the inputted digital intermediate frequency signal, or digital (lower band spectrum) signal 113. For example, the second selective filtering module 134 may comprise a bandpass filter for a band component or sub-band component of the received satellite GNSS signal in its derivative form of the digital intermediate frequency signal or quasi-base band signal. The second selective filtering module 134 outputs a resultant digital baseband signal, which in one configuration comprises a GNSS signal at a band of interest (e.g. L2 or G2 or B2), the in-GNSS-band NBI, and the noise.
In one embodiment, as illustrated in
The second harmonic-resistant frequency translator 211 or secondary downconverter represents a digital downconverter, which can operate in the digital frequency domain, to process, mix or shift an inputted digital intermediate frequency signal (e.g., quasi baseband signal) to a digital baseband signal. The output of the second harm resistant frequency translator 211 is inputted to a second digital low pass filter 213 such as a band low pass filter or bandpass low pass filter that passes the lower band or low band of the received satellite signal (e.g., received GNSS signal). The second digital low pass filter 213 or bandpass filter (BPF) is configured to reject the image band (e.g., to reduce or attenuate aliased signal components) of the output or mixer output component from the secondary downconverter.
The output of the first selective filtering module 114 is coupled to the input of the first narrowband interference (NBI) mitigation module 110. The first NBI system 110 is configured to reject an interference component that interferes with the received radio frequency signal (e.g., associated with a first sub-band or set of channels). For example, in one embodiment the first NBI mitigation system 110 comprises a first narrow band notch filter 205, such as a first adaptive narrow band notch filter that has a dynamically adjustable attenuation versus frequency response, a selectable attenuation versus frequency response (e.g., notch depth, notch width), or a library of preset frequency responses to dynamically select from. The first narrow band notch filter 205 can be selectively controlled by a controller or the output of the first wide band interference mitigation module 117 to temporarily disable the notch filter for sampling intervals or samples that correspond to blanked digital intermediate frequency signal components.
The output of the second selective filtering module 134 is coupled to the input of the second narrowband interference (NBI) mitigation system 130. The second NBI system 130 is configured to reject an interference component that interferes with the received radio frequency signal (e.g., associated with a first sub-band or set of channels). For example, in one embodiment the second NBI mitigation system 130 comprises a second narrow band notch filter 215, such as a first adaptive narrow band notch filter that has a dynamically adjustable attenuation versus frequency response, a selectable attenuation versus frequency response, or a library of preset frequency responses to dynamically select from. The second narrow band notch filter 215 can be selectively controlled by a controller or the output of the second wide band interference mitigation module 137 to temporarily disable the notch filter for sampling intervals or samples that correspond to blanked digital intermediate frequency signal components.
The first digital automatic gain control (DAGC) 207 is configured to adjust the gain to compensate for attenuation or amplification of the magnitude of received signal components, or corresponding samples, of the satellite signal. The second digital automatic gain control (DAGC) 217 is configured to adjust the gain to compensate for attenuation or amplification of the magnitude of received signal components, or corresponding samples, of the satellite signal.
To the extent a relatively strong in-GNSS-band NBI component is present in digital baseband signal 115, the GNSS receiver will tend to experience signal to noise ratio (SNR) degradation; such degradation is significantly determined by the relative location of NBI reference to the pseudorandom noise (PN) like signal in the frequency domain and the de-spreading gain that a specific PN sequence provides. For a GNSS satellite that transmits a spread spectrum signal (e.g., like code division multiple access modulation) the de-spreading gain or spread spectrum processing gain is the ratio of the spread bandwidth or total bandwidth of the electromagnetic signal (e.g., frequency spectrum of modulated radio frequency or microwave signal that contains approximately 95 percent of the transmitted energy) to the baseband bandwidth. To mitigate the impact of the in-GNSS-band NBI on the PN sequence demodulation performance, the NBI mitigation system (110, 130) is configured to adaptively reject the NBI in the received signal or digital baseband. The residual or resultant signal (116, 136) ideally, contains only the PN signal and the noise.
In one embodiment, the blanking process in WBI mitigation module (117, 137) introduces the phase discontinuity which negatively impact the stability and accuracy of the NBI mitigation system (110, 130); therefore the blanking enable signal (119, 139(passes onto the NBI mitigation system (110, 130) to make it pause the adaptive updating during the blanking period.
In
As shown in
The selected channel or appropriate sample stream will be further processed by the GNSS channel processing module 121, which typically comprises one or multiple carrier phase demodulators, the PN code generator sampled at multiple delayed phase/time offsets (e.g., or at taps of shift registers), the binary offset sub-carrier (BOC) modulator (used for modern GNSS signal such as GPS L1C, Bei Dou B1C, Galileo E1 signals etc.), and multiple accumulators, to create a bank or accumulations of in-phase (I) and quadrature (Q) measurements at an interval (e.g., of millisecond (ms) or multiple milliseconds) to drive the baseband tracking loop (e.g., to support timely, reliable, precise carrier phase and/or code phase measurements).
In one embodiment, the GNSS channel processing module 121 may comprise a baseband tracking loop module (e.g., 711 in
In one embodiment, the navigation, control and interface module 122 takes a pseudo-range measurements and carrier phase measurements and other related information from the satellites 101 to generate the positioning solution, which is used as a feedback to align the receiver crystal-grade clock (e.g., of lower temporal precision) with the satellite-based atomic grade clock (e.g., of higher temporal precision); the solution also, combined with other information, generates the in-view satellites 101 list to control the appropriate receiver resource allocation.
As used in this document, accumulators refer to registers, flip-flops, electronic memory, or other electronic data storage devices that are used to store the results of computations of an electronic data processor, such as a microprocessor, microcontroller, application specific integrated circuit (ASIC), system on chip (SOC), programmable logic array (PLA), field programmable gate array (FPGA), an arithmetic logic unit (ALU), a Boolean logic unit (BLU), digital logic circuits, a digital signal processor or another data processing device. The correction signals derived from the baseband tracking loops control the numerically controlled oscillators (NCOs) in the GNSS channel processing module 121 to maintain the synchronization between the received signal in the channel and the local replica of that channel (e.g., for decoding or demodulation of certain signal components and/or precise, pseudorange measurements, code phase measurements and/or carrier phase measurements).
As used in this document, pseudorange measurements or pseudorange estimates are based on code phase measurements and/or carrier phase measurements between a given GNSS satellite and GNSS receiver. For example, by evaluating a reference symbol or edge (e.g., leading edge, trailing edge or pulse) of received code phase signal with respect to the replica code phase signal, the GNSS rover receiver can estimate the difference between a transmit time (e.g., consistent with GNSS time standard and any available correction for clock bias) at the given GNSS satellite, a reception time (e.g., consistent with GNSS time standard and any available correction for clock bias) at the GNSS receiver, where the difference is multiplied by the speed of light (e.g., approximately 3×108 m/s).
In one embodiment, the navigation, control and interface module 122 takes the pseudorange and carrier phase measurements and other related information from the satellites to generate a positioning solution. In certain configurations, the position solution of the GNSS receiver, or its antenna, can be used as a feedback to align the GNSS receiver crystal-grade clock with the satellite-based atomic grade clock; the solution also, combined with other information, generates the in-view satellites list to control the appropriate GNSS receiver resource allocation.
To produce or operate a reliable interference resistant receiver system, the low pass filtering module (203, 213) is configured to address the potential limitations of the harm-resistant frequency translator, or to minimize the potential deleterious effects of aliasing in the processing (e.g., rotating or mixing of) digital baseband signal. After the analog-to-digital converter (112, 132) in the digital filtering, one or more digital filters, such as low pass filters (203, 213), can attenuate alias signal components, such as the second harmonic and the third harmonic, and the shifted second harmonic and shifted third harmonic which are shifted by the reference frequency of the mixer (e.g., in translator 201, 211), to the extent that the digital intermediate frequency approaches a zero frequency or DC or is spaced apart from the shifted second harmonic and shifted third harmonic. Ideally, the digital low pass filter (203, 213) such as a decimation filter or decimation filters do not affect the target or desired signal components of digital intermediate frequency or baseband, but can a attenuate one or more harmonics. Further, the decimation filter can shift an image or artifact of the received analog frequency to the negative frequency or complex frequency for attenuation, to the extent that the shifted image or artifact lies within (or wraps around to the real frequency domain to support) the attenuation or rejection band of the LPF. If interference is present, the interference can add to the inherent aliasing of the fundamental frequency or its harmonics, or provide narrow-band or wide-band components adjacent or near the fundamental frequency or its harmonics that can be addressed at proper time periods by the filtering. For example, the receiver is configured to apply adaptive and dynamic filtering to corresponding sampling intervals (e.g., epochs) of the received GNSS signal, and its derived digital baseband signal to mitigate certain WBI, NBI, or both.
The order or filter order of a digital filter, such as LPF (203, 213) may represent the maximum delay (e.g., of sampling intervals) that are used to create the filtered output from a digital input, such as the maximum number of samples (e.g., higher number of maximum sampling interval M or N) used in the difference equation in the time domain. The anti jamming capability of the receiver tends to increase of the ADC bit width (e.g., such as from 32 bit data blocks to 64 bit data blocks for faster or greater data processing throughput) of the ADC (112, 132) which consequently increases the logic complexity of the digital processing.
In one example, the low pass filter (203, 213) or other anti-aliasing filter is used at the output of the harmonic-resistant frequency translator (201, 211) to attenuate, reduce or ameliorate aliasing, mixing images, and mixing artifacts that might otherwise degrade performance. In one embodiment, the band low-pass filter (203,213) comprises an anti-aliasing digital filter that ideally provides a substantially flat or substantially uniform amplitude response versus frequency for baseband signal frequency components within the passband and complete attenuation or greater than a minimum threshold attenuation (relative to passband) of the baseband signal frequency components outside of the passband with a linear phase response over the entire passband.
For example, the anti-aliasing filter at the output of the harmonic-resistant frequency translator (201, 211) could be configured as an Nth order digital bandpass filter (BPF), which would require 2N+2 coefficients. However, the anti-aliasing filter at the output of the harmonic-resistant frequency translator (201, 211) is generally configured as the low pass filter (LPF) with equivalent bandwidth to the above BPF that only requires N/2 coefficients due to its symmetric property. To utilize this symmetric property, the harmonic-resistant frequency translator (201, 211) is configured to translate the digital intermediate frequency signal of interest into a digital baseband signal. However, within the harmonic-resistant frequency translator (201, 211), the digital mixer (rotation) is imperfect; the imperfection introduces the harmonics around its fundamental frequency (e.g., frequency of local oscillator input to the digital mixer). Further, the harmonics can introduce alias components from the interference around the harmonics into GNSS processing bandwidth. In order to reduce the logic complexity and processing burden, a harmonic-resistant translator (201, 211) may comprise a harmonic-resistant mixer (201, 211) that is designed to translate the GNSS signal to near zero frequency with sufficient suppression of the harmonic aliasing. The detailed design of the harmonic-resistant mixer is described in conjunction with the block diagram of
After the digital mixing by the harmonic-resistant frequency translator (201, 211), the digital baseband signal (202, 212) is at baseband or quasi-baseband (e.g., at, near, or approaching zero frequency). Therefore, any adjacent-GNSS-band interferences can be attenuated or mitigated using the LPF (203, 213). The resultant digital baseband signal (115, 135) that is outputted by the LPF (203, 213) comprises the in-GNSS-band interference, the GNSS signal, and the accompanying noise.
In one embodiment, a notch filter (205, 215) or another narrow band rejection filter is configured to attenuate an interference component of the baseband or digital intermediate frequency signal, wherein the narrow band rejection filter comprises: an adaptive notch filter and a controller for controlling the magnitude versus frequency response of the adaptive notch filter, where the controller has an estimator to estimate filter coefficients of the adaptive notch filter. For example, each NBI mitigation system (110, 130) comprises a corresponding narrow band notch filter (205, 215). The digital notch filter (205, 215) is configured with appropriate coefficients to produce a filter magnitude versus frequency response to mitigate the impact of the in-GNSS-band NBIs. The notch filter (205, 215) is commonly controlled or driven by an adaptive algorithm in either frequency domain or time domain. Theoretically, an Nth order notch filter (205, 215) can reject up to N interferences. However, the mismatch between the order of the notch filter (205, 215) and the number of the interference components (e.g., at different frequencies and associated interference bandwidths at those frequencies) negatively impact the accuracy.
In one embodiment, the WBI mitigation module (117, 137) uses the blanking (e.g., a blanking method), which can introduce a discontinuity or phase jump in the signal 118 or waveform. Such phase jump can reduce the stability and accuracy of the notch filter (205, 215). Therefore, in practice on a band-by-band basis (e.g., for lower band or higher band, independently) each NBI mitigation system (110, 130) or its notch filter (205, 215) uses a corresponding blanking enable signal (119, 139) to disable the adaptive updates during the period that phase jump is detected for the respective band or sub-band (e.g., lower band or higher band).
For each band or sub-band, the output signal or digital filtered baseband signal (206, 216) from the notch filter (205, 215) comprises the GNSS signal, the noise, and the residual of the interference (e.g., which is supposed to be at minimal strength). Because the filtered baseband signal (206, 216) comprises the GNSS signal and the noise, a low bit-width quantization of the ADC (112, 132) would be sufficient to ensure or promote the quantization loss less than a design threshold (e.g., 0.1 dB for low-bit width, where such design threshold may increase for a greater bit-width quantization). A lesser quantization level tends to reduce the logic complexity for the channel processing. However, the dynamic range of the filtered baseband signal (206, 216) tends to be dramatically different based on the interference signal(s) that are present at any time. For example, the digital signal (113, 133) may comprise strong NBIs at some times, whereas at other times the digital signal (113, 133) is noise-like or has the general characteristics of spread-spectrum modulated signal with PN code. For each band or sub-band, the respective digital AGC (DAGC) (207, 217) can adaptively adjust the magnitude of the respective filtered baseband signal (206, 216) which is susceptible to changes from interference, to yield resultant signal (116, 136) for a band or sub-band at a constant magnitude.
In
With respect to the lower band, similar to the above upper band operation, the second ADC 132 outputs a digital intermediate frequency signal 133 that is processed by the second WBI mitigation system 137. In the selective filtering module 134, the processed, mitigated signal 138 is translated into the digital baseband signal by the second harmonic-resistant frequency translator 211. The second LPF 213 suppress the adjacent-GNSS-band interference in the second digital baseband signal 212. The NBI component of signal 214 at the second path is mitigated through the second notch filter 215. The second notch filter 215 uses the second blanking enable signal 139, to disable the updates during the period of phase jump. The second DAGC 217 adaptively adjust the magnitude of the signal 216 to a constant level, a target level or within a target range.
In
Each harmonic-resistant frequency translator (201, 211) comprises one or more digital harmonic-resistant mixers 303 or digital rotators that receives one or more corresponding digital mixing signals 302, such as a sine or in-phase digitally generated signal and cosine or quadrature digitally generated signal from a look-up table 301. For example, the digital harmonic resistant mixer may comprise a first digital mixer and a second digital mixer, where the first digital mixer is configured to generate an in-phase component (I) of the baseband signal or near-baseband signal and where the second digital mixer is configured to generate an quadrature (Q) component of the (encoded/modulated) baseband frequency signal, where the I component and Q component collectively for a channel, set of channels, or aggregate channel (e.g., super channel) representative of the set of channels for a received GNSS signal or GNSS composite signal contains all, or a majority of, the electromagnetic signal energy.
The first mixer receives a digital (intermediate frequency) signal (113,133), derived from the received GNSS signal, and a local reference intermediate frequency signal, derived from local carrier signal or carrier replica, and generates the I component of the (encoded/modulate) baseband frequency or near-baseband frequency signal. Similarly, the second mixer receives a digital (intermediate frequency) signal (113,133), derived from the received GNSS signal, and a local reference intermediate frequency signal, derived from local carrier signal or carrier replica, and generates the I component of the (encoded/modulate) baseband frequency or near-baseband frequency signal.
In one embodiment, the look-up table comprises a first look-up table for generation of a corresponding sine or in-phase local reference intermediate frequency signal and a second look-up table for generation of a cosine or quadrature-phase (e.g., offset 90 degrees from the in-phase component) local reference intermediate frequency signal. Further, in one configuration, the look-up table 301, which comprises the first look-up table and the second look-up table, is configured to generate one or more high precision, low resolution digital mixing signals 302 whose harmonic content is sufficiently suppressed. For example, the look-up table 301 may model or store the signal in terms of sine component, or a cosine component, or both at one or more fundamental frequencies for one or more integer number of cycles. In the look-up table 301, the stored or modeled signal may be stored as magnitude and phase values in the look-up table, for instance.
In one embodiment, a local intermediate frequency (IF) signal generator, such as a carrier numerically controlled oscillator (NCO) with a frequency scaler, or an IF NCO provides an local IF signal or local reference signal as an input to the digital harmonic resistant mixer 303. Further, local IF frequency generator, IF NCO or carrier NCO (e.g., as well as the code NCO for the same channel, set of channels or aggregate set of channels), may be driven by a clock signal with a frequency that is precisely controlled by a phase locked loop or the clock tracking loop module (e.g., 730 in
In one configuration, a phase accumulator or accumulations for correlation data processing of encoded signal components may be used as an input to the clock tracking loop module (e.g., 730 in
Referring to the charts (352, 354, 356) of
In the uppermost chart 352 of
Here in the example, the strong adjacent-GNSS-band NBI component 312 is too close to the harmonic signal 322 which tends to translate the interference signal of NBI component 312 into or with the digital baseband signal 332. As a result, the demodulation of the GNSS signal in the digital baseband 331 will be disturbed by the baseband interference signal 332. This illustrative scenario exemplifies the importance to suppress the harmonic signal 322 and 323, which in turn suppresses the baseband interference signal 332.
In one configuration, the critical factors to suppress the undesired harmonic component are the length of one cycle (e.g., corresponding to the local oscillator frequency of the carrier offset frequency, intermediate frequency signal of the local oscillator or NCO for digital mixing) and the precision of each element in the look-up table (LUT) 301. The ideal phase rotation with only fundamental line spectrum can be represented by a float format vector Aejωn, where
A is the magnitude of the phasor,
ω is the normalized frequency, and
n is the index of the sample.
Referring to
In a first example of
In a second example of
In one embodiment, the receiver may use a multi-stage down-conversion process, even within the secondary downconverter. In the secondary downconverter, the digital harmonic-resistant mixer can implement a first-stage mixer process and a second stage filtering process. First, the digital harmonic-resistant mixer can use the coarse resolution (shorter cycle), high precision approach to translate the GNSS signal into the near digital baseband or quasi-baseband (e.g., near rather than at ideal 0 frequency) that is located within the passband of LPF (203, 213). As long as the interference aliasing is prevented by the harmonic-resistant frequency translator (201, 211) (e.g., to prevent aliasing that falls outside the bandwidth of the low pass filter, such as the negative frequency or imaginary frequency domain), in the second stage filtering process the LPF (203, 213) is able to reject the adjacent-GNSS-band interference. Following the LPF (203, 213, the notch filter (205, 215) can further take out the in-band, GNSS-band interference. Therefore, it is possible that during certain time periods, the resultant filtered baseband signal (206, 216) doesn't comprise any strong interference relative to the digital baseband signal; during such time periods the harmonic constraint is less strict, where in the first-stage mixer process, the harmonic-resistant frequency translator (201, 211) can operate (e.g., temporarily) with a long cycle (high resolution), such that the mixing can allow relative high harmonic components during such time periods.
In an alternate embodiment, an additional or supplemental digital mixer may follow the LPF (203, 213) to shift or translate any digital quasi-baseband signal to a true digital baseband signal at zero frequency or over a bandwidth of the baseband that is inclusive of a zero frequency component.
In
In
in
Following the low pass filtering of the LPF or bandpass filtering of selective filtering module (114, 134), the notch filter (205,215) either in finite impulse response (FIR), an infinite impulse response (IIR), or hybrid of both the FIR and IIR, is used to mitigate the impact of in-GNSS-band NBI 441. Because of the spectrum overlap between the in-GNSS-band NBI 441 and the desired GNSS signal component 311, the mitigation of the NBI signal 441 is susceptible to distortion (e.g., typically, inevitably distorts) the waveform of the desired GNSS signal component 311, which adversely impact the tracking accuracy. Therefore, a notch filter control module is configured to control the notch filter (205, 215) to facilitate advantageously enabling notch filter for sampling intervals that overweigh the disadvantage from the potential distortion. The blanking algorithm in WBI mitigation module (117, 137) provides a control signal or data message indicative of phase discontinuity which negatively affect the tracking performance of the notch filter (205, 215) for one or more sampling intervals to control, disable and enable the notch filter (205, 215). If a disable notch-filter signal or blanking enable signal (119, 139) is provided to notch filter (205, 215), the adaptive updates of notch filter are delayed, stayed or held off over the period (e.g., series of successive sampling intervals or epochs, which is a measurement time unit within the GNSS receiver) that phase discontinuity happens.
At the output of the notch filter (205, 215), the digital filtered baseband signal (206, 216) only comprises the desired GNSS signal component 311 and the noise. However, the magnitude of signal (206, 216) can be dramatically diversified between the case that in-GNSS-band interference signal 441 is present and the case that interference signal 441 is absent. The two-bit quantization loss to noise-like signal desired GNSS signal component 311 is negligible; therefore to simplify the channel processing logic, the digital AGC (DAGC) (207, 217) is configured to adaptively scale the signal 413 to a constant level, a target signal level, or within a target signal range (e.g., digital compression). The resultant signal (116,136), whose spectrum can be distorted from interference without proper or reliable interference mitigation, drives the channel tracking loops to produce the navigation solution.
The second case starts with the input signal 531 (for digital filtering), which comprises the GNSS signal (e.g., 311) encoded with the desired PN signal, either the adjacent-GNSS-band NBI signal (e.g., NBI component 312 in
In order to achieve the target gain control goal of DAGC system (207, 217),
A first comparator 547 compares the resultant signal 544 with a desired threshold signal 546. If the first comparator 547 determines that signal 544 is greater than the signal 546 and that the signal 545 is greater than 1, the counter 541 decrements or decreases the scale signal (S) 545 in response to the decrement enable input 543 that is provided by the firsts comparator 547.
The second comparator 548 compares resultant signal 544 with desired signal 546. If the second comparator 548 determines that the signal 544 is less than the signal 546 and if the scale signal (S) 545 does not exceeds the maximum value (e.g., 2X−1) of the counter 541, the second comparator 548 generates an increment signal 542, which enables the counter 541 to increase or increment the scale signal 545. Otherwise, the signal 545 maintained constant and the signal 544 stays at the corresponding desirable level that is associated with the then-constant scale signal (S) 545.
The bit extractor 549 takes the appropriate bits with sufficiently small quantization loss to produce the resultant signal (116,136) for the GNSS channel processing.
In
In an alternate embodiment, the complex channel of
In
To support dual PN demodulation (e.g., of L1C signal or an alternate BOC signal), the channel comprises two code demodulation paths. The first code enable unit 603 provides the first code enable signal 623 (e.g., first clock signal) to drive the first code phase accumulator 605. At each code enable signal 623 (e.g., first clock signal) the first code rate signal 627, generated by the code NCO 607, is accumulated by the first accumulator 605 (e.g., first code phase accumulator). If the first accumulator detects an overflow, the first accumulator 605 generates a first code advance signal 625. The first code advance signal 625 moves the PN coder 608 (e.g., first PN code generator) from the current state (e.g., current chip) to the next state (e.g., next chip) to generate the first PN sequence (e.g., which can be unique to each GPS satellite or channel on such GPS satellite).
The chip refers the unit of clock cycle in the GPS receiver that uses spread spectrum modulation, where the chipping rate of the code or PN sequence on each GNSS channel is generally known or published. If or when the received PN sequence is aligned with a replica of the PN code sequence, the code phase; hence, the associated pseudo-range between the receiver and corresponding satellite that transmitted the received PN sequence, can be precisely measured to about one chip resolution.
Because all the PN sequences on a given satellite are coherent (e.g., substantially coherent), and the frequency-dependent and frequency-independent channel perturbation is the identical for all the PN signal on the same frequency (e.g., same carrier frequency or same power spectral density function with encoded sideband components), the second code path is associated with the first code path to reflect such coherency. The first code advance signal 625 also drives the second code enable module 604 to generate the second code enable signal 638 (e.g., second clock signal), which further drives the second code phase accumulator 606.
A multiplier 690 scales the first code rate signal 627 (from the code NCO) by a multiple signal 639 to generate the second code rate signal 624. For some frequencies, the two PN code rates differ by a multiple of 10/20, for example, the combination of GPS L1-CA versus GPS L1-P as complex channel for the system
The first coder 608 generates the first set of the PN sequences with different delays, such as generating the signal 628-1 to the signal 628-m, where m is the mth signal in the first set of PN sequences. The second coder 609 generates the second set of the PN sequences with different delays, such as the generation of the signal 629-1 to the signal 629-m. The different delays may be structured as early, prompt and late delays for use (e.g., dual-use) in the system of
Each complex channel supports n accumulators, where n is a positive integer. The first PN selection unit 611-1 may comprise a multiplexor, for example. The first PN selection unit 611-1 selects or picks one PN sequence from either the first set of PN sequences 628-1 to 628-m, or from the second set of PN sequences 629-1 to 629-m based on the first selection signal 631-1 through the nth selection signal 631-n, which may be based on feedback from detection of a correlation peak and/or discriminator/envelope detector output. In one possible example, the first PN selection unit 611-1 outputs the first selected PN sequence 632-1, which the summer 691 (e.g., adder) combines or sums with the first overlay code (signal) 633-1, provided by the first code unit 612-1 (e.g., SCM[1]), to generate the first local code signal 637-1.
The first-integration-and-dump module 613-1 receives the first carrier baseband signal 636-1 (e.g., first carrier demodulated signal), the first local code signal 637-1, and the first W-code selection signal 614-1 (e.g., Wsel[1]). The first-integration-and-dump module 613-1 is configured to generate the first in-phase (e.g., I component) and quadrature-phase accumulation signal (Q-component) 634-1 of the P-code signal (e.g., P-code phase) based on the first local code signal 637-1 based on the first carrier baseband signal 636-1 (e.g., first carrier demodulated signal), the first local code signal 637-1, and the first W-code selection signal 614-1 (e.g., Wsel[1]). In one illustrative configuration, the first-integration-and-dump module 613-1 stores and makes available the first in-phase (e.g., I component) and quadrature-phase accumulation signal (Q-component) 634-1 of the P-code signal (e.g., P-code phase) for later application to a discriminator, envelope detector, or dot product evaluation of correlated signal amplitude or power spectral density to generate an error signal for code phase tracking and/or code NCO adjustment. In some configurations, data overlayed on the GNSS signal may result from encoding of the GNSS carrier signal with data, such as navigation-related data, encryption data, P(Y) codes and/or W-codes.
Similarly, the nth PN selection unit 611-n may comprise a multiplexor, for example. The nth PN section unit selects or picks one PN sequence from either the first set of PN sequences 628-1 to 628-m, or from the second set of PN sequences 629-1 to 629-m, based on the nth selection signal 631-n. In one possible example, the nth PN selection unit 611-n outputs the nth selected PN sequence 632-n, which the summer 692 (e.g., adder) combines or sums with the nth overlay code (signal) 633-n provided by the nth second code unit 612-n (e.g., SCM(n)) to generate the nth local code signal 637-n.
The nth-integration-and-dump module 613-n receives the nth carrier baseband signal 636-n (e.g., nth carrier demodulated signal), the nth local code signal 637-n, and the nth W-code selection signal 614-n (e.g., Wsel[n]). The nth-integration-and-dump module 613-n is configured to generate the first in-phase (e.g., I component) and quadrature-phase accumulation signal (Q-component) 634-n of the P-code signal (e.g., P-code phase) based on the nth local code signal 637-n based on the nth carrier baseband signal 636-n (e.g., nth carrier demodulated signal), the nth local code signal 637-n, and the nth W-code selection signal 614-n (e.g., Wsel[n]). In one illustrative configuration, the nth-integration-and-dump module 613-n stores and makes available the nth in-phase (e.g., I component) and quadrature-phase accumulation signal (Q-component) 634-n of the P-code signal (e.g., P-code phase) for later application to a discriminator, envelope detector, or dot product evaluation of correlated signal amplitude or power spectral density to generate an error signal for code phase tracking and/or code NCO adjustment. In some configurations, data overlayed on the GNSS signal may result from encoding of the GNSS carrier signal with data, such as navigation-related data, encryption data, P(Y) codes and/or W-codes.
The new L1C signal is designed for interoperability with the Galileo GNSS system and is also backward compatible with the current civil signal on L1 GPS. L1C signal will be transmitted at a higher power level and include some advanced design features associated with Galileo Binary Offset Carrier (BOC) modulation for enhanced performance.
The L2C signal is formed by multiplexing a first PN code (e.g., civil moderate code) and a second PN code (civil long code) and modulating the carrier with (e.g., by bipolar phase shift keying (BPSK)), where the civil moderate code is modulated by a navigation message that can be demodulated by civilian receivers in addition to the L1C/A navigation messages. The L2C civil long code signal can be used as a pilot, but is not representative of pilot PN sequence or pilot PRN sequence. The L2C signal supports dual-frequency civilian GPS receivers to correct the ionospheric group delay and to provide potentially faster signal acquisition. Further, the L2C signal provides enhanced reliability or greater operating range (e.g., from terrain attenuation, vegetation shading or multipath) because L2C broadcasts at a higher effective power than the legacy L1 C/A signal.
In
The first phase selection module 659-1 selects the GNSS received signal aligning with the L1P phase (plane) (as opposed to L1C or L1-C/A phase plane of the received GNSS signals for a given satellite) to generate the signal 648-1 for the first code demodulation processing; the nth phase selection 659-n selects the signal aligning with the L1P phase (plane) to generate the signal 648-n for the nth code demodulation processing. In one embodiment, each first phase selection module, 659-1 through 659n, comprises a multiplexer that is controlled by a selection signal from a channel baseband tracking loop module (e.g., code tracking loop, carrier tracking loop). Further, the channel baseband tracking loop comprises: (a) a carrier NCO that provides a local carrier frequency oscillator, or IF frequency oscillator, that is aligned with the L1P carrier phase (e.g., of its baseband I and Q components after pull-in or signal acquisition of a certain channel) of a received GNSS signal; (b) a correlator that accepts (samples of) such L1P-aligned local carrier frequency oscillator signal and (samples of) the received, evaluated GNSS signal (of the same satellite, such as a candidate L2C carrier signal) to be aligned to provide (candidate) correlations; (c) a carrier loop discriminator, alone to together with a carrier loop filter (e.g., Kalman filter), for evaluating (e.g., searching for) L1P phase (plane alignment) by estimating a peak correlation among such provided (candidate) correlations between L1P-aligned local carrier frequency oscillator signal and received, evaluated GNSS (e.g., L2C carrier signal); (d) the carrier loop discriminator configured to provide a carrier error signal for the carrier NCO to maintain such alignment (e.g., to select a phase-aligned L2C channel from the candidates).
The P code is available to the public; the P code and is the W-code are applied to an exclusive OR digital logic (e.g., logic gate or digital signal processing) to cryptographically generate a sequence or word that produces or yields the Y-code. For example, a GPS GNSS transmits the Y-code if the anti-spoofing module of the satellite is set to the “on” state by the U.S. government, or its agencies. The encrypted signal is generally referred to as the P(Y) code.
Because the L1P channel only processes or encodes one PN sequence at any given time consistent with the above coherency with the L1C signal, the first code enable unit 641 is configured to generate a first code enable signal 651 (e.g., first clock signal), or use a common code enable unit for the systems of
A code numerically controlled oscillator (code NCO) 649 is configured to generate a first code rate signal 649-1, such as a clocked, discrete-time representation of waveform (e.g., generally sinusoidal waveform or square-wave). The first code rate signal 649-1, which can comprise a derivative signal from the code rate signal 627 (of
The first P-coder 643 (e.g., first code generator) is configured to generate multiple PN sequences, such as 653-1 to 653-m, where m is the maximum number of possible PN sequences, which are offset in time or phase to provide different code phases (e.g., early, prompt, and late times), are fed into the multiple PN selection units, such as 644-1 to 664-n, where n is the maximum number (e.g., positive integer) of code phases or where n equals m. In one example, the first PN selection 644-1 comprise a multiplexor that selects the prompt or on-time PN sequence 653-3, which may represent 653-nOT (nOT=3 in this implementation or 653-3) or another PN sequence within 653-1 to 653-m, to generate the first code demodulation signal 654-1 through first code demodulation signal 654-n. For instance, the first code (demodulation) signal 654-1 through first code (demodulation) signal 654-n represent P-codes with corresponding I-components and Q-components, where prompt P-code could be selected based on correlation process (e.g., code correlator and code loop discriminator/envelope detector) as input to, or an integrated system of, the PN selection module (644-1).
Because of the encryption modulation on the L1P (e.g., L1P(Y)) signal, the first W-code accumulator 645-1 (e.g., WAcc[1]) combines the first code demodulation signal 654-1 and the first carrier (demodulated) signal 648-1 or first baseband signal from the respective phase selection module 659-1 to integrate over the encryption chip period to generate the first W-code accumulation signal 655-1. The unknown W-code may be applied to encrypt the P(Y) code at a known frequency (e.g., 500 Kilohertz (KHz)), which is less than the PN-code chip rate and which provides the above corresponding encryption chip period (e.g., inversely proportional to a fixed frequency, to the extent applicable) for integration. The W-code demodulation unit 646 squares or processes the signal 655-1 to wipe off (e.g., blindly wipe off) or remove the unknown W code modulation and to create the first W-code removed signal 656-1. Therefore, the L1P, signal I and Q (vector) components, can be demodulated or decoded from the modulated L1P signal with the carrier and P(Y) codes (e.g., for purposes code tracking and carrier tracking loops) even though the W-code is not known by the public or commercial end users.
In one embodiment, the first integration and dump module 647-1 further accumulates the W-code removed signal 656-1 to generate the multi-millisecond integration signal 657-1. For example, the first integration-and-dump module 647-1 (e.g., labeled I&D[1] in
In
For example, for the nth code demodulation signal 654-n, the nth W-code accumulator 645-n (e.g., labeled WAcc[n] in
In one embodiment, the nth integration and dump module 647-n (e.g., labeled I&D[n] in
The processing of the L2 signals in
In one embodiment, each first phase selection module, 689-1 through 689n, comprises a multiplexer that is controlled by a selection signal from a channel baseband tracking loop module (e.g., code tracking loop, carrier tracking loop). Further, the channel baseband tracking loop comprises: (a) a carrier NCO that provides a local carrier frequency oscillator, or IF frequency oscillator, that is aligned with the L1P carrier phase (e.g., of its baseband I and Q components after pull-in or signal acquisition of a certain channel) of a received GNSS signal; (b) a correlator that accepts (samples of) such L1P-aligned local carrier frequency oscillator signal and (samples of) the received, evaluated GNSS signal (of the same satellite, such as a candidate L2C carrier signal) to be aligned to provide (candidate) correlations; (c) a carrier loop discriminator, alone to together with a carrier loop filter (e.g., Kalman filter), for evaluating (e.g., searching for) L1P phase (plane alignment) by estimating a peak correlation among such provided (candidate) correlations between L1P-aligned local carrier frequency oscillator signal and received, evaluated GNSS (e.g., L2C carrier signal); (d) the carrier loop discriminator configured to provide a carrier error signal for the carrier NCO to maintain such alignment (e.g., to select a phase-aligned L2C channel from the candidates).
Because L2P channel only processes one PN sequence at a given time, the code enable signal 661 (e.g., clock signal) is synchronized with the code enable signal 623 (in
In
The second integration-and-dump unit 667-1 (e.g., labeled I&D[1] in
The nth PN selection module 664-n selects a PN sequence, from the set of signals (673-1 to 673-m, inclusive) to generate the corresponding nth code demodulation signal 674-n. For example, if n equals 2, the second PN selection module 664-2 selects a PN sequence, from the set of signals (673-2) to generate the corresponding nth code demodulation signal 674-2.
The nth W-code accumulator 665-n combines the nth code demodulation signal 674-n and the nth carrier demodulated signal 688-n to integrate over the encryption chip period to generate the nth W-code accumulation signal 675-n. For example, the second W-code accumulator 665-2 combines the second W-code demodulation signal 674-2 and the second carrier demodulated signal 688-2 to integrate over the encryption chip period to generate the nth W-code accumulation signal 675-2.
In one embodiment, the second W-code demodulation module 666 processes or combines the first W-code accumulation signal 655-1 (of
In alternate embodiments, the second W-code demodulation module 666 may process or combine W-code accumulated signals from the L1P channel or the L2P channel to create a corresponding W-code removed signal (676-1 through 676-n, inclusive) in accordance with various techniques that may be applied cumulatively or separately. Under a first technique, the second W-code demodulation module 666 may process or combine the nth W-code accumulation signal 655-n or another W-code accumulation signal (e.g., 655-1 through 655-n−1, inclusive) from L1P channel (of
The nth integration and dump unit 667-n further accumulates the W-code removed signal 676-n to generate the multi-millisecond integration signal 677-n. As previously noted, the decoded output (e.g., 677-1 through 677-n) of each integration-and-dump module (667-1 through 667-n) in
In one configuration, each NCO (771, 774, 776) may be used to form a phase-locked-loop controlled oscillator, which can receive phase alignment feedback from one or more of the following: correlators (723, 726), alone, or together with, post-correlation processing of signal strength evaluators; LOS estimation module 704; vector tracking module 705; or receiver data processing system. The system of
In the carrier tracking loop of the channel baseband tracking module 711, the carrier NCO 771 is coupled to one or more look-up tables 770, such as sine map and cosine map, to create a carrier signal (e.g., 714) for the carrier wipe-off function of the carrier demodulator 602. In one configuration, carrier signals generated by the look-up tables 770 are generally 90 degrees out of phase generally sinusoidal waveforms. Further, the carrier NCO 771 receives a stable clock signal (at a clock frequency or pulse train) from the clock tracking loop(s) 730. The clock signal may be corrected for first clock bias component of a mobile GNSS receiver clock to a respective satellite clock for each satellite, a second bias component of the respective satellite clock to a GNSS clock time for a GNSS constellation, and third bias component of the respective satellite block to a GNSS clock for a different GNSS constellation, among other things.
In the code tracking loop of the channel baseband tracking module 711, the code NCO 774 is coupled to a code generator 773, such as code generator for generating a P(Y) code or another PN code based on the inputted local code oscillator signal from the code NCO 774. In turn, the code generator 773 is coupled to shift register 772 to form prompt P code and late P code signals for input to the first correlators (723, 726) at the output of the shift register to complement the early P code signal, directly, from the code generator 773.
In the clock tracking loop module 730, the clock NCO 776 drives a waveform look-up table 775 to provide a pulse train, square-wave or other suitable clock signal at a target frequency and target phase based on the second correlators 726 for clock tracking and based on the first set of first correlators 723 for code loop tracking (e.g., code phase loop tracking) and carrier loop tracking (e.g., carrier phase look tracking), which collectively can be referred to as channel baseband tracking.
In some embodiments, an external sensor of
In certain embodiments, a skyward-facing or upward-facing imaging device, a monocular camera, a stereo vision camera, or a LIDAR (Light Detection and Raging) system, a radar (e.g., radio detection and ranging) system provides external data to the LOS estimation module 704. In turn, the LOS estimation module 704 is configured to estimate whether satellite signals from certain satellites (e.g., satellite identifiers based on an almanac of rising and setting times of the corresponding satellites and GNSS time) within the constellation are line-of-sight (LOS) propagation paths between satellite and GNSS receiver, or blocked or materially attenuated by an attenuation threshold with respect to local terrain, ground clutter, obstructions, vegetation (e.g., tree canopy), buildings or other structures.
For one or more epochs or sampling intervals, the LOS data estimation module 704 may do one or more of the following: (a) estimate or determine LOS data, motion-corrected LOS data, or Doppler-corrected LOS data; (b) determine a list of available GNSS satellite channels to be scanned, polled, surveyed or processed for (e.g., aggregate) carrier phase and code phase tracking of the GNSS receiver, such as measurement of pseudorange between the polled or scanned GNSS satellites and the GNSS receiver by eliminating blocked or materially attenuated signals of GNSS satellites that would otherwise be in view or reception range based on observations of a skyward-facing or upward-facing imaging device, a monocular camera, a stereo vision camera, radar system, a LIDAR system; (c) temporarily mask/ignore reception of, or de-weight carrier phase measurements and/or code phase measurements of GNSS satellites that are associated with obstructions that block or severely attenuate LOS GNSS signals traversing sky regions (e.g., semi-spherical zones or arc-bounded regions, such as near the horizon of the Earth or at low elevation angle, such as 20 degrees or less, with respect to the surface of the Earth) estimated by: (1) (previously) received or stored satellite (almanac/ephemeris/navigation) data indicative of excluded/blocked/attenuated satellites at or below a threshold low elevation angle, and/or (e.g., alone, or together with) (2) observed by a skyward-facing or upward-facing imaging device, a monocular camera, a stereo vision camera, radar system, or a LIDAR system that is co-located with the mobile or rover GNSS receiver; (d) temporarily exclude unavailable satellites (e.g., from navigation determinations, position estimates and raw carrier phase measurements and code phase measurements associated with occluded or blocked satellites signals or materially attenuated satellite signals) or de-weight unavailable satellites from navigation solution estimation, such as estimation of position, attitude or motion estimation solutions of the GNSS receiver, the rover/mobile GNSS receiver or reference GNSS receiver, and/or (e) temporarily exclude unavailable satellites from channel baseband processing or de-weight position data or LOS data until the LOS estimation module 704 (or its external sensor) determines that the corresponding satellite with a certain satellite identifier returns to reliable reception of GNSS satellite signal (e.g., of sufficient signal quality, signal strength, or geometric dilution of precision (GDOP), or low bit error rate of encoded data) or LOS for some minimum time period (e.g., hysteresis, transition delay window or transition dwell time).
In another embodiment, the external sensor, which is coupled to or integrated into the LOS estimation module 704, may comprise an inertial measurement unit (IMU) for providing aiding-motion data to the LOS estimation module 704. For example, the aiding-motion data may comprise any of the following: position versus time data, velocity data, and acceleration data associated with a corresponding mobile or rover GNSS receiver. The LOS estimation module 704 can estimate LOS velocity-aiding data or motion-aiding data associated with a corresponding GNSS satellite and respective GNSS rover or mobile receiver. For example, the LOS estimation module 704 may estimate Doppler shift in a pseudorange measurement based on the aiding-motion data from the IMU, or based on the applicable LOS velocity-aiding data to the signal propagation path of the received GNSS signal between a rover GNSS receiver and a corresponding GNSS satellite. Because the satellites orbit and are in constant motion, the Doppler shift may be detectable even if the reference, mobile or rover GNSS receiver is stationary, and the motion of the GNSS receiver, such as via the IMU, or one or more accelerometers or gyroscopes, can be used to estimate Doppler shift from movement of the mobile or rover GNSS receiver relative to any given GNSS satellite within view or reception range.
In an alternate embodiment, the LOS estimation module 704 receives aiding motion data from an IMU of the rover or mobile GNSS receiver; the LOS estimation module 704 estimates a compensating adjustment or time offset to the clock local oscillator (e.g., clock NCO) (e.g., of the clock tracking loop(s) 730), which is applied to the code local oscillator (e.g., code NCO) and/or the carrier local oscillator (e.g., carrier NCO) (e.g., of the channel baseband tracking loop(s) 711) to adjust the generated local code signal or code replica based on a Doppler shift of the received GNSS signal at the rover or mobile GNSS receiver.
In one embodiment, the vector tracking module 705 combines the signal tracking of a set of multiple GNSS satellite channels (e.g. all received GNSS channels) consistent with an estimated position, estimated velocity and estimated time (e.g., epoch) that is provided by a navigation, control and interface module 122 (e.g., in
In one embodiment, for a code wipe-off configuration of a phased-lock loop (PLL) with code NCO 774 (e.g., and carrier NCO 771) to process encoded data or data modulation (e.g., repeating navigation-related data message for a given satellite, until the message is updated over time) on the received GNSS signal, a bank of first correlators 723 (e.g., first correlators) for baseband carrier tracking may provide accumulations or correlations (e.g., 634-1 to 634-n, inclusive, and 735-1 to 735-m, inclusive) that indicate similarity between samples of the inputted baseband signal and corresponding samples of one or more replica code signals (e.g., prompt code replica signals or early, prompt and late code replica signals of a previously decoded or known repeating navigation data message for a given satellite until the message is updated), such as a primary correlation associated with an in-phase component and a secondary correlation associated with a quadrature-phase component (e.g., demodulated IQ signals) of a set of one or more corresponding GNSS channels (e.g., all polled, scanned or serially surveyed, available GNSS channels with reliable signal quality or signal strength).
The above code loop tracking is based on (e.g., tracking, alignment of time synchronization of) a carrier numerically controlled oscillator (NCO) 771 of the channel baseband processing module 711, where the carrier is modulated with the code for a corresponding satellite, satellite channel, or set of channels. For example, the channel baseband tracking loop 711 may further comprise: (a) a carrier loop discriminator, alone or together with a carrier loop filter, that is configured to select the correlations in which the prompt correlations have the greatest normalized magnitude or greater normalized magnitude (e.g., than associated early and prompt correlations), or (b) the carrier loop discriminator and a code loop discriminator that is configured to select the correlations in which the prompt correlations have the greatest normalized magnitude or greater normalized magnitude (e.g., relatively greater than associated early and prompt correlations for a given sampling interval, clock cycle, or chip). If a particular embodiment of the channel baseband tracking loop 711 only includes the carrier tracking loop discriminator(s), the clock tracking loops 730 may comprise the accompanying code tracking loop discriminators or clock tracking loop discriminators. Here, it is understood that the first correlators 723 may include integrators to average the products of prompt correlations and/or early, prompt, and late correlations for storage in corresponding accumulators 736.
In certain embodiments, the channel baseband tracking loop 711 may comprise a carrier loop discriminator that is configured select correlations from the carrier loop discriminator to adjust the carrier NCO 771 in the channel baseband and the clock NCO 776 within the clock tracking loop 730. Further, in some configurations the carrier NCO 771 will be adjusted, replica signals, based on the selected correlations from the carrier loop discriminator and any of the following: (a) LOS estimation data from the LOS estimation module 704, (b) position versus time data, velocity data or acceleration data (e.g., for Doppler adjustments to estimated pseudoranges between the GNSS receiver and any given GNSS satellite within view/reception range) from an external sensor of the LOS estimation module 704, such as an inertial measurement unit (IMU) or multi-axis accelerometer or gyroscope associated with the GNSS receiver, (c) clock bias data and corresponding clock bias compensating data, and (d) carrier bias data and corresponding carrier bias compensating data. For each received GNSS channel, each GNSS channel with a set of received GNSS channels, or a representative collective GNSS channel (e.g., super GNSS channel), the (Doppler) adjustment of the carrier NCO 771 is configured to adjust the time synchronization and temporal tracking of the samples of locally generated carriers or samples of carrier replicas with respect to the received intermediate frequency signal, received near-baseband signal and/or received baseband signal to produce clean carrier wipe-off accumulations (e.g., IQ accumulations) at digital baseband signal with IQ components (636-1 through 636-n) that are substantially free of unwanted images in the digital frequency domain, artifacts or residual carrier components.
In one configuration, a bank of second correlators 726 (e.g., second correlators), for clock tracking may provide accumulations or correlations (725-1, 725-m) that indicate similarity between samples of the inputted baseband signal and corresponding samples of one or more replica clock signals (e.g., delayed or shifted clock replica signals), of a set of one or more corresponding GNSS channels (e.g., all polled, scanned or serially surveyed, available GNSS channels with reliable signal quality or signal strength) based on (e.g., tracking, alignment of time synchronization of) a clock phase numerically controlled oscillator 776 of the clock tracking loop module 730. The clock tracking loop module 730 can support accurate tracking the code phase, the frequency of the code phase, or both based on precise control and adjustment of the numerically controlled oscillator (NCO) of the clock tracking loop module 730, where the clock tracking loops 730 can provide incremental code-phase adjustments to the code NCO 774 related to same GNSS channel or set of channels for given satellite.
Similarly, a bank of first correlators 723 (e.g., first correlators), for code tracking and/or carrier tracking may provide accumulations or correlations (634-1, 634-m) accumulations or correlations (e.g., in accumulators 736) that indicate similarity between samples of the inputted baseband signal and corresponding samples of one or more replica code signals (e.g., prompt code replica signals or early, prompt and late code replica signals), such as a primary correlation associated with an in-phase component of carrier signal and a secondary correlation associated with a quadrature-phase component of carrier signal (e.g., IQ signals for certain GNSS signals can be demodulated to decode navigation data). The first correlators 723 support a set of one or more corresponding GNSS channels (e.g., all polled, scanned or serially surveyed, available GNSS channels with reliable signal quality or signal strength) based on (e.g., tracking, alignment of time synchronization of) a clock phase numerically controlled oscillator 776 of the clock tracking loop module 730.
In an alternate embodiment, a channel code loop discriminator is coupled to the accumulators 736, or integrated within the accumulators 736; the channel code loop discriminator is configured to select the correlations in which the samples of prompt correlations have the greatest normalized magnitude or greater normalized magnitude (e.g., relatively greater than the magnitude associated early and prompt correlations). Further, in certain embodiments, the first correlators 736 may include integrators to average the products of prompt correlations and/or early, prompt, and late correlations for storage in corresponding accumulators 736, which are coupled to the output of first correlators 723.
Further, in an alternate embodiment, the clock tracking loop 730 may comprise a clock loop filter coupled to the output to filter the control feedback signal provided to the code NCO 774 in the code tracking loop and/or the carrier NCO 771 in the carrier tracking loop. For example, the clock tracking loop 730 with a code loop discriminator is configured to use the selected correlations at the output of the second correlators 726, the accumulators 738, or (a clock loop discriminator associated with the second correlators 726 and the accumulator 738) to adjust the code NCO 774, the carrier NCO or both. Further, in some configurations the code NCO 774 will be adjusted based on the samples of selected correlations (e.g., with prompt correlations of maximum magnitude) from the code loop discriminator and any of the following: (a) carrier accumulations and correlations (634-1 through 634-m, 735-1 to 735-m, or both) from the accumulators 736 for the carrier tracking loop module and/or code tracking loop module to derive carrier aiding data, (b) code accumulations and correlations (725-1 to 735-m) from the accumulators 738 for the code tracking loop module to derive a clock tracking error signal, a code bias signal, or control feedback signal, and (c) feedback or error data (e.g., 727-1 to 727-m) outputted from the vector tracking module 705.
In one embodiment, the vector tracking module 705 may provide feedback or error data to each numerically controlled oscillators (771, 774) of respective channel baseband tracking modules 711 (e.g., one channel baseband tracking loop per received GNSS satellite) and each clock NCO 776 of the respective clock tracking loops 730 (e.g., one clock tracking loop per received GNSS satellite) to replicate or determine a local estimation of one or more of the following: each carrier phase of corresponding GNSS channel, or set of corresponding GNSS channels, or a collective representative (e.g., super) GNSS channel that represents a set of corresponding GNSS channels; each code phase of a corresponding encoded GNSS signal or a set of corresponding GNSS channels, or a collective representative (e.g., super) GNSS channel that represents a set of corresponding GNSS channels; and each code frequency of a corresponding encoded GNSS signal or set of corresponding GNSS channels, or a collective representative (e.g., super) GNSS channel that represents a set of corresponding GNSS channels.
In an alternate embodiment, such local replicas or local estimations of carrier phase and/or code phase may be applied to an extended Kalman filter or Kalman filter within navigation, control and interface module 122 and/or the vector tracking module 705 for reduction of errors, such as phase tracking error (e.g., carrier phase tracking error, code phase tracking error), GNSS receiver clock bias, and satellite clock bias.
In some embodiments, at the output of the correlators (723, 726), the decoded or demodulated baseband signal may contain navigation-related data, such as GNSS time, almanac, ephemeris data, satellite status, data, and orbital data, encoded on the corresponding GNSS signal, which can be applied to the navigation control and interface module 122, the vector tracking module 705, or both. For example, the L1C/A is encoded with navigation data that can be fully demodulated by civilian GNSS receivers.
In one configuration, the vector tracking module 705 provides estimated carrier frequency and phase and estimated code frequency and phase (727-1 to 727-m) for each received GNSS satellite signal, consistent with input from a navigation, control and interface module 122 that may comprise an extended Kalman filter or Kalman filter that processes one or more demodulated/decoded GNSS channels from correlators (623, 726) or accumulators (736, 738) associated with the output of the correlators (623, 726) to reduce error in the decoded data of the each GNSS output channel.
The channel baseband tracking module 711 and the clock tracking loop(s) 730 may be configured in accordance with various technical approaches (e.g., embodiments) that may be applied cumulatively or separately.
Under a first technical approach, in each channel baseband tracking module 711 or each second-stage of the carrier demodulator 719, a local carrier generator, such as carrier NCO 771, for each channel or set of channels, (e.g., or aggregate local carrier generator), alone or in conjunction with a waveform look-up table 770 (e.g., sine and/or cosine look-up table), provide an estimated replica or an estimated (e.g., aggregate) local carrier signal (712 and/or 714) for a set of one or more corresponding GNSS signals to first correlators 723 or the carrier demodulator 602, such as the second-stage carrier demodulator 719). First correlators 723 may comprise one or more mixers, shift registers (e.g., to generate delayed signals), and integration-and-dump modules. Further under the first approach, in each code tracking loop, a local code generator, such as an code NCO 774, for each channel or a set of GNSS signals (e.g., aggregate local code generator) provides an estimated local replica (713) (e.g., P code, or P(Y) code) or an estimated aggregate local code signal for a corresponding set of GNSS signals.
Under a second technical approach, in each channel baseband tracking module 711, a local carrier generator, such as carrier NCO 770, for each channel or set of channels (e.g., a super channel or aggregate local carrier generator) provides an estimated replica or an estimated (e.g., an aggregate or super) local carrier signal (714, 712) for a set of one or more corresponding GNSS signals; further, in the same channel baseband tracking module 711, a local code generator, such as code NCO 774, for each channel or set of channels (e.g., aggregate or super local carrier generator) provides an estimated replica or an estimated code signal (e.g., an aggregate or super local code and local carrier signals 713) for a set of one or more corresponding GNSS signals.
Under a third approach, in the clock tracking loop(s) 730, a local clock generator, such as clock NCO 776, for each channel or set of channels (e.g., aggregate local carrier generator) provides an estimated replica or an estimated (e.g., aggregate) local clock signal (706), or incremental adjustments thereto for each clock cycle or chip, for a set of one or more corresponding GNSS signals, where each local clock signal is associated with any of the following: a corresponding respective received satellite GNSS signal, a corresponding GNSS satellite and respective satellite clock bias, and a particular reference, rover or mobile GNSS receiver and respective receiver clock bias. The local clock generator (e.g., clock NCO 776) can provide the estimated local clock signal to the second correlators 726, which comprise mixers and shift registers, and integration-and dump-modules.
The set of one or more estimated local code signals (713) (e.g., P codes) and estimated local baseband signals (636-1 to 636-n) (e.g., demodulated carrier signals) with I components and Q components are inputted to a corresponding first correlators 723 (or correlator bank) for baseband tracking of the carrier; the correlated output signal (634-1 to 634-n, inclusive), or decoded information signal is provided to each channel baseband tracking module 711 and clock tracking loop module 730 for a corresponding channel or set of channels of any given GNSS satellite. the correlated output signal or demodulated output (e.g., IQ vector signals) is stored in accumulators 736, from which it can be provided to the clock tracking loop module 730. Each clock tracking loop (730) is configured to provide a clock output signal 706 based on outputs of the correlator banks (723 or 726) that can be stored in accumulators (736, 738).
In an alternate embodiment, the carrier and code replicas for the same set of GNSS signals or same set of GNSS channels may share a common correlator bank, such as first correlators 723, while each second correlator 726 within the bank of second correlators 726 is dedicated to clock loop tracking module 730 for a corresponding GNSS satellite.
In one embodiment, the modulated or encoded intermediate frequency signal or modulated, encoded digital baseband signal (e.g., near baseband frequency signal) (e.g., 115, 135) is inputted to the demodulator 602 of
In
The second-stage carrier demodulator 719 receives input signals 714 (e.g., P codes of one or more phases, such as early, prompt and late P-codes) and 721 to remove (e.g., completely) a carrier component of the carrier signal to provide the encoded baseband signal (636-1 through 636-n) without a carrier (e.g., and with reduced carrier phase tracking error and/or reduced code phase tracking error arising from potential variance in coherence, or with reduced demodulation phase noise in the process), such as a encoded baseband signal (636-1 to 636-n, inclusive) with a removed/stripped carrier to the correlator bank of first correlators (723) for channel baseband tracking.
Meanwhile, the frequency scaling module 715 adjusts, translates, or scales the clock signal(s) 706 (e.g., clock estimation signals or incremental changes to the clock signal) from respective output(s) of a set of clock tracking loops 730 at reference frequency to the clock rate data or clock-related data 716 at a corresponding channel signal frequency or corresponding set of channels (e.g., coherent or time/phase synchronized GNSS channels), consistent with the channel tracking outputs (712, 713, 714) of the channel baseband tracking module 711. In one embodiment, the clock-related data 716, alone or together with carrier frequency data, contributes to carrier-aiding data to the code tracking loop, such as a tracking error signal component 621.
In an alternate embodiment, the frequency scaling module 715 can be coupled between the channel baseband processing module 711 and the summer 750 to scale the carrier phase error component.
The summer 750 adds or processes the scaled clock rate 716, which is derived from one or more clock signals, and the carrier frequency data 712 (or aggregate carrier phase error component and/or code phase error component) associated with the LOS data (e.g., from the same GNSS satellite) from the coherently or synchronously related channel baseband tracking module 711, alone or together with the LOS estimation module 704, to generate the tracking error signal 621 (e.g., with carrier phase tracking error component, a code phase error tracking component and clock tracking error component), which is provided to a demodulator 602 or digital mixer to remove the tracking error component (e.g., carrier tracking error signal component) from the baseband signal, among other things.
In one embodiment, a correlator bank of first correlators 723 (e.g., first correlator bank) mixes a bank of carrier-removed signals (636-1 to 636-n) or rather encoded baseband signals (e.g., modulated with a code), which are derived from the output of the second carrier demodulator 719, with a set of local code replica signals 713 to generate a bank in-phase and quadrature-phase (IQ) accumulation signal 735-1 to 735-m, inclusive (also denoted as signal 634-1 to 634-n, inclusive), which can be stored in accumulators 736. For example, the first correlators 723 may detect the modulated components or encoded components (e.g., IQ components) of the baseband signal based on product (e.g., a time-averaged, smoothed or integrated) of the inputted, encoded baseband signal and a locally generated signal or local code replica (e.g., or a set phase delayed replicas) that has some of the same characteristics as the encoded baseband signal.
In one configuration, the vector tracking module 705, the LOS data estimation module 704, or both combine the bank of correlations (735-1 to 735-m) from the first channel to the mth channel (or 641-1 to 641-n, inclusive) to produce the estimated LOS data 727-1 (e.g., for each applicable satellite and corresponding GNSS receiver) for the first channel and up to 727-m for the mth channel. Further, the external sensor (e.g., imaging system) of the LOS data estimation module 704 also independently produces the estimated LOS data (728-1 to 728-m) for the first channel and up to the mth channel (e.g., for each applicable satellite and corresponding GNSS receiver). The baseband tracking module 711 of the first channel uses correlation signal 735-1 (or 641-1), while the second through mth channel use correlations (735-1 to 735-m) to produce the residual frequency 714. For each channel or set of channels, baseband tracking module 711 is configured to use the LOS data signal (727-1 to 727-m), (728-1 to 728-m) and the derived satellite LOS to produce the LOS-affiliated carrier frequency 712 (e.g., LOS frequency or carrier frequency associated with LOS data).
The LOS carrier removed signal or local code replica signal 713, which can be generated by an matched filter with an impulse response that is reversed in time, is also used by the second correlator(s) 726 to generate the in-phase and quadrature-phase correlation 725-1 for the first channel. The clock tracking loop 730 takes, the clock correlation 725-1 and the channel correlation 735-1 of the first channel, up to the clock correlation 725-m and channel correlation 735-m of the mth channel to produce the clock frequency 706 at the reference frequency.
In
To reduce the data storage size, one narrow bandwidth signal is selected for each GMSS constellation. The multiplexer (MUX) 802 selects one of the GNSS constellation signals at digital baseband (115, 135) or an intermediate frequency to acquire and outputs selected signal 803. Because selected signal 803 may comprise the intermediate frequency component, the mixer 804 or digital mixer is used to translate the signal 803 into the digital baseband signal 806 (e.g., mixed signal).
In one embodiment, the low pass filter (LPF) 807 reduces the bandwidth of the selected signal 803 by filtering digital baseband signal 806, where the decimation of the LPF 807 is viable or capable to reduce the rate of the filtered signal 819. For example, the LPF 807 comprises a digital baseband finite impulse response (FIR) filter that decimates by ND, which tends to reduce images or unwanted artifacts in signal 819 from the digital mixer 804 that can potentially appear in the resultant signal 806.
In one embodiment, the level translation unit 808 further reduces the number of bits of each sample of the code (e.g., P(Y) encoded baseband) because the coarse acquisition doesn't need multi-bit precision as tracking requires. By the LPF 807 and the level translation unit 808 applying two successive decimations, the size of storage memory 809 can be significantly reduced.
A cold start of the GNSS receiver means initialization or turning on the GPS receiver without stored (e.g., locally stored) prior data to assist or reduce the pull-in or ambiguity resolution of carrier phase or acquisition period to attain a precise position solution (e.g., based on precise carrier delta pseudorange and code pseudo range solutions). A warm start of the GNSS receiver means initialization or turning on the GPS receiver with stored (e.g., locally stored) prior data to assist or reduce the pull-in or ambiguity resolution of carrier phase or acquisition period. The method of the acquisition process of a GNSS signal
In step S858, in preparation for or prior to starting the acquisition process of each GNSS signal (e.g., at least one GNSS signal per GNSS constellation), the pseudo-noise (PN) sequences (e.g., pilot or training PN or pseudo-random noise (PRN) sequences) from the received GNSS channel, set of GNSS channels, or aggregate GNSS channel (e.g., super channel) representative of the set of GNSS channels, are acquired are recorded, stored in, or written into, data storage 811 through the interface register 818. For example, the data storage 811 may comprise registers of an electronic data processor 827 (e.g., of the control module or GNSS receiver), nonvolatile random access electronic memory, electronic memory, magnetic storage device, a disk drive, or an optical storage device.
Further, in accordance with step S858, each GNSS channel, set of channels or aggregate channel (e.g., at least one GNSS signal per GNSS constellation) is: (a) down-converted to an intermediate frequency, baseband frequency (e.g., down-conversion with carrier wipe-off) or near-baseband frequency signal that comprises: (a) an in-phase component (I component), or (b) a quadrature component (Q-component), or both I and Q components with the encoded PN sequence (e.g., pilot or training PN sequence); (b) a correlation and code-wipe-off process is applied to the I component, the Q component, or both with the encoded PN sequence (e.g., to identify leading edges, trailing edges, pulses, or pulse transitions of the PN sequence for proper initial/preliminary alignment of carrier tracking loop and/or code tracking loop); (c) such PN sequence, which can comprise an I component, a Q component or both may represent a pilot PN sequence, or a training PN sequency for the signal acquisition process.
For example, the quadrature-component, Q-component (e.g., which can be post-correlation Q-component and prior to discrimination or envelope detection) of pseudo-noise (PN) sequences (e.g. training or pilot PN sequence) that are acquired are recorded, stored in, or written into, data storage 811 through the interface register 818. If the Q component of the received signal is used, the signal power of the Q component may predominately include noise energy when the carrier tracking is phase locked in the carrier phase locked loop; may contain greater energy in the signal components when the carrier tracking is not locked (e.g., during acquisition and searching modes).
In step S859, a sub-band or band selection module 802 (e.g., multiplexer), GNSS receiver or data processor 827 selects the GNSS signal 803 (e.g., channel, set of channels, or aggregate channel representative of a set) for acquisition or decoding from a bank of candidates, associated with the received GNSS signal(s) or (IF) digital signals (113, 133), consistent with the recorded or stored pilot PN sequences (e.g., associated with digital baseband frequency), where the GNSS signal may be associated with a GNSS channel, a set of GNSS channels, or an aggregate channel representative of the set of GNSS channels. Accordingly, consistent with the block diagram of
In step S860, based on the characteristics of the selected GNSS signal 803, the control module 813, GNSS receiver, or electronic data processor 827 is configured for the acquisition parameters, such as the coherent integration period (e.g., associated with control signal, such as enable coherent integration signal 822) and data/overlay code pattern or code specifications (e.g. associated with the bit pattern selection signal 821). The coherent integration period may vary between: (a) a lower limit of a tracking integration rate after acquisition, pull-in, or convergence of the carrier phase ambiguity solution of GNSS channel, where the tracking integration rate may depend upon the digital IF sampling rate or P-code encoding rate, and (b) an upper limit during an acquisition or search for proper data integration boundaries for corresponding encoded baseband signals. The data/overlay code pattern or code specifications (e.g., PN code specifications) may include the modulation frequency of the encoded information, such as the code chipping period or chipping rate of the P-code, reference power spectrum of the PN code plus data at baseband frequency, and autocorrelation period and autocorrelation interval of the P-code. Similarly, the code specifications for C/A (coarse-acquisition-code), publicly available specifications for encryption W-code; and autocorrelation period and autocorrelation interval for C/A code may supplement the PN code data for acquisition purposes.
In step S861, to manage the signal acquisition process, the control module 813, GNSS receiver or electronic data processor 827 generates or controls the state of the GNSS receiver to control one or more of the following control signals or commands: the bit pattern selection signal 821, coherent integration signal 822, Discrete Fourier Transform (DFT)/fast Fourier Transform (FFT) selection signal 823, LD_PNI signal 814, carrier frequency offset signal 824, and buffer memory control signal 826. The coherent integration signal 822 may support any of the following: enable integration, disable integration and integration period, lower limit of integration period and upper limit of integration period, acquisition mode and tracking mode. The DFT/FFT selection signal supports selection of a particular form of Fourier transform to simplify calculations of the dot product or evaluation of the signal power or energy of various I and Q components of the encoded baseband PN code. To provide accurate analysis, Fourier transforms may be based on an assumption that the amplifier and filter components are generally linear in a phase response versus frequency for any passband of an amplifier, filter, or both, for example.
The carrier frequency offset may support: enable, disable, and change/adjustments in the carrier frequency, intermediate frequency, frequency scaling, and associated phase of the carrier frequency or intermediate frequency. LD_PNI signal 814 (e.g., load PN sequence signal) may represent the enable, disable, load next PN sequence (or next I component and/or Q component, without or without integration or averaging, in the register, stack, memory queue or data storage. The LD-PNI or load next PN sequence also include code frequency, code phase, or code phase adjustment of output of local oscillator, such as code NCO that pertains to a corresponding or next PN sequence.
In step S862, the control module 813 is configured to provide the carrier frequency offset signal 824: (a) that selects the carrier frequency or intermediate frequency (IF) from the frequency offset lookup table(s) 810 (or look-up table 301); (b) where the look-up table(s) (810 or 301) drive or is driven by the carrier numerically controlled oscillator (NCO) (e.g., 805 or another NCO) which generates the local carrier signal 825 or local IF signal; and (c) where the mixer 804 combines the digital intermediate frequency signal or selected digital baseband signal 803 and local carrier 825 to generate the (encoded) digital baseband signal 806. For example, for certain GPS signals, the baseband signal 806 may be encoded with P(Y) codes, encryption W-codes, or both.
In step S863 in one embodiment, the LPF 807, which comprises a decimation module, reduces the rate of the baseband signal 806. The resultant filtered signal 819 is further quantized by the level translation unit 808 to reduce the quantization level, which is then stored in the buffer data storage device 809. For example, the buffer data storage device 809 may comprise electronic buffer memory for storing filtered baseband signal samples that are sampled at rate of a suitable sampling interval (e.g., millisecond(s) duration samples or chip-rate sampling duration of corresponding baseband modulated by P-code, precision code, encoded navigation-related data, encoded message or pseudo-noise code or pseudo random noise code (PRN)). A chip is a unit of clock cycle in a spread spectrum modulated GNSS received signal, such as GPS received signal that is encoded or modulated with a PN signal or PRN signal associated with navigation-related information, where the bandwidth of the modulated GNSS received signal is generally proportional to the chipping rate.
In step S864, the control module 813, GNSS receiver or data processor 827 is configured to generate buffer memory control signal 826 to align the first sample of the PN sequence (e.g., pilot or training PN sequence) in the buffer data storage 809 with the GNSS receiver (millisecond or chip-rate) clock edge or symbol transition of the clock signal or local replica of the phase code signal, such as the P(Y) encoded signal. In one embodiment, the LOS estimation module 704, alone or together with the external data source (e.g., IMU) may provide a Doppler shift, motion aiding data, or LOS data to adjust the clock frequency and/or phase of a clock oscillator (e.g., clock NCO) in one or more clock tracking loops 730, which in turn, in some configurations, may drive or reduce tracking error a carrier NCO and a corresponding code NCO, of the respective carrier tracking loop, code tracking loop, or channel tracking loop for one or more GNSS channels.
In step S865, the control module 813, GNSS receiver or data processor 827 is configured to generate (e.g., via one or more shift registers or digital delay units) the PN local signal 818 (e.g., the P(Y) encoded signal) with one or more phase offsets (e.g., early, prompt and late phase offsets or one or more chips proportional to the chipping rate of the modulation on the encoded GNSS signal channel) in data storage 811 and to transfer the scheduled PN sequence 812 into the data processing module 815 an electronic data processor 827 that communicates to associated data storage device via a data bus, where the electronic data processor 827 is configured with software instructions to provide one or more of the following: correlators, mixers, accumulators, and integrators.
In step S866, the data processing module 815, GNSS receiver or data processor 827 combines, integrates and/or correlates the PN code sequence (e.g., pilot or training PN code sequence in the buffer data storage) with the time-aligned GNSS receiver generated replica or local PN code sequences (with one or more corresponding phase offsets in the data processing module 815): (a) to generate multiple sub-millisecond integrations (e.g., from correlations or accumulations) to search for bit or word transitions in the PN sequence (e.g., pilot or training PN sequence for a corresponding GNSS channel or set of GNSS channels); (b) to evaluate the signal energy of the integrations, such as FFT/DFT signal products represent dot product power (e.g., substantially coherent dot product power) of various in-phase (I) components, and quadrature components (Q) with different time/phase offsets; and (c) based on signal energy analysis, to generate a carrier offset signal 824 or tracking error for aligning the carrier frequency/carrier phase, or change in carrier phase, with respect to the code phase of the local PN sequence or PN code replica in accordance with the carrier tracking loop, code tracking loop, and/or channel tracking loop and any accompanying tracking error signal.
For example, in accordance with control data from the control module 813, the data processing module 815, GNSS receiver, or data processor 827 can apply a Fourier transform, such as fast Fourier Transform (FFT) and Discrete Fourier Transform (DFT) to the integrations (e.g., sub-millisecond integrations or integrations over an encoding or encryption chip period). In particular, the control module 813 is configured to generate the FFT/DFT selection signal 823 to conduct the spectrum analysis (e.g., frequency versus magnitude response in the frequency domain) on the integrations (e.g., sub-millisecond integrations).
Step S866 may be accomplished by applying various techniques, cumulatively or separately. Under a first technique for executing step S866, the processing module 815, GNSS receiver or data processor 827 is configured to linearly combine, integrate, or otherwise manipulate or process selected FFT/DFT signal products, such as in-phase components and quadrature phase components, of the local replica PN code sequence and the received PN code sequence encoded on the baseband signal, in the frequency domain (e.g., based on the data bit pattern selection signal 825). For example, for an integration time (e.g., one or more chips, epochs or successive sampling intervals at a higher acquisition integration rate, rather than a lower steady-state tracking rate for carrier, code and/or clock loops), the FFT/DFT signal products may represent dot product power (e.g., substantially coherent dot product power) of various in-phase (I) components, and quadrature components (Q) with different time offsets (e.g., early, prompt and late time offsets) from processing the output of one or more correlators of carrier loop discriminators and/or code loop discriminators (e.g., delay lock loop discriminators).
Under a second technique for executing step S866, for an integration time (e.g., one or more chips, epochs or successive sampling intervals at a higher acquisition integration rate, rather than a lower steady-state tracking rate for carrier, code and/or clock loops), the control module 813 can provide a data bit pattern selection signal 821 or a control signal for the data processing module 815 to make and store (and optionally rank) various combinations (e.g., products or dot products) of DFTs or FFTs with greatest signal power (e.g., maximum power density or maximum signal magnitude over a frequency range of interest) at the current code shift, where combination comprises a product or integrated product of a sample or time-shifted sample (e.g., early, prompt or late sample) of local replica PN code sequence and sample of received PN code sequence encoded on the baseband signal. Further, to the extent that the local PN code sequence is temporally coherent, synchronized with or aligned to track the received PN code sequence on the baseband, the correlations between the local PN code sequence and the received PN sequence correspond to combinations of DFTs or FFTs with greatest signal power (e.g., maximum power density or maximum signal magnitude over a frequency range of interest) at the current code shift, such as prompt correlation, as opposed to an early (e.g., earlier) correlation or a late (e.g., later) correlation with a lower power density, or lower signal magnitude over a frequency range of interest.
Under a third technique, the control module 813 instructs the processing module 815 to remove or wipe-off codes, such as P(Y) codes, encryption, and/or encoded W-codes of the baseband signal (e.g., L1P or L2P signal with removed P(Y) codes) to detect, decode, or promote detection or decoding navigation-related data, along with bit error rate (BER), symbol error rate (SER) and other digital signal metrics to verify/confirm proper acquisition of one or more GNSS channels, or to continue searching for signal acquisition in accordance with the method of
Under a fifth technique, a GNSS channel, set of GNSS channels, or aggregate GNSS channel is not found or acquired, or has insufficient signal strength, signal-to-noise ratio, signal quality or GDOP based on insufficient received energy (e.g., for reception, decoding or demodulation) associated with the I and Q components and dot products thereof of the Fourier transform analysis.
In step S867, after making and evaluating (e.g., and optionally ranking) the combinations (e.g., products, averaged products or time-integrated products) of DFTs and FFT, the GNSS receiver, processing module 815, or data processor 827 increments or shifts by one sample (e.g., offset time or delay unit) the PN sequence (e.g., training or pilot sequence) in the process module 815, and the process or step S866 is repeated (e.g., by correlation with a prompt, early or late local code replica of the PN sequence) until a sufficient number (or all) the code phases (e.g., in the PN or PRN local sequence memory 811) have been evaluated or tried (e.g., exhaustively, completely or at least once). Then, the PN load signal 814 transfers the next PN sequence 812 into the process module 815. Process S866 to S867 repeats until all the PN sequence (e.g., pilot or training PN sequences) in the memory buffer 811 has been evaluated and iterated for a particular GNSS signal, such as a respective representative GNSS satellite signal per each constellation.
In step S868, if a sufficient threshold number of stored samples of pilot PN sequences are responsive to pull-in, initialize or establish preliminary, substantially coherent tracking of carrier phase and code phase in step S866 on a GNSS channel (e.g., that conforms to a reference substantially coherent dot product power of Fourier transforms, reference spectral energy density, minimum signal quality or minimum signal-to-noise ratio), the code search for a PN or PRN on a GNSS channel, set of GNSS channels, or aggregate GNSS channel representative of the set of the GNSS channels, is completed. Accordingly, once the data processing module 815, the GNSS receiver or the data processor 827 processes a sufficient number of stored samples of PN sequences (e.g., pilot or training PN sequences) for a corresponding GNSS channel, set of GNSS channels, or aggregate GNSS channel representative of the set of the GNSS channels, the next channel or set of channels is evaluated in according with step S869.
In step S869, if a sufficient threshold number of PN sequences (e.g., pilot PN sequences) is analyzed, evaluated or reviewed at a carrier frequency, GNSS channel, set of channels, or aggregate channel for a corresponding satellite in a constellation, the offset signal 824 can be changed or controlled for the next GNSS channel, or set of GNSS channels, or aggregate GNSS channel with a certain GNSS carrier frequency for a different GNSS satellite signal on the same GNSS constellation or a satellite of a different satellite GNSS constellation and repeats the steps S862 though S868.
In step S870, a GNSS receiver, a channel selector (e.g., sub-band selection module or multiplexer) 802, or an electronic data processor 827 selects a received GNSS signal (115, 135 in
In step 872, a control module 813 or electronic data processor 827 of the GNSS receiver is configured for signal acquisition and correlation, wherein the signal acquisition parameters comprise any of the following: a coherent integration period and data/overlay code pattern or code specifications based on the selected GNSS signal.
In step S874, a control module 813, a channel baseband tracking module 711 or carrier tracking loop module is configured to provide a carrier frequency offset signal that selects the carrier frequency or intermediate frequency (IF) from the frequency offset lookup table(s) that supports a carrier numerically controlled oscillator (NCO) 805 (e.g., carrier NCO) to generate the local carrier signal or local IF signal (e.g., with associated I and Q components as vectors), or change to the local carrier signal the local IF signal for carrier phase loop tracking; and mixing the digital intermediate frequency signal and local carrier to generate or translate a frequency of the (encoded) digital baseband signal.
In step S876, a mixer 804 or electronic data processor 827 mixes the digital intermediate frequency (IF) 803 and local IF signal 825 to generate or translated a frequency of the encoded digital baseband signal. Alternately, a mixer 804 or electronic data processor 827 mixes the received signal (e.g., 803′) and the local carrier frequency signal (e.g., 825′) to generate or translate a frequency of the (encoded) digital baseband signal.
In step S878, a low-pass filter (807, 203, 213 or 114, 134) or bandpass filter is configured to filter, or filter and decimate the digital baseband signal to reduce or eliminate aliasing or analog-to-digital conversion artifacts.
In step S880, a control module 813 or electronic processor generates a buffer memory control signal to align the first sample of the pilot PN sequence in the buffer data storage device 809 with a clock edge or symbol transition of the clock signal or local replica of the phase code signal.
In step S882, the PN coder or code generator, alone or together with one or more shift registers or digital delay units, generates the PN local signal with one or more phase offsets of one or more chips proportional to the chipping rate of the modulation on the encoded GNSS signal channel in data storage device 811 and transfers each scheduled PN sequence of the PN local signal into a data processing module 815.
In step S884, the data processing module 815 or one or more correlators (e.g., 723 in
In step S886, data processing module 815 or electronic data processor 827 generates multiple sub-millisecond integrations from correlations or accumulations in the data storage device, or the data processing module 815, or both to search for bit or word transitions in the pilot PN sequence for a corresponding GNSS channel or set of GNSS channels.
In step S888, data processing module 815 or electronic data processor 827 evaluates the signal energy of the integrations of the correlations or accumulations based on Fourier transform signal products that represent dot product power or substantially coherent dot product power of various in-phase (I) components (e.g., I vectors) and corresponding quadrature components (Q) (e.g., Q vectors) with different time/phase offsets (e.g., early, prompt (e.g., on-time) and late offsets, expressed in chips or phase).
In step S890, the data processing module 815 or electronic data processor 827 generates a tracking error (e.g., with a carrier phase error component, a code phase error component, and/or clock bias component) based on signal energy analysis where tracking error is configured to align the carrier frequency/carrier phase, or change in carrier phase, with respect to the code phase of the local PN sequence or PN code replica in accordance with the carrier tracking loop, the code tracking loop, and/or collective channel tracking loop. For example, the data processing module 815 or electronic data processor 827 generates a tracking error (e.g., with a carrier phase error component, a code phase error component, and/or clock bias component) based on signal energy analysis where tracking error is configured to align the carrier frequency/carrier phase, or change in carrier phase, with respect to the code phase of the local PN sequence or PN code replica in accordance with the carrier tracking loop, the code tracking loop, and/or collective channel tracking loop in conjunction with (or corrected for) clock bias (e.g., satellite clock bias, rover GNSS receiver bias, inter-constellation GNSS bias) incremental updates to the clock signal inputted to the carrier NCO and code NCO.
In certain configurations, the code tracking loop or code NCO can be updated (e.g., indirectly for clock bias) in the form of incremental code-phase updates that include a clock bias component or clock bias update at an update rate that is based on, at least partially, a real-time, ionosphere-propagation-corrected, pseudo-range estimate, which is susceptible to dithering reception of certain satellite signals, between a rover or mobile GNSS receiver and a corresponding satellite. Similarly, the carrier tracking loop or carrier NCO can be updated (e.g., indirectly for clock bias) in the form of incremental carrier-phase updates that include a clock bias component or clock bias update at an update rate that is based on, at least partially, a real-time, ionosphere-propagation-corrected, pseudo-range estimate, which is susceptible to dithering reception of certain satellite signals, between a rover or mobile GNSS receiver and a corresponding satellite.
In step S900, a selection module 802 selects a received GNSS signal (e.g., 102) as a channel, set of channels, or aggregate channel representative of the set, for acquisition to acquire the received GNSS (102) signal that is susceptible to Doppler frequency shifts or propagation-related frequency shifts.
In step S902, a frequency offset module 810 provides a carrier frequency offset signal to generate one or more candidates of the local carrier frequency signal or local intermediate frequency (IF) signal based on evaluation of signal energy associated with correlations. For example, each of the candidates of the local carrier frequency signal generally has relative phase offsets with respect to others of the candidates. Further, the frequency offset module 810 is configured to provide a carrier frequency offset signal to generate one or more candidates of the local carrier frequency signal or local intermediate frequency (IF) signal based on feedback. The feedback comprises evaluation of signal energy associated with correlations and a frequency hypothesis (for the true or synchronized local carrier frequency signal relative to the carrier of the GNSS received signal), where each of the candidates having relative phase offsets with respect to others of the candidates. For example, the frequency hypothesis is configured to depend upon factors such as the GNSS system, the GNSS satellite, the GNSS received signal and encoding/modulation, relative movement of the rover GNSS receiver and satellite transmitting the received GNSS signal.
In step S904, a mixer 804 mixes the received GNSS signal and the local carrier frequency signal or local carrier intermediate frequency signal to provide a baseband signal in which a carrier of the received GNSS signal is removed (e.g., wiped-off).
In step S906, a low-pass filter 807 filters (e.g., or filters and decimates) the received samples of digital baseband signal that is encoded by a received pseudo random noise code (PN) sequence. For example, the low-pass filter 807 is configured to low-pass filter 807 and decimate of the digital baseband signal to reduce a sampling rate of the baseband signal. Further, an electronic data processor or translation module 808 (e.g., L level translation module) quantizes the filtered signal to reduce possible quantization levels for (reduced or efficient) storage in a buffer data storage device or electronic memory of the GNSS receiver.
In step S908, a control module 813 generates a buffer memory control signal to attempt to align temporally one or more received samples of the received PN sequence, or a portion thereof, in a buffer data storage device with a clock edge or symbol transition of the clock signal of a set of local samples of corresponding PN local sequence, or portion thereof, of a local signal or PN replica signal. In practice, the signals can be aligned by adjusting, by a chip or fractional chip, the clock edge or symbol transition of the block signal or a set of local samples of the corresponding PN local sequence, or portion thereof, of a local signal or PN replica signal. If the PN local sequence has a data size that exceeds the respective size of registers within an electronic data processor or data processing module 815, a portion of the PN local sequence, which is commensurate with the available respective size of the registers, is processed with respect to a corresponding portion of the received PN sequence; such that the alignment of the PN local sequence to the received PN sequence may take multiple iterations of evaluations of portions of PN sequences in accordance with correlations and integrations of the signal acquisition process of
In step S910, one or more shift registers or digital delay units configured to generate the set of the local samples of the local PN sequence or local PN signal with respective one or more phase offsets. For example, the set of local samples of the local PN sequence, or portion thereof, with one or more phase offsets is stored in memory 8111.
In step S912, one or more memory devices (809, 811) transfer the received samples of each scheduled received PN sequence, or portion thereof, and the corresponding local samples of the local PN sequence, or portion thereof, into a data processing module 815 or a set of correlators.
In step S914, a set of correlators or the data processing module 815 correlates the received samples of the received PN code sequence, or portion thereof, in a buffer data storage (e.g., 809) with the respective set of local samples of the local PN sequence, or portion thereof, in memory (e.g., 811) to pursue identification of a temporally aligned (e.g., aligned integration and dump phase for code-loop discrimination/detection in signal acquisition mode), local PN code sample associated with the corresponding selected, received GNSS signal.
In step S916, a plurality of integrators or the data processing module 815 generates integrations of the correlations at millisecond or sub-millisecond intervals to identify clock edge or symbol transitions in the received PN code sequence, or portions thereof. Step S916 may be carried out in accordance with various techniques, which may be applied separately and cumulatively. In accordance with a first technique, the integrators or the data processing module 815 generates multiple sub-millisecond integrations from the correlations or accumulations in the data storage device, the data processing module 815, or both to search for bit or word transitions, as identifiable symbol transitions, in the received PN sequence, or portion thereof, for a corresponding GNSS channel or set of GNSS channels based on a data hypothesis (e.g., of the PN code or encoded information on the respective GNSS channel) consistent with publicly available specifications of the received PN code sequence. For example, the data hypothesis may comprise any of the following (e.g., publicly available technical information): the modulation type (e.g., Bipolar Phase Shift Keying (BPSK) or Binary Offset Carrier (BOC)), carrier frequency, modulation frequency or rate, PN code rate (e.g., Mega-chips/second), length or size of the PN code, navigation message data modulation rate (e.g., bits per second), code overly, code pattern, pilot PN code sequence, portion of pilot PN code sequence, selection of bit pattern from bit pattern library, symbol duration, duration of period or full cycle of the PN code before it is repeated, among other things.
In accordance with a second technique, the integrators or data processing module 815 generates multiple sub-millisecond integrations from the correlations or accumulations in the data storage device, the data processing module 815, or both to search for bit or word transitions, as identifiable symbol transitions, in the received PN sequence, or portion thereof, for a corresponding GNSS channel or set of GNSS channels based on one or more of the following associated with the selected received GNSS signal: a recorded pilot PN sequence, a stored pilot PN sequence, a coherent integration period, a data/overlay code pattern, code specifications, or specifications related to the data hypothesis, frequency hypothesis or both.
In accordance with a third technique, the GNSS received signal comprises a L1 C/A (coarse-acquisition signal) or a L2C signal that is modulated with a navigation data message.
In accordance with a fourth technique, a control module 813, an electronic data processor, or the data processing module 815 of the GNSS receiver is configured for signal acquisition and correlation, wherein the signal acquisition parameters comprise any of the following: a coherent integration period and data/overlay code pattern or code specifications based on the selected GNSS signal, or specifications related to the data hypothesis, frequency hypothesis or both.
In accordance with a fifth technique, a control module 813, an electronic data processor or the data processing module 815 is configured to manage one or more of the following control signals or commands of the control module 813 or the electronic data processor to facilitate step S916, alone or together with step S918: a bit pattern selection signal, a coherent integration signal, a Fourier transform (e.g., Fourier transform type or parameter) selection signal, a load stored PN sequence signal, a carrier frequency offset signal, and buffer memory control signal, or specifications related to the data hypothesis, frequency hypothesis or both.
In step S918, a data processing module 815 or an evaluator (e.g., discriminator) evaluates the signal energy of the integrations of the correlations between received samples and local samples for each sampling interval or epoch, where the candidate local carrier frequency or candidate local IF corresponding to the correlations with the greatest signal energy or magnitude with identifiable symbol transitions is generally indicative of acquisition of or the identification of the proper temporally aligned carrier frequency offset of the GNSS signal among the generated candidates to compensate for the Doppler frequency shifts or propagation-related frequency shifts.
Step S918 may be carried out in accordance with various examples, which may be applied separately and cumulatively. In accordance with a first example, the integrators, evaluators, discriminators or data processing module 815 evaluates the signal energy of the integrations of the correlations or accumulations based on Fourier transform signal products that represent dot product power or substantially coherent dot product power of various in-phase (I) components, and quadrature components (Q) with different time/phase offsets of the received GNSS signal.
In accordance with a second example, during or after identifying the proper temporally aligned carrier frequency offset, the integrators, evaluators, discriminators or data processing module 815 determine the Fourier transform signal products, which comprise fast Fourier transform signal products, discrete Fourier transform signal products, representative of dot product power of various in-phase (I) components, and quadrature components (Q) with different time offsets comprising early, prompt and late time offsets arising from processing the output of one or more correlators of the acquisition engine, and carrier loop discriminators and/or code loop discriminators.
In accordance with a third example, after the proper temporally aligned carrier frequency offset is identified, the evaluators, channel tracking module, or data processing module 815 tracks an error based on signal energy analysis, of the integrations of the correlations. For example, the carrier offset signal or tracking error is configured to align the carrier frequency/carrier phase, or change in carrier phase, with respect to the code phase of the local PN sequence or PN code replica in accordance with the carrier tracking loop, code tracking loop, and/or channel tracking loop and any accompanying tracking error signal.
In accordance with a fourth example, during a signal acquisition mode for an integration time at a higher acquisition integration rate than a lower steady-state tracking rate for carrier, code and/or clock loops, the evaluators, channel tracking module or data tracking module providing a data bit pattern selection signal or a control signal for the data processing module 815 to make and store various ranked combinations of products or dot products of Fourier transforms with greatest signal power, where the signal power comprises maximum power density or maximum signal magnitude over a frequency range of interest at the current code shift associated with the correlating.
In accordance with a fifth example, after making and evaluating and ranking combinations of products, or time-integrated products of discrete Fourier transforms and/or fast Fourier transforms, the GNSS receiver, shift registers, delay lines or data processing module 815 shifting by one sample the PN pilot sequence in one or more registers of the processing module 815 or data storage device. Further, in accordance with the fifth example, in conjunction with one or more of the steps S912, S914 and S916, repeating a correlation process by correlating the shifted pilot PN sequence with the local code replica of the PN sequence until a sufficient number of the code phases in the data storage device or registers have been evaluated to acquire or pull-in the carrier frequency of the selected, received GNSS signal that is compensated for Doppler frequency shift or propagation-related frequency shift.
In accordance with a sixth example, if the maximum integration (e.g., maximum magnitude of integrations), at the selected carrier frequency, is sufficiently above a threshold (e.g., signal or channel threshold), the evaluators, discriminators or data processing module 815 determines that the code search (e.g., and associated code shift) is complete for a PN on a GNSS channel, set of GNSS channels, or aggregate GNSS channel representative of the set of the GNSS channels; further, the GNSS receiver or its data processing module 815 uses the selected carrier frequency and the code shift to pull-in, initialize or establish preliminary, substantially coherent tracking of carrier phase and code phase on a GNSS channel that conforms to the threshold (e.g., signal or channel threshold), such as a reference spectral energy density, minimum signal quality or minimum signal-to-noise ratio.
In step S920, the integrators or the data processing module 815 generate multiple sub-millisecond integrations from the correlations or accumulations in the data storage device (e.g., 809. 811), the data processing module 815, or both to search for bit or word transitions, as identifiable symbol transitions, in the received PN sequence, or portion thereof, for a corresponding GNSS channel or set of GNSS channels based on a data hypothesis consistent with publicly available specifications of the received PN code sequence. In one embodiment, a data hypothesis means a data pattern, data overlay, data structure, chip rate, data modulation rate, modulation parameters, and/or data technical specifications that provide timing or other information about the encoded data or PN code(s) on the corresponding GNSS channel or set of GNSS channels, or bit or word transitions, or symbol transitions in the received PN sequence. For example, the data hypothesis may comprise any of the following (e.g., publicly available technical information): the modulation type (e.g., Bipolar Phase Shift Keying (BPSK) or Binary Offset Carrier (BOC)), carrier frequency, modulation frequency or rate, PN code rate (e.g., Mega-chips/second), length or size of the PN code, navigation message data modulation rate (e.g., bits per second), code overly, code pattern, pilot PN code sequence, portion of pilot PN code sequence, selection of bit pattern from bit pattern library, symbol duration, duration of period or full cycle of the PN code before it is repeated, among other things.
The chip rate may represent the unit of a clock cycle in the receiver or spread spectrum receiver that uses a PN sequence.
In step S922, an evaluator, discriminator, or the data processing module 815 evaluates the signal energy of the integrations of the correlations or accumulations. For example, an evaluator, discriminator, or the data processing module 815 evaluates the signal energy of the integrations of the correlations or accumulations (e.g., integration and dump processing) based on Fourier transform signal products that represent dot product power or substantially coherent dot product power of various in-phase (I) components, and quadrature components (Q) with different time/phase offsets of the received GNSS signal.
In step S924 for integration time at a higher acquisition integration rate than a lower steady-state tracking rate for carrier, code and/or clock loops, an integrator or data processing module 815 provides a data bit pattern selection signal or a control signal for the data processing module 815 to make and store various ranked combinations of received signal energy evaluations or products (e.g., dot products of Fourier transforms) with greatest signal power, where the signal power comprises maximum power density or maximum signal magnitude over a frequency range of interest at the current code shift associated with the correlating.
In step S926, after making and evaluating and ranking combinations of products, or time-integrated products of discrete Fourier transforms and/or fast Fourier transforms, the GNSS receiver or data processing module 815, shifts by one sample the PN pilot sequence (e.g., data overlay or data code pattern) in one or more registers of the processing module 815 or data storage device (e.g., 809, 811).
In step S928, the data processing module 815 repeats a correlation process by correlating the shifted pilot PN sequence (e.g., shifted data overlay or data code pattern) by a chip or fraction of a chip with the local code replica of the PN sequence until a sufficient number of the code phases in the data storage device (e.g., 809, 811) or registers have been evaluated to acquire or pull-in the carrier frequency of the selected, received GNSS signal that is compensated for Doppler frequency shift or propagation-related frequency shift.
In step S930, if the maximum integration (e.g., maximum magnitude of the integration), at the selected carrier frequency, is sufficiently above a threshold, the selected carrier frequency and the code shift is used to pull-in, initialize or establish preliminary, substantially coherent tracking of carrier phase and code phase on a GNSS channel that conforms to a reference spectral energy density, minimum signal quality or minimum signal-to-noise ratio, the code search for a PN on a GNSS channel, set of GNSS channels, or aggregate GNSS channel representative of the set of the GNSS channels, is completed.
In step S932, a frequency offset module 810, such as frequency offset look-up table, provides a carrier frequency offset signal to generate one or more candidates of the local carrier frequency signal or local intermediate frequency (IF) signal based on evaluation of signal energy associated with correlations and a frequency hypothesis, each of the candidates having relative phase offsets with respect to others of the candidates, wherein frequency hypothesis is configured to depend upon factors such as the GNSS system, the GNSS satellite, the GNSS received signal and encoding/modulation, relative movement of the rover GNSS receiver and satellite transmitting the received GNSS signal. For example, a frequency offset module 810 or frequency offset look-up table is associated with or drives a numerically controlled oscillator 805 (NCO) to generate a precise local oscillator 805 signal (e.g., cosine signal, sine signal, or both to product I and Q components at the received digital baseband sign) for input to the mixer 804.
The Global Navigation Satellite System (GNSS) receiver architecture is well-suited to mitigate multiple forms of electromagnetic (EM) interference, such as interfering microwave or radio frequency signals. To address and mitigate electromagnetic interference, the GNSS receiver is susceptible to minimal EM interference in realization of a practical digital filter configuration of reasonable logic complexity because of the inherent technical limits of finite quantization. For example, the dual-stage digital down-conversion sufficiently suppresses the harmonic-related signal components arising from the inherent technical limits of finite quantization, where the dual-stage, down-conversion design can facilitate reduced, reasonable logic complexity. In particular, the dual-stage down-conversion design, which comprises an analog primary downconverter and a digital secondary downconverter alone, or in conjunction with, analog low pass filter and digital low pass filter facilitate aliasing suppression to suppress or attenuate harmonics arising from quantization that would otherwise occur in conjunction with a direct, single stage downconverter that a digital-to-analog converter that converts an analog baseband signal to a digital baseband signal.
The GNSS receiver is well-suited to reduce, ameliorate, or mitigate integrated or combined wideband interference (WBI) and narrowband interference (NBI), such that the interference mitigation system supports improvement the GNSS receiver performance when exposed to various EM interference sources. For example, the GNSS receiver features a digital automatic gain control (DAGC) to adjust the magnitude of the samples, outputted from a filter or other interference-rejection module, to a desired level in a reduced precision system.
The GNSS receiver supports a flexible, configurable channel structure to support various carrier phase or pseudorange measurements based on the individual or combined signal tracking. For example, the GNSS receiver facilitates digital data processing of an aggregate channel, such as a super-channel bundle, to synchronize the shared parameters between different GNSS signals, such as the civil signal and the encrypted signal, while dealing with the technical differences in the GNSS signal components separately. For potentially enhanced efficiency, the signal acquisition engine can be based on batch processing to speed up the initial signal detection with a cold-receiver start and signal recovery/reacquisition associated with a warm-receiver start, where the warm-receiver can use previously stored data on a position solution, resolved integer ambiguities, or bias values to promote reacquiring time of arrival data for one or more GNSS satellite signals, or for reconvergence and establishing carrier phase lock on one or more GNSS satellite signals . . . .
In certain configurations, the LOS estimation module interfaces with an external sensor, such as an IMU, to enable the signal tacking aiding across the tracked satellites.
This document (including the drawings) claims priority and the benefit of the filing date based on U.S. provisional application No. 63/363,277, filed Apr. 20, 2022, and U.S. provisional application No. 63/268,221, filed Feb. 18, 2022, under 35 U.S.C. § 119 (e), where the provisional applications are hereby incorporated by reference herein.
Number | Date | Country | |
---|---|---|---|
63268221 | Feb 2022 | US | |
63363277 | Apr 2022 | US |