Fraction continue

Ne pas confondre avec la notion de fonction continue.

2 = 1 + 1 2 + 1 2 + 1 2 + ⋯ {\displaystyle {\sqrt {2}}=1+{\cfrac {1}{2+{\cfrac {1}{2+{\cfrac {1}{2+\,\cdots }}}}}}}

Exemple de développement infini en fraction continue.

En mathématiques, une fraction continue ou fraction continue simple ou plus rarement fraction continuée^[1] est une expression de la forme :

a 0 + 1 a 1 + 1 a 2 + 1 a 3 + … {\displaystyle a_{0}+{\cfrac {1}{a_{1}+{\cfrac {1}{a_{2}+{\cfrac {1}{a_{3}+\dots }}}}}}}

comportant un nombre fini ou infini d'étages.

On montre qu'on peut « représenter » — en un sens qui sera précisé — tout nombre réel sous forme d'une fraction continue, finie ou infinie, dans laquelle a 0 {\displaystyle a_{0}} est un entier relatif et les autres a i {\displaystyle a_{i}} des entiers strictement positifs (appelés les coefficients du développement en fraction continue, ou ses quotients partiels, ou encore ses quotients incomplets).

Comme dans la notation décimale usuelle, où chaque réel est approché par des nombres décimaux de plus en plus précisément au fur et à mesure de la donnée des décimales successives, de même chaque réel est approché par des fractions étagées de la forme ci-dessus de plus en plus précisément au fur et à mesure qu'on ajoute des étages. En outre, s'il faut une infinité de décimales pour décrire exactement un nombre non décimal, il faut un développement infini en fraction continue pour décrire exactement un nombre irrationnel.

Les fractions continues sont utiles en approximation diophantienne, notamment parce qu'elles fournissent, en un certain sens, les « meilleures » approximations des réels par des rationnels. Cette propriété est à l'origine d'algorithmes pour l'approximation de racines carrées, mais aussi de démonstrations d'irrationalité voire de transcendance pour certains nombres comme π ou e. La périodicité des fractions continues des racines carrées d'entiers strictement supérieurs à 1 et sans facteur carré a des conséquences utiles pour l'étude de l'équation de Pell-Fermat.

Déjà usitées chez les mathématiciens indiens au Moyen Âge, les fractions continues sont étudiées en Europe dès le XVII^e siècle. Elles sont maintenant généralisées à d'autres expressions, appliquées aux approximations de séries entières appelées approximant de Padé, ou encore adaptées aux applications linéaires.

La notion de fraction continue est vaste et se retrouve dans de nombreuses branches des mathématiques. Les concepts associés peuvent être relativement simples comme l'algorithme d'Euclide, ou beaucoup plus subtils comme celui de fonction méromorphe^{[note 1]}.

Il est possible, dans un premier temps, de voir une fraction continue comme une suite d'entiers qui « représente » un réel. Cette situation est un peu la même que celle du système décimal qui représente π par la suite d'entiers 3, 1, 4, 1, 5, 9… Sous forme de fraction continue, la suite est 3, 7, 15, 1, 292, 1, 1… Un premier champ d'étude consiste à étudier la relation entre la suite 3, 7, 15, 1, 292, 1, 1… et celle des nombres rationnels que propose la fraction continue, en l'occurrence 3, 22/7, 333/106, etc., il permet de savoir comment passer de la première suite à la deuxième, comment la deuxième converge et répond à d'autres questions de cette nature. Tel est essentiellement l'objet de cet article.

Les fractions continues ont une relation particulière avec les racines carrées ou plus généralement les nombres, dits irrationnels quadratiques, de la forme a + b d {\displaystyle a+b{\sqrt {d}}} où a {\displaystyle a} et b {\displaystyle b} sont des nombres rationnels, b {\displaystyle b} non nul, et d > 1 {\displaystyle d>1} un entier sans facteur carré. Les fractions continues associées sont périodiques, à partir d'un certain rang, c'est-à-dire que la suite des entiers formant la fraction continue se répète à partir d'un certain rang et jusqu'à l'infini^[2]. Cette situation est à l'image des représentations décimales infinies de nombres rationnels. Ces fractions continues permettent de résoudre un célèbre problème d'arithmétique appelé équation de Pell-Fermat^[3]. Cette question fait l'objet de l'article « Fraction continue d'un irrationnel quadratique ».

À l'image du système décimal, la fraction continue offre des nombres rationnels de plus en plus approchés de leur cible. Ces approximations sont bien meilleures que celles décimales. La deuxième approximation décimale de π, égale à 31/10 possède un dénominateur relativement proche de celui de la deuxième approximation de la fraction continue 22/7, en revanche 22/7 est plus de 30 fois plus précis que 31/10. Ce type d'approche d'un nombre réel par un nombre rationnel est appelé approximation diophantienne. Les fractions continues y jouent un grand rôle. Elles ont permis de construire les premiers nombres transcendants connus : les nombres de Liouville^[4] ou de montrer que le nombre e est irrationnel^[5]. À condition de généraliser la définition d'une fraction continue, il devient possible de montrer que π est aussi irrationnel — cette approche est traitée dans l'article « Fraction continue et approximation diophantienne ». (En fait, e et π sont même transcendants, d'après le théorème d'Hermite-Lindemann.)

Une fraction continue ne concerne pas uniquement les nombres mais aussi les fonctions. On généralise encore plus les fractions continues en remplaçant les coefficients par des polynômes^[6]. Une motivation provient de l'analyse complexe, qui a pour objet l'étude des fonctions de la variable complexe à valeurs complexes, dérivables en tant que telles. L'approche classique consiste à les définir comme séries entières donc comme limites de polynômes. Une spécificité fréquente de ce type de fonction est de posséder des pôles. Si, au lieu d'approcher la fonction par des polynômes, on utilise des quotients, on construit une suite d'approximants de Padé qui ne possède pas nécessairement cette faiblesse^[7].

D'autres propriétés ont été étudiées. À la différence du système décimal, un entier apparaissant dans une fraction continue n'est en général pas borné par 9, il peut devenir arbitrairement grand. Alexandre Khintchine s'est intéressé à la moyenne, au sens de limite des moyennes géométriques de tous ces dénominateurs. Pour presque tous les nombres, cette moyenne est la même (le mot « presque » possède ici le sens technique de la théorie de la mesure) ; cette moyenne est appelée constante de Khintchine.

Il est aussi possible de construire des développements en fractions en plaçant les barres de fraction sur le numérateur et non en dessous : on obtient un développement en série de Engel :

x 0 = 1 a 0 + 1 a 0 a 1 + 1 a 0 a 1 a 2 + ⋯ + 1 a 0 a 1 … a n + ⋯ = 1 + 1 + 1 + ⋯ a 2 a 1 a 0 . {\displaystyle x_{0}={\frac {1}{a_{0}}}+{\frac {1}{a_{0}a_{1}}}+{\frac {1}{a_{0}a_{1}a_{2}}}+\cdots +{\frac {1}{a_{0}a_{1}\ldots a_{n}}}+\cdots ={\cfrac {1+{\cfrac {1+{\cfrac {1+\cdots }{a_{2}}}}{a_{1}}}}{a_{0}}}.}

L'usage des fractions continues est ancien. Aryabhata (476-550), un mathématicien indien les utilise pour résoudre des équations diophantiennes ainsi que pour approximer précisément des nombres irrationnels^{[8][réf. incomplète]}. Brahmagupta (598-668) étudie plus en profondeur l'équation maintenant dite de Pell-Fermat, en utilisant une identité remarquable. Il cherche à résoudre l'équation x² – 61y² = 1 et trouve la plus petite solution : x = 1 766 319 049 et y = 226 153 980.

Au XII^e siècle, la méthode est enrichie par Bhāskara II. Un algorithme, la méthode chakravala, analogue à celui des fractions continues, permet de résoudre le cas général^[9]. La différence la plus marquante avec la méthode européenne ultérieure est qu'il autorise les nombres négatifs dans la fraction, permettant une convergence plus rapide^[10].

L'apparition en Europe est plus tardive et italienne. Raphaël Bombelli (1526-1572) fait usage d'un ancêtre des fractions continues pour le calcul d'approximations de la racine carrée de 13^[11]. Pietro Cataldi (1548-1626) comprend que la méthode de Bombelli s'applique pour toutes les racines carrées, il l'utilise pour la valeur 18 et écrit un petit opuscule à ce sujet^[12]. Il remarque que les approximations obtenues sont alternativement supérieures et inférieures à la racine carrée cherchée.

Un progrès décisif a lieu en Angleterre. Le 3 janvier 1657, Pierre de Fermat défie les mathématiciens européens avec plusieurs questions dont l'équation déjà résolue par Brahmagupta^[13]. La réaction des anglais, piqués au vif^[14], est rapide. William Brouncker (1620-1684) trouve la relation entre l'équation et la fraction continue, ainsi qu'une méthode algorithmique équivalente à celle des indiens pour le calcul de la solution. Il produit la première fraction continue généralisée, pour le nombre 4/π^{[note 2]}. Ces résultats sont publiés par John Wallis qui en profite pour démontrer les relations de récurrence utilisées par Brouncker et Bhāskara II. Il donne le nom de fraction continue dans la phrase : « Nempe si unitati adjungatur fractio, quae denominatorem habeat continue fractum^[15] ». À cette époque, Christian Huygens (1629-1695) utilise les approximations rationnelles données par le développement en fractions continues pour déterminer le nombre de dents des engrenages d'un automate planétaire^[16].

Quelques questions théoriques sont résolues au siècle suivant. L'usage montre que l'algorithme des fractions continues permet de résoudre l'équation de Pell-Fermat en utilisant le fait que la fraction est périodique à partir d'un certain rang. Leonhard Euler (1707-1783) montre que si un nombre possède une fraction continue périodique, alors il est solution d'une équation du second degré à coefficients entiers^[17]. La réciproque, plus subtile^[18], est l'œuvre de Joseph-Louis Lagrange (1736-1813). Durant ce siècle, Jean-Henri Lambert (1728-1777) trouve une nouvelle utilité aux fractions continues. Il les utilise pour montrer l'irrationalité de π.

Cet usage devient fréquent au XIX^e siècle. Évariste Galois (1811-1832) trouve la condition nécessaire et suffisante pour qu'une fraction continue soit immédiatement périodique^{[note 3]}. Joseph Liouville utilise le développement en fraction continue pour exhiber des nombres transcendants : les nombres de Liouville^[4]. En 1873, Charles Hermite prouve la transcendance de e. Un sous-produit de sa preuve est une nouvelle démonstration de l'expression de la fraction continue simple de e trouvée par Euler^[19]. À la fin du siècle, Henri Padé (1863-1953) développe la théorie^[20] des approximants qui portent maintenant son nom et qui sont des fractions continues de polynômes. Cette technique est utilisée par Henri Poincaré (1854-1912) pour démontrer la stabilité du système solaire^[21]. Georg Cantor (1845-1918) prouve à l'aide des fractions continues que les points d'un segment et ceux situés à l'intérieur d'un carré sont en bijection^[22]. Les fonctions de cette nature sont étudiées dans le cadre de la théorie du chaos ; elles sont discontinues sur chaque point rationnel de l'intervalle [0, 1]^[23].

On commence par rappeler le déroulement de l'algorithme dû à Euclide de recherche du PGCD, en analysant l'exemple des deux nombres entiers 15 625 et 6 842. On procède à une suite de divisions euclidiennes avec reste :

15 625 = 2 × 6 842 + 1 941 , 6 842 = 3 × 1 941 + 1 019 , 1 941 = 1 × 1 019 + 922 , 1 019 = 1 × 922 + 97 , 922 = 9 × 97 + 49 , 97 = 1 × 49 + 48 , 49 = 1 × 48 + 1. {\displaystyle {\begin{matrix}15\;625&=&2\times &6\;842&+&1\;941,&\\6\;842&=&3\times &1\;941&+&1\;019,&\\1\;941&=&1\times &1\;019&+&922,&\\1\;019&=&1\times &922&+&97,&\\922&=&9\times &97&+&49,&\\97&=&1\times &49&+&48,&\\49&=&1\times &48&+&1.&\end{matrix}}}

Une autre manière d'interpréter cet algorithme consiste à approcher par étapes le quotient 15 625 / 6 842. La partie entière de ce quotient est 2, ce qui permet d'écrire :

15 625 6 842 = 2 + 1 941 6 842 . {\displaystyle {\frac {15\;625}{6\;842}}=2+{\frac {1\;941}{6\;842}}.}

Que peut-on dire de la fraction 1 941 /6 842, à part qu'elle est plus petite que 1 ? Elle est comprise entre 1/4 et 1/3, son inverse, 6 842 / 1 941, possède comme partie entière : 3 ; et plus précisément, si l'on utilise les résultats de la deuxième division euclidienne :

1 941 6 842 = 1 3 + 1 019 1 941 . {\displaystyle {\frac {1\;941}{6\;842}}={\cfrac {1}{3+{\cfrac {1\;019}{1\;941}}}}.}

Ainsi de proche en proche :

15 625 6 842 = 2 + 1 3 + 1 019 1 941 = 2 + 1 3 + 1 1 + 922 1 019 = ⋯ = 2 + 1 3 + 1 1 + 1 1 + 1 9 + 1 1 + 1 1 + 1 48 {\displaystyle {\cfrac {15\;625}{6\;842}}=2+{\cfrac {1}{3+{\frac {1\;019}{1\;941}}}}=2+{\cfrac {1}{3+{\cfrac {1}{1+{\cfrac {922}{1\;019}}}}}}=\dots =2+{\cfrac {1}{3+{\cfrac {1}{1+{\cfrac {1}{1+{\cfrac {1}{9+{\cfrac {1}{1+{\cfrac {1}{1+{\cfrac {1}{48}}}}}}}}}}}}}}}

qui est bien une fraction continue. On utilise parfois la notation suivante, plus commode :

15 625 6 842 = [ 2 , 3 , 1 , 1 , 9 , 1 , 1 , 48 ] . {\displaystyle {\frac {15\;625}{6\;842}}=[2,3,1,1,9,1,1,48].}

On peut comparer 15 625 / 6 842 à ses réduites obtenues en tronquant successivement le nombre d'étages de la fraction continue. Le tableau suivant donne les troncatures en notation fractionnelle puis décimale, et la différence entre la réduite et le nombre 15 625 / 6 842.

La suite des erreurs est décroissante en valeur absolue et de signes alternés.

Soit r = p / q {\displaystyle r=p/q} un nombre rationnel (avec p {\displaystyle p} et q {\displaystyle q} entiers et q > 0 {\displaystyle q>0} ). On cherche pour r {\displaystyle r} un développement fini en fraction continue, c'est-à-dire une écriture de r {\displaystyle r} sous la forme [ a 0 , … , a n ] {\displaystyle [a_{0},\dots ,a_{n}]} avec n {\displaystyle n} entier naturel, a 0 {\displaystyle a_{0}} entier relatif et a , … , a n {\displaystyle a_{,}\dots ,a_{n}} entiers > 0. On applique pour cela l'algorithme d'Euclide :

On pose p 0 = p {\displaystyle p_{0}=p} , p 1 = q {\displaystyle p_{1}=q} , et l'on construit les entiers a 0 {\displaystyle a_{0}} et p 2 {\displaystyle p_{2}} par division euclidienne :

p 0 = a 0 p 1 + p 2 , a 0 ∈ Z , 0 ⩽ p 2 < p 1 . {\displaystyle p_{0}=a_{0}p_{1}+p_{2},\quad a_{0}\in \mathbb {Z} ,0\leqslant p_{2}<p_{1}.}

Puis, tant que p i {\displaystyle p_{i}} n'est pas nul, on définit les entiers a i − 1 {\displaystyle a_{i-1}} et p i + 1 {\displaystyle p_{i+1}} par :

p i − 1 = a i − 1 p j + p i + 1 , {\displaystyle p_{i-1}=a_{i-1}p_{j}+p_{i+1},}

avec a i − 1 {\displaystyle a_{i-1}} entier au moins égal à 1 (pour i > 1 {\displaystyle i>1} ) et 0 ⩽ p i + 1 < p i {\displaystyle 0\leqslant p_{i+1}<p_{i}} . L'algorithme d'Euclide s'arrête. On note n {\displaystyle n} le plus grand entier pour lequel p n + 1 {\displaystyle p_{n+1}} n'est pas nul. On sait donc que p n / p n + 1 {\displaystyle p_{n}/p_{n+1}} est égal à l'entier a n {\displaystyle a_{n}} . On a alors :

p 0 p 1 = [ a 0 , a 1 , … , a n ] . {\displaystyle {\frac {p_{0}}{p_{1}}}=[a_{0},a_{1},\dots ,a_{n}].}

En effet :

p 0 / p 1 = [ a 0 , p 1 / p 2 ] = [ a 0 , a 1 , p 2 / p 3 ] = ⋯ = [ a 0 , a 1 , … , a n ] , {\displaystyle p_{0}/p_{1}=[a_{0},p_{1}/p_{2}]=[a_{0},a_{1},p_{2}/p_{3}]=\dots =[a_{0},a_{1},\dots ,a_{n}],}

ou encore :

p n / p n + 1 = [ a n ]  et pour  j = n − 1 , n − 2 , … , 0   : p j / p j + 1 = [ a j , p j + 1 / p j + 2 ] = [ a j , a j + 1 , … , a n ] . {\displaystyle p_{n}/p_{n+1}=[a_{n}]{\text{ et pour }}j=n-1,n-2,\ldots ,0~:\quad p_{j}/p_{j+1}=[a_{j},p_{j+1}/p_{j+2}]=[a_{j},a_{j+1},\dots ,a_{n}].}

Représentation géométrique^[24] —

L'algorithme d'Euclide permet de calculer une fraction continue dans le cas des nombres rationnels. Cet algorithme admet dans ce cadre une interprétation géométrique. Soit r = p/q un nombre rationnel, on considère un rectangle de longueur p et de largeur q, et on le pave par des carrés de côté q.

Si r est un entier, le pavage comporte exactement r carrés. Sinon, soit a₀ le nombre de carrés insérés dans le rectangle, ou encore, le premier terme de la fraction continue. Il reste une bande non pavée de dimension q × b₁ avec b₁ égal à p – a₀q ; on pave cette bande avec des carrés de dimension maximale, c'est-à-dire de côté b₁. Le nombre de carrés est égal au deuxième terme a₁ de la fraction continue. En réitérant la méthode, on obtient l'intégralité des coefficients a_p.

Dans l'image ci-contre, on pave le rectangle 30 × 13 par deux carrés de côtés 13. Il reste une bande de longueur 13 et de largeur 4. En termes de fraction continue, on obtient l'égalité :

30 13 = 2 + 4 13 = 2 + 1 ( 13 4 ) {\displaystyle {\frac {30}{13}}=2+{\frac {4}{13}}={\color {Red}2}+{\cfrac {1}{\left({\cfrac {13}{4}}\right)}}} .

Ensuite, on remarque qu'il est possible de remplir la bande restante de 3 carrés de côté 4 et il reste une bande de longueur 4 et de largeur 1, ce qui permet de terminer le calcul de la fraction continue :

13 = 3 × 4 + 1 ⇒ 30 13 = 2 + 1 ( 3 × 4 + 1 4 ) = 2 + 1 3 + 1 4 = [ 2 , 3 , 4 ] {\displaystyle 13=3\times 4+1\;\Rightarrow \;{\frac {30}{13}}={\color {Red}2}+{\cfrac {1}{\left({\cfrac {3\times 4+1}{4}}\right)}}={\color {Red}2}+{\cfrac {1}{{\color {Red}3}+{\cfrac {1}{\color {Red}4}}}}=[2,3,4]} .

La même construction permet de trouver le rationnel dont on connait le développement en fraction continue. Dans l'image de gauche on peut retrouver le rationnel dont le développement est [1, 1, 2, 3]. Le dernier coefficient est égal à 3, on trouve donc 3 petits carrés de côté 1, qui donnent la taille du carré suivant (3). L'avant dernier coefficient 2 indique qu'il existe deux carrés moyens de côté 3. Ces deux côtés et le petit carré donnent la taille du carré plus grand (7). Le coefficient associé est égal à 1, il n'en existe donc qu'un unique de cette nature. Le carré plus grand (7) et le carré moyen (3) donnent le côté du dernier carré (10). Les deux derniers carrés donnent la longueur totale du rectangle (17). La fraction recherchée est égale à 17/10.

Le processus s'arrête car p et q sont commensurables c'est-à-dire qu'il existe une longueur l et deux entiers a et b tels que p = la et q = lb.

Considérons maintenant un rectangle de longueur L et de largeur l. Si le quotient L/l n'est pas rationnel, c'est-à-dire si les longueurs L et l sont incommensurables, le processus ne s'arrête pas.

Tel est le cas pour la figure de droite représentant un rectangle d'or, c'est-à-dire un rectangle dont le rapport de la longueur sur la largeur est égal à φ le nombre d'or. On ne peut placer qu'un carré dans chaque bande ce qui amène à la représentation : φ = [ 1 , 1 , 1 , ⋯ ] {\displaystyle \varphi =[1,1,1,\cdots ]} .

La suite des numérateurs et celle des dénominateurs sont de Fibonacci.

On a donc montré que pour tout rationnel r {\displaystyle r} , l'algorithme d'Euclide fournit un développement fini en fraction continue de r {\displaystyle r} (réciproquement, tout nombre qui possède un développement fini en fraction continue est évidemment rationnel). Le développement [ a 0 , … , a n ] {\displaystyle [a_{0},\dots ,a_{n}]} obtenu ainsi a la particularité que si n {\displaystyle n} est non nul, alors a n > 1 {\displaystyle a_{n}>1} . On en déduit un second développement : r = [ a 0 , … , a n − 1 , 1 ] {\displaystyle r=[a_{0},\dots ,a_{n}-1,1]} . Ce sont les deux seuls^[25],[26].

Quand on adjoint, au calcul des a i {\displaystyle a_{i}} de ce développement, le calcul des numérateurs h i {\displaystyle h_{i}} et dénominateurs k i {\displaystyle k_{i}} des différentes réduites :

h 0 k 0 = [ a 0 ] = a 0 1 , h 1 k 1 = [ a 0 , a 1 ] = a 0 + 1 a 1 = a 0 a 1 + 1 a 1 , … {\displaystyle {\frac {h_{0}}{k_{0}}}=[a_{0}]={\frac {a_{0}}{1}},\quad {\frac {h_{1}}{k_{1}}}=[a_{0},a_{1}]=a_{0}+{\frac {1}{a_{1}}}={\frac {a_{0}a_{1}+1}{a_{1}}},\quad \ldots }

cet algorithme d'Euclide devient l'algorithme d'Euclide étendu^{[27][source détournée]}. Plus précisément, la suite des couples d'entiers ( u j , v j ) {\displaystyle (u_{j},v_{j})} , fournie par l'algorithme étendu appliqué à ( p , q ) {\displaystyle (p,q)} , coïncide avec la suite des ( k i , h i ) {\displaystyle (k_{i},h_{i})} , aux signes près et à un décalage près des indices : k i = ( − 1 ) i u i + 2 {\displaystyle k_{i}=(-1)^{i}u_{i+2}} et h i = ( − 1 ) i + 1 v i + 2 {\displaystyle h_{i}=(-1)^{i+1}v_{i+2}} . Pour tout i {\displaystyle i} , les entiers h i {\displaystyle h_{i}} et k i {\displaystyle k_{i}} sont donc premiers entre eux et p i + 1 = ( − 1 ) i ( q h i − 1 − p h i − 1 ) {\displaystyle p_{i+1}=(-1)^{i}(qh_{i-1}-ph_{i-1})} . En particulier : la dernière réduite, h n / k n {\displaystyle h_{n}/k_{n}} , est la fraction p / q {\displaystyle p/q} mise sous forme irréductible et l'avant-dernière correspond à la solution particulière de l'identité de Bézout fournie par l'algorithme d'Euclide étendu : PGCD ( p , q ) = p n + 1 = ( − 1 ) n ( q h n − 1 − p h n − 1 ) {\displaystyle {\text{PGCD}}(p,q)=p_{n+1}=(-1)^{n}(qh_{n-1}-ph_{n-1})} .

Démonstration

D'après le § « L'algorithme » de l'article sur l'algorithme d'Euclide étendu, les (u_j, v_j) sont définis par la relation de récurrence

u j + 1 = − a j − 1 u j + u i j 1 , v j + 1 = − a j − 1 v j + v j − 1 {\displaystyle u_{j+1}=-a_{j-1}u_{j}+u_{ij1},\quad v_{j+1}=-a_{j-1}v_{j}+v_{j-1}}

et l'initialisation

u 0 = 1 , u 1 = 0 , v 0 = 0 , v 1 = 1. {\displaystyle u_{0}=1,\quad u_{1}=0,\quad v_{0}=0,\quad v_{1}=1.}

D'après le § « Réduites d'une fraction continue » ci-dessous, les (k_i, h_i) suivent la relation de récurrence

k i = a j k i − 1 + k j i 2 , h i = a j h i − 1 + h i − 2 , {\displaystyle k_{i}=a_{j}k_{i-1}+k_{ji2},\quad h_{i}=a_{j}h_{i-1}+h_{i-2},}

même en leur affectant les valeurs supplémentaires

k − 2 = 1 , k − 1 = 0 , h − 2 = 0 , h − 1 = 1. {\displaystyle k_{-2}=1,\quad k_{-1}=0,\quad h_{-2}=0,\quad h_{-1}=1.}

Le lien annoncé entre ces deux suites de couples d'entiers s'en déduit, par une récurrence d'ordre 2.

Une remarque permet de généraliser la méthode précédente à un réel quelconque. Pour l'illustrer, appliquons-la sur le nombre π. La première étape, dans le cas d'un rationnel, était le calcul du quotient de la division euclidienne du numérateur par le dénominateur, ce qui n'a plus de sens pour un réel, en revanche le résultat était égal à la partie entière du rationnel, or la partie entière d'un réel a un sens. La partie fractionnaire, nécessairement plus petite que 1, était inversée, ce qui est encore possible ici. On obtient :

π = 3 + 1 ( 1 π − 3 ) ≈ 3 + 1 7 , 062 513 305 931 {\displaystyle \pi =3+{\cfrac {1}{\left({\cfrac {1}{\pi -3}}\right)}}\approx 3+{\frac {1}{7,062\;513\;305\;931}}}

Comme π – 3 est plus petit que 1 (c'est une partie fractionnaire) son inverse est plus grand que 1, et n'est pas entier puisque π est irrationnel. On peut donc lui appliquer la même démarche :

1 π − 3 ≈ 7 + 0 , 062 513 305 931 ≈ 7 + 1 15 , 996 594 41 et π ≈ 3 + 1 7 + 1 15 , 996 594 41 . {\displaystyle {\frac {1}{\pi -3}}\approx 7+0,062\;513\;305\;931\approx 7+{\frac {1}{15,996\;594\;41}}\quad {\text{et}}\quad \pi \approx 3+{\cfrac {1}{7+{\cfrac {1}{15,996\;594\;41}}}}.}

La nouvelle valeur, approximativement égale à 15,997, est encore un irrationnel strictement supérieur à 1, d'où la possibilité d'une nouvelle étape, puis d'une nouvelle :

π ≈ 3 + 1 7 + 1 15 + 1 1 + 0 , 003 417 ≈ 3 + 1 7 + 1 15 + 1 1 + 1 292 + 0 , 6 . {\displaystyle \pi \approx 3+{\cfrac {1}{7+{\cfrac {1}{15+{\cfrac {1}{1+0,003\;417}}}}}}\approx 3+{\cfrac {1}{7+{\cfrac {1}{15+{\cfrac {1}{1+{\cfrac {1}{292+0,6}}}}}}}}.}

Puisque que π est irrationnel, le processus ne s'arrête jamais (en imaginant que le calcul est réalisé avec une infinité de décimales). On obtient comme suite de fractions 3 puis 22/7 ≈ 3,1428 puis 333/106 ≈ 3,14150 puis 355/113 ≈ 3,1415929 et enfin 103 993 / 33 102, proche de π avec une précision meilleure que le milliardième. Une fois encore, la suite des erreurs est décroissante en valeur absolue et de signes alternés.

• On appelle développement en fraction continue ou développement en fraction continue simple (ou arithmétique^[28])^{[note 4]} toute suite, finie ou infinie, non vide, ( a p ) {\displaystyle (a_{p})} dont le premier terme a 0 {\displaystyle a_{0}} est un entier relatif et tous les termes suivants sont des entiers strictement positifs.
• L'ensemble de ses indices est soit de la forme { 0 , 1 , … , n } {\displaystyle \left\{0,1,\dots ,n\right\}} pour un certain entier naturel n {\displaystyle n} s'il s'agit d'une suite finie, soit égal à ℕ pour une suite infinie.
• Sa réduite d'indice p est le rationnel [ a 0 , a 1 , … , a p ] {\displaystyle [a_{0},a_{1},\dots ,a_{p}]} , défini par [ a 0 ] = a 0 , [ a 0 , a 1 ] = a 0 + 1 a 1 , … [ a 0 , … , a p ] = a 0 + 1 [ a 1 , … , a p ] . {\displaystyle [a_{0}]=a_{0},\quad [a_{0},a_{1}]=a_{0}+{\frac {1}{a_{1}}},\quad \dots [a_{0},\dots ,a_{p}]=a_{0}+{\frac {1}{[a_{1},\dots ,a_{p}]}}.}

Soit ( a p ) {\displaystyle (a_{p})} un développement en fraction continue. On lui associe deux suites d'entiers ( h p ) {\displaystyle (h_{p})} et ( k p ) {\displaystyle (k_{p})} , définies par récurrence par :

h − 2 = 0 , h − 1 = 1 , h p = a p h p − 1 + h p − 2 k − 2 = 1 , k − 1 = 0 , k p = a p k p − 1 + k p − 2 . {\displaystyle {\begin{aligned}h_{-2}=0,\quad &h_{-1}=1,\quad &h_{p}=a_{p}h_{p-1}+h_{p-2}\\k_{-2}=1,\quad &k_{-1}=0,\quad &k_{p}=a_{p}k_{p-1}+k_{p-2}.\end{aligned}}}

Alors, pour tout indice p {\displaystyle p}  :

[ a 0 , a 1 , … , a p ] = h p k p ( 1 ) a p = − h p k p − 2 − h p − 2 k p h p k p − 1 − h p − 1 k p ( 2 ) h p − 1 k p − h p k p − 1 = ( − 1 ) p . ( 3 ) {\displaystyle {\begin{matrix}[a_{0},a_{1},\ldots ,a_{p}]={\frac {h_{p}}{k_{p}}}&\quad &(1)\\a_{p}=-{\frac {h_{p}k_{p-2}-h_{p-2}k_{p}}{h_{p}k_{p-1}-h_{p-1}k_{p}}}&&(2)\\h_{p-1}k_{p}-h_{p}k_{p-1}=(-1)^{p}.&&(3)\end{matrix}}}

La propriété (3) montre, par application du théorème de Bézout, que les entiers h p {\displaystyle h_{p}} et k p {\displaystyle k_{p}} sont premiers entre eux.

Ces trois propriétés se démontrent directement^[29] mais sont aussi des cas particuliers de celles des fractions continues généralisées, démontrées dans l'article correspondant. On y donne également une interprétation matricielle de la définition des h p {\displaystyle h_{p}} et k p {\displaystyle k_{p}} , dont résulte immédiatement, par transposition, une propriété duale de (1)^[30] :

[ a p , a p − 1 , … , a 0 ] = h p h p − 1  (si  a 0 > 0 ) et  [ a p , a p − 1 , … , a 1 ] = k p k p − 1  (si  p > 0 ). {\displaystyle [a_{p},a_{p-1},\ldots ,a_{0}]={\frac {h_{p}}{h_{p-1}}}{\text{ (si }}a_{0}>0{\text{) et }}[a_{p},a_{p-1},\ldots ,a_{1}]={\frac {k_{p}}{k_{p-1}}}{\text{ (si }}p>0{\text{).}}}

Dans l'algorithme d'Euclide développé précédemment, l'entier a i {\displaystyle a_{i}} est le quotient de p i {\displaystyle p_{i}} dans la division euclidienne par p i + 1 {\displaystyle p_{i+1}} . C'est donc la partie entière du réel x i {\displaystyle x_{i}} égal à p i + 1 / p i {\displaystyle p_{i+1}/p_{i}} . La partie fractionnaire x i − a i {\displaystyle x_{i}-a_{i}} de x i {\displaystyle x_{i}} est p i + 2 / p i + 1 {\displaystyle p_{i+2}/p_{i+1}} , inverse du réel x i + 1 {\displaystyle x_{i+1}} .

On peut alors définir un développement en fraction continue pour tout réel x {\displaystyle x} . Le symbole ⌊ s ⌋ {\displaystyle \lfloor s\rfloor } désignant la partie entière du nombre s {\displaystyle s} , on pose :

x 0 = x , {\displaystyle x_{0}=x,}

ainsi que la définition récurrente : tant que x_j n'est pas entier,

a j = ⌊ x j ⌋ , x j = a j + 1 x j + 1 . {\displaystyle a_{j}=\lfloor x_{j}\rfloor ,\quad x_{j}=a_{j}+{\frac {1}{x_{j+1}}}.}

Si x {\displaystyle x} est rationnel, comme on l'a vu plus haut, il existe un n {\displaystyle n} tel que x n {\displaystyle x_{n}} soit entier : on pose a n = x n {\displaystyle a_{n}=x_{n}} , l'algorithme s'arrête, et les deux suites ( a i ) {\displaystyle (a_{i})} et ( x i ) {\displaystyle (x_{i})} sont finies. Si x {\displaystyle x} est irrationnel, l'algorithme ne s'arrête jamais et les deux suites sont infinies.

• La suite ( a p ) {\displaystyle (a_{p})} est appelée le développement en fraction continue du réel x.
• Le réel x p {\displaystyle x_{p}} (strictement supérieur à 1 si p > 0 {\displaystyle p>0} ) est appelé le quotient complet de x d'indice p {\displaystyle p} .
• Sa partie entière a p {\displaystyle a_{p}} est le coefficient de ce développement, ou quotient partiel, ou encore quotient incomplet de x d'indice p {\displaystyle p} .

On peut formaliser de manière plus informatique cet algorithme :

• Donnée : un nombre x réel.
• Initialisation : on assigne la valeur x à la variable X. La liste a est vide.
• Boucle : On assigne à la variable A la partie entière de X, on concatène cette valeur à la liste a. Si X est entier, l'algorithme s'arrête. Si X n'est pas entier, on assigne à X la valeur de 1 /(X – A) et on recommence au début de la boucle.

Soit, en Python

# Initialisation
X = x
a = []

while not X.is_integer():
    # Partie entière
    A = int(X)
    # Concaténation
    a.append(A)
    # Nouveau X
    X = 1 / (X - A)

Ou encore, on peut exprimer le même algorithme récursivement :

• si x {\displaystyle x} est entier, son développement est ( x ) {\displaystyle (x)}  ;
• sinon, soit a 0 {\displaystyle a_{0}} sa partie entière ; le développement de x {\displaystyle x} est : a 0 {\displaystyle a_{0}} , suivi du développement de 1 / ( x − a 0 ) {\displaystyle 1/(x-a_{0})} .

Si x {\displaystyle x} est irrationnel, deux notations sont fréquemment utilisées dans ce contexte :

x = a 0 + 1 a 1 + 1 a 2 + 1 ⋯ = [ a 0 , a 1 , a 2 , ⋯ ] . {\displaystyle x=a_{0}+{\cfrac {1}{a_{1}+{\cfrac {1}{a_{2}+{\cfrac {1}{\cdots }}}}}}=[a_{0},a_{1},a_{2},\cdots ].}

Elles seront légitimées plus loin : on verra entre autres que la suite des réduites converge vers x {\displaystyle x} .

Soient x {\displaystyle x} un réel, ( a p ) {\displaystyle (a_{p})} son développement en fraction continue, ( h p ) {\displaystyle (h_{p})} et ( k p ) {\displaystyle (k_{p})} les suites des numérateurs et dénominateurs des réduites associées à cette fraction continue, et ( x p ) {\displaystyle (x_{p})} la suite des quotients complets de x {\displaystyle x} .

Pour tout indice p {\displaystyle p} , on dispose de l'égalité :

x = [ a 0 , a 1 , … , a p − 1 , x p ] . {\displaystyle x=[a_{0},a_{1},\dots ,a_{p-1},x_{p}].}

Or la démonstration des propriétés (1) et (2) ci-dessus des réduites d'une fraction continue reste valide si l'entier a p {\displaystyle a_{p}} est remplacé par le réel x p {\displaystyle x_{p}} . On obtient donc, respectivement :

x = x p h p − 1 + h p − 2 x p k p − 1 + k p − 2 ( 1 ′ ) x p = − x k p − 2 − h p − 2 x k p − 1 − h p − 1 . ( 2 ′ ) {\displaystyle {\begin{matrix}x={\frac {x_{p}h_{p-1}+h_{p-2}}{x_{p}k_{p-1}+k_{p-2}}}&\quad &(1')\\x_{p}=-{\frac {xk_{p-2}-h_{p-2}}{xk_{p-1}-h_{p-1}}}.&&(2')\end{matrix}}}

La propriété (2’) :

• permet de calculer l'entier a_p (partie entière de x_p) en utilisant x comme seul réel, sans utiliser la suite de réels (x_i), qui peut accumuler à chaque étape des imprécisions si l'algorithme est utilisé sur l'outil informatique et conduire ainsi à des valeurs erronées à partir d'un certain rang ;
• montre que la suite des réels | k p x − h p | {\displaystyle |k_{p}x-h_{p}|} est strictement décroissante.

La différence de deux réduites successives d'une fraction continue infinie s'écrit h p + 1 k p + 1 − h p k p = ( − 1 ) p k p k p + 1 {\displaystyle {\frac {h_{p+1}}{k_{p+1}}}-{\frac {h_{p}}{k_{p}}}={\frac {(-1)^{p}}{k_{p}k_{p+1}}}} (voir supra), ce qui constitue le point de départ du théorème ci-dessous.

La loi de Gauss-Kuzmin donne la probabilité pour que le nombre k {\displaystyle k} apparaisse dans la suite des coefficients.

p ( k ) := P ( X = k ) = − log 2 ⁡ ( 1 − 1 ( 1 + k ) 2 )   . {\displaystyle p(k):=\mathbb {P} (X=k)=-\log _{2}\left(1-{\frac {1}{(1+k)^{2}}}\right)~.}

Plus le nombre est grand, plus la probabilité est asymptotiquement faible. Cette loi n'est valable que statistiquement, si on prend un nombre au hasard, pas forcément pour un nombre donné. Par exemple, pour les irrationnels quadratiques possédant un cycle, ou les rationnels, la loi est manifestement fausse. Mais ces nombres exceptionnels forment un ensemble maigre parmi tous les réels, et ne changent donc pas statistiquement la loi générale.

Par exemple, le nombre 292, qui apparait tôt dans la fraction continue de π, n'avait qu'une chance sur 59596 (environ) d'apparaitre, si on considère π comme un nombre aléatoire.

• Développements périodiques (nombres quadratiques)^[31] :
  • √2 = [ 1 , 2 , 2 , 2 , … ] = [ 1 , 2 ¯ ] {\displaystyle [1,2,2,2,\dots ]=[1,{\overline {2}}]}  ;
  • √3 = [ 1 , 1 , 2 , 1 , 2 , … ] = [ 1 , 1 , 2 ¯ ] {\displaystyle [1,1,2,1,2,\dots ]=[1,{\overline {1,2}}]}  ;
  • nombre d'or φ = [ 1 , 1 , 1 , 1 , 1 , 1 , 1 , … ] {\displaystyle \varphi =[1,1,1,1,1,1,1,\dots ]} = [ 1 ¯ ] {\displaystyle [{\overline {1}}]} , suite A000012 de l'OEIS ;
  • nombre métallique d'indice n φ n = [ n , n , n , n , n , … ] {\displaystyle \varphi _{n}=[n,n,n,n,n,\dots ]} .
• 2 + 3 = [ 3 , 6 , 1 , 5 , 7 , 1 , 1 , 4 , 1 , 38 , 43 , … ] {\displaystyle {\sqrt {2}}+{\sqrt {3}}=[3,6,1,5,7,1,1,4,1,38,43,\dots ]} , suite A089078 de l'OEIS.
• Constante délienne 2 3 = [ 1 , 3 , 1 , 5 , 1 , 1 , 4 , 1 , 1 , 8 , 1 , ⋯ ] {\displaystyle {\sqrt[{3}]{2}}=[1,3,1,5,1,1,4,1,1,8,1,\cdots ]} , suite A002945 de l'OEIS.
• e = [ 2 , 1 , 2 , 1 , 1 , 4 , 1 , 1 , 6 , 1 , 1 , 8 , … ] {\displaystyle [2,1,2,1,1,4,1,1,6,1,1,8,\dots ]} , suite A003417 de l'OEIS ; a 3 m − 1 = 2 m , a 0 = 2 ,  sinon  a n = 1 {\displaystyle a_{3m-1}=2m,a_{0}=2,{\text{ sinon }}a_{n}=1} (voir infra).
• π = [ 3 , 7 , 15 , 1 , 292 , 1 , 1 , 1 , 2 , 1 , 3 , 1 , … ] {\displaystyle [3,7,15,1,292,1,1,1,2,1,3,1,\dots ]} , suite A001203 de l'OEIS.
• Constante d'Euler γ = [ 0 , 1 , 1 , 2 , 1 , 2 , 1 , 4 , 3 , 13 , 5 , ⋯ ] {\displaystyle \gamma =[0,1,1,2,1,2,1,4,3,13,5,\cdots ]} , suite A002852 de l'OEIS
• [ 1 , 2 , 3 , 4 , 5 , 6 , … ] = ∑ n ⩾ 0 1 n ! 2 ∑ n ⩾ 0 n n ! 2 = 1,433 12742 … {\displaystyle [1,2,3,4,5,6,\dots ]={\frac {\sum _{n\geqslant 0}{\frac {1}{n!^{2}}}}{\sum _{n\geqslant 0}{\frac {n}{n!^{2}}}}}=1{,}43312742\dots }  ; décimales données par la suite A060997 de l'OEIS.
• [ 2 , 3 , 5 , 7 , 11 , 13 , 17 , 19 , … ] = 2,313 0367364 … {\displaystyle [2,3,5,7,11,13,17,19,\dots ]=2{,}3130367364\dots } , suite des nombres premiers ; décimales données par la suite A064442 de l'OEIS.

La liste des développements en fraction continue publiés dans l'OEIS se trouve ici.

Les fractions continues sont principalement considérées comme un objet arithmétique des mathématiques. Il est pourtant relativement aisé d'en faire une simple expression géométrique si l'on représente une fraction y/x par un rectangle de hauteur y et de largeur x. Alors la pente de sa diagonale est y/x. Une fraction continue s'exprime alors comme un assemblage de rectangles.

Avec cette méthode géométrique particulièrement simple il est aisé de retrouver l'approximation 3 ≈ 265 153 {\displaystyle {\sqrt {3}}\approx {\frac {265}{153}}} telle qu'elle fut exprimée par Archimède, au IIIe siècle de notre ère, dans son traité De la mesure du Cercle. De la même façon, la construction géométrique par une série de carrés alternés donne les nombres de la suite de Fibonacci avec une pente de valeur φ = 1 + 5 2 {\displaystyle \varphi ={\frac {1+{\sqrt {5}}}{2}}}

√2

√3 et approximation d'Archimède

Nombre d'or et suite de Fibonacci

Les usages des fractions continues sont nombreux. On trouvera par exemple dans Fraction continue et approximation diophantienne les démonstrations de l'irrationalité de e ou de π, dans Fraction continue d'un irrationnel quadratique un exemple de résolution d'équation de Pell-Fermat ou dans Approximant de Padé un prolongement analytique de la série entière de la fonction tangente. L'usage donné ici ne nécessite pour sa compréhension que les propriétés décrites dans cet article.

Christian Huygens souhaite construire, à l'aide d'un mécanisme de type horlogerie, un automate représentant le mouvement des planètes autour du soleil : « Ayant trouuè et fait executer depuis peu une machine automate qui represente les mouvements des Planetes dont la construction est d'une facon particuliere et assez simple a raison de son effect, au reste d'une grande utlitè a ceux qui estudient ou observent le cours des astres^[32]. » La difficulté à laquelle il est confronté est liée au rapport de la durée d'une année terrestre et de celle de Saturne. En un an, la Terre tourne de 359° 45′ 40″ 30‴ et Saturne de 12° 13′ 34″ 18‴. Le rapport est égal à 77 708 431/2 640 858. Combien faut-il de dents sur les deux engrenages supportant respectivement la Terre et Saturne ?

Huygens sait que les fractions continues offrent le meilleur compromis, ce qu'il exprime ainsi : « Or, lorsqu'on néglige à partir d'une fraction quelconque les derniers termes de la série et celles qui la suivent, et qu'on réduit les autres plus le nombre entier à un commun dénominateur, le rapport de ce dernier au numérateur sera voisin de celui du plus petit nombre donné au plus grand ; et la différence sera si faible qu'il serait impossible d'obtenir un meilleur accord avec des nombres plus petits^[33]. »

Un calcul en fraction continue montre que :

77 708 431 2 640 858 = [ 29 , 2 , 2 , 1 , 5 , 1 , 4 , 1 , 1 , 2 , 1 , 6 , 1 , 10 , 2 , 2 , 3 ] {\displaystyle {\frac {77\,708\,431}{2\,640\,858}}=[29,2,2,1,5,1,4,1,1,2,1,6,1,10,2,2,3]}

On obtient la suite de fractions : 29/1, 59/2, 147/5, 206/7, 1 177/40 ... Les deux premières solutions ne sont guère précises, dans le premier cas, à la fin d'une rotation de Saturne, la position de la terre est fausse à près d'un demi-tour, dans l'autre cas l'erreur dépasse 4°. La cinquième est techniquement difficile, elle demande la fabrication d'une roue à plus de 1 000 dents ou plusieurs roues. La quatrième offre une précision proche de 3/1 000. C'est celle que choisit Huygens.

Si la terre fait cent tours complets, sur l'automate planétaire Saturne en fait 700/206, soit trois tours et un angle de 143° 18′. Dans la réalité, Saturne a tourné de 143° 26′. Soit une erreur de 8 minutes d'angle, largement inférieure aux imprécisions mécaniques de l'horloge. Un calcul analogue montre que la fraction 147/5 donne, dans le même contexte, une erreur supérieure à un degré, pour une mise en œuvre d'une difficulté technique comparable.

Un calcul, dans la partie introductive de l'article, montre comment déterminer la fraction continue de π. Néanmoins, chaque étape est plus pénible car elle demande une précision sur la valeur initiale de plus en plus grande. Les séries entières, convergeant vers π, offrent bien une solution théorique pour le calcul de chaque coefficient de la fraction continue, mais il est calculatoirement trop inextricable pour être utilisable. Pour cette raison, il est plus simple d'obtenir une expression en fraction continue généralisée, en autorisant des numérateurs non nécessairement égaux à 1. La première fraction de ce type fut produite par Brouncker^{[note 2]} :

π 4 = 0 + 1 1 + 1 2 2 + 3 2 2 + 5 2 2 + … . {\displaystyle {\frac {\pi }{4}}=0+{\cfrac {1}{1+{\cfrac {1^{2}}{2+{\cfrac {3^{2}}{2+{\cfrac {5^{2}}{2+\dots }}}}}}}}.}

Une démonstration de cette égalité figure dans l'article « Formule de fraction continue d'Euler », par évaluation au point 1 d'une fraction continue généralisée de la fonction Arctangente. Ainsi, une fraction continue ne s'applique pas uniquement aux nombres, mais aussi à certaines fonctions.

De même, Euler a développé la fonction exponentielle en une fraction continue de fonctions d'une forme appropriée : e 1 / s = [ 1 , s − 1 , 1 , 1 , 3 s − 1 , 1 , 1 , 5 s − 1 , 1 , 1 , 7 s − 1 , … ] , {\displaystyle {\rm {e}}^{1/s}=[1,s-1,1,1,3s-1,1,1,5s-1,1,1,7s-1,\ldots ],} de manière à obtenir la fraction continue de e :

e = [ 2 , 1 , 2 , 1 , 1 , 4 , 1 , 1 , 6 , 1 , 1 , 8 , 1 , 1 , 10 , 1 , … ] . {\displaystyle {\rm {e}}=[2,1,2,1,1,4,1,1,6,1,1,8,1,1,10,1,\ldots ].}

• (en) Jonathan Borwein, Alf van der Poorten, Jeffrey Shallit et Wadim Zudilin (de), Neverending Fractions, Cambridge University Press, 2014
• (en) C. Brezinski, History of Continued Fractions and Padé Approximants, Berlin, Springer-Verlag, 1991 (DOI 10.1007/978-3-642-58169-4, lire en ligne)
• Bulletin de l'APMEP n^o 450
• Eugène Charles Catalan, « Théorie des fractions continues », Nouvelles annales de mathématiques, 1^re série, t. 8,‎ 1849, p. 154-202 (lire en ligne)
• Roger Descombes, Éléments de théorie des nombres, PUF, 1986
• Daniel Duverney, Théorie des nombres, Dunod, 2007
• G. H. Hardy et E. M. Wright (trad. de l'anglais par F. Sauvageot), Introduction à la théorie des nombres [« An Introduction to the Theory of Numbers »], Vuibert-Springer, 2007
• (en) Doug Hensley, Continued Fractions, World Scientific, 2006 (lire en ligne)
• Jean Trignan, Introduction aux problèmes d'approximation : fractions continues, différences finies, Éd. du Choix, 1994 (ISBN 978-2-909028-16-3)
• Georges Valiron, Théorie des fonctions, Masson, Paris, 1966, Notions sur les fractions continues arithmétiques p. 17-24

• Approximant de Padé
• Factorisation par fraction continue
• Fraction continue d'un irrationnel quadratique
• Fraction continue généralisée
• Fraction continue et approximation diophantienne
• Formule de fraction continue d'Euler
• Développement en série de Pierce
• Développement en série de Engel
• Mesure d'irrationalité
• Polyèdre de Klein
• Quotients partiels restreints (en)
• Concernant la convergence :
  • Constante de Khintchine
  • Constante de Lévy
  • Problème de convergence
  • Théorème de Lochs
  • Théorème de Śleszyński-Pringsheim

• Notices d'autorité :

  • LCCN
  • Japon
  • Israël
  • Tchéquie
• Notices dans des dictionnaires ou encyclopédies généralistes :

  • Britannica
  • Den Store Danske Encyklopædi
  • Gran Enciclopèdia Catalana
  • Internetowa encyklopedia PWN
• (en) Gilles Lachaud, « Continued fractions, binary quadratic forms, quadratic fields, and zeta functions », dans Algebra and topology 1988, Taejon, Korea Inst. Tech., 1988 (lire en ligne), p. 1-56
• (en) Un calculateur en ligne de fraction continue
• (de) Une illustration graphique de l'algorithme d'Euclide itéré pour le calcul d'une fraction continue (avec Cinderella).
• Charles Walter, « L1 Maths et L1 Info, Option Arithmétique, chap. 2 : Fractions continues », sur Université de Nice, 2011

• Arithmétique et théorie des nombres