Target Tracking for Multistatic Radar with Transmitter Uncertainty∗ Sora Choi, Christian R. Berger, David Crouse, Peter Willett, and Shengli Zhou ECE Department, U-2157, University of Connecticut, Storrs CT 06269† ABSTRACT We present a target tracking system for a specific sort of passive radar, that using a Digital Audio/Video Broadcast (DAB/DVB) network for illuminators of opportunity. The system can measure bi-static range and range-rate. Angular information is assumed here unavailable. The DAB/DVB network operates in a single frequency mode; this means the same data stream is broadcast from multiple senders in the same frequency band. This supplies multiple measurements of each target using just one receiver, but introduces an additional ambiguity, as the signals from each sender are indistinguishable. This leads to a significant data association problem: as well as the usual target/measurement uncertainty there is additional “list” of illuminators that must be contended with. Our intention is to provide tracks directly in the geographic space, as opposed to a two-step procedure of formation of tracks in (bi-static) range and range-rate space to fuse these onto a map. We offer two solutions: one employing joint probabilistic data association (JPDA) based on an Extended Kalman Filter (EKF), and the other a particle filter. For the former, we explain a “super-target” approach to bring what might otherwise be a three-dimensional assignment list down to the two dimensions the JPDAF needs. The latter approach would seem prohibitive in computation even with these; as such, we discuss the use of a PMHT-like measurement model that greatly reduces the numerical load.
1. INTRODUCTION Passive radar is a bi-static system [5] that uses illuminators of opportunity to detect and track airborne targets. In a bi-static radar, sender and receiver are not co-located; in passive radar radio or television stations take the place of the sender and only the receiver is under control. Our interest is on a target tracking system using passive radar with digital broadcast signals: Digital Audio/Video Broadcast (DAB/DVB) [6,15]. It is important to recognize the difference between passive radar using DAB/DVB versus more traditional systems using, for example, commercial FM stations [12] or television [13]. In the traditional case the receiver observes a “direct blast” of signal followed by replicas representing reflections of the signal off whatever targets (and clutter) may be in the scene. Detection of these replicas is difficult because: • The observations process is continuous in time: there is no “pulse” or “scan” concept as one might find in a traditional surveillance system. • It is tempting to use the “direct” signal as a matched filter, to correlate for potential replicas. However, the signal contains both noise and reflections from the same targets at earlier times: the “known” signal in the matched filter is actually not known. • The reflections are themselves corrupted by the in-band direct-blast signal which is, essentially, noise. That is, the signal to noise ratio (SNR) is necessarily strongly negative. On the other hand, the DAB/DVB signals that we are interested in are orthogonal frequency-division multiplexed (OFDM). OFDM signals are relatively long packets of digital signals (actually, each amongst a large set of several hundred sinusoids modulates its own digital signature) and the target-reflection replicas are generally of a delay that falls well within the packet’s duration. That is, in contrast to the traditional case: ∗ This
work was supported by the US Office of Naval Research under contract N00014-09-10613.
[email protected].
† Contact:
• Each packet (or group of packets if information integration becomes necessary) gives rise to a “scan” of passive radar “hits.” • Since digital waveforms can be exactly demodulated essentially without error (television and radio would not work, otherwise) there is no concern about a dirty matched-filter template. • The reflections from targets can be thought of as forming the digital channel’s “impulse response” rather in the same way as multipath propagation would. From the perspective of OFDM these are a nuisance, but one that is easy to observe in the frequency domain as a predictable phase shift per carrier. That is, target reflections come for free as part of the OFDM demodulation process. Actually, a receiver in the OFDM-based DAB/DVB system would be able to “hear” several transmitters: the nearest, but also some others more distant. In some modulation schemes (as cellular telephony) this would be undesirable interference, and is avoided (in the cell case) by a re-use pattern of frequencies that avoids adjacent cells utilizing on the same bands. However, in OFDM as part of what one might think of as the equalization process the multiple transmissions become essentially overlaid, and, remarkably, actually enhance the SNR [18]. The European system is an almost perfect match to surveillance needs, and might ultimately be considered as deployable, in force-projection situations. The domestic (US) system uses 8VSB [19], and although there are advantages versus analog transmission, the benefits to surveillance are less clear. A passive radar can be covert, but the primary advantages in a non-hostile situation (where television and radio stations are available and working) are of cost and ubiquity: receivers are comparatively cheap, and with enough of them an air picture of considerable accuracy can be made available. But there are challenges, and here we shall focus on two of them. One is the measurements: The current receiver can give bi-static range and range-rate, but no angular measurements or only measurements of poor accuracy. Bi-static range by one illuminator indicates an ellipse on which target lies. Range-rate is dependent on location and velocity of target based on geometry of illuminator and receiver. Without angular information, to decide the location of target we need at least three illuminators. Coupled to this, another challenge is the unknown data association between measurements and illuminators. DAB/DVB operates in a so-called single frequency network (SFN) that introduces ambiguity: uncertainty of association between illuminator and signal. To address this challenge in passive radar, recent work has tried to extract additional features [11]. They proposed the use of the radar cross section (RCS). On the other hand, estimation of the trajectory in only range and range-rate has been suggested. One approach is target tracking based on the probability hypothesis density (PHD) [16]. Another approach uses a multi-stage Multi Hypothesis Tracking (MHT) as the solution for removing ghosts [8], and track initiation is incorporated in [3]. Estimation of the geo-trajectory is considerably more involved: one approach is to track in range and range-rate, and to fuse the resultant tracks. Our desire, however, is a direct geo-coordinate estimation algorithm, and we present two of these in this manuscript. The first is a JPDA/EKF. Considering the uncertainty of association between sender and signal, it becomes a 3-D association problem between measurements, targets, and illuminators. Three dimensional association implies high complexity and does not match the JPDAF model. To avoid it we will introduce the concept of a “super-target,” whereby 3-D association can be recast as 2-D association. The second is a particle filter. However, multi-target tracking in with particle filters is difficult due to data association: we need to calculate particles’ likelihoods, and in order to do this under the usual target tracking model we must consider all possible 3-D assignment tables from the target/measurement/illuminator lists. That is, of course, possible, but the numerical load is unattractive. Consequently here data association is simplified by choosing the measurement model used in probabilistic multiple hypothesis tracking (PMHT). To the PMHT, each measurement’s assignment is independent of all other measurements. In our problem, then, each delay can come from any illuminator/target (or false-alarm) pair; and the event that all measurements from the same pair is not impossible. In a sense, there is no concept of a “scan” and it is a perfectly valid PMHT implementation to update separately using each measurement, one at a time; but for us, the feature is that the particles’ likelihoods multiply across measurements.
Our manuscript is structured as follow. Section 2 establishes the model. The JPDA/EKF is described with simulation results in Section 3. In Section 4 we apply the PMHT-like particle filter to our problem with simulation results. We compare the two filters in Section 5, and summarize in Section 6.
2. SCENARIO (i) xs
There are several senders at = (xs , ys ) for i = 1, · · · , Ns and only one receiver located at xr = (xr , yr ). The receiver can measure range γ(t) and range rate γ(t). ˙ If the target is located at p(t) = (x(t), y(t)) with velocity v(t) = (x(t), ˙ y(t)), ˙ the state x(k) of the target is (x(t), x(t), ˙ y(t), y(t)). ˙ Then, (i) γ(x(t), x(i) s ) = ||p(t) − xr || + ||p(t) − xs ||
γ(x(t), ˙ x(i) s ) =
(p(t) − xr )T · v(t) + ||p(t) − xr ||
(i) (p(t) − xs )T · v(t) . (i) ||p(t) − xs ||
(1) (2)
We use the typical discrete dynamic system [2], described by x(k) = F(k)x(k − 1) + ν(k).
(3)
in which ν(k) is the sequence of zero-mean white Gaussian process noise with covariance Q(k). Although a linear plant is considered in this paper, it can be easily extended to nonlinear. The measurements at the receiver are delay and Doppler, so the observations model function h(x(k), xs ) is h(x(k), xs ) = [γ(x(k), xs ), γ(x(k), ˙ xs )]T .
(4)
Hence, the measurement equation is described by z(k) = h(x(k), x(i) s ) + ω(k)
(5)
with ω(k) being a sequence of zero-mean white Gaussian measurement noises with covariance R(k). As a reference for evaluation, the Cramer-Rao Lower bound (CRLB) of measurements based on range and range-rate for single time is considered, and is useful to identify illuminator/receiver/target geometries that are difficult. While the estimate at each time is using the previous estimated results, the measurement CRLB (MCRLB) does not reflect previous measurements, association, missed detections nor clutter and may not match closely the tracking performance. The likelihood function is µ ¶ Ns Y 1 1 (i) T −1 (i) Λ(x) = exp − (z − h(x, xs )) R (z − h(x, xs )) 2πR 2 i=1
(6)
The Fisher Information matrix, the inverse of the CRLB, is the following: F
= =
E[∇x ln Λ · ∇x ln Λ0 ] Ns X
(7)
−1 0 ∇x h(x, x(i) ∇x h(x, x(i) s )R s )
(8)
i=1
in which " ∇x h
=
∂γ(x,x(i) s ) ∂x (i) ∂ γ(x,x ˙ s ) ∂x
# (9)
and ∂γ ∂p ∂γ ∂v ∂ γ˙ ∂p ∂ γ˙ ∂v
(i)
=
p − xs ||p −
(i) xs ||
+
p − xr ||p − xr ||
(10)
= 0
(11) (i)
=
(i)
(˜ p − xs ) · (˜ v · (p − xs )) ||p −
(i) xs ||3
(i)
=
p − xs ||p −
(i) xs ||
+
+
(˜ p − xr ) · (˜ v · (p − xr )) ||p − xr ||3
p − xr ||p − xr ||
(12) (13)
˜ = (−y(t), ˜ = (y(t), −x(t)). where v ˙ x(t)) ˙ and p
3. JOINT PROBABILISTIC DATA ASSOCIATION VIA SUPER-TARGETS 3.1 The Algorithm The usual JPDA is based on measurement-to-target association probabilities (the “β’s) but does not consider illuminators as the association is assumed to be known. After introducing the general JPDA, we will extend the JPDA to incorporate an illuminator association list. The basic assumptions of the JPDA [1] are the following: 1. The number of target is known. 2. One measurement is associated with at most one target 3. Every measurement is independent on each other. Let us consider just measurement-to-target association. Association events show how all measurements are related to targets (and clutter, usually the 0th target). These events are winnowed via validation gates: If the event are feasible by gating, the event is chosen. Those events θ = [θjt ] can be expressed by event matrix: If a measurement j and a target t including clutter is associated, the entry θjt is 1, else 0. Based on feasible events, joint association event probabilities are calculated, and from these the marginal association probability for each measurement and each target can be derived. In our scenario, one more association list – illuminators – is required, and thus 3-dimensional association events need to be generated. However, if the probability of clutter does not vary with illuminator, 3-dimensional association is unnecessary. Let us consider a pair, of a target plus an illuminator, as a “super-target”. If we have three illuminators, then for each target there are three super-targets. By increasing the number of targets, 3-D association can be recast as 2-D. Note that the first of the two JPDA assumptions are satisfied under this formulation. The third is not: If three measurements come from one target via three different illuminators, those measurements cannot be independent. In a feasible joint association events θ, θj t˜ expresses association between a super-target t˜ and measurement j. The marginal association probability of each measurement and target is X βj t˜ , P {θj t˜|Z k } = P (θ|Z k ) (14) θ:θj t˜=1
Let M (k) be the number of measurements in the validation region at time k and Z k the set of measurements up to and including time k. By Bayes’ rule the joint association event probabilities P {θ(k)|Z k } are P {θ(k)|Z k }
=
1 P {Z(k)|θ(k), M (k), Z k−1 }P {θ(k)|M (k)} c
(15)
where c is a normalization constant. False measurements do not require illuminator association, and are as usual assumed to be uniform on the surveillance region of volume V , with prior pmf µF (φ) governing their number. Naturally, this must be selected with consideration of the number of illuminators. If τj is defined as an indicator to select the single measurement density, M (k)
P {Z(k)|θ(k), M (k), Z k−1 }
=
Y Y
j=1
p[zj (k)|θj t˜(k), Z k−1 ]
(16)
t˜ M (k)
= V
Y Y
−φ
j=1
{fj t˜[zj (k)]}τj
(17)
t˜
where fj t˜[zj (k)] = N [zj (k); zˆj t˜(k|k − 1), Sj t˜(k)] with the predicted measurement zˆt˜(k|k − 1) and innovation covariance St˜j (k). Let φ be the total number of false measurement in event θ, and δt˜ the indicator of supert˜ for super-target t˜, we have target detection corresponding to θ. With the probability of detection PD P {θ(k)|Z(k), M (k)}
=
Y ˜ φ! t˜ 1−δt˜ t δt˜ ) (1 − PD ) . µF (φ) (PD M (k)
(18)
t˜
as the event probability. Note that marginal association probabilities from a super-target t˜ relate in a pretabulated way as βj t˜ ↔ βjts for some target t and illuminator s. ˆ (k|k) is updated as: With the marginal association probability the estimate state x ˆ t (k|k) = x ˆ t (k|k − 1) + x
Ns X
Wst (k)νst (k)
(19)
s=1
where Wst (k) is the filter gain on sender s and target t and the combined innovation is m(k)
νst (k) =
X
βjts (k)(zj (k) − zˆst (k|k − 1)).
(20)
j=0
The covariance P (k|k) associated with the updated state is Pt (k|k)
= E{[x(k) − x ˆ(k|k)][x(k) − x ˆ(k|k)]0 |Z k } =
Ns X m X
βjts (k)E{[x(k) − x ˆ(k|k)][x(k) − x ˆ(k|k)]0 |θjts (k), Z k }
(21) (22)
s=1 j=0
=
Ns X m X
βjts (k)Pjts (k|k) +
s=1 j=0
Ns X m X
βjts (k)ˆ xjts x ˆ0jts − x ˆ(k|k)ˆ x(k|k)0
(23)
s=1 j=0
Ns X m X c = [( βjts (k))Pts (k|k) + β0ts (k)Pts (k|k − 1)] s=1
+
j=1 Ns X m X
βjts (k)ˆ xjts x ˆ0jts − x ˆt (k|k)ˆ xt (k|k)
(24)
s=1 j=0
The innovation covariance matrix is Sts (k) = R(k) + Hts (k)Pts (k|k − 1)Hts (k)0
(25)
Hts (k) = ∇x hs (k)|x=ˆxt (k|k−1) .
(26)
where
c The filter gain Wts (k) is obtained by Pts (k|k − 1)Hts (k)Sts (k)−1 . Here, Pts (k|k) is −1 Pts (k|k − 1) − Wts (k)Sts (k)Wts (k)0 = [I − Wts (k)H(k)]Pts (k|k − 1) ∗ [I − Wts (k)H(k)]0 + Wts (k)R(k)Wts (k)0 .
(27)
and the operation of the super-target JPDA is described. 35 Target 1 Target 2 Receiver Illuminator T=1
35 Target 1 Target 2 Receiver Illuminator T=1
30
30
25 25
20 20
15 15
10 10
5 5
0
5
10
15
20
25
0
5
10
30
(a) Simple case: no simultaneous crossing.
15
20 [Km]
25
30
35
40
(b) Challenging case: simultaneous crossing.
Figure 1. Example: Trajectories of targets.
3.2 Simulation Results Simulation results are presented for the super-target JPDAF using three senders and two moving targets. The number of false measurement µF (φ) follows the Poisson pmf: µF (φ) = e−λV
(λV )φ φ!
(28)
where λ is the clutter density, and the parametric JPDA is used. The targets are as in figure 1(a); note that they do not cross at the same time. The noise variance is σx2 = σy2 = 10−2 and σx2˙ = σy2˙ = 10−6 . Different probabilities of detection are explored, and 103 Monte Carlo trials are used. The RMSE of each target‡ is reported based on probabilities of detection 99% and 80%. Figures 2 and 3 show the RMSE of position and velocity. The maximum RMSE is 0.15 km for the target 1 and 0.16 km for the target 2 considering the probability of detection 99% and 80%. Especially with reference to the MCRLB, which illustrates the inadequacy of the measurements for target 1 near the end of its trajectory when it is essentially between an illuminator and the receiver, the JPDA appears to be working adequately. In figure 4 we consider a higher measurement noise level: variance σx2 = σy2 = 10−1 and σx2˙ = σy2˙ = 10−5 . The JPDAF appears to be quite sensitive to measurement noise level. ‡ We would not report RMSE for a track that is lost, which here we judge from a threshold on the number of gated measurements. It is interesting that with this value set to 10 – and a large value must be used due to the multiple illuminators – there were no declared lost tracks in this simulation.
0.2 PD 99% PD 80% MCRLB
PD 99% PD 80% MCRLB
0.18
0.3 0.16 0.25
0.14
RMS [Km]
RMS [Km]
0.12 0.2
0.15
0.1
0.08
0.06
0.1
0.04 0.05 0.02
0
10
20
30
40
50 Time [s]
60
70
80
90
0
100
10
20
30
(a) Target 1
40
50 Time [s]
60
70
80
90
100
70
80
90
100
(b) Target 2
Figure 2. RMSE of Position by JPDA, low measurement noise. −3
6
−3
x 10
8 PD 99% PD 80% MCRLB
x 10
PD 99% PD 80% MCRLB
7
5 6 4
RMS [Km/s]
RMS [Km/s]
5
3
4
3 2 2 1 1
0
0
10
20
30
40
50 Time [s]
60
70
80
90
0
100
0
10
(a) Target 1
20
30
40
50 Time [s]
60
(b) Target 2
Figure 3. RMSE of Velocity by JPDA, low measurement noise.
4. PARTICLE FILTERING VIA THE PMHT MEASUREMENT MODEL 4.1 The Algorithm A particle filter [17] approximates the posterior probability p(α(k)|Z k ) by particles α1 (k), · · · , αM (k) having probability mass w1 (k), · · · , wM (k). From Bayes’ rule, the posterior distribution is given: p(x(k)|Z k ) = R
p(Z(k)|x(k))p(x(k)|Z k−1 ) . p(Z(k)|x(k))p(x(k)|Z k−1 )dx(k)
(29)
From (29), p(α(k)|Z k ) ∝
p(Z(k)|α(k))p(α(k)|Z k−1 )
(30)
Np
=
p(Z(k)|α(k))
X j=1
p(α(k)|αj (k − 1))wj (k − 1)
(31)
−3
7 0.5
0.45
x 10
High R Low R MCRLB
High R Low R MCRLB
6
0.4
5 0.35
RMS [Km/s]
RMS [Km]
0.3
0.25
4
3
0.2
0.15
2
0.1
1 0.05
0
10
20
30
40
50 Time [s]
60
70
80
90
100
0
10
20
30
(a) Position
40
50 Time [s]
60
70
80
90
100
(b) Velocity
Figure 4. RMSE for target 1 by the JPDA based on different measurement noise levels.
PNp where Np is the number of particles. Considering j=1 p(α(k)|αj (k − 1))wj (k − 1) as a prior probability, the likelihood function p(Z(k)|α(k)) is inserted to get the posterior probability p(α(k)|Z k ). Thus we can sample from p(α(k)|Z k−1 ) by choosing particles with probability w(k) and drawing from p(α(k)|αj (k − 1)). We will use sampling/importance resampling (SIR)[9] to draw from p(α(k)|Z k ). In SIR, αj has associated weight wj p(Z(k)|αj ) wj = PNp i=1 p(Z(k)|αj )
(32)
To estimate trajectories of multiple targets, it is essential to approximate the posterior probability of multiple targets by particles. Particles can estimate all targets separately or together, which latter implies a stacked vector including all targets’ states. However, high dimension requires a prohibitive number of particles. To avoid this, the targets’ posterior probabilities are assumed independent. A key ingredient to particle filtering is the likelihood of the observations given a particle’s location. There is a significant complexity to this in our problem due to the 3-D assignment model described in Section 3. It is tempting to use the same “trick” of the super-target to simplify; this would be possible. But it must be recalled that even in the resultant 2-D problem measurements are strongly dependent, since each target can be represented by a number of measurements that may be as large as the number of illuminators, and a particle update must reflect all such information. As a result, we propose to simplify data association via the assumption 14] that all measurements are independent. The implication is that one target can give rise to more than one measurement. Given the multiplicity of illuminators this is sensible; however, under the PMHT measurement model there is very little structure, and indeed admits events in which more measurements can be ascribed to a target than there are illuminators. This is an inadequacy of the model, we recognize it and see if it matters. We have M (k)
p(Z(k)|αj , M (k)) =
Y
p(zi (k)|αj , M (k)),
(33)
p(zi (k)|αj , Ls )p(Ls |αj , M (k))
(34)
i=1
and p(zi (k)|αj , M (k)) =
Ns X s=0
where Ls means that the measurement is associated with illuminator s. If s = 0, zi (k) is from clutter and p(zi (k)|L0 ) = V +Ns P1d (Nt −1) approximately. For s = 1, · · · , Ns , p(zi (k)|αj , Ls ) = N (hs (αj ), R).
(35)
PNs PNs By the fact s=0 p(Ls |αj , M (k)) = 1, p(L0 |αj , M (k)) = 1 − s=1 p(Ls |αj , M (k)). From here, s > 0 is assumed and we use the concept of the super-target again. Since we consider the targets individually, each target/illuminator pair can be regarded as a different target. In [7], p(zi (k)|αj , M (k)) is calculated under same assumption. Then, 1 − Ns p(L1 |αj , M (k)) s=0 Pmin(M (k),Ns ) p(Ls |αj , M (k)) = (36) kξ(M (k)−k)(Nks )Pdk (1−Pd )Ns −k k=1 s 6= 0 Pmin(M (k),Ns ) Ns M (k) i=1 ξ(M (k)−i)(Nis )Pdi (1−Pd )Ns −i where ξ(k) =
(λV )k k!
exp(−λV ).
Based on the weights above, the SIR filter is implemented in the following steps: 1. Prediction: Draw new particles as αj (k|k − 1) = F(k)αj (k − 1) + v(k)
(37)
2. Weighting: Calculate the weights wj . 3. Sampling: Following wj , α1 (k), · · · , αNp (k) are resampled from α1 (k|k − 1), · · · , αNp (k|k − 1).
4.2 Simulation We will use the same setup as we explored in Section 3. With probability of detection 99%, the particle filters in figures 5 and 6 appear to work well, with maximum RMSE of position 0.075km. However in the 80% probability of detection case, the RMSE for target 1 increases§ , following a trend similar to that of the MCRLB: apparently the particle filter is more sensitive to the probability of detection than is the JPDAF. We consider a larger measurement noise, also as in Section 3. Regarding figure 7, the maximum RMSE of position is 0.075 km with the low noise and 0.1 km with the high noise. It is interesting that the particle filter appears to be sensitive to probability of detection but less so to measurement noise level, the opposite trend from the JPDAF.
5. SOME COMPARISON The particle filter and JPDAF are compared in figure 8 using the ground truth from figure 1(a). For target 1, it is hard to say which is better; however, for target 2, the particle filter is far preferable. We also investigate a more extreme case in which targets cross simultaneously, shown in figure 1(b). In figure 9 the RMSE of position is shown. Target 1, in particular, exposes the weakness of the particle filter approach. All told, it would appear that the particle filter is more accurate in less challenging situations, a result that is scarcely surprising considering that it is an approximation to the optimal nonlinear filter. However, in the case of close (crossing) tracks the particle filter’s independence assumption (and our decision not to adopt a “stacked” particle) appears to be a problem. §
There is a similar increase for target 2, but due to geometry it is less pronounced.
0.5
0.16 PD 99% PD 80% MCRLB
0.45
PD 99% PD 80% MCRLB
0.14
0.4 0.12 0.35 0.1 RMS [Km]
RMS [Km]
0.3
0.25
0.2
0.08
0.06
0.15 0.04 0.1 0.02
0.05
0
10
20
30
40
50 Time [s]
60
70
80
90
0
100
10
20
30
(a) Target 1
40
50 Time [s]
60
70
80
90
100
(b) Target 2
Figure 5. RMSE of position using the particle filter. 0.025
0.014 PD 99% PD 80% MCRLB
PD 99% PD 80% MCRLB 0.012
0.02 0.01
RMS [Km/s]
RMS [Km/s]
0.015 0.008
0.006
0.01
0.004 0.005 0.002
0
0
10
20
30
40
50 Time [s]
60
70
80
90
100
0
0
(a) Target 1
10
20
30
40
50 Time [s]
60
70
80
90
100
(b) Target 2
Figure 6. RMSE of velocity using the particle filter.
6. SUMMARY Passive radar has many advantages – covertness, propinquity, parsimony – and the version that we explore here, that based on OFDM illuminators from DAB/DVB signals, can also claim accuracy. But there are challenges. First, the measurements are of bi-static range and range-rate, and hence to decide the location of target we need at least three illuminators. And second, there is an additional layer of unknown data association between measurements and illuminators, since DAB/DVB ensures commercial radio and television end-user coverage by a network of illuminators each transmitting exactly the same waveform. Previous research has suggested that the best way to combat the latter is to form tracks in a proxy ellipse/ellipse-rate space, and thence to dispatch the former via fusion of these tracks. We would prefer a single-step solution, a tracker that operates directly in geographic coordinates. But then the challenges remain. We have suggested the two approaches: a modification on the JPDAF that proposes “super-targets” that comprise both illuminator and target; and a particle filter that follows the PMHT measurement model (that all measurements are conditionally independent). Both appear to work adequately. Our prime task for the near future is to test the approaches on real data, as in [4].
0.5
0.02
High R Low R MCRLB
0.4
0.016
0.35
0.014
0.3
0.012
0.25
0.01
0.2
0.008
0.15
0.006
0.1
0.004
0.05
0.002
0
10
High R Low R MCRLB
0.018
RMS [Km/s]
RMS [Km]
0.45
20
30
40
50 Time [s]
60
70
80
90
0
100
10
20
30
(a) Position
40
50 Time [s]
60
70
80
90
100
(b) Velocity
Figure 7. RMSE for target 1 using the particle filter, based on different measurement noise levels. 0.18 JPDA 1 PF 1 MCRLB
JPDA 2 PF 2 MCRLB
0.16
0.3 0.14 0.25
0.2
RMS [Km]
RMS [Km]
0.12
0.15
0.1
0.08
0.06 0.1 0.04 0.05 0.02
0
10
20
30
40
50 Time [s]
60
70
80
90
100
0
10
(a) Target 1
20
30
40
50 Time [s]
60
70
80
90
100
(b) Target 2
Figure 8. Comparison of approaches: the RMSR of position for the less challenging case of figure 1(a).
REFERENCES 1. Y. Bar-Shalom, X.-R. Li, Multitarget-Multisensor Tracking: Principles and Techniques, John Wiley & Sons, New York, NY, USA, 2001. 2. Y. Bar-Shalom, X.-R. Li, and T. Kirubarajan, Estimation with Application to Tracking and Navigation, YBS publishing, Storrs, CT, USA, 1995. 3. C. Berger, M. Daun, and W. Koch, “Low Complexity Track Initialization from a Small Set of Non-invertible Measurements,” EURASIP Journal on Advances in Signal Processing, Vol. 2008, DOI 756414. 4. C. Berger, B. Demissie, J. Heckenbach, S. Zhou, and P. Willett, “Signal Processing for Passive Radar using OFDM Waveforms,” submitted to the IEEE Journal on Selected Topics in Signal Processing (special issue on MIMO Radar). 5. M. Cherniakov (editor), Bistatic Radars: Emerging Technology, Wiley, 2008. 6. C. Coleman, H. Yardley, and R. Watson, “A Practical Bistatic Passive Radar System for Use with DAB and DRM Illuminators,” Proceedings of the IEEE Radar Conference, May 2008, pp. 15141519. 7. D. Crouse, M. Guerriero, and P. Willett, “A Critical Look at the PMHT,” Journal of Advances in Information Fusion, 2009.
1
0.18
0.8
0.16
0.7
0.14
0.6
0.12 RMS [Km]
RMS [Km]
0.9
0.2
JPDA 1 PF 1 MCRLB
0.5
0.1
0.4
0.08
0.3
0.06
0.2
0.04
0.1
0.02
0
10
20
30
40
50 Time [s]
60
(a) Target 1
70
80
90
100
JPDA 2 PF 2 MCRLB
0
10
20
30
40
50 Time [s]
60
70
80
90
100
(b) Target 2
Figure 9. Comparison of approaches: the RMSE of position for the more challenging case of figure 1(b).
8. M. Daun and W. Koch, “Multistatic Target Tracking for Non-Cooperative Illuminating by DAB/DVB-T,” Proceedings of OCEANS 2007. DIO 4302238. 9. N. Gordon, D. Salmond, and A. Smith, “Novel Approach to Nonlinear/Non-Gaussian Bayesian State Estimation,” IEE Proceedings-F, Vol. 140, No. 2 April 1993. 10. H. Grifths and C. Baker, “Passive Coherent Location Radar Systems. Part 1: Performance Prediction,” IEE Proceedings on Radar, Sonar & Navigation, Vol. 152, No. 3, pp. 153-159, June 2005. 11. S. Herman, and P. Moulin, “A Particle Filtering Approach to FM-Band Passive Radar Tracking and Automatic Target Recognition,” Proceedings of the IEEE Aerospace Conference, Vol. 4, pp. 1789-1808, 2002. 12. P. Howland, D. Maksimiuk and G. Reitsma, “FM Radio Based Bistatic Radar,” IEE Proceedings Radar, Sonar & Navigation, Vol. 152, No. 3, pp. 107-115, June 2005. 13. P. Howland, “Target Tracking Using Television-Based Bistatic Radar,” IEE Proceedings on Radar, Sonar & Navigation, Vol. 146, No. 3, pp. 166-174, June 1999. 14. C. Hue, J-P. Le Cadre, ”Tracking Multiple Objects with Particle Filtering,” IEEE Transactions on Aerospace and Electronic Systems, Vol. 38, No. 3 July 2002. 15. H. Kuschel, “VHF/UHF Radar, Part 1: Characteristics,” Electronics & Communication Engineering Journal, Vol. 14, No. 2, pp. 61-72, April 2002. 16. M. Tobias and A. Lanterman, ”Probability Hypothesis Density-based Multitarget Tracking with Bistatic Range and Doppler Observations,” IEE Proceedings on Radar, Sonar & Navigation, Vol. 152, No. 3, pp. 195-205, June 2005. 17. S. Arulampalam, S. Maskell, N. Gordon, and T. Clapp, “A Tutorial on Particle Filters for Online Nonlinear/Non-Gaussian Bayesian Tracking,” IEEE Transactions on Signal Processing, Vol. 50, No. 2 February 2002. 18. Z. Wang and G. Giannakis, “Wireless Multicarrier Communications: Where Fourier Meets Shannon,” IEEE Signal Processing Magazine, Vol. 17, No. 3, pp. 29-48, May 2000. 19. Y. Wu, S. Hirakawa, U. Reimers and J. Whitaker, “Overview of Digital Television Development Worldwide,” Proceedings of the IEEE, Volume 94, Issue 1, pp. 8-21, January 2006.