On the Discrete-Time Modeling of a DC-to-DC Power Converter and Control Design with Discrete-Time Sliding Modes

This work presents a novel discrete-time modeling of a boost dc-to-dc power converter by means of the symplectic Euler method. Then, on the basis of this model, a discrete-time sliding mode regulator is designed in order to force the power converter to behave as a dc-to-ac power converter. Simulation and experimental results are carried on, where the great performance of the proposed methodology is verified.


Introduction
Switched mode dc-to-dc converters as in [1,2] are mainly used as constant current sources for LED, LED flashlights, industry lighting, mobile phones, and automotive applications as in [3], among other commercial devices.With respect to dc-to-ac power conversion using the traditional switched dc-to-dc power converter topologies has recently attracted the attention of researchers.So far, various control techniques, either linear or nonlinear, to regulate these converters have been proposed, such as I/O feedback linearization [4], linear designs [5], sliding mode control [6,7], current-modecontrol [8], artificial neural networks [9], fuzzy logic control [10], and passivity-based control [11], among others.Since dcto-dc converters are absolutely nonlinear systems, for better performance, adopting nonlinear control methods is a good solution.
It is worth mentioning that all of the previously cited works are designed in the continuous-time domain and have been digitally implemented.The recent advancements in digital microprocessor technology have rendered cheaper, simpler, and more flexible digital implementation of control algorithms designed in the continuous-time setting.The main problem is the degradation of the expected performance, due to the fact that the designing of the control law does not take into account the sampled dynamics of the system.This problem motivates the following possible solutions: (1) the designing of the control law in the continuoustime setting and its implementation by means of zeroorder holder; (2) the exact sampling of the continuous-time system and the designing of the control law in the digital setting; (3) the use of much simpler discretization methods, such as those due to Euler (explicit or implicit), and the designing of the control law in the digital setting.
The first solution has the drawback of possible poor performance of the resulting sampled controller.For the second solution, the problem of exact sampling of continuoustime systems is not trivial.In fact, in general, a sampled closed representation of the sampled dynamics does not exist; whereas for linear systems, a sampled model in closed form can be easily obtained [12]; for nonlinear systems in general, the sampled data representations are given in the form of infinite series [13].Hence, in practice, one uses truncated models of desired approximation [7].The third solution [14] has the disadvantage that the accuracy of the Mathematical Problems in Engineering resulting approximate discrete-time system decreases as the sampling period increases.But, there is an alternative to these approaches that is only valid for dynamical systems modeled by Euler-Lagrange methodology, the use of variational integrators as the symplectic Euler method.The advantages of this method are explained and compared with common Euler (explicit or implicit) methods in [15], and the application in robotics has been satisfactory developed in [16] and electric motors in [17].
On the other hand, the sliding mode control [18] is a popular technique among control engineer practitioners due to the fact that introduces robustness to unknown bounded perturbations that belong to the control subspace; moreover, the residual dynamics under the sliding regime, that is, the sliding mode dynamics, can easily be stabilized with a proper choice of the sliding surface.For an insightful treatment of the sliding mode technique, we refer the reader to the book entitled "Sliding mode control in electro-mechanical systems, " [19].Electrical and electromechanical systems become vulnerable when output tracking signals present small oscillations of finite frequency known as chattering.The chattering problem is harmful because it leads to low control accuracy, high wear of moving mechanical parts, and high heat losses in power circuits.The chattering phenomenon can be caused by the deliberate use of classical sliding mode control technique.This control technique is characterized by a discontinuous control action with an ideal infinite frequency.When fast dynamics are neglected in the mathematical model such phenomenon can appear.Another situation responsible for chattering is due to implementation issues of the sliding mode control signal in digital devices operating with a finite sampling frequency, where the ideal infinite switching frequency of the control signal cannot be fully implemented, [20].
Therefore, in this work one proposes to derive a novel sampled model for the boost converter by means of a variational integrator.For that reason, a Lagrangian formulation for the dc-to-dc boost converter must be first established in order to apply the symplectic Euler method.Then, on the basis of the obtained sampled model a discrete-time sliding mode controller is designed in order to force the dc-to-dc boost converter to behave as a dc-to-ac converter.
The rest of this work is organized as follows: Section 2 describes the discrete Lagrangian mechanics for conservative and nonconservative systems.Section 3 deals with the discrete-time modeling of the boost power converter by means of its Lagrangian formulation.The control algorithm is designed in Section 4 based on a discrete-time sliding mode regulator technique.Section 5 deals with the presentation of simulation and real-time results.Finally, some comments conclude the work.

Lagrangian Mechanics and Variational Integrators
In this section, the basic principles of Lagrangian mechanics and discrete Lagrangian mechanics are reviewed from a variational point of view.

Continuous-Time Lagrangian Mechanics.
To define a mechanical system, let  be a smooth manifold, called the configuration space, and let  be its tangent bundle, called the phase space.The Lagrangian is a function  :  → R and is commonly defined as the kinetic energy  (usually, the only function of the velocity) minus the potential energy  of the system (usually, the only function of the state variable): The functional action is then introduced as the integral of  along a path () for time  ∈ [0, ]: ( With this definition, the main result of Lagrangian dynamics, Hamilton's principle, can be expressed quite simply: this variational principle states that the correct path of motion of a dynamical system is such that its action has a stationary value; that is, the integral along the correct path has the same value as within first-order infinitesimal perturbations.As an integral principle this description encompasses the entire motion of a system between two fixed times (0 and  in our setup).In more ways than one, this principle is very similar to a statement on the geometry of the path (): the action can be seen as the analog of a measure of curvature, and the path is such that this curvature is extremized (i.e., minimized or maximized).
In order to determine the path which optimizes the action we take a variation.A variation of the path  is written as  and can be thought of as an infinitesimal perturbation to the path at each point, with the important property that the perturbation is null at the endpoints of the path.Computing variations of the action induced by variations  of the path () results in where integration by parts is used in the last equality.When the endpoints of () are held fixed with respect to all variations () (i.e., (0) = () = 0), the most right term in the previous equation vanishes.Therefore, the condition of stationary action for arbitrary variations  with fixed endpoints stated in Hamilton's principle directly indicates that the remaining integrand in the previous equation must be zero for all time , yielding what is known as the Euler-Lagrange equation: For a given Lagrangian, this formula will give the equations of motion of the system.Finally to account for nonconservative forces or dissipation , the least action principle is modified as follows: ∫  0  ( () , q ())  + ∫  0  ( () , q ()) ⋅  = 0, (5) this is known as the Lagrange-d' Alembert principle.

Discrete Mechanics and Variational
Integrators.The idea of variational integrators was studied in [21,22], among others, and a general theory was developed over the subsequent decade (see [23] for a comprehensive survey).The main idea when deriving discrete geometric mechanics is to handle the variational nature of mechanics and to preserve this variational structure in the discrete setting.In fact, very few integrators have a variational nature: the simpler Euler methods as the explicit and the implicit ones are not variational, and both are characterized by demanding small sampling periods in order to perform satisfactorily.Instead of simply approximating the Euler-Lagrange equation to first (or higher) order, one can directly discretize the variational principle behind it.That is, if one designs a discrete equivalent of the Lagrangian, then, a discrete Euler-Lagrange equation can be easily derived from it by paralleling the derivations followed in the continuous case.In essence, good numerical methods will come from discrete analogs to the Euler-Lagrange equation, that is, equations that truly derive from a variational principle.
Suppose that we wish to discretize a mechanical system, whose continuous Lagrangian is  :  → R, on the time interval  ∈ [0, ].First, we partition the time interval by a finite number of points 0 =  0 , . . .,   =  where ,  ∈ N; likewise, the path  : [0, ] →  is replaced by a sequence of configurations  0 , . . .,   , where   ≈ (  ).Then the discrete Lagrangian is a function   :  ×  → R that approximates the Lagrangian  on each time interval [  ,  +1 ], with ℎ being the time interval between two samples ℎ =  +1 −   (chosen here to be constant for simplicity), so that Now, the right-hand side integral can be approximated through a generalized midpoint quadrature rule; that is, by the length of the interval times the value of the integrand is evaluated somewhere between   and  +1 and with q replaced by ( +1 −   )/ℎ: where  ∈ [0, 1] is parameter that generalizes the quadrature rule.For  = 1/2, the quadrature is second-order accurate, while any other value leads to linear accuracy.Also, one can approximate the right-hand side integral through a generalized trapezoidal quadrature rule: where again  ∈ [0,1] is parameter that generalizes the quadrature rule and for  = 1/2, the quadrature is secondorder accurate, meanwhile any other value leads to linear accuracy.
Remark 1.When the integrals are approximated by a generalized midpoint quadrature rule and  is fixed to 0 the variational integrator is known as the symplectic Euler method; when  is fixed to 1/2 the variational integrator is known as the midpoint Euler method.On the other hand, when the integrals are approximated by the generalized trapezoidal quadrature rule and  is fixed to 1/2 the variational integrator is known as the Störmer/Verlet method.
Given the discrete Lagrangian, the discrete functional action becomes simply a sum: Taking fixed-endpoint variations of this discrete action   , we obtain where  1  (resp.,  2 ) denotes the partial derivative with respect to the first (resp., second) arguments of   .Using the fixed endpoint condition  0 =   = 0, one gets Setting this variation to be equal to 0 and noting that each   is arbitrary, we arrive at the discrete Euler-Lagrange (DEL) equation Notice that this condition only involves three consecutive positions.Therefore, for two given successive positions   and  +1 , (12) defines  +2 .Thus, these equations of motion are actually the algorithm for an integrator.And since the DEL equation derives from the extremization of a discrete action, such an algorithm enforces the variational aspect of the motion numerically.
If we take initial conditions ( 0 ,  1 ), then the discrete Euler-Lagrange equations define a recursive rule for calculating the sequence {  }  =0 .Regarded in this way, they define a map (  ,  +1 ) → ( +1 ,  +2 ) which we can think of as a one-step integrator for the system defined by the continuous Euler-Lagrange equations.
We are primarily interested in discrete Lagrangian mechanics for deriving integrators.Any integrator which is the discrete Euler-Lagrange equation for some discrete Lagrangian is called a variational integrator.As we have seen previously, variational integrators can be implemented by taking two configurations  0 and  1 of the system, which should approximate  0 ≈ ( 0 ) and  1 ≈ ( 0 + ℎ) and then solve the discrete Euler-Lagrange equation ( 12) for  2 .This process can then be repeated to calculate an entire discrete trajectory.The map ( −1 ,   ) → (  ,  +1 ) defined by the discrete Euler-Lagrange equations is known as the discrete evolution map.
In mechanics, the initial conditions are typically specified as a position and a velocity or momentum rather than two positions, hence it is beneficial to write (12) in a positionmomentum form [24]. Therefore one defines the momentum at time   to be where the second equality holds due to (12).The positionmomentum form of the variational integrator discussed previously is then given by For (  ,   ) known, ( 14) is an implicit equation whose solution gives  +1 . +1 is then substituted in (15) in order to find  +1 .This provides an update rule in phase space.Now we consider nonconservative systems, those with forcing and/or dissipation; in such case, the discrete action can be modified by adding the nonconservative force term and using the discrete Lagrange-d' Alembert principle as in [23]: where  −  (  ,  +1 ) and  +  (  ,  +1 ) are discrete external forces acting, respectively, on the right of   and on the left of  +1 .In other words, can be seen as a two-point quadrature approximation of the continuous forcing term, that is, It follows that the forced discrete Euler-Lagrange equation is as follows: The forced discrete Euler-Lagrange equation can be expressed in a convenient, position-momentum form as follows: This variational treatment of energy decay, despite its simplicity, has also been proven to be superior to the usual time integration schemes that often add numerical viscosity to get stability [24].

Discrete-Time Modeling for a Boost Power Converter
In this section, a Lagrangian approach is used for deriving a novel discrete-time model for the boost power converter.The Lagrangian modeling is based on a suitable parametrization, in terms of the switch position parameter, of the EL functions describing each intervening system and subsequent application of the Lagrange formalism.

Discrete-Time Euler-Lagrange Modeling for a Boost Power
Converter.Let us define T  ( q  ) and V  (  ) as the kinetic and potential energies of the circuit, respectively.One can denote by F  ( q ) the Rayleigh dissipation function of the circuit and by V ,nc () the nonconservative potential function.These quantities are taken from [25] where they readily are found to be where   and   are the circulating charge in the inductor and the electrical charge stored in the output capacitor, respectively, with  = (  ,   )  ; moreover, q  and q  are the corresponding currents with q = ( q  , q  ), the constant parameters are  as the inductance,  as the capacitance, and  as the resistance.The vector  is defined as  = (, 0)  , where  is the dc power supply value.Figure 1 shows an electric diagram of a boost power converter circuit.Now, a nonconservative Lagrangian can be formulated as follows: where L  (, q ) is the conservative Lagrangian defined as usual: Following (7) one can obtain the discrete-time version of the conservative Lagrangian as Moreover the discrete external forces  −  (  ,  +1 ) and  +  (  ,  +1 ) appearing in (19) are as follows: where   = (,0)  and   are the sampled version of the external generalized forces of control input forces and dissipation forces, respectively.The vector of dissipation forces are obtained as where its sampled version is as follows: where again, the right-hand side integral can be approximated through a one-point quadrature, resulting as follows: After some algebraic steps and assigning  = 0 (i.e., the symplectic Euler method), the updated rule for position shown in (19) yields to where   = ( , ,  , )  and   = ( , ,  , )  ; moreover, Then, solving (28) for  +1 results in and finally in scalar form: where V  = 1 −   .Now, the updated rule for momentum equation shown in (28) Then, by replacing (30) in (32) it yields to The variable  , can be omitted since it does not enter in the other equations.In the continuous-time model [25] for the boost power converter the first state variable is defined as the derivative of the circulating charge   , that is, the input current represented as  1 = q  .The second state variable is defined as the output voltage capacitor and it is represented as  2 =   / where   is the electrical charge stored in the output capacitor.Therefore, one proposes the following novel approximated discrete-time model for the boost power converter: where x1, =  , , x2, =  , ,

Discrete-Time Control for a Boost Power Converter
There are two options for controlling dc-to-dc power converters, the indirect and direct methods.The first method is characterized by controlling the inductor current where its reference signal is selected such that a desired capacitor voltage is obtained.This has the advantage that the system is minimum phase with respect to the inductor current variable, but the problem consists of guaranteeing the convergence of the output capacitor voltage to its reference signal.In the second method, the output capacitor voltage is directly controlled, where common converter topologies as the boost, buck-boost, and cúk result to be nonminimum phase systems with the main difficulty of stabilizing the unstable inductor current variable [25].Now, the control problem is formulated.The control problem consists of directly controlling the output capacitor voltage variable x2, to track a desired reference signal x2, (  ) and at the same time to stabilize the unstable inductor current variable.

Sliding Mode Regulator for a Boost Power Converter. Let us define the following tracking error:
The reference signal is supposed to be generated by an autonomous exosystem given by with initial conditions  1,0 =  2,0 =  and  3,0 = , with the following output Then,  2, (  ) is a sinusoidal shape signal with amplitude of √ 2, a frequency equal to , and a bias value equal to ; moreover,   = ( 1, ,  2, ,  3, )  .Finally one determines that x2, (  ) = π2, (  ) =  2, (  ).
We define the steady state error as where x = (x 1, x2, )  and π (  ) = (π 1, π2, )  with π1, =  1, and π2, =  2, .Then, the dynamic equation for (38) can be obtained by taking one step ahead: Proposition 2. Consider the sliding mode function as and the sliding mode controller: where V eq, is the equivalent control calculated as usual from  +1 = 0 and the design parameter with  = 1,2; moreover,  1, is the right hand of the first equation in (39); then, the posed control problem is solved.
Proof.Condition  +1 = 0 results as follows: After making some simplifications one obtains where Equation (44) has been dealt in the work presented by [26], where the following equivalent control solution has been proposed: where Δ  =  2  − 4    .With this equivalent control, the sliding function is zeroed in one or more steps or tends asymptotically to zero.When 0 > V eq, > 1, the control action is saturated (41).In this case, for analysis purposes, the saturated control can be modified as follows: with V 0 = 1/2; moreover, after some algebraic steps, the sliding function dynamics (44) can be represented in the following form: In order to reveal the structure of the sliding function dynamics, let us represent it as follows: and the equivalent control in the case of Δ  ≥ 0 is rewritten as follows:  Substituting (50) in (47), and (47) in (49) yields to and by taking the absolute values At this point, one proposes the following candidate Lyapunov function: where the increment of the Lyapunov function results as follows: and by using condition (51), the increment of the Lyapunov function reduces to Δ < 0. Therefore, the sliding function (40) and the equivalent control (46) are decreasing monotonically, and after a finite number of steps, 0 < V eq, < 1 is achieved.Discrete-time sliding mode will take place from the next sampling point onwards.
Remark 3. When V eq, = −  /2  or V eq, = −  /  is chosen in (46), a similar Lyapunov analysis can be carried on, where the sliding function (40) can be demonstrated to decrease monotonically.
After the sliding mode occurs, one has  2, = − 1  1, (see (40)), and considering only the linear part of the first equation of (39) along with (46) as in [19], then, the motion of the linearized closed-loop system (sliding mode motion) will be governed by with  1,, as a function of higher order terms that vanish at the origin with their first derivative; then, one can assign a desired eigenvalue and solve for  1 , that is, With this choice it is easy to realize that the linear approximation of  1, in the sliding regime (56) is asymptotically stable, and under the property of center manifold, one has a locally attractive manifold on (38), that is,  1, → 0 ⇒ x1, = x1, and  2, → 0 ⇒ x2, = x2, as  → ∞.
Remark 4. Due to the complexity of the equivalent control (46), the sliding mode dynamics (56) must be determined for the three conditions presented in (46), where the design parameter  1 must stabilize the sliding mode dynamics for all cases.

Reference Signal Calculation for the Inductor Current.
The smooth mappings π1, (  ) :  0 → R and π2, (  ) :  0 → R (where  0 is an open neighborhood of   = 0), with π1,0 (0) = 0 and π2,0 (0) = 0, are such that the pair ( 1, (  ), π2, (  )) is the unique solution to the following difference equations: From ( 37) and ( 59) one can determine π2, (  ) = ( 1, +  3, ).Then, one can calculate   from (57) as follows: and substituting it in (58) yields to The finding of a solution to this difference equation results in a difficult task that can be solved by proposing an approximated solution as in [27].Thus, one proposes the following polynomial as an approximated solution for π1, (  ):   (63) Note that the bias value for π1, (  ) is the monomial with  18 as a coefficient, that is,  2 /().This fact was used when linearizing the sliding mode dynamics.

Simulations and Experimental Results
5.1.Open-Loop Simulations.Open-loop simulations were carried out in order to verify the performance of the obtained discrete-time The following two figures show the comparison of the continuous-time model of the boost power converter simulated with a sampling time of 0.0001 s and the discrete-time model with a sampling time of 0.001 s.The switch input value was set to  = 0. Figure 2 shows the performance of the inductor current, and Figure 3 shows the performance of the output capacitor voltage.In both figures can be appreciated that the discrete-time model it can accurately track the corresponding continuous-time signals.

Comparison of Two Control
Schemes.The designed controller in [28] (see the appendix) is discretized by means of the explicit Euler method and compared with the controller here designed.It is worth mentioning that the discretization of continuous-time controllers via this method is a common practice among control engineer practitioners.Figure 4 shows the performance of the output capacitor voltage with a sample of 10 s.In this figure one observes that both controllers have a good tracking performance.Figure 5 shows the simulation with a sample time of 50 s.In this figure one can observe that the proposed controller still performs well (Figure 5(b)); meanwhile, the sampled controller (Figure 5(a)) shows a difficulty for tracking the reference signal at lower peaks.Finally, Figure 6 shows simulation results where the sample time is incremented to 100 s.Here, one can appreciate the difference between the sampling of continuous-time controllers and the discrete-time control law design based on a good sampled model.There is a slight deviation in the output voltage capacitor for the proposed controller as can be seen in Figure 6(c), with the sampled controller, the output voltage capacitor signal is distorted (see Figure 6(a)).This is the intention of this simulation exercise, to show the degradation of the discretized controller with larger sampling periods and at the same time to show that the proposed controller can still perform better even with the larger sampling periods.It is worth mentioning that these discretetime control laws were simulated with the continuous-time plant model.The design parameter  1 in the sliding function (40) is fixed to 1 by means of a trial and error method.In the following, simulation results for the tracking of a sinusoidal signal with frequency of  = 1 rad/s will be shown.Figure 7 shows the sliding function where the finite-time convergence can be appreciated.

Simulations of the
Figure 8 shows the minimum phase behavior of the inductor current that is commonly known to be a nonminimum phase variable.
The output tracking result of the voltage capacitor can be appreciated in Figure 9. Now, simulation results for the tracking of a sinusoidal signal with a frequency value of  = 5 rad/s will be shown.Figure 10 shows the sliding function where once again the finite-time convergence can be appreciated.Figures 11 and 12 show the good tracking performance of the inductor current and the voltage capacitor, respectively.

Real-Time Experimentation.
The experimental setup is shown in Figure 13.Here, the VARIAC is a three-phase variable transformer fed from a three-phase voltage source.By rotating the knob of the VARIAC, the amplitude of the three-phase voltage source is regulated.These voltages are fed to the power module (Semikron) that incorporates a three-phase rectifier and a single switching transistor that is connected to the boost power converter elements; in that way, a regulated dc input voltage is available for the boost power converter.The control algorithm and PWM generation are programmed in Simulink and implemented with a DSP board  (dSPACE 1104).This board comes along with a library that easily incorporates with Simulink.Analog-to-digital converters included in the DSP board acquire the signals from the inductor current and capacitor voltage.Once the DSP executes the control algorithm in each sampling step, it generates one digital signal for switching the insulated-gate bipolar transistor (IGBT).This digital signal is TTL level and it is converted to a CMOS level of 15 V.This voltage level is the required one for switching on the IGBT.A block diagram of this setup is shown in Figure 14.
The power converter has the same parameter values of the digitally simulated boost circuit in Section 5.2.The inductor current and the output voltage capacitor are measured by means of hall-type sensors, as the HX 10-P and LV 25-P, respectively; both are manufactured by LEM.A first-order Butterworth low-pass filter having an 50 rad/s edge frequency was used for filtering the inductor current in order to attenuate the measurement noise.The sampling period has been fixed at ℎ = 70 s.The dc input voltage  was open loop (the switch position V  was set to 1) calibrated to a value of 20 V; then, the boost power converter was closed loop in order that the output voltage capacitor is forced to track a sinusoidal signal with an amplitude of 10 V and bias of 40 V.The real-time results for a frequency value of 1 rad/s are shown in Figures 15 and 16 and for a value of 5 rad/s are shown in Figures 17 and 18.In general, a good output tracking performance in both cases can be appreciated.The tracking of the inductor current shows in both cases big peak values of short duration that are admissible; moreover the first-order filter introduces the appreciated delay.
Remark 5.The tracking of the inductor is not so accurate due to the effect of the PWM switching which is correct since the system is not regulating the inductor current but the capacitor voltage.
In order to better evaluate the control performance, we consider the precision error   and the chattering effect  ℎ .The former is defined as the relative error of the output variable and it is calculated as the difference between the reference value and the average steady-state control output, divided by the reference value where   is the imposed reference and  mc is the average of the output controlled variable in steady state.The amplitude of the signal characterizes the latter, and it is given by the difference between the amplitude of the output variable and the average in steady state of the reference signal and divided by the latter signal where  mc is the maximum amplitude of the deviations of controlled signals with respect to its average value  mc .Table 1 summarizes the evaluation parameters for the two frequencies.One can observe that the proposed controller in closed loop yields to small percentiles in precision error and in chattering as well.
Remark 6.The chattering problem can be significantly reduced with the incorporation of soft computing techniques.
In particular, sliding mode control when combined with artificial neural networks can be useful to estimate the equivalent   control and magnitude of the switching control component using the conventional gradient learning algorithm.On the other hand, instead of using the well-defined sliding mode control functions, the control action can be derived from a fuzzy logic-based control table, compounded by fuzzy membership functions.In such a way, the convergence on the switching manifold and overshooting is softened, alleviating in this way the chattering problem.More directions can be found in [29].

Conclusions
A discrete-time sliding mode control scheme has been designed for the dc-to-ac boost power conversion problem.The contributions of this work are the following.
(i) A novel discrete-time model for the boost power converter circuit was obtained by means of a variational integrator scheme.(ii) The discrete-time-sliding-mode-based controller can stabilize the nonminimum phase system with a proper choice of a sliding surface, resulting in the output tracking of the capacitor voltage and the stabilization of the inductor current.(iii) The reference signal for the inductor current was proposed as a polynomial of order 3 that solves with a good accuracy the corresponding solution of the reference difference equation.Moreover, simulation and experimental results validate the good performance of the proposed dc-to-ac boost power conversion scheme.Some interesting issues, such as the robustness of the controller with respect to parameter variations, determination of the sampling period, are currently under study.

Figure 4 :
Figure 4: Comparison of discrete-time controllers with a sample time of 10 s.Output capacitor voltage (continuous line), reference signal (dashed line).(a) Output capacitor voltage obtained with sampled controller.(c) Output capacitor voltage obtained with proposed controller.(b) and (d) are zoom views of same variables on the left.

Figure 5 :
Figure 5: Comparison of discrete-time controllers with a sample time of 50 s.Output capacitor voltage (continuous line), reference signal (dashed line).(a) Output capacitor voltage obtained with sampled controller.(c) Output capacitor voltage obtained with proposed controller.(b) and (d) are zoom views of same variables on the left.

Figure 6 :
Figure 6: Comparison of discrete-time controllers with a sample time of 100 s.Output capacitor voltage (continuous line), reference signal (dashed line).(a) Output capacitor voltage obtained with sampled controller.(c) Output capacitor voltage obtained with proposed controller.(b) and (d) are zoom views of same variables on the left.
Proposed Control Algorithm.Closedloop simulations were carried out in order to verify the performance of the proposed controller.The initial conditions for the exosystem (36) are set to  = 7.071,  = 60 in order to generate a sinusoidal shape signal with amplitude value of 10 and a bias value of 60.The parameters of the boost converter are  = 0.098 H,  = 10 F,  = 241 Ω, and  = 40 V.

Figure 12 :
Figure 12: Output tracking of the voltage capacitor with  = 5 rad/s, reference signal (dashed line), real signal (solid line).

Figure 14 :
Figure 14: Block diagram of the experimental setup.

Figure 16 :
Figure 16: Output tracking of the voltage capacitor with  = 1 rad/s, reference signal (dashed line), real-time signal (solid line).
1, ≈  1 (ℎ) is the sampled version of the input current, and  2, ≈  2 (ℎ) is the sampled version of the output voltage capacitor.