mardi 23 avril 2019

Une série de Fourier exotique

Les exemples de séries de Fourier qui sont utilisés sont presque toujours des fonctions $C^1$ par morceaux dont les coefficients sont des fractions rationnelles en $n$ vérifiant une asymptotique $c_n\sim_\infty {1\over n^a}$ avec $a\in{\mathbb N}$ . Le plus souvent on peut vérifier que $a=1$ pour des fonctions discontinue  et $a\geq 2$ pour des fonctions continues, pourtant on peut avoir des séries de Fourier avec des comportements plus exotiques correspondant à des fonctions continues mais pas $C^1$ par morceaux. Je suis tombé un peu par hasard sur un tel exemple assez facilement calculable que je n'ai retrouvé nulle part avant (mais peut être l'avez vous déjà rencontré?) :




La fonction $f(x)=\sqrt{\vert \sin(x)\vert}$   est paire , $\pi$ périodique et ${\mathbb L}^2([0;\pi])$ elle admet donc un développement en séries de Fourier , au moins au sens d'une décomposition dans la base Hilbertienne des $\cos(nx)$ et $\sin(nx)$. Cette fonction est continue  mais pas $C^1$ par morceaux,  car sa dérivée  n'admet pas de limites finies à droite et à gauche aux points $(k\pi)_{k\in{\mathbb Z}}$. A cause de la singularité de $x\mapsto \sqrt{x}$ on a $\lim_{x\to k\pi^\pm}f'(x)=\pm \infty$ ce qui correspond aux tangentes verticales en ces points


Ce  développement en séries de Fourier ne contient que des termes en $\cos(nx)$  (avec seulement les $n$ pairs en fait) dont les coefficients sont exprimés par  la formule usuelle sur laquelle on va faire une intégration par parties :

\begin{align*}
a_n&={2\over\pi}\int_0^\pi \sqrt{\sin(x)}\cos(nx) \,dx\\
&= {2\over\pi}\left(\left[\sqrt{\sin(x)}{\sin(nx)\over n }\right]_0^\pi-  \int_0^\pi {\cos(x)\over 2\sqrt{\sin(x)}}{\sin(nx)\over n } \,dx\right)\\
&= -{1\over n\pi}  \int_0^\pi {\cos(x)\over \sqrt{\sin(x)}}\sin(nx) \,dx=-{1\over n\pi}  b_n\\
\end{align*}
on a posé  $b_n=-\pi na_n $ et on va calculer astucieusement (grâce aux formules d'addition des $\cos$ et $\sin$) :

\begin{align*}
b_{n+1}-b_{n-1}
&= \int_0^\pi {\cos(x)\over \sqrt{\sin(x)}}(\sin((n+1)x)-\sin((n-1)x)) \,dx\\
&= \int_0^\pi {\cos(x)\over \sqrt{\sin(x)}}(2\sin(x)\cos(nx)) \,dx\\
&= \int_0^\pi  \sqrt{\sin(x)}(2\cos(x)\cos(nx)) \,dx\\
&= \int_0^\pi  \sqrt{\sin(x)}(\cos((n+1)x)+\cos((n-1)x)) \,dx\\
&= {\pi\over 2}(a_{n+1}+a_{n-1})
\end{align*}
on a donc obtenue une équation récurrente linéaire sur les $a_n$ :
\begin{align*}
a_{n+1}+a_{n-1}=-(2n+2)a_{n+1}+(2n-2)a_{n-1}
\Leftrightarrow
a_{n+1}={2n-3\over 2n+3}a_{n-1}
\end{align*}
ceci permet de trouver facilement une forme close pour les $a_{2n}$ sous forme d'un produit (les termes impairs étant nuls) :
\begin{align*}
a_{2n}&=a_{(2n-1)+1}={4n-5\over 4n+1}a_{2(n-1)}\\
&={4n-5\over 4n+1}\times \dots\times {3\over 5}\times {-1\over 1}a_{0}\\
&=\prod_{k=1}^n {4k-5\over 4k+1}a_0=-\prod_{k=2}^n {4k-5\over 4k+1}a_0
\end{align*}
où la valeur $a_0$  s'exprime avec la fonction $\beta$
\begin{align*}
a_0&={2\over\pi}\int_0^\pi \sqrt{\sin(x)} \,dx
={2\over\pi}\beta\left({1\over 2},{3\over 4}\right)= 1.525519527003627\dots
\end{align*}

il reste à faire l'évaluation asymptotique  de $a_{2n}$ . Puisque c'est un produit  on passe par le logarithme et on est ramené à étudier une série divergente :

\begin{align*}
\ln\left(\prod_{k=2}^n {4k-5\over 4k+1}\right)
&=\ln\left(\prod_{k=1}^n {1-{5\over 4k}\over 1+{1\over 4k}}\right)\\
&=\sum_{k=2}^n \ln\left(1-{5\over 4k}\right)-\ln\left(1+{1\over 4k}\right)\\
&=\sum_{k=2}^n -{5\over 4k}-{1\over 4k}+{\cal O}\left({1\over k^2}\right)\\
&=\sum_{k=2}^n -{3\over 2k}+{\cal O}\left({1\over k^2}\right)\\
&= -{3\over 2}\ln(n)+\ln(C)+{\cal O}\left({1\over n}\right)\\
\end{align*}
où on a utilisé que $\sum_{k=1}^n{1\over k}= \ln(n)+\gamma+{\cal O}\left({1\over n}\right)$. On peut finalement conclure que
$$ a_{2n}= -{C\over n\sqrt{n}}\left(1+{\cal O}\left({1\over n}\right)\right)\sim_\infty -{C\over n\sqrt{n}}$$

Faute d'autres idées pour calculer la constante $C$  j'obtiens le résultat final :

 $$\sqrt{\vert \sin(x)\vert}=\underbrace{{2\over\pi}\beta\left({1\over 2},{3\over 4}\right)}_{a_0\approx 1.525519527}\left({1\over 2}+\sum_{n=1}^{\infty} \underbrace{\prod_{k=1}^n {4k-5\over 4k+1} }_{{a_{2n}\over a_0}\sim_\infty -{C\over n\sqrt{n}}}\cos (2n x)\right)~~~~ \forall x\in {\mathbb R}$$

Au départ l'idée m'est venue d'un exercice posé par un collègue, qui se traite avec des calculs tout à fait similaires :
$$\ln(\vert \sin(x)\vert)= -\ln(2)-\sum_{n=1}^\infty {1\over n} \cos(2nx)$$
j'ai aussi trouvé un autre exemple de ce type sur stackexchange  avec la fonction impaire définie sur $[0,\pi] $ par $f(x)=\ln (\sqrt{1+\sin (x)}+\sqrt{\sin (x)})$ alors on a:
 $$
f(x)=\sum_{n=0}^{\infty} \frac{4^{-n}(2 n) !}{(2 n+1)(n !)^{2}} \sin ((2 n+1) x)
 $$

Le calcul est basée sur la même approche, après une intégration par parties on trouve une expression qui doit permettre de retomber sur une récurrence linéaire (que je n'ai pas trouvé !) :

$$
b_{n}=\frac{2}{\pi} \int_{0}^{\pi} \ln (\sqrt{1+\sin (x)}+\sqrt{\sin (x)}) \sin (n x) d x
=\frac{2}{\pi n} \int_{0}^{\pi / 2} \frac{\cos(x) \cos (n x)}{\sqrt{\sin (x)(1+\sin (x))}} d x $$


en utilisant la formule de Stirling on obtient un comportement asymptotique comparable à celle de $\sqrt{\vert \sin(x)\vert}$ :


 $$\frac{4^{-n}(2 n) !}{(2 n+1)(n !)^{2}}
\sim_\infty
4^{-n}{\left({2n\over e}\right)^{2n}\sqrt{4\pi n}\over(2n+1) \left({n\over e}\right)^{2n}2\pi n}
={\sqrt{2}\over(2n+1) \sqrt{\pi n}}
\sim_\infty{1\over \sqrt{2\pi }}{1\over n\sqrt{ n}}
$$

si vous avez d'autres exemples de comportement exotiques pour des séries de Fourier n'hésitez pas à les partager ici même!


Aucun commentaire:

Enregistrer un commentaire

Pour écrire des formules mathématiques vous pouvez utiliser la syntaxe latex en mettant vos formules entre des "dollars" $ \$....\$ $ par exemple :
- $\sum_{n=1}^\infty {1\over n^2}={\pi^2\over 6}$ s'obtient avec \sum_{n=1}^\infty {1\over n^2}={\pi^2\over 6}
- $\mathbb R$ s'obtient avec {\mathbb R} et $\mathcal D$ s'obtient avec {\mathcal D}
- pour les crochets $\langle .,. \rangle$ dans les commentaires utilisez \langle .,. \rangle
vous pouvez écrire du html dans les commentaires :
- italique <i> ... </i> gras <b> ... </b>
- lien <a href="http://adresse "> .... </a>