项目四无穷级数跟微分方程.docVIP

PAGE PAGE 138 项目四 无穷级数与微分方程 实验1 无穷级数（基础实验） 实验目的 观察无穷级数部分和的变化趋势,进一步理解级数的审敛法以及幂级数部分和对函数的 逼近. 掌握用Mathematica求无穷级数的和, 求幂级数的收敛域, 展开函数为幂级数以及展 开周期函数为傅里叶级数的方法. 数项级数 例1.1 (教材 例1.1) (1) 观察级数的部分和序列的变化趋势. (2) 观察级数的部分和序列的变化趋势. 输入 s[n_]=Sum[1/k^2,{k,n}];data=Table[s[n],{n,100}]; ListPlot[data]; N[Sum[1/k^2,{k,Infinity}]] N[Sum[1/k^2,{k,Infinity}],40] 则输出(1)中级数部分和的变化趋势图1.3. 图1.3 级数的近似值为1.64493. 输入 s[n_]=Sum[1/k,{k,n}];data=Table[s[n],{n,50}]; ListPlot[data,PlotStyle-PointSize[0.02]]; 则输出(2)中级数部分和的的变化趋势图1.4. 图1.4 例1.2 (教材 例1.2) 画出级数的部分和分布图. 输入命令 Clear[sn,g];sn=0;n=1;g={};m=3; While[1/n10^-m,sn=sn+(-1)^(n-1)/n; g=Append[g,Graphics[{RGBColor[Abs[Sin[n]],0,1/n], Line[{{sn,0},{sn,1}}]}]];n++]; Show[g,PlotRange-{-0.2,1.3},Axes-True]; 则输出所给级数部分和的图形（图1.5）,从图中可观察到它收敛于0.693附近的一个数. 图1.5 例1.3 求的值. 输入 Sum[x^(3k),{k,1,Infinity}] 得到和函数 例1.4 (教材 例1.3) 设 求. 输入 Clear[a]; a[n_]=10^n/(n!); vals=Table[a[n],{n,1,25}]; ListPlot[vals,PlotStyle-PointSize[0.012]] 则输出的散点图（1.6）,从图中可观察的变化趋势. 输入 Sum[a[n],{n,l,Infinity}] 则输出所求级数的和. 图1.6 求幂级数的收敛域 例1.5 (教材 例1.4) 求的收敛域与和函数. 输入 Clear[a]; a[n_]=4^(2n)*(x-3)^n/(n+1); stepone=a[n+1]/a[n]//Simplify 则输出 再输入 steptwo=Limit[stepone,n-Infinity] 则输出 这里对a[n+1]和a[n]都没有加绝对值. 因此上式的绝对值小于1时, 幂级数收敛; 大于1 时发散. 为了求出收敛区间的端点, 输入 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}] 则输出 函数的幂级数展开 例1.6 (教材 例1.5) 求的6阶麦克劳林展开式. 输入 Series[Cos[x],{x,0,6}] 则输出 注:这是带皮亚诺余项的麦克劳林展开式. 例1.6 (教材 例1.6) 求在处的6阶泰勒展开式. 输入 Series[Log[x],{x,1,6}] 则输出 例1.7 (教材 例1.7) 求的5阶泰勒展开式. 输入 serl=Series[ArcTan[x],{x,0,5}]; Poly=Normal[serl] 则输出的近似多项式 通过作图把和它的近似多项式进行比较. 输入 Plot[Evaluate[{ArcTan[x],Poly}],{x,-3/2,3/2}, Pl