Polynôme de Bessel

En mathématiques, les polynômes de Bessel sont une suite de polynômes orthogonaux. Il en existe plusieurs définitions, mais toutes liées. La définition la plus courante est celle donnée par la somme[1]:

y n ( x ) = ∑ k = 0 n ( n + k ) ! ( n − k ) ! k ! ( x 2 ) k {\displaystyle y_{n}(x)=\sum _{k=0}^{n}{\frac {(n+k)!}{(n-k)!\,k!}}\,\left({\frac {x}{2}}\right)^{k}}

Une autre définition, préférée dans le traitement du signal, est parfois appelée polynômes de Bessel inverses [2],[3] :

θ n ( x ) = x n y n ( 1 x ) = ∑ k = 0 n ( n + k ) ! ( n − k ) ! k ! x n − k 2 k {\displaystyle \theta _{n}(x)=x^{n}\,y_{n}\left({\frac {1}{x}}\right)=\sum _{k=0}^{n}{\frac {(n+k)!}{(n-k)!\,k!}}\,{\frac {x^{n-k}}{2^{k}}}}

Les coefficients de la deuxième définition sont les mêmes que dans la première, mais l'ordre des monômes est inversé. On a ainsi, par exemple pour l'ordre 3 :

y 3 ( x ) = 15 x 3 + 15 x 2 + 6 x + 1 e t θ 3 ( x ) = x 3 + 6 x 2 + 15 x + 15 {\displaystyle y_{3}(x)=15x^{3}+15x^{2}+6x+1\quad \mathrm {et} \quad \theta _{3}(x)=x^{3}+6x^{2}+15x+15\,}

Cette deuxième famille est utilisée dans la conception des filtres de Bessel.

Propriétés

Définition en termes de fonctions de Bessel

Un polynôme de Bessel peut aussi être défini à partir des fonctions de Bessel, dont les polynômes tirent leur nom :

y n ( x ) = x n θ n ( 1 x ) {\displaystyle y_{n}(x)=\,x^{n}\theta _{n}\left({\frac {1}{x}}\right)\,}
y n ( x ) = 2 π x e 1 / x K n + 1 2 ( 1 x ) {\displaystyle y_{n}(x)={\sqrt {\frac {2}{\pi x}}}\,\mathrm {e} ^{1/x}K_{n+{\frac {1}{2}}}\left({\frac {1}{x}}\right)}
θ n ( x ) = 2 π x n + 1 / 2 e x K n + 1 2 ( x ) {\displaystyle \theta _{n}(x)={\sqrt {\frac {2}{\pi }}}\,x^{n+1/2}\mathrm {e} ^{x}K_{n+{\frac {1}{2}}}(x)}

où Kn(x) est une fonction de Bessel modifiée de deuxième espèce[2]:7,34. Par exemple[4]:

y 3 ( x ) = 15 x 3 + 15 x 2 + 6 x + 1 = 2 π x e 1 / x K 3 + 1 2 ( 1 x ) {\displaystyle y_{3}(x)=15x^{3}+15x^{2}+6x+1={\sqrt {\frac {2}{\pi x}}}\,\mathrm {e} ^{1/x}K_{3+{\frac {1}{2}}}\left({\frac {1}{x}}\right)}

Définition en termes de fonctions hypergéométriques

Un polynôme de Bessel peut aussi être défini comme une fonction hypergéométrique confluente[5]:

y n ( x ) = 2 F 0 ( − n , n + 1 ; ; − x / 2 ) = ( 2 x ) − n U ( − n , − 2 n , 2 x ) = ( 2 x ) n + 1 U ( n + 1 , 2 n + 2 , 2 x ) . {\displaystyle y_{n}(x)=\,_{2}F_{0}(-n,n+1;;-x/2)=\left({\frac {2}{x}}\right)^{-n}U\left(-n,-2n,{\frac {2}{x}}\right)=\left({\frac {2}{x}}\right)^{n+1}U\left(n+1,2n+2,{\frac {2}{x}}\right).}

On retrouve une expression similaire pour les polynômes de Bessel généralités (cf. infra)[2]:35:

y n ( x ; α , β ) = 2 F 0 ( − n , n + a − 1 ; ; − x / b ) = ( b x ) n + a − 1 U ( n + a − 1 , 2 n + a , b x ) . {\displaystyle y_{n}(x;\alpha ,\beta )=\,_{2}F_{0}(-n,n+a-1;;-x/b)=\left({\frac {b}{x}}\right)^{n+a-1}U\left(n+a-1,2n+a,{\frac {b}{x}}\right).}

Les polynômes de Bessel inverses peuvent s'exprimer à partir des polynômes de Laguerre généralisés :

θ n ( x ) = n ! ( − 2 ) n L n − 2 n − 1 ( 2 x ) {\displaystyle \theta _{n}(x)={\frac {n!}{(-2)^{n}}}\,L_{n}^{-2n-1}(2x)}

dont on tire une expression sous forme de fonction hypergéométrique :

θ n ( x ) = ( − 2 n ) n ( − 2 ) n 1 F 1 ( − n ; − 2 n ; 2 x ) {\displaystyle \theta _{n}(x)={\frac {(-2n)_{n}}{(-2)^{n}}}\,\,_{1}F_{1}(-n;-2n;2x)}

avec (−2n)n pour le symbole de Pochhammer.

L'inversion pour les monômes s'écrit

( 2 x ) n n ! = ( − 1 ) n ∑ j = 0 n n + 1 j + 1 ( j + 1 n − j ) L j − 2 j − 1 ( 2 x ) = 2 n n ! ∑ i = 0 n i ! ( 2 i + 1 ) ( 2 n + 1 n − i ) x i L i ( − 2 i − 1 ) ( 1 x ) . {\displaystyle {\frac {(2x)^{n}}{n!}}=(-1)^{n}\sum _{j=0}^{n}{\frac {n+1}{j+1}}{j+1 \choose n-j}L_{j}^{-2j-1}(2x)={\frac {2^{n}}{n!}}\sum _{i=0}^{n}i!(2i+1){2n+1 \choose n-i}x^{i}L_{i}^{(-2i-1)}\left({\frac {1}{x}}\right).}

Fonction génératrice

Les polynômes de Bessel, avec un décalage d'indice, ont pour fonction génératrice

∑ n = 0 ∞ 2 π x n + 1 2 e x K n − 1 2 ( x ) t n n ! = 1 + x ∑ n = 1 ∞ θ n − 1 ( x ) t n n ! = exp ⁡ ( x ( 1 − 1 − 2 t ) ) . {\displaystyle \sum _{n=0}^{\infty }{\sqrt {\frac {2}{\pi }}}x^{n+{\frac {1}{2}}}\mathrm {e} ^{x}K_{n-{\frac {1}{2}}}(x){\frac {t^{n}}{n!}}=1+x\sum _{n=1}^{\infty }\theta _{n-1}(x){\frac {t^{n}}{n!}}=\exp \left(x(1-{\sqrt {1-2t}})\right).}

En dérivant selon t et en annulant x, on obtient la fonction génératrice pour les polynômes ( θ n ) n ≥ 0 {\displaystyle \left(\theta _{n}\right)_{n\geq 0}}

∑ n = 0 ∞ θ n ( x ) t n n ! = 1 1 − 2 t e x ( 1 − 1 − 2 t ) . {\displaystyle \sum _{n=0}^{\infty }\theta _{n}(x){\frac {t^{n}}{n!}}={\frac {1}{\sqrt {1-2t}}}\mathrm {e} ^{x(1-{\sqrt {1-2t}})}.}

Une fonction génératrice similaire existe pour les ( y n ) n ≥ 0 {\displaystyle \left(y_{n}\right)_{n\geq 0}} [1]:106:

∑ n = 0 ∞ y n − 1 ( x ) t n n ! = exp ⁡ ( 1 − 1 − 2 x t x ) . {\displaystyle \sum _{n=0}^{\infty }y_{n-1}(x){\frac {t^{n}}{n!}}=\exp \left({\frac {1-{\sqrt {1-2xt}}}{x}}\right).}

En posant t = z-xz2/2, on trouve l'expression de la fonction exponentielle :

e z = ∑ n = 0 ∞ y n − 1 ( x ) ( z − x z 2 / 2 ) n n ! . {\displaystyle \mathrm {e} ^{z}=\sum _{n=0}^{\infty }y_{n-1}(x){\frac {(z-xz^{2}/2)^{n}}{n!}}.}

Formules de récurrence

Une définition utilisée plus couramment pour le calcul des valeurs des polynômes de Bessel est la formule de récurrence :

y 0 ( x ) = 1 {\displaystyle y_{0}(x)=1\,}
y 1 ( x ) = x + 1 {\displaystyle y_{1}(x)=x+1\,}
y n ( x ) = ( 2 n − 1 ) x y n − 1 ( x ) + y n − 2 ( x ) {\displaystyle y_{n}(x)=(2n\!-\!1)x\,y_{n-1}(x)+y_{n-2}(x)\,}

et

θ 0 ( x ) = 1 {\displaystyle \theta _{0}(x)=1\,}
θ 1 ( x ) = x + 1 {\displaystyle \theta _{1}(x)=x+1\,}
θ n ( x ) = ( 2 n − 1 ) θ n − 1 ( x ) + x 2 θ n − 2 ( x ) {\displaystyle \theta _{n}(x)=(2n\!-\!1)\theta _{n-1}(x)+x^{2}\theta _{n-2}(x)\,}

Équation différentielle

Les polynômes de Bessel sont les solutions polynomiales de l'équation différentielle :

x 2 d 2 y n d x 2 ( x ) + 2 ( x + 1 ) d y n ( x ) d x − n ( n + 1 ) y n ( x ) = 0 {\displaystyle x^{2}{\frac {\mathrm {d} ^{2}y_{n}}{\mathrm {d} x^{2}}}(x)+2(x\!+\!1){\frac {\mathrm {d} y_{n}(x)}{\mathrm {d} x}}-n(n+1)y_{n}(x)=0}

et

x d 2 θ n d x 2 ( x ) − 2 ( x + n ) d θ n ( x ) d x + 2 n θ n ( x ) = 0 {\displaystyle x{\frac {\mathrm {d} ^{2}\theta _{n}}{\mathrm {d} x^{2}}}(x)-2(x\!+\!n){\frac {\mathrm {d} \theta _{n}(x)}{\mathrm {d} x}}+2n\,\theta _{n}(x)=0}

Orthogonalité

Les polynômes de Bessel sont orthogonaux pour le poids e-2/x sur le cercle unité du plan complexe[1] : avec le symbole de Kronecker :

∀ ( n , m ) ∈ N 2 ,   n ≠ m ⟹ ∫ 0 2 π y n ( e i θ ) y m ( e i θ ) i e i θ d θ = 0 {\displaystyle \forall (n,m)\in \mathbb {N} ^{2},\ n\neq m\Longrightarrow \int _{0}^{2\pi }y_{n}\left(\mathrm {e} ^{\mathrm {i} \theta }\right)y_{m}\left(\mathrm {e} ^{\mathrm {i} \theta }\right)\mathrm {i} \mathrm {e} ^{\mathrm {i} \theta }\mathrm {d} \theta =0}

Généralisation

Forme explicite

Une généralisation des polynômes de Bessel ont été suggérés dans la littérature scientifique (Krall, Fink), comme suit :

y n ( x ; α , β ) := ( − 1 ) n n ! ( x β ) n L n ( − 1 − 2 n − α ) ( β x ) , {\displaystyle y_{n}(x;\alpha ,\beta ):=(-1)^{n}n!\left({\frac {x}{\beta }}\right)^{n}L_{n}^{(-1-2n-\alpha )}\left({\frac {\beta }{x}}\right),}

ce qui donne une généralisation des polynômes inverses avec :

θ n ( x ; α , β ) := n ! ( − β ) n L n ( − 1 − 2 n − α ) ( β x ) = x n y n ( 1 x ; α , β ) . {\displaystyle \theta _{n}(x;\alpha ,\beta ):={\frac {n!}{(-\beta )^{n}}}L_{n}^{(-1-2n-\alpha )}(\beta x)=x^{n}y_{n}\left({\frac {1}{x}};\alpha ,\beta \right).}

Les coefficients explicites des polynômes yn(x ; α , β) sont[1]:108:

y n ( x ; α , β ) = ∑ k = 0 n ( n k ) ( n + k + α − 2 ) k _ ( x β ) k . {\displaystyle y_{n}(x;\alpha ,\beta )=\sum _{k=0}^{n}{\binom {n}{k}}(n+k+\alpha -2)^{\underline {k}}\left({\frac {x}{\beta }}\right)^{k}.}

Par conséquent, les polynômes θn(x ; α , β) peuvent être explicitement écrits ainsi :

θ n ( x ; α , β ) = ∑ k = 0 n ( n k ) ( 2 n − k + α − 2 ) n − k _ x k β n − k . {\displaystyle \theta _{n}(x;\alpha ,\beta )=\sum _{k=0}^{n}{\binom {n}{k}}(2n-k+\alpha -2)^{\underline {n-k}}{\frac {x^{k}}{\beta ^{n-k}}}.}

Avec la fonction poids

ρ ( x ; α , β ) := 1 F 1 ( 1 , α − 1 , − β x ) {\displaystyle \rho (x;\alpha ,\beta ):=\,_{1}F_{1}\left(1,\alpha -1,-{\frac {\beta }{x}}\right)}

on retrouve l'orthogonalité :

∀ n ≠ m , ∮ c ρ ( z ; α , β ) y n ( z ; α , β ) y m ( z ; α , β ) d z = 0 {\displaystyle \forall n\neq m,\oint _{c}\rho (z;\alpha ,\beta )y_{n}(z;\alpha ,\beta )y_{m}(z;\alpha ,\beta )\mathrm {d} z=0}

avec C une courbe passant autour de l'origine.

On retrouve les polynômes de Bessel pour α = β = 2, et on obtient bien la fonction poids vue au-dessus ρ(x) = exp(−2 / x).

Formule de Rodrigues pour les polynômes de Bessel

La formule de Rodrigues pour les polynômes de Bessel comme solutions particulières des équations différentielles donnent :

B n ( α , β ) ( x ) = a n ( α , β ) x α e − β x d n d x n ( x α + 2 n e − β x ) {\displaystyle B_{n}^{(\alpha ,\beta )}(x)={\frac {a_{n}^{(\alpha ,\beta )}}{x^{\alpha }\mathrm {e} ^{-{\frac {\beta }{x}}}}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}\left(x^{\alpha +2n}\mathrm {e} ^{-{\frac {\beta }{x}}}\right)}

où a(α, β)
n sont des coefficients de normalisation.

Polynômes de Bessel généralisés

Selon cette généralisation, on trouve l'équation différentielle généralisée pour les polynômes de Bessel généralisés :

x 2 d 2 B n , m ( α , β ) ( x ) d x 2 + [ ( α + 2 ) x + β ] d B n , m ( α , β ) ( x ) d x − [ n ( α + n + 1 ) + m β x ] B n , m ( α , β ) ( x ) = 0 {\displaystyle x^{2}{\frac {\mathrm {d} ^{2}B_{n,m}^{(\alpha ,\beta )}(x)}{\mathrm {d} x^{2}}}+[(\alpha +2)x+\beta ]{\frac {\mathrm {d} B_{n,m}^{(\alpha ,\beta )}(x)}{\mathrm {d} x}}-\left[n(\alpha +n+1)+{\frac {m\beta }{x}}\right]B_{n,m}^{(\alpha ,\beta )}(x)=0}

avec 0 ≤ m ≤ n {\displaystyle 0\leq m\leq n} . Les solutions sont :

B n , m ( α , β ) ( x ) = a n , m ( α , β ) x α + m e − β x d n − m d x n − m ( x α + 2 n e − β x ) . {\displaystyle B_{n,m}^{(\alpha ,\beta )}(x)={\frac {a_{n,m}^{(\alpha ,\beta )}}{x^{\alpha +m}\mathrm {e} ^{-{\frac {\beta }{x}}}}}{\frac {\mathrm {d} ^{n-m}}{\mathrm {d} x^{n-m}}}\left(x^{\alpha +2n}\mathrm {e} ^{-{\frac {\beta }{x}}}\right).}

Zéros

En désignant les zéros de yn(x ; α , β) par α(n)
k(α , β), et ceux de θn(x ; α , β) par β(n)
k(α , β), on a les estimations suivantes[2]:82:

2 n ( n + α − 1 ) ≤ α k ( n ) ( α , 2 ) ≤ 2 n + α − 1 , {\displaystyle {\frac {2}{n(n+\alpha -1)}}\leq \alpha _{k}^{(n)}(\alpha ,2)\leq {\frac {2}{n+\alpha -1}},}

et

n + α − 1 2 ≤ β k ( n ) ( α , 2 ) ≤ n ( n + α − 1 ) 2 , {\displaystyle {\frac {n+\alpha -1}{2}}\leq \beta _{k}^{(n)}(\alpha ,2)\leq {\frac {n(n+\alpha -1)}{2}},}

pour tout α ≥ 2 {\displaystyle \alpha \geq 2} . De plus, tous ces zéros ont des parties réelles négatives.

Des résultats plus fins peuvent être donnés en utilisant des théorèmes plus puissants sur l'estimation des zéros de polynômes (comme le théorème parabolique de Saff et Varga, ou des techniques d'équations différentielles[2]:88,[6]). On a par exemple[7] :

2 2 n + α − 2 3 ≤ α k ( n ) ( α , 2 ) ≤ 2 n + α − 1 . {\displaystyle {\frac {2}{2n+\alpha -{\frac {2}{3}}}}\leq \alpha _{k}^{(n)}(\alpha ,2)\leq {\frac {2}{n+\alpha -1}}.}

Valeurs particulières

Représentation graphique des six premiers polynômes de Bessel

Les six premiers polynômes de Bessel sont :

y 0 ( x ) = 1 y 1 ( x ) = x + 1 y 2 ( x ) = 3 x 2 + 3 x + 1 y 3 ( x ) = 15 x 3 + 15 x 2 + 6 x + 1 y 4 ( x ) = 105 x 4 + 105 x 3 + 45 x 2 + 10 x + 1 y 5 ( x ) = 945 x 5 + 945 x 4 + 420 x 3 + 105 x 2 + 15 x + 1 {\displaystyle {\begin{aligned}y_{0}(x)&=1\\y_{1}(x)&=x+1\\y_{2}(x)&=3x^{2}+3x+1\\y_{3}(x)&=15x^{3}+15x^{2}+6x+1\\y_{4}(x)&=105x^{4}+105x^{3}+45x^{2}+10x+1\\y_{5}(x)&=945x^{5}+945x^{4}+420x^{3}+105x^{2}+15x+1\end{aligned}}}

Aucun polynôme de Bessel ne peut être factorisé avec des coefficients strictement rationnels[8].

Représentation graphique des six premiers polynômes de Bessel inverses

Les six premiers polynômes de Bessel inverses sont donc :

θ 0 ( x ) = 1 θ 1 ( x ) = x + 1 θ 2 ( x ) = x 2 + 3 x + 3 θ 3 ( x ) = x 3 + 6 x 2 + 15 x + 15 θ 4 ( x ) = x 4 + 10 x 3 + 45 x 2 + 105 x + 105 θ 5 ( x ) = x 5 + 15 x 4 + 105 x 3 + 420 x 2 + 945 x + 945 {\displaystyle {\begin{aligned}\theta _{0}(x)&=1\\\theta _{1}(x)&=x+1\\\theta _{2}(x)&=x^{2}+3x+3\\\theta _{3}(x)&=x^{3}+6x^{2}+15x+15\\\theta _{4}(x)&=x^{4}+10x^{3}+45x^{2}+105x+105\\\theta _{5}(x)&=x^{5}+15x^{4}+105x^{3}+420x^{2}+945x+945\end{aligned}}}

Applications

L'équation différentielle des polynômes de Bessel apparait dans l'étude de l'équation des ondes en coordonnées sphériques[9].

Articles connexes

Références

  • (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Bessel polynomials » (voir la liste des auteurs).
  1. ↑ a b c et d H. L. Krall et O. Frink, « A New Class of Orthogonal Polynomials: The Bessel Polynomials », Trans. Amer. Math. Soc., vol. 65, no 1,‎ , p. 100-115 (DOI 10.2307/1990516).
  2. ↑ a b c d et e (en) Emil Grosswald, Bessel Polynomials (Lecture Notes in Mathematics), New York, Springer, (ISBN 978-0-387-09104-4).
  3. ↑ Christian Berg et C. Vignat, « Linearization coefficients of Bessel polynomials and properties of Student-t distributions », (consulté le ).
  4. ↑ Wolfram Alpha example
  5. ↑ (en) Petre Dita et Nicolae Grama, On Adomian's Decomposition Method for Solving Differential Equations, (lire en ligne).
  6. ↑ (en) E. B. Saff et R. S. Varga, « Zero-free parabolic regions for sequences of polynomials », SIAM J. Math. Anal., vol. 7, no 3,‎ , p. 344-357 (DOI 10.1137/0507028).
  7. ↑ (en) M. G. de Bruin, E. B. Saff et R. S. Varga, « On the zeros of generalized Bessel polynomials. I », Indag. Math., vol. 84, no 1,‎ , p. 1-13.
  8. ↑ (en) Michael Filaseta et Ognian Trifinov, « The Irreducibility of the Bessel Polynomials », Journal für die Reine und Angewandte Mathematik, vol. 2002, no 550,‎ , p. 125–140 (DOI 10.1515/crll.2002.069, CiteSeerx 10.1.1.6.9538).
  9. ↑ http://bib.univ-oeb.dz:8080/jspui/bitstream/123456789/8804/1/MemoireMeradHadjer.pdf
Si ce bandeau n'est plus pertinent, retirez-le. Cliquez ici pour en savoir plus.
Si ce bandeau n'est plus pertinent, retirez-le. Cliquez ici pour en savoir plus.

Certaines informations figurant dans cet article ou cette section devraient être mieux reliées aux sources mentionnées dans les sections « Bibliographie », « Sources » ou « Liens externes » ().

Vous pouvez améliorer la vérifiabilité en associant ces informations à des références à l'aide d'appels de notes.

  • « The On-Line Encyclopedia of Integer Sequences (OEIS) », The OEIS Foundation Inc. (vois les suites OEIS A001497, OEIS A001498 et OEIS A104548)
  • (en) Leonard Carlitz, « A Note on the Bessel Polynomials », Duke Math. J., vol. 24, no 2,‎ , p. 151–162 (DOI 10.1215/S0012-7094-57-02421-3, MR 0085360)
  • (en) H. Fakhri et A. Chenaghlou, « Ladder operators and recursion relations for the associated Bessel polynomials », Physics Letters A, vol. 358, nos 5–6,‎ , p. 345–353 (DOI 10.1016/j.physleta.2006.05.070, Bibcode 2006PhLA..358..345F)
  • S. Roman, The Umbral Calculus (The Bessel Polynomials §4.1.7), New York, Academic Press, (ISBN 978-0-486-44139-9)

Liens externes

  • icône décorative Portail de l'analyse