Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Repaso Laplace

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 44

ESCUELA SUPERIOR POLITÉCNICA DEL LITORAL

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

Este es un solucionario de problemas de Ecuaciones Diferenciales correspondiente a la Segunda


Evaluación, donde constan ejercicios tipo examen. Esta obra ha sido elaborada por Roberto Cabrera
y Christian de La Rosa, ex – estudiante de la ESPOL, con el auspicio de la directiva A.E.F.I.E.C. de los
años 2006, 2007, 2008. Modificado y corregido dos veces por Roberto Cabrera.
Ecuaciones Diferenciales – II Parcial

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

= 4 L{e }+ 6 L{t }− 3L{sen(4t )} + 2 L{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

= L{t e }+ L{4te }+ L{4e }+ L{e cosh (2t )}


2 t t t −4t

= L{t e }+ 4 L{te }+ 4 L{e }+ L{e cosh (2t )}


2 t t t −4t

Aplicando el primer teorema de la traslación:


{ } { } {} {
L t 2 e t + 4 L te t + 4 L e t + L e −4t cosh (2t ) }
s+4
L{t e }+ 4 L{te }+ 4 L{e }+ L{e cosh (2t )} =
− 4t 2! 1 1
2 t t t
+4 +4 +
(s − 1)3
(s − 1) 2
s − 1 (s + 4 )2 − 4
5s 4 + 29 s 3 + 9s 2 − 21s + 20
=
(s − 1)3 (s + 2)(s + 6)

-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

• Encuentre la transformada de la primera derivada de f(t)

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

pero lim e − sP f (P ) = 0 asumiendo que f (t ) es de orden exponencial


P →∞

= sF(s) - f(0)

• Encuentre la transformada de la función tf(t)

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

Derivando ambos lados de la igualdad tenemos :



d d
F (s ) = ∫ e − st f (t )dt
ds ds 0

∂ − st
=∫ e f (t )dt
0
∂s

= ∫ − te − st f (t )dt
0

= − ∫ e − st [tf (t )]dt
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

Por serie de potencias sabemos que sen ( t )= ∑


+∞
(− 1)n t 2
n = 0 (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

• Encuentre la transformada de la integral de f(t)

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

Entonces sabemos que :


L{g ' (t )} = sL{g (t )} − g (0)
t 
L{ f (t )} = sL ∫ f (u )du 
0 
Despejando tenemos que :
t  L{ f (t )} F (s )
L ∫ f (u )du  = =
0  s s

• Encuentre la transformada f(t)/t

 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

• Demuestre el segundo teorema de la traslación


Si L{ f (t )} = F (s ), entonces L{µ (t - a ) f (t − a )} = e − as F (s )

Tenemos : L{µ (t - a ) f (t − a )} = ∫ e − st µ (t - a ) f (t − a )dt
0

Entonces : L{µ (t - a ) f (t − a )} = ∫ e − st f (t − a )dt
a

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

e − t 2 ;2 n < t < 2 n + 1


• Encuentre la transformada f (t ) =  n = 0,1,2,3,....
 0 ;2n + 1 < t < 2n + 2

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

• Encuentre la transformada de la siguiente gráfica

Tenemos que encontrar la transformada de una función periódica:


sen(t ) 0 < t < π
g (t ) =  extendida periodicamente con periodo 2π
 0 π < t < 2π

1
L{g (t )} = ∫e
− st
g (t )dt
1 − e −2πs 0
π
1
L{g (t )} =
1 − e −2πs ∫0
e − st sen(t )dt

Integro por partes : u = e - st → du = − s e - st


dv = sen(t )dt → v = − cos(t )

∫e sen(t )dt = − cos(t )e - st − s ∫ e - st cos(t )dt


− st

Integro por partes : u = e - st → du = − s e - st


dv = cos(t )dt → v = sen(t )

∫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

• Demuestre el teorema de la convolución

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

La región en el plano en donde se llevará a cabo la integración es:

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

Resolución de ecuaciones diferenciales mediante las transformada de


Laplace

Encuentre la solución de la siguiente ecuación diferencial:

• y ' ' '+4 y ' '+5 y '+2 y = 10 cos(t ), y (0) = y ' (0) = 0 y ' ' (0) = 3

Aplicando la transformada de Laplace


L{y ' ' '} + 4 L{y ' '} + 5L{y '} + 2 L{y} = 10 L{cos(t )}
Encuentro las transformadas necesarias :
L{y ' ' '} = s 3Y ( s ) − s 2 y (0) − sy ' (0) − y ' ' (0) = s 3Y ( s ) − 3
L{y ' '} = s 2Y ( s ) − sy (0) − y ' (0) = s 2Y ( s)
L{y '} = sY ( s) − y (0) = sY ( s )
L{y} = Y ( s )
s
L{cos(t )} = 2
s +1
Re emplazando las transformadas :

(s Y (s) − 3)+ 4(s Y (s))+ 5(sY (s)) + 2Y (s) = 10 s s+ 1


3 2
2

(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

3s 2 + 10s + 3 = ( A + B + D )s 4 + (2A + 3B + C + 4D + E )s 3 + (2A + 3B + 2C + 5D + 4E )s 2 + (2A + 3B + C + 2D + 5E )s + ( A + 2B + 2C + 2E )


Tenemos el siguiente sistema de ecuaciones :
A+ B + D = 0
2A + 3B + C + 4D + E = 0
2A + 3B + 2C + 5D + 4E = 3
2A + 3B + C + 2D + 5E = 10
A + 2B + 2C + 2E = 3
De donde A = -1, B = 2, C = -2, D = -1, E = 2
−1 2 2 −s+2
Y ( s) = + − +
s + 2 s + 1 (s + 1)2 s 2 + 1
 − 1 2 2 − s + 2 
y (t ) = L−1{Y ( s )} = L−1  + − + 2 
 s + 2 s + 1 (s + 1)2
s + 1 
y (t ) = −e − 2t + 2e − t − 2te − t − cos(t ) + 2sen(t )

- 20 -

Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial

Encuentre la solución de la siguiente ecuación diferencial:


d2y − 4t + 8π ;0 < t < 2π
• + 4 y = h(t ), donde h(t ) =  , y (0 ) = 2, y ' (0) = 0
dt 2
 0 ; t > 2π
Aplicando la transforma da de Laplace
L{y ' '} + 4 L{y} = L{h (t )}
Encuentro las transforma das necesarias :
L{y ' '} = s 2Y ( s ) − sy (0 ) − y ' (0 ) = s 2Y ( s ) − 2 s
L{y} = Y ( s )
L{h (t )} = L{(µ 0 (t ) − µ 2π (t ) )(− 4t + 8π )} = L{(µ 0 (t ) )(− 4t + 8π )} + 4 L{(µ 2π (t ) )(t − 2π )}
− 4 8π 1
L{h (t )} = 2
+ + 4e − 2πs 2
s s s
Re emplazando :
− 4 8π 1
s 2Y ( s ) − 2 s + 4Y ( s ) = 2 + + 4 e − 2 πs 2
s s s

(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

• Determinar la solución del siguiente problema de valor inicial:

Primero se expresa en términos de funciones escalones de la siguiente manera:

Se reemplaza en la ecuación diferencial y se procede a resolverla usando transformadas de Laplace:

Despejando Y(S):

Encontrando la solución mediante transformada inversa de Laplace:

i)

ii)

iii) Entonces

iv)

- 22 -

Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial

Encuentre la solución de la siguiente ecuación integro - diferencial:


t
t3
• ∫ y (u ) y (t − u )du = 2 y (t ) + − δ (t )
0
6
Aplicando la transforma da de Laplace
 t  t3 
L  y (u ) y (t − u ) du  = 2 L{y (t )} + L   − L{δ (t ) }

 0  6 
Encuentro las transforma das necesarias :
 t 
L  y (u ) y (t − u ) du  = L{y (t ) * y (t )} = Y 2 ( s )

 0 
L{y (t )} = Y ( s )
t3  3! 1
L  = 4 = 4
6
  6 s s
L{δ (t ) } = 1
Re emplazando :
1
Y 2 ( s ) = 2Y ( s ) + −1
s4
 s4 −1
Y 2 ( s ) − 2Y ( s ) +  4  = 0
 s 
 s4 −1 4s 4 − 4s 4 + 4
2 ± 4 − 4 4  2 ±
 s  s4
Y1, 2 ( s ) = =
2 2
1
Y1 ( s ) = 1 + 2 → y1 (t ) = δ (t ) + t
s
1
Y2 ( s ) = 1 − 2 → y 2 (t ) = δ (t ) − t
s

- 23 -

Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial

Encuentre la solución de la siguiente ecuación diferencial de coeficientes variables:


• ty ' '+ (1 − 2t ) y '−2 y = 0, y (0) = 1, y ' (0) = 2
Aplicando la transforma da de Laplace
L{ty ' '} + L{(1 − 2t ) y '} − 2 L{y} = 0
Encuentro las transforma das necesarias :

L{ty ' '} = −


d
ds
L{y ' '} = −
d 2
ds
[ ]
s Y ( s ) − sy (0 ) − y ' ( 0) = − s 2Y ' ( s ) − 2 sY ( s ) + 1

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

Encuentre la solución de la siguiente ecuación diferencial de coeficientes variables:


• ty ' '−(t + 2 )y '+3 y = t − 1
Aplicando la transforma da de Laplace
L{ty ' '} − L{(t + 2 ) y '} + 3 L{y} = L{t } − L{1}
Encuentro las transforma das necesarias :

L{ty ' '} = −


d
ds
L{y ' '} = −
d 2
ds
[ ]
s Y ( s ) − sy ( 0) − y ' ( 0 ) = − s 2Y ' ( s ) − 2 sY ( s ) + k1

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 Y ' ( s ) − 2 sY ( s ) + k ) − (− sY ' ( s ) + (2 s − 1)Y ( s ) − 2k ) + 3Y ( s ) = 1s− s


2
1 1 2

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

1) Usando el método de eliminación, resuelva el siguiente sistema de ecuaciones


diferenciales:

2 x'+ y'− x − y = e − t ( 1)



x'+ y'+2 x + y = e t (2)
Restando: (1)-(2);

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

x' = −C 1 sent + C 2 cos t

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

De la primera ecuación despejamos y;


1 1 dx
⇒ y= ( 2x )−
3 3 dt
2 x´
⇒ y= x−
3 3
Reemplazando y en la segunda ecuación:
dy 2 x´´
⇒ = x´ −
dt 3 3
2 x´´ 2 x´
⇒ x´ − = x− − 2x
3 3 3 3

Multiplicando la ecuación por 3;

⇒ 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

* Encuentre la solución particular del problema anterior dado:


x (0)=8, y (0)=3

Del ejercicio anterior:

x = C1e 4t + C 2 e − t
2
y = − C1e 4 t + C 2 e −t
3

Como x (0)=8, entonces:


8= C1+C2 1
Como y (0)=3, entonces:
2
3 = − C1 + C 2 2
3

Con 1 y 2 se obtiene un sistema de 2 ecuaciones con dos incógnitas; resolviendo el


sistema se obtiene:
C2=5, C1=3

⇒ La solución particular es: x = 3e 4t + 5e − t


y = −2e 4 t + 5e − t

- 28 -

Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial

Método de los operadores diferenciales

1) Usando el método de las operaciones diferenciales resuelva el siguiente sistema de


ecuaciones diferenciales:
 (D 2 − 4D + 4)x 1 + (D 2 + 2D )x 2 = t
 2
(D − 2D )x 1 + (D + 4D + 4)x 2 = e
2 t

 (D − 2 ) 2 x 1 + D(D + 2 )x 2 = t

D(D − 2 )x 1 + (D + 2 ) x 2 = e
2 t

Encontrando x 1 ( t ) usando la regla de Kramer se obtiene que:

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

Encontrando la solución particular x p 1:


x p 1 = a + bt + ce t ;
x'p 1 = b + ce t ;
x''p 1 = ce t ;
3 t
Reemplazando en x1''(t) − 4x1 (t) = −1 − t + e , se obtiene :
4
ce t − 4(a + bt + ce t ) = −1 − t + e t ;
3
4
3
− 4a − 4 bt − 3ce t = −1 − t + e t ;
4

- 29 -

Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial

Ahora se procede a encontrar la solución x 2 ( t ), usando la regla de Kramer :


(D − 2 ) 2 t
D( D − 2 ) e t (D − 2 ) 2 e t − D( D − 2 )t (D − 2 )[( D − 2 )e t − Dt ] (D − 2 )(e t − 2e t − 1)
x 2 (t ) = = = =
(D − 2 ) 2 D( D + 2 ) − 4( D 2 − 4 ) − 4( D 2 − 4 ) − 4(D 2 − 4)
D(D − 2 ) (D + 2 ) 2
(D − 2 )(− e t − 1) − e t + 2e t + 2 et + 2
x 2 (t ) = = = ;
− 4( D 2 − 4 ) − 4(D 2 − 4) − 4(D 2 − 4 )
− 4(D 2 − 4)x 2 ( t ) = e t + 2 ;
et 1
( D 2 − 4 )x 2 ( t ) = − − ;
4 2
et 1
x 2 ' ' ( t ) − 4x 2 ( t ) = − − ;
4 2

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

2.-) Usando el método de los operadores diferenciales resuelva el sistema:

 ( D + 2 )( x 1 ) + (D − 1)( x 2 ) = −sent ( 1)

(D − 3)( x 1 ) + (D + 2 )( x 2 ) = 4 cos t ( 2 )

Multiplico 1 por (D − 3) ∧ 2 por (D + 2 )

(D + 2 )(D − 3)x 1 + (D − 1)(D − 3)x 2 = (D − 3)( −sent )



(D + 2 )(D − 3)x 1 + (D + 2 ) 2 x 2 = (D + 2 )( 4 cot)

(D 2 − 4D + 3)x 2 − (D 2 + 4D + 4 )x 2 = − cos t + 3sent + tsent − 8 cos t


( −8D − 1)x 2 = 7 sent − 9 cos t
− 8x 2 ' − x 2 = 7 sent − 9 cos t
8x 2 ' + x 2 = 9 cos t − 7 sent ;
8x 2 ' + x 2 = 0 ;
x 2 = e rt ;
x 2 ' = re rt ;
8r + 1 = 0 ;
1
r=− ;
8
1
− t
x 2 = Ce 8 ;
x 2 = A cos t + Bsent ;
x 2 ' = −Asent + B cos t ;
8x 2 ' + x 2 = 0 ;
8(− Asent + B cos t ) + A cos t + Bsent = 0 ;
(− 8A + B )sent + (8B + A ) cos t = 9 cos t − 7 sent ;
− 8A + B = −7 ;

8B + A = 9 ;
Resolviendo el sistema :
A = 1, B = 1,
La solución particular es :
x 2 p = cos t + sent ;
1
− t
x 2 = Ce 8 + cos t + sent ;

- 31 -

Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial

Ahora procedemos a encontrar x 1 del sistema de ecuaciones:


 ( D + 2 )( x 1 ) + (D − 1)( x 2 ) = −sent ( 1)

(D − 3)( x 1 ) + (D + 2 )( x 2 ) = 4 cos t ( 2 )
x 1 '+2 x 1 + x 2 '−x 2 = −sent ; (1)
x 1 '−3x 1 + x 2 '+2 x 2 = 4 cos t (2)
Restando (1) y (2), se obtiene :
5x 1 − 3x 2 = −sent − 4 cos t ;
x 1 = 3x 2 − sent − 4 cos t;
 − 18 t 
x 1 = 3 Ce + cos t + sent  − sent − 4 cos t;
 
1
− t
x 1 = 3Ce 8
− cos t + 2sent ;
La solución del sistema es :
 1
− t
x
 1 = 3 Ce 8
− cos t + 2 sent ;
 1
x = Ce − 8 t + cos t + sent ;
 2

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

 x'−3x'+2 y = sent; x(0) = 0; 1



4 x − y '− y = cos t ; y (0) = 0; 2
Aplicando transformada de Laplace a las dos ecuaciones:
L [x'] -3 L [x ] +2 L [ y ] = L [sent ]
L [4 x] - L [ y '] - L [ y ] = L [cos 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

Sumo 1 y 2, entonces se obtiene:


s 2 − 3s − 4
− [8 + ( s − 3)( s + 1)]y ( s ) =
s2 +1
[ ]
− s 2 − 2s + 5 y( s) =
( s − 4)( s + 1)
s2 + 1
( s 2 − 3s − 4)  As + B Cs + D 
⇒ y (s) = − 2 = − 2 + 2
( s − 2s + 5)( s + 1)  s − 2s + 5 s + 1 
2

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

De la ecuación 4x-y’-y=cos(t); podemos encontrar x(t):


y ′ + y + cos( t )
x( t ) =
4
3 −t 11 11 7
y ′( t ) = − [e cos(2 t ) − e − t sen (2 t )] − [ −2 e − t sen( 2 t ) − e − t cos(2 t )] − sen( t ) + cos(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

X' '+ Y'+3X = 15e − t


2) Resolver  con las condiciones X(0)=0, X’(0)=0, Y(0)=0, Y’(0)=0.
Y' '−4 X'+3Y = 15t
Aplicando la transformada de Laplace a ambas ecuaciones:
[
£ [x''+ y'−x] = £ 15e − t ]

£ [y''−4 x'− y ] = £ [15t ]

[  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

Obteniendo x(t) aplicando transformada de Laplace inversa a X(S) :


− 1 2s − 1 s  -1  − 1 
 2s − 1  -1  s 
x(t) = 15£ - 1  + 2 2 +  = 15£   + 15£  2
-1
2  + 15£  2 ;
 s (s + 1) (s 2
+ 1 )   s  (s + 1)   + 1 )
(s
 − 1
£ - 1   = −1;
 s 
 s 
£ -1  2  = cos t ;
 (s + 1)
 2s − 1   s   1 
£ -1  2 2
= 2£ - 1  2 2
− £ -1  2 2
 (s + 1)   (s + 1)   (s + 1) 

- 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

Método de los valores y vectores propios

1) Resuelva por el método de los valores y vectores propios el siguiente sistema:


x ' = −4 x + y + z
y' = x + 5 y − z
z' = y − 3z
 x'   − 4 1 1  x 
    
 y'  =  1 5 − 1  y 
 z'   0 1 − 3  z 
    
− 4 1 1 
  det(A-λI)=0
A =  1 5 − 1
 0 1 − 3
 

−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  − 1e + 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 − λ 

(−3 − λ )[λ (λ + 1)] + 2[−1 − 2(1 + λ )] = 0


− (λ + 3)(λ )(λ + 1) + 2(−3 − 2λ ) = 0
− λ3 − 4λ2 − 3λ − 6 − 4λ = 0
λ 3 + 4λ 2 + 7 λ + 6 = 0
λ1 = −2 λ 2 = −1 + 2 i λ 3 = −1 − 2i

Se procede a encontrar el vector propio asociado al siguiente valor:


* λ1 = −2
 −1 0 2 0  −1 0 2 0  −1 0 2 0
     
 1 1 0 0 ≈  0 1 2 0 ≈  0 1 2 0
 − 2 −1 2 0  0 −1 − 2 0  0 0 0 0
     
y= -2z
x= 2z
 2 
 
v1 =  − 2 
 1 
 
Se procede a encontrar el vector propio asociado al siguiente valor λ 3 = −1 − 2i :
 − 2 + 2i 0 2 0  6 0 2 ( −2 − 2 i ) 0  6 0 2( − 2 − 2 i ) 0 
     
 1 2i 0 0 ≈  1 2i 0 0 ≈ 1 2i 0 0 ≈
   
1 + 2 i 0  − 2 −1 1 + 2 0   
 −2 −1 −3 − 1 + 2 0 
0

- 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

Entonces procedemos a encontrar la primera solución l.i. con λ1 = −2 :


 2 
−2 t
 
x1 = e  − 2 
 1 
 

Ahora procedemos a encontrar la segunda y tercera solución l.i. con λ 3 = −1 − 2i , tiene la


siguiente forma λ = α + βi , por lo tanto las otras dos soluciones son:
()
x 2 = e αt sen (β t ) a + e αt cos(βt ) b ; ()
x 3 = e αt cos(β t )(a ) − e αt
sen (βt )(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

Aplicaciones de Sistema: Masa – Resorte – Amortiguador

1) Una masa de 1 kilogramo sujeta a un resorte con una constante k = 9 m/seg


se suelta del reposo 1 metro debajo de la posición de equilibrio del sistema masa-
resorte, y empieza a vibrar. Después de π / 2 segundos, la masa es golpeada hacia
arriba por un martillo que ejerce un impulso de 3 newtons.

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

Como no hay amortiguador C=0;


En t = π / 2 segundos hay un impulso hacia arriba de 3 Newtons, por lo tanto hay una perturbación
 π
f( t ) = −3δ t −  , el signo negativo se debe a que tomamos el eje de referencia positivo hacia abajo.
 2

La ecuación diferencial que representa al sistema es:

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

2) Un sistema vibratorio compuesto de un resorte de constante k = 4N / m , un


amortiguador de c = 6Ns / m , tiene adherido una bola metálica de 20 Newton de peso.
Determine la forma en que vibra la masa si inicialmente esta en la posición de equilibrio
y sin velocidad inicial, y si desde el tiempo t = 0 actúa una fuerza perturbadora definida
así:

100t ; t ∈ [0 ,2 )
f( t ) = 
400 − 100t ; t ∈ (2 , 4 ]

La ecuación diferencial que representa al sistema es:

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

La ecuación diferencial queda expresada de la siguiente forma:


d2y dy
2 2 +6 + 4 y = 100 tu( t ) − 200(t − 2 )u( t − 2 ) + 100(t − 4 )u( t − 4);
dt dt

- 42 -

Roberto Cabrera V.
Ecuaciones Diferenciales – II Parcial

Ahora se puede proceder a resolver la ecuación diferencial mediante la transformada de Laplace:


 d2 y dy 
L 2 2
+6 + 4 y  = L [100 tu( t ) − 200(t − 2 )u( t − 2 ) + 100(t − 4 )u( t − 4 )];
 dt dt 
 d2 y dy 
L 2
+3 + 2 y  = L [50 tu( t ) − 100(t − 2 )u( t − 2 ) + 50(t − 4 )u( t − 4 )];
 dt dt 

La posición inicial del sistema es y(0)=0 metro, y la velocidad inicial es y’(0)=0:

[s y(s) − sy(0) − y'(0) + 3sy(s) − 3y(0) + 2 y(s)] = 50


2

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

Para encontrar y 1 ( t ) , se procede a usar el teorema de la integral de la transformada de Laplace:

 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

Resolviendo el sistema de ecuaciones, se obtiene:


B = 501, A = -50;
 50  −1  50 50 
L −1   =L  −  = 50e −t − 50−2 t ;
 (s + 2)(s + 1)   (s + 1) (s + 2) 

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

∫ [− 50e + 25 − 2 θ ]0 du = ∫ [− 50e − u + 25 − 2 u − (− 50 + 25 )]0 du = ∫ [− 50e − u + 25 −2 u + 25]du ;


−θ u u

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

∫ [− 50e + 25− 2 u + 25]du = 50e − t − 12.5e − 2 t + 25t − 37.5;


−u

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

Ahora y(t) es:

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

Se puede representar y(t) en como una función con regla de correspondencia:

50e − t − 12.5e −2 t + 25t − 37.5 0 ≤ t < 2;



y( t ) = 50e −t (1 + 2 e 2 ) − 12.5e − 2 t (1 + 2e 4 ) + 75t − 212.5 2 ≤ t < 4;

50e (1 + 2 e + e ) − 12.5e (1 + 2e + e ) + 100t − 350 ;


 −t −2 t
2 4 4 8
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.

La ecuación diferencial que representa al sistema es:

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

La posición inicial del sistema es y(0)=1 metro, y la velocidad inicial es y’(0)=1/3:


Reemplazando las condiciones se obtiene:

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

Aplicaciones de Circuitos Eléctricos

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 

1  10 s  1 s+6 6  −30 s  1 s+6 6 


Q(s ) = e  s − −  − e  − − 
s ( s + 6) + 64 ( s + 6) + 64 
5  ( s + 6) + 64 ( s + 6) + 64 
2 2

2 2


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

Q(5) = 0 Condensador descargado

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

Q(20) = 2.08 x10 −25 coulombs

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

10 10e −10 s 100e−10 s


s 2Q ( s ) − sQ ( 0 ) − Q ' ( 0 ) + 150 s Q ( s ) − 150 Q ( 0 ) + 5000Q ( s ) = − −
s2 s2 s
10 −10 s 10 −10 s 100
(s 2
+ 150 s + 5000 ) Q ( s ) =
s2
−e
s2
−e
s
10 10 100
( s + 50 )( s + 100 ) Q ( s ) = 2 − e−10 s 2 − e−10 s
s s s

 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.

También podría gustarte