# Tracking of Noncooperative Airborne Targets Using ADS-B Signal and Radar Sensing

Authors: Ming-Shih Huang; Ram M. Narayanan; Yan Zhang; Arthur Feinberg
Journal: International Journal of Aerospace Engineering
 (2013)

---

## Abstract

As the Automatic Dependent Surveillance-Broadcast (ADS-B) system has gained wide acceptance, additional exploitations of the radioed satellite-based information are topics of current interest. One such opportunity includes the augmentation of the communication ADS-B 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 ADS-B radar system using radar ranging techniques and particle filter algorithms. In addition, the low-update-rate measurement due to the ADS-B system specification is discussed in order to provide acceptable estimation results. Simulation results show satisfactory tracking capability up to several kilometers with acceptable accuracy.


---

## Body

## 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 self-reliance. 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 altitude-reporting transponders. However, TCAS is ineffective for noncooperative aircraft and also in low-altitude 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 ADS-B, 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 ADS-B’s wide acceptance in Europe and USA, NAV CANADA commenced operational application of ADS-B as a means of providing aircraft surveillance information to Air Traffic Controllers (ATC) [5], and AirServices Australia commissioned the ADS-B Upper Airspace Project (UAP), providing ADS-B coverage across the whole continent [6]. In addition, ADS-B is being used as the key solution for Unmanned Aerial System (UAS) integration in the National Airspace (NAS). When properly equipped with ADS-B, both pilots and controllers will see the same real-time displays of air traffic, substantially improving safety and minimizing collision probability.However, before the ADS-B implementation and operation are fully complete, there will be a transition period involving coexistence of ADS-B equipped and nonequipped aircraft. Similar to TCAS, the ADS-B 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 ADS-B 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 ADS-B system.Since the emergence of ADS-B concept, some researchers have considered the utilization of existing and installed infrastructure of the surveillance radar to combine with the satellite-based ADS-B system within the perspective of ATC. More interestingly, the use of the ADS-B signal itself to detect noncooperative targets from the ADS-B 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 ADS-B radar system lies in that the system insightfully exploits the ADS-B 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 radar-communication system specifications and problems of low-update-rate measurement due to the ADS-B system requirement have not been fully discussed. This paper provides details of the ADS-B radar system design, including the link power budget, and RF signal interference analysis. Furthermore, multitarget tracking using particle filter for the ADS-B radar system is proposed. In addition, the slow convergence in particle filter algorithm and the measurement at low-update-rate, owing to the ADS-B signal broadcast rate of once per second, is addressed. It is important to note that in the ADS-B radar system, the radar capability is independent of the ADS-B message data integrity. Therefore, the coexistence of ADS-B 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. Section2 briefly reviews the evolution of air surveillance systems. The ADS-B radar system architecture and link budget are provided in Section 3. Then, the multitarget estimation algorithm for ADS-B radar system is proposed and discussed in Section 4. Computer simulations that demonstrate the functionality of ADS-B 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 transponder-equipped 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, ground-based 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 early-to mid-1950s to serve as a last resort collision avoidance safety net. TCAS operates similarly to the ground-based SSR but independently interrogates surrounding aircraft on a 1030-MHz radio channel. The pilot will be alerted to the presence of the intruding aircraft replying to the interrogation via 1090-MHz radio frequency. Current generation TCAS II, jointly developed by the US Radio Technical Commission for Aeronautics (DO-185B) and European Organization for Civil Aviation Equipment (ED-143), 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 mid-air collisions and has been very successful since its introduction. However, the major concern with regard to either TCAS or ADS-B is that they are not required all the time, and ADS-B out is not yet mandated in most countries. Aircraft equipped with TCAS and/or ADS-B 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 transponder-required airspace.


### 2.3. ADS-B

ADS-B is redefining the paradigm of communication, navigation, and surveillance in ATM. An ADS-B equipped aircraft determines its own position using GNSS and periodically broadcasts its four-dimensional 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 ADS-B system and potential ground stations without expectation of an acknowledgement or reply. One of the most significant advantages of the ADS-B 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 ADS-B surveillance information. ADS-B system provides accurate information and improves situational awareness. Moreover, ADS-B enables a shift from a centralized, ground-based 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. ADS-B Radar

While the ADS-B system can only be considered as a decentralized cooperative surveillance ATM technique, the ADS-B radar system deals with both cooperative and noncooperative targets. Figure1 depicts the general working principle for ADS-B radar concept, as well as the comparison between the ADS-B system and the ADS-B radar system. Noncooperative targets, such as unmanned aerial vehicles (UAVs) and private jets, are blind to the ADS-B systems when not equipped with the ADS-B 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 ADS-B radar system is introduced as an insightful modification to the standard ADS-B system in the interests of a safe backup service without significantly enlarging the volume of ADS-B equipment. Since the ADS-B systems broadcast signals periodically, the reflective ADS-B 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 ADS-B equipped targets, the ADS-B radar system can suppress the creation of tracks based on ADS-B messages that contain intentionally incorrect position information.Surveillance principles and concepts for ADS-B and ADS-B radar. (a) ADS-B: communication between equipped aircraft and ATC. (b) ADS-B radar: able to detect both equipped and nonequipped aircraft.
(a)
(b)


## 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 transponder-equipped 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, ground-based 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 early-to mid-1950s to serve as a last resort collision avoidance safety net. TCAS operates similarly to the ground-based SSR but independently interrogates surrounding aircraft on a 1030-MHz radio channel. The pilot will be alerted to the presence of the intruding aircraft replying to the interrogation via 1090-MHz radio frequency. Current generation TCAS II, jointly developed by the US Radio Technical Commission for Aeronautics (DO-185B) and European Organization for Civil Aviation Equipment (ED-143), 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 mid-air collisions and has been very successful since its introduction. However, the major concern with regard to either TCAS or ADS-B is that they are not required all the time, and ADS-B out is not yet mandated in most countries. Aircraft equipped with TCAS and/or ADS-B 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 transponder-required airspace.


## 2.3. ADS-B

ADS-B is redefining the paradigm of communication, navigation, and surveillance in ATM. An ADS-B equipped aircraft determines its own position using GNSS and periodically broadcasts its four-dimensional 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 ADS-B system and potential ground stations without expectation of an acknowledgement or reply. One of the most significant advantages of the ADS-B 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 ADS-B surveillance information. ADS-B system provides accurate information and improves situational awareness. Moreover, ADS-B enables a shift from a centralized, ground-based 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. ADS-B Radar

While the ADS-B system can only be considered as a decentralized cooperative surveillance ATM technique, the ADS-B radar system deals with both cooperative and noncooperative targets. Figure1 depicts the general working principle for ADS-B radar concept, as well as the comparison between the ADS-B system and the ADS-B radar system. Noncooperative targets, such as unmanned aerial vehicles (UAVs) and private jets, are blind to the ADS-B systems when not equipped with the ADS-B 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 ADS-B radar system is introduced as an insightful modification to the standard ADS-B system in the interests of a safe backup service without significantly enlarging the volume of ADS-B equipment. Since the ADS-B systems broadcast signals periodically, the reflective ADS-B 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 ADS-B equipped targets, the ADS-B radar system can suppress the creation of tracks based on ADS-B messages that contain intentionally incorrect position information.Surveillance principles and concepts for ADS-B and ADS-B radar. (a) ADS-B: communication between equipped aircraft and ATC. (b) ADS-B radar: able to detect both equipped and nonequipped aircraft.
(a)
(b)


## 3. System Design Specification

A conceptual architecture of the ADS-B radar system is depicted in Figure2. 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 ADS-B message, as well as the target estimation method in radar processing. The fundamental principle of the ADS-B radar system is to use reflective ADS-B signals with only four additional L-band antennas on the standard ADS-B system to perform multitarget estimation and tracking for both cooperative and noncooperative targets. Because the ADS-B receivers only respond to the envelope of the incoming signal, phase modulation will not affect the normal ADS-B signal decoding for other aircraft. As a result, the standard ADS-B system will still be able to correctly decode the phase-modulated signal to obtain the communication message. The ADS-B radar system retains the communication capability and receives the standard ADS-B messages from neighboring aircraft, just as the standard ADS-B system does. Once per second, it also broadcasts the phase modulated ADS-B signal. After sending out the electromagnetic wave, the receiving antennas capture the returned signal reflected from surrounding metallic objects, and cross-correlation is performed between the transmitted burst and the returned echoes. Both the ADS-B In information and the estimated location from the radar report can be fed into Cockpit Display of Traffic Information (CDTI) [17].Figure 2
Conceptual architecture of ADS-B radar system.The design consideration is focused on the feasibility of utilizing the returned ADS-B signal as a radar echo in terms of the signal-to-noise ratio and radar performance. In addition, the 500 W peak transmit power, the 120μs pulse duration specifications of the 1090 MHz ADS-B 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 ADS-B 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 ADS-B 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 ADS-B radar design is to exam the feasibility of treating the ADS-B signal as a radar signal.In order to improve the radar performance, the major modification to the ADS-B signal is random biphase modulation introduced to each bit. Figure3 shows the simulated ADS-B waveform and the proposed ADS-B radar waveform. Both signals are of the same duration and have identical digital messages, except that the ADS-B radar signal has a random 180° phase shift. By randomizing the transmit signal, the matched filtering operation can be performed by cross-correlating the reflected signal with a time-delayed replica of the transmit waveform. The biphase modulation renders both positive and negative products, forcing the autocorrelation to be statistically zero except for zero time-lag. 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 ADS-B radar signal and the ADS-B signal, as depicted in Figure 4. It is clear that the autocorrelation of the ADS-B radar signal outperforms that of the original modulation-free ADS-B signal and shows an improved correlation peak-to-noise ratio.Table 1
Effects of random biphase modulation on correlation results.
Messages
Modulation
PPM
PPM with random biphase PM
Same phase(no phase shift)
Inverted phase(180° phase shift)
(1, 1) or (0, 0)
Positivesum-product
Positive
sum-product
Negative
sum-product
(1, 0) or (0, 1)
Zerosum-product
Zerosum-product
Zerosum-productFigure 3
Illustration of the ADS-B waveform (blue) versus the ADS-B radar waveform (red).Figure 4
Autocorrelation of standard ADS-B signal and autocorrelation of the ADS-B radar signal.Since it is very likely to have other ADS-B signals in the same channel, it is critical to analyze the interference of the ADS-B signals and possibly other ADS-B radar signals as well. Figure5 shows that a typical ADS-B radar signal is uncorrelated with both a standard ADS-B signal as well as another independent ADS-B radar signal, thereby indicating that the on-board ADS-B radar receiver will not be affected by standard ADS-B or ADS-B radar transmissions from other aircraft in the vicinity.Figure 5
Cross-correlation of a standard phase-modulated ADS-B signal with another standard ADS-B signal and another phase-modulated ADS-B radar signal.Because the modulated signal has wider bandwidth (the resultant waveform still complies with the ADS-B standard), the range resolution is improved and the range side lobes are also further suppressed [18]. The range resolution of 75 meters for the ADS-B radar system is sufficient for most collision avoidance applications.


### 3.2. System Configuration and Link Budget Analysis

The ADS-B radar system is composed of the following components: (a) ADS-B transceiver, (b) ADS-B 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 ADS-B 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:(1)Pr=PtGtGrλ2σ(4π)3R4,
where Pr and Pt are the received and transmitted power, respectively; Gt and Gr are the receiving and transmitting antennas gain, respectively; λ is the signal wavelength; σ is the target’s radar cross-section (RCS); R 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 signal-to-noise 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.Table 2
Link budget analysis.
Receiver noise floor = −110.9 dBm
Received signal power = −105.9 dBm
Noise figure
3 dB
Peak transmit power
+57 dBm
Bandwidth
1 MHz baseband
Processing gain(545 samples coherently integrated)
27.4 dB
Antenna gain(omnidirectional)
0 dB
Assumed RCS
0 dBsm
1/(4π)3
−33.0 dB
1/(Range)4 (at 4.5 km)
−146.1 dBm−4
Square of wavelength
−11.2 dBsm
Signal-to-noise ratio, SNR = 5.0 dB


## 3.1. Signal Waveform and Interference Analysis

Each of the ADS-B 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 ADS-B 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 ADS-B radar design is to exam the feasibility of treating the ADS-B signal as a radar signal.In order to improve the radar performance, the major modification to the ADS-B signal is random biphase modulation introduced to each bit. Figure3 shows the simulated ADS-B waveform and the proposed ADS-B radar waveform. Both signals are of the same duration and have identical digital messages, except that the ADS-B radar signal has a random 180° phase shift. By randomizing the transmit signal, the matched filtering operation can be performed by cross-correlating the reflected signal with a time-delayed replica of the transmit waveform. The biphase modulation renders both positive and negative products, forcing the autocorrelation to be statistically zero except for zero time-lag. 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 ADS-B radar signal and the ADS-B signal, as depicted in Figure 4. It is clear that the autocorrelation of the ADS-B radar signal outperforms that of the original modulation-free ADS-B signal and shows an improved correlation peak-to-noise ratio.Table 1
Effects of random biphase modulation on correlation results.
Messages
Modulation
PPM
PPM with random biphase PM
Same phase(no phase shift)
Inverted phase(180° phase shift)
(1, 1) or (0, 0)
Positivesum-product
Positive
sum-product
Negative
sum-product
(1, 0) or (0, 1)
Zerosum-product
Zerosum-product
Zerosum-productFigure 3
Illustration of the ADS-B waveform (blue) versus the ADS-B radar waveform (red).Figure 4
Autocorrelation of standard ADS-B signal and autocorrelation of the ADS-B radar signal.Since it is very likely to have other ADS-B signals in the same channel, it is critical to analyze the interference of the ADS-B signals and possibly other ADS-B radar signals as well. Figure5 shows that a typical ADS-B radar signal is uncorrelated with both a standard ADS-B signal as well as another independent ADS-B radar signal, thereby indicating that the on-board ADS-B radar receiver will not be affected by standard ADS-B or ADS-B radar transmissions from other aircraft in the vicinity.Figure 5
Cross-correlation of a standard phase-modulated ADS-B signal with another standard ADS-B signal and another phase-modulated ADS-B radar signal.Because the modulated signal has wider bandwidth (the resultant waveform still complies with the ADS-B standard), the range resolution is improved and the range side lobes are also further suppressed [18]. The range resolution of 75 meters for the ADS-B radar system is sufficient for most collision avoidance applications.


## 3.2. System Configuration and Link Budget Analysis

The ADS-B radar system is composed of the following components: (a) ADS-B transceiver, (b) ADS-B 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 ADS-B 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:(1)Pr=PtGtGrλ2σ(4π)3R4,
where Pr and Pt are the received and transmitted power, respectively; Gt and Gr are the receiving and transmitting antennas gain, respectively; λ is the signal wavelength; σ is the target’s radar cross-section (RCS); R 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 signal-to-noise 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.Table 2
Link budget analysis.
Receiver noise floor = −110.9 dBm
Received signal power = −105.9 dBm
Noise figure
3 dB
Peak transmit power
+57 dBm
Bandwidth
1 MHz baseband
Processing gain(545 samples coherently integrated)
27.4 dB
Antenna gain(omnidirectional)
0 dB
Assumed RCS
0 dBsm
1/(4π)3
−33.0 dB
1/(Range)4 (at 4.5 km)
−146.1 dBm−4
Square of wavelength
−11.2 dBsm
Signal-to-noise ratio, SNR = 5.0 dB


## 4. Multitarget Estimation and Tracking Algorithm

The Multiple Signal Classification (MUSIC) technique [19] was initially employed for the ADS-B radar to perform detection and estimation of single target angle-of-arrival (AOA). Due to the nature of the ADS-B 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 well-known inherent constraint that the number of the receivers has to be larger than that of the targets. With fixed number of N antennas, the maximum number of detectable targets using the MUSIC algorithm is N-1. While it is desirable to minimize the cost and amount of antennas used in the ADS-B 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], cumulant-based coherent signal subspace method [21], and focusing matrix for coherent signal subspace processing (for wide-band signals) [22], these approaches do not apply to the scenario that the ADS-B radar system encounters owing to either unavailable prior information or limited signal bandwidth from the ADS-B 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 time-of-arrival (TOA) and the particle filter approaches for multitarget estimation and tracking for the ADS-B 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 ADS-B messages with radar reports are described in [7, 8].


### 4.1. Multitarget Detection and Ranging through TOA

Lets(t) represent the transmitted signal. The distance between the target and the own aircraft can be obtained by cross-correlating the delayed transmitted signal, s(t-τd), with the received signal, s(t-τr), where τd is the internal delay and τr=2R/c is the round-trip time to the target. The target range is R, and the speed of light is c. In theory, the cross-correlation peak occurs when τd=τr, 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 ADS-B radar system, s(t) is the phase-modulated ADS-B signal. By continuously cross-correlating 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 ADS-B radar system. Furthermore, the problem of slow convergence in particle filter and measurements at low-update-rate 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 byxt whose temporal evolution is given by the state equation:
(2)xt=f(xt-1)+vt-1,
where f is the state transition function, and vt is the process noise with zero mean but not necessarily having a Gaussian distribution. In the ADS-B radar system, the components of the state vector will be target locations. At each discrete time point an observation yt, related to the state vector, can be represented as follows:
(3)yt=h(xt)+nt,
where h is the measurement function and nt is the measurement noise with zero mean but not necessarily having a Gaussian distribution. nt is assumed to be uncorrelated with vt. In our application, yt is the range information and h 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 N particles and associated weights (xt-1n,w(xt-1n))n=1,…,N 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 a posteriori density. This is accomplished using
(4)w(xti)∝w(xt-1i)p(yt∣xti),where the final weights sum to one, namely,∑i=1Nw(xti)=1.(4)
Systematic resampling: after a predetermined number of iterations, take N samples with replacement from the current particle set based on w(xti).(5)
Iteration: letting t=t+1, repeat the process until desired estimation error is achieved.For abruptly changing systems, Interacting Multiple Model (IMM) method [32–34] and Generalized Pseudo-Bayesian (GPB) [35] are widely used in the target tracking literature. To reduce the computational complexity on the ADS-B 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 ADS-B 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 two-fold 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. Figure6 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 Neff, 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, Neff 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 N~eff=[1+N-1∑i=1N[N×w(xti)-1]2]-1, which ranges from 100% to 0%. The resampling step is performed when the number of meaningful particles is less than a predetermined threshold value, Nthreshold. This enables the particle filter algorithm to simply do the resampling at appropriate times. The numerical approximation of adaptive resampling for m targets is given as follows, if m is known.Figure 6
Resampling mechanism for multiple targets.The procedure is described below.(1)
Set upNthreshold, which represents the number of the meaningful particles, to be within [0,1].(2)
Derivewnew={max(wki)∣k=1,…,m,i=1,…,N}.(3)
CalculateN~eff=[1+N-1∑i=1N[N×wki-1]2]-1,k=1,…,m.(4)
Ifmin{N~eff(w1),…,N~eff(wm)}<Nthreshold, then take N samples with replacement from the current particle set using wnew.


#### 4.2.2. Supplementary Particle Filter Algorithm

As the ADS-B radar measurement is available once per second and the PF can be completed in the order of one-tenth 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 ADS-B 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 2-3 times faster than PF, and the estimation error of SPF is about one-half of the estimation error using standard PF. The gain is significant especially during the very first few measurements.


## 4.1. Multitarget Detection and Ranging through TOA

Lets(t) represent the transmitted signal. The distance between the target and the own aircraft can be obtained by cross-correlating the delayed transmitted signal, s(t-τd), with the received signal, s(t-τr), where τd is the internal delay and τr=2R/c is the round-trip time to the target. The target range is R, and the speed of light is c. In theory, the cross-correlation peak occurs when τd=τr, 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 ADS-B radar system, s(t) is the phase-modulated ADS-B signal. By continuously cross-correlating 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 ADS-B radar system. Furthermore, the problem of slow convergence in particle filter and measurements at low-update-rate 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 byxt whose temporal evolution is given by the state equation:
(2)xt=f(xt-1)+vt-1,
where f is the state transition function, and vt is the process noise with zero mean but not necessarily having a Gaussian distribution. In the ADS-B radar system, the components of the state vector will be target locations. At each discrete time point an observation yt, related to the state vector, can be represented as follows:
(3)yt=h(xt)+nt,
where h is the measurement function and nt is the measurement noise with zero mean but not necessarily having a Gaussian distribution. nt is assumed to be uncorrelated with vt. In our application, yt is the range information and h 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 N particles and associated weights (xt-1n,w(xt-1n))n=1,…,N 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 a posteriori density. This is accomplished using
(4)w(xti)∝w(xt-1i)p(yt∣xti),where the final weights sum to one, namely,∑i=1Nw(xti)=1.(4)
Systematic resampling: after a predetermined number of iterations, take N samples with replacement from the current particle set based on w(xti).(5)
Iteration: letting t=t+1, repeat the process until desired estimation error is achieved.For abruptly changing systems, Interacting Multiple Model (IMM) method [32–34] and Generalized Pseudo-Bayesian (GPB) [35] are widely used in the target tracking literature. To reduce the computational complexity on the ADS-B 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 ADS-B 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 two-fold 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. Figure6 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 Neff, 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, Neff 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 N~eff=[1+N-1∑i=1N[N×w(xti)-1]2]-1, which ranges from 100% to 0%. The resampling step is performed when the number of meaningful particles is less than a predetermined threshold value, Nthreshold. This enables the particle filter algorithm to simply do the resampling at appropriate times. The numerical approximation of adaptive resampling for m targets is given as follows, if m is known.Figure 6
Resampling mechanism for multiple targets.The procedure is described below.(1)
Set upNthreshold, which represents the number of the meaningful particles, to be within [0,1].(2)
Derivewnew={max(wki)∣k=1,…,m,i=1,…,N}.(3)
CalculateN~eff=[1+N-1∑i=1N[N×wki-1]2]-1,k=1,…,m.(4)
Ifmin{N~eff(w1),…,N~eff(wm)}<Nthreshold, then take N samples with replacement from the current particle set using wnew.


### 4.2.2. Supplementary Particle Filter Algorithm

As the ADS-B radar measurement is available once per second and the PF can be completed in the order of one-tenth 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 ADS-B 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 2-3 times faster than PF, and the estimation error of SPF is about one-half of the estimation error using standard PF. The gain is significant especially during the very first few measurements.


## 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 byxt whose temporal evolution is given by the state equation:
(2)xt=f(xt-1)+vt-1,
where f is the state transition function, and vt is the process noise with zero mean but not necessarily having a Gaussian distribution. In the ADS-B radar system, the components of the state vector will be target locations. At each discrete time point an observation yt, related to the state vector, can be represented as follows:
(3)yt=h(xt)+nt,
where h is the measurement function and nt is the measurement noise with zero mean but not necessarily having a Gaussian distribution. nt is assumed to be uncorrelated with vt. In our application, yt is the range information and h 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 N particles and associated weights (xt-1n,w(xt-1n))n=1,…,N 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 a posteriori density. This is accomplished using
(4)w(xti)∝w(xt-1i)p(yt∣xti),where the final weights sum to one, namely,∑i=1Nw(xti)=1.(4)
Systematic resampling: after a predetermined number of iterations, take N samples with replacement from the current particle set based on w(xti).(5)
Iteration: letting t=t+1, repeat the process until desired estimation error is achieved.For abruptly changing systems, Interacting Multiple Model (IMM) method [32–34] and Generalized Pseudo-Bayesian (GPB) [35] are widely used in the target tracking literature. To reduce the computational complexity on the ADS-B 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 ADS-B 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 two-fold 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. Figure6 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 Neff, 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, Neff 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 N~eff=[1+N-1∑i=1N[N×w(xti)-1]2]-1, which ranges from 100% to 0%. The resampling step is performed when the number of meaningful particles is less than a predetermined threshold value, Nthreshold. This enables the particle filter algorithm to simply do the resampling at appropriate times. The numerical approximation of adaptive resampling for m targets is given as follows, if m is known.Figure 6
Resampling mechanism for multiple targets.The procedure is described below.(1)
Set upNthreshold, which represents the number of the meaningful particles, to be within [0,1].(2)
Derivewnew={max(wki)∣k=1,…,m,i=1,…,N}.(3)
CalculateN~eff=[1+N-1∑i=1N[N×wki-1]2]-1,k=1,…,m.(4)
Ifmin{N~eff(w1),…,N~eff(wm)}<Nthreshold, then take N samples with replacement from the current particle set using wnew.


## 4.2.2. Supplementary Particle Filter Algorithm

As the ADS-B radar measurement is available once per second and the PF can be completed in the order of one-tenth 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 ADS-B 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 2-3 times faster than PF, and the estimation error of SPF is about one-half 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 ADS-B radar system, the transmitted signal is simulated according to the ADS-B radar message format, which follows the requirements of the ADS-B signal specification. The ADS-B radar waveform, as shown in Figure7, is generated by adding random biphase modulation into each bit of the ADS-B signal. The signals received by multiple sensors are essentially the attenuated and delayed replica of the transmitted signal with measurement noise, σn. 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 ADS-B 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 ADS-B radar signal is broadcast every second, a new measurement will also arrive approximately once per second.Figure 7
Transmitted ADS-B radar signal waveform.Figure 8
Transmitted ADS-B radar signal and received signals from four sensors.Without prior estimation of motion model, additional positional uncertainty is incorporated in the state transition equation, as described in Section4.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 Nthreshold=5% 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 non-cooperative aircraft having an RCS of 20 dBsm at a range of 10 km. According to the link budget analysis, as discussed in Section3.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 root-mean-square 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.Table 3
RMSE comparison (20 MC trials).
RMSE after 5 seconds
RMSE after 10 seconds
RMSE after 15 seconds
PF
885.2 m
488.7 m
289.6 m
SPF
212.9 m
146.5 m
45.33 mFigure 9
Tracking trajectories of PF and SPF methods against true target (20 MC trials).Figure 10
Range errors during each iteration (one trial).Figure 11
RMSE for a target with constant velocity, as well as Gaussian distributed acceleration and heading direction (20 MC trials).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, Figure12 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.Figure 12
Tracking performance for a maneuvering target (20 MC trials).Figure 13
Range errors during each iteration (one trial).Figure 14
RMSE for a maneuvering target (20 MC trials).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 Figure15. 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.Figure 15
Tracking performance for multiple targets (20 MC trials).


## 6. Conclusions

The ADS-B system is considered as a cooperative surveillance technique to improve situational awareness while the ADS-B radar system is an innovative add-on implementation that exploits the standard ADS-B 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 ADS-B message. The communication message is encoded using PPM while the radar signature is embedded in the sequence of the random phase modulation.The ADS-B radar signal and system design are discussed thoroughly in this paper, and computer simulation demonstrates that the proposed approach allows the ADS-B radar system to track multiple noncooperative targets simultaneously. This research work opens a lot of opportunities for the ADS-B system to incorporate airborne sense-and-avoid capability, which can be useful for the purpose of collision avoidance for resource-limited UAVs [38, 39]. Moreover, the radar functionality of the ADS-B 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 ADS-B reports.



---

*Source: 521630-2013-02-24.xml*