REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE L'ENSEIGNEMENT SUPERIEUR

ET DE LA RECHERCHE SCIENTIFIQUE
UNIVERSITE LARBI BEN MHIDI D'OUM El BOUAGHI

FACULTE DES SCIENCES EXACTES ET SCIENCES DE LA NATURE ET DE LA VIE

DEPARTEMENT DE MATHEMATIQUES ET INFORMATIQUE
N°d’ordre :……/2019

MEMOIRE
Présenté pour obtenir le diplôme de Master en Mathématiques
Option : Mathématiques appliquées

THEME
FONCTIONS DE BESSEL ET APPLICATIONS A LA
RESOLUTION D 'UN PROBLEME DE PHYSIQUE
MATHEMATIQUE

Présenté par: MERAD Hadjer

Soutenu le : 08/06/2019

Devant le Jury:
Président : Pr N. Merazga [Link] El Bouaghi
Rapporteur : Dr I. Boussefsaf [Link] El Bouaghi
Examinateurs : Dr B. Hadjou [Link] El Bouaghi
Dr M. Saadi [Link] El Bouaghi

2018/2019
Table des matières

INTRODUCTION 6

1 NOTIONS PRELIMINAIRES 8
1.1 Quelques fonctions spéciales : ; et p Fq . . . . . . . . . . . . . . . . . . 8
1.1.1 Fonction Gamma . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.1.2 Fonction Bêta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.1.3 Fonctions Hypergéométriques . . . . . . . . . . . . . . . . . . . . . 13
1.2 Polynômes orthogonaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3 Formule de Rodrigues . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.4 Fonction génératrice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.5 Méthode de Frobenius . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

2 FONCTIONS ET POLYNOMES DE BESSEL 19

I Fonctions de Bessel 20
2.1 Fonctions de Bessel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.1.1 Résolusion de l’équation de Bessel . . . . . . . . . . . . . . . . . . . 21
2.1.2 Relations de récurrence . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.1.3 Forme intégrale des fonctions de Bessel . . . . . . . . . . . . . . . . 36
2.1.4 Fonction génératrice . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.2 Fonctions de Bessel Modi…ées . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.2.1 Relations de récurrence . . . . . . . . . . . . . . . . . . . . . . . . . 43
2.2.2 Forme intégrale des fonctions de Bessel modi…ées . . . . . . . . . . 44

1
 2.3 Autres types de fonctions de Bessel . . . . . . . . . . . . . . . . . . . . . . 47
2.3.1 Fonctions de Hankel . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.3.2 Fonction de Kelvin . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
2.4 Fonctions de Bessel sphériques . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.5 Comportement des fonctions de Bessel dans les cas limites . . . . . . . . . 50

II Polynômes de Bessel 51
2.6 Polynômes de Bessel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
2.6.1 Résolusion de l’équation de polynômes de Bessel . . . . . . . . . . . 52
2.6.2 Relation entre les polynômes et les fonctions de Bessel . . . . . . . 53
2.6.3 Orthogonalité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.6.4 Relations de récurrence . . . . . . . . . . . . . . . . . . . . . . . . 56
2.6.5 Formule de Rodrigues . . . . . . . . . . . . . . . . . . . . . . . . . 56
2.6.6 Fonction génératrice . . . . . . . . . . . . . . . . . . . . . . . . . . 56
2.7 Polynômes de Bessel généralisés . . . . . . . . . . . . . . . . . . . . . . . . 57
2.7.1 Formule de Rodrigues pour les polynômes de Bessel généralisés . . . 58
2.7.2 Relations de récurrence pour les polynômes de Bessel généralisés . . 58
2.7.3 Orthogonalité des polynômes de Bessel généralisés . . . . . . . . . . 59

3 APPLICATIONS 60
3.1 Equation de la chaleur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
3.2 Equation d’Helmholtz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
3.3 Equation du pendule simple . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3.4 Equation d’onde en coordonnées sphériques . . . . . . . . . . . . . . . . . . 64

Conclusion 67

Bibliographie 69

2
 Résumé

Ce mémoire est consacré à l’étude des fonctions et des polynômes de Bessel avec
quelques applications dans le domaine de la physique .
En premier lieu, nous avons exposé les dé…nitions et les propriétés essentielles in-
tervenant au fur et à mesure dans notre travail.
Ensuite nous avons étudié profondément les fonctions et les polynômes de Bessel
qui sont solutions de l’équation di¤érentielle de Bessel de second ordre.
En…n nous avons présenté quelques applications dans domaine de la physique pour
illustrer notre études sur les fonctions et polynômes de Bessel en question.
Mots clés :Fonctions de Bessel, Polynômes de Bessel, Relation de Récurrence,
Fonction génératrice.

3
 Abstract

This dissertation is devoted to studying Bessel functions and polynomials with some
applications in physics domain.
First, we exposed the de…nitions and essential properties used in this work.
Then we studied Bessel’s functions and polynomials, which are solutions of Bessel’s
second-order di¤erential equation.
Finally, we presented some applications in physics …elds to illustrate our studies.
Keywords: Bessel Functions, Bessel Polynomials, Recurrence Relation, Generat-
ing Function.

4
 REMERCIEMENTS

La réalisation de ce mémoire a été possible grâce au concours de plusieurs personnes

à qui je voudrais témoigner toute ma reconnaissance.
Je voudrais tout d’abord adresser toute ma gratitude à mon directeur de ce mémoire
Professeurs [Link].
Je désire aussi remercier très sincèrement les professeurs N. Merazga, B. Hadjou
et M. Saadi d’avoir accepté de juger ce travail.
En…n, je tiens également à remercier tous les professeurs qui m’ont enseigné durant
ces cinq années de formation.

5
INTRODUCTION

La théorie des fonctions spéciales est née dans les travaux de plusieurs mathémati-
ciens tels que : Euler, Gauss, Laplace, Jacobi, Riemann, Tchébychev et Bessel..., elle
est depuis longtemps une discipline classique des mathématiques, qui s’enracine en
analyse mathématique et qui possède un grand champ d’applications. Beaucoup de ces
fonctions spéciales sont liées à des questions de physique mathématique et jouent un
rôle important pour fournir des solutions de certains problèmes fondamentaux de la
physique dont la dynamique est gouvernée par des équations di¤érentielles du second
ordre à coe¢ cients variables, c’est le cas notamment de la fonctions d’Airy pour le
potentiel linéaire, la fonction d’Hermite pour l’oscillateur harmonique, la fonction de
Laguerre pour l’oscillateur de Morse , la fonction hypergéométrique pour le cas de
l’atome d’hydrogène , la fonction de Jacobi pour le potentiel de Pöschl–Teller [1] ,
la fonction de Bessel pour l’équation de la di¤usion de la chaleur dans les cylindres,
oscillation d’une membrane cylindrique, propagation des ondes, les fonction de Heun
pour les systèmes ayant des massas variables et en…n les fonctions de Romanoveski
pour le potentiel trigonométrique de Rosen–Morse en mécanique supersymétrique [2]
etc.
Dans ce mémoire on s’intéresse uniquement aux fonctions et polynômes de Bessel
qui sont un cas particulier de la grande famille des fonctions spéciales. Comme il est
connu que leurs découvertes, reviennent à l’époque du mathématicien suisse Daniel
Bernoulli et portent le nom du mathématicien allemand Friedrich Wilhelm Bessel. En
1816, Bessel a développé l’analyse de ces fonctions dans le cadre de ses études du
mouvement des planètes induit par l’interaction gravitationnelle , comme étant une
généralisation des travaux antérieurs de Bernoulli. Par ailleurs, ces fonctions de Bessel

6
INTRODUCTION

ont connu de très nombreuses applications ,envahissantes plusieurs domaines de la

physique mathématique et n’ont cessé de se développer à ce jour d’aujourd’hui, de la
physique classique macroscopique à la physique moderne microscopique, nous trouvons
par exemple le cas du propagateur de l’équation de Dirac [3] dont la solution
s’exprime en terme de ces fonctions de Bessel, la solution du propagateur libre de
l’équation de Feshbach-villars [4], la particule libre de Schro dinger à trois dimensions
en representation sphérique etc......et en…n, l’extension des fonctions de Bessel déformée
a récemment vu le jour dans le cadre de la théorie des dérivées fractionnaires [5] et la
dérivée de Dunkel [6]. En conséquence, notre but de ce mémoire est de comprendre et
maitriser ces fonctions et polynômes de Bessel .
Ce mémoire est alors articulé principalement de trois chapitres.
Dans le premier chapitre, nous donnons quelques dé…nitions essentielles et quelques
propriétés utiles que nous utiliserons par la suite. Le deuxième chapitre est divisé en
deux parties, dans la première partie, on dé…nit les fonctions de Bessel avec quelques
propriétés essentielles : comme les relations de récurrence, représentation intégrale,
fonction génératrice, fonctions de Bessel de première et de deuxième espèce, fonctions de
Bessel modi…ées et sphériques. On donnera aussi une étude et un intérêt particulier aux
fonctions de Bessel d’indice entier et demi entier. Dans la deuxième partie on dé…nit les
polynômes de Bessel et on expose quelques unes de leurs propriétés importantes, aussi
que leurs relation avec les fonctions de Bessel d’ordre demi entier. Dans le troisième
chapitre, pour illustrer notre études sur les fonctions et polynomes de Bessel, on donne
quelques applications dans le domaine de la physique telles que: le pendule simple,
l’équation de la chaleur, l’équation d’Helmholtz et l’équation d’onde en représentation
sphérique.

7
Chapitre 1

NOTIONS PRELIMINAIRES

1.1 Quelques fonctions spéciales : ; et p Fq

Dans cette section, on exposera quelques dé…nitions et rappels utiles pour la suite
de notre étude.

1.1.1 Fonction Gamma

En mathématiques, la fonction gamma, notée par la lettre grecque (x); est une fonc-
tion réelle, considérée également comme une fonction spéciale. Historiquement, c’est
Euler qui a introduit et développé cette fonction. Elle prolonge la fonction factorielle
à l’ensemble des nombres réels (à l’exception des entiers négatifs) et est dé…nie par :
Z +1
(x) = e t tx 1 dt , x > 0: (1.1.1)
0

En e¤et pour tout x 2 R, la fonction

]0; +1[ ! R

t ! e t tx 1

est positive et continue sur ]0; +1[, donc elle est localement intégrable sur ]0; +1[:
L’integrale de ([Link]) converge si et seulement si les intégrales ci dessous

8
NOTIONS PRELIMINAIRES

Z 1 Z +1
t x 1
e t dt et e t tx 1 dt
0 1

convergent .

Pour 0 < t < 1 , on a

e t tx 1
tx 1

Z 1
tx 1 1
L’intégrale tx 1 dt = x 0
= x
converge si et seulement si x > 0 et d’apres le
0
critere d’équivalence on obtient la convergence de

Z 1
e t tx 1 dt < 1,
0

pour tout x > 0

t
1
Au voisinage de +1; et pour tout x > 0; on a e t tx 1
=o t2
car lim t2 e 2 tx 1
= 0:
t!1

En déduit la convergence de l’intégrale dans ([Link]) :

Proprietes de la fonction Gamma.

On introduit ici quelques propriètés de la fonction gamma ;

1. On
(x + 1) = x (x): (1.1.2)

Pour montrer cette proprièté, on integre par parties et on obtient pour x > 0
Z b Z
t x 1 e b b x e a ax 1 b t x
e t dt = + e t dt:
a x x x a

Passant à la limite, il vient

Z b Z
t x 1 1 b tx 1
(x) = lim e t dt = lim e t dt = (x + 1)
b!+1 a b!+1 x a x
a!0 a!0

Pour x entier, il est facile d’avoir

9
NOTIONS PRELIMINAIRES

2. Pour x entier, il est facile d’avoir

(n + 1) = n!:

On a pour n = 0

Z +1
(1) = e t dt = 1
0

et d’apres ([Link]), on a

(2) = 1 (1) = 1 = 1!

(3) = 2 (2) = 2 = 2!

(4) = 3 (3) = 6 = 3!

On déduit la formule demandée

3. On a la propriété
1 (2n)! p
(n + ) = 2n ; 8n 2 N: (1.1.3)
2 2 n!

Faisons un changement de variable t = y2 dans la dé…nition de la fonction

gamma ([Link]) ; on obtient
Z +1
y 2 2x 1
(x) = 2 e y dy:
0

En particulier pour x = 21 ; on aura l’intégrale d’une fonction Gaussienne

Z +1
1 y2
p
( )=2 e = :
2 0

D’apres la proprièté ([Link]) scussissivement , on obtient

1 1 1
(n + 1 ) = (n ) (n )
2 2 2
..
.
1 3 31 1
= (n )(n )::: ( );
2 2 22 2

10
NOTIONS PRELIMINAIRES

qui s’écrit aussi comme suit

1 (2n 1)(2n 3):::3 1p
(n + ) = :
2 2n
Multiplions et divisons par (2n)(2n 2):::2, il vient
1 (2n)(2n 1)(2n 2)(2n 3):::3 2 1p
(n + ) = :
2 (2n)(2n 2):::2 2n
D’où le résultat
1 (2n)! p
(n + ) = 2n :
2 2 n!
Remarque 1.1.1 On montre aussi
p
3 3
n! (n + ) = (n + 1) (n + ) = 2n 1
(2n + 1)! (1.1.4)
2 2 2
Remarque 1.1.2 La formule de duplication de Legendre [11]
22x 1
1
(2x) = p (x) (x + ):
2
4. Fonction gamma pour les valeurs négatives

De la formule ([Link]), il est facile de montrer la fonction gamma est dé…nie aussi
pour x < 0: Pour ce faire, on voit que la partie gauche de la relation est bien dé…ni
pour x + 1 > 0 c-à-d pour x > 1; donc la formule ([Link]) est dé…nie pour x > 1
et ce processus est répété pour dé…nir (x) pour x < 0:

4.1. ( m) = 1 8m 2 N :

(x+1)
De ([Link]) , on écrit que (x) = x
, on remarque que : (0) = 1; alors on
déduit

( 1 + 1) = 1 ( 1) = 1 ) ( 2 + 1) = 2 ( 2) = [Link]

4.2. Pour n 2 N; on a
1 1
( )= ( + 1);
2 2
utilison la formule ([Link]), il vient
1 1 1 1 1
( + 1) = ( )) ( )= 2 ( );
2 2 2 2 2

11
NOTIONS PRELIMINAIRES

et d’aprés ([Link]) ; on trouve

1 p
( )= 2 ;
2
scussissivement on obtient

1 ( 1)n 22n n! p
( n+ )= :
2 (2n)!

1.1.2 Fonction Bêta

La fonction bêta est un type d’intégrale d’Euler dé…nie pour tous nombres réels x et
y strictement positifs, notée (x; y). Historiquement, elle a été étudiée par Euler et
Legendre et doit son nom à Jacques Binet, Elle est dé…nie par :

Z 1
(x; y) = tx 1 (1 t)y 1 dt , x > 0 et y > 0: (1.1.5)
0

Considérons la fonction

f : t ! tx 1 (1 t)y 1 :

Cette fonction est positive et continue sur l’intervalle ]0; 1[:

Au voisinage de 0, on a

tx 1 (1 t)y 1
tx 1

et d’après le critère de Riemann, la fonction de signe constant tx 1 est intégrable ssi

x 1> 1 c-à-d x > 0:

Au voisinage de 1, on a

tx 1 (1 t)y 1
(1 t)y 1

et d’après le critère de Riemann, la fonction de signe constant (1 t)y 1 est intégrable

ssi y > 0:

12
NOTIONS PRELIMINAIRES

Proprietes de la fonction Bêta

1. La fonction Bêta est symétrique, juste par le changement de variable u = 1 t,

on otient directement
(x; y) = (y; x):

2. On peut s’écrire en fonction des fonctions Gamma selon la formule suivante

(x) (y)
(x; y) = : (1.1.6)
(x + y)

Remarque 1.1.3 Une autre propriète de la fonction Gamma, donnée par

(x) (1 x) = ; pour tout 0 < x < 1: (1.1.7)

sin x

Preuve. Pour montrer cette dernière proprièté([Link]), on utilisera l’intégrale

suivante[11] Z 1
zx 1
dz = : (1.1.8)
0 1+z sin x
z y
En posant 1+z
=y ou z = 1 y
; l’integral précedente ([Link]) devient
Z 1
= y x 1 (1 y) x dy;
0

qui est égale à la fonction Béta

= (x; 1 x)

et en…n, utilisons la proprièté ([Link]) ; on trouve

(x) (1 x)
= = :
(x + 1 x) sin x

1.1.3 Fonctions Hypergéométriques

Les fonctions hypergéométriques sont exprimables en séries entières qui ont une forme
un peu spéciale, elles permettent également d’exprimer des fonctions spéciales et des

13
NOTIONS PRELIMINAIRES

polynômes orthogonaux comme ceux de Bessel, de Legendre, Jacobi ou de Hermite.

Nous dé…nissons la fonction hypergéométrique généralisée p Fq dépendant de plusieurs
paramètres i et j de la manière suivante :

2 3
1; 2 ; :::; p; X
1
( 1 )n ( 2 )n ; :::( p )n x
n
p Fq ( 1; 2 ; :::; p; 1; 2 ; :::; q ; x) = :p Fq
4 x5 =
( 1 )n ( 2 )n ::( q )n n!
1; 2 ; :::; q; n=0

où les ( )n sont donnés par

( )n = ( + 1):::( + n 1)
( + n)
( )n =
( )

représentent le symbole de Pochhammer.

Si i est égal à zéro ou un nombre entier négatif alors la série se termine, c-à-d p Fq
devient un polynôme.
En particulier citons des types celèbres de fonction hypergéométrique :

si p = 1 et q = 0; alors

X (a)n
1 F0 (a; x) = xn
n>0
n!
est une solution de l’équation di¤érentielle suivante :

d
(1 x) a 1 F0 (a; x) = 0:
dx

Si p = 2 et q = 0; alors
X (a)n (b)n
2 F0 (a; b; x) = xn
n>0
n!

est une solution de l’équation di¤érentielle suivante :

d22 d
x 2
+ [ 1 + (1 + b + a)x] + ab 2 F0 (a; b; :; x) = 0:
dx dx

14
NOTIONS PRELIMINAIRES

Si p = 2 et q = 1; alors

X (a)n (b)n
2 F1 (a; b; c; x) = xn
n>0
n!(c)n
est une solution de l’équation di¤érentielle :

d2 d
x(1 x) + [c (a + b + 1)x] ab 2 F1 (a; b; c; x) = 0:
dx2 dx

1.2 Polynômes orthogonaux

La théorie des polynôme orthogonaux fait partie des branches de l’analyse. Elle
a connu une grande application dans plusieurs domaines de la physique , en
particulier en mécanique quantique et n’a cessé de se développer jusqu’à ce jour.
les polynômes orthogonaux sont dotés d’une fonction (x) dite fonction poid, cette
fonction est dé…nie et intégrable sur [a; b] :

Dé…nition 1.2.1 on appelle fonction poid sur l’intervalle [a; b] ; toute fonction (x)
véri…ant les deux conditions
-la positivité, (x) > 0 8x 2 [a; b] ;
Z b
-la convergence, (x)dx < 1:
a

Dé…nition 1.2.2 Soit fPn (x)gn une suite de polynôme avec deg [ Pn (x)] = n et
soit (x) une fonction poid sur [a; b] :On dit que fPn (x)gn est une suite de
polynômes orthogonaux par rapport à (x) dans [a; b] si
Z b
Pn (x)Pm (x) (x)dx = 0 si n 6= m:
a

1.3 Formule de Rodrigues

Considérons l’équation di¤érentielle du second ordre homogène à coe¢ cients variables
suivante :
p(x)y 00 + q(x)y 0 + y = 0 (1.3.9)

15
NOTIONS PRELIMINAIRES

où p(x) et q(x) sont des polynômes tels que deg(p(x)) 2 et deg(q(x)) 1 et

est une constante.
On dé…nit la fonction w(x) par la formule suivante

Zx
1 q(t)
w(x) = exp( dt):
p(x) p(t)
x0

Alors la suite des polynômes orthogonaux , solutions l’équation di¤érentielle ([Link])

est donnée par la formule suivante dite formule de Rodrigues (anciennement appelée
formule de Ivory-Jacobi) :

cn d n
Pn (x) = (w(x)(p(x))n )
w(x) dxn
où cn est une constante de normalisation, que sera déterminée par la condition suiv-
ante Z a
Pn (x)Pn (x) (x)dx = 1:
b

Voir [16] :

1.4 Fonction génératrice

Soit une fonction G(x; t) développable en série entière de t dans un certain domaine D
comme suit :
X
1
n
G(x; t) = n (x)t :
n= 1

On dit que G(x; t) est la fonction génératrice de la suite de fonctions f n (x)gn .

1.5 Méthode de Frobenius

En analyse, la méthode de Frobenius, du nom du mathématicien allemand Ferdinand
Georg Frobenius, est une technique d’obtention du développement en série des solutions
d’une équation di¤é[Link] se focalise ici sur le cas des équations linéaires du second
ordre, homogène suivante :

16
NOTIONS PRELIMINAIRES

2
2dy dy
x + xp(x) + q(x)y = 0 (1.5.10)
dx2 dx
où p(x) et q(x) sont des fonctions développables en série entière ou des polynômes
donnés par les séries suivantes:
X
1 X
1
p(x) = pm xm ; q(x) = q m xm : (1.5.11)
m=0 m=0

Cherchons les solutions de cette équation di¤érentielle([Link]) sous forme de séries

X
1
y(x) = an xs+n avec a0 6= 0 (1.5.12)
n=0

où s est un nombre réel, appelé l’indice de la série.

Maintenant, nous étudions les cas des points suivants
1) On dit que le point x = 0 est un point régulier si les limites suivantes
8
< lim p(x)
x!0
: lim q(x)
x!0

existent et sont …nies, sinon le point x = 0 est dit singulier.

2) Si le point x = 0 est un point singulier, et les limites ci-dessous existent et sont
…nies 8
< lim xp(x)
x!0
: lim x2 q(x)
x!0

alors le point x = 0 est dit point singulier régulier, sinon il est dit point singulier
irrégulier.

Théorème 1.5.1 Fuchs

Si le point d’expansion d’une série est un point régulier et au pire des cas un point
singulier régulier, alors la méthode de Frobenius nous fournit au moins une solution et
la nature des solutions obéit à l’équation indicielle suivante

a0 s2 + (p0 1)s + q0 = 0 avec a0 6= 0 (1.5.13)

17
NOTIONS PRELIMINAIRES

Remarque 1.5.1 Cette équation est déduite par l’insertion de ([Link]) ; ([Link]) dans
([Link]) au point x = 0:

Pour obtenir la forme de la 2eme solution, qui est dépendante de la di¤érence des
racines s1 , s2 de ([Link]) , nous distinguons les 3 cas suivants.

Si les indices sont égaux, alors on obtient une solution.

Si les indices sont di¤érents par un nombre non entier, on peut obtenir deux solutions
ys1 et ys2 :

Si les indices sont di¤érents par un nombre entier, alors l’indice le plus grand nous
fournit une solution.

Remarque 1.5.2 Si la méthode de Frobenius nous fournit une seule solution alors la
2eme solution peut être obtenir via cette formule
Z x Ru
exp( P (v)dv)
y2 (x) = y1 (x) du;
y12 (u)

où y1 est la première solution de l’équation ([Link]) .

Voir [7] [17]

18
Chapitre 2

FONCTIONS ET POLYNOMES DE
BESSEL

19
 Partie I

Fonctions de Bessel

20
Fonctions de Bessel

2.1 Fonctions de Bessel

2.1.1 Résolusion de l’équation de Bessel

Historiquement, les fonctions de Bessel font partie principale de la famille de fonctions
spéciales de la physique mathématique, découvertes avant par le mathématicien Daniel
Bernoulli. En 1816, le mathématicien Friedrich Wilhelm Bessel a développé l’analyse
de ces fonctions dans le cadre de ses études sur les mouvements des planètes induits
par la théorie de la gravitation et a présenté une généralisation et une extension des
travaux de recherche de Bernoulli.
Soit une équation di¤érentielle de second ordre à coe¢ cients variables suivante :
d2 y dy
2
+ P (x) + Q(x)y = 0 (2.1.1)
dx dx
Nous rappelons que de nombreuses fonctions spéciales découlent de cette équation
([Link]), mais dans notre étude, nous nous limiterons aux équations di¤érentielles
linéaires suivante :
d2 y dy
x2 + x + (x2 2
)y = 0 (2.1.2)
dx2 dx
dont les solutions sont des fonctions de Bessel d’ordre ; où x est la variable, y(x) la
fonction inconnue et un paramètre réel.
Nous constatons que le point x = 0 est un point singulier régulier. En e¤et
8
< P (x) = x = 1
x2 x
(2.1.3)
: Q(x) = 2 x 2 2
x

et leurs limites véri…ent

8
< lim xP (x) = 1;
x!0
(2.1.4)
: lim x2 Q(x) = 2
:
x!0

A…n de résoudre cette dernière équation ([Link]), nous utilisons la méthode de

Frobenius, qui est une technique de résolution des équations di¤érentielles linéaires,
basée sur le développement en série de puissances dans x, sous la forme suivante
X
y(x) = ak xk+ ; x > 0: (2.1.5)
k>0

21
Fonctions de Bessel

Par un calcul direct, on aura,

X X
y(x) = (k + )ak xk+ 1
et y
½ (x) = (k + )(k + 1)ak xk+ 2
: (2.1.6)
k>0 k>0

Nous remplaçons dans l’équation ([Link]) on trouve :

X
2 2 2 2 +1
( )a0 x +(2 + +1 )a1 x + (k + )2 2
)ak + ak 2 xk+ = 0 (2.1.7)
k>2

et pour que y satisfasse l’équation di¤érentielle ([Link]), il faut que le coe¢ cient
de chaque puissance de x dans ([Link]) soit égal à zéro, on trouve :
8
< a0 ( 2 2
)=0
(2.1.8)
: a1 (2 + 2 + 1 2
) = 0:

Comme a0 6= 0, de la première relation ([Link]), on déduit = , puis on insère

dans la deuxième relation on aura a1 ( 2v + 1) = 0 qui s’annulera pour a1 = 0
1
ou = 2
:

Remarque 2.1.1 Il est remarquable que l’équation di¤érentielle ([Link]) reste invariante
si on change par , par conséquent, par la suite, on prendra juste > 0:

On trouve:
8
< a( 2 2
)=0
0
) = (2.1.9)
: a0 6= 0

2 2
a1 (2 + +1 )=0 (2.1.10)

1 cas si
er
= , > 0: la deuxième relation ([Link]) devient

a1 (2 + 1) = 0

ce qui implique que a1 = 0; et de la relation de récurrence apres quelques

réarrangements de ([Link])

((k + )2 2
)ak + ak 2 = 0; k > 2

22
Fonctions de Bessel

il vient, tous les coe¢ cients impairs de la série sont nuls.

En e¤et

((k + )2 2
)ak + ak 2 = 0; 8k > 2 ) ak = ak 2
(k+ )2 2 8k > 2

d’où
a1 = 0; a1 = a3 = a5 = a2k+1 = : : : = 0

On déduit alors les coe¢ cients pairs comme suit :

8
>
> a2 = a0
= a0
>
< 2(2+2 ) 22 (1+ )
a2 ( 1)2 a2
a4 = 4(4+2 = (2.1.11)
>
>
) 4
2 2!(1+ )(2+ )
>
: 1)3 a0
a6 = 6(6+2 ) = 26 3!(1+( )(2+
a4
)(3+ )

et le terme général s’écrira alors:

( 1)k a0
a2k = : (2.1.12)
22k k!(1 + )(2 + ) : : : (k + )
En injectant ([Link]) dans ([Link]), on obtient la 1ere solution de l’équation sous la
forme suivante:

!
X ( 1)k x2k
J (x) = a0 x 1+ : a0 arbitraire
k>1
22k k!(1 + )(2 + ) : : : (k + )

Maintenant si on prend x < 0 dans ([Link]), dans ce cas en posant x = t ,

t > 0 et z(t) = y(x) on aura,
8
< x2 y" + xy 0 + (x2 2
)y = 0; x < 0
(2.1.13)
: t2 z" + tz 0 + (t2 2
)z = 0; t > 0

Pour solutionner cette équation di¤érentielle ([Link]), on procède par le même

démarche que la première équation di¤érentielle ([Link]), on prend la solution sous
forme de série comme suit

X
z(t) = ak tk+ (2.1.14)
k>0

23
Fonctions de Bessel

on déduit la solution

X ( 1)k t2k
z (t) = a0 t (1 + ); t > 0 (2.1.15)
k>1
22k k!(1 + )(2 + ) : : : (k + )

on revient à l’ancienne variable x, et par un simple remplacement on trouve

X ( 1)k x2k
J (x) = a0 ( x) (1 + ); x < 0 (2.1.16)
k>1
22k k!(1 + )(2 + ) : : : (k + )

On note que on peut regrouper les deux solutions pour x 2 R , l’expression de

J s’écrira :
X ( 1)k x2k
J (x) = a0 jxj (1 + ): (2.1.17)
k>1
22k k!(1 + )(2 + ) : : : (k + )

A ce niveau, pour discuter la deuxième solution, on applique le théorème de

Fuchs : on a

8 8
>
> < est entier 2 f0; 1; 2 : : :g
>
< 2 est entier
( )=2 : est demi entier 2 1 3
; ; : : : 2n+1
>
>
2 2 2
>
: 2 n’est pas un entier ( 6= 0)

Il est préférable de discuter les cas ci-dessus .

2 n’est pas un entier

2 cas
eme
= , > 0 et 2 n’est pas un entier: de la deuxiéme relation ([Link]),
il vient
a1 (1 2 )=0 (2.1.18)

ce qui implique que a1 = 0; et de la relation de récurrence, on a :

8
< a1 = a3 = a5 = a2k+1 = : : : = 0
: ((k )2 2
)ak + ak 2 = 0; k > 2

24
Fonctions de Bessel

ak 2
ak =
k(k 2 )
La 2eme solution de l’équation sera :

X ( 1)k x2k
J (x) = a0 jxj (1 + ); x 2 R , a0 arbitraire
k>1
22k k!(1 )(2 ) : : : (k )

Maintenant, dans le but d’écrire les solutions en fonction de la fonction

Eulerienne, nous nous concentrons uniquement quand x > 0

Remarque 2.1.2 On rappelle la relation suivante de la fonction Gamma ([Link])

(s + n) = (s + n 1) (s + n 1)

= (s + n 1)(s + n 2)(s + n 3):::(s + 1)s (s)

la forme de J se met alors comme suit

x X ( 1)k x 2
J (x) = ( ) ( )2k ; a0 = (2.1.19)
2 k>0 k! (1 + k + ) 2 (1 + )

et celle de J par :

x X ( 1)k x 2
J (x) = ( ) ( )2k ; a0 = (2.1.20)
2 k>0
k! (1 + k ) 2 (1 )
(1)

Théorème 2.1.1 Lorsque n’est pas un entier J (x)et J (x) sont deux solutions linéaire-
ments indépandantes(2)
(1)
si on prend 0 on commance par la 2éme cas et on decalu la 1ere cas ci dessus de 2 n’est pas
un entier
(2)
2 n’est pas un entier implique que n’est pas un entier mais la resproque contient les cas suivants
:2 est un entier ou 2 n’est pas un entier ( est un demi entier) .

25
Fonctions de Bessel

Preuve. Calculons le Wronskien de W (J (x); J (x)) :Par un calcul direct on

obtient;

J (x) J (x) 2 sin

W (J (x); J (x)) = = 6= 0 (2.1.21)
J 0 (x) J 0 (x) x

Pour démontrer ([Link]) , on a ; J (x) et J (x) satisfont ([Link])

x2 J" (x) + xJ 0 (x) + (x2 2

)J (x) = 0 (2.1.22)

x2 J" (x) + xJ 0 (x) + (x2 2

)J (x) = 0

On fait la soustraction de la première équation multipliée par J (x) de la deuxième

équation multipliée par J (x) et on divise le résultat obtenu par x; on aura

x (J (x)J" (x) J (x)J" (x)) + J (x)J 0 (x) J (x)J 0 (x) = 0

d
) x J (x)J 0 (x) J (x)J 0 (x) =0
dx
d
) (xW (J (x); J (x))) = 0
dx
C
W (J (x); J (x)) = ; où C est une constante
x
Cherchons maintenant la valeur de C. Puisque C est une constante donc on calcule
cette dernière pour x ! 0:
Pour determiner cette constante; on prend les premiers termes de J (x) et J (x)

1 x x 1
J (x) = ( ) , J 0 (x) =
(1 + ) 2 2 ( + 1)
1 x x 1
J (x) = ( ) , J 0 (x) =
(1 ) 2 2 ( + 1)

et par un simple calcul, on trouve

1
1 x x 1 x x 1
W (J (x); J (x)) = ( ) ( )
(1 + ) 2 2 ( + 1) (1 ) 2 2 ( + 1)
2 x 1 2 x 1 2x 1
= = =
(1 + ) (1 ) ( ) (1 ) ( ) (1 )

26
Fonctions de Bessel

où l’on a utilisé la relation ([Link]) de fonction Gamma suivante :

(x) (1 x) =
sin x
2 sin
W (J (x); J (x)) =
x
Il vient alors
2 sin
C= 6= 0

Remarque 2.1.3 Nous notons ici que J (x) est une série entière(1) dont le rayon de
convergence est in…ni, en appliquant le lemme de D’alembert suivant :
( 1)k x
uk = ( )2k
k! (1 + k + ) 2
uk+1 1 x
= ( )2
uk (k + 1)(1 + k + ) 2
car quand k ! 1 le rapport tend vers 0 quel que soit x.

Remarque 2.1.4 la fonction J est appelée fonction de Bessel de 1ere espéce d’order
:

2 est un entier et est un entier (2)

Théorème 2.1.2 lorsque = n est un entier, on a

J n (x) = ( 1)n Jn (x) (2.1.23)

Preuve. Si n > 0; et quand (1 + k n) 0 d’aprés les propriètés de la fonction

Eulerienne, on a (1 + k n) = 1 donc les termes associés de J n (x) sont nuls,
la série se réduira alors à
x X ( 1)k x
n
J n (x) = ( ) ( )2k
2 k>n
(k + 1) (1 + k n) 2
X ( 1)k
(1)
J (x) = ( x2 ) x 2k
k! (1+k+ ) ( 2 )
k>0
(2)
2 est un entier ) est un entier ou est un demi entier

27
Fonctions de Bessel

x X ( 1)k+n x
n
J n (x) = ( ) ( )2k+2n
2 k>0
(k + 1 + n) (1 + k n + n) 2
x X ( 1)n ( 1)k x
= ( )n ( )2k
2 k>0 (k + 1 + n) (1 + k) 2

X ( 1)k x
J n (x) = ( 1)n ( )2k+n = ( 1) n
Jn (x)
k>0
k! (n + k + 1) 2

Si n < 0; dans ce cas, nous pouvons écrire n = p, avec p > 0.

Jp (x) = ( 1)p J p (x):

On essaie de chercher une deuxième solution linéairement indépendante de J de

l’équation de Bessel .

Théorème 2.1.3 Lorsque est un entier, les deux soulutions de Bessel linéairement
indépendantes sont :

Jn (x) (2.1.24)
(cos ) J (x) J (x)
Yn (x) = lim Y (x) = lim
!n !n sin
cos n Jn (x) J n (x)
Preuve. Nous constatons que pour = n la fonction Yn (x) = sin n

n’est pas dé…nie car

8
< cos n = ( 1)n ; J (x) = ( 1)n J (x)
n n
: et sin n = 0
et pour cela, on dé…nit Yn par

Yn = lim Y
!n

dont la limite satisfait à la règle de l’Hopital a…n de lever cette indétermination,

on a :

28
Fonctions de Bessel

Yn (x) = lim Y (x)

!n
cos
J (x) J (x)
= lim
!n sin
[@=@ (cos J (x) J (x))] =n
=
[@=@ (sin )] =n
cos n [@=@ (J (x)] =n [@=@ (J (x)] =n
=
[ cos ] =n
1 @ @
= J (x) ( 1) J (x) :
@ @ =n

A ce niveau, nous devons montrer que Yn (x) est une solution de l’équation de
Bessel et elle est linéairement indépendante de Jn (x).
On a
1)Yn (x) est une solution de l’équation de Bessel
d2 J dJ
x2 2
+x + (x2 2
)J = 0 (2.1.25)
dx dx
Tout d’abord, nous dérivons l’équation di¤érentielle ([Link]) partiellement par rap-
port à , on trouve

d2 @J d @J @J
x2 2
+x + (x2 2
) 2 J =0 (2.1.26)
dx @ dx @ @
Comme J véri…e l’équation de Bessel d’ordre , nous avons donc

d2 @J d @J @J
x2 +x + (x2 2
) 2 J =0 (2.1.27)
dx2 @ dx @ @
En multipliant ([Link]) par ( 1) et on fait la soustraction des deux équations,
on aboutit à cette expression

2 d2 @J @J d @J @J @J @J
x ( ( 1) )+x ( ( 1) )+(x2 2
)( ( 1) ) 2 (J ( 1) J ) = 0:
dx2 @ @ dx @ @ @ @
Passons à la limite, quand ! n; cette dernière équation se convertira à :

d2 d 2n
x2 2
Yn (x) + x Yn (x) + (x2 n2 )Yn (x) (Jn (x) ( 1)n J n (x)) = 0:
dx dx

29
Fonctions de Bessel

Utilisons cette proprièté des fonctions de Bessel démontrée en haut ([Link]),

J n (x) = ( 1)n Jn (x)

il vient alors

d2 d
x2 2
Yn (x) + x Yn (x) + (x2 n2 )Yn (x) = 0
dx dx
Nous remarquons e¤ectivement que Yn (x) satisfait à l’équation de Bessel d’ordre
n.
2) Yn (x) est linéairement indépendante de Jn (x).
En appliquant

J (x) Y (x)
W (J (x); Y (x)) =
J 0 (x) Y 0 (x)

il vient alors

1
W (J (x); Y (x)) = cos (J (x)J 0 (x) J 0 (x)J (x)) + J 0 (x)J (x) J (x)J 0 (x)
sin

qui se simpli…e comme suit

1 2
= W (J (x); J (x)) = 6= 0
sin x

où l’on a utilisé ([Link]) et à la limite ! n; on voit clairement que W (Jn (x); Yn (x)) 6=
0, ce qui prouve l’indépendance des solutions.

Remarque 2.1.5 si n’est pas un entier dans ce cas,Y est une solution particulière de
l’équation de Bessel et s’écrit comme combinaison linéaire de deux solutions linéairement
indépendantes.

Remarque 2.1.6 Yn (x) est appelée fonction de Bessel 2eme espèce.

Théorème 2.1.4 (Expression explicite de Yn pour n entier)

30
Fonctions de Bessel

2n x o
Yn (x) = log
+ Jn (x)
2
1X
1
1 x
( 1)k [ (k) + (n + k)] ( )2k+n
k=0
(k)!(k + n)! 2

1 X1
n
x
(n k 1) ( )2k n

k=0
2

où (k) = 1 + 12 + 13 ::: + k1 ; (0) = 0 et (1)

est la constante d’Euler .

Pour n = 0

2n o 2X
1
x 1 x
Y0 (x) = log + J0 (x) ( 1)k (k) ( )2k
2 k=0
(k)!(k)! 2

2.1.2 Relations de récurrence

Théorème 2.1.5 Les relations de récurrence pour les fonctions de Bessel de première
espèce sont données comme ci-dessous quelque soit les valeurs de :

d
1. dx
[x J (x)] = x J 1 (x)

d
2. dx
[x J (x)] = x J1+ (x)

3. J 0 (x) = J 1 (x) x
J (x)

4. J 0 (x) = x J (x) J +1 (x)

5. J 0 (x) = 12 [J 1 (x) J +1 (x)]

2
6. J 1 (x) +J +1 (x) = x
J (x)

Preuve. [Link] a
X ( 1)k x
J (x) = ( )2k+
k>0
k! (1 + k + ) 2
0
(1) (1)
= (1) '0.5772156...

31
Fonctions de Bessel

alors

d d (1) X ( 1)k x
[x J (x)] = ( )2k+2
dx dx k>0 k! (1 + k + ) 2
X 2(k + )( 1)k
= x2k+2 1

k>0
k! (1 + k + )22k+2
X ( 1)k x
= x 2k+2
( )2k+ 1
=x J 1 (x)
k>0
k! (k + )2 2

2.
d X (2k)( 1)k
x J (x) = x2k 1
dx k>0
k! (1 + k + )22k+

k =r+1)
X 2(r + 1)( 1)r+1
= 2(r+1)+
x2r+1
r>0
(r + 1)! (r + 2 + )2
X ( 1)r
= x 2r+1+
x2r+1+ = x J +1 (x)
r>0
r! (r + 2 + )2

3.
d 1
[x J (x)] = x J (x) + x J 0 (x) = x J 1 (x)
dx
et nous divisons par x :
[Link] e¤ectuons la dérivation dans le résultat (2) de [T heoreme [Link] en multi-
pliant par x :
[Link].
[Link].

Remarque 2.1.7 Ces relations [T heoreme [Link] restent valables pour la fonction Y (x)

Preuve. lorsque est non entier

de (1) et (2) du Théorème 2.1.4

32
Fonctions de Bessel

d d cos J (x) J (x)

[x Y (x)] = x
dx dx sin
[cos x J 1 (x) ( x J +1 (x))]
=
sin
x [cos f( 1) + g J 1 (x) + J +1 (x)]
=
sin f( 1) + g
x [ cos(( 1) )J 1 (x) + J +1 (x)]
=
sin(( 1) )
= x Y 1 (x)

si est un entier on prend la limite seulement.

33
Fonctions de Bessel

2 est un entier , est un demi-entier

2n+1
Théorème 2.1.6 lorsque 2 est un entier et est un demi-entier = 2
; n2N ,
les deux solutions de Bessel linéairements indépendantes sont:

r
n x 2n+1 1 d n sin x
J 2n+1 = ( 1)n 2 ( ) 2 ( ) (2.1.28)
2 2 x dx x
r
n x 2n+1 1 d n cos x
J 2n+1 = 2 ( ) 2 ( )
2 2 x dx x
1
Preuve. Remplaçons = 2
dans la 1ere solution ([Link])

x X ( 1)k x
J (x) = ( ) ( )2k
2 k>0 k! (1 + k + ) 2
on trouve

X 1 x 2k+ 1
J1 = ( 1)k 3 ( )
2 (2.1.29)
2
k>0
k! (k + 2 ) 2

Puis utilisons la proprièté ([Link]) dans ([Link])

r r
2 X ( 1)k 2k+1 2
J1 = x = sin x
2 x k>0 (2k + 1)! x

De même pour J 1 (x), on aura

2

X 1 x 2k 1
J 1 = ( 1)k 1 ( )
2
2
k>0
k! (k + 2 ) 2

en utilisant la propriété de la fonction Gamma ([Link]),

Il vient alors
X 2k 2k+ 21
k2 2 1
J 1 = ( 1) p x2k 2
2
k>0
(2k)!

qui peut s’écrire

r r
2 X 1 2k 2
= ( 1)k x = cos x
x k>0 (2k)! x

34
Fonctions de Bessel

A ce stade, cherchons la forme générale Jn+ 1 par récurrence à partir de J 1 , posons

2 2

1
n
x 2 Jn+ 1 = Un
2
3
n
x 2 Jn+ 3 = Un+1
2

et d’après la relations de récurrence (2) du [Théorème 2.1.4] pour les fonctions de

Bessel

d
x J = x J +1 (x);
dx
alors

dUn d h n 1
i
n 1
= x 2 Jn+ 1 = x 2 Jn+ 3 (x) = xUn+1
dx dx 2 2

on déduit la relation de récurrence suivante

1 d
Un+1 = Un :
x dx

Nous calculons les termes U1 ; U2 ::::

1 d
U1 = U0
x dx
1 d 2
U2 = ( 1)2 ( ) U0
x dx
1 d 3
U3 = ( 1)3 ( ) U0
x dx
..
.

d’où le terme général s’écrit :

1 d n
Un = ( 1)n ( ) U0 ;
x dx
q
1 1 2
où U0 = x 2 J1 = x 2
x
sin x, alors on obtient
2

r
n n x 2n+1 1 d n sin x
J 2n+1 = ( 1) ( ) 2 ( )
2 2 x dx x
Pour J 2n+1 on suit les mêmes démarches et on montre que les deux solutions sont
2

lineairements indépendantes

35
Fonctions de Bessel

Jn+ 1 (x) J n 1 (x) 2 sin 2n+1

2 2
W (J (x); J (x)) = 2 2
= = cos n
0
Jn+ 1 (x) J0 n 1 (x) x x
2 2

Remarque 2.1.8 J 2n+1 et J 2n+1 sont appelées Fonctions de Bessel d’indice demi-entier.
2 2

Remarque 2.1.9 On peut montrer que Y 2n+1 (x) et J 2n+1 sont dépendantes ,
2 2

Y 2n+1 (x) = ( 1)n+1 J 2n+1

2 2

En utilisant la dé…ntion ([Link])

cos(n + 12 ) Jn+ 1 (x) J n 1 (x)

2 2
Yn+ 1 (x) = 1
2 sin(n + 2
)

et par un développement des fonction cos(n + 12 ) et sin(n + 12 ), on obtient

J n 1 (x)
2
= ( 1)n+1 J 2n+1
( 1)n 2

2.1.3 Forme intégrale des fonctions de Bessel

Théorème 2.1.7 La forme intégrale des fonctions de Bessel est donnée par :
Z 1
( 21 x) 1 1
J (x) = p (1 t2 ) 2 eixt dt; > (2.1.30)
( + 1=2) 1 2
1
Pour assurer la convergence de l’intégrale nous mettons > 2

Preuve. Considérons l’intégrale I de…nie par

Z 1
1
I = (1 t2 ) 2 eixt dt (2.1.31)
1
Z 1 X (ixt)k
1
= (1 t2 ) 2 dt
1 k>0
k!
X (ix)k Z 1
1
= (1 t2 ) 2 tk dt
k>0
k! 1

36
Fonctions de Bessel

Comme l’intégrale Z 1
1
(1 t2 ) 2 tk dt
1

est nulle pour k impair, alors k doit être pair, on prend k = 2s; l’intégrale s’écrira :
Z 1 Z 1
1 1
2 k
(1 t ) 2 t dt = (1 t2 ) 2 t2s dt
1 1
Z 1
1
= 2 (1 t2 ) 2 t2s dt
0

on fait un changement de variable tel que u = t2 , on obtient :

Z 1
1 1
= (1 u) 2 us 2 du
0

qui est la fonction Bêta, et on trouve

Z 1
1 1 1 1 1 1
(1 u) 2 us 2 du = ( + ; s + ); (s > car k > 0) , >
0 2 2 2 2
1
( + 2 ) (s + 12 )
=
( + s + 1)
Remplaçons cette intégrale dans ([Link]) ; il vient,

X (ix)2s ( + 1 ) (s + 1 )
2 2
I =
s>0
2s! ( + s + 1)
1 X (i2 )s (x)2s (s + 21 )
= ( + )
2 s>0 2s! ( + s + 1)
1 X ( 1)s (x)2s (2s)! p
= ( + )
2 s>0 2s! ( + s + 1)22s s!
1 p x X ( 1)s ( x )2s+
2
= ( + )( )( )
2 2 s>0
( + s + 1)s!

qui s’écrit encore comme,

1 p x
I = ( + )( )( ) J (x)
2 2
où l’on a utilisé la dé…nition de la fonction J (x) ([Link]).
Donc la forme intégrale de J (x) sera

37
Fonctions de Bessel

1 x
J (x) = 1 p ( ) I
( + 2
)( ) 2

2.1.4 Fonction génératrice

La fonction génératrice pour les fonctions de Bessel est donnée par la formule suivante
:

1 1 X
+1
G(x; t) = exp x t = tn Jn (x)
2 t n= 1
1 1
Preuve. Nous développons exp 2
x t t
en puissances de t

1 1 1 1x
exp x t = exp xt : exp
2 t 2 2t
r +1 1x m
X
+1 1
xt X 2 2 t
=
r=0
r! m=0
m!
X
+1 1 r+m
xr+m tr m
( 1)m 2
=
r;m=0
r!m!

Maintenant on prend n = r m

X
+1 X
+1
x n+2m tn X
+1
( 1)m = tn Jn (x)
n= 1 m=0
2 m! (n + m)! n= 1

38
Fonctions de Bessel

2.2 Fonctions de Bessel Modi…ées

Dans cette section, on s’intéresse à d’autre type des fonctions de Bessel ,

Théorème 2.2.1 Soit l’equation suivante

d2 y dy
x2 2
+ x + (( x)2 2
)y = 0 (2.2.32)
dx dx
est une équation réductible à l’équation de Bessel où un paramètre complexe.
Preuve. par le changement de variable suivant: t = x; on a
dy
dx
= dy
dt
d2 y 2 d2 y
dx2
= dx2
l’équation ([Link]) se transforme à
d2 y dy
t2 2
+ t + (t2 2
)y = 0:
dt dt
qui est l’equation de Bessel dont sa solution générale formelle est donnée par:

y = AJ (t) + BY (t)

ou avec l’ancienne variable

y = AJ ( x) + BY ( x):

Maintenant, cherchons la solution de l’equation suivante

d2 y dy
x2 2
+x (x2 + 2 )y = 0 (2.2.33)
dx dx
Nous remarquons que c ’est un cas particulier de l’expression ([Link]) juste en mettant
= i; alors sa solution est

y = AJ (ix) + BY (ix)

Remarque 2.2.1 Pour éviter les arguments complexes, il est préférable d’introduire
d’autre fonction réelle de Bessel dite fonction de Bessel modi…ée de première espèce
I (x) dé…nie par:
I (x) = (i) J (ix);

39
Fonctions de Bessel

qui est une solution de l’équation ([Link]).

Preuve. Montrons que I (x) est réelle
X ( 1)k ix
(i) J (ix) = i ( )2k+
n>0
n! (1 + k + ) 2
X ( 1)k x
= i i2k+ ( )2k+
n>0
k! (1 + k + ) 2
X 1 x
= ( )2k+
k>0
k! (1 + k + ) 2

Remarque 2.2.2 Si est un entier on a cette relation

In (x) = I n (x)

Preuve. on montre aisement que

n
In (x) = (i) Jn (ix)
n n
= (i) ( 1) J n (ix)

n n n
= (i) ( 1) I n (x)(i)

= (i2 ) n
( 1) n
I n (x)

= I n (x)

Remarque 2.2.3 I (x) et I (x) sont deux solutions de l’équation de Bessel modi…ée.

Remarque 2.2.4 Lorsque est non entier,I (x) et I (x) ce sont deux solutions linéaire-
ments indépendantes, car J (ix) et J (ix) sont deux solutions linéairements indépen-
dantes.

Cherchons la 2eme solution de l’équation de Bessel modi…ée, linéairement indépen-

dante de In (x)

Théorème 2.2.2 quand est entier K (x) est une solution de l’équation de Bessel mod-
i…ée qui est linéairement indépendante de I (x)

40
Fonctions de Bessel

I (x) I (x)
Kn (x) = lim K = lim
!n 2 !n sin
Preuve. Si n’est pas un entier, K (x) s’écrit comme combinaison lineaire de
deux solutions lineairements indépendantes I (x) et I (x); on a

I (x) I (x)
K (x) = :
sin

et dans ce cas K (x)et I (x) sont deux solutions linéairements indépendantes de

l’équation de Bessel modi…ée. Par contre si est un entier, la fonction Kn (x) n’est
pas dé…nie et en conséquence , on dé…nit Kn (x) par

Kn = lim K
!n

et a…n de lever cette indétermination, on a :

@ @
@
I (x) @
I (x) =n
Kn (x) =
2 cos
n
( 1) @ @
= I (x) I (x)
2 @ @ =n

A ce stade, nous montrons que Kn (x) est une solution de l’équation de Bessel
Modi…ée et est linéairement indépendante de I (x).
1) Kn (x) est une solution de l’équation de Bessel Modi…ée,
lorsque I est un solution, on a

d2 I dI
x2 2
+x (x2 + 2
)I = 0
dx dx

de cette dernière équation, nous dérivons par rapport à , on trouve

d2 @I d @I @I
x2 2
+x (x2 + 2
) 2 I =0
dx @ dx @ @

de meme pour I ; on obtient

d2 d ( 1)
x2 2
K (x) + x K (x) + (x2 n2 )K (x) = (2 I (x) 2 I (x))
dx dx 2

41
Fonctions de Bessel

et à la limite lorsque ! n , le second membre est égal à 0; d’où Kn (x) satisfait

à l’équation de Bessel modi…ée
2) pour démontrer l’independance des solutions, calculons W (I (x); K (x)) ; on
a,

W (I (x); K (x)) = I (x)K 0 (x) I 0 (x)K (x) = W (I (x); I (x)) :

2 sin

Utilisons ([Link]) ; on trouve

1
W (I (x); K (x)) = 6= 0
x

et passons à la limite ! n, on déduit

1
W (In (x); Kn (x)) = 6= 0
x

ce qui prouve l’indepandance des deux solutions.

Remarque 2.2.5 K (x) est appelée fonction de Bessel modi…ée de deuxième espèce.

Remarque 2.2.6 On peut exprimer K (x) en fonction de J (x) et Y (x);

I (x) I (x)
K (x) = (2.2.34)
2 sin
i J (ix) i J (ix)
=
2 sin

et par de…nition Y (x) est dé…nie comme suit

Y (ix) sin = cos J (ix) J (ix)

Injectons l’expression J (ix) dans ([Link]), K (x) s’écrit

i cos J (ix) i Y (ix) sin i J (ix)

K (x) =
2 sin
2 1
+1 i cos i
= i iY (ix) + J (ix)
2 sin

et pour simpli…er la forme de K (x); on a

42
Fonctions de Bessel

2 1 2
i cos i = i cos + i:i
i
sachant que i = e 2 , alors
i
2 1 2
i cos i = i cos + i:(e 2 )
i
= i cos + i:e

= i cos + i(cos i sin )

= i sin

on obtient,
+1
K (x) = i [iY (ix) + J (ix)]
2
Remarque 2.2.7 la solution générale de l’équation peut être aussi écrite sous la forme
suivante
y = CI (x) + DK (x)

2.2.1 Relations de récurrence

Théorème 2.2.3 Quelques relations de récurrence pour les fonctions de Bessel modi…ées
de première espèce sont données comme suit :

d
1. dx
[x I (x)] = x I (x)

d
2. dx
[x I (x)] = x I1+ (x)

3. I 0 (x) = I 1 (x) x
I (x)

4. I 0 (x) = x I (x) + I +1 (x)

5. I 0 (x) = 12 [I 1 (x) +I +1 (x)]

2
6. I 1 (x) I +1 (x) = x
I (x)

Remarque 2.2.8 Ces relations de récurrence se démontrent de la même manière que

pour la fonction J (x), juste en changenant x par ix.

Remarque 2.2.9 Ces relations [T heoreme [Link] restent valables pour la fonction K (x);

43
Fonctions de Bessel

2.2.2 Forme intégrale des fonctions de Bessel modi…ées

les formes intégrales des fonctions de Bessel modi…ées sont données par les formules
suivantes :

i) pour la fonction I (x), on a

Z
( 12 x) 1
xt 1 1
I (x) = p e (1 t2 ) 2 dt, >
( + 1=2) 1 2

ii) pour la fonction K (x), on a

p Z 1
1 xt 1 1
K (x) = ( x) e (t2 1) 2 dt, > ,x>0
( + 12 ) 2 1 2

Preuve.

( 12 x) R1 1 1
i) On a I (x) = (i) J (ix) et J (x) = p
( +1=2) 1
(1 t2 ) 2 eixt dt,( > 2
); on a
alors Z
i ( 12 x) 1
1
I (x) = (i) p (1 t2 ) 2 e xt
dt
( + 1=2) 1

d’ou le résultat demandé.

ii) Avant de trouver la forme intégrale de K (x); on pose

Z 1
1 1
P (x) = x e xt (t2 1) 2 dt, > ,x>0
1 2

et nous véri…ons que P (x) satisfait l’equation de Bessel modi…ée, en premier lieu,
nous calculons
Z 1 Z 1
dP (x) 1 xt 2 1
xt 1
= x e (t 1) 2 dt + x te (t2 1) 2 dt (2.2.35)
dx 1 1

Z 1 Z 1
d2 P (x) 2 xt 2 1
1 xt 1
= ( 1)x e (t 1) 2 dt + 2 x te (t2 (2.2.36)
1) 2 dt
dx2 1 1
Z 1
1
+x t2 e xt
(t2 1) 2 dt
1

44
Fonctions de Bessel

et nous remplaçons ([Link]) et ([Link]) dans l’équation de Bessel modi…ée ([Link]),

on voit bien que P (x) satisfait l’équation de Bessel modi…ée
2
2d P (x) dP (x)
x 2
+ x (x2 + 2
)P (x)
Z dx 1
dx
1
= x ( 1) 2 xt + x2 t2 + xt x2 2 e xt
(t2 1) 2 dt
1
Z 1
+1 1 1
= x x(t2 1) 2( + )t e xt
(t2 1) 2 dt
1 2
Z 1
+1 xt + 12 1 1
= x e x(t2 1) 2te xt
( + )(t2 1) 2 dt
1 2
Z
+1
1
d h xt + 12
i
= x e (t2 1) dt
1 dt
h i1
+1 xt + 21
= x e (t2 1) =0
1

et donc la fonction P prend la forme suivante

P (x) = AI (x) + BK (x)

D’une part, on remarque que P est une fonction positive et tend vers 0 lorsque x ! 1.
xt
(t2 1) < e 2 pour x ! 1

donc
Z 1
xt
xt
P (x) < x e2e dt
Z1 1
xt
= x e 2 dt
1
1
2 xt
= x e 2
x 1
xt
1
= 2x e 2

et à la limite quand x tend vers 1; P vaut 0

Dautre part et d’après la dé…nition de la fonction de Bessel modi…ée de première espèce

I (x); on a I (x) ! 1 quand x ! 1. Alors l’expression de P s’écrira sous
la forme suivante :
P (x) = BK (x)

45
Fonctions de Bessel

ou bien
I (x) I (x)
=B
2 sin
et pour trouver la forme intégrale …nale de K (x); il reste à déterminer la constante B.

Pour cette raison, nous examinons le comportement de P (x) et K (x) lorsque x ! 0:

On a les formes I (x) et I (x) à l’ordre 1 en x sont respectivement :

1 x 1 x
( ) ; ( ) ;
(1 + ) 2 (1 ) 2

puis nous remplaçons le terme dominant dans K (x), on obtient

x
K (x) ( )
2 sin (1 ) 2
et utilisons cette propriète de la fonction Eulerienne ([Link])

on trouve
1
( )2
K (x) !
x
Pour étudier le comportement de P (x) près de x = 0;on fait le changement de
variable suivant
u
t=1+
x
P devient
Z 1
2u u2 1 1
P = x e x u( + 2 ) 2 du
0 x x x
Z 1
1 x 2x 1
= e e u (1 + ) 2 u2 1
du
x 0 u
et pour x ! 0; P se met
Z 1
1 u 2 1 1
P (x) ! e u du = (2 )
x 0 x
2x 1
x
où l’on utilisé e ' 1 et (1 + u
) 2 '1

Pour déterminer B; on compare les deux expressions suivantes

P (x) = BK (x)
1
1 ( )2
(2 ) = B
x x

46
Fonctions de Bessel

B sera d’apres la formule de duplication de Legendre [11]

22x 1
1
(2x) = p (x) (x + )
2

(2 ) ( + 12 )
B= = p
( )2 1

2.3 Autres types de fonctions de Bessel

2.3.1 Fonctions de Hankel

Les fonctions de Hankel parfois sont appelées fonctions de Bessel de troisième espèce,
sont dé…nies par :

H (1) (x) = J (x) + iY (x)

H (2) (x) = J (x) iY (x)

(1)
H (x) est appelle fonction de Hankel de première espèce

(2)
H (x) est appelle fonction de Hankel de la deuxième espèce

Remarque 2.3.1 En raison de l’indépendance linéaire de la fonction de Bessel de pre-

mière et deuxième espèce, les fonctions de Hankel fournissent une paire alternative de
solutions à l’équation di¤érentielle de Bessel

Remarque 2.3.2 On peut écrire K (x) sous la forme suivante:

+1
K (x) = i H (1) (ix)
2

2.3.2 Fonction de Kelvin

Soit l’équation di¤érentielle suivante

d2 y dy
x2 2
+x (i(kx)2 + 2
)y = 0 (2.3.37)
dx dx

47
Fonctions de Bessel

2
elle est écrite sous la forme de l’équation de Bessel modi…ée avec = ik 2 , alors sa
solution generale est

y = CI ( x) + DK ( x)
1 1
= CI (i 2 kx) + DK (i 2 kx)

comme I (x) = (i) J (ix), on déduit

3 1
y = C(i) J (i 2 kx) + DK (i 2 kx)

et nous obtenons des fonctions réelles données par les dé…nitions suivantes

3
ber (x) = Re J (i 2 x)
3
bei (x) = Im J (i 2 x)
1
ker (x) = Re(i) K (i 2 x)
1
kei (x) = Im(i) K (i 2 x)

Remarque 2.3.3 ber (x) bei (x) ker (x) et kei (x) sont appellées fonctions de Kelvin.

2.4 Fonctions de Bessel sphériques

Dans cette section, on s’intéresse à étudier l’équation di¤érentielle suivante :

d2 y dy
x2 + (1 2 )x + (( x )2 + ( 2 2 2
)y = 0 (2.4.38)
dx2 dx

elle est réductible à l’équation de Bessel où et sont des paramètres

complexes.
Preuve. En premier lieu, e¤ectuant le 1er changement suivant y = x z alors
8
< dy
=x dz
+ x 1
z
dx dx
: d2 y
=x d2 z
+2 x 1 dz
+ ( 1)x 2
z
dx2 dx2 dx

Alors

48
Fonctions de Bessel

2
+2 dz dz dz
x 2
+2 x +1 + ( 1)x z+(1 2 ) x +1 + x z + 2 2 x2 + ( 2 2 2
) x z=0
dx dx dx
et par un calcul direct, on trouve
d2 z dz
x2 2 + x + 2 2 x2 2 2
z = 0:
dx dx
en deuxième lieu, en e¤ectuant un autre changement de variable t = x , il vient
8
< dz
= dz dt
= x 1 dz
dx dt dx dt
: d2 z
= ( 1)x 2 dz
+ 2 2
x 2 d2 z
dx2 dt dt2
Aprés le remplacement, on obtient

2
2 2d
z dz dz
t
2
+ ( 1)t + t + 2 2 t 2 2
)z = 0
dt dt dt
en simpli…ant, on trouve
d2 z dz
t2 2 + t + 2 t 2
)z = 0:
dt dt
On voit clairement que cette équation n’est autre que l’équation de Bessel modi…ée
dont sa solution est donnée par la forme suivante

z = AJ ( t) + BY ( t)

y = x (AJ ( x ) + BY ( x ))

Soit l’équation suivante

d2 y dy
x2 2 + 2x + (( x)2 l(l + 1))y = 0
dx dx
on remarque, qu’elle est de la forme de l’équation precédente ([Link]), avec:
8 8
>
> 1 2 = 2 >
> = 21
>
> >
>
>
< >
<
=1 =
)
>
> 2 2
= 2 >
> =1
>
> >
>
>
: 2 >
:
2 2
= l(l + 1) = l + 12
sa solution s’écrit directement sous la forme suivante

1
y=x 2 AJl+ 1 ( x) + BYl+ 1 ( x)
2 2

49
Fonctions de Bessel

2.5 Comportement des fonctions de Bessel dans les

cas limites
Dans cette section, on étudie les comportements des fonctions de Bessel

Théorème 2.5.1 Lorsque x ! 0

1 x
1. J (x) ( +1) 2
8
>
> 1
( ) x
( > 0)
>
< 2
2. Y (x) cos 1 x
( < 0)
>
>
sin ( +1) 2
>
: 2
ln x + ( = 0)

Preuve. 1) J (x) est une série entière ,il su¢ t de prendre le premier terme
2) Substituons l’expression équivalent de de J (x) et J (x) dans Y (x) on obtient

1 2 cos 1 x
Y (x) ( ) +
x sin ( + 1) 2

pour > 0; en prenant le premier terme dominant (car x plus petit), on obtient
1 2
Y (x) ( ) x

pour < 0; le deuxième terme est dominé (x plus petit), on obtient Y (x)
cos 1 x
sin ( +1) 2

quand = 0, en utilisant la formule explicite de Y (x):

Théorème 2.5.2 Lorsque x ! 1

1
2 1
1 J (x) x
2
cos x + 2 2

1
2 1
2 Y (x) x
2
sin x + 2 2

Preuve.
voir[9; ChapVII]

50
 Partie II

Polynômes de Bessel

51
Polynômes de Bessel

2.6 Polynômes de Bessel

2.6.1 Résolusion de l’équation de polynômes de Bessel

Les polynômes de Bessel sont les solutions polynomiales de l’équation di¤érentielle
suvante :

d2 y dy
x2 2
+ 2(x + 1) n(n + 1)y = 0; n 2 N+ (2.6.39)
dx dx
On résoud cette équation di¤érentielle à l’aide des séries entières

X
y(x) = ak xk+ , a0 6= 0 (2.6.40)
k>0

Remplaçons ([Link]) dans ([Link]) et multiplions le résultat obtenu par x, il vient :

X X X X
(k+ )(k+ 1)ak xk+ +1
+2 (k+ )ak xk+ +1
+2 (k+ )ak xk+ n(n+1) ak xk+ +1
=0
k>0 k>0 k>0 k>0

En égalant à zéro les coe¢ cients de même degré, alors pour k = 0, on a :

2a0 = 0 ) =0

et y(x) se simpli…e comme :

X
y(x) = ak x k (2.6.41)
k>0

Maintenant, injectons ([Link]) dans ([Link]) ; on obtient cette relation de recurrence

entre les coe¢ cients de la série

2(k + 1)ak+1 + ak [k (k + 1) n(n + 1)] = 0

ce qui implique
[k (k + 1) n(n + 1)]
ak+1 = ak
2(k + 1)
Constatons que pour k = n; on a : ak+1 = ak+2 = ak+3 = ::: = 0:

52
Polynômes de Bessel

Pour déterminer le terme générale de la série ([Link]) ; calculons les coe¢ cients
de la série, on a alors les relations suivantes :

n(n + 1) n(n + 1)(n 1)(n + 2)

a1 = a0 ; a2 = a0 :::::::::
2 8

le terme général sera alors

(n + k)
ak+1 = a0
k!(n k)!2k

D’où la solution de la série yn (x) s’écrit :

x n(n + 1)(n 1)(n + 2) x 2 (2n)! x n

yn (x) = a0 1 + n(n + 1) + + ::: +
2 2 2 n! 2

où sous forme condensée

Xn
(n + k)! x k
= ; a0 = 1
k=0
k!(n k)! 2

A ce niveau, pour trouver la deuxième solution linéairement indépendante de l’équation

2
di¤erentiellen ([Link]), en utilisant le changement suivante y(x) = (x)e x et aprés un
calcul simple, on obtient la solution suivante

2
y(x) = e x yn ( x)

et la solution générale est donnée par

Y = Ayn (x) + By(x)

2.6.2 Relation entre les polynômes et les fonctions de Bessel

Soit l’équation di¤érentielle linéaire suivante :

d2 f df
x2 + 2x + x2 n(n + 1) f = 0; n 2 N+ (2.6.42)
dx2 dx
1
On fait le changement suivant f (x) = x 2 J(x) et par un calcul simple et direct
on trouve ,
2 2
2d J dJ 1
x 2
+x + (x2 n+ )J = 0;
dx dx 2

53
Polynômes de Bessel

qui est une équation di¤érentielle de Bessel d’ordre n + 21 dont la solution s’exprime en
terme de Jn+ 1 et J n 1 :
2 2
w(x)
Pour la même équation, en posant f (x) = x
, l’équation se convertira à

x2 (w" + w) n(n + 1) = 0

ix
puis par un autre changement de variable z = ix et w(x) = e y(ix), l’équation en
haut prendra la forme suivante
d2 y dy
z2 2 n(n + 1)y = 0
dx2 dx
et en…n par un autre changement de variable t = z1 ; il est aisé d’avoir cette équation
di¤érentielle ,
d2 y dy
2
t2
+ 2(t + 1) n(n + 1)y = 0;
dt dt
qui est l’équation di¤érentielle associée aux polynômes de Bessel.
Trouvons à ce niveau les relations reliantes les polynômes et les fonctions de Bessel
d’ordre n + 12 [12; Chap5]

1 x
1
2
h i
yn ( ) = eix i n 1
Jn+ 1 (x) + in J n 1 (x) (2.6.43)
ix 2 2 2

et on déduit
8 1
< J (x) = x 2
eix i n 1
yn ( 1
) +e ix n+1
i yn ( ix1 )
n+ 21 2 ix
1 (2.6.44)
: J 1 (x) = x 2 ix n
e i yn ( 1
) +e ix
i n
yn ( ix1 ) :
n 2 2 ix

2.6.3 Orthogonalité
Il est bien connu que les polynômes de Bessel satisfont la propriété “d’orthogonalité”
2
sur un cercle d’unité C tel que la fonction poids dans le plan complexe vaut (x) = e x

.
Nous montrons que les polynômes de Bessel forment un système orthogonal :

Théorème 2.6.1
8
Z < 0
1 m 6= n
(z)ym (z)yn (z)dz =
2 i pzp=1 : ( 1)n+1 2 m=n
2n+1

où ym (z) sont des polynômes de Bessel.

54
Polynômes de Bessel

Remarque 2.6.1 La fonction poids et le cercle d’unité ne sont pas unique, car la fonction
poid multipliée par un constante ou avec une fonction analytique additive sur une courbe
arbitraire autour de l’origine véri…e toujours la relation d’orthogonalité.

Preuve. Si m 6= n
l’équation di¤érentielle peut s’écrire sous la forme :

2 2
(z 2 e z yn0 )0 n(n + 1)e z yn = 0

En multipliant les deux membres par ym et en intégrant par partie autour du cercle
d’unité on trouve :

Z Z
2 2
n(n + 1) e z yn ym dz = yn0 )0 ym dz
(z 2 e z (2.6.45)
pzp=1 pzp=1
Z
2 2
2 0
= pzp=1 z e z yn ym z2e z yn0 ym
0
dz
pzp=1

En échangeant n par m, alors :

Z Z
2 2
m(m + 1) e z ym yn dz = (z 2 e z
0 0
ym ) yn dz
pzp=1 pzp=1
Z
2 2
2 0
= pzp=1 z e z ym yn z2e z
0 0
ym yn dz (2.6.46)
pzp=1

Faisons la soustraction de ([Link]) et ([Link]) , on obtient(1)

Z
2
[(n m)(n + m + 1)] e z ym yn dz = 0;
pzp=1

comme m 6= n Z
1 2
) e z ym (z)yn (z)dz = 0
2 i pzp=1

si m = n voir [12; Chap4]

(1)
I les 1er terme dans ([Link]) et dans ([Link]) sont nuls selon la formule de Green- Riemann
d(f ) = 0

55
Polynômes de Bessel

2.6.4 Relations de récurrence

Il est utile aussi d’exposer quelques fameuses relations de récurrence pour les polynômes
de Bessel, on a :

Théorème 2.6.2 Soit les relations de récurrence pour les polynômes de Bessel suivantes
:

1. yn+1 = (2n + 1)xyn + yn 1

2. x2 yn0 = (nx 1)yn + yn 1

3. x2 yn0 1 = yn (nx 1)yn 1

4. x(yn0 + yn0 1 ) = n(yn yn 1 )

5. (nx + 1)yn0 + yn0 1 = n2 y n

Avec y 1 (x) = y0 (x) = 1

[12; Chap 3]

2.6.5 Formule de Rodrigues

L’élément yn (x) de la suite des polynômes de Bessel fyn (x)g est déterminé
par la formule de Rodrigues suivante :

n 2 dn 2
yn (x) = 2 ex (e x x2n )
dxn
voir [12] [13] [14]

2.6.6 Fonction génératrice

La fonction génératrice G(x; t) pour les polynômes de Bessel est donnée par l’expression
suivante [13]: !
1
X tn 1 (1 2xt) 2
x
G(x; t) = yn 1 (x) =e
n>0
n!

Remarque 2.6.2 y 1 (x) = y0 (x) = 1

56
Polynômes de Bessel

Preuve. La fonction génératrice peut être obtenue à partir des formules suivantes
[9]:

X tn 2 1 1
Jn 1 (r) = ( ) 2 cos(r2 2rt) 2 (2.6.47)
n>0
2 n! r
X tn 2 1 1
J n+ 12 (r) = ( ) 2 sin(r2 2rt) 2 (2.6.48)
n>0
n! r

D’apres ([Link]) on sait que

1 1 1
h i
yn 1 ( ) = ( ) 2 r 2 eir i n
Jn 1 (r) + in 1 J n+ 21 (r) (2.6.49)
ir 2 2

tn
Nous multiplions ([Link]) par n!
et nous utilisons ([Link]) et ([Link]) , nous obtenons
X 1 tn h 1 1
i
yn 1 ( ) = eir cos(r2 2irt) 2 i sin(r2 + 2rt) 2
n>0
ir n!

ou bien
1
(r2 +2irt) 2 )
= ei(r
1
posons x = ir
, cette dernière relation devient
1
X tn i 1
ix ( i21x2 + 2tx ) 2
yn 1 (x) = e
n>0
n!
1

= e
1
x [ 1 x2xt
2 ]
2

ou
X 1 tn 1
1
(1 2xt) 2
yn 1 ( ) = e x

n>0
ir n!

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

Les polynômes de Bessel généralisés yn (x; a; b) sont les solutions de l’équation di¤éren-
tielle suivante :

d2 y dy
x2 2
+ (ax + b) n(n + a 1)y = 0 ; n 2 N+
dx dx

57
Polynômes de Bessel

Remarque 2.7.1 Dans le cas a = b = 2; on retombe sur l’équation des polynômes

de Bessel ([Link]).

Les polynômes de Bessel généralisés peuvent se representer par une série hyper-
géométrique suivante :
yn (x; a; b) =
X
n
n!(n + k + a 2)(k) x k x
= :2 F0 ( n; a 1 + n; :; ):
k=0
k!(n k)! b b
qui s’écrit aussi comme :
X
n
x k
yn (x; a; b) = Cn;k (n + k + a 2)(k)
k=0
b

où z (k) = z(z 1):::(z k + 1) et Cn;k sont les coe¢ cients binomials.

2.7.1 Formule de Rodrigues pour les polynômes de Bessel

généralisés
La formule de Rodrigues pour les polynômes de Bessel généralisés est donnée par :

n 2 a b dn 2n+a 2 b
yn (x; a; b) = b x ex (x e x )
dxn

2.7.2 Relations de récurrence pour les polynômes de Bessel

généralisés
la Relation de récurrence pour les polynômes de Bessel généralisés est donnée par

h x i
(n+a 1)(2n+a 2)yn+1 = (2n + a)(2n + a 2) +a 2 (2n+a 1)yn +n(2n+a)yn 1
b

Remarque 2.7.2 les deux premières relations de récurrence pour les polynômes de
Bessel [Théorème 2.7.2] peuvent s’obtenir en posant a = b = 2.

58
Polynômes de Bessel

2.7.3 Orthogonalité des polynômes de Bessel généralisés

Soit (x)la fonction poids pour les polynômes de Bessel généralisés, tel que :

1 X
n
(a) b
(x) =
2 i n>0 (n + a 1) x

les polynômes de Bessel généralisés véri…ent la relation d’orthogonalité suivante.

8 8
< 1 Z < 0 m 6= n
ym (z; a; b)yn (z; a; b) (z; a; b)dz =
: 2 i pzp=1 : ( 1)n+1 bn! (a)
m=n
2n+a 1 (n+a 1)

Pour plus de détails concernant ces propriétés pour les polynômes de Bessel généralisés
voir [12] :

59
Chapitre 3

APPLICATIONS

3.1 Equation de la chaleur

Comme il est connu que l’équation de la chaleur est donnée par l’équation
di¤érentielle suivante :

@ 2u z @ @u
2
= ,06 61,t>0,z2R (3.1.1)
@t @ @

ou bien

@ 2u 1 @u @ 2 u
=z + 2 (3.1.2)
@t2 @ @
A…n de solutionner cette équation ([Link]), nous utilisons la méthode de séparation
des variables, on pose

u( ; t) = W ( )G(t)

Par un simple calcul, il vient

G" 1 W0 W"
= + = ; >0
zG W W

On obtient alors un systéme d’équation di¤érentielle découplé suivant

8
< G" + zG = 0
(3.1.3)
: W " + 1W 0 + W = 0

60
APPLICATIONS
p
Pour simpli…er la 2eme équation, on pose s = , il vient
dW dW ds p dW
= = (3.1.4)
d ds d ds
2 2
dW dW
=
d 2 ds2
En remplaçant ([Link]) dans ([Link]), on obtient
d2 W dW
2
+ + W =0 (3.1.5)
ds s ds
s2
puis en multipliant cette dernière équation ([Link])par ; il vient

d2 W dW
s2 2
+s + s2 W = 0 (3.1.6)
ds ds
Qui est exactement l’équation de Bessel d’order 0 donc sa solution générale se met

W (s) = CJ0 (s) + DY0 (s)

pour la première équation ([Link]), sa soulution est donnée par

p p
G(t) = a cos( z t) + b sin( z t):

Alors l’expression …nale de la solution de l’équation de la chaleur s’écrira alors

X p p
u( ; t) = Gn (t)(CJ0 ( ) + DY0 ( ))
n

3.2 Equation d’Helmholtz

Une deuxième application est de grande importance dans le domaine de la physique,
décrivant la propagation des ondes de D’Alembert, et est connue par l’équation
d’Helmholtz (d’après le physicien Hermann von Helmholtz) suivante :

4 + k2 =0 (3.2.7)

avec une solution continue à l’intérieur d’une couronne sphérique de rayons a; b(a <
b): L’équation ([Link]) exprimée en coordonnées sphériques s’écrit :
1 @ 2@ 1
r + + k2 (r; ; ') = 0 (3.2.8)
r2 @r @r r2

61
APPLICATIONS

où est la partie angulaire du Laplacien 4 :

1 @ 1 @2
= + (3.2.9)
sin @ sin2 @'2

Par la méthode de séparation des variables de la forme (r; ; ') = R(r)F ( ; '), on
aura le système découplé suivant

1 @ 2@ F
r R + k2 r2 = =C (3.2.10)
R @r @r F

ou
8
< @ 2 @
r @r + k2 r2 C R(r) = 0
@r
(3.2.11)
: [ + C] F ( ; ') = 0
La deuxième équation di¤érentielle est connue en mécanique quantique , décrit la
dynamique angulaire du système et admet des solutions sur la sphère avec C = l (l + 1)
où l est un nombre entier quantique positif jouant le rôle de la valeur propre du
moment angulaire , sa solution est donnée par
1
(2l + 1)(l m)! 2
m
Yl ( ; ') = ( 1) m
eim' Plm (cos ); 0 6 l 6 1 , l6m6l
4 (l + m)

avec Plm sont les polynômes de Legendre[1].

Revenons à la la première équation di¤érentielle radiale du sysème ([Link]) qui
est en relation avec notre exposé, en posant x = kr; l’équation se transforme alors

d2 R 2 dR l(l + 1)
+ + 1 R = 0; (3.2.12)
dx2 x dx x2

ou bien
d2 R dR
x2 2
+ 2x + x2 l(l + 1) R = 0: (3.2.13)
dx dx
On voit que cette équation représente exactement l’équation de Bessel sphérique ([Link]),
avec les paramètres suivants
1
= 2

= =1
2 2 2
=
1
=l+ 2

62
APPLICATIONS

La solution de cette dernière équation s’écrira sous la forme suivante :

1
y = x 2 AJl+ 1 (x) + BYl+ 1 (x)
2 2
1
= (kr) 2 AJl+ 1 (kr) + BYl+ 1 (kr)
2 2

et la solution …nale cherchée du problème posé est :

1
X X
l
(r; ; ') = (kr) 2 Alm Jl+ 1 (kr) + Blm Yl+ 1 (kr) Ylm ( ; ')
2 2
l>0 m= l

3.3 Equation du pendule simple

Dans cette application, on traitera un probleme de la mecanique classique et il s’agit
des oscillations d’un pendule simple. En général, la dynamique d’un système classique
est déterminée par une fonction dite fonction Lagrangienne, qui est la di¤érence des
deux énergies cinétique et potentielle, donnée par:
1
L( ; _ ; t) = f ( ; _ ; t) = ml2 _ + mgl cos (3.3.14)
2
où _ est la vitesse angulaire donnée par _ = ddt , m est la masse du point matériel
attaché au pendule, l est la longueur du pendule dépend du temps vaut l = t et g
est la pesenteur.

Maintenant pour déterminer les équations du mouvement du système, nous

appliquons l’équation de Euler -Lagrangele déduite du principe de moindre action
, Z
f ( ; _ ; t)dt = 0

il vient les équations de mouvemet suivante

@f d @f
=0 (3.3.15)
@ dt @ _
@f @f d @f
Nous calculons @
; @_ et dt @ _
;
8
>
> @f
=mgl sin
>
< @
@f
= ml2 _ = m 2 t2 _ (3.3.16)
>
> @_
>
: d2
d @f
dt @ _
= 2m 2
t_ + m 2 2•
t avec • = dt2

63
APPLICATIONS

puis les injectons dans ([Link]), on obtient cette équation di¤érentielle non lineaire.

2 2•
m t 2m 2
t_ mgl sin = 0 (3.3.17)

Dans notre étude, on se limite seulement aux petites oscillations; en prenat sin ;
il vient

2 2•
t +2 2
t _ + gl = 0 (3.3.18)

Par un changement adéquat suivant,

8 q
< u = 2 gt
(3.3.19)
: (u) = u (u)

et par un simple et direct calcul, les dérivées …gurantes dans l’équation di¤érentielle
([Link]) s’écrivent en fonction des nouvelles variables comme suit ,
8 pg
>
> _ t = d = d du ; du =
>
< dt du dt dt t
_ t = 2g _ u (3.3.20)
>
> h i u
>
: •t = 2
•u _ d2
2g
u
1
u u
avec _ u = du
d
et •u = du 2

Notre équation ([Link]) se sipmli…er à

d2 d
u2 2
+u + (u2 1) = 0 (3.3.21)
du du

Cette dernière équation admet des fonctions de Bessel d’order 1 comme solution

(u) = AJ1 (u) + BY1 (u) (3.3.22)

et la solution en fonction de l’ancienne variable devient

r r
C gt D gt
(t) = p J1 (2 ) + p Y1 (2 )
t t

3.4 Equation d’onde en coordonnées sphériques

On consacre cette séction à une appliquation pour les polynome de Bessel, de ce fait,
on traitera l’équation d’onde suivante:

64
APPLICATIONS

1 @ 2u
4u = 2 2 (3.4.23)
c @t
où 4 est le laplacien et c est la vitesse de propagation des ondes.
Utilisons les coordonnées sphérique et la méthode de séparation des variables, on
pose :
u(r; ; ; t) = R(r) ( ) ( )T (t) (3.4.24)

Ecrivons l’operateur Laplacien en coordonnées sphérique puis insérons ([Link]) dans

l’équation ([Link]) ; on trouve les expressions des parties sépariées espace et temps
suivantes :

R" 2 R0 1 0 1 1 "
1 "
1 T"
+ + 2 + 2 + = = k2 (3.4.25)
R rR r tan r r sin2
2 c2 T
8
< R" + 2 R0 + 1 0 1
+ 1 "
+ 1 "
= k2
R r R r2 tan r2 r2 sin2
(3.4.26)
: T " + c2 T k 2 = 0
Nous remarquons que la première expression est une autre équation aux dérivées
@
partielles ; @
@r @
et @
@
; utilisons de nouveau la méthode de séparation des variable , il
vient

R" 2 R0 1 0
1 1 " "
r2 sin2 + + 2 + 2 + k2 = = q2: (3.4.27)
R rR r tan r

ou bien le système suivant

8
< "
+ q2 =0
(3.4.28)
: r2 sin2 R"
+ 2 R0
+ 1 0 1
+ 1 "
+ k2 = q2
R r R r2 tan r2

Il est clair que la solution de la première équations di¤érentielle angulaire en ( ) est

;
( ) = aq cos q + bq cos q : (3.4.29)

et pour résoudre la deuxième équation di¤érentielle du système, on a besoin d’une

autre constante de séparation,

2 R"R0 0
1 "
q2
r + 2r + r2 k 2 = + + = n(n + 1) (3.4.30)
R R tan sin2

65
APPLICATIONS

ou sous la forme suivante

8
< r2 R" + 2r R0 + r2 k 2 = n(n + 1)
h 0R R i (3.4.31)
: 1
+
"
+ q2
= n(n + 1)
tan sin2

1
En simpli…ant la première équation ([Link]), en posant, R(r) = ixe x y(x) avec x = ikr
;
l’équation radiale ([Link])devient :
2
2d y dy
x + 2(x + 1) n(n + 1)y = 0
dx2 dx

qui est l’équation des polynômes de Bessel.

66
CONCLUSION

Dans ce mémoire, nous avons principalement exposé un modeste travail sur les fonctions
et les polynômes de Bessel , qui sont considérés comme un cas particulier de la grande
famille des fonctions spéciales et que leurs applications n’ont cessé de se développer et
d’envahir plusieurs domaines et disciplines scienti…ques, en particulier la physique.
Ce mémoire est composé principalement de trois chapitres. Dans le premier chapitre,
nous avons présenté quelques dé…nitions générales, quelques propriétés et notions es-
sentielles que nous les avons utilisées, citons par exemple: les fonctions Gamma, Béta
et hypergéométrique , les polynômes orthogonaux et la méthode de Frobenius. Le
deuxième chapitre est consacré à l’étude des fonctions et polynômes de Bessel.
Pour les fonctions de Bessel, plusieurs dé…nitions et propriétés essentielles ont
été présentées : comme les relations de récurrence, représentation intégrale, fonction
génératrice, les fonctions de Bessel de première et de deuxième espèce, fonctions de
Bessel modi…ées et sphériques
Pour les polynômes de Bessel, quelques propriétés importantes ont , leurs relation
avec les fonctions de Bessel .
Dans le troisième chapitre, nous avons présenté quelques applications sur la dy-
namique de certains systèmes physiques gouvernée par des fonctions et polynômes de
Bessel, nous citons: le pendule simple , l’équation de la chaleur, l’équation d’Helmholtz
et l’équation d’onde en représentation sphérique.
En …n pour terminer cette conclusion, on mentionne que le traitement d’autres
fonctions spéciales comme, celles de Heun et de Romanoveski ainsi l’extension de ces
fonctions de Bessel dans le cadre de la théorie de la déformation, est d’une grande
importance et nécessite vraiment une étude mathématiques profonde , ce que nous

67
Conclusion

comptons réaliser.

68
Bibliographie

[1] Mécanique quantique, tome3 : Théorie non relativiste, L D Landau et E. Lifchitz

, edition mir 1967; Mécanique quantique Tome 1, C. Cohen-Tannoudji, [Link] et
[Link]ë, Editions Hermann (1997).

[2] C B Compean and M Kirchbach, J. Phys. A: Math. Gen. 39 (2006) 547.

[3] G. J. Papadopoulos and J. T. Devreese, Phys. Rev. D 13, 2227.

[4] [Link] and [Link], Reports on Mathematical Physics, 77 (2016) 69.

[5] [Link], [Link], [Link], Fractional integrals and Derivatives theory

and Applications. Gordon and Breach Science Publishers, Amsterdam (1993).

[6] M.B. Zahaf and D. Manchon, arXiv:1206.2522 [[Link]].

[7] Special functions for scientists and engineers, W. W. Bell, [Link] nostrand company
LTD, 1967.

[8] Outils mathématiques à l’usage des scienti…ques et ingénieurs, Elie Belorizky :,EDP
sciences, 2007.

[9] A Treatise on the Theory of Bessel Functions, G. N. Watson ,Cambridge University

Press 1922.

[10] Fonctions spéciales de la physique mathématique, [Link] , O¢ ce des

publications universitaires, 1983.

[11] Table of Integrals, Series and Products ,I. S. Gradshteyn and I. M. Ryzhik, Academic
Press, New York, 1979.

69
Bibliographie

[12] Bessel polynomials, Emill Grosswald,Springer-Verlag Berlin Heidelberg New York

,1978.

[13] Emill Grosswald, Les groupesde Galois des polynômes de Bessel,Séminaire

[Link]èbre et théorie des nombres, vol.18, N 17, pp 2, 15 mars 1965.

[14] [Link] et [Link] : A New class of orthogonal polynomials : the Bessel polyno-
mials, [Link], vol.65, pp 100-115, 1949.

[15] Analyse De Fourier et application aux problemes de valeurs aux limites, M.R. Spiegel,
France Mac Graw Hill (1999).

[16] [Link], Polynômes orthogonaux, Québec- Canada

[17] Mathématiques de l’ingénieur, Rémi Vaillancourt, Département de mathématiques,

Université d’Ottawa, Ottawa, ON, Canada, K1N 6N5

70