圆周率

圆周率是数学常数，等于任何圆的周长和其直径的比，一个常见的近似值等于，常用符号表示。

是无理数，不能用分数表示出来（即它的小数部分是无限不循环小数），但近似$\frac{22}{7}$等有理数。学界认为的数字序列在统计上是随机分布，但迄今未能证明。此外，还是超越数，亦即它不是任何有理系数多项式的根；化圆为方的问题不可能用尺规作图解决。

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

定义涉及圆，在三角学和几何学的许多公式，特别是广泛应用在圆形、球形或椭球形相关公式中。在近代数学分析里，改由实数系统谱性质中的特征值或周期来定义，其他数学领域如数论、统计以及几乎所有物理学领域均有出现，的广泛用途使它成为科学界内外最广为人知的数学常数。几本专门介绍的书籍经已出版，圆周率日（3月14日）和值计算突破记录也往往会成为报纸的新闻头条。此外，背诵值的世界记录已达10万位。

<image>

## 基本概念

### 名称

数学家用小写希腊字母表示圆周和其直径之比，有时也将其拼写为“Pi”，来自希腊语“περίμετρος”（周长）的首字母。英语的发音与英文单词“Pie”（p，西式馅饼）相同。π的小写字母（或其无衬线体）在数学要和表示求积符号的大写Π相区分开。

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

### 定义

<image>

常用定义为圆的周长$C$与直径$d$的比值：

$\pi = \frac{C}{d}$。

无论圆的大小如何，比值$\frac{C}{d}$为恒值。如果圆的直径变为原先的二倍，周长也变为二倍，比值$\frac{C}{d}$不变。目前的定义暗地用了欧几里得几何的一些定理，虽然圆的定义可扩展到任意曲面（即非欧几里得几何），但这些圆不符合定律$\pi = \frac{C}{d}$。

这里，圆的周长指其圆周的弧长，弧长这概念可以不依赖几何学，而是用微积分学的极限来定义。例如，若想计算笛卡儿坐标系中单位圆$x^2+y^2=1$上半部分的弧长，需要用到积分：

$$\pi = \int_{-1}^1 \frac{dx}{\sqrt{1-x^2}}.$$

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

$\pi=\int_{-\infty}^\infty\frac{dx}{1+x^2}.$

这些依赖周长、且暗地依赖积分的定义如今在文献中并不常见。雷默特解释说现代教微积分时，大学一般将微分学课程安排在积分学课程之前，所以不依赖于后者的π的定义就很有必要了。其中一种定义由理查·巴尔策提出，由爱德蒙·兰道推广，其表述如下：是两倍于能使余弦函数等于零的最小正数。余弦函数可以由独立于几何之外的幂级数定义，或者使用微分方程的解来定义。

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

$\{\dots,-2\pi i, 0, 2\pi i, 4\pi i,\dots\} = \{2\pi ki| k\in\mathbb Z\}$。

并且其中包括独特的正实数。

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

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

### 无理及正规性

是无理数，无法表示成两整数之比的形式（形如$\frac{22}{7}$的分数常用来近似表达，但是没有任何普通分数（指整数的比）可以取到的精确值）。由于是无理数，故可表示为无限不循环小数。有多种方法能证明π是无理数，这些证明也都要用到微积分学和反证法。可以用有理数来近似的程度还无法准确得知（称为无理性度量），不过估计其无理性度量比$e$或$\ln(2)$的要大，但是小于刘维尔数的无理性度量。

统计随机性检验，包括正规数检验，可验证的位数没有明显的固定模式。小数中任意固定长度的序列（如3位数000，001……999）出现几率都相同。不过有关是正规数的猜想既无证明，亦无证伪。

电脑出现后可生成大量的不同位数，并统计分析之。金田康正详细统计分析了的十进制数字，并验证了其分布正规：例如，假设检定0到9十个数的出现频率，找不到有特定重复规律的证据。根据无限猴子定理，任何任意长度、由随机内容组成的子序列看起来都有可能像不随机生成。因此，就算的小数序列通过了随机性统计测试，其中也可能有几位的数字看起来似有规律可循而非随机数，例如的十进制写法在小数第762位后开始出现了连续六个9。

### 超越性

<image>

不仅是无理数，还是超越数，即不是任何有理系数多项式的根。（比方说，试图解有限项方程$\frac{x^5}{120} - \frac{x^3}{6} + x = 0$来求的值）

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

### 连分式

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

$$\pi=3+\textstyle \frac{1}{7+\textstyle \frac{1}{15+\textstyle \frac{1}{1+\textstyle \frac{1}{292+\textstyle \frac{1}{1+\textstyle \frac{1}{1+\textstyle \frac{1}{1+\ddots}}}}}}}$$

在这连分数的任意一点截断化简，都能得到π的近似值；前四位近似值是3、$\frac{22}{7}$、$\frac{333}{106}$、$\frac{355}{113}$。这些数在历史上是最广为人知且广为使用的几个近似值。用以上方式得出的的近似值要比任何有相同或更小的整数分母的其他整数分数近似值更接近。是超越数，据定义来说它不是代数数，又因此不可能是二次无理数；是故不能表示为循环连分数。尽管的简单连分数没有表现出任何其他明显规律，数学家发现了数条广义连分数能表示，例如：

$$\pi=\textstyle \cfrac{4}{1+\textstyle \frac{1^2}{2+\textstyle \frac{3^2}{2+\textstyle \frac{5^2}{2+\textstyle \frac{7^2}{2+\textstyle \frac{9^2}{2+\ddots}}}}}} =3+\textstyle \frac{1^2}{6+\textstyle \frac{3^2}{6+\textstyle \frac{5^2}{6+\textstyle \frac{7^2}{6+\textstyle \frac{9^2}{6+\ddots}}}}} =\textstyle \cfrac{4}{1+\textstyle \frac{1^2}{3+\textstyle \frac{2^2}{5+\textstyle \frac{3^2}{7+\textstyle \frac{4^2}{9+\ddots}}}}}$$

### 近似值

圆周率近似值包括：

- 整数：3
- 分数（依准确度顺序排列）：$\frac{13}{4}$、$\frac{16}{5}$、$\frac{19}{6}$、$\frac{22}{7}$、$\frac{179}{57}$、$\frac{267}{85}$、$\frac{333}{106}$、$\frac{355}{113}$、$\frac{52163}{16604}$、$\frac{53228}{16943}$、$\frac{55358}{17621}$、$\frac{57843}{18412}$、$\frac{60328}{19203}$、$\frac{103993}{33102}$、$\frac{245850922}{78256779}$（选自 及。）
- 小数（整数后首100位）：3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 06286 20899 86280 34825 34211 70679...（另见）

其他进位制的近似值

- 十六进制（整数后首20位）：3.243F
- 六十进制（整数后首20位）：60…

### 复数与欧拉恒等式

<image>

任何复数（以$z$为例）都可以表示为一组实数对：极坐标系用实数$r$表示半径，代表复平面上复数$z$离原点的距离；实数$\varphi$则表示夹角，即这条半径（复平面上复数$z$与原点的连线）与正实轴经顺时针转动的夹角。这样一来，$z$就可写成

$z = r\cdot(\cos\varphi + i\sin\varphi)$，这里$i$代表虚数单位，即$i^2=-1$。

复分析中，欧拉公式将三角函数与复指数函数糅合在一起：

$e^{i\varphi} = \cos \varphi + i\sin \varphi$，这里数学常数e是自然对数的底数。

欧拉公式确立了$e$的复指数与复平面上以原点为圆心的单位圆上的点之间的关系，而且当$\varphi=\pi$时，欧拉公式就能改写为欧拉恒等式的形式：

e^(i⋅pi)+1。此等式亦称“最奇妙的数学公式”（the most remarkable formula in mathematics），全因它将五个最基本的数学常数简洁联系起来。

欧拉等式亦可用于求出方程$z^n=1$的$n$个不同复数根（这些根叫做$n$次单位根”），可以根据以下公式求得：

$e^{ \frac{2 \pi i k}{n}} \qquad (k = 0, 1, 2, \dots, n - 1)$。

### 谱特征

<image>

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

在许多用处中都会以特征值形式出现。例如理想的振动弦问题可以建模为函数$f$在单位区间$[0,1]$的图形，固定边界值为$f(0)=f(1)=0$。弦振动的模态会是微分方程的$f^n(x)+\lambda^2f(x)=0$，此处λ是相关的特征值。受施图姆－刘维尔理论限制，$\lambda$只能是一些特定的数值。而$\lambda=\pi$即为一个特征值，因为函数$f(x)=\sin (\pi x)$满足边界条件及微分方程$\lambda=\pi$。

<image>

是上述方程的最小特征值，也和弦振动的基本模式有关。一种让弦振动的方式是提供弦能量，能量会满足维廷格函数不等式，其中提到若函数$f:[0,1]\rightarrow \mathbb{C}$使得$f(0)=f(1)=0$，且$f$和$f'$都是平方可积函数，则以下的不等式成立：

$$\pi^2\int_0^1|f(x)|^2\,dx\le \int_0^1|f'(x)|^2\,dx,$$

此例中等号成立的条件恰好是$f$为$\sin (\pi x)$倍数的时候。因此似乎是维尔丁格不等式的最佳常数，也是最小的特征值（根据雷利商数的计算方式）

在更高维度的分析也有类似的角色，出现在其他类似问题的特征值中。就如以上所述，的一项特点是等周定理中的最佳常数：周长为$P$的平面若尔当曲线，所围面积$A$满足以下的不等式

$4\pi A\le P^2$。

$A=\pi r^2$及$P=2\pi r$，故等号成立的条件是曲线为圆形。

圆周率也和庞加莱不等式的最佳常数有关，是一维及二维的狄氏能量特征向量最佳值中最小，会出现在许多经典的物理现象中，例如经典的位势论。其一维的情形即为维廷格不等式。

圆周率也是傅里叶变换的重要常数，傅里叶变换属于积分变换，将实数线上有复数值、可积分的函数，转换为以下形式：

$$\hat{f}(\xi) = \int_{-\infty}^\infty f(x) e^{-2\pi i x\xi}\,dx.$$

傅里叶变换有几种不同的写法，但不论怎么写，傅里叶变换及反傅里叶变换中，一定会有某处出现。不过上述的定义是最经典的，因为其描述了L2空间中唯一的幺正算符，也是$L^1$空间到$L^\infty$空间的代数同态。

不确定性原理也用到。不确定性原理提出了可以将函数在空间及在频域中局部化程度的下限，用傅立叶转换的方式表示：

$\int_{-\infty}^\infty x^2|f(x)|^2\,dx\ \int_{-\infty}^\infty \xi^2|\hat{f}(\xi)|^2\,d\xi\ge \left(\frac{1}{4\pi}\int_{-\infty}^\infty |f(x)|^2\,dx\right)^2$。

物理的结果，有关量子力学中同时观测位置及动量的不确定性，见下文。傅立叶分析中出现是史东－凡纽曼定理的结果，证实了海森伯群的薛定谔表示是唯一。

### 高斯积分

<image>

高斯积分是对高斯函数$e^{-x^2}$在整条实轴上的积分，即函数下方与X轴围成的面积，其结果为$\sqrt{\pi}$。

$$\int_{-\infty}^\infty e^{-x^2} \, dx=\sqrt{\pi}$$

此积分的计算可以先计算$f(x)=e^{-x^2}$对整条实轴的积分的平方，通过转换笛卡尔坐标系为极坐标系从而求得

$$\left( \int_{-\infty}^\infty e^{-x^2} \, dx \right)^2 = \iint_{\mathbf{R}^2} e^{-(x^2+y^2)}\,dxdy = \int_0^{2\pi} \int_0^{\infty} e^{-r^2}r\,dr\,d\theta = \pi$$

其他计算方法可参阅高斯积分。高斯函数更一般的形式为$f(x)=a\exp{\frac{-(x-b)^2}{2c^2}}$，求一般形式的高斯积分均可通过换元积分法转化为求$f(x)=e^{-x^2}$的积分。

另外，当高斯函数为以下形式时，它则是平均数为$\mu$和标准差为$\sigma$的正态分布的几率密度函数：

$$f(x) = {1 \over \sigma\sqrt{2\pi} }\,\exp{\frac{-(x-\mu )^2}{2\sigma^2}}$$

这函数是概率密度函数，函数下方与X轴围成的面积必须为1，令$\mu=0$和$\sigma=1$即可变换得出$\int_{-\infty}^\infty e^{-x^2} \, dx =\sqrt{\pi}$。概率论与统计学领域经常使用正态分布来作为复杂现象的简单模型：例如科学家通常假设大多数试验观测值的随机误差都是服从正态分布。

<image>

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

$\sigma_{x}\sigma_{p} \geq \frac{\hbar}{2}$。

这里的$\sigma_{x}$与$\sigma_{p}$分别为位置与动量的标准差，$\hbar$是约化普朗克常数，而不等式的等号当且仅当粒子的波函数为高斯函数使成立。

同样地，作为唯一独特的常数使得高斯函数等于其自身的傅里叶变换，此时的高斯函数形式为$f(x)=e^{-\pi x^2}$。根据豪（Howe）的说法，建立傅里叶分析基本定理的“全部工作（whole business）”简化为高斯积分。

## 历史

### 远古时期

圆周率在远古时期（公元前一千纪）已估算至前两位（3.1）。有些埃及学家声称，远至古王国时期时期的古埃及人已经用$\frac{22}{7}$作为圆周率的约数，但这说法受到质疑。

最早有记载的对圆周率估值在古埃及和巴比伦出现，两估值都与圆周率的正确数值相差不到百分之一。巴比伦曾出土一块公元前1900至1600年的泥板，泥板上的几何学陈述暗示人们当时把圆周率视同$\frac{25}{8}$（等于3.125）。埃及的莱因德数学纸草书（鉴定撰写年份为公元前1650年，但抄自一份公元前1850年的文本）载有用作计算圆面积的公式，该公式中圆周率等于$(\frac{16}{9})^2$（≈3.1605）。

公元前4世纪的《百道梵书》的天文学运算把$\frac{339}{108}$（≈3.139，精确到99.91%）用作圆周率估值。公元前150年前其他印度文献把圆周率视为$\sqrt{10}$（≈3.1622）。

### 割圆时代

<image>

第一条有纪录、严谨计算数值的算法是用正多边形的几何算法，在公元前250年由希腊数学家阿基米德发明。这算法用了有一千年之久，因而有时亦称阿基米德常数。阿基米德的算法是在计算圆的外切正六边形及内接正六边形的边长，以此计算的上限及下限，之后再将六边形变成十二边形，继续计算边长，一直计到正96边形为止。他根据多边形的边长证明$\frac{223}{71}<\pi<\frac{22}{7}$（也就是$3.1408<\pi<3.1429$）。阿基米德得到的上限$\frac{22}{7}$也造成常见误解，认为就等于$\frac{22}{7}$。在公元前150年，希腊罗马的科学家克劳狄乌斯·托勒密在《天文学大成》一书中提到的数值是3.1416，可能来自阿基米德，也可能来自阿波罗尼奥斯。数学家在1630年利用多边形的方式计算到第39位小数，一直到1699年，其他数学家才利用无穷级数的方式打破其纪录，计算到第71位小数。

<image>

中国历史上，的数值有3、3.1547（公元前一世纪）、$\sqrt{10}$（公元前100年，数值约3.1623）及$\frac{142}{45}$（第三世纪，数值约3.1556）。大约在公元265年，曹魏数学家刘徽创立割圆术，用正3072边形计算出Π的数值为3.1416。他后来又发明了较快的算法，利用边数差两倍的正多边形，其面积的差值会形成等比数列，其公比为$\frac{1}{4}$的原理，配合96边形算出的值为3.14。祖冲之在公元480年利用割圆术计算12288边形边长，得到的值在3.1415926和3.1415927之间。他同时提出了π的约率$\frac{22}{7}$和密率$\frac{355}{113}$。在之后的八百年内，这都是Π最准确的估计值。178为纪念祖冲之对圆周率发展的贡献，日本数学家三上义夫将这推算值命名为“祖冲之圆周率”，简称“祖率”。

印度天文学家阿耶波多在公元499年的著作《阿里亚哈塔历书》中使用了3.1416的数值。斐波那契在大约1220年用独立于阿基米德多边形法，计算出3.1418。意大利作家但丁·阿利吉耶里用的数值则是$3+\frac{\sqrt{2}}{10}\approx 3.14142$。

波斯天文学家卡西在1424年利用3×228边的多边形，计算到六十进制的第9位小数，相当十进制的第16位小数。这一突破成为当时的纪录，延续了约180年。法国数学家弗朗索瓦·韦达在1579年用3×217边形计算到第9位小数，佛兰芒数学家阿德里安·范·罗门在1593年计算到第15位小数。荷兰数学家鲁道夫·范·科伊伦在1596年计算到第20位小数，他之后又计算到第35位小数（因此在二十世纪初之前，圆周率在德国会称为鲁道夫数）。荷兰科学家威理博·司乃耳在1621年计算到第34位小数，而奥地利天文学家克里斯托夫·格林伯格在1630年用1040边形计算到第38位小数，至今这仍是利用多边形算法可以达到最准确的结果。

### 无穷级数

16及17世纪时，开始改用无穷级数的方式去计。无穷级数是一组无穷数列的和。无穷级数让数学家可以计算出比阿基米德以及其他用几何方式计算的数学家更准确的结果。虽然詹姆斯·格雷果里及戈特弗里德·莱布尼茨等欧洲数学家利用无穷数列计算Π而使得该方法为大家所知，但这种方法最早是由印度科学家在大约1400到1500年之间发现。第一个记载用无穷级数计算的人是约公元1500年左右时，印度天文学家尼拉卡莎·萨默亚士在他的著作《系统汇编》中用梵语诗所记录。当时没有这数列对应的证明，而证明出现在另一本较晚的印度作品《基本原理》，年代约在公元1530年。尼拉卡莎将该数列归功于更早期的印度数学家桑加马格拉马的马德哈瓦（1350–1425）。相关的无穷级数有许多，包括有关$\sin$、$\tan$及$\cos$的，现在称为马德哈瓦数列或π的莱布尼茨公式。玛达瓦在1400年用无穷级数计算到第11位小数，但在1430年一位波斯数学家卡西利用多边形算法否定了他算的结果。

<image>

欧洲发现的第一条无穷项圆周率公式是无穷乘积（和一般用来计算的无穷级数不同），由法国科学家弗朗索瓦·韦达在1593年发现：

$$\frac2\pi = \frac{\sqrt2}2 \cdot \frac{\sqrt{2+\sqrt2}}2 \cdot \frac{\sqrt{2+\sqrt{2+\sqrt2}}}2 \cdots$$

约翰·沃利斯在1655年发现了沃利斯乘积，是欧洲发现的第二条无穷项圆周率公式：

$$\frac{\pi}{2} = \frac{2}{1} \cdot \frac{2}{3} \cdot \frac{4}{3} \cdot \frac{4}{5} \cdot \frac{6}{5} \cdot \frac{6}{7} \cdot \frac{8}{7} \cdot \frac{8}{9} \cdots$$

微积分学由英国科学家艾萨克·牛顿及德国数学家戈特弗里德·莱布尼茨在1660年代发明，许多计的无穷级数出现。牛顿自己就用反正弦（$\arcsin$）数列在1655年或1666年将近似到第15位小数，后来写到“我很羞愧告诉你我为了计算它用了多少数字，我当时没有做其他事。”

苏格兰数学家詹姆斯·格雷果里在1671年发现了马德哈瓦公式，莱布尼茨也在1674年发现：

$$\arctan z = z - \frac {z^3} {3} +\frac {z^5} {5} -\frac {z^7} {7} +\cdots$$

这公式即为格雷果里－莱布尼茨公式，在$z=1$时数值为$\frac{\pi}{4}$。1699年时英国数学家亚伯拉罕·夏普用格雷果里－莱布尼茨公式，在$z=\frac{1}{\sqrt{3}}$时计算，计算到π的第71位小数，打破由多边形算法得到的第39位小数的记录。格雷果里－莱布尼茨公式在$z=1$时非常简单，但收敛到最终值的速度非常慢，现在不会再用此公式来计π。

约翰·梅钦在1706年用格雷果里－莱布尼茨级数产生了可以快速收敛的公式：

$$\frac{\pi}{4} = 4 \, \arctan \frac{1}{5} - \arctan \frac{1}{239}$$

梅钦用这公式计到第100位小数后来其他数学家也发展了一些类似公式，现在称为梅钦类公式，创下了许多计算位数的纪录。在进入电脑时代时，梅钦类公式仍然是耳熟能详可以计算的公式，而且在约250年的时间里，很多有关位数的纪录都是梅钦类公式所得，比如在1946年时由达尼尔·弗格森（Daniel Ferguson）用这类公式计到第620位小数，是没有计算设备辅助的最佳纪录。

1844年，计算天才扎卡里亚斯·达斯在德国数学家卡尔·弗里德里希·高斯的要求下以梅钦类公式心算了的200位小数，并创下纪录。英国数学家威廉·谢克斯花了15年的时间计算到小数707位，不过第528位小数出错，后面的小数也都不正确。

#### 收敛速度

有些的无穷级数收敛的比其他级数要快，数学家一般会选用收敛速度较快的级数，可以在较少的计算量下计算，且达到需要的准确度。以下是的莱布尼茨公式：

$$\pi = \frac{4}{1} - \frac{4}{3} + \frac{4}{5} - \frac{4}{7} + \frac{4}{9} - \frac{4}{11} + \frac{4}{13} - \cdots$$

随着一项一项的值加入总和中，只要项次够多，总和最后会慢慢接近π。不过此数列的收敛速度很慢，要到50万项之后，才会精确到的第五位小数。

尼拉卡莎在15世纪发展了的另一条无穷级数，收敛速度比格雷果里－莱布尼茨公式快很多：

$$\pi = 3 + \frac{4}{2\times3\times4} - \frac{4}{4\times5\times6} + \frac{4}{6\times7\times8} - \frac{4}{8\times9\times10} + \cdots$$

以下比较两条级数的收敛速率：

无穷级数｜第1项｜前2项｜前3项｜前4项｜前5项｜收敛到

$\pi = \frac{4}{1} - \frac{4}{3} + \frac{4}{5} - \frac{4}{7} + \frac{4}{9} - \frac{4}{11} + \frac{4}{13} \cdots.$｜4.0000｜2.6666…｜3.4666…｜2.8952…｜3.3396…｜3.1415…

$\pi = {{3}} + \frac{{4}}{2\times3\times4} - \frac{{4}}{4\times5\times6} + \frac{{4}}{6\times7\times8} \cdots.$｜3.0000｜3.1666…｜3.1333…｜3.1452…｜3.1396…

计算前五项后，格雷果里－莱布尼茨级数的和跟的误差为0.2，而尼拉卡莎级数和的误差为0.002。尼拉卡莎级数收敛快很多，也甚为适合用来计的值。收敛更快的级数有梅钦类公式及楚德诺夫斯基算法，后者每计一项就可以得到14位正确的小数位。

### 无理与超越性

并非所有和有关的研究都旨在提高计算它的准确度。1735年，欧拉解决了巴塞尔问题，建立了所有平方数倒数和与的关系。之后欧拉发现了欧拉乘积公式，得到了、素数的重要关联，对日后黎曼ζ函数的研究影响深远。

$$\frac{\pi^2}{6} = \frac{1}{1^2} + \frac{1}{2^2} + \frac{1}{3^2} + \frac{1}{4^2} + \cdots$$

1761年，瑞士数学家约翰·海因里希·朗伯用正切函数的无穷连分数表达式证明了是无理数。1794年，法国数学家阿德里安－马里·勒让德证明了$\pi^2$也是无理数。1882年，德国数学家费迪南德·冯·林德曼证明了对任何非零代数数$\alpha$，$e^{\alpha}$都是超越数，该结论后来由魏尔斯特拉斯推广为林德曼－魏尔斯特拉斯定理。据此定理和欧拉公式，π只能是超越数，进而证实了勒让德和欧拉提出的π超越性猜想。哈代在其著作《数论导引》中则称此证明在提出后，经过希尔伯特、施瓦兹和其他一些人化简过。

### 引入符号

<image>

在用专指“圆周率”之前，希腊字母即已用于几何概念中。威廉·奥特雷德在1647年起在《数学之钥》（Clavis Mathematicae）就已经用及$\delta$（对应p和d的希腊字母）来表示圆的周长及直径的比例。

威廉·琼斯在他1706年出版的《新数学导论》（A New Introduction to the Mathematics）提到了，是目前已知最早专门用希腊字母表示圆周和其直径比例的人。这希腊字母第一次出现是在书中讨论一块半径1的圆时提到“其圆周长一半”。琼斯选用可能因它是希腊文“周边”一词“περιφέρεια”的首字母。不过琼斯提到，他那些有关的算式出自“真正聪明的约翰·梅钦先生”，人们推测在琼斯之前，约翰·梅钦就已开始用表示圆周率。

琼斯在1706年开始使用此希腊字母，但直到莱昂哈德·欧拉在其1736年出版的《力学 (书籍)》中开始使用之后，其他数学家才纷纷开始用指代圆周率。在此之前，数字家可能用像c或p之类的字母代表圆周率。欧拉与欧洲其他数学家间时常互相写信来往，的用法迅速传播开来。1748年欧拉在他的《无穷小分析引论》再一次提到了，写道：“简洁起见，我们将此数字写为，等于半径为1的圆周长的一半。”这表示方式之后也推展到整片西方世界。

## 现代数值近似

### 计算机时代与迭代算法

<image>

高斯－勒让德算法：一开始设定

$$\scriptstyle a_0 = 1 \quad b_0 = \frac{1}{\sqrt 2} \quad t_0 = \frac{1}{4} \quad p_0 = 1$$

迭代计算:$\scriptstyle a_{n+1} = \frac{a_n+b_n}{2} \quad \quad b_{n+1} = \sqrt{a_n b_n}$

$$\scriptstyle t_{n+1} = t_n - p_n (a_n-a_{n+1})^2 \quad \quad p_{n+1} = 2 p_n$$

则Π的估计值为

$$\scriptstyle \pi \approx \frac{(a_n + b_n)^2}{4 t_n}$$

二十世纪中期计算机技术发展、革新再次引发了计算Π位数的热潮。美国数学家约翰·伦奇及李维·史密斯在1949年用桌上型计算机计算到1120位。同年，乔治·韦斯纳（George Reitwiesner）及约翰·冯·诺伊曼带领的团队利用反三角函数（arctan）的无穷级数，用ENIAC计算到了小数后2037位，花了70小时的电脑工作时间。这纪录后来多次由其他透过arctan级数计算出的结果打破（1957年到7480位小数，1958年到第一万位数，1961年到第十万位小数），直到1973年，小数点后第一百万位小数经已算出。

1980年代有两项发明加速计算了π。第一项是发现了新的迭代法去计Π的值，计算速度比无穷级数快很多；另一项是发现了可以快速计算大数字乘积的乘法算法。电脑大部分的工作时间都是在计乘法，这类算法对现代计Π格外重要。这类算法包括嘉良对马（Karatsuba）算法、谭曲（Toom－Cook）乘法及以傅里叶变换为基础的乘法算法（傅里叶乘法）。

迭代算法最早是在1975年至1976年间分别由美国物理学家尤金·萨拉明及奥地利科学家理查·布兰特独立提出。这两条算法没有依赖无穷级数来计算。迭代会重复特定计算，将前一次的计算结果作为这一次的输入值，使得计算结果渐渐的趋近理想值。此方式的原始版本其实是在160年前由卡尔·弗里德里希·高斯提出，现在称为算术－几何平均数算法（AGM法）或高斯－勒让德算法。萨拉明及布兰特都曾修改之，这算法也称为萨拉明－布兰特算法。

迭代算法收敛速度比无穷级数快很多，在1980年代以后广为使用。无穷级数随着项次的增加，一般来说正确的位数也会增加几位，但迭代算法每计算多一次，正确位数会呈几何级数增长。例如萨拉明－布兰特算法每计算多一次，正确位数会是之前的二倍。1984年加拿大人乔纳森·波温及彼得·波温提出迭代算法，每计算多一次，正确位数会是之前的四倍，1987年时有另一条迭代算法，每计算多一次，正确位数会是之前的五倍。日本数学家金田康正使用的算法在1955年及2002年间创下了若干项纪录。不过迭代算法的快速收敛也有其代价，需要的内存明显比无穷级数多。

### 计算的意义

<image>

一般而言，值并不需要过于精确便能够满足大部分数学运算的需求。按照约·安（Jörg Arndt）及古里斯佗夫·希奴（Christoph Haenel）的计算，39位精确度已可将可观测宇宙圆周的精确度准确至一粒原子大小，足以运算绝大多数宇宙学的计算需求。尽管如此，和π有关的成就往往成为世界各地的新闻头条；部分人出于对破纪录的冲动，依然奋力算出π小数点后上千甚至上百万位。此外也有测试超级计算机、测试数值分析算法（包括）等实际好处。纯粹数学这领域也能计算π的位数评定其随机度。

### 快速收敛级数

<image>

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

$$\frac{1}{\pi} = \frac{2 \sqrt 2}{9801} \sum_{k=0}^\infty \frac{(4k)!(1103+26390k)}{k!^4(396^{4k})}.$$

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

$$\frac{1}{\pi} = \frac{12}{640320^{3/2}} \sum_{k=0}^\infty \frac{(6k)! (13591409 + 545140134k)}{(3k)!(k!)^3 (-640320)^{3k}}.$$

此公式每计算一项就能得到π的约14位数值，因而用于突破圆周率的数位的计算。利用这公式，楚德诺夫斯基兄弟于1989年算得π小数点后10亿（109）位，法布里斯·贝拉于2009年算得2.7千亿（2.7×1012）位，亚历山大·易和近藤滋在2011年算得一万亿（1013）位。类似的公式还有拉马努金-佐藤级数。

2006年，加拿大数学家西蒙·普劳夫利用PSLQ整数关系算法按照以下模版生成了几条计算π的新公式：

$\pi^k = \sum_{n=1}^\infty \frac{1}{n^k} \left(\frac{a}{q^n-1} + \frac{b}{q^{2n}-1} + \frac{c}{q^{4n}-1}\right)$。

$q$为e，$k$是奇数，$a,b,c$是普劳夫计算出的有理常数。

### 统计模拟法

布丰投针问题，多枚长度为ℓ的针随机地抛掷向平面。

统计模拟法是以概率统计理论为指导的一类非常重要的计数方法，经大量重复试验计算事件发生频率，按照大数定律（即当试验次数充分大时，频率充分接近概率）可以求得的近似值。布芬（Buffon）投针问题就是其中一项实例：长度$l$的针随机往画满间距$t\left ( l\leq t \right )$的平行线的平面上抛掷$n$次，如果针与平行直线相交$m$次，$n$充分大就可根据以下公式算出的近似值：

$$\pi \approx \frac{2n\ell}{mt}$$

用统计模拟法计的另一例子是随机往内切四分之一圆的正方形内抛掷大量点，落在四分之一圆内的点的数量与抛掷点的总量的比值会近似于$\frac{\pi}{4}$。

此外还可用随机游走试验，并用统计模拟法计算值，如抛掷一枚均匀的硬币$N$次，并记录正面朝上的次数，所得结果中，正面朝上的次数$n_N$服从二项分布且

$$\Pr(n_N = m) = \binom N m (\frac{1}{2})^m (\frac{1}{2})^{N-m}$$

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

$$W_N = \sum_{k=1}^N X_k$$

设k为“硬币正面朝上的次数”减去“硬币反面朝上的次数”，即可得到$m-\left ( N-m \right )=k$。变换式子，得$m=\frac{N+k}{2}$，因此

$\Pr(W_N = k) = \binom N \frac{N+k}{2} \frac{1}{2^N}$，其中$k=-N,-N+2,-N+4,\ldots,N-2,N$。

可证明。

$E(W_N) = 0$，$E(W_N^2) = N$，以及$E(|W_N|) = \binom{N}{\left\lceil{N/2}\right\rceil \frac{\left\lceil{N/2}\right\rceil}{2^{N-1 = \begin{cases} \frac{(N-1)

},& \text{若 }N\text{偶，} \\ \frac{N

},& \mbox{若 }N\mbox{奇。} \end{cases}$

并且当N变大时，$E\left ( \left \vert W_N \right \vert \right )$的值会渐近于$\sqrt{\frac{2N}{\pi}}$，因此当N充分大时可根据以下公式算出的近似值：

$$\pi \approx \frac{2N}{|W_N|^2}$$

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

### 阀门算法

1995年引入的两条算法开辟了研究的新途径。因为每计算出一位数字，该数就会像流过阀门的水一样不会再出现在后续的计算过程中，这种新进算法叫阀门算法。这就与无穷级数及迭代算法形成对比——无穷级数和迭代算法自始至终的每一步计算都会涉及到之前所有步骤计算出的中间值。

1995年，美国数学家斯坦·瓦格纳和斯坦利·拉比诺维茨（Stanley Rabinowitz）发明了一种简单的阀门算法，其运算速度类似arctan算法，但速度比迭代算法慢。

贝利－波尔温－普劳夫公式（BBP）是另一条阀门算法，属于一种位数萃取算法。1995年，西蒙·普劳夫等人发现

$$\pi = \sum_{k=0}^\infty \frac{1}{16^k} \left( \frac{4}{8k + 1} - \frac{2}{8k + 4} - \frac{1}{8k + 5} - \frac{1}{8k + 6}\right)$$

这公式和其他公式不同，可以计算的任何十六进小数位，而不用计算前面全部小数位。十六进数位可计算得到特定二进数位；想要得到八进制数位的话，计算一、两位十六进小数即可。目前也已发现一些这种算法的变体，不过还没有发现针对十进制、可以快速生成特定小数位的位数萃取算法。位数萃取算法的一项重要用途是用来确认声称是计算到小数位数的新纪录：若有声称是新纪录的计算结果出现，先将十进制的数值转换到十六进制，再用贝利－波尔温－普劳夫公式去确认最后一些位数（用乱数决定），若这些位数都对，就能有一定把握认为此计算结果是对的。

1998年到2000年间，分布式计算计划PiHex用贝拉公式（贝利－波尔温－普劳夫公式的一种变体）计算第1015位，结果是0。2010年9月，有雅虎员工用公司的Apache Hadoop应用程序在上千台电脑计算π在2×1015位开始往后256位，其第2×1015位刚好也是0。

### 利用伽玛函数计算

伽玛函数，$\Gamma(n) =(n-1)!$，可以被用作计算圆周率。

$$\pi= \Gamma(\frac{1}{2})^2$$

证明如下：

利用欧拉反射公式

$$\Gamma(z) * \Gamma(1-z) = \frac{\pi}{\sin( \pi z)}$$

令 $z = \frac{1}{2}$

$$\Gamma(\frac{1}{2})*\Gamma(1-\frac{1}{2}) = \frac{\pi}{\sin (\frac{\pi}{2})}$$

因为 $\sin(\frac{\pi}{2})$

$$\Gamma(\frac{1}{2})^2 = \pi$$

## 用途

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

### 几何学与三角学

<image>

出现在基于圆的几何图形（如椭圆、球、圆锥与环面）的面积、体积公式中。下面是一些用到π的常见公式：

- 半径$r$的圆周长$2\pi r$。
- 半径$r$的圆面积$\pi r^2$。
- 半径$r$的球体积$\frac{4}{3}\pi r^3$。
- 半径$r$的球面面积$4\pi r^2$。

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

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

$$\int_{-1}^1 \sqrt{1-x^2}\,dx = \frac{\pi}{2}.$$

$\sqrt{1-x^2}$的积分表示上半圆（此处的平方根由勾股定理得出），从积分上下限-1到1的积分$\int ^x_{-1}$可用来计算半圆与x 轴间的面积。

<image>

三角函数要用到角，而数学家常用弧度作角度单位。π在弧度制起重要作用，数学家将周角，即360度定义为$2\pi$度。由这条定义可得，180度＝弧度，1度＝$\frac{\pi}{180^\circ}$弧度。因此，常用的三角函数的周期为的倍数；例如，正弦和余弦周期为π，任何角度$\theta$和任何整数$k$都有

$\sin\theta = \sin\left(\theta + 2\pi k \right)$及$\cos\theta = \cos\left(\theta + 2\pi k \right)$。

### 拓扑学

<image>

常数出现在将平面微分几何及其拓扑学联系起来的高斯－博内定理中。具体来说，如果紧曲面Σ的高斯曲率为$K$，那么有

$\int_\Sigma K\,dA = 2\pi \chi(\Sigma)$。

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

$$A(S) = \int_S 1\,dA = 2\pi\cdot 2 = 4\pi$$

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

常数还出现在拓扑学的许多其他的积分公式中，特别是那些涉及通过陈－韦伊同态的特征类。

### 向量分析

<image>

向量分析是与向量场的性质有关的微积分的分支，并有许多物理用途，例如用在电磁学中。位于三维笛卡尔坐标系原点的点源$Q$的牛顿位势为

$$V(\mathbf{x}) = -\frac{k Q}{|\mathbf{x}|}$$

表示位于距原点$\left \vert \boldsymbol{x} \right \vert$的单位质量（或电荷）的势能，而$k$是维度常数。在这里由$\mathrm{E}$表示的场可以是（牛顿）引力场或（库仑）电场，是位势的负梯度：

$$\mathbf{E} = -\nabla V.$$

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

$4\pi k Q =$ ${\scriptstyle S}$ $\mathbf{E} \cdot d\mathbf{A}$

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

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

$$\Delta V(\mathbf x) = -4\pi k Q\delta(\mathbf x).$$

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

$$\Delta V(\mathbf x) = -4\pi k \rho(\mathbf x)$$

其中$\rho$是分布函数。

<image>

常数在与爱因斯坦场方程中的四维势起类似的作用，爱因斯坦方程是形成广义相对论基础的一条基本公式，并且把引力的基本相互作用描述为物质和能量引起的时空弯曲的结果：

$$R_{\mu\nu} - \frac{1}{2} R g_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8 \pi G}{c^4} T_{\mu\nu},$$

$R_{\mu v}$是里奇曲率张量，$R$是数量曲率，$g_{\mu v}$是度量张量，$\Lambda$是宇宙学常数，$G$是万有引力常数，$c$是真空中的光速，而$T_{\mu v}$是应力－能量张量。爱因斯坦方程的左边是度量张量的拉普拉斯算子的非线性模拟，并化简（reduce）至在弱域的极限，而右边是分布函数的模拟乘以$8\pi$。

### 柯西积分公式

<image>

在复分析中，沿复平面若尔当曲线的围道积分是研究解析函数的重要手段之一。简化版的柯西积分公式表明，对任何若尔当曲线$\gamma$内任一点$z_0$，以下围道积分给出$2\pi i$：

$$\oint_\gamma \frac{dz}{z-z_0} = 2\pi i.$$

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

一般形式的柯西积分公式建立了全纯函数$f(z)$在若尔当曲线$\gamma$上的值与曲线内任意点$z_0$处值的关系：

$$\oint_\gamma { f(z) \over z-z_0 }\,dz = 2\pi i f (z_{0})$$

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

$$\oint_\gamma g(z)\, dz =2\pi i \sum \operatorname{Res}( g, a_k )$$

### Γ函数与斯特灵公式

<image>

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

$$\Gamma(z) = \int_{0}^{\infty} \frac{t^{z-1}}{\mathrm{e}^t} \,{\rm{d}}t$$

再利用解析延拓可以将它的定义域扩展到除去非正整数的整块复数域。当自变量$z=n$取正整数时，$\Gamma$函数给出阶乘$\left ( n-1 \right )!$；当自变量取半整数时，计算结果含有。例如$\Gamma(\frac{1}{2}) = \sqrt{\pi}$，$\Gamma(\frac{5}{2}) = \frac {3 \sqrt{\pi}} {4}$。

根据魏尔施特拉斯分解定理，$\Gamma$函数可分解为如下的无穷乘积：

$$\Gamma(z) = e^{-\gamma z}\prod_{n=1}^\infty \frac{e^{\frac{z}{n}}}{1+\frac{z}{n}}$$

$\gamma$是欧拉－马斯刻若尼常数。利用该分解公式和$\Gamma$函数在$z=\frac{1}{2}$的值$\Gamma\left ( \frac{1}{2} \right )^2=\pi$，亦可以证明沃利斯乘积式。$\Gamma$函数和黎曼ζ函数、的恒等式存在关联，其中扮演着重要的角色。

$\Gamma$函数常用于计算$n$维欧氏空间中n 维球的体积和n 维球面的表面积。对$n$维欧氏空间中半径为$r$ $n$维球，其体积$V_n(r)$和表面积$S_{n-1}(r)$满足：

$$V_n(r) = \frac{\pi^{\frac{n}{2}}}{\Gamma(\frac{n}{2}+1)}r^n$$

$$S_{n-1}(r) = \frac{n\pi^{\frac{n}{2}}}{\Gamma(\frac{n}{2}+1)}r^{n-1}$$

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

$$2\pi r = \frac{S_{n+1}(r)}{V_n(r)}.$$

当$n$很大，用$\Gamma$函数可得到阶乘$n!$的近似公式$n! \sim \sqrt{2 \pi n} \left(\frac{n}{e}\right)^n$，称斯特灵公式，等价于：

$$\pi = \lim_{n\to\infty} \frac{e^{2n}n!^2}{2 n^{2n+1}}.$$

斯特灵近似的几何应用之一是埃尔哈特体积猜想。将$n$维欧几里得空间的单纯形记作$\Delta_n$，$\left ( n+1 \right )\Delta_n$则表示该单纯形的所有面扩大$n+1$。于是

$$\operatorname{Vol}((n+1)\Delta_n) = \frac{(n+1)^n}{n!} \sim \frac{e^{n+1}}{\sqrt{2\pi n}}.$$

这是仅含一点晶格点之凸体体积的（最佳）上界。

### 数论与黎曼ζ函数

<image>

<image>

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

$$\zeta(2) = \frac{1}{1^2} + \frac{1}{2^2} + \frac{1}{3^2} + \cdots$$

找到这无穷级数的解析解是数学界著名的“巴塞尔问题”。1735年，欧拉解决了这问题，他得到该无穷级数等于$\frac{\pi^2}{6}$。欧拉的结论可推导出数论中一项结果，即两随机整数互质（无公因数）的概率为$\frac{6}{\pi^2}$。整数可由质数$p$整除的概率为$\frac{1}{p^2}$（例如，连续7个正整数只有一个可以7整除），任取两随机整数都能以质数$p$整除的概率为$\frac{1}{p^2}$，至少有一数不能整除的概率则为$1-\frac{1}{p^2}$。又，一随机整数能否以两不同质数整除是相互独立事件，两随机整数互质的概率可以表示成关于所有质数$p$的无穷乘积

$$\begin{align} \prod_p^\infty \left(1-\frac{1}{p^2}\right) &= \left( \prod_p^\infty \frac{1}{1-p^{-2}} \right)^{-1}\\ &= \frac{1}{1 + \frac{1}{2^2} + \frac{1}{3^2} + \cdots }\\ &= \frac{1}{\zeta(2)} = \frac{6}{\pi^2} \approx 61\%. \end{align}$$

这结论可结合随机数生成器，用统计模拟法计的近似值。

巴塞尔问题的结论意味着几何导出量的数值与质数的分布有着深刻的关联。巴塞尔问题是谷山－志村定理的一种特殊情况，是安德烈·韦伊对玉河数的猜想的一项特例，即猜想一个这种形式的算术量关于所有质数$p$的无穷乘积能够等于一个几何量——某局部对称空间体积的倒易。巴塞尔问题中，这空间是双曲3-流形SL2(R)/模群。

$\zeta$函数同样满足黎曼方程的公式，其中用到了和伽玛公式：

$$\zeta(s) = 2^s\pi^{s-1}\ \sin\left(\frac{\pi s}{2}\right)\ \Gamma(1-s)\ \zeta(1-s)\!.$$

除此之外。

$\zeta$函数导数也满足

$$\exp(-\zeta'(0)) = \sqrt{2\pi}.$$

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

### 傅里叶级数

<image>

周期函数的傅里叶级数很自然出现了。周期函数即实数的小数部分所构成群$T=\frac{R}{Z}$上的函数。傅里叶分解指出，$T$上的复值函数$f$可表示为无穷多个$T$的酉特征的线性叠加之和。也就是说，$T$到圆群$U(1)$（模为1的复数组成的乘法群）的映射是连续群同态。$T$的特征都有$e_n(x)= e^{2\pi i n x}$的形式，是一条定理。

$T$有唯一的特征值，直到复共轭，那是一群同态。在圆群用哈尔测度，常数是这特征值的拉东－尼科迪姆导数值的一半。其他的特征值的导数值为$2\pi$的正整数倍。因此，常数是独特的数字，以至于配备了其哈尔测度的群$T$，有对于$2\pi$整数倍的点阵的庞特里亚金对偶性。这是泊松和公式的一维版本。

### 模形式与函数

常数与模形式和Θ函数密切相关——比如，椭圆曲线中的j变量就很大程度涉及楚德诺夫斯基算法（一种快速计算Π的方法）。

模形式是以在上半平面的全纯函数的在模群$SL_2(\mathbb Z)$（或其子群，$SL_2(\mathbb Z)$是$SL_2(\mathbb R)$的一格）下的变换特性归纳。Θ函数便是一例：

$$\theta(z,\tau) = \sum_{n=-\infty}^\infty e^{2\pi i nz + i\pi n^2\tau}$$

它是一种名为雅可比形式的模形式，有时以诺姆 (数学)$q=e^{\pi i \tau}$表达。

常数是特殊常数，它会使雅可比$\Theta$函数形成自守形式，即该函数会以特定方式变换。有若干恒等式在所有自守形式下成立。例如：

$$\theta(z+\tau,\tau) = e^{-\pi i\tau -2\pi i z}\theta(z,\tau)$$

它使得$\theta$必然在离散海森伯群下以表示（representation）变换。一般模形式和其他$\Theta$函数也包含，这也是根据史东–冯纽曼定理。

### 柯西分布与位势论

<image>

积分$\int_{-\infty }^{\infty } \frac{1}{x^2+1} \, dx = \pi$，柯西分布$g(x)=\frac{1}{\pi}\cdot\frac{1}{x^2+1}$是概率密度函数，总概率等于1。

柯西分布的香农熵等于$\log(4\pi)$，也含。

<image>

柯西分布在位势论中扮演着重要的角色因为它是最简单的福斯坦堡测度和与在半平面上做布朗运动相关联的经典泊松核。共轭谐波函数以及希尔伯特变换与泊松核的渐近线有关。希尔伯特变换$H$是由奇异积分的柯西主值给出的积分变换

$Hf(t) = \frac{1}{\pi}\int_{-\infty}^\infty \frac{f(x)}{x-t}$。

常数是唯一的（正）归一化因子因此$H$定义了一个在实数轴上的平方可积分实值函数的希尔伯特空间上的线性复结构。和傅里叶变换一样，希尔伯特变换就其在希尔伯特空间$L^2(R)$的变换特性而言可以完全特征化。直到归一化，它是唯一的与正膨胀对易且与实数轴的所有反射反对易有界线性算子。常数是唯一能使这变换幺正的归一化因子。

### 复变动态系统

<image>

大卫·波（David Boll）在1991年发现在曼德博集合分形也有Π出现。他检查在曼德博集合在$(-0.75,0)$位置的特性。若考虑坐标在“颈部”$(-0.75,\varepsilon)$的点，而$\varepsilon$趋近零，在发散之前迭代的次数和$\varepsilon$相乘，会趋近。若是在右侧尖点处附近的点$(0.25,\varepsilon)$也会有类似的特性：在发散之前迭代的次数和$\varepsilon$的平方根相乘，也会趋近。

## 数学之外的Π

### 描述物理现象

与圆以及球坐标系关系密切，即使不是物理常数，也常出现在描述宇宙的基本原则方程中。比方说，经典力学领域的简单公式给出长L的单摆小幅摆动的近似周期$T \approx 2\pi \sqrt\frac{L}{g}$，$g$为地球引力加速度常数。

海森堡不确定性原理是量子力学的基本公式，表明测量粒子时，其位置不确定度（$\Delta x$）与动量不确定度（$\Delta p$）不可能同时达到任意小（$h$为普朗克常数）：$\Delta x\, \Delta p \ge \frac{h}{4\pi}$。

近似三这特性，和电子偶素的半衰期相对较长有密切的联系。其半衰期的倒数和精细结构常数$\alpha$的关系为$\frac{1}{\tau} = 2\frac{\pi^2 - 9}{9\pi}m\alpha^{6}$，$m$为电子质量。

许多结构工程的公式也有，例如欧拉推导的挫曲公式说明了长度为$L$、截面二次轴矩为I的细长形物体，在不挫曲的条件下可以承受的最大轴向负载$F$：

$$F =\frac{\pi^2EI}{L^2}.$$

流体动力学的斯托克斯定律中也有。斯托克斯定律是半径约为$R$的小球体在黏度$\eta$的流体中以速度$v$运动时会受到的阻力满足：

$$F =6 \, \pi \, \eta \, R \, v .$$

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

### Π的记忆技巧

π文字学（或译作圆周率的语言学）是指记住的大量位值，并将其世界纪录载于吉尼斯世界纪录大全的做法。维尔·美拿（Rajveer Meena）于2015年3月21日在印度于9小时27分钟内背诵了7万位的Π，创下吉尼斯世界纪录大全认证的世界纪录。2006年，日本退休工程师原口证在千叶县于官员见证下背诵了十万位小数，但他未获吉尼斯世界纪录大全认证。

常用于记忆Π的一项技巧是背诵以单词长度代表数值的故事或诗歌：第一单词有三字母，第二单词有一字母，第三单词有四字母，第四单词有一字母，第五单词有五字母，如此类推。早期例子是英国科学家詹姆士·金斯设计的诗歌：“How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics.”这类诗歌有时在英文中称为“piem”。除了英文，用于记忆Π的诗歌亦有不同语言的版本。但是，记忆的人一般并不以诗歌记忆来创下纪录，而是用如记忆数字规律或轨迹法的方法。

有好几位作家仿照上述记忆技巧，用的数值创作了新型的约束写作方式，当中单词长度须符合的数值。《The Cadaeic Cadenza》以上述技巧包含了前3835位的值，一本标准长度的书《Not a Wake》有一万单词，其中各单词亦代表了的一位。

### 大众文化

<image>

也许因为的公式很简短而且四处可见，比其他数学常数在流行文化中更常见。

在2008年由英国公开大学及英国广播公司联合制作的记录片《数学的故事》于2008年十月由英国广播公司第四台播放。影片讲述了英国数学家马库斯·杜·索托伊在到访印度研究当地三角学的贡献时，展示出历史上Π最精确的计算公式的信息图形。

巴黎的科学博物馆发现宫有间圆形房间叫“房”，墙上刻有的707位数，数字贴在圆顶状的天花板上，由大型的木制字符组成。数值是1853年由英国数学家威廉·尚克思计算出来，但是该结果于第528位后开始出现谬误，在1946年发现，1949年修正。

卡尔·萨根的小说《接触未来》中则暗示说，宇宙的创造者在Π的数字中暗藏了一则信息。Π的数字也用在凯特·布希所出的专辑中的《Pi》的歌词里。

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

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

```
在1958年，阿尔伯特·伊格尔[[轉 (角)|提议]]将换成[[Τ|<span class="texhtml">&tau;</span>（tau）]]以便简化公式。<math>\tau</math>在此定义为的兩倍<ref name="AE1958"></ref>。然而，没有任何其他作者曾这样使用过<math>\tau</math>。有些人使用不同的值，<math>\tau=2\pi=6.283185\ldots</math>。<ref>Sequence A019692</ref>这些人称<math>\tau</math>不论是作为弧度制下圆周长的1[[轉 (角)|转]]还是作为弧长与半径的比值（而不是与直径的比值）都比自然，也能因此简化许多公式。<ref></ref><ref></ref>有媒体报道称，因为<math>\tau</math>的值大小约为6.28，現已有人在6月28日庆祝“<math>\tau</math>节”，并吃“两个派”；<ref></ref>然而，主流数学界还并未使用<math>\tau</math>。<ref></ref>
```

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

## 参考资料

### 书籍

- English translation by Catriona and David Lischka.
- , English translation by Stephen Wilson.
- . issue 3 Jan/Feb, issue 4 Mar/Apr, issue 5 May/Jun

## 延伸阅读

- 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
- 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
- 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
- 威廉·尚克斯, Contributions to Mathematics Comprising Chiefly of the Rectification of the Circle to 607 Places of Decimals, 1853, pp. i–xvi, 10
- 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
- 史丹·维根, "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

## 外部链接

- Science/Math/Recreations/Specific_Numbers/Pi/Digits/
- Wolfram Mathworld上的 “Pi”
- Wolfram Alpha上的 Representations of Pi
- Pi搜索引擎：可搜寻Π、2和e的20亿位
- Demonstration by Lambert (1761) of irrationality of Π, online and analyzed BibNum (PDF).
- 寻找π值的计划
- 前1.2百万位中的部分资料
- 将π前一万位化作音乐旋律
- SuperPI 计算π值的软件，电脑硬件玩家常用来测试电脑运算速度（日文）
- 计算圆周率
- PiFast 个人电脑上最快的计算π值软件，是个人电脑计算π值纪录保持软件。
- 用大炮求圆周率：以蒙特卡罗算法，利用圆的面积求圆周率。
- 展示将π值计算到10万位的网站

## 基本资料

- 作战范围：数表—无理数$\color{blue}\sqrt{2}$ - $\color{blue}\varphi$ - $\color{blue}\sqrt{3}$ - $\color{blue}\sqrt{5}$ - $\color{blue}\delta_S$ - $\color{blue}e$ - $\color{blue}\pi$
- nav：no
- 货币符号：$\pi$
- 面值：$\pi\approx$
- OEIS编号：A000796
- 类型：无理数、超越数
- 定义：$\pi = \frac{C}{d}.$，其中$C$为圆周长、$d$为直径、$\pi = \int_{-1}^1 \frac{dx}{\sqrt{1-x^2}}.$
- 以此为根的多项式或函数：$e^{ix}+1 = 0$
- 连分数：$\pi=3+\textstyle \frac{1}{7+\textstyle \frac{1}{15+\textstyle \frac{1}{1+\textstyle \frac{1}{292+\textstyle \frac{1}{1+\textstyle \frac{1}{1+\textstyle \frac{1}{1+\ddots}}}}}}}$
- 系列：$\pi = \sum_{k=0}^\infty\frac{4(-1)^k}{2k+1}$
- 二进制：2
- 十六进制：16
- 六十进制：60…