ZHANG et al.: MULTITARGET AOA ESTIMATION USING WIDEBAND LFMCW SIGNAL AND TWO RECEIVER ANTENNAS 7103
where A
T
, f
min
, k, T
sweep
denote the amplitude of the signal,
the starting frequency of frequency modulation, the frequency
modulated slope, and the time duration of signal, respectively.
The instantaneous frequency of the LFMCW signal is
f =
1
2π
d(2πf
min
t + πkt
2
)
dt
= f
min
+ kt. (2)
In our system, the transmitted signal will be received by the
receiver after propagation. Suppose that the signal received by
the receiver antennas is reflected by a specific reflector. Then,
the received signal can be expressed as
r
0
(t)=A
R
e
j(2πf
min
(t−τ )+πk(t−τ )
2
)
, (3)
where A
R
and τ denote the reflected signal amplitude and the
time of flight (TOF), respectively. Comparing (1) and (3), we
can see that the received signal is a scaled and delayed version
of the transmitted signal.
According to [23], the amplitude A
R
can be expressed as
A
R
= A
T
R
√
G
T
G
R
d(t)
α
, (4)
where R denotes the reflection coefficient which is determined
by the radar cross section (RCS) of the reflector, G
T
and G
R
de-
note the gain of transmitter and receiver antennas, respectively,
d(t) denotes the propagation distance from the transmitter to the
reflector and back to the receiver, α is the attenuation coefficient
which is determined by the environment of propagation.
The time of flight τ is given by
τ =
d(t)
c
, (5)
where c is the velocity of light in the space.
As shown in Fig. 1, the received signal is the superposition
of a series of signals from different reflectors perturbed by the
noise of the receiver. Under the assumption that signals scatter-
ing and reflecting from unsmooth surface of the same reflector
are merged together because of the limited bandwidth, the ag-
gregated signal r(t) can be written as
r(t)=
N
i=1
A
Ri
e
j(2πf
min
(t−τ
i
)+πk(t−τ
i
)
2
)
+ w(t), (6)
where N is the number of reflectors, i denotes the index of
the ith path of propagation and w(t) is the zero-mean additive
complex white Gaussian noise at the receiver with variance σ
2
0
.
III. T
HE PROPOSED AOA ESTIMATION ALGORITHM
In this section, we describe the proposed AOA estimation
algorithm in detail. There are mainly three steps involved in the
proposed algorithm. First, the received signal is demodulated by
mixing with the transmitted signal to transform the wideband
non-stationary signal into superposition of a series of single-
tone signals. Second, the number of reflectors is determined
using hypothesis tests. Third, signals from different reflectors
are separated using bandpass filtering and the corresponding
AOA of each reflector is estimated using the phase difference
between two receiver antennas.
Fig. 2. An illustration of the transmitted signal and reflected signal.
A. Preprocessing of the Received Signal
LFMCW signal is a wideband non-stationary signal. Com-
pared to the transmitted signal, the received signal is attenuated
on amplitude and delayed in time. The frequency relation be-
tween the transmitted signal and received signal is shown in
Fig. 2. By utilizing the time delay characteristic of the received
signal, we can mix the received signal with the transmitted sig-
nal to extract the information about the reflector. Specifically,
the mixed signal can be written as
s
d
(t)=s(t)r
∗
(t)=
N
i=1
A
T
A
Ri
e
j(2πf
min
τ
i
+2πktτ
i
−πkτ
2
i
)
+ A
T
e
j(2πf
min
t+πkt
2
)
w
∗
(t). (7)
The first part in (7) corresponds to the signal that is the super-
position of a series of single-tone signals from different reflec-
tors. In t his paper, we assume that the distances of reflectors are
different. The frequency of this signal corresponds to the TOF,
which is determined by the distance from the transmitter to a
reflector and back to the receiver. In other words, the signals re-
flected from a specific distance have the same frequency and the
signals reflected from different distances have different frequen-
cies [24][25]. Note that although the Doppler shift will affect
the frequency of the received signal, when the frequency slope
of the LFMCW signal is high, the effect of the Doppler shift is
negligible. In such a case, bandpass filter can be used to sepa-
rate the signals from reflectors at different distances. The second
part in (7) corresponds to the cross-term of signal and noise. The
mixing process transforms the wideband non-stationary signals
into stationary signals with a single frequency, due to which
we can use classic signal processing techniques such as DFT to
extract the information of the signal.
In practice, the received signal is sampled by an A/D converter
to obtain its discrete version as follows:
s
d
(mT
s
)=
N
i=1
A
T
A
Ri
e
j(2πf
min
τ
i
+2πkmT
s
τ
i
−πkτ
2
i
)
+ A
T
e
j(2πf
min
mT
s
+πk(mT
s
)
2
)
w
∗
(mT
s
), (8)
where T
s
denotes the sampling interval.