Ecuaciones Diferenciales y en Diferencias

Economı́a Matemática: Ecuaciones Diferenciales y en
Diferencias
Ivan D’Agostino
Universidad Argentina de la Empresa - 2024

Ivan D’Agostino Economı́a Matemática: Ecuaciones Diferenciales 1y/en
24Diferencias
Introducción
Para esta sección, recomiendo los capı́tulos 14, 15.1 a 15.3 y 16 del
libro de Alpha Chiang o los capı́tulos 2 (sin incluir 2.3), 3.3 y 6 del
libro de Lomelı́ y Rumbos.

24Diferencias
Introducción
En el tema anterior, verificamos que podemos realizar estática

comparada sin necesidad de resolver el modelo.
Al estudiar las derivadas parciales, entendemos cómo cambia el
equilibrio cuando se modifica alguna variable exógena del modelo.
Sin embargo, nada sabemos sobre qué sucede entre un equilibrio y
otro, cómo se llega al nuevo equilibrio, si es uno estable o inestable,
etc.
Todos estos temas se relacionan con la parte dinámica de los
modelos, y no con la estática.

24Diferencias
Introducción
Algunas consideraciones importantes
En esta unidad, comenzaremos a estudiar variables que pueden
depender del tiempo. Si este último pertenece a los reales positivos, las
llamaremos ecuaciones diferenciales. En caso de que el tiempo
pertenezca a los naturales, trabajaremos con ecuaciones en
diferencias.
Pero, para estar consideradas dentro de estos grupos, no solo pueden
depender del tiempo, si no de la tasa de cambio respecto del
tiempo.
Algunas caracterı́sticas importantes:
Si la ecuación no involucra al tiempo (t), la llamaremos autónoma.
Definiremos el orden de una ecuación dependiendo de la derivada o
variable rezagada de mayor orden que contenga.
Si no hay un término independiente (constante), diremos que se trata
de una ecuación homogénea.
Si tanto la variable como sus derivadas o rezagos aparecen de forma
lineal, la ecuación será lineal. Esto no aplica para el tiempo.
Por último, si solo aparecen derivadas o rezagos respecto de la
misma variable, la llamaremos ordinaria.
24Diferencias
Introducción
Algunos ejemplos
ẋ = 5x + t
ẋ = 2x
ẋ = t 2 x
ẍ − 3ẋ + 2x − 2 = 0
xt = 2xt−1 + 2
xt + 2xt−1 − 3xt−2 = 0
xt2 + 3xt+1 = 10
24Diferencias
Ecuaciones Diferenciales
Una ecuación diferencial de orden ”n”, lineal, no autónoma, ordinaria
es una tal que tiene la forma:
F (an x n , an−1 x (n−1) , ..., a2 ẍ, a1 ẋ, a0 x, t) = h(t)
Donde las variables con un punto arriba representan las derivadas

parciales respecto del tiempo (y la cantidad de puntos el orden de la
derivada). Los coeficientes además son constantes.
En el curso, trabajaremos con ecuaciones diferenciales ordinarias,
lineales, no homogéneas y autónomas de primer orden.
Las ecuaciones con las que trabajaremos en el curso tienen una solución
que es la suma, a su vez, de una solución general y una particular:
x = xh + xp
La solución general responde a la parte homogénea de la ecuación (es

decir, la que obtendrı́amos si la ecuación diferencial estuviese igualada
a cero) y la particular está asociada a h(t).
24Diferencias
Solución
Supongamos que tenemos la siguiente ecuación diferencial:
aẋ = cx + b
Donde a, c y b son constantes. Si b=0, la ecuación serı́a homogénea.

Como vimos en el pizarrón, la solución será:
x(t) = Ae λt + B
Donde dijimos que λ = ca , B = − bc y que, para determinar A,

necesitamos un valor inicial para x.
Dado que trabajaremos siempre con el mismo tipo de ecuaciones
diferenciales, la solución que encontramos recién se va a generalizar a
todas las que tengamos que resolver. Ası́, usaremos el método de
coeficientes indeterminados, donde tendremos que hallar el valor de
los coeficientes λ y B.
24Diferencias
Solución por Coeficientes Indeterminados

Veamos como se resuelve una ecuación con este método. Tenemos la
siguiente ecuación diferencial:
3ẋ = 9x − 5
Proponemos como solución:
x(t) = Ae λt + B
Entonces, sustituimos en la ecuación a resolver:
3(Aλe λt ) = 9(Ae λt + B) − 5
Ahora, juntamos lo que contiene el exponencial de un lado y las

constantes del otro:
3Aλe λt − 9Ae λt = 9B − 5

24Diferencias
A continuación, factorizamos del lado izquierdo:
Ae λt (3λ − 9) = 9B − 5
Seguidamente, supondremos que tanto lo que está dentro del
paréntesis del lado izquierdo del igual como todo el término de la
derecha serán iguales a 0. Eso nos garantizará encontrar un valor para
los coeficientes tal que esta igualdad se sigue manteniendo. La
utilización del cero es conveniente.
Comencemos por el lado izquierdo:
3λ − 9 = 0
λ=3
Ahora, vamos por el lado derecho:
9B − 5 = 0
5
B=
9
24Diferencias
Ası́, llegamos a la siguiente solución:
5
x(t) = Ae 3t +
9
Como dijimos anteriormente, para encontrar un valor para A,
necesitamos una condición inicial. Supongamos que x(0) = 10:
5
10 = Ae 3×0 +
9
5
10 = A +
9
5
10 − = A
9
Ası́, la solución finalmente será
5 5
x(t) = (10 − )e 3t +
9 9
Ivan D’Agostino Economı́a Matemática: Ecuaciones Diferenciales10y/en
24Diferencias
Diagramas de Fase y Estado y el estudio de la Estabilidad

En casos más complejos, muchas veces no podremos encontrar una
solución explı́cita como lo hicimos en el ejemplo anterior.
En esos casos, será útil construir el Diagrama de Fase y el Diagrama
de Estado de las ecuaciones. Para ello, primero debemos encontrar los
llamados puntos fijos o estados estacionarios del sistema.
Dada la ecuación:
ẋ = f (x)
Se llama punto fijo a un valor x ∗ que satisface f (x ∗ ) = 0
Para encontrar los puntos fijos, impondremos entonces ẋ = 0.

24Diferencias
Diagramas de Fase, Estado y el Estudio de la Estabilidad

Supongamos que tenemos la siguiente ecuación diferencial:
ẋ = −x + 2
Imponiendo ẋ = 0, encontramos el punto fijo 2:
0 = −x + 2
x =2
Notemos que siempre que elijamos valores mayores a 2, ẋ será

negativa. Eso implica que x está decreciendo. Lo contrario sucede para
valores menores a 2, donde x crece.
Ya estamos en condiciones de construir el diagrama de fase, que
ubica a ẋ en el eje ”y” y a x en el eje ”x”.
Con él, podemos entender la dinámica implı́cita en la ecuación.

24Diferencias
Figure: Diagrama de fase de la ecuación ẋ = −x + 2

24Diferencias
La otra pregunta que surge respecto de la solución es su convergencia
¿Qué sucede a medida que pasa el tiempo? ¿La ecuación converge o
diverge a infinito?
Una vez que resolvimos la ecuación, también podemos responder esta
pregunta.
Continuando con la última ecuación, si la resolvemos llegaremos a:
x(t) = Ae −t + 2
El estudio de la convergencia consiste en tomar el lı́mite de la

ecuación cuando t tiende a infinito:
lim Ae −t + 2
t→∞
−∞
lim Ae +2=2
t→∞
Notemos que este último lı́mite converge ya que tiende a 2.

24Diferencias

Del último lı́mite, deducimos lo siguiente:
Para que una ecuación diferencial converja, λ debe ser menor a 0.
De lo contrario (λ mayor a 0), el lı́mite diverge.
Si λ es igual a 0, nos mantenemos en el mismo lugar.
Esto lo podremos ver construyendo el diagrama de estado.

24Diferencias
Figure: Diagramas de Estado de ẋ = 10e λt + 2.

24Diferencias
Ecuaciones en Diferencias
Cuando el tiempo posee un dominio discreto (en los números
naturales), estaremos trabajando con ecuaciones en diferencias.
Una ecuación en diferencias de orden n tiene la forma:
F (an−1 xt+n−1 , an−2 xt+n−2 , ..., a1 xt−1 , a0 xt , t) = h(t)
Aquı́ también trabajaremos con ecuaciones en diferencias lineales,

autónomas, ordinarias, no homogéneas y de orden 1.
En estos casos, la solución también vendrá como la suma de una
solución homogénea y una solución particular:
x = xh + xp

24Diferencias
Supongamos que debemos resolver la siguiente ecuación en diferencias:
xt = xt−1 (1 + r ) + c
Para resolver la parte homogénea, suponemos c=0 e iteramos dándole
varios valores a t, hasta encontrar un patrón:
x1 = x0 (1 + r )
x2 = x1 (1 + r )
x2 = [x0 (1 + r )](1 + r )
x2 = x0 (1 + r )2
Notemos que, de seguir iterando, llegaremos a la solución homogénea:
xt = x0 (1 + r )t
Para la solución particular, propondremos una constante:
xt = a
24Diferencias
Para la solución particular, propondremos una constante:
xt = a
Usemos la propuesta en nuestra ecuación en diferencias:
a = a(1 + r ) + c
a = a + ar + c
−ar = c
c
a=−
r
Ası́, la solución general será:
c
xt = x0 (1 + r )t −
r

24Diferencias
Al igual que con las ecuaciones diferenciales, como trabajaremos con el
mismo tipo de ecuaciones en diferencias durante todo el curso,
podremos proponer una solución y resolver vı́a el método de
coeficientes indeterminados.
Notemos que la solución a proponer tendrá siempre la forma:
xt = Aλt + B
Para notar que esta solución propuesta es correcta, veamos si podemos

llegar a la misma solución que encontramos vı́a iteración:
Aλt + B = (Aλt−1 + B)(1 + r ) + c

Aλt − (Aλt−1 )(1 + r ) = B(1 + r ) − B + c
1
Aλt (1 − (1 + r )) = B + Br − B + c
λ

24Diferencias
Del lado izquierdo, nos queda:
1
1− (1 + r ) = 0
λ
1
(1 + r ) = 1
λ
1 1 −1
( )−1 = ( )
λ 1+r
λ=1+r
Del lado derecho, nos queda:
B + Br − B + c = 0
Br = −c
c
B=−
r
Si imponemos A = x0 ,la solución a la que llegamos es:
c
xt = x0 (1 + r )t −
r
24Diferencias
Al igual que para las ecuaciones diferenciales, aquı́ también podemos
estudiar la convergencia de la solución. Esto también lo haremos
tomando el lı́mite de la solución cuando t tiende a ∞.
Las reglas de convergencia, sin embargo, cambian levemente:
Si |λ| es menor a 1, la solución converge. De ser negativo, converge
haciendo una especie de zigzag.
Si |λ| es mayor o igual a 1, la solución diverge.
Aquı́ también podemos construir diagramas de fase y estado.

24Diferencias
Figure: Diagramas de Estado de xt = 10λt + 2.

24Diferencias
Ecuaciones Diferenciales y en Diferencias
Visto esto, vamos a la práctica.

24Diferencias

Ecuaciones Diferenciales y en Diferencias

Cargado por

Copyright:

Formatos disponibles

Ecuaciones Diferenciales y en Diferencias

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ecuaciones Diferenciales y en Diferencias

Cargado por

Copyright:

Formatos disponibles

Economı́a Matemática: Ecuaciones Diferenciales y en

Universidad Argentina de la Empresa - 2024

Universidad Argentina de la Empresa - 2024

Universidad Argentina de la Empresa - 2024

En el tema anterior, verificamos que podemos realizar estática

Universidad Argentina de la Empresa - 2024

F (an x n , an−1 x (n−1) , ..., a2 ẍ, a1 ẋ, a0 x, t) = h(t)

Donde las variables con un punto arriba representan las derivadas

La solución general responde a la parte homogénea de la ecuación (es

Donde a, c y b son constantes. Si b=0, la ecuación serı́a homogénea.

Donde dijimos que λ = ca , B = − bc y que, para determinar A,

Solución por Coeficientes Indeterminados

Proponemos como solución:

Entonces, sustituimos en la ecuación a resolver:

Ahora, juntamos lo que contiene el exponencial de un lado y las

Universidad Argentina de la Empresa - 2024

Diagramas de Fase y Estado y el estudio de la Estabilidad

Universidad Argentina de la Empresa - 2024

Diagramas de Fase, Estado y el Estudio de la Estabilidad

Imponiendo ẋ = 0, encontramos el punto fijo 2:

Notemos que siempre que elijamos valores mayores a 2, ẋ será

Universidad Argentina de la Empresa - 2024

Figure: Diagrama de fase de la ecuación ẋ = −x + 2

El estudio de la convergencia consiste en tomar el lı́mite de la

Notemos que este último lı́mite converge ya que tiende a 2.

Diagramas de Fase, Estado y el Estudio de la Estabilidad

Universidad Argentina de la Empresa - 2024

Diagramas de Fase, Estado y el Estudio de la Estabilidad

Figure: Diagramas de Estado de ẋ = 10e λt + 2.

Universidad Argentina de la Empresa - 2024

F (an−1 xt+n−1 , an−2 xt+n−2 , ..., a1 xt−1 , a0 xt , t) = h(t)

Aquı́ también trabajaremos con ecuaciones en diferencias lineales,

Universidad Argentina de la Empresa - 2024

Usemos la propuesta en nuestra ecuación en diferencias:

Universidad Argentina de la Empresa - 2024

Para notar que esta solución propuesta es correcta, veamos si podemos

Aλt + B = (Aλt−1 + B)(1 + r ) + c

Universidad Argentina de la Empresa - 2024

Universidad Argentina de la Empresa - 2024

Diagramas de Fase, Estado y el Estudio de la Estabilidad

Figure: Diagramas de Estado de xt = 10λt + 2.

Universidad Argentina de la Empresa - 2024

Visto esto, vamos a la práctica.

Universidad Argentina de la Empresa - 2024

También podría gustarte