NONLINEAR CONTROLLER COMPARISON ON A BENCHMARK ...
NONLINEAR CONTROLLER COMPARISON ON A BENCHMARK ...
NONLINEAR CONTROLLER COMPARISON ON A BENCHMARK ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
a total closed-loop expression, provided that the plant input, up, is de ned according<br />
to the following interconnection constraint:<br />
up = ; @Vc(qcqp)<br />
: (3.13)<br />
@qp<br />
Thus, the second and most important factabout EL systems is that the closed-loop<br />
system is itself an EL system, and the closed-loop EL parameters (TclVclFcl) are<br />
simply the sum of the parameters of the plant and the control. In other words,<br />
we can shape the energy and dissipation of the closed loop system as desired by<br />
choosing the dynamics of the control in the correct manner: it must comply with<br />
the aforementioned interconnection constraint. In [10] this idea is further re ned<br />
by showing that the injected dissipation function in the control need not necessarily<br />
be a function of the derivative of the plant variables. Rather, a dynamic system in<br />
the control is su cient, under certain passivity conditions on the original plant, to<br />
guarantee closed loop asymptotic stability. In addition, Ortega et.al. [10] add the<br />
constraint thatthere is limited actuator power available, or in other words<br />
up umax: (3.14)<br />
Following the design outlined in [15], a suitable Rayleigh dissipation function Fc was<br />
chosen:<br />
Fc = 1<br />
2ab _<br />
qc<br />
(3.15)<br />
where a and b are parameters that scale the injected damping. Also, the potential<br />
energy function of the control was constructed so as to ensure the stability of the zero<br />
equilibrium point. A su cient condition is that Vcl = Vp + Vc have a global minimum<br />
at zero. Since the potential energy of the plant is a quadratic function of the states,<br />
we chose Vc as<br />
Vc(q3qc) = 1<br />
Z Z qc+bq3<br />
q3<br />
k1 tanh(s)ds + k2 tanh(s)ds: (3.16)<br />
b 0<br />
0<br />
so that it was the sum of two strictly increasing integrals that are minimized at the<br />
zero state. Here tanh(s) was chosen as a saturation function to model the actuator<br />
20