non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.
|
|
- Damian Jordan
- 5 years ago
- Views:
Transcription
1 LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis As model equaions in and of hemselves Populaion dynamics moivaes he sudy of linear models whose coefficien marices are non-negaive or posiive. P= ( pij, i, j m P or P> means p or p > for all i, j f( x = Px R R = ij m m maps non -negaive cone The soluion of he marix, difference equaion x ( = Px (, Z x( = x remains in he non-negaive cone, which we say is (forward invarian. ij Models make simplifying assumpions. Homogeneiies Demographic homogeneiy Temporal homogeneiy Spaial homogeneiy We ll consider relaxing he firs wo assumpions. Demographic inhomogeneiy Differences among individuals include: Age Size Gender Geneics Behavioral : A he exremes: Classical models (mos oal populaion size is he sae variable (all individuals reaed as idenical Individual based models Each individual is a sae variable In beween: Srucured models Sae variables are classes of individuals (individuals wihin classes are idenical
2 A Classic Example: (Leslie Age Srucured Marix Models x x i vecor of age classes x i number of individuals in age class i,,...,m For example, age classes could be in decades: x pre-eens, x eens, x 3 wenies, For i =,3,, m xi( = τ i, i xi( where τ ii,- is he fracion ha survives from age i o age i. For he las age class i = m x ( = τ x ( τ x ( m m, m m mm m For he i = age class, we have o consider newborns If fi is he per capia offspring produced by age class i per uni ime, hen m = i = i i x ( f x ( where In marix form : f f f 3 f m L x ( = Lx ( 3 m,m mm A more general marix model allows for ransiions among all classes for newborns o fall ino more han one class. x ( = Px ( P= F T, F = f, T = τ ( ij ( ij f ij per capia j-class newborns from i-class individuals or L f f f m f m 3 m,m mm ij fracion of j-class individuals ha move o i-class m ij, ij i The asympoic dynamics of soluions depend on he eigenvalues of he ransiion marix A. Does non-negaiviy ell us anyhing abou eigenvalues?
3 D Example a b Consider P = >. c d Could be a Leslie age-srucured model: Eigenvalues : L ( ( 4 λ = a d ± a d bc ± Some algebra shows : λ < λ b Eigenvecors : v± = λ± a More algebra shows v > and v is no posiive (second componen is negaive f f Asympoic dynamics of he equaion Similarly, Noe x ( = Px ( x ( = Px x = cv c v x( x c > = Px = cav cav x( = cλ v c λ v and x( = Px( = cλ v c λ v x ( = cλ v cλ v λ λ x ( = cλ v cλ v > lim x ( = (unbounded growh < lim x ( = (exincion Moreover, x ( cλ v c λ v cv c( λ / λ v = = x ( c v c v cv c λ / λ v Therefore, λ λ x ( v lim = x ( v ( Regardless of he populaion growh, he normalized (class disribuion sabilizes! Wha made all his work mahemaically? The posiive x marix has a posiive, simple, sricly dominan eigenvalue wih an associaed posiive eigenvecor no oher independen eigenvecor is non-negaive Perron s Theorem A posiive m x m marix has a posiive, simple, sricly dominan eigenvalue wih an associaed posiive eigenvecor no oher independen eigenvecor is non-negaive 3
4 As in he x case, we have he resul---- Suppose P > is a posiive m x m marix, λ > is is sricly dominan (simple eigenvalue, and v > is an associaed eigenvecor. Then he soluion of saisfies x ( = Px (, x lim x ( = if λ < lim x ( = if λ > x ( v In eiher case : lim = x ( v When does a non - negaive marix P a b = c d have a posiive eigenvalue & eigenvecor and no oher independen eigenvecor non-negaive? No boh b and c can equal. Assume WLOG ha b >. A Drawback: In applicaions P > is rarely posiive! An invesigaion of he formulas λ ( a d ( d a 4 bc = b v = ( d a ( d a 4bc shows a necessary and sufficien condiion is bc >, A permuaion marix M is he ideniy I wih is rows reordered. MP is he marix P wih is rows reordered PM is he marix P wih is columns reordered MPM is he marix P wih is rows and columns reordered (by he same permuaion i.e., boh off diagonal enries are nonzero. Such a marix is called irreducible. 4
5 DEFINITION: P is reducible if for some permuaion marix M A MPM = B C A and B are square marices Oherwise P is irreducible. In he dynamical model, irreducible means here is flow from each class o all oher classes (in a finie number of seps. Frobenius Theorem A non-negaive, irreducible m x m marix has a posiive, simple, dominan eigenvalue: λ > λ i wih an associaed posiive eigenvecor no oher independen eigenvecor is non-negaive F. R. Ganmacher, The Theory of Marices, Volume II, Chelsea, 96 Berman & Plemmons, Nonnegaive Marices in he Mahemaical Sciences, SIAM, 994 Wha is missing? Sric dominance of he posiive eigenvalue : λ > λ ι This is necessary for convergence of he normalized disribuion : x ( v lim = x ( v DEFINITION: A non-negaive, irreducible marix is primiive if is dominan eigenvalue is sricly dominan. The Fundamenal Theorem of Demography Suppose he m x m marix P is non-negaive & irreducible. Le > be he dominan eigenvalue of P and v > be an associaed posiive eigenvecor. Le x ( be he soluion of x ( = Px ( Then λ λ λ x( = x ( x. > lim x ( = (unbounded growh < lim x ( = (exincion Suppose P is primiive. Then x ( v lim =. x ( v 5
6 How o deermine if a non-negaive marix is irreducible? and primiive? Some useful facs: An irreducible, non-negaive marix is primiive if is race is posiive. A non-negaive marix is irreducible and primiive if and only if n P > for some posiive ineger n. There are many ess for specialized marices. L Example: Leslie Marices f f fm fm τ = τ 3 τmm, τ mm fi, τ mm, < τ ii, for i =,3,, m. Suppose m is an adul class : fm >. Then L is irreducible.. If, in addiion, race L = f τ mm >, hen L is primiive. 3. Bu i is common in models ha f = (juvenile class τ = (maximum age class. m L is primiive if wo consecuive classes are ferile, i. e., if wo consecuive f i are nonzero. Example: Usher Marices f τ f f3 fm fm τ τ fi, < τ i, i τ3 τ33 for i =,3,, m U = τii, τii τi, i τ m, m for i =,,, m τmm, τ mm Commonly used as a size-srucured model: individuals eiher remain in a size class or move one ahead each ime sep.. Suppose m is an adul class : fm >. Then L is irreducible.. If, in addiion, a leas one τ ii >, hen race U > and U is primiive. x ( = Px ( P= F T, F = f, T = τ ( ij ( ij f ij per capia j-class newborns from i-class individuals ij fracion of j-class individuals ha move o i-class m ij, ij i Suppose P > is irreducible and primiive. Then he Fundamenal Theorem of Demography applies and he normalized disribuion sabilizes. 6
7 Problems: How o calculae he dominan eigenvalue λ? Deermine when λ saisfies λ < or λ >? How does λ depend on he model parameers? EXAMPLE: x Leslie Marix L f f, = τ < τ τ ( ( ( 4 λ = f τ f τ fτ f τ Rewrie he inequaliy λ < as ( ( f τ f τ 4 fτ f τ < ( f τ 4( fτ f τ < ( f τ ( f τ 4( fτ f τ < 4 4( f τ ( f τ f fτ < τ Define n= f f τ n < λ < Similar calculaions show n > λ > n = λ = n has a biological inerpreaion : f Number of offspring when of age τ f f τ Define he number of offspring when of age x he probabiliy of reaching age n= f f τ n < λ < Similar calculaions show n > λ > n = λ = n has a biological inerpreaion : τ 7
8 Define f f τ f τ τ n= f f τ he number of offspring when of age x he probabiliy of reaching age x he probabiliy of surviving one more ime uni n < λ < Similar calculaions show n > λ > n = λ = n has a biological inerpreaion : τ Define n= f f τ n < λ < Similar calculaions show n > λ > n = λ = n has a biological inerpreaion : f f τ f τ τ f τ τ τ he number of offspring when of age x he probabiliy of reaching age x he probabiliy of surviving wo more ime unis Define n= f f τ n < λ < Similar calculaions show n > λ > n = λ = n has a biological inerpreaion : τ f f τ f τ τ f τ τ f τ τ 3 he number of offspring when of age x he probabiliy of reaching age x he probabiliy of surviving hree more ime unis Define n= f f τ n < λ < Similar calculaions show n > λ > n = λ = n has a biological inerpreaion : τ f f τ f τ τ f τ τ f τ τ 3 = expeced number of offspring during lifeime = ne reproducive number 8
9 Define n= f f τ n < λ < Similar calculaions show n > λ > n = λ = n has a biological inerpreaion : f f τ f τ τ f τ τ f τ τ 3 3 ( = f fτ τ τ τ = f fτ τ τ n = ne reproducive number General Definiion of he Ne Reproducive Number The eigenvalues of a marix ( T = τ, τ, τ ij ij ij i= saisfy λ <. We assume all eigenvalues saisfy λ <. Tha is o say, he specral radius ρ(t of T saisfies ρ(t <. Then he inverse exiss. m ( 3 I T I T T T = - > DEFINITION: If F( I T has a dominan eigenvalue n, hen n is called he ne reproducive number of P = F T. THEOREM: Suppose P = F T is irreducible and le λ > denoe is dominan eigenvalue. Suppose ( F, ρ( T <, F I T has a dominan eigenvalue n >. Then one of he following holds λ < n<, < n< λ, n= λ =. This is a sligh improvemen of a heorem of JC & Z. Yicang proved by Li & Schneider (J. Mah. Biol. 44, no. 5 (, AN EXAMPLE: he X Leslie marix (revisied L f f, τ = < τ τ f f =, = ( = τ < F T ρ T τ τ f f F( I T = τ τ f fτ f = τ τ Has a sricly dominan eigenvalue n= f f τ τ. 9
10 ij Biological Inerpreaion of n e e em e e e m 3 I T = = I T T T em em emm ( e = he fracion of ime an j - class newborn spends in class i Suppose here is only one newborn class: n * * ( F I T = f f f m F = n = f e f e f e m m For an inerpreaion of n when here is more han one newborn class see he monograph erraa. Selec a disribuion of newborns (among he classes. Calculae he expeced per capia number of offspring from each newborn class and ake he minimum. n is he maximum of hese minima aken over all possible disribuions of newborns. (The max and min can be inerchanged. n and λ can be used o sudy asympoic dynamics of populaions One advanage for n is ha formulas ofen exis ha relae i o he enries in he marices F and T. This is almos never he case for λ. Examples appear in he Monograph. f τ f f3 fm fm τ τ Usher Marix : τ3 τ33 (for size srucured populaions U = τ m, m τmm, τ mm Assume τ ii <. Then m i τ j, j n= f i ( where τ = i= j= τ jj An Applicaion Hudson River sriped bass (P. H. Leslie Populaion is divided ino 6 age classes: -yr-old (eggs, -yr-old (larvae -yr-old,, 4-yr-old (juveniles 5-yr-old,, 5-yr-old (aduls Ages 3 5 are subjec o fishing
11 The Leslie marix Can calculae λ numerically (using a compuer. n= f f τ fτ τ f m τ τ τ m m 3 3 3, ( τ τ τ f ( τ τ τ = f ,5 =.999 < λ < populaion is going exinc (I urns ou ha λ = Formulas for n can make parameer sudies possible. Suppose fishing moraliy could be managed in such a way ha he wo oldes ages are no aken. Suppose his resuls in an increased survivorship of s for hese older aduls. Quesion: wha survivorship s, if any, will keep he populaion from going exinc? ( τ τ τ ( τ τ τ ( τ τ f ( τ τ τ f5 ( ττ3 τ4,3τ5,4 f6 ( ττ3 τ4,3τ5,4τ6,5 ( τ τ f ( τ τ τ f5 ( ττ3 τ4,3s f6 ( ττ3 τ4,3s n= f f ,5 = f ,3 = f ,3 3 n= s s 3 n> if s > No possible.
Vehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More information1 Review of Zero-Sum Games
COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any
More informationStationary Distribution. Design and Analysis of Algorithms Andrei Bulatov
Saionary Disribuion Design and Analysis of Algorihms Andrei Bulaov Algorihms Markov Chains 34-2 Classificaion of Saes k By P we denoe he (i,j)-enry of i, j Sae is accessible from sae if 0 for some k 0
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationApplication 5.4 Defective Eigenvalues and Generalized Eigenvectors
Applicaion 5.4 Defecive Eigenvalues and Generalized Eigenvecors The goal of his applicaion is he soluion of he linear sysems like where he coefficien marix is he exoic 5-by-5 marix x = Ax, (1) 9 11 21
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationChapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws
Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationExercises: Similarity Transformation
Exercises: Similariy Transformaion Problem. Diagonalize he following marix: A [ 2 4 Soluion. Marix A has wo eigenvalues λ 3 and λ 2 2. Since (i) A is a 2 2 marix and (ii) i has 2 disinc eigenvalues, we
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationChapter Three Systems of Linear Differential Equations
Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More information5.2. The Natural Logarithm. Solution
5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,
More informationMath 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.
Mah 36. Rumbos Spring 1 1 Soluions o Assignmen #6 1. Suppose he growh of a populaion is governed by he differenial equaion where k is a posiive consan. d d = k (a Explain why his model predics ha he populaion
More informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationChallenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k
Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,
More informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationConcourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations
Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial
More informationComparing Means: t-tests for One Sample & Two Related Samples
Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More informationIntroduction to Probability and Statistics Slides 4 Chapter 4
Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random
More information5. Stochastic processes (1)
Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly
More informationApproximation Algorithms for Unique Games via Orthogonal Separators
Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define
More informationV L. DT s D T s t. Figure 1: Buck-boost converter: inductor current i(t) in the continuous conduction mode.
ECE 445 Analysis and Design of Power Elecronic Circuis Problem Se 7 Soluions Problem PS7.1 Erickson, Problem 5.1 Soluion (a) Firs, recall he operaion of he buck-boos converer in he coninuous conducion
More informationSystem of Linear Differential Equations
Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y
More informationKEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationSPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),
SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric
More informationWe just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n
Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.
More informationSolutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore
Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationNotes for Lecture 17-18
U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationLet ( α, β be the eigenvector associated with the eigenvalue λ i
ENGI 940 4.05 - Sabiliy Analysis (Linear) Page 4.5 Le ( α, be he eigenvecor associaed wih he eigenvalue λ i of he coefficien i i) marix A Le c, c be arbirary consans. a b c d Case of real, disinc, negaive
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationChapter Floating Point Representation
Chaper 01.05 Floaing Poin Represenaion Afer reading his chaper, you should be able o: 1. conver a base- number o a binary floaing poin represenaion,. conver a binary floaing poin number o is equivalen
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationMath Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.
Mah 50-004 Week 5: Secion 7.4, mass-spring sysems. These are noes for Monday. There will also be course review noes for Tuesday, posed laer. Mon Apr 3 7.4 mass-spring sysems. Announcemens: Warm up exercise:
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationTHE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX
J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he
More informationProblem Set on Differential Equations
Problem Se on Differenial Equaions 1. Solve he following differenial equaions (a) x () = e x (), x () = 3/ 4. (b) x () = e x (), x (1) =. (c) xe () = + (1 x ()) e, x () =.. (An asse marke model). Le p()
More informationBiol. 356 Lab 8. Mortality, Recruitment, and Migration Rates
Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese
More informationOn linear ODEs with a time singularity of the first kind and unsmooth inhomogeneity
Rachůnková e al. RESEARCH On linear ODEs wih a ime singulariy of he firs kind and unsmooh inhomogeneiy Irena Rachůnková, Svaoslav Saněk, Jana Vampolová and Ewa B Weinmüller 2* * Correspondence: ewa.weinmueller@uwien.ac.a
More informationMacroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3
Macroeconomic Theory Ph.D. Qualifying Examinaion Fall 2005 Comprehensive Examinaion UCLA Dep. of Economics You have 4 hours o complee he exam. There are hree pars o he exam. Answer all pars. Each par has
More informationAge (x) nx lx. Age (x) nx lx dx qx
Life Tables Dynamic (horizonal) cohor= cohor followed hrough ime unil all members have died Saic (verical or curren) = one census period (day, season, ec.); only equivalen o dynamic if populaion does no
More information2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS
Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described
More informationBifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays
Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationT. J. HOLMES AND T. J. KEHOE INTERNATIONAL TRADE AND PAYMENTS THEORY FALL 2011 EXAMINATION
ECON 841 T. J. HOLMES AND T. J. KEHOE INTERNATIONAL TRADE AND PAYMENTS THEORY FALL 211 EXAMINATION This exam has wo pars. Each par has wo quesions. Please answer one of he wo quesions in each par for a
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More information4.1 - Logarithms and Their Properties
Chaper 4 Logarihmic Funcions 4.1 - Logarihms and Their Properies Wha is a Logarihm? We define he common logarihm funcion, simply he log funcion, wrien log 10 x log x, as follows: If x is a posiive number,
More informationOptimality Conditions for Unconstrained Problems
62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More informationDISCRETE GRONWALL LEMMA AND APPLICATIONS
DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE MAA NORTH CENTRAL SECTION MEETING AT UND 24 OCTOBER 29 Gronwall s lemma saes an inequaliy ha is useful in he heory of differenial equaions. Here is
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationTracking Adversarial Targets
A. Proofs Proof of Lemma 3. Consider he Bellman equaion λ + V π,l x, a lx, a + V π,l Ax + Ba, πax + Ba. We prove he lemma by showing ha he given quadraic form is he unique soluion of he Bellman equaion.
More informationLecture 4 Notes (Little s Theorem)
Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,
More informationOn Oscillation of a Generalized Logistic Equation with Several Delays
Journal of Mahemaical Analysis and Applicaions 253, 389 45 (21) doi:1.16/jmaa.2.714, available online a hp://www.idealibrary.com on On Oscillaion of a Generalized Logisic Equaion wih Several Delays Leonid
More information4.6 One Dimensional Kinematics and Integration
4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of
More informationOn Boundedness of Q-Learning Iterates for Stochastic Shortest Path Problems
MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 2, May 2013, pp. 209 227 ISSN 0364-765X (prin) ISSN 1526-5471 (online) hp://dx.doi.org/10.1287/moor.1120.0562 2013 INFORMS On Boundedness of Q-Learning Ieraes
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More informationThe Brock-Mirman Stochastic Growth Model
c December 3, 208, Chrisopher D. Carroll BrockMirman The Brock-Mirman Sochasic Growh Model Brock and Mirman (972) provided he firs opimizing growh model wih unpredicable (sochasic) shocks. The social planner
More informationECON 482 / WH Hong Time Series Data Analysis 1. The Nature of Time Series Data. Example of time series data (inflation and unemployment rates)
ECON 48 / WH Hong Time Series Daa Analysis. The Naure of Time Series Daa Example of ime series daa (inflaion and unemploymen raes) ECON 48 / WH Hong Time Series Daa Analysis The naure of ime series daa
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationKEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More information. Now define y j = log x j, and solve the iteration.
Problem 1: (Disribued Resource Allocaion (ALOHA!)) (Adaped from M& U, Problem 5.11) In his problem, we sudy a simple disribued proocol for allocaing agens o shared resources, wherein agens conend for resources
More informationMODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE
Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS
More informationMath 333 Problem Set #2 Solution 14 February 2003
Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial
More informationEssential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems
Essenial Microeconomics -- 6.5: OPIMAL CONROL Consider he following class of opimizaion problems Max{ U( k, x) + U+ ( k+ ) k+ k F( k, x)}. { x, k+ } = In he language of conrol heory, he vecor k is he vecor
More informationt + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that
ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationOptimal Paired Choice Block Designs. Supplementary Material
Saisica Sinica: Supplemen Opimal Paired Choice Block Designs Rakhi Singh 1, Ashish Das 2 and Feng-Shun Chai 3 1 IITB-Monash Research Academy, Mumbai, India 2 Indian Insiue of Technology Bombay, Mumbai,
More informationLars Nesheim. 17 January Last lecture solved the consumer choice problem.
Lecure 4 Locaional Equilibrium Coninued Lars Nesheim 17 January 28 1 Inroducory remarks Las lecure solved he consumer choice problem. Compued condiional demand funcions: C (I x; p; r (x)) and x; p; r (x))
More informationHypothesis Testing in the Classical Normal Linear Regression Model. 1. Components of Hypothesis Tests
ECONOMICS 35* -- NOTE 8 M.G. Abbo ECON 35* -- NOTE 8 Hypohesis Tesing in he Classical Normal Linear Regression Model. Componens of Hypohesis Tess. A esable hypohesis, which consiss of wo pars: Par : a
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationProblem set 2 for the course on. Markov chains and mixing times
J. Seif T. Hirscher Soluions o Proble se for he course on Markov chains and ixing ies February 7, 04 Exercise 7 (Reversible chains). (i) Assue ha we have a Markov chain wih ransiion arix P, such ha here
More information12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =
1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationThe Arcsine Distribution
The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More informationON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES WITH IMMIGRATION WITH FAMILY SIZES WITHIN RANDOM INTERVAL
ON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES ITH IMMIGRATION ITH FAMILY SIZES ITHIN RANDOM INTERVAL Husna Hasan School of Mahemaical Sciences Universii Sains Malaysia, 8 Minden, Pulau Pinang, Malaysia
More information