项目四无穷级数和微分方程.pdfVIP

项目四 无穷级数与微分方程 实验 1 无穷级数（基础实验） 实验目的 观察无穷级数部分和的变化趋势 , 进一步理解级数的审敛法以及幂级数部分和对函数的 逼近 . 掌握用 Mathematica 求无穷级数的和 , 求幂级数的收敛域 , 展开函数为幂级数以及展 开周期函数为傅里叶级数的方法 . 数项级数 例 1.1 ( 教材 例 1.1) 1 (1) 观察级数 的部分和序列的变化趋势 . 2 n 1 n 1 (2) 观察级数 的部分和序列的变化趋势 . n 1 n 输入 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 n 1 1 例 1.2 ( 教材 例 1.2) 画出级数 ( 1) 的部分和分布图 . n n 1 输入命令 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 例 1.3 求 的值 . 2 4n