$ \newcommand{\uni}[1]{\ \mathrm{#1}} \newcommand{\E}[1]{\times 10^{#1}} \newcommand{\ket}[1]{|#1\rangle} $

萌狸
狸群 【圆周率节】祖冲之是如何计算 π 的

yxlllc
7:02:50
今天是第30个圆周率日
栗子球
7:03:11
30年前发生了什么
yxlllc
7:03:15
3b1b不打算出个关于pi的科普吗
xyfd
7:03:54
无数人做过的东西
yxlllc
7:04:46
去年我在票圈里科普过祖冲之当时是怎么算pi的
yxlllc
7:14:45
大概是这样
悟理萌主
7:16:50
悟理萌主
7:16:56
好复杂的样子
yxlllc
7:18:52
是啊,这可是大概1500年前
yxlllc
7:19:06
就这个领先西方将近1000年
IAsimov
7:40:25
为啥霍金去世,这么多人发朋友圈和动态表示很难受
yxlllc
7:41:22
本质上是在为自己难受
悟理萌主
7:42:06
我不难受
yxlllc
7:42:20
有句话说得非常好忘了是谁说的了,这个世界只有两种人,一种是英雄,一种是崇拜英雄的人
IAsimov
7:43:57
我就挺乐观
IAsimov
7:44:48
大师留下很多财富,与其无病呻吟,不如追随他们的脚步,完成他们未竟的事业
yxlllc
7:44:53
为自己的平庸而难受
IAsimov
7:45:55
平庸和平凡还是有点区别的
栗子球
8:06:31
(站在巨人的肩膀
subsophiec
8:17:50
膜大佬
 
Probe
14:07:29
https://www.youtube.com/watch?v=bcPTiiiYDs8  
Probe
14:07:51
yxlllc
14:10:10
这波可以
Probe
15:02:51
lc 之前应该先预测一下分钟物理出的相对论
Probe
15:03:09
这下lc 的预言家身份坐实了
栗子球
15:04:52
因果律武器
yxlllc
15:05:29
不是,pi节是真的啊,确实有这么个东西
yxlllc
15:06:01
3b1b既然做数学科普,当然完全有可能来这么一波
yxlllc
15:06:12
一切都是很符合逻辑的
Probe
15:07:16
我感觉pi节是不是近年才炒起来的的啊
Probe
15:07:25
至少高中我没怎么听说过
N31Welt
15:07:47
http://mp.weixin.qq.com/s/kdPfp0o9YyBdXxrnqrXj2g
氘化氢
15:08:07
没有吧
yxlllc
15:08:19
yxlllc
15:08:27
人家谷歌都知道这个节日
N31Welt
15:09:00
在1988年3月14日那天,Larry Shaw带着博物馆的员工以及其他参与者,一起绕着博物馆的纪念碑转了22/7(π的近似值)圈,还一起吃水果派、分享与π有关的知识。
yxlllc
15:09:04
谷歌这波宣传看来是全世界都知晓了
N31Welt
15:09:52
此后,每一年的3月14日,旧金山科学博物馆都会举办π的庆祝活动。
yxlllc
15:10:05
所以中国的科普事业真的落后太多了
yxlllc
15:10:19
人家的科普已经成为文化了
栗子球
15:10:25
(我才知道
小灰机
15:10:28
呃。。
N31Welt
15:10:49
派节历史估计比群里很多人年龄大
小灰机
15:10:51
pi节不是人尽皆知吗。。
Probe
15:11:01
每年pi节mathematica还打折的
小灰机
15:11:01
比我大。
Probe
15:11:17
我真是最近两三年才听说的
yxlllc
15:11:46
我大概初中的时候就知道了
小灰机
15:12:16
初中知道+1
小灰机
15:12:30
N31Welt
15:13:44
我们创个e节吧
小灰机
15:14:12
貌似世界上并没有公认的e节诶。。
N31Welt
15:14:12
每年2月7号,下午6点28开始
yxlllc
15:14:13
所以感觉人家外国人在科普祖冲之的时候我们自己人却不知道这么一回事感觉很不好
N31Welt
15:14:30
@小灰机 所以从我们开始
小灰机
15:14:35
祖冲之,小学奥数书上随处可见
yxlllc
15:15:09
然而很少有人讲他具体是怎么算的pi
小灰机
15:15:19
有啊
小灰机
15:15:42
不就是把圆近似成多边形算面积的吗?
N31Welt
15:15:55
内接n边形,老祖非常能算
yxlllc
15:16:03
不只是怎么简单
N31Welt
15:16:14
小灰机
15:16:18
???
小灰机
15:16:22
啥意思
氘化氢
15:16:23
2月71号。
N31Welt
15:16:37
每年2月7号,下午6点28开始
yxlllc
15:16:52
比如割圆术是刘徽发明的,很多人都不知道
yxlllc
15:17:04
论贡献实际上刘徽的功劳最大
小灰机
15:17:20
但祖冲之算啊
yxlllc
15:17:32
理论体系是他构建的,祖冲之按照这个方法计算
小灰机
15:17:36
刘徽难道不讲么。。
小灰机
15:18:01
我记得《华罗庚题典》里面就有好像。。
xyfd
15:18:12
http://blog.sciencenet.cn/blog-100379-1081681.html
xyfd
15:18:23
xyfd
15:18:38
我觉得这个问题重点在于他使用了怎样的计算技巧
xyfd
15:18:46
总不能强算吧……
xyfd
15:19:55
yxlllc
15:20:02
这种明显都是上帝视角了,那个年代哪来的无穷级数
xyfd
15:20:23
怎么就没有……
xyfd
15:20:34
中国古代算数特别强
yxlllc
15:21:50
没有类似的概念
yxlllc
15:22:17
哪怕是三角函数都几乎没有
xyfd
15:22:37
小灰机
15:23:05
@xyfd 中国弱在没有完整的理论体系
xyfd
15:23:49
有没有理论体系和计算pi有啥关系……
xyfd
15:24:01
反正就是算到那个精度了
yxlllc
15:24:08
反正经过我的研究就是这么算的
xyfd
15:24:23
你估计过这个工作量了吗
氘化氢
15:24:48
很大
氘化氢
15:24:53
但是可以算
yxlllc
15:24:54
反正那个年代不可能写成sinx 这种东西
氘化氢
15:25:02
当时算数很发达的
yxlllc
15:25:07
压根就没这种概念
N31Welt
15:25:23
祖冲之割到16384边形,实在算不动了,卒
xyfd
15:25:50
yxlllc
15:26:41
我们确实可以用高等数学的方法来理解割圆术的本质
xyfd
15:26:46
注1:据说巴比伦人在公元前2000年就会算平方根,他们这样算根号二,x_{n+1} = (x_n+2/x_n)/2。这其实就是现代的牛顿法,收敛超快!我们日常计算机就应该是用这个算法。不确定祖冲之是否用的这个聪明算法。

注2:据说巴比伦人也会解一元二次方程。这可以作为一个文明发展的标志。曾经有个老师跟博主讲蒙古人如何以野蛮战胜文明,他就说成成吉思汗灭花剌子模(就是阿尔-卡西的国家)的时候,花剌子模都已经能够解一元二次方程了。其实不止能解一元二次方程,他们还把pi算到了16位有效数字。

注3:我从来没见过祖冲之
的圆周率的原始出处,据说在隋书-律历志。有这本书的同志,希望能够拍个照分享下,谢谢!

注4:博友尤明庆老师也研究过同样的问题。不过,他只做一次变换,没有对变换所得数列再作变换,所以需要割到192边形。

注5:其实一个更大的迷是,祖冲之是如何找到pi的有理近似355/113的。
yxlllc
15:27:59
至少刘徽是按我这么算的
wuyudi(1160020387)
15:28:55
拉马努金如何发现那些公式对于我永远是个谜
yxlllc
15:29:06
主流观点认为祖冲之也是采用刘徽这种方法来算的,如果承认这个那祖冲之得出的pi近似值就是我上面那个式子
小灰机
15:30:43
xyfd
15:31:19
小灰机
15:32:14
对啊,Intel的SuperPi比AMD强很多啊
yxlllc
15:38:11
对应该是这个,从正六边形开始算,我刚才那个是从正方形开始算