MATLAB Examples

Airplane Tracking Using ADS-B Signals and ADALM-PLUTO Radio

This example shows you how to track planes by processing Automatic Dependent Surveillance-Broadcast (ADS-B) signals using MATLAB® and Communications System Toolbox™. You can either use captured signals or receive signals in real time using the Communications System Toolbox Support Package for ADALM-PLUTO Radio. The example can show the tracked planes on a map, if you have the Mapping Toolbox™.

Contents

Required Hardware and Software

To run this example using captured signals, you need the following software:

To receive signals in real time, you also need the following hardware:

  • ADALM-PLUTO radio

and the following software

For a full list of Communications System Toolbox supported SDR platforms, refer to Supported Hardware section of Software Defined Radio (SDR) discovery page.

Background

ADS-B is a cooperative surveillance technology for tracking aircraft. This technology enables an aircraft to periodically broadcast its position information (altitude, GPS coordinates, heading, etc.) using the Mode-S signaling scheme.

Mode-S is a type of aviation transponder interrogation mode. When an aircraft receives an interrogation request, it sends back a transponder's squawk code. This is referred to as Mode 3A. Mode-S (Select) is another type of interrogation mode that is designed to help avoid interrogating the transponder too often. More details about Mode-S can be found in [ 1 ]. This mode is widely adopted in Europe and is being phased in for North America.

Mode-S signaling scheme uses squitter messages, which are defined as a non-solicited messages used in aviation radio systems. Mode-S has the following properties:

  • Transmit Frequency: 1090 MHz
  • Modulation: Pulse Position Modulation
  • Data Rate: 1 Mbit/s
  • Short Squitter Length: 56 microseconds
  • Extended Squitter Length: 112 microseconds

Short squitter messages contain the following information:

  • Downlink Format (DF)
  • Capability (CA)
  • Aircraft ID (Unique 24-bit sequence)
  • CRC Checksum

Extended squitter (ADS-B) messages contain all the information in a short squitter and one of these:

  • Altitude
  • Position
  • Heading
  • Horizontal and Vertical Velocity

The signal format of Mode-S has a sync pulse that is 8 microseconds long followed by either 56 or 112 microseconds of data as illustrated in the following figure.

Run the Example

Type plutoradioADSBExample in the MATLAB Command Window or click the link to run the example. You need to enter the following information when you run the example:

  1. Reception duration in seconds,
  2. Signal source (captured data or PlutoSDR),
  3. Optional output methods (map and/or text file).

The example shows the information on the detected airplanes in a tabular form as shown in the following figure.

You can also observe the airplanes on a map, if you have a valid license for the Mapping Toolbox.

Receiver Structure

The following block diagram summarizes the receiver code structure. The processing has four main parts: Signal Source, Physical Layer, Message Parser, and Data Viewer.

Signal Source

This example can use two signal sources:

  1. ''Captured Signal'': Over-the-air signals written to a file and sourced from a Baseband File Reader object at 2.4 Msps
  2. ''PlutoSDR'': ADALM-PLUTO radio at 12 Msps

If you choose ''PlutoSDR'' as the signal source, this application will search your computer for ADALM-PLUTO radios and ask you to choose one of them as the signal source.

The signal source is configured to provide enough samples to contain a maximum of 180 complete extended squitter messages. Since an extended squitter message is 288 samples long, set SamplesPerFrame property of the signal source object to 51840 samples, which is a 180 times 288. The rest of the algorithm searches for Mode-S packets in this frame of data and outputs all correctly identified packets. This type of processing is defined as batch processing. An alternative approach is to process one extended squitter message at a time. This single packet processing approach incurs 88 times more overhead than the batch processing, while it has 88 times less delay. Since the ADS-B receiver is delay tolerant, batch processing was used.

Physical Layer

The baseband samples received from the signal source are processed by the physical (PHY) layer to produce packets that contain the PHY layer header information and raw message bits. The following diagram shows the physical layer structure.

The baseband signal has a sampling rate of 12 MHz. Given the data rate is 1 Mbit/s, this corresponds to 12 samples per symbol. The rest of the receiver works on the magnitude of the interpolated complex values.

The packet synchronizer works on subframes of data equivalent to two extended squitter packets, that is, 1440 samples at 12 MHz or 120 micro seconds. This subframe length ensures that a whole extended squitter packet is contained in the subframe. The Packet synchronizer first correlates the received signal with the 8 microsecond preamble and finds the peak value. Then, it validates the synchronization point by checking if it matches the preamble sequence, [1 0 0 0 0 0 1 0 1 0 0 0 0 0 0], where a '1' represents a high value and a '0' represents a low value.

The Mode-S PPM modulation scheme defines two symbols. Each symbol has two chips, where one has a high value and the other has a low value. If the first chip is high followed by low chip, this corresponds to the symbol being a 1. Alternatively, if the first chip is low followed by high chip, then the symbol is 0. The bit parser demodulates the received chips and creates a binary message. The binary message is validated using a CRC checker. The output of bit parser is a vector of Mode-S physical layer header packets that contains the following fields:

  • RawBits: Raw message bits
  • CRCError: FALSE if CRC checks, TRUE if CRC fails
  • Time: Time of reception in seconds from start of receiver
  • DF: Downlink format (packet type)
  • CA: Capability

Message Parser

The message parser extracts data from the raw bits based on the packet type as described in [ 2 ]. This example can parse short squitter packets and extended squitter packets that contain airborne velocity, identification, and airborne position data.

Data Viewer

The data viewer shows the received messages on a graphical user interface (GUI). For each packet type, the number of detected packets, the number of correctly decoded packets, and the packet error rate (PER) is shown. As data is captured, the application lists information decoded from these messages in a tabular form.

You can also launch the map and start text file logging using the GUI.

  • Start Data Logging - Save the captured data in a TXT file. You can use the saved data for post processing.
  • Launch Map - Launch the map where the tracked flights can be viewed.

Example Code

The receiver asks for user input and initializes variables. Then it calls the signal source, physical layer, message parser, and data viewer in a loop. The loop also keeps track of the radio time using the frame duration and lost samples reported by the signal source.

The latency output of the signal source is an indication of when the samples were actually received and can be used to determine how close to real time the receiver is running. A latency value of greater than one indicates that the receiver was not able to process the samples in real time. Latency is reported in terms of the number of frames. It can be between 1 and 128. If latency is greater than 128, then samples are lost.

% Request user input from the command-line for application parameters
userInput = plutoradioAdsbUserInput;

% Calculate ADS-B system parameters based on the user input
[adsbParam,sigSrc] = plutoradioAdsbConfig(userInput);

% Create the data viewer object and configure based on user input
viewer = plutoradioAdsbViewer('LogFileName', userInput.LogFilename);
if userInput.LogData
  startDataLog(viewer);
end
if userInput.LaunchMap
  startMapUpdate(viewer);
end

% Create message parser object
msgParser = helperAdsbRxMsgParser(adsbParam);

% Start the viewer and initialize radio time
start(viewer)
radioTime = 0;

% Main loop
while radioTime < userInput.Duration
  % Receive baseband samples (Signal Source)
  rcv = sigSrc();

  % Process physical layer information (Physical Layer)
  [pkt,pktCnt] = plutoradioAdsbRxPhy(rcv, radioTime, adsbParam);

  % Parse message bits (Message Parser)
  [msg,msgCnt] = msgParser(pkt,pktCnt);

  % View results packet contents (Data Viewer)
  update(viewer, msg, msgCnt);

  % Update radio time
  radioTime = radioTime + adsbParam.FrameDuration;
end

% Stop the viewer and release the signal source
stop(viewer)
release(sigSrc)

Further Exploration

You can explore following functions and System objects for details of the physical layer implementation:

Selected Bibliography

  1. International Civil Aviation Organization, Annex 10, Volume 4. Surveillance and Collision Avoidance Systems.
  2. Technical Provisions For Mode S Services and Extended Squitter (Doc 9871)