Soluciones a Problemas de Ecuaciones Diferenciales
Soluciones a Problemas de Ecuaciones Diferenciales
Lista 1
1.1. Problema 1
**Solución al Problema de Valor Inicial (P.V.I.)**
Dado el sistema:
dP
= (10−1 − 10−7 P )P
dt
P (0) = 5000
dP P
= rP 1−
dt K
Al comparar con la ecuación dada:
dP P
= 10−1 P − 10−7 P 2 = 10−1 P 1− 6
dt 10
Identificamos: - Tasa de crecimiento r = 10−1 = 0,1 (mes−1 ) - Capacidad de carga
K = 106 = 1, 000, 000
El valor lı́mite corresponde a la capacidad de carga:
1
La población inicial es P (0) = 5,000, y buscamos t tal que P (t) = 2,500.
La solución general de la ecuación logı́stica es:
K
P (t) =
K−P0
1+ P0
e−rt
Sustituyendo los valores:
1,000,000 1,000,000
P (t) = =
1+ 1,000,000−5,000
e−0,1t 1 + 199e−0,1t
5,000
Igualando a 2,500:
1,000,000
2,500 =
1 + 199e−0,1t
Resolviendo paso a paso:
1. Invertir la ecuación:
1,000,000
1 + 199e−0,1t = = 400
2,500
2. Simplificar:
399
199e−0,1t = 399 =⇒ e−0,1t = ≈ 2,005025
199
3. Aplicar logaritmo natural:
4. Despejar t:
0,6955
t≈ ≈ −6,955 meses
−0,1
**Respuesta final:**
t ≈ −6,96 meses
2
1.2. Problema 2
Para resolver la ecuación diferencial ordinaria (EDO):
donde:
x si 0 ≤ x < 1
f (x) =
−x si x ≥ 1
2x f (x)
y′ + 2
y=
1+x 1 + x2
**Paso 2:** Calcular el factor integrante µ(x):
2x
R
dx 2
µ(x) = e 1+x2 = eln(1+x ) = 1 + x2
d
(1 + x2 )y ′ + 2xy = f (x) =⇒ (1 + x2 )y = f (x)
dx
**Paso 4:** Integrar ambos lados:
Z
2
(1 + x )y = f (x) dx + C
R
f (x) dx + C
y(x) =
1 + x2
R
**Paso 6:** Calcular f (x) dx en los intervalos dados:
- Para 0 ≤ x < 1 (f (x) = x):
x2
Z
x dx = + C1
2
- Para x ≥ 1 (f (x) = −x):
x2
Z
−x dx = − + C2
2
3
**Paso 7:** Aplicar la condición inicial y(0) = 0:
0+C
y(0) = = C = 0 =⇒ C = 0
1+0
**Paso 8:** Asegurar continuidad en x = 1:
- Para x → 1− :
1
1
y(1− ) = 2
=
2 4
- Para x → 1+ :
+− 12 + C2
y(1 ) =
2
Igualando y(1− ) = y(1+ ):
1 − 1 + C2
= 2 =⇒ C2 = 1
4 2
**Solución final:**
x2
si 0 ≤ x < 1
2(1+x2 )
y(x) =
2−x2
si x ≥ 1
2(1+x2 )
1.3. Problema 3
**Solución del Problema de Valor Inicial (P.V.I.)**
Dada la ecuación diferencial:
xy ′ + y = y 2 ln(x)
y(e) = 2
1 y 2 ln(x)
y′ + y =
x x
Esta es una ecuación de Bernoulli de la forma:
y ′ + P (x)y = Q(x)y n
4
ln(x)
donde P (x) = x1 , Q(x) = x
, y n = 2.
**Paso 2: Realizar la sustitución para convertirla en lineal**
Hacemos la sustitución v = y 1−n = y −1 . Entonces:
v ′ = −y −2 y ′
1 ln(x)
−y −2 y ′ − y −1 = −
x x
Sustituyendo v y v ′ :
1 ln(x)
v′ − v = −
x x
**Paso 3: Resolver la ecuación lineal**
Esta es ahora una ecuación lineal en v de la forma:
′1 ln(x)
v + − v=−
x x
El factor integrante µ(x) es:
R
− x1 dx 1
µ(x) = e = e− ln(x) =
x
Multiplicamos la ecuación por µ(x):
1 ′ 1 ln(x)
v − 2v = − 2
x x x
El lado izquierdo es la derivada de xv :
d v ln(x)
=− 2
dx x x
Integramos ambos lados:
Z
v ln(x)
=− dx
x x2
1
Para resolver la integral, usamos integración por partes con u = ln(x) y dv = x2
dx:
Z Z
ln(x) ln(x) 1 ln(x) 1
2
dx = − + 2
dx = − − +C
x x x x x
5
Por lo tanto:
v ln(x) 1
= + +C
x x x
Despejando v:
v = ln(x) + 1 + Cx
1
= ln(x) + 1 + Cx
y
Aplicamos la condición inicial y(e) = 2:
1 1 1
= ln(e) + 1 + Ce =⇒ = 1 + 1 + Ce =⇒ = 2 + Ce
2 2 2
Resolviendo para C:
1 3 3
Ce = − 2 = − =⇒ C = −
2 2 2e
**Paso 5: Escribir la solución final**
Sustituyendo C en la ecuación:
1 3x
= ln(x) + 1 −
y 2e
Despejando y:
1
y(x) = 3x
ln(x) + 1 − 2e
**Respuesta final:**
1
y(x) =
3x
ln(x) + 1 −
2e
6
1.4. Problema 4
**Solución de la Ecuación Diferencial Ordinaria**
Dada la ecuación diferencial:
∂M ∂
= [−(2y − x + 5)] = −2
∂y ∂y
∂N ∂
= [2x − y − 4] = 2
∂x ∂x
∂M ∂N
Como ∂y
̸= ∂x
, la ecuación no es exacta.
**Paso 3: Buscar un factor integrante**
Calculamos:
∂M ∂N
∂y
− ∂x −2 − 2 −4
= =
N 2x − y − 4 2x − y − 4
Este resultado no es función solo de x o solo de y, por lo que probamos con un factor
integrante de la forma µ(x + y).
**Paso 4: Realizar cambio de variable**
Hacemos el cambio u = x + y, y buscamos un factor integrante µ(u).
Multiplicamos la ecuación original por µ(u):
7
∂ ∂
[µ(u)(2x − y − 4)] = [−µ(u)(2y − x + 5)]
∂x ∂y
Calculamos las derivadas parciales:
∂
[µ(u)(2x − y − 4)] = µ′ (u)(2x − y − 4) + µ(u) · 2
∂x
∂
[−µ(u)(2y − x + 5)] = −µ′ (u)(2y − x + 5) − µ(u) · 2
∂y
Igualando ambas expresiones:
µ′ (u)(2x − y − 4 + 2y − x + 5) + 4µ(u) = 0
Simplificando:
µ′ (u)(x + y + 1) + 4µ(u) = 0
Como u = x + y, tenemos:
µ′ (u)(u + 1) + 4µ(u) = 0
dµ 4
=− du
µ u+1
Integrando:
C
ln |µ| = −4 ln |u + 1| + C =⇒ µ(u) =
(u + 1)4
Tomando C = 1, el factor integrante es:
1
µ(u) =
(x + y + 1)4
**Paso 5: Multiplicar por el factor integrante y resolver**
Multiplicamos la ecuación original por µ(u):
2x − y − 4 2y − x + 5
dy − dx = 0
(x + y + 1)4 (x + y + 1)4
Ahora verificamos que la ecuación es exacta. Calculamos:
8
2(x + y + 1)4 − 4(2x − y − 4)(x + y + 1)3
∂ 2x − y − 4
=
∂x (x + y + 1)4 (x + y + 1)8
∂F 2y − x + 5 ∂F 2x − y − 4
=− , =
∂x (x + y + 1)4 ∂y (x + y + 1)4
∂F
Integramos ∂y
con respecto a y:
2x − y − 4
Z
F (x, y) = dy
(x + y + 1)4
Hacemos el cambio de variable v = x + y + 1, dv = dy:
2x − (v − x − 1) − 4 3x − v + 1 − 4 3x − v − 3
Z Z Z
F (x, y) = dv = dv = dv
v4 v4 v4
Separamos la integral:
v −3 v −2
Z Z
−4 −3
F (x, y) = (3x − 3) v dv − v dv = (3x − 3) − + g(x)
−3 −2
Simplificando:
3x − 3 1 x−1 1
F (x, y) = − 3
+ 2 + g(x) = − 3 + 2 + g(x)
3v 2v v 2v
∂F
Derivamos F con respecto a x e igualamos a ∂x
:
∂F 1 3(x − 1) 2y − x + 5
=− 3 + 4
+ g ′ (x) = −
∂x v v v4
Sustituyendo v = x + y + 1 y simplificando, obtenemos:
g ′ (x) = 0 =⇒ g(x) = C
9
x−1 1
− 3
+ =C
(x + y + 1) 2(x + y + 1)2
**Respuesta final:**
x−1 1
− + =C
(x + y + 1)3 2(x + y + 1)2
1.5. Problema 5
**Demostración del lı́mite para la solución de la EDO lineal**
Consideremos la ecuación diferencial ordinaria lineal de primer orden:
Z x
−αx αt
Φ(x) = e b(t)e dt + C
0
β
lı́m Φ(x) =
x→+∞ α
**Paso 3: Reescribir la solución**
Dividimos la solución en dos partes:
Z x
−αx
Φ(x) = e b(t)eαt dt + Ce−αx
0
Observamos que:
10
Por lo tanto, el comportamiento de Φ(x) cuando x → +∞ está dominado por el primer
término.
**Paso 4: Aplicar la regla de L’Hôpital**
Consideramos el lı́mite:
Z x
−αx
lı́m e b(t)eαt dt
x→+∞ 0
∞
Este es un lı́mite del tipo ∞
, por lo que aplicamos la regla de L’Hôpital:
Rx
0
b(t)eαt dt b(x)eαx b(x) β
lı́m αx
= lı́m αx
= lı́m =
x→+∞ e x→+∞ αe x→+∞ α α
**Paso 5: Conclusión**
Por lo tanto, combinando los resultados:
β β
lı́m Φ(x) = +0=
x→+∞ α α
**Respuesta final:**
β
lı́m Φ(x) =
x→+∞ α
11
Capı́tulo 2
Lista 2
2.1. Problema 1
**Solución del circuito RLC en serie**
Dado el circuito RLC con los siguientes parámetros: - Resistencia R = 0,06 Ω - Induc-
50
tancia L = 0,01 H - Capacitancia C = 89
F - Fuente de voltaje E(t) = 0,1 sin(10t) V -
Condiciones iniciales: q(0) = 0 C, i(0) = 0 A
**Paso 1: Ecuación diferencial del circuito**
La ecuación diferencial para un circuito RLC en serie es:
d2 i di 1 dE
L 2 +R + i=
dt dt C dt
Sustituyendo los valores y derivando E(t):
d2 i di 89 d
0,01 + 0,06 + i = [0,1 sin(10t)] = cos(10t)
dt2 dt 50 dt
Multiplicando por 100 para simplificar:
d2 i di
2
+ 6 + 178i = 100 cos(10t)
dt dt
**Paso 2: Solución de la ecuación homogénea**
La ecuación caracterı́stica es:
r2 + 6r + 178 = 0
Raı́ces:
12
√ √
−6 ± 36 − 712 −6 ± −676
r= = = −3 ± 13i
2 2
Solución homogénea:
Derivadas:
Agrupando términos:
(−100A + 60B + 178A) cos(10t) + (−100B − 60A + 178B) sin(10t) = 100 cos(10t)
Simplificando:
Igualando coeficientes:
−60A + 78B = 0
Resolviendo el sistema:
13
78 60
A= ≈ 0,78, B= ≈ 0,60
100 100
**Paso 4: Solución general**
**Solución final:**
i(t) = e−3t (−0,78 cos(13t) − 0,64 sin(13t)) + 0,78 cos(10t) + 0,60 sin(10t)
2.2. Problema 2
**Solución del Problema de Valores en la Frontera (PVF)**
Dado el sistema:
y ′′ + βy = sin t
y(0) = 0,
y(π) = 0
con β ≥ 0 y β ̸= 1.
**Paso 1: Solución general de la ecuación homogénea**
La ecuación homogénea asociada es:
y ′′ + βy = 0
14
p p
yh (t) = C1 cos( βt) + C2 sin( βt)
Simplificando:
Igualando coeficientes:
A(β − 1) = 1
B(β − 1) = 0
Como β ̸= 1, obtenemos:
1
A= , B=0
β−1
Por lo tanto, la solución particular es:
sin t
yp (t) =
β−1
**Paso 3: Solución general de la ecuación completa**
La solución general es la suma de la solución homogénea y particular:
p p sin t
y(t) = C1 cos( βt) + C2 sin( βt) +
β−1
**Paso 4: Aplicación de las condiciones de frontera**
1. Para y(0) = 0:
sin 0
C1 cos(0) + C2 sin(0) + = C1 = 0
β−1
15
2. Para y(π) = 0:
p p sin π p
C1 cos( βπ) + C2 sin( βπ) + = C2 sin( βπ) = 0
β−1
**Paso 5: Análisis de casos**
√
- **Caso 1:** Si sin( βπ) ̸= 0, entonces C2 = 0, y la solución es trivial:
sin t
y(t) =
β−1
√ √
- **Caso 2:** Si sin( βπ) = 0, entonces β = n con n ∈ Z+ . Es decir:
β = n2 para algún n ≥ 1
sin t
y(t) = C2 sin(nt) +
β−1
**Conclusión:**
El PVF admite soluciones no triviales cuando β = n2 para algún entero n ≥ 1 (es
decir, cuando β es el cuadrado de un entero positivo). Para estos valores, la solución es:
sin t
y(t) = C sin(nt) +
n2 − 1
donde C es una constante arbitraria no nula.
**Respuesta final:**
Valores de β : β = n2 con n ∈ N, n ≥ 1
sin t
Solución no trivial: y(t) = C sin(nt) + (C ̸= 0)
n2 − 1
2.3. Problema 3
**Solución general de la ecuación diferencial**
Dada la ecuación diferencial no homogénea:
y ′′ − 4y ′ + 3y = sin3 (x)
16
Resolvemos primero la ecuación homogénea:
y ′′ − 4y ′ + 3y = 0
√
2 4± 16 − 12 4±2
r − 4r + 3 = 0 =⇒ r = = =⇒ r1 = 3, r2 = 1
2 2
Por lo tanto, la solución general de la ecuación homogénea es:
yh (x) = C1 e3x + C2 ex
3 sin(x) − sin(3x)
sin3 (x) =
4
Ası́, la ecuación original se convierte en:
3 sin(x) − sin(3x)
y ′′ − 4y ′ + 3y =
4
**Paso 3: Solución particular para cada componente**
Buscamos una solución particular yp (x) para la ecuación no homogénea. Consideramos
los dos términos por separado:
3
1. **Para 4
sin(x):**
Suponemos una solución de la forma:
′
yp1 (x) = A cos(x) − B sin(x)
′′
yp1 (x) = −A sin(x) − B cos(x)
3
(−A sin(x) − B cos(x)) − 4(A cos(x) − B sin(x)) + 3(A sin(x) + B cos(x)) = sin(x)
4
17
Simplificando:
3
(−A + 4B + 3A) sin(x) + (−B − 4A + 3B) cos(x) = sin(x)
4
Igualando coeficientes:
3
2A + 4B = y − 4A + 2B = 0
4
Resolviendo el sistema:
3 3 3 3
B = 2A =⇒ 2A + 4(2A) = =⇒ 10A = =⇒ A = , B=
4 4 40 20
Por lo tanto:
3 3
yp1 (x) = sin(x) + cos(x)
40 20
2. **Para − 14 sin(3x):**
Suponemos una solución de la forma:
′
yp2 (x) = 3C cos(3x) − 3D sin(3x)
′′
yp2 (x) = −9C sin(3x) − 9D cos(3x)
1
(−9C sin(3x)−9D cos(3x))−4(3C cos(3x)−3D sin(3x))+3(C sin(3x)+D cos(3x)) = − sin(3x)
4
Simplificando:
1
(−9C + 12D + 3C) sin(3x) + (−9D − 12C + 3D) cos(3x) = − sin(3x)
4
Igualando coeficientes:
1
−6C + 12D = − y − 12C − 6D = 0
4
Resolviendo el sistema:
1 1 1 1
D = −2C =⇒ −6C + 12(−2C) = − =⇒ −30C = − =⇒ C = , D=−
4 4 120 60
18
Por lo tanto:
1 1
yp2 (x) = sin(3x) − cos(3x)
120 60
**Paso 4: Solución particular total**
Sumamos las soluciones particulares obtenidas:
3 3 1 1
yp (x) = yp1 (x) + yp2 (x) = sin(x) + cos(x) + sin(3x) − cos(3x)
40 20 120 60
3 3 1 1
y(x) = yh (x) + yp (x) = C1 e3x + C2 ex + sin(x) + cos(x) + sin(3x) − cos(3x)
40 20 120 60
**Respuesta final:**
3 3 1 1
y(x) = C1 e3x + C2 ex + sin(x) + cos(x) + sin(3x) − cos(3x)
40 20 120 60
2.4. Problema 4
**Solución general de la ecuación diferencial**
Dada la ecuación diferencial lineal no homogénea de tercer orden:
y ′′′ + 2y ′′ + y ′ = 4et
y ′′′ + 2y ′′ + y ′ = 0
19
Por lo tanto, la solución general de la ecuación homogénea es:
y ′′′ + 2y ′′ + y ′ = 4et
yp (t) = Aet
Igualando coeficientes:
4A = 4 =⇒ A = 1
yp (t) = et
**Respuesta final:**
20
2.5. Problema 5
**Demostración de la fórmula para el Wronskiano**
Sean ϕ1 y ϕ2 dos soluciones de la ecuación diferencial homogénea de segundo orden:
y ′′ + a1 (x)y ′ + a2 (x)y = 0
Rx
− a1 (t)dt
W [ϕ1 , ϕ2 ](x) = e x0
W [ϕ1 , ϕ2 ](x0 )
ϕ1 (x) ϕ2 (x)
W [ϕ1 , ϕ2 ](x) = = ϕ1 (x)ϕ′2 (x) − ϕ2 (x)ϕ′1 (x)
ϕ′1 (x) ϕ′2 (x)
**Paso 2: Derivada del Wronskiano**
Calculamos la derivada del Wronskiano con respecto a x:
d
W (x) = ϕ′1 (x)ϕ′2 (x) + ϕ1 (x)ϕ′′2 (x) − ϕ′2 (x)ϕ′1 (x) − ϕ2 (x)ϕ′′1 (x) = ϕ1 (x)ϕ′′2 (x) − ϕ2 (x)ϕ′′1 (x)
dx
ϕ′′1 (x) = −a1 (x)ϕ′1 (x) − a2 (x)ϕ1 (x)ϕ′′2 (x) = −a1 (x)ϕ′2 (x) − a2 (x)ϕ2 (x)
d
Sustituyendo en la expresión para dx
W (x):
d
W (x) = ϕ1 (x) (−a1 (x)ϕ′2 (x) − a2 (x)ϕ2 (x)) − ϕ2 (x) (−a1 (x)ϕ′1 (x) − a2 (x)ϕ1 (x))
dx
Simplificando:
d
W (x) = −a1 (x) (ϕ1 (x)ϕ′2 (x) − ϕ2 (x)ϕ′1 (x)) = −a1 (x)W (x)
dx
**Paso 4: Resolución de la ecuación diferencial para W (x)**
21
Hemos obtenido:
dW
= −a1 (x)W
dx
Esta es una ecuación diferencial separable. Resolviéndola:
Z
dW R
= −a1 (x)dx ln |W | = − a1 (x)dx + CW (x) = Ce− a1 (x)dx
W
Evaluando en x = x0 :
Rx
− 0 a1 (t)dt
W (x0 ) = Ce x0
=C
Por lo tanto:
Rx
− a1 (t)dt
W (x) = W (x0 )e x0
**Respuesta final:**
Rx
− a1 (t)dt
W [ϕ1 , ϕ2 ](x) = e x0
W [ϕ1 , ϕ2 ](x0 )
22
Capı́tulo 3
Lista 3
3.1. Problema 1
**Solución del sistema de ecuaciones diferenciales y cálculo de eAt **
*** Parte 1: Resolución del sistema X ′ (t) = AX(t)
Dado el sistema:
−2 0 1
′
X (t) = −6 −2 0 X(t)
19 5 −4
**Paso 1: Encontrar los valores propios de A**
Calculamos el determinante de A − λI:
−2 − λ 0 1
−2 − λ 0 −6 −2 − λ
det(A − λI) = −6 −2 − λ 0 = (−2 − λ) −1
5 −4 − λ 19 5
19 5 −4 − λ
λ1 = −1, λ2 = −3, λ3 = −4
23
Para λ1 = −1:
−1 0 1 x
(A + I)v1 = 0 =⇒ −6 −1 0 y = 0
19 5 −3 z
1
Solución: v1 = −6
1
Para λ2 = −3:
1 0 1 x
(A + 3I)v2 = 0 =⇒ −6 1 0 y = 0
19 5 −1 z
1
Solución: v2 = 6
−1
Para λ3 = −4:
2 0 1 x
(A + 4I)v3 = 0 =⇒ −6 2 0 y = 0
19 5 0 z
−1
Solución: v3 = 3
2
**Paso 3: Solución general del sistema**
1 1 −1
X(t) = c1 e−t −6 + c2 e−3t 6 + c3 e−4t 3
1 −1 2
3.2. Problema 2
*** Parte 2: Cálculo de eAt para la matriz dada
Dada la matriz:
−4 1 1
A = 1 5 −1
0 1 −3
24
**Paso 1: Encontrar los valores propios de A**
−4 − λ 1 1
5−λ −1 1 −1 1 5−λ
det(A−λI) = 1 5−λ −1 = (−4−λ) −1 +1
1 −3 − λ 0 −3 − λ 0 1
0 1 −3 − λ
= −λ3 − 2λ2 + 6λ + 56 + λ + 4
= −λ3 − 2λ2 + 7λ + 60 = 0
Valores propios:
λ1 = −3, λ2 = 4, λ3 = 5
2. Exponencial de la matriz A:
e−3t 0 0
eAt = P 0 e4t 0 P −1
5t
0 0 e
25
3.3. Problema 3
**Solución del sistema de ecuaciones diferenciales no homogéneo**
Dado el sistema:
t
1 1 0 e
Y ′ (t) = AY (t) + F (t),
2t
con A = 1 1 0 , F (t) = e
0 0 3 te3t
1−λ 1 0
= (3 − λ) (1 − λ)2 − 1 = (3 − λ)(λ2 − 2λ) = 0
det(A − λI) = 1 1−λ 0
0 0 3−λ
Raı́ces: λ1 = 0, λ2 = 2, λ3 = 3.
**Vectores propios:** - Para λ1 = 0:
1 1 0 x 1
0 y = 0 =⇒ v1 = −1
1 1
0 0 3 z 0
- Para λ2 = 2:
−1 1 0 x 1
1 −1 0 y = 0 =⇒ v2 = 1
0 0 1 z 0
- Para λ3 = 3:
−2 1 0 x 0
1 −2 0 y = 0 =⇒ v3 = 0
0 0 0 z 1
**Solución homogénea:**
1 1 0
Yh (t) = c1 −1 + c2 e2t 1 + c3 e3t 0
0 0 1
26
t
ae
1
1. **Para et (primer componente):** Proponemos Yp1
t
= a2 e . Sustituyendo:
0
a1 et = (a1 + a2 )et + et
a1 = a1 + a2 + 1
a2 et = (a1 + a2 )et =⇒ =⇒ a1 = −1, a2 = 0
a2 = a1 + a2
0 = 0 + 0
t
−e
Yp1 = 0
0
b1 e2t
2. **Para e2t (segundo componente):** Proponemos Yp2 = b2 e2t . Sustituyendo:
0
2b1 e2t = (b1 + b2 )e2t + e2t
2b1 = b1 + b2 + 1
2b2 e2t = (b1 + b2 )e2t =⇒ =⇒ b1 = 1, b2 = 1
2b2 = b1 + b2
0 = 0 + 0
2t
e
2t
Yp2 = e
0
0
3. **Para te3t (tercer componente):** Proponemos Yp3 = . Sustituyen-
0
3t
(c1 t + c2 )e
do:
(3c1 t + 3c2 + c1 )e3t = 3(c1 t + c2 )e3t + te3t
Igualando coeficientes:
3c1 = 3c1 + 1 (No tiene solución)
Ajustamos la propuesta:
0
Yp3 =
0
(c1 t2 + c2 t)e3t
27
Derivando y sustituyendo:
1
2c1 = 1 =⇒ c1 = , c2 = 0
2
0
Yp3 = 0
1 2 3t
2
te
**Solución particular total:**
t 2t
−e + e
Yp (t) = e2t
1 2 3t
2
te
**Respuesta final:**
1 1 0 −et + e2t
Y (t) = c1 −1 + c2 e2t 1 + c3 e3t 0 + e2t
1 2 3t
0 0 1 2
te
3.4. Problema 4
**Solución del circuito eléctrico R-L-C**
Dado el circuito con los siguientes parámetros: - Resistencias: R1 = 8 Ω, R2 = 3 Ω
- Inductancias: L1 = L2 = 1 H - Fuente de voltaje: E(t) = sin(100t) V - Condiciones
iniciales: i1 (0) = 0 A, i2 (0) = 0 A
**Paso 1: Modelado del sistema**
Aplicando las leyes de Kirchhoff, obtenemos el siguiente sistema de ecuaciones dife-
renciales:
28
L1 di1 + R1 i1 + R2 (i1 − i2 ) = E(t)
dt
−11 − λ 3
det(A − λI) = = λ2 + 14λ + 24 = 0
3 −3 − λ
Raı́ces:
λ1 = −2, λ2 = −12
Vectores propios:
Para λ1 = −2:
−9 3 x 1
= 0 =⇒ v1 =
3 −1 y 3
Para λ2 = −12:
1 3 x 3
= 0 =⇒ v2 =
3 9 y −1
Solución homogénea:
i (t) 1 3
1h = c1 e−2t + c2 e−12t
i2h (t) 3 −1
**Paso 3: Solución particular**
29
Proponemos una solución particular de la forma:
i (t) A cos(100t) + B sin(100t)
1p =
i2p (t) C cos(100t) + D sin(100t)
Sustituyendo en el sistema y resolviendo, obtenemos:
300 1100
i (t) 1 3 − cos(100t) + 10009 sin(100t)
1 = c1 e−2t + c2 e−12t + 10009
300 3300
i2 (t) 3 −1 10009
cos(100t) + 10009
sin(100t)
0 = 3c1 − c2 +
300
10009
Resolviendo:
100
c1 = 0, c2 =
10009
**Solución final:**
3.5. Problema 5
**Demostración de la unicidad de la solución para un sistema diferencial lineal con
matriz diagonalizable**
**Teorema:** Sea A una matriz de orden n × n diagonalizable. Entonces, para todo
X0 ∈ Rn , el problema de valor inicial:
30
tiene una **solución única**.
—
**Demostración:**
**Paso 1: Diagonalización de A** Como A es diagonalizable, existe una matriz inver-
tible P y una matriz diagonal D tal que:
A = P DP −1 ,
X(t) = P Y (t).
Y ′ (t) = DY (t).
yi (t) = eλi t (P −1 X0 )i .
En forma matricial:
Y (t) = eDt P −1 X0 ,
31
**Paso 5: Unicidad de la solución** Dado que: 1. La construcción de eAt es única (pues
P y D están únicamente determinadas por A). 2. La solución X(t) = eAt X0 satisface tanto
la ecuación diferencial como la condición inicial X(0) = X0 . 3. Por el teorema de existencia
y unicidad para sistemas lineales, esta es la **única solución**.
—
**Conclusión:** Para toda matriz diagonalizable A y cualquier condición inicial X0 ,
el PVI:
X ′ (t) = AX(t), X(0) = X0
32
Capı́tulo 4
Lista 4
4.1. Problema 1
**Solución para determinar la transformada de Laplace de g(t)**
Dada la función:
e−2t
−2e2t cos(5t) + e4t + 1
g(t) =
58
**Paso 1: Simplificar la expresión**
Expandimos el término dentro del corchete:
1 1 1
g(t) = − cos(5t) + e2t + e−2t
29 58 58
33
1 −2t
3. Para 58
e :
1 −2t 1 1
L e = ·
58 58 s + 2
**Paso 3: Combinar los resultados**
Sumamos las transformadas de cada término:
s 1 1
L{g(t)} = − + +
29(s2 + 25) 58(s − 2) 58(s + 2)
**Paso 4: Simplificar la expresión final**
Podemos escribir el resultado como:
1 1 1 s
L{g(t)} = + − 2
58 s − 2 s + 2 29(s + 25)
**Respuesta final:**
s 1 1
L{g(t)} = − + +
29(s2 + 25) 58(s − 2) 58(s + 2)
4.2. Problema 2
**Solución para calcular la transformada de Laplace de h(t)**
Dada la función definida por tramos:
π
2
para 0 ≤ t < π2 ,
π
π − t
para ≤ t < π,
2
h(t) =
π para π ≤ t < 2π,
t − 2π
para 2π ≤ t < 3π,
34
n π π o π π
L − u t− = − e− 2 s .
2 2 2s
2. **Para π − t en [ π2 , π):**
n π o − π2 s π − π2 s π 1
L (π − t)u t − =e L{π − (t + )} = e − .
2 2 2s s2
e−πs
−πs −πs 0 1
L {−(π − t)u(t − π)} = −e L{π − (t + π)} = −e − = 2 .
s s2 s
3. **Para π en [π, 2π):**
π −πs
L {πu(t − π)} =e .
s
π
L {−πu(t − 2π)} = − e−2πs .
s
4. **Para t − 2π en [2π, 3π):**
e−2πs
L {(t − 2π)u(t − 2π)} = .
s2
e−3πs
L {−(t − 2π)u(t − 3π)} = −e−3πs L{(t + 2π − 2π)} = − .
s2
**Paso 3: Combinar todos los términos**
Sumamos todas las transformadas calculadas:
Simplificamos la expresión:
π
e−πs + e−2πs e− 2 s e−3πs
π π π π −πs
L{h(t)} = 1 − e− 2 s + e− 2 s + + e − e−2πs − 2 − 2 .
2s s2 s s s
**Nota:** La función h(t) es periódica con periodo 3π, por lo que también se puede
calcular usando la fórmula para funciones periódicas:
Z 3π
1
L{h(t)} = −3πs
h(t)e−st dt.
1−e 0
Sin embargo, el método presentado arriba es más directo para esta función especı́fica.
35
4.3. Problema 3
**Solución para determinar la transformada inversa de Laplace de F (s)**
Dada la función:
e−5s
F (s) =
(s2 + 4)(s2 − 9)
**Paso 1: Descomponer en fracciones parciales**
Primero, ignoramos el término exponencial e−5s y trabajamos con:
1
G(s) =
(s2 + 4)(s2 − 9)
1 As + B Cs + D
= + 2
(s2 + 4)(s2 − 9) s2 + 4 s −9
Igualando coeficientes:
A+C =0
B + D = 0
−9A + 4C = 0
−9B + 4D = 1
Resolviendo el sistema:
1 1
A = 0, B=− , C = 0, D=
13 13
Por lo tanto:
1 1 1 1
G(s) = − · 2 + · 2
13 s + 4 13 s − 9
**Paso 2: Aplicar la transformada inversa a G(s)**
Usamos las siguientes transformadas conocidas:
−1 1 sin(at) −1 1 sinh(at)
L 2 2
= , L 2 2
=
s +a a s −a a
36
Aplicando a G(s):
1 sin(2t) 1 sinh(3t)
g(t) = L−1 {G(s)} = − · + ·
13 2 13 3
sin(2t) sinh(3t)
g(t) = − +
26 39
**Paso 3: Incorporar el término exponencial e−5s **
La función original tiene un factor e−5s , que corresponde a un desplazamiento en el
tiempo:
F (s) = e−5s G(s) =⇒ f (t) = L−1 {F (s)} = g(t − 5)u(t − 5)
sin(2(t − 5)) sinh(3(t − 5))
f (t) = − + u(t − 5)
26 39
**Respuesta final:**
sinh(3(t − 5)) sin(2(t − 5))
f (t) = − u(t − 5)
39 26
4.4. Problema 4
**Solución del sistema de ecuaciones diferenciales usando transformada de Laplace**
Dado el circuito eléctrico con los siguientes parámetros: - Inductancias: L1 = 5 H,
L2 = 10 H - Resistencias: R1 = 30 Ω, R2 = 100 Ω - Fuente de voltaje: E(t) = 2 · u4 (t) V
(función escalón unitario en t = 4) - Condiciones iniciales: I1 (0) = 0 A, I2 (0) = 0 A
**Paso 1: Modelado del sistema**
Aplicando las leyes de Kirchhoff, obtenemos el siguiente sistema de ecuaciones dife-
renciales acopladas:
L1 dI1 + R1 (I1 − I2 ) = E(t)
dt
37
Simplificando:
dI1 + 6I1 − 6I2 = 2 u4 (t)
dt 5
Reorganizando términos:
(s + 6)L{I1 } − 6L{I2 } =
2e−4s
5s
−3L{I1 } + (s + 13)L{I2 } = 0
2e−4s 2e−4s
(s + 13) · 5s
3· 5s
L{I1 } = , L{I2 } =
∆ ∆
Simplificando:
38
Descomponemos L{I1 }:
s + 13 A B C
= + +
s(s + 4)(s + 15) s s + 4 s + 15
Resolviendo:
13 9 2
A= , B=− , C=
60 44 165
Análogamente para L{I2 }:
6 D E F
= + +
5s(s + 4)(s + 15) s s + 4 s + 15
Resolviendo:
1 3 2
D= , E=− , F =
50 110 275
**Paso 5: Aplicar la transformada inversa de Laplace**
Considerando el término e−4s , que indica un retardo en el tiempo:
2 13 9 −4(t−4) 2 −15(t−4)
I1 (t) = − e + e u(t − 4)
5 60 44 165
6 1 3 −4(t−4) 2 −15(t−4)
I2 (t) = − e + e u(t − 4)
5 50 110 275
**Respuesta final:**
13 9 −4(t−4) 4 −15(t−4)
I1 (t) = − e + e u(t − 4)
150 110 825
3 9 −4(t−4) 12 −15(t−4)
I2 (t) = − e + e u(t − 4)
125 275 1375
4.5. Problema 5
**Demostración de la propiedad de la transformada de Laplace para la integral de una
función**
Sea f : [0, +∞) → R una función seccionalmente continua y de orden exponencial c
en [0, +∞). Queremos demostrar que:
Z t
F (s)
L f (τ ) dτ =
0 s
39
donde F (s) = L{f (t)} es la transformada de Laplace de f .
—
**Demostración:**
**Paso 1: Definición de la transformada de Laplace para la integral**
Consideramos la función: Z t
g(t) = f (τ ) dτ
0
y g(0) = 0.
**Paso 3: Aplicar la transformada de Laplace a ambos lados**
Aplicamos la transformada de Laplace a la ecuación g ′ (t) = f (t):
sL{g(t)} = F (s)
F (s)
L{g(t)} =
s
Por lo tanto: Z t
F (s)
L f (τ ) dτ =
0 s
—
**Respuesta final:**
Z t
F (s)
L f (τ ) dτ =
0 s
40