割圓術 (劉徽)

劉徽割圓術原理
三國時代數學家劉徽的割圓術是中國古代數學中“一個十分精彩的算法”。在此之前，圓周率採用“徑一週三”的實驗數據。東漢科學家張衡採用和。劉徽認為過大。。孫吳天文學家王蕃採用。這些圓周率都是實驗值，都只准確到二位數字。劉徽是中國數學史上最先創造了一個從數學上計算圓周率到任意精確度的迭代程序的數學家。他自己通過分割圓為192邊形，計算出圓周率在3.141024 與 3.142704之間，取其近似，並以 表示。這個數值準確到三位數字，比前人的圓周率數值都準，但他自己次承認這個數值偏小。後來劉徽發明一種快捷算法，可以只用96邊形得到和1536邊形同等的精確度，從而得令他自己滿意的。
劉徽割圓術簡單而又嚴謹，富於程序性，可以繼續分割下去，求得更精確的圓周率。南北朝時期著名數學家祖沖之用劉徽割圓術計算11次，分割圓為12288邊形，得圓周率=3.1415926，成為此後千年世界上最準確的圓周率。

**圓面積公式**
圓面積=圓的半周長X半徑
劉徽割圓術是建立在圓面積論的基礎之上的。他首先論證，將圓分割成多邊形，分割來越細，多邊形的邊數越多，多邊形的面積就和圓面積沒有差別了。他說，將6邊形一邊的長度乘以圓半徑，再乘3，得12邊形的面積。將12邊形的一邊長乘半徑，再乘6，得24邊形面積。越割越細，多邊形和圓面積的差越小。如此割了再割，最後終於和圓合為一體，毫無差別了。
6邊形的面積顯然和圓面積相差很多。
內接正12邊形面積 = 6邊形面積+6個藍色三角形面積，向圓面積趨近了一步。
正24邊形面積=6邊形面積+6個藍色三角形面積+12個黃色三角形面積，更加接近圓面積了。
顯然：
正12邊形面積 N邊形面積N邊形面積   圓面積。
他又指出：6邊形之外，遺留了半徑的一小段d ，稱為餘徑。將餘徑d乘多邊形的一邊，所得長方形ABCD，已經越出圓周範圍之外。如果將圓周分割得很細，餘徑d趨向於0，而長方形ABCD的面積也趨向於0。
顯然，劉徽之所以研究餘徑，目的是從上限和下限兩個方面逐步逼近圓面積：
  內接 6*2N邊形面積   圓面積  內接 6*2N邊形面積+6*2N*d*L。
劉徽進一步證明圓面積=圓周/2  ×  半徑。
關於多邊形的面積，劉徽有如下公式：
2 N邊形的面積= N邊形的半周長×R。
**,
其中L為N邊形的單邊長，R為圓半徑。
此公式可用劉徽出入相補原理證明： 將內接2N邊形，分割，然後重新排列成寬為  L x N/2,  高為R的長方形；
顯然2N邊形的面積=長方形面積==N邊形的半周長 * R
當
N邊形的半周長圓的半周長
 2N邊形面積=N邊形的半周長 *  R 圓面積
所以
圓的半周長 * R = 圓面積
因此
圓周  = 2* 圓面積/R
圓周率圓周/直徑=  2* 圓面積/(R*2R)= 圓面積/R2
** 2N邊形的面積/R2

**割圓術程序**
劉徽從半徑1尺圓的內接正6邊形開始，逐次分割為12邊形，24邊形，48邊形，96邊形。反覆使用勾股定理求得各多邊形的邊長，又用劉氏多邊形面積公式求多邊形面積。

**分割6邊形為12邊形**
劉徽割圓術原理
令圓直徑為2尺，折半得半徑1尺。圓內接正6邊形的邊長也是1尺。
如圖：
半徑OA=r=1尺=10寸
6邊形單邊長AB=M=10寸
從圓心O作AB的垂直平分線OC，將AB平分為二，
:AP=BP=M/2，AP+BP=AB
垂直平分線OC和圓周相交於C，
作直線AC
AC就是12邊形的一邊，
OAP是一個直角三角形
弦=半徑=r=10寸
勾=AP=M/2=5寸
股OP 可用勾股定理求得：
令弦長=X，股長=G， 句長=M/2，則：
 平方寸

因為1寸 =100000忽
1平方寸 =10000000000平方忽
忽
APC是一個小直角三角形
令小弦AC長度為m，令小句PC長度為j
忽
用勾股定理求m：

**平方忽
12邊形的一邊長度忽
:
12邊形的一邊長度的一半忽

**分割12邊形為24邊形**
將上一輪的多邊形邊長m作為新一輪割圓的開始，
作替換M=m=12邊形的一邊長度忽
繼續將此多邊形的一邊平分，週而復始，重複使用：




由上M^2已有現成數值




:24邊形一邊長度 

**分割24邊形為48邊形**
將第二輪的多邊形邊長m作為第三輪割圓的起點，
作替換

:



開平方，得48邊形一面忽
根據劉徽多邊形面積公式：
96邊形的面積= 48邊形的半周長×半徑=,
所以96邊形的面積
平方忽
 平方寸

**分割48邊形為96邊形**
將第三輪的多邊形邊長m作為第四輪割圓的起點
作替換忽

:



開方得
96邊形的一邊忽
根據劉徽多邊形面積公式：
192邊形的面積96邊形的半周長×半徑=
所以192邊形的面積 平方忽
平方忽
 平方寸

**劉徽圓周率不等式**
劉徽圓周率不等式示意圖
劉徽利用多邊形面積差的幾何學，得出圓周率的雙邊不等式。
如圖：
黃色代表N邊形面積
黃色+綠色代表2N邊形面積
綠色代表2N邊形面積與N邊形面積之差=
長方形ABCD面積
C代表圓面積。
如下不等式成立：

或

當N=96，2N=192：
192邊形面積 
96邊形的面積 
:
192邊形面積和96邊形的面積之差（差冪）




即

劉徽認為這個面積已經超過圓面積，所以將192邊形的面積的整數部分定為圓面積：
圓面積~192邊形面積=
所以圓周率=圓面積/半徑2
這就是徽率。
實際上只要計算精確度夠高，劉徽割圓術可以計算到任何精確度，不僅限於二位小數點。

**圓周率捷法**
劉徽在得圓周率=3.14之後，將這個數值和晉武庫中漢王莽時代製造的銅製體積度量衡標準嘉量斛的直徑和容積檢驗，發現3.14這個數值還是偏小。於是繼續割圓到1536邊形，求出3072邊形的面積，得到令自己滿意的圓周率。但是劉徽卻不敘述“分割96邊形為192邊形”，“分割192邊形為384邊形”，“分割384邊形為768邊形”，“分割768邊形為1536邊形”：因為他發現了一個快捷的算法，只要利用96邊形的數據經過一次除法和一次加法，就可以獲得和計算到1536邊形同等的精確度 ，省去了4次開方計算；畢竟在三國時代用籌算進行開方相當的繁難。
劉徽圓周率捷法乃是以他素有研究的多邊形面積差為基礎的。
令表示2N邊形的面積和N邊形的面積差
，，，形成一個等比級數:



因此


其中


劉徽圓周率捷法，可以解釋如下幾個問題：time=2009-03-22T17:57:36+00:00
1)為什麼劉徽割圓術以多邊形面積為基礎，因為圓周率捷法必須用到多邊形面積差。
2)劉徽對割圓術的陳述為什麼止於96邊形。因為他發明了一個便捷的方法，只用96邊形數據，就可以算出相當於1536邊形(甚至12288邊形)的精確度。
3)晉武庫一段的作者，非劉徽莫屬，而不可能出自祖沖之。面積差法本來就是他推求不等式 的基礎。從 到一脈相承。何況九章算術中全無“祖沖之注”的痕跡，而且一字不提祖沖之密率。
劉徽的 = 後來見於印度數學中，足證古印度數學採用劉徽注《九章算術》

**割圓術迭代公式**
如令半徑=1，
從



可簡化為：






由此可得劉徽割圓術迭代公式：
:            
圓周率= 3*2^N * m

**π的連平方根表示式**
根據劉徽割圓術迭代公式：
:            
圓周率= 3*2^N * m
從半徑=1的內接6邊形開始：
各多邊形的一邊長m：





半徑=1圓形正內接多邊形面積：





**祖率**
南北朝數學家祖沖之，並沒有發明新的方法計算圓周率，而是將劉徽割圓術的計算，繼續分割到12288邊形，又用劉徽多邊形面積公式，求得24576邊形的面積：
。
再用劉徽圓周率不等式：
。
其中：

。
得不等式：
。
取八位有效數字即得祖沖之著名的圓周率不等式：
。
祖沖之算得的圓周率準確到小數點後7位，保持了世界最準確圓周率達900年之久。祖沖之熟悉何承天調日法，以3為弱率， 以4為強率，通過調日法計算7次得圓周率約率，計算23次得密率。
根據調日法計算出來的約率和密率都是強率；所謂約率只意味這個數值和圓周率的誤差較大，並無約率“小於”圓周率的意思。

**和阿基米德割圓術比較**
比劉徽早出生500年的希臘數學家阿基米德用阿基米德割圓術計算圓周率，他的論證以計算線長為依據，在推導過程中不考慮多邊形面積面積，和劉徽的以面積計算為中心的割圓術成對照。他用兩套不同的方法方法，先多次分割圓的切線，證明π>；另用內接多邊形，計算到96邊形，證明π十七次調日值　阿基米德弱值 3.140845 ＜ 劉徽弱值 3.141024 ＜ π ＜ 二十三次調日值　祖沖之密率 3.14159292035 ＜ 劉徽強值 3.142704 ＜ 七次調日值　阿基米德強值 3.142857 。
。

**參見**
割圓術 (趙友欽)
*《九章算術》

**參考**

**延伸閱讀**
吳文俊主編 《中國數學史大系》第三卷 第一章第三節 劉徽的割圓術  152-164頁  ISBN 7-303-04557-0
*傅海倫編著 《中外數學史概論》 第四章 第三節 劉徽的割圓術  49-52頁 ISBN 987-7-03-018477-1