Fis-Mat

 

A Equação de Bessel

bessel01.GIF (603 bytes)

7. Ortogonalidade das Funções de Bessel

 

 

Sejam Sigma.GIF (54 bytes) e mi.GIF (59 bytes) dois autovalores da equação

x2y" + xy’ + (lambda.GIF (59 bytes)2x2 - ni.GIF (55 bytes)2)y = 0

(17)

Nestas condições, e conforme visto no item 4,

y1 = JNi2.GIF (55 bytes) (Sigma.GIF (54 bytes)x)          e          y2 = JNi2.GIF (55 bytes) (mi.GIF (59 bytes)x)

(18)

são soluções da equação (17) para os correspondentes autovalores (lambda.GIF (59 bytes) = Sigma.GIF (833 bytes) e lambda.GIF (59 bytes) = mi.GIF (59 bytes), respectivamente). Podemos então escrever:

x2y1" + xy1' + (Sigma.GIF (833 bytes)2x2 - ni.GIF (55 bytes)2)y1 = 0

(19)

x2y2" + xy2' + (mi.GIF (59 bytes)2x2 - ni.GIF (55 bytes)2)y2 = 0

(20)

Multiplicando (19) por y2 e (20) por y1 e subtraindo membro a membro as equações (20) da (19), assim modificadas, temos:

x2(y2y1"-y1y2") + x(y2y1'-y1y2') = (mi.GIF (59 bytes)2 - Sigma.GIF (833 bytes)2)x2y1y2

Dividindo por x e manipulando os termos, podemos escrever:

bessel20a.GIF (1085 bytes)

ou

bessel21.GIF (823 bytes)

(21)

Integrando (21) no intervalo [a,b]:

Bessel21a.GIF (850 bytes).

Dividindo por (mi.GIF (59 bytes)2 - Sigma.GIF (833 bytes)2) e utilizando as igualdades (18):

Bessel21b.GIF (1743 bytes)

Para a=0 e b=R, temos:

bessel22.GIF (1680 bytes)

(22)

e para mi.GIF (59 bytes) setaD.GIF (62 bytes) Sigma.GIF (833 bytes) a expressão fica:

bessel22a.GIF (1633 bytes)

Aplicando L’Hospital e a seguir multiplicando o numerador e o denominador do segundo membro da expressão por Sigma.GIF (833 bytes):

bessel23.GIF (3401 bytes)

(23)

Mas JNi2.GIF (55 bytes)(Sigma.GIF (833 bytes)R) é solução da equação de Bessel

Sigma.GIF (833 bytes)2R2y"(R) + Sigma.GIF (833 bytes)Ry’(R) + (Sigma.GIF (833 bytes)2R2 - ni.GIF (55 bytes)2)y(R) = 0,

ou seja

Sigma.GIF (833 bytes)2R2JNi2.GIF (55 bytes)"(Sigma.GIF (833 bytes)R) + Sigma.GIF (833 bytes)RJNi2.GIF (55 bytes)(Sigma.GIF (833 bytes)R) = - (Sigma.GIF (833 bytes)2R2 - ni.GIF (55 bytes)2)JNi2.GIF (55 bytes)(Sigma.GIF (833 bytes)R).

(24)

Substituindo (24) em (23):

bessel24a.GIF (1508 bytes)

ou

bessel25.GIF (1587 bytes)

(25)

Se Sigma.GIF (833 bytes), mi.GIF (59 bytes) e R forem tais que as raízes n-ésima e m-ésima da função de Bessel de ordem ni.GIF (55 bytes) igualem, respectivamente a alfaniene.GIF (101 bytes) = Sigma.GIF (833 bytes)R e alfanieme.GIF (110 bytes) = mi.GIF (59 bytes)R, as expressões (22) e (25) podem ser agrupadas numa única expressão:

bessel26.GIF (1216 bytes) *

(26)

conhecida como relação de ortogonalidade das funções de Bessel. De (26) observa-se que as funções Raizxis.GIF (142 bytes)bessel26a.GIF (369 bytes) e Raizxis.GIF (142 bytes)bessel26b.GIF (381 bytes) são ortogonais no intervalo [0,R]; ou então, as funções bessel26a.GIF (369 bytes) e bessel26b.GIF (381 bytes) são ortogonais em relação à função peso x.

 

* deltami.GIF (71 bytes)n,m = 0 se n dif.GIF (56 bytes) m, e deltami.GIF (71 bytes)n,m = 1 se n = m. Observar que se n dif.GIF (56 bytes) m, a equação válida é a (22), e se n = m (ou mi.GIF (59 bytes) setaD.GIF (62 bytes) Sigma.GIF (833 bytes) ) recaímos na equação (25). Como estamos trabalhando com raízes, JNi2.GIF (55 bytes)(Sigma.GIF (833 bytes)R) = JNi2.GIF (55 bytes)(mi.GIF (59 bytes)R) = 0 em ambos os casos, o que nos leva facilmente à generalização expressa em (26). Voltar