试验六一元函数的极值与数据拟合.doc

实验八 无穷级数 实 验 目 的 1．会进行级数的求和与收敛性判别 2．会求函数项级数的求和及其收敛域。 3．掌握函数的幂级数展开及应用 4．了解周期函数的傅立叶级数逼近 实 验 步 骤 进行级数求和运算，必须在运算命令前打开mathematica的程序包symbolic.m，格式如下： algebra\symbolic.m 一、数项级数求和 格式： Sum[f(n),{n,a,b}] 求 例：（1） （2） （3） （4） （5） 输入： algebra\symbolic.m Sum[1/(n*(n+1)),{n,1,Infinity}] Sum[1/n^2,{n,1,Infinity}] Sum[(-1)^(n+1)/n,{n,1,Infinity}] Sum[1/n,{n,1,Infinity}] Sum[n!/n^n,{n,1,Infinity}] 输出：（1）1 （2）Pi2/6 (3) Infinity (发散于无穷大) Sum[n!/n^n,{n,1,Infinity}] （原样输出，无结果） 从上例可见：Sum[ ]是用定义求级数和的，结果为有限实数的，级数收敛于它。但有些级数仅用定义无法求出其和，甚至无法判断其收敛性，如上例中的（5）。 这时，借助Limit[ ]函数，运用比较判别法、比值判别法先来判断级数是否收敛，若收敛，再用NSum[ ]函数求出级数和的近似值。 如上例中的（5）： 输入：algebra\symbolic.m a[n_]:= n!/n^n; Limit[a[n+1]/a[n],n-Infinity] NSum[n!/n^n,{n,1,Infinity}] 输出： 比值判别法,1，收敛。 1.87985 二、函数项级数求和 例2、函数项级数求和，并求其收敛域： （1） （2） 输入：algebra\symbolic.m Sum[x^n/n!,{n,0,Infinity}] Sum[(-1)^(n+1)*x^n/n,{n,1,Infinity}] 输出：Ex ln(1+x) 以下求(2)的收敛域： 输入：Clear[f,n] f[x_]:=(-1)^(n+1)*x^n/n ; 定义函数 a[n_]:=(-1)^(n+1)/n; 定义系数项 Limit[a[n]/a[n+1],n-Infinity] Print[“R=”,Abs[%]] Sum[f[1],{n,1,Infinity}] Sum[f[-1],{n,1,Infinity}] 输出： R=1 求出收敛半径R=1 Ln2 在 x=1处收敛于ln2 Infinity 在 x=-1处发散 所以收敛域为（-1，1） 三、幂级数展开 格式： Series[f(x),{x,x0,n}] 求f(x)在x=x0处的n阶泰勒多项式及余项o[x-x0]n. Normal [Series[f(x),{x,x0,n}]] 求f(x)在x=x0处的n阶麦克劳林级数（不含余项） 例3、（1）sinx的10阶麦克劳林级数 （2）ex的6阶麦克劳林级数 输入：Clear[x] Series[Sin[x],{x,0,10}] Normal [%] Normal [Series[Exp[x], {x,0,6}]] 输出：略 *三、富里埃级数 例 将f(x)= 展开为富里埃级数，并作图显示富里埃级数逼近f(x)的图形。 输入： a0=1/Pi*( Integrate[x,{x,0,Pi}]) （*计算a0*） an=1/Pi*(Integrate[x*Cos[n*x],{x,0,Pi}]) （*计算an*） bn=1/Pi*(Integrate[x*Sin[n*x],{x,0,Pi}]) （*计算bn*） f[x_]:=Which[-3Pi=x-2Pi,0,-2Pi=x-Pi,x+2Pi, -Pi=x0,0,0=xPi,x,Pi=x2Pi,0,2Pi=x=3Pi,x-2Pi]; （*分段函数*） For[i=1,i40,i+=5,fu[x_]:=a0/2+Sum[an*Cos[n*x]+bn*Sin[n*x],{n,1,i}]; （*8个不同级富里埃级