Position Control of a Maglev System Fed by a DC/DC Buck Power Electronic Converter

In this paper, we solve the problem of position regulation in a magnetic levitation system that is fed by a DC/DC Buck power electronic converter as a power ampliﬁer. We present a formal asymptotic stability proof. Although this result is local, the merit of our proposal relies on the fact that this is the ﬁrst time that such a control problem is solved for a magnetic levitation system, a nonlinear electromechanical plant. In this respect, we stress that most works in the literature on control of electromechanical systems actuated by power electronic converters are devoted to control brushed DC motors which are well known to have a linear model. Furthermore, despite the plant that we control in the present paper is complex, our control law is simple. It is composed by four nested loops driven by one sliding mode controller, two proportional-integral controllers, and a nonlinear proportional-integral-derivative position controller. Each one of these loops is devoted to control each one of the subsystems that compose the plant: electric current through the converter inductor, voltage at the converter capacitor, electric current through the electromagnet, and position of the ball. Thus, our proposal is consistent with the simple and intuitive idea of controlling each subsystem of the plant in order to render robust the control scheme. We stress that such a solution is complicated to derive using other control approaches such as diﬀerential ﬂatness or backstepping. In this respect, our proposal relies on a novel passivity-based approach which, by exploiting the natural energy exchange between the mechanical and electrical dynamics, renders possible the design of a control scheme with the above cited features.


Introduction
One common technique that is used to supply power to electromechanical systems is pulse width modulation (PWM). However, the hard commutation that is intrinsic to PWM stresses the electromechanical system inducing abrupt changes in its dynamics which are observed as sudden variations in voltages and electric currents [1]. One manner to avoid this situation is the employment of DC/DC power electronic converters. Since these devices have embedded capacitors and inductors, they provide smooth voltages and electric currents, diminishing the effects of hard commutation in PWM-based power amplifiers. e mathematical models of some DC/DC power electronic converter-DC motor systems were proposed for the first time in [2]. Since then, many works have been reported on control of several DC/DC power electronic converter topologies and DC motors [3][4][5][6][7][8][9][10][11][12]. Among the proposed control techniques are differential flatness, proportional-integral (PI) control, generalized PI control, passivity, adaptive control, PI fuzzy control, LQR (linearquadratic regulator) control, backstepping, and hierarchical control.
e control problems that have been solved are unidirectional velocity regulation and tracking, velocity and torque control focusing on electrical transients, smooth velocity starters, and active disturbance rejection. In recent works [13][14][15], the introduction of an inverter between the DC/DC power electronic converter and the DC motor has rendered possible the bidirectional control of velocity. e approach in [16], to control the DC/DC Buck power electronic converter-DC motor system, was inspired in part by [17][18][19]. Control scheme in [16] has the advantage of including a PI loop to control voltage at converter capacitor, a PI loop to control motor armature's current, and an external PI loop to regulate motor velocity. Hence, the main components of the successful strategies employed in industry to control electromechanical systems are included in the proposal of [16]. Moreover, another internal loop is devoted to control electric current through converter inductance. is loop is driven by a sliding mode control, a common strategy for control of power electronic converters in practice. e approach is proven in experiments to be robust with respect to parametric uncertainties and external disturbances.
On the contrary, magnetic levitation systems are commonly used as benchmark problems to test novel control approaches. Among the proposed control techniques, the passivity-based approaches presented in [20][21][22] have been welcomed in the control community. In particular, the solution presented in [20] is interesting because it possesses a classical proportional-integral-derivative (PID) controller to cope with the mechanical part of the system. However, since the design is performed in terms of magnetic flux, instead of electric current, efforts are oriented to avoid the implementation of any internal loop to cope with the electrical dynamics.
is is because of the complications arising from magnetic flux measurements. In this respect, we stress that experimental results have been reported in the literature showing that such internal loop is necessary to improve performance in practice, see [23], for instance. e novel control technique known as immersion and invariance (I&I) has been employed in [24,25] to control magnetic levitation systems. Novelty in those applications is that a (small) parasitic capacitance is considered to be present at terminals of the electromagnet. e main target to use I&I in such a control problem is to extend the application of any control law, say w, that has been designed when such a parasitic capacitor is not present. However, since this requires to feedback the time derivative of w, the online computation of an important number of additional complex terms are required.
In the present paper, we extend the work in [16] to control the ball position in a magnetic levitation system which is fed by a DC/DC Buck power electronic converter.
is implies that additional inductance and capacitance with considerable values are included in the electrical circuitry of the magnetic levitation system. Since a magnetic levitation system only requires unipolar voltage, such a power converter topology is adequate and any inverter is not required. We stress that a magnetic levitation system is a complex and nonlinear system. Hence, controlling for the first time and from a theoretical point of view, a plant with these features when it is fed by a DC/ DC Buck power electronic converter represents one important contribution of the present paper.
Despite the complex and nonlinear nature of the magnetic levitation system, our proposal is simple. It is composed by a PI loop to control voltage at the converter capacitor, a PI loop to control the electromagnet electric current, and an external PID loop to regulate the ball position. As in [16], an additional sliding modes' internal loop is employed to control electric current through the converter inductance. Formalizing this intuitively simple idea to control a complex plant is another important contribution of the present paper. e key for this is a novel passivity-based approach exploiting energy ideas, i.e., we take advantage from the natural energy exchange among the several subsystems to design the control law. is represents another contribution of the present paper.
is paper is organized as follows. In Section 2, we introduce the plant to be controlled and present its dynamical model. e passivity properties of the plant are described in Section 3 where we also give some insight on the rationale behind our approach. Our main result is presented in Section 4. In Section 5, we present a simulation study and, finally, some concluding remarks are given in Section 6.

Mathematical Model
e DC/DC Buck power electronic converter-Magnetic levitation system is depicted in Figure 1(a). e DC/DC Buck power converter is composed by a transistor Q, a diode D, an inductor L c , a capacitor C, and a resistance R c . Symbols i c and υ represent electric current through inductance L c and voltage at capacitor terminals C, respectively, whereas E stands for voltage of the DC power supply. e system input is u which only takes the discrete values 0, +1 { } representing the off and on states of transistor Q, see Figure 1(b). e magnetic levitation system consists of an electromagnet, with inductance L(y) and internal resistance R, and a ball with mass m, made in a ferromagnetic material, which receives an upwards magnetic force F from the electromagnet. is force must cancel the downwards ball weight mg in order to levitate the ball in space. Electromagnet is basically a ferromagnetic core with a conductor wire wound around it. e electric voltage υ is applied at the electromagnet terminals which force an electric current i to flow through the electromagnet winding and this current produces the attractive magnetic force F � (1/2)(dL(y)/dy)i 2 on the ball. Symbol λ represents the magnetic flux produced by electric current i within the electromagnet core. Ball position, measured from the bottom of the electromagnet to the top of ball is represented by y ≥ 0. We remark that inductance of electromagnet, L(y) > 0 for all y ≥ 0, depends on the ball position y in the form shown in Figure 2. In order to understand this, recall that the magnetic flux is given as λ � L(y)i. Suppose that i remains constant and the ball approaches to electromagnet, i.e., y decreases. is reduces both the air gap and the reluctance. Hence, λ increases. Since λ � L(y)i and i remains constant, this means that L(y) must increase.
us, L(y) increases as y decreases. When y ⟶ ∞, it is obtained as the case when the ball is not 2 Complexity present, and in such a case, L(y) reaches a minimum positive value. is also means that dL(y)/dy < 0. We refer the reader to [23] for a detailed description of a magnetic levitation system as well as for precise instructions to construct one of them for experimental purposes. Furthermore, the complete procedure to obtain its dynamical model is presented and some experiments are provided to identify its parameters. Also, some controllers are designed and tested experimentally.
Important for our purposes is the following class of saturation functions. Definition 1. Given positive constants L * and M, with L * < M, a function σ: R ⟶ R: ς↦σ(ς) is said to be a strictly increasing linear saturation for (L * , M) if it is locally Lipschitz, strictly increasing, and satisfies [28] σ(ς) � ς, when |ς| ≤ L * ,

The Rationale behind Our Proposal
Consider the following slightly modified version of the mathematical model in (1)-(5): where G(y) � dP(y)/dy with P(y) is a positive semidefinite scalar function. e total energy stored in the system is given as follows: where the first term stands for electric energy stored in the capacitor of the Buck power converter, whereas the last three terms stand for the magnetic, kinetic, and potential energies stored in the electrical and the mechanical subsystems, respectively, of the magnetic levitation system. e time derivative of V e along the trajectories of system in (7) is given as follows: Notice that the cancellation of terms 1 2 represent (1) natural energy exchange between the electrical and the mechanical subsystems of the magnetic levitation system and (2) natural energy exchange between the capacitor and the electrical subsystem of the magnetic levitation system. Another cancellation of the terms involves ± G(y) _ y which represents the exchange between kinetic and potential energies in the magnetic levitation system. Hence, if we define the input i c and the output υ, then e expression in (11) proves that the model in (7) is output strictly passive [29], Definition 6.3.
In the present paper, we exploit these properties by proceeding as follows. First, we design u as a sliding mode controller to force i c to reach a desired function i d c . en, i c � i d c is employed as the control input for the sliding surface systems (2)-(4). In order to perform this step, we first obtain the error equation for this system by adding and subtracting some convenient terms (notice that these terms are not introduced using any control law), i.e., e equivalence of these expressions and those in (2)-(4) can be verified by reducing the redundant terms in the three expressions in (12). Defining V E1 (ε, I, _ y) � (1/2)Cε 2 + (1/2)L(y)I 2 + (1/2)m _ y 2 and using the three expressions in (12), we have that Hence, y and taking advantage from several natural cancellations, as in (10), we have where us, F s suitably shapes the potential energy of the mechanical subsystem to have a unique minimum at y � y * , whereas F d , i d cd , and υ d , represent the damping injection terms, as usual in standard passivity-based control [20].
In Section 4, we will show that several cross terms arising from the rectangular brackets in (15) do not cancel naturally.
is means that additional terms must be included in the control law if they are required to be cancelled. Hence, we prefer to dominate these terms instead of feeding back them in order to artificially cancel them. is allows us to design simpler control laws when compared to previous passivitybased approaches [20] where those terms must be online computed and fed back in order to be cancelled artificially.

Complexity
In this respect, we stress that it is recognized in the literature that increasing the number of online computations deteriorates performance because this increases numerical errors and the effects of noise. Moreover, we will also show that including PI and PID controllers (instead of proportional controllers as above) is straightforward. Notice that this feature is important to render robust the control scheme. e features described in this paragraph render novel and advantageous our passivity-based approach with respect to that in [20] where the natural cancellations shown in (10) are not exploited.

Main Result
Our main result is stated in the following proposition. (1)- (5) in a closed loop with the following controller:

Proposition 1. Consider the mathematical model in
where y * > 0 is a real constant standing for the desired position, h(y) � σ(y), and sat(z) � σ(z), where σ(·) is a strictly increasing linear saturation function for some (L * , M) (see Definition 1). Furthermore, it is also required that function σ(·) be continuously differentiable such that e closed-loop state evolution is assumed to be constrained to a subset D ⊂ R 7 , where dL(y) dy > c 1 , for some c 1 > 0 and c 2 > 0. Under these conditions there always exist constant scalars α, β, k p1 , k i1 , k p , k d , k i , α p , and α i , such that the closed-loop system has a unique equilibrium point which is asymptotically stable as long as At this equilibrium point, y � 0.

4.1.
Reaching the Sliding Surface. e time derivative of the positive definite and radially unbounded scalar function where [15] has been used, if (24) implies (23). From the sliding condition _ s c � 0, [1,22], we find that the equivalent control satisfies the following bound: which means that the sliding regime is possible. On the contrary, (24) ensure that the sliding surface s c � i c − i * c � 0 is reached, i.e., i c � i * c is reached. us, we only have to study the stability of dynamics (2)-(5) in closed loop with (17)- (20) when evaluated at i c � i * c .

Closed-Loop Dynamics on the Sliding Surface.
Using [16] in [2], and adding and subtracting the terms i * , where dL y * dy � dL(y) dy y�y * .
is means that the only equilibrium point of the closed- Notice that this closed-loop dynamics is autonomous because it can be written as _ ζ � f(ζ) for some nonlinear f(·) ∈ R 7 .

Stability Analysis.
e closed-loop dynamics (26)-(33) and (20) can be rewritten as follows: Notice that (35)-(38) are almost identical to the openloop dynamics in (7) if we replace y, _ y, i, υ, G, i c by y, _ y, i, Υ, G, I c . One important difference is that the resistances R c and R in (7) have been enlarged to (1/R c ) + k p1 and R + α p L(y) in (35) and (36), respectively. Moreover, we can see that suitable damping can be introduced thanks to term k d _ y in the definition of G. Another important difference is the three new equations in (38) which represent the integral terms of the PI electric current controller, the PI controller of voltage at the capacitor, and the PID position controller, which are intended to compensate for the effects of the gravity term mg. ese observations motivate the use of the following "energy" storage function for the closed-loop dynamics:

Complexity
We stress that function V m ( _ y, y, z − mg/k i ) defined in (40) is very similar to the function V(q, _ q, z − (k i ′ ) − 1 g(q d )) analyzed in [30].
us, conditions to ensure that V m ( _ y, y, z − mg/k i ) is positive definite and radially unbounded are the same as those introduced in [30]. For the sake of completeness of this work, these conditions are presented in Appendix A as (A.1), (A.2), and (A.5) and k p > 0, α > 0, and β > 0.
us, the function W( _ y, y, z − mg/k i , i, z 1 , e, ξ) qualifies as a Lyapunov function candidate because it is positive definite and radially unbounded if k i1 > 0 and α i > 0. e first two terms in W represent the electric energy stored in the converter's capacitor and "energy" stored in the integral term of the PI voltage controller. e third and fourth terms represent the magnetic energy stored in the electrical system and the "energy" stored in the integral term of the PI electric current controller. On the other hand, function V m includes the kinetic energy and the closed-loop "potential energy" as well as the "energy" stored in the integral of position through an integral of s(·). e cross terms αmh(y) _ y and αβms(z) _ y are required to provide _ W with negative quadratic terms in both h(y) and s(z). In this respect, it is easy to verify that d dt Since Υ depends on both h(y) and s(z), negative quadratic terms on both h(y) and s(z) are required to dominate some cross terms in both h(i) and i and s(z) and i.
is is the reason for the cross terms αmh(y) _ y and αβms (z) _ y (the quadratic term − (R + α p L(y))i 2 already exists).
e integral term αk d y 0 h(r)dr is intended to cancel an undesired cross term appearing in the time derivative of αmh(y) _ y. Notice that term ie, arising from the product i Υ above, cancels with term − ei, arising from (d/dt)(Ce 2 /2). Moreover, some third order terms involving i * appear from iΥ.
ese terms can be dominated by quadratic negative terms in _ y and i because |h(y)| and |s(z)| are bounded by finite constants. is is the reason to employ a PID position controller with saturated proportional and integral actions.
After some straightforward natural cancellations (i.e., not requiring to use additional terms in the control law to achieve them), which include the closed-loop equivalents of (10) (see Remark 6), we find that the time derivative of W along the trajectories of the closed-loop systems (35)-(38) is given as follows: Taking into account (19), (31), and Appendix B, it is found that _ W can be upper bounded as follows: x T � [| _ y|, |h(y)|, |s(z)|, |i|, |e|], where k d1 , k d2 , k d3 , k d4 , α p1 , α p2 , α p3 , α p4 , k p2 , and k p3 are positive constant scalars such that k d1 + k d2 + k d3 + k d4 � k d , α p1 + α p2 + α p3 + α p4 � α p , and k p2 + k p3 � k p1 . e entries of matrix Q are defined as follows: where k r � max L(y) 2F * |dL(y)/dy| 8 Complexity and constants k σ , k δ , k * , k t , and I * (0) are defined in Appendix B. Notice that the four leading principal minors of matrix Q can always be rendered positive definite by suitable selection of the controller gains k p1 , k d1 , k p , k i , and α p1 , and hence, λ min (Q) > 0. Furthermore, all of the expressions in the seven rows in (45) can be rendered negative, at least locally, using suitable gains k d2 , k d3 , k d4 , α p2 , α p3 , α p4 , k p2 , and k p3 . Hence, it is concluded that _ W ≤ 0 for all ζ ∈ D, where D is a subset of R 7 containing the origin. us, stability of the origin is concluded. Since the closed-loop system is autonomous, the LaSalle invariance principle is invoked to prove asymptotic stability. is completes the proof of Proposition 1.

Remark 1.
Notice that the ith leading principal minor of matrix Q can be rendered positive by choosing large enough Q ii entry of matrix Q and choosing small some constants such as α > 0 and β > 0. In this respect, it is important to stress that each one of the Q ii entries depend on a controller gain. us, once the ith leading principal minor is rendered positive by choosing a suitable controller gain, the (i + 1)th leading principal minor can be rendered positive by suitably enlarging the entry Q (i+1)(i+1) of matrix Q and so on. Moreover, these ideas suggest that an intuitive try-and-error tuning procedure can also be derived; choose larger controller gains in the Q ii entries as i is larger.
Remark 2. Notice that the system evolves by itself as long as the sliding surface is not reached. However, it is proven in standard books on sliding mode control that the system reaches the sliding surface in a finite time which depends on the initial conditions. Once the sliding surface is reached, the system evolves on the sliding surface starting from the state values that it has at the time where the sliding surface is reached. us, if initial conditions are close to the desired equilibrium point, then asymptotic convergence to the desired equilibrium point is ensured by Proposition 1.
On the other hand, limits in currents i c and i are imposed by the circuit parameters and the applied voltage. With the problem at hand, limits on these currents are imposed by the DC power supply E which is manipulated through the variable u representing the on-off state of transistor Q. In Section 5, we explain how the variables u and E affect the system performance. Finally, saturation of the internal PI's only might occur if the state variables are allowed to take large values. Since our stability result is local, this prohibits the state to take such large values. Figure 3, we present a block diagram of the control scheme in control scheme in Proposition 1. Notice that this controller is made up of three main loops: (1) a PI controller for voltage at the DC/DC Buck power converter output (at the capacitor terminals), (2) a PI controller for electric current through the electromagnet, and (3) a nonlinear PID (NPID) controller for the ball position. us, our proposal contains the fundamental components in industrial applications and, hence, it is expected to be robust with respect to parametric uncertainties and external disturbances. Furthermore, another internal loop is provided to control electric current through the inductance of the DC/ DC Buck power converter. is loop is driven by a sliding mode controller which constitutes one common technique to control power electronic devices in practice.

Remark 4.
It is stressed that we introduce the factor L(y) as a part of the proportional gain in (18) in order to ensure to be constant the term that is added to the integral term in (27).
is is a necessary step for the integral action of the PI controller of voltage at the capacitor. In this respect, the last two terms in (17) are included in order to cancel some terms arising in (26) because of the fact that − C _ υ must be added and subtracted to complete (26). (18) and the terms − k m e _ y 2 − k n e|i| in (17), for some positive constants k q , k f , k h , k m , and k n , to proceed as in [30] to dominate the positive terms in the five rows of (45). However, we have decided not to include the above cited terms because of several reasons. (1) e stability result would still remain to be local, as usual in magnetic levitation systems. (2) Including the above terms in (18) would result in additional complex terms that should be cancelled using the definition of i * c in (17). (3) In order to maintain the simple control law, we have decided not to proceed including the terms − k q i _ y 2 − k f i|i| − k h i|e| in (18) which renders useless to include − k m e _ y 2 − k n e|i| in (17). us, proposing a simple and robust control law performing well is one important objective of our proposal.

Remark 6.
e novel passivity-based approach that is employed in this paper has the following properties: (i) Several terms cancel naturally. is means that they cancel without requiring to compute and feedback them. is property is a direct consequence of the fact that the closed-loop dynamics (35)-(38) is almost identical to the open-loop model in (6).
Recall that the open-loop model was proven in Section 3 to be output strictly passive, and the existence of several natural cancellations is instrumental for this. is property is opposite to what happens in exact feedback linearization approaches where the plant undesired terms are online computed and fed back in order to force their cancellation.
ese natural term cancellations are very useful to obtain a simple control law. As stated above, other control approaches require to use additional terms in the control law in order to cancel these terms. (ii) A nested-loop passivity-based control approach is exploited in [20]. is means that the electric current error is first proven to converge exponentially to zero, and this allows to use this variable as a vanishing perturbation for the mechanical subsystem. is, however, requires the online computation of either the time derivative of the desired force or the time derivative of the desired electric current. Instead of that we use an approach which is similar to what was called in [20] passivity-based control with total energy shaping. Although the latter approach has been disregarded in [20] by arguing that it results in more complex controllers, we prove the opposite in the present paper. is is one important novelty of our approach. (iii) e previous features of our approach allow (1) to straightforwardly include PI internal loops and an external PID loop, which are important to improve the robustness properties of the control scheme and (2) to avoid the requirement on the exact knowledge of the electric resistance of the electromagnet, which is a parameter that changes during normal operation conditions.

Remark 7.
e present work is inspired by [16] in the sense that some PI and PID loops are included for both the electrical and the mechanical dynamics of the electromechanical system, and a sliding mode controller is devoted to control electric current through the inductor of the DC/DC Buck converter, see Figure 4. However, the control scheme in [16] is designed for unidirectional control of velocity in a brushed DC motor, whereas the controller in the present work contains several refinements that extend work in [16] to control position in a magnetic levitation system, a highly nonlinear system. ese refinements include (1) a clever selection of a nonlinear PID position controller, (2) a clever selection of internal PI controllers, and (3) proposing a suitable Lyapunov function for stability analysis. We stress that, aside from the sliding mode controller, the closed-loop system in [16] is linear, whereas both plant and controller are nonlinear in the present work.

Simulation Results
In this section, we present a numerical example to give some insight on the achievable performance when the controller in Proposition 1 is employed. To this aim, we use the numerical values of the magnetic levitation system that has been tested experimentally in [23]. In that work, the electromagnet inductance is modeled as follows: where k 0 � 36.3 × 10 − 3 H, k � 3.5 × 10 − 3 H, and a � 5.2 × 10 − 3 m, and the remaining parameters are R � 2.72 Ohm, m � 0.018 kg, and g � 9.81 m/s 2 . e practical range of input voltages is [0, +12] V and the range of electric current through the electromagnet is [0, +3] A. e Buck DC/DC power electronic system parameters were chosen as L c � 0.686 H, C � 114.4 × 10 − 6 F, and R c � 28.5 Ohm. We also use E � 12 V if 2 s < t < 3 s and E � 50 V, otherwise in order to study the response when disturbances appear in the DC power supply. ese parameters are similar to the experimental values reported in [16]. e controller gains were chosen to be k p � 8, k d � 1, k i � 2, α p � 470, α i � 1000, α � 64, β � 1, k p1 � 6000, and k i1 � 18000. Inspired by [28], we used the saturation function: where In Figure 5, we can see that the actual ball position y reaches its desired value y * in the steady state. Notice that the settling time is about 0.5 s. We observe that the position response is very damped. In this respect, we have performed several additional simulations which make us to conclude that the term (1 + αβ(k d /k i )) _ y appearing in (20) is responsible for such a damped response.
In Figures 6 and 7, we verify that electric current through the electromagnet i and voltage at the electromagnet terminals υ remain within [0, +3] A and [0, +12] V, respectively, the actual ranges of values reported experimentally in [23].
In Figure 8, we present electric current through the converter inductor which also remains within the range [0, +3] A. Finally, in Figure 9, we verify that evolution of all the above signals is achieved by applying, at the transistor input, a switching signal only taking the discrete values 1 or 0. Notice that, despite this hard switching signal, electric current through the electromagnet is smooth enough, see Figure 6. Moreover, this is achieved despite voltage at the electromagnet terminals has a small ripple, see Figure 7.
We also observe in Figure 5 that the step changes in the DC power supply, i.e., E, appearing at t � 2 s and t � 3 s have an almost imperceptible effect in the ball position. Moreover, the effects of these changes are more noticeable in Figures 7  and 8 as a simple change in ripple of both voltage at the electromagnet terminals and electric current through the converter inductor. e effects of these step changes are very small thanks to the employment of a sliding mode controller for electric current through the converter inductor and this is the reason why we employ sliding mode control in our approach.
Finally, we stress that the control signal u is constrained to only take the discrete values 0 or 1 and all of the system evolution must be controlled while u takes these discrete values. However, the designer has an additional degree-of-    Complexity freedom in this respect; a larger value for E can be chosen if the system evolution requires more control effort. We have arrived to this conclusion during the simulations that we have performed. For instance, we can observe in Figure 9 three clear zones at t � 1 s, t � 4 s, and t � 5.5 s which are evidences that control effort saturates at its extremum values, i.e., either 1 or 0, when step changes in the desired position are commanded. We also have observed that the ball escapes if larger reference changes are commanded. is is the reason why we have used E � 50 V for most of the time in the above simulations. Now, we present some simulation results when using the following classical control scheme: It is assumed that the DC/DC Buck power electronic converter is not present and a static PWM-based power amplifier is employed. Notice that position error is given as y − y * for the classical PID controller.
is is usual in control of magnetic levitation systems. A feature that is required because of the negative gain of the system which, in turn, is produced by the fact that (dL(y)/dy) < 0. e controller gains were chosen to be k pi � 4.7, k ii � 100, k d � 0.0277, k p � 0.866, and k i � 1.082. is results in two real dominant closed-loop poles located at s � − 21.5 and s � − 2.11. A third real closed-loop pole is at s � − 140 and two very fast complex conjugate poles, due to the PI control of the electrical dynamics, located at s � − 19 ± 500j.
Our intention is not to compare results obtained with the controller in Proposition 1 and result with the classical control scheme since such a comparison would be unfair. We just want to point out on some features of the classical control scheme. e simulation results are shown in Figures 10-12. e desired position y * is given, again, as in (50).
Notice that the position response has a very short rise time but a large settling time. is is produced by initial fast changes on both the applied voltage and the resulting electric current, which tend slowly to constant values afterwards. Notice that position response exhibits a large overshoot despite the dominant closed-loop poles being real. Recall that the fast complex conjugate poles due to the electrical dynamics cannot produce such a slow overshoot. As it is clearly explained in [31], reason for such a large overshoot is the open-loop unstable pole of the magnetic levitation system which is responsible of its open-loop instability. It is also demonstrated in [31] that this feature cannot be avoided when using the classical control scheme in (51) and (52).
Since classical control schemes rely on linear approximations of the plant to control, it is reasonable to wonder whether some advantages could be obtained when proposing control schemes that take into account more information of the plant nonlinear dynamical model. is is the intention of nonlinear control schemes as the one presented in Proposition 1 aside from taking into account some additional dynamics as that of the DC/DC Buck power electronic converter.
Notice that the expressions in (18)- (20) are the nonlinear versions of (51) and (52), i.e., they represent a PI electric current loop and a PID position loop. On the contrary, (16) and (17), represent a sliding mode control loop for electric current at the converter inductor and a nonlinear PI loop for voltage at the converter capacitor. is means that (16) and (17) are additional control loops that are included to cope with system components that are not present in standard magnetic levitation systems.
is is what we mean when stating that the control scheme in Proposition 1 is simple if we take into account the complexity of the nonlinear model of the plant.

Conclusions
In this paper, we have solved the position regulation problem in a magnetic levitation system when it is fed by a DC/DC Buck power electronic converter as power amplifier.  12 Complexity Contrary to previous works in the literature, which are concerned with control of DC motors, we have solved this problem for the first time for a nonlinear electromechanical system. is fact represents the merit of our proposal despite the local nature of the results. In this respect, we stress that the region of attraction might be enlarged by including additional nonlinear terms to the control law. However, we have decided not to do this because we are interested in presenting a simple control law in order to render it free of numerical errors and noise amplification.
Our proposal relies on a novel passivity-based approach which exploits the natural energy exchange existing among the mechanical and electrical subsystems that compose the plant. is allows several nonlinear terms to naturally cancel, i.e., without requiring to include additional terms in the control law, which renders simple the control law. Moreover, contrary to previous well known works in the literature, our approach does not rely on proving exponential convergence to zero of the electrical dynamics. is is good news since, otherwise, we would require to online compute and feedback both the time derivative of the desired electric current in the electromagnet and the time derivative of the desired voltage at the electromagnet terminals. Such online computations would render the control law complex and sensitive to numerical errors and noise amplification.
Finally, although the control law is simple, the stability proof may require much attention from the reader. is might be seen as a disadvantage of the approach but it is the authors' belief that it is the role of the control community to cope with mathematical analysis and other related theoretical complications. e important practical thing is that the resulting control scheme is intuitively simple to understand and to implement: a multiloop scheme driven by PI and PID controllers provided with explicit tuning guide lines. We must also say that our approach only considers the regulation task.