Lasalle invariance principle for ordinary dierential equations and applications
Table Of Contents
- Acknowledgment i
Certication ii
Approval iii
Introduction v
Dedication vi
1 Preliminaries 2
- 1.1Denitions and basic Theorems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
- 1.2Exponential of matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2 Basic Theory of Ordinary Dierential Equations 7
- 2.1Denitions and basic properties . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
- 2.2Continuous dependence with respect to the initial conditions . . . . . . . . . . . . 11
- 2.3Local existence and blowing up phenomena for ODEs . . . . . . . . . . . . . . . . 12
- 2.4Variation of constants formula . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3 Stability via linearization principle 21
- 3.1Denitions and basic results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
4 Lyapunov functions and LaSalle’s invariance principle 26
- 4.1Denitions and basic results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
- 4.2Instability Theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
- 4.3How to search for a Lyapunov function (variable gradient method) . . . . . . . . . 30
- 4.4LaSalle’s invariance principle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
- 4.5Barbashin and Krasorskii Corollaries . . . . . . . . . . . . . . . . . . . . . . . . . . 32
- 4.6Linear systems and linearization . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
5 More applications 40
- 5.1Control design based on lyapunov’s direct method . . . . . . . . . . . . . . . . . . 41
Conclusion 48
Bibliography 48
Thesis Abstract
Abstract
The LaSalle invariance principle is a fundamental tool in the analysis of stability for systems described by ordinary differential equations. This principle provides a powerful method for proving the asymptotic stability of equilibrium points or limit cycles. The essence of the LaSalle invariance principle lies in constructing a Lyapunov function that decreases along solutions of the differential equation, eventually leading to the convergence of trajectories to a region called the largest invariant set. This region typically corresponds to the set of points where the derivative of the Lyapunov function is zero. The LaSalle invariance principle has a wide range of applications in various fields, including control theory, dynamical systems, and mathematical biology. In this research project, we delve into the theoretical foundations of the LaSalle invariance principle and explore its applications in the context of ordinary differential equations. We study the conditions under which the LaSalle invariance principle can be applied to guarantee stability properties of solutions to differential equations. By analyzing the behavior of trajectories near equilibrium points or limit cycles, we aim to establish the asymptotic stability of these critical points using Lyapunov functions constructed based on the LaSalle invariance principle. Furthermore, we investigate the robustness of the LaSalle invariance principle with respect to perturbations in the system dynamics. Understanding how small changes in the system parameters affect the stability analysis is crucial for practical applications of the principle. We explore techniques for incorporating uncertainties into the stability analysis framework and develop methods for quantifying the impact of perturbations on the stability properties of the system. Moreover, we consider practical examples to illustrate the effectiveness of the LaSalle invariance principle in analyzing stability properties of real-world systems. By applying the principle to specific differential equations modeling physical, biological, or engineering systems, we demonstrate how it can be used to predict the long-term behavior of trajectories and assess the stability of the system under different operating conditions. Overall, this research project contributes to the understanding and application of the LaSalle invariance principle for ordinary differential equations, shedding light on its significance in stability analysis and its relevance to a wide range of scientific and engineering disciplines.
Thesis Overview
<p>
</p><p>PRELIMINARIES<br>1.1 Denitions and basic Theorems<br>In this chapter, we focussed on the basic concepts of the ordinary dierential equations. Also, we<br>emphasized on relevant theroems in ordinary dierential equations.<br>Denition 1.1.1 An equation containing only ordinary derivatives of one or more dependent vari-<br>ables with respect to a single independent variable is called an ordinary dierential equation ODE.<br>The order of an ODE is the order of the highest derivative in the equation. In symbol, we can<br>express an n-th order ODE by the form<br>x(n) = f(t; x; :::; x(nô€€€1)) (1.1.1)<br>Denition 1.1.2 (Autonomous ODE ) When f is time-independent, then (1.1.1) is said to be<br>an autonomous ODE. For example,<br>x0(t) = sin(x(t))<br>Denition 1.1.3 (Non-autonomous ODE ) When f is time-dependent, then (1.1.1) is said to<br>be a non autonomous ODE. For example,<br>x0(t) = (1 + t2)y2(t)<br>Denition 1.1.4 f : Rn ! Rn is said to be locally Lipschitz, if for all r > 0 there exists k(r) > 0<br>such that<br>kf(x) ô€€€ f(y)k k(r)kx ô€€€ yk; for all x; y 2 B(0; r):<br>f : Rn ! Rn is said to be Lipschitz, if there exists k > 0 such that<br>kf(x) ô€€€ f(y)k kkx ô€€€ yk; for all x; y 2 Rn:<br>Denition 1.1.5 (Initial value problem (IVP) Let I be an interval containing x0, the follow-<br>ing problem (<br>x(n)(t) = f(t; x(t); :::; x(nô€€€1)(t))<br>x(t0) = x0; x0(t0) = x1; :::; x(nô€€€1)(t0) = xnô€€€1<br>(1.1.2)<br>is called an initial value problem (IVP).<br>x(t0) = x0; x0(t0) = x1; :::; x(nô€€€1)(t0) = xnô€€€1<br>are called initial condition.<br>2<br>Lemma 1.1.6 [9](Gronwall’s Lemma) Let u; v : [a; b] ! R+ be continuous such that there exists<br>> 0 such that<br>u(x) +</p><p>x<br>a<br>u(s)v(s)ds; for all x 2 [a; b]:<br>Then,<br>u(x) e</p><p>x<br>a<br>v(s)ds<br>; for all x 2 [a; b]:<br>Proof .<br>u(x) +</p><p>x<br>a<br>u(s)v(s)ds<br>implies that<br>u(x)<br>+</p><p>x<br>a<br>u(s)v(s)ds<br>v(x):<br>So,<br>u(x)v(x)<br>+</p><p>x<br>a<br>u(s)v(s)ds<br>v(x);<br>which implies that<br>x<br>a<br>u(x)v(x)<br>+</p><p>x<br>a<br>u(s)v(s)ds<br>ds</p><p>x<br>a<br>v(x)ds:<br>So, taking exponential of both side we get<br>u(x) +</p><p>x<br>a<br>u(s)v(s)ds</p><p>x<br>a<br>u(s)v(s)ds:<br>Thus,<br>u(x) e</p><p>x<br>a<br>v(s)ds<br>; x 2 [a; b]:<br>Corollary 1.1.7 Let u; v : [a; b] ! R+ be continuous such that<br>u(x)</p><p>x<br>a<br>u(s)v(s)ds; for all x 2 [a; b]:<br>Then, u = 0 on [a; b].<br>Proof . Now,<br>u(x)</p><p>x<br>a<br>u(s)v(s)ds<br>implies that<br>u(x)</p><p>x<br>a<br>u(s)v(s)ds<br>1<br>n<br>+</p><p>x<br>a<br>u(s)v(s)ds; for all n 1:<br>So, by Gronwall’s lemma,<br>u(x)<br>1<br>n<br>e</p><p>x<br>a u(s)v(s)ds;<br>so as<br>n ! 1; u(x) ! 0:<br>Thus, u(x) = 0, since u(x) 0. Hence, u = 0 on [a; b].<br>3<br>1.2 Exponential of matrices<br>Denition 1.2.1 Let A 2 Mnn(R), then eA is an n n matrix given by the power series<br>eA =<br>1X<br>k=0<br>Ak<br>k!<br>The series above converges absolutely for all A 2 Mnn(R)<br>Proof . The n-th partial sum is<br>Sn =<br>Xn<br>k=0<br>Ak<br>k!<br>So, let n > m Then,<br>Sn ô€€€ Sm =<br>Xn<br>k=m+1<br>Ak<br>k!<br>:<br>So,<br>kSn ô€€€ Smk<br>Xn<br>k=m+1<br>kAkk<br>k!<br>:<br>So as<br>m ! 1; kSn ô€€€ Smk ! 0<br>So, (Sn)n is Cauchy. Thus, converges.<br>Theorem 1.2.2 [3](Cayley Hamilton Theorem)<br>Let A 2 Mnn(R) and () = det(I ô€€€ A) its characteristic polynomial then<br>(A) = 0:<br>Proof . Let A 2 Mnn(R);<br>() = det(I ô€€€ A) = c0 + c1 + c22 + ::: + cnn:<br>adj(A ô€€€ I) = B0 + B1 + B22 + ::: + Bnô€€€2nô€€€2 + Bnô€€€1nô€€€1;<br>where Bi 2 Mnn(R) for i = 0; 1; 2; :::; n; but, from linear algebra we have that<br>Aô€€€1 =<br>adj(A)<br>det(A)<br>;<br>where adj(A) denotes the adjugate or classical adjoint of A. So,<br>det(I ô€€€ tA)I = (I ô€€€ tA)adj(I ô€€€ tA):<br>(A ô€€€ I)(B0 + B1 + B22 + ::: + Bnô€€€2nô€€€2 + Bnô€€€1nô€€€1) = (c0 + c1 + c22 + ::: + cnn)I:<br>Observe that the entries in adj(I ô€€€tA) are polynomials in of degree at most nô€€€1. So, Bi is the<br>zero matrix for i = n. Equating the coecients of n on both sides gives<br>c0I + c1A + c2A2 + ::: + cnAn = 0:<br>Thus,<br>(A) = 0:<br>4<br>Example 1.2.3 (Application of Cayley Hamilton Theorem)<br>Find etA for A =</p><p>0 1<br>ô€€€1 0<br>!<br>Solution:<br>The characteristic equation is s2 + 1 = 0, and the eigenvalues are 1 = i, and 2 = ô€€€i. So, by<br>Theorem 1.2.2 we have that,<br>etA = 0I + 1A;<br>where we are to nd the value of 0, and 1. So,<br>eti = cos t + i sin t = 0 + 1i<br>eô€€€ti = cos t ô€€€ i sin t = 0 ô€€€ 1i<br>which implies that 0 = cos t, and 1 = sin t. So,<br>etA = cos(t)I + sin(t)A =</p><p>cos t sin t<br>ô€€€sin t cos t<br>!<br>Theorem 1.2.4 [11] Let A;B 2 Mnn(R). Then,<br>(1) If 0 denotes the zero matrix, then e0 = I, the identity matrix.<br>(2) If A is invertible, then eABAô€€€1<br>= AeBAô€€€1.<br>Proof . Recall that, for all integers s 0, we have (ABAô€€€1)s = ABsAô€€€1. Now,<br>eABAô€€€1<br>= I + ABAô€€€1 +<br>(ABAô€€€1)2<br>2!<br>+ :::<br>= I + ABAô€€€1 +<br>AB2Aô€€€1<br>2!<br>+ :::<br>= A(I + B +<br>B2<br>2!<br>+ :::)Aô€€€1<br>= AeBAô€€€1:<br>(3) If A is symmetric such that A = AT , then<br>e(AT ) = (eA)T :<br>Proof .<br>eA =<br>1X<br>k=0<br>Ak<br>k!<br>:<br>Then<br>eAT<br>=<br>1X<br>k=0<br>(AT )k<br>k!<br>=<br>1X<br>k=0<br>(Ak)T<br>k!<br>= (<br>1X<br>k=0<br>Ak<br>k!<br>)T = (eA)T :<br>(4) If AB = BA, then<br>eA+B = eAeB:<br>Proof .<br>eAeB = (I + A +<br>A2<br>2!<br>+<br>A3<br>3!<br>+ :::)(I + B +<br>B2<br>2!<br>+<br>B3<br>3!<br>+ :::)<br>= (<br>1X<br>k=0<br>Ak<br>k!<br>)(<br>1X<br>j=0<br>Bj<br>j!<br>)<br>=<br>1X<br>k=0<br>1X<br>j=0<br>(A + B)k+j<br>j!k!<br>5<br>Put m = j + k, then j = m ô€€€ k then from the binomial theorem that<br>eAeB =<br>1X<br>m=0<br>1X<br>k=0<br>AmBmô€€€k<br>(m ô€€€ k)!k!<br>=<br>1X<br>m=0<br>Am<br>m!<br>1X<br>k=0<br>m!<br>(m ô€€€ k)!<br>Bmô€€€k<br>k!<br>=<br>1X<br>m=0<br>(A + B)m<br>m!<br>= eA+B:<br>Theorem 1.2.5 [9]<br>detA<br>dt<br>= AetA = etAA; for t 2 R:<br>Proof . x(t; x0) = etAx0. Then,<br>dx(t; x0)<br>dt<br>= etAx0A =<br>1X<br>k=0<br>tkAk<br>k!<br>x0A = ( lim<br>n!1<br>Xn<br>k=0<br>tkAk<br>k!<br>)x0A<br>= lim<br>n!1<br>Xn<br>k=0<br>tkAk+1<br>k!<br>x0 = lim<br>n!1<br>Xn<br>k=0<br>AtkAk<br>k!<br>x0 = A<br>1X<br>k=0<br>tkAk<br>k!<br>x0 = AetAx0<br>So,<br>detA<br>dt<br>= AetA = etAA:<br>Proposition 1.2.6 The solution x(:; x0) of the following linear space<br>(<br>x0(t) = Ax(t); t 2 R<br>x(0) = x0 2 Rn<br>where A 2 Mnn(R), is given by<br>x(t; x0) = etAx0:<br>6</p>
<br><p></p>