Tema 4

Prof.
Jos Luis Quintero 1

TEMA 4
ECUACIONES DIFERENCIALES PARCIALES

4.1. INTRODUCCIN

En este tema se vern procedimientos para resolver ecuaciones en derivadas
parciales que surgen con frecuencia en problemas donde aparecen vibraciones,
potenciales y distribuciones de temperatura. Se llaman problemas de valores en la
frontera y se describen mediante ecuaciones en derivadas parciales de primer o
segundo orden, que son relativamente simples. En el caso de las ecuaciones de
segundo orden, lo que se hace es hallar las soluciones particulares de una ecuacin en
derivadas parciales reducindolas a dos o ms ecuaciones diferenciales ordinarias. Se
comenzar con el mtodo de separacin de variables para ecuaciones en derivadas
parciales lineales. La aplicacin de este mtodo recordar conceptos de temas
anteriores como valores propios, funciones propias y el desarrollo de una funcin en
una serie infinita de funciones ortogonales.

4.2. ESPACIO VECTORIAL CON PRODUCTO INTERIOR

Definicin 1. Un espacio vectorial real es un conjunto V (cuyos elementos se
llamarn vectores), con una operacin binaria de adicin y una multiplicacin escalar
que satisfacen las siguientes propiedades: ( , R)
a. Dado cualquier par de vectores x, y V , existe un nico vector x y V + , llamado
la suma de x y y
b. x , y V , x y z x y z ( ) ( ) + + = + +
c. x , y V , x y y x + = +
d. 0 x x 0 0 x x V V , + = + =
e. x x x x 0 ( ) V V , ( ) + =
f. Dado cualquier vector x V y cualquier nmero real , existe un nico vector
x V llamado el producto de y x
g. x y x y ( ) + = +
h. x x x ( ) + = +
i. x x ( ) ( ) =
j. x x 1. =

Definicin 2. Se dice que un subconjunto S de un espacio vectorial V, es un
subespacio de V, si x y , S y , R el vector x y S + .

Observacin 1. Es claro que { } 0 y V son subespacios de V.

Prof. Jos Luis Quintero 2
Definicin 3. Se dice que un producto interior est definido en un espacio vectorial
real V, si con cada par de vectores x e y en V se asocia un nmero real x y tal que
a. x y y x =
b. x y x y ( ) ( ) = para cada nmero real
c.
1 2 1 2
x x y x y x y ( ) + = +
d. x x 0 y x x 0 = si, y solamente si, x 0 = .

Definicin 4. Un espacio vectorial con un producto interior se conoce como un
espacio euclidiano o espacio vectorial con producto interior.

Si se aplica (a) a (b) y a (c), se ve que un producto interior tambin satisface
las propiedades x y x y ( ) ( ) = y
1 2 1 2
x y y x y x y ( ) + = + . En general
1 2 1 2 1 1 1 2 2 1 2 2
x x y y x y x y x y x y
1 2 1 2 1 1 1 2 2 1 2 2
( ) ( ) ( ) ( ) ( ) ( ) + + = + + +
donde
1 2 1
, , y
2
son escalares arbitrarios.

Ejemplo 1. Sean x
1 n
(x ,..., x ) = y y
1 n
(y ,..., y ) = vectores en
n
R , y defnase a x y
por x y
1 1 n n
x y ... x y = + + . Entonces
n
R se convierte en un espacio euclidiano y, como
tal, se llama espacio euclidiano de dimensin n. En
2
R y
3
R este producto interior
no es otro que el producto punto usado en fsica, donde la definicin generalmente
se expresa geomtricamente como el producto de la longitud de x por la longitud de y
por el coseno del ngulo que forman ambos.

Ejemplo 2. Considere el espacio C a,b (

de todas las funciones continuas en a,b (

, con
f g definido por

f g
b
a
f(x)g(x)dx =
.
Pruebe que f g define un producto interior.
Solucin.
No es difcil demostrar que satisface todas las exigencias de un producto interior: Se
tiene que

b b
a a
f(x)g(x)dx f(x)g(x)dx =

,
siendo un nmero real, y

b b b
1 2 1 2
a a a
f (x) f (x) g(x)dx f (x)g(x)dx f (x)g(x)dx + = + (

.
Finalmente, al recordar que la integral de una funcin no negativa es no negativa, y
que la integral de una funcin continua no negativa es cero si, y solamente si, la
funcin es idnticamente cero se ve que

f f
b
2
a
(f(x)) dx 0 =
,
y que f f 0 = si, y solamente si, f 0 = .

En adelante, siempre que se haga referencia a
n
R o a C a,b (

como espacios
euclidianos, se considerar que se estn utilizando los productos interiores definidos en
los ejemplos anteriores, a menos que se mencione expresamente lo contrario.

Ejemplo 3. Sea r una funcin no negativa en C a,b (

que se anula, cuando ms, en un
nmero finito de puntos en el intervalo a,b (

. Entonces

f g
b
a
f(x)g(x)r(x)dx =

tambin define un producto interior en C a,b (

.

La funcin r se llama la funcin peso para este producto interior, y se
observar que cuando r es idnticamente igual a 1 se tendr la frmula definida en el
ejemplo anterior. Se encontrar este producto interior nuevamente cuando se estudien
problemas de valor en la frontera, para ecuaciones diferenciales.

4.3. LONGITUD. MEDIDA ANGULAR. DISTANCIA

Se definir longitud, medida angular y distancia, en trminos del producto
interior de un espacio euclidiano. Cada uno de estos conceptos tiene un significado
bien definido en el espacio euclidiano, de dimensin 2, y es razonable exigir que
cualquier definicin que se adopte, se reduzca a la ya usual en
2
R .

Sea x
1 2
(x , x ) = cualquier vector en
2
R . Entonces, la longitud de x indicada por
x , es el nmero real no negativo x
2 2
1 2
x x = + . Esta expresin puede replantearse
en trminos del producto interior en
2
R como x x x = y se tiene la siguiente
definicin:

Definicin 5. La longitud (o norma) de un vector x en un espacio euclidiano, se
define como el nmero real no negativo x x x = .

As, la longitud de un vector x
1 n
(x ,..., x ) = en
n
R , es x
2 2
1 n
x ... x = + + ,
mientras que la longitud de un vector f en C a,b (

es

f
1/ 2
b
2
a
(f(x)) dx
(
=
(

.

A continuacin, se observa que si x y y son dos vectores no nulos cualesquiera,
en
2
R , la frmula
x y

x y
cos , 0
= ,
es una consecuencia inmediata de la ley de los cosenos. La expresin
x y
x y

tambin es significativa en un espacio euclidiano arbitrario, hecho que sugiere
considerarla como un candidato razonable para la definicin general del cos . Sin
embargo, antes de proceder, debe establecerse la desigualdad
x y
x y
1 1

para cada par de vectores no nulos en un espacio euclidiano, ya que, por supuesto,
cualquier definicin de cos debe satisfacer la desigualdad 1 cos 1 . Este hecho
surgir como una consecuencia del siguiente resultado importante, conocida como la
desigualdad de Schwarz o de Cauchy-Schwarz.

TEOREMA 1. (Desigualdad de Schwarz) Si x y y son dos vectores cualesquiera en
un espacio euclidiano, entonces x y x x y y
2
( ) ( )( ) .
Demostracin.
Primero se observa que esta desigualdad es inmediata si x o y, es el vector cero, ya
que entonces ambos miembros de la desigualdad son cero. As, es suficiente
considerar el caso en que x y y son no nulos. Aqu se usan propiedades del producto
interior para desarrollar x y x y ( ) ( ) donde y son nmeros reales
arbitrarios. Por la propiedad d de la definicin 3 se tiene
x y x y x x x y y y
2 2
0 ( ) ( ) ( ) 2 ( ) ( ) = + ,
de donde x y x x y y
2 2
2 ( ) ( ) ( ) + . Ahora si y y = y x x = se obtiene
x x y y x y x x y y 2 ( ) 2( )( ) que es equivalente a x y x x y y . Ahora
bien, si se eleva al cuadrado, se obtiene lo que se quera demostrar.

En
n
R la desigualdad de Schwarz toma la forma
2
n n n
2 2
i i i i
i 1 i 1 i 1
x y x y
= = =
| | | | | |
| | |
| | |
\ \ \

y en C a,b (

se convierte en

2
b b b
2 2
a a a
f(x)g(x)dx (f(x)) dx (g(x)) dx
| | | | | |
| | |
\ \ \

.

La primera de estas desigualdades es vlida para cualquier coleccin
1 n
x ,..., x ,
1 n
y ,..., y de nmeros reales, y generalmente se le llama la desigualdad de Cauchy.
Es importante recordarla ya que es frecuentemente til en la deduccin de otras
desigualdades aritmticas.

En la notacin de la definicin 5, la desigualdad de Schwarz se convierte en
x y x y , y afirma que el valor absoluto del producto interior de dos vectores no
excede el producto de las longitudes de los vectores. As
x y
x y
1
,
siempre que x y y sean no nulos.

Ejemplo 4. Sean
1 n
a ,..., a nmeros reales positivos. Use la desigualdad de Cauchy
para demostrar que
2
1 n
1 n
1 1
(a ... a ) ... n
a a
| |
+ + + +
|
\
.
Solucin.
Sean
x y
1 n
1 n
1 1
( a ,..., a ) ,...,
a a
| |
| = =
|
\
.
De modo que
x x
1 n
a ... a = + + , y y x y
1 n
1 1
... , n
a a
= + + =
y aplicando la desigualdad de Cauchy se obtiene lo que se quera demostrar.

Ejemplo 5. Sean a, b y c, nmeros reales positivos, tales que a b c 1 + + = . Utilice la
desigualdad de Cauchy para demostrar que
1 1 1
1 1 1 8
a b c
| | | | | |

| | |
\ \ \
.
Solucin.
Sean
x y
1 1 1
( a, b, c) , ,
a b c
| |
= =
|
\
.
De modo que
x x a b c 1, = + + = y y x y
1 1 1
, 3
a b c
= + + =
y aplicando la desigualdad de Cauchy
1 1 1
9
a b c
+ + .
Por lo tanto manipulando se tiene

bc ac ab bc ac ab bc ac ab abc (1 a)(1 b)(1 c)
9 1 8 8 8
abc abc abc abc
1 a 1 b 1 c 1 1 1
8 1 1 1 8
a b c a b c
+ + + + + +

| | | | | | | | | | | |

| | | | | |
\ \ \ \ \ \

Ejemplo 6. Demuestre que la desigualdad
2
2 2
1 n 1 n
a ... a a ... a
n n
+ + + + | |

|
\

persiste para cualquier coleccin de nmeros reales
1 n
a ,..., a .
Solucin.
Sean
x y
1 n
a a
(1,...,1) ,...,
n n
| |
= =
|
\
.
Entonces
x x y y
2 2
1 n
2
a ... a
n,
n
+ +
= = , x y
1 n
a ... a
n
+ +
=
y aplicando la desigualdad de Cauchy se obtiene lo que se quera demostrar.

Ejemplo 7. Sea
1
f C 0,1 (

. Pruebe que

1/ 2
1 1
2 1 1 2
2 4 15
0 0
( t)f(t)dt f '(t) dt

(
`

)

.
Solucin.
Integrando por partes:
2
t t 1
2 2 2
u f(t) du f '(t)dt , dv t v = = = = .

1 1 1 1
2 2 2
t t t t t t 1
2 2 2 2 2 2 2
0 0 0 0
( t)f(t)dt ( )f(t) ( )f '(t)dt ( )f '(t)dt = =

.
Aplicando la desigualdad de Schwarz:

1/ 2 1/ 2
1 1 1
2 2
2 2 t t t t
2 2 2 2
0 0 0
( )f '(t)dt ( ) dt (f '(t)) dt
( (

( (

.

1 1 1
2 4 3 2 5 4 3
2 t t t t t t t t 6 15 10 1 1 1 1
2 2 4 2 4 20 8 12 20 8 12 120 120
0 0 0
( ) dt ( )dt ( )
+
= + = + = + = =

.
Por lo tanto:

1/ 2
1
2
2 1/ 2 2 t t 1 1 1 2
2 2 120 4 15
2 2 15 4 15
0
( ) dt ( )
(
= = = =
(

.
En consecuencia se demuestra que:

1/ 2
1 1
2 1 1 2
2 4 15
0 0
( t)f(t)dt f '(t) dt

(
`

)

.

Definicin 6. Si x y y son vectores no nulos en un espacio euclidiano, se define el
coseno del ngulo entre ellos como
x y
x y
cos

= .

Queda finalmente pendiente por definir la distancia entre dos puntos
cualesquiera (esto es, vectores) en un espacio euclidiano. Nuevamente esto se hace
sencillamente copiando la definicin de
2
R donde la distancia entre x y y es la longitud
del vector x y . As se establece la definicin de distancia:

Definicin 7. La distancia entre dos vectores x y y en un espacio euclidiano es, por
definicin, x y x y d( , ) = .

Pero es sta una definicin razonable del trmino distancia? Con objeto de
contestar esta pregunta, primero debe decidirse qu propiedades se requieren de la
distancia, en general. En ese aspecto, los matemticos estn de acuerdo, y lo han
decidido como sigue: la distancia entre dos puntos debe ser un nmero real no
negativo que es cero si, y solamente si, los puntos coinciden. Debe ser independiente
del orden en el que se consideren los puntos y, finalmente, la desigualdad del
tringulo, famosa en la geometra plana, debe satisfacerse.

As, con el objeto de justificar el uso del trmino distancia en la definicin 7
se debe demostrar que d(x,y) es un nmero real que satisface las siguientes
propiedades:
a. x y d( , ) 0
b. x y d( , ) 0 = si, y solamente si, x y =
c. x y y x d( , ) d( , ) =
d. x y y z x z d( , ) d( , ) d( , ) + , x, y y z tres vectores cualesquiera.

Las tres primeras de estas propiedades se deducen inmediatamente de la
definicin de longitud y de los axiomas que gobiernan un producto interior. La ltima,
sin embargo, no es tan obvia. Para demostrarla, primero se establecer una
desigualdad que es de cierta importancia por s misma.

LEMA 1. Si x y y son vectores arbitrarios en un espacio euclidiano, entonces
x y x y + + .
Demostracin.
( )
x y x y x y x x x y y y
x x x x y y y y x x y y x x y y
x y
1 / 2 1/ 2
1/ 2
2 1/ 2
( ) ( ) ( ) 2( ) ( )
( ) 2 ( )
+ = + + = + + ( (

(
(
+ + = + = +
(

+

La desigualdad del tringulo sigue inmediatamente de este resultado. En
efecto, x z x y y z x y y z ( ) ( ) = + + , que precisamente es lo que se tena
que demostrar. Finalmente, se observa que la funcin distancia definida antes tambin
goza de las propiedades notables siguientes:
a. x y x y d( , ) d( , ) = para todo nmero real
b. x z y z x y d( , ) d( , ) + + =

4.4. ORTOGONALIDAD. ORTOGONALIZACIN

Definicin 8. Se dice que dos vectores en un espacio euclidiano son ortogonales o
perpendiculares, si el coseno del ngulo que forman, vale cero.

En seguida se generalizar algo ms la nocin de ortogonalidad, pero primero
se probar un teorema particularmente notable.

TEOREMA 2. (Pitgoras) Dos vectores x y y en un espacio euclidiano se dice que son
ortogonales si, y solamente si, x y x y
2 2 2
+ = + .
Demostracin.
x y x y x y x x x y y y x x y y
2 2 2
( ) ( ) 2( ) 2( ) + = + + = + + = + +
As x y x y
2 2 2
+ = + si, y solamente si, x y 0 = , como se afirm.

Ejemplo 8. Sean
1 n
x x , ..., vectores mutuamente perpendiculares en un espacio
euclidiano. Demuestre que
1 n 1 n
x x x x
2 2 2
... ... + + = + + .
Solucin. Se sabe que
1 n 1 n 1 n
i i i j i i j
x x x x x x
x x x x x x x
2
n n n n n n
2
i 1 i 1 j 1 i 1 i 1 j 1
i j i j
... ( ... ) ( ... )
( ) ( ) ( )
= = = = = =

+ + = + + + +
= + = +

As
1 n i 1 n
x x x x x
n
2 2 2 2
i 1
... ...
=
+ + = = + +

ya que
i j
x x 0, i j = .

Dicho esto, se contina la discusin, dando la siguiente definicin:

Definicin 9. Se dice que un conjunto de vectores
1 2 i
x x x , ,..., ,... en un espacio
euclidiano, es un conjunto ortogonal si
i
x 0 para todo i, y que
i j
x x 0 = siempre
que i j . Si, adems,
i i
x x 1 = para cada i, se dice que el conjunto es ortonormal.

As, un conjunto ortogonal es un conjunto de vectores no nulos perpendiculares
entre s, mientras que un conjunto ortonormal es un conjunto ortogonal en el cual cada
uno de los vectores es de longitud unitaria. Por economa de notacin, cuando se
discutan conjuntos ortonormales, las propiedades anteriores se combinarn
frecuentemente escribiendo
i j

x x

ij
0 si i j
1 si i j

= =

=
.
El smbolo
ij
que se introduce aqu, se le llama la delta de Kronecker.

TEOREMA 3. Un conjunto ortogonal de vectores { }
1 2 n
x x x , ,..., es linealmente
independiente.
Demostracin.
Se sabe de cursos anteriores de lgebra Lineal que un conjunto de vectores
{ }
1 2 n
x x x , ,..., es linealmente independiente si
1 2 n
x x x 0
1 2 n 1 2 n
... ... 0 + + + = = = = = .
Se tiene una combinacin lineal de los vectores ortogonales:
1 2 n
x x x 0
1 2 n
... + + + =
Haciendo producto escalar en la expresin anterior por un vector del conjunto
i
x se
tiene
1 i 2 i i i n i i
x x x x x x x x x
2
1 2 i n i
... ... 0 0 + + + + + = = .
De modo que { }
i
0 i 1,2,...,n = . Esto implica que el conjunto ortogonal de vectores
{ }
1 2 n
x x x , ,..., es linealmente independiente (forman una base) por lo tanto forman un
sistema completo.

Ejemplo 9. En
3
R los vectores (1, 0, 0), (0,2, 0), (0, 0, 1/ 2) forman un conjunto
ortogonal, mientras que los vectores de base cannica (1,0,0), (0,1,0) y (0,0,1)
forman un conjunto ortonormal. En forma ms general, el conjunto que consiste de los
vectores de la base en
n
R es ortonormal.

Ejemplo 10. Se define un polinomio trigonomtrico de grado 2n 1 + como una
expresin de la forma
0
1 2 n 1 2 n
a
f(x) a cos(x) a cos(2x) ... a cos(nx) b sen(x) b sen(2x) ... b sen(nx)
2
= + + + + + + + +
donde
0 n
a ,...,b son nmeros reales, y
n
a 0 o
n
b 0 o ambos. Sea
n
T el conjunto de
todos los polinomios trigonomtricos de grado 2n 1 + , junto con el polinomio cero.
Haciendo que
n
T sea un espacio euclidiano al definir la adicin y la multiplicacin
escalar de los polinomios trigonomtricos, trmino a trmino como con los polinomios
ordinarios y un producto interior dado por

f g f(x)g(x)dx
,
el conjunto de funciones 1, cos(x), sen(x), ..., cos(nx), sen(nx) es un conjunto ortogonal
en
n
T ya que, para enteros m y n no negativos

sen(mx)sen(nx)dx 0,

sen(mx)cos(nx)dx 0,

cos(mx)cos(nx)dx 0, si m n
.
Para normalizar este conjunto se observa que

1dx 2
y

2 2
sen (mx)dx cos (mx)dx , si m 0

= = >

.
Por tanto, las funciones
1 cos(x) cos(nx)
, ,..., ,
2
sen(x) sen(nx)
,...,

forman un conjunto ortonormal en
n
T .

Ejemplo 11. Se deduce del ejemplo precedente que el conjunto (infinito) 1, cos(x),
sen(x), , cos(nx), sen(nx), es ortogonal en C , (

.

En lugar de empezar con la situacin ms general, se introducir el proceso de
ortogonalizacin con dos ejemplos. El primero es extrado de
2
R donde se considera
un par de vectores linealmente independientes
1
x y
2
x . Entonces
1
x y
2
x forman
una base para
2
R y, en este caso, el problema se convierte en la sustitucin de
1
x y
2
x por una base ortogonal
1
e y
2
e construida en trminos de
1
x y
2
x de alguna
manera razonable. Una forma sencilla de resolver el problema es simplemente tomar
=
1 1
e x y luego se expresa
2
e en la forma
2 2 1
e x e = , y luego se determina de
modo que la condicin de ortogonalidad
2 1
e e 0 = quede satisfecha. Esto da la
ecuacin
2 1 1 1
x e e e ( ) 0 = , y en consecuencia, el valor de es
2 1
1 1
x e
e e
.
Con esto
2
e se ha determinado en trminos de
1 1
x e ( ) = y
2
x , y la base
1
x y
2
x en
2
R se ha ortogonalizado.

Ejemplo 12. Si
1
x (1,1) = y
2
x (0,1) = , entonces
(0,1) (1,1) 1
(1,1) (1,1) 2
= =
,
y as
1
e (1,1), =
2
e
1 1
2 2
( , ) = .

Como otro ejemplo, se ortogonaliza una base arbitraria
1
x ,
2
x y
3
x en
3
R . El
procedimiento esencialmente es el mismo usado antes en
2
R , y se inicia eligiendo
1 1
e x = . El segundo paso consiste en determinar
2
e de acuerdo con el par de
ecuaciones
2 1
e e 0, =
2 2 1
e x e = , lo que da nuevamente
2 1
1 1
x e
e e
.

Est claro que
2
e no es el vector cero, y tambin lo es que
1
e y
2
e
pertenecen al subespacio de
3
R generado por
1
x y
2
x . En consecuencia
1 2
e e ( , ) es
un subespacio de
1 2
x x ( , ) . Ms an, ya que los vectores ortogonales
1
e y
2
e son
linealmente independientes,
1 2
e e ( , ) tiene la misma dimensin que
1 2
x x ( , ) . As,
1 2 1 2
e e x x ( , ) ( , ) = .

Combinado con el hecho de que
1
x ,
2
x y
3
x forman una base para
3
R , esta
igualdad implica que
3
x no pertenece al subespacio de
3
R generado por
1
e y
2
e . El
proceso de ortogonalizacin debe completarse aceptando que
3
e sea la componente
de
3
x , perpendicular al subespacio
1 2
e e ( , ) . As, se hace
3 3 1 2
e x e e
1 2
= , y se
encuentran
1
y
2
por medio de las condiciones de la ortogonalidad
1 2 1 3
e e e e =
2 3
e e 0 = = . Ellas dan el par de ecuaciones dadas por
3 1 1 1
x e e e
1
0 ( ), =
3 2 2 2
x e e e
2
0 ( ) = de donde
3 1 3 2
1 1 2 2
x e x e

e e e e
1 2
,

= =

.

Con esto se completa el proceso de ortogonalizacin de la base
1
x ,
2
x y
3
x .
Ahora se puede disponer de la situacin general requerida para ortogonalizar un
conjunto arbitrario de vectores linealmente independientes
1 2
x x , , ... en un espacio
euclidiano. Se va a presentar un teorema que puede resumir este proceso de
generalizacin.

TEOREMA 4. Sea
1 2
x x , , ... un conjunto (finito o infinito) de vectores linealmente
independientes en un espacio euclidiano T. Entonces existe un conjunto ortogonal
1 2
e e , , ... en T tal que para cada entero n,
1 n 1 n
e e x x ( ,..., ) ( ,..., ) = . Ms an, el
n
e
puede escogerse de acuerdo con la regla
1 1
e x = , y
n 1 n 1 1 n
e x e e
1 n
...
+ +
= ,
donde
n 1 1 n 1 2 n 1 n
1 1 2 2 n n
x e x e x e

e e e e e e
1 2 n
, , ... ,
+ + +

= = =

.

El mtodo de ortogonalizacin descrito en este teorema, es conocido como el
proceso de ortogonalizacin de Gram-Schmidt.

Ejemplo 13. En este ejemplo se aplica el proceso de ortogonalizacin a los conjuntos
infinitos de vectores, linealmente independientes,
2
1, x, x , ... en el espacio de
polinomios P con producto interior definido por

p q
1
1
p(x)q(x)dx
.
La ortogonalizacin va como sigue:
1
e 1 = ,

1 1

e e
1
1
dx 2
= =
,

2

e
1
1
2
1
x , donde xdx 0
= = =
.
As,
2
e x = ,

2 2

e e
1
2 2
3
1
x dx
= =
,
3
e
2
1 2
x x, = donde

1 1
2 3 3 1 1
1 2
2 3 2
1 1
x dx , x dx 0

= = = =

.
As,
3
e
2 1
3
x = ,

3 3

e e
1
2 2 8 1
3 45
1
(x ) dx
= =
.
Continuando de este modo, se obtiene la sucesin ortogonal

2 1
3
1, x, x ,
3 4 2 3 6 3
5 7 35
x x, x x , ... +
en P. Cuando se presentan con constantes de multiplicacin adecuadas, estos
polinomios se convierten en los famosos polinomios de Legendre del Anlisis.

4.5. EL ESPACIO DE LAS FUNCIONES CONTINUAS POR
TRAMOS

Dada la importancia de las diferentes series que estn por considerarse, reviste
cierto inters expresar la discusin siguiente de tal modo que, de una categora de
funciones, se abarque tanto como sea posible. Desafortunadamente, la mayora de las
generalizaciones que van ms all del espacio de funciones continuas, son sumamente
tcnicas e inaccesibles en un curso de este nivel. Sin embargo, es posible extender las
conclusiones para incluir funciones continuas por tramos, y ya que bien vale la pena
hacer esta modesta generalizacin, esta seccin se destina a demostrar cmo el
conjunto de todas estas funciones en un intervalo determinado, puede estar dentro de
un espacio euclidiano. Por conveniencia, se empieza por recordar la definicin de
continuidad por tramos.

Definicin 10. Se dice que una funcin de valor real f es continua por tramos, en un
intervalo a,b (

si
(i) f es definida y continua en todos, excepto un nmero finito, de puntos de a,b (

y
(ii) los lmites
0 0 0 0
h 0 h 0
f(x ) lm f(x h) , f(x ) lm f(x h)
+
+

= + = existen en cada punto
0
x en
a,b (

. (Observar que slo uno de estos lmites es apropiado, si
0
x es un punto
extremo de a,b (

.)

Cuando
0
x es un punto de continuidad de f cada uno de estos lmites es igual
al valor de f en
0
x y, entonces, se tiene que
0 0 0
f(x ) f(x ) f(x )
+
= = . Con mayor
generalidad, la exigencia de que ambos lmites sean finitos en todas partes de a,b (

,
implica que las nicas discontinuidades de f son discontinuidades por salto. Ms an,
la diferencia
0 0
f(x ) f(x )
+
mide la magnitud del salto de la funcin f en
0
x .
1. Si f es continua por tramos en a,b (

, entonces

b
a
f(x)dx

existe, y es independiente de los valores que f (si los toma) en sus puntos de
discontinuidad. En particular, si f y g son idnticas por todas partes en a,b (

, salvo
en sus puntos de discontinuidad, entonces

b b
a a
f(x)dx g(x)dx =

.
2. Si f y g son continuas por tramos en a,b (

, tambin lo es su producto f.g.
3. Cada funcin continua en a,b (

es continua por tramos.

Dicho esto, ahora se fija la atencin en el problema de convertir el conjunto de
funciones continuas por tramos en a,b (

en un espacio euclidiano. En vista del hecho
de que este conjunto incluye las funciones continuas en a,b (

solo es razonable exigir
que la construccin sea concebida de tal modo que el espacio euclidiano resultante
tenga a C a,b (

como un subespacio. Esto, a su vez, sugiere que se defina a f g por la
frmula

b
a
f g f(x)g(x)dx =
.
Pero la frmula anterior produce realmente un producto interior en el conjunto de
funciones continuas por tramos en a,b (

? La respuesta desagradable es, NO. Para ver
lo que est equivocado, sea n(x) una funcin que es cero dondequiera en a,b (

excepto en un nmero finito de puntos.

Se dice que tal funcin es una funcin nula, y tiene la molesta propiedad de
que

b
2
a
n(x) dx 0 =

a pesar del hecho de que n no es la funcin cero. Esto, por supuesto, invalida el uso de
la frmula cuestionada como un producto interior, ya que, por definicin, el producto
interior de un vector no cero por s mismo, no puede ser cero.

Est perfectamente claro, sin embargo, que la dificultad anterior desaparecer
si se pasa por alto el hecho de que una funcin nula no es idnticamente cero, y se la
trata como si lo fuera. Pero entonces, para que sea consistente, deben mirarse
tambin dos funciones continuas por tramos cualesquiera como idnticas, siempre
que difieran solamente en un nmero finito de puntos, y esto es exactamente lo que se
requiere para hacer que la frmula cuestionada d un producto interior.

4.6. CONVERGENCIA DE SUCESIONES

Definicin 11. Se dice que una sucesin { } { }
k 1 2
x x x , ,... = de vectores, en un espacio
euclidiano V converge al vector x en V si, y solamente si,
k
x x
k
lim 0
= . En este
caso se dice que x es el lmite de { }
k
x , denotndose el hecho de que la sucesin { }
k
x
converge a x, escribiendo { }
k
x x
k
lim
= .

El estudio de convergencia en espacios euclidianos de dimensin finita es
esencialmente el mismo que el estudio de convergencia de sucesiones de nmeros
reales. Sin embargo, en espacios de dimensin infinita, tales como C a,b (

, esta
situacin se vuelve mucho ms compleja y, en consecuencia, ms interesante, porque
entonces el tipo de convergencia definido antes es totalmente distinto del estudiado en
clculo. En efecto se ver momentneamente, que en un espacio de funciones con
producto interior integral, la afirmacin

k

f f
1/ 2
b
2
k
k k
a
l m l m f (x) f(x) dx 0

| |
= = (
|
\

de ninguna manera es la misma que decir que la sucesin { }
k
f converge a la funcin f
en cada punto de a,b (

. En anlisis, tal convergencia se conoce como convergencia
media para enfatizar que se calcula por integracin, lo que, en un sentido, es un
proceso de promedio generalizado.

Ejemplo 14. La sucesin de funciones
{ }
2 3
x, x , x ,... converge en la media en C 1,1 (

a la funcin cero, ya que

1/ 2 1/ 2
1
k 2k
k k k
1
2
l m x 0 l m x dx l m 0
2k 1
| | | |
= = =
| |
+
\ \
.
Aunque
{ }
2 3
x, x , x ,... no converge a cero en cada punto en el intervalo 1,1 (

. En
efecto, en x 1 = la sucesin converge a 1, mientras que en x 1 = no converge (ver
figura 1).


Figura 1. Sucesin de funciones del ejemplo 14

4.7. BASES EN ESPACIOS EUCLIDIANOS DE DIMENSIN
INFINITA

En secciones anteriores se utiliz el proceso de ortogonalizacin de Gram-
Schmidt para demostrar que cada espacio euclidiano de dimensin finita tiene una
base ortonormal
1 n
e e ,..., y que cada vector en tal espacio puede expresarse
unvocamente en la forma

1 1 n n
x x e e x e e ( ) ... ( ) = + + . (1)

Pero el proceso de Gram-Schmidt puede aplicarse igualmente bien en un
espacio euclidiano de dimensin infinita donde puede usarse para dar un conjunto
ortonormal infinito
1 2
e e , ,... . Este hecho sugiere que una versin ampliada de los
resultados iniciales, podra aceptarse, y las secciones restantes del presente tema se
dedican a realizar esta sugerencia. La forma ms obvia para intentar generalizar (1) en
presencia de un conjunto ortonormal
1 2
e e , ,... en un espacio de dimensin infinita V,
es reemplazar su segundo miembro por la serie infinita

k k
x e e
k 1
( )
. (2)

Sin embargo, en la ausencia de cualquier informacin adicional, est claro que
no hay razn a priori para suponer que esta serie converge y, mucho menos, que
converge a x
*
. No obstante, es conveniente tener una notacin que exprese el hecho
de que (2) se deduce de x, aunque de modo puramente formal. La usada con ms
frecuencia es
k k
x x e e
k 1
( )
,
usndose el smbolo (el que, de paso, nada tiene que ver con el usado antes, para
relaciones equivalentes) para enfatizar que las series en cuestin podran no converger
a x. Por supuesto, en caso de que suceda, se escribe
k k
x x e e
k 1
( )
=
=
,
y se dice que la serie converge en la media a x. En uno u otro caso, los productos
interiores
k
x e son llamados las coordenadas o (generalizando) los coeficientes de
Fourier de x con respecto al conjunto ortonormal
1 2
e e , ,...

Est claro que los coeficientes de Fourier de x dependen del conjunto
ortonormal con respecto al cual se calculan. No tan claro, pero igualmente ciero, es
que (1) puede converger en la media a x para un conjunto ortonormal, pero no para
otro. Los ejemplos siguientes ilustran ambos puntos.

Ejemplo 15. Calcule los coeficientes generalizados de Fourier de la funcin f(x) x = ,
x en C , (

con respecto al conjunto ortonormal
senx sen2x
, ,

sen3x
,...

Solucin.
En este caso los coeficientes son

k

x e
1
xsenkxdx, k 1,2,...,
= =

y usando integracin por partes, se obtiene

2
, k 1,3,5,...,
1 2
k
xsenkxdx cosk
k
2
k 2, 4, 6,...
k

= =

En consecuencia, los coeficientes de Fourier de x con respecto al conjunto ortonormal
dado son
k 1
( 1) (2 / k)
, y (2) adopta la forma dada por la expresin

sen2x sen3x
x 2 senx ... .
2 3
| |
+
|
\

En prximos temas se ver que esta serie realmente converge en la media a x, de
suerte que se justifica escribir
k 1
k 1
senkx
x 2 ( 1)
k
=
=
.

Ejemplo 16. Si en el ejemplo anterior se reemplaza al conjunto ortonormal por el
conjunto dado por
1 cos x cos2x
, , ,...,
2

entonces los coeficientes de Fourier son

1 1
xdx 0 y x coskxdx 0, k 1,2,...
2

= = =

.
Pero todas estas integrales son cero, y en consecuencia (2) se convierte en x 0 .
Ahora queda claro que no se tiene una igualdad, ya que

x
1/ 2
1/ 2
3
2
2
x dx 0
3
| |
| |
= = |
|
|
\
\
.

Observacin 2. Para desarrollar f en una serie de funciones ortogonales, es necesario
que no sea ortogonal a cada elemento del conjunto ortogonal. Para evitar este
problema se supondr, en temas siguientes, que un conjunto ortogonal es completo.
Esto quiere decir que la nica funcin ortogonal a cada miembro del conjunto es la
funcin cero.

4.8. DESIGUALDAD DE BESSEL. IGUALDAD DE PARSEVAL.

TEOREMA 5. Sea
1 2
e e , ,... un conjunto ortonormal de vectores en un espacio
euclidiano V, de dimensin infinita, y sea x un vector arbitrario en V. Entonces
k
x e x
2
2
k 1
( )

(desigualdad de Bessel).
Ms an,
1 2
e e , ,... es una base para V si, y solo si,
k
x e x
2
2
k 1
( )
=
=

(igualdad de Parseval).
Demostracin. La demostracin descansa en el clculo del valor de la expresin
k k
x x e e
2
n
k 1
( )
=

para cualquier x en V y cualquier entero n. Ahora
k k k k k k
k k j j k k
x x e e x x e e x x e e
x x x e x e x e e x e e
2
n n n
k 1 k 1 k 1
n n n
k 1 j 1 k 1
( ) ( ) ( )
2 ( )( ) ( ) ( )
= = =
= = =
| | | |
| | =
| |
\ \
| | | |
| | = +
| |
\ \

(Se ha cambiado el ndice de la sumatoria en el primer factor del ltimo trmino, por
conveniencia de clculo). Pero ya que
k
e es ortonormal,
j k
e e
jk
= , se sigue que
j j k k j k j k k
x e e x e e x e x e e e x e
n n n n n
2
j 1 k 1 j 1 k 1 k 1
( ) ( ) ( )( )( ) ( ) .
= = = = =
| |
| |
|
| = =
| |
\
\

Entonces

k k k
x x e e x x e
2
n n
2
2
k 1 k 1
( ) ( )
= =
=

. (3)

Para demostrar la desigualdad de Bessel, es suficiente con observar que
k k k
x x e e x x e
2
n n
2
2
k 1 k 1
0 ( ) ( )
= =
=

para todo x en V y toda n. En consecuencia
k
x e x
n
2
2
k 1
( )
=

para toda n, y las sumas parciales de la serie
k
x e
2
k 1
( )

forman una sucesin no decreciente (acotada) de nmeros reales no negativos. Por un
teorema bien conocido del clculo se concluye que esta serie converge y que
k
x e x
n
2
2
k 1
( )
=

.
Terminando la demostracin ahora se supone que
1 2
e e , ,... es una base ortonormal
para V. Entonces, se sabe que
k k
x x e e
k 1
( )
=
=
.
Por lo tanto
k k
x x e e
2
n
n
k 1
lm ( ) 0
=
=
,
y se sigue de (3) que
k
x x e
n
2
2
n
k 1
lm ( ) 0
=
(
( =
(

.
Por tanto
k
x e x
2
2
k 1
( )
=
=
,
que es la igualdad de Parseval. Finalmente, ya que los eslabones de esta cadena de
razonamientos son reversibles, se concluye que
1 2
e e , ,... es una base siempre que la
igualdad de Parseval se satisfaga, y la demostracin est completa.

4.9. PROBLEMA DE STURM-LIOUVILLE

Ejemplo 17. y'' y 0, y(0) 0, y( ) 0 + = = = .
Solucin.
Para resolver este ejemplo simplemente se aplican las condiciones de frontera a la
solucin general
1 2
c senx c cos x + para deducir que
2
c 0 = y que
1
c es arbitraria. As,
y csenx = donde c es una constante arbitraria, es una solucin general de este
problema particular.
Se dir que las condiciones de frontera dadas, son homogneas siempre que
1 2
0 = = . Las soluciones de un problema con valor en la frontera, que considera un
operador diferencial lineal L : C a,b (

, estn ntimamente ligadas a las soluciones
de la ecuacin
Ly y = , (4)
donde es un parmetro desconocido. En este marco, se requiere encontrar todos los
valores de para los cuales (4) admite soluciones no triviales en , y luego encontrar
las soluciones correspondientes a estas . Se puede observar que (4) puede volver a
escribirse (L I)y 0 = , donde I indica la transformacin identidad que enva cada
funcin en C a,b (

sobre s misma, o como
2 1 0
a (x)y'' a (x)y' a (x) y 0 + + = (

si
2
2 1 0
L a (x)D a (x)D a (x) = + + .

As, para cada valor de , (4) es una ecuacin diferencial lineal homognea de
segundo orden, y el conjunto solucin de cualquier problema (homogneo) con valor
en la frontera que considere esta ecuacin, es el espacio nulo, en , del operador
L I .

Ejemplo 18. Resuelve el problema con valor en la frontera
y'' y 0, y(0) 0, y( ) 0 + = = = . (5)
Solucin.
Aqu, es el espacio de todas las funciones dos veces continuamente diferenciables en
0, (

que se anulan en los puntos extremos del intervalo, y L es el operador lineal de
segundo orden
2
D .(El signo menos se ha introducido nicamente para simplificar los
resultados finales. Sin l, los valores adecuados de seran negativos.) Se distinguen
tres casos, segn que 0, 0 y 0 = < > .
Caso I. 0 = . Aqu la ecuacin diferencial tiene a
1 2
c c x + como su solucin general, y
las condiciones de frontera implican que
1 2
c c 0 = = . As, (5) no tiene soluciones no
triviales cuando 0 = .
Caso 2. 0 < . Ahora,
x x
1 2
y c e c e

= + , y las condiciones de frontera dan otra vez
y 0 .
Caso 3. 0 > . Aqu, la solucin general de y'' y 0 + = es

1 2
y c sen x c cos x = + , (6)
y las condiciones de frontera dan el par de ecuaciones
2
c 0 = ,
1
c sen 0 = .

As, (5) admite soluciones no triviales si, y solamente si, sen 0 = ; esto es,
si, y solamente si, asume uno de los valores
2
n
n = , donde n 1,2,... =

Para cada uno de estos valores de , la constante
1
c en (6) permanece
arbitraria, siguindose que el espacio solucin correspondiente a
n
es el subespacio
unidimensional de C 0, (

generado por la funcin sen(nx). Los nmeros
2
n
n = que
determinan los casos en que (5) tiene soluciones no triviales son llamados los valores
propios para este problema, y cada solucin no trivial correspondiente al valor propio
n
es llamada vector propio o funcin propia perteneciente a
n
. Esta terminologa se
generalizar luego, para incluir una clase mucho ms amplia de problemas y, por el
momento, slo se pide que observe que cualquier conjunto de funciones propias, tal
como sen(x), sen(2x), sen(3x), , uno para cada valor propio, es ortogonal en C 0, (

.

Definicin 12. Un problema de Sturm-Liouville es una ecuacin diferencial ordinaria
de segundo orden homognea de la forma

d
(r(x)y') (q(x) p(x))y 0
dx
+ + = , a x b < < (7)
sujeta a las condiciones de frontera

1 1 2 2
y(a) y'(a) 0, y(b) y'(b) 0 + = + = , (8)
donde p(x), q(x), r(x) y r '(x) son funciones continuas en a,b (

. Los coeficientes en las
ecuaciones (7) y (8) se suponen reales e independientes de . Adems,
1
y
1
no
son cero ambas ni
2
y
2
no son cero ambas.

Resolver un problema de Sturm-Liouville es hallar los valores propios
(autovalores) y funciones propias (autofunciones) no triviales que satisfagan (7) y (8).
Al problema anterior se le suele llamar problema regular o problema normal de
Sturm-Liouville.

TEOREMA 6. (Propiedades del problema regular de Sturm-Liouville).
a. Existe un nmero infinito de valores propios reales que se pueden ordenar en
forma creciente,
1 2 n
... ... < < < < de modo que
n
cuando n .
b. A cada valor propio corresponde solo una funcin propia .
c. Las funciones propias que corresponden a los diversos valores propios son
linealmente independientes.
d. El conjunto de funciones propias que corresponde al conjunto de los valores
propios es ortogonal con respecto a la funcin peso p(x) en el intervalo a,b (

.
Demostracin de la propiedad (d).
Sean
m
y y
n
y las funciones propias que corresponden a los valores propios
m

y
n
, respectivamente. Entonces

'
m m m
d
(r(x)y ) (q(x) p(x))y 0
dx
+ + = (9)

'
n n n
d
(r(x)y ) (q(x) p(x))y 0
dx
+ + = (10)
Al multiplicar la ecuacin (9) por
n
y y la ecuacin (10) por
m
y y restar los
resultados se obtiene
' '
m n m n m n n m
d d
( )p(x)y y y (r(x)y ) y (r(x)y )
dx dx
= .
Se integra por partes este resultado, desde x a = hasta x b = , y se llega a

b
' ' ' '
m n m n m n n m m n n m
a
( ) p(x)y y dx r(b) y (b)y (b) y (b)y (b) r(a) y (a)y (a) y (a)y (a)
( (
=

Ahora bien, las funciones propias
m
y y
n
y deben satisfacer, ambas, las
condiciones en la frontera (8). Por ejemplo, de acuerdo a una de las condiciones se
tiene
' '
1 m 1 m 1 n 1 n
y (a) y (a) 0, y (a) y (a) 0 + = + = . Para que
1
y
1
satisfagan este
sistema, no nulos simultneamente, el determinante de los coeficientes debe ser igual
a cero:
' '
m n n m
y (a)y (a) y (a)y (a) 0 = . Se aplica un argumento semejante a la otra
condicin de (8) y se obtiene
' '
m n n m
y (b)y (b) y (b)y (b) 0 = . Como los dos miembros del
lado derecho de la ecuacin son cero, se ha establecido la relacin de ortogonalidad

b
m n m n
a
p(x)y (x)y (x)dx 0, =
.

Ejemplo 19. Resuelva el problema de valores en la frontera
y'' y 0, y(0) 0, y(1) y'(1) 0 + = = + = . (11)
Solucin.
Se puede comprobar que para 0 < y para 0 = , el problema planteado por la
ecuacin (11) slo posee la solucin trivial y 0 = . Cuando 0 > , la solucin general de
la ecuacin diferencial es
1 2
y c cos x c sen x = + . Ahora bien, la condicin y(0) 0 =
implica que en esta solucin,
1
c 0 = . Cuando
2
y c sen x = , se satisface la segunda
condicin de frontera y(1) y'(1) 0 + = , siempre que
2 2
c sen c cos 0 + = . Al
escoger
2
c 0 , se ve que esta ltima ecuacin es equivalente a tg = . Si se
hace x = , los valores propios del problema (11) son, entonces,
2
n n
x = , donde

n
x , n 1,2,3,... = son las races positivas consecutivas. En general, las funciones propias
del problema son
{ }

n
sen( x , n 1,2,3,... =

Al identificar
1 1
r(x) 1, q(x) 0, p(x) 1, 1, 0, = = = = =
2 2
1, 1 = = se ve que
la ecuacin (11) es un problema regular de Sturm-Liouville. As,
{ } n
sen( x ,
n 1,2,3,... = es un conjunto ortogonal con respecto a la funcin peso p(x) 1 = en el
intervalo 0,1 (

.

A continuacin se consideran cierto nmero de problemas de Sturm-Liouville
que se encontrarn repetidamente en trabajos posteriores y que, por conveniencia de
referencia, se resuelven aqu, de una vez por todas. En cada caso se limitar a
encontrar valores propios y funciones propias, y se pospone cualquier discusin de
cmo usar esta informacin para resolver problemas especficos con valor en la
frontera que involucren el operador considerado.

Ejemplo 20. Resuelva el sistema de Sturm-Liouville
y'' y 0, y(0) 0, y(L) 0 + = = = , 0 x L < <
Solucin.
Esta es una variante de un problema que se ha considerado muchas veces y, en esta
ocasin, se encuentra que las constantes dadas por

2 2
n
2
n
, n 1,2,...
L
= =
y las funciones

n
n x
y (x) sen , n 1,2,...,
L
| |
= =
|
\

son considerados conjuntos completos de valores propios y funciones propias. Se ver
en detalle su resolucin: Se tienen tres opciones:
2
0 k < = :
2 2 2
y'' k y 0 m k 0 m k = = =
kx kx
1 2
y(x) c e c e
= + ,
1 2 2 1
y(0) c c 0 c c = + = = .
Sustituyendo:
kx kx
1
y(x) c (e e )
= ,
kL kL
1 1 2
y(L) c (e e ) 0 c 0, c 0 y(x) 0
= = = = =
Conduce a solucin trivial.
0 : =
1 2 2 1
y'' 0 y(x) c x c y(0) c 0 y(x) c x = = + = = =
1
y(L) c L 0 = = , pero L 0 , entonces conduce a solucin trivial.
2
0 : k > = :
2
1 2
y'' k y 0 y(x) c cos(kx) c sen(kx) + = = +
Al evaluar
1 2
y(0) c 0 y(x) c sen(kx) = = =
2
n
y(L) c sen(kL) 0 sen(kL) 0 kL n k
L
= = = = = .
Al sustituir se tiene:
n 2
n
y c sen x
L
| |
=
|
\

(Solucin no trivial).
Como
2
n
n n
k
L L
| |
= =
|
\
Autovalores.

Ejemplo 21. Resuelva el sistema de Sturm-Liouville y'' y 0, y'(0) 0, y'(L) 0 + = = = .
Solucin.
Un clculo similar en todos los aspectos al que se us para resolver el ejemplo
anterior, revela que los valores propios para este problema son las constantes no
negativas

2 2
n
2
n
, n 0,1,2,...
L
= =
y que

n
n x
y (x) cos , n 0,1,2,...,
L
| |
= =
|
\

constituye un conjunto completo de funciones propias.

Ejemplo 22. Resuelva el sistema de Sturm-Liouville
y'' y 0, y(0) 0, hy(L) y'(L) 0, + = = + =
considerando que h y L son constantes positivas.
Solucin.
Como es usual, se estudia por casos, dependiendo del signo algebraico de , y
nuevamente se encuentra que no hay valores propios 0 . Por otra parte, cuando
0 > ,
1 2
y c sen x c cos x, = + y la primera condicin de frontera implica que
2
c 0 = .
As, queda pendiente (si es que queda) elegir en forma tal que la funcin
1
y c sen( x) = , con
1
c 0 , satisfaga la ecuacin hy(L) y'(L) 0 + = . Esto, a su vez,
implica que debe escogerse en forma tal que

L
sen( L) cos( L)
hL
= ; (12)
una ecuacin que vuelve a escribirse como
t g
hL
= (13)
haciendo L = . Aunque es imposible resolver (12) explcitamente para , sus
soluciones pueden visualizarse conforme surgen, a travs de (13), en los puntos de
interseccin de las grficas de las funciones tg y / hL.

Hay una infinidad de tales puntos localizados simtricamente con respecto al
origen. As, el problema considerado tiene un nmero infinito de valores propios
positivos
2
n
n
2
, n 1,2,...
L
= = . A continuacin se muestra la grfica (figura 2) para el

caso h 1, L 1, x = = = ):


Figura 2. Problema del Sturm-Liouville del ejemplo 22

con
n n 1 +
< para toda n, y
n
n
lm
= . (Tambin observe, de la geometra del caso,

se tiene que
n 1 n
n
lm( )
+
= ). Las funciones

n
n
x
y (x) sen sen( x), n 1,2,...
L
| |
= = =
|
\

constituyen un conjunto completo de funciones propias para este problema.

Ejemplo 23. Considere el problema con valores en la frontera

y'' 2y' (1 )y 0
y(0) 0, y(1) 0
+ + =
= =
.
Halle los valores propios y las funciones propias y determine una relacin de
ortogonalidad entre las soluciones.
Solucin.
Si 0 < se tiene que la solucin es
x
1 2
y(x) e c cos(kx) c sen(kx)
= + (

, donde
2
k = .
Con las condiciones iniciales se tiene
1
1 2
0 c , 0 e c sen(k)
= = . Si
2
c 0, k n = . Por lo
tanto
2 2
n = y
x
n
y (x) e sen(n x)
= . Sea la funcin de peso

2x
p(x) e = , se tiene la
siguiente relacin de ortogonalidad entre las soluciones:

1
2x x x
0
e e sen(m x)e sen(n x)dx 0

=
.

4.10. FORMA AUTOADJUNTA DE PROBLEMAS DE
STURM-LIOUVILLE

Si los coeficientes son continuos y a(x) 0 para todo x en algn intervalo,
entonces cualquier ecuacin diferencial de segundo orden
a(x)y'' b(x)y' (c(x) d(x))y 0 + + + = (14)
se puede llevar a la forma autoadjunta, ecuacin (7), si se multiplica por el factor
integrante
(b(x) a(x))dx 1
e
a(x)

.
Para comprobarlo, se observa que la ecuacin diferencial
(b / a)dx (b / a)dx (b / a)dx (b / a)dx b(x) c(x) d(x)
e y'' e y' e e y 0
a(x) a(x) a(x)
| |

+ + + =
|
\

es la misma que

(b / a)dx (b / a)dx (b / a)dx d c(x) d(x)
e y' e e y 0
dx a(x) a(x)
| |
(
+ + =
|
(

\
. (15)

Entonces
(b / a)dx
r(x) e
= ,
(b / a)dx c(x)
q(x) e
a(x)

= ,
(b / a)dx d(x)
p(x) e
a(x)

= .

Desde luego que no es necesario llevar a una ecuacin diferencial de segundo
orden como la (14) a la forma autoadjunta (7), para resolverla. Se usar la forma de
la ecuacin (15) para determinar la funcin peso p(x) que se necesita en la relacin de
ortogonalidad.

Ejemplo 24. Encuentre los valores propios y las funciones propias para el problema de
Sturm-Liouville y'' 4y' (4 9 )y 0, y(0) 0, y(a) 0 + + = = = .
Solucin.
En primer lugar, se escribir en forma autoadjunta como

4x 4x 4x
d dy
e (4e (9e ))y 0, y(0) 0, y(a) 0
dx dx
| |
+ = = =
|
\
.
As, las funciones propias correspondientes a diferentes valores propios para
este problema sern mutuamente ortogonales en el espacio euclidiano C 0, a (

con
producto interior calculado con relacin a la funcin peso
4x
p(x) 9e = .
Caso I. 0 > . Aqu la solucin general de y'' 4y' (4 9 )y 0 + + = viene dada por
( 2 3 )x ( 2 3 )x
1 2
y c e c e
+
= + ,
y las condiciones de frontera implican que
1 2
c c 0, + =
( 2 3 )a ( 2 3 )a
1 2
c e c e 0
+
+ = .
As
1 2
c c 0 = = y el problema no tiene valores propios positivos.
Caso II. 0 = . Ahora la solucin general de la ecuacin es
2a
1 2
(c c x)e ,
+ y
nuevamente se encuentra que
1 2
c c 0 = = .
Caso III. 0 < . Aqu
2x
1 2
y e (c sen(3 x) c cos(3 x))
= +
y la exigencia de que y(0) y(a) 0 = = da
2 1
c 0, c sen a 0 = = . En consecuencia el
problema tiene soluciones no triviales si, y solamente si, satisface la ecuacin
sen(3 a) 0 = , siguindose que los valores propios para este problema son

2 2
n
2
n
, n 1,2,...
9a
= =
Para encontrar un conjunto de funciones propias correspondientes, ahora se hace

1 2
c 1, c 0 = = y
n
= y se obtienen

2x
n
n x
(x) e sen , n 1,2,...
a
| |
= =
|
\
.

Ejemplo 25. Encuentre los valores propios y las funciones propias del problema de
Sturm-Liouville dado por y'' ky 0, y(0) y(2 ), y'(0) y'(2 ) + = = = .
Solucin.
Caso 1. k 0 = .
1 2
y''(x) 0 y(x) c c x = = + .
Autovalor k 0 = . Autofuncin
1 1
y c c 0 =

1 1 2 2 1 1
y(0) y(2 ) c c 2 c c 0 y(x) c , y'(0) y'(2 ) 0 0 y(x) c = = + = = = = = .
Caso 2.
2
k , 0 = .
2
3 4
y'' y 0 y(x) c cosh( x) c senh( x) = = + .
= = +
= = +

3 3 4
4 3 4
y(0) y(2 ) c c cosh(2 ) c senh(2 ) ,
y'(0) y'(2 ) c c senh(2 ) c cosh(2 )

3 4 3 4
c (cosh(2 ) 1) c senh(2 ) 0 , c senh(2 ) c (cosh(2 ) 1) 0 + = + = . Por lo tanto se
tiene un sistema de ecuaciones lineales homogneo. El valor del determinante del
sistema corresponde a
= +
=
2 2 2 2
(cosh(2 ) 1) sen h(2 ) cos h(2 ) 2cosh(2 ) 1 sen h(2 )
2 2cosh(2 ) 0

En consecuencia
3 4
c c 0 = = es la nica solucin y se tiene y(x) 0 = .
Caso 2 (Alternativo).
2
k , 0 = .
2 x x
3 4
y'' y 0 y(x) c e c e

= = + .

2 2 2 2
3 4 3 4 3 4 3 4
y(0) y(2 ) c c c e c e , y'(0) y'(2 ) (c c ) (c e c e )

= + = + = =

2 2 2 2
3 4 3 4
c (e 1) c (e 1) 0 , c (e 1) c (e 1) 0

+ = = . Por lo tanto se tiene un
sistema de ecuaciones lineales homogneo. El valor del determinante del sistema
corresponde a
2 2 2 2 2 2
2(e 1)(e 1) 2(1 e e 1) 2(2 (e e )) 0

= + = + .
En consecuencia
3 4
c c 0 = = es la nica solucin y se tiene y(x) 0 = .
Caso 3.
2
k , 0 = .
2
5 6
y'' y 0 y(x) c cos( x) c sen( x) + = = + .
= = + =
=

5 5 6
6 6 5
y(0) y(2 ) c c cos(2 ) c sen(2 ) , y'(0) y'(2 )
c (c cos(2 ) c sen(2 ))

5 6 5 6
c (cos(2 ) 1) c sen(2 ) 0 , c sen(2 ) c (1 cos(2 )) 0 + = + = . Por lo tanto se tiene
un sistema de ecuaciones lineales homogneo. El valor del determinante del sistema
corresponde a
2 2 2 2
(cos(2 ) 1) sen (2 ) cos (2 ) 2cos(2 ) 1 sen (2 ) 2(cos(2 ) 1) = + = .
cos(2 ) 1 0 cos(2 ) 1 2 2n n = = = = . Por tanto los autovalores vienen
dados por
2 2
n
k n n 1,... = = = y las autofunciones
n 5 6
y (x) c cos(nx) c sen(nx) = + .
SOLUCIN FINAL:
2
k n n 0,... = = y las autofunciones
n 5 6
y (x) c cos(nx) c sen(nx) = + .
4.11. FUNCIONES PARES Y FUNCIONES IMPARES

La tarea de calcular las integrales que surgen en el estudio de series
ortogonales, puede simplificarse a menudo, al sacar provecho de la simetra de las
funciones involucradas. Esta tcnica, generalmente se formaliza por la introduccin de
las nociones de funciones pares e impares, como sigue:

Definicin 13. Se dice que una funcin f, definida en un intervalo centrado en el
origen, es par si f( x) f(x) = para toda x en el dominio de f, e impar si f( x) f(x) = .

Esto, por supuesto, es justamente otra forma de decir que una funcin es par si
su grfica es simtrica con respecto al eje vertical, e impar si su grfica es simtrica
con respecto al origen. As, para valores enteros de n,
n
x es par, si n es par; e impar,
si n es impar. La importancia de las funciones pares e impares para este trabajo, nace
de las igualdades

a a
a 0
f(x)dx 2 f(x)dx
=

siempre que f sea par e integrable, y

a
a
f(x)dx 0

siempre que f sea impar e integrable. Una propiedad un tanto menos obvia de las
funciones pares e impares, se establece en el siguiente

TEOREMA 7. Cada funcin en el intervalo a, a (

puede expresarse justamente en
forma de una suma de una funcin par y una funcin impar.
Demostracin. Sea f una funcin arbitraria en a, a (

, y sea

E 0
f(x) f( x) f(x) f( x)
f (x) , f (x)
2 2
+
= = .

Es trivial verificar que
E
f es par,
0
f impar, y que
E 0
f f f = + . As, f tiene cuando
menos una descomposicin de la forma deseada, y queda por demostrar que es la
nica. Con este fin, se supone que tambin se tiene
E 0
f g g = + , con
E
g par y
0
g impar.
Entonces,
E 0 E 0
f f g g + = + , y
E E 0 0
f g g f = . Pero la diferencia de dos funciones pares
es par, y la diferencia de dos funciones impares es impar. As, la funcin definida por la
igualdad anterior es simultneamente par e impar, y as, debe ser la funcin cero. En
otras palabras,
E E 0 0
f g g f 0 = = , de donde se sigue que
E E
g f = y
0 0
g f = , como se
deseaba. Las funciones
E
f y
0
f son conocidas, respectivamente, como las partes par e
impar de f.

Ejemplo 26. Si
x
f(x) e = , entonces
x x
E
e e
f (x)
2

+
= y
x x
0
e e
f (x)
2

= .
As, las partes par e impar de la funcin exponencial son el coseno hiperblico y el
seno hiperblico, respectivamente.

4.12. SERIES DE FOURIER

En esta seccin se inicia el estudio de series ortogonales considerando el
desarrollo en serie relativo a las funciones 1, cos x, senx, cos(2x), sen(2x),... . Ya se
ha visto que estas funciones son mutuamente ortogonales en PC , (

, y se probar
en forma breve que tambin son una base. Garantizando lo cierto de este hecho,
puede utilizarse la frmula
k
k
k
x e
x e
e
2
k 1

para expresar cualquier funcin continua por tramos f en el intervalo , (

en la
forma
2 2 2
k 1
f 1 f cos(kx) f sen(kx)
f(x) cos(kx) sen(kx)
1 cos(kx) sen(kx)
=
(

(
= + +
(

(media)
donde la notacin (media) indica que la serie en cuestin converge en la media a f.
Pero ya que

2
1 dx 2
= =
,

2
2
cos(kx) cos (kx)dx
= =
,

2
2
sen(kx) sen (kx)dx
= =
,
entonces se tiene que

0
k k
k 1
a
f(x) (a cos(kx) b sen(kx)
2
=
= + +
(media) , (16)
donde

k
1
a f(x)cos(kx)dx , k 0,1,2,...
= =
,

k
1
b f(x)sen(kx)dx , k 1,2,...
= =

Esta representacin particular es conocida como su desarrollo en serie de
Fourier en el intervalo , (

, y las
k
a y
k
b son llamadas los coeficientes de
Fourier. Una vez ms se enfatiza que (16) debe leerse como una afirmacin de que la
serie en cuestin converge en la media a f, no que converja punto por punto en el
sentido de que
0
0 k 0 k 0
k 1
a
f(x ) (a cos(kx ) b sen(kx ))
2
=
= + +

para toda
0
x en , (

. En efecto ya que el valor de f en
0
x puede cambiarse
arbitrariamente, sin cambiar los valores de sus coeficientes de Fourier, esto sera
demasiado esperar. Pero, sorprendentemente, siempre que f se comporte
razonablemente bien, converge a f(x) para toda x.

Ejemplo 27. Encuentre el desarrollo en serie de Fourier, de la funcin
1 x 0
f(x)
1 0 x
< <
=

< <
.
Solucin.
En este caso f es una funcin impar en , (

. En consecuencia, as es f(x)cos(kx), de
modo que
k
a 0 = , para toda k. Por otra parte, f(x)sen(kx) es par y, por lo tanto, se
tiene
k
0
4
k 1,3,5,... 2 2
b sen(kx)dx (1 cos(k )) k
k
0 k 2, 4, 6,...
= = =

.
Por tanto, el desarrollo en serie de Fourier de f es
k 1
4 sen(3x) sen(5x) 4 sen(2k 1)x
f(x) sen(x) ...
3 5 2k 1
=
(
= + + + =
(

A continuacin se presentan algunas grficas (ver figuras 3, 4 y 5):
4sen(x)
y =

Figura 3. Primer desarrollo de Fourier para el ejemplo 27

4 sen(3x) sen(5x)
y sen(x)
3 5
(
= + +
(

Figura 4. Segundo desarrollo de Fourier para el ejemplo 27
4 sen(3x) sen(5x) sen(7x) sen(9x)
y sen(x)
3 5 7 9
(
= + + + +
(

Figura 5. Tercer desarrollo de Fourier para el ejemplo 27

TEOREMA 8. Sea f una funcin suave por tramos en PC , (

, por lo que se entiende
que f tiene una primera derivada continua por tramos en , (

. Entonces el desarrollo
en serie de Fourier para f converge punto por punto dondequiera en , (

, y tiene el
valor
0 0
f(x ) f(x )
2
+
+

en cada punto
0
x en el interior del intervalo, y
f( ) f( )
2
+
+

en .

Cuando se aplican estos resultados al ejemplo anterior, permiten afirmar que la
serie
4 sen(3x) sen(5x)
sen(x) ...
3 5
(
+ + +
(

converge punto por punto en el intervalo , (

a

1 si x 0,
0 si x , 0,
1 si 0 x
< <
=
< <
.

As, por ejemplo, cuando x 2 = , el valor de la serie es 1, y entonces
4 1 1 1
1 1 ...
3 5 7
(
= + +
(

,
o
1 1 1
1 ...
4 3 5 7
= + +

TEOREMA 9. El desarrollo en serie de Fourier de una funcin f suave por tramos en
PC , (

converge punto por punto en toda la recta real. Ms an, si F indica la
extensin peridica de f, entonces el valor de la serie es
0
F(x ) cuando
0
x es un punto
de continuidad de F, y
0 0
F(x ) F(x )
2
+
+

cuando
0
x es una discontinuidad por salto de F.

Ejemplo 28. Halle el desarrollo en serie de Fourier de la funcin f(x) x , x = < < .
Solucin.
En este caso, f es una funcin par en , (

. En consecuencia,
k
b 0 = para toda k,
mientras que, para k 0 ,

k
0 0
0
2
2 2
0
2 2 xsen(kx) 1
a x cos(kx)dx sen(kx)dx
k k
4
k 1,3,5,...,
2 2
cos(kx) (cos(k ) 1)
k
k k
0 k 2, 4, 6,...
(
= = (

(

= = =

=

Finalmente, cuando k 0 = , se tiene

0
0
2
a xdx
= =
,
siguiendo que
2 2 2
k 1
4 cos(3x) cos(5x) 4 cos(2k 1)x
x cos(x) ...
2 2
3 5 (2k 1)
=
| |
= + + + =
|

\

Algunas grficas (ver figuras 6 y 7):
4 cos(3x)
y cos(x)
2 9
| |
= +
|
\

4 cos(3x) cos(5x) cos(7x)
y cos(x)
2 9 25 49
| |
= + + +
|
\


Ejemplo 29. Sea g la funcin en PC , (

definida por
0 x 0
g(x)
1 0 x
< <
=

< <
.
Entonces, con f como en el ejemplo 27,
1 1 1
g (1 f) f,
2 2 2
= + = +
y se concluye que el desarrollo en serie de Fourier de g es
k 1
1 2 sen(2k 1)x
f(x)
2 2k 1
= +

.
La moraleja de este ejemplo es que una serie de Fourier puede, algunas veces,
encontrarse sin recurrir a la integracin. Otra forma de expresar la serie sera
k
k 1
1 1 1 ( 1)
f(x) sen(kx)
2 k
=

= +

(ver figuras 8, 9, 10 y 11) . Grficas:
1
1 2
k 1 s (x) sen(x)
2
= = +

3
1 2 sen(3x)
k 3 s (x) sen(x)
2 3
(
= = + +
(


7
1 2 sen(3x) sen(5x) sen(7x)
k 7 s (x) sen(x)
2 3 5 7
(
= = + + + +
(


11
1 2 sen(3x) sen(5x) sen(7x) sen(9x) sen(11x)
s (x) sen(x)
2 3 5 7 9 11
(
= + + + + + +
(


Figura 11. Cuarto desarrollo de Fourier para el ejemplo 29

4.13. SERIES DE SENOS Y COSENOS

En los ejemplos anteriores se abus del hecho de que las funciones
consideradas eran pares o impares, para simplificar la tarea de encontrar su desarrollo
en serie de Fourier. Esta tcnica puede ser explotada ms frecuentemente de lo que se
pudiera esperar, y es de suma importancia para sacarla al descubierto.

Especficamente, si f es una funcin par en PC , (

, entonces, para todos los
valores de k, f(x)cos(kx) es par, y f(x)sen(kx) es impar. As, se tiene que

0
f(x)cos(kx)dx 2 f(x)cos(kx)dx,

=

f(x)sen(kx)dx 0,

siguindose que el desarrollo en serie de Fourier de una funcin par en PC , (

incluye solamente trminos cosenos y pueden calcularse de acuerdo con la frmula
0
k
k 1
a
f(x) a cos(kx)
2
=
= +
,
donde

k
0
2
a f(x)cos(kx)dx
.
Un argumento similar muestra que el desarrollo en serie de Fourier de una funcin par
en PC , (

incluye solamente trminos senos, y se calcula con
k
k 1
f(x) b sen(kx)
=
=
,
donde

k
0
2
b f(x)sen(kx)dx
.

En las aplicaciones de la teora de las series de Fourier, frecuentemente se
necesita obtener el desarrollo en serie para una funcin f continua por tramos, que
est definida solamente en el intervalo 0, (

. Una forma de hacerlo es, extendiendo f
en todo el intervalo 0, (

(donde con esto quiere decirse que una F ests definida en
, (

de tal modo, que F coincide con f en 0, (

), y entonces F se desarrolla como
una serie de Fourier.

Lo esencial de este mtodo tiene relacin con la forma en que f es extendida a
, (

. Esto por supuesto, puede hacerse en la forma que sea (en tanto que la funcin
resultante pertenezca a PC , (

), pero las dos extensiones son las ms convenientes
e importantes. La primera es la llamada extensin par de f, indicada por
f
E , y definida
por
f
f(x) 0 x
E (x)
f( x) x 0

=

<

mientras la segunda es la extensin impar de f, indicada por
f
O , y definida por
f
f(x) 0 x
O (x)
f( x) x 0

=

<
.
Se tienen entonces los desarrollos

0
f k k
0
k 1
a 2
E (x) a cos(kx), a f(x)cos(kx)dx
2
=
= + =

y

f k k
0
k 1
2
O (x) b sen(kx), b f(x)sen(kx)dx
=
= =

.

Estas series son llamadas respectivamente, los desarrollos en serie de
Fourier de coseno y de seno de f. Un trmino, un tanto extraviado desarrollo semi-
intervalo tambin es usado en este contexto.

Ejemplo 30. Halle el desarrollo de Fourier de
2
0 x 0
f(x)
x 0 x
<
=

<

y pruebe, usando el desarrollo anterior, que
n 1 2
2
n 1
( 1)
12
n
+
=

=
.
Solucin.
Sea
0
n n
n 1
a
f(x) (a cos(nx) b sen(nx))
2
=
= + +
.

3 2
2
0
0
0
1 1 x
a x dx . .
3 3
= = =


2
2
n
0 0
0
2
2
0
0
2 2
1 1 x sen(nx) 2
a x cos(nx)dx x.sen(nx)dx
n n
1 x sen(nx) 2 x.cos(nx) sen(nx)
n n n
n
2 2 2 2
n impar : . n par : . Por tanto : a
n n n n
n n

(
(
= =
(

(
(
(
= +
(
(

( (
= = = =
( (

n
n
2
2
( 1) .
n
=

2
2
n
2
0
0
0
2 2 2
2
2 2
n 1 n
2
n
1 1 x cos(nx) 2 x 1
b x sen(nx)dx sen(nx) cos(nx)
n n n
n
1 2 2 1 4 1
n impar : n par :
n n n n n
n n
( 1) 2(1 ( 1)
Por tanto : b
n
+
(
(
(
= = + +
(
(

( (
( (
= + = = =
( (
( (

( (

=
n 1 n
2 3
) ( 1) 2(( 1) 1)
n
n n
+
(

= +
(
(

De modo que:
2 n 1 n
n
2 3
n 1
2 ( 1) 2(( 1) 1)
f(x) ( 1) cos(nx) sen(nx)
6 n
n n
+
=
( | |

= + + + ( |
|
(
\
.
Si se hace x 0 = se tiene:
2 2 n 1 2
n n
2 2 2
n 1 n 1 n 1
2 2 ( 1)
0 ( 1) ( 1) .
6 6 12
n n n

+
= = =

= + = =

Ejemplo 31. Encuentre la serie de Fourier de
0 x 0
f(x)
1 0 x
< <
=

<
.

Rta.
n
n 1
1 1 1 ( 1)
f(x) sen(nx)
2 n
=

= +
(ver figuras 12, 13, 14 y 15)

1
1 2
n 1 s (x) sen(x)
2
= = +

3
1 2 sen(3x)
n 3 s (x) sen(x)
2 3
(
= = + +
(


7
1 2 sen(3x) sen(5x) sen(7x)
n 7 s (x) sen(x)
2 3 5 7
(
= = + + + +
(


11
1 2 sen(3x) sen(5x) sen(7x) sen(9x) sen(11x)
s (x) sen(x)
2 3 5 7 9 11
(
= + + + + + +
(


Figura 15. Cuarto desarrollo de Fourier para el ejemplo 31

Ejemplo 32. Encuentre la serie de Fourier de
0 x 0
f(x)
x 0 x
< <
=

<

y usela para calcular el valor de la serie
k
2
k 1
1 ( 1)
(2k 1)
(2k 1)
=
| |
+ |
|
.
Solucin.

2 2
0
0
0
1 x
a xdx .
2 2 2

= = = =

k
k
2
0
1 ( 1) 1
a x cos(kx)dx
k

= =
.

k 1
k
0
1 ( 1)
b xsen(kx)dx
k
+
= =
.
Se tiene entonces que
k k 1
2
k 1
( 1) 1 ( 1)
f(x) cos(kx) sen(kx)
4 k
k
+
=
| |

= + + |
|

Si x 0 : =
2
2 2
k 1 k 1
2 1 1
4 8
(2k 1) (2k 1)

= =

= =

Si x :
2
=
k 1
2
k 1 k 1
k 1 k 1 k
k 1 k 1 k 1
2 1 ( 1)
cos (2k 1) sen k
4 2 k 2
(2k 1)
( 1) ( 1) ( 1)
sen k
4 k 2 4 2k 1 4 2k 1

+
= =

+ +
= = =
| | | |
= +
| |
\ \
| |
= = =
|

\

Por lo tanto:
k 2
2
k 1
1 ( 1)
(2k 1) 8 4
(2k 1)
=
| |

+ = |
|
.

Ejemplo 33. Desarrolle
0, x 0
f(x)
x, 0 x
< <
=

< <

en serie de Fourier.
Solucin.
Con p = , se tiene que

2
0
0
0
0
1 1 1 x
a f(x)dx 0dx ( x)dx x ,
2 2

(
(
= = + = =
(
(

(

0
n
0
n
2 2
0
0 0
1 1
a f(x)cos(nx)dx odx ( x)cos(nx)dx
1 sen(nx) 1 1 cos(nx) cos(n ) 1 1 ( 1)
( x) sen(nx)dx
n n n n
n n

(
= = +
(

(
+
= + ( = = =

(

En forma anloga, se obtiene que

n
0
1 1
b ( x)sen(nx)dx .
n
= =

Por lo tanto,
n
2
n 1
1 ( 1) 1
f(x) cos(nx) sen(nx)
4 n
n
=
| |

= + + |
|

4.14. CAMBIO DE INTERVALO

Hasta ahora se ha tratado exclusivamente con funciones en los intervalos , (

y 0, (

. Para muchas finalidades, sin embargo, este marco es demasiado restrictivo, y
ahora se propone generalizar los resultados a un intervalo arbitrario a,b (

. Pero ms
bien empezar desde luego con el caso ms general, resultar ms sencillo si primero
se consideran intervalos de la forma p,p (

y sus espacios euclidianos asociados
PC p,p (

. Porque as, la situacin puede resolverse mejor. En realidad, casi es obvio
que las funciones

x x 2 x 2 x
1, cos , sen , cos , sen , ...
p p p p
| | | | | | | |
| | | |
\ \ \ \

son mutuamente ortogonales en PC p,p (

. Ms an, justamente como en el caso
donde p = , puede demostrarse que estas funciones son una base para este espacio
y, en consecuencia, que sus series ortogonales asociadas (las que, de paso, son an
llamadas series de Fourier) convergen en la media. Y finalmente, teniendo en cuenta la
longitud del intervalo, todas las observaciones anteriores referentes a la convergencia
por puntos, son vlidas en este caso. Entonces
0
k k
k 1
a k x k x
f(x) a cos b sen
2 p p
=
| | | | | |
= + +
| | |
\ \ \
,
donde

p
k
p
1 k x
a f(x)cos dx
p p
| |
=
|
\
,

p
k
p
1 k x
b f(x)sen dx
p p
| |
=
|
\

para toda k. Y con esto se ha logrado.

La discusin anterior puede adaptarse fcilmente para manipular el espacio
euclidiano PC a,b (

. En efecto, si se establece 2p b a = , de manera que
a,b a, a 2p = + ( (

, las funciones anteriores tambin son una base para PC a, a 2p + (

.
Esto inmediatamente lleva a las frmulas para calcular el desarrollo en serie de Fourier
de una funcin f en PC a,b (

de la forma:
0
k k
k 1
a 2k x 2k x
f(x) a cos b sen
2 b a b a
=
| | | | | |
= + +
| | |

\ \ \
,
donde

b
k
a
2 2k x
a f(x)cos dx,
b a b a
| |
=
|

\

b
k
a
2 2k x
b f(x)sen dx,
b a b a
| |
=
|

\

para toda k.

Ejemplo 35. Halle el desarrollo en serie de Fourier, en PC 0,1 (

de la funcin f(x) x = .
Solucin.
En este caso, b a 1 = , de modo que

1 1
k k
0 0
a 2 x cos(2k x)dx, b 2 xsen(2k x)dx = =

.
Integrando por partes se tiene ahora:
0 k k
1
a 1 , a 0, k 0 , b
k
= = =
.
Por lo tanto
1 1 sen(4 x) sen(6 x)
f(x) sen(2 x) ...
2 2 3
| |
= + + +
|
\
.

Ejemplo 36. Encuentre el desarrollo en serie de Fourier de la funcin
x 2 2 x 3
f(x)
4 x 3 x 4

=

.
Solucin. Las frmulas seran

4 3 4
k
2 2 3
a f(x)cos(k x)dx (x 2)cos(k x)dx (4 x)cos(k x)dx = = +

4 3 4
k
2 2 3
b f(x)sen(k x)dx (x 2)sen(k x)dx (4 x)sen(k x)dx = = +

Aunque estas integrales pueden calcularse directamente, los clculos pueden
simplificarse considerablemente al tomar en cuenta el siguiente argumento:
Designando por F la extensin peridica de f en todo el eje x. Entonces, las funciones
F(x)cos(k x) y F(x)sen(k x) son peridicas con perodo 2, y se tiene

a 2 4 a 2 4
a 2 a 2
F(x)cos(k x)dx f(x)cos(k x)dx, F(x)sen(k x)dx f(x)sen(k x)dx
+ +
= =

para cualquier nmero real a. (En este punto se utiliza el hecho obvio de que si g es
continua por tramos en ( , ) con perodo 2p, entonces

a 2p b 2p
a b
g(x)dx g(x)dx
+ +
=

para cualquier par de nmeros reales a y b.) Ahora se hace a 1 = y se obtiene

1 1
k k
1 1
a F(x)cos(k x)dx , b F(x)sen(k x)dx

= =

Pero en el intervalo 1,1 (

, F coincide con la funcin par x . Entonces
k
b 0 =
para toda k, y

1
k
0
a 2 x cos(k x)dx =
.
As

2 2
0 k
4
k impar
a 1, a
k
0 k par, k 0
= =

,
y el desarrollo en serie de Fourier de f es
2 2 2
1 4 cos(3 x) cos(5 x)
f(x) cos( x) ...
2
3 5
| |
= + + +
|
\
.

Ejemplo 37. Encuentre la serie de Fourier de cosenos de
1 0 x 1
f(x)
2 x 1 x 2

=

.
Solucin.
x 2 2 x 1
f(x) 1 1 x 1
2 x 1 x 2
+
.

2
2
1 2
1
0
0
0 1
1
x 3 3
a dx (2 x)dx x 2x 1 2 .
2 2 2
| |
= + = + = + = |
|
\

1 2 2
n
0 1 1
n n n
a cos x dx 2 cos x dx x cos x dx
2 2 2
| | | | | |
= +
| | |
\ \ \

n 1,5, 9,13,... =

1
1
0
0
n 2 n 2
cos x dx sen x .
2 n 2 n
| | | |
= =
| |

\ \

n 3,7,11,15,... =

1
1
0
0
n 2 n 2
cos x dx sen x .
2 n 2 n
| | | |
= =
| |

\ \

n 2, 6,10,14,... =

1
1
0
0
n 2 n
cos x dx sen x 0.
2 n 2
| | | |
= =
| |
\ \

n 4, 8,12,16,... =

1
1
0
0
n 2 n
cos x dx sen x 0.
2 n 2
| | | |
= =
| |
\ \

n 1,5, 9,13,... =

2
2
n
1
1
n 4 n 4
a 2 cos x dx sen x
2 n 2 n
| | | |
= = =
| |

\ \
.
n 3,7,11,15,... =

2
2
n
1
1
n 4 n 4
a 2 cos x dx sen x
2 n 2 n
| | | |
= = =
| |

\ \
.
n 2, 6,10,14,... =

2
2
n
1
1
n 4 n
a 2 cos x dx sen x 0
2 n 2
| | | |
= = =
| |
\ \
.
n 4, 8,12,16,... =

2
2
n
1
1
n 4 n
a 2 cos x dx sen x 0
2 n 2
| | | |
= = =
| |
\ \
.
n 1,5, 9,13,... =

2
2
2 2 2 2
1
1
n 2x n 4 n 4 2
x cos x dx sen x cos x
2 n 2 2 n
n n
( | | | | | |
= + =
( | | |

\ \ \

n 3,7,11,15,... =

2
2
2 2 2 2
1
1
n 2x n 4 n 4 2
2 n 2 2 n
n n
( | | | | | |
= + = +
( | | |

\ \ \

n 2, 6,10,14,... =

2
2
2 2 2 2
1
1
n 2x n 4 n 8
2 n 2 2
n n
( | | | | | |
= + =
( | | |
\ \ \

n 4, 8,12,16,... =

2
2
2 2
1
1
n 2x n 4 n
x cos x dx sen x cos x 0
2 n 2 2
n
( | | | | | |
= + =
( | | |
\ \ \

n 1,5, 9,13,... =
n
2 2 2 2
2 4 4 2 4
a .
n n n
n n
= + + =

n 3,7,11,15,... =
n
2 2 2 2
2 4 4 2 4
a .
n n n
n n
= + + =

n 2, 6,10,14,... =
n
2 2
8
a .
n
=
n 4, 8,12,16,... =
n
a 0. =
n
2 2 2
n 1
3 4 1 (2n 1) ( 1) 1 (2n)
f(x) cos x cos x
4 2 2
(2n 1) (2n)
=
| |
| | | |
= + + |
| |
|
\ \
\

(ver figuras 16, 17 y 18)

2
2
3 4 1
s (x) cos x cos( x)
4 2 2
| | | |
= +
| |
\ \


4
2
3 4 1 1 3
s (x) cos x cos( x) cos x
4 2 2 9 2
| | | | | |
= + +
| | |
\ \ \


6
2
3 4 1 1 3 1 5 1
s (x) cos x cos( x) cos x cos x cos(3 x)
4 2 2 9 2 25 2 18
| | | | | | | |
= + + +
| | | |
\ \ \ \


Ejemplo 38. Sin recurrir a la integracin, encuentre el desarrollo en serie de Fourier
de la funcin
3
f(x) cos (x).sen(2x) = en el intervalo ( , ) .
Solucin.
Usando las identidades trigonomtricas, se puede demostrar que
1
2
sen cos sen( ) sen( ) = + + (

. Por tanto:

3 1
2
1 1
2 4
1 1
4 8
1 cos(2x)
cos (x)sen(2x) cos(x). .sen(2x) cos(x) sen(2x) cos(2x)sen(2x)
2
cos(x)sen(2x) cos(x)sen(4x)
sen(3x) sen(x) sen(5x) sen(3x)
+
= = + (

= +
= + + + ( (

3 1 1
4 8 8
sen(x) sen(3x) sen(5x) = + +

Ejemplo 39. Sea
1
10
1
10
x 0 x 10
f(x)
(20 x) 10 x 20

.
Encuentre el desarrollo en serie de Fourier de:
a. cosenos de f(x) y sela para demostrar que
2
2
k 0
1
8
(2k 1)
=
+
.
Solucin.
( (
= + = + = + + = ( (
( (

10 20 10 20
2 2
2 1 x x 1
0
200 100 2 2 100
0 10 0 10
a xdx (20 x)dx (20x ) (50 200 200 50) 1

10 20
k x k x 2
k
200 20 20
0 10
20 20 20
k x k x k x 1
100 20 20 20
0 10 10
a xcos( )dx (20 x)cos( )dx . k 1,...
xcos( )dx 20 cos( )dx xcos( )dx

(
= + = (
(

(
= + (
(

10 20
20
20x k x 400 k x 400 k x 20x k x 400 k x 1
100 k 20 2 2 20 k 20 k 20 2 2 20
10
k k
0 10
200 k 400 k 400 400 k 400 200 k 1
100 k 2 2 2 2 2 2 k 2 2 2 k
k k k
( sen( ) cos( )) sen( ) ( sen( ) cos( ))
sen( ) cos( ) sen( ) cos(k ) sen(

(
= + + + (
(

= + +

400 k
2 2 2 2
k
k k 4
2 2 2
k
) cos( )
2cos( ) 1 ( 1) . Si k 2n,
(
+
(

(
= =

2n n 2n 4 1 2
2 2 2 2 2 2 2
4n n n
2cos( ) 1 ( 1) 2cos(n ) 2 ( 1) 1

( (
= = = (

Por lo tanto
k
( 1) 1
k x 1 2
2 2 2 10
k
k 1
f(x) cos( )
=
(
= +
(

.
En particular, para x 0 = se obtiene
k
2
( 1) 1
1 2 1 4 1 4 1 1 1
2 2 2 2 2 2 2 2 2 2 8
k (2k 1) (2k 1) (2k 1)
k 1 k 1 k 0 k 0
0 0

+ + +
= = = =
(
= + = = =
(

.
b. senos de f(x) y use la identidad de Parseval para hallar el valor de
4
k 0
1
(2k 1)
=
+
.
Solucin.

(
( = + =
(

(
( = +
(

= +

10 20
k x k x 2
k
200 20 20
0 10
20 20 20
k x k x k x 1
100 20 20 20
0 10 10
10
20x k x 400 k x 400 1
100 k 20 2 2 20 k
k
0
b xsen( )dx (20 x)sen( )dx . k 1,...
xsen( )dx 20 sen( )dx xsen( )dx
( cos( ) sen( ))

(
(
+
(

(
= + + + +
(

=

20
20
k x 20x k x 400 k x
20 k 20 2 2 20
10
k
10
200 k 400 k 400 400 k 400 200 k 400 k 1
100 k 2 2 2 2 k k 2 k k 2 2 2 2
k k
400 k 1
100 2 2 2
k
cos( ) ( cos( ) sen( ))
cos( ) sen( ) cos(k ) cos( ) cos(k ) cos( ) sen( )
sen(

(
+ =
(

400 k 8 k
2 2 2 2 2 2
k k
) sen( ) sen( ).

Por lo tanto
k
sen( )
8 k x 2
2 2 20
k
k 1
f(x) sen( )
=
(
=
(
(

.
Al calcular

20
2
0
f(x) dx (

se obtiene:

20
10 20 10 20
3
3 2
(20 x) 2 2 20 1 1 1 x
100 100 100 3 3 3
0 0 10 0
10
f(x) dx x dx (20 x) dx

(
= + = = ( (

(

.
Aplicando la igualdad de Parseval:
2 2 k k
sen ( ) sen ( )
4 4
64 20 20 5 2 2 1
4 4 3 4 64.3 4 48
k k (2k 1)
k 1 k 1 k 0
.

+
= = =
= = =

4.15. ALGUNAS DEFINICIONES

Definicin 14. Una ecuacin diferencial parcial (EDP) es una relacin funcional
entre una funcin incgnita u de varias variables independientes, dichas variables
independientes y las derivadas parciales de u.

Ejemplo 40.
z
2
y
,
2 2 2
2 2 2
t t t
0
x y z

+ + =

son ecuaciones diferenciales parciales o en derivadas parciales.

Ejemplo 41. Se tienen las siguientes ecuaciones diferenciales parciales:
a.
2
tt xx
u a u 0 = (Ecuacin de onda unidimensional que caracteriza la propagacin de
ondas en algunos medios y las vibraciones mecnicas de una cuerda vibrante)
b.
xx yy zz
u u u 0 + + = (Ecuacin de Laplace que se presenta en el estudio de
potenciales magntico, elctrico, gravitatorio y en el flujo de calor)

Definicin 15. El orden de una EDP es el mayor orden de todas las derivadas
parciales que figuran en la ecuacin.
Ejemplo 42. La ecuacin
4 2
4 2
t t
0
x x

+ =

es una ecuacin diferencial parcial de orden cuatro.

Definicin 16. El grado de una EDP lo da la potencia de la derivada de mayor orden
en dicha ecuacin.

Ejemplo 43.
2
2
z z
3x 6
x
x

=

es una ecuacin diferencial parcial de orden 2 y grado 1.

Observacin 3. Si se utiliza la notacin
x
u u/ x = ,
y
u u/ y = ,
2 2
xx
u u / x = ,
2
xy
u u / x y = ,
se puede escribir una ecuacin diferencial parcial general de segundo orden para una
funcin u de dos variables independientes x e y definida en un dominio D, en la forma

x y xx xy yy
F(x, y,u,u ,u ,u ,u ,u ) 0 = , (17)

relacin que debe verificarse para todo punto (x,y) del interior R del dominio D de la
funcin u.

Definicin 17. Una solucin de (17) en R es cualquier funcin u definida y dos veces
continuamente diferenciable en R que verifica (17) idnticamente para todo (x,y) en
2
R .

4.16. DIFERENCIAS ENTRE LAS ECUACIONES DIFERENCIALES
ORDINARIAS (EDO) Y LAS ECUACIONES DIFERENCIALES
PARCIALES (EDP)

a. En cuanto al tipo de problemas donde se presentan, se puede decir que las
EDO describen problemas unidimensionales, mientras que las EDP describen
problemas de dos o ms dimensiones.
b. En cuanto al espacio de soluciones, se puede decir que en el caso de las EDO el
espacio de soluciones es de dimensin finita, y adems la dimensin de dicho
espacio coincide con el orden de la ecuacin diferencial, mientras que en el caso de
las EDP no necesariamente el espacio de soluciones es de dimensin finita y
adems no depende del orden de la ecuacin. Es bueno hacer notar que en la
mayora de los casos los espacios de soluciones para EDP son de dimensin
infinita.
c. En cuanto a la cantidad de condiciones suplementarias que hay que
agregar a la ecuacin para determinar en forma nica una solucin, en el
caso de las EDO, se necesitarn tantas como el nmero de constantes arbitrarias
que aparezcan en la solucin, o sea, tantas como el orden de la ecuacin. En el
caso de las EDP, an en el caso comparativamente ms simple de las ecuaciones
lineales la situacin es mucho ms complicada. Las soluciones de una EDP
dependen no ya de constantes arbitrarias, sino de funciones arbitrarias.

4.17. SOLUCIONES DE ALGUNAS ECUACIONES DIFERENCIALES
PARCIALES SENCILLAS

Con el objeto de obtener algunas ideas relacionadas con la naturaleza de las
soluciones de ecuaciones diferenciales parciales, considere los siguientes ejemplos:

Ejemplo 44. La ecuacin de segundo orden
2
2
u
0
y

puede ser resuelta integrando dos veces respecto de
y:
u
f(x)
y
, u yf(x) g(x) = +
donde f y g son funciones arbitrarias.

Ejemplo 45. Obtenga soluciones de la ecuacin diferencial parcial

2
2
U
6x 12y
x y
= +

. (18)
Solucin. Aqu la variable dependiente U depende de dos variables independientes x y
y. Para hallar soluciones, se busca determinar U en trminos de x y y, esto es, U(x,y)
Si se escribe (18) como
2
U
6x 12y
x y
| |
= +
|

\
,
se puede integrar con respecto a x, manteniendo y constante, para hallar

2 2
U
3x 12xy F(y)
y
= + +
(19)
donde se ha aadido la constante de integracin, la cual depende de y y as
realmente es una funcin de y denotada por F(y). Ahora se integra (19) con respecto a
y manteniendo x constante para hallar

2 3
U 3x y 4xy F(y)dy G(x) = + + +
(20)
esta vez aadiendo una funcin arbitraria de x dada por G(x). Puesto que la integral de
una funcin arbitraria de y es otra funcin arbitraria de y, se puede escribir (20) como

2 3
U 3x y 4xy H(y) G(x) = + + + . (21)

Esto se puede chequear al sustituirla en (18) y obtener una identidad. Puesto
que (18) es una ecuacin diferencial parcial de segundo orden, y (21) tiene dos
funciones arbitrarias, por analoga con las ecuaciones diferenciales ordinarias se llama
a (21) la solucin general de (18). Usando la misma analoga, se llamara por
supuesto a cualquier solucin obtenida de la solucin general (21) por selecciones
particulares de las funciones arbitrarias, tal como por ejemplo
3
H(y) y , = G(x) sen2x = ,
una solucin particular. Se construye as la siguiente definicin:

Definicin 18. Dada una ecuacin diferencial parcial de orden n, una solucin que
contenga n funciones arbitrarias se llama la solucin general y cualquier solucin
obtenida de esta solucin general por selecciones particulares de las funciones
arbitrarias se llama una solucin particular.

Observacin 4. Como en el caso de ecuaciones diferenciales ordinarias, puede
suceder que existan soluciones singulares las cuales no se puedan obtener de la
solucin general por cualquier seleccin de las funciones arbitrarias.

Como en el caso de ecuaciones diferenciales ordinarias, con frecuencia se
necesitan determinar soluciones de ecuaciones diferenciales parciales que satisfagan
condiciones dadas. Por ejemplo, suponga que se desea resolver la ecuacin diferencial
(18) sujeta a las dos condiciones

2
U(1, y) y 2y , U(x,2) 5x 5 = = . (22)

Entonces de la solucin general (21) y la primera condicin (22) se tiene
2 3 2
U(1, y) 3(1) y 4(1)y H(y) G(1) y 2y = + + + = o
2 3
H(y) y 5y 4y G(1) =
de modo que

2 3 2 3
U 3x y 4xy y 5y 4y G(1) G(x) = + + + . (23)

Si ahora se usa la segunda condicin en (22), se tiene
2 3 2 3
U(x,2) 3x (2) 4x(2) (2) 5(2) 4(2) G(1) G(x) 5x 5 = + + + =
de donde
2
G(x) 33 27x 6x G(1) = + . Usando esto en (23), se obtiene la solucin
2 3 2 3 2
U 3x y 4xy y 5y 4y 27x 6x 33 = + + + .

Se podra usar la misma terminologa de problemas de valor inicial y de
frontera para las ecuaciones diferenciales parciales como se hizo para las ecuaciones
diferenciales ordinarias. Sin embargo, debido a que generalmente hay una
combinacin de condiciones de frontera e iniciales, con frecuencia se refiere a tales
problemas como problemas de valor de frontera.

Ejemplo 46. Hallar una solucin al problema de valor de frontera

2
2
x
U U
2 , U(0, y) 0 , U (x, 0) x
x y x

= + = =

Solucin.
Escribiendo la ecuacin como
U
U 2
x y
| |
=
|

\

e integrando con respecto a x se tiene
U
U 2x F(y)
y
= +
,
la cual es una ecuacin lineal con factor integrante
y
e
. Por tanto,
y y y
(e U) 2xe e F(y)
y

= +

o
y y y
U(x, y) 2x e e F(y)dy e G(x)
= + +
,
donde G(x) es arbitraria. Escribiendo
y y
H(y) e e F(y)dy
,
se tiene

y
U(x, y) 2x H(y) e G(x) = + + . (24)

De U(0, y) 0 = se encuentra
y
H(y) G(0)e = , as (24) llega a ser
y y
U(x, y) 2x G(0)e e G(x) = + .

Diferenciando con respecto a x y colocando y 0 = , se encuentra

3
2
x
x
U (x, 0) 2 G'(x) x o G(x) 2x c
3
= + = = + + ,
de donde,
3
y y
x
U(x, y) 2x G(0)e e 2x c
3
| |
= + + + |
|
\
.
Puesto que c G(0), = entonces se tiene que
3 y
y
x e
U(x, y) 2xe 2x
3
= + .

4.18. SIGNIFICADO GEOMTRICO DE LAS SOLUCIONES
GENERAL Y PARTICULAR

En el ejemplo 45, se obtuvo la solucin general

2 3
U 3x y 4xy H(y) G(x) = + + + . (25)

Suponga ahora que se escogen funciones particulares para H(y) y G(x), y se
reemplaza U por z. Entonces (25) toma la forma z f(x, y) = , la cual se interpreta como
una superficie S en un sistema de coordenadas rectangular o xyz. La superficie est
formada por los puntos con coordenadas (x,y,z) que satisfacen z f(x, y) = .

Para funciones arbitrarias H(y) y G(x), se obtiene una familia de superficies
cada miembro de la cual corresponde a una seleccin particular de H(y) y G(x), esto
es, una solucin particular. La ecuacin diferencial que tenga esto como una solucin
se llama entonces la ecuacin diferencial de la familia de superficies.

4.19. ECUACIONES DIFERENCIALES PARCIALES QUE SURGEN
DE LA ELIMINACIN DE FUNCIONES ARBITRARIAS

Puesto que las soluciones de ecuaciones diferenciales parciales involucran
funciones arbitrarias, parece lgico que se debiera obtener ecuaciones diferenciales
parciales por el proceso inverso de eliminar tales funciones. Esta idea es til porque
ayuda a enriquecer el conocimiento de cmo se pueden resolver ecuaciones
diferenciales parciales. Se van a considerar algunos ejemplos.

Ejemplo 47. Encuentre una ecuacin diferencial parcial de primer orden que tenga
como su solucin general

2
U y F(x) 3x 4y = + , (26)
donde F(x) es una funcin arbitraria de x.
Solucin.
Si se diferencia a (26) con respecto a y, se obtiene

U
2yF(x) 4
y
= +
. (27)
Entonces eliminando F(x) entre (26) y (27), se encuentra la ecuacin deseada
U
y 2U 6x 4y
y
.
Chequeo.
2
U
y 2U y[2yF(x) 4] 2[y F(x) 3x 4y] 6x 4y
y
= + + =
.

Ejemplo 48. Encuentre una ecuacin diferencial parcial de primer orden que tenga
como solucin general
z F(3x 4y) = , (28)
donde F es una funcin arbitraria.
Solucin.
Sea u 3x 4y = . Entonces (28) llega a ser
z f(u) = (29)
Diferenciando (29) con respecto a x, se tiene

z z u
. F'(u).3 3F '(u)
x u x

= = =

. (30)
Diferenciando (29) con respecto a y, se tiene

z z u
. F'(u).( 4) 4F '(u)
y u y

= = =

. (31)
Eliminando F'(u) entre (30) y (31) produce la ecuacin
z z
4 3 0
x y

+ =

.

Ejemplo 49. Encuentre una ecuacin diferencial parcial de segundo orden que tenga
como su solucin general
U xF(y) yG(x) = + , (32)
donde F y G son funciones arbitrarias.
Solucin.
Se puede eliminar F(y) en (32) dividiendo ambos lados de (32) por x y diferenciando el
resultado con respecto a x. Entonces se encuentra
U y
F(y) G(x)
x x x x
( (
= +
( (

, esto es,
U
x U xyG'(x) yG(x)
x
,
la cual se puede escribir

U
x U y[xG'(x) G(x)]
x
. (33)
Si ahora se dividen ambos lados de (33) por y y se diferencia con respecto a y, se
encuentra

2
1 U U U U
x U 0 xy x y U 0
y y x y x x y
( | |
= + =
( |

\
,
la cual construye la ecuacin de segundo orden deseada. Note que la segunda
ecuacin tambin se puede escribir como
2 2

2
U U U U U
xy x y U 0 puesto que
x y x y x y y x

+ = =

.

Se deberan hacer algunas observaciones acerca de los resultados anteriores.

Observacin 5. Al diferenciar las funciones arbitrarias se asume por supuesto que
ellas son diferenciables. De otra manera, no se tiene derecho a diferenciar.

Observacin 6. La ecuacin diferencial obtenida en cada ejemplo representa la
ecuacin diferencial de la familia representada por la solucin general.

Observacin 7. Si una solucin tiene un nmero dado de n funciones arbitrarias, con
frecuencia es fcil escribir una ecuacin diferencial de orden mayor que n teniendo
esta solucin. Por ejemplo, es fcil ver que (32) es una solucin de
4
2 2
U
0
x y
=

pero sta es de orden 4 y no 2. Cuando se busca la ecuacin diferencial, se busca
aquella del menor orden.

4.20. EL MTODO DE SEPARACIN DE VARIABLES

A veces, para una ecuacin diferencial en derivadas parciales lineal
homognea, es posible obtener soluciones particulares en forma de producto
u(x, y) X(x)Y(y) = . El uso de este producto, llamado mtodo de separacin de
variables, permite reducir la ecuacin diferencial en derivadas parciales a varias
ecuaciones diferenciales ordinarias. Con este propsito, se hace notar que

2 2
2 2
u u u u
X' Y , XY ' , X'' Y , XY ''
x y
x y

= = = =

,
en donde las primas indican diferenciales ordinarias.

Ejemplo 50. Halle soluciones en forma de producto de la ecuacin

2
2
u u
4 .
y
x

=
(34)
Solucin.
Si u X(x)Y(y), = entonces (34) se transforma en X'' Y 4XY ' = . Despus de dividir
ambos miembros entre 4XY, se logra separar las variables:
X'' Y '
4X Y
= .
Puesto que el lado izquierdo de esta ecuacin es independiente de y y es idntico al
lado derecho, el cual es independiente de x, se concluye que ambos miembros deben
ser iguales a una constante. En la prctica es conveniente escribir esta constante real
como
2
, o bien como
2
. Se distinguen los casos siguientes:
CASO I: Si
2
ctte , = las igualdades
2
X'' Y '
4X Y
= =
conducen a las ecuaciones
2
X'' 4 X 0 = y
2
Y ' Y 0 = . Estas ltimas tienen las
soluciones
2 x 2 x
1 2
X c e c e

= + y
2
y
3
Y c e
= , respectivamente. As, una solucin

particular de (34) es
2 2 2
2 x 2 x y y 2 x y 2 x
1 2 3 1 1
u XY (c e c e )c e A e B e
+
= = + = + ,
en donde
1 1 3
A c c = y
1 2 3
B c c = .
CASO II. Si
2
ctte = , las igualdades
2
X'' Y '
4X Y
= = conducen a
2
X'' 4 X 0 + = y
2
Y ' Y 0 + = . Puesto que las soluciones de estas ecuaciones son
4 5
X c cos(2 x) c sen(2 x) = + y
2
y
6
Y c e
= ,
respectivamente, otra solucin de (34) es
2 2
y y
2 2
u A e cos(2 x) B e sen(2 x)

= + , en
donde
2 4 6
A c c = y
2 5 6
B c c = .
CASO III. Si ctte 0 = , se tiene que X'' 0 = y Y ' 0 = . En este caso
7 8
X c x c = + y
9
Y c = , de modo que
3 3
u A x B = + en donde
3 7 9
A c c = y
2 8 9
B c c = .

Se deja como ejercicio verificar que las tres soluciones obtenidas satisfacen la
ecuacin dada.

Ejemplo 51. Halle una solucin en forma de producto de

2
2
u u
k , k 0
t
x

= >

que satisfaga las condiciones u(0, t) 0, u(L, t) 0 = = .
Solucin.
Si u XT = , se puede escribir la ecuacin dada como
2
X'' T '
X kT
= =
lo que conduce a
2
X'' X 0 + = y
2
T ' k T 0 + = , de modo que se tiene
1 2
X c cos( x) c sen( x) = + y
2
k t
3
T c e

= ,
respectivamente. Ahora bien, puesto que u(0, t) X(0)T(t) 0 = = , u(L, t) X(L)T(t) 0 = = se
debe tener X(0) 0 = y X(L) 0 = . Estas son condiciones de frontera para la ecuacin
diferencial ordinaria. Aplicando la primera de tales condiciones, resulta de inmediato
1
c 0 = . Por lo tanto,
2
X c sen( x) = . La segunda condicin de frontera implica ahora que
2
X(L) c sen( L) 0 = = . Si
2
c 0 = , entonces X 0 = , de modo que u 0 = . Para obtener una
solucin no trivial u, se debe tener
2
c 0 y entonces la ltima ecuacin se satisface
cuando sen( L) 0 = . Entonces L n = o
n
L
= ,
siendo n 1,2,3,... =
Por tanto,
2 2 2 2
k t k(n / L )t
2 3 n
n
u (c sen( x)(c e ) A e sen x
L

| |
= =
|
\

satisface la ecuacin dada y ambas condiciones adicionales. El coeficiente
2 3
c c se
reescribe como
n
A para recalcar que se obtiene una solucin diferente para cada n. Se
deja como ejercicio verificar que si se usa
2
0 , no se llega a una solucin que
satisfaga u(0, t) 0 = y u(L, t) 0 = .

4.21. PRINCIPIO DE SUPERPOSICIN

TEOREMA 10. Si
1 2 n
u , u , ..., u son soluciones de una ecuacin diferencial paracial
lineal homognea, entonces la combinacin lineal
1 1 2 2 k k
u c u c u ... c u = + + + , donde los

i
c , i 1,2,...,k = son constantes, tambin es una solucin.

En la prxima seccin se supondr formalmente que cada vez que se tenga un
conjunto infinito
1 2
u , u , ... de soluciones de una ecuacin lineal homognea, an se
puede obtener otra solucin u formando la serie infinita
k
k 1
u u
=
=
.

Ejemplo 52. En virtud del principio de superposicin, la funcin definida mediante la
serie
2 2 2
k(n / L )t
n
n 1
n
u A e sen x
L

=
| |
=
|
\

debe satisfacer tambin, aunque sea formalmente, a la ecuacin del ejemplo anterior,
junto con las condiciones dadas.

4.22. PROBLEMAS DE CONDICIN DE FRONTERA

Las siguientes ecuaciones diferenciales parciales lineales

2
2
u u
k , k 0,
t
x

= >

2 2
2
2 2
u u
a
x t

=

,
2 2
2 2
u u
0
x v

+ =

,
desempean un papel importante en muchas reas de fsica e ingeniera. Las primeras
dos ecuaciones son conocidas como ecuacin del calor en una dimensin y
ecuacin de onda en una dimensin respectivamente. Una dimensin se refiere al
hecho de que x denota una dimensin espacial, en tanto que t generalmente
representa tiempo. La ltima ecuacin se llama ecuacin de Laplace. Se usar el
mtodo de separacin de variables para resolver varios problemas aplicados, cada uno
de los cuales es descrito por una de las ecuaciones anteriores adems de ciertas
condiciones adicionales. Estas condiciones adicionales consisten en:
a. Condiciones de frontera: u o u x especificada para x cons tante = ; u o u y
especificada para y cons tante = ,
b. Condiciones iniciales: u en t 0 = para la primera ecuacin, o bien, u y u t en
t 0 = para la segunda ecuacin.

La descripcin matemtica colectiva de un problema de esta naturaleza es
conocida como problema de condicin en la frontera.

En una frontera se pueden especificar los valores de una de las siguientes
funciones:

x x
(i) u , (ii) u , o bien (iii) u hu , h cons tante + .

Una condicin a la frontera del primer tipo (i), se llama condicin de
Dirichlet; del segundo tipo (ii), condicin de Neuman, y del tercer tipo (iii),
condicin de Robin.

4.23. ECUACIONES DIFERENCIALES PARCIALES LINEALES DE
SEGUNDO ORDEN

Ya se han visto cmo se resuelven algunos problemas de frontera que
involucran ecuaciones diferenciales parciales. Se enfocar ahora el estudio sobre cierto
tipo de ecuaciones diferenciales parciales importantes que surgen en varios problemas
aplicados. Hay tres tipos importantes de problemas que involucran:
a. Conduccin o difusin de calor
b. Vibraciones u oscilaciones
c. Temperatura estacionaria, potencial elctrico o gravitacional

Definicin 19. Si u representa la variable dependiente, x e y las variables
independientes, entonces la forma general de una ecuacin diferencial parcial
lineal de segundo orden con dos variables independientes x e y, es

2 2 2
2 2
u u u u u
A B C D E Fu G
x y x y
x y

+ + + + + =

, (35)
en donde A, B, C, D, E, F y G son funciones de x e y.

Observacin 8. Cuando G(x, y) 0 = , la ecuacin se llama homognea; en cualquier
otro caso es no homognea.
Ejemplo 53.
2 2
2 2
u u
0
x y

+ =

es homognea.
2
2
u u
xy
y
x

=
no es homognea.

Definicin 20. Una solucin de una ecuacin (35) en derivadas parciales con dos
variables independientes x e y es una funcin u(x,y) que posee todas las derivadas
parciales que indica la ecuacin y que la satisface en alguna regin del plano xy.

Una ecuacin en derivadas parciales, lineal de segundo orden con dos variables
independientes y con coeficientes constantes, puede pertenecer a uno de tres tipos
generales. Esta clasificacin slo depende de los coeficientes de las derivadas de
segundo orden. Naturalmente, se supone que al menos uno de los coeficientes A, B y C
no es cero.

Definicin 21. La ecuacin diferencial parcial lineal y de segundo orden

2 2 2
2 2
u u u u u
A B C D E Fu 0
x y x y
x y

+ + + + + =

,
en donde A, B, C, D, E, F y G son constantes reales es

a. parablica, si
2
B 4AC 0 = .
Ejemplo 54.
2
2
u u
3 0 (A 3,B 0, C 0)
y
x

= = = =
. Esta ecuacin modela la distribucin

de temperatura en una varilla (ecuacin de calor unidimensional).
b. hiperblica, si
2
B 4AC 0 > .
Ejemplo 55.
2 2
2 2
u u
0 (A 1,B 0, C 1)
x y

= = = =

. Esta ecuacin modela las vibraciones
en una cuerda (ecuacin de onda unidimensional).
c. elptica, si
2
B 4AC 0 < .
Ejemplo 56.
2 2
2 2
u u
0 (A 1,B 0, C 1)
x y

+ = = = =

. Esta ecuacin modela el potencial o
distribucin de temperatura estacionaria no dependiente del tiempo (ecuacin de
Laplace).

Estas ecuaciones tienen mtodos para hallar su solucin, entre los que
destacan:
a. Mtodo de separacin de variables
b. Mtodo de la transformada de Laplace
c. Mtodo de la integral de Fourier

4.24. ECUACIN DE CALOR UNIDIMENSIONAL

A continuacin se empezar formulando las ecuaciones del flujo de calor que
describen la transferencia de energa trmica (o calorfica). La energa trmica est
causada por la agitacin de las molculas. Existen dos procesos bsicos que
intervienen en la transferencia de la energa trmica: conduccin y conveccin. La
conduccin resulta de las colisiones de molculas colindantes, por las que la energa
cintica de vibracin de una molcula se transfiere a su vecina ms cercana. La
energa trmica se transporta de esta manera por conduccin incluso si las molculas
mismas no se mueven de forma apreciable de su posicin. Por otro lado, una molcula
vibrante que se traslada de una regin a otra, lleva su energa trmica con ella. Este
segundo fenmeno se conoce con el nombre de conveccin. Con objeto de comenzar el
estudio con problemas relativamente simples, se estudiar el flujo de calor slo en los
casos en que la conduccin de energa trmica predomina sobre la conveccin. ste es
el caso del flujo de calor en slidos, o incluso en fluidos (lquidos y gases) en los que la
velocidad es suficientemente pequea.

4.25. DEDUCCIN DE LA CONDUCCIN DEL CALOR EN UNA
VARILLA UNIDIMENSIONAL

Suponga que una varilla circular delgada de longitud L tiene un rea transversal
A y coincide con el eje x en el intervalo 0,L (

(ver figura 19). Suponga lo siguiente:
x
x=0
x x+ x
x=L
Q

Figura 19. Flujo de calor en una varilla delgada

a. El flujo de calor dentro de la varilla toma lugar solo en la direccin x.
b. La superficie lateral, o curva, de la varilla est aislada; es decir, el calor no escapa
de su superficie.
c. No se genera calor dentro de la varilla.
d. La varilla es homognea; es decir, su masa por unidad de volumen es una
constante.
e. El calor especfico y la conductividad trmica K del material de la varilla, son
constantes.
A fin de obtener la ecuacin diferencial parcial que se cumple para la
temperatura u(x,t), son necesarias dos leyes empricas de conduccin del calor:
a. La cantidad de calor Q en un elemento de masa m es
Q mu = , (36)
donde u es la temperatura del elemento.
b. La rapidez de flujo de calor
t
Q por la seccin transversal indicada en la figura, es
proporcional al rea A de la seccin transversal y la derivada parcial con respecto a
x de la temperatura:

t x
Q KAu = . (37)

Puesto que el calor fluye en la direccin de temperatura decreciente, el signo
menos en (37) se usa para asegurar que
t
Q sea positiva para
x
u 0 < (flujo de calor a
la derecha) y negativa para
x
u 0 > (flujo de calor a la izquierda). Si la parte circular de
la varilla mostrada en la figura entre x y x x + es muy delgada, entonces u(x,t) se
toma como la temperatura aproximada en cada punto del intervalo. Ahora la masa de
la porcin es m A x = , y de (36) se deduce que la cantidad de calor en sta es
Q A xu = . (38)
Adems, cuando el calor fluye en la direccin positiva de x, se ve de (37) que el
calor se acumula en el corte con razn neta

x x x x
KAu (x, t) KAu (x x, t) KA u (x x, t) u (x, t) + = + ( (

. (39)

Al diferenciar (38) con respecto a t, se ve que esta razn neta se expresa
tambin mediante

t t
Q xu = . (40)

Al igualar (39) y (40), se obtiene

x x
t
u (x x, t) u (x, t) K
u
x
+
=

. (41)

Por ltimo, al tomar el lmite de (41) cuando x 0 , se obtiene la ecuacin
xx t
(K / )u u = (La definicin de la segunda derivada parcial viene dada por la
expresin
x x
xx
x 0
u (x x, t) u (x, t)
u lm
x
+
=
. Normalmente se establece que k K / = , y

esta constante se llama difusividad trmica.

4.26. PROBLEMAS ILUSTRATIVOS DE LA ECUACIN DE CALOR
UNIDIMENSIONAL

Ejemplo 57. Considere una barra delgada o varilla de largo L con una distribucin
longitudinal de temperatura f(x) y cuyos extremos se mantienen a una temperatura
constante de cero grados en todo instante. Si la varilla satisface las suposiciones
enunciadas en la seccin anterior, entonces la temperatura u(x,t) en la varilla se
determina a partir del problema de valores en la frontera

2
2
u u
k , k 0, 0 x L, t 0
t
x

= > < < >
, (42)
u(0, t) 0, u(L, t) 0, t 0 = = > , (43)
u(x, 0) f(x), 0 x L = < < . (44)

La constante k es proporcional a la conductividad trmica y se llama
difusividad trmica. Este problema representa el caso de una ecuacin de calor
homognea con condicin inicial no nula y condiciones de frontera nulas.
Solucin.
Usando el producto u XT = y
2
como constante de separacin, se llega a
2
X'' T '
X kT
= =
y

2
X'' X 0, X(0) 0, X(L) 0, + = = = (45)

2
T ' k T 0 + = (46)
Ya se ha obtenido resolviendo (45) la solucin de la ecuacin (42) sujeta a las
condiciones de frontera (43). Se vi que (45) tena una solucin no trivial slo si el
parmetro adoptaba los valores

n
, n 1,2,3,...
L
= = (47)
Las correspondientes soluciones de (45) eran entonces

1
n
X c sen x , n 1,2,3,...
L
| |
= =
|
\
(48)

Los valores (47) para los cuales (45) tiene una solucin no trivial son conocidos
como valores caractersticos o como valores propios. Las soluciones (48) se
denominan funciones caractersticas o funciones propias. Se hace notar que para
un valor de distinto de los dados en (47), la nica solucin de (45) es la funcin
X 0 = . A su vez, esto implicara que una funcin que satisface (42) y (43) en u 0 = no
es una solucin del problema original de condicin en la frontera cuando f(x) 0 .
Naturalmente se supone que se cumple esta ltima condicin. Puesto que la solucin
de (46) es
2 2 2 2
k t k(n /L )t
3 3
T c e c e

= = , los productos

2 2 2
k(n /L )t
n n
n
u XT A e sen x
L

| |
= =
|
\
(49)
satisfacen la ecuacin diferencial (42) y las condiciones de frontera (43) para cada
valor del entero positivo n. Por conveniencia se ha reemplazado la constante
1 3
c c por
n
A . Para que las funciones dadas en (49) satisfagan la condicin inicial (44), se
tendra que elegir coeficientes constantes
n
A de modo que

n n
n
u (x, 0) f(x) A sen x
L
| |
= =
|
\
. (50)

En general, la condicin (50) no se satisface para cualquier f. Por lo tanto, hay
que admitir que (47) no es una solucin del problema dado. Sin embargo, por el
principio de superposicin, la funcin

2 2 2
k(n / L )t
n n
n 1 n 1
n
u(x, t) u A e sen x
L

= =
| |
= =
|
\

(51)
tambin satisface (42) y (43). Sustituyendo t 0 = en (51) se obtiene
n
n 1
n
u(x, 0) f(x) A sen x
L
=
| |
= =
|
\
.

Se ve que esta ltima expresin es el desarrollo en serie de senos en medio
intervalo. Por consiguiente si se hace la identificacin
n n
A b = , n 1,2,3,... = , se tiene
que

L
n
0
2 n
A f(x)sen x dx
L L
| |
=
|
\
.
Se concluye que la solucin del problema de condicin de frontera descrito en (42),
(43) y (44) est dada por la serie infinita

L
2 2 2
k(n / L )t
0
n 1
2 n n
u(x, t) f(x)sen x dx e sen x
L L L

=
| | | | | |
=
| | |
\ \ \

.

Ejemplo 58. Ecuacin de calor homognea con condicin inicial no nula y condiciones
de frontera no nulas y constantes.

2
2
U U
k , 0 x L, k 0, t 0
t
x

= < < > >
,

1 2
U(0, t) K , U(L, t) K , t 0 = = > ,
U(x, 0) f(x), 0 x L = < < .
Solucin.
Sea U(x, t) W(x, t) (x) = + . Se tiene
t t xx xx
U (x, t) W(x, t) , U (x, t) W (x, t) ''(x) = = + .
El problema se convierte en

2
2
W W
k ''(x) , 0 x L, k 0, t 0
t
x
| |

+ = < < > > |
|

\
,

1 2
W(0, t) K (0), W(L, t) K (L), t 0 = = > ,
W(x, 0) f(x) (x), 0 x L = < < .
Condiciones sobre (x) :
1 2
''(x) 0, (0) K , (L) K = = = .
Entonces:

2 1
1 2 1
K K
(x) x , (0) K , (L) K K , .
L
= + = = = =
El problema es:

2
2
W W
k , 0 x L, k 0, t 0
t
x
| |

= < < > > |
|

\
,
W(0, t) 0, W(L, t) 0, t 0 = = > ,

K K
2 1
1
L
W(x, 0) f(x) x K , 0 x L
= < < .
Solucin:
2 n
k( ) t
2 1 n x L
1 n
L
n 1
K K
U(x, t) x K a sen( )e
L
= + +
,
donde

L
n x
n
L
0
2
a f(x) (x) sen( )dx
L

= (

Ejemplo 59. La temperatura u(x,t) en una varilla no aislada de longitud L, con
condiciones de frontera no nulas se determina a partir del problema de valores en la
frontera

xx t
ku f(x, t) u , k 0, 0 x L, t 0 + = > < < > ,

1 2
u(0, t) (t), u(L, t) (t), t 0 = = > ,
u(x, 0) (x), 0 x L = < < .
Solucin.
En primer lugar, se puede escribir u(x, t) F(x, t) w(x, t) = + , donde F(x,t) es una funcin
que debe resolver el problema
xx 1 2
F (x, t) 0, F(0, t) (t), F(L, t) (t) = = = . De modo que
xx 1 2
F (x, t) 0 F(x, t) g (t) xg (t) = = + . Al imponer la primera condicin de borde se tiene
1 1
F(0, t) (t) g (t) = = . Al imponer la segunda condicin de borde se tiene
(t) (t)
2 1
1 2 2 1 2 2
L
F(x, t) (t) xg (t) F(L, t) (t) (t) Lg (t) g (t)

= + = = + = .
Por lo tanto se obtiene la funcin definida por
x
1 2 1
L
F(x, t) (t) (t) (t) = + (

.
Por otro lado w(x,t) es una funcin incgnita auxiliar cuya determinacin se estudiar a
continuacin. Como
1
F(0, t) (t) = y
2
F(L, t) (t) = por construccin, resulta que
w(0, t) 0 = y w(L, t) 0 = . Adems w(x, 0) u(x, 0) F(x, 0) (x) F(x, 0) = = en 0 x L < < , y
xx t t xx t t
kw f(x, t) w F (x, t) kw w F (x, t) f(x, t) + = + = + (

.
Si se denotan H(x) (x) F(x, 0) = y
t
G(x, t) f(x, t) F (x, t) = , el problema original queda
convertido como

xx t
kw G(x, t) w , k 0, 0 x L, t 0 + = > < < > ,
w(0, t) 0, w(L, t) 0, t 0 = = > ,
w(x, 0) H(x), 0 x L = < < .
Este problema representa el caso de una ecuacin de calor no homognea con
condicin inicial no nula y condiciones de frontera nulas.
Si se ensaya con una solucin de la forma
n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
,
y se expresa
n x
n
L
n 1
G(x, t) G (t)sen( )
=
=
,
con

L
n x
n
L
0
2
G (t) G(x, t)sen( )dx, n 1,...
L

= =
,
entonces
2 n n x
xx n
L L
n 1
w (x, t) w (t)( ) sen( )

=
=
,
' n x
t n
L
n 1
w (x, t) w (t)sen( )
=
=
.
Sustituyendo en la ecuacin
xx t
kw G(x, t) w + = se tiene
2 ' n n x n x n x
n n n
L L L L
n 1 n 1 n 1
k w (t)( ) sen( ) G (t)sen( ) w (t)sen( )

= = =
+ =

.
Reacomodando trminos:
Ecuacin lineal de primer orden
' 2 ' 2 n n x n
n n n n n n
L L L
n 1
w (t) k( ) w (t) G (t) sen( ) 0 w (t) k( ) w (t) G (t)

=
(
+ = + =

.
Como
n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
,
entonces de la condicin inicial del problema se tiene
n x
n
L
n 1
H(x) w (0)sen( )
=
=

(serie de senos de Fourier)
y por tanto se obtiene el coeficiente
n
w (0) como

L
n x
n
L
0
2
W (0) H(x)sen( )dx
L

=
.
De manera que hay que resolver una ecuacin lineal ordinaria no homognea de
primer orden sujeta a una condicin:

L
' 2 n n x
n n n n
L L
0
2
w (t) k( ) w (t) G (t) , W (0) H(x)sen( )dx
L

+ = =

Mtodo 1.
Solucin general:
2 2 2 n n n
k( ) t k( ) t k( ) t
L L L
n n n
w (t) e G (t)e dt C e M (t) C

(
= + = + (
(

.
Valor C:
n n
C w (0) M (0) = . Solucin particular:
2 n
k( ) t
L
n n n n
w (t) e M (t) w (0) M (0)
= + (

.
Resumen de pasos a seguir:
1. Construccin de G(x,t).
t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
3. Construccin de H(x). H(x) (x) F(x, 0) =
4. Clculo de
n
w (0)
5. Clculo de
n n
M (t) y M (0) .
6. Clculo de los
n
w (t).
7. Construccin de w(x,t).
8. Construccin de u(x,t).

Mtodo 2.
' 2 2 n n
n n n n n n n
L L
{w (t)} k( ) {w (t)} {G (t)} sw (s) w (0) k( ) w (s) G (s)

+ = + = L L L .
2 n n n n n
n n n n
L 2 2 2 n n n
L L L
G (s) w (0) G (s) w (0)
s k( ) w (s) G (s) w (0) w (s)
s k( ) s k( ) s k( )

+
(
+ = + = = +

+ + +
.
2 2 n n
k( ) t k( ) t
-1 -1 n
L L
n n n n n
2 n
L
G (s)
w (t) {w (s)} w (0)e P (t) w (0)e
s k( )

= = + = +
`
+

)
L L .

t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de { }
n n
G (t) y P (t) L .
6. Clculo de los
n
w (t).

Solucin final:
n x x
1 2 1 n
L L
n 1
u(x, t) F(x, t) w(x, t) (t) (t) (t) w (t)sen( )
=
= + = + + (

.

Ejemplo 60. (Caso particular del ejemplo 59) Ecuacin de calor homognea con
condicin inicial no nula y condiciones de frontera no nulas.
Solucin.
En este caso
t
G(x, t) F (x, t) = y se aplican cualquiera de los mtodos descritos en el
ejemplo 59.

Ejemplo 61. La temperatura u(x,t) en una varilla de longitud L con extremos aislados.
Paso 1. Planteamiento del problema.

xx t
ku u , k 0, 0 x L, t 0 = > < < > ,

x x
u (0, t) 0, u (L, t) 0, t 0 = = > ,
u(x, 0) f(x), 0 x L = < < .
Paso 2. Planteamiento y resolucin de las ecuaciones diferenciales ordinarias.
Si se ensaya con soluciones del tipo u(x, t) X(x)T(t) = con X(x) 0 , T(t) 0 se tiene
kX'' T XT ' = . Al dividir cada miembro entre kXT e igualando a una constante negativa
se obtiene
2
X'' T '
X kT
= = ,
lo que origina las dos ecuaciones diferenciales ordinarias:

2
2 2 k t
1 2 3
X'' X 0 X(x) c cos( x) c sen( x) , T ' k T 0 T(t) c e

+ = = + + = =
Paso 3. Utilizacin de las condiciones de borde.

2
X'' X 0 , X'(0) 0 , X'(L) 0 + = = = (Problema de Sturm-Liouville clsico)
Solucin.
n x
n n
L
X (x) a cos( ) n 0,1,2,...
= =
Paso 4. Aplicacin del principio de superposicin.
2 n
k( ) t
n x L
n n
L
u (x, t) b cos( )e

= ,
entonces
2 n
k( ) t
n x L
n n
L
n 0 n 0
u(x, t) u (x, t) b cos( )e

= =
= =

.
Paso 5. Utilizacin de la condicin inicial.
0 n x n x n x
n 0 n n
L L L
n 0 n 1 n 1
A
u(x, 0) f(x) b cos( ) b b cos( ) A cos( )
2

= = =
= = = + = +

,
donde

L L
n x
0 n
L
0 0
2 2
A f(x)dx , A f(x)cos( )dx
L L

= =

.
Serie de cosenos de Fourier para f(x)
Paso 6. Solucin del problema planteado.

2 n L
k( ) t
n x n x L
L L
0
n 0
2
u(x, t) f(x)cos( )dx cos( )e
L

=
| |
=
|
\

.

Ejemplo 62.

xx t
ku f(x, t) u , k 0, 0 x L, t 0 + = > < < > ,

x x
u (0, t) u (L, t) (t), t 0 = = > ,
u(x, 0) (x), 0 x L = < < .
Solucin.
xx x x
F (x, t) 0, F (0, t) F (L, t) (t) = = = . De modo que
xx x 2 1 2
F (x, t) 0 F (x, t) g (t) F(x, t) g (t) xg (t) = = = + . Al imponer la primera condicin
de borde se tiene
x 2
F (0, t) (t) g (t) = = . Al imponer la segunda condicin de borde se
tiene
1 x
F(x, t) g (t) x (t) F (L, t) (t) (t) = + = = . Por lo tanto se obtiene la funcin
definida por
1
F(x, t) g (t) x (t) = + , donde
1
g (t) es una funcin arbitraria que por
comodidad se asumir
1
g (t) 0 = , de modo que F(x, t) (t)x = . Por otro lado w(x,t) es
una funcin incgnita auxiliar cuya determinacin se estudiar a continuacin.

Como
x x
F (0, t) F (L, t) (t) = = , resulta que
x
w (0, t) 0 = y
x
w (L, t) 0 = . Adems se tiene
que w(x, 0) u(x, 0) F(x, 0) (x) F(x, 0) = = en 0 x L < < , y
xx t t xx t t

.

t
convertido como

xx t
kw G(x, t) w , k 0, 0 x L, t 0 + = > < < > ,

x x
w (0, t) 0, w (L, t) 0, t 0 = = > ,
w(x, 0) H(x), 0 x L = < < .
condicin inicial no nula y extremos aislados.

n x
n
L
n 0
w(x, t) w (t)cos( )
=
=
,
y se expresa
n x
n
L
n 0
G(x, t) G (t)cos( )
=
=
,
con

L
n x
n
L
0
2
G (t) G(x, t)cos( )dx, n 0,...
L

= =
,
se tiene
2 n n x
xx n
L L
n 0
w (x, t) w (t)( ) cos( )

=
=
,
' n x
t n
L
n 0
w (x, t) w (t)cos( )
=
=
.
xx t
2 ' n n x n x n x
n n n
L L L L
n 0 n 0 n 0
k w (t)( ) cos( ) G (t)cos( ) w (t)cos( )

= = =
+ =

.
' 2 ' 2 n n x n
n n n n n n
L L L
n 1
w (t) k( ) w (t) G (t) cos( ) 0 w (t) k( ) w (t) G (t)

=
(
+ = + =

.

Como
n x
n
L
n 0
w(x, t) w (t)cos( )
=
=
,
n x
n
L
n 0
H(x) w (0)cos( )
=
=

(serie de cosenos de Fourier)
n
w (0) como

L
n x
n
L
0
2
W (0) H(x)cos( )dx
L

=
.

L
' 2 n n x
n n n n
L L
0
2
w (t) k( ) w (t) G (t) , W (0) H(x)cos( )dx
L

+ = =

Mtodo 1.
Solucin general:
2 2 2 n n n
k( ) t k( ) t k( ) t
L L L
n n n

(
= + = + (
(

.
Valor C:
n n
C w (0) M (0) = . Solucin particular:
2 n
k( ) t
L
n n n n
w (t) e M (t) w (0) M (0)
= + (

.

t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de
n n
M (t) y M (0) .
6. Clculo de los
n
w (t).

Mtodo 2.
' 2 2 n n
n n n n n n n
L L

+ = + = L L L .
2 n n n n n
n n n n
L 2 2 2 n n n
L L L
G (s) w (0) G (s) w (0)
s k( ) w (s) G (s) w (0) w (s)
s k( ) s k( ) s k( )

+
(
+ = + = = +

+ + +
.
2 2 n n
k( ) t k( ) t
-1 -1 n
L L
n n n n n
2 n
L
G (s)
w (t) {w (s)} w (0)e P (t) w (0)e
s k( )

= = + = +
`
+

)
L L .

t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de { }
n n
G (t) y P (t) L .
6. Clculo de los
n
w (t).

Solucin final.
n x
1 0 n
L
n 1
u(x, t) F(x, t) w(x, t) (t)x w (t) w (t)cos( )
=
= + = + +
.

Ejemplo 63.

xx t
ku u , k 0, 0 x L, t 0 = > < < > ,

x
u(0, t) 0, u (L, t) 0, t 0 = = > ,
u(x, 0) f(x), 0 x L = < < .
se obtiene
2
X'' T '
X kT
= = ,

2
2 2 k t
1 2 3

+ = = + + = =

2
X'' X 0 , X(0) 0 , X'(L) 0 + = = = (Problema de Sturm-Liouville clsico)
Solucin.
(2n 1) x
n n
2L
X (x) a sen( ) n 0,1,2,...
+
= =
(2n 1) 2
k( ) t
(2n 1) x
2L
n n
2L
u (x, t) b sen( )e
+
+
= ,
(2n 1) 2
k( ) t
(2n 1) x
2L
n n
2L
n 0 n 0
u(x, t) u (x, t) b sen( )e
+
+
= =
= =

.
(2n 1) x
n
2L
n 0
u(x, 0) f(x) b sen( )
+
=
= =
, donde

L
(2n 1) x
n
2L
0
2
b f(x)sen( )dx
L
+
=
.
Serie de senos de Fourier para f(x).

(2n 1) 2 L
k( ) t
(2n 1) x (2n 1) x
2L
2L 2L
0
n 0
2
u(x, t) f(x)sen( )dx sen( )e
L
+
+ +
=
| |
=
|
\

.

Ejemplo 64.

xx t
ku f(x, t) u , k 0, 0 x L, t 0 + = > < < > ,

1 x 2
u(0, t) (t), u (L, t) (t), t 0 = = > ,
u(x, 0) (x), 0 x L = < < .
Solucin.
xx 1 x 2
F (x, t) 0, F(0, t) (t), F (L, t) (t) = = = . De modo que
xx 1 2
1 1
1 2 x 2 2
F(x, t) (t) xg (t) F (L, t) (t) g (t) = + = = . Por tanto se obtiene la funcin definida
por
1 2
F(x, t) (t) (t)x = + . Por otro lado w(x,t) es una funcin incgnita auxiliar cuya
determinacin se estudiar a continuacin.

Como
1
F(0, t) (t) = y
x 2
F (L, t) (t) = por construccin, resulta que w(0, t) 0 = y
x
w (L, t) 0 = . Adems w(x, 0) u(x, 0) F(x, 0) (x) F(x, 0) = = en 0 x L < < , y
xx t t xx t t

.

t
convertido como

xx t
kw G(x, t) w , k 0, 0 x L, t 0 + = > < < > ,

x
w(0, t) 0, w (L, t) 0, t 0 = = > ,
w(x, 0) H(x), 0 x L = < < .

condicin inicial no nula, condicin de frontera nula y extremo aislado.

(2n 1) x
n
2L
n 0
w(x, t) w (t)sen( )
+
=
=
,
y se expresa
(2n 1) x
n
2L
n 0
G(x, t) G (t)sen( )
+
=
=
,
con

L
(2n 1) x
n
2L
0
2
G (t) G(x, t)sen( )dx, n 0,...
L
+
= =
,
(2n 1) (2n 1) x 2
xx n
2L 2L
n 0
w (x, t) w (t)( ) sen( )
+ +
=
=
,
(2n 1) x '
t n
2L
n 0
+
=
=
.
xx t
(2n 1) (2n 1) x (2n 1) x (2n 1) x 2 '
n n n
2L 2L 2L 2L
n 0 n 0 n 0
k w (t)( ) sen( ) G (t)sen( ) w (t)sen( )

+ + + +
= = =
+ =

.
(2n 1) x ' 2 ' 2 n n
n n n n n n
L 2L L
n 0
w (t) k( ) w (t) G (t) sen( ) 0 w (t) k( ) w (t) G (t)
+

=
(
+ = + =

.
Como
(2n 1) x
n
2L
n 0
w(x, t) w (t)sen( )
+
=
=
,
entonces de la condicin inicial dada en el problema se tiene
(2n 1) x
n
2L
n 0
H(x) w (0)sen( )
+
=
=

n
w (0) como

L
(2n 1) x
n
2L
0
2
W (0) H(x)sen( )dx
L
+
=
.

L
(2n 1) x ' 2 n
n n n n
L 2L
0
2
w (t) k( ) w (t) G (t) , W (0) H(x)sen( )dx
L
+
+ = =
.

Mtodo 1.
Solucin general:
(2n 1) (2n 1) (2n 1) 2 2 2
k( ) t k( ) t k( ) t
2L 2L 2L
n n n
+ + +

(
= + = + ( (

(

.
Valor C:
n n
C w (0) M (0) = . Sol. particular:
(2n 1) 2
k( ) t
2L
n n n n
w (t) e M (t) w (0) M (0)
+
= + (

.

t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de
n n
M (t) y M (0) .
6. Clculo de los
n
w (t).

Mtodo 2.
(2n 1) (2n 1) ' 2 2
n n n n n n n
2L 2L
+ +
+ = + = L L L .
(2n 1) 2
n n n
2L
n n n n
n
(2n 1) (2n 1) (2n 1) 2 2 2
2L 2L 2L
s k( ) w (s) G (s) w (0)
G (s) w (0) G (s) w (0)
w (s)
s k( ) s k( ) s k( )
+
+ + +
(
+ = +

+
= = +
+ + +
(2n 1) (2n 1) 2 2
k( ) t k( ) t
-1 -1 n
2L 2L
n n n n n
(2n 1) 2
2L
G (s)
w (t) {w (s)} w (0)e P (t) w (0)e
s k( )
+ +

+

= = + = +
`
+

)
L L .

t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de { }
n n
G (t) y P (t) L .
6. Clculo de los
n
w (t).
Solucin final:
(2n 1) x
1 2 n
2L
n 0
u(x, t) F(x, t) w(x, t) (t) (t)x w (t)sen( )
+
=
= + = + +
.

Ejemplo 65.

xx t
ku u , k 0, 0 x L, t 0 = > < < > ,

x
u (0, t) 0, u(L, t) 0, t 0 = = > ,
u(x, 0) f(x), 0 x L = < < .
se obtiene
2
X'' T '
X kT
= = ,

2
2 2 k t
1 2 3

+ = = + + = =

2
X'' X 0 , X'(0) 0 , X(L) 0 + = = = (Problema de Sturm-Liouville clsico)
Solucin.
(2n 1) x
n n
2L
X (x) b cos( ) n 0,1,2,...
+
= =
(2n 1) 2
k( ) t
(2n 1) x
2L
n n
2L
u (x, t) a cos( )e
+
+
= ,
(2n 1) 2
k( ) t
(2n 1) x
2L
n n
2L
n 0 n 0
u(x, t) u (x, t) a cos( )e
+
+
= =
= =

.
(2n 1) x
n
2L
n 0
u(x, 0) f(x) a cos( )
+
=
= =
, donde

L
(2n 1) x
n
2L
0
2
a f(x)cos( )dx
L
+
=
.
Serie de cosenos de Fourier para f(x).

(2n 1) 2 L
k( ) t
(2n 1) x (2n 1) x
2L
2L 2L
0
n 0
2
u(x, t) f(x)cos( )dx cos( )e
L
+
+ +
=
| |
=
|
\

.

Ejemplo 66.

xx t
ku f(x, t) u , k 0, 0 x L, t 0 + = > < < > ,

x 1 2
u (0, t) (t), u(L, t) (t), t 0 = = > ,
u(x, 0) (x), 0 x L = < < .
Solucin.
xx x 1 2
F (x, t) 0, F (0, t) (t), F(L, t) (t) = = = . De modo que
xx 1 2
x 1 2
F (0, t) (t) g (t) = = . Al imponer la segunda condicin de borde se tiene
1 1 2 1 1 1 2 1
F(x, t) g (t) x (t) F(L, t) (t) g (t) L (t) g (t) (t) L (t) = + = = + = . Por lo tanto
se obtiene la funcin definida por
2 1
F(x, t) (t) (x L) (t) = + . Por otro lado w(x,t) es
una funcin incgnita auxiliar cuya determinacin se estudiar a continuacin.

Como
x 1
F (0, t) (t) = y
2
F(L, t) (t) = por construccin, resulta que
x
w (0, t) 0 = y
w(L, t) 0 = . Adems w(x, 0) u(x, 0) F(x, 0) (x) F(x, 0) = = en 0 x L < < , y
xx t t xx t t

.

t
convertido como

xx t
kw G(x, t) w , k 0, 0 x L, t 0 + = > < < > ,

x
w (0, t) 0, w(L, t) 0, t 0 = = > ,
w(x, 0) H(x), 0 x L = < < .

condicin inicial no nula, condicin de frontera nula y extremo aislado.

(2n 1) x
n
2L
n 0
w(x, t) w (t)cos( )
+
=
=
,
y se expresa
(2n 1) x
n
2L
n 0
G(x, t) G (t)cos( )
+
=
=
, con

L
(2n 1) x
n
2L
0
2
G (t) G(x, t)cos( )dx, n 0,...
L
+
= =
,
(2n 1) (2n 1) x 2
xx n
2L 2L
n 0
w (x, t) w (t)( ) cos( )
+ +
=
=
,
(2n 1) x '
t n
2L
n 0
w (x, t) w (t)cos( )
+
=
=
.

xx t
(2n 1) (2n 1) x (2n 1) x (2n 1) x 2 '
n n n
2L 2L 2L 2L
n 0 n 0 n 0
k w (t)( ) cos( ) G (t)cos( ) w (t)cos( )

+ + + +
= = =
+ =

.
(2n 1) x ' 2 ' 2 n n
n n n n n n
L 2L L
n 0
w (t) k( ) w (t) G (t) cos( ) 0 w (t) k( ) w (t) G (t)
+

=
(
+ = + =

.

Como
(2n 1) x
n
2L
n 0
w(x, t) w (t)cos( )
+
=
=
,
entonces de la condicin inicial dada en el problema se tiene
(2n 1) x
n
2L
n 0
H(x) w (0)cos( )
+
=
=

(serie de cosenos de Fourier)
n
w (0) como

L
(2n 1) x
n
2L
0
2
W (0) H(x)cos( )dx
L
+
=
.

L
(2n 1) x ' 2 n
n n n n
L 2L
0
2
w (t) k( ) w (t) G (t) , W (0) H(x)cos( )dx
L
+
+ = =

Mtodo 1.
Solucin general:
(2n 1) (2n 1) (2n 1) 2 2 2
k( ) t k( ) t k( ) t
2L 2L 2L
n n n
+ + +

(
= + = + ( (

(

.
Valor C:
n n
(2n 1) 2
k( ) t
2L
n n n n
w (t) e M (t) w (0) M (0)
+
= + (

.

t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de
n n
M (t) y M (0) .
6. Clculo de los
n
w (t).

Mtodo 2.
(2n 1) (2n 1) ' 2 2
n n n n n n n
2L 2L
+ +
+ = + = L L L .
(2n 1) 2
n n n
2L
n n n n
n
(2n 1) (2n 1) (2n 1) 2 2 2
2L 2L 2L
s k( ) w (s) G (s) w (0)
G (s) w (0) G (s) w (0)
w (s)
s k( ) s k( ) s k( )
+
+ + +
(
+ = +

+
= = +
+ + +
(2n 1) (2n 1) 2 2
k( ) t k( ) t
-1 -1 n
2L 2L
n n n n n
(2n 1) 2
2L
G (s)
w (t) {w (s)} w (0)e P (t) w (0)e
s k( )
+ +

+

= = + = +
`
+

)
L L .

t
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de { }
n n
G (t) y P (t) L .
6. Clculo de los
n
w (t).
Solucin final:
(2n 1) x
2 1 n
2L
n 0
u(x, t) F(x, t) w(x, t) (t) (x L) (t) w (t)cos( )
+
=
= + = + +
.

Ejemplo 67. Sea una varilla de aluminio de longitud L que inicialmente est a la
temperatura uniforme de 25 C
. Suponga que, en el instante t 0 = , el extremo x 0 =

se enfra hasta 0 C
, mientras que el extremo x L = se calienta hasta 60 C
y que, a
partir de entonces, los dos se mantienen a esas temperaturas. Demuestre que la
distribucin de temperatura u(x,t) viene dada como
2 2 2
n 0.86n t / L
n 1
60x 50 70 n x
u(x, t) ( 1) e sen
L n n L

=
| | | |
= + +
| |

\ \
.
Solucin.
Planteamiento del problema:

xx t
0.86u u , 0 x L , t 0
u(0, t) 0 , u(L, t) 60 , t 0
u(x, 0) 25 , 0 x L
= < < >
= = >
= < <
.
Sea u(x, t) v(x, t) w(x, t) = + . Se tiene entonces el siguiente problema:

xx xx t t
0.86w 0.86v w v , 0 x L , t 0
w(0, t) 0 v(0, t) , w(L, t) 60 v(L, t) , t 0
w(x, 0) 25 v(x, 0) , 0 x L
+ = + < < >
= = >
= < <
.
Exigencias sobre v(x,t): v'' 0, v(0) 0, v(L) 60 = = = .
60x
L
v(x) =
El nuevo problema sera:

xx t
60x
L
0.86w w , 0 x L , t 0
w(0, t) 0, w(L, t) 0 , t 0
w(x, 0) 25 , 0 x L
= < < >
= = >
= < <

cuya solucin viene dada por
2 2 2
0.86n t / L n x
n
L
n 1
w(x, t) c e sen( )

=
=
,

L
n x 2
n
L L
0
c f(x)sen( )dx
.
En este caso

L
60x n x 2
n
L L L
0
c (25 )sen( )dx
.
Al calcular
n
c se tiene

L L L
60x n x 50 n x 120 n x 2
n
L L L L L L L
0 0 0
c (25 )sen( )dx sen( )dx xsen( )dx

= =

n 50 50 120L 50 70L
n
n n n n n
c cos(n ) cos(n ) ( 1)

= + + = + .
Por tanto
2 2 2
n 0.86n t / L
n 1
60x 50 70 n x
u(x, t) ( 1) e sen
L n n L

=
| | | |
= + +
| |

\ \
.

Ejemplo 68. Dado el problema de valor de frontera

t xx
x
u ku hu , 0 x , t 0 , k 0 , h 0
u(0, t) u( , t) 0 , t 0,
u(x, 0) e , 0 x
= < < > > >
= =
= < <

haga
ht
u(x, t) e w(x, t)
= y obtenga su solucin.
Solucin.
ht ht ht ht x
t t xx xx
u e w , u he w e w , u e w , u(x, 0) w(x, 0) e ,

= = + = = =
ht ht
u(0, t) e w(0, t) 0 w(0, t) 0 , u(L, t) e w(L, t) 0 w(L, t) 0

= = = = = = .
En consecuencia:
ht ht ht ht
t xx t xx
he w e w ke w he w w kw

+ = = .
Problema de valor en la frontera equivalente:

t xx
x
w kw , 0 x , t 0 , k 0
w(0, t) w( , t) 0 , t 0,
w(x, 0) e , 0 x
= < < > >
= =
= < <

Sea w(x, t) X(x)T(t) = . Por tanto
xx
w X'' T = y
t
w XT ' = . De modo que kX'' T XT' = . Al
dividir por kXT 0 se obtiene
X'' T '
X kT
= .
Igualando a una constante de separacin:
X'' T '
K
X kT
= = ,
se tienen el par de ecuaciones ordinarias: X'' KX 0 , T ' kKT 0 = + = . Se estudian 3
casos:
Caso I:
2
K 0 , K > = .
Ecuaciones ordinarias:

2
2 x x 2 k t
X'' X 0 X(x) Ae Be , T ' k T 0 T(t) Ce

= = + + = = .
Condiciones de borde homogneas:
w(0, t) 0 X(0) 0 A B 0 A B
w( , t) 0 X( ) 0 Ae Be A(e e ) 0 A 0 B 0

= = + = =
= = + = = = =

Solo se obtiene solucin trivial.
Caso II: K 0 = .
Ecuaciones ordinarias: X'' 0 X(x) A Bx , T ' 0 T(t) C = = + = = .
w(0, t) 0 X(0) 0 A 0 , w( , t) 0 X( ) 0 B 0 B 0 = = = = = = = .
Caso III:
2
K 0 , K < = .
Ecuaciones ordinarias:

2
2 2 k t
X'' X 0 X(x) A cos( x) Bsen( x) , T ' k T 0 T(t) Ce

+ = = + = = .

n
w(0, y) 0 X(0) 0 A 0
w( , t) 0 X( ) 0 B sen( ) 0 n , n 1,2,3,...
= = =
= = = = =

Se tiene:

2 2
n kt n kt
n n n
n 1 n 1
w(x, t) B sen(nx)C e , w(x, t) D e sen(nx)

= =
= =

Condicin inicial:
x
n
n 1
w(x, 0) e D sen(nx)
=
= =
.
En consecuencia:

x x
n
2
0
0
2 2 n sen(nx)
D e sen(nx)dx . e cos(nx)
n
n 1
(
= =
(

+

2 2
2 n 2 n
n impar . (1 e ) n par . (1 e ).
n 1 n 1

= + =

+ +

Por tanto
n 1
n
2
2 n
D . (1 ( 1) e )
n 1
+
= +
+

Se obtiene:
2
n 1 n kt
2
n 1
2 n
w(x, t) (1 ( 1) e )e sen(nx)
n 1
+
=
= +

2
ht n 1 n kt
2
n 1
2 n
u(x, t) e (1 ( 1) e )e sen(nx)
n 1
+
=
= +
.

Ejemplo 69. La temperatura u(x,t) en una varilla no aislada de longitud L, con
condiciones de frontera no nulas, con flujo de calor en la superficie de la varilla en
estudio hacia el medio que la rodea cuya temperatura se mantiene constante, digamos
m
u , se determina a partir del problema de valores en la frontera

xx m t
ku h(u u ) u , k 0, h ctte, 0 x L, t 0 = > < < > ,

1 2
u(0, t) (t), u(L, t) (t), t 0 = = > ,
u(x, 0) (x), 0 x L = < < .
Solucin.
xx 1 2
xx 1 2
1 1
(t) (t)
2 1
1 2 2 1 2 2
L

= + = = + = .
x
1 2 1
L
F(x, t) (t) (t) (t) = + (

. Por otro
lado w(x,t) es una funcin incgnita auxiliar cuya determinacin se estudiar a
continuacin.

Como
1
F(0, t) (t) = y
2
F(L, t) (t) = por construccin, resulta que w(0, t) 0 = y
w(L, t) 0 = . Adems w(x, 0) u(x, 0) F(x, 0) (x) F(x, 0) = = en 0 x L < < , y
xx t t xx t t

.

m t
G(x, t) (h(F(x, t) u ) F (x, t)) = + , se tiene que el
problema original queda convertido como

xx t
kw G(x, t) hw w , k 0, 0 x L, t 0 + = > < < > ,
w(0, t) 0, w(L, t) 0, t 0 = = > ,
w(x, 0) H(x), 0 x L = < < .

condicin inicial no nula y condiciones de frontera nulas.

n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
,
y se expresa
n x
n
L
n 1
G(x, t) G (t)sen( )
=
=
,
con

L
n x
n
L
0
2
G (t) G(x, t)sen( )dx, n 1,...
L

= =
,
entonces
2 n n x
xx n
L L
n 1
w (x, t) w (t)( ) sen( )

=
=
,
' n x
t n
L
n 1
=
=
.
xx t
kw G(x, t) hw w + = , se tiene
2 ' n n x n x n x n x
n n n n
L L L L L
n 1 n 1 n 1 n 1
k w (t)( ) sen( ) G (t)sen( ) h w (t)sen( ) w (t)sen( )

= = = =
+ =

.
' 2 ' 2 n n x n
n n n n n n
L L L
n 1
w (t) (k( ) h)w (t) G (t) sen( ) 0 w (t) (k( ) h)w (t) G (t)

=
(
+ + = + + =

.

Como
n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
,
n x
n
L
n 1
H(x) w (0)sen( )
=
=

n
w (0) como

L
n x
n
L
0
2
W (0) H(x)sen( )dx
L

=
.

L
' 2 n n x
n n n n
L L
0
2
w (t) (k( ) h)w (t) G (t) , W (0) H(x)sen( )dx
L

+ + = =

Mtodo 1.
Solucin general:
2 2 2 n n n
(k( ) h)t (k( ) h)t (k( ) h)t
L L L
n n n

+ + + (
= + = + (
(

.
Valor C:
n n
2 n
(k( ) h)t
L
n n n n
w (t) e M (t) w (0) M (0)
+
= + (

.
m t
G(x, t) (h(F(x, t) u ) F (x, t)) = +
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de
n n
M (t) y M (0) .
6. Clculo de los
n
w (t).

Mtodo 2.
' 2 2 n n
n n n n n n n
L L
{w (t)} (k( ) h) {w (t)} {G (t)} sw (s) w (0) (k( ) h)w (s) G (s)

+ = + + = L + L L .
2 n
n n n
L
n n n n
n
2 2 2 n n n
L L L
s k( ) h w (s) G (s) w (0)
G (s) w (0) G (s) w (0)
w (s)
s k( ) h s k( ) h s k( ) h

(
+ + = +

+
= = +
+ + + + + +
2 2 n n
(k( ) h)t (k( ) h)t
-1 -1 n
L L
n n n n n
2 n
L
G (s)
w (t) {w (s)} w (0)e P (t) w (0)e
s k( ) h

+ +

= = + = +
`
+ +

)
L L .

m t
G(x, t) (h(F(x, t) u ) F (x, t)) = +
2. Clculo de los
n
G (t) .
4. Clculo de
n
w (0)
5. Clculo de { }
n n
G (t) y P (t) L .
6. Clculo de los
n
w (t).

Solucin final:
n x x
1 2 1 n
L L
n 1
=
= + = + + (

.

Ejemplo 70. Encuentre la solucin u(x,t) de la ecuacin de calor unidimensional que
satisface las condiciones de frontera dadas por

xx t
u 1 u , 0 x , t 0
u(0, t) t , u( , t) 0 , t 0
u(x, 0) x , 0 x
+ = < < >
= = >
= < <
.
Solucin.
Paso 1. Construccin de G(x,t).
x x
F(x, t) t t t(1 )

= = .
x x
G(x, t) 1 1

= + = .
Paso 2. Clculo de los
n
G (t) .

n 1 2 2 x 1 2
n
2 2 n 2 n
n 0 0
G (t) xsen(nx)dx cos(nx) sen(nx) ( 1)

+

(
= = + =
(

.
Paso 3. Construccin de H(x). H(x) x 0 x = = .
Paso 4. Clculo de
n
w (0).

n 1 2 2
n
n
0
w (0) xsen(nx)dx ( 1)
= =
.
Paso 5. Clculo de
n n
M (t) y M (0) . (Mtodo 1)
2 2
n 1 n t n 1 n t 2 2
n
n 3
n
M (t) ( 1) e dt ( 1) e
+ +
= =
,
n 1 2
n
3
n
M (0) ( 1)
+
= .
n
w (t). (Mtodo 1)
2 2 2 2
n t n 1 n t n 1 n 1 n 1 n t n t 2 2 2 2 2
n
3 n 3 3
n n n
w (t) e ( 1) e ( 1) ( 1) ( 1) e e n 1
+ + + +

( (
= + = +
( (

.
Paso 5. Clculo de { }
n n
G (t) y P (t) L . (Mtodo 2)
{ }
{ }
n 1 n 1 2 2 1
n n
n n s
G (s) G (t) ( 1) ( 1)
+ +

= = L =L .
n 1
2 1
-1 n 1 -1 n 1 -1 n s 2 2
n
n n 2 4 4 2 4
2 2
n n n n
4 4 2 4
2
2
n
n
t 2
n 1 -1 n 1 2 n t 4 4 2
n
3 3 2 2 4
n n n n
2 4
( 1)
1 1
P (t) ( 1) ( 1)
s n
s n s (s )
P (t) ( 1) ( 1) e senh( )
(s )
+
+ +

+ +

= = =
` ` `
+
+ + +
)
) )

= =
`
+
)
L L L
L

n
w (t). (Mtodo 2)
2 2
2 2 n t n t
n t n
2 2 t 2
2 2 n 1 n 1 n t n 1 n 1 n t n t e e 4 2 4 2 2 2
n
3 2 n 3 2 n
n n
2 2 2 2
n 1 n 1 n t n 1 n t n 1 n t n t 2 2 2 2 2
3 3 n 3
n n n
w (t) ( 1) e senh( ) ( 1) e ( 1) e ( 1) e
( 1) ( 1) e ( 1) e ( 1) e e n 1

+ + + +

+ + + +

(
(
= + = +
(
(

(
= + = +
(

Paso 7. Construccin de w(x,t).
n 1 2 2
( 1) n t n t 2 2
3
n
n 1
w(x, t) e e n 1 sen(nx)
=
(
= +
(

.
Paso 8. Construccin de u(x,t).
n 1 2 2
( 1) n t n t 2 x 2
3
n
n 1
u(x, t) t(1 ) e e n 1 sen(nx)
+

=
(
= + +
(

.

Ejemplo 71. Si hay flujo de calor en la superficie de la varilla en estudio hacia el
medio que la rodea cuya temperatura se mantiene constante igual a cero, la ecuacin
de calor se convierte en la expresin
xx t
ku hu u = . Considere la temperatura u(x,t) en
una varilla delgada con condiciones u(0, t) u(L, t) 0 = = y
3 x
L
u(x, 0) sen( )
= . Halle la
solucin de la ecuacin sujeta a las condiciones.
Solucin.
Paso 1. Construccin de G(x,t). F(x, t) 0 = . G(x, t) 0 = .
n
G (t) .
n
G (t) 0 = .
Paso 3. Construccin de H(x).
3 x
L
H(x) sen( )
= .
Paso 4. Clculo de
n
w (0).
L L
(3 n) x (3 n) x
3 x n x 2 1
n
L L L L L L
0 0
L
(3 n) x (3 n) x
1 L L
L 3 n L 3 n L
0
W (0) sen( )sen( )dx (cos( ) cos( ))dx
( sen( ) sen( ) 0 , n 3
+

+
+
= =
= =

L L
L
2 3 x 6 x 6 x 2 1 1 L
3
L L L L L 6 L
0
0 0
W (0) sen ( )dx (1 cos( ))dx (x sen( )) 1

= = = =

.
Paso 5. Clculo de
n n
M (t) y M (0) . (Mtodo 1)
n n
M (t) M (0) 0 = = .
n
w (t). (Mtodo 1)
2 3
(k( ) h)t
L
3
w (t) e

+
= .
Paso 5. Clculo de { }
n n
G (t) y P (t) L . (Mtodo 2) { } L =
n n
G (t) P (t) 0 = .
n
w (t). (Mtodo 2)
2 3
(k( ) h)t
L
3
w (t) e

+
= .
Paso 7. Construccin de w(x,t).
2 3
(k( ) h)t
3 x L
L
w(x, t) e sen( )
= .
Paso 8. Construccin de u(x,t).
2 3
(k( ) h)t
3 x L
L
u(x, t) e sen( )
= .

4.27. ECUACIN DE ONDA UNIDIMENSIONAL

Hasta el momento, en el estudio de ecuaciones en derivadas parciales, el nico
problema fsico que se ha introducido es la conduccin del calor. Para ampliar la
discusin se estudiar ahora las vibraciones de cuerdas y membranas perfectamente
elsticas. Se comenzar formulando las ecuaciones que rigen las vibraciones de una
cuerda a partir de principios fsicos y se demostrar que las condiciones de contorno
adecuadas son similares en sentido matemtico a las condiciones de contorno para la
ecuacin de calor. Finalmente, se resolvern algunos ejemplos por el mtodo de
separacin de variables.

4.28. DEDUCCIN DE LA ECUACIN DE ONDA PARA LA
CUERDA VIBRANTE

Considere una cuerda de longitud L tal como una cuerda de guitarra tensa ante
dos puntos sobre el eje x, por ejemplo x 0 = y x L = . Cuando comienza a vibrar,
suponga que el movimiento se lleva a cabo en el plano xu, de tal modo que cada punto
de la cuerda se mueve en direccin perpendicular al eje x (vibraciones transversales).
De modo que u(x,t) representa el desplazamiento vertical de cualquier punto de la
cuerda, medido a partir del eje x, cuando t 0 > . Adems se supone que:

x x + x
( ) +
2
T x x,t
( )
1
T x,t
+ =
1 1 2
1
x

Figura 20. Vibraciones transversales en una cuerda delgada

La cuerda es perfectamente flexible.
La cuerda es homognea, esto es, , su masa por unidad de longitud es
constante.
Los desplazamientos u son pequeos comparados con el largo de la cuerda.
La pendiente de la curva es pequea en todos sus puntos.
La tensin T acta tangente a la cuerda y su magnitud T es la misma en todos sus
puntos.
La tensin es grande en comparacin con la fuerza de gravedad.
No hay otras fuerzas externas actuando sobre la cuerda.

En la figura 20, las tensiones
1
T y
2
T son tangentes a los extremos de la curva
en el intervalo x, x x + (

. Para
1
y
2
pequeos, la fuerza vertical neta que acta
sobre el elemento s correspondiente a la cuerda es, por consiguiente,
2 1 2 1 x x
Tsen Tsen Ttg Ttg T u (x x, t) u (x, t) = = + (

,

en donde
1 2
T T T = = . Ahora, s x es la masa de la cuerda en x, x x + (

y al
aplicar la segunda ley de Newton se obtiene
x x tt
T u (x x, t) u (x, t) xu + = (

o sea
x x
tt
u (x x, t) u (x, t)
u
x T
+
=
.

Si se toma el lmite cuando x 0 , esta ltima ecuacin se transforma en
xx tt
u ( / T)u = . Si se toma
2
a / T = , se obtiene
2
xx tt
u a u = .

4.29. PROBLEMAS ILUSTRATIVOS DE LA ECUACIN DE ONDA
UNIDIMENSIONAL

Ejemplo 72. Problema tpico de condicin en la frontera es

2 2
2
2 2
u u
a , 0 x L, a 0, t 0
x t

= < < > >

, (52)
u(0, t) 0, u(L, t) 0, t 0 = = , (53)

t 0
u
u(x, 0) f(x), g(x), 0 x L
t
=
= = < <
. (54)
Solucin.
Las condiciones de frontera (53) simplemente dicen que los extremos de la cuerda
permanecen fijos en todo instante. En t 0 = , las funciones f y g dadas en (54)
especifican la configuracin inicial y la velocidad inicial de cada punto de la cuerda,
respectivamente. En el contexto est implcito que f es continua y que f(0) 0 = y
f(L) 0 = . Sea u(x, t) X(x)T(t) = . La ecuacin se convierte en
2
a X'' T T '' X = . De modo
que

2
2 2
a X'' T T '' X
, XT 0
a XT a TX
= ,
2
2
X'' T ''
X
a T
= = .
Se tienen
2
X'' X 0 X Acos( x) Bsen( x) + = = + ,
2 2
T '' a T 0 T Ccos(a t) Dsen(a t) + = = + .
Condiciones de borde: u(0, t) 0 X(0)T(t) 0 X(0) 0. = = =
Se tiene: X(0) A cos(0) Bsen(0) 0 A 0 = + = = , u(L, t) 0 X(L)T(t) 0 X(L) 0. = = =
Entonces:
n
L
X(L) Bsen( L) 0 sen( L) 0 , n 1,2,3,...
= = = = = . Por lo tanto:
n n n n
an t an t n x
u (x, t) C cos D sen B sen
L L L
| | | | | | | |
= +
| | | |
\ \ \ \
.
n n
n 1
an t an t n x
u(x, t) E cos F sen sen
L L L
=
| | | | | | | |
= +
| | | |
\ \ \ \
.
Usando condiciones iniciales:

L
n n
0
n 1
n x 2 n x
u(x, 0) f(x) E sen E f(x)sen dx
L L L
=
| | | |
= = =
| |
\ \

t n n
n 1
an an t an an t n x
u (x, t) E sen F cos sen
L L L L L
=
| | | | | | | |
= +
| | | |
\ \ \ \

L
t n n
0
n 1
an n x an 2 n x
u (x, 0) g(x) F sen F g(x)sen dx
L L L L L
=
| | | |
= = =
| |
\ \

L
t n n
0
n 1
an n x 2 n x
u (x, 0) g(x) F sen F g(x)sen dx
L L an L
=
| | | |
= = =
| |
\ \

Observacin 9. Si la cuerda se suelta en reposo (g(x) 0) = , entonces
n
F 0 = y la
solucin se expresa como
n
n 1
an t n x
u(x, t) E cos sen
L L
=
| | | |
=
| |
\ \
.
Usando la identidad trigonomtrica
sen(A B) sen(A B)
sen(A)cos(B)
2
+ +
= ,
la solucin obtenida es

n
n 1
n n
n 1 n 1
1 n n
u(x, t) E sen (x at) sen (x at)
2 L L
1 n n
E sen (x at) E sen (x at)
2 L L
1
f(x at) f(x at) (Ecuacin de d' Alembert)
2
=

= =
( | | | |
= + +
( | |
\ \
(
| | | |
( = + +
| |
\ \ (

= + + (

Observacin 10.
n n
n 1
an t n x
u(x, t) G sen sen
L L
=
| | | | | |
= +
| | |
\ \ \
donde
2 2
n n n
G E F = + y

n n
n n
n n
E F
sen( ) , cos( )
G G
= = .

Ejemplo 73.

xx tt
t
4u u , 0 x 2 , t 0
u(0, t) 0 , u(2, t) 0, t 0
u (x, 0) 0 , u(x, 0) x , 0 x 2
= < < >
= =
= = < <

Solucin. Sea u(x, t) X(x).T(t) = . Por tanto:
xx
u X'' T = y
tt
u XT '' = . De modo que:
4X'' T XT '' = .
Al dividir por 4XT 0 se obtiene
X'' T ''
X 4T
= .
Igualando a una constante de separacin:
X'' T ''
K
X 4T
= = ,
se tienen el par de ecuaciones ordinarias: X'' KX 0 , T '' 4KT 0 = = . Se estudian 3
casos:
Caso I:
2
K 0 , K > = . Ecuaciones ordinarias:

2 x x 2 2 t 2 t
X'' X 0 X(x) Ae Be , T '' 4 T 0 T(t) Ce De

= = + = = +
2 2 2 2
u(0, t) 0 X(0) 0 A B 0 A B
u(2, t) 0 X(2) 0 Ae Be A(e e ) 0 A 0 B 0

= = + = =
= = + = = = =

Caso II: K 0 = . Ecuaciones ordinarias:
X'' 0 X(x) A Bx , T '' 0 T(t) C Dt = = + = = + .
u(0, t) 0 X(0) 0 A 0 , u(2, t) 0 X(2) 0 2B 0 B 0 = = = = = = = .
Caso III:
2
K 0 , K < = . Ecuaciones ordinarias:

2 2
X'' X 0 X(x) Acos( x) Bsen( x) , T '' 4 T 0 T(t) Ccos(2 t) Dsen(2 t) + = = + + = = +

n
n
2
t
u(0, t) 0 X(0) 0 A 0
u(2, t) 0 X(2) 0 B sen(2 ) 0 , n 1,2,3,...
T '(t) 2 ( Csen(2 t) Dcos(2 t)) u (x, 0) 0 T '(0) 0 D 0.
= = =
= = = = =
= + = = =

Se tiene:
n x n x
n n n
2 2
n 1 n 1
u(x, t) B sen( )C cos(n t) u(x, t) E sen( )cos(n t)

= =
= =

Condicin de borde no homognea:
n x
n
2
n 1
u(x, 0) f(x) E sen( )
=
= =
.
En consecuencia:
( ) ( ) ( )

2
2
n x n x n x
n
2 2 2
0
0
n 1
n
2 2
E xsen dx sen x cos
n n
2 4 2 4 4
n impar 2 . n par 2 . Por tanto : E ( 1) .
n n n n n

+
(
= =
(

= = = = = ( (

Se obtiene:
n 1
n x
2
n 1
4 ( 1)
u(x, t) cos(n t)sen( ) x
n
= =
.

Ejemplo 74. Problema tpico de condicin en la frontera es

2 2
2
2 2
u u
a , 0 x L, a 0, t 0
x t

= < < > >

, (55)

x x
u (0, t) 0, u (L, t) 0, t 0 = = , (56)

t 0
u
u(x, 0) f(x), g(x), 0 x L
t
=
= = < <
. (57)
Solucin.
En t 0 = , las funciones f y g dadas en (57) especifican la configuracin inicial y la
velocidad inicial de cada punto de la cuerda, respectivamente. En el contexto est
implcito que f es continua y que f(0) 0 = y f(L) 0 = . Sea u(x, t) X(x)T(t) = . La ecuacin
se convierte en
2
a X'' T T'' X = .
De modo que
2
2 2
a X'' T T '' X
, XT 0
a XT a TX
= ,
2
2
X'' T''
X
a T
= = . Se tienen
2
X'' X 0 X Acos( x) Bsen( x) + = = + ,
2 2
T '' a T 0 T Ccos(a t) Dsen(a t) + = = + .
Condiciones de borde:
x
u (0, t) 0 X'(0)T(t) 0 X'(0) 0 = = = ,
x
u (L, t) 0 X'(L)T(t) 0 X'(L) 0. = = =
Por lo tanto:
n n n n
an t an t n x
u (x, t) C cos D sen B cos
L L L
| | | | | | | |
= +
| | | |
\ \ \ \
, n 0,1,... =
0 n n
n 1
an t an t n x
u(x, t) E E cos F sen cos
L L L
=
| | | | | | | |
= + +
| | | |
\ \ \ \
.
Usando condiciones iniciales:

0
0 n n
n 1 n 1
L L
0 n
0 0
2E n x n x
u(x, 0) f(x) E E cos E cos
L 2 L
1 2 n x
donde E f(x)dx , E f(x)cos dx
L L L

= =
| | | |
= = + +
| |
\ \
| |
= =
|
\

t n n
n 1
an an t an an t n x
u (x, t) E sen F cos cos
L L L L L
=
| | | | | | | |
= +
| | | |
\ \ \ \

L
t n n
0
n 1
an n x an 2 n x
u (x, 0) g(x) F cos F g(x)cos dx
L L L L L
=
| | | |
= = =
| |
\ \

L
t n n
0
n 1
an n x 2 n x
u (x, 0) g(x) F sen F g(x)cos dx
L L an L
=
| | | |
= = =
| |
\ \

Ejemplo 75.

2
xx tt
a u f(x, t) u , a 0, 0 x L, t 0 + = > < < > ,

1 2
u(0, t) (t), u(L, t) (t), t 0 = = > ,

t
u(x, 0) (x), u (x, 0) (x), 0 x L = = < < .
Solucin.
xx 1 2
xx 1 2
1 1
(t) (t)
2 1
1 2 2 1 2 2
L

= + = = + = .
x
1 2 1
L
F(x, t) (t) (t) (t) = + (

. Por otro
lado w(x,t) es una funcin incgnita auxiliar cuya determinacin se estudiar a
continuacin.

Como
1
F(0, t) (t) = y
2
F(L, t) (t) = por construccin, resulta que w(0, t) 0 = y
w(L, t) 0 = . Por otro lado se tiene que
w(x, 0) u(x, 0) F(x, 0) (x) F(x, 0) = = ,
t t t t
w (x, 0) u (x, 0) F (x, 0) (x) F (x, 0) = =
en 0 x L < < , y
2 2
xx tt tt xx tt tt
a w f(x, t) w F (x, t) a w w F (x, t) f(x, t) + = + = + (

.
Si se denotan
H(x) (x) F(x, 0) = ,
t
I(x) (x) F (x, 0) = ,
tt
G(x, t) f(x, t) F (x, t) = ,
el problema original queda convertido como

2
xx tt
a w G(x, t) w , a 0, 0 x L, t 0 + = > < < > ,
w(0, t) 0, w(L, t) 0, t 0 = = > ,

t
w(x, 0) H(x), w (x, 0) I(x), 0 x L = = < < .

Este problema representa el caso de una ecuacin de onda no homognea con
condiciones iniciales no nulas y condiciones de frontera nulas.
n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
,
y se expresa
n x
n
L
n 1
G(x, t) G (t)sen( )
=
=
,
con

L
n x
n
L
0
2
G (t) G(x, t)sen( )dx, n 1,...
L

= =
,
entonces
2 n n x
xx n
L L
n 1
w (x, t) w (t)( ) sen( )

=
=
,
'' n x
tt n
L
n 1
=
=
.
2
xx tt
a w G(x, t) w , + = se tiene
2 2 '' n n x n x n x
n n n
L L L L
n 1 n 1 n 1
a w (t)( ) sen( ) G (t)sen( ) w (t)sen( )

= = =
+ =

.
Ecuacin lineal de segundo orden
'' 2 '' 2 an n x an
n n n n n n
L L L
n 1
w (t) ( ) w (t) G (t) sen( ) 0 w (t) ( ) w (t) G (t)

=
(
+ = + =

.
Como
n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
,
entonces de la primera condicin inicial del problema se tiene
n x
n
L
n 1
H(x) w (0)sen( )
=
=

n
w (0) como

L
n x
n
L
0
2
w (0) H(x)sen( )dx
L

=
.
De la segunda condicin inicial del problema se tiene
' n x
n
L
n 1
I(x) w (0)sen( )
=
=

'
n
w (0) como

L
' n x
n
L
0
2
w (0) I(x)sen( )dx
L

=
.
segundo orden sujeta a dos condiciones:

L L
'' 2 ' an n x n x
n n n n n
L L L
0 0
2 2
w (t) ( ) w (t) G (t) , w (0) H(x)sen( )dx , w (0) I(x)sen( )dx
L L

+ = = =

Mtodo 1.
Solucin general:
an t an t an t an t an t an t L
n 1 2 n n
L L an L L L L
w (t) C cos( ) C sen( ) sen( ) G (t)cos( )dt cos( ) G (t)sen( )dt

(
= + +
(

an t an t an t an t L
n 1 2 n n
L L an L L
w (t) C cos( ) C sen( ) M (t)sen( ) N (t)cos( )

( = + +

.
n
1 n
N (0)L
C w (0)
an
= +
.
'
' n
2 n n
L.N (0) L
C w (0) M (0)
an an
| |
= + |
|

\
.
tt
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
3. Construccin de H(x) y I(x). H(x) (x) F(x, 0) =
t
I(x) (x) F (x, 0) =
4. Clculo de
n
w (0) y
'
n
w (0).
5. Clculo de
'
n n n n n
M (t), M (0), N (t), N (0) y N (0) .
6. Clculo de los
n
w (t).

Mtodo 2.
'' 2 2 ' 2 an an
n n n n n n n n
L L
{w (t)} ( ) {w (t)} {G (t)} s w (s) sw (0) w (0) ( ) w (s) G (s)

+ = + = L L L .
'
2 2 ' n n n an
n n n n n
L 2 2 an
L
n
2 an
L
G (s) sw (0) w (0)
s ( ) w (s) G (s) sw (0) w (0) w (s)
s ( )
G (s)
s (
+ +
(
+ = + + =

+
=
+

'
n n
2 2 2 2 2 an an
L L
sw (0) w (0)
.
) s ( ) s ( )

+ +
+ +

-1 -1 ' n an t an t L
n n n n
L an L 2 2 an
L
' an t an t L
n n n
L an L
G (s)
w (t) {w (s)} w (0)cos( ) w (0)sen( )
s ( )
P (t) w (0)cos( ) w (0)sen( ).

= = + +
`
+

)
= + +
L L

tt
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
t
I(x) (x) F (x, 0) =
4. Clculo de
n
w (0) y
'
n
w (0).
5. Clculo de { }
n n
G (t) y P (t) L .
6. Clculo de los
n
w (t).

Solucin final:
n x x
1 2 1 n
L L
n 1
=
= + = + + (

.

Ejemplo 76. Resuelva la ecuacin diferencial

2
xx t tt
1 2
t
a u 2u u u , a 0 , 0 x L , t 0
u(0, t) (t) , u(L, t) (t) , t 0
u(x, 0) (x) , u (x, 0) (x) , 0 x L
+ = > < < >
= = >
= = < <
.
Solucin.
que debe resolver la ecuacin con condiciones

xx 1 2
F (x, t) 0, F(0, t) (t), F(L, t) (t) = = = .
De modo que
xx 1 2
F (x, t) 0 F(x, t) g (t) xg (t) = = + . Al imponer la primera condicin de
borde se tiene
1 1
1 2
F(x, t) (t) xg (t) = +
(t) (t)
2 1
2 1 2 2
L
F(L, t) (t) (t) Lg (t) g (t)

= = + = . Por lo tanto
x
1 2 1
L
F(x, t) (t) (t) (t) = + (

. Por otro lado w(x,t) es una funcin auxiliar cuya
determinacin se estudiar a continuacin. Como
1
F(0, t) (t) = y
2
F(L, t) (t) = por
construccin, resulta que w(0, t) 0 = y w(L, t) 0 = . Por otro lado se tiene que
w(x, 0) u(x, 0) F(x, 0) (x) F(x, 0) = = ,
t t t t
w (x, 0) u (x, 0) F (x, 0) (x) F (x, 0) = = en
0 x L < < , y
2 2
xx t t tt tt xx t tt tt t
a w 2w 2F w F w F a w 2w w w F 2F F + + = + + = + + (

.

Si se denotan
H(x) (x) F(x, 0) = ,
t
I(x) (x) F (x, 0) = ,
tt t
G(x, t) F (x, t) 2F (x, t) F(x, t) = + (

,
el problema original queda convertido en

2
xx t tt
t
a w 2w w G(x, t) w , a 0 , 0 x L , t 0
w(0, t) 0 , w(L, t) 0 , t 0
w(x, 0) H(x) , w (x, 0) I(x) , 0 x L
+ + = > < < >
= = >
= = < <
.

Este problema representa el caso de una ecuacin de onda no homognea con
condiciones iniciales no nulas y condiciones de frontera nulas. Si se ensaya con una
solucin de la forma
n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
, y
n x
n
L
n 1
G(x, t) G (t)sen( )
=
=
, con

L
n x
n
L
0
2
G (t) G(x, t)sen( )dx, n 1,...
L

= =
,
entonces
2 n n x
xx n
L L
n 1
w (x, t) w (t)( ) sen( )

=
=
,
' n x
t n
L
n 1
=
=
,
'' n x
tt n
L
n 1
=
=
.
2
xx t tt
a w 2w w G(x, t) w + + = , se tiene
2 2 ' n n x n x n x n x
n n n n
L L L L L
n 1 n 1 n 1 n 1
'' n x
n
L
n 1
a w (t)( ) sen( ) 2 w (t)sen( ) w (t)sen( ) G (t)sen( )
w (t)sen( )

= = = =
=
+ +
=

'' ' 2 an n x
n n n n
L L
n 1
'' ' 2 an
n n n n
L
w (t) 2w (t) (( ) 1)w (t) G (t) sen( ) 0
w (t) 2w (t) (( ) 1)w (t) G (t)
(
+ + =

+ + =

Como
n x
n
L
n 1
w(x, t) w (t)sen( )
=
=
,
entonces de la primera condicin inicial del problema se tiene
n x
n
L
n 1
H(x) w (0)sen( )
=
=

n
w (0) como

L
n x
n
L
0
2
W (0) H(x)sen( )dx
L

=
.
De la segunda condicin inicial del problema se tiene
' n x
n
L
n 1
I(x) w (0)sen( )
=
=

'
n
w (0) como

L
' n x
n
L
0
2
W (0) I(x)sen( )dx
L

=
.

'' ' 2 an
n n n n
L
L L
' n x n x
n n
L L
0 0
w (t) 2w (t) (( ) 1)w (t) G (t) ,
2 2
W (0) H(x)sen( )dx , W (0) I(x)sen( )dx
L L

+ + =
= =

Mtodo 1. Solucin general:

t t an t an t
n 1 2
L L
t
2t 2t Le an t an t an t an t
n n
an L L L L
w (t) C e cos( ) C e sen( )
sen( ) G (t)e cos( )dt cos( ) G (t)e sen( )dt

= + +
(
(

t
t t an t an t Le an t an t
n 1 2 n n
L L an L L
w (t) C e cos( ) C e sen( ) sen( )M (t) cos( )N (t)

( = + +

n
1 n
N (0)L
C w (0)
an
= +
.
'
' n n
2 n n n
L(2N (0) N (0)) L
C w (0) w (0) M (0)
an an
| |
+
= |
|

\
.
tt
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
t
I(x) (x) F (x, 0) =
4. Clculo de
n
w (0) y
'
n
w (0).
5. Clculo de
'
n n n n n
M (t), M (0), N (t), N (0) y N (0) .
6. Clculo de los
n
w (t).

Mtodo 2.
'' ' 2 an
n n n n
L
2 ' 2 an
n n n n n n n
L
{w (t)} 2 {w (t)} (( ) 1) {w (t)} {G (t)}
s w (s) sw (0) w (0) 2sw (s) 2w (0) (( ) 1)w (s) G (s)
+ + =
+ + + =
L L L L

2 2 ' an
n n n n
L
' '
n n n n n n n
n
2 2 2 2 2 2 2 2 an an an an
L L L L
s 2s 1 ( ) w (s) G (s) (s 2)w (0) w (0)
G (s) (s 2)w (0) w (0) G (s) (s 1)w (0) w (0) w (0)
w (s) .
(s 1) ( ) (s 1) ( ) (s 1) ( ) (s 1) ( )

(
+ + = + +

+ +
= = + +
+ + + +
-1 -1 t ' t n an t an t L
n n n n n
L an L 2 2 an
L
t ' t an t an t L
n n n n
L an L
G (s)
w (t) {w (s)} w (0)e cos( ) (w (0) w (0))e sen( )
(s 1) ( )
P (t) w (0)e cos( ) (w (0) w (0))e sen( )

= = + +
`
+

)
= + +
L L

tt
G(x, t) f(x, t) F (x, t) =
2. Clculo de los
n
G (t) .
t
I(x) (x) F (x, 0) =
4. Clculo de
n
w (0) y
'
n
w (0).
5. Clculo de { }
n n
G (t) y P (t) L .
6. Clculo de los
n
w (t).

Solucin final:
n x x
1 2 1 n
L L
n 1
=
= + = + + (

.

4.30. ECUACIN DE LAPLACE

A continuacin, se aplicarn las tcnicas aprendidas anteriormente a un tipo
diferente de problemas. Se considera el problema estacionario del calor en un dominio
bidimensional, con condiciones no homogneas todas. En ese caso, la ecuacin en
derivadas parciales es lineal y homognea, pero las condiciones de contorno, aunque
lineales, son todas no homogneas. Si todas las condiciones de contorno no son
homogneas, no se podr aplicar el mtodo de separacin de variables a este
problema en su forma actual. Se podr superar esta dificultad mediante la
descomposicin en problemas donde el mtodo de separacin de variables, sea
apropiado.

4.31. PROBLEMAS DONDE INTERVIENE LA ECUACIN DE
LAPLACE EN COORDENADAS CARTESIANAS

Aunque no se presentar su deduccin, esta ecuacin en dos o tres dimensiones
surge en problemas independientes del tiempo que conciernen a potenciales como el
electrosttico, el gravitacional y la velocidad en mecnica de fluidos. Adems, una
solucin de la ecuacin de Laplace se puede interpretar como la distribucin de
temperatura en un estado estable. De manera que, una solucin u(x,y) podra
representar la temperatura que vara de un punto a otro, pero no con el tiempo, en
una placa rectangular. La ecuacin de Laplace, en dos y en tres dimensiones, se
abrevia en la forma
2
u 0 = , donde
2 2
2
2 2
u u
u
x y

= +

y
2 2 2
2
2 2 2
u u u
u
x y z

= + +

se llaman, respectivamente, laplaciano en dos dimensiones y laplaciano en tres
dimensiones, de una funcin u.

Ejemplo 77. Problema de Dirichlet para un rectngulo.

2 2
2 2
u u
0, 0 x a, 0 y b
x y

+ = < < < <

u(x, 0) 0, u(x,b) 0, 0 x a = = < <
u(0, y) 0, u(a, y) f(y), 0 y b = = .
Solucin.
Sea u(x, y) X(x)Y(y) = . La ecuacin se convierte en X'' Y Y '' X = . De modo que

X'' Y Y '' X
, XY 0
XY XY
= ,
2
X'' Y ''
X Y
= = .
Se tienen
2 x x
X'' X 0 X Ae Be

= = + ,
2
Y '' Y 0 Y Ccos( y) Dsen( y) + = = + .
Primeras condiciones de borde homogneas:
u(x, 0) 0 X(x)Y(0) 0 Y(0) 0. = = = u(x,b) 0 X(x)Y(b) 0 Y(b) 0. = = =
Se tiene: Y(0) Ccos(0) Dsen(0) 0 C 0 = + = =
Se tiene:
n
b
Y(b) Dsen( b) 0 sen( b) 0 , n 1,2,3,...
= = = = = .
Por lo tanto:
n x n x
b b
n n n n
n y
u (x, y) A e B e D sen
b

| |
| |
= +
|
|
\
\
.
Otra condicin de borde homognea: u(0, y) 0 X(0)Y(y) 0 X(0) 0. = = =
Se tiene:
0 0
X(0) Ae Be 0 A B = + = = . Por lo tanto:

n x n x n x n x
b b b b
n n n n n
n n
n y n y
u (x, y) A e A e D sen E e e sen
b b
n x n y n x n y
2E senh sen F senh sen
b b b b

| | | |
| | | |
=
| |
| |
\ \
\ \
| | | | | | | |

| | | |
\ \ \ \

n
n 1
n x n y
u(x, y) F senh sen
b b
=
| | | |
=
| |
\ \


b
n n
0
n 1
b
2
n
an
0 b.senh
b
an n y an 2 n y
U(a, y) f(y) F senh sen F senh f(y)sen dy
b b b b b
n y
F f(y)sen dy
b
=
| |
|
\
| | | | | | | |
= = =
| | | |
\ \ \ \
| |
=
|
\

Ejemplo 78. Placa rectangular con caras laterales aisladas.

2 2
2 2
u u
0, 0 x a, 0 y b
x y

+ = < < < <

u(x, 0) 0, u(x,b) f(x), 0 x a = = < <

x x
u (0, y) 0, u (a, y) 0, 0 y b = = .
Solucin.
Sea u(x, y) X(x)Y(y) = . La ecuacin se convierte en X'' Y Y '' X = . De modo que

X'' Y Y '' X
, XY 0
XY XY
= ,
2
X'' Y ''
X Y
= = .
Se tienen
2
X'' X 0 X Acos( x) Bsen( x) + = = + ,
2 y y
Y '' Y 0 Y Ce De

= = + .
Primeras condiciones de borde homogneas:
x
u (0, y) 0 X'(0)Y(y) 0 X'(0) 0. = = =
x
u (a, y) 0 X'(a)Y(y) 0 X'(a) 0. = = =
Se tiene: X'(0) A sen(0) B cos(0) 0 B 0 B 0 = + = = =
Se tiene:
n
a
X'(a) A sen( a) 0 0 .
= = = =
Si 0 = ,
0 0 0 0 0
u (x, y) (A B x)(C D y) = + +
ox
u (0, y) 0 X'(0)Y(y) 0 X'(0) 0 B 0 = = = =
ox
u (a, y) 0 X'(a)Y(y) 0 X'(a) 0 0 0 = = = =
Otra condicin de borde homognea:
0
u (x, 0) 0 X(x)Y(0) 0 Y(0) 0 C 0 = = = =
Por lo tanto:
0 0
u (x, y) E y = . Si
n
a
= ,
n y n y
a a
n n n n
n x
u (x, y) C e D e A cos
a

| |
| |
= + |
|
|
\
\

Otra condicin de borde homognea:
n n n
u (x, 0) 0 X(x)Y(0) 0 Y(0) 0 C D = = = =
Por lo tanto:
n y n y
a a
n n n
n x n y n x
u (x, y) e e E cos senh F cos
a a a

| |
| | | | | |
= |
| | |
|
\ \ \
\

0
n
n 1
F y n y n x
u(x, y) F senh cos
2 a a
=
| | | |
= +
| |
\ \


0
n
n 1
a a
0 0
0 0
a a
2
n n
bn
0 0 a.senh
a
F b bn n x
u(x,b) f(x) F senh cos donde
2 a a
2 2
F b f(x)dx F f(x)dx
a ab
bn 2 n x n x
F senh f(x)cos dx F f(x)cos dx
a a a a
=
| |
|
\
| | | |
= = +
| |
\ \
= =
| | | | | |
= =
| | |
\ \ \

Ejemplo 79.

2 2
2 2
u u
0, 0 x a, 0 y b
x y

+ = < < < <

u(x, 0) f(x), u(x,b) g(x), 0 x a = = < <
u(0, y) F(y), u(a, y) G(y), 0 y b = = < < .
Solucin.
PRINCIPIO DE SUPERPOSICIN. Un problema de Dirichlet para un rectngulo se
puede resolver con facilidad separando variables cuando se especifican condiciones
homogneas para dos fronteras paralelas, sin embargo, el mtodo de separacin de
variables no se aplica a un problema de Dirichlet cuando las condiciones en la frontera
en los cuatro lados del rectngulo son no homogneas. Para salvar esta dificultad se
descomponen el problema en dos problemas llamados PROBLEMA 1 y PROBLEMA 2.
Estos problemas vienen dados como

PROBLEMA 1

2 2
2 2
u u
0, 0 x a, 0 y b
x y

+ = < < < <

u(x, 0) f(x), u(x,b) g(x), 0 x a = = < <
u(0, y) 0, u(a, y) 0, 0 y b = = < < .

PROBLEMA 2

2 2
2 2
u u
0, 0 x a, 0 y b
x y

+ = < < < <

u(x, 0) 0, u(x,b) 0, 0 x a = = < <
u(0, y) F(y), u(a, y) G(y), 0 y b = = < < .

De modo que:
SOLUCIN PROBLEMA = SOLUCIN PROBLEMA 1 + SOLUCIN PROBLEMA 2.

Ejemplo 80. Ecuacin de Poisson.

xx yy
u u G(x, y), 0 x a, 0 y b + = < < < < ,
u(x, 0) 0, u(x,b) 0, 0 x a = = < < ,
u(0, y) f(y), u(a, y) g(y), 0 y b = = < < .
Solucin.
n y
n
b
n 1
u(x, y) u (x)sen( )
=
=
,
y se expresa
n y
n
b
n 1
G(x, y) G (x)sen( )
=
=
, con

L
n y
n
b
0
2
G (x) G(x, y)sen( )dy, n 1,...
L

= =
,
se tiene que
n y ''
xx n
b
n 1
u (x, y) u (x)sen( )
=
=
y
n y 2 n
yy n
L b
n 1
u (x, y) u (x)( ) sen( )
=
=
.

xx yy
u u G(x, y) + = se tiene
n y n y n y '' 2 n
n n n
b b b b
n 1 n 1 n 1
u (x)sen( ) u (x)( ) sen( ) G (x)sen( )

= = =
=

.
n y '' 2 '' 2 n n
n n n n n n
b b b
n 1
u (x) ( ) u (x) G (x) sen( ) 0 u (x) ( ) u (x) G (x)

=
(
= =

.
Como
n y
n
b
n 1
u(x, y) u (x)sen( )
=
=
,
entonces de la condicin de borde del problema se tiene
n y
n
b
n 1
f(y) u (0)sen( )
=
=

n
u (0) como

b
n y
n
b
0
2
u (0) f(y)sen( )dy
b

=
.
De la segunda condicin de borde del problema se tiene
n y
n
b
n 1
g(y) u (a)sen( )
=
=

n
u (a) como

b
n y
n
b
0
2
u (a) g(y)sen( )dy
b

=
.

b b
n y n y '' 2 n
n n n n n
b b b
0 0
2 2
u (x) ( ) u (x) G (x) , u (0) f(y)sen( )dy , u (a) g(y)sen( )dy
b b

+ = = =

Mtodo 1.
Solucin general:

n x n x
n 1 2
b b
b n x n x n x n x
n n
n b b b b
u (x) C cosh( ) C senh( )
senh( ) G (x)cosh( )dx cosh( ) G (x)senh( )dx

= + +
(
(

n x n x b n x n x
n 1 2 n n
b b n b b
u (x) C cosh( ) C senh( ) M (x)senh( ) N (x)cosh( )

( = + +

.

n
1 n
N (0)b
C u (0)
n
= +
.
n a b n a n a 1
2 n 1 n n
n a b n b b
senh( )
b
C u (a) C cosh( ) M (a)senh( ) N (a)cosh( )

( ( =

.

1. Clculo de los
n
G (x) .
2. Clculo de
n
u (0) y
n
u (a) .
3. Clculo de
n n n n n
M (x), M (a), N (x), N (0) y N (a).
4. Clculo de los
n
u (x) .
5. Construccin de u(x,y).

Mtodo 2.
'' 2 2 ' 2 n n
n n n n n n n n
b b
{u (x)} ( ) {u (x)} {G (x)} s u (s) su (0) u (0) ( ) u (s) G (s)

= = L L L .
2 2 ' n n n n
n n n n n
b 2 2 n
b
n n
2 2 n
b
G (s) su (0) u (a)
s ( ) u (s) G (s) su (0) u (0) u (s)
s ( )
G (s) su (0)
s ( ) s
+ +
(
= + + =

= +

'
n
2 2 2 2 n n
b b
u (0)
.
( ) s ( )

+

-1 -1 ' n n x b n x
n n n n
b n b 2 2 n
b
' n x b n x
n n n
b n b
G (s)
u (x) {u (s)} u (0)cosh( ) u (0)senh( )
s ( )
P (x) u (0)cosh( ) u (0)senh( ).

= = + +
`

)
= + +
L L

' ' n a b n a n n a
n n n n n n n n
b n b n a b
b.senh( )
b
u (a) P (a) u (0)cosh( ) u (0)senh( ) u (0) u (a) P (a) u (0)cosh( )

(
= + + =

1. Clculo de los
n
G (x) .
2. Clculo de
n
u (0) y
n
u (a) .
3. Clculo de { }
'
n n n
G (x) , P (x) y u (0) L .
4. Clculo de los
n
u (x) .
5. Construccin de u(x,y).

Solucin final:
n y
n
b
n 1
u(x, y) u (x)sen( )
=
=
.
Ejemplo 81. Placa semiinfinita.

2 2
2 2
u u
0, 0 x L, y 0
x y

+ = < < >

u(0, y) 0, u(L, y) 0, y 0 = = >

0
u(x, 0) u , 0 x L , u(x, y) M = < < < .
Solucin.

X'' Y ''
X'' Y XY '' ctte X'' ctteX 0 , Y '' ctteY 0
X Y
= = = = + = .

n x
n 1
L
X'' ctteX 0 , X(0) X(L) 0 X (x) C sen( ) n 1,...
= = = = =

n y n y
L L
n 2 3
Y (y) C e C e n 1,...

= + =
n y n y
n x L L
n n
L
n 1
u(x, y) A e B e sen( )

=
(
= + (
(

.

Si y + entonces
n y
L
e

+ por tanto para garantizar que la temperatura
est acotada se tiene que hacer
n
A 0 = . Usando la condicin de borde se tiene

n 1
L
0
0 n x n x
0 n n
L L
0
n 1
2u 1 ( 1)
2u
u B sen( ) B sen( )dx
L n
+

=
(
+

= = =

,
de modo que la solucin se expresa como
n 1 n y n y
1 ( 1)
2u
n x 0 n x L L
n
L n L
n 1 n 1
u(x, y) B e sen( ) e sen( )
+ (
+
(

= =
= =

.

4.32. PROBLEMAS DONDE INTERVIENE LA ECUACIN DE
LAPLACE EN COORDENADAS POLARES

Pregunta.
Al usar coordenadas polares, cmo se expresa la ecuacin de Laplace dada en
coordenadas rectangulares como
2 2
2 2
u u
0
x y

+ =

?
Respuesta.
Sigamos los siguientes pasos:
a. Teniendo en cuenta que x r cos = , y rsen = , pruebe que:
x
r x :
r
Se tiene que
x
r r.cos x
r
= =
.
y
r y :
r
Se tiene que
y
r r.sen y
r
= =
.
x
y :
Se tiene que
x
r.sen y
= =
.
y
x :
Se tiene que
y
r.cos x
= =
.
b. Emplee la regla de la cadena para funciones de dos variables y pruebe que:
u u u
r x y .
r x y

= +

Prueba.
Se tiene que
u u x u y u u
. . .cos .sen
r x r y r x y

= + = +

. De modo que:
u u u u u
r r.cos . r.sen . x y
r x y x y

= + = +

.
u u u
y x
x y

= +

.
Se tiene que
u u x u y u u
. . r.sen r.cos
x y x y

= + = +

.
Entonces
u u u
y x
x y

= +

. Suponga que la funcin u satisface las condiciones
requeridas para que se aplique la regla de la cadena.
c. Pruebe que: (Diferenciando las ecuaciones en b)
2 2 2
2
u u u u u u
r x y (cos ) (sen )
r x r y r x y
r

+ = + + +

.
Demostracin evidente.
2 2 2
2
u u u u u
y x x y
x y x y

= +

.
d. Emplee c. para obtener
2 2 2 2 2 2
2
2 2
u u u u u u u
r r xr yr y x
r x r y r x y
r

+ + = + +

.
e. Pruebe que: (Idea: sustituya
u
x
y
u
y
para u en b.)
2 2 2
2
u u u
r x y
x r x y
x

= +

. Demostracin evidente.
2 2 2
2
u u u
r x y
y r y x
y

= +

2 2 2
2
u u u
y x
x x y
x

= +

2 2 2
2
u u u
y x
y y x
y

= +

f. Emplee d. y e. para establecer
2 2 2 2
2 2
2 2 2 2
u u u u u
r r r
r
r x y
| |

+ + = + |
|

\
.
Usando d. se tiene
2 2 2 2 2 2
2
2 2
u u u u u u u
r r xr yr y x
r x r y r x y
r

+ + = + +

.
Usando e. se tiene

2 2 2 2 2 2 2 2
2 2 2 2
2 2 2 2 2
2 2
2
2
u u u u u u u u u
r r x xy y xy y xy
r x y y x x y
r x y x
U U
x xy
y x
y

+ + = + + + + +

En conclusin, en coordenadas polares, la ecuacin de Laplace toma la forma:
2 2
2 2 2
u 1 u 1 u
. . 0
r r
r r

+ + =

.

Ejemplo 82. Problema de Dirichlet para un crculo. Considere el problema de
resolver la ecuacin de Laplace en la regin circular r a < , sujeta a la condicin de
frontera u(a, ) f( ), = en donde f es una funcin dada sobre 0 2 < .
Solucin.
Para completar el planteamiento del problema se observa que, para que u(r, ) sea de
valores nicos, es necesario que u sea peridica en con perodo 2 . Adems, se
expresa explcitamente que u(r, ) debe estar acotada para r a , ya que esto ser
importante posteriormente. Para aplicar el mtodo se separacin de variables a este
problema se supone que u(r, ) R(r) ( ) = . Al sustituir en la ecuacin diferencial se
obtiene
2
1 1
R''(r) ( ) .R'(r) ( ) .R(r) ''( ) 0
r
r
+ + = ,
o bien,
2
R'' R' ''
r r K
R R

+ = =
,
en donde K es la constante de separacin. Por tanto se obtienen las dos ecuaciones
diferenciales ordinarias
2
r R'' rR' KR 0 , '' K 0 + = + = .

En este problema no hay condiciones homogneas en la frontera; sin embargo,
recuerde que las soluciones deben ser acotadas y tambin peridicas en , con
perodo 2 . Se considerarn sucesivamente los casos en que K es negativa, cero y
positiva. Si K 0, < sea
2
K = , en donde 0 > ; entonces las ecuaciones ordinarias a
resolver son:
2 2 2
r R'' rR' R 0 , '' 0 + + = = . En particular la ecuacin
2
'' 0 =
tiene como solucin ( ) Ae Be

= + , la cual es peridica solo si A B 0 = = , de modo
que K 0, < no ofrece una solucin de inters. Si K 0 = , entonces las ecuaciones
ordinarias a resolver son
2
r R'' rR' 0 , '' 0 + = = . La ecuacin '' 0 = tiene como
solucin ( ) A B = + , la cual es peridica si B 0 = , de modo que ( ) es constante.
Por otro lado la ecuacin
2
r R'' rR 0 + = es de Cauchy-Euler y tiene la solucin
R(r) C Dln(r) = + . No es posible aceptar el trmino logartmico, si U(r, ) ha de
permanecer acotada cuando r 0 ; de donde D 0 = . Por tanto, correspondiente a
K 0 = , se obtiene la solucin
0 0
U (r, ) D = .

Por ltimo, si K 0 > , se hace
2
K = , en donde 0 > . Entonces las ecuaciones
ordinarias a resolver vienen dadas por
2 2 2
r R'' rR' R 0 , '' 0 + = + = . La ecuacin
2 2
r R'' rR' R 0 + = es una ecuacin de Cauchy-Euler y su solucin es R(r) Ar Br

= + ,
mientras que la ecuacin
2
'' 0 + = tiene solucin ( ) Ccos( ) Dsen( ) = + . Para
que ( ) sea peridica con perodo 2 es necesario que sea un entero positivo n.
Con n = se deduce que debe descartarse la solucin r
, en virtud de que se vuelve

no acotada cuando r 0 . Como consecuencia B 0 = y las soluciones apropiadas son

n n
n n n
u (r, ) E r cos(n ) F r sen(n ) , n 1,2,... = + = Estas funciones, junto con
0 0
u (r, ) D = ,
forman un conjunto de soluciones fundamentales del presente problema.

De la manera acostumbrada, suponga que u(r, ) puede expresarse como una
combinacin lineal de las soluciones fundamentales; es decir,
n 0
n n
n 1
E
u(r, ) r (E cos(n ) F sen(n ))
2
=
= + +
.

Entonces, la condicin de frontera requiere que
n 0
n n
n 1
E
u(a, ) f( ) a (E cos(n ) F sen(n ))
2
=
= = + +

donde,

2
0
0
1
E f( )d
2

=
,

2
n
n
0
1
E f( )cos(n )d
a

=
,

2
n
n
0
1
F f( )sen(n )d
a

=
.

Ejemplo 83. Problema de Dirichlet para un crculo.
2 2
2 2 2
U 1 U 1 U
. . 0
r r
r r
sen( ) 0
U(1, ) f( )
0 2

+ + =

= =

.
Solucin.
n 0
n n
n 1
E
U(r, ) r (E cos(n ) F sen(n ))
2
=
= + +

donde,

2
0
0
1
E f( )d
,

2
n
n
0
1
E f( )cos(n )d
a

=
,

2
n
n
0
1
F f( )sen(n )d
a

=
.
Para el caso de inters:

0
0
1 2
E sen( )d .
= =

n
2
0
0
n 1
n
2 2 2
1 1 n 1
E sen( )cos(n )d . seny.sen(ny) cos y cos(ny)
n
n 1
1 n 1 n 2 ( 1) 1
n impar . .0 n par . . E
n
n 1 n 1 (n 1)
+
(
= = +
(

=

n
2
0
0
1 1 n 1
F sen( )sen(n )d . cos y.sen(ny) cos ysen(ny) 0
n
n 1
(
= = + =
(

n 1
n
2
n 1
1 1 ( 1) 1
u(r, ) r cos(n )
n 1
+
=

= +

4.33. PROBLEMAS PROPUESTOS

1. Encuentre x y para cada uno de los siguientes pares de vectores en
3
R .
a. x y (1/ 2,2, 1) (4, 2,3) = =
b. x y (2 / 3,1/ 2,1) ( 1/ 2, 4,2) = =

2. Encuentre f g para cada uno de los siguientes pares de vectores en C 0,1 (

.
a. f g
2
(x) x (x) 1 x = =
b. f g (x) sen( x / 2) (x) cos( x / 2) = =
c. f g
1 1 1
2 2 2
(x) x (x) x = =

3. Encuentre f g para cada uno de los siguientes pares de vectores en C 0,1 (

si el
producto interior est definido con respecto a la funcin de peso r
x
(x) e = .
a. f g
x
(x) 1 2x (x) e
= =
b. f g (x) x (x) 1 x = =

4. Sea r cualquier funcin en C a,b (

que se anula para, cuando ms, un nmero
finito de valores en el intervalo a,b (

. Demuestre que

f g
b
a
f(x)g(x) r(x) dx =

define un producto interior en C a,b (

.

5. Encuentre la longitud de los siguientes vectores en
4
R .
a. (3, 4, 3,1)
b.
1
2
( , 3, 2,1)
c.
2 1 2 1
3 2 3 2
( , , , )

6. Encuentre la distancia entre cada uno de los siguientes pares de puntos en
4
R .
a. x y (3, 0, 1,5) (2,2, 1,3) = =
b. x y
3 3 2 1 1
3 6 2 2 2
( ,1, 0,2) ( , , , ) = =

7. Calcule f para cada uno de los siguientes vectores en C 0,1 (

.
a. f
2
(x) 1 x =
b. f(x) ln(x 1) = +
c. f
x / 2
(x) e =
d. f
2
(x) sen( x)
=

8. Encuentre el ngulo entre cada uno de los siguientes pares de vectores en
3
R .
a. x y
1 1
2 2
(1,1,1) ( , 1, ) = =
b. x y (1, 1, 0) (2, 1,2) = =
c. x y ( 3, 1, 0) (1,2, 5) = =
9. Encuentre el coseno del ngulo de cada uno de los siguientes pares de vectores en
3
P si el producto interior es

p q
1
1
p(x)q(x)dx
.
a. 1, x
b.
2
x, x
c. x,1 x

10. Repita el apartado anterior usando ahora el producto interior

p q
1
0
p(x)q(x)dx =
.

11. Sean p y q polinomios en el espacio vectorial
2
P (espacio de los polinomios de
coeficientes reales de grado a lo sumo dos). Discuta si la operacin definida de la
forma p, q p( 1)q( 1) p(0)q(0) p(1)q(1) = + + es un producto interno en
2
P .

12. Utilice la ley de los cosenos para establecer que
x y
x y
cos

=
para cualquier par de vectores no nulos x y y en
2
R .

13. Demuestre que x y x y + = + si, y solamente si, y x = o x y = para algn
nmero real 0 .

14. Sean x y y vectores arbitrarios en un espacio euclidiano, y se supone x y = .
Demuestre que x y + y x y son ortogonales. Interprete este resultado de forma
geomtrica.

15. Sea (V, , ) < > un espacio producto interno real y sean v, w elementos de V. Pruebe
que v w v w v w
2 2 2 2
2( ) + + = + . Esta igualdad se conoce como ley del
paralelogramo.

16. Demuestre que las funciones dadas son ortogonales en el intervalo indicado:
a.
2
1 2
f (x) x, f (x) x , 2 x 2 = =
b.
3 2
1 2
f (x) x , f (x) x 1, 1 x 1 = = +
c.
x x x
1 2
f (x) e , f (x) xe e , 0 x 2

= =
d.
2
1 2
f (x) cos x, f (x) sen x, 0 x = =

17. Demuestre que el conjunto de funciones dado es ortogonal en el intervalo indicado.
Encuentre la norma de cada funcin del conjunto.
a.
2
senx, sen3x, sen5x, ..., 0 x

b.
n x m x
p p
{1, cos( ), sen( )} n 1,2,3,... m 1,2,3,... p x p

= =

18. Verifique mediante integracin directa que las funciones son ortogonales con
respecto a la funcin de peso indicada, en el intervalo dado.
a.
2
2 x
0 1 2
H (x) 1, H (x) 2x, H (x) 4x 2 w(x) e , x
= = = = < <
b.
2 x 1
0 1 2
2
L (x) 1, L (x) x 1, L (x) x 2x 1 w(x) e , 0 x
= = + = + = < <

19. Sea
n
{ (x)} un conjunto ortogonal de funciones en a x b tal que
0
(x) 1 = y
1
(x) x. = Demuestre que

b
n
a
( x ) (x)dx 0 + =

para n 2,3,..., = y constantes y cualesquiera.

20. Sea
n
{ (x)} un conjunto ortogonal de funciones en a x b . Demuestre que
2 2 2
m n m n
(x) (x) (x) (x) + = + , m n .

21. Halle los valores de a, b y c para que f(x) ax 2 = + ,
2
g(x) bx cx 1 = + + y h(x) x 1 =
sean mutuamente ortogonales en el intervalo 0,1 (

.

22. Demuestre que las funciones
1
f (x) 1 = y
2
f (x) x = son ortogonales, con respecto al
producto interno integral en el intervalo 1,1 (

, y determine las constantes A y B
de modo que la funcin
2
3
f (x) 1 Ax Bx = + + sea ortogonal a
1
f y
2
f .

23. Sea
1 1
2 2
C , (

con el producto interno

1/ 2
1/ 2
f, g r(x)f(x)g(x)dx
,
donde la funcin peso r est definida mediante
2
1
r(x)
1 x
=
.
Pruebe que la sucesin { }
n
n 0
f

=
, donde
n
f (x) cos(n.arccos(x)) = es ortogonal.

24. Ortogonalice cada una de las bases siguientes en
3
R .
a. (1,1, 0), ( 1,1, 0), ( 1,1,1)
b.
1
3
(2, 1,1), (3, 2, ), (1, 0,1)
c. (1, 0, 0), (0,1, 0), ( 1, 0,1)

25. Ortogonalice cada uno de los conjuntos de vectores siguientes en C 0,1 (

.
a.
x x
e , e

b.
x 2x
e , e
c.
x
1, 2x, e

26. Determine cul de las siguientes sucesiones { }
k
a converge en R, y encuentre el
lmite de cada sucesin convergente.
a.
k
k
( 1)
a
k
=
b.
k
k
a 1 ( 1) = +
c.
k
k
a 1 e
= +
d.
k 1 k
k
k
2 ( 1)
a
2
+
+
=

27. Determine cul de las siguientes sucesiones { }
k
x converge en
3
R , y encuentre el
lmite de cada sucesin convergente. (En cada caso, los vectores
1 2
e e , y
3
e son
una base ortonormal para
3
R .
a.
k 1 2 3
x e e e
1 k 1
1
k k 2
+ | |
= + + +
|
+
\

b.
k 1 2 3
x e e e
k
senk
(1 2 ) (lnk)
k
= + +

28. Calcule los coeficientes de Fourier de la funcin
2
f(x) x , x = en C , (

con respecto al conjunto ortonormal
1 cos x cos2x
, , ,...,
2

y tambin con respecto al conjunto ortonormal
senx sen2x
, ,...

29. Repita el ejercicio anterior, para la funcin
x
e , x .

30. Considere que las funciones 1, senx, cos x,..., senkx, coskx,... son una base para el
espacio euclidiano C , (

. Encuentre los coeficientes de Fourier de f en C , (

en trminos de esta base. Cul es el desarrollo en serie de f?

31. Encuentre todas las soluciones de cada uno de los problemas siguientes con valor
en la frontera.
a. y'' y 0; y(0) 1, y( ) 1 + = = =
b. y'' y 0; y(0) 0, y'( ) 0 + = = =
c. y'' 4y sen(2x); y(0) 0, y( ) 0 + = = =
d.
x
y'' 5y e ; y(0) 0, y( ) 0 + = = =
e. y'' 9y 0; y'(0) 0, y'( ) 0 + = = =
f. y'' 9y cos(2x); y'(0) 0, y'( ) 0 + = = =
g. y'' 9y x; y'(0) 0, y'( ) 0 + = = =

32. Calcule los valores propios y las funciones propias para los problemas con valores
en la frontera en los siguientes ejercicios:
a. y'' (1 )y 0 ; y(0) 0, y( ) 0 + + = = =
b. y'' 2y' (1 )y 0 ; y(0) 0, y(1) 0 + + = = =
c. 4y'' 4y' (1 )y 0 ; y( 1) 0, y(1) 0 + + = = =
d. y'' 2y' (1 )y 0 ; y'(0) 0, y'( ) 0 + + = = =

33. La ecuacin diferencial de Hermite y'' 2xy' 2ny 0, n 0,1,2,... + = = tiene
soluciones polinomiales
n
H (x). Escriba la ecuacin en su forma autoadjunta y
deduzca una relacin de ortogonalidad.

34. Cul es la funcin propia correspondiente para el problema de valores en la
frontera y'' y 0, y'(0) 0, y( / 2) 0 + = = = para 25 = ?

35. La ecuacin diferencial de Laguerre xy'' (1 x)y' ny 0, n 0,1,2,... + + = = tiene
soluciones polinomiales
n
L (x) . Escriba la ecuacin en su forma autoadjunta y
deduzca una relacin de ortogonalidad.

36. Desarrolle la funcin g(x) sen(x), 0 x = , en trminos de las funciones propias
del problema de Sturm-Liouville
y'' y 0
y(0) 0, y( ) y'( ) 0
+ =
= + =
.

37. Determine los valores propios y las funciones propias del siguiente problema de
Sturm-Liouville

2
x y'' 3xy' y 0
y(1) 0, y(e) 0
+ + =
= =
.

38. Demuestre que las funciones
x x n x n x
1, cos , sen ,..., cos , sen ,...
a a a a
| | | | | | | |
| | | |
\ \ \ \

son linealmente independientes en PC a, a (

.

39. Descomponga cada una de las siguientes funciones en sus partes par e impar:
a.
x 1
x 1
+

b.
1
x 1 +

c. xcox(x) cos(2x)
d.
0
k k k k
k 1
a
f(x) (a cos(kx) b sen(kx)) a ,b cons tantes
2
=
= + +

40. Encuentre el desarrollo en serie de Fourier de cada funcin dada:
a. f(x) x, x = < <
b.
x
f(x) e , x = < <
c.
1 1
2 2
1 1
2 2
1 1
2 2
x x
f(x) 0 x
x x
<
<

d. f(x) ( x)( x), x = +
e.
x x 0
f(x)
x 0 x
+ < <
=

< <

41. a. Encuentre el desarrollo en serie de Fourier de la funcin
0 x 0
f(x)
x 0 x
<
=

<
.
b. Use esta serie para demostrar que
2
2 2 2
1 1 1
1 ...
8
3 5 7
= + + + +

42. Encuentre el desarrollo en serie de Fourier de la funcin f(x) cos(x) , 0 x = < < y
use este resultado para deducir que
2 2 2 2
2 1 3 5 7
...
16
2 1 6 1 10 1 14 1
= + +

.

43. Encuentre el desarrollo en serie de Fourier de senos de
x
e , siendo 0 x < < .

44. Sea f una funcin en PC 0, (

con la propiedad
f x f x
2 2
| | | |
+ =
| |
\ \
, 0 x 2 < < .
Demuestre que los nicos trminos no cero en el desarrollo en serie de Fourier de
coseno de f son de la forma
k
A cos(kx) , k impar donde

2
0
k
4
f(x)cos(kx)dx k 1,3,5,7,...
A
0 k 0,2, 4, 6,...
=
=

45. Encuentre el desarrollo en serie de Fourier de la funcin en PC 2,2 (

definida por
0 2 x 1
f(x) x 1 x 1
0 1 x 2
< <
= < <
< <
.

46. Encuentre el desarrollo en serie de Fourier de la funcin
1 8 x 9
f(x)
10 x 9 x 10
< <
=

< <
.

47. Sin recurrir a la integracin, encuentre el desarrollo en serie de Fourier, x ( , ) ,
de cada una de las siguientes funciones:
a.
3
f(x) sen (x) =
b.
2
x
f(x) sen(x)cos
2
| |
=
|
\

c.
2 2 2
f(x) sen (7x)sen (x)cos (x) =

48. Demuestre que
n
k 1
1
sen n x
2 1
cos(kx)
x 2
2sen
2
=
| | | |
+
| |
\ \
+ =
| |
|
\
,
como sigue:
Usando la frmula de Euler
ix
e cos(x) isen(x) = + se deduce que
n
k 1
1
cos(kx)
2
=
+

es la parte real de la expresin
n
ikx
k 1
1
e
2
=
+
.
Evale esta expresin usando la frmula para la suma de los n primeros trminos
de una serie geomtrica, y luego encuentre su parte real.

49. En cada uno de los ejercicios siguientes, encuentre la solucin u(x,t) de la ecuacin
de calor unidimensional que satisface las condiciones de frontera dadas:
a.
x 0 x L 2
u(0, t) u(L, t) 0; u(x, 0)
L x L 2 x L
<
= = =

<

b. u(0, t) u(L, t) 0; u(x, 0) x(L x) = = =
c.
x x
x
u (0, t) u (L, t) 0 u(x, 0) ksen
L
| |
= = =
|
\
, k ctte.
d.
x
u(0, t) 100, u(L, t) 0; u(x, 0) 100cos
2L
| |
= = =
|
\

e.
xx t
L
2
L
2
kU U
U(0, t) 0, U(L, t) 0
1 0 x
U(x, 0)
0 x L
= =
< <

=

< <

f.

xx t
2
u x u , 0 x , t 0
u(0, t) t , u( , t) 0 , t 0
u(x, 0) cos x , 0 x
+ = < < >
= = >
= < <

50. La temperatura en una barra aislada delgada, de longitud L, satisface las
condiciones en los puntos extremos u(0, t) 0 = y u(L, t) 1 = , y la condicin inicial
x
u(x, 0) sen
L
| |
=
|
\
.
a. Encuentre la temperatura en la barra como una funcin de posicin y tiempo.
b. Cul es la temperatura de estado estable en la barra (esto es, la temperatura
cuando t )?

51. Para una barra aislada de longitud 1, suponga que se mantiene la temperatura fija
con valor 0 en ambos puntos extremos y la distribucin de temperatura inicial se
da segn T(x) 200sen( x) 100sen(2 x) = . Suponga que la temperatura u(x,t)
satisface
t xx
u u = .
a. Encuentre u(x,t) para todo t 0 , 0 x 1 .
b. Dnde se encuentra el punto ms caliente de la barra cuando t 0 = ?

52. Suponga que la temperatura u(x,t) satisface
t xx
u u = .Suponga que para una barra
aislada de longitud 1, se mantiene la temperatura fija en 0 en ambos puntos
extremos, y la distribucin inicial se da segn la funcin T(x) 3sen(2 x) = .
a. Si la constante k en la ecuacin de calor tiene valor numrico 0.05, encuentre
la distribucin de temperatura u(x,1) en el tiempo t 1 = .
b. Si en el tiempo t 1 = la temperatura en cada punto ha decrecido a la mitad de
su valor inicial, cul es el valor numrico de k?

53. Resuelva el problema de valores en la frontera

x
xx t
ku Ae u , 0, 0 x 1, t 0
u(0, t) 0, u(1, t) 0, t 0, u(x, 0) f(x), 0 x 1
+ = > < < >
= = > = < <

La ecuacin diferencial es una forma de la ecuacin de conduccin de calor cuando
se genera calor dentro de una varilla delgada, por ejemplo, por decaimiento
radiactivo del material.

54. Resuelva el problema de valor de frontera

t xx
u 1 u , 0 x 1, t 0
u(0, t) 0, u(1, t) 0, u(x, 0) sen( x)
= + < < >
= = =
.

55. Conduccin de calor con radiacin. Una posibilidad interesante que puede
surgir es el caso donde uno de los extremos de la varilla, digamos x 0 = , se
mantiene a alguna temperatura, por ejemplo a 0 C
, mientras que el otro extremo

x L = irradia a un medio circundante, el cual se asume que est tambin a la
temperatura 0 C
. Se puede asumir como antes que la distribucin de temperatura

inicial est especificada por (x) . Resuelva la siguiente ecuacin diferencial que
modela la situacin antes descrita:

2
t xx
x x
u a u , 0 x L , t 0
u (0, t) u (L, t) hu(L, t) 0
u(x, 0) (x)
= < < >
= + =
.

56. Encuentre la temperatura de estado estable en una varilla que satisface la
condicin en la frontera
x
u (0, t) u(0, t) 0 = en el extremo x 0 = y se mantiene a la
temperatura constante T en el extremo x L = .

57. Resuelva la siguiente ecuacin diferencial:

t xx
x 1 2
u 10u 250u , 0 x 11 , t 0
u (0, t) K , u(11, t) K
u(x, 0) g(x)
= < < >
= =
.

58. Resuelva la ecuacin diferencial parcial de la difusin de calor en una barra
delgada de longitud L, cuya difusividad trmica es , aislada trmicamente en sus
dos extremos, si su temperatura inicial tiene la forma de un tringulo issceles con
valor mximo de temperatura
0
T 0 > en el punto medio de la barra y temperatura
nula en los dos extremos.

59. Sea una varilla de aluminio de longitud L que inicialmente est a la temperatura
uniforme de 25 C
. Suponga que, en el instante t 0 = , el extremo x 0 = se enfra

hasta 0 C
, mientras que el extremo x L = se calienta hasta 60 C
y que, a partir
de entonces, los dos se mantienen a esas temperaturas.
a. Encuentre la distribucin de temperatura en la varilla en cualquier instante t.
b. Use solamente el primer trmino de la serie para la temperatura u(x,t) con el
fin de hallar la temperatura aproximada en x 5 cm = , cuando t 30 s = y
cuando t 60 s = . (L 20 cm) =
c. Use los dos primeros trminos de la serie para u(x,t) a fin de encontrar un
valor aproximado de u(5,30) . Cul es la diferencia porcentual entre las
aproximaciones con uno y dos trminos? El tercer trmino tiene algn efecto
apreciable para este valor de t? (L 20 cm) =
d. Use el primer trmino de la serie para u(x,t) para estimar el intervalo de
tiempo que debe transcurrir antes de que la temperatura en x 5 cm = llegue a
ser menos de 1% de su valor de estado estable. (L 20 cm) =

60. Si hay flujo de calor en la superficie de la varilla en estudio hacia el medio que la
rodea cuya temperatura se mantiene constante, digamos
m
u la ecuacin de calor
se convierte en
xx m t
ku h(u u ) u = . Considere la temperatura u(x,t) en una varilla
delgada con u(0, t) u(L, t) 0 = = y u(x, 0) f(x) = . En vez de estar aislado
lateralmente, la varilla pierde calor hacia el ambiente (con temperatura fija igual a
cero) a una razn proporcional a u(x,t). Demuestre que la solucin de la ecuacin
xx t
ku hu u = sujeta a las condiciones anteriores viene dada por la expresin
2 2 2
ht n t / L
n
n 1
n x
u(x, t) e c e sen
L

=
| |
=
|
\
, donde

L
n
0
2 n x
c f(x)sen dx
L L
| |
=
|
\
.

61. Suponga que se pierde calor a travs de la superficie lateral de una varilla delgada
de longitud L, el cual pasa al medio ambiente que est a la temperatura cero. Si se
aplica la ley lineal de transferencia de calor, la ecuacin adopta la forma
xx t
ku hu u = , 0 x L < < , t 0 > , donde h es una constante. Halle la temperatura
u(x,t) si la temperatura inicial en toda la varilla es f(x) y los extremos x 0 = y
x L = estn aislados.

62. Resuelva el problema de valor en la frontera

2
t xx
x 1 2
u u u , 0 x L , t 0
u (0, t) K , u(L, t) K
u(x, 0) 0
= < < >
= =
.

63. Considere la ecuacin
xx t
au bu cu 0 + = , donde a, b y c son constantes.
a. Haga
t
u(x, t) e w(x, t)
= , donde es una constante, y encuentre la ecuacin

diferencial parcial correspondiente para w.
b. Si b 0 , demuestre que puede elegirse de modo que la ecuacin que se
encontr en la parte a) no tenga trmino en w.
c. Resuelva el problema

xx t
x
u u u 0 , 0 x , t 0
u (0, t) u( , t) 0 , t 0
u(x, 0) 0 , 0 x
+ = < < >
= = >
= < <
.

64. En cada uno de los ejercicios siguientes, encuentre la solucin u(x,t) de la ecuacin
de onda unidimensional que satisface las condiciones de frontera dadas:
a.
( )

2
xx tt
2hx L
L 2
t
x L
L 2
a U U
U(0, t) 0, U(L, t) 0
0 x
U(x, 0) , U (x, 0) 0
2h 1 x L
= =
< <

= =
< <

b.

2
tt xx
x
t
u a u , 0 x L , t 0
u(0, t) u (L, t) hu(L, t) 0 , t 0
u(x, 0) f(x) , u (x, 0) g(x)
= < < >
= + = >
= =

c.

2
tt xx
t
u a u f(x) , 0 x L , t 0
u(0, t) , u(L, t)
u(x, 0) 0 , u (x, 0) 0
= + < < >
= =
= =

d.

3
xx tt
t
u tsen x u , 0 x , t 0
u(0, t) 0 , u( , t) , t 0
u(x, 0) x , u (x, 0) 0 , 0 x
+ = < < >
= = >
= = < <

e.

x
xx tt
t
u te u , 0 x , t 0
u(0, t) 0 , u( , t) 0 , t 0
u(x, 0) 1 , u (x, 0) 0 , 0 x
+ = < < >
= = >
= = < <

65. Pruebe que la solucin u(x,t) del problema

2
tt xx
t
u a u , 0 x L , t 0
u(0, t) 0 , u(L, t) 0
u(x, 0) f(x) , u (x, 0) g(x)
= < < >
= =
= =

puede escribirse como u(x, t) v(x, t) w(x, t) = + , en donde v(x,t) es la solucin del
mismo problema con g(x) 0 = y w(x,t) es la solucin del mismo problema con
f(x) 0 = . De este modo, v(x,t) representa el movimiento de una cuerda a partir del
reposo con el desplazamiento inicial f(x) y w(x,t) representa el movimiento de una
cuerda que se pone en movimiento a partir de su posicin de equilibrio con
velocidad inicial g(x). Por lo tanto, se puede obtener la solucin del problema
general al resolver dos problemas subsidiarios.

66. Cuerda pulsada. Resuelva la siguiente ecuacin diferencial:

2
xx tt
L
2
t
L
2
a u u , 0 x L , t 0
u(0, t) u(L, t) 0 , t 0
bx 0 x
u(x, 0) , u (x, 0) 0
b(L x) x L
= < < >
= = >

= =

.

67. Cuerda vibrante con amortiguamiento pequeo. Resuelva la siguiente ecuacin
diferencial:

2
xx t tt
t
a u u u , 0 x L , t 0
u(0, t) u(L, t) 0 , t 0
u(x, 0) f(x) , u (x, 0) 0
= < < >
= = >
= =
.

68. Cuerda vibrante bajo la gravedad. Muestre que la solucin de la ecuacin
diferencial

2
xx tt
t
a u g u , 0 x L , t 0
u(0, t) u(L, t) 0 , t 0
u(x, 0) 0 , u (x, 0) 0
= < < >
= = >
= =

viene dada por
2
3 2 3 3 2
4gL 1 x at 1 3 x 3 at gx(L x)
u(x, t) sen cos sen cos ...
L L L L
a 1 3 2a
| | | | | | | | | |
= + +
| | | | |
\ \ \ \ \

69. Resuelva los siguientes problemas de valor en la frontera:
a.

2
xx tt
x
t
a u u , 0 x L , t 0
u(0, t) 0 , u (L, t) K , t 0
u(x, 0) 0 , u (x, 0) 0 , 0 x L
= < < >
= = >
= = < <

b.

xx tt
2
x
2 2
2 1 1
t
2 2
u 1 u , 0 x , t 0
u(0, t) 1 , u ( , t) , t 0
u(x, 0) 1 x x , u (x, 0) 0 , 0 x L
+ = < < >
= = >
= + = < <

70. Pruebe que una solucin del problema de valores en la frontera

xx tt
t
2
2
u u u , 0 x , t 0
u(0, t) u( , t) 0 , t 0 , u (x, 0) 0 , 0 x
x 0 x
u(x, 0)
x x
= + < < >
= = > = < <
< <

=

<

es
k 1
2
2
k 1
4 ( 1)
u(x, t) sen((2k 1)x)cos( (2k 1) 1 t)
(2k 1)
+
=
= +
.

71. Al someter una cuerda vibratoria a una fuerza vertical externa, que vara en
funcin de la distancia horizontal al extremo izquierdo, la ecuacin de onda tiene la
forma
2
xx tt
a u Ax u + = , donde A es una constante. Resuelva esta ecuacin
diferencial parcial sujeta a u(0, t) u(1, t) 0, t 0 = = > , u(x, 0) 0 = ,
t
u (x, 0) 0 = ,
0 x 1 < < .

72. Resuelva el siguiente problema de valor en la frontera:

2
xx tt
t
a u u , 0 x L , t 0
u(0, t) 0 , u(L, t) 0 , t 0
u(x, 0) x(L x) , u (x, 0) 0 , 0 x L
= < < >
= = >
= = < <

73. Resuelva el siguiente problema de valor en la frontera:

2
xx tt
n t n
n 1 n 1
a u u , 0 x L , t 0
u(0, t) 0 , u(L, t) 0 , t 0
n x n x
u(x, 0) A sen , u (x, 0) B sen , 0 x L
L L

= =
= < < >
= = >
| | | |
= = < <
| |
\ \

74. Resuelva

xx yy
y y
U U 0
U(0, y) 0, U(1, y) 1 y
U (x, 0) 0, U (x,1) 0
+ =
= =
= =
.

75. Resuelva

xx yy
x x
U U 0
U(x, 0) x, U(x,1) 1 x
U (0, y) 0, U (1, y) 0
+ =
= =
= =
.

76. Resuelva
1 2
rr r
2 1
2
U r U r U 0
U(1, ) sen (2 )

+ + =
= +
.

77. En cada uno de los siguientes ejercicios, encontrar la solucin de la ecuacin de
Laplace en el rectngulo 0 x L < < , 0 y M < < que satisface las condiciones de
frontera dadas:
a. u(0, y) u(L, y) 0, u(x, 0) u(x,M) 2x(x L) = = = =
b.
L
2
2
2 3L L L
2 2 2
0 0 x
u(0, y) u(L, y) 0, u(x, 0) u(x,M)
x x x L
<
= = = =

+ <

c.
y
M
u(x, 0) u(x,M) x(x L), u(0, y) u(L, y) sen( )
= = = =

78. Deduzca la frmula para la temperatura de estado estable en una placa cuadrada
de cuyas caras estn aisladas en forma tal que el calor no fluye a travs de
ellas, bajo la suposicin de que
x x
u (0, y) u ( , y) u(x, ) 0, u(x, 0) f(x) = = = = .

79. Resuelva la siguiente ecuacin diferencial:
xx yy
u u 0, 0 x a, 0 y b
u(0, y) 0, u(a, y) 0, 0 y b
u(x, 0) 0, u(x,b) f(x), 0 x a
+ = < < < <
= = < <
= = < <

.

80. Determine la temperatura de estado estable u(r, ) en un anillo semicircular, si
u(a, ) ( ), u(b, ) 0, 0
u(r, ) 0, u(r, ) 0, a r b
= = < <
= = < <
.

81. Determine la temperatura de estado estable u(r, ) en una placa semicircular de
radio r 2 = , si
0
u 0 2
u(2, )
0 2
< <
=

< <
,
con
0
u constante, y las orillas 0 = y = estn aisladas.

82. Determine la temperatura de estado estable u(r, ) en una placa circular de radio
a; si se sabe que
r
u (a, 0) f( ) = , 0 2 < < .

RESPUESTAS A LOS PROBLEMAS PROPUESTOS

1. a. -5 b. 11/3 3. a. 0 b. 3-e 5. a. 35 b. 5/2 c.
5
6
2
7. a.
2
15
30 b. 2(1 ln2) c. e 1 d.
1
2
2
9. a. 0 b. 0 c.
1
2
10.
1 1 1
2 4 2
3, 15,
15. a. (1,1,1), (-1,1,0), (0,0,1) b.
1 2
18 83
(2, 1,1), (4, 11, 19), (5,7, 3)
c. (1,0,0), (0,1,0), (0,0,1)
16. a.
x x x
2
2
e , e e
e 1
b.
2
x 2x x
2(e e 1)
e , e e
3(e 1)
+ +
+

c.
x
1, 2x 1, e 6(e 3)x 10 4e + +
17. a. 0 b. No converge c. 1 d. 2 18. a.
1 2 3
e e e + + b. No converge
20. Respecto al primer conjunto
0
e e
c
2

;
k
k
2
( 1)
c (e e ), k 1,2,3,...
(1 k )

= =
+

Respecto al segundo conjunto
k 1
k
2
( 1) k
c (e e ), k 1,2,3,...
(1 k )
+

= =
+

31. a. y cos x csenx = + c. No hay solucin e. y c.cos(3x) = g. No hay solucin
32. a.
2
n n
n 1 , (x) sen(nx), n 1,2,... ortogonal en C 0, = = = (

b.
2 2 x
n n
n , (x) e sen(n x), n 1,2,3,...;
= = = ortogonal en C 0,1 (

con respecto
a la funcin de peso
2x
e .
c.

x / 2
2 2
n n
x / 2
n x
e sen , n 2, 4, 6,...
2
n , (x)
n x
e cos , n 1,3,5,...
2
| |
=
|
\
= =

| |
=
|
\
ortogonal en C 1,1 (

respecto a
la funcin peso
x
e
.
d.
2 x
n n
n , (x) e (ncos(nx) sen(nx)), n 1,2,3,...
= = + =

39. a.
2
E 0
2 2
x 1 2x
f (x) , f (x)
x 1 x 1
+
= =

b.
E 0
2 2
1 x
f (x) , f (x)
1 x 1 x
= =

c.
E 0
f (x) cos(2x) , f (x) x cos(x) = =
d.
0
E k 0 k
k 1 k 1
a
f (x) a cos(kx) , f (x) b sen(kx)
2

= =
= + =

40. a.
k 1
k 1
( 1)
2 sen(kx)
k
+
=
b.
k
2
k 1
e e 1 ( 1)
(cos(kx) ksen(kx))
2
1 k

=
(

( +
+ (

c.
n 1
2
n 1
1 1 1 n 2 1 n 2 n
sen cos(nx) ( 1) cos sen sen(nx)
4 2 n 2 n n 2 2
n
+
=
( | | | | | |
+ + +
`
( | | |

\ \ \ )

d.
k 1
2
2
k 1
2 ( 1)
4 cos(kx)
3
k
+
=
e.
k 1
1
2 sen(kx)
k

41. a.
k k 1
2
k 1
( 1) 1 ( 1)
cos(kx) sen(kx)
4 k
k
+
=
| |

+ + |
|
5.
2
k 1
8 k
sen(2kx)
4k 1

43.
k 1
2
k 1
2 k
1 ( 1) e sen(kx)
1 k
+
=
(
+

45.
k
k 1
1 k x
f(x) A cos
4 2
=
| |
= +
|
\
, donde
2 2
2 2
k
2 2
2 4
k 1,5, 9,...
k
k
8
k 2, 6,10,...
A
k
2 4
k 3,7,11,...
k
k
0 k 4, 8,12,16,...
.
46.
k
2 2
k 1
3 ( 1) 1 1
cos(k x) sen(k x) , 8 x 10
4 k
k
=
(

+ + < <
(

47.a.
3 1
sen(x) sen(3x)
4 4
b.
1 1
sen(x) sen(2x)
2 4
+
c.
1 1 1 1 1
cos(4x) cos(10x) cos(14x) cos(18x)
16 16 32 16 32
+ +
49.a.
2 2 2 2
(2m 1) t / L a
2 2
m 0
4L 1 (2m 1) x
u(x, t) e sen
L
(2m 1)
(
+
(

=
+ | |
=
|
+ \

b.
2 2 2 2 2
(2m 1) t / L a
3 3
m 0
8L 1 (2m 1) x
u(x, t) e sen
L
(2m 1)
(
+
(

=
+ | |
=
|
+ \

c.
2 2 2 2
4m t / L a
2
m 1
2k 1 2m x
u(x, t) 1 2 e cos
L
4m 1
(

(

=

| |
=
`
|
\

)

d.
2 2 2 2
n t / L a
2
n 1
x 200 1 n x
u(x, t) 100 1 e sen
L L
n(4n 1)
(

(

=
| | | |
= +
| |
\ \

e.
n
2 2 2
k(n / L )t 2 n x
L
n 1
cos( ) 1
2
u(x, t) e sen( )
n

=
+
=

50.a.
2 2 2 2 n
n t / L a
n 1
x 2 ( 1) n x
u(x, t) e sen
L n L
(

(

=
| |
= +
|
b.
x
u(x, t)
L
=
51.a.
2 2
t 4 t
u(x, t) 200e sen( x) 100e sen(2 x)

= b. x 2 3 =
52.b.
2
(ln2) (4 )
53.
2 2
kn t
n
n 1
u(x, t) (x) A e sen(n x)

=
= +
donde

x
2
A
(x) e (e 1)x 1
k

(
= + +

y

L
n
0
A 2 f(x) (x) sen(n x)dx = (

54.
2 2
n t
3 3
n 1
1 2(1 cos(n ))
u(x, t) x(1 x) 1 e sen(n x)
2
n

=
(
= +
(

55.

2 2
L
2 2
a t n
n
n n
2 2
0
n n 1
h
u(x, t) 2 ( )cos( )d e cos( x)
L(h ) h

=
(
+
=
(
+ +
(

donde
n
son las raices positivas de tg( L) h = .
Nota. Pruebe que

2 2
L
2
2 n
n n
2 2
0
n
L(h ) h
sen( x) sen ( x)dx
2(h )
+ +
= =
+

56. v(x) T(1 x) /(1 L) = + +
59.a.
2 2 2
0.86n t / L
n 1
60x 50 70 n x
u(x, t) cos(n ) e sen
L n n L

=
| | | |
= + +
| |

\ \

b. u(5,30) 12.6 C ; u(5, 60) 13.7 C

c. u(5,30) 14.1 C ; 12% ; tercer tr mino 0.005 C

d. t 160 s
61.

L L
2 2 2
ht n t / L
0 0
n 1
1 2 n x n x
u(x, t) e f(x)dx f(x)cos dx e cos
L L L L

=
(
| | | | | |
( = +
| | |
\ \ ( \

64.a.
n
2 n at n x
L L 2 2
n 1
sen( )
8h
u(x, t) cos( )sen( )
n

=
=

b.
n n n n n
n 1
u(x, t) A cos(a t) B sen(a t) sen(a x)
=
= + (

donde

L L
n n n n
2 2
0 0
n n n
1 1
A f(x)sen( x)dx B g(x)sen( x)dx
sen( x) a sen( x)
= =

n
son las raices positivas de htg( L) = .
e.
n 1
n
2
n 1
2 ( 1) 1
M
u(x, t) cos(nt) t sen(nx)
n
n
+
=
( (
+

(
= +
(

66.
2 2
n 1
4bL 1 n n at n x
u(x, t) sen cos sen
2 L L
n
=
| | | | | |
=
| | |
\ \ \

69.a.
n 1
2KL 1 (2n 1) x (2n 1) at
u(x, t) Kx sen cos
2n 1 L L
=
| | | |
= +
| |

\ \

b.
2
3
n 1
x 4 (2n 1) cos(n ) 2
u(x, t) 1 sen((2n 1)x)cos((2n 1)t)
2
(2n 1)
=
(
+
=
(

71.
n
3
2 2 3 3
n 1
A 2A ( 1)
u(x, t) (x x ) cos(n at)sen(n x)
6a a n
= +

72.
2
n 1
3 3
n 1
4L 1 n at n x
u(x, t) 1 ( 1) cos sen
L L
n
+
=
| | | |
(
= +
| |

\ \

73.
n n
n 1
n at L n at n x
A cos B sen sen
L n a L L
=
( | | | | | |
+
( | | |
\ \ \

74.
n
2 2
n 1
1 2 1 ( 1)
U(x, y) x senh(n x)cos(n y)
2
n senh(n )
=

= +

75.

n n
n y n y
n n n n
2 2 n n 2 2 n n
n impar
1 4 4e 4e 4
U(x, y) (A e B e )cos(n x) A , B
2
n (e e ) n (e e )

+ +
= + + = =

76.
4 1
2
U(r, ) 1 r cos(4 ) =
77.a.

m
m 0
2
m
3 3
(2m 1) (M y) (2m 1) y (2m 1) x
u(x, y) C senh senh sen
L L L
16L
donde C
(2m 1) senh (2m 1) M/ L
=
( + + + | | | | | |
= +
( | | |
\ \ \
=
+ + (

b.

n
n 1
2
n
n
3 3
n (M y) n y n x
u(x, y) C senh senh sen
L L L
L n n
donde C n sen 4 ( 1) cos
2 2
n senh(n M/ L)
=
( | | | | | |
= +
( | | |
\ \ \
( | | | |
= +
`
( | |
\ \ )

c.
{ }

m
m 0
2
m
3 3
senh( x / M) senh (L x) /M sen( y / M)
u(x, y)
senh( L / M)
(2M 1) (M y) (2m 1) y (2m 1) x
C senh senh sen
L L L
8L
donde C
(2m 1) senh (2m 1) M/ L
=
+ (

= +
( + + + | | | | | |
+
( | | |
\ \ \
=
+ + (

78.

0
n n
0
n 1
a senh(n( y))cos(nx) 2
u(x, y) ( y) a donde a f(x)cos(nx)dx
2 senh(n )
=

= + =

79.
n
n 1
n y n x
u(x, y) A sen sen
a a
=
| | | |
=
| |
\ \
donde
( )

a
n
n b
0
a
2 n x
A f(x)sen dx
a
asenh

| |
=
|
\

80.
n
n 2n 2n
3 2n 2n
n 1
4 1 ( 1) r b a
u(r, ) . sen(n )
r
n a b
=
| |
=
|

81.
( )
n n
2 0 0
n 1
sen
u 2u r
u(r, ) cos(n )
2 n 2

=
| |
= +
|

Tema 4

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

Tema 4

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tema 4

Cargado por

Copyright:

Formatos disponibles

Prof.

Jos Luis Quintero 1

, y (2) adopta la forma dada por la expresin

= = . A continuacin se muestra la grfica (figura 2) para el

= . (Tambin observe, de la geometra del caso,

= . Sea la funcin de peso

(ver figuras 12, 13, 14 y 15)

= , respectivamente. As, una solucin

. Esta ecuacin modela la distribucin

. Normalmente se establece que k K / = , y

. Suponga que, en el instante t 0 = , el extremo x 0 =

, mientras que el extremo x L = se calienta hasta 60 C

+ = > < < >

+ + = > < < >

, en virtud de que se vuelve

+ = > < < >

= = > = < <

, mientras que el otro extremo

. Se puede asumir como antes que la distribucin de temperatura

= < < >

. Suponga que, en el instante t 0 = , el extremo x 0 = se enfra

, mientras que el extremo x L = se calienta hasta 60 C

= < < >

= , donde es una constante, y encuentre la ecuacin

= < < >

= + < < >

+ = < < >

+ = < < >

= < < >

= < < >

= < < >

= < < >

= < < >

+ = < < >

= + < < >

= = > = < <

= < < >

= < < >

También podría gustarte