数学实验 Mathematic实验九 无穷级数.doc

天水师范学院数学与统计学院 实验报告 实验项目名称 无穷级数 所属课程名称 数学实验 实 验 类 型 微积分实验 实 验 日 期 2011.11.16 班 级 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 1.掌握用Mathematica求无穷级数的和; 2.用Mathematica求幂级数的收敛域; 3.用Mathematica展开函数为幂级数以及展开周期函数为傅里叶级数的方法. 【实验原理】 1.求无穷和命令Sum. Sum命令也可以用来求无穷和. Sum[1/n^2，{n，1，Infinity}]Series[y[x]，{x，0，5}]展开成幂级数后，有时为了近似计算或作图，需要把余项去掉，只要使用Normal命令，例如输入 Series[Exp[x]，{x，0，6}] Normal[%]fx=Normal[Series[Exp[x]，{x，0，3}]] Plot[fx，{x，-3，3}]Plot[Evaluate[fx]，{x，-3，3}]ListPlot[Table[j^2，{j，16}]，PlotStyle?PointSize[0.012]]Clear[g，gf]； g[x_]:=x/；0?x1 g[x_]:=-x/；-1?x0 g[x_]:=g[x-2]/；x?1 gf=Plot[g[x]，{x，-1，6}]…(表达式)/；…(条件)”来定义周期性分段函数更方便些.用Plot命令可以作出分段函数的图形，但用Mathematica命令求分段函数的导数或积分时往往会有问题.用which定义的分段函数可以求导，但不能积分.Mathematica内部函数中有一些也是分段函数，如：Mod[x，1]，Abs[x]，Floor[x]和Unitstep[x].其中只有单位阶跃函数Uniltstep[x]可以用Mathematica命令来求导和求定积分,因此在求分段函数的傅里叶系数时，对分段函数的积分往往要分区间来积,在被积函数可以用单位阶跃函数Unitstep的四则运算和复合运算表达时，计算傅里叶系数就比较方便了. 【实验环境】 Mathematic 4 二、实验内容： 【实验方案】 1.数项级数; 2.求幂级数的收敛域; 3.函数的幂级数展开; 4.傅里叶级数. 【实验过程】（实验步骤、记录、数据、分析） 1.数项级数 Sum[1/(4n^2+8n+3)，{n，1，Infinity}]的和. 例9.2　求级数的和. Sum[x^(3k)，{k，1，Infinity}]，求. Clear[a]； a[n_]=10^n/(n!)； vals=Table[a[n]，{n，1，25}]； ListPlot[vals，PlotStyle?PointSize[0.012]]Sum[a[n]，{n，1，Infinity}]收敛域与和函数. Clear[a]； a[n_]=4^(2n)*(x-3)^n/(n+1)； stepone=a[n+1]/a[n]//Simplify steptwo=Limit[stepone，n?Infinity]ydd=Solve[steptwo?1，x] zdd=Solve[steptwo?-1，x]Simplify[a[n]/.x?(49/16)] Simplify[a[n]/.x?(47/16)] Sum[4^(2n)*(x-3)^n/(n+1)，{n，0，Infinity}]的6阶麦克劳林展开式. Series[Cos[x]，{x，0，6}]在处的6阶泰勒展开式. Series[Log[x]，{x，1，6}]的5阶麦克劳林展开式. ser1=Series[ArcTan[x]，{x，0，5}]； poly=Normal[ser1]Plot[Evaluate[{ArcTan[x]，poly}]，{x，-3/2，3/2}，PlotStyle?{Dashing[{0.01}]，GrayLevel[0]}，AspectRatio?1]在处的8阶泰勒展开，并通过作图比较函数和它的近似多项式. Clear[f]； f[x_]=Exp[-(x-1)^2*(x+1)^2]； poly2=Normal[Series[f[x]，{x，1，8}]] Plot[Evaluate[{f[x]，poly2}]，{x，-1.5，1.5}，PlotRange?{-2，3/2}，Plo