24无穷级数.pptVIP

2.4 无穷级数 2.4.1 级数部分和的变化 1、符号求和函数symsum 1)、symsum(S)：其中S为待求和的级数的通项表达式.其功能是求出通项为S的级数关于系统默认变量的有限和. 2)、symsum(S，v)：其中S为待求和的级数的通项表达式，v为求和变量.其功能是求出通项为S的级数关于变量v的有限和(例如v从0到k-1的有限和). 3)、symsum(S，a，b)或symsum(S，v，a，b)：求从a到b的级数的和.其中b可以取有限数，也可以取无穷.此命令既可以用于求级数的部分和，也可以用于判别级数的收敛性. 2、数值求和函数sum 3、vpa(f,n) 指定f 的运算精度 2.4.2 常数项级数敛散性 1、符号求和函数symsum； 2、求极限函数指令limit； 3、绝对值指令abs; 2.4.3 幂级数 1、泰勒(Taylor)幂级数展开指令 taylor( ),具体功能及应用实例见§2.2节； 2、求极限指令limit； 3、求级数和指令symsum； 4、多变量函数的泰勒(Taylor)幂级数展开, 主要以调用Maple语言的mtaylor( )函数来直接求取多变量函数的泰勒(Taylor)幂级数展开. F=maple(mtaylor,f,[x1,x2,…,xn],k)：根据原点展开； F=maple(mtaylor,f,[x1=a1,…,xn=an],k)：根据(a1,…,an)点展开. 其中，k-1为展开的最高阶次，f为原多变量函数. 2.4.4 傅里叶级数 1、积分指令int( )； 2、矩阵的点乘积； 3、了解在设计Matlab绘图时的图形标注函数及其功能. * 解： 1）先用数值计算方法求值 format long;sum(2.^[0:64]) 由于数值计算中使用了double数据类型，至多只能保留16位有效数字，因此结果并不很精确。若利用符号求和指令可以求出精确的结果。 syms k; symsum(2^k,0,64) 2) n=0:50;S1=sum(1./(2.^(2*n-1).*(2*n-1))); format long;S1 syms n; S2=symsum(1/(2^(2*n-1)*(2*n-1)),0,50) 3）syms n x;S50=symsum((-1)^(n+1)*x/(n*(2*n-1)),n,1,50) 解：本题是求解级数与极限的综合问题，在命令窗口中直接输入下一语句求解。 syms k n; lim=limit(symsum(1/k,k,1,n)-log(n),n,inf) 得出的结果为Euler常数。其前100位数字为 vpa(lim,100) 解：1）输入语句 syms n; symsum(1/((2*n-1)*(2*n+1)),1,inf) 2）输入语句 syms n; symsum(1/n^2,1,inf) 解：对于级数，当n无限增大时，它的一般项不趋于零， 则级数发散。 1）syms n; limit((n/(2*n-1))^2,n,inf) ans = 1/4 该级数发散。 2） syms n; limit((n^(n+1/n))/(n+1/n)^n,n,inf) ans = 1 该级数发散。 解：由于所运用的代码比较简单，没有必要编写M文件. 1）在命令窗口输入 syms n; symsum(1/n^(3/2),1,inf) ans = zeta(3/2) 故级数收敛，且其和为zeta(3/2)。 该值约为：vpa(zeta(3/2),10) ans = 2.612375349 2) 首先输入 syms n; symsum(sin(1/n),1,inf) ans = sum(sin(1/n),n = 1 .. Inf) 从结果中并不能看出敛散性，若用比较判别法， 输入程序： limit(sin(1/n)/(1/n),n,inf) 得到以下结果 ans = 1 所以是发散的 *