Ecuación Diferencial de Bessel y El Método de Frobenius

Ecuacin diferencial de Bessel y el mtodo de
Frobenius
Prof. Dr. Francisco Jos Snchez Sesma

Ayte. Julio Csar Pardo Daino
Matemticas Aplicadas
15 de octubre de 2015
Resumen
Se obtendr la funcin de Bessel de primera especie y de orden p a
partir de la ecuacin diferencial de Bessel de orden p, por medio del
mtodo de Frobenius.
Introduccin
Una clase amplia de problemas en fsico-matemticas, conduce a ecuaciones diferenciales lineales de segundo orden y de coeficientes variables.
dy
d2 y
+ Q( x ) + R( x )y = 0
2
dx
dx
Ejemplo de estas ecuaciones son la ecuacin de Airy, Bessel, Chebyshev,
Hermite, Laguerre y de Legendre por mencionar algunas. Es por ello que
se han desarrollado diversos mtodos para darles solucin. Uno de los
mtodos ms famosos es la solucin de este tipo de ecuaciones por medio
de series de potencia. Esto es que exista una solucin de la forma:
P( x )
y( x ) =
a n ( x x0 ) n
n =0
Debemos considerar tres aspectos, el primero es demostrar que la ecuacin diferencial estudiada admite una solucin como la mostrada anteriormente, el segundo es determinar formalmente los coeficientes an , de tal
manera que y( x ) satisfaga la ecuacin diferencial y por ltimo determinar
si la serie determinada converge y, si es as, para que valores de x x0 .
1
Algunas definiciones y teoremas

Para poder establecer el mtodo de solucin de la ecuacin diferencial
de Bessel de orden p se deben dar un par de definiciones, as como de
teoremas para tener un esquema matemtico que permita desarrollar la
solucin de dicha ecuacin. Para empezar consideremos nuevamente la
ecuacin diferencial lineal de segundo orden y de coeficientes variables:
d2 y
dy
+ a1 ( x ) + a2 ( x ) y = 0
2
dx
dx
Dividamos la ecuacin (1) entre a0 ( x ) obteniendo as:
a0 ( x )
dy
d2 y
+ P( x ) + Q( x )y = 0
dx2
dx
Donde P( x ) =
a1 ( x )
a0 ( x )
y Q( x ) =
(1)
(2)
a2 ( x )
a0 ( x )
Definicin 1. Considere la ecuacin (2). Se dice que x0 es un punto ordinario,

si los coeficientes P( x ) y Q( x ) son analticos en x0 . Si uno de ellos o ambos
no son analticos en x0 se dice que ese punto es singular.
Por analiticidad de los coeficientes P( x ) y Q( x ) en x0 entiendase que
existe una serie de potencias centrada en x0 que converge a P( x ) y Q( x )
respectivamente para toda x de algn intervalo abierto que incluya a x0 .
Definicin 2. Considere la ecuacin (2) y suponga que por lo menos uno de
los coeficientes P( x ) y Q( x ) de la ecuacin no es analtico en x0 , de modo x0
es un punto singular. Si las funciones definidas por
( x x0 ) P ( x ) y ( x x0 )2 Q ( x )
son ambas analticas en x0 , entonces el punto x0 se llama punto singular
regular; si una o ambas no son analticas en x0 entonces al punto x0 se le
llama punto singular irregular.
Ahora que tenemos este par de definiciones, enunciaremos un teorema
sin demostracin que ser la base para desarrollar nuestro mtodo.
Teorema 3. Considere la ecuacin (2) y considere que existe un punto singular x0
en la ecuacin diferencial. Entonces, la ecuacin diferencial tiene por lo menos una
solucin no trivial de la forma
( x x0 )r
Cn (x x0 )n
(3)
n =0
donde r es una constante definida que se puede determinar, esta solucin es valida
alrededor de x0 .
2
Ahora que ya tenemos la seguridad de que por lo menos existe una

solucin y es de la forma (3) alrededor de un punto singular regular x0 de la
ecuacin diferencial, debemos preguntarnos cmo debemos proceder para
encontrar el nmero r y los coeficientes Cn de esta soucin? El procedimiento
fue desarrollado por el matemtico alemn Frobenius (1849-1917), la teora
que justifica el mtodo es algo complicada, por lo cul no ahondaremos en
ella. A continuacin se describen los pasos a seguir para la utilizacin del
mtodo de Frobenius.
1. Sea x0 un punto singular regular de la ecuacin (2) entonces, existe
una solucin y es de la forma:
y( x ) =
Cn (x x0 )n+r
(4)
n =0
con C0 6= 0
2. Suponiendo que la derivacin es vlida trmino a trmino se tiene
dy
= (n + r )Cn ( x x0 )n+r1
dx
n =0
(5)
d2 y
=
(n + r )(n + r 1)Cn ( x x0 )n+r2
dx2
n =0
(6)
Teniendo esto, se sustituyen las series(4), (5) y (6) en la ecuacin (1).

3. Se procede a simplificar la expresin resultante a fin de que tome la
forma
K0 ( x x 0 ) r + k + K1 ( x x 0 ) r + k +1 + K2 ( x x 0 ) r + k +2 + . . . = 0
(7)
donde k es un cierto nmero entero y los coeficientes Ki son funciones

de r y ciertos coeficientes Cn .
4. A fin de que (7) sea valida para toda x se debe hacer:
K0 = K1 = K2 = . . . = 0
(8)
5. Al igualar a cero el coeficiente K0 de la mnima potencia, se obtiene una

ecuacin cuadrtica en r, llamada ecuacin indicial. Las dos raices de
esta ecuacin cuadrtica son los nicos valores posibles de la constante
r de la solucin supuesta.
3
6. Ahora, se igualan a cero los coeficientes restantes K1 , K2 , . . .. Esto conduce a un conjunto de condiciones que contienen la constante r, que
deben ser satisfechas por los diferentes coeficientes Cn de la serie.
7. De las dos raices obtenidas r1 y r2 , considerando que Re(r1 ) Re(r2 ),
se sutituye la raiz r1 en las condiciones obtenidas en el paso 6,y se
escogen los Cn para satisfacer estas condiciones. Ntese que si r1 y r2
son reales y desiguales, entonces r1 es la raiz ms grande.
8. Si r2 6= r1 , se puede repetir el procedimiento del paso 7, utilizando la
raiz r2 . De esta manera, se puede obtener una segunda solucin de la
forma deseada; ntese que si r1 y r2 son reales y desiguales, entonces
r2 es la raiz ms pequea. Sin embargo, la segunda solucin obtenida
en este paso, no puede ser linealmente independiente de la solucin
que se obtuvo del paso 7. Tambin en caso de que r1 y r2 sean reales e
iguales, a solucin obtenida en este paso, es claramente idntica a la
que se obtuvo en el paso 7.
Ecuacin diferencial de Bessel de orden p

Consideremos ahora la ecuacin de Bessel de orden p, donde p < y
adems p > 0:
dy
d2 y
+ x + ( x 2 p2 ) y = 0
(9)
2
dx
dx
Esta ecuacin es una de las ms importantes en la fsico-matemtica,
hace aparicin en los problemas de tipo cilindricos, ejemplo de ello es la
propagacin de ondas o la propagacin de calor en una placa entre otros.
Reescribamos la ecuacin de Bessel:
x2
d2 y 1 dy x2 p2
+
+
y=0
(10)
dx2
x dx
x2
Observamos que x0 = 0 es un punto singular y adems regular de la
ecuacin (10), entonces por el Teorema 3 tenemos que existe una solucin
no trivial de la forma:
y( x ) =
Cn xn+r
(11)
n =0
Entonces podemos sustituir las derivadas de dicha solucin en la ecuacin (10), teniendo as:
n =0
n =0
n =0
(n + r)(n + r 1)Cn xn+r + (n + r)Cn xn+r + (x2 p2 ) Cn xn+r = 0

(12)
Dicha ecuacin puede ser reescrita cmo:
n =0
n =0
{(r + n)2 p2 }Cn xn+r + Cn xr+n+2 = 0
(13)
Sacando los dos primeros trminos de la primer suma tenemos:
(r2 p2 )C0 xr + {(r + 1)2 p2 }C1 xr+1 +

{(r + n)2 p2 }Cn + Cn2 x n+r = 0
n =2
(14)
Del primer coeficiente K0 al igualarlo a 0 obtenemos la ecuacin indicial
representada por r2 p2 = 0, ya que como sabemos C0 6= 0. De aqu que las
raices de r sean r1 = p y r2 = p. Como se estableci que p > 0, entonces
de aqu tenemos que r1 > r2 , entonces se usar r1 para el anlisis. Igualando
el coeficiente K1 = 0 obtenemos:
{(r + 1)2 p2 }C1 = (2p + 1)C1 = 0
(15)
De la ecuacn (15) vemos que que forzosamente C1 = 0. Ahora igualemos a 0 el coeficiente n-simo:
{(r + n)2 p2 }Cn + Cn2 = 0 n 2
(16)
Esta ecuacin la podemos reescribir cmo:

Cn =
Cn2
n(n + 2p)
(17)
La ecuacin (17) es la relacin de recurrencia para las constantes Cn con

n = 0, 1, 2, . . .. De ellas podemos ver claramente que para el caso de n impar
las constantes son igual a cero (C2n+1 = 0). Para el caso de las constantes
pares, es fcil demostrar que se cumple la siguiente relacin de recurrencia:
C2n =
(1)n C0
[2 4 6 . . . 2n] [(2 + 2p)(4 + 2p)(6 + 2p) . . . (2n + 2p)]
(18)
Consideremos que p Z, podemos reescribir la ecuacin (18) en trminos del factorial, teniendo as:
C2n =
(1)n C0 p!
+ p)!
22n n!(n
5
(19)
Ahora sustituimos la ecuacin (19) en la solucin (11) y tenemos:
y( x ) =
(1)n C0 p! 2n+ p
x
2n
n=0 2 n! ( n + p ) !
(20)
La cual podemos expresar como:

y( x ) = C0 p!2 p
(1)n x 2n+ p
n=0 n! ( n + p ) ! 2
(21)
1
Ahora considerando que C0 = p!2
p , obtenemos la llamada funcin de
Bessel de primera clase de orden p entero, expresada cmo:
y( x ) = J p ( x ) =
(1)n x 2n+ p
n=0 n! ( n + p ) ! 2
(22)
Hasta ahora slo se consider que p Z, ahora veamos cuando p <,

en este caso, extendemos el uso del factorial de nmeros enteros a nmeros
reales. Esto se logra con la funcin Gamma, cuya propiedad para n Z es
(n + 1) = n!; tenemos as finalmente la funcin de Bessel de primera clase
de orden p:
Jp (x) =
x 2n+ p
(1)n
n =0 ( n + 1 ) ( n + p + 1 ) 2
(23)
Falta verificar el caso para r2 = p, pero el procedimiento es anlogo; la

solucin completa vendra dada por:
y( x ) = AJ p ( x ) + BJ p ( x )
p 6 Z
(24)
Como vemos este resultado limita el uso a los p que no sean enteros,
donde la funcin J p y J p son linealmente independientes, es por ello que con
otros procedimientos se obtiene otra funcin que garantice la independencia
lineal de las funciones para los p <, esta funcin es conocida como funcin
de Bessel de segunda clase de orden p (Yp ). La demostracin de esta funcin
sale de los alcances del presente trabajo.

Ecuación Diferencial de Bessel y El Método de Frobenius

Cargado por

Copyright:

Formatos disponibles

Ecuación Diferencial de Bessel y El Método de Frobenius

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ecuación Diferencial de Bessel y El Método de Frobenius

Cargado por

Copyright:

Formatos disponibles

Ecuacin diferencial de Bessel y el mtodo de

Prof. Dr. Francisco Jos Snchez Sesma

Algunas definiciones y teoremas

Definicin 1. Considere la ecuacin (2). Se dice que x0 es un punto ordinario,

Ahora que ya tenemos la seguridad de que por lo menos existe una

Teniendo esto, se sustituyen las series(4), (5) y (6) en la ecuacin (1).

donde k es un cierto nmero entero y los coeficientes Ki son funciones

5. Al igualar a cero el coeficiente K0 de la mnima potencia, se obtiene una

Ecuacin diferencial de Bessel de orden p

(n + r)(n + r 1)Cn xn+r + (n + r)Cn xn+r + (x2 p2 ) Cn xn+r = 0

{(r + n)2 p2 }Cn xn+r + Cn xr+n+2 = 0

Sacando los dos primeros trminos de la primer suma tenemos:

(r2 p2 )C0 xr + {(r + 1)2 p2 }C1 xr+1 +

{(r + 1)2 p2 }C1 = (2p + 1)C1 = 0

{(r + n)2 p2 }Cn + Cn2 = 0 n 2

Esta ecuacin la podemos reescribir cmo:

La ecuacin (17) es la relacin de recurrencia para las constantes Cn con

Ahora sustituimos la ecuacin (19) en la solucin (11) y tenemos:

La cual podemos expresar como:

Hasta ahora slo se consider que p Z, ahora veamos cuando p <,

Falta verificar el caso para r2 = p, pero el procedimiento es anlogo; la

También podría gustarte