Spectral estimation of radar time-of-arrival periodicities

• US 5,396,250 A
• Filed: 01/27/1994
• Issued: 03/07/1995
• Est. Priority Date: 12/03/1992
• Status: Expired due to Fees

• Current Assignee: United States Of America As Represented By The Secretary Of The Air Force
• Original Assignee: United States Of America As Represented By The Secretary Of The Air Force
• Inventors: Sharpin, David L., Shaw, Rudy L., Tsui, James B. Y.
• Primary Examiner(s): Sotomayor, John B.
• Application Number:
  US08/187,356
• Time in Patent Office: 404 Days
• Field of Search: 342/13, 342/196
• US Class Current: 342/13
• CPC Class Codes: G01S 7/021 Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
• CPC Class Codes: G01S 7/2923 based on data belonging to a number of consecutive radar periods

Abstract

Disclosed is a method to extract the periodic properties (or Pulse Repetition Intervals (PRIs)) of radar signals whose time-of-arrival at an airborne platform have been time tagged by an Electronic Warfare (EW) receiver. The PRIs are determined using a modified Discrete Fourier Transform (DFT) written as ##EQU1## where k represents the frequency components, the T_i are the individual Time-Of-Arrival (TOA) data values and N is the last TOA value collected.

Disclosed are three possible methods of reducing the number of computations involved with the modification of the Discrete Fourier Transform (DFT) equation to facilitate its use with Radar Time-Of-Arrival (TOA) to extract periodic properties of radar signals. This is particularly applicable to radar pulse trains which are interleaved in time and where each individual pulse train may be staggered or jittered in time so that a time domain deinterleaving scheme must be employed. A spectral estimator can be computed by taking the modulus of the DFT of the time series data and extending the periodogram estimator by making a slight modification of the DFT equation.

Claims

1. Apparatus to extract the Pulse Repetition Intervals (PRIs) of radar signals comprising a train of radar pulses, wherein said apparatus comprises:

an Electronic Warfare (EW) receiver at an airborne platform which includes means for receiving said radar signals and time tagging the time-of-arrival (TOA) of each radar pulse as a time data value.

a memory in which is stored M values of a quantity _R(k)=e -j2π

k/N for values of k from 0 to (N-1) inclusive;

where k represents the frequency components, the T_i are the individual time-of-arrival (TOA) data values and N is the last TOA value collected; and

math means having inputs from the memory supplying said M values and input of the time date values from said Electronic Warfare receiver effective for each time data value T_i for performing the operation ##EQU7## on each of the M values R(k) from the memory and supplying the M results to M summing means respectively;

wherein the M summing means include means to accumulate the values of R(k) for all values of T_i in the train of radar pulses, and to supply the resulting M values to further processing means;

wherein the further processing means includes to determine the Pulse Repetition Intervals using a Discrete Fourier Transform (DFT) equation written as ##EQU8##

2. Apparatus according to claim 1, including means wherein, approximate values of pulse repetition frequencies being known a priori, only certain k values are computed and searched.

3. Apparatus according to claim 1, including means for first performing said Discrete Fourier Transform (DFT) equation on a coarse set of TOA values and finding dominant sets of k values, and means for then performing said Discrete Fourier Transform (DFT) equation on fine resolution TOA values only around the dominant k values.

4. Apparatus according to claim 1, including means for calculating said Discrete Fourier Transform (DFT) equation starting from a minimum value kmin, through some relatively low k value (k₁), means for finding the lowest k value in this range, then using this value to remove TOA pulses that correspond to the peak k value'"'"'s PRF;

said Discrete Fourier Transform (DFT) equation being again calculated, starting from the previous k₁ and the process of removal being repeated until the last TOA pulse train is removed.

5. A method to extract the Pulse Repetition Intervals (PRIs) of radar signals comprising a train of radar pulses, wherein said method comprises:

using an Electronic Warfare (EW) receiver at an airborne platform for receiving said radar signals and time tagging the time-of-arrival (TOA) of each radar pulse as a time data value;

storing in a memory M values of a quantity _R(k)=e -j2π

k/N for values of k from 0 to (N-1) inclusive;

where k represents the frequency components, the T_i are the individual time-of-arrival (TOA) data values and N is the last TOA value collected;

using inputs from the memory supplying said M values and input of the time data values from said Electronic Warfare receiver effective for each time data value T_i and performing the operation ##EQU9## on each of the M values R(k) from the memory and supplying the M results to M summing means respectively; and

using the M summing means to accumulate the values of R(k) for all values of T_i in the train of radar pulses, and supplying the resulting M values to further processing means;

wherein in the further processing means the Pulse Repetition Intervals are determined using a Discrete Fourier Transform (DFT) equation written as ##EQU10##

6. A method according to claim 5, wherein approximate values of pulse repetition frequencies being known a priori, only certain k values are computed and searched.

7. A method according to claim 5, comprising first performing said Discrete Fourier Transform (DFT) equation on a coarse set of TOA values and finding dominant sets of k values, and then performing said Discrete Fourier Transform (DFT) equation on fine resolution TOA values only around the dominant k values.

8. A method according to claim 5, comprising calculating said Discrete Fourier Transform (DFT) equation starting from a minimum value k_min, through some relatively low k value (k₁), finding the lowest k value in this range, then using this value to remove TOA pulses that correspond to the peak k value'"'"'s PRF;

said Discrete Fourier Transform (DFT) equation being again calculated, starting from the previous k₁ and the process of removal being repeated until the last TOA pulse train is removed.