Semana 5:
Tiempo discreto: Ecuaciones en Diferencias de
segundo orden
Chiang, Alpha. 1996. Métodos Fundamentales Economía
Matemática -3ra. Edición. MC Graw Hill.
Aldo Ramírez
aframire@[Link]
2019
ECUACIÓN EN DIFERENCIAS DE SEGUNDO ORDEN
Sabemos que Δyt, es:
yt yt 1 yt
Ahora definamos una segunda diferencia para “y”:
• Δ2yt = Δ(Δyt)= Δ(yt+1 -yt)
• Δ2yt = (yt+2 –yt+1) - (yt+1 -yt)
• Δ2yt = yt+2 - 2yt+1+yt
Así, una diferencia segunda de yt puede transformarse en una
suma de términos que incluyen un retraso de dos períodos.
Ecuación en diferencias de segundo orden
• Consideremos la ecuación en diferencias de segundo orden.
𝑦 𝑡 +2 +𝑎 1 𝑦 𝑡 +1+ 𝑎 2 𝑦 𝑡 =𝑐
Donde “a1”, “a2” y “c” son constantes; cuando “c” es cero
tenemos el caso homogéneo.
Recordemos que si yc es la ecuación complementaria e yp la
solución particular, entonces:
yt y c y p
Es la solución general de la ecuación completa.
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
La Integral particular
En el caso de coeficientes y término constantes, el caso más sencillo es asumir
yt=k, con lo que yt+1=yt+2=k, entonces:
k+a1k+a2k=c
k= c / (1+a1+a2)
yp=k= c / (1+a1+a2)
Donde debe cumplirse que (1+a1+a2) sea diferente de cero
Ejemplo: Halle la yp de la siguiente ED
yt 2 3 yt 1 4 yt 6
Si yt=k, entonces yt+1=yt+2=k, con lo que:
k-3k+4k=6
k=3
yp=3
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
La Integral particular
Ejemplo:
yt 2 yt 1 2 yt 12
Si hacemos yt=k…
No tendremos solución, entonces…
Habrá que suponer yt=kt, entonces yt+1=k(t+1), yt+2=k(t+2) con lo
que:
k(t+2) + k(t+1) - 2kt =12
k=4, entonces
yp=kt=4t
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
La Función Complementaria
La FC es la solución general de la ecuación en forma reducida (u homogénea):
a1 yt 1 a2 yt
Sabemos por la semanayt4 2que: 0
yt=Abt
Por lo tanto:
yt+1=Abt+1 , yt+2=Abt+2
La ecuación se transforma en:
Abt 2 a1 Abt 1 a 2 Abt 0 o también:
Abt b 2 a1b a 2 0
Si escogemos los valores de “b” que satisfacen la ecuación, la conjetura que
dice que y=Abt funcionará adecuadamente.
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
La Función Complementaria
Para que la última ecuación resuelva la ED, debemos hacer que A=0 o
encontrar las raíces de:
b 2 a1b a 2 0
Para esto nos sirve la fórmula cuadrática, que nos dará (en este caso) dos
soluciones (raíces) posibles:
2
a1 a1 4a2
b1 , b2
2
Donde debe cumplirse que:
b1+b2=-a1, y también b1.b2=a2
Tenemos dos soluciones como parte de la yc:
t t
yc y1 y2 A1b1 A2b2
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
Caso 1: Raíces reales distintas
Cuando:
a12>4a2
Entonces b1 y b2 son dos números reales diferentes y la FC se mantiene como:
t t
yc y1 y2 A1b1 A2b2
Ejemplo, resolver (cuaderno), con las condiciones iniciales y0=4 e y1=5:
yt 2 yt 1 2 yt 12
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
Caso 2: Raíces reales repetidas
Cuando:
a12=4a2
Entonces b1 y b2 son dos números reales iguales y la FC puede ser escrita
como:
yc y1 y2 A1b1 A2b2 A1 A2 b t A3b t
t t
Sin embargo, con sólo una raíz real necesitamos otra expresión linealmente
independiente (li) de A3bt para poder encontrar la función primitiva.
Usemos A4tbt que es li para poder encontrar la sg. (“t” es la clave)
Entonces:
yc A3b t A4tb t
Ejemplo: resolver
yt 2 6 yt 1 9 yt 4
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
Raíces complejas
Cuando a12<4a2 , entonces con la raíz:
a12 4a2
Sabemos que podemos hacer:
1 4a2 a12 i 4a2 a12
Y ahora:
a1 4a2 a12
h v
2 2
Entonces las dos raíces pueden ser expresadas como un par de números
complejos conjugados:
b1b2 h vi
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
Análisis de Raíces Complejas
Con:
v h
sen cos
R R
Y con:
R a2 h vi t R t Cost iSent
Tenemos que:
yc A1 h vi A2 h vi
t t
yc A1 R t cos t isen t A2 R t cos t isen t
yc R t A1 A2 cos t A1 A2 isen t
yc R t A5 cos t A6 sent
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO ORDEN CON
COEFICIENTES Y TÉRMINO CONSTANTES
Análisis de Raíces Complejas
Recordemos:
Resolver:
yt 2 4 yt 1 16 yt 0
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
Convergencia de la trayectoria temporal
Igual que en las ED de primer orden, la convergencia depende de que y c tienda a
cero cuando t tiende a infinito. Las propiedades de las raíces “b” se mantienen.
Con raíces reales diferentes, si ambas raíces “b” cumplen que:
b1 1 b2 1
Entonces, la solución converge.
Si una de las raíces es mayor a la otra en valor absoluto, la convergencia depende
de que la mayor de ellas en valor absoluto sea menor a uno. A la raíz mayor en
valor absoluto se le llama “Raíz dominante”.
Una Trayectoria temporal será convergente sí y solo sí la raíz dominante es menor
a uno en valor absoluto.
ECUACIONES DIFERENCIALES LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES Y TÉRMINO CONSTANTES
Convergencia de la trayectoria temporal
Con raíces reales repetidas, tenemos dos términos en la yc, que son:
yc A3b t A4tb t
Entonces, la solución converge si:
b 1
La Trayectoria temporal será convergente sí y solo sí la raíz es menor a uno en
valor absoluto.
A pesar de que “t” va creciendo, el efecto combinado con un b en valor absoluto
menor a uno es decreciente.
Ejemplo, analizar la convergencia de:
t
yt 3 2 4t
Modelo de aplicación
• MODELO DE INTERACCIÓN DE MULTIPLICADOR CON
ACELERADOR DE SAMUELSON
{
𝑌 𝑡 = 𝐶𝑡 + 𝐼 𝑡 + 𝐺0
𝐶 𝑡= 𝛾 𝑌 𝑡 −1 , 0<𝛾 <1
𝐼 𝑡 = 𝛼 ( 𝐶 𝑡 − 𝐶 𝑡 −1 ) ,𝛼 >0
INTERACCIÓN DE MULTIPLICADOR CON ACELERADOR DE
SAMUELSON
ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES CONSTANTES Y TÉRMINO
VARIABLE
Ecuaciones en diferencias de segundo orden con término variable
Veamos ahora:
yt 2 a1 yt 1 a2 yt cm t
Donde “m” es una función de “t”.
De nuevo, nuestra función complementaria yc no se altera pues es
la ecuación cuando cmt es cero.
Pero yp debe calcularse de manera diferente.
Método de Coeficientes indeterminados
Primero veamos el caso en que “c” y “m” son constantes.
ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES CONSTANTES Y TÉRMINO
VARIABLE
Método de Coeficientes indeterminados
Veamos el siguiente ejemplo:
yt 2 yt 1 3 yt 7 t
Aquí tenemos c=1 y m=7. Debemos asegurarnos que el término variable 7 t sea
susceptible de diferenciación, así su primera diferencia será:
7 7
t t 1
7 t
7 1 7 6 7
t
t
Su segunda diferencia:
7 t 6 7 t 67 67 67 17 t 367
t 1 t t
Todas las diferencias sucesivas estarán en función de 7 t, por lo tanto la solución
probable a la integral particular podrá ser yt=B(7)t. Por lo tanto usamos esta yt y sus
versiones para t+1 y t+2:
B7 B7 3B7 7 t
t 2 t 1 t
ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES CONSTANTES Y TÉRMINO
VARIABLE
Método de Coeficientes indeterminados
Sigamos con:
B7 B7 3B7 7 t
t 2 t 1 t
Entonces:
B 7 2 7 3 7 7 t
t
Con lo que obtenemos B:
B=1/(49+7-3)=1/53
yp=B(7)t=(1/53) (7)t
Este es un equilibrio móvil.
ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES CONSTANTES Y TÉRMINO
VARIABLE
Método de Coeficientes indeterminados
Veamos ahora el caso ctn, uno tal como:
yt 2 5 yt 1 2 yt t 2
Tomemos las tres primeras diferencias de t2:
t 2 t 1 t 2 2t 1
2
2t 2 t 2 2t 1 2t 1 2t 2(t 1 t ) 2
3t 2 2t 2 2 0
Como a partir de la tercera diferencia el resultado es cero, la forma funcional
buscada parece ser un polinomio de grado 2, tal como:
yt B0 B1t B2t 2
ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES CONSTANTES Y TÉRMINO
VARIABLE
Método de Coeficientes indeterminados
Sigamos con:
yt B0 B1t B2t 2
Entonces:
yt 1 B0 B1 (t 1) B2 (t 1) 2 B0 B1 B2 B1 2 B2 t B2t 2
yt 2 B0 B1 (t 2) B2 (t 2) 2 B0 2 B1 4 B2 B1 4 B2 t B2t 2
Sustituyendo en
yt 2 5 yt 1 2 yt t 2
Obtenemos finalmente:
8B0 7 B1 9 B2 8B1 14 B2 t 8B2t 2 t 2
ECUACIONES EN DIFERENCIAS LINEALES DE SEGUNDO
ORDEN CON COEFICIENTES CONSTANTES Y TÉRMINO
VARIABLE
Método de Coeficientes indeterminados
Igualando los términos tenemos:
8B0 7 B1 9 B2 0
8B1 14 B2 0
8 B2 1
Hallamos los valores:
B0=13/256, B1=-7/32 y B2=1/8
Obtenemos finalmente la yp:
y p B0 B1t B2t 2
13 7 1
yp t t2
256 32 8