Notes on the stability of dynamic systems and the use of Eigen Values.

Similar documents
. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

Solution in semi infinite diffusion couples (error function analysis)

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

Chapter Lagrangian Interpolation

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

( ) () we define the interaction representation by the unitary transformation () = ()

2/20/2013. EE 101 Midterm 2 Review

Comb Filters. Comb Filters

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

Density Matrix Description of NMR BCMB/CHEM 8190

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

Let s treat the problem of the response of a system to an applied external force. Again,

Density Matrix Description of NMR BCMB/CHEM 8190

Advanced time-series analysis (University of Lund, Economic History Department)

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

On One Analytic Method of. Constructing Program Controls

First-order piecewise-linear dynamic circuits

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

Chapter 6: AC Circuits

CHAPTER 10: LINEAR DISCRIMINATION

Response of MDOF systems

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

CS286.2 Lecture 14: Quantum de Finetti Theorems II

Mechanics Physics 151

Mechanics Physics 151

( ) [ ] MAP Decision Rule

Mechanics Physics 151

Control Systems. Mathematical Modeling of Control Systems.

UNIVERSITAT AUTÒNOMA DE BARCELONA MARCH 2017 EXAMINATION

Robust and Accurate Cancer Classification with Gene Expression Profiling

FI 3103 Quantum Physics

Part II CONTINUOUS TIME STOCHASTIC PROCESSES

Lecture VI Regression

Born Oppenheimer Approximation and Beyond

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

TSS = SST + SSE An orthogonal partition of the total SS

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

by Lauren DeDieu Advisor: George Chen

Graduate Macroeconomics 2 Problem set 5. - Solutions

Linear Response Theory: The connection between QFT and experiments

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Consider a system of 2 simultaneous first order linear equations

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

Motion in Two Dimensions

Robustness Experiments with Two Variance Components

Volatility Interpolation

GORDON AND NEWELL QUEUEING NETWORKS AND COPULAS

Normal Random Variable and its discriminant functions

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

Summary: Solving a Homogeneous System of Two Linear First Order Equations in Two Unknowns

Relative controllability of nonlinear systems with delays in control

ln 2 1 ln y x c y C x

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

Chapter 2 Linear dynamic analysis of a structural system

Lecture 6: Learning for Control (Generalised Linear Regression)

COMPUTER SCIENCE 349A SAMPLE EXAM QUESTIONS WITH SOLUTIONS PARTS 1, 2

Chapters 2 Kinematics. Position, Distance, Displacement

PHYS 1443 Section 001 Lecture #4

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Observer Design for Nonlinear Systems using Linear Approximations

Comparison of Differences between Power Means 1

3. OVERVIEW OF NUMERICAL METHODS

FTCS Solution to the Heat Equation

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

Scattering at an Interface: Oblique Incidence

PHYS 705: Classical Mechanics. Canonical Transformation

Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system Y X (2) and write EOM (1) as two first-order Eqs.

GMM parameter estimation. Xiaoye Lu CMPS290c Final Project

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

Lecture 2 M/G/1 queues. M/G/1-queue

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

Track Properities of Normal Chain

OP = OO' + Ut + Vn + Wb. Material We Will Cover Today. Computer Vision Lecture 3. Multi-view Geometry I. Amnon Shashua

Department of Economics University of Toronto

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

CHAPTER 5: MULTIVARIATE METHODS

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

Existence of Time Periodic Solutions for the Ginzburg-Landau Equations. model of superconductivity

Cubic Bezier Homotopy Function for Solving Exponential Equations

Pendulum Dynamics. = Ft tangential direction (2) radial direction (1)

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

II The Z Transform. Topics to be covered. 1. Introduction. 2. The Z transform. 3. Z transforms of elementary functions

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION

Advanced Machine Learning & Perception

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment

A Tour of Modeling Techniques

[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Displacement, Velocity, and Acceleration. (WHERE and WHEN?)

Chapter 5. The linear fixed-effects estimators: matrix creation

Linear Dynamic Models

P R = P 0. The system is shown on the next figure:

Should Exact Index Numbers have Standard Errors? Theory and Application to Asian Growth

Transcription:

Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon (994) Tme Seres nalss chaper & Mahemacal Revew ppend. Prepared b rel Vale 3 Le he frs-order lnear non-lnear non-homogeneous ssem of dfference equaons be n general form: z () z f where s a n n mar of me nvaran coeffcens and z and f are n vecor of daed varables. The varable z s pcall he sae vecor of he ssem a me ; f s a vecor of (possbl me-dependen) forcng erms ofen hough of as eogenous. For llusrave purposes we wll work wh a ssem of he form: a c b d () or, a c b d. (3) We can wre he general soluon o an lnear ssem lke n () n he form of a complemenar homogeneous ssem plus he parcular soluon, or: z z z g C P (4) Ths resul (4), allows us o solve for he equlbrum of he ssem n () n wo sages: we sar wh he homogeneous ssem z z, and complee he ask b fndng he parcular soluon o he full ssem. Parcular soluons are eas o fnd when he lernavel he analss can be done n erms of dfferenal equaons, changng noaon (). Recall b Kolmogorov s connu heorem, ha for he dscree { } process and T > f here es posve consans α β [ s ] D s ; s, T E α, β, D s.., hen here ess a connuous verson of, ().

forcng erms are me nvaran because n he sead sae and. Thus, I I where d c b a (5) regroupng erms, ( ) Mar Iden I where (6) and, ( ) I (7) ha s ff ( s nverble,.e. non-sngular. Pluggng (6) no (3): ) (8) ( ) I Regroupng erms, we ge he homogeneous lnear ssem of frs-order dfference equaons: (9) ssume we know he nal condon, hen. () Hence, he ssem converges o he sead sae ff. Ths can be deermned hrough a specral or egenvalue represenaon of. To do hs we need o use Jordan decomposon, whch requres ha mar has lm s n lnearl ndependen egenvecors,.e. dsnc egenvalues (could have some bu never all egenvalues repeaed). The Jordan decomposon of s: () Λ e e

where e s he mar of egenvecors and Λ s he mar of egenvalues, a dagonal mar wh he egenvalues along he prncpal dagonal and zeros elsewhere. For he ssem: λ Λ λ () and because Λ s dagonal, mus be: Λ λ, where λ (3) λ λ The egenvalues are obaned as follows. Whou loss of general drop he me superscrp n (), hen mus be ha: e Λ e defne λ as he vecor of egenvalues and I an (4) n n den mar, hen Λ λ I, and ( λ I ) e λ I e e (5) Thus, ff e s non-sngular, hen mus be ha: ( I ) de λ (6) (6) s known as he characersc equaon from where we oban he egenvalues. Noce, ha for he ssem: a b λ de c d λ ( λ I ) de a λ b de (7) c d λ ( a λ)( d λ) bc ad aλ λ dλ bc λ ( a d ) λ ( ad bc) ( ) T where a d Trace( ) and ( ad bc) de( ) D. Then, he characersc equaon can be re-wren n erms of he T and D as follows: λ Tλ D (8)

Solvng he quadrac equaon for he egenvalues: T T 4D T T 4D λ, and λ (9) λ s he larges egenvalue and λ s he smalles egenvalue. lso, λ λ T and λ λ λ D. Recall ha, Λ, hen we can assess he local sabl properes λ of he planar ssem from he dscrmnan T 4D (a parabola) of he characersc equaon: ) Iff T > 4D hen boh egenvalues are real and dsnc. ) Iff T 4D hen boh egenvalues are equal and mar canno be dagonalzable and we need anoher mehod o fnd he egenvalues. 3) Iff T < 4D, hen boh egenvalues are dsnc and comple conjugaes: T 4D T T 4D T λ, and λ () lso: ) Iff λ <, he modulus of he egenvalues le whn he un crcle (SINK),.e. he sead sae s sable as λ as. B) Iff λ >, he modulus of he egenvalues le ousde he un crcle (SOURCE),.e. he sead sae s unsable as λ ± as. C) Iff s.. λ >, here ess some egenvalue whch modulus les ousde he un crcle (bu no all) (SDDLE), hen he sead sae s neher sable nor unsable. Onl a ver specfc se of nal condons wll ake he ssem o he sead sae along he saddle pah (sable manfold) drven b he sable egenvalues. Fnall facorzng he characersc equaon p, we oban he sragh lnes p() p( ), whch jonl wh he dscrmnan parabola dvdes he plane n each regons for he characerzaon of he sabl properes of he ssem n erms of he egenvalues (roos), or equvalenl from he race (T) and he deermnan (D) of he coeffcen mar (see graph below):

D T < - T > D > Comple Roos D > D > 4 5 8 p() p(-) Comple Roos D < 6 7a 7c - - T (,) & D (, ) T Real Roos Real Roos T < 7b Real Roos T > 3 - Real Roos D < - Real egenvalues zones: Regon : Saddle, Regon : Source, oscllaor dvergence, Regon 3: Saddle, Regon 4: Source, monoone dvergence, Regon 7: Snk, 7a & 7c monoone convergence, 7b oscllaor convergence, Regon 8: Source, monoone dvergence, Comple egenvalues (conjugae pars) zones: Regon 5: Source, oscllaor dvergence, Regon : Snk, oscllaor convergence.

FINL REMRK: For he nonlnear case, we lnearzed he ssem usng Talor s epanson and fndng he Jacoban (J) evaluaed a he sead sae. Then we can proceed wh he prevous analss changng for J.