WO2005040848A1 - Cross-correlation removal of carrier wave jamming signals - Google Patents

Cross-correlation removal of carrier wave jamming signals Download PDF

Info

Publication number
WO2005040848A1
WO2005040848A1 PCT/US2004/035350 US2004035350W WO2005040848A1 WO 2005040848 A1 WO2005040848 A1 WO 2005040848A1 US 2004035350 W US2004035350 W US 2004035350W WO 2005040848 A1 WO2005040848 A1 WO 2005040848A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
carrier wave
wave jamming
spread spectrum
phase
Prior art date
Application number
PCT/US2004/035350
Other languages
French (fr)
Inventor
Paul A. Underbrink
Henry D. Falk
Charles P. Norman
Original Assignee
Sirf Technology, Inc.
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
Priority claimed from US10/689,565 external-priority patent/US7822105B2/en
Priority claimed from PCT/US2004/028542 external-priority patent/WO2005022187A2/en
Application filed by Sirf Technology, Inc. filed Critical Sirf Technology, Inc.
Publication of WO2005040848A1 publication Critical patent/WO2005040848A1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/21Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
    • G01S19/215Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service issues related to spoofing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/246Acquisition or tracking or demodulation of signals transmitted by the system involving long acquisition integration times, extended snapshots of signals or methods specifically directed towards weak signal acquisition
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/35Constructional details or hardware or software details of the signal processing chain
    • G01S19/37Hardware or software details of the signal processing chain
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0205Details
    • G01S5/0221Receivers
    • G01S5/02213Receivers arranged in a network for determining the position of a transmitter
    • G01S5/02216Timing or synchronisation of the receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details 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/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7097Interference-related aspects
    • H04B1/71Interference-related aspects the interference being narrowband interference
    • H04B1/7101Interference-related aspects the interference being narrowband interference with estimation filters
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04KSECRET COMMUNICATION; JAMMING OF COMMUNICATION
    • H04K3/00Jamming of communication; Counter-measures
    • H04K3/20Countermeasures against jamming
    • H04K3/22Countermeasures against jamming including jamming detection and monitoring
    • H04K3/224Countermeasures against jamming including jamming detection and monitoring with countermeasures at transmission and/or reception of the jammed signal, e.g. stopping operation of transmitter or receiver, nulling or enhancing transmitted power in direction of or at frequency of jammer
    • H04K3/228Elimination in the received signal of jamming or of data corrupted by jamming
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30Acquisition or tracking or demodulation of signals transmitted by the system code related

Definitions

  • This invention relates generally to signal processing in a radio receiver and more particularly to removal of undesired signals received at the radio receiver.
  • RF carrier wave
  • CW jamming is a source of interference in spread spectrum systems, such as CDMA cellular telephone systems and satellite positioning systems.
  • Spread spectrum communication systems use lower power signals spread across the frequency spectrum and are subject to interference from carrier waves used in other communication systems.
  • the problem of CW jamming is further complicated because of the geographical area covered by spread spectrum system may include the whole Earth, for example the United States' Global Position System
  • Systems consistent with the present invention provide a receiver that is capable of receiving a spread spectrum signal that contains a CW jamming signal along with a weak signal.
  • the signal is processed with a crosscorrelator that enables the CW jamming signal to be identified, tracked, and reproduced.
  • the replicated CW jamming signal is subtracted from the received signal after demodulation, thus enabling the weak signal to be processed.
  • FIG. 1 is a flow diagram depicting strong signal cancellation in a weak spread spectrum signal using crosscorrelation.
  • FIG. 2 illustrates a signal-processing diagram for identifying and removing CW jamming signals where the cancelled signal of FIG. 1 maybe the CW jamming signal.
  • FIG. 3 is a block diagram of electrical components of FIG. 2.
  • FIG. 4 is a flow diagram of the identifying and removing CW jamming signals.
  • FIG. 1 a flow diagram 100 depicting strong signal cancellation in a weak spread spectrum signal using crosscorrelation is shown.
  • a strong/weak or near/far signal isolation provided by a spread spectrum, pseudo random number (PRN) code family such as used in CDMA spread spectrum systems is dependent upon the crosscorrelation between the various code members of the family.
  • PRN pseudo random number
  • the isolation of two signals at the same frequency (or multiples of the code repetition rate, in this case 1 KHz) is about 21 to 23 dB. If the relative strengths of two signals differ by more than this limit, the weaker signal cannot be discriminated using only the spreading code.
  • a method of removing the effects of the stronger signal may be applied if the weaker signal is to be tracked.
  • the crosscorrelation effect is at its maximum when the relative Doppler frequency offset between the relatively strong and weak signals is an integer multiple of 1 KHz in the case of coarse acquisition code (C/A) in the GPS signals.
  • C/A coarse acquisition code
  • a general solution to the problem of tracking a weak signal spread spectrum signal in the presence of a stronger spread spectrum signal is based on the premise that all aspects of the strong signal's interference can either be measured or calculated in order to remove it from the weaker signal.
  • the solution can be implemented in any multi-channel receiver having the ability to control a channel's frequency and phase as well as selecting the desired spreading code and setting that code's phase position.
  • the receiver typically employs two channels, one to track the weak signal and one to track the interfering strong signal.
  • the channel that is used to track the strong signal is not required if the characteristics such as power, code phase and frequency of the strong signal can be obtained or accurately estimated by alternate means.
  • the procedure starts 102 with the strong signal being acquired 104, such as by tracking the strong signal in a first channel of the receiver.
  • the channel provides a measurement of the signal strength of the strong signal along with the phase of the ( carrier signal and the spreading code. Additional channels may be used to track additional strong signals (not shown in FIG. 1).
  • the code phase of the spreading code of the weaker signal are predicted 106 based on the 50Hz navigation data code data (D) by methods known in the art.
  • a second channel in the receiver is dedicated to receiving the compound carrier signal and tracking 108 the predicted weak signal component.
  • the second receiver channel correlates the incoming signal with the second code at the predicted frequency and signal phase.
  • the resulting in-phase and quadrature (I,Q) measurements contain both the weak signal and the strong signal, each spread by their unique code.
  • the product Code2R*Code2 is the autocorrelation of the received code 2 and the replica code 2.
  • the autocorrelation function has a value of 1 if the replica code is aligned with the received code. This crosscorrelation of replica code 2 with
  • Strong X (Code2R*CodeX) is next computed 110 to be removed from the compound signal.
  • Codel and Code2 are both members of a PRN code family, and their autocorrelation and crosscorrelation properties are known. It is therefore possible to calculate the crosscorrelation of the two codes at their respective phases by simply multiplying each bit of Codel by the corresponding (in time) bit of Code2 to produce their crosscorrelation value. Since there may be a relative Doppler frequency offset between the two codes, the phase of the codes will process past one another over time and create a new crosscorrelation function. For the GPS system the greatest delta code Doppler typically encountered is about plus or minus 9KHz which is equivalent to six code chips per second (1540 carrier cycles per code chip), and thus the maximum recalculation rate of the crosscorrelation value is roughly 6 times per second.
  • the maximum crosscorrelation occurs at a frequency offset of zero with peaks occurring at intervals of 1000Hz. There is an attenuation of the crosscorrelation as the frequency offset moves away from zero. This attenuation follows the well known sin(x)/x curve. If 10ms measurements are used for tracking or acquisition, the attenuation factor would be equal to sin( ⁇ freq* ⁇ /100 Hz)/( ⁇ freq* ⁇ /100 Hz). This produces an attenuation of -lOdB at about a 75Hz delta frequency. Other local peaks in the sin(x)/x curve (i.e. locally minimum attenuation) occur at 150Hz and 250Hz with attenuations of -13.5dB and -18dB, respectively. This implies that for a desired strong signal suppression of lOdB, only the first lobe of the sin(x)/x function need be considered; however, should additional suppression be desired, the entire curve may be considered.
  • the next step entails computing 112, for each strong signal, the product of the strong signal amplitude and the calculated frequency and time domain (code phase) crosscorrelation.
  • the weak signal is finally extracted by subtracting 114 this product from the compound signal and processing is complete 116.
  • the weak signal thus extracted is subsequently processed in the receiver circuitry as known in the art.
  • the in-phase (I) and quadrature amplitude (Q) of each strong signal is obtained by measurement in each strong signal's own individual receiver channel or by estimation through independent means. Because the strong signal is being actively tracked by the receiver's phase lock loops, the phase of the strong signal is presumed to be near zero radians and thus nearly all the signal power is in the in-phase portion.
  • a signal comprising a strong signal SI modulated with a first code Codel summed with a weak signal w2 modulated with a second code Code2 produces (Sl*codel+w2*Code2).
  • the sum of the two signals is correlated with a replica of the second code Code2R to produce ⁇ Code2R*(Sl*Codel+w2*Code2) ⁇ , where the sum ⁇ includes all chips of the PRN code used to modulate the weak signal w2.
  • the autocorrelation of a code with itself is 1 so the preceding equation can be rewritten as ⁇ Sl*Codel*Code2+w2 ⁇ . It can see that in order to obtain w2 Sl*Codel*Code2 is removed.
  • the code dependent portion of the crosscorrelation factor is computed from the known relative states of the PRN code generators to predict the crosscorrelation between a strong signal of unit power and zero frequency offset, and a weak signal. This factor is multiplied by the amplitude of the corresponding strong signal and adjusted for frequency attenuation before it is subtracted from the composite signal.
  • the various Gold codes used to modulate the PRN signals are all derived from a two code sequence Gl and G2 where the bits of the two code sequences are combined through an
  • Gold code selected. It is known that an XOR operation using binary numbers is mathematically equivalent to multiplication of ⁇ 1. This allows expressing the equations below in term of products of ⁇ 1 while in reality the implementation could be with binary numbers with XORs.
  • I Summation index ranges from 0 to 1022.
  • SatlGl(I) Value of satellite l's Gl coder chip at state I. Possible values are ⁇ 1.
  • SatlG2(I) Value of satellite l's G2 coder chip at state I. Possible values are ⁇ 1.
  • Sat2Gl(I) Value of satellite 2's Gl coder chip at state I. Possible values are ⁇ 1.
  • Sat2G2(I) Value of satellite 2's G2 coder chip at state I. Possible values are ⁇ 1.
  • offset time difference between the satellite 1 and 2 in units of chips ⁇ *Phase change per chip between satellite 1 and 2 in radians.
  • the 1023 states of Gl and G2 may be stored linearly in permanent memory. Thus, it is possible to quickly gather 8, 16, 32 or some other convenient number of bits with a single controller load instruction by computing the address of the desired chip and the shift required to align it. Thirty-two bits is a particularly convenient number because 31 divides 1023 evenly. If a single 1024 bit correlator match filter is employed, the 32 bits of the PN code are derived at a time from the G1/G2 states and loaded into a 1024 wide reference shift register. The current implementation thus reads 32 bits at a time and uses 31 of them at a time for each of 33 intervals that span the 1023 chips of the C/A code. The 31 bit sums are broken into four parts of 8, 8, 8, and 7 bits, and each 7 or 8 bit sum is multiplied by e "j ⁇ I where I changes by 7.75 chips for each part. The form of the sum is:
  • I Outer index ranges from 0 to 32
  • J Inner index ranges from 0 to 7 for the first three sums and from 0 to 6 for the last sum.
  • the inner sums are computed in parallel by using a 32 bit word that contains all 31 bits and using bitwise XOR to perform the multiplications and shifting and adding to sum the 1 bit products.
  • Note that all of the multiplications of the Gl and G2 codes in the above equation are implemented by bit- wise XOR instructions. The above algorithm is in error by at most -17 dB from an exact computation, and requires about 6000 controller operations to complete.
  • the code dependent crosscorrelation factors are computed for all strong and weak signal pairs with small frequency differences, i.e. frequency differences that could cause strong-weak crosscorrelation interference.
  • strong signals are those with C/N 0 >40 dB and weak signals are those with C/N 0 >30 dB.
  • the maximum "significant" frequency difference is 90Hz.
  • the code dependent cross correlation factor for each possibly interfering pair of signals is computed for each of the measurements that might potentially be used by the tracking and signal processing algorithms. For example, if early, punctual and late measurements are used by the tracking loops, the correlation factors for each of these code alignments is computed and stored in the tables.
  • ⁇ F Frequency difference between a strong and weak signal in Hz
  • Mod modulo offset to give a range of -500 Hz to +500 Hz
  • the attenuation only needs to be recomputed if the frequency difference changes by more than 5Hz.
  • phase and amplitude of the strong signal is required to remove the crosscorrelation.
  • the method used in the preferred embodiment is to track the strong signal on its own dedicated channel and collect the I, Q measurements output over the exact same interval that the weak signal I, Q samples are taken.
  • the known phase and frequency of the replica signal that is used to track the strong signal is an excellent approximation of the actual phase and frequency of the strong signal.
  • the magnitude of the I measurement provides a good approximation of the amplitude of the strong signal.
  • the bi-phase modulation of the strong signal data bits D may cause the phase of the strong signal to rotate 180 degrees whenever the data bits transition from a 1 to 0 or from a 0 to 1.
  • the phase of the strong signal is corrected by adding 180 degrees to the phase of the replica signal whenever the sign of the I measurement for the strong signal is negative.
  • SecondPhase SecondCorrelationPhase+DeltaPhase
  • FirstMag FirstCorrelationMag+FirstCodeOffsetFraction* Strongl *FrequencyAttentuation
  • SecondMag SecondCorrelationMag*(l -FirstCodeOffsetFraction) * Strong Frequency Attenuation
  • Co ⁇ rectedWeakIQ WeakIQ-FirstMag*e- jFirstPhase -SecondMag*e- jSecondphase
  • StrongDoppler Doppler of last output to the strong signals channel
  • ⁇ T Difference in time between outputs to the weak and strong channels
  • TableEntryOCodeState Code state difference of the first element of the crosscorrelation table.
  • WeakCarrierPhase Carrier phase angle of last output to the weak signal channel
  • StrongCarrierPhase Carrier phase angle of last output to the strong signal channel
  • DeltaKHz Nearest integer multiple of 1 KHz of the difference between the weak and strong channels Doppler. In units of KHz.
  • FirstCorrelationPhase Phase entry in the crosscorrelation table for the chip indicated by FirstCodeOffset
  • SecondCorrelationPhase Phase entry in the crosscorrelation table for the chip indicated by FirstCodeOffset +1 chip.
  • FirstCorrelationMag Magnitude entry in the crosscorrelation table for the chip indicated by FirstCodeOffset.
  • SecondCorrelationMag Magnitude entry in the crosscorrelation table for the chip indicated by FirstCodeOffset+1 chip.
  • CorrectedWeakIQ IQ correlation corrected for crosscorrelation from the strong signal.
  • CorrectedWeaklQ is computed for the early, on time, and late correlators by shifting the FirstCodeOffset appropriately, such as by half a chip each. These modified correlations are then used normally in the carrier and code tracking software for the weak signal. The algorithm attenuates the crosscorrelation by at least 10 dB without attenuating the weak signal, and is repeated for each strong signal that may be interfering with the weak signal.
  • FIG. 1 has described crosscorrelation to cancel a strong signal having a PRN from a weak signal having different PRN
  • the procedure may be modified to cancel a CW jamming signal and enhance a weak signal having a PRN.
  • FIG. 2 a signal-processing diagram for identifying and removing CW jamming signals where the cancelled signal of FIG. 1 may be a CW jamming signal is shown.
  • the receiver 200 receives a spread spectrum signals at antenna 202.
  • the crosscorrelator 204 is placed in a mode to identify CW jamming signals. In that mode, the crosscorrelator employs all ones for a PRN code.
  • the crosscorrelator 204 may have a signal processor 206 and a match filter 208 that despreads the spread spectrum signal.
  • the resultant signal is the CW jamming signal (the strong signal) and is tracked in the tracker 210.
  • An example of the tracker 210 in the current implementation may be a phase lock loop circuit.
  • the identified CW jamming signal is then used to generate a replica CW jamming signal 212.
  • the replica CW jamming signal may be created using a signal processor 214 and a match filter 216.
  • other approaches to generating a desired signal may be employed. Examples of other approaches include, but are not limited to voltage controlled oscillators, digital signal processing, analog signal processing.
  • the spread spectrum signal received by receiver 200 is then processed by a crosscorrelator 214 with dispreading codes rather than the all ones.
  • the received weak signal in the spread spectrum signal is demodulated by the signal processor in the crosscorrelator 218.
  • the generated CW jamming signal is canceled from the weak signal by a canceller 220.
  • the resultant signal is then tracked by a tracker 222.
  • the CW jamming signal that is being tracked by tracker 210 is further processed to remove the desired weak signal from the tracked CW jamming signal resulting in a clean replica of the unwanted signal on the demodulation process of the desired signal in block 224.
  • FIG. 3 a block diagram of electrical components of FIG. 2 is shown.
  • the spread spectrum signal containing the weak signal and the CW jamming signal are received at the receiver 200 via antenna 202.
  • the spread spectrum signal is demodulated and the CW jammer signal is filtered by the jammer filter 302.
  • the resulting CW jamming signal is tracked by jammer tracker 304.
  • the tracked jamming signal is then replicated by a jammer replica wave generator 306 as a replicated CW jamming signal.
  • the CW signal has the characteristic of a constant phase, therefore the replicated CW jamming signal is scaled and rotated 308 to an appropriate phase to be subtracted from the demodulated weak signal.
  • the weak signal is demodulated by a demodulator 310.
  • the phase and magnitude replicated CW jamming signal is then subtracted from the weak signal from the demodulator 302 by a signal canceller 312.
  • the resulting weak signal is then tracked by a track control circuit 314.
  • the track control circuit 314 outputs the weak signal to the demodulator 310 and two signal combiners 316 and 318.
  • the weak signal is subtracted from the CW jammer signal by the signal combiner 316 and the resulting signal is used to scale and rotate the replicated CW jamming signal.
  • the weak signal is subtracted from the CW jamming signal by signal combiner 318 in order to provide a more accurate CW jamming signal to be replicated by the jammer replica wave generator 306.
  • FIG. 4 a flow diagram 400 of the identifying and removing CW jamming signals is shown.
  • the procedure starts 402 with a receiver 200 receiving a spread spectrum signal 404.
  • the signal may be filtered and the CW jamming signal is identified 406 using a crosscorrelator.
  • the tracked CW jamming signal is then replicated 410 by a jammer replica wave generator 306.
  • the replicated CW jammer signal is then subtracted from the received signal 412.
  • the CW jamming signal has been removed from the received signal and processing is complete 414.

Abstract

A receiver capable of receiving a spread spectrum signal and having a crosscorrelator that enables a carrier wave (CW) jamming to be identified, tracked, replicated and removed from the received spread spectrum signal after demodulation of a weak signal has occurred.

Description

CROSS-CORRELATION REMOVAL OF CARRIER WAVE JAMMING SIGNALS
CROSS REFERENCE TO RELATED APPLICATIONS
[001] This application claims the benefit of U.S. Non-Provisional Patent Application No.
10/689,565, filed on October 20, 2003, and titled "CROSS-CORRELATION REMOVAL OF CARRIER WAVE JAMMING SIGNALS." This application also claims priority to International Patent Application No. PCT/US04/28542, filed on September 2, 2004, and titled "CONTROL AND FEATURES FOR SATELLITE POSITIONING SYSTEM RECEIVERS" and to International Patent Application No. PCT/US04/28926, filed on September 2, 2004, and titled "SIGNAL PROCESSING SYSTEM FOR SATELLITE POSITIONING SIGNALS", which claim priority to United States Provisional Patent Application No. 60/546,816, filed on February 23, 2004, titled "CONTROL AND FEATURES FOR SATELLITE POSITIONING SYSTEM RECEIVERS", and United States Provisional Patent Application No. 60/547,385, filed on February 23, 2004, titled "OVERALL SYSTEM ARCHITECTURE AND RELATED FEATURES", which claim priority to US Provisional Patent Application No. 60/499,961, filed on September 2, 2003 and titled "A GPS SYSTEM" which is hereby incorporated by reference in its entirety.
BACKGROUND OF THE INVENTION [002] 2. Field of the Invention.
[003] This invention relates generally to signal processing in a radio receiver and more particularly to removal of undesired signals received at the radio receiver.
[004] 3. Related Art.
[005] Numerous communication systems rely on receiving one or more radio frequency
(RF) signals. As the available frequency bands become more congested, the numerous radio frequencies used by the different communication systems start to interfere with each other. One type of interference encountered in the different communication systems is carrier wave interference or more commonly called carrier wave (CW) jamming.
[006] CW jamming is a source of interference in spread spectrum systems, such as CDMA cellular telephone systems and satellite positioning systems. Spread spectrum communication systems use lower power signals spread across the frequency spectrum and are subject to interference from carrier waves used in other communication systems. The problem of CW jamming is further complicated because of the geographical area covered by spread spectrum system may include the whole Earth, for example the United States' Global Position System
(GPS).
[007] Attempts to eliminate CW jamming signals have occurred by regulating frequency use and by adding dedicated circuitry to radio receivers. The added circuitry being included in radio receivers often result in additional power consumption and expense. Some examples of jamming signal detection approaches and signal performance improvement approaches may be found in United States Patents 6,448,925, 6,577,271, 6,606,349, and patent application
09/562,749.
[008] Therefore, there is a need for methods and systems for identifying and removing CW jamming signals that overcomes the disadvantages set forth above, and others previously experienced.
SUMMARY [009] Systems consistent with the present invention provide a receiver that is capable of receiving a spread spectrum signal that contains a CW jamming signal along with a weak signal. The signal is processed with a crosscorrelator that enables the CW jamming signal to be identified, tracked, and reproduced. The replicated CW jamming signal is subtracted from the received signal after demodulation, thus enabling the weak signal to be processed. [010] Other systems, methods, features and advantages of the invention will be or will become apparent to one with skill in the art upon examination of the following figures and detailed description. It is intended that all such additional systems, methods, features and advantages be included within this description, be within the scope of the invention, and be protected by the accompanying claims.
BRIEF DESCRIPTION OF THE FIGURES [Oil] The components in the figures are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention. In the figures, like reference numerals designate corresponding parts throughout the different views.
[012] , FIG. 1 is a flow diagram depicting strong signal cancellation in a weak spread spectrum signal using crosscorrelation.
[013] FIG. 2 illustrates a signal-processing diagram for identifying and removing CW jamming signals where the cancelled signal of FIG. 1 maybe the CW jamming signal. [014] FIG. 3 is a block diagram of electrical components of FIG. 2.
[015] FIG. 4 is a flow diagram of the identifying and removing CW jamming signals.
DETAILED DESCRIPTION [016] Turning first to figure 1, a flow diagram 100 depicting strong signal cancellation in a weak spread spectrum signal using crosscorrelation is shown. A strong/weak or near/far signal isolation provided by a spread spectrum, pseudo random number (PRN) code family such as used in CDMA spread spectrum systems is dependent upon the crosscorrelation between the various code members of the family. In the case of a satellite positioning system, such as GPS, the isolation of two signals at the same frequency (or multiples of the code repetition rate, in this case 1 KHz) is about 21 to 23 dB. If the relative strengths of two signals differ by more than this limit, the weaker signal cannot be discriminated using only the spreading code. A method of removing the effects of the stronger signal may be applied if the weaker signal is to be tracked. [017] The crosscorrelation effect is at its maximum when the relative Doppler frequency offset between the relatively strong and weak signals is an integer multiple of 1 KHz in the case of coarse acquisition code (C/A) in the GPS signals. A general solution to the problem of tracking a weak signal spread spectrum signal in the presence of a stronger spread spectrum signal is based on the premise that all aspects of the strong signal's interference can either be measured or calculated in order to remove it from the weaker signal. The solution can be implemented in any multi-channel receiver having the ability to control a channel's frequency and phase as well as selecting the desired spreading code and setting that code's phase position. The receiver typically employs two channels, one to track the weak signal and one to track the interfering strong signal. However, the channel that is used to track the strong signal is not required if the characteristics such as power, code phase and frequency of the strong signal can be obtained or accurately estimated by alternate means. [018] As shown generally in FIG. 1, the procedure starts 102 with the strong signal being acquired 104, such as by tracking the strong signal in a first channel of the receiver. The channel provides a measurement of the signal strength of the strong signal along with the phase of the ( carrier signal and the spreading code. Additional channels may be used to track additional strong signals (not shown in FIG. 1).
[019] The code phase of the spreading code of the weaker signal, along with its received frequency and signal phase, are predicted 106 based on the 50Hz navigation data code data (D) by methods known in the art. A second channel in the receiver is dedicated to receiving the compound carrier signal and tracking 108 the predicted weak signal component.
[020] The second receiver channel correlates the incoming signal with the second code at the predicted frequency and signal phase. The resulting in-phase and quadrature (I,Q) measurements contain both the weak signal and the strong signal, each spread by their unique code. Correlation by multiplication of the replica code for the second signal, Code2R, with the incoming signal yields the product Code2R*(weak2*Code2+StrongX*CodeX+ . . .) where weak2 is the power of weak signal 2, Code2 is the actual code for satellite 2 broadcasting the weak signal 2, StrongX (X=l, 3, 4, . . . ) is the power of strong signal X, and CodeX is the actual code for satellite X contained in the signal. The product Code2R*Code2 is the autocorrelation of the received code 2 and the replica code 2. The autocorrelation function has a value of 1 if the replica code is aligned with the received code. This crosscorrelation of replica code 2 with
Strong X (Code2R*CodeX) is next computed 110 to be removed from the compound signal.
[021] Codel and Code2 are both members of a PRN code family, and their autocorrelation and crosscorrelation properties are known. It is therefore possible to calculate the crosscorrelation of the two codes at their respective phases by simply multiplying each bit of Codel by the corresponding (in time) bit of Code2 to produce their crosscorrelation value. Since there may be a relative Doppler frequency offset between the two codes, the phase of the codes will process past one another over time and create a new crosscorrelation function. For the GPS system the greatest delta code Doppler typically encountered is about plus or minus 9KHz which is equivalent to six code chips per second (1540 carrier cycles per code chip), and thus the maximum recalculation rate of the crosscorrelation value is roughly 6 times per second. [022] The maximum crosscorrelation occurs at a frequency offset of zero with peaks occurring at intervals of 1000Hz. There is an attenuation of the crosscorrelation as the frequency offset moves away from zero. This attenuation follows the well known sin(x)/x curve. If 10ms measurements are used for tracking or acquisition, the attenuation factor would be equal to sin(Δfreq*π/100 Hz)/(Δfreq*π/100 Hz). This produces an attenuation of -lOdB at about a 75Hz delta frequency. Other local peaks in the sin(x)/x curve (i.e. locally minimum attenuation) occur at 150Hz and 250Hz with attenuations of -13.5dB and -18dB, respectively. This implies that for a desired strong signal suppression of lOdB, only the first lobe of the sin(x)/x function need be considered; however, should additional suppression be desired, the entire curve may be considered.
[023] The next step entails computing 112, for each strong signal, the product of the strong signal amplitude and the calculated frequency and time domain (code phase) crosscorrelation. The weak signal is finally extracted by subtracting 114 this product from the compound signal and processing is complete 116. The weak signal thus extracted is subsequently processed in the receiver circuitry as known in the art.
[024] The in-phase (I) and quadrature amplitude (Q) of each strong signal is obtained by measurement in each strong signal's own individual receiver channel or by estimation through independent means. Because the strong signal is being actively tracked by the receiver's phase lock loops, the phase of the strong signal is presumed to be near zero radians and thus nearly all the signal power is in the in-phase portion.
[025] A signal comprising a strong signal SI modulated with a first code Codel summed with a weak signal w2 modulated with a second code Code2 produces (Sl*codel+w2*Code2). The sum of the two signals is correlated with a replica of the second code Code2R to produce Σ{Code2R*(Sl*Codel+w2*Code2)}, where the sum Σ includes all chips of the PRN code used to modulate the weak signal w2. The autocorrelation of a code with itself is 1 so the preceding equation can be rewritten as Σ{Sl*Codel*Code2+w2}. It can see that in order to obtain w2 Sl*Codel*Code2 is removed. Since we know Codel and Code2, we can easily calculate their crosscorrelation. This leaves us to estimate the value of SI that can be done by independently tracking the strong signal on a separate channel, or by any other convenient means. This computed value of Sl*Codel *Code2 would be sufficient if the strong signal SI and the weak signal w2 where at the same frequency. The two signals are received at different frequencies, however, due to the Doppler effect as well as other factors enumerated previously. [026] We know that strength of the crosscorrelation varies with the difference between these frequencies in a sin(x)/x relationship. We may therefore calculate an attenuation factor based on the difference in frequency between the strong and the weak signal and apply it to the computed crosscorrelation. Furthermore, if more than one strong signal is present, an attenuation factor is computed for each strong signal.
[027] The code dependent portion of the crosscorrelation factor is computed from the known relative states of the PRN code generators to predict the crosscorrelation between a strong signal of unit power and zero frequency offset, and a weak signal. This factor is multiplied by the amplitude of the corresponding strong signal and adjusted for frequency attenuation before it is subtracted from the composite signal.
[028] The various Gold codes used to modulate the PRN signals are all derived from a two code sequence Gl and G2 where the bits of the two code sequences are combined through an
XOR operation after G2 has been offset some number of bits relative to Gl depending on the
Gold code selected. It is known that an XOR operation using binary numbers is mathematically equivalent to multiplication of ±1. This allows expressing the equations below in term of products of ±1 while in reality the implementation could be with binary numbers with XORs.
[029] The correlation between two C/A codes can generally be expressed as:
[030] ΣSatlGl(I)*SatlG2(I)*Sat2Gl(I-offset)*Sat2G2(I-offset)*e"jΔΘI
[031] where
[032] I=Summation index ranges from 0 to 1022.
[033] SatlGl(I)= Value of satellite l's Gl coder chip at state I. Possible values are ±1.
[034] SatlG2(I)= Value of satellite l's G2 coder chip at state I. Possible values are ±1.
[035] Sat2Gl(I)= Value of satellite 2's Gl coder chip at state I. Possible values are ±1.
[036] Sat2G2(I)= Value of satellite 2's G2 coder chip at state I. Possible values are ±1.
[037] offset=time difference between the satellite 1 and 2 in units of chips ΔΘ*Phase change per chip between satellite 1 and 2 in radians.
[038] It should be noted that when the difference I-offset is less than 0, 1023 is added to the difference to maintain the value in the range of 0 to 1022. In other words, the domain of the functions returning coder chip states is limited to the range of 0 to 1022.
[039] The computation time required to compute the 1023 bit-by-bit correlations can be accelerated by making use of standard controller instructions that perform 8, 16 or 32 bit- wise XORs with a single controller instruction. The following will demonstrate the method of computing eight chips in parallel. Those skilled in the art will immediately recognize that the scheme can be easily modified to accommodate some other convenient number of bits per controller XOR operation. Such an approach may employ a single 1024 bit correlator match filter that performs the correlation for 1023 G1/G2 states concurrently.
[040] The 1023 states of Gl and G2 may be stored linearly in permanent memory. Thus, it is possible to quickly gather 8, 16, 32 or some other convenient number of bits with a single controller load instruction by computing the address of the desired chip and the shift required to align it. Thirty-two bits is a particularly convenient number because 31 divides 1023 evenly. If a single 1024 bit correlator match filter is employed, the 32 bits of the PN code are derived at a time from the G1/G2 states and loaded into a 1024 wide reference shift register. The current implementation thus reads 32 bits at a time and uses 31 of them at a time for each of 33 intervals that span the 1023 chips of the C/A code. The 31 bit sums are broken into four parts of 8, 8, 8, and 7 bits, and each 7 or 8 bit sum is multiplied by e"jΔΘI where I changes by 7.75 chips for each part. The form of the sum is:
[041] Σ(eiΔΘI*31*Σ(SatlGl(I*31+J)*SatlG2(I*31+J)*Sat2Gl(I*31+Joffset)*Sat2G2(P31*
J-offset))+e-jΔΘ(P31+7-75)^
*Sat2G2(I*31+J+8-offset))+e"jΔΘ(I*31+15-5)*Σ(SatlGl(I*31+J+16)*SatlG2(I*31+J+16)
*Sat2Gl(I*31+J+16-offset)*Sat2G2(I*31+J+16-offset))+e-jΔΘ(P31+23-25)*Σ(SatlGl(P31+J+24)
*Satl G2(I*3 l+J+24)*Sat2Gl (1*3 l+J+24-offset)*Sat2G2(I*3 l+J+24-offset)))
[042] where
[043] I=Outer index ranges from 0 to 32 [044] J=Inner index ranges from 0 to 7 for the first three sums and from 0 to 6 for the last sum. The inner sums are computed in parallel by using a 32 bit word that contains all 31 bits and using bitwise XOR to perform the multiplications and shifting and adding to sum the 1 bit products. [045] Note that all of the multiplications of the Gl and G2 codes in the above equation are implemented by bit- wise XOR instructions. The above algorithm is in error by at most -17 dB from an exact computation, and requires about 6000 controller operations to complete. [046] Periodically, as needed, the code dependent crosscorrelation factors are computed for all strong and weak signal pairs with small frequency differences, i.e. frequency differences that could cause strong-weak crosscorrelation interference. In the current implementation strong signals are those with C/N0 >40 dB and weak signals are those with C/N0 >30 dB. Because 10ms integrations of I and Q measurements are used by the code and phase tracking loops, the maximum "significant" frequency difference (modulo lOOOHz) is 90Hz. In the preferred embodiment the code dependent cross correlation factor for each possibly interfering pair of signals is computed for each of the measurements that might potentially be used by the tracking and signal processing algorithms. For example, if early, punctual and late measurements are used by the tracking loops, the correlation factors for each of these code alignments is computed and stored in the tables.
[047] These tables only need be updated at a 10Hz rate because the maximum Doppler difference is less than 9KHz or less than six chips per second. In addition to maintaining the crosscorrelation table, the frequency attenuation of the crosscorrelation due to the frequency difference is computed at the 10Hz rate. The attenuation can be expressed as: [048] Frequency Attenuation=sin (ΔF mod 1000*π/l 00)/(ΔF mod 1000*π/l 00) [049] where
[050] ΔF=Frequency difference between a strong and weak signal in Hz Mod=modulo offset to give a range of -500 Hz to +500 Hz [051] The attenuation only needs to be recomputed if the frequency difference changes by more than 5Hz.
[052] An estimate of the phase and amplitude of the strong signal is required to remove the crosscorrelation. The method used in the preferred embodiment is to track the strong signal on its own dedicated channel and collect the I, Q measurements output over the exact same interval that the weak signal I, Q samples are taken. The known phase and frequency of the replica signal that is used to track the strong signal is an excellent approximation of the actual phase and frequency of the strong signal. Furthermore, because the strong signal is in phase lock, the magnitude of the I measurement provides a good approximation of the amplitude of the strong signal. Finally, the bi-phase modulation of the strong signal data bits D may cause the phase of the strong signal to rotate 180 degrees whenever the data bits transition from a 1 to 0 or from a 0 to 1. In the current implementation, the phase of the strong signal is corrected by adding 180 degrees to the phase of the replica signal whenever the sign of the I measurement for the strong signal is negative.
[053] Every 10ms a new set of I, Q correlation data is available from the channel assigned to track the weak signal. The tables of crosscorrelation factors are checked to predict the presence of any interfering strong signals. If strong signals are predicted, the following subtraction is performed to remove the strong signal crosscorrelations: [054] FirstCodeOffset=WeakCodeState-StrongCodeState-StrongDoppler*ΔT- TableEntryOCodeState [055] DeltaPhase=WeakCarrierPhase-StrongCarrierPhase-StrongDoppler* ΔT +DeltaKHz*StrongCodeState [056] FirstPhase=FirstCorrelationPhase+DeltaPhase
[057] SecondPhase=SecondCorrelationPhase+DeltaPhase
[058] FirstMag=FirstCorrelationMag+FirstCodeOffsetFraction* Strongl *FrequencyAttentuation [059] SecondMag=SecondCorrelationMag*(l -FirstCodeOffsetFraction) * Strong Frequency Attenuation [060] CoιrectedWeakIQ=WeakIQ-FirstMag*e-jFirstPhase-SecondMag*e-jSecondphase
[061] where {
[062] WeakCodeState=Code state of last output to the weak signals channel
[063] StrongCodeState=Code state of last output to the strong signals channel
[064] StrongDoppler=Doppler of last output to the strong signals channel ΔT=Difference in time between outputs to the weak and strong channels [065] TableEntryOCodeState=Code state difference of the first element of the crosscorrelation table. [066] WeakCarrierPhase=Carrier phase angle of last output to the weak signal channel
[067] StrongCarrierPhase=Carrier phase angle of last output to the strong signal channel [068] DeltaKHz=Nearest integer multiple of 1 KHz of the difference between the weak and strong channels Doppler. In units of KHz. [069] FirstCorrelationPhase=Phase entry in the crosscorrelation table for the chip indicated by FirstCodeOffset
[070] SecondCorrelationPhase=Phase entry in the crosscorrelation table for the chip indicated by FirstCodeOffset +1 chip.
[071] FirstCorrelationMag=Magnitude entry in the crosscorrelation table for the chip indicated by FirstCodeOffset.
[072] SecondCorrelationMag=Magnitude entry in the crosscorrelation table for the chip indicated by FirstCodeOffset+1 chip.
[073] FirstCodeOffsetFraction^Fraction of a chip in FirstCodeOffset.
[074] StrongI= Absolute value of I correlation from the strong channel.
[075] FrequencyAttenuation=Attenuation due to frequency offset.
[076] WeakIQ=IQ correlation from the weak signal's channel.
[077] CorrectedWeakIQ=IQ correlation corrected for crosscorrelation from the strong signal. CorrectedWeaklQ is computed for the early, on time, and late correlators by shifting the FirstCodeOffset appropriately, such as by half a chip each. These modified correlations are then used normally in the carrier and code tracking software for the weak signal. The algorithm attenuates the crosscorrelation by at least 10 dB without attenuating the weak signal, and is repeated for each strong signal that may be interfering with the weak signal.
[078] While FIG. 1 has described crosscorrelation to cancel a strong signal having a PRN from a weak signal having different PRN, the procedure may be modified to cancel a CW jamming signal and enhance a weak signal having a PRN. In FIG. 2, a signal-processing diagram for identifying and removing CW jamming signals where the cancelled signal of FIG. 1 may be a CW jamming signal is shown. The receiver 200 receives a spread spectrum signals at antenna 202. The crosscorrelator 204 is placed in a mode to identify CW jamming signals. In that mode, the crosscorrelator employs all ones for a PRN code. The crosscorrelator 204 may have a signal processor 206 and a match filter 208 that despreads the spread spectrum signal. The resultant signal is the CW jamming signal (the strong signal) and is tracked in the tracker 210. An example of the tracker 210 in the current implementation may be a phase lock loop circuit.
[079] The identified CW jamming signal is then used to generate a replica CW jamming signal 212. The replica CW jamming signal may be created using a signal processor 214 and a match filter 216. In other implementations, other approaches to generating a desired signal may be employed. Examples of other approaches include, but are not limited to voltage controlled oscillators, digital signal processing, analog signal processing.
[080] The spread spectrum signal received by receiver 200 is then processed by a crosscorrelator 214 with dispreading codes rather than the all ones. The received weak signal in the spread spectrum signal is demodulated by the signal processor in the crosscorrelator 218. The generated CW jamming signal is canceled from the weak signal by a canceller 220. The resultant signal is then tracked by a tracker 222. The CW jamming signal that is being tracked by tracker 210 is further processed to remove the desired weak signal from the tracked CW jamming signal resulting in a clean replica of the unwanted signal on the demodulation process of the desired signal in block 224.
[081] The identification of the CW jamming signal has been described as occurring with three crosscorrelators 204, 218 and 212. But in practice, only one crosscorrelator may be used allowing the crosscorrlator to distinguish between PRN codes and to identify and remove CW jamming signals. Thus, the generation of a replica of the unwanted signal (CW jamming signal) occurs during the demodulation of the wanted signal.
[082] Turning to FIG. 3, a block diagram of electrical components of FIG. 2 is shown. The spread spectrum signal containing the weak signal and the CW jamming signal are received at the receiver 200 via antenna 202. The spread spectrum signal is demodulated and the CW jammer signal is filtered by the jammer filter 302. The resulting CW jamming signal is tracked by jammer tracker 304. The tracked jamming signal is then replicated by a jammer replica wave generator 306 as a replicated CW jamming signal. The CW signal has the characteristic of a constant phase, therefore the replicated CW jamming signal is scaled and rotated 308 to an appropriate phase to be subtracted from the demodulated weak signal. The weak signal is demodulated by a demodulator 310. The phase and magnitude replicated CW jamming signal is then subtracted from the weak signal from the demodulator 302 by a signal canceller 312. The resulting weak signal is then tracked by a track control circuit 314. The track control circuit 314 outputs the weak signal to the demodulator 310 and two signal combiners 316 and 318. The weak signal is subtracted from the CW jammer signal by the signal combiner 316 and the resulting signal is used to scale and rotate the replicated CW jamming signal. Similarly, the weak signal is subtracted from the CW jamming signal by signal combiner 318 in order to provide a more accurate CW jamming signal to be replicated by the jammer replica wave generator 306.
[083] Even though the current implementation has been described with CW jamming signals in general, one skilled in the art would recognize that specific known signals such as IS- 95 pilot signals, Galileo signals, or other known carriers may be identified, tracked and remove from a demodulated desired signal. Unlike other known methods of identifying and removing CW jamming signals prior to signal processing or at the front end of a receiver, the current approach identifies and removes the CW jamming signal at the back end. This is advantageous when desired signals are received at low signal strength, such as in CDMA or satellite positioning system. The signal strength is not sufficient to support preprocessing of the signal while maintaining the desired weak signals.
[084] Turning to FIG. 4, a flow diagram 400 of the identifying and removing CW jamming signals is shown. The procedure starts 402 with a receiver 200 receiving a spread spectrum signal 404. The signal may be filtered and the CW jamming signal is identified 406 using a crosscorrelator. Once the CW jamming signal is identified 406, it is tracked 408. The tracked CW jamming signal is then replicated 410 by a jammer replica wave generator 306. The replicated CW jammer signal is then subtracted from the received signal 412. The CW jamming signal has been removed from the received signal and processing is complete 414. [085] The foregoing description of an implementation has been presented for purposes of illustration and description. It is not exhaustive and does not limit the claimed inventions to the precise form disclosed. Modifications and variations are possible in light of the above description or may be acquired from practicing the invention. For example, the described implementation includes software but the invention may be implemented as a combination of hardware and software or in hardware alone. Note also that the implementation may vary between systems. The claims and their equivalents define the scope of the invention.

Claims

CLAIMS What is claimed is:
1. A radio receiver apparatus in receipt of a spread spectrum radio signal having a first signal tracking channel and a second signal tracking channel, comprising: a demodulator that demodulates a first signal in the spread spectrum radio signal into the first signal tracking channel and a second signal in the spread spectrum radio signal into the second signal tracking channel; a crosscorrelator connected to the first tracking channel and the second tracking channel; a signal processor that identify a carrier wave jamming signal with the crosscorrelator; a tracker that tracks the carrier wave jamming signal; and a signal canceller subtracts the carrier wave jamming signal from the spread spectrum signal.
2. The radio receiver apparatus of claim 1 , where the signal canceller further including: a signal generator that generates a replica carrier wave jamming signal having a phase from the carrier wave jamming signal having another phase and subtracts the replica carrier wave jamming signal from the spread spectrum signal to cancel the carrier wave jamming signal.
3. The radio receiver apparatus of claim 2, further includes: a signal rotator that rotates the phase of the replica carrier wave jamming signal.
4. The radio receiver apparatus of claim 3, where the signal rotator adjusts the phase of the replica carrier wave jamming signal to match the other phase of the carrier wave jamming signal in the spread spectrum signal.
5. The radio receiver apparatus of claim 1, where the crosscorrelator has a code of all ones for a pseudo random number (PRN) code.
6. The radio receiver apparatus of claim 1, wherein the spread spectrum radio signal is a position signal.
7. The radio receiver apparatus of claim 1, wherein the crosscorrelator is at least a 1024 bit wide correlator.
8. The radio receiver apparatus of claim 7, where the crosscorrelator further includes: an at least a 1024 bit wide match filter.
9. A method of removing a carrier wave jamming signal from a spread spectrum signal having a first signal tracking channel and a second signal tracking channel, comprising: demodulating a first signal in the spread spectrum radio signal into the first signal tracking channel and a second signal in the spread spectrum radio signal into the second signal tracking channel; correlating the first tracking channel and the second tracking channel with a crosscorrelator; computating a product of the first signal fracking channel and the second signal tracking channel to obtain a carrier wave jamming signal; tracking the carrier wave jamming signal; and canceling the carrier wave jamming signal from the spread spectrum signal.
10. The method of claim 9, where canceling further includes: generating a replica carrier wave jamming signal having a phase from the carrier wave jamming signal having another phase; and subtracting the replica carrier wave jamming signal from the spread spectrum signal to cancel the carrier wave jamming signal.
11. The method of claim 10, further includes: rotating the phase of the replica carrier wave jamming signal.
12. The method of claim 11 , further includes : adjusting the phase of the replica carrier wave jamming signal to match the other phase of the carrier wave jamming signal in the spread spectrum signal.
13. The method of claim 9, further includes: configuring the crosscorrelator with a code of all ones for a pseudo random number (PRN) code.
14. The method of claim 9, wherein the spread spectrum radio signal is a position signal.
15. The method of claim 9, wherein the crosscorrelator is at least a 1024 bit wide correlator.
16. The method of claim 15, where the crosscorrelator further includes: filtering with an at least a 1024 bit wide match filter.
17. A receiver in receipt of a spread spectrum radio signal having a first signal fracking channel and a second signal tracking channel, comprising: demodulation means for demodulating a first signal in the spread spectrum radio signal into the first signal fracking channel and a second signal in the spread spectrum radio signal into the second signal tracking channel; correlation means for correlating the first tracking channel and the second tracking channel; computation means for computing a product of the first signal tracking channel and the second signal tracking channel to obtain a carrier wave jamming signal; means for tracking the carrier wave jamming signal; and canceling means that cancels the carrier wave jamming signal from the spread spectrum signal.
18. The receiver of claim 17, where cancellation means further includes: generating means for generation of a replica carrier wave jamming signal having a phase from the carrier wave jamming signal having another phase; and subtracting means for subtraction of the replica carrier wave jamming signal from the spread spectrum signal to cancel the carrier wave jamming signal.
19. The receiver of claim 18 , further includes : means for rotating the phase of the replica carrier wave jamming signal.
20. The receiver of claim 19, further includes : means for adjusting the phase of the replica carrier wave jamming signal to match the other phase of the carrier wave jamming signal in the spread spectrum signal.
21. The receiver of claim 17, further includes : means for configuring the correlation means with a code of all ones for a pseudo random number (PRN) code.
22. The radio receiver apparatus of claim 17, wherein the spread spectrum radio signal is a positioning signal.
23. The radio receiver apparatus of claim 9, wherein the correlation means is at least a 1024 bit wide correlator.
24. The radio receiver apparatus of claim 23, where the correlation means further includes: an at least a 1024 bit wide match filter.
PCT/US2004/035350 2003-10-20 2004-10-20 Cross-correlation removal of carrier wave jamming signals WO2005040848A1 (en)

Applications Claiming Priority (10)

Application Number Priority Date Filing Date Title
US10/689,565 US7822105B2 (en) 2003-09-02 2003-10-20 Cross-correlation removal of carrier wave jamming signals
US10/689,565 2003-10-20
US54681604P 2004-02-23 2004-02-23
US54738504P 2004-02-23 2004-02-23
US60/546,816 2004-02-23
US60/547,385 2004-02-23
PCT/US2004/028542 WO2005022187A2 (en) 2003-09-02 2004-09-02 Control and features for satellite positioning system receivers
USPCT/US04/28542 2004-09-02
PCT/US2004/028926 WO2005047923A2 (en) 2003-09-02 2004-09-02 Signal processing system for satellite positioning signals
USPCT/US04/28926 2004-09-02

Publications (1)

Publication Number Publication Date
WO2005040848A1 true WO2005040848A1 (en) 2005-05-06

Family

ID=34528385

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2004/035350 WO2005040848A1 (en) 2003-10-20 2004-10-20 Cross-correlation removal of carrier wave jamming signals

Country Status (1)

Country Link
WO (1) WO2005040848A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009158049A1 (en) * 2008-06-27 2009-12-30 Sirf Technology, Inc. Method and apparatus for mitigating the effects of cw interference via post correlation processing in a gps receiver
CN103954977A (en) * 2014-05-12 2014-07-30 武汉大学 Method and system for sensing GNSS deception jamming
EP2930536A3 (en) * 2014-04-09 2015-10-28 The Mitre Corporation Positioning, navigation, and timing device interference and spoofing detector with timing mitigation
CN112180407A (en) * 2020-09-14 2021-01-05 北京宏德信智源信息技术有限公司 Method for eliminating mutual interference of strong and weak signals

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0884855A1 (en) * 1996-11-29 1998-12-16 Mitsubishi Denki Kabushiki Kaisha Shared wave removing device
WO2000055992A1 (en) * 1999-03-17 2000-09-21 Ericsson, Inc. Synchronization and cell search method and apparatus for wireless communications
US6282231B1 (en) * 1999-12-14 2001-08-28 Sirf Technology, Inc. Strong signal cancellation to enhance processing of weak spread spectrum signal

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0884855A1 (en) * 1996-11-29 1998-12-16 Mitsubishi Denki Kabushiki Kaisha Shared wave removing device
WO2000055992A1 (en) * 1999-03-17 2000-09-21 Ericsson, Inc. Synchronization and cell search method and apparatus for wireless communications
US6282231B1 (en) * 1999-12-14 2001-08-28 Sirf Technology, Inc. Strong signal cancellation to enhance processing of weak spread spectrum signal

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009158049A1 (en) * 2008-06-27 2009-12-30 Sirf Technology, Inc. Method and apparatus for mitigating the effects of cw interference via post correlation processing in a gps receiver
GB2473792A (en) * 2008-06-27 2011-03-23 Sirf Tech Inc Method and apparatus for mitigating the effects of cw interference via post correlation processing in a GPS receiver
GB2473792B (en) * 2008-06-27 2012-11-07 Sirf Tech Inc Method and apparatus for mitigating the effects of cw interference via post correlation processing in a GPS receiver
EP2930536A3 (en) * 2014-04-09 2015-10-28 The Mitre Corporation Positioning, navigation, and timing device interference and spoofing detector with timing mitigation
US10509130B2 (en) 2014-04-09 2019-12-17 The Mitre Corporation Positioning, navigation, and timing device interference and spoofing detector with timing mitigation
CN103954977A (en) * 2014-05-12 2014-07-30 武汉大学 Method and system for sensing GNSS deception jamming
CN112180407A (en) * 2020-09-14 2021-01-05 北京宏德信智源信息技术有限公司 Method for eliminating mutual interference of strong and weak signals

Similar Documents

Publication Publication Date Title
US7822105B2 (en) Cross-correlation removal of carrier wave jamming signals
EP1238485B1 (en) Strong signal cancellation to enhance processing of weak spread spectrum signal
EP1832005B1 (en) Method for multi-path detection in spread spectrum receivers
US5596600A (en) Standalone canceller of narrow band interference for spread spectrum receivers
RU2339051C2 (en) SATELLITE "GALILEO" SIGNAL AltBOC PROCESSING EQUIPMENT
US5268927A (en) Digital adaptive transversal filter for spread spectrum receivers
US6711219B2 (en) Interference cancellation in a signal
KR940002993B1 (en) Null processing receiver apparatus and method
US6879913B1 (en) Indoor GPS clock
US5949815A (en) Code multipath reduction using weighted correlation channels
US20110254732A1 (en) System and method for multi-correlation with modulation-adapted filter for the fast acquisition and the tracking of a radio navigation signal in the presence of jamming
KR100810481B1 (en) Methods and apparatuses for processing of global positioning system using a matched filter
EP1321774B1 (en) A method for performing reacquisition in a positioning receiver, and corresponding electronic device
Glennon et al. A review of GPS cross correlation mitigation techniques
Mongrédien et al. Testing GPS L5 acquisition and tracking algorithms using a hardware simulator
JP2004515962A (en) Method for splitting a signal and signal processing circuit and device utilizing such a method
US6327298B1 (en) Post-correlation temporal nulling
EP2188907B1 (en) Method and device for multipath mitigation
WO2005040848A1 (en) Cross-correlation removal of carrier wave jamming signals
US20070076786A1 (en) Methods and apparatuses for processing spread spectrum signals
Shanmugam Improving GPS L1 C/A code correlation properties using a novel multi-correlator differential detection technique
TWI712270B (en) System and method for mitigating interference in the presence of interference
Mao et al. Adaptive multipath mitigation tracking system for GPS receiver
Iswariya et al. FFT based acquisition techniques of GPS L2C signals
Bose GPS Satellite Signal Acquisition and Tracking

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

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: A1

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
122 Ep: pct application non-entry in european phase