圆周率

圆的周长与直径的比值

圆周率 (英文 名称为: Pi ),常用符号π来表示。圆周率的定义为:圆形的 周长 与直径之比,公式为π= C / D , D 为圆形的直径; C 为圆形的周长 。π既是一个 无理数 也是一个 超越数 ,并且他可以表示为无穷连分式的形式 。

中国 数学家 刘徽 是中国最早运用科学方法计算π值的人,约263年,他用割圆术得出π=3.14。此外,中国数学家 祖冲之 也求得 圆周 率3.1415926<π<3.1415927,他算出π的8位可靠数字,保持世界记录九百多年 。至1949年, 英国 数学家列维·史密斯和雷恩奇算出π的1121位,是人工算π的最高记录 。 1949年,人类第一次使用 电子 计算机 计算π ,至2024年,π值已计算到 小数点 后约105万亿位。

π是数学上的一个重要 常数 ,与许多数学问题、数学公式有关 ,例如初等数学中的 圆周长 、圆面积、椭圆面积、球体积计算等 ,在数学的无穷 、 级数 、积分 的计算中也常常使用π ,很多函数中也可以看到π的存在,如: 柯西分布 、伽玛( gaussian )函数等 。在数学之外的 介子 理论、交流电路等方面同样引入π来进行计算推导 。

定义

几何定义

圆周率定义为圆形的 周长 与直径之比,它也等于圆形之面积与半径平方之比 。

即:

为圆形的直径; 为圆形的周长 。

积分定义

19世纪, 奥古斯丁-路易·柯西 、 卡尔·魏尔施特拉斯 给出了极限、收敛的精确定义,确立了以极限理论为基础的 数学分析 体系,使微积分有了严格的理论基础 。 英国 数学家 约翰·沃利斯 对单位圆的研究,得到 ,这个式子表示的是一个单位圆的面积,后来, 瑞士 数学家 莱昂哈德·欧拉 给出了通式: 。

复数定义

复数 ,可以理解为一个 实数 和一个 虚数 复合而成的“复合数” ,通常称 的表达式为复数,其中 和 为任意实数, 为 虚数单位 。

由欧拉 恒等式 可知, 是定义为一个单位长度在圆弧上旋转半圈后得到的数,该数就是 ,所以 。其中 为 自然对数 的 底数 。

符号

圆周率符号为 (英文名称为: ),取自 希腊语 “周围—πúρωαπó”的第一个字母 。1706年,威尔士数学家威廉·琼斯(Wiliam Jones,1675—1749年)引入了符号 代表圆周率 。

简史

实验获取阶段

公元前950年前后, 基督教 《圣经》中最早记载了圆周率为“3”, 巴比伦 、 印度 、 中国 等也长期使用这个粗略而简单实用的数值,约公元前250年, 古巴比伦 的一块石碑上写着,圆与它的内接正六边形的 周长 之比为1:0.96,这意味着圆周率 的值为3.125 。中国的 刘徽 提出的“圆径一而周三”曾广泛流传,《 周髀算经 》中,就记载有这一结论,木工师傅有两句从古流传下来的口诀“周三径一,方五斜七”,意思是说,直径为1的圆,周长大约是3,边长为5的正方形, 对角线 之长约为7。这正反映了早期人们对圆周率 和 这两个 无理数 的粗略估计。 东汉 时期还 明文 规定圆周率取3为计算面积的标准,后人称之为“古率” 。

几何算法阶段

圆内接、外切多边形

阿基米德 ( Archimedes ,公元前287—公元前212年)是科学地研究圆周率这一 常数 的第一个人,开创了几何计算 的阶段,真正使圆周率计算建立在科学的基础上,他提出了一种能够借助数学过程,把 的值精确到任意精度的方法。由图1知, 圆周长 大于内接正四边形而小于外切正四边形,因此 < < ,在他的一篇论文《圆的测定》之中, 阿基米德 第一次使用上、下界来确定 的近似值,他用几何方法证明了“圆周长与圆直径之比小于 而大于 ”到150年左右,希腊天文学家托勒密(Ptolemy,90—168年)得出 3.1416,这是自 阿基米德 以来取得的巨大进步 。

割圆术的应用

割圆术是 中国 数学家 刘徽 提出的,他是中国最早运用科学方法计算 值的人,在263年前后,他就用此方法得出 =3.14,后人称之为“徽率”。割圆术仅用内接 正多边形 就确定出了圆周率的上、下界,比 阿基米德同时用内接和外切正多边形的方法简捷得多,割圆术中刘徽提供了一种精加工方法,通过计算圆内接正3072边形的周长,得出圆周率π的近似值约为3.1416。此外,中国 的另一位数学家 祖冲之 (429—500年 )也求得圆周率3.1415926< <3.1415927,同时得到 的两个近似分数,即 约率 为 , 密率 为 ,他算出的 的8位可靠数字,不但在当时是最精密的圆周率,而且保持世界记录九百多年,被命名为“ 祖率 ” 。

1150年, 印度 数学家婆什迦罗第二(1114—1185年) 则算出 的近似值为 3.1416 。约1424年, 阿拉伯 数学家 卡西 经过演算805306368个内接与外切 正多边形 的 周长 ,最终得到 3.14159265358979325,这个 值有17个 有效数字 ,首次突破由 祖冲之 所创造的记录 。

分析算法阶段

解析表达式的发现

1579年, 法国 数学家韦达(Francois Viete,1540—1630年)在《数学定律,应用于三角形》中,通过计算圆的正 393216边形,得出3.1415926535< <3.1415926537,同时他利用分析式和 级数 乘积来刻画 :

cos cos cos

=

这个式子给科学家们指出了一个崭新的计算 的思路,这是分析法计算圆周率时代的第一个解析表达式。1621年, 德国 — 荷兰 物理学家、数学家斯涅耳(Willebrod Snell,1580—1626年)对原来古典方法作出了一些三角上的改进,利用这一改进,只需算到 边形,就可得到34位 值。1630年, 罗马 数学家格林贝格(Peter Andreas Grünberg,1939—不详)利用斯涅耳的方法得到40位 值 。

反正切函数表达式

在1706年 英国 数学家、天文学家梅钦发现了Machin公式:

actan arctan

arctan + - + ,利用此公式计算得到了π的100位数值,具体为:3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679。

1734—1735年,在前人的基础上, 莱昂哈德·欧拉 (Leonhard Euler,1707—1783年)得出著名的欧拉公式 :

+ + + +

到了1737年,欧拉利用反正切函数arctan ,根据格雷戈里的展开式将右边展开,得到了著名的欧拉级数 :

1844年,德国汉堡的数学家斯特拉斯尼茨基和达什使用施瓦兹公式计算到小数点后200位。1949年, 英国 数学家列维·史密斯和雷恩奇算出 的1121位,是人工算 的最高记录 。

概率方法的使用

18世纪 法国 科学家 蒲丰 (Buffon,1707—1788年)创造的一种几何与分析思想之外的方法求 的值,即“投针求的 概率论 方法” 。

将长为 的匀质细针随机地掷于画了等距 平行线 族的平面上,相邻两平行线距离为 ,求针与平行线 相交 的概率。

解法:设 为针 中点 到平行线的距离,则 ,设 为针与平行线的夹角,则 。

由图(1)知,当且仅当 sin , 时,时针与 平行线 相交,于是所求概率:

,于是 。1901年, 意大利 数学家拉泽里尼投针3408次,求得 3.141592 。

电子计算机的使用

1949年,人类第一次在 美国 马里兰州 阿伯丁的军队 弹道 研究实验室里用电子 计算机 算 ,包括 约翰·冯·诺依曼 在内,用ENIAC和梅钦公式计算 值到2048位,突破了千位数。1959年, 美国国际 商业机器公司(IBM)制成第二代电子计算机—世界上的一台 晶体管 电子计算IBM—7090 。1961年,丹尼尔·尚克斯和雷恩奇用IBM—7090分别使用 挪威 斯托默 公式:

arctan arctan arctan

和 高斯 公式:

arctan arctan arctan

将 计算到 小数点后100256位,达到10万位高峰 。

1985年,乔纳森·波尔文和彼得·波尔文发表波尔文算法:

初值: ,

重复计算: 和 ,

最后得到 ,它四次收敛于圆周率 。

2002年, 日本 的 金田康正 利用

初值: , , ,

重复计算: , , , ,

最后计算: = 在 日立 制作所 的超并行超级电脑系统日立—SR—8000—MPP下,净耗时400多小时,算出12411亿位值,这一成绩刷新了2061亿位的 世界纪录 。

2020年,美国北 阿拉巴马 慈善计算组织的创始人蒂莫西·穆利肯使用个人电脑,耗时303天,将圆周率的值计算到 小数点 后约50万亿位。2021年, 瑞士 格劳宾登州 应用科学大学的研究人员宣布,他们借助研究所的一台 计算机 ,耗时108天,将 值计算到 小数点 后约62.8万亿位,创下截止到2020年 值计算的最精确纪录 。

几何定义的推倒过程

的推导过程如下:

任意取一个直径是 的圆,作这个圆的内接 正多边形 ,作内接正 边形,并且把正 边形的边长记做 , 周长 记做 。那么把这个正 边形的边数加倍,且无限制地继续加倍的过程,就得到一系列边数递增的圆内接正多边形,这些正多边形的周长组成一个无限 数列 :

(1)

设圆的周长为 ,因为三角形两边的长度的和大于第三边的长度,所以内接正 边形的边长 小于内接正 边形的边长 的两倍,就是 ,因此 …所以有 :

(2)

对于一切的 ,得到 :

(3)

根据(1)式和(2)式可以推导出 数列 (2)递增数列;另一方面,这个数列的每一项都小于一个常量,这时,数列(1)就有极限,当圆内接 正多边形 的边数无限增加的时候,内接正多边形的 周长 就接近于一个确定的值,这个值叫做圆的周长 :

(4)

这里 为 正多边形 的边数, 是它的周长,设圆的直径是 ,它的周长是 ,圆内接正 边形的边长、周长分别是 、 , 如图 2,得到 :

,因此 ,取极限后就得到:

这就说明了圆的周长和直径的比值是一个 常数 。把这个常数记作 就得到 :

性质

也是一个常数,古希腊 欧几里得 《 几何原本 》(约公元前3世纪初)中提到圆周率是常数, 中国 古算书《 周髀算经 》(约公元前2世纪)中有“径一而 周三 ”的记载,也认为圆周率是常数 。

无理性

是一个 无理数 ,也就是说它是一个无限不 循环小数 。1767年, 德国 物理学家约翰·海恩里希·兰伯特(Johann Heinrich Lambert,1728—1777年)率先证明了 是一个无理数,根本找不到任何分数表达式 。证明 是无理数时, 兰伯特 首先证明了将tan 写成连分式的展开形式,当 是一个非零的正整数时,tan 是一个无理数。 因为tan ,这就意味着 是一个无理数 。

超越性

1882年, 德国 数学家 费迪南德·冯·林德曼 证明了 是 超越数 ,既 不可能是任何整数 多项式 的根 。当一个数可以被写成含有理系数的多项式 方程 的根的形式时,不管这个数是 实数 还是 复数 ,则这个数都可以被定义为 代数数 ,否则,就是超越数。这就是说,如果存在非零的 有理数 使得方程 成立,我们就说式中的 是一个 代数数 。而当 为一个 超越数 时,这个数就不是任何一个含非零的 有理数 系数的 多项式 方程 的根。 费迪南德·冯·林德曼 就证明了圆周率 就是这样一个超越数 。

的 无理数 证明,通俗易懂的是 美国 数学家尼云(Ivan Morton Niven 1915—1999年)在1947年给出的如下证法:

设 ,其中 是 互素 的正整数,构造函数 以及

计算可知

由于

因此,

对 求直至 阶的 偏导数 ,注意到 和 都是整数,可推得 和 也都是整数,由于 时 sin ,因此

sin ,这说明 是一个正整数,但是 时有 sin ,根据 夹逼 定理 ,当 时,有

sin sin

这显然与 是一个正整数相矛盾,因此 是 无理数 。

连分式

圆周率 也可以用无限 级数 的嵌套分数的连分式形式来表示,表达式为 :

即:

计算方法举例

利用级数

幂 级数 在 时,收敛于 ,即

= ( )(1)

对此逐项积分得:

arctan ( )(2)

,因右端 级数 当 时收敛,故上式在 时也成立,当 时,由(2)得

(3)

此式又被称为Π的莱布尼茨公式,是史上首条π的精确 无穷级数 公式。 然而,此式右端的 交错级数 是 条件收敛 而非 绝对收敛 ,若用其前 项和作近似,产生的误差将介于 与 之间,因此该 级数 的 收敛速度 是不能令人满意的 。

为了在计算时收敛速度快一些,可取较小的 正数 为 ,并代入(2)式,一般来说, 的 绝对值 越小,该 幂级数 的收敛速度越快,而 的绝对值越大,该级数的收敛速度越慢 。如在(2)式的两端令 ,得等式:

= (4)

(4)式中 级数 的 收敛速度 显然要优于(3)式中级数的收敛速度,但是 级数 中含有的根式不易计算。17世纪后的许多著名数学家想到,将 分成两个或两个以上小角度的和,每个小角度用反正切表示,这样就可利用(2)式计算每个小角度进而算出 。一般地,取tan 为某个真分数,由

计算出

就可得到 恒等式 ,例如:令tan ,因

于是有 :

arctan arctan

用类似的方法,还可进一步将上面的arctan 分成两个更小的角度之和 。找到恒等式后,利用(2)式便可计算 了,例如:由上面的等式可得:

这样计算 就更简单, 收敛速度 也较(3)快 。

利用正多边形的面积

当圆的半径为1时,其面积恰好为 ,因此 必介于单位圆的内接正多边形面积与外切正多边形面积之间,因单位圆的内接正 边形面积为 sin ,外切正 边形面积为 tan ,故 sin tan ,下面利用上式来求 的近似值,为简便起见,令 ( ),则有 sin < < tan ,因此当用 sin 去近似 时,误差不超过 tan - sin ,记 = sin , = tan ,则有:

……,……,

=3.141592517 =3.141592722

=3.141592645 =3.141592658

=3.141592653 =3.141592653

……,……,

由此可见, 随 的增大而增大, 随 的增大而减小,且两者越来越接近 。

数学中的应用

是数学上的一个重要 常数 。随着近代数学不断深入的发展,逐步发现 与许多数学问题、数学公式有关 。

初等数学

在很多初等数学公式中出现,例如 圆周长 , 为半径 ;圆面积 ;椭圆面积 ;球体积 等 。

无穷

与 “无穷 ” 关系密切,其中 的无穷表达式主要指无限连分式 (数 )、 无穷乘积 、无穷 级数 、反正切式等 。

无穷乘积表达式

1650年, 英国 数学家Waliss得到 的无穷乘积表达式 :

反正切式

的无穷级数非常之多,其中反正切展开式使用率较高。

在两角和的 正切 公式tan =(tan tan )/( tan tan ), 中,令tan =a,tan =b,则 arctan , arctan 再取反正切得:arctan arctan arctan (1),设 , ,则(1)式为:

arctan arctan arctan (2),在(2)式中若令 ,则得

arctan arctan (3),若令 ,则得

arctan arctan arctan ,再由(3)式得

2arctan arctan ,如此递推下去,可得许多 的反正切式 。

级数

级数的定义为:给个一个数列 ,将其各项依次相加,得到的表达式

称为数项的级数或者无穷级数(简称为级数) 。

傅里叶级数

大多数周期信号都可以用 正弦 和 余弦 级数的展开式来表示,一个周期信号函数的傅里叶级数可表示为:

( cos sin )

cos

sin

式中: 为信号 基频 ,单位为Hz; 为周期信号的周期; 和 之间的关系可表示为 , 角频率 和 之间的关系可表示为:

利用傅里叶 级数 ,周期信号可以展开成无限多个 正弦 项和 余弦 项之和 。

其它级数

(1) (Leibnig级数)

(2) (Fousies级数)

数据来源:

积分

概率积分

(Euler- Poisson 积分)

Fejer积分

Fejer积分是指带Fejer核 的积分,这类积分在理论证明及工程应用上非常广泛,许多重要 定理 的证明、 习题 的解答都能看到它们的作用 。

柯西积分

若函数 在简单正向闭曲线 所围成的区域 内解析,在区域 的边界 上连续, 是区域 内任意一点,则:

函数

函数的定义为:设有两个变量 和 和一个实数集的 子集 ,若对于 中的每个值 ,变量 按照一定的法则有唯一一个确定的值 与之对应,则称变量 是变量 的函数,记作 ,其中数集 称为函数的定义域,相应的函数值的全体 称为函数的 值域 。

柯西分布

如果一个 随机变量 的密度函数形如:

则称该随机变量为服从参数为 的柯西分布 。

伽玛(高斯)函数

反常积分 ,在 时,是收敛的,因此当 取大于零的不同值时,它就有不同的值,所以它确定是一个 的单值函数,这个函数叫做伽玛函数,记作 ,或高斯 函数,记作 ,即:

伽玛函数在 工程技术 问题的研究中有广泛的应用 。

胡尔维茨zeta函数

的函数 方程 是 的函数方程在 为 有理数 时的一个特殊情形,如果 和 都是整数, ,那么对所有的 ,则有 :

cos

黎曼zeta函数

阿道夫·胡尔维茨 函数 方程 中,当 时,得到 黎曼 函数方程,即 :

cos 或者,等价地有:

sin 。

其他应用

概率论

7世纪中叶,概率论在 帕斯卡 和 皮耶·德·费玛 两人间书信来往中开始得到发展。它形成一门学问是1812年 皮埃尔-西蒙·拉普拉斯 发表的解析概率论,距今约180年的历史。把 与概率联系在一起,可以从小球命中圆、 蒲丰 氏问题(投针问题)着手 。

小球命中圆

在边长 的正方形中有个直径为 的圆, 投掷小球 ,问小球命中圆的概率是多少?其中命中正方形四周边线或命中 圆周 都算是命中圆 。

解:

正方形的面积是 ,正中圆的面积是 。故小球命中圆的概率 %,因为投球太容易,这样提出问题不太适合,因为瞄准当 中投 ,相当多的人是百发百中的,即概率应是100%,但为理论上即从数学上的考虑是 或28.26% 。

投针问题

法国 的科学家 蒲丰 (1707—1788年)于1777年提出了这样的一个问题,在画有距离为 的 平行线 的平面上,任意投一枚针,设针长为 ,试求针与任一条直线 相交 的概率 。

设 为针的 中点 到最近一条平行线的距离,显然易知: ,由几何概率定义,它等于阴影区域的面积与矩形面积之比,即 :

显然,当二平行线距离 与针长度 相等 时,针与直线相交叉的概率为 ,又设 为投针次数; 为 相交 次数,得: ,由实验结果就可得值 的近似值 。在历史上也有很多人用“蒲氏原理”进行过试验求出 的不同的近似值 。

介子理论

1935年, 日本 物理学家 汤川秀树 提出介子理论,他认为:就像 光子 是和 电磁场 及电磁力联系在一起的一样,U—场也存在一个“量子”;汤川认为 核力 有一个 — 厘米的作用球,在 — 厘米处,核力突然下降,核力可以看成是一个核子发射出一个量子,另一个核子吸收同一个量子,根据 不确定性原理 ,可以得出: ,其中 , ,代入公式得到:

当 取 厘米时,可求出此量子的质量约为电子质量的200倍,因为 中子 和 质子 、质子和质子以及中子和中子之间的作用是 相等 的,所以,这种粒子可以以三种形式出现:中性、带正电荷或负电荷,并且电荷大小等于质子的电荷 。

交流电路

欧拉公式 的几何意义是:当一个 向量 在复平面上以匀 角速度 绕原点转动时,矢量 端点 在实轴和虚轴上的投影分别作简谐振动,在 正弦 交流电路中 电流 可以表示为 sin ,由欧拉公式 cos sin ,把等式各项中的 代换为 ,得 :

cos sin

在交流电路理论中用电流 相量 来代表相应的旋转 向量 ,把 三角函数 运算转化为 复数 运算,从而简化计算过程,设三相四线制供电系统中的线电流为 :

sin ,

sin ,

sin ,

对三个以相同 角速度 沿逆时针方向旋转的 向量 ,可以利用 平行四边形 法则求出它们的合矢量,合矢量在坐标轴上的投影就代表合成的简谐振动;求出三个电流的合 向量 后,中线上中合成的 电流 也就确定了 :

电阻 元件、 电感元件 和 电容 元件上电压的合成结果可用同样的方法得到 。

国际圆周率日

2011年, 国际数学联盟 正式宣布,将每年的3月14日设为国际圆周率日。2019年, 联合国 教科文组织宣布3月14日为国际数学日。2020年3月,中国数学界以网络科普讲座的形式庆祝了第一个国际数学日 。

相关事件

2021年, 瑞士 科学家利用 超级计算机 历时108天,将圆周率计算出 小数点 后62.8万亿位数。

2022年, 谷歌 将圆周率计算到小数点后100万亿位,创下当时的 世界纪录 。 2023年4月, Solidigm 公司追平这一纪录。

2024年3月14日,总部位于 加利福尼亚州 的 计算机 存储公司Solidigm发布声明称,该公司已将圆周率Pi(π)计算到小数点后约105万亿位,打破此前100万亿位的世界纪录。本次计算历时75天,利用了100万GB数据,需要的计算能力与数十万部智能手机相当。

2025年12月,StorageReview团队使用单台戴尔PowerEdgeR7725服务器,历时4个月将圆周率计算至 小数点 后314万亿位,此次计算的能耗效率较高,每万亿位数字消耗13.70千瓦时,该能耗水平仅为此前300万亿位规模计算集群能耗的一小部分。

参考资料

  1. 圆周率与国际数学日 宜昌市科学技术局 2023-11-10
  2. 新纪录!这个数已精确到小数点后105万亿位 腾讯网 2024-03-16
  3. 历史人物 保定地方志 2023-11-14
  4. StorageReview Sets New Pi Record: 314 Trillion Digits on a Dell PowerEdge R7725 Storage Review 2026-06-17