Parameters Identiﬁcation of a Permanent Magnet DC Motor: A Review

: Since permanent magnet direct current (DC) motors are mainly used in various industrial automation applications, the demand for electric motors is increasing rapidly. However, in the mass production of electric motors, often, only random inspections are used to check the speciﬁcations and performance of electric motors. For manufacturing or engineering application staff to have a more thorough understanding of the characteristics of the motor, it is necessary to conduct a full or quick inspection during the production process to ensure the quality of the electric motor. Based on this, this literature review reveals several methods and algorithms often used to estimate DC motor parameters, given the importance of knowing the parameters of the DC motor and the lack of research on estimating the parameters of the DC motor.


Introduction
Permanent magnet direct current (DC) motors (also known as brushed DC motors) have been widely employed in numerous industries and applications, ranging from robotics and automation to automotive systems. The global brushed DC motor market is expected to reach USD 8742.20 million by 2029 at a CAGR of 5.30% during the forecast period, as shown in Figure 1 [1]. The brushed DC motor is an excellent choice for low-torque applications because it can change pace or speed with variable speed control options, including below-and above-rated speed [2]. Examples of brushed DC motors are car robots used for exploration, surveying, or mapping purposes, such as MER-B and MER-1 for exploitation on the planet Mars; robot vacuum cleaners; and toy robots. Brushed motors are often used in industries with high torque requirements at low speeds, such as a printing press, spinning drive, or agitator. Nevertheless, age and accompanying wear have changed the characteristics of DC motors. As a result, adjusting the motor characteristics as the motor matures can improve the DC motor model's accuracy [3].

Introduction
Permanent magnet direct current (DC) motors (also known as brushed DC motors) have been widely employed in numerous industries and applications, ranging from robotics and automation to automotive systems. The global brushed DC motor market is expected to reach USD 8742.20 million by 2029 at a CAGR of 5.30% during the forecast period, as shown in Figure 1 [1]. The brushed DC motor is an excellent choice for lowtorque applications because it can change pace or speed with variable speed control options, including below-and above-rated speed [2]. Examples of brushed DC motors are car robots used for exploration, surveying, or mapping purposes, such as MER-B and MER-1 for exploitation on the planet Mars; robot vacuum cleaners; and toy robots. Brushed motors are often used in industries with high torque requirements at low speeds, such as a printing press, spinning drive, or agitator. Nevertheless, age and accompanying wear have changed the characteristics of DC motors. As a result, adjusting the motor characteristics as the motor matures can improve the DC motor model's accuracy [3].

DC Motor Dynamic Response and Parameter Estimation
The dynamic response of a DC motor is its ability to modify its speed or torque in reaction to changes in the input signal fast and precisely. The electrical and mechanical properties of a DC motor govern its dynamic reaction. The motor's dynamic reaction is influenced mainly by its electrical time constant, mechanical time constant, and inertia. The time it takes for the electrical circuit of the motor to stabilize at 63.2% of its final value following a step change in the input signal is known as the electrical time constant. The resistance and inductance of the motor's windings determine it. The time required for the mechanical system of the motor to reach 63.2% of its final value following a step Electronics 2023, 12, 2559 3 of 24 change in the output torque is known as the mechanical time constant. The rotor and load inertia of the motor and any friction or damping in the system all contribute to this. The motor's inertia also has an impact on its dynamic responsiveness. A motor with low inertia will react more rapidly, whereas a motor with high inertia will take longer to establish a steady-state speed or torque.
The advantage of the dynamic response is that it allows the motor to respond to changes in input or load quickly and with high precision. This allows the motor to change speed or torque instantly and accurately, making it essential in applications requiring rapid changes, high control, and adaptability to changing operating conditions. The dynamic response also allows for more in-depth measurement and analysis of the motor's response to changes, which is essential in maintaining, repairing, and optimizing DC motor performance. Maintaining system stability and avoiding oscillations or overshoots in the motor's output need a quick and precise dynamic reaction. Additionally, a well-tuned control system can aid in lowering energy usage and enhancing the motor's general effectiveness.
The dynamic relationships between the DC motor equations are briefly discussed in this section. The dynamic reaction of the motor may be divided into three stages: acceleration, steady state, and deceleration, as shown in Figure 2 [45]. The variables that make up the DC motor mathematical model are described in Table 1.
Electronics 2023, 12, x FOR PEER REVIEW 3 of 25 resistance and inductance of the motor's windings determine it. The time required for the mechanical system of the motor to reach 63.2% of its final value following a step change in the output torque is known as the mechanical time constant. The rotor and load inertia of the motor and any friction or damping in the system all contribute to this. The motor's inertia also has an impact on its dynamic responsiveness. A motor with low inertia will react more rapidly, whereas a motor with high inertia will take longer to establish a steady-state speed or torque. The advantage of the dynamic response is that it allows the motor to respond to changes in input or load quickly and with high precision. This allows the motor to change speed or torque instantly and accurately, making it essential in applications requiring rapid changes, high control, and adaptability to changing operating conditions. The dynamic response also allows for more in-depth measurement and analysis of the motor's response to changes, which is essential in maintaining, repairing, and optimizing DC motor performance. Maintaining system stability and avoiding oscillations or overshoots in the motor's output need a quick and precise dynamic reaction. Additionally, a well-tuned control system can aid in lowering energy usage and enhancing the motor's general effectiveness.
The dynamic relationships between the DC motor equations are briefly discussed in this section. The dynamic reaction of the motor may be divided into three stages: acceleration, steady state, and deceleration, as shown in Figure 2 [45]. The variables that make up the DC motor mathematical model are described in Table 1.    The acceleration state of a DC motor refers to its ability to change its speed in response to a change in the input signal. It specifically relates to how quickly the motor speeds up when the input signal moves from a lower value to a higher value. A DC motor's dynamic response and acceleration state are closely connected, and the motor's electrical and mechanical time constants, inertia, and load torque all play a role. As the input signal changes, a motor in a fast acceleration condition can react swiftly and quickly attain its new steady-state speed.
The steady state of a DC motor refers to the point at which the motor's speed and torque have stabilized and are no longer changing in response to changes in the input signal. A DC motor's dynamic response governs how rapidly it may regain a steady state following a change in the input signal. A motor with a slower reaction will take more time to stabilize, whereas one with a better dynamic response can quickly attain a steady-state speed or torque following a change in the input signal. Furthermore, a well-tuned control system can aid in reducing oscillations and overshoot of the motor's output during the transition to a steady state, resulting in a smoother and more stable operation.
The deceleration state of a DC motor refers to its ability to slow down and stop in response to a change in the input signal. It specifically refers to the rate at which the input signal's value shifts from higher to lower, causing the motor's speed to drop. A DC motor's dynamic reaction is directly tied to its deceleration state, determined by its electrical and mechanical time constants, inertia, and load torque. As the input signal changes, a motor with a fast deceleration condition may react swiftly and halt quickly.

DC Motor Mathematical Model of DC Motor
DC motors are hybrid systems made up of both mechanical and electrical components. Therefore, electrical and mechanical equations allow for observing the motor's dynamic behavior, as seen in Equations (1) and (2).
The two connection equations discussed in Equations (3) and (4) are connected to the previous set of equations.
By substituting Equations (3) and (4) into Equations (1) and (2) and assuming that the motor is without load (T L = 0), we obtain Equations (5) and (6), and their representation can be observed in Figure 3 [40]. According to the results of Equations (5) and (6), Figure 4 displays the motor model's block diagram.

Least Squares Method
The least squares method was developed in 1795 by Karl Friedrich Gauss, who also used it for astronomical calculations. He proposed that the most likely values are also the most suitable for the unknowable but desired parameters. According to his definition, "The value of the unknown quantities that has the greatest chance of being correct is the one for which the accuracy-measuring numbers are multiplied by the sum of the squares of the discrepancies between the observed and calculated values". Since then, Electronics 2023, 12, 2559 5 of 24 several technical issues have been solved using the least squares technique. Numerous analyses of its qualities were conducted, numerical methods for various applications were suggested, and the technique was changed to meet the required specifications. The same attempt can be seen in the controlled system parameter estimation area. In recent years, the issue of recognizing a dynamic process has drawn much attention. The numerous methods created for collecting and evaluating input-output data vary from the most basic deterministic procedures to beautiful numerical and statistical strategies based on the findings of optimum estimation theory. As a result, the least squares approach has significantly improved and grown in popularity in parameter estimation.

Least Squares Method
The least squares method was developed in 1795 by Karl Friedrich Gauss, who also used it for astronomical calculations. He proposed that the most likely values are also the most suitable for the unknowable but desired parameters. According to his definition, "The value of the unknown quantities that has the greatest chance of being correct is the one for which the accuracy-measuring numbers are multiplied by the sum of the squares of the discrepancies between the observed and calculated values". Since then, several technical issues have been solved using the least squares technique. Numerous analyses of its qualities were conducted, numerical methods for various applications were suggested, and the technique was changed to meet the required specifications. The same attempt can be seen in the controlled system parameter estimation area. In recent years, the issue of recognizing a dynamic process has drawn much attention. The numerous methods created for collecting and evaluating input-output data vary from the most basic deterministic procedures to beautiful numerical and statistical strategies based on the findings of optimum estimation theory. As a result, the least squares approach has significantly improved and grown in popularity in parameter estimation.
Ref. [11] used the Simulink design optimization method to estimate unknown PMDC motor parameter values using pattern search and nonlinear least squares. The PMDC motor dynamic equation was used to create a PMDC motor Simulink model and calculate the initial estimate of the motor parameters, which were used as initial values in the Sim-

Least Squares Method
The least squares method was developed in 1795 by Karl Friedrich Gauss, who also used it for astronomical calculations. He proposed that the most likely values are also the most suitable for the unknowable but desired parameters. According to his definition, "The value of the unknown quantities that has the greatest chance of being correct is the one for which the accuracy-measuring numbers are multiplied by the sum of the squares of the discrepancies between the observed and calculated values". Since then, several technical issues have been solved using the least squares technique. Numerous analyses of its qualities were conducted, numerical methods for various applications were suggested, and the technique was changed to meet the required specifications. The same attempt can be seen in the controlled system parameter estimation area. In recent years, the issue of recognizing a dynamic process has drawn much attention. The numerous methods created for collecting and evaluating input-output data vary from the most basic deterministic procedures to beautiful numerical and statistical strategies based on the findings of optimum estimation theory. As a result, the least squares approach has significantly improved and grown in popularity in parameter estimation.
Ref. [11] used the Simulink design optimization method to estimate unknown PMDC motor parameter values using pattern search and nonlinear least squares. The PMDC motor dynamic equation was used to create a PMDC motor Simulink model and calculate the initial estimate of the motor parameters, which were used as initial values in the Sim- Ref. [11] used the Simulink design optimization method to estimate unknown PMDC motor parameter values using pattern search and nonlinear least squares. The PMDC motor dynamic equation was used to create a PMDC motor Simulink model and calculate the initial estimate of the motor parameters, which were used as initial values in the Simulink design optimization method.
Using Equations (1)-(4), the relationship between K e and K t can be determined. The back emf constant of the motor and the torque constant of the motor will be described in Equations (7) and (8).
By assuming that the electromagnetic losses are equal to zero, the connection between these constants may be found. This shows that the mechanical power and the electrical power that the back emf voltage dissipates in the armature are equivalent. Equation (9) defines the relationship between mechanical power (P m ) and electrical power (P e ): From Equations (7) and (8) above, we obtain Equation (10). The result implies K e = K t .
The motor parameter values that result from the use of the nonlinear least square algorithm are considered accurate results because the output response from the Simulink modal is almost symmetrical with the actual measured response of the motor, which makes the results obtained using the nonlinear least squares algorithm more accurate, as described in Table 2. Ref. [12] describes estimating dc motor parameters in steady state and deceleration. It is explained that the value of K e = K t , so there are four unknown parameters in the steadystate interval Armature resistance (R), the mechanical constants equal values (K t ), torque loss (T 0 ), and Friction coefficient (B). In this article, two sets of voltage inputs, V ss1 and V ss2 , are used as measurements to improve the accuracy of parameter estimation. Any unknown parameters can be solved using the least squares method described in Equation (11).
In the deceleration interval, the current equals zero, and the motor slowly stops. At this time, the motor inertia calculation formula can be explained by the Equation (12).
This article uses the open source Python program to develop a parameter measurement system for permanent magnet DC motors. Among them, the signal measurement system built on the Python platform captures the motor drive's voltage, current, and speed signals and integrates the signal data. Enter MATLAB R2021a to estimate the acceleration parameters and the least squares method and propose a method to judge the acceleration and deceleration intervals. The results of the acceleration method are compared with the measurement using the least squares method. The result is that using the least squares method has more stable and complete results for each DC motor parameter described in Table 3.

Metaheuristics Algorithms
Over the past few years, there has been a surge in the utilization of metaheuristic algorithms. The primary reason for this is the ease of implementation and adaptability to diverse problems as expounded in reference [24]. Hence, metaheuristic algorithms based on the population are frequently employed in the academic literature. The categorization of metaheuristic algorithms can be delineated into four principal groups based on the inspiratory source utilized to generate novel solutions. These four categories are evolutionary algorithms (EA), swarm intelligence algorithms, physics-based algorithms, and human-based algorithms, as detailed in reference [46]. It is commonly posited within the scholarly literature that the impetus behind the advancement of Evolutionary Algorithms (EA) can be attributed to the influential framework provided by Darwin's seminal theory of evolution, which emphasizes the notion of the preservation of the "fittest" individuals in the maintenance of species survival. This study employs genetic procedures, specifically, crossover, mutation, and selection, in order to generate improved quality offspring solutions. The swarm intelligence algorithm is rooted in the emulation of the aggregate actions exhibited by animals or insects, more specifically, in their pursuit of locating food sources and securing potential partners for the propagation of their species. Intelligent behavior, characterized by decentralized decision making and informed by the examination of local information and interaction with the search environment, can be exhibited by members of this flock. The tracing operators of algorithms that are based on human activity have the ability to imitate various human behaviors such as thinking, learning, speaking, and teaching. A taxonomy of metaheuristic algorithms based on their source of inspiration and notable examples associated with each branch of MSA is presented in Figure 5. In the present study, a review was conducted with a focus on three specific population-based metaheuristic algorithms, namely, differential evolution, particle swarm optimization, and cuckoo search. The selection of these algorithms was based on their pertinence to examining the implications of utilizing dynamic connections within the context of population-based metaheuristic algorithms.

Differential Evolution (DE)
DE is a population-based direct stochastic optimization technique initially developed by Storn and Price [70]. It is easy, powerful, and simple, qualities that make it appealing for numerical optimization over continuous search space. Compared to evolutionary algorithms, DE takes a more greedy and less stochastic approach to problem resolution. DE combines straightforward arithmetic operators with the traditional crossover, mutation, and selection operators to evolve from a randomly generated beginning population to a final solution. In addition, the way that DE implements the mutation process differs significantly from other evolutionary algorithms. The vector differentials between the members of the current population are used in DE mutation operation to determine the level and direction of perturbation that will be applied to each subject of the mutation operation.

Differential Evolution (DE)
DE is a population-based direct stochastic optimization technique initially developed by Storn and Price [70]. It is easy, powerful, and simple, qualities that make it appealing for numerical optimization over continuous search space. Compared to evolutionary algorithms, DE takes a more greedy and less stochastic approach to problem resolution. DE combines straightforward arithmetic operators with the traditional crossover, mutation, and selection operators to evolve from a randomly generated beginning population to a final solution. In addition, the way that DE implements the mutation process differs significantly from other evolutionary algorithms. The vector differentials between the members of the current population are used in DE mutation operation to determine the level and direction of perturbation that will be applied to each subject of the mutation operation. Scholarly article [18] enumerates several methodologies employed in optimization, particularly the Genetic Algorithm (GA), Differential Evolutions (DE) with two strategies (DE/rand/1/exp, DE/best/1/bin), Teaching-Learning-Based Optimization (TLBO), and Artificial Bee Colony (ABC). This article presents an analysis of the variation in simulation time of the motor model and memory assistance, utilizing three distinct approaches to minimize the overall computational time. The diagrammatic representation of the DC motor is depicted in Figure 6, along with the operative device that may or may not have a correlation with the motor. The potential association between the item in question and the motor system is uncertain and warrants further investigation. The specification of the drive's parameter requirements is contingent upon various factors. The proposed methodology can be applied to DC motors featuring distinct excitation systems such as parallel excitation or permanent magnet configurations. minimize the overall computational time. The diagrammatic representation of the DC motor is depicted in Figure 6, along with the operative device that may or may not have a correlation with the motor. The potential association between the item in question and the motor system is uncertain and warrants further investigation. The specification of the drive's parameter requirements is contingent upon various factors. The proposed methodology can be applied to DC motors featuring distinct excitation systems such as parallel excitation or permanent magnet configurations. In this article, the motor stroke response can be simulated using Equations (1) and (2). Equation (2) is converted into Equation (13) which is described below: Equations (1) and (2) are combined because the induced voltage relies on the motor speed and the motor torque ( ) depends on the motor current. In Equations (14) and (15), and ( ) are written.
where is the DC motor constant and the magnetic flux is considered to be constant in the instance of the particular motor, is taken to be constant (referred to as the motor constant in the continuation). The load given in (2) can be broken down further, as shown in (16).
The components and , respectively, stand in for Coulomb and viscous friction when there is no load from the working machine, while stands in for any air resistance in the fan along the motor axis. In the event that the working machine creates a load, the variables , , and indicate friction, air resistance at the fan at the motor axis, and the combined working machine load. Equations (17) and (18) are created by putting the following Equations (14)-(16) into (1) and (2). In this article, the motor stroke response can be simulated using Equations (1) and (2). Equation (2) is converted into Equation (13) which is described below: Equations (1) and (2) are combined because the induced voltage e relies on the motor speed and the motor torque τ(t) depends on the motor current. In Equations (14) and (15), e and τ(t) are written. e = c m ω (14) Since c m = k m Φ, where k m is the DC motor constant and the magnetic flux Φ is considered to be constant in the instance of the particular motor, c m is taken to be constant (referred to as the motor constant in the continuation). The load T L given in (2) can be broken down further, as shown in (16).
The components T la and T lb , respectively, stand in for Coulomb and viscous friction when there is no load from the working machine, while T lc stands in for any air resistance in the fan along the motor axis. In the event that the working machine creates a load, the variables T la , T lb , and T lc indicate friction, air resistance at the fan at the motor axis, and the combined working machine load. Equations (17) and (18) are created by putting the following Equations (14)- (16) into (1) and (2).
Results show that Differential Evolution (DE)/rand/1/exp is the most powerful method described in Table 4. The results are the same for each run within the scope of the individual test data. TLBO is only slightly worse, but computation time is two times longer. The results obtained using DE/best/1/bin and ABC are also acceptable. Only the results obtained using GA are too bad. Based on that, DE/rand/1/exp can be the best method for the problem discussed. The weakness of this article is that it does not consider the dynamic properties of torsion, mass with torsional stiffness, and damping effects. The behavior of a swarm of living objects, such as birds and other flies, is mimicked by the bio-inspired metaheuristic search approach known as particle swarm optimization [71]. The standard particle swarm optimization is known to attract suboptimal solutions that are only sometimes the best solution [71,72]. As a result, the search process quickly converges. The particles continue searching within a confined and smaller area of the search space after they are caught at such local optimal solutions. All particles seek locally in certain areas, making it difficult or impossible to escape from such local optimum search ranges. Once the particles have prematurely converged, they may continue to do so until they reach a solution that is very close to the others. As a result, the personal and overall best may be limited to that tiny area of the search field. As a result, the algorithm will not be able to execute the exploration process, which will improve the exploitation process. A typical PSO randomly starts a swarm of particles in the solution search space with a normal distribution and specified upper and lower limit ranges. Each particle has a velocity vector that directs its subsequent travel while being measured against a predetermined goal function to obtain the optimal answer [71,73]. Each particle's velocity is adjusted following its personal best solution (pBest) and overall best solution (gBest). A particle finds the best solution is pBest, while the entire swarm finds the best solution is gBest [72,73].
A particle may quickly become stuck in a local optimum in the typical PSO, which might converge to a decent solution but has the issue of premature convergence. This is especially true in a significant problem search space with numerous local optimal solutions. This limits the particles' ability to discover new regions of the search space. The introduction and manipulation inertia weight factor, the inclusion of diversity in the swarm, and the hybridization of an algorithm with another metaheuristic algorithm, such as genetic algorithms, are just a few of the solutions that have been suggested in the literature [74]. These techniques enable the particles to avoid being trapped in a local optimum and improve the global searching process.
Ref. [26] compared the modified PSO method to the original PSO algorithm. PSO adaptations include Chaos-Initialized PSO (CI PSO), Adaptive Inertia Weight Factor PSO (AIWF PSO), and constricted with linearly decreasing inertia weight PSO (CIW PSO). Computer modeling and system parameter estimation were performed using MATLAB/Simulink software. This article presents a DC motor system with the two-mass model (2MM) described in the Figure 6 and will be simplified using Equations (19)- (21).
In Figure 7, the motor is linked to the load inertia through the implementation of a flexible shaft. The variable J 1 signifies the inertia of the motor, while J 2 denotes the inertia of the load. Flexible couplings are frequently employed in order to facilitate the connection of loads, thus permitting a particular degree of shaft misalignment. This particular mechanism is observable in various systems such as wind turbines and robot arms that employ flexible joints, assuming that the flexible coupling exhibits greater stiffness relative to the flexible center shaft. In this scenario, a simplified model of the system may be employed, in which solely two masses are interlinked by a pliant shaft, with the consideration of their inherent damped effect. Equation (22) explains the transfer function of the system, which takes motor torque as input and produces the inertial angular velocity of the motor as an output.
= − + ( − ) − (20) In Figure 7, the motor is linked to the load inertia through the implementation of a flexible shaft. The variable signifies the inertia of the motor, while denotes the inertia of the load. Flexible couplings are frequently employed in order to facilitate the connection of loads, thus permitting a particular degree of shaft misalignment. This particular mechanism is observable in various systems such as wind turbines and robot arms that employ flexible joints, assuming that the flexible coupling exhibits greater stiffness relative to the flexible center shaft. In this scenario, a simplified model of the system may be employed, in which solely two masses are interlinked by a pliant shaft, with the consideration of their inherent damped effect. Equation (22) explains the transfer function of the system, which takes motor torque as input and produces the inertial angular velocity of the motor as an output. The Equation (23) represents the task of identifying the system model parameters for the 2MM as an optimization problem utilizing a minimal cost function.
The result is that the Constricted with Linearly Decreasing PSO (CIW PSO) outperforms all algorithms with its ability to get out of the local optimal, converge to values very close to the real ones, and obtain the best cost function, as described in Table 5. The Equation (23) represents the task of identifying the system model parameters for the 2MM as an optimization problem utilizing a minimal cost function.
The result is that the Constricted with Linearly Decreasing PSO (CIW PSO) outperforms all algorithms with its ability to get out of the local optimal, converge to values very close to the real ones, and obtain the best cost function, as described in Table 5. In reference to Ref. [27], a novel approach for determining the parameters of independently excited DC motors is presented utilizing the Chaotic Initialized Particle Swarm Optimization (CIPSO) algorithm. The process of parameter estimation is transformed into an optimization problem through the utilization of an objective function. The strategy that has been presented holds significant importance in the precise estimation of motor parameters compared to the standard particle swarm optimization technique. This is evidenced by the low mean squared error observed between the actual speeds and the estimated speeds. The DC motor is represented by means of a transfer function in the context of academic writing. In this study, the optimal estimation of five important parameters of a DC motor, namely, the moment of inertia, viscous friction, electromotive force constant, resistance, and inductance, was performed through the application of CIPSO. The generation of the initial population swarms was accomplished through the implementation of a stochastically generated marquee map. The estimated parameters are subjected to a comparison with the true values alongside the parameters obtained through the standard particle swarm optimization technique. The CIPSO method offers the distinct advantages of providing relatively precise parameter estimates with an insignificant mean squared error. In this article, the findings presented in Table 6 indicate the superior performance of the CIPSO algorithm, compared to the SPSO algorithm, in relation to both mean square error (MSE) values and the accuracy of the estimations for the motor parameters. The MSE value attained by CIPSO, which is much lower than the MSE achieved by SPSO, is 1.399 × 10 −16 . The CIPSO-estimated parameters and the real motor parameters are quite near to each other.

Cuckoo Search Optimization (CSO)
The CSO algorithm exhibits a comparatively complex nature, unlike the remaining two algorithms. In contrast to the aforementioned algorithms, similarities to the Genetic Algorithm (GA) are evident in this approach, as with the previous one. Notably, however, convergence is demonstrated in a reduced number of iterations. The CSO method constitutes a bioinspired algorithm that draws inspiration from the peculiar reproductive strategy of cuckoo birds, which involves laying eggs in the nests of other bird species for subsequent care and rearing. In contrast to the remaining two algorithms, the CSO method comprises specific parameters. The determination of certain parameters is a complex subject of inquiry in academia. The body of scholarly works have recommended a proportion of 25% [75].
Ref. [39] introduced a modified cuckoo search algorithm as an effective parametric estimation tool for optimizing motor performance. The charge function is derived from the discrepancy between the actual and desired current and velocity measurements, which arise upon applying a step input voltage to the motor. The present study entailed a comparison of the modified metaheuristic cuckoo search algorithm with the Steiglitz-McBride method and the standard cuckoo search technique.
By utilizing Equations (1)-(4), it is possible to formulate a differential equations system predicated upon the motor's current and speed. In this manner, a load-free motor configuration (T L = 0) is considered, and substitutions are made by utilizing (3) and (4) on (1) and (2), correspondingly. Equations (24) and (25) represent the consequent outcome of this procedure.
Equation (24) represents the current equation, whilst Equation (25) pertains to the velocity equation. The two variables under consideration, namely, velocity and current, are deemed to be quite feasible in terms of measurement. Nonetheless, the values of the six undetermined variables are contingent upon the type of motor utilized. Conversely, it may be posited that electrical constants exhibit magnitudes that are typically commensurate with those of mechanical constants [2]. Consequently, the issue may be simplified to acquiring five key parameters. This article exclusively employs the steady-state equation for determining the parameters. A steady state refers to a condition in which the variables characterizing the system remain constant over time. In regard to regulation, the state is defined as stabilized when the magnitude of the output remains within 2% of the ultimate value of the response to the step input. Thus, in order to derive distinct correlations, an analysis of the system's response to the step input should be conducted. By solely considering the steady state, Equations (7) and (8) can be utilized to derive the subsequent association.
In the present context, the parameter denoted by K serves as a measure of the magnitude of both the electrical and mechanical constants. I ss denotes the steady-state current, while V ss represents the voltage exerted in steady-state conditions, also known as the step magnitude. Additionally, ω ss signifies the velocity attained during steady-state conditions. It is contended that Equations (26) and (27) are applicable exclusively to the motor reaction to step inputs in a state of equilibrium.
Numerous scholarly investigations have demonstrated that the Integral Absolute Error (IAE) along the trajectory is a fitness function of notable efficacy. The motor is explicated through a set of differential equations. Hence, it is crucial to consider the inaccuracies present in both the current and velocity variables. Prior studies have demonstrated the suitability of employing the Euclidean distance as a form of fitness function when dealing with two vectors. Specifically, the original proposed fitness function, denoted as (28), involves the use of estimated values for I ss and ω ss .
The outcomes demonstrate a noteworthy enhancement over the conventional cuckoo search algorithm and exhibit superior performance compared to the Steiglitz-McBride technique. The efficacy of the modified cuckoo search algorithm in comparison with the standard cuckoo search algorithm has been demonstrated through results tabulated in Table 7. Notably, the modified cuckoo search algorithm adheres to the steady-state relationship, which is frequently invalidated by the application of random values estimated by the standard cuckoo search algorithm. The tabulated data for the performance of the two algorithms is presented in Table 7. The root means square error (RMSE) for velocity and current is computed for each methodology. The algorithm developed by Steiglitz and McBride demonstrates an RMSE value of 2.1542 for velocity and an RMSE value of 0.0296. The modified cuckoo search method exhibits an RMSE of 0.8562 for velocity and an RMSE of 0.0242 Table 7. Comparing nominal parameters, parameters acquired using the original cuckoo method, and parameters obtained using the modified cuckoo algorithm.

Parameter
Nominal In contrast, it offers marginally superior outcomes than the Steiglitz-McBride algorithm. The modified cuckoo search's shortcoming, despite the fact that this approach is an enhanced heuristic and metaheuristic method, is that it is only possible to use the modified cuckoo search with DC motors.
Ref. [40] employed a test that illustrates the use of dynamic response relationships as search constraints in a metaheuristic algorithm employed as a parametric estimator, and it was conducted with three algorithms (Gray Wolf Optimizer, Jaya Algorithm, and cuckoo search algorithm). This article compares the three original algorithms with three modified algorithms. Not only are they compared with the original algorithms, but the calculations in this article also use two motors for comparison.
The steady state is employed as a search constraint in article [40], which is similar to article [39] in that it estimates the DC motor parameters solely in the transient state. As a result of this connection, the search algorithm's random parameters are used less frequently, and it requires fewer iterations to provide usable results.
The system is considered to be in a steady state when its output has stabilized and its reaction varies over time by less than 2%. The current and speed remain consistent in this regard. As a result, both the velocity derivative and the current derivative are zero. The similarity in magnitude between the electrical and mechanical constants is a well-known phenomenon. K is therefore utilized exactly in both cases.
When the system's output varies with respect to time, this is known as a transient. Even though it is known that the current has a direct proportionate impact on the velocity derivative, the relationship in this phase is not readily visible. This relationship causes the maximum current value to cause the current derivative's maximum value, as explained in Equation (29).
Since B is already defined as a function of K, removing J only leaves it specified as a function of K, as can be shown from Equation (29); the parameter J can be a function of both B and K.
In this article, the fitness/cost function is the Euclidean distance between the sum of the errors in flow and the sum of the errors in speed, explained in Equation (31).
The absolute percentage error for each test concerning the nominal parameters was calculated to assess the algorithm performance's correctness. Tables 8 and 9 for Motor 1 and Motor 2 describe the findings. The numerical and tabular results demonstrate improvements when comparing the modified algorithms to the originals. The focus of the findings that follow will be the changed algorithms. The observed findings show that the updated method converges more quickly and adjusts the cost curves in the initial rounds. Given that the adjustment speed had been preset, the search ended when the cost function reached a value of less than 0.1. The drawback of this article is that it was only tested in an accelerated state and a steady state, even though measuring a DC motor according to changes in speed has three stages: acceleration, steady state, and deceleration.
The beginning values are those that each algorithm uses to determine the outcomes and the improvement. An exact improvement cannot be measured because the method starts with random numbers. Nonetheless, the findings demonstrate a clear improvement in each algorithm, particularly in the first iterations, where the application of dynamic response relations causes the error to be decreased in fewer iterations. As a result, all algorithms show decreased error and quicker convergence. The pattern persists across different motors, although the GWO m performs better with Motor 1 than with Motor 2.

Discussion
This section discusses in detail the comparison of each method and algorithm for estimating DC motor parameters and comparing the results of each study.

Least Squares Method
The six measurements above can be divided into two parts: measurements using simulations in MATLAB Simulink and direct measurements using real DC motors, which are described in Table 10. 2. Improvement and Implementation of Model Identification for Permanent Magnet DC Motors [12] Ref. [11] has the drawback of not having comparative data between the value manufacturer and parameter estimates, which are fully explained in Ref. [12]. With this value comparison, the reader can find out how effectively this method is used to solve the problem of estimating DC motor parameters. Article [12] also explains in detail how to obtain DC motor parameter estimates.

Differential Evolution
Article [18] examines measurements for estimating DC motor parameters using differential evolution, using MATLAB Simulink to create DC motor simulations. Both have the same advantages by comparing parameter estimation results with several algorithms. The advantages of differential evolution in Ref. [18] obtain the best results but have a weakness. Namely, the calculation time is twice as long as the other algorithms. Another drawback in Ref. [18] is that it does not consider dynamic torsional properties, mass with torsional stiffness, or damping effects, which could be the topic of the following journal. The strength of this article is that it considers the constant, linear, and quadratic load sections, although it does not cover all the physical phenomena during motor starting. Future research will likely use different evolution to estimate DC motor parameters due to the many articles cited each year, described in Figure 8 [76]. The application of this algorithm hybridization is also widely used, among others, in Table 11. Ref. [11] has the drawback of not having comparative data between the value manufacturer and parameter estimates, which are fully explained in Ref. [12]. With this value comparison, the reader can find out how effectively this method is used to solve the problem of estimating DC motor parameters. Article [12] also explains in detail how to obtain DC motor parameter estimates.

Differential Evolution
Article [18] examines measurements for estimating DC motor parameters using differential evolution, using MATLAB Simulink to create DC motor simulations. Both have the same advantages by comparing parameter estimation results with several algorithms. The advantages of differential evolution in Ref. [18] obtain the best results but have a weakness. Namely, the calculation time is twice as long as the other algorithms. Another drawback in Ref. [18] is that it does not consider dynamic torsional properties, mass with torsional stiffness, or damping effects, which could be the topic of the following journal. The strength of this article is that it considers the constant, linear, and quadratic load sections, although it does not cover all the physical phenomena during motor starting. Future research will likely use different evolution to estimate DC motor parameters due to the many articles cited each year, described in Figure 8 [76]. The application of this algorithm hybridization is also widely used, among others, in Table 11.

Particle Swarm Optimization
The two articles [26,27] regarding measurements for estimating DC motor parameters using particle swarm optimization use MATLAB Simulink to identify DC motor parameters. The two articles explain the same flow and how to obtain the estimated DC motor parameters, depicted in Figure 9.
Electronics 2023, 12, x FOR PEER REVIEW 1 Figure 9. The process of estimating model system parameters with particle swarm optimizat Even though it uses the same flow, article [27] obtains slightly better results th ticle [26]. This difference can be due to the two articles using different DC motors. In [26], the authors modified particle swarm optimization (PSO), whereas, in an articl the authors used chaotic initialized particle swarm optimization (CIPSO) describe finite nonlinear system with a deterministic dynamic behavior that has ergodic an chastic properties [99].
Considering that there is little research on the estimation of DC motor param using particle swarm optimization here, the authors suggest that in the future, w combine chaotic initialized particle swarm optimization (CIPSO) with constricted p swarm optimization. Because chaotic initialized particle swarm optimization can im the global exploration process for the parameter search space by allowing the part fluctuate above the search space, using inertia weight methods will improve the ex tion process.

Cuckoo Search Optimization
From a review of two articles [39,40], both use Simulink for DC motor simu The two articles have their respective advantages in terms of results. The advant article [39] is that there is a more detailed measurement of the DC motor parameter mation results and it compares the RMSE of the current and velocity between mo tions of CSO, standard CSO, and Steiglitz-McBride. The strength of article [40] is the ification of three algorithms (Gray Wolf Optimizer, Jaya Algorithm, and cuckoo algorithm) to prove that the modifications made to the cuckoo search algorithm in Figure 9. The process of estimating model system parameters with particle swarm optimization.
Even though it uses the same flow, article [27] obtains slightly better results than article [26]. This difference can be due to the two articles using different DC motors. In article [26], the authors modified particle swarm optimization (PSO), whereas, in an article [27], the authors used chaotic initialized particle swarm optimization (CIPSO) described as a finite nonlinear system with a deterministic dynamic behavior that has ergodic and stochastic properties [99].
Considering that there is little research on the estimation of DC motor parameters using particle swarm optimization here, the authors suggest that in the future, we can combine chaotic initialized particle swarm optimization (CIPSO) with constricted particle swarm optimization. Because chaotic initialized particle swarm optimization can improve the global exploration process for the parameter search space by allowing the particle to fluctuate above the search space, using inertia weight methods will improve the exploration process.

Cuckoo Search Optimization
From a review of two articles [39,40], both use Simulink for DC motor simulation. The two articles have their respective advantages in terms of results. The advantage of article [39] is that there is a more detailed measurement of the DC motor parameters estimation results and it compares the RMSE of the current and velocity between modifications of CSO, standard CSO, and Steiglitz-McBride. The strength of article [40] is the modification of three algorithms (Gray Wolf Optimizer, Jaya Algorithm, and cuckoo search algorithm) to prove that the modifications made to the cuckoo search algorithm in article [39] can be extended to various metaheuristics. The modified algorithm achieved better results than the three standard algorithms above to obtain estimates of DC motor parameters. The drawback of the two articles above is that the modified algorithm can only determine the estimated DC motor parameters and cannot determine other types of motors. Measurements are only carried out in the steadystate position.
In the future, it is hoped that many will apply the cuckoo search optimization algorithm considering that there are still many variants of cuckoo search optimization that have not been used to estimate the DC motor parameters depicted in Figure 10 [100], and we can combine cuckoo search optimization with other algorithms such as hybrid cuckoo search optimization with a genetic algorithm.

Quantitative Comparison of Computational Costs
The computational cost of parameter estimation in DC motors is affected by several factors, including the dimension of the parameter space to be estimated, as a more significant number of parameters increases the search space and generally leads to an increase in computational cost for the optimization algorithm; the evaluation time required to calculate the fitness or objective function, which involves simulating motor behavior and can be time consuming; and the convergence speed of the optimization algorithm, as faster convergence reduces the number of iterations required and, consequently, lowers the overall computational cost.

Least Squares Method
The computational cost of the least squares method depends on the optimization algorithm used, such as Gauss-Newton or Levenberg-Marquardt. The least squares method is often efficient for parameter estimation in smaller-sized problems. The computational cost is relatively low compared to the metaheuristic optimization algorithms.

Differential Evolution (DE)
DE generally has moderate to high computational costs. It utilizes a complex search strategy involving trial vector generation and crossover operations. DE may require more function evaluations and iterations to converge, especially for higher-dimensional parameter spaces. The computational cost increases with the dimension of the problem and the number of iterations required.

Particle Swarm Optimization (PSO)
PSO usually has a moderate computational cost. It relies on the movement of particles through the search space to find optimal parameter values. PSO can converge relatively quickly, but the computational cost increases with larger swarm sizes and longer simulation times to evaluate the fitness function. The number of iterations required affects the computational cost.

Cuckoo Search Optimization (CSO)
CS often has a lower computational cost compared to DE and PSO. It uses a random search strategy which can efficiently explore the parameter space. However, CS may require more iterations to converge than PSO and DE, which can impact the overall computational cost. The evaluation time of the fitness function remains a significant factor. Electronics 2023, 12, x FOR PEER REVIEW 19 of 25

Quantitative Comparison of Computational Costs
The computational cost of parameter estimation in DC motors is affected by several factors, including the dimension of the parameter space to be estimated, as a more significant number of parameters increases the search space and generally leads to an increase in computational cost for the optimization algorithm; the evaluation time required to calculate the fitness or objective function, which involves simulating motor behavior and can be time consuming; and the convergence speed of the optimization algorithm, as faster convergence reduces the number of iterations required and, consequently, lowers the overall computational cost.

Least Squares Method
The computational cost of the least squares method depends on the optimization algorithm used, such as Gauss-Newton or Levenberg-Marquardt. The least squares method is often efficient for parameter estimation in smaller-sized problems. The computational cost is relatively low compared to the metaheuristic optimization algorithms.

Conclusions
This paper reviews several methods and algorithms used to estimate the excess parameters of DC motors from the least squares method. This method has more applications for estimating DC motor parameters than metaheuristic algorithms due to the ease of calculating using the least squares method. Judging from the two articles that use this method, article [12] has advantages compared to similar articles using the least squares method because article [12] uses actual motor measurements and reasonably explains how to obtain each parameter estimate.
The advantage of differential evolution is that there are many studies on this algorithm, as seen from the fact that its use continues to increase yearly. There are also many studies on this algorithm in combination with several other algorithms to achieve optimal results.
The advantages of particle swarm optimization are the speed in estimating compared to other algorithms and the ease of its relation with low limitations on the environment and the objective function. There is still enough space for research on estimating DC motor parameters.
The advantage of using cuckoo search is its simplicity, and in the future, more variants should be used, as described in Figure 10, to estimate DC motor parameters.
The advantage of using complex metaheuristic methods is that they are necessary for handling large search spaces, global optimization, non-differentiable or discrete functions, flexibility, limited problem-specific information, and balancing exploration-exploitation trade-offs in complex optimization problems, surpassing the limitations of simple methods. In the future, this metaheuristic algorithm will be even more used and developed, considering that metaheuristic algorithms can be hybridized with various methods and algorithms.
The three metaheuristic algorithms above also show good results because they are close to the parameter values of the DC motor manufacturers. Hybridizing these three metaheuristic algorithms to estimate DC motor parameters can optimize the results and existing iterations, considering that each algorithm has advantages and disadvantages.