FORMAS NORMALES

ECUACIONES DIFERENCIALES PARCIALES

E. SAEZ

La idea general de la Teora de Formas Normales consiste en encontrar un sistema

de coordenadas donde una determinada expresion se reduce a una expresion equiva-
lente que contiene solo los terminos mas relevantes, por ejemplo, es conocido en el
estudio de las cuadraticas reales en dos variables reales con coeficientes constantes,
que existen traslaciones y rotaciones en el plano que permiten definir nuevas coorde-
nadas donde la ecuacion de la cuadratica no contiene el termino mixto y no contiene
al menos un termino lineal ya que son irrelevantes para la forma cualitativa de la
grafica de la cuadratica. Dichas ecuaciones son llamadas Canonicas. Otro ejemplo
conocido es en el estudio de las matrices cuadradas, donde se demuestra que existen
bases adecuadas en los Espacios Vectoriales respectivos tales que mediante transfor-
maciones lineales es posible representar la matriz en terminos de los valores propios
que son los n umeros realmente relevantes de la matriz. Dichas matrices se dicen de
Jordan.
La expresion general de una Ecuacion Diferencial Parcial Lineal Real de Segundo
Orden en dos variables, definida en un dominio del plano en coordenadas cartesianas
es de la forma
2 2 2
a(x, y) xu2 + 2b(x, y) xy
u
+ c(x, y) yu2 + d(x, y) u + e(x, y) u +
(1) x y
f (x, y)u + g(x, y) = 0
donde los coeficientes a, b, c : R son funciones en dos variables que admiten
desarrollos de Taylor convergentes y no se anulan simultaneamente en . Los coefi-
cientes d, e, f, g : R son funciones en dos variables y continuas.
Como la ecuacion (1) es de segundo orden, veremos en lo que sigue que siempre es
posible reducir los coeficientes de las derivadas de segundo orden a constantes muy
simples mediante un cambio de coordenadas definidas por sistema de ecuaciones de
la forma
= (x, y)

(2) (,)
= (x, y) con (x,y)
6= 0 , (x, y)
tal que (1) en las nuevas coordenadas es equivalente a una de los siguientes tipos
de ecuaciones mas sencillas, llamadas Formas Normales , o bien , formas Formas

Departamento de Matem atica, UTFSM

email: [email protected].
1
2
E. SAEZ

onicas de (1).
Can

2u 2u

1.a) 2

2
+ T.O.I. =0
2u

1.b) + T.O.I. =0




2u
(3) FormasNormales : 2.a) 2
+ T.O.I. =0
2u




2.b) 2
+ T.O.I. =0
3) 2u 2u
+ 2 + T.O.I. =0

2

donde T.O.I., designa los terminos de orden inferior al efectuar el cambio de coorde-
nadas (2) a (1) para obtener (3).
Definicion. Diremos que la ecuacion (1) es de tipo Hiperb olica,
o Parabolica,
o bien, Elptica, si y solo si existe un cambio de coordenadas tal que la ecuacion se
puede escribir en la Forma Normal 1., o 2., o bien 3., respectivamente.
Notese de (3) que no existe unicidad de Formas Normales. Ademas, las formas 1.),
2.) y 3.) son por definicion de tipo Hiperbolicas, Parabolicas y Elptica, respectiva-
mente.
En Fsica e Ingeniera es frecuente encontrar las E.D.P.
2
u 2u
x2 t2 = 0 Ecuacion de Onda Unidimensional

2u u
x 2 t = 0 Ecuacion de Calor Unidimensional
u2 + 2 u2 = 0 Ecuacion de Laplace Bidimensional
2
x y

Estas ecuaciones son ejemplos inmediatos de Ecuaciones de tipo Hiperbolica, Parabolica

y Elptica, respectivamente.
Teorema
La E.D.P. (1) es reducible en a la forma:
1) Hiperbolica si y solo si b(x, y)2 a(x, y)c(x, y) > 0 , (x, y)
2) Parabolica si y solo si b(x, y)2 a(x, y)c(x, y) = 0 , (x, y)
3) Elptica si y solo si b(x, y)2 a(x, y)c(x, y) < 0 , (x, y)
Demostraci
on
Sea : R tal que (x, y) = b(x, y)2 a(x, y)c(x, y). Demostraremos
en primer lugar que en , el signo del discriminante (x, y) es invariante bajo un
cambio de coordenadas del tipo (2). Por (2) es inmediato que existe una cadena de
dependencia de las variables del tipo:

x
u
y
 FORMAS NORMALES DE ECUACIONES ... 3

Usando la Regla de la Cadena y cambiando de notacion se obtiene que las derivadas

parciales respecto de x e y, estan dadas por las expresiones:


ux = u x + u x
uy = u y + u y


uxx = u x2 + 2u x x + u x2 + T.O.I.


uxy = u y x + u y x + u y x + u x y + T.O.I.
2 2

yy = u y + 2u y y + u y + T.O.I.
u

Reemplazando las derivadas parciales anteriores en (1) se obtiene que la forma general
de la Ecuacion Parcial en las cordenadas (, ) es del tipo
(4) Au + 2Bu + Cu + T.O.I. = 0
donde los coeficientes de las derivadas de segundo orden estan dadas por:

A = ax2 + 2bx y + cy2
(5) B = ax x + b(x y + y x ) + cy y
C = a 2 + 2b + c 2
x x y y

Con un calculo operatorio simple se verifica la identidad:

 2
2 2 2 2 (, )
B AC (b ac)(x y y x ) (b ac)
(x, y)
La identidad anterior, demuestra que el signo de la expresion = b2 ac es invariante
bajo cambio de coordenadas pues
 2
2 2 (, )
Sing(B AC) = Sing(b ac) si y solo si 6= 0
(x, y)

(). Es inmediato de (3) que las Formas Normales : Hiperbolica, Parabolica y

Elptica implican, respectivamente: (x, y) 1(resp. 41 ), (x, y) 0 y (x, y)
1, de donde los signos de (x, y) en son positivos, cero y negativo ya que son
invariantes bajo cambio de coordenadas.
() Demostraremos a continuacion que dependiendo del signo de , existe un
cambio de coordenadas (2) que transforma (1) a una de las respectivas Formas Nor-
males en (3).
Caso 1.- Supongamos que (x, y) > 0 en y busquemos funciones = (x, y), =
(x, y) tales que anulen simultaneamente las expresiones de A y C en (5). Si tales
1
funciones existen entonces B 6= 0 en (5) pues B 2 AC > 0 y dividiendo (4) por 2B
se tiene la Forma Normal 1.b) de (3).
La demostracion se reduce a analizar la consistencia del sistema de ecuaciones
diferenciales de primer orden
ax2 + 2bx y + cy2 = 0
(6)
ax2 + 2bx y + cy2 = 0
4
E. SAEZ

Como las ecuaciones del sistema son similares, el estudio de la primera de ellas es
valido para la segunda ecuacion cambiando el rol de por .
La primera ecuacion del sistema, dividida por y2 tiene la forma
 2
x x
(7) a + 2b + c = 0
y y
Si a 6= 0 en cada punto de , se tiene una ecuacion de segundo grado donde el
cuociente xy juega el rol de incognita. Sus soluciones son:

x b + b2 ac x b b2 ac
= , o bien, =
y a y a
Sean las ecuaciones diferenciales ordinarias

dy b b2 ac dy b + b2 ac
(8) = , y, =
dx a dx a
Si (x, y) = cte es una solucion general de la segunda ecuacion anterior, entonces es
dy
inmediato que x dx + y dy = 0 y en consecuencia dx = xy , de donde el cuociente
x
y
es solucion de la ecuacion de segundo grado (7) y la funcion satisface la primera
ecuacion del sistema (6). Analogamente si (x, y) = cte es solucion general de la
primera ecuacion de (8), entonces la funcion satisface la segunda ecuacion del
sistema (6).
Para demostrar que la transformacion

= (x, y)
= (x, y)
define un cambio coordenadas en , basta demostrar que el jacobiano de la transfor-
macion no se anula en . En efecto las funciones y , por (7) y (8) satisfacen
(
2
xy = b+ ba ac

xy = b b2 ac
a

Luego

x x b2 ac x y y x
y
y
= 2 a
6= 0 , de donde, y y
6= 0
(,)
Entonces, x y y x = (x,y)
6= 0 , x, y
Si a 0, entonces b(x, y) 6= 0 en cada punto de pues por hipotesis (x, y) > 0. Si
c 0, basta dividir (1) por 2b para tener la forma normal de la ecuacion diferencial.
Sea entonces c(x, y) 6= 0 en cada punto de . El sistema (6) se reduce a la expresion
2bx y + cy2 = 0
2bx y + cy2 = 0
 FORMAS NORMALES DE ECUACIONES ... 5

Dividiendo la primera ecuacion del sistema por x2 se tiene la ecuacion de segundo

grado con incognita el cuociente xy ;
y y
2b + c( )2 = 0
x x
o bien factorizando ,
y y
(2b + c ) = 0
x x
Sus soluciones son :
y y 2b
= 0 , o bien , =
x x c
Sean las ecuaciones diferenciales ordinarias
dx dx 2b
=0 , y, =
dy dy c
Sean x = cte , y , (x, y) = cte las respectivas soluciones generales de las ecuaciones
ordinarias, entonces

= x
= (x, y)
satisface el sistema (6). El sistema como una transformacion define un cambio de
(,)
coordenadas pues (x,y) = y 6= 0 , ya que, b 6= 0. Lo anterior demuestra solo la parte
1) del Teorema.
Caso 2.- Supongamos que (x, y) 0 en . Analogamente a lo anterior busque-
mos ahora funciones = (x, y), = (x, y) tales que anulen simultaneamente las
expresiones de A y B en (5).
Si a 0 en (1) entonces b 0 pues b2 ac 0, luego c 6= 0 en cada punto de .
Basta dividir (1) por c para tener la forma Parabolica 2.b).
Si a 6= 0 en cada punto de en (1), del caso 1.- se sabe que la funcion = (x, y)
obtenida de la segunda ecuacion diferencial en (8), (que se reduce a la forma mas
dy
simple dx = ab ), tiene la propiedad de anular la expresion de A . Pero B 2 AC 0
pues el sig((x, y)) es invariante bajo cambios de coordenadas, entonces la misma
funcion = (x, y) anterior, tambien anula B. Por otro lado basta tomar (x, y) = y
ya que en el caso c 6= 0 en cada punto de , no anula la expresion de C (si c 0
el teorema es inmediato pues (1) dividida por a se encuentra en la forma normal
Parabolica 2.b))
Para terminar la demostracion de este caso falta demostrar que la transformacion
definida por el sistema:

= (x, y)
= y , define un cambio de coordenadas en .
6
E. SAEZ

Pero
(, ) x y x b
= = x , y =
(x, y) 0 1 y a
Si b 0 en , entonces c 0 y la demostracion es inmediata pues basta dividir (1)
por a para tener la forma normal Parabolica 2.b).
Si b 6= 0 en cada punto de , entonces x 6= 0 y el Jacobiano anterior no se anula
en lo que termina la demostracion del caso 2.-.
Caso 3.- Supongamos que (x, y) < 0 en . Busquemos funciones = (x, y), =
(x, y) tales que simultaneamente anule la expresiones de B y se tenga la identidad
A C en (5). Entonces de (5) y de las identidades en , A C 0 y B 0 se tiene
respectivamente el sistema:
a(x2 x2 ) + 2b(x y x y ) + c(y2 y2 ) = 0

(9)
ax x + b(x y + y x ) + cy y = 0
Multiplicando la segunda ecuacion del sistema por 2i y sumando con la primera se
obtiene:
a(x + ix )2 + 2b(x + ix )(y + iy ) + c(y + iy )2 = 0
Dividiendo la ecuacion por (y + iy )2 se tiene:
x + ix 2 x + ix
a( ) + 2b( )+c=0
y + iy y + iy
x +ix
Esta expresion significa que el cuociente y +iy
es solucion de la ecuacion de segundo
grado aX 2 + 2bX + c = 0, es decir,

x + ix b i ac b2
(10) =
y + iy a
Notese que a 6= 0 en cada punto de pues b2 ac < 0 (el caso a = 0 no es posible
pues b2 < 0 es una contradiccion .)
Consideremos la ecuacion diferencial ordinaria

dy b + i ac b2
= , y sea , (x, y) = (x, y) + i(x, y) = c1 + ic2 C
dx a
su solucion general, donde c1 + ic2 es una constante compleja arbitraria.
dy
Pero la diferencial de (x, y) = cte es d = x dx + y dy = 0, entonces dx
=
y
x

y en consecuencia se tiene (10).

El calculo anterior significa que las funciones

= (x, y)
(11)
= (x, y)
satisfacen el sistema (9) y en consecuencia A C B 0.
 FORMAS NORMALES DE ECUACIONES ... 7

Falta demostrar que el sistema (11) define una transformacion que es un cambio
(,)
de coordenadas, o bien , que (x,y) 6= 0 en cada punto de .
Consideremos la parte imaginaria

x + ix y x x y ac b2
Im( )= = 6= 0 , en cada punto de
y + iy y2 + y2 a
Entonces
(, )
= x y y x 6= 0 , en cada punto de
(x, y)
lo que concluye la demostracion.
Comentario 1. El teorema anterior se generaliza a E.D.P. lineales de segundo
orden definidas en un dominio Rn . Por cambios de coordenadas la forma general
de la E.D.P. se puede escribir en la forma normal
n
X 2u
i 2 + T.O.I. = 0 , y se tiene la siguiente clasificacion :
i=1
i

1) Elptica i, i > 0 (resp. i < 0)

2) Parabolica i, i = 0
3) Hiperbolica (!)i, i > 0 (resp. i < 0) , j < 0 (resp. j > 0) si i 6= j
4) Ultrahiperbolica Existe mas de un i > 0 y mas de un i < 0
Comentario 2. Sea la funcion asociada a una misma E.D.P. de la forma (1)
: R2 R, tal que , (x, y) = b2 (x, y) a(x, y)c(x, y)
Consideremos la particion de en los tres conjuntos
= 1 2 3 donde 1 = 1 (, 0), 2 = 1 (0) y 3 = 1 (0, )
Entonces, una misma ecuacion puede ser de los tres tipos: Elptica , Parabolica, o
bien, Hiperbolica si se restringe a los conjuntos 1 , 2 , o bien , 3 respectivamente.
Es decir, el tipo de ecuacion depende del dominio de definicion de los coeficientes de
la E.D.P.
Ejercicio interesante
Encontrar la solucion de la ecuacion de Onda Unidimensional
2u 2
2 u
=k donde k > 0 es una constante,
t2 x2
tal que satisfaga las condiciones iniciales de Cauchy:
u(x, 0) = f (x)
(12) u con f, g C 1
t
(x, 0) = g(x)

Solucion: Si el coeficiente k 6= 1 en la ecuacion de Onda, entonces dicha ecuacion

no esta escrita en forma normal por la presencia de este coeficiente diferente de uno .
8
E. SAEZ

Es inmediato que, (x, y) = b2 ac = k 2 > 0 y en consecuencia la ecuacion de Onda

es de tipo Hiperbolico.
Sean las ecuaciones diferenciales ordinarias

dt b + b2 ac 1 dt b b2 ac 1
= = , = =
dx a k dx a k
Entonces las respectivas soluciones generales son:
x kt = c1 , x + kt = c2 , donde c1 , c2 son constantes arbitrarias.
Luego existe el cambio de coordenadas
= x kt
, donde la Ec. de Onda es de la forma Hiperbolica
= x + kt
De (5) se tiene que B = 2k 2 y la ecuacion en las nuevas coordenadas es de la forma
2u
2B + T.O.I = 0. Es simple verificar que los terminos T.O.I 0, dividiendo la
2
u
ecuacion anterior por 4k 2 se tiene la forma hiperbolica = 0. Es claro que cualquier
funcion que depende solo de unas de las variables es anulada por la derivada mixta
(se supone de clase C 2 ). Entonces
u = () + () , con , C 2 funciones arbitrarias
son soluciones de la forma hiperbolica. En consecuencia, regresando a las variables
originales por el cambio de coordenadas ,
(13) u(x, t) = (x kt) + (x + kt)
son soluciones arbitrarias de la ecuacion de Onda.
Consideremos la siguiente interpretacion de la solucion u = (xkt) . Supongamos
que la grafica de en el instante inicial t = 0 es como la curva 0 de la figura. Entonces
la grafica de u = (x kt) es una curva dinamica t que como funcion del tiempo, se
desplaza s = kt unidades hacia la derecha a velocidad k, pues kt > 0. Analogamente
la grafica de (x + kt) es una curva que como funcion del tiempo se desplaza s = kt
unidades a velocidad k hacia la izquierda pues kt < 0. Por estas interpretaciones
, se dice que las soluciones de la ecuacion de Onda, u(x, t) = (x kt) + (x + kt)
son sumas de ondas viajeras como indica la figura:
u
0 t

kt
x
 FORMAS NORMALES DE ECUACIONES ... 9

Busquemos ahora las funciones , tales que la(s) solucion(es) satisface(en) las condi-
ciones iniciales de Cauchy, es decir, la(s) solucion(es) satisfacen las condiciones (12),
entonces tomando t = 0 en (13) se tiene que , satisfacen el sistema:

f (x) = (x) + (x)
g(x) = k (x) + k (x)
Derivando la primera ecuacion del sistema anterior respecto de x se tiene que las
derivadas de , satisfacen el sistema:

(x) + (x) = f (x)

k (x) + k (x) = g(x)
Despejando del sistema las derivadas de , tenemos:
(
(x) = kf (x)g(x)
2k

(x) = kf (x)+g(x)
2k

Integrando respecto de x se tiene las identidades:

Rx
(x) f (x) 1

2
2k R0
g()d
f (x) 1 x
(x) 2 + 2k 0 g()d
Reemplazando x x kt y x x + kt en la primera y segunda identidad respecti-
vamente se obtiene
( R xkt
(x kt) f (xkt)
2
1
2k g()d
f (x+kt) 1
R0x+kt
(x + kt) 2
+ 2k 0 g()d
Finalmente, sumando las dos identidades anteriores encontramos la solucion llamada
de DAlambert del problema de la Ecuacion de Onda con condiciones iniciales de
Cauchy.

1 1
R x+kt
u(x, t) = 2
[f (x + kt) + f (x kt)] + 2k xkt
g()d

RESUMEN
Caso 1.-
b2 (x, y) a(x, y)c(x, y) > 0, (x, y) E.D.P. Hiperbolica.
i) Si a 6= 0, sean las Ecuaciones Diferenciales Ordinarias:

dy b + b2 ac dy b b2 ac
= , =
dx a dx a

(x, y) = c1
y las respectivas soluciones generales :
(x, y) = c2
10
E. SAEZ

Entonces el cambio de coordenadas:

2u

= (x, y)
, Transforma (1) a la forma , 2B + T.O.I. = 0
= (x, y)
Dividiendo la expresion anterior por 2B , se obtiene la Forma Normal
Hiperbolica,

2u
(14) + T.O.I. = 0

ii) Si a = 0 y c 6= 0 , en cada punto de , sean las Ecuaciones Diferenciales

Ordinarias:
dx dx 2b
=0 , =
dy dy c

x = c1
y las respectivas soluciones generales :
(x, y) = c2
Entonces el cambio de coordenadas:
2u

= x
, Transforma (1) a la forma , 2B + T.O.I. = 0
= (x, y)
Dividiendo la expresion anterior por 2B , se obtiene la Forma Normal
Hiperbolica ().
iii) Si a = 0 y c = 0 , en cada punto de , basta dividir (1) por 2b para obtener
la Forma Normal Hiperbolica en las coordenadas cartesianas originales,

2u
+ T.O.I. = 0
xy

Caso 2.-
b2 (x, y) a(x, y)c(x, y) 0, (x, y) E.D.P. Parabolica.
i) Si a = 0, en cada punto de , basta dividir (1) por c para obtener la
Forma Normal Parabolica en las coordenadas cartesianas originales,

2u
+ T.O.I. = 0
y 2

ii) Si a 6= 0 y c 6= 0 , en cada punto de , sea la Ecuacion Diferencial

Ordinaria
dy b
= ,
dx a
y la respectiva solucion general: (x, y) = c1
 FORMAS NORMALES DE ECUACIONES ... 11

Entonces el cambio de coordenadas:

2u

= (x, y)
, Transforma (1) a la forma , C 2 + T.O.I. = 0
= y
Dividiendo la expresion anterior por C , se obtiene la Forma Normal
Hiperbolica ,

2u
+ T.O.I. = 0
2

iii) Si a 6= 0 y c = 0 , en cada punto de , basta dividir (1) por a para obtener

la Forma Normal Parabolica en las coordenadas cartesianas originales,

2u
+ T.O.I. = 0
x2

Caso 3.-
b2 (x, y) a(x, y)c(x, y) < 0, (x, y) E.D.P. Elptica.
Sea la Ecuacion Diferencial Ordinaria Compleja:

dy b + i ac b2
= ,
dx a
y la respectiva solucion general : (x, y) = (x, y) + i(x, y) = c1 + ic2 C
Entonces el cambio de coordenadas:
2u 2u

= (x, y)
, Transforma (1) a la forma , A 2 + A 2 + T.O.I. = 0
= (x, y)
Dividiendo la expresion anterior por A , se obtiene la Forma Normal Elptica,

2u 2u
+ 2 + T.O.I. = 0
2

Ejercicios propuestos:
1) Clasificar y reducir a su forma normal:
i) uxx + 2uxy + uyy + ux uy = 0
ii) uxx + 2uxy + 5uyy + 3ux + u = 0
2 2
iii) y 2 uxx + 2xyuxy 3x2 uyy = yx ux 3 xy uy
iv) x2 uxx y 2 uyy = 0
2) Demostrar que terminos de primer orden de una E.D.P., se pueden anular
bajo la sustitucion: u = veAx+By .
12
E. SAEZ

3) Cual es la solucion de la E.D.P.

uxx uyy = 3x2 2y
u(x, 0) = 41 x2

tal que satisface las condiciones: u
y
(x, 0) = x
4) Sea u = u(x, t) y considere la E.D.P.
2u u x 2u
x + 2 = , x 6= 0 , c > 0.
x2 x c2 t2
i) Introduzca en la E.D.P. la sustitucion v(x, t) = xu(x, t).
ii) Escribir la E.D.P. obtenida en i) en su forma normal.
ii) Encontrar la solucion de la 
E.D.P planteada, tal que satisface las condi-
u(x, 0) = x
ciones iniciales de Cauchy: u
t
(x, 0) = 2c.
5) Sea la ecuacion de Onda Unidimensional
2z 2z
= , con < x < , t > 0
t2 x2
i) Cual es la solucion que satisface las condiciones iniciales de Cauchy ?
z(x, 0) = 1 x2
z
t
(x, 0) = 0
ii) En el plano xz haga un bosquejo de la solucion para t = 0 y t = 1.
iii) En el plano xz. Para que tiempo t > 0 la solucion encontrada en i)
pasa por el punto (x, z) = (0, 4) ?. Haga un bosquejo de la onda.
iv) Existe algun instante t tal que la solucion encontrada en i) pasa, en el
plano xz, por el punto (x, z) = (1, 21 ).
6) Sea el problema de Cauchy
2u 2u
t2
= x2 
; x, t R
sen x si x [0, ]
u(x, 0) =
0 si otro caso
u
t
(x, 0) 0
Haga un bosquejo en el plano ux, de la onda solucion del problema, en el
instante t = .

References
[1] A. R. Castro F. Curso b
asico de ecuaciones en derivadas parciales . Addison-Wesley Iberoame-
ricana . Wilmington, Delaware. E.U.A. 1997.
[2] I. Peral A. Primer curso de ecuaciones en derivadas parciales . Addison-Wesley / Universidad
Autonoma de Madrid. Wilmington, Delaware. E.U.A. 1995.