Métodos para obtener soluciones numéricas de los modelos DSGE

El método de Blanchard y Khan (1980)

Un ejemplo
Un ejemplo en Matlab

Clase 5B: Métodos de solución de modelos

DSGE
Macrodinámica I

Hamilton Galindo

Junio - Agosto
2015

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Contenido

1 Métodos para obtener soluciones numéricas de los modelos DSGE

2 El método de Blanchard y Khan (1980)

3 Un ejemplo

4 Un ejemplo en Matlab

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Fundamentos I

1 Métodos para obtener soluciones numéricas de los modelos DSGE

2 Estos métodos son usados para obtener soluciones estables para cual-
quier sistema de ecuaciones en diferencia estocástica no lineal.
3 En la solución del modelo nos tenemos que asegurar que la solución
sea estable; es decir, una solución en la cual las variables per cápita
no explote rápidamente. Para ello hay que imponer condiciones de
estabilidad sobre la solución.
4 Solución numérica: es un conjunto de series de tiempo, una para ca-
da variable relevante en el modelo económico, que en cada periodo
satisface todas las condiciones en el modelo.
5 Simulación: es un procedimiento por el cual una solución numérica
es encontrada para cada realización del vector estocástico del choque
exógeno que afecta a la economı́a.
6 Tres métodos para resolver el sistema log-lineal (encontrar la solución
numérica):

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Fundamentos II

Coeficientes indeterminados (Uhlig, xxx)

Blanchar y Kahn ()
Método de descomposición eigenvalore-eigenvectores (C.A. Sims,
xxxx)
Nota: estos métodos son para sistemas lineales (o log-lineales). Los
métodos para sistemas no lineales son dos: [1] Método de expectativas
parametrizadas (den Hann-Marcet): usan las condiciones de primer or-
den, pero apróxima las expectativas condicionales (que aparecen en la
condiciones de primer orden) por polinomios exponeniales y [2] Clase
de nétodos de proyección: que parametrizan las reglas de decisión o
las ecuaciones de control como funciones polinomiales de las variables
de estado.
7 La principal diferencia de estos métodos es el camino por el cual las
condiciones de estabilidad son impuestas sobre la solución numérica.
8

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Método de Blanchard y Khan (1980)

1 Este método está basado en la idea de resolver el sistema lineal por

medio de encontrar la senda estable del sistema no lineal.
2

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Forma estado-espacio I

[1] Forma estado-espacio generalizado

La forma estado-espacio generalizada para sistemas de ecuaciones no li-

neales, dinámicas y estocásticas con expectativas es:

A0 Et Xt+1 = A1 Xt + B0 vt+1 (1)

En la literatura existe varios métodos para resolver este tipo de mo-

delos (ver DeJong y Dave, XXXX):
El método de Blanchard y Kahn (1980) usa la descomposición de
Jordan cuando la matriz A0 es invertible. Cuando no es invertible usa
la descomposición QZ.

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Forma estado-espacio II

[2] Forma estado-espacio alternativa

A0 Et Xt+1 = A1 Xt + B0 vt+1
Et Xt+1 = A−1 A1 Xt + A−1 B0 vt+1
| 0{z } | 0{z }
=A =B
Et Xt+1 = AXt + Bvt+1 (2)

En la ecuación [2] la variable vt+1 puede ser choques iid con: media
igual a cero (E (vt ) = 0) y autocorrelación nula. Alternativamente,
vt+1 puede comportarse como un proceso AR(1), la cual depende de
choques exógenos iid.

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Forma estado-espacio III

En este sistema de ecuaciones se puede definir dos tipos de variables:

[1] wt es el vector de variables backward looking (predeterminadas o
variables de estado), y [2] yt es el vector de variables forward looking
(variables de control):  
w
Xt = t
yt
Introduciendo lo anterior en la ecuación [2] se tiene:
   
wt+1 w
= A t + Bvt+1 (3)
Et yt+1 yt

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Transformando el sistema de ecuaciones I

[1] Descomposición de Jordan de A

   
wt+1 w
= A t + Bvt+1 (4)
Et yt+1 yt

La matriz A puede ser expresada (por la descomposición de Jordan)

de la siguiente manera:
A = PΛP −1
Donde P son los eigenvectores y Λ es la matrix diagonal de eigenva-
lores.

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Transformando el sistema de ecuaciones II

Al introducir la descomposición de A en la ecuación [4] se tiene:

   
wt+1 w
= PΛP −1 t + Bvt+1 (5)
Et yt+1 yt
Ordenando los términos:
   
wt+1 w
P −1 = ΛP −1 t + P −1 Bvt+1 (6)
Et yt+1 yt
Donde: P −1 B = R

[2] Partición del modelo

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Transformando el sistema de ecuaciones III

La matriz Λ, que representa los valores propios, se puede expresar de

la siguiente manera:  
Λ
Λ= 1 (7)
Λ2
En la cual Λ1 es el eigenvalor estable y Λ2 es el inestable.
En base a esta partición la matriz de eigevectores se particiona
también:  
−1 P11 P12
P =
P21 P22
Asimismo se hace con matrix R:
 
R1
R=
R2

[3] Sistema de ecuaciones transformado

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Transformando el sistema de ecuaciones IV

Al introducir las particiones descritas en la ecuación [6] se tiene:

        
P11 P12 wt+1 Λ P11 P12 wt R
= 1 + 1 vt+1
P21 P22 Et yt+1 Λ2 P21 P22 yt R2
(8)
Se define dos nuevas variables: yet y w
et
    
P11 P12 wt w
= t (9)
e
P21 P22 Et yt yet
Lo anterior implica que:

w
et = P11 wt + P12 yt (10)
yet = P21 wt + P22 yt (11)

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Transformando el sistema de ecuaciones V

[4] Desacoplamiento de las ecuaciones

De la ecuación [9] se podrı́a obtener dos ecuaciones; es decir, el

sistema se puede desacoplar:
    
P11 P12 wt w
= t (12)
e
P21 P22 Et yt yet
Esta ecuación implica dos ecuaciones:

Ecuación estable
w
et+1 = Λ1 w
et + R1 vt+1 (13)
Ecuación inestable
yet+1 = Λ2 w
et + R2 vt+1 (14)

La ventaja de este desacoplamiento es que cada ecuación puede ser

resuelta por separado.
Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE
 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Solución I

[1] Estrategia de solución

Resolver la ecuación inestable transformada.

Resolver la ecuación estable transformada.
Trasladar la solución anterior al problema original.

[2] Solución de la ecuación inestable

La ecuación inestable es:

yet+1 = Λ2 w
et + R2 vt+1 (15)

Resolviendo esta ecuación hacia adelante (t + j):

Et yet+j = (Λ2 )j yet (16)

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Solución II

Como kΛ2 k > 1, la única solución estable es yet = 0, ∀t. Entonces se

tiene que:
yet = P21 wt + P22 yt = 0
y por tanto:
−1
yt = −P22 P21 wt (17)
Esta ecuación expresa las variables de control (forward looking) en
función de las variables de estado (predeterminadas). Por tanto, esta
ecuación representa la función de polı́tica.

[3] Solución de la ecuación estable

La ecuación estable es:

w
et+1 = Λ1 w
et + R1 vt+1 (18)

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Solución III

Resolviendo esta ecuación hacia adelante (t + j):

et+j = (Λ1 )j w
Et w et (19)

Como kΛ1 k < 1, no hay problemas de inestabilidad. Entonces al

resolver simultaneamente la ecuación [17] que se obtiene del
−1
componente inestable (yt = −P22 P21 wt ) y la definición de w
et
(w
et = P11 wt + P12 yt ) se tiene que:

−1
et = (P11 − P12 P22
w P21 )wt (20)
Reemplazando esta última expresión en la ecuación estable [17]:

w
et+1 = Λ1 w
et + R1 vt+1
Considerando
−1 −1
et = (P11 − P12 P22
w et+1 = (P11 − P12 P22
P21 )wt w P21 )wt+1

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Solución IV

Se tiene que:

−1 −1
wt+1 = (P11 − P12 P22 P21 )−1 Λ1 (P11 − P12 P22 P21 )wt (21)
−1
+(P11 − P12 P22 P21 )−1 R1 vt+1

Esta ecuación expresa la dinámica de las variables de estado.

[4] Solución total

La solución total del sistema está descrito por la ecuación [21] y [17]:
−1
yt = −P22 P21 wt
−1 −1
wt+1 = (P11 − P12 P22 P21 )−1 Λ1 (P11 − P12 P22 P21 )wt
−1
+(P11 − P12 P22 P21 )−1 R1 vt+1

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Solución V

Todas las variables están en función de las variables de estado (predetermi-

nadas). Cabe mencionar que este sistema tiene una estructura recursiva.

[5] Formulación recursiva

Para el cálculo de IRF y momentos se sigue los siguientes pasos:

1 Empezar del valor de estado estacionario: w0 = 0 (recordar que las
variables son log-lineales).
2 Simular choques {vt } de la distribución normal.
3 Simular {wt } de {vt } recursivamente usando la ecuación de la
dinámica de las variables de estado:
−1 −1
wt+1 = (P11 − P12 P22 P21 )−1 Λ1 (P11 − P12 P22 P21 )wt
−1
+(P11 − P12 P22 P21 )−1 R1 vt+1

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Solución VI

4 Calcular {yt } de {wt } usando la función de polı́tica:

−1
yt = −P22 P21 wt

5 Finalmente, calcular IRF y momentos.

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Un ejemplo I
Modelo de Brock y Mirmam (xxxx)

el modelo está en la pag 203

[1] Condiciones de primer orden

[2] Sistema log-lineal

pag 210

[3] Solución del sistema: método de BK

pag 212

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Un ejemplo en Matlab
Modelo de Brock y Mirmam (xxxx)

ver el programa Blanchard Kahn.m

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE

 Métodos para obtener soluciones numéricas de los modelos DSGE
El método de Blanchard y Khan (1980)
Un ejemplo
Un ejemplo en Matlab

Qué método usa Dynare

El método de perturbación

Hamilton Galindo Clase 5B: Métodos de solución de modelos DSGE