Design and Implementation of FMCW Surveillance Rad
Design and Implementation of FMCW Surveillance Rad
6, 2018
1 Abstract—We propose and implement a low-complexity architecture for acquisition and a FPGA-based processing
dual-chirp FMCW radar system for surveillance applications. platform of a radar sensor with a single transmitter and
The FMCW radar is suitable for the detection of various multiple receivers. However, these algorithms usually
positions of multiple targets for the monitoring of concealed
perform full-dimensional FFT when detecting in order to
humans. For a surveillance FMCW radar system, it is necessary
to distinguish between stationary and moving targets while distinguish between stationary and moving targets and thus
maintaining a low level of complexity. However, conventional require high complexity levels. The most important issue in
FMCW radar systems are very complex with many chirps to relation to surveillance radar is a quick identification of the
distinguish between stationary and moving targets. Even in presence of a moving target. Therefore, these algorithms are
hardware with low complexity levels, in order to distinguish not suitable for surveillance radar systems due to their high
between a stationary and a moving target, the proposed
complexity.
algorithm employs only dual chirps. Experimental results show
that the proposed algorithm can effectively distinguish between In this paper, we propose and implement a low-cost
moving and stationary targets despite its low complexity and FMCW radar algorithm with very low complexity for
low-level hardware. surveillance applications. The proposed radar algorithm
recognizes a moving target with an inexpensive chip and a
Index Terms—FMCW; Radar; Surveillance; Low simple RF. In order to distinguish between stationary and
complexity. moving targets, we initially employ only dual chirps. Due to
the characteristics of FMCW radar signals, when the target
I. INTRODUCTION moves, the phase according to the chirp changes due to the
Recently, there have been several studies of radar sensors Doppler effect. That is, the proposed algorithm recognizes
for surveillance applications [1]–[10]. These areas are the existence of the target by observing the phase change
popular because radar sensors have many advantages between the dual chirp signals. Furthermore, this system is
compared to other sensors. For example, radar sensors can implemented as a 24 GHz module with low complexity in
detect at longer distances than an ultrasonic sensor with a order to verify that the proposed algorithm is suitable for
short detection distance and can detect targets even in surveillance applications.
environments in which a camera sensor cannot operate easily, The structure of the paper is as follows. In Section II, we
such as at night or in fog [1]–[3]. introduce and define the system model and the main notations
Meanwhile, in the various radar systems used for used in order to ensure a proper understanding of the FMCW
surveillance applications, the frequency-modulated algorithm, such as the TX and RX signals used in FMCW
continuous wave (FMCW) radar system represents a radar systems. Section III shows the conventional FMCW
promising technology due to advantages such as lower radar detection algorithm using the 3D FFT. In Section IV,
complexity and cost levels [6]–[15]. Compared to pulse radar we present the structure of the proposed algorithm. In
systems, FMCW radar systems have relatively low Sections V and VI, simulation and experimental results for
transmitted peak power levels and can meet certain range and various cases are provided. Finally, Section VII concludes
velocity resolution requirements with relatively simple this paper.
hardware and architectures [6].
In earlier work [10]–[12], FMCW radar algorithms with II. SYSTEM MODELS
multiple channels for surveillance were designed and This section discusses the different models of FMCW
implemented. In one study [10], the authors investigated the radar systems and the notations used in them. We consider
implementation of 24 GHz FMCW radar for surveillance and reflected signals of FMCW radar in order to detect the
detected the range, velocity and angle of targets. Meanwhile, ranges, the velocities, and the azimuth distributions of targets
in other words [11]–[13], researchers presented a scalable [6]. To facilitate an easier understanding of the concepts, we
initially consider the case of a single target, after which we
Manuscript received 14 December, 2017; accepted 10 August, 2018. extend this to the case of multiple targets.
This work was supported in DGIST was supported in part by the DGIST Consider a uniform linear array (ULA) of K receiving
R&D Program of the Ministry of Science, ICT and Future Planning, Korea
(18-IT-01).
(RX) antenna arrays with distance d between adjacent arrays
60
ELEKTRONIKA IR ELEKTROTECHNIKA, ISSN 1392-1215, VOL. 24, NO. 6, 2018
x0 (t ) exp j 2 f 0t t 2 / 2 , (1)
where a0 is the complex amplitude, w l ,k is the additive white
where 0 t T , f0 is the carrier frequency, T is the Gaussian noise (AWGN) signal, fD is the Doppler frequency
frequency modulation period during which the TX frequency due to the moving of the target and is the wavelength. By
sweeps over the entire bandwidth B, and is the chirp rate, denoting the velocity term as vl f D exp j 2 lf DT and
i.e., B / T . All L chirp signals are transmitted during a
the DOA term as k exp j 2 dk sin / , xl ,k (t )
frame; i.e., TF = TL. The signal reflected from the target is
received with delay and direction of arrival (DOA) at can simply be represented as follows
the k-th array, as shown in Fig. 2. Here, the l-th chirp slot is
denoted by xl,k(t) and is expressed as xl ,k (t ) a0 x0 t vl f D k w l ,k (t ). (3)
a a0 exp j 2 f 0 / 2 , 2
(5)
Yl ,k u yl ,k [n]WNqn ,
n 0
RFFT
(8)
61
ELEKTRONIKA IR ELEKTROTECHNIKA, ISSN 1392-1215, VOL. 24, NO. 6, 2018
operation; i.e., W exp j 2 / N . That is, the range the range and velocity map with N RFFT N CFFT . Finally,
domain NRFFT-point FFT output vector at the l-th chirp and with an NAFFT-point FFT operation on the range and velocity
the k-th array Yl ,k Yl ,k [1], Yl ,k [2],..., Yl ,k [ N RFFT ] is map with N RFFT N CFFT in the array domain, the 3D map
which reflects the range, velocity and DOA is obtained.
obtained. Subsequently, the N CFFT point FFT operation on
However, in order to perform a 3D FFT, very high
Yl , k in the chirp domain is conducted where N CFFT is the complexity is required, as is well known. This causes the
number of points for the FFT in the chirp domain. This gives processing time and cost to increase.
Meanwhile, in order to reduce the computational two beat signals among all L beat signals according to the
complexity of multi-dimensional FFT operations, FFT velocity of interesting target; that is, the p-th and the q-th beat
algorithms with reduced dimensions have been proposed [16]. signals, yp,k[n] and yq,k[n], are selected. We then subtract the
In these algorithms, the range of the FFT output is initially two signals, with this action denoted by d pq,k [n] , as follows
determined, after which the peak of the magnitude of the FFT
output is detected. The authors use the characteristic by
d pq,k [n] y p ,k [n] yq ,k [n]
which the targets exist only in the range bin with the peak of
the magnitude of the FFT output. They conduct a
chirp-domain FFT operation on only the range bins in which
axb [nts ] k p q w p,k [n] wq,k [n]. (9)
62
ELEKTRONIKA IR ELEKTROTECHNIKA, ISSN 1392-1215, VOL. 24, NO. 6, 2018
Here, we properly select two chirps, i.e., p and q. This is In order to estimate the range of the target, we perform
done because we want to avoid a case in which the proposed FFT on d pq ,k [n] with respect to the sample domain n, as in
algorithm does not work properly. For example, if the phase
(8); thus, the NRFFT-FFT output of d pq,k [n] denoted by
changes at the two chirps are symmetrical, the output after
subtraction may be zero, as if there is no moving target. D pq ,k [u ] is obtained. We then obtain the index ipeak
However, the duration of one frame is relatively short corresponding to the beat frequency of d pq ,k [n] through the
compared to the moving time of the target. By properly
selecting two chirps, we can avoid cases in which the phase peak detection of | D pq,k [u ] | for 0 u N RFFT 1 . Figure
changes of two chirps are symmetrical. 6 shows an example of the result of the peak detection of
| D pq ,k [u ] | when ipeak 35 . From the obtained ipeak , the
estimated range R̂ is calculated as [16]
c ˆ c fˆb c ipeak f s
Rˆ , (11)
2 2 2 N RFFT
0.8
0.6
|Dpq,k[u ]|
0.4
ipeak =35
0.2
0
Fig. 4. The structure of range detection of the proposed algorithm. 15 20 30 2535 40 45 50 55
sample index, n
Figure 5 shows the waveforms of beat signals and their Fig. 6. Example of the result of peak detection of | D pq,k [u ] | .
difference according to a sample index for stationary and
moving targets, respectively, when l = 1 and 3. As shown in Furthermore, we use only two RX antennas in order to
Fig. 5(a), in the case of the stationary target, it can be seen ensure a low cost and low complexity. The proposed
that the difference between the two beat signals becomes algorithm does not perform array-domain FFT in order to
almost zero. In the other hand, in the case of the moving detect the DOA information. In the proposed algorithm, the
target, the difference between the two beat signals follows a DOA is estimated using the phase change between two arrays
sinusoidal waveform with the beat frequency as shown in instead of an array-domain FFT, as follows [17]
Fig. 5(b).
ˆ sin 1 Y p,1[ipeak ]Y p,2 [ipeak ]
1 y1,1 [n ]
sin 1 exp j 2 d sin / , (12)
Amplitude
y3,1 [n ]
0
d 13,1 [n ] where denotes the term apart from except DOA term in
-1
(3). Of course, this simple DOA estimation is only possible in
20 40 60 80 100 120 140 160 180 200 the condition with a single target. However, even with
sample index, n
(a)
multiple targets, if the ranges or velocities of the targets are
not identical, the outcome will be similar to that of a
1 single-target condition in a FMCW radar system [10].
Amplitude
y1,1 [n ]
0
y3,1 [n ] V. SIMULATION RESULTS
-1 d 13,1 [n ]
This section addresses the simulation results in order to
20 40 60 80 100 120 140 160 180 200
observe the feasibility of the proposed algorithm. The
sample index, n parameters for the simulation of FMCW radar system are
(b)
shown in Table I. In order to consider the multiple targets
Fig. 5. Difference between beat signals with l = 1 and 3 according to sample
index in the cases of stationary and moving targets: (a) stationary target; (b) condition, the number of targets is set to four and the 4th
moving target. target is set to stationary.
63
ELEKTRONIKA IR ELEKTROTECHNIKA, ISSN 1392-1215, VOL. 24, NO. 6, 2018
0.5 Target #3
0.25
0
0 10 20 30 40 50
Range (m)
Fig. 7. Snapshot of the result of detection of the proposed scheme.
1.5
estimated range and DOA of the proposed algorithm which is
1 experimented in the anechoic chamber. For this experiment,
out of three targets, two are stationary and only one target is
0.5 moving. The two stationary targets are located at 3.3 m and
6.2 m, respectively. The moving target moves between the
0 two stationary targets. As shown in Fig. 12, the power
0 5 10 15 20
SNR [dB] spectral density (PSD) of two stationary targets is
Fig. 8. RMSE results according to SNR. insignificant. On the other hand, for the moving target, the
64
ELEKTRONIKA IR ELEKTROTECHNIKA, ISSN 1392-1215, VOL. 24, NO. 6, 2018
peak of the PSD is located at approximately 5 m. From this slow mode. These results imply that the proposed algorithm
result, the proposed algorithm can detect only the moving can detect a moving target well in condition which contain
target. Fig. 13 shows that the proposed algorithm accurately clutter and stationary targets despite its low complexity.
detects the DOA of moving target.
Fig. 14. Photo of surveillance radar experiment for actual moving with
single target.
10
Fig. 11. Photo of back end modulation (BEM) [18].
range (m)
1
5
Moving target (5.2 m)
0.75 Stationary target (3.3 m)
Stationary target (6.2 m)
0
0 10 20 30 40 50 60
frame index
PSD
0.5
(a)
10
0.25
range (m)
5
0
0 5 10 15 20
Range (m)
Fig. 12. Experiment result of range detection of the proposed algorithm. 0
0 10 20 30 40 50 60
frame index
1 (b)
Fig. 15. Experiment results of range according to frame index (time) by our
0.8 developed surveillance radar: (a) slow mode; (b) fast mode.
65
ELEKTRONIKA IR ELEKTROTECHNIKA, ISSN 1392-1215, VOL. 24, NO. 6, 2018
the dashed and solid arrows show the movement of the target REFERENCES
1 and target 2, respectively. Figure 18 shows that according [1] M. Richards, M. A. Richards, J. A. Scheer, W. A. Holm, Principles of
to the elapsed time the movement of two targets is detected Modern Radar: Basic Principles. Institution of Engineering and
Technology. 2010. DOI: 10.1049/sbra021e.
well. It implies that the proposed algorithm can detect not
[2] J. Li, P. Stoica, “MIMO radar with collocated antennas”, IEEE Signal
only a moving target well but also even in multiple targets Process. Mag., vol. 24, no. 5, pp. 106–114, 2007. DOI:
under many clutter and stationary obstacle conditions. 10.1109/msp.2007.904812.
[3] H. H. Meinel, “Evolving automotive radar: From the very beginnings
1 into the future”, in Proc. European Conf. Antennas and Propagation,
Hague, Nethelands, 2014, pp. 3107–3114. DOI: 10.1109/eucap.2014.
Target 1 6902486.
0.75 [4] G. Reina, D. Johnson, J. Underwood, “Radar sensing for intelligent
vehicles in urban environments”, Sensors, vol. 15, no. 6,
pp. 14661–14678, 2015. DOI: 10.3390/s150614661.
[5] A. Natarajan, A. Komijani, X. Guan, A. Babakhani, A. Hajimiri, “A
PSD
0.5
77-GHz phased-array transceiver with on-chip antennas in silicon: [3]
Target 2 Transmitter and local LO-path phase shifting”, IEEE J. Solid-state
circuits, vol. 41, no. 12, pp. 2807–2819, 2006. DOI: 10.1109/jssc.
0.25 2006.884817.
[6] C. Hu, Y. Liu, H. Meng, “Randomized switched antenna array FMCW
radar for automotive applications”, IEEE Trans. on Vehi. Tech., vol. 63,
0 no. 8, 2014. DOI: 10.1109/tvt.2014.2308895.
0 5 10 15 20 25 30
[7] A. Tessmann, S. Kudszus, T. Feltgen, M. Riessle, C. Sklarczyk,
range (m) W. H. Haydl, “Compact single-chip W-band FMCW radar modules for
Fig. 17. Photo of surveillance radar experiment for actual moving with commercial high-resolution sensor applications”, IEEE Trans. on
multiple targets. Microwave Theory and Tech., vol. 50, no. 12, pp. 2995–2997, 2002.
DOI: 10.1109/tmtt.2002.805162.
10 [8] D. C. W. Low, K. W. Chang, R. Lin, E. W. Lin, H. Wang,
Target 1 Target 2 M. Biedenbender, G. S. Dow, B. R. Allen, “A single-chip W-band
8 transceiver with front-end switching receiver for FMCW radar
applications”, IEEE Microwave Millimeter-Wave Monolithic Circuits
Symp. Dig., pp. 225–228, 1995. DOI: 10.1109/mcs.1995.470952.
range (m)
66