利用泰勒级数计算无理数和以与其他任意无理数近似值.docVIP

利用泰勒级数计算无理数和以及其他任意无理数的近似值 适用章节：第十二章第四、五节 幂级数，泰勒级数（同济大学数学系《高等数学（第六版）》） 一、问题提出： 和是两个常用的无理数。是一个在 HYPERLINK /wiki/%E6%95%B8%E5%AD%B8 \o 数学 数学及 HYPERLINK /wiki/%E7%89%A9%E7%90%86%E5%AD%B8 \o 物理学 物理学普遍存在的 HYPERLINK /wiki/%E6%95%B0%E5%AD%A6%E5%B8%B8%E6%95%B0 \o 数学常数 数学常数，其定义为 HYPERLINK /wiki/%E5%9C%93%E5%BD%A2 \o 圆形 圆形之 HYPERLINK /wiki/%E5%91%A8%E9%95%BF \o 周长 周长与 HYPERLINK /wiki/%E7%9B%B4%E5%BE%91 \o 直径 直径之比，是精确计算圆周长、圆面积、 HYPERLINK /wiki/%E7%90%83_(%E6%95%B0%E5%AD%A6) \o 球 (数学) 球体积等量的关键值，的近似值是3.141592654。是高等数学中常用的无理数，的值约为 2.718282。它们是怎样计算出来的，大家可能不是很清楚，我们将利用幂级数以及泰勒级数展开法借助数学软件Mathmatica来计算它们的近似值。问： （1）计算的近似值。 （2）计算的近似值。 （3）怎样利用余项判断近似值的误差。 （4）考虑一下其它无理数的计算，比如等。 二、涉及知识点 幂级数，泰勒级数 三、所使用的软件和关键语句 软件：Mathmatica 5 关键语句：Series[f{,,}] %将函数在处展开到阶幂级数。 四、实现的过程和结果 1．提示：， 。 2．实验过程： （1）计算的近似值 在Mathmatica 5工作页面输入如下命令： n = 20 Series[ArcTan[x], {x, 0, n}] %在展开到=20阶 f[x_] = 4 Normal[%] %去掉幂级数余项且乘4倍 NumberForm[f[1.] , 10] %令得的近似值,输出10位有效数字 wucha = N[1/(2 n + 3), 10] %误差精度，小于该数值，输出10位数字 （结果与上面程序每一行对应，其中ArcTan[x]表示函数，Series[]级数展开函数） 3.041839619 0.02325581395 %精确到整数位 计算结果不太精确，我们增大的阶数，如下 n=100； Series[ArcTan[x], {x, 0, n}]; % 展开到100阶 f[x_] = 4 Normal[%]; NumberForm[f[1.] , 10] wucha = N[1/(2 n + 3), 10] 输出结果： 3.12159 0.004926108374 %精确到小数点后一位 n=1000； Series[ArcTan[x], {x, 0, n}]; % 展开到1000阶 f[x_] = 4 Normal[%]; NumberForm[f[1.] , 10] wucha = N[1/(2 n + 3), 10] 输出结果： 3.139592656 0.0004992511233 %精确到小数点后一位 n=10000； Series[ArcTan[x], {x, 0, n}]; % 展开到10000阶 f[x_] = 4 Normal[%]; NumberForm[f[1.] , 10] wucha = N[1/(2 n + 3), 10] 输出结果： 3.141392654 0.00004999250112 %精确到小数点后3位 n=100000； Series[ArcTan[x], {x, 0, n}]; % 展开到100000阶 f[x_] = 4 Normal[%]; NumberForm[f[1.] , 10] wucha = N[1/(2 n + 3), 10] 输出结果： 3.1415726