Algorithme de Liu Hui pour π

La méthode de Liu Hui pour calculer l'aire d'un cercle

L'algorithme de Liu Hui pour Modèle:MathPi a été inventé vers le Modèle:S- par Liu Hui, mathématicien chinois du royaume de Cao Wei. Avant lui, le rapport de la circonférence d'un cercle à son diamètre était souvent pris expérimentalement à trois en Chine, tandis que Zhang Heng (78 – 139) le prenait à 3,1724 (rapport du cercle céleste au diamètre de la terre, Modèle:Formule) ou encore à $\sqrt{10}\approx 3,162$ . Liu Hui , non satisfait de cette dernière valeur, a fait remarquer qu'elle était trop grande. Un autre mathématicien : Wang Fan (219 – 257) a proposé Modèle:Formule . Toutes ces valeurs empiriques de Modèle:MathPi sont exactes à deux chiffres (c'est-à-dire une décimale). Liu Hui a été le premier mathématicien chinois à fournir un algorithme rigoureux pour le calcul de Modèle:MathPi pour une précision quelconque. Le calcul personnel de Liu Hui à partir d'un polygone à 96 côtés a fourni une précision de cinq chiffres : Modèle:Formule .

Liu Hui a remarqué dans son commentaire des Neuf chapitres sur l'art mathématique que le rapport de la circonférence d'un hexagone inscrit au diamètre du cercle était de trois, donc Modèle:MathPi doit être supérieur à trois. Il a ensuite fourni une description détaillée étape par étape d'un algorithme itératif pour calculer Modèle:MathPi pour toute précision voulue basée sur des polygones bissecteurs ; il a estimé Modèle:MathPi entre 3,141024 et 3,142708 à partir d'un 96-gone ; il a suggéré que 3,14 était une approximation assez bonne et a exprimé Modèle:MathPi sous la forme 157/50 ; il a admis que ce nombre était un peu petit. Plus tard, il a inventé une méthode rapide pour l'améliorer et a obtenu Modèle:Formule avec seulement un 96-gone, un niveau de précision comparable à celui d'un 1536-gone. Sa contribution la plus importante dans ce domaine a été son algorithme itératif simple de calcul de Modèle:MathPi.

Aire d'un cercle

L'aire de l'intérieur d'un cercle est égale au rayon multiplié par la moitié de la circonférence, ou Modèle:Formule.

Liu Hui a écrit :

Apparemment, Liu Hui avait déjà maîtrisé le concept de limite

$\underset{N\to \mathrm{\infty }}{\lim }\text{aire du }N\text{-gone}=\text{aire du cercle}.$

De plus, Liu Hui a prouvé que l'aire d'un cercle est la moitié de sa circonférence multipliée par son rayon. Il a dit : Modèle:Citation bloc

Dans le diagramme Modèle:Formule = rayon excédentaire. En multipliant Modèle:Formule par un côté, on obtient Modèle:Formule oblong qui dépasse la limite du cercle. Si un côté du polygone est petit (c'est-à-dire qu'il y a un très grand nombre de côtés), alors le rayon en excès sera petit, donc la zone en excès sera petite.

Comme dans le diagramme, lorsque Modèle:Formule, Modèle:Formule et Modèle:Formule .

Lorsque Modèle:Formule, la moitié du périmètre du Modèle:Formule-gone s'approche d'un demi-cercle, donc une demi-circonférence d'un cercle multipliée par son rayon est égale à l'aire du cercle. Liu Hui n'a pas expliqué en détail cette déduction. Cependant, cela va de soi en utilisant le "principe du complément entrée-sortie" de Liu Hui qu'il a fourni ailleurs dans Les Neuf chapitres sur l'art mathématique : Modèle:Citation

Ainsi, le réarrangement des six triangles verts, des trois triangles bleus et des trois triangles rouges en un rectangle de largeur vaut Modèle:Formule, et la hauteur Modèle:Formule montre que l'aire du dodécagone veut Modèle:Formule .

En général, la multiplication de la moitié de la circonférence d'un Modèle:Formule-gone par son rayon donne l'aire d'un Modèle:Formule-gone. Liu Hui a utilisé ce résultat de manière répétitive dans son algorithme pour Modèle:MathPi.

Inégalité de Liu Hui pour Modèle:MathPi

Inégalité de Liu Hui pour Modèle:MathPi

Liu Hui a prouvé un encadrement de Modèle:MathPi en considérant l'aire des polygones inscrits à Modèle:Formule et Modèle:Formule côtés.

Dans la figure, la zone jaune est celle d'un Modèle:Formule-gone, d'aire $A_N$, et la zone jaune plus la verte est celle d'un Modèle:Formule-gone, d'aire $A_{2N}$. Par conséquent, la zone verte représente la différence entre celles du Modèle:Formule-gone et du N-gone :

$D_{2N}=A_{2N}-A_N.$

La zone rouge a même aire que la verte, et est donc aussi d'aire $D_{2N}$. Donc

Zone jaune + zone verte + zone rouge = $A_{2N}+D_{2N}.$

Soit $A_C$ l'aire du disque. Alors

$A_{2N}<A_C<A_{2N}+D_{2N}.$

Si le rayon du cercle est pris égal à 1, alors nous avons l'inégalité de Liu Hui pour Modèle:MathPi :

$A_{2N}<\pi <A_{2N}+D_{2N}.$

Algorithme itératif

Algorithme de Liu Hui pour Modèle:MathPi

Liu Hui a commencé par un hexagone inscrit. Soit Modèle:Formule la longueur d'un côté Modèle:Formule de l'hexagone, Modèle:Formule est le rayon du cercle.

La droite (Modèle:Formule) est la médiatrice de [[[:Modèle:Formule]]] ; [[[:Modèle:Formule]]] devient un côté du dodécagone (12-gone), soit Modèle:Formule sa longueur . Soit Modèle:Formule la longueur de [[[:Modèle:Formule]]] et Modèle:Formule la longueur de [[[:Modèle:Formule]]].

Modèle:Formule, Modèle:Formule sont deux triangles rectangles. Liu Hui utilise plusieurs fois le théorème de Pythagore :

$G^2=r^2-{\left(\frac{M}{2}\right)}^2$

$G=\sqrt{r^2-\frac{M^2}{4}}$

$j=r-G=r-\sqrt{r^2-\frac{M^2}{4}}$

$m^2={\left(\frac{M}{2}\right)}^2+j^2$

$m=\sqrt{{\left(\frac{M}{2}\right)}^2+j^2}$

$m=\sqrt{{\left(\frac{M}{2}\right)}^2+{(r-G)}^2}$

$m=\sqrt{{\left(\frac{M}{2}\right)}^2+{(r-\sqrt{r^2-\frac{M^2}{4}})}^2}$

À partir de là, il existe une technique pour déterminer Modèle:Formule à partir de Modèle:Formule, qui donne la longueur du côté d'un polygone ayant le double d'arêtes. En partant d'un hexagone, Liu Hui a pu déterminer la longueur du côté d'un dodécagone à l'aide de cette formule. On itère pour déterminer la longueur du côté d'un tétraicosagone (24 côtés) compte tenu de la longueur du côté d'un dodécagone. On pourrait le faire récursivement autant de fois que nécessaire. Sachant déterminer l'aire de ces polygones, Liu Hui a pu alors approcher Modèle:MathPi.

Avec $r=10$ unités, il a obtenu

aire d'un 96-gone $A_{96}=313+\frac{584}{625}$

aire d'un 192-gone $A_{192}=314+\frac{64}{625}$

Différence entre un 96-gone et un 48-gone :

$D_{192}=314+\frac{64}{625}-(313+\frac{584}{625})=\frac{105}{625}$

par l'inégalité de Liu Hui pour Modèle:MathPi :

$A_{2N}<A_C<A_{2N}+D_{2N}.$

Puisque Modèle:Formule = 10, $A_C=100\times \pi$

donc:

$314+\frac{64}{625}<100\times \pi <314+\frac{64}{625}+\frac{105}{625}$

$314+\frac{64}{625}<100\pi <314+\frac{169}{625}$

$3,141024<\pi <3,142704.$

Il n'a jamais pris Modèle:MathPi comme la moyenne de la borne inférieure 3,141024 et de la borne supérieure 3,142704. Au lieu de cela, il a suggéré que 3,14 était une approximation suffisamment bonne pour Modèle:MathPi, et l'a exprimé sous la forme d'une fraction Modèle:Sfrac ; il a souligné que ce nombre est légèrement inférieur à la valeur réelle de Modèle:MathPi.

Liu Hui a effectué son calcul par calcul de la tige et a exprimé ses résultats avec des fractions. Cependant, la nature itérative de l'algorithme de Liu Hui pour Modèle:MathPi est assez claire :

$2-m^2=\sqrt{2+(2-M^2)},$

où Modèle:Formule est la longueur d'un côté du polygone d'ordre suivant coupé en deux par Modèle:Formule . Le même calcul est fait à plusieurs reprises, chaque étape ne nécessitant qu'une addition et une extraction de racine carrée.

Méthode rapide

Le calcul des racines carrées des nombres irrationnels n'était pas une tâche aisée au Modèle:S- avec les réglettes. Liu Hui a découvert un raccourci en comparant les différentiels d'aire des polygones et a découvert que la proportion de la différence d'aire des polygones d'ordre successif était d'environ 1/4.

Soit Modèle:Formule_{Modèle:Formule} la différence des aires de Modèle:Formule-gones et (Modèle:Formule)-gones

$D_N=A_N-A_{N/2}$

Il a trouvé:

$D_{96}\approx \frac{1}{4}{D}_{48}$

$D_{192}\approx \frac{1}{4}{D}_{96}$ Modèle:Refl

Ainsi:

$\begin{array}{cc}{D}_{384}& \approx \frac{1}{4}{D}_{192}\\ D_{768}& \approx {\left(\frac{1}{4}\right)}^2{D}_{192}\\ D_{1536}& \approx {\left(\frac{1}{4}\right)}^3{D}_{192}\\ D_{3072}& \approx {\left(\frac{1}{4}\right)}^4{D}_{192}\\ & ⋮\end{array}$

Aire du cercle de rayon unitaire =

$\pi =A_{192}+D_{384}+D_{768}+D_{1536}+D_{3072}+\cdots \approx A_{192}+F\cdot D_{192}.$

Dans lequel

$F=\frac{1}{4}+{\left(\frac{1}{4}\right)}^2+{\left(\frac{1}{4}\right)}^3+{\left(\frac{1}{4}\right)}^4+\cdots =\frac{\frac{1}{4}}{1-\frac{1}{4}}=\frac{1}{3}.$

C'est-à-dire que toutes les zones excédentaires ultérieures s'élèvent à un tiers de $D_{192}$

aire du cercle unitaire $=\pi \approx A_{192}+\frac{1}{3}{D}_{192}\approx \frac{3927}{1250}=3,1416.$ Modèle:Refl

Liu Hui était assez satisfait de ce résultat car il avait obtenu le même résultat avec le calcul pour un 1536-gone, obtenant la surface d'un 3072-gone. Cela répond à quatre questions :

1. Pourquoi s'est-il arrêté net à Modèle:Formule₁₉₂ dans la présentation de son algorithme ? Parce qu'il a découvert une méthode rapide pour améliorer la précision de Modèle:MathPi, obtenant le même résultat obtenu par un 1536-gone avec seulement un 96-gone. Après tout, le calcul des racines carrées n'était pas une tâche simple avec le calcul de la tige. Avec la méthode rapide, il lui suffisait d'effectuer une soustraction supplémentaire, une division supplémentaire (par 3) et une addition supplémentaire, au lieu de quatre extractions de racine carrée supplémentaires.
2. Pourquoi a-t-il préféré calculer Modèle:MathPi en calculant les aires au lieu des circonférences des polygones successifs ? Car la méthode rapide nécessitait des informations sur la différence des aires des polygones successifs.
3. Qui est le véritable auteur du paragraphe contenant le calcul de $\pi =\frac{3927}{1250}.$
4. Ce célèbre paragraphe commençait par Modèle:Citation. De nombreux érudits, parmi lesquels Yoshio Mikami et Joseph Needham, pensaient que le paragraphe "conteneur en bronze de la dynastie Han" était l'œuvre de Liu Hui et non de Zu Chongzhi comme d'autres le pensaient, en raison de la forte corrélation des deux méthodes par le calcul de la surface, et parce que il n'y avait pas un seul mot mentionnant le résultat de Zu Modèle:Nobr obtenu par un 12288-gone.

Développements ultérieurs

Liu Hui a établi un algorithme solide pour le calcul de Modèle:MathPi avec n'importe quelle précision.

• Zu Chongzhi connaissait le travail de Liu Hui et a obtenu une plus grande précision en appliquant son algorithme à un 12288-gone.

D'après la formule de Liu Hui pour un Modèle:Formule-gone :

$A_{2N}=m_N\times r$

Pour un 12288-gone inscrit dans un cercle de rayon unitaire :

$A_{24576}=3,14159261864<\pi$ .

De l'inégalité de Liu Hui pour Modèle:MathPi :

$A_{24576}<\pi <A_{24576}+D_{24576}$

Dans lequel $D_{24576}=A_{24576}-A_{12288}=0,0000001021$

$A_{24576}=3,14159261864<\pi <3,14159261864+0,0000001021$ .

Donc

$3,14159261864<\pi <3,141592706934$

Tronqué à huit chiffres significatifs :

$3,1415926<\pi <3,1415927$ .

C'est la fameuse inégalité de Zu Chongzhi pour Modèle:MathPi.

Zu Chongzhi a ensuite utilisé la formule d'interpolation de He Chengtian (何承天, 370-447) et a obtenu une fraction approximative : $\pi \approx \frac{355}{113}$ .

Cependant, cette valeur de Modèle:MathPi a disparu dans l'histoire chinoise pendant une longue période (par exemple le mathématicien de la dynastie Song Qin Jiushao a utilisé $\pi =\frac{22}{7}$ et $\pi =\sqrt{10})$), jusqu'à ce que le mathématicien de la dynastie Yuan Zhao Yuqin travaille sur une variante de l'algorithme de Liu Hui, en coupant en deux un carré inscrit et obtient à nouveau $\pi \approx \frac{355}{113}$ .

Signification de l'algorithme de Liu Hui

L'algorithme π de Liu Hui a été l'une de ses contributions les plus importantes aux mathématiques chinoises anciennes. Il était basé sur le calcul de l'aire d'un Modèle:Formule -gone, contrairement à l'algorithme d'Archimède basé sur la longueur de la circonférence du polygone. Avec cette méthode, Zu Chongzhi a obtenu le résultat à huit chiffres : 3,1415926 < Modèle:MathPi < 3,1415927, qui a détenu le record mondial de la valeur la plus précise de Modèle:MathPi pendant des siècles, jusqu'à ce que Madhava de Sangamagrama calcule 11 chiffres au Modèle:S- ou Jamshid al -Kashi calcule 16 chiffres en 1424 ; les meilleures approximations de Modèle:MathPi connues en Europe n'étaient précises qu'à 7 chiffres jusqu'à ce que Ludolph van Ceulen calcule 20 chiffres en 1596.

Voir aussi

• Méthode d'exhaustion (Modèle:-s-)
• Algorithme de Zhao Youqin pour Modèle:MathPi (Modèle:Sp-)
• Modèle:En Proof of Newton's Formula for Pi (Modèle:S-)

Remarques

Modèle:Refa La valeur correcte est de 0,2502009052

Modèle:Refa Les valeurs correctes sont :

$A_{192}=3,1410319509$

$D_{192}=0,0016817478$

$\pi \approx A_{192}+\frac{1}{3}{D}_{192}\approxeq 3,1410319509+0,0016817478/3$

$\pi \approx 3,1410319509+0,0005605826$

$\pi \approx 3,1415925335.$

La méthode rapide de Liu Hui était potentiellement capable de fournir presque le même résultat que l'algorithme classique avec un 12288-gone (3,141592516588) avec seulement un 96-gone.

BIbliographie

• Modèle:Ouvrage.
• Wu Wenjun ed, History of Chinese Mathematics Vol III (in Chinese) Modèle:Isbn