《数学实验报告》样例[精选].doc

数学实验报告 实验序号:1 日期:2010年6月18 日星期五 班级 07303 姓名 王 杨 学号 0730352 实验 名称 无穷级数与函数逼近 问题背景描述： 无穷级数在函数性质研究与工程近似计算中有着重要的理论和应用价值。函数逼近就是用一系列简单函数无限接近给定的函数，它又分为局部和整体逼近。在传统的课堂教学中，这两种概念很难讲清楚。 实验目的： 用Mathematica显示无穷数项级数前项部分和的变化趋势，理解无穷级数的概念。学习Mathematica关于函数展开成幂级数的命令及方法。研究幂级数的收敛域。观察Fourier级数对周期函数的逼近情况。 实验原理与数学模型： 1、无穷数项级数前项部分和称为该级数的部分和，若存在，称该级数收敛且极限值为该级数的和。 2、在收敛域内幂级数收敛于一个和函数，但反之，任给函数是否可以看成某个幂级数的和函数，即函数幂级数的展开问题。 3、Fourier级数经常用于谐波分析，在物理和电工学中有重要应用。 实验所用软件及版本： Mathematica 4.1 主要内容(要点)： 级数的部分和的变化趋势； 将函数展开成泰勒级数； 幂级数收敛域； Fourier级数逼近函数的演示。 实验过程记录(含基本步骤、主要程序清单及异常情况记录)： 级数部分和的变化趋势 在实验一中，我们通过图形能够清楚地显示极限的变化趋势，而级数的和就是部分和序列的极限，下面我们采用散点图来观察级数部分和序列的变化趋势。 例1 观察级数的部分和序列的变化趋势，并求和。 解：（1）可以用以下两种方法，从图形来观察级数的敛散性。 （a）利用“Table”命令生成部分和数列的数据点集后作点图，输入语句如下： s[n_]:=Sum[(-1)^(k-1)/(2k-1),{k,1,n}]; data=Table[s[n],{n,1,400}] ;ListPlot[data] 图1 图2 运行后见图1。从图中可以看到级数收敛，级数和大约为0.786。 （b）将级数的所有部分和用竖直线段画出，得到类似条形码的图形，通过这种图形来看出级数的收敛性。输入命令如下： sn=0;n=1;h={};m=3;While[1/n10^(-m),sn=sn+(-1)^(n-1)/(2n-1);h=Append[h,Graphics[{RGBColor[Abs[Sin[n]],0,1/n],Line[{{sn,0},{sn,1}}]}]];n++]; Show[h,PlotRange-{-0.2,1.3},Axes-True] 运行后见图2。从图中可以看出，级数的和在0.785与0.786之间，并且可以通过改变m的值来提高观察到的和的精度。 （2）求和。 对于这个级数，可以通过基本输入模板利用求和符号来直接求和，只要输入命令Sum[(-1)^(n-1)/(2n-1),{n,1,Infinity}]，运行后得到级数和的精确值：。但并不是所有的级数都能如此求和，但可以利用下面的命令来求得和的近似值： NSum[(-1)^(n-1)/(2n-1),{n,Infinity}] 运行后得结果均为：0.785398。为了得到更高的精度，还可以用“N”命令来求函数“Sum”所得的值。大家可以比较一下下面两条命令输出结果的差异： N[Sum[(-1)^(n-1)/(2n-1),{n,Infinity}],30] N[NSum[(-1)^(n-1)/(2n-1),{n,Infinity}],30] 函数的幂级数展开 如果函数在邻域内具有任意阶导数，则函数可以展开为处的幂级数 ， 称之为泰勒级数。特别地，当时，称为麦克劳林级数。 将函数展开为的幂级数，并利用图形考察幂级数的部分和逼近函 数的情况。 解：根据幂级数的展开公式，若能展开成的幂级数，其展示为，因此首先定义函数，再计算点的阶导数，最后构成和式。不妨设，输入命令如下： m=-1;f[x_]:=(1+x)^m;x0=0; g[n_,x0_]:=D[f[x],{x,n}]/.x-x0; s[n_,x_]:=Sum[g[k,x0]/k!*(x-x0)^k,{k,0,n}]; 命令中的函数s[n_,x_]表示的是函数在处的n阶泰勒级数。下面我们通过以下命令观察幂级数的部分和逼近函数的情况： t=Table[s[n,x],{n,20}]; p1=Plot[Evaluate[t],{x,-1/2,1/2}]; p2=Plot[(1+x)^m,{x,-1/2,1/2},Plot