本实验的目的是用Mathematica显示级数部分和的变化趋.doc

级 数本实验的目的是:用Mathematica显示级数部分和的变化趋势；学会如何利用幂级数的部分和对函数的逼近以及进行函数值的近似计算；展示傅里叶级数对周期函数的逼近情况。级数部分和的变化趋势在实验一中，我们通过图形能够清楚地显示极限的变化趋势，而级数的和就是部分和序列的极限，下面我们采用散点图来观察级数部分和序列的变化趋势。例1 观察级数的部分和序列的变化趋势，并求和。解：（1）可以用以下两种方法，从图形来观察级数的敛散性。（a）利用“Table”命令生成部分和数列的数据点集后作点图，输入语句如下： 图1 图2运行后见图1。从图中可以看到级数收敛，级数和大约为0.786。（b）将级数的所有部分和用竖直线段画出，得到类似条形码的图形，通过这种图形来看出级数的收敛性。输入命令如下： 运行后见图2。从图中可以看出，级数的和在0.785与0.786之间，并且可以通过改变m的值来提高观察到的和的精度。（2）求和。对于这个级数，可以通过基本输入模板利用求和符号来直接求和，只要输入命令，运行后得到级数和的精确值：。但并不是所有的级数都能如此求和，但可以利用下面的命令来求得和的近似值： 运行后得结果均为：0.785398。为了得到更高的精度，还可以用“N”命令来求函数“Sum”所得的值。大家可以比较一下下面两条命令输出结果的差异：函数的幂级数展开 如果函数在邻域内具有任意阶导数，则函数可以展开为处的幂级数 ，称之为泰勒级数。特别地，当时，称为麦克劳林级数。将函数展开为的幂级数，并利用图形考察幂级数的部分和逼近函数的情况。解：根据幂级数的展开公式，若能展开成的幂级数，其展示为，因此首先定义函数，再计算点的阶导数，最后构成和式。不妨设，输入命令如下：命令中的函数s[n_,x_]表示的是函数在处的n阶泰勒级数。下面我们通过以下命令观察幂级数的部分和逼近函数的情况：运行结果如图3。从图中可以看到，当n越大时，幂级数越逼近函数。 图33、傅里叶级数设是以2T为周期的周期函数，在任一周期内，除在有限个第一类间断点外都连续，并且只有有限个极值点，则可以展开为傅里叶级数：，其中，且傅里叶级数在任一点处收敛于。下面通过例子来研究傅里叶级数的生成以及级数部分和逼近函数的情况。画图观察周期为、振幅为1的方波函数展成的傅里叶级数的部分和逼近的情况。解：根据傅里叶系数公式可得：，，，故输入以下命令，从输出的图形（图4）观察傅里叶级数的部分和逼近的情况： 图4从图表可以看出，n越大逼近函数的效果越好，还可以注意到傅里叶级数的逼近是整体性的。将函数展开成周期为3的傅里叶级数。解：在下面的Mathematica程序中，先定义了一个产生傅里叶级数的前k项部分和的函数，在后面的“For”循环中连续调用该函数，输出了的傅里叶级数的前项部分和函数，并画出了图形及同一坐标下与其傅里叶级数前7项和函数的图形（如图5）。 图5实验习题观察级数的部分和序列的变化趋势，并求和改变例2中m及的数值来求函数的幂级数及观察其幂级数的逼近函数的情况。观察函数展成的傅里叶级数的部分和逼近的情况。