Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Aplicaciones Edo

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

Ecuaciones Diferenciales Ordinarias

Aplicaciones
Karina Malla Buchhorsts
Departamento de Matematicas ? UCN
marzo de 2013

Indice
1. Aplicaciones: Mecanica
2. Aplicaciones: razon de cambio
2.1. Mezclas . . . . . . . . . . .
2.2. Crecimiento y decrecimiento
2.3. Crecimiento de poblacion . .
2.4. Enfriamento . . . . . . . . .

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

.
.
.
.

5
5
6
7
9

1.

Aplicaciones: Mecanica

Consideraremos en esta seccion algunas aplicaciones de las Leyes de Newton, que nos permitiran aplicar nuestros
conocimientos de ecuaciones diferenciales de primer orden.
La cantidad de movimiento es igual al producto de la masa del cuerpo por su velocidad, esto es, m(t)v(t). En la
mayora de las aplicaciones la masa permanece constante, luego la cantidad de movimiento es mv(t). La mecanica clasica
trata del movimiento de objetos grandes comparados con un a tomo y lento comparado con la velocidad de la luz.
Leyes de Newton
1. Cuando un cuerpo no esta sujeto a ninguna fuerza externa, se mueve a velocidad constante.
2. Cuando un cuerpo esta sujeto a una o mas fuerzas externas, la razon de cambio con respecto al tiempo de la cantidad
de movimiento del cuerpo es igual a la suma vectorial de las fuerzas externas F que actuan sobre e l. Es decir,
dx
d
(mv(t)) = F(t, x, )
dt
dt

(1)

dx
es su velocidad.
dt

donde x(t) es la posicion del cuerpo en el instante t y

3. Cuando dos cuerpos actuan recprocamente, la fuerza que el primero ejerce sobre el segundo es igual en magnitud,
pero de direccion opuesta, a la fuerza que el segundo cuerpo ejerce sobre el primero.
La segunda Ley de Newton, que se aplica en un marco de referencia en el que un cuerpo no perturbado se mueve a
velocidad constante, permite formular las ecuaciones del movimiento de un cuerpo. Si a = dv/dt es la aceleracion del
cuerpo en el instante t, entonces la ecuacion (1) es
ma = F
(2)
La ecuacion (1) o (2) es una ecuacion diferencial de primer orden en v siempre que la fuerza F no dependa de x. En
esta seccion se usaran los sistemas de unidades: el sistema ingles, el sistema centmetro-gramo-segundo (cgs) y el sistema
metro-kilogramo-segundo (mks). las diferentes unidades de los tres sistemas se presentan en la siguiente tabla:
x
fuerza
masa
distancia
tiempo
aceleracion

Sistema ingles
libra
unidad de masa
pie
segundo
pie/seg 2

Sistema cgs
dina
gramo
centmetro
segundo
cm/seg 2

Sistema mks
newton
kilogramo
metro
segundo
m/seg 2

Recordemos que la fuerza de atraccion gravitacional que ejerce la tierra sobre un cuerpo se llama peso del cuerpo y se
representa por w. Para un cuerpo en cada libre (sin resistencia del aire), la u nica fuerza que actua sobre e l es su peso, de
modo que e sta es la fuerza resultante. La aceleracion es la debida a la gravedad, se representa por g y es aproximadamente
32 pies/s2 en el sistema ingles, 980 cm/s2 en el sistema cgs y 9,8 m/s2 en el sistema mks. La segunda ley de Newton
F = ma genera la relacion w = mg, es decir
w
(3)
m=
g
Para la cada de un cuerpo no se conoce una ley general que exprese la resistencia R del aire. En algunos casos la magnitud
de la resistencia del aire es R = kv, mientras que en otros es R = kv 2 , donde la constante de proporcionalidad depende
de las circunstancias.
Para un cuerpo en movimiento rectilneo, es decir, en movimiento a lo largo de una recta L como se muestra en la
figura (), escogemos un punto de referencia fijo como el origen O, una direccion fija como positiva y una unidad de
distancia, entonces la coordenada x de la posicion del cuerpo C a partir del origen O indica la distancia o desplazamiento
de C. La velocidad instantanea de C es v = dx/dt y la aceleracion instantanea de C es
a=

dv dx
dv
dv
=
=v
dt
dx dt
dx

Las fuerzas, desplazamientos, velocidades y aceleraciones en el sentido positivo de L son cantidades positivas; mientras que las que estan en sentido negativo son cantidades negativas. De esta forma la segunda ley de Newton ma = F para
el movimiento de C es una de las siguientes expresiones:
dv
dt
d2 x
m 2
dt
dv
mv
dx
m

(4)

(5)

(6)

donde F es la fuerza resultante que actua sobre el cuerpo.


Para resolver un problema de mecanica u otra aplicacion con ecuaciones diferenciales, es necesario tener en cuenta
dos etapas: una primera etapa de formulacion de un modelo matematico que interprete la situacion y una segunda de
resolucion de la ecuacion diferencial involucrada en el problema.
Ejemplo 1 Consideremos un objeto de masa m que cae bajo la influencia de la fuerza de gravedad y con una velocidad
inicial v0 dirigida hacia abajo. Si la fuerza de gravedad es constante y la resistencia del aire es proporcional a la
velocidad del ,objeto, determine la ecuacion del movimiento del objeto.
Formulacion del problema:
El movimiento del objeto es a lo largo de un eje vertical, sobre el cual se escoge como origen el punto donde el objeto
fue soltado y se denota con x(t) la distancia que el objeto a recorrido en el instante t, como se muestra en la siguiente
figura:
O
2

C
1

Las fuerzas que actuan sobre el objeto son: la fuerza de gravedad F1 = mg, donde g es la constante gravitacional cerca
de la tierra; la fuerza opuesta por la resistencia del aire F2 = kv(t) = kx(t), donde k es la constante de proporcionalidad.
Por lo tanto, la fuerza resultante que actua sobre el objeto es
F = F1 + F2 = mg kv(t).
Entonces, aplicando la segunda ley de Newton y considerando que la velocidad inicial es v0 , se obtiene el modelo matematico siguiente
dv
m
= mg kv, v(0) = v0 .
(7)
dt
Solucion. El modelo matematico (7) del problema planteado es un problema de valor inicial, con una ecuacion diferencial
de primer orden en que corresponde a una ecuacion de variables separables puesto que se puede escribir en la forma


m
dv = dt
mg kv
integrando y considerando la condicion inicial, se obtiene la solucion
mg 
mg  kt/m
e
v(t) =
+ v0
k
k

(8)

La ecuacion (8) nos da la velocidad del objeto en el instante t, si queremos conocer la distancia recorrida en un tiempo t,
integramos dx/dt = v(t) para obtener
Z
mg
m
mg  kt/m
e
x(t) = v(t)dt =
t
v0
+C
k
k
k
y considerando la condicin x(0) = 0, resulta la constante
C=

mg 
m
v0
k
k
3

de donde se obtiene la ecuacion del movimiento del objeto


mg
mg 
m
x(t) =
v0
(1 ekt/m )
t+
k
k
k

(9)

Observacion 1 El Problema anterior se puede resolver para diferentes valores de las constante m, k y v0 . Ademas, se
puede responder a una variedad de preguntas sobre la cada de un objeto con resistencia del aire proporcional a su
velocidad. Por ejemplo, cuando t tiende a la velocidad (8) tiende a llamada velocidad lmite del objeto.
t=0
x(t)
F2=kv(t)
t
F1=mg

Ejemplo 2 Un objeto cae desde el reposo hacia el mar. El peso del objeto es de 64 lb y la resistencia del aire igual a 1/2
de su velocidad v. Si 4 segundos despues cae al agua y la resistencia de e sta es igual a 1/4 del cuadrado de la velocidad,
encuentre la velocidad del objeto:
1. En el aire y en el momento de tocar el agua.
2. En el agua.
Solucion. Tal como el problema anterior, se toma el eje vertical positivo hacia abajo y elegimos el origen en el punto
en que comienza la cada. El problema tiene dos partes: a) Cuando el objeto esta en el aire, y b) cuando el objeto esta en
el agua.
1. Antes de caer al agua, las fuerzas que actuan sobre el objeto son: F1 = mg = 64 lbs, que actua hacia abajo y, por
lo tanto, es positiva. F2 = 1/2v que actua hacia arriba y, por lo tanto, es negativa, es decir, F2 = 1/2v.
Tomando g = 32, se tiene que m = 2 y de acuerdo a la segunda ley de Newton, con F = F1 + F2 , y la condicion
inicial v(0) = 0, se obtiene el PVI
dv
1
2
= 64 v,
v(0) = 0
(10)
dt
2
La solucion de la ecuacion diferencial de variables separables del PVI (10), es
1
t
v(t) = 128 + C e 4
con la condicion inicial v(0) = 0 obtenemos c = 128, luego la velocidad del objeto en el aire es

1
t
v(t) = 128 1 e 4 , 0 t 4
En particular, cuando t = 4, obtenemos
v(4) = 128(1 e1 ),
que es la velocidad en el instante en que el objeto toca el agua.
4

2. En el agua, las fuerzas que actuan sobre el objeto son: su peso F1 = 64 libras, que actua hacia abajo, por lo tanto,
1
es positiva y F2 = v 2 , que actua hacia arriba, por lo tanto, es negativa. Considerando que en el instante t = 4 la
4
velocidad v(4) = 128(1 e1 ) es aproximadamente igual a 81, de la segunda ley de Newton, se tiene el PVI
2

dv
1
= 64 v 2 ,
dt
4

v(4) = 81

(11)

La ecuacin diferencial del PVI (11), separando las variables, es


1
dv
= dt
v 2 16
8
integrando, por fracciones parciales, se obtiene
v(t) =

16(1 + C e4t
),
1 C e4t

1 ce Finalmente, haciendo uso de la condicion inicial, obtenemos que la constante C es aproximadamente 0,8, por
lo tanto, la velocidad del objeto en el agua es
v(t) =

2.

16(1 + 0,8 e4t


),
1 0,8 e4t

Aplicaciones: razon de cambio

dx
representa la razon de cambio de
Si x(t) representa la magnitud de una cantidad presente en el tiempo t, entonces
dt
esa cantidad con respecto del tiempo. Consideraremos en esta seccion algunas situaciones que corresponden a problemas
de razon de cambio, esto es: mezclas, crecimiento y decrecimiento, crecimientos de poblacion enfriamiento.

2.1.

Mezclas

Supongamos que una sustancia S fluye hacia una cierta mezcla de un recipiente, con una cierta rapidez, y la mezcla
se mantiene uniforme por agitacion. Ademas, la mezcla uniforme debera salir del recipiente y pasar a otro. Se busca
determinar la cantidad de sustancia S presente en la mezcla en el tiempo t.
dx
Representaremos por x(t) a la cantidad de S en el tiempo t y as
corresponde a la razon de cambio de la cantidad
dt
de sustancia, de esta manera se tiene la ecuacion
dx
= Rapidez de ENTRADA Rapidez de SALIDA,
dt

(12)

que sirve como modelo para resolver los siguientes problemas.


Ejemplo 3 Un deposito contiene 50 galones de salmuera en el que hay disueltas 75 libras de sal. A partir de un instante
t = 0, comienza a fluir en e l salmuera que contiene una libra de sal por galon, a una velocidad de 2 galones por minuto.
Simultaneamente la mezcla sale del deposito a razon de un galon por minuto.
1. Que cantidad de sal hay en el deposito en un tiempo t?.
2. En que instante habra 125 libras de sal disueltas en el deposito?.
Solucion
Sea x(t) la cantidad de sal presente en el deposito, como dijimos anteriormente, la rapidez de entrada es
ENTRADA = (2 gal/min)(1 lib/gal) = 2 lib/min.
Para encontrar la rapidez de salida, debemos considerar que siendo 50 galones la cantidad original de salmuera, t minutos
despues aumenta a 50 + (2 1)t = 50 + t galones, por lo que la concentracion en t minutos es
x(t)
lib/gal,
50 + t
5

en consecuencia, la rapidez de salida es



SALIDA = (1 gal/min)

x(t)
lib/gal
50 + (2 1)t


=

x(t)
lib/min.
50 + t

y por (12), tenemos la ecuacion


dx
x(t)
=2
dt
50 + t
como el deposito contena originalmente 75 libras de sal, tenemos la condicion inicial x(0) = 75, y con ello el modelo
matematico del problema es el
dx
x(t)
=2
,
x(0) = 75
dt
50 + t
1. La ecuacion diferencial del PVI, es decir
dx
x(t)
+
=2
dt
50 + t
es lineal, luego su solucion se obtiene aplicando metodos adecuados, en consecuencia
x(t) = 50 + t +

c
50 + t

con c una constante arbitraria. Haciendo uso de la condicion inicial, obtenemos c = 1250. De esta manera se tiene
que la cantidad de sal en un instante t cualquiera es
x(t) = 50 + t +

1250
50 + t

2. El deposito contendra 125 libras de sal si x(t) = 125, luego se tiene la ecuacion en t
125 = 50 + t +

1250
50 + t

Las posibles soluciones de esta ecuacion son: t = 64,04 y t = 39,04, como t debe ser positivo, conclumos que
el deposito contiene 125 libras de sal al cabo de 64,04 minutos.

2.2.

Crecimiento y decrecimiento

Muchos fenomenos que involucran crecimiento o decrecimiento (descomposicion) tienen como modelo matematico
al problema de valor inicial
dx
= kx(t),
x(t0 ) = x0
(13)
dt
Por ejemplo. En fsica (13) es un modelo para aproximar la cantidad restante de una substancia que se desintegra radiactivamente o para determinar la temperatura de un cuerpo que se enfra; en biologa se observa que la rapidez con que
ciertas bacterias se multiplican es proporcional al numero de bacterias presentes en un instante t, es decir, se puede utilizar
(13); finalmente, en qumica (13), sirve de modelo para describir la cantidad restante de una substancia durante ciertas
reacciones.
Ejemplo 4 La rapidez con la que se desintegran nucleos radiactivos es proporcional al numero de nucleos que estan
presentes en una muestra dada. La mitad del numero original de nucleos radiactivos ha experimentado la desintegracion
en un periodo de 1500 anos.
1. Que porcentaje de nucleos radiactivos originales continuaran despues de 4500 anos?.
2. En cuantos anos quedara solamente un decimo del numero original?.
Solucion En primer lugar buscaremos un modelo matematico del problema. Si denotamos por x(t) a la cantidad de
dx
nucleos radiactivos presentes despues de t anos, entonces
representa la rapidez con la que se desintegran los nucleos.
dt
Dado que esta desintegracion es proporcional a la cantidad presente, tenemos la ecuacion
dx
= kx(t),
dt
6

dx
donde k es una constante de proporcionalidad. Como x es positiva y decreciente,
< 0, entonces la ecuacion anterior
dt
se puede escribir en la forma
dx
= kx(t),
dt
con k > 0. Si denotamos por x0 la cantidad inicialmente presente, entonces tenemos la condicion inicial x(0) = x0 y con
ello tenernos que el modelo matematico del problema es el PVI, del tipo (13),
dx
= kx(t),
dt

x(0) = x0

(14)

La ecuacion del PVI (14) es de variables separables y lineal homogenea; aplicando las tecnicas conocidas, por ejemplo (),
se obtiene
x(t) = c ekt
y al aplicar la condicion inicial se encuentra c = x0 , de donde la solucion del PVI () es
x(t) = x0 ekt
Para determinar la constante de proporcionalidad k utilizamos la condicion dada en el problema x(1500) =
 1
encuentra ek = 21 1500 . Reemplazando ek en la solucion de la ecuacion se tiene que

1
2 x0 ,

se

"  1 #t
1 1500
x(t) = x0
2
de manera que el numero de nucleos radiactivos presentes en el tiempo t es
t
  1500
1
x(t) = x0
2
1. Despues de 4500 anos quedaran x(4500) = 18 x0 nucleos radiactivos, o bien, 12.5 % de los nucleos originales.
2. Para saber en cuantos anos quedara un decimo del numero original de nucleos, hacemos x(t) =
1
=
10

1
10 x0 ,

de donde

t
  1500
1
2

Resolviendo para t se obtiene


t = 1500

ln 10
ln 2

que es aproximadamente 4985 anos.

2.3.

Crecimiento de poblacion

Consideremos una poblacion (humana, de bacterias, etc.) y representemos por x(t) a la cantidad de individuos en el
tiempo t. Si suponemos que la tasa de crecimiento es proporcional a la poblacion presente, que la tasa de mortalidad es
cero y que en un tiempo t0 la poblacion es x0 , entonces un modelo matematico para el problema es
dx
= kx(t),
dt

x(t0 ) = x0 ,

(15)

dx
donde k > 0 es la constante de proporcionalidad, que es positiva dado que la poblacion x y su crecimiento
son
dt
positivos.
La ecuacion diferencial del PVI (15) es lineal y de variables separables, luego es facil de resolver. En efecto, su
solucion es
x(t) = c ekt ,
y haciendo uso de la condicion inicial, obtenemos c = x0 ek(tt0 ) de manera que la solucion del PVI es
x(t) = x0 ek(t0 )
Esta ley de crecimiento de poblacion, bastante aceptable para intervalos cortos de tiempo, se llama ley Maltusiana.
7

(16)

Ejemplo 5 Suponga que la poblacion de una cierta comunidad se duplica cada 5 anos y que crece como se describio anteriormente.
1. Cuanto demorara en triplicarse?.
2. Si despues de 3 anos la poblacion es de 10000 habitantes. Cual era la poblacion original?
Solucion.
Si la poblacion original de la comunidad mencionada es x0 , entonces la condicion inicial del problema es x(0) = x0
y el modelo matematico para la poblacion es
dx
= kx(t),
dt

x(0) = x0

De acuerdo con (16) la solucion a este PVI es x(t) = x0 ekt . Si la poblacion se duplica en 5 anos, entonces x(5) = 2x0 ,
es decir
2x0 = x0 e5k
de donde obtenemos la constante de proporcionalidad k =

ln 2
5

de manera que la poblacion en un ano t cualquiera es

x(t) = x0 e

ln 2
5 t

1. La poblacion se triplica si x(t) = 3x0 , es decir si


x0 e
resolviendo para t, obtenemos t =

5 ln 3
ln 2

ln 2
5 t

= 3x0

= 7,9. Esto es, la poblacion se triplica en 7,9 anos.

2. Si en 3 anos la poblacion es de 10000 habitantes, entonces tenemos la condicion x(3) = 10000. Reemplazando esta
condicion en x(t), se tiene
3 ln 2
10000 = x0 e 5
de donde obtenemos x0 = 6597,539. Luego, la poblacion original era de 6597 habitantes.
Una representacion mas realista del crecimiento de una poblacion se logra con la ecuacion diferencial
dx
= ax bx2 ,
dt

(17)

donde a y b son constantes positivas. En la mayora de los casos la constante b es muy pequena comparada con a,
en consecuencia, para una poblacion x suficientemente pequena, el termino ax predomina y as la poblacion crece
en forma muy rapida respecto del tiempo. Por otro lado, cuando x llega a ser suficientemente grande, predomina
el termino bx2 y se tiene una disminucion de la rapidez de crecimiento. La ecuacion (17) se llama ecuacion
logstica de crecimiento. Es claro que esta ecuacion es una ecuacion diferencial de variables separables y tambien
una ecuacion de Bernoulli, por lo que su resolucion no trae mayores dificultades.
Ejemplo 6 Suponga que un estudiante portador de un virus regresa a un campus universitario aislado, que tiene 1000
estudiantes. Si la rapidez con que el virus se propaga es proporcional al producto del numero de estudiantes contagiados
por los no contagiados y se observa que a los 4 das hay 50 contagiados, entonces determine el numero de estudiantes
contagiados a los 6 das.
Solucion.
Si denotamos por x(t) al numero de alumnos contagiados con el virus, entonces 1000 x corresponde a los no
contagiados y la rapidez con que el virus se propaga es dada por la ecuacion
dx
= kx(1000 x)
dt
y como el portador es uno, se tiene la condicion inicial x(0) = 1. Luego, el modelo matematico del problema es
dx
= kx(1000 x),
dt

x(0) = 1

Para resolver este PVI consideramos que la ecuacion es de variables separables, separando las variables e integrando, se
obtiene
ln x ln(1000x) = 1000tk + c,
de la condicion inicial, con t = 0 y x = 1, obtenemos la solucion
x(t) =

1000 e1000tk
999 + e1000tk

Para determinar la constante de proporcionalidad consideramos que x(4) = 50 y obtenemos k = 0,0009906. Finalmente,

reemplazando k en la solucion se obtiene la funcion que indica el nmero


de contagiados en un tiempo t, esto es
x(t) =

1000 e0,9906k
999 + e9906k

El numero de contagiados a los 6 das es x(6) = 276 estudiantes.

2.4.

Enfriamento

La ley de Newton del enfriamiento dice que, en un cuerpo que se esta enfriando, la rapidez con que cambia la temperatura T (t) es proporcional a la diferencia entre la temperatura del cuerpo y la temperatura constante M , del medio que
lo rodea. Es decir,
dT
= K(T M )
(18)
dt
donde K es una constante de proporcionalidad.
Ejemplo 7 El plasma sanguneo se almacena a una temperatura de 40o y antes de ser usado debe estar a 90o de temperatura. Si se coloca el plasma dentro de un horno a 120o , demora 45 minutos en alcanzar los 90o . Suponga que la ley de
enfriamiento de Newton es aplicable y que la constante K de enfriamiento es independiente de la temperatura del horno.
Cuanto demorara el plasma en estar en condiciones de ser usado si la temperatura del horno es: 1) 100o y 2) 140o ?.
Solucion Si el plasma se almacena a una temperatura de 40o , T (t) es la temperatura en un instante t y la temperatura
ambiental (en el horno) es M = 120o , entonces, de acuerdo con (18) el modelo matematico que rige el fenomeno es dado
por el PVI
dT
= K(T 120),
T (0) = 40.
dt
La ecuacion diferencial del PVI es de variables separables, pues se puede escribir en la forma
dT
= Kdt
T 120
en consecuencia, integrando se obtiene ln |T 120| = Kt + c y usando la condicion inicial T (0) = 40, se obtiene la
constante c = ln 80, luego
ln |T 120| = Kt + ln 80
El valor de la constante K es independiente de la temperatura dci horno, luego consideramos T (45) = 90 para obtener
K = 0,0218 y reemplazando en la expresion anterior se logra la solucion del PVI
ln |T 120| = 0,0218t + ln 80
De donde obtenemos que la temperatura del plasma a una temperatura ambiental M es
ln |T M | = 0,0218t + ln 80
1. Si la temperatura del horno es 100o , entonces la temperatura del plasma en un instante t es dada por ln |T 100| =
0,0218t + ln 80.
Ahora, el plasma logra los 90o , si y solo si, T (t) = 90, de donde ln |90100| = 0,0218t+ln 80, luego t = 82,19.
Es decir, el plasma esta en condiciones de ser usado a los 82,19 miriutos
2. Si la temperatura del horno es 140o , entonces la temperatura del plasma en un instante t es dada por ln |T 140| =
0,0218t + ln 80. Ahora, el plasma logra los 90, para ser usado, si y slo si ln |90 140| = 0,0218t + ln 80, de
donde se obtiene t = 31,79. Es decir, el plasma esta en condiciones de ser usado a los 31,79 minutos.

También podría gustarte