幂级数在近似计算中的应

1、论文4 幂级数在近似计算中的应用谢文清 江权霞 (指导老师:陈引兰)数学与统计学院1001班摘要：形如的函数项级数称为幂级数,幂级数可以看成是一个“无限次多项式”，它无论在理论上还是实践上都是一个有力的工具.本文主要运用幂级数的展开式,对无理数等，利用计算机相关软件，进行近似计算. 关键词：幂级数、近似计算1.理论依据以某个幂级数展开式为基础，然后把所需要求的量表达成级数的和，并依据要求，选取部分和作这个量的近似值，误差用余项估计.我们先给出一些基本初等函数的幂级数展开式及它们对应的余项2的近似计算 本节利用两个函数的幂级数展开式来近似计算，在相同的误差条件下，取不同的，若取级数的前n项和作为

2、的近似值，对应的n值不一样，这就为幂级数在近似计算中的应用提供了很大的空间.由函数的幂级数展开式知若取时， （1） 等式的右端是一个交错级数且是收敛的，实际计算时，我们只能使用有限项。如果取级数前n项之和作为的近似值即，其误差为，为了保证误差不超过，就要取级数（1）的前20000项进行计算，计算量之大可以想象.它的收敛速度很慢.对于展开式而言，当越小收敛越快，恰恰在端点收敛最慢. 以下取的求和的级数相应它的收敛速度要稍快些.现若取带入展开式得 （2） 若取级数的前n项和作为的近似值，其误差为 下面实现(2)式的计算，若要求误差小于(计算的程序见附录1)当n=8时，现取，显见，记，而，所以，就是

3、 (3)下面实现(3)的计算，若要求误差小于(计算的程序见附录2)当n=7时， 对于，误差一样（如要求误差小于），取不同的，对应部分和的项数n与近似计算的值如下表1n20000873.141673.14156对于的展开式而言，取 (4)下面实现(4)的计算，若要求误差小于(计算的程序见附录3)当n=4时，综上，知当误差确定时，对相同的幂级数展开式，x的取值不同，所取部分和的项数不同，近似计算的值也不同，对不同的幂级数展开式结果亦然.当然，当误差改变时，我们同样可以利用幂级数展开式估算出的值，其精确度更高.3.数的近似计算当=1时，所以取作为近似值，则误差为.例如:精确到，则需要(见附录4).扩

4、广:利用幂级数推导是无理数.反证法:假设是有理数，则等式左边是一个整数，右端第一项是整数，而k是小数；即右端不是整数，矛盾.故是无理数.4.对数的计算 利用对数的幂级数展开式，作对数的近似计算。根据对数的特征，只要计算出正整数的特征，那么由对数的运算，其它有理数的对数也就知道了.以ln(1+x)的麦克劳林级数作为出发点当时,当取前n项作为其近似值，其误差如要精确到就要截取一万项来计算，另外上面的展开式的收敛域为,这就不能直接用它来计算其它整数的对数.下面用一个收敛较快的幂级数来计算利用的幂级数展开式令,即带入(5),有估计余项如下如要精确到，即使,只要 (见附录5）拓展:令，有 这是一个递推公

5、式，所以据此可求任何正整数的对数，相应的也可求有理数的对数.如:当N=2时，即，有(的结果见附录6) 当N=4时，即，有(的结果见附录7)如此进行下去，可得ln6,ln7,的值利用上述计算方法，通过换底公式，我们可以计算得到了的一些近似计算结果并与数学用表中值进行比较（见表）表 的幂级数近似计算结果与数学用表中数值的比较12345678910幂级数 算00.301030.477060.602060.69870.778090.845040.900900.954121数学用 表00.30100.47710.60210.69900.77820.84510.90310.95421通过此表，知幂级数作为

6、近似计算的工具，结果与真实值很相近.参考文献1 董延闿.级数M.上海:上海科学技术出版社,1982.2华东师范大学数学系.数学分析.M.北京:高等教育出版社,19993周晓阳.数学实验与Matlab.武汉:华中科技大学出版社,2002附录1.s=0;n=1;ps=pi;while abs(s-ps)1e-4 s=(-1)(n-1)*2*3(1/2)/(2*n-1)*3(n-1)+s; n=n+1;ends,n程序所得结果为s=3.14167431n = 8 即为使计算结果精确到小数后第四位，只需求对应级数前7项的和利用Matlab软件算得 syms k symsum(-1)(k-1)*x(2*

7、k-1)/(2*k-1),k,1,8)ans = x-1/3*x3+1/5*x5-1/7*x7+1/9*x9-1/11*x11+1/13*x13-1/15*x15syms k f=6*(-1)(k-1)*(1/sqrt(3)(2*k-1)/(2*k-1)symsum(f,k,1,7)结果为ans =3.141674312. s=0;n=1;ps=pi;while abs(s-ps)1e-4 s=4*(-1)(n-1)/(2*n-1)*1/2(2*n-1)+1/3(2*n-1)+s; n=n+1;ends,n计算结果为s =3.14156158n = 73. s=3;n=1;ps=pi;while abs(s-ps)1e-4 s=(2*n-1)!/(2*n)!*(2*n+1)*2(2*n+1)+s; n=n+1;ends,n计算结果为s=3.14115 n=44. ff=sym(n*n!=107);solve(ff,n)ans = 10先算syms k nsymsum(1/sym(!),1,10)ans = 1.7182818则5. ff=sym(4*(2*n+1)*3(2*n-1)=104);solve(ff,n)ans = 4syms k nsymsum(2/3(2*-1)*(2*-1),1,4)a