实验1无穷级数基础实验.docx

1、实验1无穷级数基础实验项目四 无穷级数与微分方程实验1 无穷级数（基础实验）实验目的观察无穷级数部分和的变化趋势,进一步理解级数的审敛法以及幂级数部分和对函数的逼近. 掌握用Mathematica求无穷级数的和, 求幂级数的收敛域, 展开函数为幂级数以及展开周期函数为傅里叶级数的方法.基本命令1. 求无穷和的命令Sum 该命令可用来求无穷和. 例如,输入 Sum1/n2,n,l,Infinity则输出无穷级数的和为 命令Sum与数学中的求和号相当. 2. 将函数展开为幂级数的命令Series 该命令的基本格式为 Seriesfx,x,x0,n它将展开成关于的幂级数. 幂级数的最高次幂为余项用表

2、示. 例如,输入 Seriesyx,x,0,5则输出带皮亚诺余项的麦克劳林级数3. 去掉余项的命令Normal在将展开成幂级数后, 有时为了近似计算或作图, 需要把余项去掉. 只要使用Normal命令. 例如,输入 SeriesExpx,x,0,6 Normal%则输出 4. 强制求值的命令Evaluate如果函数是用Normal命令定义的, 则当对它进行作图或数值计算时, 可能会出现问题. 例如,输入fx=NormalSeriesExpx,x,0,3Plotfx,x,-3,3则只能输出去掉余项后的展开式而得不到函数的图形. 这时要使用强制求值命令Evaluate, 改成输入 PlotEval

3、uatefx,x,-3,3则输出上述函数的图形.5. 作散点图的命令ListPlot ListPlot 为平面内作散点图的命令, 其对象是数集,例如,输入ListPlotTablej2,j,16,PlotStyle-PointSize0,012则输出坐标为的散点图(图1.1).图1.16. 符号“/；”用于定义某种规则,“/；”后面是条件. 例如,输入Clearg,gf;gx_:=x/;0=x1gx_:=-x/;-1=x=1则得到分段的周期函数再输入 gf=Plotgx,x,-1,6则输出函数的图形1.2.图1.2注:用Which命令也可以定义分段函数, 从这个例子中看到用“(表达式)/; (

4、条件)”来定义周期性分段函数更方便些. 用Plot命令可以作出分段函数的图形, 但用Mathematica命令求分段函数的导数或积分时往往会有问题. 用Which定义的分段函数可以求导但不能积分. Mathematica内部函数中有一些也是分段函数. 如:Modx,1,Absx,Floorx和UnitStepx.其中只有单位阶跃函数UnitStepx可以用Mathematica命令来求导和求定积分. 因此在求分段函数的傅里叶系数时, 对分段函数的积分往往要分区来积. 在被积函数可以用单位阶跃函数UnitStep的四则运算和复合运算表达时, 计算傅里叶系数就比较方便了.实验举例 数项级数例1.1

5、(教材例1.1)(1)观察级数的部分和序列的变化趋势.(2) 观察级数的部分和序列的变化趋势.输入sn_=Sum1/k2,k,n;data=Tablesn,n,100;ListPlotdata;NSum1/k2,k,InfinityNSum1/k2,k,Infinity,40则输出(1)中级数部分和的变化趋势图1.3.图1.3级数的近似值为1.64493.输入sn_=Sum1/k,k,n;data=Tablesn,n,50;ListPlotdata,PlotStyle-PointSize0.02;则输出(2)中级数部分和的的变化趋势图1.4.图1.4例1.2 (教材例1.2) 画出级数的部分和

6、分布图.输入命令Clearsn,g;sn=0;n=1;g=;m=3;While1/n10-m,sn=sn+(-1)(n-1)/n;g=Appendg,GraphicsRGBColorAbsSinn,0,1/n,Linesn,0,sn,1;n+;Showg,PlotRange-0.2,1.3,Axes-True;则输出所给级数部分和的图形（图1.5）,从图中可观察到它收敛于0.693附近的一个数.图1.5例1.3 求的值.输入 Sumx(3k),k,1,Infinity得到和函数例1.4 (教材例1.3)设 求. 输入Cleara;an_=10n/(n!);vals=Tablean,n,1,25

7、;ListPlotvals,PlotStyle-PointSize0.012则输出的散点图（1.6）,从图中可观察的变化趋势. 输入 Suman,n,l,Infinity则输出所求级数的和.图1.6求幂级数的收敛域 例1.5 (教材例1.4)求的收敛域与和函数. 输入Cleara;an_=4(2n)*(x-3)n/(n+1);stepone=an+1/an/Simplify则输出再输入 steptwo=Limitstepone,n-Infinity则输出这里对an+1和an都没有加绝对值. 因此上式的绝对值小于1时, 幂级数收敛; 大于1时发散. 为了求出收敛区间的端点, 输入ydd=Solv

8、esteptwo=1,xzdd=Solvesteptwo=-1,x则输出由此可知,当时,级数收敛,当或时,级数发散. 为了判断端点的敛散性, 输入 Simplifyan/.x-(49/16)则输出右端点处幂级数的一般项为因此,在端点处,级数发散. 再输入 Simplifyan/.x-(47/16)则输出左端点处幂级数的一般项为因此,在端点处, 级数收敛. 也可以在收敛域内求得这个级数的和函数. 输入 Sum4(2n)*(x-3)n/(n+1),n,0,Infinity则输出函数的幂级数展开 例1.6 (教材例1.5)求的6阶麦克劳林展开式.输入 SeriesCosx,x,0,6则输出注:这是带

9、皮亚诺余项的麦克劳林展开式.例1.6 (教材例1.6)求在处的6阶泰勒展开式.输入 SeriesLogx,x,1,6则输出例1.7 (教材例1.7) 求的5阶泰勒展开式.输入serl=SeriesArcTanx,x,0,5;Poly=Normalserl则输出的近似多项式 通过作图把和它的近似多项式进行比较. 输入PlotEvaluateArcTanx,Poly,x,-3/2,3/2,PlotStyle-Dashing0.01,GrayLevel0,AspectRatio-l则输出所作图形（图1.7）, 图中虚线为函数,实线为它的近似多项式.图1.7例1.9 求在处的8阶泰勒展开, 并通过作图

10、比较函数和它的近似多项式. 输入Clearf;fx_=Exp-(x-1)2*(x+1)2;poly2=NormalSeriesfx,x,1,8PlotEvaluatefx,poly2,x,-1.75,1.75,PlotRange-2,3/2,PlotStyle-Dashing0.01,GrayLevel0则得到近似多项式和它们的图1.8.图1.8例1.10 求函数在处的阶泰勒展开, 通过作图比较函数和它的近似多项式, 并形成动画进一步观察. 因为所以输入DoPlotSum(-1)j*x(2j+1)/(2j+1)!,j,0,k,Sinx,x,-40,40,PlotStyle-RGBColor1,

11、0,0,RGBColor0,0,1,k,1,45则输出为的3阶和91阶泰勒展开的图形. 选中其中一幅图形,双击后形成动画. 图1.9是最后一幅图.图1.9例1.11 利用幂级数展开式计算(精确到).因为根据在处的展开式有故前项部分和为输入命令sn_=3(1-1/(5*34)-SumProduct5i-1,i,1,k-1/(5k k!3(4k),k,2,n-1);rn_=Product5i-1,i,1,n-1/5n/n!3(4n-5)/80;delta=10(-10);n0=100;DoPrintn=,n,sn=,Nsn,20;Ifrndelta,Break;Ifn=n0,Printfailed

12、,n,n0则输出结果为傅里叶级数例1.12 (教材例1.8) 设是以为周期的周期函数,它在的表达式是将展开成傅里叶级数. 输入Clearg;gx_:=-1/;-Pi=x0gx_:=1/;0=xPigx_:=gx-2Pi/;PiRGBColor0,1,0; 则输出的图形 (图1.10).图1.10因为是奇函数, 所以它的傅里叶展开式中只含正弦项. 输入b2n_:=b2n=2 Integrate1*Sinn*x,x,0,Pi/Pi;fourier2n_,x_:=Sumb2k*Sink*x,k,1,n;tun_:=Plotgx,Evaluatefourier2n,x, x,-Pi,5 Pi,Plot

13、Style-RGBColor0,1,0,RGBColor1,0.3,0.5,DisplayFunction-Identity; (*tun是以n为参数的作图命令*)tu2=Tabletun,n,1,30,5; (*tu2是用Table命令作出的6个图形的集合*)toshow=Partitiontu2,2; (*Partition是对集合tu2作分割, 2为分割的参数*)ShowGraphicsArraytoshow (*GraphicsArray是把图形排列的命令*)则输出6个排列着的图形（图1.11）,每两个图形排成一行. 可以看到越大, 的傅里叶级数的前项和与越接近. 图1.11实验习题1.求下列级数的和:(1) (2) (3) (4) 2. 求幂级数的收敛域与和函数. 3. 求函数的6阶麦克劳林多项式. 4. 求的6阶麦克劳林多项式.5. 设,求的5阶和10阶麦克劳林多项式,把两个近似多项式和函数的图形作在一个坐标系内.6. 设在一个周期内的表达式为, 将它展开为傅里叶级数(取6项), 并作图.7. 设在一个周期内的表达式为, 将它展开为傅里叶级数(取8项), 并作图. 8. 求级数的和的近似值.