ATLAS Ouvrir dans l'atlas →

Définition 1Distributions et EDP

Fonctions test et noyau régularisant

Soit Ω\Omega un ouvert non vide de ℝn{\mathbb R^n}. On note ∣x∣=∥x∥2|x| = \|x\|_2 la norme euclidienne (Normes), B(a,r)B(a, r) et Bˉ(a,r)\bar B(a, r) les boules ouverte et fermée qu'elle définit, et λ\lambda la mesure de Lebesgue de ℝn\mathbb R^n (Lebesgue sur ℝⁿ). On écrit ∫Af\int_A f pour ∫Af dλ\int_A f\,d\lambda.

Un multi-indice est un α=(α1,…,αn)∈ℕn\alpha = (\alpha_1, \dots, \alpha_n) \in {\mathbb N^n}, de longueur ∣α∣=α1+⋯+αn|\alpha| = \alpha_1 + \cdots + \alpha_n. Pour ff de classe C∞C^\infty sur Ω\Omega, on pose ∂αf=∂1α1⋯∂nαnf\partial^\alpha f = \partial_1^{\alpha_1} \cdots \partial_n^{\alpha_n} f ; l'ordre des dérivations est indifférent, par le théorème de Schwarz (Schwarz, (i)). On écrit β⩽α\beta \leqslant \alpha si βi⩽αi\beta_i \leqslant \alpha_i pour tout ii, et eje_j est le multi-indice dont seule la jj-ième composante vaut 11.

Le support d'une fonction continue f:Ω→ℝf : \Omega \to {\mathbb R}, noté supp⁡f\operatorname{supp} f, est l'adhérence dans Ω\Omega de {f≠0}\{f \neq 0\}. Une fonction test sur Ω\Omega est une fonction de classe C∞C^\infty sur Ω\Omega dont le support est compact ; elles forment l'ensemble 𝒟(Ω){\mathcal D}(\Omega).

Pour un compact K⊂ΩK \subset \Omega, on note 𝒟K(Ω)\mathcal D_K(\Omega) l'ensemble des φ∈𝒟(Ω)\varphi \in \mathcal D(\Omega) de support contenu dans KK, et, pour m∈ℕm \in \mathbb N, pm(φ)=max⁡∣α∣⩽m sup⁡x∈Ω∣∂αφ(x)∣,p_m(\varphi) = \max_{|\alpha| \leqslant m}\ \sup_{x \in \Omega} |\partial^\alpha \varphi(x)|, qui est fini, puisque chaque ∂αφ\partial^\alpha\varphi est continue et nulle hors du compact supp⁡φ\operatorname{supp} \varphi (Compacts de ℝ, (iii)).

Une suite (φk)(\varphi_k) de 𝒟(Ω)\mathcal D(\Omega) converge vers φ\varphi dans 𝒟(Ω)\mathcal D(\Omega) s'il existe un compact K⊂ΩK \subset \Omega tel que supp⁡φk⊂K\operatorname{supp}\varphi_k \subset K pour tout kk, et si pm(φk−φ)→0p_m(\varphi_k - \varphi) \to 0 pour tout mm : les supports restent dans un même compact, et chaque dérivée converge uniformément.

Pour ε>0\varepsilon > 0, on pose Ωε={x∈Ω:Bˉ(x,ε)⊂Ω}\Omega_\varepsilon = \{x \in \Omega : \bar B(x, \varepsilon) \subset \Omega\} ; si Ω=ℝn\Omega = \mathbb R^n, Ωε=ℝn\Omega_\varepsilon = \mathbb R^n.

Soit enfin θ\theta la fonction du nœud Fonctions plateau, nulle sur ]−∞,0]]{-\infty}, 0] et strictement positive sur ]0,+∞[]0, +\infty[. On pose σ(t)=θ(t)θ(t)+θ(1−t),ρ~(x)=θ(1−∣x∣2),cn=∫ℝnρ~,\begin{aligned} \sigma(t) &= \frac{\theta(t)}{\theta(t) + \theta(1 - t)}, \\ \tilde\rho(x) &= \theta(1 - |x|^2), \qquad c_n = \int_{\mathbb R^n} \tilde\rho, \end{aligned} puis ρ=ρ~/cn\rho = \tilde\rho / c_n et, pour ε>0\varepsilon > 0, ρε(x)=ε−nρ(x/ε)\rho_\varepsilon(x) = \varepsilon^{-n}\rho(x/\varepsilon). La famille (ρε)(\rho_\varepsilon) est le noyau régularisant.

Énoncé

(i) Opérations. 𝒟(Ω)\mathcal D(\Omega) est un espace vectoriel, stable par chaque ∂α\partial^\alpha et par multiplication par toute fonction aa de classe C∞C^\infty sur Ω\Omega. Plus précisément, ∂α(aφ)\partial^\alpha(a\varphi) est une somme de termes c ∂βa ∂α−βφc\,\partial^\beta a\,\partial^{\alpha - \beta}\varphi, où β⩽α\beta \leqslant \alpha et où les entiers c⩾0c \geqslant 0 ne dépendent que de α\alpha et β\beta.

Si φk→φ\varphi_k \to \varphi dans 𝒟(Ω)\mathcal D(\Omega), alors ∂αφk→∂αφ\partial^\alpha\varphi_k \to \partial^\alpha\varphi et aφk→aφa\varphi_k \to a\varphi dans 𝒟(Ω)\mathcal D(\Omega). Si ω⊂Ω\omega \subset \Omega est ouvert, toute φ∈𝒟(ω)\varphi \in \mathcal D(\omega), prolongée par 00, appartient à 𝒟(Ω)\mathcal D(\Omega) et y a le même support.

(ii) Noyau régularisant. cnc_n est un réel strictement positif, et ρ\rho est de classe C∞C^\infty sur ℝn\mathbb R^n, positive, paire, de support Bˉ(0,1)\bar B(0, 1) et d'intégrale 11. Pour ε>0\varepsilon > 0, ρε\rho_\varepsilon a les mêmes propriétés, avec le support Bˉ(0,ε)\bar B(0, \varepsilon), et ρε(x)=hε(∣x∣)\rho_\varepsilon(x) = h_\varepsilon(|x|), où hε:ℝ→ℝh_\varepsilon : \mathbb R \to \mathbb R est de classe C∞C^\infty, décroissante sur [0,+∞[[0, +\infty[ et nulle sur [ε,+∞[[\varepsilon, +\infty[.

(iii) Troncatures. σ\sigma est de classe C∞C^\infty sur ℝ\mathbb R, croissante, à valeurs dans [0,1][0, 1], nulle sur ]−∞,0]]{-\infty}, 0] et égale à 11 sur [1,+∞[[1, +\infty[. Si KK est un compact contenu dans un ouvert OO de ℝn\mathbb R^n, il existe χ∈𝒟(O)\chi \in \mathcal D(O), à valeurs dans [0,1][0, 1], égale à 11 sur un ouvert contenant KK.

(iv) Voisinages des compacts. Pour ε>0\varepsilon > 0, Ωε\Omega_\varepsilon est ouvert. Tout compact K⊂ΩK \subset \Omega est contenu dans Ωε\Omega_\varepsilon pour ε\varepsilon assez petit. Si K⊂ΩεK \subset \Omega_\varepsilon est compact, l'ensemble K+Bˉ(0,ε)K + \bar B(0, \varepsilon) des points k+vk + v, où k∈Kk \in K et ∣v∣⩽ε|v| \leqslant \varepsilon, est un compact contenu dans Ω\Omega.

Démonstration

(i) Stabilité. Le support d'une somme est contenu dans la réunion des supports, et ∂αφ\partial^\alpha\varphi, comme aφa\varphi, est nulle sur l'ouvert Ω∖supp⁡φ\Omega \setminus \operatorname{supp}\varphi : leurs supports sont des fermés de Ω\Omega contenus dans le compact supp⁡φ\operatorname{supp}\varphi, donc compacts (Parties compactes, (i)). Les sommes, les dérivées partielles et les produits de fonctions de classe C∞C^\infty le sont (Classe Cᵏ, (ii) et (iii)).

Produits. Pour ff et gg de classe C1C^1, ∂j(fg)=(∂jf) g+f ∂jg\partial_j(fg) = (\partial_j f)\,g + f\,\partial_j g, par la dérivée d'un produit en la variable xjx_j (Opérations, (ii)). Si ∂α(aφ)\partial^\alpha(a\varphi) est une somme de termes c ∂βa ∂α−βφc\,\partial^\beta a\,\partial^{\alpha - \beta}\varphi, l'application de ∂j\partial_j à chacun donne c ∂β+eja ∂α−βφ+c ∂βa ∂α+ej−βφ,c\,\partial^{\beta + e_j}a\,\partial^{\alpha - \beta}\varphi + c\,\partial^\beta a\,\partial^{\alpha + e_j - \beta}\varphi, deux termes de la forme voulue pour α+ej\alpha + e_j : la récurrence sur ∣α∣|\alpha| conclut.

Suites. Si supp⁡φk⊂K\operatorname{supp}\varphi_k \subset K, la fonction φ\varphi, limite uniforme des φk\varphi_k, est nulle hors de KK, et supp⁡φ⊂K\operatorname{supp}\varphi \subset K, puisque KK est fermé. Les supports de ∂αφk\partial^\alpha\varphi_k et de aφka\varphi_k sont aussi dans KK. D'une part pm(∂αφk−∂αφ)⩽pm+∣α∣(φk−φ)p_m(\partial^\alpha\varphi_k - \partial^\alpha\varphi) \leqslant p_{m + |\alpha|}(\varphi_k - \varphi). D'autre part, si NmN_m est la somme des entiers cc qui interviennent pour les ∣α∣⩽m|\alpha| \leqslant m et Am=max⁡∣β∣⩽msup⁡K∣∂βa∣A_m = \max_{|\beta| \leqslant m} \sup_K |\partial^\beta a|, fini par compacité, pm(aφk−aφ)⩽NmAm pm(φk−φ).p_m(a\varphi_k - a\varphi) \leqslant N_m A_m\, p_m(\varphi_k - \varphi).

Prolongement. La fonction prolongée coïncide avec φ\varphi sur l'ouvert ω\omega et est nulle sur Ω∖supp⁡φ\Omega \setminus \operatorname{supp}\varphi, ouvert puisque le compact supp⁡φ\operatorname{supp}\varphi est fermé dans ℝn\mathbb R^n (Parties compactes, (ii)). Ces deux ouverts recouvrent Ω\Omega, et la classe C∞C^\infty est locale (Partition de l'unité, (i)). Le support, fermé de ℝn\mathbb R^n, ne change pas.

(ii) Régularité et signe. x↦1−∣x∣2x \mapsto 1 - |x|^2 est polynomiale, et θ\theta est de classe C∞C^\infty (Fonctions plateau, (i)) : ρ~\tilde\rho l'est par composition. Elle est strictement positive exactement sur B(0,1)B(0, 1), et paire. Tout point xx tel que ∣x∣=1|x| = 1 est limite des points (1−1k)x(1 - \frac1k)x de B(0,1)B(0, 1), et Bˉ(0,1)\bar B(0, 1) est fermée : le support est Bˉ(0,1)\bar B(0, 1).

La constante. Sur ]0,+∞[]0, +\infty[, θ′(x)=θ(x)/x2⩾0\theta'(x) = \theta(x)/x^2 \geqslant 0 (Fonctions plateau, (i), avec P1(u)=u2P_1(u) = u^2), et θ′=0\theta' = 0 sur ]−∞,0[]{-\infty}, 0[ : θ\theta est croissante (Signe de la dérivée, (i)). Donc 0⩽ρ~⩽θ(1)0 \leqslant \tilde\rho \leqslant \theta(1), et ρ~\tilde\rho est nulle hors de [−1,1]n[-1, 1]^n, d'où cn⩽2nθ(1)c_n \leqslant 2^n\theta(1).

Comme ∣x∣⩽n∥x∥∞|x| \leqslant n\|x\|_\infty (Espaces métriques, (ii)), le cube {∥x∥∞<12n}\{\|x\|_\infty < \frac1{2n}\} est contenu dans B(0,12)B(0, \frac12), où ρ~⩾θ(34)>0\tilde\rho \geqslant \theta(\frac34) > 0. Ce cube a pour mesure n−nn^{-n} (Lebesgue sur ℝⁿ, (iv)), donc cn⩾θ(34) n−n>0c_n \geqslant \theta(\frac34)\,n^{-n} > 0.

Dilatées. ρε\rho_\varepsilon est de classe C∞C^\infty, positive, paire, de support Bˉ(0,ε)\bar B(0, \varepsilon). La formule du changement de variables (Changement de variables, (ii)), pour le difféomorphisme y↦εyy \mapsto \varepsilon y de ℝn\mathbb R^n, de jacobien εn\varepsilon^n, donne ∫ℝnρε(x) dx=∫ℝnρε(εy) εn dy=∫ℝnρ=1.\begin{aligned} \int_{\mathbb R^n} \rho_\varepsilon(x)\,dx &= \int_{\mathbb R^n} \rho_\varepsilon(\varepsilon y)\,\varepsilon^n\,dy \\ &= \int_{\mathbb R^n} \rho = 1. \end{aligned}

Enfin hε(s)=ε−ncn−1 θ(1−s2/ε2)h_\varepsilon(s) = \varepsilon^{-n}c_n^{-1}\,\theta\big(1 - s^2/\varepsilon^2\big) est de classe C∞C^\infty, décroissante sur [0,+∞[[0, +\infty[ comme composée de la fonction croissante θ\theta et d'une fonction décroissante, et nulle pour s⩾εs \geqslant \varepsilon.

(iii) La fonction σ\sigma. Le dénominateur est strictement positif : θ(1−t)>0\theta(1 - t) > 0 si t<1t < 1, et θ(t)>0\theta(t) > 0 si t>0t > 0. Comme t↦1/tt \mapsto 1/t est de classe C∞C^\infty sur ]0,+∞[]0, +\infty[ (preuve de Partition de l'unité, (iii)), σ\sigma est de classe C∞C^\infty. Elle est nulle pour t⩽0t \leqslant 0, égale à 11 pour t⩾1t \geqslant 1, puisque θ(1−t)=0\theta(1 - t) = 0, et σ′(t)=θ′(t) θ(1−t)+θ(t) θ′(1−t)(θ(t)+θ(1−t))2,\sigma'(t) = \frac{\theta'(t)\,\theta(1 - t) + \theta(t)\,\theta'(1 - t)}{(\theta(t) + \theta(1 - t))^2}, qui est positive, puisque θ\theta et θ′\theta' le sont.

Une fonction bosse. La fonction τ(u)=σ(4−4u) σ(4+4u)\tau(u) = \sigma(4 - 4u)\,\sigma(4 + 4u) est de classe C∞C^\infty, à valeurs dans [0,1][0, 1], égale à 11 si ∣u∣⩽34|u| \leqslant \frac34 et nulle si ∣u∣⩾1|u| \geqslant 1.

Troncature. Si KK est vide, χ=0\chi = 0 convient. Sinon, soit δ(x)=min⁡(1,inf⁡a∉O∥x−a∥∞)\delta(x) = \min\big(1, \inf_{a \notin O}\|x - a\|_\infty\big), et δ=1\delta = 1 si O=ℝnO = \mathbb R^n. Elle est 11-lipschitzienne pour ∥⋅∥∞\|\cdot\|_\infty, et B∞(x,δ(x))⊂OB_\infty(x, \delta(x)) \subset O pour x∈Ox \in O : la preuve de Recouvrement par des cubes, (i), vaut pour tout ouvert. Elle atteint sur le compact KK un minimum δ0>0\delta_0 > 0 (Compacts de ℝ, (iii)) ; posons η=δ0/2\eta = \delta_0/2.

Pour m∈ℤnm \in {\mathbb Z^n}, la fonction gm(x)=∏lτ((xl−ηml)/η)g_m(x) = \prod_l \tau\big((x_l - \eta m_l)/\eta\big) est de classe C∞C^\infty, à valeurs dans [0,1][0, 1], égale à 11 sur le cube {∥x−ηm∥∞⩽34η}\{\|x - \eta m\|_\infty \leqslant \frac34\eta\} et nulle hors de {∥x−ηm∥∞<η}\{\|x - \eta m\|_\infty < \eta\}. L'ensemble MM des mm tels que le cube {∥x−ηm∥∞⩽η2}\{\|x - \eta m\|_\infty \leqslant \frac\eta2\} rencontre KK est fini, KK étant borné. On pose χ=1−∏m∈M(1−gm),\chi = 1 - \prod_{m \in M}(1 - g_m), de classe C∞C^\infty et à valeurs dans [0,1][0, 1].

Pour x∈Kx \in K, soit mlm_l le plus grand entier tel que ml⩽xlη+12m_l \leqslant \frac{x_l}\eta + \frac12 : alors ∥x−ηm∥∞⩽η2\|x - \eta m\|_\infty \leqslant \frac\eta2, donc m∈Mm \in M, et xx appartient à l'ouvert {∥y−ηm∥∞<34η}\{\|y - \eta m\|_\infty < \frac34\eta\}, où gm=1g_m = 1 et χ=1\chi = 1. La réunion de ces ouverts, pour m∈Mm \in M, contient KK.

Hors du compact F=⋃m∈M{∥y−ηm∥∞⩽η}F = \bigcup_{m \in M}\{\|y - \eta m\|_\infty \leqslant \eta\}, tous les gmg_m sont nuls, et χ=0\chi = 0. Si y∈Fy \in F, il existe m∈Mm \in M et x∈Kx \in K tels que ∥y−ηm∥∞⩽η\|y - \eta m\|_\infty \leqslant \eta et ∥x−ηm∥∞⩽η2\|x - \eta m\|_\infty \leqslant \frac\eta2 ; alors ∥y−x∥∞⩽32η<δ0⩽δ(x)\|y - x\|_\infty \leqslant \frac32\eta < \delta_0 \leqslant \delta(x), et y∈Oy \in O. Ainsi supp⁡χ⊂F⊂O\operatorname{supp}\chi \subset F \subset O, et χ∈𝒟(O)\chi \in \mathcal D(O).

(iv) Compacts. Soit K⊂ΩK \subset \Omega compact non vide, et δ0>0\delta_0 > 0 le minimum sur KK de la fonction δ\delta de (iii), pour O=ΩO = \Omega. Si ε<δ0\varepsilon < \delta_0, x∈Kx \in K et ∣w−x∣⩽ε|w - x| \leqslant \varepsilon, alors ∥w−x∥∞⩽∣w−x∣<δ(x)\|w - x\|_\infty \leqslant |w - x| < \delta(x) (Espaces métriques, (ii)), et w∈Ωw \in \Omega : K⊂ΩεK \subset \Omega_\varepsilon.

Ouverture. Soit x∈Ωεx \in \Omega_\varepsilon. La boule K=Bˉ(x,ε)K = \bar B(x, \varepsilon) est fermée et bornée, donc compacte (L'espace ℝⁿ, (iii)), et contenue dans Ω\Omega ; soit δ0\delta_0 comme ci-dessus. Si ∣x′−x∣<δ0|x' - x| < \delta_0 et ∣w−x′∣⩽ε|w - x'| \leqslant \varepsilon, le point z=x+(w−x′)z = x + (w - x') est dans KK, et ∥w−z∥∞⩽∣x′−x∣<δ(z)\|w - z\|_\infty \leqslant |x' - x| < \delta(z) : w∈Ωw \in \Omega. Donc B(x,δ0)⊂ΩεB(x, \delta_0) \subset \Omega_\varepsilon.

Somme. K×Bˉ(0,ε)K \times \bar B(0, \varepsilon) est fermé et borné dans ℝ2n\mathbb R^{2n}, donc compact, et son image par l'addition, linéaire donc continue, est compacte (Parties compactes, (iii)). Elle est contenue dans Ω\Omega par définition de Ωε\Omega_\varepsilon.

□
Remarque

Pour n=1n = 1, un calcul numérique donne c1=∫−11θ(1−x2) dx≈0,444c_1 = \int_{-1}^1 \theta(1 - x^2)\,dx \approx 0{,}444, d'où ρ(0)=θ(1)/c1≈0,829\rho(0) = \theta(1)/c_1 \approx 0{,}829. La fonction ρε\rho_\varepsilon, de même aire 11 sur un intervalle de longueur 2ε2\varepsilon, a pour maximum ρ(0)/ε\rho(0)/\varepsilon : elle se resserre et s'élève à mesure que ε\varepsilon décroît, et c'est ce qui lui permettra d'approcher la masse de Dirac au nœud Exemples de distributions.

Remarque

L'erreur la plus répandue est d'oublier la condition sur les supports dans la convergence de 𝒟(Ω)\mathcal D(\Omega). Sur ℝ\mathbb R, les fonctions φk(x)=1k ρ(x−k)\varphi_k(x) = \frac1k\,\rho(x - k) vérifient pm(φk)=pm(ρ)/k→0p_m(\varphi_k) = p_m(\rho)/k \to 0 pour tout mm, mais leurs supports [k−1,k+1][k - 1, k + 1] s'enfuient : elles ne convergent pas dans 𝒟(ℝ)\mathcal D(\mathbb R).

Ce n'est pas un détail. La forme linéaire φ↦∑j⩾1j φ(j)\varphi \mapsto \sum_{j \geqslant 1} j\,\varphi(j), somme finie pour chaque fonction test, prend la valeur ρ(0)≈0,829\rho(0) \approx 0{,}829 sur chaque φk\varphi_k ; avec la seule convergence uniforme des dérivées, elle ne serait pas continue, alors qu'elle sera une distribution au sens du nœud Distributions.

Remarque

La fonction ρ\rho est une fonction plateau au sens du nœud Fonctions plateau, mais construite sur la norme euclidienne plutôt que coordonnée par coordonnée. Ce caractère radial ne sert à rien pour la régularisation elle-même ; il sera décisif au nœud Régularité des fonctions harmoniques, où une moyenne pondérée par ρε\rho_\varepsilon se décompose en moyennes sur des boules.

Ouvrir dans l'atlas interactif →