Fonction de Bessel
En mathématiques, et plus précisément en analyse, les fonctions de Bessel, appelées aussi quelquefois fonctions cylindriques<ref>Albert Wangerin, « Fonctions cylindriques ou fonctions de Bessel », dans Encyclopédie des sciences mathématiques pures et appliquées. Tome II. Cinquième volume, J. Molk (éd.), Paris, Gauthier-Villars, 1912, p. 209.</ref>, découvertes par le mathématicien suisse Daniel Bernoulli, portent le nom du mathématicien allemand Friedrich Wilhelm Bessel. Bessel développa l'analyse de ces fonctions en 1816 dans le cadre de ses études du mouvement des planètes induit par l'interaction gravitationnelle, généralisant les découvertes antérieures de Bernoulli. Ces fonctions sont des solutions canoniques y(x) de l'équation différentielle de Bessel :
- <math>x^2 \frac{\mathrm d^2 y}{\mathrm dx^2} + x \frac{\mathrm dy}{\mathrm dx} + (x^2 - \alpha^2)y = 0</math>
pour tout nombre réel ou complexe α. Le plus souvent, α est un entier naturel (alors appelé l'ordre de la fonction), ou un demi-entier.
Il existe deux sortes de fonctions de Bessel :
- les fonctions de Bessel de première espèce, Jn, solutions de l'équation différentielle ci-dessus qui sont définies en 0 ;
- les fonctions de Bessel de seconde espèce, Yn, solutions qui ne sont pas définies en 0 (mais qui ont une limite infinie en 0).
Les représentations graphiques des fonctions de Bessel ressemblent à celles des fonctions sinus ou cosinus, mais s'amortissent comme s'il s'agissait de fonctions sinus ou cosinus divisées par un terme de la forme Modèle:Sqrt.
Les fonctions de Bessel sont aussi connues sous le nom de fonctions cylindriques, ou d'harmoniques cylindriques, parce qu'elles font partie des solutions de l'équation de Laplace en coordonnées cylindriques (intervenant, par exemple, dans la propagation de la chaleur dans un cylindre).
Elles interviennent dans beaucoup de problèmes physiques présentant une symétrie cylindrique:
- les ondes électromagnétiques ou les ondes acoustiques dans un guide cylindrique (antenne ou tuyau) ;
- les modes de vibration d'une fine membrane circulaire ou annulaire ;
- l'étude d'instruments d'optique comme les fibres optiques constituées d'un cœur et d'une gaine optique concentriques;
- le pendule de Bessel ;
- les phénomènes de diffraction par une fente circulaire ;
- l'étude de la modulation de fréquence en télécommunications.
Expression des fonctions de Bessel
Pour les valeurs entières de Modèle:Math, les fonctions de Bessel de première espèce Modèle:Mvar sont définies par la série entière (de rayon de convergence infini) suivante<ref>Modèle:Note autre projet</ref> :
- <math>J_n(x)=\sum_{p=0}^\infty {(-1)^p \over p! (n+p)!} {\left({x \over 2}\right)}^{2p+n}</math>.
Plus généralement, pour α non entier, on a le développement analogue
- <math> J_\alpha(x) = \sum_{p=0}^\infty \frac{(-1)^p}{p! \, \Gamma(p+\alpha+1)} {\left({x \over 2}\right)}^{2p+\alpha}</math>
où Modèle:Math est la fonction gamma, généralisant la fonction factorielle à des valeurs non entières.
Les fonctions de Bessel de deuxième espèce, également appelées fonctions de Neumann ou encore fonctions de Weber-Schläfli, sont définies par :
- <math>Y_n(x)=\lim_{\lambda \to n} {J_\lambda(x) \cos(\lambda \pi) - J_{-\lambda}(x) \over \sin(\lambda \pi)}</math>.
Intégrales de Bessel
Pour les valeurs entières de Modèle:Math, les fonctions de Bessel peuvent être représentées par des intégrales :
- <math>J_n(x) = \frac1{\pi} \int_0^\pi \cos (n \tau - x \sin \tau) \,\mathrm{d}\tau</math>
ou encore par :
- <math>J_n (x) = \frac1{2 \pi} \int_{-\pi}^\pi{\rm e}^{-\mathrm i\,(n \tau - x \sin \tau)} \,\mathrm d\tau</math>.
C'est la définition qu'en donna Bessel, et qui lui servit à obtenir de nombreuses propriétés de ces fonctions (à commencer par l'équation différentielle, qui en découle par différentiation sous le signe d'intégration, suivie d'une intégration par parties). Cette définition peut s'étendre au cas α non entier (pour Re(x) > 0), en ajoutant un autre terme<ref>Modèle:Lien web.</ref>,<ref>Modèle:Lien web.</ref>,<ref name=Watson/>,<ref name="Arfken & Weber"/>:
- <math>J_\alpha(x) =
\frac1{\pi} \int_0^\pi \cos(\alpha\tau- x \sin\tau)\,{\rm d}\tau - \frac{\sin(\alpha\pi)}{\pi} \int_0^\infty {\rm e}^{-x \sinh(t) - \alpha t} \, {\rm d}t.</math>
Relation avec les séries hypergéométriques
Les fonctions de Bessel peuvent également s'exprimer sous forme de série hypergéométrique comme
- <math>J_\alpha(x)=\frac{(x/2)^\alpha}{\Gamma(\alpha+1)}\;_0F_1 (\alpha+1; -\tfrac{x^2}4)</math>.
Cette expression est liée au développement des fonctions de Bessel à l'aide de la Modèle:Lien.
Relation avec les polynômes de Laguerre
Notant Modèle:Mvar le k-ième polynôme de Laguerre, les fonctions de Bessel peuvent être exprimées ainsi<ref>{{#invoke:Langue|indicationDeLangue}} Gábor Szegő, Orthogonal Polynomials, Modèle:4e éd., Providence, RI, AMS, 1975.</ref> :
- <math>\frac{J_\alpha(x)}{\left( \frac{x}{2}\right)^\alpha}= \frac{{\rm e}^{-t}}{\alpha!} \sum_{k=0}^\infty \frac{1}Modèle:K+ \alpha \choose k L_k^{(\alpha)}\left( \frac{x^2}{4 t}\right) \frac{t^k}{k!}</math>,
où l'expression de droite ne dépend pas de t et demande, pour être généralisée au cas α non entier, l'utilisation de dérivées fractionnaires.
Propriétés des Modèle:Math
- Relations de récurrence<ref name=AbramowitzStegun>Modèle:Abramowitz et Stegun, Modèle:P..</ref> :
- <math>J_{\alpha-1}(z)+J_{\alpha+1}(z)={2\alpha\over z} J_\alpha(z)</math>,
- <math>J_{\alpha-1}(z)-J_{\alpha+1}(z)=2J'_\alpha(z)</math>,
- <math>J'_\alpha(z)=J_{\alpha-1}(z)-{\alpha\over z}J_\alpha(z)=-J_{\alpha+1}(z)+{\alpha\over z}J_\alpha(z)</math>.
- On en déduit :
- <math>J_1(x)=-J_0'(x)</math>,
- <math>\frac{\mathrm d}{\mathrm dx}(x^n J_n(x))=x^n J_{n-1}(x)</math>.
- Orthogonalité :
- <math>\lambda_i</math> et <math>\lambda_j</math> étant deux zéros distincts de <math>J_n</math>, on a : <math>\int_0^1x J_n(\lambda_i x) J_n(\lambda_j x)\,\mathrm dx = 0</math>.
Jn est souvent défini par l'intermédiaire d'une série de Laurent, correspondant à la fonction génératrice :
- <math>{\rm e}^{(x/2)(t-1/t)} = \sum_{n=-\infty}^\infty J_n(x) t^n</math> ;
cette approche est celle de Peter Andreas Hansen en 1843. Elle peut se généraliser à des ordres n non entiers, par l'intermédiaire, par exemple, d'intégrales de contour.
Des développements analogues; mais utilisant des séries trigonométriques, sont dus à Jacobi et Anger ; on a<ref>{{#invoke:Langue|indicationDeLangue}} Annie Cuyt, Vigdis Petersen, Brigitte Verdonk, Haakon Waadeland et William B. Jones, Handbook of Continued Fractions for Special Functions, Springer, 2008, p. 344.</ref>
- <math>{\rm e}^{{\rm i}z \cos \phi} = \sum_{n=-\infty}^\infty{\rm i}^n J_n(z){\rm e}^{{\rm i}n\phi},</math>
et
- <math>{\rm e}^{{\rm i}z \sin \phi} = \sum_{n=-\infty}^\infty J_n(z){\rm e}^{{\rm i}n\phi}.</math>
Développements asymptotiques
Les fonctions de Bessel ont les formes asymptotiques suivantes (pour Modèle:Math). Près de 0 (et plus précisément pour <math>0 < x \ll \sqrt{\alpha + 1}</math>), on a<ref name="Arfken & Weber">{{#invoke:Langue|indicationDeLangue}} Modèle:Lien et Hans J. Weber, Mathematical Methods for Physicists, Modèle:6e éd., Harcourt, San Diego, 2005 Modèle:ISBN.</ref> :
- <math>J_\alpha(x) \approx \frac1{\Gamma(\alpha+1)} \left( \frac x2\right)^\alpha </math>
- <math>Y_\alpha(x) \approx \begin{cases}
\displaystyle \frac2{\pi} \left[ \ln \left( \frac x2 \right) + \gamma \right]& \text{si } \alpha=0 \\ \displaystyle -\frac{\Gamma(\alpha)}{\pi} \left( \frac2x\right) ^\alpha & \text{si } \alpha > 0 \end{cases} </math> où Modèle:Math est la constante d'Euler-Mascheroni (0,577…) et Modèle:Math est la fonction gamma. Pour x tendant vers l'infini (et plus précisément pour <math>x \gg |\alpha^2 - 1/4|</math>), ces développements deviennent<ref name="Arfken & Weber" /> :
- <math>J_\alpha(x)\approx \sqrt{\frac2{\pi x}}\cos \left( x-\frac{\alpha\pi}2- \frac{\pi}4\right)</math>
- <math>Y_\alpha(x) \approx \sqrt{\frac2{\pi x}}\sin \left( x-\frac{\alpha\pi}2- \frac{\pi}4\right)</math>.
La forme asymptotique ci-dessus pour Modèle:Mvar est aussi un équivalent pour Modèle:Mvar (complexe non nul) fixé, quand Modèle:Mvar tend vers Modèle:Math. Autrement dit<ref>{{#invoke:Langue|indicationDeLangue}} NIST Digital Library of Mathematical Functions, §10.19(i).</ref> :
- <math>J_\alpha(x)\underset{\overset{\alpha\to+\infty}{}}\sim\frac1\sqrt{2\pi\alpha}\left(\frac{\mathrm e\,x}{2\alpha}\right)^\alpha
</math>.
Zéros
Comme Modèle:Mvar est une solution non nulle d'une équation différentielle linéaire d'ordre 2, ses zéros (à l'exception éventuelle de Modèle:Math) sont simples<ref>Modèle:Ouvrage.</ref>, donc (cf. Relations de récurrence ci-dessus) différents de ceux de Modèle:Math.
La conjecture de Bourget
Bessel avait démontré que pour n entier positif, Jn(x) admet une infinité de zéros<ref>{{#invoke:Langue|indicationDeLangue}} F. Bessel, « Untersuchung des Theils der planetarischen Störungen », Berlin Abhandlungen, 1824, article 14.</ref>. Cependant, les graphes de Jn semblent montrer que ces zéros sont distincts pour différentes valeurs de n, en dehors de Jn(0) = 0. Ce phénomène est appelé la conjecture de Bourget<ref>Modèle:Article.</ref> ; elle fut démontrée par Carl Siegel en 1929<ref name=Watson>{{#invoke:Langue|indicationDeLangue}} G. N. Watson, A Treatise on the Theory of Bessel Functions, Modèle:2e éd., 1995, Cambridge University Press, Modèle:Lire en ligne, p. 176 et 484-485.</ref>.
Transcendance
Siegel a également démontré en 1929 que lorsque Modèle:Mvar est rationnel et Modèle:Mvar est un nombre algébrique non nul, Modèle:Math, les nombres Modèle:Math et Modèle:Math sont transcendants<ref>Modèle:Ouvrage.</ref>, de même que la valeur en Modèle:Mvar de la fonction de Bessel modifiée Modèle:Mvar<ref>Modèle:Article.</ref>. On sait aussi que toutes les racines des dérivées d'ordre supérieur <math>J_\alpha^{(n)}</math> pour Modèle:Math sont transcendantes, sauf les cas particuliers <math>J_1^{(3)}(\pm\sqrt3)=0</math> et <math>J_0^{(4)}(\pm\sqrt3)=0</math><ref name=lorch>Modèle:Article.</ref>.
Notes et références
Modèle:Traduction/Référence Modèle:Références
