Loading AI tools
histoire des techniques de calcul d'une valeur approchée de π De Wikipédia, l'encyclopédie libre
Dans l'histoire des mathématiques, les approximations de la constante π ont atteint une précision de 0,04 % de la valeur réelle avant le début de notre ère (Archimède). Au Ve siècle, des mathématiciens chinois les ont améliorées jusqu'à sept décimales.
De grandes avancées supplémentaires n'ont été réalisées qu'à partir du XVe siècle (Al-Kashi). Les premiers mathématiciens modernes ont atteint une précision de 35 décimales au début du XVIIe siècle (Ludolph van Ceulen) et 126 chiffres au XIXe siècle (Jurij Vega), dépassant la précision requise pour toute application concevable en dehors des mathématiques pures.
Le record de l'approximation manuelle de π est détenu par William Shanks, qui a calculé 527 décimales correctes vers 1873. Depuis le milieu du XXe siècle, l'approximation de π est effectuée sur ordinateurs par des logiciels spécifiques.
Le , le record est établi avec cent mille milliards de décimales par Emma Haruka Iwao, travaillant sur Google Cloud durant 157 jours[1].
Les approximations les plus connues de π datant avant l'ère commune étaient exactes à deux décimales ; cela a été amélioré avec les mathématiques chinoises en particulier au milieu du premier millénaire, à une précision de sept décimales. Après cela, aucun progrès n'a été réalisé jusqu'à la fin de la période médiévale.[pas clair]
Les mathématiques babyloniennes évaluaient généralement π à 3, ce qui suffisait pour les projets architecturaux de l'époque (notamment dans la description du temple de Salomon dans la Bible hébraïque). Les Babyloniens savaient qu'il s'agissait d'une approximation, et une vieille tablette mathématique babylonienne exhumée près de Suse en 1936 (datée entre le XIXe et le XVIIe siècle avant notre ère) apporte une meilleure approximation de π, 25⁄8 = 3,125, soit proche de 0,5% de la valeur exacte[2],[3],[4],[5].
À peu près au même moment, le papyrus Rhind (daté de la Deuxième Période intermédiaire égyptienne, soit 1600 avant notre ère) donne l'approximation 256⁄81 ≈ 3,16 (précis à 0,6 %) en calculant l'aire d'un octogone régulier approximant un cercle[6].
Les calculs astronomiques dans le Shatapatha Brahmana (vers le VIe siècle avant notre ère) utilisent l'approximation fractionnaire 339/108 ≈ 3,139[7].
Au IIIe siècle avant notre ère, Archimède a prouvé les inégalités strictes 223⁄71 < π < 22⁄7, au moyen d'un 96-gone (précision de 2 × 10−4 et 4 × 10−4, respectivement). Il est aussi possible de prouver que 22 / 7 dépasse π grâce à un calcul d'intégrales élémentaire[8], mais cette technique n'existait pas à l'époque d'Archimède.
Au IIe siècle, Ptolémée a utilisé la valeur 377⁄120, la première approximation connue précise à trois décimales[9].
Le mathématicien chinois Liu Hui, en 263, a calculé que π serait compris entre 3,141024 et 3,142708, valeurs correspondant respectivement aux constructions d'un 96-gone et d'un 192-gone ; la moyenne de ces deux valeurs est 3,141866 (précision de 9 × 10−5). Il a également suggéré que 3,14 était une approximation suffisante pour des applications pratiques. Il a souvent été crédité d'un résultat plus précis π ≈ 3927/1250 = 3,1416 (précision de 2 × 10−6), bien que certains chercheurs pensent plutôt que celui-ci est dû au mathématicien chinois Zu Chongzhi (429-500)[10]. Zu Chongzhi est connu pour avoir donné l'encadrement de π entre 3,1415926 et 3,1415927, correct à sept décimales près. Il a donné deux autres approximations de π : π ≈ 22/7 et π ≈ 355/113. La dernière fraction est la meilleure approximation rationnelle possible de π en utilisant moins de cinq chiffres décimaux au numérateur et au dénominateur. Le résultat de Zu Chongzhi restera l'approximation la plus juste pendant près d'un millénaire.
À l'époque de l'empire Gupta (VIe siècle), le mathématicien Aryabhata, dans son traité astronomique Āryabhaṭīya, a calculé la valeur de π à cinq chiffres significatifs (π ≈ 62832/20000 = 3,1416)[11]. Il l'a utilisée pour calculer une approximation de la circonférence de la Terre[12]. Aryabhata a affirmé que son résultat donnait « approximativement » (āsanna) la circonférence d'un cercle. Nilakantha Somayaji (école du Kerala) a fait remarquer, au XVe siècle, que le mot signifie non seulement que c'est une approximation, mais que la valeur est « incommensurable » (c.-à-d. irrationnelle)[13].
Au Ve siècle, π était connu à environ sept décimales dans les mathématiques chinoises, et à environ cinq dans les mathématiques indiennes. Les progrès suivants n'ont été réalisés que plus d'un millénaire après, au XIVe siècle, lorsque le mathématicien et astronome indien Madhava de Sangamagrama, fondateur de l'école du Kerala d'astronomie et de mathématiques, a découvert une série infinie pour π, maintenant connue sous le nom de formule de Madhava-Leibniz[14],[15], et a donné deux méthodes pour calculer la valeur de π. Ce faisant, il a obtenu la série infinie
et a utilisé les 21 premiers termes pour calculer les 11 premières décimales de π.
L'autre méthode qu'il utilisait était d'ajouter un terme restant à la série originale de π. Il a utilisé le terme restant
dans l'extension de série infinie de π⁄4 pour améliorer l'approximation de π à 13 décimales, lorsque n = 75.
En 1424, Al-Kashi, un astronome et mathématicien perse, a calculé correctement 2π à 9 chiffres sexagésimaux[16]. Convertie en base dix, cette approximation est donnée par 17 chiffres :
ce qui équivaut à
Il a atteint ce niveau de précision en calculant le périmètre d'un polygone régulier à 3 × 228 côtés[17].
Dans la seconde moitié du XVIe siècle, le mathématicien français François Viète a découvert un produit infini qui converge vers π, connu sous le nom de formule de Viète.
Le mathématicien germano-néerlandais Ludolph van Ceulen (vers 1600) a calculé les 35 premières décimales de π à l'aide d'un 262-gone. Il était si fier de ce résultat qu'il l'a fait inscrire sur sa pierre tombale.
Dans Cyclometricus (1621), Willebrord Snell affirme que le périmètre du polygone inscrit approche la circonférence deux fois plus vite que le périmètre du polygone circonscrit correspondant. La preuve sera apportée par Christiaan Huygens en 1654. Snell a obtenu 7 décimales de π avec un 96-gone (en)[18].
En 1789, le mathématicien slovène Jurij Vega a calculé les 140 premières décimales de π, cependant les 14 dernières étaient incorrectes[19] ; ce résultat a néanmoins été le record du monde pendant 52 ans, jusqu'en 1841, lorsque William Rutherford (en) a calculé 208 décimales dont les 152 premières étaient correctes. Vega a ensuite amélioré la formule de John Machin à partir de 1706 et sa méthode est encore mentionnée aujourd'hui.
Le mathématicien amateur anglais William Shanks a passé plus de 20 ans à calculer les 707 premières décimales de π. Lorsqu'il eut terminé sa tâche en 1873, les 527 premières décimales étaient correctes. Il calculait de nouvelles décimales toute la matinée et passait l'après-midi à vérifier le travail de la matinée. Il s'agissait de l'approximation la plus précise de π jusqu'à l'avènement de l'ordinateur, un siècle plus tard.
En 1910, le mathématicien indien Srinivasa Ramanujan a trouvé plusieurs séries infinies de π à convergence rapide, comme
qui calcule à chaque nouvelle itération huit décimales supplémentaires de π. Ses séries sont maintenant la base des algorithmes les plus rapides actuellement utilisés pour calculer π.
À partir du milieu du XXe siècle, tous les calculs de π ont été effectués à l'aide de calculatrices ou d'ordinateurs.
En 1944, D. F. Ferguson, à l'aide d'une calculatrice mécanique, a constaté que William Shanks avait commis une erreur à la 528e décimale et que toutes les décimales suivantes étaient incorrectes.
Dans les premières années de l'ordinateur, un développement de π à 100 000 décimales[20] a été calculé par le mathématicien de l'université du Maryland Daniel Shanks (n'ayant aucun rapport avec William Shanks mentionné ci-dessus) et son équipe au Naval Research Laboratory à Washington. En 1961, Shanks et son équipe ont utilisé deux séries différentes pour calculer les décimales de π. Pour l'une, on savait que toute erreur produirait une valeur légèrement trop élevée, et pour l'autre, on savait que toute erreur produirait une valeur légèrement trop faible par rapport à la valeur réelle ; par conséquent, tant que les deux séries produisaient les mêmes chiffres, il n'y avait pratiquement aucun doute sur la justesse des décimales calculées. Les premiers 100 265 chiffres après la virgule de π ont été publiés en 1962[20].
En 1989, les frères Chudnovsky ont calculé correctement π à 4,8 × 108 décimales grâce au superordinateur IBM 3090 en utilisant la variation suivante de la série infinie de Ramanujan :
En 1999, Yasumasa Kanada et son équipe de l'université de Tokyo ont calculé π à 2,061 584 3 × 1011 décimales, grâce au superordinateur Hitachi SR8000/MPP en utilisant une autre variante de la série infinie de Ramanujan.
Certains records depuis 2009 ont été accomplis par des calculs sur des ordinateurs personnels en utilisant l'algorithme de Chudnovski
Le 6 décembre 2002, Yasumasa Kanada et une équipe de 9 autres personnes, utilisant le Hitachi SR8000, un supercalculateur avec 1 téraoctet de mémoire, ont calculé π à environ 1,241 1 × 1012 décimales en environ 600 heures[21],[22].
Le 29 avril 2009, Daisuke Takahashi et al., de l'université de Tsukuba, avec un supercalculateur T2K, ont plus que doublé le record précédent, en calculant π à plus de 2,57 × 1012 décimales, en moins de 74 heures au total (calcul et vérification)[23].
Le 31 décembre 2009, Fabrice Bellard a utilisé un ordinateur personnel pour calculer un peu moins de 2,7 × 1012 décimales de π. Les calculs ont été effectués en binaire, puis le résultat a été converti en base 10. Les étapes de calcul, de conversion et de vérification ont pris au total 131 jours[24].
De 2010 à 2020, tous les records ont été obtenus à l'aide du logiciel « y-cruncher » d'Alexander Yee :
Archimède, dans De la mesure du cercle, a créé le premier algorithme pour le calcul de π, basé sur l'idée que le périmètre de n'importe quel polygone inscrit dans un cercle est inférieur à la circonférence du cercle qui, à son tour, est inférieur au périmètre de tout polygone circonscrit à ce cercle. Il a commencé par des hexagones réguliers inscrits et circonscrits, dont les périmètres sont facilement déterminés. Il montre alors comment calculer les périmètres des polygones réguliers ayant deux fois plus de côtés qui sont inscrits et circonscrits autour du même cercle. Il s'agit d'une procédure récursive qui serait décrite aujourd'hui comme suit :
Soit pk et Pk les périmètres de polygones réguliers à k côtés, inscrits et circonscrits autour du même cercle, respectivement. Alors,
Archimède a utilisé le processus décrit pour calculer successivement P12, p12, P24, p24, P48, p48, P96 et p96[39]. Il a obtenu l'encadrement suivant en partant des trois dernières valeurs
On ne sait pas pourquoi Archimède s'est arrêté à un polygone de 96 côtés. Héron rapporte dans son Metrica qu'Archimède a continué le calcul dans un livre maintenant perdu, mais lui attribue une valeur incorrecte[40].
Archimède n'a pas utilisé la trigonométrie dans ce calcul et la difficulté d'appliquer la méthode réside dans l'obtention de bonnes approximations pour les racines carrées des calculs impliqués. La trigonométrie a probablement été utilisée par Claudius Ptolémée pour obtenir la valeur de π donnée dans l'Almageste (environ 150 après J.C.)[41].
La formule de Viète, publiée par François Viète en 1593, a été obtenue en utilisant une méthode polygonale, avec des aires (plutôt que des périmètres) de polygones dont les nombres de côtés sont des puissances de 2[42].
Une amélioration trigonométrique par Willebrord Snell (1621) a permis d'obtenir de meilleures limites à partir de la méthode du polygone. Ainsi, des résultats plus précis ont été obtenus, nécessitant moins de côtés[43]. La dernière tentative majeure de calculer π par cette méthode a été effectuée en 1630 par Christopher Grienberger, qui a calculé 39 décimales en utilisant l'amélioration de Snell.
Pour des calculs rapides, on peut utiliser des formules similaires à celle de John Machin :
avec le développement en série de Taylor de la fonction arc tangente. Cette formule se vérifie plus facilement à l'aide de coordonnées polaires de nombres complexes, produisant :
((x, y) = (239, 132) est une solution de l'équation de Pell x2 − 2y2 = −1.)
Les formules de ce type sont connues sous le nom de formules du type de Machin. Elles ont bien été utilisées dans l'ère informatique pour calculer un nombre record de décimales de π mais plus récemment, d'autres formules similaires ont été utilisées.
Par exemple, Shanks et son équipe ont utilisé la formule suivante de Machin en 1961 pour calculer les 100 000 premiers chiffres de π :
ainsi que
servant de vérification.
Le record enregistré en décembre 2002 par Yasumasa Kanada (voir supra) utilisait les formules suivantes de Machin :
K. Takano (1982).
F. C. W. Störmer (1896).
D'autres formules qui ont été utilisées pour calculer les estimations de π sont entre autres :
Liu Hui (voir aussi Formule de Viète) :
Madhava :
Euler :
Transformation de Newton/Euler[44] :
où (2k + 1)!! désigne le produit des entiers impairs jusqu'à 2k + 1.
Le travail de Ramanujan est à la base de l'algorithme de Chudnovsky, les algorithmes les plus rapides utilisés pour calculer π.
Les développements décimaux extrêmement longs de π sont généralement calculés avec des algorithmes tels que l'algorithme de Gauss-Legendre et l'algorithme de Borwein (en). Ce dernier, trouvé en 1985 par Jonathan et Peter Borwein, converge très rapidement :
Pour et
où , la suite converge de manière quartique vers π et donne 171 décimales en trois étapes[45],[46], donc plus de 1012 décimales en 20 étapes. Cependant, l'algorithme de Chudnovsky (qui converge linéairement) est plus rapide[réf. nécessaire] que ces formules itératives.
Le premier million de décimales de π et 1/π sont disponibles grâce au projet Gutenberg (voir les liens externes ci-dessous).
Historiquement, la base 60 a été utilisée pour les calculs de π. Dans cette base, π peut être approché à huit chiffres (décimaux) significatifs avec le nombre 3:8:29:4460, c'est-à-dire :
(Le chiffre sexagésimal suivant est 0, ce qui permet de tronquer ici pour obtenir une approximation relativement bonne.)
En outre, les expressions suivantes peuvent être utilisées pour estimer π :
Provenant de Ramanujan, celui-ci a affirmé que la déesse de Namagiri lui était apparue dans un rêve et lui avait dit la vraie valeur de π[51].
Cette approximation curieuse résulte de l'observation que la 193e puissance de 1/π, produisant la suite 1122211125[52]...[précision nécessaire]
π peut être obtenu à partir d'un cercle, si son rayon et son aire sont connus, en utilisant la relation :
Si on considère un cercle de rayon r avec son centre situé en (0, 0), tout point dont la distance à l'origine est inférieure à r sera situé dans le disque. Le théorème de Pythagore donne la distance de n'importe quel point (x, y) au centre :
Les carrés dont le centre se trouve à l'intérieur ou exactement à la limite du cercle peuvent ensuite être comptés en testant si, pour chaque point (x, y),
Le nombre total de point satisfaisant cette condition se rapproche ainsi de l'aire du cercle, qui peut ensuite être utilisée pour calculer une approximation de π.
Cette formule s'écrit :
Autrement dit, on commence par choisir une valeur pour r. En considérant tous les points (x, y) dans lesquels x et y sont des entiers compris entre −r et r. Il suffit alors de diviser la somme représentant l'aire d'un cercle de rayon r, par r2 pour trouver une approximation de π. Par exemple, si r est égal à 5, alors les points considérés sont :
(−5,5) | (−4,5) | (−3,5) | (−2,5) | (−1,5) | (0,5) | (1,5) | (2,5) | (3,5) | (4,5) | (5,5) |
(−5,4) | (−4,4) | (−3,4) | (−2,4) | (−1,4) | (0,4) | (1,4) | (2,4) | (3,4) | (4,4) | (5,4) |
(−5,3) | (−4,3) | (−3,3) | (−2,3) | (−1,3) | (0,3) | (1,3) | (2,3) | (3,3) | (4,3) | (5,3) |
(−5,2) | (−4,2) | (−3,2) | (−2,2) | (−1,2) | (0,2) | (1,2) | (2,2) | (3,2) | (4,2) | (5,2) |
(−5,1) | (−4,1) | (−3,1) | (−2,1) | (−1,1) | (0,1) | (1,1) | (2,1) | (3,1) | (4,1) | (5,1) |
(−5,0) | (−4,0) | (−3,0) | (−2,0) | (−1,0) | (0,0) | (1,0) | (2,0) | (3,0) | (4,0) | (5,0) |
(−5,−1) | (−4,−1) | (−3,−1) | (−2,−1) | (−1,−1) | (0,−1) | (1,−1) | (2,−1) | (3,−1) | (4,−1) | (5,−1) |
(−5,−2) | (−4,−2) | (−3,−2) | (−2,−2) | (−1,−2) | (0,−2) | (1,−2) | (2,−2) | (3,−2) | (4,−2) | (5,−2) |
(−5,−3) | (−4,−3) | (−3,−3) | (−2,−3) | (−1,−3) | (0,−3) | (1,−3) | (2,−3) | (3,−3) | (4,−3) | (5,−3) |
(−5,−4) | (−4,−4) | (−3,−4) | (−2,−4) | (−1,−4) | (0,−4) | (1,−4) | (2,−4) | (3,−4) | (4,−4) | (5,−4) |
(−5,−5) | (−4,−5) | (−3,−5) | (−2,−5) | (−1,−5) | (0,−5) | (1,−5) | (2,−5) | (3,−5) | (4,−5) | (5,−5) |
Les 12 points (0, ±5), (±5, 0), (±3, ±4), (±4, ±3) sont exactement sur le cercle, tandis que 69 points sont à l'intérieur. Par conséquent, l'aire est 81, et π est approximé à environ 3,24 car 81/52 = 3,24. Les résultats pour certaines valeurs de r sont présentés dans le tableau ci-dessous :
r | Aire | Approximation de π |
---|---|---|
2 | 13 | 3,25 |
3 | 29 | 3,22222 |
4 | 49 | 3,0625 |
5 | 81 | 3,24 |
10 | 317 | 3,17 |
20 | 1257 | 3,1425 |
100 | 31417 | 3,1417 |
1000 | 3141549 | 3,141549 |
π a de nombreuses représentations en fraction continue généralisée générées par une règle simple[54]:
est la série entière pour arctan(x) en posant x = 1. Elle converge trop lentement pour avoir un intérêt pratique.
Sachant que arcsin(1) = π/2, le développement en série entière de arc sin peut être simplifié pour obtenir :
mais avec une convergence presque aussi lente que la précédente (10 termes supplémentaires donnent trois décimales de plus).
En notant que
nous avons
avec une convergence telle que pour cinq termes supplémentaires, on obtient au moins trois décimales de plus.
La formule Bailey-Borwein-Plouffe (BBP) pour calculer π a été découverte en 1995 par Simon Plouffe. En utilisant le système hexadécimal, la formule peut calculer n'importe quelle décimale de π sans avoir à calculer les décimales situées avant elle.
En 1996, Simon Plouffe en a déduit un algorithme pour extraire le n-ième chiffre décimal de π, en base 10, en un temps de calcul amélioré d'un facteur O(n3(log n)3). L'algorithme ne nécessite presque aucune mémoire de stockage[55].
Le temps de calcul de la formule de Plouffe a été amélioré de O(n2) par Fabrice Bellard, qui a démontré une formule alternative pour calculer π[56].
Au cours des années, plusieurs logiciels ont été écrits pour calculer π sur les ordinateurs personnels. Ils impliquent généralement une vérification et un échange de disque efficace pour faciliter des calculs extrêmement coûteux en ressources.
Seamless Wikipedia browsing. On steroids.
Every time you click a link to Wikipedia, Wiktionary or Wikiquote in your browser's search results, it will show the modern Wikiwand interface.
Wikiwand extension is a five stars, simple, with minimum permission required to keep your browsing private, safe and transparent.