1

MIMO UWB Radar System with Compressive
Sensing

arXiv:1604.04291v1 [cs.IT] 14 Apr 2016

Xia Li, Student Member, IEEE, Zhen Hu, Robert C. Qiu, Fellow, IEEE,

Abstract—A multiple input multiple output ultra-wideband
cognitive radar based on compressive sensing is presented in this
letter. For traditional UWB radar, high sampling rate analog to
digital converter at the receiver is required to meet Shannon
theorem, which increases hardware complexity. In order to
bypass the bottleneck of ADC or further increase the radar
bandwidth using the latest wideband ADC, we propose to exploit
CS for signal reconstruction at the receiver of UWB radar for
the sparse targets in the surveillance area. Besides, the function
of narrowband interference cancellation is integrated into the
proposed MIMO UWB radar. The field demonstration proves
the feasibility and reliability of the proposed algorithm.

Target 2

Target 1

AWG

RX antenna

RX antenna
TX antenna

TX antenna

Index Terms—MIMO, UWB, radar, CS, NBI, testbed.
Transmitter
based on
FPGA

Trigger signal
Receiver

I. I NTRODUCTION
USB 2.0

In the radar system, the unprecedented radio bandwidth
provides advantages such as the high precision of range
estimation. However, the extremely high sampling rate of
analog to digital converter (ADC) becomes a major challenge
in the ultra-wideband (UWB) radar systems. According to
Nyquist sampling theorem, the sampling rate should be at least
twice the bandwidth of the signal. For better system performance, oversampling is often required and results even higher
sampling rate. Compressed sensing (CS) gives an opportunity
to overcome this challenge. CS [1] theory specifies a new
signal acquisition approach, in a high probability, allowing
the acquisition of signals at a much lower data rate than the
Nyquist sampling rate.
Some previous work of our group [2] applied CS to the
communication application. In the communication system, the
received waveforms are synchronized by some mechanism
like preamble. However, in the MIMO radar application, the
received signal cannot be considered as synchronized since
the arrival time of signal depends on the over the air delay,
which is determined by the target distance to the radar. The
unknown delay leads to the randomness of the measurement
matrix construction. A cyclic prefix (CP) is used in waveform
design to make sure the measurement matrix is constant and
unrelated to the time delay. Therefore, each receiver antenna
can construct the basis matrix locally. Based on the analysis
of the system, the architecture of the testbed is proposed. This
testbed is based on our previous results of the UWB time reversal experiment [3] and cognitive UWB radar experiment [4]–
[6]. Both of them are on the full sampling rate. However, the
new system satisfies the specialty of the compressed sensing
X. Li, Z. Hu, and R. C. Qiu are with the Department of Electrical and Computer Engineering, Center for Manufacturing Research, Tennessee Technological University, Cookeville, TN, 38505, e-mail: xli43@students.tntech.edu;
zhu@tntech.edu; rqiu@tntech.edu.

Gigabyte
Ethernet
MATLAB
computing
engine

Fig. 1.

Block diagram of the system architecture

structure. The sampling rate at receiver can be much lower
than the Nyquist rate used in the two previous systems.
The experiment results based on measured data are provided
to demonstrate the effectiveness of the algorithm. The performance curves of the algorithm under different NBI and SNR
conditions are provided in simulation.
II. BACKGROUND AND P ROBLEM F ORMULATION
Since the transmitter and receiver of the our system are
synchronized by a trigger signal as shown in Fig. 1, ToAs
between transmitter antennas and receiver antennas can be
estimated from the recovered channel information, which are
used to estimate the target positions [4].
In traditional channel estimation theory, the channel vector
h ∈ Cn can be estimated from the observations of the linear
convolution, Xh = y ∈ Cn ., where the observation matrix
X ∈ Cn×n is a partial Toeplitz matrix.
By using CS theory, a sparse, high-dimensional vector
h ∈ Cn can be accurately recovered from a small number of
linear observations, RΩ Xh = ȳ ∈ Cm . Consider an arbitrary
index set Ω ⊂ {0, 1, . . . , n − 1} whose cardinality |Ω| = m.
We define the operator RΩ ∈ Cm×n that restricts a vector to
the entries listed in Ω. In our testbed, the ȳ can be viewed as
the sub-sampling of y. A delay in ȳ means the change of RΩ ,
which means the measurement matrix RΩ X is different. Since
the delay is unknown to the receiver, channel information
cannot be recovered at receiver locally.

2

A cyclic extension of time length TG is used to change
the linear convolution to circulant convolution. If the channel
delay is within TG , the system can be described in matrix as
y = Xh + n, where X is a circulant matrix as


xn−1 xn−2 · · · x1
x0
 x0
xn−1 . . . x2
x1 


(1)
X= .
.
.
.
..  .
..
..
..
 ..
. 
xn−2 xn−3 · · · x0 xn−1
The delay in dictionary nd can be represented by the
channel latency. This is proved by (2). It holds when the delay
nd is smaller than the length of CP.
yc [n − nd ]

=

=

=

=

N
−1
X
m=0
N
−1
X

of the received waveform with full sampling rate.
 
 
 
h11
h21

 


h12
h22
ȳ1 ȳ2 = X1 X2 
| {z } | {z }
h1
h2


+ n̄1 n̄2
(5)
The corresponding partial circulant matrix is defined as X̄ =
RΩ X. The partial random circulant matrix is generated by a
Rademacher sequence. It is proved in [8], the sth restricted
isometry constant δs has
( 3
)
r
3
s
s2
E [δs ] 6 C1 max
log 2 n,
log s log n ,
(6)
m
m

h [m] x [hn − (nd + m)iN ]

where C1 > 0 is a universal constant.
(6) means that for a given δ ∈ (0, 1), we have E [δs ] 6 δ
provided by
o
n
3
3
(7)
m > C2 max δ −1 s 2 log 2 n, δ −2 slog2 nlog2 s ,

h [m0 − nd ] x [hn − m0 iN ]

where C2 > 0 is another universal constant.
For 0 6 λ 6 1 and C3 > 0, we can have

h [m] x [hn − nd − miN ]

m=0
N
−1
X
m0 =0
N
−1
X

h [m − nd ] x [hn − miN ]

P (δs > E [δs ] + λ) 6 e
(2)

m=0

The vector h is the sampled signal of the channel impulse
response h (t). h (t) P
can be treated as a time limited pulse train
of the form h (t) = αm δ (t − τm Ts ), where the amplitudes

2

m

If the delay τm is an integer, then all the energy from αm is
mapped to tap hτm . However, if τm is not an integer, its energy
will leak to all taps of hk . However, most of the energy is kept
in the neighborhood of the original pulse locations. When the
target number is small, the number of wireless paths is also
small. In this case, h can be considered as a sparse signal.
For a 2 × 2 MIMO system shown in Fig. 1, the system
equation can be described as
 
 
 
h11
h21

 


h12
h22
y1 y2 = X1 X2 
| {z } | {z }
h1
h2


+ n1 n2 ,
(4)
where y1 and y2 are the received waveforms from receiver
antennas 1 and 2. For the transmitted waveforms x1 and x2 ,
X1 and X2 are the respective Toeplitz matrices. n1 and n2
are the i.i.d. complex Gaussian variables for the two receiver
antennas. hji is the impulse response of the channel between
transmitter antenna i to receiver antenna j.
For the CS-based radar system, the sub-sampled received
waveform can be considered as the result of down-sampling

,

(8)

s
C3 m
log2 slog2 n.

where σ =
In radar systems, the number of multipath is equivalent to
the number of targets. The channel information is sparse to
the most cases. The channel information can be recovered by
using the following convex optimization program
min (kh1 k1 + kh2 k1 )

Xh1 − ȳ1
subject to
Xh2 − ȳ2

αm of channel taps are complex valued. The samples of h are
modeled as [7]
π
sin (πτm )
1 X
 . (3)
hk = √
αm e−j N (k+(N −1)τm )
π
n m
sin N
(τm − k)

−λ2
σ2

2
2

6 ε1 ,
6 ε2

(9)

where k·k1 is the `1 norm and k·k2 is the `2 norm. εj is
chosen such that it bounds the amount of noise in the measurement from jth receiver antenna. The convex optimization
problem (9) is a second-order cone problem and can be solved
efficiently by utilizing the Basis Pursuit Denoising (BPDN)
algorithm [9].
III. NBI CANCELLATION USING CS THEORY
The system is modified to achieve the robustness under
NBI. In (4), NBI signal is not taken into account in algorithm,
NBI is treated as noise. However, the sparsity of the NBI can
be used to recover its information from sub-sampled received
waveform.


ȳ1

ȳ2



=X



h1

h2



+



n̄1

n̄2



=X



h1

h2



+



n̄1

n̄2



+



j̄1

+



FH jf 1

j̄2



FH j f 2


(10)

where FH is a randomly chosen rows of the discrete Fourier
matrix and can be written as FH = RΩ FH . jj is a vector of
the NBI signal received at antenna j. jf j is a vector of the
NBI in frequency domain. A recent research [10] has proved

3

that this measurement matrix has small restricted isometry
constants with very high probability.
For the receiver antenna j, we can have


 hj1

ȳj = X1 X2
+ F−1 jf j + n̄j
hj2


i hj1
h
(11)
=
FH  hj2  +n̄j .
X1 X2
|
{z
}
jf j
| {z }
Ā
sj

Because both the channel information and the NBI interference information are sparse, the interference signal and
channel information can be recovered simultaneously by using
the following convex optimization program
min (ks1 k1 + ks2 k1 )

As1 − ȳ1
subject to
As2 − ȳ2

=



|

X1

X2

X1

X2
{z
X̄





hj1
hj2



2
2

6
6

,

(12)

+ FH jf j + n̄j − FH j̃fj

Local Oscillator 4.0 GHz

I
D/A Conveter
Dual- core
1Gsps/14- bit

Xilinx Virtex-5 FPGA

Q

Power
Amplifier

Quadrature Modulator

Block diagram of the system architecture
Quadrature Demodulator
400~6000 MHz

I

A/D Converter
DDR2
Memory
Low Noise RF Bandpass
Amplifier
filter

A/D Converter

Q

FPGA
Array

Local Oscillator
4.0 GHz

Host Computer
Xilinx
Chipscope Pro
Analyzer

I
DDR2
Memory

A/D Converter

Low Noise RF Bandpass
Amplifier
filter
Q
Quadrature Demodulator
400~6000 MHz

Fig. 3.

A/D Converter

4 Channels 3Gsps Digitizer with FPGA Array In NI PXIe Chassis

Block diagram of the system architecture

traditional system. Due to the CS algorithm, the sampling rate
of the received signal can be reduced to much lower than
the Nyquist rate. As an example, we set the receiver ADC
sampling rate 13 of the full signal bandwidth, which is 166.66
Msps, as shown in Fig. 3. The reduced sampling rate in the
receiver can greatly reduce the cost and power consumption of
the receiver. The analog RF down-converter is the same as the
traditional system. A bandpass filter is used to remove noise
outside of signal bandwidth.
V. E XPERIMENT R ESULTS AND P ERFORMANCE A NALYSIS



(13)

hj

min (kh1 k1 + kh2 k1 )

Xh1 − ȳ10
subject to
Xh2 − ȳ20
where

Power Splitter
1:2

Reference Clock 20 MHz
For Local Oscillator

Baseband
Processing

Fig. 2.

Power
Amplifier

Q

JTAG

The channel information can be recovered by using the
following convex optimization program

ε00j

1Gsps/8-bit
2 Cores both with I&Q

Host
computer

Quadrature Modulator



hj1
+w̄j + n̄j
} hj2
| {z }


Digital Arbitrary
Waveform Generators

Power Splitter
1:2

ε01
ε02

where s1 contains both the channel information and interference position in frequency domain. ε0j bounds the energy of
n̄j .
NBI signal is always sparse in frequency domain due to
its narrow band feature. The number of taps in channel
information equals the number of targets in the surveillance
area and suffers the influence of the environment like clutters.
It is more likely NBI has better sparsity structure than the
channel information. The precision of the recovery algorithm
relies on the signal sparsity. The NBI signal recovery should be
more accurate. Thus, the estimated interference signal j̃f can
be used to cancel the interference in the received waveform.
By subtracting the interference signal jf using the estimated
interference signal j̃f , the remaining signal w̄ is obtained.
ȳj0 =

I
D/A Conveter
Dual- core
1Gsps/14- bit

2
2

6 ε001 ,
6 ε002

(14)

bounds the energy of noise signal w̄j + n̄j .
IV. CS- BASED MIMO RADAR SYSTEM

Based on the analysis on the previous section, a CS-based
cognitive MIMO radar system is proposed. The transmitter
architecture is the same as our previous work [3], as shown
in Fig. 2. In the experiment, the carrier frequency is set to
4 GHz. The bandwidth of the radar waveform is 500 MHz,
making the system an ultra-wideband system.
To recover the original signal, the sampling rate of the
radar received waveform must be higher than 500 MHz in

With one target in surveillance area, the experiment result is
shown in Fig. 4. The starting time of recording received signal
is the starting time of waveform transmitting. The position of
radar waveform in the received signal can be detected by using
the correlation relationship of CP. The channel information
should have only one path in this situation. However, in real
system, the channel delay is not an integer in model (3). Fig. 4
(b) is the result by solving (14). Most energy is mapped to
tap 3 while the neighborhood of the tap holds the residue.
Since sparsity is considered, the values of all other taps are
close to zero. Fig. 4 (c) Interference information is obtained
by solving (11). All of the energy is mapped to tap 45.
In order to analyze the performance of the proposed algorithms, several simulations are conducted. For each point,
1 × 104 simulations are performed. Fig. 5 shows the effect of
single tone NBI to compressed sensing system when SNR =
20 dB. Rhombused curve is obtained by (9) while stared curve
is obtained by recovering channel information using (12). The
vertical axis is the probability of detection, which means the
probability when correct channel tap location is obtained. The
result shows when SIR goes below 15 dB, system performance
would be severely influenced by NBI.

4

20
100
0

50

0

1

0

0.9

-10

0.8

-20

0.7

-30

-40

0.6
Pd

Amplitude (dB)

Amplitude (dB)

Amplitude

-20

10

-40

0.5

-50

0.4

-60

With second step, single tone NBI
With only one step, single tone NBI
With second step, dual tone NBI
With only first step, dual tone NBI

0.3

-60
-50

-70

0.2

-80

0.1
-20

-80
-90

-100

0

4

2
4 x 10
Time index

-100

-10
SNR (dB)

-5

0

-100
0
20
40
Channel tap index

0

20
40
60
Frequency bin

Fig. 4. Experiment data: (a) Received sub-sampled baseband signal; (b)
Recovered channel information; (c) Recovered interference information
1
0.9

Fig. 6.

Performance curve of simulation

algorithm in real-world situation. Simulation results show that
the proposed algorithms are able to estimate the locations
of the jamming interference and provide a performance gain
compared to the system without NBI cancellation.

0.8

ACKNOWLEDGMENT

0.7

This work is funded by National Science Foundation
through two grants (ECCS-0901420 and ECCS-0821658), and
Office of Naval Research through two grants (N00010-10-10810 and N00014-11-1-0006). This work is partly funded by
an AFOSR subcontract though the prime contractor (RNET
Technologies, Inc.) on the CIRE (Center for Innovative Radar
Engineering) contract FA8650- 10-D-1750 with AFRL/RY.

Pd

0.6
0.5
0.4
0.3
Without NBI cancellation
With NBI cancellation

0.2
0.1
12

Fig. 5.

-15

12.5

13

13.5

14
SIR (dB)

14.5

15

15.5

16

The effect of NBI cancellation

With a strong interference (SIR = 0 dB), Fig. 6 gives system
performance under difference additive white Gaussian noise
(AWGN) levels. Three targets situation is considered. Both
single tone and dual tone NBIs results are shown. Two stared
curves are for (14). Two rhombused curves are for (12). With
the second step of channel estimation, the accuracy becomes
much higher. In the case of dual-tone, the NBI signal becomes
less sparse than the single tone case. Because the CS algorithm
highly relies on the sparsity of the signal, the performance is
worse than the single tone case.
VI. C ONCLUSION
In this letter, we have proposed a CS-based MIMO UWB
cognitive radar. An algorithm is designed to get the ToA information from the sub-sampled echoed radar waveform. The
algorithm is based on circulant convolution. It can estimate
the channel information. Based on the equations of the whole
system, a hardware architecture has been proposed to fit into
the special structure of the CS system. NBI signal in the
surveillance area is also considered. An algorithm to estimate
and cancel the NBI interference is proposed. The hardware
testbed is used to verify the effectiveness of the proposed

R EFERENCES
[1] E. Candès, J. Romberg, and T. Tao, “Robust uncertainty principles: Exact
signal reconstruction from highly incomplete frequency information,”
Information Theory, IEEE Transactions on, vol. 52, no. 2, pp. 489–509,
2006.
[2] P. Zhang, Z. Hu, R. Qiu, and B. Sadler, “A compressed sensing based
ultra-wideband communication system,” in Communications, 2009.
ICC’09. IEEE International Conference on. IEEE, 2009, pp. 1–5.
[3] Z. Chen, C. Zhang, F. Lin, J. Yu, X. Li, Y. Song, R. Ranganathan,
N. Guo, and R. C. Qiu, “Towards a large-scale cognitive radio network:
Testbed, intensive computing, frequency agility and security,” in Computing, Networking and Communications (ICNC), 2012 International
Conference on. IEEE, 2012, pp. 556–562.
[4] X. Li, Z. Hu, R. Qiu, N. Guo, G. Scalzi, and M. Wicks, “Waveform
diversity for closed-loop wide-band cognitive sensing system,” in IEEE
Waveform Diversity and Design Conference, vol. 2012, 2012.
[5] X. Li, Z. Hu, R. Qiu, Z. Wu, J. Browning, and M. Wicks, “Demonstration
of cognitive radar for target localization under interference,” Aerospace
and Electronic Systems, IEEE Transactions on, vol. 50, no. 4, pp. 2440–
2455, 2014.
[6] X. Li, Z. Hu, R. C. Qiu, and M. Wicks, “Experimental demonstration
of cognitive radar for target localization under strong interference,”
in MILITARY COMMUNICATIONS CONFERENCE, 2012-MILCOM
2012. IEEE, 2012, pp. 1–6.
[7] J. Van de Beek, O. Edfors, M. Sandell, S. Wilson, and P. Borjesson,
“On channel estimation in ofdm systems,” in Vehicular Technology
Conference, 1995 IEEE 45th, vol. 2. IEEE, 1995, pp. 815–819.
[8] H. Rauhut, J. Romberg, and J. A. Tropp, “Restricted isometries for partial random circulant matrices,” Applied and Computational Harmonic
Analysis, vol. 32, no. 2, pp. 242–254, 2012.
[9] S. S. Chen, D. L. Donoho, and M. A. Saunders, “Atomic decomposition
by basis pursuit,” SIAM review, vol. 43, no. 1, pp. 129–159, 2001.
[10] E. J. Candes, Y. C. Eldar, D. Needell, and P. Randall, “Compressed
sensing with coherent and redundant dictionaries,” Applied and Computational Harmonic Analysis, vol. 31, no. 1, pp. 59–73, 2011.

