Funzioni di Airy

In matematica le funzioni di Airy sono due funzioni speciali indicate rispettivamente con A i ( x ) {\displaystyle \mathrm {Ai} (x)} e B i ( x ) {\displaystyle \mathrm {Bi} (x)} che traggono il nome da quello dell'astronomo inglese George Biddell Airy (1801-1892). Esse costituiscono le soluzioni dell'equazione differenziale ordinaria, detta "di Airy",

f x f = 0. {\displaystyle f''-xf=0.}

Questa è la più semplice equazione differenziale lineare del secondo ordine dotata di un punto in cui il carattere delle soluzioni passa da oscillatorio a esponenziale. Spesso con il nome di "funzione di Airy" si intende la sola A i ( x ) {\displaystyle \mathrm {Ai} (x)} . Tale funzione può sorgere per esempio dall'equazione di Helmholtz in una sola dimensione (ordinaria):

f + k ( x ) 2 f = 0 , {\displaystyle f''+k(x)^{2}f=0,}

nel caso in cui la componente del vettore d'onda dipenda dalla radice della direzione:

k ( x ) 2 = x . {\displaystyle k(x)^{2}=x.}

Storia

La funzione di Airy prende il nome dall'astronomo inglese George Biddell Airy, che la incontrò nei suoi studi di ottica (Airy 1838). La notazione A i ( x ) {\displaystyle \mathrm {Ai} (x)} fu introdotta da Harold Jeffreys. Airy è diventato l'astronomo reale inglese nel 1835, e tenne il posto fino al suo pensionamento nel 1881.

Definizioni

Grafico di A i ( x ) {\displaystyle \mathrm {Ai} (x)} in rosso e di B i ( x ) {\displaystyle \mathrm {Bi} (x)} in blu

Per valori reali della x {\displaystyle x} , la funzione di Airy A i ( x ) {\displaystyle \mathrm {Ai} (x)} viene definita dal seguente integrale improprio:

A i ( x ) := 1 π 0 cos ( t 3 3 + x t ) d t 1 π lim b 0 b cos ( t 3 3 + x t ) d t . {\displaystyle \mathrm {Ai} (x):={\frac {1}{\pi }}\int _{0}^{\infty }\cos \left({\frac {t^{3}}{3}}+xt\right)dt\equiv {\dfrac {1}{\pi }}\lim _{b\to \infty }\int _{0}^{b}\cos \left({\dfrac {t^{3}}{3}}+xt\right)\,dt.}

L'integrale quando b {\displaystyle b\to \infty } converge anche se l'integrando non si annulla a causa delle rapide oscillazioni, per il lemma di Riemann-Lebesgue (la loro presenza può essere verificata effettuando una integrazione per parti).

Derivando sotto il simbolo di integrale, si ottiene che f = A i ( x ) {\displaystyle f=\mathrm {Ai} (x)} soddisfa l'equazione differenziale di Airy:

f x f = 0. {\displaystyle f''-xf=0.}

Questa equazione ha due soluzioni linearmente indipendenti. A meno di una costante moltiplicativa, A i ( x ) {\displaystyle \mathrm {Ai} (x)} è la soluzione soggetta alla condizione f 0 {\displaystyle f\to 0} se x + {\displaystyle x\to +\infty } . La scelta standard per l'altra soluzione è la funzione di Airy del secondo tipo, indicata con B i ( x ) {\displaystyle \mathrm {Bi} (x)} . Questa soluzione ha la stessa ampiezza di oscillazione di A i ( x ) {\displaystyle \mathrm {Ai} (x)} per x {\displaystyle x\to -\infty } , ma sfasata di π / 2 {\displaystyle \pi /2} .

B i ( x ) = 1 π 0 [ exp ( t 3 3 + x t ) + sin ( t 3 3 + x t ) ] d t . {\displaystyle \mathrm {Bi} (x)={\frac {1}{\pi }}\int _{0}^{\infty }\left[\exp \left(-{\tfrac {t^{3}}{3}}+xt\right)+\sin \left({\tfrac {t^{3}}{3}}+xt\right)\,\right]dt.}

Proprietà

I valori di A i ( x ) {\displaystyle \mathrm {Ai} (x)} e B i ( x ) {\displaystyle \mathrm {Bi} (x)} e delle loro derivate per x = 0 {\displaystyle x=0} sono dati da

A i ( 0 ) = 1 3 2 3 Γ ( 2 3 ) , A i ( 0 ) = 1 3 1 3 Γ ( 1 3 ) , B i ( 0 ) = 1 3 1 6 Γ ( 2 3 ) , B i ( 0 ) = 3 1 6 Γ ( 1 3 ) . {\displaystyle {\begin{aligned}\mathrm {Ai} (0)&{}={\frac {1}{3^{\frac {2}{3}}\Gamma ({\tfrac {2}{3}})}},&\quad \mathrm {Ai} '(0)&{}=-{\frac {1}{3^{\frac {1}{3}}\Gamma ({\tfrac {1}{3}})}},\\\mathrm {Bi} (0)&{}={\frac {1}{3^{\frac {1}{6}}\Gamma ({\tfrac {2}{3}})}},&\quad \mathrm {Bi} '(0)&{}={\frac {3^{\frac {1}{6}}}{\Gamma ({\tfrac {1}{3}})}}.\end{aligned}}}

Qui Γ {\displaystyle \Gamma } denota la funzione Gamma. Segue che il Wronskiano di A i ( x ) {\displaystyle \mathrm {Ai} (x)} e B i ( x ) {\displaystyle \mathrm {Bi} (x)} per x = 0 {\displaystyle x=0} vale 1 / π {\displaystyle 1/\pi } .

Quando x {\displaystyle x} è positivo, A i ( x ) {\displaystyle \mathrm {Ai} (x)} è positiva, concava e decrescente esponenzialmente a zero, mentre B i ( x ) {\displaystyle \mathrm {Bi} (x)} è positiva, convessa e crescente esponenzialmente. Quando x {\displaystyle x} è negativo, A i ( x ) {\displaystyle \mathrm {Ai} (x)} e B i ( x ) {\displaystyle \mathrm {Bi} (x)} oscillano intorno a zero con frequenza crescente e ampiezza decrescente. Questo è ottenibile dalle sottostanti formule asintotiche delle funzioni di Airy.

Le funzioni di Airy sono ortogonali,[1] nel senso che

A i ( t + x ) A i ( t + y ) d t = δ ( x y ) . {\displaystyle \int _{-\infty }^{\infty }\mathrm {Ai} (t+x)\mathrm {Ai} (t+y)dt=\delta (x-y).}

Argomenti complessi

Si può estendere la definizione di funzione di Airy al piano complesso definendo

A i ( z ) = 1 2 π i exp ( t 3 3 z t ) d t , {\displaystyle \mathrm {Ai} (z)={\frac {1}{2\pi i}}\int \exp \left({\frac {t^{3}}{3}}-zt\right)\,dt,}

dove l'integrale è definito su un percorso che inizia in un punto all'infinito con argomento π / 3 {\displaystyle -\pi /3} e finisce in un punto all'infinito con argomento π / 3 {\displaystyle \pi /3} . Alternativamente, possiamo usare l'equazione differenziale f x f = 0 {\displaystyle f''-xf=0} per estendere A i ( x ) {\displaystyle \mathrm {Ai} (x)} e B i ( x ) {\displaystyle \mathrm {Bi} (x)} a funzioni intere sul piano complesso.

Grafici

[ A i ( x + i y ) ] {\displaystyle \Re \left[\mathrm {Ai} (x+iy)\right]} [ A i ( x + i y ) ] {\displaystyle \Im \left[\mathrm {Ai} (x+iy)\right]} | A i ( x + i y ) | {\displaystyle |\mathrm {Ai} (x+iy)|\,} arg [ A i ( x + i y ) ] {\displaystyle \arg \left[\mathrm {Ai} (x+iy)\right]\,}
[ B i ( x + i y ) ] {\displaystyle \Re \left[\mathrm {Bi} (x+iy)\right]} [ B i ( x + i y ) ] {\displaystyle \Im \left[\mathrm {Bi} (x+iy)\right]} | B i ( x + i y ) | {\displaystyle |\mathrm {Bi} (x+iy)|\,} arg [ B i ( x + i y ) ] {\displaystyle \arg \left[\mathrm {Bi} (x+iy)\right]\,}

Formule asintotiche

A i {\displaystyle Ai} (in blu) e la forma asintotica sinusoidale/esponenziale di A i {\displaystyle Ai} (in violetto)
B i {\displaystyle Bi} (in blu) e la forma asintotica sinusoidale/esponenziale di B i {\displaystyle Bi} (in violetto)

Il comportamento asintotico delle funzioni di Airy con | z | {\displaystyle |z|} tendente all'infinito mantenendo costante il valore di arg ( z ) {\displaystyle \arg(z)} dipende da quest'ultimo: questo è chiamato il fenomeno di Stokes. Per | arg ( z ) | < π {\displaystyle |\arg(z)|<\pi } si ha la seguente stima asintotica per A i ( z ) {\displaystyle \mathrm {Ai} (z)} :[2]

A i ( z ) e 2 3 z 3 2 2 π z 1 4 [ n = 0 ( 1 ) n Γ ( n + 5 6 ) Γ ( n + 1 6 ) ( 3 4 ) n 2 π n ! z 3 n / 2 ] . {\displaystyle \mathrm {Ai} (z)\sim {\dfrac {e^{-{\frac {2}{3}}z^{\frac {3}{2}}}}{2{\sqrt {\pi }}\,z^{\frac {1}{4}}}}\left[\sum _{n=0}^{\infty }{\dfrac {(-1)^{n}\Gamma (n+{\frac {5}{6}})\Gamma (n+{\frac {1}{6}})\left({\frac {3}{4}}\right)^{n}}{2\pi n!z^{3n/2}}}\right].}

Se ne ha una uguale per B i ( z ) {\displaystyle \mathrm {Bi} (z)} , ma applicabile solo quando | arg ( z ) | < π / 3 {\displaystyle |\arg(z)|<\pi /3} :

B i ( z ) e 2 3 z 3 2 π z 1 4 [ n = 0 Γ ( n + 5 6 ) Γ ( n + 1 6 ) ( 3 4 ) n 2 π n ! z 3 n / 2 ] . {\displaystyle \mathrm {Bi} (z)\sim {\frac {e^{{\frac {2}{3}}z^{\frac {3}{2}}}}{{\sqrt {\pi }}\,z^{\frac {1}{4}}}}\left[\sum _{n=0}^{\infty }{\dfrac {\Gamma (n+{\frac {5}{6}})\Gamma (n+{\frac {1}{6}})\left({\frac {3}{4}}\right)^{n}}{2\pi n!z^{3n/2}}}\right].}

Delle formule più accurate per A i ( z ) {\displaystyle \mathrm {Ai} (z)} e per B i ( z ) {\displaystyle \mathrm {Bi} (z)} quando π / 3 < | arg ( z ) | < π {\displaystyle \pi /3<|\arg(z)|<\pi } o, equivalentemente, per A i ( z ) {\displaystyle \mathrm {Ai} (-z)} e B i ( z ) {\displaystyle \mathrm {Bi} (-z)} quando | arg ( z ) | < 2 π / 3 {\displaystyle |\arg(z)|<2\pi /3} ma non zero, sono:[3]

A i ( z ) sin ( 2 3 z 3 2 + π 4 ) π z 1 4 B i ( z ) cos ( 2 3 z 3 2 + π 4 ) π z 1 4 . {\displaystyle {\begin{aligned}\mathrm {Ai} (-z)&{}\sim {\frac {\sin \left({\frac {2}{3}}z^{\frac {3}{2}}+{\frac {\pi }{4}}\right)}{{\sqrt {\pi }}\,z^{\frac {1}{4}}}}\\[6pt]\mathrm {Bi} (-z)&{}\sim {\frac {\cos \left({\frac {2}{3}}z^{\frac {3}{2}}+{\frac {\pi }{4}}\right)}{{\sqrt {\pi }}\,z^{\frac {1}{4}}}}.\end{aligned}}}

Segue dal loro comportamento asintotico che sia A i ( x ) {\displaystyle \mathrm {Ai} (x)} e B i ( x ) {\displaystyle \mathrm {Bi} (x)} hanno un'infinità di zeri nell'asse reale negativo. La funzione A i ( x ) {\displaystyle \mathrm {Ai} (x)} non ha altri zeri nel piano complesso, mentre la funzione B i ( x ) {\displaystyle \mathrm {Bi} (x)} ha anche un'infinità di zeri nel settore { z C : π / 3 < | arg ( z ) | < π / 2 } {\displaystyle \{z\in \mathbb {C} :\pi /3<|\arg(z)|<\pi /2\}} .

Quando arg ( z ) = 0 {\displaystyle \arg(z)=0} , cioè per i numeri reali, queste sono buoni approssimazioni ma non sono asintotiche poiché il rapporto fra A i ( z ) {\displaystyle \mathrm {Ai} (-z)} o B i ( z ) {\displaystyle \mathrm {Bi} (-z)} e l'approssimazione sovrastante tende ad infinito ogni volta che il seno o il coseno si annullano. Stime asintotiche per questi limite sono comunque disponibili e sono elencate in (Abramowitz and Stegun, 1954) e (Olver, 1974).

Relazioni con altre funzioni speciali

Per argomenti positivi, le funzioni di Airy sono collegate alle funzioni di Bessel modificate:

A i ( x ) = 1 π 1 3 x K 1 / 3 ( 2 3 x 3 / 2 ) , {\displaystyle \mathrm {Ai} (x)={\frac {1}{\pi }}{\sqrt {{\frac {1}{3}}x}}\,K_{1/3}\left({\frac {2}{3}}x^{3/2}\right),}
B i ( x ) = 1 3 x [ I 1 / 3 ( 2 3 x 3 / 2 ) + I 1 / 3 ( 2 3 x 3 / 2 ) ] . {\displaystyle \mathrm {Bi} (x)={\sqrt {{\frac {1}{3}}x}}\left[I_{1/3}\left({\frac {2}{3}}x^{3/2}\right)+I_{-1/3}\left({\frac {2}{3}}x^{3/2}\right)\right].}

Dove I ± 1 / 3 {\displaystyle I_{\pm 1/3}} e K 1 / 3 {\displaystyle K_{1/3}} sono soluzioni di

x 2 f + x f ( x 2 + 1 / 9 ) f = 0. {\displaystyle x^{2}f''+xf'-(x^{2}+1/9)f=0.}

La derivata prima della funzione di Airy è

A i ( x ) = x π 3 K 2 3 ( 2 3 x 3 2 ) . {\displaystyle \mathrm {Ai'} (x)=-{\frac {x}{\pi {\sqrt {3}}}}\,K_{\frac {2}{3}}\left({\tfrac {2}{3}}x^{\frac {3}{2}}\right).}

Per argomenti negativi, le funzioni di Airy sono collegate alle funzioni di Bessel:

A i ( x ) = 1 3 x [ J 1 / 3 ( 2 3 x 3 / 2 ) + J 1 / 3 ( 2 3 x 3 / 2 ) ] , {\displaystyle \mathrm {Ai} (-x)={\frac {1}{3}}{\sqrt {x}}\left[J_{1/3}\left({\frac {2}{3}}x^{3/2}\right)+J_{-1/3}\left({\frac {2}{3}}x^{3/2}\right)\right],}
B i ( x ) = x 3 [ J 1 / 3 ( 2 3 x 3 / 2 ) J 1 / 3 ( 2 3 x 3 / 2 ) ] . {\displaystyle \mathrm {Bi} (-x)={\sqrt {\frac {x}{3}}}\left[J_{-1/3}\left({\frac {2}{3}}x^{3/2}\right)-J_{1/3}\left({\frac {2}{3}}x^{3/2}\right)\right].}

Dove J ± 1 / 3 {\displaystyle J_{\pm 1/3}} sono soluzioni di

x 2 f + x f + ( x 2 1 / 9 ) f = 0. {\displaystyle x^{2}f''+xf'+(x^{2}-1/9)f=0.}

Le funzioni di Scorer, che risolvono l'equazione f x f = 1 / π {\displaystyle f''-xf=1/\pi } , possono anche essere espresse in termini delle funzioni di Airy:

G i ( x ) = B i ( x ) x A i ( t ) d t + A i ( x ) 0 x B i ( t ) d t , {\displaystyle \mathrm {Gi} (x)=\mathrm {Bi} (x)\int _{x}^{\infty }\mathrm {Ai} (t)\,dt+\mathrm {Ai} (x)\int _{0}^{x}\mathrm {Bi} (t)\,dt,}
H i ( x ) = B i ( x ) x A i ( t ) d t A i ( x ) x B i ( t ) d t . {\displaystyle \mathrm {Hi} (x)=\mathrm {Bi} (x)\int _{-\infty }^{x}\mathrm {Ai} (t)\,dt-\mathrm {Ai} (x)\int _{-\infty }^{x}\mathrm {Bi} (t)\,dt.}

Trasformata di Fourier

Usando la definizione di funzione di Airy A i ( x ) {\displaystyle \mathrm {Ai} (x)} , è semplice mostrare che la sua trasformata di Fourier è data da

F ( A i ) ( k ) := A i ( x )   e 2 π i k x d x = e i 3 ( 2 π k ) 3 . {\displaystyle {\mathcal {F}}(\mathrm {Ai} )(k):=\int _{-\infty }^{\infty }\mathrm {Ai} (x)\ e^{-2\pi ikx}\,dx=e^{{\frac {i}{3}}(2\pi k)^{3}}.}

Applicazioni

La funzione di Airy è la soluzione dell'equazione di Schrödinger per una particella confinata in una buca di potenziale triangolare e per una particella in un campo unidimensionale uniforme di forze. Per lo stesso motivo, questa funzione serve a fornire un'approssimazione uniforme vicino ad un punto di svolta nell'approssimazione WKB, dove il potenziale può essere approssimato localmente da una funzione lineare della posizione. La soluzione della buca di potenziale triangolare è direttamente rilevante per la comprensione di molti dispositivi a semiconduttore.

La funzione di Airy inoltre sottolinea la forma dell'intensità vicino a una caustica ottica direzionale, come quella dell'arcobaleno. Storicamente, fu questo problema matematico che portò Airy a sviluppare questa funzione speciale.

Funzione Zeta di Airy

La funzione Zeta di Airy, studiata da Crandall (1996), è una funzione analoga alla funzione zeta di Riemann e relativa agli zeri della funzione A i ( x ) {\displaystyle \mathrm {Ai} (x)} .

Detta a 1 {\displaystyle a_{1}} , a 2 {\displaystyle a_{2}} , ... la successione degli x {\displaystyle x} in cui A i ( x ) = 0 {\displaystyle \mathrm {Ai} (x)=0} , ordinati in base al loro valore assoluto, la funzione Zeta di Airy è definita dalla serie

ζ A i ( s ) = i = 1 1 | a i | s . {\displaystyle \zeta _{\mathrm {Ai} }(s)=\sum _{i=1}^{\infty }{\frac {1}{|a_{i}|^{s}}}.}

Questa serie converge quando la parte reale di s {\displaystyle s} è maggiore di 3 / 2 {\displaystyle 3/2} e può essere estesa per prolungamento analitico ad altri valori di s {\displaystyle s} . Come la funzione Zeta di Riemann, il cui valore ζ ( 2 ) = π 2 / 6 {\displaystyle \zeta (2)=\pi ^{2}/6} è la soluzione al problema di Basilea, la funzione Zeta può essere valutata esattamente in s = 2 {\displaystyle s=2} :

ζ A i ( 2 ) = i = 1 1 a i 2 = 3 5 / 3 Γ 4 ( 2 3 ) 4 π 2 , {\displaystyle \zeta _{\mathrm {Ai} }(2)=\sum _{i=1}^{\infty }{\frac {1}{a_{i}^{2}}}={\frac {3^{5/3}\Gamma ^{4}({\frac {2}{3}})}{4\pi ^{2}}},}

dove Γ {\displaystyle \Gamma } è la funzione Gamma, una variante continua del fattoriale. Valutazioni simili sono anche possibili per valori di s {\displaystyle s} più grandi. È stato congetturato che il prolungamento analitico della funzione Zeta di Airy valutato in s = 1 {\displaystyle s=1} valga

ζ A i ( 1 ) = 3 2 / 3 Γ ( 2 3 ) Γ ( 4 3 ) . {\displaystyle \zeta _{\mathrm {Ai} }(1)={\frac {-3^{-2/3}\Gamma ({\frac {2}{3}})}{\Gamma ({\frac {4}{3}})}}.}

Note

  1. ^ David E. Aspnes, Physical Review, 147, 554 (1966)
  2. ^ Abramowitz & Stegun (1970, pag. 448), Equaz 10.4.59 e 10.4.63
  3. ^ Abramowitz & Stegun (1970, pag. 448), Equaz 10.4.60 e 10.4.64

Bibliografia

  • (EN) Milton Abramowitz e Irene A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 1954. (Vedi §10.4). National Bureau of Standards.
  • (EN) George B. Airy (1838): On the intensity of light in the neighbourhood of a caustic, Transactions of the Cambridge Philosophical Society, 6, pp. 379–402
  • (EN) Frank Olver, Asymptotics and Special Functions, Academic Press, 1997, ISBN 978-15-68-81069-0. (Cap. 11)
  • (EN) Frank Olver (2008): Chapter AI: Airy and related functions, Chapter 9 della Digital Library of Mathematical Functions.
  • (EN) Richard E. Crandall, On the quantum zeta function, in Journal of Physics A: Mathematical and General, vol. 29, n. 21, 1996, pp. 6795–6816, DOI:10.1088/0305-4470/29/21/014, ISSN 0305-4470 (WC · ACNP), MR 1421901.

Voci correlate

Altri progetti

Altri progetti

  • Wikimedia Commons
  • Collabora a Wikimedia Commons Wikimedia Commons contiene immagini o altri file su funzioni di Airy

Collegamenti esterni

  • Airy, funzioni di, in Enciclopedia della Matematica, Istituto dell'Enciclopedia Italiana, 2013. Modifica su Wikidata
  • (EN) Eric W. Weisstein, Funzioni di Airy, su MathWorld, Wolfram Research. Modifica su Wikidata
  • (EN) Funzioni di Airy, su Encyclopaedia of Mathematics, Springer e European Mathematical Society. Modifica su Wikidata
  • Funzioni di Airy in MathWorld
  • funzione di Airy Ai(x) in functions.wolfram.com
  • funzione di Airy Bi(x) in functions.wolfram.com
  • derivata della funzione di Airy Ai(x) in functions.wolfram.com
  • derivata della funzione di Airy Bi(x) in functions.wolfram.com
  • Funzione Zeta di Airy in MathWorld
Controllo di autoritàThesaurus BNCF 33538 · LCCN (EN) sh85003008 · GND (DE) 4225959-9 · BNF (FR) cb12404112c (data) · J9U (ENHE) 987007292907605171
  Portale Matematica: accedi alle voci di Wikipedia che trattano di matematica