UNIVERSIDAD DEL MAGDALENA

FACULTAD DE ING. ELECTRONICA

ECUACIONES DIFERENCIALES PARCIALES

TALLER 2°SEGUIMIENTO

DOCENTE: deud soto palomino

PRESENTADO POR:
MARÍA BEATRIZ PINZÓN MEJÍA (2018119088)

GRUPO: 3
2020
R/ En este caso por ser una cuerda infinita, con posicion inicial y velocidad inicial dadas,
podemos usar la formula de d´Alambert;
1
u(x ,t )= ¿
2

con f ( x )=e−|x|, g ( x ) cos 4 x , c=2

x +2t
1
u(x ,t )= ¿+e−|x+2 t|¿ + 1 ∫ cos ( 4 α ) d α
2 2c x−2 t

1 1
u(x ,t )= ¿+e−|x+2 t|¿ + sen 4 α| x+2 t
2 16 x −2t

La solución buscada es:

1
u(x ,t )= ¿ +
2

1
e−|x+2 t|¿ + ( sen ( 4 x +8 t )−sen ( 4 x−8 t ) )
16

Solucion:

Para este caso, cuerda vibrante de longitud L, con extremos fijos, posicion inicial dada y

velocidad inicial 0, tenemos la solución utilizando la formula obtenida por variables

separables:
 ∞
y ( x ,t )=∑ An sen
n=1
( nπxl ) cos ( nπxl )
con An

L
2 nπ
An = ∫ f ( x ) sen
L0 L ( )
x dx

Para este caso C=3 , f ( x )=1−cos ( x ) L=2 π

2π
2 nπ
An = ∫ (1−cosx ) sen
2π 0 2π
x dx ( )
8 (cos ( πn )−1)
An = ó An =8 ¿ ¿ ¿
n ( n2−4 )

Por lo tanto :

∞
y ( x ,t )=∑ An sen
n=1
( n2 x) cos( 32 nt )
∞
y ( x ,t )=∑ 8 ¿ ¿ ¿ ¿
n=1

16 2 n−1 3(2 n−1)

{ 3
¿ ( 2 n−1 ) −4
sen
2
0
(
x cos
2 ) (
t ) sin es impar

sin es par

∞
16 2n−1
y ( x ,t )=∑
n=1
3
( 2 n−1 ) −4
sen
2( )
x cos ( 6 n−3 ¿¿¿ 2 t ) Es solución al problema.
Solucion :

transformada de fourier a ambos lados

∂2 y ∂2 y
F { }
∂t
=144 F
∂x { }
∂2 ^ 2
2
y =−144 w ^y
∂t

∂2 ^ 2
2
y +144 w ^y =0
∂t

Es una E . D . P homogenea lineal de segundo ordencon coeficientes constantes ,

y la ecuacion caracteristica es :

m2 +144 w 2=0

m=± 12 wi

^y ( w, t )=c 1 sen ( 12 wt ) +c 2 cos ( 12 wt )

Transformando las condiciones iniciales

F { ∂∂ yf ( x , 0) }=F {0 }
∂y
( w , 0 ) =0
∂t

F { y (x , 0) }=F { e−5|x|}

^y ∞
iwx
( w, 0) =¿ ∫ e−5|x| e dx ¿
−∞

Mediante el uso de software

^y 10
( w, 0 ) =¿ ¿
w2 +25

10
Luego usando ^y (w , 0)= 2
w + 25

10
^y ( w, 0)=c 1 sen ( 12w 0 ) +c 2 cos (12 w 0 )= 2
w +25

10
c 2= 2
w +25

Luego

10
^y ( w, t )=c 1 sen ( 12 wt ) + 2
cos ( 12 wt )
w +25

∂y
Usando ( w , 0 )=0
∂t

∂y 120 w
( w , 0 ) =c 1 12 w cos ( 12 w 0 )− sen ( 12 w 0 ) =0
∂t w+25

12 w c 1=0=¿ c 1=0

10
^y ( w, t )= 2
cos ( 12 wt )
w +25
Volvendomos a las variables originales

10
F−1 { ^y (w ,t ) }=F−1
{ 2
w +25
cos ⁡(12 wt )
}
∞
1 10
y ( x ,t )= ∫ 2
cos ( 12 wt ) e−iwx dw Es solucion del problema
2 π −∞ w +25

Usamos la transformada de fourier en coseno pues conocemos y (0 , t ) y es una cuerda semi infinita

∂2 y ∂2 y
Fs { } { }
∂ t2
=9 F s
∂ x2

∂2 ^ ( 2
y w ,t )=−9 w ^y s ( w , t ) + wy ( 0 , t ) para y ( 0 , t )=0
2 s
∂t

∂2 ^ ( 2
y w ,t ) +9 w ^y s ( w ,t )=0
2 s
∂t

Es una E . D . P , y su ecuacioncaracteristica es ;

m 2 +9 w 2=0

m=± 3 wi
Luego ^y s=c 1 sen ( wt ) +c 2 cos ( 3 wt )

Ahorausamos las condiciones iniciales para hallar las constantes pero antes≤¿
aplicamos la transformada

∂2 y ∂
Fs { ∂t 2 }
( x , 0 ) =F s { 0 }=¿ ^y s ( w , o ) =0
∂t

y F s { y ( x , 0 ) }=F s { f ( x ) }

^y s ( w , o )=F s f ( x )= x (1−x ) , 0 ≤ x ≤ 1
{ { 0 , x >1 }}
1 ∞
^y s ( w , o )=∫ x ( 1−x ) sen ( wx ) dx +∫ 0 sen ( wx ) dx
0 1

−1
^y s ( w , o )= ( wsen ( w ) +2 cos ( w )−2 )
w3

Entonces haciendouso de lascondiciones iniciales

∂
^y =3 w c1 cos ( 3 wt ) −3 w c 2 sen( 3 wt )
∂t s

∂
^y ( w , 0 )=3 w c1 cos ( 0 ) −3 w c 2 sen ( 0 ) =0
∂t s

3 w c 1=0=¿ c 1=0

^y s (w , t)=c2 cos( 3 wt )

Luego
∂ −1
^y s ( w , 0 )=c 2 cos ( 3 w 0 )= 3 ( wsen ( w )+ 2cos ( w )−2 )
∂t w

1
c 2=¿− ( wsen ( w ) +2 cos ( w )−2 )
w3

−cos ( 3 wt )
^y s ( w , t )= ( wsen ( w )+ 2cos ( w )−2 )
w3

Devolviendo a las variables originales

cos (3 wt )
F−1 { ^y (w ,t ) }=F−1 { w3
( 2−2cos ( w )−wsen(w) ) }
Por lo tanto la solucion buscada es :

∞
2 cos ( 3 wt ) sen ( wx )
y ( x , t )= ∫ [ 2−2cos ( w ) −wsen ( w ) ] dw
η0 w3

∂2 y 2
2 ∂ y
=c 0< x< L, t> 0
∂ t2 ∂ x2

y ( x , 0 ) =Asen ( πxL ) , 0< x< L

∂y
( x , 0 ) =0 para 0 < x< L
∂t

y ( 0 , t )= Asen ( πL0 )=0 ,Y ( L , t )= Asen ( π )=0

Aplicandotransformada de laplas a amboslados

∂2 y ∂2 y
ʆ { } { }
∂ t2
=c
2
ʆ
∂ x2

2 ∂2 ( 2 ∂y(
c 2
Y x , s )=s Y ( x , s )−sy ( x , 0 ) x , 0)
∂x ∂t

∂2 s2 −sA πx
∂x 2
Y −
c2
Y=
c
sen
L ( )
Enuna E . D . P lineal de segundo orden, hallando la Yh

2 s2
m− =0
c2

s
m=±
c

s −s
x x
Y h=c 1 e c + c 2 e c

Ahora hallandoYp mediante elmetodo de coeficientes indeterminados

Y p=Bcos ( πxL )+ Lsen( πxL )

−π 2 πx π 2 πx
Y ´ ´=
L 2
Bcos
L L ( )
− 2 csen
L ( )
Remplazando en la E . D. P

−π 2 πx π 2 πx s2 πx s2 πx sA πx
L 2
Bcos
L
−
L 2
csen ( )
L
−
L 2
Bcos
L ( )
−
L 2
csen
L C( )
= sen
L ( ) ( )
πx π 2 B S2 B πx π 2 c S 2 B sA πx
cos ( )(
L L 2
+ 2 + sen
L ) ( )(
L L 2
+ 2 = sen
L C L) ( )
 π 2+ s 2
B ( L2 )
=0 B=0

π 2+ s 2 sA
c( L2
=
c )
sA L2
c=
c ( π 2+ s 2 )

s
A L2 s
−s
x x πx
Y =Y h+ Y p=c 1 e c + c2 e c
+ 2 2
c(π +s ) ( )
sen
L

Para tener una solucion acotada necesitamos que C 1=0 debido a que
s
x
c
e =¿ ∞ cuando S=¿ ∞ por lo tanto

A L2 s
−s
x πx
Y ( x , s )=c2 e c
+ 2 2
c (π + s )
sen
L ( )
Usando la condicion de frontera

Y ( 0 , t )=0=¿ ʆ { y (0 , t) }= ʆ { 0 }=¿Y ( 0 , s )=0

0 A L2 s
Y 0 , s =c 1 e +
( ) sen ( 0 )=0
c ( π 2+ s 2 )

c 1=0

Luego

A L2 s πx
Y ( x , s )= 2 2
c (π + s )
sen
L ( )
Devolviendo a la variable original
 A L2 πx −1 s
ʆ −1 {Y ( x , s) }=
c
sen
L ( ) {
ʆ
( π + s2 )
2 }
A L2 πx
y ( x , t )=
c
sen ( )
L
cos ( πt ) Es solucion

Aplicanoransformada de laplace

∂2 u ∂2 u
ʆ { } { }
∂t
2
=a ʆ
∂ x2

2
2 ∂u ( 2 ∂
s U ( x , t ) −s a ( x , 0 )− x , 0 )=a U (x ,s)
∂t ∂ x2

2 ∂2 2
a 2
U −s U =0
∂x

∂2 s2
U − 2 U =0
∂ x2 a

Es una E . D . P de segundo ordenhomogenea y su ecuacion caracteristica es :

2 s2 s
m− 2
=0=¿ m=±
a a
 s −s
x x
U ( x , s )=c 1 e a +c 2 e a

Transformando las condiciones de frontera

ʆ {u ( 0 , t ) }= ʆ { f ( x ) }=¿ U ( 0 , s )=F ( s )

ʆ lim u ( x ,t ) = ʆ { 0 }=¿ lim U ( x , s )=0

{ x→ ∞ } x →∞

Usando lascondiciones de frontera

s −s
x x
lim U ( x , s ) =lim c 1 e a +c 2 e a
=0
x→ ∞ x →∞

s
x
Concluimos que C 1=0 porque si x=∞ , e a → ∞

−s
x
a
U ( x , s )=c 2 e

Luego

U ( x , 0 ) =c 2 e0=F ( s )

c 2=F ( s )

−s
x
a
U ( x , s )=F ( s ) e

Devolviendonos a las variables originales

−x
ʆ −1
{ U ( x , s ) }= ʆ −1
{ ʆ {f ( x ) }e } a
s

U ( x , t )=f t− ( ax ) H (t− ax )
o bien
 ( xa ) para t ≥ ax
u ( x , t )=
{
f t−

0 para t<
x
a
}