Performance Comparison of Meta-Heuristics Applied to Optimal Signal Design for Parameter Identification

This paper presents a comparative study that explores the performance of various meta-heuristics employed for Optimal Signal Design, specifically focusing on estimating parameters in nonlinear systems. The study introduces the Robust Sub-Optimal Excitation Signal Generation and Optimal Parameter Estimation (rSOESGOPE) methodology, which is originally derived from the well-known Particle Swarm Optimization (PSO) algorithm. Through a real-life case study involving an Autonomous Surface Vessel (ASV) equipped with three Degrees of Freedom (DoFs) and an aerial holonomic propulsion system, the effectiveness of different meta-heuristics is thoroughly evaluated. By conducting an in-depth analysis and comparison of the obtained results from the diverse meta-heuristics, this study offers valuable insights for selecting the most suitable optimization technique for parameter estimation in nonlinear systems. Researchers and experimental tests in the field can benefit from the comprehensive examination of these techniques, aiding them in making informed decisions about the optimal approach for optimizing parameter estimation in nonlinear systems.


Introduction
The design of input signals is undoubtedly one of the key elements for the success of system identification [1].Well-designed signals can stimulate and highlight the main dynamics of the system, allowing these behaviors to be accurately estimated and represented by parametric models [2].Consequently, they constitute a fundamental element for obtaining accurate models and, in turn, for tuning model-based controllers [3], designing state observers [4], and developing simulators [5], among other applications [6].
However, the task of Optimal Signal Design (OSD) is not a trivial problem.As presented in the literature, especially in [7][8][9], this task is defined as a mixed-integer nonlinear optimization problem.It requires competent approaches to simultaneously address signal design and parametric estimation problems.
Regardless of the methodology adopted, there is a consensus that the approach must generate signals that provide effective stimuli to the main dynamics of the system [10].These signals must exhibit appropriate excitation persistence, covering a wide range of amplitudes and frequencies.Furthermore, it is essential that the signals be feasible for the actuator system and adhere to safety conditions during the identification experiment development [11].
Meeting all these requirements is challenging, especially in nonlinear systems, where signals with high persistence and detailed encoding are necessary.In such cases, there can be a significant increase in computational complexity, leading to time-consuming computational solutions [12].
To address the Optimal Signal Design (OSD) problem, it is common to formulate the optimization problem using the inverse of the Fisher Information Matrix (FIM).This approach allows for measuring signal excitation persistence, irrespective of the type of estimation algorithm implemented.
Regarding the solution to the optimization problem, various deterministic and probabilistic optimization algorithms are employed [13].These include Steepest-descent, Conjugategradient, Quasi-linearization, and Newton-Raphson, among others.Additionally, bio-inspired meta-heuristics are classic examples of probabilistic methods [14].
A comprehensive survey of the literature reveals several algorithms in this context.Concerning solutions based on deterministic methods, works such as [15,16] serve as classic examples of applications in the field of OSD.In these approaches, signals are optimized using deterministic algorithms that target metrics based on the Fisher Information Matrix (FIM), including criteria like A-Optimality, E-Optimality, and D-Optimality Another approach found in the literature combines deterministic methods with optimal control theory.Examples of these combined approaches can be found in [17][18][19][20].
Specifically, in the work by [17], an identification methodology is applied to aerial vehicles.The optimal signals are generated through a combination of dynamic programming and optimal control, which enhances traditional identification maneuvers.It is worth noting the significant reliance on the optimal control input signal.While these approaches yield satisfactory results, it is evident that further improvements to the entry maneuvers could lead to a substantial increase in the problem's spatial complexity, rendering it intractable for the techniques presented.
In terms of methodologies, it is worth highlighting the works of [7,12,21], which focus on designing multilevel pseudo-random sequence signals.The work in [7] employed the PSO algorithm, ref. [12] utilized the Ant Colony Optimization (ACO) algorithm, and [21] opted for the Genetic Algorithm (GA) to address the OSD problem.What sets these works apart is their ability to design signals with advanced coding, reasonable processing times, and high-quality solutions.In the case of [7], the signal design also optimized the experiment time, a feature not found in other works.Additionally, the authors developed a dual-layer optimization methodology, using their own metrics that cater to specific needs.
Considering all the mentioned works, it becomes apparent that Optimal Signal Design can be highly complex, particularly when dealing with signals featuring multiple levels of amplitude or intricate coding.In such cases, the number of potential combinations in the solution space can become prohibitively large, potentially rendering the project infeasible or leading to extremely high computational processing times, as exemplified in [12].
In this context, the application of meta-heuristics has shown promise, as demonstrated in the works of [7,12,21].However, a more comprehensive investigation, involving a comparison of various techniques, is still needed.It is worth noting that this area of research is not yet extensive in the literature, but it holds promise for OSD and deserves further attention.
In light of this context, this paper presents a comparative study of meta-heuristics applied to the OSD problem for estimating parameters of Nonlinear Dynamic Systems (NDSs).Several meta-heuristics will be investigated using the methodology proposed by [22], referred to as rSOESGOPE.While the PSO algorithm [23] served as the optimizer, this paper will also explore other meta-heuristics to enhance the performance of rSOESGOPE.The traditional algorithms such as GA [24], Bat Algorithm (BA) [25], Whale Optimization Algorithm (WOA) [26], Grey Wolf Optimizer (GWO) [27], Salp Swarm Algorithm (SSA) [28], Arithmetic Optimization Algorithm (AOA) [29], and Multi-trial vector-based Differential Evolution (MTDE) [30] will be analyzed.
This comparative study aims to explore and assess the potential of these approaches for solving this class of optimization problems.To this end, real case studies involving an ASV with three DoFs and an aerial holonomic propulsion system are presented and discussed.
After this summary, the contributions of this paper can be highlighted as follows: • a performance comparison of eight meta-heuristics applied in the rSOESGOPE concept for parametric estimation, not yet explored in the literature; • a description of the rSOESGOPE concept that allows obtaining robust models from the use of multiple projection excitation signals; • a real case application involving a 3 DoF ASV, which presents a problem with a complex solution space for signal generation and a challenging parametric estimate due to the hydrodynamic phenomena involved; • an in-depth analysis of the ASV's identification process, providing valuable insights into the modeling and generation of identification signals.
This paper is organized as follows: Section 2 presents the OSD problem statement, demonstrating its characteristics and complexities; Section 3 introduces the rSOESGOPE method; Section 4 elaborates on the fundamental aspects for the application of the algorithms in the rSOESGOPE concept; Section 5 presents the results of applying the meta-heuristics to a real case study problem; Section 7 presents some concluding remarks.

Problem Formulation
Consider a nonlinear dynamic system R(Γ) which can be satisfactorily approximated by the nonlinear model M(Γ): where f and h are the system-dependent nonlinear functions on the state vector x ∈ R n , the output vector y ∈ R m , the input vector u ∈ R p , and the set of parameters of the model Γ ∈ R r .Additionally, it is admitted that the R(Γ) has some restrictions that need to be respected during its operation, represented by where x and x ∈ R n represent, respectively, the upper and lower vector states' operating limits, y and y ∈ R m express the lower and upper vessels' output vector operating limits, respectively, while u and u ∈ R p represent the upper and lower input signal operating limits, respectively.It is also assumed that Γ represents the unknown optimal set of the system's parameters, and M(Γ) best represents R(Γ) for any given signal u ∈ U, where U represents the entire possible input signal domain.
In this context, designing an optimal identification signal involves finding the best signal u ⊕ ∈ U that, when applied to the real system R(Γ), enables the subsequent parametric estimation Γ+ , satisfying the following relationship: M( Γ+ , u) ≈ R(Γ, u).Therefore, designing an optimal identification signal requires discovering u ⊕ through an optimization search defined by

S(M(
where Γ− is the initial parameter estimation Γ− and S is an arbitrary mixed-integer optimization algorithm that searches for the best u ⊕ ∈ U. Finally, the optimal parameter set is defined by where P is an arbitrary optimization procedure that uses the mathematical model M and the signal u ⊕ to estimate the parameter set Γ+ in which M( Γ+ , u) best represents R(Γ, u ⊕ ).Additionally, in this notation, it is possible to replace [x, y] by M(Γ, u) or R(Γ, u) without loss of generalization.

Robust SOESGOPE
The Robust SOESGOPE (rSOESGOPE) is a methodology developed for designing identification signals and promoting robust parametric estimation.
This approach is a derivation of the original SOESGOPE method developed by [7] for designing single optimal signals for optimal parametric estimation of NDS.
However, as demonstrated by [8], the SOESGOPE method does not yield satisfactory results in situations where the initial estimate has a high level of uncertainty regarding the real system.
To address these uncertainties, the rSOESGOPE method proposes the use of multiple identification signals Then, the same U ⊕ will also be able to excite R(Γ) and consequently obtain Γ+ that respects M( Γ+ , U ⊕ ) ≈ R(Γ, U ⊕ )." These concepts can be better understood in Figure 1, which illustrates the trust region surrounding Γ, the designed signals, and the stages of parametric estimation: Thinking in stages, the rSOESGOPE method proposes a two-step approach: • Generate the benchmark parameter sets P ⊕ = [ Γp 1 , Γp 2 , . . ., Γp n ] and use it as a parameter to find each signal of Apply U ⊕ to the real system to find Γ+ .
Thus, according to Section 2, the following steps can be defined: To tackle the two-step approach, a two-layer optimization strategy was developed to solve S. Two well-known optimization algorithms were employed: the PSO algorithm [31] and the Interior-Points Algorithm (IPA) [32].
The PSO algorithm was dedicated to the external optimization layer, specifically for the design of u(Ξ) and, consequently, the parameterization of Ξ, which represents the swarm's individuals.
In the proposed method, the projected signals are AutoRegressive exogenous inputs with an Amplitude-Modulated Pseudo-Random Binary Signal (APRBS).The internal layer is assigned to the IPA, which is responsible for the parametric estimation required for evaluating signal quality.This same IPA is used to solve P and find Γ+ .
To assess the signal quality generated by S, an objective function is proposed, consisting of a weighted sum of three metrics.
To ensure the analysis of excitation persistence, the following procedure was employed: 1.
PRECISION OUTPUT METRIC f o (•): This metric is responsible for measuring the difference between the states and outputs obtained from M( Γ+ , u) and M( Γp , u).It can be mathematically described as or, in the case of variables with significantly different dimensions, where ( x+ (k), ŷ+ (k)) and ( xp (k), ỹp (k)) are the states and outputs of M( Γ+ ) and M( Γp ), respectively.N represents the number of samples k in the experiment (t u ), and Q x and Q y are weightings of the state and output variables, respectively.These weightings allow for the treatment of variables based on their degrees of uncertainty.

2.
RECOVERABILITY METRIC f δ(•): This metric measures the capability to estimate the correct parameter set given u and Γp .Mathematically, f δ(•) is represented by the sum of the relative errors between the final estimation Γ+ and the reference Γp : where Γ+ i and Γp i represent the i-th parameter of Γ+ and Γp , respectively.This metric evaluates the similarity given the results provided by the previous metric f o (•).
It is important to mention that different parameter sets can generate the same output, and therefore, this metric also assesses how closely the estimated parameters match the benchmark set ( Γp ).
To ensure compliance with the constraints, Θ(•) was used to penalize a given signal u(Ξ) if it drives the system out of the desired operational restrictions.Mathematically, this can be expressed as follows: where x i and x i represent the lower and upper bounds of the i-th state, respectively, and y i and y i represent the lower and upper limits of the i-th output, respectively.Using the metrics described above, ref. [22] proposed the following objective function to be minimized.This objective function is composed of the weighted sum of the three metrics presented: where Ξ are the parameters that encode the signal, and k o , k δ, and k Θ ∈ R ≥0 are constant weightings related to metrics and established according to the priority.

Meta-Heuristic Application
This section presents fundamental aspects for the application of meta-heuristics in the rSOESGOPE concept.It is divided into two parts.Sections 4.1 and 4.2 present the signal representation and the adjustment criteria for meta-heuristics, respectively.

APRBS Representation and Settings
A fundamental step to applying any meta-heuristic in an OSD problem is the signal definition and its representation.These definitions determine the composition of the individuals in the technique and also the search space dimension.
In the rSOESGOPE method, the chosen signal is APRBS, which allows the incorporation of operational or safety restrictions together with the excitation persistence optimization [7,21,33].It is a fact that, in principle, unplanned signals cannot be guaranteed.
The APRBS strategy generates a coding vector Ξ k ∈ R 2×k , defined by k stages or pairs of amplitude and time interval, mathematically represented by [21] where [t 1 , t 2 , . . ., t k ] are the interval times and [A 1 , A 2 , . . ., A k ] are the respective amplitudes.Therefore, this APRBS signal with six stages yields an encoded signal Ξ 6 ∈ R 12 or Ξ 6 = [A 1 , . . ., A 6 , t 1 , . . ., t 6 ]. Figure 2 shows this example, where the realization of these parameters in u(Ξ 6 ) is marked as black: Thus, the proposed challenge for the meta-heuristics is the determination of the elements Ξ k , whose dimension depends on the number of stages k, where defining the number of stages is not a trivial task.As the object of study can vary, it is difficult to determine or establish a general rule.
Furthermore, it is known that the greater the number of signal stages, the greater the chance of producing signals with a wide range of amplitude and frequency.However, the larger the solution space is, the greater the effort required by the meta-heuristic, and therefore, the greater the number of individuals in the population should be.To address this situation, it is advisable to conduct systematic tests that analyze different configurations of the stages and population numbers, as well as the study of the necessary signals for U ⊕ .Therefore, by seeking a feasible trade-off relationship, this procedure ensures more accurate and appropriate definitions for the problem.

Meta-Heuristic Parameters
To ensure a fair comparison between the techniques, the algorithms were implemented with the same number of iterations and population size.
By default, the number of iterations is set to 100, while the population size was determined based on the study presented in Section 5.2.This study aimed to identify the optimal configuration, especially regarding the number of signal stages, which significantly impacts the solution space.
The remaining control parameters were defined in accordance with the original algorithm guidelines.You can find these parameters listed in Table 1.

Real Case Study: AERO4River
This section presents a study of metaheuristics applied in the rSOESGOPE concept.A real case study of the AERO4River ASV modeling is presented and discussed.
For this purpose, Section 5.1 introduces the AERO4River ASV and its mathematical model.Section 5.2 describes the procedures for applying the rSOESGOPE method and generating robust signals.Section 5.3 presents the results and implications of using multiple signals for each investigated optimization technique.Section 5.4 compares the results of each obtained model M( Γ+ ) with the real ASV, i.e., R(Γ) for different input signals.

AERO4River ASV
The AERO4River ASV is a vessel developed to perform autonomous hydrological measurements in environments with underwater obstacles and shallow and rapid water flows.It is based on a catamaran-type vessel and features an innovative air propulsion system with azimuth control, resulting in a highly maneuverable ASV with three overactuated DoFs, capable of meeting the required operating conditions.Figure 3 illustrates it: For more constructive information, it is recommended to consult the research in [22,34,35].

Kinematics and Dynamics Modeling
In the context of surface vessels, modeling requires the development of a description with three Degrees of Freedom (DoFs), namely, translation movements along the x and y axes ( surge and sway) and rotation around the z-axis ( yaw).Other Degrees of Freedom (DoFs) are disregarded due to their minimal influence on the vessel's dynamics [36].
According to marine vehicle nomenclature, it is customary to express state vectors as follows: η = [x, y, ψ] T , representing inertial positions (x, y) and angular position (ψ) in the vehicle's inertial frame; ν = [u, v, r] T , representing linear velocities (u, v) and angular velocity (r) in the rigid-body frame [37].
The general ASV dynamics and kinematics modeling without disturbance is represented as follows [38]: where M RB ∈ R 3×3 is the rigid-body Inertia Matrix; C RB (ν) ∈ R 3×3 is the Coriolis and centripetal rigid-body Matrix; J(ψ) is the Jacobian matrix that relates the velocities in the rigid-body and inertial frame; the matrices M A ∈ R 3×3 and C A (ν) ∈ R 3×3 represent the added mass phenomenon; D(ν) ∈ R 3×3 represents the hydrodynamic damping; τ ∈ R 3×1 is the vector of generalized forces and moments applied by the propellant system.The matrices M RB , C RB (ν), and J(ψ) are represented by [38] where m is the vehicle total mass, I z is the z-axis moment of inertia, and x g is the gravity center displacements on the x-axis.
As for the M A and C A (ν) matrices (which are related to the added mass phenomenon), they are defined by [38] where X u, Y v, N ṙ represent the decoupling parameters and the symmetry of the vessel, while the other parameters express the vessel's asymmetries.
Each element of M A represents the mass added in a specific direction.For instance, Y ṙ corresponds to the added mass affecting the y-axis due to the acceleration ṙ.Furthermore, it is defined as Finally, the hydrodynamic damping D(ν) is expressed as a combination of linear and nonlinear damping, which is directly related to quadratic damping.Thus, D(ν) is given by [38] where the elements of linear damping X u , Y v , and N r represent the uncoupled linear parameters and the vessel's symmetry, while X v , X r , Y u , Y r , N u , and N v represent its asymmetries.
For large-scale boats, the asymmetries can be disregarded in the formulation, simplifying the linear damping relationship for a given direction of movement.Additionally, in the nonlinear damping, X |u|u , Y |v|v , and N |r|r represent the nonlinear symmetrical and uncoupling damping of quadratic order in the directions of the x, y, and z axes, respectively.The remaining parameters represent the asymmetries and couplings.
Assuming that M RB and C RB (ν) are related to previously known variables, the parameter identification problem boils down to finding M A , C A (ν), and D(ν).Therefore, the solution to the parameter identification problem exists within a R 27 dimensional space or finding of the hydrodynamic derivatives.

rSOESGOPE Application
To apply the rSOESGOPE method, the first step is to obtain an initial estimate of the vessel's dynamics.Based on this estimation, benchmark parameters are designed for later use in designing identification signals, applying them to the vessel, and ultimately estimating the posteriori parameters.
Thus, the application of the method follows a flowchart with four well-defined steps: (1) Prior Estimation, (2) Generate Benchmark Parameters, (3) Signal Design, and (4) Parametric Estimation.These steps are also illustrated in Figure 4:

Priori Estimation
Two approaches were employed to obtain the initial estimation of the vessel's dynamics.First, the operational and inertial parameters were obtained through direct measurements or Computational Fluid Dynamics (CFD) simulations.The measured parameters are presented in Table 2, while the computationally simulated ones can be found in Table 3.In these tables, [F x , F y , τ n ] represent the propulsion system capacity for surge, sway, and yaw, respectively.Additionally, [t F x , t F y , t τ n ] denote the time constants related to the dynamics of surge, sway, and yaw, respectively.Subsequently, to obtain the initial values of the main hydrodynamic parameters for the simpler uncoupled model with nine variables, an experimental test was conducted using an APRBS-based signal [39].The results of this test are presented in Table 4:

Generation of Benchmark Parameters
To define it is first necessary to determine the number of characters n that make up U ⊕ .This number should be sufficient to minimize uncertainties about Γ− and ensure the subsequent robustness of the method.
For the presented case study, the use of up to five signals in the composition of U ⊕ will be investigated.Therefore, P ⊕ is also composed of five parameter sets Γp , which were generated by five random samples belonging to the uncertainty region of 50% around Γ− .The equation used is shown below: where σ ∈ R represents the level of parametric uncertainty in the range of [0, 100], and r ∈ R 9×1 is a vector of random numbers with uniform distributions ([−1, 1]).Table 5 shows the generated and required P ⊕ for the design of U ⊕ according to Section 3.

Designing Signals
Firstly, it is necessary to determine the number of stages in the signal and the metaheuristic's population size.
As shown in Section 4, the number of stages affects the elements of the individuals, i.e., the dimension of Ξ, and also the excitation's persistence property.Moreover, the number of individuals also affects the metaheuristic's ability to search for signals with high excitation persistence.
Therefore, in order to find the best configuration, a study was developed using the PSO, which includes a single-signal design.For a better signal persistence analysis and its relationship with the number of stages and the swarm size, no parametric constraints are considered.With this in mind, the following definitions were used: ERROR WEIGHTINGS: Only the velocities in the rigid-body frame [u, v, r] were considered for the optimization process, weighted by the inverse of the variance.Therefore, Q x was disregarded, and Q y was defined as follows: where σu = 0.05 [m/s], σv = 0.05 [m/s], and σr = 0.01 [rad/s] are estimated standard deviations of the measurement errors of surge, sway, and angular rate, respectively.• SIGNAL CODING: The OSD problem consists of the determination of u(Ξ k ), where Ξ k is defined as follows: where k is the stage number, [F k x , F k y , τ k n ] are the APRBS amplitudes, and [t k F x , t k F y , t k τ n ] are the APRBS swap times.Table 6 presents the simulation results for 50 iterations (referred to as '50 signals') for each configuration under study.These results are represented by the average (av) and standard deviation (std) of f (Ξ).Additionally, all simulations utilized the same sets of initial parameters Γ− (refer to Table 4) and reference Γp (as shown in Equation ( 21)).Furthermore, they employed identical initial conditions for generating the initial signals.Table 6 demonstrates an improvement in signal persistence as the number of particles increases.This improvement is noticeable from 5 to 20 particles, with significant enhancements.Additionally, it is observed that there is an enhancement in signal quality up to six stages.Beyond this threshold, the signal's persistence decreases, indicating that further increasing the signal encoding is no longer beneficial.
Consequently, following the proposed framework, the optimal configuration for this problem involves using 20 particles and six stages.
Based on these findings, five signals were generated for each metaheuristic.Essentially, each signal . Therefore, Equation ( 5) can also be applied without loss of generality as follows: where S represents the double-layer optimization methodology (referred to as "Metaheuristic + IPA") whose metaheuristics will be investigated.Additionally, parametric constraints have been applied to generate achievable signals that closely resemble the ASV behavior, resulting in k θ = 2 × 10 5 .Figures 5-12 display the generated APRBS signals.These figures exhibit the typical characteristics of projected signals, including multiple amplitudes and varying excitation frequencies.It is worth noting that, for the analyzed vessel, the total experiment time does not exceed 32 s.This observation is noteworthy because there is often a natural assumption that longer excitation signals are superior to shorter ones.However, the results suggest that this is not an absolute truth.Experiments with excessively long duration may not necessarily enhance parameter learning but rather contribute to the richness of excitation in the projected signals.Achieving this richness is possible through the optimization process, even in shorter experiments, as long as the quality of excitation is maintained.sample posed the greatest challenges for the methods, particularly for the GA, which struggled to satisfy the boundary conditions.Additionally, the last column of the table displays the total sum of f (Ξ) for all five signals.In this regard, it can be observed that the GA and AOA techniques yielded the poorest results.Conversely, GWO, PSO, BA, and MTDE produced very similar results, with PSO achieving the best outcome.

Parametric Estimation
After designing ⊕ , the final step is to apply it to the real denoted as R(Γ), and PE.In this case, PE problem aims to determine Γ ∈ R 27 by exploring the capabilities of the rSOESGOPE method, as specified in Section 5.2.It is possible to analyze the PE results using up to five optimal signals.Therefore, for each U ⊕ n , a posteriori estimate Γ+ is computed from where n represents the number of signals of U ⊕ n = [u ⊕ 1 , . . ., u ⊕ n ] and P represents the IPA responsible for solving the PE problem, as shown in Section 3.
For comparison purposes, the Γ+ obtained by each method is evaluated against the scenario provided by R(Γ, U ⊕ ).The results are presented in Table 8, where the Root-Mean Square Error (RMSE) measures the performance.

Table 8. PE performance analysis using U ⊕
n .
- The results in Table 8 reveal that using more than one signal in u ⊕ led to improvements in the final estimate of Γ+ .This improvement was observed for the GWO, WOA, PSO, SSA, BA, and MTDE methods, where three or more signals were used, resulting in enhancements.However, the same trend was not observed for the AOA and GA techniques, as the addition of signals worsened their performance with each addition to U ⊕ n .This phenomenon can be attributed to the difficulties these techniques face in the design phase, as indicated in Table 7.
Based on these results, it is also evident that the use of multiple signals led to improvements in six out of the eight metaheuristics.It is important to note that the addition of signals should be carried out cautiously and in moderation, as seen in the case of WOA, where using more than four signals was no longer beneficial.These findings suggest that adding signals without caution can either amplify already identified features or introduce new dynamics to the estimated model.
Tables 9 and 10 present the best estimates obtained based on the from Table 8.Consequently, the GWO technique employs five signals, Optimization Algorithm three signals, and PSO employs five signals.The AOA GA techniques utilize only one signal, while SSA, BA, and MTDE utilize four signals.

Model Validation
To evaluate the quality of the models obtained from each meta-heuristic, three different maneuver signals were employed ( All experimental data, including the ASV's states and signals, are available in the dataset [40]).In scenario 1 (C1), the ASV undergoes frontal motion with small yaw corrections.In scenario 2 (C2), the vessel experiences the traditional zigzag motion.Finally, in the last scenario (C3), the ASV is subjected to coupled motions.
Table 11 shows the models' performance, while the signals used and the behavior of the estimated models can be seen in Figures 13 and 14.In addition, the performance of each model M( Γ+ ) is measured by RMSE.Analyzing the experimental results and their performance in Table 11, it can be concluded that the majority of metaheuristics yielded good results.The exceptions are the AOA and BA methods, as evident in Figure 14.In this figure, it is apparent that the AOA failed to capture the primary dynamics of the ASV, resulting in a significant divergence in 3 DoFs.On the other hand, the Bat Algorithm (BA) method demonstrated strong generalization abilities for sway and yaw motions but exhibited a high degree of divergence in the surge DoF during Track 3.
As for the other techniques, it is noteworthy that the estimated models exhibited behaviors closely resembling the ASV in 3 DoFs, as illustrated in Figures 13 and 14.Among these techniques, it is worth mentioning that the GWO and MTDE methods delivered the best performances and outperformed the PSO.According to Table 11, GWO reduced the RMSE by approximately 6.21%, while MTDE reduced it by 0.69% when compared to PSO.
Another fact showing that the results obtained were satisfactory is found in Table 12, which shows the RMSE of each DoF over all experiments.
Considering that the experiments were conducted in an open environment exposed to wind, waves, and other external disturbances, the results remain consistent, except for the AOA and BA techniques.For instance, in the case of GWO, an RMSE of 0.24 [m/s] was obtained for surge, and 0.20 [m/s] for sway, which correspond to errors of 12% and 10%, respectively, when considering a linear speed operation of 2.00 [m/s].Regarding yaw movement, the RMSE error of 0.34 [rad/s] would equate to 11.85% if the ASV operated at full capacity.Hence, both results are acceptable under the presented conditions, and the same holds true for methods using multiple signals, as demonstrated in Figures 13 and 14.
From Table 12, it can be observed that the methods using multiple signals (GWO, WOA, PSO, SSA, BA, and Multi-trial vector-based Differential Evolution (MTDE)) yield better results compared to the estimation models using a single signal (AOA and GA).For example, when comparing the results of GA with those of GWO, a reduction in RMSE of 14.28% is observed in the surge DoF, and similar improvements are observed in other DoFs.

Qualitative Analysis of Algorithm Performance
Taking a closer look at the algorithmic results, it is worthwhile to compare the techniques that achieved success, referred to as Group 1 (GWO, WOA, PSO, SSA, BA, and MTDE) and Group 2 (AOA and GA), compared to techniques that did not perform well in the studied problem.
In the case of Group 2, the GA exhibited low performance.Despite its wide applicability in complex problems, it is known that the GA can face challenges related to parameterization and slow convergence, especially in complex and multidimensional problems.These challenges became evident in the final model of the algorithm, which was developed using only one optimal excitation signal.Given the complexity of the problem, this resulted in a model with limited overall learning and reduced performance.
Similarly, within Group 2, the AOA also failed to achieve satisfactory results.Premature convergence and difficulties in exploring the search space were observed.These issues can be attributed to the use of only one signal in constructing the model.Similar to the GA, this led to a model with limited overall learning.Contributing factors include the dependence on the exploration-exploitation trade-off of two parameters with fixed values, the absence of information sharing and guidance towards the best global solution, and the modification of only the best solution based on arithmetic operators during the exploration and exploitation stages of the algorithm.
Turning to the results from Group 1, it is evident that all meta-heuristics generated models using three or more sub-optimal signals.This indicates that these algorithms did not encounter issues with premature convergence or difficulties in exploring the solution space.Furthermore, they allowed for the development of models with increased learning through the addition of new signals.
Among the six techniques in Group 1, three stood out as highly efficient: GWO, and the meta-heuristics MTDE and PSO.These techniques yielded very similar results and demonstrated efficiency in exploring the search space and addressing the complexities of OSD problems.GWO achieved the best performance, which can be attributed to (i) its lower sensitivity to parameters, making it easier to configure and fine-tune compared to other algorithms, and (ii) its ability to explore a wide range of solutions within its search space before converging towards solutions closer to the global optimum.
It is important to emphasize that each algorithm possesses individual characteristics and may excel in different scenarios.In this study, GWO demonstrated superiority, even outperforming the originally used PSO in the rSOESGOPE method.Finally, for a better understanding, Table 13 below provides a comprehensive comparison of the entire study conducted in this article.

Conclusions
This paper presented a comparison of metaheuristic performance when applied to Optimal Signal Design for estimating parameters of a nonlinear system.In particular, various metaheuristics were investigated and applied in the context of rSOESGOPE, which originally utilized the PSO algorithm as the primary optimizer.To assess the performance of these metaheuristics, a modeling problem of a catamaran-type ASV was employed.From the results, several key points can be emphasized:

•
ASV was chosen because it presents several difficulties in the scope of modeling and identification, which must deal with geometric and mechanical asymmetries, coupling between the DoFs, a high level of uncertainty about the initial estimates, and other aspects such as environmental disturbances.• Furthermore, these challenges also include the design of signals in large search spaces and the parametric estimation of many parameters.

•
The techniques were analyzed at each stage of the rSOESGOPE methodology, from configuration, signal design, and parameter estimation to the model validation obtained in different identification scenarios.

•
Considering the Optimal Signal Design phase, it has been proposed and demonstrated that conducting the study to define the characteristics of the APRBS signal plays a fundamental role in the projected signal quality.This study shows the importance of defining the trade-off between the persistence of the excitation and the problem solvability.Although long multi-stage APRBS signals provide richer excitation characteristics, the nonlinear solution space grows exponentially, which can make it difficult, costly, or even intractable.

•
In general, the performance of the estimated model improved with the insertion of signals, minimizing the uncertainties of the initial estimates, as initially proposed.
Exceptions occurred only for the GA and AOA models, reaffirming the difficulties in the design stage regarding compliance with the established constraint restrictions and improvement of the persistence of excitation.There is a high similarity between the estimated models and the real system.In these scenarios, the GWO method presents the best RMSE among all the methods used, followed by the MTDE and PSO methods, respectively.It can also be observed that the AOA method could not find a solution that satisfies the system, as its results were about 5 times larger than the other ones.Furthermore, regarding the RMSE for each DoF, the GWO obtained the best results considering surge and sway motions and is surpassed in the yaw behavior only by the SSA method, which is still ranked only in fifth place among the metaheuristics' results.
Finally, it is possible to appreciate the validity of the rSOESGOPE concept, which has proven its effectiveness with various algorithms.Regardless of the adopted metaheuristic, it was possible to mitigate the impact of initial uncertainties on the system and signal design, leading to the estimation of robust models.This achievement opens the door to the development of robust control strategies and fault-tolerant systems for the use of ASV in more challenging and hostile environments.
As part of future work and research, the aim is to conduct a more comprehensive assessment of the algorithms' performance under various environmental conditions and in harsher scenarios (e.g., fast-flowing water, waves, and/or winds) to gain a more comprehensive understanding of their applicability and limitations.
Furthermore, a future project involves the creation of a hybrid metaheuristic model that combines the strengths of the top-performing algorithms.This could lead to a more robust and efficient solution, further enhancing the capabilities of our models and their adaptability to a broader range of challenges.
A schematic video demonstrating the application of the technique is available at https://www.youtube.com/watch?v=OG0gr67YMNg, accessed on 11 October 2023.

Figure 3 .
Figure 3. Real view of the developed catamaran.

Table 1 .
Parameter values for the comparative algorithms.

Table 6 .
Analysis of number of stages and individuals.

Table 7
presents the properties of the projected signals across all studied metaheuristics.

Table 7
presents the performance of the generated signals for each Γp i .It is also evident that the reference parameters used impose varying degrees of difficulty.Notably, the Γp

Table 12 .
Analysis of Root-Mean Square Error (RMSE) results for each Degree of Freedom (DoF).