Les cas où on peut calculer exactement la transformation de Fourier d'une fonction ne sont pas fréquents, quand on fait le tour des ressources disponibles sur le net on se rend compte qu'on retombe toujours sur les mêmes exemples : fonction porte, Gaussienne, Lorentzienne ... alors pour sortir un peu de l'ordinaire voici le calcul de la transformée de Fourier de la fonction ${\rm sech}(x)= {1\over \cosh(x)}$ :
$$ \widehat{\rm sech}(\xi)=\int_{-\infty}^{+\infty} e^{-ix\xi} {1\over \cosh(x)}dx=\pi \,{\rm sech}\left({\pi\xi\over 2}\right)~~~~(1)$$Pour commencer il faut remarquer que la fonction ${\rm sech}(x)$ est $C^\infty({\mathbb R})$ et à décroissance rapide puisque :
$${\rm sech}(x)\sim_{\pm\infty} {e^{-\vert x\vert }\over 2}$$
sa transformée de Fourier est donc bien définie au sens des fonctions de la classe de Schwarz ${\mathcal S}({\mathbb R})$ avec l'intégrale $(1)$. La manière la plus simple de calculer cette transformée de Fourier repose sur l'utilisation du théorème des résidu, ce qui explique peut être que ce résultat ne soit pas très populaire dans les cours d'analyse de Fourier qui sont souvent séparés des cours d'analyse complexe! Petit rappel sur le théorème des résidus :
Théorème des résidus: Soient $U$ un ouvert simplement connexe de $\mathbb C$ , $F$ une fonction holomorphe sur $U\setminus\{z_1,\dots, z_n\}$ et $\Omega\subset U$ un ouvert tel que $\gamma=\partial \Omega$ soit un lacet simple $C^1$ par morceau ne passant par aucun des points $z_k$ alors
$$\oint_\gamma F(z) dz = 2i\pi\sum_{z_k\in\Omega} {\rm res}(F,z_k)$$
où ${\rm res}(F,z_k)=a_{-1}$ dans le développement de Laurent $F(z)=\sum_{p\in{\mathbb Z}}a_p (z-z_k)^p$
$$\oint_\gamma F(z) dz = 2i\pi\sum_{z_k\in\Omega} {\rm res}(F,z_k)$$
où ${\rm res}(F,z_k)=a_{-1}$ dans le développement de Laurent $F(z)=\sum_{p\in{\mathbb Z}}a_p (z-z_k)^p$
Ici on va donc poser
$$ F(z)=e^{-iz\xi} {1\over \cosh(z)}=e^{-iz\xi} {2\over e^z + e^{-z}}=e^{-iz\xi} {2e^z\over e^{2z}+1}$$
cette fonction est clairement analytique sauf aux points où son dénominateur s'annule :
$$e^{2z}+1=0\Leftrightarrow e^{2z}=-1=e^{i\pi}\Leftrightarrow {z}=(2k+1)i{\pi\over 2},~~ k\in{\mathbb Z}$$
on va intégrer $F(z)$ sur un contour encerclant le pôle $z_0=i{\pi\over 2}$ composé de quatre segments
$$\gamma=[-R,R]\cup [R,R+i\pi]\cup [-R+i\pi,R+i\pi]\cup [-R+i\pi,-R]=\gamma_1\cup\gamma_2\cup\gamma_3\cup\gamma_4$$
et orienté dans le sens direct.
En utilisant le DSE usuel de $e^x$ avec $x=2z$ et en $z_0=i{\pi\over 2}$
$$e^{2z}=\sum_{p\in{\mathbb Z}}2^pe^{2z_0}(z-z_0)^p=\sum_{p\in{\mathbb Z}}-2^p(z-z_0)^p
\Rightarrow e^{2z}+1=-2(z-z_0)+\dots$$
on en déduit le développement en série de Laurent de $F$ (ou au moins les premiers termes) :
$$F(z)=e^{-iz\xi} {2e^z\over e^{2z}+1}= e^{-iz_0\xi} {2e^{z_0}\over -2(z-z_0)}+\sum_{p\geq 0}a_p(z-z_0)^p$$
ce qui donne le résidu de $F$ en $z_0$ :
$${\rm res}(F,z_0)=e^{-iz_0\xi} {2e^{z_0}\over -2}=e^{\pi\xi\over 2} {2e^{i\pi\over 2}\over -2}=-ie^{\pi\xi\over 2}$$
on peut maintenant appliquer le théorème des résidus le long de $\gamma$ pour obtenir que :
$$\oint_\gamma F(z) dz = 2\pi e^{\pi\xi\over 2}$$
il reste à calculer $\oint_\gamma$ et à faire tendre $R\to\infty$ en décomposant le calcul intervalle par intervalle :
sur $\gamma_1$ on a tout simplement :
$$\int_{\gamma_1} F(z) dz =\int_{-R}^{+R} e^{-ix\xi} {1\over \cosh(x)}dx\mathop{\longrightarrow}_{R\to\infty}\widehat{\rm sech}(\xi)$$
sur $\gamma_3$ on a $z=x+i\pi$ donc :
$$\cosh(x+i\pi)={e^{x+i\pi}+e^{-x-i\pi}\over 2}={e^{x} e^{i\pi}+e^{-x}e^{-i\pi}\over 2}=-\cosh(x)$$
comme pour $\gamma_1 $ on trouve
$$\int_{\gamma_3} F(z) dz =\int_{R}^{-R} e^{-ix\xi+\pi\xi} {1\over -\cosh(x)}dx\mathop{\longrightarrow}_{R\to\infty} e^{\pi\xi}\widehat{\rm sech}(\xi)$$
sur les segments $\gamma_2$ et $\gamma_4$ on a $z=\pm R+it$ donc :
$$\vert F(\pm R+it)\vert = e^{t\xi} {2\over \vert e^{R\pm it}+e^{-R\mp it}\vert}\leq 2{e^{t\xi}\over e^{R}}\mathop{\longrightarrow}_{R\to\infty}0$$
on en déduit que pour $j=2$ et $4$ que
$$\left\vert \int_{\gamma_j}F(z) dz\right\vert\leq \int_0^\pi 2{e^{t\xi}\over e^{R}}dt\leq 2\pi {e^{\pi\xi}\over e^{R}}\mathop{\longrightarrow}_{R\to\infty}0$$
au final en faisant tendre $R\to\infty$ on obtient :
$$\oint_\gamma F(z) dz\mathop{\longrightarrow}_{R\to\infty} \left(1+e^{\pi\xi}\right)\widehat{\rm sech}(\xi)+0 = 2\pi e^{\pi\xi\over 2}$$
il ne reste plus qu'à simplifier :
$$\widehat{\rm sech}(\xi) = {2\pi e^{\pi\xi\over 2}\over 1+e^{\pi\xi}}= {2\pi \over e^{\pi\xi\over 2}+e^{-{\pi\xi\over 2}}}=\pi\, {\rm sech}\left({\pi\xi\over 2}\right)$$
On remarquera aussi que en utilisant la formule de changement d'échelle
$${\mathcal F}(f(ax))(\xi)={1\over a}\widehat{f}\left({\xi\over a}\right)$$
avec $a=\sqrt{\pi\over 2}$ on obtient que la fonction $f(x)={\rm sech}\left(\sqrt{\pi\over 2}x\right)$ vérifie $\widehat{f}(\xi)=\sqrt{2\pi}f(\xi)$ comme la fonction Gaussienne $\exp(-x^2/2)$ ...
Merci, très intéressant!!!
RépondreSupprimer