Research Article  Open Access
MingShih Huang, Ram M. Narayanan, Yan Zhang, Arthur Feinberg, "Tracking of Noncooperative Airborne Targets Using ADSB Signal and Radar Sensing", International Journal of Aerospace Engineering, vol. 2013, Article ID 521630, 12 pages, 2013. https://doi.org/10.1155/2013/521630
Tracking of Noncooperative Airborne Targets Using ADSB Signal and Radar Sensing
Abstract
As the Automatic Dependent SurveillanceBroadcast (ADSB) system has gained wide acceptance, additional exploitations of the radioed satellitebased information are topics of current interest. One such opportunity includes the augmentation of the communication ADSB signal with a random biphase modulation for concurrent use as a radar signal. This paper addresses the formulation and analysis of a suitable noncooperative multitarget tracking method for the ADSB radar system using radar ranging techniques and particle filter algorithms. In addition, the lowupdaterate measurement due to the ADSB system specification is discussed in order to provide acceptable estimation results. Simulation results show satisfactory tracking capability up to several kilometers with acceptable accuracy.
1. Introduction
During the past half century, air traffic management (ATM) has relied heavily on centralized ground stations using primary surveillance radar (PSR) and secondary surveillance radar (SSR) systems. With dramatic increase in the number of aircraft in the future, PSR and SSR will soon reach their limitations in their capability to monitor aerospace traffic patterns and aircraft locations in a timely manner to avoid collisions. Several decentralized approaches have drawn more attention lately thanks to their potential for quicker response and selfreliance. Independent of air traffic control (ATC), the widely used traffic alert and collision avoidance system (TCAS) monitors the airspace around an aircraft and communicates with other aircraft equipped with a corresponding transponder. TCAS is effective and successful in Class A, B, and C airspaces and Class E airspace above 10,000 feet, since aircraft in this airspace are required to be equipped with altitudereporting transponders. However, TCAS is ineffective for noncooperative aircraft and also in lowaltitude Class E and uncontrolled airspace. To enable the transformation of the ATM to a new paradigm that can meet the demand for the next 20 years and beyond, several developmental programs are underway, such as Single European Sky Air traffic Research System (SESAR) in Europe [1], Next Generation Air Transportation System (NextGen) in U.S.A. [2], and Collaborative Actions for Renovation of Air Traffic Systems (CARATS) in Japan [3]. The philosophy is to move away from legacy ground based technologies to a new and more dynamic satellite based technology. A key element of SESAR and NextGen is ADSB, which uses the Global Navigation Satellite System (GNSS) signals to provide air traffic controllers and pilots with precise position information in space, in contrast to the traditional surveillance radar derived data [4]. Aircraft transponders receive GNSS signals and use them to determine the aircraft’s precise location in the sky, which is combined with other relevant data and broadcast to other aircraft and air traffic control facilities via a digital data link. Besides ADSB’s wide acceptance in Europe and USA, NAV CANADA commenced operational application of ADSB as a means of providing aircraft surveillance information to Air Traffic Controllers (ATC) [5], and AirServices Australia commissioned the ADSB Upper Airspace Project (UAP), providing ADSB coverage across the whole continent [6]. In addition, ADSB is being used as the key solution for Unmanned Aerial System (UAS) integration in the National Airspace (NAS). When properly equipped with ADSB, both pilots and controllers will see the same realtime displays of air traffic, substantially improving safety and minimizing collision probability.
However, before the ADSB implementation and operation are fully complete, there will be a transition period involving coexistence of ADSB equipped and nonequipped aircraft. Similar to TCAS, the ADSB system works only for cooperative aircraft equipped with a corresponding system and is, thus, blind to noncooperative targets that are not similarly equipped. It has been pointed out that the limited use of ADSB as the sole means of surveillance may lead to a reduction of the integrity of the entire ATC system [7, 8]. Furthermore, localized problems, such as lack of connectivity to the required four visible satellites, will confuse not only aircraft pilots but also ATC [9, 10]. Hence, it is desired to find a way to cope with the noncooperative targets while retaining the benefits of the ADSB system.
Since the emergence of ADSB concept, some researchers have considered the utilization of existing and installed infrastructure of the surveillance radar to combine with the satellitebased ADSB system within the perspective of ATC. More interestingly, the use of the ADSB signal itself to detect noncooperative targets from the ADSB message and from the radar processing, as an onboard collision avoidance system, was first described in a patent disclosure [11], and the concept subsequently developed further [12–15]. The novelty of the ADSB radar system lies in that the system insightfully exploits the ADSB out signal, which is primarily designed for communication purposes, as a radar signal to perform multiple target estimation and tracking, thereby creating a multifunctional waveform. However, details of such radarcommunication system specifications and problems of lowupdaterate measurement due to the ADSB system requirement have not been fully discussed. This paper provides details of the ADSB radar system design, including the link power budget, and RF signal interference analysis. Furthermore, multitarget tracking using particle filter for the ADSB radar system is proposed. In addition, the slow convergence in particle filter algorithm and the measurement at lowupdaterate, owing to the ADSB signal broadcast rate of once per second, is addressed. It is important to note that in the ADSB radar system, the radar capability is independent of the ADSB message data integrity. Therefore, the coexistence of ADSB and radar sensing achieves increased reliability and redundancy, because it becomes possible to identify malfunctioning sensors and removes false alarms through performing conformance checks by comparing communication and radar reports.
The remainder of the paper is organized as follows. Section 2 briefly reviews the evolution of air surveillance systems. The ADSB radar system architecture and link budget are provided in Section 3. Then, the multitarget estimation algorithm for ADSB radar system is proposed and discussed in Section 4. Computer simulations that demonstrate the functionality of ADSB radar system for multitarget estimation performance are provided in Section 5. Conclusions are drawn in Section 6.
2. Evolution of Surveillance Technologies
For the sake of completeness and continuity, a brief and comparative review of the existing and new surveillance technologies is provided in this section. The major advantages and disadvantages for each approach are pointed out with an emphasis on the noncooperative surveillance capability. In addition, it can also be seen through the transition of the surveillance technologies that ATM is shifting from centralized systems to decentralized systems as the density of the air traffic continues to increase.
2.1. PSR and SSR
For the past few decades, PSR and SSR have been the main two components of an ATC station. PSR has the capability to detect large metal objects, including cooperative and noncooperative targets, while SSR works only for transponderequipped aircraft. SSR relies on aircraft with corresponding transponder, but the report provides aircraft identification. Both PSR and SSR were designed for low and medium traffic situations. Due to the insufficient range ambiguity, they are unable to meet the challenge of high traffic volume. Moreover, groundbased surveillance radar has limitations at low altitudes and may be influenced by atmospheric conditions.
2.2. TCAS
Due to continuing growth in air traffic, TCAS or other similar devices have been in various stages of research and development since the earlyto mid1950s to serve as a last resort collision avoidance safety net. TCAS operates similarly to the groundbased SSR but independently interrogates surrounding aircraft on a 1030MHz radio channel. The pilot will be alerted to the presence of the intruding aircraft replying to the interrogation via 1090MHz radio frequency. Current generation TCAS II, jointly developed by the US Radio Technical Commission for Aeronautics (DO185B) and European Organization for Civil Aviation Equipment (ED143), issues two types of advisories: the resolution advisory (RA), which identifies an intruder that is considered a collision threat, and the traffic advisory (TA), which identifies an intruder that may soon cause an RA. According to the predicted closest point of approach (CPA), TCAS produces a TA at approximately 45 seconds and an RA at 35 seconds to CPA. TCAS is designed to reduce the incidence of midair collisions and has been very successful since its introduction. However, the major concern with regard to either TCAS or ADSB is that they are not required all the time, and ADSB out is not yet mandated in most countries. Aircraft equipped with TCAS and/or ADSB are still exposed to danger of collisions in low altitude of Class E and uncontrolled airspace owing to their inability to detect noncooperative targets and unawareness of any illegal intruder in transponderrequired airspace.
2.3. ADSB
ADSB is redefining the paradigm of communication, navigation, and surveillance in ATM. An ADSB equipped aircraft determines its own position using GNSS and periodically broadcasts its fourdimensional position (latitude, longitude, altitude, and time), track and ground speed, aircraft or vehicle identification, and other additional relevant data as appropriate, for example, intended trajectories [16], to nearby aircraft also equipped with the ADSB system and potential ground stations without expectation of an acknowledgement or reply. One of the most significant advantages of the ADSB system is that it minimizes radio frequency (RF) spectral congestion as would be generated by TCAS. Any user, either aircraft or ground stations within broadcasting range, may receive and process ADSB surveillance information. ADSB system provides accurate information and improves situational awareness. Moreover, ADSB enables a shift from a centralized, groundbased ATM system to a decentralized network involving pilots and aeronautical operational control centers. It might eventually allow pilots to use onboard instruments and electronics to maintain a safe separation and to reduce their reliance on ground controllers.
2.4. ADSB Radar
While the ADSB system can only be considered as a decentralized cooperative surveillance ATM technique, the ADSB radar system deals with both cooperative and noncooperative targets. Figure 1 depicts the general working principle for ADSB radar concept, as well as the comparison between the ADSB system and the ADSB radar system. Noncooperative targets, such as unmanned aerial vehicles (UAVs) and private jets, are blind to the ADSB systems when not equipped with the ADSB system, but they pose equal collision danger if not being detected. Moreover, an inexpensive backup surveillance is necessary in the case of lost, incorrect, jammed, or compromised satellite information possibly caused by localized problem or device malfunction. The ADSB radar system is introduced as an insightful modification to the standard ADSB system in the interests of a safe backup service without significantly enlarging the volume of ADSB equipment. Since the ADSB systems broadcast signals periodically, the reflective ADSB electromagnetic energy could be viewed and extended for use as radar echoes. The returned echoes will be exploited to detect and track targets, including noncooperative ones, and hence, the pilot can be alerted without relying on guidance from the ground station to maintain aviation safety. Furthermore, by comparing the communication and radar reports for the ADSB equipped targets, the ADSB radar system can suppress the creation of tracks based on ADSB messages that contain intentionally incorrect position information.
(a)
(b)
3. System Design Specification
A conceptual architecture of the ADSB radar system is depicted in Figure 2. This differs from the corresponding figure in Zhang and Qiao [11] by pointing out the possibility to fuse the information from the radar processing and the ADSB message, as well as the target estimation method in radar processing. The fundamental principle of the ADSB radar system is to use reflective ADSB signals with only four additional Lband antennas on the standard ADSB system to perform multitarget estimation and tracking for both cooperative and noncooperative targets. Because the ADSB receivers only respond to the envelope of the incoming signal, phase modulation will not affect the normal ADSB signal decoding for other aircraft. As a result, the standard ADSB system will still be able to correctly decode the phasemodulated signal to obtain the communication message. The ADSB radar system retains the communication capability and receives the standard ADSB messages from neighboring aircraft, just as the standard ADSB system does. Once per second, it also broadcasts the phase modulated ADSB signal. After sending out the electromagnetic wave, the receiving antennas capture the returned signal reflected from surrounding metallic objects, and crosscorrelation is performed between the transmitted burst and the returned echoes. Both the ADSB In information and the estimated location from the radar report can be fed into Cockpit Display of Traffic Information (CDTI) [17].
The design consideration is focused on the feasibility of utilizing the returned ADSB signal as a radar echo in terms of the signaltonoise ratio and radar performance. In addition, the 500 W peak transmit power, the 120 μs pulse duration specifications of the 1090 MHz ADSB standard, and the broadcasting rate impose constraints on the system design, requiring a detailed analysis.
3.1. Signal Waveform and Interference Analysis
Each of the ADSB Extended Squitter is 120 μs long, including 8 μs preamble and 112 μs data block, and is similar to TCAS message format, except for longer 56 μs ADS information. Besides, the ADSB message is digitally encoded using pulse position modulation (PPM). When a communication signal is utilized and treated as a radar signal, many issues need to be addressed. First, the signal must be suitably modulated for performing matched filtering in the receiver. The product of the PPM signal waveforms for two different bits can only be either zero or positive, which translates into a raised correlation noise floor. Second, the constraints of pulse duration, peak transmit power, signal bandwidth, and coherent returned signals degrade the radar performance. Therefore, the primary task for the ADSB radar design is to exam the feasibility of treating the ADSB signal as a radar signal.
In order to improve the radar performance, the major modification to the ADSB signal is random biphase modulation introduced to each bit. Figure 3 shows the simulated ADSB waveform and the proposed ADSB radar waveform. Both signals are of the same duration and have identical digital messages, except that the ADSB radar signal has a random 180° phase shift. By randomizing the transmit signal, the matched filtering operation can be performed by crosscorrelating the reflected signal with a timedelayed replica of the transmit waveform. The biphase modulation renders both positive and negative products, forcing the autocorrelation to be statistically zero except for zero timelag. The effect of the sum of the products in the autocorrelation function is provided in Table 1. The significance of the random biphase modulation (0° or 180°) can be seen clearly in the autocorrelation of the ADSB radar signal and the ADSB signal, as depicted in Figure 4. It is clear that the autocorrelation of the ADSB radar signal outperforms that of the original modulationfree ADSB signal and shows an improved correlation peaktonoise ratio.

Since it is very likely to have other ADSB signals in the same channel, it is critical to analyze the interference of the ADSB signals and possibly other ADSB radar signals as well. Figure 5 shows that a typical ADSB radar signal is uncorrelated with both a standard ADSB signal as well as another independent ADSB radar signal, thereby indicating that the onboard ADSB radar receiver will not be affected by standard ADSB or ADSB radar transmissions from other aircraft in the vicinity.
Because the modulated signal has wider bandwidth (the resultant waveform still complies with the ADSB standard), the range resolution is improved and the range side lobes are also further suppressed [18]. The range resolution of 75 meters for the ADSB radar system is sufficient for most collision avoidance applications.
3.2. System Configuration and Link Budget Analysis
The ADSB radar system is composed of the following components: (a) ADSB transceiver, (b) ADSB encoder/decoder, (c) radio frequency (RF) electronics with up and downfrequency conversion, crosstalk cancellation, and filtering capabilities, (d) a phase modulator, and (e) four additional antennas operating at 1090 MHz. In particular, (d) and (e) are additional to the original ADSB system to process the reflected radar signals and to estimate target locations.
The link power budget is analyzed based on the radar range equation as follows: where and are the received and transmitted power, respectively; and are the receiving and transmitting antennas gain, respectively; is the signal wavelength; is the target’s radar crosssection (RCS); is the range to the target. To improve range resolution, common radar signal pulse durations are as short as tens of nanoseconds. However, since the signal energy is emitted continuously, integration gain can be achieved by utilizing the 545 sinusoidal cycles in each bit. As shown in Table 2, for a target having an RCS of 0 dBsm (1 square meter) at a distance of 4.5 km, the signaltonoise ratio (SNR) is 5.0 dB, a desirable value. For a large airliner with a higher RCS of 20 dBsm (100 square meters), the operational range for the same SNR value could extend to 14 km, which allows more than one minute of reaction time.

4. Multitarget Estimation and Tracking Algorithm
The Multiple Signal Classification (MUSIC) technique [19] was initially employed for the ADSB radar to perform detection and estimation of single target angleofarrival (AOA). Due to the nature of the ADSB radar signal, there are several important factors that limit the extension to multitarget estimation, primarily (a) the number of detectable targets and (b) signal coherency. The MUSIC algorithm has the wellknown inherent constraint that the number of the receivers has to be larger than that of the targets. With fixed number of antennas, the maximum number of detectable targets using the MUSIC algorithm is . While it is desirable to minimize the cost and amount of antennas used in the ADSB radar system, the detection capability becomes strictly limited. To further elaborate the signal coherency problem, the returned signals bouncing off from different targets are highly correlated because all the reflected radar echoes are simply scaled replicas of the transmitted signal but slightly frequency shifted by the Doppler shift, which is about 1 kHz for the relative speed of 150 m/s. Although there has been research dealing with correlated signals, such as constrained MUSIC [20], cumulantbased coherent signal subspace method [21], and focusing matrix for coherent signal subspace processing (for wideband signals) [22], these approaches do not apply to the scenario that the ADSB radar system encounters owing to either unavailable prior information or limited signal bandwidth from the ADSB signal specification. Other methods, like the least squares and the unscented Kalman filter approaches, require prior information of the motion model. A hybrid approach that combines the timeofarrival (TOA) and the particle filter approaches for multitarget estimation and tracking for the ADSB radar system is proposed in this paper. In the particle filter algorithm, the unavailable target motion can be dealt with using the concept of jittering [23], by introducing a noise model to cover the uncertainty of the target location. After the radar report is obtained, multisensor data fusion is a discrete problem, and attempts to combine ADSB messages with radar reports are described in [7, 8].
4.1. Multitarget Detection and Ranging through TOA
Let represent the transmitted signal. The distance between the target and the own aircraft can be obtained by crosscorrelating the delayed transmitted signal, , with the received signal, , where is the internal delay and is the roundtrip time to the target. The target range is , and the speed of light is . In theory, the crosscorrelation peak occurs when , from which the range can be determined. In our simulation experiment, the peak location may vary a little due to the measurement noise and phase coding scheme. In the ADSB radar system, is the phasemodulated ADSB signal. By continuously crosscorrelating the received signal with the delayed transmitted signal replicas, the number of peaks above a predetermined threshold indicates the number of detected targets. Moreover, the ranges can be calculated based on the delay times as described above. To correctly associate the ranges with the correct targets can be a difficult problem [24, 25]. However, because the sensors are located only several meters apart, the returned echoes from different targets usually arrive at the sensors in the same sequence.
4.2. Multitarget Estimation
The particle filter algorithm is briefly reviewed in this section. In addition, a simple multitarget estimation method based on particle filter is proposed for the application of the ADSB radar system. Furthermore, the problem of slow convergence in particle filter and measurements at lowupdaterate is addressed.
4.2.1. Particle Filter Algorithm
Particle filter [26–28], also known as sequential Monte Carlo method, bootstrap filtering [29], and the condensation algorithm in computer vision [30], is very suitable for nonlinear and nonGaussian estimation and tracking applications as encountered in the real world. A particle filter is essentially composed of three stages: prediction, update, and resampling. The prediction stage uses the system model to predict the state probability density function (PDF) forward from one measurement time to the next. Since the state is usually subject to unknown disturbances, prediction generally spreads the state PDF. The update operation takes the latest measurement to modify the prediction PDF using Bayes’ formula. A resampling step was introduced by Gordon et al. [31] in order to discard the particles with very low weights. Particle filters work by providing a Monte Carlo approximation to the PDF which can be easily updated to incorporate new information as it arrives.
A discrete time estimation problem is considered here. The state vector is denoted by whose temporal evolution is given by the state equation: where is the state transition function, and is the process noise with zero mean but not necessarily having a Gaussian distribution. In the ADSB radar system, the components of the state vector will be target locations. At each discrete time point an observation , related to the state vector, can be represented as follows: where is the measurement function and is the measurement noise with zero mean but not necessarily having a Gaussian distribution. is assumed to be uncorrelated with . In our application, is the range information and is the process to obtain the target ranges through the TOA technique.
A numerical integration method for the particle filter algorithm is listed below.(1)Particle generation: within the field of view (FOV), create particles and associated weights according to the uniform distribution.(2)Prediction: particles propagate according to evolution equation (1).(3)Measurement update: use the available measurements to compute the likelihood of each particle and update the weights of all particles with posteriori density. This is accomplished using where the final weights sum to one, namely, .(4)Systematic resampling: after a predetermined number of iterations, take samples with replacement from the current particle set based on .(5)Iteration: letting , repeat the process until desired estimation error is achieved.
For abruptly changing systems, Interacting Multiple Model (IMM) method [32–34] and Generalized PseudoBayesian (GPB) [35] are widely used in the target tracking literature. To reduce the computational complexity on the ADSB radar system, these approaches are not adopted. In fact, we do not even use the motion equation, since the dynamic model of the object cannot be obtained. In the prediction step, particles are simply scattered with a large variance that is able to cover abrupt motion changes.
Assume there are two targets around the own aircraft equipped with ADSB radar system. Each particle will be assigned two weights according to the state PDF. In the resampling step, these two weights will be merged by taking the largest weight for each particle. This process has twofold implication: (a) the peak selection from the posterior distribution and (b) partition of the particles to represent each target. The index of the maximum value out of the two weights for one particle will be used to mark which target this particle is associated with and is more likely to be close to. Since no two targets should and will share the same location, the larger weight of the partitioned particle is significant and retained while the smaller weight is discarded. After the resampling procedure, all particles whose both weights are low will be discarded and only those particles with at least one high weight will survive. Moreover, parts of the particles will be representing Target 1 if their first weights are clearly larger than the second ones. Figure 6 depicts the idea of such resampling mechanism for multiple targets. The resampling step removes particles that are improbable to be any of the targets, but it requires additional processing latency. Hence, we would like to intelligently pick up the right time to do the resampling step. The parameter , the effective sample size, can be used to measure the degeneracy of the algorithm [36, 37]. Once the particles spread out all over the place and only a certain amount of particles are meaningful, will become small. This is the right time to spend extracomputational expense to discard low weight particles. A good estimate of the effective sample size is given by , which ranges from 100% to 0%. The resampling step is performed when the number of meaningful particles is less than a predetermined threshold value, . This enables the particle filter algorithm to simply do the resampling at appropriate times. The numerical approximation of adaptive resampling for targets is given as follows, if is known.
The procedure is described below.(1)Set up , which represents the number of the meaningful particles, to be within .(2)Derive .(3)Calculate .(4)If , then take samples with replacement from the current particle set using .
4.2.2. Supplementary Particle Filter Algorithm
As the ADSB radar measurement is available once per second and the PF can be completed in the order of onetenth of a second, most of the time the system is idling and waiting for the new measurement. To improve the estimation performance, some sort of upsampling process is desirable and can be realized through piecewise constant interpolation between successive measurements. We name the algorithm using interpolated measurement on PF as the supplementary particle filter (SPF). Before the new measurement arrives, SPF performs iteratively using the current observation as if the target were static. Standard PF is vulnerable to sample impoverishment, because a finite number of particles are used to approximate a continuous distribution. The benefit of SPF over standard PF is that SPF minimizes the estimation error when the sample size is not sufficiently large. SPF provides improved estimation accuracy because particles will be redistributed to high likelihood areas over iterations although the same measurement information is reused. The sample resolution is essentially enhanced during SPF iterations, thereby improving the estimation performance.
SPF not only improves the estimation accuracy between successive measurements but also benefits further the estimation result when new measurement arrives, because more particles are already allocated to local mode of the posterior density. Through extensive simulation analysis, it has been noticed that PF takes many iterations to converge. While the ADSB message is available only once a second, a few dozen seconds may be needed to obtain accurate estimated target locations. However, for the pilots to react on imminent collisions, it is critical to improve the convergence rate and estimation accuracy in spite of the low measurement update rate. Overall, the convergence rate of SPF is about 23 times faster than PF, and the estimation error of SPF is about onehalf of the estimation error using standard PF. The gain is significant especially during the very first few measurements.
5. Simulation Results
In order to ascertain the effectiveness of the algorithm for the ADSB radar system, the transmitted signal is simulated according to the ADSB radar message format, which follows the requirements of the ADSB signal specification. The ADSB radar waveform, as shown in Figure 7, is generated by adding random biphase modulation into each bit of the ADSB signal. The signals received by multiple sensors are essentially the attenuated and delayed replica of the transmitted signal with measurement noise, . The sensor locations are not restricted, and in the simulation setup, the sensors are circularly positioned and 30 meters apart, using the maximum available distance on an airplane. The transmitted ADSB radar signal and the returned signals received by the four sensors are shown in Figure 8. The time differences between the transmitted and received signals are used to calculate the ranges from the target to each of the sensors. Since the ADSB radar signal is broadcast every second, a new measurement will also arrive approximately once per second.
Without prior estimation of motion model, additional positional uncertainty is incorporated in the state transition equation, as described in Section 4.2.1, to allow the particles to have the potential to move around and to compensate the unknown maneuver while searching for the best solution. As long as one or a few particles are able to “follow” the target, the calculated weights will be high, and subsequently, in the next resampling step, a large amount of particles will be drawn to the neighboring regions of the particles that have high weights. In the simulation, the positional uncertainty is set to 300 m after considering the relative aircraft speed and the finite amount of particles. By taking into account the detection range of 10 km and the computational load for multitarget tracking, the number of particles is chosen to be 10,000 and for the following three simulation scenarios: (1) constant velocity with random acceleration and direction noise, (2) basic flight maneuvers, and (3) multiple targets.
To test the tracking capability, we first assume a noncooperative aircraft having an RCS of 20 dBsm at a range of 10 km. According to the link budget analysis, as discussed in Section 3.2, the SNR is 11.1 dB. The target is moving at a constant velocity of 150 m/s with a Gaussian distributed acceleration and heading direction. Figure 9 depicts the tracking performance comparison of PF and SPF against the true target. In addition, the range errors of one trial and the rootmeansquare error (RMSE) of twenty Monte Carlo (MC) trials are plotted in Figures 10 and 11. The plots indicate that SPF estimated location converged to the true target location much faster, especially during the first few seconds. The range RMSEs of SPF and PF estimates at the 5th, 10th, and 15th seconds are listed in Table 3. Note that the difference between the RMSE values of the two tracking methods increases during the first few measurements, because each new SPF estimate intelligently exploits prior estimates closer to the true target location which is not done when using PF. After convergence, SPF also performs better than PF in terms of RMSE.

In the following simulation, the benefit of the additive positional uncertainty will be shown through tracking a maneuvering target even though the motion model is unknown. For a noncooperative and maneuvering target 5 km away, Figure 12 shows that both PF and SPF are capable to track the target with decent accuracy. Again, 20 MC trials are performed, and the range error for each iteration and the RMSE results are depicted in Figures 13 and 14, respectively. Compared to the target at 10 km away in the first example, the error in the position estimate rolls off much faster when the target is closer to the sensors. Nevertheless, the improvement of SPF over PF is still noticeable in terms of convergence rate and estimation accuracy.
A multitarget scenario is simulated as well to test the capability to track multiple surrounding targets simultaneously. Since this paper focuses on the noncooperative targets, we assume two noncooperative targets located relative to own aircraft at −4000 m, −3000 m and −5000 m, −5000 m, respectively. The relative speeds are 100 m/s and 150 m/s, and the trajectories are drawn in Figure 15. Even though the two target trajectories cross over, the proposed method is still able to provide satisfactory estimated trajectories. However, similar to PSR, the target identification is not included in the radar report and trajectory, for each aircraft needs to be constructed using data association methods, which is another discrete problem.
6. Conclusions
The ADSB system is considered as a cooperative surveillance technique to improve situational awareness while the ADSB radar system is an innovative addon implementation that exploits the standard ADSB system with only a few additional integrated devices for the detection of noncooperative targets. It is significant to point out the independence of the communication and radar capability, because the radar performance will not be deteriorated in the event of an incorrect ADSB message. The communication message is encoded using PPM while the radar signature is embedded in the sequence of the random phase modulation.
The ADSB radar signal and system design are discussed thoroughly in this paper, and computer simulation demonstrates that the proposed approach allows the ADSB radar system to track multiple noncooperative targets simultaneously. This research work opens a lot of opportunities for the ADSB system to incorporate airborne senseandavoid capability, which can be useful for the purpose of collision avoidance for resourcelimited UAVs [38, 39]. Moreover, the radar functionality of the ADSB radar system has the potential to serve as a backup surveillance in the event of loss of GNSS function and make the communication system spoof resistant to incorrect ADSB reports.
References
 SESAR Consortium, The ATM Target Concept: D3, Document Number DLM06120010200, 108 pages, 2007, http://www.eurocontrol.int/sesar/gallery/content/public/docs/DLM06120010200.pdf.
 Joint Planning and Development Office, Concept of Operations for the Next Generation Air Transportation System, Version 3.2, 153 pages, 2010, http://jpe.jpdo.gov/ee/docs/conops/NextGen_ConOps_v3_2.pdf.
 Study Group for the Future Air Traffic Systems, Collaborative Actions for Renovation of Air Traffic Systems (CARATS): Longterm Vision for the Future Air Traffic Systems, 56 pages, 2010, http://www.mlit.go.jp/common/000128185.pdf.
 P. Brooker, “SESAR and NextGen: investing in new paradigms,” Journal of Navigation, vol. 61, no. 2, pp. 195–208, 2008. View at: Publisher Site  Google Scholar
 G. Wright, “NAV Canada implements ADSB,” in Proceedings of the Integrated Communications, Navigation and Surveillance Conference (ICNS '09), pp. 1–9, Arlington, Va, USA, 2009. View at: Publisher Site  Google Scholar
 L. Purton, H. Abbass, and S. Alam, “Identification of ADSB system vulnerabilities and threats,” in Proceedings of the 33rd Australasian Transport Research Forum (ATRF '10), pp. 1–16, Canberra, Australia, September 2010, http://www.patrec.org/web_docs/atrf/papers/2010/1916_222%20%20Purton%20Abbass%20Alam.pdf. View at: Google Scholar
 J. A. Besada, J. Garcia, G. de Miguel, F. J. Jimenez, G. Gavin, and J. R. Casar, “Data fusion algorithms based on radar and ADS measurements for ATC application,” in Proceedings of the IEEE International Radar Conference, pp. 98–103, Alexandria, Va, USA, May 2000. View at: Google Scholar
 J. L. R. Da Silva, J. F. B. Brancalion, and D. Fernandes, “Data fusion techniques applied to scenarios including ADSB and radar sensors for air traffic control,” in Proceedings of the 12th International Conference on Information Fusion (FUSION '09), pp. 1481–1488, Seattle, Wash, USA, July 2009. View at: Google Scholar
 C. Mensing, S. Sand, and A. Dammann, “Hybrid data fusion and tracking for positioning with GNSS and 3GPPLTE,” International Journal of Navigation and Observation, vol. 2010, Article ID 812945, 12 pages, 2010. View at: Publisher Site  Google Scholar
 A. Grant, P. Williams, N. Ward, and S. Basker, “GPS jamming and the impact on maritime navigation,” Journal of Navigation, vol. 62, no. 2, pp. 173–187, 2009. View at: Publisher Site  Google Scholar
 Y. Zhang and J. Qiao, “ADSB radar system,” US Patent 7414567, Intelligent Automation, Rockville, Md, USA, 2008. View at: Google Scholar
 M.S. Huang, R. M. Narayanan, and A. Feinberg, “Multiple targets estimation and tracking for ADSB radar system,” in Proceedings of the IEEE/AIAA 27th Digital Avionics Systems Conference (DASC '08), pp. 3C11–3C110, St. Paul, Minn, USA, 2008. View at: Publisher Site  Google Scholar
 C. P. Lai, Y. J. Ren, and C. Lin, “ADSB based collision avoidance radar for unmanned aerial vehicles,” in 2009 IEEE MTTS International Microwave Symposium (IMS '09), pp. 85–88, Boston, Mass, USA, June 2009. View at: Publisher Site  Google Scholar
 S. Wang, Y. Zhang, A. Huston, Z. Li, M. Mallo, and J. Fagan, “Phase modulated waveforms for transponderbased radar sensing: signal optimization and experiments,” IEEE Transactions on Aerospace and Electronic Systems, vol. 47, no. 3, pp. 1733–1753, 2011. View at: Publisher Site  Google Scholar
 M.S. Huang and R. M. Narayanan, “Noncooperative collision avoidance concept for Unmanned Aircraft System using satellitebased radar and radio communication,” in Proceedings of the 30th Digital Avionics Systems Conference—Closing the Generation Gap: Increasing Capability for Flight Operations among Legacy, Modern and Uninhabited Aircraft (DASC '11), pp. 5C21–5C29, Seattle, Wash, USA, 2011. View at: Publisher Site  Google Scholar
 R. Barhydt and A. W. Warren, “Development of intent information changes to revised minimum aviation system performance standards for automatic dependent surveillance broadcast (RTCA/DO242A),” NASA Technical Memorandum TM 2002211669, 2002. View at: Google Scholar
 K. Owusu and G. Dunstone, “Development of cockpit display of traffic information (CDTI),” in Proceedings of the 3rd Meeting of Automatic Dependent Surveillance—Broadcast (ADSB) Study and Implementation Task Force (ADSB TF '05), Bangkok, Thailand, March 2005. View at: Google Scholar
 S. R. J. Axelsson, “Noise radar using random phase and frequency modulation,” IEEE Transactions on Geoscience and Remote Sensing, vol. 42, no. 11, pp. 2370–2384, 2004. View at: Google Scholar
 R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986. View at: Google Scholar
 T. Akiyama, T. Yamaoka, and N. Hamada, “2D angle estimation using the constrained MUSIC with circular array,” in Proceedings of the 5th International Symposium on Signal Processing and Its Applications, pp. 877–880, Brisbane, Australia, August 1999. View at: Google Scholar
 Z. Saidi and S. Bourennane, “Cumulantbased coherent signal subspace method for bearing and range estimation,” EURASIP Journal on Advances in Signal Processing, vol. 2007, Article ID 84576, 2007. View at: Publisher Site  Google Scholar
 H. Hung and M. Kaveh, “Focussing matrices for coherent signalsubspace processing,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 36, no. 8, pp. 1272–1281, 1988. View at: Publisher Site  Google Scholar
 P. Fearnhead, Sequential Monte Carlo methods in filter theory [Ph.D. thesis], University of Oxford, Oxford, UK, 1998.
 Y. BarShalom and T. E. Fortmann, Tracking and Data Association, Academic Press, New York, NY, USA, 1988.
 I. Hwang, H. Balakrishnan, K. Roy, and C. Tomlin, “Multipletarget tracking and identity management with application to aircraft tracking,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 3, pp. 641–653, 2007. View at: Publisher Site  Google Scholar
 F. Gustafsson, F. Gunnarsson, N. Bergman et al., “Particle filters for positioning, navigation, and tracking,” IEEE Transactions on Signal Processing, vol. 50, no. 2, pp. 425–437, 2002. View at: Publisher Site  Google Scholar
 M. S. Arulampalam, S. Maskell, N. Gordon, and T. Clapp, “A tutorial on particle filters for online nonlinear/nonGaussian Bayesian tracking,” IEEE Transactions on Signal Processing, vol. 50, no. 2, pp. 174–188, 2002. View at: Publisher Site  Google Scholar
 B. Ristic, S. Arulampalam, and N. Gordon, Beyond the Kalman Filter: Particle Filters for Tracking Applications, Artech House, Norwood, Mass, USA, 2004.
 N. Gordon, D. Salmond, and C. Ewing, “Bayesian state estimation for tracking and guidance using the bootstrap filter,” Journal of Guidance, Control, and Dynamics, vol. 18, no. 6, pp. 1434–1443, 1995. View at: Google Scholar
 M. Isard and A. Blake, “CONDENSATION—conditional density propagation for visual tracking,” International Journal of Computer Vision, vol. 29, no. 1, pp. 5–28, 1998. View at: Google Scholar
 N. J. Gordon, D. J. Salmond, and A. F. M. Smith, “Novel approach to nonlinear/nongaussian Bayesian state estimation,” IEE Proceedings, Part F, vol. 140, no. 2, pp. 107–113, 1993. View at: Google Scholar
 E. Mazor, A. Averbuch, Y. BarShalom, and J. Dayan, “Interacting multiple model methods in target tracking: a survey,” IEEE Transactions on Aerospace and Electronic Systems, vol. 34, no. 1, pp. 103–123, 1998. View at: Google Scholar
 Z. Khan, T. Balch, and F. Dellaert, “An MCMCbased particle filter for tracking multiple interacting targets,” in Proceedings of the 8th European Conference on Computer Vision (ECCV '04), vol. 3024 of Lecture Notes in Computer Science, pp. 279–290, 2004. View at: Google Scholar
 H. A. P. Blom, “An efficient filter for abruptly changing systems,” in Proceedings of the 23rd IEEE Conference on Decision and Control, pp. 656–658, Las Vegas, Nev, USA, December 1984. View at: Google Scholar
 J. K. Tugnait, “Detection and estimation for abruptly changing systems,” in Proceedings of the 20th IEEE Conference on Decision and Control, pp. 1357–1362, San Diego, Calif, USA, December 1981. View at: Google Scholar
 J. S. Liu, “Metropolized independent sampling with comparisons to rejection sampling and importance sampling,” Statistics and Computing, vol. 6, no. 2, pp. 113–119, 1996. View at: Google Scholar
 A. Kong, J. S. Liu, and W. H. Wong, “Sequential imputations and Bayesian missing data problems,” Journal of the American Statistical Association, vol. 89, no. 425, pp. 278–288, 1994. View at: Google Scholar
 B. A. Kumar and D. Ghose, “Radarassisted collision avoidance/guidance strategy for planar flight,” IEEE Transactions on Aerospace and Electronic Systems, vol. 37, no. 1, pp. 77–90, 2001. View at: Publisher Site  Google Scholar
 S. C. Han, H. Bang, and C. S. Yoo, “Proportional navigationbased collision avoidance for UAVs,” International Journal of Control, Automation and Systems, vol. 7, no. 4, pp. 553–565, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 MingShih Huang 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.