Transformation de Landen

John Landen (-)

Les transformations de Landen transforment une intégrale elliptique (ou une fonction elliptique jacobienne) en une autre intégrale elliptique (ou une autre fonction elliptique jacobienne) de même forme, mais avec des variables différentes. Les transformations croissantes augmentent le module des fonctions elliptiques tandis que les transformations décroissantes le diminuent. Elles permettent d'effectuer plus rapidement les calculs numériques des fonctions elliptiques.

Historique

Cette transformation est l'œuvre initiale du mathématicien anglais John Landen (1719-1790) qui proposa en 1775 un changement de variable très réussi pour les intégrales et fonctions elliptiques.

Il a pu montrer que la longueur d'un arc d'hyperbole pouvait être exprimée par les longueurs des arcs de deux ellipses différentes. Chacune d’elles présente une excentricité qui peut être identifiée au module elliptique. La relation particulière entre les excentricités des deux ellipses que Landen découvrit plus tard porte son nom.

Cette transformation a été redécouverte indépendamment par Carl Friedrich Gauss. La forme actuelle de la transformation de Landen a été développée par Joseph-Louis Lagrange, Adrien-Marie Legendre et Gauss. En utilisant la transformation de Landen, Gauss a calculé la longueur de la lemniscate[réf. souhaitée]. En particulier, les travaux de Legendre ont joué un rôle majeur pour Niels Henrik Abel et Charles Gustave Jacob Jacobi dans leur développement des fonctions elliptiques.

Landen a découvert une nouvelle façon de calculer, et pas seulement les fonctions elliptiques. Son idée principale, selon laquelle la fonction calculée peut être représentée par une fonction de la même forme mais avec d'autres paramètres qui tendent vers certaines limites lors de la récursion, a ensuite été largement utilisée en mathématiques computationnelles. Cette transformation modulaire joue un rôle important dans les mathématiques modernes.

Changement de variable et nouveau module

Au sens large, la transformation de Landen désigne une transformation reposant sur le principe selon lequel la fonction calculée peut être représentée par une fonction de la même forme mais avec d'autres paramètres qui tendent vers certaines limites lors de la récursion ; tandis que la transformation de Landen à proprement parler désigne la transformation qui utilise le changement de variable que Landen a proposé[réf. souhaitée].

On pose :

{k0=kk'0=k′θ0=θ.

Les transformations changent le module k0 en un autre module k1 ou k−1 en changeant la variable d'intégration θ0 en une nouvelle variable θ1 ou θ−1 définie ainsi :

sin⁡(2θn+1−θn)=knsin⁡θn⇔sin⁡θn=21+knsin⁡θn+1cos⁡θn+11−(2kn1+kn)2sin2θn+1=(1+kn′+1′)sin⁡θn+1cos⁡θn+11−kn+12sin2θn+1
sin⁡θn+1=(1+kn)sin⁡θn1+knsin2θn⇔sin⁡θn=21+knsin⁡θn+11+1−(2kn1+kn)2sin2θn+1=(1+kn′+1′)sin⁡θn+11+1−kn+12sin2θn+1

Ainsi, pour la transformation de Landen comme pour la transformation gaussienne, on a les relations suivantes :

{kn+1=2kn1+kn=21−kn'241+1−kn'2,kn′+1′=1−kn1+kn=1−1−kn'21+1−kn'2⇔kn−1=1−1−kn21+1−kn2=1−kn′1+kn′,kn′−1′=21−kn241+1−kn2=2kn′1+kn′

Une transformation qui change kn en kn+1 est une transformation croissante et une transformation qui change kn en kn−1 est une transformation décroissante parce que :

{kn⩽1⇒kn+1⩾kn⇒kn+1⩾knkn−1=kn2(1+kn′)2⇒kn−1⩽kn

On peut aussi exprimer ainsi la relation de récurrence :

kn−1=tan2arcsin⁡kn2=th2argth⁡kn2

Transformations de Landen des intégrales elliptiques de première espèce

Transformation de Landen

Démonstration

On commence par vérifier la réciprocité de :

sin⁡(2θn+1−θn)=knsin⁡θn⇔sin⁡θn=21+knsin⁡θn+1cos⁡θn+11−(2kn1+kn)2sin2θn+1

On a :

sin⁡θn=sin⁡(2θn+1)1+kn2+2kncos⁡(2θn+1),cos⁡θn=kn+cos⁡(2θn+1)1+kn2+2kncos⁡(2θn+1),tan⁡θn=sin⁡(2θn+1)kn+cos⁡(2θn+1)

En multipliant cette dernière égalité par cos⁡θn[kn+cos⁡(2θn+1)], on a bien : sin⁡(2θn+1−θn)=knsin⁡θn. Ce changement de variable permet que l'angle transformé devienne plus petit que l'angle d'origine : kn⩽1⇒2θn+1−θn⩽θn⇒θn+1⩽θn.

Les grandeurs apparaissant dans l'intégrale elliptique de première espèce sont :

1−kn2sin2θn=1−sin2(2θn+1−θn)=cos⁡(2θn+1−θn)=1+kncos⁡(2θn+1)1+kn2+2kncos⁡(2θn+1)
dθn=cos2θndtan⁡θn=2[1+kncos⁡(2θn+1)]1+kn2+2kncos⁡(2θn+1)dθn+1

Ainsi :

dθn1−kn2sin2θn=21+kndθn+11−kn+12sin2θn+1

Amplitude

On calcule l'amplitude φn+1L (on pose : φ0L=φ) :

sin2φnL=1−cos2(2φn+1L)1+kn2+2kncos⁡(2φn+1L)
⇒cos2(2φn+1L)+2knsin2φnLcos⁡(2φn+1L)+kn2sin2φnL−1+sin2φnL=0
⇒cos⁡(2φn+1L)=−knsin2φnL±kn2sin4φnL−kn2sin2φnL+1−sin2φnL∈[−1;+1]
⇒cos⁡(2φn+1L)=−knsin2φnL+cos⁡φnL1−kn2sin2φnL
⇒φn+1L=arcsin⁡1+knsin2φnL−cos⁡φnL1−kn2sin2φnL2

On calcule l'amplitude φnL. Puisque φn+1L∈[φnL2;φnL]⇔φnL∈[φn+1L;2φn+1L], on a :

{φn+1L∈[0;π2−arccos⁡kn2]⇔φnL∈[0;π2]φn+1L∈[π2−arccos⁡kn2;π2]⇔φnL∈[π2;π]
φnL=arccos⁡kn+cos⁡(2φn+1L)1+kn2+2kncos⁡(2φn+1L)=arccos⁡1−(1+k'n+1)sin2φn+1L1−kn+12sin2φn+1L

Si on utilise la notation de l'intégrale elliptique avec un point-virgule, on est obligé de décomposer l'intégrale en tronçons ne dépassant pas π/2 puisque cette notation ne permet que d'exprimer des intégrales elliptiques dont l'amplitude ne dépasse pas π/2 : si φnL∈]π/2;π], F(φnL,kn)=2K(kn)−F(π−φnL,kn)=2K(kn)−F(sin⁡φnL;kn).

Autre formulation

Si kn, kn′+1′, φnL et φn+1L sont tels que (1+kn)(1+kn′+1′)=2 et tan⁡(φnL−φn+1L)=kn′+1′tan⁡φn+1L, alors la transformation de Landen stipule que :

F(φnL,kn)=21+knF(φn+1L,kn+1)=(1+kn′+1′)F(φn+1L,kn+1)

La transformation de Landen peut donc être exprimée soit en fonction de son module elliptique k, soit en fonction de son comodule k′.

En effet :

  • En réécrivant sin⁡(2φn+1L−φnL)=knsin⁡φnL :
sin⁡φn+1Lcos⁡(φn+1L−φnL)+cos⁡φn+1Lsin⁡(φn+1L−φnL)=knsin⁡(φnL−φn+1L)cos⁡φn+1L+kncos⁡(φnL−φn+1L)sin⁡φn+1L
et en divisant par cos⁡φn+1Lcos⁡(φn+1L−φnL), on a :
tan⁡φn+1L−tan⁡(φnL−φn+1L)=kntan⁡(φnL−φn+1L)+kntan⁡φn+1L
ce qui donne :
tan⁡(φnL−φn+1L)=kn′+1′tan⁡φn+1L⇒tan⁡φnL−tan⁡φn+1L1+tan⁡φnLtan⁡φn+1L=kn′+1′tan⁡φn+1L⇒tan⁡φnL=(1+kn′+1′)tan⁡φn+1L1−kn′+1′tan2φn+1L
  • De plus, on a :
(1+kn)(1+kn′+1′)=2

Transformation gaussienne

Démonstration

On commence par vérifier la réciprocité de :

sin⁡θn+1=(1+kn)sin⁡θn1+knsin2θn⇔sin⁡θn=21+knsin⁡θn+11+1−(2kn1+kn)2sin2θn+1

On a bien :

sin⁡θn=2sin⁡θn1+knsin2θn1+1−(2kn1+kn)2(1+kn)2sin2θn(1+knsin2θn)2=2sin⁡θn1+knsin2θn1+1−knsin2θn1+knsin2θn=sin⁡θn

Les grandeurs apparaissant dans l'intégrale elliptique de première espèce sont :

1−kn2sin2θn=(1+kn−kncos2θn)2−(1+kn)2+(1+kn)2cos2θncos⁡θn=(1+knsin2θn)2−(1+kn)2sin2θncos⁡θn
dθn+1=dsin⁡θn+1cos⁡θn+1=(1+knsin2θn)(1+kn)cos⁡θn−2kn(1+kn)sin2θncos⁡θncos⁡θn+1(1+knsin2θn)2dθn=(1−knsin2θn)(1+kn)cos⁡θn(1+knsin2θn)(1+knsin2θn)2−(1+kn)2sin2θndθn

On a alors :

dθn1−kn2sin2θn=(1+knsin2θn)dθn+1(1−knsin2θn)(1+kn)=dθn+1(1+kn)1−4knsin2θn(1+knsin2θn)2=dθn+1(1+kn)1−(2kn1+kn)2sin2θn+1

Amplitude

On calcule l'amplitude φn+1 (on pose : φ0G=φ) :

φn+1G=arcsin⁡(1+kn)sin⁡φnG1+knsin2φnG

L'amplitude φnG vaut :

φnG=arcsin⁡2sin⁡φn+1G1+kn+1+kn2+2kncos⁡(2φn+1G)=arcsin⁡(1+k'n+1)sin⁡φn+1G1+1−kn+12sin2φn+1G

Transformation AGM quartique

En plus des changements de variable landenien et gaussien, il en existe d'autres, par exemple celui-ci :

tan⁡θn+1=1+sin⁡θn1−k24cos⁡θn[réf. à confirmer]

Itération

Quel que soit kn, φnL varie de 0 à π lorsque φn+1L varie de 0 à π/2.
Transformation de Landen
[0;π]→[0;π/2]φnL↦φn+1L
Cliquer pour voir et modifier le graphique.
Quel que soit kn, φnG varie de 0 à π/2 lorsque φn+1G varie de 0 à π/2.
Transformation de Gauss
[0;π/2]→[0;π/2]φnG↦φn+1G
Cliquer pour voir et modifier le graphique.
φn+1L(φnL) et φn+1G(φnG) pour kn∈[0;1] par pas de 0,05. La courbe rouge correspond au cas kn=0,8.

En répétant plusieurs fois de suite la transformation de Landen ou gaussienne, on aura k∞→1 si on utilise la transformation croissante, et k−∞→0 si on utilise la transformation décroissante. φ évolue ainsi : kn↗⇒φnL/G↘ et k−n↘⇒φ−nL/G↗. Lorsque le module est égal à 0 ou 1, l'intégrale elliptique peut être calculée analytiquement :

F(φ,0)=φ
F(φ,1)=∫0φdθ1−sin2θ=artanh⁡sin⁡φ=12ln⁡1+sin⁡φ1−sin⁡φ=ln⁡1+sin⁡φcos⁡φ=ln⁡cos⁡φ2+sin⁡φ2cos⁡φ2−sin⁡φ2=ln⁡1+tan⁡φ21−tan⁡φ2=ln⁡tan⁡(φ2+π4)

Si l'on part d'un module k et d'une amplitude φ arbitraires, une intégrale elliptique générale de première espèce peut être calculée numériquement ainsi :

  • Pour la transformation de Landen :
F(φ,k)=21+kF(φ1L,k1)=21+k21+k1F(φ2L,k2)=21+k0⋯21+kn−1F(φnL,kn)=k1kF(φ1L,k1)=k1kk2k1F(φ2L,k2)=knk0⋯kn−1kF(φnL,kn)=1+k−12F(φ−1L,k−1)=1+k−121+k−22F(φ−2L,k−2)=1+k−12⋯1+k−n2F(φ−nL,k−n)=k−1kF(φ−1L,k−1)=k−1kk−2k−1F(φ−2L,k−2)=k−nkk−1⋯k−nF(φ−nL,k−n)
F(φ,k)=ln⁡tan⁡(φ∞L2+π4)∏n=0∞21+kn=ln⁡tan⁡(φ∞L2+π4)∏n=0∞knk=limn→∞φ−nL∏i=1n1+k−i2=limn→∞φ−nLk−nk∏i=1nk−i
  • Pour la transformation gaussienne :
F(φ,k)=11+kF(φ1G,k1)=11+k11+k1F(φ2G,k2)=11+k0⋯11+kn−1F(φnG,kn)=k12kF(φ1G,k1)=k12kk22k1F(φ2G,k2)=knk0⋯kn−12nkF(φnG,kn)=(1+k−1)F(φ−1G,k−1)=(1+k−1)(1+k−2)F(φ−2G,k−2)=(1+k−1)⋯(1+k−n)F(φ−nG,k−n)=2k−1kF(φ−1G,k−1)=2k−1k2k−2k−1F(φ−2G,k−2)=2nk−nkk−1⋯k−nF(φ−nG,k−n)
F(φ,k)=limn→∞ln⁡tan⁡(φnG2+π4)∏i=0n11+ki=limn→∞ln⁡tan⁡(φnG2+π4)∏i=0nkn2nk=φ−∞G∏n=1∞(1+k−n)=φ−∞Glimn→∞2nk−nk∏i=1nk−i

où φ∞L, φ−∞L, φ∞G et φ−∞G sont les valeurs asymptotiques de l'amplitude transformée. On remarquera que φ−nL double à chaque itération lorsque n→∞.

Les convergences de kn et k−n sont quadratiques : −lg⁡(1−kn) et −lg⁡k−n doublent à peu près à chaque itération, ce qui signifie que peu d'itérations suffisent : k=0,1⇒1−k5<10−14, k=10−100⇒1−k11<10−14, k=0,9⇒k−5<10−14 et k=1−10−14⇒k−8<10−14.

L'intégrale elliptique complète de première espèce est :

K(k)=F(π2,k)=1+k−12F(π,k−1)=(1+k−1)K(k−1)=(1+k−1)(1+k−2)K(k−2)=π2∏n=1∞(1+k−n)

Transformation de Landen du nom elliptique

On a :

K(kn+1)=(1+kn)K(kn)⇒K(2kn1+kn)=(1+kn)K(kn)⇒K(2kn′1+kn′)=(1+kn′)K(kn′)⇒K(kn′−1′)=(1+kn′)K(kn′)⇒
K(kn′)=(1+kn′+1′)K(kn′+1′)
⇒K(k'n)K(kn)=2K(kn′+1′)K(kn+1)
⇒q(kn)2=q(kn−1)

Transformations de Landen des intégrales elliptiques de deuxième espèce

Transformation de Landen

Démonstration

On pose (avec a>b), :

{F(a,b,φ)=∫0φdθa2cos2θ+b2sin2θE(a,b,φ)=∫0φa2cos2θ+b2sin2θdθ et {an−1=an+bn2bn−1=anbncn−1=an−bn2

Le module devient kn=1−(bnan)2=cnan, si bien que kn−1=an−bnan+bn. On a an2=bn2+cn2. a−1 et b−1 sont respectivement la moyenne arithmétique et géométrique de a et b. Si la transformation est itérée plusieurs fois, alors les paramètres a−i et b−i convergent très rapidement vers une valeur commune, même s’ils sont initialement d’ordres de grandeur différents. La valeur limite est appelée moyenne arithmético-géométrique de a et b et notée AGM⁡(a,b) ou M(a,b). On a alors a−∞=b−∞=AGM⁡(a,b). On a :

21+kn−1dθn1−kn2sin2θn=dθn−11−kn−12sin2θn−1⇒2dθnan2cos2θn+bn2sin2θn=dθn−1an−12cos2θn−1+bn−12sin2θn−1

et :

cos⁡(2θnL)=−kn−1sin2θn−1L+cos⁡θn−1L1−kn−12sin2θn−1L
⇒{2cos2θnL=1−kn−1sin2θn−1L+cos⁡θn−1L1−kn−12sin2θn−1L2sin2θnL=1+kn−1sin2θn−1L−cos⁡θn−1L1−kn−12sin2θn−1L
⇒2an2cos2θnL+2bn2sin2θnL=(an2+bn2)(cos2θn−1L+sin2θn−1L)−(an2−bn2)kn−1sin2θn−1L+(an2−bn2)cos⁡θn−1L1−kn−12sin2θn−1L=(an2+bn2)cos2θn−1L+2anbnsin2θn−1L+(an−bn)cos⁡θn−1L(an+bn)2−(an−bn)2sin2θn−1L=4(an−12cos2θn−1L+bn−12sin2θn−1L)−2bn−12+4cn−1cos⁡θn−1Lan−12cos2θn−1L+bn−12sin2θn−1L

En intégrant, on obtient :

F(an,bn,φnL)=12F(an−1,bn−1,φn−1L)=12mF(an−m,bn−m,φn−mL)
E(an,bn,φnL)=E(an−1,bn−1,φn−1L)−bn−122F(an−1,bn−1,φn−1L)+cn−1sin⁡φn−1L
⇒E(φnL,kn)=an−1anE(φn−1L,kn−1)−bn−122anan−1F(φn−1L,kn−1)+cn−1ansin⁡φn−1L
⇒E(φnL,kn)=1+kn′2E(φn−1L,kn−1)−kn′1+kn′F(φn−1L,kn−1)+1−kn′2sin⁡φn−1L
⇒E(φnL,kn)=11+kn−1E(φn−1L,kn−1)−1−kn−12F(φn−1L,kn−1)+kn−11+kn−1sin⁡φn−1L

On fera attention comme précédemment si on utilise les notations avec un point-virgule.

Itération

E(a,b,φ)−a2F(a,b,φ)=E(a−1,b−1,φ−1L)−a−12F(a−1,b−1,φ−1L)+(a−12−a22−b−122)F(a−1,b−1,φ−1L)+c−1sin⁡φ−1L=E(a−1,b−1,φ−1L)−a−12F(a−1,b−1,φ−1L)−a−1c−1F(a−1,b−1,φ−1L)+c−1sin⁡φ−1L=E(a−1,b−1,φ−1L)−a−12F(a−1,b−1,φ−1L)−2a−1c−1F(a,b,φ)+c−1sin⁡φ−1L=E(a−2,b−2,φ−2L)−a−22F(a−2,b−2,φ−2L)−2a−1c−1F(a,b,φ)−4a−2c−2F(a,b,φ)+c−1sin⁡φ−1L+c−2sin⁡φ−2L=AGM⁡(a,b)φ−∞L−AGM⁡(a,b)φ−∞L−∑n=1∞2na−nc−nF(a,b,φ)+∑n=1∞c−nsin⁡φ−nL=−∑n=1∞2na−nc−nF(a,b,φ)+∑n=1∞c−nsin⁡φ−nL
E(a,b,φ)=[a2−∑n=1∞2na−nc−n]F(a,b,φ)+∑n=1∞c−nsin⁡φ−nL
E(φ,k)=[1−∑n=1∞2na−nc−na2]F(φ,k)+∑n=1∞c−nasin⁡φ−nL=[1−2a−1c−1a2−2a−1c−1a22a−2c−2a−1c−1−2a−1c−1a22a−2c−2a−1c−12a−3c−3a−2c−2−⋯]F(φ,k)+c−1asin⁡φ−1L+a−1ac−2a−1sin⁡φ−2L+a−1aa−2a−1c−3a−2sin⁡φ−3L+⋯=[1−k22(1+k−12+k−1k−24+k−1k−2k−38+⋯)]F(φ,k)+k−11+k−1sin⁡φ−1L+11+k−1k−21+k−2sin⁡φ−2L+11+k−111+k−2k−31+k−3sin⁡φ−3L+⋯=[1−k(∑n=0∞∏i=0nk−i2n)]F(φ,k)+∑n=1∞k−n∏i=1n(1+k−i)sin⁡φ−nL=[1−k(k2+kk−14+kk−1k−28+kk−1k−2k−316+⋯)]F(φ,k)+k−1k2k−1sin⁡φ−1L+k−2k2k−1k−12k−2sin⁡φ−2L+k−3k2k−1k−12k−2k−22k−3sin⁡φ−3L+⋯=[1−k(k2+kk−14+kk−1k−28+kk−1k−2k−316+⋯)]F(φ,k)+k[k−12sin⁡φ−1L+k−1k−24sin⁡φ−2L+k−1k−2k−38sin⁡φ−3L+⋯]=[1−k(∑n=0∞∏i=0nk−i2n)]F(φ,k)+k∑n=1∞∏i=1nk−i2nsin⁡φ−nL

L'intégrale elliptique complète de deuxième espèce est :

E(k)=[1−k(∑n=0∞∏i=0nk−i2n)]K(k)

Moyenne arithmético-géométrique

On remarque une formule de Gauss :

F(a,b,π/2)=12F(a−1,b−1,π)=F(a−1,b−1,π/2)=F(a+b2,ab,π/2)

Donc, on a :

F(a,b,π/2)=π2AGM⁡(a,b)=1aK(1−(ba)2)=2a+bK(a−ba+b)⇔AGM⁡(a,b)=πa2K(1−(ba)2)=π(a+b)4K(a−ba+b)

En opérant a→1 et b→k, on a :

K(k)=π2AGM⁡(1,k′)=π2AGM⁡(1+k,1−k)

Transformation de Gauss

On a aussi une transformation de Gauss des intégrales elliptiques de deuxième espèce.

E(φn+1G,kn+1)=(1+kn′+1′)E(φnG,kn)−kn′+1′F(φn+1G,kn+1)+(1−1−kn+12sin2φn+1)cot⁡φn+1

ce qui donne :

E(φnG,kn)=(1+kn′)E(φn−1G,kn−1)−kn′F(φnG,kn)+(1−1−kn2sin2φn)cot⁡φn

Transformations de Landen des intégrales elliptiques de troisième espèce

On a aussi une transformation de Landen et de Gauss des intégrales elliptiques de troisième espèce,.

Transformation de Landen

Transformation de Gauss

Transformation de Landen des fonctions elliptiques jacobiennes

Le module elliptique initial est 0<k<1. La transformée de Landen croissante augmente le paramètre et la transformée de Landen décroissante diminue le paramètre. En répétant la transformation de Landen croissante, le paramètre converge vers 1 et la fonction elliptique se rapproche d'une fonction hyperbolique. En répétant la transformation de Landen décroissante, le paramètre converge vers 0 et la fonction elliptique se rapproche d'une fonction trigonométrique. Cette propriété rend la transformée de Landen utile pour les calculs numériques de fonctions elliptiques. On utilise ici la transformation de Landen pour la transformation croissante des fonctions elliptiques jacobiennes et la transformation de Gauss pour leur transformation décroissante.

Démonstration de la transformation croissante

sn⁡(un,kn)=21+knsn⁡(un+1,kn+1)cn⁡(un+1,kn+1)1−(2kn1+kn)2sn2(un+1,kn+1)=21+knsn⁡(1+kn2un,2kn1+kn)cn⁡(1+kn2un,2kn1+kn)dn⁡(1+kn2un,2kn1+kn)
cn⁡(un,kn)=1−sn2(un,kn)=1−21+knsn2(1+kn2un,2kn1+kn)dn⁡(1+kn2un,2kn1+kn)=dn2(1+kn2un,2kn1+kn)−1−kn1+kn2kn1+kndn⁡(1+kn2un,2kn1+kn)
dn⁡(un,kn)=1−kn2sn2(un,kn)=1−2kn1+knsn2(1+kn2un,2kn1+kn)dn⁡(1+kn2un,2kn1+kn)=dn2(1+kn2un,2kn1+kn)+1−kn1+kn21+kndn⁡(1+kn2un,2kn1+kn)

Démonstration de la transformation décroissante

En opérant les transformations un→un1+kn et kn→1−kn′1+kn′ au sein de sn⁡[(1+kn)un,2kn1+kn]=sn⁡(un+1,kn+1)=(1+kn)sn⁡(un,kn)1+knsn2(un,kn), on a :

sn⁡(un,kn)=21+kn′sn⁡(1+kn′2un,1−kn′1+kn′)1+1−kn′1+kn′sn2(1+kn′2un,1−kn′1+kn′)
cn⁡(un,kn)=1−sn2(un,kn)=cn⁡(1+kn′2un,1−kn′1+kn′)dn⁡(1+kn′2un,1−kn′1+kn′)1+1−kn′1+kn′sn2(1+kn′2un,1−kn′1+kn′)
dn⁡(un,kn)=1−kn2sn2(un,kn)=1−1−kn′1+kn′sn2(1+kn′2un,1−kn′1+kn′)1+1−kn′1+kn′sn2(1+kn′2un,1−kn′1+kn′)=dn2(1+kn′2un,1−kn′1+kn′)−2kn′1+kn′21+kn′−dn2(1+kn′2un,1−kn′1+kn′)

Conversion de nombres imaginaires

La transformation de Landen croissante et la transformation de Landen décroissante sont alternées par la transformation imaginaire.

sn⁡(iu,k′)=isc⁡(u,k)=isn⁡(u,k)cn⁡(u,k)

En utilisant la transformation de Landen croissante, on a :

isn⁡(u,k)cn⁡(u,k)=2i1+ksn⁡(1+k2u,2k1+k)cn⁡(1+k2u,2k1+k)dn⁡(1+k2u,2k1+k)dn2(1+k2u,2k1+k)−1−k1+k2k1+kdn⁡(1+k2u,2k1+k)=4ik(1+k)2sn⁡(1+k2u,2k1+k)cn⁡(1+k2u,2k1+k)dn2(1+k2u,2k1+k)−1−k1+k

En utilisant la transformation imaginaire, on a :

sn⁡(iu,k′)=4k(1+k)2sc⁡(1+k2iu,1−k1+k)nc⁡(1+k2iu,1−k1+k)dc2(1+k2iu,1−k1+k)−1−k1+k=4k(1+k)2sn⁡(1+k2iu,1−k1+k)dn2(1+k2iu,1−k1+k)−1−k1+kcn2(1+k2iu,1−k1+k)=21+ksn⁡(1+k2iu,1−k1+k)1+1−k1+ksn2(1+k2iu,1−k1+k)

En opérant iu→u et k′→k, on a :

sn⁡(u,k)=21+k′sn⁡(1+k′2u,1−k′1+k′)1+1−k′1+k′sn2(1+k′2u,1−k′1+k′)

Il s'agit d'une transformation de Landen décroissante.

Transformations de degré supérieur

Les transformations de modules suivantes peuvent être effectuées à l'aide des Amplitudinis Sine et des Amplitudinis Delta :

K(k)=3sn⁡[13K(k),k]2−sn⁡[13K(k),k]K(k3sn⁡[13K(k),k]4)=51+2dn⁡[25K(k),k]+2dn⁡[45K(k),k]K(k5sn⁡[15K(k),k]4sn⁡[35K(k),k]4)=71+2dn⁡[27K(k),k]+2dn⁡[47K(k),k]+2dn⁡[67K(k),k]K(k7sn⁡[17K(k),k]4sn⁡[37K(k),k]4sn⁡[57K(k),k]4)

Ici, la fonction sinus de Jacobi sn⁡[13K(k),k] est solution de l'équation k2x4−2k2x3+2x−1=0.

Globalement, la formule suivante s'applique à toutes les valeurs n∈ℕ et 0⩽k⩽1 :

K(k)=n∑a=1ndn⁡[2anK(k),k]K(kn∏a=1nsn⁡[2a−1nK(k),k]2)
F(arctan⁡x,k)=1∑a=1ndn⁡[2anK(k),k]F(∑a=1narctan⁡{dn⁡[2anK(k),k]x},kn∏a=1nsn⁡[2a−1nK(k),k]2)

Deux exemples de la transformation du troisième degré susmentionnée sont présentés ci-dessous :

(368)2(23)4−2(368)2(23)3+2(23)−1=0⇔x=sn⁡[13K(368),368]=23

K(368)=3sn⁡[13K(368),368]2−sn⁡[13K(368),368]K{k3sn⁡[13K(368),368]4}=32K[(368)3(23)4]=32K(616)

(51016)2(45)4−2(51016)2(45)3+2(45)−1=0⇔x=sn⁡[13K(51016),(51016)]=45

K(51016)=3sn⁡[13K(51016),51016]2−sn⁡[13K(51016),51016]K{k3sn⁡[13K(51016),51016]4}=2K[(51016)3(45)4]=2K(108)

En les comparant selon un motif en forme de X, des modules elliptiques complémentaires de Pythagore deviennent visibles des deux côtés de l'échelle de l'équation :

K(368)=32K(616)

K(51016)=2K(108)

Les modules de même couleur sont complémentaires les uns des autres en termes pythagoriciens.

Récapitulatif

Modules parent et enfant

Modules parent et enfant
kn+1=2kn1+kn=21−kn'241+1−kn'2 kn′+1′=1−kn1+kn=1−1−kn'21+1−kn'2
kn−1=1−1−kn21+1−kn2=1−kn′1+kn′ kn′−1′=21−kn241+1−kn2=2kn′1+kn′

Intégrales elliptiques de première espèce

Transformations de Landen des intégrales elliptiques incomplètes de première espèce
Transformation croissante Transformation décroissante
Landen F(φ,k)=21+kF(arcsin⁡1+ksin2φ−cos⁡φ1−k2sin2φ2,2k1+k) F(φ,k)=11+k′F(arccos⁡1−(1+k′)sin2φ1−k2sin2φ,1−k′1+k′)
Gauss F(φ,k)=11+kF(arcsin⁡(1+k)sin⁡φ1+ksin2φ,2k1+k) F(φ,k)=21+k′F(arcsin⁡(1+k′)sin⁡φ1+1−k2sin2φ,1−k′1+k′)
F(φ,k)=ln⁡tan⁡(φ∞L2+π4)∏n=0∞21+kn=ln⁡tan⁡(φ∞L2+π4)∏n=0∞knk=limn→∞φ−nL∏i=1n1+k−i2=limn→∞φ−nLk−nk∏i=1nk−i
F(φ,k)=limn→∞ln⁡tan⁡(φnG2+π4)∏i=0n11+ki=limn→∞ln⁡tan⁡(φnG2+π4)∏i=0nkn2nk=φ−∞G∏n=1∞(1+k−n)=φ−∞Glimn→∞2nk−nk∏i=1nk−i
Transformations de Landen des intégrales elliptiques complètes de première espèce
Transformation croissante Transformation décroissante
Landen K(k)=21+kF(arcsin⁡1+k2,2k1+k) K(k)=21+k′K(1−k′1+k′)
Gauss K(k)=11+kK(2k1+k) K(k)=21+k′K(1−k′1+k′)


Représentation par produit infini
Transformation décroissante
Gauss K(k)=π2∏n=0∞21+k−′n′=π2∏n=1∞(1+k−n)








Quotient K′/K
K(k'n)K(kn)=2K(kn′+1′)K(kn+1)

Intégrales elliptiques de deuxième espèce

Transformations de Landen des intégrales elliptiques incomplètes de deuxième espèce
Transformation décroissante
Landen E(φ,k)=1+k′2E(arccos⁡1−(1+k′)sin2φ1−k2sin2φ,1−k′1+k′)−k′1+k′F(arccos⁡1−(1+k′)sin2φ1−k2sin2φ,1−k′1+k′)+k2sin⁡(2φ)41−k2sin2φ
E(φ,k)=[1−k22(∑n=0∞∏i=1nk−i2n)]F(φ,k)+k∑n=1∞∏i=1nk−i2nsin⁡φ−nL
Transformations de Landen des intégrales elliptiques complètes de deuxième espèce
Transformation décroissante
Landen E(k)=(1+k′)E(1−k′1+k′)−2k′1+k′K(1−k′1+k′)

Fonctions elliptiques jacobiennes de base

Transformations de Landen d'une fonction elliptique
Transformation croissante (Landen) Transformation décroissante (Gauss)
sn⁡(u,k)=21+ksn⁡(1+k2u,2k1+k)cn⁡(1+k2u,2k1+k)dn⁡(1+k2u,2k1+k) sn⁡(u,k)=21+k′sn⁡(1+k′2u,1−k′1+k′)1+1−k′1+k′sn2(1+k′2u,1−k′1+k′)
cn⁡(u,k)=dn2(1+k2u,2k1+k)−1−k1+k2k1+kdn⁡(1+k2u,2k1+k) cn⁡(u,k)=cn⁡(1+k′2u,1−k′1+k′)dn⁡(1+k′2u,1−k′1+k′)1+1−k′1+k′sn2(1+k′2u,1−k′1+k′)
dn⁡(u,k)=dn2(1+k2u,2k1+k)+1−k1+k21+kdn⁡(1+k2u,2k1+k) dn⁡(u,k)=1−k′1+k′−[1−dn2(1+k′2u,1−k′1+k′)]1−k′1+k′+[1−dn2(1+k′2u,1−k′1+k′)]

Références

A : Elliptic Integrals - Harris Hancock, John Wiley & Sons, New York (1917)

Autres références :

  1. ↑ An investigation of a general theorem for finding the length of any arc of any conic hyperbola, by means of two elliptic arcs with some other new and useful theorems deduced therefrom DOI 10.1098/rstl.1775.0028
  2. ↑ C. F. Gauss et Nachlass, « Arithmetisch geometrisches Mittel, Werke, Bd. 3 », Königlichen Gesell. Wiss., Göttingen,‎ , p. 361–403
  3. ↑ (en) G.M. Scarpello, D. Ritelli et A. Scimone, The hyperbola rectification from Maclaurin to Landen and the Lagrange-Legendre transformation for the elliptic integrals (lire en ligne)
  4. ↑ H. McKean and V. Moll, Elliptic Curves: Function Theory, Geometry, Arithmetic; Cambridge University Press, Cambridge (1997). (ISBN 0-521-65817-9).
  5. ↑ (en) A.L. Baker, Elliptic Functions - An Elementary Textbook for Students of Mathematics, John Wiley & Sons, New York, (lire en ligne)
  6. ↑ M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover Publications, New York (1964). (ISBN 0-486-61272-4).
  7. 1 2 Ces deux équations sont un travail partiellement inédit, notamment les membres de ces équations qui contiennent limn→∞ n'apparaissent pas dans les ouvrages. Ces membres-là ne sont peut être pas si utiles.
  8. ↑ G.N. Watson, The Marquis and the Land-Agent; A Tale of the Eighteenth Century, Mathematical Gazette 17, 5-17 (1934). Foredrag for the British Mathematical Association, 1933. Gjengitt i J.J. Berggren, J. Borwein and P. Borwein, Pi: A Source Book, Springer-Verlag, New York (2004). (ISBN 978-1-4419-1915-1).
  9. ↑ Elliptic Functions: Landen's Transformation
  10. ↑ La démonstration sans an, bn et cn, mais avec kn paraît plus difficile. Les intégrales de Bulirsch paraissent avoir un avantage. De plus, à notre époque, les formes symétriques de Carlson ont supplanté les formes de Legendre.
  11. 1 2 Transformations de Landen et de Gauss des intégrales elliptiques des trois espèces (en) « Transformation de Landen », sur NIST Handbook of Mathematical Functions, Cambridge University Press, (ISBN 978-0521192255)
  12. ↑ Louis V. Osmania University et Digital Library Of India, On The Numerical, At The University Press, (lire en ligne)
  13. ↑ Pi & the AGM: A Study in Analytic Number Theory and Computational Complexity, frères Borwein

Liens externes