Aplicaciones
Aplicaciones
Aplicaciones
Ecuaciones Diferenciales
Gloria Aguilar
Natalia Boal
Carmelo Clavero
Francisco Gaspar
Universidad de Zaragoza
1
Selección de ejemplos obtenidos de los libros de Nagle-Saff, Braun, Finzio-Ladas, Campbell,
entre otros.
Trayectorias ortogonales
dx ∂F
∂y
ii) Usando que la pendiente de una curva ortogonal (perpendicular) a una curva es la inversa
de la pendiente de la curva dada, demuestra que las curvas ortogonales a la familia F (x, y) = k
satisfacen la ecuación diferencial
∂F ∂F
(x, y)dx − (x, y)dy = 0.
∂y ∂x
iii) Utilizando la ecuación diferencial precedente, demuestra que las trayectorias ortogonales
de la familia de circunferencias x2 + y2 = k son rectas que pasan por el origen.
2. Usa el método del problema anterior para encontrar las trayectorias ortogonales de cada
familia de curvas dada. Dibuja conjuntamente la familia de curvas y sus trayectorias ortogonales.
Problemas de mezclas
3. En los problemas de mezclas se quiere calcular la cantidad de una sustancia, x(t), que hay
en un tanque en cada instante de tiempo t. Usando que la derivada de x respecto a t expresa la
razón de cambio de la sustancia presente en el tanque, se cumple la relación
dx
= velocidad de entrada - velocidad de salida .
dt
Dada la velocidad a la que un fluido que contiene la sustancia entra en el tanque y la concen-
tración de la sustancia, se cumple la relación
Por ejemplo, en un tanque de 1000 l (litros) una solución salada de salmuera empieza a fluir a
velocidad constante de 6 l/min. La solución se mantiene bien agitada y fluye hacia el exterior
con velocidad de 6 l/min. Sabiendo que la concentración de salmuera que entra es de 1 Kg/l,
2
calcula cuándo la concentración de sal será de 1/2 Kg/l. Si la velocidad de salida es 5 l/min,
determina la concentración de sal en función del tiempo.
4. La corriente sanguı́nea lleva un medicamento hacia el interior de un órgano a razón de
3cm3 /seg y sale de él a la misma velocidad. Se sabe que el volumen del órgano es 125 cm3
y la concentración del medicamento que entra es de 0.2 g/cm3 . ¿Cuál es la concentración del
medicamento en el órgano si inicialmente no habı́a vestigio alguno del medicamento? ¿Cuándo
la concentración será de 0.1 g/cm3 ?
5. El agua del rı́o Aguadulce fluye hacia el lago Magdalena a razón de 300 gal/min (1 galón
es una medida inglesa de capacidad que equivale a 4.5459 litros; en Estados Unidos equivale a
3.7853 litros). El lago Magdalena contiene aproximadamente 100 millones de galones de agua.
La fumigación de los naranjales cercanos ha ocasionado que la concentración de plaguicidas en
el lago llegue a ser de 35 partes por millón. Si se suspende la aplicación de plaguicidas, ¿cuánto
tiempo transcurrirá antes de que la concentración de los mismos en el lago esté por debajo de
10 partes por millón? (Se supone que el rı́o no contiene plaguicidas y que el volumen del lago
es constante).
6. El modelo malthusiano de crecimiento de una población p(t) supone que la tasa de creci-
miento es proporcional a la población presente. Sabiendo que la población de Estados Unidos
en 1790 era de 3.93 millones y en 1800 de 5.31 millones, usa el modelo anterior para conocer
la población en función del tiempo. Este modelo supone que la tasa de mortalidad es nula,
que desde luego es errónea. Parece natural pensar que la tasa de mortalidad natural también
es proporcional al tamaño de la población. No obstante, debido a otros factores de mortalidad
(desnutrición, enfermedades, crı́menes violentos, etc), se puede suponer que la tasa de morta-
lidad es proporcional al número de interacciones bipartitas. Para una población de tamaño p,
existen p(p − 1)/2 interacciones de este tipo. Prueba que con esta hipótesis, el PVI que rige el
modelo tiene la forma
dp
= ap − bp2, p(0) = p 0 .
dt
Esta ecuación se conoce con el nombre de ecuación logı́stica. Calcula la población en función
del tiempo. En el caso particular de que la población en 1790 sea de 3.93 millones, en 1840 de
17.07 millones y en 1890 de 62.95 millones, determina la solución usando el modelo logı́stico.
7. Un modelo de crecimiento de población que se utiliza en predicciones actuariales se basa en
la ecuación de Gompertz
dp
= p(a − b ln p),
dt
con a y b constantes reales.
i) Resuelve la ecuación diferencial para calcular p(t).
ii) Si p(0) = p0 > 0, proporciona una fórmula para p(t) en términos de a, b, p0 y t.
iii) Describe el comportamiento de p(t) cuando t → ∞, considerando los casos b > 0 y b < 0.
8. Una bola de nieve se derrite de forma que la razón de cambio de su volumen es proporcional
al área de su superficie. Si el diámetro inicial de la bola es de 4 pulgadas y al cabo de 30
min. es de 3 pulgadas, ¿cuándo será de 2 pulgadas?, ¿cuando desaparecerá la bola de nieve?
Si suponemos que la bola se derrite tal que la razón de cambio de su diámetro es proporcional
al área de su superficie, con los mismos datos, ¿cuándo será su diámetro de 2 pulgadas?; en
términos matemáticos ¿cuándo desaparecerá la bola de nieve?.
9. Uno de los métodos más precisos para determinar la edad de restos arqueológicos es el
método del carbono 14 (C 14 ), basado en que para cualquier organismo vivo una proporción
3
constante de átomos de carbono está formada por el isótopo radiactivo C 14 . La proporción
permanece prácticamente constante durante toda la vida y cuando el organismo muere el C 14
sigue su proceso de desintegración, con lo cual la proporción disminuye. Un modelo simple para
describir el fenómeno supone que la cantidad de átomos que se desintegran es proporcional a la
cantidad de átomos presentes, es decir, se cumple la ecuación diferencial
dN
= −λN (t),
dt
siendo N (t) la cantidad de C14 en una muestra en el tiempo t. Suponiendo que N (0) = N0
calcula la solución del PVI. A R(t) = λN (t) se le llama tasa de desintegración y R(0) es
la tasa original de desintegración que coincide con la tasa de desintegración de la materia viva.
Sabiendo que la edad media del C14 (la edad media es el tiempo que debe transcurrir para que
se desintegre la mitad de la sustancia) es aproximadamente de 5600 años, resuelve los siguientes
casos:
i) El nivel de carbón vegetal extraı́do en las grutas de Lascaux (Francia) en 1950 dio una
medida de 0.91 desintegraciones por minuto y gramo, mientras que la materia viva da 6.68
desintegraciones. Calcula la época en que las grutas estuvieron habitadas.
ii) En la excavación en Ninpur (Babilonia), en 1950 el carbón vegetal de una viga dio
4.09 desintegraciones por minuto y gramo. Suponiendo que este carbón se formó en la época
Hamurabi, calcula una fecha probable de la sucesión de Hamurabi.
iii) En una cueva de Sudáfrica se encontró un cráneo humano junto a los restos de una fogata.
Los arqueólogos creen que la edad del cráneo es igual a la de la fogata. Sabiendo que solo un
2% de la cantidad original de C 14 queda en la madera, calcula la edad aproximada del cráneo.
10. Sea T (t) la temperatura interior de un edificio para el que la razón de cambio de temperatura
es la diferencia entre la razón a la que aumenta y la razón a la que disminuye. Suponemos
que afectan tres factores en la temperatura. El primero (calor producido por las personas,
luces y máquinas) incrementa la temperatura a razón de H(t). El segundo es el calentamiento
(enfriamiento) producido por la calefacción (aire acondicionado) a razón de U (t). El tercer factor
es el efecto de la temperatura exterior M (t), para el que la razón de cambio es proporcional a
la diferencia entre la temperatura exterior y la interior (ley de Newton del enfriamiento).
Recogiendo todo, la ecuación diferencial que modela el fenómeno es
dT
= k(M (t) − T (t)) + H(t) + U (t),
dt
con k constante real que depende de las propiedades del edificio (número de puertas y ventanas,
aislamiento, material, etc); al valor 1/k se le llama constante de tiempo del ediftcio. Resuelve
la ecuación diferencial para calcular la temperatura en función del tiempo.
i) Si la temperatura exterior es constante M = 0 y se cumple H = U = 0, escribe la solución
sabiendo que T (t0) = T0 (esto refleja como var´ıa la temperatura del edificio).
ii) Supongamos que H es constante H0 , U = 0 (no hay calefacción o enfriamiento) y la
temperatura exterior está dada por M (t) = M0 − B cos ωt, donde B es una constante real
positiva y ω = π/12 (onda senoidal de perı́odo 24 horas, con mı́nimo en t = 0 (medianoche) y
máximo en t = 12 (mediodı́a)). Calcula la solución sabiendo que a medianoche la temperatura
es T0 (esto refleja como varı́a la temperatura en primavera u otoño cuando no hay calefacción
ni aire acondicionado).
iii) En el apartado anterior ii) supongamos que hay un termostato para controlar la tem-
peratura del edificio en relación con la temperatura deseada Td . Suponiendo que el calen-
tamiento o enfriamiento suministrado es proporcional a la diferencia de temperatura, es decir,
4
U (t) = ku(Td T− (t)), calcula T (t) siendo de nuevo T (0) = T0 (esto refleja como var´ıa la tem-
peratura en verano (aire acondicionado) y en invierno (calefacción)).
11. En una mañana, mientras las personas trabajan, la calefacción mantiene la temperatura
interior a 210 C. A mediodı́a se apaga y se sabe que la temperatura exterior es de 120 C durante
la tarde. Sabiendo que la constante de tiempo del edificio es de 3 horas, ¿cuándo la temperatura
será de 160 C? Si se dejan abiertas algunas ventanas, la constante de tiempo se reduce a 2 horas;
¿cuándo se alcanzarán ahora los 160 C?.
12. La ley de Stefan de la radiación establece que la razón de cambio de la temperatura de
un cuerpo que se encuentra a T grados Kelvin en un medio a M grados Kelvin, es proporcional
a M 4 − T 4 . Plantea y resuelve la ecuación diferencial. Explica por qué las leyes de Newton
y Stefan son aproximadamente iguales cuando T se acerca a M y M es constante. (Ayuda:
factoriza M 4 − T 4).
13. Un vino tinto se saca de la bodega, que es un lugar fr´ıo a 100C, y se deja reposar en un
cuarto con temperatura de 230 C. ¿En qué momento la temperatura del vino será de 180 C si
transcurren 10 minutos para alcanzar los 150C?
Mecánica Newtoniana
14. A un objeto de masa m se le aplica una velocidad inicial dirigida hacia abajo v0 y cae bajo
la acción de la fuerza de la gravedad, que se supone constante, y la fuerza debida a la resistencia
del aire que es proporcional a la velocidad. Calcula la ecuación del movimiento del objeto.
15. Un paracaidista cuya peso (masa) es de 100 Kg se deja caer desde un helicóptero situado a
3000 metros. Suponiendo que la resistencia del aire es proporcional a la velocidad con constante
k1 = 20 kg/seg si el paracaı́das está cerrado y k1 = 100 kg/seg si está abierto y que el paracaı́das
se abre transcurridos 30 segundos, ¿cuándo llegará a la superficie?. Si el paracaı́das se abre al
cabo de 1 minuto, ¿cuándo llegará el paracaidista a la superficie?.
16. Vuelo de un cohete. Un cohete con masa inicial de m0 Kg se lanza verticalmente desde
la superficie de la Tierra. El cohete expele gas a razón de α Kg/seg y a una velocidad constante
de β m/seg relativa al cohete. Suponiendo que el campo gravitacional es constante de g kg/seg2,
la segunda ley de Newton da lugar a la ecuación
dv
(m0 − αt) − αβ = −g(m0 − αt),
dt
donde v = dx/dt es la velocidad del cohete, x es su altura respecto de la superficie de la Tierra
y m0 αt− es la masa del cohete a los t segundos del lanzamiento. Sabiendo que la velocidad
inicial es cero, resuelve la ecuación anterior para calcular la velocidad y altura del cohete para
0 ≤ t ≤ m0/α.
5
donde k es la constante de reacción y a(0), b(0) son conocidos. Resuelve en el caso que
a(0) = 3, b(0) = 1 mol/l y k = 1 l/mol.s. Interpreta el resultado.
18. Supongamos que dos tanques A, B de volúmenes V1 , V2 están separados por una frontera,
a través de la cual una sustancia fluye proporcionalmente (con constante k) de la diferencia de
concentraciones C1 − C2 . Calcula la concentración en función del tiempo en cada tanques.
19. Consideremos un sistema mecánico que consiste de un resorte (muelle) en espiral suspendido
de un soporte rı́gido con una masa sujeta al extremo (ver figura 1). Para analizar el fenómeno
usamos dos leyes: ley de Hooke y la segunda ley de Newton. La ley de Hooke establece
que el resorte ejerce una fuerza de restitución opuesta a la dirección de alargamiento del resorte
(peso=kl donde k es la constante de restitución y l es el alargamiento hasta la posición de
equilibrio). Para aplicar la ley de Newton, hay que tener en cuenta que sobre la masa m actúa
la fuerza de la gravedad (mg, g vale 9.8m/seg 2 ), la fuerza de restitución (− kx − mg donde x
es el alargamiento), la fuerza de amortiguación que es proporcional a la velocidad de la masa
(− b(dx/dt) con b la constante de amortiguación) y las fuerzas externas (f (t)). Ası́, la ecuación
diferencial que resulta es
d 2x dx
m 2 + b + kx = f (t).
dt dt
Si b = 0 el sistema es no amortiguado; en otro caso se dice amortiguado. Si f (t) = 0 se dice
libre; en otro caso se dice forzado. Sabiendo esto, resuelve los siguientes casos.
i) Una masa de 5 Kg se sujeta de un resorte y ocasiona un estiramiento de 0.5 m. A
continuación se tira de la masa 0.1 m abajo del punto de equilibrio y se le aplica una velocidad
dirigida hacia arriba de 0.1 m/seg. Calcula la ecuación del movimiento armónico simple de la
masa. ¿Cuando alcanzará la masa por primera vez su mı́nima altura después de haberse puesto
en movimiento?
ii) Una masa de 2 kg estira un resorte 49 cm al llegar al reposo en equilibrio. La constante
√
de amortiguación es 8 5N ·seg/m. Si la masa se tira 10 cm hacia abajo del punto de equilibrio
y se le aplica una velocidad de 2 m/seg dirigida hacia abajo, ¿cuál es el desplazamiento máximo
que alcanzará a partir de la posición de equilibrio?
iii) Una masa de 2 kg estira un resorte 20 cm al llegar al reposo en equilibrio. La constante
de amortiguación es 5N · seg/m. En t = 0 la masa se desplaza 5 cm hacia abajo y se suelta. En
el mismo instante se aplica una fuerza externa f (t) = 0.3 cos t N al sistema. Calcula la ecuación
de movimiento de la masa.
Sabiendo que la frecuencia de resonancia del sistema es γr /2π donde
.
k b2
γr = − ,
m 2m2
determina la frecuencia de resonancia del sistema. Se dice que un sistema está en resonancia
cuando es estimulado por una fuerza externa cuya frecuencia coincide con la frecuencia de
resonancia del sistema. En el caso de que la constante de amortiguamiento sea pequeña, si la
fuerza externa tiene una frecuencia cercana a la frecuencia del sistema, el sistema está sujeto a
grandes oscilaciones. Estas grandes vibraciones son las que preocupan a los ingenieros, puesto
que pueden ocasionar que las alas de los aviones se rompan, los puentes se desplomen o (menos
importante) que las copas de vidrio se hagan añicos.
20. Se considera un servomecanismo que modela a un piloto automático. Si y(t) es la dirección
real de la nave en el instante t y g(t) la dirección deseada, entonces el error o desviación es
6
e(t) = y(t) − g(t). La segunda ley de Newton, expresada en términos de momentos de las fuerzas
proporciona la ecuación
Iyjj(t) = −ke(t),
donde I es el momento de inercia del eje del timón y k es una constante positiva. Calcula el
error cuando g(t) = at con a constante, sabiendo que el eje del timón se encuentra inicialmente
en reposo. Supongamos ahora que para controlar las oscilaciones se suministra al eje del timón
una fuerza giratoria proporcional a ej (t), pero de signo opuesto. Prueba que ahora la ecuación
diferencial es
Iyjj(t) = −ke(t) − µej(t),
donde µ es una constante positiva. Determina ahora el error sabiendo que g(t) = a con a
constante y de nuevo el eje del timón se encuentra inicialmente en reposo.
21. Consideremos un circuito eléctrico elemental que consta de una fuerza electromotriz (baterı́a
o generador), una resistencia, una bobina y un condensador en serie (ver figura 2). Estos circuitos
se rigen por dos principios fı́sicos: la conservación de la carga eléctrica y la conservación de la
energ´ıa. Las leyes que formulan estos principios son las leyes de Kirchoff que establecen:
i) La corriente I que pasa a través de cada elemento de un circuito en serie debe ser la misma.
ii) La suma algebraica de los cambios instantáneos de potencial (caı́da de voltaje) en un
circuito cerrado debe ser cero.
En cada elemento la ca´ıda de voltaje es la siguiente:
a) Según la ley de Ohm, la caı́da de voltaje en una resistencia es proporcional a la corriente
I que pasa por ella, es decir, ER = RI donde R es la resistencia.
b) Según las leyes de Faraday y Lenz, la caı́da de voltaje en una bobina es proporcional
dI
a la razón de cambio instantánea de la corriente I que pasa por ella, es decir, EL = L donde
dt
L es la inductancia.
c) La caı́da de voltaje en un condensador es proporcional a la carga eléctrica q del conden-
1 1
sador, es decir, E = q donde C es la capacitancia y es la elastancia.
C
C C
dq
Se sabe que la carga y la corriente están relacionadas por I = . Ası́, la ecuación diferencial
dt
que rige el modelo es
d2q dq 1
L 2 +R + q = E(t),
dt dt C
donde E(t) es el voltaje aplicado por la fuerza electromotriz., o en términos de la corriente
d 2I dI + 1 I = dE(t)
L 2 + R .
dt dt C dt
A partir de esto calcula la corriente (amperios) y la carga (culombios) en los siguientes casos;
en todos ellos analiza la solución.
i) La resistencia es de 0.02 ohmios, la inductancia es de 0.001 henrios, la capacitancia es de
2 faradios, la fuerza electromotriz es E(t) = sen(100t) voltios y la corriente y carga inicial son
cero.
ii) La resistencia es de 120 ohmios, la inductancia es de 4 henrios, la capacitancia es de
2200−1 faradios, la fuerza electromotriz es E(t) = 10cos(2t) voltios y la corriente y carga inicial
son cero.
iii) La resistencia es de 10 ohmios, la inductancia es de 4 henrios, la capacitancia es de 0.01
faradios, la fuerza electromotriz es E(t) = E0cosγt voltios con E0, γ constantes y la corriente y
carga inicial son cero.
22. En el diseño de una planta de tratamiento de aguas residuales, se propuso el siguiente PVI
60 − H = 77.8Hjj + 19.42(Hj)2, H(0) = Hj(0) = 0,
7
donde H(t) es el nivel del fluido en una cámara de eyección y t es el tiempo en segundos. Usa
el método de Euler para aproximar H(t) en [0, 5].
8
donde y1 es la concentración de células, y2 es la concentración de gluco..., y3 es la concentracion de
ácido glucónico, y4 es la concentración de glucosa y bi , i = 1, · · ·, 5 son parámetros. Sabiendo que
b1 = 0.949, b2 = 3.439, b3 = 18.72, b4 = 37.51, b5 = 1.169, y1(0) = y4(0) = 5, y2(0) = y3(0) = 0,
resuelve el modelo en el intervalo [0, 9] usando el método de Euler e interpreta los resultados.
27. Sobre una superficie lisa horizontal se sujeta una masa de 2 Kg a una superficie vertical por
un medio de un resorte cuya constante es de 4 N/m. Otra masa de 1 kg se conecta al primer
objeto mediante un resorte con constante 2 N/m (ver figura 3). Sabiendo que los objetos se
desplazan 3 m hacia la derecha de sus posiciones de equilibrio y luego se sueltan, comprueba
que el sistema diferencial que rige el modelo es
d 2x
2 + 6 x − 2y = 0,
2
2dt
dy
+2 2=0
y− x ,
dt2
dx dy
x(0) = y(0) = 3, (0) = = 0.
dt dt
Calcula un sistema diferencial de primer orden equivalente al dado. Resuelve el problema,
representa la solución e interpreta los resultados.
28. Se considera el circuito eléctrico dado en la figura 4. Demuestra que el sistema diferencial
que rige el modelo es
1 t
I + 64I2jj = −2 sen ,
9 1 24
1 jj jj
64 I + 9I − 64I =
=3 I2 + I3,
I1 0,
3 2
donde I1.I2, I3 son las intensidades que circulan por las distintas ramas de la red.
i) Elimina I1, I3 del sistema para probar que se cumple
(92)(642)I(4) + (82)(64)Ijj + I2 = 0.
2 2
9
A. Manteniendo los restantes datos iguales, resuelve el problema, representa las soluciones e
interpreta los resultados.
31. Un famoso experimento de fı́sica llevado a cabo por Thomson para determinar la proporción
de carga respecto a la masa m para un electrón, usa las ecuaciones diferenciales
d2x + He dy
m = Ee,
dt2 2 dt
dx
dy
m 2 − He = 0,
dt dt
donde H es la intensidad del campo magnético y E es la intensidad del campo eléctrico que se
han aplicado.
i) Transforma el sistema anterior en un sistema equivalente de primer orden.
ii) Resuelve para el caso particular H = 4, e = 1, m = 2, E = 6.
iii) Resuelve usando transformada de Laplace, suponiendo que H = 1/4, e = 4, m = 1, E = 2
con condiciones iniciales x(0) = y(0) = 0, xj(0) = yj(0) = 1
10
reduce a la mitad suponiendo que la contaminación ha cesado, es decir, las afluencias que no
provienen de ninguno de los lagos son de agua limpia. Suponiendo que todos los lagos tienen la
misma concentración original de contaminación p, ¿cuánto tardará en reducirse la contaminación
al 5 por ciento del nivel original?.
donde c es una constante que depende del rozamiento y k es una constante que depende de la
elasticidad y tensión de la cuerda.
i) Resuelve (3) cuando c < k y f (x) = sen πx + sen 3πx cos πx.
ii) Si c = 1/2, k = 1 calcula T tal que
1
"u(x, t)"L2(0,1) ≤ "f (x)"L2(0,1), ∀ t ≥ T.
2
36. La transmisión de impulsos eléctricos en un cable largo con una distribución de capacitancia,
inductancia y resistencia, viene modelada por la ecuación de los telegraftstas:
11
∂ 2u ∂u 2
2∂ u
+2
a + bu − c = 0, a, b, c ∈ R, c > 0, 0 < x < π, t > 0.
∂t2 ∂t ∂x2
12
Calcula la solución formal del problema anterior cuando las condiciones iniciales y de contorno
son
u(0, t) = u(π, t) = 0, t ≥ 0,
u(x, 0) = 0, 0 ≤ x ≤ π,
∂u
(x, 0) = g(x), 0 ≤ x ≤ π.
∂t
13