= ( ) ) or a system of differential equations with continuous parametrization (T = R

Similar documents
( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

Math 333 Problem Set #2 Solution 14 February 2003

Chapter 2. First Order Scalar Equations

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Predator - Prey Model Trajectories and the nonlinear conservation law

SOLUTIONS TO ECE 3084

Two Coupled Oscillators / Normal Modes

3.6 Derivatives as Rates of Change

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Problem Set on Differential Equations

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Solutions Problem Set 3 Macro II (14.452)

Some Basic Information about M-S-D Systems

Solutions to Assignment 1

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

System of Linear Differential Equations

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

This document was generated at 7:34 PM, 07/27/09 Copyright 2009 Richard T. Woodward

Echocardiography Project and Finite Fourier Series

Economics 8105 Macroeconomic Theory Recitation 6

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

Macroeconomic Theory Ph.D. Qualifying Examination Fall 2005 ANSWER EACH PART IN A SEPARATE BLUE BOOK. PART ONE: ANSWER IN BOOK 1 WEIGHT 1/3

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Appendix 14.1 The optimal control problem and its solution using

5.1 - Logarithms and Their Properties

Fishing limits and the Logistic Equation. 1

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

KINEMATICS IN ONE DIMENSION

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

Seminar 4: Hotelling 2

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Matlab and Python programming: how to get started

EXERCISES FOR SECTION 1.5

Unit Root Time Series. Univariate random walk

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

The average rate of change between two points on a function is d t

15. Vector Valued Functions

Differential Equations

10. State Space Methods

Phys1112: DC and RC circuits

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Math 2214 Solution Test 1A Spring 2016

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

Math 10B: Mock Mid II. April 13, 2016

This document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC

t 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...

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

4.1 - Logarithms and Their Properties

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

Cooperative Ph.D. Program in School of Economic Sciences and Finance QUALIFYING EXAMINATION IN MACROECONOMICS. August 8, :45 a.m. to 1:00 p.m.

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Math 2214 Solution Test 1B Fall 2017

The Brock-Mirman Stochastic Growth Model

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

Solutionbank Edexcel AS and A Level Modular Mathematics

Math 115 Final Exam December 14, 2017

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

Chapter 7: Solving Trig Equations

13.3 Term structure models

Chapter 3 Boundary Value Problem

Announcements: Warm-up Exercise:

1. An introduction to dynamic optimization -- Optimal Control and Dynamic Programming AGEC

3.1 More on model selection

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

Today: Falling. v, a

Basic Circuit Elements Professor J R Lucas November 2001

ES.1803 Topic 22 Notes Jeremy Orloff

4.6 One Dimensional Kinematics and Integration

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

CHAPTER 12 DIRECT CURRENT CIRCUITS

EE100 Lab 3 Experiment Guide: RC Circuits

Lecture Notes 3: Quantitative Analysis in DSGE Models: New Keynesian Model

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

BU Macro BU Macro Fall 2008, Lecture 4

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Instructor: Barry McQuarrie Page 1 of 5

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

Y 0.4Y 0.45Y Y to a proper ARMA specification.

The expectation value of the field operator.

From Complex Fourier Series to Fourier Transforms

Notes 04 largely plagiarized by %khc

AP Chemistry--Chapter 12: Chemical Kinetics

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

Math Final Exam Solutions

Transcription:

XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of ime in one of wo formas: x= or x x x = () The variable belongs o some index se T, usually he nonnegaive inegers or nonnegaive reals. When we alk abou how hese sysems evolve wih ime, we have a sysem of difference equaions when he paramerizaion is discree ( T = Z + ) or a sysem of differenial equaions wih coninuous paramerizaion (T = R + ). These look like, respecively: x = + f( x) and: dx d = g( x) Economiss ofen prefer o work wih difference equaions. In your firs-year macro classes, you will see many of hese. We look a how variables in he economy progress from one year o he nex, which we migh describe wih he sysem of equaions: c f c k, k = + ( ) = + g( c, k ) Nex year s consumpion can be deermined by a funcion of oday s consumpion and oday s capial sock. The same goes for nex year s capial sock. A firs order difference equaion conains only variables of one ime period on he righ-hand side; in conras, a second order difference equaion conains wo periods: c f c, k, c, k = ( ) + Similarly, an n-h order difference equaion has he sae of variables for he pas n periods deermining he nex value. When I wrie diff. eq.s (pronounced diffy queues ) I am saing somehing abou boh difference and differenial equaions. Diff. eq.s are also known as laws of moion or evoluionary laws. Firs-order linear differenial equaions are he simples class of difference equaions o work wih. They look like: x+ = ax + b The differenial equaion is called homogeneous when b and ime-invarian when a = a, b = b are consans. Yeserday we looked a a sysem of homogeneous, ime-invarian, firs-order linear differenial equaions. Here s a simple one: k = k k = + δ ( δ) k Nex year s capial sock is he same as his year s, minus some consan depreciaion. Then we can forecas he capial sock wo years in he fuure: k = δ k δ k ( ) = ( ) + + In fac, we can say ha he capial sock n years ino he fuure will be: n ( ) = ( ) + n k+ n= δ k δ k0 Summer 00 mah class noes, page 9

where k 0 is he value of he capial sock a he beginning of ime. Firs-order diff. eq.s can be solved subjec o one iniial value. Tha is, provided we know wha he capial sock was a =0, we can predic i for he res of eerniy. The given value doesn necessarily have o be iniial we could have been given he value of any k and we could solve he sysem a all. Given he general firs-order difference equaion, he general soluion is: x = as x as br 0 + s r s r = The soluion is much easier o express if we have consans for a and b, bu his isn always he case. Example: I sar my life wih A 0 in asses. Each year, asses grow wih a a rae of ρ. Each year I make invesmens (or wihdrawals, if negaive) of i, which we ake as exogenous (no depending on asses): A = + + ( ρ) A + i Solve his in erms of A 0 and i. Sysems of linear differenial equaions can be expressed in marix form, in one of hese manners: x + a a x b x a a + = x + x ( + ) a a x b or: b x a a ( + ) = () () () () () x() + () b () eiher of which I will wrie compacly as: v v v x = x + + A b The soluion here is: v v v x A A KA A A x A A KA A A b = ( ) + ( ) 0 0 s+ s+ s s s= 0 Example: The capial sock evolves according o some equaion: k+ = ( + ρ n) k + i Whereas populaion growh can be described as: N = + + ( n) N + δk Solve his sysem in erms of k 0, N 0, and i. If he i, wha can you say abou sabiliy? A firs-order linear difference equaion is he bes ha you can hope for. There are ricks for making hem homogeneous when hey are no. I m no going o ge ino hose. Many equaions in macro will be second-order difference equaions. In fac, you ll become very well acquained wih his one: ( ) β ( ( ) + )= ( ( ) ) f k U f k k U f k k Summer 00 mah class noes, page 93

This defines k + as a funcion (a leas implicily) of k and k. A is bes, his equaion will ake he form: φβ( φk k)= φk k+ k+ = φ( + β) k φ βk Second-order differenial equaions are a bich o solve. Even his one, a secondorder, ime-invarian, homogenous linear difference equaion (as nice as hey ge), is ricky enough. Second-order diff. eq.s require wo iniial values o solve (ofen one of hese will be he acual iniial value and he oher he erminal value). We solve hem using a highly scienific mehod known as guess and verify, or more euphemisically, as he mehod of undeermined coefficiens. The idea is ha, given a difference equaion like: x+ = f( x, x ), we can guess ha here is some funcional form g so ha x+ = g( x) and x = g( x ). We inver he second of hese, x g x, sick i ino f, = ( )?? = ( )= ( ( ))= ( ) x f x, x f x, g x g x + and see wheher his guess gives us back ha x g x has he same funcional form. If i does, our guess is good. Le s go back o our nices possible second-order equaion: k+ = φ( + β) k φ βk My guess is going o be ha we can wrie each k as jus a consan imes he preceding one: k = σk This means ha k = σ k, which we subsiue back ino our difference equaion. k+ = φ( + β) k φ β σ k k+ = ( φ( + β) φ β σ ) k If his is rue, hen all he suff in fron of k mus also equal σ. σ = φ( + β φβ σ ) σ σφ φβσ + φ β = ( σ φβ) ( σ φ ) We have deermined ha wo roos possibly saisfy his equaion; his effecively confirms ha he linear guess was correc. We should sill check ha hese are sensible soluions. This was he paricular funcional form ha solved his paricular problem. For oher difference equaions, you need o use oher guesses. If you are familiar wih he paricular funcional forms involved, you migh are more likely o guess correcly. If you re uncerain, check a book like Sydsæer, Srøm, and Berck for some recommendaions. + = ( ) In difference equaions, we re considering funcions ha essenially (afer some rearranging) look like: x x = + f( x) x Summer 00 mah class noes, page 94

Le s call he lef hand side of his equaion x, and define he righ hand side as some funcion g. Then a difference equaion can always be expressed as: x = g( x)( ) The change in ime is = ( + ) =, so muliplying he righ hand side by his changes nohing. This looks like a lo like a firs-order approximaion, which is exacly wha a difference equaion is. Taking he limi of his as 0, we can wrie his as a coninuous ime differenial equaion: dx = g( x)( d) dx d = g( x) Differenial equaions are always expressed a he derivaive of a variable wih respec o ime. Someimes his is denoed by dx d, someimes by x ( ), and very ofen by «x. As wih difference equaions, he simples a firs-order, linear, ime-invarian, homogeneous equaion is he simples differenial equaion possible o solve. These equaions look like: «x = ax x ( )= ax() Wha funcion is is own derivaive? Tha s e y. Acually, i s Ce y, where C is any consan. Wha funcion imes a consan (say, γ ) is is own derivaive? Well, ha s Ce γ y. And ha s exacly he soluion o his differenial equaion. a a «x = ax x = Ce = x 0 e (You should probably remember his formula.) C is called a consan of inegraion, which equals he iniial value x 0 a =0. If you don know an iniial value, you should leave C in as an undeermined coefficien. More generally, given any linear firs order differenial equaion of he form: «x = ax + b where a and b are consans ha do no depend on x, he soluion o his problem should saisfy: x = C a( s) ds b r a s ds dr + ( ) ( ) exp exp 0 0 r Higher order and nonlinear differenial equaions are messy o solve (bu generally, more solvable han difference equaions), and sysems of equaions require some ricks. There are formulas for a few soluions, bu ofen you have o employ a bunch of fun echniques o ge hem o work ou. Go look hem up if you ever need hem (and hopefully you won ). Solving diff. eq.s is difficul, bu forunaely, we re ofen no ineresed in he acual soluions. Consider an economy, he capial sock and consumpion of which can be described by he pair of differenial equaions: «α k = Ak ( δ + n) k c U c c«= c [ ( ) ] α A( α) k ( δ + n) θ [ ] cu ( c ) Summer 00 mah class noes, page 95

Being able o solve hese would allow us o say, given ha we know Albania s consumpion and capial sock oday, in he year 086 hey will be such-and-such. This may be imporan in iself. However, no being able o solve hem explicily, we can sill answer many quesions abou he growh (or decline) of he variables. Growh of capial sock and sandards of living are fundamenal issues in developmen. Growh of money socks and real money balances are cenral o moneary heory. A sock variable or sae is like he variable x. In conras, he variable «x is called a flow variable. Think of a pond being filled wih waer. The sock a any ime is how much waer is in he pond (gallons or cubic meers of waer); he flow is he rae a which he seam feeds ino i (measured in gallons per minue). When looking a a pair of differenial equaions in which each sock is an inpu ino he flow of each variable (as wih capial and consumpion, above), we ofen consruc phase diagrams in order o idenify regions in which each is growing or declining. There are a number of seps in his.. Draw a wo-dimensional graph, wih each sock variable on one axis.. For each variable, idenify under wha condiions here is no change in he sock. You do his by seing «x and solving he equaion for one variable in erms of he oher. Plo his funcion (or correspondence). This is called an x-consan locus. 3. The x-consan locus or loci separae he graph ino several regions. In each of hese regions, figure ou wheher x is increasing or decreasing (all you need do is evaluae a one poin in each region). In he region, skech in some arrows parallel o he x-axis poining in he appropriae direcion. 4. Overlay he wo graphs. You ve now picked ou regions where x is growing and y is growing, x is growing and y is declining, and such. Le s look a he example from he previous page. Firs we consruc he graph wih capial and consumpion on is axes. Then we pick ou capial, and se he law of moion for «k equal o zero. This solves easily for c as a funcion of k, and we skech in hese poins a which he capial sock is unchanging. [Elsewhere, we do he same hing for he law of moion for «c, and we idenify he locus along which consumpion c ( ) k«α c = Ak δ + n k is unchanging.] Now we need o pick ou a poin above he k-consan locus and one below he locus o ell us he direcion of changes in he capial sock for each he k Summer 00 mah class noes, page 96

c wo regions. I m going o pick one jus off he c axis and one jus off he k axis, since ha ll allow me o essenially se one of he variables equal o zero in each case. Evaluaing he funcion: «k Ak α n k c = ( δ + ) a c = c > 0 and k = ε 0, I find ha k«c = <0. Everywhere above he locus, capial sock is decreasing. I pencil in some lile arrows going o he lef. Then I evaluaion he funcion a c = ε 0 and k = k. I find ha «α k Ak n k = δ +, ( ) k which looks like i should be posiive, as long as k isn oo big (presumable, oo big is where i inersecs he k-axis a second ime. I skech in my arrows. In he oher picure, I do a similar hing for consumpion. Finally, I combine he wo picures, giving me somehing ha looks like his: c «k «k «c And here you have idenified he direcions of change in he variables, so you can also race ou pahs for how hey follow, from a paricular saring poin. Knowing Albania s curren saes, we migh pu hem in he boom lef region, and hen predic ha hey follow a pah of growing consumpion and capial sock for a while bu hen hey cross over he op of he k-seady locus, and consumpion coninues o grow, bu no k he capial sock (maybe hey are selling off heir asses o buy bigger and beer cars). This is he usefulness of phase diagrams. All you need is a pair of differenial equaions o skech hem. Where hose differenial equaions come from is a differen maer. I don ge o each you guys abou opimal conrol heory, so you ll jus have o wai for macro. Summer 00 mah class noes, page 97