Développement en cotangente continue de Lehmer

En mathématiques, le développement en cotangente continue d'un nombre réel est une écriture de ce nombre utilisant une suite de nombres entiers et la fonction cotangente. Il a été découvert par Derrick Lehmer et permet de retrouver des propriétés équivalentes à celles du développement en fraction continue, tout en permettant d'approcher le réel limite de façon plus efficace.

Définition et propriétés[modifier | modifier le code]

En 1938, Lehmer remarque que la plupart des développements de nombres se basent sur l'itération suivante :

x=f\left(x_{0},f(x_{1},f(x_{2},\ldots ))\right)

Par exemple,

$f(x,y)=x+{\frac {y}{q}}$ donne le développement en base $q$
$f(x,y)=x+{\frac {1}{y}}$ donne le développement en fraction continue
$f(x,y)={\frac {1}{x}}(1+y)$ donne le développement en série de Engel

Il considère donc la fonction

f(x,y)={\frac {xy+1}{y-x}}=\cot \left(\operatorname {arccot}(x)-\operatorname {arccot}(y)\right).

Le développement en cotangente continue d'un nombre réel $x$ est la suite de nombres entiers $(n k)$ tel que :

x=\cot \left(\sum _{k=0}^{+\infty }(-1)^{k}\operatorname {arccot}(n_{k})\right)

Développements réguliers et réduits[modifier | modifier le code]

Le développement en cotangente continue est dit régulier si :

pour tout $k$ , $n k$ est un entier positif
la suite vérifie

\forall k\in \mathbb {N} ,\ n_{k+1}\geqslant n_{k}^{2}+n_{k}+1

Dans le cas où la suite est finie, l'inégalité est stricte pour le dernier

n k

non nul.

Un développement en cotangente continue est dit réduit s'il est interrompu à partir d'un certain terme. On en déduit une approximation rationnelle du nombre, pas forcément sous sa forme réduite.

Lehmer établit aussi que tout nombre réel $x$ admet un développement en cotangente continue régulier unique et convergent, fini si et seulement si $x$ est rationnel.

Liens avec le développement en fraction continue généralisée[modifier | modifier le code]

À partir des entiers du développement en cotangente continue d'un nombre, on peut construire son développement en fraction continue généralisée. En effet, on a, puisque $f(x,y)=x+{\frac {x^{2}+1}{y-x}}$ :

\forall x,x=\cot \left(\sum _{k=0}^{+\infty }(-1)^{k}\operatorname {arccot}(n_{k})\right)=n_{0}+{\frac {n_{0}^{2}+1}{n_{1}-n_{0}+{\dfrac {n_{1}^{2}+1}{n_{2}-n_{1}+{\frac {n_{2}^{2}+1}{n_{3}-n_{2}+\ddots }}}}}}

Construction du développement[modifier | modifier le code]

Développement pour un nombre rationnel[modifier | modifier le code]

Pour $x = p / q$ rationnel, la suite du développement s'obtient à partir d'un algorithme similaire à l'algorithme d'Euclide pour $p$ et $q$ :

Initialisation

On pose $p 0 = p$ et $q 0 = q$ .

Récurrence

À l'étape $k$ :

$n k$ est le quotient de la division euclidienne de $p k$ par $q k$
$q k +1$ est le reste de la division euclidienne de $p k$ par $q k$
$p k +1 = p k n k + q k +1$

L'algorithme s'arrête quand le reste est nul.

On l'applique pour $x = 37/25$

Étape $k$	Dividende $p_{k}$	Diviseur $q_{k}$	Équation $p_{k}=q_{k}n_{k}+q_{k+1}$	Quotient $n_{k}$	Reste $q_{k+1}$	Itération $p_{k}n_{k}+q_{k+1}=p_{k+1}$
0	37	25	37 = 25 × 1 + 12	1	12	37 × 1 + 25 = 62
1	62	12	62 = 12 × 5 + 2	5	2	62 × 5 + 12 = 322
2	322	2	322 = 2 × 161 + 0	161	0	Fin de l'algorithme

Soit :

{\frac {37}{25}}=\cot(\operatorname {arccot}(1)-\operatorname {arccot}(5)+\operatorname {arccot}(161))

Développement pour un nombre irrationnel[modifier | modifier le code]

Pour $x$ irrationnel, la suite du développement s'obtient par récurrence :

Initialisation

On pose $x 0 = x$ .

Récurrence

À l'étape $k$ :

n_{k}=\lfloor x_{k}\rfloor

x_{k+1}={\frac {x_{k}\lfloor x_{k}\rfloor +1}{x_{k}-\lfloor x_{k}\rfloor }}

On l'applique pour $x = e$

Étape $k$	Racine $x_{k}$	Partie entière $n_{k}$	Itéré $x_{k+1}$
0	$e$ ≈ 2,7182818...	2	8,961055...
1	8,961055956...	8	75,62293911...
2	75,62293911...	75	8949,669392...
3	8949,669392...	8949	119646723,6...

Soit :

\mathrm {e} =\cot(\operatorname {arccot}(2)-\operatorname {arccot}(8)+\operatorname {arccot}(75)-\operatorname {arccot}(8949)+\operatorname {arccot}(119646723)-\dots )

suite A002668 de l'OEIS

D'autres suites pour des constantes classiques sont données sur le site de l'OEIS :

{\sqrt {2}}=\cot(\operatorname {arccot}(1)-\operatorname {arccot}(5)+\operatorname {arccot}(36)-\operatorname {arccot}(3406)+\operatorname {arccot}(14694817)-\dots )

suite A002666 de l'OEIS

\pi =\cot(\operatorname {arccot}(3)-\operatorname {arccot}(73)+\operatorname {arccot}(8599)-\operatorname {arccot}(400091364)+\operatorname {arccot}(371853741549033970)-\dots )

suite A002667 de l'OEIS

{\begin{aligned}\varphi &=\cot(\operatorname {arccot}(1)-\operatorname {arccot}(4)+\operatorname {arccot}(76)-\operatorname {arccot}(439204)+\operatorname {arccot}(84722519070079276)-\dots )\\&=\cot \left(\sum _{k=0}^{+\infty }(-1)^{k}\operatorname {arccot}(L_{3^{k}})\right)\end{aligned}}

suite A006267 de l'OEIS

où $\varphi$ est le nombre d'or et $L p$ désigne le nombre de Lucas d'indice $p$ .

Constante de Lehmer[modifier | modifier le code]

Dans sa comparaison entre les développements en fraction continue et ceux en cotangente continue, Lehmer cherche à déterminer les cas où le développement converge le plus lentement. Le développement en fraction continue à la convergence la plus lente correspond à la valeur :

{\frac {-1+{\sqrt {5}}}{2}}=0+{\frac {1}{1+{\frac {1}{1+\ddots }}}}=[0,1,1,\ldots ]

Le développement en fraction cotangente continue à la convergence la plus lente correspond au cas :

n_{0}=0,\,\forall k\in \mathbb {N} ,\ n_{k+1}=n_{k}^{2}+n_{k}+1

soit

\xi =\cot(\operatorname {arccot}(0)-\operatorname {arccot}(1)+\operatorname {arccot}(3)-\operatorname {arccot}(13)+\operatorname {arccot}(183)-\operatorname {arccot}(33673)+\operatorname {arccot}(1133904603)-\ldots )\approx 0,59263271.

Ce nombre est appelée constante de Lehmer. Il a affirmé que ce nombre était transcendant mais sa preuve n'est pas valide.

Applications[modifier | modifier le code]

Les développements en cotangente continue sont une alternatives aux fractions continues en approximation diophantienne.

Le développement en cotangente continue, grâce à sa convergence très rapide, permet de construire des formules du type de Machin très efficaces. En effet, en reprenant l'exemple du développement de $37/25$ vu supra :

{\frac {37}{25}}=\cot(\operatorname {arccot}(1)-\operatorname {arccot}(5)+\operatorname {arccot}(161))

permet d'écrire :

\operatorname {arccot} \left({\frac {37}{25}}\right)=\operatorname {arccot}(1)-\operatorname {arccot}(5)+\operatorname {arccot}(161)

soit

\operatorname {arccot}(1)={\frac {\pi }{4}}=\arctan \left({\frac {25}{37}}\right)+\arctan \left({\frac {1}{5}}\right)-\arctan \left({\frac {1}{161}}\right)

Références[modifier | modifier le code]

(en) Derrick H. Lehmer, « A cotangent analogue of continued fractions », Duke Mathematical Journal, vol. 4, n^o 2,‎ juin 1938, p. 323-340 (DOI 10.1215/S0012-7094-38-00424-7, lire en ligne)
Tanguy Rivoal, « Propriétés diophantiennes du développement en cotangente continue de Lehmer », Monatshefte für Mathematik, vol. 150,‎ 2007, p. 49–71 (DOI 10.1007/s00605-006-0415-7, lire en ligne)
(en) Jeffrey Shallit, « Predictable Regular Continued Cotangent Expansions », Journal of Research of the National Bureau of Standards- B. Mathematical Sciences, vol. 80B, n^o 2,‎ 1976 (lire en ligne)

Liens externes[modifier | modifier le code]

(en) Eric W. Weisstein, « Lehmer Cotangent Expansion », sur MathWorld
(en) Eric W. Weisstein, « Lehmer's Constant », sur MathWorld

Portail de l'analyse