Repaso Laplace
Repaso Laplace
Repaso Laplace
Solucionario de
problemas de Ecuaciones
Diferenciales
2do Parcial (3ra versión)
• Resolución de ecuaciones diferenciales alrededor de puntos singulares
• Transformada de Laplace
• Resolución de ecuaciones diferenciales mediante la transformada de
Laplace
• Resolución de sistemas de Ecuaciones diferenciales
• Aplicaciones de las ecuaciones diferenciales de segundo orden
• Series de Fourier
• Ecuaciones en Derivadas Parciales
Roberto Cabrera V.
dcabrera@fiec.espol.edu.ec
06/02/2009
Transformada de Laplace
Halle:
• {
L 4e 5t + 6t 3 − 3 sen(4t ) + 2 cos(2t ) }
Por la propiedad de linealidad tenemos que:
{
L 4e 5t + 6t 3 − 3 sen(4t ) + 2 cos(2t ) }
L{4e 5t
+ 6t 3 − 3 sen(4t ) + 2 cos(2t )} = L{4e }+ L{6t }+ L{− 3 sen(4t )} + L{2 cos(2t )}
5t 3
1 3! 4 s
=4 +6 4 −3 2 +2 2
s−5 s s + 16 s +4
4 36 12 2s
= + − +
s − 5 s 4 s 2 + 16 s 2 + 4
Halle
• {
L (t + 2 ) e t + e −4t cosh (2t )
2
}
Por la propiedad de linealidad tenemos que:
{
L (t + 2) e t + e −4t cosh (2t )
2
}
L{(t + 2) e cosh (2t )} = L{(t + 2) e }+ L{e }
cosh (2t )
2 −4t 2 − 4t
t
+e t
{( ) } { cosh(2t )}
= L t + 4t + 4 e + L e
2 t − 4t
-8-
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Demuestre:
• Demuestre el primer teorema de la traslación
{
Si L{ f (t )} = F (s ), entonces L e at f (t ) = F (s − a ) }
∞
Tenemos : L{ f (t )} = ∫ e − st f (t )dt = F (s )
0
∞
{ }
Entonces : L e at f (t ) = ∫ e − st e at f (t )dt
0
∞
= ∫ e −(s −a )t f (t )dt → si s = s − a
0
()
∞
= ∫ e − st f (t )dt = F s = F (s − a )
0
Halle:
• {
L senh 3 (2t ) cos(t ) }
Por la propiedad de linealidad tenemos que:
{
L senh 3 (2t ) cos(t ) }
e 2t − e −2t 3
{
L senh (2t ) cos(t ) = L
3
} 2
cos(t )
1
{(
= L e 6t − 3e 2t + 3e −2t − e −6t cos(t )
8
) }
1
[{ } { } { } {
= L e 6t cos(t ) + L − 3e 2t cos(t ) + L 3e −2t cos(t ) + L − e −6t cos(t )
8
}]
1
[{ } { } {
= L e 6t cos(t ) − 3L e 2t cos(t ) + 3L e − 2t cos(t ) − L e −6t cos(t )
8
} { }]
Aplicando el primer teorema de la traslación:
1
8
[{ } { } {
L e 6t cos(t ) − 3L e 2t cos(t ) + 3L e −2t cos(t ) − L e −6t cos(t ) } { }]
1 s−6 s−2 s+2 s+6
= −3 +3 −
8 (s − 6 ) + 1
2
(s − 2) + 1 (s + 2) + 1 (s + 6)2 + 1
2 2
= 2
(
48 s 4 − 46 s 2 + 185 )
( )( )(
s − 12 s + 37 s 2 − 4 s + 5 s 2 + 4 s + 5 s 2 + 12 s + 37 )( )
-9-
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Si L{ f (t )} = F (s ), entonces L{ f' (t )} = sF (s ) − f (0 )
∞ P
Tenemos : L{ f' (t )} = ∫ e − st f ' (t )dt = lim ∫ e − st f ' (t )dt
P →∞
0 0
− st
Integrando por partes : u = e → du = -s e - st dt
dv = f ' (t )dt → v = f (t )
P
P
f ' (t )dt = lim f (t )e − st + s ∫ e − st f (t )dt
P
lim ∫ e − st
P →∞ P →∞ 0
0 0
− sP P
= lim e f (P ) − f (0 ) + s ∫ e − st f (t )dt
P →∞
0
∞
= s ∫ e − st f (t )dt − f (0 ) + lim e − sP f (P )
P →∞
0
= sF(s) - f(0)
d
Si L{ f (t )} = F (s ), entonces L{tf (t )} = − F (s )
ds
∞
Tenemos : L{ f (t )} = ∫ e − st f (t )dt = F (s )
0
d
→ L{tf (t )} = − F (s )
ds
- 10 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
• {
L t 2 cos(at ) }
Por la propiedad de la derivada de la transformada tenemos que:
{
L t 2 cos(at ) }
L{t cos(at )} = (− 1)
2 d2 2
F (s)
ds 2
d2 s
= 2 2
ds s + a 2
d a 2 − s 2
=
(
ds s 2 + a 2 2 )
=
(
− 2s s 2 + a 2 )
2
( ) ( )
− 2 s 2 + a 2 2s a 2 − s 2
(
s2 + a2 )
2
=
(
2s s 2 − 3a 2 )
(
s2 + a2
3
)
- 11 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Halle:
•
L
cos t
( )
t
Usando la propiedad de la transformada de la derivada
L
( )
cos t
t
( )
Si f (t ) = sen t , entonces f' (t) =
cos t( ), además f(0) = 0
2 t
L{ f ' (t )} = sF (s ) − f (0)
( )
cos t
L = sL sen t { ( )}
2 t
( )
cos t
L = 2 s L sen t { ( )}
t
Encuentro la transformada de sen t ( )
2 n +1
( t) + ( t) − ( t)
3 5 7
3 5 7
sen( t ) = t −
1
t2 t2 t2
+ .... = t − + − + .... 2
3! 5! 7! 3! 5! 7!
{ ( )}
Γ 3
L sen t = 3 2 −
Γ( ) ( ) ( ) ( )
5
2 +
Γ 7
2 −
Γ 9
2 + ....
5 7 9
s 2
3! s 2
5! s 2
7! s 2
π 1
= 3 1 − 2 +
1 2
2 s
2
−
( ) ( )
1 2
2 s
3
+ ...
2s
2 2 s 2! 3!
π
1
−
= 3
e 4s
2s 2
L
( )
cos t π − 4s
= 2s 3 e
1
π − 4s
= 1 e
1
t 2s 2 s 2
- 12 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
t F (s )
Si L{ f (t )} = F (s ), entonces L ∫ f (u )du =
0 s
t
Si g (t ) = ∫ f (u )du , entonces g' (t) = f(t) y g(0) = 0
0
f (t )
∞
Si L{ f (t )} = F (s ), entonces L = ∫ F (u )du
t s
f (t )
Si g (t ) = , entonces f (t) = t g(t)
t
Entonces sabemos que :
L{ f (t)} = L{t g(t)}
d
L{ f (t)} = −L{g (t)}
ds
Integrando ambos lados tenemos que :
s ∞
L{g (t)} = − ∫ f (u )du = ∫ f (u )du
∞ s
f (t )
∞
L = ∫ f (u )du
t s
- 13 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Halle:
4t t 1 − 4θ
• L e t ∫ e sen(3θ )dθ
0
θ
t
1
L e 4t t ∫ e − 4θ sen(3θ )dθ
0
θ
Por el primer teorema de la traslación tenemos que :
{ }
L e 4t g (t ) = G (s − 4 )
t 1
Debo encontrar G(s) que es L t ∫ e − 4θ sen(3θ )dθ ,
0θ
por el teorema de la derivada de la transformada sabemos que :
d
L{t h(t )} = −H ( s)
ds
t 1 − 4θ M ( s) 1
Encuentro H(s) = L ∫ e sen(3θ )dθ = , si M ( s ) = L e − 4θ sen(3θ )
0 θ s θ
∞
x( θ )
De donde hallamos M(s) = L = ∫ X (u )du
θ s
{ }
X (u ) = L e − 4θ sen(3θ ) que por el primer teorema de traslación es :
X (u ) = L{e sen(3θ )} =
− 4θ 3 3
=
(u + 4) 2
+9 u + 8u + 25
2
∞ ∞ ∞
x( θ ) 3 u + 4 π s + 4
M(s) = L = ∫ X (u )du = ∫ 2 du = arctan = − arctan
θ s s u + 8u + 25 3 s 2 3
M (s) 1 π s + 4
H(s) = = − arctan
s s 2 3
s + 4
arctan
d 1 π s + 4 π 3 3
G ( s ) = L{t h(t )} = − − arctan = 2 + −
ds s 2 3 2 s s s 2 + 8s + 25 s2 ( )
s
arctan
t
1 π 3 3
L e 4t t ∫ e − 4θ sen(3θ )dθ = G (s − 4 ) = + −
0
θ 2( s − 4 )
2
( 2
)
(s − 4) (s − 4) + 8(s − 4 ) + 25 (s − 4)2
- 14 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Si t = u + a → u = t - a y dt = du
Cuando t = a → u = 0 y t = ∞ → u = ∞
∞
L{µ (t - a ) f (t − a )} = ∫ e − s (u + a ) f (u )du
0
∞
L{µ (t - a ) f (t − a )} = e − as
∫e
− su
f (u )du = e −as F (s )
0
t
[µ 0 (t ) − µ1 (t ) + µ 2 (t ) − µ 3 (t ) + µ 4 (t ) − µ 5 (t )....]
−
f (t ) = e 2
1
Por el primer teorema de la traslación tenemos que L{ f(t)} = G s +
2
G ( s ) = L{µ 0 (t ) − µ1 (t ) + µ 2 (t ) − µ 3 (t ) + µ 4 (t ) − µ 5 (t )...}
1 e − s e − 2 s e −3 s e − 4 s
G (s) =
s
−
s
+
s
−
s
+
s
1
s
(
− ... = 1 − e − s + e − 2 s − e −3s + e − 4 s − ... )
+∞ n
1 1 1 1 es
G (s) = ∑ − s = =
s n=0 e s
1+ s
1 s es +1 ( )
e
s+ 1
1 e 2
L{ f(t)} = G s + =
2 1 s+ 1
s + e 2 + 1
2
- 15 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
sen(3t )
• L sen(t )µ π (t ) + δ (t )
4
t
sen(3t )
L sen(t )µ π (t ) + δ (t )
4
t
sen(3t ) sen(3t )
L sen(t )µ π (t ) + δ (t ) = L sen(t )µ π (t ) + L δ (t )
4
t 4 t
Para la primera transformada utilizo el segundo teorema de la traslación :
π
π − s
L f t − µ π (t ) = e 4 F ( s )
4 4
Pero debo desplazar la función que multiplica al escalón :
π π π π π π 2 π π
sen t − + = cos t − sen + sen t − cos = cos t − + sen t −
4 4 4 4 4 4 2 4 4
2 π π 2 π π
L cos t − + sen t − µ π (t ) = L cos t − µ π (t ) + L sen t − µ π (t )
2 4 4 4 2 4 4 4 4
π π
2 −4s s 1 2 −4s s +1
= e 2 + 2 = e 2
2 s + 1 s + 1 2 s + 1
Para la segunda transformada utilizo la función impulso :
sen(3t ) sen(3t )
L δ (t ) = e −0 s lim = 1(3) = 3
t t →0 t
sen(3t )
π
2 −4s s +1
L sen(t )µ π (t ) + δ (t ) = e 2 +3
4
t 2 s + 1
- 16 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
∫e
− st
(
sen(t )dt = − cos(t )e - st − s e - st sen(t ) + s ∫ e - st sen(t )dt )
(s 2
)
+ 1 ∫ e − st sen(t )dt = e - st (− s sen(t ) − cos(t ) )
e − st (− s sen(t ) − cos(t ) )
∫ e sen(t )dt =
− st
s2 +1
Re emplazando :
π
1 e − st (− s sen(t ) − cos(t ) )
L{g (t )} =
1 − e −2πs s2 +1 0
1 1 + e −πs 1
L{g (t )} = − 2πs =
1− e s +1 1− e
2 −πs
s2 +1 ( )( )
- 17 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
t
Si L {F(s)} = f (t ) y L {G(s)} = g (t ), entonces L{F ( s )G ( s )} = ∫ f (u ) g (t − u )du = f (t ) * g (t )
-1 -1
t
L ∫ f (u ) g (t − u )du = F ( s )G ( s )
0
donde F(s) = L{ f(t)}, G(s) = L{g(t)} por lo que :
t ∞ t ∞ t
L ∫ f (u ) g (t − u )du = ∫ e − st ∫ f (u ) g (t − u )du dt = ∫ ∫e
− st
f (u ) g (t − u )du dt = lim S M
M →∞
0 t =0 u =0 t =0 u =0
M t
∫ ∫e
− st
donde S M = f (u ) g (t − u )du dt
t =0 u =0
t-u=v
0 M 0 M
Luego de hacer el cambio t-u=v la región cambia, por lo que el integral se transforma en:
∂ (u , t )
S M = ∫∫ e − st f (u ) g (t − u )du dt = ∫∫ e − s (u + v ) f (u ) g (v) du dv
Rtu Ruv
∂ (u , v )
Donde el Jacobiano de la transformación es :
∂u ∂u
∂ (u , t ) ∂u ∂v = 1 0 = 1
J= =
∂ (u , v ) ∂t ∂t 1 1
∂u ∂v
M M −v
− s (u + v )
De donde S M = ∫ ∫e
v =0 u =0
f (u ) g (v)du dv
e − s (u + v ) f (u ) g (v) u + v ≤ M
Definamos otra función K(u,v) =
0 u+v > M
M M ∞ ∞
SM = ∫
v =0 u =0
∫ K (u, v)du dv , entonces lim S M = M →∞ ∫ ∫ K (u, v)du dv
v =0 u =0
=
∞
∫ ∫e
v =0 u =0
∞
− s (u + v )
f (u ) g (v) du dv = {∫ e
0
∞
− su
f (u ) du }{∫ e
0
∞
− sv
}
g (v) dv = F ( s )G ( s )
- 18 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Halle:
s
• L−1 2
( )
s 2 + a 2
s
L−1
( ) 2
s 2 + a 2
Usando el integral de convolución tenemos que :
s 1 1
L−1 2 = cos(at ) * sen(at )
( 2 2
)(
s +a s +a
2
)
a
s 1 t sen(a(t − u ))
L−1 2 = ∫ cos(au ) du
( 2 2
)(
s +a s +a
2
)
0 a
t
1
= ∫ cos(au )[sen(at ) cos(au ) − cos(at )sen(au )]du
a0
t t
1 1
= sen(at )∫ cos 2 (au )du − cos(at )∫ sen(au ) cos(au )du
a 0
a 0
1 + cos(2au ) sen(2au )
t t
1 1
= sen(at )∫ du − cos(at )∫ du
a 0
2 a 0
2
1 t sen(2at ) 1 1 − cos(2at )
= sen(at ) + − cos(at )
a 2 4a a 4a
1 t sen(at ) cos(at ) 1 sen 2 (at )
= sen(at ) + − cos (at )
a 2 2a a 2 a
t sen(at )
=
2a
- 19 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
• y ' ' '+4 y ' '+5 y '+2 y = 10 cos(t ), y (0) = y ' (0) = 0 y ' ' (0) = 3
(s + 4s + 5s + 2)Y (s) − 3 = 10 s s+ 1
3 2
2
2
+ 10s + 3
(s + 1)2 (s + 2 )Y ( s) = 3s
s2 +1
3s 2 + 10 s + 3 A B C Ds + E
Y ( s) =
(s + 1)(s + 1) (s + 2)
2 2
= + +
s + 2 s + 1 (s + 1)2
+ 2
s +1
→ 3s 2
+ 10 s + 3 = A(s + 1) (s + 1)+ B (s + 1)(s + 1)(s + 2) + C (s + 1)(s + 2 ) + (Ds + E )(s + 1) (s + 2)
2 2 2 2 2
- 20 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
(s 2 + 4 Y (s) =) 2 s 3 + 8πs − 4
s 2
4
+ e − 2 πs 2
s
2 s + 8πs − 4
3
4
Y (s) = + e − 2 πs 2 2
2 2
s s +4 ( ) (
s s +4 )
Encuentro fracciones parciales :
2 s 3 + 8πs − 4 A B Cs + D
• = + 2 + 2
(
s2 s2 + 4 s s )
s +4
( ) (
2 s + 8πs − 4 = As s 2 + 4 + B s 2 + 4 + (Cs + D )s 2
3
)
2 s 3 + 8πs − 4 = ( A + C )s 3 + (B + D )s 2 + (4 A )s + (4 B )
A+C = 2
B+D=0
4 A = 8π
4 B = −4
Re solviendo el sistema tenemos que : A = 2π , B = -1, C = 2 - 2π , D = 1
4 A B Cs + D
• 2 2 = + 2 + 2
s s +4 ( s s )
s +4
( ) ( )
4 = As s + 4 + B s 2 + 4 + (Cs + D )s 2
2
4 = ( A + C )s + (B + D )s + (4 A )s + (4 B )
3 2
A+C =0
B+D=0
4A = 0
4B = 4
Re solviendo el sistema tenemos que : A = 0 , B = 1, C = 0, D = −1
2π 1 (2 - 2π )s + 1 1 1
Y (s) = − 2 + 2
+ e − 2πs 2 − 2
s s s +4 s s +4
sen ( 2t ) sen (2(t − 2π ))
y (t ) = 2π − t + (2 - 2π ) cos( 2t ) + + µ 2π (t ) (t − 2π ) −
2 2
- 21 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Despejando Y(S):
i)
ii)
iii) Entonces
iv)
- 22 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
- 23 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
d
L{(1 − 2t ) y '} = L{y '} − 2 L{ty '} = (sY ( s ) − y (0 ) ) + 2 [sY ( s ) − y (0 ) ]
ds
L{(1 − 2t ) y '} = (sY ( s ) − 1) + 2 (Y ( s ) + sY ' ( s ) ) = 2 sY ' ( s ) + (s + 2 )Y ( s ) − 1
L{y} = Y ( s )
Re emplazando :
(− s Y ' ( s ) − 2 sY ( s ) + 1) + (2 sY ' ( s ) + (s + 2 )Y ( s ) − 1) − 2Y ( s ) = 0
2
(− s + 2 s )Y ' ( s ) + (− 2 s + s + 2 − 2 )Y ( s ) = 0
2
− s (s − 2 )Y ' ( s ) − sY ( s ) = 0
− s (s − 2 )Y ' ( s ) = sY ( s )
Y ' (s) s
=
Y ( s ) − s (s − 2 )
Y ' (s) ds
∫ Y (s)
=− ∫ (s − 2 )
ln (Y ( s ) ) = − ln (s − 2 ) + ln( K )
K
Y (s) = → y (t ) = Ke 2 t
s−2
y ( 0 ) = Ke 2 ( 0 ) = 1 → K = 1
y (t ) = e 2 t
- 24 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
d
L{(t + 2 ) y '} = L{ty '} + 2 L{y '} = − [sY ( s ) − y (0 ) ] + 2 (sY ( s ) − y (0 ) )
ds
L{(t + 2 ) y '} = 2 (sY ( s ) − k1 ) − (Y ( s ) + sY ' ( s ) ) = − sY ' ( s ) + (2 s − 1)Y ( s ) − 2 k1
L{y} = Y ( s )
1 1 1− s
L{t } − L{1} =− = 2
s2 s s
Re emplazando :
(− s 2
)
+ s Y ' ( s ) + (− 2 s − 2 s + 1 + 3 )Y ( s ) + (k1 + 2 k1 ) =
1− s
s2
1 − s − 3k1 s 2
− s (s − 1)Y ' ( s ) − 4 (s − 1)Y ( s ) =
s2
2 − 1 + s + 3k1 s 2
Y ' ( s) + Y (s) =
s s 3 (s − 1)
2
u (s ) = e
∫ s ds = e 2 ln (s ) = (s )2
− 1 + s + 3k1 s 2
u (s )Y ( s ) = u (s ) ∫ ds
s (s − 1)
3
2
(s )2 Y ( s ) = − 1 + s + 3k1s ds = 1 + 3k1 1 ds
∫ ∫
s (s − 1) s s −1
s 2Y ( s ) = ln (s ) + 3k1 ln (s − 1) + k 2 = ln s (s − 1) ( 3 k1
)+ k 2
Y (s) =
ln s (s − 1)( k
3 k1
)
+ 22 → y (t ) = L−1 {Y ( s )}
2
s s
ln s (s − 1)
y (t ) = L−1
(
3 k1
k
+ 22 = L−1
)
ln s (s − 1)3 k1 ( ) + L −1 k 2
2
2
s s s2 s
L−1
(
ln s (s − 1)3 k1 ) = L −1
(
ln s (s − 1)
3 k1
) s1 = L −1
{F ( s )G ( s )} = f (t ) * g (t )
s2
2
−1 (3k1 + 1)s − 1
{(
f (t ) = L−1 ln s (s − 1)
3 k1
)}→ tf (t ) = − L −1 d
( 3 k1
ln s (s − 1) = − L )
ds s ( s − 1)
(3k1 ) 1
tf (t ) = − L−1 + = − 3k1e t + 1 → f (t ) = (
− 3k1e t + 1
) ( )
( s − 1) s t
1
g (t ) = L−1 2 = t
s
(
− 3k1e t + 1 ) t
(
− 3k1e u + 1 )
(t − u )du + k 2t
y (t ) =
t
* t + k 2 t → y (t ) = ∫
0
u
- 25 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Método de eliminación
Se obtiene:
x'−3x − 2 y = e −t − e t
Despejando y :
x' 3x e t − e −t
y= − +
2 2 2
Reemplazando y en (1):
'
x' 3x e t e − t x' 3x e t e − t
2 x'+ − + − − x − − + − = e − t
2 2 2 2 2 2 2 2
t −t t −t
x' ' 3x' e e x' 3x e e
⇒ 2 x'+ − + + −x− + − + = e −t
2 2 2 2 2 2 2 2
x' ' x
⇒ + =0 ⇒ x' '+ x = 0
2 2
⇒ si x = e rt
⇒ e rt [r 2 + 1] = 0
⇒ r2 + 1 = 0
⇒ r1 , 2 = ±i
⇒ x = C 1 cos t + C 2 sent
x' 3x e t e − t
⇒ y= − + −
2 2 2 2
C1 C 3 e t e −t
⇒ y=− sent + 2 cos t − (C 1 cos t + C 2 sent ) + −
2 2 2 2 2
t −t
C 3C 2 C 2 3C 1 e e
y = − 1 − sent + − cos t + −
144
2 24243 1242
4 43 24 2 2
K1 K2
e e−t
t
e t − e −t
⇒ y = k 1 sent + k 2 cos t + − ; Pero = senh t
2 2 2
Solución:
x = C 1 cos t + C 2 sent
y = k 1 sent + k 2 cos t + senht
- 26 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
2) Utilice el método de eliminación para encontrar la solución general del sistema lineal
dado, donde x´, y´, z´ denotan diferenciación con respecto a t.
dx
= 2x − 3y 1
dt
dy
= y − 2x
dt 2
⇒ 2 x´ − x´´ = 2 x − x´ − 6 x
⇒ x´´ −3 x´ −4 x = 0
Obtenemos una ecuación diferencial de coeficientes constantes:
Resolviendo la ecuación 3 con x=ert;
[
⇒ e rt r 2 − 3 r − 4 = 0 ]
Ecuación Característica
⇒ r 2 − 3r − 4 = 0
( r − 4 )( r + 1 ) = 0
⇒ r1 = 4 , r2 = −1
⇒ x1 = e 4t , x 2 = e −t
⇒ x = c1e 4t + c 2 e −t ;
Ahora encontremos y:
2
y= ( x ) − 1 ( x ')
3 3
x = C1e + C 2 e −t
4t
x' = 4C1e 4t − C 2 e −t
- 27 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
⇒ Reemplazando x, y x’ en y:
2
3
[ 1
] [
⇒ y = C1e 4t + C 2 e −t − 4C1e 4t − C 2 e −t
3
]
2
y = − C1e 4 t + C 2 e −t
3
x = C1e 4t + C 2 e − t
2
y = − C1e 4 t + C 2 e −t
3
- 28 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
(D − 2 ) 2 x 1 + D(D + 2 )x 2 = t
D(D − 2 )x 1 + (D + 2 ) x 2 = e
2 t
t D(D + 2 )
et ( D + 2 )2 ( D + 2 )2 t − D( D + 2 )e t (D + 2 )[(D + 2 )t − De t ]
x1 (t ) = = =
(D − 2 ) 2
D( D + 2 ) (D − 2 ) ( D + 2 ) − D(D + 2 )D( D − 2 ) (D + 2 )(D − 2 )[(D + 2 )(D − 2 ) − D 2 ]
2 2
D(D − 2 ) ( D + 2 )2
(D + 2 )(1 + 2 t − e t ) 2 − e t + 2 + 4t − 2e t 4 + 4 t − 3e t
x1 (t ) = = =
(D − 4)(D − 4 − D
2 2 2
) − 4(D 2 − 4 ) − 4( D 2 − 4 )
4 + 4 t − 3e t
x1 (t ) = ;
− 4(D 2 − 4 )
− 4(D 2 − 4 )x 1 ( t ) = 4 + 4 t − 3e t ;
3 t
(D 2 − 4 )x 1 ( t ) = −1 − t + e ;
4
3
x 1 ' ' ( t ) − 4 x 1 ( t ) = −1 − t + e t ;
4
Encontrando la solución hom ogénea:
x1 ''( t) − 4x1 (t ) = 0 ;
x 1 ( t ) = e rt ;
e rt [r 2 − 4 ] = 0 ;
r 2 − 4 = 0; r1 , 2 = ±2 ;
x 1 h ( t ) = C 1e + C 2 e −t ;
2t
- 29 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
r1 , 2 = ±2 ;
x 2 h (t ) = C 1 e 2 t + C 2 e −2 t ;
Encontrand o la solución particular :
x 2 p = a + be t ;
x' 2 p = be t ;
x' ' 2 p = be t ;
et 1
Reemplazando x 2p en x 2 ' ' (t) − 4x 2 (t) = − − :
4 2
et 1
be t − 4(a + be t ) = − − ;
4 2
et 1
− 3be t − 4a = − − ;
4 2
1 1
b= ; a= ;
12 8
1 et
x 2p = + ;
8 12
1 et
x 2 (t ) = C 1 e 2 t + C 2 e − 2 t + + ;
8 12
La solución es:
−t 1 1 1 t
x 1 (t) = C 1 e + C 2 e + 4 + 4 t − 4 e ;
2t
t
x ( t ) = C e 2 t + C e − 2 t + 1 + e ;
2 1 2
8 12
- 30 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
( D + 2 )( x 1 ) + (D − 1)( x 2 ) = −sent ( 1)
(D − 3)( x 1 ) + (D + 2 )( x 2 ) = 4 cos t ( 2 )
- 31 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
- 32 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Método de Laplace
1) Utilice el método de las transformadas de Laplace para resolver el problema de valor
inicial dado. Aquí x’, y’, etc. denotan diferenciación con respecto a t.
1
4
5 x( s ) − x(0) − 3x( s ) + 2 y ( s ) = s 2 + 1 ; 1
(−4)( s − 3) x ( s ) − 8 y ( s ) = − s 2 + 1 ;
≈
4 x( s ) − sy ( s ) − y (0) − y ( s ) = s ; 4( s − 3) x ( s ) − ( s − 3)( s + 1) y ( s ) = s ( s − 3) ;
s2 +1 2
s2 +1
1
( s − 3) x( s ) + 2 y ( s ) = 2
(−4) s +1
( s − 3) s
4 ( ) − ( s + 1) y ( s ) = 2
x s s +1
s 2 − 3s − 4 ( A + C )( s 3 ) + ( B + D − 2C ) s 2 + ( A − 2 D + 5C ) s + ( B + 5D)
⇒ 2 =
( s − 2s + 5)( s 2 + 1) ( s 2 − 2s + 5)( s 2 + 1)
A + C = 0 A = 11 / 10 ;
B + D − 2C = 1 B = −1 / 2 ;
⇒ Resolviendo el sistema:
A − 2D + 5C = −3 C = −11 / 10 ;
B + 5D = −4 D = −7 / 10 ;
11 1 11 7
10 s − 2 − s−
! y(s) = − 2 + 102 10
s − 2 s + 5 s + 1
- 33 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
1 11 11 7
− s s+
y(s ) = 2 10 + 10 2 10
[
(s − 1) + 4
2
] s +1
1 1 11 s 11 s 7 1
y(s ) = ⋅ − ⋅ + ⋅ 2 + ⋅ 2
[ 2
] [ 2
]
2 (s − 1) + 4 10 (s − 1) + 4 10 s + 1 10 s + 1
1 2 11 (s − 1 + 1) 11 s 7 1
y(s ) = ⋅ − ⋅ + ⋅ 2 + ⋅ 2
[ 2
] [2
]
4 (s − 1) + 4 10 (s − 1) + 4 10 s + 1 10 s + 1
1 2 11 (s − 1) 11 2 11 s 7 1
y(s ) = ⋅ − ⋅ − ⋅ + ⋅ 2 + ⋅ 2
[ ] [ ] [
4 (s − 1) + 4 10 (s − 1) + 4 20 (s − 1) + 4 10 s + 1 10 s + 1
2 2 2
]
! Aplicando transformada inversa de Laplace a y(s):
L−1 [ y ( s )] = y (t ) ;
1 2 11 − 1 (s − 1) 11 − 1 2 11 − 1 s 7 − 1 1
y( t ) = L− 1 − L − L + L 2 + L 2
[(s − 1) + 4] 10 [(s − 1) + 4] 20 [(s − 1) + 4] 10 s + 1 10 s + 1
2 2 2
4
1 −t 11 11 11 7
y( t ) = e sen (2 t ) − e − t cos(2 t ) − e − t sen (2 t ) + cos(t ) + sen (t )
4 10 20 10 10
3 −t 11 − t 11 7
y( t ) = − e sen (2 t ) − e cos(2 t ) + cos(t ) + sen (t )
10 10 10 10
y ′( t ) 1 − t 5 11 7
= e cos( 2 t ) + e − t sen (2 t ) − sen( t ) + cos(t )
4 5 8 40 40
y( t ) 3 11 − t 11 7
= − e − t sen( 2 t ) − e cos(2 t ) + cos( t ) + sen (t )
4 40 40 40 40
La solución:
3 −t 11 −t 1 7
x( t ) = − 40 e cos( 2 t ) + 20 e sen (2 t ) − 10 sen( t ) + 10 cos(t )
y( t ) = − 3 e −t sen (2 t ) − 11 e − t cos(2 t ) + 11 cos(t ) + 7 sen (t )
10 10 10 10
- 34 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
[ 2
]
s X(s) − sx(0) − x'(0) + [sY( S ) − y(0 )] − X( S ) =
15
s+1
;
[ ]
s 2 Y(s) − sy(0 ) − y'(0 ) − 4[sX( S ) − x(0 )] − Y( S ) = 15 ;
s2
[ 2
] 15
s X(s) + [sY( S )] − X( S ) = s + 1 ;
[ ]
s 2 Y(s) − 4[sX( S )] − Y( S ) = 15 ;
s2
2
(s − 1)X(s) + sY(S ) = 15 ;
s+1
− 4sX( S ) + (s − 1)Y(s) = 15 ;
2
s2
Aplicando la regla de Kramer :
15
s
s+1
15
(s2 − 1) 15(s − 1) − 152s 15(s − 1)(s + 1) − 15 15(s − 1) − 15 15(s − s ) − 15
2 2
2
X( S ) = 2s = s2 + 1 2 s = s+1 s = 1 s = s
(s − 1) s (s − 1) + 4s 2
s + 2s + 1
4 2
(s + 1)
2 2
(s + 1)2
2
− 4s (s 2 − 1)
15(s 2 − s ) − 15 15(s 2 − s − 1) A Bs + C Ds + E
X( S ) = = = 15 + 2 + 2
s(s + 1)(s + 1) s(s + 1)(s + 1)
2 2 2 2
s (s + 1) (s + 1)
2
Expresando X(s) como la suma de fracciones parciales se obtiene que los valores de los coeficientes son :
A = -1, B = 2, C = -1, D = 1, E =0
Por lo tanto X(s) lo expresamos como :
− 1 2s − 1 s
X( S ) = 15 + 2 + 2
s (s + 1) (s + 1)
2
- 35 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
s -1 s
L-1 =£ 2 = (sent * cos t );
(s + 1) (s + 1)(s + 1)
2 2 2
t
sen( t ) + sen( 2 u − t ) usen( u ) cos( 2 u − t )
(sent * cos t ) = ∫0sen( u ) cos(t − u )du = ∫0
t t
du = −
2 2 4 0
t
tsen( t ) cos( t ) cos( −t ) tsent
(sent * cos t ) = − −0+ = ;
2 4 4 0 2
s tsent
L-1 = (sent * cos t ) = ;
(s + 1)
2
2
2
1 1 t t cos( 2 u − t ) − cos( t )
L-1 2 2
= £ -1 2 = (sent * sent ) = ∫ sen( u )sen( t − u )du = ∫ du ;
(s + 1)(s + 1)
(s + 1)
2 0 0
2
1 sen( 2 u − t ) u cos t t sen( t ) t cos t sen( −t ) sent − t cos t
L-1 2 2
= − = − − =
(s + 1) 4 2 0 4 2 4 2
1 sent − t cos t
L-1 2 2
=
(s + 1) 2
2s − 1 tsent sent − t cos t sent t cos t
L-1 2 2
= 2 − = tsent − +
(s + 1) 2 2
2 2
−1 2s − 1 s
x(t) = 15£ -1 + 15£ -1 2 2
+ 15£ -1 2
(s + 1)
s (s + 1)
sent t cos t
x(t) = −15 + 15 tsent − + + 15[cos t ]
2 2
x( t ) = −15 + 15tsent − 7.5sent + 7.5t cos t + 15 cos t
Ahora encontremos y(t) usando una ecuación del sistema:
X' '+ Y'+3X = 15e − t
y' = 15e − t − x' '−3x ;
x = −15 + 15tsent − 7.5sent + 7.5t cos t + 15 cos t ;
x' = 15t cos t + 15sent − 7.5 cos t − 7.5tsent + 7.5 cos t − 15sent ;
x' = 15t cos t − 7.5tsent ;
x' ' = −15tsent + 15 cos t − 7.5t cos t − 7.5sent ;
y' = 15e − t + 15tsent − 15 cos t + 7.5t cos t + 7.5sent − 3(− 15 + 15tsent − 7.5sent + 7.5t cos t + 15 cos t );
y' = 15e − t − 30tsent − 15t cos t − 60 cos t + 30sent ;
y= (∫ 15e −t
− 30tsent − 15t cos t − 60 cos t + 30sent )
y = −15e − t + 30t cos t − 30sent − 15tsent − 15 cos t − 60sent − 30 cos t + C ;
y(0) = 0 , entonces C = 45;
La solucion es :
x( t ) = −15 + 15tsent − 7.5sent + 7.5t cos t + 15 cos t
−t
y( t ) = −15e + 30t cos t − 15tsent − 90sent − 45 cos t + 45;
- 36 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
−4−λ 1 1
det(A − λI ) = 1 5−λ −1 = ( −4 − λ )[( 5 − λ )( −3 − λ ) + 1] − 1[( −3 − λ ) − 1] = 0
0 1 −3−λ
(λ + 4)(λ − 5)(λ + 3) = 0
λ1 = −4 λ2 = 5 λ3 = −3
λ = −4
x 10z 10
0 1 1 x 0 y = −z
y + z = 0 y ⇒ − z ⇒ − 1 z
1 9 − 1 y = 0 …. z z 1
x + 9 y − z = 0 x = 10 z
0 1 1 z 0
10
υ1 = − 1
1
λ =5
− 9 1 1 x 0 x z 1
− y + 8 z = 0 y = 8z
1 0 − 1 y = 0 ….. y ⇒ 8z ⇒ 8 z
0 1 − 8 z 0 x − z = 0 x=z z z 1
1
υ2 = 8
1
λ = −3
− 1 1 1 x 0 x z 1
− x + y + z = 0
1 8 − 1 y = 0 y ⇒ 0 ⇒ 0 z
0 1 0 z 0 x − z = 0 z z 1
- 37 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
1
υ3 = 0
1
x(t ) = c1υ1e λ1t + c 2υ 2 e λ2t + c3υ 3 e λ3t
10 1 1
− 4t 5t
x(t )c1 − 1e + c 2 8 e + c3 0 e −3t
1 1 1
2) Resolver el sistema
− 3 0 2
X' = 1 − 1 0 X
− 2 − 1 0
− 3 − λ 0 2
1 −1− λ 0 =0
−2 −1 − λ
- 38 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
( )
6x = 4 + 2 2 i z ;
3 y = (− 1 + 2 i )z ;
z ∈ Re ales.
Entonces:
2 2i
x = + z ; .
3 3
1 2i
y = − + z ;
3 3
z ∈ Re ales .
Entonces podemos concluir que el vector propio complejo asociado a este valor de
λ 3 = −1 − 2i es:
2
+ 2 i z
3 3
x
1 2i
v = y = − + z ,
z 3 3
z
Podemos usar un vector propio que tenga la forma de v, si z = 3 :
( )
x 2 + 2 i 2 2
( )
v = y = − 1 + 2 i = − 1 + 2 i ,
z 3
3 0
12 3 123
a b
2 2 2 2
( ) ( ) ( ) ( )
x 2 = e −t cos − 2 t − 1 − e − t sen − 2 t 2 = e − t cos 2 t − 1 + e −t sen 2 t 2 ;
3 3
0 0
- 39 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
2 2 2 2
−t −t
−t
( ) ( )
−t
( )
x 3 = e sen − 2 t − 1 + e cos − 2 t 2 = −e sen 2 t − 1 + e cos 2 t 2 ; ( )
3 3
0 0
Por lo tanto la solución general es:
x = C 1 x 1 +C 2 x 2 + C 3 x 3 ;
2 2 2 2 2
−t
−2 t
x = C 1 e − 2 +C e cos 2( )
t − 1 + e −t
sen ( )
2 t 2 + C 3 − e −t
( )
sen 2 t − 1 + e −t
( )
cos 2 t 2 ;
1
0
2
3 3 0
- 40 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
a) Determine una función que defina la posición ‘’y’’ de la masa en cualquier instante
‘’t’’.
b) Halle la posición de la masa en los tiempos t= π / 4 segundos y t= π segundos.
d2 y dy
m 2 +C + Ky = f ( t )
dt dt
d2 y π
2
+ 9Ky = −3δ t − ;
dt 2
Para resolver esta ecuación diferencial aplicamos la transformada de Laplace a ambos lados de la ecuación:
d2 y ϖ
L 2 + 9 y = L − 3δ t − ;
dt 2
π
− s
s 2 Y(s) − sy(0 ) − y' (0 ) + 9 y(s) = −3e 2 ;
La posición inicial del sistema es y(0)=1 metro, y la velocidad inicial es y’(0)=0:
π
− s
s 2 y(s) − s + 9 y( s) = −3e 2 ;
π
(s + 9 )y(s) = s − 3e
− s
2 2
;
π
− s π
− s
s − 3e
2
s 3e 2
y(s) = 2 = 2 − 2 ;
s +9 s +9 s +9
− s
π
− π2 s
-1 s 3e 2 -1 s -1 3e ;
y( t ) = L − =L 2 −L
s2 + 9 s2 + 9 s + 9 s2 + 9
π π
y( t ) = cos 3t − sen 3 t − u t − ;
2 2
- 41 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
a)
π
cos 3t ; t<
2
y( t ) =
cos 3t − sen 3 t − π t≥
π
2 2
b)
2
y( π / 4 ) = cos( 3π / 4 ) = − m,
2
y( π) = cos( 3π) − sen 3( π − π / 2 ) = −1 − ( −1) = 0m
100t ; t ∈ [0 ,2 )
f( t ) =
400 − 100t ; t ∈ (2 , 4 ]
d2 y dy
m 2 +c + Ky = f( t );
dt dt
Asumiendo que la gravedad es 10m / s 2 :
w 20
m= = = 2 Kg.
g 10
d2y dy
2 2 +6 + 4 y = f( t );
dt dt
Antes de resolver la ecuación diferencial aplicando la transformada de Laplace, se recomienda expresar la
función f(t) en términos de de funciones escalones multiplicadas por las funciones que se encuentran en cada
uno de los intervalos mostrados en la regla de correspondencia:
f(t ) = 100tu( t) − 100tu(t − 2) + (400 − 100t )u( t − 2) − (400 − 100t )u(t − 4);
f(t ) = 100tu( t) − 100(t )u(t − 2) + (400)u( t − 2) − 100tu( t − 2) − (400)u(t − 4) + 100tu(t − 4);
f(t ) = 100tu( t) − 200(t )u(t − 2) + (400)u(t − 2) − (400)u(t − 4) + 100tu( t − 4);
f(t ) = 100tu( t) − 200(t − 2 + 2 )u(t − 2) + (400)u( t − 2) − (400)u(t − 4) + 100(t − 4 + 4)u(t − 4);
f(t ) = 100tu( t) − 200(t − 2 )u(t − 2) − 400u(t − 2) + (400)u(t − 2) − (400)u( t − 4) + 100(t − 4)u( t − 4) + 400u( t − 4);
f(t ) = 100tu( t) − 200(t − 2 )u(t − 2) + 100(t − 4)u(t − 4);
- 42 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
s
−
100
s
e 2 2
−2 s
+
50 − 4 s
s2
e
50 100 − 2 s 50 − 4 s
s 2 y( s ) + 3sy( s) + 2 y(s ) = − 2 e + 2e ;
s2 s s
y( s )[s 2 + 3s + 2 ] = 2 − 2 e − 2 s + 2 e − 4 s ;
50 100 50
s s s
50 100 50
y( s ) = 2 2 − 2 2 e −2s + 2 2 e −4s ;
s (s + 3s + 2 ) s (s + 3s + 2 ) s (s + 3s + 2 )
50 100 50
y( s ) = 2 − 2 e −2s + 2 e − 4s ;
s (s + 2 )( s + 1) s (s + 2 )( s + 1) s (s + 2 )( s + 1)
1442443 14442444 3 144424443
y1( s ) y2( s ) y 3( s )
50
y 1 ( t ) = L −1 2
s (s + 2 )(s + 1)
50 −1 50 tu
Si L −1 = f ( t ) , entonces L 2 = ∫ ∫ f(θ)dθdu
(s + 2 )(s + 1) s (s + 2 )(s + 1) 0 0
50 −1 A B A(s + 1) + B(s + 2 )
L −1 =L + = L −1 ;
(s + 2 )(s + 1) (s + 2 ) (s + 1) (s + 2 )(s + 1)
A + B = 0
A + 2B = 50
Entonces:
50 tu t u
L −1 ∫∫
= f ( θ )dθ = ∫ ∫ (50e −θ − 50 −2 θ )dθdu ;
s 2
(s + 2 )( s + 1) 00 0 0
- 43 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
t t t
0 0 0
t
∫ [− 50e + 25 − 2 u + 25]du = [50e − u − 12.5e − 2 u + 25u]0 = [50e −t − 12.5e − 2 t + 25t − (50 − 12.5)];
−u t
0
t
Por lo tanto:
50
L −1 −t −2 t
= 50e − 12.5e + 25t − 37.5 ;
s2
(s + 2 )( s + 1 )
y 1 ( t ) = 50e − 12.5e − 2 t + 25t − 37.5 ;
−t
100 50
y 2 ( t ) = L −1 2 e − 2 s = 2L −1 2 e −2 s
s (s + 2 )(s + 1) s (s + 2 )(s + 1)
50
y 2 ( t ) = 2L −1 2 e − 2 s = 2 (50e −(t − 2 ) − 12.5e −2 ( t −2 ) + 25(t − 2 ) − 37.5)u( t − 2 );
s (s + 2 )(s + 1)
e − 4 s = (50e −(t − 4 ) − 12.5e −2 ( t − 4 ) + 25(t − 4 ) − 37.5 )u( t − 4);
50
y 3 ( t ) = L −1 2
s (s + 2 )(s + 1)
y( t ) = y 1 ( t ) + y 2 ( t ) + y 3 ( t );
y( t ) = 50e − t − 12.5e − 2 t + 25 t − 37.5 + 2 ((50e −( t − 2 ) − 12.5e − 2 ( t − 2 ) + 25(t − 2 ) − 37.5 )u( t − 2 ))
+ ((50e −( t − 4 ) − 12.5e − 2 ( t − 4 ) + 25(t − 4 ) − 37.5 )u( t − 4 ));
- 44 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
3) Una masa de 5kg se sujeta a un resorte suspendido del techo y ocasiona que el resorte se
estire 2 metros al llegar al reposo en equilibrio. Se eleva luego la masa 1 metro sobre el punto
de equilibrio y se le aplica una velocidad dirigida hacia arriba de 1/3 m/seg. Determine:
a) La ecuación del movimiento armónico simple de la masa.
π
b) La posición del objeto en t = segundos
4
d2 y dy
m 2
+C + Ky = f ( t )
dt dt
a)
Como no hay amortiguador C=0, además no existe fuerza perturbadora que se aplique al sistema por lo tanto
f(t)=0, la posición inicial de la masa es 1 metro sobre la posición de equilibrio por lo tanto si tomamos el eje de
referencia positivo hacia arriba la posición inicial de la masa sera 1 metro. Y la velocidad es 1/3 m/seg.
d2 y
5 2 + ky = 0 ;
dt
Se debe encontrar el valor de k:
Como la masa es 5kg y si se asume la gravedad 10m/seg 2 , el peso será de 50 Newton, al sujetar el resorte la
masa se estira 2 metros, lo que me indica de manera implícita la constante del resorte que se la puede calcular
mediante:
F = k∆l , donde F es el peso del objeto y ∆l la longitud del estiramiento. Despejando k se obtiene k=25N/m.
Para resolver esta ecuación diferencial aplicamos la transformada de Laplace a ambos lados de la ecuación:
d2 y
L 5 2 + 25 y = L [0];
dt
5s Y(s ) − 5sy(0 ) − 5 y' (0 ) + 25 y(s) = 0 ;
2
- 45 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
5
5s 2 y(s ) − 5s − + 25y(s) = 0 ;
3
(5s 2 + 25)y(s) = 5s + 5
3
(s 2 + 5)y(s) = s + 1
3
s 1
y(s ) = 2 +
(s + 5) 3(s + 5) ;
2
s 1 s 1
y( t ) = L-1 2 + = L-1 2 + L-1 2 ;
(s + 5 ) 3(s + 5) (s + 5) 3(s + 5)
2
1
y( t ) = cos 5t + sen 5(t )
15
b)
π / 4 segundos es:
La posición del objeto en
π π 1 π
y = cos 5 − sen 5
4 4 15 4
π 2 1 2 2 1 2 16 8 2
y = − − =− 1 + = − =−
4 2 15 2 2 15 2 15 15
- 46 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
1) Un circuito LRC con R=12 ohmios, L=1, C=0.01 faradios se conecta a una batería
que transmite un voltaje de 20 voltios. Si el interruptor esta inicialmente apagado y se lo
enciende después de 10 segundos, permaneciendo conectada por un lapso de 20
segundos y luego desconectada definitivamente. Si inicialmente no hay carga en el
condensador y la corriente inicial es cero, determine:
a) La carga acumulada en el condensador en los tiempos t=5s, y t=20s.
b) La intensidad de corriente que atraviesa el circuito en los tiempos t=8s, y t=40s.
20
10 30 t
1
LQ ' '+ RQ'+ Q = ε (t ) =20u(t-10)-20u(t-30)
C
1
l[LQ' ' ] + l[RQ' ] + l[ Q ] = l[ε( t )]
C
e −10 s − e −30 s
s 2 Q( s ) + 12sQ( s ) + 100Q( s ) = 20
s
e −10 s e −30 s
( s 2 + 12s + 100)Q( s ) = 20 −
s s
e −10 s e −30 s
Q( s ) = 20 2 −
s ( s + 12s + 100) s ( s + 12s + 100)
2
- 47 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
A Bs + C
+ 2 ⇒ As 2 + 12 As + 100 A + Bs 2 + Cs = 1
s s + 12s + 100
1
A=
100
1
B=−
100
12
C=−
100
1
100 − 1 s + 12
100 2
s s + 12 s + 100
Q(t ) = l −1 [Q( s )]
1 − 6 ( t − 10 ) 3
Q( t ) = 1− e cos 8( t − 10) − e −6( t −10 )sen 8( t − 10) U 10 ( t )
5 4
1 3
− 1 − e −6( t − 30 ) cos 8( t − 30) − e −6( t − 30 )sen 8( t − 30 ) U 30 ( t )
5 4
Cuando t=5s
Cuando t=20s
1 1 −6( t −10 ) 3
Q(t ) = − e cos 8(t − 10) − e −6 (t −10) sen8(t − 10)
5 5 20
1 1 −60 3
Q(20) = − e cos 80 − e −60 sen80
5 5 20
1 1 −60 3
Q(20) = − e (−0.110) − e −60 (−0.993)
5 5 20
- 48 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
2) Un circuito LRC con R=150 ohmios, L=1 Henrio, C=0.0002 faradios en t=0 se le
aplica un voltaje que crece linealmente de 0 a 100 voltios, durante 10 segundos, para
luego cesar por tiempo indefinido. Si inicialmente no hay carga en el condensador y la
corriente inicial es cero, determine:
a) La carga en cualquier instante de tiempo
b) La corriente del circuito en t=20s.
R = 150r Q ( 0) = 0
L = 1H Q ' ( 0) = 0
C = 2 × 10−4 F
V(t)
100
0 10 t
LQ ''+ RQ '+ 1/ C Q = V ( t )
Q ''+ 150Q '+ 5000Q = 10t ( µ 0 ( t ) − µ10 ( t ) )
Q ''+ 150Q '+ 5000Q = 10t µ 0 ( t ) − 10t µ10 ( t ) + 100 µ10 ( t ) − 100 µ10 ( t )
Q ''+ 150Q '+ 5000Q = 10t µ 0 ( t ) − 10 ( t − 10 ) µ10 ( t ) − 100 µ10 ( t )
Encontrando la transformada:
A = 1/ 500
B = −3 / 50000
10 A B C D
" = 2+ + + ⇒
s ( s + 50 )( s + 100 ) s
2
s s + 50 s + 100 C = 1/12500
D = −1/ 50000
"
- 49 -
Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial
Q(20segundos)=0
∂Q ( t )
i (t ) =
∂t
i(20segundos)=0
- 50 -
Roberto Cabrera V.