0% encontró este documento útil (0 votos)
4 vistas40 páginas

Soluciones a Problemas de Ecuaciones Diferenciales

Cargado por

Sweet
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
4 vistas40 páginas

Soluciones a Problemas de Ecuaciones Diferenciales

Cargado por

Sweet
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Capı́tulo 1

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

donde t se mide en meses.


1. Valor lı́mite de la población**
Reconocemos la ecuación como del tipo logı́stico:

 
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:

lı́m P (t) = 1,000,000


t→∞

2. Tiempo para alcanzar la mitad de la población inicial**

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:

−0,1t = ln(2,005025) ≈ 0,6955

4. Despejar t:

0,6955
t≈ ≈ −6,955 meses
−0,1
**Respuesta final:**

t ≈ −6,96 meses

*Nota: El resultado se redondeó a dos decimales.*

2
1.2. Problema 2
Para resolver la ecuación diferencial ordinaria (EDO):

(1 + x2 )y ′ + 2xy = f (x), y(0) = 0

donde:

x si 0 ≤ x < 1

f (x) =
−x si x ≥ 1

**Paso 1:** Reescribir la EDO en forma estándar:

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

**Paso 3:** Multiplicar la EDO por µ(x) y simplificar:

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

**Paso 5:** Resolver para y(x):

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

**Paso 1: Reescribir la ecuación en forma estándar**


Dividimos ambos lados por x (para x > 0):

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 ′

Multiplicamos la ecuación original por −y −2 :

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

**Paso 4: Revertir la sustitución y aplicar la condición inicial**


Recordando que v = y1 :

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:

(2x − y − 4)dy − (2y − x + 5)dx = 0

**Paso 1: Reescribir la ecuación en forma estándar**


Expresamos la ecuación en la forma M (x, y)dx + N (x, y)dy = 0:

−(2y − x + 5)dx + (2x − y − 4)dy = 0

donde: - M (x, y) = −(2y − x + 5) - N (x, y) = 2x − y − 4


**Paso 2: Verificar si es exacta**
Calculamos las derivadas parciales:

∂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):

µ(u)(2x − y − 4)dy − µ(u)(2y − x + 5)dx = 0

Para que la nueva ecuación sea exacta, debe cumplirse:

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

Esta es una ecuación diferencial para µ(u):

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

−2(x + y + 1)4 + 4(2y − x + 5)(x + y + 1)3


 
∂ 2y − x + 5
− =
∂y (x + y + 1)4 (x + y + 1)8
Simplificando, se comprueba que las derivadas son iguales, por lo que la ecuación es
exacta.
**Paso 6: Encontrar la solución general**
Existe una función F (x, y) tal que:

∂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

Por lo tanto, la solución general es:

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:

y ′ + αy = b(x) con α > 0

donde b(x) es una función continua en [0, +∞) y satisface:

lı́m b(x) = β con β ∈ R


x→+∞

**Paso 1: Solución general de la EDO**


La solución general de la ecuación y ′ + αy = b(x) está dada por:

Z x 
−αx αt
Φ(x) = e b(t)e dt + C
0

donde C es una constante determinada por la condición inicial.


**Paso 2: Análisis del lı́mite cuando x → +∞**
Queremos demostrar que:

β
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:

lı́m Ce−αx = 0 ya que α > 0


x→+∞

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:

ih (t) = e−3t (C1 cos(13t) + C2 sin(13t))

**Paso 3: Solución particular**


Buscamos una solución particular de la forma:

ip (t) = A cos(10t) + B sin(10t)

Derivadas:

i′p (t) = −10A sin(10t) + 10B cos(10t)

i′′p (t) = −100A cos(10t) − 100B sin(10t)

Sustituyendo en la ecuación diferencial:

(−100A cos(10t)−100B sin(10t))+6(−10A sin(10t)+10B cos(10t))+178(A cos(10t)+B sin(10t)) = 100 co

Agrupando términos:

(−100A + 60B + 178A) cos(10t) + (−100B − 60A + 178B) sin(10t) = 100 cos(10t)

Simplificando:

(78A + 60B) cos(10t) + (−60A + 78B) sin(10t) = 100 cos(10t)

Igualando coeficientes:

78A + 60B = 100

−60A + 78B = 0

Resolviendo el sistema:

13
78 60
A= ≈ 0,78, B= ≈ 0,60
100 100
**Paso 4: Solución general**

i(t) = e−3t (C1 cos(13t) + C2 sin(13t)) + 0,78 cos(10t) + 0,60 sin(10t)

**Paso 5: Aplicación de condiciones iniciales**


Para t = 0, i(0) = 0:

0 = e0 (C1 cos(0) + C2 sin(0)) + 0,78 cos(0) + 0,60 sin(0) =⇒ C1 = −0,78

Derivando i(t) y evaluando en t = 0:

i′ (0) = −3C1 + 13C2 + 6 = 0 =⇒ 2,34 + 13C2 + 6 = 0 =⇒ C2 ≈ −0,64

**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

Cuya solución general es:

14
p p
yh (t) = C1 cos( βt) + C2 sin( βt)

**Paso 2: Solución particular de la ecuación no homogénea**


Buscamos una solución particular de la forma:

yp (t) = A sin t + B cos t

Sustituyendo en la ecuación original:

(−A sin t − B cos t) + β(A sin t + B cos t) = sin t

Simplificando:

(−A + βA) sin t + (−B + βB) cos t = sin t

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

En este caso, C2 puede ser cualquier constante, y la solución general es:

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)

**Paso 1: Solución de la ecuación homogénea asociada**

16
Resolvemos primero la ecuación homogénea:

y ′′ − 4y ′ + 3y = 0

La ecuación caracterı́stica es:


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

**Paso 2: Simplificación del término no homogéneo**


Utilizamos la identidad trigonométrica para expresar sin3 (x):

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 sin(x) + B cos(x)

Calculamos sus derivadas:


yp1 (x) = A cos(x) − B sin(x)

′′
yp1 (x) = −A sin(x) − B cos(x)

Sustituyendo en la ecuación diferencial:

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) = C sin(3x) + D cos(3x)

Calculamos sus derivadas:


yp2 (x) = 3C cos(3x) − 3D sin(3x)

′′
yp2 (x) = −9C sin(3x) − 9D cos(3x)

Sustituyendo en la ecuación diferencial:

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

**Paso 5: Solución general**


Finalmente, la solución general de la ecuación diferencial es la suma de la solución
homogénea y la solución particular:

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

**Paso 1: Solución de la ecuación homogénea asociada**


Resolvemos primero la ecuación homogénea:

y ′′′ + 2y ′′ + y ′ = 0

La ecuación caracterı́stica es:

r3 + 2r2 + r = 0 =⇒ r(r2 + 2r + 1) = 0 =⇒ r(r + 1)2 = 0

Las raı́ces son: - r1 = 0 (raı́z simple) - r2 = −1 (raı́z doble)

19
Por lo tanto, la solución general de la ecuación homogénea es:

yh (t) = C1 + C2 e−t + C3 te−t

**Paso 2: Solución particular de la ecuación no homogénea**


Buscamos una solución particular para:

y ′′′ + 2y ′′ + y ′ = 4et

Observamos que el término no homogéneo es 4et . Como et no es solución de la ecuación


homogénea, proponemos una solución particular de la forma:

yp (t) = Aet

Calculamos sus derivadas:

yp′ = Aet yp′′ = Aet yp′′′ = Aet

Sustituyendo en la ecuación diferencial:

Aet + 2Aet + Aet = 4et =⇒ 4Aet = 4et

Igualando coeficientes:
4A = 4 =⇒ A = 1

Por lo tanto, la solución particular es:

yp (t) = et

**Paso 3: Solución general**


La solución general de la ecuación diferencial completa es la suma de la solución
homogénea y la solución particular:

y(t) = yh (t) + yp (t) = C1 + C2 e−t + C3 te−t + et

**Respuesta final:**

y(t) = C1 + C2 e−t + C3 te−t + et

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

definidas en un intervalo I que contiene a x0 . Queremos demostrar que el Wronskiano


W [ϕ1 , ϕ2 ](x) satisface:

Rx
− a1 (t)dt
W [ϕ1 , ϕ2 ](x) = e x0
W [ϕ1 , ϕ2 ](x0 )

**Paso 1: Definición del Wronskiano**


El Wronskiano de ϕ1 y ϕ2 está dado por:

ϕ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

**Paso 3: Sustitución usando la ecuación diferencial**


Como ϕ1 y ϕ2 son soluciones de la ecuación diferencial, tenemos:

ϕ′′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 − λ

= (−2 − λ)[(−2 − λ)(−4 − λ) − 0] − 1[(−6)(5) − (−2 − λ)(19)]

= (−2 − λ)(λ2 + 6λ + 8) − (−30 + 19λ + 38)

= −λ3 − 8λ2 − 20λ − 16 + 8 − 19λ

= −λ3 − 8λ2 − 39λ − 8 = 0

Resolviendo la ecuación cúbica, encontramos los valores propios:

λ1 = −1, λ2 = −3, λ3 = −4

**Paso 2: Encontrar los vectores propios correspondientes**

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 − λ

= (−4 − λ)[(5 − λ)(−3 − λ) + 1] − 1[−3 − λ] + 1[1]

= (−4 − λ)(λ2 − 2λ − 14) + λ + 4

= −λ3 − 2λ2 + 6λ + 56 + λ + 4

= −λ3 − 2λ2 + 7λ + 60 = 0

Valores propios:
λ1 = −3, λ2 = 4, λ3 = 5

**Paso 2: Calcular eAt usando la forma canónica de Jordan**


Dado que A es diagonalizable (todos los valores propios son distintos), podemos ex-
presar:
eAt = P eDt P −1

donde D es la matriz diagonal de valores propios y P es la matriz de vectores propios.


**Respuestas finales:**
1. Solución del sistema:
     
1 1 −1
     
X(t) = c1 e−t −6 + c2 e−3t  6  + c3 e−4t  3 
     
     
1 −1 2

2. Exponencial de la matriz A:
 
e−3t 0 0
 
eAt = P  0 e4t 0  P −1
 
 
5t
0 0 e

*Donde P es la matriz cuyas columnas son los vectores propios de A.*

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

**Paso 1: Solución del sistema homogéneo Yh (t)**


Encontramos los valores propios de A:

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

**Paso 2: Solución particular Yp (t)**


Buscamos Yp (t) según la forma de F (t):

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:

(2c1 t + c2 + 3c1 t2 + 3c2 t)e3t = 3(c1 t2 + c2 t)e3t + te3t

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

**Paso 3: Solución general**


Combinando las soluciones homogénea y particular:
       
t 2t
1 1 0 −e + e
       
Y (t) = c1 −1 + c2 e2t 1 + c3 e3t 0 +  e2t 
       
       
1 2 3t
0 0 1 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

L2 di2 + R2 (i2 − i1 ) = 0



dt

Sustituyendo los valores dados:



 di1 + 11i1 − 3i2 = sin(100t)

dt

 di2 − 3i1 + 3i2 = 0



dt

**Paso 2: Solución del sistema homogéneo**


La matriz del sistema homogéneo es:
 
−11 3
A= 
3 −3
Valores propios:

−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 300 3300


A=− , B= , C= , D=
10009 10009 10009 10009
**Paso 4: Solución general**
Combinando las soluciones homogénea y particular:

       
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)

**Paso 5: Aplicación de condiciones iniciales**


Para t = 0:

0 = c1 + 3c2 −
 300
10009

0 = 3c1 − c2 +
 300
10009

Resolviendo:

100
c1 = 0, c2 =
10009
**Solución final:**

300 −12t 300 1100


i1 (t) = e − cos(100t) + sin(100t)
10009 10009 10009
100 −12t 300 3300
i2 (t) = − e + cos(100t) + sin(100t)
10009 10009 10009

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:

X ′ (t) = AX(t), X(0) = X0

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 ,

donde D = diag(λ1 , λ2 , . . . , λn ) contiene los valores propios de A, y las columnas de P


son los vectores propios correspondientes.
**Paso 2: Cambio de variable** Definimos una nueva variable Y (t) mediante:

X(t) = P Y (t).

Sustituyendo en la ecuación diferencial:

P Y ′ (t) = AP Y (t) = P DP −1 P Y (t) = P DY (t).

Multiplicando ambos lados por P −1 :

Y ′ (t) = DY (t).

**Paso 3: Solución del sistema desacoplado** El sistema Y ′ (t) = DY (t) es diagonal,


por lo que cada componente yi (t) satisface:

yi′ (t) = λi yi (t), yi (0) = (P −1 X0 )i .

La solución única para cada componente es:

yi (t) = eλi t (P −1 X0 )i .

En forma matricial:
Y (t) = eDt P −1 X0 ,

donde eDt = diag(eλ1 t , eλ2 t , . . . , eλn t ).


**Paso 4: Transformación de vuelta a X(t)** Aplicando el cambio de variable inverso:

X(t) = P Y (t) = P eDt P −1 X0 = eAt X0 ,

donde eAt = P eDt P −1 es la exponencial de la matriz A.

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

tiene como **solución única**:


X(t) = eAt 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:

e−2t  2 e2t e−2t


−2e2t cos(5t) + e4t + 1 = − cos(5t) +

g(t) = +
58 58 58 58

Simplificamos los coeficientes:

1 1 1
g(t) = − cos(5t) + e2t + e−2t
29 58 58

**Paso 2: Aplicar la transformada de Laplace a cada término**


Utilizamos las siguientes propiedades de la transformada de Laplace: 1. L{cos(at)} =
s 1
s2 +a2
2. L{ebt } = s−b

Aplicamos a cada término de g(t):


1
1. Para − 29 cos(5t):
 
1 1 s
L − cos(5t) = − · 2
29 29 s + 25
1 2t
2. Para 58
e :  
1 2t 1 1
L e = ·
58 58 s − 2

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π,

y extendida periódicamente con periodo T = 3π.


**Paso 1: Expresar h(t) en términos de funciones escalón unitario**
Usamos la función escalón unitario u(t−a) para describir los cambios en cada intervalo:
πh  π i h  π i
h(t) = u(t) − u t − +(π−t) u t − − u(t − π) +π [u(t − π) − u(t − 2π)]+(t−2π) [u(t − 2π
2 2 2
**Paso 2: Aplicar la transformada de Laplace a cada término**
Calculamos la transformada de Laplace para cada tramo, utilizando las propiedades
de linealidad y desplazamiento:
π
1. **Para 2
en [0, π2 ):**
nπ o π
L u(t) = .
2 2s

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:

e−πs π e−2πs e−3πs


 
π π −πs − π π 1 π
L{h(t)} = − e 2 +e 2 s
− 2 + 2 + e−πs − e−2πs + 2 − 2 .
2s 2s 2s s s s s s s

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

**Paso 4: Expresión final**


La transformada de Laplace de h(t) es:
π
π e−πs + e−2πs π −πs −2πs
 e− 2 s + e−3πs
L{h(t)} = + + e − e − .
2s s2 s s2

**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)

Descomponemos G(s) en fracciones parciales:

1 As + B Cs + D
= + 2
(s2 + 4)(s2 − 9) s2 + 4 s −9

Multiplicamos ambos lados por el denominador común:

1 = (As + B)(s2 − 9) + (Cs + D)(s2 + 4)

Expandiendo y agrupando términos:

1 = (A + C)s3 + (B + D)s2 + (−9A + 4C)s + (−9B + 4D)

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)

Donde u(t − 5) es la función escalón unitario.


**Paso 4: Expresión final de la transformada inversa**

 
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

L2 dI2 + R2 I2 + R1 (I2 − I1 ) = 0



dt

Sustituyendo los valores dados:



5 dI1 + 30(I1 − I2 ) = 2u4 (t)

dt

10 dI2 + 100I2 + 30(I2 − I1 ) = 0



dt

37
Simplificando:

 dI1 + 6I1 − 6I2 = 2 u4 (t)

dt 5

 dI2 − 3I1 + 13I2 = 0



dt

**Paso 2: Aplicar la transformada de Laplace**


Aplicamos la transformada de Laplace a ambas ecuaciones, considerando las condicio-
nes iniciales nulas:

sL{I1 } + 6L{I1 } − 6L{I2 } = 2 e−4s
·

5 s

sL{I2 } − 3L{I1 } + 13L{I2 } = 0


Reorganizando términos:

(s + 6)L{I1 } − 6L{I2 } =
 2e−4s
5s

−3L{I1 } + (s + 13)L{I2 } = 0

**Paso 3: Resolver el sistema algebraico**


Expresamos el sistema en forma matricial:
    
2e−4s
s+6 −6 L{I1 }
   =  5s 
−3 s + 13 L{I2 } 0
Calculamos el determinante de la matriz de coeficientes:

∆ = (s + 6)(s + 13) − 18 = s2 + 19s + 60

Las soluciones son:

2e−4s 2e−4s
(s + 13) · 5s
3· 5s
L{I1 } = , L{I2 } =
∆ ∆
Simplificando:

2(s + 13)e−4s 6e−4s


L{I1 } = , L{I2 } =
5s(s2 + 19s + 60) 5s(s2 + 19s + 60)
**Paso 4: Descomposición en fracciones parciales**
Factorizamos el denominador:

s2 + 19s + 60 = (s + 4)(s + 15)

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

**Paso 2: Derivada de g(t)**


Observamos que:
g ′ (t) = f (t)

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):

L{g ′ (t)} = L{f (t)}

Utilizando la propiedad de la transformada de Laplace para derivadas:

sL{g(t)} − g(0) = F (s)

Dado que g(0) = 0, simplificamos:

sL{g(t)} = F (s)

**Paso 4: Despejar L{g(t)}**


Dividiendo ambos lados por 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

También podría gustarte