Linear Dynamic Models

Similar documents
Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

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

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

Two Coupled Oscillators / Normal Modes

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

Let us start with a two dimensional case. We consider a vector ( x,

BU Macro BU Macro Fall 2008, Lecture 4

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

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

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

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

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

Math 334 Fall 2011 Homework 11 Solutions

10. State Space Methods

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

LAPLACE TRANSFORM AND TRANSFER FUNCTION

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

Theory of! Partial Differential Equations!

Chapter 2. First Order Scalar Equations

Lecture Notes 2. The Hilbert Space Approach to Time Series

Theory of! Partial Differential Equations-I!

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

t + 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 II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y

The expectation value of the field operator.

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

Chapter 6. Systems of First Order Linear Differential Equations

Second Order Linear Differential Equations

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

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

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

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

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.

After the completion of this section the student. Theory of Linear Systems of ODEs. Autonomous Systems. Review Questions and Exercises

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

The Brock-Mirman Stochastic Growth Model

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

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

Vehicle Arrival Models : Headway

6. 6 v ; degree = 7; leading coefficient = 6; 7. The expression has 3 terms; t p no; subtracting x from 3x ( 3x x 2x)

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

System of Linear Differential Equations

Announcements: Warm-up Exercise:

MATH 2050 Assignment 9 Winter Do not need to hand in. 1. Find the determinant by reducing to triangular form for the following matrices.

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

Introduction to Physical Oceanography Homework 5 - Solutions

T. J. HOLMES AND T. J. KEHOE INTERNATIONAL TRADE AND PAYMENTS THEORY FALL 2011 EXAMINATION

Instructor: Barry McQuarrie Page 1 of 5

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

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

Some Basic Information about M-S-D Systems

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

INSTANTANEOUS VELOCITY

INDEX. Transient analysis 1 Initial Conditions 1

Radical Expressions. Terminology: A radical will have the following; a radical sign, a radicand, and an index.

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

Let ( α, β be the eigenvector associated with the eigenvalue λ i

13.3 Term structure models

Introduction to DSGE modelling. Nicola Viegi. University of Pretoria

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

5. Stochastic processes (1)

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Distance Between Two Ellipses in 3D

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

KINEMATICS IN ONE DIMENSION

SOLUTIONS TO ECE 3084

Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.

Math Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.

Section 4.4 Logarithmic Properties

ME 391 Mechanical Engineering Analysis

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

1 Differential Equation Investigations using Customizable

Fishing limits and the Logistic Equation. 1

14 Autoregressive Moving Average Models

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

An Introduction to Malliavin calculus and its applications

STATE-SPACE MODELLING. A mass balance across the tank gives:

Properties of Autocorrelated Processes Economics 30331

Predator - Prey Model Trajectories and the nonlinear conservation law

Chapter 3 Boundary Value Problem

Complete solutions to Exercise 14(b) 1. Very similar to EXAMPLE 4. We have same characteristic equation:

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

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

The Contradiction within Equations of Motion with Constant Acceleration

Solutions Problem Set 3 Macro II (14.452)

Ordinary Differential Equations

Application 5.4 Defective Eigenvalues and Generalized Eigenvectors

Exercises: Similarity Transformation

A Dynamic Model of Economic Fluctuations

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand

T 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

Robert Kollmann. 6 September 2017

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Transcription:

Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference equaions. Analsis of he ssem, focusing on he Samuelson model

. Placing he Samuelson model in sae space form Consumpion funcion c Invesmen funcion = α i = β ( c c ) Implicaion = α( + β ) αβ + g 3 An approach o analzing ssems of difference equaions We are going o use a sandard framework for suding dnamic ssems. I is comprised of wo groups of equaions Equaions describing how observables relae o saes. Each of hese akes he form of a linear relaionship, so ha we have a marix ssem of he form (if we absrac from consan and rend erms) 4

Y S S S j = π j + π j +... + π jn n Y Y π π π n S Y π π π n S = =...... Ym πm πm πmn Sn 5 And he saes evolve as S = M S +... M S + G e +... G e j j, jm m, j jn n S S M M M, S S M M M m, = =...... S m Mm Mm Mmm Sm, S m + Ge SGZ [ more Macro 008: Lecure on 5A his par laer] 6

General framework Can be wrien in simple manner Y =Π S S = MS + Ge Is eas o sud and undersand d Can be used for man of differen problems Highlighs imporan consideraions 7 Solving his ssem recursivel Recursive subsiuion from iniial S, given pah of e s S = MS + Ge 0 S = MS + Ge = M[ MS + Ge ] + Ge 0 = + + M S0 MGe Ge 8

coninuing S = M S0 + M Ge + M Ge +... + MGe + Ge [ check : S = M S0 + MGe+ Ge] 9 So, where we are a depends on The iniial sae S 0 The series of e s ha ake place along he pah beween and How imporan hese e s are for S (G) Various powers of M [ver similar o firs order difference equaion, excep in vecors ] 0

Samuelson s s model (wihou g): sae evoluion = α( + β) αβ S = α( + β) αβ = = 0 MS Samuelson s s model (wihou g): oher variables 0 Y = c = 0 α =ΠS i α i = c = α

Samuelson s model (wih g) Assume g = g + e Implies S = g 3 S Samuelson s model (wih g): saes = g α ( + β ) αβ = 0 0 0 + e 0 0 g = MS + Ge 4

Samuelson s s model (wih g): oher variables Y g = S = c g i g 0 0 0 0 = c 0 α 0 g i α SGZ Macro 008: Lecure 5A 5. Analzing he ssem Qualiaive dnamics of ssem (as in Samuelson s regions) are governed b M. Characerisic roos of M are he soluions o he polnomial Iz-M =0 0, whereiisan ideni marix wih same dimension as M and X means deerminan of X. 6

Samuelson s model (wihou g) S = α ( + β) αβ = = 0 MS 7 Finding roos in Samuelson s model Calculaing Iz-M, we find he same resul as in Samuelson s foonoe #: roos are soluion o quadraic equaion 0 ( + α) β αβ 0 = 0 z 0 z ( + α) β αβ = z = z α( + β) z+ αβ 8

This is a quadraic equaion, so There are wo roos There can eiher be wo real roos or wo complex roos The roos can be of varie of sizes in erms of absolue value: Less han one (sable); Equal o one (uni, borderline); Greaer han one (unsable) Higher order polnomials also have roos wih las wo properies (real or complex, small or large) 9 Complex roos Will no be a special focus of our discussion, alhough he were a ke par of Samuelson s analsis (one definiion of business ccles ) Lead o oscillaor naure of dnamic responses (as in Samuelson s able ) Can eiher be sable (damped) or unsable (explosive) 0

3. Sochasic Ssems such as vecor auoregressions Suppose ha we rea he e in he sae ssem above as a series of unpredicable random variables, so ha our expecaion (forecas) of e +j given informaion a is zero. We wrie his as E e +j =0 Then, i is eas o produce a forecas (compue a condiional expecaion) And his forecas akes a simple form In words, he forecas depends on where we sar (iniial S) and he dnamic forces capured b M. ES j + j = M S + E M + j { Ge+... + MGe+ j + Ge+ j} j = M S

Forecasing oher variables We know ha Y is relaed o S b a simple expression, so we ge forecass ver easil Y =ΠS + j + j E Y =Π ES =ΠM S j + j + j 3 Srucure above absracs from consan erms (normal values) bu hese are eas o add, using eiher of wo mahemaicall equivalen approaches APPROACH : "previous ssem is deviaions from normal values" Y Y =Π( S S) S S = M( S S) + Ge APPROACH : "add consan erms o ever equaion, hen work ou normal values" Y A S S B MS = +Π = + + Ge S = ( I M) B; Y= A+ΠS 4

Puing he ssem o work Empirical models (like VARs) Linear raional expecaions models Relaed heoreical device Soluion of LRE model occurs in sae space form 5 Appendix. Analical framework for dnamic linear ssem. Eigenvecor-Eigenvalue decomposiion of M 3. Link o naure of soluion o second order difference equaion (eg Samuelson model) 6

. Analical quesion Recursive soluion akes form S = M S0 + M Ge + M Ge +... + MGe + Ge Wha is behavior of M raised o he j h power? 7 Ke properies from linear algebra An marix M which has disinc (no repeaed) roos can be decomposed ino M=PμP -, where P is a marix of eigenvecors (characerisic vecors) and μ is a diagonal marix wih he eigenvalues on he diagonal and P - is he inverse of P so ha PP - =I. Furher, M =(PμP - )(PμP - ) = (Pμ P - ) and more generall M j = (Pμ j P - ). Since μ is diagonal, μ j is a diagonal marix wih he eigenvalues raised o he jh power on he diagonal. 8

Samuelson model w/o g ( ) m m M α + β αβ = = 0 0 9 P Eigenvecors and Eigenvalues μ μ μ - = P = μ μ μ μ 0 μ = 0 μ μ = [ m+ m + 4 m] μ = [ m m + 4 m] 30

Powers of M j j M = ( Pμ P ) j + j + j + j + μ μ μμ + μμ = j j j j μ μ μ μ μμ + μμ = [ 0 ] j + j M S μ μ μ μ + μ μ = + j+ j+ j+ j+ μ μ μ μ 3 Soluion Hence, he soluion of he difference equaion is of he form described in he main lecure, = θ μ + θ μ θ = ( μ μ μ )/( μ μ ) 0 θ = ( μ + μ μ )/( μ μ ) 0 3