第三章 基础实验
基础实验方法—基础实验示范：函数与简单函数表示
第一部分 实验指导书
一、实验目的
1．理解Taylor公式的意义；
2．认识Taylor公式的地位和作用； 3．了解较复杂函数的简单函数表示。
二、实验使用的软件
Mathematica 5.0或以上版本.
三、实验的基本理论及方法
1．Taylor公式
1．1带皮亚诺余项的Taylor公式
设函数f(x)在x0处n阶可导, 则
f(x)?
?k?0nf(k)(x0)(x?x0)k??((x?x0)n). k!特别地x0?0，即得Maclaurin公式
f(x)?
?k?0nf(k)(0)kx??(xn). k!1．2带拉格朗日余项的Taylor公式
设函数f(x)?C[a, b],且f(x)?C(a, b),x,x0?[a,b], 则
(n)(n?1)f(x)??k?0nf(k)(x0)f(n?1)(?)k(x?x0)?(x?x0)n?1 k!(n?1)!其中?介于x与x0之间.
特别地x0?0，即得Maclaurin公式
f(x)??k?0nf(k)(0)kf(n?1)(?)n?1x?x k!(n?1)!其中?介于x与0之间.
2．幂级数展开
给定函数f(x)及任意一点x0是否能找到一个幂级数
?an?0?n(x?x0)，在其收敛区间内
的和函数恰好就是给定的函数f(x)呢？如果能找到这样的幂级数，我们就说f(x)在x0能展开成幂级数，而该幂级数就称为f(x)的在该点处的幂级数展开式。
3．傅里叶级数展开
对波的研究在物理学和工程技术中显得非常重要，它反映了物质作周期运动的运动规律，我们常常用一个以T为周期的周期函数f(t)?f(t?T)来描述它。而简谐振动是最简
?t??)，其中y表示动点的位置，t表示时单的一种周期运动，其运动规律为y?Asin(间，A表示振幅，?是初相，?为角频率.那么其它的波能否用无穷多个简谐波的叠加来表示是傅里叶级数所要解决的问题。
若函数f(x)是以2?为周期的周期函数，且在区间[??,?]上连续或只有有限个第一类间断点，而且只有有限个极值点（上述条件称为狄里克雷充分条件），则有 （1）当x是f(x)的连续点时，
a0? f(x)???(ancosnx?bnsinnx) (1)
2n?1其中的系数an,bn 由式（2）确定
1??a?f(x)cosnxdx,n?0,1,2,???n???? （2） ???b?1f(x)sinnxdx,n?1,2,?n??????其中，式(1)的右端称为函数f(x)的傅立叶级数；式（2）称为傅立叶系数公式。
（2）当x是f(x)的间断点时，傅立叶级数收敛于
1[f(x?0)?f(x?0)] 2四、实验的内容与步骤
1.编写Mathematica程序，从图象上观察多项式与函数的接近或逼近
2x在同一坐标系里分别作出多项式函数y?x，y?x?，y?x2，y?x?x3，
2x3x3x3x54，y?x?y?x??x，y?x??，…和函数y?sinx的图象．观察这些
3!5!3!3!多项式函数的图象向y?sinx的图象逼近的情况．
函数f(x)、g(x)?sinx在区间[a,b]上图象可用如下Mathematica程序画出
f[x_]=expr ;
g[x_]=Sin[x];
Plot[{f,g}，{x，a，b},{PlotStyle->{RGBColor[1,0,0], RGBColor[0,1,0]}}] 思考：哪些多项式函数能与y?sinx逼近？在什么范围内逼近？其它函数？
2.构造多项式与函数逼近
n设多项式函数pn(x)?a0?a1x???anx与函数f(x)逼近，则
f(k)(0)ak?,k?0,1,?,n
k!Mathematica计算程序如下
n=n0;
f[x_]=expr ;
a[x,k_]=D[f,{x,k}]/k!; Table[a[0,k],{k,0,n}];
p[x_,n]=Sum[a[0,k]*x^k,{k,0,n}]
先对sinx分别构造一阶、二阶、…、十五阶Maclaurin多项式，并从图象观察逼近程度与范围。
取x0?1对sinx分别构造一阶、二阶、…、十五阶Taylor多项式，并从图象观察逼近程度与范围。
再取x0??/3对sinx分别构造一阶、二阶、…、十五阶Taylor多项式，并从图象观察逼近程度与范围。
当n??时，Maclaurin(Taylor)多项式函数趋向于什么函数?
3．傅立叶级数
分别取n?10,20，画出函数y??1sin?2k?1?x在区间??3?,3??上的图象．当2k?1k?1nn??时，这个函数趋向于什么函数?
Mathematica程序是：
输出结果：
0.750.50.25-6-4-2-0.25-0.5-0.7524 图3.1.1
Mathematica程序是：
输出结果：
图3.1.2
Mathematica没有专门的命令将一个周期函数进行傅里叶级数展开，但我们可以通过下列的程序将一个以2?为周期的周期函数展开成有限阶不带任何余项的傅里叶级数
n=Input[“n=”]; f[x_]=
Input[＂f[x]=＂]
L=(1.0/Pi*NIntegrate[f[x],{x,-Pi,Pi}]; For[i=1,i<=n,i++,
L=L+(1.0/Pi)*NIntegrate[f[x]*Cos[i*x],{x,-Pi,Pi}]*Cos[i*x]+(1.0/Pi) *NIntegrate[f[x]*Sin[i*x],{x,-Pi,Pi}]*Sin[i*x];]; L
例 编辑一个程序从图形上演示傅里叶级数逐步逼近锯齿波 f(x)??的过程。
解
n＝Input[“n=”];f[x_]=
Which[x>=-Pi&&x<0,x+Pi,x>=0&&x<=Pi,x];
?x??,???x?0
x,0?x???L=(1.0/Pi)*NIntegrate[f[x],{x,-Pi,Pi}];
For[i=1,i<=n,i++,L=L+(1.0/Pi)*NIntegrate[f[x]*Cos[i*x],{x,-Pi,Pi}]*Cos[i*x]+1.0/P
i*NIntegrate[f[x]*Sin[i*x],{x,-Pi,Pi}]*Sin[i*x];
Plot[L,{x,-2.0*Pi,2.0*Pi},Axes->True];]
图3.1.3
图3.1.4
图3.1.5
图 3.1.3、图 3.1.4及图3.1.5分别是锯齿波的 10阶、20阶傅里叶和式及其本身的图形。从这些图形我们可以观察出，随着n的增大，傅里叶和式的图形越来越接近锯齿波的图形，因此傅里叶和式的图形在n??时的极限状态即为锯齿波。
作业：
1．解读实验指导书。（做几件事？目的？） 2．制订实验计划。（实验程序；实验思路。）