Fis-Mat

 

A Equação de Bessel

bessel01.GIF (603 bytes)

2. Solução da Equação (1) em série de potências

 

 

//////Seja bessel02b.GIF (343 bytes)uma possível solução de (1). Ao admitirmos esta possibilidade, na realidade o que estamos fazendo é assumir y, dado pela equação acima, como solução de (1), e propondo-nos a determinar os coeficientes an, bem como gama.GIF (62 bytes). Nestas condições temos

bessel02c.GIF (491 bytes)

e

bessel02d.GIF (643 bytes)

//////Substituindo estes valores de y, y’ e y" em (1) obtemos:

bessel02e.GIF (1077 bytes)

ou

bessel02f.GIF (826 bytes)

ou, ainda,

bessel02g.GIF (1344 bytes)

//////A equação (3) pode ser desmembrada em três igualdades:

bessel0456.GIF (1287 bytes)

//////As igualdades (4) e (6) são chamadas, respectivamente, equação indicial e equação de recorrência.

//////Da equação indicial, e assumindo ao dif.GIF (56 bytes) 0, vem

gama.GIF (62 bytes) = maismenos.GIF (61 bytes) ni.GIF (55 bytes)

o que implica, de (5), em:

a1 = 0.

A equação de recorrência (6) pode ser escrita na forma:

bessel06a.GIF (401 bytes)

e como a1 = 0, segue-se que todos os coeficientes an com índice ímpar são nulos: a3 = a5 = a7 = ... = 0. Por outro lado, os coeficientes an com índice par são dados por:

bessel06b.GIF (1607 bytes)

Generalizando, temos:

bessel06c.GIF (871 bytes)

(n par)

 

//////Efetuando a mudança de variável n setaD.GIF (62 bytes) m tal que n = 2m, os coeficientes an podem ser expressos por:

bessel06d.GIF (926 bytes)

ou

bessel06e.GIF (921 bytes)

ou, finalmente:

bessel07.GIF (938 bytes)

 

/
/
//////
Lembrando que bessel02b.GIF (343 bytes), n = 2m e gama.GIF (62 bytes) = maismenos.GIF (61 bytes) ni.GIF (55 bytes), temos:

bessel08.GIF (1293 bytes)

/
//////
Substituindo (7) em (8) as soluções y1(x) e y2(x) ficam:
/

bessel09.GIF (2320 bytes)

//////As expressões (9) podem ainda ser escritas, em termos de funções gama (gamaM.GIF (65 bytes)) [4], da seguinte maneira:

bessel10.GIF (2233 bytes)

//////Determinadas duas soluções y1 e y2 para a equação (1), podemos, uma vez comprovada a existência destas soluções, bem como sua independência linear, escrever a solução geral como

y(x) = C1y1(x) + C2y2(x).

(11)

No item a seguir veremos que (11), com y1 e y2 dados por (10), não é solução geral válida para ni.GIF (55 bytes) = 0, 1, 2, ...

Notas e Referências:

4. A função gama goza da seguinte propriedade: gamaM.GIF (65 bytes)(z+1) = zgamaM.GIF (65 bytes)(z). Portanto: gamaM.GIF (65 bytes)(m+p+1) = (m+p)gamaM.GIF (65 bytes)(m+p) = (m+p)(m-1+p)gamaM.GIF (65 bytes)(m-1+p) = (m+p)(m-1+p)(m-2+p)gamaM.GIF (65 bytes)(m-2+p) = (m+p)(m-1+p).........(1+p)gamaM.GIF (65 bytes)(1+p). Então, (m+p)(m-1+p)......(1+p) = FGama01.GIF (308 bytes). Analogamente, mostra-se que (m-p)(m-1-p).......(1-p) = FGama02.GIF (306 bytes). Voltar