Target Detection
Target Detection
Radar Fundamentals
1.1 INTRODUCTION
This chapter introduces some of the major radar concepts that will be described
in greater detail in subsequent chapters. One key concept, however, is covered
primarily in this chapter: the radar range equation (RRE). Although specific
forms of the RRE are discussed in other chapters, its derivation and definitions
of all major terms are covered here. The remainder of the chapter provides an
overview of other significant radar concepts.
• Search:
– To survey a volume of space and report locations of targets
– To measure target position in two or three dimensions (i.e., 2-D or 3-D)
1
2 Chapter 1
• Tracking:
– Obtain more accurate target location by “smoothing” measurements
– Estimate the “state vector” of the target (i.e., positions, rates, and possibly
accelerations) to predict the target state vector for a future time
• Track-While-Scan:
– Combine the search and tracking functions into one radar mode
– Use data processing to initiate and maintain tracks, while simultaneously
searching for new targets
– Add tracking capability without using additional radar resources.
The reflected radar energy (or waveform) is subject to additive thermal noise,
primarily due to active electronics in the receive chain. It is in this “noisy” envi-
ronment within which targets must be sensed. This radar return or “echo” is de-
picted in Figure 1.1.
Two key target detection attributes are the “probability of detection” (Pd) and
“probability of false alarm” (Pfa). The radar’s design objective is to maximize Pd
while maintaining a small (and possibly constant) Pfa.
Detectability (or Pd) depends on the ratio of reflected target energy to average
thermal noise power, which is defined as the signal-to-noise ratio (SNR). SNR is
an important radar performance metric. The performance of nearly all radar
functions depends on the SNR. For example, Pd increases monotonically with in-
creasing SNR.
Two approaches to increasing SNR are to use higher-energy (i.e., larger ampli-
tude or longer duration) waveforms or to add (or “integrate”) multiple return
pulses, either coherently (in-phase) or non-coherently (summing magnitudes or
without phase coherence).
A target-present decision when no target (only noise) is present is called a
“false alarm.” The converse to this error is declaring “no target present” when a
target is actually present. This latter error is referred to as a “missed detection.”
Multiple pulses can be transmitted as a pulse “train” or “burst.” This type of
waveform is depicted in Figure 1.2. If each pulse in a burst has energy E, then
the waveform illustrated in the figure has total energy NE.
The minimum “non-self eclipsing” range is based on the pulse duration τ and
is given by:
, (1.1)
where c is the speed of light. This phenomenon arises from the fact that recep-
tion cannot begin until transmission is complete. The maximum “un-eclipsed”
range is therefore defined as:
T −τ ,
R max = c (1.2)
uneclipsed 2
4 Chapter 1
where T is defined as the pulse repetition interval (PRI). Alternatively, the pulse
repetition frequency (PRF) is defined as:
1 1 . (1.3)
PRF = =
T PRI
1. (1.4)
B =
τ
Figure 1.3 shows an ideal transmitted rectangular radar pulse. Two targets (of
equal amplitude) separated by much greater than τ in time are easily “resolved.”
The inherent range resolution capability of the waveform is defined as:
c .
δR = (1.5)
2B
The reflected energy received by the radar may originate from unintended ob-
jects such as land or sea backscatter, weather (e.g., rain), or man-made objects
such as buildings or other structures. This type of radar return is referred to as
“clutter.” When the radar is used for weather detection, clutter is the desired tar-
get.
A key discriminating feature used to differentiate targets from clutter is the
observed or measured speed (or range-rate). This method is based on exploiting
the Doppler-shifting effect, that is, the phenomenon that returns from higher
, (1.6)
where Ṙ is the target or clutter range-rate and λ is the radar operating wave-
length.
A second key discriminating feature used to separate targets from clutter is
the “polarization” of the reflected electromagnetic (EM) signal. This is especially
true for clutter arising from rain. The return from spherical rain drops when us-
ing a circularly polarized waveform is reduced compared to the target response
when a linearly polarized waveform is employed. Different polarizations will
reflect (or refract) distinctly from different shaped targets; for example, sharp
edges tend to spread energy into the primary and orthogonal polarizations,
whereas smooth reflectors result in a single reflected polarization. Polarization is
an important feature for separating certain classes of targets.
waveform generator, exciter and transmitter, receiver, signal processor, data pro-
cessor, and operator displays and controls.
Note that for reflector-antenna radars, the phased-array is replaced with a re-
flector antenna. Also for mechanically scanned radars, the beam steering genera-
tor is replaced with antenna pedestal control. In the case of solid-state transmit-
receive configuration, the transmitter is replaced with T/R modules (this type of
radar is referred to as active-aperture radar). Most often in modern radars, the
signal processing is performed in software, and the signal processor and data
processor are combined into a single signal/data processor (i.e., computer or
multiprocessor).
Figure 1.6 is a block diagram of a phased-array antenna (also known more
generally as an electronically-steered array).
A phased-array antenna is a “sampled aperture,” where the antenna consists
of N small antennas or antenna elements. Each is “phased-steered” to receive or
transmit from an angle θ steer using phase shifters. The phase shifters are usu-
ally digitally selected to obtain quantized phase (e.g., 0, 22.5°, 45°, 67.5°, 90°).
This type of antenna is used for narrowband operation. When wideband opera-
tion is required, some combination of time-delay and phase-shift steering is
needed.
Radar Fundamentals 9
Pt Gt
Transmit Power Density = , (1.7)
4π R 2 Lt
where Pt Gt is the product of the peak transmitter power and the transmit an-
tenna gain, and R2 Lt is the product of range of the target from the radar squared
and the total transmit losses. For a target with an RCS of σ, the echo or reflected
power density is given by:
10 Chapter 1
Pt Gt σ
Echo Power Density = . (1.8)
( 4π R )
2
2
Lt
The received power at the radar antenna aperture of area Ar and receive loss of
Lr is then:
Pt Gt σ Ar
Power at Radar Aperture = . (1.9)
( 4π R )
2
2
Lt Lr
When the thermal front-end noise can be modeled as having a “white” spectrum
with power spectral density of kTs, that is, the product of Boltzmann’s constant
and the system noise temperature, and the radar has a processing (or noise)
bandwidth of B, the noise at the output of the matched filter is given by:
Pt Gt Ar σ
SNR = . (1.11)
( 4π R2 )
2
k Ts B Lt Lr
4π Ar , (1.12)
Gr =
λ2
then solving for the receive aperture and substituting into equation (1.11) yields
the commonly used “sensitivity” form of the RRE:
Pt Gt Gr λ 2 σ
SNR = . (1.13)
( 4π ) 3 4
k Ts B R Lt Lr
S S 1
SIR = ≈ = . (1.15)
I N+ I ( SNR ) −1
+ ( SIR )
−1
12 Chapter 1
Therefore, when the interference level is much lower than thermal noise level,
the equations (1.11) and (1.13) are a useful performance figure-of-merit (FoM) for
radars. However, when the interference level approaches the noise level and ex-
ceeds it, the SIR defined in equation (1.15) becomes the performance FoM.
As can be seen from equation (1.15), when the thermal noise is exceeded by
the thermal noise-plus-interference, target detection sensitivity is reduced and is
now limited by the SIR. Mathematically, this can be expressed as:
2 2
P J G J G rj λ P J G J G rj λ
kT s B → kT s B + -------------------------
- ≈ -------------------------
-, (1.16)
2 2
( 4πR J ) ( 4πR J )
where PJ GJ is the product of the jammer power and antenna gain (referred to as
the effective radiated power of the jamming source), RJ is the range of the jam-
mer to the radar, Grj is the radar receive antenna gain in the jammer direction,
and jammer bandwidth ≥ B. Note that the jamming energy impinging on the ra-
dar antenna aperture decays as an inverse-range-squared law, compared with
the target return, which goes as inverse-range to the fourth power. This enables
relatively low-powered jamming sources to degrade radar performance, which
is the inherent advantage of the jammer over radars.
Substituting equation (1.16) into (1.13) yields the RRE for the case where jam-
ming levels far exceed thermal noise:
Pt Gt Gr λ 2 σ RJ2
SIR = . (1.17)
( 4π ) 2
PJ GJ R 4 Lt Lr
1.6.2.1 Radar Range Equation for Volume Search. For volumetric searches as de-
picted in Figure 1.4, the RRE is adapted in the following manner. First, note that
the antenna beamwidth is related to the antenna gain via the relationship:
14 Chapter 1
λ λ π
θ3 ≈ = = 2
A λ2G G . (1.18)
4π
π .
Ab ≈ θ3 AZ θ3 EL = 4 (1.19)
G
Also, the average transmitted power can be defined as:
Pt PRF . (1.20)
PAVE = Pt DF = Pt τ PRF =
B
where DF is the radar duty factor (percentage of time allowed to transmit) and
PRF is the pulse repetition frequency. Noting that:
Ψ ΨG
= = PRF (or beams / s) , (1.21)
θ3 AZ θ3 EL TSC 4π TSC
where Ψ is the area to be searched in radians2 and TSC is the scan or “frame”
time for the search, and substituting equation (1.20) into (1.21) yields:
ΨG B PAVE ,
= PRF = (1.22)
4π TSC Pt
or:
PAVE 4π TSC B .
Pt = (1.23)
ΨG
Now substituting equation (1.23) into the RRE of equation (1.13) results in:
Equation (1.24) is proportional to the PAVE Ar product and is not a function of op-
erating frequency. Also note in contrast that equation (1.11), the sensitivity form
of the RRE, is proportional to the product Pt Gt Ar . Theoretically then, radars at
any operating frequency can search a volume equally well if they possess the
same PAVE Ar .
In practice, however, lower operating frequencies are usually employed for
search radars since the number of beams required to search a given volume, for
a fixed antenna aperture size, is many fewer than at higher frequencies. At
higher frequencies this can lead to timeline occupancy issues since the larger
number of beams required will demand that longer portions of the timeline be
scheduled. When these scheduling periods exceed the specified frame time, the
radar is said to be “occupancy-limited.”
1.6.2.2 Radar Range Equation for Horizon Fence Search. Early warning (EW) mis-
sile surveillance radars and ballistic missile defense (BMD) radars often employ
horizon fence searches to detect and acquire targets. Figure 1.9 depicts a typical
search fence that covers ± 60 degrees of azimuth.
The fence concept is based on the fact that if radars have adequate detection
range on the targets of interest, then any ascending ballistic target must fly
through a fence and will be detected. Therefore, instead of performing a radar
resource-intensive volume search (i.e., due to the much larger number of
beams), a single row of beams at or above the horizon is sufficient for missile
search and acquisition applications.
Starting with the RRE of equation (1.24) for volume search, modifications can
be made to take advantage of the horizon fence characteristics. First note that for
missiles flying through the fence:
θ3 θ3 R
TSC = = , (1.25)
E T N vT
where Ė T is the target elevation rate, R is the target range, νT is the vertical tar-
get velocity, and N is the number of looks required for detection. Noting the
beamwidth relationship in equation (1.18), (1.25) can be expressed as:
16 Chapter 1
2R π .
TSC = (1.26)
N vT G
σ PAVE Ar
SNR = . (1.27)
( 2 π ) k Ts R3 Ψ N vT Lt Lr Gr
As can be seen from equation (1.27), the RRE for horizon fence search is weakly
a function of operating frequency due to the G term and is also follows an in-
verse R3 law rather than R4 as do the RREs in equations (1.13) and (1.24).
1.6.2.3 Radar Range Equation for Tracking. The primary driving requirement on
tracking is angle accuracy given by:
θ3
σθ = , (1.28)
km 2 SNR η
where km and η are the monopulse slope and number of independent measure-
ments smoothed by the tracking filter. The SNR is defined by the track sensitiv-
ity form of the RRE given by equation (1.11). Substituting average power for
peak power as defined in equation (1.20) yields:
PAVE Gt Ar σ
SNR = , (1.29)
)
PRF ( 4π R 2 2 k Ts Lt Lr
Radar Fundamentals 17
θ32 PRF ( 4π ) R k Ts Lt Lr .
24
σ 2
= (1.30)
θ
2 km2 η PAVE Ar Gt σ
Now substituting for antenna beamwidth from equation (1.18) and noting that η
= PRF Tt where Tt is the time in track (and PRF refers to the track update rate)
yields:
σ 2
=
( 4π )
3
k Ts R 4 Lt Lr . (1.31)
θ
2 km2 Tt PAVE Ar Gt Gr σ
1.6.2.4 Radar Range Equation Summary. Tables 1.1 and 1.2 provide a summary of
the RREs for the radar applications discussed in this section.
This topic is covered in detail in Chapter 2 for detection in noise, clutter, and
jamming environments, respectively. However, some of the basic concepts are
introduced in this chapter.
In general, detection performance is a function of SNR (or in cases of clutter or
jamming, signal-to-clutter ratio (SCR) or SIR, respectively) and typically em-
ploys threshold tests to declare target detections. References [6] and [7] address
the theory for this critical radar function. Most modern radars use matched-filter
receivers to maximize the SNR at the processed output prior to detection deci-
sions.
18 Chapter 1
Radar
Radar Appropriate Form of the Parametric
Application Radar Range Equation Factor
σ TSC
SNR = PAVE Ar
( 4π ) k Ts R4 Ψ Lt Lr
Volume search PAVE A
σ PAVE Ar PAVE A
SNR =
Horizon fence search
( 2 π ) k Ts R Ψ N vT Lt Lr
3
Gr G
Pt Gt Ar σ
SNR = Pt AG
( 4π R )
Track sensitivity 2
2
k Ts B Lt Lr
Track accuracy σ 2
=
( 4π ) 3
k Ts R 4 Lt Lr
PAVE AG 2
θ
2 k Tt PAVE Ar Gt Gr σ
2
m
Radar
Radar Appropriate Form of the Parametric
Application Radar Range Equation Factor
σ TSC
Volume search SNR = PAVE Ar PAVE A
16 k Ts R 4 Ψ Lt Lr
πσ PAVE Ar PAVE A
Horizon fence search SNR =
16 k Ts R3 Ψ N vT Lt Lr Gr G
Pt Gt Ar σ
SNR = Pt AG
( 4π R )
Track sensitivity 2
2
k Ts B Lt Lr
Track accuracy σ 2
=
( 2π ) 4
k Ts R 4 Lt Lr
PAVE AG 2
θ
2 k Tt PAVE Ar Gt Gr σ
2
m
Radar Fundamentals 19
H1
s(t) > VT , (1.32)
< VT
H0
where s(t) is the output of the optimal matched filter and VT is the detection
threshold. The detection threshold is usually based upon an assumption of
Rayleigh-distributed random noise magnitudes.
The target’s RCS is an approximation devised to account for the portion of scat-
tered energy reflected back toward the radar from the target. The target RCS
may vary in amplitude (or fluctuate) versus time. The RCS depends on the RF
operating frequency of the radar (e.g., ultra-high frequency [UHF], L-band, X-
band) and on the reflecting structure of the particular target. Complex scattering
surfaces such as an aircraft are composed of many individual RF scatterers. Sim-
pler target shapes, like a conical missile warhead, may present only one or two
scattering sources.
The interpretation of the reflected energy scattered by the target then becomes
driven by the waveform range resolution:
c cτ eff . (1.33)
δR = =
2B 2
20 Chapter 1
When δR is ≥ target “length,” a single return occurs at the radar that is a sum of
the complex-valued individual scattering components. RCS fluctuations arise
due to constructive and destructive combination of individual scattering center
components within the resolution-limited waveform bandwidth.
However, when δR << target “length,” multiple scattering returns are re-
solved in range (resulting in little or no RCS fluctuation). Models of target RCS
fluctuation (e.g., Swerling Models I-IV, log-normal) were developed to allow
detection analysis to be performed without having actual RCS data versus
viewing angle.
The probability of false alarm affects the number of “false target responses” that
must be evaluated by the radar processor. If a search volume such as that in Fig-
ure 1.10 consists of Nb antenna beams, where Nb = Nazimuth × Nelevation, and Nr
range gates, then the average number of false alarms per frame or search volume
is:
N FA = PFA ⋅ N b ⋅ N r . (1.34)
Typically, the search PFA is selected to result in ≤ NFA/second. Therefore, for ex-
ample, if NFA/second = 5 false alarms/second, and Nb = 100 and Nr = 1000, then
for a scan or frame time (TSC) of 2 seconds, the PFA allowed is:
N FA ⋅ TSC ( 5) ( 2 )
PFA = = = 10 −4 . (1.35)
Nb ⋅ Nr 10 5
fX ( x ) dx ,
∞
PD = ∫ VT
(1.36)
Radar Fundamentals 21
where fx (x) is the RCS fluctuation density, and VT is the selected detection
threshold, and:
VT2
∞ −
PFA = ∫ fn ( n ) dn = e 2σ 2 , (1.37)
VT
(
VT2 = − n Pfa 2σ 2 . ) (1.38)
Substituting (1.38) into (1.36) allows the calculation of PD for any analytical RCS
model.
Calculating PD for a specified PFA and target RCS model requires performing the
integration defined in (1.36). For exponentially distributed RCS (e.g., Swerling I),
the integral has a closed-form solution:
PD = (P )
FA
1 + SNR . (1.39)
22 Chapter 1
n PFA
SNR REQ = − 1. (1.40)
n PD
For probability densities in (1.36) that do not have a closed-form solution, nu-
merical integration can be used to evaluate PD. PD has been evaluated for many
target models and is documented in texts such as references [6] and [7].
VT 2 = − 2 σˆ 2 ln PFA . (1.41)
Radar Fundamentals 23
Due to the variance on the total noise estimate, use of equation (1.41) in place
of equation (1.38) incurs a loss in detectability referred to as a CFAR loss. For this
reason, CFAR tends to be used where the additional signal processing loss (in-
cluded in the denominator of the RRE as a component of Lr) can be tolerated.
Therefore, it is usually used for tracking but not for functions such as long-range
search where the additional loss is undesirable.
where:
2 ∆R , (1.42)
∆ fd =
λ
and ∆Ṙ and λ are the difference between the target or clutter range-rates and the ra-
dar operating wavelength, respectively. A typical case is illustrated in Figure 1.12.
For example, land clutter is usually distinguishable from moving targets
since the mean velocity of the land reflections is zero. However, Doppler shifts
for rain in wind conditions or sea clutter at high sea state can overlap with
those of slow-moving targets (especially surface targets like ships or tanks). As
a benefit to target detection, certain RF waveform polarizations, for example,
circular, can result in reduced reflections from rain clutter as compared with
those of the target.
There are two basic approaches to mitigating clutter effects to better detect
moving targets:
MTI cancellers operate by subtracting sequential radar returns with the objec-
tive that near-stationary clutter will be cancelled, but moving targets will not.
Figure 1.12 Doppler Frequency Separation of Target from Land and Sea
Clutter
Radar Fundamentals 25
The simplest MTI canceller is the single-delay (or two-pulse) canceller depicted
in Figure 1.13. The clutter appearing at the canceller output is reduced substan-
tially by the attenuation of low-frequency clutter. Optimum performance is ob-
tained when target Doppler shift is near PRF/2 (i.e., least attenuation of target
response).
Pulse-Doppler waveforms are coherent bursts that consist of N pulses spaced
by a uniform delay of PRI (pulse repetition interval). A pulse-Doppler pulse
train waveform is shown in Figure 1.14.
The matched filter to a Pulse-Doppler waveform consists of sub-pulse
matched filtering (i.e., range processing), with the results stored for M range
cells, followed by N Doppler filtering. Clutter (and its “alias”) will appear in
lower and higher Doppler filters. This is depicted in Figure 1.15.
c
Range : δ R =
2B
Angle : δθ = θ3 , (1.43)
1
Doppler : δ f =
T
where:
τ = pulse length
B = pulse bandwidth
c = speed of light
θ3 = 3 dB antenna beamwidth
T = integration time.
δR
Range accuracy ≅ = σR
2 SNR
δθ
Angle accuracy ≅ = σ θ ; km ≈ 1.6 . (1.44)
km 2 SNR
δf
Doppler accuracy ≅ =σf
2 SNR
Tracking radars predict the future position of targets by estimating range and
angular rates. These predictions are used to position the antenna beam for sub-
sequent radar transmit and receive actions.
Monopulse is a technique to measure target angles by separating two-dimen-
sional receive antennas into azimuth and elevation “quadrants”; and combining
them in a specific way to enable the estimation of angular positions. An “error”
pattern is formed as e ( θ ) = ( θ ) ⁄ Σ ( θ ) as shown in Figure 1.17. The amplitude
and sign of e ( θ ) indicates the distance of the target from the antenna pointing
angle.
• Usually, the previous two or three scans of data are used to perform track initi-
ation (TI)
• Subsequently, tracks are maintained (or updated) on each new scan (re-visit)
of a particular angular region.
1.11 REFERENCES
[1] R. Nitzberg, Radar Signal Processing and Adaptive Systems, Artech House, 1999
[2] D. K. Barton, Modern Radar System Analysis, Artech House, 1988
[3] D. R. Wehner, High Resolution Radar, Artech House, 1987
30 Chapter 1
1.12 PROBLEMS
1. Consider a radar with a peak power of 100 kW, 45 dB transmit antenna gain,
a 0 dBsm target, a 10 square meters antenna aperture, a target range of
500 km, with transmit and receive losses of 3 dB each. Estimate the received
power at the radar aperture using the relationship below.
t tσ Ar
PG
Prec = .
( 4 π R2 )
2
Lt Lr
2. Consider a radar being degraded by barrage noise jamming. For a radar re-
ceiver with system noise of –143 dBm, a jammer effective radiated power
(PJGJ) of 10 W, a 0.01 square meter effective radar antenna aperture (equiva-
lent to a –30 dB sidelobe), with the jammer at 1,000 km:
PJ GJ Ar'
N effective = kTs Br + .
4π RJ
How much degradation in equivalent noise power does the jammer pro-
duce? Using the radar from problem 1, estimate the effective signal-to-inter-
ference ratio (SIR). Is this good or bad for the radar’s performance?
3. Consider a radar with 2 false alarms per second, 1,000 range cells, 16 beams,
and a 10-second search frame time. What probability of false alarm does this
require? What if the search requires 32 beams and a 5-second frame time?
N FA ⋅ Τ
PFA = .
Nb ⋅ Nr
Radar Fundamentals 31
n Pfa
SNR REQ = − 1
1
PD = ( PFA ) 1 + SNR
n Pd
5. If a probability of detection of 0.95 is required for the system, what SNRs are
required to achieve these for the probabilities of false alarm calculated in
problem 3? Are these SNRs reasonable for a radar? Why or why not?
6. Consider a radar with range, angle (both azimuth and elevation), and Dop-
pler resolutions of 15 meters, 20 milliradians, and 100 Hz. If an SNR of 10 dB
is available, calculate the measurement accuracies that can be expected. If
the original target was a fighter aircraft with a 10 dBsm RCS, what accuracies
would be achieved for a missile with a 0 dBsm RCS and for a 747 aircraft
with a 25 dBsm RCS? Are these reasonable?
δR
Range accuracy ≅ = σR
2 SNR
θ3
Angle accuracy ≅ = σθ
km 2 SNR
δf
Doppler accuracy ≅ =σf
2 SNR
7. If an LFM or “chirp” waveform with a 10 µsec pulse length and a 1,000 MHz
bandwidth is used, what is the effective range resolution? What range accu-
racies would be achieved for the three target types in problem 5?