SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Open Badges Research

Adaptive highly localized waveform design for multiple target tracking

Ioannis Kyriakides1*, Darryl Morrell2 and Antonia Papandreou-Suppappola1

Author Affiliations

1 School of Electrical, Computer and Energy Engineering, Arizona State University, Tempe, AZ, USA

2 Department of Engineering, Arizona State University East, Mesa, AZ, USA

For all author emails, please log on.

EURASIP Journal on Advances in Signal Processing 2012, 2012:180  doi:10.1186/1687-6180-2012-180

The electronic version of this article is the complete one and can be found online at: http://asp.eurasipjournals.com/content

Received:10 January 2012
Accepted:24 July 2012
Published:21 August 2012

© 2012 Kyriakides et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


When tracking multiple targets, radar measurements from weak targets are often masked by the ambiguity function (AF) sidelobes of the measurements from stronger targets. This results in deteriorated tracking performance and lost tracks. In this study, we consider the design of configurable waveforms whose AF sidelobes can be positioned to unmask weak targets. Specifically, we construct multicarrier phase-coded (MCPC) waveforms based on Björck constant amplitude zero-autocorrelation (CAZAC) sequences. The MCPC CAZAC waveforms exhibit wide regions in their AF surface without sidelobes and allow for selective positioning of sidelobes. We apply these waveforms in the context of a target tracker by selecting waveform parameters that minimize the expected tracking error. We show that this is accomplished by selecting the position of AF sidelobes to unmask weak targets. The target tracker is based on an independent partitions likelihood particle filter that is capable of processing the high-resolution measurements resulting from the Björck CAZAC sequences and tracks a fixed and known number of targets. Using simulations, we demonstrate the improvement in tracking performance when we adaptively select the MCPC CAZAC waveforms over tracking using non-adaptive waveform configurations or single-carrier phase-coded CAZAC waveforms.


When tracking multiple targets using radar sensors, weak targets are often difficult to observe in the presence of strong targets. This is because the ambiguity function (AF) sidelobes of measurements from strong targets are higher than the AF mainlobe of measurements originating from weak targets. As a result, the joint tracking performance of a multitarget tracker, expressed either in terms of mean-squared error (MSE) or percentage of lost tracks, is poor. The location and magnitude of the AF measurement sidelobes in the delay-Doppler plane are directly related to the location and magnitude of the AF sidelobes of the transmitted waveform. The AF is in turn defined by the type of transmitted signal and its parameters [1-3]. Therefore, there is a need to design configurable radar waveforms and develop an adaptive radar sensor configuration technique to position sidelobes from strong target returns away from the predicted locations of weak targets.

In [2,3], the processing of the return signal was performed by partitioning the delay-Doppler plane into resolution cells with fixed locations. These cells were constructed in such a way as to approximate a probability of detection contour that depends on the signal type and its parameters, which were assumed to be fixed. A detection in a resolution cell was declared based on the thresholded output of a matched filter placed on the centroid of the cell. However, the shape of the probability of detection contour is often not well approximated by a tessellating shape, resulting in measurement errors. Instead of using a fixed waveform, adaptive waveform techniques were used to minimize either the tracking error or validation gate volume in [4,5]. Moreover, in [6,7] waveform parameter adaptation was used to minimize track loss in the presence of clutter. In [8], the probability of track loss and a function of estimation error covariance was minimized by selecting both waveform parameters and detection thresholds for range and range rate tracking in clutter. In [9], the time to detect new targets was minimized by posing the problem as a partially observed Markov decision problem. In addition, in [10], the one step ahead expected information from the target kinematic model was maximized using the appropriate waveform selection.

The methods mentioned above rely on linear observation models that do not accurately represent physical systems with nonlinear characteristics. In [11,12], iterative adaptive waveform techniques were developed for nonlinear system models with a single target and using frequency-modulated waveforms. These adaptive waveform techniques assume that the measurements are processed using matched filters on a fixed grid.

In this study, we develop a tracker that selectively collects measurements based on the predicted target state (instead of a fixed grid for linear systems) to track a known and fixed number of targets. We implement the tracker using a new method that does not require the collection of measurements exhaustively on a fixed grid in the AF plane. The new independent partitions likelihood particle filter (IP-LPF) tracker adaptively configures multicarrier phase-coded (MCPC) waveforms [1] so that their AF sidelobes can be positioned in such a way as to not mask weak targets in the presence of strong targets. We also develop an adaptive configuration strategy to select the MCPC waveform parameters based on the relative positioning of the targets in the delay-Doppler plane. In contrast to previous methods [11,12], the proposed adaptive waveform selection method is not iterative; in contrast, waveform parameters are directly selected based on the state information on the weak target relative to the strong targets. The new IP-LPF tracker uses a proposal distribution that is based on the independent partitions algorithm [13,14] and the likelihood particle filter [15]. The particle filtering framework accommodates the propagation of resolution cells to locations of interest in the delay-Doppler plane based on prior information on the target state. It can also use the exact shape of the probability of detection contour as a resolution cell instead of the simplified tessellating shape. With this approach, we do not need to exhaustively collect measurements on all points of a fixed grid. Instead, the matched filter is matched to locations that are likely, given the belief about the target state; these locations are represented by each of the particles of the particle filter.

We employ MCPC waveforms in our work as their AFs can exhibit both wide regions with zero magnitude as well as non-zero sidelobes that can appropriately be positioned based on how the waveform parameter values are chosen. In order to construct these MCPC waveforms, we use multicarrier modulation and equal-length Björck constant amplitude zero-autocorrelation code (CAZAC) sequences [16-18] that are cyclic-shifted versions of one another. A Björck CAZAC sequence provides higher-resolution measurements than a linearly frequency-modulated chirp [2,3] due to its highly concentrated AF. The high concentration in the AF plane results in improved tracking performance as we demonstrated in [19] for the single target case. For the multitarget case, measurements from single-carrier phase-coded (SCPC) CAZAC sequences exhibit sidelobes that are spread in the delay-Doppler plane and could mask weak targets. Our proposed use of configurable MCPC CAZACs, on the other hand, can adaptively position the waveform sidelobes to not mask weak targets.

We configure the MCPC CAZAC parameters at each time step of the tracking scenario to minimize the predicted tracking error. The waveform parameters selected are the cyclic shift in frequency used to generate the waveform, the number of CAZAC sequences used, and the length of the CAZAC sequences. We present a computationally feasible method of selecting parameters to position sidelobes that accounts for the predicted relative positioning of the targets. Furthermore, our simulations demonstrate that the minimization of the predicted tracking error, achieved by selecting the MCPC waveform parameters, can be achieved by positioning the AF sidelobes such that a weak target is not masked in the presence of strong targets.

The rest of the article is organized as follows. In the following section, we present the MCPC CAZAC waveforms and investigate the properties of their AFs. In Section “IP-LPF algorithm”, we provide a detailed description of the IP-LPF algorithm and its application in minimizing the predicted tracking error. In Section “Adaptive waveform selection”, we integrate the IP-LPF with a waveform configuration algorithm and demonstrate its performance in tracking multiple targets in Section “Simulation results”.

MCPC CAZAC sequences

Björck CAZAC sequences

A CAZAC sequence ξ(m) with finite length M, has constant magnitude, |ξ(m)| = 1, m = 0,…,M − 1, and zero-autocorrelation, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M1">View MathML</a>, for n ≠ 0, where the addition is modulo M[17,20]. An example of a CAZAC sequence with quadratic phase is the Björck CAZAC sequence. For prime length M = 1, mod 4, it is given by [16,17]

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M2">View MathML</a>


where m, mod M (or m modulo M) is the remainder of the division m/M, and [ (m/M)] is the Legendre symbol that is given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M3">View MathML</a>

Björck CAZAC sequences are an attractive choice for target tracking with radar [17] as their constant amplitude allows for continuous transmission of peak power and can thus lead to increases in signal-to-noise ratio (SNR). They also exhibit very tight localization in the delay-Doppler plane that can enhance the range resolution and range-rate resolution of the measurements. The discrete AF of a Björck CAZAC sequence is given by [21]

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M4">View MathML</a>


where n and νare the discrete delay and Doppler parameters, respectively. Specifically, the AF exhibits a large spike at the origin (nν)=(0,0) of the discrete delay-Doppler plane, with very small sidelobes. An example of the AF of a Björck CAZAC of length M = 1,741 is shown in Figure 1.

thumbnailFigure 1. AF surface plot of a Björck CAZAC of length M = 1,741.

MCPC Björck CAZAC sequences

As the AF of a Björck CAZAC sequence is very highly localized, we want to exploit its properties to position AF sidelobes for multiple target tracking. In particular, we use the fact that a cyclic frequency-shifted CAZAC is also a CAZAC [20] and also that a sum of cyclic frequency-shifted sequences has an AF surface whose sidelobe locations depend on the difference in cyclic frequency shift, number of sequences, and sequence length. Note that, although cyclic permutations of CAZACs are possible in both time and frequency, we restrict our attention to frequency shifts as they result in wide zero regions in the AF plane and can better facilitate the adaptive positioning of the AF sidelobes.

The MCPC scheme combines multiple waveforms that are modulated by orthogonal carriers; the carriers are separated in frequency using orthogonal frequency division multiplexing [1]. The phase coding is required to reduce the bandwidth of the CAZAC sequence so that it meets transmission requirements. We use this scheme to form the MCPC CAZAC waveform by combining Q cyclically permuted Björck CAZAC sequences. Specifically, if we cyclic frequency shift ξ(m) in (1) using the frequency shift ζq(assumed mod M) to obtain the qth SCPC cyclic frequency-shifted CAZAC waveform

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M5">View MathML</a>


then the MCPC CAZAC waveform, modulated with carrier frequency ζc, is given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M6">View MathML</a>


where m = 0,…,MQ−1, and ⌊·⌋ denotes rounding down to the nearest integer. Note that we restrict ζq = qζ (mod M) in (3) as this selection of cyclic frequency shift causes the positioning of the sidelobes of the AF to depend on ζ, thus facilitating adaptive waveform configuration in our proposed algorithm. Thus, Θ = (QMζ) in (4) defines the three parameters of the MCPC CAZAC waveform.

When processing the SCPC CAZAC waveform in (3) using the AF in (2), the narrowband assumption is used that states that the transmitted waveform does not experience any time scale changes due to target motion. This assumption is valid since the time-bandwidth product of the waveform can be shown to be much less than <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M7">View MathML</a> as the speed of propagation in the air c is large, where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M8">View MathML</a> is the target range rate ( [22], p. 241). As we restrict the time-bandwidth product of an SCPC sequence and an MCPC sequence to be the same, the narrowband assumption also holds for MCPC CAZACs. Note also that we double the number of possible AFs by taking the Fourier transform (FT) of each of the MCPC CAZAC waveforms that we construct. The AF of the transformed waveform is equal to the AF of the original waveform with the delay and Doppler variables interchanged. This offers a convenient method of producing additional sidelobe positioning options with little effort.

AF surface of MCPC CAZAC waveforms

The AF surface of the unmodulated MCPC CAZAC waveform in (4) is given by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M9">View MathML</a>. Using (2), (3), and (4), the AF is given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M10">View MathML</a>


where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M11">View MathML</a> is the energy of sΘ(m) that is normalized to have the same energy as ξ(m). Next, we consider two separate cases of cyclic frequency shifts: ζ=0 and ζ > 0.

Zero cyclic frequency-shift

When ζ = 0, two of the exponential terms in (5) cancel out. We can also simplify the summations <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M12">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M13">View MathML</a> where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M14">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M15">View MathML</a> are integers (see [23] for the derivation details). The resulting AF of the MCPC CAZAC with Θ = (QM,0) becomes

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M16">View MathML</a>


We can see that the AF in (6) is non-zero only if <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M17">View MathML</a> is a multiple of Q, thus resulting in zero AF surface regions of width Q. Although these regions can be used to reveal weak targets at selected areas in the AF plane, we also need to reduce the sidelobes near the origin of the AF. The area in the delay-Doppler measurement plane near the AF origin is the area that is most commonly interrogated by the IP-LPF tracker when accurately tracking a target, as we will see in Section “IP-LPF algorithm”. Since we already have zero AF surface regions in the interval n = 1,…,Q − 1, we need to investigate the shape of the AF surface along the Doppler axis ν at n = 0. Setting <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M18">View MathML</a> in (6), we obtain the AF surface as

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M19">View MathML</a>

Since |ξ(m)| = 1 for all m, we conclude that the AF surface is non-zero only when νis an integer multiple of M. Therefore, the location of the sidelobes when n = 0 can also be chosen by adjusting the value of M. An example of this is shown in Figure 2 that depicts the AF surface of an MCPC CAZAC waveform with Θ = (130,13,0); all non-zero sidelobes exist when n is an integer multiple of Q = 130.

thumbnailFigure 2. AF surface of an MCPC Björck CAZAC with parameters: (a) Θ = (Q,M, ζ) = (130,13,0) and (b) Θ = (130,13,1).

Positive cyclic frequency-shift

When ζ > 0, we obtain higher diversity in the locations of the AF sidelobes. Specifically, when ζ ≠ 0 in (5), we observe that the terms ej2Π(⌊(m − n) / Q)⌋)qζ/Mand <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M20">View MathML</a> repeat multiple times in the summation in Equation (5). This is due to the summation of these terms over q = 1,…,Q − 1 and the modulo M / ζ effect of the two exponential functions which is due to M / ζ < Q. Therefore, we can factor the repeating terms and rewrite (5) as

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M21">View MathML</a>


Note that q and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M22">View MathML</a> now vary from 0 to ⌊Q/β⌋ − 1. We choose Q, M and ζ such that β = ⌊(M−1)/ζ⌋ + 1 is approximately a multiple of Qfor most choices of ζ = 1,…,M − 1. This eliminates the summation terms for q and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M23">View MathML</a> that fall between the values of (β(⌊Q/β⌋−1) + β−1) and Q − 1. These terms were omitted in (7), which explains the use of the approximation symbol, as they only cause a negligible variation of sidelobes in the AF surface compared to the exact expression. The accuracy of the above approximation can be verified using a numerical-based analysis, i.e., the generation of the AF surface using the Matlab code used in this work which is available to the reader upon request. The summation with respect to <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M24">View MathML</a> can be simplified using <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M25">View MathML</a>, where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M26">View MathML</a> is an integer. We then let <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M27">View MathML</a>, where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M28">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M29">View MathML</a>. Also, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M30">View MathML</a>, where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M31">View MathML</a> is an integer. We let <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M32">View MathML</a> with <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M33">View MathML</a> an integer and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M34">View MathML</a>. This simplifies Equation (7) to

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M35">View MathML</a>


This expression shows that non-zero values of the AF exist for <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M36">View MathML</a> integer and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M37">View MathML</a>, i.e., <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M38">View MathML</a> for integer <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M39">View MathML</a>. This provides for controlled size valleys in the AF surface.

We also examine what happens along the Doppler axis ν at zero delay and ζ > 0. If we set n = 0 or <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M40">View MathML</a> in (8), then the term <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M41">View MathML</a> reveals that the AF surface <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M42">View MathML</a> has sidelobes that periodically repeat with period βM. Evaluating the above expression at the in-between intervals, we can obtain the AF surface sidelobe values. We then choose to use only waveforms with parameters Q, M, and ζ with relatively low sidelobe levels in their AF surface with n = 0.

When ζ = 1, β = M, larger valleys appear in the AF surface. Specifically, the AF in (8) becomes:

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M43">View MathML</a>

Using <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M44">View MathML</a> since <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M45">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M46">View MathML</a> and, therefore, having <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M47">View MathML</a> above expression becomes

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M48">View MathML</a>

Then we note that the factor <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M49">View MathML</a> can only take the values of 0 if <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M50">View MathML</a> and −1 if <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M51">View MathML</a> since both m and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M52">View MathML</a> take values less than β. This implies that non-zero values of the AF surface <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M53">View MathML</a> only exist at delay locations such that <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M54">View MathML</a> or <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M55">View MathML</a> are multiples of β. For <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M56">View MathML</a> which restricts <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M57">View MathML</a> to be a multiple of βwith <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M58">View MathML</a>, and since <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M59">View MathML</a>, M<Q then indices of the waveform in the AF expression summation are very limited compared to the waveforms’ length MQ (i.e., m<M). Therefore, the case where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M60">View MathML</a> in the AF expression appears in a very small number of additive terms and is omitted in the following analysis. Letting <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M61">View MathML</a>, and using <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M62">View MathML</a>, we obtain (for details, see [23])

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M63">View MathML</a>

Since this implies that <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M64">View MathML</a>, thus <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M65">View MathML</a> with β = M, it follows that non-zero sidelobes of the AF surface appear at intervals of Q + (Q/M) in the delay. This is demonstrated in the AF surface of the MCPC CAZAC waveform with Θ = (130,13,1), as shown in Figure 2b.

In summary, the possibility of choosing the parameters Θ = (Q,M,ζ) of an MCPC CAZAC waveform, and also rotating the entire AF surface by choosing to take the FT of the waveform, enables us to position sidelobes in order to minimize the predicted MSE, as we will show in Section “Adaptive waveform selection”.

IP-LPF algorithm

Tracking model

We consider L targets moving in a two-dimensional (2D) plane, where the number of targets is fixed and known. The target dynamics are modeled by a linear, constant velocity model [24] given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M66">View MathML</a>


where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M67">View MathML</a> is the state vector for the lth target at time k, T denotes vector transpose, xl,k, yl,k and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M68">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M69">View MathML</a> are the position and velocity in Cartesian coordinates, respectively, the matrix F is given by F =[1 Δt 0 0;0 1 0 0;0 0 1 Δt;0 0 0 1] (with each row in square brackets), Δt is the time difference between observations, and vl,k is a zero-mean, additive white Gaussian process with diagonal covariance matrix <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M70">View MathML</a> that models target deviations from constant velocity. The model in (9) can be used to determine the kinematic prior probability distribution function <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M71">View MathML</a> for the lth target. The multitarget state vector is expressed in terms of the state vectors of each target as <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M72">View MathML</a>. Following [13], we refer to each component xl,k of Xk as a partition. Since we assume that the targets move independently, the multitarget kinematic prior distribution is given by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M73">View MathML</a>.

A radar sensor collects information on the range and range rate of the targets in the scene relative to the sensor by transmitting pulses and processing the returns after they are reflected by the targets. The return waveform provides range information, in the form of time delays, and range rate information, in the form of frequency shifts of the return waveforms relative to the transmitted waveform. Assuming point targets, the range and range rate for partition l at time step k, relative to the uth sensor, u = 1,…,U are given, respectively, by [11]<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M74">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M75">View MathML</a>, where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M76">View MathML</a> are the Cartesian coordinates of the location of the uth sensor, and the U sensors are assumed to transmit and receive waveforms independently. The discrete time shift value and discrete Doppler shift value at the u sensor, due to the lth target, are given by [11]<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M77">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M78">View MathML</a>, respectively, where round (·) transforms the real number to the nearest integer, c is the velocity of propagation in the medium, fc is the carrier frequency, Ts is the sampling period, and M is the total number of waveform samples.

Matched filter statistic

At every time step k, a signal s(m), m = 0,…,M − 1, is simultaneously transmitted from each sensor in different frequency bands to avoid interference. The received signal at the uth sensor after demodulation is a linear combination of the reflections from all L targets, and it is given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M79">View MathML</a>

Here, ζc = fcTs is the discrete carrier frequency of the transmitted waveform. The sum of random complex returns, Al,k, from many different target scatterers on target l are zero-mean, complex Gaussian with known variance <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M80">View MathML</a> and follow the Swerling I model [25]. Each target is assumed to have a different radar cross section (RCS) [26] and thus its return signal has a different strength that is represented by the variance of Al,k. It is also assumed that the return signal strength depends only on the target RCS and not on the distance between the sensor and the target; the distance is compensated for by amplifying returns that arrive later in time. The noise terms vu,k(m), u = 1,…,U, are assumed to be zero-mean complex Gaussian with variance 2 N0 and independent for each sensor. The SNR is given by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M81">View MathML</a>[2], where lw is the index of the weakest target and Es is the energy of the transmitted waveform from each sensor. When no target is present, du,k(m) = vu,k(m).

At the receiver, the return signal is matched filtered with a signal representing returns from Θ targets, at different time shifts <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M82">View MathML</a> and different frequency shifts <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M83">View MathML</a>, λ = 1,…,Θ. These time-frequency shifts are derived from the belief in target state using a particle filtering approach. The matched filter output is thus given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M84">View MathML</a>


where Md > M should be large enough to accommodate a maximum delay in the signal due to a reflection from the target. Moreover, Θ in (10) equals 1 when independently proposing one partition and Θ = L if particles of L partitions are proposed. The matched filter statistic that we will use for estimation is <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M85">View MathML</a>, and it is written in terms of the AF of s(m) in Equation (2).

Measurement likelihood

The statistical properties of the matched filter statistic yu,k depend on the fact that Al,k and vu,k(m) are independent, zero-mean, and complex Gaussian. It can be shown that <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M86">View MathML</a> in (10) is also complex Gaussian with zero mean and yu,k is exponentially distributed both under hypothesis H0 (not target is assumed present) and under hypothesis H1(L targets are assumed present). The two hypothesis formulation for the measurement likelihood is thus given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M87">View MathML</a>


where (see [23] for derivation)

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M88">View MathML</a>


When using the MCPC CAZAC waveforms to compute the measurement likelihoods in (11), we can reduce the computational complexity by approximating the above variance expressions. In particular, since the AF sidelobes of MCPC CAZAC waveforms are zero at the locations where, according to the belief in target state, the targets are expected to be we can set AFs(n,ν) to be 0 for n ≠ 0 and ν ≠ 0 in the above expressions. In addition, for the SCPC waveforms the above holds only approximately due to their non-zero, however, very low AF surface sidelobes. Also, using the fact that AFs(0,0) = 1, we can let <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M89">View MathML</a> for all l (where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M90">View MathML</a> is a nominal value that we choose, since we assume the target strength to be unknown). Based on this, we can then set <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M91">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M92">View MathML</a> in (12).

Likelihood partition sampling

The highly concentrated AF of a Björck CAZAC sequence provides a highly concentrated likelihood proposal distribution and a high measurement accuracy. However, the proposal process needs to be modified to sample particles from the likelihood instead of the kinematic prior since the former is much more localized than the latter. To achieve this, we use a likelihood particle filter [15], where the importance density depends on the measurements rather than the kinematic prior.

We propose to integrate the use of the likelihood proposal with the independent partition (IP) particle filtering [13,14] concept to efficiently propose particles. In the IP, we propose individual partitions of the multitarget state vector, each representing the state of a single target. We then combine the more accurate partition proposals into particles. The IP algorithm is an approximation to the joint multitarget probability density particle filter [14]; the approximation is accurate when the targets are well separated in the observation space. When targets are close in measurement space, their partitions cannot be independently proposed as described above. Due to our use of the Björck CAZAC sequences that have sharply peaked AFs, the measurements can be well approximated as independent [27]. Our resulting algorithm, the IP-LPF, belongs to the class of sequential partition algorithms [28]. Algorithms of this class propose partitions sequentially and then combine them into particles.

Specifically, we first independently evaluate likelihood values at discrete delay-Doppler bins for each partition. Using these values, we create histograms and sample partition states. Note that we narrow our bin selection to a region of probability of almost one in the kinematic prior partition sample. This is necessary to ensure a minimum number of bins to build the histogram and to ensure that the sample from the measurements is consistent with the kinematic prior. We then evaluate partition weights by combining measurements from the different sensors using the kinematic prior. Using the normalized partition weights, we independently sample values for each partition. We combine the sampled partitions into particles, compute the weights of the particles, and estimate and resample the particles.

Note that range information from two sensors, combined with kinematic prior knowledge, can provide adequate information to estimate the position of a single target [19]. Specifically, geometry is used to find the intersection between two circles whose radii (in Cartesian coordinates) are the ranges of the sensors. The two intersections of these circles provide two coordinate locations, one of which can be selected that agrees with the kinematic prior information. For multiple targets, there are multiple of these circles for each sensor and each target, and multiple intersection points that do not correspond to true coordinate locations. Therefore, the use of three sensors can help clear the ambiguity by providing fewer intersections of three circles. In order to avoid complicated geometry, we first process the returns of two sensors and sample Cartesian coordinate target locations using the likelihood. We then weight the sampled locations with measurements from the third sensor. Our method is for a general number of sensors equal or greater than three; in this work, we kept the number of sensors to a minimum of three.

The partitions sampling based on the likelihood is performed in two stages. In Stage 1, we utilize information from only two of the sensors in order to propose a preliminary set of partitions. This avoids the complex geometry required to sample Cartesian locations from range and range rate information obtained from three or more sensors. In Stage 2, we refine our partitions selection by sampling from the preliminary set of partitions created in the first stage using information from all the sensors.

Stage 1: partitions sampling

We start by propagating each state partition without noise. We let λdenote the proposed partition at time k and l denote the partition that represents the true state of the lth target. Assuming that we use a sequential importance resampling particle filter [15], the ith state particle, i=1,…,N, is given by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M93">View MathML</a>. Using <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M94">View MathML</a>, we obtain

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M95">View MathML</a>

We want to determine a region of delay-Doppler bins that could contain observations if the true state is <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M96">View MathML</a>. This region is obtained from the spread of the kinematic prior that determines the possible states of partition λ. If we assume for simplicity that the variances <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M97">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M98','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M98">View MathML</a> of the kinematic prior in the 2D (x,y) dimensions are equal, then with probability of almost one, the proposed particle will fall within 3σx from <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M99">View MathML</a> and within 3σy from <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M100">View MathML</a>. The maximum and minimum possible sampled x and y coordinates will then yield the maximum and minimum range. That is, if we assume that the target is at angle Π/2 with the sensor, then <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M101">View MathML</a>. The range would then increase/decrease by the amount <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M102">View MathML</a><a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M103">View MathML</a>. The delay would also assume minimum and maximum index values given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M104">View MathML</a>

Similarly, the minimum and maximum values can be obtained for the range rate and thus the Doppler

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M105">View MathML</a>

where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M106">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M107">View MathML</a>.

We form all combinations of indices for delay and Doppler that lie within the minimum and maximum delay and Doppler values using

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M108">View MathML</a>


<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M109','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M109">View MathML</a>


where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M110','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M110">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M111','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M111">View MathML</a>. Then, we evaluate the matched filter output at each of these values and for sensor u,

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M112','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M112">View MathML</a>


Note that we have used only one delay-Doppler pair <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M113','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M113">View MathML</a> in the template signal representing a single partition λ. The single partition likelihood for this delay-Doppler bin is given by:

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M114','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M114">View MathML</a>


We evaluate the likelihood ratio for each delay-Doppler bin as

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M115','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M115">View MathML</a>


We then obtain

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M116','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M116">View MathML</a>


and the normalized distribution

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M117','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M117">View MathML</a>


from which we sample <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M118','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M118">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M119','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M119">View MathML</a><a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M120','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M120">View MathML</a>, jn = <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M121','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M121">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M122','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M122">View MathML</a>, for each particle i and each sensor u = 1,2. The resulting sampled range and range rate and the bias are, respectively, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M123','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M123">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M124">View MathML</a>(2 fcTs), <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M125','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M125">View MathML</a>. The values of <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M126','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M126">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M127','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M127">View MathML</a>, in turn, yield proposed state values <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M128','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M128">View MathML</a>. This is accomplished by taking the intersection of two circles in the 2D Cartesian plane and choosing the intersection point that mostly agrees with the kinematic prior information. This process is illustrated in Figure 3. We note that the sampled partitions <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M129','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M129">View MathML</a> from Stage 1 are based on information provided from only two sensors. Therefore, some of these partitions may be incorrect, as previously explained. However, the value of <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M130','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M130">View MathML</a> can be used in Stage 2 to evaluate likelihoods for three sensors in order to sample proposal partitions and help remove partitions that have incorrectly been sampled.

thumbnailFigure 3. Schematic of the likelihood proposal process. Each particle <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M131','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M131">View MathML</a> is deterministically propagated forward (left top figure, arrows 1,2), the observation points for sensors 1 and 2 are defined (right top and right bottom), one point is sampled from each observation set of each sensor (arrows 3,4), and one of the two states <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M132','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M132">View MathML</a> formed by the sampled ranges in the Cartesian coordinates that agrees more with the prior is selected (left bottom).

Stage 2: partitions sampling

During Stage 1, we propose partitions <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M133','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M133">View MathML</a>, i = 1,…,N, from delay-Doppler bins associated with sensors u = 1,2. In Stage 2, we utilize the return signals transmitted by U ≥ 3 sensors to refine our choice of partitions and to more accurately represent the target state. We first describe the complexity in calculating the partition weights and the approximation we use to make the computation tractable before we provide the details on the sampling process.

After matched filtering, and using the locations in the delay-Doppler plane derived from the proposed partitions, the measurements from sensors u = 1,…,U are given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M134','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M134">View MathML</a>


where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M135','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M135">View MathML</a> is the delay-Doppler pair that corresponds to the state <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M136','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M136">View MathML</a> and the uth sensor, and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M137','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M137">View MathML</a> is the true target state xl,k. Therefore, the single partition likelihood function for each proposed partition λ of particle ɨ =1,…,N is given by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M138','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M138">View MathML</a>. Here, the hypothesis of particle ɨ and partition λ is that the partition state equals <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M139','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M139">View MathML</a> and not <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M140','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M140">View MathML</a> for i ≠ ɨ, while <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M141','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M141">View MathML</a>, i = 1,…,N, u = 1,…,U are the measurements obtained from matched filters at the delay-Doppler location defined by the particle proposed target state vectors <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M142','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M142">View MathML</a>, i = 1,…,N. However, each likelihood for sensor u is a multivariate exponential distribution [29] that grows in dimensionality as the number of particles N increases. We approximate the likelihood for each partition to be <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M143','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M143">View MathML</a>, where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M144','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M144">View MathML</a> denotes the likelihood that a target exists at <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M145','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M145">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M146','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M146">View MathML</a> denotes the likelihood that a target does not exist at <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M147','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M147">View MathML</a>. In [27], we show that the covariance between measurements <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M148','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M148">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M149','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M149">View MathML</a>, ɨ ≠ i, depends on the filter proximity (i.e., the closeness of <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M150','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M150">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M151','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M151">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M152','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M152">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M153','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M153">View MathML</a>) relative to the AF spread. Therefore, the measurement independence approximation for the Björck CAZAC is reasonable due to its concentrated AF. Using this approximation, the weights for partition λ of particle ɨ = 1,…,Nare

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M154','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M154">View MathML</a>


If we divide the right-hand side by the constant <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M155','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M155">View MathML</a>, and we use (17)–(19), we obtain

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M156','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M156">View MathML</a>

where the likelihood probability functions are given in (11) for a single target. This is then normalized

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M157','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M157">View MathML</a>



<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M158','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M158">View MathML</a>


We finally perform partition resampling, where we sample a partition index <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M159','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M159">View MathML</a>, ɨ =1,…,N, from the distribution of <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M160','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M160">View MathML</a> with replacement. The resulting selected partition has value <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M161','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M161">View MathML</a> and selection probability <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M162','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M162">View MathML</a>.

Particle weighting

After partition resampling, we assemble particles from the sampled partitions as <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M163','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M163">View MathML</a>. We weigh these particles with weights that incorporate prior and measurement information. To find the weight equation, we start by defining a measurement matrix Yk that is composed of measurements from <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M164','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M164">View MathML</a> and contains measurements from each of the U sensors. Specifically,

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M165','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M165">View MathML</a>


<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M166','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M166">View MathML</a>


The likelihood function (for a single partition case) for each proposed particle ɨ =1,…,N is given by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M167','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M167">View MathML</a>, i = 1,…,N, u = 1,…,U. Here, the hypothesis of particle ɨ is that the state equals <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M168','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M168">View MathML</a>, while <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M169','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M169">View MathML</a> are the measurements obtained from matched filters at the delay-Doppler location defined by the particle proposed target state vectors <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M170','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M170">View MathML</a>. Note that this likelihood is a multivariate exponential distribution [29]. Using similar arguments as for the single partition case, we approximate the likelihood for each particle to be <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M171','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M171">View MathML</a>, where <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M172','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M172">View MathML</a> is the likelihood given that the target state equals <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M173','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M173">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M174','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M174">View MathML</a> is the likelihood given that no targets exist having state <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M175','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M175">View MathML</a>.

The weight of particle ɨ [15] using the aforementioned assumptions is given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M176','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M176">View MathML</a>

Dividing by the constant <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M177','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M177">View MathML</a>, using the likelihood in (11), and normalizing the weights by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M178','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M178">View MathML</a>, we obtain the normalized weighs

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M179','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M179">View MathML</a>


The state estimate is thus given by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M180','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M180">View MathML</a>. The algorithm is outlined next.

IP-LPF algorithm

■ For each partition λ = 1,…,Λ and for each particle i = 1,…,N

Stage 1: Likelihood Partition Sampling

✶ Let <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M181','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M181">View MathML</a>

✶ For each sensor u = 1,…,U

✻ For <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M182','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M182">View MathML</a> and for <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M183','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M183">View MathML</a>

♦ Form <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M184','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M184">View MathML</a> using (13) and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M185','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M185">View MathML</a> using (14)

♦ Evaluate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M186','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M186">View MathML</a> using (15) and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M187','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M187">View MathML</a> using (19)

♦ Sample <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M188','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M188">View MathML</a>

♦ Let <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M189','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M189">View MathML</a>

♦ Let <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M190','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M190">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M191','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M191">View MathML</a>

✶ Calculate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M192','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M192">View MathML</a> from <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M193','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M193">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M194','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M194">View MathML</a>

Stage 2: Likelihood Partition Sampling

✶ For each particle ɨ =1,…,N

✻ Evaluate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M195','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M195">View MathML</a> using (20) and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M196','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M196">View MathML</a> using (22)

✻ Sample <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M197','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M197">View MathML</a>

✶ Let <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M198','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M198">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M199','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M199">View MathML</a>

Particle Weighting

■ For each particle i = 1,…,N

✶ Assemble particles <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M200','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M200">View MathML</a>

✶ Evaluate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M201','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M201">View MathML</a> using (24)

✶ For each particle ɨ =1,…,N

✻ Evaluate particle weights: <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M202','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M202">View MathML</a><a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M203','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M203">View MathML</a>

■ Estimate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M204','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M204">View MathML</a>

■ Increment k by 1

Adaptive waveform selection

In order to further improve tracking performance, we adaptively select the parameters of the MCPC CAZAC transmit waveform at each time step k so that we can minimize the predicted tracking root mean-squared error (RMSE). The three MCPC CAZAC parameters we consider are <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M205','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M205">View MathML</a> in (4), where Qk is the number of cyclically permuted Björck CAZAC sequences at time k, Mk is the length of the sequences at time k, and ζk is a parameter that controls the cyclic frequency shift at time k. The expected RMSE is given by the cost function

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M206','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M206">View MathML</a>


where the weighting matrix Cmakes the units of the cost function consistent by compensating for the differing units of the state vector. The subscript in the expectation operator E·[·] shows the dependance of the expected RMSE on the random target strength vector Ak, the random noise matrix vk, the unknown true target state Xk, and the estimate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M207','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M207">View MathML</a>, given the multitarget state estimate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M208','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M208">View MathML</a> at k−1 and the choice of Θk.

Next, we identify the set of values that the multitarget state estimate <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M209','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M209">View MathML</a> can take in terms of the delay-Doppler locations associated with it. As described in Section “IP-LPF algorithm”, in order to propose particles, we have considered a discrete finite set of delay-Doppler locations for each partition and each particle. This set corresponds to Cartesian coordinate locations that are most likely to occur according to the kinematic prior and the set of particles <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M210','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M210">View MathML</a>, i = 1,…,N generated at the previous time step k − 1. The set is given in (13) and (14) as <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M211','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M211">View MathML</a>, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M212','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M212">View MathML</a> and <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M213','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M213">View MathML</a>, for λ = 1,…,Λ, u = 1,…,2 and i = 1,…,N. We use index ȷ to denote a member of the set <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M214','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M214">View MathML</a>, of cardinality <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M215','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M215">View MathML</a>, consisting of the N particles that could be sampled by the IP-LPF proposal and subsequently weighted. Therefore, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M216','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M216">View MathML</a> is a large set including all combinations of possible delay-Doppler locations from two of the sensors for each target and particle. The process of forming partitions from sampled delay-Doppler locations is explained in Section “Stage 1: partitions sampling” and illustrated in Figure 3. Subsequently, one possible outcome of the likelihood sampling process and particle weighting is N weight-particle pairs <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M217">View MathML</a> corresponding to delay-Doppler locations <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M218','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M218">View MathML</a>, λ = 1,…,Λ, u = 1,…,U, i = 1,…,N. Similarly, based on the target motion model, we can identify a discrete finite set of possible true target states Xk. Each possible true target state <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M219','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M219">View MathML</a> with index <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M220','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M220">View MathML</a> is a member of the set <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M221','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M221">View MathML</a>, of cardinality <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M222','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M222">View MathML</a>. <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M223','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M223">View MathML</a> is related to corresponding delay-Doppler locations <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M224','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M224">View MathML</a>, l = 1,…,L, u = 1,…,U.

From the above, we may rewrite the cost function in (26) as

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M225','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M225">View MathML</a>

where the probability distributions <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M226','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M226">View MathML</a><a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M227','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M227">View MathML</a>, and p(vk) are defined in the context of the motion and measurement models in Sections “Tracking model” and “Matched filter statistic”.

In order to minimize the cost function, we need to minimize the probability <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M228','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M228">View MathML</a> and the particle weights <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M229','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M229">View MathML</a> in (25) for particles (i) such that <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M230','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M230">View MathML</a>. As the ȷth set of particles <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M231','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M231">View MathML</a> results from sampling by the IP-LPF, we will follow the sampling process of the IP-LPF and identify the selection probability for each partition of particles <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M232','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M232">View MathML</a>. According to Section “Stage 1: partitions sampling”, we obtain values <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M233','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M233">View MathML</a> for each partition λ = 1,…,Θ and each particle i = 1,…,N by sampling delay-Doppler bins from sensors u = 1,2 with probability <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M234','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M234">View MathML</a> given by (19). The values <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M235','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M235">View MathML</a> allow us to evaluate the likelihoods for the U sensors in order to sample partitions. In Section “Stage 2: partitions sampling”, we obtain partitions <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M236','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M236">View MathML</a> with selection probability <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M237','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M237">View MathML</a> given by (22). These sampled partitions are combined into particles into particles <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M238','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M238">View MathML</a> in Section “Particle weighting”. From the sampling process of each particle <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M239','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M239">View MathML</a>, we conclude that the probability of each particle being selected is <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M240','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M240">View MathML</a>. Therefore, any set of particles <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M241','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M241">View MathML</a> appears with probability <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M242','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M242">View MathML</a>. Furthermore, from (21), (22), for <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M243','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M243">View MathML</a> and from (17), (19) for <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M244','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M244">View MathML</a> we observe that the above sampling probabilities depend on the single partition likelihood ratio which using (16) is proportional to <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M245','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M245">View MathML</a>. Since <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M246','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M246">View MathML</a>, the selection probability monotonically increases with the matched filter statistic <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M247','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M247">View MathML</a>. Therefore, in order to minimize the above selection probability the matched filter statistic needs to be minimized for the delay-Doppler values in (13) and (14) with the additional constraint that <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M248','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M248">View MathML</a> for all partitions λ, particles i and sensors u. These sets of delay-Doppler locations correspond to the belief on target state as explained previously and only include delay-Doppler locations that imply erroneous target states <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M249','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M249">View MathML</a> (i.e., AF sidelobes). Since the matched filter statistic is a random variable it is minimized by minimizing its variance, given in (12) with Λ = 1, with respect to the waveform parameters.

Next, we observe that the particle weights <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M250','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M250">View MathML</a> in (25) contain the likelihood ratio both in the numerator and denominator. This, together with the fact that the prior has a wide spread compared to the likelihood, makes the particle weights nearly constant. Therefore, particle weights cannot be significantly reduced by adjusting the waveform parameters.

Therefore, the focus is on minimizing the matched filter statistic variance in (12) with respect to the waveform parameters specifically for the delay-Doppler values in (13) and (14) and such that <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M251','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M251">View MathML</a>. Since the matched filter statistic variance depends on the AF of the waveform, and since the above-mentioned set of delay-Doppler values refer to AF locations where the target is expected to be at time step k, excluding the AF mainlobe, the problem of minimizing the RMSE reduces to the problem of reducing AF sidelobes in the area where the target is expected to exist. This is a very well defined area in the delay-Doppler plane that is given by the sequential tracking process of the particle filter as explained in Section “IP-LPF algorithm”. Configuring the waveform so that zero sidelobes appear in selected areas of the AF surface was described in Section “AF surface of MCPC CAZAC waveforms”. Therefore, at each time step of the scenario, the parameters of the waveform to be transmitted at the next time step are selected such as to achieve low AF sidelobes in the areas where the weak target is expected to be found, resulting to a minimization of the expected RMSE. This is computationally efficient compared to iterative methods of waveform parameter selection [11,12]. However, the entire multitarget particle filtering method proposed is still associated with a large computational load which is not expected to reach real time operation with state-of-the art hardware which also limits the number of targets that can be tracked.

It is noted that this method works well only if the number of weak targets is low. The AF surface valleys created by these waveforms are, as expected, of limited size since the uncertainty cannot be entirely eliminated. Therefore, if multiple weak targets happen to be relatively positioned such that AF surface valleys cannot be configured to contain them then these targets will be masked. The problem of unmasking a larger number of weak targets is, therefore, an open problem and a limitation of the proposed method. Moreover, there is a prediction error associated with each target location which is estimated based on the Bayesian methodology employed. In this study, the prediction error is minimized as targets are highly localized when using the high resolution, high AF surface peaked CAZAC-based waveforms. The AF surface valleys designed are then large enough to contain this uncertainty and guarantee the unmasking of a weak target.

Furthermore, it is noted that the weak targets need to have a signal strength that is well above the noise level so that they are observable. In this study, it is assumed that what keeps weak targets masked are in fact the sidelobes from stronger measurement returns and not the noise. In order to initially detect weak targets, when no prior tracking information on their state is available, a sequential selective positioning of the sidelobes over different regions of the field of view is necessary. Once a weak target is detected and the tracking process begins then the selective positioning of the sidelobes based on prior tracking information described in this study is possible to take place.

Simulation results

We consider two simulation scenarios to demonstrate the performance of tracking multiple targets. The first scenario consists of one weak target and two strong targets; the second scenario consists of three targets of equal strength. Three different types of waveforms will be used: (a) SCPC Björck CAZAC, (b) MCPC Björck CAZAC with fixed parameters, and (c) MCPC Björck CAZAC with adaptively configured parameters. Three targets move in a 2D plane. The motion is completed in 199 time steps. Three sensors located at χ1 = −1000 m, ψ1 = 500 m, χ2 = 2500 m, ψ2 = 500 m, and χ3 = 500 m, ψ3 = 0 m collect range and range rate measurements. The trajectory of the target and the location of the sensors are shown in Figure 4. The target is assumed to move according to a nearly constant velocity model with covariance matrix Q = diag (225 64 225 64).

thumbnailFigure 4. Target trajectory and sensor location.

In the first scenario, the weak target l = 2 has a cross-sectional area such that, for SNR varying as 5, 10, 12, 15, 17, 20 dB, <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M252','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M252">View MathML</a>. The strong targets are characterized by <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M253','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M253">View MathML</a>. The noise variance is N0 = 1 and the waveform energy is Es = 1. In the second scenario, all three targets are observed with SNR that varies as 5,10,12,15,17,20 dB. The SCPC Björck CAZAC waveform has length M = 1,741. The choice of parameters of the MCPC waveforms was limited to combinations of values {M,Q} = {7,245},{11,154},{13,130}, and ζ = 0,1 in order to reduce computational expense in the adaptive selection process. The FT of the waveform was also used to introduce another degree of freedom by rotating the AF. All waveforms are sampled at 8 MHz and frequency modulated by fc = 40 GHz. The speed of propagation of the waveforms is c = 2.997925×108 m/s. For the simulations, we used N = 300 particles, initialized by drawing from a Gaussian distribution with mean the true initial target position and covariance Q0 = diag (1000 100 1000 100). The results were averaged over 300 Monte Carlo runs. The parameters for the adaptively configured MCPC waveform are selected at each time step as described in Section “Adaptive waveform selection”, while for the fixed MCPC the parameters were selected randomly at the beginning of the scenario.

For the first scenario, the RMSE tracking performance is shown in Figure 5a for different values of SNR and for all waveforms. The percentage of lost tracks is shown for each waveform and SNR value in Figure 5b. A lost track is declared if, for 6 consecutive steps, the tracking error exceeded 300 m. We observe that the MCPC waveform with adaptive configuration (indicated as AMCPC in all figures) clearly outperformed the SCPC Björck CAZAC and fixed MCPC waveform when considering the number of lost tracks. In terms of the RMSE, the SCPC Björck CAZAC appears to have similar performance as the adaptive MCPC case since both have the same measurement resolution. Performing well in RMSE is, however, not useful if it is accompanied with a high number of lost tracks. The non-adaptive MCPC waveform case has the lowest performance rating due to its high sidelobes that are not avoided during measurement. When there are no successful tracks, the RMSE value is shown as zero in Figure 5a,b.

thumbnailFigure 5. (a) RMSE versus SNR with 95% confidence intervals and (b) percentage of lost tracks versus SNR for three waveforms when tracking one weak and two strong targets.

The corresponding results from the second scenario are demonstrated in Figures 5a and 6a. We can observe that, if the targets have equal strength, then the adaptive MCPC and SCPC CAZAC perform similarly as their AF sidelobes do not mask weak targets. The large sidelobes of the non-adaptive MCPC, however, still result in large errors. Another observation is that in the first scenario, in the adaptive MCPC case, the results are improved compared to the second scenario. This is because in the first scenario, two of the targets have higher SNR values than in the second scenario.

thumbnailFigure 6. (a) RMSE versus SNR with 95% confidence intervals and (b) percentage of lost tracks for three waveforms when tracking targets of equal strengths.


We developed the IP-LPF algorithm, a particle filtering method based on the IPs approach and the likelihood particle filter, to track a fixed and known number of targets. A particle filter selects measurements based on the belief on the target state instead of collecting measurements exhaustively on a fixed grid. Moreover, the likelihood particle filter is capable of processing measurements resulting from the use of waveforms with high-resolution properties such as Björck CAZACs. In addition, we developed MCPC waveforms whose AF sidelobes can be adaptively positioned. We outlined a configuration strategy for selecting the parameter values of MCPC waveforms to position AF sidelobes such that weak targets are unmasked by minimizing the predicted MSE. We demonstrated with simulations that when tracking targets with different strengths using single Björck CAZAC or fixed parameter MCPC waveforms results in deteriorated tracking performance. On the other hand, the use of adaptively configured MCPC waveforms enables the tracking of weak targets in the presence of strong targets and offers significant tracking performance improvements.

Competing interests

The authors declare that they have no competing interests.


This study was supported under MURI Grant No. AFOSR FA9550-05-1-0443.


  1. N Levanon, E Mozeson, Radar Signals (Wiley, New York, 2004)

  2. C Rago, P Willett, Y Bar-Shalom, Detection-tracking performance with combined waveforms. IEEE Trans. Aerosp. Electron. Syst 34(2), 612–624 (1998). Publisher Full Text OpenURL

  3. R Niu, P Willett, Y Bar-Shalom, Tracking considerations in selection of radar waveform for range and range-rate measurements. IEEE Trans. Aerosp. Electron. Syst 38(2), 467–487 (2002). Publisher Full Text OpenURL

  4. DJ Kershaw, RJ Evans, Optimal waveform selection for tracking systems. IEEE Trans. Inf. Theory 40(5), 1536–1550 (1994). Publisher Full Text OpenURL

  5. DJ Kershaw, RJ Evans, Waveform selective probabilistic data association. IEEE Trans. Aerosp. Electron. Syst 33(4), 1180–1188 (1997)

  6. S-M Hong, RJ Evans, HS Shin, Control of waveforms and detection thresholds for optimal target tracking in clutter. IEEE Conference on Decision and Control, vol. 4 ((2000), pp), . 3906–3907 OpenURL

  7. SM Hong, RJ Evans, HS Shin, Optimization of waveform and detection threshold for target tracking in clutter. Proceedings of the 40th SICE Annual Conference ((2005), pp), . 42–47 OpenURL

  8. SM Hong, RJ Evans, HS Shin, Optimization of waveform and detection threshold for range and range-rate tracking in clutter. IEEE Trans. Aerosp. Electron. Syst 41(1), 17–33 (2005). Publisher Full Text OpenURL

  9. BF La Scala, W Moran, RJ Evans, Optimal adaptive waveform selection for target detection. International Conference on Radar ((2003), pp), . 492–496 OpenURL

  10. SD Howard, S Suvorova, W Moran, Waveform libraries for radar tracking applications. International Conference on Waveform Diversity and Design (2004)

  11. SP Sira, A Papandreou-Suppappola, D Morrell, Dynamic configuration of time-varying waveforms for agile sensing and tracking in clutter. IEEE Trans. Signal Process 7(55), 3207–3217 (2007)

  12. SP Sira, A Papandreou-Suppappola, D Morrell, Time-varying waveform selection and configuration for agile sensors in tracking applications. IEEE International Conference on Acoustics, Speech, and Signal Processing, vol. 5 ((2005), pp), . 881–884 OpenURL

  13. M Orton, W Fitzgerald, A Bayesian approach to tracking multiple targets using sensor arrays and particle filters. IEEE Trans. Signal Process 50(2), 216–223 (2002). Publisher Full Text OpenURL

  14. C Kreucher, K Kastella, AO Hero III, Tracking multiple targets using a particle filter representation of the joint multitarget probability density. SPIE Int. Symp. on Optical Science and Technology, vol. 5204 ((2003), pp), . 258–259 OpenURL

  15. MS Arulampalam, N Gordon, S Maskell, T Clapp, A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking. IEEE Trans. Signal Process 2(50), 174–188 (2002)

  16. G Björck, Functions of modulus one on <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M254','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/180/mathml/M254">View MathML</a> whose Fourier transforms have constant modulus, and cyclic n-roots. Proc. NATO Adv. Study Inst. on Recent Adv. in Fourier Analysis and its Applications, ed. by JS Byrnes, JL Byrnes ((1990), pp), . 131–140 OpenURL

  17. J Benedetto, I Konstantinidis, K Okoudjou, A Bourouihiya, Concatenating codes for improved ambiguity behavior. Int. Conf. on Electromagnetics in Advanced Applications ((2007), pp), . 464–467 OpenURL

  18. JJ Benedetto, I Konstantinidis, J Donatelli, C Shaw, A Doppler statistic for zero autocorrelation waveforms. Conf. on Inf. Sciences and Systems ((2006), pp), . 1403–1407 OpenURL

  19. I Kyriakides, D Morrell, JJ Benedetto, I Konstantinidis, A Papandreou-Suppappola, Target tracking using particle filtering and CAZAC sequences. Waveform Design and Diversity Conference ((2007), pp), . 367–371 OpenURL

  20. JJ Benedetto, JJ Donatelli, Ambiguity function and frame-theoretic properties of periodic zero-autocorrelation waveforms. IEEE J. Sel. Topics Signal Process 1(1), 6–20 (2007)

  21. A Kebo, JJ Benedetto, MR Dellomo, JM Sieracki, I Konstantinidis, Ambiguity and sidelobe behavior of CAZAC coded waveforms. IEEE Radar Conference ((2007), pp), . 99–103 OpenURL

  22. HLV Trees, in Detection Estimation and Modulation Theory, Part III (Wiley, New York, 1971) OpenURL

  23. Kyriakides I, On the use of Monte Carlo techniques for integrated sensing and processing (PhD thesis, Arizona State University, 2008), .

  24. GJ Foster, T Phan, JJ Petruzzo III, Track filtering of boosting targets. Proceedings of the 35th Southeastern Symposium on System Theory, vol. 35 ((2003), pp), . 450–454 OpenURL

  25. MI Skolnik, in Introduction to Radar Systems (McGraw-Hill, New York, 1980) OpenURL

  26. EF Knott, JF Shaeffer, MT Tuley, in Radar Cross Section, ((SciTech Publishing Inc), . , Raleigh NC, 2004) OpenURL

  27. Kyriakides I, On the validity of the measurement independence approximation when using single and MCPC waveforms based on Björck CAZAC sequences in multiple target radar tracking with a particle filter (Technical Report EEE-5-30-2008-1, Department of Electrical Engineering, Arizona State University, 2007, http://www), . ikyriakides.net webcite

  28. I Kyriakides, D Morrell, A Papandreou-Suppappola, Sequential Monte Carlo methods for tracking multiple targets with deterministic and stochastic constraints. IEEE Trans. Signal Process 56(3), 937–948 (2008)

  29. R Mallik, On multivariate Rayleigh and exponential distributions. IEEE Trans. Inf. Theory 49(6), 1499–1515 (2003). Publisher Full Text OpenURL