Transformation de Weierstrass

Weierstrass a utilisé cette transformation dans sa preuve originale du théorème d'approximation de Weierstrass. Elle est également appelée transformation de Gauss ou de Gauss–Weierstrass du nom de Carl Friedrich Gauss et aussi transformation de Hille d'après Einar Carl Hille qui l'a étudié en profondeur. La généralisation W_t mentionnée supra est connue en traitement du signal comme un filtre de Gauss et en traitement de l'image (quand définie sur ${\displaystyle \mathbb {R} ^{2}}$) comme un flou gaussien.

Toute fonction constante est également à sa transformée de Weierstrass.

La transformée de Weierstrass d'un polynôme est un polynôme de même degré. En effet, si H_n désigne le polynôme d'Hermite (physique) de degré n, alors la transformée de Weierstrass de H_n(x/2) est égale à xⁿ. On peut le voir en utilisant le fait que la fonction génératrice des polynômes d'Hermite est proche du noyau gaussien utilisée dans la définition de la transformée de Weierstrass.

La transformée de Weierstrass de la fonction e^ax (avec a une constante arbitraire) est e^a2 e^ax. Les fonctions exponentielles de la forme e^ax sont donc des fonctions propres de la transformation de Weierstrass. (Cette propriété est, en fait, vraie pour toutes les transformations basées sur une convolution.)

En posant a=bi avec i l'unité imaginaire, et grâce à l'identité d'Euler, on peut vérifier que les transformées de Weierstrass de fonction cos(bx) est e^−b2 cos(bx) et celle de sin(bx) vaut e^−b2 sin(bx).

La transformée de Weierstrass de la fonction e^ax2 est

${\displaystyle {\frac {1}{\sqrt {1-4a}}}\mathrm {e} ^{\frac {ax^{2}}{1-4a}}}$ vraie si a < 1/4 et non définie sinon.

En particulier, en choisissant a négatif, il est évident que la transformée de Weierstrass d'une gaussienne est une autre gaussienne, mais plus large.

La transformation de Weierstrass qui associe une fonction f à une fonction F est linéaire et invariante par translation : la transformée de f(x + a) est F(x + a).

Si la transformée F(x) est bien définie pour les réels x = a et x = b, alors il existe pour toutes les valeurs comprises entre elles et y forme une fonction analytique ; plus encore, F(x) existe pour tout complexe x tel que a ≤ Re(x) ≤ b et forme une fonction holomorphe sur cette bande du plan complexe. C'est la traduction précise de la propriété de "lissage" de la transformation.

Si f est intégrable sur toute la droite réelle (i.e. f ∈ L¹(ℝ)), sa transformée de Weierstrass F l'est également, de plus, si f est positive, F aussi et leurs intégrales sur ℝ sont égales. On retrouve ainsi la propriété physique que l'énergie thermique totale est conservée par l’équation de la chaleur, ou que la quantité totale de matériau diffusant est conservé par l'équation de diffusion.

Par la suite, il vient que pour 0 < p ≤ ∞ et f ∈ L^p(ℝ), alors F ∈ L^p(ℝ) et ||F||_p ≤ ||f||_p. La transformée de Weierstrass permet de définir un opérateur borné W : L^p(ℝ) → L^p(ℝ).

Si f est assez régulière, alors la transformée de Weierstrass de la dérivée d'ordre k de f est égale à la dérivée d'ordre k de la transformée de Weierstrass de f.

On peut relier la transformation de Weierstrass W et la transformée bilatérale de Laplace L. En notant

${\displaystyle g(x)=\mathrm {e} ^{-{\frac {x^{2}}{4}}}f(x)}$

alors

${\displaystyle W[f](x)={\frac {1}{\sqrt {4\pi }}}\mathrm {e} ^{-x^{2}/4}L[g]\left(-{\frac {x}{2}}\right).}$

La formule qui suit, proche de la transformée de Laplace d'une gaussienne, et un analogue réel de la transformation de Hubbard-Stratonovich (en), est assez simple à établir :

${\displaystyle \mathrm {e} ^{u^{2}}={\frac {1}{\sqrt {4\pi }}}\int _{-\infty }^{\infty }\mathrm {e} ^{-uy}\mathrm {e} ^{-y^{2}/4}\,\mathrm {d} y.}$

En remplaçant u par l'opérateur formel de dérivation D=^d⁄_dx et en utilisant l'opérateur de décalage de Lagrange

${\displaystyle \mathrm {e} ^{-yD}f(x)=f(x-y)}$,

(cette égalité se déduit de la série de Taylor de l'exponentielle), on a

${\displaystyle {\begin{aligned}\mathrm {e} ^{D^{2}}f(x)&={\frac {1}{\sqrt {4\pi }}}\int _{-\infty }^{\infty }\mathrm {e} ^{-yD}f(x)\mathrm {e} ^{-y^{2}/4}\,\mathrm {d} y\\&={\frac {1}{\sqrt {4\pi }}}\int _{-\infty }^{\infty }f(x-y)\mathrm {e} ^{-y^{2}/4}\,\mathrm {d} y=W[f](x)\end{aligned}}}$

ce qui permet de déduire une expression formelle de la transformation de Weierstrass W,

${\displaystyle W=\mathrm {e} ^{D^{2}}}$

où l'opérateur est à comprendre comme agissant sur f (x) comme

${\displaystyle \mathrm {e} ^{D^{2}}f(x)=\sum _{k=0}^{\infty }{\frac {D^{2k}f(x)}{k!}}~.}$

L'égalité précédente dépend de propriétés de convergence non précisées, et l'égalité W=e^D2 n'est donc pas toujours vérifiée ; ainsi, il existe des fonctions f dont la transformée de Weierstrass est bien définie, mais pour lesquelles e^D2f(x) n'a aucun sens.

Elle est cependant utile et permet de trouver rapidement les transformées de Weierstrass de fonctions usuelles simples, comme les polynômes, les fonctions exponentielles et trigonométriques.

L'opérateur inverse formel de la transformation de Weierstrass est donc donné par

${\displaystyle W^{-1}=\mathrm {e} ^{-D^{2}}~.}$

Une fois encore, la formule est vérifiée sous réserve de propriétés de convergence et ne peut que servir de guide. Elle est vraie pour certaines classes de fonctions si l'opérateur donné dans le membre droit est bien défini[2].

On peut, par un autre moyen, essayer d'inverser la transformée de Weierstrass : d'une fonction analytique donnée

${\displaystyle F(x)=\sum _{n=0}^{\infty }a_{n}x^{n}~,}$

on applique W⁻¹ pour obtenir

${\displaystyle f(x)=W^{-1}[F(x)]=\sum _{n=0}^{\infty }a_{n}W^{-1}[x^{n}]=\sum _{n=0}^{\infty }a_{n}H_{n}(x/2)}$

avec les polynômes d'Hermite (physiques) H_n.

Cette formule de f(x) reste un indicateur formel, vraie si la série à droite converge. Mais si, par exemple, f ∈ L²(ℝ), alors il suffit de connaître les valeurs des dérivées de F en x = 0 pour obtenir les coefficients a_n, et donc reconstruire f comme une série de polynômes d'Hermite.

Une troisième méthode d'inversion consiste à utiliser l'égalité faisant intervenir la transformation de Laplace donnée supra, dont on connait la formule d'inversion. Le résultat est donné dans la suite pour les distributions.

La transformation de Weierstrass généralisée permet de donner une forme exacte de la solution de l'équation de la chaleur (ici en 2D).

Le noyau de convolution peut être remplacé par le noyau gaussien ${\displaystyle {\frac {1}{\sqrt {4\pi t}}}\mathrm {e} ^{-{\frac {x^{2}}{4t}}}}$ (pour t > 0) au lieu de ${\displaystyle {\frac {1}{\sqrt {4\pi }}}\mathrm {e} ^{-{\frac {x^{2}}{4}}}}$, qui permet de définir l'opérateur W_t de la transformation généralisée de Weierstrass.

Pour les petites valeurs de t, W_t[f] est très proche de f, mais régulière. Plus t est grand, plus l'opérateur va moyenner et changer les variations de f. Physiquement, W_t revient à suivre l'évolution de l'équation de la chaleur (ou de diffusion) avec t correspondant au temps, avec une propriété additive,

${\displaystyle W_{s}\circ W_{t}=W_{s+t},}$

qui se traduit par « diffuser pendant une durée t, puis pendant une durée s, est équivalent à diffuser pendant une durée s + t ». On peut étendre l'idée en définissant à t = 0 l'opérateur W₀ qui est l'opérateur identité (le noyau devenant alors une distribution de Dirac), et construisant ainsi un semi-groupe d'opérateurs.

Le noyau ${\displaystyle {\frac {1}{\sqrt {4\pi t}}}\mathrm {e} ^{-{\frac {x^{2}}{4t}}}}$ de la transformation généralisée de Weierstrass est parfois appelé noyau de Gauss–Weierstrass, et est la fonction de Green de l'équation de diffusion sur ℝ.

W_t peut être écrit à partir de W : pour une fonction f(x), on définit une nouvelle fonction f_t(x) = f(x√t), puis on pose W_t[f](x) = W[f_t](x/√t), une conséquence de la règle de substitution.

La transformation de Weierstrass peut également être définie pour certaines classes de distributions ou "fonctions généralisées"[3]. Par exemple, la transformée de Weierstrass d'une fonction delta de Dirac est la gaussienne ${\displaystyle {\frac {1}{\sqrt {4\pi }}}\mathrm {e} ^{-x^{2}/4}}$.

Dans ce contexte, des formules d'inversion plus rigoureuses peuvent être prouvées :

${\displaystyle f(x)=\lim _{r\to \infty }{\frac {1}{i{\sqrt {4\pi }}}}\int _{x_{0}-ir}^{x_{0}+ir}F(z)\mathrm {e} ^{\frac {(x-z)^{2}}{4}}\,\mathrm {d} z}$

où x₀ est un réel donnée pour lequel F(x₀) existe, l'intégrale s'étend autour de la ligne verticale de plan complexe de partie réelle x₀, et cette limite est à voir au sens des distributions.

Partant de là, la transformation de Weierstrass peut être définie pour des fonctions à valeurs réelles ou complexes définies sur ℝⁿ. Il suffit de poser la convolution sur ℝⁿ et de considérer le terme (x − y)² comme le carré de la norme euclidienne du vecteur x − y; le facteur constant devant l’intégrale est également à ajuster pour retrouver une gaussienne d'intégrale égale à 1.

Plus généralement, la transformation de Weierstrass peut être définie sur une variété riemannienne : l'équation de la chaleur peut y être formulée (par l'opérateur de Laplace-Beltrami sur la variété), et la transformée de Weierstrass W[f] est alors donnée en suivant la solution de l'équation de la chaleur dans le temps, avec comme solution initiale une "distribution de température" f.

En considérant la convolution avec le noyau 1/(π(1 + x²)) au lieu d'une gaussienne, on trouve la transformation de Poisson qui a des propriétés de régularisation et moyennage similaires.

• Flou gaussien
• Filtre de Gauss
• Q-représentation Husimi

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Weierstrass transform » (voir la liste des auteurs).