ABP – Ecuaciones
Diferenciales
Integrantes:
Benjamín Lorenzo Soto Mayta - Darcie Sol Ayala Tapia
PROBLEMA 4
PLANTEAMIENTO
Donde:
Segunda ley de Kirchhoff (LVK): • La caída de voltaje del inductor.
• Inductor: almacenan energía en
En cada instante de tiempo la suma algebraica de las forma de campo magnético.
tensiones alrededor de cualquier contorno cerrado del
circuito es cero. • La caída de voltaje del resistor.
• Resistor: limita o regula el flujo
de corriente eléctrica en los
circuitos eléctricos
• La caída de voltaje del capacitor.
• Capacitor: almacenan energía en
Ley de mallas: forma de campo eléctrico.
Es un método que nos permite hallar el valor de las
corrientes. Donde, una malla es una trayectoria cerrada o
lazo de un circuito que no contiene ninguna otra
trayectoria cerrada dentro de ella.
De acuerdo con Martínez et al. (2006) “Si el circuito
tiene N mallas se necesitan N ecuaciones simultaneas e
independientes, que surgirán de aplica la segunda ley de
Kirchhoff a cada una de las mallas”. (p. 36)
PLANTEAMIENTO
Condiciones iniciales y el contexto
asumido:
𝑀𝑎𝑙𝑙𝑎 1 : 𝐿 1′ + 𝑅 1 𝐼 = 𝐸 …(1)
𝐶𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠 : 𝐼 1 ( 0 )=0 𝐴 , 𝐼 2 ( 0 ) =0 𝐴
𝑞2
𝑀𝑎𝑙𝑙𝑎 2 : − 𝑅1 𝐼 + + 𝑅 2 𝐼 2= 0 …( 2)
𝐶
Asumimos que el circuito estaba abierto y en
t=0 se cierra instantáneamente.
2
′ 𝑅1 𝐸
𝐼1 + ( 𝐼 1 − 𝐼 2 )=
𝐿 𝐿
𝐼2
( 𝑅 1+ 𝑅 2 ) 𝐼 2′ − 𝑅 1 𝐼 1′ + =0
𝐶
RESOLUCIÓN 1: MÉTODO DE
ELIMINACIÓN
( 𝐷+25 ) 𝐼 1 −25 𝐼 2=50 …(1)
( 3 𝐷 +5 ) 𝐼 2 − 2 𝐷𝐼 1 =0 … (2)
Aplicando el operador a (1) y multiplicando por a (2), luego sumamos:
( 3 𝐷+5 )( 𝐷+25 ) 𝐼 1− ( 3 𝐷+5 ) 25 𝐼 2=250 … (3)
25 ( 3 𝐷+5 ) 𝐼 2− 25 × 2 𝐷𝐼 1= 0 …( 4)
¿ Hallando las raíces, por
fórmula general:
( 3 𝐷+5 )( 𝐷+25 ) 𝐼 1− 50 𝐷𝐼 1=250
( 3 𝐷 2+30 𝐷+125 ) 𝐼 1=250 … (5)
−5 𝑡
𝐼 1 ( 𝑡 ) =𝐶 1 𝑒 𝑠𝑖𝑛 ( )
5 √6
3
−5 𝑡
𝑡 +𝐶 2 𝑒 𝑐𝑜𝑠( )
5 √6
3
𝑡
Solución particular: 𝐼 1𝑝= 𝐴
Reemplazando en (5) : 125 𝐴=250 → 𝐴=2
c
𝐷𝑒 ( 1 ) : ( 𝐷+25 ) 𝐼 1 −50=25 𝐼 2
+ cos … (7)
𝑈𝑠𝑎𝑛𝑑𝑜 𝑙𝑎𝑠 𝑐𝑜𝑛𝑑𝑖𝑐𝑖𝑜𝑛𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠: 𝐼 1 ( 0 )=0 𝐴, 𝐼 2 ( 0 )=0 𝐴
𝐶 2=− 2
Reemplazando los valores de y en la funciones (6) y (7)
𝐼 1 ( 𝑡 )=4 √ 6 𝑒 −5 𝑡 𝑠𝑒𝑛 ( 5 √6
3 ) (
𝑡 −2 𝑒 −5 𝑡 c 𝑜𝑠 √ 𝑡 +2
5 6
3 )
𝐼 2 (𝑡)=
10 √ 6
3
5 6
3 (
𝑒 − 5 𝑡 𝑠𝑒𝑛 √ 𝑡 )
SEGUNDO MÉTODO:
TRANSFORMADA DE LAPLACE
• Segunda malla:
• Primera malla:
𝐼 1 ( 𝑠 ) 𝑆 − 𝑓 ( 0)+25 𝐼 1 ( 𝑠 ) − 25 𝐼 2 ( 𝑠 )=50 𝐿 { 1 }
( 50+25 ) 𝐼 2 ( 𝑠 ) 𝑆− ( 50+25 ) 𝐼 2 ( 0 ) −50 𝐼1 ( 𝑠 ) 𝑆+50 𝐼 1 ( 0 )+125 𝐼 2 ( 𝑠 )=0
50 3 𝐼 2 ( 𝑠 ) 𝑆 − 2 𝐼 1 ( 𝑠 ) 𝑆+ 5 𝐼 2 ( 𝑠 )=0
𝐼 1 ( 𝑠 ) [ 𝑆+25 ] −0 − 25 𝐼 2 ( 𝑠 )= … (𝛼)
𝑆 𝐼 2 ( 𝑠 ) [ 3 𝑆+5 ] − 2 𝐼 1 ( 𝑠 ) 𝑆=0 … ( 𝛽 )
Luego, para tener una ecuación respecto a I2: 100
𝐼 2 ( 𝑠) =
[ 50
]
2
3 𝑆 +30 𝑆+125
(𝟐 𝑺)× 𝐼 1 ( 𝑠 ) [ 𝑆+ 25 ] − 25 𝐼 2 ( 𝑠 ) =
𝑆
… (𝛼 )
(𝑺 +𝟐𝟓) × [ 𝐼 2 ( 𝑠 ) [ 3 𝑆+5 ] − 2 𝑆 𝐼 1 ( 𝑠 ) =0 ] … ( 𝛽 )
¿
Se dividirá entre 3 y se
completará cuadrados.
𝐼 2 ( 𝑠 ) [ − 25 × ( 2 𝑆 ) ] + 𝐼 2 ( 𝑠 ) [ ( 3 𝑆+5 ) × ( 𝑆+25 ) ] =100
100
𝐼 2 ( 𝑠 ) [ 3 𝑆 2+5 𝑆+75 𝑆 − 50 𝑆+125 ] =100 3
𝐼 2 ( 𝑠) =
125
𝑆2 +10 𝑆+
𝐼 2 ( 𝑠 ) [ 3 𝑆 +30 𝑆+125 ] =100
2
3
100
3 100
𝐼 2 ( 𝑠) =
2
𝑆 +10 𝑆+5 +
2 12 5
−5
2
𝐼 2 ( 𝑠) =
3
×sin (√ 503 𝑡 )× 𝑒 −5𝑡
√
3 50
100 3
3
𝐼 2 ( 𝑠) =
Racionalizamos:
( 𝑆+5 )2 + 50
3
Se multiplicara y dividirá por un factor 𝐼 2 ( 𝑠) =
10 √ 6 −5 𝑡
3
𝑒 sin(5 √6
3
𝑡 )
𝐼 2 ( 𝑠) =
100
3
×
√ 50
3
√
2 50 50
( 𝑆+5 ) +
3 3
100
𝐼 2 ( 𝑠) =
3
×si n (√ ) 50
𝑡 ×𝑒 −5 𝑡
√
50 3
3
Luego, para tener una ecuación respecto a I1:
[
( 𝟑 𝑺 +𝟓 ) × 𝐼 1 ( 𝑠 ) [ 𝑆+25 ] −25 𝐼 2 ( 𝑠 )=
50
𝑆 ]… (𝛼)
¿
(𝟐𝟓)× [ 𝐼 2 ( 𝑠 ) [ 3 𝑆+5 ] − 2 𝐼 1 ( 𝑠 ) 𝑆= 0 ] … ( 𝛽 )
250
𝐼 1 ( 𝑠 ) [ ( 𝑆+25 ) × ( 3 𝑆+5 ) ] −50 𝑆 𝐼 1 ( 𝑠 )=150+
𝑆
250
𝐼 1 ( 𝑠 ) [ 3 𝑆 +7 5 𝑆+5 𝑆+125 −50 𝑆 ] =150+
2
𝑆
250
𝐼 1 ( 𝑠 ) [ 3 𝑆 +30 𝑆+125 ] =150+
2
𝑆
1
150 250 3
𝐼 1 ( 𝑠 )= + ×
( 3 𝑆 +30 𝑆 +125 )
( ) 1
2
125
𝑆 3 𝑆 2 +1 0 𝑆 +
3 3
250
50 3
𝐼 1 ( 𝑠 )= +
( ) 𝑆 ( 𝑆 +30 𝑆 +125 )
2
2 125
𝑆 +10 𝑆 +
3
250 250
50 3 50 3
𝐼 1 ( 𝑠 )= + = +
( (
)+30 𝑆+125 )
( ) 𝑆 (( 𝑆+5 ) + 503 )
2
125
𝑆 +10 𝑆+2 𝑆 𝑆 ( 𝑆+5 ) + 50
2 2
3 3
250
3 𝐴 𝐵𝑆 +𝐶
= +
( 2
𝑆 ( 𝑆 +5 ) +
50
3
𝑆
)( 𝑆 +5 ) 2+ 50
3 ( ) Se aplicará
250
3
2
(
= 𝐴 × ( 𝑆+5 ) +
50
3
+ ( 𝐵𝑆+ 𝐶 ) × ( 𝑆 ) ) fracciones parciales
Para S = 0 Para S = -5
250
3 (
= 𝐴 × ( 0+5 )2 +
50
3 )
250
3 (
= 𝐴 × ( − 5+5 )2 +
50
3 )
+ ( 𝐵× ( −5 ) +𝐶 ) × ( − 5 )
250
= 𝐴×
125
250
3
=(2) × ( 50
3 )+ ( 25 𝐵 − 5 𝐶 )
3 3 5 0= ( 25 𝐵 − 5 𝐶 )
𝐴=2 1 0=5 𝐵− 1𝐶 … ( 𝛿 )
Para S = -4
250
3 (
= 𝐴 × ( − 4 +5 )2 +
50
3 )
+ ( 𝐵× ( − 4 ) +𝐶 ) × ( − 4 )
250
3
=(2) × 1+
50
3( )
+ ( 16 𝐵 − 4 𝐶 )
4 8=( 16 𝐵 − 4 𝐶 )
1 2=4 𝐵 −1 𝐶 …(𝛾)
Luego por medio de un sistema de ecuaciones despejamos las variables:
1 2=4 𝐵 −1 𝐶 …(𝛾) 1 2=4 ( − 2 ) −1 𝐶 …(𝛾)
(−)
1 0=5 𝐵− 1𝐶 … ( 𝛿 ) 2 0=−1 𝐶
2= − 𝐵
C =− 20
𝐵=− 2
50 2 −2 𝑆 − 20
𝐼 1 ( 𝑠 )= + +
( 2
𝑆 +10 𝑆+
125
3 ) 𝑆
(
( 𝑆+5 )2+
50
3 )
{( }
50 2 2𝑆 20
𝐿
−1
{𝐼 1 ( 𝑠 ) }= 𝐿− 1 + − −
( 𝑆+5 )2+
50
3 ) 𝑆
( ( 𝑆+5 )2+
50
3 ) ( ( 𝑆+5 )2+
50
3 )
{( √
} {( √
} {( √
}
50 50 50
−1
𝐼1 ( 𝑡 ) = 𝐿
50
×
3
+2𝐿
−1 1
{ }
−2 𝐿
−1 𝑆
+
5
−
5
×
3
−𝐿
−1 20
×
3
) √ ) ( ) ( ) √ ) √
2 50 50 𝑆 2 50 2 50 2 50 50 2 50 50
( 𝑆+5 ) + ( 𝑆+5 ) + ( 𝑆+5 ) + ( 𝑆+5 ) + ( 𝑆+5 ) +
3 3 3 3 3 3 3 3
𝐼1 ( 𝑡 ) =
50
𝑠𝑖𝑛 (√ ) 50
𝑡 ×𝑒 −5 𝑡 + 2− 2 𝑐𝑜𝑠 (√ ) 50
𝑡 × 𝑒− 5 𝑡 +
2× 5
𝑠𝑖𝑛 (√ ) 50
𝑡 ×𝑒 −5 𝑡 −
20
𝑠𝑖𝑛 (√ )50
𝑡 ×𝑒 −5 𝑡
√ √ √
50 3 3 50 3 50 3
3 3 3
Racionalizamos:
𝐼 1 ( 𝑡 ) =2 −2 𝑒− 5 𝑡 𝑐𝑜𝑠 ( 5 √6
3 )
𝑡 + 4 √ 6 𝑒 −5 𝑡 𝑠𝑖𝑛 √ 𝑡
5 6
3 ( )
ÁNALISIS DE
LOS
RESULTADOS
RESOLUCIÓN : MÉTODO DE
ELIMINACIÓN PARA EL CASO DE
VARIACIÓN PARAMETRO
𝐿=2 h𝑒𝑛𝑟𝑖𝑒𝑠
𝑅1=25 𝑜h𝑚𝑠 𝐶=0,008 𝑓𝑎𝑟𝑎𝑑𝑠
𝑅2 =25 𝑜h𝑚𝑠 𝑉 =100 𝑣𝑜𝑙𝑡𝑠
Malla 2:
Malla 1:
𝑑 𝐼2 𝑑 𝐼1
𝑑 𝐼1 ( 25+25 ) − 25 +125 𝐼 2=0
2 +25 𝐼 1 −25 𝐼 2= 100 𝑑𝑡 𝑑𝑡
𝑑𝑡
2𝐼1 (𝑠 ) 𝑆−2𝐼1 ( 0 )+25𝐼 1 (𝑠 ) − 25𝐼2 ( 𝑠 )=100𝐿 {1}… (𝑎) (25+25 ) 𝐼 2 (𝑠 ) 𝑆− (25+25 ) 𝐼 2 (0) −25𝐼 1( 𝑠) 𝑆+25𝐼1 (0)+125𝐼 2( 𝑠)=0…(𝑏)
2 𝐼 1 ( 𝑠 ) 𝑆+25 𝐼 1 ( 𝑠 ) −25 𝐼 2 ( 𝑠 )=100 𝐿 {1 } …(𝑎)
−25 𝐼 2 ( 𝑠 ) 𝑆+ 𝐼 2 ( 𝑠 ) [ ( 2 𝑆+5 ) × ( 2 𝑆+ 25 ) ] =100
( 50 ) 𝐼 2 ( 𝑠 ) 𝑆 −25 𝐼 1 ( 𝑠 ) 𝑆+125 𝐼 2 ( 𝑠 )=0 …(𝑏)
𝐼 2 ( 𝑠 ) [ 4 𝑆 +35 𝑆 +125 ] =100
2
Factorizamos e en y : 100
𝐼 2 ( 𝑠) =
100 2
4 𝑆 +35 𝑆 +125
𝐼 1 ( 𝑠 ) [ 2 𝑆+25 ] − 25 𝐼 2 ( 𝑠 ) = …(𝑎)
𝑆 Factorizamos y completamos cuadrados:
𝐼 2 ( 𝑠 ) [ 50 𝑆+125 ] −25 𝐼 1 ( 𝑠 ) 𝑆=0 …(𝑏) 25 ×
5 √ 31
100 8 1
𝐼 2 ( 𝑠) = = ×
𝐼 2 ( 𝑠 ) [ 2 𝑆+5 ] − 𝐼 1 ( 𝑠 ) 𝑆=0 … (𝑏)
( ) ) ( ) 5 √ 31
2 2
35 775
( 5 √ 31
2
4 𝑆+ + 35
𝑆+ +
8 16 8 8 8
Multiplicamos por factores que permitan
dejar todo en una sola incógnita: Aplicamos Laplace inversa:
[
( 𝑆 ) × 𝐼 1 ( 𝑠 ) [ 2 𝑆+25 ] − 25 𝐼 2 ( 𝑠 )=
100
𝑆 ]…(𝑎)
¿ 𝐼 2(𝑡)=
40 √ 31
31
×𝑒
− 35 𝑡
8
sen ( 5 √31
8
𝑡 )
( 2 𝑆+25 ) × [ 𝐼 2 ( 𝑠 ) [ 2 𝑆+5 ] − 𝐼 1 ( 𝑠 ) 𝑆=0 ] …(𝑏)
−25 𝐼 2 ( 𝑠 ) 𝑆+ 𝐼 2 ( 𝑠 ) [ ( 2 𝑆+5 ) × ( 2 𝑆+ 25 ) ] =100
Se usará (b) para encontrar el valor de I1:
[ 2 𝑆+5 ]
𝐼 2 ( 𝑠 ) [ 2 𝑆+5 ] − 𝐼 1 ( 𝑠 ) 𝑆=0 …(𝑏) 𝐼 1 ( 𝑠 ) 𝑆=𝐼 2 ( 𝑠 ) [ 2 𝑆+5 ] 𝐼 1 ( 𝑠 )= 𝐼 2 (𝑠 )
𝑆
[ 2 𝑆+5 ] 100 𝐴 𝐵𝑆+𝐶
𝐼 1 ( 𝑠 )= × = +
𝑆 4 𝑆 +35 𝑆+125 𝑆 4 𝑆 + 35 𝑆+125
2 2
200 𝑆+500= 𝐴 × ( 4 𝑆 + 35 𝑆+125 ) + ( 𝐵𝑆+𝐶 ) × 𝑆
2
𝐴= 4 𝐵=− 16 𝐶=−16
(( ) (( )
4 −16 𝑆+60 4 16 𝑆 60 1
𝐼 1 ( 𝑠 )= + = − +
) ( ) ) ( )
𝑆 4 𝑆 +35 𝑆+125 𝑆
2
4 2
4 2
+ √ + √
2 2
35 5 31 35 5 31
𝑆+ 𝑆+
8 8 8 8
( ( 𝑆 + 35 )
) ((
35
−
)
4 8 8 15
𝐼 1 ( 𝑠 )= −4 +
+( ) ) +( )
𝑆 2 2
(𝑆 + 35 ) 5 √ 31 √ 31
2 2
35 5
𝑆+
8 8 8 8
(( ( )
)) (
35 35
𝑆+ ×− 4
( )
4 8 8 15
𝐼 1 ( 𝑠 )= −4 − +
) ( ) ( ) ( ) ( )
𝑆 35 2
5 √ 31
2
35 2 5 √ 31
2
35 2
5 √ 31
2
𝑆+ + 𝑆+ + 𝑆+ +
8 8 8 8 8 8
(( ( )
)
35 5 √ 31
𝑆+
(( )
4 8 65 1 8
𝐼 1 ( 𝑠 )= −4 + × ×
5 √ 31
) +( ) ) ( )
𝑆 35
2
5 √ 31
2
2 35
2
5 √ 31
2
𝑆+ 𝑆+ +
8 8 8 8 8
( ( 𝑆 + 35
8 )
) (( )
5 √ 31
4 52 √ 31 8
𝐼 1 ( 𝑠 )= −4 + ×
( 8) ( 8 ) ) ( )
𝑆 35 5 √ 31
2 2
31 35
2
5 √ 31
2
𝑆 + + 𝑆+ +
8 8
Finalmente usamos:
( ) ( )
35 35
−
8
𝑡 5 √ 31 52 √ 31 − 8
𝑡 5 √ 31
𝐼 1 ( 𝑠 )= 4 − 4 𝑒 cos 𝑡 + 𝑒 sin 𝑡
8 31 8
Corriente I2 sin variación de parámetro
Corriente I2 con variación de parámetro
ÁNALISIS
DE LOS
RESULTAD
OS
Corriente I1 con variación de parámetro
ÁNALISIS
DE LOS
RESULTAD
Corriente I1 sin variación de parámetro OS
REFERENCIAS
Martínez, A., Domínguez, E., & Díaz, A. (2006). Sistemas y circuitos (2.a ed., Vol. 1).
Universidad Politécnica de Cartagena.
https://repositorio.upct.es/bitstream/handle/10317/7661/syc.pdf?
sequence=1&isAllowed=y
D. Young, H., & A. Roger, F. (2013). Física universitaria con física moderna (13.a ed.,
Vol. 2). Pearson.
GRAC
IAS
Por su atención