Next Article in Journal
Ulam Stability for Boundary Value Problems of Differential Equations—Main Misunderstandings and How to Avoid Them
Previous Article in Journal
LAMBERT: Leveraging Attention Mechanisms to Improve the BERT Fine-Tuning Model for Encrypted Traffic Classification
Previous Article in Special Issue
Dynamical Sphere Regrouping Particle Swarm Optimization: A Proposed Algorithm for Dealing with PSO Premature Convergence in Large-Scale Global Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparative Analysis of Metaheuristic Algorithms for Enhanced Parameter Estimation on Inverted Pendulum System Dynamics

by
Daniel Sanin-Villa
1,
Miguel Angel Rodriguez-Cabal
1,
Luis Fernando Grisales-Noreña
2,
Mario Ramirez-Neria
3,* and
Juan C. Tejada
4
1
Department of Mechatronics and Electromechanics, Instituto Tecnológico Metropolitano, Medellín 050034, Colombia
2
Department of Electrical Engineering, Faculty of Engineering, Universidad de Talca, Curicó 3340000, Chile
3
InIAT Instituto de Investigación Aplicada y Tecnología, Universidad Iberoamericana, Ciudad de México 01219, Mexico
4
Artificial Intelligence and Robotics Research Group (IAR), Universidad EIA, Envigado 055428, Colombia
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(11), 1625; https://doi.org/10.3390/math12111625
Submission received: 24 April 2024 / Revised: 11 May 2024 / Accepted: 15 May 2024 / Published: 22 May 2024

Abstract

:
This research explores the application of metaheuristic algorithms to refine parameter estimation in dynamic systems, with a focus on the inverted pendulum model. Three optimization techniques, Particle Swarm Optimization (PSO), Continuous Genetic Algorithm (CGA), and Salp Swarm Algorithm (SSA), are introduced to solve this problem. Through a thorough statistical evaluation, the optimal performance of each technique within the dynamic methodology is determined. Furthermore, the efficacy of these algorithms is demonstrated through experimental validation on a real prototype, providing practical insights into their performance. The outcomes of this study contribute to the advancement of control strategies by integrating precisely estimated physical parameters into various control algorithms, including PID controllers, fuzzy logic controllers, and model predictive controllers. Each algorithm ran 1000 times, and the SSA algorithm achieved the best performance, with the most accurate parameter estimation with a minimum error of 0.01501 N m and a mean solution error of 0.01506 N m. This precision was further underscored by its lowest standard deviation in RMSE (1.443 99 × 10−6 N m), indicating remarkable consistency across evaluations. The 95% confidence interval for error corroborated the algorithm’s reliability in deriving optimal solutions.

1. Introduction

The inverted pendulum symbolizes the equilibrium between stability and instability within control systems. It finds real-world uses in robotics, transportation, and automation industries [1]. However, effectively managing inverted pendulum systems entails gauging their characteristics, a task made difficult by their intricate and nonlinear nature. Conventional approaches to parameter estimation may not always yield results [2].

1.1. Metaheuristic Algorithms

Over the years, these algorithms have emerged to address this problem. Inspired by phenomena and social dynamics, these algorithms offer search methods that cover solution spaces. Through applying optimization principles, metaheuristic algorithms support parameter estimation for managing inverted pendulums, leading to stability, robustness, and effectiveness. Exploring the domain of algorithms for parameter estimation in controlling inverted pendulums requires a combination of control theory, optimization tactics, and computational intelligence [3].
Different kinds of algorithms can be found that are centered on the optimization of mechanical and engineering design problems. Due to the nature of these types of problems, which are nonlinear non-convex mathematical models, the main kinds of algorithms used are the metaheuristic techniques such as the particle swarm optimization (PSO) [4], genetic algorithm (GA) [5], and salp swarm algorithms (SSA) [6], among others. For instance, the work developed by Moraes et al. uses an enhanced PSO algorithm for parameter estimation of water in oil emulsification in a duct. The methodology implemented can obtain good-quality solutions, where the objective function describes multiple variables in a non-convex mathematical model—showing the PSO as an efficient tool [7]. Other works, in which the PSO was used, have been carried out by Liu et al. [4], where the authors present a novel variation of the PSO using the Chebyshev functions to improve the algorithm’s efficiency. The PSO was tested with EC2017 and CEC2022 benchmark functions, where the objective is the dynamic parameter estimation of different mathematical problems related to control operations. Additionally, Rodriguez et al. implemented a CGA algorithm to minimize the weight of a drive shaft, solving the mathematical model of the nonlinear nonconvex type [8]. On the other hand, Ab Wahab et al. implemented a CGA algorithm for mobile robot global route planning, where the authors developed a set of constrictions related to blockages in the routes to allow the algorithm to find the optimal path for each route [9].

1.2. Applications and Challenges in Control Systems

Inverted pendulum systems find applications across automotive, aerospace, manufacturing, and robotics sectors. Within the field of robotics, these systems play roles in navigation, manipulation, and surveillance tasks. Utilizing the dynamics of inverted pendulums, mobile robots and humanoid platforms can attain agility, maneuverability, and adaptability when operating in environments [10]. Inverted pendulum control systems are pivotal in maintaining vehicle stability, ensuring safe navigation, and improving overall performance and safety in autonomous vehicles, especially in challenging terrains or adverse weather conditions. Accurate parameter estimation techniques such as mass distribution, moment of inertia, and friction coefficients enable precise control and trajectory tracking [11]. In industrial automation, inverted pendulum systems optimize the operation of conveyor belts, robotic arms, and other automated machinery to improve efficiency, throughput, and reliability in production processes [12]. The aerospace industry also utilizes this type of control system in designing and operating unmanned aerial vehicles (UAVs) and spacecraft, enhancing UAVs’ maneuverability, stability, and payload capacity for tasks such as surveillance, reconnaissance, and aerial photography. In rehabilitation robotics, inverted pendulum-based systems assist individuals with mobility impairments or neurological disorders. By providing dynamic balance support and gait assistance, these systems aid in rehabilitation therapy, helping patients regain mobility, improve their motor control, and enhance their overall quality of life. Accurate parameter estimation is crucial for achieving optimal performance and functionality across diverse applications of inverted pendulum systems [13,14].

1.3. Relevance of Torque Prediction in Mechanical Systems

Torque is a fundamental parameter in determining a mechanical system’s dynamic behavior and stability. An accurate understanding of torque through dynamic modeling is essential for developing effective control strategies. An inadequate torque prediction can compromise a system’s efficiency, reliability, and safety. Therefore, dynamic modeling enables engineers to understand how a system responds to external inputs and disturbances. Accurately predicting the torque on mechanical systems is extremely important in control engineering. This prediction significantly affects the control system’s design, operation, and performance in various applications [15]. However, without precise torque predictions, the dynamic model may mischaracterize the system’s behavior, leading to ineffective control strategies. Control algorithms such as PID (Proportional–Integral–Derivative) controllers, model predictive controllers, and adaptive control schemes rely heavily on accurate torque predictions [16]. These predictions help generate control actions that stabilize the system and achieve the desired performance objectives. Estimating the torque allows controllers to efficiently manage the system performance, adjusting for disruptions and uncertainties as they occur. This involves reducing energy usage, lessening strain, and optimizing effectiveness [17]. Engineers can optimize control strategies by accurately estimating the torque, ensuring the system operates safely while attaining the desired performance metrics. Where energy-efficient control schemes can be implemented, they prolong the mechanical components’ lifespan [18]. In safety-critical applications such as automotive systems, robotics, and industrial machinery, accurate torque prediction is imperative for ensuring the safety and reliability of the controlled system. Overestimation or underestimation of the torque can lead to unexpected system behavior, potentially resulting in equipment failure, operational accidents, or injury to personnel. Control systems can implement appropriate safety measures and ensure reliable operation by accurately predicting the torque under varying conditions. Optimizing performance metrics like speed, accuracy, and precision requires a deep understanding of torque dynamics. Accurate torque prediction equips control systems to achieve optimal performance while adhering to operational constraints and performance specifications [19].
Researchers in control engineering continually strive to enhance the accuracy and efficiency of torque prediction in mechanical systems for control purposes. Advanced modeling techniques and optimization algorithms have emerged as promising avenues for improving torque prediction accuracy and optimizing control strategies. However, existing studies have identified limitations and challenges associated with current approaches, necessitating further exploration and refinement of methodologies. A torque-based nonlinear predictive control approach has been proposed for automotive powertrains, highlighting the significance of dynamic torque control [20]. This approach underscores the importance of accurate torque estimation for optimizing powertrain performance and efficiency. Similarly, sensorless control techniques have been developed for synchronous motor drives to achieve resilient torque control, particularly in the presence of parameter errors [21,22]. However, limitations such as applicability under large parameter errors and accuracy constraints at lower speeds have been noted, underscoring the need for robust torque estimation methods. Nonlinear model predictive control (NMPC) has garnered attention for internal combustion engines, offering enhanced performance through advanced modeling and optimization [23]. Yet, challenges remain, including potential tracking errors induced by predictive model inaccuracies. Similarly, real-time model predictive control (MPC) of drive systems with elastic transmission has shown promise but faces limitations in addressing severe torsional vibrations and load disturbances [24]. Metaheuristic optimization techniques have been explored for tuning proportional–integral controllers in drive systems, demonstrating their potential to improve the control efficiency [25]. However, further investigation is necessary to explore their applicability across different control scenarios and systems. Furthermore, the optimization-based motion prediction of mechanical systems has been investigated, highlighting the benefits of sensitivity analysis for optimizing lifting motion. The ongoing quest to enhance torque prediction accuracy and control efficiency underscores the interdisciplinary nature of control engineering. By addressing the limitations and challenges identified in the current research, there is a way to develop more robust and effective control strategies for various mechanical systems. This study aims to contribute to this topic by investigating novel torque prediction using optimization methodologies for parameter estimations, drawing upon insights from advanced modeling techniques, optimization algorithms, and practical applications in various domains. In the pursuit of advancing torque prediction and control methodologies, several seminal studies have laid the groundwork for current research. For instance, a torque-based nonlinear predictive control approach has been proposed for automotive powertrains, emphasizing the importance of dynamic torque control. Despite its potential, challenges remain, including the risk of tracking errors due to predictive model inaccuracies. Innovative approaches leveraging metaheuristic optimization techniques have been explored for optimizing control strategies in drive systems. While these approaches can improve control efficiency, their applicability across different scenarios and systems warrants further investigation [26]. Integrating advanced modeling techniques, optimization algorithms, and practical applications represents a fertile ground for advancing torque prediction and control strategies. Using these strategies, engineers can enhance performance and efficiency across various mechanical systems [27].

1.4. Contributions and Research Objective

The findings of this study have implications for advancing control systems engineering by providing solutions to improve the performance and reliability of inverted pendulum-based applications across different fields. Using algorithms, we aim to address the challenges related to parameter estimation in a mathematical model of inverted pendulum systems, for more resilient and adaptable control strategies in real-world settings. This research is based on a range of literature exploring methods and approaches to find the physical parameters to estimate the torque in the mechanism, and three optimization techniques are proposed: the PSO, CGA, and SSA, where a statistical evaluation was carried out to obtain the best evaluation and response of each technique attained to the dynamic methodology.

2. Methodology

This section presents the methodological strategies employed in estimating the model parameters for a mechanism using experimental data. The model parameters were effectively ascertained through the experimental measurements from an actual implementation of the mechanism with a metaheuristic algorithm.

2.1. Inverted Pendulum System

The inverted pendulum system comprises a planar double inverted pendulum. A DC motor drives the initial link, whereas the subsequent link is an underactuated simple pendulum. θ 1 and θ 2 represent the angular positions of the links, u stands for the control torque input, while m 1 and m 2 are the masses of the links. The lengths of the links are denoted by l 1 and l 2 , and l c 1 and l c 2 indicate the distances to the center of masses. Finally, I 1 and I 2 stand for the inertia of the links. The system is shown in Figure 1.

Mathematical Model

Underactuated Euler–Lagrange systems of fourth order can be generally represented as [28]
P ( q ( t ) ) u ( t ) = M ( q ( t ) ) q ¨ ( t ) + C ( q ( t ) , q ˙ ( t ) ) q ˙ ( t ) + g ( q ( t ) ) .
Here, q ( t ) = q 1 ( t ) q 2 ( t ) denotes the generalized coordinates, M ( q ( t ) ) R 2 × 2 represents the inertia matrix, and C ( q ( t ) , q ˙ ( t ) ) R 2 × 2 describes the Coriolis and centrifugal forces. The vector g ( q ( t ) ) R 2 represents the gravitational forces, and P ( q ( t ) ) R 2 maps the external forces. Additionally, u ( t ) R signifies the control input.
Such systems, described by (1), can be alternatively expressed in state space as
x ˙ a ( t ) = x b ( t ) x ˙ b ( t ) = f ( x a ( t ) , x b ( t ) ) + ζ ( x a ( t ) ) u ( t ) ,
where x a ( t ) represents the joint positions vector q 1 ( t ) q 2 ( t ) , x b ( t ) denotes the articular velocities vector q ˙ 1 ( t ) q ˙ 2 ( t ) , and
f ( x a ( t ) , x b ( t ) ) = M ( q ( t ) ) 1 C ( q ( t ) , q ˙ ( t ) ) q ˙ ( t ) + g ( q ( t ) )
ζ ( x a ( t ) ) = M ( q ( t ) ) 1 P ( q ( t ) ) .
Utilizing the Euler–Lagrange formalism, the dynamic model of the system can be represented as
M ( q ( t ) ) q ¨ ( t ) + C ( q ( t ) , q ˙ ( t ) ) q ˙ ( t ) + g ( q ( t ) ) = P τ ( t ) ,
where
q ( t ) = θ 1 ( t ) θ 2 ( t )
P = 1 0
M ( q ) = β 1 + β 2 + 2 β 3 cos ( θ 2 ( t ) ) M 1 , 2 ( θ ) β 2 + β 3 cos ( θ 2 ( t ) ) β 2
C ( q , q ˙ ) = β 3 θ ˙ 1 ( t ) sin ( θ 2 ( t ) ) β 3 C 1 , 2 ( θ ) β 3 θ ˙ 1 ( t ) sin ( θ 2 ( t ) ) 0
g ( q ) = β 4 g cos ( θ 1 ( t ) ) + β 5 g cos ( θ 1 ( t ) + θ 2 ( t ) ) β 5 g cos ( θ 1 ( t ) + θ 2 ( t ) ) ,
where
β 1 = m 1 l c 1 2 + m 2 l c 2 2 + I 1
β 2 = m 2 l c 2 2 + I 2
β 3 = m 2 l 1 l c 2
β 4 = m 1 l c 1 + m 1 l 1
β 5 = m 2 l c 2 .

2.2. Experimental Setup

The experimental pendulous prototype is displayed in Figure 2. It has a DC motor model NC5475, manufactured by NISCA, Tokio, Japan; to drive the first link. The angular positions of both links are monitored with incremental encoders of 10,000 counts per revolution, consistent with the earlier example. Power amplification is facilitated by the amplifier model VoltPAQ-X2, manufactured by Quanser, Markham, Canada. The control strategy is implemented in the Matlab–Simulink platform, with a sampling time of 0.001   s .
The actual physical parameters of the prototype are outlined below and were measured before the assembly:
  • Link inertias: I 1 = 0.00053 [kg·m2] and I 2 = 0.00077 [kg·m2].
  • Masses of the links: m 1 = 0.210 [kg] and m 2 = 0.1 [kg].
  • Lengths of the links: l 1 = 0.15 [m] and l 2 = 0.3 [m].
  • Distances to the center of mass: l c 1 = 0.12 [m] and l c 2 = 0.15 [m].
  • Armature resistance: k τ = 0.0724 [ Ω ].
  • Torque constant: τ m = 2.983 [ Ω ].
The provided values do not encompass any joints, glue, screws, or any accessories that are not part of the model. To refine the model, it is necessary to compute the equivalent values of mass, center of mass location, and inertia for each link. To assess the dynamic model, a trajectory is suggested with the initial conditions set as x 1 ( 0 ) = π 2 [rad] and x 3 ( 0 ) = 0 [rad]. The trajectory involves a seamless transition from rest to rest.
y f * ( t ) = β 2 + β 3 β 2 x δ 1 * ( t ) + x δ 3 * ( t ) y f * ( t ) = β 2 + β 3 β 2 ( θ 1 * ( t ) π 2 ) + θ 2 * ( t )
Initially, at t = 0 , the trajectory starts from y f * ( 0 ) = 0 , setting θ 1 * ( t ) = π 2 and θ 2 * ( t ) = 0 . At t = 4.5 , it moves to y f * ( 6 ) = β 2 + β 3 β 2 π 6 π 6 with θ 1 * ( t ) = 3 2 π and θ 2 * ( t ) = π 6 over a duration of 2.5 seconds. Subsequently, at t = 9.5 [s], it transitions to y f * ( 13.5 ) = β 2 + β 3 β 2 π 6 + π 6 with θ 1 * ( t ) = π 3 and θ 2 * ( t ) = π 6 within 4 s. Finally, at t = 17.5 [s], it returns to the initial position until the test concludes. Figure 3 illustrates the rest-to-rest positions of the pendulum.
The velocity and angular acceleration for each link will be estimated from the position values derived from the smooth trajectory. These kinematic features will inform the evaluation of the prototype’s dynamic behavior, facilitating the acquisition of the experimental torque data necessary to proceed with the parametric estimation.
The experimental torque is estimated based on the motor’s angular velocity ( ϕ ( t ) ˙ ) and the supply voltage ( V ( t ) ), expressed as
τ = k 1 R a V ( t ) k 1 k 2 R a ϕ ˙ .
Here, k 1 represents the torque constant, k 2 denotes the back-EMF constant, and R a stands for the winding resistance.

2.3. Algorithm Implementation

To solve the optimization problem addressed, it is necessary to determine the limits of the optimization parameters, which define the solution space to the algorithms. These upper and lower limits ensure the algorithm produces accurate and relevant results. By defining these limits, the algorithm can effectively narrow down the range of potential solutions, resulting in faster and more efficient problem solving. The data presented in Table 1 provide the upper and lower limits for the model parameters to be estimated, with consideration of the physical constraints. These limits were based on the estimated actual parameters and were determined by ensuring that the center of mass distance was below the length of the link.
To realize the optimization problem, the root mean square error is used as a fitness function (Equation (18)); this calculates the square root over the average of the squared differences between the predicted values of the torque and the actual (experimental) values. This allows evaluation of how closely the model’s predictions match the real-world data [29], where the main objective is to minimize this function, which means that the mathematical model with the estimated parameters matches the experimental values. To quantify the results obtained by the algorithms, the uncertainty associated with the mean solution errors is evaluated, where the confidence interval of 95% is calculated for each algorithm’s mean error across the 1000 simulation runs. The confidence interval provides a range in which the true mean error is expected with a 95% probability.
The mean solution is calculated by a set of RMSE measurements (Equation (18)) for each algorithm. n is the number of observations, τ ^ i is the predicted torque value for the i-th observation, and τ i is the experimental value for the i-th observation.
R M S E = 1 n i = 1 n ( τ ^ i τ i ) 2
The RMSE is a useful fitness function for assessing the effectiveness of a developed model with the estimated parameters. A lower RMSE implies that the model’s predictions are closer to the actual values, which is what we desire. Conversely, a higher RMSE indicates significant discrepancies between the predictions and actual values, resulting in poor parameter estimation on the model performance [30]. Therefore, a lower RMSE denotes more accurate predictions, reflecting a better parameter estimation.
The standard deviation s, is calculated with Equation (19):
s = 1 n 1 i = 1 n ( R M S E i R M S E ¯ ) 2 .
The standard error of the mean (SEM) is then given by Equation (20):
S E M = s n .
We use the standard error and the Z-score corresponding to the desired confidence level to calculate the 95% confidence interval for the mean error. The Z-score is a statistical measure that indicates how many standard deviations an element is away from the mean of a distribution. For a two-tailed test at a 95% confidence level, the Z-score is approximately 1.96.
Thus, the confidence interval (CI) is calculated as follows:
C I = R M S E ¯ ± ( Z · S E M ) .
In this study, we applied this methodology to compute the confidence intervals for the mean error obtained after applying the three proposed algorithms: SSA, CGA, and PSO. The computed intervals provide insight into the precision of the mean error estimates and allow us to ascertain the statistical significance of the differences observed between the algorithms’ performances in determining the parameters. The algorithms that determine the parameters of the model are presented in Table 1. Their objective is to minimize the RMSE, finding the specific parameters that close the gap between the predictive torque and the experimental measurement, as shown in Equation (18). Each metaheuristic algorithm proposed in this work is described below. Each algorithm setting was optimized using a stochastic PSO algorithm to obtain the best parameters and improve the performance for solving the optimization problem.

2.4. Particle Swarm Optimization

The PSO is a metaheuristic technique for optimizing non-convex nonlinear mathematical models. It is a bio-inspired algorithm that simulates the way flocks of birds sweep a terrain in search of food, where each animal is modeled as a particle. Each particle in the swarm has a cognitive and social component, allowing the swarm to follow the best particle that finds the best solution in the current iteration. Adding a velocity that changes in each iteration depending on its position allows the exploration of the entire solution space [4,31]. Algorithm 1 shows the process to solve an optimization problem. The RMSE is often employed as a fitness function in optimization scenarios to reduce the disparity between the predicted and observed values. The selection of the RMSE as the fitness metric does not directly impact how the optimization techniques perform; rather, it consistently gauges the error these algorithms strive to diminish. For instance, the PSO algorithm utilizing the RMSE offers a continuous error gauge, steering the particles toward the possible solution. The algorithm fine-tunes particles’ positions to minimize the RMSE. This choice does not inherently influence how well these optimization techniques function; instead, they leverage the RMSE to direct their quest for a solution. The efficacy of each technique in minimizing the RMSE depends on its mechanisms and parameter configurations.
Algorithm 1: Particle swarm optimization procedure.
Mathematics 12 01625 i001

2.5. Continuous Genetic Algorithm

The genetic algorithm is a classic optimization technique widely used to solve continuous problems with nonlinear mathematical models. The algorithm starts by creating an initial population of individuals, each evaluated based on the objective function and a set of constraints representing the problem’s physics [5]. These enable the evaluation of the response viability, allowing the algorithm to consider non-feasible points as potential solutions. This facilitates better exploration and exploitation of the solution space, preventing the algorithm from getting stuck in local optima. Subsequently, descendant populations are generated using recombination, selection, and mutation techniques. The objective function and set of restrictions are evaluated in each iteration, allowing for advancement through the solution space to find an adequate and viable solution to the analyzed problem [8]. Algorithm 2 briefly details the CGA process to solve an optimization problem.
Algorithm 2: Continuous genetic algorithm procedure.
Mathematics 12 01625 i002

2.6. Salp Swarm Algorithm

The salp swarm algorithm is a bio-inspired metaheuristic technique with key factors to avoid being trapped in local optima, allowing the exploration and exploitation of the solution space [6,30]. This algorithm is based on swarm intelligence, which simulates the behavior of salps in the ocean, which move in nature in the form of a chain, where the first individual is taken as the leader particle, and the rest of the salps are the followers, where in each iteration, the search is carried out for the area in which the food is found, which represents the optimal zone; then the movement of the leader particle towards the zone is carried out, solving the mathematical model of the various points where the followers are located, to update the positions of the leader particle, where the objective is that all the particles approach the area of the best food, finding the best solution to the problem [32]. Algorithm 3 briefly details the SSA process to solve an optimization problem.
Algorithm 3: Salp swarm algorithm procedure.
Mathematics 12 01625 i003

3. Results

Figure 4 presents the experimental torque results correlated with the position, velocity, and acceleration data outlined in Section 2.2.
To provide a clear and direct comparison of the results on the parameter estimation for the mechanism, Figure 5 shows the plots of the torque predicted using the mathematical model with the parameters found by the different algorithms. It also compares the experimentally measured torque with the theoretical model alongside the predictions made by the SSA, CGA, and PSO algorithms. A zoomed view box gives a detailed comparison for a specific time interval between 5 s and 6.6 s, allowing for closer scrutiny of the algorithms’ performance in relation to the experimental data.
Table 2 shows the parameter estimation found by each algorithm; those values correspond to the minimum RMSE over the 1000 runs.
The proposed mathematical model was validated by comparing it with the experimental torques. The overlapping lines show that the model closely aligns with the experimental observations, which serves as a benchmark for the algorithms. The estimated torques of the SSA, CGA, and PSO algorithms appear as distinct lines. The SSA estimation adheres most closely to the experimental and theoretical benchmarks, indicating that the SSA’s parameter optimization can accurately calculate the system dynamics over the observed time. The CGA and PSO estimations generally follow the same pattern; however, their divergence from the experimental and theoretical lines in different points indicates less precision in the model fitting.
The performance of the SSA algorithm is consistent, as indicated by its proximity to the experimental and theoretical lines shown in Figure 5. This aligns with the low mean RMSE and small confidence interval mentioned earlier, which confirms that it is a robust optimization tool. On the other hand, the CGA and PSO lines exhibit greater variance from the expected torque values, particularly in areas where sharp changes in torque occur. This behavior could be attributed to these algorithms’ exploration and exploitation mechanisms, which may not be as well-tuned as the SSA algorithm to capture the intricacies of the modeled system.
A comparative result for the performance of the proposed algorithms, with value numbers, is presented in Table 3 and Table 4. The SSA algorithm achieved the best parameter estimation for the mechanism with a minimum error at 0.015015147   N m compared to the experimental results. This algorithm achieved a mean solution error of 0.015066207   N m , the lowest among the tested algorithms. The standard deviation RMSE for SSA was also the lowest, 1.44399 × 10−6 N m, indicating highly consistent results across the 1000 runs. The 95% confidence interval for the error was narrow, further confirming the reliability of the SSA in finding the optimal solutions. This algorithm demonstrated a moderate mean computing time of 916.5614   s , with a standard deviation of 18.7228   s and a closed confidence interval that suggests consistent performance in terms of time efficiency.
The CGA also offered interesting results, with the best solution RMSE recorded at 0.015039607   N m . However, the error of the mean solution was slightly higher at 0.015490795   N m , with the standard deviation of error being more considerable, implying less consistency in obtaining the best-fit parameters compared to the SSA. The error’s confidence interval was wider than the SSA’s, implying a larger CGA performance variation. The algorithm also had a higher mean computing time of 954.1851 s, with an identical standard deviation to the SSA (18.7228 s). The PSO, although the fastest with a mean computing time of 888.4771   s and a relatively low standard deviation of 17.4958   s , had a higher mean RMSE solution, the largest among the three algorithms, and the confidence interval for the mean RMSE was notably broad. These values indicate a significant variation in PSO’s results’ accuracy, making it less reliable for parameter estimation despite its computational speed.
The bar chart in Figure 6 provides a clear visual comparison of the average computing time required by the three metaheuristic algorithms. The PSO demanded the least amount of computing time on average, indicating a more time-efficient approach in parameter estimation, with the average time being noticeably less than that of the other two algorithms. The confidence limits in the bars suggest consideration of the variability in computing time, indicating the reliability and consistency of each algorithm’s performance over the 1000 runs. While the PSO led in terms of time efficiency, the SSA algorithm followed as a close second, displaying a marginally higher average computing time. According to the analysis, the average CGA algorithm took the longest time to compute, which may not be ideal when time is critical. Moreover, the analysis indicates an inverse relationship between the precision of the algorithms discussed in detail in Table 4. This implies that there may be a trade-off situation where the CGA algorithm, despite its longer computation time, may not necessarily produce better solutions, as demonstrated by its mean solution error compared to the SSA and PSO algorithms.
Figure 7 presents a comparative analysis of the average RMSE for the applied metaheuristic algorithms. The bar chart reveals that the SSA algorithm has the lowest average RMSE, indicating the highest level of precision in fitting the experimental measurements to the mathematical model; it also has a comparatively lower height of the SSA confidence interval limits with a confidence of 95%. These error bars indicate variability in the accuracy of the algorithms from one run to another, with the SSA showing the least variability, followed by the CGA and PSO. The CGA algorithm, while surpassing the performance of the PSO, exhibits a higher average RMSE than the SSA, as visualized by the taller height of its corresponding bar. This suggests the CGA is less precise than the SSA but more so than the PSO. The error bars for the CGA and PSO are noticeably longer, particularly for the PSO, which signals a greater spread in the RMSE values and, consequently, less consistent performance. The PSO algorithm’s bar, which is the tallest and the longest error bar, indicates that it has the poorest response for this problem.

4. Conclusions

A parameter estimation via dynamic modeling of an inverted pendulum system was carried out, where the comparison between the mathematical model and the experimental measurements of the torque express a correlation that shows the accuracy of the methodology described.
To solve the mathematical model for the parameter estimation, three algorithms were implemented, where the best solution was achieved by the SSA, which was 0.163% better than the CGA and 0.2% than the PSO. Regarding the repeatability of the solution, the SSA also presented the best standard deviation, reflecting a better quality of solution compared with the CGA and the PSO. Finally, in terms of the computing time, the PSO algorithm took 7.013% less time to obtain the response compared with the SSA and CGA. The SSA was selected as the best algorithm, since, for the problem addressed, the quality of the response has more influence.
This research introduces an innovative methodology and empirically validates dynamic modeling for parameter estimation in inverted pendulum systems. By integrating metaheuristic algorithm techniques and a comparative analysis of their performance, this study improves the precision of the modeling of complex dynamic systems. In future work, it is proposed that these parametric estimation methods be used to determine gains and control parameters of various linear and nonlinear techniques to control dynamic systems like mobile robots, drones and other robotic systems. Furthermore, a comparison of the effectiveness of metaheuristic methodologies with other algorithms to obtain nonlinear models will be studied, such as in MPC.

Author Contributions

The authors contributions to the achievement of the research are as follows: Conceptualization, D.S.-V.; methodology, L.F.G.-N.; software, M.A.R.-C. and J.C.T.; validation, M.R.-N., J.C.T. and M.A.R.-C.; formal analysis, D.S.-V.; investigation, D.S.-V., M.A.R.-C. and J.C.T.; resources, M.R.-N.; data curation, J.C.T.; writing—original draft preparation, J.C.T. and M.A.R.-C.; writing—review and editing, D.S.-V. and L.F.G.-N.; visualization, D.S.-V.; supervision, M.A.R.-C.; project administration, D.S.-V.; funding acquisition, M.R.-N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by InIAT Instituto de Investigación Aplicada y Tecnología Universidad Iberoamericana Ciudad de México under Project DINVP-051 and Universidad EIA under Project INV-CO-049-2023.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, M. Ramirez-Neria ([email protected]), upon reasonable request.

Acknowledgments

Thanks to Instituto Tecnológico Metropolitano, Universidad de Talca, InIAT Instituto de Investigación Aplicada y Tecnología-Universidad Iberoamericana and IAR Research Group-Universidad EIA.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nayak, A.; Banavar, R.N. Almost-global tracking of the unactuated joint in a pendubot. IFAC Pap. OnLine 2018, 51, 137–142, 6th IFAC Workshop on Lagrangian and Hamiltonian Methods for Nonlinear Control LHMNC 2018. [Google Scholar] [CrossRef]
  2. Yan, Z.; Lai, X.; Meng, Q.; Wu, M.; She, J.; Iwasaki, M. Modeling, analysis, and adaptive neural modified-backstepping control of an uncertain horizontal pendubot with double flexible joints. Control Eng. Pract. 2023, 139, 105647. [Google Scholar] [CrossRef]
  3. Anh, H.P.H.; Son, N.N.; Van Kien, C.; Ho-Huu, V. Parameter identification using adaptive differential evolution algorithm applied to robust control of uncertain nonlinear systems. Appl. Soft Comput. 2018, 71, 672–684. [Google Scholar] [CrossRef]
  4. Liu, H.; Wang, W.; Cheng, X.; Zheng, H. Particle swarm optimization with Chebychev functional-link network model for engineering design problems. Appl. Soft Comput. 2022, 131, 109819. [Google Scholar] [CrossRef]
  5. Milenković, B.; Jovanović, D.; Krstić, M. Application of Particle Swarm Optimization for Classical Engineering Problems. IJEEC Int. J. Electr. Eng. Comput. 2020, 5, 42–49. [Google Scholar] [CrossRef]
  6. Medjahed, S.A.; Ouali, M. A new hybrid SSA-TA: Salp Swarm Algorithm with threshold accepting for band selection in hyperspectral images. Appl. Soft Comput. 2020, 95, 106534. [Google Scholar] [CrossRef]
  7. Moraes, A.O.; Mitre, J.F.; Lage, P.L.; Secchi, A.R. A robust parallel algorithm of the particle swarm optimization method for large dimensional engineering problems. Appl. Math. Model. 2015, 39, 4223–4241. [Google Scholar] [CrossRef]
  8. Rodriguez-Cabal, M.A.; Grisales-Norea, L.F.; Maŕn, J.A.; Montoya, O.D. Optimal design of transmission shafts: A continuous genetic algorithm approach. Stat. Optim. Inf. Comput. 2019, 7, 802–815. [Google Scholar] [CrossRef]
  9. Ab Wahab, M.N.; Nazir, A.; Khalil, A.; Ho, W.J.; Akbar, M.F.; Noor, M.H.M.; Mohamed, A.S.A. Improved genetic algorithm for mobile robot path planning in static environments. Expert Syst. Appl. 2024, 249, 123762. [Google Scholar] [CrossRef]
  10. Turrisi, G.; Carlos, B.B.; Cefalo, M.; Modugno, V.; Lanari, L.; Oriolo, G. Enforcing Constraints over Learned Policies via Nonlinear MPC: Application to the Pendubot. IFAC Pap. OnLine 2020, 53, 9502–9507, 21st IFAC World Congress. [Google Scholar] [CrossRef]
  11. Cychowski, M.; Szabat, K. Model predictive speed control with optimal torque constraints handling of drive systems with elastic transmission. In Proceedings of the 2009 IEEE International Electric Machines and Drives Conference, Miami, FL, USA, 3–6 May 2009; pp. 251–258. [Google Scholar] [CrossRef]
  12. Khalil, W.; Boyer, F. An efficient calculation of computed torque control of flexible manipulators. In Proceedings of the 1995 IEEE International Conference on Robotics and Automation, Nagoya, Japan, 21–27 May 1995; IEEE: New York, NY, USA; pp. 609–614. [Google Scholar] [CrossRef]
  13. Du, C.; Yan, F.; Yan, Y.B.; Yang, P.L. Methods of engine torque estimation for control algorithms. Trans. Csice 2008, 26, 446–451. [Google Scholar]
  14. Leuer, M.; Bocker, J. Self-optimizing Model Predictive Direct Torque Control for electrical drives. Proccedings of the 2015 IEEE 24th International Symposium on Industrial Electronics (ISIE), Buzios, Brazil, 3–5 June 2015; Volume 6, pp. 1046–1051. [Google Scholar] [CrossRef]
  15. Miranda, H.; Cortes, P.; Yuz, J.I.; Rodriguez, J. Predictive Torque Control of Induction Machines Based on State-Space Models. IEEE Trans. Ind. Electron. 2009, 56, 1916–1924. [Google Scholar] [CrossRef]
  16. Bilal, A.; Waheed, A.; Shah, M.H. A Comparative Study of Machine Learning Algorithms for Controlling Torque of Permanent Magnet Synchronous Motors through a Closed Loop System. Proccedings of the 2019 Second International Conference on Latest Trends in Electrical Engineering and Computing Technologies (INTELLECT), Karachi, Pakistan, 13–14 November 2019; pp. 1–6. [Google Scholar] [CrossRef]
  17. Novak, M.; Xie, H.; Dragicevic, T.; Wang, F.; Rodriguez, J.; Blaabjerg, F. Optimal Cost Function Parameter Design in Predictive Torque Control (PTC) Using Artificial Neural Networks (ANN). IEEE Trans. Ind. Electron. 2021, 68, 7309–7319. [Google Scholar] [CrossRef]
  18. Xiang, M.S.Y.; Mutschler, S.; Brix, D.I.N.; Brach, C.; Geimer, M. Optimization of Hydrostatic-Mechanical Transmission Control Strategy by Means of Torque Control. Proccedings of the 12th International Fluid Power Conference, Dresden, Germany, 9–11 March 2020. [Google Scholar]
  19. Chang, S.; Gordon, T.J. Model-based predictive control of vehicle dynamics. Int. J. Veh. Auton. Syst. 2007, 5, 3. [Google Scholar] [CrossRef]
  20. Atabay, O.; Ötkür, M.; Ereke, İ.M. Model based predictive engine torque control for improved drivability. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2018, 232, 1654–1666. [Google Scholar] [CrossRef]
  21. Varatharajan, A.; Pellegrino, G.; Armando, E.; Hinkkanen, M. Sensorless Control of Synchronous Motor Drives: Accurate Torque Estimation and Control Under Parameter Errors. IEEE J. Emerg. Sel. Top. Power Electron. 2021, 9, 5367–5376. [Google Scholar] [CrossRef]
  22. Hajian, M.; Soltani, J.; Markadeh, G.; Hosseinnia, S. Adaptive Nonlinear Direct Torque Control of Sensorless IM Drives with Efficiency Optimization. IEEE Trans. Ind. Electron. 2010, 57, 975–985. [Google Scholar] [CrossRef]
  23. Kang, M.; Shen, T. Nonlinear model predictive torque control for IC engines. In Proceedings of the 11th World Congress on Intelligent Control and Automation, Shenyang, China, 29 June–4 July 2014; pp. 804–809. [Google Scholar] [CrossRef]
  24. Szabat, K. Efficient real-time model predictive control of the drive system with elastic transmission. IET Control Theory Appl. 2010, 4, 37–49. [Google Scholar]
  25. Costa, B.L.G.; Graciola, C.L.; Angélico, B.A.; Goedtel, A.; Castoldi, M.F. Metaheuristics optimization applied to PI controllers tuning of a DTC-SVM drive for three-phase induction motors. Appl. Soft Comput. 2018, 62, 776–788. [Google Scholar] [CrossRef]
  26. Elaziz, M.A.; Elsheikh, A.H.; Oliva, D.; Abualigah, L.; Lu, S.; Ewees, A.A. Advanced Metaheuristic Techniques for Mechanical Design Problems: Review. Arch. Comput. Methods Eng. 2022, 29, 695–716. [Google Scholar] [CrossRef]
  27. Oladipo, S.; Sun, Y.; Wang, Z. Optimization of PID Controller with Metaheuristic Algorithms for DC Motor Drives: Review. Int. Rev. Electr. Eng. (IREE) 2020, 15, 352. [Google Scholar] [CrossRef]
  28. Ochoa-Ortega, G.; Villafuerte-Segura, R.; Luviano-Juárez, A.; Ramírez-Neria, M.; Lozada-Castillo, N. Cascade Delayed Controller Design for a Class of Underactuated Systems. Complexity 2020, 2020, 2160743. [Google Scholar] [CrossRef]
  29. Sanin-Villa, D.; Montoya, O.D.; Grisales-Noreña, L.F. Material property characterization and parameter estimation of thermoelectric generator by using a master–slave strategy based on metaheuristics techniques. Mathematics 2023, 11, 1326. [Google Scholar] [CrossRef]
  30. Sanin-Villa, D.; Montoya, O.D.; Gil-González, W.; Grisales-Noreña, L.F.; Perea-Moreno, A.J. Parameter Estimation of a Thermoelectric Generator by Using Salps Search Algorithm. Energies 2023, 16, 4304. [Google Scholar] [CrossRef]
  31. Grisales-Noreña, L.; Cortés-Caicedo, B.; Montoya, O.D.; Sanin-Villa, D.; Gil-González, W. Integration of BESS in grid connected networks for reducing the power losses and CO2 emissions: A parallel master-stage methodology based on PDVSA and PSO. J. Energy Storage 2024, 87, 111355. [Google Scholar] [CrossRef]
  32. Mirjalili, S.; Gandomi, A.H.; Mirjalili, S.Z.; Saremi, S.; Faris, H.; Mirjalili, S.M. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. Adv. Eng. Softw. 2017, 114, 163–191. [Google Scholar] [CrossRef]
Figure 1. Schematics of the Pendubot system.
Figure 1. Schematics of the Pendubot system.
Mathematics 12 01625 g001
Figure 2. Experimental pendulum prototype.
Figure 2. Experimental pendulum prototype.
Mathematics 12 01625 g002
Figure 3. Inverted pendulum trajectory. (a) t = 0 s. (b) t = 7 s. (c) t = 13.5 s.
Figure 3. Inverted pendulum trajectory. (a) t = 0 s. (b) t = 7 s. (c) t = 13.5 s.
Mathematics 12 01625 g003
Figure 4. Experimental results for inverted pendulum.
Figure 4. Experimental results for inverted pendulum.
Mathematics 12 01625 g004
Figure 5. Comparison of average RMSE for algorithms.
Figure 5. Comparison of average RMSE for algorithms.
Mathematics 12 01625 g005
Figure 6. Comparison of average computing times for algorithms.
Figure 6. Comparison of average computing times for algorithms.
Mathematics 12 01625 g006
Figure 7. Comparison of average RMSE for algorithms.
Figure 7. Comparison of average RMSE for algorithms.
Mathematics 12 01625 g007
Table 1. Upper and lower limits for the parameters to be estimated on the model.
Table 1. Upper and lower limits for the parameters to be estimated on the model.
Parameter m 1 [ kg ] lc 1 [ m ] I 1 [ kgm 2 ] m 2 [ kg ] lc 2 [ m ] I 2 [ kgm 2 ]
Upper0.40.150.0010.40.30.001
Lower0.0010.0010.00010.0010.10.0001
Table 2. Comparative results of parameter estimation.
Table 2. Comparative results of parameter estimation.
ParametersExperimentalSSACGAPSO
m 1 0.200000.205790.149580.25497
l c 1 0.120000.150000.011940.15000
I 1 0.000530.000590.000730.00100
m 2 0.140000.134960.204610.13557
l c 2 0.150000.120010.104020.10000
I 2 0.000770.000510.000120.00010
Table 3. Comparative results of metaheuristic algorithms for parameter estimation—error metrics.
Table 3. Comparative results of metaheuristic algorithms for parameter estimation—error metrics.
AlgorithmsBest Solution ErrorMean Solutions ErrorStandard Deviation of ErrorConfidence Interval for Error
SSA0.0150151470.015066207 1.44399 × 10 6 [0.01506, 0.01506]
CGA0.0150396070.015490795 7.34505 × 10 6 [0.01548, 0.01549]
PSO0.0150450760.016549584 1.47122 × 10 5 [0.01654, 0.01655]
Table 4. Comparative results of metaheuristic algorithms for parameter estimation—time metrics.
Table 4. Comparative results of metaheuristic algorithms for parameter estimation—time metrics.
AlgorithmsMean TimeStandard Deviation of TimeConfidence Interval for Time
SSA916.561418.7228[912.89167, 920.23103]
CGA954.185118.7228[950.51546, 957.85481]
PSO888.477117.4958[885.04796, 891.90631]
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

Sanin-Villa, D.; Rodriguez-Cabal, M.A.; Grisales-Noreña, L.F.; Ramirez-Neria, M.; Tejada, J.C. A Comparative Analysis of Metaheuristic Algorithms for Enhanced Parameter Estimation on Inverted Pendulum System Dynamics. Mathematics 2024, 12, 1625. https://doi.org/10.3390/math12111625

AMA Style

Sanin-Villa D, Rodriguez-Cabal MA, Grisales-Noreña LF, Ramirez-Neria M, Tejada JC. A Comparative Analysis of Metaheuristic Algorithms for Enhanced Parameter Estimation on Inverted Pendulum System Dynamics. Mathematics. 2024; 12(11):1625. https://doi.org/10.3390/math12111625

Chicago/Turabian Style

Sanin-Villa, Daniel, Miguel Angel Rodriguez-Cabal, Luis Fernando Grisales-Noreña, Mario Ramirez-Neria, and Juan C. Tejada. 2024. "A Comparative Analysis of Metaheuristic Algorithms for Enhanced Parameter Estimation on Inverted Pendulum System Dynamics" Mathematics 12, no. 11: 1625. https://doi.org/10.3390/math12111625

APA Style

Sanin-Villa, D., Rodriguez-Cabal, M. A., Grisales-Noreña, L. F., Ramirez-Neria, M., & Tejada, J. C. (2024). A Comparative Analysis of Metaheuristic Algorithms for Enhanced Parameter Estimation on Inverted Pendulum System Dynamics. Mathematics, 12(11), 1625. https://doi.org/10.3390/math12111625

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