圆周率的由来

数学史期末作业（论文）No

1

题目：圆周率的由来、应用及历史作用

姓名：

班级：

学号：

数学史期末作业（论文）No

2

目录

引言

1圆周率的由来.........................................................4

1.1古希腊求值.......................................................5

1.2古中国求值.......................................................5

1.3伊斯兰求值.......................................................5

1.4现代求值.........................................................6

2圆周率的应用..........................................................6

2.1用圆周率来测试计算机性能...........................................6

2.2圆周率在C语言中的应用.............................................6

3圆周率的历史作用....................................................10

3.1通过找出各种表达式..............................................10

3.2通过计算圆的面积和周长..........................................10

3.3用来进行一些函数的定义，积分的计算，指数的构成..................10

数学史期末作业（论文）No

3

引言

众所周知，圆周率一般以来表示，是一个在数学及物理学普遍存在的数学常数．它

定义为圆形之周长与直径之比．它也等于圆形之面积与半径平方之比．是精确计算圆周长、

圆面积、球体积等几何形状的关键值．圆周率是一个常数（约等于3.1415926），是代表圆

周长和直径的比例．它是一个无理数，即是一个无限不循环小数．圆周率在生产实践中应

用非常广泛，在科学不很发达的古代，计算圆周率是一件相当复杂和困难的工作．俗话说

得好，“有理走遍天下，无理寸步难行”圆周率好比这个“理”．有了圆周率不仅解决

了困惑众多数学家的三大著名几何问题之一的化圆为方的不可能性更为后续的数学研究

奠定了基础．因此，圆周率的理论和计算在一定程度上反映了一个国家的数学水平．

本文通过对圆周率各个时期由来的认识，深刻的理解到圆周率的历史价值，包括通过

π找出各种表达式，通过计算圆的面积和周长，一些函数的定义，积分的计算，指数的

构成等都要用到；还介绍了圆周率在测试计算机性能上的应用和圆周率在C语言上的应

用，最后还详述了圆周率的历史作用。

数学史期末作业（论文）No

4

1圆周率的由来

很早以前,人们看出,圆的周长和直经的比是个与圆的大小无关的常数,并称之为圆周

率.1600年,英国威廉.奥托兰特首先使用表示圆周率,因为是希腊之“圆周”的第一个

字母,而δ是“直径”的第一个字母,当δ=1时,圆周率为.1706年英国的琼斯首先使用

.1737年欧拉在其著作中使用.后来被数学家广泛接受,一直没用至今.

是一个非常重要的常数.一位德国数学家评论道：“历史上一个国家所算得的圆周率

的准确程度,可以做为衡量这个这家当时数学发展水平的重要标志.”古今中外很多数学家

都孜孜不倦地寻求过值（如图1所示）的计算方法.

图1

1.1古希腊求值

公元前200年间古希腊数学家阿基米德首先从理论上给出值的正确求法.他用圆外

切与内接多边形的周长从大、小两个方向上同时逐步逼近圆的周长,巧妙地求得．

公元前150年左右,另一位古希腊数学家托勒密用弦表法(以1的圆心角所对弦长乘以

360再除以圆的直径)给出了的近似值3.1416.

1.2古中国求值

公元200年间,我国数学家刘徽提供了求圆周率的科学方法----割圆术]1[（如图2所示）,

体现了极限观点.刘徽与阿基米德的方法有所不同,他只取“内接”不取“外切”.利用圆

面积不等式推出结果,起到了事半功倍的效果.而后,祖冲之在圆周率的计算上取得了世界

领先地位,求得“约率”和“密率”(又称祖率)得到3.1415926<<3.1415927.可惜,祖冲

数学史期末作业（论文）No

5

之的计算方法后来失传了.人们推测他用了刘徽的割圆术,但究竟用什么方法,还是一个

谜.

正六边形正十二边形正二十四边形正四十八边形

图2

1.3伊斯兰求值

15世纪,伊斯兰的数学家阿尔.卡西通过分别计算圆内接和外接正32边形周长,把

值推到小数点后16位,打破了祖冲之保持了上千年的记录.

1.4现代求值

本世纪50年代以后,圆周率的计算开始借助于电子计算机,从而出现了新的突破.目

前有人宣称已经把计算到了亿位甚至十亿位以上的有效数字.

人们试图从统计上获悉的各位数字是否有某种规律.竞争还在继续,正如有人所说,

数学家探索中的进程也像这个数一样:永不循环,无止无休……

2圆周率的应用

2.1用圆周率来测试计算机性能

它现在可以被人们用来测试或检验超级计算机的各项性能，特别是运算速度与计算过

程的稳定性．这对计算机本身的改进至关重要．

就在几年前，当Intel公司推出奔腾

(Pentium)

时，发现它有一点小问题，这问题正是通

过运行的计算而找到的．这正是超高精度的计算直到今天仍然有重要意义的原因之一．

2.2圆周率在C语言中的应用

2.2.1简单技术公式计算圆周率—掌握C语言的循环]2[

C语言课程盗了循环章节，它的魅力就逐渐显现了出来，许多小程序和算法都可以让

学生去尝试时限．而我们在高等数学泰勒公式章节学习的圆周率公式这个时候就派上了用

场．

中学时我们就已经知道：

1

4

tan



，从而4arctan1

，如果我们应用泰勒公式将

arctan

展开，就可以得到

数学史期末作业（论文）No

6











12

)1(

753

arctan

12

1

753

n

xxxx

xx

n

n

把1x代入，就得到







12

1

)1(

7

1

5

1

3

1

1

4

1

n

n



这个简单的公式就可以用循环实现来实现圆周率的近似计算（虽然它并不真正适用于计算

圆周率），下面来分析一下．

我们用C语言中的double类型的变量来存放最终的结果，而上面公式中的n作循环变

量再适合不过了，在实际编码中我们取

n

从1变化到30000.由于double类型本身的精度限

制，

n

其实也没有必要取太大的值．在循环中，

n

每次增加2，而要注意类型转换和在循环

累加中项的符号变化，这里为了清楚起见，我们增加了变量s，由它产生了正负符号的变

化．下面来看代码：



()

.

#main

hstdioinclude

;1,1intsi

double;0d

)30000(iwhile



);1*2)(/(1*idoublesd

;1*s

;i



);*4,"%("intdnIfpifpr

该程序在VC++6.0中运行结果3.141559，达到预期要求．

作为C语言实验课课后练习，可尝试进行下面公式的算法实现：



17

1

13

1

11

1

7

1

5

1

1

3





17

1

13

1

11

1

7

1

5

1

1

6

3



数学史期末作业（论文）No

7



222

2

5

1

3

1

1

1

8





222

2

3

1

2

1

1

1

6



这样锻炼了C语言能力，同时也加深了数学能力的培养，一举两得．

2.2.2圆周率的概率算法—掌握C语言的函数

随机模拟方法就是这样一种计算圆周率的方法．虽然这种方法并不能“非常好”的逼

近圆周率，但是在其他领域，随机模拟具有非常大的作用．这种算法也称为蒙特卡洛算法．

蒙特卡洛算法简单描述：在树枝积分法中，利用求单位圆的

4

1

的面积来求得

4



，从而

得到．单位圆的

4

1

面积是一个扇形，它是边长为1单位正方形的一部分．只要求出扇形面

积

1

S

在正方形面积S中占的比例

S

S

K1

就立即能得到

1

S

，从而得到的值．办法是在正方

形中随机投入很多个点落在扇形内．将落在扇形内的点数m与所投点的总数n的比

n

m

作为

k的近似值．

P

落在扇形内的充要条件是

122yx

．

代码：







hmathinclude

hstdlibinclude

hstdioinclude

.

#.

#.

#void

()main



doublei;int

;,,,nkjl

;0n

);30000;0(iiifor



;0.32767/()randj

;0.32767/()randk

);**(kkjjsqrtl

;)1(nlif



);/*4,"%("intinnIfpifpr

数学史期末作业（论文）No

8



运行结果：π=3.126167

由于没初始化随机种子，所以该寒暑在执行过程中并不具有“随机性”，由于在TC和

在VC++6.0中初始化随机种子的函数并不相同，本程序为保持兼容性，未加入那些内容．

2.2.3“外星人计算



的程序？”—C语言阅读理解

在网上流传很久的“精简”代码，如果你进入谷歌的主页，并且搜索“外星人计算



的

程序”这个词汇，很容易看到这样一个程序：

d%aed/a),e,4d"printf("%.14,c-2;*cg0),dfor(;a/5;])f[bc;-bfor(;){Main(

g;f[2801],e,d,2800,cb,10000,aInt

stdio.h"include"

#



0;returnb;d*-b;--,-gd/-g,-d%f[b]a,*[b])fdc;For(b

该程序当然违背了C语言的代码规范，不过如果你运行它，你会惊奇的发现它得到了

圆周率小数点后面800位数字．这段代码来自“国际C语言混乱代码大赛”（The

InternationalObfuscatedCCodeContest）．

3圆周率的历史作用

3.1通过找出各种表达式

1579年法国韦达发现了关系式，首次摆脱了几何学的陈旧方法，寻求到了的解析

表达式．1650年瓦里斯把表示成无穷乘积，无穷连分数，无穷级数等各种值表达式纷纷

出现，值计算精度也迅速增加．稍后，莱布尼茨发现接着欧拉证明了这些公式的计算量都

很大．尽管形式非常简单，值的计算方法的最大突破是找到了它的反正切函数表达

式．1706年英国数学家麦欣首先发现了其计算速度远远超过方典算法．

3.2通过计算圆的面积和周长

某个古代文牍员以不同长度的半径画了一些圆，他取了每个圆的直径（将半径加倍）

只是为了好玩．他决定以每个圆的直径为单位长度在圆周上丈量．令人惊奇的是，不管圆

的大小如何，圆周总是直径的3倍多一点．由于与圆的特殊关系，故数学家设计用来计

算出圆的面积和周长的新方法．

例已知一个圆形花坛的直径是4米，沿它的外侧铺一条1米宽的小路，求这条小路的面

积。（精确到0.1平方米）

数学史期末作业（论文）No

9

解花坛半径是:224（米）

15.7223.14-1)(21)(23.14

（平方米）

3.3一些函数的定义，积分的计算，指数的构成等都要用到

例如，1777年，法国数学家蒲丰研究投针问题，将一根长为l的的针任意投到画有间

距为l)a(a的平行线的平面上，他得到得结论是：该针与任一平行线相交的概率是

a

p

21

，圆周率与随机现象产生了密切联系即在概率中也有作用．在数学中还有一个重

要公式2)

1

1log(4

i

i

i





将圆周率与虚数单位

i

联系起来．

随着数学的不断发展，π的应用不再局限于求圆的面积和周长，椭圆，萁舌线，旋

轮线等面积公式中也都出现了π值。此外，一些函数的定义，积分的计算，指数的构成等

都要用到π。例如，1777年，法国数学家蒲丰研究投针问题，将一根长为l的的针任意投

到画有间距为a(a>l)的平行线的平面上，他得到得结论是：该针与任一平行线相交的概

率是p=2l/aπ，圆周率与随机现象产生了密切联系即π在概率中也有作用。在数学中还有

一个重要公式π=4log(1-i/1+i)^i/2将圆周率与虚数单位i联系起来。1740年，欧拉进

一步得到关系式e^iπ+1=0，将数学中5个重要的数学最重要的两个运算符号统一在一个

公式中，令人拍案叫绝！在数论中，法国人沙特尔1904年得到一个定理：任一写下两个

整数，则它们互素的概率是6、π，一个简单的圆周率π几乎无所不在。

背诵圆周率能够锻炼人的记忆力，我国桥梁专家茅以升年轻时就能背诵圆周率锻炼记

忆力。晚年时仍能轻松地背出圆周率的100位数值。

可见圆周率π不仅与我们身边的数学紧密相连更与我们的生活息息相关。俗话说得

好，“有理走遍天下，无理寸步难行”圆周率π就好比这个“理”。有了圆周率π不仅解决

了困惑众多数学家的三大著名几何问题之一的化圆为方的不可能性更为后续的数学研究

奠定了基础。