INSTITUTO TECNOLÓGICO DE CELAYA

G
Ecuaciones diferenciales

PROFESOR
Ramos Beltrán José Alfredo

ALUMNO
Gaviña Vázquez Byron Eduardo
21030323

César Eduardo Martínez Arredondo

21030441

Fecha
Tabla de contenido
2. Introducción.....................................................................................................................................................3

3. Marco Teórico.......................................................................................................................................................6
4. Desarrollo
5. Conclusiones........................................................................................................................................................10
Conclusión de Gaviña Vázquez Byron Eduardo..................................................................................................10
Conclusión de César Eduardo Martinez Arredondo...........................................................................................10

6........................................................................................................................... Referencias Bibliográficas

..............................................................................................................................................................................10
 1. Introducción
Las ecuaciones diferenciales son herramientas fundamentales para modelar fenómenos
dinámicos en ciencias e ingeniería. En el contexto de circuitos eléctricos, las leyes de Kirchhoff y
la Ley de Ohm permiten establecer relaciones entre corriente, voltaje, y resistencia, las cuales
dan lugar a sistemas de ecuaciones diferenciales que describen el comportamiento temporal de
las cargas, corrientes y voltajes en los componentes.

El presente trabajo aborda el análisis de un circuito RLC en paralelo, planteando el modelo

matemático que lo describe mediante las leyes antes mencionadas. El propósito principal es
resolver este modelo utilizando transformadas de Laplace y métodos numéricos asistidos por
software para obtener las soluciones específicas. El uso de herramientas computacionales,
como Python o MATLAB, es esencial para resolver ecuaciones diferenciales complejas de
manera eficiente y para visualizar las soluciones, proporcionando una comprensión más
profunda del fenómeno físico.
 2. Marco Teórico

El modelo del circuito RLC en paralelo se basa en las siguientes leyes fundamentales:
1. Ley de Kirchhoff de Corrientes (KCL): En un nodo, la suma de las corrientes
entrantes es igual a la suma de las corrientes salientes.

2. Ley de Kirchhoff de Tensiones (KVL): La suma algebraica de las tensiones en

un lazo cerrado es igual a cero.

3. Ley de Ohm: Relaciona la corriente y el voltaje en un resistor mediante V =IR

En este circuito, se modelan las siguientes componentes:

 Capacitor (C): Almacena energía eléctrica en forma de carga. Su relación es
q=CV , donde q es la carga.

 Inductor (L): Almacena energía en forma de campo magnético, relacionado

por V =IR

 Resistores (R1 y R2): Disipan energía según V =IR

El sistema se describe por dos ecuaciones diferenciales derivadas de las leyes

mencionadas, considerando las corrientes y tensiones en cada rama. Estas ecuaciones
vinculan el comportamiento de la carga q(t), la corriente i(t), y el voltaje aplicado E(t).

3. Planteamiento del Problema

a) Aplicando las leyes de Ohm y Kirchhoff muestre que el sistema de ecuaciones

diferenciales que modela la carga en el capacitor (q)t y la corriente i 3 ( t ) en la red
eléctrica esta dado por
dq 1
R1 + q+ R 1 i 3 ( t )= E ( t )
dt C
di 3 −1
L1 + R 2 i 3 ( t )= q=0
dt C
Paso 1: Planteamiento del sistema de ecuaciones diferenciales
A partir de las leyes de Kirchhoff:
1. En la rama R1 y C:

dq q 1
R1 + q=E(t)
dt C

2. En la rama R2 y L:
2
d q dq q
L 2
+ R2 − =0
dt dt C

Paso 2: Resolución con Transformadas de Laplace

1. Tomando la transformada de Laplace:

dq q 1
 Para R1 + =E ( t ) :R1 sQ ( s ) + Q(s )=E ( s )
dt C C
2
d q dq q
 Para L 2
+ R2 − =0
dt dt C
2 1
L s Q( s)+ R2 sQ (s )− Q(s)=0
C

2. Resolviendo para Q(s), obtenemos una función de transferencia que depende de

E(s). Luego se aplica la transformada inversa para encontrar q(t).

3. Paso 3: Resolución numérica

Se resolvió el sistema con parámetros dados
(L=1 H , R 1=5 Ω , R 2=10 Ω ,C=0.01 F , E (t)) usando métodos numéricos. Los
resultados obtenidos se presentaron gráficamente para i3(t) en el intervalo
0 ≤ t ≤ 10 y se construyó una tabla con los valores de q ( t ) e i 3 (t) para los primeros
5 segundos.
Usando las leyes de Kirchhoff y Ohm:
Para el circuito paralelo, la corriente en cada rama suma hasta la corriente total:
i 1+ i2 +i 3=i
Aplicamos la Ley de Kirchhoff de tensiones:
 En la rama del capacitor (R1 y C):
q
R1 i1 + =E (t)
C
dq
Donde q(t) es la carga en el capacitor y i 1=
dt
En la rama del inductor (R2 y L):
di 3 q
L + R2 i 3=
dt C

dq
Sustituyendo i 1= y reorganizando:
dt
 Primera ecuación
dq q
R1 + =E(t)
dt C

 Segunda Ecuación
2
d q dq q
L 2
+ R2 − =0
dt dt C

b) Determinación de la carga q(t)

 Con los valores dados:

 L=1 H
 R1=5 Ω
 C=0.01 F

{
0 0 ≤ t<1
 E ( t ) 120 sin ( t ) 1 ≤t <5
0 5≤ t
La solución requiere:
1. Resolver el sistema diferencial usando condiciones iniciales (q ( 0 )=0 , i ( 0 )=0)
2. Usar métodos analíticos (Transformada de Laplace) o numéricos (CAS).

Vamos a resolver el problema utilizando Transformadas de Laplace.

Paso 1: Sistema de ecuaciones diferenciales
De las ecuaciones del circuito, tenemos:
dq q
R1 + =E(t)
dt C
2
d q dq q
L 2
+ R2 − =0
dt dt C
Reescribimos la ecuación (1) en términos de q(t):
dq q
=i 1 ,⇒ R 1 i 1 + =E(t )
dt C
De la ecuación (2):
2
d q dq q
L 2
+ R2 − =0
dt dt C

Paso 2: Transformada de Laplace

Aplicamos la Transformada de Laplace a las ecuaciones
Propiedades:

L
[ dq
dt
=sQ ( s )−q (0)
]
L
[ d2q 2
dt 2
dq
=s Q ( s )−s q ( 0 )− (0)
dt ]
Ecuación 1:

R1 L
[ dq 1
+ L [ q (t) ]=L [ E( t) ]
dt C ]
Sustituyendo:
1
R1 ( sQ ( s )−q ( 0 ) )+ Q ( s )=E (s )
C
Como q (0) = 0 queda
1
R1 sQ ( s ) + Q ( s )=E(s)
C
Ecuación 2:

Sustituyendo:

( 2
L s Q ( s )−sq ( 0 ) −
dq
dt )
( 0 ) + R2 ¿

dq
Con q ( 0 )=0 y ( 0 )=0queda
dt
2 1
L s Q ( s ) + R2 s Q ( s ) − Q ( s )=0
C

Paso 3: Resolviendo Q(s):

Del primera sistema:

Del segundo sistema:

Sustituimos E(s) en el dominio de Laplace. E(t) tiene la forma:

Su transformada es:
120
E ( s )=L [ 120 sin ⁡( t) ] = 2
, para 1 ≤ t<5
s +1
Sustituyendo:

Paso 4: Transformada inversa para obtener q(t)

Para simplificar la anti transformada, descomponemos en fracciones parciales:
120 A B
= +
1 s +1 R s+ 1
2
( s2 +1)(R1 s+ ) 1
C C
Resolviendo para A y B, integramos y usamos tablas de transformadas inversas.
Finalmente:

q ( t )=L
−1
[ Q(s)]
Realizamos la descomposición:

Expresión:
120
Q ( s )=
1
(s 2+1)(R1 s + )
C
 Planteamos la descomposición:
A B
Q(s)= +
s +1 R s + 1
2
1
C

Paso 2: Multiplicación por denominador común

2 1
Multiplicamos ambos lados por el denominador completo (s +1)(R1 s+ )
C

(
120=A R1 s +
1
C )
+ B( s 2+1)

Expandiendo:
A 2
120=A R 1 s+ +B s + B
C

Paso 3: Comparación de coeficientes

Comparando los coeficientes de cada potencia de s obtenemos:
1. Coeficiente de s2:
B=0.
2. Coeficiente de s:
120
A R1=120 ⟹ A=
R1
3. Coeficiente constante:
A
+ B=120
C
120
4. Sustituyendo A= y B=0
R1
120
R1
=120
C
Resolviendo:
120
=120 ⇒ R1 C=1
R1
Paso 4: Solución de fracciones parciales
La ecuación queda:
 120
R1
Q ( s )= 2
s +1
Transformada inversa:

[ ]
120
R1
q ( t )=L−1
s 2 +1

Sabemos que:

L−1
[ ]
2
1
s +1
=sin ⁡(t )

Entonces la solución final para q(t):

120
q ( t )= sin ⁡(t)
R1

Implementación de código para resolver el problema del circuito RLC en paralelo

Código:

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# Parámetros del circuito

L = 1 # Henrios
R1 = 5 # Ohmios
R2 = 10 # Ohmios
C = 0.01 # Faradios

# Función E(t)
def E(t):
if 0 <= t < 1:
return 0
elif 1 <= t < 5:
return 120 * np.sin(t)
else:
 return 0

# Sistema de ecuaciones diferenciales

def circuito_RLC(t, y):
q, dq = y # y[0] = q(t), y[1] = dq/dt
i1 = dq
i3 = (q / C - R2 * i1) / L
dqdt = (E(t) - R1 * i1 - q / C) / L
return [dq, dqdt]

# Solución numérica
t_span = (0, 10) # Intervalo de tiempo
t_eval = np.linspace(0, 10, 1000) # Puntos de evaluación
y0 = [0, 0] # Condiciones iniciales q(0) = 0, dq/dt(0) = 0

# Resolver el sistema de ecuaciones diferenciales

sol = solve_ivp(circuito_RLC, t_span, y0, t_eval=t_eval, args=())

# Extraer solución
t = sol.t
q = sol.y[0] # Carga q(t)
dqdt = sol.y[1] # Derivada de la carga dq/dt
i3 = dqdt # i3(t) es equivalente a dq/dt

# Generar gráfica para i3(t)

plt.figure(figsize=(10, 6))
plt.plot(t, i3, label='$i_3(t)$', linewidth=2)
plt.title('Corriente $i_3(t)$ en el circuito RLC', fontsize=14)
plt.xlabel('Tiempo (s)', fontsize=12)
plt.ylabel('Corriente $i_3$ (A)', fontsize=12)
plt.grid(True)
plt.legend()
plt.show()
 # Construir tabla para los primeros 5 segundos
indices_5s = np.where(t <= 5)[0]
tabla_5s = {
"Tiempo (s)": t[indices_5s],
"Carga q(t) (C)": q[indices_5s],
"Corriente i3(t) (A)": i3[indices_5s]
}

import pandas as pd

# Crear tabla y mostrarla

tabla_5s_df = pd.DataFrame(tabla_5s)
print(tabla_5s_df.head()) # Mostrar primeras filas de la tabla

# Opcional: Exportar tabla a un archivo CSV

tabla_5s_df.to_csv("resultados_rlc.csv", index=False)

Explicación:
1. Definición de parámetros:
 Se definen los valores para L, R1, R2, y C.
 Se define la función E(t) como una pieza por partes.
2. Sistema de ecuaciones:
 Se define la función circuito_RLC para representar el sistema de
ecuaciones diferenciales.
3. Resolución numérica:
 Se usa solve_ivp para integrar el sistema de ecuaciones en el
intervalo t=[0,10].
4. Gráfica:
 Se genera una gráfica para la corriente i3(t) en función del tiempo.
5. Tabla de resultados:
 Se crea una tabla con los valores de q(t) e i3(t) en los primeros 5
segundos y se exporta opcionalmente como archivo CSV
 4. Conclusiones
Conclusión de Gaviña Vázquez Byron Eduardo
El análisis del circuito RLC en paralelo permitió comprender el comportamiento dinámico de la
carga y la corriente en sus componentes. Las ecuaciones diferenciales derivadas de las leyes de
Kirchhoff modelaron correctamente el sistema, proporcionando una descripción matemática de
su evolución temporal.
Conclusión de César Eduardo Martínez Arredondo

5.Referencias Bibliográficas

 Dorf, R. C., & Svoboda, J. A. (2010). Introduction to Electric Circuits. Wiley.

 Kreyszig, E. (2011). Advanced Engineering Mathematics. Wiley.
 Franklin, G. F., Powell, J. D., & Emami-Naeini, A. (2015). Feedback Control of Dynamic Systems. Pearson.