Voir la source de Algorithme de Liu Hui pour π

{{titre mis en forme|Algorithme de Liu Hui pour {{mathPi}}}}
[[Fichier:Cutcircle2.svg|droite|vignette| La méthode de Liu Hui pour calculer l'aire d'un cercle]]
'''L'algorithme de Liu Hui pour {{mathPi}}''' a été inventé vers le {{s-|III}} par [[Liu Hui]], mathématicien chinois du royaume de [[Royaume de Wei|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 &#x2013; 139) le prenait à 3,1724 (rapport du cercle céleste au diamètre de la terre, {{Formule|92/29}}) ou encore à <math> \sqrt{10} \approx 3,162</math> . Liu Hui , non satisfait de cette dernière valeur, a fait remarquer qu'elle était trop grande. Un autre mathématicien : [[Wang Fan]] (219 &#x2013; 257) a proposé {{Formule|π ≈ 142/45 ≈ 3,156|§=π ≈ 142/45 ≈ 3,156}} <ref>{{article |lang=en |auteur=Herman C. Schepler |année=1950 |titre=The Chronology of Pi |périodique=Mathematics Magazine |volume=23 |numéro=3 |pages=165–170 |ISSN=0025-570X}}.</ref>. Toutes ces valeurs empiriques de {{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 {{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 : {{Formule|π ≈ 3,1416|§=π ≈ 3,1416}} .

Liu Hui a remarqué dans son commentaire des ''[[Les Neuf Chapitres sur l'art mathématique|Neuf chapitres sur l'art mathématique]]'' <ref>Needham, Volume 3, 66.</ref> que le rapport de la circonférence d'un hexagone inscrit au diamètre du cercle était de trois, donc {{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 {{MathPi}} pour toute précision voulue basée sur des polygones bissecteurs ; il a estimé {{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é {{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 {{Formule|π ≈ 3,1416}} 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 {{MathPi}}.

== Aire d'un cercle ==
[[Fichier:Cutcircle.svg|droite|vignette| L'aire de l'intérieur d'un cercle est égale au rayon multiplié par la moitié de la circonférence, ou {{Formule|''A'' {{=}} ''r'' × ''C''/2 {{=}} ''r'' × ''r'' × π}}.]]
Liu Hui a écrit :

{{citation bloc|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 <ref>D'abord remarqué par le mathématicien japonais [[Yoshio Mikami]]</ref>

: <math>\lim_{N \to \infty}\text{aire du }N\text{-gone} = \text{aire du cercle}. \, </math>

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 :
{{citation bloc|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 {{Formule|''d''}} = rayon excédentaire. En multipliant {{Formule|''d''}} par un côté, on obtient {{Formule|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 {{Formule|''N'' → ∞}}, {{Formule|''d'' → 0}} et {{Formule|ABCD → 0}} .

{{citation bloc|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 {{Formule|''N'' → ∞}}, la moitié du périmètre du {{Formule|''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'' : {{citation|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 {{Formule|3''L''}}, et la hauteur {{Formule|''R''}} montre que l'aire du dodécagone veut {{Formule|3''RL''}} .

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

== Inégalité de Liu Hui pour {{MathPi}} ==
[[Fichier:Liuhui_Pi_Inequality.svg|droite|vignette| Inégalité de Liu Hui pour {{MathPi}}]]
Liu Hui a prouvé un encadrement de {{MathPi}} en considérant l'aire des polygones inscrits à {{Formule|''N''}} et {{Formule|2''N''|§=2''N''}} côtés.

Dans la figure, la zone jaune est celle d'un {{Formule|''N''}}-gone, d'aire <math>A_N</math>, et la zone jaune plus la verte est celle d'un {{Formule|2''N''|§=2''N''}}-gone, d'aire <math>A_{2N}</math>. Par conséquent, la zone verte représente la différence entre celles du {{Formule|2''N''|§=2''N''}}-gone et du ''N''-gone :

: <math>D_{2N}=A_{2N}-A_N.</math>

La zone rouge a même aire que la verte, et est donc aussi d'aire <math>D_{2N}</math>. Donc
: Zone jaune + zone verte + zone rouge = <math>A_{2N} + D_{2N}.</math>

Soit <math>A_{C}</math> l'aire du disque. Alors
: <math>A_{2N} < A_{C} < A_{2N} + D_{2N}.</math>

Si le rayon du cercle est pris égal à 1, alors nous avons l'inégalité de Liu Hui pour {{MathPi}} :
: <math> A_{2N} < \pi < A_{2N} + D_{2N}.</math>

== Algorithme itératif ==
[[Fichier:Liuhui_geyuanshu.svg|droite|vignette| Algorithme de Liu Hui pour {{MathPi}} ]]
Liu Hui a commencé par un hexagone inscrit. Soit {{Formule|M|§=''M''}} la longueur d'un côté {{Formule|AB}} de l'hexagone, {{Formule|''r''}} est le rayon du cercle.

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

{{Formule|APO}}, {{Formule|APC}} sont deux triangles rectangles. Liu Hui utilise plusieurs fois le [[théorème de Pythagore]] :

: <math>{} G^2 = r^2 - \left(\tfrac{M}{2}\right)^2</math>
: <math>{} G = \sqrt{r^2- \tfrac{M^2}{4}}</math>
: <math>{} j = r - G = r - \sqrt{r^2- \tfrac{M^2}{4}}</math>
: <math>{} m^2 = \left(\tfrac{M}{2}\right)^2 + j^2</math>
: <math>{} m = \sqrt{\left(\tfrac{M}{2}\right)^2 + j^2}</math>
: <math>{} m = \sqrt{\left(\tfrac{M}{2}\right)^2 + \left(r - G\right)^2}</math>
: <math>{} m = \sqrt{\left(\tfrac{M}{2}\right)^2 + \left(r - \sqrt{r^2- \tfrac{M^2}{4}}\right)^2}</math>

À partir de là, il existe une technique pour déterminer {{Formule||§=''m''}} à partir de {{Formule|M|§=''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 {{MathPi}}.

Avec <math>r = 10</math> unités, il a obtenu

: aire d'un [[Polygone|96-gone]] <math>{}A_{96} = 313 + {584 \over 625} </math>
: aire d'un 192-gone <math>{}A_{192} = 314 + {64 \over 625} </math>
: Différence entre un  96-gone et un 48-gone :

: <math>{}D_{192} = 314+\frac{64}{625} - \left(313+\frac{584}{625}\right) = \frac{105}{625}</math>

: par l'inégalité de Liu Hui pour {{MathPi}} :

: <math> A_{2N} < A_{C} < A_{2N} + D_{2N}.</math>

: Puisque {{Formule|''r''}} = 10, <math>A_{C} = 100 \times \pi</math>
: donc:
: <math>{}314+\frac{64}{625}<100 \times \pi <314 +\frac{64}{625} +\frac{105}{625}</math>
:: <math>{}314+\frac{64}{625}<100 \pi <314+ \frac{169}{625}</math>
: <math>{} 3,141024 < \pi < 3,142704.</math>

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

Liu Hui a effectué son calcul par [[Technique de la multiplication en Chine antique|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 {{MathPi}} est assez claire :

: <math> 2-m^2 =\sqrt{2+(2-M^2)}\, ,</math>

où {{Formule|''m''}} est la longueur d'un côté du polygone d'ordre suivant coupé en deux par {{Formule|''M''}} . 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 [[Nombre irrationnel|nombres irrationnels]] n'était pas une tâche aisée au {{S-|III}} avec les [[Baguettes à calculer|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<ref>Yoshio Mikami: Ph.D. Dissertation 1932</ref>.

Soit {{Formule|D}}<sub>{{Formule|N}}</sub> la différence des aires de {{Formule|''N''}}-gones et ({{Formule|''N''/2}})-gones

: <math>D_N = A_N - A_{N/2}\,</math>

Il a trouvé:

: <math>D_{96} \approx \tfrac{1}{4} D_{48}</math>
: <math>D_{192} \approx  \tfrac{1}{4} D_{96}</math> {{refl|1}}

Ainsi:

: <math>
\begin{align}
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{align}
</math>

Aire du cercle de rayon unitaire =

: <math>{}\pi = A_{192} + D_{384} + D_{768}+D_{1536}+D_{3072} + \cdots \approx A_{192} + F \cdot D_{192}. \, </math>

Dans lequel

: <math> 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}.</math>

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

: aire du cercle unitaire <math> {}=\pi \approx A_{192} + \frac{1}{3}D_{192} \approx {3927 \over 1250} = 3,1416.\,</math> {{refl|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 :

# Pourquoi s'est-il arrêté net à {{Formule|''A''}}<sub>192</sub> dans la présentation de son algorithme ? Parce qu'il a découvert une méthode rapide pour améliorer la précision de {{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 [[Technique de la multiplication en Chine antique|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.
# Pourquoi a-t-il préféré calculer {{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.
# Qui est le véritable auteur du paragraphe contenant le calcul de <math>\pi = {3927 \over 1250}.</math>
# Ce célèbre paragraphe commençait par {{citation|Un conteneur en bronze de la dynastie Han dans l'entrepôt militaire de la [[Dynastie Jin (265-420)|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 {{nobr|3,1415926 < {{MathPi}} < 3,1415927}} obtenu par un 12288-gone.

== Développements ultérieurs ==
Liu Hui a établi un algorithme solide pour le calcul de {{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 {{Formule|2''N''}}-gone :
:: <math> A_{2N} = m_{N} \times r</math>
: Pour un 12288-gone inscrit dans un cercle de rayon unitaire :
:: <math> A_{24576}=3,14159261864 < \pi </math> .

: De l'inégalité de Liu Hui pour {{MathPi}} :
:: <math> A_{24576} < \pi < A_{24576} +D_{24576}</math>
: Dans lequel <math> D_{24576}=A_{24576}- A_{12288}=0,0000001021</math>
:: <math> A_{24576}=3,14159261864 < \pi <3,14159261864 +0,0000001021</math> .
: Donc
:: <math>3,14159261864 < \pi <3,141592706934</math>

Tronqué à huit chiffres significatifs :

:: <math> 3,1415926 < \pi <3,1415927</math> .

C'est la fameuse inégalité de Zu Chongzhi pour {{MathPi}}.

Zu Chongzhi a ensuite utilisé la formule d'interpolation de [[He Chengtian]] ([[:zh:何承天 (南朝)|何承天]], 370-447) et a obtenu une fraction approximative : <math>\pi \approx {355 \over 113}</math> .

Cependant, cette valeur de {{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é <math>\pi = {22\over 7}</math> et <math>\pi=\sqrt {10})</math>), 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 <math>\pi \approx {355 \over 113}</math> <ref>Yoshio Mikami said about the work of Zhao Yu Xin:"The sides and consequently the perimeters of these polygons are successively calculated in such a manner as followed by Liu Hui of old", p136, Development of Mathematics in China and Japan</ref>.

== 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 {{Formule|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 < {{MathPi}} < 3,1415927, qui a détenu le record mondial de la valeur la plus précise de {{MathPi}} pendant des siècles<ref>Robert Temple, The Genius of China, a refined value of pi, p144-145, {{ISBN|1-85375-292-4}}</ref>, jusqu'à ce que [[Madhava de Sangamagrama]] calcule 11 chiffres au {{S-|XIV}} ou [[Al-Kashi|Jamshid al -Kashi]] calcule 16 chiffres en 1424 ; les meilleures approximations de {{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]] ({{-s-|V}})
* Algorithme de Zhao Youqin pour {{mathPi}} ({{Sp-|XIII|-|XIV|s}})
* {{en}} [https://proofwiki.org/wiki/Newton%27s_Formula_for_Pi Proof of Newton's Formula for Pi] ({{S-|XVII}})

== Remarques ==

: {{Refa|1|1}} La valeur correcte est de {{formatnum:0.2502009052}}
: {{Refa|2|2}} Les valeurs correctes sont :
:: <math>A_{192}= 3,1410319509</math>
:: <math>D_{192}=0,0016817478</math>
:: <math>\pi \approx A_{192}+ \frac{1}{3} D_{192}\approxeq 3,1410319509 +0,0016817478/3</math>
::<math>\pi \approx 3,1410319509 +0,0005605826</math>
::<math>\pi \approx 3,1415925335.</math>
::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.

== Références ==
{{Traduction/Référence|en|Liu Hui's π algorithm|1144603846}}
{{Références|taille=36}}

== BIbliographie ==

*{{ouvrage |lang=en |auteur1=Joseph Needham |année=1986 |titre=Science and Civilization in China |sous-titre=Volume 3, Mathematics and the Sciences of the Heavens and the Earth |lieu=Taipei |éditeur=Caves Books, Ltd}}.
* Wu Wenjun ed, ''History of Chinese Mathematics'' Vol III (in Chinese) {{isbn|7-303-04557-0}}

{{Portail|géométrie}}

[[Catégorie:Mathématiques chinoises]]
[[Catégorie:Pi]]