Next Article in Journal
Combination of the Single-Valued Neutrosophic Fuzzy Set and the Soft Set with Applications in Decision-Making
Next Article in Special Issue
Analysis of Chaotic Response of Frenkel-Kontorova-Tomlinson Model
Previous Article in Journal
Shortest Path Solution of Trapezoidal Fuzzy Neutrosophic Graph Based on Circle-Breaking Algorithm
Previous Article in Special Issue
Anatomical Considerations and Study of the Fractal Dimension around the Posterior Superior Alveolar Artery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrete-Inverse Optimal Control Applied to the Ball and Beam Dynamical System: A Passivity-Based Control Approach

by
Oscar Danilo Montoya
1,2,*,
Walter Gil-González
3 and
Carlos Ramírez-Vanegas
4
1
Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Carrera 7 No. 40B-53, Bogotá, D.C 11021, Colombia
2
Laboratorio Inteligente de Energía, Universidad Tecnológica de Bolívar, km 1 vía Turbaco, Cartagena 131001, Colombia
3
Grupo GIIEN, Facultad de Ingeniería, Institución Universitaria Pascual Bravo, Campus Robledo, Medellín 050036, Colombia
4
Facultad de Ciencias Básicas, Universidad Tecnológica de Pereira, Pereira 660003, Colombia
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(8), 1359; https://doi.org/10.3390/sym12081359
Submission received: 5 August 2020 / Accepted: 12 August 2020 / Published: 14 August 2020
(This article belongs to the Special Issue Advances in Nonlinear, Discrete, Continuous and Hamiltonian Systems)

Abstract

:
This express brief deals with the problem of the state variables regulation in the ball and beam system by applying the discrete-inverse optimal control approach. The ball and beam system model is defined by a set of four-order nonlinear differential equations that are discretized using the forward difference method. The main advantages of using the discrete-inverse optimal control to regulate state variables in dynamic systems are (i) the control input is an optimal signal as it guarantees the minimum of the Hamiltonian function, (ii) the control signal makes the dynamical system passive, and (iii) the control input ensures asymptotic stability in the sense of Lyapunov. Numerical simulations in the MATLAB environment allow demonstrating the effectiveness and robustness of the studied control design for state variables regulation with a wide gamma of dynamic behaviors as a function of the assigned control gains.

1. Introduction

The ball and beam dynamical system is a classical and well-known nonlinear dynamical system that attracts much attention in the control area [1]. The main challenge in this plant is to regulate all the state variables at the origin, taking into account that it corresponds to a fourth-order dynamical system ( n = 4 ) with only one control input ( m = 1 ) [2]. Additionally, the system contains strong nonlinearities such as products between variables and trigonometric functions [3]. Analyzing the ball and beam system from the control analysis point of view is interesting as its model can be used to understand transportation systems, communication, or power system dynamics [4]. In the literature, the problem of control in ball and beam systems has been addressed with linear and nonlinear techniques with some of these are presented as follows; proportional-integral control [5], feedback of state variables [6], linear matrix inequalities [7] exact feedback linearization [8], control Lyapunov functions [9,10], adaptive control design [2], passivity-based control [11], fuzzy logic [12,13], and artificial neural networks [14], among others.
The main contributions of our approach, different from previous works, can be summarized as follows [15].
The application of the discrete-inverse optimal control to regulate all the state variables of the ball and beam system guaranteeing passivity, stability, and optimality properties.
The numerical validation via simulations by working the discrete equivalent nonlinear model of the system without any special assumption on the open- or closed-loop dynamics.
The robustness and effectiveness of the discrete-inverse optimal control design when parametric variations affect the discrete dynamical model.
It is important to mention that after an exhaustive revision of the literature about the ball and beam dynamical system modeling and control, we identify that the discrete-inverse optimal control has not been studied in this system. This indicates that it is a clear opportunity for research that this article tries to fill.
The remainder of this document is organized as follows. Section 2 presents the general dynamical formulation for the ball and beam dynamical system on the continuous and discrete domains. Section 3 presents the general theory about discrete-inverse optimal control applied on nonlinear systems by highlighting its passivity, stability, and optimality properties. Section 4 shows all the numerical simulations and their corresponding analysis and discussion. Section 5 presents the main conclusions derived from this research study.

2. Dynamical Model and Discretization

The ball and beam system generates a fourth-order nonlinear set of ordinary nonlinear differential equations [2]. The main objective is to regulate all the state variables around the origin of coordinates guaranteeing stability in closed-loop. Figure 1 illustrates the physical form of the ball and beam system.
In this figure the variables and parameters have the following interpretation; θ and r are the beam angle and the ball position, respectively. Moreover, k 1 is the steady-state gain; τ is the time-constant, L beam is the length of the beam; and m and J b are the mass and moment of inertia of the ball, respectively. In addition, R is the radius of the ball, g is the acceleration due to gravity, r arm is the distance between screw and motor gear, and V m is the input of the system [2]. The dynamical formulation of the ball and beam system depicted in Figure 1 is presented below.
r ¨ = m r arm g R 2 L beam m R 2 + J b sin θ m J b R 2 + m r θ ˙ 2 ,
θ ¨ = 1 τ θ ˙ + k 1 τ v m ,
To reach the state space representation of the dynamical model (1), let us define k b b = m r arm g R 2 L beam m R 2 + J b , h = m J b R 2 + m , and x 1 , x 2 , x 3 , x 4 T = r , r ˙ , θ , θ ˙ T , which allows finding the following fourth-order nonlinear dynamical model,
x ˙ 1 = x 2 ,
x ˙ 2 = k b b sin x 3 h x 1 x 4 2 ,
x ˙ 3 = x 4 ,
x ˙ 4 = 1 τ x 4 + k 1 τ u ,
where u = v m corresponds to the control input. Note that the equilibrium point of the dynamical system (2) is x 1 🟉 , x 2 🟉 , x 3 🟉 , x 4 🟉 T = r 0 , 0 , 0 , 0 T , being r 0 an arbitrary point between the extremes of the beam; nevertheless, the unique physical solution possible is when r 0 = 0 , due to the gravity force in the real system will make any other point unstable as can be seen in Figure 1.
To design the proposed controller it is needed to represent the continuous dynamical system (2) into a discrete equivalent. For doing so, the classical forward difference is applied [16], i.e.,
x ˙ = f x x k + 1 = Δ k f x k + x k ,
where subscript k is the current sample and Δ k represents the discretization time. If we apply the forward difference in the set of Equations (2), the the following discrete system is reached. l
x 1 k + 1 = Δ k x 2 k + x 1 k ,
x 2 k + 1 = Δ k k b b sin x 3 k h x 1 k x 4 k 2 + x 2 k ,
x 3 k + 1 = Δ k x 4 k + x 3 k ,
x 4 k + 1 = Δ k 1 τ x 4 k + k 1 τ u k + x 4 k ,

3. Inverse Optimal Control Design

In this section three main aspects of the inverse optimal control design will be explored for general nonlinear discrete systems [15,17]. For doing so, let us define the general structure of the system under analysis as follows [18,19].
Definition 1.
A nonlinear dynamical system in the discrete domain with the form,
x k + 1 = f x k + g x k u k ,
y k = h x k + j x k u k ,
fulfills passivity properties, is globally asymptotically stable, and also there is a control law with the form u k = y k , such that a functional cost function is minimized, i.e., u k is an optimal control law. Note that in (4), y k is the output of the system and h x k and j x k take the followings structures,
h x k = g T x k P f x k ,
j x k = 1 2 g T x k P g x k
being Q a symmetry positive definite matrix, i.e., P = P T 0 .
To demonstrate each one of the properties presented in Definition 1, let us consider a candidate Lyapunov function with a quadratic form as follows,
V x k = 1 2 x k T P x k ,
which is positive definite for all x k 0 and zero only for x k = 0 . In addition, let us define a general form for the control input u k as follows,
u k = β x k + v k ,
where v k is the new input and β x k can be defined as presented below,
β x k = I + j ( x k ) 1 h ( x k ) ,
being I an identity matrix with appropriate dimensions.
Definition 2.
The dynamical system (4) exhibits passivity properties if there is a matrix P such that the following inequality is held.
f x k + g x k β x k T P f x k + g x k β x k x k T P x k ,

3.1. Passivity

To demonstrate passivity properties in the the dynamical discrete system consider Lemma 1 as presented below [20,21,22].
Theorem 1.
The dynamical system in (4) is a feedback passive system for the output y ˜ k . The control input is defined as (7), where y ˜ k takes the following form,
y ˜ k = h ˜ x k + j x k v k .
where
h ˜ x k = g T x k P f ˜ x k ,
f ˜ x k = f x k + g x k β x k
Proof. 
To proof the feedback passivity properties of the dynamical system (4), consider the variation of the Lyapunov function for the current and the future states as follows,
Δ V = V x k + 1 V x k .
Note that using (4) and (7) in (12), we have
Δ V = 1 2 f x k + g x k β x k T P f x k + g x k β x k 1 2 x k T P x k + f x k + g x k β x k T P g x k v k + 1 2 v k T g T x k P g x k v k .
From (13), we can observe that
f x k + g x k β x k T P g x k v k = h ˜ T x k v k ,
v k T g T x k P g x k v k = 2 v k T j T x k v k .
Now, if we consider Definition 2 and expressions in (14) to be replaced in (13), then we have
Δ V y ˜ k T v k ,
which confirms that the discrete system is passive from the output y ˜ k to the new input v k and the proof about passivity is completed. □

3.2. Stability

To demonstrate passivity properties the stability properties in the sense of Lyapunov for closed-loop operation, let us consider the following Lemma.
Theorem 2.
The system (4) is asymptotically stable in the sense of Lyapunov with the control input (7) if v k is defined as
v k = y ˜ k = I + j ( x k ) 1 h ˜ ( x k ) .
Proof. 
To proof stability in the sense of Lyapunov, we can transform the dynamical system (4) with the control input (7) as an equivalent system with the following structure
x k + 1 = f ˜ x k + g x k w k .
Now, if we consider the difference between the current and the next step of the Lyapunov function defined in (12), we have
Δ V = 1 2 f ˜ x k + g x k w k T P f ˜ x k + g x k w k 1 2 x k T P x k = f ˜ T x k P g x k v k + 1 2 v k T g T x k P g x k v k + 1 2 f ˜ T x k P f ˜ x k x k T P x k .
From (18), we can note that
f ˜ T x k P g x k v k + 1 2 v k T g T x k P g x k v k = y ˜ k T v k ,
in addition, from Lemma 2, we know that v k = y ˜ k , which implies that in conjunction with (19), the expression (18) takes the following form
Δ V = 1 2 f ˜ T x k P f ˜ x k x k T P x k v k 2 < 0 ,
which allows to conclude that the system (4) is globally asymptotically stable in x k = 0 as the candidate Lyapunov function V ( x k ) = 1 2 x k T P x k is radially unbounded. This completes the proof. □

3.3. Optimality

Theorem 3.
The inverse control law (7) is considered optimal since it stabilizes the dynamical system as presented in Section 3.2, and it minimizes the following functional cost
F = k = 0 L x k , β x k ,
where L x k , β x k is the LaGrangian function of the system that can be written as
L x k , β x k = l x k + β T x k β x k .
being l x k defined as
l x k = x k T P x k f ˜ T x k P f ˜ x k 2 ,
Note that the optimal solution for the functional cost is F 🟉 = V x 0 , being it x 0 the initial condition for the dynamical system (4).
Proof. 
To demonstrate the control law β x k , that is, an optimal function, let us consider the Hamiltonian of the system as
H x k , u k = L x k , β x k + V x k + 1 V x k ,
which has the global minimum as H x k , u k u k = 0 .
To minimize this Hamiltonian function, we can rewrite (23) considering (22) as follows,
min β x k l x k + β T x k β x k + V x k + 1 V x k = 0 .
The solution of the minimization function (25) considering (23) and the variation of the candidate Lyapunov function (12) as presented in [18], it is taken the following form,
h T x k + 2 β T x k j x k + f T x k y k T g T x k P g x k = 0 ,
in addition, if we consider (4b) and (5), then we can simplify (27) as presented below,
β T x k j x k + h T x k j x k + β T x k j T x k j x k = 0 .
It is important to mention that the solution of (27) for β x k takes the following structure,
β x k = I + j x k 1 h x k ,
which confirms control function initially defined in (8) as an optimal control law as it minimizes the functional cost (21).
In order to determine the optimal value for the Lagrangian function (21), let us consider that the interval of analysis 0 , N , being N a natural number with the following result,
k = 0 L x k , β x k = V x N + V x 0 + k = 0 H x k , β x k ,
In the case of the optimal control law β x k , this is optimal if it makes zero the Hamiltonian function H x k , β x k demonstrated in [18]; in addition, we know based on the stability properties of the inverse optimal control that when N the Lyapunov function V x N 0 for any initial condition x 0 , which implies that F 🟉 = V x 0 . □

3.4. General Commentaries

In the application of the studied inverse optimal control it is worthy to mention the following.
To stabilize a nonlinear discrete dynamical system with the form defined in (4) it is used the optimal control law u k = β x k guaranteeing passivity, stability, and optimallity properties.
The application of the inverse optimal control design is subject to the fact that the dynamical system be zero detectable, which can be expressed as presented in Definition 3.
Definition 3.
A system (4) is locally zero-state observable (locally zero-state detectable) if there is a neighborhood Z of x k = 0 R n such that for all x 0 Z
y k | u k = 0 = h ϕ k , x 0 , 0 = 0 k x k = 0 ,
where ϕ k , x 0 , 0 = f k x k is the trajectory of the unforced dynamics x k + 1 = f x k with initial condition x 0 . If Z = R n , the system is zero-state observable (respectively zero-state detectable).

4. Numerical Validation

To demonstrate the effectiveness and robustness of the studied discrete-inverse optimal control for regulation variables in the ball and beam dynamical system, we use the parametric information reported in [2] as presented in Table 1.
With the parameters in Table 1, the coefficients k b b and h are 0.4183 and 0.7143, respectively. The discretization time Δ k is assigned as 1 × 10 3 s. To evaluate the dynamical performance of the proposed approach, the initial condition of the state variables is 15 cm for the ball position, i.e., for x 1 0 , and zero for the rest of state variables.
The parametrization of the control design is based on the structure of the P matrix, which for the ball and beam system can take the following form,
P = r 0 0 j 1 0 r 0 j 2 0 0 r j 3 j 1 j 2 j 3 r
note that P is a positive definite matrix, as it can be rewritten as P = R + J , where R is a positive diagonal and J is skew-symmetry, for r > 0 and j i > 0 .
From the structure of the P matrix, we can observe that there exist four parameters that need to be adjusted to reach the desired dynamical performance, i.e., r and j 1 to j 4 . Here, they are tuned via a trial and error procedure by running multiple simulations in the MATLAB/OCTAVE programming environments.

4.1. Regulation of the State Variables

In this simulation we present the ability of the proposed control to regulate all the state variables using the anti-symmetry nature of the P matrix, where the control gains were assigned as follows, r = 10.37 , j 1 = 16 , j 2 = 20 , and j 3 = 5 .
Figure 2 presents the numerical behavior of each one of the state variables of the ball and beam system after its numerical simulations in MATLAB/OCTAVE software.
From Figure 2 we can observe the following.
All the state variables are regulated when have passed 4000 samples, i.e., ~4 s, which implies that the discrete-inverse optimal control fulfill the control objective when the control input (28) is applied to the discrete equivalent system (3).
The ball position exhibits a smooth dynamical behavior from the initial position (like a second-order system), i.e., x 1 0 = 15 cm, to the origin with a minimum overpass, which implies that the selection of the control gains was appropriate. Nevertheless, this behavior can also be improved (smoothing) if an optimization procedure over gains in P is made as recommended in [13].
The control input u k presented in Figure 2c reaches the zero value when all the state variables are in the origin of coordinates, which is a natural behavior as it is a nonlinear function of all these variables working as a proportional controller that reduces its amplitude when the regulated variables are near to the origin of coordinates.
It is worth mentioning that the proposed controller can regulate all the state variables for the ball and beam system even if this system is strong nonlinear with unique control input and four dynamical equations with a basic implementation in the discrete domain. Figure 3 provides the MATLAB/OCTAVE code implementation that generates all the dynamical behaviors reported in Figure 2.
Note that the main idea of providing the MATLAB/OCTAVE code in Figure 3 is the possibility of using it for possible comparisons in future works by engineering students and researches related to nonlinear control areas.

4.2. Dynamical Performance for Different Control Gains

This subsection explores the effect that has the variations in the control gains in the matrix P regarding the dynamical performance of the ball position. Figure 4 presents some dynamical outputs of the ball position for variations in the control gain j 3 . These variations begin with j 3 = 4 and ends with j 3 = 6 .
The results in Figure 4 demonstrate that the dynamical behavior of the ball position is strongly related to the control gain j 3 , as values of approximately 4 produce higher oscillations in its position that vanish after 6500 samples, i.e., 6.5 s approximately. When this gain increases until 6, the system has low oscillations, and the settling time is ~4.5 s.
It is important to point out that similar behaviors on the ball position can be reached if r or j 2 are also modified because, as aforementioned, the control input (28) is a nonlinear function of the state variables, which implies that these allow governing the closed-loop behavior of the ball and beam system.

4.3. Effect of the Parameter Variations

To evaluate the ability of the proposed discrete-inverse optimal control to deal with parametric uncertainties, we consider that parameters h and k b b have a 20 % of error in their calculations. This implies that in the dynamical model, these parameters are different from the user to design the controller. Figure 5 presents the dynamical performance of the ball position for different combinations of the parameters h and k b b .
From Figure 5, we can observe that the variation in the parameters associated with the physical quantities of the ball and beam system, i.e., h and k b b , has effects on the dynamical behavior of the ball position. Nevertheless, this does not compromise the stability properties ensured by the discrete-inverse optimal control law proposed in this research. In addition, from this simulation, it is possible to note that the setting time remains constant at approximately 4.5 s, which is the same time reported when all the parameters are entirely known. This confirms the robustness of the proposed control approach to possible parametric uncertainties.

5. Conclusions and Future Works

In this paper, the nonlinear discrete-inverse optimal control approach to regulate all the state variables in the discrete version of the ball and beam system was proposed. Numerical results confirm that the proposed controller allows guaranteeing passivity, stability, and optimally properties during closed-loop operation. In addition, numerical simulations show that the dynamical behavior of the system is highly dependent on the control gains contained in the P matrix as one or more of then can produce similar dynamical performances with small oscillations around the origin of coordinates. The control design robustness was tested by varying the parameters k b b and h from 80 % to 120% of their nominal rates, where simulations demonstrate that in all of these cases the setting time remains constant at approximately 4.5 s with smooth oscillations around the control objective, i.e., the origin of coordinates.
As future works it will be possible to develop the following. (i) The experimental validation of the discrete-inverse optimal control and their comparisons with nonlinear control methodologies, (ii) the application of the interconnection and damping assignment control to regulate all the state variables by taking advantage on the open-loop model of the ball and beam via Hamiltonian functions, and (iii) the extension of the proposed discrete-inverse optimal control to the four-order model of the synchronous machine system to regulate the angular deviation when short-circuit scenarios appear in the power system.

Author Contributions

Conceptualization, O.D.M., and W.G.-G.; methodology, O.D.M., and W.G.-G.; formal analysis, O.D.M., and W.G.-G.; investigation, O.D.M., and W.G.-G.; resources, O.D.M., W.G.-G., and C.R.-V.; writing—original draft preparation, O.D.M., and W.G.-G.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the National Scholarship Program Doctorates of the Administrative Department of Science, Technology, and Innovation of Colombia (COLCIENCIAS), by calling contest 727-2015.

Acknowledgments

The authors want to thank to Vicerrectoria de Investigación, Innovación y Extensión from Universidad Tecnológica de Pereira by support given in this investigation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kagami, R.M.; da Costa, G.K.; Uhlmann, T.S.; Mendes, L.A.; Freire, R.Z. A Generic WebLab Control Tuning Experience Using the Ball and Beam Process and Multiobjective Optimization Approach. Information 2020, 11, 132. [Google Scholar] [CrossRef] [Green Version]
  2. Koo, M.S.; Choi, H.L.; Lim, J.T. Adaptive nonlinear control of a ball and beam system using the centrifugal force term. Int. J. Innov. Comput. Inf. Control 2012, 8, 5999–6009. [Google Scholar]
  3. Ye, H.; Gui, W.; Yang, C. Novel Stabilization Designs for the Ball-and-Beam System. IFAC Proc. Vol. 2011, 44, 8468–8472. [Google Scholar] [CrossRef] [Green Version]
  4. Rahmat, M.F.; Wahid, H.; Wahab, N.A. Application of intelligent controller in a ball and beam control system. Int. J. Smart Sens. Intell. Syst. 2010, 3, 45–60. [Google Scholar] [CrossRef] [Green Version]
  5. Meenakshipriya, B.; Kalpana, K. Modelling and Control of Ball and Beam System using Coefficient Diagram Method (CDM) based PID controller. IFAC Proc. Vol. 2014, 47, 620–626. [Google Scholar] [CrossRef]
  6. Ding, M.; Liu, B.; Wang, L. Position control for ball and beam system based on active disturbance rejection control. Syst. Sci. Control. Eng. 2019, 7, 97–108. [Google Scholar] [CrossRef] [Green Version]
  7. Qi, X.; Li, J.; Xia, Y.; Wan, H. On stability for sampled-data nonlinear ADRC-based control system with application to the ball-beam problem. J. Franklin Inst. 2018, 355, 8537–8553. [Google Scholar] [CrossRef]
  8. Chen, C.C.; Chien, T.L.; Wei, C.L. Application of Feedback Linearization to Tracking and Almost Disturbance Decoupling Control of the AMIRA Ball and Beam System. J. Optim. Theory Appl. 2004, 121, 279–300. [Google Scholar] [CrossRef] [Green Version]
  9. Aguilar-Ibañez, C.; Suarez-Castanon, M.S.; de Jesús Rubio, J. Stabilization of the Ball on the Beam System by Means of the Inverse Lyapunov Approach. Math. Probl. Eng. 2012, 2012, 1–13. [Google Scholar] [CrossRef]
  10. Li, E.; Liang, Z.Z.; Hou, Z.G.; Tan, M. Energy-based balance control approach to the ball and beam system. Int. J. Control 2009, 82, 981–992. [Google Scholar] [CrossRef]
  11. Kelly, R.; Sandoval, J.; Santibáñez, V. A Novel Estimate of The Domain of Attraction of an IDA-PBC of a Ball and Beam System. IFAC Proc. Vol. 2011, 44, 8463–8467. [Google Scholar] [CrossRef] [Green Version]
  12. Chang, Y.H.; Chang, C.W.; Tao, C.W.; Lin, H.W.; Taur, J.S. Fuzzy sliding-mode control for ball and beam system with fuzzy ant colony optimization. Expert Syst. Appl. 2012, 39, 3624–3633. [Google Scholar] [CrossRef]
  13. Castillo, O.; Lizárraga, E.; Soria, J.; Melin, P.; Valdez, F. New approach using ant colony optimization with ant set partition for fuzzy control design applied to the ball and beam system. Inf. Sci. 2015, 294, 203–215. [Google Scholar] [CrossRef]
  14. Zavala, S.J.; Yu, W.; Li, X. Synchronization of two ball and beam systems with neural compensation. IFAC Proc. Vol. 2008, 41, 12781–12786. [Google Scholar] [CrossRef]
  15. Ornelas, F.; Loukianov, A.G.; Sanchez, E.N. Discrete-time Robust Inverse Optimal Control for a Class of Nonlinear Systems. IFAC Proc. Vol. 2011, 44, 8595–8600. [Google Scholar] [CrossRef] [Green Version]
  16. Atkinson, C.; Osseiran, A. Discrete-space time-fractional processes. Fract. Calc. Appl. Anal. 2011, 14. [Google Scholar] [CrossRef]
  17. Carrasco-Gutierrez, C.E.; Sosa, W. A discrete dynamical system and its applications. Pesqui. Oper. 2019, 39, 457–469. [Google Scholar] [CrossRef]
  18. Sanchez, E.N.; Ornelas-Tellez, F. Discrete-Time Inverse Optimal Control for Nonlinear Systems; CRC Press Taylor and Francis Group: Boca Raton, FL, USA, 2017. [Google Scholar]
  19. Galor, O. Discrete Dynamical Systems; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar] [CrossRef]
  20. Gil-González, W.; Serra, F.M.; Montoya, O.D.; Ramírez, C.A.; Orozco-Henao, C. Direct Power Compensation in AC Distribution Networks with SCES Systems via PI-PBC Approach. Symmetry 2020, 12, 666. [Google Scholar] [CrossRef] [Green Version]
  21. la Sen, M.D. On the Passivity and Positivity Properties in Dynamic Systems: Their Achievement under Control Laws and Their Maintenance under Parameterizations Switching. J. Math. 2018, 2018, 1–16. [Google Scholar] [CrossRef] [Green Version]
  22. Navarro-López, E.; Fossas-Colet, E. Dissipativity, passivity, and feedback passivity in the nonlinear discrete-time setting. IFAC Proc. Vol. 2002, 35, 257–262. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic representation of the ball and beam dynamical system.
Figure 1. Schematic representation of the ball and beam dynamical system.
Symmetry 12 01359 g001
Figure 2. Time-domain behavior of the state variables and control input in the ball and beam system: (a) position and and speed of the ball, (b) angular position and angular speed of the beam, and (c) control input.
Figure 2. Time-domain behavior of the state variables and control input in the ball and beam system: (a) position and and speed of the ball, (b) angular position and angular speed of the beam, and (c) control input.
Symmetry 12 01359 g002
Figure 3. MATLAB/OCTAVE implementation of the inverse optimal controller for the ball and beam system.
Figure 3. MATLAB/OCTAVE implementation of the inverse optimal controller for the ball and beam system.
Symmetry 12 01359 g003
Figure 4. Dynamical behavior of the ball position for different values of the gain j 3 .
Figure 4. Dynamical behavior of the ball position for different values of the gain j 3 .
Symmetry 12 01359 g004
Figure 5. Dynamical behavior of the ball position when parametric variations are experimented in h and k b b parameters.
Figure 5. Dynamical behavior of the ball position when parametric variations are experimented in h and k b b parameters.
Symmetry 12 01359 g005
Table 1. Parameter information of the ball and beam system.
Table 1. Parameter information of the ball and beam system.
ParameterValueUnity
L beam 42.55cm
r arm 2.54cm
R1.27cm
m64mg
g9.81m/s2
J b 4.1290 × 10 6 kgm2
k 1 1.76rad/sv
τ 28.5ms

Share and Cite

MDPI and ACS Style

Danilo Montoya, O.; Gil-González, W.; Ramírez-Vanegas, C. Discrete-Inverse Optimal Control Applied to the Ball and Beam Dynamical System: A Passivity-Based Control Approach. Symmetry 2020, 12, 1359. https://doi.org/10.3390/sym12081359

AMA Style

Danilo Montoya O, Gil-González W, Ramírez-Vanegas C. Discrete-Inverse Optimal Control Applied to the Ball and Beam Dynamical System: A Passivity-Based Control Approach. Symmetry. 2020; 12(8):1359. https://doi.org/10.3390/sym12081359

Chicago/Turabian Style

Danilo Montoya, Oscar, Walter Gil-González, and Carlos Ramírez-Vanegas. 2020. "Discrete-Inverse Optimal Control Applied to the Ball and Beam Dynamical System: A Passivity-Based Control Approach" Symmetry 12, no. 8: 1359. https://doi.org/10.3390/sym12081359

APA Style

Danilo Montoya, O., Gil-González, W., & Ramírez-Vanegas, C. (2020). Discrete-Inverse Optimal Control Applied to the Ball and Beam Dynamical System: A Passivity-Based Control Approach. Symmetry, 12(8), 1359. https://doi.org/10.3390/sym12081359

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop