ECUACIONES DIFERENCIALES A DERIVADAS PARCIALES
DEFINICIÓN: Se llama ecuación diferencial a derivadas parciales EDP a una ecuación de la forma
𝐹 𝑥1 , … , 𝑥𝑛 , 𝑢, 𝑢𝑥1 , … , 𝑢𝑥𝑛 , … , 𝑢𝑥 𝑘1 𝑥 𝑘2 … 𝑥𝑘𝑛 = 0
1 2 𝑛
𝑚
donde 𝐹: Ω ⊂ ℝ → ℝ, con Ω abierto, 𝑥1 , … , 𝑥𝑛 ∈ Ω son las variables independientes y 𝑢 = 𝑢(𝑥1 , … , 𝑥𝑛 ) es la
variable dependiente.
DEFINICIÓN: El orden de la EDP viene dado por la derivada de mayor orden que aparece en la ecuación.
EJEMPLOS. en la siguiente tabla se presentan algunas EDP junto a su correspondiente orden:
ECUACIÓN EXPRESIÓN ORDEN
Laplace 𝑢𝑥𝑥 + 𝑢𝑦𝑦 = 0 2
Fourier 𝑢𝑡 − 𝛼 2 𝑢𝑥𝑥 = 0 2
Onda 𝑢𝑥𝑥 − 𝑢𝑦𝑦 = 0 2
Euler-Bernoulli 𝑢𝑡𝑡 + 𝛽 2 𝑢𝑥𝑥𝑥𝑥 = 0 4
Schrodinguer 𝑢𝑡 − 𝑖𝛽𝑢𝑥𝑥 = 0 2
Del calor unidimensional 1 2
𝑢𝑥𝑥 = 2 𝑢𝑡𝑡
𝑣
Hemos visto que la solución general de una ecuación diferencial ordinaria de orden 𝑛 contiene 𝑛 constantes
arbitrarias, y que, los valores de esas constantes quedan determinados cuando imponemos apropiadas condiciones
iniciales. La solución general de una ecuación diferencial a derivadas parciales de orden n, por otro lado, contiene
frecuentemente 𝑛 funciones arbitrarias.
EJEMPLO: La ecuación unidimensional del calor
1
𝑢𝑥𝑥 = 2 𝑢𝑡𝑡 ,
𝑣
tiene como solución a
𝑢 𝑥, 𝑡 = 𝑓 𝑥 + 𝑡𝑣 ,
donde 𝑓 es una función arbitraria de una variable 𝑟. Pues,
𝑑2 𝑓 𝑑 2𝑓
𝑢𝑥𝑥 = 2 𝑥 + 𝑡𝑣 𝑢𝑡𝑡 = 𝑣 2 2 𝑥 + 𝑡𝑣 .
𝑑𝑟 𝑑𝑟
La solución general de la ecuación dada es:
𝑢 𝑥, 𝑡 = 𝑓 𝑥 + 𝑡𝑣 + 𝑔 𝑥 + 𝑡𝑣
donde 𝑓 y 𝑔 son funciones arbitrarias.
Método de separación de variables
En dos dimensiones y haciendo cierta suposición sobre sus soluciones es posible resolver algunas clases de EDP. Este
procedimiento es el llamado método de separación de variables donde se supone que la función buscada 𝑢 𝑥, 𝑦 es
de la forma
𝑢 𝑥, 𝑦 = 𝑓 𝑥 𝑔(𝑦)
Con esta suposición se cumple
𝑢𝑥 = 𝑓′ 𝑔, 𝑢𝑦 = 𝑓 𝑔′
𝑢𝑥𝑥 = 𝑓 ′′𝑔, 𝑢𝑥𝑦 = 𝑓 ′ 𝑔′ 𝑢𝑦𝑦 = 𝑓 𝑔′′
⋮
𝑢𝑥…𝑥𝑦…𝑦=𝑓𝑘 𝑔𝑛−𝑘
Sustituyendo estas expresiones en la EDP, en algunas ocasiones es posible reducir una EDP a un sistema de EDO con
dos ecuaciones que puede resolverse con los métodos habituales.
EJEMPLO: Resuelve la siguiente EDP de segundo orden por el método de separación de variables
𝑢𝑥𝑥 = 4𝑢𝑦
Asumiremos por tanto que la solución 𝑢(𝑥, 𝑦) es de la forma
𝑢 𝑥, 𝑦 = 𝑓 𝑥 𝑔 𝑦
Entonces,
𝑢𝑥 (𝑥, 𝑦) = 𝑓 ′ 𝑥 𝑔 𝑦 , 𝑢𝑥𝑥 (𝑥, 𝑦) = 𝑓 ′′ 𝑥 𝑔 𝑦 , 𝑢𝑦 (𝑥, 𝑦) = 𝑓 𝑥 𝑔′(𝑦)
y sustituyendo en la EDP obtendremos
𝑓 ′′ 𝑥 𝑔 𝑦 = 4𝑓 𝑥 𝑔′ 𝑦
Si ahora se agrupan de forma independiente los términos en 𝑥 e 𝑦 obtenemos
𝑓′′(𝑥) 𝑔′(𝑦)
=4
𝑓(𝑥) 𝑔(𝑦)
El primer miembro es una función que sólo depende de 𝑥 y el segundo miembro es una función que sólo depende
de 𝑦, así que la única forma de que se cumpla la igualdad es que ambos valores sean iguales a una constante 𝑐 ∈ ℝ,
que se denomina constante de separación:
𝑓′′(𝑥) 𝑔′(𝑦)
=𝑐=4
𝑓(𝑥) 𝑔(𝑦)
Igualando cada fracción a la constante de separación obtenemos dos ecuaciones diferenciales ordinarias
𝑓 ′′ 𝑥 = 𝑐 𝑓 𝑥 , 4𝑔′ 𝑦 = 𝑐𝑔 𝑦
que se resuelven de forma independiente, la primera es una ecuación EDO lineal de segundo orden y la segunda es
una EDO lineal de primer orden.
Ecuación del calor
La ecuación del calor describe la variación de la temperatura en una región a lo largo del tiempo. En el caso de la
ecuación del calor en una dimensión, describiría la temperatura en una barra de longitud 𝐿 con el tiempo 𝑡 y se
expresa como
𝑢𝑡 = 𝛼𝑢𝑥𝑥 𝑡 > 0, 𝑥 ∈ (0, 𝐿)
൞𝑢 0, 𝑥 = 𝑓 𝑥 , 𝑥 ∈ (0, 𝐿)
𝑢 𝑡, 0 = 𝑢 𝑡, 𝐿 = 0, 𝑡 > 0.
Las condiciones
𝑢 𝑡, 0 = 𝑢 𝑡, 𝐿 = 0
son las condiciones de contorno e indican que la temperatura en los extremos de la barra es constante y es igual a 0.
Mientras que la condición
𝑢 0, 𝑥 = 𝑓 𝑥
es una condición inicial e indica la distribución de temperatura en la barra en el instante inicial. El valor 𝛼 es la
difusividad térmica y depende del material que forma la barra.
Vamos a utilizar el método de separación de variables para resolver esta ecuación, para ello supondremos que la
solución 𝑢 = 𝑢(𝑥, 𝑡) puede ponerse como producto de dos funciones en cada una de las variables independientes:
𝑢 𝑥, 𝑡 = 𝐹 𝑡 𝐺(𝑥)
Derivando y sustituyendo en la ecuación obtenemos
𝐹 ′ 𝑡 𝐺 𝑥 = 𝛼 2 𝐹 𝑡 𝐺′′(𝑥)
Obviamente la función nula 𝑢 = 0 es solución de la ecuación en derivadas parciales, sin embargo sólo será solución
del problema si 𝑓 = 0 que es el caso trivial; así que buscaremos soluciones alternativas y por tanto supondremos
que 𝐺 𝑥 ≠ 0 y 𝐹 𝑡 ≠ 0, por tanto:
1 𝐹′(𝑡) 𝐺′′(𝑥)
=
𝛼 2 𝐹(𝑡) 𝐺(𝑥)
Un lado de la igualdad depende sólo de 𝑡 y el otro depende sólo de 𝑥, por tanto, para que se de la igualdad ambos
miembros deben ser constantes:
1 𝐹′(𝑡) 𝐺′′(𝑥)
= 𝜆 =
𝛼 2 𝐹(𝑡) 𝐺(𝑥)
con 𝜆 ∈ ℝ, la constante de separación y donde hemos elegido el signo menos por convención.
De las últimas igualdades obtenemos dos ecuaciones diferenciales:
𝐹 ′ 𝑡 − 𝜆𝛼 2 𝐹 𝑡 = 0
𝐺 ′′ 𝑥 − 𝜆𝐺 𝑥 = 0
Las condiciones de contorno se transforman en
𝑢 𝑡, 0 = 𝐹 𝑡 𝐺 0 = 0
𝑢 𝑡, 𝐿 = 𝐹 𝑡 𝐺 𝐿 = 0
y siendo 𝑡 es arbitraria, se deduce que:
𝐺 0 =𝐺 𝐿 =0
Tenemos así el siguiente problema de valores iniciales
𝐺 ′′ 𝑥 − 𝜆𝐺 𝑥 = 0
ቐ 𝐺 0 =0
𝐺 𝐿 =0
cuya solución dependerá del valor del parámetro de separación 𝜆.
1. Caso 𝝀 = 𝟎. En este caso la ecuación diferencial sería
𝐺 ′′ 𝑥 = 0
cuya solución general se obtiene integrando dos veces respecto a 𝑥,
𝐺 𝑥 = 𝐴𝑥 + 𝐵
Si utilizamos las condiciones de contorno
𝐴0+𝐵 =0 𝐵=0
൝ ⇒ቊ
𝐴𝐿+𝐵 =0 𝐴=0
pues 𝐿 ≠ 0, entonces obtenemos la solución nula, 𝐺 𝑥 = 0.
2. Caso 𝝀 > 𝟎. Supongamos ahora que 𝜆 es positivo. La ecuación diferencial sería
𝐺 ′′ 𝑥 − 𝜆𝐺 𝑥 = 0
que tiene por solución general
𝐺 𝑥 = 𝐴𝑒 𝜆𝑥 + 𝐵𝑒 − 𝜆𝑥
Utilizamos las condiciones iniciales para determinar los valores de 𝐴 y 𝐵:
𝜆0 + 𝐵𝑒 − 𝜆0 = 0 𝐴 = −𝐵 𝐴=0
൝𝐴𝑒 ⇒൝ ⇒ ቊ
𝐴𝑒 𝜆𝐿 + 𝐵𝑒 − 𝜆𝐿 = 0 𝐴 = −𝐵𝑒 −2 𝜆𝐿 𝐵=0
pues 𝑒 −2𝜇𝐿 ≠ 1, que nos da para el PVI de nuevo la solución nula 𝐺(𝑥) = 0.
3. Caso 𝝀 < 𝟎. Supongamos ahora que 𝜆 es positivo, y por tanto lo podemos poner de la forma 𝜆 = − 𝜆 . La
ecuación diferencial sería
𝐺 ′′ 𝑥 + 𝜆2 𝐺 𝑥 = 0
que tiene por solución general
𝐺 𝑥 = 𝐴𝑐𝑜𝑠 𝜆 𝑥 + 𝐵𝑠𝑒𝑛 𝜆 𝑥
Utilizamos las condiciones de contorno para encontrar los valores de 𝐴 y 𝐵
𝐴𝑐𝑜𝑠 𝜆 0 + 𝐵𝑠𝑒𝑛 𝜆0=0 𝐴=0
ቐ ⇒൝
𝐴𝑐𝑜𝑠 𝜆 𝐿 + 𝐵𝑠𝑒𝑛 𝜆𝐿=0 𝐵𝑠𝑒𝑛 𝜆 𝐿 = 0
𝑛2 𝜋 2
Como no queremos la solución nula debe ocurrir 𝐵 ≠ 0 y 𝑠𝑒𝑛 𝜆𝐿=0⇔𝜆= − 𝐿2 . Luego, para cada valor de
n ∈ ℕ, tendremos una posible solución de la EDO,
𝑛2 𝜋 2
𝐺𝑛 𝑥 = 𝐵𝑛 𝑠𝑒𝑛 𝑥 .
𝐿2
Los únicos valores que proporciona soluciones distintas de la solución nula son
𝑛2 𝜋 2
𝜆𝑛 = − 2
𝐿
Para estos valores y utilizando la otra ecuación diferencial
′
𝑛2 𝜋 2 2
𝐹 𝑡 + 2 𝛼 𝐹 𝑡 =0
𝐿
cuya solución para cada n ∈ ℕ es de la forma
𝛼 2 𝑛 2 𝜋2
− 𝑡
𝐹𝑛 𝑡 = 𝐴𝑛 𝑒 𝐿2
Y la solución de la EDP será, para cada n ∈ ℕ, de la forma
𝛼 2 𝑛2 𝜋 2 𝑛𝜋
− 𝑡
𝑢𝑛 𝑥, 𝑡 = 𝐶𝑛 𝑒 𝐿2 𝑠𝑒𝑛 𝑥
𝐿
donde hemos puesto 𝐶𝑛 = 𝐴𝑛 𝐵𝑛 .
Utilizando algunos herramientas de la teoría de series de Fourier, la condición inicial 𝑢 0, 𝑥 = 𝑓 𝑥 determinará los
valores de 𝐶𝑛 para todo n ∈ ℕ.
