WO2013184933A1 - Signal receiver with group delay and amplitude distortion compensation - Google Patents
Signal receiver with group delay and amplitude distortion compensation Download PDFInfo
- Publication number
- WO2013184933A1 WO2013184933A1 PCT/US2013/044541 US2013044541W WO2013184933A1 WO 2013184933 A1 WO2013184933 A1 WO 2013184933A1 US 2013044541 W US2013044541 W US 2013044541W WO 2013184933 A1 WO2013184933 A1 WO 2013184933A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- signals
- received signals
- signal receiver
- frequency
- Prior art date
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
- H04B17/21—Monitoring; Testing of receivers for calibration; for correcting measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
- G01S19/235—Calibration of receiver components
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2662—Symbol synchronisation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
- H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
- H04L27/38—Demodulator circuits; Receiver circuits
- H04L27/3809—Amplitude regulation arrangements
Definitions
- the disclosed embodiments relate generally to the field of signal processing at signal receivers, and in particular to a system and method for compensating for group delay distortions and, optionally, amplitude distortions at a signal receiver.
- Signal receivers often include apparatus and circuitry (such as, analog signal processing circuitry) to condition received signals to have desired amplitude, phase and/or frequency characteristics.
- the apparatus and circuitry at the signal receivers may also be designed to eliminate or reduce the effects of interfering noise and other distortions introduced into the signals along the signal propagation path, prior to being received at the signal receivers.
- the signal receivers may introduce additional signal distortions or undesired characteristics to the received signals.
- phase distortion (sometimes herein called group-delay distortion).
- Phase distortions introduced by the components of the signal receiver result from a non-linear phase response of the respective components of the signal receiver. Phase distortions are observed as non-linear phase vs. frequency characteristics, or a non-uniform group delay over the frequency pass band of interest. Phase distortions limit the ability of the signal receiver to combine information from signals having different frequencies, thereby limiting the accuracy of measurements derived from wideband signals obtained from combining multiple signals or multi- frequency signals.
- Amplitude distortions introduced by the components of the signal receiver result from a non-symmetric magnitude response of the respective components of the signal receiver. Amplitude distortions are observed as non- symmetric magnitude vs. frequency characteristics over the frequency pass band of interest.
- Some embodiments provide a signal receiver that includes an antenna interface for receiving signals from an antenna, analog signal processing circuitry coupled to the antenna interface for processing the received signals to produce filtered signals, sampling circuitry to sample the filtered signals so as to produce digitized received signals, a compensation mechanism to receive the digitized received signals and compensate for nonuniform group delay and amplitude distortion introduced by the analog signal processing circuitry to produce compensated digitized received signals, and a digital processor to process the compensated digitized received signals so as to produce a result.
- Some embodiments provide a method of signal compensation performed at a signal receiver.
- the method includes configuring a compensation mechanism at the signal receiver in accordance with processed digital signals received by the compensation mechanism in response to injection of a calibration signal at the signal receiver by a calibration signal generator, the configuring including configuring the compensation mechanism to compensate for non-uniform group delay and amplitude distortion introduced by analog signal processing circuitry of the signal receiver.
- the method further includes receiving signals at the signal receiver, processing the received signals at the signal receiver to produce digitized received signals, and compensating the digitized received signals using the compensation mechanism to produce compensated signals.
- Some embodiments provide a method of signal compensation performed at a signal receiver.
- the method includes receiving signals from an antenna; processing the received signals with analog signal processing circuitry to produce filtered signals; sampling the filtered signals with sampling circuitry to produce digitized received signals;
- the compensation mechanism is configured by injecting a calibration signal at the signal receiver; transforming digital signals, corresponding to the injected calibration signal, received by the compensation mechanism from time-domain to frequency-domain to produce a complex- valued frequency-domain representation of the processed digital signals; and extracting magnitude and phase values from the complex- valued frequency-domain representation of the processed digital signals.
- Configuring the compensation mechanism further includes computing phase residuals by subtracting the extracted phase values from target frequency-dependent phase information, computing magnitude response residuals by determining ratios of the extracted magnitude values with target magnitude response information, generating a set of complex values from the computed phase residuals and magnitude response residuals, and transforming the set of complex values from frequency-domain to time-domain to produce complex- valued time- domain filter coefficients corresponding to a finite impulse response filter.
- Figure 1 illustrates an example application (a GNSS satellite navigation system) for a signal receiver with group delay compensation, in accordance with some embodiments.
- Figure 2 is a block diagram illustrating a signal receiver with group delay compensation, in accordance with some embodiments.
- Figures 3A-3F include block diagrams illustrating components of a signal receiver with group delay compensation, in accordance with some embodiments.
- Figure 4A includes a flow diagram illustrating configuring of a compensation mechanism of a signal receiver with group delay compensation, in accordance with some embodiments.
- Figures 4B-4C include flow diagrams illustrating configuring of a
- Figures 5A-5D include flow charts illustrating methods of signal
- first means "first,” “second,” etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another.
- a first contact could be termed a second contact, and, similarly, a second contact could be termed a first contact, without changing the meaning of the description, so long as all occurrences of the "first contact” are renamed consistently and all occurrences of the second contact are renamed consistently.
- the first contact and the second contact are both contacts, but they are not the same contact.
- Figure 1 includes a block diagram illustrating an example of a system (e.g., a
- GNSS Global Navigation Satellite System
- Signal Receiver 120 a signal receiver with group delay and amplitude compensation, according to some embodiments.
- GNSS 100 comprises one or more GNSS satellite(s) 110 and Signal Receiver
- GNSS satellite(s) 110 transmits signals (e.g., signals containing navigation information) to be received by Signal Receiver 120.
- the signals transmitted by GNSS satellite(s) 110 travel through various propagation media (e.g., layers of the atmosphere) prior to reaching Signal Receiver 120.
- the propagation media include undesired Noise and Distortion Source(s) 130 which introduce noise and/or distortions into the transmitted signals, causing the signals transmitted by GNSS satellite(s) 110 to be impacted detrimentally.
- Some detrimental effects of noise and distortions introduced by Noise and Distortion Source(s) 130 include degradation of signal to noise ratio (e.g., by the introduction of undesired noise power with frequency content outside— and often within— the frequency bands of interest), and undesirable and often non-uniform, non-deterministic changes to the amplitude, frequency and phase content of the transmitted signals.
- Signal Receiver 120 typically includes analog and digital circuitry to at least partially mitigate or counteract the detrimental effects of Noise and Distortion Source(s) 130 on the received signal.
- Signal Receiver 120 includes signal conditioning elements (e.g., filters and amplifiers) in the analog signal processing circuitry that selectively emphasize signals having frequencies of interest, and reject or attenuate signals that do not have frequencies within the bandwidth of interest.
- signal conditioning elements e.g., filters and amplifiers
- FIG. 2 is a block diagram illustrating a Signal Receiver 120 in accordance with some embodiments.
- Signal Receiver 120 receives signals from Antenna 202.
- Signal Receiver 120 includes Antenna Interface 204, Receiver(s) 206 and Compensation Mechanism(s) 212 (also herein called digital compensator(s)) for every frequency band of interest, and Digital Processor 214.
- Signal Receiver 120 also optionally includes Calibration Signal Generator 230, Clock 240, Housing 250, and Circuit Board 260.
- Antenna Interface 204 receives signals from an antenna (e.g., Antenna 202) and optionally from a calibration signal generator (e.g., Calibration Signal Generator 230).
- an antenna e.g., Antenna 202
- a calibration signal generator e.g., Calibration Signal Generator 230
- Receiver(s) 206 includes Analog Signal Processing Circuitry 208 and
- Analog Signal Processing Circuitry 208 is coupled to Antenna Interface 204 for processing the received signals to produce filtered signals.
- Analog Signal Processing Circuitry 208 includes various frequency, amplitude and phase conditioning components, such as, one or more analog filters and/or one or more gain (amplification) stages.
- Analog Signal Processing Circuitry 208 corresponds to or includes a low noise amplifier.
- Analog Signal Processing Circuitry 208 includes a (quadrature) demodulator to down-convert the received signals to produce orthogonal baseband signals (for example, orthogonal I and Q signal components).
- Sampling Circuitry 210 samples the filtered signals from Analog Signal Processing Circuitry 208 so as to produce digitized received signals.
- circuitry for producing the digitized received signals further includes quantization circuitry and digitization circuitry. It may be noted that in some implementations, Signal Receiver 120 includes separate Receivers 206 (and consequently, separate Analog Signal Processing Circuitry 208 and Sampling Circuitry 210) for two or more frequency bands of interest. Optionally, Signal Receiver 120 includes separate Receivers 206 for every frequency band of interest.
- Receivers 206-a (and consequently, Analog Signal Processing Circuitry 208-a and Sampling Circuitry 210-a) is configured to operate at the LI frequency band (e.g., 1575.42 ⁇ 16 MHz; or 1559 MHz - 1591 MHz)
- Receiver 206-b (and consequently, Analog Signal Processing Circuitry 208-b and Sampling Circuitry 210-b) is configured to operate at the L2 frequency band (e.g., 1227.6 ⁇ 16 MHz; or 1211 MHz - 1243 MHz)
- Receiver 206-c (and consequently, Analog Signal Processing Circuitry 208- c and Sampling Circuitry 210-c) is configured to operate at the L5 frequency band (e.g., 1176.45 ⁇ 16 MHz; or 1160 MHz - 1192 MHz).
- Compensation Mechanism(s) 212 receives the digitized received signals and compensates for non-uniform group delay and, in some implementations, amplitude distortion introduced by the analog signal processing circuitry 208 to produce compensated digitized received signals.
- the compensation mechanism e.g., Compensation Mechanism 212 compensates for signal distortions in the digitized received signals produced by the analog signal processing circuitry in one or more predefined ranges of frequencies.
- the signal distortions are phase distortions introduced by the non- linear phase response of Analog Signal Processing Circuitry 208.
- the non- linear phase response of Analog Signal Processing Circuitry 208 causes undesirable variations in group delay when processing signals that span a desired range of frequencies.
- the signal distortions include amplitude distortions introduced by a non- symmetric magnitude response of Analog Signal Processing Circuitry 208.
- the non- symmetric magnitude response of Analog Signal Processing Circuitry 208 causes undesirable effects on signals that span a desired range of frequencies.
- Compensation Mechanism(s) 212 is implemented in a digital signal processor, sometimes called a DSP. In some other implementations,
- Compensation Mechanism(s) 212 is implemented in circuitry, such as an FIR filter and one or more registers. Alternatively, or in addition, Compensation Mechanism(s) 212 is implemented in software executed by a general purpose processor. It may be noted that the term "non-uniform," as used herein, indicates having different delays for different frequencies in a predefined range of frequencies. In some implementations, Signal Receiver 120 is configured to include a separate Compensation Mechanism 212 for every frequency band of interest.
- Compensation Mechanism 212-a is configured to operate at the LI frequency band (e.g., a frequency band that includes 1559 MHz - 1591 MHz)
- Compensation Mechanism 212-b is configured to operate at the L2 frequency band (e.g., a frequency band that includes 1211 MHz - 1243 MHz)
- the L2 frequency band e.g., a frequency band that includes 1211 MHz - 1243 MHz
- Compensation Mechanism 212-c is configured to operate at the L5 frequency band (e.g., a frequency band that includes 1160 MHz - 1192 MHz).
- Digital Processor 214 processes compensated digitized received signals so as to produce a Result 220.
- the result e.g., Result 220
- the result includes a range to satellite, ranges to multiple satellites, navigation result(s), geographical location(s), and/or satellite time value(s).
- Digital Processor 214 is implemented using one or more microprocessors or other programmable processors.
- Digital Processor 214 is further described herein with reference to Figure 3F.
- Digital Processor 214 is implemented within or close to the baseband system; as such, Digital Processor 214 is configured to operate on baseband signals.
- Digital Processor 214 optionally includes GNSS Processor 216 and Microprocessor 218.
- GNSS Processor 216 typically includes circuitry, such as correlators, for analyzing signals received from GNSS satellites and thereby assisting Microprocessor 218 to perform navigation functions and optionally other functions.
- Digital Processor 214 includes and executes control instructions for controlling when Calibration Signal Generator 230, described below, injects a calibration signal at Antenna Interface 204 and for controlling configuration of Compensation Mechanism 212 (for example, through control signals provided by Microprocessor 218) in accordance with digitized received signals received in response to the injection of the calibration signal at Antenna Interface 204.
- Calibration Signal Generator 230 is coupled to
- Antenna Interface 204 to inject a calibration signal at Antenna Interface 204.
- Clock 240 provides synchronized clock timing signals to Calibration Signal Generator 230 and Sampling Circuitry 210.
- Antenna Interface 204 Analog Signal Processing
- Circuitry 208 Sampling Circuitry 210, Compensation Mechanism 212, Digital Processor 214 and Calibration Signal Generator 230 are all contained within Housing 250.
- Antenna Interface 204 Analog Signal Processing
- Circuitry 208, Sampling Circuitry 210, Compensation Mechanism 212, Digital Processor 214 and Calibration Signal Generator 230 are mounted on a single circuit board (e.g., Circuit Board 260).
- Antenna Interface 204 is not mounted on the circuit board on which the other components are mounted.
- Circuit Board 260 is contained within Housing 250.
- the signal receiver e.g., Signal Receiver 120
- the signal receiver is a satellite signal receiver.
- FIG. 3A is a block diagram illustrating an Antenna Interface 204 in accordance with some embodiments.
- Antenna Interface 204 includes one or more filters (e.g., Filter(s) 304-a and Filter(s) 304-b) to limit the frequencies of Received Signals 302 to frequencies of interest.
- Filter(s) 304-a and Filter(s) 304-b include filters with fixed or variable (e.g., tunable) properties.
- Antenna Interface 204 also includes one or more Amplifiers 306 for amplifying or strengthening signals of interest.
- Amplifiers 306 may include one or more amplifiers with fixed or variable (e.g., tunable) properties.
- FIG. 3B is a block diagram illustrating a Compensation Mechanism 212 and an optional Configuring Mechanism 320.
- Compensation Mechanism 212 receives digitized received signals (e.g., Digitized Received Signals 312) and compensates for non-uniform group delay and, in some implementations, amplitude distortion introduced by Antenna Interface 204 and Analog Signal Processing Circuitry 208 to produce digitized received signals (e.g., Digitized Received Signals 312) and compensates for non-uniform group delay and, in some implementations, amplitude distortion introduced by Antenna Interface 204 and Analog Signal Processing Circuitry 208 to produce
- digitized received signals e.g., Digitized Received Signals 312
- amplitude distortion introduced by Antenna Interface 204 and Analog Signal Processing Circuitry 208
- compensated digitized received signals e.g., Compensated Digitized Received Signals 316.
- Compensation Mechanism 212 is configured (for example, using Configuring Mechanism 320) to process Digitized Received Signals 312 using a finite impulse response filter or FIR filter (e.g., Complex FIR 314) using a set of complex filter coefficients (e.g., Complex FIR coefficients 334) to compensate for nonuniform group delay and, in some implementations, amplitude distortion introduced by at least the analog signal processing circuitry (e.g., Analog Signal Processing Circuitry 208).
- a finite impulse response filter or FIR filter e.g., Complex FIR 314
- a set of complex filter coefficients e.g., Complex FIR coefficients 3344
- the compensation mechanism e.g., Compensation
- Mechanism 212) is configured (e.g., by Configuring Mechanism 320) in accordance with signals received by the compensation mechanism in response to injection of the calibration signal at the antenna interface, for example by Calibration Signal Generator 230.
- Configuring Mechanism 320 performs a method of configuring the compensation mechanism as described with reference to Figures 5A-5D below.
- Coherent Sample Accumulator 322 (shown in more detail in Figure 3E), which is used to store and coherently accumulate (across multiple calibration cycles) samples of the digitized received signals received in response to injection of the calibration signal, and thus outputs coherently accumulated samples.
- coherent accumulation is performed by coherently combining samples of the digitized received signals received in response to injection of the calibration signal— across multiple calibration cycles— using methods such as summing, averaging, filtering (e.g., by the use of Kalman filters or Bayesian estimators) or any combination thereof.
- Coherent Sample Accumulator 322 may be implemented in hardware (e.g., accumulation registers, clock dividers, etc.) or in software executed by a digital signal processor or general purpose microprocessor.
- An FFT 324 (e.g., N-Point FFT 324-a, shown in Figure 3E) performs a complex Fourier Transform on the coherently accumulated (e.g., summed or averaged) samples to produce a set of complex values for a plurality of frequency components.
- the output of FFT 324 is, or includes, N complex values, each having a real component and an imaginary component.
- FFT 324 may be implemented in software, hardware (e.g., on an FFT chip), or on a digital signal processor.
- the output of the FFT 324 is converted to produce the measured phase response (which is called the Actual Phase Response 402 in Figure 4A) and, in some implementations, the measured magnitude response (which is called the Actual Magnitude Response 412 in Figure 4B and Actual Magnitude Response 422 in Figure 4C) of Antenna Interface 204 and Analog Signal Processing Circuitry 208.
- the measured phase response represents the measured phases of the filtered signal (e.g., the digitized received signal received in response to the injected calibration signal) at various frequencies present in the filtered signal.
- the measured magnitude response represents the measured magnitude of the filtered signal (e.g., the digitized received signal received in response to the injected calibration signal) at various frequencies present in the filtered signal.
- the measured phase response (which is called the Actual Phase Response 402 in Figure 4A) is obtained from a polynomial fit (e.g., polynomial smoothing) to the measured phases of the filtered signal (e.g., the digitized received signal received in response to the injected calibration signal) at various frequencies.
- the measured phase response is represented (e.g., approximated) by a polynomial having a numerical order; the numerical order being selected so as to eliminate undesired noise in the polynomial representation of the measured phases of the filtered signal, while preserving (e.g., approximating) the effects of group delay distortions on the measured phases of the filtered signal.
- the polynomial is then sampled at a plurality of frequencies (e.g., corresponding to frequency bins of N-point FFT 324-a) to produce smoothed measured phases of the filtered signal at the plurality of frequencies.
- smoothed measured phases are generated (e.g., computed) from the phase values (e.g., Actual Phase Response 402) of corresponding frequency components of the digitized received signals using interpolation methods (e.g., curve- fitting, polynomial interpolation, spline interpolation, Gaussian interpolation, regression-based methods, etc.).
- the measured magnitude response (which is called the Actual Magnitude Response 412 in Figure 4B and Actual Magnitude Response 422 in Figure 4C) is obtained from the measured magnitudes of the filtered signal at various frequencies by similar methods as those described above for obtaining the measured phase response.
- Target Response Curve 326 corresponds to a desired set of frequency-dependent phase values and/or frequency-dependent magnitude values stored previously or (optionally) generated from Actual Phase Response 402 and/or Actual
- the target phase values (e.g., Target Response Curve 326 or Target Response 404) are generated (e.g., computed) from a straight line fit (e.g., Straight Line 404-a) to the phase values (e.g., Actual Phase Response 402) of corresponding frequency components of the digitized received signals.
- the target phase values are generated (e.g., computed) from the phase values (e.g., Actual Phase Response 402) of corresponding frequency components of the digitized received signals using interpolation methods (e.g., curve- fitting, polynomial interpolation, spline interpolation, Gaussian interpolation, regression-based methods, etc.).
- the target phase values are fixed values determined in advance, independent of the measured phase response (e.g., Actual Phase Response 402).
- the target magnitude values are obtained from magnitude values (e.g., Actual Magnitude Response 412, shown in Figure 2C) of corresponding frequency components of the digitized received signals by similar methods as those described above for obtaining the target phase values.
- the target magnitude values are fixed values determined in advance, independent of the measured magnitude response (e.g., Actual Magnitude Response 412).
- residual phase differences e.g., Residuals 328 or Phase
- Residuals 406) are computed by subtracting the phase values (e.g., Actual Phase Response 402) of corresponding frequency components of the digitized received signals from the target phase values (e.g., Target Response Curve 326 or Target Response 404).
- phase values e.g., Actual Phase Response 402
- target phase values e.g., Target Response Curve 326 or Target Response 404.
- residual magnitude differences are computed by determining ratios of the extracted magnitude values (e.g., Actual Magnitude Response 412) of corresponding frequency components of the digitized received signals with the target magnitude response information (e.g., Target Response Curve 326 or Target Response 404).
- An inverse Fast Fourier Transform (e.g., iFFT 330) is performed on a representation of the residual phase differences (e.g., Residuals 328 or Phase Residuals 406) and, in some implementations, magnitude response residuals (e.g., Magnitude Residuals 416) to produce an initial set of complex values corresponding to FIR coefficients.
- iFFT 330 may be implemented using software executed by a digital signal processor or general purpose microprocessor, or hardware (e.g., on an FFT circuit). In some
- FFT 324 and iFFT 330 are implemented using the same FFT circuit, or the same program(s) executed by a digital signal processor or general purpose microprocessor.
- a windowing operation (Windowing 332) is performed on initial set of complex values obtained from iFFT 330, which conditions the results produced by iFFT 330 by reducing the magnitude of the edge components of the results of iFFT 330 by multiplying the results of iFFT 330 by a windowing function (e.g., a Tukey window, or even more specifically a Tukey 0.5 window).
- Windowing 332 is implemented using software executed by a digital signal processor or general purpose microprocessor, or hardware. "Edge components" of results of iFFT 330 are typically components not in a predefined central portion of the results.
- the center components include at least components 32 to 96
- the edge components include components 1 to 31 and 97 to 128 or a subset of components 1 to 31 and 97 to 128.
- any of a number of windowing functions are used to window the set of complex values obtained from iFFT 330 to produce the complex FIR coefficients, such as windowing functions described in G. Heinzel, A. Rudiger, R. Schilling, "Spectrum and spectral density estimation by the Discrete Fourier transform (DFT), including a comprehensive list of window functions and some new flat-top windows," Internal Report, Max-Planck-Institut fur Gravitationsphysik, Hannover, 2002; and FJ Harris, “On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform," Proceedings of IEEE, vol. 66, pp. 51-83, 1978.
- DFT Discrete Fourier transform
- Configuring Mechanism 320 The method performed by Configuring Mechanism 320 is explained further in relation to Figures 5A-5D.
- FIG. 3C is a block diagram illustrating a Finite Impulse Response Filter 314 as used in Compensation Mechanism 212.
- FIR Filter 352 e.g., Complex FIR 314, shown in Figure 3B
- Complex Filter Coefficients 354 e.g., Complex FIR coefficients 334, shown in Figure 3B.
- Complex Filter 354 e.g., Complex FIR coefficients 334, shown in Figure 3B.
- Coefficients 354 for FIR Filter 352 are stored in a Coefficients Register 356.
- Coefficients Register 356 is implemented internally within FIR Filter 352.
- Complex Filter Coefficients 354 are generated (e.g., by Configuring Mechanism 320) by processing calibration signals injected at antenna interface 204 during one or more calibration periods.
- Filter Coefficient Register 356 are obtained from Filter Coefficients Library 358.
- Filter Coefficients Library 358 stores multiple sets of filter coefficients as shown in Figure 3C.
- a multiplexer 359 is used to select a set of filter coefficients from the Filter Coefficients Library 358, for example, based on operating conditions (e.g., temperature) in or near the signal receiver 120.
- FIR Filter 352 filters digitized received signals using the selected set of set of filter coefficients from the Filter Coefficients Library 358.
- the multiple sets of filter coefficients stored in Filter Coefficients Library 358 are filter coefficients generated at a plurality of different operating conditions (e.g., a plurality of measured temperatures).
- the process of generating the multiple sets of filter coefficients in Filter Coefficients Library 358 is sometimes called factory calibration, as the coefficients may be generated "at the factory” during or shortly after manufacturing the signal receiver.
- filter coefficients generated through factory calibration are obtained using an external calibration signal (i.e., a calibration signal generated external to the receiver), and in such embodiments the signal receiver need not include Calibration Signal Generator 230.
- Figure 3D is a block diagram illustrating a Calibration Signal Generator 230.
- Calibration Signal Generator 230 is coupled to the Antenna Interface 204, to inject a calibration signal at Antenna Interface 204.
- the calibration signal produced by Calibration Signal is the calibration signal produced by Calibration Signal
- the Generator 230 includes a time domain signal corresponding to orthogonal frequency components in a predefined range of frequencies.
- the duration (or period) of the time domain signal is equal to, or an integral multiple of, the inverse of the frequency spacing between the orthogonal frequency components.
- the calibration signal includes a step carrier signal. In some embodiments, the calibration signal includes a swept carrier signal.
- the calibration signal includes evenly spaced frequency components, such as frequency components from XI MHz to X2 MHz in Y MHz increments (e.g., 1 MHz increments), where XI and X2 MHz correspond to the frequency span or bandwidth of the received signal (received at Antenna Interface 204), sometimes called the passband of interest, and Y MHz defines the (typically uniform) spacing (i.e., frequency interval) between consecutive frequency components.
- some implementations provide a duration (or period) of the calibration signal equal to, or an integral multiple of, the inverse of the frequency interval (i.e., spacing - e.g., Y MHz) between the consecutive frequency components.
- the calibration signal is a multi-tone (sometimes called a multi- carrier) Calibration Signal 368 that is an OFDM
- the calibration signal corresponds to a multi-tone signal that spans the passband of interest.
- the multi-tone signal is generated, for example, by a method employed by OFDM modulators.
- a control signal e.g., Enable, Disable, Timing reference
- Control Signal 360 is used to control the timing and/or duration of injection of the calibration signal.
- Digital Processor 214 controls, at least in part through Control Signal 360, Calibration Signal Generator 230.
- Digital Processor 214 determines the timing and/or duration of injection of the calibration signal, including the interval between injection instances and the duty cycle associated with the calibration signal.
- the timing and/or duration of injection of the calibration signal is equal to, or an integral multiple of, the inverse of the frequency spacing between consecutive (and orthogonal) frequency components that constitute Calibration Signal 368.
- Control Signal 360 controls Clock Generator 361. In some embodiments,
- Clock Generator 361 provides a timing reference to Quadrature Modulator 366 and to a Store 362 that stores a quadrature time series for generating orthogonal tones, also called multi-tone time-domain digital signals.
- the multi-tone time-domain digital signals are read from Store 362 and converted into analog signals by one or more Digital-to- Analog converter(s) (e.g., DAC 364-a and DAC 364-b), which provide the resulting analog signals to Quadrature Modulator 366.
- Quadrature Modulator 366 modulates the analog signals onto a carrier signal to produce Calibration Signal 368.
- Store 362 stores a time-domain representation of a target set of multi-tone signals, which are converted into an analog signal by one or more DAC's 364 and then modulated onto a carrier signal to produce the Calibration Signal 368.
- the calibration signal generator e.g., Calibration Signal
- the calibration signal injects the calibration signal at Antenna Interface 204 at a duty cycle less than a predefined fraction (e.g., at a duty cycle less than one percent).
- the calibration signal is injected at a low duty cycle (e.g., a duty cycle of less than 10 per cent, 1 one percent, 0.1 percent or 0.02 percent).
- calibration signals are added for an insertion period (sometimes also called a calibration period) once per duty period (e.g., a period of time between 20 milliseconds and 2 seconds).
- the insertion period of the calibration signal has a duration of no less than one "symbol period" of the Calibration Signal 368 (e.g., the period of time required to maintain orthogonality of all the frequency components in the calibration signal). Typically, the insertion period of the calibration signal has a duration of two or more symbol periods.
- a very low duty cycle calibration signal minimally impacts the GNSS signal carrier to noise power density ratio (CNo).
- the calibration signal generator injects the calibration signal at the antenna interface as a pseudorandom or random signal, to reduce any detrimental impact on the various GNSS signal structures.
- the calibration signal has an amplitude at least a predefined amount (e.g., 40 dB) greater than an amplitude of the received signals from the antenna.
- the average power of the calibration signal is less than the broadband noise floor for the signal receiver.
- the amplitude and duration of the calibration signal is determined in accordance with a criterion that the average power of the calibration signal is maintained lower than a value which would impact gain settings of a gain control circuit (e.g., an Automatic Gain Control circuit) of the signal receiver.
- the gain control circuit corresponds to the AGC circuit of the analog signal processing circuitry of the signal receiver.
- the calibration signal is injected after fixed or variable intervals of time, such as 10 milliseconds, 10 seconds, 10 minutes, 20 minutes, etc. Maintaining intervals of time between the injection instances of the calibration signal facilitates keeping the average calibration signal power lower than a value (corresponding to the broadband noise floor of the signal receiver) that would cause the gain control circuit (e.g., the AGC circuit) to change its gain settings.
- the gain control circuit e.g., the AGC circuit
- the AGC circuit may be switched off (or prevented from changing its state) during injection of the calibration signal.
- the calibration signal generator e.g., Calibration Signal
- Calibration Signal Generator 230 is controlled, at least in part, by the digital processor (e.g., Digital Processor 214) of the signal receiver 120.
- the digital processor would determine the timing and/or duration of injection and potentially the interval between calibration signal injection instances.
- Calibration Signal Generator 230 under control of Digital Processor 214, injects different calibration signals for each receiver 206, with the calibration signal for each receiver 206 having a distinct set of multi-tone time-domain digital signals in Store 362, a respective duty cycle, and distinct timing from the calibration signal injected for the other receivers 206.
- FIG. 3E is a block diagram illustrating Accumulated FFT 370, according to some embodiments.
- Accumulated FFT 370 includes two or more registers (for example, N-Stage Shift Register 322-a and N-Stage Accumulator 322 -b), each register to store and coherently accumulate (for example, using methods described elsewhere in this document) samples from a respective orthogonal baseband signal of two or more orthogonal baseband signals generated from the received signals (e.g., Digitized Received Signals 312).
- a coherent accumulation (e.g., sum or average) of the samples is calculated by performing a coherent (e.g., time-synchronized) summation of samples over M calibration cycles and optionally then by dividing the resulting values by the number (e.g., M) of calibration cycles (e.g., implemented using a shift register) to produce a coherent sum (and optionally, a coherent average) of samples across the M calibration cycles.
- a coherent accumulation e.g., an average or summation
- the coherent accumulation across calibration cycles allows for better calibration accuracy.
- Accumulated FFT 370 includes Coherent Sample Accumulator
- Accumulated FFT 370 receives Digitized RF Samples 372 (for example, obtained from Digitized Received Signals 312), Sample Clock 374 and Symbol Time Enable (duty cycle) 376. Accumulated FFT 370 processes Digitized RF Samples 372 to produce N Complex Frequency Domain Samples over M-Symbols 378, as described below.
- Coherent Sample Accumulator 322 is used to store and coherently accumulate
- Digitized RF samples 372 e.g., Digitized Received Signals 312 received by the
- Digitized RF samples 372 are shifted and stored using an N-Stage Shift Register 322-a and coherently accumulated (e.g., combined by summation or filtering approaches), and optionally, averaged using N-Stage Accumulator 322-b.
- N-Stage Shift Register 322-a coherently accumulated (e.g., combined by summation or filtering approaches), and optionally, averaged using N-Stage Accumulator 322-b.
- FFT 324 (e.g., N-point FFT (over M symbols) 324-a) performs a complex
- a timing reference for the signal samples (e.g., Digitized RF Samples 372) is provided, in some implementations, using Sample Clock 374 (for example, obtained from or synchronized with Clock 240).
- Sample Clock 374 for example, obtained from or synchronized with Clock 240.
- Another timing reference, Symbol Time Enable (duty cycle) 376 enables operation of the Coherent Sample Accumulator 322 and FFT 324 at times during or corresponding to injection of the calibration signal.
- clock dividers Divide by N 374-a and Divide by M 374-b, divide the frequency of Sample Clock 374 by predefined values (e.g., N and M, respectively) to produce clock timing signals with reduced frequencies, as required by N-Stage Accumulator 322-b and N-point FFT (over M symbols) 324-a, respectively.
- FIG. 3F is a block diagram illustrating a Digital Processor 214, according to some embodiments.
- Digital Processor 214 can be used in a Global Navigation Satellite System (GNSS) application shown in Figure 1.
- Digital Processor 214 typically includes one or more processing units (CPU's) 602, one or more network or other communications interfaces 604, memory 610, optional supplemental receivers 670 (e.g., one or more receivers to receive navigation aiding signals, such as wide-area differential satellite positioning system signals), and one or more communication buses 609 for interconnecting these components.
- the communication buses 609 may include circuitry (sometimes called a chipset) that interconnects and controls communications between system components.
- Digital Processor 214 is coupled to a satellite positioning system receiver 674 (e.g., signal receiver 120 excluding Digital Processor 214) to receive digitized (and typically compensated) satellite signals for processing.
- Digital Processor 214 optionally includes a user interface 605 that includes a display device 606 and input devices 608 (e.g., keyboard, mouse, touch screen, keypads, etc.).
- Memory 610 includes high-speed random access memory, such as DRAM, SRAM, DDR RAM or other random access solid state memory devices; and optionally (but typically) includes non- volatile memory, such as one or more magnetic disk storage devices, optical disk storage devices, flash memory devices, or other non-volatile solid state storage devices.
- Memory 610 optionally includes one or more storage devices remotely located from CPU(s) 602.
- Memory 610 or alternately the non- volatile memory device(s) within memory 610, comprises a non-transitory computer readable storage medium.
- memory 610 or the computer readable storage medium of memory 610 stores the following programs, modules and data structures, or a subset thereof:
- an operating system 612 that includes procedures for handling various basic system services and for performing hardware dependent tasks
- a communication module 614 that is used for connecting Digital Processor 214 to other computers, wide-area differential satellite positioning systems, local positioning systems, etc., via the one or more communication interfaces 604 (wired or wireless) and one or more communication networks, such as the Internet, other wide area networks, local area networks, metropolitan area networks, and so on;
- a user interface module 616 that receives commands from the user via one or more input devices 608 of user interface 605, generates user interface objects in display device 606, and/or displays maps, coordinates, routes, etc., related to the position of Signal Receiver 120;
- a navigation module 618 that produces navigation results (e.g., a range to satellite, ranges to multiple satellites, geographic positioning, location information, and/or a time value) by processing digitized satellite signals received from satellite positioning system receiver 674;
- a calibration module 620 that controls when calibration signal generator 230 ( Figure 2) injects a calibration signal (e.g., the timing and/or duration of injection of the calibration signal) at the antenna interface 204 and (optionally) for controlling configuration of the compensation mechanism(s) 212; and
- stored calibration sets 622 that include multiple sets of filter coefficients (stored, for example, in a Filter Coefficients Library 358); optionally, stored calibration sets 622 correspond to previously performed calibrations (e.g., factor calibrations) for different operating conditions (e.g., for different temperature conditions).
- Each of the above identified elements may be stored in one or more of the previously mentioned memory devices, and corresponds to a set of instructions for performing a function described above.
- the set of instructions can be executed by one or more processors (e.g., the CPUs 602).
- the above identified modules or programs i.e., sets of instructions
- memory 610 may store a subset of the modules and data structures identified above.
- memory 610 may store additional modules and data structures not described above.
- Figure 3F shows a "Digital Processor"
- Figure 3F is intended more as functional description of the various features which may be present in a mobile satellite navigation receiver than as a structural schematic of the embodiments described herein.
- items shown separately could be combined and some items could be separated.
- FIG. 4A is a flow diagram illustrating the Configuring Mechanism 320 described previously in relation to Figure 3B.
- the phase values of corresponding frequency components of the digitized received signals represent the measured phases (e.g., Actual Phase Response 402) of the filtered signal at various frequencies present in the calibration signal.
- a smoothed version of the measured phase response (Actual Phase Response 402, Figure 4A) is obtained from a polynomial fit (e.g., polynomial smoothing) to the measured phases of the filtered signal (e.g., the digitized received signal received in response to the injected calibration signal) at various frequencies.
- the smoothed measured phase response is represented (e.g., approximated) by a polynomial of a numerical order selected to eliminate undesired noise in the polynomial representation of the measured phases of the filtered signal, while
- target phase values are previously stored.
- the target phase values e.g., Target Response Curve 326 or Target Response 404
- the measured phases e.g., Actual Phase Response 402
- a straight line fit e.g., Straight Line 404-a
- the residual phase differences e.g., Residuals 328 or
- Phase Residuals 406 are computed by subtracting the phase values of corresponding frequency components of the digitized received signals from the target phase values. As shown in Figure 4A, Phase Residuals 406 are computed by subtracting Actual Phase
- Figures 4B-4C are flow diagrams illustrating operation of the Configuring
- the magnitude values of corresponding frequency components of the digitized received signals represent the measured magnitudes (e.g., Actual Magnitude Response 412 shown in Figure 4B and Actual Magnitude Response 422 shown in Figure 4C) of the filtered signal at various frequencies present in the calibration signal.
- target magnitude values are previously stored.
- the target magnitude values e.g., Target Response Curve 326 or Target Magnitude Response 414 shown in Figure 4B or Target Magnitude Response 424 shown in Figure 4C
- the measured magnitudes e.g., Actual Magnitude Response 412
- the magnitude response residuals are computed by determining amplitude ratios between one or more target amplitude values and amplitude values of corresponding frequency components of the digitized received signals.
- Magnitude Residuals 406 and Magnitude Residuals 416 are computed by dividing Target Magnitude Response 414 (or 424, Figure 4C) by Actual Magnitude Response 412 (or 422, Figure 4C).
- the complex filter coefficients (e.g., Complex FIR coefficients 334 or FIR Coefficients 408, FIR Coefficients 418, or FIR Coefficients 428) of the FIR filter (e.g., Complex FIR 314) correspond to residual phase differences (e.g., Residuals 328 or Phase Residuals 406) between target phase values (e.g., Target Response Curve 326 or Target Response 404) and phase values (e.g., Actual Phase Response 402) of corresponding frequency components of the digitized received signals (e.g., Digitized Received Signals 312), and residual magnitude ratios (e.g., Residuals 328 or Magnitude Residuals 416 or Magnitude Residuals 426) between target magnitude values (e.g., Target Response Curve 326 or Target Magnitude Response 414 or Target Magnitude Response 424) and magnitude values (e.g., Actual Magnitude Response 412 or Actual Magnitude Response 422) of corresponding frequency components
- residual phase differences e.g.,
- An inverse FFT (e.g., iFFT 330) is performed on the residual phase differences (e.g., Residuals 328 or Phase Residuals 406) and magnitude response residuals (e.g., Magnitude Residuals 416 or Magnitude Residuals 426), which results in the complex filter coefficients (e.g., Complex FIR coefficients 334 or FIR Coefficients 408) of the FIR filter (e.g., Complex FIR 314).
- the iFFT produces an initial set of complex filter coefficients, which are then multiplied by a windowing function (e.g., a Tukey window) to produce the complex filter coefficients used by an FIR filter in the Compensation
- Mechanism 212 to filter digitized received signals to produce compensated digitized received signals.
- Figures 5A-5C and Figure 5D include a flowchart representing method 500 and method 570 respectively for signal compensation performed at a signal receiver (e.g., Signal Receiver 120), according to certain embodiments of the invention.
- Method 500 and method 570 may be governed by instructions that are stored in a computer readable storage medium and that are executed by one or more processors of one or more servers.
- Each of the operations shown in Figures 5A-5D may correspond to instructions stored in a computer memory or computer readable storage medium.
- the computer readable storage medium may include a magnetic or optical disk storage device, solid state storage devices such as Flash memory, or other non-volatile memory device or devices.
- the computer readable instructions stored on the computer readable storage medium are in source code, assembly language code, object code, or other instruction format that is interpreted by one or more processors.
- a signal receiver e.g., a signal receiver
- Signal Receiver 120 receives (510) signals from an antenna (e.g., Antenna 202).
- the signal receiver processes (512) the received signals with analog signal processing circuitry (e.g., Analog Signal Processing Circuitry 208) to produce filtered signals.
- the signal receiver samples (514) the filtered signals with sampling circuitry (e.g., Sampling Circuitry 210) to produce digitized received signals (e.g., Digitized Received Signals 312).
- the signal receiver compensates (515) the digitized received signals, with a compensation mechanism (e.g., Compensation Mechanism 212), for non-uniform group delay and amplitude distortion introduced by the analog signal processing circuitry to produce compensated digitized received signals (e.g., Compensated Digitized Received Signals 316).
- a compensation mechanism e.g., Compensation Mechanism 212
- compensated digitized received signals e.g., Compensated Digitized Received Signals 316.
- compensating the digitized received signals includes (516) compensating for signal distortions in a predefined range of frequencies.
- the predefined range of frequencies optionally includes one or more of: (LI, El) 1559 MHz - 1591 MHz, 1559 MHz to 1610 MHz, (L2, L5) 1211 MHz - 1243 MHz, and 1160 MHz - 1192 MHz.
- the amplitude distortion comprises (517) deviation of a magnitude response of the analog signal processing circuitry from a uniform magnitude response across a predefined frequency range or a magnitude response curve that is symmetric with respect to a predefined center frequency of the predefined frequency range.
- target magnitude response curve e.g., Target Magnitude Response 414 and Target Magnitude Response 424 are symmetric with respect to a predefined center frequency of the predefined frequency range.
- the signal receiver processes (518) the compensated digitized received signals so as to produce a result (e.g., Result 220).
- the result includes (519) a navigation result.
- the result includes a range to satellite, ranges to multiple satellites, navigation result(s), geographical location(s), and/or satellite time value(s).
- the signal receiver configures (520) the compensation mechanism to process the digitized received signals using a finite impulse response filter (FIR filter), such as Complex FIR 314 ( Figure 3B), with complex filter coefficients (e.g., Complex FIR coefficients 334 shown in Figure 3B) to compensate for non-uniform group delay and amplitude distortion introduced by at least the analog signal processing circuitry.
- FIR filter finite impulse response filter
- the set of complex filter coefficients of the FIR filter correspond to residual phase differences between target phase values and phase values of corresponding frequency components of the digitized received signals and to amplitude ratios between one or more target amplitude values and amplitude values of corresponding frequency
- FIR Coefficients 408 are obtained from an inverse transform (iFFT 330) of Phase Residuals 406 and Magnitude Residuals 416 (or Magnitude Residuals 426).
- Phase Residuals 406 are the differences between (i.e., obtained from subtracting) Target Response 404 and Actual Phase Response 402.
- Magnitude Residuals 416 are the ratios between (i.e., obtained from dividing) Target Magnitude Response 414 and Actual
- Signal Receiver 120 performs additional steps for configuring (524) Compensation Mechanism 212 as described in steps 526-550 described below.
- the signal receiver injects (526) a calibration signal at the antenna interface
- the calibration signal generator injects (528) the calibration signal at the antenna interface at a duty cycle of less than a predefined fraction (e.g., one percent).
- the calibration signal includes (530) a time domain signal corresponding to orthogonal frequency components in a predefined range of frequencies.
- the calibration signal includes (532) a step carrier signal.
- the calibration signal includes (534) a swept carrier signal.
- the calibration signal has (536) an amplitude at least a predefined amount (e.g., 40 dB) greater than an amplitude of the received signals from the signal receiver's antenna (e.g., Antenna 202).
- a predefined amount e.g. 40 dB
- an amplitude of the received signals from the signal receiver's antenna e.g., Antenna 202).
- the amplitude and duration of the calibration signal is determined (538) in accordance with a criterion that the average power of the calibration signal is maintained lower than a value which would impact gain settings of a gain control circuit of the signal receiver. Operation of the calibration signal generator is further described with reference to Figure 3D above.
- the signal receiver provides (540) synchronized clock timing signals from a clock (e.g., Clock 240), to the signal receiver's calibration signal generator and sampling circuitry.
- a clock e.g., Clock 240
- Signal Receiver 120 configures (542)
- Compensation Mechanism 212 in accordance with signals received by Compensation Mechanism 212 in response to injection of the calibration signal at Antenna Interface 204.
- configuring Compensation Mechanism 212 includes (543) coherently accumulating digitized received signals corresponding to the injected calibration signals.
- Coherent Sample Accumulator 322 (including two or more registers) is used to coherently accumulate digitized received signals corresponding to the injected calibration signals.
- Signal Receiver 120 stores (544) control instructions at
- control instructions include (546) instructions for controlling when Calibration Signal Generator 230 injects the calibration signal at Antenna Interface 204.
- control instructions include (548) instructions for controlling configuration of Compensation Mechanism 212 in accordance with the digitized received signals received in response to the injection of the calibration signal at Antenna Interface 204.
- Signal Receiver 120 configures (550) Compensation
- Mechanism 212 by performing additional, or alternative, steps 552-564 described next, for example using Configuring Mechanism 320 (shown and described above with reference to Figure 3B).
- Signal Receiver 120 transforms (552) the digitized received signals, received in response to injection of the calibration signal, from time-domain to frequency-domain to produce a complex-valued frequency-domain representation of the digitized received signals. It may be noted that a real- valued representation is a special case of a complex- valued representation.
- the transformation from time- domain to frequency-domain is a Discrete Fourier Transform implemented using, for example, a Fast Fourier Transform (FFT) implementation. As described in relation to the example embodiment in Figures 3B and 3E, an FFT 324 is be used to produce N Complex Frequency Domain Samples over M-Symbols 378.
- FFT Fast Fourier Transform
- Signal Receiver 120 extracts (554) magnitude (e.g.,
- phase (e.g., Actual Phase Response 402, Figure 4A) of the digitized received signals are eliminated, for example by methods such as "phase unwrapping" (e.g., by the addition or subtraction of integer multiples of 2 ⁇ radians).
- phase unwrapping e.g., by the addition or subtraction of integer multiples of 2 ⁇ radians.
- the phase (e.g., Actual Phase Response 402, Figure 4 A) of the digitized received signals is converted to a continuous function of frequency by eliminating discontinuities introduced, for example, by phase-wrapping.
- Signal Receiver 120 obtains (556) target frequency-dependent phase information and target magnitude response information.
- target frequency-dependent phase information (Target Response Curve 326 in Figure 3B or Target Response 404 in Figure 4A) corresponds to a desired set of frequency-dependent phase values stored previously or generated from phase values (e.g., Actual Phase Response 402, Figure 4A) of corresponding frequency components.
- the target (also referred to herein as "desired") frequency-dependent phase information is a polynomial (for example, linear) relationship between phase and frequency.
- a polynomial relationship between phase and frequency is obtained, for example, using a polynomial fit (e.g., a straight line fit such as Straight Line 404-a shown in Figure 4A) to the (actual) extracted phase information (e.g., Actual Phase Response 402, Figure 4A) or by using an arbitrary polynomial fit (e.g., a straight line fit such as Straight Line 404-a shown in Figure 4A) to the (actual) extracted phase information (e.g., Actual Phase Response 402, Figure 4A) or by using an arbitrary
- target frequency-dependent phase information is obtained either by generating this information using curve fitting (or the like) on phase information extracted during calibration, or by retrieving predefined information, such as a "target" polynomial or straight line.
- target magnitude response information (Target Response Curve 326 in Figure 3B or Target Magnitude Response 414 in Figure 4B) corresponds to a desired set of frequency-dependent magnitude values stored previously or generated from magnitude values (e.g., Actual Magnitude Response 412, Figure 4B) of corresponding frequency components.
- the target (also referred to herein as "desired") magnitude response information is a polynomial (for example, symmetric) relationship between phase and frequency.
- Signal Receiver 120 computes (558) phase residuals by subtracting the extracted phase values from the target frequency-dependent phase
- Phase Residuals 406 are computed by subtracting Actual Phase Response 402 from Target Response 404.
- Signal Receiver 120 computes (559) magnitude response residuals by determining ratios of the extracted magnitude values with the target magnitude response information.
- Magnitude Residuals 416 are computed by dividing Actual
- Signal Receiver 120 generates (560) a set of complex values from the computed phase residuals and magnitude response residuals. For example, if a computed phase residual is ' ⁇ ' and the computed magnitude response residual is 'R' for a given frequency, then the complex value generated for that frequency corresponding to magnitude 'R' and phase residual ' ⁇ ' would be of the form R x (cos ⁇ + j sin ⁇ ).
- Signal Receiver 120 transforms (562) the set of complex values (e.g., complex filter response) from frequency-domain to time-domain to produce complex- valued time-domain filter coefficients (e.g., FIR Coefficients 408, Figure 4A) corresponding to a finite impulse response filter (e.g., Complex FIR 314, Figure 3B).
- the transformation from frequency-domain to time-domain is an Inverse Discrete Fourier Transform, implemented, for example, through an inverse Fast Fourier Transform (iFFT) implementation (e.g., iFFT 330, Figure 3B).
- iFFT inverse Fast Fourier Transform
- Signal Receiver 120 multiplies (564) initial complex- valued time-domain filter coefficients with a windowing function (e.g., a Tukey window) to produce the complex-valued time-domain filter coefficients.
- the step of transforming includes multiplying an initial set of complex- valued time-domain filter coefficients by a windowing function (e.g., Windowing 332, Figure 3B) to produce the complex-valued time-domain filter coefficients (e.g., FIR Coefficients 408, Figure 4A) corresponding to a finite impulse response filter (e.g., Complex FIR 314, Figure 3B).
- a windowing function e.g., Windowing 332, Figure 3B
- the complex-valued time-domain filter coefficients e.g., FIR Coefficients 408, Figure 4A
- the resulting coefficients are used by the FIR filter during normal operation to compensate for group delay distortions at the signal receiver.
- a signal receiver receives (572) signals.
- the signal receiver processes (574) the received signals at the signal receiver to produce digitized received signals.
- the signal receiver configures (576) a compensation mechanism (for example, but not limited to Compensation Mechanism 212) at the signal receiver in accordance with digital signals received by the compensation mechanism in response to injection of a calibration signal at the signal receiver by a calibration signal generator (for example, but not limited to Calibration Signal Generator 230).
- the configuring includes (578) configuring the compensation mechanism to compensate for non-uniform group delay and amplitude distortion introduced by analog signal processing circuitry of the signal receiver.
- the amplitude distortion comprises (579) deviation of a magnitude response of the analog signal processing circuitry from a uniform magnitude response across a predefined frequency range or a magnitude response curve that is symmetric with respect to a predefined center frequency of the predefined frequency range.
- target magnitude response curve e.g., Target Magnitude Response 414 and Target Magnitude Response 424.
- configuring includes (580) coherently accumulating digitized received signals corresponding to the injected calibration signals.
- the signal receiver configures (582) the compensation mechanism by performing additional steps (e.g., all or a subset of steps 552-564, Figure 5C) described above with regard to method 500.
- the step (576) of configuring the compensation mechanism is performed at various instances in time (e.g., intermittently and/or periodically), for example based on control signals provided by a digital processor at the signal receiver.
- the signal receiver compensates (582) the digitized received signals using the compensation mechanism to produce compensated signals.
- compensating the digitized received signals includes (517) compensating for signal distortions in a predefined range of frequencies that includes (LI, El) 1559 MHz - 1591 MHz, or 1559 MHz to 1610 MHz, or (L2, L5) 1211 MHz - 1243 MHz, or 1160 MHz - 1192 MHz.
- the signal receiver processes (586) the compensated signals so as to produce a result.
- the result includes (588) a navigation result.
- the result includes a range to satellite, ranges to multiple satellites, navigation result(s), geographical location(s), and/or satellite time value(s).
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Noise Elimination (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Radio Relay Systems (AREA)
- Circuits Of Receivers In General (AREA)
Abstract
Description
Claims
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
BR112014025626-8A BR112014025626B1 (en) | 2012-06-08 | 2013-06-06 | SIGNAL RECEIVER |
EP13801209.1A EP2859660B1 (en) | 2012-06-08 | 2013-06-06 | Signal receiver with group delay and amplitude distortion compensation |
CN201380023270.8A CN104380610B (en) | 2012-06-08 | 2013-06-06 | There is group delay and the signal receiver of amplitude distortion compensation |
JP2015516215A JP2015526936A (en) | 2012-06-08 | 2013-06-06 | Signal receiver with group delay and amplitude distortion compensation |
AU2013271536A AU2013271536B2 (en) | 2012-06-08 | 2013-06-06 | Signal receiver with group delay and amplitude distortion compensation |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261657650P | 2012-06-08 | 2012-06-08 | |
US61/657,650 | 2012-06-08 | ||
US13/892,210 US8837654B2 (en) | 2012-06-08 | 2013-05-10 | Signal receiver with group delay and amplitude distortion compensation |
US13/892,210 | 2013-05-10 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013184933A1 true WO2013184933A1 (en) | 2013-12-12 |
Family
ID=51390271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2013/044541 WO2013184933A1 (en) | 2012-06-08 | 2013-06-06 | Signal receiver with group delay and amplitude distortion compensation |
Country Status (6)
Country | Link |
---|---|
EP (1) | EP2859660B1 (en) |
JP (1) | JP2015526936A (en) |
CN (1) | CN104380610B (en) |
AU (1) | AU2013271536B2 (en) |
BR (1) | BR112014025626B1 (en) |
WO (1) | WO2013184933A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108141275A (en) * | 2015-09-15 | 2018-06-08 | 高通股份有限公司 | For the power control scheme based on duty ratio of satellite communication |
CN112004263A (en) * | 2019-05-27 | 2020-11-27 | 中兴通讯股份有限公司 | Group delay fluctuation calibration method and device |
CN116136602A (en) * | 2023-04-14 | 2023-05-19 | 福建福大北斗通信科技有限公司 | Device and method for in-band spectrum amplitude and time delay consistency of Beidou anti-interference channel |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107306241A (en) * | 2016-04-18 | 2017-10-31 | 中兴通讯股份有限公司 | A kind of method, device and the signal transmitter of signal precompensation |
DE102016117587B3 (en) * | 2016-09-19 | 2018-03-01 | Infineon Technologies Ag | CIRCUIT ARRANGEMENT WITH OPTIMIZED FREQUENCY TRANSITION AND METHOD FOR CALIBRATING A CIRCUIT ARRANGEMENT |
KR102077000B1 (en) * | 2018-01-29 | 2020-04-07 | 주식회사 만도 | Apparatus and Method for compensating refection loss of antenna for Radar, and Radar Apparatus using the same |
CN109217843A (en) * | 2018-09-28 | 2019-01-15 | 西安空间无线电技术研究所 | A kind of asymmetric FIR distortion compensation filter design method of satellite launch channel dual domain |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999060701A1 (en) * | 1998-05-18 | 1999-11-25 | Technische Universiteit Delft | Method and device for filtering a digital signal with fractional delay |
US20020137510A1 (en) | 2001-01-10 | 2002-09-26 | Hughes Electronics | Method and system for adaptive equalization for receivers in a wide-band satellite communications system |
US20020168951A1 (en) * | 2001-01-12 | 2002-11-14 | Tod Paulus | Notch filter for DC offset reduction in radio-frequency apparatus and associated methods |
EP2086194A2 (en) * | 2008-02-04 | 2009-08-05 | Her Majesty the Queen in Right of Canada, as represented by the Minister of Industry, through The Communications Research Centre Canada | Distortion Compensation Circuit and Method |
WO2010097349A2 (en) | 2009-02-27 | 2010-09-02 | Astrium Limited | Compensation apparatus |
US20110280570A1 (en) * | 2010-05-14 | 2011-11-17 | Verizon Patent And Licensing Inc. | Compensating for end-to-end group delay ripples |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7333423B2 (en) * | 2004-03-31 | 2008-02-19 | Intel Corporation | Transceiver with calibrated I and Q paths and methods for deconvolved calibration |
JP4648948B2 (en) * | 2005-06-24 | 2011-03-09 | パナソニック株式会社 | Wireless receiver |
CA2576778C (en) * | 2006-02-07 | 2014-09-02 | Xinping Huang | Self-calibrating multi-port circuit and method |
CN101964449A (en) * | 2010-08-27 | 2011-02-02 | 中国科学院上海微系统与信息技术研究所 | On-track correcting device of space-borne phased array transmitting antenna |
-
2013
- 2013-06-06 JP JP2015516215A patent/JP2015526936A/en active Pending
- 2013-06-06 EP EP13801209.1A patent/EP2859660B1/en active Active
- 2013-06-06 BR BR112014025626-8A patent/BR112014025626B1/en active IP Right Grant
- 2013-06-06 WO PCT/US2013/044541 patent/WO2013184933A1/en active Application Filing
- 2013-06-06 AU AU2013271536A patent/AU2013271536B2/en active Active
- 2013-06-06 CN CN201380023270.8A patent/CN104380610B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999060701A1 (en) * | 1998-05-18 | 1999-11-25 | Technische Universiteit Delft | Method and device for filtering a digital signal with fractional delay |
US20020137510A1 (en) | 2001-01-10 | 2002-09-26 | Hughes Electronics | Method and system for adaptive equalization for receivers in a wide-band satellite communications system |
US20020168951A1 (en) * | 2001-01-12 | 2002-11-14 | Tod Paulus | Notch filter for DC offset reduction in radio-frequency apparatus and associated methods |
EP2086194A2 (en) * | 2008-02-04 | 2009-08-05 | Her Majesty the Queen in Right of Canada, as represented by the Minister of Industry, through The Communications Research Centre Canada | Distortion Compensation Circuit and Method |
WO2010097349A2 (en) | 2009-02-27 | 2010-09-02 | Astrium Limited | Compensation apparatus |
US20110280570A1 (en) * | 2010-05-14 | 2011-11-17 | Verizon Patent And Licensing Inc. | Compensating for end-to-end group delay ripples |
Non-Patent Citations (1)
Title |
---|
See also references of EP2859660A4 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108141275A (en) * | 2015-09-15 | 2018-06-08 | 高通股份有限公司 | For the power control scheme based on duty ratio of satellite communication |
CN108141275B (en) * | 2015-09-15 | 2020-12-01 | 高通股份有限公司 | Satellite communication power control method and device based on duty ratio |
CN112004263A (en) * | 2019-05-27 | 2020-11-27 | 中兴通讯股份有限公司 | Group delay fluctuation calibration method and device |
CN112004263B (en) * | 2019-05-27 | 2023-05-30 | 中兴通讯股份有限公司 | Group delay fluctuation calibration method and device |
CN116136602A (en) * | 2023-04-14 | 2023-05-19 | 福建福大北斗通信科技有限公司 | Device and method for in-band spectrum amplitude and time delay consistency of Beidou anti-interference channel |
Also Published As
Publication number | Publication date |
---|---|
EP2859660A4 (en) | 2016-01-20 |
BR112014025626A2 (en) | 2017-06-20 |
BR112014025626B1 (en) | 2022-08-16 |
AU2013271536A1 (en) | 2014-08-07 |
EP2859660A1 (en) | 2015-04-15 |
JP2015526936A (en) | 2015-09-10 |
AU2013271536B2 (en) | 2017-08-10 |
EP2859660B1 (en) | 2018-11-21 |
CN104380610B (en) | 2016-08-17 |
CN104380610A (en) | 2015-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2013271530B2 (en) | Signal receiver with group delay compensation | |
AU2013271536B2 (en) | Signal receiver with group delay and amplitude distortion compensation | |
US8837654B2 (en) | Signal receiver with group delay and amplitude distortion compensation | |
CN104185974B (en) | Apparatus and method for the real-time I/Q disequilibrium regulatings of broadband RF receiver | |
US9118401B1 (en) | Method of adaptive interference mitigation in wide band spectrum | |
US8649464B2 (en) | Quadrature receiver and method of compensating for I/Q imbalance using a calibration signal | |
US9356732B2 (en) | Quadrature error detection and correction | |
US20140348018A1 (en) | Self-interference cancellation | |
WO2014145350A1 (en) | Quadrature error detection and correction | |
CN102098258B (en) | Method for removing narrow-band interference and self-adapting filter | |
EP2894823B1 (en) | Coefficient estimation for digital IQ calibration | |
US8836552B1 (en) | Direct input radio frequency complex analog to digital converter with corrective factors | |
JP2005348236A (en) | Array antenna transmitter and receiver, and calibration method | |
CN104065598B (en) | Broadband IQ disequilibrium regulatings method, apparatus and system | |
WO2014163640A1 (en) | Vswr estimation using spectral analysis to suppress external interference | |
Soltanian et al. | Utilization of multi-rate signal processing for GNSS-SDR receivers | |
US20240187025A1 (en) | Feedforward distortion cancellation apparatus and method | |
EP2978183A1 (en) | Filtering a modulated signal | |
SG185841A1 (en) | Method and system for performing complex sampling of signals by using two or more sampling channels and for calculating time delays between these channels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13801209 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2013271536 Country of ref document: AU Date of ref document: 20130606 Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2013801209 Country of ref document: EP |
|
ENP | Entry into the national phase |
Ref document number: 2015516215 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: 112014025626 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: 112014025626 Country of ref document: BR Kind code of ref document: A2 Effective date: 20141014 |