FACULTAD DE CIENCIAS ECONÓMICAS

“Año de la Lucha contra la Corrupción y la Impunidad”

CURSO: ANÁLISIS MATEMÁTICO PARA ECONOMISTAS III

CICLO: 3er
AULA/TURNO: 54-E Tarde
DOCENTE: GARCÍA BONILLA, Quicio Alberto
INTEGRANTES:

 ALCALÁ USUCACHI, María Isabel

 CASTRO CABRERA, Angie Betsabé
 VARA CARHUAMACA, Dana Samira
 VILLAR MACHACA, Lisbeth Marisol
 CONCHA ALEJO, Nadia Jimena
 SALVADOR ANTÓN, Mario elíseo

Lima, agosto 2019

 INDICE

I. INTRODUCCIÓN

1. ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO ORDEN CON

COEFICIENTES CONSTANTES Y TÉRMINO CONSTANTE

a) La solución particular

b) La función complementaria

c) La convergencia de la trayectoria de tiempo

2. MODELO DE INTERACCIÓN DE MULTIPLICADOR CON ACELERADOR DE

SAMUELSON

a) Marco de referencia

b) La solución

c) Convergencia contra divergencia

d) Un resumen gráfico

3. LA INFLACIÓN Y EL DESEMPLEO EN TIEMPO DISCRETO

a) El modelo

b) La ecuación en diferencias en p

c) La trayectoria de tiempo de p

d) El análisis de U

e) La relación de Phillips de largo plazo

4. GENERALIZACIONES A ECUACIONES CON TÉRMINOS VARIABLES Y DE

ORDEN SUPERIOR

a) El término variable con forma de Cm^t

b) El término variable con forma de Ct^n

 c) Ecuaciones en diferencias lineales de orden superior

d) La convergencia y el teorema de Schur

5. MODELO DE INVENTARIO DE METSLER

II. CONCLUSIONES

III. BIBLIOGRAFÍA
 I. INTRODUCCIÓN

En el presente trabajo consideraremos sistemas de ecuaciones en los que cada variable está

indexada en el tiempo t = 0, 1, 2,..., y variables correspondientes a tiempos distintos están

relacionadas de una manera no trivial. Tales sistemas se denominan sistemas de ecuaciones en

diferencias y son muy útiles para describir sistemas dinámicos en tiempo discreto

El estudio de modelos dinámicos en economía es importante dado que permite eliminar la hipótesis

(estática) de que el proceso de ajuste es instantáneo e inevitablemente da lugar a un equilibrio. En

un contexto dinámico, esta propiedad de estabilidad tiene que ser comprobada y no puede ser

asumida a priori.

Veremos cómo algunos modelos económicos incluyen ecuaciones en diferencias que relacionan

entre sí a 𝑃𝑡 con 𝑃𝑡−1 . Como el valor de P en un periodo puede determinar en forma única al valor

de P en el siguiente, la trayectoria de tiempo de P se hace totalmente determinada una vez que se

especifica un valor inicial de 𝑃0 .

Sin embargo, suele ocurrir que el valor de una variable económica para el periodo t (digamos 𝑦𝑡 )

depende no sólo de 𝑦𝑡−1 sino también de 𝑦𝑡−2 . Esta situación va a dar lugar a una ecuación en

diferencias de segundo orden.

1. ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO ORDEN CON

COEFICIENTES CONSTANTES Y TÉRMINO CONSTANTE

Las ecuaciones en diferencias de segundo orden suelen adoptar la siguiente forma:

𝑦𝑡+2 + 𝑎1 𝑦𝑡+1 + 𝑎2 𝑦𝑡 = 𝐶

Reconocemos que esta ecuación es lineal, no homogénea, con coeficientes contantes y término

constante.

a) La Solución Particular

Podemos esperar que la solución de la ecuación anterior tenga dos componentes: una solución

particular 𝑦𝑝 , que representa el nivel de equilibrio intertemporal de 𝑦, y una función

complementaria 𝑦𝑐 que específica para cada periodo la desviación respecto al equilibrio. La

solución particular, que se define como cualquier solución de la ecuación completa, puede algunas

veces encontrarse simplemente al intentar una solución de la forma 𝑦𝑡 = 𝑘. Al sustituir este valor

constante de 𝑦 en la ecuación lineal planteada inicialmente, obtenemos:

𝑦𝑡+2 + 𝑎1 𝑦𝑡+1 + 𝑎2 𝑦𝑡 = 𝐶

𝑘 + 𝑎1 𝑘 + 𝑎2 𝑘 = 𝐶

𝑘(1 + 𝑎1 + 𝑎2 ) = 𝐶

𝐶
𝑘=
1 + 𝑎1 + 𝑎2

Entonces siempre que (1 + 𝑎1 + 𝑎2 ) ≠ 0, la integral particular será:

𝐶
𝑦𝑝 = 𝑘 =
1 + 𝑎1 + 𝑎2

En caso de 𝑎1 + 𝑎2 ≠ −1.
En caso de 𝑎1 + 𝑎2 = −1. La solución de prueba 𝑦𝑡 = 𝑘 fallaría y podemos intentar en su lugar

𝑦𝑡 = 𝑘𝑡, por lo tanto al sustituir tendríamos:

𝑦𝑡+2 + 𝑎1 𝑦𝑡+1 + 𝑎2 𝑦𝑡 = 𝐶

𝑘(𝑡 + 2) + 𝑎1 𝑘(𝑡 + 1) + 𝑎2 𝑘𝑡 = 𝐶

𝑘(𝑡(1 + 𝑎1 + 𝑎2 ) + 𝑎1 + 2) = 𝐶

𝐶
𝑘=
𝑡(1 + 𝑎1 + 𝑎2 ) + 𝑎1 + 2

𝐶
𝑘=
𝑎1 + 2

𝐶
𝑦𝑝 = 𝑘𝑡 = 𝑡
𝑎1 + 2

En caso de 𝑎1 + 𝑎2 = −1, 𝑎1 ≠ −2.

b) La Función Complementaria

Comenzaremos con una ecuación reducida:

𝑦𝑡+2 + 𝑎1 𝑦𝑡+1 + 𝑎2 𝑦𝑡 = 0

Por ecuaciones diferenciales de primer orden intentaremos una solución de la forma:

𝑦𝑡 = 𝐴𝑏 𝑡

𝑦𝑡+1 = 𝐴𝑏 𝑡+1

𝐴𝑏 𝑡+2 + 𝑎1 𝐴𝑏 𝑡+1 + 𝑎2 𝐴𝑏 𝑡 = 0

𝐴𝑏 𝑡 (𝑏 2 + 𝑎1 𝑏 + 𝑎2 ) = 0
 𝑏 2 + 𝑎1 𝑏 + 𝑎2 = 0

Esta ecuación cuadrática posee dos raíces características:

−𝑎1 ± √𝑎12 − 4𝑎2

𝑏1 , 𝑏2 =
2

Cada una de las cuales es aceptable en la solución 𝐴𝑏 𝑡 . De hecho, tanto 𝑏1 como 𝑏2 deben aparecer

en la solución general de la ecuación en diferencia homogénea porque, al igual que en el caso de

las ecuaciones diferenciales, esta solución general debe consistir de dos partes linealmente

independientes, cada una con su propia constante arbitraria multiplicativa.

Podemos encontrar tres situaciones posibles respecto a las raíces características, dependiendo

de la expresión con raíz cuadrada.

 Caso 1 (Raíces Reales Diferentes)

Cuando 𝑎12 > 4𝑎2 , la raíz cuadrada es un número real, y 𝑏1 𝑦 𝑏2 son reales y diferentes. En ese

caso, 𝑏1𝑡 𝑦 𝑏2𝑡 son linealmente independientes y la función complementaria puede escribirse

simplemente como una combinación lineal de estas expresiones; es decir:

𝑦𝑐 = 𝐴1 𝑏1𝑡 + 𝐴2 𝑏2𝑡

 Caso 2 (Raíces Reales Repetidas)

Cuando 𝑎12 = 4𝑎2 se anula la raíz cuadrada y las raíces características se repiten:

𝑎1
𝑏(= 𝑏1 = 𝑏2 ) = −
2

La función complementaria se puede expresar en un solo termino:

𝑦𝑐 = 𝐴1 𝑏1𝑡 + 𝐴2 𝑏2𝑡 = (𝐴1 + 𝐴2 )𝑏 𝑡 ≡ 𝐴3 𝑏 𝑡

Para suministrar la constante que nos falta, debemos multiplicar 𝑏 𝑡 por la variable 𝑡. Por lo tanto,

el nuevo termino tendrá la forma de:

𝐴4 𝑡𝑏 𝑡

Finalmente, la función complementaria para el caso de raíces repetidas nos queda de la siguiente

forma:

𝑦𝑐 = 𝐴3 𝑏 𝑡 + 𝐴4 𝑡𝑏 𝑡

 Caso 3 (Raíces Complejas)

Cuando 𝑎12 < 4𝑎2 ,las raíces características son complejas conjugadas, están en la forma:

𝑎1 √4𝑎2 − 𝑎12
𝑏1 , 𝑏2 = − ±
2 2

√4𝑎2 −𝑎12
𝑎1
Donde – =ℎy = 𝑣, entonces la función complementaria tomara la forma de:
2 2

𝑦𝑐 = 𝐴1 𝑏1𝑡 + 𝐴2 𝑏2𝑡

𝑦𝑐 = 𝐴1 (ℎ + 𝑣)𝑡 + 𝐴2 (ℎ − 𝑣)𝑡

Para facilitar la solución de la Función Complementaria, usaremos el Teorema de Moivre que nos

permite transformar la Función Complementaria en términos trigonométricos, tomara la siguiente

forma:

(ℎ ± 𝑣)𝑡 = 𝑅 𝑡 (𝐶𝑜𝑠𝜃𝑡 ± 𝑆𝑒𝑛𝜃𝑡)

Donde el valor de R siempre será considerado positivo y tomará la siguiente forma:

 𝑎12 + 4𝑎2 − 𝑎12
𝑅 = √ℎ 2 + 𝑣 2 = √ = √𝑎2
4

Y 𝜃 es la medida en radianes del ángulo en el intervalo [0;2π), que satisface las condiciones:

ℎ −𝑎1
𝐶𝑜𝑠 𝜃 = =
𝑣 2√𝑎2

𝑣 𝑎12
𝑆𝑒𝑛 𝜃 = = √1 −
𝑅 4𝑎2

La función complementaria finalmente toma la forma de:

𝑦𝑐 = 𝐴1 𝑅 𝑡 (𝐶𝑜𝑠𝜃𝑡 + 𝑆𝑒𝑛𝜃𝑡) + 𝐴2 𝑅 𝑡 (𝐶𝑜𝑠𝜃𝑡 − 𝑆𝑒𝑛𝜃𝑡)

𝑦𝑐 = 𝑅 𝑡 [(𝐴1 + 𝐴2 )𝐶𝑜𝑠𝜃𝑡 + (𝐴1 − 𝐴2 )𝑆𝑒𝑛𝜃 ]

c) La convergencia de la Trayectoria del Tiempo

Como en el caso de las ecuaciones diferenciales de primer orden, la convergencia de la

trayectoria de tiempo 𝑦𝑡 depende solo de si 𝑦𝑐 tiene a cero cuando 𝑡 → ∞.

A diferencia de las ecuaciones diferenciales de primer orden, tendremos que considerar dos

raíces características en lugar de una. Tendremos los siguientes casos de las raíces reales

diferentes:

 Si |𝑏1 | > 1 y |𝑏2 | > 1 , entonces ambos términos componentes en la función

complementaria 𝐴1 𝑏1𝑡 y 𝐴2 𝑏2𝑡 serán explosivos y, por lo tanto, 𝑦𝑐 debe ser divergente.

 Si |𝑏1 | < 1 y |𝑏2 | < 1 , entonces ambos términos de 𝑦𝑐 van a convergir a cero a medida

que 𝑡 aumenta indefinidamente, así como también lo hará 𝑦𝑐 .

  Si |𝑏1 | > 1 y |𝑏2 | < 1 , el termino 𝐴2 𝑏2𝑡 tiende que disminuir, mientras que el otro termino

tiende a desviarse cada vez mas de cero, la trayectoria será divergente.

A la raíz con el valor absoluto más alto le llamaremos raíz dominante.

Una trayectoria de tiempo será convergente si y solo si la raíz dominante es menor que 1 en el valor

absoluto.

2. MODELO DE INTERACCIÓN DE MULTIPLICADOR CON ACELERADOR DE

SAMUELSON

Como una ilustración del uso de las ecuaciones en diferencias de segundo orden en economía,

citemos un trabajo clásico del profesor Paul Samuelson, el primer economista que gano el premio

Novel. Nos referimos a su modelo clásico de interacción, que busca explorar el proceso dinámico

de la determinación del ingreso cuando el principio de aceleración está en operación junto al

multiplicador Keynesiano. Entre otras cosas, ese modelo sirve para demostrar que la sola

interacción del multiplicador y del acelerador tiene la capacidad de generar fluctuaciones cíclicas

en forma endógena.

a) El marco de referencia

Suponga que el ingreso nacional 𝑌𝑡 está compuesto de tres componentes de corriente de gasto: el

consumo 𝐶𝑡 , la inversión 𝐼𝑡 y el gasto del gobierno 𝐺𝑡 . El consumo se concibe como una función

no del ingreso presente sino del ingreso del periodo anterior 𝑌𝑡−1 ; por simplicidad, se supone que

𝐶𝑡 es estrictamente proporcional a 𝑌𝑡−1 . La inversión, que es de la variedad “inducida”, es una

función de la tendencia prevaleciente del gasto del consumidor. Por supuesto que es a través de

esta inversión inducida como el principio de aceleración entra en el modelo. Específicamente,

vamos a suponer que 𝐼𝑡 conserva una relación fija con el incremento del consumo ∆𝐶𝑡−1 ≡ 𝐶𝑡 −
𝐶𝑡−1 . Por otro lado, el tercer componente 𝐺𝑡 se considera como exógeno; de hecho, vamos a

suponerlo como constante y vamos a denotarlo simplemente como 𝐺0 .

Estas hipótesis pueden traducirse en el siguiente conjunto de ecuaciones:

𝑌𝑡= 𝐶𝑡 + 𝐼𝑡 + 𝐺0

𝐶𝑡 = 𝑦𝑌𝑡−1 (0 < 𝑦 < 1) (18.13)

𝐼𝑡 = 𝛼(𝐶𝑡 − 𝐶𝑡−1 ) (𝛼 > 0)

Donde 𝑦 (la letra griega gamma) representa la propensión marginal al consumo y 𝛼 representa el

acelerador (abreviatura de coeficiente de aceleración). Observe que si se cancela la inversión

inducida en el modelo nos queda una ecuación en diferencias de primer orden que abarca el proceso

del multiplicador dinámico. Sin embargo, si se incluye la inversión inducida, tenemos una ecuación

en diferencias de segundo orden que ilustra la interacción del multiplicador y el acelerador.

En virtud de la segunda ecuación, podemos expresar 𝐼𝑡 en términos de ingreso como sigue:

𝐼𝑡 = 𝛼(𝑦𝑌𝑡−1 − 𝑦𝑌𝑡−2 ) = 𝛼𝑦(𝑌𝑡−1 − 𝑌𝑡−2 )

Al sustituir esto y la ecuación de 𝐶𝑡 en la primera ecuación en (18.13) y reordenando, el modelo

puede condensare en la ecuación individual

𝑌𝑡 − 𝑦(1 + 𝛼)𝑌𝑡−1 + 𝛼𝑦𝑌𝑡−2 = 𝐺0

O en forma equivalente (después de desplazar los subíndices hacia adelante por dos periodos),

𝑌𝑡+2 − 𝑦(1 + 𝛼)𝑌𝑡+1 + 𝛼𝑦𝑌𝑡 = 𝐺0 (18.14)

Dado que ésta es una ecuación en diferencia lineal de segundo orden con coeficientes constantes y

término constante, podemos resolverla por el método recién aprendido.

 b) La solución

Como solución particular tenemos, por (18.2),

𝐺0 𝐺0
𝑌𝑝 = =
1 − 𝑦(1 + 𝛼) + 𝛼𝑦 1 − 𝑦

Podemos observar que la expresión 1/(¡ −𝑦) es simplemente el multiplicador que prevalecería en

ausencia de la inversión inducida. Por lo tanto 𝐺0 /(1 − 𝑦) -el elemento exógeno de gasto por el

multiplicador- debe darnos el ingreso de equilibrio “𝑌 ∗ ” en el sentido de que este nivel de ingreso

satisface la condición de equilibrio “ingreso nacional = gasto total”. Sin embargo, siendo la

solución particular del modelo, también nos da el ingreso 𝑌̅ de equilibrio intemporal.

Respecto a la función complementaria, hay tres casos posibles. El caso 1 (𝑎12 > 4𝑎2), en el presente

contexto, se caracteriza por

𝑦 2 (1 + 𝛼)2 > 4𝛼𝑦 𝑜 𝑦(1 + 𝛼)2 > 4𝛼

O sea

4𝛼
𝑦>
(1 + 𝛼)2

En forma similar, para caracterizar los casos 2 y 3, solo necesitamos cambiar el signo > de la última

desigualdad por = 𝑦 <, respectivamente. En la figura 18.2 hemos dibujado la gráfica de la ecuación

𝑦 = 4𝛼/(1 + 𝛼)2 . De acuerdo con la discusión anterior, los pares ( 𝛼, 𝑦) que se localizan

exactamente sobre esta curva pertenecen al caso 2. Por otro lado, los pares (𝛼, 𝑦) situados arriba

de esta curva (que implica valores mas altos que 𝑦) tienen que ver con el caso 1 y los situados

debajo de la curva con el caso 3.

Esta clasificación es gráfica en la figura 18.2, es de interés porque revela claramente las condiciones

bajo las cuales las fluctuaciones cíclicas pueden surgir en forma endógena de la interacción del

multiplicador con el acelerador. Pero esto no nos dice nada acerca de la convergencia o la

divergencia de la trayectoria de tiempo de 𝑌. Por lo tanto, a nosotros nos toca distinguir, para cada

caso, entre los sub casos amortiguado y explosivo. Por supuesto que podríamos tomar el camino

fácil ilustrando simplemente estos sub casos al citar ejemplos numéricos específicos. Pero

intentemos la tarea más gratificante, aun cuando más ardua, de encontrar las condiciones generales

bajo las cuales prevalecen la convergencia y la divergencia.

c) Convergencia contra divergencia

La ecuación en diferencias (18.14) tiene la ecuación característica

𝑏 2 − 𝑦(1 + 𝛼)𝑏 + 𝛼𝑦 = 0

Que suministra las dos raíces

 𝑦(1 + 𝛼) ± √𝑦 2 (1 + 𝛼)2 − 4𝛼𝑦
𝑏1 , 𝑏2=
2

Puesto que la cuestión de la convergencia contra la divergencia depende de los valores de 𝑏1 y 𝑏2 ,

y como 𝑏1 y 𝑏2 a su vez dependen de los valores de los parámetros 𝛼 y 𝑦, las condiciones para

convergencia y divergencia deben ser expresables en términos de los valores de 𝛼 y 𝑦. Para hacer

esto, podemos usar el hecho de que -por (16.2)- las dos raíces características siempre se relacionan

entre sí por las siguientes ecuaciones:

𝑏1 + 𝑏2 = 𝑦(1 + 𝛼) (18.15)

𝑏1 𝑏2 = 𝛼𝑦 (18.15´)

Basándonos en estas dos ecuaciones, podemos observar que

(1 − 𝑏1 )(1 − 𝑏2 ) = 1 − (𝑏1 + 𝑏2 ) + 𝑏1 𝑏2

= 1 − 𝑦(1 + 𝛼) + 𝛼𝑦 = 1 − 𝑦 (18.16)

En vista de la especificación del modelo de que 0 < 𝑦 < 1, se hace necesario imponer sobre las

dos raíces la condición

0 < (1 − 𝑏1 )(1 − 𝑏2 ) < 1 (18.17)

Examinemos ahora la cuestión de la convergencia para el caso 1, en el cual las raíces son reales y

diferentes. Dado que por hipótesis 𝛼 y 𝑦 son ambos positivos, (18.15´) nos dice que 𝑏1 𝑏2 > 0, lo

que implica que 𝑏1 y 𝑏2 poseen el mismo signo algebraico. Aún más, ya que 𝑦(1 + 𝛼) > 0 (18.15)

indica que tanto 𝑏1 como 𝑏2 deben ser positivos. Así, la trayectoria de tiempo 𝑌𝑡 no puede tener

oscilaciones en el caso 1.
Aun cuando los signos 𝑏1 y 𝑏2 se conocen ahora, para el caso 1 existen en realidad cinco

combinaciones posibles de valores de (𝑏1 , 𝑏2 ), cada uno con su propia implicación respecto a los

valores correspondientes de 𝛼 y 𝑦:

𝑖) 0 < 𝑏2 < 𝑏1 < 1 → 0 < 𝑦 < 1; 𝛼𝑦 < 1

𝑖𝑖) 0 < 𝑏2 < 𝑏1 = 1 → 𝑦 = 1

𝑖𝑖𝑖) 0 < 𝑏2 < 1 < 𝑏1 → 𝑦 > 1

𝑖𝑣) 1 = 𝑏2 < 𝑏1 → 𝑦 = 1

𝑣) 1 < 𝑏2 < 𝑏1 → 0 < 𝑦 < 1; 𝛼𝑦 > 1

La posibilidad 𝑖, para la cual 𝑏1 y 𝑏2 son fracciones positivas, satisface debidamente la condición

(18.17) y, por lo tanto, está de acuerdo con la especificación del modelo 0 < 𝑦 < 1. El producto

de las raíces también debe ser una fracción positiva bajo esta posibilidad y mediante (18.15´),

implica que 𝛼𝑦 < 1. En contraste, las tres siguientes posibilidades violan la condición (18.17) y

son valores inadmisibles de 𝑦 (vea el ejercicio 18.2-3). Por lo tanto, deben descartarse. Pero la

posibilidad 𝑣 puede ser aceptable. Con 𝑏1 y 𝑏2 , mayores que uno, (18.17) todavía puede satisfacerse

si (1 − 𝑏1 )(1 − 𝑏2 ) < 1. Pero esta vez tenemos 𝛼𝑦 > 1 (en vez de < 1) de (18.15´). La conclusión

es que existe solo dos sub casos admisibles para el caso 1. El primero –la posibilidad 𝑖—incluye

las raíces fraccionarias 𝑏1 y 𝑏2 ,y por lo tanto suministra una trayectoria de tiempo convergente de

𝑌. El otro sub caso –la posibilidad 𝑣 – se caracteriza por raíces mayores que uno, y por lo tanto

produce una trayectoria de tiempo divergente. Sin embargo, por lo que toca a los valores de 𝛼 y 𝑦,

el asunto de convergencia y la divergencia depende solo de si 𝛼𝑦 < 1 o 𝛼𝑦 > 1. Esta información

se resume en la parte superior a la tabla 18.1, donde el sub caso convergente se rotula como 1𝐶, y

el sub caso divergente como 1𝐷.

El análisis del caso 2, con raíces repetidas, es de naturaleza similar. Las raíces son ahora 𝑏 = 𝑦(1 +

𝛼)/2, con un signo positivo porque 𝛼 y 𝑦 son positivos. Entonces, nuevamente no hay oscilación.

Esta vez podemos clasificar el valor de 𝑏 solo de acuerdo con tres posibilidades:

𝑣𝑖) 0 < 𝑏 < 1 → 𝑦 < 1 ; 𝛼𝑦 < 1

𝑣𝑖𝑖) 𝑏 = 1 → 𝑦 = 1

𝑣𝑖𝑖𝑖) 𝑏 > 1 → 𝑦 < 1 ; 𝛼𝑦 > 1

Para la posibilidad 𝑣𝑖, 𝑏(= 𝑏1 = 𝑏2 )es una fracción positiva; entonces las implicaciones relativas

a 𝛼 y 𝑦 son totalmente idénticas con las de la posibilidad 𝑖 para el caso 1. De manera análoga, la

posibilidad 𝑣𝑖𝑖𝑖, con 𝑏(= 𝑏1 = 𝑏2 )mayor que uno, puede satisfacer a (18.17) solo si 1 < 𝑏 < 2;

si es así, suministra los mismos resultados que la posibilidad 𝑣. Por otro lado, la posibilidad 𝑣𝑖𝑖

viola (18.17) y debe descartarse. Entonces de nuevo hay solo dos sub casos admisibles. El primero

– la posibilidad 𝑣𝑖—suministra una trayectoria de tiempo convergente, mientras que el otro –la

posibilidad 𝑣𝑖𝑖𝑖—da una divergente. En términos de 𝛼 y 𝑦, los sub casos convergente y divergente

se asocian de nuevo, respectivamente, con 𝛼𝑦 < 1 y 𝛼𝑦 > 1. Estos resultados se listan en la parte

media de la tabla 18.1, donde los dos sub casos se rotulan 2𝐶 (convergente) y 2𝐷 (divergente).

TABLA 18.1 Casos y sub casos del modelo de Samuelson

Caso Sub caso Valores de 𝜶 y 𝒚 Trayectoria de tiempo 𝒀𝒕

[Link] reales

diferentes
 𝟒𝜶 1𝐶: 0 < 𝑏2 < 𝑏1 < 1 𝛼𝑦 < 1 No oscilatoria y no fluctuante
𝒚>
(𝟏 + 𝜶)𝟐
1𝐷: 1 < 𝑏2 < 𝑏1 𝛼𝑦 > 1

[Link]íces reales

repetidas

𝟒𝜶 2𝐶: 0 < 𝑏 < 1 𝛼𝑦 < 1 No oscilatoria y no fluctuante

𝒚=
(𝟏 + 𝜶)𝟐
2𝐷: 𝑏 > 1 𝛼𝑦 > 1

[Link]

complejas

𝟒𝜶 3𝐶: 𝑅 < 1 𝛼𝑦 < 1 Con fluctuación escalonada

𝒚<
(𝟏 + 𝜶)𝟐
3𝐷: ≥ 1 𝛼𝑦 ≥ 1

Finalmente, en el caso 3, con raíces complejas, tenemos fluctuación escalonada y por lo tanto ciclos

de negocios endógenos. En este caso, debemos ver el valor absoluto 𝑅 = √𝑎2 vea (18.2) para la

pista de convergencia y la divergencia, donde 𝑎2 es el coeficiente del término 𝑦𝑡 en la ecuación en

diferencias (18.1). En el presente modelo, tenemos 𝑅 = √𝛼𝑦 , lo que da lugar a las tres siguientes

posibilidades:

𝑖𝑥) 𝑅 < 1 → 𝛼𝑦 < 1

𝑥) 𝑅 = 1 → 𝛼𝑦 = 1

𝑥𝑖) 𝑅 > 1 → 𝛼𝑦 > 1

Aun cuando todos estos resultan ser admisibles (vea el ejercicio 18.2-4), solo la posibilidad 𝑅 < 1

implica una trayectoria de tiempo convergente y califica como el sub caso 3𝐶 en la tabla 18.1. Las

otras dos se rotulan colectivamente como el sub caso 3𝐷.

En resumen, de la tabla 18.1 podemos concluir que puede presentarse una trayectoria de tiempo

convergente si y solo si 𝛼𝑦 < 1.

d) Un resumen gráfico

El análisis anterior nos ha conducido a una clasificación algo compleja de casos y sub casos. Nos

ayudaría a tener una representación visual del esquema de clasificación. Esto se ofrece en la figura

18.2.

El conjunto de todos los pares admisibles (𝛼, 𝑦) en el modelo se muestran en la figura 18.2 por el

área rectangular sombreada de varias maneras. Dado que se excluyen los valores de 𝑦 = 0 y 𝑦 =

1, así como el valor 𝛼 = 0, el área sombreada es una especie de rectángulo sin lados. Ya hemos

graficado la ecuación 𝑦 = 4𝛼/(1 + 𝛼)2 para delimitar los tres casos principales de la tabla 18.1:

los puntos sobre esa curva pertenecen al caso 2; los puntos situados al norte de la curva (que

representan valores más altos de 𝑦) pertenecen al caso 1; los situados al sur (con valores más bajos

de 𝑦) son del caso 3. Para distinguir entre los sub casos convergente y divergente, ahora añadimos

la gráfica de 𝛼𝑦 = 1(una hipérbola rectangular) como otra línea de demarcación. Los puntos

situados al norte de la hipérbola rectangular satisfacen la desigualdad 𝛼𝑦 > 1, mientras que los

ubicados por debajo de aquella corresponden a 𝛼𝑦 < 1. Entonces podemos demarcar fácilmente

los sub casos. Para el caso 1, la región sombreada con líneas discontinuas, que está por debajo de

la hipérbola, corresponde al sub caso 1𝐶,pero la región sombreada con líneas continuas se asocia

con el sub caso 1𝐷. Para el caso 2, que se relaciona con los puntos situados sobre la curva 𝑦 =
4𝛼/(1 + 𝛼)2 , el sub caso 2𝐶 cubre la parte con pendiente hacia arriba de esa curva, y el sub caso

2𝐷las partes con pendiente hacia abajo. Finamente, para el caso 3, la hipérbola rectangular sirve

para separar la región sombreada llana (sub caso 3𝐶) de la región sombreada y granulada (sub caso

3𝐷). Debe observar que esta última también incluye loas puntos localizados sobre la hipérbola

rectangular misma, debido a la desigualdad débil en la especificación 𝛼𝑦 ≥ 1.

Como la figura 18.2 es la depositaria de todas las conclusiones cualitativas del modelo, dado

cualquier par ordenado (𝛼, 𝑦),siempre podemos encontrar gráficamente al sub caso correcto al

graficar el par ordenado en el diagrama.

Por analogía con los casos de oscilación uniforme y fluctuación uniforme, podemos denominar

esta situación como “fluctuación escalonada uniforme”. Sin embargo, no podemos esperar en

general el aspecto de uniformidad en este último caso sea perfecto porque, en forma similar a lo

que se hizo en la figura 1.18, podemos aceptar solo aquellos puntos sobre una curva seno o coseno

que correspondan a valores enteros de 𝑡1 pero estos pueden inducir sobre un conjunto de puntos

totalmente diferentes sobre la curva para cada periodo de fluctuación.

Ejemplo

1. Si el acelerador es 0.8 y la propensión marginal al consumo es 0.7, ¿Qué tipo de trayectoria de

tiempo de interacción va a resultar? El par ordenado (0.8, 0.7) se ubica en la región de gris liso,

sub caso 3𝐶; entonces la trayectoria de tiempo se caracteriza por la fluctuación escalonada

amortiguada.

2. ¿Qué tipo de interacción está implicada con 𝛼 = 2 y 𝑦 = 0.5? el par ordenado (2,0.5) esté

situado exactamente sobre la hipérbola rectangular, bajo el sub caso 3𝐷. La trayectoria de tiempo

de 𝑌 nuevamente exhibe fluctuación escalonada, pero ésta no va ser ni explosiva ni amortiguada.

Ejercicio

1. Consulte la figura 18.2 y encuentre los sub casos a los que pertenecen los siguientes

conjuntos de valores de 𝛼 y 𝑦, y describa de forma cualitativa la trayectoria de tiempo de

interaccion.

𝑎) 𝛼 = 3.5; 𝑦 = 0.8 𝑐) 𝛼 = 0.2; 𝑦 = 0.9

𝑏) 𝛼 = 2; 𝑦 = 0.7 𝑑) 𝛼 = 1.5; 𝑦 = 0.6

2. De los valores de 𝛼 y 𝑦 dados en las partes 𝑎) y 𝑐) del problema 1, encuentre los valores

numéricos de las raíces características en cada caso, y analice la naturaleza de la trayectoria

de tiempo. ¿concuerdan sus resultados con los obtenidos anteriormente?

3. Verifique las posibilidades 𝑖𝑖, 𝑖𝑖𝑖 y 𝑖𝑣 en el caso 1 implican valores inadmisibles de 𝑦.

4. Muestre que en el caso 3 nunca podemos encontrar 𝑦 ≥ 1.

3. LA INFLACIÓN Y EL DESEMPLEO EN TIEMPO DISCRETO

La interacción de la inflación y del desempleo, estudiada antes en el marco de referencia de tiempo

continuo, también puede abordarse en el tiempo discreto. Usando esencialmente las mismas

hipótesis económicas, ilustraremos en esta sección la forma en que puede reformularse este modelo

como un modelo de ecuación en diferencias.

𝑤 = 𝑓(𝑈) 𝑝 =𝑤−𝑇 𝑝 = 𝛼 − 𝑇 − 𝛽𝑈

a) El modelo

La formulación anterior en tiempo continuo de la interacción de la inflación y el desempleo

consistía en tres ecuaciones diferenciales:

 Relación de Phillips con expectativas:

𝑝 = 𝛼 − 𝑇 − 𝛽𝑈 + 𝑔𝜋 (0 < 𝑔 ≤ 1)

 Expectativas adaptativas:

𝑑𝜋
= 𝑗(𝑝 − 𝜋) (0 < 𝑗 ≤ 1)
𝑑𝑡

 Política monetaria:

𝑑𝑈
= −𝑘(𝑚 − 𝑝) (𝑘 > 0)
𝑑𝑡

Están presentes tres variables endógenas: 𝑝 (tasa real de inflación), π (tasa esperada de inflación)

y 𝑈 (tasa de desempleo), Aparecen hasta seis parámetros en el modelo; entre estos, el parámetro 𝑚

la tasa de crecimiento del dinero nominal (o sea la tasa de expansión monetaria) difiere de los

demás en que su magnitud se establece como una decisión de política. W (tasa de crecimiento del

salario monetario)

Cuando se vierte en el molde del análisis de periodos, la relación de Phillips con expectativas

simplemente se transforma en:

𝑝𝑡 = 𝛼 − 𝑇 − 𝛽𝑈𝑡 + 𝑔𝜋𝑡 (𝛼, 𝛽 > 0 ; 0 < 𝑔 ≤ 1)

En la ecuación de las expectativas adaptativas, la derivada debe reemplazarse por una expresión de

diferencias:

𝜋𝑡+1 − 𝜋𝑡 = 𝑗(𝑝𝑡 − 𝜋𝑡 ) (0 < 𝑗 ≤ 1)

Por la misma razón, la ecuación de política monetaria debe cambiarse a…

𝑈𝑡+1 − 𝑈𝑡 = −𝑘(𝑚 − 𝑝𝑡+1 ) (𝑘 > 0)

Estas tres ecuaciones constituyen la nueva versión del modelo de inflación – desempleo.

b) La ecuación en diferencias en 𝑝
Como primer paso en el análisis de este nuevo modelo, una vez más tratamos de condensar el

modelo en una sola ecuación con una sola variable. Sea esa variable 𝑝 . De acuerdo con esto,

centraremos nuestra atención en el molde de análisis de periodos de la relación de Phillips. Sin

embargo, ya que el molde de análisis de periodos de la relación de Phillips a diferencia de las otras

dos ecuaciones no describe por sí misma un patrón de cambio, depende de nosotros la creación de

este patrón. Esto se logra al 𝑜𝑏𝑡𝑒𝑛𝑒𝑟 𝑙𝑎𝑠 𝑑𝑖𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎𝑠 𝑑𝑒 𝑝𝑡 , es decir, al tomar la primera

diferencia de 𝑝𝑡 de acuerdo con la definición.

∆𝑝𝑡 ≡ 𝑝𝑡+1 − 𝑝𝑡

Intervienen dos pasos en esto. Primero, desplazamos hacia delante los subíndices de tiempo en el

molde de análisis de periodos de la relación de Phillips un periodo, para obtener.

𝑝𝑡+1 = 𝛼 − 𝑇 − 𝛽𝑈𝑡+1 + 𝑔𝜋𝑡+1

Entonces restamos el molde de análisis de periodos de la relación de Phillips de la formula anterior,

para obtener la primera diferencia de 𝑝𝑡 queda el patrón deseado de cambio:

𝒑𝒕+𝟏 − 𝒑𝒕 = −𝜷(𝑼𝒕+𝟏 − 𝑼𝒕 ) + 𝒈(𝝅𝒕+𝟏 − 𝝅𝒕 )

Reemplazando:

𝑈𝑡+1 − 𝑈𝑡 = −𝑘(𝑚 − 𝑝𝑡+1 ) 𝜋𝑡+1 − 𝜋𝑡 = 𝑗(𝑝𝑡 − 𝜋𝑡 )

𝒑𝒕+𝟏 − 𝒑𝒕 = 𝜷𝒌(𝒎 − 𝒑𝒕+𝟏 ) + 𝒈𝒋(𝒑𝒕 − 𝝅𝒕 )

Observe que en 𝑝𝑡+1 − 𝑝𝑡 = 𝛽𝑘(𝑚 − 𝑝𝑡+1 ) + 𝑔𝑗(𝑝𝑡 − 𝜋𝑡 ) los patrones de cambio de las otras

dos variables tal como están dadas por 𝑈𝑡+1 − 𝑈𝑡 = −𝑘(𝑚 − 𝑝𝑡+1 ) y

𝜋𝑡+1 − 𝜋𝑡 = 𝑗(𝑝𝑡 − 𝜋𝑡 ) Se han incorporado al patrón de cambio de la variable 𝑝 . Entonces

𝑝𝑡+1 − 𝑝𝑡 = 𝛽𝑘(𝑚 − 𝑝𝑡+1 ) + 𝑔𝑗(𝑝𝑡 − 𝜋𝑡 ) Incorpora ahora toda la información del presente

modelo.

Sin embargo, el término 𝜋𝑡 es ajeno al estudio de 𝑝 y es necesario eliminarlo de

𝑝𝑡+1 − 𝑝𝑡 = 𝛽𝑘(𝑚 − 𝑝𝑡+1 ) + 𝑔𝑗(𝑝𝑡 − 𝜋𝑡 ) . Con este fin, hacemos uso del hecho de que…

𝒈𝝅𝒕 = 𝒑𝒕 − (𝜶 − 𝑻) + 𝜷𝑼𝒕 Mediante 𝒑𝒕 = 𝜶 − 𝑻 − 𝜷𝑼𝒕 + 𝒈𝝅𝒕

Al sustituir esto en 𝑝𝑡+1 − 𝑝𝑡 = 𝛽𝑘(𝑚 − 𝑝𝑡+1 ) + 𝑔𝑗(𝑝𝑡 − 𝜋𝑡 ) y reducir términos, obtenemos:

(𝟏 + 𝜷𝒌)𝒑𝒕+𝟏 − [𝟏 − 𝒋(𝟏 − 𝒈)]𝒑𝒕 + 𝒋𝜷𝑼𝒕 = 𝜷𝒌𝒎 + 𝒋(𝜶 − 𝑻)

Pero ahora aparece un término 𝑈𝑡 que debe eliminarse. Para hacerlo, diferenciamos lo último para

obtener un término (𝑈𝑡+1 − 𝑈𝑡 ) y luego usar 𝑈𝑡+1 − 𝑈𝑡 = −𝑘(𝑚 − 𝑝𝑡+1 ) para eliminar este

último. Solo después de este proceso más bien engorroso de sustituciones, obtendremos la ecuación

en diferencias deseada solo en la variable 𝑝 , la cual, cuando se normaliza apropiadamente, adopta

la forma…

1 + 𝑔𝑗 + (1 − 𝑗)(1 + 𝛽𝑘) 1 − 𝑗(1 − 𝑔) 𝑗𝛽𝑘𝑚

𝑝𝑡+2 − 𝑝𝑡+1 + 𝑝𝑡 =
1 + 𝛽𝑘 1 + 𝛽𝑘 1 + 𝛽𝑘

a1 a2 a3

c) La trayectoria de tiempo de p

El valor de equilibrio intertemporal de p, dado por la integral particular de este último, es

𝑐 𝑗𝛽𝑘𝑚
𝑝̅ = = =𝑚
1 + 𝑎1 + 𝑎2 𝛽𝑘𝑗

Mediante:
Siempre que (1 + 𝑎1 + 𝑎2 ) ≠ 0, la integral particular será:

𝐶
𝑦𝑝 = 𝑘 =
1 + 𝑎1 + 𝑎2

En caso de 𝑎1 + 𝑎2 ≠ −1

Por lo tanto, al igual que en el modelo de tiempo continuo, la tasa de inflación de equilibrio es

exactamente igual a la de expansión monetaria.

En cuanto a la función complementaria, pueden surgir ya sea raíces reales diferentes (caso 1), o

raíces repetidas (caso 2), o raíces complejas (caso 3), dependiendo de la magnitud relativa de 𝑎12

y 4𝑎2 . En el presente modelo,

≥ ≥
𝑎12 < 4𝑎2 𝑠𝑖 𝑦 𝑠𝑜𝑙𝑜 𝑠𝑖 [1 + 𝑔𝑗 + (1 − 𝑗)(1 + 𝛽𝑘)]2 <
4[1 − 𝑗(1 − 𝑔)](1 + 𝛽𝑘)

1 1 1
Si 𝑔 = ,𝑗 = y 𝛽𝑘 = 5 , por ejemplo, entonces 𝑎12 = (5 )2 mientras que 4𝑎2 = 20 ; por lo
2 3 6

tanto resulta el caso 1. Pero si g = j = 1, entonces 4𝑎2 = 4(1 + 𝛽𝑘) > 4 , y tenemos el caso 3.

Sin embargo, en vista del mayor número de parámetros en el presente modelo, no es factible

construir una gráfica de clasificación como en el modelo de Samuelson.

Sin embargo, todavía puede proceder el análisis de convergencia a lo largo de la misma línea que

en el modelo de interacción de multiplicador con acelerador de Samuelson. Específicamente,

−𝑎1 +√𝑎12 −4𝑎2 −𝑎1 −√𝑎12 −4𝑎2 −2𝑎1

recordamos 𝑟1 + 𝑟2 = + = = −𝑎1 que las dos raíces
2 2 2

características b1 y b2 deben satisfacer las dos siguientes relaciones:

1 + 𝑔𝑗
𝑏1 + 𝑏2 = −𝑎1 = +1−𝑗 >0
1 + 𝛽𝑘
 1 − 𝑗(1 − 𝑔)
𝑏1 𝑏2 = 𝑎12 = ∈ (0,1)
1 + 𝛽𝑘

Aún más, en el presente modelo tenemos

𝛽𝑗𝑘
(1 − 𝑏1 )(1 − 𝑏2 ) = 1 − (𝑏1 + 𝑏2 ) + 𝑏1 𝑏2 = >0
1 + 𝛽𝑘

Ahora considere el caso 1, donde las dos raíces b1 y b2 son reales y diferentes. Como el producto

𝑏1 𝑏2 es positivo, b1 y b2 deben tener el mismo signo. Aun mas, como la suma es positiva, b1 y b2

deben ser ambas positivas, lo que implica que no puede ocurrir oscilación. A partir de

𝛽𝑗𝑘
(1 − 𝑏1 )(1 − 𝑏2 ) = 1 − (𝑏1 + 𝑏2 ) + 𝑏1 𝑏2 = > 0 Podemos inferir que ni b1 ni b2 pueden
1+𝛽𝑘

ser iguales a uno; ya que de otro modo (1 − 𝑏1 )(1 − 𝑏2 ) seria cero, violando la desigualdad

indicada. Esto significa que, en términos de las diferentes posibilidades de las combinaciones de

(b1, b2) enumeradas en el modelo de Samuelson, las posibilidades ii y iv no pueden surgir aquí.

También es inaceptable tener una raíz mayor que uno y la otra raíz menor que uno; puesto que de

otra manera (1 − 𝑏1 )(1 − 𝑏2 ) sería negativo. Entonces también se descarta la posibilidad iii. Se

sique que b1 y b2 deben ser ya sea los dos mayores que uno, o ambos menores que uno. Si b1 > 1

1−𝑗(1−𝑔)
y b2 > 1 (posibilidad v) sin embargo, se violaría 𝑏1 𝑏2 = 𝑎12 = ∈ (0,1) . Por eso lo
1+𝛽𝑘

único viable es la posibilidad i, con b1 y b2 siendo ambas fracciones positivas, de modo que la

trayectoria de tiempo de p es convergente.

El análisis de caso 2 no es muy diferente. Mediante un razonamiento prácticamente idéntico,

podemos concluir que la raíz repetida b solo puede ser una fracción positiva en este modelo; es

decir, la posibilidad vi es factible, pero no las posibilidades vii y viii. La trayectoria de tiempo de

p en el caso 2 es nuevamente no oscilatoria y convergente.

Para el caso 3, la convergencia requiere que R (el valor absoluto de las raíces complejas) sea menor

𝑎12 +4𝑎2 −𝑎12

que uno. Por 𝑅 = √ℎ2 + 𝑣 2 = √ = √𝑎2 , R = √𝑎2 . Siempre que 𝑎2 sea una fracción
4

1−𝑗(1−𝑔)
positiva 𝑏1 𝑏2 = 𝑎12 = ∈ (0,1) ciertamente tenemos R < 1. Entonces, la trayectoria de
1+𝛽𝑘

tiempo de p en el caso 3 también es convergente, aunque esta vez habrá una fluctuación escalonada.

d) El Análisis De U

Para analizar esta vez la trayectoria de tiempo del desempleo, podemos tomar como punto

de partida:

𝑈𝑡+1 − 𝑈𝑡 = −𝑘(𝑚 − 𝑝𝑡+1 ) (18.20)

Para deshacernos de p en esta ecuación sustituimos primero:

𝑝𝑡+1 =∝ −𝑇 − 𝛽𝑈𝑡+1 + 𝑔𝜋𝑡+1 (18.18´)

Para obtener:

(1 + 𝛽𝑘)𝑈𝑡+1 − 𝑈𝑡 = 𝑘(𝛼 − 𝑇 − 𝑚) + 𝑘𝑔𝜋𝑡+1 (18.28)

En seguida para prepararse para la sustitución de la otra ecuación, al tomar la diferencia

encontramos:

(1 + 𝛽𝑘)𝑈𝑡+2 − (2 + 𝛽𝑘)𝑈𝑡+1 + 𝑈𝑡 = 𝑘𝑔(𝜋𝑡+2 − 𝜋𝑡+1 ) (18.29)

En vista de la presencia de una expresión en diferencias en 𝜋 a la derecha; podemos sustituirla

por una versión desplazada

hacia delante de la ecuación de expectativas adaptativas.

(1 + 𝛽𝑘)𝑈𝑡+2 − (2 + 𝛽𝑘)𝑈𝑡+1 + 𝑈𝑡 = 𝑘𝑔𝑗(𝑝𝑡+1 − 𝜋𝑡+1 ) (18.30)

Se debe eliminar las variables p y 𝜋 antes de que surja una ecuación en diferencias propias en U.

𝑘𝑝𝑡+1 = 𝑈𝑡+1 − 𝑈𝑡 + 𝑘𝑚 (18.31)

Aún más al multiplicar g𝜋𝑡 = 𝑝𝑡 − (𝛼 − 𝑇) + 𝛽𝑈𝑡 por (-kj) y al desplazar los subíndices de

tiempo; obtenemos:

−𝑘𝑗𝑞𝜋𝑡+1 = −𝑘𝑗𝑝𝑡+1 + 𝑘𝑗(𝛼 − 𝑇) − 𝛽𝑘𝑗𝑈𝑡+1

= −𝑗(𝑈𝑡+1 − 𝑈𝑡 + 𝑘𝑚) + 𝑘𝑗(𝛼 − 𝑇) − 𝛽𝑘𝑗𝑈𝑡+1

= −𝑗(1 + 𝛽𝑘)𝑈𝑡+1 + 𝑗𝑈𝑡 + 𝑘𝑗(𝛼 − 𝑇 − 𝑚)

Estos dos resultados expresan a 𝑝𝑡+1 y 𝜋𝑡+1 en términos de la variable U y por lo tanto nos

permite obtener la ecuación en diferencias deseada solo en la variable U.

1 + 𝑔𝑗 + (1 − 𝑗)(1 + 𝛽𝑘) 1 − 𝑗(1 − 𝑔)

𝑈𝑡+2 − 𝑈𝑡+1 + 𝑈𝑡
1 + 𝛽𝑘 1 + 𝛽𝑘

𝑘𝑗[𝛼−𝑇−(1−𝑔)𝑚]
=
1+𝛽𝑘

e) La Relación De Phillips De Largo Plazo

La tasa de equilibrio intertemporal del desempleo es:

 1
̅=
𝑈 [𝛼 − 𝑇 − (1 − 𝑔)𝑚]
𝛽

Pero como se halló que la tasa de inflación es igual a la tasa de expansión monetaria 𝑝̅ =m,

̅ 𝑐𝑜𝑛 𝑝̅ mediante la ecuación:

podemos enlazar a 𝑈

̅ = 1 [𝛼 − 𝑇 − (1 − 𝑔)𝑝̅ ]
𝑈 (18.34)
𝛽

Dado que esta ecuación tiene que ver solo con las tasas de desempleo y de inflación de equilibrio,

se dice que ilustra la relación de Phillips a largo plazo.

Un caso especial de (18.34) ha recibido mucha atención entre los economistas: el caso de g = 1. Si

g = 1, el término 𝑝̅ tendrá coeficiente cero y por tanto sale de la escena. En otras palabras, U va a

transformarse en una función constante de 𝑝̅ .

En el diagrama estándar de Phillips, donde la tasa de desempleo se gráfica en el eje horizontal,

este resultado da lugar a una curva vertical de Phillips de largo plazo. El valor de U en este caso,

denominado tasa natural de desempleo, es consistente con cualquier tasa de inflación de equilibrio,

con la notable implicación de política de que, a largo plazo, no hay intercambios entre los demonios

gemelos de la inflación y del desempleo tal como sucede a corto plazo.

PERO, ¿QUÉ OCURRE SI g < 1?

En ese caso, el coeficiente de 𝑝̅ de (18.34) será negativo. Entonces, la curva de Phillips de largo

plazo tendrá pendiente hacia abajo, suministrando con ello todavía una relación de intercambio

entre inflación y desempleo.

Ya sea que la curva de Phillips de largo plazo tenga pendiente vertical o negativa, depende en forma

crítica por lo tanto del valor del parámetro g, el cual, de acuerdo con la relación de Phillips con

expectativas, mide hasta qué punto la tasa de inflación esperada logra influir en la estructura de los

salarios y en la tasa de inflación real.

Generalizaciones a ecuaciones con términos variables y de orden superior

El termino variable con forma de Cmt.

Cuando el termino constante c se reemplaza por un término variable el único efecto recaerá sobre

la solución particular. Para encontrar la nueva solución particular, podemos aplicar el método de

coeficientes indeterminados. Para este método se requiere que el término variable y sus derivadas

sucesivas juntas tomen solo un número finito de tipos diferentes de expresiones, aparte de las

constantes multiplicativas.

Ejemplo 1

𝑦𝑡+2 + 𝑦𝑡+1 − 3𝑦𝑡 = 7𝑡

Aquí tenemos c=1 y m=7. Evaluemos primero si el término variable 7𝑡 suministra un número

finito de tipos de expresiones de diferencias sucesivas. De acuerdo con la regla para tomar

diferencias (∆𝑦𝑡 = 𝑦𝑡+1 − 𝑦𝑡 )

La primera diferencia del término es:

∆7𝑡 = 7𝑡+1 − 7𝑡 = (7 − 1)7𝑡 = 6(7)𝑡

La segunda diferencia, ∆2 (7𝑡 ), puede expresarse como:

 ∆(∆7𝑡 ) = ∆6(7𝑡 ) = 6(7)𝑡+1 − 6(7)𝑡 = 6(7 − 1)7𝑡 = 36(7)𝑡

Todas las diferencias sucesivas serán algún múltiplo de 7𝑡 como la primera y la segunda. Como

solamente hay un tipo de expresión, podemos intentar una solución 𝑦𝑡 = 𝐵(7)𝑡 para la solución

particular, donde B es un coeficiente indeterminado.

Sustituyendo la solución de prueba y sus versiones correspondientes los periodos (𝑡 + 1) y

(𝑡 + 2) en la ecuación en diferencia dada, obtenemos:

𝐵(7)𝑡+2 + 𝐵(7)𝑡+1 − 3𝐵(7)𝑡 = 7𝑡 o 𝐵(72 + 7 − 3)(7)𝑡 = 7𝑡

Despejando B

1 1
𝐵= =
49 + 7 − 3 53

Y podemos escribir la solución participar como

1
𝑦𝑝 = 𝐵(7)𝑡 = (7)𝑡
53

Esto puede tomarse como un equilibrio móvil. Se puede verificar la corrección de la solución al

sustituirla en la ecuación en diferencias y comprobar que va a obtener una identidad 7𝑡 = 7𝑡

Generalizaciones

El resultado obtenido en el ejemplo 1 puede generalizarse fácilmente del termino variable 7𝑡 al de

𝑐𝑚𝑡 . De nuestra experiencia, esperamos que todas las diferencias sucesivas de 𝑐𝑚𝑡 adopten la

misma forma de expresión, es decir, 𝐵𝑚𝑡 , donde B es una constante multiplicativa.

Entonces podemos intentar una solución 𝑦𝑡 = 𝐵𝑚𝑡 para la solución particular, cuando se da la

ecuación en diferencias

𝑦𝑡+2 + 𝑎1 𝑦𝑡+1 − 𝑎2 𝑦𝑡 = 𝑐𝑚𝑡

Usando la solución de prueba 𝑦𝑡 = 𝐵𝑚𝑡 , que implica que 𝑦𝑡+1 = 𝐵𝑚𝑡+1 , etc., podemos expresar

la ecuación como

𝐵𝑚𝑡+2 + 𝑎1 𝐵𝑚𝑡+1 + 𝑎2 𝐵𝑚𝑡 = 𝑐𝑚𝑡

Es decir,

𝐵(𝑚2 + 𝑎1 𝑚 + 𝑎2 )𝑚𝑡 = 𝑐𝑚𝑡

Entonces, el coeficiente B en la solución de prueba debe ser

𝑐
𝐵=
𝑚2 + 𝑎1 𝑚 + 𝑎2

Y la solución particular deseada desde de la ecuación puede escribirse como

𝑐
𝑦𝑝 = 𝐵𝑚𝑡 = 𝑚2 +𝑎 𝑚𝑡 (𝑚2 + 𝑎1 𝑚 + 𝑎2 ≠ 0)
1 𝑚+𝑎2

No se permite que el denominador de B sea cero. SI lo es, entonces debemos usar en su lugar la

solución de prueba 𝑦𝑡 = 𝐵𝑚𝑡 ; o si esta también falla, 𝑦𝑡 = 𝐵𝑡 2 𝑚𝑡 .

Ejercicios propuestos
Encuentre la solución particular de:

a) 𝒚𝒕+𝟐 − 𝟐𝒚𝒕+𝟏 + 𝟓𝒚𝒕 = 𝟑𝒕

Donde c=1 y m=3 y los coeficientes son: 𝑎0 = 1, 𝑎1 = −2, 𝑎2 = 5

Aplicando la fórmula:

1
𝑦𝑝 = 𝐵𝑚𝑡 = 2
3𝑡
3 + (−2)(3) + 5

1
𝑦𝑝 = 3𝑡
8

El termino variable con forma de 𝒄𝒕𝒏

Tenemos el término variable 𝑐𝑡 𝑛 , donde c es cualquier constante y n es un entero positivo.

Ejemplo 2

Encuentre la solución de

𝑦𝑡+2 + 5𝑦𝑡+1 + 2𝑦𝑡 = 𝑡 2

Las tres primeras diferencias de 𝑡 2 (un caso especial de 𝑐𝑡 𝑛 con c=1 y n=2) se encuentran como

sigue.

∆𝑡 2 = (𝑡 + 1)2 − 𝑡 2 = 2𝑡 + 1

∆2 𝑡 2 = ∆(∆𝑡 2 ) = ∆(2𝑡 + 1) = ∆2𝑡 + ∆1 = 2(𝑡 + 1) − 2𝑡 + 0 = 2 [∆ 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 = 0]

∆3 𝑡 2 = ∆(∆2 𝑡 2 ) = ∆2 = 0
Como las diferencias adicionales solo van a dar cero, hay en conjunto tres tipos diferentes de

expresiones 𝑡 2 (del termino variable mínimo), t y una constante (a partir de las diferencias

sucesivas).

Por lo tanto intentemos la solución

𝑦𝑡 = 𝐵0 + 𝐵1 𝑡 + 𝐵2 𝑡 2

Para la solución particular, con coeficientes indeterminados 𝐵0 , 𝐵1 y 𝐵2. Observe que esta

solución implica

𝑦𝑡+1 = 𝐵0 + 𝐵1 (𝑡 + 1) + 𝐵2 (𝑡 + 1)2 = (𝐵0 + 𝐵1 + 𝐵2 ) + (𝐵1 + 2𝐵2 )𝑡 + 𝐵2 𝑡 2

𝑦𝑡+2 = 𝐵0 + 𝐵1 (𝑡 + 2) + 𝐵2 (𝑡 + 2)2 = (𝐵0 + 2𝐵1 + 4𝐵2 ) + (𝐵1 + 4𝐵2 )𝑡 + 𝐵2 𝑡 2

Cuando esto se sustituye en la ecuación en diferencias, tenemos

(8𝐵0 + 7𝐵1 + 9𝐵2 ) + (8𝐵1 + 14𝐵2 )𝑡 + 8𝐵2 𝑡 2 = 𝑡 2

Al igualar los dos lados término a término, vemos que se requiere que los coeficientes

indeterminados satisfagan las siguientes ecuaciones simultáneas:

8𝐵0 + 7𝐵1 + 9𝐵2 = 0

8𝐵1 + 14𝐵2 = 0

8𝐵2 = 1

Los valores obtenidos son:

13
𝐵0 =
256
 7
𝐵1 = −
32

1
𝐵2 =
8

Dándonos la solución particular

13 7 1
𝑦𝑝 = − 𝑡 + 𝑡2
256 32 8

Generalizaciones

En la nueva solución de prueba, obviamente debe estar el término 𝐵𝑛 𝑡 𝑛 , que corresponda al

término variable dado. Aun más, ya que la diferenciación sucesiva del término suministra las

expresiones diferentes 𝑡 𝑛−1 , 𝑡 𝑛−2 , … , 𝑡 y 𝐵0 (constante), la nueva solución de prueba para el caso

del término variable 𝑐𝑡 𝑛 debe escribirse como

𝑦𝑡 = 𝐵0 + 𝐵1 𝑡 + 𝐵2 𝑡 2 + ⋯ + 𝐵𝑛 𝑡 𝑛

Esta solución de prueba también puede fallar. En ese caso, el truco nuevamente es multiplicar la

solución original de prueba por una potencia suficientemente alta de t, es decir, podemos intentar

en su lugar

𝑦𝑡 = 𝑡(𝐵0 + 𝐵1 𝑡 + 𝐵2 𝑡 2 + ⋯ + 𝐵𝑛 𝑡 𝑛 )

Ecuaciones en diferencias lineales de orden superior

El orden de una ecuación en diferencias indica la diferencia de mayor orden presente en la

ecuación; pero también indica el número máximo considerado de periodos de desfasamiento de

tiempo. Una ecuación en diferencias lineal de orden n (con coeficientes constantes y término

constante) puede entonces escribirse en general como

𝑦𝑡+𝑛 + 𝑎1 𝑦𝑡+𝑛−1 + ⋯ + 𝑎𝑛−1 𝑦𝑡+1 + 𝑎𝑛 𝑦𝑡 = 𝑐

Las ecuaciones en diferencia lineales se pueden clasificar en:

• Homogéneas si c = 0

• Completas si c ≠ 0

• De coeficientes constantes si 𝑎𝑖

• De coeficientes no constante si 𝑎𝑖 (𝑡)

Enunciamos una serie de teoremas que nos guían en la búsqueda de la solución general de este

tipo de ecuaciones.

Teorema 1 (Teorema de existencia y unicidad) Dada la ecuación y dados n números reales k0,

k1, . . . , kn−1 existe una única solución que verifica

𝑦(0) = 𝑘0 , 𝑦(1) = 𝑘1 , … 𝑦(𝑛 − 1) = 𝑘𝑛−1

Demostración: Definimos la sucesión 𝑦(𝑡) de la siguiente forma

𝑦0 = 𝑘0

𝑦1 = 𝑘1

.
.
.
 𝑦𝑛−1 = 𝑘𝑛−1

𝑦𝑛 = −𝑎1 (0)𝑦𝑛−1 − … − 𝑎𝑛 (0)𝑦0

𝑦𝑛+1 = −𝑎1 (1)𝑦𝑛 − … − 𝑎𝑛 (1)𝑦1

.
.
.
𝑦𝑡 = −𝑎1 (𝑡 − 𝑛)𝑦𝑡−1 − ⋯ − 𝑎𝑛 (𝑡 − 𝑛)𝑦𝑡−𝑛

Es decir, {y (t)} se define mediante una ley de recurrencia. Además {y (t)} es solución de la

ecuación y satisface las condiciones iniciales y además es única, ya que si existe otra solución que

verifique las condiciones iniciales, entonces coinciden en todos los puntos pues la propia ley de

recurrencia determina los valores posteriores de la solución.

Para cada conjunto de condiciones iniciales {k0, k1, . . . , kn−1} hay una única solución, pero al

variar las condiciones iniciales varia la solución, por tanto una ecuación en diferencias lineal

homogénea tiene infinitas soluciones dependiendo de las condiciones iniciales dadas. De ahí, la

importancia del siguiente teorema.

Teorema 2: Toda combinación lineal de soluciones de una ecuación en diferencias lineal

homogénea de orden n es también solución.

Hallando la Solución General

Para encontrar la solución particular podemos intentar 𝑦𝑡 = 𝑘 (el caso de equilibrio intertemporal

estacionario). Si esto falla podemos intentar 𝑦𝑡 = 𝑘𝑡 o 𝑦𝑡 = 𝑘𝑡 2 , etc., en ese orden.

La función complementaria va ser una ecuación polinomial de grado n:

𝑏 𝑛 + 𝑎1 𝑏 𝑛−1 + ⋯ + 𝑎𝑛−1 𝑏 + 𝑎𝑛 = 0
Habrá n raíces características 𝑏𝑖 (𝑖 = 1,2, … , 𝑛), todo lo cual debe entrar en la función

complementaria, por lo tanto:

𝑦𝑐 = ∑ 𝐴𝑖 𝑏𝑖𝑡
𝑖=1

Suponiendo que las raíces son todas reales y diferentes.

 En caso de que haya raíces reales repetidas (digamos, 𝑏1 = 𝑏2 = 𝑏3 ) entonces los tres
primeros términos de la suma de la ecuación deben modificarse como
𝐴1 𝑏1𝑡 + 𝐴2 𝑡𝑏1𝑡 + 𝐴3 𝑡 2 𝑏1𝑡

 Si existe un par de raíces complejas conjugadas—digamos, 𝑏𝑛−1 , 𝑏𝑛 , entonces los dos

últimos términos de la suma de la ecuación deben combinarse según la expresión
𝑅 𝑡 (𝐴𝑛−1 cos 𝜃𝑡 + 𝐴𝑛 𝑠𝑒𝑛 𝜃𝑡)

También puede asignarse una expresión similar a cualquier otro par de raíces complejas. Sin

embrago, en caso de dos pares repetidos, a uno de los dos debe dársele un factor multiplicativo de

𝑡𝑅 𝑡 en lugar de 𝑅 𝑡 .

Después de haber encontrado 𝑦𝑝 y 𝑦𝑐 la solución general de la ecuación en diferencias completa

se obtiene mediante una suma como la siguiente:

𝑦𝑡 = 𝑦𝑝 + 𝑦𝑐

La convergencia y el teorema de Schur

Cuando se tiene una ecuación en diferencias de orden superior que no se resuelva fácilmente,

podemos determinar la convergencia de la trayectoria de tiempo relevante en forma cualitativa sin

tener que batallar con la solución cuantitativa real.

La trayectoria de tiempo puede convergir si y solo si cada una de las raíces de la ecuación

característica es menor que 1 en valor absoluto. El siguiente teorema llamado también como

Teorema de Schur, se aplica en ese caso.

Las raíces de la ecuación polinomial de grado n

𝑎0 𝑏 𝑛 + 𝑎1 𝑏 𝑛−1 + ⋯ + 𝑎𝑛−1 𝑏 + 𝑎𝑛 = 0

Serán todas menores que la unidad en valor absoluto si y solo los siguientes n determinantes

SON TODOS POSITIVOS

El teorema de Schur es una contraparte perfecta de la ecuación en diferencias del teorema de Routh

introducido antes en el marco de las ecuaciones diferenciales.

Se explica de la mejor manera con ayuda de las líneas segmentadas que forman las participaciones

de cada determinante para formar cuatro áreas. Cada área del determinante k-esimo, ∆𝑘 , siempre

consiste de un subdeterminante kxk. El área superior izquierda tiene solo 𝑎0 en la diagonal, ceros
por arriba de la diagonal, y subíndices progresivamente mayores para los coeficientes sucesivos en

cada columna que está debajo de los elementos de la diagonal. Cuando transponemos los elementos

del área superior izquierda, obtenemos el área inferior derecha. Regresando al área superior

derecha, colocamos ahora el coeficiente 𝑎𝑛 solo en la diagonal, con ceros debajo de la diagonal, y

subíndices progresivamente más pequeños para los coeficientes sucesivos a medida que subimos

por cada colina desde la original. Cuando se transponen los elementos de esta área, obtenemos el

área inferior izquierda.

La aplicación de este teorema es directa. Dado que los coeficientes de la ecuación característica de

la ecuación son los mismos que los que aparecen en el lado izquierdo de la ecuación en diferencias

original, podemos introducirnos directamente en los determinantes citados.

Ejercicios donde se aplica el teorema de Schur:

Ejercicio 1

Observe que, en nuestro contexto, siempre tenemos 𝑎0 = 1

1 1
Pruebe la convergencia de la trayectoria de 𝑦𝑡+2 + 6 𝑦𝑡+1 − 6 𝑦𝑡 = 2 mediante el teorema de Schur.

1 1
Aquí los coeficientes son 𝑎1 = 6 , 𝑎2 = − 6 (con n=2). Entonces tenemos

1
1 −
𝑎0
∆1 = |𝑎
𝑎2
| = | 6| = 35 > 0
2 𝑎0 1 36
− 1
6
 1 1
1 0 −
6 6
𝑎2 𝑎1 | 1 |
𝑎0 0 1
1 0 −
𝑎
∆2 = |𝑎1
𝑎0 0 𝑎2
| = | 6 6| = 1176 > 0
2 0 𝑎0 𝑎1 1 1296
𝑎1 𝑎2 0 𝑎0 − 0 1
| 6 1 1 6 |
1 − 0
6 1
6

Estos realmente satisfacen la condición necesaria y suficiente para la convergencia.

Ejercicio 2

Prueba la convergencia de las soluciones de las siguientes ecuaciones en diferencias por el teorema

de Schur:

1
a) 𝑦𝑡+2 − 2 𝑦𝑡+1 − 𝑦𝑡 = 0

1
Los coeficientes son: 𝑎0 = 1, 𝑎1 = − 2 , 𝑎2 = −1

𝑎0 𝑎2 1 −1
∆1 = |𝑎 𝑎0 | = |−1 1 | = 0
2

1
1 0 −1 −
𝑎0 0 𝑎2 𝑎1 2
| 1 |
𝑎 𝑎0 0 𝑎2
∆2 = |𝑎1 | = −2 1 0 −1 = 1 > 0
2 0 𝑎0 𝑎1
𝑎1 𝑎2 0 𝑎0 | −1 0 1
1 − |
1 2
− −1 0
2 1

De este ejercicio se puede apreciar que la ecuación no satisface la condición necesaria ya que los

determinantes no son estrictamente positivos.

4. MODELO DE CICLO DE INVENTARIO DE METZLER

Metzler fue el primero en investigar formalmente las consecuencias de los esfuerzos de los

empresarios para mantener sus niveles de existencias, a través de variaciones apropiadas en la

producción. Aquí vamos a exponer las características básicas de los modelos que construyó con

esta finalidad para después pasar a examinar uno de sus modelos de segundo orden.

El producto total del periodo es la suma de los bienes de consumo y de inversión producidos, siendo

el volumen de estos últimos una constante exógena. El volumen total de los bienes de consumo se

compone de dos partes:

(1) La producción que normalmente debe venderse, de acuerdo con las expectativas de ventas de

los productores.

(2) La producción que se destina a mantener las existencias a su nivel deseado.

Este segundo componente puede, desde luego, ser negativo, lo cual vendría a significar que los

vendedores quieren producir menos de lo que esperan vender, cubriendo la diferencia mediante la

reducción deseada de existencias.

Desde luego, estas expectativas pueden no ser exactas, es decir, las ventas reales pueden ser

diferentes de las ventas esperadas, suponiendo esta diferencia una variación no deseada de

existencias. Obsérvese que las ventas reales coinciden con la demanda de consumo del periodo

que no debe ser confundida con la producción de bienes de consumo, ya que en el modelo de

Metzler la producción total (tal como arriba se define) y la demanda de bienes de consumo, pueden

ser diferentes.
A partir de estas características generales se puede obtener un gran número de modelos distintos

según los supuestos concretos que se formulen en torno a la formación de expectativas y a cómo

se determina el nivel deseado de existencias.

Vamos aquí a suponer que las expectativas son del tipo “ingenuo”, es decir, que las ventas

corrientes esperadas Ut son iguales a las ventas realizadas en el periodo precedente; en ese caso

Ut = Ct - 1 (1.1)

Donde Ut es por tanto el componente que arriba llamamos (1) y Ct-1 la demanda de consumo del

periodo anterior, aunque no se presume la existencia de ningún tipo de desfase en la función de

consumo, de modo que

Ct = bYt 0<b<1 (1.2)

Y por lo tanto

Ut = bYt – 1 (1.3)

Para especificar el componente (2), suponemos que los productores desean mantener una cierta

relación constante entre sus existencias y sus ventas: llamemos a esta relación k (que será el

“acelerador de las existencias”). Puesto que las ventas reales sólo serán conocidas ex post, los

productores aplicarán esta relación a las ventas esperadas para calcular el nivel deseado de

existencias, Qt, que deberá ser alcanzado en el periodo:

Qt = kUt k>0 (1.4)

El componente (2) – la inversión en existencias – viene dado pues por

Qt - Qt - 1 = kUt - Qt – 1 = k bYt – 1 - Qt – 1 (1.5)

Donde Qt – 1 es el nivel de existencias en el momento en que los productores hacen sus planes (al

comienzo del periodo t, es decir, al final de t -1). Ahora bien, Qt – 1 es el nivel de existencias que

los productores habían planeado para el periodo t -1, es decir kUt – 1 =k bYt – 2 menos

la variación involuntaria de existencias (si es que ésta tuvo lugar) en ese periodo entre las ventas

realizadas Ct – 1 = bYt – 1, y las ventas esperadas, Ut – 1= bYt – 2. Así pues, tenemos que

Qt – 1 = kbYt – 2 – b (Yt – 1 – Yt – 2 ) (1.6)

Como decíamos al principio, el producto (la renta) nacional viene dado por

Yt = Ut + (Qt - Qt – 1) + I0 (1.7)

donde I0 es la inversión autónoma.

Tras realizar una serie de sencillas sustituciones obtenemos la siguiente ecuación, no homogénea

de segundo orden

Yt – (2+k)bYt – 1 + (1+k)bYt – 2 = I0 (1.8)

Para ésta, se obtiene una solución particular probando con 𝑌̅𝑡 = constante, lo que da

1
𝑌̅𝑡 = 𝑡−𝑏 I0 (1.9)

Es decir, que el equilibrio estacionario viene dado como siempre por la aplicación del

multiplicador al gasto autónomo constante.

La ecuación homogénea correspondiente tiene por ecuación característica

λ2 – (2+k)bλ + (1+k)b = 0 (1.10)

Aplicando ahora las condiciones de estabilidad tenemos que, puesto la propensión al consumo es

menor que la unidad, la condición esencial es que

 1
𝑏 < 1+𝑘 (1.11)

La sucesión de signos de los coeficientes de la ecuación (5.10) es + - + de forma que no puede

haber raíces negativas. El discriminante de (1.10) es

Δ = (2+k)2b2 – 4(1+k)b ,

de forma que

> > 4(1+𝑘)

Δ<0 si b < (2+𝑘)2 (1.12)

Y puesto que la desigualdad

4(1+𝑘) 1
> 1+𝑘 (1.13)
(2+𝑘)2

Se cumple para cualquier k > 0, nos encontramos con los siguientes tipos de movimiento

(referidos a la situación de equilibrio estacionario):

1
0 < 𝑏 < 1+𝑘 , oscilante y convergente;

1
𝑏 = 1+𝑘 , oscilante y de amplitud constante;

1 4(1+𝑘)
< 𝑏 < (2+𝑘)2 , oscilante y divergente;
1+𝑘

4(1+𝑘)
𝑏 ≥ (2+𝑘)2 , monótono y divergente.
Es interesante observar que no puede existir un movimiento monótono estable; la estabilidad sólo

puede adoptar la forma de oscilaciones atenuadas. Obsérvese también que, ceteris paribus, cuando

menor sea la relación que los fabricantes desean mantener entre existencias y ventas, más probable

es que tenga lugar un movimiento estable. El siguiente gráfico nos muestra como siempre las

diferentes regiones que el lector puede tratar de interpretar, a modo de ejercicio

Gandolfo, G. (1976). Métodos y modelos matemáticos de la dinámica económica. [Figura 5.6]

Ejemplo:

1. En el periodo inicial a la renta se encuentra en equilibrio en el valor 2.000; la producción (y

las ventas) de bienes de consumo es de 1.000; las existencias se encuentran al nivel deseado y

la inversión autónoma es de 1.000. En el periodo 1 la inversión autónoma pasa a ser de 1.100

y se mantiene en este nivel en los periodos sucesivos. Dado que la propensión al consumo es

de 0,5 y que la relación deseada existencias/ ventas es de 0,2, calcúlense los valores de Yt para

un cierto número de periodos, y averígüese matemáticamente el tipo de movimiento que

presenta Yt.
Solución:

Los cálculos iterativos se presentan en el cuadro. Los encabezamientos de las columnas no

requieren ninguna aclaración. Los valores de las columnas señaladas con un asterisco se suman

para formar Yt; el resto de las columnas se presentan para que se puedan seguir los cálculos.

En el periodo 0 el sistema se encuentra en equilibrio; en el periodo 1 la inversión autónoma pasa

a ser de 1.100; la producción de bienes de consumo sigue en 1.000, la renta es de 2.100 y el

consumo es de 1.050, mayor por tanto que el consumo esperado (la producción de bienes de

consumo para la venta, Ut): la diferencia origina una reducción de existencias.

En el periodo 2 la producción de bienes de consumo es de 1.050 para hacer frente al consumo

esperado más 60 para llevas las existencias a su nivel deseado. La renta pasa así a ser de 2.210, el

consumo de 1.105, de nuevo superior a su nivel esperado, dando lugar a nuevas reducciones de

existencias, etc. Los valores de Yt se han calculado hasta t=17 aplicando directamente la ecuación

Yt =1,1Yt-1 – 0,6Yt-2 + 1.100

Datos:

 b=0.5 𝑄̂ t=kUt Ct = bYt-1

 k=0.2 Ut = Ct-1
t I0(*) Ut (*) 𝑄̂ 𝑡 𝑄̂ 𝑡 – Qt-1(*) Yt Ct Ct-Ut Qt
0 1000 1000 200 0 2000 1000 0 200
1 1100 1000 200 0 2100 1050 50 150
2 1100 1050 210 60 2210 1105 55 155
3 1100 1105 221 66 2271 1135.5 30.5 190.5
4 1100 1135.5 227.1 36.6 2272.1 1136.05 0.55 226.55
5 1100 1136.05 227.21 0.66 2236.71 1118.355 -17.695 244.905
6 1100 1118.355 223.671 -21.234 2197.12 1098.56 -19.795 243.466
7 1100 1098.56 219.712 -23.754 2174.81 1087.405 -11.155 230.867
8 1100 1087.405 217.481 -13.386 2147.02 1073.51 -13.895 231.376
9 1100 1073.51 214.702 -16.674 2186.54 1093.27 19.76 194.942
10 1100 1093.27 218.654 23.712 2200.78 1100.39 7.12 211.534
11 1100 1100.39 220.078 8.544 2208.93 1104.465 4.075 216.003
12 1100 1104.465 220.893 4.89 2209.36 1104.68 0.215 220.678
13 1100 1104.68 220.936 0.258 2204.94 1102.47 -2.21 223.146
14 1100 1102.47 220.494 -2.652 2199.82 1099.91 -2.56 223.054
15 1100 1099.91 219.982 -3.072 2196.84 1098.42 -1.49 221.472
16 1100 1098.42 219.684 -1.788 2196.63 1098.315 -0.105 219.789
17 1100 1098.315 219.663 -0.126 2198.19 1099.095 0.78 218.883

Determinación del tipo de movimiento:

1
0 < 0.5< 0 < 0.5 < 0.83
1+0.2

Entonces nos encontramos ante un movimiento oscilante y convergente, es decir, que la

producción total en los primeros períodos presenta fuerte oscilaciones que con el transcurso del

tiempo irán disminuyendo convergiendo así a una cantidad que en nuestro caso particular estará

próximo a 2.200
Veamos ahora qué pasa con el movimiento; evidentemente es oscilante y atenuado en torno al

valor de equilibrio Yt = [1/ (1 – 0,5)] X 1.100 = 2.200. A modo de comprobación puede verse

con facilidad que las raíces de la ecuación característica

Conclusión:

No puede existir un movimiento monótono estable; la estabilidad solo puede adoptar la forma, de

oscilaciones atenuadas. Obsérvese también que cuanto menos sea la relación que los fabricantes

desean mantener entre existencias y ventas, más probable es que tenga un lugar movimiento

estable.

La solución de equilibrio de inventario de Metzler está dado por:

1
𝑌̅𝑡 = 𝑡−𝑏 I0

.
 II. CONCLUSIONES

 Para los tres casos de las raíces características, la trayectoria de tiempo va a convergir en el
equilibrio intertemporal, si y sólo si el valor absoluto de cada una de las raíces es menor de
1, independientemente de cuáles puedan ser las condiciones iniciales.
 No puede existir un movimiento monótono estable; la estabilidad solo puede adoptar la
forma, de oscilaciones atenuadas. Obsérvese también que cuanto menos sea la relación que
los fabricantes desean mantener entre existencias y ventas, más probable es que tenga un
lugar movimiento estable.

III. BIBLIOGRAFIA

 Chiang, A. (2006) Métodos Fundamentales de Economía Matemática, cuarta edición

(p.568, 583, 584, 585)

 Gandolfo, G (1976) Métodos y modelos matemáticos de la dinámica económica (p.85 – 90)