Dans notre système universitaire, l'enseignement des outils de comparaison asymptotique se limite en général à l'étude des développements limités (DL). Si cet outil est incontournable son apprentissage se réduit bien souvent à connaître par cœur une liste de DL usuels (les séries géométrique et exponentielle ou la formule du binôme étendue aux exposants non entiers) et quelques techniques de combinaison. La question des développements asymptotiques (DA) est souvent bien plus difficile à aborder, c'est bien dommage, car elle fourni des problèmes pas forcément plus difficiles techniquement mais plus profond conceptuellement et qui illustrent parfaitement cette citation de Bertrand Russell :
Je vous propose donc d'étudier le problème suivant , familier de ceux qui connaisse le théorème de raréfaction des nombres premier :
image wikipédia |
«trouver un développement asymptotique de la solution
$x$ de l'équation $x\ln(x)=t$ quand $t\to\infty$ »
$x$ de l'équation $x\ln(x)=t$ quand $t\to\infty$ »
Pour résoudre de manière asymptotique ce type d'équation il faut d'abord s'assurer que le problème a bien une solution unique, par exemple en faisant l'étude de la fonction :
$$f(x)=x\ln(x)\Rightarrow f'(x)=\ln(x)+1>0 \text{ si } x\geq 1$$
on en déduit que $f$ est bijective (car monotone) de $ [1,+\infty[$ dans $ f( [1,+\infty[)=[0,+\infty[$ on est donc certain que quand $t\to\infty$ on a bien une unique solution $x(t)\geq 1$, mais on peut aussi en déduire que $x(t)\to\infty$ quand $t\to\infty$. Cette simple remarque va permettre d'enclencher une stratégie dite «d'amorçage/pompage»* pour trouver un développement asymptotique précis de $x=x(t)$. Pour résumer, cette stratégie se décompose en 2 étapes :
- l'amorçage : on cherche un premier équivalent simple de la solution voir une domination asymptotique du reste
- le pompage : on réinjecte l'estimation asymptotique dans l'équation de départ pour en sortir une nouvelle estimation avec un terme de reste plus précis
$$x\ln(x)=t\Rightarrow x={t\over \ln(x)}$$
pour connaitre l'équivalent de $x$ il nous faut donc trouver un équivalent de $\ln(x)$, pour ça on passe au logarithme dans l'équation de départ:
$$ x\ln(x)=t
\Rightarrow\ln(x)+\underbrace{\ln(\ln(x))}_{=o(\ln(x))}=\ln(t)
\Rightarrow\ln(x)\sim_\infty \ln(t) $$
on en déduit un premier équivalent de la solution :
$$ x={t \over \ln(x)}\sim_\infty {t \over \ln(t)}$$
Pour obtenir un DA on va maintenant chercher à mettre en facteur $t/\ln(t)$ dans l'expression de $x$ :
$$ x={t \over \ln(x)}= {t \over \ln({t \over \ln(x)})}
= {t \over \ln(t)-\ln(\ln(x))}= {t \over \ln(t)}\times {1\over 1-{\ln(\ln(x))\over \ln(t)}}$$
comme $\ln(x)\sim_\infty \ln(t)$ on a aussi $\ln(\ln(x))\sim \ln(\ln(t))$ (voir le billet sur le passage au logarithme dans les équivalents) et on peut poser $u={\ln(\ln(x))\over \ln(t)}\sim_\infty {\ln(\ln(t))\over \ln(t)}\to0$ quand $t\to \infty$ pour majorer le reste via un DL de ${1\over 1-u}=1+O(u)$ quand $u\to0$ :
$$ x={t \over \ln(t)}\times {1\over 1-{\ln(\ln(x))\over \ln(t)}}
={t \over \ln(t)}\times \left(1+O\left({\ln(\ln(t))\over \ln(t)}\right)\right)
={t \over \ln(t)}+O\left({t\ln(\ln(t))\over \ln(t)^2}\right)
$$
On a finit l'étape d'amorçage! Il ne reste "plus" qu'à recommencer pour remplacer le terme ${\ln(\ln(x))\over \ln(t)}$ par un équivalent et améliorer l'estimation du reste. D'abord on passe au logarithme :
$$\begin{align*}
{\ln(\ln(x))}&={\ln(\ln(t/\ln(x)))} ={\ln(\ln(t)-\ln(\ln(x)))}\\
&= {\ln\left(\ln(t)\left(1-{\ln(\ln(x))\over\ln(t)}\right)\right)}\\
&= {\ln(\ln(t))+\ln\left(1-{\ln(\ln(x))\over\ln(t)}\right)}\\
&= {\ln(\ln(t))}+O\left({\ln(\ln(x))\over\ln(t)}\right)
= {\ln(\ln(t))}+O\left({\ln(\ln(t))\over\ln(t)}\right)\\
\end{align*}
$$
où on a utilisé que $\ln(1-u)=O(u)$ quand $u={\ln(\ln(x))\over\ln(t)}\sim_\infty {\ln(\ln(t))\over\ln(t)}\to 0$ . Ensuite on divise par $\ln(t)$
$$\begin{align*}
{\ln(\ln(x))\over\ln(t)}
&= {\ln(\ln(t))\over\ln(t)}+O\left({\ln(\ln(t))\over\ln(t)^2}\right)
\end{align*}$$
et on utilise le DL ${1\over 1-u}=1+u+O(u^2) $ en $u={\ln(\ln(x))\over \ln(t)}\to 0$ :
$$\begin{align*}
{1\over 1-{\ln(\ln(x))\over \ln(t)}}
&= 1+{\ln(\ln(x))\over \ln(t)} +O\left(\left({\ln(\ln(x))\over \ln(t)}\right)^2\right)\\
&= 1+{\ln(\ln(t))\over \ln(t)} +O\left({\ln(\ln(t))\over\ln(t)^2}\right) +\underbrace{O\left(\left({\ln(\ln(t))\over \ln(t)}\right)^2\right)}_\text{terme dominant}\\
\end{align*}$$
enfin on multiplie les 2 résultats :
$$\begin{align*}x&={t \over \ln(t)}\times {1\over 1-{\ln(\ln(x))\over \ln(t)}}&={t \over \ln(t)}+{t\ln(\ln(t))\over \ln(t)^2} +O\left({t\ln(\ln(t))^2\over\ln(t)^3}\right)
\end{align*}
$$
En continuant on trouvera une série de termes en ${t\ln(\ln(t))^k\over \ln(t)^n}$ avec $1\leq k< n$ et un reste en $O(t\ln(\ln(t))^n/ \ln(t)^{n+1})$ ... Je vous laisse réitérer pour trouver par pompage les termes suivant du DA :-)
$$\begin{align*}x&=
{t \over \ln(t)}+{t\ln(\ln(t))\over \ln(t)^2}+{t\ln(\ln(t))^2\over \ln(t)^3}-{t\ln(\ln(t))\over \ln(t)^3} +O\left({t\ln(\ln(t))^3\over\ln(t)^4}\right)
\end{align*}$$
vous pouvez facilement juger de la précision du DA limité à 4 termes en prenant un exemple numérique:
-->x=1000 x = 1000. -->t=x*log(x) t = 6907.7553 -->x1=t/log(t) x1 = 781.38492 -->x2=x1+t*log(log(t))/(log(t)^2) x2 = 974.01164 -->x3=x2+t*(log(log(t))^2)/(log(t)^3) x3 = 1021.4979 -->x4=x3-t*(log(log(t)))/(log(t)^3) x4 = 999.70854
* je me souviens qu'on m'a introduit ce genre de méthode sous ce nom «d'amorçage/pompage» mais je ne trouve aujourd'hui aucune trace de cette terminologie ...
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>