Notas de EDP

J. Duoandikoetxea - L. Escauriaza

29 de noviembre de 2013
Capı́tulo 1

Introducción

1.1. ¿Qué son las EDP’s?

Una EDP es una ecuación con una incógnita u que es función de x =
(x1 , . . . , xn ), a la que se pide que satisfaga
F (x, u, Du, D2 u, . . . , Dm u) = 0.
Aquı́ Dj u representa todas las derivadas de u de orden j y
n
∂ |α| u X
Dα u = α1 αn , |α| = αi , si α = (α1 , . . . , αn ) ∈ Nn .
∂x1 . . . ∂xn
i=1

Terminologı́a
Orden de la ecuación. Es el mayor orden de derivación de u que aparece
en la ecuación.
Ecuación lineal. Aquella en la que los términos en los que aparecen u o
sus derivadas son lineales.
Todas las ecuaciones que estudiamos con detalle en el curso son lineales
y de orden menor que tres.
Ecuación cuasilineal. La misma definición que la anterior pero ahora
referida sólo a los términos que contienen a las derivadas de orden igual al
orden de la ecuación.
Los siguientes conceptos y teorema serán muy últiles en este curso:
∂f ∂f
∇f = ( ∂x 1
, . . . , ∂x 1
).
La divergencia de un campo vectorial, F = (F1 , . . . , Fn ), y el Lapla-
ciano de una función f definidos en un abierto Ω ⊂ Rn
∂F1 ∂Fn
∇·F = + ··· + .
∂x1 ∂x1
4f = ∇2 f = ∇ · (∇f ).

1
Teorema 1 (Teorema de la divergencia). Ω es una abierto acotado de Rn
con frontera regular, F : Ω −→ Rn es una campo vectorial en Ω, F ∈ C 1 (Ω),
ν es el vector normal exterior a ∂Ω. Entonces,
Z Z
∇ · F (x) dx = F · ν dσ.
Ω ∂Ω

1.2. Ejemplos de ecuaciones

Los tres primeros ejemplos son las ecuaciones que se estudian en el curso.

1. Ecuación de ondas. Se busca u(x, t) con x ∈ Rn , tal que

utt − c2 4u = F (x, t).
Aquı́ c > 0 es una constante, F una función dada y 4u = ux1 x1 +
· · · + uxn xn , es el laplaciano de u en las coordenadas espaciales (t es el
tiempo). Si n = 1, se reduce a, utt − c2 uxx = F .
2. Ecuación del calor. Se busca u(x, t) con x ∈ Rn tal que
ut − 4u = F (x, t).
Si n = 1, se reduce a ut − c2 uxx = F .
3. Ecuación del potencial. Se busca u(x) con x ∈ Rn tal que
4u = F (x).
Si n = 2, se reduce a uxx + uyy = F .
Una ecuación que no estudiaremos es la siguiente.
4. Ecuación de Schrödinger. Se busca u(x, t), con x ∈ Rn y tal que
iut − 4u = F (x, t).

En lugar de ecuaciones se pueden considerar sistemas.

5. Ecuaciones de Cauchy-Riemann. Se buscan u(x, y), v(x, y) definidas
en R2 tales que (
ux = vy ,
uy = −vx .

6. Ecuaciones de Navier-Stokes. Encontrar, u : R3 × [0, +∞) −→ R3 y

p : R3 × [0, +∞) −→ R, tal que se verifiquen las ecuaciones

∂t u − 4u + u · ∇u = −∇p,

∇ · u = 0,

u(0) = a,


2
 en R3 × (0, +∞), donde a : R3 −→ R3 es una campo vectorial con
divergencia nula, ∇ · a = 0.

7. Sistema de elasticidad. Encontrar u : R3 −→ R3 , tal que se verifiquen

las ecuaciones
µ4u + (µ + λ)∇ (∇ · u) = 0
en Ω ⊂ R3 , donde µ, λ > 0 son las constantes de elasticidad o de
Lamé del sólido Ω.

1.3. Cambio de variables

Repaso de la regla de la cadena para funciones de varias variables. Sea
v(y) una función de clase C m definida en D ⊂ Rn y ϕ : D0 ⊂ Rn → D,
biyectiva y de clase C m ; se define u(x) = v(ϕ(x)). Entonces
n
∂u X ∂v ∂ϕi
(x) = (ϕ(x)) · (x).
∂xj ∂yi ∂xj
i=1

Con la función inversa de ϕ se puede escribir la relación en la forma v(y) =

u(ϕ−1 (y)) y escribir todas las derivadas de v en términos de las de u.

Ejemplos
1. La ecuación de primer orden, ux + uy = 0, se transforma en vξ = 0, por
el cambio de variables independientes
(
ξ = x + y,
η = x − y.

2. Se define en el plano el cambio de variable

( (
ξ = x + ct , x = (ξ + η)/2 ,
cuya inversa es
η = x − ct , t = (ξ − η)/2c .

u(x, t) y v(ξ, η) se relacionan por u(x, t) = v(x + ct, x − ct). Escribir la

expresión utt − c2 uxx en términos de las derivadas de v. Se debe llegar a
−4c2 vξη .
3. El laplaciano en coordenadas polares. Se define en el plano el cambio de
variable a coordenadas polares habitual x = r cos θ, y = r sin θ, donde r > 0
y −π < θ < π. Si v(r, θ) = u(r cos θ, r sin θ), se verifica que
1 1
4u = uxx + uyy = vrr + vr + 2 vθθ .
r r

3
4. El laplaciano en coordenadas cilı́ndricas y esféricas. Las coordenadas
cilı́ndricas son x = r cos θ, y = r sin θ y z = z, donde r > 0, 0 < θ < 2π y
−∞ < z < +∞. Si v(z, θ, z) = u(r cos θ, sin θ, z),
1 1
4u = uxx + uyy + uzz = vrr + vr + 2 vθθ + vzz .
r r
Las coordenadas esféricas son (ρ, θ, ϕ) y su relación con las cartesianas
es x = ρ cos θ sin ϕ, y = ρ sin θ sin ϕ y z = ρ cos ϕ, ρ > 0, 0 < θ < 2π y
0 < ϕ < π. En este caso, si v(ρ, θ, ϕ) = u(ρ cos θ sin ϕ, ρ sin θ sin ϕ, ρ cos ϕ),
se verifica que
2 1 1 cot ϕ
4u = vρρ + vρ + 2
vθθ + 2 vϕϕ + 2 vϕ .
ρ (ρ sin ϕ) ρ ρ

Cambio de función incógnita

Si en la ecuación ut − uxx − αu = 0 hacemos el cambio v(x, t) =
e−αt u(x, t), queda la ecuación del calor para v.
Si en la ecuación, ut = u2x + uxx , hacemos el cambio v = eu , desaparece
el término no lineal y v es solución de vt − vxx = 0.
La siguiente fórmula es muy útil:

4(f g) = f 4g + g4f + 2∇f · ∇g.

Si v = e−ϕ u, calcular e−ϕ 4u en términos de v.

1.4. Problema de Cauchy. Condiciones iniciales y

de contorno. Problema bien planteado
Una ecuación puede tener en principio muchas soluciones, ¿qué datos se
necesitan para determinar una de ellas? En EDOs, podı́a ser el valor de la
función y algunas derivadas en un punto; ahora eso no es suficiente.
Por ejemplo, sea ux = 0 en el plano. La solución es u(x, y) = ϕ(y), para
alguna función ϕ; por tanto:
- si nos dan u(a, b), tenemos u(x, b) para todo x;
- si nos dan u(a, b) y u(a0 , b) y no son iguales, no hay solución y si son
iguales, uno de los datos es innecesario.
En general, si a(x, y), b(x, y) y c(x, y) son funciones regulares, el proble-
ma de Cauchy asociado a una ecuación de primer orden consiste en encontrar
una solución local u de la ecuación

aux + buy = c (1.1)

que verifique la condición de Cauchy

u(f (s), g(s)) = h(s) (1.2)

4
a lo largo de una curva S = {(f (s), g(s))} y cerca de un punto P0 =
(f (s0 ), g(s0 )) de S. Veremos que esto es posible si sabemos a priori que
la curva S es no “caraterı́stica”en P .
Una curva S es caraterı́stica en uno de sus puntos P si su vector tan-
gente en P es paralelo al vector (a(P ), b(P )). En caso contrario, se dice no
caraterı́stica en P . Es decir, se ha de verificar que
 
a(f (s), g(s)) b(f (s), g(s))
 6= 0. (1.3)
f 0 (s) g 0 (s)

 

Las curvas caracterı́sticas asociadas a la ecuación (1.1) que pasan por

(f (s), g(s)) se obtienen resolviendo el sistema de EDO’s con condiciones
iniciales 
∂x
 ∂t = a(x, y),

∂y
∂t = b(x, y),
(1.4)

x(0, s) = f (s), y(0, s) = g(s).


y si u es solucion de clase C 1 de la ecuación lineal en un entorno de P y

z(t, s) = u(x(t, s), y(t, s)), se verifica que
∂z
= c(x(t, s), y(t, s))
∂t
y Z t
z(t, s) = h(s) + c(x(τ, s), y(τ, s))dτ. (1.5)
0
En particular,
Z t
u(x(t, s), y(t, s)) = h(s) + c(x(τ, s), y(τ, s))dτ,
0

y los valores de u están determinados a lo largo de una curva caracterı́stica

que pasa por S por el valor de u en la intersección de la caracterı́stica con
S y por los valores de c a lo largo de la caracterı́stica.
Recı́procamente, si a, b, c y (f (s), g(s), h(s)) son de clase C 1 en entornos
de P y s0 respectivamente, la teorı́a de las EDO’s asegura que el sistema
(1.4) admite una solución de clase C 1 para (t, s) cerca de (0, s0 ). Es decir,
existe un δ > 0 tal que la transformación, x = x(t, s), y = y(t, s), está bien
definida y es de clase C 1 en (−δ, δ) × (s0 − δ, s0 + δ). Además, si S es no
caraterı́stica en P , el Jacobiano de la transformación
(
x = x(t, s),
(1.6)
y = y(t, s),

en (0, s0 ) es  
a(f (s0 ), g(s0 )) b(f (s0 ), g(s0 ))
 6= 0
f 0 (s0 ) g 0 (s0 )

 

5
y el teorema de la función inversa garantiza que el sistema (1.6) se puede
invertir para (x, y) cerca de P y (t, s) cerca de (0, s0 ). Si definimos entonces

u(x, y) = z(t(x, y), s(x, y)),

donde z se define mediante (1.5), la regla de la cadena muestra que u es una

función de clase C 1 en un entorno de P . Como z(t, s) = u(x(t, s), y(t, s)),
haciendo t = 0, vemos que u verifica la condición de Cauchy (1.2) en S cerca
P . Derivando la misma identidad respecto a t, obtenemos que
∂z ∂x ∂y
= ux (x(t, s), y(t, s)) + uy (x(t, s), y(t, s))
∂t ∂t ∂t
y como, ∂z∂t = c(x(t, s), y(t, s)), se sigue que u verifica la ecuación (1.1) en
un entorno de P .
Por lo tanto, el problema de Cauchy asociado a la ecuación lineal de orden
uno tiene una única solución local de clase C 1 en un entorno de P si S no
es una curva caracterı́stica en P , a(x, y), b(x, y), c(x, y) y (f (s), g(s), h(s))
son funciones de clase C 1 cerca de P y s0 respectivamente.
Otro ejemplo es el asociado a la ecuación

ux − uy = 0.

La solución general en este caso es de la forma

u(x, y) = ϕ(x + y),

donde ϕ es cualquier función de clase C 1 . Las soluciones son constantes en

las rectas x + y = k y las curvas S en la que se dan los datos de Cauchy
deben contener un único punto de cada una de estas rectas si queremos que
el problema de Cauchy esté bien planteado para todos los datos.
Las rectas y = k, en el primer ejemplo y las rectas x + y = k, en el
segundo, son caracterı́sticas y ponen limitaciones al lugar en que se deben
dar los datos para que el problema de Cauchy asociado tenga siempre sentido
(poder dar datos arbitrarios sobre S). Las soluciones de ambas ecuaciones
son constantes a lo largo de las caracterı́sticas y los datos h no pueden ser
arbitrarios a lo largo de las mismas.
Se puede comprobar, que la condición de que S sea no caracterı́stica en
P para la ecuación lineal, equivale a pedir que las ecuaciones
(
aux + buy = c
u(f (s), g(s)) = h(s)

determinen unı́vocamente el valor de todas las derivadas en S y cerca de P

de una posible solucion u de la ecuación en un entorno de P . Es decir, que
 
a(f (s), g(s)) b(f (s), g(s))
 6= 0. (1.7)
f 0 (s) g 0 (s)

 

6
 En general, el problema de Cauchy asociado a una ecuación de primer
orden consiste en encontrar en un entorno de P una solución u de clase C 1
del problema (
F (x, y, u, ux , uy ) = 0
.
u(f (s), g(s)) = h(s)
Se dice que la curva en el espacio, (f (s), g(s), h(s)), por la que queremos
pasar una superficie, z = u(x, y), que verifique la ecuación es no caraterı́stica,
si todas las derivadas parciales de u están univocamente determinadas en P
por la ecuación y la condición de Cauchy en S.
En las ecuaciones no lineales el ser no caraterı́stica depende en general
del dato de Cauchy h que damos sobre S y no sólo de S. Por ejemplo, para
la ecuación cuasilineal

a(x, y, u)ux + b(x, y, u)uy = c(x, y, u), (1.8)

la curva, (f (s), g(s), h(s)), es no caracterı́stica si

 
a(f (s), g(s), h(s)) b(f (s), g(s), h(s))
 6= 0, para s cerca de s0 ,
f 0 (s) g 0 (s)

 

condición que depende de los valores de h en S. En este caso, si a, b, c y

(f (s), g(s), h(s)) son respectivamente de clase C 1 en entornos de P y s0 , la
teorı́a de las EDO’s asegura que el sistema

∂x
∂t = a(x, y, z),



 ∂y = b(x, y, z),

∂t (1.9)
∂z
∂t = c(x, y, z),




x(0, s) = f (s), y(0, s) = g(s), z(0, s) = h(s),

tiene una solución de clase C 1 para (t, s) cerca de (0, s0 ). Además, si S es

no caraterı́stica en P , el Jacobiano de la transformación
(
x = x(t, s),
(1.10)
y = y(t, s),

en (0, s0 ) es
 
a(f (s0 ), g(s0 ), h(s0 )) b(f (s0 ), g(s0 ), h(s0 ))
 6= 0
f 0 (s0 ) g 0 (s0 )

 

y el teorema de la función inversa garantiza otra vez que el sistema (1.10)

se puede invertir para (x, y) cerca de P y (t, s) cerca de (0, s0 ). Si definimos

u(x, y) = z(t(x, y), s(x, y)),

7
la regla de la cadena muestra que u es de clase C 1 en un entorno de P y
haciendo t = 0 en la relación, z(t, s) = u(x(t, s), y(t, s)), vemos que u verifica
la condición de Cauchy (1.2) en S y cerca de P . Derivando otra vez la misma
identidad respecto a t, concluimos que
∂z ∂x ∂y
= ux (x(t, s), y(t, s)) + uy (x(t, s), y(t, s)) ,
∂t ∂t ∂t
que combinado con las identidades que verifican x(t, s), y(t, s) y z(t, s) en
(1.9), implica otra vez que u es solución de (1.8) en un entorno de P .
En las ecuaciones de segundo orden puede haber caracterı́sticas pero el
concepto es algo más complicado. En este caso, el problema de Cauchy con-
siste en encontrar una función u(x, y) de clase C 2 que verifique la ecuación
de segundo orden en un entorno de P y cuyas derivadas de orden menor
que dos sobre la curva S son conocidas . Por ejemplo, si S = {y = 0}, nos
podemos plantear encontrar u(x, y) de clase C 2 , definida cerca de (0, 0) y
tal que 
uyy = F (x, y, u, ux , uy , uxx , uxy ),

u(x, 0) = ϕ0 (x),

uy (x, 0) = ϕ1 (x).


Se comprueba fácilmente que y = 0 no es una curva caracterı́stica para

este último problema de Cauchy:
Todas las derivadas de una posible solución están unı́vocamente deter-
minadas en la curva, y = 0, por las derivadas del dato de Cauchy, ϕ0 , ϕ1 y
las derivadas de F .

1.5. Clasificación de las ecuaciones lineales de se-

gundo orden.
La ecuación cuasilineal de segundo orden en dos variables es de la forma
auxx + 2buxy + cuyy = d, (1.11)
donde a, b, c y d dependen de x, y, u, ux y uy . El problema de Cauchy
consiste en encontrar una solución de (1.11) con valores compatibles de u,
ux y uy sobre una curva S. Si S = {(f (s), g(s))} el dato de Cauchy es
u(f (s), g(s)) = h(s), ux (f (s), g(s)) = φ(s), uy (f (s), g(s)) = ψ(s) (1.12)
y debe satisfacer la condicion de compatibilidad
h0 (s) = φ(s)f 0 (s) + ψ(s)g 0 (s). (1.13)
Esto equivale a dar sobre S los valores de u y los de su derivada normal
−ux g 0 + uy f 0
u(f (s), g(s)) = h(s), p 0 = ℵ(s). (1.14)
f 2 + g02

8
 Si u(x, y) es solución de (1.11), la ecuación (1.11) y las condiciones de
compatibilidad sobre S para los valores de ux y uy implican que

auxx + 2buxy + cuyy = d, f 0 uxx + g 0 uxy = φ0 , f 0 uxy + g 0 uyy = ψ 0 ,

que es un sistema con coeficientes que dependen de s y del dato de Cauchy

y que tiene por incógnitas las derivadas segundas de u en S. Éste sistema
tendrá solución única si
 
 a 2b c 
0 0
4 = f g 0 0  = ag 2 − 2bf 0 g 0 + cf 2 6= 0.
 0 
 0 f 0 g0

Si 4 = 6 0 en (f (s0 ), g(s0 ), h(s0 ), φ(s0 ), ψ(s0 )) y derivando la ecuación

(1.11), se concluye que todas las derivadas de una solución de (1.11), (1.12) (o
de (1.11) y (1.14)), quedan unı́vocamente determinadas a lo largo de la curva
(f (s), g(s), h(s), φ(s), ψ(s)) y cerca de (f (s0 ), g(s0 ), h(s0 ), φ(s0 ), ψ(s0 )) por
los datos de Cauchy y los valores de los coeficientes de la ecuación (1.11) a
lo largo de la curva.
Las curvas caracterı́sticas para (1.11), (1.12) son las curvas en R5

(f (s), g(s), h(s), φ(s), ψ(s))

tales que 4 ≡ 0 a lo largo de su recorrido. Si la ecuación es lineal, la condición

de ser caracterı́stica es independiente del dato de Cauchy (h(s), φ(s), ψ(s))
y sólo depende de la proyección de la curva en el plano xy, a(x, y), b(x, y) y
c(x, y), como ocurre con las ecuaciones de primer orden; es decir, de la traza
de S = {(f (s), g(s))}.
En particular, en el caso lineal y si S está dada implicitamente por la
ecuación, φ(x, y) = cte. ((f 0 , g 0 ) es paralelo al vector (φy , −φx )), la condición
de ser curva caraterı́stica es

∇φ 6= 0 y a(x, y)φ2x + 2b(x, y)φx φy + c(x, y)φ2y = 0. (1.15)

Recı́procamente, si φ verifica (1.15), las curvas φ(x, y) = cte. son curvas

caracterı́sticas de la ecuación.
La ecuación lineal de segundo orden es

a(x, y)uxx + 2b(x, y)uxy + c(x, y)uyy = d(x, y, u, ux , uy )

y para encontrar una forma equivalente y más sencilla de la ecuación alre-

dedor de un punto P = (x0 , y0 ), se puede consider un cambio de variables
abstracto (
ξ = α(x, y),
η = β(x, y),

9
con funciones regulares α, β con Jacobiano no nulo. En tal caso, se puede
mostrar que u, como función de las variables (ξ, η), es solución de la ecuación
˜ η, u, uξ , uη ),
ãuξξ + 2b̃uξη + c̃uηη = d(ξ,

donde

ã = aαx2 + 2bαx αy + cαy2 = A∇α · ∇α,

b̃ = aαx βx + b (αx βy + αy βx ) + cαy βy = A∇α · ∇β,
c̃ = aβx2 + 2bβx βy + cβy2 = A∇β · ∇β,

y A es la matriz simétrica  
a b
.
b c
Ejemplo 1. Las caracterı́sticas de las ecuaciones, uyy − c2 uxx = 0 y de
uxy = 0, son respectivamente las curvas, x ± cy = const. y x = const.,
y = const.
Si a y c son nulas en un entorno de P y b(P ) 6= 0, la ecuación (1.11)
se reduce dividiendo por b y cerca de P a una ecuación del tipo (1.16).
Supuesto que c > 0, la ecuación caracterı́stica (1.15) tiene dos soluciones
independientes definidas en un entorno de P , α(x, y) y β(x, y), si E = ac −
b2 < 0 cerca de P ; es decir, por cada punto en el entorno pasan dos curvas
caracterı́sticas no paralelas. En este caso el operador se dice hiperbólico y
con el cambio de variable de variables, ξ = α(x, y), η = β(x, y), la ecuación
original se transforma en

uξη = F (ξ, η, u, uξ , uη ). (1.16)

Si además hacemos el cambio, X = 21 (ξ + η), T = 1

2 (ξ − η), la última
ecuación se transforma en una ecuación de tipo ondas

uT T − uXX = H(X, T, u, uX , uT )

y la ecuación de ondas se considera la forma canónica de las ecuaciones de

tipo hiperbólico.
Ejemplo 2. La ecuación uxx + uyy = 0 no tiene curvas caracterı́sticas.
Si E > 0 en un entorno de P , la ecuación se dice de tipo elı́ptico. En
este caso, no hay curvas caraterı́siticas y (1.15) tiene una solución a valores
complejos φ(x, y). Si
1 1
α = (φ + φ̄) , β(x, y) = (φ − φ̄),
2 2i
resulta que la ecuación original se transforma en

uξξ + uηη = G(ξ, η, u, uξ , uη ),

10
y la ecuación de Laplace se considera como la forma canónica de las ecua-
ciones de tipo elı́ptico.
Ejemplo 3. Las caracterı́sticas de la ecuación uxx −uy = 0 son las curvas
y = const.
Si E ≡ 0 en un entorno de P , la ecuación se dice de tipo parabólico. Por
cada punto pasa una única caracterı́stica y la ecuación caracterı́stica (1.15)
tiene una solución α(x, y) con gradiente no nulo en P . Además, E ≡ 0
implica que A∇α ≡ 0 en un entorno de P . En este caso, tomando una
función arbitraria β(x, y) (si αx 6= 0, considerar β(x, y) = y y si αy 6= 0,
hacer β(x, y) = x) que complete el cambio de variables, la nueva ecuación
es
uηη = M (ξ, η, u, uξ , uη ).
La ecuación del calor se considera como la forma canónica de las ecua-
ciones de tipo parabólico.

1.6. Teorema de Cauchy-Kowalevsky

Consideramos la versión para n = 2 en ecuaciones de primer orden con
datos de Cauchy en y = 0 y con uy despejada.

Teorema 2. Se considera el problema de Cauchy

(
uy = F (x, y, u, ux ) ,
(P1 )
u(x, 0) = h(x)

y supongamos que h es analı́tica en un entorno de 0 y que F lo es en un

entorno de (0, 0, h(0), h0 (0)). Entonces, existe una única solución analı́tica
de (P1 ) en un entorno de (0, 0).

Consideramos la versión para n = 2 en ecuaciones de segundo orden con

datos de Cauchy en y = 0 y con uyy despejada.

Teorema 3. Se considera el problema de Cauchy


uyy = F (x, y, u, ux , uy , uxx , uxy ) ,

(P2 ) u(x, 0) = ϕ0 (x) ,

uy (x, 0) = ϕ1 (x) .


y supongamos que ϕ0 y ϕ1 son funciones analı́ticas en un entorno de 0 y

que F lo es en un entorno de

(0, 0, ϕ0 (0), ϕ00 (0), ϕ1 (0), ϕ000 (0), ϕ01 (0)).

Entonces, existe una única solución analı́tica de (P2 ) en un entorno de (0, 0).

11
 De la ecuación y de las condiciones de Cauchy, todos los valores de

∂ i+j u
cij = (0, 0),
∂xi ∂y j
quedan determinados por los datos y la ecuación. Con estos, se puede escribir
formalmente la serie de Taylor de una posible solución u en (0, 0), como
+∞
X cij i j
u(x, y) = xy .
i!j!
i,j=0

Lo interesante, que no haremos en clase, es comprobar que la serie con-

verge en un entorno de (0, 0) y define efectivamente una solución de (P1 ) o
(P2 ) en un entorno de (0, 0).
El problema de Cauchy
(
F (x, y, u, ux , uy ) = 0,
(1.17)
u(f (s), g(s)) = h(s),

es equivalente por medio del cambio de variables

(
x = f (s) + tg 0 (s),
y = g(s) − tf 0 (s),

que transforma localmente la curva S = {(f (s), g(s))} en el plano (x, y) en

la curva t = 0 en el plano (t, s), a un problema del tipo
(
G(s, t, u, us , ut ) = 0,
(1.18)
u(s, 0) = h(s)

y que S sea no caracterı́stica en P = (f (s0 ), g(s0 )) para el problema de C

auchy (1.17) es equivalente a que t = 0 no lo sea para (1.18) en (s0 , 0); es
decir, que (1.18) se pueda escribir como
(
ut = H(s, t, u, us ),
u(s, 0) = h(s).

Además, como la inversa o la composición de funciones analı́ticas de va-

riable real son también analı́ticas de variable real, el teorema de Teorema
de Cauchy-Kowalevsky garantiza que si F (x, y, z, p, q) en (1.17) es analı́tica
en (x, y, z, p, q), f , g, h son analı́ticas en s0 y S es no caracterı́stica en P ,
entonces (1.17) tiene una única solución analı́tica
+∞
X cij
u(x, y) = (x − f (s0 ))i (y − g(s0 ))j
i!j!
i,j=0

12
en algún entorno de P .
Un enunciado análogo se puede dar para el problema de Cauchy asociado
a la ecuación de segundo orden

F (x, y, u, ux , uy , uxx , uxy , uyy ) = 0,

u(f (s), g(s)) = h(s),

 ∂u
∂ν (s, 0) = χ(s),

si F , f , g, h, χ son analı́ticas de variable real y S no es caracterı́stica para

este problema en P .
La unicidad en el teorema de Cauchy-Kowalevsky se refiere sólo a solu-
ciones analı́ticas y no se excluye la existencia de posibles soluciones que no
sean analı́ticas. Las demostraciones de estos resultados las encontraréis en
el libro de F. John, Partial Differential Equations.

Condiciones iniciales o condiciones de contorno

Estos conceptos los entenderemos con algunos ejemplos asociado a la
ecuación de Laplace y a la Ecuación del calor:

El problema de contorno o problema de Dicihlet asociado a la ecuación

de Laplace en un dominio Ω ⊂ Rn y con dato de contorno o de Dirichlet
ϕ en C(∂Ω), consiste en encontrar u en C 2 (Ω) ∩ C(Ω), que verifique
(
4u = 0 en Ω,
u=ϕ en ∂Ω.

El problema de condiciones iniciales asociado a la ecuación del calor en

Rn y con dato inicial ϕ en C(Rn ), consiste en encontrar u en C 2,1 (Rn ×
(0, +∞)) ∩ C(Rn × [0, +∞)) tal que
(
4u − ∂t u = 0 en Rn × (0, +∞),
u(x, 0) = ϕ(x) en Rn .

El problema de condiciones iniciales con condición de contorno nula de

tipo Dirichlet para la ecuación del calor en Ω ⊂ Rn y con dato inicial f
en C(Ω), consiste en encontrar u en C 2,1 (Ω×(0, +∞))∩C(Ω×[0, +∞))
que verifique

4u − ∂t u = 0
 en Ω × (0, +∞),
u(x, 0) = f (x) en Ω,

u=0 en ∂Ω × [0, +∞).


13
 La condición de contorno se dice de tipo Neumann nula, cuando la
condición en la frontera lateral es:
∂u
= ∇u · ν = 0, en ∂Ω × (0, +∞),
∂ν
donde ν es el vector normal unitario exterior a ∂Ω.

Problema bien planteado

J. Hadamard introdujo el concepto de problema bien planteado en los
términos siguientes: un problema está bien planteado si existe solución, es
única y depende continuamente de los datos iniciales.
Un ejemplo de Hadamard : El problema de Cauchy para la ecuación
(
uxx + uyy = 0,
u(x, 0) = uy (x, 0) = 0.

está mal planteado en el sentido de Hadamard, pues si resolvemos


uxx + uyy = 0,

u(x, 0) = 0,

uy (x, 0) = n−1 sin nx,


su solución es
1
un (x, y) = sin nx sinh ny,
n2
y no hay dependencia continua porque aunque los datos iniciales sean muy
pequeños (n muy grande), la solución tiende a infinito con n.
Otro ejemplo:
La ecuación, 4u = 0, en un dominio Ω con condiciones de tipo Neumann
∂u
= χ en ∂Ω,
∂ν
no tiene solución si la integral de χ sobre ∂Ω no es nula, pues por el teorema
de la divergencia,
Z Z Z
∂u
χ dσ = dσ = 4u dx = 0;
∂Ω ∂Ω ∂ν Ω

y si lo es, la solución no es única, se pueden añadir constantes.

Una aclaración importante: que un problema esté bien planteado o no
depende de la elección del espacio en que se busquen las soluciones. Un
problema puede tener más de una solución pero sólo una si nos restringimos
a una determinada clase de funciones. La dependencia continua siempre se
referirá a la manera en que se mida la continuidad, dependerá del espacio

14
en el que se midan los datos y de aquél donde se mida el tamaño de las
soluciones que se encuentren.
Los ejemplos anteriores muestran que el problema de Cauchy no está en
general bien planteado: no hay en general dependencia continua de los da-
tos, como muestra el ejemplo de Hadamard. Sin embargo, los problemas
con condiciones iniciales y condiciones de contorno si están en general bien
planteados para ciertas elecciones de los espacios donde se eligen los datos
y en donde se encuentran las soluciónes.

15
Capı́tulo 2

Ondas en una dimensión

2.1. Deducción de la ecuación

Tenemos una cuerda que se mueve y dos variables (x, t): (x, u(x, t)) es
la posición de la cuerda en tiempo t ≥ 0, ρ(x) es la densidad de la cuerda,
que se supone independiente del tiempo, T (x, t) denota la magnitud de la
tensión de la cuerda en (x, u(x, t)), una fuerza que actúa en la dirección de
la tangente a la posición de la cuerda en tiempo t y F (x, t) es la densidad
de la magnitud de una fuerza vertical que actúa sobre la cuerda.
Consideramos un trozo de cuerda, de x a x + h. Por la segunda ley
de Newton, la fuerza total que actúa sobre el segmento de cuerda entre
(x, u(x, t)) y (x + h, u(x + h, t)), es igual a la masa por la aceleración:

(1, ux (x + h, t)) (1, ux (x, t))

T (x + h, t) p − T (x, t) p
2
1 + ux (x + h, t) 1 + u2x (x, t)
Z x+h Z x+h
+ (0, 1) ρ(s)F (s, t) ds = (0, 1) ρ(s)utt (s, t) ds.
x x

La primera componente de esta ecuación muestra que

T (x, t)
p
1 + u2x (x, t)

es una función constante de x que supondremos igual a T (t). La segunda

ecuación implica que

T (t)
utt − c2 (x, t)uxx = F, si c2 (x, t) = .
ρ(x)

Si suponemos que los extremos de la cuerda de longitud l permanecen

fijos en todo tiempo en las posiciones (0, 0) y (l, 0), que la posición inicial de
la cuerda, u(x, 0) = f (x) y la velocidad inicial de la cuerda, ut (x, 0) = g(x),

16
son funciones conocidas, concluimos que u es la solución del problema de
Cauchy con condiciones de contorno de tipo Dirichlet



 utt − c2 (x, t)uxx = F, si 0 < x < l y t > 0,

u(x, 0) = f (x), si 0 ≤ x ≤ l,


 ut (x, 0) = g(x), si 0 ≤ x ≤ l,

u(0, t) = u(l, t) = 0, si t ≥ 0.

Si suponemos que c2 (x, t) ≡ c2 es constante, la ecuación asociada es la

ecuación de ondas unidimensional:

utt − c2 uxx = F.

2.2. Solución de d’Alembert

Buscamos la solución general de la ecuación homogénea

utt − c2 uxx = 0. (2.1)

Haciendo el cambio de variable ξ = x + ct , η = x − ct (ya hecho en clase;

sección 1.3) se llega a la ecuación

−4c2 vξη = utt − c2 uxx = 0

La solución general de esta ecuación es v(ξ, η) = ϕ(ξ)+ψ(η). Esto implica

que la de (2.1) es
u(x, t) = ϕ(x + ct) + ψ(x − ct).
Las gráficas de ϕ(x + ct) y ψ(x − ct) se dibujan desplazando las de ϕ(x)
y ψ(x) a la izquierda y derecha, respectivamente, en un cantidad ct. Se dice
que u es la suma de una onda regresiva (ϕ) y de una onda progresiva (ψ); c
es la velocidad de propagación.

2.3. Problema de Cauchy para la ecuación de on-

das
Buscamos la solución del problema de valores iniciales

2
utt − c uxx = 0 , x ∈ R, t > 0 ,

u(x, 0) = f (x) , x ∈ R ,

ut (x, 0) = g(x) , x ∈ R .


Hay que determinar ϕ y ψ de la solución general de la ecuación u(x, t) =

ϕ(x + ct) + ψ(x − ct), haciendo t = 0 en u y en ut . Tendremos, ϕ(x) + ψ(x) =

17
f (x), ϕ0 (x)+ψ 0 (x) = f (x) y ϕ0 (x)−ψ 0 (x) = 1c g(x). Resolviendo este sistema
de ecuaciones:
1 x 1 x
Z Z
1 1
ϕ(x) = f (x) + g(s) ds + k, ψ(x) = f (x) − g(s) ds − k,
2 2c 0 2 2c 0

para alguna constante k que podemos suponer siempre que es igual a cero,
lo que da

f (x + ct) + f (x − ct) 1 x+ct

Z
u(x, t) = + g(s) ds. (2.2)
2 2c x−ct

El razonamiento anterior implica la uncidad de la solución si suponemos

a priori que u ∈ C 2 (R × (0, +∞)) ∩ C 1 (R × [0, +∞)).
Si f es par o impar respecto a un punto x = l, entonces

f (x + ct) + f (x − ct)
2
es respectivamente par o impar respecto al mismo punto x = l. Lo mismo
ocurre con
1 x+ct
Z
g(s) ds,
2c x−ct
cuando g es par o impar respecto a x = l.

2.4. Dominios de dependencia y de influencia

Dependencia. La solución en (x, t) depende del valor de f (posición ini-
cial) en dos puntos, x − ct y x + ct y del valor de g (velocidad inicial) en el
intervalo, [x − ct, x + ct].
Influencia. El valor de f en x0 interviene en el valor de la solución sobre
los puntos que están en las rectas x − ct = x0 y x + ct = x0 . El valor de g
en x0 interviene en el valor de la solución sobre los puntos que están en la
región x − ct < x0 < x + ct (cono de luz ).
Suponer que f vive en el intervalo [a, b] y que g es idénticamente nula;
mostrar dónde vive la solución. Lo mismo, pero cambiando f y g.

2.5. Ecuación no homogénea

Vamos a resolver

2
utt − c uxx = F ,
 x ∈ R, t > 0 ,
u(x, 0) = 0 , x ∈ R,

ut (x, 0) = 0 , x ∈ R.


18
 Con el mismo cambio de variables,
 
ξ+η ξ−η
v(ξ, η) = u , ,
2 2c

queda  
2 ξ+η ξ−η
−4c vξη = F , .
2 2c
De u(x, 0) = 0 se deduce v(ξ, ξ) = 0 y derivando vξ (ξ, ξ) + vη (ξ, ξ) = 0. De
ut (x, 0) = 0 se deduce que vξ (ξ, ξ) − vη (ξ, ξ) = 0. Por tanto,

v(ξ, ξ) = vξ (ξ, ξ) = vη (ξ, ξ) = 0

y el dato de Cauchy de v a lo largo de la recta diagonal, ξ = η, es nulo.

Por el Teorema Fundamental del cálculo, si (ξ, η) está en el semiplano
u > v, dibujando en un plano de coordenadas (u, v) el triángulo de vértices
(ξ, η), (ξ, ξ) y (η, η), que denotamos T (ξ, η), obtenemos que
Z η
vξ (ξ, η) = vξη (ξ, v) dv,
ξ
Z ξ Z ξ Z η  Z
v(ξ, η) = vξ (u, η) du = vξη (u, v)dv du = − vξη (u, v)dudv
η η u T (ξ,η)
 
u+v u−v
Z
1
v(ξ, η) = 2 F , dudv.
4c T (ξ,η) 2 2c

Con el cambio de variables de integración, u = r + cs, v = r − cs y

recordando que, ξ = x + ct, η = x − ct, el triángulo T (ξ, η) se transforma en
el triángulo caracterı́stico T (x, t) de vértices (x, t), (x − ct, 0) y (x + ct, 0),
dudv = 2cdrds y
!
1 t
Z Z Z x+c(t−s)
1
u(x, t) = F (r, s) drds = F (r, s) dr ds,
2c T (x,t) 2c 0 x−c(t−s)

El triángulo T (x, t) se denomina triángulo caracterı́stico porque sus lados

son las rectas caracterı́sticas que pasan por (x, t), r+cs = x+ct, r−cs = x−ct
y el eje real s = 0. Este triángulo es el dominio de dependencia de u(x, t)
con respecto a los valores de F .
El valor de F en (x, t) interviene (dominio de influencia) en el cálculo de
los valores de la solución u en los puntos de la región comprendida entre las
caracterı́sticas que pasan por (x, t) y situada encima de la recta horizontal
s = t y se denomina dominio de influencia de los valores de F en (x, t).
Si F ( · , t) es par o impar respecto a l ∈ R, para todo t ≥ 0, entonces
u( · , t) es respectivamente par o impar respecto a l ∈ R, para todo t ≥ 0.

19
 La solución completa del problema de Cauchy no homogéneo

2
utt − c uxx = F (x, t) ,
 x ∈ R, t > 0 ,
u(x, 0) = f (x) , x ∈ R,

ut (x, 0) = g(x) , x ∈ R,


es
x+ct
f (x + ct) + f (x − ct)
Z
1
u(x, t) = + g(s) ds
2 2c x−ct
Z t Z x+c(t−s)
1
+ F (r, s) drds.
2c 0 x−c(t−s)

2.6. Soluciones generalizadas o débiles.

Es fácil comprobar que la solución u es C 2 (R × [0, +∞)) si f ∈ C 2 (R),
g ∈ C 1 (R) y F ∈ C 1 (R × [0, +∞)). Sin embargo, aparecen de forma natural
datos iniciales que no tienen esa regularidad, por ejemplo, cuando f es lineal
a trozos (una cuerda que tiene esquinas). En esos casos, la solución es regular
sólo en ciertas regiones separadas por rectas en las que no es C 2 (ni siquiera
C 1 a veces). Estas singularidades de u aparecen a partir de las singularidades
de f , g y F y se propagan a lo largo de las rectas caracterı́sticas que pasan
por las singularidades de f , g y F .
Estudiar la regularidad de u si f = |x|, g ≡ 0 y F ≡ 0.

2.7. Ondas en una semirrecta

1. Condición de contorno de tipo Dirichlet. Buscamos la solución del
problema 
2
utt − c uxx = 0 , x > 0, t > 0 ,



u(x, 0) = f (x) , x > 0 ,


 ut (x, 0) = g(x) , x > 0 ,

u(0, t) = h(t) , t > 0.
La solución general es u(x, t) = ϕ(x + ct) + ψ(x − ct) y como antes, las
condiciones iniciales determinan ϕ(x) y ψ(x) para x > 0. Esto da la solución
en la región x − ct > 0 con la misma expresión que en (2.2).
Haciendo x = 0, tenemos ϕ(ct) + ψ(−ct) = h(t) si t ≥ 0. Esto determina
ψ(t) para t < 0:
ψ(t) = h(−t/c) − ϕ(−t),
de donde obtenemos que para x − ct < 0,
x+ct
f (x + ct) − f (ct − x)
Z
1  x
u(x, t) = + g(s) ds + h t − .
2 2c ct−x c

20
 La regularidad de u en x − ct > 0 y en x − ct < 0 se deduce de la
de f , g y h. La regularidad sobre la recta x = ct requiere condiciones de
compatibilidad, por ejemplo, u es continua sobre la recta si f (0) = h(0). Se
demuestra que u es C 2 en el primer cuadrante si f y h son C 2 ([0, +∞)), g
es C 1 ([0, +∞)) y se verifican las condiciones de compatibilidad f (0) = h(0),
g(0) = h0 (0) y h00 (0) = c2 f 00 (0).
Cuando h es nula, la solución es la misma que la se obtiene al resolver
un problema en toda la recta con los datos iniciales f y g extendidos a
toda la recta de forma impar en 0 (Método de reflexiones) y considerando
la restricción de dicha solución al primer cuadrante; es decir, extendiendo f
y g a R de forma que f (x) = −f (−x) y g(x) = −g(−x).
2. Condición de contorno de tipo Neumann. Ahora estudiamos el pro-
blema 


utt − c2 uxx = 0 , x > 0, t > 0 ,

u(x, 0) = f (x) , x > 0 ,


ut (x, 0) = g(x) , x > 0 ,

u (0, t) = h(t) , t > 0 .
x

No hay diferencias en x − ct > 0. En x = 0 sale ϕ0 (ct) + ψ 0 (−ct) = h(t) si

t > 0, de donde, ϕ(ct) − ψ(−ct) = cH(t) − β donde H 0 (t) = h(t) y β es
constante. Luego
ψ(t) = ϕ(−t) − cH(−t/c) + β para t < 0.
La solución del problema en x − ct < 0 es
Z x+ct Z ct−x 
f (x + ct) + f (ct − x) 1
u(x, t) = + g(s) ds + g(s) ds
2 2c 0 0
Z t−x/c
−c h(s) ds + 2k + β.
0
Para que u sea continua en x = ct hay que escoger necesariamente k y
β de forma que, 2k + β = 0.
Observación. Cuando h es nula la solución es la misma que la que se
obtiene al resolver un problema en toda la recta con los datos iniciales f y g
extendidos a R como funciones pares respecto al 0 (método de reflexiones).

2.8. Ondas en una cuerda finita

Situamos una cuerda de longitud l entre los puntos 0 y l. El problema
de la vibración de la cuerda con condiciones de contorno de tipo Dirichlet es



 utt − c2 uxx = 0 , 0 < x < l, t > 0 ,

u(x, 0) = f (x) , 0 < x < l,


 ut (x, 0) = g(x) , 0 < x < l,

u(0, t) = h (t) , u(l, t) = h (t), t > 0 .
0 1

21
 Otra vez, u(x, t) = ϕ(x + ct) + ψ(x − ct) y las condiciones iniciales dan
ϕ(s) y ψ(s) cuando 0 < x < l. Con estos valores podemos escribir la solución
del problema en el triángulo {(x, t) : t > 0 y 0 < x − ct < x + ct < l}.
Usando las condiciones de contorno tenemos

ϕ(ct) + ψ(−ct) = h0 (t), t > 0,

ϕ(l + ct) + ψ(l − ct) = h1 (t), t > 0,

o bien,

ψ(s) = −ϕ(−s) + h0 (−s/c), s < 0,

ϕ(l + s) = −ψ(l − s) + h1 (s/c), s > 0.

Se ve que empezando con ϕ y ψ en el intervalo (0, l), estas dos ecuaciones

permiten definir ϕ(s) para todo s > 0 y ψ(s) para todo s < l, lo que
determina u.
Observación. Cuando h0 y h1 son nulas, la solución es la misma que la
que se obtiene al resolver un problema en toda la recta con datos iniciales
extendidos de forma impar en 0 y en l respectı́vamente (método de refle-
xiones); es decir, extendiendo f y g a R de forma que f (x) = −f (−x) y
f (l + x) = −f (l − x), g(x) = −g(−x) y g(l + x) = −g(l − x).
Una función impar respecto a 0 y l es automáticamente periódica de
periodo 2l.
Cuando las condiciones de contorno son de tipo Neumann hay que hacer
las mismas modificaciones que en el caso de la semirrecta. Si las condiciones
son nulas, las extensiones requeridas son pares en 0 y en l.
Ahora hay dos extremos y pueden haber dos condiciones de contorno
distintas, una de tipo Dirichlet y otra de tipo Neumann. Si son u(0, t) = 0
y ux (l, t) = 0, se construye la solución extendiendo f como función impar
respecto a x = 0 y par respecto a x = l, g como función par respecto a x = l
e impar respecto a x = 0 y resolviendo el problema en toda la recta.
Una función par respecto a 0 y l es automáticamente periódica de periodo
2l.
Una función par respecto a 0 e impar respecto a l es automáticamente
periódica de periodo 4l.

2.9. Conservación de la energı́a

La función Z l
E(t) = c2 u2x (x, t) + u2t (x, t) dx
0

22
es constante en [0, +∞), si u ∈ C 2 ([0, l] × [0, +∞)) es solución de

2
utt − c uxx = 0 , 0 < x < l, t > 0 ,

u(x, 0) = f (x) , 0 < x < l ,

ut (x, 0) = g(x) , 0 < x < l


que verifica condiciones de contorno de Dirichlet o Neumann nulas en la

frontera lateral del cilindro [0, l] × [0, +∞).
Este resultado implica la unicidad de solución al problema de Cauchy
asociado a la ecuación de ondas, con condiciones de contorno de tipo Diri-
chlet o tipo Neumann y en la clase de soluciones C 2 ([0, l] × [0, +∞)).
Este resultado también es cierto si la cuerda vibrante se remplaza por
una membrana, es decir, para otras dimensiones n del espacio euclı́deo en
el que se encuentre la membrana Ω ⊂ Rn . En este caso se utiliza el teorema
de la divergencia para mostrar que
Z
E(t) = c2 |∇u|2 + u2t dx
Ω

es constante si utt − c2 ∆u = 0 en Ω × [0, +∞) ⊂ Rn+1 , u verifica condiciones

de contorno de Dirichlet o Neumann nulas en la frontera lateral de Ω y
u ∈ C 2 (Ω × [0, +∞)).

Demostración.
Z
0
E (t) = 2 c2 ∇u · ∇ut + ut utt dx
ZΩ Z
2 2 2


= 2 c ∇ · (ut ∇u) + ut utt − c 4u dx = 2c ut ∇u · ν dσ,
Ω ∂Ω

que es cero si u o su derivada normal son nulas la frontera lateral ∂Ω ×

(0, +∞).

23
Capı́tulo 3

La ecuación del calor en la

recta

3.1. Deducción de la ecuación

Sea Ω ⊂ R3 un sólido y u(x, t) la temperatura en uno de sus puntos x
en tiempo t > 0. Si D es una parte pequeña de Ω, la cantidad de calor en
D en tiempo t es Z
Q(t) = ρu(x, t) dx ,
D
donde ρ es la densidad de Ω. La variación total del calor en D es
Z
0
Q (t) = ρut (x, t) dx.
D

Las leyes de la termodinámica dicen que el flujo del calor M , un campo

vectorial, va de las regiones más calientes a las más frias y que es proporcional
al gradiente de la temperatura: M = −A∇u, donde A es el calor especı́fico.
La variación total del calor en D coincide con la cantidad neta de calor
que entra y sale de D, junto los cambios de calor que son generados directa-
mente sobre D por cambios de temperatura externos con densidad F (x, t);
es decir, si ν es la normal unitaria exterior a ∂D,
Z Z Z
ρut (x, t) dx = − M · ν dσ + ρF (x, t) dx,
D ∂D D
Z Z Z
ρut (x, t) dx = A ∇u · ν dσ + ρF (x, t) dx
D ∂D D
y por el teorema de la divergencia
Z
[ρut (x, t) − A4u(x, t) − ρF (x, t)] dx.
D

24
Como D es arbitrario, la temperatura verifica la ecuación
A
ut − c2 4u = F en D × (0, +∞), c2 =
ρ
y con el cambio de variables, v(x, t) = u(cx, t), la ecuación se transforma en

vt − 4v = G(x, t), G(x, t) = F (cx, t),

que se denomina la ecuación del calor.

Queremos resolver el problema de valores iniciales no homogéneo:
(
ut − uxx = F, x ∈ R, t > 0,
(3.1)
u(x, 0) = f (x), x ∈ R.

“Lo natural” es buscar la solución en

C 2 (R × (0, +∞)) ∩ C(R × [0, +∞)) ∩ L∞ (R × (0, +∞)),

si f ∈ C(R) ∩ L∞ (R) y F ∈ C ∞ (R × [0, +∞)).

3.2. Soluciones autosemejantes

Para resolver el problema de valores inciales no homogéneo (3.1), se
empieza por intentar resolver el problema de valores iniciales homogéneo
(
uxx − ut = 0, en R × (0, +∞),
(3.2)
u(x, 0) = f (x), en R.

Para ello, buscamos primero las soluciones autosemejantes de la ecua-

ción, ut − uxx = 0 y consideraremos “sumas”de estas soluciones, para obte-
ner más soluciones. Después, intentamos determinar los “coeficientes de la
suma” para que se verifiquen las condiciones iniciales.
Si u(x, t) verifica, ut − uxx = 0 en R × (0, +∞), también lo verifican
uλ (x, t) = u(λx, λ2 t), para todo λ real y u(x + x0 , t + t0 ), para todo (x0 , t0 )
en R2 .
¿Existen soluciones autosemejantes, es decir, tales que

u(λx, λ2 t) = λα u(x, t),

√
para algún α y todo
√ λ > 0? Si esto ocurre, haciendo λ = 1/ t, tenemos que
u(x, t) = tα/2 u(x/ t, 1), es decir,
√
u(x, t) = tα/2 Φ(x/ t), si Φ(x) = u(x, 1).

A las soluciones de este tipo se las llama soluciones autosemejantes de

la ecuación.

25
 √
Si u(x, t) = tα/2 Φ(x/ t) es una solución autosemejante, Φ tiene que
verificar
s α
Φ00 (s) + Φ0 (s) − Φ(s) = 0 (3.3)
2 2
y si pedimos a la solución que conserve el calor total, es decir, que u( · , t)
sea integrable en R para todo t > 0 y que la integral en x de u no dependa
de t > 0, Z Z
α+1
u(x, t) dx = t 2 Φ(y) dy,
R R
resulta que α debe ser igual a −1. En este caso, (3.3) se reduce a
s
(Φ0 + Φ)0 = 0,
2
cuya solución general es
2
Z s τ2 s2
− s4
Φ(s) = αe e 4 dτ + βe− 4 , si α, β ∈ R.
0

De estas sólo la segunda es integrable en R,

Z +∞ 2
Z s τ2
Z +∞ Z 1 s2
− s4 2
e e 4 dτ ds = se− 4 (1−τ ) dτ ds
0 0 0 0
Z 1 Z +∞ Z 1
s2 2 dτ
= se− 4 (1−τ ) dsdτ = 2 = +∞
0 0 0 1 − τ2
2 /4
y nos quedamos con la solución, Φ(s) = e−s . Ası́
2 /4t
u(x, t) = βt−1/2 e−x

y elegimos la constante β para que la integral en x de v(x, t) sea 1,

1 −x2 /4t
u(x, t) = √ e .
4πt
Esta función verifica, ut − uxx = 0, en R × (0, +∞). Si la trasladamos en
la dirección espacial en una dirección dada, y ∈ R, es decir, consideramos
u(x − y, t), también verifica la ecuación del calor en R × (0, +∞).
La función
2
(
√ 1 e−|x| /4t , si t > 0,
G(x, t) = 4πt
0, si t ≤ 0,
se denomina el núcleo de Gauss o núcleo Gaussiano.
Escogiendo para cada y un valor f (y) y “sumando sobre todos los valores
de y”, obtenemos Z
u(x, t) = G(x − y, t)f (y) dy, (3.4)
R

26
que es una suma formal de soluciones y quizás una solución del problema
de valores iniciales homogéneo (3.2).
Cada término del integrando en la fórmula anterior es solución de la
ecuación del calor en las variables (x, t). Para que u sea solución, necesitamos
primero, que la fórmula (3.4) tenga sentido y después, que se pueda derivar
bajo el signo integral. Usando el teorema de la convergencia dominada se
puede probar que esto es ası́ si f ∈ L∞ (R). Lo mismo es cierto si f ∈ Lp (R)
y 1 ≤ p ≤ ∞. Lo último se deduce fácilmente a partir de las acotaciones en
el Teorema 4.

3.3. La solución al problema de valores iniciales

¿Cuánto vale
√ la solución anterior si t = 0? Haciendo el cambio de varia-
bles, x − y = t z, resulta que

1
Z +∞ √ y2
u(x, t) = √ f (x − t y)e− 4 dy.
4π −∞
De esta fórmula, de la acotación,
√ y2 y2
|f (x − t y) e− 4 | ≤ kf kL∞ (R) e− 4

y del teorema de la convergencia dominada es fácil deducir que

lı́m u(x, t) = f (x),

t→0+

si f es continua en x y acotada en R.

Teorema 4. Si f ∈ C(R) ∩ L∞ (R), la función definida por (3.4) verifica

que, u ∈ C ∞ (R×(0, +∞))∩C(R×[0, +∞))∩L∞ (R×[0, +∞)) y es solución
de (3.2).
1 s 2
Demostración. Como G(x, t) = t− 2 Φ( √xt ), Φ(s) = √1 e− 4
4π
y ∂t G − ∂x2 G =
0, en R × (0, +∞), se verifica que

∂xα ∂tβ G(x, t) = ∂xα+2β G(x, t)

1+α+2β x 1+α+2β x x2
= t− 2 Φ(α+2β) ( √ ) = t− 2 Pα+2β ( √ )e− 4t ,
t t
donde Pα+2β es un polinomio de grado α + 2β. De estas fórmulas, se deduce
que para todo  > 0 y R > 0 existe N = N (, R, α, β) > 0 tal que
2 /N
|∂xα ∂tβ G(x − y, t)| ≤ N e−y , si |x| ≤ R y  ≤ t ≤ R.

27
Esto y el teorema de convergencia dominada muestran que u está en C ∞ (R×
(0, +∞)) si f ∈ L∞ (R). Además, como
Z Z
G(x − y, t) dy = G(z, t) dz = 1, si t > 0, (3.5)
R R

se verifica que
|u(x, t)| ≤ kf kL∞ (Rn )
y u está acotada en el semiplano superior.
Para mostrar la continuidad de u en R×{0} se utiliza la continuidad de f
en R, el decaimiento en el infinito y la propiedad (3.5) del nucleo Gaussiano:
fijado x0 ∈ R y dado  > 0, existe un δ > 0 tal que |f (y) − f (x0 )| ≤ , si
|y − x0 | < δ y
Z
|u(x, t) − f (x0 )| ≤ G(x − y, t)|f (y) − f (x0 )| dy
Bδ (x0 )
Z
+ G(x − y, t)|f (y) − f (x0 )| dy
Rn \Bδ (x0 )
Z
≤  + 2kf k∞ G(x − y, t) dy, (3.6)
Rn \Bδ (x0 )

para todo (x, t) en el semiplano superior. Por otro lado, δ ≤ |y−x0 | ≤ 2|y−x|,
si |x − x0 | ≤ δ/2 y |y − x0 | ≥ δ, de donde
1 |x−y|2 |y−x|2 1 δ2 |y−x|2
G(x − y, t) ≤ (4πt)− 2 e− 8t e− 8t ≤ (4πt)− 2 e− 32t e− 8t

y Z
δ2
G(x − y, t) dy ≤ N e− 32t , si |x − x0 | ≤ δ/2 . (3.7)
Rn \Bδ (x0 )

De (3.6) y (3.7) se deduce que

δ2
|u(x, t) − f (x0 )| ≤  + 2N kf kL∞ (R) e− 32t , si |x − x0 | ≤ δ/2

y que u es continua en (x0 , 0).

Algunas propiedades de la solución obtenida:

• Si m ≤ f (x) ≤ M , entonces m ≤ u(x, t) ≤ M para todo x ∈ R y t > 0,
es decir, las temperaturas máxima y mı́nima originales no se pueden superar
ni por arriba ni por abajo.

Demostración.
R Escribir la fórmula para u(x, t), utilizar que m ≤ f ≤ M y
que R G(z, t) dz = 1, si t > 0.

28
 R R
• Conservación del calor : R u(x, t) dx = R f (x) dx, para todo t > 0, si
además f está en L1 (R).

Demostración. Si f tiene soporte compacto, u( · , t) y sus derivadas decaen

hacia cero si x tiende hacia ±∞ y
Z Z
∂t u(x, t) dx = ∂x2 u(x, t) dx = 0.
R R

El caso general se obtiene aproximando f con funciones de soporte com-

pacto. También se puede utilizar el siguiente razonamiento que justifica el
teorema de Fubini-Tonelli:
Z Z Z  Z
u(x, t) dx = f (y) G(x − y, t) dx dy = f (y) dy.
R R R R

• Si además f ∈ L2 (R), entonces

Z Z T Z Z
2 2
u (x, T ) dx + 2 |∂x u(x, t)| dxdt = f 2 (x) dx, (3.8)
R 0 R R

para todo T > 0.

Demostración. Si f tiene soporte compacto, u( · , t) y sus derivadas decaen

hacia cero si x tiende hacia ±∞ y

0 = u ∂t u − ∂x2 u = 12 ∂t u2 + (∂x u)2 − ∂x (u∂x u).

La identidad es consecuencia de integrar la fórmula anterior en R×(0, T ).

• El análogo de la identidad anterior para cilindros acotados es la si-

guiente fórmula: si u ∈ C 2 (Ω × (0, +∞)) ∩ C(Ω × [0, +∞)) verifica

∂t u − 4u = 0, en Ω × (0, +∞),

u(x, 0) = f (x) en Ω,

u o ∂u
∂ν = 0, en ∂Ω × (0, T ).


Entonces,
Z Z T Z Z
2 2
u (x, T ) dx + 2 |∇u| dxdt = f 2 (x) dx, (3.9)
Ω 0 Ω Ω

para todo T > 0. Las identidades (3.9) y (3.8) se conocen como la identidad
de energı́a para la ecuación del calor.

29
Demostración. Se integra sobre Ω × (0, T ) la identidad

0 = u (∂t u − 4u) = 12 ∂t u2 + |∇u|2 − ∇ · (u∇u) .

y se aplica el teorema de la divergencia para concluir que

Z Z
∇ · (u∇u) dx = u ∂u
∂ν dσ = 0
Ω ∂Ω
y
Z Z T Z 
2 2
∂T u (x, T ) dx + 2 |∇u(x, t)| dxdt = 0, para todo T > 0.
Ω 0 Ω

3.4. La unicidad
• La función
x 2
∂x G(x, t) = − √ e−x /4t
2 4πt 3

cumple la ecuación del calor en t > 0 y lı́mt→0+ ∂x G(x, t) = 0, para todo

x. Si entendemos la condición inicial en este sentido, tenemos una solución
no nula del problema (3.2) y con f nula. Por tanto, no habrı́a unicidad.
Ahora bien, si entendemos la condición inicial como un lı́mite en las dos
variables, esta solución ya no valdrı́a, pues al acercarnos al (0, 0) a lo largo
de la parábola x2 = t, el lı́mite es infinito.
• Ejemplo de Tychonoff. La función
∞
X x2k 2
u(x, t) = g (k) (t) , g(t) = e−1/t , g(0) = 0, (3.10)
(2k)!
k=0

verifica la ecuación del calor en R2 , u ∈ C ∞ (R2 ) y u( · , 0) ≡ 0. Además,

existe N ≥ 1 tal que
N x2
− 12
|u(x, t)| ≤ N e t Nt , si x ∈ R y t > 0; (3.11)

es decir, u crece “demasiado” en el infinito en cada tiempo fijo, t > 0.

Demostración. Con la fórmula de representación de Cauchy para las deri-

vadas de una función analı́tica mostraremos que existe N ≥ 1 tal que
1
|g (k) (t)| ≤ N 1+k k!t−k e− N t2 , si t > 0 y k ≥ 0. (3.12)

Esta acotación y la fórmula de Stirling,

k!
lı́m √ = 1,
k k
k→+∞


2πk e

30
implican que el término general de la serie (3.10) verifica
2k  2 k
   
x N N x 1
e− N t2 , si x ∈ R, t > 0 y k ≥ 0,
 (k)
g (t) ≤
 (2k)!  k! t

que la serie (3.10) converge uniformemente sobre compactos de R × (0, +∞),

que su suma u verifica (3.11), es continua en R2 y u( · , 0) ≡ 0. Lo mismo es
cierto para cualquier derivada de u.
Para ver que (3.12) es cierto, dado t > 0 y  en (0, 1),
2
e−1/z
Z
(k) k!
g (t) = dz,
2πi |z−t|=t (z − t)k+1

por la fórmula de Cauchy y elegimos  > 0 pequeño de forma que Bt (t)
esté inscrito en un cono con vértice en (0, 0) y con ángulo medio de apertura
igual a π/8. Entonces, existe N ≥ 1 tal que si z = reiθ , se verifica
1
< 1/z 2 = r−2 cos 2θ ≥ , si |z − t| = t, t > 0
N t2
y
1
e− N t2
Z
k! 1
|g (k)
(t)| ≤ |dz| ≤ N 1+k k!t−k e− N t2 .
2π |z−t|=t (t)k+1

De lo anterior concluimos que la ecuación del calor en R × (0, +∞) no

tiene solución única sin exigir condiciones adicionales a la solución.
En los resultados siguientes veremos algunas condiciones adicionales que
garantizan la unicidad de solución para el problemas de valores iniciales
sobre un intervalo de R, un dominio acotado de Rn , la recta real R y el
espacio euclı́deo Rn .

Teorema 5 (Principio del máximo débil). Sean ΩT = (a, b) × (0, T ] y ΓT =

[a, b] × {0} ∪ {a, b} × (0, T ], la frontera parabólica del cilindro ΩT y 0 < T <
+∞. Si u ∈ C(ΩT ) ∩ C 2 (ΩT ) y ∂xx u − ∂t u ≥ 0 en ΩT , entonces

máx u ≤ máx u.
ΩT ΓT

Análogamente hay un principio del mı́nimo débil: aplicar el Teorema 5

a −u, si ∂xx u − ∂t u ≤ 0 en ΩT y concluir que

mı́n u ≥ mı́n u.
ΩT ΓT

Demostración. Suponemos primero que uxx − ut > 0 en ΩT . Si u tiene

un máximo local en algún (x0 , t0 ) que es interior a ΩT , se verifica que
ut (x0 , t0 ) = 0 y uxx (x0 , t0 ) ≤ 0. Si hubiera un máximo local en (x0 , T ),

31
para algún a < x0 < b, se tendrı́a que, ut (x0 , T ) ≥ 0 y uxx (x0 , T ) ≤ 0 y
ambos casos son incompatibles con la hipótesis, uxx − ut > 0 en ΩT .
Si uxx −ut ≥ 0 en ΩT , definimos u (x, t) = u(x, t)+x2 y por el resultado
anterior
máx u ≤ máx u ≤ máx u ≤ máx u +  máx {a2 , b2 },
ΩT ΩT ΓT ΓT

lo que implica el resultado si hacemos  tender a 0.

• Lo mismo es cierto en Rn si ΩT = Ω×(0, T ], ΓT = Ω×{0}∪∂Ω×(0, T ],
4u−∂t u ≥ 0 en ΩT , Ω es un abierto acotado de Rn y n ≥ 2. En este caso, en
un posible máximo interior (x0 , t0 ) a ΩT , la matriz Hessiana de u en (x0 , t0 )
es semidefinida negativa y su traza, 4u(x0 , t0 ), es menor o igual que cero.
• Consecuencia: unicidad de solución al problema de valores iniciales
asociado a la ecuación del calor con condiciones de Dirichlet en la frontera
lateral cuando la base es un intervalo o un abierto acotado en Rn y para
posibles soluciones en C 2 (Ω × (0, T ]) ∩ C(ΩT ): para verificarlo aplicar el
principio del máximo débil a ± (u1 − u2 ) en ΩT , si u1 y u2 son soluciones en
C 2 (Ω × (0, T ]) ∩ C(ΩT ) de

4u − ∂t u = F, en Ω × (0, T ],

u(x, 0) = f (x), en Ω,

u = 0 en ∂Ω × (0, T ],


para concluir que u1 = u2 en Ω × [0, T ].

Teorema 6 (Principio del máximo débil para R). Sea u ∈ C(R × [0, T ]) ∩
C 2 (R × (0, T ]) y tal que ∂xx u − ∂t u ≥ 0 en R × (0, T ], u(x, t) ≤ M en
R × [0, T ], para algún M ≥ 0 y u(x, 0) = f (x) en R. Entonces
u(x, t) ≤ sup f, si x ∈ R y 0 ≤ t ≤ T.
R
Demostración. Fijado y ∈ R y 0 < µ ≤ 1, consideramos la función
vµ (x, t) = u(x, t) − µ |x − y|2 + 2t .

Se verifica que, ∂xx vµ − ∂t vµ ≥ 0 en R × (0, T ] y por el principio del

máximo débil para cilindros acotados
vµ (y, t) ≤ máx vµ , si Γρ = [y − ρ, y + ρ] × {0} ∪ {y − ρ, y + ρ} × [0, T ]
Γρ

y 0 ≤ t ≤ T . En la parte lateral de Γρ
|x − y|2 + 2t ≥ ρ2 y vµ (x, t) ≤ M − µρ2 ≤ sup f,
R

si ρ ≥ ρ(M, supR f, µ). Por tanto,

vµ (y, t) = u(y, t) − 2µt ≤ sup f, si 0 ≤ t ≤ T
R

y haciendo, µ → 0+ , obtenemos el resultado.

32
 Lo mismo funciona en Rn × (0, T ] si cambiamos vµ por

vµ (x, t) = u(x, t) − µ |x − y|2 + 2nt .

• Consecuencia: unicidad de solución para el problema de condiciones

iniciales asociado a la ecuación del calor en R × (0, +∞) en la clase de
soluciones C 2 (R × (0, +∞) ∩ C(R × [0, +∞)) ∩ L∞ (R × [0, +∞). Si f ∈
L∞ (R) ∩ C(R), la solución de (3.2) generada con el núcleo de Gauss está en
C 2 (R × (0, +∞) ∩ C(R × [0, +∞)) ∩ L∞ (R × [0, +∞) y es única.
• La conclusión del Teorema 6 también es válida si la condición, u(x, t) ≤
2
M , se cambia por, u(x, t) ≤ M eM |x| en R × [0, T ], para algún M > 0
(Comparar esto con las soluciones de Tychonoff).

3.5. Ecuación no homogénea: método de Duhamel

Buscamos una solución para el problema no homogéneo
(
ut − uxx = F (x, t), x ∈ R, t > 0,
(3.13)
u(x, 0) = 0, x ∈ R.

La ecuación diferencial ordinaria no homogénea con condición inicial

(
y 0 + ay = F (t),
y(0) = 0,

se resuelve por el método de variación de las constantes, que es equivalente

al método de Duhamel que describimos a continuación: la solución y(t) se
puede escribir como Z t
y(t) = γ(t, s) ds.
0
Entonces,
Z t Z t
0 0
y (t) = γ(t, t) + γt (t, s) ds, y + ay = γ(t, t) + γt (t, s) + aγ(t, s) ds
0 0

y si queremos que y sea solucion, esto se cumplirá si elegimos γ de forma

que (
γt (t, s) + aγ(t, s) = 0, en (s, +∞)
γ(s, s) = F (s), en t = s,
que es un problema homogéneo que sabemos resolver. Si w(t, s) es la solución
de la ecuación homogénea
(
wt (t, s) + aw(t, s) = 0, en (0, +∞)
w(0, s) = F (s), en t = 0,

33
w(t, s) = e−at F (s) y γ(t, s) = w(t − s, s); es decir.
Z t
y(t) = w(t − s, s) ds.
0

Para resolver (3.13), procedemos de forma similar

Teorema 7. Si w( · , · , s) es la solución de
(
wt − ∂x2 w = 0, x ∈ R, t > 0,
w(x, 0, s) = F (x, s), x ∈ R,

entonces
Z t Z tZ
u(x, t) = w(x, t − s, s) ds = G(x − y, t − s)F (y, s) dyds (3.14)
0 0 R

es solución de (3.13).

Demostración. Para deducir tal fórmula, escribimos

Z t
u(x, t) = γ(x, t, s) ds
0

y un cálculo formal muestra que

Z t
∂t u − ∂x2 u = γ(x, t, t) + ∂t − ∂x2 γ(x, t, s) ds.


0

Lo lógico es elegir γ de forma que

(
∂t − ∂x2 γ(x, t, s) = 0,


en R × (s, +∞),
γ(x, s, s) = F (x, s), en R × {s}

y como la ecuación es invariante por traslaciones,

γ(x, t, s) = w(x, t − s, s).

Finalmente, Z
w(x, t, s) = G(x − y, t)F (y, s) dy.
R

El principio de Duhamel se puede aplicar también a la ecuación de ondas

no homogénea.

34
Teorema 8. Sea w( · , · , s) la solución de

2 2 2
∂t w − c ∂x w = 0,
 x ∈ R, t > 0,
w(x, 0; s) = 0, x ∈ R;

∂t w(x, 0; s) = F (x, s), x ∈ R,


entonces Z t
u(x, t) = w(x, t − s, s) ds (3.15)
0
es solución de la ecuación de ondas

2 2 2
∂t u − c ∂x u = F (x, t),
 x ∈ R, t > 0,
u(x, 0) = 0, x ∈ R, (3.16)

ut (x, 0) = 0, x ∈ R.


Podéis comprobar que esto coincide con los resultados obtenidos ante-
riormente para esta ecuación.

Demostración. El lector puede comprobrar que (3.15) es formalmente una

solución de (3.16). Para deducir la fórmula, escribimos
Z t
u(x, t) = γ(x, t, s) ds
0

y un cálculo formal muestra que

Z t
∂t2 u c2 ∂x2 u ∂t2 − c2 ∂x2 γ(x, t, s) ds,


− = ∂t γ(x, t, t) +
0

si suponemos que γ(x, t, t) = 0. Lo lógico es elegir γ de forma que

(
∂t2 − c2 ∂x2 γ(x, t, s) = 0,


en R × (s, +∞),
γ(x, s, s) = 0, ∂t γ(x, s, s) = F (x, s), en R × {s}

y como la ecuación es invariante por traslaciones,

γ(x, t, s) = w(x, t − s, s).

3.6. La ecuación del calor en un cilindro

• La solución del problema
(
ut − uxx = F (x, t), en R × (0, +∞),
u(x, 0) = f (x), en R,

35
es par o impar respecto a x = l en todo tiempo t > 0, si la misma propiedad
es cierta para f y F ( · , t), para todo t > 0. Por ello, podemos utilizar el
método de reflexiones para encontrar la solución del problema

ut − uxx = F (x, t), x > 0, t > 0,

u(x, 0) = f (x), x > 0,

u(0, t) = 0, t > 0,


resolviendo un problema en todo R que tiene por datos las extensiones im-
pares a R de f y F ( · , t). Si F ≡ 0, la solución es
Z ∞
1 h 2 2
i
u(x, t) = √ f (y) e−(x−y) /4t − e−(x+y) /4t dy.
4πt 0
• Si la condición de contorno es de tipo Neumann, ux (0, t) = 0, la ex-
tensión debe ser par.
• El método de reflexiones también funciona para la ecuación del calor
con condiciones de Dirichlet, Neumann o con condiciones mixtas sobre un
intervalo finito. Para ello, se calculan las correspondientes extensiones de f y
F ( · , t) a todo R como funciones pares o impares respecto a los extremos del
intervalo, se resuelve el problema asociado a las extensiones de f y F ( · , t)
sobre R y se restringe la solución calculada al cilindro asociado al intervalo.
• El método de Duhamel funciona de forma análoga en un cilindro de
base finita o de media recta.

36
Capı́tulo 4

Separación de variables

4.1. Calor en una barra finita: método de Fourier

• Consideramos el problema

ut − uxx = 0 ,
 0 < x < π,t > 0,
u(x, 0) = f (x) , 0 < x < π, (4.1)

u(0, t) = u(π, t) = 0 , t > 0.


El método de Fourier (separación de variables) consiste en buscar solu-

ciones no nulas en variables separadas de la ecuación, u(x, t) = X(x)T (t),
que cumplan las condiciones de contorno laterales. Después se intenta que
una “combinación lineal”o suma infinita de éstas verifique la condición ini-
cial y sea la solución de la ecuación.
Si u(x, t) = X(x)T (t) es una solución no nula en [0, π] de ut − uxx = 0,
se debe cumplir que, X(x)T 0 (t) − X 00 (x)T (t) = 0; es decir,

X 00 (x) T 0 (t)
= ,
X(x) T (t)

son funciones constantes y podemos suponer que para algún λ ∈ R,

X 00 (x) T 0 (t)
= = −λ.
X(x) T (t)

Las condiciones de contorno exigen que X(0) = X(π) = 0, de modo

que necesitamos encontrar las soluciones no nulas del problema (llamado de
Sturm-Liouville) (
X 00 + λX = 0 en [0, π],
(4.2)
X(0) = X(π) = 0.
Es fácil comprobar que (4.2) tiene soluciones no nulas si y sólo si, λ = k 2
(valores propios), k = 1, 2, . . . y que éstas soluciones son múltiplos de sin kx

37
(funciones propias). Para estos valores de λ, las correspondientes soluciones
no nulas de la ecuación,
T 0 + λT = 0,
2t
son múltiplos de e−k y las soluciones de variables separadas de (4.1), son
2t
uk (x, t) = e−k sin (kx), si k ≥ 1.

Para resolver (4.1) y suponiendo que no hubiera problemas para la con-

vergencia y la conmutación de las derivadas con la suma infinita, la ecuación
y las condiciones de contorno se cumplen si tomamos por solución de (4.1),
la suma infinita
∞
2
X
u(x, t) = bk e−k t sin kx. (4.3)
k=0

Para la condición inicial necesitarı́amos (formalmente) que

∞
X
bk sin kx = f (x) (4.4)
k=1

y la cuestión que surge es si existen coeficientes bk que hagan posible la

identidad.
• Una familia ortogonal de funciones {ϕk : k ≥ 1} es completa en L2 (a, b)
si para toda f en L2 (a, b) existen números {αk } tales que,
N
X
kf − αk ϕk kL2 (a,b)
k=1

tiende a cero, si N → +∞. En tal caso,

Z b Z b Z b
f ϕl = (f − TN f ) ϕl dx + TN f ϕl dx
a a a
Z b Z b
= (f − TN f ) ϕl dx + αl ϕ2l ,
a a
PN
si N ≥ l ≥ 1 y TN f = k=1 αk ϕk . Por la desigualdad de Hölder
Z b 
 

 (f − TN f ) ϕl dx ≤ kf − TN f kL2 (a,b) kϕl kL2 (a,b)
a

que tiende a cero si N crece. Es decir,

Z b Z b
αl = f ϕl dx/ ϕ2l dx, si l ≥ 1. (4.5)
a a

38
 Si {ϕk } es completa, la serie +∞
P
k=1 αk ϕk se denomina serie de Fourier de
f en la base {ϕk }. La ortogonalidad de {ϕk } y (4.5) implican la desigualdad
de Bessel para {ϕk }
N
X N
X Z b
kf − αk ϕk k2L2 (a,b) = kf k2L2 (a,b) − αk2 ϕ2k dx ≥ 0, si N ≥ 1 (4.6)
k=1 k=1 a

y pasando al lı́mite
+∞
X Z b
αk2 ϕ2k dx ≤ kf k2L2 (a,b) .
k=1 a

La completitud de {ϕk } y (4.6) dan la identidad de Parseval o de Plan-

cherel para {ϕk },
+∞
X Z b
kf k2L2 (a,b) = αk2 ϕ2k dx, si f ∈ L2 (a, b).
k=1 a

4.2. Criterio de Weierstrass

Teorema 9 (Criterio de Weierstrass). Si {fk } es una sucesión de funciones,
fk : Ω −→ R, Ω es una abierto en Rn y para cada K ⊂ Ω compacto de Ω
existe una sucesión {Mk } de números positivos tales que
+∞
X
|fk | ≤ Mk en K y Mk < +∞.
k=1
P+∞
Entonces, la serie de funciones, k=0 fk , converge uniformemente sobre
compactos de Ω. Además, si fk ∈ C(Ω), para todo k ≥ 1, entonces f ∈ C(Ω).
Demostración. Si SN (x) = +∞
P
k=1 fk (x) y 1 ≤ N < M ,

M
X
|SN (x) − SM (x)| ≤ Mk , si x ∈ K,
k=N +1

suma que será más pequeña que  > 0, si N ≥ N . Es decir, {SN } es una
sucesión de Cauchy en L∞ (K) y la serie converge uniformemente sobre K.
La suma de la serie será continua en Ω, por ser lı́mite uniforme de funciones
continuas sobre compactos de Ω, si cada fk ∈ C(Ω).

Teorema 10 (Derivabilidad de una serie de funciones). Si además, fk ∈

C 1 (Ω), para todo k ≥ 1 y para cada K ⊂ Ω compacto existe una sucesión
Nk , tal que
+∞  
X  ∂fk 
Nk < ∞ y   ≤ Nk en K,
∂xi 
k=0

39
para cada i = 1, · · · , n y k ≥ 1. Entonces, f ∈ C 1 (Ω) y
+∞
∂f X ∂fk
= , para i = 1, . . . , n.
∂xi ∂xi
k=0

La demostración la podéis encontrar en el libro J.E. Marsden y M. J.

Hoffman, Análisis Clásico Elemental.

4.3. Series de Fourier

Serie trigonométrica de periodo 2π
• Es toda serie de la forma
∞
a0 X
+ (ak cos kx + bk sin kx).
2
k=1

• Al haber un paso al lı́mite, que es la suma de la serie, esta puede

entenderse de distintas formas: convergencia puntual, uniforme, en normas
variadas. . .

Coeficientes de Fourier
• De la ortogonalidad de la familia {1, cos kx, sin kx} en [−π, π], si f
coincide con la suma de su serie de Fourier y la serie converge uniformente
(o en L2 (−π, π)) a f , se obtiene que

1 π 1 π
Z Z
ak = f (x) cos kx dx y bk = f (x) sin kx dx, si k ≥ 0
π −π π −π

y se denominan los coeficientes de Fourier de f en [−π, π]. Estos tienen

sentido y se pueden calcular si f está en L1 (−π, π).
• La serie trigonométrica construida con estos coeficientes se llama la
serie de Fourier de f .
• Las propiedades elementales de los coeficientes de Fourier son

Linealidad: ak (f + g) = ak (f ) + ak (g), bk (f + g) = bk (f ) + bk (g).

Acotación: |ak |, |bk | ≤ kf k1 /π, si k ≥ 0.

Lema de Riemann Lebesgue: Si f ∈ L1 (−π, π), entonces

lı́m ak = lı́m bk = 0.
k→∞ k→∞

40
Demostración. Si f = χ[a,b] , [a, b] ⊂ [−π, π],
1 b sin kb − sin ka
Z
ak (f ) = cos kx dx = ,
π a k
que tiende a cero, si k → +∞. Si f ∈ C([−π, π]), para cada  > 0 existe
una partición, −π
P =−1x0 < x1 < . . . < xN = π, de forma que la función
escalonada, S = N i=0 f (xi )χ[xi ,xi+1 ) verifica, kf − SkL∞ (−π,π) ≤ 2 y

|ak (f )| ≤ |ak (f − S)| + |ak (S)| ≤  + |ak (S)|, para todo k ≥ 1.

Como lı́mk→+∞ ak (S) = 0, podremos encontrar k tal que |ak (f )| ≤ 2,
si k ≥ k . En general, si f ∈ L1 (−π, π), existe g ∈ C([−π, π]), tal que
kf − gkL1 (−π,π) ≤ π y
|ak (f )| ≤ |ak (f − g)| + |ak (g)| ≤  + |ak (g)|, para todo k ≥ 1.
Como lı́mk→+∞ ak (g) = 0, podremos encontrar k tal que, |ak (f )| ≤ 2, si
k ≥ k .
• Si f ∈ C 1 ([−π, π]) y f (π) = f (−π), la relación entre los coeficientes
de Fourier de f y f 0 es la siguiente:
ak (f 0 ) = kbk (f ) , bk (f 0 ) = −kak (f ), si k = 1, 2, 3 . . .
• Si f ∈ C ∞ ([−π, π]) y sus derivadas son periódicas de periodo 2π, los
coeficientes de Fourier de f verifican
2
|ak (f )| + |bk (f )| ≤kf (N ) kL1 (−π,π) , para todo k ≥ 0 y N ≥ 1.
kN
Demostración. Integrar por partes.

Núcleo de Dirichlet
• Si SN (f )(x) es la N -ésima suma parcial de la serie de Fourier de una
función f en L1 (−π, π),
N
" #
1 π
Z
1 X
SN (f )(x) = f (t) + cos kx cos kt + sin kx sin kt dt
π −π 2
k=1
N
" #
1 π
Z
1 X
= f (t) + cos k(x − t) dt.
π −π 2
k=1

La fórmula, 2 cos x sin y = sin (x + y) − sin (x − y), muestra que

N
( )
1 X y 
2 + cos ky sin
2 2
k=1
y  X N        
1 1 1
= sin + sin k + y − sin k − y = sin N + y
2 2 2 2
k=1

41
y
Z π
1 sin(N + 1/2)x
SN (f )(x) = f (t)DN (x − t) dt, donde DN (x) = .
π −π 2 sin (x/2)
La función DN se denomina núcleo de Dirichlet. Es una función regular
y periódica de periodo 2π.
• Extendiendo f a todo R como una función periódica de periodo 2π y
haciendo el cambio de variables natural, podemos escribir
1 π 1 π
Z Z
SN (f )(x) = f (x − t)DN (t) dt = f (x + t)DN (t) dt. (4.7)
π −π π −π
• Como SN (1) ≡ 1, (4.7) implica que
1 π
Z
DN (t) dt = 1. (4.8)
π −π

Resultados de convergencia
A partir del lema de Riemann-Lebesgue y de la representación (4.7) con
el núcleo de Dirichlet podemos probar el siguiente resultado de convergencia
puntual de la serie de Fouriter:
Teorema 11. Si f : [−π, π] −→ R verifica, f (π) = f (−π) y para algún
α > 0 y M > 0,
|f (x) − f (y)| ≤ M |x − y|α , si x e y ∈ [−π, π]. (4.9)
Entonces, lı́mN →∞ SN (f )(x) = f (x) en [−π, π].
Demostración. La extensión periódica de f a todo R verifica (4.9) y por
(4.8), podemos escribir
1 π
Z
SN (f )(x) − f (x) = [f (x − t) − f (x)] DN (t) dt
π −π
1 π f (x − t) − f (x)
Z
= sin (N + 21 )t dt = b 1 (gx ),
π −π 2 sin ( 2t ) N+ 2

si
f (x − t) − f (x)
gx (t) = .
2 sin ( 2t )
La condición de regularidad (4.9) muestra que gx ∈ L1 (−π, π),
M |t|α
|gx (t)| ≤ . |t|α−1 , si |t| ≤ π
| sin ( 2t )|
y por el Lema de Riemann-Lebesgue,
lı́m b 1 (gx ) = 0.
N →+∞ N + 2

42
Desigualdad de Bessel
Teorema 12. Si f ∈ L2 (−π, π), entonces
∞ π
a20 X 2
Z
1
+ (ak + b2k ) ≤ |f (x)|2 dx. (4.10)
2 π −π
k=1

Demostración. La familia de funciones {1, cos kx, sin kx} es ortogonal en

L2 (−π, π), Z π Z π
2
cos kx dx = sin2 kx dx = π, si k ≥ 1
−π −π

y la desigualdad de Bessel asociada es (4.10).

Un resultado de convergencia uniforme

Sabemos que si una sucesión de funciones continuas converge uniforme-
mente, el lı́mite es una función continua y por tanto, si una serie de Fourier
converge uniformemente, su suma es una función continua.

Teorema 13. Sea f ∈ C 1 ([−π, π]) o C 1 a trozos y tal que f (π) = f (−π).
Entonces, su serie de Fourier converge uniformemente a f en [−π, π].

Demostración. Una tal f verifica las condiciones del teorema 11 con α = 1

y M = kf 0 kL∞ (−π,π) , y su serie de Fourier converge puntualmente a f en
[−π, π]. Al mismo tiempo,

|ak | + |bk | ≤ k −2 + k 2 a2k + b2k , si k ≥ 1.

(4.11)

La relación entre los coeficientes de Fourier de f y los de su derivada, la

desigualdad de Bessel y (4.11) implican que
+∞ +∞ Z π
X X
0 2 1
0 2
2
a2k b2k f 02 (x) dx


k + = ak (f ) + bk (f ) ≤
π −π
k=1 k=1

y la convergencia uniforme de la serie de Fourier se sigue del criterio de

Weierstrass y de la desigualdad

|ax cos kx + bk sin kx| ≤ |ak | + |bk |, si k ≥ 1.

Una vez conocido este resultado, podemos demostrar la identidad de

Parseval para {1, cos kx, sin kx} con un argumento de densidad y deducir
que dicha familia forma un sistema completo en L2 (−π, π).

Teorema 14. La serie de Fourier de una función f ∈ L2 (−π, π) converge

a f en L2 (−π, π) y la desigualdad de Bessel es una identidad.

43
Demostración. La desigualdad de Bessel es equivalente a que

N
! 12
√ a20 X
kSN (f )kL2 (−π,π) = π + (a2k + b2k ) ≤ kf kL2 (−π,π) , (4.12)
2
k=1

si N ≥ 1 y f ∈ L2 (−π, π), que junto con la desigualdad triangular implica

que

kf − SN (f )kL2 (−π,π) ≤ kf − ϕkL2 (−π,π) + kϕ − SN (ϕ)kL2 (−π,π)

 
1
+kSN (ϕ−f )kL2 (−π,π) ≤ 1 + √ kf −ϕkL2 (−π,π) +kϕ−SN (ϕ)kL2 (−π,π) .
π

En la desigualdad anterior, podemos elegir ϕ ∈ C0∞ (−π, π) que verifique

kf − ϕkL2 (−π,π) ≤ 

y como SN (ϕ) converge uniformemente a ϕ en [−π, π],

kϕ − SN (ϕ)kL2 (−π,π)

tiende a cero, si N → +∞ y concluimos que también SN f converge a f en

L2 (−π, π). Éste resultado y la identidad
N
!
a 2 X
kf − SN (f )k2L2 (−π,π) = kf k2L2 (−π,π) − π 0
+ a2k + b2k ,
2
k=1

muestran que la desigualdad de Bessel (4.12) es de hecho una identidad,

llamada identidad de Plancherel o Parseval,
+∞ π
a20 X 2
Z
1
+ (ak + b2k ) = f 2 (x) dx, si f ∈ L2 (−π, π),
2 π −π
k=1

y por tanto, la familia de funciones {1, cos kx, sin kx} es completa en L2 (−π, π).

Funciones pares e impares

• Series de Fourier de cosenos/senos. Dada una función en [0, π], se
llama serie de Fourier de cosenos de f a la serie de Fourier de su extensión
par al intervalo [−π, π]. La serie de Fourier de senos de f es la serie de
Fourier de su extensión impar a [−π, π].
• Las funciones {sin kx : k ≥ 1} y {cos kx : k ≥ 0} forman sistemas
completos en L2 (0, π) respectivamente.

44
Demostración. Sea fe la extensión par a [−π, π] de f en L2 (0, π). Los coefi-
cientes de Fourier de fe son respectivamente
1 π e 2 π
Z Z
ak =
e f (x) cos kx dx = f (x) cos kx dx, si k ≥ 0
π −π π 0
y Z π
ebk = 1 fe(x) sin kx dx = 0, si k ≥ 1.
π −π
Como
N
!
a0 X
kf − ak cos kx kL2 (−π,π)
e
e + e
2
k=1
tiende a cero y [0, π] ⊂ [−π, π], la serie de Fourier de cosenos
+∞
2 π
Z
a0 X
+ ak cos kx, ak = f (x) cos kx dx, si k ≥ 0
2 π 0
k=1

converge a f en L2 (0, π).

Si fe es la extensión impar de f , su serie de Fourier sólo tiene senos. En
este caso
2 π
Z
bk = ebk = f (x) sin kx dx, si k ≥ 1
π 0
y
XN N
X
kf − bk sin kxkL2 (0,π) ≤ kfe − ebk sin kxkL2 (−π,π) ,
k=1 k=1
que tiende a cero si N → +∞.

Cambio de periodo
• Si f es periódica de periodo 2l, su serie de Fourier en (−l, l) se escribe
∞    
a0 X kπx kπx
+ ak cos + bk sin ,
2 l l
k=1
con
Z l   Z l  
1 kπx 1 kπx
ak = f (x) cos dx, bk = f (x) sin dx
l −l l l −l l
y las funciones    
kπx kπx
{1, cos , sin }
l l
forman un sistema completo en L2 (−l, l).
• Las familias de funciones
{1, cos kπx , sin kπx } , {cos kπx kπx





l l l : k ≥ 0} , {sin l : k ≥ 1}
son sistemas completos en L2 (0, l), si l > 0.

45
Demostración. El cambio de variable, g(x) = f ( lx 2
π ), transforma L (−l, l) en
2
L (−π, π) y la información sobre la serie de Fourier clásica de g en [−π, π],
en información para la nueva serie de Fourier de f en (−l, l).

4.4. Vuelta al problema del calor en una barra fi-

nita
Los coeficientes bk en (4.4) son los de la serie de Fourier de senos de f
en (0, π):
2 π
Z
bk = f (x) sin kx dx, si k ≥ 1.
π 0
Esta sucesión de números está acotada, si f ∈ L1 (0, π) y la serie
∞
2t
X
bk e−k sin kx,
k=1

converge uniformemente junto con sus derivadas en R × [, +∞) para todo
 > 0, es una solución de la ecuación del calor y verifica la condición de
contorno en [0, π] × (0, +∞).
2
Demostración. Si α, β ∈ N, k α+2β e−k  ≤ C(α, β, ), para todo k ≥ 1,
 2
 2
|∂xα ∂tβ bk e−k t sin kx | . kf kL1 (0,π) k α+2β e−k t
 k
1
≤ C(α, β, , f )  , si x ∈ R y t ≥ 2
e
y la afirmación es consecuencia del criterio de Weierstrass y de la conver-
gencia de la serie de término general rk , 0 < r < 1.

Convergencia al dato inicial

Teorema 15. Si f ∈ L2 (0, π), entonces u( · , t) converge a f en L2 (0, π).
Demostración. Sabemos que {sin kx : k ≥ 1} es completa en L2 (0, π) y por
la identidad de Plancherel para esta familia
+∞
X 2 πX 2 2
ku(·, t) − f k2L2 (0,π) = k bk (e−k t − 1) sin kxk2L2 (0,π) = bk (1 − e−k t )2
2
k=1 k

y el lı́mite a la derecha es cero por el Teorema de Convergencia Dominada

o la siguiente acotación:
N
2 2
X X
b2k (1 − e−k t )2 ≤ b2k (1 − e−k t )2 +  ,
k k=1
P+∞ 2
si N se elige tal que, k=N +1 bk ≤ .

46
 Si f ∈ C([0, π]), f (0) = f (π) = 0 y fe es la extensión impar de f a
R respecto de 0 y π, fe ∈ C(R) ∩ L∞ (R). En tal caso, la solución de (4.1)
que obtenemos por el método de reflexiones está en C ∞ ([0, π] × (0, +∞)) ∩
C([0, π]×[0, +∞)). También se puede mostrar que la solución generada por el
método de separación de variables es, en tal caso, continua en [0, π]×[0, +∞),
y por el principio del máximo débil para un cilindro de base acotada, las
dos soluciones son iguales (Ver la sección titulada Extensión de la validez
de esas soluciones en el libro de Hans F. Weinberger pero no intentar su
lectura hasta haber leı́do el último tema en este curso).
Es fácil mostrar que la solución generada por el método de separación de
variables es continua en la clausura del cilindro [0, π] × R, si f ∈ C 1 ([0, π])
y f (0) = f (π) = 0.

Demostración. Los coeficietes de la serie de Fourier de senos de f verifican

Z π Z π
2 0 2 ak
bk = − f (x) (cos kx) dx = f 0 (x) cos kx dx = , si k ≥ 1,
πk 0 πk 0 k
donde ak es el coeficiente de Fourier de la serie de cosenos de f 0 en [0, π].
Esto implica que
1
|bk | ≤ 2 + a2k , si k ≥ 1, (4.13)
k
y (4.13), la identidad de Plancherel para la serie de cosenos de f 0 en [0, π] y
el criterio de Weiestrass, muestran que la serie (4.3) converge uniformemente
hacia una función continua en [0, π] × R.

Otras condiciones de contorno

• Si las condiciones de contorno son, ux (0, t) = ux (π, t) = 0 (de tipo
Neumann), llegamos a

X 00 + λX = 0, con X 0 (0) = X 0 (π) = 0,

que tiene soluciones no nulas para λ = 0 (las constantes) y λ = k 2 (múltiplos

de cos kx), si k ≥ 1. La solución se escribe como
∞
a0 X 2
u(x, t) = + ak e−k t cos kx,
2
k=1

y los ak se obtienen como los coeficientes de la serie de cosenos de f .

Métodos similares a los de la sección anterior muestran que u ∈ C ∞ ([0, π]×
(0, +∞)) si f es integrable, que u es continua en [0, π] × [0, +∞) si f ∈
C 1 ([0, π]) y que u ∈ C 1 ([0, π]×[0, +∞), si f ∈ C 2 ([0, π]) y f 0 (0) = f 0 (π) = 0.
• Si las condiciones de contorno son, u(0, t) = 0 y ux (π, t) = 0, estas
conducen al problema de Sturm-Liouville

X 00 + λX = 0, con X(0) = X 0 (π) = 0,

47
que tiene como valores propios λk = (k + 1/2)2 y como funciones propias a
Xk = sin(k + 1/2)x, si k ≥ 0. Ası́,
∞
2
X
u(x, t) = ck e−(k+1/2) t sin(k + 1/2)x
k=0

y para encontrar la solución por el método de la separación de variables

necesitamos saber que
∞
X
f (x) = ck sin(k + 1/2)x,
k=0

en algún sentido ¿Es {sin (k + 1/2)x : k ≥ 0} una familia completa en

L2 (0, π)?
Para resolver este problema procedemos como en el método de reflexiones
y de forma análoga a como probamos que las funciones asociadas a las
series de cosenos y senos son familias completas en L2 (0, π): si fe denota la
extensión de f a toda la recta real como función impar respecto a x = 0 y
par respecto a x = π, fe(2x) es periódica de periodo 2π y su serie de Fourier,
se reduce a una serie de senos impares. Como f (x) = fe( x2 ), resulta que f
admite un desarrollo con una serie de Fourier del tipo anterior con
2 π
Z
ck = f (x) sin(k + 1/2)x dx
π 0

que convege a f en L2 (0, π).

4.5. Ecuación no homogénea

El método de separación de variables requiere condiciones de contorno
homogéneas. Si no lo son, podemos escoger una función v que cumpla las
condiciones de contorno del problema, restarla de u y resolver un nuevo
problema para w = u − v. Esto puede afectar al segundo miembro de la
ecuación y a la condición inicial.

4.5.1. Segundo miembro no nulo

Consideramos ahora el problema no homogéneo

ut − uxx = F (x, t) ,
 0 < x < π,t > 0,
u(x, 0) = 0 , 0 < x < π,

u(0, t) = u(π, t) = 0 , t > 0.


que también se puede resolver por separación de variables: imitamos el méto-

do de variación de las constantes (o de Duhamel) y se escribe u como una

48
serie de funciones de variables separadas con nuevas funciones de la variable
temporal t que multiplican a las soluciones del problema de Sturm-Liouville
asociado al problema homogéneo; es decir, escribimos
∞
X
u(x, t) = Tk (t) sin kx,
k=1

y se requiere que u cumpla la ecuación y la condición inicial (el método

asegura las condiciones de contorno). Desarrollando F según las mismas
funciones propias, podremos escribir
∞
X
F (x, t) = Ck (t) sin kx
k=1

y las funciones Tk (t) se determinan como las soluciones de la ecuación dife-

rencial ordinaria (
Tk0 (t) + k 2 Tk (t) = Ck (t),
Tk (0) = 0,
que resolvemos por el método de variación de las constantes.

4.6. Separación de variables para la ecuación de

ondas
Sea el problema de ondas en una cuerda finita

2
utt − c uxx = 0 ,

 0 < x < π,t > 0,

u(0, t) = u(π, t) = 0 , t > 0 ,


u(x, 0) = f (x) , 0 < x < π,

u (x, 0) = g(x) ,
t 0 < x < π.

Buscamos soluciones de la forma X(x)T (t) que cumplan la ecuación y las

condiciones de contorno, X(0) = X(π) = 0. Ası́ se llega al mismo problema
de Sturm-Liouville que antes: los valores propios son k 2 , k ≥ 1 y las funciones
propias son sin kx. La ecuación diferencial ordinaria para Tk con cada uno
de éstos valores es
Tk00 (t) + c2 k 2 Tk (t) = 0,
que tiene como solución general, Tk (t) = A cos kct + B sin kct. Por tanto,
una posible solución se escribe como
∞
X
u(x, t) = (Ak cos kct + Bk sin kct) sin kx
k=1

y se determinan los coeficientes con las condiciones iniciales.

49
 Ası́ queda que
2 π 2 π
Z Z
Ak = f (x) sin kx dx y ckBk = g(x) sin kx dx.
π 0 π 0
En este caso es más difı́cil estudiar la convergencia de la serie a partir del
tamaño de los coeficientes. Se puede hacer indirectamente usando fórmulas
trigonométricas y observando que la serie que representa a u es la que co-
rresponde a la solución construida utilizando por el método de d’Alembert
o de las caracterı́sticas.
Las variantes con distinto tipo de condiciones de contorno o con segundo
miembro a la derecha son como en el caso de la ecuación del calor. En
particular, para resolver el problema no homogéneo



 utt − c2 uxx = F (x, t) , 0 < x < π , t > 0 ,

u(0, t) = u(π, t) = 0 , t > 0 ,
(4.14)


 u(x, 0) = 0 , 0 < x < π,

u (x, 0) = 0 ,
t 0 < x < π,

y en analogı́a con el metodo de variación de las constantes o de Duhamel,

intentamos escribir la solución de (4.14) como
∞
X
u(x, t) = Tk (t) sin kx. (4.15)
k=1

Sustituyendo en la ecuación y si F (x, t) = ∞

P
k=1 Ck (t) sin kx, se debe verifi-
car que (
Tk00 (t) + c2 k 2 Tk (t) = Ck (t),
Tk (0) = Tk0 (0) = 0,
para todo k ≥ 1. Como conocemos las soluciones de la ecuación homogénea,
por el método de variación de las constantes escribimos Tk (t) como

Ak (t) cos kct + Bk (t) sin kct,

donde las funciones Ak y Bk verificarán


0 0
Ak (t) cos kct + Bk (t) sin kct = 0,

−A0k (t) sin kct + Bk0 (t) cos kct = Ck (t)/kc, (4.16)

Ak (0) = Bk (0) = 0,


si k ≥ 1.
Las ecuaciones (4.16) implican formalmente que u es solución de (4.14).
El análisis de la convergencia de la serie (4.15) y de la serie de sus derivadas
parciales se puede hacer con el criterio de Weierstrass para la convergencia
uniforme de series de funciones.

50
4.7. Separación de variables en otros dominios
Si Ω es un abierto acotado de Rn con frontera regular y queremos resolver

4u − ∂t u = 0,
 en Ω × (0, +∞),
u = 0, en ∂Ω × [0, +∞),

u(x, 0) = f (x), en Ω,


podemos empezar por buscar soluciones de variables separadas, X(x)T(t),

que cumplan las condiciones de contorno. Esto da lugar a encontrar los λ ∈ R
tal que existe una solución no nula de
(
4X + λX = 0, en Ω,
X = 0, en ∂Ω,

y e−λt X es una tal solución. Es conocido que la solución al problema de

Sturm-Liouville es una familia numerable {ek } de funciones en C ∞ (Ω) que
verifican (
4ek + λ2k ek = 0, en Ω,
k = 1, 2, . . .
ek = 0, en ∂Ω,
donde
0 < λ1 < λ2 ≤ λ3 ≤ . . . λk ≤ . . . , lı́m λk = +∞,
k→+∞

y que {ek } es un sistema completo en L2 (Ω). Estas funciones se denominan

las autofunciones del operador de Laplace en Ω con condiciones de Dirichlet
nulas. Ası́ podemos construir la solución como
+∞ Z
−λ2k t
X
u(x, t) = ak e ek , si ak = f ek dx.
k=1 Ω

El problema

2
4u − ∂t u = 0,
 en Ω × (0, +∞),
u = 0, en ∂Ω × [0, +∞),

u(x, 0) = f (x), ut (x, 0) = g(x) en Ω,


se resuelve escribiendo
+∞ 
X 
bk
u(x, t) = ak cos λk t + λk sin λk t ek ,
k=1

con Z Z
ak = f ek dx y bk = gek dx.
Ω Ω

51
Capı́tulo 5

La ecuación del potencial en

el plano

5.1. Los problemas de Dirichlet y Neumann.

• En fı́sica aparecen campos incompresibles (conservativos) e irrotacio-
nales: campos magnéticos, flujos de temperatura, velocidades de fluidos, etc.
Si V = P i + Qj es un tal campo en Ω ⊂ R2 , lo anterior corresponde a que
su divergencia y rotacional

∇ · V = ∂x P + ∂y Q , ∇ × V = (∂x Q − ∂y P )k,

sean nulos en Ω. Un campo irrotacional es el gradiente de una función u :

Ω ⊂ R2 −→ R, V = ∇u, que es única excepto por contantes y que se
denomina el potencial del campo V . Como

∇ · V = ∇ · (∇u) = 4u = 0, en Ω,

el potencial es una función armónica.

En general se pueden hacer mediciones de V en la frontera de Ω y conocer
la componente normal de V (V · ν) o su componente tangential (V · t) en ∂Ω,
donde ν y t son los vectores normal y tangencial unitarios a ∂Ω. Esto equivale
a conocer la derivada normal de u en ∂Ω (V ·ν = ∇u·ν) o los valores de u en
∂Ω. Para entender el segundo caso, si u ∈ C 1 (Ω) y ∂Ω = {σ(t) : t ∈ [0, 1]},
donde σ : [0, 1] −→ R2 es una parametrización periódica del borde de Ω, se
verifica que
Z t
u(σ(t)) = u(σ(0)) + ∇u(σ(τ )) · σ̇(τ ) dτ
0
Z t Z
= u(σ(0)) + V (σ(τ )) · t(τ )|σ̇(τ )| dτ = u(σ(0)) + P dx + Q dy,
0 γt

52
donde γt es el arco orientado en ∂Ω que conecta σ(0) con σ(t). Por tanto,
los valores de V en Ω quedan determinados por el valor del gradiente de
cualquiera de las soluciones de los problemas:
( (
4u1 = 0, en Ω, 4u2 = 0, en Ω,
∂u1
o
∂ν = h, en ∂Ω, u2 = f, en ∂Ω,

donde h y f corresponden a las mediciones que se tengan de V en la frontera

de Ω.
• El primer problema se llama problema de Neumann y el segundo de
Dirichlet. Si Ω es acotado los dos problemas tienen solución única en C 2 (Ω)∩
C 1 (Ω).
∂u
Demostración. Si u ∈ C 2 (Ω) ∩ C 1 (Ω) es armónica en Ω y u o ∂ν es cero en
∂Ω,
∇ · (u∇u) = u4u + |∇u|2 = |∇u|2
y por el teorema de la divergencia,
Z Z Z
|∇u|2 dx = ∇ · (u∇u) dx = u ∂u
∂ν dσ = 0.
Ω Ω ∂Ω

• Si Ω es un dominio simplemente conexo y Ψ : B −→ Ω es la trans-

formación conforme de Riemann de la bola unidad B en Ω, las funciones
Ui = ui ◦ Ψ, i = 1, 2, son armónicas en B y verifican
( (
4U1 = 0, en B, 4U2 = 0, en B,
∂U1 0
,
∂ν = h ◦ Ψ|Ψ |, en ∂B, U2 = f ◦ Ψ, en ∂B.

• Esto justifica que intentemos resolver el problema de Dirchlet en el

cı́rculo (
4u = 0, en B,
(5.1)
u = f, en ∂B,

y encontraremos una solución u en C ∞ (B) ∩ C(B), si f ∈ C(∂B).

Utilizando coordenadas polares (r, θ), f se convierte en una función con-
tinua y periódica de periodo 2π,

ϕ(θ) = f (cos θ, sin θ) y u(r, θ) = u(r cos θ, r sin θ),

tendrá que estar en C ∞ ([0, 1) × [−π, π]) ∩ C([0, 1] × [−π, π]) y ser, junto con
sus derivadas parciales, periódica de periodo 2π, para cada 0 ≤ r ≤ 1.

53
 En la nuevas coordenadas, el problema de Dirichlet (5.1) se reduce a
encontrar u ∈ C ∞ ([0, 1) × [−π, π]) ∩ C([0, 1] × [−π, π]), periódica de periodo
2π y tal que
(
urr + 1r ur + r12 uθθ = 0, en [−π, π] × [0, 1),
(5.2)
u(1, θ) = f (θ), en [−π, π].

En estas coordenadas es posible resolver (5.2) por separación de varia-

bles. Las soluciones en variables separadas no nulas, R(r)Θ(θ), nos lleva a
las ecuaciónes
r2 R00 + rR0 Θ00
− = = −λ.
R Θ
Ası́ tendremos el problema de Sturm-Liouville

Θ00 + λΘ = 0, Θ y sus derivadas son periódicas de periodo 2π,

que equivale a encontrar los λ ∈ R tales que

(
Θ00 + λΘ = 0, en [−π, π],
(5.3)
Θ(−π) = Θ(π), Θ0 (−π) = Θ0 (π),

tiene solutiones no nulas. Esto ocurre si λ = 0 con Θ ≡ 1 o λ = k 2 y

Θk = A cos kθ + B sin kθ, si k ≥ 1.

Demostración. Si λ = 0, Θ = Aθ + B, con A, B en R y para ser periódica

A deber√ ser cero. Es
√ decir, Θ0 = 1 es solución con λ0 = 0. Si λ < 0,
Θ = Ae −λ θ + Be− −λ θ y que Θ y Θ0 sean periódicas nos lleva al sistema
de ecuaciones
( √ √
A sinh ( −λ π) − B sinh ( −λ π) = 0,
√ √
A sinh ( −λ π) + B sinh ( −λ π) = 0,
√ √
cuya única solución es la nula. Si λ > 0, Θ = A cos ( λ θ) + B sin ( λ θ) y
las condiciones de contorno implican que
( √ √ √ √
A cos ( λ π) + B sin ( λ π)B = A cos ( λ π) − B sin ( λ π)B,
√ √ √ √
−A sin ( λ π) + B cos ( λ π)B = A sin ( λ π) + B cos ( λ π)B,
√
que tiene solución no nula si sin ( λ π) = 0; es decir, para λk = k 2 , k ≥ 1 y
con Θk = Ak cos kθ + Bk sin kθ.

Resolvemos la EDO para R según los valores hallados de λ

 
00 d dR
r2 R + rR0 − λR = r r − λR = 0.
dr dr

54
 d d
La ecuación es de tipo Euler y por el cambio de variables r = es , r dr = ds ,
es equivalente a
d2 R
− λR = 0.
ds2
Si λ = 0, tenemos, R(r) = A log r + B y si λ = k 2 , R(r) = Ar−k + Brk .
Para que u sea continua en 0 debemos eliminar las soluciones de varia-
bles separadas que no lo son, las asociada a log r y r−k y con las restantes
escribimos formalmente
∞
a0 X k
u(r, θ) = + r (ak cos kθ + bk sin kθ). (5.4)
2
k=1

Si pedimos a u que cumpla la condición de contorno deducimos que ak

y bk deben ser los coeficientes de Fourier de f en L2 (−π, π).

Regularidad
Si f ∈ L1 (−π, π), los coeficientes ak y bk están acotados y la serie con-
verge absoluta y uniformemente en la bola B(0, 1 − ρ), para todo 0 < ρ < 1
y se puede derivar dentro de la suma tantas veces como se quiera: u está en
C ∞ (B) y 4u = 0 en el interior de B.

Demostración. Por Weierstrass y la acotación

|k l rk (ak cos kθ + bk sin kθ)| ≤ Cl kf kL1 (−π,π) (1 − ρ2 )k ,

si r ≤ 1 − ρ, k, l ≥ 0 y θ ∈ R.

Núcleo de Poisson
Sustituyendo ak y bk por la fórmula de los coeficientes:
∞
" #
1 π
Z
1 X k
u(r, θ) = f (φ) + r cos k(θ − φ) dφ
π −π 2
k=1
1 π
Z
= f (φ)Pr (θ − φ) dφ,
π −π

donde
∞
1 X k 1 − r2
Pr (θ) = + r cos kθ =
2 2(1 − 2r cos θ + r2 )
k=1

es el núcleo de Poisson.

55
Demostración. La fórmula anterior es consecuencia de las fórmulas
∞
X
[r2 + 1 − 2r cos θ] rk cos kθ
1
∞
X
= {[rk+2 + rk ] cos kθ − rk+1 [cos (k + 1)θ + cos (k − 1)θ]}
1
∞
X ∞
X ∞
X ∞
X
k+2 k k
= r cos kθ + r cos kθ − r cos kθ − rk+2 cos kθ
1 1 2 0
= r cos θ − r2 ,
∞
X r cos θ − r2
rk cos kθ = ,
r2 + 1 − 2r cos θ
1

y de sumar 1/2 a la suma infinita anterior.

Si x = (x1 , x2 ), x1 = r cos θ, x2 = r sin θ y si |y| = 1, y = (cos φ, sin φ),

la longitud de arco en ∂B es dσ = dφ y

1 − r2 1 − |x|2
= ,
1 − 2r cos (θ − φ) + r2 |x − y|2

de donde
1 − |x|2
Z
1
u(x) = f (y) dσ(y). (5.5)
2π ∂B |x − y|2
Si en el razonamiento anterior hacemos f ≡ 1, sale u ≡ 1; es decir,

1 − |x|2
Z
1
dσ(y) = 1, si |x| < 1.
2π ∂B |x − y|2

Continuidad hasta el borde

Hemos visto que u ∈ C ∞ (B), si f ∈ L1 (∂B). Si además f ∈ C(∂B),
entonces u ∈ C(B) y u = f en ∂B.

Demostración. Si y0 ∈ ∂B y  > 0, existe δ > 0 tal que |f (y) − f (y0 )| < ,

si |y − y0 | < δ, |y| = 1. Entonces,

1 − |x|2
Z
1
u(x) − u(y0 ) = (f (y) − f (y0 )) dσ(y)
2π ∂B |x − y|2
1 − |x|2
Z
1
= (f (y) − f (y0 )) dσ(y)
2π ∂B∩Bδ (y0 ) |x − y|2
1 − |x|2
Z
1
+ (f (y) − f (y0 )) dσ(y).
2π ∂B\Bδ (y0 ) |x − y|2

56
La primera integral es menor que . La segunda, está acotada por

1 − |x|2
2kf kL∞ (∂B) , si |x| ≤ δ/2
(δ/2)2

y tiende a cero, si x ∈ B se acerca a y0 .

Cı́rculo de radio R
Si x0 ∈ R2 , R > 0 y ϕ ∈ C(∂BR (x0 )). Una función u en C 2 (BR (x0 )) ∩
C(B R (x0 )) es solución del problema de Dirichlet
(
4u = 0, en BR (x0 ),
(5.6)
u = ϕ, en ∂BR (x0 ),

si v(ξ) = u(x0 + Rξ) es solución en C 2 (B) ∩ C(B) de (5.1) con f (ξ) =

ϕ(x0 + Rξ). Para el cı́rculo unidad sabemos como construir una solución de
(5.5) y deshaciendo el cambio de variables, obtenemos una solución de (5.6).
De (5.5),

1 − |ξ|2
Z
1
u(x0 + Rξ) = ϕ(x0 + Ry) dσ(y), si |ξ| < 1.
2π ∂B |ξ − y|2

y escribiendo, x = x0 + Rξ,

R2 − |x − x0 |2
Z
1
u(x) = ϕ(x0 + Ry) dσ(y), si |x − x0 | < R.
2π ∂B |x − (x0 + Ry)|2

Finalmente, parametrizando ∂B como y = eiθ , −π ≤ θ ≤ π, x0 + Reiθ

es una parametrización de ∂BR (x0 ) y

R2 − |x − x0 |2
Z
1
u(x) = ϕ(y) dσ(y), si |x − x0 | < R, (5.7)
2πR ∂BR (x0 ) |x − y|2

es una solución de (5.6) que hemos construido utilizando el método se sepa-

ración y cambios de variables.

Algunas propiedades de esta solución.

La solución (5.7) verifica las siguientes propiedades:

•Principio del máximo y del mı́nimo.

mı́n f ≤ u(x) ≤ máx f, si x ∈ BR (x0 ).

∂BR (x0 ) ∂BR (x0 )

57
Demostración. Observar que

R2 − |x − x0 |2
Z
1
dσ(y) = 1, si |x − x0 | < R.
2πR ∂BR (x0 ) |x − y|2

•Propiedad del valor medio.

Z
1
u(x0 ) = f dσ.
2πR ∂BR (x0 )

Demostración. Hacer x = x0 en (5.7).

•Desigualdad de Harnack. Si u ≥ 0 en BR (x0 ), entonces

R − |x − x0 | R + |x − x0 |
u(x0 ) ≤ u(x) ≤ u(x0 ),
R + |x − x0 | R − |x − x0 |

para todo x ∈ BR (x0 ). Además,

sup u≤9 ı́nf u.

BR/2 (x0 ) BR/2 (x0 )

Demostración. Observar que

R − |x − x0 | R2 − |x − x0 |2 R + |x − x0 |
≤ 2
≤ ,
R + |x − x0 | |x − y| R − |x − x0 |

si x ∈ BR (x0 ) e y ∈ ∂BR (x0 ). Para la segunda desigualdad utilizar (5.7) y

que
R + |x − x0 | R − |x − x0 | 1
≤ 3, ≥ ,
R − |x − x0 | R + |x − x0 | 3
si x ∈ B R (x0 ), y ∈ ∂BR (x0 ).
2

5.2. Principio del máximo

Teorema 16 (Principio del máximo débil). Sea Ω un abierto acotado del
plano y u ∈ C 2 (Ω) ∩ C(Ω). Si 4u ≥ 0 en Ω, entonces

máx u = máx u.
Ω ∂Ω

Demostración. Si 4u > 0 en Ω y u alcanza su máximo en Ω en algún x en

Ω, la matriz Hessiana de u en x es semidefinida negativa y su traza 4u, es
menor o igual que 0. Si 4u ≥ 0, aplicar lo anterior a u = u + |x|2 y hacer
tender  a cero.

58
 El principio del máximo implica que el problema de Dirichlet para el
operador de Laplace tiene a lo más una solución u ∈ C 2 (Ω) ∩ C(Ω).
Teorema 17 (Propiedad del valor medio). Sea Ω un abierto del plano,
u ∈ C 2 (Ω) armónica en Ω, x0 ∈ Ω y R > 0 tal que B R (x0 ) ⊂ Ω. Entonces,
Z
1
u(x0 ) = u(y) dσ.
2πR ∂BR (x0 )

Demostración. Por la unicidad de solución para problema de Dirichlet, u

coincide en BR (x0 ) con la solución v construida en (5.7) para
(
4u = 0, en BR (x0 ),
v = u, en ∂BR (x0 ),

y esta verifica la propiedad del valor medio en BR (x0 ).

El recı́proco también es cierto.

Teorema 18. Si u en C(Ω) cumple la propiedad del valor medio, entonces
u es C ∞ en Ω y armónica en Ω.
Demostración. La función
− 1
(
N e 1−|x|2 , si |x| < 1,
ρ(x) =
0, si |x| ≥ 1,

es C ∞ en R2 , tiene integral igual a 1 si N se elige de forma adecuada. Si

definimos Z
u (x) = −2 ρ( x−y
 )u(y) dy,
Ω
u esC∞ en R2 . Por la propiedad del valor medio e integración en coorde-
nadas polares,
Z Z Z
−2 1 r
u (x) =  ρ( x−y
 )u(y) dy = 2
ρ(  )u(y) dσdr
B (x) 0 ∂Br (x) 
Z  Z
2πr r
= u(x) 2
ρ(  ) dr = u(x) ρ(y) dy = u(x),
0 R2

si d(x, ∂Ω) > . Por tanto, u es C ∞ en Ω. Para ver que u es armónica en Ω,

desarrollamos por Taylor hasta orden 2 alrededor de x en Ω,

u(y) = u(x) + ∇u(x) · (y − x) + 12 D2 u(x)(y − x) · (y − x) + O(|y − x|3 ),

y como
Z Z
1 1 r2 δij
2πr (yi − xi ) dσ = 0 , 2πr (yi − xi )(yj − xj ) dσ = 2 ,
∂Br (x) ∂Br (x)

59
si i, j = 1, . . . , n, tendremos que
Z
r2
u(x) = 2πr 1
u(y) dσ = u(x) + 4 4u(x) + O(r3 );
∂Br (x)

es decir,
0 = 14 4u(x) + O(r),
y hacemos r tender hacia 0.

Teorema 19 (Principio del máximo fuerte). Sea Ω es un abierto conexo en

el plano y u en C 2 (Ω) ∩ C(Ω) una función armónica en Ω. Si u alcanza en
Ω un máximo o mı́nimo global, entonces u es constante en Ω.

Demostración. Sea M el valor del máximo de u en Ω. Entonces,

H = {x ∈ Ω : u(x) = M }

es no vacı́o, cerrado y por la propiedad del valor medio también es abierto

en Ω: si y ∈ H
Z Z
1
M = u(y) = u(x) dσ , (M − u(x)) dσ = 0,
2πR ∂BR (y) ∂BR (y)

si R < d(y, ∂Ω). Pero M − u ≥ 0 en Ω y concluimos que u ≡ M en

B(y, d(y, ∂Ω)). Por tanto, H = Ω y u es constante en Ω.

Teorema 20 (Un teorema de tipo Liouville). Si u es armónica en R2 y

acotada inferiormente (o superiormente), entonces u es constante en R2 .

Demostración. Si u ≥ 0, fijar R > 0, escribir la desigualdad de Harnack

para u en BR (0) y hacer tender R a infinito. Esto implica que u(x) = u(0),
para todo x. Si u ≥ −M , trabajar con u + M y si u ≤ M , trabajar con
M − u.

5.3. El problema no homogéneo en el cı́rculo

El problema no homogéneo
(
4u = F, en B,
u = 0, en ∂B,

es equivalente en coordenadas polares a encontrar u(r, θ), periódica de pe-

riodo 2π junto con sus derivadas y tal que
(
urr + 1r ur + r12 uθθ = F (r, θ), en [−π, π] × [0, 1),
u(1, θ) = 0, en [−π, π],

60
si F (r, θ) = F (r cos θ, r sin θ). Escribimos la posible solución como
∞
a0 (r) X
u(r, θ) = + ak (r) cos kθ + bk (r) sin kθ.
2
k=1

Esta función será solución del problema no homogéneo si

 2


 a00k + 1r a0k + kr2 ak = Fk (r), en [0, 1),
b00 + 1 b0 + k2 b = G (r),

en [0, 1),
k r k r2 k k


 ak (1) = bk (1) = 0,

a , b son acotadas en [0, 1],
k k

y
∞
F0 (r) X
F (r, θ) = + Fk (r) cos kθ + Gk (r) sin kθ,
2
k=1

y las ecuaciones se resuelven por el método de variación de las constantes.

5.4. El problema de Dirichlet en otros dominios

circulares
Cuando se trabaja en dominios como un anillo o el exterior de un cı́rculo,
donde todavı́a es posible separar las variables al trabajar en coordenadas
polares, se plantean los problemas de encontrar u de clase C 2 en BR2 \ B R1 ,
0 < R1 < R2 < ∞, continua en la clausura de BR2 \ BR1 y tal que

4u = 0, en BR2 \ B R1 ,

u = f2 , en ∂BR2 ,

u = f1 , en ∂BR1 ,


o de encontrar u de clase C 2 en Rn \ B R , R > 0, continua y acotada en

Rn \ BR y tal que (
4u = 0, en Rn \ B R ,
u = f, en ∂BR .
En coordenadas polares los problemas son respectivamente equivalentes
a encontrar u de clase C ∞ en (R1 , R2 ) × [−π, π], periódica de periodo 2π
junto con sus derivadas, continua en [R1 , R2 ] × [−π, π] y tal que

1 1
urr + r ur + r2 uθθ = 0,
 en (R1 , R2 ) × [−π, π],
u(R2 , θ) = ϕ2 (θ) = f2 (R2 cos θ, R2 sin θ), en [−π, π],

u(R1 , θ) = ϕ1 (θ) = f1 (R1 cos θ, R1 sin θ), en [−π, π],


61
o de encontrar u de clase C ∞ en (R, +∞) × [−π, π], periódica de periodo 2π
junto con sus derivadas, continua y acotada en [R, +∞) × [−π, π] y tal que
(
urr + 1r ur + r12 uθθ = 0, en (R, +∞) × [−π, π],
u(R, θ) = ϕ(θ) = f (R cos θ, R sin θ), en [−π, π].

En el primer caso, lo razonable es escribir u como

+∞  k
1 X r
u(r, θ) = (a0 + b0 log r) + (ak cos kθ + bk sin kθ)
2 R2
k=1
+∞ 
R1 k
X 
+ (dk cos kθ + ck sin kθ) , (5.8)
r
k=1

donde los coeficientes ak , bk , dk y ck se determinan por comparación de

series de Fourier e imponiendo que u(R1 , θ) = ϕ1 (θ), u2 (R2 , θ) = ϕ2 (z). En
el segundo, descartamos las soluciones de variables separadas que no son
acotadas en (R, +∞) × [−π, π] y escribimos u como
∞  k
a0 X R
u(r, θ) = + (ak cos kθ + bk sin kθ), (5.9)
2 r
k=1

donde ak y bk son los coeficientes de Fourier de ϕ.

Con un análisis similar al que hemos hecho con la solución generada
por separación de variables para el problema de Dirichlet en el interior de
una bola, se puede mostrar que las soluciones (5.8) y (5.9) son continuas y
acotadas en la clausura de sus dominios cuando su dato frontera es continuo
en la frontera del dominio.
El problema de Dirichlet en el exterior de una bola BR tiene a lo más una
solución en la clase de funciones C 2 (R2 \ BR ) ∩ C(R2 \ BR ) ∩ L∞ (R2 \ BR ):
basta con considerar
u = ± u −  log Rr ,
aplicar el principio del máximo débil a u en BR0 \ BR , donde R0 > R se
elige grande para que u ≤ 0 en ∂BR0 y se hace tender  a cero para concluir
que
± u(x) ≤ ∓ máx u, si x ∈ R2 \ BR .
∂BR

La condición de estar en L∞ (R2 \ BR ) es necesaria para la unicidad,

pues u = log R − log r es armónica en R2 \ B R , u = 0 en ∂BR y u no es
idénticamente nula.
También se puede usar separación de variables para resolver el problema
de Dirichlet en la región limitada por las rectas θ = α, θ = β y la circun-
ferencia r = R y con datos de Dirichlet nulos en los lados rectilı́neos del

62
cono: 
1 1
urr + r ur + r2 uθθ = 0,
 0 < r < R, α < θ < β,
u(r, α) = u(r, β) = 0, 0 ≤ r < 1, (5.10)

u(1, θ) = f (θ) α ≤ θ ≤ β.


Las soluciones en variables separadas no nulas, R(r)Θ(θ), nos lleva a las

ecuaciónes
r2 R00 + rR0 Θ00
− = = −λ
R Θ
y a buscar los λ ∈ R tales que
(
Θ00 + λΘ = 0, en (α, β),
Θ(α) = Θ(β) = 0,
tiene solución no nula. La solución son números positivos
λ21 < λ22 < . . . < λ2k . . .
con soluciones asociadas Θ1 , Θ2 , . . . , Θk . . . que forman un sistema completo
en L2 (α, β). Las soluciones de variables separadas acotadas son rλk Θk y la
solución de (5.10) se escribe como
+∞
X
u(r, θ) = ak (r/R)λk Θk .
k=1

La separación de variables también se puede utilizar para resolver el

problema de Dirichlet en un anillo limitado por un cono (θ = α y θ = β)
y la circunferencias r = R1 y r = R2 , 0 ≤ R1 < R2 ≤ ∞, con datos de
Dirichlet nulos en los lados rectilı́neos del cono y dos funciones dadas sobre
los arcos de las circunferencias dentro del cono. En este caso las soluciones
de variables separadas son r±λk Θk y la solución se puede escribir como
+∞
X
u(r, θ) = ak (R1 /r)λk Θk + bk (r/R2 )λk Θk .
k=1

5.5. El Problema de Neumann

Queremos encontrar u ∈ C ∞ (B) ∩ C 1 (B) tal que
(
4u = 0, en B,
∂u
∂ν = f, en ∂B,
si f ∈ C(∂B). Todo es como en el problema de Dirichlet hasta que llegamos a
(5.4). En ∂B el vector normal exterior es (x, y)/r y ∂u ∂u
∂ν = ∂r y por separación
de variables, escribimos la posible solución como
∞
a0 X k
u(r, θ) = + r (ak cos kθ + bk sin kθ).
2
k=1

63
La condición de contorno de tipo Neumann obliga a que
∞
X
k(ak sin kθ + bk cos kθ) = f (θ),
k=1

que sólo tiene solución si Z

f dσ = 0.
∂B
Si esto se cumple, kak y kbk deben ser los coeficientes de Fourier de f y la
solución queda determinada excepto por la constante a0 .
La condición de integral nula para f es necesaria para la existencia de
solución. En efecto, por el teorema de la divergencia
Z Z
∂u
4u dx = dσ.
B ∂B ∂ν

5.6. La ecuación del potencial en un rectángulo

El problema de Dirichlet en el rectángulo Ω = [0, a] × [0, b] para la
ecuación de Laplace
(
uxx + uyy = 0, en Ω,
(5.11)
u = ϕ, en ∂Ω,

ϕ ∈ C(∂Ω), se reduce a resolver el siguiente problema:


uxx + uyy = 0,
 0 < x < a, 0 < y < b,
u(0, y) = u(a, y) = 0, 0 < y < b, (5.12)

u(x, 0) = f1 (x), u(x, b) = f2 (x), 0 < x < a.


Para resolverlo buscamos soluciones no nulas de variables separadas,

X(x)Y (y), que verifiquen X(0) = X(a) = 0. El problema de Sturm-Liouville
es como los del capı́tulo anterior: los valores propios son λ = (kπ/a)2 , k ≥ 1
y las funciones propias asociadas, múltiplos de sin (kπx/a). La EDO para Y
con estos valores propios tiene por solución general
   
kπy kπy
A cosh + B sinh
a a

y la solución del problema se puede escribir como

∞       
X kπy kπy kπx
u(x, y) = Ak cosh + Bk sinh sin .
a a a
k=1

64
 Sustituyendo y por 0 y b e igualando a f1 (x) y f2 (x) respectivamente,
se determinan los coeficientes por comparación de series Fourier. Para de-
terminar Ak y Bk , hay que resolver un sistema lineal de dos ecuaciones con
dos incógnitas, para cada k ≥ 1.
El problema (5.11) se resuelve descomponiendo ϕ en suma de dos fun-
ciones que “viven” en los lados opuestos del rectángulo Ω, resolviendo el
análogo del problema (5.12), pero cambiando los lados horizontales del cua-
drado por los verticales y escribiendo la solución de (5.11) como la suma de
la solución del último problema descrito y la de (5.12).
Para el problema no homogéneo

uxx + uyy = F,
 0 < x < a, 0 < y < b,
u(0, y) = u(a, y) = 0, 0 < y < b,

u(x, 0) = 0, u(x, b) = 0, 0 < x < a,


escribimos la posible solución como

+∞  
X kπx
u(x, y) = Yk (y) sin ,
a
k=1

donde las funciones Yk son las soluciones de

(
2
Yk00 − kπ
a Yk = Fk (y), 0 < y < b,
Yk (0) = Yk (b) = 0
y
+∞  
X kπx
F (x, y) = Fk (y) sin , en [0, a] × [0, b]
a
k=0

es la serie de Fourier de F ( · , y) en [0, a].

El problema en un cilindro infinito

Se considera el problema siguiente: encontrar u acotada tal que

uxx + uyy = 0,
 0 < x < a, y > 0,
u(0, y) = u(a, y) = 0, y > 0,

u(x, 0) = f (x), 0 < x < a.


La primera parte es igual que en el caso anterior: para que la solución

sea acotada en el cilindro [0, a] × [0, +∞), excluimos ekπy/a y queda
∞
X
Bk e−kπy/a sin kπx


u(x, y) = a . (5.13)
k=1

65
 Los coeficientes Bk se determinan a partir de la serie de Fourier de f .
• La solución es única en la clase de funciones

C 2 ((0, a) × (0, +∞)) ∩ C([0, a] × [0, +∞)) ∩ L∞ ([0, a] × [0, +∞)).

• Se puede comprobar que (5.13) verifica estas condiciones, si f ∈ C([0, a])

y f (0) = f (a).

5.7. Ecuación del potencial en un semiplano

Consideramos el problema
(
uxx + uyy = 0, x ∈ R, y > 0,
(5.14)
u(x, 0) = f (x), x ∈ R.

• La ecuación homogénea tiene soluciones no nulas que se anulan en la

frontera, u(x, y) = y, por lo que no hay unicidad sin hipótesis adicionales.
El siguiente principio del máximo para el semiplano superior, muestra
que (5.14) tiene a lo más una solución u ∈ C 2 (R2+ ) ∩ C(R2+ ) que es acotada.

Teorema 21. Si u ∈ C 2 (R2+ ) ∩ C(R2+ ) verifica u ≤ M para algún M ≥ 0 y

4u ≥ 0 en R2+ . Entonces,
sup u ≤ sup f.
R2+ R

Demostración. Sean ΩT = BT (x0 ) × (0, T ), si T > 0, x0 ∈ R y

y
uT = u − 1
3 (x − x0 )2 + 3 (y − 2T ).
T 2
T2
Se verifica que
4uT ≥ 0, en ΩT
y por el principio del máximo débil para abiertos acotados

sup uT ≤ sup uT .
ΩT ∂ΩT
√
En la frontera lateral de
√ Ω T , uT ≤ M − T ≤ supR f , si T es grande.
En la superior, uT ≤ M − T ≤ supR f , si T es grande. Entonces,

y
uT (x0 , y) = u(x0 , y) + 3 (y − 2T ) ≤ sup f, si 0 ≤ y ≤ T
T2 R

y hacemos T tender a +∞.

66
 • Para resolver (5.14) buscamos soluciones autosemejantes: si u(x, y) es
armónica también lo es u(λx, λy) y si es homogénea, u(λx, λy) = λα u(x, y),
para todo λ > 0,
   
α x α x
u(x, y) = y u ,1 = y Φ , si y > 0.
y y

Si tenemos una solución, sus trasladadas en la variable x también son

soluciones
y α Φ(y −1 (x − z))
y las usamos para escribir u como “suma” de soluciones
Z ∞   Z ∞
α x−z α+1
u(x, y) = y Φ c(z) dz = y Φ(z)c(x − zy) dz, (5.15)
−∞ y −∞

donde c(z) debe ser elegido para que se cumpla la condición de contorno.
• Para que el lı́mite cuando y tiende a cero en (5.15) no sea 0 o ∞
necesitamos que α sea −1 y que Φ sea integrable en R. Para que el lı́mite
sea f elegimos c = f , si Φ tiene integral igual a 1 en R. Para que el laplaciano
(5.15) sea nulo, Φ debe verificar la EDO

(1 + s2 )Φ00 (s) − 2(α − 1)sΦ0 (s) + α(α − 1)Φ(s) = 0,

que cuando α = −1 es
00

(1 + s2 )Φ = 0.
Las soluciones integrables de esta ecuación son Φ(s) = N (1 + s2 )−1 y la
que tiene integral 1 corresponde a N = 1/π.
• Ası́ obtenemos que un candidato a solución de (5.14) es
Z
1 y
u(x, y) = f (y) dz. (5.16)
π R (x − z)2 + y 2

Teorema 22. La función u definida por (5.16) está en C ∞ (R2+ ) ∩ C(R2+ ) ∩

L∞ (R2+ ), si f ∈ C(R) ∩ L∞ (R) y es solución de (5.14).

Demostración. El nucleo
1 y
P (x, y) =
π x + y2
2

verifica
C(α, β, , R)
|∂xα ∂yβ P (x, y)| ≤ , si x ∈ R,  ≤ y ≤ R
1 + x2
y
C(α, β, , R)
|∂xα ∂yβ P (x − z, y)| ≤ , si |x| ≤ R,  ≤ y ≤ R,
1 + z2

67
que junto con el teorema de convergencia dominada muestra que u es regular
para y > 0.
Para acotar u, recordamos que
Z Z
1 y 1 1
2 2
dz = dz = 1, si y > 0,
π R (x − z) + y π R 1 + z2
y  Z 
1 y
|u(x, y)| ≤ kf kL∞ (R) dz ≤ kf kL∞ (R) .
π R (x − z)2 + y 2
Para la continuidad de u en (x0 , 0), dado  > 0 existe un δ > 0 tal que,
|f (x) − f (x0 )| ≤ , si |x − x0 | ≤ δ. Entonces,
Z
|u(x, y) − f (x0 )| ≤ P (x − z, y)|f (z) − f (x0 )| dz
Bδ (x0 )
Z
+ P (x − z, y)|f (z) − f (x0 )| dz
R\Bδ (x0 )
Z
≤  + 2kf kL∞ (R) P (x − z, y) dz.
R\Bδ (x0 )

Pero, P (x − z, y) ≤ N P (x0 − z, y), si |x − x0 | ≤ 2δ , |z − x0 | ≥ δ y

Z Z Z
N 1
P (x − z, y) dz ≤ N P (x0 − z, y) dz = π dz,
R\Bδ (x0 ) R\Bδ (x0 ) |z|≥ yδ 1 + z2

que tiende a cero, si y desciende hacia 0.

• Si f es par o impar respecto a l, u( · , y) también es par o impar respecto

a l para todo y > 0 y los problemas

uxx + uyy = 0, en (0, +∞) × (0, +∞),

u(0, y) = 0, en [0, +∞)

u(x, 0) = f (x), en [0, +∞),



uxx + uyy = 0,
 en (0, a) × (0, +∞),
u(0, y) = u(a, y) = 0, en [0, +∞)

u(x, 0) = f (x), en [0, a],


se pueden resolver por el método de reflexiones.

• La solución de estos dos problemas es única en la clase de funciones
que son C 2 en el interior del dominio, continuas y acotadas en la clausura.
• Cuando f verifica la condición de compatibilidad, la solución generada
por el método de reflexiones está en la clase de unicidad.

68