EDO Unidad 2 ING

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 41

ECUACIONES

DIFERENCIALES
RUBÉN ZÁRATE ROJAS
Ingeniero Civil - Ingeniero Industrial
Magíster en Ingeniería Civil ©
Magíster en Ingeniería Industrial

E-mail: [email protected]
Ecuaciones diferenciales de
Primer Orden y Primer Grado

P ( x,y ) dx + Q ( x,y ) dy =
0

2
Variables Separadas
P ( x,y ) dx + Q ( x,y ) dy =
0
P ( x,y ) = f1 ( x )
Q ( x,y ) = f2 ( y )
f1 ( x ) dx + f2 ( y ) dy = 0

∫ f1 ( x ) dx + ∫ f2 ( y ) dy =
0
F1 ( x ) + F2 ( y ) =
C
3
Variables Separables
P ( x,y ) dx + Q ( x,y ) dy =
0
P ( x,y ) = f1 ( x ) g1 ( y )
Q ( x,y ) = f2 ( x ) g2 ( y )
f1 ( x ) g1 ( y ) dx + f2 ( x ) g2 ( y ) dy =
0

4
Variables Separables (II)
f1 ( x ) g1 ( y ) f2 ( x ) g2 ( y )
dx + dy =
0
f2 ( x ) g1 ( y ) f2 ( x ) g1 ( y )
f1 ( x ) g2 ( y )
dx + dy =
0
f2 ( x ) g1 ( y )

F(x) + G(y) =
C
5
Homogéneas
P ( x,y ) dx + Q ( x,y ) dy = 0
f ( λx, λy ) =λ f ( x,y )
n

 y n  y
P ( x,y ) x=
= n
P  1,  x f  
 x x
 y y
Q ( x,y ) x=
= n n
Q  1,  x g  
 x x
n  y y
x f   dx + x g   dy =
n
0
x x 6
Homogéneas (II)
y y
f   dx + g   dy = 0
x x
y dy dz
z = ⇒ y = zx = x +z
x dx dx
f ( z ) dx + g ( z ) dy =
0 dy
= −
f(z)
dx g(z)

7
Homogéneas (III)
f(z) dz
− =x +z
g(z) dx
dx dz
=
x ϕ (z) − z
E.D. a variables separadas

8
Reducibles a Homogéneas

(a1x + b1y + c1 ) dx + (a2x + b2y + c 2 ) dy =


0
a1b2 − a2b1 ≠ 0
x= u+α ; y = v+β
a1α + b1β +=
c1 0 ; a2α + b2 β + =
c2 0
(a1u + b1v ) du + (a2u + b2v ) dv =
0
9
Reducibles a Homogéneas

(a1x + b1y + c1 ) dx + (a2x + b2y + c 2 ) dy =


0
a1b2 − a2b1 =
0
a1 b1
= = k
a2 b2
=
z a1x + b1y
E.D. a variables separables
10
Exactas
P ( x,y ) dx + Q ( x,y ) dy =
0
Py ( x,y ) = Qx ( x,y )
x y

∫ P ( x,y ) dx + ∫ 0
Q ( x ,y ) dy =
0
x0 y0

F ( x,y ) = C
11
Factor Integrante
P ( x,y ) dx + Q ( x,y ) dy = 0
Px ( x,y ) ≠ Qx ( x,y )
Py − Qx ∫
= A ( x ) ⇒= = f (x)
A( x )dx
z e
Q
Py − Qx ∫
= B (y) ⇒= = g( y)
B( y )dy
− z e
P
z P ( x,y ) dx + z Q ( x,y ) dy =
0
E.D. Exacta 12
Lineal
b(x)
y '+ a ( x ) y =
z ( x ) y '+ z ( x ) a ( x ) y =
z (x) b(x)
z '(x) = z (x) a(x)
z ( x ) y '+ z ' ( x ) y =
z (x) b(x)

(z (x) y)' = z (x) b(x)


13
Lineal
z ( x ) y = ∫ z ( x ) b ( x ) dx

z (x) = e ∫ a( x ) dx

z '(x) = a(x) e ∫ a( x ) dx

z '(x) = z (x) a(x)

14
Lineal

∫ a( x ) dx ∫ a( x ) dx
b ( x ) dx + C
e y ∫ e

∫− a( x )dx  b(x) e ∫ a( x )dx dx + C 


 ∫
y e


15
Bernoulli
y '+ a ( x ) y =
b(x) y n

n ≠ 0,n ≠ 1
k −1
y = z y' = k z z'
k

kz z '+ a ( x ) z =
k −1 k
b(x) z
kn

a(x) b ( x ) kn−k +1
z '+ z= z
k k
16
Bernoulli
kn − k + 1 =0
1
k=
1− n
z'+ (1 − n ) a ( x ) z = (1 − n ) b ( x )
E.D. Lineal

− (1−n) a dx (1 − n) b e(1−n) ∫ a dx dx + C 

1−n
y =
e
 
17
Riccati
y ' = a ( x ) y +b ( x ) y + c ( x )
2

y = yo ( x ) Solución particular
=
y yo + z Solución general

y=' y o '+ z '


z=' ( 2a yo + b ) z + a z 2

E. D. de Bernoulli

18
Clairaut
y x y '+ f ( y ' )
=
y' = p
=y c x + f ( c ) I. G.
=y p x + f (p )

 df ⇒( F )x,y =
0 I. S.

 x + =0
 dp

19
Lagrange
= y x f ( y ') + g ( y ')
p y x f (p ) + g (p )
y ' ==
 df dg  dp
p =f ( p ) +  x + 
 dp dp  dx
dx
= a ( p ) x + b ( p ) E.D. lineal
dp
20
Lagrange (II)
x = F ( p,c )
y x f (p ) + g (p )
=
 (⇒ G x,y,c =
0)
 x = F ( p,c )
I. G.

21
22
Lineal
Particular

− a( x )dx ∫

a( x )dx
yp = e b(x) e dx
Representa el nivel de equilibrio de y(t)

23
Lineal
Complementaria

yc = C e ∫− a( x )dx
dx
Representa la desviación del equilibrio

Si lim y c = 0
x →∞

y(t) es dinámicamente estable

24
1. Modelo de ajuste de precios EVANS

Establece que el precio P de un cierto


artículo varía de forma que su razón de
cambio respecto al tiempo es proporcional
a la demanda en exceso (Qd– Qs)

25
1. Modelo de ajuste de precios EVANS
dP
= m ( Qd − Qs )
dt
(Qd – Qs): exceso de la demanda

Qd= c + d P : función lineal de la demanda

Qs= g + h P: función lineal de la oferta

m: constante positiva

26
1. Modelo de ajuste de precios EVANS
En el equilibrio:
Qd = Q s ; c + b P = g + h P
c−g
El precio de equilibrio es: P=
h−b

= m ( ( c + bP ) − ( g + hP ) )
dP
dt
dP
= m ( b − h ) P + m ( c − g)
dt

27
1. Modelo de ajuste de precios EVANS

dP
= AP + B
dt
ED lineal: A = m (b - h) ; B = m (c - g)
∫ A dt  − ∫ A dt 
P(t) e  ∫ Be + C
 
∫ − A dt B − At
∫ A dt = A t ∫ Be ∫ Be
− At
= dt = − e
A
 B −A t  B
P(t) =e  − e + C =
At
− + Ce At

 A  A
28
1. Modelo de ajuste de precios EVANS
B B
P(0) =− +C⇒C =P(0) +
A A
 B  At B
P(t) =P(0) +  e −
 A A
En condiciones normales: h > b

 c − g  −m(h−b)t c − g
P(t) =
P(0) −  e +
 h −b h−b


P(t) =P(0) − P  e − m(h −b)t
+ P
 
29
2. El modelo de DOMAR
Trata de determinar la trayectoria
cronológica a lo largo de la que puede
crecer una economía, manteniendo una
utilización plena de su capacidad
productiva.
Si son constantes la propensión
marginal al ahorro (s) y la razón de
capital marginal / producción (k),
determinar la función de inversión que
se necesita para el crecimiento deseado.

30
2. El modelo de DOMAR
a. El cambio de la demanda agregada
es igual al cambio de la inversión por el
multiplicador 1/s:

dY 1 dI
=
dt s dt

31
2. El modelo de DOMAR
b. El cambio de la capacidad productiva
es igual al cambio de las existencias de
capital por el multiplicador 1/k:

dQ 1 dK
=
dt k dt

32
2. El modelo de DOMAR
c. El cambio de las existencias de capital
es igual a la inversión:

dK
=I
dt
dQ 1
= I
dt k

33
2. El modelo de DOMAR
d. Para la capacidad plenamente
empleada:

dQ dY
=
dt dt
1 1 dI
I=
k s dt

MATEMÁTICA PARA CONTADORES - Ing. Rubén Zárate Rojas, MSc 34


2. El modelo de DOMAR
dI s
− dt =
0
I k
s
ln I − t =
C
k
s
t
I =C e k

Para t=0, I(0)= C


s
t
I = I (0) e k

35
3.- El modelo de SOLOW
Examina las trayectorias de crecimiento en
equilibrio, tanto del capital como de la
mano de obra, con empleo completo.

a. La producción es una función


linealmente homogénea del capital de la
mano de obra, con beneficios, o
rendimientos a escalas constantes:

Y=f(K , L)…(1)

36
3.- El modelo de SOLOW

b. Una proporción constante s de la


producción se ahorra e invierte:
dK
= sY …(2)
dt

37
3.- El modelo de SOLOW

c. La oferta de mano de obra crece a un


índice constante r:

L = Lo e rt
…(3)

38
Deducir la ecuación diferencial en función
de la variable z=K/L
Y = f(K,L) …(1)
dK
= sY …(2)
dt
dK
= s f(K,L) …(3)
dt
L = Lo e rt
…(4)
dK
= s f(K,Loe ) …(5)
rt

dt 39
Deducir la ecuación diferencial en función
de la variable z=K/L

K = zL
K = zLo e rt

dK rt dz
= zr Lo e + L o e
rt

dt dt
dK rt  dz 
= Lo e  zr + 
dt  dt 
40
Deducir la ecuación diferencial en función
de la variable z=K/L
 dz 
( =
s f K,L oe

)
Lo e  zr + 
rt rt

dt 

 K  rt  dz 
Lo e s f  =
rt
,1 Lo e  zr + 
 dt 
rt
 Lo e 
 K  dz
sf  rt
,1=
 zr +
 Lo e  dt

dz
= s f ( z,1) − zr
dt
41

También podría gustarte