lOMoARcPSD|25287928
ECUACIONES DIFERENCIALES EJE NUMERO 4 TAREA
calcula integral (Fundación Universitaria del Área Andina)
Escanea para abrir en Studocu
Studocu no está patrocinado ni avalado por ningún colegio o universidad.
Descargado por Jonathan Caro Mora (
[email protected])
lOMoARcPSD|25287928
ECUACIONES DIFERENCIALES
TRANSFORMADAS DE LAPLACE EN LA SOLUCIÓN DE ECUACIONES
DIFERENCIALES
ACTIVIDAD EVALUATIVA EJE 4
MICHAEL NICOLAS GOMEZ CURREA
FABIAN ALEXIS LEUDO RAMOS
JOSÉ JULIÁN NIÑO PRADA
SERGIO ANDRÉS RAMÍREZ ACOSTA
SILVIA REBECA VEGA RIAÑO
Docente
FUNDACIÓN UNIVERSITARIA DEL AREANDINA
2021
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
INTRODUCCION
En la actualidad la aplicación de circuitos eléctricos resulta común, dado que muchos de
los aparatos que nos rodean en nuestra vida, como el celular, la televisión y la computadora,
entre muchos otros, requieren de circuitos eléctricos simples, combinados y complejos para
su funcionamiento. (MAXIMILIANO, 2013)
En la presente actividad tenemos dos situaciones problema las cuales se pretenden resolver
mediante la transformada de lapace, esta es una forma de dar solución a este tipo de
ecuaciones de una manera sencilla en comparación con las demás.
La transformada de Laplace es de gran importancia en la ingeniería ya que permite reducir
ecuaciones diferenciales ordinarias con coeficientes constantes a simples expresiones
algebraicas de sencilla resolución.
La transformada de Laplace de una función f(t) definida para todos los números positivos
t ≥ 0 es la función F(s), definida por:
∝
𝑭(𝒔) = 𝑳 {𝒇(𝒕)} = ∫ 𝒆−𝒔𝒕 𝒇(𝒕) 𝒅𝒕.
𝟎
siempre y cuando la integral este definida.
La función original f(t) se conoce como transformada inversa o inversa de F(s), es decir:
𝑳−𝟏 {𝑭(𝒔)} = 𝒇(𝒕)
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
Situación problémica 1
Circuito eléctrico De los circuitos usados en los contextos y aplicaciones de la Ingeniería
y otras disciplinas, el circuito RLC (figura 1) aparece con frecuencia, ya que cada circuito
real tiene una cierta resistencia finita. Un circuito RLC se compone de los elementos pasivos:
resistencia, bobina y condensador.
Un circuito RLC en serie tiene una fuente de voltaje dada por
𝑽(𝒕) = 𝒔𝒆𝒏 𝟏𝟎𝟎𝒕 ,
𝑢𝑛 𝑟𝑒𝑠𝑖𝑠𝑡𝑜𝑟 𝑑𝑒 𝟎. 𝟎𝟐 Ω,
𝑢𝑛 𝑖𝑛𝑑𝑢𝑐𝑡𝑜𝑟 𝑑𝑒 𝟎. 𝟎𝟎𝟏 𝑯
𝑢𝑛 𝑐𝑎𝑝𝑎𝑐𝑖𝑡𝑜𝑟 𝑑𝑒 𝟐 F
Si la corriente y la carga iniciales en el capacitor son iguales a cero, determinar la corriente
en el circuito para 𝒕 > 𝟎.
La ecuación para la corriente es:
𝑑𝐼 𝑑 2 𝐼 𝐼 𝑑𝑉
𝑅 +𝐿 2+ =
𝑑𝑡 𝑑𝑡 𝐶 𝑑𝑡
𝑑𝐼 𝑑2𝐼
(0,02) (
+ 0,001 ) + 0,5𝐼 = 100cos (100𝑡)
𝑑𝑡 𝑑𝑡 2
𝑑2𝐼 𝑑𝐼
(0,001) + ( 0,02 ) + 0,5𝐼 = 100cos (100𝑡)
𝑑𝑡 2 𝑑𝑡
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
𝑑2𝐼 (0,02) 𝑑𝐼 0,5 100 cos (100𝑡)
+ + 𝐼=
𝑑𝑡 2 (0,001) 𝑑𝑡 (0,001) (0,001)
𝑑2𝐼 𝑑𝐼
2
+ 20 + 500𝐼 = 100000cos (100𝑡)
𝑑𝑡 𝑑𝑡
Las condiciones iniciales son 𝐼 (0) = 0 𝑦 𝐼′ (0) = 0.
Usando transformada de Laplace:
𝑠
𝑠 2 𝐼 (𝑠) − 𝑠𝐼 (0) − 𝐼 ′ (0) + 20𝑠𝐼 (𝑠) − 20𝐼(0) + 500𝐼(𝑠) = 100000( )
𝑠2 + 1002
𝑠
𝑠 2 𝐼(𝑠) + 20𝑠𝐼 (𝑠) + 500𝐼(𝑠) = 100000 ( )
𝑠2 + 1002
𝑠
(𝑠 2 + 20𝑠 + 500)𝐼(𝑠) = 100000 ( )
𝑠2 + 1002
100000 𝑠
𝐼(𝑠) = ( 2 )
(𝑠 2 + 20𝑠 + 500) 𝑠 + 1002
𝑠 2 + 20𝑠 + 500 = 𝑠 2 + 20𝑠 + 400 + 100 = (𝑠 + 10)2 + 400
100000 𝑠 𝐴𝑠 + 𝐵 𝐶𝑠 + 𝐷
𝐼(𝑠) = ( 2 )= + 2
((𝑠 + 10) + 400) 𝑠 + 100
2 2 (𝑠 + 10) + 400 𝑠 + 1002
2
Entonces,
(𝐴𝑠 + 𝐵 )(𝑠 2 + 1002 ) + (𝐶𝑠 + 𝐷)((𝑠 + 10)2 + 400) = 100000
𝐴+𝐶= 0
𝐵 + 20𝐶 + 𝐷 = 0
10000𝐴 + 500𝐶 + 20𝐷 = 100000
10000𝐵 + 500𝐷 = 0
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
3800 4000 3800 80000
𝐴= ;𝐵 = − ;𝐶 = − ;𝐷 =
377 377 377 377
De esta forma,
3800 4000 3800 80000
𝑠− − 𝑠+
𝐼(𝑠) = 377 2 377 + 377 377
(𝑠 + 10) + 400 𝑠 2 + 1002
3800 4000 20
𝑠 − 377 3800 𝑠 − 19
𝐼1 (𝑠) = 377 = ( )
(𝑠 + 10)2 + 400 377 (𝑠 + 10)2 + 400
20
3800 𝑠 + 10 − 10 − 19
𝐼1 (𝑠) = ( )
377 (𝑠 + 10)2 + 400
210
3800 𝑠 + 10 − 19
𝐼1 (𝑠) = ( + )
377 (𝑠 + 10)2 + 400 (𝑠 + 10)2 + 400
3800 𝑠 + 10 3800 210 1 20
𝐼1 (𝑠) = ( ) − ( ) ( ) ( ) ( )
377 (𝑠 + 10)2 + 400 377 19 20 (𝑠 + 10)2 + 400
3800 𝑠 + 10 2100 20
𝐼1 (𝑠) = ( )−( )( )
377 (𝑠 + 10) + 400
2 377 (𝑠 + 10)2 + 400
3800 −1 𝑠 + 10 2100 −1 20
𝐼1 (𝑠) = 𝐿 { }−( )𝐿 { }
377 (𝑠 + 10) + 400
2 377 (𝑠 + 10)2 + 400
3800 −10𝑡 2100 −10𝑡
𝐼1 (𝑠) = 𝑒 𝑐𝑜𝑠20𝑡 − 𝑒 𝑠𝑒𝑛20𝑡
377 377
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
3800 80000 3800 80000
− 𝑠+ − 𝑠+
𝐼2 (𝑠) = 377 377 = 377 377
𝑠 2 + 1002 (𝑠 − 100𝑖 )(𝑠 + 100𝑖 )
𝑎 = 100𝑖 = 𝛼 + 𝑖𝛽 → 𝛼 = 0; 𝛽 = 100
𝑎̅ = −100𝑖
3800 80000
− 𝑠+
𝑄(𝑠) = 377 377
(𝑠 − 100𝑖 )
3800 80000
− 377 (100) + 377 1900 400𝑖
𝑄(𝑎) = 𝑄(100𝑖 ) = =− − = 𝑄1 + 𝑄2
(100𝑖 + 100𝑖 ) 377 377
1900 400
𝑄1 = − ; 𝑄2 = −
377 377
3800 80000
− 377 𝑠 + 377
−1
𝐼2 (𝑡) = 𝐿 { } = 2𝑒 𝛼𝑡 (𝑄1 𝑐𝑜𝑠𝛽𝑡 − 𝑄2 𝑠𝑒𝑛𝛽𝑡)
(𝑠 2 + 1002 )
Entonces,
1900 400 3800 800
𝐼2 (𝑡) = 2𝑒 0𝑡 (− 𝑐𝑜𝑠100𝑡 + 𝑠𝑒𝑛100𝑡) = − 𝑐𝑜𝑠100𝑡 + 𝑠𝑒𝑛100𝑡
377 377 377 377
Por lo tanto, a solución es
3800 −10𝑡 2100 −10𝑡 3800 800
𝐼(𝑡) = 𝑒 𝑐𝑜𝑠20𝑡 − 𝑒 𝑠𝑒𝑛20𝑡 − 𝑐𝑜𝑠100𝑡 + 𝑠𝑒𝑛100𝑡
377 377 377 377
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
Situación problémica 2
Esta aplicación de las transformadas de Laplace en la solución de ecuaciones diferenciales
de segundo orden se presenta con el propósito de determinar el movimiento de oscilación de
dos resortes acoplados de forma horizontal mediante una ecuación lineal de segundo grado y
transformadas de Laplace, considerando también conceptos básicos de la Ley de Hooke y las
leyes de Newton. En una superficie horizontal suave, una masa
𝒎𝟏 = 𝟏 𝒌𝒈
está unida a una pared fija mediante un resorte con constante de resorte
𝒌𝟏 = 𝟐 𝑵/𝒎 .
Por su parte, otra masa
𝒎𝟐 = 𝟐 𝒌𝒈
está unida al primer objeto mediante un resorte con constante de resorte
𝒌𝟐 = 𝟒 𝑵/𝒎.
Los objetos están alineados en forma horizontal, de modo que los resortes tengan su
longitud natural.
Si ambos objetos se desplazan 𝟑 𝒎𝒆𝒕𝒓𝒐𝒔 a la derecha de sus posiciones de equilibrio y
luego se liberan, ¿Cuáles son las ecuaciones de movimiento de los dos objetos?
𝑭𝟏 = −𝒌𝟏 𝒙 𝑭𝟐 = 𝒌𝟐 (𝒚 − 𝒙)
𝑭𝟑 = −𝒌𝟐 (𝒚 − 𝒙)
Aplicando la 2da ley de Newton a los objetos:
𝑑2𝑥
𝑚 1 2 = 𝐹1 + 𝐹2 = −𝑘1 𝑥 + 𝑘2 (𝑦 − 𝑥 )
𝑑𝑡
𝑑2𝑥
𝑚1 + (𝑘1 + 𝑘2 )𝑥 − 𝑘2 𝑦 = 0
𝑑𝑡 2
y,
𝑑2𝑦
𝑚2 = 𝐹3 = −𝑘2 (𝑦 − 𝑥 )
𝑑𝑡 2
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
𝑑2𝑦
𝑚2 + 𝑘2 𝑦 − 𝑘2 𝑥 = 0
𝑑𝑡 2
Resolviendo el sistema,
𝑑2𝑥
+ 6𝑥 − 4𝑦 = 0
𝑑𝑡 2
𝑑2𝑦 𝑑2𝑦
2 + 4𝑦 − 4𝑥 = 0 → + 2𝑦 − 2𝑥 = 0
𝑑𝑡 2 𝑑𝑡 2
Condiciones iniciales: 𝑥 (0) = 3, 𝑥′(0) = 0; 𝑦(0) = 3, 𝑦′(0) = 0.
La transformada de Laplace para cada ecuación es:
s 2 X(s) − sx(0) − x ′ (0) + 6X(s) − 4Y(s) = 0
s 2 Y(s) − sy(0) − y′(0) + 2Y(s) − 2X(s) = 0
Sustituyendo las condiciones iniciales:
s 2 X(s) − 3s + 6X(s) − 4Y(s) = 0 → (𝑠 2 + 6)X(s) − 4Y(s) = 3s
s 2 Y(s) − 3s + 2Y(s) − 2X(s) = 0 → −2𝑋(𝑠) + (𝑠 2 + 2)𝑌(𝑠) = 3𝑠
𝑠2 + 6 −4 |
∆= | = 𝑠 4 + 8𝑠 2 + 4
−2 𝑠2 + 2
3𝑠 −4
| | 3
3𝑠 𝑠 2 + 2 = 3𝑠 + 18𝑠
𝑋(𝑠) =
∆ 𝑠 4 + 8𝑠 2 + 4
2
|𝑠 + 6 3𝑠
| 3
𝑌(𝑠) = −2 3𝑠 = 3𝑠 + 24𝑠
∆ 𝑠 4 + 8𝑠 2 + 4
3𝑠 3 + 18𝑠 𝐴𝑠 + 𝐵
𝑋(𝑠) = 4 2
=
𝑠 + 8𝑠 + 4 4
(𝐴𝑠 + 𝐵 )(𝑠 2 − 2) + 4(𝐶𝑠 + 𝐷) = 2𝑠 3 + 42𝑠
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
𝐴=2
𝐵 =0
2𝐶 − 𝐴 = 42
2𝐷 − 𝐵 = 0
𝐴 = 2; 𝐵 = 0; 𝐶 = 22; 𝐷 = 0
Entonces;
2𝑠 22𝑠
𝑋 (𝑠) = + 2
4 𝑠 −2
2𝑠 22𝑠 𝛿(0)
𝑥 (𝑡) = 𝐿−1 { + 2 }= + 11𝑒 −√2𝑡 + 11𝑒 √2𝑡
4 𝑠 −2 2
12𝑠 3 + 36𝑠 𝐴𝑠 + 𝐵 𝐶𝑠 + 𝐷
𝑌(𝑠) = = + 2
4( 𝑠 2 − 2) 4 𝑠 −2
(𝐴𝑠 + 𝐵 )(𝑠 2 − 2) + 4(𝐶𝑠 + 𝐷) = 12𝑠 3 + 36𝑠
𝐴 = 12
𝐵 =0
2𝐶 − 𝐴 = 36
2𝐷 − 𝐵 = 0
𝐴 = 12; 𝐵 = 0; 𝐶 = 24; 𝐷 = 0
Entonces;
12𝑠 24𝑠
𝑌(𝑠) = + 2
4 𝑠 −2
12𝑠 24𝑠
𝑦 (𝑡) = 𝐿−1 { + 2 } = 3𝛿(0) + 12𝑒 −√2𝑡 + 12𝑒 √2𝑡
4 𝑠 −2
Descargado por Jonathan Caro Mora ([email protected])
lOMoARcPSD|25287928
BIBLIOGRAFIA
➢ G. James, Matemáticas Avanzadas para Ingeniería, Pearson Educación, 2002, pp.
130-135.
➢ R.A. Serway y R.J. Beichner, Física II Para Ciencias e Ingeniería, McGraw-Hill,
2001, pp. 877-878
➢ G. Calandrini, Guía de Definiciones y Teoremas estudiados en el curso de
funciones de Variable Compleja.1er. cuatrimestre 2011, pp. 56-63.
Descargado por Jonathan Caro Mora ([email protected])