TheoremComplete

LaSalle's Invariance Principle

LaSalle's invariance principle extends Lyapunov's method to handle cases where VΛ™\dot{V} is only negative semi-definite, by analyzing the invariant structure of the set where VΛ™=0\dot{V} = 0.


Statement

Theorem7.9LaSalle's invariance principle

Let Ξ©βŠ‚Rn\Omega \subset \mathbb{R}^n be a compact set that is positively invariant for xβ€²=F(x)\mathbf{x}' = \mathbf{F}(\mathbf{x}) (i.e., trajectories starting in Ξ©\Omega remain in Ξ©\Omega for all tβ‰₯0t \geq 0). Let V:Ξ©β†’RV: \Omega \to \mathbb{R} be C1C^1 with VΛ™(x)≀0\dot{V}(\mathbf{x}) \leq 0 on Ξ©\Omega. Define E={x∈Ω:VΛ™(x)=0}E = \{\mathbf{x} \in \Omega : \dot{V}(\mathbf{x}) = 0\} and let MM be the largest invariant subset of EE.

Then every solution starting in Ξ©\Omega approaches MM as tβ†’βˆžt \to \infty:

lim⁑tβ†’βˆžd(x(t),M)=0.\lim_{t \to \infty} d(\mathbf{x}(t), M) = 0.


Proof

Proof

Since VV is non-increasing along trajectories and bounded below on Ξ©\Omega (compact), V(x(t))V(\mathbf{x}(t)) converges to some limit cc as tβ†’βˆžt \to \infty.

The Ο‰\omega-limit set Ο‰(x)\omega(\mathbf{x}) is nonempty (by compactness), contained in Ξ©\Omega, and invariant. For any yβˆˆΟ‰(x)\mathbf{y} \in \omega(\mathbf{x}), there exists tnβ†’βˆžt_n \to \infty with x(tn)β†’y\mathbf{x}(t_n) \to \mathbf{y}, so V(y)=lim⁑V(x(tn))=cV(\mathbf{y}) = \lim V(\mathbf{x}(t_n)) = c.

Therefore VV is constant (=c= c) on Ο‰(x)\omega(\mathbf{x}). Since Ο‰(x)\omega(\mathbf{x}) is invariant and VV is constant there, VΛ™=0\dot{V} = 0 on Ο‰(x)\omega(\mathbf{x}). Thus Ο‰(x)βŠ‚E\omega(\mathbf{x}) \subset E. Since Ο‰(x)\omega(\mathbf{x}) is invariant, Ο‰(x)βŠ‚M\omega(\mathbf{x}) \subset M (the largest invariant subset of EE). β– \blacksquare

β– 

Applications

ExampleDamped oscillator

xβ€²β€²+bxβ€²+kx=0x'' + bx' + kx = 0 (b,k>0b, k > 0). As a system: xβ€²=(y,βˆ’kxβˆ’by)T\mathbf{x}' = (y, -kx-by)^T. Use V=12(kx2+y2)V = \frac{1}{2}(kx^2 + y^2) (energy).

VΛ™=kxy+y(βˆ’kxβˆ’by)=βˆ’by2≀0\dot{V} = kxy + y(-kx-by) = -by^2 \leq 0.

E={y=0}E = \{y = 0\}. On EE: xβ€²=y=0x' = y = 0 (so x=constx = \text{const}) and yβ€²=βˆ’kxy' = -kx (so x=0x = 0 for invariance). Therefore M={(0,0)}M = \{(0,0)\}.

By LaSalle: all trajectories approach (0,0)(0,0). Asymptotic stability is established even though VΛ™\dot{V} is only semi-definite.

ExampleNonlinear damping

x′′+g(x′)+h(x)=0x'' + g(x') + h(x) = 0 where h(0)=0h(0) = 0, xh(x)>0xh(x) > 0 for x≠0x \neq 0, and g(0)=0g(0) = 0, yg(y)>0yg(y) > 0 for y≠0y \neq 0.

V=∫0xh(ΞΎ) dΞΎ+12y2>0V = \int_0^x h(\xi)\,d\xi + \frac{1}{2}y^2 > 0. Then VΛ™=h(x)y+y(βˆ’g(y)βˆ’h(x))=βˆ’yg(y)≀0\dot{V} = h(x)y + y(-g(y)-h(x)) = -yg(y) \leq 0.

E={y=0}E = \{y = 0\}, and on EE: yβ€²=βˆ’h(x)y' = -h(x), so invariance requires h(x)=0h(x) = 0, hence x=0x = 0. So M={(0,0)}M = \{(0,0)\} and the origin is asymptotically stable.

RemarkComparison with strict Lyapunov

Finding a strict Lyapunov function (VΛ™<0\dot{V} < 0) can be very difficult. LaSalle's principle allows the use of natural energy functions (which typically satisfy only V˙≀0\dot{V} \leq 0) and compensates by analyzing the invariant set structure. This makes it far more practical than requiring strict negative definiteness.