Fonction de Bessel

Chen Thomas
[Link].thomas1[at][Link]
17 mai 2024

Attention
1. Ce document contient certainement des coquilles. N’hésitez pas à me le signaler. De même si vous avez
une question.
2. Pour les recasages, ce sont les miens mais ce développement se case peut-être ailleurs et je n’y ai pas
réfléchi.
3. Il se peut que ce développement dure plus de 15 minutes. J’ai essayé de le découper pour faire des
recollements personnalisés.

Leçons
• 221 : Equations différentielles linéaires. Systèmes d’équations différentielles linéaires. Exemples et ap-
plications.
• 243 : Séries entières, propriétés de la somme. Exemples et applications.

Références
[1] C. Deschamps, F. Moulin, etc. Maths Tout-en-un MP/MP*-MPI/MPI*. Dunod, 2022.
[2] S. Francinou, H. Gianella, S. Nicolas. Exercices de mathématiques. Oraux X-ENS, Analyse 4. Cassini, 2012.

Tout est dans [2] et on peut quelque temps s’aider de [1].

Soit (E) l’équation différentielle xy ′′ + y ′ + xy = 0 sur R.

Proposition 1. Il existe une unique solution de (E) développable en série entière qui vaut 1 en 0. Toute
autre solution définie sur ]0, α[ pour α > 0 non proportionnelle à celle-ci est non bornée. On a d’ailleurs une
expression intégrale de la fonction développable en série entière.
X
Démonstration. Soit f = an xn une série entière de rayon de convergence R > 0. Soit x ∈] − R, R[.
n≥0

+∞
xf ′′ (x) =
X X
an n(n − 1)xn−1 = an+1 n(n + 1)xn .
n=2 n=1

+∞
f ′ (x) =
X
an+1 (n + 1)xn .
n=0
+∞
X
xf (x) = an−1 xn .
n=1

1
 Chen Thomas

Analyse : f est solution de (E) si, et seulement si

+∞
X
(an+1 n(n + 1) + an+1 (n + 1) + an−1 ) xn + a1 x0 = 0
n=1

ce qui revient à demander par unicité des coefficients d’une série entière

∀n ∈ N∗ , an+1 (n + 1)2 + an−1 = 0, a1 = 0.

Par récurrence, on a donc immédiatement

−a2n (−1)n+1 a0
∀n ∈ N, a2n+1 = 0 ; a2n+2 = = · · · = .
[2(n + 1)]2 22(n+1) (n + 1)!2

Si on impose f (0) = 1 à savoir a0 = 1, on a donc

+∞
X (−1)n 2n
f (x) = x .
n=0
4n (n!)2

Synthèse : soit f série entière définit ci-avant. Montrons que son rayon de convergence est strictement positif.

(−1) p
 si n = 2p
Pour cela, le lemme d’Abel nous dit que R = +∞. En effet, en notant bn = 4p (p!)2 , la
0 si n = 2p + 1

suite (bn Rn )n∈N est bornée pour tout R > 0. Il suffit d’écrire

(−1)p 2p
∀p ∈ N, n = 2p, R −→ 0.
4p (p!)2 n→+∞

Par les calculs antérieurs à l’analyse, f est solution de (E) et f (0) = 1 : c’est ce qu’il fallait démontrer.
f g
Soit g telle que (f, g) soit libre et g solution de (E). Soit W ( wonskien de f, g, c’est-à-dire, W = ′ ′ =
f g
f g ′ − f ′ g. Alors W (f, g) est dérivable de dérivée W ′ = f ′ g ′ + f g ′′ − f ′′ g − f ′ g ′ donc xW ′ = f (−g ′ − xg) −
(−f ′ − xf )g = −f g ′ + f ′ g = −W . Ainsi, (xW )′ est constant sur ]0, α[ donc il existe une constante réelle A tel
que ∀x ∈]0, α[, xW (x) = A. A est non nulle : sinon, W est nulle ce qui contredit la liberté de (f, g). g étant
bornée, je la prolonge en 0. Ainsi,

A = xf (x)g ′ (x) − xf ′ (x)g(x) −→ xf (x)g ′ (x).

x→0+
Z x0
1
Puisque f (x) −→ 1, on a donc g ′ (x) A /x. Soit x0 ∈]0, α[. Puisque dt diverge, par théorème d’inté-
x→0+ x→0+ 0 t
gration d’équivalents, on a
g(x0 ) − g(x) ∼ A(ln(x0 ) − ln(x)).
x→0

Cela contredit le caractère bornée de g. Ainsi, par l’absurde, on a g non bornée.

Z π/2
Soit J : x ∈ R 7→ cos(x sin(t))dt. Soit f : (x, t) ∈ R × [0, π/2] 7→ cos(x sin(t)) ∈ R.
0
∂f ∂2f
• ∀t ∈ [0, π/2], x 7→ f (x, t) ∈ C 2 (R) ; ∀t ∈ [0, π/2], ∀x ∈ R, (t, x) = − sin(t) sin(x sin(t)), 2 (t, x) =
∂x ∂x
− sin2 cos(x sin(t)).
• Pour tout x réel, ces dernières applications sont mesurables sur [0, π/2].
• Domination : chacune de ces applications sont continues sur le segment [0, π/2] donc dominée sur [0, π/2].
Par le théorème de dérivation C 2 sous le signe intégral, J est de classe C 2 et
Z π/2 Z π/2
J ′ (x) = − sin(t) sin(x sin(t))dt , J ′′ (x) = − sin2 cos(x sin(t))dt.
0 0

[Link].thomas1[at][Link] Page 2
 Chen Thomas

Ainsi,
Z π/2 Z π/2
′′ 2
x(J (x) + J(x)) = x (1 − sin (t)) cos(x sin(t))dt = cos(t)(x cos(t)) cos(x sin(t))dt.
0 0

On réalise une intégration par parties, on est en présence de deux fonctions C 1 : on intègre t 7→ (x cos(t)) cos(x sin(t))
dont une primitive est t 7→ sin(x sin(t)) et on dérive cos. On obtient alors
Z π/2
π/2
x(J ′′ (x) + J(x)) = [sin(x sin(t)) cos(t)]0 + sin(t) sin(x sin(t))dt = −J ′ (x).
| {z } 0
=0

On a donc J solution de (E). De plus, J est bornée sur R. Par précédent, J = λf avec λ ∈ R. Or, J(0) = 1 =
λf (0) = λ donc J = f .

[Link].thomas1[at][Link] Page 3