Radar Doppler DSP
Radar Doppler DSP
ESIEE, Paris
September 1996
SPRA342
IMPORTANT NOTICE
Texas Instruments (TI) reserves the right to make changes to its products or to discontinue any
semiconductor product or service without notice, and advises its customers to obtain the latest version of
relevant information to verify, before placing orders, that the information being relied on is current.
TI warrants performance of its semiconductor products and related software to the specifications applicable
at the time of sale in accordance with TI’s standard warranty. Testing and other quality control techniques
are utilized to the extent TI deems necessary to support this warranty. Specific testing of all parameters of
each device is not necessarily performed, except those mandated by government requirements.
Certain application using semiconductor products may involve potential risks of death, personal injury, or
severe property or environmental damage (“Critical Applications”).
Inclusion of TI products in such applications is understood to be fully at the risk of the customer. Use of TI
products in such applications requires the written approval of an appropriate TI officer. Questions concerning
potential risk applications should be directed to TI through a local SC sales office.
In order to minimize risks associated with the customer’s applications, adequate design and operating
safeguards should be provided by the customer to minimize inherent or procedural hazards.
TI assumes no liability for applications assistance, customer product design, software performance, or
infringement of patents or services described herein. Nor does TI warrant or represent that any license,
either express or implied, is granted under any patent right, copyright, mask work right, or other intellectual
property right of TI covering or relating to any combination, machine, or process in which such
semiconductor products or services might be or are used.
Other brands and names are the property of their respective owners.
Contents
Abstract ........................................................................................................................... 7
Product Support ............................................................................................................. 9
World Wide Web......................................................................................................... 9
EE students are motivated by technical contests ..................................................... 10
A Tracking Radar is a good support for DSP teaching ............................................. 11
Axir = Automatic X-band Instrumentation Radar....................................................... 13
Aerial, RF head and pedestal ................................................................................... 14
Front-end processor for monopulse radar (FEPMR)................................................. 16
Tracker...................................................................................................................... 16
Radar Manager ......................................................................................................... 17
Simulator and Built-In-Test-Equipment ...................................................................... 18
AXIR applications ..................................................................................................... 18
Conclusion .................................................................................................................... 20
Figures
Figure 1. General architecture of AXIR............................................................................ 14
Figure 2. Air security........................................................................................................ 19
Figure 3. Weather hazards detection and warning .......................................................... 20
Digital Monopulse Doppler Radar and
DSP Teaching
Abstract
In 1995 several teams of EFREI students submitted a project to
the Texas Instruments DSP Challenge. This paper describes the
adventure of one project which passed with success the European
semi-finals to go and participate to the worldwide final in May,
1996 in ATLANTA.
The winning project deals with the Digital Signal Processing
(hardware and software) associated with a new class of a Digital
Monopulse Doppler Radar. The study has been an interesting
practical application of the Radar and DSP courses.
AXIR (Automatic X-Band Instrumentation Radar) is a very
innovative design. The analog RF head is reduced to the
minimum. The principal tasks are realized digitally by seven TMS
320 DSPs. Nearly all the circuits are located at the rear of the
printed planar monopulse antenna.
The digital processors are organized on four PCBs (Printed Circuit
Board)
1) Front-End Processor (Digital Pulse Compression and Doppler
Filtering)
2) Tracker ( Three Alpha-Beta evolutive filters smoothing the
cartesian coordinates)
3) Radar Manager (Automatic determination of the waveform and
of the coefficients)
4) Simulator and Built-In-Test-Equipment (BITE)
This low-cost tracking radar has many applications related to civil
aviation, meteorological surveillance or general instrumentation.
Product Support
Tracker
The Tracker is a second PCB comprising one DSP TMS320C40.
This design was studied by a second EFREI team and was
submitted to the TI challenge. The software tasks of this board
are:
1) to compute difference and sum channels in azimuth, elevation
and range,
2) to derive angular and range deviations as quotients,
Radar Manager
The Radar Manager is a third PCB comprising one DSP
TMS320C40. This board was studied by an other EFREI team to
be submitted to the TI challenge. This design is the most
intelligent part of the radar AXIR A software similar to an expert
system is in charge of piloting all the important parameters of the
radar.
q The class of the target is determined as a function of the
Signal To Noise Ratio, the range, the altitude and the velocity
along the trajectory.
q The waveform is regularly optimized, on a per burst basis. The
Digital Pulse Compression transmit sequence is selected as a
function of the Range and of the Radar Cross Section
estimated for the target. The Pulse Compressor code is
computed as a function of the Doppler velocity. The Pulse
Repetition Interval is also selected in a collection of pre-
defined values, via a cost-function optimization, avoiding
automatically Range Doppler ambiguities.
q The Doppler complex coefficients are derived first from the 3
Doppler filters outputs, during the acquisition process and the
preliminary tracking. Later during the confirmed tracking stage,
the Doppler filters are associated to the smoothed trajectory
radial velocity.
q The Alpha-Beta coefficients for the X, Y, Z tracker are selected
as a function of the target class.
q The operational sequencing, missions and constraints defined
by the operator are acting on the Radar Manager program.
AXIR applications
AM is a light and low-cost Instrumentation Radar designed fix civil
applications (no ECCM features). It is basically a Monopulse
Doppler Tracking Radar. The design of this sensor is very
innovative, using a minimum of components, affording thus a low
serial production cost and an excellent reliability.
Concerning the technical performances AXIR is able to track one
target up to 100 kin, while surveying other possible targets in the
beam and in 14 remaining range gates. The acquisition range
depends of the target RCS and of the a priori information about
the trajectory and the velocity. A typical value is around 20 km.
An optional Tracking TV equipment could be associated with the
Radar. allowing the operator to get a look on the tracked target,
depending of the range and of weather conditions. Data about the
trajectory (spherical and/or Cartesian coordinates. derivatives) are
available for Graphical or Text displays and for recording.
AXIR is capable of filling many jobs, principally in the area of civil
aviation (figures 2 and 3) :
q Short range air-control for small airfields, air-clubs or piloting
schools, public or private
q Air-security for parallel runways or difficult access airports
q Runways and taxiways ground control of one aircraft or
vehicle, with warning of around intruders
q Weather hazards detection and warning ( wind-shears, air
gusts, tornadoes ...) while tracking a plane
q Tracking of atmospheric sounders
q Trajectory control, monitoring and recording for aerobatics
sports events
q Cloud back-scattering analysis (via scanning)
Conclusion
AXIR has been an excellent application for a total of 16 EFREI
students, to realize practically hardware breadboards and
software. The four Processors were simulated and tested in High
Level Language. Each processor was tested independently due to
the limited time of this experience. EFREI is looking for a contract
to develop completely the pedestal, the RF head and to finish a
global integration.
EFRIE, France
December 1995
SPRA299
IMPORTANT NOTICE
Texas Instruments (TI) reserves the right to make changes to its products or to discontinue any
semiconductor product or service without notice, and advises its customers to obtain the latest version of
relevant information to verify, before placing orders, that the information being relied on is current.
TI warrants performance of its semiconductor products and related software to the specifications applicable
at the time of sale in accordance with TI’s standard warranty. Testing and other quality control techniques
are utilized to the extent TI deems necessary to support this warranty. Specific testing of all parameters of
each device is not necessarily performed, except those mandated by government requirements.
Certain application using semiconductor products may involve potential risks of death, personal injury, or
severe property or environmental damage (“Critical Applications”).
Inclusion of TI products in such applications is understood to be fully at the risk of the customer. Use of TI
products in such applications requires the written approval of an appropriate TI officer. Questions concerning
potential risk applications should be directed to TI through a local SC sales office.
In order to minimize risks associated with the customer’s applications, adequate design and operating
safeguards should be provided by the customer to minimize inherent or procedural hazards.
TI assumes no liability for applications assistance, customer product design, software performance, or
infringement of patents or services described herein. Nor does TI warrant or represent that any license,
either express or implied, is granted under any patent right, copyright, mask work right, or other intellectual
property right of TI covering or relating to any combination, machine, or process in which such
semiconductor products or services might be or are used.
Other brands and names are the property of their respective owners.
CONTACT INFORMATION
Abstract
FEPMR is standing for Front-End and Doppler Processing for
Monopulse Doppler Radar. It is a PCB at the frontier between the
RF and the Digital Circuits for the radar AXIR. There are four
channels, corresponding to the four squinted beams of the
antenna ( Right, Left, Up and Down ). Each beam channel is
processed in phase (I) and quadrature ( Q) real baseband signals.
The eight video analog signals are amplified and then converted
by 8 ADC ( Analog to Digital Converters ). Four Texas Instruments
(TI) TMS320C50 dgital signal processors (DSPs) ( one for each
beam channel ) are operating two major algorithms in cascade.
A flexible FIR ( Finite Impulse Response ) Digital Filter with
complex coefficients is first applied as a Digital Pulse
Compressor to provide 16 range bins. AXIR transmitter is
using a quadriphase sequence (1 up to 31 sub-pulses).
Coefficients for the reference code are optimised by the Radar
Manager ( project C ) to reduce the range side lobes.
The second part implements three IIR ( Infinite Impulse
Response ) Digital Filters, using also complex coefficients as a
Doppler filter. These coefficients are dynamically computed by
the Radar Manager ( project C ) as a function of the tracked
target range and Doppler. An output time decimation is
provided in accordance with the coherent and non-coherent
integration process.
For each beam, each range cell and each Doppler, the I and Q
signals are combined to form magnitudes to be delivered to the
TRACKER ( project B ) via 4 dual-port RAM. At the end of the
digital calculations, there are 4 * 16 * 3 bins, corresponding
respectively to the 4 beams, the 16 ranges and the 3 Doppler
filters. The output rate is fixed to 2 ms, independently of the
flexible Pulse Repetition Interval ( PRI ) controlled dynamically
by the Radar Manager.
This project is implemented on a Printed Circuit Board ( PCB ),
using four TMS320C50. In a close future, a single TMS320C82
could probably do the same job within the same real time
specifications because of its improvements.
This document was an entry in the 1995 DSP Solutions
Challenge, an annual contest organized by TI to encourage
students from around the world to find innovative ways to use
DSPs. For more information on the TI DSP Solutions Challenge,
see TI’s World Wide Web site at www.ti.com.
Introduction
This project deals with the study of a Front-End Processing for a
Monopulse Doppler Radar ( FEPMR ), in charge of Pulse
Compression and Doppler Digital Filtering. AXIR ( Automatic
Xband Instrumentation Radar ) is a concept of a low cost radar,
specially designed for the TEXAS INSTRUMENTS DSP
SOLUTIONS CHALLENGE. The idea was to use the
performances of the TMS320C4x and TMS320C5x families to
design a performant and intelligent radar processing, made of
multiprocessors DSP chips.
The objective product cost is under 80.000 $ for the complete
radar, including the pedestal, the antenna, the transmitter, the
R.F. Front-End processors and the display. AXIR is basically a
Monopulse Doppler Ground Based Radar with many civilian
applications:
Tracking of Meteo sounders.
Cloud Doppler analysis for Meteo purposes.
Short range air control for parallel runways or difficult access
airports.
Wind shear and burst detection.
Trajectory control for sportive aerobatics events or air clubs.
General low cost instrumentation radar ( radar cross section
evaluation, tutorial radar for universities, private air-
tracking…).
The radar processor has been divided into four separate projects
submitted to the TI challenge, respectively in charge of the
following tasks:
B. DIGITAL TRACKER.
C. RADAR MANAGER.
D. TARGET SIMULATOR.
The pedestal uses robotics low cost and modern technology. The
antenna box is free to rotate 360°, both in azimuth and elevation. A
standard PC, under MICROSOFT WINDOWS / MSDOS is used as
an operator remote control, graphics display, monitoring and
recording. A radio link could be managed between the sensor and
the operator desk. This project corresponds to the Printed Circuit
Board number 1 ( PCB #1 ).
Three undergraduates students from the Ecole FRancaise
d’Electronique et d’Informatique have elaborated the hardware
and the software of this project, under the supervision of an
advising Professor.
Video Amplification
There are four channels, corresponding to the four squinted
beams of the antenna ( Right, Left, Up and Down ). Each beam
channel is processed in phase ( I ) and quadrature ( Q ) real
baseband signals.
These eight signals are amplified by eight monolithic operational
video amplifiers NE 5539 ( Philips Semiconductors RF
Communications Products ). This model was selected for its wide
bandwidth and its high slew rate. They provide a true differential
input impedance device, with proper external compensation to
design operations over a wide range of closed-loop gains. The
amplifiers are providing 30 dB gain, a 2 MHz , bandwidth and a
(
low noise 4nV/ Hz . )
Figure 1. Frequency Response of the Amplifier
As the DSP is running with interruption IT1, the signal which will
activate the interruption during at least 3 consecutive machine
cycles have to be chosen. The input of this interruption is a low
active state and as the frequency of new samples is 2 Mhz, a
combination between the signal PRI and the 2 Mhz clock, which
are synchroned because issued from the same Master Clock (
project C ), can be used as IT1.
The 16 range cells given by the DPC algorithm are the elementary
informations taken from the received signal on a PRI period. They
will be associated with the Doppler filtering which will give the
information necessary to track the detected target.
After having simulated the algorithm in C language and in
comparison with the arrays of calculation, the reason for ordering
the data in memory and the importance of the shift implemented
by the LTD instruction can be seen. To optimise real time and
speed calculation, the for loop can be replaced by a shift of the
data in memory (LTD). This instruction loads the T register with
the contents of the address, adds the result of the previous
multiply to the accumulator and shifts the data to the next higher
address in data memory. Associated with MPY ( multiplication of
the contents of the T register with the contents of the address ),
the calculations are optimised for real time specifications.
ZAC ;Zero the accumulator
LTD VideoQ ;CoefR . Doppler
MPY LxQ
LTS Videol ;CoefQ . DopplerQ - CoefI . DopplerI
MPY Lxl
DMOV Videol
These instructions can be repeated nLx times with the RPTB
instruction, improving the time for calculations.
Doppler Processing
The Doppler Processing is used to reduce the bandwidth (
coherent integration ) and to reject fixed echoes. It is applied on
the 16 complex range bins calculated by the Pulse Compressor. It
corresponds to a complex coefficients BandPass Filter ( BPF ),
centered around a specific Doppler frequency.
This filter is used 3 times, with 3 different center frequencies, on
the 16 * 2 inputs ( real and imaginary parts ). So, the outputs
contained 16 * 2 * 3 values. These 3 filters are Digital order one
filters, based on a relationship between the input sequence
DPC_Result ( n ) and the output sequence Doppler ( n ).
Calculations
The Z transfer of the three Doppler filters is:
K
H (z) =
1-C.z-1
Each filter uses 2 coefficients : K and C.
K corresponds to a gain and is common to the 3 filters.
1- C
ENBW = where C is the complex pole.
1+ C
The calculations are done for the 16 range cells ( real and imaginary
parts ). The results give 6 arrays of 16 elements, loaded between
0400h and 045Fh in the data memory.
The initial conditions are set at the beginning of the tracking
phase:
1 1
ENBW = ⇒ C =
3 2
1
fcentral = (C0 )
2
flower = (C - 1)
1
6
fupper= (C1)
5
6
Figure 5. Poles of 3 Narrow Filters Around the Central Doppler Frequency During
the Tracking Phase (Example)
C Program
The program in C language is specified before writing the ASM
program and is given in the Appendix F. It can help to optimise the
ASM program by taking advantages ( for execution time and data
memory ) of the TMS C50 features. To maximise the real time
calculations, the same instructions than in the DPC program are
used : LTD, MPY, DMOV and RPTB.
Approximate Modulus
Magnitudes of previous results are calculated and transmitted to
the Digital Tracker.
The following approximate formula is used:
1
M= I + Q + I − Q
2
Architecture
The signals CLKR ( Receive Clock Signal ), DR ( Receive serial
Data signal ), FSR ( Receive Framing synchronisation Signal ) of
our DSP's are connected to the Radar Manager' DSP ( project C ).
Configuration
The use of the serial port receiver needs to program three
memory registers :
Hexadecimal address
SPC Serial Port Control register 22h
DRR Data Receive Register 20h
IMR Interrupt Mask Register 4h
IFR Interrupt Flag Register 6h
‘i’, ‘c’, ‘s’ :integers ( 16 bits ) Real part ( α ) for the inferior,
central and superior
frequencies.
Transmission Protocol
Two different methods have been studied:
The Radar Manager (project C) transmits only the variables
that need to be changed (for example : K, 1200, i, 1300, C,
2400). The main advantage of this method consists in the fact
that the Radar Manager doesn't need to transmit the eight
variables. But it needs a lot of cycles at the reception for the
tests to recognise the variables.
An other method consists to transmit the 8 words (16 bits).
When the receive buffer is full, the integer is placed into the
corresponding variable. This method has been retained
because of its speed : all the variables are sent every 2 ms.
The starting reception is detected by an interrupt signal. This
signal is generated when the first integer is received. Then the
procedure global_variables is executed.
Appendices
General description of the Front-End Processing. A
Global diagram of the 4 beams. B
Functional diagram for one antenna channel. C
Electrical diagram for one channel. D
Memory allocation. E
C program for simulation. F
C and ASM programs for the serial communication. J
All the addresses are given in hexadecimal and refer to the DATA
memory and not to the program memory.
{ int i;
switch ( nlx )
{ case 3 : LxQ [1] = 1; LxQ [2] = -1;
break;
case 5 : LxQ [1] = 1; LxQ [2] = 1; LxQ [3] = -1; LxQ [4] = 1;
break;
case7 : LxQ [1] = 1; LxQ [2] = 1; LxQ [3] = -1; LxQ [4] = -1 LxQ [5] = 1;
LxQ [6] = -1;
break;
case 11: LxQ [1] = 1; LxQ [2] = 1; LxQ [3] = -1 LxQ [4] = -1; LxQ [5] = -1;
LxQ [6] = 1; LxQ [7] = -1; LxQ [8] = -1; LxQ [9] = 1; LxQ [10] = -1;
break;
case l3: LxQ [1] = 1; LxQ [2] = 1; LxQ [3] = 1; LxQ [4] = 1; LxQ [5] = -1;
LxQ [6] = -1; LxQ [7] = 1; LxQ [8] = 1; LxQ [9] = -1; LxQ [10] = 1;
LxQ [11] = -1; LxQ [12] = 1;
break;
case 15: LxQ [1] = 1; LxQ [2] = 1; LxQ [5] = 1; LxQ [9] = -1; LxQ [12] = 1;
LxQ [13] = -1; LxQ [14] = 1; LxQ [3] = 1; LxQ [4] = 1; LxQ [6] = -1;
LxQ [7] = -1; LxQ [8] = 1; LxQ [10] = -1; LxQ [11] = 1;
break;
case 2l: LxQ [1] = 1; LxQ [2] = -1; LxQ [3] = -1; LxQ [4] = -1; LxQ [5] = -1;
LxQ [6] = -1; LxQ [7] = -1; LxQ [8] = -1; LxQ [9] = 1; LxQ [10] = 1;
LxQ [11] = -1; LxQ [12] = -1; LxQ [13] = 1; LxQ [14] = -1; LxQ [15] = 1;
LxQ [16] = -1; LxQ [17] = 1; LxQ [18] = -1; LxQ [19] = -1; LxQ [20] = 1
break;
case 3l: LxQ [1] = 1; LxQ [3] = -1; LxQ [7] = 1; LxQ [8] = 1; LxQ [9] = 1;
LxQ [11] = 1; LxQ [14] = 1; LxQ [16] = 1; LxQ [19] = -1; LxQ [21] = -1;
LxQ [22] = 1; LxQ [23] = -1; LxQ [27] = 1; LxQ [29] = -1; LxQ [30] = 1;
LxQ [2] = 1; LxQ [4] = -1; LxQ [5] = -1; LxQ [6] = -1; LxQ [10] = -1;
LxQ [12] = 1; LxQ [13] = -1; LxQ [15] = 1; LxQ [17] = -1; LxQ [18] = -1;
LxQ [20] = 1; LxQ [24] = 1; LxQ [25] = -1; LxQ [26] = 1; LxQ [28] = -1;
}
}
}
}
{ int i;
{ int i;
clrscr ( );
for ( i =0 ; i < ntxt – 1 ; i++)
{ printf ( "%c", texte [i] );
}
printf ( "%c \n", texte [ntxt – 1 ] );
for ( i = 0 ; i <nvideo ; i++ )
{ printf ( "%2d : Valeur %4d %4d \n", i, LxQ [i], LxI [i] );
}
I = getch ( );
}
/* Calculation of Doppler coefficients from the frequencies and the amplitude delivered by the project C via the serial port.
=> normally designed by the projet C, but here for the simulation. */
void calcul_coef ( float fzero, float fmoins, float fplus, float amp, int coefQ [ ], int coefI [ ] )
{ int i;
/* Main program. */
int main ( )
doppler1Q [i] = K*DPC_ResultQ [i] + coefQ [1]*doppler1Q [i-1] - coefI [1]*doppler1l [I-1];
doppler1l [i] = K*DPC_ResultI [i] + coefI [1]*doppler1Q [i-1] + coefQ[1]*doppler1I [I-1];
doppler2Q [i] = K*DPC_ResultQ [i] + coefQ [2]*doppler2Q [i-1] - coefI [2]*doppler2I [i 1];
doppler2l [i] = K*DPC_ResultI [i] + coefI [2]*doppler2Q [i-1] + coefQ [2]*doppler2I [I-1];
}
affiche ( nDPC_Result, doppler0Q, doppler0I, 7, "doppler0" );
affiche ( nDPC_Result, doppler1Q, doppler1l, 7, "doppler1" );
affiche ( NPC_Result, doppler2Q, doppler2l, 7, "doppler2");
// Constant //
#define nb_word_transfered 8
#define begin_adr Ox9OOOh //9000h->9016h//
#define adr_serial_buffer 32
// Global variables //
L, K, I, c, s, I, C, S :integer;
// Reception of interruption Rint signales that the first byte is being transmitted //
// Execution of procedure void change_global_variables (void) //
// Local variables //
integer *buffer_DRR;
integer *register;
{
register=04h; // adress of the register IMR //
*register=0h; // inhibit interruptions //
register=022h; // adress of the register SPC II
buffer_DRR = adr_serial_buffer;
K= *buffer_DRR;
repeat until (*register & 04 00h == 1) // wait until bit RRDY==1; an integer is received //
L= *buffer_DRR;
repeat until (*register & 04 00h == 1) // wait until bit RRDY==1; an integer is received //
i = *buffer_DRR;
repeat until (*register & 04 00h == 1) // wait until bit RRDY==1; an integer is received //
c= *buffer_DRR;
repeat until (*register & 04 00h == 1) // wait until bit RRDY==1; an integer is received //
s= *buffer_DRR;
repeat until (*register & 04 00h == 1) // wait until bit RRDY==1; an integer is received //
I= *buffer_DRR;
repeat until (*register & 04 00h == 1) // wait until bit RRDY==1; an integer is received //
C= *buffer_DRR;
repeat until (*register & 04 00h == 1) // wait until bit RRDY==1; an integer is received //
S= *buffer_DRR;
register=04h; // adress of the register IMR //
*register = 0l0h; // allow RINT interruption II
}
CHANGE_GLOBAL_VARiABLES
LACL DRR ; Load received value
SACL *+ ; Write to memory block
TEST1 BIT *,5,SPC ; wait transition from 0 to 1 of bit 10 (RRDY) of SPC
BCND TEST1, NTC ; TC=0, return to TEST1
NEXT LACL DRR ; Load received value
SACL *+ ; Write to memory block
….
RETE ; Clears INTM (enable interrupt) and execute RETI,
; return from interrupt
Glossary
ADC : Analog to Digital Converter.
AXIR : Automatic Xband Instrumentation Radar.
CLKR : Receive CLocK
DMA : Double Memory Access.
DR : Receive serial Data signal
DRR : Data Receive Register
DSP : Digital Signal Processor.
FIR : Finite Impulse Response.
FSR : Receive Framing Synchronisation signal
IFR : Interrupt Flag Register
IIR : Infinite Impulse Response.
IMR : Interrupt Mask Register
MSB : Most Significant Bit
PRF : Pulse Repetition Frequency.
PRI : Pulse Repetition Interval.
R.F. : Radio Frequency.
SPC : Serial Port Control register