Next Article in Journal
Aerodynamic Shape Optimization of Subsonic/Supersonic Flows Integrating Variable-Fidelity Longitudinal Trim Analysis
Previous Article in Journal
Conceptual Design of Layered Distributed Propulsion System to Improve Power-Saving Benefit of Boundary-Layer Ingestion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms

by
Raymundo Peña-García
1,†,
Rodolfo Daniel Velázquez-Sánchez
1,†,
Cristian Gómez-Daza-Argumedo
2,†,
Jonathan Omega Escobedo-Alva
2,†,
Ricardo Tapia-Herrera
3,† and
Jesús Alberto Meda-Campaña
1,*,†
1
Instituto Politécnico Nacional SEPI ESIME Zacatenco, Ciudad de México 07738, Mexico
2
Instituto Politécnico Nacional SEPI ESIME Ticomán, Ciudad de México 07340, Mexico
3
CONAHCYT—Instituto Politécnico Nacional SEPI ESIME Zacatenco, Ciudad de México 07738, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Aerospace 2024, 11(2), 142; https://doi.org/10.3390/aerospace11020142
Submission received: 14 December 2023 / Revised: 25 January 2024 / Accepted: 31 January 2024 / Published: 8 February 2024

Abstract

:
This research introduces a physics-based identification technique utilizing genetic algorithms. The primary objective is to derive a parametric matrix, denoted as A, describing the time-invariant linear model governing the longitudinal dynamics of an aircraft. This is achieved by proposing a fitness function based on the properties of the transition matrix and taking advantage of some of the capabilities of the genetic algorithm, mainly those of restricting the search ranges of the unknowns. In this case, such unknowns are related to the type of aircraft and flight conditions that are considered during the identification process. The proposed identification method is validated with a reliable nonlinear model that can be found in the literature, as well as with the calculation of the trim condition and linearization generally used in aircraft dynamics. In summary, this study suggests that the genetic algorithm provided with the adequate fitness function could be an appealing alternative for aircraft model identification, even when limited data are available. Furthermore, in some cases, linearization using a genetic algorithm can be more efficient than classical methods.

1. Introduction

1.1. Motivation

One of the fundamental pursuits in the scientific domain is the development of mathematical models that accurately represent physical systems. A fundamental criterion for evaluating such models is their usefulness. This may involve the ability to predict crucial aspects of the behavior of a physical system or to provide critical information to the design process through model parameters. In all cases, the mathematical model must achieve a balance between simplicity to be of practical use, and complexity to describe the most important characteristics of the system being modeled [1,2].
A long-standing need in aeronautics has been that of obtaining reliable and accurate models for aircraft dynamics. This goes beyond aircraft design and flight simulation; it involves the implementation and validation of control systems. In the current era, because of the impressive development of unmanned aerial vehicles (UAVs), the search for control systems has become more demanding, both in civil and military applications. The arrival of high-reliability, low-cost sensors and sophisticated control systems has broken down previous barriers. Small aircraft can now reach higher altitudes, higher speeds, and longer ranges. This provides the opportunity to implement cutting-edge algorithms, which are both smarter and more efficient [3].
However, the effort to identify an optimal model for advanced control system design objectives remains a costly and complex task. Given this reality, it is imperative to reevaluate the identification process, seeking an innovative approach that is simple, fast, and effective. It is crucial to avoid heuristic models and costly experiments, recognizing that the simplest option is often the most prudent. Several compelling and direct examples can be found in [4]. Over the past two decades, identification in aeronautics, particularly for manned airplanes, has been vigorously explored by prominent civil military and government entities such as the Deutsches Zentrum für Luft-und Raumfahrt (DLR) and NASA [5].
A correct identification system incorporates experimental data and undergoes validation in flight. Nowadays, a huge amount of information on experimental identification and validation processes can be obtained from various technical reports and manuals [6], affording insights into many facets of authentic identification procedures. Recent research [5] offers an excellent overview of the state of the art in this field. Tools for the time domain are comprehensively examined in [7], with additional contributions from the German Aerospace Center featured in [8]. Furthermore, system identification techniques for aircraft dynamics in the frequency domain have been provided in [9].
In [5], it is noted that advances in machine learning have yet to yield substantial benefits for aircraft system identification, primarily due to the absence of a physics-based model. This limitation becomes evident when employing neural networks to represent dynamics [10] or aircraft dynamics as a black box [11], or when resorting to genetic algorithms in pursuit of a heuristic model [12]. There are also published works where two techniques are combined; for instance, in [13], genetic algorithms are part of a hybrid identification system originally based on optimization, and, in [14], the blend of neural networks and fuzzy systems is considered to obtain the longitudinal model of an aircraft.
When the integration of artificial intelligence is judiciously combined with an understanding of physics-based aircraft models, along with the aforementioned state-of-the-art knowledge, it enables the imposition of parametric constraints that facilitate the use of swift convergence algorithms. This approach circumvents the conventional issue of infinite solutions in computational algorithms and yields more cost-effective results compared to relying solely on artificial intelligence or physics-based algorithms [15]. As demonstrated, this robust identification methodology proves invaluable not only for linear dynamics identification but also for the linearization under diverse trim conditions. It offers a simple, effective, and tailored solution to the formidable challenges encountered in today’s aerospace industry.
It is generally recognized that the identification of nonlinear models poses a greater challenge compared to linear models. For example, successful identification through neural networks requires an appropriate determination of the number of layers, neurons per layer, and learning factor, all depending on the characteristics of the nonlinear system under consideration. Similarly, for fuzzy systems, an optimal approximation depends on the selection of an appropriate number of fuzzy rules, together with a judicious choice of membership functions. Even in the case of SINDY, customization is often imperative to address the unique complexities associated with the nonlinear system being identified [16]. On the other hand, when the mentioned techniques are applied to identify linear systems, the problem reduces to parametric identification [17].
While a large number of approaches are available in the literature to successfully identify both linear and nonlinear systems, this paper focuses specifically on the utilization of a practical and effective technique based on genetic algorithms to linearize the longitudinal dynamics of an aircraft by introducing a new fitness function.

1.2. Contribution

A notable disadvantage of the application of genetic algorithms in obtaining mathematical models from real or simulated data is that, during the evaluation of the fitting function, the model obtained must be simulated in each iteration to determine the error between the data and the response of the intermediate model. This causes a very high computational cost that may be out of the reach of researchers with limited budgets. To reduce this computational cost, and considering that the linear system is sufficient to describe the longitudinal dynamics of an aircraft, in this work, the transition matrix [18,19] is considered instead of the complete simulation of the intermediate model in each iteration of the genetic algorithm.
Therefore, the main goal of the current work is that of presenting a fitness function based on the transition matrix, such that, when the fitness function is included in the genetic algorithm, the matrix, A, involved in the linear approximation of the aircraft, is promptly obtained, and relatively few data are needed.

1.3. Manuscript Organization

The rest of this paper is organized as follows. In Section 2, the nonlinear model of the aircraft is presented. Section 3 describes the trim condition considered. In Section 4, the traditional linearization of the nonlinear model is performed. Below, Section 5 presents a concise analysis of the parameters included in matrix A. The main result, i.e., the linearization achieved by the genetic algorithm, is discussed in detail in Section 6. Finally, in Section 7, some conclusions are presented.

2. Nonlinear 6-DOF Dynamic Model of a Research Aircraft

The primary function of a nonlinear model is to represent and capture complex relationships between variables that do not follow a linear pattern, allowing greater flexibility in modeling intricate patterns and dependencies. It is important to mention that, although in the current work only the longitudinal dynamics are excited, it is necessary to pay attention to the complete dynamics of the aircraft. Next, the nonlinear model used in the analysis and validation of the identification process of the genetic algorithm is presented synthetically. This complete and validated model was extracted from [20], which is an excellent summary of [21].
In accordance with rigid body dynamics and classical approximation from Newton–Euler, the following analysis can be carried out:
Translation equation on body axes:
V ˙ c m / e b = 1 m F b ω b / e b × V b ,
where m is the mass, the sub-index c m / e means from the center of mass to Earth axis, and super-index b indicates body frame; this nomenclature comes from kinematics equations deduction (see [22]), with  V ˙ c m / e b = [ u ˙ v ˙ w ˙ ] T as the linear acceleration vector, V c m / e b = [ u v w ] T as the linear velocity vector, F b as the external forces, ω ˙ b / e b = [ p ˙ q ˙ r ˙ ] T as the angular acceleration, and  ω b / e b = [ p q r ] T as the angular rate.
Rotational equations on body axes:
ω ˙ b / e b = ( J b ) 1 M b ω b / e b × J b ω b / e b
where J b is the inertial matrix, and  M b are the external moments.
Kinematic Euler Equation:
At this point, it is important to consider a transformation of angular rates from body axes to the general frame Euler system (inertial), as follows:
Φ ˙ = H ( Φ ) ω ˙ b / e b ,
where Φ ˙ = [ ϕ ˙ θ ˙ ψ ˙ ] T is the vector of Euler angles rate and
H ( Φ ) = 1 s e n ( ϕ ) t a n ( θ ) c o s ( ϕ ) t a n ( θ ) 0 c o s ( ϕ ) s e n ( ϕ ) 0 s e c ( θ ) s e n ( ϕ ) c o s ( ϕ ) s e c ( θ ) .
The full nonlinear model is considered nonautonomous and depends on the full state (even forces and moments are direct functions of the states) and control inputs, resulting in the following:
X ˙ = f ( X , U ) ,
with the state vector
X = x 1 x 12 T = u v w p q r ϕ θ ψ P N P E h ¯ T ,
where P N , P E , and P D are the x y z aircraft displacements on North–East–Down navigation frame system and h ¯ = P D .
As mentioned above, external forces and moments are functions of states and control inputs, resulting in the following:
F b ( X , U ) ,
M b ( X , U ) ,
which ultimately leads to
x ˙ 1 x ˙ 2 x ˙ 3 x ˙ 4 x ˙ 5 x ˙ 6 x ˙ 7 x ˙ 8 x ˙ 9 x ˙ 10 x ˙ 11 x ˙ 12 = 1 m F b ( X , U ) x 4 x 5 x 6 × x 1 x 2 x 3 ( J b ) 1 M b ( X , U ) x 4 x 5 x 6 × J b x 4 x 5 x 6 H ( Φ ) x 4 x 5 x 6 C S ( ϕ , θ , ψ ) u v w .
In Equation (5), C S ( ϕ , θ , ψ ) is the body to navigation axis transformation:
c θ c ψ ( c ϕ s ψ + s ϕ s θ c ψ ) ( s ϕ s ψ + c ϕ s θ c ψ ) c θ s ψ ( c ϕ c ψ + s ϕ s θ s ψ ) ( s ϕ c ψ + c ϕ s θ s ψ ) s θ s ϕ c θ c ϕ c θ ,
where c ( · ) and s ( · ) are c o s i n e and s i n e function, respectively.

The Aircraft

The original airplane is a twin-engine research civil aircraft with a weight mass of 120,000 kg and inertial moment I y of 7.68 × 10 6 kg · m 2 , a mean aerodynamic chord of 6.6 m, and a wing platform area of 260 m 2 . The dynamics of the control input and the allowable limits on the elevator are specified as 15 δ e ˙ 15 degrees per second and 25 δ e 10 degrees, respectively. For the two motors, the rate and end limits are 1.6 δ t h ˙ 1.6 degrees per second and 0.5 δ t h 10 degrees.

3. Trim Condition

In the context of identification, the equilibrium condition will be a steady-state flight with wings level and no sideslip [22]. Determining the equilibrium point in a coupled, nonlinear, multivariable system, such as Equation (9), is nontrivial and requires a specialized algorithm [23]. Before linearization, a fitting condition for the system must be selected in the general form X ˙ = f ( X , U ) :
f ( X e , U e ) = 0 .
In accordance with Equation (5), navigation variables X e R 9 and U e R 5 are neglected, i.e.,  u ˙ = v ˙ = w ˙ = p ˙ = q ˙ = r ˙ = 0 for steady-state flight. At this point, the methodology presented in [22] is considered, whose problem statement is the following:
min X Ω f ^ ( X ) = f ( X ) + C T H C ,
where H is a weight diagonal matrix, and  C R 9 + 1 is the restriction vector. The process is then performed with the Nelder and Mead Simplex minimization search algorithm, using the command fminsearch() MATLAB, which is a nonlinear programming solver to find the minimum of an unconstrained multivariable function using the derivative-free method. The cost function and initial lookup values are required. In this case, all functions in C that are equal to zero are
C = x ˙ 1 x ˙ 9 ( | V c m / e b | V e ) γ x 2 x 7 x 9 ;
in Equation (9), V e is the required operation speed in trim condition, | V c m / e b | = x 1 2 + x 2 2 + x 3 2 , and γ = x 8 a t a n ( x 3 / x 1 ) is the flight path angle. The weight matrix was defined as H = d i a g ( 3 ) after a trial and error process that provided valid results. The final function cost is f ^ ( X ) = C T H C .
Using fminsearch() with the cost function defined above, after 9 iterations and with a minimization error of 1 × 10 19 , the clipping condition X e obtained is x 1 = 110 m / s , x 3 = 2.0882 m / s , x 8 = 0.0634 rad , and all remaining states are zero.

4. Linearization

For the coefficient identification process, it is necessary to obtain the time-invariant linear system around the equilibrium point X e . To this end, one may consider Equation (4) with x ˙ R n , x R n , u R m , and the mapping function f : R n + m R n . This expression can be rewritten in an implicit form as follows:
F ( x ˙ , x , u ) = x ˙ f ( x , u ) = 0 ,
where f : R 2 n + m 0 n . By linearizing Equation (10) through a Taylor expansion series around the equilibrium point and considering only low-order terms (linear), the following results are obtained:
F ( x ˙ , x , u ) = F ( x ˙ 0 , x 0 , u 0 ) + F x ˙ ( x ˙ 0 , x 0 , u 0 ) ( x ˙ x ˙ 0 ) + F x ( x ˙ 0 , x 0 , u 0 ) ( x x 0 ) + F u ( x ˙ 0 , x 0 , u 0 ) ( u u 0 ) .
Now, the equilibrium condition,
F ( x ˙ 0 , x 0 , u 0 ) = 0 ,
must be evaluated. Then, from Equations (10) and (12), x ˙ = x ˙ 0 , x = x 0 , u = u 0 , and define the partial derivatives as
F x ˙ ( x ˙ 0 , x 0 , u 0 ) = E ,
F x ( x ˙ 0 , x 0 , u 0 ) = A ^ ,
F u ( x ˙ 0 , x 0 , u 0 ) = B ^ ,
and, by simplifying Equation (11), one has
0 = E δ x ˙ + A ^ δ x + B ^ δ u .
Thus, solving for δ x ˙ yields the following:
δ x ˙ = A δ x + B δ u ,
where
A = A ^ E ,
B = B ^ E .
It is important to mention that the LTI system x ˙ = A x + B u with matrices Equations (18) and (19) will provide values of perturbations ( δ x ) from an initial equilibrium condition X e , as in the case of the control inputs.
The partial derivatives appearing in Equations (13)–(15) can be obtained by using a numerical algorithm that employs a finite differential approximation, a feasible approach in simulation where one has control over the reliability and integrability of the data. While the process is relatively straightforward, careful consideration must be afforded to numerical errors, which can be mitigated by avoiding excessively small values of h. In the literature, it is widely acknowledged that central differentiation is more accurate compared to one-sided difference formulas [24,25]. Then, for any partial derivatives, one has
F ( t k ) x ( t k ) = F ( x ( t k + 1 ) ) F ( x ( t k 1 ) ) 2 h ,
where h = x ( t k + 1 ) x ( t k ) . This method is considered to determine E ^ , A ^ , and B ^ column by column, i.e.,
E i , j = F i ( x ˙ j ( t k + 1 ) ) F i ( x ˙ j ( t k 1 ) ) 2 h j , A ^ i , j = F i ( x j ( t k + 1 ) ) F i ( x j ( t k 1 ) ) 2 h j , B ^ i , j = F i ( u j ( t k + 1 ) ) F i ( u j ( t k 1 ) ) 2 h j .
To obtain data and calculate the matrices in Equation (21), the initial condition and the input signals are used as perturbations. For instance, to produce a response in all states of the nonlinear model with no input signal, x ( 0 ) 1 = 10 m / s and x ( 0 ) 2 = 5 m / s are chosen as initial conditions; these values are additional to the trim condition. Then, after an increment of the state variable x j from h j = 1 × 10 12 , the function F i is computed from each state x j to finally calculate the partial derivatives that make up matrices as in Equation (21). According to Equations (18) and (19), the state-space matrices for the 6-degree-of-freedom linear system are determined. Such a system is described by the differential equation x ˙ = A x + B u , where x R 9 and u R 5 . However, the longitudinal dynamics can be separated as in the following reduced system:
A = 0.052 0.0011 6.852 9.7903 0.2333 0.9104 107.9602 0.6215 0.0044 0.0431 1.4537 0 0 0 1 0 ,
B = 0 0.8022 0 9.8100 9.8100 0 12.6387 0 0 0 0 5.0424 0 0.3924 0.3924 0 0 0 0 0 .
The matrix A is directly related to Equation (24) and is considered parametric, implying a physical interpretation.

5. Parameters

The matrix coefficient in a state-space form has a well-known analytic expression in the field of flight dynamics. In this particular context, only longitudinal dynamics are considered. The matrix A is square and corresponds to the vector state x = [ u w q θ ] T . Each coefficient is assigned a specific nomenclature, delineating moments and forces as functions of the state variables, i.e.,
A = X u X w X q X θ Z u Z w Z q Z θ M u M w M q 0 0 0 1 0 .
In Equation (24), X θ and Z θ represent the linearized force projections due to gravity on body axes and are direct functions of the trim condition defined by θ e ; M refers to moments on body axes, specifically controlling pitch. Data for the common coefficient range parameters and their physical interpretation can be found in various sources, including [22,26,27,28]. These significant values are contingent upon the specific aircraft and prevailing flight conditions.

6. Main Result: Linearization by Means of the Genetic Algorithm (GA)

In this section, the nonlinear model introduced in Section 2 is employed as the benchmark for the evaluation of the performance of the linearization methods outlined in Section 4, including the genetic algorithm-based linearization approach. The maneuver considered corresponds to a cruise flight condition wherein the aircraft is traveling at a speed of 122.8357 m / s at sea level. The specified state variables include an angle of attack of 1.1354 degrees, a pitch angle of −3.7701 degrees, and all remaining states set to zero. The trim elevator input is −5.9505 degrees, and both motors operate at the 67.5% of total power. Utilizing the nonlinear model defined by Equations (1)–(6), the longitudinal behavior of the aircraft is simulated, deviating from the trim condition with perturbations of Δ u 0 n l = 10 m / s , Δ α 0 n l = 2.3309 degrees, Δ q 0 n l = 12 degrees per second, and Δ θ 0 n l = 0  degrees over a duration of t = 180 s with a sampling time of T = 0.05 s .
In order to establish a reference for comparing the genetic algorithm (GA) linearization, the linear system with matrix A, derived through the linearization procedure detailed in Section 4 and represented by Equation (22), is simulated under identical conditions. The numerical results are illustrated in Figure 1, considering only perturbation data from the trimmed flight condition above. The mean squared error (MSE) between the linear states and their corresponding nonlinear counterparts is computed as 0.0034 .
Next, the nonlinear system is linearized using a genetic algorithm. Before presenting the results, a concise description of the genetic algorithm is provided.

6.1. Genetic Algorithm (GA)

A genetic algorithm (GA) is an optimization technique that mimics the processes of evolution and natural selection. Usually, it is applied to solve complex optimization problems.
In a genetic algorithm, a potential solution to the problem is represented as a chromosome, typically as a binary string, but other representations are also possible, depending on the problem at hand. A population is a collection of such chromosomes. The size of the population is a user-defined parameter. The fitness function evaluates how well each chromosome in the population solves the problem. It assigns a fitness score to each chromosome based on the goal of the problem. The fitness function is problem-specific and provided by the user.
Chromosomes with higher fitness scores have a greater chance of being selected for reproduction. The most common selection methods include roulette wheel selection and tournament selection. Crossover is the process of combining two parent chromosomes to create one or more offspring. Various crossover operators are used, such as one-point, two-point, and uniform crossover. The crossover rate (probability of crossover) is a user-defined parameter. Mutation involves randomly changing some bits in a chromosome with a low probability. It introduces genetic diversity into the population. The mutation rate (probability of mutation) is another user-defined parameter. Genetic algorithms terminate when a stopping criterion is met. Common termination conditions include a maximum number of generations, reaching a satisfactory fitness level, or a time limit.
A typical genetic algorithm operates in a loop, which includes selection, crossover, mutation, and replacement of the old population with the new one. This process continues for a predefined number of generations [29,30,31,32,33]. It is important to mention that the GA has been successfully used to solve a number of problems in aeronautics. For instance, a methodology involving a multiobjective genetic algorithm for exploring parametrized microstructures is presented in [34]. In [35], the authors investigate the aerodynamics of desert locust tandem wings through computational fluid dynamics (CFD) simulations. Using 2D and 3D Navier–Stokes equations, they explore wing interactions and corrugation effects. Optimization with genetic algorithms and Nash game theory improves gliding performance by at least 77%. In [36], an automated framework for developing interpolatable aeroelastic reduced-order models (AE ROMs) across diverse flight conditions is presented. By combining system identification, state-consistency enforcement (SCE), and a genetic algorithm (GA), the approach addresses the issue of state inconsistency in ROMs.
In MATLAB, optimization via genetic algorithm (GA) can be performed using the function ga, which is part of the Global Optimization Toolbox. The syntax is as follows:
x = g a ( f u n , n v a r s , A , b , A e q , b e q , l b , u b , n l c o n , o p t s ) ,
where f u n is the fitness function to be optimized; x is the vector of variables included in the fitness function; n v a r s is the dimension of the vector x; the pair ( A , b ) defines the inequality restrictions in the form of A · x b ; the pair ( A e q , b e q ) defines the equality restrictions in the form of A e q · x = b e q ; the vectors l b and u b , both of dimension n v a r s , define the lower and upper limits for the search of x, respectively; n l c o n defines the nonlinear constraints; o p t s allows the user to modify parameters such as tolerance, the type of plot to be depicted, the maximal number of generations, etc. For more details, please refer to [37].
An essential challenge when working with GA is the need to define an appropriate fitness function. This is a crucial factor as the fitness function plays a critical role in evaluating how effectively a solution addresses the specific problem. The selection of the fitness function can significantly influence the overall performance of the algorithm and the quality of the solutions it produces.
With this in mind, a simple yet practical fitness function is introduced in the following sections. The function is designed to derive a parametric matrix A that characterizes the longitudinal dynamics of an aircraft. It utilizes data obtained from the simulation of the highly reliable nonlinear aircraft model described in Section 2.

6.2. Linearization via GA

In this section, the objective is to find a matrix A in the form of Equation (24) for the linearization of the longitudinal dynamics of the aircraft described in Equations (1)–(6) using genetic algorithms (GA). To accomplish this, a fitness function is derived in a practical manner. The unknown x s values to be determined by the GA are organized in matrix A G A as follows:
A G A = x ( 1 ) x ( 4 ) x ( 7 ) 9.7903 x ( 2 ) x ( 5 ) 107.9602 x ( 9 ) x ( 3 ) x ( 6 ) x ( 8 ) 0 0 0 1 0 ,
where the constants appearing in Equation (26) are fixed values dependent on the characteristics of the modeled aircraft. In other words, only nine of the sixteen elements of the matrix regarding Equation (26) need to be determined.
It is proposed that the fitness function evaluates the matrix described in Equation (26) in each iteration, comparing the linear states with their corresponding nonlinear states at arbitrary sample instants. To minimize computational load, the inclusion of the transition matrix is considered in this process.

6.2.1. Transition Matrix

The transition matrix, denoted as Φ ( t t 0 ) = e A ( t t 0 ) , describes the evolution of a time-invariant linear system over time, considering a given initial condition x ( t 0 ) . It shows how the state of the system at an initial time t 0 is transformed to the state at a later instant t. For continuous-time systems described by x ˙ ( t ) = A x ( t ) , the transition matrix can be used to obtain x ( t ) from the initial condition x ( t 0 ) (with t 0 = 0 ) as follows:
x ( t ) = Φ ( t , 0 ) x ( 0 ) .
Please refer to [18,19] for further details. In MATLAB, the transition matrix can be computed using the function e x p m ( ) , and its application in this work is explained below.

6.2.2. Fitness Function

For optimized computational resource utilization, the proposal is to apply the fitness function to a reduced set of samples. Rather than utilizing the entire set of simulated data, consisting of 3601 data points collected over a duration of 180 s with a sampling interval of T = 0.05 s , a smaller subset must be considered. Such a subset is suggested to comprise the region where the linear approximation is desired to closely match the nonlinear data.
Taking this into account, from the simulation of the nonlinear model Equations (1)–(6) and based on the transition matrix mentioned above, the following algorithm describes the fitness function used to linearize the nonlinear dynamics of the aircraft.
Explanation:
  • Line 1: The name of the function is f i t n e s s , and the argument x R 9 is the vector of the variables to be found by the GA in order to minimize this function.
  • Line 2: Global variables n d a t a , x N L , x 0 , and t are used in this function, where n d a t a is the number of samples considered during the linearization, x N L R n d a t a × 4 is a matrix containing the values of the nonlinear states u , w , q , θ , from the numerical simulation at desired instants, x 0 represents the vector of initial conditions, and  t R n d a t a is a vector with the corresponding time instants at which matrix x N L has been obtained.
  • Line 3: Matrix A G A is constructed from elements of vector x from a previous iteration.
  • Line 4: Matrix x L R n d a t a × 4 is initialized as a matrix of zeros with n d a t a rows (number of samples) and four columns (states u , w , q , θ ).
  • Lines 5–7: A loop iterates over the elements of t, using the exponential matrix operation to compute the states of x L at each one of the time instants included in vector t.
  • Lines 8–11: Squared errors s q e i for each state, with  i = 1 , , 4 .
  • Line 12: The total error f is computed as the square root of the sum of squared errors.
  • Line 13: The function returns f as the result of the fitness evaluation.

6.3. Numerical Simulations

In this section, the results of two linearizations using the genetic algorithm (GA) are presented. These experiments were performed on a personal computer with the following specifications: Intel(R) Core(TM) i5-10500H CPU @ 2.50 GHz, 16 GB of RAM, running MATLAB 2019b.

6.3.1. Simulation 1: GA Linearization Using 31 Samples

In this case, only the samples obtained during the first 3 s of the nonlinear simulation with a sampling interval of T = 0.1 s are considered. This involves selecting every other sample from the nonlinear simulation, resulting in a total of 31 samples. In this example, according to the fitness function outlined in Algorithm 1, one has n d a t a = 31 , x N L R 31 × 4 , x 0 = [ 1050.20940 ] T , and  t R 31 .
Algorithm 1 Function f i t n e s s
  1:
procedure Fitness(x)
  2:
    global  n d a t a , x N L , x 0 , t
  3:
     A G A x ( 1 ) x ( 4 ) x ( 7 ) 9.7903 x ( 2 ) x ( 5 ) 107.9602 x ( 9 ) x ( 3 ) x ( 6 ) x ( 8 ) 0 0 0 1 0
  4:
     x L zeros ( size ( t ) , 4 )
  5:
    for  k = 1  to  n d a t a  do
  6:
         x L ( k , : ) expm ( A G A · t ( k ) ) · x 0
  7:
    end for
  8:
     s q e 1 ( x N L ( 1 ) x L ( 1 ) ) T · ( x N L ( 1 ) x L ( 1 ) )
  9:
     s q e 2 ( x N L ( 2 ) x L ( 2 ) ) T · ( x N L ( 2 ) x L ( 2 ) )
10:
     s q e 3 ( x N L ( 3 ) x L ( 3 ) ) T · ( x N L ( 3 ) x L ( 3 ) )
11:
     s q e 4 ( x N L ( 4 ) x L ( 4 ) ) T · ( x N L ( 4 ) x L ( 4 ) )
12:
     f s q e 1 + + s q e 4
13:
    return f
14:
end procedure
During the execution of the GA in MATLAB, neither inequality nor equality restrictions nor nonlinear constraints are considered. However, due to the limitations of the available computer resources, the decision was made to search for the solution around the values obtained by the linearization procedure explained in Section 4.
Therefore, the GA in MATLAB was executed using the following three lines of code:
o p t s = o p t i m o p t i o n s ( g a , P l o t F c n , @ g a p l o t b e s t f ) f u n = @ f i t n e s s x = g a ( f u n , 9 , [ ] , [ ] , [ ] , [ ] , l b , u b , [ ] , o p t s )
with
l b = [ 0.06 , 0.3 , 0.005 , 0.002 , 1 , 0.05 , 5.5 , 2 , 0 ] u b = [ 0 , 0 , 0 , 0 , 0 , 0 , 7 , 0 , 0.7 ] ,
It is worth mentioning that the options chosen in Equation (28) are primarily intended for monitoring the GA’s behavior during execution, which may lead to a higher consumption of computer resources. A more refined selection of these parameters could potentially yield more efficient results. For a comprehensive description of MATLAB’s function g a , interested readers are referred to [37].
The elapsed time is 355.863055 s, and the matrix obtained by the GA is
A G A 1 = 0.0572 0.0020 5.5 9.7903 0.2367 0.9198 107.9602 0.0001 0.005 0.05 1.5393 0 0 0 1 0 .
Then, matrix Equation (30) is used to simulate the linear system for 180 s with a sampling time of T = 0.05 s in order to compare the states of the linear identified system with the nonlinear ones. The results are presented in Figure 2. The overall mean square error (MSE) between the nonlinear states and their linear approximations is 0.0167.
Notice that the MSE obtained in this experiment is greater than the one obtained by the linearization of Section 4, but the GA has been executed considering only 31 samples. In the next experiment, a larger set of samples is considered.

6.3.2. Simulation 2: GA Linearization Using 66 Samples

Now, the samples that occurred during the initial 3 s of the nonlinear simulation with a sampling time of T = 0.1 s , along with samples taken from the remainder of the signal with T = 5 s , are considered, resulting in a total of 66 samples. Therefore, one has n d a t a = 66 , x N L R 66 × 4 , x 0 = [ 1050.20940 ] T , and t R 66 . The execution of the GA in MATLAB is exactly the same as presented in Equations (28) and (29). The elapsed time is 1304.814475 s, and the matrix obtained by the GA is
A G A 2 = 0.0549 0.0019 6.1047 9.7903 0.2209 0.9005 107.9602 0.0001 0.0046 0.0499 1.5690 0 0 0 1 0 .
As before, the linear system with state matrix Equation (31) is simulated for 180 s with a sampling time of T = 0.05 s , and the comparisons with the nonlinear states are depicted in Figure 3. The overall mean square error (MSE) between the nonlinear states and their linear approximations is notably smaller: 4.3984 × 10 4 , but the execution time is larger than in the previous case.

6.4. Linearization via Least Squares (LS) Method

In contrasting the potential of the genetic algorithm (GA), a presentation of a simple linear regression utilizing the least squares method is provided. It is widely recognized that perfect identification results when applying this method to data generated by a linear time-invariant (LTI) system without noise or external perturbations. However, complications in identification may arise when applying this method to data from a nonlinear system [7]. The method was, nonetheless, applied using the same nonlinear simulation data employed in the GA algorithm. Initially, the analysis was conducted using a limited dataset comprising 31 samples, with a final time of 3 s and samples taken at 0.1 s intervals. The matrix obtained is presented in Equation (32), and the results of the numerical simulation are depicted in Figure 4. The overall MSE between the nonlinear states and their linear approximations is 6.2725 × 10 6 .
A L Q 1 = 0.003 0.103 10.4008 14.6247 0.0805 1.0208 118.5216 1.5095 0.0032 0.0489 1.5903 0.2304 0 0 1 0 .
Finally, the entire dataset produced by the nonlinear simulation was considered, spanning 180 s with a sample interval of 0.05 s. In this case, the matrix obtained is shown in Equation (33), while the numerical results are illustrated in Figure 5. In this case, the overall MSE is 235.2556 .
A L Q 2 = 0.001 0.0095 2.3994 0.21 0.0056 0.037 3.3085 0.2203 0 0.004 0.0231 0.0026 0 0.0001 0.0018 0.0049 .
This analysis is summarized in Table 1.

7. Conclusions

The traditional linearization algorithm is based on the idealized assumption of perfectly clean and complete data, often considering h = 1 × 10 12 , rendering it impractical for real-world experimental data. Nonetheless, it has been demonstrated in this study that the genetic algorithm (GA) not only generates matrices with clear physical interpretations but also exhibits the essential robustness necessary for reliable experimental identification. Furthermore, the automation of the linearization process is achievable.
In situations characterized by limited computational resources, where applying GA from scratch to linearize nonlinear dynamics is impeded by unknown search ranges, the GA approach based on the proposed transition matrices proves valuable as a post-processing tool. This is particularly advantageous when measurements of the nonlinear system are available, as demonstrated in this study, showcasing its capacity to enhance results obtained through traditional linearization methods.
Most identification methodologies, including the one presented here, can demonstrate robustness against noise when appropriately designed. With this proposal, assurance is provided that matrix A possesses a physical interpretation, and the presence of nonlinearity does not impede the linear identification process.
In contrast to other methodologies that face challenges in deriving a matrix A with a clear physical interpretation, this approach excels in this regard. Furthermore, the model obtained through this methodology proves functional across a broader operational zone compared to alternative methods.
Considering these aspects, this approach emerges as an alternative for identifying linear time-invariant (LTI) equations from experimental data. The exploration of this potential, particularly in broader experimental scenarios, constitutes an essential avenue for future research.

Author Contributions

Conceptualization, J.O.E.-A. and J.A.M.-C.; formal analysis, R.P.-G., R.D.V.-S., C.G.-D.-A., J.O.E.-A. and R.T.-H.; investigation, J.O.E.-A., R.T.-H. and J.A.M.-C.; methodology, J.O.E.-A., R.T.-H. and J.A.M.-C.; project administration, J.A.M.-C.; writing—original draft, J.O.E.-A. and J.A.M.-C.; writing—review and editing, J.A.M.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

This work was supported in part by Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCYT), through the Scholarship Sistema Nacional de Investigadores (SNI); in part by Instituto Politécnico Nacional through Research Project under Grant 20230023 and 20231126; in part by the Scholarship Estímulo al Desempeño de los Investigadores (EDI); in part by the Scholarship Comisión de Operación y Fomento de Actividades Académicas (COFAA); and in part by the Scholarship Beca de Estímulo Institucional de Formación de Investigadores (BEIFI).

Conflicts of Interest

The authors declare no conflicts of interest.

Notation and Definitions

AMatrix for LTI longitudinal dynamics
A ^ Jacobian matrix of nonlinear equations with respect to the state variables
BInput Matrix for LTI system
B ^ Jacobian matrix of nonlinear equations with respect to inputs
c ( · ) Cosine function
CMatrix for output definition in state form
EJacobian matrix of nonlinear equations with respect to the first-order state variables
F b Force vector in body frame
F ^ Jacobian matrix of nonlinear equations with respect to the first-order state variables
hStep size for numerical derivative, h = x ( t k + 1 ) x ( t k )
h ¯ Aircraft flying height
HWeight diagonal matrix for optimization criteria
H ( Φ ) Transformation matrix of angular rates from body axes to the inertial system
J b Inertial tensor of body center of gravity
mAircraft mass
M b Moment vector in body axes
M q Moment due to pitch rate
M u Moment due to horizontal velocity in body axes
M w Moment due to vertical speed in body axis
pRoll rate in x body axis
p ˙ Angular acceleration in the x-axis in body frame
P E Aircraft displacement to the east direction (navigation frame)
P N Aircraft displacement to north direction (navigation frame)
qPitch rate in y body axis
q ˙ Angular acceleration in the y-axis in body frame
rYaw rate in body axes
r ˙ Angular acceleration on the z-axis in body frame
s ( · ) Sine function
TStep size for genetic algorithm
uFlight velocity in x-axes in body frame from the equilibrium state (cruise flight)
UInput vector control
u ˙ Flight acceleration in x-axes in the body frame and from the equilibrium state (cruise flight)
U e Input vector control for equilibrium condition
u 0 Input for equilibrium condition in the linear system
vLateral component for flight speed from equilibrium condition, in body frame
V c m / e b Aircraft velocity from the center of mass to Earth axis, super-index b indicates body frame
V e Flight velocity in trim condition
v ˙ Lateral component for acceleration from equilibrium point, in body frame
wVertical component for flight velocity from equilibrium condition, in body frame
w ˙ Vertical component for acceleration from equilibrium condition, in body frame
XTwelve dimensional state vector
xx displacement on body frame
x 0 n l Initial condition for non linear model
X e Vector state in trim condition
X θ Forces in x body axis due to pitch angle
x 0 State values in equilibrium point for the linear system definition
x 0 ˙ First-order state values in equilibrium point for the linear system definition
X ˙ First-order vector state
x N L Matrix containing the values of the nonlinear states
X q Forces in x body axis due to pitch rate
X u Forces in the x-body axis due to velocity in the x-axis of the body frame
X w Forces in the x-body axis due to velocity in the z-axis of the body frame
yLateral displacement in body axis
zVertical displacement in body axis
Z θ Forces in z-body axis due to pitch angle
Z q Forces in z-body axis due to pitch rate
Z u Forces in the z-body axis due to velocity in the x-axis of the body frame
Z w Forces in the z-body axis due to velocity in the z-axis of the body frame
γ Longitudinal flight path angle
θ Pitch angle
θ e Pitch angle for trim condition
θ ˙ Pitch rate
Φ ˙ Vector of Euler angle rates, angular rates of aircraft respect to the inertial system
ϕ Roll angle
ϕ ˙ Roll rate
Φ ( t t 0 ) Transition matrix Φ ( t t 0 ) = e A ( t t 0 )
Φ ( t , 0 ) Transition matrix with initial condition x ( 0 ) = 0
ψ Yaw angle
ψ ˙ Yaw rate
ω b / e b Angular rate vector in body axis
ω ˙ b / e b Angular acceleration vector in body axis

References

  1. Hanks, H.T.; Tran, H.T. Mathematical and Experimental Modeling of Physical and Biological Processes; Chapman and Hall/CRC: Boca Raton, FL, USA, 2009. [Google Scholar]
  2. Nelles, O. Nonlinear System Identification: From Classical Approaches to Neural Networks and Fuzzy Models; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 2001. [Google Scholar]
  3. Martinez, O.A.; Cardona, M. State of the Art and Future Trends on Unmanned Aerial Vehicle. In Proceedings of the 2018 International Conference on Research in Intelligent and Computing in Engineering (RICE), San Salvador, El Salvador, 22–24 August 2018; pp. 1–6. [Google Scholar] [CrossRef]
  4. Brunton, S.L.; Proctor, J.L.; Kutz, N.J. Discovering Governing Equations From Data by Sparse Identification of Nonlinear Dynamical Systems. Proc. Natl. Acad. Sci. USA 2016, 113, 3932–3937. [Google Scholar] [CrossRef] [PubMed]
  5. Grauer, J.A.; Morelli, E.A. Introduction to the Advances in Aircraft System Identification from Flight Test Data Virtual Collections. J. Aircr. 2023, 60, 1329. [Google Scholar] [CrossRef]
  6. Hosseini, B.; Steinert, A.; Hofmann, R.; Fang, X.; Steffensen, R.; Holzapfe, F. Advancements in the Theory and Practice of Flight Vehicle System Identification. J. Aircr. 2023, 6, 1419–1436. [Google Scholar] [CrossRef]
  7. Jategaonkar, R.V. Flight Vehicle System Identification: A Time-Domain Methodology; American Institute of Aeronautics and Astronautics, Inc.: Reston, VA, USA, 2015. [Google Scholar]
  8. Deiler, C.; Mönnich, W.; Seher-Wei, S.; Wartmann, J. Retrospective and Recent Examples of Aircraft and Rotorcraft System Identification at DLR. J. Aircr. 1996, 33, 1–26. [Google Scholar] [CrossRef]
  9. Tischler, M.B.; Remple, R. System Identification. Engineering Methods with Flight Test Examples; AIIA Education: Richmond, VA, USA, 2006. [Google Scholar]
  10. Cedillo, S.G.T.; Bonello, P. Empirical identification of the inverse model of a squeeze-film damper bearing using neural networks and its application to a nonlinear inverse problem. J. Vib. Control 2018, 24, 357–378. [Google Scholar] [CrossRef]
  11. Kirkpatrick, K.; May, J.; Valasekz, J. Aircraft System Identification Using Artficial Neural Networks. In Proceedings of the 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Grapevine, TX, USA, 7–10 January 2013; pp. 1–12. [Google Scholar]
  12. Ulinowicz, M.; Narkiewicz, J. Aircraft parameter Identification Using Genetic Algorithm. In Proceedings of the 29th Congress of the International Council of the Aeronautical Sciences, St. Petersburg, Russia, 7–12 September 2014; pp. 317–334. [Google Scholar]
  13. Bruce, P.; Kellett, M. Maximum likelihood identification of linear aircraft dynamics using a hybrid genetic algorithm. In Proceedings of the 36th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 12–15 January 1998; pp. 357–378. [Google Scholar]
  14. Roy, A.G.; Peyada, N.K. Longitudinal Aircraft Parameter Estimation Using Neuro-Fuzzy and Genetic Algorithm Based Method. In Proceedings of the AIAA Atmospheric Flight Mechanics Conference, Grapevine, TX, USA, 9–13 January 2017; pp. 1–12. [Google Scholar]
  15. Anderson, M.B. Genetic Algorithms in Aerospace Design: Substantial Progress, Tremendous Potential; Technical Report; Sverdrup Technology Inc.: Allen Park, MI, USA; TEAS Group: Manchester, UK, 2002. [Google Scholar]
  16. Brunton, S.L.; Kutz, J.N. Data Driven Science and Engineering Machine Learning, Dynamical Systems, and Control; Department of Applied Mathematics, University of Washington: Philadelphia, PA, USA, 2017. [Google Scholar]
  17. Miri, N.; Mohammadzaheri, M.; Chen, L. An enhanced physics-based model to estimate the displacement of piezoelectric actuators. J. Intell. Mater. Syst. Struct. 2015, 26, 1442–1451. [Google Scholar] [CrossRef]
  18. Kailath, T.; Sayed, A.H.; Hassibi, B. Linear Estimation; Prentice Hall: Englewood Cliffs, NJ, USA, 2000. [Google Scholar]
  19. Ogata, K. Modern Control Engineering; Pearson Education: Saddle River, NJ, USA, 2010. [Google Scholar]
  20. Lambrechts, P.; Bennani, S.; Looye, G.; Moormann, D. The RCAM design challenge problem description. In Proceedings of the Robust Flight Control, Toulouse, France, 7 March 1997; Magni, J.F., Bennani, S., Terlouw, J., Eds.; Springer: Berlin/Heidelberg, Germany, 1997; pp. 149–179. [Google Scholar]
  21. Lambrechts, P.F. Technical Publication TP-088-3, GARTEUR-FM(AG08); Technical Report; Group for Aeronautical Research and Technology in Europe (GARTEUR): Madrid, Spain, 1997. [Google Scholar]
  22. Stevens, B.L.; Lewis, F.L.; Johnson, E.N. Aircraft Control and Simulation: Dynamics, Controls Design, and Autonomous Systems, 3rd ed.; Wiley: Hoboken, NJ, USA, 2016. [Google Scholar]
  23. Hong-Gi, L. Linearization of Nonlinear Control Systems; Springer: Singapore, 2022. [Google Scholar]
  24. Kiusalaas, J. Numerical Methods in Engineering with MATLAB®, 3rd ed.; Cambridge University Press: Cambridge, UK, 2015. [Google Scholar] [CrossRef]
  25. Miller, G. Numerical Analysis for Engineers and Scientists; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar] [CrossRef]
  26. Roskam, J. Airplane Flight Dynamics and Automatic Flight Controls; DAR Corporation: Lawrence, KS, USA, 2003. [Google Scholar]
  27. Blakelock, J.H. Automatic Control of Aircraft and Missiles; Wiley: New York, NY, USA, 1991. [Google Scholar]
  28. McLean, D. Automatic Flight Control Systems; Prentice Hall: New York, NY, USA, 1990. [Google Scholar]
  29. Holland, J.H. Adaptation in Natural and Artificial Systems; University of Michigan Press: Cambridge, MA, USA, 1975. [Google Scholar]
  30. Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning; Addison-Wesley: Boston, MA, USA, 1989. [Google Scholar]
  31. Mitchell, M. An Introduction to Genetic Algorithms; MIT Press: Cambridge, MA, USA, 1998. [Google Scholar]
  32. Haupt, R.L.; Haupt, S.E. Practical Genetic Algorithms; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  33. López-González, A.; Meda-Campaña, J.A.; Hernández-Martínez, E.G.; Paniagua-Contro, P. Multi robot distance based formation using Parallel Genetic Algorithm. Appl. Soft Comput. 2020, 86, 105929. [Google Scholar] [CrossRef]
  34. Liu, S.; Acar, P. Parameter Space Exploration of Cellular Mechanical Metamaterials Using Genetic Algorithms. AIAA J. 2023, 61, 3633–3643. [Google Scholar] [CrossRef]
  35. Isakhani, H.; Yue, S.; Xiong, C.; Chen, W. Aerodynamic Analysis and Optimization of Gliding Locust Wing Using Nash Genetic Algorithm. AIAA J. 2021, 59, 4002–4013. [Google Scholar] [CrossRef]
  36. Shu, J.I.; Wang, Y.; Brown, A.; Kaminsky, A. Genetic-Algorithm-Guided Development of Parametric Aeroelastic Reduced-Order Models with State-Consistence Enforcement. AIAA J. 2023, 61, 3976–3994. [Google Scholar] [CrossRef]
  37. MathWorks. Genetic Algorithm Options (Global Optimization Toolbox); Mathworks: Natick, MA, USA, 2023. [Google Scholar]
Figure 1. Simulated data with an initial condition located away from the equilibrium point. The linearization process, as described in Section 4, is applied to obtain the corresponding linearized representation.
Figure 1. Simulated data with an initial condition located away from the equilibrium point. The linearization process, as described in Section 4, is applied to obtain the corresponding linearized representation.
Aerospace 11 00142 g001
Figure 2. Simulated data with initial condition located away from the equilibrium point. Linearization achieved via GA using the first three seconds of nonlinear simulated data, sampled at T = 0.1 s .
Figure 2. Simulated data with initial condition located away from the equilibrium point. Linearization achieved via GA using the first three seconds of nonlinear simulated data, sampled at T = 0.1 s .
Aerospace 11 00142 g002
Figure 3. Simulated data with initial condition located away from the equilibrium point. Linearization achieved using a GA using the initial three seconds of nonlinear simulated data, sampled at T = 0.1 s , followed by T = 5 s for the rest of the nonlinear data.
Figure 3. Simulated data with initial condition located away from the equilibrium point. Linearization achieved using a GA using the initial three seconds of nonlinear simulated data, sampled at T = 0.1 s , followed by T = 5 s for the rest of the nonlinear data.
Aerospace 11 00142 g003
Figure 4. Linearization achieved using linear regression with least squares method using the initial three seconds of nonlinear simulated data, sampled at T = 0.1 s .
Figure 4. Linearization achieved using linear regression with least squares method using the initial three seconds of nonlinear simulated data, sampled at T = 0.1 s .
Aerospace 11 00142 g004
Figure 5. Linearization accomplished through linear regression employing the least squares method, utilizing the complete dataset of the nonlinear simulated data sampled at T = 0.05 s .
Figure 5. Linearization accomplished through linear regression employing the least squares method, utilizing the complete dataset of the nonlinear simulated data sampled at T = 0.05 s .
Aerospace 11 00142 g005
Table 1. MSE and execution time.
Table 1. MSE and execution time.
Linearization viaMSEExecution Time
Section 4 0.0034 —–
GA (31 samples) 0.0167 355.86 s
GA (66 samples) 4.3984 × 10 4 1304.81 s
LS (31 samples) 6.2725 × 10 6 —–
LS (full dataset) 235.2556 —–
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Peña-García, R.; Velázquez-Sánchez, R.D.; Gómez-Daza-Argumedo, C.; Escobedo-Alva, J.O.; Tapia-Herrera, R.; Meda-Campaña, J.A. Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms. Aerospace 2024, 11, 142. https://doi.org/10.3390/aerospace11020142

AMA Style

Peña-García R, Velázquez-Sánchez RD, Gómez-Daza-Argumedo C, Escobedo-Alva JO, Tapia-Herrera R, Meda-Campaña JA. Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms. Aerospace. 2024; 11(2):142. https://doi.org/10.3390/aerospace11020142

Chicago/Turabian Style

Peña-García, Raymundo, Rodolfo Daniel Velázquez-Sánchez, Cristian Gómez-Daza-Argumedo, Jonathan Omega Escobedo-Alva, Ricardo Tapia-Herrera, and Jesús Alberto Meda-Campaña. 2024. "Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms" Aerospace 11, no. 2: 142. https://doi.org/10.3390/aerospace11020142

APA Style

Peña-García, R., Velázquez-Sánchez, R. D., Gómez-Daza-Argumedo, C., Escobedo-Alva, J. O., Tapia-Herrera, R., & Meda-Campaña, J. A. (2024). Physics-Based Aircraft Dynamics Identification Using Genetic Algorithms. Aerospace, 11(2), 142. https://doi.org/10.3390/aerospace11020142

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop