A Primer on Continuous-time Economic Dynamics

Similar documents
Section 2.3. Matrix Inverses

Matrix- System of rows and columns each position in a matrix has a purpose. 5 Ex: 5. Ex:

Numbers and indices. 1.1 Fractions. GCSE C Example 1. Handy hint. Key point

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

SIMPLE NONLINEAR GRAPHS

50 AMC Lectures Problem Book 2 (36) Substitution Method

18.06 Problem Set 4 Due Wednesday, Oct. 11, 2006 at 4:00 p.m. in 2-106

Reference : Croft & Davison, Chapter 12, Blocks 1,2. A matrix ti is a rectangular array or block of numbers usually enclosed in brackets.

Eigenvectors and Eigenvalues

MCH T 111 Handout Triangle Review Page 1 of 3

Factorising FACTORISING.

Differentiation of Polynomials

Lecture 6: Coding theory

CS 491G Combinatorial Optimization Lecture Notes

6.5 Improper integrals

Matrices SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics (c) 1. Definition of a Matrix

Solutions for HW9. Bipartite: put the red vertices in V 1 and the black in V 2. Not bipartite!

For a, b, c, d positive if a b and. ac bd. Reciprocal relations for a and b positive. If a > b then a ab > b. then

F / x everywhere in some domain containing R. Then, + ). (10.4.1)

1 PYTHAGORAS THEOREM 1. Given a right angled triangle, the square of the hypotenuse is equal to the sum of the squares of the other two sides.

Surds and Indices. Surds and Indices. Curriculum Ready ACMNA: 233,

for all x in [a,b], then the area of the region bounded by the graphs of f and g and the vertical lines x = a and x = b is b [ ( ) ( )] A= f x g x dx

M344 - ADVANCED ENGINEERING MATHEMATICS

LINEAR ALGEBRA APPLIED

AP Calculus BC Chapter 8: Integration Techniques, L Hopital s Rule and Improper Integrals

Linear Algebra Introduction

Calculus AB. For a function f(x), the derivative would be f '(

Logarithms LOGARITHMS.

Chapter Gauss Quadrature Rule of Integration

Let s divide up the interval [ ab, ] into n subintervals with the same length, so we have

Mid-Term Examination - Spring 2014 Mathematical Programming with Applications to Economics Total Score: 45; Time: 3 hours

22: Union Find. CS 473u - Algorithms - Spring April 14, We want to maintain a collection of sets, under the operations of:

Solutions to Assignment 1

Proportions: A ratio is the quotient of two numbers. For example, 2 3

Counting Paths Between Vertices. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs. Isomorphism of Graphs

Now we must transform the original model so we can use the new parameters. = S max. Recruits

Algebra 2 Semester 1 Practice Final

Project 6: Minigoals Towards Simplifying and Rewriting Expressions

CIT 596 Theory of Computation 1. Graphs and Digraphs

Math 32B Discussion Session Week 8 Notes February 28 and March 2, f(b) f(a) = f (t)dt (1)

I 3 2 = I I 4 = 2A

ax bx c (2) x a x a x a 1! 2!! gives a useful way of approximating a function near to some specific point x a, giving a power-series expansion in x

The Double Integral. The Riemann sum of a function f (x; y) over this partition of [a; b] [c; d] is. f (r j ; t k ) x j y k

GRAND PLAN. Visualizing Quaternions. I: Fundamentals of Quaternions. Andrew J. Hanson. II: Visualizing Quaternion Geometry. III: Quaternion Frames

The Ellipse. is larger than the other.

Here we consider the matrix transformation for a square matrix from a geometric point of view.

5. Every rational number have either terminating or repeating (recurring) decimal representation.

April 8, 2017 Math 9. Geometry. Solving vector problems. Problem. Prove that if vectors and satisfy, then.

Calculus Module C21. Areas by Integration. Copyright This publication The Northern Alberta Institute of Technology All Rights Reserved.

Worksheet #2 Math 285 Name: 1. Solve the following systems of linear equations. The prove that the solutions forms a subspace of

x ) dx dx x sec x over the interval (, ).

Lesson 2.1 Inductive Reasoning

Lesson 2.1 Inductive Reasoning

Section 4: Integration ECO4112F 2011

Green s Theorem. (2x e y ) da. (2x e y ) dx dy. x 2 xe y. (1 e y ) dy. y=1. = y e y. y=0. = 2 e

Sturm-Liouville Theory

Integration. antidifferentiation

Introduction to Algebra - Part 2

] dx (3) = [15x] 2 0

AP Calculus AB Unit 4 Assessment

y z A left-handed system can be rotated to look like the following. z

( ) { } [ ] { } [ ) { } ( ] { }

CHAPTER 4: DETERMINANTS

University of Sioux Falls. MAT204/205 Calculus I/II

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

Solids of Revolution

Exercise sheet 6: Solutions

H (2a, a) (u 2a) 2 (E) Show that u v 4a. Explain why this implies that u v 4a, with equality if and only u a if u v 2a.

Particle Physics. Michaelmas Term 2011 Prof Mark Thomson. Handout 3 : Interaction by Particle Exchange and QED. Recap

The Trapezoidal Rule

QUADRATIC EQUATION. Contents

Lecture 1 - Introduction and Basic Facts about PDEs

September 13 Homework Solutions

Topics Covered AP Calculus AB

K 7. Quadratic Equations. 1. Rewrite these polynomials in the form ax 2 + bx + c = 0. Identify the values of a, b and c:

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Topic 1 Notes Jeremy Orloff

Chapter 16: Magnetic Materials and a Few Devices

APPENDIX. Precalculus Review D.1. Real Numbers and the Real Number Line

Mathematics Number: Logarithms

Chapter Five - Eigenvalues, Eigenfunctions, and All That

sec x over the interval (, ). x ) dx dx x 14. Use a graphing utility to generate some representative integral curves of the function Curve on 5

MULTIPLE CHOICE QUESTIONS

Calculus Cheat Sheet. Integrals Definitions. where F( x ) is an anti-derivative of f ( x ). Fundamental Theorem of Calculus. dx = f x dx g x dx

INTEGRATION. 1 Integrals of Complex Valued functions of a REAL variable

Math 211A Homework. Edward Burkard. = tan (2x + z)

Chapter 6 Techniques of Integration

21.1 Using Formulae Construct and Use Simple Formulae Revision of Negative Numbers Substitution into Formulae

Logarithms. Logarithm is another word for an index or power. POWER. 2 is the power to which the base 10 must be raised to give 100.

12.4 Similarity in Right Triangles

Lesson 55 - Inverse of Matrices & Determinants

Applications of Definite Integral

Data Structures LECTURE 10. Huffman coding. Example. Coding: problem definition

Partial Differential Equations

Bases for Vector Spaces

y = c 2 MULTIPLE CHOICE QUESTIONS (MCQ's) (Each question carries one mark) is...

Core 2 Logarithms and exponentials. Section 1: Introduction to logarithms

Mathematics. Area under Curve.

Lecture 2: Cayley Graphs

16z z q. q( B) Max{2 z z z z B} r z r z r z r z B. John Riley 19 October Econ 401A: Microeconomic Theory. Homework 2 Answers

Transcription:

Eonomis 205A Fll 2008 K Kletzer A Primer on Continuous-time Eonomi Dnmis A Liner Differentil Eqution Sstems (i) Simplest se We egin with the simple liner first-orer ifferentil eqution The generl solution is ẋ (0) 0 0 e t n the initil onition is stisfie the prtiulr solution The growth rte of is given n 0 e t Net we solve the two-imensionl sstem given ẋ (0) 0 ẏ (0) 0 It is useful to write these in mtri form s ẋẏ 0 (0) 0 (0) 0 0 The solution for this two-imensionl sstem is simpl 0 e t 0 e t Eh of these prolems the one-imensionl n the two-imensionl re emples of initilvlue prolems This sstem possesses single ste stte 0 0n 0 0 A phse igrm plots out s funtion of for ifferent possile initil vlues of n The resulting lous relting to t eh time given n prtiulr initil vlues 0 n 0 epits trjetor long whih n move over time Arrows re usull rwn to epit the iretion of motion in the -plne s t inreses Consier three ses In the first the prmeters n re othpositive For ninitilpoint

( 0 0 ) tht is not the ste stte trjetor moves w from the origin This is the unstle se In the seon n re oth negtive n n trjetor onverges to the ste stte This is the stle se In this se trjetories smptotill pproh the ste stte prllel to the -is if <<0 n onversel In the thir se one root or is positive n the other is negtive For the emple in whih >0 > n trjetor suh tht 0 0onverges to the ste stte t the rte An trjetor suh tht 0 0moves w from the ste stte n onverges towrs the -is smptotill This is the sle-pth stle se Consier n emple of sle-pth nmis suh tht n The solution is given 0 et 0 e t As t grows without oun pprohes zero s grows without oun s 0e t (ii) Generl se Let z e n n-imensionl olumn vetor over the rels n ż Az where z (0) z 0 The mtri A isn n with onstnt rel oeffiients Impose the restrition tht A is nonsingulr For emple onsier the two-imensionl ifferentil eqution sstem ẋẏ (0) (0) We hnle this prolem hnging oorintes We wnt to fin hnge of sis from new sis suh tht the initil-vlue prolem eomes ŷ ẋ t â 0 0 whih hs the solution given in prt (i) ŷ ŷ ât 0 e t ŷ 0 e (0) ŷ (0) 0 0 0 This is one fining nonsingulr 2 2 mti M with onstnt oeffiients suh tht â 0 M n M M We n lws o this given the ssumptions me so fr The ifferentil eqution sstems re ientil ut epresse in ifferent oorintes ŷ 0 ŷ 0 2 to

3 The solution for the initil-vlue prolem ẋẏ (0) (0) 0 0 is given M ât 0 e So the onl issue left is fining M Strt oserving tht M les to where t ŷ 0 e â 0 M 0 This n e rewritten s two equtions m m 2 M where M â âm âm 2 m m 2 m 2 m 22 m m 2 â 0 ŷ 0 M m2 m22 0 0 M M m m 2 0 0 â 0 0 n m 2 m 22 m 2 m 22 These impl tht in the generl se we seek solutions to 0 0 m 2 m 22 A λi ν 0 where I is the ientit mtri of sme imension s A λ is slr n ν is olumn vetor For A nonsingulr this onl hs solutions suh tht et A λi 0 whih is n th orer polnomil You just solve for the n vlues of λ λ λ n n then sustitute to fin the orresponing n olumn vetors ν j M is the mtri me up of the n olumn vetors The roots λ λ n n e omple numers If root is omple then its omple onjugte is lso root (omple roots re isusse in susetion (iii)) For the two-imensionl se we hve tht â λ λ 2 ν The generl solutions for n re given α ν e λ t + α 2 ν 2 e λ2t m m 2 n ν 2 m 2 m 22

with the strighforwr etension to the n imensionl se Here α n α 2 re solutions to (0) (0) α ν + α 2 ν 2 We will use the onventionl terminolog: λ is lle n eigenvlue n ν is lle n eigenvetor et A λi 0is lle the hrteristi polnomil An emple: let 6 A 2 The eigenvlues n ssoite eigenvetors re λ 4 λ 2 ν ν 2 Note tht n non-zero multiple of these eigenvetors is lso n eigenvetor The generl solutions for n re α 2 2 e 4t + α2 3 e t The initil onitions re stisfie hoies of α n α 2 tht stisf + α2 whih implies (iii) Comple roots α 2 α 2 0 6 0 5 5 3 0 0 n α 2 3 0 5 + 6 0 5 The roots of the hrteristi polnomil for the n-imensionl liner initil-vlue prolem n e omple numers For emple onsier the sstem whih hs the hrteristi polnomil with the two roots ẋ ẏ 0 0 λ 2 +0 λ i n λ 2 i The eigenvetors ssoite with eh of these eigenvlues re i ν n ν 2 The generl solution is given the rel prt of α e it + α2 i i i e it 3 4

5 where e iθt os θt + i sin θt Sustituting we hve the generl solution α sin t os t + α2 os t sin t n the prtiulr solution sin t os t 0 + 0 os t sin t whih ou n hek setting t to zero in the generl solution For this emple ll trjetories re ounterlokwise irles roun the origin prts In the generl se in two imensions we might hve roots tht re omple with non-zero rel The generl solution will e the rel prt of whih equls λ γ + θi n λ 2 γ θi α ν e (γ+θi)t + α 2 ν 2 e (γ θi)t e γt α ν e θit + α 2 ν 2 e θit The prt insie the squre rkets is simpl irulr motion (lokwise or ounterlokwise epening on the oeffiients of the mtri A) The eponentil in front epns or ontrts these les epening on whether γ is positive or negtive B Non-liner nmil sstems In mn eonomis pplitions the nmis will not e represente liner ifferentil equtions In tpil ses we n onl prtill hrterize the nmis of our moel lthough in some speil instnes ou n integrte the sstem An emple of non-liner moel is the stnr one-setor optiml growth moel with iminishing returns to pitl The nmis re given k f (k) nk n ċ σ () (f (k) n ρ)

for σ () u () The initil onitions re etermine k u () 0 n the trnsverslit onition 6 lim u ( t e ρt t ) k t 0 We nnot integrte the eqution for ċ unless oth the intertemporl elstiit of sustitution n mrginl proutivit of pitl re onstnt One pproh to hrterizing the nmis of pitl umultion is to emine these ner ste stte We impose the In onitions so tht ste stte with positive k n eists First the sstem n e pproimte ner the ste stte tking first-orer Tlor series epnsion out the ste stte: n ċ f (k) σ ( ) k f (k) n(k k ) ( ) (k k )+(f (k ) n ρ) σ ( ) σ ( ) σ 2 ( ) ( ) We will el with the prtiulr solution tht stisfies the initil n trnsverslit onitions fter we hrterize the generl solutions to the two-imensionl linerize ifferentil eqution sstem n Notie tht this emple simplifies to k ρ (k k ) ( ) ċ f (k ) σ ( ) (k k ) when we impose the ste-stte onitions ċ 0n k 0 In mtri form the pproimte sstem is ċ k 0 f (k ) ρ σ( ) ( ) (k k ) where the oeffiients in the mtri A re ll onstnt The hrteristi polnomil is whih hs solutions In this se λ ± 2 λ (ρ λ) + f (k ) ( ρ ± σ ( ) 0 ρ 2 4 f (k ) σ ( ) λ+ >ρ n λ < 0 n the linerize sstem is sle-pth stle in neighorhoo of the ste stte ( k )There is n eigenvetor ssoite with eh of these eigenvlues The one ssoite with the negtive )

7 eigenvlue will give us trjetor for ( k k ) tht onverges to the ste stte This trjetor will e the onl one tht stisfies the trnsverslit onition n the itionl fesiilit onition tht k t 0 for ll t Wht is the slope of this eigenvetor? Our hnge of oorintes ove le us to the eqution Aν j λ j ν j 0 for eh eigenvlue n ssoite eigenvetor So we hve the equtions n where f (k ) ν 2 λ ν σ ( ) ν + ρν 2 λ ν 2 ν For λ λ these impl tht the slope of the stle sle pth ν /ν 2 ρ λ whihis positive n steeper thn the slope of the k 0lous t the ste stte in the phse igrm The slope of the unstle sle pth ρ λ + is negtive euse λ + >ρ It is importnt to know tht the solution to the linerize ifferentil eqution sstem pproimtes the nmis of the originl non-liner sstem ner the ste stte A si theorem in nmil sstems proves tht in neighorhoo of ste stte the trjetories for the generl solution of the lineriztion of the ifferentil eqution sstem out the ste stte re ritrril lose (in the epsilon-elt sense) to the trjetories of the non-liner sstem uner generl onitions For ontinuousl ifferentile funtions u () n f (k) the trjetories of the linerize n non-liner sstem re tngent to eh other t the ste stte This mens tht if ou fin the slope of the stle sle pth for the liner pproimtion (the slope of the stle root eigenvetor) ou hve foun the slope of the sle pth for the non-liner sstem t the ste stte It lso mens tht the growth rtes of ( ) n (k k ) re pproimtel equl to λ long the stle sle pth in neighorhoo of the ste stte The most importnt onition for this theorem to hol is tht the eigenvlues of the linerize sstem ll hve non-zero rel prts If some of the eigenvlues re omple with zero rel prt then the linerize sstem is not n pproimte sstem n its nmis m hve little in ommon with the nmis of the originl non-liner sstem This is lso n importnt theorem in nmil ν ν 2

8 sstems In generl non-liner sstem n hve mn ste sttes When ou n pproimte its nmis lineriztion of the ifferentil eqution sstem (the ommon se) ou will e oing so roun prtiulr ste stte The liner sstem will hve unique ste stte so ou will hve ifferent liner pproimtion for eh ste stte when there is more thn one When the linerize sstem is stle we ll the originl sstem loll stle When it is unstle we ll the non-liner sstem loll unstle n so forth for lol sle-pth stilit Sometimes we n s something or even lot out the nmis of non-liner sstem w from the ste stte ut we nnot o so using lol pproimtions