MathStudio使用方法入门6连分数.ppt

MathStudio for iPad 使用方法入门 （6） cFrac Convergents 连分数 渐近分数 2014年3月14日 ;什么是连分数？连分数有什么用？ 初次看到连分数时，觉得它太麻烦，有些望而生畏；后来看到华罗庚的 《从祖冲之的圆周率谈起》，深入浅出地讲了连分数的用途和算法，觉得它很有 用，也很好玩，只是觉得用辗转相除法算起来很费力费时，因而就敬而远之了。 近来，发现MathStudio里的内置函数cFrac可以把实数转化为连分数，手指轻点几 下，就代替了大量纸笔之劳，太方便了！ 什么是连分数呢？ 若a0为整数，a1、a2、a3、a4… 皆为正整数，则下式就是简单连分数;连 分 数 的 基 本性 质 任何一个有理数一定可以唯一地写成有限连分数； 任何一个无理数一定可以唯一地写成无限连分数； 利用实数的连分数表达式的逐次截断值可以求出该实数的近似值—渐近分数。 连分数的逐次截断值从左、右两个方向交叉地逐次逼近真值。 a0 a2 a4 …… a …… a5 a3 a1 上式中，a的下标是偶数的项是递增逼近真数，下标是奇数的项是递减逼近真数 连 分 数 有 什 么 用？ 对上述基本性质我们不究其理，求尽其用。对函数连分式也不作探讨。 连分数的上述性质，和我们日常生活有关吗？ 有！比如说： 一. π值的计算，我们怎么得到祖冲之的约率、密率？ 二. 现在的历法为什么要四年一闰？一百年要少一闰？四百年又要加一闰？ 三. 阴历的大小月和闰月是怎么安排的？ 四. 日月蚀的发生有什么规律？ 五. 火星大冲、火星“逆行”现象有什么规律？ 诸如此类，凡是用有理数最佳逼近实数的问题以及自然界声、光、电、水、气等 周期性波动现象的研究，连分数都可作为计算工具来运用。 ; 一. π 值 的 计 算 古之九数，圆周率三，圆径率一，其术疏舛。自刘歆、张衡、刘徽、王蕃、皮延宗之徒，各设新率，未臻折衷。宋末，南徐州从事史祖冲之，更开密法，以圆径一亿为一丈，圆周盈数三丈一尺四寸一分五厘九毫二秒七忽，朒数三丈一尺四寸一分五厘九毫二秒六忽，正数在盈朒二限之间。密率，圆径一百一十三，圆周三百五十五。约率，圆径七，周二十二。又设开差冪，开差立，兼以正圆参之。指要精密，算氏之最者也。所著之书，名为《缀术》，学官莫能究其深奥，是故废而不理。 《隋书》卷十六 志第十一 律历上 求算圆周率的值是数学中一个非常重要也是非常困难的研究课题。中外古代许多数学家都致力于圆周率的计算。公元466年，祖冲之将圆周率算到小数点后7位的精确度，这一纪录在世界上保持了一千年之久。密率355/113 （3.1415929...）则是一个易于记忆（三个连续奇数：113355）的分数近似值，它是分母小于16604的分数中最接近π的。（引自维基百科） ;MathStudio的cFrac可以得到实数的连分式， 渐近分数怎么算？ 下式引自华罗庚《从祖冲之圆周率说起》 ;已知连分数的部分商，计算渐近分数 左图 编一个小程序，计算π连分数的前10个截断值 数组a(1)～a(10) cFrac计算出的π的10个部分商 p(1)、q(1) 、p(2) 、q(2)分别是前2个截断值3/1、 22/7的分子和分母 c(1)～c(10) 计算出π的10个截断值，受小数位的 限制，后六位的数值都相同了 从eval(p)和eval(q),得到 π的10个渐近分数： 3 22/7 333/106 355/113 103993/33102 104348/33215 208341/66317 312689/99532 833719/265381 1146418/364913 ;二. 为什么阳历要四年一闰？ 一百年要少一闰？四百年又要加一闰？ 历法是一种推算年、月、日的时间长度和它们之间的关系，制定时间的序列的方法，它是长时间的记时系统。 历法能使人类确定每一日在无限的时间序列中的确切位置并记录历史。 回归年（地球绕太阳公转一周） 是365日5时48分46秒余。比历法年365天多了5时48分46秒，即0.242199…天。 左图是回归年与历法年（平年）的