EJERCICIOS DE E.D.

P LINEAL DE VALERIA IRIO (PAG-45)

Seccin 1:

1. Resuelva :

II) = ( ) , > 0 ,
u(x, 0) = x , >0

Solucin:

= ( ) (, ) = ( ) + ()

(, 0) = 0 + () = () = 2

La solucin general es: (, ) = ( ) + 2 , > 0

2. Verifique si las siguientes problemas tienen solucin y en este caso ,halle las
soluciones

II) = , (, ) 2
( 2 , ) = 2
Solucin:
Parametrizando (utilizaremos el teorema 2.1)
: t ((t), (t)) . Definido () = ( 2 , )
((t), (t)) = 0 , ((t), (t)) = 1. Entonces
0 1
| | = 2 0 0 . Esto implica que existe infinitas soluciones para
2 1
0
Resolvemos el sistema para hallar una solucin:
(, ) = 0 , (, 0) = 2 (, ) = = (, ) = 2
{ { {
(, ) = 1 , (, 0) = (, ) = + = (, ) = +
Cambio de variable:
= (, ) =
{ , 0 Tomemos (, ) = (, )
= (, ) =
Entonces (, ) = 2 + (, ) = 2 + + ()
Volviendo a las variables originales
(, ) = 2 + () .Evaluando en las condiciones iniciales tenemos: () =
0 , 0
La solucin general ser:
(, ) = 2 , (, ) = {(, )/ 0 , }
Observacin:
 Si restringimos el abierto a = {(, )/ 0 , + }
= {(, )/ 0 , } la solucin de la E.D.P ser nica (estos anlisis se
hacen graficando)

IV) = 2 , (, ) 2
(, 2 ) = 1
Solucin:
Parametrizando (utilizaremos el teorema 2.1 para garantizar la existencia y
unicidad)
: t ((t), (t)) . Definido () = (, 2 )
((()), (())) = (1,0) . Entonces
1 0
| | = 2 0 0 . Esto implica que existe infinitas soluciones para
1 2
0
Resolvamos el sistema para hallar la solucin:
(, ) = 0 , (, 0) = (, ) = + = (, ) = +
{ { {
(, ) = 1 , (, 0) = 2 (, ) = = (, ) = 2
Cambio de variable:
= (, ) =
{ , 0 Tomemos (, ) = (, )
= (, ) =
Entonces (, ) = 2( + ) 2 (, ) = 2 2 + 2 3 + ()
Volviendo a las variables originales
(, ) = ( )2 + 2( )()3 + () .Evaluando en las condiciones
iniciales tenemos: () = 0 , 0
La solucin general ser:
(, ) = 2 + 1 2 , (, ) = {(, )/ 0 , }
Observacin:
Si restringimos el abierto a = {(, )/ 0 , + }
= {(, )/ 0 , } la solucin de la E.D.P lineal ser nica

Seccin 2:
1. Halle las curvas caractersticas planas
II) 5 + 4 = 3 + 1 + 2 3
Solucin:
Paramerizando
() = ((), ()) = 5 () = 5 + 1 . (1)
() = ((), ()) = 4 () = 4 + 2 . (2)
 De (1) y (2) tenemos:

()1 ()2
= 4(() 1 ) = 5(() 2 ) 4() 5() = 41 52
5 4
Entonces las ecuaciones caractersticas son:
4 5 =

2. Resuelva, indicando la regin del plano donde la solucin es vlida:

() + = 3 + 3, (0, ) = ,

Solucin:

*) Hallamos las curvas caractersticas:

2
= = =
2

Cambiando de variable:

= =
2
{ { 2
= =+
2 2

Tomando (, ) = (, )

(, ) (, ) (, )
= +

5 5
= + = 3 + 3 = ( 3 + 3)
2 2

1 3 1 3 2
(, ) = 4 + 2 + () (, ) = 4 + 2 + ( )
8 2 8 2 2

2
Sea = = 0 =
2
 2 2
(0, ) = () = : ( ) = 2
2

Por tanto:

2
1 3
(, ) = 8 4 + 2 2 + 2 , (, ) 2

() 2 + 2 = , (1, ) = , > 0. ( )

Solucin:

Hallamos las curvas caractersticas:

2
= 2 2 = 2 =

Cambiando de variable:

= =
{ = { =
1

Tomando (, ) = (, )

(, ) (, ) (, )
= +

2 2
1 1
= + ( ) 2 = 2 + ( ) = ( )
1 1 1

1
= ( )
1

(, ) = ln(1 ) + () (, ) = ln ( ) + ( )

 1 1
(1, ) = ln ( ) + ( )
1
1 1
Sea = = 1+

+
() = + ln(1 ) ( )=+ ln( )

Por tanto:

+
(, ) = ln ( ) + , , 0, >
2 1

3.- Resolver los problemas a continuacin, indicando la regin donde la solucin

es vlida:

() 2 + 2 = 2 , (, 2) = 1, > 0.

Solucin:

La curva () = (, 2, 1) es regular, pues los vectores ( (), ()) = (1,2) y

(((), (), ()), ((), (), ())) = ( 2 , 4 2 ) nunca son paralelos, > 0

La solucin es obtenida integrndose a lo largo de las caractersticas que

intersectan la curva . Esto equivale a resolver el sistema:

= = 2 , (0, ) = () =
{ = = 2 , (0, ) = () = 2
= = 2 , (0, ) = () = 1
1 1 1
1) 2 = = +, pero: (0, ) = = = = 1

1 1 1 2
2) = = , pero: (0, ) = 2 = = =
2 + 2 12

1 1 1
3) 2 = = + , : (0, ) = 1 = = 1 = 1

Luego de 1) tenemos: = 1+, y de 2) se tiene: = 2(+1)

Entonces: 1+ = 2(+1)
+ 2
=

Por lo tanto:
2
(, ) = ,
+ 2 1

SECCIN 3:

1. Halle la solucin general de cada uno de los siguientes EDPs:

x
ii). U x sen (Ecuacin Diferencial Parcial lineal de primer grado no homognea)
y

u x
sen
x y

x
y dx
a u sen

1 x
u y sen dx
y y

x
u y cos y
y

iv). u x 3u y senx cos y (Ecuacin Diferencial Parcial lineal de primer grado no homognea)

a. Calcular la solucin general Un de la ecuacin homognea

ux -3uy = 0
1 a=1 b = -3 c=0
2 Resolver
 dy b 3
3
dx a 1

dy 3dx
y 3x

y y
3nn
x x
3 Tomar = x
4 La solucin de
c
d
u , n n e a , como c = 0

u , n n e
0

u , n f n

u , n f n

y
5 u x, y
x
b. Calculo de la solucin particular up por inspeccin: ux = cos x
seny
uy
3
1
ya que cos x 3 seny cos x seny Por lo tanto la solucin particular
3
1
deber ser: up senx cos y
3
La solucin general de la EDP no homo es:
y
u x, y senx cos y
1
x 3
vi) yux xuy = 0 (Ecuacin Diferencial Parcial lineal de primer grado no homognea)

1 a = y, b = -x, c=0
2 Resolver
dy b x

dx a y

dy x

dx y

ydy xdx
y2 x2 K

y2 x2 K

n x, y K y 2 x 2
3 Tomar = x
4 La solucin es
c
a d , como c = 0
u , n n e
u , n f n
5
ux, y x 2 y 2
(Ecuacin Diferencial Parcial casi lineal de primer orden no
viii) ux + xuy = x3+3xy homogneo)

A. Calculando la solucin general un de la ecuacin homognea.

1 a = 1, b = 0, c=x
2 Resolver
 dy 0

dx 1

dy fodx
ykn y

3 Tomar = x
4 La solucin es
c
a d
u , n n e
c x
d 1 dx Xdx x
2
f
a

u , n f n e x
2

u , n f y e x
2

B. Hallando up

u x2 2 3y

ux 2x

u x xu 2 x x( x 2 2 3 y ) 2 x x 3 2 x 3xy x 3 3xy

u p x, y x 2 2 3 y
La solucin general es ux, y x 2 2 3 y y e x
2

5u x 4u y u x 3 1 2e.3 y
A. Calculando la solucin general uM de la ecuacin homognea:
1. a = 5, b = 4; c=1
2. Resolver
 dy 4

dx 5

4
d y dx
5

4 4
y x K y x
5 5

3. x
c
a d
4. La solucin es u , n n e
c 1 1 1
d d 5d 5d x
5
x

1
x
u , n n e 5

4 x
1
u , n y 5 x e 5

B. Hallando la solucin particular.

Por inspeccin: u x 3x 2 30 x 150

30 3 y
uy e
65

u x, y x 3 15 x 2 150 x 749
10 3 y
e
65

120 3 y
5u x 4u y u 15 x 2 150 750 e
65

10 3 y
x 3 15 x 2 150 x 749 e
65

x 3 1 2e 3 y
 u p x, y x 3 15 x 2 150 x 749
10 3 y
e
65

La solucin general es:

4 x
1
u x, y
10 3 y
y 5 x e 3 x 15 x 150 749 65 e
3 2

xii) (x+)x + (y+)y +=0 (, , )

Solucin:

La ecuacin general tiene la forma:

(, ) + (, ) + (, ) + (, ) = 0
= x + , = y + , = , = 0

Resolver la ecuacin:
+
= = +
+ +

Integrando se obtiene:

ln (y+) = ln(x+)

ln (y+) - ln(x+) = lnk

+
ln ( + ) = lnk

+
=k
+

+
tomemos t = ; como ty o, podemos escoger s = x
+
obteniendo entonces el cambio variable.
+ =
s = x; t = + { = ( + ) = ( + )

y + = t (s+)

x+=s+

Tomando (x, y)= V(s, t) satisface la EDP.

 (, ) (, ) (, )
= . + .

Vs = x . Xs + y.ys

VS = x (1) + y (t)

(S+) VS = (S + ) x + (S + )t y

(S+) Vs = v
Vs
= +
V

Integrando con respecto a "" , obtenemos

ln v = - ln (s+)

V= (s +)-

Pasando (, ) (, ); ya que (x, y) = V(s, t)

(x, y) = (x+)-

Ejercicio 3. ii) x + y =

Solucin:

Veamos si ! solucin (por Teorema):

Hallemos las curvas caractersticas de la EDP. Para ello definamos la curva
suave de la E.D.P. tal que admite la parametrizacin.

C : I CR R3

t C (t) = ((t), (t), (t))

tal que:

(t) = a ((t), (t), (t)) = 1

(t) = b ((t), (t), (t)) = 1

(t) = c ((t), (t), (t)) =

Tomemos la curva inicial parametrizada (auxiliar)

 : I R R3

t (t) = (t, -t + 1; (t))

(t) = (1,-1, (t)) 0

es regular

Probar que los vectores (a ((t)), b ((t)) y ( (t) , (t)) nunca son paralelos.

1 1
| |=| | = -2 0
1 1

Los vectores nunca son paralelos

Los vectores no son paralelos y como es regular, entonces existe una nica
solucin para el problema.

Resolver los sistemas caractersticos con condiciones iniciales:

Xs (s, t) = a (x (s, t), y(s, t), V(s, t)) = 1 , X (0, t) = t

Ys (s, t) = b (x (s, t), y(s, t), V(s, t)) = 1 , Y (0, t) = -t + 1

Xs (s, t) = c (x (s, t), y(s, t), V(s, t)) = V , V (0, t) = (t)

De esto tenemos que:

Xs = 1

Integrando con respecto a "" , obtenemos:

X(s, t) = s+c

Pero X (0, t)

X (s, t) = 0+c = t c = t

X(s, t) = s + t

ys = 1

Integrando con respecto a "" , obtenemos

Y (s, t) = S + C2

Pero y (0, t) = -t +1 = C2

Y (s, t) = S + 1 t
 Tomando: (, ) = (, )
(,) (,) (,)
= . + . (por la regla de la cadena)

Vs = x. xs + y. ys

Vs = x (1) + y. (1)

Vs =V
Vs
=1
V

Integrando con respecto a "" , obtenemos

lnV(s, t) = s + c + f (t)

V (s, t) = ef(t) .es+c = k .es.ef(t)

V (s, t) = k es+f (t)

Pero x = s + t y = s+1-t

xt=s y + t -1 = s

x-t=y+t1
xy+1
=t
2

x+y1
Como x = s + t s = 2

Pasando (, ) (, ); ya que (x, y) = V(s, t)

V (s, t) = k es+f (t)

+1 x-y+1
(x, y) = k 2 . f(
2
)
EJERCICIOS DE EDP NO LINEAL DE CAMURRI

(VAMOS A RESOLVERLO UTILIZANDO EL SISTEMA DE CHARPIT)

2. (, ) = 1
Solucin:
Ecuaciones de charpit

= = = =
+ 0 0
De la ltima igualdad se obtiene: = 0 = 0 = =
Tomando: =
1 = 0 1
Resulta el sistema: { = () =
=0

+ = +

La solucin completa de (, ) es: 2 + = +

4) (, ) = 2 2 + 2 1

Solucin:
Ecuaciones de charpit

2
= = 2 2 2
= 3 = 2
2 2 2 + 2 2 2

De la ltima igualdad se obtiene:

= (, , , ) =

2 2 + 2 1 = 0 1
Resultando el sistema {
= 0 2

Reemplazando 2 1 se obtiene: 2 2 2 + 2 1 = 0 2 =
1
2 2 +1
1 1
1 2 1 2
=( 2 2 ) = ( )
+1 2 2 + 1

Reemplazando en la ecuacin: + =
1 1
1
1 2 1 2
Se obtiene (2 2 +1) + (2 2 +1) = + = (2 2 + 1)2

1
Integrando se obtiene: + = (2 2 + 1)2

Resolviendo la integral:

2 2 + 1
() = 2 () = ; sec() =
1 1
1 3 ()
(2 2 + 1)2 =

Utilizaremos la siguiente frmula:

1 2
() = 2 ()() + 2 ()
1 1

Se obtiene:
 3 () 1 1 1
= ( sec() () + () )
2 2
1 1 1
(2 2 + 1)2 = 2 2 + 1 + | + 2 2 + 1|
2 2

La solucin completa de (, ) es:

2( + + ) = 2 2 + 1 + | + 2 2 + 1|

6) = +
Solucin:

=

Integrando en ambos miembros, obtenemos:

ln() ln() = ln() ln ( ) = ln() = = ,
(, , ) = = 0
{
(, , , ) = = 0

Multiplicando a la 1 ecuacin por , y a la 2 ecuacin por ( 1); y luego

restando, se tiene:

1+
+ ( 1) = 0 + 1 = 0 = 1 + =

1+
=

Ahora reemplazando los valores de y en la ecuacin: + =

1+ 1+
+ =

1+
(1 + ) + ( ) =

Integrando, tenemos:
1+ 2
(1 + ) + ( ) = 0
2

Ejercicio 8:
 x 1 p y 2q2

Solucin:

1. Resolver el sistema de Charpit

F x 1 p y 2 q 2 , Las ecuaciones de Charpit son:

dx dy du dp dq

Fp Fq pFp qFq Fx pFu Fy qFu

Calculando:

F p x 1

Fq 2 y 2 q

pFp qFq px 1 q 2 y 2 q 2
px 1 2 y 2 q 2

y 2q2 2 y 2q2

pFp qFq y 2 q 2

Fx p

Fu 0

Fx pFu p p0 p

Fy 2yq 2

Fy qFu 2yq 2

dx dy du dp dq

x 1 2 y q y q
2 2 2
p 2 yq 2

dx dy d dp dq
2u 2
x 1 2y q
2
y q p 2 yq 2

I II III IV V
Igualando II y V

dy dq
2

2y q 2 yq 2

2 yd y qd q
2

2y q2

dy dq

y q

ln y ln q

ln y 1 ln q nc

ln y 1 ln cq

c
y 1 cq q
q

gualando I y IV

dx dp

x 1 p

ln x 1 ln p

1
ln ln p
x 1

a2
p
x 1

Luego:
 pd x qd y d u

a 2 d x ad y
du
x 1 y

a 2 ln x 1 a ln y u

a 2 ln x 1 a ln y u

PROBLEMAS DE ECUACIONES DE LA FSICA MATEMTICA SMIRNOV

Mtodos de las caractersticas

Hallar la solucin general de la siguiente ecuacin

EJERCICIO 21
2 2 + 2 + + = 0
Solucin:

1) Veamos que tipo de EDP es:

= 2

(, ) = 2 ; (, ) = ; (, ) = 2 ; (, ) = ; (, )
= ; (, ) = 0 ; (, ) = 0

= 2 2 2 2 = 0

2) Hallemos las curvas caractersticas

= = 2 = = + = =

3) Hallemos la solucin general

Hallamos su forma cannica , tomando (, ) = (, )
+ = 0
Escoger , de tal manera que satisfaga uno de ellos : {
+ = 0
Tomando (, ) = el cual cumple con la primera ecuacin
 Como la EDP es del tipo parablico su solucin en la forma cannica ser:
(, ) + (, , , , ) = 0
Tomando (, ) =
1
(, ) = 2 ( ) = 1

1 1
(, , , , ) = (2 + + ) + ( 2 ( ) + ( )) = 0
2
Esto implica que:
= 0 (, ) = 1 () (, ) = 1 () + 2 ()
Volviendo a las variables originales
(, ) = 1 () + 2 () (Solucin general)

Ejercicio 25:

Indicar el tipo de E.D.P., luego reducir a la forma cannica y encuentre solucin:

2u 2u 2 u
2
2u 2 u
2
2u
x 2
2 xy y 2 yz z 2 xz 0
x 2 yx y 2 yz z 2 zx

1 caso: i) Para las variables x e y

Tenemos los coeficientes

a1 x, y, z x 2
b1 x, y, z xy
c1 xx, y, z y 2

El discriminante es:

b12 x, y, z a1 x, y, z c1 x, y, z

x2 y2 x2 y2 0

La ecuacin es de tipo parablico

ii) Hallando las caractersticas:

dy b1 x, y, z x y y
2 m1 x, y, z
dx a1 x, y, z x x

dy y

dx x

dy dx

x y

ln y ln k1 x

y k1 x

y
k 1
x

Tenemos x, y, z
y
z

2 caso: i) Para las variables x e z

x 2 u xx 2 zxu xz z 2 u zz

Identificando coeficientes

a 2 x, y , z x 2
b2 x, y, z zx
c 2 xx, y, z z 2

El discriminante es:

b 2 a2 c2 x 2 z 2 x 2 z 2 0

La ecuacin es de tipo parablico

ii) Hallando las caractersticas:

 dy b2 z x z
2 m2 x, y, z
dx a2 x x

dz z

dx x

dz d
z
x
x

ln z ln k 2 x

z k2 x

z
k 2
x

Tomando x, y, z
z
x

u x, y, z v , ,

u x v . x v . x v . x

u xx v x v x v x v x v v
' ' ' ' ' '
x x v xx v . xv . x v

u xx v x v . x v . x x v xx v . xv . x v . x x
u xx v x v x x v x x
2

v xx v . x x v . x x v x v . xx
2

u xx v x v x v x 2v x x 2v x x 2v x x
2 2 2

2v . xx v . xx v . xx
 y z y z y z
2 2

u xx v 2 v 2 v 0 2v 2 2 2v 2 0 2v 2 0
x x x x x x

2y 2z
v . 3 v . 3 v .0
x x

y
2 2
z 2 yz 2y 2z
u xx 2 v 2 v 4 v 3 v 3 v
x x x x x

y2 z2 2 yz 2y 2z
x 2 u xx v
2
v 2 v
2
v v
x x x x x

u yy v y v y v y 2v y y 2v y y 2v y y
2 2 2

v . yy v . yy v . yy

2
1 1 1
u yy v v 0 v 1 2v 0 2v 1 2v 0 1
2 2

x x x

v .0 v .0 v .0

1 2
u yy v 2 v v
x x

y2 2y2
y 2 u yy v yv v
x2 x

u zz v z v z v z 2v z z 2v z z 2v z z
2 2 2

v . zz v . zz v . zz

2
1 1 1
u zz v 0 v v 1 2v 0 2v 01 2v 1
2 2

x
x x

v .0 v .0 v .0
 1 1
u zz v 2 v 2v
x x

z2 2z 2
z 2u zz v z 2
v v
x2 x

u x v x v x v x

u xy v x v x v x v x v x v x
' ' ' ' ' '

u xy v y v . y v . y x v xy v . yv . y v . y x

v xy v . yv . y v . y x v . xy

u xy v x y v x y v x y v xy v y x v x y v x y

v xy v . y x v . y x v x y v . xy

u xy v x y v x y v x y v x y y x v x y y x

v . xy v . xy v . xy v x y y x

u y v . y v . y v . y

u yz v y v y v y v y v y v y
' ' ' ' ' '

u yz v z v . z v . z y v yz v . zv . z v . z y

v yz v . zv . z v . z y v . yz

u yz v y z v y z v z y v yz v z y v y z v y z

v yz v . z y v . z y v y z v . yz
u yz v y z v y z v y z v y z z y v y z z y

v y z z y v . yzv . yz v yz

y 1 z y 1 z
u xy v 2 v 2 0 v 0 1v 2 0 2
x x x x x x

y 1 z 1
v 2 1 0 v 2 1 00 v 2 v 0 v 0
x x x x

y z y z 1
u xy v 3 3 v 2 v v 2 2 v
x x x x x

2y2 2 yz 2y2 2 yz 2y
2 xyuxy 2
v 2
v v v v
x x x x x

u yz v z y v y z v y z v y z z y v y z z y

v y z z y v . yz v . yz v . yz

1 1 1 1
u yz v 0 v 0 v 11 v 00
x x x x

1 1
v 1 0 1 v 01 1 v 0 v 0 v 0
x x

1 1 1
u yz v v v v
2
x x x

2 yz 2 yz 2 yz
2 xyu yz 2 yzv 2
v v v
x x x

u zx v z x v z x v z x v z x x z v z x x z

v z x x z v . zx v . zx v . zx
 1 z z y 1
u zx v 10 v 2 v 10 v 0 2 2
x x x x x

y 1 z 1
v 00 2 1 v x 0 x 2 1 v 0 v x 2 v 0
x

z y y z 1
u zx v 3 3 v 2 v 2 v 2 v
x x x x x

2z 2 2 yz 2 yz 2z 2 2z
2 zxuzx v v v v v
x2 x2 x x x

Sumando:

y2 z2 2 yz 2y 2z
v
2
v 2 v
2
v v
x x x x x

y2 2y2
v y 2
v v
x2 x

z2 2z 2
v z 2
v v
x2 x

2y2 2 yz 2y2 2 yz 2y
2
v 2
v v v v
x x x x x

2 yz 2 yz 2 yz
2 yzv 2
v v v
x x x
2z 2 2 yz 2 yz 2z 2 2z
2
v 2
v v v v
x x x x x

2y 2z 2z 2 2 yz
v v y 2 v z 2 v v v
x z x x

2y 2 yz 2z 2 2z
v 2 yzv v v v 0
x x x x

y 2 v z 2 v 2 yzv 0

y z 2 v 0

2 v 0

v 0

d dv
0
d d

dv
d
d
0d

,
dv
d

d v , d

, , ...x

como u x, y, z v , ,

y z
y , , z y
x x
 x es :

y z y z
u x, y, z z y , ,
x x x x

PROBLEMAS DE FSICA MATEMATICA DE TIJANOV

Ejercicio 53
Una cuerda infinita se excita mediante una desviacin inicial localizada con
forma de parbola cuadrtica
a) Hallar las formulas que representan el perfil de la cuerda para > 0

-c c

Como se trata de una cuerda entonces la ecuacin es de onda

representada
= 2 ,
Vamos a determinar la forma de la cuerda en el instante = 0 ( (, 0) =
() )

Hallando ()
() = 1 2 2 + = 0
{ 22 = 0 2 = 0
() = 1 2 + 2 + = 0

Entonces: () = 1 2 + (0) = = () = 1 2 + = 0

1 = 2
2
Esto implica: () = 2 2 + () = (1 2 ) esta funcin
representa la forma que toma la onda en el instante = 0
 Ahora vamos hallar la funcin que plantea la distribucin de las velocidades
( = () ) en el instante = 0 , donde se deduce () = 0
Dndole una forma matemtica a nuestro problema tenemos:
= 2 ,
0 , < <
2
(, 0) = () = (1 2 ) , < <

{ 0 , < < +
{ (, 0) = 0
Su solucin est dada por la formula de DAlambert:

1
(, ) = [( ) + ( + )] 1
2
Donde:
( ):
( + ):

Para obtener las formulas en la condicin del problema estudiaremos la

particin del plano por las caractersticas de 1 trazadas en los extremos del
intervalo (, ) donde el desplazamiento inicial es diferente de cero
() 0 )

Daremos en primer lugar las formulas que determinan el perfil de la cuerda

cuando t=cte. limitndonos a dos casos :

0 < < < < +

Si , t = cte , 0 < < +
Para < <
< < 2
{ ( ) = 0 ( + ) = 0 1 (, )
< + <
=0

Para < < +

2 < <
{ ( ) = 0 ( + ) =
< + < + 2
(+)2
(1 )
2
( + )2
2 (, ) = (1 )
2 2
Para + < <
 < < 2 ()2
{ ( ) = (1 2 ) ( +
+ 2 < + <
(+)2
) = (1 )
2

2 + 2 2
3 (, ) = (1 )
2
Para < < +
2 < < ()2
{ ( ) = (1 2 ) ( + ) = 0
< + < + 2

( )2
4 (, ) = (1 )
2 2
Para + < < +
< < +
{ ( ) = 0 ( + ) = 0
+ 2 < + < +

5 (, ) = 0
Donde resulta:
0 , < <
( + )2
(1 ) , < < +
2 2
2 + 2 2
(, ) = (1 ) , + < <
2
( )2
(1 ) , < < +
2 2
{ 0 , + < < +

b) Hallar las formulas que representan la ley del movimiento de los puntos de
la cuerda con diferentes abscisas para t>0

Ahora daremos las formulas que determinan u(x,t) para x=cte. Y representa la
ley del movimiento del punto de la cuerda con la abscisa fija. Escogeremos un
valor fijo de x en cada uno de los intervalos: < < , < < 0 ,
0< < , < <+ y estudiaremos como varia la expresin de la
solucin donde t varia de 0 hasta + (es anlogo a lo que se hizo anterior)

Para < <

+ + 2 < <
0<< { ( ) = 0 ( + )
< + <
=0
1 (, ) = 0
 + 2 < < 2 +
<< { ( ) = 0 ( + )
< + <
( + )2 ( + )2
= (1 ) 2 (, ) = (1 )
2 2 2

< < 2
< < + { ( ) = 0 ( + ) = 0
< + < +
2 (, ) = 0

Agrupando:
+
0 , 0<<
(+)2
(, ) = (1 ) , << , < <
2 2

{ 0 , < < +

Anlogamente se obtienen para los dems intervalos

2 +2 2 +
(1 ) , 0<<
2
(+)2 +
(, ) = (1 ) , << , < < 0
2 2

{ 0 , < < +

2 +2 2
(1 ) , 0<<
2
()2 +
(, ) = (1 ) , << , 0 < <
2 2
+
{ 0 , < < +

+
0 , 0<<

( )2 + +
(, ) = (1 ) , < < , < < +
2 2
+
{ 0 , < < +

EJERCICIO DE FSICA MATEMTICA N 58

Resolver el problema de la propagacin de las oscilaciones elctricas en un cable

infinito con la condicin
 GL = CR

Donde G, L, C, R son la conductancia de fugas, la autoinduccin, la capacidad y la

resistencia por unidad de longitud del cable*.La tensin y la intensidad de la
corriente en el cable en el instante inicial se suponen conocidas.

*Esta condicin garantiza la posibilidad de que las ondas pasen por el cable sin
variar su forma .En lo sucesivo, a un cable que cumple esta condicin lo
llamaremos: cable de transmisin sin distorsin.

Aplicarem
os:

La ley de voltajes de kirchoff:

Ley de ohm: La diferencia de potencial existente entre dos extremos de un

conductor atravesado por una corriente elctrica es igual al producto de la
intensidad de esta por la resistencia entre ambos puntos.

La diferencia de potencial en la bobina de inductancia es:

VL = -L


=

Veamos:

=

= ( ) + () + + |


0 = + ( ) + () ++


+ (R dx)i + (Ldx) = 0

Esto lo podemos escribir como:

Vx + Ri + Lit = 0

Ahora hallemos la segunda ecuacin:

Aplicaremos:

La diferencia de potencial entre las placas de un condensador = carga

almacenada /capacitancia

= /

Ley de Ohm:

Q = IT

R = 1/G

de Q = IT, tenemos

= + =

()
= = =

= /

de

I=i3 +i2 i3= I-i2

i2 = (Gdx) V i3= I-(Gdx) V ..()

Adems: i3 = i4 + I +

i3 = i4 + I + ..()

Pero: i4 = (Cdx)

En ():

i3 = (Cdx) + I + ..()

De () y () tenemos:

I-(Gdx) V = (Cdx) + I +

0 = + () + ()

Esta ecuacin la podemos escribir como:

Ix + CVt + GV=0 (2)

Diferenciando (1) con respecto a x obtenemos:

+ + = 0

= (3)
Diferenciando (2) con respecto a t obtenemos:

Ixt + GVt + CVtt = 0. (4)

LIxt + LGVt + LCVtt = 0.. (4) X L

Reemplazando (3) en la ecuacin anterior:

+LGVt + LCVtt = 0 (5)

+ CVt + GV = 0.. De (2)

CVt + GV = -. (6)

(6) en (5):

-Vxx + R CVt + GV + LGVt + LCVtt = 0

Ordenando: VXX = LCVtt + RC + LGVt + RGV = 0 (7)

Si G=R=0:

VXX = LCVtt

Por comodidad:

1/k2 =

VXX = 1/k2 Vtt

Entonces tenemos el siguiente problema:

VXX = 1/k2 Vtt

(, 0) = ()
1
Vt(x, t) = g(x) = ( + )
2

Hallemos la solucin de la ecuacin caracterstica

2 2
a(x, t)dt -2b(x,t)dxdt +c(x, t) dx = 0

2
donde: a = 1; b= 0; c = - 1/k

(dt)2 - 1/k2 (dx)2 = 0 ( )2 = k2 = k

x = kt + C1 x = -kt + C2
 x - kt = C1 x + kt = C2

Ahora definamos el cambio de coordenadas:

x- kt = x + kt =

Ahora despejemos x e y:

Definamos V(x, t) = V*(, )

Calculemos las derivadas:

Vxx = V* (x )2+ 2V* x x + V* (x )2 + V* xx + V*xx =

V* (1 )2+ 2V* (1) (1) + V* (1 )2 + V* (0) + V*(0)

Vxx = V* + 2V* + V*

Vtt = V* (t )2+ 2V* t t + V* (t )2 + V* tt + V*tt =

V* (-k )2+ 2V* (-k) (k) + V* (k )2 + V* (0) + V*(0)

Vtt = k2V* - 2k2V* + k2V*

Reemplazamos en la ecuacin:

VXX = 1/k2 Vtt

V* + 2V* + V* = 1/k2 ( k2V* - 2k2V* + k2V* )

V* + 2V* + V* -V* + 2V* -V* =0

4V* =0

En estas nuevas coordenadas la ecuacin de onda es V* = 0

Esta es la forma canonca de la ecuacin de la onda, resolvamos esta ecuacin:

Primero (V*) = 0

De esto tenemos que V* es independiente de, luego: V* = h(n)

Integrando obtenemos: V*(,) = h(n)dn + F ()

F (): es una constante de integracin de la derivada parcial respecto a , h(n)dn :

es solo otra funcin de , la que escribiremos como G(). As:
 V*(,) = G () + F ()

Donde F y G deben ser funciones de una variable con segunda derivada continua,
por lo dems son arbitrarias. Por tanto pasamos a las coordenadas originales
tenemos:

V (,) = F(x-kt) + G(x+kt)

Esta ecuacin se llama la solucin DAlambert de la ecuacin de onda.

Utilizando las condiciones iniciales:

V(x, 0) = G(x) + F(x) = f(x)..... (#)

Vt(x, t) = -kF(x-kt) + kG(x+kt)

Vt(x, 0) = -kF(x) + kG(x) = g(x). (##)

Integrando la ltima ecuacin:

x
0 Vt (,0) d = -k0 ()d + k0 ()d = 0 () d

-kF(x)-F(0) + kG(x) G(0) = 0 ()d

1
-F(x) + G(x) = 0 () d - F (0) + G (0)

1
F(x) -G(x) = 0 () d + F (0) -G (0)

Sumando esta ecuacin con la ecuacin (#)

1
2F(x) = f (x) - 0 () d +F (0) - G (0)

De donde:
1 1 1 1
F(x) = f (x) - 2 0 () d + 2 F (0) - 2 G (0)
2

Tambin de la ecuacin (#)

1 1 1 1
G(x) = f(x)-G(x) = 2 f (x) +2 0 () d - 2 F (0) + 2 G (0)

Ahora reemplazando en la solucin, tenemos que

v(x, t) = F (x-kt) + G(x + kt)

 1 1 1 1
= 2 f (x-kt) - 2 0 () d + 2 F (0) - 2 G (0)

1 1 + 1 1
+ 2 f(x+kt) + 2 0 () d - 2 F (0) + 2 G (0)

1 1 +
V(x, t) = 2 (f(x-kt) + f (x + kt)) + 2 () d

Pero R= G= 0

1
Entonces : g(x) = ( + )=0
2

1 1 +
V(x, t) = (f(x-kt) + f (x + kt)) + 0d
2 2

1
V(x, t) = (f(x-kt) + f (x + kt))
2

Hallemos la solucin de la ecuacin para la tensin usando la ecuacin

VXX = LCVtt + RC + LGVt + RGV = 0 (7)

Sujeta a las condiciones

(, 0) = ()
1
Vt(x, t) = g(x) = ( + )
2

Y usando la condicin:

=

Otra forma de escribir la ecuacin (7) es la siguiente:

1
Vxx LC + ( + ) 2 (, ) + ( )2 V = 0 ( )
2 4

Comprobando:

1
Vxx LC Vtt +( + ) + 4 ( + )2 V + ( )2 V =0
4

+
Vxx LC Vtt LC( ) ( + )2 V+ ( )2 V =0
4 4

Vxx LC Vtt ( + ) = 0

En ( ) reemplazamos la ecuacin : =

1
Vxx LC + ( + ) 2 (, ) = 0
2

1
Vxx LC 2 /t2 + ( + ) /t + 4 ( + )2 V(x, t) = 0

1
Vxx LC Vtt + ( + ) Vt + 4 ( + )2 V = 0

Vxx LCVtt ( + ) Vt - ( + )2 V = 0
4

1 1
- Vxx + Vtt + ( + ) Vt + 4 ( + )2 V = 0 ( )

tomando en ( ) una solucin de la ecuacin de la forma :

1
V(x,t) = 2( + ) .U(x,t)
1 1
1
Vt = ( + ) . 2( + ) .U(x,t) + 2( + ) .Ut(x,t)
2

1 2 1
Vtt = ( + ) . 2( + ) .U(x,t)
4

1
1
2 ( + ) . 2( + ) . Ut (x,t)

1 1
1
- ( + ) . 2( + ) . Ut (x,t) + 2( + ) . Utt (x,t)
2

1
( + )
Vx = 2 . Ux(x, t)

1
( + )
Vxx = 2 . Uxx(x, t)

Reemplazamos en ( ), obtenemos:
1
- Uxx(x, t) + Utt(x, t) = 0

1
Utt(x, t) - Uxx(x, t) = 0 ( )

1
Ya que V(x, t) = 2( + ) .U(x, t)
1
V(x, 0) = 2( + )0.U(x, 0) = 0 .U(x, 0) = U(x, 0)

V(x, 0) = U(x, 0)
 f (x) = U(x,0)

1 1
1
Vt(x, t) = ( + ) . 2( + ) .U(x, t) + 2( + ) .Ut(x, t)
2

1 1
1
Vt(x, 0) = ( + ) . 2( + ).0 .U(x, t) + 2( + ).0.Ut(x, t)
2

1
g(x) = ( + ).U(x, 0) +Ut(x, 0)
2

1 1
2 ( + )f(x) = 2 ( + ) f(x) + Ut(x,0)

Ut(x, 0) = 0

Ahora tenemos un nuevo problema:

1
Utt - Uxx = 0 . de ( )

U(x, 0) = f (x)

Ut(x, 0) = 0

Hallemos la solucin:

Consideremos:

1
k2 =

Utt - k2 Uxx = 0

hallemos la ecuacin caracterstica:

2 2
a(x, t) (dt) -2b(x, t) dxdt +c(x,t) (dx) = 0

donde: a = - k2 , b= 0, c = - 1

- k2 (dt)2 +0+ (dx)2 = 0 (dx)2 = k2 (dt)2

x = kt + C1 x = -kt + C2

x - kt = C1 x + kt = C2

Ahora definamos el cambio de coordenadas:

x- kt = x + kt =

Definamos U (x, t) = V (, )
Calculemos las derivadas:

Uxx = V (x)2+ 2V x x + V (x)2 + V xx + Vxx =

V (1)2+ 2V (1) (1) + V (1)2 + V (0) + V(0)

Uxx = V + 2V + V

Utt = V (t )2+ 2V t t + V (t )2 + V tt + Vtt =

V (-k )2+ 2V (-k) (k) + V (k )2 + V (0) + V(0)

Utt = k2V - 2k2V + k2V

Reemplazamos en la ecuacin :

Utt - k2 Uxx = 0

( k2V - 2k2V + k2V) - k2 (V + 2V + V ) = 0

-4 k2 V = 0

En estas nuevas coordenadas la ecuacin de onda es V = 0

Esta es la forma canonca de la ecuacin de la onda, resolvamos esta ecuacin:

Primero (V) = 0

De esto tenemos que V es independiente de, luego: V = h(n)

Integrando obtenemos: V (,) = h(n) dn + F ()

F (): es una constante de integracin de la derivada parcial respecto a , h(n)dn :

es solo otra funcin de , la que escribiremos como G(). As:

V (,) = G () + F ()

Donde F y G deben ser funciones de una variable con segunda derivada continua,
por lo dems son arbitrarias. Por tanto pasamos a las coordenadas originales
tenemos:

U(x, t) = F(x-kt) + G(x+kt)..()

Esta ecuacin se llama la solucin DAlambert de la ecuacin de onda.

Utilizando las condiciones iniciales:

U(x, 0) = G(x) + F(x) = f(x)..... (1)

Ut(x, t) = -kF(x-kt) + kG(x+kt)

Ut(x, 0) = -kF(x) + kG(x) = 0. (2)

Integrando la ltima ecuacin:

x
0 Ut (,0) d = -k0 ()d + k0 ()d = 0 0d

-kF(x)-F (0) + kG(x) G (0) = 0

-F(x) + G(x) = - F (0) + G (0).()

Sumando esta ecuacin con la ecuacin (1)

U(x, 0) = G(x) + F(x) = f(x)

2G(x) = f (x) -F (0) + G (0)

1 1 1
G(x) = f (x) - 2 F (0) + 2 G (0)
2

1 1 1
F(x) = f (x) + F (0) - G (0)
2 2 2

En () : U(x,t) = F(x-kt) + G(x+kt)

1 1 1 1 1 1
= f (x-kt) + 2 F (0) - 2 G (0)+ f (x+kt) - 2 F (0) + 2 G (0)
2 2

1
U(x,t) = f (x-kt) + f (x+kt)
2

1
( + )
Pero V(x,t) = 2 .U(x,t)

1
1 1 1
(, ) = 2
2( + ) f(x -
t) + f(x +
t)

Adems: =

1 1
1
(

+ ) f( x t ) + f( x + t)
= 2

2
 1 1
f( x t ) + f( x + t)
V(x, t) = ( )

2

Hallemos la solucin de la ecuacin para la intensidad de corriente, usando las

ecuaciones siguientes:

Vx + RI + LIt = 0 (1)

Ix + CVt + GV=0 (2)

Y las condiciones:

I(x, 0) = F(x)

1
It(x, t) = ( + ) F(x)
2

Adems : =

Diferenciando (1) con respecto a t y (2) con respecto a x

Vxt + RIt + LItt = 0 (1)

Ixx + CVtx + GVx=0 (2)

CVxt + RCIt +LCItt= 0 (1C) ()

CVtx = - Ixx GVx de (2)

En (): - Ixx GVx + RCIt +LCItt = 0. ()

de (1): -Vx = RI + LIt ()

() en ():

- Ixx + G (RI + LIt) + RCIt +LCItt= 0

Ixx = LCItt + (LG + CR) It +RGI

Otra forma de escribir la ecuacin anterior es la siguiente:

 2 1
LC + ( + ) 2 (, ) + ( )2 I(x, t) = 0
2 2 4

Usando la condicin: =

2 2 1
LC 2 + ( + ) + 4 ( + )2 I(x, t) = 0
2

2 2 1
LC 2 + ( + ) + 4 ( + )2 I(x, t) = 0
2

IXX LCItt LC ( + )It - ( + )2 I(x, t) = 0
4

1 1
- Ixx + Itt + ( + ) It + 4 ( + )2 I(x,t) = 0 (*)

tomando en ( ) una solucin de la ecuacin de la forma:

1
I(x, t) = 2( + ) .V(x, t)
1 1
1
It = ( + ) . 2( + ) .V(x, t) + 2( + ) .Vt(x, t)
2

1 2 1
Itt = ( + ) . 2( + ) .V(x,t)
4

1
1
2 ( + ) . 2( + ) . Vt (x, t)

1 1
1
2 ( + ) . 2( + ) . Vt (x,t) + 2( + ) . Vtt (x, t)

1
Ix = 2( + ) . Vx(x, t)
1
( + )
Ixx = 2 . Vxx(x, t)

Reemplazamos en la (*) obtenemos:

1
vtt(x, t) - vxx(x, t) = 0. ()

1
Ya que I(x, t) = 2( + ) .V(x, t)
1

I(x, 0) = 2( + )0.V(x, 0) = 0 .V(x, 0) = V(x, 0) = F(x)

V(x, 0) = F(x)

 1 1
1
It(x, t) = ( + ) . 2( + ) .V(x, t) + 2( + ) .Vt(x, t)
2

1 1
1
It(x, 0) = ( + ) . 2( + ).0 .V(x, t) + 2( + ).0.Vt(x, t)
2

1
It(x, 0) = ( + ).V(x, t) + Vt(x, t)
2

1 1
2( + ) F(x) = 2 ( + ).V(x, 0) + Vt(x, 0)

1 1
2( + ) F(x) = 2 ( + ) . F(x) + Vt(x, 0)

Vt(x, 0) = 0

ahora tenemos un nuevo problema :

1
Vtt - Vxx = 0 ..de ()

V(x, 0) = f (x)

Vt(x, 0) = 0

Hallemos la solucin:

1
Consideremos: k2 =

Vtt - k2 Vxx = 0

Hallemos la ecuacin caracterstica:

2 2
a(x,t)(dt) -2b(x,t)dxdt +c(x,t) (dx) = 0

donde: a = - k2 , b= 0, c =1

- k2 (dt)2 +0+ (dx)2 = 0 (dx)2 = k2 (dt)2

x = kt + C1 x = -kt + C2

x - kt = C1 x + kt = C2

Ahora definamos el cambio de coordenadas:

x- kt = x + kt =

Definamos V(x, t) = V*(, )

Calculemos las derivadas:

Vxx = V* (x)2+ 2V* x x + V* (x)2 + V* xx + V*xx =

V* (1)2+ 2V* (1) (1) + V* (1)2 + V* (0) + V*(0)

Vxx = V* + 2V* + V*

Vtt = V* (t )2+ 2V* t t + V* (t )2 + V* tt + V*tt =

V* (-k )2+ 2V* (-k) (k) + V* (k )2 + V* (0) + V*(0)

Vtt = k2V* - 2k2V* + k2V*

Reemplazamos en la ecuacin:

Vtt -k2 Vxx = 0 obtenemos

-4k2V* =0

V* = 0

En estas nuevas coordenadas la ecuacin de onda es V* = 0

Esta es la forma canonca de la ecuacin de la onda, resolvamos esta ecuacin:

Primero (V*) = 0

De esto tenemos que V* es independiente de, luego: V* = h(n)

Integrando obtenemos: V*(,) = h(n) dn + H ()

F (): es una constante de integracin de la derivada parcial respecto a , h(n)dn :

es solo otra funcin de , la que escribiremos como G(). As:

V*(,) = G () + H ()

Donde F y G deben ser funciones de una variable con segunda derivada continua,
por lo dems son arbitrarias. Por tanto pasamos a las coordenadas originales
tenemos:

V(x, t) = H(x-kt) + G(x+kt)

Esta ecuacin se llama la solucin DAlambert de la ecuacin de onda.

Utilizando las condiciones iniciales:

Utilizando las condiciones iniciales:

V(x, 0) = G(x) + H(x) = F(x) ..... (1)

Vt(x, t) = -kH(x-kt) + kG(x+kt)

Vt(x, 0) = -kH(x) + kG(x) = 0=g(x) . (2)

Integrando la ltima ecuacin:

x
0 Vt (,0) d = -k0 ()d + k0 ()d = 0 ()d

-kH(x)-H(0) + kG(x) G(0) = 0 ()

1
-H(x) + G(x) = 0 () - H (0) + G (0)

Sumando esta ecuacin con la ecuacin (1)

G(x) + H(x) = F(x)

1
Obtenemos: 2G(x) = 0 () - H (0) + G (0) + F(x)

1 1 1 1
G(x) = 2 0 () - 2 H (0) + 2 G (0) + 2 F(x)

1 1 1 1
H(x) = F(x) ()+ 2 H (0) - 2 G (0)
2 2 0

Pero V(x, t) = H(x-kt) + G(x+kt)

1 + 1 1 1
V(x, t) =2 0 () d - 2 H (0) + 2 G (0) + 2 F(x+kt) +

1 1 1 1
F(x-kt) () + 2 H (0) - 2 G (0)
2 2 0

1 1 +
V(x,t)= 2 F(x+kt) + F(x-kt) + 2 ()d

1
Pero V(x, t) = 2( + ) .I(x,t)
1
1 1 +
V(x, t) = 2 F(x+kt) + F(x-kt) + 2 ()d = 2( + ) .I(x,t)
Como g(x) = 0

1
1 1 +
F(x+kt) + F(x-kt) + 2 0d = 2( + ) .I(x,t)
2

1
1
F(x+kt) + F(x-kt) = 2( + ) .I(x, t)
2

1
1
I(x, t) = F(x+kt) + F(x-kt) 2( + )
2

1
1 1 1 1
= 2( + ) F(x + t) + F(x - t) . k =
2

Adems: =

1 1
(
1
+ ) f( x+ t ) + f( x t)
= 2

2

1 1
f( x+ t ) + f( x t)
I(x, t) = ( )

2

5. Una cuerda homognea de longitud l. sujeta a los extremos x =0 y x= 1,

vibra bajo la accin de una fuerza exterior armnica F x, t f x semcut

, calculador por unidad de longitud. Hallar la desviacin de ux, t .

Investigar la posibilidad de resonancia y, el dado el caso, hallar la

solucin.
Por los datos del problema el modelo a seguir el de la propagacin de la
onda en un intervalo, con trmino exterior y condiciones homogneas en
la frontera.
 u tt c 2 u xx g x, t ,0 x L, t 0

u x,0 x

u t x,0 x

u 0, t 0, u l , t 0

Primer dato: una cuerda homognea de longitud L, es decir 0 < x < L.

Segundo dato: como f(x, t) es fuerza exterior armnica, es decir:

g x, t F x, t f x semcut

Tercer dato: ux,0 x mide la desviacin de ux, t en cualquier

punto x que esta entre 0 y L, pero cuando t = 0.

Cuarto dato: u t x,0 x mide la velocidad inicial en cualquier punto

de la cuerda.

Quinto dato: u0, t 0, ul , t 0 por que la cuerda esta sujeta a los

extremos.

u tt c 2 u xx F x, t

u x,0 x

u t x ,0 x

Quedando:

u0, t 0, ul , t 0

Solucin:

La solucin de este problema esta dado por:

 ux, t vx, t x, t .......I

Donde vx, t es solucin del siguiente sistema:

u tt c 2 u xx 0.......1

u x,0 x .........*

u t x,0 x

u 0, t 0, u l , t 0

Y x, t es la solucin del sistema siguiente:

tt c 2 xx F x, t

x,0 0.........* *

t x ,0 0

0, t 0, l , t 0

1) solucionando el sistema (*)

x, t T t X x

vtt xt T ' ' t

v xx x ' ' T t

Sustituyendo en (1) y separando las variables tenemos:

T ' ' t x' ' x

x
c T t
2
X x

 T ' ' t c 2T t 0 X ' ' x xx 0...2

1) caso: 0

p2 , p 0

Reemplazando en (2)

T ' ' t c 2T t 0

Es una ecuacin diferencial de segundo

orden su ecuacin caracterstica es:

2 p 2c 2 0

pc

T t c1e pct c2 e pct

Usando la condicin de frontera u0, t ul , t 0

v0, t T t X x 0..........3

T t X 0 0

X 0 0 , pues se busca ux, t 0

vL, t T t X L 0..........4

X L 0 , pues se busca ux, t 0

Resolviendo
 X ' ' x X x 0

X ' ' x p 2 X x 0

Su ecuacin caracterstica es:

2 p2 0

X x c1e pt c2 e pt

v x, t c1e pct c2 e pct b1e px b2 e px

v0, t c1e pct c2 e pct b1 b2

b1 b2 0 b1 b2

vL, t c1e pct c2 e pct b1e pL b2 e pL
b1e pL b2 e pL 0

b2 e pL b2 e pL 0

b2 e pL e pL 0
0

b2 0

Se tendra ux, t 0, se busca una solucin ux, t 0

2) caso: p 2 0
 En X ' ' x p 2 X x 0

X ' ' x 0

r 2 0 r 0m 2

X x a bx Satisface la condicin (3) y (4)

0 x0 a a 0

0 X L a Lb Lb 0 b 0

Luego F x 0 (no es til) porque ux, t 0 y

Se busca ux, t 0

2) caso: p 2 0

X ' ' x p 2 X x 0

Su ecuacin es caracterstica es:

r 2 p 2 0 r pi

X x A cos px Bi sec mx

Hacer A a cos , B a sec

X x a cos px cos asenxsen

a cos px

De manera anloga procedemos para T t

T t b cosct p
 vx, t a cos px b cosct

vx, t d cos px cosct , donde s = ab de las

condiciones de frontera u0, t 0, uL, t 0

d cos ct 0 , para t > 0

d cos L 0 cosct 0, t 0

Se cumple para cos = 0 = /2

Se cumple para cos pL 0

pL 2n 1
2

pL n
2 2

pL n

n
p
L

n n
v n x, t d cos cos ct
L 2 L

n n
d cos x cos sen xsen
L 2 L 2
n n
cos L c cos sen L ctsen

n n n
sen x d cos cos ct dsensen ct
L L L
Donde cn dsen kn d cos

n n n
vn x, t sen x cn sen ct k n cos ct
L L L

Por la superposicin de soluciones:

Solucin general

n n n
vn x, t sen x cn sen ct k n cos ct
n 1 L L L

Con condicin inicial ux,0 x

n
u x,0 k n sen x x
n 1 L

n
x sen
2 L
kn
L 0 L
xdx

Con la condicin ut x,0 x

n n
ut x,0 cn csen x x
n 1 L L

n
x sen
2 L
cn xdx
nc 0 L

Luego hallando la solucin x, t es solucin particular dada por:

la cual se halla haciendo de manera anloga a lo anterior)

n
x, t bn t sen x
n 1 L

bn t Son coeficientes que deben determinarse si se asume que

 b0 0 , bn 0 0 n 1,2,...

Con condiciones de frontera se satisfacen:

b 0senp
n 1
n n x 0, 0 x L

b' 0senp
n 1
n n x 0, ...0 x L

Para hallar bn t calculamos:

t b' n t senp n x

tt b' ' n t senp n x

x bn t n cos p n x

xx bn t n2 senp n x

b' ' n n2 c 2 bn t senpn x F x, t

Ftt t b' ' n t p n2 c 2 bn t F x, t sen n xdx

2 L
l 0

Est es una ecuacin diferencial ordinaria que se puede

resolver por variacin de parmetros:

b' ' n t pn2 c 2 bn t F x sentsen n xdx

2 L

l
0

p F x sentsen n xdx
2 L
(I)
l 0
 sent f x sentsenn xdx
2p
l

n
f t senp ct d
1 t
bn n n pn
pn c 0 L

Quedando:

n n n n
u x, t sen x c n sen ct k n cos ct bn t sen x
n 1 L L L n 1 L

Seccin: mtodo de propagacin de Ondas

(Mtodo de DALAMBERT)

Ejercicio 60 (Ecuaciones de la fsica Matemtica SAMARSKY)

A una barra semi infinita 0 x con el extremo x 0 libre se

le comunica una velocidad axial igual a V0 en el segmento [c, 2c] e
igual a cero fuera de este segmento.

El valor del desplazamiento longitudinal de las secciones

transversales de la barra se puede representar, para mayor claridad,
en una direccin perpendicular al eje x, es decir se puede proceder
como en el caso de la cuerda. Hallar la ecuacin del desplazamiento
de la barra semi infinita.

Solucin
 Ecuacin dif. De segundo grado de la propagacin
unidireccional de una cuerda.

u tt a 2 u xx , 0 x, t

El desplazamiento de la cuerda en el instante t = 0 (parte de

reposo) como no nos dicen nada de que es diferente

ux,0 0, 0 x

Con extremo x = 0 libre, significa que la tensin de la cuerda

en x = 0 es cero, es decir no hay fuerzas extremas.

u
T x, t T 0, t K 0
x x 0

u x 0, t 0 0t

La velocidad axial inicial igual a V0 en el segmento [c, 2c] e

igual a cero fuera de este segmento.

0,0 x c

ut x,0 V0 , c x 2c
0,2 x x

ux,0 x 0

Desarrollo por el mtodo de DALAMBERT

u x, t x x

u t x,0 x x

0,0 x
x x 0
0, x 0
 x , x 0
x
x , x 0

0,0 x c
v , c x 2c x[0,+]
0

0,2c x 0
x
0,c x 0
v0 ,2c x c


0, x 2c

0, x 2c
v ,2c x c


0

x 0,c x c
v , c x 2c
0

0,2c x

Se obtiene la solucin de la ecuacin de oscilaciones:

x at x at
u x, t x at x at
2

1 x at
u x, t x at x at d
2a x at

u x, t x at x at
Ejercicio 9(Piter Oneil), Pgina 256.
Ecuacin de Onda

Resuelva el problema con valores en la frontera usando la separacin de variables

2 2
= 8 , 0 < < 4 , > 4
2 2
(0, ) = (4, ) = 0 , 0
(, 0) = 2 ( 4),

(, 0) = 1, 0 4

Solucin
Aplicando el mtodo de separacin de variables:
(, ) = (). ()
(, ) = (). ,, ()
(, ) = ,, (). ()
Reemplazando en la ecuacin de onda:
2 2
=8 2
2

(). ,, () = 8 ,, (). ()
,, () 1 ,, ()
= =
() 8 ()
,, () = () ,, () = 8()
) = 2 > 0 , 0
) ,, () = 2 ()
,, () 2 () = 0
2 2 = 0
=+

() = 1 + 2

) ,, () = 82 ()
,, () 82 () = 0
2 82 = 0
 +
= 8
() = 3 8 + 4 8
:
(, ) = (1 + 2 )(3 8 + 4 8 )

(0, ) = (1 + 2 )(3 8 + 4 8 ) = 0
1 + 2 = 0 1 = 2

(4, ) = (1 4 + 2 4 )(3 8 + 4 8 ) = 0
1 4 + 2 4 = 0
1 4 1 4 = 0
1 ( 4 1 4 ) = 0
1 = 0, 2 = 0
(, ) = 0
) = 0 :
,, () = ()
,, () () = 0
2 = 0 = 0 ( = 2)
() = 1 0 + 2 0
() = 1 + 2

(0, ) = 0 (0). () = 0 (0) = 0 1 = 0

(4, ) = 0 (4). () = 0 (4) = 0 2 . 4 = 02 = 0
() = 1 + 2 = 0 +0.x=0
(, ) = 0

) = 2 < 0 , 0

) ,, () = 2 ()
,, () + 2 () = 0
 2 + 2 = 0
=+

() = 1 cos + 2 sin
: 1 = cos ; 2 = asin
() = cos cos + asin sin
:

() = cos( + )

) ,, () = 82 ()
,, () + 82 () = 0
2 + 82 = 0
+
= 8
() = 3 cos(8) + 4 sin(8)
: 3 = cos ; 4 = bsin
() = cos . cos(8) bsin . sin(8)
:

() = cos( + 8)

, :

(, ) = cos( + ) cos( + 8)

(, ) = cos( + ) cos( + 8) ; = .

(, ) = cos cos( + 8) ; 0 ; cos( + 8) 0
cos = 0

=
2

(4, ) = cos( + 4) cos( + 8) = 0

 0 ; cos( + 8) 0

cos( + 4) = 0

+ 4 = +
2

+ 4 = +
2 2

=
4

Por tanto una solucin de (, )

(, ) = cos ( + ) cos ( + 8 )
2 4 4

(, ) = (cos cos ( ) sin . sin ( ))
2 4 2 4

.(cos cos (8 ) sin . sin (8 ) )
4 4

= (sin ( )) (cos cos (8 ) sin . sin (8 ) )
4 4 4

= sin ( ) ( cos cos (8 ) + sin . sin (8 ) )
4 4 4

sin = cos =

(, ) = sin ( ) ( cos (8 ) + sin (8 ) )
4 4 4

Por superposicin se tiene:

(, ) = sin ( ) [ cos (8 ) + sin (8 )] . . ()
4 4 4
=1

Usando condiciones iniciales para hallar y



(, 0) = sin ( ) [ cos(0) + sin(0)] = 2 ( 4)
4
=1

sin ( ) = 2 ( 4)
4
=1

Esta es una funcin en serie de Fourier, la cual es impar en el intervalo (-

4,4).su coeficiente est dado por:

1 4 2
= ( 4) sin ( )
2 0 4

1 4 3 1 4 2
= sin ( ) 4 sin ( )
2 0 4 2 0 4

Calculando:

4

1). = 3 sin ( )
0 4

4
= 3 = 3 2 ; = sin ( ) = cos ( )
4 4

4
4 3 4 4 2
=[ cos ( )] + 3 cos ( )
4 0
0 4

256 12 4 2
= cos() + cos ( )
0 4
4
256
12 2 4 8 4
= (1) + [ sin ( )] sin ( ) ]
4 0 0 4

4
256 +1
12 8 4 4 4
= (1) + [ [ cos ( )] + cos ( ) ]
4 0 0 4

4
256 +1
12 8 16 4 4
= (1) + [ [ cos() + sin ( )] ]
4 0
 256 12 8 16 16
= (1)+1 + [ ( (1) + 2 2 (0))]

256 12(8)(16)(1)
= (1)+1 +
3 3

4 4

1 = 4 sin ( ) = 4 2 sin ( )
2
0 4 0 4

4
4 2 4 4
1 = 4 [[ cos ( )] + 2 cos ( ) ]
4 0 0 4

4
64 8 4 4 4
1 = 4 [ cos() + [ sin ( )] sin ( ) ]
4 0 0 4

4
64 +1
8 16 4 4
1 = 4 [ (1) + [ sin(0) [ cos ( )] ]]
4 0

64 8 16 16
1 = 4 [ (1)+1 + [ 2 2 cos() + 2 2 ]]

64 128 128
1 = 4 [ (1)+1 + 3 3 (1) + 3 3 ]

Reemplazando y 1 en :

16 16 +1
96(1) 1 (1)+1 2
2
= [ (1) + ] + (4)(64) [[ + (1) + ]]
2 3 3 2 3 3 3 3

(1)+1 6(1) (1)+1 2 2

= 128 [ + 3 3 ] + 128 [[ + 3 3 (1) + 3 3 ]]

2(1)+1 8(1) 2
= 128 [ + 3 3 + 3 3]

Usando la otra condicin



(, ) = sin ( ) [8 cos (8 ) 8 sin (8 )]
4 4 4 4 4
=1

(, 0) = sin ( ) 8 =1
4 4
=1

2 4
8 = sin ( )
4 4 0 4

4
1
= sin ( )
28 4 0 4

4
4 2 8 8
= [ cos ( )] = cos() +
8 4 0 8 8

8 8
= (1)+1 +
8 8

Reemplazando y en (), se tiene la solucion general

2(1)+1 8(1)
(, ) = sin ( ) [[128 [ + 3 3
4
=1

2 8 8
+ ]] cos (8 ) + [ (1)+1 + ] sin (8 )]
3 3 4 8 8 4

EJERCICIO 21.
Considere el problema con valores en la frontera:

= 9 + cos(), 0 < < 4, >0

{ (0, ) = (4, ) = 0, 0
(, 0) = (4 ), (, 0) = 0, 04
 a) Escriba una solucin en serie:

1) Hallemos la solucin para el problema homogneo

= 9 , 0 < < 4, >0

{ (0, ) = (4, ) = 0, 0
(, 0) = (4 ), (,
0) = 0, 04

Solucin por separacin de variables

(, ) = ()()

(0, ) = (0)() = 0 (0) = 0

(4, ) = (4)() = 0 (4) = 0

(, 0) = ()(0) = 0 (0) = 0

() 1 ()
= =
() 9 ()

() + () = 0 , (0) = (4) = 0

Para = 0, = 2 < 0, (, )

Para = 2 > 0:

() + 2 () = 0 2 + 2 = 0 = () = 1 + 2

(0) = 1 = 0 () = 2 (4) = 2 4

Tomemos: 2 0 4 = 0 4 = , = 1,2,

2
2 2
= =
16

() = 2 ( ) (1 )
4

92 2 3
() + 9 () = 0 , (0) = 0 2 + 16
= 0 = 4
 3 3
() = 1 cos ( ) + 2 ( )
4 4
3 3 3 3 3
() = 1 ( ) + 2 cos ( ) (0) = 2 = 0
4 4 4 4 4

2 = 0

3
: () = 1 cos ( ) (2 )
4

3
De (1 ) (2 ) tenemos: (, ) = cos (4 ) ( 4 )

Estas funciones satisfacen la ecuacin de onda, las condiciones de frontera y

condicin (, 0) = 0, excepto la condicin: (, 0) = (4 )
3
Por el principio de sustitucin tenemos: (, ) =
=1 cos (4 ) ( 4 )

(, 0) = ( ) = (4 )
4
=1

2 4
Utilizando serie de Fourier: = 4 0 (4 ) ( 4 )

2 4 2 4
= 4 ( ) 2 ( )
4 0 4 4 0 4
4
2 4
= 2 ( ) 2 ( )
0 4 4 0 4

Integracin por partes:

4
0 ( 4 )

= =
{ = ( ) = 4 cos ( )

4 4

4
64
( ) = (1)
0 4
 4
0 2 ( 4 )

= 2 = 2
{ = ( ) = 4 cos ( )
4 4

4 2 8
2 ( ) = cos ( ) + cos( )
4 4 4

= =
{ = ( ) = 4 sen ( )
4 4

4 2 32 128
2 ( ) = cos ( ) + 2 2 ( ) + 3 3 cos( )
4 4 4 4
4
4 2 32 128 4
2 ( ) = cos ( ) + 2 2 ( ) + 3 3 cos( )|
0 4 4 4 4 0

4
64 128 128
2 ( ) = (1) + 3 3 (1) 3 3
0 4

Por tanto:

64
= [1 (1) ]
3 3
La solucin para el problema homogneo ser:

64
3
(, ) = [1 (1) ] cos ( ) ( )
3 3 4 4
=1

128 3 (2 1)
(, ) = cos ( (2 1)) ( )
(2 1)3 3 4 4
=1
2) Hallemos la solucin para el problema No homogneo
= 9 + cos(), 0 < < 4, >0
{ (0, ) = (4, ) = 0, 0
(, 0) = 0, (, 0) = 0, 0 4
Por separacin de variable se obtiene:


(, ) = () ( ) (3 )
4
=1

4

(, 0) = (0) ( ) = 0 (0) = 0. ( ) = 0 (0) = 0
4 0 4
=1

Anlogamente se obtiene: (0) = 0

Derivando (3 ) para luego remplazarlo en la ecuacin de onda no homognea, de

eso se obtiene:

9 2 2
[ () + ()] ( ) = cos()
16 4
=1

Utilizando la serie de Fourier se obtiene:

9 2 2 2 4
() + () = cos() ( )
16 4 0 4

1 41 1
= [ ( + ) + ( )]
2 0 2 4 2 4

9 2 2 2 2
() + () = ( + ) [1 (1) ]
16 (4 + ) ( 4)

Vamos a solucionarlo por variacin de parmetros:

:
() = + , {
:

1) Hallando
 9 3 3
() + 16 2 2 () = 0 = 1 cos ( ) + 2 ( )
4 4

2) Hallando

Recordemos:

= 1 1 + 2 2 + + , : + 1 1 + 2 2 + + = ()

1 , 2 , ,

() 1

= , () = | |
0 () 1 1 1

() = (0) 1(0)

Para nuestro caso:

3 3
= 1 1 + 2 2 1 = cos ( ) , 2 = ( )
4 4

3 3
cos ( ) ( ) 1 0 3 3
4 4
() = | 3 3 3 3
| (0) = |0 3 | = () =
4 4
( ) ( ) 4
4 4 4 4

()1 2 2
1 = () = ( + ) [1 (1) ]
() (4 + ) ( 4)

3
0 ) ( 3
1 () = (1)1+2
| 4 | = ( )
3 3 4
1 ( )
4 4
4 2 2 3
1 = ( + ) [1 (1) ] ( )
3 (4 + ) ( 4) 4

4 2 2 3
1 = ( + ) [1 (1) ] ( )
0 3 (4 + ) ( 4) 4

16 2 2 3
1 = 2 2
( + ) [1 (1) ][cos ( ) 1]
9 (4 + ) ( 4) 4
 3
()2 cos ( ) 0 3
2 = , 2 () = (1)2+2 | 4 | = ( )
() 3 3 4
( ) 1
4 4

4 2 2 3
2 = ( + ) [1 (1) ] ( )
0 3 (4 + ) ( 4) 4

16 2 2 ]
3
2 = ( + ) [1 (1) ( )
92 2 (4 + ) ( 4) 4

Luego:

= 1 1 + 2 2

16 2 2 ]
3 3
= { ( + ) [1 (1) [cos ( ) 1]} cos ( )
92 2 (4 + ) ( 4) 4 4
16 2 2 3 3
+{ 2 2( + ) [1 (1) ] ( )} ( )
9 (4 + ) ( 4) 4 4

() = +
3 3
= 1 cos ( ) + 2 ( )
4 4
16 2 2 3
+ { 2 2 ( + ) [1 (1) ] [cos ( )
9 (4 + ) ( 4) 4
3
1]} cos ( )
4
16 2 2 3 3
+{ 2 2( + ) [1 (1) ] ( )} ( )
9 (4 + ) ( 4) 4 4

(0) = 1 = 0 1 = 0

(0) = 0 2 = 0

La solucin general de la EDO es:

 16 2 2 3 3
() = { 2 2
( + ) [1 (1) ] [cos ( ) 1]} cos ( )
9 (4 + ) ( 4) 4 4
16 2 2 3 3
+{ 2 2( + ) [1 (1) ] ( )} ( )
9 (4 + ) ( 4) 4 4

La solucin del problema no homogneo seria:

16 2 2 3
(, ) = ({ 2 2
( + ) [1 (1) ] [cos ( )
9 (4 + ) ( 4) 4
=1
3
1]} cos ( )
4
16 2 2
+{ 2 2( + ) [1
9 (4 + ) ( 4)
3 3
(1) ] ( )} ( )) ( )
4 4 4

Esto implicara que la solucin de nuestro problema es:

(, ) = (, ) + (, )

LIBRO DE PETER ONEIL (PAG. 258)

EJERCICIO 23.
Considere el problema con valores en la frontera :
= 9 + , 0 < < 1, >0
{ (0, ) = (1, ) = 0, 0
(, 0) = 0, (, 0) = , 0 1

Escriba una solucin en serie

1) Hallemos la solucin para el problema homogneo

= 9 , 0 < < 1, >0

{ (0, ) = (1, ) = 0, 0
(, 0) = 0, (,
0) = , 0 1

Solucin por separacin de variables

(, ) = ()()
 (0, ) = (0)() = 0 (0) = 0

(1, ) = (1)() = 0 (1) = 0

(, 0) = ()(0) = 0 (0) = 0

() 1 ()
= =
() 9 ()

() + () = 0 , (0) = (1) = 0

Para = 0, = 2 < 0 (, )

Para = 2 > 0

() + 2 () = 0 2 + 2 = 0 = () = 1 + 2

(0) = 1 = 0 () = 2 (1) = 2 . 1

Tomemos: 2 0 = 0 = , = 1,2,

= 2 = 2 2 . () = 2 (). 1

() + () = 0 , (0) = 0 2 + 92 2 = 0 = 3

() = 1 cos(3 ) + 2 (3 ) (0) = 1 = 0 () = 2 (3 ) 2

De 1 2 tenemos: (, ) = (3 ) ()

Estas funciones satisfacen la ecuacin de onda, las condiciones de frontera y

Condicin inicial, excepto (, 0) =

Por el principio de sustitucin tenemos: (, ) =

=1 (3 ) ()

(, 0) = 3 ( ) =
=1

1
Utilizando serie de Fourier: 3 = 2 0 ()

Integracin por partes:

= =
{ = () = 1 cos()

 1 2
3 = 2 [ (1) ] = 2 2 (1)
3
La solucin para el problema homogneo ser:

2
(, ) = 2 2
(1) (3 ) ()
3
=1

2) Hallemos la solucin para el problema No homogneo

= 9 + , 0 < < 1, >0
{ (0, ) = (1, ) = 0, 0
(, 0) = 0, (, 0) = 0,
Por separacin de variable se obtiene:

(, ) = ()() 3
=1

1
(, 0) = (0)() = 0 (0) = 0. () = 0 (0) = 0
=1 0

Anlogamente se obtiene: (0) = 0

Derivando 3 para luego remplazarlo en la ecuacin de onda no homognea, de

eso se obtiene:

[ () + 92 2 ()] () =
=1

Utilizando la serie de Fourier se obtiene:

1
2(1)
() + 92 2 () = 2 () =
0

2 2
2(1)
() + 9 () =

vamos a solucionarlo por variacin de parmetros

:
() = + , {
:
1) Hallando

() + 92 2 () = 0 = 1 cos(3 ) + 2 (3 )

2) Hallando

Recordemos

= 1 1 + 2 2 + + , : + 1 1 + 2 2 + + = ()

1 , 2 , ,

() 1

= , () = | |
0 () 1 1 1

() = (0) 1(0)

Para nuestro caso

= 1 1 + 2 2 1 = cos(3 ) , 2 = (3 )

cos(3 ) (3 ) 1 0
() = | | (0) = | | = 3 () =
3(3 ) 3(3 ) 0 3
3

()1 2(1) 0 (3 )
1 = () = 1 () = (1)1+2 | |
() 1 3(3 )
= (3 )
2(1) 2(1) 2(1)
1 = 32 2 (3 ) 1 = 0 32 2 (3 ) = (cos(3 )
93 3
1) 1

()2 cos(3 ) 0
2 = , 2 () = (1)2+2 | | = (3 )
() 3(3 ) 1

2(1)
2 = (3 ) 2
32 2

2(1) 2(1)
= (3 ) = (3 ) 2
0 32 2 93 3

1 2

2(1) 2(1)
= (cos(3 ) 1) cos(3 ) (3 )(3 )
93 3 93 3
 2(1)
= ((6 ) cos(3 ))
93 3
Luego:

() = +
2(1)
= 1 cos(3 ) + 2 (3 ) + ((6 ) cos(3 ))
93 3
2(1) 2(1)
(0) = 1 + (1) = 0 1 =
93 3 93 3
2(1) 4(1)
(0) = 32 + 6 = 0 2 =
93 3 93 3
La solucin general de la EDO es:

2(1) 4(1) 2(1)

() = cos(3 ) (3 ) + ((6 ) cos(3 ))
93 3 93 3 93 3

4(1)
() = [cos(3 ) (3 ) + ((6 ) (3 ) )]
93 3
La solucin del problema no homogneo seria:

4(1)
(, ) = [cos(3 ) (3 ) + ((6 )
93 3
=1
(3 ) )]()

Esto implicara que la solucin de nuestro problema es:

(, ) = (, ) + (, )

Ejercicio 3(Peter Onill), Pgina 281.

Ecuacin de Onda

Determinar las siguientes caractersticas de la onda para el problema

tt = 72 xx, para - < x < , t > 0

(x, 0) = f(x) = cos(x)

t (x,0) = g (x) = 1 x2

Solucin

Como la ecuacin va a ser hiperblica, su ecuacin caracterstica es:

(dx)2- 49 (dt)2 = 0 ( )2 = 49 = 7

x = 7t + C1 x = -7t + C2

x - 7t = C1 x + 7t = C2

Ahora definamos el cambio de coordenadas:

x- 7t = x + 7t =

Ahora despejemos x e y:
1 1
X= 2 (+), y = 14 (- +)

Definamos U (,) = u (x (+), y (-+))

Calculemos las derivadas:

uxx = U x + U x + U x + U x = U + 2U + U

utt = U (-7) + U (7)

y utt = -7U (-7) + U (7) + 7U (-7) + U (7) = 49U - 2(49)U + 49U

Entonces

utt 49uxx = -4(49)U

En estas nuevas coordenadas la ecuacin de onda es U = 0

Esta es la forma canonca de la ecuacin de la onda, resolvamos esta ecuacin:

Primero (U) = 0

De esto tenemos que U es independiente de, luego: U = h(n)

Integrando obtenemos: U (,) = h(n) dn + F () es una constante de integracin
de la derivada parcial respecto a , adems h(n)dn es solo otra funcin de , la
que escribiremos como G(). As:

U () = G () + F ()

Donde F y G deben ser funciones de una variable con segunda derivada continua,
por lo dems son arbitrarias. Por tanto pasamos a las coordenadas originales
tenemos:

(x, t) = F(x-7t) + G(x+7t)

Esta ecuacin se llama la solucin DAlambert de la ecuacin de onda.

Utilizando las condiciones iniciales:

(x,0) = G(x) + F(x) = f(x)..... (1)

t(x, t) = -7F(x-7t) + 7G(x+7t)

t(x,0) = -7F(x) + 7G(x) = g(x). (2)

Integrando la ltima ecuacin:

x
0 ut (,0) d = -7 0 ()d + 70 ()d = 0 () d

-7F(x)-7(0) + 7G(x) G (0) = 0 () d

1
-7F(x) + G(x) = 7 0 () d - F (0) + G(0)

Sumando esta ecuacin con la ecuacin (1)

1
2G(x) = f (x) +7 0 () d - F (0) G (0)

1 1 1 1
De donde. G(x) = 2 f(x) + 14 0 () d - 2 F (0) + 2G (0)

Tambin de la ecuacin (1)

1 1 1 1
F(x) = f(x)-G(x) = 2 f (x) - 14 0 () d+2 F (0) - 2 G (0)

Ahora reemplazando en la solucin, tenemos que

U(x,t) = F (x-7t) + G(x + 7t)

1 1 7 1 1
= 2 f (x-7t) - 14 0 () d + 2 F(0) - 2 G(0)
 1 1 +7 1 1
+ 2 f(x+7t) + 14 0 ()d - 2 F(0) + 2 G(0)

1 1 7 1 1 +7
u(x,t) = 2 f(x-7t) - 14 0 ()d + 2 f (x+7t)+ 14 0 ()d

1 1 +7
u(x,t) = 2 f(x-7t) + f (x + 7t)) 14 7 ()d

Ahora como f(x) = cos(x) y g(x) = 1-X2

1 1 +7
u(x, t) = 2 (cos((x-7t))+ cos((x+7t))) + 14 7 (1-2)d

1 1 1
u(x, t) = 2 (cos((x-7t))+ cos((x+7t))) + 14(|+7 - 3 |+7
7 3 7
)

1 1 1
u(x,t) = (cos((x-7t))+cos((x+7t)))+14(14t-3(x3+ x2.7t+x.49t2+73.t3x3+ x2.7t
2

x.49t2+73.t3))

1 1 14 2 686 3
u(x,t) = 2 (cos((x-7t))+ cos((x+7t))) + 14 (14t - x t.- t)
3 3

1 x2t 49
u(x,t) = 2 (cos((x-7t))+ cos((x+7t))) + t - t
3 3

Pg. 281 seccin 4.4

Ejercicio 10
Determine las caractersticas de la ecuacin de onda para el problema
= 4 , < < , >0
{ 3
(, 0) = 1 , (, 0) = ( ), < <
Y escribir la solucin de DAlambert
Solucin:
La ecuacin de onda es de tipo Hiperblico
1) Hallemos sus caractersticas
(, ) = 1 , (, ) = 0 , (, ) = 4
2 = 1
= = 2 {
2 + = 2
Vamos a tomar de tal manera que cumpla :
+ 2 = 0 (, ) = 2
{ { , (, ) = (, )
2 = 0 (x, t) = 2 +
Como la EDP es del tipo Hiperblica su forma canonca ser:
 (, ) = 0
2(, ) + (, , , , ) = 0 {
(, , , , ) = 0
Entonces : 16 (, ) = 0 (, ) = 1 () + 2 ()
Volviendo a las variables originales: (, ) = 1 (2 ) + 2 (2 + ) (1)
La recta = 0 no es tangente a las curvas caractersticas de esto se deduce que
el problema tiene una nica solucin
Determinaremos la solucin, para ello se debe verificar que la solucin (1)
satisface las condiciones inciales
(, 0) = 1 (2) + 2 (2) = x 3 1
derivando: 21 (2) + 22 (2) = 3x 2 , 1 (2) + 2 (2) = ( )
21 (2) + 22 (2) = 3x 2 3x2 +2( )
De esto se tiene: { 2 (2) =
21 (2) + 22 (2) = 2( ) 4
2 3
3x + 2( ) 1
2 (2) = 2 = + ( ) +
4 2
3 (2 + )3 1
1 2 + 1
2 () = + ( ) + 2 (2 + ) = + ( ) +
16 2 16 2

3 1 (2 )3 1 2 1
2 () = ( ) 2 (2 ) = + ( )
16 2 16 2
La solucin de DAlambert ser:

1
(, ) = [(2 + 1)3 + (2 1)3 ]
16

Libro Peter Oneil Pg. 281

Ejercicio 19: Resuelva el problema

, para - < x < , x > 0
u tt c 2 u xx F x, t

u x ,0 f x

u t x,0 g x , para - < x <

Donde c= 2, f (x) = X senx g(x) = e-x

F(x, t) = xt
 Solucin
, para - < x < , t > 0
utt 4u xx xt

u x,0 xsenx
, para - < x < (I)

ut x,0 ,epara
x - < x <

La solucin de (I) es: ux, t vx, t x, t

En la que vx, t es la solucin del problema de cauchy de la ecuacin

homognea.

vtt 4v xx 0

v x,0 x xsenx (II)

v t x ,0 x e x

Y x, t es la solucin del problema de cauchy con ecuacin no

homognea.
tt 4v xx g x, t xt

x,0 0 (III)

t x,0 0

La solucin de problema de cauchy II

 x, t x ct 1 xct d
x ct
v x, t
1
2 2c

x 2t senx 2t x 2t senx 2t 1 x2t e d

x 2t
v x, t
1
2 4

x 2t

v x, t x 2t sen x 2t x 2t sen x 2t e
1 1
2 4 x 2t

v x, t
1
2 4

x 2t senx 2t x 2t senx 2t 1 e x2c e x2c

x 2t senx 2t x 2t senx 2t 1 e x2t 1 e x 2t

v x, t
1
2 4 4
Para la solucin del problema de Cauchy I se cumple la misma
transformacin:

x 2t x 2t

Del problema (II), tomando como el dominio del triangulo R formado

por las caractersticas que pasan por (x, t).

El cambio de variables se reducen a la ecuacin no homognea a la

forma

4c 2 g ,
2 2

El dominio R se transforma su vez en el dominio D integrado la

ecuacin anterior en D se obtiene:

' ' ' '

v ,
1
4c 2 g
D
2
,
2
d ' d '

Volviendo a las variables originales:

 x, t g x' , t 'dx' dt ' , ya que se cumple con las condiciones
1
2c
R
x,0) = 0 tx,0) = 0
1 t
x' t ' dx' dt '
x 2 t 2 y '
x, t
0
4 x 2t 2t '

x2 x 2 t 2 t '

x, t
1 t
0 2
t ' dt '
x 2 t 2 t '
4

x, t
1 t
8 0

t ' x 2t 2t ' x 2t 2t ' dt '
2 2

1 t x 4t 4t ' 4 xt 4 xt'8tt '

2 2 2

x, t
8 0
2 dt '

t '
x 4t 2
4t ' 2
4 xt 4 xt '8tt '

x, t t '8 xt 8tt 'dt '

1 t
8 0

x, t
1 t
8 0

8 x t '8 xt'2 dt '

x, t
1 t
8 0

t ' xtt' xt'2 dt '

x, t xtt' xt' dt '

t
2
0

t t
xtt'2 xt'3 xtt2 xt3 xt 3 xt3
x, t
2 0
3 0
2 3 2 3

3 xt3 2 xt3 xt 3
x, t
6 6

xt3
x, t
6

Sabiendo que la solucin de I es:

 ux, t vx, t x, t

1
x 2t senx 2t x 2t senx 2t
2

1 x 2t 1 x 2t xt 3
e e
4 4 6

Ejercicio 19 de Libro Peter Oneil Pg. 316

Resuelva el problema:

u 2
9 u2
t x

, para 0 < x < 5 (I)

u x,0 0

u x, t 0, u 5, t ,3para t 0

Constryase la funcin
Donde T1 = 0, T2 = 3, L = 5
u x T1 T2 T1 x
e

v x 0 3 0 x
5

v x
3
x
5

Y consideramos el problema:
, para 0 < x < l, t > 0
t 9 xx

3,xpara 0 < x < l . (II)

x,0 0 v x
5

x, t 0, l , t
, t e> 0
La solucin de (II) es: ux, t vx x, t donde x, t es solucin de II

Hallando x, t por el mtodo de separacin de variables:

x, t f xg t

Luego:

x, t 2
f x g ' t y 2 f ' ' x g t
t x

Reemplazando estos en t 9 xx

f xg ' t y9 f ' ' xg t

g ' t f ' ' x f ' ' x g ' t

9 g t g x f x 9 g t

A. Tomando la 1 ecuacin:

f ' ' x f x 0,0 x 5 (1)

Como:

0, t f 0 g t 0 f 0 0
5, t f 5g t 0 f 5 0 (2)

Resolviendo (1) se analiza los casos:

i) si 0 la solucin es:

F x c1e x
c2 e x

F 0 c1 c2 c1 0 c2 0

F x 0, esta es la solucin que nos interesa

ii) si 0 reemplazamos en 1:
(Multiplicidad 2)
 f ' ' x 0 r 2 0 r 0

La solucin general es: f x c1e 0 x c2 xe0 x

F x c1 c2 x

F 0 c1 c1 0

F L c1 c2 L 0

0 L.c 2 0

c2 0

F x 0 , esta es la solucin que nos interesa.

iii) si 0 hacemos 2 , remplazamos en (1):

f ' ' x f x 0

f ' ' x 2 f x 0

f ' ' x 2 f x 0

Su ecuacin caracterstica es:

r2 2 0

Entonces la solucin es:

f x c1 cos wx c2 senwx

Luego de (2)
 f 0 c1 0

f 5 c1 cos 5 c2 sen5 0

0cos 5 c2 sen5 0

Si c2 0 f x 0 y esta es la solucin que no se quiere.

n
sen5 0 5 n , n 1,2,...
5

Los valores propios son:

n
n
5

n
Haciendo c 2 1 f n x sen x (3)
5

B. Considerando ahora la ecuacin: g ' t 9g t 0

Su ecuacin caracterstica es: R 9 0 R 9

Su solucin es g t c2 e 9 .(4)

De (3) y (4), tenemos la solucin:

9 n 2 2
n
n x, t cn e
t
25
sen x (5)
5

Por superposicin de soluciones a (5):

9 n 2 2

n
x, t c n e
t
25
sen x (6)
n 1 5

Usando la condicin:
 x,0
3x Para 0 x 5
5

En (6)

n
x,0 c n sen
3
x x
n 1 5 5

2 5 3x n
cn sen
5 5
0 5
xdx

6 5 n
cn
5 0
xsen
5
xdx

u x du d x

n n 5 n n
d v sen xd x v sen xdx sen xdx
5 5 n 5 5
5 n
v cos x
n 5

6 5 5 5 5
c n
5 n
x cos
n
x
n cos
n
sdx

n
5
6 5 5 25
xcons x 2 2 sen x
5 n n n 5 0

6 5
5cons.n 25 sen.n
5 n n
2 2

6 5
cn cos n
6
1n
25 n n
 9 n 2 2

n n
x, t 1 e
6 25
t
sen x
n 1 n 5

9 n 2 2
n n
ux, t x 1 e
3 6 25
t
sen x
5 n 5

Pg.316 Problema de Calor

Ejercicio 6
Hallar la solucin del problema con valores de frontera

= 4 , 0 < < 3, >0

{ (0, ) = (3, ) = 0, 0
2
(, 0) = , 03

Hallemos la solucin por separacin de variables

(, ) = ()()

(0, ) = (0)() = 0 (0) = 0

(3, ) = (3)() = 0 (3) = 0

() 1 ()
= =
() 4 ()

1) () + () = 0 , (0) = (3) = 0

Para = 0, = 2 < 0 (, )

Para = 2 > 0

() + 2 () = 0 2 + 2 = 0 = () = 1 + 2

(0) = 2 = 0 () = 1 (3) = 1 3 = 0

Tomemos: 1 0 3 = 0 = , = 1,2,
3

2 2
= 2 = . 1 = 1 () = ( ) . 1
9 3
 2 2
2 2 () 2 2
2) () + () = 0 () = () = 9 2
9 9

2 2

De 1 2 tenemos: (, ) = ( ) 9
3

Estas funciones satisfacen la ecuacin de calor y las condiciones de frontera

excepto la condicin inicial, (, 0) = 2
2 2

Por el principio de sustitucin tenemos: (, ) =
=1 ( ) 9
3

(, 0) = ( ) = 2
3
=1

2 3
Utilizando serie de Fourier: = 3 0 2 ( )
3

Integracin por partes:

= 2 = 2
{ = ( ) = 3 sen ( )
3 3

= =
4 3 3
= sen ( ) { = ( ) = cos ( )
0 3 3 3

36
= (1)
2 2
La solucin para este problema ser:

36
2 2
(, ) = (1) ( ) 9
2 2 3
=1

15) Use la idea del problema 13 para resolver:

 = ( + 6 ) 0 < < 4, > 0
{ (0, ) = (4, ) = 0 0
(, 0) = 1

Solucin:
(Por ejercicio 13):
La ecuacin diferencial parcial puede ser transformada en una ecuacin de calor
estndar eligiendo apropiadamente y y haciendo (, ) = + (, )

(0, ) = (0, ) = 0 (0, ) = 0

(4, ) = 4+ (4, ) = 0 (4, ) = 0
(, 0) = (, 0) = 1 (, 0) =

= + + +
= + + +
= 2 + + + + + + +

+ + + =
2 + + + + + + + + 6 + + 6 +

(1 ) = + ( + 6) + ( 2 + 6 )
= 6, =0
1
=
7

1
= 0 < < 4, > 0
{ 7
(0, ) = (4, ) = 0
(, 0) =

Resolvemos por separacin de variables:

Sustituimos (, ) = ()() en la ecuacin de calor, as obtenemos:
7 () ()
=
() ()
()
+ () = 0
2
Para = 0, = , la solucin es trivial.
Para = 2 > 0, se tiene:
() = 1 cos() + 2 ()
(0) = 2 = 0 () = 2 ()

(4) = 2 (4) = 0, supongamos que 2 0 (4) = = 4
2 2
= () = 2 ( 4 )
16

() + () = 0
2 2
() = 1 112

2 2
(, ) = 112 ( )
4
Por el principio de superposicin tenemos:

2 2

(, ) =
=1 112 ( 4 )


(, 0) = ( ) = 6
4
=1
Por series de Fourier tenemos:
1 4 6
= ( )
2 0 4
6 6
Integrando por partes: = = 6
4
= ( ) = cos( )
4 4

2
= (1 24 (1) )
2 2
576

2 2 2
24 (1)
(, ) = (1 ) 112 ( )
2 2 576 4
=1
Esto implica que:

2 2 2
6 24 (1)
(, ) = (1 ) 112 ( )
2 2 576 4
=1
 EJERCICIO SOBRE TRASNFORMADA DE FOURIER

2 y c 2 2 y
, x , t 0
t 2 x 2

y t 2 , para 2 t 2
y x,0 f t
0, parat 2, parat 2

y
x,0 0, para xz
t

Solucin: mediante transformada

Debido a que x vara sobre toda la recta, podemos intentar la transformada de

Fourier en la variable x.

Para esto transformaremos y(x, t) como una funcin de x, dejando t como

parmetro.

2 y 2 y
F 2 c 2 F 2 ............ *
t x

Calculando

2 y y
2
F 2 x, t e ix dx
t t 2

2 y
2
y x, t e ix dx
t

2 y
F 2 2 y, t ...........................( I )
t t

Luego:
 2 y y
2
F 2 x, t e ix dx
t t
2

Por la propiedad de la derivada de la transformada.

F F n i f
n

Quedando (*) al reemplazar (I) y (II)

y , t 2 y , t
t 2

y , t 2 y , t 0.............................. * *
t 2

Su ecuacin caracterstica de (**)

Transformada de Forier de onda y calor

Newman Dirichit

R 2 c 2 2 0

R ci

La Solucin de (**) es:

y, t a cosct b senct

y ,0 a cos0 b sen0

y ,0 a

y ,0 F y x,0 f

 4 t 2 , para 2 t 2
Calculando f de f t
0, t 2, t 2

f f t e it dt 0e it dt 4 t 2 e it dt 0e it dt

2

2

0 4 t 2 e it dt
2

2 2
4e it dt t 2 e it dt
2 2

4 2 2

i 2
ie it dt t 2 e it dt
2

2
4 (***)
e it (* * *)
i 2

Calculando (***)

t e it dt
2
Por integracin por partes

u t 2 , du 2tdt

1 it
dv e it , v e it dt e
i

t 2 it 2 it
t
it
2
e dt e te dt
i i

t 2 it 2

i
e
i te it dt .(*****)

(****)
 2 it t 2 it 1 t it t it

it
t e dt e e dt e e ............(1)
i i i 2

u t , du dt
t it
dv e it v e it e
i

2
4 it t 2 it 2 t it
f
1
e e e 2 e it
i i i i 2

2
4 it t 2 it 2t it 2
e e e 3 e it
i i i i 2

4 it t 2 it 2t it 2
e e e 3 e it
i i i i

4 2 e it t 2 2 e it 2t i e it 2e it
i 3

4 2

t 2 2 2ti 2 it
e
2

i 3 2

4 2

4 2 4i 2 4 2 4 2 2i 2 2it
e e 2 i
i 3

4i 2 4i 2 e 2i e 2i
i 3

8i 2i

i
2

e e 2 i
2

4 e 2 i e 2 i 4
2 sen2 cos2
2 2

f 2 sen2 cos2
4


y ,0 f

Sabiendo que:

y
x,0 0
t

y
x, t ca senct cb senct
t

y
x, t cb
t

y
F x,0 cb
t

F 0 cb

0 cb b 0

f , t f cosct

Aplicando la inversa de y, t

y x, t f cosct e ix d w
1
2

sen2 cos2 cosct e ix d w

y x, t
1 4
2 3

sen2 cos2
cosct cos x i cos x d w

y x, t
1
2 3

Como la funcin desplazamiento forma valores reales, debemos

tomas la parte real de esta integral

sen2 cos2
y x, t cosct cos wx d w
1
2
3

PAG 337. PROBLEMA DE DIRICHLET PARA EL RECTNGULO

Ejercicio 3

+ = 0, 0 < < 1, 0 < < 4

{ (0, ) = (1, ) = 0, 04

(, 0) = 0, (, 4) = ( ) , 01
2
Por separacin de variable:

(, ) = ()()

(0, ) = (0)() = 0 (0) = 0

(1, ) = (1)() = 0 (1) = 0

(, 0) = ()(0) = 0 (0) = 0
 () ()
= =
() ()

1) () + () = 0 , (0) = (1) = 0

Para = 0, = 2 < 0 (, )

Para = 2 > 0

() + 2 () = 0 2 + 2 = 0 = () = 1 + 2

(0) = 1 = 0 () = 2 (1) = 2 = 0

Tomemos: 2 0 = 0 = , = 1,2,

= 2 = 2 2 . 1 = 1 () = (). 1

2) () 2 2 () = 0 , (0) = 0 2 2 2 = 0 =

() = 1 + 2 (0) = 1 + 2 = 0 ()

= 21 ( ) = 21 ()
2

() = () 2

De 1 2 tenemos: (, ) = ( ) ()

Estas funciones armnicas satisfacen la ecuacin de La place y las condiciones de

frontera excepto, (, 4) = ( 2 )

Por el principio de sustitucin tenemos: (, ) =

=1 ( ) ()

(, 4) = ( ) (4) = ( )
2
=1

1
Utilizando serie de Fourier: (4) = 2 0 ( 2 ) ( ) =

Para hacer fcil la integral vamos a utilizar la transformacin trigonomtrica

+ 1
( ) . ( ) = [ ]
2 2 2
Sea
 +
= + = (1+2) (12)
2 2
{ = ; =
2 2
= = 2
2

1 1 (1+2) (12)
Esto implicara: = 2 0 [ ]
2 2 2

1
(1 + 2)
1 =
0 2
= =
{ = (1 + 2) 2 (1 + 2)
= cos
2 (1 + 2) 2

4(1)
1 =
2 (1 + 2)2
1 (12) 4(1)
Anlogamente 2 = 0 = 2 (12)2
2

32(1)+1
1 2 = =
2 (12)2 (1+2)2
32(1)+1
2 (4)(12)2 (1+2)2

La solucin general ser:

32(1)+1
(, ) = ( ) ()
2 (4)(1 2)2 (1 + 2)2
=1

PAG 363. PROBLEMA DE NEUMANN PARA UN RECTNGULO

Ejercicio 3 [hemos igualado la condicin (, 0) = 0 ]
+ = 0, 0 < < , 0 < <
{ (, 0) = 0; (, ) = 6 3, 0
(0, ) = 0, (, ) = 0 , 0

Por separacin de variable:

(, ) = ()()

(0, ) = (0)() = 0 (0) = 0

(1, ) = ()() = 0 () = 0
 (, 0) = ()(0) = 0 (0) = 0

() ()
= =
() ()

1) () + () = 0 , (0) = () = 0

Para = 0, = 2 < 0 (, )

Para = 2 > 0

() + 2 () = 0 2 + 2 = 0 = () = 1 + 2

(0) = 2 = 0 () = 1 () = 1 = 0

Tomemos: 1 0 = 0 = , = 0,1,2,

= 2 = 2 . 1 = 1 () = (). 1

2) () 2 () = 0 , (0) = 0 2 2 = 0 =

() = 1 + 2 () = 1 2 = 0 (0) = (1 2 ) = 0

Para = 0 () =
+
Para 1 1 = 2 () = 21 ( ) = 21 () () =
2
(). 2

De 1 2 tenemos: (, ) = ( ) ()

Estas funciones son armnicas las cuales satisfacen la ecuacin de Laplace y las
condiciones de frontera excepto, (, ) = 6 3

Por el principio de superposicin tenemos: (, ) = 0 +

=1 ( ) ()

(, ) = ( ) () = 6 3
=1

2
Utilizando serie de Fourier: () = 0 (6 3) ( ) =

2 2
= 6( ) 3() = 1 2
0 0
 2 6
1 = 6() = 2 ((1) 1) ; 2 = 0
0
12 12
Esto implicara: () = 2 ((1) 1) = 3 () ((1) 1)

La solucin general ser:

12
(, ) = 0 + ((1) 1)( ) ()
3 ()
=1