5 - Funciones de Bessel

E cu a ci o n e s D i f e r e n ci a l e s d e O r d en Su p er i o r
Pa r t e V
Funciones de Bessel
Ing. Ramn Abascal
Profesor Titular de Anlisis de Seales y Sistemas y Teora de los Circuitos II en la UTN, Facultad Regional Avellaneda Buenos Aires, Argentina 2006
Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel.
5.3
Funciones de Bessel.
5.1 - Ecuacin diferencial de Bessel: Segn se vio en la Tercera Parte que una ecuacin diferencial de segundo orden, con coeficientes variables responde a la frmula general ( x - a ) 2 . y" + ( x - a) . ( x ) . y' + ( x ) . y = 0
P P
(3.2)
La ecuacin x y" + x y' + ( x - ) y = 0

2
P
2
P
2
P
Ecuacin de Bessel
(5.1)
que se conoce como Ecuacin Diferencial de Bessel, y en la que es un nmero real no negativo, constituye un caso particular entre las ecuaciones diferenciales de segundo orden. Comparando entre si las ecuaciones (3.2) y (5.1), surgen a primera vista las siguientes relaciones: a = 0 (x) = 1, y ( x ) = x 2 - 2
P P P P
Solucin de la Ecuacin de Bessel: Cualquier funcin y ( x ) que satisfaga la ecuacin (5.1) ser una solucin particular de la Ecuacin de Bessel. Adems, si y1 ( x ) e y2 ( x ) son ambas soluciones de aquella, entonces tambin lo son A y1 ( x ) y B y2 ( x ), donde A y B son dos constantes arbitrarias:
B B B B B B B B
x 2 A y1 " ( x ) + x A y 1 ' ( x ) + A ( x 2 - 2 ) y 1 ( x ) = 0
P P B B B B P P P P B B
x 2 B y2 " ( x ) + x B y 2 ' ( x ) + B ( x 2 - 2 ) y 2 ( x ) = 0
Sumando miembro a miembro estas dos igualdades, obtenemos: x2 ( A y 1 " + B y 2 " ) + x ( A y 1 ' + B y2 ' ) + ( x 2 - 2 ) ( A y 1 + B y2 ) = 0
P P B B B B B B B B P P P P B B B B
(5.2) (5.3)
Si llamamos: y = A y1 + B y2
B B B B
entonces: y' = A y1' + B y2'

B B B B
y" = A y1" + B y2"

B B B B
Lo que demuestra que la (5.3) es tambin solucin de la ecuacin de Bessel. Vamos a tratar ahora de hallar una solucin particular de la ecuacin, es decir, vamos a definir cmo debe ser la funcin y ( x ) para que sea solucin de dicha ecuacin.
R. Abascal - Anlisis de Seales y Sistemas
Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel
5.4
A tal fin, comenzaremos por sustituir en la misma la solucin general (3.4), es decir: y(x) =
k=0
ck ( x - a ) k+
B B P P
Por razones de orden prctico, hemos cambiado el nombre del ndice m por k. Hecho lo cual, reemplazamos el valor de y ( x ) en la ecuacin de Bessel, con lo que obtenemos la ecuacin:
(k + ) (k + -1) ck x
B B P
k+
P
+ (k + ) ck x
B B P
k+
P
k=0
ck x
B B P
k++2
P
-
P
2
P
k=0
ck xk+ = 0
B B P P P
k=0
k=0
(5.4) Esta ecuacin debe cumplirse para todo x, lo que a su vez exige que sean nulos todos los coeficientes que multiplican a cada una de las potencias de x. Esta conclusin nos conduce a encontrar ciertas ecuaciones a partir de las cuales podremos determinar el valor de tales coeficientes. Empecemos por el caso en que k = 0. El coeficiente de x est dado precisamente por la suma de los coeficientes que multiplican a x en cada una de las sumatorias, es decir:
P P P P
( - 1 ) c o + c o - 2 c o = 0
B B B B P P B B
Ecuacin Indicial
Esta igualdad se conoce como Ecuacin Indicial, porque podemos a partir de ella calcular los coeficientes ck (Ver apartado 3.1).
B B
Para resolver la Ecuacin Indicial es necesario que se cumpla que co 0, pues de lo contrario como puede verse por simple inspeccin de la ecuacin, si co = 0, nos encontraramos en un callejn sin salida. Simplificando, hallamos que:
B B B B B B
2 2 = 0
P P P P
Las soluciones de esta ecuacin algebraica de segundo grado son: 1 = +

B B
2 = -
B B P P
De modo similar, si k = 1, la suma de los coeficientes de x1+ , que tambin debe ser nula, es: ( + 1 ) c1 + ( + 1 ) c1 - 2 c1 = ( 2 + 2 + 1 - 2 ) c1 = 0
B B B B P P B B P P P P B B
Reemplazando en ella la primera de las dos soluciones posibles: = 1 = +

B B
tenemos: ( 2 + 2 + 1 - 2 ) c 1 = ( 2 + 1 ) c 1 = 0
P P P P B B B B
Como establecimos que es un nmero no negativo, entonces c1 deber ser necesariamente nulo.
B B
5.5
Nota: Incluso si hubiramos aceptado la segunda solucin: = 2 = -

B B
habramos obtenido: ( 2 - 2 + 1 - 2 ) c 1 = ( - 2 + 1 ) c 1 = 0
P P P P B B B B
por lo que c1 debera ser tambin igual a cero, salvo en el caso puntual en que fuera igual a 0,5(1). Es decir que para la generalidad de los casos se cumple la condicin
B B TP PT
c1 = 0
B B
Finalmente, para k = 2, 3, 4, . . . , obtenemos la forma general: ( k + ) ( k + - 1 ) ck + ( k + ) ck + ck -2 - 2 ck = 0

B B B B B B P P B B
( k + )2 ck - ( k + ) ck + ( k + ) ck - 2 ck + ck -2 = 0
P P B B B B B B P P B B B B
Al simplificarla, resulta [ ( k + )2 - 2 ] ck + ck -2 = 0
P P P P B B B B
Despejando ck obtenemos la siguiente frmula de recurrencia que permite obtener todos los coeficientes, pares e impares, a partir del conocimiento de los dos primeros:
B B
ck =
B B
- ck -2
B B P
(k + )2 - 2
P P
De aqu, como c1 es igual a cero se deduce que todos los coeficientes de ndice impar son nulos:
B B
c1 = c3 = c5 = c7 = . . . = c2s+1 = 0
B B B B B B B B B B
Para determinar los coeficientes de ndice par, sabiendo que = , empezaremos por reemplazar este valor en la frmula de recurrencia: ck =
B B
- ck -2
B B
=
P P P P P P
- ck -2
B B
=
P P P P
- ck -2
B B
(k + )2 - 2
k 2 + 2 k + 2 - 2
k(k+ 2)
P
Como k es un nmero par, podemos hacer: k = 2 m, con m = 0 , 1 , 2, 3, . . .
con lo cual, las frmulas anteriores se modifican como sigue: c2m =

B B
- c2 ( m - 1)
B B
- c2 ( m - 1)
B B
2m(2m+2)
4m(m+)
P
(1)
TP PT
Por ahora dejamos de lado el caso particular en el que = 0,5, que veremos oportunamente. R. Abascal - Anlisis de Seales y Sistemas
5.6
Esta frmula nos permitir calcular todos los coeficientes de ndice 2m, a partir del conocimiento de co, pero como ste no est definido, se acostumbra adoptar para el mismo el siguiente valor, que conduce a aplicaciones interesantes de las funciones de Bessel:
B B
co =
B B
1 2 !
P P
=
P
1 2 ( + 1)
P
A partir de aqu, apelando a la frmula de recurrencia, podemos obtener los dems coeficientes: Para c2 ( En cuyo caso, m = 1 ),
B B
c2 =
B B
- co
B B
- 1 22 ( + 1 ) 2 ( + 1 )
P P P P
4(1 + )
Si recordamos que (n+1) = n. (n ) la expresin del coeficiente c2 se puede reformular como sigue:
B B
c2 =
B B
- 1 22+ ( + 1 ) ( + 1 )
P P B B
=
P P
- 1 22+ ( + 2 )
De modo similar, para c4 ( m = 2 ) ser: c4 =

B B
- c2
B B
=
P P
1 8 . 22+ ( + 2 ) ( + 2 )
=
P P
1 24+ 2! ( + 3 )
8(2 + )
Por reiteracin, esta ltima frmula conduce a la ecuacin siguiente, que permite calcular en forma directa cualquier coeficiente de ndice par: c2m =
B B
(-1)m
P P
(5.5)
2
P
2m+
P
m! ( + m + 1 )
5.2 - Funciones de Bessel: Aplicando para los coeficientes de la ecuacin (3.6) el valor dado por la ecuacin anterior, y recordando que hemos elegido = , se obtiene una solucin particular de la Ecuacin de Bessel, que se conoce como Funcin de Bessel de primera clase y que se designa como J ( x ):
B B
J ( x ) = x
B B P
m=0
( - 1 ) m x2m
P P P
2
P
2m+
P
m! ( + m + 1 )
Funcin de Bessel de 1 Clase

P
(5.6)
El exponente 2m en el numerador tiene en cuenta que los coeficientes de ndice impar son nulos.
5.7
Si en lugar de = hubiramos tomado la segunda solucin de la ecuacin indicial, es decir = , habramos llegado a la solucin:
P P
J ( x ) = x
B B
m=0
( - 1 ) m x2m
P P P P
22m- m! ( m + 1 )
P P
Si no es un entero, cualquier combinacin lineal de estas dos ltimas ecuaciones es una Solucin General de la ecuacin diferencial de Bessel. Por ejemplo: y ( x ) = A J ( x ) + B J ( x )
B B B B
donde A y B son nmeros complejos arbitrarios. En efecto, esta ecuacin no es otra que la (5.3), en la que y1 = J ( x )
B B B B
y2 = J ( x )
B B B B B B B B
Por el contrario, si es un entero, = n, se demuestra que Jn ( x ) y J-n ( x ) son linealmente dependientes entre s, y la solucin anterior no es vlida. En efecto, las ecuaciones son ahora:
Jn ( x ) = x n
B B P P
m=0
P
( - 1 )m x2m
P P P P
(5.7)
22m+n m! ( n + m ) !
P
J-n ( x ) = x -n
B B P P
m= n
( - 1 )m x2m
P P P P
(5.8)
22m-n m! ( m - n ) !
P P
En la segunda de ellas, los ndices comienzan a partir de m = n, porque, como el factorial de un nmero entero negativo es infinito, todos los trminos para los que m es menor que n son nulos, al estar divididos por una cantidad infinita. Ahora demostraremos la dependencia lineal entre ambas soluciones. Empezaremos llamando: m = n+s Por lo tanto, m - n = s, y 2m - n = 2n + 2s - n = 2s + n
Reemplazando en la (5.8) resulta:

J-n ( x ) =
B B
m=n
P
( - 1 ) m x2m-n
P P P P
s=0
( - 1 ) n+s x2s+n
P P P P P P
22m-n m! ( m - n )!
P
22s+n ( n + s )! s!
Es decir
5.8
J-n ( x ) = ( - 1 )n xn
B B P P P P
s=0
P
( - 1 ) s x 2s
P P P P P
22s+n s! ( n + s )!
P P B B
Si cambiamos el nombre del ndice s por m, esta ecuacin es igual a ( - 1 ) n . Jn ( x ). Queda por consiguiente probado que Jn ( x ) y J-n ( x ) son linealmente dependientes entre s.
B B B B
Como conclusin, cuando es un nmero entero, para encontrar una solucin general de la ecuacin diferencial se hace necesario recurrir a otro tipo de funciones Bessel: Las llamadas funciones de Bessel de segunda clase que estudiaremos en el pargrafo 5.9. A continuacin representamos las funciones Jo ( x ) y J1 ( x ):
B B B B
1,2 1 0,8 0,6 0,4 0,2 0 -0,2 -0,4 -0,6 Jo (x) J1(x)
5.3 - Expresion integral de las funciones de Bessel: Hasta ahora hemos expresado las funciones de Bessel a travs de desarrollos en series de potencias de x. Sin embargo, como veremos en este apartado, dichas funciones fueron introducidas histricamente como integrales que involucran relaciones entre funciones trigonometricas de x, de la forma siguiente: Sea una funcin f(z) analtica en una corona con centro en el orgen de coordenadas, a la que desarrollaremos en serie alrededor del punto z0 = (0, 0).
B B
Empezaremos por recordar el desarrollo en serie de Laurent de una funcin f ( z ), en zo = (0,0):

B B
f(z) =
n=0
an z
B B P
n
P
n=1
bn z -n
B B P P
5.9
Los coeficientes an y bn estn dados por las frmulas:

B B B B
an =
B B
1 2i
f(z) z
P
dz
n+1
P
bn =
B B
1 2i
f(z) z
P
dz
-n+1
P
Para nuestro caso vamos a llamar: an = n

B B B B
bn = -n
B B B B
As, la serie de Laurent que define la funcin f ( z ), puede escribirse:

f(z) =
n=0
n z n +
P B B P
n=1
-n z -n =
B B P P
n z n
P B B P
n=-
Donde
n =
B B
1 2i
f ( z ) dz C z n+1
P P
Demostraremos seguidamente que si la funcin f ( z ) es f() = e siendo = ei

P P B B
x ( 1) 2
B P P
es un complejo de mdulo unitario, entonces n es igual a la funcin de Bessel de primera clase y orden n. Antes de continuar, trataremos de averigar qu forma adopta n en este caso. Derivando respecto de ,
B B
d = i e i d
P P
y reemplazando en la integral, obtenemos:

B
n =
B B
1 2i
C
2
x ( e i e i) 2
P P
e in . e i
P P P P
i e i d =
P P
1 2
e ix sen d 0 e in
P P P P
2
B
y finalmente:
B B
n =
B B
1 2
e i (x sen
P
n)
P
Funcin de Bessel de primera clase y orden n
(5.9)
5.10
En la Seccin 5.4 trataremos de relacionar esta forma con la definida en el pargrafo 5.2. De hecho, demostraremos que n = Jn ( x )
B B B B
Por otro lado, como e i (x sen

P
n)
P
= cos ( x sen - n ) + i sen ( x sen - n )
la funcin de Bessel se puede escribir tambin de la manera siguiente:

B
2
B
n =
B B
1 2
cos ( x sen - n ) d +
sen ( x sen - n ) d
(5.10)
2 0
En la ltima integral el integrando es el seno de una diferencia. Por tanto, dicha integral puede desdoblarse as:
B
2
B
sen ( x sen - n ) d =
2
B
0
B B
=
B B
0
B
sen ( x sen ) cos n d
0
B B
cos ( x sen ) sen n d
Llamemos: f1 ( ) = sen ( x sen )

B
y
2
B B
f2 ( ) = cos ( x sen )
2
Entonces:
B
2
B
sen ( x sen - n ) d = f1 ( ) cos n d

0
B B
f2 ( ) sen n d
B B
(5.11)
0
B
Tratndose de un seno, f1 ( ) es una funcin impar, lo mismo que f2 ( ), por tratarse de un coseno, es una funcin par.
B B
Si comparamos la frmula de los coeficientes an en el desarrollo de una funcin en serie de Fourier,

B B B
2
B
an =
B B B B
1
B
f1 ( t ) cos n t d t
B B
con la primer integral del segundo miembro de la (5.11), vemos que, si an es el coeficiente del desarrollo de Fourier de f1 ( ), la misma es igual a an :
B B B B B B
2
B
f1 ( ) cos n d = an
B B B B
0
B
5.11
Pero siendo f1 ( ) impar, los coeficientes an, y por tanto, la integral, debern ser todos nulos. De modo similar, como f2 ( ) es una funcin par, la segunda integral en la (5.10) es tambin nula:
B B B B B B B
2
B
cos ( x sen ) sen n d =
f2 ( ) sen n d = bn = 0
B B B B
0
B
En resumen, la (5.11) queda reducida a la igualdad:

2
n =
B B
1 2
cos ( x sen n ) d
Expresin Integral de la Funcin de Bessel de 1 Clase
Siendo el integrando una funcin coseno, de perodo 2 , puede ponerse, alternativamente:

2
n =
B B
1 2
cos ( x sen n ) d =
cos ( x sen n ) d
(5.12)
5.4 - Demostracin de la equivalencia entre las formas polinmica e integral: A continuacin, a partir de los desarrollos en serie de las funciones elementales, probaremos que las formas (5.6) (5.7) de las funciones de Bessel y la expresin (5.9) son equivalentes. La forma exponencial el nmero complejo , de mdulo unitario, es: = 1 . e i
P P
de lo cual se deducen las relaciones siguientes: = cos + i sen , n = cos n + i sen n

P P
-1 = cos - i sen ,
P P
-n = cos n - i sen n
P P
(5.13)
Restando entre s las dos primeras se advierte que: - -1 = 2 i sen

P P
de donde: e
P
i sen
P
= e
P
2
P
. e
P
2
P
Si elevamos ambos miembros a la potencia x, obtenemos: e

P
i x sen
P
x ( - -1 )
B
x
2
P P
-x
e
P
2
PB PB P
= e
. e
P
2
P
Al desarrollar las dos exponenciales en serie, encontramos una nueva igualdad:

5.12
2 2 3 3 e i x sen = 1 + x + 1 x + 1 x + . . . 2 2! 4 3! 8
P P P P P P P P P P
2 1- x + 1 x - 1 2 2! 42 3!
P P P P
x3 + . . . 83
P P P P
Haremos seguidamente el producto de cada uno de los trminos de una serie por todos los de la otra. Agrupando los correspondientes a potencias de de igual orden, resulta: e i x sen = 0
P P P P
2 x4- 1 x6 + ... 1- x + 1 4 (2!)2 16 (3!)2 64

P P P P P P P P P P
x - 1 2 2! 1
P P P P
x3 + 1 23 2! 3!
P P P P P P
x5 - 1 x7 + . . . 25 3! 4! 27
P P P P P P P P P P
+ 2
2! + -1
P P
x2 - 1 4 3!
P P
x4 + 1 x6 - . . . 24 2! 4! 26
P P P P P P P
+ ... +
3 x5 + 1 x7 - . . . - x + 1 x - 1 2 2! 23 2! 3! 25 3! 4! 27
P P P P P P P
+ -2
P P
x2 - 1 2! 4 3! 1
P P
x4 + 1 x6 - . . . 2! 4! 26 24
P P P P P P P P
+ ...
(5.14)
Recordemos la frmula (5.7) de la funcin Bessel de 1 clase:
n
P P
Jn ( x ) = x
B B
m=0
( - 1 )m x2m
P P P P
22m+n m! ( n + m ) !
P P
para valores sucesivos de n, esta funcin vale, respectivamente: Para n = 0: (-1)0 x0 + x 0 (-1) x2 + x 0 (-1)2 x4 + x 0 (-1)3 x6 + . . . Jo ( x ) = x 0 20 0! 0! 22 1! 1! 24 2! 2! 26 3! 3!
P P P P P P P P P P P P P P P P B B P P P P P P P P P P P P P P P P
2 x4 - 1 x6 + . . . Jo ( x ) = 1 - x + 1 4 (2!)2 16 (3!)2 64
P P P P P P B B P P P P
Para n = 1: (-1)0 x0 + x 1 (-1) x2 + x 1 (-1)2 x4 + x 1 (-1)3 x6 + . . . J1 ( x ) = x 1 2 0! 1! 23 1! 2! 25 2! 3! 27 3! 4!

P P P P P P P P P P P P P P P P B B P P P P P P P P P P P P P P P P
J1 ( x ) =
B B
x - 1 2 2!
x3 + 1 23 2! 3!
P P P P
x5 - 1 x7 + . . . 25 3! 4! 27
P P P P P P P P
5.13
De modo similar: J2 ( x ) =
B B
1 2!
x2
P P
1 3!
22
P P
x4 - 1 x6 + . . . 24 2! 4! 26
P P P P P P P P
... J-1 ( x ) = B B
x + 1 2 2! 1 2!
P P
x3 - 1 x5 + 1 x7 - . . . 2! 3! 25 3! 4! 27 23
P P P P P P P P P P P P
J-2 ( x ) =
B B
x2 - 1 3! 22
P P
x4 + 1 x6 - . . . 24 2! 4! 26
P P P P P P P P
... Podemos ver que cada uno de los parntesis en la (5.14) se corresponde exactamente con una de las funciones Jn ( x ). Reemplazando en la mencionada igualdad, se obtiene:
B B
e i x sen = Jo ( x ) + J1 ( x ) + 2 J2 ( x ) + 3 J3 ( x ) + . . . +
P P B B B B P P B B P P B B
+ -1 J-1 ( x ) + -2 J-2 ( x ) + -3 J-3 ( x ) + . . .

P P B B P P B B P P B B
(5.15)
Vimos por otra parte que si n es un entero, entonces: J-n ( x ) = ( - 1 ) n . Jn ( x ),

B B P P B B
Sustituyendo esta igualdad en la (5.15), y agrupando trminos, hallamos: ei x sen = Jo ( x ) + ( - -1 ) J1 ( x ) + ( 2 + -2 ) J2 ( x ) + ( 3 - -3 ) J3 ( x ) + . . . (5.16) Recurriremos ahora a las ecuaciones (5.13),
P P B B P P B B P P P P B B P P P P B B
n = cos n + i sen n
P P
-n = cos n - i sen n
P P
(5.13)
Sumndolas o restndolas entre s, deducimos respectivamente: n + -n = 2 cos n

P P P P
n - -n = 2 i sen n
P P P P
Entonces, reemplazando en la (5.16), hallamos: ei x sen = Jo ( x ) + 2 [ J2 ( x ) cos 2 + J4 ( x ) cos 4 + . . . ] +

P P B B B B B B
+ 2 i [ J1 ( x ) sen + J3 ( x ) sen 3 + J5 ( x ) sen 5 + . . . ]

B B B B B B
(5.17) (5.18)
Pero tambin, como ei x sen = cos ( x sen ) + i sen ( x sen )

P P
igualando las partes reales e imaginarias, se obtiene:
5.14
cos ( x sen ) = Jo ( x ) + 2 [ J2 ( x ) cos 2 + J4 ( x ) cos 4 + . . . ]

B B B B B B
(5.19) (5.20)
sen ( x sen ) = 2 [ J1 ( x ) sen + J3 ( x ) sen 3 + J5 ( x ) sen 5 + . . . ]

B B B B B B
Al multiplicar ambos miembros de la primera de estas igualdades por cos n e integrar luego entre 0 y , encontramos:
cos n cos ( x sen ) d = Jo ( x ) cos n d + 2 J2m ( x ) cos 2m cos n d

B B B B
m=1
En el segundo miembro podemos observar que, con la nica excepcin del trmino para el cual n es igual a 2m, todos los dems son nulos, porque estn multiplicados por integrales de funciones coseno(1) entre los lmites 0 y . Por lo tanto, la ecuacin se reduce as:
TP PT
cos n cos ( x sen ) d = 2 Jn ( x )

B B
cos2 n d
P P
donde n es entero y par (n = 2m). Ahora slo queda resolver la ltima integral:

2
P P
cos n d =
1 + cos 2n 2
d = 1 d + 0 2 0
cos 2n 2
d = 2
Por consiguiente Jn ( x ) =
B B
cos n . cos ( x sen ) d
(5.21)
0 De idntica forma, multiplicando la (5.20) por sen n, e integrando, se obtiene, para n impar: 1
Jn ( x ) =
B B
sen n . sen ( x sen ) d
(5.22)
Es posible agrupar estas dos ecuaciones, (5.21) y (5.22) en una sola. En efecto, como la primera de ellas es igual a cero para n impar y la segunda lo es igualmente cuando n es par, podemos decir que, para cualquier valor de n entero y positivo, se satisface la igualdad: Jn ( x ) =
B B
[ cos n . cos ( x sen ) + sen n . sen ( x sen ) ] d
0 O tambin, como el integrando es igual al coseno de la diferencia entre n y x sen :

(1)
TP PT
Recurdese que cos a cos b = 1/2 [( cos a + b ) - cos ( a - b )]
5.15
Jn ( x ) =
B B
cos ( n - x sen ) d
0 Pero esta ltima integral es la misma que la que aparece en la frmula (5.12), con lo que queda demostrado que n = Jn ( x )
B B B B
5.5 - Funcin de Bessel para = 0,5: Un caso particular, interesante por las aplicaciones que se deducen del mismo, es aquel para el cual es igual a 0,5. Trataremos en este apartado de hallar la funcin de Bessel correspondiente. Para ello, reemplazando en la (5.6) obtenemos:
J0,5 ( x ) = x 0,5
B B P P
m=0
P
( - 1 )m x2m
P P P P
22m+0,5 m! ( m + 0,5 + 1)
P
Que se puede escribir tambin: J0,5 ( x ) =

B B
x
P
0,5
P
m=0
P
( - 1 )m x2m
P P P P
22m m! ( m + 1,5 )
P
En el denominador aparece la funcin ( m + 1,5 ), que analizaremos a partir de la igualdad: (n) = (n-1) (n-1) En este caso: J0,5 ( x ) =
B B
(5.23)
x
P
0,5
P
m=0
( - 1 )m x2m
P P P P
(5.24)
2
P
2m
P
m! ( m + 0,5) ( m + 0,5 )
Volviendo a aplicar la misma relacin, ( n ) = ( n - 1 ) . ( n - 1 ), esta vez a ( m + 0,5 ), y en forma reiterada, ( m + 0,5 ) = ( m + 0,5 - 1 ) . ( m + 0,5 - 1 ) = = ( m + 0,5 - 1 ) ( m + 0,5 - 2 ) ( m + 0,5 - 3 ) . . . ( 0,5 ) = ( m - 0,5 ) ( m - 1,5 ) ( m - 2,5 ) . . . ( 0,5 ) Recordemos que el valor de la funcin Gamma de 0,5 (Ver Parte I) es: ( 0,5 ) = Entonces, para m = 0 tenemos:
5.16
1 2 Para m = 1: 3 2 Para m = 2 5
( 0,5 ) =
3 2
. 1 2 2
( 0,5 ) =
3 4
2 2 .......
. 3 . 1 2 2 2
( 0,5 ) =
15 8
Notar que por razones de mayor claridad estamos utilizando indistintamente las formas fraccionaria o decimal, segn nos convenga en cada oportunidad. Reemplazando en la (5.24), y efectuando operaciones, se obtiene: J0,5 ( x ) =
B B
x
P
0,5
P
4 x2
P P
+
P
8x4
P P P
...
2 Si sacamos factor comn , queda: J0,5 ( x ) =

B B
4 . 3 .
16 . 2 . 15 .
x
P
0,5
P
2
0,5
P P
x2 + x4 ... 3 60
P P P P P P
A continuacin sacaremos nuevamente factor comn, esta vez, 2/x, J0,5 ( x ) =

B B
x 2 2
P
2 x
x 3 + x5 . . . 6 120
P P P P P P
0,5
P
J0,5 ( x ) =
B B
x3
P P
+
P
x5
P P
...
3!
5!
El ltimo factor es igual al desarrollo en serie de sen x. Reemplazando, se obtiene finalmente la funcin de Bessel para = 0,5: J0,5 ( x ) =
B B
2
P
0,5
P
sen x
De modo similar se demuestra que el valor de la funcin de Bessel para = - 0,5 es: J0,5 ( x ) =
B B
2
P
0,5
P
cos x
5.17
5.6 - Propiedades de las funciones de Bessel: Estudiaremos en este apartado y los siguientes algunas propiedades de las funciones de Bessel. En primer lugar, demostraremos la igualdad: d dx La demostracin ms sencilla consiste en trabajar con ambos miembros en forma separada, hasta llegar a una misma expresin, lo que nos permitir, por carcter transitivo, confirmar la validez de la ecuacin. Empezaremos por el miembro de la izquierda. Para ello, tomemos la (5.7), y multipliqumosla por x -n , con lo cual obtenemos:
P P
x-n Jn ( x )
P P B B
= - x-n Jn+1 ( x )
P P B B
(5.25)
x
P
-n
P
Jn ( x )
B B
m=0 2
P
( - 1 )m x2m
P P P P
2m+n
P
m! ( n + m )!
que ahora debemos derivar respecto de x:
d dx
x
P
-n
P
Jn ( x )
B B
m=0
( - 1 )m 2 m x2m-1
P P P P P P
(5.26)
2 2m+n m! ( n + m )!
B B
Pasaremos ahora a analizar la funcin de Bessel Jn+1 ( x ), que aparece en el miembro de la derecha en la (5.25), y cuya forma es la siguiente:
Jn+1 ( x ) = xn+1
B B P P
m=0
( - 1 )m x2m
P P P P
2
P
2m+n+1
P
m! ( n + m + 1 )!
Mulitiplicaremos ambos miembros por x-n:

P P
x-n Jn+1 ( x ) = x
P P B B
m=0
P
( - 1 )m x2m
P P P P
2 2m+n+1 m! ( n + m + 1 )!
P
A continuacin vamos a recurrir al artificio de reemplazar: m + 1 = k y por tanto, m = k 1 Entonces:
x
P
-n
P
Jn+1 ( x ) = x .
B B
( - 1 )k-1 x2k-2
P P P P
(5.27)
k=1 2
P
2k-2+n+1
P
( k - 1 )! ( n + k )!
5.18
Al llegar aqu introduciremos en la sumatoria el trmino correspondiente a k = 0, aprovechando que el mismo es nulo, con lo cual no altera el valor de la sumatoria. Efectivamente, como el producto factorial de un nmero entero negativo, en este caso (-1)!, es infinito, entonces: ( - 1 )k-1 x2k-2
P P P P
=
k=0
P
( - 1 ) -1 x -2
P P P P P
= 0
2 2k-2+n+1 ( k - 1 )! ( n + k )!
P P
2 -2+n+1 ( - 1 )! ( n )!
Ello nos permite completar la sumatoria en la (5.27), agregndole el trmino mencionado sin que la misma se altere en absoluto. Es decir:
x-n Jn+1 ( x ) = x
P P B B
k=0
P
( - 1 )k-1 x2k-2
P P P P P
2 2k+n-1 ( k - 1 )! ( n + k )!
Si introducimos la variable x dentro de la sumatoria, y multiplicamos numerador y denominador por 2 k, obtenemos:

x-n Jn+1 ( x ) =
P P B B
k=0 2
( - 1 )k-1 . x . 2 . k . x2k-2
P P P P
( - 1 )k-1 . 2 k . x2k-1
P P P P
2k+n-1
P P
2 k ( k -1 )! ( n + k )!
k=0 2
P
2k+n
P
( k )! ( n + k )!
Apelando tambin a la igualdad: ( - 1 ) k-1 = ( - 1 ) k . ( - 1 ) -1 = - ( - 1 ) k,

P P P P P P P P
llegamos por fin a la ecuacin:
x
P
-n
P
Jn+1 ( x ) =
B B
k=0
( - 1 )k . 2 k . x2k-1
P P P P P P
2 2k+n ( k )! ( n + k )!
cuyo segundo miembro es idntico al de la (5.26), salvo porque en la sumatoria hemos reemplazado el nombre del ndice m por k. Queda por lo tanto demostrada la propiedad (5.25). Ms sencillo resulta demostrar la propiedad siguiente: d dx En efecto, en este caso, a partir de (5.7) tenemos:

xn J n ( x )
P P B B
= xn Jn-1 ( x )
P P B B
(5.28)
x Jn ( x )
P P B B
= x
P
2n
P
m=0
P
( - 1 )m x2m
P P P P
=
P
m=0
P
( - 1 )m x2(m+n)
P P P P
2 2m+n m! ( m + n )!
P
2 2m+n m! ( m + n )!
P
Derivando respecto de x:
d dx
x Jn ( x )
P P B B
m=0
( - 1 )m 2 ( m + n ) x2(m+n)-1
P P P P
2 2m+n m! ( m + n )!
P P
5.19
Dividimos ahora numerador y denominador por 2 y por ( m + n ), y operando, se obtiene en forma sucesiva:
d dx d dx d dx
xn J n ( x )
P P B B
m=0
P
( - 1 )m x2(m+n)-1
P P P P P
2 2m+n-1 m! ( m + n - 1 )! ( - 1 )m x2m+(n-1)
P P P P
x Jn ( x )
P P B B
= x
P
n
P
m=0
P
2 2m+(n-1) m! ( m + n - 1 )!
P
x Jn ( x )
P P B B
= x .x
P P P
n-1
P
m=0
( - 1 ) m x2m
P P P P
= xn Jn-1 ( x )
P P B B
2
P
2m+(n-1)
P
m! ( m + n - 1 )!
Con lo cual queda demostrada la (5.28).
5.7 - Propiedades de las races. En este apartado veremos algunas propiedades relacionadas con las races de las funciones de Bessel. 1 - Alternancia de las races de Jn-1 ( x ) y Jn ( x ).
B B B B
De las ecuaciones (5.25) y (5.28): d dx y d dx surge que entre cada dos ceros de Jn ( x ) cae un cero de Jn+1 ( x ) y a la inversa. Lo mismo vale para las funciones Jn-1 ( x ) y Jn ( x ). Y por extensin a todas las funciones de la familia.
B B B B B B B B
xn J n ( x )
P P B B
= xn Jn-1 ( x )
P P B B
x-n Jn ( x )
P P B B
= - x-n Jn+1 ( x )
P P B B
Notemos en primer lugar que la funcin xn Jn ( x ) tiene los mismos ceros que Jn ( x ), ms otros n ceros ubicados en el origen de coordenadas.
P P B B B B
Jn ( x )
B B
x-n Jn ( x )
P P B B
5.20
Por su parte, la funcin x-n Jn ( x ) tendr tambin los mismos ceros que Jn ( x ) , ms otros n ceros ubicados en el infinito.
P P B B B B
Es decir que si y son dos ceros consecutivos, no nulos ni impropios, de Jn ( x ), tambin lo son de xn Jn ( x ) y de x-n Jn ( x ).
B B P P B B P P B B
Adems, entre y debe existir un valor para el cual xn Jn ( x ) presente un mximo o un mnimo. Dicho valor es en consecuencia un cero de la funcin
P P B B
d dx
xn J n ( x )
P P B B
= xn Jn-1 ( x )
P P B B
y por lo tanto, por un razonamiento similar, tambin lo es de Jn-1 ( x ).

B B
Es decir:
Jn-1 ( ) = 0
B B
( 0)
Veamos esto en forma grfica: Jn1 ( x )

B B
Jn ( x )
B B
Razonando en forma inversa que antes, podemos decir que si es un cero no nulo de xn Jn-1 ( x )
P P B B
n Jn-1 ( ) = 0
P P B B
( 0)
B B P P
es necesariamente tambin un cero de Jn-1 ( x ). Esto es lgico, puesto que los n ceros de xn estn todos en el origen. De modo similar, entre y debe existir un valor, , para el cual x-n Jn ( x ) presente un mximo o un mnimo. Dicho valor es en consecuencia un cero de la funcin
P P B B
d dx
x -n Jn ( x )
P P B B
= - x -n Jn+1 ( x )
P P B B
Por lo tanto, si es un cero no infinito de x-n Jn+1 ( x ), es decir

P P B B
-n Jn+1 ( ) = 0
P P B B
( )
B B
entonces es necesariamente un cero de Jn+1 ( x ).
5.21
2 - Las funciones de Bessel tienen infinitos ceros. Es posible demostrar que las funciones de Bessel tienen infinitos ceros. O, lo que es igual, que la ecuacin: Jn ( x ) = 0
B B
para todo n, tiene infinitas races. La representacin de las funciones J0 ( x ) y J1 ( x ), al final del pargrafo 5.2, no deja dudas al respecto. Al aplicar a las mismas la propiedad 1, deducimos que tambin J2 ( x ) tendra infinitas races, y por induccin, lo mismo debera ocurrir con todas las funciones de Bessel de primera clase y coeficiente n, entero.
B B B B B B
Volveremos ms adelante sobre este punto. 3 - Ecuaciones diferenciales de segundo orden con coeficientes variables diferentes: Sean las dos ecuaciones diferenciales siguientes: u" ( x ) + g ( x ) . u ( x ) = 0 y v" ( x ) + h ( x ) . v ( x ) = 0 h(x) g(x) (5.29) (5.30)
con la condicin:
Esta condicin es necesaria, puesto que si h y g fueran iguales, estaramos en presencia de una misma ecuacin. Despejando en las (5.29) y (5.30) hallamos: u" ( x ) = - u . g y v" ( x ) = - v . h
Si multiplicamos ambos miembros por v en la primera ecuacin y por u en la segunda, es decir: u" v = - u v . g y v" u = - u v . h
y restamos miembro a miembro, resulta: v u" - u v" = ( h - g ) . u v Si y son dos races consecutivas de la ecuacin u ( x ) = 0, lo que equivale a decir que: u() = 0 y u ( ) = 0, (5.31)
como el segundo miembro de la (5.31) es igual a cero, debe cumplirse necesariamente que
5.22
v ( ) . u"( ) - u ( ) . v"( ) = 0 o tambin: v ( ) . u"( ) - u ( ) . v"( ) = 0 De aqu: v ( ) . u"( ) = u ( ) . v"( ) = 0 y v ( ) . u"( ) = u ( ) . v"( ) = 0 (5.33) Supongamos por un momento que elegimos una funcin v ( x ) que satisface la ecuacin diferencial (5.30), pero tal que ni ni sean races suyas. En tal caso, es evidente que para que se cumplan las relaciones (5.32) y (30,.33), la derivada segunda de la funcin u ( x ) debe ser nula en y en , respectivamente. La explicacin anterior nos fuerza a aceptar que u ( x ) y sus derivadas primera y segunda deben comportarse aproximadamente como muestra la figura siguiente: (5.32)
u'(x) u" (x)
u (x)
Ahora procederemos a integrar la ecuacin (5.31) entre y :

( v u" - u v" ) dx =
( h - g ) . u v . dx
Y calcularemos la primera integral por partes. Para lo cual, teniendo en cuenta que u" dx = y lo mismo: v" dx = dv' la procederemos a modificar previamente as: du' dx = du' dx
5.23
( v du' - u dv' )

( h - g ) . u v . dx
v du'
u dv' =
( h - g ) . u v . dx

v u'
v' u' dx - u v'
u' v' dx
( h - g ) . u v . dx
Simplificando las dos primeras integrales, iguales entre s, resulta
es decir:
( v u' - u v' )
( h - g ) . u v . dx
(5.34)
v ( ) u' ( ) - u ( ) v' ( ) - v ( ) u' ( ) + u ( ) v' ( ) = Pero como u ( ) = 0 y u ( ) = 0, simplificando nuevamente queda: v ( ) u' ( ) - v ( ) u' ( ) =
( h - g ) u v dx
( h - g ) u v dx
(5.35)
Supongamos ahora que la funcin v ( x ), que como dijimos no es nula ni en ni en , se mantiene positiva (negativa) en todo el intervalo entre y . Hagamos aqu tambin un anlisis grfico mostrando la situacin en el intervalo : u (x) v(x) u' ( x )
La figura muestra en primer lugar que u' ( ) tiene distinto signo que u' ( ).
5.24
Tambin observamos en la ecuacin (5.35) que si v ( x ) se mantiene positiva (o negativa) en todo el intervalo , el primer trmino de la ecuacin no puede ser igual a cero en nngn punto del mismo, por lo que la integral no puede ser nula en el intervalo, lo que implica que, necesariamente, debe ser h(x) g(x) Esto es obvio, porque si fueran iguales, la integral sera nula. 4 - Races de Jo ( x ):
B B
Volvamos a la ecuacin de Bessel (5.1): x2 y" + x y' + ( x2 - n2 ) y = 0

P P P P P P
y hagamos en ella el cambio de variable: u = y x y = u x -1/2

P P
(5.36)
La derivada primera de esta funcin es y' = u' x -1/2 P P
1 2
u x -3/2
P P
y la derivada segunda: y" = u" x-1/2 1 2

P P
u' x-3/2 1 u' x-3/2 + 2

P P P P
3 4
u x-5/2
P P
Al reemplazar en la ecuacin diferencial de Bessel obtenemos: x2

P P
u" - u' + 3 u x xx 4 x2 x
P P
+x
u' - u x 2x x
+ ( x 2 - n2 )
P P P P
u x
= 0
Ecuacin que, una vez simplificada, queda as: x2 u" - x u' +

P P
3 u + x u' - u 4 2
+ x2 u - n 2 u = 0
P P P P
Simplificando de nuevo y agrupando los coeficientes de las potencias de igual grado de la variable u, resulta x2 u" + u
P P
x2 - n 2 + 1 4
P P P P
= 0
Dividiendo por x2 y reagrupando trminos, esta ecuacin se puede modificar como sigue:
P P
5.25
u" + u
2 1 - 4n -1 4 x2
P P P
= 0
(5.37)
En el caso en que n = 0, la ecuacin anterior se simplifica ms an, quedando: u" + u Sintetizando: Esta ecuacin es la misma que la (5.29), con g = 1 + Su solucin es: u = y x Adems, como se trata de la ecuacin de Bessel, con n = 0, tambin debe ser y = Jo .
B B
1 +
1 4 x2
P
= 0
(5.38)
1 4 x2
P
(5.39)
Por tanto: u = y x = Jo ( x ) x
B B
(5.40)
Verificaremos a continuacin esta afirmacin. Si en la ecuacin (5.7),
Jn ( x ) = x n
B B P P
m=0
( - 1 )m x2m
P P P P
22m+n m! ( n + m ) !
P P
hacemos n = 0, entonces:
Jo ( x ) =
B B
m=0
P P
( - 1 )m x2m
P P P P
22m ( m! ) 2
P P P P P P
Es decir:
2 x4 - 1 x6 + . . . Jo ( x ) = 1 - x + 1 4 (2!)2 16 (3!)2 64
P P B B P P P P
Reemplazado en la (5.40), nos queda: u = x Jo ( x ) = x1/2 B B P P
1 4
x 5/2 +
P P
1 64
x 9/2 P P
1 2304
x 13/2 + . . .
P P
Derivaremos ahora para obtener u': u' = 1 2 x -1/2 - 5 x 3/2 + 9 x 7/2 - 13 x 11/2 + . . . 8 128 4608
P P P P P P P P
5.26
Y derivando nuevamente, u" = - 1 4

Tambin:
x -3/2 - 15 16
P P P P
x 1/2 +
P P
63 256
x 5/2 P P
143 9216
x 9/2 + . . . . .
P P
u 4x
P
=
2
P
1 x -3/2 - 1 x 1/2 + 1 x 5/2 - 1 x 9/2 + . . . 4 16 256 9216

P P P P P P
Basta reemplazar estos valores en la (5.38), para verificar que efectivamente la (5.40) es una solucin de la misma. En efecto, de tal reemplazo resulta: u" + u + u 4x
P
= 2
P
1 x -3/2 - 15 x 1/2 + 63 x 5/2 - 143 x 9/2 + . . . + 4 16 256 9216

P P P P P P P P P P
+ x1/2 1 4
1 x 5/2 + 1 x 9/2 - 1 x 13/2 + . . . + 4 64 2304

P P P P P P
x -3/2 - 1 x 1/2 + 1 x 5/2 - 1 x 9/2 + . . . 16 256 9216

P P P P P P P P
Y al agrupar en esta ecuacin los coeficientes de las potencias de x de igual grado, hallamos: u" + u + u 4x
P
=
2
P
- 1 4 4
P P
x -3/2 +
P P P
16 - 15 16 16
P
1 16
x 1/2 +
P P P
63 - 64 + 1 256 256 256
143 - 144 + 1 x 9/2 + . . . 9216 9216 9216

P P
Como todos los coeficientes resultan nulos, comprobamos que, efectivamente, la (5.40) es solucin de la (5.38). Volvamos ahora a la ecuacin (5.30) y hagamos en ella h ( x ) = 1. Este valor satisface la condicin h ( x ) g ( x ). En efecto, para cualquier x , g = 1 + 1 4x
P
se verifica que: >

2
h = 1
Al reemplazar en (5.30) el valor adoptado para h, la misma queda as: v" + v = 0 Una solucin de esta ecuacin es:
5.27
v = sen ( x - a ) En efecto: v' = cos ( x - a ) y v" = - sen ( x - a )
con a = cte.
Consideremos ahora dos races consecutivas de v ( x ): a y a + sen ( x - a )
x a a+ a+2
Las races consecutivas de las funciones de Bessel, a partir de la segunda, estn separadas entre s por un valor muy prximo a . Efectivamente, en el caso de J0, las primeras races caen en:
B B
x = 2,40,
B B
5,52,
8,6557,
11,7915, 14,9298, etc.
Por su parte, las de J1 estn ubicadas en x = 0, 3,83, 7,015, 10,17, 13,32, etc. Estos valores confirman, por una parte, la alternancia entre las races de una y otra funcin, y por la otra, que la separacin entre races consecutivas es muy prxima a . En efecto, en el primer caso, la separacin es respectivamente: 3,120, 3,1357, 3,1358, 3,1383, . . . O, tratndose de J1,
B B
3,185, 3,155, 3,150, . . . De todo lo visto, podemos extraer las siguientes conclusiones: A partir de la ecuacin: v u" - u v" = ( h - g ) . u v Si es una raz de u ( x ), es decir, u ( ) = 0, entonces v() = 0 v() = 0 u" ( ) = 0 u" ( ) = 0 Por idntico razonamiento, si es una raz de u ( x ), (5.31)
Por otra parte, a partir de la igualdad:

5.28
v ( ) u' ( ) - v ( ) u' ( ) = podemos dedudcir que si h g, entonces
( h - g ) . u ( x ) . v ( x ) dx
(5.35)
v ( ) u' ( ) - v ( ) u' ( ) 0 O, lo que es igual: v ( ) u' ( ) v ( ) u' ( ) De aqu podemos deducir tambin que, si v ( ) = 0, entonces v() 0
Supongamos ahora que elegimos el punto a coincidente con . Como dos races cualesquiera de la funcin Jo ( x ), que son a su vez races de u ( x ), difieren entre s un valor muy aproximado, pero no exactamente igual, a , nos encontraremos con una situacin aproximada a la que representa la figura siguiente:
B B
u (x) u' (x) a a +2 x
v (x)
La figura muestra que entre cada dos races de v ( x ) cae una raz de Jo ( x ). Ntese que la funcin v ( x ) conserva su signo entre y , lo que satisface el postulado que acompaa a la ecuacin (5.35):
B B
Si v ( x ) se mantiene positiva (o negativa) en todo el intervalo , el primer trmino de la ecuacin no puede ser igual a cero en nngn punto del mismo, por lo que la integral no puede ser nula en el intervalo, lo que implica que, necesariamente, debe ser h(x) g(x)
5.29
5.8 - Ortogonalidad de las Funciones de Bessel. En esta seccin probaremos que las funciones de Bessel forman un conjunto ortogonal en el intervalo entre x = 0 y x = 1. Para ello en la ecuacin diferencial de Bessel, (5.1), reemplazaremos la variable x por = ax, donde a es una constante arbitraria.
La ecuacin de Bessel queda entonces as: 2 . Jn " ( ) + . Jn ' ( ) + ( 2 - n 2 ) . J n ( ) = 0

(5.41)
Derivando Jn ( ) respecto de x obtenemos:

B B
d Jn ( )
B B
dx
B B
= d Jn ( ) d d dx
B B
= a
d Jn ( )
B B
d J 'n ( )
B B
d Jn ( ) = d
d Jn ( ) = a dx
B B
Derivando nuevamente, la derivada segunda es: d2 J n ( ) = d 2

P P B B P P B B B B
J"n ( )
B B
a2
P P
Jn' ( ) y Jn" ( ) son la derivadas respecto de x. Reemplazando en (5.41), la misma queda as: x2 . J n " ( ) + x . J n ' ( ) + ( a 2 x2 - n 2 ) J n ( ) = 0
P P B B B B P P P P P P B B
o tambin, dividiendo por x2:

P P
Jn " ( ) +
B B
1 J '() + n x
B B
a2 P P
n2
P P P
x2
P
Jn ( ) = 0
B B
(5.42)
Tambin aqu ensayaremos la solucin: u = Jn x

B B
Jn = u x-1/2
B B P P
Derivando respecto de x: Jn' = u' x-1/2 B B P P
1 2
u x-3/2
P P
Y derivando una vez ms: Jn" = u" x -1/2

B B P P
1 2
u' x -3/2
P P
1 2
u' x -3/2 +
P P
3 4
u x -5/2 =
P P
5.30
= u" x -1/2 u' x -3/2 +

P P P P
3 4
u x -5/2
P P
Por fin, al remplazar en la (5.42), la ecuacin diferencial queda ahora formulada as: u" x -1/2 u' x -3/2 + 3 u x -5/2 + u' x-3/2 1 u x -5/2 + a2 u x -1/2 n2 u x -5/2 = 0 4 2
P P P P P P P P P P P P P P P P P P
A continuacin, simplificaremos los dos trminos en u', multiplicaremos el resto por x agruparemos los coeficientes de u; con lo cual la ecuacin se reduce a la forma siguiente:
P
1/2
P
, y
u" +
1 4 x2
P
+ a2
P P P
n2
P P
u = 0
x2
P P
Reagrupando de nuevo, la ecuacin puede ahora sintetizarse como vemos aqu abajo: u" + a2 P P
4 n2 - 1
P P
u = 0
(5.43)
4x
P
Si en lugar del reemplazo = ax hubiramos hecho = bx, siguiendo idnticos pasos habramos llegado a una ecuacin totalmente similar, excepto por el nombre de la variable: v" + b2 P P
4 n2 - 1
P P
v= 0
(5.44)
4x
P
Las dos ltimas ecuaciones se pueden simbolizar as: u" + G ( x ) . u = 0 y donde v" + H ( x ) . v = 0 G ( x ) = a2 P P
4 n2 - 1
P P
4 x2
P
2 H ( x ) = b2 - 4 n - 1 4 x2
P P P P P
Si nos detenemos un instante para analizar lo hecho hasta aqu, podemos ver que estas ecuaciones son idnticas a las (5.29) y (5.30), a partir de las cuales obtuvimos en su momento la (5.34). Siguiendo ahora los mismos pasos que en aquella ocasin, pero integrando esta vez entre los lmites 0 y x, llegamos a la igualdad siguiente:
x x
( v u' - u v' )
0
( H - G ) . u v . dx
5.31
Si restamos G ( x ) de H ( x ), es decir: H - G = b2 - a 2
P P P P
y reemplazamos en la integral anterior, nos queda:

x x
y como
( v u' - u v' )
0
= (b - a )
P P P P
u v . dx
(5.45)
u = x y v = x
. Jn ( ax )
B B
. Jn ( bx )
B B
obtenemos finalmente:
x x
( v u' - u v' )
0
= ( b 2 - a2 )
P P P P
x . Jn ( ax ) . Jn ( bx ) . dx
B B B B
Las derivadas de u y v son, respectivamente: u' = 1 2 y Por tanto: v u' - u v' = 1 2 1 2 v u' - u v' = x a J'n ( ax ) Jn ( bx ) - b J'n ( bx ) Jn ( ax )
P P B B B B P P B B B B
x-1/2 Jn ( ax ) + a x1/2 J'n ( ax )

P P B B P P B B
v' =
1 2
x-1/2 Jn ( bx ) + b x1/2 J'n ( bx )

P P B B P P B B
Jn ( ax ) Jn ( bx ) + a x J'n ( ax ) Jn ( bx ) B B B B P P B B B B
Jn ( bx ) Jn ( ax ) - b x J'n ( bx ) Jn ( ax )
B B B B P P B B B B
Reemplazamos ahora en (5.45), y obtenemos

x
( b2 - a 2 )
P P P P
x Jn ( ax ) Jn ( bx ) x = x
B B B B
a J'n ( ax ) Jn ( bx ) - b J'n ( bx ) Jn ( ax )
B B B B B B B B
Al establecer el intervalo de integracin entre 0 y 1, la ecuacin anterior se transforma en la siguiente:

1
( b2 - a2 ) x Jn ( ax ) Jn ( bx ) dx = a J'n ( a ) Jn ( b ) - b J'n ( b ) Jn ( a )
P P P P B B B B B B B B B B B B
(5.46)
5.32
A continuacin analizaremos una por una las diferentes alternativas que pueden darse respecto de esta ecuacin. Primer caso: a y b son ceros de Jn ( x ) Entonces, el segundo miembro de la (5.46) es igual a cero:
B B
( b2 - a 2 )
P P P P
x Jn ( ax ) . Jn ( bx ) dx = 0
B B B B
Estamos en presencia de dos factores cuyo producto es cero. Existen dos posibilidades: Si ambas races son distintas, a b, o la integral es nula, o bien a = - b. Obviamente entonces, para que la integral
1
x Jn2 ( ax ) dx
B PB P
(5.47)
sea distinta de cero, la nica posibilidad es que: a2 = b2

P P P P
Esta es por tanto la condicin de ortogonalidad. Segundo caso: Una segunda alternativa consiste en que a y b sean ambas races de la ecuacin J'n ( x ) = 0,
B B
observando la (5.46) vemos que se reproduce una situacin idntica a la del primer caso. Es decir, que la condicin necesaria para que la integral (5.47) sea no nula es la misma que antes: a2 = b2
P P P P
Tercer caso: a y b son ambas races de la ecuacin Jn+1 ( x ) = 0

B B
(5.48)
De acuerdo con la (5.25), sabemos que d dx Al efectuar esta derivada llegamos a la igualdad: - n x -n-1 Jn ( x ) + x -n J'n ( x ) = - x -n Jn+1 ( x )
x -n Jn ( x )
P P B B
= - x -n Jn+1 ( x )
P P B B
Simplificando x -n en los dos miembros de la misma, obtenemos

P P
- n x -1 Jn ( x ) + J'n ( x ) = - Jn+1 ( x )
P P B B B B B B
5.33
Jn+1 ( x ) + J'n ( x ) = n x -1 Jn ( x )
B B B B P P B B
Pero a y b son races de la (5.48), es decir: Jn+1 ( a ) = 0

B B
Jn+1 ( b ) = 0,
B B
entonces: J'n ( a ) = n a-1 Jn ( a )

B B P P B B
o bien: a J'n ( a ) = n Jn ( a )
B B B B
De modo similar obtendramos: b J'n ( b ) = n Jn ( b )

B B B B
Multiplicando estas dos ltimas igualdades entre s, en forma cruzada, obtenemos: b J'n ( b ) n Jn ( a ) = a J'n ( a ) n Jn ( b )
B B B B B B B B
(5.49)
De aqu, simplificando n en ambos miembros y pasando todo al primero, resulta: b J'n ( b ) Jn ( a ) - a J'n ( a ) Jn ( b ) = 0
B B B B B B B B
El primer miembro de esta ltima igualdad coincide exactamente con el ltimo miembro de la (5.46). Por tanto, nuevamente deber ser nula la integral que aparece en ella, salvo que fueran a2 = b2
P P P
La conclusin obvia es que esta tercera alternativa es equivalente a los dos anteriores. Cuarta posibilidad: a y b son dos races, diferentes, de Jn-1 ( x ) = 0
B B
(5.50)
Segn la (5.28), sabemos que d dx Efectuando esta derivada llegamos a la expresin: n xn-1 Jn ( x ) + xn J'n ( x ) = xn Jn-1 ( x )
xn J n ( x )
P P B B
xn Jn-1 ( x )
P P B B
Simplificamos aqu xn antes de pasar todo al primer miembro:

P P
n x -1 Jn ( x ) + J'n ( x ) =
P P B B B B
Jn-1 ( x )
B B B B
n Jn ( x ) + x J'n ( x ) - x Jn-1 ( x ) = 0
B B B B
(5.51)
Como a y b son races de la (5.50), es decir:
5.34
Jn-1 ( a ) = 0
B B
Jn-1 ( b ) = 0,
B B
Al reemplazar cada una de estas igualdades en la ecuacin (5.51), sta queda, respectivamente igual a: n Jn ( a ) = - a J'n ( a )
B B B B
b J'n ( b ) = - n Jn ( b )
B B B B
Multiplicando miembro a miembro verificamos que: b J'n ( b ) n Jn ( a ) = a J'n ( a ) n Jn ( b )

B B B B B B B B
Pero ya habamos obtenido esta igualdad, en la ecuacin (5.41), por lo que tambin la cuarta posibilidad conduce a un resultado idntico a los tres primeros supuestos.
Quinta posibilidad: Consideremos ahora la ecuacin: x J'n ( x ) + C Jn ( x ) = 0

B B B B
en la que C es una constante arbitraria, en tanto que a y b son dos races distintas de la ecuacin, es decir: a J'n ( a ) + C Jn ( a ) = 0 y b J'n ( b ) + C Jn ( b ) = 0 O bien: a J'n ( a ) = - C Jn ( a ) y b J'n ( b ) = - C Jn ( b )
B B B B B B B B B B B B B B B B
Si como en los casos anteriores, multiplicamos miembro a miembro en forma cruzada, tendremos: a J'n ( a ) . C Jn ( b ) = b J'n ( b ) . C Jn ( a )
B B B B B B B B
Simplificando C, y pasando todo al primer miembro, obtenemos nuevamente una ecuacin que ya nos est resultando familiar: b J'n ( b ) . Jn ( a ) - a J'n ( a ) . Jn ( b ) = 0
B B B B B B B B
La conclusin es, una vez ms, que este caso es tambin similar a los anteriores. En resumen, si se cumple que a y b son races de cualesquiera de las cinco ecuaciones supuestas, la condicin necesaria para que la integral
1
x Jn ( ax ) Jn ( bx ) dx
B B B B
sea distinta de cero es en todos los casos la misma:
a 2 = b2
P P P P
5.35
Conclusin: Condiciones de ortogonalidad de las funciones de Bessel, en el intervalo {0;1} Recordemos que dos funciones, fj ( x ) y fk ( x ), de un conjunto son ortogonales entre s cuando
B B B B
= 0, fj ( x ) fk ( x ) dx =
B B B B
si j k si j = k
0,
En ocasiones, la ortogonalidad requiere la existencia de una funcin "complemento" o "peso", g ( x ), de la siguiente forma:
b
= 0, g ( x ) fj ( x ) fk ( x ) dx =
B B B B
si j k si j = k
0,
Tal es el caso de las funciones de Bessel, donde: El complemento es g ( x ) = x, como puede verse en los cinco casos expuestos. Las funciones fj ( x ) y
B B B B B B
fk ( x ) son, respectivamente:
B B B B B B
fj ( x ) = Jn ( ax ) = Jn ( j x ) y fk ( x ) = Jn ( bx ) = Jn ( k x )
B B B B B B B B
Partiendo de esto, si llamamos C al conjunto de las funciones Bessel de primera clase: C = [ J n ( i x ) ]

B B B B
constituye un Conjunto Ortogonal de funciones en el intervalo {0, 1}, cuyo peso es g ( x ) = x. Es decir que si: J n ( j x ) C
B B B B
Jn ( k x ) C
B B B B
entonces:
1
= 0,
B B B B B B B B
si j 2 k 2
B B P P B B P
x Jn ( j x ) Jn ( k x ) dx = 0, si j 2 = k 2
B B P P B B P P
5.9 - Desarrollo de funciones en serie de Funciones de Bessel: En el apartado anterior hemos demostrado que las funciones de Bessel Jn ( x ),
B B
con n = 0, 1, 2, 3, 4 . . .
constituyen un conjunto ortogonal de funciones en el intervalo 0 x 1, y que lo mismo ocurre en el caso de las familias de funciones siguientes
5.36
Jn ( i x ), Jn+1 ( i x ), Jn-1 ( i x ), J'n ( i x ), x J'n ( i x ) + C Jn ( i x ),

B B B B B B B B B B B B B B B B B B B B B B B B
Por consiguiente es posible en principio desarrollar una funcin f ( x ) como una suma constituida por los infinitos trminos de una cualquiera de las familias mencionadas, donde el parmetro i vara de cero a infinito. Por ejemplo: f(x) =
i=0
c i J n ( i x ) = c o J n ( o x ) + c 1 J n ( 1 x ) + c 2 J n ( 2 x ) + . . .
B B B B B B B B B B B B B B B B B B B B B B B B
En los textos especializados se encuentran tablas de los valores de Jn ( i x ).

B B B B
Esta sumatoria constituye de hecho una serie generalizada de Fourier, f(x) = en la cual: j ( x ) = J n ( j x )
B B B B B B
n =1
c n . n ( x ) = c 1 1 ( x ) + c 2 2 ( x ) + . . . + c n . n ( x ) + . . .
B B B B B B B B B B B B B B B B
(5.52)
k ( x ) = J n ( k x )
B B B B B B B B
Si multiplicamos ambos miembros de la serie (5.52) por k ( x ) y a continuacin integramos en un perodo que en forma genrica vamos a designar como a, b, resultar:
b
B B B B
f ( x ) . k ( x ) . d x =
b
cn
B
.
B
a b
k ( x ) . n ( x ) . d x =
B B B B
n=1
= c1
B B
a b
1 ( x ) . k ( x ) d x + c 2
B B B B B B
2 ( x ) . k ( x ) d x + . . . +
B B B B
+ ck
B B
k ( x ) . k ( x ) d x + . . .
B B B B
Pero
b
= Mk
B B
si k = n (5.53)
si k n
B B
k ( x ) . n ( x ) . d x
B B B B
= 0
De aqu llegamos a la siguiente igualdad, que nos permite definir los coeficientes ck.
b
f ( x ) . k ( x ) . d x =
B B B B
ck . M k
B B B B
(5.54)
B
5.37
Al introducir este resultado en la sumatoria podemos escribir, cambiando el nombre de la variable de integracin para evitar confusiones:
b
B B
f ( x ) = k ( x ) . 1 k =1 Mk
B B
f ( ) . k ( ) . d
B B B B
El segundo miembro en esta ecuacin constituye lo que se conoce como una Serie generalizada de Fourier. Volviendo a la ecuacin (5.52), de modo similar tenemos:
b
= Mk
B B
si k = j si k j
j ( x ) . k ( x ) . d x
B B B B
= 0
haciendo los reemplazos adecuados al caso presente y teniendo en cuenta el peso, x, hallamos la igualdad:
1
0
B B
x . Jn2 ( j x ) . d x = Mk
B PB P B B B B
Obsrvese que Mk es la Norma del conjunto ortogonal de funciones. Por su parte, efectuando los reemplazos correspondientes en la ecuacin (5.54)
b
f ( x ) . k ( x ) . d x =
B B B B
ck . M k
B B B B
(5.54)
para la serie de funciones de Bessel tenemos:

1
f ( x ) . x . J n ( j x ) . d x =
B B B B B B
ck . M k
B B B B
De aqu podemos deducir el valor de los coeficientes:

1
ck =
B B
1 Mk
B B
x . f ( x ) . J n ( j x ) . d x
B B B B B B
5.10 - Funciones de Bessel de segunda clase. En el apartado 5.2 se vio que cuando es un nmero entero n, las funciones Jn ( x ) y J-n ( x ) mantienen una relacin de linealidad entre s, por lo que para resolver la ecuacin diferencial respectiva se haca necesario recurrir a las llamadas Funciones de Bessel de Segunda Clase.
B B B B
5.38
Es decir que cuando es un nmero entero, la solucin general de la ecuacin de Bessel no puede obtenerse como una suma ponderada de Jn ( x ) y J-n ( x ), como hacamos en el caso contrario(1).
B B B B TP PT
En este apartado trataremos la situacin particular mencionada. Veamos primero el caso en que n = 0. En tal supuesto, la ecuacin diferencial de Bessel, x2 y" + x y' + ( x2 - n2 ) y = 0
P P P P P P
se reduce a x y + y + x y = 0 La ecuacin indicial ( - 1 ) c o + co - n2 co = 0

B B B B P P B B
se reduce a su vez a: 2 - + = 0
P P
1 = 2 = 0
B B B B
Es decir que la ecuacin indicial tiene una raz doble igual a cero. Esto nos obliga a buscar un camino diferente para encontrar la solucin. En el caso general, para que la ecuacin (3.6)
2 = - ( o - 1 ) ( o - 1 ) - 4 o 2
B B B B P P B
tenga una raz real doble, es necesario y suficiente que se cumpla la igualdad: ( o - 1 ) 2 = 4 o
B B P P B
De aqu, se deduce que = - ( o - 1 ) 2

B B
= 1 - o 2
B B
Una solucin particular de la ecuacin diferencial de segundo orden con coeficientes variables puede ser, como hemos visto, de la forma (3.7): y1 ( x ) = x ( co + c1 x + c2 x2 + . . . ) (2)
B B P P B B B B B B P P TP PT
La siguiente es asimismo una segunda solucin particular de la ecuacin y2 ( x ) = u ( x ) . y 1 ( x )

B B B B
(5.55)
(1)
TP
Por supuesto, esta conclusin no invalida todo lo visto hasta aqu sobre las funciones de Bessel de subndice n, entero.
PT
(2)
TP
Como regla general, descomponer la solucin en el producto de una cierta potencia de x, convenientemente elegida, por una serie geomtrica en x, facilita grandemente el trabajo posterior
PT
5.39
en donde u (x ) es una funcin a determinar(3).

TP PT
Para obtener la funcin u ( x ) debemos calcular las derivadas primera y segunda de y2 ( x ):

B B
y2 ( x ) = u( x ) . y1 ( x ) + u ( x ) . y1 ( x )
B B B B B B
(5.56)
B B
y2 ( x ) = u( x ) . y1 ( x ) + 2 u ( x ) . y1 ( x ) + u ( x ) . y1 ( x )
B B B B B B
Reemplazando estos valores en la ecuacin general (3.2), con a = 0, x2 y" + ( x ) . y' + ( x ) . y = 0

P P
(5.57)
tenemos: x2 [u( x ) . y1 ( x ) + 2 u ( x ) . y1 ( x ) + u ( x ) . y1 ( x )] +
P P B B B B B B
+ ( x ) [u( x ) . y1 ( x ) + u ( x ) . y1 ( x )] + ( x ) . u ( x ) . y1 ( x ) = 0
B B B B B B
Reagrupando adecuadamente los trminos de esta ecuacin, x2 [u( x ) . y1 ( x ) + 2 u ( x ) . y1 ( x )] + ( x ) . u( x ) . y1 ( x ) +

P P B B B B B B
+ u ( x ) [ x2 . y1 ( x ) + ( x ) . y1 (x ) + ( x ) . y1 ( x )] = 0
P P B B B B B B
vemos que la expresin que aparece multiplicando a u ( x ) es justamente el primer miembro de la (5.57) , por lo que la misma es igual a cero. Al simplificarla, queda solamente: x2 [u( x ) . y1 ( x ) + 2 u ( x ) . y1 ( x )] + ( x ) . u( x ) . y1 ( x ) = 0
P P B B B B B B
Desarrollemos ahora ( x ) en serie de potencias ( x ) = ( o x + 1 x2 + 2 x3 + 3 x4 + . . . )

B B B PB P P P B B P P B B P P
y dividamos todo por x2 y1, para obtener en forma explcita las derivadas de u ( x ):
P P B B
u ( x ) + ( 2 y1 + o + 1 + 2 x + 3 x2 + . . . ) u ( x ) = 0 y1 x
B B B B B PB P B B B B P P B B
(5.58)
Por otro lado, la derivada primera de y1 ( x ) = x ( c o + c 1 x + c 2 x2 + . . . )

B B P P B B B B B B P P
es:
y1 ( x ) = x1 ( co + c1 x + c2 x2 + . . . ) + x ( c1 + 2 c2 x + . . . )
B B P P B B B B B B P P P P B B B B
Esta ecuacin se puede simplificar unificando las potencias de x dentro y fuera de los parntesis, asi: y1 ( x ) = x1 ( co + c1 x + c2 x2 + . . . ) + x1 ( c1 x + 2 c2 x2 + . . . )
B B P P B B B B B B P P P P B B B B P P
Esto nos permite sacar factor comn x1, con lo cual tenemos:
P P
y1 ( x ) = x1 [ ( co + c1 x + c2 x2 + . . . ) + ( c1 x + 2 c2 x2 + . . . )]
B B P P B B B B B B P P B B B B P P
(3)
TP
La determinacin de la funcin u (x) requiere un proceso sencillo pero largo. Si se desea obviarlo, sugerimos avanzar a la pgina 5.41, donde damos una apretada sntesis del mismo.
PT
5.40
A continuacin vamos a dividir ambos miembros por y1 ( x ):

B B
y1 ( x ) = y1 ( x )
B B B B
x1 [ co + c1 x + c2 x2 + . . . ) + ( c1 x + 2 c2 x2 + . . . )]
P P B B B B B B P P B B B B P P
x ( c o + c 1 x + c 2 x2 + . . . )
P P B B B B B B P P
Simplificando: y1 ( x ) = y1 ( x )
B B B B
1 x
( c o + c 1 x + c 2 x2 + . . . ) + ( c 1 x + 2 c 2 x2 + . . . )
B B B B B B P P B B B B P P
( c o + c1 x + c2 x + . . . )
B B B B B B P P
Pero la serie que multiplica a en el numerador es igual a la que aparece en el denominador; por lo que, simplificando nuevamente, obtenemos: y1 ( x ) = y1 ( x )
B B B B
( c 1 x + 2 c 2 x2 + . . . )
B B B B P P
+ (x)
x ( c o + c1 x + c2 x + . . . )
B B B B B B P P
En esta ecuacin, el significado de ( x ) es obvio. Reemplazando en la (5.58) el resultado anterior, se llega a la igualdad: u + [ 2 Llamando: ( x ) = ( x ) + 1 + 2 x + 3 x2 + . . . +
B PB P B B B B P P
o
B B
+ 1 + 2 x + 3 x2 + . . . + ( x ) ] . u = 0
B PB P B B B B P P
alcanzamos una simplificacin an mayor, que nos permite escribir: u + [ 2 x + o

B B
+ ( x ) ] . u = 0
(5.59)
Hemos hallado oportunamente que cuando la ecuacin indicial deriva en una raz doble para la constante , el valor de sta ltima es: = 1 - 0
B B
2 + 0 = 1
B B
2 Reemplazando en la ecuacin (5.59): u + [ 1 x Dividamos ambos miembros por u'. Si pasamos los trminos entre parntesis al segundo miembro, resulta: u u = 1 x + (x) + ( x ) ] . u = 0
Al llegar aqu, podemos hacer:
5.41
u u
d u' u'
Entonces, integrando, resulta: ln u = - ln x + ( x ) dx Y por tanto: u = 1 x Como ( x ) es un polinomio en x, la integral que aparece en el exponente es asimismo un polinomio en x; ello permite desarrollar la exponencial en esta ecuacin como una serie de potencias de x. Si llamamos hi a los coeficientes de la misma, podemos escribir simblicamente:
B B
e
P
(x) dx
u =
P
1 x
( 1 + ho x + h1 x2 + h2 x3 + . . . ) =
B B B B P P B B P P
1 x
+ ho + h1 x + h2 x2 + . . .
B B B B B B P P
Integrando nuevamente, podemos tambin poner: u = ln x + k1 x + k2 x2 + . . .

B B B B P P
(5.60)
El significado de los coeficientes kj es tambin suficientemente explcito:

B B
k1 = h o ,
B B B B
k2 =
B B
h1
B B
, etc.
2 ______________________________________________________________________________ Recapitulacin: La ecuacin de Bessel x2 y" + x y' + ( x2 - n2 ) y = 0

P P P P P P
cuando n es igual a cero se reduce a x y + y + x y = 0 La ecuacin indicial no nos sirve en este caso para avanzar hacia la solucin general. Por tanto, recurriremos a la ecuacin diferencial de segundo orden con coeficientes variables. Una solucin particular de la misma es: y1 ( x ) = x ( c o + c 1 x + c 2 x2 + . . . )
B B P P B B B B B B P P
O bien
y2 ( x ) = u ( x ) . y1 ( x )
B B B B
Donde la funcin u ( x ), como se ha visto, es de la forma u = ln x + k1 x + k2 x2 + . . .

B B B B P P
(5.60)
5.42
Al reemplazar este valor en la ecuacin anterior, obtenemos y2 ( x ) = y1 ( x ) . ln x + x ( k1 x + k2 x2 + . . . ) . ( co + c1 x + c2 x2 + . . . )

B B B B P P B B B B P P B B B B B B P P
Ahora estamos en condiciones de seguir en busca de la solucin y2 ( x ).

B B
Al efectuar el producto de las dos series entre parntesis, aparece una nueva serie de potencias en x. Llamando Am al coeficiente de la potencia emsima, arribaremos a la expresin simblica siguiente:
B B
y2 ( x ) = y1 ( x ) . ln x + x
B B B B P
m=1
m xm
B B P P
Cuando = 0, la ecuacin anterior se simplifica as:
y2 ( x ) = y1 ( x ) . ln x +
B B B B
m xm
B B P P
m=1
Aplicaremos ahora este resultado a la solucin de la ecuacin de Bessel: y2 ( x ) = Jo ( x ) . ln x +

B B B B
m xm
B B P P
(5.61)
m=1
Derivando sucesivamente, obtenemos:
y'2 ( x ) = J'o ( x ) . ln x + Jo ( x ) + m m xm-1 x m=1

B B B B B B B B P P B B B B B B B B B B
y"2 ( x ) = J"o ( x ) . ln x + J'o ( x ) + x J'o ( x ) - Jo ( x ) + m (m - 1) m xm-2 x x2 m=1

B B P P P P
y"2 ( x ) = J"o ( x ) . ln x + 2 J'o ( x ) x

B B B B B B
Jo ( x ) + m (m - 1) xm-2 m x2 m=1
B B B B P P P P
Como para n = 0 la ecuacin de Bessel es: x y" + y' + x y = 0 y por su parte, Jo ( x ) es solucin de la misma, reemplazando en la ltima igualdad las sucesivas derivadas de y ( x ) hallamos:
B B
( x Jo + J'o + x J"o ) ln x + Am xm+1 +

B B B B B B B B P P
Jo
B B
m Am xm-1 + 2 J'o B B P P B B
Jo
B B
m=1
m=1
m ( m - 1 ) Am xm-1 = 0
B B P P
m=1
5.43
Al ser Jo una solucin particular, la suma

B B
x Jo + J'o + x J"o
B B B B B B
es necesariamente nula. Por lo tanto, luego de simplificar y ordenar los trminos, resulta:
2 J'o ( x ) +
B B
Am x
B B P
m+1
P
m Am x
B B P
m-1
P
m ( m - 1 ) Am xm-1 = 0
B B P P
m=1
m=1
m=1
Como veremos enseguida, nos conviene agrupar las dos ltimas sumatorias, que slo difieren entre s en el factor multiplicador de cada trmino; obtenemos: 2 J'o ( x ) +
B B
Am xm+1
B B P P
( m + m2 - m ) Am xm-1 = 0
P P B B P P
m=1
m=1
Efectivamente, ahora podemos simplificar ms y menos m:
2 J'o ( x ) +
B B
Am x
B B P
m+1
P
m2 Am xm-1 = 0
P P B B P P
(5.62)
m=1
m=1
Hemos visto (5.7) que J ( x ) = x

B B P P
m=0
( - 1 )m x2m
P P P P
22m+ m! ( + m + 1 )
P P
Si = 0, ser:
J0 ( x ) =
B B
m=0
( - 1 )m x2m
P P P P
22m m! ( m + 1 )
P P
Cuya derivada respecto de x es:
J0 ( x ) =
B B
m=0
( - 1 )m 2m x2m-1
P P P
22m m! ( m + 1 )
P P
Con el objeto de simplificar esta ecuacin, le introduciremos previamente ciertas modificaciones. Comenzaremos por recordar que: ( m + 1 ) = m! = m ( m-1 )! Entonces
J0 ( x ) =
B B
m=0
P
( - 1 )m m x2m-1
P P P
22m-1 m ( m 1 )! m!
P
Ahora s, simplificaremos "m" en numerador y denominador:

5.44
J0 ( x ) =
B B
m=0
( - 1 )m x2m-1
P P P P
22m-1 ( m 1 )! m!
P P
Adems, como el primer trmino de la sumatoria, es decir, el que corresponde al subndice m = 0, es nulo porque en su denominador aparece el factorial de - 1, cuyo valor es infinito, podemos hacer tambin:
J0 ( x ) =
B B
m=1
( - 1 )m x2m-1
P P P P
(5.63)
2
P
2m-1
P
( m 1 )! m!

Al introducir este valor de la derivada en la (5.62), obtenemos:
m=1
( - 1 )m x2m-1
P P P P
Am x
B B P
m+1
P
m2 Am xm-1 = 0
P P B B P P
(5.64)
P
2
P
2m-2
P
( m 1 )! m!
m=1
m=1
Como esta expresin debe tener validez general, cualquiera sea el valor de x, la suma de los trminos correspondientes a cada potencia de dicha variable x deber ser siempre igual a cero. Ntese que como m, sbndice, es un entero, las potencias de x son todas impares. Clculo de los coeficientes Am:
B B
A partir de la condicin mencionada en el prrafo anterior podemos calcular los sucesivos coeficientes Am, como veremos ms abajo. Pero previamente, como para un dado valor de m, los exponentes de x son diferentes en cada sumatoria, para evitar confusiones cambiaremos el nombre de los subndices, as:
B B
m=1 2
P
( - 1 )m x2m-1
P P P P
h=1
Ah x
B B P
h+1
P
k2 Ak xk-1 = 0
P P B B P P
(5.65)
P
2m-2
P
( m 1 )! m!
k=1
Coeficientes de x1:
P P
Los coeficientes de la potencia primera de x en cada una de las sumatorias son, respectivamente: 2m - 1 = 1 h + 1 = 1 k - 1 = 1 m = 1
Luego h = 0 Luego k = 2
Reemplazando estos valores en la (5.62), obtenemos: - x

P P P P
+ A0 x + 4 A2 x = - x + A0 x + 4 A2 x = 0
B B P P B B P P B B P P B B P P
2 ( 0 )! ( 1 )!
P P
Pero, como las sumatorias en la (5.64) comienzan todas a partir del subndice m = 1, A0 es nulo. Entonces
B B
5.45
4 A2 = 1
B B P P
A2 = 1 / 4 = 0,25
B B
Coeficientes de x3:
P P
2m - 1 = 3 h + 1 = 3 k - 1 = 3
m = 1 h = 2 k = 4
Reemplazando en la (5.65), obtenemos en este caso:

P
x3
P P
+ A2 x3 + 16 A4 x3 = ( 1/8 + 1/4 + 16 A4 ) x3 = 0
B B P P B B P P B B P P
2 ( 1 )! ( 2 )!
P P
O sea, 16 A4 = - 1/8 - 1/4 = - 3/8

B B P P
A4 = - 3 / 128
B B
Coeficientes de x5:
P P
2m - 1 = 5 h + 1 = 5 k - 1 = 5 Reemplazando en la (5.65):
P
m = 3 h = 4 k = 6
- x5
P P
A4 x5 + 36 A6 x5 = 0
B B P P B B P P
2 ( 2 )! ( 3 )!
P P
De aqu, 36 A6 = 1/192 + 3/128 = 4/192 = 1/48

B B P P
A6 = 11 / 13824,
B B
etc.
Observemos que los nicos coeficientes que aparecen son aquellos cuyo ndice es par, es decir: A2, A4, A6, A8, etc.
B B B B B B B B
Reemplazando sus valores en la solucin propuesta, y2 ( x ), sta ltima queda as:

B B
y2 ( x ) = Jo ( x ) . ln x +
B B B B
1 4
x2
P P
3 128
x4 +
P P
11 13824
x6 . . .
P P
(5.66)
que trataremos de expresar como una sumatoria. Para ello, teniendo en cuenta que, como acabamos de ver, m es siempre par, pondremos m = 2r Con lo cual se obtiene para los coeficientes la siguiente frmula general: Am =
B B
1 +
1 2
1 3
+ ... +
1 r
( 1 ) r -1
P P P P P
22r ( r! ) 2
P
5.46
Verificaremos a continuacin su validez. En efecto, si m = 2 (Y por tanto, r = 1): A2 =

B B P
( 1 ) 0
P P
=
2
P P
1 4
2 ( 1! )
P
Si m = 4, r = 2: A4 =
B B
1 +
1 2
P
( 1 ) 2 ( 2! )
P P
=
2
P
3 2 . 16 . 4
3 128
Si m = 6, r = 3: A6 =
B B
1 +
1 2
1 3
( 1 ) 2
P P
=
2
P P
11 6 . 64 . 36
11 13824
2 ( 3! )
P P
etc. Ahora podemos reemplazar el valor de los coeficientes Am en la (5.66):

B B
y2 ( x ) = Jo ( x ) . ln x +
B B B B
r=1
1 +
1 + ... + 1 2 r
( 1 ) r -1
P P
x 2r
P P
2
P
2r
P
( r! )
P
2
P
Expresando en forma simblica la suma que aparece entre parntesis, esta ecuacin se puede expresar tambin como
r
y2 ( x ) = Jo ( x ) . ln x +
B B B B
r=1
p =1
1
p
( 1 ) r -1
P P
x 2r
P P P
2
P
2r
P
( r! )
P
Como Jo ( x ) e y2 ( x ) son linealmente dependientes entre s, es necesario recurrir a una solucin independiente, como por ejemplo:
B B B B
Yo ( x ) = [ y2 ( x ) + Jo ( x ) ]
B B B B B B
Se conviene en hacer: = 2/ y = - ln 2
siendo la llamada Constante de Euler, cuya definicin y valor son: = lim
k=1
1
k
- ln = 0,5772156 ...
r
Yo ( x ) =
B B
Jo ( x ) ln x + - ln 2 +
B B
r=1
p =1
1 p
( 1 ) r -1
P P
x 2r
P P
2
P
2r
P
( r! )
P
5.47
Yo ( x ) =
B B
Jo ( x )
B B
ln x + 2
r=1
( 1 ) r -1
P P
x 2r
P P P
p =1 p
2 ( r! )
P P P
2r
Esta ecuacin es la que se conoce como Funcin de Bessel de Segunda Clase, o Funcin de Neumann, de orden cero. Otras funciones usuales son: La Funcin de Bessel de Segunda Clase y orden , o Funcin de Neumann de Orden , que se define as: Y ( x ) =
B B
1 sen
J ( x ) cos - J ( x )
B B B B
Funcin de Bessel de 2 Clase (Funcin de Neumann)
Y la de orden n: Yn ( x ) = lim Y ( x )
B B B B
5.11 - Solucin general de la Ecuacin diferencial de Bessel: Una solucin general de la ecuacin de Bessel para cualquier valor de , es la siguiente: y ( x ) = C1 J ( x ) + C2 Y ( x )
B B B B B B B B
Donde C1 y C2 son constantes complejas arbitrarias.

B B B B
Otra solucin general est basada en las funciones de Bessel de tercera clase, o funciones de Hankel: H1, ( x ) = J ( x ) + i Y ( x )
B B B B B B
Y H2, ( x ) = J ( x ) - i Y ( x )
B B B B B B
Funciones de Hankel
En este caso, la solucin completa es de la forma: y ( x ) = A . H1, ( x ) + B . H2, ( x )

B B B B
en donde A y B son tambin constantes complejas arbitrarias. Las funciones de Hankel son de aplicacin para resolver problemas de la Fsica en los cuales aparecen magnitudes en cuadratura.
5.48
5.12 - Cuadro Sinptico. Ecuacin Diferencial de Bessel x2 y" + x y' + ( x2 - 2 ) y = 0

P P P P P P
Frmula general
, nmero real no negativo J ( x ) , funcin de Bessel de primera clase

B B
Solucin particular
J ( x ) = x
B B P
m=0
( - 1 ) m x2m
P P P P
22m+ m! ( + m + 1 )
P P B B B B
Solucin general Caso particular: = n, , nmero entero Solucin general
y ( x ) = A J ( x ) + B J ( x ) A y B, nmeros complejos arbitrarios x2 y" + x y' + ( x2 - n2 ) y = 0

P P P P P P
No aplicable cuando es un entero
y ( x ) = A J ( x ) + B Y ( x )
B B B B
Y ( x ) , Funcin de Bessel de segunda clase

B B
Funciones de Bessel
Funcin de Bessel de J ( x ) = x primera clase

B B P P
m=0
( - 1 ) m x2m
P P P P
22m+ m! ( + m + 1 )
P P
Idem, expresin integral n, entero no negativo Caso particular: = 0,5
Jn ( x ) =
B B
cos ( n - x sen ) d
0,5
P P
0 J0,5 ( x ) =
B B
sen x
x
B B
J- 0,5 ( x ) =
B B
2
P
0,5
cos x
P
Funcin de Bessel de 1 Y ( x ) = segunda clase o sen Funcin de Neumann

B B
[ J ( x ) cos - J ( x ) ]
B B
Funciones de Hankel
H1, ( x ) = J ( x ) + i Y ( x )
B B B B B B
H2, ( x ) = J ( x ) - i Y ( x )
B B B B B B
Utiles cuando en las ecuaciones aparecen trminos imaginarios.
5.49
El presente captulo 5, dedicado a la Ecuacin diferencial de Bessel y a su solucin, no pretende ser exhaustivo en el estudio de las Funciones de Bessel. Todo lo contrario, la multitud de aplicaciones de las mismas como solucin de una cantidad de ecuaciones que se presentan en los diferentes dominios de la Fsica, ha llevado al desarrollo de otras alternativas de una riqueza tal, en cantidad y en propiedades, que excede largamente los propsitos de este tratado, pero que pueden ser consultadas tanto en los textos especializados como tambin en muchos libros dedicados al Clculo Superior.
5.13 - Problemas. Las Ecuaciones de Bessel aceptan en general soluciones comunes, lo que las hace particularmente interesantes para resolver muchos problemas de la Fsica. Y aunque dichas soluciones pueden parecer difciles o largas en su aplicacin, existen en general tablas que ayudan grandemente a su evaluacin. En este apartado veremos primero cmo resolver algunas Ecuaciones de Bessel, y a continuacin, cmo transformar ciertas ecuaciones diferenciales en Ecuaciones de Bessel, para de esta forma simplificar el problema de hallar la solucin.
5.13.1 - Resolver las ecuaciones siguientes: a) x2 y" + x y' + ( x2 - 1 ) y = 0 4 Solucin:

P P P P
Aqu es: = 0,5 y por lo tanto, la solucin general es: y ( x ) = A J1/2 ( x ) + B J1/2 ( x ) = A
B B B B
2 x
sen x + B
2 x
cos x
Si hacemos: A = 2 i y B = 2 ,
la ecuacin anterior se transforma en la siguiente:

-
y(x) = x
P
1 2
( cos x + i sen x ) = x
P
1 2
. e ix
P P
b) x2 y" + x y' + ( x2 - 25 ) y = 0 4 R: y ( x ) = A J2,5 ( x ) + B J- 2,5 ( x )

P P P P B B B B
5.50
c) x2 y" + x y' + ( x2 - 9 ) y = 0
P P P P
R: y ( x ) = A J3 ( x ) + B Y3 ( x )
B B B B
d) x2 y" + x y' + x2 y = 0
P P P P
Solucin: En este caso es igual a cero:
= 0
2 4 6 R: y ( x ) = J0 ( x ) = 1 - x + x - x + . . . 4 64 2304
P P P P P P B B
5.13.2 - Modificar las ecuaciones siguientes para transformarlas en ecuaciones de Bessel, y hallar la solucin: a) x2 y" + 5 x y' + x2 y = 0
P P P P
Solucin: Hagamos
P P
y = z.x
P P P
y ' = x -1 z + z ' x
y" = ( - 1 ) x -2 z + 2 x -1 z' + x z"

P P P P P P
Reemplazando en la ecuacin obtenemos: ( - 1 ) x z + 2 x +1 z' + x + 2 z" + 5 x z + 5 z ' x + 1 + z . x + 2 = 0

P P P P P P P P P P P P
Podemos simlificar x , que aparece en todos los sumandos:

P P
( - 1 ) z + 2 x z' + x 2 z" + 5 z + 5 z ' x + z . x 2 = 0

P P P P
A continuacin ordenaremos esta ecuacin: x 2 z" + ( 2 + 5 ) x z ' + [ ( - 1 ) + 5 + x 2 ] z = 0

P P P P
De aqu: Y tambin:
2+5 = 1
= -2
(-1)+ 5 = -4
De donde, reemplazando, obtenemos: x 2 z" + x z ' + ( x 2 - 4 ) z = 0

P P P P
Que es una ecuacin de Bessel, en z, y en la cual, n es igual a 2. La solucin es: z ( x ) = A J 2 ( x ) + B Y2 ( x )

B B B B
Finalmente, como y = z . x , la solucin general, expresada como funcin de y, es:

P P P P
y ( x ) = A x J2 ( x ) + B x Y2 ( x )
P P B B P P B B
Donde A y B son dos constantes complejas arbitrarias. b) x2 y" + 6 x y' + x2 y = 0

P P P P
5.51
c) x2 y" + 3 x y' + x2 y = 0 2
P P P P
5.13.3 - Resolver las ecuaciones siguientes transformndolas previamente en ecuaciones de Bessel: a) x2 y" + x y' + m2 ( x2 - 2 ) y = 0
P P P P P P P P
Hagamos el reemplazo: t = mx Ahora, y ( x ) es una funcin de t: y ( x ) = y ( m -1 t ).

P P
Las derivadas de y respecto de t son: dy dt y

P P
= dy dx
. dx dt
= m -1 . y ' ( x )
P P
d2 y = d y ' . d x d t2 dx dt
P P
= m -2 . y " ( x )
P P
Reemplazamos ahora x, y '( x ) e y " ( x ) en la ecuacin original, con lo cual esta queda expresada como una funcin de t: m2 t 2 y" + m t y' + ( t 2 - 2 m 2 ) y = t 2 d 2 y + t d y + ( t 2 - 2 ) y = 0 m2 m dt2 dt
P P P P P P P P P P P P P P P P P P P P P P
Hemos obtenido as una ecuacin de Bessel, en la cual la variable es t, y el valor de : = m La solucin general es, por lo tanto: y ( t ) = A J ( t ) + B Y ( t )
B B B B
b) x2 y" + 5 x y' + 2 x2 y = 0
P P P P
El reemplazo adecuado en este caso es:
y = zx
P P
Con esta sustitucin, la ecuacin diferencial resulta (Ver problema 5.2 a): x 2 z" + ( 2 + 5 ) x z ' + [ ( - 1 ) + 5 + 2 x 2 ] z = 0
P P P P
De aqu:
2+5 = 1
= 2
Reemplazando en la ltima ecuacin, tenemos: x 2 z" + x z ' + 2 ( x 2 - 2 ) z = 0

P P P P
5.52
Esta ltima ecuacin es similar a la del problema 5.3 a, y por tanto se resuelve de manera idntica a aquella.
c) x y" + m y ' + n y = 0 Esta ecuacin aunque diferente de la ecuacin diferencial de Bessel, para ciertos valores de m y n acepta como solucin a la funcin de Bessel de orden cero: y ( t ) = J0 ( t )
B B
En efecto, hagamos el reemplazo: x = t 2. Como

P P
2 4 6 y(t) = 1 - t + t - t + ... 4 64 2304 entonces 2 x3 + ... y(x) = 1- x + x 4 64 2304 Derivando: 2 y ' ( x) = - 1 + x - x + ... 4 32 768 Tambin 2 x . y"( x ) = x - x + ...
P P P P P P P P P P P P P P P P
32
384
Reemplazando en la ecuacin original podemos arribar a las siguientes conclusiones: Para que el trmino independiente sea nulo, deberan ser m = 4 o bien, n = 1/4. Para que los trminos en x se anulen, deberan ser: m = 1 y n = 1/4. Estos ltimos valores anulan tambin el trmino en x2 y sucesivos. En definitiva, la ecuacin dada, en la que
P P
m =1
n = 1 4 donde t = x
tiene como solucin y ( t ) = J0 ( t )

B B
5.53
Problemas resueltos con Matlab: Calcular el valor de la funcin de Bessel de primera clase y orden 2, para x = 3. Comparar el resultado con la suma de los cuatro primeros trminos de la serie:
% Utilizaremos la funcin Matlab: BESSELJ (Nu, x): BESSELJ(2, 3) ans = 0.4861
A partir de los primeros cuatro trminos de la serie:
J2 ( 3 ) = 3
B B P
2
P
m=0
P P
( - 1 )m x2m
P P P P
22m+2 m! ( m + 2 ) !
P P
J2 ( 3 ) =
B B P
32
4
P
32 . 32
P P P P
+
P
32 . 34
P P P P
32 . 36
P P P P
= 0,4836
2 2!
B
2 3!
P P B
2 2! 4!
P
2 3! 5!
P P
Calcular el valor de la funcin J0,3 ( x ) para los primeros valores enteros de x. Hallar por aproximaciones sucesivas los primeros cuatro ceros, y determinar la separacin entre ellos.
BESSELJ (0.3, 1) ans = 0.7402 BESSELJ (0.3, 2) ans = 0.4257 BESSELJ (0.3, 3) ans = -0.0673 BESSELJ (0.3, 4) ans = - 0.3638 BESSELJ (0.3, 5) ans = -0.2968 BESSELJ (0.3, 6) ans = 0.0058 BESSELJ (0.3, 7) ans = 0.2567 BESSELJ (0.3, 8) ans = 0.2538 BESSELJ (0.3, 9) ans = 0.0317 BESSELJ (0.3, 10) ans = - 0.1946 BESSELJ (0.3, 11) ans = - 0.2289 BESSELJ (0.3, 12) ans = -0.0589 BESSELJ (0.3, 13) ans = 0.1495 Ubicacin ceros: Distancia entre ceros
BESSELJ (0.3, 2.855) ans = -4.2812e-004
3,127
BESSELJ (0.3, 5.982) ans = -7.2278e-005
3,137
BESSELJ (0.3, 9.119) ans = 8.9611e-005
3,140
BESSELJ (0.3, 12.259) ans = 6.4857e-005
5.54
Utilizando exclusivamente Matlab, verificar, para x = /2, el cumplimiento de la frmula: J0,5 ( x ) =

B B
2
P
0,5
P
sen x
0,5
P P
x 4
P
J0,5 ( ) = 2
B B
2
P
sen 2
BESSELJ (0.5, pi/2) ans = 0.6366
Verificacin en el segundo miembro de la igualdad:

a b J a = (pi)^2 = 9.8696 b= 4/a = 0,4053 J = b^0.5 = 0.6366
% Recordar que sen / 2 = 1
Repetir el problema anterior, siendo x = /3, para verificar la frmula: J- 0,5 ( x ) =

B B
2
P
0,5
P
cos x
BESSELJ (-0.5, pi/3) ans = 0.3898
Verificacin:
x = (pi)^3 x = 1.0472 y = cos (x) y = 0.5000 a = 2*3/(pi^2) a = 0.6079 b = a^0,5 b = 0.7097 J = b*y J = 0.3898
Calcular la funcin de Bessel de segunda clase y orden cero, para x = :

BESSELY (0, pi) ans = 0.3284
5.55
Construccin de tablas de funciones de Bessel. Ejemplos:

NU = [0] NU = 0 X = [0 0.2 0.4 0.6 0.8 1] X= 0 0.2000 0.4000 0.6000 0.8000 1.0000 BESSELJ(NU, X) ans = 1.0000 0.9900 0.9604 0.9120 0.8463 0.7652 De otra forma: BESSELJ (0, (0 : .2 : 1)) ans = 1.0000 0.9900 0.9604 0.9120 0.8463 0.7652
Construmos ahora una tabla idnticapero para otros valores de :

BESSELJ (1, (0 : .2 : 1)) ans = 0 0.0995 0.1960 0.2867 0.3688 0.4401
% Tabla de Funciones Bessel para = 2 BESSELJ (2, (0 : .1 : 1)) ans = Columns 1 through 7 0 0.0012 0.0050 0.0112 0.0197 Columns 8 through 11 0.0588 0.0758 0.0946 0.1149 Observemos tambin la forma siguiente:
BESSELJ (1, (0 : .2 : 1) ' ) ans = 0 0.0995 0.1960 0.2867 0.3688 0.4401
0.0306
0.0437
Ntese que en este caso, a la inversa de los ejemplos anteriores, los valores para constante aparecen en las columnas. La modificacin se debe al agregado del apstrofe entre los dos
parntesis de cierre. La utilidad del cambio la veremos en el ejemplo siguiente, que permite construir una tabla (O matriz) bidimensional:
5.56
Construccin de la tabla de valores para la funcin de Bessel de primera clase, para: Rango de la variable : De 0 a 0,5, en pasos de 0,1. Rango de la variable x: De 0 a 1, en pasos de 0,2. BESSELJ (0 : .1 : .5, (0 : .2 : 1)')
ans = 1.0000 0.9900 0.9604 0.9120 0.8463 0.7652 0 0.8274 0.8626 0.8573 0.8248 0.7708 0 0.6815 0.7633 0.7932 0.7902 0.7615 0 0.5542 0.6666 0.7237 0.7458 0.7402 0 0.4455 0.5753 0.6524 0.6949 0.7094 0 0.3545 0.4913 0.5816 0.6399 0.6714
Construccin de la tabla de valores para la funcin de Bessel de primera clase, para: Rango de la variable : De 0 a 5, en pasos de 1. Rango de la variable x: De 0 a 2, en pasos de /3.
BESSELJ (0 : 1 : 5, (0 : pi/3: 2*pi)') ans = 1.0000 0 0 0 0 0 0.7441 0.4550 0.1250 0.0223 0.0030 0.0003 0.1698 0.5689 0.3734 0.1443 0.0401 0.0087 -0.3042 0.2846 0.4854 0.3335 0.1514 0.0521 -0.3781 -0.1348 0.3137 0.4344 0.3085 0.1547 -0.0979 -0.3446 -0.0337 0.3188 0.3991 0.2909 0.2203 -0.2124 -0.2879 0.0291 0.3157 0.3728
Construccin de la tabla de valores para la funcin de Bessel de segunda clase, para: Rango de la variable : De 0 a 0,5, en pasos de 0,1. Rango de la variable x: De 0 a 1, en pasos de 0,2.
BESSELY (0 : 1 : 5, (0 : pi/3: 2*pi)') ans = -Inf -Inf -Inf -Inf 0.1242 -0.7411 -1.5396 -5.1396 0.5180 -0.0547 -0.5703 -1.0344 0.3284 0.3589 -0.0999 -0.4861 -0.0896 0.3700 0.2663 -0.1157 -0.3339 0.0667 0.3594 0.2078 -0.2291 -0.2391 0.1530 0.3365
-Inf -Inf -27.9080 -208.0617 -2.3930 -8.1063 -0.8284 -1.6235 -0.4321 -0.7095 -0.1212 -0.3930 0.1683 -0.1222
5.57
Alternancia de los ceros de las funciones de Bessel de ndice consecutivo:

% Comprobar en forma grfica la alternancia entre los ceros % de las funciones J3(x) y J2(x) % Rango de x: De 0 a 20 en pasos de 0.2: J2x = besselj(2,(0:0.2:20)); 0.5 J3x = besselj(3,(0:0.2:20)); x = [0 : 0.2 : 20]; 0.4 J3(x) plot(x,J2x,x,J3x)
B B
0.3 0.2 0.1 0 -0.1 -0.2 -0.3 -0.4
J2(x)
B B
10
12
14
16
18
20
Propiedades de las funciones de Besssel:

% Verificar la propiedad J-n(x) = (-1)^n*Jn(x). % Si n es par, entonces: J-n(x) = Jn(x). Sea por ejemplo n = 2: syms x J_2x = besselj(-2,x) J_2x = besselj(2,x) % Por el contrario, si n es impar: J_3x = besselJ (-3,x) J_3x = - besselj(3,x)
Verificar con la ayuda de su representacin grfica, la frmula: d/dx (1/x^n*Jn(x) = - 1/x^n*Jn+1(x): syms x % Sea por ejemplo n = 2 x^(-2)*besselj(2,x) ans = 1/x^2*besselj(2,x) F1 = diff (1/x^2*besselj(2,x)) F1= - 2/x^3*besselj(2,x)+1/x^2*(besselj(1,x) -2/x*besselj(2,x)) F2 = 1/x^2*besselj(3,x) F2 = 1/x^2*besselj(3,x)
5.58
% Si representamos ambas funciones, verificamos que las mismas son iguales y opuestas: ezplot (F1) hold on 1/x^2*besselj(3,x) 0.04 ezplot (F2) grid 0.03
0.02 0.01 0 -0.01 -0.02 -0.03 -0.04 -5 -4 -3 -2 -1 0 x 1 2 3 4 5
Races de Jo(x): Con ayuda de la grfica de la funcin, hallar por aproximaciones sucesivas las dos primeras races y calcular la separacin entre ellas:
% Empecemos por representar la funcin 1 syms x J0 = besselj(0,x); ezplot (J0) 0.5 % En el grfico vemos que las races caen en % x1 = 2,4 y x2 = 5,5 aproximadamente. besselj(0,(5.52 : 0.001 : 5.525)) 0 ans = -0.0000 0.0003 0.0007 0.0010 0.0013 0.0017 % Una raz cae en x = 5,52. La otra en: -0.5 -6 -4 -2 0 2 besselj(0,(2.40 : 0.001 : 2.405)) x ans = 0.0025 0.0020 0.0015 0.0009 0.0004 -0.0001 % Es decir, en 2.405 aprox. La diferencia entre ambas es: % Distancia entre ceros: 5.52 - 2.405 ans = 3.1150 % Vemos nuevamente que la diferencia entre races es una cantidad prxima a pi.
besselj(0,x)
Resolver la ecuacin del problema 5.13.1.a syms y t Y = dsolve('t^2*D2y+t*Dy+(t^2-.25)*y = 0') Y = (C1*sin(t)+C2*cos(t))/t^(1/2) % Observar que esta ecuacin coincide con la solucin del problema, % con la condicin de asignar a C1 y C2 los valores adecuados.
5.59
Solucin de ecuaciones de Bessel o transformables Bessel: % Resolver la ecuacin de Bessel: y " *x^2 + y ' *x + (x^2 - n^2)*y = 0 % Ejemplo 1: n = 0 syms y t Y = dsolve('t^2*D2y+t*Dy+t^2*y=0') Y = C1*besselj(0,t)+C2*bessely(0,t) % Ejemplo 2: n = 2 Y = dsolve('t^2*D2y+t*Dy+(t^2-4)*y=0') Y = C1*besselj(2,t)+C2*bessely(2,t) % Ejemplo 3: n = 5 Y = dsolve('t^2*D2y+t*Dy+(t^2-25)*y=0') Y = C1*besselj(5,t)+C2*bessely(5,t) % etc. syms y t % Resolver la ecuacin y " *x^2 + 6*y ' *x + x^2 *y = 0 Y = dsolve('t^2*D2y+6*t*Dy+t^2*y=0')
Y = (C1*cos(t)*t^2-3*C1*cos(t)-3*C1*sin(t)*t+C2*sin(t)*t^2-3*C2*sin(t) +3*C2*cos(t)*t)/t^5
% Resolver la ecuacin y " *x + y ' + x^2 * y/4 = 0 Y = dsolve('t*D2y + Dy + t^2*y/4 = 0') Y = C1*bessely(0,1/3*t^(3/2))+C2*besselj(0,1/3*t^(3/2)) % Resolver la ecuacin y " *x^2 + 3/2*y ' *x + x^2 *y = 0 Y = dsolve('t^2*D2y+3/2*t*Dy+t^2*y=0') Y = (C1*besselj(1/4,t)+C2*bessely(1/4,t))/t^(1/4) % Resolver la ecuacin: y " * x^2 + 6*y ' * x + x^2 * y = 0 Y = dsolve('t^2*D2y+5*t*Dy+2*t^2*y=0') Y = (C1*besselj(2,2^(1/2)*t)+C2*bessely(2,2^(1/2)*t))/t^2
5.60
Verificacin de la solucin( Grfico):

syms y t Y = dsolve('t^2*D2y+t*Dy+t^2*y=0') Y = C1*besselj(0,t)+C2*bessely(0,t) % Hagamos C1 = 1 y C2 = 1: y = besselj(0,t)+bessely(0,t) y = besselj(0,t)+bessely(0,t) d1 = diff(y,t) d1 = -besselj(1,t)-bessely(1,t) d2 = diff(d1,t) d2 = -besselj(0,t)+1/t*besselj(1,t)-bessely(0,t)+1/t*bessely(1,t) F = t^2*d2+t*d1+t^2*y F = t^2*(-besselj(0,t)+1/t*besselj(1,t)-bessely(0,t)+1/t*bessely(1,t))+ + t*(-besselj(1,t)-bessely(1,t))+t^2*(besselj(0,t)+bessely(0,t)) ezplot(F) axis([0 6 .1 .1])
t^2*(-besselj(0,t)+1/t ~~~ elj(0,t)+bessely(0,t)) 0.1 0.08 0.06 0.04 0.02 0 -0.02 -0.04 -0.06 -0.08 -0.1
3 t

5 - Funciones de Bessel

Cargado por

Copyright:

Formatos disponibles

5 - Funciones de Bessel

Cargado por

Información del documento

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

5 - Funciones de Bessel

Cargado por

Copyright:

Formatos disponibles

E cu a ci o n e s D i f e r e n ci a l e s d e O r d en Su p er i o r

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel.

La ecuacin x y" + x y' + ( x - ) y = 0

entonces: y' = A y1' + B y2'

y" = A y1" + B y2"

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel

Las soluciones de esta ecuacin algebraica de segundo grado son: 1 = +

Reemplazando en ella la primera de las dos soluciones posibles: = 1 = +

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel.

Nota: Incluso si hubiramos aceptado la segunda solucin: = 2 = -

Finalmente, para k = 2, 3, 4, . . . , obtenemos la forma general: ( k + ) ( k + - 1 ) ck + ( k + ) ck + ck -2 - 2 ck = 0

Como k es un nmero par, podemos hacer: k = 2 m, con m = 0 , 1 , 2, 3, . . .

con lo cual, las frmulas anteriores se modifican como sigue: c2m =

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel

De modo similar, para c4 ( m = 2 ) ser: c4 =

Funcin de Bessel de 1 Clase

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel.

Reemplazando en la (5.8) resulta:

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel

Empezaremos por recordar el desarrollo en serie de Laurent de una funcin f ( z ), en zo = (0,0):

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel.

Los coeficientes an y bn estn dados por las frmulas:

Para nuestro caso vamos a llamar: an = n

As, la serie de Laurent que define la funcin f ( z ), puede escribirse:

Demostraremos seguidamente que si la funcin f ( z ) es f() = e siendo = ei

y reemplazando en la integral, obtenemos:

Funcin de Bessel de primera clase y orden n

R. Abascal - Anlisis de Seales y Sistemas

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel

Por otro lado, como e i (x sen

= cos ( x sen - n ) + i sen ( x sen - n )

la funcin de Bessel se puede escribir tambin de la manera siguiente:

sen ( x sen ) cos n d

cos ( x sen ) sen n d

Llamemos: f1 ( ) = sen ( x sen )

sen ( x sen - n ) d = f1 ( ) cos n d

Si comparamos la frmula de los coeficientes an en el desarrollo de una funcin en serie de Fourier,

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel.

cos ( x sen ) sen n d =

En resumen, la (5.11) queda reducida a la igualdad:

Expresin Integral de la Funcin de Bessel de 1 Clase

Siendo el integrando una funcin coseno, de perodo 2 , puede ponerse, alternativamente:

de lo cual se deducen las relaciones siguientes: = cos + i sen , n = cos n + i sen n

Restando entre s las dos primeras se advierte que: - -1 = 2 i sen

Si elevamos ambos miembros a la potencia x, obtenemos: e

Al desarrollar las dos exponenciales en serie, encontramos una nueva igualdad:

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel

2 x4- 1 x6 + ... 1- x + 1 4 (2!)2 16 (3!)2 64

Recordemos la frmula (5.7) de la funcin Bessel de 1 clase:

Para n = 1: (-1)0 x0 + x 1 (-1) x2 + x 1 (-1)2 x4 + x 1 (-1)3 x6 + . . . J1 ( x ) = x 1 2 0! 1! 23 1! 2! 25 2! 3! 27 3! 4!

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel.

+ -1 J-1 ( x ) + -2 J-2 ( x ) + -3 J-3 ( x ) + . . .

Vimos por otra parte que si n es un entero, entonces: J-n ( x ) = ( - 1 ) n . Jn ( x ),

Sumndolas o restndolas entre s, deducimos respectivamente: n + -n = 2 cos n

Entonces, reemplazando en la (5.16), hallamos: ei x sen = Jo ( x ) + 2 [ J2 ( x ) cos 2 + J4 ( x ) cos 4 + . . . ] +

+ 2 i [ J1 ( x ) sen + J3 ( x ) sen 3 + J5 ( x ) sen 5 + . . . ]

Pero tambin, como ei x sen = cos ( x sen ) + i sen ( x sen )

igualando las partes reales e imaginarias, se obtiene:

R. Abascal - Anlisis de Seales y Sistemas

Ecuaciones Diferenciales de Orden Superior. Parte 5 Funciones de Bessel