ON IDENTIFICATION OF DYNAMICAL SYSTEM PARAMETERS FROM EXPERIMENTAL DATA
|
|
- Brenda West
- 5 years ago
- Views:
Transcription
1 ON IDENTIFICATION OF DYNAMICAL SYSTEM PARAMETERS FROM EXPERIMENTAL DATA M Shatalov*, I Fedotov** * CSIR Manufacturing and Materials PO Box 305, Pretoria 0001, SCIR, South Africa and Department of Mathematics and Statistics, PBX680, Pretoria 0001 FIN Tshwane University of Technology, South Africa, mshatlov@csircoza **Department of Mathematics and Statistics, PBX680, Pretoria 0001 FIN Tshwane University of Technology, South Africa, fedotovi@tutacza Key words: Dynamical system, quadrature rules, least squares method, constraints Abstract In this paper we present a new method of numerical evaluation of unknown coefficients of a dynamical system having available information about unknown phase trajectories at some time values The method consists in the direct integration of given dynamical system with posterior application of a quadrature rules Using the least square method and possible Constraints we obtain a linear system for determining an unknown coefficient A numerical example illustrates the method INTRODUCTION The problem of determination of dynamical systems coefficients from experimental data is solved for a practically important case of the linear systems with respect to their unknown coefficients The problem often occurs in interpretation of experimental data in mathematical biology, ecology, medicine, chemical kinetics, economy etc One of the first example of such problems (the predator-prey equations) consists of a pair of first order, non-linear differential equations frequently used to describe the dynamics of biological systems in which two species interact They were proposed independently by Alfred J Lotka [1] and Vito Volterra in 1926 [2] This system can be written in the form x 1(t) = x 1 (a 11 a 12 x 2 ) x 2(t) = x 2 (ηa 12 x 1 a 22 ) When solved for x 1 and x 2 the above system of equations yields and x 1 = 0, x 1 = 0 1
2 x 1 = a 22, x 1 = a 11 ηa 12 a 12 hence there are two equilibria The solution in the neighborhood of the first (saddle) fixed point does not have any essential physical meaning The second (center) fixed point represents neighbourhood at which both populations uphold their current, non-zero numbers The level of population at which this equilibrium is achieved depends on the chosen values of the parameters a 11, a 12, a 22, η The value of these parameters are, generally speaking, unknown and determining of its values is confronted with serious difficulties In general case the structure of interactions between species in the dynamical system proceeds from a real physical problem and is supposed to be known The method of finding the unknown parameters is based on integration of both parts of equations of the dynamical system and on applying regression methods to the obtained overdetermined system of linear algebraic equations with constraints The least squares method is used for solution of this problem An example of a classical model predator- prey is considered The proposed method could be used for approximate determination of the dynamical system parameters The method can be easily generalized for systems with non-linear dependence on coefficients, for example for systems [3] which describe a deterministic mathematical model for the transmission dynamics of HIV infection in the presence of a preventive vaccine EVALUATION OF COEFFICIENTS Let x (t) =(x 1 (t), x 2 (t),, x m (t)) be a vector in R m depending on variable t where x k (t) : [t 0, t N ] R and f(t, x(t)) = (f 1 (t, x), f 2 (t, x),, f n (t, x)) a vector in R n where f j (t, x) : [t 0, t N ] R m R Consider a dynamical system x k(t) = or, in the matrix form with initial conditions a k,j f j (t, x (t)), (k = 1, 2,, m) (1) x (t) = Af(t, x(t)) (1 ) x (t 0 ) = x 0 = (x 1,0, x 20,, x m, 0 ) (2) The entries a jk of matrix A are unknown but we suppose that there is information about values of functions x k (t) at some points t j = t 0 + jh j It may be, for example, some statistical data of the form t 0 t 1 t i t M x k (t 0 ) = x k,0 x k (t 1 ) = x k,1 x k (t i ) = x k,i x k (t N ) = x k,m, 2
3 where (k = 1, 2,, m) For the sake of simplicity we suppose that h j = h is a constant We suppose that system (1) (3) has on interval [t 0, t N ] a unique solution Moreover, we suppose that the mapping F j : [t 0, t N ] R where F j (t) = f j (t, x(t)) is bounded and has a bounded second derivative with respect to t on the interval [t 0, t N ], that is for all j there exist constants C such that F j (t) C and Fj (t) C for all t [t 0, t N ] To avoid further use of triple subscripts we introduce the following notations using the columns of the transpose a 1,1 a 2,1 a m,1 A = a 1,2 a 2,2 a m,2 = ( ) α 1 α 2 α n a 1,n a 2,n am,n of matrix A, that is we introduce vectors α k = ( a k,1 a k,2 a k,n ) def = ( αn(k 1)+1 α n(k 1)+2 α kn ), (k = 1, 2,, m) and also we denote by α the N 1 column α = ( α 1 α 2 α m ) where N = mn Now system (1) can be rewritten in the form x k(t) = (α k, f(t, x(t)) = α (k 1)n+j f j (t, x), ( k = 1, 2,, m) (3) Integration with respect to t from t 0 to t i, (i = 1, 2,, M) gives t i t 0 x k(t)dt = α (k 1)n+j t i t 0 f j (t, x)dt (4) The left hand side of (5) after integration can be written as t i t 0 x k(t)dt = x k (t i ) x k (t 0 ) = x k,i x k,0 def = k,i (5) We evaluate the integrals in the right hand side of (4) using a quadrature rule, for example the adaptive trapezoidal rule in the form We denote here t i t 0 f j (t, x (t))dt = T i (f j ) def = q i,j (6) q i,j def = I i (f j ) + R i (f j ) def = I i,j + R ij (7) 3
4 where I i,j = h 2 i [f j (t, x (t l )) + f j (t, x (t l 1 ))] (8) and the error functional R i (f j ) can be evaluated by R ij 1 12 C(t N t 0 )h 2 (9) Thus, we obtain the following, generally speaking, overdetermined system for finding the unknown numbers α 1, α N : α (k 1)n+j q ij ki = 0, k = 1, 2,, m; i = 12,, N (10) In many cases the numbers α 1, α N satisfy certain Constraints Frequently these Constraints are linear We limit ourselves to this case Thus system (10) must be resolved providing that the following Constraints are satisfied: N c i,j α j = b i, (11) and system (11) in matrix form can be written as Cα = b (12) where C is an N 1 N matrix and N 1 mn: c 1,1 c 1,2 c 1,N c 2,1 c 2,2 c 2,N C = c N1,1 c N1,2 c N1,N The values of entries of A can be found using method of least squares by minimization of the functional L = 1 2 m α (k 1)n+j q i,j k,i k=1 2 N 1 To minimize functional F we write all M equations L L = 0,, = 0 α 1 α M λ i b i N c i,j α j (13) 4
5 This is equivalent to the system ( N ) N 1 α (k 1)n+j q i,j q i,l + λ i c i,(k 1)n+l = k,i q i,l, (l = 1, 2,, n), (k = 1,, m) (14) and we have to solve this system provided that (11) is satisfied Introducing the notations we obtain the system q i,j q i,l = T j,l ; k,i q i,l = p k,l (15) M 1 T j,l α (k 1)n+j + λ i c i,(k 1)n+l = p k,l ( l = 1, 2,, n), (k = 1,, m) (16) This system together with Constraints can be written in the following matrix form: where T C 1 0 T 0 0 C T C m C 1 C 2 C m 0 α 1 α 1 α m λ p 1 p 2 = p m b (17) T 1,1 T 2,1 T n,1 T 1,2 T 2,2 T n,2 T = ; (18) T 1,n T 2,n T n,n is an n n matrix The matrix C 1 is the n N 1 matrix consisting of first n rows of C, analogously, C 2 consists of next n rows of C etc so that C = ( C 1 C 2 C m ) λ = ( ) λ 1 λ 2 λ N1 ; p1 = ( ) p 11 p 12 p 1n ; ; p m = ( ) ( ) p m,1 p m,2 p m,n ; b = b1 b 2 b N1 Relations (16) are exact The idea of our method is to replace of T j,l with I i,j = h i [f j (t, x (t l )) + f j (t, x (t l 1 ))] It is possible to prove that the error 2 5
6 after this replacement will be of order h meaning that it will be small if h is small Remark 1 If the assumption of boundedness of F is replaced by assumption that F is Lipschitzian or is continuous with bounded variation then inequality (9) must be replaced by R ij C(t N t 0 )h (see for example [4]) and in this case we can not guarantee that the error will be small In what follows, we suppose that q i,j defined by (7) is replaced with r i,j = I i (f j ) (defined by (8)) and p k,l is replaced by s k,l = have to replace T j,l by I j,l = following form: where is an k,i r i,l At the same time we N r i,j r i,l Thus our system can be written in the I C 1 0 I 0 0 C I C m C 1 C 2 C m 0 n n matrix and α 1 α 1 α m λ s 1 s 2 = s m b (19) I 1,1 I 2,1 I n,1 I 1,2 I 2,2 I n,2 I = ; (20) I 1,n I 2,n I n,n s 1 = ( s 11 s 12 s 1n ) ; ; sm = ( s m,1 s m,2 s m,n ) ; EXAMPLE (Predator-Prey Model) x 1(t) = a 11 x 1 (t) a 12 x 1 (t)x 2 (t) + 0x 2 x 2(t) = 0x 1 + ηa 12 x 1 (t)x 2 (t) a 22 x 2 (t) Here m = 2, n = 3 and according to our notations ( ) x1 x = ; f = x 1 ( ) x x 1 x 2 a11 a ; A = 12 0 ; 2 0 ηa 12 a 22 x 2 A = a 11 0 a 12 ηa 12 0 a 22 ; } (21) 6
7 α 1 = a 11 a 12 0 = α 1 α 2 α 3 ; α 2 = 0 ηa 12 a 22 Constraints (Cα = b) can be written as = ηα 2 + α 5 = 0, α 3 = 0, α 4 = 0 In matrix form 0 η ( ) α1 = α Thus, C = η and to find α we must solve the system I 0 C 1 0 I C 2 α 1 α 2 C 1 C 2 0 λ where s 1 α 4 α 5 α 6 α 5 + ηα 2 α 4 α 3 ; α = ( α1 α 2 ) = 0, C 1 = η 0 0, C 2 = (x 1,i x 1,0 ) I i (f 1 ) (x 1,i x 1,0 ) I i (f 2 ) ; s 2 = (x 1,i x 1,0 ) I i (f 3 ) I = Ii 2(f 1) I i (f 1 )I i (f 2 ) I i (f 1 )I i (f 3 ) = I i (f 2 )I i (f 1 ) Ii 2(f 2) I i (f 2 )I i (f 3 ) s 1 s 2 0 (x 2,i x 2,0 ) I i (f 1 ) (x 2,i x 2,0 ) I i (f 2 ) (x 2,i x 2,0 ) I i (f 3 ) I i (f 3 )I i (f 1 ) I i (f 3 )I i (f 2 ) Ii 2(f 2) 7
8 and I i (f 1 ) = h 2 I i (f 3 ) and = h 2 ix (x 1;l + x 1;l 1 ) ; I i (f 2 ) = h 2 ix (x 21;l + x 2;l 1 ) ix (x 1;l x 2;l + x 1;l 1 x 2;l 1 ) ; I i (f 1 ) = h 2 i (x 1,l + x 1,l 1 ), I i (f 2 ) = h 2 i (x 1,l x 2,l + x 1,l 1 x 2,l 1 ), Numerical realizationwas accomplished using MathCad13 i (x 21,l + x 2,l 1 ) Let us generate 120 I i (f 3 ) = h points of solution of the system (21) with the coe cients a 11 = 1:8; a 12 = 0:25; 2 a 22 = 0:7; = 0:1 (a 21 = 0:025) by means of direct solution of the initial value problem with x (0) = 100; y (0) = 10 by the Runge- Kutta method These points are considered further as "experimental data" Application of the proposed method gives us the following set of estimated parameters: ea 11 = 1:839; ea 12 = 0:256; ea 21 = 0:026; ea 22 = 0:715: We see that the values of coe cients are close enough to original values Further, the initial value problem was solved with the new coe cients and old initial conditions Comparison of the "experimental Data" with the estimated solutions are presented in gures 1 and 2 One can see that the estimated values of predators and preys are close enough to the original "experimental data" Divergence of the graphs is stipulated by the errors of the trapezoidal rule used at the stage of numerical integration Application of the quadrature formulae of higher orders allows to achieve more accurate results Numerical realization was accomplished using MathCad13 Let us generate 120 points of solution of the system (21) with the coefficients a 11 = 18, a 12 = 025, a 22 = 07, η = 01 (a 21 = 0025) by means of direct solution of the initial value problem with x (0) = 100, y (0) = 10 by the Runge- Kutta method These points are considered further as experimental data Application of the proposed method gives us the following set of estimated parameters: ã 11 = 1839, ã 12 = 0256, ã 21 = 0026, ã 22 = 0715 We see that the values of coefficients are close enough to original values Further, the initial value problem was solved with the new coefficients and old initial conditions Comparison of the experimental Data with the estimated solutions are presented in figures 1 and 2 One can see that the estimated values of predators and preys are close enough to the original experimental data Divergence of the graphs is stipulated by the errors of the trapezoidal rule used at the stage of numerical integration Application of the quadrature formulae of higher orders allows to achieve more accurate results 8
9 CONCLUSION 1 The Method of identification of unknown parameters of a given dynam- CONCLUSION ical system from experimental data is formulated on the basis of least squares 1 The method Method of identi cation of unknown parameters of a given dynamical system from 2 The experimental error functionaldata is estimated is formulated as a value of onorder thehbasis - step of least in squares statistical data method 3 Example of predator-prey model identification is considered as an exam- of general functional theory andis it estimated is shown that as obtained a value parameters of order give accurate h - the step in 2 Thepleerror statistical data interpretation of the experimental data 3 Example REFERENCES of predator-prey model identi cation is considered as an example of general [1] theory Lotka, AJ, and Elements it is of shown Mathematical that obtained Biology, 1956, parameters Dover, New York give accurate interpretation[2] of Volterra the experimental V, Variations anddata fluctuations of the number of individuals in animal species living together In Animal Ecology McGraw-Hill, 1931 Translated from 1928 edition by R N Chapman REFERENCES [3]Gumel AB, Moghadas SM, Mickens RE Effect of a preventive vaccine on the dynamics of HIV transmission, Communications in Nonlinear Science [1] Lotka, AJ, Elements of Mathematical Biology, 1956, Dover, New York and Numerical Simulations, 9 (2004), pp [2] Volterra [4] SS V, Variations Dragomir and I and Fedotov, uctuations A Grüss type of inequality the number for mapping of individuals of in animal species bounded living variation together and application In Animal to numerical Ecology Analysis, McGraw-Hill, Nonlinear Functional 1931 Translated from 1928 edition by R N Chapman Analysis and applications, Vol6,No3, 2001, pp [3]Gumel AB, Moghadas SM, Mickens RE E ect of a preventive vaccine on the dynamics of HIV transmission, Communications in Nonlinear Science and Numerical Simulations, 9 (2004), pp [4] SS Dragomir and I Fedotov, A Grüss type inequality for mapping of bounded variation and application to numerical 9 Analysis, Nonlinear Functional Analysis and applications, Vol6,No3, 2001, pp
Parameter Identification for the Gumel-Mickens HIV Transmission Model under Missing Observable Data
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 11 (2017), pp. 7759 7770 Research India Publications http://www.ripublication.com/gjpam.htm Parameter Identification for
More informationOBSERVABILITY AND OBSERVERS IN A FOOD WEB
OBSERVABILITY AND OBSERVERS IN A FOOD WEB I. LÓPEZ, M. GÁMEZ, AND S. MOLNÁR Abstract. The problem of the possibility to recover the time-dependent state of a whole population system out of the observation
More informationA NOVEL METHOD OF INTERPOLATION AND EXTRAPOLATION OF FUNCTIONS BY A LINEAR INITIAL VALUE PROBLEM
A OVEL METHOD OF ITERPOLATIO AD EXTRAPOLATIO OF FUCTIOS BY A LIEAR IITIAL VALUE PROBLEM Michael Shatalov Sensor Science and Technolog o CSIR Manuacturing and Materials, P.O.Bo 395, Pretoria, CSIR and Department
More informationSufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations
/. Inst. Maths Applies (1977) 19, 377-383 Sufficiency Conditions for Finite Escape Times in Systems of Quadratic Differential Equations W. M. GETZ AND D. H. JACOBSON National Research Institute for Mathematical
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationANALYTIC SOLUTIONS TO A CLASS OF TWO-DIMENSIONAL LOTKA-VOLTERRA DYNAMICAL SYSTEMS
PROCEEDINGS OF THE INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 18 21, 2000, Atlanta, USA pp. 137 142 ANALYTIC SOLUTIONS TO A CLASS OF TWO-DIMENSIONAL LOTKA-VOLTERRA DYNAMICAL
More informationWeek 15-16: Combinatorial Design
Week 15-16: Combinatorial Design May 8, 2017 A combinatorial design, or simply a design, is an arrangement of the objects of a set into subsets satisfying certain prescribed properties. The area of combinatorial
More informationMatrix Multiplication
3.2 Matrix Algebra Matrix Multiplication Example Foxboro Stadium has three main concession stands, located behind the south, north and west stands. The top-selling items are peanuts, hot dogs and soda.
More informationCartesian Tensors. e 2. e 1. General vector (formal definition to follow) denoted by components
Cartesian Tensors Reference: Jeffreys Cartesian Tensors 1 Coordinates and Vectors z x 3 e 3 y x 2 e 2 e 1 x x 1 Coordinates x i, i 123,, Unit vectors: e i, i 123,, General vector (formal definition to
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationMatrices: 2.1 Operations with Matrices
Goals In this chapter and section we study matrix operations: Define matrix addition Define multiplication of matrix by a scalar, to be called scalar multiplication. Define multiplication of two matrices,
More informationBasic Concepts in Linear Algebra
Basic Concepts in Linear Algebra Grady B Wright Department of Mathematics Boise State University February 2, 2015 Grady B Wright Linear Algebra Basics February 2, 2015 1 / 39 Numerical Linear Algebra Linear
More informationLinear Systems and Matrices
Department of Mathematics The Chinese University of Hong Kong 1 System of m linear equations in n unknowns (linear system) a 11 x 1 + a 12 x 2 + + a 1n x n = b 1 a 21 x 1 + a 22 x 2 + + a 2n x n = b 2.......
More informationReview of Basic Concepts in Linear Algebra
Review of Basic Concepts in Linear Algebra Grady B Wright Department of Mathematics Boise State University September 7, 2017 Math 565 Linear Algebra Review September 7, 2017 1 / 40 Numerical Linear Algebra
More informationThe Application of the Poincart-Transform to the Lotka-Volterra Model
J. Math. Biology 6, 67-73 (1978) Journal of by Springer-Verlag 1978 The Application of the Poincart-Transform to the Lotka-Volterra Model S. B. Hsu Department of Mathematics, University of Utah, Salt Lake
More informationNonstandard Finite Difference Methods For Predator-Prey Models With General Functional Response
Nonstandard Finite Difference Methods For Predator-Prey Models With General Functional Response Dobromir T. Dimitrov Hristo V. Kojouharov Technical Report 2007-0 http://www.uta.edu/math/preprint/ NONSTANDARD
More informationMULTIPLIERS OF THE TERMS IN THE LOWER CENTRAL SERIES OF THE LIE ALGEBRA OF STRICTLY UPPER TRIANGULAR MATRICES. Louis A. Levy
International Electronic Journal of Algebra Volume 1 (01 75-88 MULTIPLIERS OF THE TERMS IN THE LOWER CENTRAL SERIES OF THE LIE ALGEBRA OF STRICTLY UPPER TRIANGULAR MATRICES Louis A. Levy Received: 1 November
More informationElementary Row Operations on Matrices
King Saud University September 17, 018 Table of contents 1 Definition A real matrix is a rectangular array whose entries are real numbers. These numbers are organized on rows and columns. An m n matrix
More informationCHAPTER 6. Direct Methods for Solving Linear Systems
CHAPTER 6 Direct Methods for Solving Linear Systems. Introduction A direct method for approximating the solution of a system of n linear equations in n unknowns is one that gives the exact solution to
More informationA Study of the Variational Iteration Method for Solving. Three Species Food Web Model
Int. Journal of Math. Analysis, Vol. 6, 2012, no. 16, 753-759 A Study of the Variational Iteration Method for Solving Three Species Food Web Model D. Venu Gopala Rao Home: Plot No.159, Sector-12, M.V.P.Colony,
More informationMATH 250 Homework 4: Due May 4, 2017
Due May 4, 17 Answer the following questions to the best of your ability. Solutions should be typed. Any plots or graphs should be included with the question (please include the questions in your solutions).
More informationComputational Aspects of Aggregation in Biological Systems
Computational Aspects of Aggregation in Biological Systems Vladik Kreinovich and Max Shpak University of Texas at El Paso, El Paso, TX 79968, USA vladik@utep.edu, mshpak@utep.edu Summary. Many biologically
More information5 Alfred Lotka, Vito Volterra, and Population Cycles
5 Alfred Lotka, Vito Volterra, and opulation Cycles Dr. Umberto D Ancona entertained me several times with statistics that he was compiling about fishing during the period of the war and in periods previous
More informationDirect Methods for Solving Linear Systems. Matrix Factorization
Direct Methods for Solving Linear Systems Matrix Factorization Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011
More informationDYNAMICS OF THE ZEROS OF FIBONACCI POLYNOMIALS M. X. He Nova Southeastern University, Fort Lauderdale, FL. D, Simon
M. X. He Nova Southeastern University, Fort Lauderdale, FL D, Simon Nova Southeastern University, Fort Lauderdale, FL P. E. Ricci Universita degli Studi di Roma "La Sapienza," Rome, Italy (Submitted December
More informationApplied Differential Equation. November 30, 2012
Applied Differential Equation November 3, Contents 5 System of First Order Linear Equations 5 Introduction and Review of matrices 5 Systems of Linear Algebraic Equations, Linear Independence, Eigenvalues,
More informationMATRICES. a m,1 a m,n A =
MATRICES Matrices are rectangular arrays of real or complex numbers With them, we define arithmetic operations that are generalizations of those for real and complex numbers The general form a matrix of
More informationMatrix Arithmetic. j=1
An m n matrix is an array A = Matrix Arithmetic a 11 a 12 a 1n a 21 a 22 a 2n a m1 a m2 a mn of real numbers a ij An m n matrix has m rows and n columns a ij is the entry in the i-th row and j-th column
More informationDEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix to upper-triangular
form) Given: matrix C = (c i,j ) n,m i,j=1 ODE and num math: Linear algebra (N) [lectures] c phabala 2016 DEN: Linear algebra numerical view (GEM: Gauss elimination method for reducing a full rank matrix
More informationTHE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A 2 2 SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS
THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS Maria P. Skhosana and Stephan V. Joubert, Tshwane University
More information10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )
c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side
More informationTheory of Ordinary Differential Equations. Stability and Bifurcation I. John A. Burns
Theory of Ordinary Differential Equations Stability and Bifurcation I John A. Burns Center for Optimal Design And Control Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute
More informationMath 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
More informationDirect Methods for Solving Linear Systems. Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le
Direct Methods for Solving Linear Systems Simon Fraser University Surrey Campus MACM 316 Spring 2005 Instructor: Ha Le 1 Overview General Linear Systems Gaussian Elimination Triangular Systems The LU Factorization
More informationPersistence and global stability in discrete models of Lotka Volterra type
J. Math. Anal. Appl. 330 2007 24 33 www.elsevier.com/locate/jmaa Persistence global stability in discrete models of Lotka Volterra type Yoshiaki Muroya 1 Department of Mathematical Sciences, Waseda University,
More informationMath 1280 Notes 4 Last section revised, 1/31, 9:30 pm.
1 competing species Math 1280 Notes 4 Last section revised, 1/31, 9:30 pm. This section and the next deal with the subject of population biology. You will already have seen examples of this. Most calculus
More informationChaos and adaptive control in two prey, one predator system with nonlinear feedback
Chaos and adaptive control in two prey, one predator system with nonlinear feedback Awad El-Gohary, a, and A.S. Al-Ruzaiza a a Department of Statistics and O.R., College of Science, King Saud University,
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2017/18 Part 2: Direct Methods PD Dr.
More informationMath 216 Final Exam 14 December, 2012
Math 216 Final Exam 14 December, 2012 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationLagrangian data assimilation for point vortex systems
JOT J OURNAL OF TURBULENCE http://jot.iop.org/ Lagrangian data assimilation for point vortex systems Kayo Ide 1, Leonid Kuznetsov 2 and Christopher KRTJones 2 1 Department of Atmospheric Sciences and Institute
More information6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities
6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities 1 Outline Outline Dynamical systems. Linear and Non-linear. Convergence. Linear algebra and Lyapunov functions. Markov
More informationStage-structured Populations
Department of Biology New Mexico State University Las Cruces, New Mexico 88003 brook@nmsu.edu Fall 2009 Age-Structured Populations All individuals are not equivalent to each other Rates of survivorship
More information1 Determinants. 1.1 Determinant
1 Determinants [SB], Chapter 9, p.188-196. [SB], Chapter 26, p.719-739. Bellow w ll study the central question: which additional conditions must satisfy a quadratic matrix A to be invertible, that is to
More informationA Moment Closure Method for Stochastic Chemical Reaction Networks with General Kinetics
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 70 2013 785-800 ISSN 0340-6253 A Moment Closure Method for Stochastic Chemical Reaction Networks with General
More informationNUMERICAL MATHEMATICS & COMPUTING 7th Edition
NUMERICAL MATHEMATICS & COMPUTING 7th Edition Ward Cheney/David Kincaid c UT Austin Engage Learning: Thomson-Brooks/Cole wwwengagecom wwwmautexasedu/cna/nmc6 October 16, 2011 Ward Cheney/David Kincaid
More informationII. Determinant Functions
Supplemental Materials for EE203001 Students II Determinant Functions Chung-Chin Lu Department of Electrical Engineering National Tsing Hua University May 22, 2003 1 Three Axioms for a Determinant Function
More information1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )
Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical
More informationApplications in Biology
11 Applications in Biology In this chapter we make use of the techniques developed in the previous few chapters to examine some nonlinear systems that have been used as mathematical models for a variety
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
More informationIntroduction to Dynamical Systems
Introduction to Dynamical Systems Autonomous Planar Systems Vector form of a Dynamical System Trajectories Trajectories Don t Cross Equilibria Population Biology Rabbit-Fox System Trout System Trout System
More informationMatrix Differentiation
Matrix Differentiation CS5240 Theoretical Foundations in Multimedia Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (NUS) Matrix Differentiation
More informationarxiv: v1 [math.ds] 11 Feb 2011
The role of space in the exploitation of resources Y. Kang and N. Lanchier arxiv:0.83v [math.ds] Feb 0. Introduction Abstract In order to understand the role of space in ecological communities where each
More informationEquilibrium points, stability and numerical solutions of fractional-order predator prey and rabies models
J. Math. Anal. Appl. 325 (2007) 542 553 www.elsevier.com/locate/jmaa Equilibrium points, stability and numerical solutions of fractional-order predator prey and rabies models E. Ahmed a, A.M.A. El-Sayed
More informationMatrix Algebra & Elementary Matrices
Matrix lgebra & Elementary Matrices To add two matrices, they must have identical dimensions. To multiply them the number of columns of the first must equal the number of rows of the second. The laws below
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationExistence of Positive Periodic Solutions of Mutualism Systems with Several Delays 1
Advances in Dynamical Systems and Applications. ISSN 973-5321 Volume 1 Number 2 (26), pp. 29 217 c Research India Publications http://www.ripublication.com/adsa.htm Existence of Positive Periodic Solutions
More informationChapter 1 Matrices and Systems of Equations
Chapter 1 Matrices and Systems of Equations System of Linear Equations 1. A linear equation in n unknowns is an equation of the form n i=1 a i x i = b where a 1,..., a n, b R and x 1,..., x n are variables.
More informationVector analysis and vector identities by means of cartesian tensors
Vector analysis and vector identities by means of cartesian tensors Kenneth H. Carpenter August 29, 2001 1 The cartesian tensor concept 1.1 Introduction The cartesian tensor approach to vector analysis
More informationIDENTIFICATION OF MULTI-DEGREE-OF-FREEDOM NON-LINEAR SYSTEMS UNDER RANDOM EXCITATIONS BY THE REVERSE PATH SPECTRAL METHOD
Journal of Sound and Vibration (1998) 213(4), 673 78 IDENTIFICATION OF MULTI-DEGREE-OF-FREEDOM NON-LINEAR SYSTEMS UNDER RANDOM EXCITATIONS BY THE REVERSE PATH SPECTRAL METHOD Acoustics and Dynamics Laboratory,
More informationMULTIPOINT BOUNDARY VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS. Katalin Károlyi Department of Applied Analysis, Eötvös Loránd University
MULTIPOINT BOUNDARY VALUE PROBLEMS FOR ORDINARY DIFFERENTIAL EQUATIONS Katalin Károlyi Department of Applied Analysis, Eötvös Loránd University HU-1117 Budapest, Pázmány Péter sétány 1/c. karolyik@cs.elte.hu
More informationMathematics 13: Lecture 10
Mathematics 13: Lecture 10 Matrices Dan Sloughter Furman University January 25, 2008 Dan Sloughter (Furman University) Mathematics 13: Lecture 10 January 25, 2008 1 / 19 Matrices Recall: A matrix is a
More informationA Discrete Model of Three Species Prey- Predator System
ISSN(Online): 39-8753 ISSN (Print): 347-670 (An ISO 397: 007 Certified Organization) Vol. 4, Issue, January 05 A Discrete Model of Three Species Prey- Predator System A.George Maria Selvam, R.Janagaraj
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationBoolean Inner-Product Spaces and Boolean Matrices
Boolean Inner-Product Spaces and Boolean Matrices Stan Gudder Department of Mathematics, University of Denver, Denver CO 80208 Frédéric Latrémolière Department of Mathematics, University of Denver, Denver
More informationLocal Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)
Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two
More informationDiscrete Math, Spring Solutions to Problems V
Discrete Math, Spring 202 - Solutions to Problems V Suppose we have statements P, P 2, P 3,, one for each natural number In other words, we have the collection or set of statements {P n n N} a Suppose
More informationObserver design for a general class of triangular systems
1st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 014. Observer design for a general class of triangular systems Dimitris Boskos 1 John Tsinias Abstract The paper deals
More informationNotes on Blackwell s Comparison of Experiments Tilman Börgers, June 29, 2009
Notes on Blackwell s Comparison of Experiments Tilman Börgers, June 29, 2009 These notes are based on Chapter 12 of David Blackwell and M. A.Girshick, Theory of Games and Statistical Decisions, John Wiley
More informationADAPTIVE CONTROL AND SYNCHRONIZATION OF A GENERALIZED LOTKA-VOLTERRA SYSTEM
ADAPTIVE CONTROL AND SYNCHRONIZATION OF A GENERALIZED LOTKA-VOLTERRA SYSTEM Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi, Chennai-600
More informationMarkov Chains. Sarah Filippi Department of Statistics TA: Luke Kelly
Markov Chains Sarah Filippi Department of Statistics http://www.stats.ox.ac.uk/~filippi TA: Luke Kelly With grateful acknowledgements to Prof. Yee Whye Teh's slides from 2013 14. Schedule 09:30-10:30 Lecture:
More informationThe Role of Space in the Exploitation of Resources
Bull Math Biol (2012) 74:1 44 DOI 10.1007/s11538-011-9649-1 ORIGINAL ARTICLE The Role of Space in the Exploitation of Resources Y. Kang N. Lanchier Received: 29 September 2010 / Accepted: 25 February 2011
More informationThis operation is - associative A + (B + C) = (A + B) + C; - commutative A + B = B + A; - has a neutral element O + A = A, here O is the null matrix
1 Matrix Algebra Reading [SB] 81-85, pp 153-180 11 Matrix Operations 1 Addition a 11 a 12 a 1n a 21 a 22 a 2n a m1 a m2 a mn + b 11 b 12 b 1n b 21 b 22 b 2n b m1 b m2 b mn a 11 + b 11 a 12 + b 12 a 1n
More informationVectors and matrices: matrices (Version 2) This is a very brief summary of my lecture notes.
Vectors and matrices: matrices (Version 2) This is a very brief summary of my lecture notes Matrices and linear equations A matrix is an m-by-n array of numbers A = a 11 a 12 a 13 a 1n a 21 a 22 a 23 a
More informationOn graph differential equations and its associated matrix differential equations
Malaya Journal of Matematik 1(1)(2013) 1 9 On graph differential equations and its associated matrix differential equations J. Vasundhara Devi, a, R.V.G. Ravi Kumar b and N. Giribabu c a,b,c GVP-Prof.V.Lakshmikantham
More informationRomberg Integration. MATH 375 Numerical Analysis. J. Robert Buchanan. Spring Department of Mathematics
Romberg Integration MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Spring 019 Objectives and Background In this lesson we will learn to obtain high accuracy approximations to
More informationNONSTANDARD NUMERICAL METHODS FOR A CLASS OF PREDATOR-PREY MODELS WITH PREDATOR INTERFERENCE
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 67 75. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationOutline. Input to state Stability. Nonlinear Realization. Recall: _ Space. _ Space: Space of all piecewise continuous functions
Outline Input to state Stability Motivation for Input to State Stability (ISS) ISS Lyapunov function. Stability theorems. M. Sami Fadali Professor EBME University of Nevada, Reno 1 2 Recall: _ Space _
More informationLecture 20 : Markov Chains
CSCI 3560 Probability and Computing Instructor: Bogdan Chlebus Lecture 0 : Markov Chains We consider stochastic processes. A process represents a system that evolves through incremental changes called
More informationChoose three of: Choose three of: Choose three of:
MATH Final Exam (Version ) Solutions July 8, 8 S. F. Ellermeyer Name Instructions. Remember to include all important details of your work. You will not get full credit (or perhaps even any partial credit)
More informationChapter 4. Systems of ODEs. Phase Plane. Qualitative Methods
Chapter 4 Systems of ODEs. Phase Plane. Qualitative Methods Contents 4.0 Basics of Matrices and Vectors 4.1 Systems of ODEs as Models 4.2 Basic Theory of Systems of ODEs 4.3 Constant-Coefficient Systems.
More informationLotka-Volterra Models Nizar Ezroura M53
Lotka-Volterra Models Nizar Ezroura M53 The Lotka-Volterra equations are a pair of coupled first-order ODEs that are used to describe the evolution of two elements under some mutual interaction pattern.
More informationSolve Linear System with Sylvester s Condensation
International Journal of Algebra, Vol 5, 2011, no 20, 993-1003 Solve Linear System with Sylvester s Condensation Abdelmalek Salem Department of Mathematics University of Tebessa, 12002 Algeria And Faculty
More informationGeneralized Riccati Equations Arising in Stochastic Games
Generalized Riccati Equations Arising in Stochastic Games Michael McAsey a a Department of Mathematics Bradley University Peoria IL 61625 USA mcasey@bradley.edu Libin Mou b b Department of Mathematics
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationSTUDY OF THE DYNAMICAL MODEL OF HIV
STUDY OF THE DYNAMICAL MODEL OF HIV M.A. Lapshova, E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper is devoted to the study of the dynamical model of HIV. An application
More informationOn pairwise comparison matrices that can be made consistent by the modification of a few elements
Noname manuscript No. (will be inserted by the editor) On pairwise comparison matrices that can be made consistent by the modification of a few elements Sándor Bozóki 1,2 János Fülöp 1,3 Attila Poesz 2
More informationInput-State Incidence Matrix of Boolean Control Networks and Its Applications
Input-State Incidence Matrix of Boolean Control Networks and Its Applications Yin Zhao, Hongsheng Qi, Daizhan Cheng The Key Laboratory of Systems & Control Institute of Systems Science, Academy of Mathematics
More informationMath 215/255 Final Exam, December 2013
Math 215/255 Final Exam, December 2013 Last Name: Student Number: First Name: Signature: Instructions. The exam lasts 2.5 hours. No calculators or electronic devices of any kind are permitted. A formula
More informationBIOS 5970: Plant-Herbivore Interactions Dr. Stephen Malcolm, Department of Biological Sciences
BIOS 5970: Plant-Herbivore Interactions Dr. Stephen Malcolm, Department of Biological Sciences D. POPULATION & COMMUNITY DYNAMICS Week 10. Population models 1: Lecture summary: Distribution and abundance
More informationDESIGNING CNN GENES. Received January 23, 2003; Revised April 2, 2003
Tutorials and Reviews International Journal of Bifurcation and Chaos, Vol. 13, No. 10 (2003 2739 2824 c World Scientific Publishing Company DESIGNING CNN GENES MAKOTO ITOH Department of Information and
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 6, Number 1, Winter 1997 ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM J. CHATTOPADHYAY, E. BERETTA AND F. SOLIMANO ABSTRACT. The present paper deals with
More informationInvertible Matrices over Idempotent Semirings
Chamchuri Journal of Mathematics Volume 1(2009) Number 2, 55 61 http://www.math.sc.chula.ac.th/cjm Invertible Matrices over Idempotent Semirings W. Mora, A. Wasanawichit and Y. Kemprasit Received 28 Sep
More informationIntroduction to Quantitative Techniques for MSc Programmes SCHOOL OF ECONOMICS, MATHEMATICS AND STATISTICS MALET STREET LONDON WC1E 7HX
Introduction to Quantitative Techniques for MSc Programmes SCHOOL OF ECONOMICS, MATHEMATICS AND STATISTICS MALET STREET LONDON WC1E 7HX September 2007 MSc Sep Intro QT 1 Who are these course for? The September
More informationPhysics 116A Determinants
Physics 116A Determinants Peter Young (Dated: February 5, 2014) I. DEFINITION The determinant is an important property of a square n n matrix. Such a matrix, A say, has n 2 elements, a ij, and is written
More informationNotes on Mathematics
Notes on Mathematics - 12 1 Peeyush Chandra, A. K. Lal, V. Raghavendra, G. Santhanam 1 Supported by a grant from MHRD 2 Contents I Linear Algebra 7 1 Matrices 9 1.1 Definition of a Matrix......................................
More informationAlgorithmic Game Theory and Applications. Lecture 4: 2-player zero-sum games, and the Minimax Theorem
Algorithmic Game Theory and Applications Lecture 4: 2-player zero-sum games, and the Minimax Theorem Kousha Etessami 2-person zero-sum games A finite 2-person zero-sum (2p-zs) strategic game Γ, is a strategic
More information/Q~ =/V~ (b, + ~. a~inj), j=l
Bulletin of Mathematical Biology, Vol. 44, No. 4, pp. 579-585, 1982. Printed in Great Britain. 0092-8240[82/0405794)7503.00/0 Pergamon Press Ltd. 9 1982 Society for Mathematical Biology NOTE GRAPH THEORETCAL
More informationStability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations
Stability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations Sreelakshmi Manjunath Department of Electrical Engineering Indian Institute of Technology Madras (IITM), India JTG Summer
More informationLinear Algebra (Review) Volker Tresp 2018
Linear Algebra (Review) Volker Tresp 2018 1 Vectors k, M, N are scalars A one-dimensional array c is a column vector. Thus in two dimensions, ( ) c1 c = c 2 c i is the i-th component of c c T = (c 1, c
More information