Semana 02
Curso: Ecuaciones Diferenciales Parciales
Tema: La transformada de Laplace y las ecuaciones diferenciales en derivadas parciales
Un capítulo importante de la transformada viene dado por su uso en la solución de ecuaciones
en derivadas parciales. Por ejemplo, supongamos que
T ( x ,t )=x +t
Y decidimos extraer la transformada respecto a la variable t, con x fijo. Tenemos
x 1
L ( T ( x , t ) ) ( s )= + 2
s s
Al extraer la transformada respecto a t se hace desaparecer justamente la variable t, i.e
L(T(x , t)) (s) resulta ser solo función de x y del parámetro s , que en adelante
denotaremos T́ s ( x ) . Por ello
∂ d
L ( T ( x , t ) )es en verdad L ( T ( x , t ) ) , de ahí que en el ejemplo en discusión tengamos
∂x dx
d 1
T́ ( x )=
dx s s
En general, suponiendo que la transformada se extrae respecto a t ,
∞
L ( T ( x , t ) ) (s)=∫ e−st T ( x , t ) dt .
0
Aceptando el intercambio entre derivada parcial e integral indefinida tendremos
∞
∂ ∂
L ( T ( x , t ) ) ( s )=∫ e−st T ( x ,t ) dt .
∂x 0 ∂ x
d ∂
i.e.
dx (
L ( T ( x , t ) ) ( s ) =L
∂x )
T ( x , t ) (s)
(1)
De otro lado, evidentemente
L ( ∂∂tT ( x , t ) ) ( s )=sL(T ( x , t )) ( s) −T ( x , 0)
Si se hubiera tomado transformada respecto a la primera variable, hubiéramos llegado a
L( ∂∂x T ( x , t )) ( s )=s L (T ( x , t ) ) ( s )−T (0 , t)
∂ d
L ( T ( x , t )) ( s )= L ( T ( x , t ) ) ( s )
∂t dt
En forma análoga, para derivadas de 2º orden
∂2 T d2
L (
∂ x2 )
( x , t ) ( s) = 2 L ( T ( x , t ) ) ( s )
dx
∂2 T ∂T
L ( ∂t 2 )
( x , t ) ( s ) =s 2 L ( T ( x ,t ) ) ( s )−s T ( x , 0 )−
∂t
(x , 0)
cuando se extrae transformada respecto a la variable t, mientras que
∂2 T ∂T
L ( ∂x 2 )
( x , t ) ( s ) =s 2 L ( T ( x ,t ) ) ( s )−s T ( 0 , t )−
∂x
(0 , t)
∂2 T d2
L ( ∂t 2
(x ,t ) )( )
s=
dt 2
L ( T ( x ,t ) ) ( s ) (2)
cuando se extrae transformada respecto a la variable x. Notar la gran simetría existente. Las
formulas (1) . (2) constituyen la clave del uso práctico de la transformada en relación a las EDP
Semana 02
En resumen
∂u
1) L [ ]
∂t
=s u ( x , s )−u(x , 0)
∂ u du
2) L [ ] =
∂ x dx
∂2 u
3) L
∂t 2[ ]
=s 2 u ( x , s ) −s u ( x , 0 )−ut ( x , 0)
∂2 u d 2 u
4) L
∂ t2[ ]
= 2
dx
∂ u(x ,0)
Nota: ut ( x , 0 )=
∂t
du
=¿ Derivada (ordinaria) de u respecto a x
dt
∂u
=u x =¿ Derivada parcial de u respecto a x
∂x
Problemas Resueltos
Prob(1) Halle la solución de la ecuación diferencial
∂u ∂u
− =1−e−t
∂x ∂t
u ( x , 0 ) =x
u ( x , t )acotada ; 0<x<1 ; t>0
Solución
∂u ∂u
L [ ] [ ]
∂x
−L
∂x
= L [ 1 ] −L [ e−t ]
du 1 1
−( s u−u ( x , 0 ) )= −
dx s s+ 1
1 1
u' −s u+ x= −
s s+ 1
1 1
u' −s u=−x+ −
s s +1
1 1
→ u=k esx +e sx ∫ e−sx −x+ −
⏟ s s +1[dx ] (α )
I
x 1
dx= e− sx + 2 e−sx
−sx
pero *) ∫ −x e
s s
1 − sx −1 −sx
*) ∫ e dx= e
s s2
1 −sx 1
*) −∫ e dx= e−sx
s +1 s( s+ 1)
x 1 1 1
→ I = e− sx + 2 e−sx − 2 e−sx + e− sx
s s s s(s+1)
x −sx 1 −sx 1 − sx
¿ e + e − e
s s s +1
En ( α )
x −sx 1 −sx 1 −sx
u=k e sx + esx [ s
e + e −
s s +1
e ]
x 1 1
¿ k e sx + + −
s s s +1
x 1 1
L−1 [ u ] =k L−1 [ e sx ] + L−1 [] [ ]
s
+ L−1 −
s s
U ( x , t )=k L−1 [ e sx ] + x +1−e−t
Además L−1 [ e sx ]=e sx δ (t)
U ( x , t )=k e sx δ ( t ) + x +1−e−t
U ( x , 0 ) =k e sx δ ( 0 )+ x+1−1=x
→ k=0
∴ U ( x , t ) =x+1−e−t ∎
Recordar: L [ δ ( t ) ] =1
L [ δ ( t−t 0 ) ] =e−s t ↔ L−1 [ e−s t ] =δ (t−t 0)
0 0
δ ( t ) : Delta de Dirac
Prob(2) Halle la solución de
∂U ∂U
=2 +U
∂x ∂t
U ( x , 0 ) =6 e−3 x
U acotada para x>0 , t>0
Solución
∂U ∂U
L [ ] [ ]
∂x
=2 L
∂t
+ L[U ]
du
=2 [ su−U (x , 0) ] +u
dx
du
−u−2 su=−2 U (x , 0)
dx
du
−( 1+2 s ) u=−12 e−3 x
dx
→ u=k e(1 +2 s) x + e(1 +2 s) x∫ e−(1+ 2 s) x (−12 e ¿¿−3 x )dx ¿
¿ k e (1+ 2 s) x −12 e (1+2 s) x ∫ e−(4 +2 s) x dx
1
¿ k e (1+ 2 s) x −12 e (1+2 s) x . e− (4 +2 s ) x
( 4 +2 s )
6 −3 x
u=k e (1+2 s ) x + 6
s+2
1
L−1 [ u ( x , t ) ] =L−1 [ k e (1 +2 s) x ] +6 L−1 [ ]
s+ 2
e−3 x
6 −3 x
U ( x , t )=k L−1 [ e(1+2 s )x ]+ e
s+2
U ( x , 0 ) =k L−1 [ e( 1+2 s ) x ] +6 e−2t e−3 x
6 e−3 x =k L−1 [ e(1 +2 s) x ] +6 e−3 x
k¿
→ U ( x , t ) =6 e−(2 t +3 x) ∎
Prob(3)
Resolver
∂U ∂2 U
=2
∂t ∂ x2
U ( 0 , t ) =0 ; U (5 , t )=0
U ( x , 0 ) =10 sen ( 4 πx )
Solución
∂U ∂2 U
L [ ] [ ]
∂t
=2 L
∂ x2
d2u d2 u
su−U ( x , 0 )=2 → s u−10 sen 4 π x=2
dx 2 dx 2
d2u
→2 −s u=−10 sen 4 π x
dx2
i.e 2 u' ' −s u=−10 sen 4 π x (M)
Cálculo de uh
s s
Sea p ( r )=2 r 2−s=0 →r 2= →r =±
2 √ 2
s s
→ uh=c 1 e
−
√ 2
x −
+ c2 e √ 2
x
Cálculo de u p
u p =A sen ( 4 πx ) + B cos (4 πx)
u p ' =4 Aπ cos ( 4 πx )−4 Bπ sen( 4 πx)
u p ' ' =−16 A π 2 sen ( 4 πx )−16 B π 2 cos (4 πx )
En <M>
−32 A π 2 sen ( 4 πx )−32 B π 2 cos ( 4 πx )−s ( Asen ( 4 πx )+ Bcos ( 4 πx ) )=−10 sen( 4 πx)
*) (−32 A π 2−sA ) sen ( 4 πx ) + (−32 B π 2−sB ) cos 4 πx=−10 sen (4 πx )
10
↔−32 A π 2−sA =−10 → A=
s+ 32 π 2
*) −32 B π 2−sB=0 → B=0
10
∴ u p= sen( 4 πx)
s +32 π 2
s s
u=c e √ +c e√
−
2
x
2
x 10
1 2 + sen(4 πx)
s+ 32 π 2
Transformando las condiciones
*) U ( 0 , t ) =0 → L ¿
10
→ 0=c1 +c 2 + → c 1+ c 2=0 → c2 =−c 1
s+ 32 π 2
*) U ( 5 , s )=0 → L [ U (5 , t) ] =L [ 0 ] → u ( 5 , s )=0
s s s s
→ 0=c1 e
5
√ +c e √ =¿ c ( e √ −e √ )=0 → c =0
2
2
−5
2
1
5
2
−5
2
1
c 2=0
En (M)
10
u ( x , t )= 2
sen 4 πx
s +32 π
1
→ L−1 [ u( x , t) ] =10 sen ( 4 πx ) L−1
[ s+ 32 π 2 ]
2
U ( x , t )=10 e−32 π t sen 4 πx ∎
Prob(4) Resolver
∂U ∂2 U
=3
∂t ∂ x2
U ( 0 , t ) =0 ,U ( π , t ) =0
U ( x , 0 ) =30 sen 5 x
Solución
∂U ∂2 U
L [ ] [ ]
∂t
=3 L
∂ x2
d2u
s u−U ( x , 0 )=3 → s u−30 sen 5 x=3 u' '
dx 2
→ 3u ' ' −su=−30 sen 5 x <M>
Cálculo de uh
s s
Sea p ( r )=3 r 2−s=0 → r 2= → r=±
3 √ 3
s s
u =c e √ +c e√
− x x
3 3
h 1 2
Cálculo de u p
u p =A sen 5 x+ Bcos 5 x
u'p =5 A cos 5 x−5 B sen 5 x
u'p' =−25 A sen 5 x−25 B cos 5 x
En <M>
−75 A sen 5 x−75 B cos 5 x−s ( A sen 5 x+ Bcos 5 x )=−30 sen 5 x
(−75 A−sA ) sen 5 x+ (−75 B−sB ) cos x=−30 sen 5 x
30
↔−75 A−sA =−30→ A=
75+ s
−75 B−sB=0→ B (75+ s )=0 → B=0
30
up= sen 5 x
s+75
s s
u=c e √ +c e√
−
2
x
2
x 30
1 2 + sen 5 x
s+ 75
Sabemos que L [ U ( 0 , t ) ] =L[0]
u ( 0 , s )=0
→ u ( 0 , s ) =c 1+ c 2=0 → c1=−c 2
De otro lado L [ U ( π , t ) ] = L[0 ]
u ( π , s ) =0
s s
→ u ( π , s )=c 1 e
−
√ + c e√
2
2
2
π
=0
s s
→ c1 e( √ −
2
π
+e √ )=0 →c =0
2
π
1
c 2=0
30
→ u p= sen 5 x →
s +75
1
U ( x , t )=30 L−1 [ ] s+75
sen 5 x
U ( x , t )=30 ( sen 5 x ) e−75t ∎
Prob(5)
Resolver la ecuación diferencial parcial
∂2 U ∂ 2 U
{
=
∂ x 2 ∂t 2
U ( 0 ,t )=0 ,U ( 1 ,t )=0
U ( x ,0 )=sen 2 π x
∂U
( x , 0 )=πsen 3 π x ; 0< x <1
∂t
Solución
∂2 U ∂2 U
1) L
[ ] [ ]
∂ x2
= L
∂ t2
Haciendo L [ U ( x , t ) ] =u (s ,t )
d2u 2 ∂U
2) 2
=s u−sU ( x , 0 )− (x , 0)
dx ∂t
¿ s2 u−s sen 2 πx−πsen 3 πx
→ u' ' −s 2 u=−s sen 2 πx −π sen 3 πx (α )
como U ( 0 , t ) =0 → L [ U ( 0 ,t ) ] =L [ 0 ] → u ( 0 , s )=0
U ( 1 ,t )=0 → L [ U ( 1 ,t ) ]=L [ 0 ] → u ( 1 , s ) =0
Cálculo de uh
Sea p ( r )=r 2 −s 2=0 → r=± s
→ uh ( s , t )=c 1 e−sx +c 2 e sx
uh ( 0 , t ) =c 1+ c 2=0 → c1=−c 2
uh ( s ,t )=c1 e−s +c 2 e s=0
→ un ( 1 , t )=c 1 e− s−c 1 e s=0 → c1 (e ¿ ¿−s−e s )=0 ¿
→ c 1=0
∴ c 2=0
Cálculo de u p
u p =Acos 2 πx+ Bsen 2 πx + M sen 3 πx + Ncos 3 πx
u p ' =−2 πA sen 2 πx +2 πB cos 2 πx+ 3 Mπ cos 3 πx−3 πN sen 3 πx
u p ' ' =−4 π 2 A cos 2 π x−4 π 2 Bsen 2 πx−9 M π 2 sen 3 πx−9 N π 2 cos 3 πx
En (α )
¿
−s ¿
¿−s sen 2 πx −π sen 3 π x
→ (−4 π 2 A−sA ) cos 2 πx + (−4 π 2 B−sB ) sen 2 πx
+ (−9 M π 2−sM ) sen 3 πx+ ( −9 N π 2−sN ) cos 3 πx=−s sen 2 πx−π sen 3 πx
− A ( 4 π 2+ s )=0 → A=0
s
−B ( 4 π 2 + s ) =−s → B=
s+ 4 π 2
π
−M ( 9 π 2+ s )=−π → M =
s+9 π 2
−N ( 9 π 2+ s )=0 → N=0
s π
up= 2
sen 2 πx+ sen 3 π x
s+ 4 π s+ 9 π 2
s π
L
−1
[ u p ]=sen 2 πx L−1 [ s+4 π 2]+ sen 3 πx L
[
−1
s+9 π 2 ]
1
U ( x , t )=cos 2 πt sen 2 πx+ sen 3 πt sen 3 πx
3
Prob(5) Resolver la ecuación diferencial
∂2U ∂2 U
{
−4 +U=16 x +20 sen x
∂ t2 ∂ x2
U ( 0 , t )=0 , U ( π , t )=16 π
∂U
( x , 0 )=0 , U ( x , 0 )=16 x
∂t
Solución
∂2 U ∂2 U
1) L
[ ] [ ]
∂ t2
−4 L
∂ x2
+ L [ U ] =L [ 16 x ] +20 L [ senx ]
Haciendo L [ U ( x , t ) ] =u (s ,t )
∂U d2u 16 x 20
[ 2
→ s u−s U ( x , 0 )−
∂t
(x , 0) −4 2 +u=
dx ] s
+ sen x
s
16 x 20
→(s¿¿ 2 u−s ( 16 x )−0)−4 u ' '+u= + senx ¿
s s
16 x 20
s2 u−16 sx−4 u' ' +u= + senx
s s
s2 u '' 1 4x 5
→ −4 sx−u + u= + senx
4 4 s s
'' 1 s2 4x 5
u − u− u=−4 sx− − senx
4 4 s s
1 4 5
u' ' − (s¿¿ 2+1)u=(−4 s− ) x− senx ¿ <M>
4 s s
como u ( 0 , t ) =0 → L [ U (0 , t) ] =L [ 0 ] → u ( 0 , s )=0
16 π
L [ U ( π , t ) ] = L [ 16 π ] → u ( π , s )=
s
Cálculo de uh
1 1
P ( r )=r 2 − ( s 2+ 1 )=0 → r 2= (s 2+1)
4 4
1 2
→ r=±
2
√ s +1
uh ( 0 , s )=c 1 +c 2=0 → c 1=c 2
−1 2 1
2
√ s +1 π √ s2+1 π
uh ( π , s ) =c 1 e +c 2 e 2 =0
−1 2 1
2
√ s +1 π √ s2+1 π
→ c1 e +c 2 e 2 =0
−1 2 1 2
c (e )=0 → c =0
2
√ s +1 π 2
√ s +1 π
1 −e 1
c 2=0
−1 2 1
2
√ s +1 x √ s2+1 x
→ uh=c 1 e +c 2 e 2
Cálculo de u p
u p =Ax +B sen x+ C cos x
u'p =A + B cosx−c senx
u ' 'p =−Bsen x−c cos x
En <M>
1 4 5
(−B senx −c cosx )− ( s 2 +1 ) ( Ax +Bsen x +c cosx )= −4 s− x − senx
4 s s ( )
1 1 A
[ ] [
↔ −B− ( s2 +1 ) B senx+ − c− ( s2 +1 ) c cos x− ( s2 +1 ) x
4 4 4 ]
4 5
(
¿ −4 s−
s )
x− senx
s
1 5 s2 5
↔ *) −B(1+ ( s2 +1 ) )=
4
−5
s
↔B + =
4 4 s ( )
5×4
→ B=
20
s (5+s ¿ ¿ 2)= ¿
s (s 2+ 5)
→ B=4
( 1s − s s+5 )
2
*) −C ¿
1 4
→= ( s 2+1 ) A=−4 s−
4 s
s 2 +1 s2 +1 16
( ) ( ) 4
A=4
s
↔ A=
s
16 x 1 s
∴ u=
s ( s s +5 ) senx
+4 − 2
16 x
L [ u ] =L [
->
−1
s ]
−1
+4 L
[ 1s − s s+5 ] senx−1
2
U ( x , t )=16 x+ 4 ( 1−cos √ 5 t ) senx
¿ 16 x+ 4 senx−4 senx cos ( √ 5 t) ∎
Curso: Ecuaciones diferenciales parciales
Tema: La transformada de Laplace y las ecuaciones en derivadas parciales
Practica Dirigida Nro 02
Problema: Resolver las siguientes ecuaciones en derivadas parciales:
U ( 0 ,t )=0
∂U ∂2 U
1.- =2 ; U (5 ,t )=0
∂t ∂x
2
U ( x , 0 )=10 sen 4 π x
U ( 0 ,t )=0
∂U ∂2 U
2.- =2 ; U ( 5 ,t )=0
∂t ∂x
2
U ( x , 0 )=10 sen ( 4 π x )−5 sen(6 πx )
∂U ∂U U ( x ,0 )=x
3.- − =1−e−t ;
∂ x ∂t U ( x , t ) acotada ; 0< x <1 ,t >0
U ( 0 , t )=0
∂U ∂2 U
4.- =3 ; π ; U ( x , 0 ) =30 cos 5 x
∂t ∂ x2 U
2 ( )
,t =0
∂U ( x , t) ∂ U ( x , t)
5.- +4 =0
∂x ∂t
U ( o , t )=2 ; U ( x , 0 )=0
∂U ∂U
6.- =2 + U ; U ( x , 0 )=6 e−3 x
∂x ∂t
U ( x , 0 )=0
∂U ∂2 U
7.- =9 ; ∂U
∂t ∂ x 2 U ( 3 , t )=2 ; ∂ t ( 0 ,t )=0
2 U ( 0 , t )=0
∂U ∂ U
8.- = 2 ; U ( x , 0 )=3 sen ( 2 πx ) ; U ( 1 , t ) =0
∂t ∂x
0< x <1
2
∂2 U 2 ∂ U ( x,t)
9.- ( x , t ) =a −g
∂t 2 ∂ x2
∂U ( x , 0) ∂U (x , t)
U ( x , 0) = =0 ; U ( 0 , t ) =0 ; lim ¿ x → ∞ =0
∂t ∂x
∂2 U (x , t) ∂2 U ( x , t)
10.- =9 ;
∂ t2 ∂ x2
U ( 0 , t ) =0 ,U ( 2 , t )=0 ; U ( x , 0 )=20 sen ( 2 πx )−10 sen ( 5 πx )
∂U ∂2 U
11.- = 2 −4 U , U ( 0 ,t )=0 , U ( π ,t )=0
∂t ∂x
U ( x , 0 ) =6 senx −4 sen 2 x
∂2 U ∂2 U ∂ U (
12.- =16 , 0 , t )=0 , U (3 , t )=0 , U ( x , 0 ) =0
∂t 2 ∂ x2 ∂ x
∂U (
x ,0 )=12 cos πx +16 cos 3 πx−4 cos 5 πx
∂t
∂2 U ∂2 U
13.- = , x >0 , t> 0
∂t 2 ∂ x 2
∂U (
U ( 0 , t ) =10 sen 2 t , U ( x , 0 )=0 , x , 0 )=0 , y
∂t
lim U ( x , t )=0
x→∞