0% encontró este documento útil (0 votos)
154 vistas13 páginas

EDP Semana 02

1) La transformada de Laplace se utiliza para resolver ecuaciones en derivadas parciales. Al aplicar la transformada a una variable, esta desaparece y la función resulta solo depender de las otras variables y del parámetro s. 2) Al derivar bajo el signo de transformada de Laplace, surgen relaciones entre la transformada de la derivada y la derivada de la transformada. Esto permite resolver ecuaciones diferenciales en derivadas parciales. 3) Se presentan fórmulas clave para transformadas de derivadas primeras y segundas que constituy
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 DOCX, PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
154 vistas13 páginas

EDP Semana 02

1) La transformada de Laplace se utiliza para resolver ecuaciones en derivadas parciales. Al aplicar la transformada a una variable, esta desaparece y la función resulta solo depender de las otras variables y del parámetro s. 2) Al derivar bajo el signo de transformada de Laplace, surgen relaciones entre la transformada de la derivada y la derivada de la transformada. Esto permite resolver ecuaciones diferenciales en derivadas parciales. 3) Se presentan fórmulas clave para transformadas de derivadas primeras y segundas que constituy
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 DOCX, PDF, TXT o lee en línea desde Scribd

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→∞

También podría gustarte