A Comprehensive Study Guide for the Radar & Satellite Engineering Course
ECE 521E • 5th Year BSc. Electrical & Communication Engineering • Masinde Muliro University of Science and Technology
Unit Overview
This study guide supports Unit: Moving Target Indicator (MTI) Radar of ECE 521E. Every surveillance radar that looks toward the earth's surface must extract aircraft, ship, or vehicle echoes from echoes that are vastly stronger — returns from hills, buildings, sea waves, rain and chaff. This undesired return is clutter. MTI processing is the classical, and still indispensable, technique for rejecting zero-Doppler clutter while passing moving targets.
Learning Outcomes
After studying this unit and completing the tutorials you should be able to:
Compute Doppler shifts and interpret their sign for approaching/receding targets.
Explain coherent and non-coherent MTI with a block diagram.
Derive the frequency responses of single and double delay-line cancellers.
Compute blind speeds and design staggered-PRF waveforms to raise them.
Design binomial FIR and recursive IIR MTI filters and sketch their notches.
Define and compute the improvement factor and subclutter visibility.
Select amplitude/statistical models for land, sea and weather clutter.
Critically assess MTI limitations and the migration path to MTD/STAP.
Indicative Study Time
Activity
Time
Guided reading (Sections 1–9)
4–5 hours
Worked examples reworked by hand
2 hours
Tutorial problems (Section 10)
3 hours
Revision of key equations & glossary
1 hour
Prerequisite knowledge: radar range equation, pulsed-radar timing (PRI, PRF, duty cycle), unambiguous range, Fourier analysis of sampled signals, basic random processes (mean, variance, Gaussian PDF).
Section 1 — Doppler Foundations of MTI
1.1 The Doppler Shift
Relative radial motion between radar and target shifts the echo frequency by the two-way Doppler amount:
fd = 2vr/λ = 2f0vr/c
vr > 0 toward the radar → positive shift (echo frequency rises).
vr < 0 receding → negative shift.
Stationary clutter on a stationary platform → fd = 0.
Example 1.1 — A 9.4 GHz (X-band, λ = 0.0319 m) radar sees a target approaching at 250 m/s.
fd = 2(250)/0.0319 ≈ +15.7 kHz. A stationary hill at the same range contributes 0 Hz. MTI exists to separate these two.
1.2 Phase Progression Pulse-to-Pulse
The echo phase is φ = 4πR/λ. Between pulses spaced Tr = 1/fr apart the target moves ΔR = vrTr, so
Δφ = 4πvrTr/λ = 2πfdTr
MTI filtering is fundamentally the measurement of this pulse-to-pulse phase rotation at a fixed range cell. If Δφ is a multiple of 2π, successive echoes are identical and the target masquerades as clutter — the blind-speed condition (Section 4).
1.3 Sampling and Doppler Ambiguity
Pulsed radar samples the target's motion at fr. Hence Doppler is only observable modulo fr: frequencies fd and fd + nfr are indistinguishable. The first blind speed corresponds exactly to fd = fr.
Check Your Understanding 1
Why does the Doppler expression contain a factor of 2?
An S-band radar (λ = 0.1 m) views a receding target producing fd = −800 Hz. Find vr.
Show that a target whose echo phase advances by exactly 2π between pulses produces identical complex samples.
Answers
The wave travels to the target and back, so range (and phase) accumulates at twice the radial rate.
vr = fdλ/2 = (−800)(0.1)/2 = −40 m/s (receding).
x(n) = A e^{j(2πfdnTr+φ)}; if fdTr = 1, the exponent changes by 2πn → x(n) identical for all n → cancels in a differencer.
Section 2 — The MTI Concept
2.1 Signal Model at the Receiver
x(n) = s(n) + c(n) + w(n)
s(n) — target echo: phase rotates by 2πfdTr per pulse.
The delay-and-subtract element is the heart of MTI: identical (stationary) echoes cancel, phase-rotated echoes survive.
Interpretation: MTI is a change detector across pulses. It answers "did anything move in this range cell?" — it does not directly measure how fast (that is pulse-Doppler/MTD territory, Section 9).
Check Your Understanding 2
Why does non-coherent MTI fail when the clutter in a cell is momentarily weak?
What two quantities must the delay line preserve for perfect cancellation?
Answers
The clutter itself supplies the phase reference; weak/absent clutter → no reference → poor subtraction of what little clutter exists and loss of target phase coherence.
Amplitude and phase — equivalently the full complex (I and Q) sample.
Section 3 — Delay-Line Cancellers
3.1 Single Canceller
y(n) = x(n) − x(n−1) ⇔ H(z) = 1 − z⁻¹
|H(f)| = 2|sin(πf/fr)|
Deep null at f = 0; nulls repeat at every multiple of fr (blind speeds).
Response of the single and double cancellers. The double canceller's sin² shape gives a wider, deeper clutter notch; both share nulls at fd = n fr (blind speeds).
3.3 Why the Double Canceller Wins Against Real Clutter
Real clutter is not perfectly stationary — vegetation sways, sea moves. Its spectrum has finite width σc. Near DC, the power response of the double canceller falls as f⁴ versus f² for the single canceller, so far more clutter power is suppressed for the same notch centre (see Fig. C in Section 8).
3.4 Binomial N-Pulse Canceller
y(n) = Σk=0N (−1)k C(N,k) x(n−k)
Weights from Pascal's triangle give a maximally flat notch; N = 2 (three pulses) is the most common practical choice; N = 3 used when weather/chaff is severe. Costs: more pulses per dwell, longer transient, hardware taps.
3.5 Implementation History
Period
Delay element
Notes
1940s–50s
Mercury/quartz acoustic delay
Analog video; fixed PRI
1960s–70s
SAW / ultrasonic lines
Improved stability
1980s–present
Digital RAM at complex baseband
Staggered PRF, adaptive weights, easy redesign
Example 3.1 — PRF = 1200 Hz. Double-canceller response at fd = 300 Hz:
|H| = 4 sin²(π·300/1200) = 4 sin²(π/4) = 4(0.5) = 2.0. The same target at 600 Hz (fr/2, the "optimum" Doppler) gives |H| = 4 — maximum response.
Check Your Understanding 3
Derive |H(f)| = 2|sin(πfTr)| from H(z) = 1 − z⁻¹ by substituting z = e^{j2πfTr}.
A double canceller has weights {1, −2, 1}. Write the four weights of the triple (N = 3) binomial canceller.
Where does the double canceller achieve maximum response, and what Doppler does that correspond to?
Answers
H = 1 − e^{−j2πfTr} = e^{−jπfTr}(e^{jπfTr} − e^{−jπfTr}) = 2j e^{−jπfTr} sin(πfTr) → magnitude as stated.
{1, −3, 3, −1}.
Max at f = fr/2 → fd = PRF/2, the "blind-quadrature" point midway between blind speeds.
Section 4 — Blind Speeds
4.1 Definition and Derivation
Blind speeds are target velocities at which the MTI response vanishes — the moving target cancels exactly like clutter:
vn = n λ fr / 2, n = 0, 1, 2, …
Physical picture: the target moves ΔR = vnTr = nλ/2 between pulses — exactly n half-wavelengths. The echo is then identical on every pulse.
4.2 The Range–Doppler Dilemma
v1Ru = λc/4
With one PRF, pushing blind speeds up (high fr) shortens the unambiguous range Ru = c/(2fr). Both cannot be large simultaneously — the fundamental motive for staggered PRF (Section 5) and PRF diversity.
Safety consequence: an airliner transiting a blind speed disappears from the MTI display. ATC waveform design must guarantee that no realistic aircraft speed coincides with a blind speed — hence staggering is effectively mandatory in surveillance radars.
Example 4.1 — λ = 0.15 m, fr = 600 Hz: v1 = (0.15)(600)/2 = 45 m/s; v2 = 90 m/s. Ru = c/(2fr) = 250 km — comfortable, which is why low-PRF S-band surveillance sets trade-off acceptably without stagger for some ground radars, but airborne/high-PRF sets cannot.
Check Your Understanding 4
Does increasing PRF raise or lower blind speeds? What is sacrificed?
Compute v1 for λ = 0.03 m, fr = 2 kHz, and state the unambiguous range.
Answers
Raises them (vn ∝ fr); unambiguous range shrinks as 1/fr.
v1 = 30 m/s; Ru = 75 km.
Section 5 — Staggered PRF
5.1 Idea
Alternate PRIs T1, T2, … cyclically. Blind speeds of the individual intervals no longer coincide, so no velocity (except zero) is blind in every interval. The composite response for a two-interval double canceller:
Staggering (f₁ = 1000 Hz, f₂ = 1250 Hz, λ = 0.1 m) lifts the first blind speed from 50 m/s to 250 m/s, at the cost of response ripple in the passband.
5.2 Choosing the Stagger Ratio
Ratio T₁:T₂
First blind speed
Passband ripple
Comment
4:5
Moderate
Low
Common choice
7:8
Higher
Moderate
Low ripple for the gain
3:4
High
High
Rarely used alone
5.3 Engineering Costs
Range-gate alignment: target echoes arrive at different times per PRI → digital interpolation before differencing.
Cancellation loss from timing jitter and misregistration.
CFAR normalisation must cope with unequal dwell statistics.
Example 5.1 — λ = 0.1 m, f1 = 900 Hz, f2 = 1200 Hz (ratio 4:3):
vb = (0.1)(900)(1200)/(2·300) = 180 m/s — versus 45 m/s uniform. Ratio 4:3 is aggressive; 5:4 would give 250 m/s? Check: (0.1)(1000)(1250)/(2·250) = 250 m/s with ratio 5:4 — gentler stagger, more practical ripple.
Check Your Understanding 5
Why does staggering not create a new blind speed below vb = λ/[2|T₁−T₂|]?
What single hardware block makes digital staggered MTI easy today?
Answers
The product of responses is zero only where either factor is zero; the first common zero occurs at the difference frequency 1/|T₁−T₂|.
Digital memory (RAM) storing complex baseband samples — delay becomes a read-address shift, trivially variable.
Section 6 — MTI Filter Design
6.1 Transversal (FIR) Structure
H(z) = Σk=0N−1 ak z−k
Design objectives: deep notch at DC, maximally flat or equiripple rejection band, controlled passband gain. The binomial family is the simplest; Chebyshev/equiripple and optimum (Wiener) weights — matched to the measured clutter covariance — give the best subclutter visibility.
6.2 Recursive (IIR) Structure
H(z) = (Σ ak z⁻ᵏ) / (1 + Σ bm z⁻ᵐ)
Sharper notch per pulse; must be checked for stability (poles inside unit circle) and short-dwell transients. Practical radars cascade a wide-notch FIR (bulk weather/sea) with a narrow-notch IIR (residual ground clutter).
6.3 Adaptive and Clutter-Locked MTI
Clutter locking: estimate the mean clutter Doppler and recenter the notch on it — vital when the platform moves (shipborne, airborne).
TACCAR (Time-Averaged-Clutter Coherent Airborne Radar): the analog ancestor of clutter locking; a phase shifter driven by a feedback loop centres the clutter spectrum at 0 Hz.
Fully adaptive MTI: per-cell weights from the sample covariance matrix — the scalar-time ancestor of STAP (Section 9).
6.4 Blind Phases and I/Q Processing
Single-channel (envelope) differencing can cancel a target whose two successive echo magnitudes coincide — "blind phase." Coherent I/Q samples preserve phase, so cancellation depends only on fd; blind phases vanish. This is why all modern MTI is implemented at complex baseband.
Check Your Understanding 6
Give one advantage and one risk of an IIR MTI filter.
Why is clutter locking essential for an airborne MTI radar flying at 200 m/s toward its coverage?
Answers
Advantage: very narrow, deep notch with few pulses. Risk: instability/transient ringing, and poor behaviour in short dwells.
Own velocity imparts a bulk Doppler to all ground clutter (up to ~2v/λ ≈ 27 kHz at λ=0.03 m); without recentering, the notch sits at 0 Hz while clutter sits at 27 kHz → no cancellation.
Section 7 — Performance Measures
7.1 Definitions
Clutter attenuation: CA = Cin/Cout
MTI improvement factor: I = (Sout/Cout)/(Sin/Cin) = CA × (Sout/Sin)
Subclutter visibility (SCV): target-to-clutter ratio (dB) at which detection is still achieved — practically SCV ≈ I − (3–6 dB losses)
7.2 Improvement Factor of the Canceller (Gaussian clutter)
With clutter rms spectral width σc, normalised σn = σc/fr:
The σₙ⁻⁴ scaling of the double canceller explains its dominance: for narrow clutter, each canceller order multiplies rejection dramatically.
Example 7.1 — σc = 2 Hz (wind-blown trees), fr = 1000 Hz → σₙ = 0.002.
I₂ ≈ 2/(3 × 1.6×10⁻¹¹) ≈ 4×10¹⁰ ≈ 106 dB ideal; real systems reach only 45–60 dB — the gap is set by hardware instabilities (Section 9), not filter theory.
7.3 What Limits I
Limiter
Mechanism
Mitigation
Transmitter phase noise
Clutter echo decorrelates pulse-to-pulse
Low-noise STALO, stable amplifier chains
Internal clutter motion
Spectrum widens beyond notch
Higher-order cancellers, notch widening
Platform motion
Clutter shifted off 0 Hz
TACCAR / clutter lock / STAP
Quantization
Finite I/Q bits
≥ 12–14 bit ADCs, dither
Stagger misregistration
Sample timing error
Digital interpolation
Check Your Understanding 7
A radar achieves CA = 55 dB but the target loses 5 dB in the filter. What is I?
Why does SCV always fall short of I?
Answers
I = CA × Sout/Sin → 55 dB − 5 dB = 50 dB.
Detection needs margin above threshold after CFAR and integration losses; SCV ≈ I minus those losses.
Section 8 — Radar Clutter and Its Models
8.1 Clutter Taxonomy
Type
Physics
σ°/η range
Spectral width
Land (surface)
Soil, rock, vegetation, buildings
−40…−5 dB (grazing-angle & freq. dependent)
≈ 0–2 Hz (wind)
Sea (surface)
Wind-driven capillary/gravity waves
−50…+5 dB (sea-state dependent)
1–100 Hz
Weather (volume)
Rain, snow; also chaff
η ≈ 10⁻⁵–10⁻³ m²/m³
Tens of Hz (m/s spread)
Biological
Birds, insects ("angels")
Highly variable
Target-like
8.2 Amplitude Statistics
Model
PDF character
Best suited for
Rayleigh
One-parameter, light tail
Homogeneous land clutter, many scatterers
Log-normal
Heavy tail, log-domain Gaussian
Spiky low-grazing-angle terrain
Weibull
Two-parameter shape flexibility
Ground & sea clutter broadly
K-distribution
Compound (speckle × texture)
High-resolution sea clutter
The double canceller's f⁴ notch (red) rejects far more of the finite-width Gaussian clutter spectrum than the single canceller's f² notch (blue) — the quantitative reason two-pulse cancelling dominates practice.
Check Your Understanding 8
Which clutter type sets the widest notch requirement for a coastal surveillance radar?
Why do heavy-tailed models force higher CFAR thresholds?
Answers
Sea clutter (moving waves) plus rain — broad spectra, demanding wide notches and higher-order cancellation.
Large-amplitude clutter spikes occur more often; without a raised threshold these spikes cross the detector and flood the display with false alarms.
Section 9 — Limitations and the Path to MTD
9.1 Limitation Summary
Blind speeds (mitigated, never eliminated, by stagger).
Scalloping — gain varies with fd, causing fluctuating target returns.
9.2 From MTI to MTD and STAP
Generation
Key addition
Benefit
MTI (this unit)
Pulse-to-pulse cancellation
Clutter rejection, detection
MTD (Moving Target Detector)
FFT Doppler filter bank + per-bin CFAR + clutter maps
True velocity, better SCV, constant false alarms
STAP
Joint space (array) × time (Doppler) adaptation
Airborne early-warning: detect targets in mainlobe clutter
Exam tip: when asked to "compare MTI and pulse-Doppler", the crisp answer is: MTI = notch filter, binary moving/not-moving decision; pulse-Doppler = filter bank, unambiguous velocity measurement at the cost of PRF/ambiguity management.
Section 10 — Tutorial Problems
Attempt all problems before consulting the solutions. Values: c = 3×10⁸ m/s.
An L-band MTI radar (f₀ = 1.25 GHz) uses PRF = 900 Hz. Find (a) λ, (b) the first three blind speeds, (c) the unambiguous range.
A single delay-line canceller is used with fr = 1 kHz. Compute |H| at fd = 100, 250 and 500 Hz. Repeat for the double canceller.
Show from first principles that v₁Ru = λc/4. Hence choose fr for λ = 0.15 m so that the first blind speed exceeds 100 m/s, and state the resulting Ru.
A two-PRF stagger uses f₁ = 950 Hz and f₂ = 1050 Hz at λ = 0.1 m. Compute the first blind speed and the stagger ratio.
Ground clutter has σc = 0.5 Hz and fr = 800 Hz. Estimate the ideal improvement factors (dB) of single and double cancellers.
Define subclutter visibility. A system quotes I = 48 dB and total processing losses of 6 dB — what SCV should the data sheet claim?
Explain why an airborne early-warning radar flying at 250 m/s cannot use a fixed zero-Doppler notch at X-band, and name two techniques that rescue it.
A coastal radar must reject both sea clutter (σc = 20 Hz) and rain (σc = 40 Hz) at fr = 1200 Hz. Sketch (describe) a suitable cascaded filter architecture.
Given I/Q samples x(n) = {3+j2, −1+j4, −2−j3, 1−j1}, compute the output of a double canceller for n = 2 and n = 3.
Write a one-page comparison of coherent and non-coherent MTI with a table of four differences and two application scenarios for each.
Selected solutions
λ = 0.24 m; v₁ = 108 m/s; blind speeds 108, 216, 324 m/s; Ru = 166.7 km.
Own velocity shifts ground clutter by up to 2v/λ ≈ 16.7 kHz (λ = 0.03 m); fixed notch at 0 Hz misses it → clutter breakthrough. Remedies: TACCAR/clutter-locked notch; STAP.
Wide-notch binomial FIR (N = 3) first stage to suppress the broad sea/rain spectra, then narrow-notch IIR second stage for residual ground clutter; adaptive per-cell gains from measured covariance.
Space-Time Adaptive Processing — joint adaptive spatial and Doppler filtering for airborne radars.
Section 12 — References and Further Reading
M. I. Skolnik, Introduction to Radar Systems, 3rd ed., McGraw-Hill, 2001 — Ch. 4 (MTI and pulse-Doppler radar).
M. I. Skolnik (ed.), Radar Handbook, 3rd ed., McGraw-Hill, 2008 — MTI, clutter, and Doppler processing chapters.
N. Levanon, Radar Principles, Wiley-Interscience, 1988.
B. R. Mahafza, Radar Systems Analysis and Design Using MATLAB, 3rd ed., CRC Press, 2013 — canceller simulation examples.
F. E. Nathanson, Radar Design Principles, 2nd ed., McGraw-Hill, 1991 — clutter statistics and MTI hardware limits.
Academic honesty note: equation forms for improvement factor differ slightly between references depending on the clutter spectral model assumed. In examinations, state the model you use (e.g., Gaussian spectrum of width σc).