MATHEMATICAL MODELS - Vol. I - Controllability, Observability, and Stability of Mathematical Models - Abderrahman Iggidr

Similar documents
Stability Analysis of the Euler Discretization for SIR Epidemic Model

Castiel, Supernatural, Season 6, Episode 18

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Existence of viscosity solutions with asymptotic behavior of exterior problems for Hessian equations

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

THE SI AND SIR EPIDEMICS ON GENERAL NETWORKS DAVID A L D O U S (BERKELEY)

K. Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Appendix: The Laplace Transform

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

MATHEMATICS. The assessment objectives of the Compulsory Part are to test the candidates :

Infinite Sequences and Series

Chapter 4 : Laplace Transform

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Stability of fractional positive nonlinear systems

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

Theory of Ordinary Differential Equations. Stability and Bifurcation II. John A. Burns

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

Axioms of Measure Theory

ECE 308 Discrete-Time Signals and Systems

Principle Of Superposition

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Chapter 10: Power Series

Properties of Fuzzy Length on Fuzzy Set

Numerical Method for Blasius Equation on an infinite Interval

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

MAT1026 Calculus II Basic Convergence Tests for Series

Generating Functions for Laguerre Type Polynomials. Group Theoretic method

DISTRIBUTION LAW Okunev I.V.

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. XII - Feedback Linearization of Nonlinear Systems - Alberto Isidori and Claudio De Persis

Simulation of Discrete Event Systems

Estimation for Complete Data

CHAPTER 10 INFINITE SEQUENCES AND SERIES

Riesz-Fischer Sequences and Lower Frame Bounds

(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3

7 Sequences of real numbers

Subject: Differential Equations & Mathematical Modeling -III. Lesson: Power series solutions of Differential Equations. about ordinary points

Sequences and Series of Functions

An Algorithm for Maximizing the Controllable Set for Open-Loop Unstable Systems under Input Saturation

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Course Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.

A Note on Matrix Rigidity

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Section 5.5. Infinite Series: The Ratio Test

Continuous Functions

sin(n) + 2 cos(2n) n 3/2 3 sin(n) 2cos(2n) n 3/2 a n =

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Decoupling Zeros of Positive Discrete-Time Linear Systems*

Linear Regression Demystified

CDS 101: Lecture 5.1 Controllability and State Space Feedback

Scientific Research of the Institute of Mathematics and Computer Science

Zeros of Polynomials

Several properties of new ellipsoids

THE KALMAN FILTER RAUL ROJAS

Lyapunov Stability Analysis for Feedback Control Design

CHAPTER 5. Theory and Solution Using Matrix Techniques

Ω ). Then the following inequality takes place:

The Differential Transform Method for Solving Volterra s Population Model

You may work in pairs or purely individually for this assignment.

Section 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations

lim za n n = z lim a n n.

Math 312 Lecture Notes One Dimensional Maps

FUZZY ALTERNATING DIRECTION IMPLICIT METHOD FOR SOLVING PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS IN THREE DIMENSIONS

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients.

Dynamics of Piecewise Continuous Functions

MATHEMATICS. The assessment objectives of the Compulsory Part are to test the candidates :

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Math 155 (Lecture 3)

BETWEEN QUASICONVEX AND CONVEX SET-VALUED MAPPINGS. 1. Introduction. Throughout the paper we denote by X a linear space and by Y a topological linear

6.3 Testing Series With Positive Terms

TEACHER CERTIFICATION STUDY GUIDE

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4

6.003 Homework #3 Solutions

Product measures, Tonelli s and Fubini s theorems For use in MAT3400/4400, autumn 2014 Nadia S. Larsen. Version of 13 October 2014.

3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials

Markscheme May 2015 Calculus Higher level Paper 3

ECON 3150/4150, Spring term Lecture 3

Discrete-Time Signals and Systems. Discrete-Time Signals and Systems. Signal Symmetry. Elementary Discrete-Time Signals.

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

CHAPTER I: Vector Spaces

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

TENSOR PRODUCTS AND PARTIAL TRACES

Summary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram.

Expectation and Variance of a random variable

A NOTE ON BOUNDARY BLOW-UP PROBLEM OF u = u p

PRELIM PROBLEM SOLUTIONS

Algebra of Least Squares

A Negative Result. We consider the resolvent problem for the scalar Oseen equation

MAS111 Convergence and Continuity

Modified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations

A Note on the Symmetric Powers of the Standard Representation of S n

OPTIMAL STOPPING AND EXIT TIMES FOR SOME CLASSES OF RANDOM PROCESSES. Vladyslav Tomashyk

Question1 Multiple choices (circle the most appropriate one):

Teaching Mathematics Concepts via Computer Algebra Systems

CHAPTER 3 A REFINED MODEL FOR GONORRHEA DYNAMICS

b i u x i U a i j u x i u x j

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

Kinetics of Complex Reactions

Transcription:

CONTROLLABILITY, OBSERVABILITY, AND STABILITY OF MATHEMATICAL MODELS Abderrahma Iggidr INRIA (Ur Lorraie) ad, Uiversity of Metz, Frace Keywords: accessibility, asymptotic stability, attractivity, chemostat, closed-loop, cotrol-lability, differetial equatio, fiite dimesioal systems, iput, Lie algebra, Lotka-Volterra systems, Lyapuov fuctios, oliear systems, observability, observer, output, stabilizatio, state. Cotets 1. Itroductio 2. Cotrollability 2.1. What is Cotrollability? 2.2. Cotrollability of Liear Systems 2.3. Accessibility Criteria for Noliear Systems 3. Stability 3.1. Geeral Defiitios 3.2. Stability of Liear Systems 3.3.Liearizatio ad Stability of Noliear Systems 3.4. Lyapuov Fuctios 3.5. Limit Cycle 3.6. Stabilizatio 3.6.1. Sufficiet Stabilizability Coditios 4. Observability 4.1. Observability of Liear Systems 4.2. Observability of Noliear Systems 4.3. Examples from Life Support Systems 5. Observers 5.1. Observers for Liear Systems 5.2. Some Noliear Observers 5.2.1. System with o Iput 5.2.2. Affie Cotrol Systems 5.2.3. Observers for a Class of No-affie Cotrol Systems 5.2.4. Asymptotic Observers Ackowledgemets Glossary Bibliography Biographical Sketch Summary This chapter presets a overview of three fudametal cocepts i Mathematical System Theory: cotrollability, stability ad observability. These properties play a promiet role i the study of mathematical models ad i the uderstadig of their behavior. They costitute the mai research subject i Cotrol Theory. Historically the

tools ad techiques of Automatic Cotrol have bee developed for artificial egieerig systems but owadays they are more ad more applied to atural systems. The mai objective of this chapter is to show how these tools ca be helpful to model ad to cotrol a wide variety of atural systems. 1. Itroductio The mai goal of this chapter is to develop i some details some otios of Cotrol Theory itroduced i Basic priciples of mathematical modelig. It cocers more specifically three structural properties of cotrol systems: cotrollability, stability ad observability. Based o the refereces listed i the ed of this chapter, we give a survey of these three properties with applicatios to various LSS examples. By a cotrol system we mea dyamical system evolvig i some state space ad that ca be cotrolled by the user. More precisely we are iterested i the study of systems that ca be modeled by differetial (respectively differece) equatios of the form dx x () t = = X ( x() t, u() t ), y() t = h( x() t ), where the variable t represets the time, the vector x(t) is the state of the system at time t, the vector u(t) is the iput or the cotrol, i.e., the actio of the user or of the eviromet ad the vector y() t is the output of the system that is, the available iformatio that ca be measured or observed by the user. The dyamics fuctio X idicates how the system chages over time. We shall address the followig problems: Whe the whole state x(t) is ot available for measuremet, how is it possible to use the iformatio provided by y( t ) together with the dyamics (1) i order to get a good estimatio of the real state of the system? This turs out to be a observability problem. How to use the cotrol u(t) i order to meet some specified eeds? This is the problem of cotrollability ad stabilizability. (1) To illustrate these various cocepts, we explai them through a simple LSS system: a epidemic model for the trasmissio of a ifectious disease. A homogeeous populatio is divided ito four classes S, E, I ad T accordig to the health of its idividuals. Let S(t) deote the umber of idividuals who are susceptible to the disease, i.e., who are ot yet ifected at time t. E(t) deotes the umber of members at time t who are exposed but ot yet ifected. I(t) deotes the umber of ifected idividuals, that is,

who are ifectious ad able to spread the disease by cotact with idividuals who are susceptible. T(t) is the umber at time t of treated idividuals. The total populatio size is deoted N = S + E + I + T. The dyamic of the disease ca be described by the followig differetial system: ds SI = bn μs β, N de SI = β ( μ + ε) E, N di = εe ( r + d + μ) I, dt = ri μt, dn = ( b μ ) N di, where the parameter b is the rate for atural birth ad μ that of atural death. The parameter β is the trasmissio rate, d is the rate for disease-related death, ε is the rate at which the exposed idividuals become ifective ad r is the per-capita treatmet rate. The parameters b, μ, β, d ad ε are assumed to be costat. The per-capita treatmet rate may vary with time. (2) Fig 1: This diagram illustrates the dyamical trasfer of the populatio. The system (2) ca be writte by usig the fractios s = S N, e = E N, i = I ad τ = T N of the classes S, E, I ad T i the populatio. These fractios satisfy the system of differetial equatios:

ds = b bs ( β d) si, de βsi ( b ε) e dei, = + + di εe r d b i di = + + + dτ = ri bτ + diτ. ( ) 2, (3) Sice s+ e+ i+ τ = 1, it is sufficiet to cosider the followig system ds = b bs ( β d) si, de = βsi ( b + ε) e + dei, di 2 = εe ( r+ d + b) i+ di. If we suppose that the proportio of ifected idividuals ca be measured ad that we ca cotrol the populatio by choosig the treatmet rate the the above system ca be see as a cotrol system of the form (1) with: the state of the system is x t s t, e t, i t u t = r t ad the measurable output is () = ( () () ()), the cotrol is ( ) ( ) () = i() t. The state space is y t { x 3 : 0 s 1, 0 e 1, 0 i 1, s e i 1} Ω= + +. The aim of Cotrol Theory is to give aswers to the followig questios: 1. Give two cofiguratios x 1 ad x 2 of the system (3), is it possible to steer the x system from 1 to x 2 by choosig a appropriate treatmet strategy r(t)? For a x give iitial state 0 = ( s0, e0, i0) ad a give time T, what are all the possible states that ca be reached i time T by usig all the possible cotrol fuctios r(t)? 2. Does the system (3) possess equilibrium poits? Are they stable or ustable? How ca the cotrol be chose i order to make a give equilibrium (for istace, the disease-free equilibrium) globally asymptotically stable? 3. Sice oly some variables ca be measured (for istace, the proportio i(t) of ifected idividuals i the above epidemic example), is it possible to use them together with the dyamics (3) of the system i order to estimate the o measurable variables s(t) ad e(t) ad how ca this be doe? (4) Questios 1 ivolve the cotrollability of cotrol systems. We shall give i Sectio 2 a geeral defiitio of this otio as well as some simple useful criteria that allow us to

study the cotrollability of some oliear systems. Questios 2 ivoke the stability ad the stabilizatio of oliear system, this will be the subject of Sectio 3 where some classical stability results are exposed. Questios 3 are coected to the observability ad the costructio of observers for dyamical systems. These problems will be addressed i Sectio 4 ad Sectio 5. I each sectio we shall apply the differet tools to various models such as predatorprey systems, fisheries ad bioreactors. Cocerig the epidemic model (2-4), the divers problems metioed above are still uder ivestigatio. We ca give here just some partial aswers. It ca be show that if the per-capita treatmet rate satisfies βε r r0 = d b ε + b the the system has a uique equilibrium poit i the domai Ω. This equilibrium is the disease free equilibrium s = 1, e = i = τ = 0 ad it is globally asymptotically stable, that is, if ( s( 0, ) e( 0, ) i( 0) ) Ω is ay iitial coditio the the solutio of the system (4) startig from this state will coverge to the disease free equilibrium, i.e., s() t 1, e() t 0adi() t 0ast. Therefore, if the treatmet satisfies the coditio (5) the the disease will die out. If r < r0 the the disease free equilibrium becomes ustable ad i this case, there is aother equilibrium which belogs to the iterior of the domai Ω. This equilibrium is the edemic equilibrium ad it is globally asymptotically stable withi the iterior of Ω provided that r < r0. The disease will be edemic i this case. 2. Cotrollability I this sectio, we preset a elemetary overview of a importat property of a system, amely that of cotrollability. This cocept has bee briefly itroduced ad explaied i Basic priciples of mathematical modelig. Here, we deal with the problem of testig cotrollability of systems of the type () = ( (), ()) x t X x t u t. (6) The state vector x(t) belogs to the state space M which will always be here ope coected subset of t u t (5) or a. The cotrol fuctio ( ) are defied o [ 0, ) ad take values i a coected subset U of. These cotrol fuctios are assumed to belog to a admissible cotrol set U ad. This admissible cotrol set is geerally specified by the problem cosidered. Here it is assumed that U ad cotais the set of piecewise costat fuctios as a dese subset, i.e., ay admissible cotrol fuctio ca be approximated by piecewise costat fuctios. The dyamics fuctio X : M x, u. For each cotrol value u U, is assumed to be aalytic o ( ) m

we deote by all x M. u u X the vector field ( vector fuctio ) defied by X ( x) X ( x, u) = for A real-valued fuctio f defied o some ope set O is said to be aalytic if for each x0 O, there exists a positive real umber r > 0 such that f(x) is the sum of a power series for all x with x x0 < r: 0 0. = 0 ( ) ( ) f x = a x x for all x satisfyig x x < r For a aalytic fuctio, the coefficiets a ca be computed as k a = ( ) ( a) f! I a similar way a fuctio f : V of several variables is said to be aalytic if it is locally give by power series. A vector fuctio X : V k (,..., ) X1 x1 xk x = ( x1,..., xk ) X ( x) = X x1 xk is aalytic if all its compoets i (,..., ) X are aalytic. u. For each x M ad each cotrol fuctio u (.) U ad, we deote by Xt x the u(). u solutio of (6) satisfyig X0 ( x) = x. We assume moreover that (). Xt ( x ) is defied for all t [ 0, ). For istace, if X is a liear vector fuctio, that is, X ( x, u) = Ax + Bu, with A ad B beig matrices with appropriate dimesios, the u(). ta t ( ) ( t s) A X x = e x+ e Bu s ds. - - - t 0 ( ) ( ) ( ) TO ACCESS ALL THE 70 PAGES OF THIS CHAPTER, Visit: http://www.eolss.et/eolss-sampleallchapter.aspx

Bibliography Bacciotti, A. (1992). Local stabilizability of oliear cotrol systems, 202 p. Series o Advaces i Mathematics for Applied Scieces. 8. Sigapore, World Scietific. [A well writte ad easy-to-read itroductio to the stabilizatio of cotrol systems]. Basti G. ad Dochai D. (1990). O-lie Estimatio ad Adaptive Cotrol of Bioreactors, 379 p. Elsevier Publisher. [Bioreactor models have specific structures which are thoroughly cosidered i this book] Berard O., Sallet G. ad Sciadra A. (1998). Noliear Observer for a Class of Biological Systems: Applicatio to Validatio of a Phytoplaktoic Growth Model. IEEE Tras. Automat. Cotrol, 43(8), pp 1056-1067. [This paper presets a applicatio of the observers theory to a biological system i order to estimate some state variables that ca ot be measured]. Borard G., Celle-Couee F. ad Gilles G. (1993). Observabilité et observateurs. Systèmes o liéaires, 1. modélisatio estimatio (ed. A.J. Fossard ad D. Normad-Cyrot), pp 177-221. Masso. [This paper (i Frech) is a readable itroductio to the observability ad observers theory]. Brockett R.W. (1983). Asymptotic stability ad feedback stabilizatio. i Differetial Geometric Cotrol Theory (Ed. Brockett R.W., Millma R.S. ad Sussma H.J.), pp. 181-191. Birkhäuser. [This semial article gives some ecessary coditios for the existece of a stabilizig feedback for a oliear system]. Gauthier J.-P., Hammouri H. ad Othma S. (1992). A simple observer for oliear systems, applicatios to bioreactors. IEEE Tras. Automat. Cotrol, 37(6), pp 875-880. [This article provides a expoetial observer for the chemostat example with aother absorbig respose]. Gauthier, J.-P., ad Kupka, I. (2001). Determiistic Observatio Theory ad Applicatios, Cambridge Uiversity Press. [This self cotaied book presets the theoretical foudatios of observability. It requires a good backgroud i differetial geometry]. Jurdjevic, V. (1997). Geometric cotrol theory, 492 p. Cambridge Studies i Advaced Mathematics. 52. Cambridge Uiversity Press. [This book presets the essetial aspects of oliear cotrol theory icludig optimal cotrol]. Nijmeijer H. ad Va der Schaft A.J. (1990). Noliear Dyamical Cotrol Systems, 467 pp. New York, Spriger-Verlag. [This book presets may techiques for the study of iput-output oliear systems]. Touzeau S. ad Gouzé J.-L. (1998). O the stock-recruitmet relatioships i fish populatio models. Evirometal modelig ad assessmet, 3:87-93. [This article gives a stage-structured model of a fish populatio]. Biographical Sketch Abderrahma Iggidr, was bor i Rabat (Morocco). He received the D.E.A i 1989 i applied mathematics ad the Ph.D. degree i 1992, i mathematics ad automatic cotrol, from the Uiversity of Metz (Frace). Sice 1993, he has bee researcher at INRIA (The Frech Natioal Istitute for Research i Computer Sciece ad Cotrol). His research iterests are the field of oliear systems, cotrollability, observability, stabilizatio, observers, ad stabilizatio by estimated state feedback from observers. His research iterests also iclude the applicatio of cotrol theory to the modelig ad the regulatio of biological systems.