Ecuaciones Diferenciales Ordinarias
Ecuaciones Diferenciales Ordinarias
x = f (r , w , p , K , t ) (I )
Sin embargo, nosotros no estudiaremos este tipo de ecuaciones diferenciales en
este manual ya que escapa a nuestros objetivos.
Ejemplos:
∂x ∂x
1. − +5 = 4 ⇒ x = f (w , r )
∂w ∂r
3 2
∂2x 3 4
2. − + ln r ∂ x + senz ∂x = 0 ⇒ x = f (w , r , z )
∂w 2 ∂r 3 ∂z
Ejemplos:
dy dz dy dz
1. − F y , , z, ; r = ln y + senz − 10 r = 0 ≡ ln ydy + senzdz = 10 rdr
dr dr dr dr
dy dw dy dw
2. − F y , , w, ; r = y +w − r 2 = 0 ≡ ydy + wdw = r 2 dr
dr dr dr dr
No obstante, en este manual sólo abordaremos las EDO que contienen derivadas
totales de una única variable dependiente con respecto a una sola variable
independiente, y que pueden expresarse en forma implícita:
dx d 2 x d (n x
F x , , ,K, ,r = 0 (III)
dr dr 2 dr n
Donde:
F : Ψ ⊆ ℜ n +2 → ℜ
dx (n d (n −1 x d 2 x dx
= f , KK , , , x, r (IV )
dr n dr n −1 dr 2 dr
Dónde:
f : Φ ⊆ ℜ n +1 → ℜ
Es decir, nosotros estudiaremos aquellas EDOs que describan una relación entre una
función desconocida “x” y sus derivadas/diferenciales totales. A la ecuación (IV) se
le denomina ecuación diferencial ordinaria normal. La solución de este tipo de EDOs
es una función desconocida “ x (r ) ” definida en un dominio “D” que admite
derivadas totales hasta de orden “n”, también definidas en “D”, tales que esta función
y sus derivadas satisfacen (IV) como una identidad cuando r ∈ D. Por tanto, la
solución general de (IV) será:
x:D⊆ℜ→ ℜ
(V)
r → x (r )
238
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Ejemplos:
dx dx
1. − F x , , r = 4 r + 6 x − 2 = 0 ≡ 4 rdx + (6 x − 2 )dr = 0 ⇒ x = x (r )
dr dr
dx d 2 x d 2 x 3
2. − F x , ,
dr dr 2
;r =
dx
+ 5 x ( )
+ 8 2 + r 2 x = 0 ⇒ x = x (r )
dr 2 dr
3
dx d 2 x d 3 x dx
4
d2x d3x
3. − F x , , , ; s = e x +4
+ ln s = 0 ⇒ x = x (s )
ds ds 2 ds 3 ds 2 ds 3
ds
dx d 2 x d ( n x
F x , , ,K, ,t = 0 (VI )
dt dt 2 dt n
Dónde:
F : Ψ ⊆ ℜ n+2 → ℜ
dx (n d (n −1 x d 2 x dx
= f , KK , , , x, t (VII )
dt n dt n −1 dt 2
dt
Dónde:
f : Φ ⊆ ℜ n +1 → ℜ
x : D ⊆ ℜ+ → ℜ
(VIII )
t → x (t )
La solución general de una EDO contendrá “n” constantes arbitrarias, por lo que
una EDO tendrá infinitas soluciones.
239
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
V.3 Existencia y unicidad de una solución
A la ecuación (VII) junto con las siguientes “n” condiciones iniciales:
x (t 0 ) = x 0
x ' (t 0 ) = x '0
(IX )
M
x (n −1 (t 0 ) = x 0( n −1
Ejemplos:
1.- Determinar la solución del siguiente problema de valor inicial:
x '' ( t ) = 10
x (0) = 5 (1)
'
x (0 ) = 4
dx ' (t )
∫ x '' (t )dt = 10dt ⇒
∫ ∫ dt = 10dt ⇒ x ' (t ) = 10 t + A
∫ (2 )
dt
1
Si junto a la ecuación (VII) se dan condiciones (denominadas condiciones de frontera) que corresponden
a valores distintos de la variable independiente se dice que se tiene un problema de valores en la frontera.
240
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Integrando (2) respecto de “t” resulta:
dx (t )
∫ x (t )dt = ∫ (10 t + A )dt ⇒ ∫ ∫ (10 t + A )dt
'
dt =
dt
x (t ) = 5 t 2 + At + B (3)
x ' (0 ) = A = 4
(4 )
x (0 ) = B = 5
x (t ) = 5 t 2 + 4 t + 5
Figura 1
241
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
2.- Determinar la solución del siguiente problema de valor inicial:
x '''(t) = 120
x(0) = 10
' (5)
x (0) = 4
''
x (0) = 8
dx ''(t)
∫ x '''(t)dt = 120dt ⇒
∫ ∫ dt = 120dt ⇒ x ''(t) = 120t + A
∫ (6)
dt
dx '(t)
∫ ∫ (120t + A)dt ⇒ x (t) = 60t (7)
' 2
dt = + At + B
dt
dx(t)
∫ x (t)dt = ∫ (60t ) ∫ ∫ (60t )
' 2 2
+ At + B dt ⇒ dt = + At + B dt
dt
A
x (t ) = 20 t 3 + t 2 + Bt + C (8 )
2
x '' (0 ) = A = 8
x (0 ) = B = 4 (9 )
'
x (0 ) = C = 8
x (t ) = 20 t 3 + 4 t 2 + 4 t + 8
242
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
En la figura 2 se han representado algunas de las curvas que forman parte de
esta familia de soluciones de x '''(t) = 120. Entre ellas, se ha representado la
solución particular que resuelve el problema de valores iniciales (5).
Figura 2
Ejemplos:
3
d2x
+ 3x 3 t 4 − 7 = 0 : EDO de segundo orden y tercer grado.
dt 2
2
= e (5+ 3 t )
dx 2
5
: EDO de primer orden y segundo grado.
dt
2
∂ 3u
∂r 3
( )
− u 2 (r )3 − sen r 3 u 2 + 5 = 0 : EDP de tercer orden y segundo grado.
243
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Clasificación de las ecuaciones diferenciales según el método de solución
Ecuaciones diferenciales
Ordinarias Parciales
Lineales No lineales
Homogéneas No homogéneas
244
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
1. EDOs lineales de primer orden con coeficientes constantes
La forma general de este tipo de ecuaciones puede expresarse como:
a0 f (t)
a1x ' + a 0 x = f (t) ⇒ Si (a1 ≠ 0) ⇒ x ' + x=
a1 a1
Resultando que:
( ) = [x e
d xe wt ' wt
]
+ wxe wt = q(t)e wt
dt
( )
d xe wt
∫ dt = ∫ q(t)e ∫ q(t)e
wt
dt ⇒ xe wt = wt
dt + C
dt
x (t ) = e − wt [∫ q(t )e wt
dt + C ] (XIII)
Para el caso homogéneo, q(t) = 0, entonces la versión homogénea de (XII)
sería:
x (t ) = Ce − wt (XV )
q q
x (t ) = e − wt e wt + C = + Ce − wt (XVII )
w w
245
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Para el caso homogéneo, q = 0, entonces la versión homogénea de (XVI) sería:
x (t ) = x c (t ) = Ce − wt (XIX )
A la solución de la ecuación diferencial homogénea (XVIII), que viene dada
por (XIX), se le suele denominar solución complementaria: x c (t ).
( )
x ' = g x p = g (x E ) = q − wx p = 0 ⇒ x p = x E =
q
(XX )
w
q
Donde a x p = x E = se le suele denominar valor de equilibrio, punto fijo o
w
solución particular de (XVI).
q
x (t ) = + Ce − wt = x p + x c (t ) (XXI )
w
( )
x (t ) = x p + x 0 − x p e − wt = x p + x c (0 ) ⋅ e − wt = x p + x c (t ) (XXIII)
La ecuación (XXIII) nos permite interpretar a la solución particular
x p = x E = q w como el equilibrio fijo (ya que no depende del tiempo) de
(
x (t ) y a la solución complementaria x c (t ) = x 0 − x p e − wt como el desvío )
de la trayectoria x (t ) de su nivel de equilibrio a lo largo del tiempo. Este
desvío crecerá exponencialmente en el tiempo a una tasa “w” si w < 0. Esto
es, si w < 0, conforme t → +∞ ⇒ x c (t ) = Ce − wt → +∞, por lo que en el
largo plazo (LP), si t → + ∞ ⇒ x (t ) → +∞. En términos formales:
x LP = lím x (t ) → +∞ (XXIV )
t → +∞
246
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Este caso, con w < 0 se dice que es inestable ya que el lím x (t ) → +∞.
t → +∞
q
x LP = lím x (t ) = = xp = xE (XXV )
t → +∞ w
x (t ) x (t )
x0 > xp ; w < 0
x0 x0 > xp ; w > 0
x0
xp xp
x0
x0 < xp ; w > 0
x0
x0 < xp;w < 0
t t
(a ) (b )
Figura 3
Ejemplos:
1.- Beneficio e inversión (Tu, 1994): En una economía donde el beneficio
“ π ” es una función decreciente del stock de capital “K” y la inversión
(I ≡ K ) es una función creciente del beneficio,
'
π = −θK
I = K ' = σπ + I A
247
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Donde θ, σ son coeficientes constantes positivos e “IA” es la inversión
autónoma, asumida constante. Se pide hallar el comportamiento temporal
del stock de capital y su valor de largo plazo, si se sabe que K(0) = K 0 > 0.
σθK = I A ⇒ K E = K p = I A σθ
IA
K (t ) = + Ce − σθt
σθ
IA IA
K(0) = + C = K0 ⇒ C = K0 −
σθ σθ
IA I A − σθt
K (t ) = + K 0 − e
σθ σθ
IA IA −σθt IA
K LP = lím K (t ) = lím + K 0 − e
= = KE
t → +∞ t → +∞ σθ σθ σθ
K (t ) = K E + (K 0 − K E )e − σθt
Q D (t + dt ) = α + βp (t + dt ) (función de demanda )
Q S (t + dt ) = θ + γp ( t ) (función de oferta )
Q D (t + dt ) = Q S (t + dt ) (Ecuación de equilibrio )
248
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Donde α, θ, γ son parámetros positivos, β < 0 , y se asume que los
productores esperan que en el siguiente instante (t + dt) el precio sea igual
al precio actual. Determine el comportamiento temporal del precio de
equilibrio (aquel que vacía el mercado) asumiendo que dt = 1, y que
p(0) = p 0 > 0. Asimismo, determine las condiciones para que el precio de
p (t + dt ) − p (t ) = p ' (t ) ⇒ p (t + dt ) = p (t ) + p ' (t )
[
Q D (t + dt ) = α + β p (t ) + p ' (t )
]
Q S (t + dt ) = θ + γp (t ) (10 )
Q D (t + dt ) = Q S (t + dt )
[ ]
α + β p (t ) + p ' (t ) = θ + γp (t ) ⇒ p ' (t ) +
β−γ
β
p (t ) =
θ−α
β
(11)
(β − γ )
θ−α − t
p (t ) = + Ce β
(12 )
β−γ
θ−α θ−α
p(0) = + C = p0 ⇒ C = p0 −
β−γ β−γ
249
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando “C” en (12), se tiene:
(β − γ )
θ−α θ−α − t
p (t ) = + p0 − e β
(13)
β−γ β−γ
(β − γ)
θ − α θ − α − β t
p LP = lím p(t) = lím + p0 − e
t → +∞ t → +∞ β − γ β − γ
β−γ
>0
β
(β − γ)
θ − α θ − α − β t θ−α
p LP = lím p(t) = lím + p0 − e
=
t → +∞ t → +∞ β − γ β−γ β−γ
θ−α
>0
β−γ
(β − γ ) (θ − α ) θ−α
p= ⇒ pE = = p LP
β β β−γ
250
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Por lo que (13) podría escribirse como:
(β − γ )
− t
p (t ) = p E + (p 0 − p E )e β
(14 )
Derivando (14) respecto del tiempo obtenemos la tasa de cambio
instantánea del precio:
(β − γ )
β−γ − t
p (t ) = −
'
(p 0 − p E )e β
(15 )
β
θβ − αγ
Q D (t + dt ) = Q S (t + dt ) = + γ [p ( t ) − p E ]
β−γ
θβ − αγ
Q LP = lím [Q D (t + dt )] = lím [Q S (t + dt )] =
t → +∞ t → +∞ β−γ
a 0 (t ) f (t )
a 1 (t )x ' + a 0 ( t )x = f (t ) ⇒ Si (a 1 (t ) ≠ 0 ∀t ) ⇒ x ' + x=
a 1 (t ) a 1 (t )
x ' + w (t )x = q (t ) (XXVI )
Su solución general la obtendremos multiplicando (XXVI) por el factor
w ( t )dt
e∫ , denominado factor de integración, de donde tenemos que:
251
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Resultando que:
d xe ∫
w ( t )dt
= x ' e ∫ w ( t )dt + w (t )xe ∫ w (t )dt = q (t )e ∫ w ( t )dt
dt
w ( t )dt
d xe ∫
w ( t )dt w ( t )dt w ( t )dt
∫ q (t )e ∫ dt ⇒ xe ∫ ∫ q (t )e ∫
∫
dt
dt =
= dt + C
− w ( t )dt
∫ ∫ w (t )dt dt + C
x (t ) = e q (t )e
∫ (XXVII )
x ' + w (t )x = 0 (XXVIII)
Cuya solución la obtenemos reemplazando q(t) = 0 en (XXVII):
− w ( t )dt
∫
x (t ) = Ce (XXIX )
Ejemplos:
1.- Crecimiento económico (Sydsæter, 2005): considere el siguiente modelo
de crecimiento económico en un país en vías de desarrollado:
252
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Determine el comportamiento temporal de “ K(t) ” si
M(t) = M 0e , K(0) = K 0 y βα ≠ θ. Encuentre cómo evoluciona en el tiempo
θt
K(t) = e ∫
βα dt
∫
θt
M 0e e
∫
− βαdt
[∫
dt + C ⇒ K(t) = eβαt M 0eθt e −βα t dt + C ]
[ ] M0
K(t) = eβα t M 0 e(θ −βα)t dt + C ⇒ K(t) = eβα t
∫
e(θ−βα)t + C
θ − βα
M0
K(t) = eθt + Ceβα t
θ − βα
M0 M0
K(0) = + C = K0 ⇒ C = K0 −
θ − βα θ − βα
M0 M 0 βα t
K(t) = eθt + K 0 − e
θ − βα θ − βα
M 0α M0 βα t
Y (t ) = e θt + α K 0 − e
θ − βα θ − βα
Y (t ) M 0α α M0 (βα − γ )t
y (t ) = = e ( θ − γ )t + K0 − e
N (t ) N 0 (θ − βα )
N0 θ − βα
253
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
2.- Modelo de bonos y tasas de interés: (Lomelí, 2003): Sea B(t ) el valor
(precio) de un bono al instante “t”. Se supone que la tasa de cambio relativo
instantánea del valor de este bono, la tasa de interés instantánea del bono al
tiempo “t”, viene dada por:
B ' (t )
r (t ) = ⇒ B ' (t ) = r (t )B(t ) ⇒ B ' (t ) − r (t )B(t ) = 0 (16 )
B (t )
Supóngase que a lo largo del periodo de vida del bono se realizan varios
depósitos (o retiros), y que el valor de la inversión en bonos al instante “t”
viene dada por:
Y (t ) = Z (t )B(t ) (17 )
Donde Z(t ) representa la cantidad de bonos que se tienen en la inversión
al instante “t”. Estando el cambio marginal de Z(t ) dado por:
Z ' (t ) = δ (t ) (18 )
Bajo todos los supuestos antes hechos, ahora vamos a calcular la variación
marginal del valor de la inversión en bonos en función de la variación
marginal de la cantidad de bonos que se tienen en la inversión y de la
variación marginal del valor de un bono, al instante “t”. Para ello, vamos a
derivar (8) respecto de “t”, de donde obtenemos:
254
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Para ello, primero resolveremos la ecuación (16) asumiendo que
B(T ) = B T , y tomando como instante inicial “T” y como instante final “t”.
t
− r ( τ )dτ
∫
Multiplicando (16) por e T se tiene:
t
− ∫ r(τ)dτ
d B(t)e
T
t t
− r(τ)dτ
∫ − r(τ)dτ
∫
B'(t)e T − r(t)B(t)e T =0⇒ =0
dt
t t
− r ( τ )dτ
∫ ∫ r (τ )dτ
B(t )e T = C ⇒ B(t ) = Ce T (22 )
∫ r (τ )dτ
B(T ) = Ce T = BT ⇒ C = BT (23)
t
− r ( τ )dτ
∫
Ahora, para resolver (21) multiplicaremos dicha ecuación por e T :
t t t
− r ( τ )dτ
∫ − r ( τ )dτ
∫ − r ( τ )dτ
∫
Y ' (t )e T − r (t )Y (t )e T = δ (t )B(t )e T
t
− ∫ r ( τ )dτ
d Y (t )e T
t
− ∫ r ( τ )dτ
= δ (t )B(t )e T (25 )
dt
t
− ∫ r ( τ )dτ
d Y (t )e T
= δ (t )B T (26 )
dt
255
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Integrando (26) obtenemos:
t
− ∫ r(τ)dτ
dY(t)e T
t t
∫ dτ
dτ = ∫ BTδ(τ)dτ
T T
t t
− r ( τ )dτ
∫ t − r ( τ )dτ
∫ t
Y (t )e T − Y (T ) = B T δ (τ )dτ ⇒ Y (t )e
∫ T = Y (T ) + B T δ (τ )dτ
∫
T T
t T
− r(τ)dτ
∫ − r(τ)dτ
∫ t
Y(t)e T − Y(T)e T = BT δ(τ)dτ ∫
T
t ∫ r ( τ )dτ
Y (t ) = Y (T ) + B T δ (τ )dτ e T
∫ (27 )
T
∫ r (τ )dτ
Reemplazando Y (T ) = YT y el valor de e T que aparece en (24)
obtenemos que:
Y t
Y (t ) = + δ (τ )dτ B(t )
∫ (28 )
T
BT
T
1
r ( t ) = r0 −
t +1
T +1
δ (t )B(t ) = − (29 )
t +1
B(T ) = 1
Y (T ) = 1
256
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Dado que:
t ≥ 0 ⇒ t + 1 > 0 ⇒ t + 1 = t + 1
T > t ⇒ T > 0 ⇒ T + 1 > 0 ⇒ T + 1 = T + 1
Por tanto:
t +1
r0 ( t − T )− ln
(T + 1)e r0 (t −T )
B(t ) = e [r0 t − ln ( t +1)− r0T + ln (T +1)] = e T +1
= (30 )
t +1
(T + 1)e r0 ( t −T ) T +1
δ (t ) =− ⇒ δ (t ) = − e r0 (T − t ) (31)
t +1 t + 1
e r0 (T − t ) − 1 (T + 1)e r0 ( t − T )
Y (t ) = 1 + ⋅
r0 t +1
k ' (t ) = f (k (t )) − c (t ) − (n + g )k (t )
' f ' (k (t )) − ρ − θg
c (t )
(33)
c ( t ) =
θ
Con las condiciones iniciales: k(0) = k 0, c(0) = c0. Donde el stock de capital
por mano de obra efectiva, k(t), y el consumo por mano de obra efectiva,
c(t), son considerados como variables endógenas. Siendo “n” la tasa de
crecimiento instantáneo de la población, “g” la tasa instantánea de
progreso tecnológico, “ρ” es la tasa de descuento intertemporal, y “θ” es el
coeficiente de aversión relativa al riesgo.
f (k (t )) = ak (t ); a >0 (34 )
257
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Se pide determinar la evolución temporal del stock de capital k(t).
f ' (k (t )) = a ; (35 )
k ' (t ) = (a − n − g )k (t ) − c (t )
' a − ρ − θg (36 )
c ( t ) =
c (t )
θ
a − ρ − θg
c ' (t ) − c (t ) = 0
(37 )
θ
a − ρ − θg
t
c (t ) = Ce θ (38 )
a −ρ − θg
t
c (t ) = C = c 0 ⇒ c (t ) = c 0 e θ (39 )
Reemplazando (39) en la primera ecuación de (36) resulta:
a −ρ − θg
t
k (t ) = (a − n − g )k (t ) − c 0
'
e θ
a −ρ − θg
t
k (t ) − (a − n − g )k (t ) = − c 0
'
e θ (40 )
La solución de (40) la podemos obtener a partir de (XXVII):
a − ρ − θg
t
(a − n −g )dt − (a − n − g )dt
k (t ) = e ∫ ∫
∫ − c0 e θ e dt + C
a − ρ − θg
− (a − n − g ) t
k ( t ) = e ( a − n − g )t − c 0 e
θ
∫ dt + C
258
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
a − ρ − θg
− ( a − n − g ) t
θ
( a − n − g )t − c 0 e
k (t ) = e + C
a − ρ − θg
− (a − n − g )
θ
a −ρ − θg
t
− c0 e θ
( a − n − g )t
k (t ) = + Ce (41)
a − ρ − θg
− (a − n − g )
θ
−c 0 c0
k (0 ) = + C = k0 ⇒ C = k0 + (42 )
a −ρ − θg a −ρ − θg
− (a − n − g ) − (a − n − g )
θ θ
a − ρ − θg
t
− c0 e θ c0 e (a − n − g )t
k (t ) = + k 0 +
a − ρ − θg a −ρ − θg
− (a − n − g ) − (a − n − g )
θ θ
259
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
La solución complementaria puede encontrarse realizando la combinación
lineal de “n” soluciones linealmente independientes. Es decir:
n
x c (t ) = ∑ kixi =k 1 x 1 + k 2 x 2 + K + k n x n (XXXII )
i =1
x1 x2 K xi K xn
x1' x '2 K xi' K x 'n
M M M M M M
W(t) = ≠0 con (0 ≤ i ≤ n)
x1i) xi2) K xii) K xin)
M M M M M M
x1n −1) x n2 −1) L xin −1) K x nn −1)
P (r ) = r n + a n −1 r n −1 + K + a 2 r 2 + a 1 r + a 0 = 0 (XXXIII)
n1 n2 n3
x c (t ) = ∑ A i e ri t + ∑ B i t (i −1) e rt + ∑ e p t [ C i cos (q i t ) + G i sen (q i t ) ] (XXXIV )
i
i =1 i =1 i =1
260
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Dónde:
n W (t )
x p (t ) = ∑ x i ∫ Wi(t ) dt (XXXVI )
i =1
x1 x2 K 0 K xn
x1' x'2 K 0 K x 'n
M M M M M M
Wi(t) =
x1k) x k2) K 0 K x kn)
M M M M M M
x1n −1) x n2 −1) n −1)
L f (t) K x n
c (t4
644444444444x7 ) 44444444448
n1 n2 n3
x(t) = ∑ Ai e ri t + ∑ Bi t (i −1)e rt + ∑ ep t[ Ci cos (qi t) + G isen(qi t) ] +
i
i =1 i =1 i =1
(XXXVII )
n Wi (t)
+ ∑ x i ∫ W(t) dt
i =1442443
1
x p (t)
q
x ' = q − wx E = g (x E ) = 0 ⇒ x E = (XXXVIII)
w
261
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Un punto de equilibrio x E es estable si para cada ε > 0 hay un “δ” tal que cada
trayectoria x (t ) con x (0 ) − x E < δ satisface x (t ) − x E < ε ∀t ≥ 0. En otras
palabras, si una trayectoria que empieza “cerca” a un punto de equilibrio permanece
cerca a este punto durante todo el tiempo futuro, entonces el punto de equilibrio se
dice que es estable.
Una condición necesaria y suficiente para que la ecuación (XXX) sea estable
(globalmente asintóticamente estable) es que todas las raíces características
(autovalores) de (XXXIII), tengan partes reales negativas.
262
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Es claro que la solución general de (XXX) poseerá un equilibrio fijo si tanto
su solución complementaria como su solución particular convergen hacia un
valor real. La convergencia de la solución complementaria se garantizará si la
parte real de las raíces características es negativa, con lo cual se cumple que
el lím x c(t) = 0, mientras que la convergencia de la solución particular se
t → +∞
garantizará si se cumple que el lím x p (t) = x E ∈ ℜ .
t → +∞
Ejemplos:
Resuelva las siguientes ecuaciones diferenciales:
a b−p
1. x '' − x = . Si se sabe que a, A, b y p son constantes positivas y que
A 2A
las condiciones iniciales son: x(0) = x 0 y x(T) = x T.
a
r1 = >0
a A
( ) 2
Pr = r − =0⇒
A a
r2 = − <0
A
x c (t ) = A 1 e (43)
a At − a At
+ A 2e
x 2(t) = e ⇒ x '2(t) = − a A e
− a At − a At
El Wronsquiano será:
a At − a At
W(t) = e e = −2 a A ≠ 0
a At − a At
a Ae − a Ae
− a At
0 e p−b
W1(t) = b − p
− a At
− a At = e
− a Ae 2A
2A
263
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
a At
e 0 b−p
W2(t) =
a At
a At b−p = e
a Ae 2A
2A
Por tanto:
p−b
x (0 ) = A 1 + A 2 + = x0
2a
p−b
A1 + A 2 = x 0 − (46 )
2a
p−b
x (T ) = A 1 e
a AT − a AT
+ A 2e + = xT
2a
p−b
A1e
a AT
+ A2e
− a AT
= xT − (47 )
2a
p−b p−b a AT
x0 − − xT − e
2a
A1 = 2a
2 a A T
1 − e
264
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
a A T p−b p−b a A T
e x T − − x0 −
e
2a 2a
A2 =
2 a AT
1 − e
p−b p−b a AT
x 0 − − xT −
e
2a 2a
x (t ) =
a At
e +
2 a AT
1 − e
a AT p−b p−b a AT
e x T − − x0 −
e
2a 2a − a At p−b
+ e +
2 a AT 2a
1 − e
5 1
2.- x '' + x ' − 2x = 5t . Si se sabe que: x (0 ) = − y x ' (0 ) = .
4 2
r1 = 1
P(r) = r 2 + r − 2 = (r − 1)(r + 2) = 0 ⇒ .
r2 = −2
x c ( t ) = A 1 e t + A 2 e −2 t
Dónde:
x 1 (t ) = e t ⇒ x 1' (t ) = e t
x 2 (t ) = e −2 t ⇒ x '2 (t ) = −2 e −2 t
265
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Por tanto, el Wronsquiano será:
et e −2 t
W (t ) = −2t
= −3e − t < 0 ⇒ W (t ) ≠ 0
et − 2e
Mientras que:
0 e −2 t
W1 (t ) = = −5 te − 2 t
5 t − 2e −2 t
et 0
W2 (t ) = = 5 te t
et 5t
Por tanto:
W1 (t ) − 5 te −2 t 5 5
∫ dt = ∫ dt = ∫ te
−t
dt = − (t + 1)e − t
W (t ) − 3e −t 3 3
W1 (t ) 5 5
x 1 (t ) ∫ dt = e t − (t + 1)e − t = − (t + 1)
W (t ) 3 3
W 2 (t ) 5 te t 5 5
∫ dt = ∫ − 3e − t dt = − 3 ∫ te
2t
dt = − e 2 t (2 t − 1)
W (t ) 12
W 2 (t ) 5 2t 5
x 2 (t ) ∫ dt = e − 2 t − e (2 t − 1) = − (2 t − 1)
W (t ) 12 12
Dónde:
5 5
x p (t ) = − (t + 1) − (2 t − 1)
3 12
5 5
x (t ) = A 1 e t + A 2 e − 2 t − (t + 1) − (2 t − 1)
3 12
266
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Simplificando:
5
x (t ) = A 1 e t + A 2 e − 2 t − (2 t + 1)
4
5
x ' (t ) = A 1 e t − 2 A 2 e − 2 t −
2
5 5
x (0 ) = A 1 + A 2 − =− (48)
4 4
5 1
x ' (0 ) = A 1 − 2 A 2 − = (49 )
2 2
A 1 = 1 y A 2 = −1 .
Figura 4
267
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
135
3. x'' ' + 13x' − 34x = e−2t . Si las condiciones iniciales son: x(0) = ,
68
103 1
x '(0) = y x' '(0) = − .
34 17
r1 = 2
(
P (r ) = r 3 + 13r − 34 = (r − 2 ) r 2 + 2 r + 17 = 0 ⇒ ) p1 = −1.
r = −1 ± 4 j ⇒
q1 = 4
Donde:
x1' (t) = 2e 2t
x1(t) = e 2t ⇒
x1''(t) = 4e 2t
e 2t e −t cos(4t) e −t sen(4t)
W(t) = 2e 2t
− e [cos(4t) + 4sen(4t)]
−t
e [4 cos(4t) − sen(4t)] = 100 ≠ 0
−t
Es decir, x1(t) = e2t, x 2(t) = e−t cos(4t) y x3(t) = e−tsen(4t) son soluciones de
xc(t) que son linealmente independientes.
Mientras que:
268
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
e2t 0 e−tsen(4t)
W2(t) = 2e2t 0 e− t[4 cos(4t) − sen(4t)] = e− t[3sen(4t) − 4 cos(4t) ]
4e2t e − 2t
− e− t[8 cos(4t) + 15sen(4t)]
Por tanto:
W1(t) 4e−4t 1 1
∫ W(t) ∫ 100 dt = 25 ∫ e
− 4t
dt = dt = − e− 4t
100
W1(t) 2t 1 − 4t 1 − 2t
x1(t)
∫ W(t) dt = e
− 100 e = − 100 e
W2(t) 1
x 2(t) cos(4t) [19sen(4t) + 8 cos(4t) ]
∫ W(t) dt = − 1700 e
− 2t
W3(t) 1
x3(t) sen(4t) [19 cos(4t) − 8sen(4t) ]
∫ W(t) dt = 1700 e
− 2t
Donde:
Por tanto, la solución general de x'' ' + 13x ' − 34x = e−2t es:
269
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
e−2t e−2t cos(4t) [19sen(4t) + 8 cos(4t) ]
x(t) = A1e2t + e− t [C1 cos(4t) + G1sen(4t)] − − +
100 1700
1 8
x(t) = A1e2t + e− t[C1 cos(4t) + G1sen(4t)] − e− 2t − e− 2t
100 1700
1
x(t) = A1e2t + e− t[C1 cos(4t) + G1sen(4t)] − e− 2t
68
1
x '(t) = 2A1e2t + e− t[(4G1 − C1) cos(4t) − (G1 + 4C1) sen(4t)] + e− 2t
34
1
x ' '(t) = 4A1e2t + e− t[(8C1 − 15G1) sen(4t) − (8G1 + 15C1) cos(4t)] − e− 2t
17
1 135
x (0 ) = A 1 + C1 − = ⇒ A 1 + C1 = 2 (50 )
68 68
1 103
x ' (0 ) = 2 A 1 + 4G 1 − C1 + = ⇒ 2 A 1 + 4G 1 − C1 = 3 (51)
34 34
1 1
x '' (0 ) = 4 A 1 − 8G 1 + 15C1 − =− ⇒ 4 A 1 − 8G 1 + 15C1 = 0 (52 )
17 17
A1 = 4, C1 = −2 y G1 = −7 / 4
7 1 − 2t
x(t) = 4e2t + e− t − 2 cos(4t) − sen(4t) − e
4 68
270
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Figura 5
Figura 6
271
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
4.- Modelo IS-LM de la Economía (Tu, 1994): Considere una economía
cerrada descrita por las siguientes ecuaciones:
Y' = h(D − S)
'
[
r = m L(Y) − M ]
D = C+ I
C = cY; (0 < c < 1)
(53)
I = − ar; (a > 0)
S = Y
L(Y) = kY; (k > 0)
s ≡ 1 − c
La primera ecuación en (53) nos indica que la renta nacional “Y” aumenta
en respuesta al exceso de demanda agregada “D”, la segunda ecuación nos
indica que la tasa de interés “r” aumenta en respuesta al exceso de demanda
de dinero (diferencia entre la demanda y oferta de dinero: “ L (Y ) − M ”), la
tercera ecuación nos indica que la demanda agregada es igual a la suma del
consumo y de la inversión, la cuarta y la quinta ecuación nos indican que es
el consumo y la inversión son funciones lineales de la renta y de la tasa de
interés respectivamente, la sexta ecuación nos señala que la oferta agregada
no es más que la producción nacional, la sexta ecuación nos indica que la
demanda de dinero es una función lineal creciente de la renta (es decir, el
dinero demandado únicamente con el propósito de realizar transacciones), la
séptima ecuación no es más que la definición de la propensión marginal a
ahorrar “s” en términos de la propensión marginal al consumo “c”. Se asume
que la oferta de dinero “ M ” es realizada por el Banco Central, y que las
velocidades de ajuste “h” y “m” son constantes e iguales a 1.
Y' = C + I − Y
(54 )
r' = kY − M
272
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Y' = cY − ar − Y = −sY − ar
' (56 )
r = kY − M
( )
Y '' = − sY ' − ar ' = −sY ' − a kY − M ⇒ Y '' + sY ' + akY = aM (57 )
−s+ s2 − 4ak −s+ ∆
r1 = =
P(r) = r 2 + sr + ak = (r − r1)(r − r2) = 0 ⇒ 2 2
−s− s2 − 4ak −s− ∆
r2 = =
2 2
r1 ∈ ℜ ∧ r1 < 0
Caso I: ∆ = s 2 − 4ak > 0 ⇒ r2 ∈ ℜ ∧ r2 < 0
r ≠ r ∧ r < r
1 2 2 1
Yc (t ) = A 1 e r1t + A 2 e r2 t (58)
Dónde:
er1t er2t
W(t) = = (r2 − r1)e(r1 + r2 )t = − ∆ e(r1 + r2 )t ≠ 0
r1er1t r2er2 t
Mientras que:
er2t
W1(t) =
0
= −aMer2t
aM r2er2 t
273
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
er1t
W2(t) =
0
= aMer1t
r1er1t aM
Por tanto:
W1(t) aM aM
Y1(t) ∫ W(t) dt = −e
r1t
⋅ e− r1t = −
r1 ∆ r1 ∆
W2(t) aMer1t aM − r2 t aM
∫ W(t) dt = ∫− dt = − ∫e dt = e− r2t
∆ e(r1 + r2 )t ∆ r2 ∆
W2(t) aM aM
Y2(t) ∫ W(t) dt = e
r2 t
⋅ e− r2 t =
r2 ∆ r2 ∆
Dónde:
aM aM aM 1 1 aM r1 − r2 aM r1 − r2 M
Yp(t) = − + = − = = =
r1 ∆ r2 ∆ ∆ r2 r1 ∆ r ⋅r
1 2 ∆ r ⋅r
1 2 k
M
Y (t ) = A 1 e r1t + A 2 e r2 t + (59 )
k
M
YLP = lím Y (t ) =
t → +∞ k
s
Caso II: ∆ = s 2 − 4 ak = 0 ⇒ r1 = r2 = r = − ∈ℜ ∧ r < 0
2
Yc (t ) = B1 e rt + B 2 te rt (60 )
Dónde:
Y1 (t ) = e rt ⇒ Y1' (t ) = re rt
Y2 (t ) = te rt ⇒ Y2' (t ) = e rt + tre rt = e rt (1 + tr )
274
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Por tanto el Wronsquiano será:
e rt te rt
W (t ) = = e 2 rt ≠ 0
re rt e (1 + tr )
rt
Mientras que:
0 te rt
W1 (t ) = = − aMte rt
aM e rt (1 + tr )
e rt 0
W2(t) = = aMe rt
rert aM
Por tanto:
W1 (t ) − aMte rt aM 1
∫ W (t ) dt = ∫ dt = − aM te − rt dt =
∫ e − rt t +
2 rt
e r r
W1 (t ) aM 1 aM 1
Y1 (t ) ∫ W (t ) dt = e
rt
⋅ e − rt t + = t +
r r r r
W2 ( t ) aMe rt aM
∫ dt = ∫ dt = aM e − rt dt = −
∫ e − rt
W (t ) e 2 rt
r
W2 ( t ) aM aM
Y2 (t ) ∫ dt = − te rt ⋅ e − rt = − t
W (t ) r r
Dónde:
aM 1 aM aM
Yp (t ) = t + −
t=
r r r r2
aM
Y (t ) = B1 e rt + B 2 te rt + (61)
r2
s
Reemplazando r = − < 0 en (61) tenemos:
2
aM
Y (t ) = B1 e −(s 2 )t + B 2 te −(s 2 )t + 4 (62 )
s2
275
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
s
Ya que r = − < 0, entonces (62) es estable. En este caso, la renta
2
nacional convergerá en el largo plazo a:
aM
YLP = lím Y (t ) = 4
t → +∞ s2
r = p1 ± q1 j
Caso III: ∆ = s 2 − 4ak < 0 ⇒
p1 ∈ ℜ; q1 ∈ ℜ − {0}
Dónde:
Es decir, Y1(t) = e p1t cos(q1t) y Y2(t) = e p1t sen(q1t) son soluciones de Yc (t)
que son linealmente independientes.
Mientras que:
sen(q1t)e p1t
W1(t) =
0
= −aMsen(q1t)e p1t
aM [p1sen(q1t) + q1 cos(q1t)]ep t 1
cos(q1t)e p1t
W2(t) =
0
= aM cos(q1t)ep1t
[p1 cos(q1t) − q1sen(q1t)]e p1t
aM
Por tanto:
276
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
W1(t) aM
∫ W(t) dt = q (p 2 + q 2 ) e − p1t [p1sen(q1t) + q1 cos(q1t)]
1 1 1
W1(t) aM
Y1(t)
∫ W(t) dt = e
p1t
cos(q1t) ⋅ e − p1t [p1sen(q1t) + q1 cos(q1t)]
(
q1 p12 + q12 )
W1(t) aM
Y1(t)
∫ W(t) dt = q (p 2 + q 2 ) cos(q1t) ⋅ [p1sen(q1t) + q1 cos(q1t)]
1 1 1
W2(t) aM
∫ W(t) dt = q (p 2 + q 2 ) e − p1t [q1sen(q1t) − p1 cos(q1t)]
1 1 1
W2(t) aM
Y2(t) ∫ W(t) dt = e
p1t
sen(q1t) ⋅ e − p1t [q1sen(q1t) − p1 cos(q1t)]
(
q1 p12 + q12 )
W2 (t) aM
Y2 (t) ∫ W(t) dt = sen(q1t) ⋅ [q1sen(q1t) − p1 cos(q1t)]
(
q1 p12 + q12 )
Dónde:
aM
Yp (t) = cos(q1t) ⋅ [p1sen(q1t) + q1 cos(q1t)] +
(
q1 p12 + q12 )
aM
+ sen(q1t) ⋅ [q1sen(q1t) − p1 cos(q1t)]
(
q1 p12 + q12 )
aM
Yp (t) =
(
p12 + q12 )
Por tanto, la solución general de (57) en este caso será:
aM
Y(t) = e p1t [C1 cos(q1t) + G1sen(q1t)] + (63)
(p 2
1 + q12 )
Pero en este caso, dado que ∆ = s 2 − 4ak < 0, entonces tenemos que:
277
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
r = − s ± ∆
=−
s
±
(
− 4ak − s 2 ) =−
s
±
(4ak − s ) ⋅2
−1
2 2 2 2 2 2
s
p1 = − < 0 (64)
s (4ak − s ) 2
2
r = − ± j = jp1 + q1 j ⇒
2 2
q1 =
4ak − s 2
>0
( )
2
s2
p12 =
4
2
⇒ p12 + q12 = ak (65)
2 4ak − s
q1 =
4
Y(t) = e −(s 2)t C1 cos
(4ak − s ) t + G sen (4ak − s ) t + M
2 2
(66)
1
2 2 k
2
El valor absoluto del coeficiente de aversión absoluta al riesgo aproximadamente mide la variación
relativa de la utilidad marginal del agente ante un incremento unitario del dinero, para una cantidad “x” de
dinero. Aquellas funciones de utilidad que poseen un coeficiente de aversión absoluta al riesgo constante
suelen ser conocidas como funciones de utilidad CARA (Constant Absolute Risk Averse) o como funciones
de utilidad exponenciales.
278
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Dado que en este caso se requiere que el coeficiente de aversión absoluta
al riesgo sea constante, entonces:
r1 = 0
p(r) = r 2 + λ a r = 0 ⇒ p(r) = r(r + λ a ) = 0 ⇒ (70)
r2 = −λ a < 0
lím u(x) = A1
t → +∞
dp (t ) t
dt
= k⋅ ∫ [D(p (υ )) − S(p (υ ))]dυ (k > 0 ) (72 )
−∞
dp(t) e t
= k ⋅ [D(p(υ)) − S(p(υ))]dυ + [D(p(υ)) − S(p(υ))]dυ
∫ ∫ (73)
dt
− ∞ e
279
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Derivando (73) respecto de “t”:
d e t
p' '(t) = k ⋅ [D(
∫p(υ)) − S(p(υ))]dυ + [D(p(υ)) − S(p(υ))]dυ
∫
dt − ∞
e
d e t
p' '(t) = k ⋅ [D(p(υ)) − S(p(υ))]dυ + [D(p(υ)) − S(p(υ))]dυ
∫ ∫
dt − ∞
e
d e d t
[D(p(υ)) − S(p(υ))]dυ + [D(p(υ)) − S(p(υ))]dυ
p' '(t) = k ⋅ ∫ ∫
dt
dt − ∞ e
644 47 0 444 8 67 08
e ∂[D(p(υ)) − S(p(υ))] dt de
p' '(t) = k ⋅ ∫ dυ + [D(p(t)) − S(p(t))] − [D(p(e)) − S(p(e))] +
− ∞ ∂t dt dt
644 47 0 444 8
t
∂[D(p(υ)) − S(p(υ))]
+ ∫ ∂t
dυ
e
p1 = 0
P (r ) = r 2 + (h − b )k = 0 ⇒ r = ± (h − b )k j = p1 ± q1 j ⇒
q1 = (h − b )k >0
3
Ver apéndice.
280
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Dónde:
Mientras que:
0 sen (q1t )
W1 (t ) = = k (c − a )sen (q1t )
k (a − c ) q1 cos( q1t )
cos (q1t ) 0
W2 (t ) = = k (a − c ) cos (q1t )
− q1sen ( q1t ) k (a − c )
Por tanto:
W1 (t ) k (c − a )sen (q1t ) k (c − a ) k (a − c )
∫ W (t ) dt = ∫ q1
dt =
q1
∫ sen (q1t )dt = q12
cos (q1t )
W1 (t ) k (a − c ) k (a − c )
p1 (t ) ∫ W (t ) dt = cos(q1t ) ⋅ cos (q1t ) = cos 2 (q1t )
q12 q12
W2 (t ) k (a − c ) cos (q1t ) k (a − c ) k (a − c )
∫ W (t ) dt = ∫ q1
dt =
q1
∫ cos(q1t )dt = q12
sen (q1t )
W2 (t ) k (a − c ) k (a − c )
p 2 (t )∫ W (t ) dt = sen (q1t ) ⋅ sen (q1t ) = sen 2 (q1t )
q12 q12
Dónde:
k (a − c ) k (a − c ) k (a − c )
p p (t ) = cos 2 (q1t ) + sen 2 (q1t ) =
q12 q12 q12
k (a − c )
p (t ) = C1 cos (q1t ) + G1sen ( q1t ) + (76 )
q12
281
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando q1 = (h − b )k en (76) resulta:
a−c
p (t ) = C1 cos (h − b )k ⋅ t + G1sen ( (h − b )k ⋅ t) + (77 )
h−b
Figura 7
dy F(t )
y ' (t ) = = ; G (y ) ≠ 0 (XXXIX )
dt G (y )
du ' du ' dx
xu ''(x) = −λ r u '(x) ⇒ x = −λ r u ' ⇒ '
= −λ r (80)
dx u x
du ' dx
∫ '
= −λ r ∫ ⇒ ln u ' = −λ r ln x + C1 (81)
u x
Dado que:
4
El valor absoluto del coeficiente de aversión relativa al riesgo aproximadamente mide la variación
porcentual de la utilidad marginal del agente ante un incremento respecto al nivel de dinero “x” del 1%.
Aquellas funciones de utilidad que poseen un coeficiente de aversión relativa al riesgo constante suelen ser
conocidas como funciones de utilidad CRRA (Constant Relative Risk Averse) o como funciones de utilidad
de elasticidad de sustitución constante (funciones de utilidad isoelásticas).
283
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando (82) en (81) resulta:
u'
ln u ' = −λ r ln x + C1 ⇒ ln u ' = ln+ C1 ⇒ ln −λ =C (83)
x r 1
u'
ln − λ
x r u'
e = eC1 ⇒ −λ r
= eC1 ⇒ u ' = eC1 x −λ r ⇒ u ' = A1x −λ r (84)
x
du
u' = = A1x − λ r ⇒ du = A1x − λr dx (85)
dx
A1
x 1−λ r + A 2 ; λ r ≠ 1
∫ du = ∫ A1 x
−λ r
dx ⇒ u (x ) = 1 − λ r (86 )
A 1 ln x + A 2 ; λr = 1
dy
=
(
y1 − α y ρ ) (87)
dx x
1 αyρ −1 1
+ = (88)
( )
ρ
y 1 − αy
y1 − α y ρ
(
y = β x −ρ + α ) −1 ρ
(89)
284
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Luego, haciendo x = K L , y = Y L , A = (α + β )−1 ρ , y δ = β (α + β ) ,
demuestre que a partir de (89) se llega a la siguiente forma especial de la
función de producción CES:
[
Y = A δK ρ + (1 − δ )L− ρ ]
−1 ρ
(90 )
dy dx
= (91)
(
y 1 − αy ρ ) x
dy dx
∫ = ∫ (92)
(
y1 − α y ρ ) x
1 ρ −1
+ αy dy = dx
(93)
∫ y 1 − αy ρ ∫ x
1 y 1ρ
ln y − ln 1 − αy ρ = ln x + C ⇒ ln = C + ln 1 − αy ρ (94 )
ρ x
y 1 ρ
ρ ρ
C + ln 1− αy
ln y 1ρ y
e x
= e ⇒ = e C 1 − αy ρ ⇒ = e ρ C 1 − αy ρ (95 )
x x
Dado que:
y > 0 ⇒ y = y
x > 0 ⇒ x = x (96 )
ρC
e > 0 ⇒ e ρC = e ρC
ρ ρ ρ ρ
y y y y y y
= >0⇒ = > 0 ⇒ = (97 )
x x x x x x
285
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando (96) y (97) en (95) se obtiene:
ρ ρ
y
x
y
= e ρC 1 − αy ρ = e ρC 1 − αy ρ ⇒ ( )
= e ρC 1 − α y ρ
( ) (98)
x
ρ
y
( )
= θ 1 − αy ρ ⇒ y ρ = θx ρ − αθy ρ x ρ ⇒ y ρ 1 + αθx ρ = θx ρ
x ( )
(
y = θ1 ρ x 1 + αθx ρ ) −1 ρ
(
⇒ y = θ −1 x − ρ + α ) −1 ρ
(99 )
(
y = βx − ρ + α )
−1 ρ
(100 )
x = K L
y = Y L
A = (α + β )−1 ρ (101)
β α
δ = ⇒ 1− δ =
α+β α+β
Se tiene:
−1 ρ
K −ρ
Y
= β + α
(
⇒ Y = βK −ρ + αL−ρ )
−1 ρ
L L
[
Y = δ (α + β )K −ρ + (1 − δ )(α + β )L−ρ ]
−1 ρ
[
Y = (α + β )−1 ρ δK − ρ + (1 − δ )L− ρ ]−1 ρ
[
Y = A δK ρ + (1 − δ )L− ρ ]
−1 ρ
lqqd
286
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
3.- Supóngase que Y = (a + αK ) t + 1 denota la producción en función del
capital “k”, donde el factor t + 1 se debe al progreso tecnológico.
Supóngase que una fracción constante 0 < s < 1 es ahorrada, y que la
acumulación del capital es igual a los ahorros, de manera que tengamos la
siguiente ecuación diferencial:
dK
= sY = s (a + αK ) t + 1
dt (102)
K (0 ) = K
0
dK dK
a + αK
= s t + 1 dt ⇒ ∫ a + αK = ∫ s t + 1 dt
2
ln a + αK = s (t + 1)3 2
+C (103)
3
a + αK = a + αK (104 )
2
s ( t + 1)
3 2
2 + C
ln (a + αK )
ln (a + αK ) = s (t + 1) 3 2
+C⇒e = e 3
3
2 2
s ( t + 1) 3 2
eC s ( t + 1) 3 2
a
a + αK = e ⋅ e C 3 ⇒ K (t ) = ⋅e 3 − (105 )
α α
2
eC s a eC a
K (0 ) = ⋅e 3 − = K0 ⇒ = K 0 + e − 2 s 3
(106 )
α α α α
a
2s
[(t +1) 3 2
−1 ] a
K (t ) = K 0 + ⋅e
3 −
α α
287
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
2. EDOs con coeficientes homogéneos
Sea una EDO de primer orden en su forma diferencial:
P ( x , t ) dx + Q (x , t ) dt = 0 (XLII)
P (αx , αt ) = α γ P (x , t )
(XLIII)
Q (αx , αt ) = α γ Q (x , t )
dt p(t x) dt
x γ p(t x) dx + x γ q(t x) dt = 0 ⇒ =− ⇒ = f (t x) (XLV)
dx q(t x) dx
t
= θ ⇒ t = θx (XLVI)
x
dt dθ
dt = θdx + xdθ ⇒ =θ+x (XLVII)
dx dx
dθ dθ f (θ) − θ dx dθ
f (θ) = θ + x ⇒ = ⇒ = (XLVIII)
dx dx x x f (θ) − θ
288
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Ejemplos:
1.- (Weber, 1982): Suponga que la tasa de incremento en el costo “x” de
realizar un pedido y supervisarlo, a medida que aumenta el tamaño del
pedido a suministrar, es igual a la suma de los cuadrados del costo y el
tamaño del pedido, dividida entre el doble del producto del costo y el
tamaño del pedido. Determinar la relación entre el costo de realizar y
supervisar, y el tamaño del pedido si x = 6 cuando t = 1.
dx x2 + t2 1 1
= = + (t x) = f (t x) (107)
dt 2tx 2 (t x)
Se puede verificar que las funciones P(x, t) y Q(x, t) de (108) son funciones
homogéneas de grado dos:
[x 2
] ) ( ) (
+ (θx)2 (θdx + xdθ) − 2x(θx)dx = 0 ⇒ θ θ2 − 1 x 2dx + 1 + θ2 x 3dθ = 0
θ (θ − 1)dx + (1 + θ ) xdθ = 0 ⇒
2 2 dx
+
(1 + θ ) dθ = 0 2
θ(θ − 1)
x 2
dx (1 + θ ) dθ = C ⇒ ln x + 2θ − 1 dθ = C
2
∫ + ∫ ∫
θ(θ − 1) ( θ − 1) θ
x 2 2
θ dθ
ln x + 2 ∫ dθ − ∫ =C
(θ − 1)
2 θ
ln x + ln θ2 − 1 − ln θ = C ⇒ ln
(
x θ2 − 1 ) =C
θ
289
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
( )
x θ −1
2
ln
e
θ
= eC ⇒
(
x θ2 − 1 ) =e C
>0⇒
(
x θ2 − 1 )=e C
θ θ
(
x θ2 − 1 )=e C
⇒x=
θ
eC (109)
2
θ θ −1
(t x) C
x=
(t x)2 − 1
e ⇒ x = ± t t − eC ( ) (110)
(
x = t t − eC ) (111)
Reemplazando las condiciones de frontera en (110) tenemos:
6 = 1 − eC ⇒ eC = −35 (112)
Reemplazando (112) en (111) tenemos:
Figura 8
290
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
2.- (Weber, 1982): La relación entre el costo promedio C , y el número de
unidades producidas “q”, es tal que el cambio en el costo promedio a
medida que crece el número de unidades es igual a la diferencia del
número de unidades y el costo promedio, dividida dicha diferencia entre el
número de unidades. Determine la relación entre el costo promedio y el
número de unidades producidas si C = 9 cuando q = 2. Grafique la
relación obtenida.
dC q−C 1
= =1− = f (q C) (114)
dq q (q C)
Expresando (114) en su forma diferencial se tiene:
({q) dC + (1
C2 q)dq = 0
−3 (115)
( )
P C,q Q C,q( )
Se puede verificar que las funciones P C, q y Q C, q de (115) son ( ) ( )
funciones homogéneas de grado uno:
( ( ) )
P α C, αq = αq = αP C, q
C dC dθ
= θ ⇒ C = θq ⇒ d C = θdq + qdθ ⇒ =θ+q (116)
q dq dq
q − θq dθ dq dθ dq dθ
q
=θ+q
dq
⇒
q
+
2θ − 1
=0⇒ ∫ q
+ ∫ 2θ − 1 = K
1 2
ln q + ln 2θ − 1 = K ⇒ 2 ln q + ln 2θ − 1 = 2K ⇒ ln q + ln 2θ − 1 = 2K
2
ln q 2 (2θ −1)
ln q 2 + ln 2θ − 1 = 2K ⇒ ln q 2(2θ − 1) = 2K ⇒ e = e 2K > 0
ln q 2 (2θ −1)
e = e 2K ⇒ q 2(2θ − 1) = e 2K (117)
291
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando (116) en (117) tenemos:
C
q 2 2
q
− 1 = e 2K
(118)
e 2K = 32 (119)
C 1 16
q 2 2 − 1 = 32 ⇒ C(q) = q + (120)
q 2 q
Figura 9
3. EDOs exactas
Supongamos que se tiene una función implícita F(x , t ) = c cuyo diferencial
total es:
∂F(x, t) ∂F(x, t)
dF(x, t) = dx + dt = 0 (XLIX)
∂x ∂t
292
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Ahora, procediendo de manera inversa que en el razonamiento anterior,
supongamos que tenemos la siguiente ecuación diferencial:
∂F(x, t)
= P(x, t)
∂x
(LI)
∂F(x, t)
= Q(x, t)
∂t
∂F(x, t) ∂F(x, t)
P(x, t) dx + Q(x, t) dt = dx + dt = 0 (LII)
∂x ∂t
Una ecuación diferencial ordinaria que tiene la forma general dada por (LII)
será exacta si cumple el siguiente criterio de exactitud:
∂P(x, t) ∂Q(x, t)
= (LIII)
∂t ∂x
∂F(x, t)
∂
∂P(x, t) ∂x ∂ 2 F(x, t)
= = (LIV)
∂t ∂t ∂x∂t
∂F(x, t)
∂
∂Q(x, t) ∂t ∂ 2 F(x, t)
= = (LV)
∂x ∂x ∂t∂x
Técnica de solución
Buscamos una función F(x, t) que satisfaga (LI). Para determinar dicha
función, deberemos proceder de la siguiente forma:
2.- Integrar la primera ecuación que aparece en (LI) respecto a “x”, colocando
en lugar de la usual constante de integración una función genérica que
dependa de “t”: f (t) . En el proceso de integración se deberá considerar a
“t” constante.
∂f (t)
4. Integrar con respecto a “t” para obtener f (t ) .
∂t
∂f (t)
∫ dt = f (t) (LVIII)
∂t
∂Q(x, t) ∂P(x, t)
−
f (t) dt
a) Si ∂x ∂t = f (t) el factor de integración será: µ(t) = e ∫ .
P(x, t)
∂P ( x , t ) ∂Q ( x , t )
−
∫ g(x) dx
b) Si ∂t ∂x = h (x ) el factor integración será: µ(x) = e .
Q (x , t )
Ejemplos:
1.- Elasticidad de la función de costos respecto a la cantidad producida
(Larson, 1995): Si C = C(q) es el costo de producir “q” unidades de
determinado bien, la elasticidad del costo se define como:
20q − C
εq = (122)
2C − 10q
q
C
⋅
dC
dq
=
20q − C
2C − 10q
⇒1 (
4− 10
2Cq
42 43q 2 dC + 1
4 )
C 242
− 20 (
Cq dq = 0
43 ) (123)
P(C,q) Q(C,q)
295
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Siguiendo los pasos de la técnica de resolución de una EDO exacta,
primero derivaremos parcialmente P(C, q) y Q(C, q) respecto a “q” y
a “C” respectivamente para verificar si se cumple o no del criterio
de exactitud dado por (LIII):
∂P(C, q) ∂Q((C, q)
= 2C − 20q = (124)
∂q ∂C
∂F ( C, q )
= P ( C, q )
∂C
(125)
∂F ( C, q )
= Q ( C, q )
∂q
F( C, q ) = ∂F( C, q ) = P ( C, q ) ∂C + f (q )
∫ ∫ (126 )
∂F( C, q ) ∂f (q ) ∂f (q )
= C 2 − 20Cq + = Q (C, q ) = C 2 − 20Cq ⇒ = 0 (128 )
∂q ∂q ∂q
f (q ) = k (129 )
Finalmente, reemplazando (129) en (127) e igualando F(C, q) a una
constante C1 , la solución de (123) resulta:
10q 2 ± 100 q 4 − 4C 2q
qC 2 − 10q 2C + C 2 = 0 ⇒ C(q) = (130)
2q
5 q 2 ± q 4 − 1000000 q
C(q) = (132)
q
5 q 2 + q 4 − 1000000 q
C(q) = ; q ≥ 100
q
dx xq + 12q
=− (133)
dq q2 + 8
∂P(x, q) ∂Q(x, q)
= 2q ≠ =q
∂q ∂x
297
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Dado que no se satisface el criterio de exactitud, ahora vamos a
buscar si existe un factor de integración adecuado para transformar
(134) en una ecuación diferencial exacta.
∂P (x , q ) ∂Q (x , q )
−
2q − q g ( x )dx
= h (x ) ⇒ µ (x ) = e ∫
∂q ∂x 1
= =
Q (x , q ) q (x + 12 ) x + 12
dx
∫
µ(x) = e x +12 =e
ln x +12
= x + 12 ⇒ como x ≥ 0 ⇒ µ(x) = x + 12
(1xq44
2
+ 84
x2+ 12
4q44 3)dx + (1
+ 96 2
x4q+
44 qx4+4
242 144
43q)dq = 0 2
(135)
P1(x, q) Q1(x, q)
∂P1(x, q) ∂Q1(x, q)
= 2qx + 24q = = 2qx + 24q
∂q ∂x
∂F(x, q)
= P1(x, q)
∂x
(136)
∂F(x, q)
= Q1(x, q)
∂q
∫
F(x, q) = ∂F(x, q) = ∫ P1(x, q)∂x + f (q) (137)
Reemplazando P1(x, q) en (137) resulta:
∫ (xq )
+ 8x + 12q 2 + 96 ∂x + f (q)
2
F(x, q) =
298
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
F(x, q) = (1
x42)q4
2
( (x 4
4+284 +4 3) + f (q)
24 ) (138)
g(x, t)
∂F ( x , q ) ∂f (q )
= x 2 q + 24qx + = Q 1 (x , q ) = x 2 q + 24 qx + 144 q
∂q ∂q
∂f (q)
= 144q (139)
∂q
F( x , q ) =
( )
x q 2 + 8 (x + 24 )
+ 72q 2 = C (141)
2
C = 25668 (142 )
F( x , q ) =
(
x q 2 + 8 (x + 24 )) + 72 q 2 = 25668
2
52488
x = −12 ±
q2 + 8
52488
x = −12 + ; q ∈ [0, +∞ ) (143)
q2 + 8
299
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Note que este problema también podemos resolverlo probando con:
∂Q(x, q) ∂P(x, q)
−
q − 2q −q f (q)dq
= f (q) ⇒ µ(q) = e ∫
∂x ∂q
= =
P(x, q) q +82 2
q +8
−qdq
∫ q 2 +8 −1 2 ln q 2 + 8 −1 2
µ(q) = e =e = q2 + 8 ⇒ como q 2 + 8 > 0
(
µ (q ) = q 2 + 8 )−1 2
=
1
q2 + 8
2 xq + 12 q
q + 8 dx + dq = 0 (144 )
14243 q2 + 8
142 4
P2 (q , x ) 4 43
Q 2 (q , x )
∂P2 ( x , q ) q ∂Q 2 ( x , q ) q
= = =
∂q 2
q +8 ∂x 2
q +8
∂F ( x , q )
= P2 ( x , q )
∂x
(145 )
∂F ( x , q )
= Q 2 (x, q)
∂q
300
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
F( x , q ) = ∫ q 2 + 8 ∂x + f (q )
F( x , q ) = q 2 + 8 x + f (q ) (147 )
144244 3
g (x ,t )
∂F ( x , q ) qx ∂f (q ) xq + 12 q
= + = Q 2 (x , q ) =
∂q q2 + 8 ∂q q2 + 8
∂f (q ) 12 q
= (148 )
∂q q2 + 8
f (q ) = 12 q 2 + 8 (149 )
Posteriormente, reemplazando (149) en (147) e igualando F(x, q) a
una constante C , la solución de (133) resulta:
F( x , q ) = q 2 + 8 x + 12 q 2 + 8 = C
F( x , q ) = q 2 + 8 (x + 12 ) = C (150 )
C = 162 2 (151)
Reemplazando (151) en (150) obtenemos:
F( x , q ) = q 2 + 8 (x + 12 ) = 162 2
162 2 52488
x= − 12 = −12 + ; q ∈ [0, +∞ )
q2 + 8 q2 + 8
301
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
3.- Modelo macroeconómico (Sydsæter 2005): Sean C(t), I(t), e Y(t) el
consumo, la inversión y la renta nacional de un país en el instante
“t”. Supongamos que:
C ( t ) + I ( t ) = Y ( t )
dC ( t )
I (t ) = k (152 )
dt
C (t ) = aY (t ) + b
Donde “a”, “b” y “k” son constantes positivas, con a < 1. Se pide
obtener la ecuación diferencial para Y(t) transformándola en una
EDO exacta a través de un adecuado factor de integración si se sabe
que Y(0) = Y0. Compare el resultado obtenido con el que obtendría
resolviendo la ecuación diferencial para Y(t) utilizando la fórmula
para resolver las EDOs lineales de primer orden con coeficientes
constantes y término fijo: ecuación (XXI). Una vez resuelta la
ecuación diferencial para Y(t) se pide determinar I(t). Finalmente, se
pide calcular el siguiente límite:
Y(t) b
lím si Y0 ≠
t → +∞I(t) 1− a
dC (t )
= aY ' (t ) (153)
dt
∂P ( Y , t ) ∂Q ( Y , t )
=0≠ = a −1
∂t ∂Y
302
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Dado que no se satisface el criterio de exactitud, ahora vamos a
buscar si existe un factor de integración adecuado para transformar
(156) en una ecuación diferencial exacta.
Probamos con:
∂P ( Y , t ) ∂Q ( Y , t )
−
1− a g ( Y )dY
∂t ∂Y = = h (Y ) ⇒ µ (Y ) = e ∫
Q( Y, t ) (a − 1)Y + b
En consecuencia, el factor de integración será:
1− a
∫ (a −1)Y + b dY − ln ( a − 1)Y + b 1
µ (Y ) = e =e =
(a − 1)Y + b
1
µ (Y ) = ±
(a − 1)Y + b
Aun cuando podemos trabajar con cualquiera de los factores de
integración que acabamos de calcular, por simplicidad vamos a
trabajar con el factor positivo. Multiplicando (156) por el factor de
integración µ (Y ) se tiene:
(ka )
dY + ({ 1) dt = 0 (157 )
1(4
a − 1)Y + b
42443 Q 1 (Y , t )
P1 ( Y , t )
∂P1 ( Y , t ) ∂Q 1 ( Y , t )
=0= =0
∂t ∂Y
∂F (Y , t )
= P1 (Y , t )
∂Y
(158 )
∂F (Y , t ))
= Q 1 (Y , t )
∂t
F( Y , t ) = ∂F ( Y , t ) = P1 ( Y , t ) ∂Y + f (t )
∫ ∫ (159 )
303
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando P1 ( Y , t ) en (159) resulta:
(ka )
F( Y , t ) = ∫ (a − 1)Y + b ∂Y + f (t )
ka
F( Y , t ) = ln (a − 1)Y + b + f ( t ) (160 )
a−
14144 424444 3
g(x , t )
∂F ( Y , t ) ∂f (t )
= = Q 1 (Y , t ) = 1
∂t ∂t
∂f (t )
=1 (161)
∂t
f (t ) = t (162 )
Subsiguientemente, reemplazando (162) en (160) e igualando
F( Y , t ) a una constante C , la solución de (156) resulta:
ka
F( Y , t ) = ln (a − 1)Y + b + t = C (163)
a −1
ka
C= ln (a − 1)Y0 + b (164 )
a −1
ka ka
F( Y , t ) = ln (a − 1)Y + b + t = ln (a − 1)Y0 + b (165 )
a −1 a −1
ka ka
ln (a − 1)Y + b + t = ln (a − 1)Y0 + b
a −1 a −1
ka (a − 1)Y + b (a − 1)Y + b a −1
ln = − t ⇒ ln = − t
a −1 (a − 1)Y0 + b (a − 1)Y0 + b ka
304
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
a −1 a −1
(a − 1)Y + b − t
(a − 1)Y + b − t
= e ka >0⇒ = e ka
(a − 1)Y0 + b (a − 1)Y0 + b
a −1 a −1
(a − 1)Y + b − t
[(a − 1)Y0 + b ] − t
b
= e ka ⇒ Y (t ) = e ka −
(a − 1)Y0 + b a −1 a −1
1− a
b t
ka b
Y (t ) = Y0 − e
+ (166 )
1− a 1− a
a −1 b
Y ' + Y = −
(167 )
ka ka
1− a 1− a
t t
− b ka b
Y (t ) = + Ce ka = + Ce ka (168 )
(a − 1) ka 1− a
b
C = Y0 − (169 )
1− a
1− a
b b t
ka
Y (t ) = + Y0 − e
(170 )
1− a 1− a
1− a
1− a b t
ka
Y ' (t ) = Y0 −
e
(171)
ka 1− a
305
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando (171) en (154) resulta:
1− a
b t
ka
I (t ) = (1 − a ) Y0 − e
(172 )
1− a
1− a
b b ka t
+ Y0 − e
Y (t ) 1 − a 1 − a
lím = lím
t → +∞ I ( t ) 1− a
t → +∞ b ka t
(1 − a ) Y0 − e
1 − a
1− a
Y (t ) 1 b b − t
ka
lím = ⋅ lím Y0 − e + 1
t → +∞ I (t )
1− a t → +∞ 1− a 1− a
1− a
1− a − t
0 < a <1 ∧ k > 0 ⇒ > 0 ⇒ si t → +∞ ⇒ e ka →0
ka
Por tanto:
Y (t ) 1
lím =
I ( t ) 1 − a
t → +∞
4. Ecuación de Bernoulli
Una ecuación diferencial del tipo Bernoulli tiene la siguiente forma general:
1
y = x 1− n ⇒ y =
n −1
(LXI )
x
x' 1
= ⋅ f (t ) + g (t ) (LXII)
n n −1
x x
306
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Derivando (LXI) respecto de “t” tenemos:
dy dy dx dy dx
= ⋅ ⇒ = (1 − n )x − n ⋅
dt dx dt dt dt
y'
= yf (t ) + g (t ) ⇒ y ' = (1 − n )yf (t ) + (1 − n )g (t )
1− n
Ejemplos:
1.- Modelo de crecimiento económico de Solow-Swan (Vinogradov, 1999): El
modelo de Solow-Swan es el modelo básico de crecimiento económico y
éste descansa en la siguiente ecuación diferencial, que explica la dinámica
del capital por unidad de trabajo efectivo:
k ' (t ) = sf (k (t )) − (n + g + δ )k (t ); k (0 ) = k 0 (173)
Dónde:
k' 1
α
+
α −1
(n + g + δ ) = s (175 )
k k
307
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Haciendo el siguiente cambio de variable:
1
y = k 1−α ⇒ y =
α −1
(176 )
k
dy dy dk dy dk
= ⋅ ⇒ = (1 − α )k − α ⋅
dt dk dt dt dt
y'
+ (n + g + δ )y = s ⇒ y ' + (1 − α )(n + g + δ )y = (1 − α )s (178 )
1− α
s
y (t ) = + Ce −[(1−α )(n +g + δ )]t (179 )
(n + g + δ )
Igualando (176) y (179) se tiene:
s
k 1−α = + Ce −[(1−α )( n + g +δ )]t
(n + g + δ )
1 (1 − α )
s
k (t ) = + Ce − [(1 − α )(n + g + δ )]t (180 )
(n + g + δ )
s
C = k 10− α − (181)
(n + g + δ )
Sustituyendo (181) en (180) se tiene:
1 (1 − α )
s s − [(1 − α )( n + g + δ )]t
k (t ) = + k 10− α − e (182 )
(n + g + δ ) (n + g + δ )
Note que gracias a que gracias a que n+g+δ > 0 y 1 − α > 0 , entonces:
308
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Por tanto:
1 (1 − α )
s
lím k (t ) = = k LP
t → +∞
(n + g + δ )
K ' = γ 1 bK α + γ 2 K (183)
K' γ2
α
− = γ1b (184 )
K K α −1
1
y = K1− α ⇒ y =
α −1
(185 )
K
dy dy dK dy dK
= ⋅ ⇒ = (1 − α )K −α ⋅
dt dK dt dt dt
y'
− γ 2 y = γ 1 b ⇒ y ' − (1 − α )γ 2 y = (1 − α )γ 1 b (187 )
1− α
γ 1b
y (t ) = Ce [(1−α )γ 2 ]t − (188 )
γ2
309
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Igualando (185) y (188) se tiene:
γ1b
K 1 − α = Ce [(1 − α )γ 2 ]t −
γ2
1 (1− α )
γ 1b
K (t ) = Ce [(1−α )γ 2 ]t − (189 )
γ 2
γ1b
C = K 10− α + (190 )
γ2
x 1' (t ) = a 11 (t )x 1 (t ) + a 12 (t )x 2 (t ) + K + a 1n (t )x n (t ) + b 1 (t )
x '2 (t ) = a 21 (t )x 1 (t ) + a 22 (t )x 2 (t ) + K + a 2 n (t )x n (t ) + b 2 (t )
(LXV)
M
x 'n (t ) = a n1 (t )x 1 (t ) + a n 2 (t )x 2 (t ) + K + a nn (t )x n (t ) + b n (t )
x 1' (t ) a 11 (t ) a 12 (t ) K a 1n (t ) x 1 (t ) b 1 (t )
'
x 2 (t ) = a 21 (t ) a 22 (t ) a 2 n (t ) x 2 (t ) b 2 (t )
⋅ + (LXVI)
M M M M M M
x 'n (t ) a n1 (t ) a n 2 (t ) K a nn (t ) x n (t ) b n (t )
1424 3 1444442444443 1 4v24 3 1 4v24 3
v
X ' (t ) A(t ) X (t ) b (t )
310
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Donde A (t ) es una matriz de dimensión n × n cuyos coeficientes variables en el
v
tiempo son continuos en a < t < b , y b (t ) es un vector de dimensión n × 1 cuyas
componentes son variable en el tiempo y continuas en a < t < b. El caso de
v v
coeficientes constantes surge como un caso particular en el que A (t ) = A y b (t ) = b
son constantes (no dependen del tiempo). El sistema (LXVI) será homogéneo si
v v v v
b (t ) = 0 y no homogéneo si b (t ) ≠ 0. La solución completa de (LXVI) viene dada por
la combinación de la solución complementaria y de la solución particular o solución
de equilibrio:
v v v
X (t ) = X c (t ) + X p (t ) (LXVII)
v
Donde X c (t ) es la solución de la parte homogénea de (LXVI), esto es, es la solución
v v v
de X ' (t ) = A (t )X (t ); y X p (t ) es la solución que fija a (LXVI).
v v v
En general, si cada uno de los vectores X 1 (t ), X 2 (t ), K , X n (t ) son soluciones de
(LXVI), también lo será su combinación lineal, esto es:
v v v v
X (t ) = c1X1 (t ) + c 2 X 2 (t ) + K + c n X n (t ) (LXVIII)
x 11 (t ) x 12 (t ) K x 1n (t )
x (t ) x 22 (t ) x 2 n (t )
[ v1 v2 vn
χ (t ) = X (t ) X (t ) K X (t ) = 21
M
] M M
(LXX)
x n1 (t ) x n 2 (t ) K x nn (t ) (n×n )
Siendo:
x11 (t ) x12 (t ) x1 n ( t )
v1 x (t ) v x 22 (t ) ; K ; X n (t ) =
v
x 2 n (t ) ⇒
X ( t ) = 21 ; X 2 (t ) =
M M M
x n1 ( t ) x n 2 (t ) x nn (t )
x 1 j (t )
vj x 2 j (t )
X (t ) = ( j = 1, 2 K , n ) (LXXI)
M
x nj (t )
311
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Al determinante de la matriz χ (t ) 6 se le conoce como Wronskiano, W (χ ), esto es:
v v v
W (χ ) = χ (t ) = X 1 (t ) X 2 (t ) K X n (t ) (LXXII)
x 1' (t ) a 11 a 12 K a 1n x 1 (t )
'
x 2 (t ) = a 21 a 22 a 2 n x 2 (t ) v v
M M ⋅ ⇒ X ' (t ) = A ⋅ X (t ) (LXXV)
M M M
x 'n (t ) a n1 a n 2 K a nn x n (t )
1 4v24 3 1 4 44 4 2 444 4 3 1
4 v2 43
X ' (t ) A X (t )
v v
Si A ≠ 0, el único punto de equilibrio7 es X* = 0 :
312
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Por otro lado, para sistemas lineales no homogéneos tales como (LXXIII), el
equilibrio puede encontrarse a partir de:
x1' (t ) a11 a12 K a1n x1* b1 0
'
x 2 (t ) = a 21 a 22 a 2 n x *2 b 2 0
⋅ + = ⇒ A ⋅
v* v v
X +b=0 (LXXVII)
M M M M M M M
x 'n (t ) a n1 a n 2 K a nn x *n b n 0
1 4v24 3 1444 424444 3 123 { v { v
v
X ' (t ) A X* b 0
v
v* v v* −1
v adjA ⋅ b
A ⋅ X = −b ⇒ X = A ⋅ b = ⇔ A ≠0 (LXXVIII)
A
Dónde:
v v v
D (t ) = X (t ) − X * (LXXX)
x 1 (0 )
v x (0 )
Donde X 0 = 2 .
M
x n (0 )
313
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
v v
Decimos que un vector X (t ) es solución de (LXXXIII), si X (t ) es derivable, satisface
el sistema de ecuaciones diferenciales y las condiciones iniciales.
p (λ ) = A − λI = 0 (LXXXIV)
v v v n
v
X (t ) = e At X 0 = Pe Ψt P −1 X 0 = ∑ cieλ t vi i (LXXXV)
i =1
Dónde:
t2 tn ∞ tk
e At = I + tA + A2 + K + An + K = ∑ Ak (t ∈ ℜ) (LXXXVI)
2! n! k =0 k!
v v v
P = [v1 v 2 K v n ](n × n ) (LXXXVII)
λ1 0 K 0 λk1 0 K 0
0 λ2 K K 0 λk2
⇒ ψ k =
K K
ψ= (LXXXVIII)
K K O K K K O K
0 0 K λ n 0 0 K λkn
314
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
∞ λk1 t k
∑
k = 0 k!
0 K
0
eλ1t
0 K 0
∞ ∞ λk2 t k
ψk tk K = 0
λ2t
eψ t = ∑ = 0 ∑ K
e K K
(LXXXIX)
k! K k =0 k! K K O K
k =0
K O K
∞ λk t k 0 0 K e λ n t
n
0 0 K ∑
k = 0 k!
Ahora vamos a verificar que la primera expresión a la derecha del signo de igualdad
de (LXXXV) es solución de (LXXXIII). Para ello vamos a derivar (LXXXVI)
respecto a “t”:
( ) = A + tA
d e At 2
+K+
t n −1
A n + K = A I + tA +
t2 2
A +K+
t n n −1
A
+ K
dt (n − 1)! 2! n!
( ) = Ae
d e At At
(XC)
dt
v'
X (t ) =
( v
d e At X 0
=
)
d e At v
X0
( ) (XCI)
dt dt
Ahora vamos a demostrar que la segunda expresión a la derecha del signo de igualdad
de (LXXXV) es solución de (LXXXIII).
315
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Teniendo en cuenta (LXXXVII) y (LXXXVIII), la igualdad dada por (XCIV) se
puede expresar de forma equivalente en forma matricial de la siguiente forma:
λ1 0 K 0
v v v v v v 0 λ2 K K
A ⋅ [v 1 v2 K v n ] = [v 1 v2 K vn ]⋅
K K O K
0 0 K λ n
Es decir:
AP = Pψ ⇒ ψ = P −1 AP o A = PψP −1 (XCV)
( )( ) ( )
A 2 = A ⋅ A ⇒ A 2 = PψP −1 ⋅ PψP −1 = Pψ P −1Pψ P −1 = Pψ (ψ )P −1 = Pψ 2 P −1
3 2
A = A ⋅A ⇒ A 3
= (Pψ P ) ⋅ (PψP ) = Pψ (P Pψ )P
2 −1 −1 2 −1 −1
= Pψ 2 (ψ )P −1 = Pψ 3 P −1
M
A k = P ψ k P −1 (XCVI)
∞ tk ∞ tk ∞ tk
e At = ∑ Ak = ∑ Pψ k P −1 =P ψ k P −1 = Peψt P −1 (XCVII)
∑
k! k!
k =0 k =0 k = 0 k!
Por último vamos a verificar que la tercera expresión a la derecha del signo de
igualdad de (LXXXV) es solución de (LXXXIII). Para ello, con el propósito de
desacoplar el sistema (LXXXIII) realizaremos el siguiente cambio de variables:
v v
X (t ) = PY (t ) (XCVIII)
316
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
de (LXXXVIII)
64 447444 8
y1' (t) λ1 0 K 0 y1(t) y1' (t) λ1y1(t)
' '
v v y (t) 0 λ2 K K y 2(t) y (t) λ y (t)
Y '(t) = ψY(t) ⇒ 2 = ⋅ ⇒ 2 = 2 2 (CI)
M K K O K M M M
y1(t) 0 0
' K λ n y n (t) y1(t) λ n y n (t)
'
y1 (t ) c1e λ1 t
λ2t
v y (t )
Y (t ) = 2 = 2
c e
(CII)
M M
y n (t ) c n e λ n t
c1eλ1t c1eλ1t
λ de (LXXXVII)
v c e 2 6 t v 44 v47444 v8 c 2eλ 2 t n
v
X(t) = P ⋅ 2 = [v v K v ] ⋅ M = ∑ ci eλ t v i
i
M
1 2 n
i =1
c n eλ n t c n eλ n t
Dónde:
v n v
v
X (0 ) = ∑ ci vi = X 0 (CIII)
i =1
Ejemplo:
Resolver el siguiente sistema:
x ' 1 1 x
' = ⋅
y − 2 4 y
(191)
v
()
X0 =
()
x 0 1
=
()
y 0 2
1− λ λ 1 = 3
()
p λ = A − λI =
1
−2 4−λ
= λ2 − 5 λ + 6 = 0 ⇒ (192)
λ 2 = 2
Para λ 1 = 3 :
v −2 1 a 0
[A − λ 1 I ]vv 1 =0⇒ ⋅ = ⇒ −2 a + b = 0 ⇒ b = 2 a (193)
− 2 1 b 0
317
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Si hacemos a = k ⇒ b = 2 k , entonces:
v k 1 v 1
v 1 = = k ⇒ si k = 1 ⇒ v 1 = (194)
2 k 2 2
Para λ 2 = 2 :
v −1 1 c 0
[A − λ 2 I ]vv 2 =0⇒ ⋅ = ⇒ −c + d = 0 ⇒ c = d (195)
− 2 2 d 0
Si hacemos c = d = k , entonces:
v k 1 v 1
v 2 = = k ⇒ si k = 1 ⇒ v 2 = (196)
k
1 1
−1
v 1 1 e 3 t 0 1 1 1 1 1 e 3 t 0 − 1 1 1
X (t ) = ⋅ ⋅
2 t 2 1
⋅ = ⋅ ⋅ ⋅
2 1 0 e 2 2 1 0 e 2 t 2 − 1 2
v e3 t 1 1
X (t ) = 3 t = e 3 t − 0e 2 t (197)
2e 2 1
Consideremos el caso en el que el sistema (LXXXIII) posee una matriz “A” que tiene
autovalores repetidos, digamos algún “ λ i = λ ” repetido “m” veces. En caso “m”
coincida con el número “ε” de autovectores linealmente independientes asociados a
“ λ i = λ ”, entonces la solución de (LXXXIII) se hallará con (LXXXV). En caso “m”
sea mayor a “ε”, la solución de (LXXXIII) se hallará con:
*Si m = 2 :
v
() v v
( v
X t = c 1 e λ t v 1 + c 2 e λ t tv 1 + v 2 ) (CIV)
*Si m = 3 :
v
() v v
( v
) v v v
X t = c 1 e λt v 1 + c 2 e λt tv 1 + v 2 + c 3 e λt t 2 v 1 + 2 tv 2 + 3 v 3
(CV)
Dónde:
(A − λI )vv i v
= v i −1 (i = 1, 2, K , m ) (CVI)
318
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Ejemplo:
Resolver el siguiente sistema:
x ' 1 1 x
' = ⋅
y − 1 3 y
(198)
v x 0 2
X0 =() ()
=
()
y 0 3
1− λ
()
p λ = A − λI =
1
−1 3− λ
= λ2 − 4 λ + 4 = λ − 2 ( )2 {
= 0 ⇒ λ1 = λ 2 = λ = 2 (199)
Para λ = 2 :
[A − λI ]vv 1 = 0 ⇒ −− 11 1 a 0
v
⋅ = ⇒ −a + b = 0 ⇒ b = a (200)
1 b 0
Si hacemos a = b = k , entonces:
v k 1 v 1
v 1 = = k ⇒ si k = 1 ⇒ v 1 = (201)
k
1 1
−1 1 c 1
[A − λI ]vv 2 v
= v1 ⇒ ⋅ = ⇒ − c + d = 1 (202)
− 1 1 d 1
v 0
v2 = (203)
1
()
x t 2 t 1
2 t 1 0
= c1 e + c 2 e t + (204)
()
y t 1
1 1
319
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Teniendo en cuenta las condiciones iniciales obtenemos:
()
x 0 2 1 0 c 1 = 2
= = c1 + c 2 ⇒ (205)
()
y 0
3
1
1 c 2 = 1
v v v
() ( )
X t = e αt h 1 cos βt + h 2 sen βt ( ) (CVII)
v v
Donde h 1 = c 1 v 1 + c 2 v 2 y h 2 = i (c 1 v 1 − c 2 v 2 ) son vectores de componentes reales.
v v v v
Ejemplos:
Resolver el siguiente sistema:
x ' − 3 4 x
' = ⋅
y − 2 1 y
(207)
v x (0 ) 1
X (0 ) = =
y (0 ) 4
−3 4 −3 − λ 4
A= ; A = −3 + 8 = 5; A − λI = = λ2 + 2 λ + 5 = 0
− 2 1 − 2 1 − λ
α = −1
λ 1 = −1 + 2i y λ 2 = −1 − 2i ⇒ (208)
β = 2
−2 − 2i 4 a 0
[A − λ 1 I ]vv 1 = ⋅ = ⇒ a = b (1 − i )
−2 2 − 2i b 0
Si hacemos a = 2 ⇒ b = 1 + i , entonces:
v 2
v1 =
1 + i
Para λ 2 = −1 − 2i :
−2 + 2i 4 c 0
[A − λ 2 I ]vv 2 = ⋅ = ⇒ c = d (1 + i )
−2 2 + 2i d 0
Si hacemos c = 2 ⇒ d = 1 − i, entonces:
v 2
v2 =
1 − i
v v v
X (t ) = e − t h 1 cos (2 t ) + h 2 sen (2 t ) (209)
Con:
v v v 2 2 2c 1 + 2 c 2
h 1 = c1 v1 + c 2 v 2 = c1 + c2 = (210)
1 + i 1 − i (1 + i )c 1 + (1 − i )c 2
v v v 2 2 2i (c 1 − c 2 )
h 2 = i (c 1 v 1 − c 2 v 2 ) = ic 1 − ic 2 = (211)
1 + i 1 − i 1i [c (1 + i ) − c 2 (1 − i )]
1 − 7i
c1 =
v 1 v 2c1 + 2c 2
X (0 ) = = h 1 = ⇒
4 (212)
4 (1 + i )c1 + (1 − i )c 2 c = 1 + 7i
2 4
v 1 v 7
h1 = y h 2 = (213)
4
3
321
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Reemplazando (213) en (209) se tiene:
Modelo Keynesiano IS-LM (Tu, 1994): En este modelo la renta “Y” responde al
exceso de la cantidad demandada (esto es, al exceso de la inversión “I” respecto al
ahorro “S”), y la tasa de interés “r” responde al exceso de la demanda de dinero
(preferencia de liquidez) “ L (Y , r ) ” respecto a la oferta de dinero “ M ” determinada
exógenamente. El modelo matemáticamente se puede expresar de la siguiente forma:
Dónde:
(
Y ' = I 0 − µr − s Y − T − T − G ) ( )
' (216)
r = θY − γr − M
Y ' − s − µ Y I 0 + (s − 1)T + G
'= ⋅ + (217)
1r θ − γ { r −M
2 3 14243 v 1 4 4 4 2v 4 4 43
v A
X' X b
Para transformar el sistema (217) en uno homogéneo vamos a determinar el punto de equilibrio
v*
X del sistema (217):
Y ' − s − µ Y * I 0 + (s − 1)T + G 0
'= ⋅ + =
θ42−
r 1 4
γ r *
3
−M 0
1 2v3 12 v3
14442 v 4443
X' A X* b
322
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
− γ µ − I0 − (s − 1)T − G
⋅
Y* − s − µ −1 − I0 − (s − 1)T − G − θ − s M
=
* ⋅ =
{ 1
r θ
42 − γ M
44424443 − s − µ
v 4 43 4 1
X* A −1 (−bv) θ −γ
v * Y * 1
[
µM + γ I 0 + (s − 1)T + G
]
X = *=
r sγ + µθ
[
θ I 0 + (s − 1)T + G − sM
]
[
µM + γ I 0 + (s − 1)T + G ] [
* µM + γ I 0 + (s − 1)T + G
Y =
]
v * Y * sγ + µθ sγ + µθ
X = *= ⇒ (218)
r
[
θ I 0 + (s − 1)T + G − sM ]
[
θ I 0 + (s − 1)T + G − sM ]
r * =
sγ + µθ sγ + µθ
v d v v Y (t ) − Y *
D (t ) = 1 = X (t ) − X * (t ) = *
(219)
d 2 r (t ) − r
v v d ' − s − µ d1
D ' (t ) = AD (t ) ⇒ 1' = ⋅ (220)
d 2 θ − γ d 2
Ahora podemos resolver el sistema homogéneo dado por (220). Para ello vamos a
calcular el polinomio característico:
−s − λ −µ
p (λ ) = = λ2 + (s + γ )λ + (sγ + µθ ) = λ2 − (trA )λ + A = 0 (221)
θ −γ−λ
Dónde:
A = sγ + µθ > 0
trA = − (s + γ ) < 0
323
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Ahora vamos a asignar valores a los parámetros del problema de manera que podamos
encontrar soluciones explícitas para los tres casos analizados en la sección anterior.
λ 1 ∈ ℜ
Caso 1: ∆ = (trA )2 − 4 A > 0 ⇒ (trA )2 > 4 A ⇒ λ 2 ∈ ℜ
λ ≠ λ
1 2
trA = −3,3
A = 0,9625
∆ = 7,04 (223)
*
Y = 4,3117
*
r = 0,02597
− 3,3 + 7,04
λ1 = ≈ −0,3234
2
− 3,3 − 7,04
λ2 = ≈ −2,9766
2
Para λ 1 :
v 10,7
v1 =
1
324
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Para λ 2 :
v 1
v2 =
10,7
Por tanto:
c1 ≈ −0,22477
(227)
c 2 ≈ 0,09335
325
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Caso II: ∆ = (trA ) − 4 A = 0 ⇒ (trA )2 = 4 A ⇒ {λ 1 = λ 2 = λ ∈ ℜ
2
trA = −1
A = 0,25
∆ = 0 (228)
*
Y = 1,84
*
r = 3,28
{λ1 = λ 2 = λ = −1 2
Para λ 1 = λ 2 = λ :
Si hacemos a = −1 ⇒ b = 3, entonces:
v −1
v1 =
3
v 0
v2 = (231)
10
326
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
En consecuencia, de (65) tenemos:
1 1
d1 − t − 1 − t − 1 0
2 t
= c1e + c 2e +
2 (232)
d 2 3 3 10
1
− t
1
− t
d1 − c e 2 − tc e 2
= 1
1
1
2
1 (233)
d 2 − t − t − t
3c1e 2 + 3tc2e 2 + 10c 2e 2
1
− t
1
− t
Y(t) − Y* − c e 2 − tc e 2
*
= 1 2
1
r(t) − r
1 1
− t − t − t
3c1e 2 + 3tc 2e 2 + 10c 2e 2
1
− t
1
− t
Y(t) Y *
− c e 2 − tc e 2
= 1 2
1
r(t) *
1 1
− t − t − t
r + 3c1e 2 + 3tc 2e 2 + 10c 2e 2
1
− t
1
− t
Y(t) = Y* − c e 2 − tc e 2
1 2
(234)
1 1 1
− t − t − t
r(t) = r* + 3c1e 2 + 3tc 2e 2 + 10c 2e 2
1
− t
1
− t
Y(t) = 1,84 − c e 2 − tc e 2
1 2
(**)
1 1 1
− t − t − t
r(t) = 3,28 + 3c1e 2 + 3tc 2e 2 + 10c 2e 2
c1 = −0,16
(235)
c 2 = −0,2
1
− t
1
− t
Y(t) = 1,84 + 0,16e 2 + 0,2te 2
1 1 1
− t − t − t
r(t) = 3,28 − 0,48e 2 − 0,6te 2 − 2e 2
327
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
λ1 = α + βi
Caso III: ∆ = (trA )2 − 4 A < 0 ⇒ (trA )2 < 4 A ⇒
λ 2 = λ1 = α − βi
trA = −1
A = 1
∆ = −3 (236)
*
Y = 1,5
*
r = 1
1 3 1 3 α = − 1 2
λ1 = − + i y λ2 = − − i⇒ (238)
2 2 2 2 β = 3 2
1 3
Para λ1 = − + i:
2 2
3
Si hacemos b = 1 ⇒ a = i, entonces:
2
3
v i
v1 = 2
1
328
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
1 3
Para λ 2 = − − i:
2 2
3 2 i − 0,75
c 0 3
[A − λ 2I]vv 2 =
⋅ = ⇒ c = − id
1
d
3 2 i 0 2
Si hacemos d = 1 ⇒ c = − 3 2 i, entonces:
3
v − i
v2 = 2
1
− t v
1
d1 3 v 3
2 h cos
t
=e t + h 2 sen (239)
1
d
2 2 2
Con:
3
v
i(c − c )
v v 1
h1 = c1v1 + c 2 v 2 = 2
2
(240)
c1+ c 2
3
v −
(c + c )
v v
h 2 = i(c1v1 − c 2 v 2 ) = 2
1 2
(241)
i(c1 − c 2 )
Reemplazando (240) y (241) en (239) tenemos:
t
1
−
d1 = e 2 3 (c − c ) i cos 3 t − 3 (c + c ) sen 3 t
2 1 2
2 1 2
2
2
(242)
t
1
−
d2 = e 2 (c + c ) cos 3 t + (c − c ) i sen 3 t
1 2
2 1 2
2
t
1
−
Y(t) = Y* + e 2 3 (c − c ) i cos 3 t − 3 (c + c ) sen 3 t
2 1
2
2 1 2
2 2
(243)
t
1
−
r(t) = r* + e 2 (c + c ) cos 3 t + (c − c ) i sen 3 t
1
2
2 1 2
2
329
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
* *
Reemplazando Y y r en el sistema (243) tenemos:
t
1
−
Y(t) = 1,5 + e 2 3 (c − c ) i cos 3 t − 3 (c + c ) sen 3 t
2
2
1 2 1 2
2 2
(***)
t
1
−
r(t) = 1 + e 2 (c + c ) cos 3 t + (c − c ) i sen 3 t
1 2
2 1 2
2
Y(0) 2 1,5 + 3 (c − c ) i
=
= 1 2
r(0) 0,8 1 +2c + c
1 2
− 0,6 − 3 i
c1 =
6
(244)
− 0,6 + 3 i
c 2 =
6
t 1
1
− 3
Y(t) = 1,5 + e 2 cos 3
t +
3
sen
t
2 10 2
2
(245)
t
1
− 2 3 3 3
r(t) = 1 + e 2 − cos t + sen t
10 2 3 2
330
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Diagramas de fase y soluciones cualitativas
Aunque muchas ecuaciones diferenciales ordinarias (EDO) no lineales de primer
orden, que suelen aparecer con relativa frecuencia en modelos económicos, no pueden
resolverse analíticamente, las propiedades cualitativas de sus soluciones pueden
algunas veces ser descritas gráficamente a través de los denominados diagramas de
fase. Supongamos que tenemos la forma general de una ecuación diferencial
autónoma8 de primer orden (lineal o no lineal respecto a “y”):
dy
= f (y ) (CVIII)
dt
En esta sección estudiaremos algunos de los posibles diagramas de fase de una única
variable que pueden surgir en el análisis cualitativo de ecuaciones diferenciales
autónomas de primer orden, aquellas que satisfacen la ecuación (CVIII). Asociados a las
curvas de fase representadas en las figuras 10, 11,…,14, y 15 tenemos sus respectivas
trayectorias, sendas u órbitas temporales (figuras 16, 17,…, 20, y 21).
dy dt
f (y )
y 0B
y
y E y 0A
Figura 10
8
Una ecuación diferencial ordinaria (EDO), tal como la ecuación (CVIII), en la que el tiempo “t” no
aparece explícitamente como argumento de “f”, se denomina autónoma.
331
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
dy dt f (y )
yE
y
y 0B y 0A
Figura 11
dy dt
D
C E
f(y)
B F
y y0 y1 y2 y
y
−
A G
H
Figura 12
dy dt f (y )
yE
y
y 0B y 0A
Figura 13
332
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
dy dt
f (y )
yE
y
y 0B y 0A
Figura 14
dy dt
f (y )
y
yE
Figura 15
Ahora veremos cómo, una vez conocida la curva de fase correspondiente a la ecuación
diferencial (CVIII), obtener información cualitativa importante de la trayectoria temporal
y(t ).
333
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
En las figuras 10, 11, 12, 13, 14 y 15 se puede observar que encima del eje horizontal el
signo de dy dt es positivo, esto es dy dt > 0, por lo que la variable “y” aumentará con el
transcurso del tiempo. Es por esta razón que los puntos que se encuentran encima del eje “y”
en las curvas de fase de las figuras 10, 11, 12, 13, 14 y 15 se desplazan de izquierda a
derecha, tal como muestran las flechas dibujadas sobre dichas curvas. De manera análoga,
podemos observar que los puntos que se encuentran debajo del eje “y” en las curvas de fase
se desplazarán de derecha a izquierda, ya que en esta dirección la dy dt < 0, por lo que “y”
disminuirá con el paso del tiempo. Asimismo, es importante resaltar que los resultados antes
vistos son independientes del signo de la variable “y”. Es decir, aun cuando los diagramas de
fase de las figuras 10, 11, 12, 13, 14 y 15 se encontrasen a la izquierda del eje vertical
(y < 0 ), el sentido de las flechas en dichas curvas de fase no se vería afectado.
La estabilidad dinámica del equilibrio: la convergencia de la trayectoria
temporal
y 0A
yE
y 0B
Figura 16
y (t )
y 0A
yE
y 0B
t
Figura 17
En la figura 12 podemos apreciar una curva de fase que forma un “lazo cerrado”. Esto
se debe a que f (y ) es una relación. Cuando el diagrama de fase es un lazo cerrado, un
movimiento oscilatorio de amplitud constante ocurrirá siempre que se verifiquen dos
condiciones: i) una parte de la curva de fase debe permanecer encima y otra parte
debajo del eje “y”, de modo que haya una etapa en la que “y” crece y otra etapa en la
que “y” decrece, ii) la curva de fase debe presentar pendiente no definida en los
puntos de intersección con el eje “y”. En la figura 12 se puede apreciar que en los
puntos “B” y “F” de la curva de fase se cumple que la dy dt = 0, y sin embargo estos
puntos no representan puntos de equilibrio. Más bien, son extremos (mínimos y
máximos) relativos de la trayectoria temporal correspondiente a esta curva de fase, tal
como se aprecia en la figura 18. En esta figura, se aprecia que la trayectoria temporal
de “y” oscila periódicamente entre el valor correspondiente al punto B (mínimo
relativo) y el valor correspondiente al punto F (máximo relativo).
De los casos antes vistos podemos concluir que una condición necesaria, pero no
suficiente, que debe satisfacer todo punto de equilibrio y E es que la dy (t E ) dt = 0.
335
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Figura 18
Tabla I
Del análisis realizado en esta sección podemos observar que la estabilidad dinámica
del equilibrio, o equivalentemente la convergencia de la trayectoria temporal y (t ),
dependerá del signo que tenga la pendiente de la curva de fase en la vecindad al punto
de intersección de dicha curva con el eje “y”. Por ejemplo, en la figura 10 se aprecia
que en la vecindad al punto y E la curva de fase presenta una d[f (y )] dy > 0, lo cual
da lugar a la inestabilidad dinámica. Mientras que en la figura 11 se aprecia que en la
vecindad al punto y E la curva de fase presenta una d[f (y )] dy < 0, lo cual da lugar a
la estabilidad dinámica de la variable “y”. Sin embargo, en la figura 12 se aprecia que
en los puntos “B” y “F”, que como ya hemos dicho no son puntos de equilibrio (son
sólo las cotas de una trayectoria temporal fluctuante), la curva de fase presenta una
d [f ( y )] dy que no está definida.
336
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Asimismo, si la curva de fase es tangencial al eje horizontal, permaneciendo siempre
a un lado de este (es decir, si el punto de tangencia es un máximo o un mínimo de la
curva de fase), el punto de equilibrio es estable por uno de sus lados e inestable por el
otro9. Este es el caso del punto y E de la figura 13, estable por izquierda (para valores
menores a y E ) e inestable por derecha (para valores mayores a y E ). En la figura 19
observamos las posibles trayectorias temporales correspondientes a la curva de fase
de la figura 13. Se observa que si y (0 ) = y 0 < y E entonces la trayectoria temporal de
“y” converge en el largo plazo hacia el valor y E . Sin embargo, si y (0 ) = y 0 > y E
entonces la trayectoria temporal de “y” diverge del valor y E (rama de la trayectoria
temporal que se encuentra a la izquierda de la asíntota vertical y encima del valor y E ).
Además, si el punto de tangencia corresponde a un punto de inflexión horizontal de la
curva de fase (figura 14), entonces el punto de equilibrio es estable (inestable) si la
curva de fase permanece encima (debajo) del eje “y” a la izquierda (derecha) del punto
de inflexión. En la figura 20 podemos ver la trayectoria temporal correspondiente al
diagrama de fase de la figura 14.
Por otro lado, es importante resaltar que pueden existir movimientos oscilatorios
(convergentes o divergentes) de amplitud no constante cuando f (y ) es una relación. En
este caso la curva de fase no será un lazo cerrado como la figura 12 (aunque seguirá
satisfaciendo las condiciones i) y ii) de ocurrencia de un ciclo de la página 334) sino
una espiral (divergente) como muestra la figura 15. La posible trayectoria temporal
correspondiente a este diagrama de fase aparece en la figura 21.
9
A este tipo de punto se le denomina shunt. Un shunt es un punto de equilibrio alrededor del cual el
sistema evoluciona sin invertir su sentido (las flechas del eje “y” apuntan en un mismo sentido).
337
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
y (t )
y 0A
yE
y 0B
Figura 19
y (t )
y 0A
yE
y 0B
t
Figura 20
Figura 21
338
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
dk (t ) dt
t
k 1* = 0 k *2 = a
Figura 22
Ejemplos:
Realice un análisis cuantitativo y/o cualitativo de las siguientes ecuaciones
diferenciales:
P ' = dP dt = k [a + bP − α − βP ] = k (a − α ) + k (b − β )P (246)
P ' − k (b − β )P = k (a − α ) (247)
339
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Para realizar el análisis cualitativo de la ecuación diferencial necesitamos construir
el diagrama de fase. Dado que la ecuación (246) representa una recta en el plano
P ' -P, necesitaremos dos puntos para poder construir dicha recta. Para ello,
primero determinaremos el punto de equilibrio, intersección con el eje “P”,
igualando la ecuación (246) a cero:
α−a
P ' = k (a − α ) + k (b − β )PE = 0 ⇒ PE =
b−β
P ' = k (a − α ) > 0
dP dt
k (a − α )
PE P0A
P
P0B
Figura 23
r − k (b − β ) = 0 ⇒ r = k (b − β )
PC = Ae k ( b −β )t
α−a
− k (b − β )PP = k (a − α ) ⇒ PP = (248)
b−β
α−a
P (t ) = Ae k (b − β )t + (249)
b −β
α−a α−a
P0 = A + ⇒ A = P0 −
b−β b−β
α − a k ( b − β )t α−a
P (t ) = P0 − e + (250)
b − β b−β
α−a
lím P (t ) = = PE > 0.
t →∞ b −β
P (t )
P0A
PE
P0B
t
Figura 24
dY
= k (D − Y ) (251)
dt
D = C + I0 + G 0
D = C (Y ) + I 0 + G 0
dY
Y' = = k [C (Y ) + I 0 + G 0 − Y ] ≡ kf (Y ) (252)
dt
Donde f (Y ) = C (Y ) + I 0 + G 0 − Y .
342
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Para obtener el diagrama de fase de la ecuación (252), primeramente vamos a
calcular la pendiente de la curva de fase en el plano Y ' -Y10:
dY ' df (Y )
= k [dC (Y ) dY − 1] ≡ k ⋅ (253)
dY dY
df (Y ) dC (Y )
Donde = − 1.
dY dY
d2Y' d 2 C (Y ) d 2 f (Y )
=k⋅ ≡ k⋅ (254)
dY 2 dY 2 dY 2
d 2 f (Y ) d 2 C (Y )
Donde = .
dY 2 dY 2
10
No tiene sentido trabajar con valores negativos de “Y” ya que la renta es una variable económica. Por
tanto, estaremos interesados en averiguar el signo de la pendiente de la curva de fase únicamente en el
primer y cuarto cuadrantes.
343
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Por otro lado, si la función de consumo es estrictamente cóncava respecto a la
renta, d 2 C (Y ) dY 2 < 0, por (254) la curva de fase también será estrictamente
cóncava respecto a la renta, d 2 f (Y ) dY 2 < 0. Además, si f (Y ) es estrictamente
decreciente con la renta, entonces df (Y ) dY = dC (Y ) dY − 1 < 0. Esto significaría
que la propensión marginal al consumo 0 < dC (Y ) dY < 1. Asimismo, al ser el
intercepto de la curva de fase con el eje “ Y ' ” k ⋅ f (0 ) = k ⋅ [C (0 ) + I 0 + G 0 ] > 0, y
al ser df (Y ) dY < 0 y d 2 f (Y ) dY 2 < 0, habrá algún valor de “Y” perteneciente al
intervalo [0, ∞ ) en el que f (Y ) = 0, con lo cual se garantizará que Y ' = 0.
f (Y E ) = C (Y E ) + I 0 + G 0 − Y E = 0
YE = f −1 (0 )
dY dt
kf (Y ) : f ' (Y ) > 0
kf (0 )
YE = f −1 (0 )
Y
kf (Y ) : f ' (Y ) < 0
Figura 25
344
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Plano de fase y retratos de fase de sistemas dinámicos autónomos
dx
= x ' = f (x , y )
dt
(CIX)
dy = y ' = g (x , y )
dt
v v x x (t )
Donde X = X (t ) = = es una solución del sistema (CX)12. Las componentes
y
y ( t )
v
de X (t ) pueden entenderse como un par de ecuaciones paramétricas de forma que
v v
para cada instante “t” se tiene un punto X = X (t ) ∈ ℜ 2 . Las ecuaciones paramétricas
x (t ) y y (t ) son funciones diferenciables respecto a “t” que satisfacen el sistema (CX)
v
sobre algún intervalo abierto “I”. Una solución X (t ) describe una curva o senda
(curva o senda de fase) en el plano x-y (plano de fase) que consta de todos los puntos
{(x (t ), y (t )) t ∈ I}. Al conjunto de todas las posibles sendas de fase13 se le denomina
retrato de fase.
11
Un sistema de ecuaciones diferenciales se denomina autónomo cuando la variable “t” no aparece
explícitamente en el sistema de ecuaciones; en caso contrario se dice que el sistema es no autónomo.
12 v
Téngase en cuenta que si X (t ) es solución de (CIX), al ser (CIX) y (CX) sistemas equivalentes, también
será solución de (CX). Por esta razón, sólo haremos referencia a uno de los dos sistemas, el sistema (CX),
en el resto de esta sección.
13
Es importante resaltar que la senda de fase no debe ser confundida con el diagrama de fase analizado en
la sección precedente. El análogo a la senda de fase para una única ecuación diferencial sería el eje “y”.
345
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Es decir, el campo vectorial de (CX) es un conjunto de vectores en el plano x-y, tal
v
que la pendiente del vector F(x , y ) en el punto (x , y ) coincide con la pendiente de la
tangente a la senda de fase que pasa por el punto (x , y ).
v
La pendiente del vector F(x , y ) en el punto (x , y ) viene dada por:
dy dy dt y'
= = ⇔ x' ≠ 0 (CXII)
dx dx dt x'
En la figura 26 se aprecia una senda de fase que es una solución particular del sistema
(CX), tal que en el instante inicial t = 0 debe satisfacer la condición inicial
v
X 0 = ( x (0 ), y (0 )) = (x 0 , y 0 ). En esta figura podemos apreciar que en el punto inicial
v
X0 los signos de las razones de cambio de “x” y “y” son
v
( ) v
( )
f X 0 = x ' (0 ) > 0 y g X 0 = y ' (0 ) > 0, entonces según se incremente el tiempo, el
v
sistema se moverá desde el punto X 0 hacia la derecha y hacia arriba. Asimismo, se
v
aprecia que en un punto genérico en el instante “t” tal como X = (x (t ), y (t )) los signos
de las razones de cambio de “x” y “y” son f (x , y ) = x ' (t ) > 0 y g (x , y ) = y ' (t ) < 0,
v
entonces según transcurra el tiempo, el sistema se moverá desde el punto X hacia la
v
derecha y hacia abajo. En el punto X = (x (t ), y (t )) se aprecia que la velocidad de
movimiento (razón o tasa de cambio instantánea respecto al tiempo) está dada por la
v v
longitud del vector F X . ( )
y
v
v X = (x (t ), y (t ))
X 0 = (x (0 ), y (0 ))
y ' (t ) v
( )
F(x , y ) = x ' (t ), y ' (t )
x (t )
'
Figura 26
14
Una órbita, a diferencia de una senda de fase, no nos da información sobre el sentido del movimiento.
346
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Para ilustrar la dinámica del sistema (CX), en principio, podemos dibujar tales
vectores en cada punto del plano x-y. Tal familia de vectores es llamada campo
vectorial15. En la práctica podemos dibujar sólo una pequeña muestra representativa
de estos vectores. Sobre la base del campo vectorial, teniendo en cuenta que los
vectores del campo vectorial son tangentes a las sendas de fase, podemos dibujar las
sendas de fase para el sistema y por tanto exhibir el retrato de fase del sistema.
x ' = x + 2 y = f (x , y ) v x ' 1 2 x
' ⇒ X ' = ' = ⋅ (CXIII)
y = 3 x + 2 y = g (x , y ) y
1 3 2 y
23 { v
A X
Para obtener algunos vectores del campo vectorial de este sistema podemos escoger
v
arbitrariamente algunos puntos del plano x-y para luego reemplazarlos en el vector X
que aparece en la ecuación (CXIII).
v x 3 v x ' 1 2 3 3
X = = ⇒ X ' = ' = ⋅ =
y 0 y 3
2 0 9
v x 4 v ' x' 1 2 4 8
X = = ⇒ X = ' = ⋅ =
y 2
y
3 2 2 16
v x 0 v x' 1 2 0 4
X = = ⇒ X ' = ' = ⋅ =
y 2
y
3 2 2 4
v x − 1 v x ' 1 2 − 1 3
X = = ⇒ X ' = ' = ⋅ =
y 2 y 3
2 2 1
v x − 3 v x' 1 2 − 3 − 3
X = = ⇒ X ' = ' = ⋅ =
y 0
y
3 2 0 − 9
v x − 1 v x' 1 2 − 1 − 5
X = = ⇒ X ' = = ⋅ =
y − 2
y'
3 2 − 2 − 7
v x 0 v x ' 1 2 0 − 10
X = = ⇒ X ' = = ⋅ =
y − 5 y' 3
2 − 5 − 10
v x 2 v x ' 1 2 2 − 8
X = = ⇒ X ' = = ⋅ =
y − 5 y' 3
2 − 5 − 4
Podemos observar que, por ejemplo, en el punto (3, 0 ) del plano x-y tendremos un
vector que apunta en la dirección (3,9 ). En el punto (2, −5 ) tendremos un vector que
apunta en la dirección (−8, −4 ). Repitiendo el proceso anterior para un gran número de
puntos del plano x-y se obtendrá el campo vectorial de la figura 2716.
15
Para obtener el campo vectorial de sistemas de ecuaciones diferenciales autónomos podemos utilizar
algunos programas matemáticos como Derive, Maple, Matlab, etc.
16
Es importante resaltar que la longitud de los vectores en la figura 27 ha sido proporcionalmente reducida
de manera que los vectores no interfieran el uno con el otro. Pero la longitud de cada vector aún sugiere la
velocidad de movimiento.
347
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Figura 27
Una vez que se ha obtenido el campo vectorial, podemos bosquejar algunas sendas de
fase teniendo en cuenta que los vectores de la figura 27 son tangentes a las sendas de
fase y que la dirección de los vectores nos da la dirección de la senda según se
incrementa el tiempo. Con el propósito de mostrar la dependencia temporal de la
solución, se han colocado flechas sobre las sendas de fase. En la figura 28 se muestra
el retrato de fase correspondiente al campo vectorial de la figura 27. Usualmente los
retratos de fase sólo incluyen algunas de las sendas de fase y no la totalidad de ellas.
Asimismo, frecuentemente el retrato de fase no va acompañado del campo vectorial,
aunque algunas veces por cuestiones didácticas se presentarán ambos bosquejados en
el plano x-y.
Figura 28
348
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Puntos fijos y estabilidad
v
(
Dado el sistema (CIX), si X * = x * , y * ) es un punto en el plano para el cual
simultáneamente f (x , y ) = 0 y g (x , y ) = 0 , entonces resulta que x ' = 0 y y ' = 0. Esto
significa que ni “x” ni “y” cambian con el tiempo: el sistema tiene un punto fijo, o
tiene un punto de equilibrio. Una vez que hemos encontrado un punto de equilibrio, lo
que nos interesa es saber si este punto es estable o inestable.
ε
( )
v
Bδ X*
v v
X0 X*
( )
v
Bε X*
Figura 29
349
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Por lo antes dicho, podemos introducir la siguiente definición de estabilidad.
v
( )
Definición I: Un punto de equilibrio X* = x * , y * es estable si para cualquier δ > 0
v v v v
existe un ε > 0 tal que si X 0 − X * < ε entonces X (t ) − X * < δ para todo “t”.
ε v
( )
Bδ X *
δ
v
X*
( )
v
Bε X *
v
X0
x
Figura 30
v
Note que la senda de fase puede alejarse del punto X * mientras permanece dentro de
( )
v
la bola Bε X * y aproximarse al punto de equilibrio en el límite. Un punto que es
estable pero que no es asintóticamente estable suele denominársele como neutral o
marginalmente estable. La figura 30 muestra un punto de equilibrio neutralmente
estable.
350
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
La estabilidad asintótica es más fuerte que la estabilidad. Esto es claro ya que si un
punto de equilibrio es asintóticamente estable, entonces debe ser estable. La condición
v
límite por sí sola no es suficiente ya que un sistema puede empezar “cerca de” X* (es
decir, dentro de la bola de radio “δ”) y aproximarse al punto de equilibrio en el límite,
pero divergir considerablemente (ir más allá) de bola de radio “ε” en determinado
periodo.
a) Es estable;
v v
b) Existe un α>0 tal que siempre que X 0 − X* < α entonces
v v
lím X (t ) − X * = 0.
t →∞
v
Si un sistema tiene un punto de equilibrio X* que es asintóticamente estable, y si
cada senda de fase se aproxima al punto de equilibrio (es decir, tanto para puntos
cercanos al punto de equilibrio como para puntos lejanos de éste), entonces el punto
de equilibrio se dice que es globalmente estable. Otra forma de considerar esto es
establecer el conjunto inicial de condiciones para las cuales el punto de equilibrio
dado sea asintóticamente estable, esto es, la bola más grande a partir de la cual
cualquier trayectoria entrante converja asintóticamente al punto de equilibrio. Este
conjunto de condiciones iniciales es llamado fuente de atracción. Un punto de
equilibrio es localmente asintóticamente estable si existe una fuente de atracción,
( )
v
Bε X * , dentro de la cual todas las sendas de fase entrantes a esta bola eventualmente
v
se aproximan al punto X* . Si la fuente de atracción es todo el plano x-y, entonces el
sistema es globalmente asintóticamente estable sobre el punto de equilibrio X* . (v )
Por lo antes dicho, podemos introducir la siguiente definición de estabilidad.
v
Definición III: Sea X* un punto de equilibrio asintóticamente estable del sistema
v
( )
v v
t →∞
v
(CIX), entonces el conjunto B X = X ∈ ℜ 2 lím X (t ) − X * = 0 es el dominio o
B(X ) = ℜ 2
v v
fuente de atracción de X. Si (o, en todo caso, si coincide con el plano de
v
fase) entonces se dice que X es globalmente asintóticamente estable. Si la estabilidad
v
únicamente se mantiene en una vecindad de X , se dice que es localmente
asintóticamente estable.
Algunas propiedades de las sendas de fase de sistemas autónomos son las siguientes:
1.- No más de una senda de fase pasa por un punto del plano x-y;
2.- Una senda de fase que empieza en un punto que no es un punto de equilibrio
únicamente alcanzará un punto de equilibrio en un periodo infinito;
3.- Ninguna senda de fase puede atravesarse así misma a menos que sea una curva
cerrada. Si la senda de fase es una curva cerrada entonces la solución es periódica.
351
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Isoclinas y líneas de fuerza en el plano de fase:
x (t ) x' f (x , y )
v
Sea X =
x v
= . Dado el sistema X ' = = ( )
v v
= F X , donde F X
y ' g (x , y )
( )
v v
y (t ) y
puede ser lineal o no lineal, los lugares geométricos de ℜ2 tales que
v x ' f (x , y ) a
X ' = ' = = , donde “a” y “b” son constantes, se denominan isoclinas.
y g (x , y ) b
Las isoclinas son de gran utilidad ya que nos permiten obtener la dirección de los
vectores del campo vectorial del sistema, lo que a su vez nos sirve para bosquejar las
sendas de fase.
En particular, las isoclinas en las que “a” y “b” son nulas (ceroclinas),
v x ' f (x , y ) 0
X ' = ' = = , nos dan los puntos para los que ya no hay ajuste
y g (x , y ) 0
dinámico para “x” y para “y”. Es decir, en las intersecciones de las ceroclinas
encontraremos los puntos de equilibrio del sistema.
y
x' = 0
x' > 0
x' < 0
x' > 0
x' < 0
x
Figura 31
352
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
La ceroclina y ' = g (x , y ) = 0 también divide el plano x-y en dos regiones, una en la
que y ' = g (x , y ) > 0 (donde “y” crece conforme transcurre el tiempo: en esta región
se dibujarán líneas de fuerza que apuntarán hacia arriba indicando el crecimiento de
“y” conforme se incrementa el tiempo) y la otra en la que y ' = g (x , y ) < 0 (donde “y”
decrece conforme transcurre el tiempo: hacia abajo indicando el decrecimiento de “y”
conforme se incrementa el tiempo). En la figura 32 se aprecia el movimiento en el
plano x-y descrito líneas arriba. En esta figura también puede notarse que en las dos
regiones en que la ceroclina y ' = g (x , y ) = 0 divide al plano x-y las líneas de fuerza
son opuestas (aquellas flechas que tienen el mismo color).
Es importante resaltar que las sendas de fase del sistema que se tracen en el plano x-y
cortarán verticalmente a la ceroclina x ' = f (x , y ) = 0 y horizontalmente a la ceroclina
y ' = g (x , y ) = 0.
y
y' > 0
y' < 0
y' < 0
x
Figura 32
Los puntos de equilibrio del sistema se encontrarán en las intersecciones de las
ceroclinas. Al superponer las ceroclinas de las figuras 31 y 32 en un mismo gráfico
(ver figura 33), el plano x-y quedará dividido en cuatro regiones en las que será
posible conocer la evolución temporal de “x” y de “y” a través de las líneas de
fuerza17. En la figura 33, sobre un punto de la ceroclina x ' = f (x , y ) = 0 se ha trazado
arbitrariamente el gradiente de f (x , y ) en dirección noroeste (de color azul). Dado que
la dirección del ∇f (x , y ) sobre cualquier punto de la ceroclina apunta hacia la
dirección donde f (x , y ) = x ' > 0, en consecuencia, las líneas de fuerza que se
encuentran arriba de la ceroclina x ' = f (x , y ) = 0 apuntarán hacia la derecha (este) ya
que en esa región, al ser x ' > 0, “x” aumentará conforme transcurra el tiempo. En la
región que se encuentra debajo de la ceroclina x ' = f (x , y ) = 0, al ser x ' < 0, las
líneas de fuerza apuntarán hacia la izquierda (oeste) indicando el decrecimiento de
“x” conforme aumenta el tiempo.
17
Es importante resaltar que las ceroclinas x ' = 0 y y ' = 0 se pueden interceptar en más de un punto de
equilibrio y, por tanto, dividir el plano x-y en más de cuatro regiones.
353
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
y
' x' = 0
y =0
∇f (x , y )
∇g (x , y )
E
x
Figura 33
y
y' = 0 x' = 0
x
Figura 34
354
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
En la figura 34 se han agregado algunas sendas de fase a la figura 33 en función de las
direcciones de las líneas de fuerza, las cuales sirven para prever el movimiento
dinámico del sistema a partir de cualquier punto inicial del plano x-y. Note que las
sendas de fase que cortan la ceroclina y ' = g (x , y ) = 0 tienen una pendiente nula en el
punto de corte (líneas punteadas que cortan a la ceroclina y ' = g (x , y ) = 0
horizontalmente), mientras que las sendas de fase que cortan la ceroclina
x ' = f (x , y ) = 0 tienen una pendiente infinita en el punto de corte (líneas punteadas
que cortan a la ceroclina x ' = f (x , y ) = 0 verticalmente).
v
v2
∇g ( x , y )
∇f ( x , y )
v
x
y' = 0 E v1
x' = 0
Figura 35: odo propio repulsor
18
Cuando las sendas de fase entran a un punto de equilibrio (sea un nodo o cualquier otro tipo de punto
crítico) se le suele denominar sumidero y cuando salen de él se le suele denominar fuente.
355
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
En las figuras 36-a y 36-b se muestran ejemplos de un nodo impropio.
f (x , y ) = 0
∇g (x , y ) g(x , y ) = 0 E
x
v v y' = 0
v2 v1
∇f ( x , y )
x' = 0
f (x , y ) = 0
y
y' = 0
v
v2
v
v1
E
x
∇g (x , y )
∇f ( x , y )
g (x , y ) = 0
x' = 0
Figura 36-b: odo impropio repulsor
356
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Un punto de silla es un punto de equilibrio que se caracteriza por ser estable en
algunas direcciones e inestable en otras. Específicamente, un punto de silla tiene un
par de ramas estables en las que las sendas de fase convergen hacia el punto de
equilibrio, y un par de ramas inestables en el que las sendas de fase se alejan del
punto de equilibrio. El resto de sendas de fase se mueven hacia el punto de silla
inicialmente pero luego se alejan de él. Debido a que la estabilidad sólo se observa en
el par de ramas estables, y ésta evidentemente no se obtiene como algo usual, en
general, al punto de silla se le suele clasificar como un punto de equilibrio inestable.
No obstante, en el caso muy poco probable, en el que en el instante inicial el estado
inicial del sistema se encontrase sobre alguna de las ramas estables, entonces el
sistema se movería a lo largo de esta rama hasta el punto de equilibrio. Es por esta
razón que algunos autores clasifican al punto de silla como un punto de equilibrio
condicionalmente estable19.
f (x , y ) = 0
∇f (x , y )
v
v1
x
E v
v2
x' = 0
∇g (x , y )
y' = 0
Figura 37: punto de silla
Los focos son puntos de equilibrio caracterizados por sendas de fase en forma de espiral
o que se aproximan cíclicamente a él (foco estable) o que se alejan cíclicamente de éste
(foco inestable). En la figura 38 se muestra un foco inestable.
19
Para un estudio de sistemas dinámicos condicionalmente estables recurra a Gandolfo, G. Economic
Dynamics. Study Edition. Springer-Verlag. Berlin. 1997.
357
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
y
y' = 0
∇f (x, y )
x
E
x' = 0
∇g ( x , y )
Figura 38: foco inestable
y' = 0
∇g (x , y )
E
x' = 0
∇f ( x , y ) x
358
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Puntos de equilibrio y criterios de estabilidad para sistemas lineales (de “n”
variables) de ecuaciones diferenciales ordinarias de primer orden con
coeficientes constantes
v v v
Hemos visto que para el sistema de la ecuación (LXXIV), X ' (t ) = A ⋅ X (t ) + b , si
A ≠ 0, entonces su único punto de equilibrio lo obtenemos de la ecuación
v
v v adjA ⋅ b
(LXXVIII), X* = A −1 ⋅ b = . Además, si A ≠ 0 y el sistema es
A
v v
( ) v v
homogéneo b = 0 , entonces el único punto de equilibrio sería el origen, X * = 0.
Teorema de Routh-Hurwitz:
Dado el sistema (LXXIV), las partes reales de todas las raíces (autovalores) del
polinomio característico de grado “n”:
a1 a3 a5 a7 a9 K 0
a 0 a2 a4 a6 a8 K 0
0 a1 a3 a5 a7 K 0
0 a0 a2 a4 a6 K 0
0 0 a1 a3 a5 K 0
0 0 a0 a2 a4 K 0
M M M M M O M
0 0 0 0 0 K a n
359
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
a1 a3 a5 a7 a9 K 0
a0 a2 a4 a6 a8 K 0
0 a1 a3 a5 a7 K 0
a1 a3 0 a0 a2 a4 a6 K 0
a 1 = a 1 > 0; > 0; KK ; >0
a0 a2 0 0 a1 a3 a5 K 0
0 0 a0 a2 a4 K 0
M M M M M O M
0 0 0 0 0 K an
a11 a12
a 11 = a 11 < 0 , > 0,
a 21 a 22
Criterio II: Una condición de estabilidad necesaria y suficiente para la matriz “A” es
que los menores principales dominantes de la matriz simétrica “B” alternen de signo,
empezando con signo negativo. Dónde:
1 1
a11 (a12 + a 21 ) K (a1n + a n1 )
2 2
1
(a 2 n + a n 2 )
1
B=
2
( 2
)
A + A = (a 21 + a12 )
1 T a 22 K
2
M M O M
1 1
(a n1 + a1n ) (a n 2 + a 2 n ) K a nn
2 2
20
Para una revisión más exhaustiva de estos criterios puede recurrir a Gandolfo, G. Economic Dynamics.
Study Edition. Springer-Verlag. Berlin. 1997.
360
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Criterio III: Sea a ij ≥ 0, i ≠ j, (es decir, todos los elementos que no pertenecen a la
diagonal principal de “A” son no negativos). Entonces, en este caso, las condiciones
de estabilidad necesarias y suficientes para la matriz “A” son las mismas que las del
criterio I. Note que estas condiciones implican que a ii < 0. Una matriz con a ii < 0 y
a ij ≥ 0, i ≠ j, es denominada “matriz Metzleriana”.
o
n
a ii < 0, a ii > ∑ a ji Dominación por columnas
j =1
j≠ i
dx
= x ' = a 11 x + a 12 y
dt x ' a a 12 x v' v
⇒ ' = 11 ⋅ ⇒ X (t ) = A ⋅ X (t ) (CXV)
dy = y ' = a x + a y y a 21 a 22 y
dt 21 22
v 0
Si A ≠ 0, entonces posee como único punto crítico a X * = , es decir, el origen.
0
El polinomio característico de (CXV) viene dado por:
a 11 − λ a 12
p (λ ) = = λ2 − (a 11 + a 22 )λ + (a 11 ⋅ a 22 − a 12 ⋅ a 21 ) = 0
a 21 a 22 − λ
p (λ ) = λ2 − (trA )λ + A = 0 (CXVI)
361
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
La estabilidad del sistema (CXV) quedará determinada por el hecho que el polinomio
característico (CXVI) tenga sus dos autovalores con parte real negativa.
}(− )
a1 a3 − trA 0
= > 0 ⇒ − trA ⋅ A > 0 ⇒ A > 0
a0 a2 1 A
trA = λ 1 + λ 2 y A = λ 1 ⋅ λ 2 (CXVIII)
Por lo que si los autovalores son reales y A > 0, entonces λ 1 y λ 2 deben tener el
mismo signo, mientras que si A < 0, los autovalores λ 1 y λ 2 deben tener signos
opuestos.
El estudio de los posibles estados del sistema (CXV) se realizará teniendo en cuenta
que:
362
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
A = λ1 ⋅ λ 2
∆=0 ∆ = (trA ) 2 − 4 A = 0
∆>0 ∆<0 IV VI V
I II
trA = λ 1 + λ 2
III
Figura 40
Debajo de la parábola tenemos que ∆ = (trA )2 − 4 A > 0, por lo que los autovalores
son reales y distintos; encima de la parábola se tiene que ∆ = (trA )2 − 4 A < 0, por lo
que los autovalores son complejo conjugados; mientras que a lo largo de la parábola
tenemos que ∆ = (trA )2 − 4 A = 0, por lo que los autovalores son reales e iguales.
363
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
de c 1 y c 2 . En este caso, en el largo plazo, las sendas de fase se aproximarán al
v
origen tangencialmente a la línea de acción del autovector V1 . En este caso se tendrá
un nodo impropio estable tal como el de la figura 36-a. De particular interés es el
v
hecho que si el punto inicial se encuentra justo sobre V1 , entonces c 2 = 0 y el
v
sistema se mueve a lo largo de la línea de acción del autovector V1 y se aproxima al
origen conforme transcurre el tiempo. De forma similar, si el punto inicial se
v
encuentra justo sobre V2 , entonces c1 = 0 y el sistema se mueve a lo largo de la
v
línea de acción de V2 , aproximándose al origen en el límite. Por tanto, el punto de
equilibrio es un nodo impropio asintóticamente estable. En segundo lugar, si
0 < λ 2 < λ1 , se verificará que cuando t → ∞ entonces e λ1 t >> e λ 2 t . Por lo que
v
cuando t → ∞ la solución tenderá a alinearse con c1e λ1 t V1 , esto es, en el largo plazo
v v
tendremos que X (t ) ≈ c1e λ1 t V1 . En este caso, las sendas de fase se alejarán del punto
v
de equilibrio tangencialmente a la línea de acción del autovector V1 conforme
transcurra el tiempo. Esto se debe a que tanto x (t ) e y (t ) crecen exponencialmente
con el tiempo.
t →∞
v
t →∞
( v
)
v
t →∞
[ (v v
lím X (t ) = lím c 1 e λ1t V1 + c 2 e λ 2 t V2 = lím e λt c 1 V1 + c 2 V2 )] v
→ 0.
364
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Caso III: autovalores reales distintos de signos opuestos
x ' = αx + βy x ' α β x
' ⇒ ' = (CXIX)
y = −βx + αy y − β α y
x = R cos θ
(CXX)
y = R senθ
21
Como ya se dijo, salvo que el estado inicial del sistema se encuentre sobre una de las dos ramas estables,
un punto de silla es inestable.
365
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Realizando operaciones elementales con las ecuaciones del sistema (CXX) se tiene
que:
R 2 = x 2 + y 2
y (CXXI)
tgθ =
x
Derivando respecto al tiempo las dos ecuaciones del sistema (CXXI) se tiene:
( )
RR ' = x (αx + β y ) + y (− βx + αy ) = α x 2 + y 2 = αR 2
1 ' x (− βx + αy ) − y (αx + βy )
2
⋅θ =
cos θ R 2 cos 2 θ
R ' = αR
R ' = αR
'
θ =
2
(
−βx + y 2 ⇒
)
θ ' = −β
(CXXIII)
R 2
R = ce αt
(CXXIV)
θ = θ 0 − βt
Donde “c” es una constante y θ(0 ) = θ 0 . Como β > 0, entonces “θ” decrece con el
tiempo, y por tanto el movimiento es en el sentido de las agujas del reloj22. Además, si
t → ∞ entonces:
R → 0 si α < 0
R → ∞ si α > 0
Por tanto, las sendas de fase son espirales que se aproximan o se alejan del origen
dependiendo del signo de “α”. En la figura 38 se muestra una espiral divergente.
22
En caso β < 0, el ángulo “θ” aumentaría con el transcurrir del tiempo, por lo que el sentido de giro sería
antihorario.
366
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Caso V: autovalores complejos con parte real nula
En la región VI, que es el eje “y” encima del origen, α = 0 y el punto de equilibrio
tiene un centro (vórtice) con una curva cerrada como trayectoria.
x ' = βy x ' 0 β x
' ⇒ ' = (CXXV)
y = −βx y − β 0 y
R ' = 0 R = c
' ⇒ (CXXVI)
θ = −β θ = θ 0 − βt
Donde “c” y θ(0 ) = θ 0 son constantes. Esto significa que las sendas de fase son
curvas cerradas (círculos o elipses) con centro en el origen. Si β > 0 el movimiento
será en el sentido de las agujas del reloj mientras que si β < 0 el movimiento será
antihorario. Un circuito completo alrededor del origen denota la fase del ciclo, que es
2 π β . En la figura 39 se muestra un centro.
367
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
A = λ1 ⋅ λ 2 Autovalores Punto de
D Estabilidad
trA = λ 1 + λ 2 de A equilibrio
Reales del
A >0 Nodo
∆>0 mismo signo: Estable
trA < 0 λ1 < λ 2 < 0 impropio
Reales del
A >0 Nodo
∆>0 mismo signo: Inestable
trA > 0 λ1 > λ 2 > 0
impropio
Tabla II
368
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Ejemplos:
Resuelva los siguientes sistemas dinámicos y realice su análisis cualitativo: represente
algunas de sus sendas de fase en el plano de fase x − y.
x ' = − y v x (0 ) 5 x ' 0 − 1 x v 5
1.- con X 0 = = ⇒ ' = ⋅ con X 0 = .
y ' = x y (0 ) 2
y 1 0 y 2
A = 1 > 0
⇒ ∆ = 0 2 − 4 (1) = −4 < 0 (255)
trA = 0
− λ −1 λ 1 = i α = 0
p (λ ) = = λ2 + 1 = 0 ⇒ ⇒ (256)
1 −λ λ
2 = − i β = 1
Para λ 1 = i :
−i −1 a 0
[A − λ 1 I ]vv 1 = ⋅ = ⇒ b = − ai
1 − i b 0
Si hacemos a = 1 ⇒ b = − i, entonces:
v 1
v1 =
− i
Para λ 2 = −i :
i −1 c 0
[A − λ 2 I ]vv 2 = ⋅ = ⇒ d = ci
1 i d 0
Si hacemos c = 1 ⇒ d = i , entonces:
v 1
v2 =
i
v v v v v
X (t ) = e 0 h 1 cos (t ) + h 2 sen (t ) = h 1 cos (t ) + h 2 sen (t ) (257)
369
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Con:
v v v 1 1 c + c 2
h 1 = c1 v1 + c 2 v 2 = c1 + c 2 = 1 (258)
−
i i − ic 1 + ic 2
v v v 1 1 i (c − c 2 )
h 2 = i (c 1 v 1 − c 2 v 2 ) = ic 1 − ic 2 = 1 (259)
− i i (c 1 + c 2 )
5 + 2i
c1 =
v 5 v c1 + c 2
X (0 ) = = h 1 = ⇒
2 (260)
2 − ic1 + ic 2 c = 5 − 2i
2 2
v 5 v −2
h1 = y h 2 = (261)
2
5
v 5 − 2 x (t ) = 5 cos (t ) − 2sen (t )
X (t ) = cos (t ) + sen (t ) ⇒ (262)
2
5 y (t ) = 2 cos (t ) + 5sen (t )
Para encontrar las órbitas del sistema en el plano de fase x − y vamos a utilizar
la ecuación (CXII):
dy y' x y2 x2
dx
=
x'
=
−y
∫
⇔ y ≠ 0 ⇒ − ydy = ∫ xdx + k1 ⇒ −
2
= k1 +
2
x2 y2
+ = k1 ⇒ x 2 + y 2 = 2 k1 ⇒ x 2 + y 2 = k (263)
2 2
En consecuencia, las órbitas del sistema vienen dadas por (263), que representa
una familia de circunferencias con centro en el origen (el punto de equilibrio es
un centro o vórtice) y radio k tal como se aprecia en la figura 41. En el
instante inicial, t = 0, se tiene que x 2 + y 2 = 29. El sentido del movimiento de la
senda de fase que satisface las condiciones iniciales se obtiene evaluando el
vector tangente a dicha senda (vector de velocidad de movimiento) en el instante
inicial. En el instante inicial, t = 0, tendríamos que el vector tangente es:
v v v
( )
F(x (0 ), y (0 )) = F0 = F(5, 2 ) = x ' (0 ), y ' (0 ) = (− 2,5 ). En consecuencia el movimiento
de la senda de fase sería contrario al sentido de giro de las agujas del reloj.
370
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
y
y' = 0
v
F0
v
2 X0
E
'
x =0 ∇g ( x , y ) 5
∇f ( x , y ) x
Figura 41
Se puede apreciar que en este caso las ceroclinas coinciden con los ejes
coordenados:
Por esta razón los gradientes de “f” y “g” se encuentran también sobre los ejes
coordenados:
0
∇ f ( x , y ) =
− 1
1
∇ g ( x , y ) =
0
En consecuencia, por debajo del eje “x” las líneas de fuerza horizontales
apuntarán hacia la derecha, y por encima del eje “x” apuntarán hacia la izquierda.
Asimismo, las líneas de fuerza verticales apuntarán hacia arriba a la derecha del
eje “y”, y hacia abajo a la izquierda del eje “y”. Note que todas las soluciones del
sistema dinámico, salvo el origen, son circunferencias. Asimismo, note que en las
intersecciones de las circunferencias con el eje “x”, y = 0, la dy dx no está
definida.
371
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
x ' = − x − y v x (0 ) 2 x ' − 1 − 1 x v 2
2.- con X 0 = = ⇒ ' = ⋅ con X 0 = .
'
y = x − y y (0 ) 1
y 1 − 1 y 1
A = 2 > 0
⇒ ∆ = (− 2 ) 2 − 4 (2 ) = −4 < 0 (264)
trA = −2
Para λ1 = −1 + (−1)i :
i −1 a 0
[A − λ1I]vv 1 = ⋅ = ⇒ b = ai
1 i b 0
Si hacemos a = 1 ⇒ b = i , entonces:
v 1
v1 =
i
Para λ 2 = −1 − (−1)i :
−i −1 c 0
[A − λ 2 I]vv 2 = ⋅ = ⇒ d = − ci
1 − i d 0
Si hacemos c = 1 ⇒ d = −i , entonces:
v 1
v2 =
− i
v v v v v
X (t ) = e − t h1 cos (− t ) + h 2 sen (− t ) = e − t h1 cos (t ) − h 2 sen (t ) (266)
23
En este punto es importante darse cuenta que, de acuerdo al sistema dado por (CXIX), β = −1 < 0.
372
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Con:
v v v 1 1 (c + c 2 )
h1 = c1v1 + c 2 v 2 = c1 + c 2 = 1 (267)
i
− i i (c1 − c 2 )
v v v 1 1 i (c − c 2 )
h 2 = i (c1v1 − c 2 v 2 ) = ic1 − ic 2 = 1 (268)
i − i − (c1 + c 2 )
v 2 v (c + c 2 ) c1 = (2 − i ) 2
X (0 ) = = h1 = 1 ⇒ (269)
1
i
1(c − c )
2 c 2 = (2 + i ) 2
v 2 v 1
h1 = y h 2 = (270)
1 − 2
Para encontrar las órbitas del sistema en el plano de fase x − y vamos a utilizar
la ecuación (CXII):
dy y' x−y 1 − (y x )
= = = ⇔ y ≠ −x (272)
dx x '
−x−y − 1 − (y / x )
dy 1− u u −1 u −1 u −1
= = ⇒ dy = dx ⇒ udx + xdu =
u + 1 dx
dx −1− u u +1 u + 1
u −1 u2 +1
udx − dx = − xdu ⇒ dx = − xdu ⇒ dx = u + 1 du
2
u +1 u +1 x u +1
dx u +1 1 x 2 + y2 y
∫ =
2 ∫ du ⇒ ln x + k = − ln
2
− tg −1 ; (x ≠ 0 )
x u +1 2 x x
373
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
1 x 2 + y2 y
F(x , y ) = ln x + k + ln + tg −1 = 0 ; (x ≠ 0 ) (274)
2
2 x x
1 x 2 + y2 y
F(x , y ) = ln x − 1, 2683 + ln + tg −1 = 0 ; (x ≠ 0 ) (275)
2
2 x x
La ecuación (275) nos representa una espiral convergente al origen que satisface
las condiciones iniciales del sistema dinámico. Para determinar el sentido del
movimiento de la senda de fase evaluaremos el vector tangente a la senda de fase
(vector de velocidad de movimiento) por ejemplo en el instante t = 0,
v v v
(
F(x (0 ), y (0 )) = F0 = F (2,1) = x ' (0 ), y ' (0 ) = (− 3,1).)
Por otro lado, del sistema (271) es fácil verificar que se cumple que:
x 2 + y 2 = 5e −2 t (276)
La ecuación (276) corresponde a una espiral, que puede interpretarse como una
circunferencia cuyo radio 5 e − t tiende a cero cuando t → ∞.
R 2 = x 2 + y 2 = 5e − 2 t ⇒ R = ± x 2 + y 2 = ± 5 e − t ⇒ c = ± 5 (277)
y (0 ) 1
tg (θ0 ) = = ⇒ θ0 = tg −1 (1 2 ) (278)
x (0 ) 2
θ( t ) = tg −1 (1 2 ) + t (279)
374
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
y
x' = 0 y' = 0
∇f (x , y ) ∇g (x , y )
v
F0
v
X0
2 x
Figura 42
Se puede apreciar que en este caso las ceroclinas vienen dadas por:
x ' = f (x , y ) = − x − y = 0 ⇒ y = − x
'
y = g (x , y ) = x − y = 0 ⇒ y = x
− 1
∇ f ( x , y ) =
− 1
( 1
∇ g x , y ) =
− 1
A = 1 > 0
⇒ ∆ = (2 ) 2 − 4 (1) = 0 (280)
trA = 2
375
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
El polinomio característico de la matriz “A” es:
1− λ 0
p(λ) = = λ2 − 2λ + 1 = 0 ⇒ {λ1 = λ2 = λ = 1 (281)
0 1− λ
Utilizando λ 1 = 1 :
0 0 a 0 a ∈ ℜ
[A − λ 1 I ]vv 1 = ⋅ = ⇒ 0⋅a +0⋅b = 0 ⇒
0 0 b 0 b ∈ ℜ
v a
v1 =
b
Utilizando λ 2 = 1 :
0 0 c 0 c ∈ ℜ
[A − λ 2 I ]vv 2 = ⋅ = ⇒ 0⋅c +0⋅d = 0 ⇒
0 0 d
0 d ∈ ℜ
v c
v2 =
d
v v v
Ya que [A − λI ]v = 0, con λ = 1, se satisface para cualquier autovector v ,
entonces deberemos escoger arbitrariamente cualquier par de autovectores
linealmente independientes para los autovalores λ 1 = λ 2 = λ = 1. Permítase que
dichos autovectores sean:
v 1 v 0
v1 = y v 2 =
0 1
v 1 0 x (t ) = c 1 e t
X (t ) = c1 e t + c 2 e t ⇒ (282)
0 1 y (t ) = c 2 e t
v 1 x (0 ) = c 1 e 0 = 1 ⇒ c 1 = 1 x ( t ) = e t
X (0 ) = ⇒ ⇒ (283)
1 y (0 ) = c 2 e 0 = 1 ⇒ c 2 = 1 y (t ) = e t
Para encontrar las órbitas del sistema en el plano de fase x − y vamos a utilizar
la ecuación (CXII):
dy y' y dy dx dy dx
dx
=
x '
=
x
⇔x≠0⇒
y
=
x
⇒ ∫ y
= ∫ x
+ k ⇒ ln y = ln x + k1
376
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
y y
ln = k1 ; (x ≠ 0 ∧ y ≠ 0 ) ⇒ = e k 1 = c > 0; (x ≠ 0 ∧ y ≠ 0 )
x x
y y
= c; (x ≠ 0 ∧ y ≠ 0) ⇒ = c; (x ≠ 0 ∧ y ≠ 0)
x x
y = cx; (x ≠ 0 ∧ y ≠ 0) (c ∈ ℜ) (284)
y = x ; (x ≠ 0 ∧ y ≠ 0 ) (285)
Se puede apreciar que en este caso las ceroclinas coinciden con los ejes
coordenados:
1
∇ f ( x , y ) =
0
0
∇ g ( x , y ) =
1
377
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
y
v
F0
1 v
∇g ( x , y ) X0
x
y' = 0 ∇f ( x , y ) 1
x' = 0
Figura 43
x ' = 2 y v x (0 ) − 2 x ' 0 2 x v − 2
4.- con X 0 = = ⇒ ' = ⋅ con X 0 = .
y ' = x y (0 ) 1 y 1 0 y 1
A = −2 < 0
⇒ ∆ = 0 2 − 4 (− 2 ) = 8 > 0 (286)
trA = 0
−λ 2 λ = 2
1
p (λ ) = = λ2 − 2 = 0 ⇒ (287)
1 −λ λ 2 = − 2
Para λ1 = 2 :
2 a 0 2
[A − λ1I ]vv1 = − 2
⋅ = ⇒ b = a
1 − 2 b 0 2
378
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Si hacemos a = 1 ⇒ b = 2 2 , entonces:
v 1
v1 =
2 2
Para λ 2 = − 2 :
2 2 c 0 2
[A − λ 2 I ]vv 2 = ⋅ = ⇒d = − c
1 2 d 0 2
Si hacemos c = 1 ⇒ d = − 2 2 , entonces:
v 1
v2 =
− 2 2
v 1 1 −
X (t ) = c1 e
2t
+ c2 e
2t
2 2 − 2 2
x (t ) = c1e + c 2e − 2 t
2t
(288)
y (t ) = 2 2 c1e 2 t − 2 2 c 2 e − 2t
− 2 +1
c1 + c 2 c1 =
v − 2
X (0 ) = = ⇒ 2
(289)
1 2 2
c1 − c2 − 2 −1
c 2 =
2 2
2
− 2 −1
x (t ) = − 2 + 1 e 2t
+
− 2t
e
2 2
(290)
− 2 +1 2 + 1 −
y (t ) =
e
2t
+ e
2t
2 2
379
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Para encontrar las órbitas del sistema en el plano de fase x − y vamos a utilizar
la ecuación (CXII):
dy y' x x2
dx
=
x'
=
2y
⇔ y ≠ 0 ⇒ 2 ydy = ∫ ∫ xdx + k ⇒ y 2 =
2
+k
y2 x2
− =1 (291)
k 2k
En consecuencia, las órbitas del sistema vienen dadas por (291), que representa
una familia de hipérbolas con vértices en el eje “y” para valores de k > 0, y con
vértices sobre el eje “x” para valores de k < 0. Independientemente del signo de
“k”, las ecuaciones de las asíntotas a dichas hipérbolas son y = ± 2 2 x . Note
que las pendientes de las rectas asíntotas coinciden con las tangentes de los
v v
ángulos de inclinación de los autovectores v1 y v 2 . Asimismo, note que en las
intersecciones de las sendas de fase con el eje horizontal, y = 0, la dy dx no está
definida. El punto de equilibrio es un punto de silla, tal como se aprecia en la figura
44. En el instante inicial, t = 0, se tiene que (1)2 1 − (− 2 )2 2 = k = −1. Por tanto,
la senda de fase que satisface las condiciones iniciales viene dada por:
x2 y2
− =1
2 1
Se puede apreciar que en este caso las ceroclinas coinciden con los ejes
coordenados:
x ' = f (x , y ) = 2 y = 0 ⇒ y = 0 : eje " x "
'
y = g (x , y ) = x = 0 ⇒ x = 0 : eje " y"
Por esta razón los gradientes de “f” y “g” se encuentran también sobre los ejes
coordenados:
0
∇ f ( x , y ) =
2
( 1
∇ g x , y ) =
0
En consecuencia, por encima del eje “x” las líneas de fuerza horizontales
apuntarán hacia la derecha, y por debajo del eje “x” apuntarán hacia la izquierda.
Asimismo, las líneas de fuerza verticales apuntarán hacia arriba a la derecha del
eje “y”, y hacia abajo a la izquierda del eje “y”. Note que todas las soluciones del
sistema dinámico, salvo el origen y las asíntotas, son hipérbolas.
380
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
y
∇f ( x , y ) v
v
1 v1
-2 F0
'
x =0 v x
v2
∇g ( x , y )
y' = 0
Figura 44
A = 2 > 0
⇒ ∆ = 0 2 − 4 (2 ) = −8 < 0 (292)
trA = 0
−λ 2 λ = 2 i α = 0
1
p (λ ) = = λ2 + 2 = 0 ⇒ ⇒ (293)
−1 − λ λ 2 = − 2 i β = 2
Para λ1 = 2 i :
a 0 2
[A − λ1I ]vv1 = − 2i 2
⋅ = ⇒ b = ai
− 1 − 2 i b 0 2
2
Si hacemos a = 1 ⇒ b = i , entonces:
2
1
v
v1 = 2
i
2
381
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Para λ 2 = − 2 i :
2i 2 c 0 2
[A − λ 2 I ]vv 2 = ⋅ = ⇒ d = − ci
− 1 2 i
d 0 2
2
Si hacemos c = 1 ⇒ d = − i, entonces:
2
1
v
v 2 = 2
− i
2
v v v
X (t ) = h1 cos 2 t + h 2sen 2 t (294)
Con:
v 1 1 (c1 + c 2 )
v v
h1 = c1v1 + c 2 v 2 = c1 2 + c 2
2 =
(295)
2
i − i i (c1 − c 2 )
2 2 2
v 1 1 i (c1 − c 2 )
v v
h 2 = i (c1v1 − c 2 v 2 ) = ic1 2 − ic 2 2 =
(296)
2
i − i − (c + c )
2 2 2
1 2
(c1 + c 2 ) c = 2 − 2 i 2
v 2 v 1
X (0 ) = = h1 = 2 ⇒ (297)
1 i (c − c ) c = 2 + 2 i 2
2
1 2
2
v 2 v 2
h1 = y h 2 = (298)
1 − 2
v 2 2
X (t ) = cos 2 t + sen 2 t
1 − 2
382
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
x (t ) = 2 cos 2 t + 2 sen 2 t
(299)
y (t ) = cos 2 t − 2 sen 2 t
Para encontrar las órbitas del sistema en el plano de fase x − y vamos a utilizar
la ecuación (CXII):
dy y' −x
dx
=
x '
=
2y
⇔ y ≠ 0 ⇒ 2 ydy = − xdx ⇒ 2 ydy = − xdx + k∫ ∫
x2 x2 y2
y2 = − +k⇒ + = 1; (k > 0 ) (300)
2 2k k
En consecuencia, las órbitas del sistema vienen dadas por (300), que representa
una familia de elipses con centro en el origen (el punto de equilibrio es un centro
o vórtice) tal como se aprecia en la figura 45.
x2 y2
+ =1
6 3
Se puede apreciar que en este caso las ceroclinas coinciden con los ejes
coordenados:
Por esta razón los gradientes de “f” y “g” se encuentran también sobre los ejes
coordenados:
0
∇ f ( x , y ) =
2
( − 1
∇ g x , y ) =
0
383
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
En consecuencia, por encima del eje “x” las líneas de fuerza horizontales
apuntarán hacia la derecha, y por debajo del eje “x” apuntarán hacia la izquierda.
Asimismo, las líneas de fuerza verticales apuntarán hacia arriba a la izquierda del
eje “y”, y hacia abajo a la derecha del eje “y”. Note que todas las soluciones del
sistema dinámico, salvo el origen, son elipses. Asimismo, note que en las
intersecciones de las elipses con el eje “x”, y = 0, la dy dx no está definida.
v
X0
1 v
∇f ( x , y ) F0
E
2 x
∇g ( x , y )
Figura 45
A = 9 > 0
⇒ ∆ = (− 6 )2 − 4 (9 ) = 0 (301)
trA = −6 < 0
−4 − λ −1
p (λ ) = = λ2 + 6 λ + 9 = 0 ⇒ {λ1 = λ 2 = λ = −3 (302)
1 −2−λ
Para λ = −3 :
−1 −1 a 0
[A − λI ]vv1 = ⋅ = ⇒ a = −b
1 1 b 0
v 1
v1 =
− 1
384
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Como apreciamos sólo existe un autovector asociado a λ = −3. Para encontrar un
v
segundo autovector que sea linealmente independiente de v 1 vamos a utilizar la
expresión (CVI):
−1 −1 c 1
[A − λI ]vv 2 v
= v1 ⇒ ⋅ = ⇒ c + d = −1 (303)
1 1 d − 1
v 0
v2 =
− 1
x (t ) −3t 1 −3t 1 0
= c1e + c 2 e t +
y (t ) − 1 − 1 − 1
x (t ) = c1e −3 t + c 2 te −3 t
(304)
y (t ) = − c1e − 3 t − c 2 (t + 1)e − 3 t
x (0 ) = c1 = 2
(305)
y (0 ) = − c1 − c 2 = 2 ⇒ c 2 = −4
x (t ) = 2 e −3 t − 4 te −3 t
(306)
y (t ) = −2 e − 3 t + 4 (t + 1)e − 3 t
Para encontrar las órbitas del sistema en el plano de fase x − y vamos a utilizar
la ecuación (CXII):
dy y' x − 2y 2(y x ) − 1
= = = ⇔ y ≠ −4 x (307)
dx x '
− 4x − y (y / x ) + 4
Si hacemos el siguiente cambio de variable:
dy 2u − 1 2u − 1 2u − 1
= ⇒ dy = dx ⇒ udx + xdu =
u + 4 dx
dx u+4 u+4
385
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
2u − 1 (u + 1)2 dx u+4
udx − dx = − xdu ⇒
dx = − xdu ⇒ = − du
u+4 u+4
x (u + 1)2
dx u+4 3
∫ =− ∫ du ⇒ ln x + k = − ln u + 1 + ; (u ≠ −1)
x (u + 1)2 u +1
y 3 y
ln x + k = − ln +1 + ; ≠ −1
x y x +1 x
x+y 3x
ln x + k = − ln + ; (y ≠ − x )
x x+y
3x
F(x , y ) = ln x + y − + k = 0; (y ≠ − x ) (309)
x+y
3x
F(x , y ) = ln x + y − + 0,1137 = 0; (y ≠ − x ) (310)
x+y
La ecuación (310) nos representa una senda de fase convergente al origen que
satisface las condiciones iniciales del sistema dinámico. Para determinar el
sentido del movimiento de la senda de fase evaluaremos el vector tangente a la
senda de fase (vector de velocidad de movimiento) por ejemplo en el instante
v v v
( )
t = 0, F(x (0 ), y (0 )) = F0 = F (2, 2 ) = x ' (0 ), y ' (0 ) = (− 10, −2 ).
Se puede apreciar que en este caso las ceroclinas vienen dadas por:
x ' = f (x , y ) = −4 x − y = 0 ⇒ y = −4 x
'
y = g (x , y ) = x − 2 y = 0 ⇒ y = x 2
− 4
∇ f ( x , y ) =
− 1
1
∇ g ( x , y ) =
− 2
386
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
En consecuencia, a la izquierda de la ceroclina x ' = 0 las líneas de fuerza
horizontales apuntarán hacia la derecha, y a la derecha de ella apuntarán hacia la
izquierda. Asimismo, las líneas de fuerza verticales apuntarán hacia arriba por
debajo de la ceroclina y ' = 0, y hacia abajo por encima de ella. Note que en las
intersecciones de las sendas de fase con la isoclina x ' = 0 (y = −4 x ), la dy dx
no está definida. El punto de equilibrio es un nodo impropio estable.
x' = 0
y
∇f ( x , y )
v y' = 0
2 X0
v
F0
v 2 x
v1
v
v2
∇g (x , y )
Figura 46
Realice el análisis cualitativo del modelo Keynesiano IS-LM (Tu 1994), dado por el
sistema (215), para el caso III.
trA = −1
A = 1
γ = 0,5; s = 0,5; µ = 0,75; θ = 1
⇒ ∆ = −3
I0 = T = G = M = 1 *
Y = 1,5
*
r = 1
387
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Teniendo en cuenta la ecuación (219) y los valores de equilibrio, se tiene:
0,5u − 1 0,5u − 1
d1 ⋅ du + u ⋅ d(d1) = ⋅ d(d1) ⇒ − u ⋅ d(d1) = d1 ⋅ du
0,5 + 0,75u 0,5 + 0,75u
2 3 3
ln d1 − ln + tg −1 u + k = 0
3 2
4 + 3u 2
2 3 3 d
ln d1 − ln + tg −1 2 + k = 0
2 3 2 d1
d2
4 + 3
d
1
2 3 3 r − 1
ln Y − 1,5 − ln + tg −1 + k = 0
2 3 2 Y − 1,5
r −1
4 + 3
Y − 1,5
Donde Y ≠ 1,5.
2 3 3 r − 1
ln Y − 1,5 − ln + tg −1 + 0,8290 = 0
2 3 2 Y − 1,5
r −1
4 + 3
Y − 1,5
La ecuación anterior nos representa una espiral convergente al punto de equilibrio que
satisface las condiciones iniciales del sistema dinámico. Para determinar el sentido del
movimiento de la senda de fase evaluaremos el vector tangente a la senda de fase
(vector de velocidad de movimiento) por ejemplo en el instante t = 0,
v v v
( )
F(Y(0), r(0)) = F0 = F(2, 0,8) = Y '(0), r '(0) = (− 0,1, 0,6).
' 2
Y = f (Y, r) = −0,5Y − 0,75r + 1,5 = 0 ⇒ r = 2 − Y
3
'
r = g(Y, r) = Y − 0,5r − 1 = 0 ⇒ r = 2Y − 2
− 0,5
∇f (Y, r) =
− 0,75
1
∇g(Y, r) =
− 0,5
389
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
En consecuencia, a la izquierda de la ceroclina Y ' = 0 las líneas de fuerza
horizontales apuntarán hacia la derecha, y a la derecha de ella apuntarán hacia la
izquierda. Asimismo, las líneas de fuerza verticales apuntarán hacia arriba por debajo
de la ceroclina r ' = 0, y hacia abajo por encima de ella. Note que en las
intersecciones de las sendas de fase con la isoclina Y ' = 0 (r = 2 − 2 3 Y), la dr dY no
está definida.
r d2
r' = 0
∇g(Y, r)
∇f (Y, r)
v
F0
v d1
X0
r* Y' = 0
Y*
Figura 47
Análisis cualitativo de la estabilidad para sistemas autónomos de dos ecuaciones
diferenciales ordinarias de primer orden no lineales
En este apartado estudiaremos la estabilidad local de un sistema autónomo de dos
ecuaciones diferenciales ordinarias de primer orden no lineales, tal como el dado por
la ecuación (CIX), a través del procedimiento de linealización, que consiste en
expandir las funciones “f” y “g” en series de Taylor alrededor de algún punto de
equilibrio del sistema, sin considerar los términos de orden superior (aquellos
términos de la serie de Taylor que poseen derivadas parciales de orden mayor a uno).
Luego se analizan las propiedades de estabilidad del sistema linealizado para
finalmente ser extendidas al sistema no lineal original.
v
(
Para examinar si un punto de equilibrio X * = x * , y * del sistema (CIX):)
x ' f ( x , y )
' =
y g (x , y )
24
Ver apéndice.
390
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
(
x ' = f (x , y ) ≈ f x * , y* + ) (
∂f x * , y*
x − x* +
)(∂f x * , y*
)y − y*
( )( )
∂x ∂y
'
= ( ) ≈ (* *
+ ) (
∂g x * , y*
− * )(
+
∂g x * , y*
)y − y*
( )( )
y g x , y g x , y x x
∂x ∂y
Pero ya que:
(
x ' f x * , y*
=
) = 0
(
y ' g x * , y*
) 0
Entonces:
x ' = f (x , y ) ≈
(
∂f x * , y* )(
x − x* + )
∂f x * , y* (
y − y*
)( )
∂x ∂y
' (
∂g x * , y* )( )
∂g x * , y* ( )( )
y = g (x , y ) ≈ x − x* + y − y*
∂x ∂y
∂f
(x , y )
∂f x * , y *
* *
( )
x f ( x , y ) ( ) )((xy −− xy ))
) (
'
∂x ∂y x − x* * *
*
' = ≈
y g (x , y ) ∂g
* *
x ,y ( * *
)
∂g x , y y − y* (
= J x ,y
) (
*
(CXXVII)
∂x ∂y
Si hacemos:
v z x − x*
Z = 1 = *
z 2 y − y
v
( v
Z ' = J x * , y* Z ) (CXXVIII)
Tabla III
392
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Teorema de Lyapunov (1955): Sean f (x , y ) y g (x , y ) funciones de clase uno, C1, y sea
v
( )
X * = x * , y * algún punto de equilibrio tanto del sistema (CIX) así como también del
v
sistema (CXXVIII). Para el sistema lineal (CXXVIII), X * = x * , y * será ( )
asintóticamente globalmente estable si y solo si los dos autovalores de J (x *
, y* )
(
tienen partes reales negativas, o equivalentemente, si y solo si J x , y tiene traza * *
)
negativa y determinante positivo. Dado que (CXXVIII) se “comporta”
v
( )
aproximadamente como (CIX) cerca de X * = x * , y * , entonces deducimos que en
v
( )
este caso X * = x * , y * es un punto de equilibrio asintóticamente localmente estable
del sistema (CIX).
( ) ( )
f x x * , y * ⋅ g y x * , y * ≠ 0 en todo ℜ 2 o f y * * x(x , y ) ⋅ g (x , y ) ≠ 0
* *
en todo ℜ 2 .
v*
( )
Entonces X = x * , y * es un punto de equilibrio asintóticamente globalmente estable.
Aplicaciones económicas
f (K , C ) = K ' = AK − BK 2 − C
(CXXVIII)
g (K , C ) = C ' = D (A − 2 BK )C
Solución:
Para determinar los puntos de equilibrio de este sistema vamos a calcular las
ceroclinas, esto es:
f (K , C ) = K ' = AK − BK 2 − C = 0 ⇒ C = AK − BK 2
(317)
g (K , C ) = C ' = D (A − 2 BK )C = 0 ⇒ C = 0 ∨ K = A 2 B
393
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Se aprecia que en el plano K-C la ceroclina K ' = 0 es una parábola
( )
C = AK − BK 2 , y que C = 0 (eje “K”) y K = A 2 B (recta perpendicular al eje
“K”) son parte de la ceroclina C ' = 0. El gráfico de estas ceroclinas puede
apreciarse en la figura 48. En la intersección de las ceroclinas se encontrarán los
puntos
v
(
X * = K * , C* ) de equilibrio del sistema (CXXVIII):
v v v
(
X1* = (0, 0 ), X *2 = (A B , 0 ) y X *3 = A 2 B , A 2 4 B . )
C
(III) (IV )
v
X *3
A 2 4B
(I ) (II )
K' = 0
C' = 0
v v
X1* A 2B X *2 K
Figura 48
Tabla IV
394
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
C
∇g (x , y )
v
X *3
∇f ( x , y )
∇ g ( x , y ) ∇f ( x , y )
v v
X1* A 2B X *2 K
∇g (x , y )
Figura 49
Para posteriormente determinar los gradientes de “f” y “g” sobre cualquier punto
de las ceroclinas f (K , C ) = K ' = 0 y g (K , C ) = C ' = 0, primero vamos a calcular
los gradientes de “f” y “g” para cualquier punto (K , C ) del plano de fase K-C:
f K (K , C ) A − 2 BK
∇ f (K , C ) = =
f C (K , C ) − 1
g K (K , C ) − 2 BDC
∇ g (K , C ) = =
g C (K , C ) D (A − 2 BK )
Para K ≥ A 2 B :
f K (K , C ) A − 2 BK ≤ 0
∇ f (K , C ) = =
f C (K , C ) − 1 < 0
Para K < A 2 B :
f K (K , C ) A − 2 BK > 0
∇ f (K , C ) = =
f C (K , C ) − 1 < 0
395
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Para K > A 2 B ∧ C = 0 :
g K (K , C ) 0 = 0
∇ g (K , C ) = =
g C (K , C ) D (A − 2 BK ) < 0
Para K < A 2 B ∧ C = 0 :
g K (K , C ) 0 = 0
∇ g (K , C ) = =
g C ( K , C ) D ( A − 2 BK ) > 0
Para K = A 2 B ∧ C > 0 :
g K (K , C ) − 2 BDC < 0
∇ g (K , C ) = =
g C (K , C ) 0 = 0
En la figura 49 se han dibujado los gradientes de “f” y “g” sobre las ceroclinas
f (K , C ) = K ' = 0 y g (K , C ) = C ' = 0 respectivamente. El gradiente de “f” apunta
hacia las región donde f (K , C ) = K ' > 0 y el gradiente de “g” apunta hacia la
región donde g (K , C ) = C ' > 0 . Esto corrobora los resultados obtenidos en la tabla
IV.
Ahora para realizar el análisis de estabilidad vamos a calcular el jacobiano del
sistema de ecuaciones diferenciales (CXXVIII) y lo evaluaremos en sus puntos
de equilibrio:
f K (K , C ) f C (K , C ) A − 2 BK −1
J (K , C ) = =
g
K ( K , C ) g C ( K , C ) − 2 BDC D ( A − 2 BK )
v
Para X1* = (0, 0 ), se tiene:
A −1
J (0,0 ) =
0 DA
( )
trJ K * , C * = trJ (0, 0 ) = A + DA = A (1 + D ) > 0
( )
J K * , C * = J (0, 0 ) = DA 2 > 0
( )
∆ = [A (1 + D )]2 − 4 DA 2 = [A (1 − D )]2 ≥ 0
v
Por tanto, X1* = (0, 0 ) es asintóticamente localmente inestable. Si D = 1 ⇒ ∆ = 0,
v
entonces X1* = (0, 0 ) sería un nodo propio asintóticamente localmente inestable,
v
pero si D ≠ 1 ⇒ ∆ > 0, entonces X1* = (0, 0 ) sería un nodo impropio
asintóticamente localmente inestable.
396
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
v
Para X *2 = (A B , 0 ), se tiene:
−A −1
J (A B , 0 ) =
0 − DA
( )
trJ K * , C * = trJ (A B ,0 ) = − A − DA = − A (1 + D ) < 0
J (K *
)
, C * = J (A B , 0 ) = DA 2 > 0
( )
∆ = [− A (1 + D )]2 − 4 DA 2 = [A (1 − D )]2 ≥ 0
v
Por tanto, X *2 = (A B , 0 ) es asintóticamente localmente estable. Si
v
D = 1 ⇒ ∆ = 0, entonces X *2 = (A B , 0 ) sería un nodo propio asintóticamente
v
localmente estable, pero si D ≠ 1 ⇒ ∆ > 0, entonces X *2 = (A B , 0 ) sería un
nodo impropio asintóticamente localmente estable.
v
( )
Para X *3 = A 2 B , A 2 4 B , se tiene:
J (A 2 B , A )
2 0 −1
4B = 2
− DA 2 0
( ) (
trJ K * , C * = trJ A 2 B , A 2 4 B = 0 )
( ) ( )
J K * , C * = J A 2 B , A 2 4 B = − DA 2 2 < 0
( )
∆ = (0 )2 − 4 − DA 2 2 = 2 DA 2 > 0
v
( )
Por tanto, X*3 = A 2B , A2 4B es asintóticamente localmente inestable. En este
caso, al ser el determinante del jacobiano negativo y el discriminante mayor a
v
( )
cero, resulta que X*3 = A 2B , A2 4B es un punto de silla normalmente
asintóticamente localmente inestable.
Se pide:
Solución:
a) Para resolver este apartado primero vamos a derivar respecto del tiempo a“k”
y a “kh” respectivamente:
k h = K H AL ⇒ k 'h = H
(
K ' (AL) − K H A 'L + AL' )
(AL)2
' (
K '(AL) − K A 'L + AL' )
k = K AL ⇒ k =
(AL)2
K'
k 'h = H − k h (η + γ)
AL
(319)
' K'
k = − k(η + γ)
AL
k 'h = λy − k h (η + γ)
' (321)
k = δy − k(η + γ)
y = k β k θh (322)
f (k h , k) = k 'h = λk βk θh − k h (η + γ)
(323)
g(k h , k) = k ' = δk βk θh − k(η + γ)
b) Para determinar los puntos de equilibrio de este sistema vamos a calcular las
ceroclinas, esto es:
647 a 48
1 1− θ 1− θ
η + γ β β
k h = 0 ⇒ k =
' ⋅ kh
= ak h β > 0 ya que : a > 0 ∧ k h > 0
λ
(324)
6474 b 48 4
1 θ θ
' δ 1−β − β
k = 0 ⇒ k = 1
⋅ kh = bk h β > 0 ya que : b > 0 ∧ k h > 0
1−
η + γ
1 1
1−β β
1−β − θ λθδ1− θ 1−β − θ
v*
( )
X = k *h , k* =
λ δ
η+γ
,
η+γ
(325)
v
Reemplazando X* en (322) resulta:
1
λθδβ 1−β −θ
y* = (326)
(η + γ)β+ θ
399
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
(v )
Se aprecia que los valores de equilibrio X*, y* dependen directamente de las
tasas de ahorro “λ” y “δ”. Por tanto, un incremento en cualquiera de estas
tasas significaría tener en el largo plazo un mayor stock de capital humano
por unidad de trabajo calificado, un mayor stock de capital físico por unidad
de trabajo calificado y un mayor ingreso por unidad de trabajo calificado.
Asimismo, ante una mayor tasa de progreso tecnológico “η” o ante una
mayor tasa de crecimiento poblacional “γ”, en el largo plazo, menores serán
los valores de equilibrio por unidad de trabajo calificado.
c) Ahora vamos a dibujar las ceroclinas del sistema (324), que se interceptan en
el punto de equilibrio del sistema:
1 1 β(1−β) θβ
1−β β
1−β − θ λθδ1− θ 1−β − θ
v*
( )
X = k *h , k* =
λ δ
,
η+γ
a a
= θ +β −1 , b θ+β −1
b b
η + γ
1 1− θ 1− θ
η+ γ β
Para la ceroclina, k 'h = 0, k = ⋅ kh β = ak h β se tiene:
λ
1− θ −β
dk (1 − θ) β
=a kh >0 ∀ kh > 0
dk h β
1− θ − 2β
2
d k = a(1 − θ)(1 − θ − β) k β > 0 ∀ kh > 0
dk 2 β2
h
h
1 θ θ
δ 1−β
Para la ceroclina, g(k h , k) = k ' = 0, k = ⋅ k h1−β = bk h1−β se tiene:
η+ γ
dk θb 1
= ⋅ 1−θ −β > 0 ∀ k h > 0
dk h (1 − β)
k h 1−β
2 θ+ 2(β −1)
d k bθ(θ + β − 1)
2 = k h 1−β < 0 ∀ k h > 0
dk h (1 − β)2
400
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
En consecuencia, la ceroclina k ' = 0 es una función estrictamente creciente
y estrictamente cóncava en k h > 0.
1− θ −β
dk (
1 − θ) 1−β
lím = lím a kh →0
k h →0 dk β
h f (k h , k)= k 'h =0 k h → 0
dk
θb 1
lím = lím ⋅ 1− θ −β → +∞
k h →0 dk
g(k h , k)= k ' = 0 k h →0 (1 − β)
h
k h 1−β
k' = 0
k* v
X*
∇g (k h , k )
∇f (k h , k )
k *h kh
Figura 50
g k (k h , k) δθk βk θh −1
∇g(k h , k) = h =
g k (k h , k) δβk k h − (η + γ)
β −1 θ
401
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
Se puede apreciar que para cualquier punto de la ceroclina f (k h , k) = k 'h = 0,
con k h > 0 y k > 0, se tiene que:
(η7 )
+ γ4
f k h (k h , k) 64 8
β θ −1 λk βk θh −1(θ − 1) < 0
∇f (k h , k) = β θ −1
= λθk k h − λk k h = β −1 θ
f k (k h , k) β −1 θ λβ k k h > 0
λβ k k h
g k h (k h , k)
δθk βk θh −1 β θ −1 > 0
∇g(k h , k) = δθk k
= δβk k h − δk k h = β−1 θ h
β −1 θ β −1 θ
g k (k h , k) 1424
(
3
) δk k h (β − 1) < 0
η + γ
v
Para X* se tiene:
λβ(η + γ)
(θ − 1)(η + γ)
(
J k *h , k* ) =
δθ(η + γ)
δ
(β − 1)(η + γ)
λ
( )
trJ k*h , k * = (η + γ)(θ + β − 2) < 0 en todo ℜ2+
( )
J k *h , k* = (η + γ)2[1 − (θ + β)] > 0 en todo ℜ2+
402
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
v
Por tanto, gracias al teorema de Lyapunov, X* es un nodo impropio
asintóticamente localmente estable del sistema (323). Además, gracias a que
(−8
) 67 (+)8 67 (−8
) 67 (+)8
v
( ) ( )
v 67
se verifica que f k h X* ⋅ g k X* = (θ − 1)(η + γ) ⋅ (β − 1)(η + γ) > 0 en todo ℜ+2 ,
v
por el teorema de Olech, X* también es un nodo impropio asintóticamente
globalmente estable del sistema (323).
(
K ' = K sK α −1 − γ
'
) (327)
P = K δ − θP
Donde 0 < s < 1, 0 < α < 1, γ > 0, θ > 0, y δ > 1. Encontrar el punto de equilibrio
(K , P ) en el primer cuadrante abierto, y verificar (si es posible) la estabilidad
* *
Para resolver el problema, primero vamos a determinar las ceroclinas. Para ello
vamos a igualar a cero las dos ecuaciones que aparecen en (327), de donde
resulta:
1
( ) s
f (K, P) = K ' = K sK α −1 − γ = 0 ⇒ K = 0 ∨ K =
γ
1− α
(328)
1 1
g(K, P) = P = K − θP = 0 ⇒ K = θ P
' δ δ δ
(s γ)δ (1−α)
(K , P ) = (s γ) (
* * 1 1− α)
,
θ
403
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
De donde se obtiene:
( )
trJ K*, P* = (1
α2 1) γ{ − θ{ < 0
−3
(−) (+) (+)
( )
J K*, P* = −(1
α2 1) γ{ θ{ > 0
−3
(−) (+)(+)
K'
K ' = sK α − γK = 0 ⇒ K ' + γK = sK α ⇒ + γK1−α = s (329)
Kα
y = K1−α (330)
dy dy dK K' y'
y' = = ⋅ = (1 − α)K −α K ' ⇒ = (331)
dt dK dt Kα (1 − α)
y'
+ γy = s ⇒ y' + (1 − α)γy = (1 − α)s (332)
(1 − α)
404
MATEMÁTICAS PARA EL AÁLISIS ECOÓMICO
Por (XXI), la solución general de (332) viene dada por:
[
y(t) = e − γ(1− α)t A + e γ(1−α)t (1 − α)sdt
∫ ]
s s
y(t) = e − γ(1−α)t A + e γ(1−α)t = Ae− γ(1− α)t + (333)
γ γ
1
s 1− α
K(t) = Ae− γ(1−α)t + (334)
γ
1
s 1−α s
K 0 = A + ⇒ A = K10−α − (335)
γ γ
1
s s 1−α
K(t) = K10−α − e − γ(1− α)t + (336)
γ γ
1
s 1−α
lím K(t) = = K*
t → +∞ γ
405
CIRO BAZÁ ECUACIOES DIFERECIALES ORDIARIAS
406