În matematică, prin funcții Bessel se înțeleg soluțiile canonice Z(z) ale ecuației diferențiale a lui Bessel (cu z real sau complex):

pentru o valoare arbitrară α reală sau complexă, numită ordinul funcției Bessel. Cele mai comune și mai importante cazuri fiind acelea în care α are o valoare întreagă n.

De altfel, α și −α produc aceeași ecuație diferențială, convențional definindu-se funcții Bessel diferite pentru cele două ordine, dar cel mai adesea sunt alese ca funcții netede de α. De asemenea funcțiile Bessel sunt cunoscute ca funcții cilindrice sau cilindrice armonice deoarece ele se regăsesc în soluția ecuației Laplace în coordonate cilindrice.

Ele au fost definite prima dată de Daniel Bernoulli și generalizate de Friderich Bessel, de unde și denumirea lor.

Aplicații ale funcțiilor Bessel Prin aplicarea metodei separării variabilelor pentru soluționarea ecuației Laplace și a ecuației Helmholz, în coordonate cilindrice sau coordonate sferice, se obține ecuația lui Bessel, din care se obțin funcțiile Bessel.

Rezolvând ecuația în sistemul de coordonate cilindrice, obținem funcții Bessel de ordin întreg (α = n); rezolvând ecuația în sistemul de coordonate sferice, obținem funcții Bessel de ordin fracționar (α = n +1/2).

Importanța funcțiilor Bessel rezultă din faptul că soluționează multe probleme de potențal static și de propagare a undelor, de exemplu: unde electromagnetice în ghiduri de undă cilindrice, transferul de căldură în obiecte cilindrice, modurile de vibrații ale unei membrane circulare subțiri, probleme de difuzie, soluția ecuației Schrödinger radiala pentru o particulă liberă, în coordinate sferice, soluționarea unor probleme de divergentă aeroelastică și flutter în Aeroelasticitate, semnale de proces în electronică precum sinteza FM sau filtre Bessel.

Definiții Deoarece ecuația lui Bessel este o ecuație diferențială ordinară de ordinul doi, aceasta va avea două soluții liniar independente, iar datorită diverselor formulări ale funcției Bessel, în serie sau integrală, se alege forma cea mai convenabilă pentru problema care se soluționează.

Funcțiile Bessel de speța I-a : Jα Funcțiile Bessel de speța I-a, notate Jα(z), sunt soluții ale ecuației diferențiale a lui Bessel, care au valoare finită în origine z = 0 pentru valori α întregi nenegative și valoare infinită în origine pentru valori α negative diferite de întregi. Tipul de soluție, întreagă sau nu, și normalizarea funcției Jα(z) sunt definite de proprietațile de mai jos.

Funcția Bessel de speța I-a este definită de următoarea serie Taylor în jurul originii z = 0:

unde Γ(z) este funcția Gamma a lui Euler, care reprezinta generalizarea funcției factorial pentru valori z diferite de întregi. Graficul funcției Bessel oscilează ca cel al funcției sinus sau cosinus, diferența fiind aceea că funcția Bessel descrește proporțional cu spre infinit, precum și faptul că rădăcinile nu sunt în general periodice, cu excepția celor asimptotice pentru valori mari ale lui z.

Pentru valori α diferite de întregi, funcțiile Jα(z) și J-α(z) sunt liniar independente, reprezentând cele două soluții ale ecuației diferențiale. Pe de altă parte, pentru α de ordin întreg, este valabilă următoarea relație (de notat că funcția Gamma devine infinită pentru argumente întregi negative):

acest lucru arătând că cele două soluții nu sunt liniar independente. În acest caz, a doua soluție liniar independentă este dată de funcția Bessel de speța a II-a.

Integrale Bessel O altă modalitate de definire a funcțiilor Bessel este cea a reprezentărilor integrale.

Pentru n întreg avem reprezentarea:

Acesta a fost modul de reprezentare folosit de Bessel, iar din această definiție au derivat mai multe proprietăți ale funcției.

Adăugând un termen suplimentar integralei de mai sus, definiția poate fi extinsă și pentru valori α diferite de întregi, reprezentarea ei fiind dată de:

O altă reprezentare integrală este și:

Relația cu seria hipergeometrică Funcția Bessel poate fi exprimată în termenii seriei hipergeometrice a lui Gauss astfel:

Această expresie se referă la dezvoltarea funcției Bessel în termenii funcției Bessel-Clifford.

Relația cu polinoamele lui Laguerre În termenii polinoamelor Laguerre, pentru orice parametru t, funcția Bessel se poate exprima astfel:

Szegö, G. Orthogonal Polynomials, 4th ed. Providence, RI: Amer. Math. Soc., 1975.

Funcțiile Bessel de speța a II-a : Yα

Funcțiile Bessel de speța a II-a, notate prin Yα(z), sunt de asemenea soluții ale ecuației diferențiale a lui Bessel. Ele au o singularitate infinită în origine (z = 0).

Funcția Yα(z) este denumită și funcția Neumann, ocazional fiind notată și cu Nα(z).

Pentru valori α diferite de întregi funcția Bessel de speța a II-a se scrise în funcție de Jα(z) sub forma:

În cazul în care α are o valoare întregă n, funcția se defineste ca limită de α → n:

scriindu-se sub formă integrală:

iar sub formă de serie:

în care ψ este funcția digamma.

În cazul în care α are o valoare diferită de întreg, funcția Yα(z) este inutilă (putând fi înlocuită oricând cu J-α(z)). Pe de altă parte, când α este un întreg n, Yn(z) este a doua soluție liniar independentă a ecuației lui Bessel. Mai mult, este valabilă o relație similară cu cea pentru funcția de speța I-a, adică:

Ambele funcții, Jα(z) and Yα(z), sunt funcții olomorfe de z în planul complex cu tăietură de-a lungul axei reale negative. Când α este un întreg, funcțiile Bessel J sunt funcții întregi de z. Dacă z este fizat, atunci funcțiile Bessel sunt funcții întregi de α.

Funcțiile Hankel : Hα Un alt mod de formulare a două soluții liniar independente ale ecuației lui Bessel sunt funcțiile Hankel , Hα(1)(z) și Hα(2)(z), definite prin:

unde i este unitatea imaginară. Aceaste combinații liniar independente sunt cunoscute și sub numele de funcții Bessel de speța a III-a. Funcțiile Hankel de prima și a doua speță sunt folosite la exprimarea soluției propagării undelor cilindrice, respectiv spre exterior sau interior, în funcție de convenția de semn aleasă pentru frecventă. Ele au fost denumite după numele lui Hermann Hankel.

Folosind relațiile de definiție de mai sus ele pot fi scrise și sub forma:

în care, dacă α este un întreg, este necesară calcularea limitei expresiilor de mai sus. Sunt valide următoarele relații indiferent dacă α este sau nu întreg:

Funcțiile Bessel modificate : Iα, Kα Funcțiile Bessel sunt definite și pentru argumente complexe ale lui z, iar un caz special important este acela al argumentului pur imaginar. În acest caz, soluțiile ecuației lui Bessel se numesc funcții Bessel modificate (sau câteodată funcții Bessel hiperbolice) de prima și a doua speță, fiind definite prin oricare din următoarele relații echivalente:

Acestea au fost alese astfel încât să aibă valori reale pentru argumente z reale și pozitive. Astfel, seria obținută pentru Iα(z) este similară cu Jα(z), dar fără a avea factorul alternant (−1)m.

Iα(z) și Kα(z) sunt cele două soluții liniar independente ale ecuației modificate a lui Bessel:

Spre deosebire de funcțiile Bessel ordinare, care oscilează ca funcții de argument real, Iα(z) și Kα(z), sunt funcții care cresc și descresc exponențial. Ca și funcția Jα(z), funcția Iα(z) are valoarea zero în punctul (z = 0) pentru α > 0 și valoarea unu pentru α = 0. În mod analog, Kα(z) este divergentă în punctul (z = 0) și tinde spre zero când z → ∞.

Funcțiile Bessel sferice : jn, yn

Când se rezolvată ecuația lui Helmhotz în coordonate sferice prin separarea variabilelor, ecuația radială capătă forma:

Cele doua soluții liniar independente ale acestei ecuații se numesc funcțiile Bessel sferice jn(z) și yn(z), iar când acestea sunt scrise cu ajutorul funcțiilor Jn(z) și Yn(z) capată formele:

Abramowitz and Stegun, p. 437, 10.1.1.

De asemenea yn(z) mai este notată cu nn(z) sau ηn(z). Unii autori numesc aceste relații funcțiile Neumann sferice.

Funcțiile Bessel sferice se mai pot scrie și sub forma:

Funcția Bessel sferică de speța I-a, j0(z), mai este cunoscuta și sub numele de funcția sinc.

În continuare se dau câteva funcții Bessel sferice:

Abramowitz and Stegun, p. 438, 10.1.11.

și

Abramowitz and Stegun, p. 438, 10.1.12;

Identitatea generală fiind:

Relații diferențiale În următoarea ecuație, fn, poate fi oricare din funcțiile jn, yn, hn(1), hn(2), unde n = ±1,±2,… .

Funcțiile Hankel sferice : h n De asemenea există funcții Hankel sferice:

De fapt, ele sunt simple expresii ale funcțiilor Bessel de ordinul (n+1/2) în termenii funcțiilor trigonometrice. În particular, pentru valori n întregi nenegative avem expresia:

iar hn(2) este funcția complex conjugată a acesteia pentru z real.

Funcțiile Riccati-Bessel : Funcțiile Riccati-Bessel diferă puțin de funcțiile Bessel sferice, fiind date de formulele:

Ele satisfac ecuația diferențială:

Acestă ecuație diferențială și soluția ei Riccati-Bessel apar în problema împrăștierii undelor elecromagnetice printr-o sferă, cunoscută ca împrăștierea Mie.

Câteodata se folosesc și notațiile ψn, χn în loc de Sn, Cn.

Formele asimptotice Funcțiile Bessel au următoarele forme asimptotice pentru valori α nenegative. Pentru valori mici ale argumentelor , obtinem:Arfken & Weber.

unde γ este constanta Euler–Mascheroni (0.5772...).

Pentru argumente mari , acestea devin:

De reținut că, pentru α=1/2 aceste formule sunt exacte (vezi funcțiile Bessel sferice de mai sus). În continuare se dau formele asimptotice și pentru alte tipuri de funcții Bessel. De exemplu, pentru valori , funcțiile Bessel modificate devin:

în timp ce pentru argumente mici , ele devin:

Pentru funcțiile Bessel de prima, a doua și a treia speță, avem următoarele dezvoltări asimptotice: pentru : pentru : pentru : pentru : în care: și , iar: Abramowitz and Stegun, p. 364, 9.2.9. și Abramowitz and Stegun, p. 364, 9.2.10.

Dacă ν este real nenegativ, iar z pozitiv, restul obținut după sumarea a n termeni din expresia P(ν,z), nu depășește termenul (n+1) în valoare absolută având același semn cu el, facându-se specificația că .

Același lucru este adevarat și pentru Q(ν,z), dar cu .

Proprietăți Pentru α = n (întreg), Jn este folosită adesea prin intermediul seriei Laurent pentru funcția generatoare:

care poate fi generalizată pentru orice α prin metoda integrării pe contur sau altă metodă. O altă relație importantă pentru ordine întregi este expansiunea Jacobi–Anger:

folosită pentru a dezvolta o undă plană ca o sumă de unde cilindrice, sau pentru a găsi seria Fourier a unui semnal modulat în frecvență.

Mai general, o serie:

este numită dezvoltarea Neumann de `f`. Pentru ν = 0 are forma explicită: unde este polinomul lui Neumann Abramowitz and Stegun, p. 363, 9.1.82 ff..

Unele funcții admit reprezentarea specială:

în care :

datorită relației ortogonale:

Mai general, dacă f are un punct de ramificație în jurul originii de așa natură încât:

atunci:

sau

unde este transformata Laplace a lui f.E. T. Whittaker, G. N. Watson, A course in modern Analysis p. 536

Un alt mod de definire a funcțiilor Bessel este dat de formula lui Poisson:

unde k > -1/2, iar z este un număr complex.I.S. Gradshteyn (И.С. Градштейн), I.M. Ryzhik (И.М. Рыжик); Alan Jeffrey, Daniel Zwillinger, editors. Table of Integrals, Series, and Products, seventh edition. Academic Press, 2007. ISBN 978-0-12-373637-6. Equation 8.411.10 Acestă formulă este folositoare în special când lucrăm cu transformata Fourier.

Funcțiile Jn, Jn, Hα(1) și Hα(2) satisfac următoarele relații de recurență:

unde Z poate fi oricare din funcțiile J, Y, H(1), sau H(2). Aceste două relații sunt adesea combinate prin adunare sau scădere, pentru obținerea altor relații. De exemplu, putem calcula funcții Bessel de ordin înalt (sau derivate de ordin înalt) dând valori de ordin scăzut (sau derivate de ordin scăzut).

În particular avem:

Funcțiile Bessel modificate se scriu prin relații similare:

și

Relațiile de recurență se scriu: