Funció de Bessel

De la Viquipèdia, l'enciclopèdia lliure
(S'ha redirigit des de: Funcions de Bessel)
La part radial dels modes de vibració d'un tambor circular segueixen la funció de Bessel.

Les funcions de Bessel són les solucions canòniques de l'equació diferencial de Bessel:[1]

que tenen com a punt singular regular i una singularitat essencial a . El paràmetre és un nombre donat que es pot considerar positiu sense cap pèrdua de generalitat.

Aplicacions de la funció de Bessel[modifica]

L'equació de Bessel sorgeix quan es busquen solucions separables a l'equació de Laplace i l'equació de Helmholtz en coordenades cilíndriques o esfèriques. Les funcions de Bessel són, per tant, especialment importants per molts problemes de propagació d'ones i potencials estàtics. En la resolució de problemes en sistemes de coordenades cilíndriques, s'obtenen funcions de Bessel d'ordre enter (α = n); en problemes esfèrics, se n'obtenen d'ordre semi-enter (α = n+1/2). Per exemple:

  • Ones electromagnètiques en una guia d'ones cilíndriques
  • Amplituds de pressió dels fluxos rotatacionals no viscosos
  • La conducció tèrmica en un objecte cilíndric
  • Modes de vibració d'una membrana acústica prima i circular (o annular) (com un tambor o un altre membranòfon)
  • Problemes de difusió en enreixats
  • Solucions de l'equació de Schrödinger radial (en coordenades esfèriques i cilíndriques) per una partícula lliure
  • Solució de patrons de radiació acústica
  • Fricció depenent de la freqüència en canoncades circulars
  • Dinàmica de cossos flotants
  • Resolució angular

Les funcions de Bessel també apareixen en altres problemes, com en el processament de senyal (per exemple, vegeu finestra de Kaiser o filtre de Bessel).

Definicions[modifica]

Com que es tracta d'una equació diferencial de segon ordre, hi ha d'haver dues solucions linealment independents. Depenent de les circumstàncies, tanmateix, són convenients formulacions diverses d'aquestes solucions. En la següent taula, es resumeixen les diferents variacions, i es descriuen en les seccions següents.

Tipus Primer tipus Segon tipus
Funcions de Bessel Jα Yα
Funcions de Bessel modificades Iα Kα
Funcions de Hankel Hα(1) = Jα + iYα Hα(2) = Jα - iYα
Funcions de Bessel esfèriques jn yn
Funcions de Hankel esfèriques hn(1) = jn + iyn hn(2) = jn - iyn

Les funcions de Bessel de segon tipus i les funcions de Bessel esfèriques de segon tipus sovint s'anoten com a Nn i nn, respectivament, més que no tant Yn i yn.[2][3]

Funcions de Bessel de primer tipus: Jα[modifica]

Les funcions de Bessel de primer tipus, denotades com a Jα(x), són solucions de l'equació diferencial de Bessel que són finites a l'origen (x = 0) per valors d'α enters o positius, i divergeixen a mesura que x tendeix a 0 per valors d' α no enters. Es pot definir la funció a través de la seva expansió en sèrie al voltant de x = 0, que es pot trobar aplicant el mètode de Frobenius a l'equació de Bessel:[4]

on Γ(z) és la funció gamma, una generalització canviada de la funció factorial per valors no enters. La funció de Bessel de primer tipus és una funció entera si α és un enter, altrament és una funció multivaluada amb singularitat en el zero. La representació gràfica de les funcions de Bessel s'assemblen força funcions sinusoidals o cosinusoidals que decauen proporcionalment a 1/√x (vegin-se també les seves formes asimptòtiques més endavant), tot i que les seves arrels no són generalment periòdiques, excepte asimptòticament amb valors grans de x. (Les sèries mostren que −J1(x) és la derivada de J0(x), tal com −sin(x) és la derivada de cos(x); més generalment, la derivada de Jn(x) es pot expressar en termes de Jn±1(x) per les identitats que es mostren més endavant.)

Representació gràfica de la funció de Bessel de primer tipus, J α(x), per ordres d'enters α = 0, 1, 2

Per valors de α no enters, les funcions J α(x) i Jα(x) són linealment dependents, i són per tant les dues solucions de l'equació diferencial. D'altra banda, per α d'ordre enter, la següent relació és vàlida (noti's que la fucnió gamma té pols simples per cadascun dels enters no positius):[5]

Això significa que les dues solucions ja no són linealment independents. En aquest cas, la segona de les solucions linealment independent a la primera es troba com la funció de Bessel de segon tipus, com es comentarà més endavant.

Integrals de Bessel[modifica]

Una altra definició de la funció de Bessel, per valors enters de n, és possible mitjançant l'ús de la representació integral:[6]

Una altra representació integral és:[6]

Aquesta és l'aproximació que va fer Bessel, i a partir d'aquesta definció va derivar diverses propietats de la funció. La definició es pot estendre a ordres no enteres per una de les integrals de Schläfli, per :[6]

[7][8][9][10]

Relació amb les sèries hipergeomètriques[modifica]

Les funcions de Bessel es poden expressa en termes de les sèries hipergeomètriques generalitzades com:[11]

Aquesta expressió està relacionada al desenvolupament de les funcions de Bessel en termes de la funció de Bessel-Clifford.

Relació amb els polinomis de Laguerre[modifica]

En termes dels polinomis de Laguerre Lk i el paràmetre elegit arbitràriament t, la funció de Bessel es pot expressar com:[12]

Funcions de Bessel de segon tipus: Yα[modifica]

Les funcions de Bessel de segon tipus, anotades com a Yα(x), i ocasionalment anotades com a Nα(x), són solucions de l'equació diferencial de Bessel que té una singularitat a l'origen (x = 0) i són multivaluades. De vegades se les anomena funcions de Weber, ja que van ser introduïdes per H. M. Weber (1873), o també funcions de Neumann en honor de Carl Neumann.[13]

Representació gràfica de la funció de Bessel de segon tipus, Yα(x), per ordres enters α = 0, 1, 2.

Per valors d'α no enters, Yα(x) està relacionada amb Jα(x) com:

En el cas d'ordre enter n, la funció es defineix prenent el límit quan α no enter tendeix a n,

Hi ha també una fórmula integral corresponent (per Re(x) > 0),[14]

Yα(x) és necessàriament una segona solució linealment indepdendent de l'equació de Bessel quan α és un enter. Però Yα(x) pot adoptar altres significats més enllà d'aquest. Es pot considerar un company 'natural' de Jα(x). Vegeu també el subapartat de les funcions de Hankel més endavant.

Quan α és un enter, a més, com ho era similarment en el cas de funcions del primer tipus, la següent relació és vàlida:

Tant J α(x) com Y α(x) són funcions holomorfes de x en el pla complex tallat al llarg de l'eix real negatiu. Quan α és un enter, les funcions de Bessel J són funcions enteres d'x. Si x es manté constant en un valor no zero, llavors les funcions de Bessel són funcions enteres d' α.

Les funcions de Bessel del segon tipus quan α és un enter és un exemple del segon tipus de solució en el teorema de Fuchs.

Funcions de Hankel:Hα(1), Hα(2)[modifica]

Una altra formulació important de les dues solucions linealment independents a l'equació de Bessel són les funció Hankel de primer i segon tipus, H(1)
α
(x)
i H(2)
α
(x)
, definides com:[15][16]

on i és la unitat imaginària. Aquestes combinacions lineals també són conegudes com a funcions de Bessel de tercer tipus; són dues solucions linealment independents de l'equació diferencial de Bessel. Duen el nom de Hermann Hankel.

La importància de les funcions de Hankel de primer i segon tipus rau més en llur desenvolupament teòric més que no en la seva aplicació. Aquestes formes de combinació lineal satisfan diverses propietats d'aspecte simple, com les fórmules asimptòtiques o les representacions integrals. Aquí, "simple" significa una aparença del factor de forma eif(x). Es pot entendre, doncs, que la funció de Bessel de segon tipus apareix de manera natural com la part imaginària de les funcions de Hankel.

Les funcions de Hankel de primer i segon tipus són usades per representar les solucions d'ones entrants i sortints d'una equació d'ones en simetries cilíndriques respectivament (o viceversa depenent de la convecció de signe de la freqüència).

Usant les relacions prèvies, es poden expressar com:

Si α és un enter, s'ha de calcular el límit the limit. Les següents relacions són vàlides, tant si α és enter com si no ho és:[17]

En particular, si α = m + 1/2 amb m enter no negatiu, les relacions superiors impliquen directament que:

Això és útil en el desenvolupament de les funcions esfèriques de Bessel (veure més a baix).

Les funcions de Hankel admeten les següents representacions integrals per Re(x) > 0:[18]

on els límits de les integrals indiquen que la integració és al llarg d'un contorn que es pot triar segons: de −∞ a 0 al llarg de l'eix real negatiu, de 0 a ±iπ al llarg de l'eix imaginari, i de ±iπ a +∞ ± iπ al llarg d'un contorn paral·lel a l'eix real.[19]

Funcions de Bessel modificades : Iα, Kα[modifica]

Les funcions de Bessel són vàlides fins i tot amb arguments complexos x, i un cas especialment particular és quan l'argument és imaginari pur, és a dir, quan la part real de l'argument és zero. En aquest cas, les solucions de l'equació de Bessel reben el nom funcions de Bessel modificades (o també funcions hiperbòliques de Bessel) de primer i segon tipus i estan definides com:[20]

quan α no és un enter; quan α és un enter, llavors s'utilitza el límit. Es tria que tinguin valors reals per arguments reals i positius x. L'expansió en sèrie per Iα(x) és doncs similar a la de Jα(x), però sense el factor alternant (−1)m.

Si −π < arg xπ/2, Kα(x) poden ser expressades com a funcions de Hankel de primer tipus:

i si π/2 < arg x ≤ π, es pot expressar com una funció de Hankel de segon tipus:

Es poden expressar les funcions primera i segona de Bessel en termes de les funcions modificades de Bessel (són vàlides si −π < arg zπ/2):[21]

Iα(x) i Kα(x) són les dues solucions linealment independents de l'equació modificada de Bessel:[22]

A diferència del cas de les funcions ordinàries de Bessel, que oscil·len com a funció d'un argument real, Iα i Kα són funcions exponencialment creixents i decreixents respectivament. Com en la funció de Bessel ordinària Jα, la funció Iα val zero a x = 0 per α > 0 i és finita a x = 0 per α = 0. De forma anàloga, Kα divergeix a x = 0 amb una singularitat de tipus logarítmica.[23]

Funcions de Bessel modificades del primer tipus, Iα(x), per α = 0, 1, 2, 3
Funcions de Bessel modificades de segon tipus, Kα(x), per α = 0, 1, 2, 3


Dues fórmules integrals per les funcions de Bessel modificades són (per Re(x) > 0):[24]

En alguns càlculs de física, pot ser útil saber el valor de la següent relació:

Es pot demostrar mostrant la igualtat amb la definició integral de més amunt per K0. Això es fa integrant una corba tancada en el primer quadrant del pla complex.

Es poden representar les funcions de Bessel modificades K1/3 i K2/3 en termes de les integrals ràpidament convergents:[25]

Les funcions de Bessel modificades de segon tipus han rebut també els següents noms (actualment en desús):

  • Funció de Basset function per Alfred Barnard Basset
  • Funció de Bessel modificada de tercer tipus
  • Funció de Hankel modificada[26]
  • Funció de Macdonald per Hector Munro Macdonald

Funcions de Bessel esfèriques: jn, yn[modifica]

Funcions de Bessel esfèriques de primer tipus, jn(x), per n = 0, 1, 2
Funcions de Bessel esfèriques de segon tipus, yn(x), per n = 0, 1, 2

Quan se soluciona l'equació de Helmholtz en coordenades esfèriques amb separació de variables, l'equació radial té la forma:

Les dues solucions linealment independents a aquest equació s'anomenen funcions de Bessel esfèriques jn i yn, i estan relacionades amb les funcions de Bessel ordinàries Jn i Yn segons:[27]

yn també es denota nn o ηn; alguns autors anomenen aquestes funcions the funcions de Neumann esfèriques.

Les funcions de Bessel esfèriques també es poden escriure com a (fórmules de Rayleigh)[28]

La primera funció de Bessel esfèrica j0(x) també es coneix com la funció sinc (no normalitzada). Les primeres funcions de Bessel esfèriques són:[29]

i:[30]

Funcions generadores[modifica]

Les funcions esfèriques de Bessel tenen les funcions generadores:[31]

Relacions diferencials[modifica]

En la següent expressió, fn és qualsevol de jn, yn, h(1)
n
, h(2)
n
per n = 0, ±1, ±2, ...[32]

Funcions de Hankel esfèriques:hn(1), hn(2)[modifica]

Les funcions esfèriques de Hankel es defineixen de forma anàloga a les no esfèriques:

De fet, això implica que existeixen expressions tancades de les funcions de Bessel d'ordre semienter en terme de funcions trigonomètriques i, per tant, també de les funcions esfèriques de Bessel. D'això es dedueix que, per n enter no negatiu es té:

i és la funció complexa conjugada d'aquesta (per a real). D'aquesta fórmula es poden deduir les formes tancades de les funcions esfèriques de Bessel ordinàries, per exemple, i , i en general per a qualsevol argument n.

Funcions de Riccati-Bessel: Sn, Cn, ξn, ζn[modifica]

Les funcions de Riccati-Bessel són una petita modificació de les funció de Bessel esfèriques:

Aquestes funcions satisfan la següent equació diferencial:

Aquesta equació difeencial i les seves solucions, les equacions de Riccati-Bessel, s'usen per resoldre el problema de la dispersió d'ones electromagnètiques per una esfera, un problema conegut com a difusió de Mie, per la publicació per primer cop d'aquests resultat per Mie, l'any 1908. Vegeu, per exemple, Du (2004).[33]

Segons Debye (1909) s'utilitza de vegades la notació en lloc de .

Funció de Bessel per a ν enter[modifica]

El cas especial més habitual i important és aquell en què ν és un nombre enter; llavors ν es coneix com l'ordre de la funció de Bessel. Les funcions de Bessel de primer ordre són les solucions de l'equació diferencial de Bessel que són finites a l'origen per a enters ν no negatius i divergeixen quan x s'aproxima a zero per a ν negatius no enters. La solució general és:

on J es defineix mitjançant la seva sèrie de Taylor al voltant de zero:

i és la funció gamma.

Referències[modifica]

A Wikimedia Commons hi ha contingut multimèdia relatiu a: Funció de Bessel
  1. E. Boyce, William; C. DiPrima, Richard. Elementary Differential Equation and Boundary Value Problem (en anglès). 10a. JohnWiley & Sons, Inc, 2012, p. 296. ISBN 978-0-470-45831-0. 
  2. Weisstein, Eric W. "Spherical Bessel Function of the Second Kind." MathWorld - A Wolfram Web Resource.
  3. Weisstein, Eric W. "Bessel Function of the Second Kind." MathWorld - A Wolfram Web Resource.
  4. Abramowitz and Stegun, p. 360, 9.1.10.
  5. Abramowitz and Stegun, p. 358, 9.1.5.
  6. 6,0 6,1 6,2 Temme, Nico M. Special functions : an introduction to the classical functions of mathematical physics. 2. print.. Nova York [u.a.]: Wiley, 1996, p. 228–231. ISBN 0471113131. 
  7. Watson, p. 176
  8. Gerlach, Ulrich. «Properties of Hankel and Bessel Functions» (en anglès), 2007. Arxivat de l'original el 23 de setembre 2010. [Consulta: 31 gener 2023].
  9. Lund, Mette. «Integral representations of the Bessel function» (en anglès), 2000. Arxivat de l'original el 3 d’octubre 2022. [Consulta: 31 gener 2023].
  10. Arfken & Weber, exercise 11.1.17.
  11. Abramowitz and Stegun, p. 362, 9.1.69.
  12. Szegö, G. Orthogonal Polynomials, 4th ed. Providence, RI: Amer. Math. Soc., 1975.
  13. [enllaç sense format] http://www.mhtlab.uwaterloo.ca/courses/me755/web_chap4.pdf
  14. Watson, p. 178.
  15. Abramowitz and Stegun, p. 358, 9.1.3, 9.1.4.
  16. Greenberg, Michael D. Advanced Engineering Mathematics (en anglès). 2a. Upper Saddle River, Nova Jersey: Prentice Hall, 1998, p. 236. ISBN 0-13-321431-1. 
  17. Abramowitz and Stegun, p. 358, 9.1.6.
  18. Abramowitz and Stegun, p. 360, 9.1.25.
  19. Watson, p. 178.
  20. Abramowitz and Stegun, p. 375, 9.6.2, 9.6.10, 9.6.11.
  21. Abramowitz and Stegun, p. 375, 9.6.3, 9.6.5.
  22. Abramowitz and Stegun, p. 374, 9.6.1.
  23. Greiner, Walter; Reinhardt, Joachim. Quantum Electrodynamics. Springer, 2009, p. 72. ISBN 978-3-540-87561-1. 
  24. Watson, p. 181.
  25. Khokonov, M. Kh. «Cascade Processes of Energy Loss by Emission of Hard Photons». Journal of Experimental and Theoretical Physics, 99, 4, 2004, pàg. 690–707. DOI: 10.1134/1.1826160.. Derived from formulas sourced to I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Fizmatgiz, Moscow, 1963; Academic Press, New York, 1980).
  26. Referred to as such in: Teichroew, D. «The Mixture of Normal Distributions with Different Variances». The Annals of Mathematical Statistics, 28, 2, 1957, pàg. 510–512. DOI: 10.1214/aoms/1177706981.
  27. Abramowitz and Stegun, p. 437, 10.1.1.
  28. Abramowitz and Stegun, p. 439, 10.1.25, 10.1.26.
  29. Abramowitz and Stegun, p. 438, 10.1.11.
  30. Abramowitz and Stegun, p. 438, 10.1.12.
  31. Abramowitz and Stegun, p. 439, 10.1.39.
  32. Abramowitz and Stegun, p. 439, 10.1.23, 10.1.24.
  33. Hong Du, "Mie-scattering calculation," Applied Optics 43 (9), 1951–1956 (2004)