Algorithme de Liu Hui pour π

L'aire de l'intérieur d'un cercle est égale au rayon multiplié par la moitié de la circonférence, ou A = r × C/2 = r × r × π.

Liu Hui a écrit :

« Multipliez un côté d'un hexagone par le rayon (de son cercle circonscrit), puis multipliez-le par trois, pour obtenir l'aire d'un dodécagone ; si nous coupons un hexagone en dodécagone, multipliez son côté par son rayon, puis multipliez à nouveau par six, nous obtenons l'aire d'un 24-gone ; plus nous coupons finement, plus la perte par rapport à l'aire du cercle est petite, donc en continuant coupe après coupe, l'aire du polygone résultant coïncidera et ne fera plus qu'un avec le cercle ; il n'y aura pas de perte. »

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

${\displaystyle \lim _{N\to \infty }{\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 :

« Entre un polygone et un cercle, il y a un excès de rayon. Multipliez le rayon en excès par un côté du polygone. La zone résultante dépasse la limite du cercle. »

Dans le diagramme d = rayon excédentaire. En multipliant d par un côté, on obtient ABCD 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 N → ∞, d → 0 et ABCD → 0 .

« Multipliez le côté d'un polygone par son rayon, et l'aire double ; multipliez donc la moitié de la circonférence par le rayon pour obtenir l'aire du cercle. »

Lorsque N → ∞, la moitié du périmètre du N-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 : « Découpez une forme géométrique en parties, réarrangez les parties pour former une autre forme, la zone de les deux formes seront identiques. »

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

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

Inégalité de Liu Hui pour π

Liu Hui a prouvé un encadrement de π en considérant l'aire des polygones inscrits à N et 2N côtés.

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

${\displaystyle D_{2N}=A_{2N}-A_{N}.}$

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

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

Soit ${\displaystyle A_{C}}$ l'aire du disque. Alors

${\displaystyle 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 π :

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

Algorithme de Liu Hui pour π

Liu Hui a commencé par un hexagone inscrit. Soit M la longueur d'un côté AB de l'hexagone, r est le rayon du cercle.

La droite (OPC) est la médiatrice de [AB] ; [AC] devient un côté du dodécagone (12-gone), soit m sa longueur . Soit j la longueur de [PC] et G la longueur de [OP].

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

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

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

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

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

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

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

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

À partir de là, il existe une technique pour déterminer m à partir de M, 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 π.

Avec ${\displaystyle r=10}$ unités, il a obtenu

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

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

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

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

par l'inégalité de Liu Hui pour π :

${\displaystyle A_{2N}<A_{C}<A_{2N}+D_{2N}.}$

Puisque r = 10, ${\displaystyle A_{C}=100\times \pi }$

donc:

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

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

${\displaystyle {}3,141024<\pi <3,142704.}$

Il n'a jamais pris π 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 π, et l'a exprimé sous la forme d'une fraction 157/50 ; il a souligné que ce nombre est légèrement inférieur à la valeur réelle de π.

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 π est assez claire :

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

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

Le calcul des racines carrées des nombres irrationnels n'était pas une tâche aisée au III^e siècle 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[4].

Soit D_N la différence des aires de N-gones et (N/2)-gones

${\displaystyle D_{N}=A_{N}-A_{N/2}\,}$

Il a trouvé:

${\displaystyle D_{96}\approx {\tfrac {1}{4}}D_{48}}$

${\displaystyle D_{192}\approx {\tfrac {1}{4}}D_{96}}$ ^[1]

Ainsi:

${\displaystyle {\begin{aligned}D_{384}&{}\approx {\tfrac {1}{4}}D_{192}\\D_{768}&{}\approx \left({\tfrac {1}{4}}\right)^{2}D_{192}\\D_{1536}&{}\approx \left({\tfrac {1}{4}}\right)^{3}D_{192}\\D_{3072}&{}\approx \left({\tfrac {1}{4}}\right)^{4}D_{192}\\&{}\ \ \vdots \end{aligned}}}$

Aire du cercle de rayon unitaire =

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

Dans lequel

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

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

aire du cercle unitaire ${\displaystyle {}=\pi \approx A_{192}+{\frac {1}{3}}D_{192}\approx {3927 \over 1250}\approx 3,1416.\,}$ ^[2]

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 à A₁₉₂ dans la présentation de son algorithme ? Parce qu'il a découvert une méthode rapide pour améliorer la précision de π, 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 π 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 ${\displaystyle \pi ={3927 \over 1250}.}$
4. Ce célèbre paragraphe commençait par « Un conteneur en bronze de la dynastie Han dans l'entrepôt militaire de la dynastie Jin... ». 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 3,1415926 < π < 3,1415927 obtenu par un 12288-gone.

Liu Hui a établi un algorithme solide pour le calcul de π 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 2N-gone :

${\displaystyle A_{2N}=m_{N}\times r}$

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

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

De l'inégalité de Liu Hui pour π :

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

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

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

Donc

${\displaystyle 3,14159261864<\pi <3,141592706934}$

Tronqué à huit chiffres significatifs :

${\displaystyle 3,1415926<\pi <3,1415927}$ .

C'est la fameuse inégalité de Zu Chongzhi pour π.

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

Cependant, cette valeur de π a disparu dans l'histoire chinoise pendant une longue période (par exemple le mathématicien de la dynastie Song Qin Jiushao a utilisé ${\displaystyle \pi ={22 \over 7}}$ et ${\displaystyle \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 ${\displaystyle \pi \approx {355 \over 113}}$ [5].

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 N -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 < π < 3,1415927, qui a détenu le record mondial de la valeur la plus précise de π pendant des siècles[6], jusqu'à ce que Madhava de Sangamagrama calcule 11 chiffres au XIV^e siècle ou Jamshid al -Kashi calcule 16 chiffres en 1424 ; les meilleures approximations de π connues en Europe n'étaient précises qu'à 7 chiffres jusqu'à ce que Ludolph van Ceulen calcule 20 chiffres en 1596.

• Méthode d'exhaustion (V^e siècle av. J.-C.)
• Algorithme de Zhao Youqin pour π (XIII^e – XIV^e siècles)
• (en) Proof of Newton's Formula for Pi (XVII^e siècle)

1. 1

La valeur correcte est de 0,2502009052

1. 2

Les valeurs correctes sont :

${\displaystyle A_{192}=3,1410319509}$

${\displaystyle D_{192}=0,0016817478}$

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

${\displaystyle \pi \approx 3,1410319509+0,0005605826}$

${\displaystyle \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.

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Liu Hui's π algorithm » (voir la liste des auteurs).

• (en) Joseph Needham, Science and Civilization in China : Volume 3, Mathematics and the Sciences of the Heavens and the Earth, Taipei, Caves Books, Ltd, 1986.
• Wu Wenjun ed, History of Chinese Mathematics Vol III (in Chinese) (ISBN 7-303-04557-0)