ATLAS Ouvrir dans l'atlas →

Proposition 3Distributions et EDP

Dérivée d'une distribution

Soit T∈𝒟′(Ω)T \in {\mathcal D'}(\Omega). Sa dérivée partielle ∂jT\partial_j T est la forme ⟨∂jT,φ⟩=−⟨T,∂jφ⟩.\langle \partial_j T, \varphi\rangle = -\langle T, \partial_j\varphi\rangle. Pour un multi-indice α\alpha, on pose ∂αT=∂1α1⋯∂nαnT\partial^\alpha T = \partial_1^{\alpha_1} \cdots \partial_n^{\alpha_n}T, si bien que ⟨∂αT,φ⟩=(−1)∣α∣⟨T,∂αφ⟩\langle \partial^\alpha T, \varphi\rangle = (-1)^{|\alpha|}\langle T, \partial^\alpha\varphi\rangle. Le laplacien est ΔT=∑j∂j2T{\Delta}T = \sum_j \partial_j^2 T ; pour une fonction uu de classe C2C^2, Δu=∑j∂j2u\Delta u = \sum_j \partial_j^2 u. Sur ℝ{\mathbb R}, on écrit T′T' pour ∂1T\partial_1 T.

La fonction de Heaviside est H=1]0,+∞[H = \mathbf 1_{]0, +\infty[}, et sgn⁡\operatorname{sgn} la fonction égale à −1-1 sur ]−∞,0[]{-\infty}, 0[, à 11 sur ]0,+∞[]0, +\infty[ et à 00 en 00.

Énoncé

(i) ∂jT\partial_j T est une distribution ; ∂i∂jT=∂j∂iT\partial_i\partial_j T = \partial_j\partial_i T ; et si Tk→TT_k \to T dans 𝒟′(Ω)\mathcal D'(\Omega), alors ∂jTk→∂jT\partial_j T_k \to \partial_j T.

(ii) Compatibilité. Si ff est de classe C1C^1 sur Ω\Omega, ∂jTf=T∂jf\partial_j T_f = T_{\partial_j f}. Par suite, pour ff de classe CkC^k et ∣α∣⩽k|\alpha| \leqslant k, ∂αTf=T∂αf\partial^\alpha T_f = T_{\partial^\alpha f}, et ΔTu=TΔu\Delta T_u = T_{\Delta u} pour uu de classe C2C^2.

(iii) Formule des sauts. Soit f:ℝ→ℝf : \mathbb R \to \mathbb R dont les restrictions à ]−∞,0[]{-\infty}, 0[ et à ]0,+∞[]0, +\infty[ sont de classe C1C^1, et telle que ff et f′f' aient des limites finies f(0±)f(0^\pm) et f′(0±)f'(0^\pm) à droite et à gauche de 00. Alors ff et f′f', définie hors de 00, sont localement intégrables, et (Tf)′=Tf′+(f(0+)−f(0−)) δ.(T_f)' = T_{f'} + \big(f(0^+) - f(0^-)\big)\,\delta.

(iv) Exemples. H′=δH' = \delta ; ∣x∣′=sgn⁡|x|' = \operatorname{sgn} et ∣x∣′′=2δ|x|'' = 2\delta ; (ln⁡∣x∣)′=vp⁡1x(\ln|x|)' = \operatorname{vp}\frac1x.

Démonstration

(i) Si φk→φ\varphi_k \to \varphi dans 𝒟(Ω)\mathcal D(\Omega), alors ∂jφk→∂jφ\partial_j\varphi_k \to \partial_j\varphi (Fonctions test, (i)), et ⟨T,∂jφk⟩→⟨T,∂jφ⟩\langle T, \partial_j\varphi_k\rangle \to \langle T, \partial_j\varphi\rangle : la forme linéaire ∂jT\partial_j T est une distribution. Le théorème de Schwarz pour φ\varphi (Schwarz, (i)) donne la symétrie. Enfin ⟨∂jTk,φ⟩=−⟨Tk,∂jφ⟩→−⟨T,∂jφ⟩\langle \partial_j T_k, \varphi\rangle = -\langle T_k, \partial_j\varphi\rangle \to -\langle T, \partial_j\varphi\rangle.

(ii) La fonction ∂jf\partial_j f est continue, donc localement intégrable. L'intégration par parties (Intégrales à paramètre, (ii)) donne ⟨∂jTf,φ⟩=−∫Ωf ∂jφ=∫Ω(∂jf) φ.\langle \partial_j T_f, \varphi\rangle = -\int_\Omega f\,\partial_j\varphi = \int_\Omega (\partial_j f)\,\varphi. La suite s'obtient par récurrence sur ∣α∣|\alpha|.

(iii) Intégrabilité. ff et f′f' sont boréliennes, continues sur chaque demi-droite ouverte, et bornées sur tout borné grâce aux limites en 00 : elles sont localement intégrables.

Prolongement en 00. Soit f+f_+ la fonction égale à ff sur ]0,+∞[]0, +\infty[ et à f(0+)f(0^+) en 00. Pour x>0x > 0, l'égalité des accroissements finis (Égalité des accroissements finis) donne cx∈ ]0,x[c_x \in \,]0, x[ tel que f+(x)−f+(0)x=f′(cx)\frac{f_+(x) - f_+(0)}x = f'(c_x), qui tend vers f′(0+)f'(0^+) : f+f_+ est de classe C1C^1 sur [0,+∞[[0, +\infty[. De même à gauche.

Intégration par parties. Soit φ∈𝒟(ℝ)\varphi \in \mathcal D(\mathbb R) de support contenu dans [−R,R][-R, R]. Le théorème fondamental (Théorème fondamental) appliqué à f+φf_+\varphi sur [0,R][0, R] donne −∫0Rfφ′=f(0+) φ(0)+∫0Rf′φ,-\int_0^R f\varphi' = f(0^+)\,\varphi(0) + \int_0^R f'\varphi, et de même −∫−R0fφ′=−f(0−) φ(0)+∫−R0f′φ-\int_{-R}^0 f\varphi' = -f(0^-)\,\varphi(0) + \int_{-R}^0 f'\varphi. La somme est la formule annoncée, évaluée sur φ\varphi.

(iv) Pour HH, f′=0f' = 0 hors de 00 et le saut vaut 11. Pour ∣x∣|x|, f′=sgn⁡f' = \operatorname{sgn} hors de 00, sans saut ; pour sgn⁡\operatorname{sgn}, f′=0f' = 0 et le saut vaut 22, d'où ∣x∣′′=sgn⁡′=2δ|x|'' = \operatorname{sgn}' = 2\delta.

Le logarithme. La fonction ln⁡∣x∣\ln|x| (Puissances réelles) est localement intégrable : pour 0<ε<10 < \varepsilon < 1, ∫ε1(−ln⁡x) dx=1−ε+εln⁡ε\int_\varepsilon^1 (-\ln x)\,dx = 1 - \varepsilon + \varepsilon\ln\varepsilon, et la convergence monotone (Convergence monotone) donne ∫01∣ln⁡x∣ dx=1\int_0^1 |\ln x|\,dx = 1, l'intégrale de Riemann d'une fonction continue sur un segment étant celle de Lebesgue (Riemann et Lebesgue).

Soit φ\varphi de support dans [−R,R][-R, R], et ε=1k\varepsilon = \frac1k. La convergence dominée donne −∫ln⁡∣x∣ φ′(x) dx=lim⁡kIk-\int \ln|x|\,\varphi'(x)\,dx = \lim_k I_k, où IkI_k est l'intégrale sur {ε⩽∣x∣⩽R}\{\varepsilon \leqslant |x| \leqslant R\}. Sur chacun des deux segments, l'intégration par parties donne Ik=ln⁡ε (φ(ε)−φ(−ε))+(∫−R−ε+∫εR)φ(x)x dx.\begin{aligned} I_k &= \ln\varepsilon\,\big(\varphi(\varepsilon) - \varphi(-\varepsilon)\big) \\ &\quad + \Big(\int_{-R}^{-\varepsilon} + \int_\varepsilon^R\Big)\frac{\varphi(x)}x\,dx. \end{aligned} Le premier terme est majoré par 2ε∣ln⁡ε∣ p1(φ)2\varepsilon|\ln\varepsilon|\,p_1(\varphi). Comme ln⁡z⩽z−1\ln z \leqslant z - 1 pour z⩾1z \geqslant 1, par les accroissements finis, on a ∣ln⁡ε∣=2ln⁡ε−1/2⩽2ε−1/2|\ln\varepsilon| = 2\ln\varepsilon^{-1/2} \leqslant 2\varepsilon^{-1/2}, donc ε∣ln⁡ε∣⩽2ε→0\varepsilon|\ln\varepsilon| \leqslant 2\sqrt\varepsilon \to 0. Le second terme tend vers ⟨vp⁡1x,φ⟩\langle \operatorname{vp}\frac1x, \varphi\rangle (Exemples de distributions, (iv)).

□
Remarque

L'erreur la plus répandue est de confondre la dérivée au sens des distributions avec la dérivée presque partout. La fonction HH est dérivable hors de 00, de dérivée nulle : sa dérivée presque partout est 00. Sa dérivée au sens des distributions est δ\delta, qui n'est pas nulle : elle garde la trace du saut, que la dérivée presque partout oublie.

x x x |x| sgn 1 −1 2 2δ sgn′
En haut, la fonction ∣x∣|x| ; au milieu, sa dérivée sgn⁡\operatorname{sgn}, qui saute de −1-1 à 11 en 00 ; en bas, la dérivée de sgn⁡\operatorname{sgn}, qui est nulle hors de 00 et égale à 2δ2\delta, dessinée par une flèche rouge de hauteur 22, la hauteur du saut. Chaque dérivation au sens des distributions transforme un coin en saut, puis un saut en masse de Dirac.
Remarque

Toute distribution est indéfiniment dérivable, alors qu'une fonction continue peut n'être dérivable nulle part. C'est la force de la définition : la dérivée est transférée sur la fonction test, qui en a autant qu'on veut. Le prix est que la dérivée d'une fonction n'est plus, en général, une fonction.

Remarque

Pour φ(x)=e−x2(1+x)\varphi(x) = e^{-x^2}(1 + x), qui n'est pas à support compact mais décroît assez vite pour que les calculs aient un sens, on trouve −∫ℝln⁡∣x∣ φ′(x) dx=∫0∞φ(x)−φ(−x)x dx=2∫0∞e−x2 dx=π≈1,7725-\int_{\mathbb R}\ln|x|\,\varphi'(x)\,dx = \int_0^\infty \frac{\varphi(x) - \varphi(-x)}x\,dx = 2\int_0^\infty e^{-x^2}\,dx = \sqrt\pi \approx 1{,}7725, valeur vérifiée par un calcul numérique des deux membres. La partie paire de φ\varphi ne compte pas : la valeur principale ne voit que la partie impaire.

Ouvrir dans l'atlas interactif →