RÉPUBLIQUE DÉMOCRATIQUE DU CONGO
MINISTÈRE DE L’ENSEIGNEMENT SUPÉRIEUR ET UNIVERSITAIRE
INSTITUT SUPÉRIEUR ET PÉDAGOGIQUE DE BOMA
SECTION : SCIENCES & TECHNOLOGIE
DÉPARTEMENT : P.T.A.E.
OPTION : PHYSIQUE & ÉLECTRICITÉ / PHYSIQUE-TECHNOLOGIE
PROMOTION : PREMIÈRE LICENCE
TP EXAMEN
DE
MATHÉMATIQUES APPLIQUÉES
Dirigé par : le CT. NGOMA MBUKU Maurice Texer.
ÉTUDE & RECHERCHES SUR LES
ÉQUATIONS DE BESSEL.
Présenté par : -
• UMBA di KIPHUNDI Christophe
• MALONDA LUVUNU Joseph
• NGOMA VANGU Léopold
• MUANDA DIFUENI François
• KHONDE KHONDE robert
• VANGU KAMBA Jean Claude
• MVUELA DISEKELE Polydore
• MALELE MANGONGO Gérard
OCTOBRE 2023
copyright@2023 licence de physique /ISP-BOMA
1
Table des matières
I. INTRODUCTION ........................................................................................................................................... 2
II. RÉSOLUTION DE L’ÉQUATION DE BESSEL. ......................................................................................... 3
C’EST L’ÉQUATION DE BESSEL MODIFIÉE .............................................................................................................. 3
III. RESOLUTION PAR LA MÉTHODE DES SÉRIES. .................................................................................. 4
IV. L’ÉQUATION INDICIELLE DE L’EDO DE BESSEL............................................................................... 4
V. APPLICATION DES THÉORÈMES DE FROBENIUS DANS LA RÉSOLUTION : .................................. 5
VI. SOLUTION DE L’ÉQUATION DIFFÉRENTIELLE DE BESSEL ............................................................ 8
VII. L’ÉQUATION DE BESSEL DE 2EME ESPECE ......................................................................................... 9
VIII. APPLICATIONS A LA RESOLUTION D 'UN PROBLEME DE ............................................................ 9
PHYSIQUE MATHEMATIQUE ....................................................................................................................... 9
IX. CONCLUSION ........................................................................................................................................... 17
X. RÉFÉRENCES............................................................................................................................................. 18
2
I. INTRODUCTION
Les fonctions de Bessel apparaissent NATURELLEMENET comme fonctions propres
d’opérateurs différentiels liés à l’études de certaines équations de la physique Mathématique. Elles
font partie des principales familles des fonctions spéciales de la physique.
Cependant, leur développement en série intervient aussi de façon importante en combinatoire
énumérative.
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.En1816,
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.
Dans cette étude, nous comptons montrer comment arriver aux solutions de l’équation
différentielle de Bessel en aboutissant à ses équations différentielles de premier ordre comme celle
de deuxième ordre pour ainsi donner la voix à ceux qui voudraient les utiliser pour savoir d’où
proviennent ces formules d’équations qui semblent très complexe de prime abord.
Nous nous y sommes attelés en neuf points dont un sur une application concrète de ces équations
sur son utilité dans le quotidien d’un physicien.
Loin de nous la prétention d’avoir traité tous les aspects de ce problème ; pourtant nous avons
essayés, autant que faire ce peut, d’épuiser la petite partie qui nous concernait à ce sujet suivant
les recommandations de notre enseignant soit : montrer la résolution en série des équations
différentielles de Bessel y compris à l’aide de la méthode de Frobenius.
3
II. RÉSOLUTION DE L’ÉQUATION DE BESSEL.
Soit une équation différentielle de second ordre à coefficients variables suivante :
𝑥 " 𝑦" + 𝑥𝑦 ' + (𝑥 " − 𝑎" )𝑦 = 0. ( 𝐈𝐈. 𝟏)
Cette équation est dite « Forme standard de l’équation différentielle de Bessel » qui est une
équation différentielle linéaire du deuxième ordre; son intérêt en ingénierie pour les problèmes de
chaleur, de communication et de phénomènes vibratoires est d’importance capitale.
2
Posons : x = iz. Et. i2 = - 1 donc 𝑧 = = - ix
3
En appliquant la règle de dérivation des fonctions composées à y' de l’équation (II.1) on a :
45 67 69 𝐝𝐲 69
𝑦 ' = 42 = 48 42
= −𝒊 𝒅𝒛
, car 42
= −𝑖
et
𝑑" 𝑦 d 𝑑𝑦 dz d 𝒅𝒚
𝑦 '' = "
= F G = F −𝒊 G (−𝑖)
𝑑𝑥 𝑑𝑧 𝑑𝑥 𝑑𝑥 𝑑𝑧 𝒅𝒛
𝑑" 𝑦 d 𝑑𝑦 dz d 𝑑" 𝑦
𝑦 '' = = F G = (−𝑖)(−𝑖) = −1
𝑑 𝑥" 𝑑𝑧 𝑑𝑥 𝑑𝑥 𝑑𝑧 𝑑 𝑧"
Or x2= i2z2. donc x2= -z2
L’équation standard de Bessel: 𝑥 " 𝑦" + 𝑥𝑦 ' + (𝑥 " − 𝑎" )𝑦 = 0. Devient alors :
−𝑧 " [−𝑦 '' ] + 𝑖𝑧[−𝑦 ' ] + ( −𝑧 " − 𝑎" )𝑦 = 0. (𝐈𝐈. 𝟐)
Où dans. 𝑥 " 𝑦"; y = y(z) et dans (𝑥 " − 𝑎" )𝑦 ; y = y(x)
è z2 . y" + z. y' – (x2 – a2)y = 0
II.2 peut s’écrire alors :
𝒛𝟐 𝒚'' + 𝒛𝒚' − ( 𝒙𝟐 + 𝒂𝟐 )𝒚 = 𝟎
C’est l’équation de Bessel modifiée
4
III. RESOLUTION PAR LA MÉTHODE DES SÉRIES.
Soit une équation différentielle de second ordre à coefficients variables suivante :
𝑥 " 𝑦" + 𝑥𝑦 ' + (𝑥 " − 𝑎" )𝑦 = 0. ( 𝐈𝐈. 𝟏)
Écrivons la sous la forme suivante :
𝑦" + 𝑝(𝑥)𝑦 ' + 𝑞(𝑥)𝑦 = 0
2 U 2 T V WT
Avec p(x) = = et q(x) = = on a :
2T 2 2T
1 ' 𝑥2 − 𝑎2
𝑦 '' + 𝑦 +X Z𝑦 = 0
𝑥 𝑥2
U
p(x) =
2
𝑥2 − 𝑎2
q(x) = F 𝑥2
G
On constate que les fonctions p et q sont singulières en x0 = 0 ; par conséquent nous allons
développer l’expression en cherchant sa série de Taylor :
U
(x-x0).p(x) = x = 1 + 0𝑥 + 0𝑥2 + ⋯. Avec une Cste P0= 1
2
U 2 T V WT
(x-x0)2.q(x) = x = 𝑥2 \ ] = −𝑎2 + 0𝑥 + 1𝑥2 + ⋯ Avec une Cste q0= −𝑎"
2 2T
IV. L’ÉQUATION INDICIELLE DE L’EDO DE BESSEL.
Qui a La forme : 𝑟 " + (p0 − 1)𝑟 + q0 = 0
L’équation indicielle de Bessel peut s’écrire en vertu des valeurs de P0= 1 et et q0= −𝑎"
𝑟 " + (1 − 1)𝑟 + (−𝑎" ) = 0
𝑟 " − 𝑎" = 0
Équation qui admet 2 racines ou 2 solutions qu’on va noter :
5
r1 = a ³ 0 positif
r2 = - a £ 0 négatif
quelle que soit la nature de la différence entre r1 et r2 ; cette équation de Bessel
admet une solution série qui s’écrit :
h
y = b Cn . 𝑥 efgU
eij
V. APPLICATION DES THÉORÈMES DE FROBENIUS DANS LA
RÉSOLUTION :
Appliquons le théorème de Frobenius : comme r1 vaut a , nous cherchons la solution en série de
l’équation de Bessel qui a la forme :
h
y = b Cn . 𝑥 efW
eij
Il faut donc connaître les dérivées du premier ordre et du deuxième ordre de y , ce qui donne :
Et si on prend la dérivée de y par rapport à x on a :
y′ = b(n + a)Cn . 𝑥 efWVU
eij
et en dérivant la dérivée première, on a :
h
y′′ = leij(n + a)(n + a − 1) Cn. 𝑥 efWV"
`
En remplaçant ce terme dans l’équation standard de Bessel :
𝑥 " 𝑦" + 𝑥𝑦 ' + (𝑥 " − 𝑎" )𝑦 = 0 ; cela donne :
h h
y′′ = leij(n + a)(n + a − 1) Cn. 𝑥 efW + leij(n + a) Cn. 𝑥 efW
h
+ leij Cn. 𝑥 efWf" - h
leij 𝒂𝟐 Cn. 𝑥 efW = 0
6
On peut simplifier cette équation, en regroupant tous les termes qui ont la puissance xn+a. ; on
aura donc :
h h
leij[(𝑛 + 𝑎)" − 𝑎" ] Cn. 𝑥 efW + leij Cn. 𝑥 efWf" = 0
or
(n+a)2 - a2. est une identité remarquable (a2-b2) = (a - b)(a + b) ; ce qui donne n(n+2a) on a
donc :
h h
leij 𝑛(𝑛 + 2𝑎) Cn. 𝑥 efW + leij Cn. 𝑥 efWf" = 0
Développons ce terme jusqu’à l’ordre 2
Si n=0 è 0
h
Si n = 1 è 0 + 1(1+2a) C1. 𝑥 UfW + lei" 𝑛(𝑛 + 2𝑎) Cn. 𝑥 efW + ….
Pour le deuxième terme on pose k = n+2 è n = k-2 et on remplace « n » par k-2 ; ça devient :
h h
Si n = 1 è 0 + 1(1+2a) C1. 𝑥 UfW + lei" 𝑛(𝑛 + 2𝑎) Cn. 𝑥 efW + loV"ij Ck − 2. 𝑥 ofW = 0
On voit que k-2 donne « k = 2 » comme « n = 2 » ; k est donc une variable muette ; on aurait pu
dire k=2 à ¥ comme j=2 à ¥ ou n=2 à ¥ ou p=2 à ¥ en tant que variable muette ; on va donc
remplacer la lettre « k » par « n » ; ce qui donne
h
Si n = 1 è 0 + 1(1+2a) C1. 𝑥 UfW + l𝒏i𝟐 𝒏(𝒏 + 𝟐𝒂) 𝐂𝐧. 𝒙𝒏f𝒂 + ∑h
𝒏i𝟐 𝐂𝐧 − 𝟐. 𝒙
𝒏f𝒂
=0
On peut donc combiner ces deux séries en base « n = 2 »
h
0 + 1(1+2a) 𝐂𝟏. 𝑥 UfW + lei"[𝑛(𝑛 + 2𝑎) Cn + Cn − 2]. 𝑥 efW = 0
Pour que cette identité soit valide, il faut que le paramètre C1 soit égale à zéro puisque à droite il
n y a pas de terme en 𝒙𝟏f𝒂 ; il est nécessaire que C1 soit égale à zéro mais il faut bien sûr que
n (n + 2a) Cn + Cn – 2 soit aussi égale à zéro.
C1 = 0
* n (n + 2a) Cn + Cn – 2 = 0. "n³2
* relation de récurrence.
On a donc C1 = 0 ; on peut tirer Cn en vertu de la relation de récurrence précédente
n (n + 2a) Cn + Cn – 2 = 0 ; on aura :
𝐶1 = 0. ; n=1
VU
𝐶𝑛 = \e(ef"W)] 𝐶𝑛 − 2; n ³ 2
7
Et on voit immédiatement que tous les termes impairs sont nuls soit
C1 = C3 = C5 = C7 =…= 0
Il n y a que les termes pairs qui sont non-nuls ; alors essayons de calculer quelques-uns de ces
termes qui sont non-nuls
−1 (−1) 1
𝐶2 = F G 𝐶0 = X 2 Z 𝐶0
2(2 + 2𝑎) 2 (1 + 𝑎)
−1 (−1) 2
𝐶4 = F G 𝐶2 = v 4 w 𝐶0
4(4 + 2𝑎) 2 . 1.2(1 + 𝑎)(2 + 𝑎)
−1 (−1) 3
𝐶6 = F G 𝐶4 = v 6 w 𝐶0
6(6 + 2𝑎) 2 . 1.2.3(1 + 𝑎)(2 + 𝑎)(3 + 𝑎)
La forme générale des coefficients pairs va s’écrire :
(−1) 𝑘
𝐶2𝑘 = v 2𝑘 w 𝐶0
2 . 𝑘! ∏𝑘𝑖=1(𝑖 + 𝑎)
Si on pose C0, qui est une constante arbitraire tel que
1
𝐶0 = 𝑎
2 G(𝑎 + 1)
Alors la forme C2k peut s’écrire tout simplement par :
(−1) 𝑘
𝐶2𝑘 =
22𝑘 . 𝑘! (1 + 𝑎)(2 + 𝑎) … (𝑘 + 𝑎) G(𝑎 + 1) 2𝑎
À titre de rappel :
x G(𝑥) = G(𝑥 + 1) et donc x G(𝑘 + 𝑎) = G(𝑘 + 𝑎 + 1) ; C2k devient alors :
(VU) 𝑘
𝐶2𝑘 = 2𝑘+𝑎 (C2k pairs). " k = 0,1,2,3, 4….¥
2 .𝑘! G(𝑎+1)
𝐶2𝑘 + 1 = 0 (C2k impairs)
8
VI. SOLUTION DE L’ÉQUATION DIFFÉRENTIELLE DE BESSEL
Elle s’écrit alors :
h
y = leij Cn . 𝑥 efgU
où r1= a. ; on peut décomposer cette somme pour les entiers pairs et impairs :
h
y = b Cn . 𝑥 efW
eij
y = ∑h h
ei~W3g•( ) + ∑ei3€~W3g•( ) or le 2
ème
terme = 0 ; donc
h h
(−1) • 2 T•‚ƒ
y= loij C2k . 𝑥 "ofW è y=b T•‚ƒ
oij 𝐾! G(…f†fU) . "
on peut simplifier en écrivant :
h
(−1) o 𝑥 "ofW
y=‡ ˆ ‰
𝐾! G(k + a + 1) 2
oij
cette somme en série est appelée Ja (x)
Ja(x) = fonction de Bessel de première espèce d’odre ou d’indice a
Une solution de l’équation de Bessel standard !
9
VII. L’ÉQUATION DE BESSEL DE 2ème ESPECE
Construisons cette équation à partir de la première espèce ; soit
h
(−1) o 𝑥 "ofW
Ja(x) = ‡ ˆ ‰
𝐾! G(k + a + 1) 2
oij
Lorsque a = 0 ; cette fonction vaut alors :
h
(−1) • 2 "o
J0(x) = b T
ˆ ‰ car G(k + 1) = 𝑘!
oij (k!) "
La fonction de Bessel de deuxième espèce et d’ordre a est construite à partir de la
fonction de Bessel de première espèce et de même ordre, on la note :
𝐽𝑎(𝑥 )cos (𝜋𝑎) − 𝐽 − 𝑎(𝑥)
𝑌𝑎 (𝑥 ) =
sin (𝜋𝑎)
Ya(x) = fonction de Bessel de deuxième espèce d’odre ou d’indice a ou fonction de
Neumann
Une solution de l’équation de Bessel standard !
VIII. APPLICATIONS A LA RESOLUTION D 'UN PROBLEME DE
PHYSIQUE MATHEMATIQUE
Exemple 1 : Propagation de la chaleur dans un cylindre infini de rayon R.
Il faut déterminer la température à l’intérieure du cylindre, si on a la distribution de
la température au moment initial et sur la surface du cylindre.
On prend que la température U dépend de la distancer. Soit la condition initiale :
10
11
12
13
Exemple 2 : cylindre infini
14
15
16
17
IX. CONCLUSION
Dans ce travail pratique, nous avons principalement exposé un modeste travail sur
les équations 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 scientifiques,
en particulier la physique.
Nous y avons montré que les équations de Bessel peuvent être modifiées ou trouver
solution avec une première et une deuxième espèce d’ordre ou d’indice « a ».
Ensuite deux exemples concrets sont placés en dernier pour en montrer l’utilité.
Nous laissons la latitude aux autres étudiants d’avenir pour approfondir le sujet et
aboutir jusqu’aux espèces non démontrés ici.
18
X. RÉFÉRENCES
[1] ELIE BELORIZKY, Outils mathématiques à l’usage des scientifiques et
ingénieurs, EDP sciences, 2007.
[2] A. NIKIFOROV, V. OUVAROV, Fonctions spéciales de la physique
mathématique, , Office des Publications universitaires,1983.
[3] J.P. PELLETIER, Techniques Numériques Appliquées au Calcul Scientifique,
Masson et Cie., Paris, 1971.
[4] A. POINSOT, "Exemple de résolution numérique de l'équation caractéristique
d'une cavité hyperfréquence contenant un échantillon diélectrique", Revue de
Physique Appliquée, Tome 12, p. 1055, juillet 1977.
[5] M.R. SPIEGEL, théories et applications de l’analyse, série Schaum-Rensselaer
Polytechnic Institute, McGraw-Hill, Paris, 1973
[6] Y. LEROYER, P. TESSON, mathématiques pour l’ingénieur, Dunod, Paris,
2009
[7] M.CHOSSAT, Y.PRIVAT, aide-mémoire mathématiques de l’ingénieur,
Dunod, Paris, 2017
[8] P. COLMEZ, éléments d’analyse et d’algèbre (et de théorie des nombres), Ed.
Polytechniques, Paris, 2009
[9] F. GASCON, fonctions de Bessel et combinatoires, universités de Québec à
Montréal-thèse de doctorat en maths, Montréal, 2000
[10] H. MERAD, fonctions de Bessel et applications a la résolution d 'un problème
de physique mathématique- thèse de master en maths, université Larbi ben Mhidi
d’Oum el Bouaghi, Alger, 2019
[11] M. PINARD, convertisseurs et électronique de puissance, Dunod, Paris, 2007
