Introduction to Systems

Simon Wagner

Fraunhofer FHR Cognitive Radar Department

CoSIP Workshop Berlin, 07.12.2016

© Fraunhofer FHR What is Radar?

 Radar is a palindrome  Indicates the basic send/echo idea  Electromagnetic wave is transmitted and reflected by the target  Echo is received after a time proportional to the distance of the target

from www.radartutorial.eu

© Fraunhofer FHR RADIO DETECTION AND RANGING (RADAR)

 Is an active system (apart from passive )  A microwave radiation is emitted by the radar. The reflected radiation, known as the echo, is backscattered from the surface and received some time later.  Passive radars use foreign sources of illumination (e.g. FM radio, DVB-T, foreign remote sensing satellites …)  Independent on the day light, operation at day and night  Operates at microwave frequencies  Usual wavelength between 10 m – 0.1 mm  Clouds, fog, smoke, dust and other materials can be penetrated

© Fraunhofer FHR

RADAR Two basic concepts

Radar sensor for object detection Imaging radar and positioning  Position measurements over time  Generation of a quasi optical image allow target tracking (SAR, ISAR)  The resolution cells (range,  Resolution cells much smaller than direction, Doppler) are greater or target dimension equal to the object dimensions  Classification with range profiles,  Classification by signal strength radar images (RCS), Doppler , polarisation, dynamics of motion…

© Fraunhofer FHR Radar Transmitter and Receiver (and other wireless transmission systems)

 Signals that appear inside a radar are real valued voltages  Only real valued signals can be transmitted as electromagnetic waveform Why is radar signal processing normally complex valued?

 Due to physical reason (antenna size, transmission medium), the waveform is mixed with a carrier signal with frequency f0

 The spectrum of the radar signal is shifted to f0

Picture from Ohm, J.R.; Lüke, H.D.; Signalübertragung, Springer-Verlag, Berlin, 2010 – in German

© Fraunhofer FHR Radar Transmitter and Receiver (and other wireless transmission systems)

 If the transmitted signal is real, its spectrum is Hermitian

∗ yields ∗ 푠푅퐹 푡 = 푠푅퐹 푡 푆푅퐹 푓 = 푆푅퐹 −푓  All information is contained in one half of the spectrum  How to create a low pass signal with the same information content?

Pictures from Ohm, J.R.; Lüke, H.D.; Signalübertragung, Springer-Verlag, Berlin, 2010 – in German

© Fraunhofer FHR Radar Transmitter and Receiver (and other wireless transmission systems)

 Equivalent low pass signal s 푡 is complex valued  Transmitted signal can be describes as real part of the product of equivalent low pass signal and a complex carrier (amplitude and phase modulation)

푗2휋푓0푡 푗휙(푡) 푠푅퐹 푡 = 푅푒 푠 푡 푒 = 푅푒 퐴(푡)푒  Real and imaginary parts of low pass signal are separated in two channels  Real part is the In-phase (I)-channel  Imaginary part is the Quadrature (Q)-channel

Picture from Ohm, J.R.; Lüke, H.D.; Signalübertragung, Springer-Verlag, Berlin, 2010 – in German © Fraunhofer FHR

COHERENT RADAR Quadrature demodulator The QDM transfers a real valued RF signal to a complex baseband signal and recovers the complex envelope.

Reference frequency (RF)

Complex envelope, base band signal

Real valued RF signal

© Fraunhofer FHR COHERENT RADAR Quadrature modulator

The QM transfers the complex baseband signal to a real valued RF signal.

© Fraunhofer FHR COHERENT RADAR Generic radar system

* Point target r Distance to a point scatterer

r c0 Velocity of light t Traveling time Antenna N(t) White noise s(t;r) Received waveform T/R Switch a complex amplitude

QM QDM  Traveling time f0 as(t;r)  Traveling distance s(t) N(t)  Received signal Z(t) Figure: Radar system with baseband signals

© Fraunhofer FHR COHERENT RADAR Received waveform

 Wave length and wave number:

c0 0  f0 Complex envelope 2 k  0   Received waveform 0

R 2f0t  2f0 c0 2R  0

 k0 R

© Fraunhofer FHR COHERENT RADAR QDM

Optimum receive filter f0 as(t;r) N(t)

Y(t) h(t) Z(t)

i.e.

© Fraunhofer FHR COHERENT RADAR

The pulse response of the optimum filter is equal to the time-inverted, complex conjugated signal

The maximum SNR is . This filter is called matched filter.

Received signal

Replica

Response of matched filter

© Fraunhofer FHR COHERENT RADAR Point spread function

© Fraunhofer FHR COHERENT RADAR Matched filter, point spread function

Matched filtering means correlation with the transmit signal. The point spread function is the reaction of the receive filter to the transmit signal. The point spread function is equal to the of the transmit signal, if a matched filter is used. In this case it is the Fourier back transform of the magnitude-squared of the signal spectrum.

Point spread function

Reflectivity of three point targets Output of the matched filter

© Fraunhofer FHR COHERENT RADAR Coherent radar Pulse repetition frequency: PRF (~ 100 Hz ... 10 kHz)

Intrapulse sampling frequency: fs (~ 10 MHz ... 1 GHz) 1 F  PRF  s T 1 f  s t

T= 1ms: Covered range =150 km t= 1ns: Range sampling = 15 cm Figure: Two time scales for pulse radar

© Fraunhofer FHR Doppler Effect

 Up to now, the radar target was a non-moving isotropic point scatterer  Now it becomes a moving isotropic point scatterer  Movement of target implies Doppler shift on radar signal proportional to the targets speed

푠(푡)푒푗2휋푓0푡 푠(푡)푒푗2휋 푓0+푓퐷 푡

 Output of matched filter is no longer the autocorrelation of the transmit signal  Cross correlation between transmitted signal and Doppler shifted receive signal

∗ ∗ 푗2휋푓퐷휏 푝 푡 = 푠 휏 푠 푡 + 휏 푑휏 푝 푡, 푓퐷 = 푠 휏 푠 (푡 + 휏)푒 푑휏

© Fraunhofer FHR

Doppler Effect –

 The accuracy with which target and Doppler (velocity) can be estimated is given by the ambiguity function

∗ 푗2휋푓퐷푡 휒 휏, 푓퐷 = 푠 푡 푠(푡 + 휏)푒 푑푡

 Basic properties of the ambiguity function  Conservation of ambiguity volume

Volume under 휒 휏, 푓퐷 depends only on signal energy, not on the shape of the waveform Radar uncertainty principle: choosing a waveform that narrows the surface in one dimension will cause it to widen in the other dimension

© Fraunhofer FHR Doppler Effect – Ambiguity Function

 To determine the range resolution, the frequency domain expression of the ambiguity function is needed

∗ 푗2휋푓휏 휒 휏, 푓퐷 = 푆 푓 푆(푓 + 푓퐷)푒 푑푓

 For the range resolution a stationary target (푓퐷 = 0) is considered

휒 휏, 0 = 푆(푓) 2푒푗2휋푓휏푑푓

 Theoretical optimal range resolution is obtained by a Dirac delta function  Such a signal would have infinite energy  Functions broadly supported in the Fourier domain improve the resolution

© Fraunhofer FHR COHERENT RADAR

The solution is to expand the bandwidth by modulation of the pulse. The Rayleigh range resolution of a waveform with a rectangular spectrum S(f) of bandwidth b is given by

without direct dependence on the pulse length.

© Fraunhofer FHR Pulse compression – waveform

 Achieve high bandwidth with pulse modulation

푇푃 푡 − 휇 2 2 푗2휋 푓0푡+ 푡 푠푅퐹 푡 = 푅푒 푟푒푐푡 푒 2 푇푃

 Frequency modulation over time is given by 1 푑휙(푡) 푓 푡 = = 푓 + 휇푡 2휋 푑푡 0  휇 – sweep rate

 푇푃 – Pulse length

 푏 ≈ μ푇푃

 Time-bandwidth product = 푇푃푏

© Fraunhofer FHR Pulse compression – Chirp waveform

Compression Chirp result

tac t t

F F-1

t|.|2 f f Power Spectrum Spectrum

© Fraunhofer FHR Pulse compression – Chirp waveform

Figure: Ambiguity function of a chirp with Gaussian envelope

© Fraunhofer FHR SYNTHETIC APERTURE RADAR Imaging radar (SAR, ISAR) is based on ..

Measurement of range change via phase (-> cross range resolution)

 Relative aspect change of the scene, the object necessary  SAR: via motion of the platform  ISAR: via motion of the target

Measurement of range (pulse compression)

© Fraunhofer FHR SYNTHETIC APERTURE RADAR Applications: Analysis of floods

SAR-image of the Elbe river near Dömitz at times of a disastrous flooding. Left: 29 August 2002, Right: 22 October 2003.  FHR

© Fraunhofer FHR Radar imaging with turn table Range histories of r isolated scatterers



Radar T

360° Beam

Turn table

© Fraunhofer FHR Radar imaging with turn table

 How to achieve resolution perpendicular to the line of sight (LOS) of the radar?  Three point scatterer rotate clockwise with angular rate 휔 and rotation angle 휃

 Range of scatterer three (푅0 ≫ 푑3)

푟3 ≈ 푅0 − 푑3 sin 휃3 = 푅0 − 푑3 sin 휔푡푠푙표푤  Received signal is replica of transmitted signal delayed by travelling time τ = 2푟/푐

2푟3 푍 푡 ~푠 푡 − 휏 = 푅푒 푒푗2휋푓0 푡 − 푐

2푟 (푡 ) 푗2휋 푓 푡 − 3 푠푙표푤 = 푅푒 푒 0 푓푎푠푡 휆 with 푐 = 휆푓

2푅 2푑 sin 휔푡 푗2휋 푓 푡 − 0 + 푠푙표푤 = 푅푒 푒 0 푓푎푠푡 휆 휆

More details in Mensa, Dean L.; High Resolution Radar Cross-Section Imaging, Artech House, Boston, 1991

© Fraunhofer FHR Radar imaging with turn table

 Modulation of the signal from pulse to pulse

1 푑휙(푡푓푎푠푡, 푡푠푙표푤) 2휔푑 2휔푥 푓 푡푠푙표푤 = = cos 휔푡푠푙표푤 = 2휋 푑푡푠푙표푤 휆 휆

 The modulation frequency is proportional to the cross range position  Resolution in cross range is achieved via observations over time  The spectra of the scatterers of previous example indicate position

Scatterer 1 Scatterer 2 Scatterer 3

More details in Mensa, Dean L.; High Resolution Radar Cross-Section Imaging, Artech House, Boston, 1991

© Fraunhofer FHR Radar imaging with turn table

Sample field with compressed ISAR image after cross range range profiles FFT

© Fraunhofer FHR Questions…

or Ideas?

© Fraunhofer FHR