Edps
Edps
Edps
Una ecuacion en derivadas parciales (EDP) es una ecuacion diferencial cuya incognita es una funcion
que depende de mas de una variable. El orden de una EDP es el orden de la derivada parcial mas alta.
En este tema vamos a estudiar algunas EDPs lineales de segundo orden a coeficientes constantes.
1Por ejemplo, un transistor se calienta cuando una corriente electrica circula por el.
1
2 Depositado en http://www.ma1.upc.edu/ed/edp.pdf
Equilibrios elasticos y termicos 1D del caso homogeneo. Buscamos los equilibrios elasticos en
ausencia de fuerzas externas y los equilibrios termicos en ausencia de focos/sumideros de calor internos:
F (x, t) 0. Equilibrio significa que el estado del cuerpo se mantiene estacionario en el tiempo, luego
buscamos soluciones u = u(x) que no dependan del tiempo y as desaparecen las derivadas parciales
ut y utt . En tal caso, las EDPs utt = c2 uxx y ut = k 2 uxx se reducen a la EDO lineal de segundo orden
u00 = 0, cuyas unicas soluciones son las funciones lineales de la forma u(x) = ax + b, con a, b R.
Queda probado pues que los unicos equilibrios elasticos de una cuerda vibrante no sometida a
fuerzas externas y los unicos equilibrios termicos de una barra sin focos ni sumideros de calor internos
son los estados (desplazamiento o temperatura) lineales.
Las versiones multidimensionales. Antes de dar las versiones multidimensionales de las ecuaciones
anteriores, recordamos que el Laplaciano de una funcion u : Rn R que depende de una variable
vectorial x = (x1 , . . . , xn ) Rn es la funcion
n
X 2u
u = div grad u = = ux1 x1 + + uxn xn .
j=1
x2j
Por ejemplo, si la funcion u depende de una unica variable x R, entonces u = uxx . En cambio, si
depende de tres variables: (x, y, z) R3 , entonces u = uxx + uyy + uzz . Ademas, cuando la funcion
dependa de la posicion x = (x1 , . . . , xn ) y el tiempo t, interpretaremos que el Laplaciano solo afecta
a las variables de posicion; es decir, el Laplaciano no incluye el termino utt .
Las versiones n-dimensionales de las ecuaciones anteriores son las siguientes.
La ecuacion de ondas que modela el movimiento ondulatorio de un cuerpo elastico Rn es
utt c2 u = F (x, t), u = u(x, t), x = (x1 , . . . , xn ) , t R.
La ecuacion del calor que modela la evolucion de la temperatura en un cuerpo Rn es
ut k 2 u = F (x, t), u = u(x, t), x = (x1 , . . . , xn ) , t > 0.
Desde un punto de vista fsico, solo interesan los casos 1D, 2D o 3D. Es decir, n 3. Al igual que en
las versiones 1D estamos suponiendo que el cuerpo es completamente homogeneo. La funcion F (x, t)
representa la accion de una fuerza exterior en el caso de la ecuacion de ondas o los focos y sumideros
de calor internos en el caso de la ecuacion de calor. Se supone que F (x, t) es un dato conocido.
En cambio, la ecuacion del calor es de primer orden en el tiempo, luego basta fijar la temperatura
inicial : u(x, 0) = f (x) para x . Y, para acabar, la ecuacion de Laplace/Poisson es estatica, luego
no tiene sentido fijar el estado inicial del objeto, ya que ese estado es justamente la incognita del
problema. Sera como preguntar de que color es el caballo blanco de Santiago.
Ejemplo 1. Sea W R3 un cuerpo de frontera S = W . Un PVI de calor en este cuerpo 3D sin focos
ni sumideros de calor interno con condiciones de frontera de tipo Neumann consiste en las ecuaciones
ut = k 2 u (x, y, z) W t > 0
u(x, y, z, 0) = f (x, y, z) (x, y, z) W
u
n (x, y, z, t) = h(x, y, z, t) (x, y, z) S t > 0
Ejemplo 2. Un PVI de calor 1D en una barra de longitud L sin focos ni sumideros de calor interno
con condiciones de frontera de tipo Neumann consiste en las ecuaciones
2
ut = k uxx
x (0, L) t > 0
u(x, 0) = f (x) x (0, L)
ux (0, t) = h i (t) t>0
ux (L, t) = hd (t) t>0
donde la temperatura inicial f : [0, L] R y los flujos hi , hd : [0, +) R son funciones conocidas.
Analogamente, si imponemos que las condiciones de frontera sean periodicas, el problema queda as:
2
ut = k uxx
x (0, L) t > 0
u(x, 0) = f (x) x (0, L)
u(0, t) = u(L, t) t>0
ux (0, t) = ux (L, t) t>0
Interpretacion del flujo en la ecuacion del calor. Para entender que es flujo de calor a traves
de la frontera, explicaremos una ley de conservacion referente a la evolucion de la temperatura en un
cuerpo 3D o una barra 1D sin focos ni sumideros de calor internos.
Empezamos por el caso 3D. Consideramos un cuerpo W R3 y sea S = W su frontera. Sea
u(x, y, z, t) una solucion del problema considerado en el ejemplo 1 e introducimos la funcion
Z
1
T (t) = u(x, y, z, t) dx dy dz
Vol(W ) W
que mide la temperatura promedio del cuerpo en el instante t. Su derivada es
k2
Z Z
1
T 0 (t) = ut (x, y, z, t) dx dy dz = u(x, y, z, t) dx dy dz
Vol(W ) W Vol(W ) W
k2 k2
I I
u
= (x, y, z, t) dS = h(x, y, z, t) dS.
Vol(W ) S n Vol(W ) S
Las propiedades que hemos usado son: derivada bajo el signo de la integral (primera igualdad), la
ecuacion del calor homogenea (segunda igualdad), el teorema de la divergencia de Gauss (tercera
igualdad, mirar el ejemplo 33 de los apuntes de Calculo Vectorial ) y las condiciones de frontera de
tipo Neumann (cuarta igualdad).
Por tanto, la tasa de variacion de la temperatura promedio es proporcional a la integral del flujo
de calor h sobre la superficie cerrada S. En realidad, la funcion h(x, y, z, t) determina exactamente a
que velocidad se escapa/entra el calor a traves de cada punto (x, y, z) S en cada instante t 0. La
temperatura promedio se mantiene constante cuando h 0 (es decir, H cuando el cuerpo esta termi-
camente aislado y ni entra ni sale calor por la frontera) o cuando S h dS 0 (es decir, cuando el
calor que entra por un lado se compensa exactamente con el calor H que sale por otro). Por contra, la
temperatura promedio aumenta/disminuye cuando la integral S h dS es positiva/negativa (es decir,
cuando hay una entrada/salida neta de calor).
A continuacion, estudiamos el caso 1D; o sea, la evolucion de la temperatura en una barra de
longitud L. Sea u(x, t) una solucion del problema considerado en la primera parte del ejemplo 2 e
introducimos la funcion
1 L
Z
T (t) = u(x, t) dx
L 0
que mide la temperatura promedio de la barra en el instante t. Su derivada es
1 L k2 L k2 k2
Z Z
x=L
T 0 (t) =
ut (x, t) dx = uxx (x, t) dx = ux (x, t) x=0 = hd (t) + hi (t) .
L 0 L 0 L L
Las propiedades que hemos usado son: derivada bajo el signo de la integral (primera igualdad), la
ecuacion del calor (segunda igualdad), el teorema fundamental del calculo (tercera igualdad) y las
condiciones de frontera (cuarta igualdad). Por tanto, la suma hd (t) + hi (t) nos dice cual es la tasa de
variacion de la temperatura promedio T (t). En otra palabras, las funciones hd (t) y hi (t) nos dicen a
que velocidad fluye el calor por el extremo derecho e izquierdo de la barra, respectivamente. Cuando
son positivas/negativas, tenemos una entrada/salida de calor por el extremo correspondiente.
Superposicion. Consideramos los dos PVIs de calor 1D en una barra de longitud L sin focos ni
sumideros de calor internos dados por
vt = k 2 vxx x (0, L) t > 0
wt = k 2 wxx x (0, L) t > 0
v(x, 0) = f (x) x (0, L) w(x, 0) = 0 x (0, L)
, .
v x (0, t) = 0 t > 0
w x (0, t) = hi (t) t>0
vx (L, t) = 0 t>0 wx (L, t) = hd (t) t>0
Ambos problemas tienen condiciones de frontera de tipo Neumann. La diferencia estriba en que el
primero tiene una unica condicion no homogena: la temperatura inicial, mientras que el segundo tiene
dos: las condiciones de frontera en los extremos de la barra.
Entonces, dadas dos soluciones cualesquiera v(x, t) y w(x, t) de estos problemas, su superposicion
(suma) u(x, t) = v(x, t) + w(x, t) es una solucion del PVI de calor 1D presentado en la primera parte
del ejemplo 2, que tiene tres condiciones no homogeneas.
En general, podemos trocear cualquier problema lineal en varios subproblemas de forma que cada
subproblema tenga pocas (quiza incluso solo una) ecuaciones/condiciones no homogeneas, siendo, por
tanto, mas simple que el problema original. En tal caso, si conseguimos resolver todos los subproblemas,
la superposicion (suma) de sus soluciones cumplira el problema original.
Homogeneizacion. Este truco es similar al anterior, pero en vez de trocear el problema original
en varios subproblemas simples, ahora queremos simplificarlo mediante un cambio de variables astuto.
Para fijar ideas, consideramos el PVI de calor 1D en una barra de longitud L = 1 sin focos ni
sumideros de calor internos con condiciones de frontera de tipo Dirichlet constantes
2
ut = k uxx 2 x (0, 1) t > 0
u(x, 0) = x x (0, 1)
.
u(0, t) = 1 t>0
u(1, t) = 2 t>0
La funcion v(x) = x + 1 cumple las condiciones de frontera: v(0) = 1 y v(1) = 2. Por tanto, si
realizamos el cambio de variables w(x, t) = u(x, t) v(x), el problema original se transforma en
wt = k 2 wxx x (0, 1) t > 0
w(x, 0) = x2 x 1 x (0, 1)
w(0, t) = 0
t>0
w(1, t) = 0 t>0
que es un problema bastante mas simple pues hemos homogeneizado las dos condiciones de frontera,
sin deshomogeneizar la EDP.
donde la fuerza externa F (x, t), el desplazamiento inicial f (x) y la velocidad inicial g(x) son funciones
conocidas. Este PVI tiene una unica solucion que viene dada por
Z x+ct Z t (Z x+c(ts) )
1 1 1
u(x, t) = f (x + ct) + f (x ct) + g(y) dy + F (y, s) dy ds.
2 2c xct 2c 0 xc(ts)
6 Depositado en http://www.ma1.upc.edu/ed/edp.pdf
La siguiente demostracion es opcional, no forma parte del temario de examen, pero deben entenderse
las consecuencias fsicas que se derivan de la formula.
Demostracion. Solo consideramos el caso F (x, t) 0; es decir, cuando no actua ninguna fuerza externa
sobre la cuerda. La idea principal consiste en realizar el cambio de variables
= x + ct, = x ct
para simplificar la EDP. Para eso debemos relacionar las derivadas parciales de la funcion transformada
v(, ) = u(x, t)
con las derivadas parciales de la funcion original u(x, t). Aplicamos repetidamente la regla de la cadena:
u v v
ux = = + = v + v
x x x
u v v
ut = = + = cv cv
t t t
ux v v v v v v
uxx = = + = + + + = v + 2v + v
x x x x x
ut v v v v v v
= c2 v 2v + v .
utt = =c c =c +c
t t t t t
v = 0
v (, ) = r() para alguna funcion r : R R arbitraria
v(, ) = p() + q() para algunas funciones p, q : R R arbitrarias
u(x, t) = p(x + ct) + q(x ct) para algunas funciones p, q : R R arbitrarias.
Es decir, la solucion general de la EDP de la cuerda vibrante infinita posee infinitas soluciones,
las cuales dependen de dos funciones arbitrarias, de la misma manera que la solucion general de una
EDO lineal de segundo orden dependa de dos constantes libres. Por tanto, para hallar la solucion del
PVI planteado, utilizaremos la misma estrategia seguida con las EDOs: determinar las dos funciones
libres imponiendo las dos condiciones iniciales. As pues, imponemos que
f (x) = u(x, 0) = p(x) + q(x), g(x) = ut (x, 0) = cp0 (x) cq 0 (x).
Derivando la primera ecuacion y multiplicando por c, se obtiene la relacion cp0 (x) + cq 0 (x) = cf 0 (x).
Combinando esta ultima relacion con la segunda ecuacion resulta que
1 0 1
p0 (x) = f (x) + g(x).
2 2c
Integrando esta ultima igualdad, se obtiene que
Z x Z x
1 1 1 1
p(x) = f (x) + g(y) dy + k, q(x) = f (x) p(x) = f (x) g(y) dy k,
2 2c 0 2 2c 0
con lo cual las funciones libres p(x) y q(x) quedan determinadas salvo una constante de integracion
comun k R. Finalmente,
Z x+ct
1 1
u(x, t) = p(x + ct) + q(x ct) = f (x + ct) + f (x ct) + g(y) dy,
2 2c xct
pues las dos constantes de integracion se cancelan entre si.
Depositado en http://www.ma1.upc.edu/ed/edp.pdf 7
Separacion de variables
El metodo de separacion de variables es un metodo para resolver problemas con simetras que poseen
una unica condicion (inicial o de frontera) no homogenea. No desarrollaremos una teora general, sino
que lo aplicaremos a tres o cuatro ejemplos concretos.
PVFs lineales de segundo orden. El susodicho metodo de separacion de variables requiere resolver
ciertos problemas de valores en la frontera (PVFs) asociados a EDOs lineales homogeneas a coeficientes
constantes de segundo orden mediante el metodo del polinomio caracterstico.
Definicion. P (m) = m2 + a1 m + a0 es el polinomio caracterstico de la EDO x00 + a1 x0 + a0 x = 0.
La importancia del polinomio caracterstico radica en que si x(t) = emt , entonces
x00 (t) + a1 x0 (t) + a0 x(t) = m2 emt + a1 memt + a0 emt = P (m)emt .
Usando esta relacion podemos expresar la solucion general de la EDO x00 + a1 x0 + a0 x = 0 en terminos
de las races de su polinomio caracterstico.
Lema. Sea xh (t) la solucion general de la EDO x00 + a1 x0 + a0 x = 0, donde a1 , a2 R. Entonces:
P (m) tiene dos races reales diferentes m1 , m2 R xh (t) = c1 em1 t + c2 em2 t ;
P (m) tiene una raz real doble m R xh (t) = em t (c1 + c2 t); y
P (m) tiene races complejas conjugadas m = i 6 R xh (t) = et c1 cos t + c2 sin t .
Demostracion. En la asignatura de Calculo 2 se explico que el conjunto de soluciones de una EDO
lineal homogenea de orden n es un subespacio vectorial de dimension n. Por tanto, si encontramos dos
soluciones linealmente independientes x1 (t) y x2 (t) de la EDO x00 + a1 x0 + a0 x = 0, podemos asegurar
que su solucion general es xh (t) = c1 x1 (t) + c2 x2 (t), con c1 , c2 R libres.
Si P (m) tiene dos races reales diferentes m1 y m2 , entonces x1 (t) = em1 t y x2 (t) = em2 t son
dos soluciones linealmente independientes.
Si P (m) tiene una raz real doble m , entonces x1 (t) = em t es una solucion y aplicando el
metodo de reduccion de orden (ver Calculo 2 ) obtenemos que x2 (t) = tem t es otra solucion.
Si P (m) tiene dos races complejas conjugadas m = i 6 R, entonces las funciones
x (t) := e( i)t = et et i = et cos t i sin t ,
8 Depositado en http://www.ma1.upc.edu/ed/edp.pdf
donde los instantes t1 6= t2 y los coeficientes ij y ai () son datos del problema. Los coeficientes a1 ()
y a0 () son constantes, pero dependen de un parametro . La funcion x(t) 0 siempre es una solucion
de estos PVFs. Es la llamada solucion trivial. Queremos saber para que valores del parametro R
existen soluciones no triviales.
Definicion. Estos valores son los valores propios (VAPS) y las soluciones no triviales son las funciones
propias (FUPS) del PVF. Un VAP es simple/doble cuando la dimension del subespacio vectorial
formado por sus FUPS es uno/dos.
Seguiremos los siguientes pasos para calcular los VAPs y sus FUPS:
1. Expresar la solucion general de la EDO homogenea en funcion del parametro R; es decir,
2. Imponer las condiciones de frontera, para as obtener un sistema lineal homogeneo de la forma
c1
A c = 0, c= .
c2
Proposicion. Listamos los VAPs y las FUPs de los siguientes PVFs clasicos.
1. PVF con condiciones de frontera tipo Dirichlet: x00 = x, x(0) = x(L) = 0.
VAPs: n = (n/L)2 , FUPs: xn (t) = sin(nt/L), para n 1.
2. PVF con condiciones de frontera tipo Neumann: x00 = x, x0 (0) = x0 (L) = 0.
VAPs: n = (n/L)2 , FUPs: xn (x) = cos(nt/L), para n 0.
3. PVF con condiciones de frontera mixtas: x00 = x, x(0) = x0 (L) = 0.
VAPs: n = (n + 1/2)2 2 /L2 , FUPs: xn (t) = sin (n + 1/2)t/L , para n 0.
4. PVF con condiciones de frontera periodicas: x00 = x, x(L) = x(L) y x0 (L) = x0 (L).
0 = 0 es un VAP simple de FUP c0 (t) 1.
n = (n/L)2 es VAP doble de FUPs cn (t) = cos(nt/L) y sn (t) = sin(nt/L), para n 1.
Depositado en http://www.ma1.upc.edu/ed/edp.pdf 9
Demostracion. Ninguno de estos cuatro PVFs tiene VAPs positivos. Por ejemplo, si x(t) es una FUP
de VAP , integrando por partes, usando que x00 = x y cualquiera de los tres primeros tipos de
condiciones de frontera anteriores (Dirichlet, Neumann o mixtas), vemos que
Z L Z L Z L RL 0
2 t=L 2 (x (t))2 dt
x0 (t) dt = x(t)x00 (t) dt x(t)x0 (t) t=0 = x(t) dt = R0L
0.
0 0 0 (x(t)) 2 dt
0
En el caso periodico se obtiene el mismo resultado, pero integrando en el intervalo [L, L].
Ademas, = 0 x0 (t) 0. Es decir, = 0 es un VAP si y solo si la funcion x(t) 1 cumple las
condiciones de frontera, lo cual solo sucede con las condiciones de tipo Neumann y las periodicas.
Finalmente, basta calcular los VAPs negativos: = 2 < 0 y sus FUPs. Recordamos que en tal
caso la solucion general de la EDO x00 = x es la combinacion lineal de funciones trigonometricas
xh (t) = c1 cos t + c2 sin t, c1 , c2 R.
1. Al imponer que xh (t) cumpla las condiciones de frontera tipo Dirichlet, obtenemos el sistema
1 0 c1
A c = 0, A = , c= .
cos L sin L c2
n 2
Por tanto, sin L = det[A ] = 0 = n = n 2
L , n Z = n = n = L , n 1.
En la ultima equivalencia hemos pasado de n Z a los enteros n 1. Las dos razones para
hacerlo son que n esta elevado al cuadrado (por tanto, los enteros negativos son superfluos) y
que la eleccion n = 0 no tiene sentido (pues estamos en el caso < 0).
Para calcular las FUPs de VAP = n = 2n , debemos resolver el sistema indeterminado
c1 = c1 cos 0 + c2 sin 0 = x(0) = 0 c1 = 0
= .
(1)n c1 = c1 cos n + c2 sin n = x(L) = 0 c2 R libre
Por tanto, xn (t) = sin(n t) = sin(nt/L) es una FUP de VAP n = (n/L)2 , para n 1.
2. Al imponer que xh (t) cumpla las condiciones de frontera tipo Newmann, obtenemos el sistema
0 1 c1
A c = 0, A = , c= .
sin L cos L c2
n 2
Por tanto, sin L = det[A ] = 0 = n = n 2
L , n Z = n = n = L , n 1.
Hemos pasado de n Z a los enteros n 1 por el mismo motivo que antes.
Para calcular las FUPs de VAP = n = 2n , debemos resolver el sistema indeterminado
= x0 (0) = 0
n c2 = c1 n sin 0 + c2 n cos 0 c1 R libre
= .
(1)n n c2 = c1 n sin n + c2 n cos n = x0 (L) = 0 c2 = 0
Por tanto, xn (t) = cos(n t) = cos(nt/L) es una FUP de VAP n = (n/L)2 , para n 0.
Hemos escrito n 0, en vez de n 1, para incluir que x0 (t) 1 es una FUP de VAP 0 = 0.
3. Al imponer que xh (t) cumpla las condiciones de frontera mixtas, obtenemos el sistema
1 0 c1
A c = 0, A = , c= .
sin L cos L c2
2 2
Y cos L = det[A ] = 0 = n = (n+1/2) L , n Z = n = 2n = (n+1/2)
L2
, n 0.
Hemos pasado de n Z a los enteros n 0 por lo de siempre.
Para calcular las FUPs de VAP = n = 2n , debemos resolver el sistema indeterminado
c1 = c1 cos 0 + c2 sin 0 = x(0) = 0 c1 = 0
= .
(1)n+1 n c1 = c1 n sin n L + c2 n cos n L = x0 (L) = 0 c2 R libre
Por tanto, xn (t) = sin(n t) = sin (n + 1/2)t/L es una FUP de VAP n = (n + 1/2)2 2 /L2 ,
para todo n 0.
10 Depositado en http://www.ma1.upc.edu/ed/edp.pdf
4. Al imponer que xh (t) cumpla las condiciones de frontera periodicas, obtenemos el sistema
0 2 sin L c1
A c = 0, A = , c= .
2 sin L 0 c2
n 2
Luego 4 sin2 L = det[A ] = 0 = n = n 2
L , n Z = n = n = L , n 1.
Hemos pasado una vez mas de n Z a los enteros n 1.
Para calcular las FUPs de VAP = n = 2n , debemos resolver el sistema indeterminado
0 = 2c2 sin(n L) = x(L) x(L) = 0 c1 R libre
= .
0 = 2n c1 sin(n L) = x0 (L) x0 (L) = 0 c2 R libre
Por tanto, cn (t) = cos(n t) = cos(nt/L) y sn (t) = sin(n t) = sin(nt/L) son dos FUPs
linealmente independientes del VAP doble n = (n/L)2 , para todo n 1.
Con esto hemos acabado la prueba de la proposicion.
Ejercicio. Resolver el PVF con condiciones mixtas x00 = x, x0 (0) = x(L) = 0.
Separacion de variables en la ecuacion de ondas 1D. Consideramos una ecuacion de ondas 1D
homogenea con condiciones de contorno de tipo Neumann homogeneas. Para simplificar supondremos
que la cuerda tiene longitud L = y que la soltamos, sin impulso, con un desplazamiento inicial
f (x) = 1 2 cos(3x). Tambien supondremos que no actua ninguna fuerza externa. Notamos por c la
velocidad a la que viajan las ondas por la cuerda. Las ecuaciones que modelan este problema son
utt = c2 uxx x (0, ) t R
u(x, 0) = 1 2 cos(3x) x (0, )
(1) ut (x, 0) = 0 x (0, ) .
u (0, t) = 0 t
x R
ux (, t) = 0 tR
La idea basica del metodo consiste en buscar soluciones en forma de variables separadas
u(x, t) = X(x)T (t)
de la parte homogenea del problema a resolver. En el caso anterior, todas las condiciones y ecuaciones
son homogeneas, salvo la referente al desplazamiento inicial, luego su parte homogenea es
utt = c2 uxx x (0, ) t R
ut (x, 0) = 0 x (0, )
(1)h .
ux (0, t) = 0 tR
ux (, t) = 0 tR
Ahora nos centramos en el problema asociado a la funcion T (t), pero teniendo en cuenta que
= n = n2 . En particular, la solucion general de la EDO T 00 + n2 c2 T = 0 es
Al imponer la condicion inicial T 0 (0) = 0, vemos que c2 = 0 y c1 R queda libre. Tras tomar c1 = 1,
que es la opcion mas simple, obtenemos la familia de funciones
As pues, hemos obtenido que todas las funciones de variables separadas de la familia
son soluciones del problema homogeneo (1)h . Estas funciones reciben el nombre de modos normales
y describen la forma en que vibra la cuerda. Concretamente, debido a la linealidad del problema
homogeneo (1)h , cualquier vibracion de la cuerda que estamos estudiando es una superposicion (suma)
de estos infinitos modos normales. En otra palabras, la solucion general del problema homogeneo (1)h
viene dada, al menos a nivel formal, por la serie
X X
u(x, t) = an un (x, t) = an cos(nx) cos(cnt),
n0 n0
se obteniene por inspeccion directa que a0 = 1, a3 = 2 y las demas amplitudes son nulas, luego
es una solucion del problema original. (En realidad es la unica, pero no lo probaremos.) As pues, en
este caso la vibracion de la cuerda es la superposicion de dos modos normales: el cero y el tres.
Observacion. La solucion anterior se puede reescribir en forma de dos ondas superpuestas viajando
en sentidos opuestos a velocidad c. Efectivamente, pues
u(x, t) = 1 2 cos(3x) cos(3ct) = 1 cos 3(x + ct) cos 3(x ct) = p(x + ct) + q(x ct),
con p(x) = 1/2 cos(3x) = q(x). (Hemos usado la relacion 2 cos a cos b = cos(a + b) + cos(a b).)
Ejercicio. Leer la entrada inglesa de Wikipedia sobre standing waves; es decir, sobre ondas estaciona-
rias. Ver alguno de los muchos videos que existen en Youtube sobre standing waves, en los cuales se
visualizan experimentalmente los primeros modos normales de vibracion de una cuerda.
Tambien se pueden ver algunos modos normales de vibracion de una membrana elastica rectangular
en un video de Youtube titulado Science fun. El experimento consiste en derramar sal encima de una
membrana negra que vibra por el sonido que emite un altavoz situado debajo para comprobar que los
modos normales cambian con la frecuencia del sonido.
Ejercicio. Escribir las dos EDOs que se obtienen al imponer que la funcion u(x, t) = X(x)T (t) cumpla
la EDP utt = kut + c2 uxx , escogiendo la opcion que proporciona una EDO lo mas simple posible
para la funcion X(x). Esta EDP recibe el nombre de ecuacion de la cuerda vibrante con friccion, pues
el termino kut proviene de una fuerza de friccion proporcional (y opuesta) a la velocidad.
12 Depositado en http://www.ma1.upc.edu/ed/edp.pdf
Desarrollos de Fourier. En el ultimo paso del ejemplo anterior, hemos conseguido determinar todos
los coeficientes libres por inspeccion directa. Cuando eso no sea posible, utilizaremos las siguientes
formulas, ya vistas en la asignatura Calculo 2, para calcular desarrollos de Fourier.
El desarrollo de Fourier completo de una funcion f : [L, L] R es
1 L
Z
a = f (x) cos(nx/L) dx,
n
a0 X
L L
f (x) + an cos(nx/L) + bn sin(nx/L), Z L
2 1
n1 bn =
f (x) sin(nx/L) dx.
L L
El desarrollo de Fourier en cosenos de una funcion f : [0, L] R es
2 L
Z
a0 X
f (x) + an cos(nx/L), an = f (x) cos(nx/L) dx.
2 L 0
n1
7. Hemos probado que cuando el tiempo tiende a infinito, la temperatura tiende al equilibrio
termico consistente en que la temperatura viene dada por la recta que une las temperaturas en
los extremos. Algo acorde con nuestra experiencia fsica, la cual nos ensena que el calor tiende
a distribuirse de la forma mas uniforme posible.
Ejercicio. Conectarse al enlace http://www-math.mit.edu/daimp/ y entender el applet de JAVA ti-
tulado Heat Equation que ejemplifica este fenomeno fsico.
Ejercicio. Probar que si substituimos las dos condiciones tipo Dirichlet constantes por dos condiciones
tipo Neumann homogeneas, entonces se cumple que
1 L
Z
lm u(x, t) = f (x) dx.
t+ L 0
La interpretacion fsica de este resultado es la siguiente. Las condiciones tipo Neumann homogeneas
equivalen a la existencia de un aislamiento termico perfecto en los extremos que impide que el calor
escape o entre, luego tan solo puede redistribuirse internamente. Por tanto, la temperatura tiende a
un valor constante y este valor debe coincidir con el promedio de la temperatura inicial.
Separacion de variables en la ecuacion de Poisson 2D en dominios rectangulares.
Objetivo. En este ultimo ejemplo del metodo de separacion de variables, vamos a resolver una ecuacion
de Poisson 2D en un dominio rectangular con condiciones de contorno de tipo Dirichlet. Dos de las
cuatro condiciones de contorno son no homogeneas. Antes de separar variables, homogeneizaremos
tanto la ecuacion de Poisson (es decir, la transformaremos en una ecuacion de Laplace) como una
condicion de contorno, mediante un cambio de variables.
2. Imponer que la funcion w(x, y) = X(x)Y (y) cumpla la parte homogenea del problema (4).
Escribir el PVF asociado a la funcion Y (y) y el problema asociado a la funcion X(x).
3. Resolver el PVF asociado a la funcion Y (y).
4. Teniendo en cuenta los VAPs del PVF anterior, resolver el problema asociado a X(x).
5. Calcular la solucion general de la parte homogenea del problema (4).
6. Resolver el problema original.
Al imponer que la funcion w(x, y) = u(x, y) v(x, y) cumpla las condiciones de contorno
correspondientes a los lados inferior, superior e izquierdo queda
0 = w(0, y) = u(0, y) v(0, y) = 0 v(0, y) = v(0, y) = 0
0 = w(x, 0) = u(x, 0) v(x, 0) = 0 v(x, 0) = v(x, 0) = 0
0 = w(x, 2) = u(x, 2) v(x, 2) = 2x2 v(x, 2) = v(x, 2) = 2x2
Necesitamos una funcion v(x, y) que cumpla estas cuatro condiciones. Para simplificar los calcu-
los, buscamos esta funcion en forma de variables separadas: v(x, y) = X(x)Y (y). Entonces, las
cuatro condiciones anteriores equivalen a
X 00 (x)Y (y) + X(x)Y 00 (y) = 2y, X(0) = 0, Y (0) = 0, X(x)Y (2) = 2x2 .
Una posible solucion es tomar X(x) = x2 y Y (y) = y. Es decir, v(x, y) = x2 y, luego
g(y) = w(, y) = u(, y) v(, y) = 1 2 y.
2. Al imponer que la funcion w(x, y) = X(x)Y (y) cumpla:
La ecuacion de Laplace wxx + wyy = 0 se obtiene que X 00 (x)Y (y) + X(x)Y 00 (y) = 0, luego
X 00 (x) Y 00 (y)
= = R.
X(x) Y (y)
La condicion de contorno w(0, y) = 0, se obtiene que X(0) = 0.
La condicion de contorno w(x, 0) = 0, se obtiene que Y (0) = 0.
La condicion de contorno w(x, 2) = 0, se obtiene que Y (2) = 0.
Por tanto, obtenemos dos problemas separados:
00 00
X (x) + X(x) = 0 Y (y) = Y (y)
(a) (b) .
X(0) = 0 Y (0) = 0 = Y (2)
El problema (b) es un PVF con condiciones de tipo Dirichlet asociado a la funcion Y (y).
3. Ya vimos que los VAPs y las FUPs del PVF (b) son:
VAPs: = n = n2 /4
n 1.
FUPs: Y (y) = Yn (y) = sin(ny/2)
4. La EDO X 00 (x) + n X(x) = 0 es lineal, homogenea y a coeficientes
constantes. Su polinomio
caracterstico es P (m) = m2 + n y sus races son m1,2 = n = n/2. Por tanto, la
solucion general de esta ecuacion es
X(x) = c1 enx/2 + c2 enx/2 , c1 , c2 R.
Al imponer la condicion adicional 0 = X(0) = c1 + c2 obtenemos que c2 = c1 , luego
X(x) = c1 (enx/2 enx/2 ), c1 R.
Tomando c1 = 1/2, obtenemos la familia de funciones
enx/2 enx/2
Xn (x) = = sinh(nx/2), n 1.
2
5. As pues, los modos normales (las FUPs) de la parte homogenea del problema (4) son
wn (x, y) = Xn (x)Yn (x) = sinh(nx/2) sin(ny/2), n 1.
En particular, resulta que, por linealidad, todas las series de la forma
X X
w(x, y) = n wn (x, y) = n sinh(nx/2) sin(ny/2)
n1 n1
6. En el paso anterior los coeficientes n haban quedado libres, pero ahora los determinamos
para as obtener la solucion final del problema (4) imponiendo la unica condicion no
homogenea del problema; a saber, la condicion de contorno en el lado derecho del rectangulo:
X X X
g(y) = w(, y) = n wn (, y) = n sinh(n/2) sin(ny/2) = bn sin(ny/2)
n1 n1 n1
donde hemos notado bn = n sinh(n/2). En la seccion sobre desarrollos de Fourier vimos que
1 2 (1)n 2 1 (1)n
Z
bn = (1 2 y) sin(ny/2) dy = 4 2 + , n1
0 n n
son los coeficientes de Fourier del desarrollo en senos de la funcion g(y) = 1 2 y en el intervalo
[0, 2]. Y deshaciendo el cambio de variables w(x, y) = u(x, y) v(x, y), la solucion final es
X bn
u(x, y) = v(x, y) + w(x, y) = x2 y + sinh(nx/2) sin(ny/2).
sinh(n/2)
n1