Search Images Maps Play YouTube News Gmail Drive More »
Sign in
Screen reader users: click this link for accessible mode. Accessible mode has the same essential features but works better with your reader.

Patents

  1. Advanced Patent Search
Publication numberUS20030156057 A1
Publication typeApplication
Application numberUS 10/303,057
Publication dateAug 21, 2003
Filing dateNov 25, 2002
Priority dateNov 24, 2001
Also published asDE60108403D1, DE60108403T2, EP1314997A1, EP1314997B1, US6686871
Publication number10303057, 303057, US 2003/0156057 A1, US 2003/156057 A1, US 20030156057 A1, US 20030156057A1, US 2003156057 A1, US 2003156057A1, US-A1-20030156057, US-A1-2003156057, US2003/0156057A1, US2003/156057A1, US20030156057 A1, US20030156057A1, US2003156057 A1, US2003156057A1
InventorsHermann Rohling, Wilhelm Guener, Hans Hommel, Dieter Nagel
Original AssigneeHermann Rohling, Wilhelm Guener, Hans Hommel, Dieter Nagel
Export CitationBiBTeX, EndNote, RefMan
External Links: USPTO, USPTO Assignment, Espacenet
Method for HPRF-radar measurement
US 20030156057 A1
Abstract
In a method for HPRF-radar measurement of the range and Doppler frequency of at least one target, a transmit signal is generated which consists of two pulse sequences that are interleaved on a pulse to pulse basis, and have the same pulse repetition frequency PRF and the same transmit frequency. The pulses of a first one of the two pulse sequences have a linearly increasing phase value with a fixed phase difference φ1n (greater than zero) from pulse to pulse with φ1n≧0, while the pulses of the second pulse sequence have a linear increasing phase value with a fixed phase difference φ2n which differs from φ1n. The two received base band signals of each individual pulse sequence are Fourier transformed, and the amplitude peaks of the resulting two Fourier spectra are determined. The amplitude peaks of the two Fourier spectra are then clustered to determine pairs of spectral lines of a Doppler frequency shift of Δ f = ϕ 2 n - ϕ 1 n 2 π PRF .
The target Doppler frequency is calculated based on two Doppler frequencies of a previously determined spectral line pair, and the target range is calculated based on the phase difference between the two spectral lines of a previously determined spectral line pair.
Images(4)
Previous page
Next page
Claims(2)
What is claimed is:
1. A method for HPRF-radar measurement of range and Doppler frequency for at least one target, said method comprising:
generating a transmit signal which includes first and second pulse sequences that are interleaved on a pulse to pulse basis, and have the same pulse repetition frequency PRF and the same transmit frequency, wherein
pulses of the first pulse sequence have a linearly increasing phase value, with a fixed phase difference φ1n from pulse to pulse, with φ1n≧0;
pulses of the second pulse sequence have a linear increasing phase value with a fixed phase difference φ2n from pulse to pulse, with φ1n≠φ2n and φ1n≧0;
Fourier transforming received base band signals of each individual pulse sequence;
determining amplitude peaks of resulting two Fourier spectra;
clustering the amplitude peaks of the two Fourier spectra to determine pairs of spectral lines of a Doppler frequency shift of
Δ f = ϕ 2 n - ϕ 1 n 2 π . PRF ;
calculating target Doppler frequency based on the two Doppler frequencies of a previously determined spectral line pair; and
calculating target range based on a phase difference between spectral lines of a previously determined spectral line pair.
2. An apparatus for HPRF-radar measurement of range and Doppler frequency for at least one target, comprising:
means for generating a transmit signal which includes first and second pulse sequences that are interleaved on a pulse to pulse basis, and have the same pulse repetition frequency PRF and the same transmit frequency, wherein
pulses of the first pulse sequence have a linearly increasing phase value, with a fixed phase difference φ1n from pulse to pulse with φ1n≧0;
pulses of the second pulse sequence have a linear increasing phase value with a fixed phase difference φ2n from pulse to pulse, with φ1n≠φ2n and φ1n≧0;
means for Fourier transforming received base band signals of each individual pulse sequence;
means for determining amplitude peaks of resulting two Fourier spectra;
means for clustering the amplitude peaks of the two Fourier spectra to determine pairs of spectral lines of a Doppler frequency shift of
Δ f = ϕ 2 n - ϕ 1 n 2 π . PRF ;
means for calculating target Doppler frequency based on the two Doppler frequencies of a previously determined spectral line pair; and
means for calculating target range based on a phase difference between spectral lines of a previously determined spectral line pair.
Description
BACKGROUND AND SUMMARY OF THE INVENTION

[0001] This application claims the priority of 011 27 990.8, filed Nov. 24, 2001, the disclosure of which is expressly incorporated by reference herein.

[0002] The invention relates to a method for High Pulse Repitition Frequency (HPRF) radar measurement of target range and Doppler frequency, such as may be used, for example, in an airborne radar system.

[0003] The use of a well known radar waveform with two slightly different transmit frequencies is often referred to in the literature as a frequency shift keying (FSK) technique. Using such a technique, target range and velocity can be measured unambiguously even in multi target situations. FIG. 1 shows the two carrier frequencies f1 and f2 as a function of time. Each frequency represented by a vertical line corresponds to a coherent processing interval (CPI), and includes a pulse sequence of N-2048 pulses each.

[0004] In each CPI of length T, a received pulsed signal is directly down converted with the individual instantaneous transmit frequency, and shows as a-single information the Doppler frequencies of the observed targets. The two base band signals are Fourier transformed, and target line left out will be observed in two adjacent CPI's in the same Doppler bin, but with different phase information. For a simultaneous range measurement, these two phases φ12 in the same Doppler bin are considered and the phase difference Δφ=φ2−φ1 is calculated. The latter is directly related to the target range R in accordance with the following equation: Δφ = 2 π Δ f 2 R c Eq . ( 1 )

[0005] The range measurement is unambiguous if Δφ<2π. In order to measure target ranges up to Rmax=150 km, the frequency shift in the FSK wave form must be: Δ f c 2 R = 1 kHz .

[0006] The basic idea is applicable to the extent that the assumptions of a purely stationary observation are fulfilled. But each target maneuver inside the measurement time 2T will have a large influence on the phase measurement in the adjacent CPI's, and will cause large range errors in the range estimation procedure in accordance with Eq. 1.

[0007] One object of the present invention is to provide a method for simultaneously and unambiguously measuring target range and target Doppler frequency for highly mobile targets, even in multiple target situations.

[0008] This and other objects and advantages are achieved by the present invention, in which a transmit oscillator generates a single unmodulated frequency. The HPRF pulse sequence is split into two interleaved parts, which are distinguished and modulated with a digital phase shift keying (PSK) modulation scheme as follows:

[0009] The pulses of the first pulse sequence have the same phase state or alternatively have a linearly increasing phase value with a fixed phase difference φ1n from pulse to pulse; and

[0010] The pulses of the second pulse sequence have a linearly increasing phase value with a fixed phase difference φ2n from pulse to pulse with φ1n≠φ2n and φ1n>0.

[0011] The latter step is important because it renders the down conversion process simple. The PSK technique guarantees unambiguous range measurement and the interleaved pulse sequences make the measurement robust against target maneuver. Before this background all important requirements can be fulfilled with the PSK waveform according to the invention.

[0012] Other objects, advantages and novel features of the present invention will become apparent from the following detailed description of the invention when considered in conjunction with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

[0013]FIG. 1 is a graphic representation of the transmit frequency for FSK-radar as a function of time;

[0014]FIGS. 2.a and 2.b are constellation diagrams which show the phase steps from pulse to pulse for the transmit signal according to the invention;

[0015]FIG. 3 shows the transmit waveform according to the invention, with the corresponding pulses to the phase steps of FIG. 2;

[0016]FIG. 4 is a schematic diagram of a radar transmit/receive system for carrying out the method according to the invention;

[0017]FIG. 5 is a constellation diagram of the transmit and receive signal (base band domain) according to the invention;

[0018]FIG. 6 shows frequency spectra of the received base band signals containing three targets; and

[0019]FIG. 7 illustrates a sampling scheme for signal processing of the received base band signal, wherein m1,c(t), m2,c(t) are the sampling sequences for demultiplexing.

DETAILED DESCRIPTION OF THE DRAWINGS

[0020] The digital PSK waveform can be described easily by the constellation diagram of the phase states of the different pulses in the pulse sequence. In a preferred embodiment a distinction is made between a fixed phase modulation of the first sequence, in which all pulses have the same phase state (for example zero phase) for each pulse, as shown in FIG. 2.a. In this case, the second pulse sequence will have a linearly increasing phase modulation with a fixed phase difference from pulse to pulse, as shown in FIG. 2.b.

[0021] For this digital PSK waveform the transmit signal consists of two interleaved pulse sequences with a fixed transmit frequency, but one having a constant phase (first pulse sequence), and the other being modulated by a linear phase modulation technique (second pulse sequence), as shown in FIG. 3.

[0022] The phase modulation from pulse to pulse can be realized for example by using the quadrature modulation technique which is well known in digital radio transmission. The amplitude envelope of the waveform is constant, so that no distortion due to nonlinearity of the individual amplifier need be considered. The linear phase modulation of the second pulse sequence is carried so that the phase shift from pulse to pulse satisfies the condition that in the complete pulse sequence the phase rotates exactly (for example) 20 times in the constellation diagram. This condition leads to a shift in the Doppler frequency domain of exactly 20 FFT bins.

[0023]FIG. 4 shows the transmitter/receiver block diagram of the radar system on the basis of a quadrature (de)modulator.

[0024] The transmit signal st(t) can be described analytically as follows, wherein p(t) is the modulator function and s(t) is the complex video signal:

p(t)=e j2πf0t

s(t)=s 1(t)+js Q(t)=ae jφ(t)

s t(t)=p(t)s(t)

s t(t)=ae f(2πf0t+φ(t)) , f 0 transmit frequency

with

φ(t)=φ1(t)=0 for 2nPRI12 ≦t<(2n+1)PRI12 , n Pulse number  Eq. (2)

φ(t)=φ2(t)=int(t/PRI−0.5)φ2n for (2n+1)PRI12 ≦t<(2n+2)PRI12  Eq. (3)

[0025] The received signal is down converted with the fixed instantaneous carrier frequency f0 of the transmit oscillator signal p(t) that is also used as the demodulator signal. The demondulator signal is then the complex conjugate of p(t): S r ( t ) = b a s 1 ( t ) * p * ( t ) = be f ( 2 π f 0 r + ϕ ( t ) ) * - j2π 0 2 = b ( t ) = b cos ϕ ( t ) + jb sin ϕ ( t )

[0026] wherein for φ(t) the above mentioned Eq's 2 and 3 still apply and b refers to the considerably lower amplitude of the receiving signal.

[0027] The advantage of the waveform according to the invention is that the phase constellation diagram of the receive signal is shifted compared with the transmit signal only due to the time delay τ or the target range R respectively, as shown in FIG. 5.

[0028] The demodulated signal is A/D converted as shown in FIG. 4, and the received base band signal is split into the two signal sequences m1, m2 (FIG. 4) corresponding to first and second transmit pulse sequence. The received base band signal m1, m2 of each individual pulse sequence is Fourier transformed, and each spectral line is represented by a complex value in amplitude and phase. The Doppler frequency shift is estimated after target detection. Due to the linear phase modulation in the second pulse sequence a single target will be observed after Fourier transformation with the same Doppler frequency but with an additional frequency shift of Δf (20 FFT bins in the considered example).

[0029] An example representing a 3 target situation is shown in FIG. 6 which contains the frequency spectra of the received base band signals of first (dark gray) and second (light gray) pulse sequences. As can be easily seen from that drawing, for each target there exists a pair of spectral lines (one from the spectrum of the first pulse sequence and one from the spectrum of the second pulse sequence) with a certain frequency shift Δf.

[0030] In order to measure target range and speed simultaneously, the following signal processing steps are performed:

[0031] calculate Fourier transformation of both receive signals;

[0032] search the maximum amplitude peaks in both spectra;

[0033] cluster the maxima of both spectra which have a frequency shift of Δf (Doppler bin pair);

[0034] calculate Doppler frequency (or target speed) for each Doppler bin pair; and

[0035] calculate phase difference for each Doppler bin pair and calculate target range.

[0036] These steps will now be described in detail based on a preferred embodiment of the invention.

[0037] Transmit signal st(t) generation and the mixing process are well known and require no further explanation. The received demodulated signal sr(t) is sampled with fa=4 MHz sampling frequency which leads to a range gate width RGW of 37.5 m. The total sampling process of the base band signal is depicted in FIG. 7.

[0038] The receive sampled sequence inside a single range gate is considered. The sample values have a time distance of PRF12. This signal sequence is split into two different interleaved sequences filled with zeros as shown in FIG. 7. The described signal processing in FIG. 7 is for demultiplexing the two received interleaved pulse sequences m1(t) and m2(t), where it is not known if the sequence m1(t) contains the echo signal of the transmit sequence with constant phase states and m2(t) the echo signal of the transmit sequence with linear phase modulation or vice versa. Thus, it is necessary to distinguish the base band signal for each range gate into two cases described by Equations (4) and (5). m 1 ( t ) = { b exp ( j [ 2 π f dop t + ϕ 1 ( t - τ ) = 0 - 2 π f 0 τ ] ) , 2 n t PRI 12 < 2 n + 1 0 , 2 n + 1 t PRI 12 < 2 n + 2 m 2 ( t ) = { b exp ( j [ 2 π f dop + ϕ 2 ( t - τ ) = 2 π ( f dop + Δ f ) t + div ( τ , PRI ) - 2 π f o τ ] ) , 2 n + 1 t / PRI 12 < 2 n + 2 0 , 2 n t / PRI 12 < 2 n + 1 ( Eq . ( 4 )

[0039] The amplitude b of the target echo signal is calculated in accordance with the radar cross section (RCS) and in accordance to the radar equation. m 1 ( t ) = { b exp ( j [ 2 π f dop t + ϕ 2 ( t - τ ) = 2 π ( f dop + Δ f ) t + div ( τ , PRI ) - 2 π f o τ ] ) , 2 n t / PRI 12 < 2 n + 1 0 , 2 n + 1 t / PRI 12 < 2 n + 2 m 2 ( t ) { b exp ( j [ 2 π f dop t + ϕ 1 ( t - τ ) = 0 - 2 π f 0 τ ] ) , 2 n + 1 t PRI 12 < 2 n + 2 0 , 2 n t PRI 12 < 2 n + 1 Eg.(5)

[0040] Inside each individual range gate the two received and demultiplexed sequences m1(t) and m2(t) are first weighted using, for example, a Blackmann or a Hanning window, denoted as w(n) and Fourier transformed by an FFT. The resulting Doppler frequency spectra are denoted by M1(t) and M2(t).

M 1(f, rangegate)=FFT(m 1(nT, rangegate)w(n))

M 2(f, rangegate)=FFT(m 2(nT, rangegate)w(n)).

[0041] For each range gate spectrum M1(f) a threshold is calculated, using a CA-CFAR algorithm (CA-CFAR: cell averaging-constant false alarm rate) such as disclosed in German patent document DE 196 00 779 A1, which is well known to those skilled in the art. The Doppler frequency of the k-th detected target is denoted by f1,k (number of FFT bin), the measured phase information on this FFT bin by φ1,k and the range gate number by rgk. This information is needed for further processing.

[0042] For all detected target amplitude peaks of M1(f) there are existing corresponding frequency peaks in the same range gates but in the spectrum M2(f) . These frequencies are shifted due to the linear phase modulation by 20 FFT bins in the considered examples. The absolute frequency shift is given by: Δ f = ϕ 2 n - ϕ 1 n = 0 , PRF

[0043] A binary factor mpos is introduced, which describes the two situations, in which m1(t) contains the echo signal of the constant phase modulation and m2(t) the echo signal of the linear phase modulation (mpos=1) or vice versa (mpos=−1).

[0044] Frequency peaks in both spectra M1(f) and M2(f) for the same range gate are correlated, if the detected frequency peaks are around 20 FFT bins apart. This situation can be described analytically by the following equation. f 2 , k , rg - f 1 , k , rg Δ f = m pos , k { - 1 , 1 }

[0045] In the case of a frequency pair f1,k,rg, f2,k,rg and the related phases φ1,k, φ2,k target speed and range can be measured unambiguously and simultaneously even in multi target situations by the following simple equations.

[0046] The target Doppler frequency can be calculated by the measured frequencies based on the knowledge of which receive pulse sequences contains the zero phase stages (mpos=1, or mpos=−1) f dop , k , rg = { f l , k , rg m pos , k = 1 f 2 , k , rg m pos , k = - 1

[0047] The radial velocity of the detected target is in this case: v k , rg = - f dop , k , rg 2 c f 0

[0048] The target range R can now be calculated simply by the phase difference directly and the given range gate information:

Δφk2,k−φ1,k

[0049] The target range ambiguity factor tramb, which has a integer value, can be calculated by: tr amb , k = int ( Δφ k ϕ 2 n - ϕ 1 n ) = int ( Δ φ k ϕ 2 n )

[0050] Using the target range ambiguity factor and the range gate information, the unambiguous target range is given by: R k = { ( rg - 0.5 ) RGW + tr amb , k R amb | m pos , k = 1 ( rg - 0.5 ) RGW + ( tr amb , k - 0.5 ) R amb | m pos , k = - 1

[0051] with rg being the number of range gate, where the target has been detected

[0052] RGW being the width of range gate

[0053] Ramb being the unabiguous range measurement interval given by PRI.

[0054] System Parameters

[0055] The table below shows typical system parameters for an airborne radar system in which the method according to the invention can be used.

Carrier
Frequency r 8.5 . . . 10.7 GHz FFT-Length 2048
Relative −1500 m/s . . . Unambiguous 1.5 km
Velocity νrel  +300 m/s Range Ramb
Doppler  90 kHz Range Gate Width 37.5 m
Frequency (−18 kHz . . . 72 RGW
Bandwidth kHz)
Distance 0 km . . . 150 km Noise Figure F 5
Measurement
Range R
Antenna 32 dB Noise Bandwidth 4 MHz
Gain G
Sampling 4 MHz Noise Temperature 293 K
Frequency α
Pulse 1/(4 MHz) = 250 ns Transmit Power 10 kW
Width Tpw
Number of 2048 Target RCS 5 m2
Pulses per
Sequence
(CPI)
PRF 100 kHz PRF12 200 kHz

[0056] The foregoing disclosure has been set forth merely to illustrate the invention and is not intended to be limiting. Since modifications of the disclosed embodiments incorporating the spirit and substance of the invention may occur to persons skilled in the art, the invention should be construed to include everything within the scope of the appended claims and equivalents thereof.

Referenced by
Citing PatentFiling datePublication dateApplicantTitle
US7652618 *Jul 1, 2008Jan 26, 2010Deere & CompanyMethod and system for determining velocity by using variable or distinct sampling rates
US8705953 *Mar 10, 2010Apr 22, 2014Carl Zeiss AgDistance measurement
US20130004154 *Mar 10, 2010Jan 3, 2013Carl eiss AGDistance measurement
US20130278455 *Apr 23, 2012Oct 24, 2013Massachusetts Institute Of TechnologyEfficient pulse doppler radar with no blind ranges, range ambiguities, blind speeds, or doppler ambiguities
DE102004054198A1 *Nov 10, 2004May 18, 2006Eads Deutschland GmbhVerfahren zur Detektion und Entfernungsbestimmung von Zielen mittels einer HPRF-Radaranlage
DE102004054198B4 *Nov 10, 2004Jan 18, 2007Eads Deutschland GmbhVerfahren zur Detektion und Entfernungsbestimmung von Zielen mittels einer HPRF-Radaranlage
Classifications
U.S. Classification342/109, 342/135, 342/131, 342/111, 342/196, 342/137, 342/115, 342/110
International ClassificationG01S7/292, G01S13/524, G01S13/26, G01S7/288
Cooperative ClassificationG01S13/26, G01S2007/2883, G01S7/2923, G01S13/524, G01S7/288
European ClassificationG01S7/292C, G01S7/288, G01S13/26, G01S13/524
Legal Events
DateCodeEventDescription
Jul 29, 2011FPAYFee payment
Year of fee payment: 8
Jul 26, 2007FPAYFee payment
Year of fee payment: 4
Apr 28, 2003ASAssignment
Owner name: EADS DEUTSCHLAND GMBH, GERMANY
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ROHLING, HERMANN, PROF. DR. DR.;GRUENER, WILHELM DR.;HOMMEL, HANS;AND OTHERS;REEL/FRAME:014009/0378;SIGNING DATES FROM 20030114 TO 20030127
Owner name: EADS DEUTSCHLAND GMBH WILLI-MESERSCHMITT-STRASSE85