V(x)
x\ne0
x=0
u(x,t)
The theory and application of control-Lyapunov functions were developed by Zvi Artstein and Eduardo D. Sontag in the 1980s and 1990s.
Consider an autonomous dynamical system with inputswhere
x\inRn
u\inRm
x*\inRn
D\subsetRn
x*=0
x* ≠ 0
Definition. A control-Lyapunov function (CLF) is a function
V:D\toR
V(x)
x\inD
x=0
x\inRn(x ≠ 0),
u\inRm
V |
(x,u):=\langle\nablaV(x),f(x,u)\rangle<0,
\langleu,v\rangle
u,v\inRn
The last condition is the key condition; in words it says that for each state x we can find a control u that will reduce the "energy" V. Intuitively, if in each state we can always find a way to reduce the energy, we should eventually be able to bring the energy asymptotically to zero, that is to bring the system to a stop. This is made rigorous by Artstein's theorem.
Some results apply only to control-affine systems—i.e., control systems in the following form:where
f:Rn\toRn
gi:Rn\toRn
i=1,...,m
Eduardo Sontag showed that for a given control system, there exists a continuous CLF if and only if the origin is asymptotic stabilizable.[5] It was later shown by Francis H. Clarke, Yuri Ledyaev, Eduardo Sontag, and A.I. Subbotin that every asymptotically controllable system can be stabilized by a (generally discontinuous) feedback.[6] Artstein proved that the dynamical system has a differentiable control-Lyapunov function if and only if there exists a regular stabilizing feedback u(x).
It is often difficult to find a control-Lyapunov function for a given system, but if one is found, then the feedback stabilization problem simplifies considerably. For the control affine system, Sontag's formula (or Sontag's universal formula) gives the feedback law
k:Rn\toRm
(m=1)
k(x)=\begin{cases}\displaystyle-
LfV(x)+\sqrt{\left[LfV(x)\right]2+\left[LgV(x)\right]4 | |
LfV(x):=\langle\nablaV(x),f(x)\rangle
LgV(x):=\langle\nablaV(x),g(x)\rangle
V
f
g
For the general nonlinear system, the input
u
u*(x)=\underset{u}{\operatorname{argmin}}\nablaV(x) ⋅ f(x,u)
Here is a characteristic example of applying a Lyapunov candidate function to a control problem.
Consider the non-linear system, which is a mass-spring-damper system with spring hardening and position dependent mass described by
| |||
m(1+q |
+K0q+K
3=u | |
1q |
qd
q
e=qd-q
r
r=e |
+\alphae
r\mapstoV(r):=
1 | |
2 |
r2
r\ne0
Now taking the time derivative of
V
V | =r |
r |
V | =( |
e |
+\alphae)(\ddot{e}+\alpha
e |
)
The goal is to get the time derivative to be
V |
=-\kappaV
V
Hence we want the rightmost bracket of
V |
(\ddot{e}+\alpha
e |
)=(\ddot{q}d-\ddot{q}+\alpha
e |
)
(\ddot{q}d-\ddot{q}+\alpha
e |
)=-
\kappa | ( | |
2 |
e |
+\alphae)
\ddot{q}
\left(\ddot{q} | ||||||||||||||||||
|
+\alpha
e |
\right)=-
\kappa | ( | |
2 |
e |
+\alphae)
u
u=
2)\left(\ddot{q} | |
m(1+q | |
d |
+\alpha
e | + |
\kappa | |
2 |
r\right)+K0q+K
| |||
1q |
\kappa
\alpha
This control law will guarantee global exponential stability since upon substitution into the time derivative yields, as expected
V |
=-\kappaV
V=V(0)\exp(-\kappat)
And hence the error and error rate, remembering that
V= | 1 | ( |
2 |
e |
+\alphae)2
If you wish to tune a particular response from this, it is necessary to substitute back into the solution we derived for
V
e
rr | =- |
\kappa | |
2 |
r2
r | =- |
\kappa | |
2 |
r
r=r(0)\exp\left(- | \kappa |
2 |
t\right)
e |
+\alphae=(
e |
(0)+\alphae(0))\exp\left(-
\kappa | |
2 |
t\right)