数学实验一new.docVIP

毕节学院实验报告 姓名 杨永福 学号 24121101039 实验项目 微积分基础 实验日期 2013年9月1日 实验目的 了解数学软件Mathematica的基本功能,通过实验来进一步学习使用该软件， 能用Mathematica验证或观察得出微积分学的几个基本结论。完成简单的数学实验。 实验环境 Windows环境下 数学软件Mathematica 实验基本理论和方法 画一些简单函数的图像，求函数的极限，求函数的导数，函数逼近法 实验步骤和内容 实验内容 1． 分别取n=5,50,500,在同一坐标系中画出区间[-4П，4П]上函数f(x)=sinx与 的图象。观察当n增加时向sinx逼近的现象 2. 对用Nsum语句计算“大和”和“小和”及它们的平均值，观察它们的变化趋势。得出S(2)的近似值。再用求数值积分的语句 NIntegrate[1/x,{x,1,2}] 求S(2)，将所得的结果与前面用定义求出的结果相比较。 (2)用语句 S[x_]:= NIntegrate[1/t,{t,1,x}] 定义函数。再用画图语句 Plot[S[x],{x,0.1,10}] 画出这个函数在区间[0。1，10]上的图像（也可以换成别的区间）。观察图像的形状，它好像是什么函数？ 观察出(2)中画出的图像可能是对数函数的图像。现在来求这个对数的底 ，它应满足条件。从图像上看比3稍小一些。以3作为初始值，用练习2所介绍的FindRoot语句求所满足的方程的根的较精确的值。也可以用练习2介绍的牛顿切线法，利用递推关系式 由根的近似求出更好的近似值。这里。迭代语句为 g[a_]:=a-(S([a]-1)a NestList[g,3,4] 计算结果发现恰是自然对数的底e。在同一坐标系中作出函数S(x)和的图像，发现二者完全重合。 实验步骤： 1. 分别取n=5,50,500,在同一坐标系中画出区间[-4П，4П]上函数f(x)=sinx与 的图象。观察当n增加时向sinx逼近的现象 语句为： fgsin=Plot[Sin[x],{x,-4pi,4pi}, PlotStyle-{RGBColor[1,0,0]}]; p[x-,n-]:=x*Product[1-x^2/((k*pi)^2),{k,1,n}] fgproduct=Plot[p[x,50],{x,-4pi,4pi}] Show[fgsin,fgproduct] 图像为： 将实验中的p[x,50]换成p[x,500]的实验结果为： fgsin=Plot[Sin[x],{x,-4pi,4pi}, PlotStyle-{RGBColor[1,0,0]}]; p[x-,n-]:=x*Product[1-x^2/((k*pi)^2),{k,1,n}] fgproduct=Plot[p[x,500],{x,-4pi,4pi}] Show[fgsin,fgproduct] 实验结论 实验结论 通过实验观察，我们发现随着n的增加Pn(x)的函数图象逼近于sinx的函数图象，当n=500时函数图象已经非常接近，所以Pn(x)是sinx的良好逼近 2.对用Nsum语句计算“大和”和“小和”及它们的平均值，观察它们的变化趋势。得出S(2)的近似值。再用求数值积分的语句 NIntegrate[1/x,{x,1,2}] 求S(2)，将所得的结果与前面用定义求出的结果相比较。 (2)用语句 S[x_]:= NIntegrate[1/t,{t,1,x}] 定义函数。再用画图语句 Plot[S[x],{x,0.1,10}] 画出这个函数在区间[0。1，10]上的图像（也可以换成别的区间）。观察图像的形状，它好像是什么函数？ (3)观察出(2)中画出的图像可能是对数函数的图像。现在来求这个对数的底 ，它应满足条件。从图像上看比3稍小一些。以3作为初始值，用练习2所介绍的FindRoot语句求所满足的方程的根的较精确的值。也可以用练习2介绍的牛顿切线法，利用递推关系式 由根的近似求出更好的近似值。这里。迭代语句为 g[a_]:=a-(S([a]-1)a NestList[g,3,4]