Figure 1.
Framework of the proposed CTG. Firstly, the complete trajectory is divided into backward, connecting, and forward parts; secondly, a new connecting part is generated to connect the backward and forward parts from different trajectories to obtain a new synthetic trajectory.
Figure 1.
Framework of the proposed CTG. Firstly, the complete trajectory is divided into backward, connecting, and forward parts; secondly, a new connecting part is generated to connect the backward and forward parts from different trajectories to obtain a new synthetic trajectory.
Figure 2.
Schematic diagram of training data preparation for connecting part generation model. Training data are collected by utilizing three partially overlapping sliding windows: the backward and forward windows, which are used to obtain input data, and the connecting window, which is used to obtain output data. The red, yellow, and blue boxes refer to the part of the backward, connecting, and forward windows that do not overlap with other windows, respectively. The orange and green boxes refer to the overlapping part of two adjacent windows.
Figure 2.
Schematic diagram of training data preparation for connecting part generation model. Training data are collected by utilizing three partially overlapping sliding windows: the backward and forward windows, which are used to obtain input data, and the connecting window, which is used to obtain output data. The red, yellow, and blue boxes refer to the part of the backward, connecting, and forward windows that do not overlap with other windows, respectively. The orange and green boxes refer to the overlapping part of two adjacent windows.
Figure 3.
Connection-based trajectory generation (CTG) in an optimization loop. and are two decision variables that describe the switch point and the alternative trajectory for the original trajectory, respectively. The CTG takes , , and the original trajectory as inputs and outputs a modified trajectory to update the optimization result.
Figure 3.
Connection-based trajectory generation (CTG) in an optimization loop. and are two decision variables that describe the switch point and the alternative trajectory for the original trajectory, respectively. The CTG takes , , and the original trajectory as inputs and outputs a modified trajectory to update the optimization result.
Figure 4.
Schematic diagram of majorization–minimization-based adversarial training. In the majorization step, a surrogate network is trained to approximate the flight performance constraints; in the minimization step, a generator is trained to minimize the output of the surrogate to consider the flight performance constraints.
Figure 4.
Schematic diagram of majorization–minimization-based adversarial training. In the majorization step, a surrogate network is trained to approximate the flight performance constraints; in the minimization step, a generator is trained to minimize the output of the surrogate to consider the flight performance constraints.
Figure 5.
Schematic diagram of the majorization step. The connecting parts generated by the generator are labeled according to the flight performance constraints and then added to the training dataset together with connecting parts that cannot meet the constraints in previous iterations. The labeled training data are used to train the supervised learning-based surrogate network, and the loss function is cross-entropy. The red, yellow, and blue boxes refer to the part of the backward, connecting, and forward trajectories that do not overlap with others, respectively. The orange and green boxes refer to the overlapping part of two adjacent trajectories.
Figure 5.
Schematic diagram of the majorization step. The connecting parts generated by the generator are labeled according to the flight performance constraints and then added to the training dataset together with connecting parts that cannot meet the constraints in previous iterations. The labeled training data are used to train the supervised learning-based surrogate network, and the loss function is cross-entropy. The red, yellow, and blue boxes refer to the part of the backward, connecting, and forward trajectories that do not overlap with others, respectively. The orange and green boxes refer to the overlapping part of two adjacent trajectories.
Figure 6.
A schematic diagram of the minimization step. The generator’s loss function includes reconstruction loss and a flight performance penalty. The reconstruction loss is the weighted mean square error between the real and generated connecting parts, and the flight performance penalty is the mean value of the surrogate network’s output.
Figure 6.
A schematic diagram of the minimization step. The generator’s loss function includes reconstruction loss and a flight performance penalty. The reconstruction loss is the weighted mean square error between the real and generated connecting parts, and the flight performance penalty is the mean value of the surrogate network’s output.
Figure 7.
Weight assignment for different trajectory points.
Figure 7.
Weight assignment for different trajectory points.
Figure 8.
Historical trajectories and six entry points.
Figure 8.
Historical trajectories and six entry points.
Figure 9.
Flight performance of different aircraft categories. Each row in the figure represents a wake turbulence category, while each column in the figure represents one or two flight performance indicators. A trajectory meets the flight performance constraints if all its points lie in the feasible area.
Figure 9.
Flight performance of different aircraft categories. Each row in the figure represents a wake turbulence category, while each column in the figure represents one or two flight performance indicators. A trajectory meets the flight performance constraints if all its points lie in the feasible area.
Figure 10.
The improvement made by VampPrior TCVAE. The commonly used fully connected layers in the encoder and decoder are replaced with temporal convolutional networks to better extract the spatial–temporal features from the trajectories, and the widely used standard Gaussian is replaced with the variational mixture of posteriors.
Figure 10.
The improvement made by VampPrior TCVAE. The commonly used fully connected layers in the encoder and decoder are replaced with temporal convolutional networks to better extract the spatial–temporal features from the trajectories, and the widely used standard Gaussian is replaced with the variational mixture of posteriors.
Figure 11.
Flyability changes throughout the training epochs. The MMAT technique is introduced to the experiment group after 600,000 epochs, and the flyable rate of the experiment group first sharply declines and then rapidly rebounds, and finally exceeds the control group.
Figure 11.
Flyability changes throughout the training epochs. The MMAT technique is introduced to the experiment group after 600,000 epochs, and the flyable rate of the experiment group first sharply declines and then rapidly rebounds, and finally exceeds the control group.
Figure 12.
Trajectory generation with the well-trained generator. (a) We randomly select some forward trajectories and connect them with the backward trajectory with connecting parts. (b) We roughly eliminate unflyable trajectories with the surrogate network. (c) We precisely eliminate unflyable trajectories according to the flight performance constraints.
Figure 12.
Trajectory generation with the well-trained generator. (a) We randomly select some forward trajectories and connect them with the backward trajectory with connecting parts. (b) We roughly eliminate unflyable trajectories with the surrogate network. (c) We precisely eliminate unflyable trajectories according to the flight performance constraints.
Figure 13.
Flight performance of trajectories generated by different models. Each row in the figure represents a trajectory generation model, while each column in the figure represents one or two flight performance indicators. A trajectory meets the flight performance constraints if all its points lie in the feasible area.
Figure 13.
Flight performance of trajectories generated by different models. Each row in the figure represents a trajectory generation model, while each column in the figure represents one or two flight performance indicators. A trajectory meets the flight performance constraints if all its points lie in the feasible area.
Figure 14.
Comparison with other trajectories. (a) Real northbound trajectories often pass through the published waypoints and only deviate from the published routes in specific regions. For example, an aircraft entering the TMA from ATAGA or IGONO tends to make a dog-leg in the yellow area. (b) Trajectories generated by VamPrior TCVAE ignore the published waypoints and could deviate from the published routes in any region. (c) Trajectories generated by CTG-MMAT can pass through the published waypoints and deviate at proper regions just like the real trajectories.
Figure 14.
Comparison with other trajectories. (a) Real northbound trajectories often pass through the published waypoints and only deviate from the published routes in specific regions. For example, an aircraft entering the TMA from ATAGA or IGONO tends to make a dog-leg in the yellow area. (b) Trajectories generated by VamPrior TCVAE ignore the published waypoints and could deviate from the published routes in any region. (c) Trajectories generated by CTG-MMAT can pass through the published waypoints and deviate at proper regions just like the real trajectories.
Figure 15.
Generating trajectories with different maneuvers. (a) Connecting to trajectory , which joins the base leg earlier than trajectory , to form a parallel-offset maneuver. (b) Connecting to trajectory , which has hold-on patterns, to form a hold-on maneuver. (c) Connecting to trajectory , which joins the final leg earlier than trajectory , to form a short-cut maneuver. (d) Connecting to trajectory , which has a dog-leg pattern, to form a dog-leg maneuver.
Figure 15.
Generating trajectories with different maneuvers. (a) Connecting to trajectory , which joins the base leg earlier than trajectory , to form a parallel-offset maneuver. (b) Connecting to trajectory , which has hold-on patterns, to form a hold-on maneuver. (c) Connecting to trajectory , which joins the final leg earlier than trajectory , to form a short-cut maneuver. (d) Connecting to trajectory , which has a dog-leg pattern, to form a dog-leg maneuver.
Figure 16.
Rerouting path planning for an aircraft entering from IGONO. At different positions, the ownship can switch to different historical trajectories.
Figure 16.
Rerouting path planning for an aircraft entering from IGONO. At different positions, the ownship can switch to different historical trajectories.
Figure 17.
Information of the models. Although the models for comparison have close trainable variable sizes, they have different FLOPs. The higher the FLOPs, the more complex the model.
Figure 17.
Information of the models. Although the models for comparison have close trainable variable sizes, they have different FLOPs. The higher the FLOPs, the more complex the model.
Figure 18.
Time consumption for each epoch.
Figure 18.
Time consumption for each epoch.
Figure 19.
The variation in reconstruction error over several epochs.
Figure 19.
The variation in reconstruction error over several epochs.
Figure 20.
The variation in flyable trajectory numbers over several epochs.
Figure 20.
The variation in flyable trajectory numbers over several epochs.
Table 1.
Adapted architecture of generator.
Table 1.
Adapted architecture of generator.
Layer Type | Layer Parameters | Output Size |
---|
Input layer | - | |
Flatten layer | - | 144 |
Dense layer | units = 100, activation = sigmoid | 100 |
Dense layer | units = 100, activation = sigmoid | 100 |
Dense layer | units = 36, activation = sigmoid | 36 |
Reshape layer | target shape = (12, 3) | |
Table 2.
Adapted architecture of discriminator.
Table 2.
Adapted architecture of discriminator.
Layer Type | Layer Parameters | Output Size |
---|
Input layer | - | |
Feature construction layer 1 | - | |
Flatten layer | - | 99 |
Dense layer | units = 330, activation = sigmoid | 330 |
Dense layer | units = 1, activation = sigmoid | 1 |