You are accessing a machine-readable page. In order to be human-readable, please install an RSS reader.
Continue Cancel clearAll articles published by MDPI are made immediately available worldwide under an open access license. No special permission is required to reuse all or part of the article published by MDPI, including figures and tables. For articles published under an open access Creative Common CC BY license, any part of the article may be reused without permission provided that the original article is clearly cited. For more information, please refer to https://www.mdpi.com/openaccess.
Feature papers represent the most advanced research with significant potential for high impact in the field. A Feature Paper should be a substantial original Article that involves several techniques or approaches, provides an outlook for future research directions and describes possible research applications.
Feature papers are submitted upon individual invitation or recommendation by the scientific editors and must receive positive feedback from the reviewers.
Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world. Editors select a small number of articles recently published in the journal that they believe will be particularly interesting to readers, or important in the respective research area. The aim is to provide a snapshot of some of the most exciting work published in the various research areas of the journal.
Original Submission Date Received: .
Journals Active Journals Find a Journal Journal Proposal Proceedings Series Topics Information For Authors For Reviewers For Editors For Librarians For Publishers For Societies For Conference Organizers Open Access Policy Institutional Open Access Program Special Issues Guidelines Editorial Process Research and Publication Ethics Article Processing Charges Awards Testimonials Author Services Initiatives Sciforum MDPI Books Preprints.org Scilit SciProfiles Encyclopedia JAMS Proceedings Series About Overview Contact Careers News Press Blog Sign In / Sign Up Submit Journals Mathematics Volume 11 Issue 19 10.3390/math11194040 Submit to this Journal Review for this Journal Propose a Special Issue ► ▼ Article Menu Article Menu Academic Editor Junyong Zhai Subscribe SciFeed Recommended Articles Related Info Link Google Scholar More by Authors Links on DOAJ Almobaied, M. Al-Nahhal, H. S. Arrieta, O. Vilanova, R. on Google Scholar Almobaied, M. Al-Nahhal, H. S. Arrieta, O. Vilanova, R. on PubMed Almobaied, M. Al-Nahhal, H. S. Arrieta, O. Vilanova, R. /ajax/scifeed/subscribe Article Views Citations - Table of Contents Altmetric share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up ... Endorse textsms ... Comment Need Help? SupportFind support for a specific problem in the support section of our website.
Get Support FeedbackPlease let us know what you think of our products and services.
Give Feedback InformationVisit our dedicated information section to learn more about MDPI.
Get Information clear JSmol Viewer clear first_page Download PDF settings Order Article Reprints Font Type: Arial Georgia Verdana Font Size: Aa Aa Aa Line Spacing: Column Width: Background: Open AccessArticle Design a Robust Proportional-Derivative Gain-Scheduling Control for a Magnetic Levitation System by Moayed AlmobaiedMoayed Almobaied SciProfiles Scilit Preprints.org Google Scholar 1,*, Hassan S. Al-NahhalHassan S. Al-Nahhal SciProfiles Scilit Preprints.org Google Scholar 1, Orlando ArrietaOrlando Arrieta SciProfiles Scilit Preprints.org Google Scholar 2,3 and Ramon VilanovaRamon Vilanova SciProfiles Scilit Preprints.org Google Scholar 3,* 1 Electrical Engineering Department, Islamic University of Gaza, Gaza 108, Palestine 2 Instituto de Investigaciones en Ingeniería, Facultad de Ingeniería, Universidad de Costa Rica, San Jose 11501-2060, Costa Rica 3 Departament de Telecomunicació i d’Enginyeria de Sistemes, Escola d’Enginyeria, Universitat Autònoma de Barcelona, Bellaterra, 08193 Barcelona, Spain * Authors to whom correspondence should be addressed. Mathematics 2023, 11(19), 4040; https://doi.org/10.3390/math11194040 Submission received: 8 September 2023 / Revised: 19 September 2023 / Accepted: 20 September 2023 / Published: 23 September 2023 (This article belongs to the Section Engineering Mathematics) Download keyboard_arrow_down Download PDF Download PDF with Cover Download XML Download Epub Browse Figures Versions NotesAbstract: This study focuses on the design of a robust PD gain-scheduling controller (PD-GS-C) for an unstable SISO (single-input, single-output) magnetic levitation system with two electromagnets (MLS2EM). Magnetic levitation systems offer various advantages, including friction-free, reliable, fast, and cost-effective operations. However, due to their unstable and highly nonlinear nature, these systems require sophisticated feedback control techniques to ensure optimal performance and functionality. To address these challenges, in this study, we derive the nonlinear state-space mathematical model of the MLS2EM and linearize it around five different operating points. The PD-GS-C controller aims to stabilize the system and improve steady-state control error. The strategy for obtaining the PD controller gains involves a parameter space technique, which specifies performance requirements. This technique results in ranges of proportional ( K P ) and derivative ( K D ) gains that are used by the PD-GS-C structure. To optimize the controller’s performance further, we utilize the big bang–big crunch optimization technique (BB-BC) to determine the optimal PD gains within the specified ranges. The optimization process focuses on achieving optimal performance in terms of a specific performance index function. This function quantifies the system’s time-domain step response criteria, which include minimizing overshoot percentage, settling time, and rising time. The index function is inversely proportional to the desired performance criteria, meaning that the goal is to maximize the index function to optimize the system’s performance. To validate the effectiveness and viability of the proposed strategy, we conducts MATLAB simulations and real-time experiments. The simulations and experimental findings serve to demonstrate the controller’s performance and verify its capabilities in stabilizing the MLS2EM magnetic levitation system. Keywords: magnetic levitation system; proportional-derivative controller; robust PD gain-scheduling controller (PD-GS-C); parameter space technique; Big Bang–Big Crunch optimization technique (BB-BC); real-time experiment MSC: 93B51 1. IntroductionA magnetically levitated mechanism (MAGLEV), also called magnetic suspension, is a naturally occurring phenomenon whereby an antigravity ferromagnetic item is maintained without any other sources except an electromagnetic field with no interaction with the external field [1]. Considering the sophistication of electromagnetic theory, this phenomenon has a broad range of applications, such as high-speed trains, frictionless melting, and magnetic bearings. With the rapid advancement of electronics and advanced control theory, magnetic levitation can be applied for a variety of practical purposes [2]. MAGLEV is a nonlinear, open-loop, unstable system that is challenging to model and control. Magnetic levitation systems are controlled as SISO system architectures, and just one electromagnet in the simplest and most common arrangement [3]. Such systems (MLS1EM) utilize magnetic forces created by providing the electromagnet and one-axis distance sensor with the necessary voltage [4]. In more complex arrangements (MLS2EM), a secondary electromagnet can be used to produce an external gravitational pull that is distinct from that produced by the first [5].For magnetic levitation systems, various well-known linear and nonlinear control techniques have been developed. Highly nonlinear operating differential equations make nonlinear controllers more desirable. To provide excellent productivity, many nonlinear control schemes require precise understanding of nonlinear plant effects. Practical implementations of nonlinear controllers are made challenging by the modeling of uncertain parameters in the MAGLEV model. On the other hand, the approximating linear model obtained by adjusting the system behavior around a preferred operating point is used in linear controller design, and as a result, the controllers are typically valid around the operating point. Linear controller performance is enhanced by using a gain-scheduling process to accommodate different operating points and enlarge the operating region [6].Numerous efforts have recently been published in the literature to control magnetic levitation systems. Input–output and input-state linearization approaches are among the feedback linearization-based controllers that have been developed [7,8,9,10,11,12], in addition to the design of linear state feedback controls [13], observer-based controls [14,15], cognitive strategies [16,17,18], slide mode control [19,20,21], backstepping control [22], model-based predictive control (MPC) [23,24], adaptive H-infinity tracking control [25], and proportional–integral–derivative (PID) control [4,26]. In addition, there are several adaptive control schemes that have already been presented in the past few years for positioning control schemes and magnetic levitation technologies. Some of these constructions can be found in [27,28,29,30,31,32,33,34,35]. An interesting method to consider for magnetic levitation positioning is that of finite time control, such as those mentioned in [36].Gain-scheduling controlled systems are common in today’s adaptive control applications, and they are succinctly examined as follows. Slide mode control and fuzzy-based gain scheduling of exact feed-forward linearization control for a magnetic levitation system are described in [37]. The authors in [38] discusse a high-gain adaptive output feedback control for a maglev apparatus. Using a Smith predictor, The authors of [39] provided a proportional–integral–derivative (PID) gain-scheduling control with second-order linear parameter variation, which eliminates time-varying delay. Additional intriguing adaptive gain-scheduling control mechanisms for practical, applications are discussed in [40,41,42,43,44]. The authors of [45] presented the design of a proportional–integral gain-scheduling control for position control of a sphere in a laboratory magnetic levitation system. The unstable and nonlinear mathematical model of the process is linearized at seven operating points. The authors of [46] presented a gain-scheduling control design procedure for classical proportional–integral–derivative controllers (PID-GS-C) for a positioning system. The method was applied to a laboratory magnetic levitation system with two electromagnets (MLS2EM), which allowed for several experimental verifications of the proposed solution.The authors of [47] provided a coherent approach to the tuning of PID controllers. By using the model-reference robust tuning (MoReRT) method, the authors proposed a way to achieve robust and effective PID controller tuning. The authors in [48] delved into deterministic multiobjective optimization methods, which are used to solve control system design problems that involve multiple conflicting objectives. The authors developed algorithms for two-degree-of-freedom PID control. An optimal indirect internal model control-proportional derivative (iIMC-PD) controller was designed as a two-loop control structure to regulate the dynamics of integrating processes with dead time [49]. The authors in [50,51] focused on analyzing optimal controller settings for controllers with a proportional–integral–derivative (PID) structure. Analysis was conducted by considering two key factors: the operating mode of the control loop (either servo or regulation mode) and the tuning mode of the controller. The goal was to determine controller settings that lead to optimal performance based on the specific operating and tuning modes.The position control of a ferromagnetic sphere in a laboratory magnetic levitation system with two electromagnets (MLS2EM) is the topic of this research. The control input is provided only to the upper electromagnet, and the bottom electromagnet can be subjected to the generated external disturbance. Two control architectures are developed to stabilize the system: proportional-derivative control (PD-C) and proportional-derivative gain-scheduling control (PD-GS-C). To enable the use of straightforward control techniques, the given mathematical model of the MLS2EM nonlinear system is linearized around five operating points. This linearization simplifies the control design process by allowing for the application of conventional control methods. For the PD-C controller, a parameter space technique is employed to design five robust PD controllers. The goal is to guarantee system stabilization and improve the steady-state control error. The parameter space approach provides ranges of proportional ( K P ) and derivative ( K D ) gains, which are crucial for the PD-GS-C structure. In the PD-GS-C architecture, the control gains are adaptively adjusted to match changes in the operating points. This adaptability is important because the bottom electromagnet of the MLS2EM system may be subjected to external disturbances, which could affect the system’s behavior and performance. To optimize the PD-GS-C controller’s performance, we use the big bang–big crunch optimization technique (BB-BC) to determine the optimal PD gains within the specified ranges obtained by the parameter space approach. The goal of optimization is to achieve the best possible performance for MLS2EM magnetic levitation. Finally, this research aims to demonstrate effective position control for a ferromagnetic sphere in an MLS2EM system using both PD-C and PD-GS-C controllers. By considering external disturbances and adapting the control gains to changing operating conditions, the proposed control architectures seek to stabilize the system and ensure accurate positioning of the ferromagnetic sphere. The use of linearization, parameter space techniques, and optimization methods contributes to the achievement of robust and efficient control of an MLS2EM magnetic levitation system.The structure of this article is described as follows. Potential models for the MAGLEV system are described in Section 2, including an explanation of the physical principles of magnetic levitation and how the system can be mathematically modeled. The parameter space approach, which sets the boundaries for the stabilization of PD gains, is described in Section 3. The utilization of the big bang–big crunch optimization technique (BB-BC) to obtain an optimal set of PD gains is also discussed in the same section. The PD gain-scheduling control (PD-GS-C) is presented in detail in Section 4, and the concept and architecture of the controller are explained. Section 5 contains illustrations and discussions of the experimental outcomes obtained through simulations, including graphs and charts that showcase the performance of the proposed PD-GS-C controller in the magnetic levitation system. Finally, Section 6 presents the key observations and findings of this research; we summarize the achievements of the study, discuss the effectiveness of the proposed PD-GS-C controller, and highlight any potential areas for further improvement or research. 2. Magnetic Levitation and ModelingFigure 1 illustrates the ML2SEM magnetic levitation equipment, which serves as an experimental benchmark for the theoretical validation of the magnetic levitation principle by suspending a ball of steel in the air while applying magnetic force. Due to its significant nonlinearity, this system is suited for use in university laboratories and is safe for assessment of the effectiveness of various controller types [52].A schematic diagram of the ML2SEM-controlled plant is shown in Figure 2, where EM1 and EM2 are the top and bottom electromagnets, respectively; m is the sphere’s mass; F e m 1 and F e m 2 are the electromagnetic forces; and F g is the force of gravity [52].The state-space mathematical model of the ML2SEM’s is expressed as: x 1 ˙ = x 2 x 2 ˙ = − F e m 1 ( x 1 , x 3 ) m + g + F e m 2 ( x 1 , x 4 ) m x 3 ˙ = 1 f i ( x 1 ) ( k i u 1 + c i − x 3 ) x 3 ˙ = 1 f i ( x d − x 1 ) ( k i u 2 + c i − x 4 ) The nonlinear functions in Equation (1) can be described as follows. F e m 1 ( x 1 , x 3 ) = x 3 2 F e m p 1 F e m p 2 e x p ( − x 1 F e m p 2 ) F e m 2 ( x 1 , x 4 ) = x 4 2 F e m p 1 F e m p 2 e x p ( − d − x 1 F e m p 2 ) f i ( x 1 ) = f i p 1 f i p 2 e x p ( − − x 1 f i p 2 ) where, x 1 ∈ [ 0 , 0.016 m ] is the sphere position, x 2 ∈ ℜ is the speed of the sphere, and x 3 and x 4 ∈ [ 0.03884 , 2.38 ] are the currents in the electromagnetic coil in amperes. The control signals that are applied to EM1 and EM2 are u 1 and u 2 ∈ [ 0.00498 , 1 ] , respectively; d denotes the spacing between electromagnets minus the diameter of the ball; f g = m g is the gravitational force; m is the ball mass; and g is the gravitational acceleration. The parameters k 1 and c i are related to the actuator nonlinear dynamic, and f i ( x 1 ) is a function of x 1 for both actuators. Table 1 shows the numerical values for the variables in Equations (1) and (2).The ferromagnetic ball is depicted in Figure 2 between EM1 and EM2. The basic objective is to achieve levitation of the ball at the location specified by the input reference signal. The stabilization and tracking of the input reference signal are the control system’s first two objectives [53]. multi-input–multioutput (MIMO) control systems can be created while both electromagnets are in use by adding an additional force from the control signal applied to EM2. Robust applications can benefit from this capability as well. As a disturbance input, the control signal ( u 2 ) can considered as the opposite. As a result, the third goal of the control system is to guarantee stability with regard to this kind of disturbance.Due to the nonlinearity of the controlled plant illustrated in Equations (1) and (2), establishing control systems that can achieve these three control objectives is rather difficult. In order to enable low-cost automation solutions, we overcome this challenge by linearizing around five operating points (o.p.s), utilizing coordinates ( x 1 ( j ) , x 2 ( j ) , x 3 ( j ) , x 4 ( j ) ) , where j is the operating point’s index and j = 1 … 5 . The general linearized state-space mathematical model can be applied in this situation by accepting u 2 = 0 , as follows [46]: Δ x ˙ ( j ) = A ( j ) Δ x ( j ) + B ( j ) Δ u 1 ( j ) Δ y ( j ) = C T ( j ) Δ x ( j ) where the differences between variables u 1 ( j ) and y ( j ) with regard to their numbers corresponding to the equilibrium point ( u 10 ( j ) and y 0 ( j ) , respectively) are represented by the formulas Δ u 1 ( j ) = u 1 ( j ) − u 10 ( j ) and Δ y ( j ) = y ( j ) − y 0 ( j ) ; y ( j ) = x 1 ( j ) is the output; the state vector is Δ x ( j ) = [ Δ x 1 ( j ) , Δ x 2 ( j ) , Δ x 3 ( j ) , Δ x 4 ( j ) ] T ; and matrix transposition is indicated by superscript T [53].Equation (3) can be expressed in general matrix form as follows: A ( j ) = 0 1 0 0 a 21 ( j ) 0 a 23 ( j ) a 24 ( j ) a 31 ( j ) 0 a 33 ( j ) 0 a 41 ( j ) 0 0 a 44 ( j ) , B ( j ) = 0 0 b 3 ( j ) b 4 ( j ) , C T ( j ) = 1 0 0 0 , Equation (5) describes the elements of matrices A ( j ) and B ( j ) in Equation (4). a 21 ( j ) = x 30 2 ( j ) m F e m p 1 F e m p 2 2 e x p ( − x 10 ( j ) F e m p 2 ) + x 40 2 ( j ) m F e m p 1 F e m p 2 2 e x p ( − d − x 10 ( j ) F e m p 2 ) , a 23 ( j ) = − 2 x 30 ( j ) m F e m p 1 F e m p 2 e x p ( − x 10 ( j ) F e m p 2 ) , a 24 ( j ) = 2 x 40 ( j ) m F e m p 1 F e m p 2 e x p ( − d − x 10 ( j ) F e m p 2 ) , a 31 ( j ) = − ( k i u 1 ( j ) + c i − x 30 ( j ) ) ( x 10 ( j ) / f i p 2 ) f i − 1 ( x 10 ( j ) ) , a 33 ( j ) = − f i ( − 1 ) ( x 10 ( j ) ) , a 41 ( j ) = − ( k i u 1 ( j ) + c i − x 40 ( j ) ) ( x 10 ( j ) / f i p 2 ) f i − 1 ( d − x 10 ( j ) ) , a 44 ( j ) = − f i ( − 1 ) ( d − x 10 ( j ) ) , b 3 ( j ) = k i f i ( − 1 ) ( x 10 ( j ) ) , b 4 ( j ) = k i f i ( − 1 ) ( d − x 10 ( j ) ) , The operating points in Equation (6) are selected in order to cover the typical operating conditions within the steady-state area of the ball position feedback input–output map and eliminate the input–output map’s extreme values due to potential instability [46]. P ( 1 ) ( 0.0063 , 0 , 1.22 , 0.39 ) , P ( 2 ) ( 0.007 , 0 , 1.145 , 0.39 ) , P ( 3 ) ( 0.0077 , 0 , 1.07 , 0.39 ) , P ( 4 ) ( 0.0098 , 0 , 0.89 , 0.39 ) , P ( 5 ) ( 0.009 , 0 , 0.9345 , 0.39 ) , 3. Stability AnalysisThe proportional-derivative gain-scheduling Control system (PD-GS-C) solution based on the control system structure is shown in Figure 3, which can be considered the stabilizing control solution starting with the linearized state-space model in Equation (4).Where the reference input is defined by r, the PD gain parameters are K p and K D , and the ML2SEM model is described by the plant, which also contains the actuators and sensors. The minimum and maximum values of the PD gain controllers are typically set to limit the control action within certain bounds. By setting these limits, the control system can prevent excessive control efforts that could potentially destabilize the system or cause other undesirable effects. These values are used in equations of the proportional-derivative gain-scheduling control system (PD-GS-C), as shown in Section 4.2. The applied technique is based on a performance specification in terms of a parameter-space technique in order to determine the maximum and minimum of the PD gain controllers ( K p m a x , K D m a x , K p m i n , and K D m i n ), as shown in the next subsections. Stabilizing Controller GainsA parameter-space approach can be defined by the method of drawing lines that indicate the boundaries of the stabilizing gain set using the frequency response of the plant [54]. G ( s ) is a single-input–single-output plant’s transfer function with a frequency response that can be illustrated by the following equation: G ( j ω ) = a ( ω ) + j b ( ω ) , ω ∈ R where a , b are real-valued functions.Equation (8) shows the transfer function of the proportional-derivative controller plant G p d ( s ) . G p d ( j ω ) = K p + j ω K D , ω ∈ R The characteristic equation of a closed-loop system with a PD controller is expressed as: 1 + ( K p + j ω K D ) ( a ( ω ) + j b ( ω ) ) = 0 Equation (9) can be rewritten in matrix form as: a − ω b b ω a K P K D = − 1 0 By solving Equations (9) or (10), K p and K D can be expressed by the following equation: K p = − a a 2 ( ω ) + b 2 ( ω ) K D = b ω ( a 2 ( ω ) + b 2 ( ω ) ) When ω changes, the parameter plane can be separated into regions by drawing Equation (11) and the root real boundary line (at ω = 0 ). The number of unstable poles in each region is guaranteed to be kept fixed by the boundary-crossing theorem [55]. This characteristic allows us to calculate the stability set at each operating point [54,56,57]. The K P − K D planes for operating points p ( 1 ) , p ( 2 ) , p ( 3 ) , p ( 4 ) , and p ( 5 ) are illustrated in Figure 4a–e, respectively. Table 2 demonstrates the stability of each zone for operating points p ( 1 ) , p ( 2 ) , p ( 3 ) , p ( 4 ) , and p ( 5 ) by picking one polynomial for each region and testing its stability. The numbers (0, 1, 2, and 3) that appear in figures represent the number of unstable poles in the zone. The green zones in all figures represent the stable regions with zero unstable poles. In this research, the robust PD controller for p ( j ) is involved by computing all stabilizing PD controller gains for the characteristic polynomial using the parameter-space approach, as shown in Figure 4. 4. Controller DesignThe results of the parameter-space approach are ranges of PD gains ( K P , K D ) . In the next subsection, the optimal PD gains were chosen using the Big Bang–Big Crunch Optimization Algorithm [58] to achieve optimal performance for ML2SEM magnetic levitation. 4.1. The Big Bang–Big Crunch Optimization AlgorithmThe Big Bang–Big Crunch Optimization technique (BB-BC) is a population-oriented optimized algorithm that derives motivation from the cosmological events of the Big Bang and the Big Crunch. As a nature-inspired method for resolving challenging optimization issues, it was initially proposed by Erol and Eksin in 2006 [58]. The technique begins with a random population of candidate solutions that indicate potential solutions to the optimization challenge. Every potential solution is represented in the search space by a point. Depending on the scenario, the search space can be multidimensional. The BB-BC algorithm can be separated into two phases: the Big Bang phase and the Big Crunch phase.The Big Bang Phase: During the Big Bang phase, the candidate solutions are randomly spread throughout the search space, comparable to the exploded atoms during the Big Bang in cosmology. This phase is in control of exploration, allowing the solutions to cover a large portion of the search space.The Big Crunch Phase: Following the Big Bang phase, the algorithm initiates the Big Crunch phase, which depicts the contraction of matter in the universe, resulting in the construction of structures. During this phase, potential solutions are drawn to promising parts of the search space, identical to the development of galaxies and other cosmic structures. This phase focuses on exploitation, refining solutions, and convergence toward optimal solutions.The output of the Big Crunch phase can be defined as the center of mass. The center of mass is denoted by x c → , which can be expressed as follows: x c → = ∑ i = 1 N 1 f i x i → ∑ i = 1 N 1 f i where x i → and x c → are vectors in an n-dimensional search space, f i is the value of this point’s objective function, and N is the population size in the Big Bang phase.The following stage generates new points that are used in the Big Bang phase following the Big Crunch phase, generating the center of mass ( x c → ). The newly generated points are redistributed in all directions around the center of mass ( x c → ): x n e w = x c → + r α ( x m a x − x m i n ) k where r is a random number, α is a parameter that limits the size of the search space, and k is the number of iterations. [58,59,60,61,62].Iterations during the Big Crunch phase drive candidate solutions closer to the best solutions in their area. Attractive forces guide the movements, which are determined by the fitness values of the solutions. Based on the performance function, the fitness function evaluates the level of accuracy of each solution.The optimization approach is presented via a specific performance index function that is inversely proportional to a dynamical system’s time-domain step response criteria (small overshoot percentage ( O S ) with significant minimization of both settling ( T s ) and rising ( T r ) times), which is demonstrated by Equation (14). A flow chart of the BB-BC optimization algorithm is shown in Figure 5. J = 100 0.05 O S + 6 T s + 12 T r The search domains for [ K P K D ] are shown in Figure 4 for each operating point. Table 3 demonstrates that the performance varies in terms of overshoot ( O S ), rise time ( T r ), settling time ( T s ), peak time ( T P ), and steady-state error ( S S E ) for the five operating points. In addition, the optimal values of the PD gain controller, as well as the maximum and minimum gains for each of the five operating points, are shown in Table 3. Furthermore, Figure 6 shows the step response of the optimal PD gains controllers obtained by employing the BB-BC algorithm for each of the five operating points. 4.2. The Proportional-Derivative Gain-Scheduling Controller (PD-GS-C)According to [63], swapping from one PD control scheme to another is a sensible way to guarantee improved performance, although process nonlinearities depend on the operating point as well. This method is defined as the PD gain schedule. Figure 3 depicts the gain-scheduling PD voltage controller used in magnetic levitation. It is mathematically represented by Equation (15). Δ u 1 ( t ) = K P e ( t ) + K D d e ( t ) d t where the proportional gain ( K P ( t ) ) and the derivative gain ( K D ( t ) ) are functions of e ( t ) such that e ( t ) = r ( t ) − y ( t ) . The following mathematical formula can be used to define the proportional gain ( K P ( t ) ) and the derivative gain ( K D ( t ) ) as a function of e(t) [46]: K P ( t ) = K P m a x − ( K P m a x − K P m i n ) e x p ( − ( α ( t ) | e ( t ) | ) ) , K D ( t ) = K D m a x − ( K D m a x − K D m i n ) e x p ( − ( α ( t ) | e ( t ) | ) ) . where 0 ≤ α ( t ) ≤ 1 and K P m a x , K P m i n , K D m a x , and K D m i n are the maximum and minimum values of the proportional gain ( K P ) and the derivative gain ( K D ), respectively.Equation (16) demonstrates that the exponential term approaches zero ( e x p ( − ( α ( t ) | e ( t ) | ) ) → 0 ) when e ( t ) is large. Consequentially, K P ( t ) = K P m a x . Similarly, when the error ( e ( t ) ) is minimal, the exponential term approaches one ( e x p ( − ( α ( t ) | e ( t ) | ) ) → 1 ) , and as a result, K P ( t ) = K P m i n .In this modification, K P and K D are evaluated by the control scheme so that when the error ( e ( t ) ) is large, we predict maximum proportional and derivative gains ( K P m a x and K D m a x ) to enhance the control signal’s voltage up to the dynamic response, and when the error ( e ( t ) ) is small, we anticipate minimal proportional derivative gains ( K P m i n and K D m i n ) to overcome the unfavorable issue of overshoot. These fluctuation in the gains are generated by the variable “ α ( t ) ” to determine how quickly K P ( t ) and K P ( t ) change between the maximum and minimum values [63].In order to transfer from one PD controller to another using a smooth and steady method, the parameter 0 ≤ α ( t ) ≤ 1 is included. To determine this parameter’s value, the following equation is applied: α ( t ) = tanh ( η β ( t ) ) = ( e x p ( 2 η β ( t ) ) − 1 ) ( e x p ( 2 η β ( t ) ) + 1 ) where the factor η was selected to guarantee specific dynamics of the fluctuation of α ( t ) and specifies the rate at which α ( t ) varies from 0 to 1. The variable β ( t ) is defined in terms of | e ( t ) | and r ( t ) [46]: β ( t ) = 1 , for | e ( t ) | > ξ 0 , for | e ( t ) | ≤ ξ , ξ = 0.9 . r ( t ) where the parameter ξ is set according to [46]. 5. Simulation ResultsSimulink in MATLAB(R2015a) was used to simulate the two control mechanisms: the optimal PD-C and PD-GS-C. Figure 7 demonstrates the control mechanism of PD-GS-C by Simulink(R2015a). The reference input in each scenario was a signal that was set to 0.008 m from the upper coil, and the responses of the controls were examined during a time period of 5 s. Graphs of the top electromagnet’s responses to the ball location, ball velocity, coil current, and control signal are presented.A PD controller with five operating points was simulated in MATLAB(R2015a). The simulation results of the control system with a PD controller developed specifically for the two o.p.s ( ( 1 ) and ( 4 ) ) are presented in Figure 8a, which shows the control signal over time. The time-varying coil current through EM1 is displayed in Figure 8b, the ball velocity is shown in Figure 8c, and Figure 8d demonstrates the controlled output position of ball over time.The results show that the system’s nonlinear effects generate oscillations at the start of transient responses.The specifications of the PD-GS-C component in Figure 3 have the variables K P m a x , K P m i n , K D m a x , and K D m i n . By using the optimal values from Table 3, η = ∈ [ 0.001 , 0.1 ] . These values were obtained using the big bang–big crunch optimization technique (BB-BC) in the context of a tradeoff for overshoot for the five operating points.The conditions listed below can be taken into consideration when designing the PD-GS-C structure: K P m i n and K D m i n are enabled during the steady state to achieve a minimal value of steady-state error ( | e ( t ) | ) to resolve the unacceptable challenge of overshoot as shown in Figure 9. K P m a x and K D m a x are generated when the process steady-state error | e ( t ) | is large in order to provide a significant control signal and to mitigate unfavorable fluctuation and error, as illustrated in Figure 9.Figure 10 illustrates responses to various variables estimated in the proposed scheme. The simulation results of the control system using the PD-GS controller are presented in Figure 10a, which shows the control signal over time. The time-varying coil current through EM1 is displayed in Figure 10b, and the ball velocity is shown in Figure 10c. Moreover, Figure 10d shows the controlled output over time.Due to the PD-GS controller’s optimized dynamic response, which is characterized by fewer oscillations and less settling time, its results are superior to those shown in Figure 8. MATLAB simulation and nonlinear laboratory equipment were used to test and validate the PD-C and PD-GS-C control structures in real-time experiments. Figure 11 shows the MATLAB simulation results of a robust proportional-derivative gain-scheduling controller (PD-GS -C).Figure 12 displays the real-time experimental results of a robust Proportional-derivative gain-scheduling controller (PD-GS-C). Significant agreement was observed between theoretical and practical observations. 6. ConclusionsTwo control schemes (PD-C and (PD-GS-C) were been designed to control the orientation of a ball in an MLS2EM system. The PD-C (proportional-derivative controller) is a straightforward control scheme that uses the proportional and derivative terms to control the system’s orientation. The PD controller was designed based on a linearized model of MLS2EM around five operating points, allowing for the application of conventional control techniques. The PD-GS-C (proportional-derivative with gain scheduling controller) was designed to switch between multiple PD controllers based on specific conditions. To achieve this, a parameter-space approach was employed to design five robust PD controllers with varying gains ( K P and K D ). The ranges of PD gains obtained using the parameter-space approach were used in the PD-GS-C structure, allowing the controller to adapt its gains according to different operating conditions. To find the optimal PD gains for the PD-GS-C controller, we utilized the big bang–big crunch optimization technique (BB-BC). This optimization method was applied to improve the system’s performance in terms of time-domain step response criteria, which include minimizing overshoot percentage, settling time, and rising time. The performance index function was defined to quantify the desired performance criteria, and the optimization process aims to maximize this index function, thereby optimizing the system’s performance. To evaluate the effectiveness of the proposed control schemes, both PD-C and PD-GS-C configurations were tested using Simulink/MATLAB simulations, as well as real-time experiments. By comparing the PD-GS controller with the conventional PD controller, this research demonstrates the efficacy of the PD-GS-C scheme in achieving optimal control performance for an MLS2EM magnetic levitation system. As the gain-scheduling approach provided good results, we intend to extend this idea to other controller configurations, such as PI and PID controls. Therefore, in future research, a comparison of gain-scheduled PI/PD/PID controllers will be conducted under the same scenario of operating points and controller scheduling. An interesting approach to be considered in future research is that of finite-time control. Finite-time control strategies aim to achieve desired control objectives within a specified time frame, which can be particularly beneficial for systems wherein rapid and precise control is crucial. Magnetic levitation systems often require precise control to maintain stable levitation and positioning of objects, making them a suitable candidate for such research. Author ContributionsConceptualization, M.A., R.V. and H.S.A.-N.; methodology, M.A. and H.S.A.-N.; software, M.A. and H.S.A.-N.; validation, R.V., O.A. and M.A.; formal analysis, M.A. and H.S.A.-N.; investigation, O.A. and R.V.; resources, M.A.; data curation, O.A.; writing—original draft preparation, M.A. and H.S.A.-N.; writing—review and editing, R.V.; visualization, M.A. and R.V.; supervision, M.A., O.A. and R.V.; project administration, M.A.; funding acquisition, R.V. All authors have read and agreed to the published version of the manuscript.FundingThis research received no external funding.Data Availability StatementNot applicable.AcknowledgmentsThis work received support from the Catalan Government under Project 2023 SGR 00197 and from the Spanish Government under MICINN projects TED2023-129134B-I00 and PID2023-105434RBC33 cofunded by the European Union ERDF funds. Financial support was also provided by the University of Costa Rica under grant 731-B9-265.Conflicts of InterestThe authors declare no conflict of interest.AbbreviationsThe following abbreviations are used in this manuscript: PD-CProportional-derivative controllerPD-GS-CProportional-derivative gain-schedulingcontrollerBB-BCBig Bang–Big Crunch optimizationReferencesSinha, P.K. Electromagnetic Suspension Dynamics & Control; Savoy Place: London, UK, 1987. [Google Scholar]Czerwiński, K.; Ławryńczuk, M. Identification of discrete-time model of active magnetic levitation system. In Proceedings of the Trends in Advanced Intelligent Control, Optimization and Automation: Proceedings of KKA 2017—The 19th Polish Control Conference, Kraków, Poland, 18–21 June 2017; pp. 599–608.Sun, Z. Magnetic Levitation Based on Switched Reluctance Actuator. Ph.D. Thesis, Hong Kong Polytechnic University, Hong Kong, China, 2011. [Google Scholar]Almobaied, M.; Al-Nahhal, H.S.; Issa, K.B. Computation of stabilizing PID controllers for magnetic levitation system with parametric uncertainties. In Proceedings of the 2023 International Conference on Electric Power Engineering–Palestine (ICEPE-P), Gaza, Palestine, 23–24 March 2023; pp. 1–7. [Google Scholar]Bojan-Dragos, C.A.; Stinean, A.I.; Precup, R.E.; Preitl, S.; Petriu, E.M. Model predictive control solution for magnetic levitation systems. In Proceedings of the 2015 20th International Conference on Methods and Models in Automation and Robotics (MMAR), Miedzyzdroje, Poland, 24–27 August 2015; pp. 139–144. [Google Scholar]Eroğlu, Y.; Ablay, G. Cascade sliding mode-based robust tracking control of a magnetic levitation system. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2016, 230, 851–860. [Google Scholar] [CrossRef]El Hajjaji, A.; Ouladsine, M. Modeling and nonlinear control of magnetic levitation systems. IEEE Trans. Ind. Electron. 2001, 48, 831–838. [Google Scholar] [CrossRef]Barie, W.; Chiasson, J. Linear and nonlinear state-space controllers for magnetic levitation. Int. J. Syst. Sci. 1996, 27, 1153–1163. [Google Scholar] [CrossRef]Beltran-Carbajal, F.; Valderrabano-Gonzalez, A.; Rosas-Caro, J.; Favela-Contreras, A. Output feedback control of a mechanical system using magnetic levitation. ISA Trans. 2015, 57, 352–359. [Google Scholar] [CrossRef] [PubMed]Morales, R.; Sira-Ramírez, H. Trajectory tracking for the magnetic ball levitation system via exact feedforward linearisation and GPI control. Int. J. Control 2010, 83, 1155–1166. [Google Scholar] [CrossRef]Nielsen, C.; Fulford, C.; Maggiore, M. Path following using transverse feedback linearization: Application to a maglev positioning system. Automatica 2010, 46, 585–590. [Google Scholar] [CrossRef]Zhang, J.; Tao, T.; Mei, X.; Jiang, G.; Zhang, D. Non-linear robust control of a voltage-controlled magnetic levitation system with a feedback linearization approach. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2011, 225, 85–98. [Google Scholar] [CrossRef]El Rifai, O.M.; Youcef-Toumi, K. Achievable performance and design trade-offs in magnetic levitation control. In Proceedings of the AMC’98-Coimbra—1998 5th International Workshop on Advanced Motion Control, Coimbra, Portugal, 29 June–1 July 1998; pp. 586–591. [Google Scholar]Hasirci, U.; Balikci, A.; Zabar, Z.; Birenbaum, L. A novel magnetic-levitation system: Design, implementation, and nonlinear control. IEEE Trans. Plasma Sci. 2010, 39, 492–497. [Google Scholar] [CrossRef]Baranowski, J.; Piątek, P. Observer-based feedback for the magnetic levitation system. Trans. Inst. Meas. Control 2012, 34, 422–435. [Google Scholar] [CrossRef]Al-Araji, A.S. Cognitive non-linear controller design for magnetic levitation system. Trans. Inst. Meas. Control 2016, 38, 215–222. [Google Scholar] [CrossRef]Benitez-Pérez, H.; Ortega-Arjona, J.; Cardenas-Flores, F.; Quiñones-Reyes, P. Reconfiguration control strategy using Takagi–Sugeno model predictive control for network control systems—A magnetic levitation case study. Proc. Inst. Mech. Eng.-Part I J. Syst. Control Eng. 2010, 224, 1022. [Google Scholar] [CrossRef]Chen, S.Y.; Lin, F.J.; Shyu, K.K. Direct decentralized neural control for nonlinear MIMO magnetic levitation system. Neurocomputing 2009, 72, 3220–3230. [Google Scholar] [CrossRef]Cho, D.; Kato, Y.; Spilman, D. Sliding mode and classical controllers in magnetic levitation systems. IEEE Control Syst. Mag. 1993, 13, 42–48. [Google Scholar]Al-Muthairi, N.; Zribi, M. Sliding mode control of a magnetic levitation system. Math. Probl. Eng. 2004, 2004, 93–107. [Google Scholar] [CrossRef]Lin, F.J.; Chen, S.Y.; Shyu, K.K. Robust dynamic sliding-mode control using adaptive RENN for magnetic levitation system. IEEE Trans. Neural Networks 2009, 20, 938–951. [Google Scholar] [PubMed]Lin, F.J.; Teng, L.T.; Shieh, P.H. Intelligent adaptive backstepping control system for magnetic levitation apparatus. IEEE Trans. Magn. 2007, 43, 2009–2018. [Google Scholar] [CrossRef]Bächle, T.; Hentzelt, S.; Graichen, K. Nonlinear model predictive control of a magnetic levitation system. Control Eng. Pract. 2013, 21, 1250–1258. [Google Scholar] [CrossRef]Iplikci, S.; Bahtiyar, B. A field-programmable gate array implementation of a real-time nonlinear Runge–Kutta model predictive control. Trans. Inst. Meas. Control 2016, 38, 555–564. [Google Scholar] [CrossRef]Javadi, A.; Pezeshki, S. A new model-free adaptive controller versus non-linear H∞ controller for levitation of an electromagnetic system. Trans. Inst. Meas. Control 2013, 35, 321–329. [Google Scholar] [CrossRef]Kumar, E.V.; Jerome, J. LQR based optimal tuning of PID controller for trajectory tracking of magnetic levitation system. Procedia Eng. 2013, 64, 254–264. [Google Scholar] [CrossRef]Wang, B.; Liu, G.P.; Rees, D. Networked predictive control of magnetic levitation system. In Proceedings of the 2009 IEEE International Conference on Systems, Man and Cybernetics, San Antonio, TX, USA, 11–14 2009; pp. 4100–4105. [Google Scholar]An, S.; Ma, Y.; Cao, Z. Applying simple adaptive control to magnetic levitation system. In Proceedings of the 2009 Second International Conference on Intelligent Computation Technology and Automation, Zhangjiajie, China, 10–11 October 2009; Volume 1, pp. 746–749. [Google Scholar]Ziętkiewicz, J. Constrained predictive control of a levitation system. In Proceedings of the 2011 16th International Conference on Methods & Models in Automation & Robotics, Miedzyzdroje, Poland, 22–25 August 2011; pp. 278–283. [Google Scholar]Chauhan, S.; Nigam, M. Model predictive controller design and perturbation study for magnetic levitation system. In Proceedings of the 2014 Recent Advances in Engineering and Computational Sciences (RAECS), Chandigarh, India, 6–8 March 2014; pp. 1–6. [Google Scholar]King, R.; Stathaki, A. Fuzzy gain-scheduling control of nonlinear processes. In Proceedings of the WSES/IEEE/IMACS World Multiconference on Circuits, Systems, Communications and Computers CSCC’99, Athens, Greece, 4–8 July 1999. [Google Scholar]Huang, Y.W.; Tung, P.C. Design of a fuzzy gain scheduling controller having input saturation: A comparative study. J. Mar. Sci. Technol. 2009, 17, 249–256. [Google Scholar] [CrossRef]Elsodany, N.M.; Rezeka, S.F.; Maharem, N.A. Adaptive PID control of a stepper motor driving a flexible rotor. Alex. Eng. J. 2011, 50, 127–136. [Google Scholar] [CrossRef]Bianchi, F.D.; Peña, R.S.S. Interpolation for gain-scheduled control with guarantees. Automatica 2011, 47, 239–243. [Google Scholar] [CrossRef]Zhao, Z.Y.; Tomizuka, M.; Isaka, S. Fuzzy gain scheduling of PID controllers. IEEE Trans. Syst. Man Cybern. 1993, 23, 1392–1398. [Google Scholar] [CrossRef]Ye, Z.; Zhang, D.; Deng, C.; Yan, H.; Feng, G. Finite-time resilient sliding mode control of nonlinear UMV systems subject to DoS attacks. Automatica 2023, 156, 111170. [Google Scholar] [CrossRef]Lashin, M.; Elgammal, A.T.; Ramadan, A.; Abouelsoud, A.; Assal, S.F.; Abo-Ismail, A. Fuzzy-based gain scheduling of Exact FeedForward Linearization control and sliding mode control for magnetic ball levitation system: A comparative study. In Proceedings of the 2014 IEEE International Conference on Automation, Quality and Testing, Robotics, Cluj-Napoca, Romania, 22–24 May 2014; pp. 1–6. [Google Scholar]Michino, R.; Tanaka, H.; Mizumoto, I. Application of high gain adaptive output feedback control to a magnetic levitation system. In Proceedings of the 2009 ICCAS-SICE, Fukuoka City, Japan, 18–21 August 2009; pp. 970–975. [Google Scholar]Puig, V.; Bolea, Y.; Blesa, J. Robust gain-scheduled Smith PID controllers for second order LPV systems with time varying delay. IFAC Proc. Vol. 2012, 45, 199–204. [Google Scholar] [CrossRef]Bianchi, F.D.; Sánchez-Peña, R.S.; Guadayol, M. Gain scheduled control based on high fidelity local wind turbine models. Renew. Energy 2012, 37, 233–240. [Google Scholar] [CrossRef]Yang, Y.; Yan, Y. Attitude regulation for unmanned quadrotors using adaptive fuzzy gain-scheduling sliding mode control. Aerosp. Sci. Technol. 2016, 54, 208–217. [Google Scholar] [CrossRef]Bedoud, K.; Ali-rachedi, M.; Bahi, T.; Lakel, R. Adaptive Fuzzy Gain Scheduling of PI Controller for Control of the Wind Energy Conversion Systems. Energy Procedia 2015, 74, 211–225. [Google Scholar] [CrossRef]Dounis, A.I.; Kofinas, P.; Alafodimos, C.; Tseles, D. Adaptive fuzzy gain scheduling PID controller for maximum power point tracking of photovoltaic system. Renew. Energy 2013, 60, 202–214. [Google Scholar] [CrossRef]Xie, Y.; Shi, H.; Alleyne, A.G.; Yang, B. Feedback shape control for deployable mesh reflectors using gain scheduling method. Acta Astronaut. 2016, 121, 241–255. [Google Scholar] [CrossRef]Bojan-Dragos, C.A.; Precup, R.E.; Preitl, S.; Hergane, S.; Hughiet, E.G.; Szedlak-Stinean, A.I. proportional–integral gain-scheduling control of a magnetic levitation system. In Proceedings of the 2016 20th International Conference on System Theory, Control and Computing (ICSTCC), Sinaia, Romania, 13–15 October 2016; pp. 1–6. [Google Scholar]Bojan-Dragos, C.A.; Precup, R.E.; Tomescu, M.L.; Preitl, S.; Tanasoiu, O.M.; Hergane, S. proportional–integral-derivative gain-scheduling control of a magnetic levitation system. Int. J. Comput. Commun. Control 2017, 12, 599–611. [Google Scholar] [CrossRef]Alfaro, V.M.; Vilanova, R. Model-Reference Robust Tuning of PID Controllers; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]Rojas, J.D.; Arrieta, O.; Vilanova, R. Industrial PID Controller Tuning; Springer: Berlin/Heidelberg, Germany, 2023. [Google Scholar]Aryan, P.; Raja, G.L.; Vilanova, R. Optimal iIMC-PD double-loop control strategy for integrating processes with dead-time. In Proceedings of the APCA International Conference on Automatic Control and Soft Computing, Caparica, Portugal, 6–8 July 2023; pp. 521–531. [Google Scholar]Arrieta, O.; Vilanova, R. Servo/Regulation tradeoff tuning of PID controllers with a robustness consideration. In Proceedings of the 2007 46th IEEE Conference on Decision and Control, New Orleans, LA, USA, 12–14 December 2007; pp. 1838–1843. [Google Scholar]Arrieta, O.; Vilaxova, R. Performance degradation analysis of controller tuning modes: Application to an optimal PID tuning. Int. J. Innov. Comput. Inf. Control 2010, 6, 4719–4729. [Google Scholar]Magnetic Levitation System 2EM (MLS2EM) User’s Manual; Inteco Ltd.: Sankt-Petersburg, Russia, 2008.Dragoş, C.A.; Precup, R.E.; Preitl, S.; Petriu, E.M.; Rădac, M.B. Control Solutions, Simulation and Experimental Results for a Magnetic Levitation Laboratory System. Available online: https://www.eurosim.info/fileadmin/user_upload_eurosim/EUROSIM_OA/Congress/2010/data/papers/155.pdf (accessed on 7 September 2023).Saeki, M. Properties of stabilizing PID gain set in parameter space. IEEE Trans. Autom. Control 2007, 52, 1710–1715. [Google Scholar] [CrossRef]Bhattacharyya, S.P.; Keel, L.H. Robust control: The parametric approach. In Advances in Control Education; Elsevier: Amsterdam, The Netherlands, 1995; pp. 49–52. [Google Scholar]Almobaied, M.; Al-Nahhal, H.S.; Issa, K.B. Robus t-Pro porti onal-Integ ral-D eriva tive Controller Design for Magnetic Levitation System Using Big Bang–Big Crunch Algorithm. Electrica 2023, 23, 270–280. [Google Scholar] [CrossRef]Almobaied, M.; Al-Nahhal, H.S.; Issa, K.B. Robust-PID controller design for magnetic levitation system using parameter space approach. World J. Adv. Eng. Technol. Sci. 2023, 8, 135–151. [Google Scholar] [CrossRef]Erol, O.K.; Eksin, I. A new optimization method: Big bang–big crunch. Adv. Eng. Softw. 2006, 37, 106–111. [Google Scholar] [CrossRef]Kumbasar, T.; Eksin, I.; Guzelkaya, M.; Yesil, E. Adaptive fuzzy model based inverse controller design using BB-BC optimization algorithm. Expert Syst. Appl. 2011, 38, 12356–12364. [Google Scholar] [CrossRef]Yılmaz, S.; Gökaşan, M. Optimal trajectory planning by big bang-big crunch algorithm. In Proceedings of the 2014 International Conference on Control, Decision and Information Technologies (CoDIT), Metz, France, 3–5 November 2014; pp. 557–561. [Google Scholar]Dincel, E.; Genc, V.I. A power system stabilizer design by big bang-big crunch algorithm. In Proceedings of the 2012 IEEE International Conference on Control System, Computing and Engineering, Penang, Malaysia, 23–25 November 2012; pp. 307–312. [Google Scholar]Almobaied, M.; Eksin, I.; Guzelkaya, M. Design of LQR controller with big bang-big crunch optimization algorithm based on time domain criteria. In Proceedings of the 2016 24th Mediterranean Conference on Control and Automation (MED), Athens, Greece, 21–24 June 2016; pp. 1192–1197. [Google Scholar]Sedaghati, A. A PI controller based on gain-scheduling for synchronous generator. Turk. J. Electr. Eng. Comput. Sci. 2006, 14, 241–251. [Google Scholar] Figure 1. ML2SEM magnetic levitation equipment. Figure 1. ML2SEM magnetic levitation equipment. Figure 2. Block diagram of the ML2SEM-controlled plant [52]. Figure 2. Block diagram of the ML2SEM-controlled plant [52]. Figure 3. Block diagram of the PD-GS-C system structure. Figure 3. Block diagram of the PD-GS-C system structure. Figure 4. (a) K P − K D plane with boundaries of P ( 1 ) ; (b) K P − K D plane with boundaries of P ( 2 ) ; (c) K P − K D plane with boundaries of P ( 3 ) ; (d) K P − K D plane with boundaries of P ( 4 ) ; (e) K P − K D plane with boundaries of P ( 5 ) . Figure 4. (a) K P − K D plane with boundaries of P ( 1 ) ; (b) K P − K D plane with boundaries of P ( 2 ) ; (c) K P − K D plane with boundaries of P ( 3 ) ; (d) K P − K D plane with boundaries of P ( 4 ) ; (e) K P − K D plane with boundaries of P ( 5 ) . Figure 5. Flow chart of the Big Bang–Big Crunch optimization algorithm for the proposed method. Figure 5. Flow chart of the Big Bang–Big Crunch optimization algorithm for the proposed method. Figure 6. Step response of the ML2SEM system for optimal proportion-derivative ( K P , K D ) values achieved by utilizing the Big Bang–Big Crunch algorithm for each of the five operating points. Figure 6. Step response of the ML2SEM system for optimal proportion-derivative ( K P , K D ) values achieved by utilizing the Big Bang–Big Crunch algorithm for each of the five operating points. Figure 7. Block diagram of the PD-GS-C by Simulink(R2015a). Figure 7. Block diagram of the PD-GS-C by Simulink(R2015a). Figure 8. (a) The PD control signal over time; (b) PD controller’s current iEM1 versus time; (c) ball velocity over time for the PD controller; (d) ball position over time for the PD controller. Figure 8. (a) The PD control signal over time; (b) PD controller’s current iEM1 versus time; (c) ball velocity over time for the PD controller; (d) ball position over time for the PD controller. Figure 9. (a) K P gain as a function of steady–state error ( | e ( t ) | ); (b) K D gain as a function of steady-state error | e ( t ) | . Figure 9. (a) K P gain as a function of steady–state error ( | e ( t ) | ); (b) K D gain as a function of steady-state error | e ( t ) | . Figure 10. (a) Control signal ( u E M 1 ) over time for the PD–GS controller; (b) PD–GS controller’s current i E M 1 versus time; (c) ball velocity over time for the PD–GS controller; (d) ball position over time for the PD–GS controller. Figure 10. (a) Control signal ( u E M 1 ) over time for the PD–GS controller; (b) PD–GS controller’s current i E M 1 versus time; (c) ball velocity over time for the PD–GS controller; (d) ball position over time for the PD–GS controller. Figure 11. MATLAB simulation results of a robust proportional-derivative gain-scheduling controller (PD-GS-C). Figure 11. MATLAB simulation results of a robust proportional-derivative gain-scheduling controller (PD-GS-C). Figure 12. Real–time experimental results of a robust proportional-derivative gain-scheduling controller (PD–GS–C); (a) ball position over time for the PD–GS controller; (b) ball velocity over time for the PD–GS controller; (c) PD–GS controller’s current i E M 1 versus time; (d) the pulses applied to coil2 as disturbance. Figure 12. Real–time experimental results of a robust proportional-derivative gain-scheduling controller (PD–GS–C); (a) ball position over time for the PD–GS controller; (b) ball velocity over time for the PD–GS controller; (c) PD–GS controller’s current i E M 1 versus time; (d) the pulses applied to coil2 as disturbance. Table 1. ML2SEM maglev model parameters [52]. Table 1. ML2SEM maglev model parameters [52]. ParameterParameter ValueUnitm 0.0571 [kg]g 9.81 [m/s2] F e m p 1 1.7521 × 10 − 2 [H] F e m p 2 5.8231 × 10 − 3 [m] f i p 1 1.4142 × 10 − 4 [ms] f i p 2 4.5626 × 10 − 3 [m] c i 0.0243 [A] k i 2.5165 [A] x d 0.075 [m] Table 2. The stability of each zone for operating points p ( 1 ) , p ( 2 ) , p ( 3 ) , p ( 4 ) , and p ( 5 ) . Table 2. The stability of each zone for operating points p ( 1 ) , p ( 2 ) , p ( 3 ) , p ( 4 ) , and p ( 5 ) . p ( j ) Test Point [ K P K D ] System PolesNumber of Unstable Poles in a Specific AreaStability of the Area [−100 5] 1 × 10 2 * −2.4885 + 0.0000i−0.6406 + 2.0984i−0.6406 − 2.0984i0.3144 + 0.0000i1unstable area P ( 1 ) [ 100 5 ] 1 × 10 2 *−2.4722 + 0.0000i−0.4420 + 2.0257i−0.4420 − 2.0257i−0.0991 + 0.0000i0stable area [ 1000 5 ] 1 × 10 2 *0.3379 + 2.2962i0.3379 − 2.2962i−2.2768 + 0.0000i−1.8542 + 0.0000i2unstable area [ − 100 5 ] 1 × 10 2 *−0.8189 + 1.9817i−0.8189 − 1.9817i−2.0321 + 0.0000i0.3107 + 0.0000i1unstable area P ( 2 ) [ 100 5 ] 1 × 10 2 *−2.0265 + 0.0000i−0.6212 + 1.8939i−0.6212 − 1.8939i−0.0901 + 0.0000i0stable area [ 1000 5 ] 1 × 10 2 *0.2069 + 2.1753i0.2069 − 2.1753i−1.8864 + 0.1139i−1.8864 − 0.1139i2unstable area [ − 100 5 ] 1 × 10 2 *−1.0710 + 1.8602i−1.0710 − 1.8602i−1.5108 + 0.0000i0.3105 + 0.0000i1unstable area P ( 3 ) [ 100 5 ] 1 × 10 2 *−0.8812 + 1.7439i−0.8812 − 1.7439i−1.5089 + 0.0000i−0.0710 + 0.0000i0stable area [ 1000 5 ] 1 × 10 2 *0.0458 + 2.0269i0.0458 − 2.0269i−1.9027 + 0.0000i−1.5311 + 0.0000i2unstable area [ − 100 5 ] 1 × 10 2 *−1.9547 + 1.8215i−1.9547 − 1.8215i0.3581 + 0.0000i−0.2211 + 0.0000i1unstable area P ( 4 ) [ 2000 5 ] 1 × 10 2 *−3.2066 + 0.0000i−0.2722 + 2.1504i−0.2722 − 2.1504i−0.0213 + 0.0000i0stable area [ 2000 − 5 ] 1 × 10 2 *−4.7157 + 0.0000i0.4825 + 1.7053i0.4825 − 1.7053i−0.0217 + 0.0000i2unstable area [ 1000 − 5 ] 1 × 10 2 *−4.4271 + 0.0000i0.3200 + 1.2076i0.3200 − 1.2076i0.0146 + 0.0000i3unstable area [ − 100 5 ] 1 × 10 2 *−1.6553 + 1.7938i−1.6553 − 1.7938i−0.5373 + 0.0000i0.3268 + 0.0000i1unstable area P ( 5 ) [ 200 5 ] 1 × 10 2 *−1.4543 + 1.5878i−1.4543 − 1.5878i−0.4857 + 0.0000i−0.1269 + 0.0000i0stable area [ 200 − 5 ] −377.1251−60.766567.240418.53702unstable area Table 3. For the given system with alternative values of operating points, the performance varies in terms of overshot, rising time, settling time, peak time, and steady-state error. Table 3. For the given system with alternative values of operating points, the performance varies in terms of overshot, rising time, settling time, peak time, and steady-state error. P ( j ) Lower Bounds [ K P K D ] Upper Bounds [ K P K D ] Optimal Values of [ K P K D ] Value of the CostFunction T r (sec) T P (sec) T s (sec) OS SSE P ( 1 ) [ 500 50 ] [ 2000 100 ] [ 568 50 ] 23.26 0.0018 0.0048 0.1512 58.7 0.1085 P ( 2 ) [ 500 50 ] [ 2000 100 ] [ 716 50 ] 26.56 0.0018 0.0048 0.0980 63.1 0.0805 P ( 3 ) [ 500 50 ] [ 2000 100 ] [ 867 50 ] 30.77 0.0019 0.0049 0.0717 56.0 0.0846 P ( 4 ) [ 4000 50 ] [ 8000 100 ] [ 6934 66 ] 23.22 0.0020 0.0048 0.4870 27.2 0.2142 P ( 5 ) [ 500 50 ] [ 4000 100 ] [ 1012 50 ] 47.46 0.0230 0.0053 0.1650 28.8 0.1455 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. © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Share and Cite MDPI and ACS Style
Almobaied, M.; Al-Nahhal, H.S.; Arrieta, O.; Vilanova, R. Design a Robust Proportional-Derivative Gain-Scheduling Control for a Magnetic Levitation System. Mathematics 2023, 11, 4040. https://doi.org/10.3390/math11194040
AMA StyleAlmobaied M, Al-Nahhal HS, Arrieta O, Vilanova R. Design a Robust Proportional-Derivative Gain-Scheduling Control for a Magnetic Levitation System. Mathematics. 2023; 11(19):4040. https://doi.org/10.3390/math11194040
Chicago/Turabian StyleAlmobaied, Moayed, Hassan S. Al-Nahhal, Orlando Arrieta, and Ramon Vilanova. 2023. "Design a Robust Proportional-Derivative Gain-Scheduling Control for a Magnetic Levitation System" Mathematics 11, no. 19: 4040. https://doi.org/10.3390/math11194040
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here. Article Metrics No No Article Access Statistics For more information on the journal statistics, click here. Multiple requests from the same IP address are counted as one view. Zoom | Orient | As Lines | As Sticks | As Cartoon | As Surface | Previous Scene | Next Scene Cite Export citation file: BibTeX | EndNote | RIS MDPI and ACS StyleAlmobaied, M.; Al-Nahhal, H.S.; Arrieta, O.; Vilanova, R. Design a Robust Proportional-Derivative Gain-Scheduling Control for a Magnetic Levitation System. Mathematics 2023, 11, 4040. https://doi.org/10.3390/math11194040
AMA StyleAlmobaied M, Al-Nahhal HS, Arrieta O, Vilanova R. Design a Robust Proportional-Derivative Gain-Scheduling Control for a Magnetic Levitation System. Mathematics. 2023; 11(19):4040. https://doi.org/10.3390/math11194040
Chicago/Turabian StyleAlmobaied, Moayed, Hassan S. Al-Nahhal, Orlando Arrieta, and Ramon Vilanova. 2023. "Design a Robust Proportional-Derivative Gain-Scheduling Control for a Magnetic Levitation System" Mathematics 11, no. 19: 4040. https://doi.org/10.3390/math11194040
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here. clear Mathematics, EISSN 2227-7390, Published by MDPI RSS Content Alert Further Information Article Processing Charges Pay an Invoice Open Access Policy Contact MDPI Jobs at MDPI Guidelines For Authors For Reviewers For Editors For Librarians For Publishers For Societies For Conference Organizers MDPI Initiatives Sciforum MDPI Books Preprints.org Scilit SciProfiles Encyclopedia JAMS Proceedings Series Follow MDPI LinkedIn Facebook Twitter © 1996-2024 MDPI (Basel, Switzerland) unless otherwise stated Disclaimer 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. Terms and Conditions Privacy Policy We use cookies on our website to ensure you get the best experience. Read more about our cookies here. Accept Share Link Copy clear Share https://www.mdpi.com/2492382 clear Back to TopTop