Le funzioni di Bessel

Le funzioni di Bessel, studiate da Daniel Bernoulli (1732) e sistemate da Friedrich Bessel (1824), sono soluzioni dell'equazione differenziale di Bessel:

x^2 y''(x) + x y'(x) + (x^2 - n^2) y(x) = 0

L'equazione appare quando si risolve l'equazione di Laplace o Helmholtz in coordinate cilindriche, rendendo Bessel fondamentale in fisica matematica, elettromagnetismo, acustica, meccanica quantistica e ingegneria.

Le quattro famiglie principali

• J_n(x): prima specie, regolare in x=0. Valori iniziali: J_0(0)=1, J_n(0)=0 per n>=1.
• Y_n(x): seconda specie (Neumann/Weber), singolare in x=0. Y_n(x) -> -inf per x -> 0+.
• I_n(x): modificata prima specie, cresce esponenzialmente come e^x.
• K_n(x): modificata seconda specie (Macdonald), decade esponenzialmente come e^(-x).

Relazioni: I_n(x) = i^(-n) J_n(i*x), K_n legate a I_n via simmetria. Le funzioni di Hankel H_n^(1), H_n^(2) = J_n +/- i*Y_n usate per onde cilindriche entranti/uscenti.

Serie di potenze

La definizione standard via serie di Frobenius:

J_n(x) = Sigma_{k=0}^{inf} ((-1)^k / (k! * Gamma(k+n+1))) * (x/2)^(2k+n)

I_n(x) = Sigma_{k=0}^{inf} (1 / (k! * Gamma(k+n+1))) * (x/2)^(2k+n)

Queste serie convergono per tutti x reali, ma per x grande sono numericamente instabili. Si usano allora formule asintotiche.

Sviluppo asintotico di Hankel

Per x -> infinito, i comportamenti asintotici sono:

• J_n(x) ~ sqrt(2/(pi*x)) * cos(x - n*pi/2 - pi/4)
• Y_n(x) ~ sqrt(2/(pi*x)) * sin(x - n*pi/2 - pi/4)
• I_n(x) ~ e^x / sqrt(2*pi*x)
• K_n(x) ~ sqrt(pi/(2*x)) * e^(-x)

J_n e Y_n oscillano smorzati come 1/sqrt(x), I_n e K_n sono monotoni. Le formule di Hankel includono correzioni O(1/x), O(1/x^2), ..., permettendo precisione arbitraria per x grande.

Ricorrenza a tre termini

Tutte le funzioni di Bessel soddisfano ricorrenze simili:

(2n/x) * J_n(x) = J_{n-1}(x) + J_{n+1}(x)

(2n/x) * Y_n(x) = Y_{n-1}(x) + Y_{n+1}(x)

I_{n-1}(x) - I_{n+1}(x) = (2n/x) * I_n(x)

Queste permettono calcolo in O(n) tempo partendo da valori iniziali. Per Bessel, occorre usare ricorrenza all'indietro (Miller) per stabilità' numerica con n grande.

Le sei modalita' del calcolatore

1. J_n(x): prima specie, valutazione via serie o asintotica.
2. Y_n(x): seconda specie, con limite per n intero.
3. I_n(x): modificata crescente.
4. K_n(x): modificata decrescente.
5. Zeri di J_n: i primi j_{n,k} per quadrature e soluzioni ai limiti.
6. Asintotica: confronto serie vs formula di Hankel.

Zeri di J_n

J_n(x) ha infiniti zeri positivi j_{n,1} < j_{n,2} < j_{n,3} < .... Alcuni valori:

• J_0: 2.405, 5.520, 8.654, 11.792, 14.931, 18.071
• J_1: 3.832, 7.016, 10.173, 13.324, 16.471, 19.616
• J_2: 5.136, 8.417, 11.620, 14.796, 17.960, 21.117
• J_3: 6.380, 9.761, 13.015, 16.223, 19.409, 22.583

Asintoticamente: j_{n,k} ~ (k + n/2 - 1/4) * pi + O(1/k). Gli zeri sono critici in problemi di autovalore per equazioni cilindriche e in quadratura di Gauss-Bessel.

Applicazioni in fisica

• Elettromagnetismo: modi di propagazione in guide d'onda cilindriche, radiazione di dipoli.
• Acustica: vibrazioni di tamburi circolari (membrane), risonatori cilindrici.
• Meccanica quantistica: particella in buca cilindrica, atomo di idrogeno (parte radiale involve Laguerre, ma Bessel in coordinate cilindriche).
• Termodinamica: conduzione del calore in cilindri, diffusione transiente.
• Ottica: diffrazione di Fraunhofer di apertura circolare (disco di Airy usa J_1).
• Ingegneria strutturale: vibrazioni di piastre circolari, analisi di colonne cilindriche.

Applicazioni matematiche

• Trasformate di Hankel: analogo cilindrico della trasformata di Fourier.
• Quadratura di Gauss-Bessel: nodi dati dagli zeri j_{n,k}, pesi relativi.
• Teoria delle probabilità': distribuzione chi-quadrata non centrale, distribuzioni radiali.
• Analisi numerica: soluzione di EDP con metodi spettrali cilindrici.
• Statistica: test di goodness-of-fit con funzioni caratteristiche cilindriche.
• Crittografia: analisi di sicurezza di schemi basati su distribuzioni radiali.

Storia delle funzioni di Bessel

Le funzioni prendono il nome dall'astronomo e matematico tedesco Friedrich Wilhelm Bessel, che le studiò sistematicamente fra il 1817 e il 1824 analizzando le perturbazioni gravitazionali nel moto dei pianeti, in relazione all'equazione di Keplero. Tuttavia, casi particolari erano già comparsi in precedenza: Daniel Bernoulli le incontrò studiando le oscillazioni di una catena pesante appesa, ed Eulero nell'analisi delle vibrazioni di una membrana circolare. Fu Bessel a sviluppare la teoria generale, introducendo la famiglia di funzioni indicizzate dall'ordine n e studiandone le proprietà. Da allora le funzioni di Bessel sono diventate onnipresenti nella fisica matematica: ogni volta che un problema ha simmetria cilindrica o sferica, esse emergono naturalmente come soluzioni delle equazioni differenziali associate.

Le funzioni di Bessel e la membrana circolare vibrante

L'esempio fisico più intuitivo delle funzioni di Bessel è la vibrazione di una membrana circolare, come la pelle di un tamburo. Risolvendo l'equazione delle onde in coordinate polari, la parte radiale della soluzione è data proprio dalle funzioni di Bessel J_n. I modi di vibrazione del tamburo corrispondono agli zeri di queste funzioni: la condizione che il bordo resti fermo impone che la funzione di Bessel si annulli lì, selezionando le frequenze permesse. A differenza di una corda vibrante, le cui frequenze sono multipli interi, quelle di un tamburo sono legate agli zeri di J_n e non sono armoniche, il che spiega il suono caratteristico e non melodico dei tamburi. Lo stesso principio governa la propagazione delle onde nelle guide d'onda cilindriche, la conduzione del calore in un cilindro e i modi dei laser.

Aggiornato a luglio 2026