Demostracion Practica Bessel 1

Captulo 10
Otras ecuaciones relacionadas con la

de Bessel
10.1.
Funciones de Bessel hiperb

olicas
La ecuacion
d2y
dy
+ x (x2 + 2 )y = 0
(10.1)
2
dx
dx
es parecida a la de Bessel, pero tiene un signo cambiado. Se le llama ecuacion de Bessel
hiperbolica (o modificada). Podemos empezar de cero el proceso de Frobenius para
encontrar la soluci
on o, alternativamente, podemos hacer un valiente cambio de variable
t = ix, siendo i = 1, la unidad imaginaria. Esto lleva facilmente a

x2
2d
2y
dt 2
+t
dy
+ (t 2 2 )y = 0
dt
(10.2)
que es de Bessel. La solucion es, por tanto, y(t) = c1 J (t) + c2Y (t) o y(t) = c1 J (t) +
c2 J (t) si #= 0, 1, 2, 3 . . . .
10.1.1. Si reemplazamos t = ix y escribimos la serie de J (ix) a partir de su

definici
on obtenemos J (ix) = i I (x) siendo
! x "2m+
1
I (x) =
m=0 m!(m + 1 + ) 2
(10.3)
una serie real.

A I (x) se le conoce como funcion de Bessel hiperbolica (o modificada) de primera especie de orden y es una solucion (real) de la ecuacion (10.1). Si #= 0, 1, 2, . . . la otra
solucion es, sencillamente, I (x). Si = n = 0, 1, 2, 3, . . . es facil ver que In (x) = In (x),
por lo que tenemos que encontrar todava la segunda solucion. Otra vez podramos hacer todo el proceso que hicimos en la ecuacion de Bessel, pero es mas sencillo partir
de las soluciones de la funcion de Bessel Jn (t), Yn (t) con t = ix. La solucion es entonces y(x) = c1 Jn (ix) + c2Yn (ix). Aqu conviene separar la parte real de la parte imaginaria
72
(1)
dandose cuenta de que la llamada 1a funcion de Hankel1 H (z) = J (z) + iY (z) cumple
(1)
(1)
que i+1 H (ix) es real, lo que lleva a definir la funcion real K (x) = 2 i+1 H (ix), la
funcion de Bessel hiperbolica (o modificada) de segunda especie, que es una segunda
solucion de la ecuacion de Bessel hiperbolica (o modificada) valida tambien en el caso
= n = 0, 1, 2, 3, . . . .
10.1.2. Una definicion alternativa es:

K (x) =
I (x) I (x)
2
sin()
con el correspondiente lmite lmn cuando sea necesario.

La solucion de la ecuacion de Bessel hiperbolica (o modificada) es pues, y(x) =
c1 I (x) + c2 I (x) si #= 0, 1, 2, . . . , o y(x) = c1 I (x) + c2 K (x) valida siempre y, en
particular, cuando = n = 0, 1, 2, 3 . . . . Estas

funciones son un poco mas aburridas que
las de Bessel:
In[179]:=
Plot@8BesselI@0, xD, BesselI@2, xD, BesselI@4, xD<, 8x, -5, 5<D

15
10
Out[179]=
-4
In[180]:=
-2
Plot@8BesselI@1, xD, BesselI@3, xD, BesselI@5, xD<, 8x, -5, 5<D

In[184]:=
2.0
Out[180]=
-4
-2
Plot@8BesselK@0, xD, BesselK@1, xD, BesselK@2, xD<, 8x, 0, 5<D
1.5
Out[184]= 1.0
-2
0.5
-4
1
(2)
Tambien se define la 2a funcion de Hankel: H (z) = J (z) iY (z).
10.2 Ecuaciones reducibles a la de Bessel
10.2.
73
Ecuaciones reducibles a la de Bessel
Tomemos la ecuacion de Bessel para la funcion w(z):

z2
d2w
dw
+
z
+ (gz2 2 )w = 0
2
dz
dz
Si g = 1 es la ecuacion de Bessel usual; si g = 1 es la ecuacion de Bessel hiperbolica. La solucion para g = 1 es w(z) = c1 J (z) + c2 K (z) o w(z) = c1 J (z) + c2 J (z)
si #= n = 0, 1, 2, 3, . . . . Para g = 1 la solucion es w(z) = c1 I (z) + c2 K (z) o w(z) =
c1 I (z)+ c2 I (z) si #= n = 0, 1, 2, 3, . . . . Hagamos ahora el cambio de variables z = bxc ,
w = y/xa .
10.2.1. La ecuacion se transforma en:

x2
con
dy
d2y
+ x + ( x + )y = 0
2
dx
dx
=
=
=
=
(10.4)
1 2a
gb2 c2
2c
a2 c2 2
Alternativamente, si tenemos la ecuacion (10.4), podemos encontrar a, b, c, g, :

a
g
c
b
=
=
=
=
2 =
(1 )/2
signo( )
#
/2
2 | |/
(1)2 4
4c2
Si (1 )2 4 > 0, entonces el cambio anterior de variables la reduce a una ecuacion

de Bessel (g = 1) o de Bessel hiperbolica (g = 1). La solucion es, por tanto, para g = 1,
y(x) = c1 xa J (bxc ) + c2 xa K (bxc ), etc.
10.2.2. La ecuacion de Airy y%% = xy se puede escribir en la forma (10.4) con
= 0, = 1, = 3, = 0. La soluci
on es y(x) = c1 x1/2 I1/3 ( 23 x3/2 )+c2 x1/2 I1/3 ( 23 x3/2 ).
Otro ejemplo lo ofrece la ecuacion:
x2 y%% + 2xy% + (x2 n(n + 1))y = 0,
(10.5)
con n N. Parece de Bessel, pero no lo es (el problema es el termino 2xy% ). Es del tipo
(10.4), con = 2, = 1, = 2, = n(n + 1).
74
10.2.3.
La soluci
on es y(x) = c1 jn (x) + c2 yn (x), siendo jn (x), yn (x) las funciones
de Bessel esf
ericas:
$
J 1 (x)
2x n+ 2
$
Y 1 (x)
yn (x) =
2x n+ 2
jn (x) =
(10.6)
(10.7)
Estas funciones se pueden escribir en terminos de funciones trigonometricas. Se demostrara por induccion que:
& %
&
1 d n sin x
jn (x) = (1) x
x dx
x
%
&n !
cos x "
n+1 n 1 d
yn (x) = (1) x
x dx
x
n n
(10.8)
(10.9)
En particular:
sin x
j0 (x) =
x
2
sin x
x
$
cos x
2
y0 (x) =
Y1/2 (x) =
cos x
x
x
10.3.
J1/2 (x) =
(10.10)
Comportamiento asint
otico de las funciones de
Bessel
Veremos ahora como se comportan las funciones de Bessel en los lmites de x 0

y x . Para x 0 es muy sencillo, basta con tomar el termino dominante en la serie
de potencias que define a cada funcion. El termino dominante es el de la potencia mas
baja. Si hay un termino ln x este domina sobre las potencias positivas x pero no sobre
10.3 Comportamiento asint

otico de las funciones de Bessel
75
las negativas. Con todo esto es muy sencillo demostrar que:

J (x)
=
x0
Jn (x)
=
x0
Y (x)
=
x0
Y0 (x)
=
x0
I (x)
=
x0
In (x)
=
x0
K (x)
=
x0
K0 (x)
=
x0
jn (x)
=
x0
yn (x)
=
x0
'
(
x
1 + O(x2 ) ,
#= n = 1, 2, 3, . . .
2 ( + 1)
(
(1)n xn '
1 + O(x2 ) ,
n = 1, 2, 3, . . .
n
2 n!
(
2 () '
x
1 + O(x2 ) ,
>0
2
(ln x + ln 2) + O(x2 ),
'
(
x
1 + O(x2 ) ,
#= n = 1, 2, 3, . . .
2 ( + 1)
(
xn '
2
1
+
O(x
)
,
n = 1, 2, 3, . . .
2n n!
(
21 () '
2
1
+
O(x
)
,
>0
x
ln x + ln 2 + O(x2 ),
'
(
xn
1 + O(x2 ) ,
n = 0, 1, 2, . . .
(2n + 1)!!
'
(
(2n 1)!!xn1 1 + O(x2 ) ,
n = 0, 1, 2, . . .
(10.11)
(10.12)
(10.13)
(10.14)
(10.15)
(10.16)
(10.17)
(10.18)
(10.19)
(10.20)
Vamos ahora al comportamiento cuando x . Aqu las cosas no son tan rigurosas,
pero no son demasiado complicadas. Se trata de hacer el cambio de variables y(x) =
v(x)x1/2 para obtener la ecuacion:
%
&
2 1/4
%%
v + 1
v = 0.
(10.21)
x2
Si queremos estudiar el lmite x podemos reducir la ecuacion a v%% + v = 0, de
solucion v(x) = A cos(x + ), siendo A y constantes de integracion. Esto sugiere que
para x la funcion de Bessel se comporta como J (x) = Ax1/2 cos(x + ). Un calculo
auxiliar que utiliza una representacion integral de la funcion de Bessel y la tecnica del
punto de silla de integracion, da algo mas concreto:
$
!
2
"
J (x) =
cos x (2 + 1)
.
(10.22)
x
4
x
Para = 0, por ejemplo, tenemos J0 (x)
cion para x ! 2, como se ve en la figura.
2
x cos
*
+
x 4 , que es suficiente aproxima-
10.3.1. Usando la definicion (9.5) de la funcion de Weber Y (x) se llega a:

Y (x)
=
x
!
2
"
sin x (2 + 1)
.
x
4
(10.23)
76
Otra vez, una aproximacion razonable si x ! 2.

Por ultimo, mencionaremos los comportamientos asintoticos de K e I , validos para
cualquier valor de :
PlotB:BesselJ@0, xD,
In[189]:=
I (x)
=
x
K (x)
=
x
2
px
CosBx -
p
4
(
x'
e 1 + O(x2 ) ,
2x
$
(
2 x '
e 1 + O(x2 )
x
F>, 8x, 0, 12<F
In[190]:=
(10.24)
(10.25)
PlotB:BesselY@0, xD,
2
px
SinBx -
p
4
F>, 8x, 0, 12<F
0.5
1.0
0.8
0.6
0.4
Out[189]=
10
12
Out[190]=
0.2
-0.5
2
10
12
-0.2
-1.0
-0.4
In[248]:=
PlotB:BesselI@0, xD,
1
2px
In[250]:=
p
2x
Exp@- xD>, 8x, 0, 6<F
1.0
0.8
50
0.6
40
Out[248]=
PlotB:BesselK@0, xD,
Exp@xD>, 8x, 0, 6<F
Out[250]=
30
0.4
20
0.2
10
1
1

Demostracion Practica Bessel 1

Cargado por

Copyright:

Formatos disponibles

Demostracion Practica Bessel 1

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

Demostracion Practica Bessel 1

Cargado por

Copyright:

Formatos disponibles

Captulo 10

Otras ecuaciones relacionadas con la

Funciones de Bessel hiperb

t = ix, siendo i = 1, la unidad imaginaria. Esto lleva facilmente a

10.1.1. Si reemplazamos t = ix y escribimos la serie de J (ix) a partir de su

una serie real.

10.1.2. Una definicion alternativa es:

con el correspondiente lmite lmn cuando sea necesario.

particular, cuando = n = 0, 1, 2, 3 . . . . Estas

Plot@8BesselI@0, xD, BesselI@2, xD, BesselI@4, xD<, 8x, -5, 5<D

Plot@8BesselI@1, xD, BesselI@3, xD, BesselI@5, xD<, 8x, -5, 5<D

Plot@8BesselK@0, xD, BesselK@1, xD, BesselK@2, xD<, 8x, 0, 5<D

Tambien se define la 2a funcion de Hankel: H (z) = J (z) iY (z).

10.2 Ecuaciones reducibles a la de Bessel

Ecuaciones reducibles a la de Bessel

Tomemos la ecuacion de Bessel para la funcion w(z):

10.2.1. La ecuacion se transforma en:

Alternativamente, si tenemos la ecuacion (10.4), podemos encontrar a, b, c, g, :

Si (1 )2 4 > 0, entonces el cambio anterior de variables la reduce a una ecuacion

10.2.2. La ecuacion de Airy y%% = xy se puede escribir en la forma (10.4) con

Veremos ahora como se comportan las funciones de Bessel en los lmites de x 0

10.3 Comportamiento asint

las negativas. Con todo esto es muy sencillo demostrar que:

10.3.1. Usando la definicion (9.5) de la funcion de Weber Y (x) se llega a:

Otra vez, una aproximacion razonable si x ! 2.

F>, 8x, 0, 12<F

F>, 8x, 0, 12<F

Exp@- xD>, 8x, 0, 6<F

Exp@xD>, 8x, 0, 6<F

También podría gustarte