级数部分和的变化.docVIP

实验一 级数部分和的变化 实验的目的 1、理解掌握无穷级数的基本概念和性质； 2、掌握利用Matlab求无穷级数的部分和的基本方法； 实验的基本理论与方法 1、常数项级数： 一般地，如果给定一个数列，则由这个数列构成的表达式 4-1) 叫做(常数项)级数，记为即 = 其中第n项叫做级数的一般项。 作级数(4.1)的前n项的和 4-2) 当n依次取1，2，3，…时，它们构造成一个新的数列： ，，，…，，…。 如果级数4-1)的部分和而成的数列有极限，即，则称级数4-1)收敛，这时极限叫做级数4-1)的和，并写作；如果没有极限，则称级数4-1)发散。 2、级数的基本性质 性质1 如果级数收敛于和，则它的各项同乘以一个常数所得的级数也收敛，其和为。 性质2 设有两个级数： 则级数 也收敛，且其和为。 性质3 在级数的前面部分去掉或加上有限项，不会影响级数的收敛性或发散性，不过在收敛时，一般来说级数的和是要改变的。 性质4 收敛级数加括弧后所成的级数仍然收敛于原来的和。 3、级数收敛的必要条件 对于级数，它的一般项与部分和有如下关系 。 假设这级数收敛于和，则 因此级数收敛的必要条件：当n无限增大时，它的一般项趋于零，即 实验使用的函数与命令 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 实验指导 例1 求级数的下列部分和 1) ；2) ；3) ；4) ； 解： 1）先用数值计算方法求值。由于操作比较简单，现在命令窗口中输入以下代码，并按回车运行，得到结果。 format long;sum(2.^[0:64]) ans = 3.689348814741910e+019 由于数值计算中使用了double数据类型，至多只能保留16位有效数字，因此结果并不很精确。若利用符号求和指令可以求出精确的结果。 syms k; symsum(2^k,0,64) ans = 36893488147419103231 2) 先用数值计算方法求值，然后再利用符号求和。 %数值计算 n=0:50;S1=sum(1./(2.^(2*n-1).*(2*n-1)));format long;S1 S1 = -1.45069385566594 %符号求和 syms n; S2=symsum(1/(2^(2*n-1)*(2*n-1)),0,50);S2 S2 = -91060996588581738673777918971180129917604150160442795459686302540617/62770650218808522028345476945481102394317880334936508870080174489600 3）利用分别用symsum进行两种符号求和。 syms k n;S3=symsum(3^(k+1)/2^k,k);S4=symsum(3^(k+1)/2^k,k,0,n-1); S3 S4 输出结果 S3 = 6*(3/2)^k S4 = 6*(3/2)^n-6 4）输入程序 syms n x;S50=symsum((-1)^(n+1)*x/(n*(2*n-1)),n,1,50) S50 = 6118312683686597643267023636391449051611/6972037522971247716453380893531230355680*x 该题主要利用了级数的基本性质进行运算。 例2 求极限 解：本题是求解级数与极限的综合问题，在命令窗口中直接输入下一语句求解。 syms k n;lim=limit(symsum(1/k,k,1,n)-log(n),n,inf) lim = eulergamma 得出的结果为Euler常数。其前100位数字为 vpa(lim,100) ans = 0.577215664901532860606512090082402431042159335939923598805767234884867726777664