MINISTRY OF EDUCATION AND TRAINING QUY NHON UNIVERSITY DAO THI HAI YEN ON THE STATE ESTIMATION PROBLEM FOR SOME CLASSES OF DYNAMICAL SYSTEMS AND ITS APPLICATION DOCTORAL THESIS IN MATHEMATICS Binh Din[.]
MINISTRY OF EDUCATION AND TRAINING QUY NHON UNIVERSITY DAO THI HAI YEN ON THE STATE ESTIMATION PROBLEM FOR SOME CLASSES OF DYNAMICAL SYSTEMS AND ITS APPLICATION DOCTORAL THESIS IN MATHEMATICS Binh Dinh - 2023 MINISTRY OF EDUCATION AND TRAINING QUY NHON UNIVERSITY DAO THI HAI YEN ON THE STATE ESTIMATION PROBLEM FOR SOME CLASSES OF DYNAMICAL SYSTEMS AND ITS APPLICATION Speciality: Mathematical Analysis Code: 46 01 02 Reviewer 1: Prof Dr Vu Ngoc Phat Reviewer 2: Assoc Prof Dr Nguyen Dinh Phu Reviewer 3: Assoc Prof Dr Pham Quy Muoi Board of Supervisors: Assoc Prof Dr Dinh Cong Huong Binh Dinh - 2023 Declaration This thesis was completed at the Department of Mathematics and Statistics, Quy Nhon University under the guidance of Assoc Prof Dinh Cong Huong I hereby declare that the results presented in here are new and original All of them were published in peer-reviewed journals and conferences For using results from joint papers I have gotten permissions from my co-authors Binh Dinh, 2023 PhD student Dao Thi Hai Yen i Acknowledgments First, I would like to express my gratitude to the two offices Phu Yen University, my workplace, and Quy Nhon University, where I have learnt and studied for a long time I would also like to express my gratitude to all people working at the Department of Postgraduate Training and the Department of Mathematics and Statistics of Quy Nhon University Next, the thesis was carried out during the years I have been a PhD student at the Department of Mathematics and Statistics, Quy Nhon University On the occasion of completing the thesis, I would like to express the deep gratitude to Assoc Prof Dr Dinh Cong Huong not only for his teaching and scientific leadership, but also for the helping me access to the academic environment through the workshops, mini courses that assist me in broadening my thinking to get the entire view on the related issues in my research Finally, I wish to acknowledge my father, my parents in law for supporting me in every decision And, my enormous gratitude goes to my husband and sons for their love and patience during the time I was working intensively to complete my PhD program Finally, my sincere thank goes to my mother for guiding me to math and this thesis is dedicated to her ii Contents Table of Notations List of drawings Introduction Preliminaries 1.1 1.2 Some basic concepts 1.1.1 Stability criteria of some classes of dynamic systems 1.1.2 Additional lemmas 14 State observer design problems 15 1.2.1 Full order state observer 15 1.2.2 Reduced order state observer 17 1.2.3 Linear functional state observer 17 1.2.4 Interval observer 19 A new method for designing observers of a nonlinear time-delay Glucose-Insulin system 21 2.1 A novel state transformation 22 2.2 Application to the GI model 28 2.2.1 State transformation for the GI model 28 2.2.2 State observer design for the GI model 29 2.2.3 Simulation results 30 A new observer for interconnected time-delay systems and its applications to fault detection problem 34 iii 3.1 New state transformation 37 3.2 Fault detection observer 45 3.3 A numerical example 46 Distributed functional interval observer design for large-scale networks impulsive systems 52 4.1 Designing distributed linear functional interval observers 53 4.2 Existence conditions of distributed linear functional interval observers 55 4.3 Solving unknown matrices 57 4.4 Numerical examples 60 iv Table of Notations Rn+ Rn×m Z+ A () A, B ∈ Rn×m , A () B In 1n : The set of all n-dimensional real non-negative vectors : The set of all real n × m dimensional matrices : Z ∩ R+ : Real part of complex number s : Transpose matrix of the matrix A : Matrix inverse of matrix A : The Moore-Penrose-inverse of A : Rank of the matrix A : The determinant of the matrix A : The matrix A is a positively definite symmetric matrix (non-negative) : All the elements of the matrix A are positive (non-negative) : aij > (>) bij , i ∈ {1, , n}, j ∈ {1, , m} : Identity matrix of size n × n : Vector in Rn with all elements equal to one T i.e 1 | {z } n ? Lr∞ M = max(A, B) A+ A− |A| x(t¯i+1 ) n S>0 P K K∞ L KL , Nullity A : Symmetric term in a symmetric matrix : The set of all inputs v(t) with the property ||v|| < ∞ : The matrix where each entry is mij = max(aij , bij ) : max(A, 0) : A+ − A : A+ + A− : The left-sided limit of x(t) for t → ti+1 : The set of n × n positive definite matrices : {γ : R+ → R+ | γ is continuous, γ(0) = 0, γ(r) > 0, r > 0} : {γ ∈ P | γ is strictly increasing} : {γ ∈ K | γ is unbounded} : {γ : R+ → R+ | γ is continuous and strictly decreasing with limt→∞ γ(t) = 0} : {β : R+ × R+ → R+ | β is continuous, β(., t) ∈ K ∀t > 0, β(r, ) ∈ L ∀r > 0} : Equal to by definition : The dimension of the N (A) = {x ∈ Rn |Ax = 0} List of drawings Trang Figure 1.1 Responses of x(t) and its estimation 11 Figure 1.2 Responses of x(t), x− (t) and x+ (t) 18 Figure 2.1 Responses of x2 (t) and x ˆ2 (t − 3) 33 Figure 2.2 Responses of x ˆ2 (t − 3) and x2 (t − 3) 34 Figure 3.1 Residual generator using third-order observer effectively triggers fault in the 1-th subsystem 53 Residual generator using third-order observer effectively triggers fault in the 2-th subsystem 53 Figure 4.1 Responses of x11 (t), x12 (t), x13 (t), x21 (t), x22 (t) and x23 (t) 68 Figure 4.2 Responses of z1 (t) = x13 (t), z1− (t) and z1+ (t) 69 Figure 4.3 Responses of z2 (t) = x23 (t), z2− (t) and z2+ (t) 69 Figure 4.4 Responses of xi1 (t) and xi2 (t) 70 Figure 4.5 Responses of zi (t) = xi1 (t), zi − (t) and zi + (t) 71 Figure 3.2 Introduction Many real-world physical systems such as electricity, water, communication networks, etc are often modeled by systems of differential equations (see [4], [19], [20], [50], [57], [58]) The qualitative theory of these dynamical systems will bring many useful applications In recent decades, along with the development of qualitative theory, the control theory has become one of important research directions To meet the actual requirements, control systems (biology, medicine, electricity, water, communication networks ,) are becoming more and more diverse and complex These control systems are often very large, in which the controls, devices, and subsystems are linked together through transmission lines and connection bands Therefore, the information received about the state of the system is often delayed (i.e no instantaneous state information is received, but only delayed information is received), missing (i.e not received full information of the whole system but only some information is received), lost (i.e information is provided intermittently), interference (information received is incorrect) Consequences inevitably face failures, disturbances, and many other complications The operation of control systems is usually based on the state information of the systems However, in some cases, these outputs often not provide enough state information for the operation of systems Therefore, the problem of designing functional state observes or designing interval observers for functions of state variables to estimate some more unknown states, which are enough to operate the system at a cost of operation, has become known as an important control problem This problem has many applications in other problems such as fault detection problems, fault estimation problems, identification problems, and signal processing problems, In control theory, a state observer is an auxiliary dynamical system that mirrors the behavior of a physical system, and it is driven by input and output measurements of the physical system in order to provide an estimate of the internal states of the physical system The primary consideration in the design of an observer is that the estimate of the states should be close to the actual value of the system states The next requirement is to provide conditions for the existence of the state observers Some direct applications of the state observer problem include: using estimated state information to design controllers for the considered dynamic system; using estimation information to solve problems of detecting, isolating, and estimating actuator faults as well as sensor faults Since the 60s of the twentieth century, along with the development of control problem, the state estimation problem has attracted a lot of research attention due to the increase in its practical applications (see, for example, [60], [12], [16], [22], [13], [79], [38], [42], [43], [57], [10], [48]) In particular, Luenberger (1966) (see [60]) studied the following control system with input u(t) ∈ Rm and output y(t) ∈ Rp : x(t) ˙ = Ax(t) + Bu(t), t ≥ 0, (1) y(t) = Cx(t), (2) where x(t) ∈ Rn is the state vector, A ∈ Rn×n , B ∈ Rn×m and C ∈ Rp×n are known constant matrices He proposed the following full order state observers for this system: x ˆ˙ (t) = Aˆ x(t) + Bu(t) + L(y(t) − yˆ(t)), t ≥ 0, yˆ(t) = Cx ˆ(t), where L ∈ Rn×p is the observer gain matrix, x ˆ(t) ∈ Rn is an estimate of the state x(t) Darouach et al (1994) (see [12]) developed full order state observers reported in [60] to deal with system (1)-(2), where an extra term Dω(t) is added in the equation (1) The authors proposed the following observer: z(t) ˙ = N z(t) + Ly(t) + Gu(t) t ≥ 0, x ˆ(t) = z(t) − Ey(t), where z(t) ∈ Rn is the observer state vector, N , L, G and E are the observer gain matrices and x ˆ(t) ∈ Rn is an estimate of the state x(t) Deng et al (2004) (see [16]) extended the results of Darouach et al (1994) to a class of systems of differential equations with two unknown inputs Eω(t) and Dω(t) appear in equations (1) and (2), respectively (i.e., noise appears in the state and output information of control system) The authors transformed the systems (1) and (2) (in case there are two more unknown inputs) to the following form: x ˆ˙ (t) = Aˆ x(t) + B¯ u(t) + E1 ω1 (t), t ≥ 0, = Cx ˆ(t) y1 (t) Then, they proposed state observers as follows: ˙ ζ(t) = N ζ(t) + Ly1 (t) + (K − M C)¯ u(t) t ≥ 0, zˆ(t) = ζ(t) + M y1 (t), where ζ(t) ∈ Rr , zˆ(t) is the estimated vector of z(t) = Kx(t), N , L and M are gain observer matrices Fairman et al (1986); Darouach et al (1999) (see [22], [13]) considered the following class of delayed linear systems: x(t) ˙ = Ax(t) + Ad x(t − τ ) + Bu(t), t ≥ 0, (3) x(θ) = φ(θ), θ ∈ [−τ, 0], (4) y(t) = (5) Cx(t), where φ(θ) is the initial condition function, τ > is the time delay, x(t) ∈ Rn is the state vector, u(t) ∈ Rm is the control vector, y(t) ∈ Rp is the output vector, matrices A, Ad , B and C are constant and of appropriate dimensions The authors proposed the following reduced order observer without time delays to estimate the (n − p) components of the state vector: x ˆ(t) = Dw(t) + Ey(t), w(t) ˙ = N w(t) + Jy(t) + Jd y(t − τ ) + Hu(t), where x ˆ(t) is the estimated vector of x(t), and D, E, N , J, Jd and H are gain observer matrices Trinh et al (2010) (see [79]) considered the time-delayed term in the equation (3) as an unknown input, i.e, Ad x(t − τ ) = W ω(t), and moved back to the state estimation problem for systems without