圓周率

来自testwiki
imported>Wolfch2025年3月21日 (五) 00:15的版本 撤销207.174.6.43讨论)的版本86512036)
(差异) ←上一版本 | 最后版本 (差异) | 下一版本→ (差异)
跳转到导航 跳转到搜索

Template:NoteTA {{#invoke:TemplateVariadicArgumentSingle|build_template |_core_template=Template:Infobox number/core |_core_args=lang |_core_insert_code= | lang$ = {{{lang$|}}} | lang$ symbol = {{{lang$ symbol|}}} }} Template:Numbers 圓周率数学常数,等於任何周长和其直径,一個常見的近似值等於Template:計算,常用符号π表示。

π无理数,不能用分数表示出来(即它的小数部分是无限不循环小数),但近似227等有理数。學界認為π的数字序列在统计上随机分布,但迄今未能证明。此外,π还是超越数,亦即它不是任何有理系数多项式化圆为方的问题不可能用尺规作图解决。

几个文明古国很早就須计算出π的精确值以便于生产的计算。西元5世纪,中國劉宋数学家祖冲之用几何方法将圆周率计算到小数点后7位。大约同时,印度数学家也将圆周率计算到小数点后5位。史上首條π的精确无穷级数公式(即π的莱布尼茨公式)直到约1000年后才由印度数学家发现。[1][2]微積分出現,π的位數很快計到數百位,足以滿足任何科學工程的計算需求。在20和21世纪,计算机技术快速发展,π的计算精度急速提高。截至2024年3月,π的十进制精度已達105万亿位。[3]几乎所有科学研究对π的精度要求都不超过几百位,当前计算π的值主要都为打破记录、测试超级计算机的计算能力和高精度乘法算法[4]Template:Rp[5]

π的定义涉及圆,在三角学几何学的许多公式,特别是广泛应用在圆形、球形或椭球形相關公式中。[6]在近代數學分析裡,π改由實數系統譜性質中的特征值週期來定義,其他數學領域數論統計以及幾乎所有物理學領域均有出現,π的广泛用途使它成为科学界内外最广为人知的数学常数。几本专门介绍π的书籍经已出版,圆周率日(3月14日)和π值计算突破记录也往往会成为报纸的新闻头条。[7]此外,背诵π值的世界记录已达10萬位。[8]

Template:TOC limit

直徑為一的圓的周長是π(3.14159265...)

基本概念

名称

数学家用小写希腊字母π表示圆周和其直径之比,有时也将其拼写为“Template:Lang”,来自希腊语“Template:Lang”(周长)的首字母。[9]英语π的发音与英文单词“Pie”(Template:IPAc-en西式馅饼)相同。[10]π的小写字母(或其无衬线体)在数学要和表示连乘积的大写Π相区分开。

关于选择符号π的原因,请参见引入π符号一节。

定义

A diagram of a circle, with the width labeled as diameter, and the perimeter labeled as circumference
圆周长略大于其直径的三倍;精确的比例称为π

π常用定义为周长C直径d的比值:[4]Template:Rp

π=Cd

无论圆的大小如何,比值Cd为恒值。如果圆的直径变为原先的二倍,周长也变为二倍,比值Cd不变。π目前的定义暗地用了欧几里得几何的一些定理,虽然圆的定义可扩展到任意曲面(即非欧几里得几何),但这些圆不符合定律π=Cd[4]

这里,圆的周长指其圆周的弧长,弧长这概念可以不依赖几何学,而是用微积分学的极限来定义。[11]例如,若想计算笛卡儿坐标系中单位圆x2+y2=1上半部分的弧长,需要用到积分[12]

π=11dx1x2.

上述积分是由卡尔·魏尔斯特拉斯于1841年对π的积分定义。[13]

π这些依赖周长、且暗地依赖积分的定义如今在文献中并不常见。雷默特(Template:Harvtxt)解释说现代教微积分時,大学一般将微分学课程安排在积分学课程之前,所以不依赖于后者的π的定义就很有必要了。其中一种定义由Template:Le提出,[14]愛德蒙·蘭道推广,[15]其表述如下:π是两倍于能使余弦函数等于零的最小正数。[4][12][16]余弦函数可以由独立于几何之外的幂级数[17]定义,或者使用微分方程的解来定义。[16]

在相似的启发下,π可以用关于复变量z复指数函数exp(z)来定义。复指数类似余弦函数,可用多种方式定义。令函数exp(z)值为一的复数集合是如下所示的(虚)等差數列:

{,2πi,0,2πi,4πi,}={2πki|k}

并且其中包括独特的正实数π[12][18]

基于同样想法但更抽象的定义运用了精巧的拓扑学代数学概念,用以下定理描述:[19]存在唯一的从加法模数整数组成的实数群R/Z到绝对值为1的复数组成的乘法群的连续同态(拓扑学概念,指在拓扑空间之间的一种态射)。数字π定义为此同态派生的模的一半。[20]

周长固定,圆会围成最大面积,π同樣表述为等周不等式中出现的常数(乘四分之一)。此外,在很多其他紧密相关的方程中,π作为某些几何或者物理过程的特征值出现;详见下文

无理及正规性

π是无理数,无法表示成两整数之比的形式(形如227的分数常用来近似表达π,但是没有任何普通分数(指整数的比)可以取到π的精确值)。[4]Template:Rp由于π是无理数,故可表示为无限不循环小数。有多种方法能证明π是无理数,这些证明也都要用到微积分学反证法π可以用有理数来近似的程度還無法準確得知(稱為無理性度量),不過估計其無理性度量比Template:MathTemplate:Math的要大,但是小於刘维尔数的無理性度量[21]

Template:Tsl检验,包括正规数检验,可验证π的位數沒有明顯的固定模式。π的小数中任意固定长度的序列(如3位數000,001……999)出現機率都相同[22]。不過有關Template:Pi正规数的猜想既無證明,亦無-{zh-hans:证伪;zh-hant:否證;}-[4]Template:Rp[22]

電腦出現後可生成大量Template:Pi的不同位数,并統計分析之。金田康正詳細統計分析了Template:Pi的十進制數字,并验证了其分布正规:例如,假設檢定0到9十個數的出現頻率,找不到有特定重复规律的證據[4]Template:Rp。根據無限猴子定理,任何任意長度、由隨機內容組成的子序列看起來都有可能像不隨機生成。因此,就算Template:Pi的小数序列通過了隨機性統計測試,其中也可能有幾位的數字看起來似有规律可循而非隨機数,例如Template:Pi的十進制写法在小數第762位后开始出现了連續六個9[4]Template:Rp

超越性

A diagram of a square and circle, both with identical area; the length of the side of the square is the square root of pi
由于Template:Pi超越數,不能利用尺规作图化圓為方

π不仅是无理数,还是超越数,即π不是任何有理系数多项式。(比方说,试图解有限项方程x5120x36+x=0来求π的值)[23]Template:NoteTag

π的超越性衍生出一些重要的结果:π不能经有限次四则运算和开平方运算有理数来获得,因此不是规矩数。换言之,尺规作图作不出长度为π的线段,也就不可能用尺规方法做出与已知圆面积相等的正方形。后者即为有名的化圓為方问题,该问题早在古典时代即已提出,曾困扰人数千年之久[24][25]。直至今天,依然有民间数学爱好者声称他们解决了这问题[26]

连分式

π像所有无理数一样无法表示成分数,但π等全部无理数都能表示成一系列叫连分数的连续分数形式:

π=3+17+115+11+1292+11+11+11+

在这连分数的任意一点截断化简,都能得到π的近似值;前四位近似值是3、227333106355113。这些数在历史上是π最广为人知且广為使用的几个近似值。用以上方式得出的π的近似值要比任何有相同或更小的整数分母的其他整数分数近似值更接近π。[27]π是超越数,据定义来说它不是代數數,又因此不可能是二次無理數;是故π不能表示为循环连分数。尽管π的简单连分数没有表现出任何其他明显规律,[28]数学家發現了数條广义连分数能表示π,例如:[29]

π=41+122+322+522+722+922+=3+126+326+526+726+926+=41+123+225+327+429+

近似值

圆周率近似值包括:

其他进位制的近似值

复数与欧拉恒等式

在复平面上以原点为圆心的单位圆内,一条射线从圆心出发至圆的边上,以此射线与圆的边的交点作与x轴的垂线并标注了夹角Template:Math和Template:Math、Template:Math函数
欧拉公式给出了Template:Math的复指数与复平面上以原點为圆心的单位圆上的之间的关系。

任何复数(以z为例)都可以表示为一组实数对:极坐标系用实数r表示半径,代表复平面上复数z离原點的距离;实数φ则表示夹角,即这条半径(复平面上复数z与原点的连线)与正实轴经顺时针转动的夹角。这样一来,z就可写成[32]

z=r(cosφ+isinφ),这里i代表虛數單位,即i2=-1。

复分析中,欧拉公式三角函数与复指数函数糅合在一起[33]

eiφ=cosφ+isinφ,这里[[E (数学常数)|数学常数Template:Math]]是自然對數的底数。

欧拉公式确立了e的复指数与复平面上以原点为圆心的单位圆上的点之间的关系,而且当φ=π时,欧拉公式就能改写为歐拉恆等式的形式:

Template:計算結果。此等式亦稱“最奇妙的数学公式”(Template:Lang-en),全因它将五个最基本的数学常数简洁联系起来[33][34]

欧拉等式亦可用于求出方程zn=1n个不同复数根(这些根叫做n单位根[35]),可以根据以下公式求得:

e2πik/n(k=0,1,2,,n1)

谱特征

震盪弦的泛音是二次微分的本徵函數,會形成泛音列。對應的本徵值會形成由Template:Pi整數倍組成的等差数列

π出现在有关几何的问题中。然而,不少和几何无关的问题也可看到π的身影。

π在許多用處中都會以特征值形式出現。例如理想的Template:Tsl問題可以建模為函數f在單位區間[0,1]的圖形,固定邊界值f(0)=f(1)=0。弦振動的模態會是微分方程fn(x)+λ2f(x)=0,此處Template:Math是相關的特徵值。受施图姆-刘维尔理论限制,λ只能是一些特定的數值。而λ=π即為一個特征值,因為函數f(x)=sin(πx)滿足邊界條件及微分方程λ=π[36]

依照第一代开尔文男爵威廉·汤姆森所述的一篇傳說,古迦太基城的外形是等周長問題的一項解Template:Harv。這些包圍著海的區域由迦太基女王狄多所圍,城不靠海的邊界須用指定大小的牛皮圍住,後來是將牛皮剪成小段

π是上述方程的最小特征值,也和弦振動的Template:Tsl有關。一種讓弦振動的方式是提供弦能量,能量會滿足維廷格函數不等式[37],其中提到若函數f:[0,1]使得f(0)=f(1)=0,且ff都是平方可積函數,則以下的不等式成立:

π201|f(x)|2dx01|f(x)|2dx,

此例中等號成立的條件恰好是fsin(πx)倍數的時候。因此π似乎是維爾丁格不等式的最佳常數,也是最小的特征值(根據雷利商數的計算方式)

π在更高維度的分析也有類似的角色,出現在其他類似問題的特徵值中。就如以上所述π的一項特點是等周定理中的最佳常數:周長為P的平面若尔当曲线,所圍面積A滿足以下的不等式

4πAP2

A=πr2P=2πr,故等號成立的條件是曲线為圓形[38]

圓周率Template:Pi也和庞加莱不等式的最佳常數有關[39]π是一維及二維的Template:Le特征向量最佳值中最小,會出現在許多經典的物理現象中,例如經典的位势论[40][41][42]。其一維的情形即為維廷格不等式。

圓周率Template:Pi也是傅里叶变换的重要常數,傅里叶变换屬於积分变换,將實數線上有複數值、可積分的函數,轉換為以下形式:

f^(ξ)=f(x)e2πixξdx.

傅里叶变换有幾種不同的寫法,但不論怎麼寫,傅里叶变换及反傅里叶变换中,一定會有某處出現π。不過上述的定義是最經典的,因為其描述了[[Lp空间|Template:Math空間]]中唯一的幺正算符,也是L1空間到L空間的代數同態[43]

不确定性原理也用到π。不确定性原理提出了可以將函數在空間及在頻域中局部化程度的下限,用傅立葉轉換的方式表示:

x2|f(x)|2dx ξ2|f^(ξ)|2dξ(14π|f(x)|2dx)2

物理的結果,有關量子力学中同時觀測位置及動量的不確定性,見下文傅立葉分析中出現Template:PiTemplate:Tsl的結果,證實了海森伯群Template:Tsl是唯一[44]

高斯积分

高斯函数f(x)=ex2的图像,函数下方与X轴围成的阴影部分面积为π

高斯积分是对高斯函数ex2在整条实轴上的积分,即函数下方与X轴围成的面积,其结果为π

ex2dx=π

此积分的计算可以先计算f(x)=ex2对整条实轴的积分的平方,通过转换笛卡尔坐标系极坐标系从而求得

(ex2dx)2=𝐑2e(x2+y2)dxdy=02π0er2rdrdθ=π

其他计算方法可参阅高斯积分。高斯函数更一般的形式为f(x)=aexp(xb)22c2,求一般形式的高斯积分均可通过换元积分法转化为求f(x)=ex2的积分。

另外,当高斯函数为以下形式时,它则是平均数μ標準差σ正态分布機率密度函數[45]

f(x)=1σ2πexp(xμ)22σ2

这函数是概率密度函数,函数下方与X轴围成的面积必须为1,令μ=0σ=1即可变换得出ex2dx=π概率论统计学领域经常使用正态分布来作为复杂现象的简单模型:例如科学家通常假设大多数试验观测值的随机误差都是服从正态分布[46]

由一维布朗运动的反正弦定律,可以通过试验正信号相对于负信号领先权过零点的分布反过来推算Template:Pi

概率论与统计学中的中心极限定理解释了正态分布以及π的核心作用,这定理本质上是联系着π谱特征海森堡不确定性原理相关的特征值,并且在不确定性原理中有

σxσp/2

这里的σxσp分別為位置與動量的標準差約化普朗克常数,而不等式的等号当且仅当粒子的波函数为高斯函数使成立[47]

同样地,π作为唯一独特的常数使得高斯函数等于其自身的傅里叶变换,此时的高斯函数形式为f(x)=eπx2[48]。根据豪(Template:Lang)的说法,建立傅里叶分析基本定理的“全部工作(whole business)”简化为高斯积分。

历史

Template:Main article Template:See also

远古时期

圓周率在远古时期(西元前一千纪)已估算至前两位(3.1)。有些埃及學家聲稱,遠至古王國時期時期的古埃及人已經用227作為圓周率的約數[49]Template:NoteTag,但這說法受到質疑。[50][51][52][53]

最早有記載的对圓周率估值在古埃及巴比伦出现,兩估值都与圆周率的正确数值相差不到百分之一。巴比伦曾出土一塊西元前1900至1600年的泥板,泥板上的幾何學陳述暗示人们当时把圓周率視同258(等於3.125)。[4]Template:Rp埃及的莱因德数学纸草书(鉴定撰寫年份為西元前1650年,但抄自一份西元前1850年的文本)載有用作計算圓面積的公式,该公式中圓周率等于(169)2(≈3.1605)。[4]Template:Rp

西元前4世紀的《Template:Link-en》的天文學運算把339108(≈3.139,精确到99.91%)用作圓周率估值[54]。西元前150年前其他印度文獻把圓周率視為10(≈3.1622)[4]Template:Rp

割圆时代

圖中有圓的外切五邊形、內接五邊形、外切六邊形及內接六邊形
Template:Pi可以透過計算圓的外切多邊形及內接多邊形周長來估算

第一條有紀錄、嚴謹計算Template:Pi數值的演算法是用正多邊形的幾何算法,在西元前250年由希臘數學家阿基米德發明。[4]Template:Rp這算法用了有一千年之久,因而有時Template:Pi亦稱阿基米德常數。[4]Template:Rp阿基米德的算法是在計算圓的外切正六邊形及內接正六邊形的邊長,以此計算π的上限及下限,之後再將六邊形變成十二邊形,繼續計算邊長,一直計到正96邊形為止。他根據多邊形的邊長證明22371<π<227(也就是3.1408<π<3.1429[55]。阿基米德得到的上限227也造成常見誤解,認為π就等於227[4]Template:Rp。在西元前150年,希臘羅馬的科學家克劳狄乌斯·托勒密在《天文学大成》一書中提到Template:Pi的數值是3.1416,可能來自阿基米德,也可能來自阿波罗尼奥斯[4]Template:Rp[56]數學家在1630年利用多邊形的方式計算Template:Pi到第39位小數,一直到1699年,其他數學家才利用無窮級數的方式打破其紀錄,計算到第71位小數[57]

獨自研究圖形的阿基米德
阿基米德發展了用多邊形近似Template:Pi的計算方式

中国历史上,π的數值有3[58]、3.1547(西元前一世紀)、10(西元前100年,數值約3.1623)及14245(第三世紀,數值約3.1556)[4]Template:Rp。大約在西元265年,曹魏數學家刘徽創立割圆术,用正3072邊形計算出Template:Pi的數值為3.1416。[59][4]Template:Rp他後來又發明了較快的算法,利用邊數差兩倍的正多邊形,其面積的差值會形成等比數列,其公比為14的原理,配合96邊形算出π的值為3.14。[59]祖冲之在西元480年利用割圆术計算12288邊形邊長,得到π的值在3.1415926和3.1415927之间。他同时提出了π的约率227和密率355113。在之後的八百年內,這都是Template:Pi最準確的估計值。[4]Template:Rp為紀念祖沖之對圓周率發展的貢獻,日本數學家三上義夫將這推算值命名為“祖沖之圓周率”,簡稱“祖率”。[60]

印度天文學家阿耶波多在西元499年的著作《阿里亞哈塔曆書》中使用了3.1416的數值。[4]Template:Rp斐波那契在大約1220年用獨立於阿基米德多邊形法,計算出3.1418[4]Template:Rp。義大利作家但丁·阿利吉耶里用的數值則是3+2103.14142[4]Template:Rp

波斯天文學家卡西在1424年利用3×228邊的多邊形,計算到六十進制的第9位小數,相當十進制的第16位小數。[61][62]這一突破成為當時的紀錄,延續了約180年。[63]法國數學家弗朗索瓦·韦达在1579年用3×217邊形計算到第9位小數[63],佛蘭芒數學家阿德里安·范·羅門在1593年計算到第15位小數[63]。荷蘭數學家鲁道夫·范·科伊伦在1596年計算到第20位小數,他之後又計算到第35位小數(因此在二十世紀初之前,圓周率在德國會稱為鲁道夫數)。[4]Template:Rp荷蘭科學家威理博·司乃耳在1621年計算到第34位小數[4]Template:Rp,而奧地利天文學家Template:Tsl在1630年用1040邊形計算到第38位小數[64],至今這仍是利用多邊形算法可以達到最準確的結果[4]Template:Rp

无穷级数

Template:Comparison pi infinite series.svg 16及17世紀時,開始改用無窮级数的方式去計π。無窮级数是一組無窮數列的和[4]Template:Rp。無窮级数讓數學家可以計算出比阿基米德以及其他用幾何方式計算的數學家更準確的結果。[4]Template:Rp雖然詹姆斯·格雷果里戈特弗里德·莱布尼茨等歐洲數學家利用無窮數列計算Template:Pi而使得该方法为大家所知,但这种方法最早是由印度科學家在大約1400到1500年之間發現。[4]Template:Rp[65]第一個记载用無窮级数計算Template:Pi的人是约西元1500年左右时,印度天文學家Template:Link-en在他的著作《Template:Tsl》中用梵語詩所記錄。[66]當時沒有這數列對應的證明,而證明出現在另一本較晚的印度作品《Template:Tsl》,年代約在西元1530年。尼拉卡莎將該數列歸功於更早期的印度數學家Template:Link-en(1350–1425)。[66]相關的無窮级数有許多,包括有關sintancos的,現在稱為Template:Link-enπ的莱布尼茨公式[66]。瑪達瓦在1400年用無窮级数計算Template:Pi到第11位小數,但在1430年一位波斯數學家卡西利用多邊形算法否定了他算的結果[67]

長髮艾萨克·牛顿的畫像
艾萨克·牛顿利用無窮级数計算Template:Pi到第15位,後來寫道:「我很羞愧的告訴你我為了這個計算用了多少個數字。」[68]

歐洲發現的第一條無窮項圓周率公式無窮乘積(和一般用來計算Template:Pi的無窮級數不同),由法國科學家弗朗索瓦·韦达在1593年發現[4]Template:Rp[69]

2π=222+222+2+22

約翰·沃利斯在1655年發現了沃利斯乘积,是歐洲發現的第二條無窮項圓周率公式[4]Template:Rp

π2=2123434565678789

微积分学由英國科學家艾萨克·牛顿及德國數學家戈特弗里德·莱布尼茨在1660年代發明,許多計Template:Pi的無窮級數出現。牛頓自己就用反正弦arcsin)數列在1655年或1666年將Template:Pi近似到第15位小數,後來寫到「我很羞愧告訴你我為了計算它用了多少數字,我當時沒有做其他事。」[68]

蘇格蘭數學家詹姆斯·格雷果里在1671年發現了馬德哈瓦公式,莱布尼茨也在1674年發現:[4]Template:Rp[70]

arctanz=zz33+z55z77+

這公式即為格雷果里-莱布尼茨公式,在z=1時數值為π4[70]1699年時英國數學家亚伯拉罕·夏普用格雷果里-莱布尼茨公式,在z=13時計算,計算到π的第71位小數,打破由多邊形算法得到的第39位小數的记录。[4]Template:Rp格雷果里-莱布尼茨公式在z=1時非常簡單,但收斂到最終值的速度非常慢,現在不会再用此公式來計π。[4]Template:Rp

約翰·梅欽在1706年用格雷果里-莱布尼茨級數產生了可以快速收斂的公式:[4]Template:Rp

π4=4arctan15arctan1239

梅欽用這公式計到π第100位小數[4]Template:Rp後來其他數學家也發展了一些類似公式,現在稱為梅欽類公式,創下了許多計算π位數的紀錄。[4]Template:Rp在進入電腦時代時,梅欽類公式仍然是耳熟能详可以計算π的公式,而且在约250年的时间里,很多有關π位數的紀錄都是梅欽類公式所得,比如在1946年時由達尼爾·弗格森(Template:Lang)用這類公式計到第620位小數,是沒有計算設備輔助的最佳紀錄。[4]Template:Rp

1844年,計算天才Template:Tsl在德國數學家卡爾·弗里德里希·高斯的要求下以梅欽類公式心算了π的200位小數,並創下紀錄。[4]Template:Rp英國數學家Template:Tsl花了15年的時間計算Template:Pi到小數707位,不過第528位小數出錯,後面的小數也都不正確。[4]Template:Rp

收敛速度

有些Template:Pi的無窮級數收斂的比其他級數要快,數學家一般會選用收斂速度較快的級數,可以在較少的計算量下計算π,且達到需要的準確度[71][4]Template:Rp。以下是Template:Pi莱布尼茨公式[4]Template:Rp

π=4143+4547+49411+413

隨著一項一項的值加入總和中,只要項次夠多,總和最後會慢慢接近π。不過此數列的收斂速度很慢,要到50萬項之後,才會精確到π的第五位小數[72]

尼拉卡莎在15世紀發展了π的另一條無窮級數,收斂速度比格雷果里-萊布尼茨公式快很多:[73]

π=3+42×3×444×5×6+46×7×848×9×10+

以下比較兩條級數的收斂速率:

π的無窮級數 第1項 前2項 前3項 前4項 前5項 收斂到
π=4143+4547+49411+413. 4.0000 2.6666… 3.4666… 2.8952… 3.3396… 3.1415…
π=3+42×3×444×5×6+46×7×8. 3.0000 3.1666… 3.1333… 3.1452… 3.1396…

計算前五項後,格雷果里-萊布尼茨級數的和跟π的誤差為0.2,而尼拉卡莎級數和的誤差為0.002。尼拉卡莎級數收斂快很多,也甚為適合用來計π的值。收斂更快的級數有梅欽類公式楚德诺夫斯基算法,後者每計一項就可以得到14位正確的小數位[71]

无理与超越性

Template:See also 并非所有和π有关的研究都旨在提高计算它的准确度。1735年,欧拉解决了巴塞尔问题,建立了所有平方数倒数和与π的关系。之后欧拉发现了欧拉乘积公式,得到了π、素数的重要關聯,對日後黎曼ζ函數的研究影響深遠。[74]

π26=112+122+132+142+

1761年,瑞士数学家约翰·海因里希·朗伯正切函数的无穷连分数表达式证明了π是無理數[4]Template:Rp[75]1794年,法国数学家阿德里安-马里·勒让德证明了π2也是无理数。1882年,德国数学家费迪南德·冯·林德曼证明了对任何非零代数数αeα都是超越数,该结论后来由魏尔斯特拉斯推广为林德曼-魏尔斯特拉斯定理。据此定理和欧拉公式,π只能是超越數,進而证实了勒让德和欧拉提出的π超越性猜想。[4]Template:Rp[76]哈代在其著作《数论导引》中则称此证明在提出後,經過希尔伯特施瓦兹和其他一些人化简过。[77]

引入Template:Pi符号

萊昂哈德·歐拉在他1736年到1748年的作品中開始用希臘字母Template:Pi表示圓周率,數學界也開始廣為使用

在用Template:Pi专指“圆周率”之前,希腊字母即已用於幾何概念中[4]Template:Rp威廉·奥特雷德在1647年起在《數學之鑰》(Clavis Mathematicae)就已經用πδ(對應p和d的希臘字母)來表示圓的周長及直徑的比例。

威廉·琼斯在他1706年出版的《新數學導論》(Template:Lang)提到了π,是目前已知最早专门用希臘字母π表示圓周和其直徑比例的人[78]。這希臘字母第一次出现是在书中討論一塊半徑1的圓時提到「其圓周長一半(π)」。琼斯選用π可能因它是希臘文“周边”一词“Template:Lang”的首字母[79]。不過琼斯提到,他那些有關π的算式出自「真正聰明的約翰·梅欽先生」,人们推測在瓊斯之前,約翰·梅欽就已开始用π表示圓周率[4]Template:Rp

瓊斯在1706年開始使用此希臘字母,但直到萊昂哈德·歐拉在其1736年出版的《Template:Tsl》中開始使用之后,其他数学家才纷纷开始用π指代圆周率。在此之前,數字家可能用像cp之類的字母代表圓周率[4]Template:Rp。歐拉與歐洲其他數學家间时常互相写信来往,π的用法迅速傳播开来[4]Template:Rp。1748年歐拉在他的《无穷小分析引论》再一次提到了π,写道:「簡潔起見,我們將此數字寫為ππ等於半徑為1的圓周長的一半。」这表示方式之後也推展到整片西方世界[4]Template:Rp

现代数值近似

计算机时代与迭代算法

一位穿著西裝男士的照片
约翰·冯·诺伊曼所屬的團隊是用數位計算機ENIAC來計Template:Pi的第一隊

Template:Quote box

二十世紀中期计算机技术发展、革新再次引发了計算Template:Pi位數的熱潮。美國數學家约翰·伦奇及李維·史密斯在1949年用桌上型計算機計算到1120位[4]Template:Rp。同年,喬治·韋斯納(George Reitwiesner)及约翰·冯·诺伊曼帶領的團隊利用反三角函数(arctan)的無窮級數,用ENIAC計算到了小數後2037位,花了70小時的電腦工作時間[80]。這紀錄後來多次由其他透過arctan級數计算出的結果打破(1957年到7480位小數,1958年到第一萬位數,1961年到第十萬位小數),直到1973年,小数点后第一百萬位小數經已算出[4]Template:Rp

1980年代有两项發明加速計算了π。第一项是發现了新的迭代法去计Template:Pi的值,計算速度比無窮級數快很多;另一项是發现了可以快速計算大數字乘積的乘法演算法[4]Template:Rp。電腦大部分的工作時間都是在計乘法,這類演算法對現代計Template:Pi格外重要[4]Template:Rp。這類演算法包括嘉良對馬(Karatsuba)算法譚曲(Toom-Cook)乘法及以傅里叶变换為基礎的乘法演算法(傅里叶乘法)[4]Template:Rp

迭代演算法最早是在1975年至1976年间分别由美國物理學家Template:Le及奧地利科學家Template:Le独立提出[4]Template:Rp。這两條演算法没有依赖無窮級數來計算。迭代會重覆特定計算,将前一次的計算結果作为這一次的輸入值,使得計算結果漸漸的趨近理想值。此方式的原始版本其實是在160年前由卡爾·弗里德里希·高斯提出,現在稱為算术-几何平均数算法(AGM法)或高斯-勒让德算法[4]Template:Rp。薩拉明及布蘭特都曾修改之,这算法也稱為薩拉明-布蘭特演算法。

迭代演算法收斂速度比無窮級數快很多,在1980年代以後廣為使用。無窮級數隨著項次的增加,一般來說正確的位數也會增加幾位,但迭代演算法每計算多一次,正確位數會呈几何级数增长。例如薩拉明-布蘭特演算法每計算多一次,正確位數會是之前的二倍。1984年加拿大人Template:TslTemplate:Tsl提出迭代演算法,每計算多一次,正確位數會是之前的四倍,1987年時有另一條迭代演算法,每計算多一次,正確位數會是之前的五倍[81]。日本數學家金田康正使用的演算法在1955年及2002年間創下了若干項紀錄[82]。不過迭代演算法的快速收斂也有其代價,需要的記憶體明顯比無窮級數多[82]

计算π的意义

當數學家發現新的算法、電腦變得普及时,Template:Pi的已知小數位急剧增加。注意垂直坐标使用了对数坐标

一般而言,π值并不需要过于精确便能够满足大部分数学运算的需求。按照約·安(Jörg Arndt)及古里斯佗夫·希奴(Template:Lang)的计算,39位精確度已可将可觀測宇宙圆周的精确度準確至一粒原子大小,足以運算絕大多數宇宙学的计算需求[83]。尽管如此,和π有關的成就往往成為世界各地的新聞頭條;部分人出于對破紀錄的冲动,依然奋力算出π小数点后上千甚至上百萬位[4]Template:Rp[84][85]。此外也有測試超级计算机、測試数值分析算法(包括Template:Tsl)等實際好處。純粹數學這领域也能计算π的位数评定其隨機度[4]Template:Rp

快速收敛级数

一位男士的肖像
斯里尼瓦瑟·拉马努金的肖像,他在印度独立工作时提出了许多计算π的新颖数列。

现代计算π的程序不仅局限于迭代算法。20世纪80与90年代,出现了可用来计算π的新无穷级数,其收敛速度可与迭代算法媲美,而又有着复杂度、内存密集度更低的优势。[82]印度数学家斯里尼瓦瑟·拉马努金是这方面的先驱,他在1914年发表了许多与π相关的公式,这些公式十分新颖,极为优雅而又颇具数学深度,收敛速度也非常快。[4]Template:Rp下式即为一例,其中用到了模方程

1π=229801k=0(4k)!(1103+26390k)k!4(3964k).

这无穷级数收敛速度远快于绝大多数反正切数列,包括梅钦公式。[4]Template:Rp第一位使用拉马努金公式计算π并取得进展的是Template:Tsl,他在1985年算得了小数点后一千七百万位。[4]Template:Rp拉马努金公式开创了现代数值近似算法的先河,此后波尔文兄弟和Template:Tsl进一步发展了这类算法。[4]Template:Rp后者于1987年提出了楚德诺夫斯基公式,如下所示:

1π=126403203/2k=0(6k)!(13591409+545140134k)(3k)!(k!)3(640320)3k.

此公式每计算一项就能得到π的约14位数值[86],因而用於突破圆周率的数位的计算。利用这公式,楚德诺夫斯基兄弟于1989年算得π小数点后10亿(109)位,法布里斯·贝拉于2009年算得2.7千亿(2.7×1012)位,亚历山大·易和近藤滋在2011年算得一万亿(1013)位。[4]Template:Rp[87][88]类似的公式还有Template:Tsl

2006年,加拿大数学家西蒙·普勞夫利用PSLQTemplate:Tsl[89]按照以下模版生成了几條计算π的新公式:

πk=n=11nk(aqn1+bq2n1+cq4n1)

q为[[E的π次方|Template:Mathπ]],k奇数a,b,c是普勞夫计算出的有理常数。[90]

統計模擬法

Template:Multiple image 統計模擬法是以概率统计理论为指导的一类非常重要的计数方法,經大量重复试验计算事件发生频率,按照大数定律(即当试验次数充分大时,频率充分接近概率)可以求得π的近似值[91]布芬(Buffon)投針問題就是其中一項實例:长度l的针随机往画满间距t(lt)的平行线的平面上抛掷n次, 如果针与平行直线相交m次,n充分大就可根据以下公式算出π的近似值[92]

π2nmt

用統計模擬法计π的另一例子是随机往内切四分之一圆的正方形内抛掷大量点,落在四分之一圆内的点的数量与抛掷点的总量的比值会近似于π4[4]Template:Rp[93]

此外还可用随机游走试验,并用統計模擬法计算π值,如抛掷一枚均匀的硬币N次,并记录正面朝上的次数,所得结果中,正面朝上的次数nN服从二項分佈

Pr(nN=m)=(Nm)(12)m(12)Nm

因为硬币均匀,所以Template:Math次试验中每次试验结果相互独立。由此可定义一系列独立的随机变量Xk(k=1,2,),当抛掷结果为正面时Xk=1否则为-1,且Xk=±1且取何值有相同概率(即,正面朝上和背面朝上的概率相同)。对随机变量Xk(k=1,2,,N)求和可得

WN=k=1NXk

Template:Math为“硬币正面朝上的次数”减去“硬币反面朝上的次数”,即可得到m(Nm)=k。变换式子,得m=N+k2,因此

Pr(WN=k)=(NN+k2)12N,其中k=N,N+2,N+4,,N2,N

可证明[94]

E(WN)=0E(WN2)=N,以及E(|WN|)=(NN/2N/22N1)={(N1)!!(N2)!!,若 N偶,N!!(N1)!!,若 N奇。

并且当N变大时,E(|WN|)的值会渐近于2Nπ,因此当N充分大时可根据以下公式算出π的近似值:[95]

π2N|WN|2

和其他计算π值的方法相比,蒙特卡洛方法收敛速度很慢,而且无论实验多少次,都无从得知π的估值已经精确到第几位。因此,当追求速度或精度时,蒙特卡洛方法不适合用来估计π[4]Template:Rp[96]

阀门算法

1995年引入的兩條算法开辟了研究π的新途径。因为每计算出一位数字,該數就会像流过阀门的水一样不会再出现在后续的计算过程中,这种新進算法叫Template:Le[4]Template:Rp[97]这就与无穷级数及迭代算法形成对比——无穷级数和迭代算法自始至终的每一步计算都会涉及到之前所有步骤计算出的中间值。[4]Template:Rp

1995年,美國數學家Template:Link-en和斯坦利·拉比諾維茨(Template:Lang)发明了一种簡單的阀门算法[97][4]Template:Rp[98],其運算速度類似arctan演算法,但速度比迭代算法慢[4]Template:Rp

贝利-波尔温-普劳夫公式(BBP)是另一條阀门算法,屬於一种Template:Tsl。1995年,西蒙·普勞夫等人發現[4]Template:Rp[99]

π=k=0116k(48k+128k+418k+518k+6)

這公式和其他公式不同,可以計算π的任何十六进小數位,而不用計算前面全部小數位[4]Template:Rp。十六进数位可计算得到特定二进数位;想要得到八进制数位的话,计算一、两位十六进小數即可。目前也已發現一些這種演算法的變體,不過還沒有发现針對十進制、可以快速生成特定小數位的位數萃取演算法[100]。位數萃取演算法的一項重要用途是用來確認聲稱是計算到π小數位數的新紀錄:若有聲稱是新紀錄的計算結果出現,先將十進制的數值轉換到十六進制,再用贝利-波尔温-普劳夫公式去確認最後一些位數(用亂數決定),若這些位數都對,就能有一定把握认为此計算結果是对的[88]

1998年到2000年間,分布式计算計畫Template:Tsl貝拉公式(贝利-波尔温-普劳夫公式的一種變體)計算π第1015位,結果是0[4]Template:Rp[101]。2010年9月,有雅虎員工用公司的Apache Hadoop應用程式在上千台電腦計算π在2×1015位开始往后256位,其第2×1015位剛好也是0[102]


利用伽瑪函數計算

伽瑪函數,Γ(n)=(n1)!,可以被用作計算圓周率。

π=Γ(12)2

證明如下:

利用歐拉反射公式

Γ(z)*Γ(1z)=πsin(πz)

z=12

Γ(12)*Γ(112)=πsin(π2)

因為 sin(π2)

Γ(12)2=π






用途

π与圆密切相关,出现在许多几何学和三角学的公式中(特别是与圆、椭圆和球体相关的那些)。 此外,π也出现在其他学科的重要公式中,比如统计学、物理学,傅立叶分析和数论的公式。

几何学与三角学

圆右上四分之一處覆蓋在正方形下的图。
圆的面积等于π乘以阴影部分面积。

π出现在基于圆的几何图形(如椭圆圆锥环面)的面积、体积公式中。下面是一些用到π的常见公式:[6]

  • 半径r的圆周长2πr
  • 半径r圆面积πr2
  • 半径r的球体积43πr3
  • 半径r的球面面积4πr2

上述公式是n维球的体积与其边界(n−1)维球的球面)的表面积的特殊情况,具体将在后文给出解释。

描述由圆生成的图形的周长、面积或体积的定积分常涉及π。例如,表示半径为1的半圆的面积的积分为[103]

111x2dx=π2.
1x2

的积分表示上半圆(此处的平方根勾股定理得出),从-1到1的积分

1x

可用来计算计算半圆与[[笛卡儿坐标系|Template:Math 轴]]间的面积。

函数图象
正弦余弦函数的重复周期为 2Template:Pi

三角函数要用到角,而数学家常用弧度作角度单位。π在弧度制起重要作用,数学家将周角,即360度定义为2π度。[104]由这条定义可得,180度=π弧度,1度=π180弧度。[104]因此,常用的三角函数的周期为π的倍数;例如,正弦和余弦周期为π,[105]任何角度θ和任何整数k都有

sinθ=sin(θ+2πk)cosθ=cos(θ+2πk)[105]

拓扑学

Template:Tsl单值化亏格为3且欧拉特征值为−4的面,作为双曲面Template:Tsl对称群PSL(2,7)的商。根据高斯-博内定理,基本域的双曲面积为8Template:Pi.

常数π出现在将Template:Tsl及其拓扑学联系起来的高斯-博内定理中。具体来说,如果曲面Template:Math高斯曲率K,那么有

ΣKdA=2πχ(Σ)

其中χ(Σ)是该曲面的欧拉示性数,是整数。[106]例如,曲率为1(也就是说其曲率半径也为1,对于球面而言此时的曲率半径与半径重合)的球面S的表面积。球面的欧拉特征数可以通过其同源组计算,其结果为2。于是,便得出

A(S)=S1dA=2π2=4π

即为半径为1的球面的表面积公式。

常数π还出现在拓扑学的许多其他的积分公式中,特别是那些涉及通过陈-韦伊同态的特征类[107]

向量分析

向量分析的方法可以通过分解成球谐函数来理解(图示)

向量分析是与向量場的性质有关的微积分的分支,并有许多物理用途,例如用在电磁学中。位于三维笛卡尔坐标系原点的点源QTemplate:Le[108]

V(𝐱)=kQ|𝐱|

表示位于距原点|𝒙|的单位质量(或电荷)的势能,而k是维度常数。在这里由E表示的场可以是(牛顿)引力場或(库仑)電場,是位势的负梯度

𝐄=V.

特殊情况有库仑定律牛顿万有引力定律高斯定律表明,通过包含原点的任何平滑、简单、封闭、可定向曲面S的场的向外通量等于4πkQ

4πkQ=\oiintS𝐄d𝐀

标准形式会将4π的这因子吸收到常数k中,但这种说法表明了它必须出现在“某处”。此外,4π是单位球面的表面积,但並没有假设S是球面。然而,作为散度定理的结果,由于远离原点的区域是真空(无源的),只有R3{0}中的表面S同调类与计算积分有关,因此可以由相同同调类中的任何方便的表面代替,特别是球形,因为球面坐标可以用于计算积分。

高斯定律的结果之一是位势V的负拉普拉斯算子等于狄拉克δ函数4πkQ倍:

ΔV(𝐱)=4πkQδ(𝐱).

通过卷积就能得到物质(或电荷)的更一般分布,给出泊松方程

ΔV(𝐱)=4πkρ(𝐱)

其中ρ是分布函数。

爱因斯坦方程表明,时空的曲率是由其中的物质能量生成

常数π在与爱因斯坦场方程中的四维势起类似的作用,爱因斯坦方程是形成廣義相對論基础的一條基本公式,并且把引力基本相互作用描述为物质能量引起的时空弯曲的结果:[109]

Rμν12Rgμν+Λgμν=8πGc4Tμν,

Rμv里奇曲率張量R数量曲率gμv度量张量Λ宇宙學常數G万有引力常数c是真空中的光速,而Tμv應力-能量張量。爱因斯坦方程的左边是度量张量的拉普拉斯算子的非线性模拟,並化簡(reduce)至在弱域的極限,而右边是分布函数的模拟乘以8π

柯西积分公式

複雜的解析函數可以以一系列的流綫和等電位綫(許多以直角相交的曲綫)視覺化,圖中是伽瑪函數的複數對數。

复分析中,沿复平面若尔当曲线围道积分是研究解析函数的重要手段之一。简化版的柯西積分公式表明,对任何若尔当曲线γ内任一点z0,以下围道积分给出2πi[110]

γdzzz0=2πi.

该命题是柯西积分定理的直接推论,后者表明上述围道积分在围道的同伦变换下保持不变,因而沿任一曲线的积分和沿以z0为圆心的圆周积分的结果相同。更为一般地,该公式对不通过z0点的任意可求长曲线都成立,但等式右边要乘以曲线关于该点的卷绕数

一般形式的柯西積分公式建立了全纯函数f(z)在若尔当曲线γ上的值与曲线内任意点z0处值的关系:[111][112]

γf(z)zz0dz=2πif(z0)

柯西积分定理是留数定理的一項特例。根据留数定理,在区域内除去有限个解析的亚纯函数g(z)在边界上的围道积分与函数在这些点的留数之和满足:

γg(z)dz=2πiRes(g,ak)

Γ函数与斯特灵公式

Template:Tsl将三塊球面霍普夫纤维化,下方是富比尼-施图迪度量黎曼球面與其富比尼-施图迪度量(如图所示的三塊平行曲面)。恒等式Template:Math可以确定一條数列

阶乘函数n!的值等于所有小于等于n的正整数之积,它的定义域只包含非负整数。Γ函数则是階乘的推广。它在复平面的右半平面定义为:

Γ(z)=0tz1etdt

再利用解析延拓可以将它的定义域扩展到除去非正整数的整塊复数域。当自变量z=n取正整数时,Γ函数给出阶乘(n1)!;当自变量取半整数时,计算结果含有π。例如Γ(1/2)=πΓ(5/2)=3π4[113]

根据魏尔施特拉斯分解定理Γ函数可分解为如下的无穷乘积:[114]

Γ(z)=eγzn=1ez/n1+z/n

γ歐拉-馬斯刻若尼常數。利用该分解公式和Γ函数在z=12的值Γ(12)2=π,亦可以证明沃利斯乘积式。Γ函数和黎曼ζ函數Template:Tsl的恒等式存在关联,其中π扮演着重要的角色

Γ函数常用于计算n维欧氏空间中n 维球的体积和n 维球面的表面积。对n维欧氏空间中半径为rn维球,其体积Vn(r)和表面积Sn1(r)满足:[115]

Vn(r)=πn/2Γ(n2+1)rn
Sn1(r)=nπn/2Γ(n2+1)rn1

两者还满足如下的关系式:

2πr=Sn+1(r)Vn(r).

n很大,用Γ函数可得到阶乘n!的近似公式n!2πn(ne)n,稱斯特靈公式[116],等价于:

π=limne2nn!22n2n+1.

斯特灵近似的几何应用之一是Template:Tsl。将n维欧几里得空间的单纯形记作Δn(n+1)Δn则表示该单纯形的所有面扩大n+1。于是

Vol((n+1)Δn)=(n+1)nn!en+12πn.

这是仅含一點晶格点之凸体体积的(最佳)上界[117]

数论与黎曼ζ函数

全部质数都有其关联的Template:Tsl,即圆的算数定域。分析数论的L函數也定域在每个质数p
基于Template:Tsl的巴塞尔问题的解:ζ(2) 的数值是Template:Tsl中一个基本域的双曲面积的2π倍。

黎曼ζ函数ζ(s) 在数学的许多领域均有应用。当自变量s=2 ,可写作

ζ(2)=112+122+132+

找到这无穷级数的解析解是数学界著名的“巴塞尔问题”。1735年,欧拉解决了这问题,他得到该无穷级数等于π26[74]。欧拉的结论可推导出数论中一項结果,即两随机整数互质(无公因数)的概率为6π2 [4]Template:Rp[118]。整数可由质数p整除的概率為1p2(例如,连续7个正整数只有一个可以7整除),任取两随机整数都能以质数p整除的概率为1p2,至少有一數不能整除的概率则为11p2。又,一随机整数能否以两不同质数整除是相互独立事件,两随机整数互质的概率可以表示成关于所有质数p的无穷乘积[119]

p(11p2)=(p11p2)1=11+122+132+=1ζ(2)=6π261%.

这结论可结合随机数生成器,用統計模擬法π的近似值。[4]Template:Rp

巴塞尔问题的结论意味着几何导出量π的数值与质数的分布有着深刻的关联。巴塞尔问题是谷山-志村定理的一種特殊情况,是Template:Tsl的一項特例,即猜想一个这种形式的算术量关于所有质数p的无穷乘积能够等于一个几何量——某Template:Tsl体积的倒易。巴塞尔问题中,这空间是Template:TslTemplate:Math[120]

ζ函数同样满足黎曼方程的公式,其中用到了π和伽玛公式:

ζ(s)=2sπs1 sin(πs2) Γ(1s) ζ(1s).

除此之外, ζ函数导数也满足

exp(ζ(0))=2π.

最终的结果是π可以从谐振子Template:Tsl中求得。这泛函行列式可以无穷乘积展开式计算,而且这种方法等价于沃利斯乘积公式。[121]这种方法可用于量子力学,尤其是玻尔模型中的变分[122]

傅里叶级数

创建缩略图出错:
Template:Pi出现在P进数中的表示(如图),它們是Template:Tsl的元素。Template:Tsl很大程度地利用了這系統。[123]

周期函数傅里叶级数很自然出现了π。周期函数即实数的小数部分所构成群T=RZ上的函数。傅里叶分解指出,T上的复值函数f可表示为无穷多个TTemplate:Tsl的线性叠加之和。也就是说,T圓群U(1)(模为1的复数组成的乘法群)的映射是连续群同態T的特征都有en(x)=e2πinx的形式,是一條定理。

T有唯一的特征值,直到复共轭,那是一群同态。在圆群用哈尔测度,常数π是这特征值的拉东-尼科迪姆导数值的一半。其他的特征值的导数值为2π的正整数倍。[20]因此,常数π是独特的数字,以至于配备了其哈尔测度的群T,有对于2π整数倍的点阵的庞特里亚金对偶性[124]。这是Template:Tsl的一维版本。

模形式与Template:Theta函数

常數π模形式Θ函數密切相关——比如,椭圆曲线中的Template:Tsl就很大程度涉及楚德诺夫斯基算法(一种快速计算Template:Pi的方法)。

模形式是以在上半平面全純函數的在Template:TslSL2()(或其子群,SL2()SL2()的一格)下的變換特性歸納。Θ函數便是一例:

θ(z,τ)=n=e2πinz+iπn2τ

它是一種名為Template:Tsl的模形式,[125]有時以Template:Tslq=eπiτ表達。

常數π是特殊常數,它會使雅可比Θ函數形成自守形式,即該函數會以特定方式變換。有若干恆等式在所有自守形式下成立。,例如:

θ(z+τ,τ)=eπiτ2πizθ(z,τ)

它使得θ必然在離散海森伯群下以表示(representation)變換。一般模形式和其他Θ函數也包含π,這也是根據Template:Tsl[125]

柯西分布与位势论

箕舌线,英文名来自玛利亚·阿涅西(1718–1799),柯西分布的一幅几何构筑图

积分1x2+1dx=π柯西分布g(x)=1π1x2+1概率密度函数,总概率等于1。

柯西分布的香农熵等于log(4π),也含π

创建缩略图出错:
柯西分布控制做布朗运动的粒子通过膜的通道

柯西分布在位势论中扮演着重要的角色因为它是最简单的Template:Tsl和与在半平面上做布朗运动相关联的经典泊松核[126]Template:Tsl以及希尔伯特变换与泊松核的渐近线有关。希尔伯特变换H是由奇异积分柯西主值给出的积分变换

Hf(t)=1πf(x)xt

常数π是唯一的(正)归一化因子因此H定义了一个在实数轴上的平方可积分实值函数的希尔伯特空间上的Template:Tsl[127]。 和傅里叶变换一样,希尔伯特变换就其在希尔伯特空间L2(R)的变换特性而言可以完全特征化。直到归一化,它是唯一的与正膨胀对易且与实数轴的所有反射反对易有界线性算子[128]。常数π是唯一能使这变换幺正的归一化因子。

複變動態系統

一個複數平面下,曼德博集合的黑色圖案,背景為藍色的
可以從曼德博集合中計算Template:Pi,計算方式和計算從(−0.75, ε)點開始,一直到發散之前的次數有關

大衛·波(David Boll)在1991年發現在曼德博集合分形也有Template:Pi出現[129]。他檢查在曼德博集合在(0.75,0)位置的特性。若考慮坐標在「頸部」(0.75,ε)的點,而ε趨近零,在發散之前迭代的次數和ε相乘,會趨近π。若是在右側尖點處附近的點(0.25,ε)也會有類似的特性:在發散之前迭代的次數和ε的平方根相乘,也會趨近π[129][130]

数学之外的Template:Pi

描述物理现象

π与圆以及球坐标系关系密切,即使π不是物理常数,也常出现在描述宇宙的基本原则方程中。比方说,经典力学领域的简单公式给出长L的单摆小幅摆动的近似周期T2πLgg为地球引力加速度常数。[131]

海森堡不確定性原理是量子力学的基本公式,表明测量粒子时,其位置不确定度(Δx)与动量不确定度(Δp)不可能同时达到任意小(h普朗克常数):[132]ΔxΔph4π

π近似三这特性,和电子偶素的半衰期相對較長有密切的联系。其半衰期的倒數和精细结构常数α的關係為[133]1τ=2π299πmα6m為電子質量。

許多結構工程的公式也有π,例如歐拉推導的挫曲公式說明了長度為L截面二次轴矩Template:Math的細長形物體,在不挫曲的條件下可以承受的最大軸向負載F[134]

F=π2EIL2.

流體動力學斯托克斯定律中也有π。斯托克斯定律是半径约为R的小球體在黏度η流體中以速度v運動時會受到的阻力满足[135]

F=6πηRv.

在理想状态下,河的曲折程度(河道本身的长度与源头到入海口的比值)随着时间的推移逐渐趋向于π。河流外边缘的快速水流弯曲会使河流内边缘加倍侵蚀,河道变得更弯曲,整條河弯折更厉害。然而,这股弯折劲儿最终会导致河流折回一开始弯折的地方,导致“短路”,并形成河迹湖。这两种相反因素使河道长度与源头到入海口的比值的平均值为Template:Pi[136][137]

Template:Pi的记忆技巧

Template:Main article π文字學(或譯作圆周率的语言学)是指記住π的大量位值[4]Template:Rp,并将其世界紀錄載於健力士世界紀錄大全的做法。維爾·美拿(Rajveer Meena)於2015年3月21日在印度於9小時27分鐘內背誦了7萬位的Template:Pi,创下健力士世界紀錄大全認證的世界紀錄。[138]2006年,日本退休工程師原口證在千葉縣於官員見證下背誦了十萬位小數,但他未獲健力士世界紀錄大全認證。[139]

常用於記憶Template:Pi的一項技巧是背誦以單詞長度代表π數值的故事或詩歌:第一單詞有三字母,第二單詞有一字母,第三單詞有四字母,第四單詞有一字母,第五單詞有五字母,如此類推。早期例子是英國科學家詹姆士·金斯設計的詩歌:「Template:Lang[4]Template:Rp這類詩歌有時在英文中稱為「Template:Lang」。除了英文,用於記憶Template:Pi的詩歌亦有不同語言的版本[4]Template:Rp。但是,記憶π的人一般並不以詩歌記憶π來創下紀錄,而是用如記憶數字規律或Template:Tsl的方法。[140]

有好幾位作家仿照上述记忆技巧,用π的數值創作了新型的Template:Tsl方式,當中單詞長度須符合π的數值。《Template:Tsl》以上述技巧包含了π前3835位的值[141],一本標準長度的書《Not a Wake》有一萬單詞,其中各單詞亦代表了π的一位。[142]

大众文化

Pi Pie at Delft University
Template:Pi-{zh-cn:派;zh-hk:批;zh-tw:派;}-;圓形西式餡餅是常見的Template:Pi雙關語(英語圓周率和-{zh-cn:派;zh-hk:批;zh-tw:派;}-同音)

也許因為π的公式很簡短而且四處可見, π比其他數學常数在流行文化中更常見Template:NoteTag

在2008年由英國公開大學英國廣播公司聯合制作的記錄片《Template:Tsl》于2008年十月由英國廣播公司第四台播放。影片讲述了英國數學家马库斯·杜·索托伊在到訪印度研究當地三角學的貢獻時,展示出歷史上Template:Pi最精確的計算公式的信息图形[143]

巴黎的科学博物馆發現宮有間圓形房間叫「π房」,牆上刻有π的707位數,數字貼在圓頂狀的天花板上,由大型的木製字符組成。數值是1853年由英國數學家Template:Tsl計算出來,但是该结果於第528位后開始出現謬誤,在1946年發現,1949年修正。[144][4]Template:Rp

卡尔·萨根的小说《接觸未來》中则暗示说,宇宙的创造者在Template:Pi的数字中暗藏了一则信息。[145]Template:Pi的数字也用在凱特·布希所出的专辑Template:Tsl中的《Pi》的歌词里。[146]

美国人在3月14日庆祝圓周率日,此节日在学生中很流行。[147]一些自称“数学极客”的人常常用π与其数位来创作一些数学或技术圈内人士才能领会到的笑话麻省理工学院则有几句包含“3.14159”的大学Template:Tsl[148]2015年的圆周率日格外重要,按照美式写法,当天的日期时间3/14/15 9:26:53较其他圆周率日包含更多位数的π[149]

北电网络于2011年举行的技术专利拍卖会上,谷歌用了一些包含π在内的数学或科学常数來竞价。[150]

Template:Anchor 在1958年,Template:Tsl提议π换成[[Τ|Template:Tau(tau)]]以便简化公式。τ在此定义为π的兩倍[151]。然而,没有任何其他作者曾这样使用过τ。有些人使用不同的值,τ=2π=6.283185[152]这些人称τ不论是作为弧度制下圆周长的1还是作为弧长与半径的比值(而不是与直径的比值)都比π自然,也能因此简化许多公式。[153][154]有媒体报道称,因为τ的值大小约为6.28,現已有人在6月28日庆祝“τ节”,并吃“两个派”;[155]然而,主流数学界还并未使用τ[156]

1897年,有业余美国数学家试图藉印第安纳州议会来通过後世所謂印第安纳圆周率法案的法案。这法案试图以法律命令强制规定数学常数之值而臭名远播。该法案描述化圆为方的方法,并间接提到了π的错误值,例如3.2。该法案通过了印第安纳州众议院的表决,但参议院否决之。[4]Template:Rp[157][158]

注释

Template:NoteFoot

參考資料

書籍

Template:Refbegin

Template:Refend

引用

Template:Reflist

延伸閱讀

Template:Refbegin

  • Template:Cite book
  • Template:Cite journal
  • Template:Cite journal
  • Chudnovsky, David V. and Chudnovsky, Gregory V., "Approximations and Complex Multiplication According to Ramanujan", in Ramanujan Revisited (G.E. Andrews et al. Eds), Academic Press, 1988, pp 375–396, 468–472
  • Cox, David A., "The Arithmetic-Geometric Mean of Gauss", L' Ensignement Mathematique, 30(1984) 275–330
  • Delahaye, Jean-Paul, "Le Fascinant Nombre Pi", Paris: Bibliothèque Pour la Science (1997) ISBN 2902918259
  • Template:Cite journal
  • Euler, Leonhard, "On the Use of the Discovered Fractions to Sum Infinite Series", in Introduction to Analysis of the Infinite. Book I, translated from the Latin by J. D. Blanton, Springer-Verlag, 1964, pp 137–153
  • Hardy, G. H. and Wright E. M., An Introduction to the Theory of Numbers first published 1938, fifth Edition 1979 with additions 2000, Clarendon Press, Oxford U.K.
  • Heath, T. L., The Works of Archimedes, Cambridge, 1897; reprinted in The Works of Archimedes with The Method of Archimedes, Dover, 1953, pp 91–98
  • Huygens, Christiaan, "De Circuli Magnitudine Inventa", Christiani Hugenii Opera Varia I, Leiden 1724, pp 384–388
  • Template:Cite journal
  • Template:Cite journal
  • Template:Cite journal
  • Niven, Ivan, "A Simple Proof that pi Is Irrational", Bulletin of the American Mathematical Society, 53:7 (July 1947), 507
  • Ramanujan, Srinivasa, "Modular Equations and Approximations to π", Quarterly Journal of Pure and Applied Mathematics, XLV, 1914, 350–372. Reprinted in G.H. Hardy, P.V. Seshu Aiyar, and B. M. Wilson (eds), Srinivasa Ramanujan: Collected Papers, 1927 (reprinted 2000), pp 23–29
  • Template:Tsl, Contributions to Mathematics Template:Sic Chiefly of the Rectification of the Circle to 607 Places of Decimals, 1853, pp. i–xvi, 10
  • Template:Cite journal
  • Tropfke, Johannes, Geschichte Der Elementar-Mathematik in Systematischer Darstellung (The history of elementary mathematics), BiblioBazaar, 2009 (reprint), ISBN 978-1-113-08573-3
  • Viete, Francois, Variorum de Rebus Mathematicis Reponsorum Liber VII. F. Viete, Opera Mathematica (reprint), Georg Olms Verlag, 1970, pp 398–401, 436–446
  • Template:Tsl, "Is Pi Normal?", The Mathematical Intelligencer, 7:3(1985) 65–67
  • Wallis, John, Arithmetica Infinitorum, sive Nova Methodus Inquirendi in Curvilineorum Quadratum, aliaque difficiliora Matheseos Problemata, Oxford 1655–6. Reprinted in vol. 1 (pp 357–478) of Opera Mathematica, Oxford 1693
  • Zebrowski, Ernest, A History of the Circle: Mathematical Reasoning and the Physical Universe, Rutgers University Press, 1999, ISBN 978-0-8135-2898-4

Template:Refend

外部連結

Template:Commons category

Template:- Template:幾何術語 Template:無理數導航

Template:Authority control Template:Good article

  1. Template:Cite book
  2. Template:Cite journal
  3. Template:Cite news
  4. 4.00 4.01 4.02 4.03 4.04 4.05 4.06 4.07 4.08 4.09 4.10 4.11 4.12 4.13 4.14 4.15 4.16 4.17 4.18 4.19 4.20 4.21 4.22 4.23 4.24 4.25 4.26 4.27 4.28 4.29 4.30 4.31 4.32 4.33 4.34 4.35 4.36 4.37 4.38 4.39 4.40 4.41 4.42 4.43 4.44 4.45 4.46 4.47 4.48 4.49 4.50 4.51 4.52 4.53 4.54 4.55 4.56 4.57 4.58 4.59 4.60 4.61 4.62 4.63 4.64 4.65 4.66 4.67 4.68 4.69 4.70 4.71 4.72 4.73 4.74 4.75 4.76 4.77 4.78 Template:Harvnb
  5. Template:Citation
  6. 6.0 6.1 Template:Harvnb
  7. Template:Cite web
  8. Template:Cite web
  9. Template:Cite journal
  10. Template:Cite web
  11. Template:Citation。 Page 102: "From a logical point of view, this is unsatisfactory at the present stage because we have not yet discussed the concept of arc length." Arc length is introduced on page 529.
  12. 12.0 12.1 12.2 Template:Citation
  13. Template:Harvtxt。魏尔斯特拉斯使用的积分实为π=dx1+x2.
  14. Template:Citation
  15. Template:Citation
  16. 16.0 16.1 Template:Cite book, p 183.
  17. Template:Cite book, p 2.
  18. Template:Citation
  19. Template:Citation, §VIII.2
  20. 20.0 20.1 Template:Citation, §II.3.
  21. Template:Cite journal
  22. 22.0 22.1 Template:Cite news
  23. Template:Cite web
  24. Template:Harvnb
  25. Template:Harvnb
  26. Template:Harvnb
    Template:Cite book,第185页
  27. 27.0 27.1 Template:Harvnb
  28. Template:SloanesRef Retrieved 12 April 2012.
  29. Template:Cite journal
  30. Template:Citation克劳狄乌斯·托勒密使用了60进制下的三位小数去近似,随后卡西将其扩展到了九位小数。参见Template:Citation
  31. 引用错误:<ref>标签无效;未给name(名称)为base 60 OEIS的ref(参考)提供文本
  32. Template:Harvnb
  33. 33.0 33.1 Template:Harvnb
  34. Maor, Eli, E: The Story of a Number, Princeton University Press, 2009, p 160, ISBN 978-0-691-14134-3 ("five most important" constants)
  35. Template:Cite mathworld
  36. Template:Citation
  37. Template:Citation
  38. Template:Citation
  39. Template:Citation, Chapter 7.
  40. Template:Citation
  41. Template:Citation
  42. Template:Citation
  43. Template:Citation
  44. Template:Harvnb
  45. Template:Harvnb
  46. Feller, W. An Introduction to Probability Theory and Its Applications, Vol. 1, Wiley, 1968, pp 174–190.
  47. Template:Citation; Section 2.7
  48. Template:Citation; Theorem 1.13.
  49. Petrie, W.M.F. Wisdom of the Egyptians (1940)
  50. Rossi, Corinna Architecture and Mathematics in Ancient Egypt, Cambridge University Press. 2007. ISBN 978-0-521-69053-9.
  51. Legon, J. A. R. On Pyramid Dimensions and Proportions (1991) Discussions in Egyptology (20) 第25-34页 [1] Template:Wayback
  52. "We can conclude that although the ancient Egyptians could not precisely define the value of Template:Pi, in practice they used it". Template:Cite journal,第70页。
    Template:Cite journal,第30页。
    参见Template:Cite journal.
    参见Template:Cite journal
  53. Egyptologist: Rossi, Corinna, Architecture and Mathematics in Ancient Egypt, Cambridge University Press, 2004, 第60-70页,200, ISBN 9780521829540.
    对此的怀疑:Template:Tsl, The Skeptic Encyclopedia of Pseudoscience, ABC-CLIO, 2002, 第407-408页, ISBN 9781576076538.
    参见Fagan, Garrett G., Archaeological Fantasies: How Pseudoarchaeology Misrepresents The Past and Misleads the Public, Routledge, 2006, ISBN 9780415305938.
    若需要一个没有Template:Pi参与的解释的列表,请见Template:Cite book
  54. Chaitanya, Krishna. A profile of Indian culture. Template:Wayback Indian Book Company (1975). 第133页。
  55. Template:Cite web
  56. Template:Harvnb
  57. Template:HarvnbGrienberger在1630年已計算到39位小數,Sharp在1699年計算到71位小數
  58. 周髀算經》注中, 趙爽指出「圓徑一而周三,方徑一而匝四」。
  59. 59.0 59.1 Template:Harvnb
  60. Template:Cite book
  61. Template:Cite journal
  62. Template:Cite web
  63. 63.0 63.1 63.2 Template:Harvnb
  64. Template:Cite book。其計算結果是3.14159 26535 89793 23846 26433 83279 50288 4196 < Template:Pi < 3.14159 26535 89793 23846 26433 83279 50288 4199
  65. Template:Harvnb
  66. 66.0 66.1 66.2 Template:Harvnb
  67. Template:Harvnb
  68. 68.0 68.1 Template:Harvnb。牛顿由Arndt引用。
  69. Template:OEIS2C
  70. 70.0 70.1 Template:Harvnb
  71. 71.0 71.1 Template:Cite journal
  72. Template:Cite journal
  73. Template:Harvnb, (formula 16.10). Note that (n − 1)n(n + 1) = n3 − n.
    Template:Cite book
  74. 74.0 74.1 Template:Harvnb
  75. Lambert, Johann, "Mémoire sur quelques propriétés remarquables des quantités transcendantes circulaires et logarithmiques", reprinted in Template:Harvnb
  76. Hardy and Wright 1938 and 2000:177 footnote §11.13-14 references Lindemann's proof as appearing at Math. Ann. 20 (1882), 213-25.
  77. cf Hardy and Wright 1938 and 2000:177 footnote §11.13-14. The proofs that e and π are transcendental can be found on pages 170-176. They cite two sources of the proofs at Landau 1927 or Perron 1910; see the "List of Books" at pages 417-419 for full citations.
  78. Template:Harvnb. A facsimile of Jones' text is in Template:Harvnb
  79. 參考Template:Harvnb威廉·奥特雷德用字母Template:Pi來表示一個圓的周長
  80. Template:Harvnb。参见 Template:Harvnb.
  81. Template:Harvnb.
    具体算法情参见Template:Harvnb
  82. 82.0 82.1 82.2 Template:Cite web
  83. Template:Harvnb“39 digits of π are sufficient to calculate the volume of the universe to the nearest atom.”(中文:39位π已足以计算宇宙到其最近的原子的体积了。)
    鉴于人们需要一些额外的数位来解决捨入誤差的问题,阿恩特称几百位小数足以应付任何科学计算。
  84. Template:Cite news
  85. Template:Cite news
  86. Template:Harvnb
  87. Template:Cite web
  88. 88.0 88.1 Template:Cite web
  89. PSLQ即Partial Sum of Least Squares,偏最小二乘和。
  90. Template:Cite web
  91. Template:Citation
  92. Template:Citation
  93. Template:Harvnb
  94. Template:Cite mathworld
  95. Template:Citation
  96. Template:Harvnb
  97. 97.0 97.1 Gibbons, Jeremy, "Unbounded Spigot Algorithms for the Digits of Pi" Template:Wayback, 2005. Gibbons produced an improved version of Wagon's algorithm.
  98. Template:Cite journal A computer program has been created that implements Wagon's spigot algorithm in only 120 characters of software.
  99. Template:Cite journal
  100. Template:Harvnb。普勞夫有找到十進制的位數萃取演算法,但其速度比完整計算之前所有位數要慢。
  101. Bellards formula in: Template:Cite web
  102. Template:Cite news
  103. Template:Cite mathworld
  104. 104.0 104.1 Template:Harvnb
  105. 105.0 105.1 Template:Harvnb
  106. Template:Citation; Chapter 6.
  107. Template:Citation; Chapter XII Characteristic classes
  108. H. M. Schey (1996) Div, Grad, Curl, and All That: An Informal Text on Vector Calculus, ISBN 0-393-96997-5.
  109. Yeo, Adrian, The pleasures of pi, e and other interesting numbers, World Scientific Pub., 2006, p 21, ISBN 978-981-270-078-0.
    Ehlers, Jürgen, Einstein's Field Equations and Their Physical Implications, Springer, 2000, p 7, ISBN 978-3-540-67073-5.
  110. Template:Citation
  111. Template:Cite mathworld
  112. Joglekar, S. D., Mathematical Physics, Universities Press, 2005, p 166, ISBN 978-81-7371-422-1.
  113. Template:Harvnb
  114. Template:Citation
  115. Template:Citation
  116. Template:Harvnb
  117. Template:Citation
  118. 此理论由Ernesto Cesàro于1881年证明. For a more rigorous proof than the intuitive and informal one given here, see Hardy, G. H., An Introduction to the Theory of Numbers, Oxford University Press, 2008, ISBN 978-0-19-921986-5, theorem 332.
  119. Ogilvy, C. S.; Anderson, J. T., Excursions in Number Theory, Dover Publications Inc., 1988, pp. 29–35, ISBN 0-486-25778-9.
  120. Template:Citation
  121. Template:Citation
  122. Template:Cite journal
  123. Template:Citation
  124. Template:Citation; Chapter 4
  125. 125.0 125.1 Template:Citation
  126. Template:Citation
  127. * Template:Citation.
  128. Template:Citation; Chapter II.
  129. 129.0 129.1 Template:Cite journal
  130. Peitgen, Heinz-Otto, Chaos and fractals: new frontiers of science, Springer, 2004, pp. 801–803, ISBN 978-0-387-20229-7.
  131. Halliday, David; Resnick, Robert; Walker, Jearl, Fundamentals of Physics, 5th Ed., John Wiley & Sons, 1997, p 381, ISBN 0-471-14854-7.
  132. Template:Cite web
  133. C. Itzykson, J-B. Zuber, Quantum Field Theory, McGraw-Hill, 1980.
  134. Low, Peter, Classical Theory of Structures Based on the Differential Equation, CUP Archive, 1971, pp 116–118, ISBN 978-0-521-08089-7.
  135. Batchelor, G. K., An Introduction to Fluid Dynamics, Cambridge University Press, 1967, p 233, ISBN 0-521-66396-2.
  136. Template:Cite journal
  137. Template:Harvnb
  138. Template:Cite web
  139. Template:Cite news
  140. Template:Cite journal
  141. Template:Cite web
  142. Template:Cite book
  143. Template:Citation, showing a visualization of the historically first exact formula, starting at 35 min and 20 sec into the second part of the documentary.
  144. Template:Harvnb
  145. Template:Harvnb这部分情节在改编的电影《超時空接觸》中省略。
  146. Template:Cite journal
  147. Template:Cite web
  148. Template:Cite web
  149. Template:Cite web
  150. Template:Cite news
  151. Template:Cite book
  152. Sequence Template:OEIS2C
  153. Template:Cite journal
  154. Template:Cite journal
  155. Template:Cite web
  156. Template:Cite journal
  157. Template:Harvnb
  158. Template:Cite journal