FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing.

Industrial 1

MATEMATICA D1
MODULO 1- UNIDAD 3: ECUACIONES DIFERENCIALES EN
DERIVADAS PARCIALES DE SEGUNDO ORDEN.
MÉTODO DE SEPARACIÓN DE VARIABLES

Introducción: Una ecuación diferencial en derivadas parciales lineal en dos

variables es de la forma:

Auxx + Buxy + Cuyy + Dux + Euy + F u = G

donde A, B, C, .., G son funciones de las variables x e y.

Generalmente en la teorı́a y aplicación de las ecuaciones diferenciales en derivadas
parciales relativas a una función u(x, y) se exige que u satisfaga no solamente la
ecuación diferencial en algún dominio de las variables independientes, sino también
alguna condición en los contornos del aquel dominio.
Las ecuaciones que representan tales condiciones de contorno pueden contener valores
de las derivadas de u, como también de la propia u, en puntos del contorno. Además
se exigen ciertas condiciones de continuidad tanto de u como de sus derivadas dentro
del dominio y en el contorno de éste.
La unión de todas las condiciones constituyen un “problema de contorno para la
función u”

1. Método de separación de variables

El objetivo es encontrar un conjunto de soluciones particulares de un problema de

contorno con condiciones homogéneas utilizando ecuaciones diferenciales ordinarias.
Se determinarán soluciones del tipo y(x, t) = X(x)T (t) que satisfagan las condiciones
de contorno.
Ejemplo: Trataremos de encontrar una fórmula para los desplazamientos trans-
versales y(x, t) de una cuerda que se encuentra estirada entre los puntos (0, 0) y (c, 0)
si la cuerda se encuentra desplazada inicialmente a la posición f (x) y luego se la deja
libre. Suponemos que a lo largo de la cuerda no actúan fuerzas exteriores.
El problema consiste en encontrar una función y(x, t) que satisfaga la ecuación de
ondas:

ytt (x, t) = a2 yxx (x, t) 0<x<c t>0 (1)

que satisfaga las condiciones de contorno:

y(0, t) = 0 y(c, t) = 0 yt (x, 0) = 0 (2)

y(x, 0) = f (x) 0<x<c (3)

donde la función f (x) es continua en el intervalo 0 ≤ x ≤ c y f (0) = f (c) = 0
FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing. Industrial 2

Buscaremos y(x, t) = X(x)T (t), luego yxx (x, t) = X 00 (x)T (t) y ytt (x, t) = X(x)T 00 (t)
Si y satisface la ecuación (1). entonces

X(x)T 00 (t) = a2 X 00 (x)T (t)

Dividiendo por a2 XT para separar las variables:

X 00 (x) T 00 (t)
= 2 (4)
X(x) a T (t)

Puesto que el primer miembro depende solo de x no puede variar con t, sin embargo
es igual a una función que depende solamente de t y que no puede variar con x. Por lo
tanto ambos miembros deben tener un valor constante común que podemos designar
por λ, luego:
X 00 (x) = λX(x) T 00 (t) = λa2 T (t) (5)
Puesto que y debe satisfacer las condiciones (2), entonces:

y(0, t) = 0 = X(0)T (t) y(c, t) = 0 = X(c)T (t) yt (x, 0) = 0 = X(x)T 0 (0)

En cada una de las ecuaciones anteriores no puede ser nulas ni T (t) ni X(x) puesto
que se obtiene una solución trivial. Resulta, entonces que y satisface las condiciones
(1) y (2) cuando X y T satisfacen los dos problemas:

X 00 (x) = λX(x) X(0) = 0 X(c) = 0 (6)

T 00 (t) = λa2 T (t) T 0 (0) = 0 (7)

donde el parámetro λ toma el mismo valor en ambos problemas. Para hallar soluciones
no triviales para los dos problemas, se puede observar el el problema (7) tiene una
sola condición de contorno y, por lo tanto puede tener muchas soluciones para cada
valor de λ, en cambio, el problema (6) tiene dos condiciones de contorno y puede ser
que para determinados valores de λ tenga soluciones no triviales.
Si λ = 0 la ecuación diferencial del problema (6) queda: X 00 (x) = 0; cuya so-
lución es X(x) = Ax + B y como X(0) = 0 = X(c), tiene que ser B = 0 y A = 0 y
X(x) = 0

Si λ = −α2 < 0, la ecuación diferencial del problema (6) queda:

X 00 (x) + α2 X(x) = 0; cuya solución es X(x) = C1 sen(αx) + C2 cos(αx)
Entonces X(0) = 0 si C2 = 0 y X(c) = 0 si C1 sen(αc) = 0, es decir:
nπ
α= n = 1, 2, ...
c
nπx
luego X(x) = C1 sen n = 1, 2, ...
c
Si λ = α2 > 0, la ecuación diferencial del problema (6) queda:
X 00 (x) − α2 X(x) = 0; cuya solución es X(x) = C1 eαx + C2 e−αx que solo verifica las
FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing. Industrial 3

condiciones de (6) cuando C1 = C2 = 0.

n2 π 2
Como λ = −α2 = − ; el problema (7) queda:
c2

00 n 2 π 2 a2
T (t) + T (t) = 0 T 0 (0) = 0
c2
nπat
que tiene solución T (t) = cos
c
Por lo tanto, son solución todas las funciones de la sucesión:
nπx nπat
yn (x, t) = Cn sen cos n = 1, 2, ... (8)
c c
Se puede demostrar que si todas las funciones del tipo de las de (8) son soluciones de
un problema de contorno homogéneo como el del ejemplo, las combinaciones lineales
de esas funciones también son solución .

2. Principio de superposición de soluciones

Teorema 1 Si y1 , y2 , ..., yn son soluciones de una ecuación diferencial parcial ho-
mogénea, entonces cualquier combinación lineal: C1 y1 + C2 y2 + ... + Cn yn ; donde las
Ci son constantes arbitrarias satisface la ecuación diferencial .

Teorema 2 Si cada una de las funciones de la sucesión infinita y1 , y2 , y3 ... satisfa-

ce la ecuación diferencial parcial homogénea y las condiciones de contorno lineales,
entonces la función y = ∞
P
n=1 Cn yn también la satisface si las constantes Cn y las
funciones yn son tales que la serie converge y es derivable hasta la derivada segunda.

Para el ejemplo desarrollado de la ecuación de la onda:

nπx nπat
Puesto que yn (x, t) = Cn sen cos n = 1, 2, ... son soluciones:
c c
∞ ∞
X X nπx nπat
Cn yn (x, t) = Cn sen cos
n=1 n=1 c c

es solución de la ecuación diferencial (1) y de las condiciones (2) pero debemos en-
contrar las constantes Cn de modo que se verifique (3):
∞ ∞
X X nπx
y(x, 0) = Cn yn (x, 0) = Cn sen = f (x) 0<x<c
n=1 n=1 c

que es la expresión del desarrollo en serie de Fourier de f (x) si la consideramos

como una función periódica, impar, de periodo 2c y en este caso:
∞
X nπx
f (x) = bn sen
n=1 c
FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing. Industrial 4
2Z c nπx
bn = f (x)sen dx
c 0 c
Luego, Cn = bn y:
∞ Z c
2X nπx nπx nπat

y(x, t) = f (x) sen dx sen cos
c n=1 0 c c c
Ejemplos: para dos distintas f (x):
1)
ytt (x, t) = 25yxx (x, t) 0<x<3 t>0 (9)
que satisfaga las condiciones de contorno:
y(0, t) = 0 y(3, t) = 0 yt (x, 0) = 0 (10)
y(x, 0) = 10sen(2πx) 0 < x < 3 (11)
nπx nπ5t
La solución y(x, t) = ∞
P
n=1 Cn sen cos ; como
3 3
P∞ nπx
y(x, 0) = n=1 Cn sen = 10sen(2πx)
3
si n = 6, C6 sen(2πx) = 10 sen(2πx)
Entonces C6 = 10 y Cn = 0 para cualquier n 6= 6
y(x, t) = 10 sen(2πx) cos(10πt)
2)
ytt (x, t) = 25yxx (x, t) 0 < x < 3 t > 0 (12)
que satisfaga las condiciones de contorno:
y(0, t) = 0 y(3, t) = 0 yt (x, 0) = 0 (13)
(
0 0<x<1
y(x, 0) = (14)
x 1<x<3
∞ Z 3
2X nπx nπx nπ5t

y(x, t) = x sen dx sen cos
3 n=1 1 3 3 3

3. Algunos problemas de contorno en problemas

de la Fı́sica

3.1. Ecuación del calor

El problema consiste en conocer la temperatura de un sólido que ocupa una región

R del espacio en un tiempo t: u(x, y, z, t). Teniendo en cuenta el material del
sólido, su densidad, calor especı́fico y conductividad térmica la ecuación diferencial
que representa la temperatura es:
ut = k(uxx + uyy + uzz ) (1)
donde k es el coeficiente conductividad térmica de la sustancia.
FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing. Industrial 5

1. En el caso unidimensional, tenemos la ecuación del calor para el caso de una va-
rilla extendida a los largo del eje x de sección transversal uniforme A (pequeña)
de modo que la temperatura es constante en cada sección.

a) los extremos están a temperatura constante cero grados

El problema de contorno resulta:

 ut
 = kuxx 0 < x < L , t > 0
u(0, t) = u(L, t) = 0


u(x, 0) = f (x) 0 < x < L

b) los extremos están aislados

El problema de contorno resulta:

 ut
 = kuxx 0 < x < L , t > 0
ux (0, t) = ux (L, t) = 0


u(x, 0) = f (x) 0 < x < L

2. El otro caso de la ecuación del calor que consideraremos es cuando las tempe-
raturas están en estado estacionario o sea cuando u no varı́a con el tiempo. En
el caso que el sólido esté en el plano, la ecuación (1) queda:

uxx + uyy = 0 Ecuación de Laplace

Por Ejemplo: Supongamos que se quiere obtener la temperatura u(x, y) corres-

pondiente al estado estacionario en una placa rectangular con las condiciones
de frontera como se indica en la figura.


y6
 uxx
 + uyy = 0 0 < x < a 0 < y < b
u = f (x)
ux (0, y) = 0 ux (a, y) = 0 0 < y < b


u(x, 0) = 0 u(x, b) = f (x) 0 < x < a

(a, b)

aislado aislado
-
u=0 x

3.2. Ecuación de la cuerda vibrante

Es el caso de una cuerda tensamente estirada con posición de equilibrio en algún

intervalo del eje x, que vibra en el plano xy, donde no actúan fuerzas exteriores
y en cada punto la cuerda tiene un desplazamiento pequeño en la dirección del
eje y.
La función y(x, t) (desplazamiento transversal) satisface la ecuación:

ytt = a2 yxx
FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing. Industrial 6

Ejemplo: Un problema tı́pico de contorno es este caso es:

ytt = a2 yxx , 0 < x < L , t > 0



y(0, t) = y(L, t) = 0, t ≥ 0 (los extremos de la cuerda







permanecen fijos en todo instante)


 y(x, 0) = f (x) , yt (x, 0) = g(x) , 0 < x < L
(en t = 0 las funciones f y g representan la posición inicial





y la velocidad inicial de cada punto de la cuerda)


Para resolverlo puede dividirse en dos problemas, cada uno de ellos tiene una condición
de frontera no homogénea, es decir, un problema donde f (x) = 0 y otro
donde g(x) = 0. Si encontramos una solución para cada uno de estos problemas, la
suma de estas soluciones verificará el problema completo.

4. Ejercicios

Para tener en cuenta en la resolución de los problemas mediante el método de

separación de variables:
- Cuando las condiciones de contorno involucran derivadas es posible que exista so-
lución para λ = 0
- Conviene usar todas las condiciones homogéneas antes de usar el principio de su-
perposición.
- Si en algún trozo de la frontera, u es justamente una función que involucra los senos
o cosenos correspondientes a la serie de Fourier, los coeficientes se pueden encontrar
por comparación debido a la unicidad de dicha serie.

2 π2 t
1. Las funciones Xn (x).Tn (t) donde Xn (x) = An sen(nπx) y Tn (t) = Bn e−n
con n ≥ 1, son soluciones del problema de contorno
uxx = ut 0 < x < 1, t > 0
u(0, t) = u(1, t) = 0.

Encontrar u(x, t) aplicando el principio de superposición tal que verifique la

condición u(x, 0) = x, 0 < x < 1.

2. Resolver los problemas de contorno de difusión de calor: calcular la distribu-

ción de la temperatura u(x, t) en una varilla de espesor mı́nimo y superficie en
perfecto aislamiento, en la cual el flujo de calor es unidimensional.
a) 4uxx = ut 0 < x < 3, t > 0
u(0, t) = u(3, t) = 0, t > 0
u(x, 0) = sen( 32 πx) − sen(πx), 0 < x < 3.

b) uxx = ut 0 < x < 2, t > 0

u(0, t) = u(2, t) = 0, t > 0
FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing. Industrial 7
x, 0<x<1

u(x, 0) = .
0, 1<x<2

3. Encontrar la expresión u(x, t) del desplazamiento de una cuerda elástica de lon-

gitud fija en sus extremos.
a) uxx = utt 0 < x < π, t > 0,
u(0, t) = u(π, t) = 0 t > 0,
u(x, 0) = 10−2 sen(3x), ut (x, 0) = 0, 0 < x < π.

b) 9uxx = utt 0 < x < 3, t > 0

u(0, t) = u(3, t) = 0 t > 0,
0, 0<x<1


u(x, 0) =  2, 1 < x < 2 , ut (x, 0) = 0, 0 < x < 3,
0, 2 < x < 3
c) uxx = utt 0 < x < 4, t > 0
u(0, t) = u(4, t) = 0 t > 0,
x, 0 < x < 2

u(x, 0) = sen π4 x, 0 < x < 4, ut (x, 0) = .
4 − x, 2 < x < 4

4. Problemas de la ecuación de Laplace o ecuación del potencial. Dada una lámina

delgada de metal, cuyas caras están perfectamente aisladas, calcular la distribu-
ción en ella de la temperatura de régimen permanente (o del potencial eléctrico)
u(x, y) que satisface las condiciones de contorno dadas.

a) uxx + uyy = 0, 0 < x < 1, 0 < y < 1

ux (0, y) = ux (1, y) = 0 0 < y < 1
u(x, 0) = 0, u(x, 1) = 1 − x, 0<x<1

b) uxx + uyy = 0, 0 < x < 5, 0 < y < ∞

u(0, y) = u(5, y) = 0
u(x, 0) = 10x
lı́my→∞ u(x, y) = 0

c) uxx + uyy = 0, 0 < x < 10, 0 < y < 30

u(0, y) = u(10, y) = u(x, 30) = 0
u(x, 0) = 2x

d) uxx + uyy = 0, 0 < x < ∞, 0 < y < 6

u(x, 0) = u(x, 6) = 0
u(0, y) =sen (πy)
lı́mx→∞ u(x, y) = 0
FI. UNLP. Matemática D1. Módulo II. Unidad 3 primer semestre 2019 Ing. Industrial 8

5. Hallar una función u(x, y) armónica en el interior del rectángulo del plano xy:
{(x, y) : 0 < x < 3, 0 < y < 2} y cuyos valores en el contorno son:
u(0, y) = 0 u(3, y) = 0 0<y<2
u(x, 0) = cos x u(x, 2) = 0 0<x<3

Bibliografı́a:

Churchill, R.V. Series de Fourier y Problemas de Contorno. McGraw-Hill.

Edwards, C.H. Jr., Penney, D. E. Ecuaciones Diferenciales Elementales y Proble-

mas con Condiciones en la Frontera. Prentice-Hall Hispanoamericana S.A.

Zill, D. G. Ecuaciones Diferenciales con Aplicaciones. Grupo Editorial Ibéroaméri-

ca.