Aplicaciones ED de Oren 1 MATE

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

APLICACIONES DE LAS ECUACIONES DIFERENCIALES DE PRIMER ORDEN

DINÁMICA DE UNA POBLACIÓN

Sea 𝑃 = 𝑃(𝑡) el número de habitantes en una población en el tiempo t, el modelo poblacional para calcular el
crecimiento exponencial comienza con el supuesto de que el crecimiento poblacional es proporcional al número
presente, es decir, 𝑃´(𝑡) = 𝐾𝑃, para cierta 𝐾 > 0. En este modelo, la tasa de crecimiento relativa o específica, es
definida por

𝑃´(𝑡)
𝑃
se asume como una constante 𝐾 > 0, si P aumenta, entonces la tasa relativa de crecimiento disminuye, o si P disminuye
la tasa de crecimiento aumenta.

El supuesto de que la tasa a la cual una población aumenta o disminuya depende solo del número presente y no de
mecanismos dependientes del tiempo tales como los fenómenos estacionales, por consiguiente, esta tasa de
crecimiento relativo lo podemos representar por

𝑃´(𝑡)
= 𝑓 (𝑃 ) o 𝑃´(𝑡) = 𝑃𝑓(𝑃) (1 ∗)
𝑃
Esta ecuación diferencial, por lo general se asume en los modelos de poblaciones de animales, se denomina hipótesis
de dependencia de la densidad poblacional.

EJEMPLO

Si la población de un país se duplica en 50 años ¿en cuántos años será el triple, suponiendo que la velocidad de aumento
es proporcional al número de habitantes?

Solución:

IDENTIFICANDO LAS VARIABLES

𝑡 el tiempo en años

𝑦 = 𝑦(𝑡) el número de habitantes en t años

𝑦0 = 𝑦(0) número inicial de habitantes en 𝑡 = 0 años

𝑦(50) = 2𝑦0

𝑦(𝑡) = 3𝑦0 ¿ 𝑡?

𝑦´ = 𝑦´(𝑡) velocidad de crecimiento de la población en 𝑡 años

𝑦´(𝑡) = 𝑘 𝑦(𝑡) crecimiento de la población es proporcional al número de habitantes, k constante de proporcionalidad.

1
𝑦´ = 𝑘 𝑦 ⟹ 𝑦´𝑑𝑡 = 𝑘𝑦𝑑𝑡 ⟹ 𝑑𝑦 = 𝑘𝑑𝑡 𝑉𝑆
𝑦

MODELO MATEMÁTICO

1
𝑑𝑦 = 𝑘𝑑𝑡
𝑦
𝑃𝑉𝐹: 𝑦(0) = 𝑦0
𝑦(50) = 2𝑦0
{𝑦(𝑡) = 3𝑦0 , ¿ 𝑡?

RESOLVIENDO EL MODELO MATEMÁTICO

1
𝑑𝑦 = 𝑘𝑑𝑡, ecuación en variables separables.
𝑦
Integrando

1
∫ 𝑑𝑦 = 𝑘 ∫ 𝑑𝑡
𝑦

𝑙𝑛|𝑦| = 𝐾𝑡 + 𝐶 ⟹ |𝑦| = 𝑒 𝐾𝑡+𝐶 = 𝑒 𝐶 𝑒 𝐾𝑡 ⟹ 𝑦 = ±𝑒 𝐶 𝑒 𝑘𝑡 ⟹ 𝑦 = 𝐶𝑒 𝐾𝑡 , ¿ 𝐶, 𝐾 constantes?

Hallando 𝐶

Si 𝑡 = 0, 𝑦 = 𝑦0 ⟹ 𝑦0 = 𝐶 ⟹ 𝑦 = 𝑦0 𝑒 𝐾𝑡

Ahora hallamos 𝐾

𝑙𝑛 2
Si 𝑡 = 50, 𝑦 = 2𝑦0 ⟹ 2𝑦0 = 𝑦0 𝑒 50𝐾 ⟹ 2 = 𝑒 50𝐾 ⟹ 50𝐾 = 𝑙𝑛 2 ⟹ 𝐾 =
50
𝑙𝑛 2 𝑡 𝑡
𝑦 = 𝑦0 𝑒 𝑡 50 = 𝑦0 𝑒 ln 250 ⟹ 𝑦(𝑡) = 𝑦0 250 ↝ población en 𝑡 años

Ahora hallamos el tiempo 𝑡 para que la población sea el triple


𝑡 𝑡 𝑡 ln 3
Si 𝑡 = 𝑡, 𝑦 = 3𝑦0 ⟹ 3𝑦0 = 𝑦0 250 ⟹ 3 = 250 ⟹ ln 3 = ln 2 ⟹ 𝑡 = 50 ≈ 79 años.
50 ln 2
Esto indica que en 79 años el número de habitantes de la población de un país será el triple de la inicial ∎

EJEMPLO

La reproducción de células cancerosas en un tejido obedece a un modelo de crecimiento exponencial. Suponga que un
tejido contiene inicialmente 10 células cancerosas y que después de una hora el número se ha duplicado.

A) Identifique las variables del problema para establecer el modelo matemático (ecuación diferencial con sus
condiciones) que represente el crecimiento de las células cancerosas en un tejido.
B) Resuelva el modelo matemático para determinar el número de células cancerosas en un tejido durante un día.
C) Interprete la solución.

Solución:

A) Identificando las variables, en este problema solamente hay dos variables, el tiempo como variable
independiente, las células cancerosas en un tejido como variable dependiente.
𝑡 el tiempo en horas.

𝑦 = 𝑦(𝑡) el número de células cancerosas en un tejido en 𝑡 horas.

𝑦(0) = 10 número inicial de células cancerosas en 𝑡 = 0 horas.

𝑦(1) = 20 células cancerosas

𝑦(24) =? ? ?

𝑦´ = 𝑦´(𝑡) velocidad de crecimiento de las células cancerosas en el tejido en 𝑡 horas

𝑦´ = 𝑘 𝑦 velocidad de crecimiento de las células cancerosas es proporcional al número de células cancerosas, k


constante de proporcionalidad.

1
𝑦´𝑑𝑡 = 𝑘 𝑦 𝑑𝑡 ⟹ 𝑑𝑦 = 𝑘𝑑𝑡 ecuación en variables separables
𝑦

1
𝑑𝑦 = 𝑘𝑑𝑡 representa el crecimiento de células cancerosas en el tejido
𝑦
MODELO MATEMÁTICO

1
𝑑𝑦 = 𝑘𝑑𝑡
𝑦
𝑃𝑉𝐹: 𝑦(0) = 10
𝑦(1) = 20
{𝑦(24) =? ? cuántas hay en 24 horas?

RESOLVIENDO EL MODELO MATEMÁTICO

Integrando

1
𝑑𝑦 = 𝑘𝑑𝑡, ecuación en variables separables
𝑦

1
∫ 𝑑𝑦 = 𝑘 ∫ 𝑑𝑡
𝑦

𝑙𝑛|𝑦| = 𝑘𝑡 + 𝑐 ⟹ |𝑦| = 𝑒 𝑘𝑡+𝑐 = 𝑒 𝑐 𝑒 𝑘𝑡 ⟹ 𝑦 = ±𝑒 𝑐 𝑒 𝑘𝑡 ⟹ 𝑦 = 𝐶𝑒 𝑘𝑡 , ¿ 𝐶, 𝑘 constantes?

Hallando C

Si 𝑡 = 0, 𝑦 = 10 ⟹ 10 = 𝐶 ⟹ 𝑦 = 10𝑒 𝑘𝑡

Ahora hallamos 𝐾
𝑡
Si 𝑡 = 1, 𝑦 = 20 ⟹ 20 = 10𝑒 𝑘 ⟹ 2 = 𝑒 𝑘 ⟹ 𝑘 = 𝑙𝑛 2 ⟹ 𝑦 = 10𝑒 𝑡 𝑙𝑛2 = 10𝑒 𝑙𝑛2 = 10 2𝑡

𝑦(𝑡) = 10 2𝑡 esto indica el número de células cancerosas en el tejido en t horas.

Ahora hallamos el número de células cancerosas en un día

𝑡 = 24 ⟹ 𝑦(24) = 10 224 = 10(16 777 216) = 167 772 160 ⟹ 𝑦(24) = 167 772 160

INTERPRETACIÓN DE LA SOLUCIÓN

Esto indica que, en un día, es decir que en 24 horas hay 167 772 160 células cancerosas en el tejido ∎

ECUACIÓN LOGÍSTICA

Suponga que un entorno es capaz de dar sustento solamente a un número fijo de


individuos K. K representa la cantidad de soportes del entorno (ambiente). Por
tanto, la función 𝑓(𝑃) dada en la ecuación anterior

𝑃´(𝑡)
= 𝑓 (𝑃 ) o 𝑃´(𝑡) = 𝑃𝑓(𝑃) (1 ∗)
𝑃
tenemos 𝑓(𝐾 ) = 0 y asumimos que 𝑓 (0) = 𝑟 habitantes. La figura muestra tres
funciones que satisfacen estas dos condiciones. El supuesto más simple que
podemos hacer es que 𝑓(𝑃) sea una línea recta, es decir,

𝑓´(𝑃) = 𝑐1
𝑃𝑉𝐹: { 𝑓 (0) = 𝑟
𝑓 (𝐾 ) = 0

𝑓(𝑃) = 𝑐1 𝑃 + 𝑐2

Usando las condiciones


𝑟
𝑓(0) = 𝑟 ⟹ 𝑐2 = 𝑟, 𝑓 (𝐾 ) = 𝑐1 𝐾 + 𝑟 = 0 ⟹ 𝑐1 = −
𝐾
𝑟
Por consiguiente, la función 𝑓 (𝑃) toma la forma de 𝑓(𝑃) = − 𝐾 𝑃 + 𝑟 y la ecuación (1 ∗), se convierte en

𝑟
𝑃´(𝑡) = 𝑃𝑓 (𝑃) = 𝑃 (𝑟 − 𝑃)
𝐾
𝑟
𝑃´(𝑡) = 𝑃 (𝑟 − 𝑃) (2 ∗)
𝐾
La ecuación (2 ∗) la podemos escribir nuevamente por

𝑃´(𝑡) = 𝑃(𝑎 − 𝑏𝑃) (3 ∗ )


𝑟
donde 𝑎 y 𝑏 son constantes positivas, 𝑎 = 𝑟, 𝑏 = 𝐾

Alrededor de 1840, el matemático y biólogo belga P. F. Verhulst se interesó en los modelos matemáticos para predecir
las poblaciones humanas de diferentes países, estudió la ecuación (3 ∗) que se conoce como la ecuación logística y su
solución se denomina función logística y la gráfica de una ecuación logística se denomina curva logística.

La ecuación diferencial 𝑃´(𝑡) = 𝐾𝑃 no ofrece un modelo muy preciso de la población cuando la población es muy
grande. Las condiciones de sobrepoblación, con sus resultantes efectos perjudiciales para el entorno, tales como la
contaminación y las demandas excesivas y competitivas por el alimento y el combustible, pueden tener un efecto
inhibidor sobre el crecimiento poblacional.

La ecuación (3 ∗) presenta un acotamiento en su solución cuando 𝑡 ⟶ +∞ , esta ecuación la podemos escribir por
𝑃´(𝑡) = 𝑃𝑎 − 𝑏𝑃2 , al término no lineal 𝑏𝑃2 , puede interpretarse como un término de inhibición o competencia. En la
mayoría de las aplicaciones, la constante 𝑎 es mucho mayor que la constante 𝑏.

Las curvas logísticas han demostrado ser mas muy precisas para pronosticar los patrones de crecimiento de ciertos tipos
de bacterias, protozoarios, pulgas acuáticas y moscos de las frutas, pero en un espacio limitado.

SOLUCIÓN DE LA ECUACIÓN LOGÍSTICA

Para la solución, podemos emplear el método de separación de variables e integración de fracciones parciales, veamos

1 1 1
𝑃´ = 𝑃𝑎 − 𝑏𝑃2 ⟹ 𝑃´ = 1 ⟹ 𝑃´𝑑𝑡 = 𝑑𝑡 ⟹ 𝑑𝑃 = 𝑑𝑡
𝑃𝑎 − 𝑏𝑃2 𝑃𝑎 − 𝑏𝑃2 𝑃𝑎 − 𝑏𝑃2
1/𝑎 𝑏/𝑎 1 1 𝑃 𝑃
( + ) 𝑑𝑃 = 𝑑𝑡 ⟹ ln|𝑃| − ln|𝑎 − 𝑏𝑃| = 𝑡 + 𝐶 ⟹ ln | | = 𝑎𝑡 + 𝑎𝐶 ⟹ | | = 𝑒 𝑎𝑡+𝑎𝐶
𝑃 𝑎 − 𝑏𝑃 𝑎 𝑎 𝑎 − 𝑏𝑃 𝑎 − 𝑏𝑃
𝑃 𝑃
= ± 𝑒 𝑎𝐶 𝑒 𝑎𝑡 = 𝑐𝑒 𝑎𝑡 ⟹ = 𝑐𝑒 𝑎𝑡 ⟹ 𝑃 = (𝑎 − 𝑏𝑃)𝑐𝑒 𝑎𝑡 ⟹ 𝑃 = 𝑎𝑐𝑒 𝑎𝑡 − 𝑏𝑐𝑃𝑒 𝑎𝑡
𝑎 − 𝑏𝑃 𝑎 − 𝑏𝑃
𝑎𝑐𝑒 𝑎𝑡 𝑎𝑐 𝑎𝑐
𝑃 (𝑡 ) = 𝑎𝑡
= −𝑎𝑡 ⟹ 𝑃(𝑡) = −𝑎𝑡
1 + 𝑏𝑐𝑒 𝑒 + 𝑏𝑐 𝑒 + 𝑏𝑐
𝑎 𝑃0
Si 𝑃(0) = 𝑃0, 𝑃0 ≠ ⟹ 𝑐=
𝑏 𝑎 − 𝑏𝑃0

𝑎𝑃0
𝑃 (𝑡 ) = solución de la ecuación logística (4 ∗)
(𝑎 − 𝑏𝑃0 )𝑒 −𝑎𝑡 + 𝑏𝑃0
GRÁFICA DE LA SOLUCIÓN DE LA ECUACIÓN LOGÍSTICA

La variable 𝑡 representa el tiempo y rara vez nos interesamos cuando


𝑡 < 0, aunque es muy importante incluir este intervalo cuando
presentamos las diferentes gráficas de la curva logística 𝑃

De (4 ∗) podemos deducir

𝑎𝑃0
𝑃 (𝑡 ) = solución de la ecuación logística (4 ∗)
(𝑎 − 𝑏𝑃0 )𝑒 −𝑎𝑡 + 𝑏𝑃0
𝑎 𝑎
lim 𝑃(𝑡) = ⟹ 𝑃(𝑡) = ↝ asíntota horizontal
𝑡⟶+∞ 𝑏 𝑏
lim 𝑃(𝑡) = 0 ⟹ 𝑃 (𝑡) = 0 ↝ asíntota horizontal (eje t)
𝑡⟶−∞

Veamos los puntos de inflexión de esta gráfica, para esto derivemos


simplemente la ecuación 𝑃´(𝑡) = 𝑃𝑎 − 𝑏𝑃2 y apliquemos la regla de la
cadena

𝑃´´(𝑡) = 𝑃´𝑎 − 2𝑏𝑃𝑃´ = 𝑃´(𝑎 − 2𝑏𝑃), donde 𝑃´ = 𝑃𝑎 − 𝑏𝑃2

𝑃´´(𝑡) = (𝑃𝑎 − 𝑏𝑃2 )(𝑎 − 2𝑏𝑃) = 𝑃(𝑎 − 𝑏𝑃 )(𝑎 − 2𝑏𝑃)

𝑃´´(𝑡) = 𝑃(𝑎 − 𝑏𝑃 )(𝑎 − 2𝑏𝑃)

𝑃´´(𝑡) = 𝑃(𝑏𝑃 − 𝑎 )(2𝑏𝑃 − 𝑎)


𝑎 𝑎
𝑃´´(𝑡) = 2𝑏2 𝑃 (𝑃 − ) (𝑃 − )
𝑏 2𝑏
𝑎
Para halar los puntos de inflexión, hacemos 𝑃´´(𝑡) = 0, descartamos los puntos críticos, 𝑃 = 0, 𝑃 = 𝑏, porque son
𝑎
asíntotas horizontales, entonces el único punto crítico de segunda especie es el punto de inflexión 𝑃 = 2𝑏 .

Analicemos la concavidad de la curva logística, en el eje P(t).

Valor de 𝑃 𝑎 𝑎 𝑎 𝑎 𝑎 𝑎
0 0<𝑃< <𝑃<
2𝑏 2𝑏 2𝑏 2𝑏 𝑏 𝑏
𝑎 𝑃´´(𝑡) ↝ + 𝑃´´(𝑡) ↝ −
𝑃=
4𝑏
CONCLUSIÓN La curva es cóncava hacia arriba La curva es cóncava hacia abajo
𝑎
Observamos que (𝑡, 2𝑏) es un punto de inflexión para la curva logística

𝑎 𝑎
Ahora analicemos el valor inicial 𝑃0 que puede estar debajo de o encima de . En nuestro cuadro veamos
2𝑏 2𝑏

Valor de 𝑃 𝑎 𝑎 𝑎 𝑎 𝑎 𝑎
0 0 < 𝑃0 < 𝑃 < < 𝑃 < 𝑃0 <
2𝑏 2𝑏 2𝑏 2𝑏 𝑏 𝑏
𝑎 𝑃´´(𝑡) ↝ + 𝑃´´(𝑡) ↝ −
𝑃=
4𝑏
CONCLUSIÓN La curva es cóncava hacia arriba, La curva es cóncava hacia abajo, la
la curva toma la forma de una S curva toma la forma de una S y el
y el punto de inflexión está en el punto de inflexión está en el
primer cuadrante 𝑡 > 0 segundo cuadrante 𝑡 < 0
𝑎
Observamos que (𝑡, 2𝑏) es un punto de inflexión para la curva logística
PROPAGACIÓN DE UNA ENFERMEDAD CONTAGIOSA Y APLICACIÓN DE LA ECUACIÓN LOGÍSTICA

Una enfermedad contagiosa se difunde en una comunidad por medio del contacto físico entre las personas.

Si 𝑥 = 𝑥(𝑡) indica el número de personas que han contraído la enfermedad (personas enfermas) en el tiempo 𝑡 y
𝑦 = 𝑦(𝑡) indica el número de personas que no han contraído la enfermedad (personas sanas) en el tiempo 𝑡,

𝑥´ = 𝑥´(𝑡) indica la razón a la que se propaga la enfermedad.

Suponga que la propagación de la enfermedad es proporcional al número de encuentros o interacciones ente estos dos
grupos de gentes, es decir 𝑥´ = 𝑘 𝑥 (𝑡) 𝑦(𝑡), además suponga que la comunidad es pequeña con 𝑛 personas y que ha
llegado a la comunidad una persona infectada, 𝑥(𝑡) y 𝑦(𝑡) se encuentran relacionados por 𝑥 + 𝑦 = 𝑛 + 1, entonces
𝑦 = 𝑛 + 1 − 𝑥, esto representa el número de personas no infectadas. Entonces el modelo del problema de valor inicial
para este problema es

𝑥´ = 𝑘 𝑥 (𝑛 + 1 − 𝑥 )
{
𝑥 (0) = 1

DRENADO DE TANQUES

EJEMPLO

Un tanque para agua tiene la forma de un cilindro recto, está lleno de agua, tiene un radio de 8 cm y una altura de 10
1
cm. Si en el fondo del tanque tiene un orificio circular de radio 12
cm, el orificio está tapado con un tapón. Cuando se
cm
saca el tapón el agua fluye por el orificio con una velocidad de 𝑣 = 4,8√ℎ , donde ℎ es el nivel del agua en el
seg
tanque.

A) Identifique las variables del problema para establecer el modelo matemático (ecuación diferencial con sus
condiciones) para determinar el cambio de nivel del agua en el tanque, usando la Ley de Torricelli.
B) Resuelva el modelo matemático determinar el tiempo en que el tanque estará vacío.
C) Interprete la solución.

Solución:

A) Identificando las variables, en este problema solamente hay


dos variables, el tiempo como variable independiente y el nivel
del agua en el tanque como variable dependiente.

Sea 𝑡 el tiempo en segundos.

ℎ = ℎ(𝑡) el nivel del agua en el tanque en el tiempo 𝑡.

∆ℎ altura del cilindro de pérdida de agua en el tanque.

∆𝑠 altura del cilindro de agua de salida por el orificio.

𝑅 = 8 radio del volumen del cilindro.

𝑟 radio del orificio por donde fluye el agua.

𝑉𝑃 volumen de pérdida de agua en el tanque,

𝑉𝑃 = 𝜋 82 ∆ℎ = 64𝜋∆ℎ.
1
𝑉𝑆 volumen de salida de agua en el tanque, cilindro circular de radio 𝑟 = 12 cm y altura ∆𝑠 (espesor del tanque).

𝜋 𝜋
𝑉𝑆 = 2
∆𝑠 = ∆𝑠
12 144
Por la Ley de Torricelli: 𝑉𝑃 = − 𝑉𝑆

𝜋 1
64𝜋∆ℎ = − ∆𝑠 ⟹ 64∆ℎ = − ∆𝑠
144 144
Dividiendo por ∆𝑡 y tomando el límite ∆𝑡 ⟶ 0 tenemos

∆ℎ 1 ∆𝑠 1 1 1
64 =− ⟹ 64 ℎ´(𝑡) = − 𝑠´(𝑡) = − 𝑣 ⟹ ℎ´(𝑡) = − 𝑣
∆𝑡 144 ∆𝑡 144 144 64(144)

Donde 𝑠´(𝑡) = 𝑣 = 4,8√ℎ es la velocidad por donde fluye el agua por el orificio, entonces

1 1 1 1 1
ℎ´(𝑡) = − 4,8√ℎ ⟹ ℎ´(𝑡)𝑑𝑡 = − 𝑑𝑡 ⟹ ℎ− 2 𝑑ℎ = − 𝑑𝑡
64(144) √ℎ 1 920 1 920
1 1
ℎ− 2 𝑑ℎ = − 𝑑𝑡 ecuación en variables seapradas
1 920
1 1
ℎ− 2 𝑑ℎ = − 𝑑𝑡 representa el cambio de nivel de agua en el tanque
1 920
Modelo Matemático
1 1
ℎ−2 𝑑ℎ = − 𝑑𝑡
1 920
𝑃𝑉𝐼: {
ℎ(0) = 10
ℎ(𝑡) = 0 ¿ 𝑡?

B) Resolviendo el Modelo Matemático


Integrando
1 1
∫ ℎ−2 𝑑ℎ = − ∫ 𝑑𝑡
1 920
1 1 𝐶 1
2√ℎ = − 𝑡 + 𝐶 ⟹ √ℎ = − 𝑡+ ⟹ √ℎ = − 𝑡+𝐶
1 920 1 920(2) 2 3 840

Hallando C

Si 𝑡 = 0, ℎ = 10 ⟹ √10 = 0 + 𝐶 ⟹ 𝐶 = √10

1
√ℎ = − 𝑡 + √10 ↝ nivel del agua en el tanque
3 840
Hallando el tiempo 𝑡 para vaciar el tanque

1
Si 𝑡 = 𝑡, ℎ = 0 ⟹ − 𝑡 + √10 = 0 ⟹ 𝑡 = 3 840 √10 𝑠𝑒𝑔.
3 840
Convirtiendo el tiempo en horas

𝑡 = 3 840√10 = 12 143. 14622 seg = 12 143. 14622 ÷ 3 600 = 3.373096171 horas

𝑡 = 3 horas + 0.373096171 horas = 3 horas + 0.373096171(60) minutos

𝑡 = 3 horas + 0.373096171
⏟ (60) minutos
22 𝑚𝑖𝑛𝑢𝑡𝑜𝑠

INTRERPRETACIÓN DE LA SOLUCIÓN

Para vaciar el tanque se demora aproximadamente 3 horas con 22 minutos ∎


LEY DE NEWTON DE ENFRIAMIENTO Y CALENTAMIENTO

Ley de enfriamiento de Newton, la velocidad a la que se enfría una sustancia al aire libre es proporcional a la diferencia
entre la temperatura de la sustancia y la del aire. Es decir

𝑇´(𝑡) = −𝑘(𝑇(𝑡) − 𝑇𝑎 )

donde 𝑇 = 𝑇(𝑡) es la temperatura de la sustancia en el tiempo 𝑡, 𝑇𝑎 es la temperatura del medio y 𝑘 es la constante


de proporcionalidad, 𝑇´ = 𝑇´(𝑡) velocidad de enfriamiento de la sustancia en el tiempo t.

EJEMPLO

Si la temperatura del medio es de 300 y la sustancia se enfría de 1000 a 700 en 15 minutos ¿en qué tiempo la
temperatura de la sustancia será de 400 ?

Solución:

IDENTIFICANDO LAS VARIABLES

𝑡 tiempo en minutos

𝑇 = 𝑇(𝑡) temperatura de la sustancia en el tiempo 𝑡

𝑇𝑎 = 30 temperatura del medio

𝑇(0) = 100 temperatura inicial

𝑇(15) = 70

𝑇(𝑡) = 40 ¿ 𝑡?

𝑇´ = 𝑇´(𝑡) velocidad de enfriamiento de la sustancia en el tiempo t

De acuerdo a la ley de Newton de enfriamiento

𝑇´(𝑡) = −𝑘(𝑇(𝑡) − 𝑇𝑎 )

1 1
𝑇´𝑑𝑡 = −𝑘𝑑𝑡 ⟺ 𝑑𝑇 = −𝑘𝑑𝑡 𝑉𝑆
𝑇 − 30 𝑇 − 30
MODELO MATEMÁTICO

1
𝑑𝑇 = −𝑘𝑑𝑡
𝑇 − 30
𝑃𝑉𝐹: 𝑇(0) = 100
𝑇(15) = 70
{ 𝑇(𝑡) = 40 ¿ 𝑡?

RESOLVIENDO EL MODELO MATEMÁTICO

Integrando

1
∫ 𝑑𝑇 = −𝐾 ∫ 𝑑𝑡
𝑇 − 30

ln|𝑇 − 30| = −𝐾𝑡 + 𝐶 ⟹ |𝑇 − 30| = 𝑒 −𝑘𝑡+𝑐 = 𝑒 𝑐 𝑒 −𝑘𝑡 ⟹ 𝑇 − 30 = ±𝑒 𝑐 𝑒 −𝑘𝑡 = 𝐶𝑒 −𝑘𝑡

𝑇 − 30 = 𝐶𝑒 −𝑘𝑡 temperatura de la sustancia en el tiempo t minutos.

Hallando C

Usando las condiciones iniciales hallamos las constantes C y K

Si 𝑡 = 0 , 𝑇 = 100 ⟹ 100 − 30 = 𝐶 ⟹ 𝐶 = 70
𝑇 − 30 = 70 𝑒 −𝑘𝑡

Hallando K

Si 𝑡 = 15 , 𝑇 = 70 ⟹ 70 − 30 = 70 𝑒 −15 𝑘 ⟹ 40 = 70 𝑒 −15 𝑘 ⟹ 4 = 7 𝑒 −15 𝑘

4 4 1 4 𝑡
ln
4
= 𝑒 −15 𝑘 ⟹ ln = −15𝐾 ⟹ 𝐾=− ln ⟹ 𝑇 − 30 = 70 𝑒 15 7
7 7 15 7
𝑡 𝑡 𝑡
4 15 4 15 4 15
𝑇 − 30 = 70 𝑒 ln(7) = 70 ( ) ⟹ 𝑇 − 30 = 70 ( )
7 7
𝑡
4 15
𝑇 − 30 = 70 ( ) temperatura de la sustancia en el tiempo t minutos.
7
Hallando el tiempo t cuando 𝑇 = 40, entonces
𝑡 𝑡 𝑡 𝑡
4 15 4 15 4 15 1 4 15
40 − 30 = 70 ( ) ⟹ 10 = 70 ( ) ⟹ 1 = 7( ) ⟹ =( )
7 7 7 7 7
1 1
1 𝑡 4 ln ln
ln = ln ⟹ 7= 𝑡 ⟹ ⟹ 𝑡 = 15 7 = 15
ln 7
= 15 (3.47722525169)
7 15 7 4 15 4 7
ln ln ln
7 7 4
𝑡 = 52.1583787754 minutos.

𝑡 ≈ 52 minutos, indica que en 52 minutos la sustancia tendrá una temperatura de 400 ∎

CAÍDA LIBRE

EJEMPLO: Un paracaidista está cayendo con una velocidad de 176 pies/seg = 53.65 m/seg cuando se abre su
paracaídas. Si la resistencia del aire es Wv 2 /256 lib, donde W es el peso total del hombre y del paracaídas. Establecer
la ecuación diferencial con sus condiciones y luego solucione para hallar la velocidad 𝑣 = 𝑣 (𝑡) del paracaidista en
función del tiempo t después de abrirse el paracaídas.

Solución:

IDENTIFICANDO LAS VARIABLES

𝑡 tiempo en segundos

𝑣 (0) = 176 pies/seg velocidad inicial con la desciende el paracaidista

𝑣 = 𝑣 (𝑡) del paracaidista en 𝑡 segundos

𝐹𝑊 = 𝑊 es el peso del hombre y del paracaídas, dirigida hacia abajo +

𝐹𝑎 = −Wv 2 /256 resistencia del aire, dirigida hacia arriba −

En este sistema actúan dos fuerzas: 𝐹𝑊 es el peso del hombre y del paracaídas
dirigida hacia abajo + y la resistencia del aire 𝐹𝑎 dirigida hacia arriba −

Por la segunda ley de Newton

∑ 𝐹 = 𝑚𝑎 ⟹ 𝑚𝑎 = peso del hombre + resistencia del aire

Wv 2 𝑊 𝑝𝑖𝑒
𝑚𝑎 = 𝑊 − , pero 𝑎 = 𝑣´ (𝑡), donde 𝑚 = , 𝑔 = 32
256 𝑔 𝑠𝑒𝑔 2
𝑊 Wv 2 𝑊 Wv 2 1 v2 256 − 𝑣 2 256 − 𝑣 2
𝑣´ = 𝑊 − ⟹ 𝑣´ = 𝑊 − ⟹ 𝑣´ = 1 − = ⟹ 𝑣´ =
𝑔 256 32 256 32 256 256 8

1 1 1 1 1 1
𝑣´ = ⟹ 𝑣´𝑑𝑡 = 𝑑𝑡 ⟹ 𝑑𝑣 = 𝑑𝑡 ⟹
256 − 𝑣 2 8 256 − 𝑣 2 8 256 − 𝑣 2 8
1 1
𝑑𝑣 = − 𝑑𝑡 cambio de velocidad con la que descendiende el paracaidista
𝑣 2 − (16)2 8

MODELO MATEMÁTICO

1 1
𝑑𝑣 = − 𝑑𝑡
𝑣 2 − (16)2 8
𝑃𝑉𝐼:
𝑣 (0) = 176
{ 𝑣 = 𝑣 (𝑡 ) ?

RESOLUCIÓN DEL MODELO MATEMÁTICO

Integrando en 𝑡

1 1 1 𝑣 − 16 1
∫ 𝑑𝑣 = ∫ 𝑑𝑡 ⟹ ln | | =− 𝑡+𝐶
𝑣 2 − (16)2 8 32 𝑣 + 16 8
𝑣 − 16
ln | | = −4𝑡 + 𝐶
𝑣 + 16
Para hallar C usamos la condición, para 𝑡 = 0 , 𝑣 = 176

176 − 16 160 5
ln | | = 𝐶 ⟹ ln | | = 𝐶 ⟹ 𝐶 = ln
176 + 16 192 6
𝑣 − 16 5 𝑣 − 16 5 5 𝑣 − 16 5 −4 𝑡
ln | | = −4𝑡 + ln ⟹ | | = 𝑒 −4𝑡+ln6 = 𝑒 −4 𝑡 ⟹ = 𝑒
𝑣 + 16 6 𝑣 + 16 6 𝑣 + 16 6

6𝑣 − 16(6) = (5𝑣 + 16(5)) 𝑒 −4 𝑡 ⟹ 6𝑣 − 5𝑣𝑒 −4 𝑡 = 16(5)𝑒 −4 𝑡 + 16(6)

6 + 5𝑒 −4 𝑡
(6 − 5𝑒 −4𝑡 )𝑣 = 16(6 + 5𝑒 −4𝑡 ) ⟹ 𝑣 (𝑡) = 16 velocidad del paracaídista en el tiempo 𝑡.
6 − 5𝑒 −4 𝑡
6 + 5𝑒 −4 𝑡
𝑣 (𝑡) = 16 ∎
6 − 5𝑒 −4 𝑡
EJEMPLO

Un cuerpo de masa 𝑚 unidades técnicas partiendo del reposo cae en un medio donde la resistencia (lb) es proporcional
al cuadrado de la velocidad (pies/seg). Si la velocidad terminal es de 150 pies/seg = 45,72 m/seg. A) Halle la
velocidad del cuerpo al cabo de 2 seg. B) Halle el tiempo necesario para que la velocidad sea de
100 pies/seg = 30,48 m/seg.

Solución:

DENTIFICANDO LAS VARIABLES

𝑡 tiempo en segundos

𝑣 (0) = 0 velocidad inicial, porque el cuerpo parte del reposo

lim 𝑣 = 150 velocidad terminal


𝑡⟶+∞

𝑣 = 𝑣(𝑡) velocidad del cuerpo en 𝑡 segundos

𝐹𝐶 = 𝑚𝑔 es peso del cuerpo con la que está cayendo, dirigida hacia abajo +

𝑅𝑚 = −𝑘 𝑣 2 resistencia del medio con la que el cuerpo está cayendo, dirigida hacia arriba −
En este sistema actúan dos fuerzas: peso del cuerpo 𝐹𝐶 = 𝑚𝑔 (dirigida hacia abajo +) y la resistencia del medio (dirigida
hacia arriba −) 𝑅𝑚 = −𝑘𝑣 2

Por la segunda ley de Newton

𝑝𝑖𝑒
∑ 𝐹 = 𝑚𝑎 ⟹ 𝑚𝑎 = peso del cuerpo + resistencia del medio, donde , 𝑔 = 32 , 𝑎 = 𝑣´(𝑡)
𝑠𝑒𝑔2
𝐾 2 𝐾
𝑚𝑎 = 𝑚𝑔 − 𝐾𝑣 2 ⟺ 𝑎 = 𝑔 − 𝑣 ⟺ 𝑣´ = 32 − 𝑣 2
𝑚 𝑚
1 1 1
𝑣´ = 1 ⟹ 𝑣´𝑑𝑡 = 𝑑𝑡 ⟹ 𝑑𝑣 = 𝑑𝑡
𝐾 𝐾 2 𝐾
32 − 𝑚 𝑣 2 32 − 𝑚 𝑣 32 − 𝑚 𝑣 2

1 1
𝑑𝑣 = − 𝑑𝑡 cambio de velocidad con la que descendiende el cuerpo
𝐾𝑣 2 − 32𝑚 𝑚
MODELO MATEMÁTICO

1 1
𝑑𝑣 = − 𝑑𝑡
𝐾𝑣 2 − 32𝑚 𝑚
𝑣 (0) = 0
𝑃𝑉𝐼: lim 𝑣 = 150
𝑡⟶+∞
𝑣 (2) = ¿ ?
{ 𝑣 𝑡) = 100 ¿ 𝑡?
(

1 1 1 𝑘 1 𝑘
∫ 𝑑𝑣 = − ∫ 𝑑𝑡 ⟹ ∫ 𝑚 𝑑𝑣 = − 𝑚 ∫ 𝑑𝑡 ⟹ ∫ 2 𝑑𝑣 = − ∫ 𝑑𝑡
𝐾𝑣 2 − 32𝑚 𝑚 2
𝑣 − 32 𝑚 𝑚
𝐾 2
𝑣 − (√32 𝐾 )

𝑚
𝑣 − √32 𝐾
√𝐾 𝐾
ln || || = − 𝑡 + 𝐶, para simplificar los cálculos en el resultado hagamos 𝐾 = 2𝑚𝑘 2
2√32 𝑚 𝑚 𝑚
𝑣 + √32 𝐾

𝑚 16
𝑣 − √32 2𝑚𝑘 2 𝑣−√ 2
√2𝑚 𝑘 2𝑚𝑘 2 √2𝑚 𝑘 𝑘 |
ln || || = − 𝑡+𝐶 ⟹ ln || | = −2𝑘 2 𝑡 + 𝐶
2√32 𝑚 𝑚 𝑚 8√2𝑚 16
𝑣 + √32 𝑣+√ 2
2𝑚𝑘 2 𝑘
4
𝑘 𝑣−𝑘 𝑘 𝑘𝑣 − 4 𝑘𝑣 − 4 8𝐶
ln | | = −2𝑘 2 𝑡 + 𝐶 ⟹ ln | | = −2𝑘 2 𝑡 + 𝐶 ⟹ ln | | = −16𝑘𝑡 + = −16𝑘𝑡 + 𝐶
8 4 8 𝑘𝑣 + 4 𝑘𝑣 + 4 𝑘
𝑣+𝑘

𝑘𝑣 − 4 𝑘𝑣 − 4 𝑘𝑣 − 4
ln | | = −16𝑘𝑡 + 𝐶 ⟹ | | = 𝑒 −16𝑘𝑡+𝐶 ⟹ = ±𝑒 −16𝑘𝑡+𝐶 = ±𝑒 𝐶 𝑒 −16𝑘𝑡 = 𝐶𝑒 −16𝑘𝑡
𝑘𝑣 + 4 𝑘𝑣 + 4 𝑘𝑣 + 4
𝑘𝑣 − 4
= 𝐶𝑒 −16𝑘𝑡
𝑘𝑣 + 4
Hallando C

Para hallar la constante C usamos la condición inicial 𝑡 = 0, 𝑣 = 0 ⟹ 𝐶 = −1

𝑘𝑣 − 4
= − 𝑒 −16𝑘𝑡
𝑘𝑣 + 4
Para hallar la constante K usamos la condición de la velocidad terminal es de 150 𝑝𝑖𝑒𝑠/𝑠𝑒𝑔

𝑘𝑣 − 4 𝑘150 − 4 4 2 2
lim = − lim 𝑒 −16𝑘𝑡 = 0 ⟹ = 0 ⟹ 𝑘150 − 4 = 0 ⟹ 𝑘 = = ⟹ 𝑘=
𝑡⟶+∞ 𝑘𝑣 + 4 𝑡⟶+∞ 𝑘150 + 4 150 75 75
2 1
75 𝑣 − 4 = −𝑒 −16 75
2 𝑣−2 2 𝑣 − 150 32 32
𝑡
⟹ 75 = −𝑒 −16 75 𝑡 ⟹ = − 𝑒 −75 𝑡 , donde = 0.4266 ≈ 0.43
2 1 𝑣 + 150 75
75 𝑣 + 4 75 𝑣 + 2
𝑣 − 150
= − 𝑒 −0.43 𝑡 ⟹ 𝑣 − 150 = −𝑣 𝑒 −0.43 𝑡 −150 𝑒 −0.43 𝑡 ⟹ (1 + 𝑒 −0.43 𝑡 )𝑣 = 150(1 − 𝑒 −0.43 𝑡 )
𝑣 + 150
150(1 − 𝑒 −0.43 𝑡 )
𝑣= ↝ velocidad del cuerpo en el tiempo 𝑡 𝑠𝑒𝑔.
1 + 𝑒 −0.43 𝑡
A) Para 𝑡 = 2 𝑠𝑒𝑔

150(1 − 𝑒 −0.43 (2) ) 150(1 − 𝑒 −0.86 ) 𝑝𝑖𝑒 𝑝𝑖𝑒


𝑣= −0.43 (2)
= −0.86
≈ 61 ⟹ 𝑣 (2) ≈ 61
1+ 𝑒 1+ 𝑒 𝑠𝑒𝑔 𝑠𝑒𝑔

INTERPRETACIÓN DE LA SOLUCIÓN

𝑝𝑖𝑒
En 2 segundos la velocidad del cuerpo es de 61 ∎
𝑠𝑒𝑔

B) Si 𝑣 = 100 ¿ 𝑡?

𝑣 − 150 100 − 150 50 1


= − 𝑒 −0.43 𝑡 ⟹ = − 𝑒 −0.43 𝑡 ⟹ − = −𝑒 −0.43 𝑡 ⟹ = 𝑒 −0.43 𝑡
𝑣 + 150 100 + 150 250 5
1 1
ln 5 = −0.43 𝑡 ⟹ 𝑡 = 0.43 ln 5 ≈ 3.7 𝑠𝑒𝑔 ⟹ 𝑡 ≈ 3.7 𝑠𝑒𝑔

INTERPRETACIÓN DE LA SOLUCIÓN
𝑝𝑖𝑒𝑠
En 𝑡 ≈ 3.7 𝑠𝑒𝑔 el cuerpo tiene una velocidad de 100 𝑠𝑒𝑔 ∎

CIRCUITOS ELÉCTRICOS

La ley que gobierna en los circuitos eléctricos es la ley de Kirchhoff, la cual describe la distribución de la corriente 𝐼.

El circuito eléctrico más simple es el circuito en serie, el cual está compuesto de una fuerza electromotriz (fem) que
actúa como una fuente de energía tal como una batería o generador, y una resistencia que usa energía como por
ejemplo un foco, plancha o cualquier artefacto eléctrico.

La ley de Ohm

En un circuito simple compuesto por una fuerza electromotriz E y una resistencia R: La fuerza electromotriz E es
directamente proporcional a la corriente instantánea 𝐼, es decir E = RI, donde R es la resistencia.

Leyes Experimentales

1. La caída de voltaje ER a través de una resistencia R es proporcional a la


corriente 𝐼 que pasa a través de la resistencia, es decir 𝐸𝑅 = 𝑅𝐼.
2. La caída de voltaje EL a través de un inductor es proporcional a la tasa de
tiempo instantánea de cambio de la corriente I , es decir 𝐸𝐿 = 𝐿𝐼´(𝑡), donde 𝐿 es
la inductancia o coeficiente de inductancia.
3. La caída de voltaje EC a través de un condensador es proporcional a la carga
1 1
eléctrica instantánea Q en el condensador, es decir EC = 𝐶 Q, donde es la
𝐶
constante de proporcionalidad, C es el coeficiente de capacitancia o simplemente
capacitancia.
La ley de Kirchhoff

La suma algebraica de todas las caídas de voltaje alrededor de un circuito eléctrico es cero.

Otra manera de enunciar: el voltaje suministrado por la fuerza electromotriz fem es igual a la suma de todas las caídas
de voltaje.

Simbología usada en los circuitos

Cantidad Símbolo Unidad


Voltaje, fem o potencial EoV Voltio
Resistencia R Ohmio
Inductancia L Henrio
Capacitancia C Faradio
Corriente I Amperio
Carga Q Culombio

Convención: la corriente fluye del lado positivo (+) de la batería a través del circuito hacia el lado negativo (–)

Por la ley de Kirchhoff, en la figura 1, la fem suministrada E es igual a la caída de voltaje a través del inductor 𝐿 es
LI´(t) más la caída de voltaje a través de la resistencia 𝑅 es RI, tenemos la ecuación diferencial para el circuito 1

𝑑𝐼
𝐿 + 𝑅𝐼 = 𝐸 (𝑡) ⟺ 𝐿 𝐼´ + 𝑅 𝐼 = 𝐸 (𝑡)
𝑑𝑡
donde 𝐿 es la inductancia (henrios), 𝑅 es la resistencia (ohmios), 𝐼 es la corriente (amperios) y 𝐸 es la fuerza
electromotriz f.e.m. (voltios), para nuestro caso 𝐿 y 𝑅 constantes, esta es una ecuación lineal.

Para el circuito en la figura 2, la fem suministrada E es igual a la caída de voltaje a través de la resistencia 𝑅 es 𝑅𝐼 más
1
la caída de voltaje a través del condensador 𝐶 es 𝐶 Q, tenemos la ecuación para el circuito 2

1
𝑅𝐼 + Q = 𝐸 (𝑡 )
𝐶
𝑑𝑄
pero 𝐼 = 𝑄´(𝑡) = , la ecuación anterior se convierte en
𝑑𝑡
𝑑𝑄 1 1
𝑅 + Q = 𝐸 (𝑡) ⟺ 𝑅 𝑄´ + Q = 𝐸 (𝑡)
𝑑𝑡 𝐶 𝐶
la ecuación diferencial lineal para la carga instantánea.

Pero si en un circuito, se tiene las caídas de voltaje en la


resistencia 𝑅 es 𝑅𝐼, inductor 𝐿 es LI´ y en el condensador 𝐶 es
1
Q, por consiguiente, por la ley d de Kirchhoff, en un circuito
𝐶
eléctrico cerrado, la fem 𝐸 es la suma algebraica de todas las
caídas de voltaje en el circuito, es decir

1
𝐸 (𝑡) = 𝐿𝐼´(𝑡) + 𝑅𝐼 + 𝑄
𝐶
Pero 𝐼(𝑡) = 𝑄´(𝑡) ⟹ 𝐼´(𝑡) = 𝑄´´ (𝑡), entonces la ecuación
anterior se convierte en una ecuación lineal no homogénea de segundo orden con coeficientes constantes

1
𝐸 (𝑡) = 𝐿𝑄´´(𝑡) + 𝑅𝑄´(𝑡) + 𝑄(𝑡)
𝐶
resolviendo esta ecuación se halla la carga 𝑄 = 𝑄(𝑡) en el circuito.

Derivando esta ecuación con respecto al tiempo 𝑡,

𝑑 𝑑 1
𝐸 (𝑡) = (𝐿𝑄´´(𝑡) + 𝑅𝑄´(𝑡) + 𝑄(𝑡))
𝑑𝑡 𝑑𝑡 𝐶
1
𝐸´(𝑡) = 𝐿𝑄´´´(𝑡) + 𝑅𝑄´´(𝑡) + 𝑄´(𝑡)
𝐶
y reemplazando 𝐼 (𝑡) = 𝑄´(𝑡) ⟹ 𝐼´(𝑡) = 𝑄´´ (𝑡) ⟹ 𝐼´´ (𝑡) = 𝑄´´´(𝑡) en la ecuación anterior se convierte en una
ecuación lineal no homogénea de segundo orden con coeficientes constantes

1
𝐸´(𝑡) = 𝐿 𝐼´´(𝑡) + 𝑅 𝐼´(𝑡) + 𝐼(𝑡)
𝐶
resolviendo esta ecuación se halla la intensidad de corriente 𝐼 = 𝐼(𝑡) en el circuito.

Si un circuito eléctrico consta de un inductor L, resistencia R y un condensador C,


en este caso se supone que 𝐸 (𝑡) = 0, entonces por la ley de Kirchhoff, se tiene

1
𝐿𝑄´´(𝑡) + 𝑅𝑄´(𝑡) + 𝑄 (𝑡 ) = 0
𝐶
representa una ecuación lineal homogénea de segundo orden con coeficientes
constantes.

EJEMPLO

Un generador con una fem de 100 voltios se conecta en serie con una resistencia de 10 ohmios y un inductor de 2
henrios. Si el interruptor K se cierra en el tiempo 𝑡 = 0, establecer el modelo matemático para la corriente y determine
la corriente en el tiempo 𝑡.

Solución:

Por la ley de Kirchhoff, la fem suministrada E=100 voltios es igual a la caída de voltaje a través del inductor
𝐿 = 2 henrios es 2I´(t) más la caída de voltaje a través de la resistencia 𝑅 = 10 ohmios es 10I, tenemos la ecuación
diferencial para el circuito
𝑑𝐼
2 + 10𝐼 = 100 ⟺ 𝐼´ + 5 𝐼 = 50
𝑑𝑡
Como el interruptor se cierra en 𝑡 = 0, indica que no hay paso de corriente en el circuito, es decir 𝐼(0) = 0

Modelo Matemático

𝐼´ + 5 𝐼 = 50
{ 𝐼(0) = 0
𝐼 (𝑡 ) = ¿ ?

𝐼𝑒 ∫ 5𝑑𝑡 = ∫ 50𝑒 ∫ 5𝑑𝑡 𝑑𝑡 + 𝐶

𝐼𝑒 5𝑡 = ∫ 50𝑒 5𝑡 𝑑𝑡 + 𝐶

𝐼𝑒 5𝑡 = 10𝑒 5𝑡 + 𝐶

Usando la condición inicial 𝐼(0) = 0, para hallar la constante C

𝑡 = 0, 𝐼 = 0 ⟹ 𝐶 = −10

𝐼𝑒 5𝑡 = 10𝑒 5𝑡 − 10 ⟹ 𝐼(𝑡) = 10 − 10 𝑒 −5𝑡 = 10(1 − 𝑒 −5𝑡 ) intensidad de corriente en el circuito en el tiempo 𝑡.

𝐼 (𝑡) = 10(1 − 𝑒 −5𝑡 ) ∎

EJEMPLO

Una fuerza fem decayente 𝐸 = 200𝑒 −5𝑡 se conecta en serie el circuito con una resistencia de 20 ohmios y un
condensador de 0,001 faradios. Suponga que la carga 𝑄 = 0 en 𝑡 = 0, encuentere la carga y la corriente en cualquier
tiempo.

Solución:

Por la ley de Kirchhoff, la fem suministrada 𝐸 = 200𝑒 −5𝑡 voltios


es igual a la caída de voltaje a través de la resistencia 𝑅 = 20
ohmios es 20𝐼 más la caída de voltaje a través del condensador
1
𝐶 = 0,01 faradios es 𝑄 = 100𝑄 , tenemos la ecuación para
0,01
el circuito

1
𝑅𝐼 + Q = 𝐸 (𝑡) ⟹ 20𝐼 + 100 Q = 200𝑒 −5𝑡
𝐶
𝑑𝑄
Pero 𝐼 = 𝑄´(𝑡) = , entonces se tiene la ecuación diferencial lineal
𝑑𝑡

20𝑄´(𝑡) + 100 Q = 200𝑒 −5𝑡 ⟹ 𝑄´(𝑡) + 5 Q = 10𝑒 −5𝑡 ecuación lineal

MODELO MATEMÁTICO

𝑄´ + 5 𝑄 = 10𝑒 −5𝑡
𝑷𝑽𝑰: { 𝑄(0) = 0
𝐼(𝑡), 𝑄(𝑡) = ¿ ?

𝑄𝑒 ∫ 5𝑑𝑡 = ∫ 10𝑒 −5𝑡 𝑒 ∫ 5𝑑𝑡 𝑑𝑡 + 𝐾 = ∫ 10𝑒 −5𝑡 𝑒 5𝑡 𝑑𝑡 + 𝐾 = ∫ 10 𝑑𝑡 + 𝐶 = 10𝑡 + 𝐾

𝑄𝑒 5𝑡 = 10𝑡 + 𝐾

Usando la condición inicial 𝑄(0) = 0, para hallar la constante K:

𝑡 = 0, 𝑄=0 ⟹ 𝐾=0
𝑄𝑒 5𝑡 = 10𝑡

𝑄(𝑡) = 10𝑡𝑒 −5𝑡 la carga en el circuito en cualquier tiempo 𝑡.

Hallando la intensidad de corriente en el circuito

𝑄(𝑡) = 10𝑡𝑒 −5𝑡 ⟹ 𝑄´(𝑡) = 10𝑒 −5𝑡 − 50𝑡𝑒 −5𝑡

Pero 𝐼 = 𝑄´(𝑡) = 10𝑒 −5𝑡 − 50𝑡𝑒 −5𝑡

𝐼 (𝑡) = 10𝑒 −5𝑡 − 50𝑡𝑒 −5𝑡 intensidad de corriente en el circuito.

Por consiguiente, la intensidad de corriente y la carga en el circuito es

𝐼 (𝑡) = 10𝑒 −5𝑡 − 50𝑡𝑒 −5𝑡 , 𝑄(𝑡) = 10𝑡𝑒 −5𝑡 ∎

También podría gustarte