CORRIGES TYPE (Deuxième Devoir / INSPEI / Abomey) : 2ème méthode : Une double intégration par parties fournit,

𝟏 𝑠ℎ(𝜆𝑥) 𝜋 𝒏 𝜋
Exercice 1 : 𝑎𝑛 (𝑔) = 𝝅 ([ cos(𝑛𝑥)] + 𝜆 ∫−𝜋 𝑠ℎ(𝜆𝑥)sin(𝑛𝑦)𝑑𝑥)
𝜆 −𝜋
1-) Calculons la série de Fourier trigonométrique réelle 2(−1)𝑛 𝑠ℎ(𝜆𝜋) 𝑛2
= − 𝜆2 𝑎𝑛 (𝑔) Et donc ∀ 𝑛 𝜖 ℕ,
𝜆𝜋
de chacune des fonctions 2𝜋 − périodique suivantes :
𝟐(−𝟏)𝒏 𝒔𝒉(𝝀𝝅) 𝝀𝟐 𝟐𝝀𝒔𝒉(𝝀𝝅) (−𝟏)𝒏
 𝒇(𝒙) = (𝒙 − 𝝅)𝟐 sur [𝟎 , 𝟐𝝅] 𝒂𝒏 (𝒈) = 𝝀𝝅
∙ 𝒏𝟐 +𝝀𝟐 = 𝝅
∙ 𝒏𝟐 +𝝀𝟐
On remarque que la fonction 𝒇 est paire, de sorte que 𝒃𝒏 (𝒇) = 𝟎 pour tout 𝒏. 1 2𝜋 1 𝜋 𝟐𝒔𝒉(𝝀𝝅)
Or 𝒂𝟎 (𝒈) = 𝜋 ∫0 𝑐ℎ(𝜆𝑥) 𝑑𝑥 = 𝜋 ∫−𝜋 𝑐ℎ(𝜆𝑥) 𝑑𝑥 =
𝝅 𝝀𝝅
𝟏 𝟐𝝅 𝟏 𝝅 𝟏 𝒚𝟑 𝟐𝝅𝟐
Par ailleurs, 𝒂𝟎 (𝒇) = ∫𝟎 (𝒙 − 𝝅)𝟐 𝒅𝒙 = ∫−𝝅 𝒚𝟐 𝒅𝒚 [ ] = 𝒔𝒉(𝝀𝝅) 𝟐𝝀𝒔𝒉(𝝀𝝅) (−𝟏)𝒏
𝝅 𝝅 𝝅 𝟑 𝟑
−𝝅 D’où 𝑺𝑮(𝒈)(𝒙) = + ∑+∞
𝒏=𝟏 𝒄𝒐𝒔(𝒏𝒙)
𝝀𝝅 𝝅 𝒏𝟐 +𝝀𝟐
et, pour tout 𝒏 ≥ 𝟏,
1 2𝜋 1 𝜋  Exponentielle complexe de 𝒉(𝒙) = 𝒆𝒙 sur ] − 𝝅 , 𝝅]
2
𝑎𝑛 (𝑓) = ∫ (𝑥 − 𝜋) cos(𝑛𝜋)𝑑𝑥 = ∫ 𝑦 2 cos(𝑛𝑦 + 𝑛𝜋)𝑑𝑦 𝜋
𝜋 0 𝜋 −𝜋
1 𝜋 𝑥 −𝑖𝑛𝑥 1 𝜋 (1−𝑖𝑛)𝑥 1 𝑒 (1−𝑖𝑛)𝑥
(−𝟏)𝒏 𝝅 (−𝟏)𝒏 𝟒 𝟒 𝑐𝑛 (ℎ) = ∫ 𝑒 𝑒 𝑑𝑥 = ∫ 𝑒 𝑑𝑥 = [ ]
𝒂𝒏 (𝒇) = 𝝅
∫−𝝅 𝒚𝟐 𝐜𝐨𝐬(𝒏𝒚)𝒅𝒚 = 𝝅
∙ 𝒏𝟐 (−𝟏)𝟐 𝝅 = 𝒏𝟐 2𝜋 −𝜋 2𝜋 −𝜋 2𝜋 1 − 𝑖𝑛 −𝜋

où l'on a effectué deux intégrations par parties. La série de Fourier de s'écrit 1 𝑒 (1−𝑖𝑛)𝜋 − 𝑒 −(1−𝑖𝑛)𝜋 1 𝑒 𝜋 − 𝑒 −𝜋 𝑠ℎ𝜋 (−1)𝑛
= ( )= ((−1)𝑛 )= ∙
2𝜋 1 − 𝑖𝑛 2𝜋 1 − 𝑖𝑛 𝜋 1 − 𝑖𝑛
𝝅𝟐 𝐜𝐨𝐬(𝒏𝒙)
donc : 𝑺𝑭(𝒇)(𝒙) = + 𝟒 ∑+∞ 𝑠ℎ𝜋
𝟑 𝒏=𝟏 𝒏𝟐 𝑺𝑯(𝒉)(𝒙) = 𝒄𝟎 (𝒉) + ∑+∞
𝒏≥𝟏[𝒄𝒏 (𝒉)𝒆
𝒊𝒏𝒕
+ 𝒄−𝒏 (𝒉)𝒆−𝒊𝒏𝒕 ] avec 𝑐0 (ℎ) = 𝜋

 𝒈(𝒙) = 𝒄𝒉(𝝀𝒙) sur ] − 𝝅 , 𝝅] 2-) - Déterminons le déplacement des points de cette corde par rapport à l’axe
La fonction 𝒈 est 2𝜋 −périodique, continue par morceaux sur ℛ et paire. des abscisses si les vitesses initiales sont nulles.
𝟏 𝝅
Pour 𝑛 𝜖 ℕ∗ , 𝒃𝒏 (𝒈) = 𝟎, puis pour 𝑛 𝜖 ℕ, 𝒂𝒏 (𝒈) = ∫ 𝒄𝒉(𝝀𝒙)𝐜𝐨𝐬(𝒏𝒙) 𝒅𝒙 La 𝜕2 𝑈 𝜕2 𝑈
𝝅 −𝝅 solution de l’équation = 𝑎2 𝜕𝑥 2 (1) constitue le déplacement des
𝜕𝑡 2
1ère Méthode : Soit 𝑛 𝜖 ℕ, 𝜕𝑈
points de la corde. | = 0 = 𝜑(𝑥) ; 𝑈(0, 𝑡) = 𝑈(𝐿, 𝑡) = 0 ; 𝑈(𝑥, 0) =
1 𝜋 1 𝜋 𝜋 𝜕𝑡 𝑡=0
𝑎𝑛 (𝑔) = ℜ𝑒(∫−𝜋 𝑐ℎ(𝜆𝑥)𝑒 𝑖𝑛𝑥 ) = ℜ𝑒(∫−𝜋 𝑒 (𝜆+𝑖𝑛)𝑥 𝑑𝑥 + ∫−𝜋 𝑒 (−𝜆+𝑖𝑛)𝑥 𝑑𝑥)
𝜋 2𝜋
𝑓(𝑥). Par la méthode de séparation des variables, on a :
(𝜆+𝑖𝑛)𝜋 −(𝜆+𝑖𝑛)𝜋 (−𝜆+𝑖𝑛)𝜋 −(−𝜆+𝑖𝑛)𝜋
1 𝑒 −𝑒 𝑒 −𝑒 𝑇 ′′ (𝑡) 𝑋 ′′ (𝑥)
= ℜ𝑒 ( + ) 𝑈(𝑥, 𝑡) = 𝑋(𝑥) ∙ 𝑇(𝑡) (2). (2) dans (1) donne =
2𝜋 𝜆 + 𝑖𝑛 −𝜆 + 𝑖𝑛 𝑎2 𝑇(𝑡) 𝑋(𝑥)

(−1)𝑛 2𝑠ℎ(𝜆𝜋) −2𝑠ℎ(𝜆𝜋) 𝟐𝒔𝒉(𝝀𝝅) (−𝟏)𝒏 Ces conditions ne peuvent être vérifiées que si elles sont égales
= ℜ𝑒 ( + )= ∙
2𝜋 𝜆+𝑖𝑛 −𝜆+𝑖𝑛 𝝅 𝒏𝟐 +𝝀𝟐
à une constante. Posons cette constante égales à – 𝜆2 (𝜆 > 0) on
Page 1 sur 6
 𝑇 ′′ (𝑡) + 𝑎2 𝜆2 𝑇(𝑡) = 0 𝐿 𝐿
a :{ 2 𝑛𝜋 2 4ℎ 𝑛𝜋
𝑋 ′′ (𝑥) + 𝜆2 𝑋(𝑥) = 0 𝐴𝑛 = ∫ 𝑓(𝑥)𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 = ∫ 2 𝑥(𝐿 − 𝑥)𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥
𝐿 𝐿 𝐿 𝐿 𝐿
0 0
𝑟 2 + 𝑎2 𝜆2 = 0 ⟹ 𝑟 = ± 𝑖𝑎𝜆
Equations caractéristique : { 𝐿 𝐿
𝑟 2 + 𝜆2 = 0 ⟹ 𝑟 = ± 𝑖𝜆 8ℎ 𝑛𝜋 8ℎ 𝑛𝜋
𝑇(𝑡) = 𝐴𝑐𝑜𝑠(𝑎𝜆𝑡) + 𝐵𝑠𝑖𝑛(𝑎𝜆𝑡) 𝐴𝑛 = 2 ∫ 𝑥𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥 − 3 ∫ 𝑥 2 𝑠𝑖𝑛 ( 𝑥) 𝑑𝑥
Les solutions sont sous la forme : { 𝐿 𝐿 𝐿 𝐿
𝑋(𝑥) = 𝐶𝑐𝑜𝑠(𝜆𝑥) + 𝐷𝑠𝑖𝑛(𝜆𝑥) 0 0

Une double intégration par partie donne :

𝑈(0, 𝑡) = 𝑋(0) = 𝑈(𝐿, 𝑡) = 𝑋(𝐿) = 0 alors 𝐶 = 0 et 𝐷𝑠𝑖𝑛(𝜆𝐿) = 0
16ℎ 16ℎ
𝑛𝜋 𝑛𝜋 𝐴𝑛 = 𝜋3 𝑛3 [1 − cos(𝑛𝜋)] = 𝜋3 𝑛3 [1 − (−1)𝑛 ] D’où
⟹ 𝜆𝐿 = 𝑛𝜋 (𝑛 ∈ ℕ) donc 𝜆 = soit 𝑋(𝑥) = 𝐷𝑠𝑖𝑛 ( 𝐿 𝑥)
𝐿
∞
𝜕𝑇(𝑡) 𝟏𝟔𝒉 [𝟏 − (−𝟏)𝒏 ] 𝒏𝝅𝒂 𝒏𝝅
= −𝑎𝜆𝐴𝑠𝑖𝑛(𝑎𝜆𝑡) + 𝑎𝜆𝐵𝑐𝑜𝑠(𝜆𝑥) 𝑼(𝒙, 𝒕) = 𝟑 ∑ 𝒄𝒐𝒔 ( 𝒕) 𝒔𝒊𝒏 ( 𝒙)
𝜕𝑈 𝜕𝑇 𝜕𝑡 𝝅 𝒏𝟑 𝑳 𝑳
| = | = 0 ⟹ { 𝜕𝑇(𝑡) 𝒏=𝟏
𝜕𝑡 𝜕𝑡
𝑡=0 𝑡=0 | = 𝑎𝜆𝐵 = 0 ⟹ 𝐵 = 0
𝜕𝑡 𝑡=0  Pour 𝑛 = 2𝑘 on a : 𝐴2𝑘 = 0 ⟹ 𝑼(𝒙, 𝒕) = 𝟎
𝑛𝜋𝑎
Soit 𝑇(𝑡) = 𝐴𝑐𝑜𝑠 ( 𝑡)  Pour 𝑛 = 2𝑘 + 1 on a : 𝐴2𝑘+1 = 𝜋3 (2𝑘+1)3 ⟹
32ℎ
𝐿

La fonction 𝑈𝑛 (𝑥, 𝑡) est donc égale à : ∞

𝟑𝟐𝒉 𝟏 (𝟐𝒌 + 𝟏)𝝅𝒂 (𝟐𝒌 + 𝟏)𝝅
𝑛𝜋𝑎 𝑛𝜋 𝑼(𝒙, 𝒕) = 𝟑 ∑ 𝟑
𝒄𝒐𝒔 ( 𝒕) 𝒔𝒊𝒏 ( 𝒙)
𝑈𝑛 (𝑥, 𝑡) = 𝐴𝑛 𝑐𝑜𝑠 ( 𝑡) 𝑠𝑖𝑛 ( 𝐿 𝑥) 𝝅 (𝟐𝒌 + 𝟏) 𝑳 𝑳
𝐿 𝒌=𝟎

Si 𝑈𝑛 (𝑥, 𝑡) est solution alors la somme ∑∞

𝑛=0 𝑈𝑛 (𝑥, 𝑡) est aussi - Ecrivons jusqu’à sept termes la série obtenue.
solution, donc : 𝑈(𝑥, 𝑡) = ∑∞
𝑛=0 𝑈𝑛 (𝑥, 𝑡) Dans ce cas, il suffit d’écrire l’expression sommative de 𝑼(𝒙, 𝒕)
∞
𝒏𝝅𝒂 𝒏𝝅 pour 𝒌 = {𝟎; 𝟏; 𝟐; 𝟑; 𝟒; 𝟓; 𝟔}.
𝑼(𝒙, 𝒕) = ∑ 𝑨𝒏 𝒄𝒐𝒔 ( 𝒕) 𝒔𝒊𝒏 ( 𝒙)
𝒏=𝟏
𝑳 𝑳 3-) En utilisant la transformation de Fourier, calculons 𝝋(𝒙, 𝒕) en admettant
𝑛𝜋 4ℎ qu’à l’instant initial toutes les particules diffusantes sont concentrées à
𝑈(𝑥, 0) = 𝑓(𝑥) ⟹ 𝑓(𝑥) = ∑∞
𝑛=0 𝐴𝑛 𝑠𝑖𝑛 ( 𝐿 𝑥) = 𝑥(𝐿 − 𝑥)
𝐿2
l’origine soit 𝝋(𝒙, 𝟎) = 𝜹(𝒙), où 𝜹(𝒙) est la fonction de Dirac.
4ℎ
La fonction 𝑓(𝑥) = 𝑥(𝐿 − 𝑥) étant continue positive et dérivable 𝜕𝜑 𝜕2 𝜑 ̂
𝜕𝜑
𝐿2 Posons 𝑇𝐹(𝜑) = 𝜑̂ donc = 𝐷 𝜕𝑥 2 ⟹ = (𝑖𝑘)2 𝐷𝜑̂ = −𝑘 2 𝐷𝜑̂
𝜕𝑡 𝜕𝑡
donc elle peut être développée en série de Fourier. D’où ̂
𝜕𝜑 2 𝑡𝐷
⟹ ̂
= −𝑘 2 𝐷𝜕𝑡 soit ln(𝜑̂) = −𝑘 2 𝑡𝐷 + 𝑐𝑠𝑡 ⟹ 𝜑̂(𝑘, 𝑡) = 𝑐𝑠𝑡𝑒 −𝑘
𝜑
𝟐 𝒕𝑫
or à 𝑡 = 0, 𝜑̂(𝑘, 0) = 𝑐𝑠𝑡 alors 𝝋 ̂(𝒌, 𝟎)𝒆−𝒌
̂(𝒌, 𝒕) = 𝝋
Page 2 sur 6
 −𝑥2 −𝑘2 𝜹2 2 ℓ 𝑛𝜋
Si 𝑓(𝑥) = 𝑒 2𝜹𝟐 ⟹ 𝑓(𝑘) = 𝛿𝑒 2 ; ⟹ 𝐵𝑛 = ℓ ∫0 𝑓(𝑥)𝑠𝑖𝑛 ( ℓ 𝑥) 𝑑𝑥
2 𝑡𝐷 1 2 𝑡𝐷 ℓ
comme 𝑒 −𝑘 = (√2𝑡𝐷𝑒 −𝑘 ) donc 2 𝑛𝜋 2
⟹ 𝐵𝑛 = ℓ ∫02 𝑓(𝑥)𝑠𝑖𝑛 ( ℓ 𝑥) 𝑑𝑥 + ℓ ∫ℓ 𝑓(𝑥)𝑠𝑖𝑛 ( ℓ 𝑥) 𝑑𝑥
ℓ 𝑛𝜋
√2𝑡𝐷
2
−𝑘 2 𝑡𝐷 1 −𝑥 2
à𝑒 ⟶ et à 𝜑̂(𝑘, 0) ⟶ 𝜑(𝑥, 0) = 𝛿(𝑥) or 2
ℓ
𝑛𝜋 2 ℓ 𝑛𝜋
√2𝑡𝐷 4𝑡𝐷 ⟹ 𝐵𝑛 = ℓ ∫0 𝑥𝑠𝑖𝑛 ( ℓ 𝑥) 𝑑𝑥 + ℓ ∫ℓ (ℓ − 𝑥)𝑠𝑖𝑛 ( ℓ 𝑥) 𝑑𝑥
2

1 +∞ 2
𝑓̂(𝑘)𝑔̂(𝑘) ⟶ ∫ 𝑓(𝑦 − 𝑥)𝑔(𝑥)𝑑𝑥 donc
√2𝜋 _∞ 4ℓ 𝑛𝜋
Une double intégration par partie donne : 𝐵𝑛 = (𝑛𝜋)2 𝑠𝑖𝑛 ( 2 )
−(𝑥−𝑦)2
1 1 +∞
𝜑(𝑥, 𝑡) = ∫ 𝑒 4𝑡𝐷 𝛿(𝑦)𝑑𝑦 ou ∞
√2𝜋 √2𝑡𝐷 _∞
𝟒𝓵 𝒏𝝅 −𝑫𝒕(𝒏𝝅)𝟐 𝒏𝝅
−𝑦2 𝝋(𝒙, 𝒕) = ∑ 𝒔𝒊𝒏 ( ) 𝒆 𝓵 𝒔𝒊𝒏 ( 𝒙)
𝜑(𝑥, 𝑡) =
1 1 +∞
− 𝑦)𝑒 4𝑡𝐷 𝑑𝑦 et comme (𝒏𝝅)𝟐 𝟐 𝓵
∫ 𝛿(𝑥
√2𝜋 √2𝑡𝐷 _∞ 𝒏=𝟏

+∞
∫_∞ 𝛿(𝑎 − 𝑦)𝑓(𝑦)𝑑𝑦 = 𝑓(𝑎).  Pour 𝑛 = 2𝑘 on a : 𝐵2𝑘 = 0 ⟹ 𝝋(𝒙, 𝒕) = 𝟎
4ℓ(−1)𝑘
𝟏 −𝒙𝟐  Pour 𝑛 = 2𝑘 + 1 on a : 𝐵2𝑘+1 = 𝜋2(2𝑘+1)2 ⟹
D’où 𝝋(𝒙, 𝒕) = 𝒆 𝟒𝒕𝑫 avec 𝑡 > 0
√𝟒𝝅𝒕𝑫
∞
𝟒𝓵(−𝟏)𝒌 𝒏𝝅 𝟐
−𝑫𝒕( ) 𝒏𝝅
4-) Trouvons la solution de cette équation de diffusion : 𝝋(𝒙, 𝒕) = ∑ 𝟐 𝟐
𝒆 𝓵 𝒔𝒊𝒏 ( 𝒙)
𝝅 (𝟐𝒌 + 𝟏) 𝓵
𝝏𝝋 𝝏𝟐 𝝋 𝒌=𝟎
= 𝑫 𝝏𝒙𝟐 satisfaisant aux conditions initiales et aux conditions
𝝏𝒕
Exercice 2:
aux limites données.
1-) Résolvons l’équation d’Euler : 𝒙𝟐 𝒚′′ + 𝒙𝒚′ − 𝒌𝟐 𝒚 = 𝟎 (𝑘 ∈ ℜ) par la
En procédant toujours la méthode de séparation des variables,
méthode de Frobenius pour 𝒌 ≠ 𝟎 et pour 𝒌 = 𝟎.
on trouve la solution générale comme suit:
 Solution pour 𝒌 ≠ 𝟎
∞
𝒏𝝅 𝟐 𝒏𝝅
𝝋(𝒙, 𝒕) = ∑ 𝑩𝒏 𝒆−𝑫( 𝓵 ) 𝒕 𝒔𝒊𝒏 ( 𝒙) La solution de cette équation est sous la forme :𝑦(𝑥) = ∑∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟
𝓵
𝒏=𝟏 ⟹ 𝑥𝑦′ = ∑∞
𝑛=0 𝑐𝑛 (𝑛 + 𝑟)𝑥
𝑛+𝑟
et 𝑥 2 𝑦′′ = ∑∞
𝑛=0 𝑐𝑛 (𝑛 + 𝑟)(𝑛 + 𝑟 − 1)𝑥
𝑛+𝑟

𝑛𝜋
𝜑(𝑥, 0) = ∑∞
𝑛=1 𝐵𝑛 𝑠𝑖𝑛 ( ℓ 𝑥) = 𝑓(𝑥) Soit ∑∞ 2 𝑛+𝑟
𝑛=0 𝑐𝑛 [(𝑛 + 𝑟)(𝑛 + 𝑟 − 1) + (𝑛 + 𝑟) − 𝑘 ] 𝑥 =0
La fonction 𝑓(𝑥) étant continue positive et dérivable donc elle Pour l’équation indicielle, on a : 𝑛 = 0 ⟹ 𝐶 [𝑟(𝑟 − 1) + 𝑟 − 𝑘 2 ] = 0
0
peut être développée en série de Fourier. D’où
𝐶0 = 1 et 𝑟 2 − 𝑘 2 = 0 ⟹ 𝑟 = ±𝑘

Page 3 sur 6
 La solution générale de l’équation d’Euler est donc : 1
𝐸 ∑∞
𝑛=2 𝑐𝑛−2 𝑥
𝑛+𝑟−1
− 4 𝐹 ∑∞
1
𝑛=3 𝑐𝑛−3 𝑥
𝑛+𝑟−1
=0
2
𝒚(𝒙) = 𝑨𝒙𝒌 + 𝑩𝒙−𝒌 avec 𝐴 et 𝐵 des constants réel. 𝟏 𝟏 𝟏
⟹ 𝒄𝒏 = 𝒏(𝒏+𝒎) (−𝜶𝒄𝒏−𝟏 − 𝟐 𝑬𝒄𝒏−𝟐 + 𝟒 𝑭𝒄𝒏−𝟑 ) or
 Solution pour 𝒌 = 𝟎
𝑚2 𝑚
Pour 𝒌 = 𝟎, l’équation différentielle devient : 𝑐0 𝑟 2 −= 0 𝑐𝑜𝑚𝑚𝑒 𝑟 = 𝑑𝑜𝑛𝑐 𝑐0 = 1
4 2
𝑥 2 𝑦 ′′ + 𝑥𝑦 ′ = 0 ⟹
𝑑
(𝑥𝑦′) = 0 donc la solution générale de cette −𝛼
𝑑𝑥 𝑐1 =
𝑚+1
équation d’Euler est : 𝒚(𝒙) = 𝑫𝒍𝒏𝒙 + 𝑪 avec 𝐷 et 𝐶 des constants réel. 1 𝛼2 𝐸
⟹ 𝑐2 = 2(𝑚 + 2) (𝑚 + 1 − 2 )
2-) Trouvons la solution dans le voisinage du point 𝑥 = 0 sous la forme
1 1 1
d’une série (donnons les quatre premiers coefficients non nuls). 𝑐3 = (−𝛼𝑐2 − 𝐸𝑐1 + 𝐹𝑐0 )
{ 3(𝑚 + 3) 2 4
𝑑 1 𝑚2 1
(𝑥𝑦′) + [ 𝐸𝑥 + 𝛼 − − 4 𝐹𝑥 2 ] 𝑦 = 0 ⟹ 𝒎
𝑑𝑥 2 4𝑥 D’où 𝒚(𝒙) = (𝟏 + 𝒄𝟏 𝒙 + 𝒄𝟐 𝒙𝟐 + 𝒄𝟑 𝒙𝟑 + ⋯ ) 𝒙 𝟐
1 1 1 𝑚2 1 𝑚2
𝑦 ′′ + 𝑥 𝑦 ′ + [2 𝐸 + 𝑥 𝛼 − 4𝑥 2 − 4 𝐹𝑥] 𝑦 = 0 ⟹ 𝑎0 = 1 et 𝑏0 = − 3-) Trouvons les solutions générales des équations de Bessel suivantes :
4
𝑚2 𝟗 3
L’équation indicielle 𝑟(𝑟 − 1) + 𝑎0 𝑟 + 𝑏0 = 0 s’écrit : 𝑟 2 − = 0 soit a-) 𝒙𝟐 𝒚′′ + 𝒙𝒚′ + (𝒙𝟐 − 𝟒) 𝒚 = 𝟎 C’est l’équation de Bessel d’ordre .
4 2

𝑟=±
𝑚
donc la plus grande des racines de l’équation D’où la solution générale est :
2
𝑚 𝒚(𝒙) = 𝑨𝑱𝟑 (𝒙) + 𝑩𝑱− 𝟑 (𝒙) avec 𝑱𝟑 (𝒙) la fonction de Bessel de 1ère
déterminante est 𝑟 = 2 . Posons 𝑦(𝑥) = ∑∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟
donc 𝟐 𝟐 𝟐

3
𝑦′ = ∑∞
𝑛=0 𝑐𝑛 (𝑛 + 𝑟)𝑥
𝑛+𝑟−1
et 𝑦′′ = ∑∞
𝑛=0 𝑐𝑛 (𝑛 + 𝑟)(𝑛 + 𝑟 − 1)𝑥
𝑛+𝑟−2
espèce d’ordre 2 ; 𝐴 et 𝐵 des constantes réelles.
L’équation devient : b-) 𝒙𝟐 𝒚′′ + 𝒙𝒚′ + 𝒙𝟐 𝒚 = 𝟎 ⟹ 𝒙𝟐 𝒚′′ + 𝒙𝒚′ + (𝒙𝟐 − 𝟎)𝒚 = 𝟎 C’est
1
∑∞ 2 𝑛+𝑟−1
𝑛=0 𝑐𝑛 (𝑛 + 𝑟) 𝑥 + 2 𝐸 ∑∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟+1
+ 𝛼 ∑∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟
− l’équation de Bessel d’ordre 0 . Ici, 𝑝 = 0 entier. D’où la solution
𝑚2 1 générale est : 𝒚(𝒙) = 𝑨𝑱𝟎 (𝒙) + 𝑩𝒀 𝟎 (𝒙) avec 𝑱𝟎 (𝒙) la fonction de Bessel
∑∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟−1
− 4 𝐹 ∑∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟+2
=0
4
1 de 1ère espèce d’ordre 0 et 𝒀 𝟎 (𝒙) la fonction de Bessel de 2ème espèce
⟹ ∑∞ 2 𝑛+𝑟−1
𝑛=0 𝑐𝑛 (𝑛 + 𝑟) 𝑥 + 2 𝐸 ∑∞
𝑛=2 𝑐𝑛−2 𝑥
𝑛+𝑟−1
+ 𝛼 ∑∞
𝑛=1 𝑐𝑛−1 𝑥
𝑛+𝑟−1
−
d’ordre 0 ; 𝐴 et 𝐵 des constantes réelles.
𝑚2 1
∑∞
𝑛=0 𝑐𝑛 𝑥
𝑛+𝑟−1
− 𝐹 ∑∞
𝑛=3 𝑐𝑛−3 𝑥
𝑛+𝑟−1
=0 𝟏 1
4 4 c-) 𝒙𝟐 𝒚′′ + 𝒙𝒚′ + (𝒙𝟐 − 𝟒) 𝒚 = 𝟎 C’est l’équation de Bessel d’ordre 2 .
𝑚2
⟹ ∑∞ 2
𝑛=0 𝑐𝑛 [(𝑛 + 𝑟) − ] 𝑥 𝑛+𝑟−1 + 𝛼 ∑∞
𝑛=1 𝑐𝑛−1 𝑥
𝑛+𝑟−1
+
4
Page 4 sur 6
 (𝑎−1)2 −4𝑏 1
D’où la solution générale est : 𝒚(𝒙) = 𝑨𝑱𝟏 (𝒙) + 𝑩𝑱− 𝟏 (𝒙) avec devient 𝑡 2 𝑢′′ + 𝑡𝑢′ + (𝑡 2 − 𝑝2 )𝑢 = 0 avec 𝑝2 = = 9 soit
𝟐 𝟐 𝑚2
1 1
2 2 𝑡 2 𝑢′′ + 𝑡𝑢′ + (𝑡 2 − 9) 𝑢 = 0. Ici, 𝑝 = 3 non entier. D’où la solution
𝐽1 (𝑥) = √𝜋𝑥 sin 𝑥 et 𝐽−1 (𝑥) = √𝜋𝑥 cos 𝑥 ; 𝐴 et 𝐵 des constantes réelles.
2 2
générale est : 𝒚(𝒙) = 𝑨𝑱𝟏 (𝟐𝒙) + 𝑩𝑱− 𝟏 (𝟐𝒙) avec 𝑱𝟏 (𝟐𝒙) la fonction de
4-) Ramenons chacune des équations suivantes à l’équation 𝟑 𝟑 𝟑

de Bessel et trouvons leurs solutions générales : 1

Bessel de 1ère espèce d’ordre 3 ; 𝐴 et 𝐵 des constantes réelles.
𝒙𝟐 𝒚′′ − 𝟑𝒙𝒚′ + (𝒙𝟒 − 𝟏𝟐)𝒚 = 𝟎. Cette équation doit avoir la forme:
Exercice 3:
2 ′′ ′
𝑥 𝑦 + 𝑎𝑥𝑦 + (𝑏 + 𝑐𝑥 𝑚 )𝑦
= 0. Par identification, on a : 𝑎 = −3 ;
1-) Montrons que la fonction de Bessel 𝑱𝒑 (𝒙) vérifie les relations de
−𝛼 1
𝑡 𝛽 𝑡 𝛽
𝑏 = −12 ; 𝑐 = 1 et 𝑚 = 4 . Posons 𝑦 = 𝑢 (𝛾) ; 𝑥 = (𝛾) où 𝛼 = récurrence :
𝟏 𝒅 𝒎
𝑎−1
= −2 ; 𝛽 =
𝑚
= 2 et 𝛾 =
2√𝑐 1
= 2 donc on a : 𝑦 = 2𝑡𝑢(𝑡) et 𝑥 = a-) (𝒙 𝒅𝒙) [𝒙𝒑 𝑱𝒑 (𝒙)] = 𝒙𝒑−𝒎 𝑱𝒑−𝒎 (𝒙)
2 2 𝑚
1
(2𝑡)2 . Après des dérivés suivant les changements de variables,
𝟏 𝒅 𝒎
(𝑎−1)2 −4𝑏 b-) (𝒙 𝒅𝒙) [𝒙−𝒑 𝑱𝒑 (𝒙)] = (−𝟏)𝒎 𝒙−𝒑−𝒎 𝑱𝒑+𝒎 (𝒙)
l’équation devient 𝑡 2 𝑢′′ + 𝑡𝑢′ + (𝑡 2 − 𝑝2 )𝑢 = 0 avec 𝑝2 = =
𝑚2

4 soit 𝑡 2 𝑢′′ + 𝑡𝑢′ + (𝑡 2 − 4)𝑢 = 0. Ici, 𝑝 = 2 entier.

c-) Si 𝒓𝟐 = 𝒙𝟐 + 𝒚𝟐 , alors 𝑱𝟎 (𝒓) = ∑+∞
𝒌=−∞ 𝑱𝟐𝒌 (𝒙) ∙ 𝑱𝟐𝒌 (𝒚)
𝒙𝟐 𝒙𝟐
D’où la solution générale est : 𝒚(𝒙) = 𝒙𝟐 [𝑨𝑱𝟐 ( 𝟐 ) + 𝑩𝒀 𝟐 ( 𝟐 )] avec

𝒙𝟐 𝒙𝟐
𝑱𝟐 ( 𝟐 ) la fonction de Bessel de 1ère espèce d’ordre 2 et 𝒀 𝟐 ( 𝟐 ) la fonction 𝒙 𝟏
d-) 𝑬𝒙𝒑 [𝟐 (𝒕 − 𝒕 )] = ∑+∞ 𝒏
𝒑=−∞ 𝒕 𝑱𝑷 (𝒙)

de Bessel de 2ème espèce d’ordre 2 ; 𝐴 et 𝐵 des constantes réelles.

𝟏
𝒙𝟐 𝒚′′ + 𝒙𝒚′ + (𝟒𝒙𝟐 − 𝟗) 𝒚. De façon analogue, on a : 𝑎 = 1 ; 𝑏 = 2-) a-) Calculons les intégrales : 𝐼 = ∫ 𝑥 5 𝐽2 (𝑥)𝑑𝑥 = ∫ 𝑥 2 [𝑥 3 𝐽2 (𝑥)]𝑑𝑥

1 𝑡
−𝛼
𝛽 𝑡 𝛽
1
𝑎−1 𝐼 = ∫ 𝑥 2 𝑑[𝑥 3 𝐽3 (𝑥)] = 𝑥 5 𝐽3 (𝑥) − 2 ∫ 𝑑[𝑥 4 𝐽4 (𝑥)]. Donc
− 9 ; 𝑐 = 4 et 𝑚 = 2. Posons 𝑦 = 𝑢 (𝛾) ;𝑥= (𝛾) où 𝛼 = =0 ;
2
𝑰 = 𝒙𝟓 𝑱𝟑 (𝒙) − 𝟐𝒙𝟒 𝑱𝟒 (𝒙)
𝑚 2√𝑐 𝑡
𝛽= = 1 et 𝛾 = = 2 donc on a : 𝑦 = 𝑢(𝑡) et 𝑥 = 2. Après des
2 𝑚 𝐾 = ∫ 𝑥 6 𝐽3 (𝑥)𝑑𝑥 = ∫ 𝑥 2 [𝑥 4 𝐽3 (𝑥)]𝑑𝑥 = ∫ 𝑥 2 𝑑[𝑥 4 𝐽4 (𝑥)]
dérivés suivant les changements de variables, l’équation
𝐾 = 𝑥 6 𝐽4 (𝑥) − 2 ∫ 𝑑[𝑥 5 𝐽5 (𝑥)] Donc 𝑲 = 𝒙𝟔 𝑱𝟒 (𝒙) − 𝟐𝒙𝟓 𝑱𝟓 (𝒙)
Page 5 sur 6
 𝐿 = ∫ 𝑥 5 𝐽0 (𝑥)𝑑𝑥 = ∫ 𝑥 4 [𝑥𝐽0 (𝑥)]𝑑𝑥 = ∫ 𝑥 4 𝑑[𝑥𝐽1 (𝑥)] 2 2
𝐽1 (𝑥) = √𝜋𝑥 sin 𝑥 et 𝐽−1⁄ (𝑥) = √𝜋𝑥 cos 𝑥
2
𝐿 = 𝑥 5 𝐽1 (𝑥) − 4 ∫ 𝑥 4 𝐽1 (𝑥)𝑑𝑥 De même ∫ 𝑥 4 𝐽1 (𝑥)𝑑𝑥 = 𝑥 4 𝐽2 (𝑥) − 2𝑥 3 𝐽3 (𝑥) 2

𝟐 𝟑−𝒙𝟐 𝟑
D’où 𝑳 = 𝒙𝟓 𝑱𝟏 (𝒙) − 𝟒𝒙𝟒 𝑱𝟐 (𝒙) + 𝟖𝒙𝟑 𝑱𝟑 (𝒙) 𝑱𝟓 (𝒙) = √𝝅𝒙 [( ) 𝐬𝐢𝐧 𝒙 − 𝒙 𝐜𝐨𝐬 𝒙]
𝒙𝟐
D’où 𝟐

b-) Exprimons chacun des résultats trouvés en fonction de 𝐽0 et 𝐽1 . 𝟐 𝒙𝟒 −𝟒𝟓𝒙𝟐 +𝟏𝟎𝟓 𝟏𝟎𝒙𝟐 −𝟕𝟎
𝑱𝟗 (𝒙) = √𝝅𝒙 [( ) 𝐬𝐢𝐧 𝒙 + ( ) 𝐜𝐨𝐬 𝒙]
2𝑝 { 𝟐
𝒙𝟒 𝒙𝟑
Par définition : 𝐽𝑝+1 (𝑥) = 𝐽𝑝 (𝑥) − 𝐽𝑝−1 (𝑥). Variation de p conduit à :
𝑥
2
𝐽2 (𝑥) = 𝐽 (𝑥) − 𝐽0 (𝑥)
𝑥 1
8 − 𝑥2 4
𝐽3 (𝑥) = ( 2 ) 𝐽1 (𝑥) − 𝐽0 (𝑥)
𝑥 𝑥
2
48 − 8𝑥 𝑥 2 − 24
𝐽4 (𝑥) =( ) 𝐽1 (𝑥) + ( ) 𝐽0 (𝑥)
𝑥3 𝑥2
𝑥 4 − 72𝑥 2 + 384 12𝑥 2 − 192
𝐽5 (𝑥) = ( ) 𝐽1 (𝑥) + ( ) 𝐽0 (𝑥)
{ 𝑥4 𝑥3

⟹ 𝑰 = (𝟐𝟒𝒙𝟑 − 𝒙𝟓 − 𝟗𝟔𝒙)𝑱𝟏 − (𝟔𝒙𝟒 − 𝟒𝟖𝒙𝟐 )𝑱𝟎

⟹ 𝑲 = (𝟏𝟗𝟐𝒙𝟑 − 𝟏𝟎𝒙𝟓 − 𝟕𝟔𝟖𝒙)𝑱𝟏 + (𝒙𝟔 − 𝟒𝟖𝒙𝟒 + 𝟑𝟖𝟒𝒙𝟐 )𝑱𝟎

⟹ 𝑳 = (𝒙𝟔 − 𝟏𝟔𝒙𝟑 + 𝟔𝟒𝒙)𝑱𝟏 + (𝟒𝒙𝟒 − 𝟑𝟐𝒙𝟐 )𝑱𝟎

3-) Exprimons chacune des fonctions 𝑱𝟓 (𝒙) et 𝑱𝟗 (𝒙) en fonction de sin(𝑥) et

𝟐 𝟐

cos(𝑥).
2𝑝
Par définition : 𝐽𝑝+1 (𝑥) = 𝐽𝑝 (𝑥) − 𝐽𝑝−1 (𝑥). Variation de p conduit à :
𝑥
3−𝑥 2 3
𝐽5 (𝑥) = ( ) 𝐽1 (𝑥) − 𝑥 𝐽−1 (𝑥)
2 𝑥2 2 2
{ 𝑥 4 −45𝑥 2 +105 10𝑥 2 −70
𝐽9 (𝑥) = ( ) 𝐽1 (𝑥) + ( ) 𝐽−1 (𝑥)
2 𝑥4 2 𝑥3 2

1
La solution générale de l’équation de Bessel d’ordre 2 s’écrit :
Page 6 sur 6