的两种新的快速逼近算法.docVIP

【标题】?的两种新的快速逼近算法 【作者】熊 俪 【关键词】?圆外切正多边形??圆内接正多边形???边长??弧长??? 【指导老师】邓 淙 【专业】数学与应用数学 【正文】1引言众所周知，圆周率?是数学中最重要和神奇的常数之一。2000多年来，人们对它进行了深入研究，发现了它的许多奇妙性质?。我国古代对?的研究，曾经作出过辉煌的贡献。追溯?的研究历史，从古自今不外两个方面，一是从理论上探索?的各种神奇性质，二是计算?的数值，现在人们已将?算到了小数后上万亿位，而且还在继续。当然，两个方面不可能截然脱离，而是互相联系，互相促进的，但又各有所向。例如，对?的数值计算就必须以理论得到的公式和方法为支撑。我国三国时期著名数学家刘徽曾创立?“割圆术”，以圆内接正多边形面积逼近圆面积，提出了计算?值的一种基本方法。祖冲之当年很可能就是利用刘徽的割圆术算到圆内接24576边形，得到了3.1415926＜?＜3.1415927?，领先于西方约1000年。随着数学的发展，人们还提出了计算?值的许多新方法，如“解析法”?、“概率法”?、“物理法”?等，真是令人耳目一新。到了上世纪70年代末，复旦大学俞文鮆教授等提出了“抛物线割圆术”?，从圆内接及外切正六边形算起，算到圆内接及外切正384边形就得到了?≈3.1415926。到90年代，邓淙教授?利用泰勒公式构造了以圆内接正多边形的边长逼近?的近似公式，提出以高精度计算?值的新方法，算到圆内接正24边形，就得到?≈3这些结果说明，现代数学赋予“割圆术”新的活力。同时也似乎在提示我们，如果类似地考虑圆外切正多边形，也应有类似的成效。受这些成果启发，本文探索了用圆外切正多边形边长逼近?值，和同时应用圆内接、外切正多边形边长逼近?值的方法，也得到了?的两种新的快速逼近算法。2已知成果的概述2.1 Ch.Huygens利用弦长求圆弧弧长的近似公式[11]如右图,圆弧?的长为?，?是?所对的弦长。而?是半弧所对的弦长。设弧长?可近似地表为???????????????????? 此处A、B是待定系数.若R是圆的半径，而?是弧?所对的圆心角，则有?（?）类似的，把?换成?，就有????（?）由此???????????? 可是?，自然必须选取A及B使?????????????????????????????因为这样一来，在所考察的公式中，左端与右端的差将仅只是含有?的项了.由上二式解得A及B的值为?，于是便得到Ch.Huygens利用弦长求圆弧长的近似公式??很易看出，其误差?可估计为：??????????????????????????????2.2?对Ch.Huygens弧长近似公式的推广1993年，文[10]将Ch.Huygens的弧长近似公式推广到了一般情形，简要过程如下(如图2)设圆弧?的长为L，半径为R，?所对的圆心角为?。图中?????????????????=2?=???=2??以?（i=1，2，???，m）表示?所对的弦长，?为已知，假设弧长L可以用???近似地线性表示，即???????????????????????????????????????????????????????????????????? L≈?=?????????????????(1)其中?为待定系数，再应用泰勒公式将弦长?展开，并代入（1）式得????????????????????????（2）用?代替?的误差为????????????????????????????????????????????（3）取常数?，使（2）中?的系数为1，而其余?的系数为0，并把得到的方程组化简得????????????????????????????????????????????（4）方程组（4）的系数行列式实际是一个?阶的范德蒙行列式，可算出常数?为????????????????????????（i=1,2,?,m）???????????????????????(5)把由（5）式所得的诸?代入（1）式，我们便得到利用弦长?,?,??求弧长的一般近似公式??????????????????????????????????????????????????????(6)???(5)式算出的诸A?都是既约分数,为了应用方便,还可以把它们化为同分母的分数.事实上,(5)式可改写为?,?? i=1,2,…..,m????????????????(7)（7）中的诸A?有共同的分母?，故（6）式可以改写为??????????????????????????????????????（8）（7）和（8）中遇到的积式?定义为1。