Capitulo3 Espsacios Hermitianos
Capitulo3 Espsacios Hermitianos
Capitulo3 Espsacios Hermitianos
Operadores Hermitianos
En la actualidad la teora de series de Fourier puede presentarse usando los conceptos y
metodos del analisis funcional, los espacios de Hilbert (extension a dimension infinita de
la nocion de espacio Eucldeo) y los operadores compactos y autoadjuntos (extension, a
dimension infinita, de los endomorfismos de R definidos por matrices simetricas). Ello
permite, por una parte, comprender mejor los metodos de Fourier, por otra, se consigue,
sin apenas esfuerzo, una gran generalidad. Los operadores Hermitianos son u
tiles porque
dan lugar a funciones ortogonales que pueden usarse como base.
Definici
on 3.1 (Operador Hermitiano) Un operador O : M 7 M se dice Hermitiano si
para toda f M cumple
hf, Ogi = hOf, gi .
(3.1)
Teorema 3.1 Los valores propios de un operador Hermitiano O son reales.
n.
Demostracio
Entonces tenemos
hf, Of i = hf, f i = ||f ||.
Por otro lado, como es Hermitiano,
hf, Of i = hOf, f i = hf, f i = ||f ||,
de donde = .
Teorema 3.2 Sea O Hermitiano y f (~x), g(~x) dos funciones propias correspondientes a
los valores propios y respectivamente. Si 6= entonces hf, gi = 0, es decir son
ortogonales.
n.
Demostracio
Entonces, si 6= tenemos
3.1
hf, gi = hf, gi ,
( ) hf, gi = 0.
hf, gi = 0.
Problema de Sturm-Liouville
Definici
on 3.2 (Operador de Sturm-Liouville) Llamamos operador de Sturm-Liouville
(y lo denotamos T ) al operador que act
ua sobre la funcion real f (x), con = (a, b),
seg
un
df (x)
1 d
R(x)
,
T f (x) =
(x) dx
dx
donde : 7 R+ y R : 7 R+ .
Teorema 3.3 El operador de Sturm-Liouville T es Hermitiano si (x) es el peso w(x)
del producto escalar y R(a) = R(b).
d
tenemos
dx
Z
w(x)
y
Z
w(x)
hT f, gi =
dx
(R(x)f (x)) g(x).
(x)
Denotando
donde las distintas funciones propias fn (x) asociadas a los distintos valores propias reales
n son ortogonales entre s con peso w(x) = (x), es decir
hfm , fn i = mn (||fm ||)2 mn hm .
3.2
(3.3)
Vamos a estudiar algunas propiedades generales de los polinomios ortogonales, que son
un conjunto de polinomios de orden n (n = 0, 1, 2, . . .) linealmente independientes entre
s en una cierta region del espacio. Para encontrar estos polinomios se pueden usar muy
diversas tecnicas. La mas habitual es partir de unas ecuaciones diferenciales concretas.
Aqu emplearemos un camino diferente (y complementario a ese) mas general que
el anterior, ya que nos permite determinar los casos que usualmente se estudian a partir de ecuaciones diferenciales concretas a partir de caractersticas en cierto modo mas
universales.
Algunas de las funciones especiales dadas aqu no estan normalizadas de la forma mas
usual. Conviene comprobar la normalizacion a la hora de comparar con otras fuentes.
Antes de estudiar la forma de las soluciones al problema de valores propias y funciones
propias expondremos algunas propiedades generales de los polinomios ortogonales.
P
(n)
Denotaremos a partir de ahora a los polinomios ortogonales como n = ni=0 i xi
y consideramos que no estan normalizados, es decir
hm , n i = mn hm .
(3.4)
(3.5)
Partiendo de la forma
para el problema de Sturm-Liouville se puede demostrar que existe un conjunto de soluciones n (x), n N tal que cada funcion n (x) es un polinomio de orden n. Esto se
cumple (a grandes rasgos) si Q(x) es como mucho un polinomio cuadratico y L(x) una
funcion lineal (observese que en ese caso cada termino de la ecuacion es un polinomio y
el grado de estos es consistente). De este modo tenemos
(q0 + q1 x + q2 x2 )n (x) + (l0 + l1 x)n (x) + n n (x) = 0.
De la expresion polinomial de n (x) obtenemos
n
X
(n)
n (x) =
j xj ,
j=0
n (x) =
n1
X
j=0
(n)
(j + 1)j+1 xj ,
n (x)
n2
X
(n)
(j + 1)(j + 2)j+2 xj ,
j=0
y por lo tanto
n2
n1
n
X
X
X
(n) j
(n) j
(n)
(q0 + q1 x + q2 x )
(j + 1)(j + 2)j+2 x + (l0 + l1 x)
(j + 1)j+1 x + n
j xj = 0.
2
j=0
j=0
j=0
(3.6)
3.3
Polinomios representativos
Para evitar la degeneracion de los valores propias es necesario que la funcion L(x) sea lineal. Dependiendo de si esta funcion es efectivamente cuadratica, lineal o simplemente una
constante tendremos tres tipos de polinomios distintos, con sus caractersticas peculiares.
Aqu desarrollaremos los polinomios de tipo Hermite, Laguerre y Gegenbauer a partir
de la exposicion general dada hasta ahora. En la tabla 3.1 mencionamos algunos de los
polinomios ortogonales correspondientes a distintas funciones de peso usuales.
Peso
Soporte
Denominacion
Notacion
(1 x2 )1/2
[1, 1]
Tn (x)
(1 x2 )1/2
[1, 1]
Un (x)
[1, 1]
Legendre
Pn (x)
[1, 1]
Gegenbauer (o ultraesfericos)
Pn (x)
[1, 1]
Jacobi
[0, )
Laguerre
Ln (x)
Hermite
Hn (x)
(1 x2 )1/2 ,
(1 x) (1 + x) ,
> 1/2
, > 1
ex x , > 1
ex
()
(,)
Pn
(x)
()
3.3.1
(3.7)
que son los llamados polinomios de tipo Hermite. Como Q(x) = q0 y L(x) = l0 + l1 x
tenemos que
l
1
A 1 (l + l1 x)x
(l + 1 x)x
w(x) = e q0 0 2 ,
(3.8)
R(x) = A e q0 0 2 ,
q0
donde A > 0 es una constante de integracion que fijamos en 1 por comodidad. Vemos
que la funcion peso esta definida positiva para todo el intervalo < x < , que sera
el intervalo en el que los polinomios tipo Hermite estaran bien definidos.
Con este intervalo de validez (y dado que l1 6= 0) podemos hacer una traslacion y una
1
dilatacion x q02 (x
forma
l0
)
l1
(3.9)
que es la forma general para la ecuacion de este tipo de polinomios. De este modo
l1 2
w(x) = e 2 x y tenemos que exigir que l1 < 0 para que el producto interior de un valor
finito.
Ademas se tiene que R() = R() con lo que se cumplen las condiciones para
que los polinomios Hn (x) sean ortogonales con valor propio real, siendo este valor propio
n = l1 n > 0.
(3.10)
(n)
j<n =
(j + 1)(j + 2) (n)
j+2
(j n)l1
del poli-
(3.11)
(n)
1
2l1
k
n!
cn ,
(n 2k)!k!
(3.12)
de donde
Hn (x) =
n
X
p=0
p(n) xp
[n/2]
X
k=0
c2
2l1
k
n!
(cx)n2k .
(n 2k)!k!
(3.13)
Con esta u
ltima ecuacion se puede obtener la funcion generatriz de los polinomios. As,
tenemos
[n/2] 2 k
X
X c
tn
n!
tn =
(cx)n2k
n!
2l1
(n 2k)!k!
n!
n=0
k=0
!
k ! X
X
(ct)2
1
(cxt)n
=
2l1
k!
n!
n=0
k=0
c2 t2
.
= exp cxt +
2l1
X
Hn (x)
n=0
(3.14)
Con este resultado podemos obtener la formula de Rodrigues, que es una formula que nos
permite obtener el polinomio tipo Hermite de grado n sin mas que derivar n veces. Como
f (x) =
X
f (n) (0)
n=0
entonces
n!
2 2
cxt+ c2lt
x y e
X
Hn (x)
n=0
n!
tn ,
(3.15)
n
2
c2 l1 x
l1 x 2
n cxt+ c2l2 t2
+t
)
(
2
1
e 2l1 c
Hn (x) = n e
=e
n
t
t
t=0
t=0
n
n
n
c2 2
l1 x 2
l1 x 2
l1 x 2
c
z
2
2l1
2
= e 2
e
=
e
e
.
l1 x
z n
l1
xn
z=
(3.16)
Como se cumplen las condiciones del teorema para que el operador de Sturm-Liouville sea
Hermitiano , sabemos que los polinomios de distinto grado son ortogonales entre si. La
norma cuadrada (hn ) se puede calcular facilmente a partir de la formula de Rodrigues:
n Z
Z
n l1 x 2
c
2
Hn (x) n e 2 dx,
hn =
Hn (x) w(x) dx =
(3.17)
l1
x
que al integrar por partes n veces (diferenciando el polinomio tipo Hermite e integrando
la exponencial, que se anula en los lmites) obtenemos
n Z
l1 x 2
c
Hn(n) (x) e 2 dx.
hn =
(3.18)
l1
(n)
(n)
Como la u
nica potencia n-esima del polinomio es n = cn tenemos que Hn (x) = cn n!
y nos queda una integral sencilla de una gaussiana, quedando la norma como
2 n r
2
c
n! .
(3.19)
hn =
l1
l1
Con todo lo expuesto hasta ahora se tiene que la relacion de recurrencia es
Hn+1 (x) cxHn (x)
c2
nHn1 (x) = 0.
l1
(3.20)
Ecuacion diferencial:
Funcion de peso:
w(x) = e
Valores propios:
Funcion generatriz:
Relacion de ortogonalidad:
n = n
P
=
n=0
Hn (x) n
t
n!
x2
x2
n
Hn (x) = (1)n e 2 d n e 2
dx
x2
Hm (x)Hn (x) e 2 dx = n! 2mn
Formula de Rodrigues:
Z
2
xt t2
x2
2
Relacion de recurrencia:
Polinomios de Hermite fsicos
En fsica se toman unos valores distintos para los dos parametros que tenemos libres
c = 2,
l1 = 2,
lo cual hace que ahora tengamos:
Hn (x) 2x Hn (x) + n Hn (x) = 0
Ecuacion diferencial:
w(x) = ex
Funcion de peso:
Valores propios:
e2xtt
Funcion generatriz:
Formula de Rodrigues:
Relacion de ortogonalidad:
n = 2n
P
Hn (x) n
=
n=0 n! t
n
2
2
Hn (x) = (1)n ex d n ex
dx
2
Hm (x)Hn (x) ex dx = 2n n! mn
Relacion de recurrencia:
3.3.2
(3.21)
(3.23)
(3.24)
(n)
Ahora los j solo dependen del termino j + 1 y, mediante induccion, se llega facilmente
a que en general se tiene
(n)
n + n n!
(n)
k =
,
(3.25)
k + k! l1nk
(n)
n
X
k=0
(n)
k xn
n X
n
n + n!
c
(l1 x)k .
=
l1
k!
k
+
k=0
(3.26)
n
X
X
Ln (x) n X X n + (l1 x)k
n + (l1 x)k ct
ct
=
t =
n!
k!
l
k!
l1
k+
k
+
1
n=0
n=0 k=0
n=0 k=0
p+q
X (l1 x)p X q + p + ct p+q
X (l1 x)p X q + p +
ct
=
=
p!
l
p! q=0
l1
p
+
q
1
p=0
q=p
p=0
!p
X
X
1
1
(cxt)p
1
cxt
=
ct p++1 =
ct +1
ct
p! (1 l1 )
p!
(1 l1 )
1 l1
p=0
p=0
cxt
exp 1
ct
l1
.
(3.27)
=
(1 ct
)+1
l1
La formula de Rodrigues se puede obtener de la siguiente manera
n X
n
n!
(n + + 1)
c
(l1 x)p
Ln (x) =
l1
(n p)!p! (p + + 1)
p=0
n
n
X
c
n p xl1 (n + + 1) p+
xl1
l e
x e
x
=
l1
(p + + 1)
p 1
p=0
n
np
n
X
d p xl1
c
d
n
xl1
n+
=
x e
e
x
l1
dxp
dxnp
p
p=0
n
c
dn
=
x exl1 n exl1 xn+ ,
l1
dx
(3.28)
(3.29)
lo que nos permite, del mismo modo que en el caso de los polinomios de tipo Hermite,
calcular la norma cuadrada de estos polinomios (de nuevo, como el operador de SturmLioville es Hermitiano, tenemos que los polinomios de distinto orden son ortogonales entre
s).
n Z
2n
Z
(n + + 1)
c
c
(n)
l1 x n+
2
n!
Ln (x) e x
dx =
,
Ln (x) w(x) dx =
hn =
l1
l1
(l1 )+1
0
0
(3.30)
y por lo tanto, la relacion de recurrencia para este tipo de polinomios es
2
c
Ln+1 (x) (l1 x + 2n + + 1)Ln (x) +
n(n + )Ln1 (x) = 0.
(3.31)
l1
Polinomios de Laguerre generalizados
De nuevo llamamos la atencion acerca de que la notacion estandar para los polinomios de
()
Laguerre generalizados es Ln (x) (con el caso = 0 siendo los polinomios de Laguerre a
c = l1 = 1,
Ecuacion diferencial:
()
()
L()
n n (x) + n Ln (x) = 0
Funcion de peso:
w(x) = x ex
Valores propios:
n = n
P
exp( xt
1t )
=
+1
n=0
(1t)
Funcion generatriz:
n
()
Ln (x) = x ex d n (ex xn+ )
dx
Formula de Rodrigues:
Relacion de ortogonalidad:
Relacion de recurrencia:
3.3.3
()
Ln (x) n
t
n!
()
Ln+1 (x)
()
x
L()
dx = n!(n + + 1)mn
m (x)Ln (x)x e
()
()
Por u
ltimo, tenemos el caso en el que
Q(x) = q0 + q1 x + q2 x2 ,
(3.32)
(3.33)
sin mayor dificultad, de modo que tenemos un intervalo de validez para los polinomios de
tipo Jacobi x (1, 1), donde
R(x) = A
(1 + x)
(1 x)
l0 l1
2q
l0 +l1
2q
w(x) = A
(1 + x)
(1 x)
l0 l1
1
2q
l0 +l1
+1
2q
(3.34)
Como siempre vamos a fijar A = 1 y ademas haremos una redefinicion de las constantes,
escribiendo l0 /q = y l1 /q = ( + + 2) obtenemos
R(x) = (1 + x)+1 (1 x)+1 ,
w(x) = (1 + x) (1 x) .
(3.35)
La ecuacion en su forma mas general para los polinomios de tipo Jacobi queda como
q(1 x2 ) Jn (x) + [ ( + + 2)x]Jn (x) + n Jn (x) = 0.
(3.36)
Se puede ver que R(1) = R(1) con lo que se cumplen las condiciones para que los
polinomios Jn (x) sean ortogonales con valor propio real (siempre y cuando , > 1),
siendo los valores propios
n = n (q(n 1) + + + 2) .
(3.37)
que son los denominados polinomios de tipo Gegenbauer (Cn (x)), esto lo hacemos ya que
todos los polinomios que vamos a nombrar son casos particulares de este. Como en el
caso de los polinomios tipo Laguerre, obviaremos el ndice
en la notacion (y de hecho
no volveremos a escribirlo barrado) durante toda la exposicion. La ecuacion diferencial
para este tipo de polinomios es
(1 x2 ) Cn (x) (2 + 1)x Cn (x) + n Cn (x) = 0.
(3.38)
n = n (n + 2) .
(3.39)
(n)
n2p
(n)
Y
n!
n
=
(n 2p)! k=1 (n 2k)(n 2k + 2) n(n + 2)
(n)
(1 n )
n! n
,
= p
4 p!(n 2p)! (p + 1 n )
(3.40)
de donde
n
X
[n/2]
(n)
n! n
(1 n )
xn2p
p p!(n 2p)! (p + 1 n )
4
p=0
k=0
1
1n n
(n) n
= n x F
, , 1 n , 2 ,
2
2
x
Cn (x) =
(n)
k xk
(3.41)
X
Cn
n=0
n!
t =
=
[n/2]
X
tn X
n=0
n!
n,p=0
p=0
(ct)n
n! cn (n + )
(1 n )
xn2p
2p
2 p!(n 2p)! () (p + 1 n )
(n + )
(1 n )
xn2p
2p)! () (p + 1 n )
22p p!(n
(L + p + ) (1 L p ) Lp
x
p)! () (+1 L )
L,p=0
p
L
X
X
(1 )
L!
ct
ct
(2x)Lp
=
2
(1 + L + ) L! p=0 p!(L p)! 2
L=0
( + L + p) (1 L p)
Lp
(1)
(1 ) ()
L
L
X
ct
(1 )
ct
=
2x (1)L
2
(1
+
L
+
)
L!
2
L=0
=
(ct)L+p
1
(1 cxt +
22p p!(L
c2 2
t)
4
(3.42)
En esta serie de igualdades hemos realizado los siguientes pasos: en la primera lnea
simplemente escribimos el polinomio de tipo Gegenbauer como una serie, en la segunda
extendimos la suma en p hasta infinito, lo cual esta permitido gracias a al factor (n2p)! =
(1 + n 2p) en el denominador, en la tercera lnea se realizo el cambio n = L p y en
la cuarta se reordenaron las series, siendo el termino entre corchetes igual a 1 debido a
propiedades de las funciones gamma, en la pen
ultima lnea resumamos la serie en p (que
no es mas que un binomio de Newton) y escribimos de forma explcita un factor 1 que
nos hace ver claramente que nuestra serie final es otro binomio.
Ahora daremos la f
ormula de Rodrigues. No la vamos a deducir porque es un calculo
con expresiones largas. En muchos textos se asigna como un ejercicio para el deducirla
para alg
un caso particular sencillo, como puede ser el caso = 1/2. La formula es
c n + n+1 + n+1 1
n
2
1/2 d
2
2
2
Cn (x) =
(x2 1)n+1/2 , (3.43)
(x
1)
2
dxn
() + n + 21
donde x es la funcion suelo (el mayor entero menor o igual a x, es decir, equivalente a
la parte entera [x]) y x la funcion techo (el menor entero mayor o igual a x).
La norma cuadrada se calcula como hemos hecho para los polinomios de tipo Hermite
y Laguerre:
Z 1
hn =
Cn (x)2 w(x) dx
1
n+1 1
2 n
Z
+
2 (n + ) 1
+ n+1
c
2
2
=
(1 x2 )n+1/2 dx
2
()
() + n + 12
1
n+1 1
n+1
2 n
+ 2
+ 2 2
c
,
(3.44)
=
2
(n + )()2
y la relacion de recurrencia tiene los coeficientes
an = c(n + ),
cn1 =
bn = 0,
hn
n+
,
hn1 n + 1
(3.45)
Polinomios de Legendre
Tomando valores particulares de tendremos distintos tipos de polinomios que aparecen
de forma bastante com
un a la hora de resolver problemas fsicos (el valor de c es solo una
cuestion de convenio). De este modo, tomando los valores
c = 2,
= 12 ,
definimos los polinomios de Legendre:
Ecuacion diferencial:
Funcion de peso:
w(x) = 1
Valores propios:
n = n(n + 1)
X
1
Pn (x) n
=
t
(1 2xt + t2 )1/2
n!
n=0
Funcion generatriz:
Formula de Rodrigues:
Relacion de ortogonalidad:
Pn (x) =
Z
d n (x2 1)n
dxn
Relacion de recurrencia:
1
2n
2 n!
mn
2n + 1
Ecuacion diferencial:
Funcion de peso:
w(x) = (1 x2 )1/2
Valores propios:
n = n(n + 2)
X
1
Un (x) n
=
t
2
(1 2xt + t ) n=0 n!
Funcion generatriz:
Formula de Rodrigues:
Relacion de ortogonalidad:
Relacion de recurrencia:
Un (x) =
Z
(n+1)! (1)n dn
(1
(2n+1)!! 1x2 dxn
x2 )n+1/2
1
()
dx = n!mn
Um
(x)Un() (x)
2
2
1x
1