- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
International Journal of Antennas and Propagation
Volume 2014 (2014), Article ID 830869, 15 pages
Frequency Diverse Array Radar Cramér-Rao Lower Bounds for Estimating Direction, Range, and Velocity
1School of Communication and Information Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
2Department of Electrical and Electronic Engineering, Imperial College London, London SW7 2AZUK, UK
Received 27 January 2014; Accepted 19 February 2014; Published 8 April 2014
Academic Editor: Frankie KitWing Chan
Copyright © 2014 Yongbing Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Different from phased-array radar, frequency diverse array (FDA) radar offers range-dependent beampattern and thus provides new application potentials. But there is a fundamental question: what estimation performance can achieve for an FDA radar? In this paper, we derive FDA radar Cramér-Rao lower bounds (CRLBs) for estimating direction, range (time delay), and velocity (Doppler shift). Two different data models including pre- and postmatched filtering are investigated separately. As the FDA radar has range-angle coupling, we use a simple transmit subaperturing strategy which divides the whole array into two subarrays, each uses a distinct frequency increment. Assuming temporally white Gaussian noise and linear frequency modulated transmit signal, extensive simulation examples are performed. When compared to conventional phased-array radar, FDA can yield better CRLBs for estimating the direction, range, and velocity. Moreover, the impacts of the element number and frequency increment are also analyzed. Simulation results show that the CRLBs decrease with the increase of the elements number and frequency increment.
The main task of a radar is to detect the existence of targets and estimate their unknown parameters, for example, range, velocity, and direction of arrival (DOA) [1–4]. All the elements in a phased-array antenna transmit the same signal at a fixed carrier frequency. By controlling the linear phase progression across elements, the beam can be steered to the desired direction precisely. But the beam steering is fixed in an angle for all the ranges and cannot mitigate undesirable range-dependent interferences. For this reason, the range and angle information of targets cannot be directly obtained for a conventional phased-array radar, which offers only the angle information. However, in some applications, it is desired that the beam pointing can be steerable to the range of interest.
Recently, frequency diverse array (FDA) is proposed in [5–7] as a method to achieve range-dependent beamforming, which can be used to suppress range-dependent interferences and clutter  or improve detection performance . The most important difference of FDA antenna from a conventional phased-array antenna is that a small amount of frequency increment compared to the carrier frequency is used across the array elements. The periodicity of FDA beampattern in range, angle, and time is studied in . In , the FDA is investigated from a simulation perspective and a low-cost FDA is designed. The FDA using chirp (linear frequency modulation) waveforms with different starting frequencies is analyzed in  to characterize the associated range-dependent beampattern. The application of FDA in synthetic aperture radar (SAR) is investigated in [13, 14]. Additional studies to exploit the range-dependent beampattern characteristics have been reported in .
FDA offers a range-angle-dependent beampattern, which is of great importance as this provides a potential for range-angle localization of targets. Therefore, this paper considers a fundamental question: what estimation performance can be achieved for an FDA radar? Referring to conventional phased-array range radar, the Cramér-Rao lower bound (CRLB) expressions for estimating velocity, direction, and Doppler shift are derived in [16, 17]. By assuming a narrowband signal model prior to matched filtering, the CRLBs for estimating the time delay, direction, and Doppler shift with a phased-array radar are detailed in [18, 19]. Furthermore, the CRLB for jointly estimating the attributes of a moving target using MIMO radar is derived in [20–22].
Different from the above literatures, this paper derives the CRLBs for jointly estimating the range, direction, and Doppler shift with an FDA radar. Two typical data models, namely, pre- and postmatched filtering are analyzed separately for static and moving targets, respectively. Since the range and angle of targets cannot be estimated directly due to the range-angle coupling, we use a simple transmit subaperturing strategy for the FDA radar. This method divides the whole array elements into two subarrays; each uses a different frequency increment. In doing so, the range and angle of targets can be solely estimated. The corresponding CRLBs for estimating direction, range, and velocity are derived. Furthermore, the impacts of the elements number and frequency increment on the CRLBs are also investigated.
The rest of the paper is organized as follows. The FDA is briefly introduced in Section 2. In Section 3, we consider the data model after matched filtering and derive the CRLB expressions for estimating the range and direction. In order to decouple the range-angle coupling response, we divide the whole array into two subarrays. The corresponding CRLBs are derived. In Section 4, the data model prior to matched filtering is investigated and the CRLB expressions for estimating the direction, range, and Doppler shift are derived, where temporally white Gaussian noise and chirp pulse signal are assumed. Finally, extensive simulation results are provided in Section 5. This paper is concluded in Section 6.
2. Frequency Diverse Array (FDA)
In conventional phased-array radars, it is assumed that an identical waveform is radiated from each element, excluding the amplitudes and phases. Different from phased-array radars, FDA elements can be either excited by the same waveform or different waveforms. For simplicity and without loss of generality, we consider a uniform linear array (ULA) FDA and assume that the waveforms radiated from each antenna element are identical but with a frequency increment of Hz. That is, the radiation frequency of the th element is where is the radar carrier frequency and is the number of the elements.
Supposing a narrowband transmit signal and taking the first element as the reference for the array, the steering vector is given by the following equation :where is the direction, is the range, is the element spacing, is the speed of light, and is the transpose. Since and , in an amplitude sense an approximation can be taken by ignoring the phase term because it has ignorable impacts. In this case, the steering vector can be simplified to
The FDA characteristics can be summarized as follows [7–12]. (1) If the frequency offset is fixed, the beam direction will vary as a function of the range ; that is, it is a range- dependent beam. (2) If the range is fixed, the beam direction will vary as a function of . This means that the FDA is frequency-increment-dependent (3) If the frequency increment across the array is not applied (i.e., ), the corresponding FDA is just a conventional ULA phased array.
3. CRLB with Signal Model after Matched Filtering
3.1. Data Models
Suppose that a target is located at . After matched filtering, the baseband equivalent of the complex valued signals at the receiver can be expressed as where is a constant for a given target, with , and is a zero-mean complex Gaussian white noise with spatial covariance: where is the expectation operator, is the conjugate transpose, is the spatial noise covariance matrix, is the noise power, and is an identity matrix.
3.2. CRLBs of Angle and Range Estimations
In the following, we discuss the CRLBs for several different cases separately.
3.2.1. Range Is Known and Angle Is Unknown
In this case, . The corresponding Fisher information matrix (FIM) is : where is the real part of the signal, is the th element of , and the is  Under the signal model (4), the FDA FIM with respect to is expressed as where is the signal-to-noise ratio (SNR). Accordingly, the CRLB is where denotes the inverse matrix. In particular, when , it is simplified to
Proof. See Appendix A.
It can be easily proved that That is to say, the FDA radar has better CRLB for angle estimation than phased-array radar.
3.2.2. Angle Is Known and Range Is Unknown
As phased-array radar has range-independent beam, here we only calculate the range estimation for FDA radar. The FIM with respect to can be expressed as The corresponding CRLB is
Proof. See Appendix B.
3.2.3. Both Angle and Range Are Unknown
In this case, the and CRLBs can be similarly determined as
Proof. See Appendix C.
When the target’s range and angle are estimated jointly, the range and angle CRLBs will be significantly degraded due to the range-angle coupling (which will be further investigated in Section 5). Consequently the range and angle of targets cannot be estimated directly by a standard FDA radar. To overcome this problem, we present a transmit subaperturing method for the FDA radar.
3.2.4. CRLBs of Transmit Subaperturing FDA Radar
To decouple the range and angle peaks and estimate both the range and angle of target, we divide the whole array into two equal subarrays [25, 26]. Suppose that the number of elements is ; each subarray has elements; namely, . The first subarray uses the frequency increment of , and the second subarray uses the frequency increment of , as shown in Figure 1. The resulting system is referred to as transmit subaperturing FDA (TS-FDA) radar.
Taking the first element of the first subarray as the reference, the new steering vector can be represented as the following equation:where
The corresponding FIM can be derived as where
The angle and range CRLBs are the first two diagonal elements of the inverse of the FIMwhere is the element at the th row and th column of the matrix.
4. CRLB with Signal Model Prior to Matched Filtering
4.1. Data Models
Unlike the first data model, the data model prior to matched filtering allows to estimate the Doppler shift (velocity). Suppose an -element antenna array receives a time delayed and Doppler-shifted echo of the transmitted signal , where is the center frequency. Knowing the time delay and Doppler shift (assuming a target with constant radial velocity), the range and radial component of velocity can be determined by and = . We denote the continuous-time signal as = , where is the sampling interval. Correspondingly, the time delay and Doppler shift in the sampled signal domain are and = , respectively .
After converting to baseband and sampling, the received signal at time becomes where is complex amplitude of the signal and is the additive noise .
We assume that the snapshots taken at cover the whole of a coherent processing interval (CPI). Therefore, the time duration of the CPI is = . Under the data model of (21), the complex amplitude is assumed to be an unknown deterministic constant during the CPI. To model the Doppler effect with a frequency shift, the radial component of the target velocity needs to be much smaller than the propagation speed (i.e., ). Then the time-bandwidth product of the complex envelope should be larger than . In addition, it is assumed that the propagation time of the signal across the array is much smaller than the reciprocal of the signal bandwidth, which is the narrowband array assumption in array processing.
Define the vector of unknown target parameters as . Stacking all samples into a single vector, (21) can be rewritten as where denotes the Kronecker product, . The noise is assumed to be zero-mean Gaussian, spatially and temporally correlated with spatiotemporal covariance  where is the temporal noise covariance matrix.
Under the above data model, the signal and noise parameters are disjoint and satisfy the space-time separability .
4.2. CRLBs of Angle, Range, and Doppler Shift
Suppose the received signal is completely covered by the observations and the sampling is dense (i.e., 0). The is assumed to be constant in the frequency band , where denotes the frequency in the continuous-time domain. The vector of Doppler shift and time delay in the continuous-time domain is defined as
For and with being an identity matrix. We define the signal power . The CRLBs expression for and follows from Appendix D as the following:where
The terms and are proportional to the root mean squared (RMS) bandwidth and RMS duration of , respectively. Note that the decoupling is a consequence of the assumed space-time separability of signal and noise models and the assumption of the complex amplitude as an unknown deterministic constant. In addition, since the signal and additive noise parameters are disjoint, the above CRLB expressions hold regardless of whether the spatiotemporal noise covariance is known or unknown.
4.3. CRLBs When Chirp Pulse Signal Is Employed
In this section, we derive the CRLB expressions for the range, Doppler shift, and direction estimation when is a chirp pulse signal with a large time-bandwidth product [31, 32]. It can be expressed as where is the pulse repetition interval, is the number of chirp pulses, and is expressed as  where is the chirp pulse duration, is the chirp bandwidth, and is the Heaviside step function.
Assume the time-bandwidth product of the pulse is . Using the signal given in (29) and (30), in continuous-time domain, we obtain the signal power , , , and . Thus, the CRLB expressions of and for FDA are derived as the following equation:
Specially, for the phased array (i.e., ), the CRLB is
When only one pulse () is used, the CRLB in (32) will be infinite because the model is not identifiable and the range and Doppler shift cannot be uniquely estimated , since phase-array radar has no range identity capability. In contrast, (31) does not have such a problem.
According to the previous discussion, the angle and range of targets cannot be estimated jointly due to the range-angle coupling. Using the similar transmit subaperturing approach presented in Section 3 and steering vector (15), the CRLBs of angle, range, and Doppler shift are derived as in the following equation:
5. Simulation and Verification
In this section, we consider several numerical examples that compare the CRLBs in different signal and noise models. Consider an X-band FDA radar with the carrier frequency GHz. We assume a ULA of elements used for transmitting. The array elements are spaced half of the wavelength apart from each other; namely, . One target of interest is supposed to reflect a plane wave that impinges on the array from direction of angle .
Under the signal model after matched filtering, Figures 2 and 3 compare the CRLBs according to (9) and (13), respectively. It can be noticed that the CRLBs are improved when a larger number of elements are employed. However, it has no significant difference when different frequency increments are used. This is because ; thus the frequency increment has a small impact on the CRLBs. In , a frequency offset selection strategy is derived to precisely steer the beam toward a fixed range with a desired angle.
Figure 4 shows the CRLBs of angle and range when both the angle and range are unknown. The CRLBs are significantly degraded due to the range-angle coupling. Consequently the range and angle of targets cannot be estimated directly by the FDA radar. However, the CRLBs decrease as the increase of the number of elements and frequency increment still holds. Moreover, generally more elements mean that better CRLBs performance can be achieved for the FDA radar.
To overcome the problem that the range and angle of targets cannot be estimated directly by the FDA radar, we use the transmit subaperturing strategy on the transmit frequency increments. Figure 5 shows the corresponding CRLBs where is assumed. It can be noticed that, to obtain a lower CRLB, the and should have inverse signs; that is, one is passive and the other is negative. One reason is that, in this case, the FDA radar has a wider system bandwidth. Figure 6 shows the CRLBs of angle and range when is employed. It can be noticed that the CRLBs performance improves with the increase of the sensor number.
Under the data model prior to matched filtering, we suppose the following signal parameters: bandwidth MHz, repetition period ms, and pulse duration μs. In this case, the approximate expressions given in (31) are valid because the transmitted has a large time-bandwidth product (). Figure 7 shows the CRLBs for direction, range, and Doppler shift as a function of SNR. Note that, when dB and are employed, we can get m and rad/s that corresponds to the CRLB for velocity is 3.27 cm/s (since and ). Since , the frequency increment has a small impact on the CRLBs. In addition, observe that the CRLB for and is block-diagonal (see (32)) and therefore decoupled; that is, remains the same whether or not is known, and similarly, is the same whether or not is known. The decoupling is a consequence of the assumed space-time separability of signal and noise models and the assumption of the complex amplitude as an unknown deterministic constant.
Figure 8 shows that the CRLBs versus SNR for different combinations of and . Comparing Figures 8 and 7, the CRLBs have been significantly improved. Likewise, comparing Figures 8 and 9, the CRLBs performance improves as the number of elements increases.
In this paper, we derive the CRLB to jointly estimate the attributes of a moving target using FDA radar and compute the corresponding CRLB expressions. First, we briefly introduce the FDA concept and make a summary on the FDA characteristics. Then, we consider two different data models, namely, pre- and postmatched filtering. Under different signal and noise models, we compute the CRLB expressions for estimating the range, direction, and Doppler shift. We demonstrate that the FDA radar beamforming is coupled in range and angle and that the target’s range and angle cannot be estimated directly by the FDA radar. To overcome this problem, this paper proposes a transmit subaperturing strategy for the FDA radar. In doing so, the range and angle of targets are estimated from the transmit-receive beamforming output. Moreover, we also specialize the CRLB results to the case of temporally white noise and a chirp pulse signal. Extensive simulation results verify the correctness of the derived CRLBs. It is shown that the CRLBs decrease with the increase of the number of elements and frequency increment. The CRLBs can be further improved through three aspects: increasing the number of elements, enhancing the system bandwidth by employing a larger frequency increment, and using transmit subaperturing strategy with more subarrays.
A. Derive the CRLB for Angle When Range Is Known
To derive the CRLB, we start with a well-known expression for the FIM under the data model in Section 3. We define the spatial noise covariance matrix as and signal-to-noise ratio (SNR) as . Suppose the target range is known, the FIM of is
For a phased-array radar, there is . We then have where, and
The FIM of the phased-array radar is
Similarly, for the FDA radar, there is . The derivation of with respect to is
Accordingly, the FIM of for the FDA can be expressed as
B. Derive the CRLB for Range When Angle Is Known
Under the data model in Section 3, when the direction is known, the parameter to be estimated is . The FIM of is
The derivation of with respect to for FDA is
The FIM of is thus given by
C. Derive the CRLB for Range and Angle
Under the data model in Section 3, both the angle and range are unknown. The range and angle of targets are estimated jointly. The FIM for parameters and can be expressed as We then have Since , then . We can get
D. General CRLB Results
Under the data model in Section 4, we derive the continuous CRLB expressions for temporally and spatially white noise and denote with being the power. We also start with FIM: where is the th element of and = . Consider For clarity, we rewrite Fisher’s information matrix as where
According to the matrix inversion lemma, the inverse matrix of is where the submatrix corresponds to the CRLB of which is nuisance and corresponds to the CRLB of and which are of interest. where
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The work described in this paper was supported in part by the National Natural Science Foundation of China under Grant 41101317, the Program for New Century Excellent Talents in University under Grant NCET-12-0095, Sichuan Province Science Fund for Distinguished Young Scholars under Grant 2013JQ0003, and Fundamental Research Fund for the Central Universities.
- F. Bandiera, M. Mancino, and G. Ricci, “Localization strategies for multiple point-like radar targets,” IEEE Transactions on Signal Processing, vol. 60, no. 12, pp. 6708–6712, 2012.
- D. R. Fuhrmann, J. P. Browning, and M. Rangaswamy, “Signaling strategies for the hybrid MIMO phased-array radar,” IEEE Journal on Selected Topics in Signal Processing, vol. 4, no. 1, pp. 66–78, 2010.
- S. Sen and A. Nehorai, “Adaptive OFDM radar for target detection in multipath scenarios,” IEEE Transactions on Signal Processing, vol. 59, no. 1, pp. 78–90, 2011.
- M. Compagnoni, P. Bestagini, F. Antonacci, A. Sarti, and S. Tubaro, “Localization of acoustic sources through the fitting of propagation cones using multiple independent arrays,” IEEE Transactions on Audio, Speech, and Language Processing, vol. 20, no. 7, pp. 1964–1975, 2012.
- P. Antonik, M. C. Wicks, H. D. Griffiths, and C. J. Baker, “Frequency diverse array radars,” in Proceedings of the IEEE Radar Conference (RADAR '06), pp. 215–217, Verona, NY, USA, April 2006.
- P. Antonik, M. C. Wicks, H. D. Griffiths, and C. J. Baker, “Multi-mission multi-mode waveform diversity,” in Proceedings of the IEEE Radar Conference (RADAR '06), pp. 580–582, Verona, NY, USA, April 2006.
- P. Antonik, H. D. Griffiths, and C. J. Baker, “Range dependent beamforming using element level waveform diversity,” in Proceedings of the International Waveform Diversity and Design Conference, pp. 1–4, Las Vegas, Nev, USA, January 2006.
- P. Baizert, T. B. Hale, M. A. Temple, and M. C. Wicks, “Forward-looking radar GMTI benefits using a linear frequency diverse array,” Electronics Letters, vol. 42, no. 22, pp. 1311–1312, 2006.
- B. W. Jung, R. S. Adve, and J. Chun, “Frequency diversity in multistatic radars,” in Proceedings of the IEEE Radar Conference (RADAR '08), pp. 1–6, Rome, Italy, May 2008.
- M. Secmen, S. Demir, A. Hizal, and T. Eker, “Frequency diverse array antenna with periodic time modulated pattern in range and angle,” in Proceedings of the IEEE Radar Conference (RADAR '07), pp. 427–430, Boston, Mass, USA, April 2007.
- J. Huang, K. F. Tong, and C. J. Baker, “Frequency diverse array with beam scanning feature,” in Proceedings of the IEEE Antennas and Propagation Conference (AP-S '08), pp. 1–4, San Diego, Calif, USA, July 2008.
- T. Higgins and S. D. Blunt, “Analysis of range-angle coupled beamforming with frequency-diverse chirps,” in Proceedings of the International Waveform Diversity and Design Conference (WDD '09), pp. 140–144, Orlando, Fla, USA, February 2009.
- J. Farooq, M. A. Temple, and M. A. Saville, “Application of frequency diverse arrays to synthetic aperture radar imaging,” in Proceedings of the International Conference on Electromagnetics in Advanced Applications (ICEAA '07), pp. 447–449, Torino, Italy, September 2007.
- J. Farooq, M. A. Temple, and M. A. Saville, “Exploiting frequency diverse array processing to improve SAR image resolution,” in Proceedings of the IEEE Radar Conference (RADAR '08), pp. 1–5, Rome, Italy, May 2008.
- W. Q. Wang, “Phased-MIMO radar with frequency diversity for rangedependent beamforming,” IEEE Sensors Journal, vol. 13, no. 8, pp. 1320–1328, 2013.
- A. L. Swindlehurst and P. Stoica, “Maximum likelihood methods in radar array signal processing,” Proceedings of the IEEE, vol. 86, no. 2, pp. 421–441, 1998.
- J. Ward, “Cramér-Rao bounds for target angle and Doppler estimation with space-time adaptive processing radar,” in Proceedings of the 29th Asilomar Conference on Signals, Systems and Computers, pp. 1198–1202, Pacific Grove, Calif, USA, November 1995.
- A. Dogandzic and A. Nehorai, “Cramér-Rao bounds for estimating range, velocity, and direction with an active array,” IEEE Transactions on Signal Processing, vol. 49, no. 6, pp. 1122–1137, 2001.
- A. Dogandzic and A. Nehorai, “Estimating range, velocity, and direction with a radar array,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '99), pp. 2773–2776, Phoenix, Ariz, USA, March 1999.
- J. J. Zhang, G. Maalouli, A. P. Suppappola, and D. Morrell, “Cramér-Rao lower bounds for the joint estimation of target attributes using MIMO radar,” in Proceedings of the International Waveform Diversity and Design Conference (WDD '09), pp. 103–107, Orlando, Fla, USA, February 2009.
- A. Hassanien, S. A. Vorobyov, and A. B. Gershman, “Moving target parameters estimation in noncoherent MIMO radar systems,” IEEE Transactions on Signal Processing, vol. 60, no. 5, pp. 2354–2361, 2012.
- Q. He, R. S. Blum, and A. M. Haimovich, “Noncoherent MIMO radar for location and velocity estimation: more antennas means better performance,” IEEE Transactions on Signal Processing, vol. 58, no. 7, pp. 3661–3680, 2010.
- S. M. Kay, Fundamentals of Statistical Signal Processing: Estimation Theory, vol. 1, Pearson, London, UK, 2nd edition, 2011.
- H. L. van Trees, Optimum Array Processing, John Wiley & Sons, New York, NY, USA, 2002.
- D. Wilcox and M. Sellathurai, “On MIMO radar subarrayed transmit beamforming,” IEEE Transactions on Signal Processing, vol. 60, no. 4, pp. 2076–2081, 2012.
- W. Q. Wang and H. Z. Shao, “Range-angle localization of targets by a double-pulse frequency diverse array radar,” IEEE Journal on Selected Topics in Signal Processing, vol. 8, no. 1, pp. 106–114, 2014.
- S. Gogineni and A. Nehorai, “Target estimation using sparse modeling for distributed MIMO radar,” IEEE Transactions on Signal Processing, vol. 59, no. 11, pp. 5315–5325, 2011.
- T. Li and A. Nehorai, “Maximum likelihood direction finding in spatially colored noise fields using sparse sensor arrays,” IEEE Transactions on Signal Processing, vol. 59, no. 3, pp. 1048–1062, 2011.
- S. Sen, “OFDM radar space-time adaptive processing by exploiting spatio-temporal sparsity,” IEEE Transactions on Signal Processing, vol. 61, no. 1, pp. 118–130, 2013.
- J. J. Blanz, A. Papathanassiou, M. Haardt, I. Furió, and P. W. Baier, “Smart antennas for combined DOA and joint channel estimation in time-slotted CDMA mobile radio systems with joint detection,” IEEE Transactions on Vehicular Technology, vol. 49, no. 2, pp. 293–306, 2000.
- P. N. Pathirana, S. C. K. Herath, and A. V. Savkin, “Multitarget tracking via space transformations using a single frequency continuous wave radar,” IEEE Transactions on Signal Processing, vol. 60, no. 10, pp. 5217–5229, 2012.
- P. H. Leong, T. D. Abhayapala, and T. A. Lamahewa, “Multiple target localization using wideband echo chirp signals,” IEEE Transactions on Signal Processing, vol. 61, no. 16, pp. 4077–4089, 2013.
- L. Zhuang, X. Liu, and W. Yu, “Precisely beam steering for frequency diverse arrays based on frequency offset selection,” in Proceedings of the International Radar Conference (RADAR '09), pp. 1–4, Bordeaux, France, December 2009.