Nonlinear H∞ synthesis is developed to solve the tracking control problem restricted to
a two degrees-of-freedom (DOF) underactuated mechanical manipulator where position measurements are the only available information for feedback. A local H∞ controller is derived
by means of a certain perturbation of the differential Riccati equations, appearing in solving
the H∞ control problem for the linearized system. Stabilizability and detectability properties
of the control system are thus ensured by the existence of the proper solutions of the unperturbed differential Riccati equations, and hence the proposed synthesis procedure obviates an extra verification work of these properties. Due to the nature of the approach, the resulting controller additionally yields the desired robustness properties against unknown but bounded
external disturbances. The desired trajectory is centered at the upright position where the
manipulator becomes a nonminimum-phase system. Simulation results made for a double
pendulum show the effectiveness of the proposed controller.
1. Introduction
The focus of this paper is to solve the tracking control problem for a 2-DOF underactuated mechanism via nonlinear H∞-control for time-varying systems [1] where joint position measurements are the only available information for feedback. Further research applications in the control of underactuated systems have gone in many directions, for example, fully actuated robots where it is required that motion continues in spite of a failure of any of its actuators. Other typical examples are the systems where its desired operation mode is oscillatory such as biped walking robots where a periodic trajectory is required to produce a coordinated motion (see, e.g., [2]); hopping robots where thrust, decompression, flight, and compression phases are also governed by a periodic motion (see, e.g., [3]); tracking control in drive systems with backlash where usually the position sensor is placed on the side of the motor instead of the side of the load (see, e.g., [4] and [5, page 456]); juggling systems [6]; among others.
Objective
In the present paper, we address the output tracking control problem in nonminimum-phase underactuated mechanical system. Representative works in this topic include orbital stabilization of underactuated systems by means of reference models as generator of limit cycles (see, e.g., [7–10]). In particular, this paper is devoted to the solution of a periodic balancing problem for a two-link underactuated mechanical manipulator introduced in [7], whose first link is not actuated whereas the second joint is actuated.
Contribution
For nonlinear mechanical systems, tracking control problem is known to be more difficult than stabilization mainly for underactuated systems whose initial conditions are close to an unstable equilibrium point. The central problem in nonminimum-phase underactuated systems, solved here, is the specification and design of output feedback inner-tracking controllers to drive the output (joint position) to a nontrivial reference trajectory in spite of external disturbances.
The prior work on the tracking control of nonminiminum-phase systems includes, among others, the results of Consolini and Tosques [11] and Berkemeier and Fearing [7] who developed an exact tracking control via state-feedback. Wang [12] partially addressed the above problem by considering the regulation problem in linear systems. A unified treatment of the control of such systems via output feedback can be found in [13]. In the present paper, the nonlinear H∞ control approach is extended for the time-varying nonlinear nonminimum-phase systems applied to tracking control problems for underactuated mechanical systems.
Methodology
The method we use for defining a desired trajectory for underactuated system is based on the work of Berkemeier and Fearing [7]. The method was successfully applied to derive a set of exact trajectories for the nonlinear equation which involve inverted periodic motion. This method was selected because the desired trajectories are at least twice-differentiable satisfying the smoothness assumption imposed on the system for the development of H∞ control theory (see, e.g., [14]).
The above problem is locally resolved within the framework of nonlinear H∞-control methods from [1, 14–16]. Those methods do not admit a straightforward application to the problem in question because in contrast to the standard case, a partial state stabilization (i.e., asymptotic stabilization of the output of the system) is only required provided that the complementary variables remain bounded. Their modification developed in the present paper is of the same level of simplicity, and it follows the common practice of proper solution to corresponding differential Riccati equations which is performed numerically.
The aforementioned H∞ synthesis took its origins from game-theoretic approach from Basar and Bernhard [15], and the ℒ2-gain analysis from Isidori and Astolfi [14]. It followed the line of reasoning, used in Orlov et al. [17], where the corresponding Hamilton-Jacobi-Isaacs expressions were required to be negative definite rather than semidefinite.
In contrast to the standard ℒ2-gain analysis from Isidori and Astolfi [14] and Van der Shaft [18] the resulting H∞ design procedure imposed the nonstabilizability-detectability conditions on the control systems. Under appropriate assumptions the existence of suitable solutions of Riccati differential equations, appearing in solving the H∞ control problem for the linearized system, was shown to be necessary and sufficient condition for a local solution of the H∞ control problem to exist. This mean that the verification of stabilizability and detectability conditions will be not required. A local solution was then derived by means of a certain perturbation of the Riccati equations when these unperturbed equations had bounded positive-semidefinite solutions. Thus, the local stabilizability and detectability properties of the control system were ensured by the existence of the proper solutions of the unperturbed Riccati equations, and hence the H∞ synthesis obviated any extra work on verification of these properties.
Organization of the Paper
The paper is organized as follows. Background materials on time-varying H∞-control synthesis are presented in Section 2. The tracking control problem of a 2-DOF underactuated system and its state equations are introduced in Section 3 while desired trajectory synthesis procedure is also discussed. A nonlinear H∞-output control for time varying systems is also constructed. Performance issues of this controller are illustrated in a simulation study in Section 4. Finally, Section 5 presents conclusions.
2. Background Material on Nonlinear H∞-Control of Time-Varying Systems2.1. Basic Assumptions and Problem Statement
Consider a nonlinear system of the form
ẋ=f(x,t)+g1(x,t)w+g2(x,t)u,z=h1(x,t)+k12(x,t)u,y=h2(x,t)+k21(x,t)w,
where x∈ℝn is the state space vector, t∈ℝ is the time, u∈ℝm is the control input, w∈ℝr is the unknown disturbance, z∈ℝl is the unknown output to be controlled, and y∈ℝp is the only available measurement on the system. The following assumptions are assumed to hold.
The functions f(x,t), g1(x,t), g2(x,t), h1(x,t), h2(x,t), k12(x,t), and k21(x,t) are piecewise continuous in t for all x and locally Lipschitz continuous in x for all t.
These assumptions are made for technical reasons. Assumption (A1) guarantees the well-posedness of the above dynamic system, while being enforced by integrable exogenous inputs. Assumption (A2) ensures that the origin is an equilibrium point of the nondriven (u=0) disturbance-free (w=0) dynamic system (1). Assumption (A3) is a simplifying assumption inherited from the standard H∞-control problem.
A causal dynamic feedback compensator
u=𝒦(ξ,t),ξ̇=η(ξ,y,t)
with internal state ξ∈ℝs, is said to be globally (locally) admissible controller if the closed-loop systems (1)-(2) are globally (uniformly) asymptotically stable when w=0.
Given a real number γ>0, it is said that systems (1), (2) haveℒ2-gain less than γ if the response z, resulting from w for initial state x(t0)=0, ξ(t0)=0, satisfies
∫t0t1∥z(t)∥2dt<γ2∫t0t1∥w(t)∥2dt
for all t1>t0 and all piecewise continuous functions w(t).
The time-varying H∞-control problem is to find a globally admissible controller (2)-(3) such that ℒ2-gain of the closed-loop systems (1), (2), (3) is less than γ. In turn, a locally admissible controller (2), (3) is said to be a local solution of the H∞-control problem if there exists a neighborhood U of the equilibrium such that inequality (4) is satisfied for all t1>t0 and all piecewise continuous functions w(t) for which the state trajectory of the closed-loop system starting from the initial point (x(t0),ξ(t0))=(0,0) remains in U for all t∈[t0,t1].
2.2. Local State-Space Solution
Assumptions (A1)–(A3) allow one to linearize the corresponding Hamilton-Jacobi-Isaacs inequalities from [1] that arise in the state feedback and output-injection design thereby yielding a local solution of the time-varying H∞-control problem. The subsequent local analysis involves the linear time-varying H∞-control problem for the system
ẋ=A(t)x+B1(t)w+B2(t)u,z=C1(t)x+D12(t)u,y=C2(t)x+D21(t)w,
where
A(t)=∂f∂x(0,t),B1(t)=g1(0,t),B2(t)=g2(0,t),C1(t)=∂h1∂x(0,t),D12(t)=k12(0,t),C2(t)=∂h2∂x(0,t),D21(t)=k21(0,t).
Such a problem is now well understood if the linear system (A.1) is stabilizable and detectable from u and y, respectively. Under these assumptions, the following conditions are necessary and sufficient for a solution to exist (see, e.g., [16]).
There exists a bounded positive semidefinite symmetric solution of the equation
-Ṗ(t)=P(t)A(t)+AT(t)P(t)+C1T(t)C1(t)+P(t)[1γ2B1(t)B1T(t)-B2(t)B2T(t)]P(t),
such that the system
ẋ=[A-(B2B2T-γ-2B1B1T)P](t)x(t)
is exponentially stable. (Throughout, a time-dependent n×n-matrix P(t) is positive semidefinite if and only if xTP(t)x≥0 for all n-vectors x and all time instants t whereas P(t) is positive definite if and only if xTP(t)x≥mxTx for all x and t, and some constant m>0. Respectively, P(t) is bounded if and only if ∥P(t)∥≤m0 for all t and some constant m0>0. )
There exists a bounded positive semidefinite symmetric solution to the equation
Ż(t)=Ã(t)Z(t)+Z(t)ÃT(t)+B1(t)B1T(t)+Z(t)×[1γ2P(t)B2(t)B2T(t)P(t)-C2T(t)C2(t)]Z(t)
specified with Ã(t)=A(t)+γ-2B1(t)B1T(t)P(t), such that the system
ẋ=[A-Z(C2TC2-γ-2PB2B2TP)](t)x(t)
is exponentially stable.
According to the time-varying bounded real lemma [17], conditions (C1) and (C2) ensure that there exists a positive constant ε0 such that the system of the perturbed differential Riccati equations
-Ṗε=Pε(t)A(t)+AT(t)Pε(t)+C1T(t)C1(t)+Pε(t)[1γ2B1(t)B1T(t)-B2(t)B2T(t)]Pε(t)+εI,Żε=Ã(t)Zε(t)+ZÃT(t)+B1(t)B1T(t)+Zε(t)×[1γ2P(t)B2(t)B2T(t)P(t)-C2T(t)C2(t)]Zε(t)+εI,
has a unique positive definite symmetric solution (Pε(t),Zε(t)) for each ε∈(0,ε0) where Ã(t)=A(t)+γ-2B1(t)B1T(t)Pε(t).
Differential equations (11) and (12) are subsequently utilized to derive a local solution of the nonlinear H∞-control problem for (1). The following resuls is extracted from [1].
Theorem 1.
Let conditions (C1) and (C2) be satisfied, and let (Pε(t),Zε(t)) be the corresponding positive solution of (11), (12) under some ε>0. Then the output feedback
ξ̇=f(ξ,t)+[1γ2g1(ξ,t)g1T(ξ,t)-g2(ξ,t)g2T(ξ,t)]×Pε(t)ξ+Zε(t)C2T(t)[y-h2(ξ,t)],u=-g2T(ξ,t)Pε(t)ξ
is a local solution of the H∞-control problem.
In what follows, Theorem 1 is used to design an H∞ tracking controller for the underactuated system.
3. H∞-Control of Underactuated System3.1. Problem Statement
Consider the equation of motion of an underactuated mechanical system given by the Lagrange equation
M(q)q̈+N(q,q̇)=Bτ+wx(t),
where q=[q1,q2]∈ℝ2 is a vector of generalized coordinates where q1 and q2 are the unactuated and actuated joints, respectively; τ∈ℝ is the vector of applied joint torques; B=[0,1]T is the input matrix that maps the torque input τ to the joint of coordinates space; wx(t)∈ℝ2 is the unknown disturbance vector to account for destabilizing model discrepancies due to hard-to-model nonlinear phenomena such as friction and backlash, t∈ℝ is the time; M(q)∈ℝ2×2 is the symmetric positive-definite inertia matrix; N(q,q̇)=[N1(q,q̇),N2(q,q̇)]T∈ℝ2 is the vector that contains the Coriolis, centrifugal, and gravity torques. Appendix A presents the dynamic model of the double pendulum.
The control objective is to design a nonlinear H∞ tracking controller that ensures
limt→∞∥q(t)-qd(t)∥=0
to be achieved asymptotically, while also attenuates the influence of external disturbances. Here, qd(t)∈ℝ2 is a continuously differentiable desired trajectory.
3.2. The Desired Trajectory
We point out that the present formulation is different from typical formulation of output tracking and regulation [1, 19], where the set point or the reference trajectory is a priori given because underactuated systems are not feedback or input-state linearizables due to its complexity. Therefore, special attention is required in the selection of the desired trajectory for the system under study. There are a few procedures to find desired trajectories for underactuated systems in literature [7, 8, 10, 20, 21], and under reasonable hypotheses all of them can be used to obtain a desired trajectory. The methodology from [7] is used here, where a set of exact trajectories is derived for the nonlinear equation of motion which involves inverted periodic motion. To this end, let us consider the desired trajectory which is solution of
ddt[qdq̇d]=[q̇dM-1(qd)[Bτ̂(qd,q̇d)-N(qd,q̇d)]],
where qd(t)∈ℝ2,q̇d(t)∈ℝ2 are the desired joint positions and velocities, respectively, and
τ̂=N2(qd,q̇d)-2M22(qd)-M12(qd)2M12(qd)-M11(qd)N1(qd,q̇d)
is the control input that makes the desired virtual output
yd(t)=2qd1(t)+qd2(t)-ϕ
remains at zero for all t≥0 when yd(t) starts at yd(0)=ẏd(0)=0,ϕ is a constant parameter that parameterizes the equilibrium manifold of the pendulum, and the oscillations given by (16)–(18) are around this manifold. Throughout, we confine our research interest in desired oscillations around the upright position of the pendulum which correspond to the more difficult case due that the open-loop system has an unstable zero dynamics. Toward this end, we choose ϕ=π for all t≥0 in (18). It was shown in [7] that (16) and (17) generate a set of exact periodic trajectories given by
q̈d1=c4sin(qd1)+c5sin(ϕ-qd1)c1-c2,
that can be interpreted as the zero dynamics of the system (16) with respect to the output yd(t). Time evolution of the desired trajectory is illustrated in Figure 1 where the value of ϕ is modified along the time:
ϕ={0,ift<5,0.025,if5≤t<12,0.1,if12≤t≤20,
where t∈ℝ is given in seconds. Notice that frequency and amplitude of oscillations change according to variations in ϕ. Figure 2 shows the profile of the frequency and amplitude of oscillations for several values of ϕ.
Plot of desired trajectories for Acrobot by selecting several values of ϕ.
Profile of the frequency and amplitude of oscillations for several values of ϕ.
3.3. The Task
Our objective is to design a controller of the form
τ=τ̂(qd,q̇d)+u
with internal state ξ(t)∈ℝ4, that ensures (15). Thus, the controller to be constructed consists of the trajectory compensator (17) and a disturbance attenuator u given in (2), (3) internally stabilizing the closed-loop system around the desired trajectory. In the sequel, we confine our investigation to the H∞ tracking problem, where
the output to be controlled is given by
z=ρ[02q1+q2-π]+[10]u(t)
with a positive weight coefficient ρ;
the joint position vector q(t)∈ℝ2 is the only available measurement, and this measurement is corrupted by the error vector wo(t)∈ℝ2, that is,
y=q(t)+wo(t).
The H∞ control problem in question is thus stated as follows. Given the system representation (14)–(23), the desired trajectory qd(t)∈ℝ2, and a real number γ>0, it is required to find (if any) a causal dynamic feedback controller (2), (3) such that the undisturbed closed-loop system is uniformly asymptotically stable around the origin, and its ℒ2-gain is locally less than γ, that is, inequality (4) is satisfied for all t1>t0 and all piecewise continuous functions w(t)=[wx(t),wo(t)]T for which the corresponding state trajectory of the closed-loop system, initialized at the origin, remains in some neighborhood of this point.
3.4. H∞ Synthesis
To begin with, let us introduce the state deviation vector x=(x1,x2)T∈ℝ4 where x1=(q1-qd1,q2-qd2)T and x2=(q̇1-q̇d1,q̇2-q̇d2)T. After that let us rewrite the system (14), the output to be controlled (22), and the output (23) in terms of the state vector x:
ddt[x1x2]=[x2M-1(x1+qd)N(x1+qd,x2+q̇d)-q̈d]+[0M-1(x1+qd)]wx+[0M-1(x1+qd)B]τ,z=ρ[02x11+x12+2qd1+qd2-π]+[10]u(t),y=x1+qd(t)+wo(t).
Clearly, the above H∞ tracking control problem is nothing else than a standard nonlinear H∞ control problem from [1] stated for a time-varying nonlinear system (1) specified as
f(x,t)=[x2M-1(x1+qd)N(x1+qd,x2+q̇d)-q̈d],g1(x,t)=[02×202×2M-1(x1+qd)02×2],g2(x,t)=[02×1M-1(x1+qd)B],h1(x,t)=ρ[02x11+x12+2qd1+qd2-π],h2(x,t)=x1+qd,k12(x,t)=[10],k21(x,t)=[02×2I2].
Now by applying Theorem 1 to system (1) thus specified, we derive a local solution of the H∞ tracking control problem.
Thus, the output feedback controller (13), specified according to (26), locally solves the H∞ tracking control problem (4)–(24). Stabilizability and detectability properties of the control systems are ensured by the existence of the proper solutions of the unperturbed differential Riccati equations, and hence the corresponding synthesis procedure obviates an extra work (formidable in the nonlinear case) on verification of these properties.
4. Simulation Results
The controller performance was studied in simulation by applying the exposed ideas to the Acrobot, depicted in Figure 3, which is a two-link planar robot with no actuator at the shoulder (link 1) and actuator at the elbow (link 2). In the simulation, performed with MATLAB, the Acrobot was required to move from [q1(0),q2(0)]=[-0.07,3.3] to the desired trajectory qd(t)∈ℝ2 and ϕ=π. The initial velocity q̇(0)∈ℝ2 and the initial compensator state ξ(0)∈ℝ4 were set to zero for all the simulations. The matrices M(q) and N(q,q̇) for the Acrobot are given in Appendix A. We seek for orbital stabilization of the unactuated link q2 around the equilibrium point q*=(0,π).
Schematic diagram of the acrobot where l1 and l2 denote the length of the nonactuated and actuated links, respectively; and m1 and m2 are the masses of each link.
The control goal was achieved by implementing the nonlinear H∞ controller with a weight parameter ρ=1 on the Acrobot. By iterating on γ, we found the infimal achievable level γ⋆≃250. However, in the subsequent simulations γ=2000 was selected to avoid an undesirable high-gain controller design that would appear for a value of γ close to the optimum. With γ=2000 we obtained that for ε=0.1 the corresponding differential Riccati equations (11)-(12) with
A(t)=[00100001A31(t)A32(t)A33(t)A34(t)A41(t)A42(t)A43(t)A44(t)],B1(t)=[02×202×2M-1(qd)02×2],B2(t)=[02×1M-1(qd)],C1(t)=ρ[00002100],C2(t)=[10000100],D12(t)=[10],D22(t)=[02×2I2]
have positive-definite solutions. These solutions can be numerically found with MATLAB. Matrix A(t) is given in Appendix B. It should be pointed out that the constant ϕ=π does not appear in C1(t) due to straightforward calculation of (A.2), but it is definitely required in (22) to improve the selection of γ which affects inequality (4) thus avoiding the synthesis of a high-gain controllers. Resulting trajectories is depicted in Figure 4. This figure demonstrates that the H∞ controller does asymptotically stabilize the system motion around the desired trajectory. In addition, the H∞ controller was successfully applied to the Acrobot under external disturbances
wxi(t)=bicos(0.1t)exp(-2t),i=1,2,
where b1=1×10-3 [N·m] and b2=2×10-3 [N·m] are the disturbance levels at first and second joints, respectively. Resulting trajectories are depicted in Figure 5. Figures 6 and 7 show the time evolution of the determinant of minors of matrices P(t)∈ℝ4×4 and Z(t)∈ℝ4×4 denoted as
Pm1(t)=P11(t),Pm2(t)=(P11(t)P12(t)P12(t)P22(t)),Pm3(t)=(P11(t)P12(t)P13(t)P12(t)P22(t)P23(t)P31(t)P32(t)P33(t)),Zm1(t)=Z11(t),Zm2(t)=(Z11(t)Z12(t)Z12(t)Z22(t)),Zm3(t)=(Z11(t)Z12(t)Z13(t)Z12(t)Z22(t)Z23(t)Z31(t)Z32(t)Z33(t)).
These Figures highlight that matrices P(t) and Z(t), which are solution of (11) and (12), respectively, are bounded and positive definite for all t≥0.
Phase portrait of the first joint trajectory and desired trajectory (+) for the unperturbed case (a) and perturbed case (b).
Time evolution of the output (continuous line) following the desired trajectory (dashed line) under perturbed torques.
Time evolution of the determinants of the principal minors of the matrix P∈ℝ4×4Pm1, Pm2, Pm3, Pm4.
Time evolution of the determinants of the principal minors of the matrix Z∈ℝ4×4Zm1, Zm2, Zm3, Zm4.
5. Conclusions
The output feedback Nonlinear H∞ tracking control problem is locally solved for an underactuaded mechanical system. The desired periodic orbit is centered at the upright position where the open-loop plant becomes a nonminimum-phase system. The developed controller drives the trajectories of the robot into a set of inverted exact desired trajectories governed by its zero dynamics. Simulation studies, made for the Acrobot, showed the effectiveness of the controller. The design of methods to generate reference trajectories evolving more frequencies and amplitudes in the upright position is in progress, and few results have been published for double-pendulums in [20, 22]. In future work there are two extensions of the result of the paper. First, one would like to synthesize the H∞ control taking into account reference trajectories derived from alternative methods. The other would also like to extend the result of the paper for the nonsmooth case.
AppendicesA. Dynamic Model of Acrobot
The equation motion of Acrobot, described by (14), was specified by applying the Euler-Lagrange formulation [23] where
M(q)=[M11(q)M12(q)M12(q)M22(q)],N(q,q̇)=[N1(q,q̇)N2(q,q̇)]
with
M11(q)=c1+c2-2c3cos(q2),M12(q)=c2-c3cos(q2),M22(q)=c2,N1(q,q̇)=c3sin(q2)q̇1q̇2+c3q̇2(q̇1+q̇2)sin(q2)-c4sin(q1)+c5sin(q1+q2),N2(q,q̇)=-c3sin(q2)q̇12+c5sin(q1+q2),
where the values of ci (i=1,…,5), given in Table 1, were taken from the experimental Acrobot provided in [7].
Parameter values for the Acrobot.
c1
c2
c3
c4
c5
0.0043
0.0051
0.0034
0.0494
0.038
B. Matrix A(t)
In this appendix we provide the computed matrix A(t) for the Acrobot which was used in the solution of differential Riccati equations (11) and (12):
A(t)=[00100001A31(t)A32(t)A33(t)A34(t)A41(t)A42(t)A43(t)A44(t)],
where
A31(t)=c2ΔM(qd)-1[c5cos(qd1+qd2)-c4cos(qd1)]-ΔM(qd)-1[c5cos(qd1+qd2)]M12(qd),A32(t)=c2ΔM(qd)-1[c3cos(qd2)q̇d1q̇d2+c3q̇d2(q̇d1+q̇d2)cos(qd2)+c5cos(qd1+qd2)]+-2c2c3ΔM(qd)-2sin(qd2)cos(qd2)N1(qd,q̇d)-c3ΔM(qd)-1N2(qd,q̇d)sin(qd2)+2c3ΔM(qd)-2sin(qd2)cos(qd2)M12(qd)N2(qd,q̇d)-ΔM(qd)-1[-c3cos(qd2)q̇d1+c5cos(qd1+qd2)]×M12(qd),A33(t)=2c2c3ΔM(qd)-1sin(qd2)q̇d2+2c3ΔM(qd)-1M12(qd)sin(qd2)q̇d1,A34(t)=2c2c3ΔM(qd)-1[sin(qd2)q̇d1+sin(qd2)q̇d2],A41(t)=ΔM(qd)-1[c5cos(qd1+qd2)]M11(qd)-ΔM(qd)-1[c5cos(qd1+qd2)-c4cos(qd1)]×M12(qd),A42(t)=2c3ΔM(qd)-1N2(qd,q̇d)sin(qd2)+ΔM(qd)-1×[c5cos(qd1+qd2)-c3cos(qd2)q̇d12]M11(qd)-2c3ΔM(qd)-2M11(qd)N2(qd,q̇d)sin(qd2)cos(qd2)-c3ΔM(qd)-1N1(qd,q̇d)sin(qd2)-ΔM(qd)-1×[2c3cos(qd2)q̇d1q̇d2+c3cos(qd2)q̇d22+c5cos(qd1+qd2)]M12(qd)+2c3ΔM(qd)-2M12(qd)N1(qd,q̇d)sin(qd2)cos(qd2),A43(t)=-2c3ΔM(qd)-1[M11(qd)+M12(qd)]sin(qd2)q̇d1,A44(t)=-2c3ΔM(qd)-1[q̇d1+q̇d2]sin(qd2)M12(qd),
where ΔM(qd)=M11(qd)M22(qd)-M122(qd).
AguilarL. T.OrlovY.AchoL.Nonlinear H∞-control of nonsmooth time-varying systems with application to friction mechanical manipulators20033991531154210.1016/S0005-1098(03)00148-1WesterveltE. R.GrizzleJ.ChevallereauC.ChoJ.MorrisB.2007London, UKTaylor & Francis/CRCM'CloskeyR. T.BurdickJ. W.VakakisA. F.On the periodic motions of simple hopping robotsProceedings of the IEEE International Conference on Systems, Man and Cybernetics (ICSMC '90)November 1990Los Angeles, Calif, USA77177710.1109/ICSMC.1990.142225AguilarL. T.OrlovY.CadiouJ. C.MerzoukiR.Nonlinear H∞-output regulation of a nonminimum phase servomechanism with backlash2007129454454910.1115/1.2719774BrogliatoB.19992ndLondon, UKSpringerBrogliatoB.RioA. Z.On the control of complementary-slackness juggling mechanical systems200045223524610.1109/9.839946BerkemeierM. D.FearingR. S.Tracking fast inverted trajectories of the underactuated Acrobot199915474075010.1109/70.782028GrizzleJ. W.MoogC. H.ChevallereauC.Nonlinear control of mechanical systems with an unactuated cyclic variable200550555957610.1109/TAC.2005.847057OrlovY.AguilarL. T.AchoL.Model orbit robust stabilization (MORS) of pendubot with application to swing up control6Proceedings of the 44th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC '05)December 2005Seville, Spain61646169ShiriaevA.PerramJ. W.Canudas-de-WitC.Constructive tool for orbital stabilization of underactuated nonlinear systems: virtual constraints approach20055081164117610.1109/TAC.2005.852568ConsoliniL.TosquesM.A morphing method for exact tracking control of nonminimum phase non linear systemsProceedings of the 46th IEEE Conference on Decision and Control (CDC '07)December 2007New Orleans, La, USA5240524510.1109/CDC.2007.4434826WangW.Output feedback control of unstable non-minimum phase systemsProceedings of the 25th Chinese Control Conference (CHICC '06)August 2006Harbin, China78178410.1109/CHICC.2006.280740IsidoriA.A tool for semiglobal stabilization of uncertain non-minimum-phase nonlinear systems via output feedback200045101817182710.1109/TAC.2000.880972IsidoriA.AstolfiA.Disturbance attenuation and H∞-control via measurement feedback in nonlinear systems19923791283129310.1109/9.159566BasarT.BernhardP.1990Boston, Mass, USABirkhäuserRaviR.NagpalK. M.KhargonekarP. P.H∞ control of linear time-varying systems: a state-space approach19912961394141310.1137/0329071OrlovY.AchoL.SolisV.Nonlinear H∞-control of time-varying systems4Proceedings of the 38th IEEE Conference on Decision and Control (CDC '99)December 1999Phoenix, Ariz, USA3764376910.1109/CDC.1999.827940van der ShaftA. J.L2-gain analysis of nonlinear systems and nonlinear state feedback H∞ control199237677078410.1109/9.256331OrlovY.yorlov@cicese.mxAguilarL. T.Non-smooth H∞-position control of mechanical manipulators with frictional joints200477111062106910.1080/0020717042000273087AguilarL. T.laguilar@citedi.mxBoikoI.i.boiko@ieee.orgFridmanL.lfridman@servidor.unam.mxIriarteR.rafael@servidor.unam.mxOutput excitation via continuous sliding-modes to generate periodic motion in underactuated systemsProceedings of the 45th IEEE Conference on Decision and Control (CDC '06)December 2006San Diego, Calif, USA1629163410.1109/CDC.2006.377713FreidovichL.RobertssonA.ShiriaevA.JohanssonR.Periodic motions of the Pendubot via virtual holonomic constraints: theory and experiments200844378579110.1016/j.automatica.2007.07.011AguilarL. T.BoikoI.FridmanL.IriarteR.Periodic motion of underactuated mechanical systems self-generated by variable structure controllers: design and experimentsProceedings of the European Control Conference (ECC '07)July 2007Kos, Greece37963801CraigJ.1989Reading, Mass, USAAddison-Wesley