Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 3820935, 11 pages
https://doi.org/10.1155/2018/3820935

Research Article

On the Output Regulation Problem: The Generalized Second-Order Underactuated Linear System Case

Carlos Aguilar-Ibanez,¹ Jesus A. Meda-Campana,² Miguel S. Suarez-Castanon,³ Jose de Jesus Rubio,⁴ and Nareli Cruz-Cortes¹

¹Instituto Politecnico Nacional, Centro de Investigacion en Computacion, Av. Juan de Dios Batiz s/n, U.P.A.L.M, Ciudad de Mexico 07738, Mexico
²Instituto Politecnico Nacional, SEPI-ESIME Zacatenco, Av. Instituto Politecnico Nacional s/n, U.P.A.L.M, Ciudad de Mexico 07738, Mexico
³Instituto Politecnico Nacional, Escuela Superior de Computo, Av. Juan de Dios Batiz s/n, U.P.A.L.M., Ciudad de Mexico 07738, Mexico
⁴Instituto Politecnico Nacional, ESIME Azcapotzalco, Av. de las Granjas 682, Ciudad de Mexico 02250, Mexico

Correspondence should be addressed to Miguel S. Suarez-Castanon; sasuarez@prodigy.net.mx

Received 7 March 2018; Revised 10 September 2018; Accepted 14 October 2018; Published 30 October 2018

Academic Editor: Konstantinos Karamanos

Copyright © 2018 Carlos Aguilar-Ibanez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this work, we aimed to solve the control regulation problem for a generalized second-order underactuated linear system in order to induce a periodic or chaotic behavior or to cancel the external perturbations, generated by an exogenous system, in the nonactuated coordinate. Further, we showed that, in some cases, it is possible to bring to zero the regulation output errors of the underactuated linear plant, depending on the structure of the plant itself and the exogenous system. In the first stage, the solution was developed for the ideal scenario, in which the whole states of the plant and of the exogenous system were available. Secondly, we showed that in some cases it was possible to solve the regulation output problem when only the observable plant output was measurable. That is, the whole plant state and the exogenous signal could be recovered, if some assumptions were fulfilled. The Lyapunov method was used to perform the stability analysis. The proposed solution was assessed through numerical simulations.

1. Introduction

In actual control applications, it is common to deal with external perturbations that affect the input of the system or the plant to be controlled. Usually, these perturbations are generated by a different external source, a so-called exogenous system. Due to the importance of this issue, a theory to deal with it has been developed. In this scenario, when the plant has to follow a trajectory, the problem is referred to as the output regulation control problem. Particularly, the exogenous system must exhibit an oscillatory behavior. In more rigorous terms, the regulation problem consists of finding a state or error feedback controller, such that asymptotic stability of the closed-loop system can be achieved, if the external exogenous signal is absent; however, in the presence of the perturbation signal generated by the exogenous system, error tracking convergence to zero can be accomplished (for a detailed treatment of this topic, please refer to [1]).

We outline here the most remarkable works related to this topic. The pioneering work [2] introduced a solution to the linear regulation problem. The same paper introduced the important notion of the linear robust regulator (a controller capable of preserving the regulation property, even when the system parameters are time variant). The authors also showed that, for this case, the linear robust regulation may be achieved, if the controller contains an internal model of the exogenous system. In [3, 4], the regulation problem was studied for the class of nonlinear systems with constant exogenous signal. Isidori and Byrnes, in [5], studied the output regulation problem of nonlinear systems with time-varying exogenous signals. Additionally, in this work, the authors probed whether a nonlinear regulator can be obtained on the basis of a set of partial differential equations, known as the FIB equations, named after Francis, Isidori, and Byrnes. A similar work was presented by Huang and Rugh in [6], where multivariable nonlinear plants with measured and unmeasured disturbance signals where considered. The following are some of the most important works related to robust output regulation: [7–12]. Finally, several control strategies for the output regulation problem can be found in the literature; however, many of them work with mathematical model approximations or switching the local models using techniques like fuzzy control [1, 13, 14]. A full review of this matter is beyond the scope of this work; however, we suggest to the interested reader the book [15].

The control regulation problem for a generalized second-order underactuated linear system perturbed by an external signal is studied in this work from the perspective of the regulation theory. The output regulation problem has been solved in the case when the system is fully actuated. Here, we are interesting in the instance where the system is not fully actuated (i.e., where there are fewer actuators than degrees of freedom). Salient to this work we found that not-fully actuated systems can be partially regulated, depending on the nature of both the plant and the exogenous system, and the interconnection matrix. Even more, we provide the necessary conditions to accomplish the output regulation for this kind of system. In other words, we were able to assure that the regulation errors and , where and are the generalized measurable position and the generalized velocity, respectively, is the exogenous system signal, and is a suitable constant matrix, susceptible to being designed with some degree of freedom. We underscore that the perturbation signal has a periodic, quasi-periodic, or chaotic behavior. We must note that, if the regulation errors go to zero, we can induce, evidently, either a chaotic or a periodic behavior, and in some cases we can cancel the effect of in any of the coordinates, depending on the nature of the matrix . The corresponding stability analysis was accomplished using the Lyapunov theory. It is worth mentioning that the obtained results can be extended to underactuated nonlinear systems, if they are locally controllable and the exogenous signal is sufficiently small. We numerically verified the effectiveness of the proposed solution using, as testbeds, the Vibrating Mechanical System (VMS) and the Inverted Pendulum Cart System (IPCS), the latter perturbed by small exogenous signals, whose time derivatives were also small. As outcomes we were able to induce a periodic and a chaotic behavior and cancel the perturbation signal in the nonactuated coordinate, even when the perturbation acted solely in the nonactuated acceleration.

At this point, it is important to remark that, in [16], a nonlinear regulator designed on the basis of fuzzy models was applied on the pendubot, which is an underactuated system. However, in that work the regulator was designed supposing that the plant and exosystem’s whole state was available (full information case). Besides, perturbations were not considered, and the dimension of the input was equal to the dimension of the output. Therefore, the results obtained in the present work may help to increase the applicability range of the regulation theory on underactuated systems.

The remainder of this work is organized as follows. In Section 2, we establish the output regulation for a generalized second-order underactuated linear system. In Section 3, we solve the established control problem in its ideal case and introduce the main proposition of this work; we also give an extension related to the main proposition, where we deal with the instance in which the state is partially known. The numerical simulations for assessing the effectiveness of the proposed solution are presented in Section 4. Concluding remarks are given in Section 5.

2. Formulation Control Problem

Let us consider the following generalized second-order underactuated linear system, modeled as

where and are the states; is the control input with ; , , , and are fixed constant matrices; and is the external signal generated by a chaotic or quasi periodic exogenous system, defined by

where is a smooth function. The control objective consists of designing a controller , formed by

where constants , , and , such that the regulation errors

converge to zero for any initial condition, whereby the constant matrix is conveniently fixed.

The aforementioned problem is closely related to the linear regulator theory (LRT), widely studied in [5, 17–19]. The main differences between these two problems are as follows: (i) the LRT was formulated assuming that the signals and have the same dimensions. In our case, we studied the instance when the dimension of is strictly less than the dimension of . That is, ; (ii) the LRT was formulated under the assumption that the exogenous system is linear (2). Conversely, in our approach the exosystem can be any attractor belonging in a compact set.

Motivation: from our point of view, this control problem has at least two relevant applications in engineering.

First application: for inducing a chaotic or quasi periodic behavior.

It is an easy way to induce oscillations, either quasi-periodic or chaotic, in some class of underactuated flat systems, or for locally controllable nonlinear underactuated systems. As we seen from equation (4), if and converge to zero, then, in the steady state, the position and velocity vectors are forced to exhibit a quasi-periodic or chaotic behavior, according to the characteristics of the signal .

Second application: for attenuating disturbances.

It consists of reducing the effect of the external signal through the control signal . In other words, we want to attenuate, as much as possible, the effect of in the system (1) through the action of the controller (3), avoiding the necessity of using either a high-gain controller or a sliding-mode controller. In fact, the control action of forces the convergence to zero of the regulation error:

keeping the system stable. Notice that if , then, in the steady state, the following inequality is fulfilled:

Hence, selecting the smallest value of all possible values of the norm of , we can assure that the effect of will be smaller on the states and . For instance, if , with , then it is easy to see that the following inequalities hold:

with and . That is, the smaller is, the smaller and are. In the developments below we give some conditions to find the matrix that assures that .

Please note that, in model (1), it is not possible to completely eliminate the effect produced by due to two facts: first, the matrix is not invertible, since ; second, is an unmatched perturbation.

We end this section by introducing some useful assumptions:

(A1) The par is stabilizable.

(A2) The constant matrices of system (1) and the structure of are known.

(A3) For the exogenous system (2), we have that , , and are well defined and bounded for all future time.

3. Problem Development

In this section, we first solve the above-established control problem in its ideal case, in which the states are always available for measurements. Afterwards, we show that this problem can also be solved using our scheme, when only the state is measurable and can always be estimated. Even more, if the interconnection matrix is not singular, it is possible to recover the states , provided that only is available for measurement. Notice that is the positions vector of the linear underactuated system.

3.1. Solving the Output Regulation Problem for the Underactuated Linear System

First of all, the closed-loop system, defined by (1) and (3), can be rewritten as

Therefore, the time derivative of the regulation error , defined in (4), satisfies the following equation:

with

Substituting the values and , obtained from definition (4), into (11), we obtain:

where

To assure that asymptotically and exponentially converges to zero, we first select and , such that becomes a Hurwitz matrix. This fact can always be fulfilled according to (A1). Next, the variables and must be selected, such that

Indeed, if (13) is fulfilled and is Hurwitz, then the error goes to zero. Next, to derive the conditions to obtain and needed to satisfy the above equation, we multiply both sides of (13) by the following nonsingular matrix:

where is the left annihilator of (that is, ). The latter leads to the following implications:

and

which leads to the following equation:

where . We pointed out that must be chosen, such that (15) is fulfilled for any . Now, if we found such , then can be straightforwardly obtained through (17), allowing us to claim that the regulation error, , converges to zero. Finally, according to relation (4), the following holds:

In other words, in the steady state the following inequalities also holds:

We summarize the previous discussion in the form of the following proposition.

Proposition 1. Consider the interconnected systems (1) and (2), in closed-loop with (3), and under the assumptions (A1), (A2), and (A3). Then, the regulation error can be accomplished if there exists a matrix that satisfies the matching equation (15).

Comment 1. Concerning (15), we pointed out that, if , with being a constant matrix and , then equation (15) leads to

It is worth mentioning that, in general, given

for systems (1) and (2), there does not always exist a matrix that satisfies the matching equation (15); besides, it is not unique. Then, there are two possibilities, which we analyze below.

(A)   can be selected to solve the disturbances attenuation problem, adding the following condition:

where is the admissible set, defined by

(B)   is accomplished provided that the nonactuated coordinate of converges to zero. That is, it is possible to completely eliminate the effects of the exogenous signal over the nonactuated coordinate of the state .

Comment 2. If there exists a matrix different from zero, then, according to the inequality (19), it can be possible to induce in the states and either a chaotic or a periodic behavior, depending on the nature of the exogenous system (2).

4. Solving the Control Problem When the State Is Available

In the previous discussion, we solved the output regulation problem for the underactuated linear system (1), when all the states were available for measurements. Now, we consider the case when only the state is available, which entails two cases: (i) the matrix is singular, is available for measurements, and the state is recovered by means of a simple Luenberger observer; (ii) the matrix is not singular, and the states and are not available. To solve the latter case, a Luenberger observer and a high-gain sliding-mode observer that work simultaneously are designed to estimate asymptotically the states and .

First Case. is singular, and and are available. To solve this case, we use the following controller:

with evolving according to the following Luenberger observer:

where is the observation error of the measurable position, and , and

Defining the observation error as

it is easy to see from (1) and (25) that their dynamics fulfill the following equation:

Observe that is Hurwitz, as long as and fulfill the following inequalities:

Therefore, according to (28a), we can claim that exponentially and asymptotically converges to zero. Hence, if and are selected according to equations (15) and (17), and following the procedure used to obtain (15), it is easy to see that the closed-loop system, defined by (1) and (24), is equivalent to the following equation:

Now, since the error exponentially converges to zero and is Hurwitz, we concluded that also converges exponentially to zero. That is, the output tracking error goes to zero when the state is estimated through the Luenberger observer (25).

Second Case. is nonsingular and is measurable. Here we propose the following controller:

where and will be proposed, such that and , for , with finite. To implement this controller, we introduce the following assumption: (A4) is invertible. Hence, the Luenberger observer is proposed as

with and fulfilling the conditions in (29). Once again, from (1) and (32), we obtain the following dynamic observation error:

Evidently, if is Hurwitz and is bounded by (A3), then the solution of (33) satisfies the inequality , for all . In other words, if (29) is fulfilled, we can assure that the observation error is ultimately bounded and stable. Based on the Luenberger observer (32), we give formulas to explicitly recover and . To this end, we have that (33) is equivalent to the following:

and

From the above, it is easy to see that (35) can be rewritten in terms of as

Notice that, from the last equation, once again are bounded. Now, if assumption (A4) holds, signal can be algebraically obtained as

Similarly, using (34), the corresponding algebraic formula for is given by

Since the time derivatives and are not available for measurement in (37) and (38), we need to propose the following reconstructors to recover them:

and

where the signals and are proposed, such that

In the next section, we use a high-order sliding-mode exact differentiator to obtain the signals and (42).

Comment 3. If and fulfill the condition (41), then we also have that and , for , with finite, assuring the exact recovering of the nonavailable signals and . On the other hand, if assumption (A4) is not fulfilled, we can only exponentially recover the state , through the Luenberger observer (25).

The sliding-mode high-order differentiator: based on Levant’s work [20], we propose the following exact differentiator:

where

Under the assumption that there exists a constant , with , for , we have that, for some , the following equalities hold:

as long as the gains are selected as

Comment 4. According to the inequalities in (29), (36), and the fact that and are both bounded, we have that , for , always holds. In fact, the actual value of cannot be accurately computed, but we can always select some sufficiently large to overcome this inconvenience.

We end this section summarizing the above discussion in the form of the following lemma.

Proposition 2. Consider system (36), with measurable output , under assumption (A3), with and fulfilling (29). Then, for sufficiently large , and selected according to (45), we have that the sliding-mode high-order differentiator (42) assures the conditions in (44). Consequently, the estimation given in (39) converges to in finite time. Complementarily, if assumption (A4) is fulfilled, then the estimation of also converges to in finite time.

Comment 5. The closed-loop system (9) can be seen as a second-order linear system, perturbed by a bounded signal , with exponentially stable. Evidently, as long as the estimations of and converge to their corresponding actual values in finite time, all the solutions of the closed-loop system, in combination with the Luenberger observer (32) and the exact differentiator (42), are bounded and well defined, for all future, fulfilling the separation principle. Consequently, the output regulation problem is accomplished.

5. Numerical Examples

To illustrate the performance of the proposed solution to the output regulation problem, we designed three numerical experiments. The first two use the underactuated Vibrating Mechanical System (VMS), and the third one uses the traditional Cart-Pendulum system (CPS). For consistency, we consider that both systems are affected by two external vibration signals, and , generated through

where is the vector of the exogenous signals.

5.1. Description of VMS for the First Two Experiments

The VMS, depicted in Figure 1, is made up of two coupled subsystems. The secondary subsystem has an active undamped dynamic vibration absorber and is coupled to the primary perturbed system. The dynamics of this system are described by the following set of differential equations:

where the variables and are, respectively, the position and the velocity of mass , in the same way that and are, respectively, the position and velocity the of mass . The springs’ linear stiffness values are denoted by and , and is the viscous damping; is a constant interconnection matrix; is the vector of the exogenous perturbation signals, generated by some target system and proposed in each experiment; and represents the system input or control action. Summarizing, , , and are, respectively, the positions vector, the velocities vector, and the exogenous signals vector. To simplify the following developments, we fixed the system parameters as . Hence, according to (1), we have that the corresponding characteristic matrices for the VMS are given as

Notice that is the actuated coordinate and is the nonactuated coordinate.

Figure 1: The underactuated VMS.

First Experiment. For this experiment, we assume that an external vibration signal, , is present in the acceleration of the nonactuated coordinate, (evidently, for all ). The goal of this experiment consists of canceling the vibration in the nonactuated coordinate, assuming that the vector position, , and are available. In other words, we consider that and , where is generated by the following system:

Evidently, for this case, the interconnection matrix is given by

Therefore, the matrix must satisfy, according to (20), the following equation:

After substituting the corresponding values of , given in (48), (49), and (50), into equation (51), we obtain

where the parameters and must be fixed, such that they assure that and . That is, we want to cancel the effect of the exogenous signal in the nonactuated coordinates and . This can be achieved by selecting (see equation (4)). To obtain the nonavailable variable , we use the observer (25), where and assure the condition (29). The initial condition for the system is given by and ; the exogenous system was initialized in ; and the corresponding control gains were set as and . In Figure 2, we show the outcome of this experiment, for , , and . In this figure, we can see the closed-loop behavior, the evolution of both the actuated and nonactuated coordinates, and the corresponding errors, and . We can also see that converges to after seconds, while the errors and converge to zero after seconds. In addition, we can see that the nonactuated coordinates remain oscillating in steady state after seconds. We can see that the control reaches a steady state at seconds.

Figure 2: Cancelation of the exogenous signal in the perturbed nonactuated coordinate of the VMS.

Second Experiment. This experiment consists of inducing a chaotic behavior in both the actuated and the nonactuated coordinates of the VMS, through the exogenous signals and , where is generated through the Genesio-Tesi attractor (GTA). This control task will be carried out assuming that the state is available. The GTA [21, 22] is described by the following set of equations:

where , , and are strictly positive constants, such that , and the interconnection matrix is fixed as

Evidently, matrix is not invertible and, therefore, does not fulfill assumption (A4). Consequently, the state cannot be recovered; hence, it has to be measured. It is easy to see that the GTA has the following Jacobian matrix:

For simplicity, we fixed the parameters

To find the needed matrix , defined as

we substitute the expressions (54), (48), and (57) into the single equation (15), leading to the following nonlinear function:

where

To make (58) hold for any , we must assure that , . The latter implies the following:

where is a free parameter. To carry out the second experiment, the gains , , and were set as in the first experiment. The parameters of the interconnection matrix were fixed as and , for . The initial condition for the GTA was initialized as , while the VMS was initialized as in the first experiment. The outcome of this experiment is shown in Figure 3, which shows the tracking errors of the actuated and nonactuated coordinates. In the same figure, we show the phase portraits of both coordinates. From this figure, we can see that the error goes very close to zero after seconds for both the actuated and the nonactuated coordinates. Additionally, we can see, from the phase portraits, that the VMS system effectively exhibits a chaotic behavior, even when the simulation time took seconds. It is worth mentioning that the experiments were carried out using a computer simulation environment that uses a sixth-order Runge-Kutta algorithm, with an integration step of . This configuration helped us see the chaotic behavior of the GTA. In fact, if we use a greater integration step, we cannot appreciate the chaotic behavior of this system.

Figure 3: Inducing a chaotic behavior in the actuated and nonactuated coordinates of the VMS, through the Genesio-Tesi system.

5.2. Description of PCS for the Third Experiment

To carry out this experiment, we consider the normalized partially feedback model of the PCS ([23]), described by the following set of differential equations:

where is the actual angular position of the pendulum with respect to the vertical, and is the normalized horizontal displacement of the cart, while their corresponding angular and horizontal velocities are and , respectively. The virtual controller is defined as , and has the same definition as in (46), provided that is sufficiently small to make the system to oscillate around to the unstable equilibrium point .

Third Experiment. This experiment consists of inducing a limit cycle using solely measurements of the pendulum and the cart positions. The limit cycle is produced by the linear exogenous system, defined as

where can be any positive constant, and the interconnection matrix is fixed as . In this case, we assume that is available for measurements and that will be recovered using the Luenberger observer defined in (32), while will be recovered using the high-order exact differentiator (42) and the reconstructor (40). To make our methodology adequate, we use the linearized model of the system (61), given by

Mimicking the steps as in the first experiment and, according to the matching equation (20), we have that

where

with defined in (62) and . Substituting , , , and into equation (64), we obtain the following:

where and can be fixed as desired. To accomplish the goal of this experiment, which is to assure that and , we use the following set up. The gains , , and were fixed as in the first experiment, while , , and . The initial condition of the exogenous system was set as , and the PCS was initialized at and . Because is nonsingular, the PCS fulfills the conditions of the assumption (A4). Therefore, we propose the parameter , for the sliding-mode high-order differentiator (42). The obtained results can be seen in Figure 4, where we present the tracking errors of the actuated and the nonactuated coordinates. In the same figure, we show the phase portraits of both coordinates. From this figure, we can see that the tracking error goes very close to zero after seconds, for both the actuated and the nonactuated coordinates. Additionally, we can see, from the phase portraits, that the PCS system effectively exhibits a limit cycle in both the pendulum angular coordinate and the cart horizontal coordinate, even when the controller uses the estimated velocity vector and the reconstructed exogenous perturbation vector .

Figure 4: Tracking errors for the angular and horizontal positions for the PCS and corresponding phase portraits of these positions.

6. Concluding Remarks

In this work, we studied the control regulation problem for a generalized second-order underactuated linear system perturbed by an external signal, generated by a chaotic or quasi-periodic exogenous system. We underscore that the well-known problem can be solved provided that the number of controllers coincides with the number of degrees of freedom. In this study, we analyzed the case in which we have fewer actuators than degrees of freedom. Under this strong assumption, we conclude that, in some cases, it is possible to assure that the regulation errors , introduced in (4), can be brought to zero relying on the system structure, which in turn depends on the nature of the exogenous system and the interconnection matrix , as we established in Proposition 1. To obtain the solution to the aforementioned problem, we worked with the ideal case, in which all the state variables and the exogenous signal were available for measurements. Then, we relaxed this condition assuming that solely the position vector, , and the exogenous signal, , were available, while the velocity vector, , can be recovered by a Luenberger observer. Besides, if the interconnection matrix is not singular, we can recover both and , only based on measurements of . It is worth mentioning that this control scheme, which is based on the linear regulation theory, can be applied to induce either periodic or chaotic behaviors in the underactuated linear system. Also, the scheme allows us to cancel in some cases the effect of the exogenous signal in the nonactuated coordinate, through appropriate control actions. This was illustrated in the first numerical simulation, where we used the VMS, which was subject to a sinusoidal exogenous signal in the nonactuated coordinate, and where the control task consisted of canceling this effect. The corresponding stability analysis was carried out applying the Lyapunov method. Finally, we performed three convincing numerical simulations, where we showed, first, that our control scheme effectively canceled the undesirable exogenous perturbations in the nonactuated coordinate and, second, that a periodic or chaotic behavior was induced.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the Secretaria de Investigacion y Posgrado del IPN, under Research Grants 20181776 and 20181101.

References

J. A. Meda-Campaña, J. C. Gómez-Mancilla, and B. Castillo-Toledo, “Exact output regulation for nonlinear systems described by Takagi-Sugeno fuzzy models,” IEEE Transactions on Fuzzy Systems, vol. 20, no. 2, pp. 235–247, 2012.

B. A. Francis and W. M. Wonham, “The internal model principle of control theory,” Automatica, vol. 12, no. 5, pp. 457–465, 1976.

E. J. Davison, “The Robust Control of a Servomechanism Problem for Linear Time-Invariant Multivariable Systems,” IEEE Transactions on Automatic Control, vol. 21, no. 1, pp. 25–34, 1976.

B. A. Francis and W. M. Wonham, “The internal model principle for linear multivariable regulators,” Applied Mathematics and Optimization. An International Journal with Applications to Stochastics, vol. 2, no. 2, pp. 170–194, 1975.

A. Isidori and C. I. Byrnes, “Output regulation of nonlinear systems,” IEEE Transactions on Automatic Control, vol. 35, no. 2, pp. 131–140, 1990.

J. Huang and W. J. Rugh, “On a nonlinear multivariable servomechanism problem,” Automatica, vol. 26, no. 6, pp. 963–972, 1990.

C. I. Byrnes, F. Delli Priscoli, A. Isidori, and W. Kang, “Structurally stable output regulation of nonlinear systems,” Automatica, vol. 33, no. 3, pp. 369–385, 1997.

F. Delli Priscoli, “Robust tracking for polynomial plants,” in Proceedings of the European Control Conference, pp. 369–373, 1993.

J. Huang, “Asymptotic tracking and disturbance rejection in uncertain nonlinear systems,” Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 40, no. 6, pp. 1118–1122, 1995.

J. Huang and C.-F. Lin, “On a robust nonlinear servomechanism problem,” Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 39, no. 7, pp. 1510–1513, 1994.

J. Huang and C.-F. Lin, “Internal model principle and robust control of nonlinear systems,” in Proceedings of the 32nd IEEE Conference on Decision and Control. Part 2 (of 4), pp. 1501–1506, December 1993.

H.-K. Khalil, “Robust servomechanism output feedback controllers for feedback linearizable systems,” Automatica, vol. 30, no. 10, pp. 1587–1599, 1994.

B. Castillo-Toledo and J. A. Meda-Campaña, “The fuzzy discrete-time robust regulation problem: a LMI approach,” in Proceedings of the 41st IEEE Conference on Decision and Control, vol. 2, pp. 2159–2164, IEEE, Las Vegas, Nev, USA, December 2002.

N. R. Cazarez-Castro, L. T. Aguilar, and O. Castillo, “Fuzzy logic control with genetic membership function parameters optimization for the output regulation of a servomechanism with nonlinear backlash,” Expert Systems with Applications, vol. 37, no. 6, pp. 4368–4378, 2010.

J. C. Doyle, B. A. Francis, and A. R. Tannenbaum, Feedback Control Theory, Courier Corporation, 2013.

J. A. Meda-Campaña, B. Castillo-Toledo, and V. Zúñiga, “On the nonlinear fuzzy regulation for under-actuated systems,” in Proceedings of the 2006 IEEE International Conference on Fuzzy Systems, pp. 2195–2202, Canada, July 2006.

B. A. Francis, “The linear multivariable regulator problem,” SIAM Journal on Control and Optimization, vol. 15, no. 3, pp. 486–505, 1977.

L. E. Ramos, B. Castillo-Toledo, and J. Alvarez, “Nonlinear regulation of an underactuated system,” in Proceedings of the 1997 IEEE International Conference on Robotics and Automation, ICRA. Part 4 (of 4), pp. 3288–3293, April 1997.

B. Castillo-Toledo, S. Celikovsky, and S. Di Gennaro, “Generalized immersion and nonlinear robust output regulation problem,” Kybernetika, vol. 40, no. 2, pp. 207–220, 2004.

A. Levant, “Higher-order sliding modes, differentiation and output-feedback control,” International Journal of Control, vol. 76, no. 9-10, pp. 924–941, 2003.

A. Luviano-Juarez, J. Cortes-Romero, and H. Sira-Ramirez, “Synchronization of chaotic oscillators by means of generalized proportional integral observers,” International Journal of Bifurcation and Chaos, vol. 20, no. 5, pp. 1509–1517, 2010.

A. G. Soriano-Sanchez, C. Posadas-Castillo, M. A. Platas-Garza et al., “Coupling strength computation for chaotic synchronization of complex networks with multi-scroll attractors,” Applied Mathematics and Computation, vol. 275, pp. 305–316, 2016.

C. Aguilar-Ibañez, J. C. Martínez-García, A. Soria-López et al., “On the Stabilization of the Inverted-Cart Pendulum Using the Saturation Function Approach,” Mathematical Problems in Engineering, vol. 2011, Article ID 856015, 14 pages, 2011.