WO2005018087A2 - Procedure pour estimer un parametre de decalage de frequence d'un oscillateur - Google Patents

Procedure pour estimer un parametre de decalage de frequence d'un oscillateur Download PDF

Info

Publication number
WO2005018087A2
WO2005018087A2 PCT/US2004/015205 US2004015205W WO2005018087A2 WO 2005018087 A2 WO2005018087 A2 WO 2005018087A2 US 2004015205 W US2004015205 W US 2004015205W WO 2005018087 A2 WO2005018087 A2 WO 2005018087A2
Authority
WO
WIPO (PCT)
Prior art keywords
oscillator
correlation functions
signal
frequency
shift parameter
Prior art date
Application number
PCT/US2004/015205
Other languages
English (en)
Other versions
WO2005018087A3 (fr
Inventor
Norman F. Krasner
Original Assignee
Qualcomm Incorporated
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qualcomm Incorporated filed Critical Qualcomm Incorporated
Priority to BRPI0413053-7A priority Critical patent/BRPI0413053A/pt
Publication of WO2005018087A2 publication Critical patent/WO2005018087A2/fr
Publication of WO2005018087A3 publication Critical patent/WO2005018087A3/fr
Priority to IL173264A priority patent/IL173264A/en

Links

Classifications

    • 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
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/16—Circuits
    • H04B1/26—Circuits for superheterodyne receivers
    • H04B1/28—Circuits for superheterodyne receivers the receiver comprising at least one semiconductor device having three or more electrodes
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06—Receivers
    • H04B1/16—Circuits
    • H04B1/30—Circuits for homodyne or synchrodyne receivers

Definitions

  • This invention relates to oscillators, and, more specifically, to procedures for estimating a frequency shift parameter of an oscillator or enhancing the frequency stability of an oscillator.
  • one or more local oscillators are used to down-convert a received signal to baseband or near-baseband frequencies.
  • a single local oscillator having a frequency set to the carrier frequency of the received signal is used to perform the down-conversion process.
  • a first local oscillator is used to down-convert the received signal to an intermediate frequency signal
  • a second local oscillator is used to down-convert the received signal to a baseband or near-baseband signal.
  • the frequency of the first local oscillator is set to the difference in frequency between the carrier frequency and the intermediate frequency.
  • the frequency of the second local oscillator is set to the intermediate frequency.
  • the one or more local oscillators may experience frequency instability and rapid frequency shifts, due to various factors such as power on/off conditions, voltage fluctuations and the like. Depending on the application, these frequency changes can cause significant problems in processing the down-converted signal.
  • frequency shifts experienced by the one or more local oscillators used to perform down-conversion will manifest themselves as phase changes in the down-converted signal. These phase changes are undesirable because they limit the size of the coherent integration time, which in turn limits the sensitivity of the GPS position determination system and its ability to detect weak GPS satellite signals.
  • the local oscillator may experience a linear frequency rate of change.
  • This linear frequency rate of change may be expressed as a sweep rate k in units of Hz/sec. Since phase is the integral of frequency, a sweep rate k manifests itself as a quadratic rate of change in phase, i.e., on the order of k 2 . If this sweep rate is known, it can be controlled and corrected for. Unfortunately, for mass market consumer electronic devices such as subscriber stations, it is too expensive and time-consuming to attempt to estimate k using a maximum likelihood approach of checking every possible hypothesis of the parameter k. Other possible techniques, such as homodyne spectral analysis, are ineffective at the low SNRs which characterize many GPS satellite transmissions.
  • a system for estimating a frequency shift parameter of an oscillator is described.
  • a received signal is down-converted using the oscillator.
  • First and second correlation functions are derived from the down-converted signal.
  • the first and second correlation functions represent the likelihood of one or more frequency hypotheses at first and second times, respectively.
  • An estimate of the frequency shift parameter is derived from the first and second correlation functions.
  • the first and second correlation functions are circularly cross-correlated with one another to form a circular cross-correlation function.
  • the peak of the circular cross-correlation function is located, and an estimate of the frequency shift parameter is derived from the location of the peak of the circular cross- correlation function.
  • the frequency shift parameter is the sweep rate, and the first and second times are displaced from one another by a delta value.
  • the sweep rate is estimated from the location of the peak of the circular cross-correlation function and the delta value.
  • the received signal comprises a composite signal representing one or more signal sources.
  • First and second correlation functions are derived from the down-converted signal for each of one or more of the signal sources.
  • the first and second correlation functions represent the likelihood of one or more frequency hypotheses at first and second times, respectively.
  • One or more of the first correlation functions are selectively combined using a quality of signal metric, e.g., SNR, to form a first combined correlation function.
  • One or more of the second correlation functions are selectively combined, again using a quality of signal metric, to form a second combined correlation function.
  • An estimate of the frequency shift parameter is derived from the first and second combined correlation functions.
  • Figure 1 is a block diagram of a system for estimating a frequency shift parameter of an oscillator.
  • Figure 2 is a block diagram of s system for estimating a frequency shift parameter of a local oscillator in a position determination system.
  • Figure 3 is an example of a two-dimensional correlation function as produced by a subscriber station in a position determination system.
  • Figure 4 is a one-dimensional correlation function derived from a cross- section of the two-dimensional correlation function of Figure 3, and representing the likelihood of various code phase hypotheses
  • Figure 5 is a one-dimensional correlation function derived from a cross- section of the two-dimensional correlation function of Figure 3, and representing the likelihood of various Doppler frequency hypotheses .
  • Figure 6 A is an example of first and second correlation functions representing the likelihood of various frequency hypotheses at first and second times, respectively, and assuming a sweep rate of 50 Hz/sec.
  • Figure 6B illustrates a circular cross-correlation function derived from the first and second correlation functions of Figure 6 A.
  • Figure 7A is an example of first and second correlation functions representing the likelihood of various frequency hypotheses at first and second times, respectively, and assuming a sweep rate of 100 Hz/sec.
  • Figure 7B illustrates a circular cross-correlation function derived from the first and second correlation functions of Figure 7 A.
  • Figure 8 is a series of plots illustrating the expected RMS error in sweep rate estimates for various input signal levels and estimated sweep rates.
  • Figure 9 is a block diagram of one embodiment of a subscriber station.
  • Figure 10 is a flowchart of one embodiment of a method for estimating a frequency shift parameter of an oscillator.
  • Figure 11 is a flowchart of one implementation of a method for estimating a frequency shift parameter of an oscillator in which the received signal is a composite signal representing a plurality of signal sources, and the plurality of correlation functions corresponding to the various signal sources are selectively combined in the course of estimating the frequency shift parameter.
  • Figure 12 is a block diagram of an embodiment of a system for estimating a frequency shift parameter of an oscillator in which a processor executes software instructions held in a memory.
  • the term "software” includes source code, assembly language code, binary code, firmware, macro-instructions, micro-instructions, or the like, or any combination of two or more of the foregoing.
  • memory refers to any processor-readable medium, including but not limited to RAM, ROM, EPROM, PROM, EEPROM, disk, floppy disk, hard disk, CD-ROM, DVD, or the like, or any combination of two or more of the foregoing, on which may be stored a series of software instructions executable by a processor.
  • processor or “CPU” refer to any device capable of executing a series of instructions and includes, without limitation, a general- or special-purpose microprocessor, finite state machine, controller, computer, digital signal processor
  • logic refers to implementations of functionality in hardware, software, or any combination of hardware and software.
  • GPS satellite includes space vehicles (SVs).
  • position determination system means any system for determining the position of a subscriber station in a wireless communications system, and includes position determination systems overlaid onto or integrated into wireless communications systems. Examples include the GPS, mobile assisted GPS,
  • GLONASS GALILEO
  • EOTA enhanced observed time difference of arrival
  • TOA uplink time of arrival
  • Figure 1 illustrates a system 100 for estimating a frequency shift parameter of an oscillator 106.
  • a signal is received over antenna 102.
  • a down-converter 110 down-converts the received signal using the oscillator 106.
  • the down-converter 110 comprises mixer 104 in combination with oscillator 106, but it should be appreciated that other examples are possible.
  • Logic 108 derives first and second correlation functions from the down-converted signal 112.
  • the first and second correlation functions represent the likelihood of one or more hypothesized frequency parameters at first and second times, respectively.
  • Logic 108 derives an estimate of the frequency shift parameter from the first and second correlation functions.
  • the signal received over antenna 102 comprises a composite signal representing one or more signal sources.
  • Down-converter 110 down- converts the signal, and logic 108 derives, for each of one or more of the signal sources, first and second correlation functions.
  • the first and second correlation functions again represent the likelihood of one or more hypothesized frequency parameters at first and second times, respectively.
  • Logic 108 estimates the frequency shift parameter from the first and second correlation functions for one or more of the signal sources.
  • logic 108 combines one or more of the first correlation functions to form a first combined correlation function, and combines one or more of the second correlation functions to form a second combined correlation function.
  • Logic 108 derives the frequency shift estimate from the first and second combined correlation functions.
  • logic 108 selectively combines one or more of the first correlation functions, using a quality of signal metric, such as SNR, to form the first combined correlation function, and selectively combines one or more of the second correlation functions, using a quality of signal metric, again such as SNR, to form the second combined correlation function.
  • Figure 2 illustrates an embodiment 200 of a system for estimating a frequency shift parameter of an oscillator in a position determination system.
  • the system is embodied or integrated within a subscriber station, which in turn is part of a GPS position determination system, but it should be appreciated that other examples are possible.
  • the system 200 functions in two modes of operation, a normal mode, and a calibration mode.
  • the normal mode will be explained first followed by the calibration mode.
  • a signal is received over antenna 202.
  • the received signal is a composite signal representing one or more position determination signal sources visible to the subscriber station.
  • the one or more signal sources are each GPS satellites 204a, 204b, and 204c visible to the subscriber station's receiver in the presence of noise, but it should be appreciated that other examples are possible.
  • Each of these satellites transmits a signal marked with a repeating PN code of 1,023 chips uniquely identifying the satellite.
  • the PN codes for the satellites are shifted versions of one another.
  • the 1,023 chips repeat every millisecond.
  • the signal is also modulated with data bits, where each data bit has a 20 ms duration in the modulated signal.
  • the received signal is down-converted to an intermediate frequency by a down-converter 206.
  • the down-converter 206 comprises a mixer 208 in combination with a local oscillator 210, with the frequency of the local oscillator 210 set to the difference in frequency between the received signal and the intermediate frequency.
  • the down-converted signal 212 is passed through low pass filter 214 to eliminate the high order component of the down-converted signal 212, and then digitized by A/D converter 216.
  • the digitized signal 218 is passed through quadrature demodulator 228, which divides the digitized signal 218 into in-phase (I) and quadrature (Q) components.
  • the quadrature modulator 228 comprises mixer 220, local oscillator 222, 90° phase shifter 224, and mixer 226, but it should be appreciated that other examples are possible.
  • the frequency of the local oscillator 222 is set at about the intermediate frequency of the signal 218 to allow the de-modulator 228 to down-convert the signal 218 to baseband or near baseband frequencies.
  • Mixer 220 down-converts the signal 218 to baseband or near baseband frequencies using the output of the local oscillator 222, and low pass filter 230 filters out the high order component of this signal to form the I component 234.
  • Mixer 226 down- converts the signal 218 to baseband or near baseband frequencies using the 90° phase- shifted version of the output of local oscillator 222, and low pass filter 232 filters out the high order component of this signal to form the Q component 236.
  • PN generator 238 generates a PN code for one of the satellites. XOR gates
  • the de-spread components are each input to signal processing logic 244.
  • Signal processing logic 244 separately and coherently integrates N c successive values of each component, and non-coherently combines M successive coherent integrations to form a correlation value C.
  • Successive correlation values C are provided to microprocessor 246.
  • Microprocessor 246 forms two-dimensional correlation functions from the correlation values C provided to it from signal processing circuitry 244. As indicated by phantom line 248, microprocessor 246 controls the PN generator 238 to shift the PN code for a particular satellite within a search window W. As indicated by phantom lines 250 and 252, it also controls either or both local oscillators 210 and 222 to vary the frequency applied by either over a plurality of discrete frequency bins. It forms the correlation function by associating the correlation values which result by independently varying these two parameters.
  • Figure 3 is an example of a two dimensional correlation function which may be generated by this process.
  • the particular shape of this correlation function is implementation dependent, so nothing in this particular example should be taken as limiting.
  • the correlation function in this particular example is the energy in dB which results by varying the code phase, i.e., the PN code shift in Vi chip increments, over a search window W, and varying the frequency applied by either or both local oscillators over a plurality of discrete Doppler frequency bins.
  • the variation in code phase which is possible defines a code phase axis or dimension 302
  • the variation in Doppler frequency which is possible defines a Doppler frequency axis or dimension 304.
  • the correlation function forms three lobes, identified with numerals 306b, 308a and 308b, respectively.
  • the Doppler frequency dimension, the correlation function also forms three lobes, identified with numerals 30ba, 306b, and 306c, respectively.
  • Figure 4 illustrates a cross-section of the correlation function along axis 4-
  • the three lobes 306b, 308a, 308b in Figure 3 manifest themselves in Figure 4 as peaks 402a, 402b, and 402c, respectively.
  • Peak 402a is the important peak for position determination purposes.
  • Figure 5 illustrates a cross-section of the correlation function along axis 5-5.
  • the three lobes 306a, 306b, 306c manifest themselves in Figure 5 as peaks 502a, 502b, 502c. Peak 502a is the important peak for position determination purposes.
  • the microprocessor 246 locates the peaks of the correlation function along both the code phase and Doppler frequency bin dimensions. It derives an estimate of the time of arrival for the satellite from the location of the peak along the code phase dimension. It derives an estimate of any Doppler frequency shift (due to relative movement between the subscriber station and the satellite) from the location of the peak along the Doppler frequency bin dimension.
  • the microprocessor 246 repeats this process for one or more additional satellites by controlling PN generator 238 to apply the PN code for these additional satellites.
  • the microprocessor 246 derives three or more time of arrival measurements for three or more of the satellites. It then accesses (or has another network entity access) an almanac (not shown) to determine the precise locations of each of the satellites. It then estimates (or has the network entity estimate) the location of the subscriber station from this information using known trilateration or triangulation techniques.
  • the microprocessor 246 requires four or more time or arrival measurements, which allows determination of the x, y and z coordinates of the position of the subscriber station and system time t. Microprocessor 246 also estimates or has estimated any Doppler frequency shift due to relative movement between the subscriber station and one or more of the satellites.
  • microprocessor 246 is configured to derive correlation functions as before, except that the axis 304 no longer represents hypothesized Doppler frequency bins. Instead, the axis 304 represent hypothesized frequencies for the I and Q signals 234 and 236 in Figure 2.
  • One or both of the local oscillators 210, 222 are assumed to be susceptible to unwanted frequency shifts or rates of change, and the purpose of the special calibration mode of operation is to estimate any such shifts or rates of change.
  • the microprocessor 246 forms a two-dimensional correlation function for a particular satellite. Then it derives from the two-dimensional correlation function a one-dimensional correlation function along the frequency dimension. In one implementation, referring to Figure 4, the microprocessor 246 derives the one-dimensional correlation function by locating the code phase associated with the peak of the main lobe 402a, and then deriving the one-dimensional correlation function at this code phase along the frequency dimension. In one example, the microprocessor 246 locates the non-interpolated peak of the correlation function along the code phase dimension using, and then refines this peak location known quadratic interpolation techniques. The one-dimensional correlation function represents the likelihood of one or more frequency hypotheses at a time ti .
  • the microprocessor 246 also utilizes this same process to derive a second one-dimensional correlation function for the same satellite.
  • This second correlation function represents the likelihood of one or more frequency hypotheses at a time t 2 displaced from the time ti by a delta value. Both of these one-dimensional correlation functions may be represented as one or more discrete correlation values corresponding to one or more discrete frequency bins.
  • the microprocessor 246 circularly cross-correlates the one-dimensional correlation functions corresponding to the first and second times ti and t 2 . The location of the peak of this circular cross-correlation function is then located. In one example, the non-interpolated peak is located, and then known quadratic interpolation techniques are used to refine the location of this peak.
  • a sweep rate for the one or more local oscillators 210, 222 is then estimated from the location of this peak and the delta value , separating the times ti and t 2 .
  • the sweep rate estimate is derived from the ratio of the peak location (expressed in frequency) divided by the delta value separating the times ti and t 2 .
  • the result is an estimate of the sweep rate in terms of Hz/sec.
  • Figures 6A-6B, 7A-7B illustrate two examples of the process of estimating a sweep rate. Both assume a -130 dBm input signal, a 3dB noise figure, and a 2 dB implementation loss.
  • Figure 6A illustrates first and second correlation functions, identified with numerals 602 and 604 respectively, corresponding respectively to times ti and t 2 , with t being one second after ti.
  • Figure 6B illustrates the circular cross- correlation function which results from circularly cross-correlating the functions 602 and 604. Assuming each frequency bin corresponds to 25 Hz, the location 608 of the peak 606 along the frequency axis, 2, is 50 Hz, consistent with an estimated 50 Hz/sec. sweep rate.
  • Figure 7A illustrates first and second correlation functions, identified with numerals 702 and 704 respectively, corresponding respectively to times ti and t 2 , with t 2 again being one second after ti.
  • Figure 7B illustrates the circular cross-correlation function which results from circularly cross-correlating the functions 702 and 704. Assuming again each frequency bin corresponds to 25 Hz, the location 708 of the peak 706 along the frequency axis, 4, is 100 Hz, consistent with an estimated sweep rate of 100 Hz/sec.
  • Figure 8 illustrates simulation results showing the expected RMS error in the sweep rate estimate for different values of input signal level and estimated sweep rate. For example, assuming an input signal level of -144 dBm and an estimated sweep rate of 100 Hz/sec, the expected RMS error in the estimate of the sweep rate is 4 Hz/ sec.
  • the sweep rate may be stored and the calibration mode terminated.
  • the system 200 may then be allowed to enter or resume the normal mode of operation.
  • the microprocessor 246, armed with the estimated sweep rate may then control one or both of the local oscillators 210, 222 to counteract, in whole or in part, this sweep rate. More specifically, both are assumed to have a variable frequency under the control of the microprocessor 246 as indicated in phantom by the lines 250 and 252.
  • Microprocessor 246 may use the estimated sweep rate to adjust the frequency of one or both of these local oscillators until the effective sweep rate is zero or reduced to a predetermined amount.
  • microprocessor 246 is configured to derive first and second correlation functions for each of one or more GPS satellites.
  • the first correlation functions represent the likelihood of one or more frequency hypotheses at a time t ls and the second correlation functions represent the likelihood of one or more frequency hypotheses at a time t 2 displaced from ti by a delta value.
  • the microprocessor 246 then derives an estimate of the frequency shift parameter from the first and second correlation functions for one or more of the GPS satellites.
  • microprocessor 246 is configured to combine one or more of the first correlation functions to form a first combined correlation function, combine one or more of the second correlation functions to form a second combined correlation function, and then derive the estimate of the frequency shift parameter from the first and second combined correlation functions.
  • the microprocessor 246 is configured to selectively combine one or more of the correlation functions using a quality of signal metric, such as SNR.
  • the microprocessor may arrange the first correlation functions in descending order of SNR, and discard any waveforms whose non- interpolated peak differs from that of the first waveform by more than +/- 1 bin.
  • a weighted sum of the remaining correlation functions may then be formed, with the weighting applied to a given waveform being determined by the difference between the SNRs of the waveforms (e.g., a 3 dB difference in SNR implies a 2:1 weighting).
  • the same process may be employed to selectively combine the second correlation functions.
  • the first and second combined correlation functions may then circularly cross- correlated, and the frequency shift parameter derived from the circular cross-correlation function as described earlier.
  • the system 200 may be embodied or integrated within a subscriber station.
  • FIG. 9 is a block diagram of a subscriber station in which the system 200 may be embodied or incorporated.
  • Radio transceiver 906 is configured to modulate baseband information, such as voice or data, onto an RF carrier, and demodulate a modulated RF carrier to obtain baseband information.
  • Antenna 910 is configured to transmit a modulated RF carrier over a wireless communications link and receive a modulated RF carrier over a wireless communications link.
  • Baseband processor 908 is configured to provide baseband information from an input device within user interface 916 to transceiver 906 for transmission over a wireless communications link. Baseband processor 908 is also configured to provide baseband information from transceiver 906 to an output device within user interface 916.
  • User interface 916 comprises a plurality of devices for inputting or outputting user information such as voice or data. The devices typically included within the user interface include a keyboard, a display screen, a microphone, and a speaker.
  • GPS receiver 912 is configured to receive and down-convert GPS satellite transmissions, and provide the down-converted information to correlator 918.
  • Correlator 918 is configured to derive GPS correlation functions from the mformation provided to it by GPS receiver 912, and provide the same to CPU 902.
  • correlator 920 produces a correlation function which is defined over a code phase dimension and a frequency dimension.
  • the frequency dimension represents Doppler frequency hypotheses.
  • the frequency dimension represents hypothesized frequencies of the baseband or near baseband signals down-converted by the GPS receiver 912.
  • the correlation function is defined over a range of code phases which define a search window W.
  • the correlation function is defined over a plurality of frequency bins.
  • Correlator 918 may also be configured to derived pilot-related correlation functions from information relating to pilot signals provided to it by transceiver 906. This information is used by the subscriber station to acquire wireless communications services.
  • Channel decoder 920 is configured to decode channel symbols provided to it by baseband processor 908 into underlying source bits.
  • the channel decoder is a Viterbi decoder.
  • the channel decoder 920 is a turbo decoder.
  • Memory 904 in configured to hold software instructions for execution by
  • CPU 902. Through execution of these instructions, CPU 902 derives times and frequency measurements from the correlation functions generated by the correlator 908.
  • blocks 902, 904, 912, and 918 embody or incorporate the system 200 of Figure 2.
  • FIG 10 is a flowchart illustrating a first embodiment 1000 of a method of estimating a frequency shift parameter of an oscillator.
  • step 1002 comprises down-converting a received signal using the oscillator.
  • step 1004 comprises deriving first and second correlation functions from the down-converted signal. These first and second correlation functions represent the likelihood of one or more hypothesized frequency parameters at first and second times, respectively.
  • Step 1006 comprises deriving an estimate of the frequency shift parameter from the first and second correlation functions.
  • the method may be performed in a position determination system such as a
  • the method is performed in a subscriber station.
  • the received signal may be a composite signal representing one or more position determination signal sources, such as GPS satellites, visible to the subscriber station's receiver.
  • the first and second correlation functions in this implementation may correspond to the same position determination signal source, e.g., GPS satellite.
  • the one or more hypothesized frequency parameters comprise one or more discrete hypothesized frequency bins.
  • the deriving step comprises circularly cross-correlating the first and second correlation functions to form a circular cross-correlation function, locating a peak of the circular cross-correlation function, and deriving an estimate of the frequency shift parameter from the location of the peak of the circular cross-correlation function.
  • the oscillator may be a local oscillator in a down-converter within a receiver in the subscriber station. In one configuration, the oscillator has a variable frequency.
  • the estimated frequency shift parameter is determined and then stored.
  • the frequency of the oscillator is controllably varied using the stored estimated frequency shift parameter to counteract in whole or in part the effects of a frequency shift or sweep rate being experienced by the oscillator.
  • the first and second times are displaced from one another by a delta value
  • the estimated frequency shift parameter is sweep rate.
  • the deriving step comprises estimating the sweep rate of the oscillator from the location of the peak of the circular cross-correlation function and the delta value.
  • the received signal comprises a composite signal representing one or more signal sources.
  • first and second correlation functions are derived from the down-converted signal.
  • the first and second correlation functions represent the likelihood of one or more hypothesized frequency parameters at first and second times, respectively.
  • An estimate of the frequency shift parameter is then derived from the first and second correlation functions for one or more of the signal sources.
  • the method is performed in a subscriber station, and the one or more signal sources are GPS satellites visible to the subscriber station's receiver in the presence of noise.
  • the deriving step may comprise combining first correlation functions for each of one or more of the signal sources to form a first combined correlation function, combining second correlation functions corresponding to one or more of the signal sources to form a second combined correlation function, and then deriving an estimate of the frequency shift parameter from the first and second combined correlation functions.
  • FIG. 11 is a flowchart of one implementation 1100 of a method of estimating a frequency shift parameter of an oscillator.
  • Step 1102 comprises down- converting a received signal comprising a composite signal representative of one or more signal sources.
  • Step 1104 comprises deriving first and second correlation functions for each of one or more of the signal sources at first and second times, respectively.
  • Step 1106 comprises selectively combining the first correlation functions of one or more of the signal sources to form the first combined correlation function using a quality of signal metric.
  • Step 1108 comprises selectively combining the second correlation functions of one or more of the second signal sources to form the second combined correlation function using a quality of signal metric.
  • Step 1110 comprises deriving an estimate of the frequency shift parameter from the first and second combined correlation functions.
  • the quality of signal metric used to combine the first and second correlation functions is SNR.
  • the first and second correlation functions may each be arranged in decreasing order of SNR. Those waveforms in which non-interpolated peak differs from that of the initial waveform by more than +/- 1 bin are excluded.
  • the remaining first correlation functions are combined through a weighted averaging approach, in which the relative weight applied to a waveform in relation to another waveform depends on the difference in SNRs between the two waveforms.
  • the remaining second correlation functions are combined using the same process.
  • the deriving step may comprise circularly cross-correlating the first and second combined correlation functions to form a circular cross-correlation function, locating a peak of the circular cross-correlation function, and deriving the estimate of the frequency shift parameter from the location of the peak.
  • a memory holds software instructions embodying any of the foregoing methods, or embodiments, implementations, or examples thereof which have been described or suggested.
  • This memory may be part of a system such as that illustrated in Figure 12.
  • system 1200 comprises a processor 1202 and the above-described memory, identified in the figure with numeral 1204.
  • processor 1202 is configured to access the memory 1204 and execute the software instructions held therein.

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
  • Superheterodyne Receivers (AREA)

Abstract

L'invention concerne un système et un procédé pour estimer un paramètre de décalage de fréquence d'un oscillateur. Un signal reçu est converti vers le bas. Une première et une deuxième fonctions de corrélation sont dérivées du signal reçu. Les première et deuxième fonctions de corrélation représentent la probabilité d'une ou plusieurs hypothèses de fréquence à un premier et à un deuxième moments, respectivement. Une estimation d'un paramètre de décalage de fréquence est dérivée à partir des première et deuxième fonctions de corrélation.
PCT/US2004/015205 2003-07-28 2004-05-14 Procedure pour estimer un parametre de decalage de frequence d'un oscillateur WO2005018087A2 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
BRPI0413053-7A BRPI0413053A (pt) 2003-07-28 2004-05-14 processo para estimar um parámetro de deslocamento de freqüencia de um oscilador
IL173264A IL173264A (en) 2003-07-28 2006-01-19 Procedure for estimating a frequency shift parameter of an oscillator

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US49082103P 2003-07-28 2003-07-28
US60/490,821 2003-07-28

Publications (2)

Publication Number Publication Date
WO2005018087A2 true WO2005018087A2 (fr) 2005-02-24
WO2005018087A3 WO2005018087A3 (fr) 2005-04-14

Family

ID=34193091

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2004/015205 WO2005018087A2 (fr) 2003-07-28 2004-05-14 Procedure pour estimer un parametre de decalage de frequence d'un oscillateur

Country Status (4)

Country Link
KR (1) KR20060036104A (fr)
BR (1) BRPI0413053A (fr)
IL (1) IL173264A (fr)
WO (1) WO2005018087A2 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7275537B2 (en) 2000-08-10 2007-10-02 Meridica Limited Device for delivering physiologically active agent in powdered form
WO2012097358A1 (fr) * 2011-01-14 2012-07-19 Qualcomm Incorporated Mises en œuvre de traitement de signal sans fil

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4426712A (en) * 1981-05-22 1984-01-17 Massachusetts Institute Of Technology Correlation system for global position receiver

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4426712A (en) * 1981-05-22 1984-01-17 Massachusetts Institute Of Technology Correlation system for global position receiver

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7275537B2 (en) 2000-08-10 2007-10-02 Meridica Limited Device for delivering physiologically active agent in powdered form
WO2012097358A1 (fr) * 2011-01-14 2012-07-19 Qualcomm Incorporated Mises en œuvre de traitement de signal sans fil

Also Published As

Publication number Publication date
WO2005018087A3 (fr) 2005-04-14
KR20060036104A (ko) 2006-04-27
IL173264A (en) 2010-12-30
BRPI0413053A (pt) 2006-10-17
IL173264A0 (en) 2006-06-11

Similar Documents

Publication Publication Date Title
JP4855677B2 (ja) 位置決定システムのための信号サーチ手続き
US7002516B2 (en) Signal detector employing correlation analysis of non-uniform and disjoint sample segments
US7280071B2 (en) Method and apparatus for computing position and velocity of receiver using independently-derived sets of satellites
CN106291614B (zh) 用于在多路径环境中跟踪卫星无线电导航信号的装置
KR101398719B1 (ko) 조합된 내비게이션 신호의 처리 방법
US7127011B2 (en) Procedure for jammer detection
US7693211B2 (en) Fast fourier transform based phase locked loop for navigational receivers
EP1321774B1 (fr) Un procédé de reacquisition dans un récepteur de positionnement, et dispositif électronique correspondant
JP2004501352A (ja) 不均一で離散的なサンプル・セグメントを相関させるコヒーレント累積システムを採用する信号検出器及び方法
WO2009064812A1 (fr) Suppression des images fantômes pour un signal sps reçu
US7903026B2 (en) Positioning apparatus and positioning apparatus control method
EP1244225B1 (fr) Système et procédé pour la determination des bords d'un élément d'information
US7096132B2 (en) Procedure for estimating a parameter of a local maxima or minima of a function
WO2005018087A2 (fr) Procedure pour estimer un parametre de decalage de frequence d'un oscillateur
JP5051011B2 (ja) サーチ周波数補正方法、測位方法、プログラム及び受信装置
US20060058027A1 (en) Method and apparatus for carrier frequency estimation and correction for GPS
JP4243817B2 (ja) スペクトラム拡散信号復調装置および逆拡散ループの制御方法
Casandra et al. Performance Evaluation of a Tracking Algorithm for Galileo E1 Signals
Bright GPS L2 C signal survey and the development of the emergent MATLAB L2 C (EMAL2) receiver

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A2

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DPEN Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed from 20040101)
WWE Wipo information: entry into national phase

Ref document number: 1020067001732

Country of ref document: KR

WWE Wipo information: entry into national phase

Ref document number: 0335/CHENP/2006

Country of ref document: IN

Ref document number: 335/CHENP/2006

Country of ref document: IN

122 Ep: pct application non-entry in european phase
ENP Entry into the national phase

Ref document number: PI0413053

Country of ref document: BR