《基础数学实验示例实验指导书.docx》由会员分享，可在线阅读，更多相关《基础数学实验示例实验指导书.docx（8页珍藏版）》请在三一办公上搜索。
1、基础数学实验示例实验指导书第三章 基础实验 基础实验方法基础实验示范：函数与简单函数表示 第一部分 实验指导书 一、实验目的 1理解Taylor公式的意义； 2认识Taylor公式的地位和作用； 3了解较复杂函数的简单函数表示。 二、实验使用的软件 Mathematica 5.0或以上版本. 三、实验的基本理论及方法 1Taylor公式 11带皮亚诺余项的Taylor公式 设函数f(x)在x0处n阶可导, 则 f(x)= k=0nf(k)(x0)(x-x0)k+o(x-x0)n). k!特别地x0=0，即得Maclaurin公式 f(x)= k=0nf(k)(0)kx+o(xn). k!12带
2、拉格朗日余项的Taylor公式 设函数f(x)Ca, b,且f(x)C(a, b),x,x0a,b, 则 (n)(n+1)f(x)=k=0nf(k)(x0)f(n+1)(x)k(x-x0)+(x-x0)n+1 k!(n+1)!其中x介于x与x0之间. 特别地x0=0，即得Maclaurin公式 f(x)=k=0nf(k)(0)kf(n+1)(x)n+1x+x k!(n+1)!其中x介于x与0之间. 2幂级数展开 给定函数f(x)及任意一点x0是否能找到一个幂级数an=0n(x-x0)，在其收敛区间内的和函数恰好就是给定的函数f(x)呢？如果能找到这样的幂级数，我们就说f(x)在x0能展开成幂级
3、数，而该幂级数就称为f(x)的在该点处的幂级数展开式。 3傅里叶级数展开 对波的研究在物理学和工程技术中显得非常重要，它反映了物质作周期运动的运动规律，我们常常用一个以T为周期的周期函数f(t)=f(t+T)来描述它。而简谐振动是最简wt+j)，其中y表示动点的位置，t表示时单的一种周期运动，其运动规律为y=Asin(间，A表示振幅，j是初相，w为角频率.那么其它的波能否用无穷多个简谐波的叠加来表示是傅里叶级数所要解决的问题。 若函数f(x)是以2p为周期的周期函数，且在区间-p,p上连续或只有有限个第一类间断点，而且只有有限个极值点，则有 当x是f(x)的连续点时， a0 f(x)=+(an
4、cosnx+bnsinnx) (1) 2n=1其中的系数an,bn 由式确定 1pa=f(x)cosnxdx,n=0,1,2,Lnp-p pb=1f(x)sinnxdx,n=1,2,Lnp-p其中，式(1)的右端称为函数f(x)的傅立叶级数；式称为傅立叶系数公式。 当x是f(x)的间断点时，傅立叶级数收敛于 1f(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的图象观察
5、这些3!5!3!3!多项式函数的图象向y=sinx的图象逼近的情况 函数f(x)、g(x)=sinx在区间a,b上图象可用如下Mathematica程序画出 fx_=expr ; gx_=Sinx; Plotf,g，x，a，b,PlotStyle-RGBColor1,0,0, RGBColor0,1,0 思考：哪些多项式函数能与y=sinx逼近？在什么范围内逼近？其它函数？ 2.构造多项式与函数逼近 n设多项式函数pn(x)=a0+a1x+L+anx与函数f(x)逼近，则 f(k)(0)ak=,k=0,1,L,n k!Mathematica计算程序如下 n=n0; fx_=expr ; ax,
6、k_=Df,x,k/k!; Tablea0,k,k,0,n; px_,n=Suma0,k*xk,k,0,n 先对sinx分别构造一阶、二阶、十五阶Maclaurin多项式，并从图象观察逼近程度与范围。 取x0=1对sinx分别构造一阶、二阶、十五阶Taylor多项式，并从图象观察逼近程度与范围。 再取x0=p/3对sinx分别构造一阶、二阶、十五阶Taylor多项式，并从图象观察逼近程度与范围。 当n时，Maclaurin(Taylor)多项式函数趋向于什么函数? 3傅立叶级数 分别取n=10,20，画出函数y=1sin(2k-1)x在区间-3p,3p上的图象当2k-1k=1nn时，这个函数趋
7、向于什么函数? Mathematica程序是： 输出结果： 0.750.50.25-6-4-2-0.25-0.5-0.7524图3.1.1 Mathematica程序是： 输出结果： 图3.1.2 Mathematica没有专门的命令将一个周期函数进行傅里叶级数展开，但我们可以通过下列的程序将一个以2p为周期的周期函数展开成有限阶不带任何余项的傅里叶级数 n=Input“n=”; fx_= Inputfx= L=(1.0/Pi*NIntegratefx,x,-Pi,Pi; Fori=1,i=-Pi&x=0&x=Pi,x; x+p,-px0x,0xpL=(1.0/Pi)*NIntegratefx,x,-Pi,Pi; Fori=1,iTrue; 图3.1.3 图3.1.4 图3.1.5 图 3.1.3、图 3.1.4及图3.1.5分别是锯齿波的 10阶、20阶傅里叶和式及其本身的图形。从这些图形我们可以观察出，随着n的增大，傅里叶和式的图形越来越接近锯齿波的图形，因此傅里叶和式的图形在n时的极限状态即为锯齿波。 作业： 1解读实验指导书。 2制订实验计划。
