MENUCours D'introduction à l'analyse numérique

Introduction

Les intégrales elliptiques complètes de première et seconde espèce sont respectivement définies par \[ K(k)=\int_0^{\pi/2}\frac{\mathrm{d}\varphi}{\sqrt{1-k^2\sin^2\varphi}} \quad\text{et}\quad E(k)=\int_0^{\pi/2}\sqrt{1-k^2\sin^2\varphi}\,\mathrm{d}\varphi \] avec \(k\in[0;1[\).

Elles ont été étudiées par Adrien-Marie Legendre en 1792 et interviennent dans de nombreux problèmes de physique comme l'étude du pendule simple, le champ magnétique créé par une spire, la mutuelle inductance entre deux spires, la perturbation des planètes, etc. Ces intégrales ne s’expriment pas en termes de fonctions simples. Bien sûr on peut utiliser une méthode numérique de calcul d’intégrale pour en donner une valeur approchée, mais la convergence n’est pas aussi rapide que celle de la méthode que l’on propose ici. En effet, ces deux intégrales sont étroitement liées à la moyenne arithmético-géométriquearithmetic-geometric mean en anglais, souvent noté AGM. ce qui permet en quelques lignes de code et quelques itérations de calculer \(K(k)\) et \(E(k)\) avec une grande précision.

Moyenne arithmético-géométrique

Considérons les suites réelles \((a_{n})\) et \((b_{n})\) définies par les relations \[ \left\{ \begin{array}{ccc} a_{n} &=&\frac{1}{2}\left(a_{n-1}+b_{n-1}\right) \\[3mm] b_{n} &=&\sqrt{a_{n-1}b_{n-1}} \\ \end{array} \right. \quad\text{avec}\quad \left\{ \begin{array}{ccc} a_{0} &=&a >0 \\[3mm] b_{0} &=&b <a \\ \end{array} \right. \] Ces suites, comme on le constate dans le très viteOn peut montrer que \[a_{n+1}-b_{n+1}\leq \frac{1}{8b}(a_{n}-b_{n})^{2}\] d'où une convergence quadratique. vers une limite commune \(M_{a,b}\) dite moyenne arithmético-géométrique.

Intégrales elliptiques et AGM

Intégrale de première espèce

Considérons l'intégrale suivante : \[ I(a,b)=\int_{0}^{\pi/2}\frac{\mathrm{d} \varphi}{\sqrt{a^{2}\cos^{2}\varphi+b^{2}\sin^{2}\varphi}} \] Cette intégrale présente la propriété remarquable d'être invariante par la transformation \(a\mapsto (a+b)/2\) et \(b\mapsto \sqrt{ab}\).

Par conséquent, si l'on appelle \((a_n)\) et \((b_n)\) les termes de la suite arithmético-géométrique générés par les termes initiaux \(a\) et \(b\), l'intégrale \(I(a,b)\) vérifie la propriété suivante \[ I(a,b)=I(a_{1},b_{1})=\ldots=I(a_{n},b_{n})=\ldots=I(M_{a,b},M_{a,b}) \] cette dernière intégrale se calculant sans difficulté : \[ I(a,b)= I(M_{a,b},M_{a,b})=\int_{0}^{\pi/2}\frac{\mathrm{d} \varphi}{M_{a,b}}=\frac{\pi}{2M_{a,b}} \] Il nous reste plus qu'à relier l'intégrale elliptique de première espèce à \(I(a,b)\). Il est facile de vérifier que \[ K(k)=I(1,\sqrt{1-k^{2}}) \quad\text{car}\quad \cos^{2}\varphi+(1-k^{2})\sin^{2}\varphi=1-k^{2}\sin^{2}\varphi \] Retenons le résultat suivant

Intégrale de seconde espèce

En plus de \(I(a,b)\), considérons deux nouvelles intégrales : \[ J(a,b)=\int_0^{\frac{\pi}{2}}\sqrt{a^2\cos^2\varphi+b^2\sin^2\varphi}\, \mathrm{d}\varphi \quad\text{et}\quad L(a,b)=\int_0^{\frac{\pi}{2}}\frac{\cos^2\varphi\, \mathrm{d}\varphi}{\sqrt{a^2\cos^2\varphi+b^2\sin^2\varphi}} \] À l'aide des deux changements de variable rencontrés dans la section précédente, on peut montrer [5] la propriété suivante : \begin{equation} {c_0}^2 L(b,a)=S I(a,b) \quad\text{avec}\quad \left\lbrace \begin{array}{rcl} {c_n}^2 &=&{a_n}^2-{b_n}^2\\ S &=&\sum_{n=0}^{\infty}2^{n-1}{c_n}^2\\ \end{array} \right. \end{equation} où \(a_n\) et \(b_n\) sont les termes de la suite arithmético-géométrique initiée par \(a\) et \(b\). Explicitons le premier terme \({c_0}^2L(b,a)\) : \[ {c_0}^2=a^2-b^2 \quad\text{et}\quad L(b,a)=\int_0^{\pi/2}\frac{\cos^2\varphi\, \mathrm{d}\varphi}{\sqrt{b^2\cos^2\varphi+a^2\sin^2\varphi}} \] Le changement de variable \(\phi=\frac{\pi}{2}-\varphi\) transforme \(L(b,a)\) comme suit : \[ L(b,a)=\int_0^{\pi/2}\frac{\sin^2\phi\, \mathrm{d}\phi}{\sqrt{a^2\cos^2\phi+b^2\sin^2\phi}} \] de sorte que \({c_0}^2L(b,a)\) vaut \begin{equation} \begin{split} {c_0}^2L(b,a)=\int_0^{\frac{\pi}{2}}\frac{a^2\sin^2\phi-b^2\sin^2\phi}{\sqrt{a^2\cos^2\phi+b^2\sin^2\phi}}\, \mathrm{d}\phi= \int_0^{\frac{\pi}{2}}\frac{a^2\, \mathrm{d}\phi}{\sqrt{a^2\cos^2\phi+b^2\sin^2\phi}}\\ -\int_0^{\frac{\pi}{2}}\sqrt{a^2\cos^2\phi+b^2\sin^2\phi}\,\mathrm{d}\phi=a^2 I(a,b)-J(a,b) \end{split} \notag \end{equation}

En reprenant la relation (1), on aboutit à \begin{equation} \boxed{\hspace{0.5em} J(a,b)=(a^2-S)I(a,b) \hspace{0.5em}} \notag \end{equation}

Il suffit maintenant de choisir \(a=1\) et \(b=\sqrt{1-k^2}\) pour obtenir la relation entre les intégrales elliptiques complètes de première et seconde espèce :

Méthode numérique

Comme on l'a vu, les suites \((a_n)\) et \((b_n)\) convergent vers la moyenne \(M_{a,b}\) à une vitesse prodigieuse. Quant à la suite \(({c_n}^2)\), elle converge vers 0 tout aussi rapidement (Tab. 2).

On peut donc obtenir la moyenne \(M_{a,b}\) et \(S\) avec une grande précision grâce à quelques itérations de la suite arithémitco-géométrique.

Voici un premier algorithme qui calcule \(K(k)\) avec un certain seuil de précision, noté \(\epsilon_0\):

Pour le calcul de \(E(k)\), on doit, au fur et à mesure des itérations, calculer \(S\). C'est sur cette somme qu'un critère de précision sera appliqué. D'autre part, si \(k=1\), il n'est pas nécessaire d'utiliser la suite arithmético-géométrique car l'on connaît la valeur théoriqueContrairement à \(K(k)\), \(E(k)\) est défini pour \(k=1\). En effet, \(E(1)=\int_0^{\pi/2}\cos\varphi\,\mathrm{d}\varphi=1\)