mathematica 数学实验报告材料 实验一

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数学实验报告

数学与统计学院

信息与计算科学(1)班

郝玉霞

201171020107

数学实验一

一、实验名:微积分基础

二、实验目的:学习使用Mathematica的一些基本功能来验证或观察得出微积分学的几个基本理论。

三、实验环境:学校机房,工具:计算机,软件:Mathematica。

四、实验的基本理论和方法:利用Mathematica作图来验证高中数学知识与大学数学容。

五、实验的容和步骤及结果

容一、验证定积分

dt

t

s

x

⎰=

1

1

与自然对数

x

b ln=

是相等的。

步骤1、作积分

dt

t

s

x

⎰=

1

1

的图象;

语句:S[x_]:=NIntegrate[1/t,{t,1,x}]

Plot[S[x],{x,0.1,10}]

实验结果如下:

图1

dt

t

s

x

⎰=

1

1

的图象

步骤2、作自然对数

x

b ln=

的图象

语句:Plot[Log[x],{x,0.1,10}] 实验结果如下:

2

1

图2

x

b ln=

的图象

步骤3、在同一坐标系下作以上两函数的图象

语句:Plot[{Log[x],S[x]},{x,0.1,10}] 实验结果如下:

2

1

图3

dt

t

s

x

⎰=

1

1

x

b ln=

的图象

容二、观察级数与无穷乘积的一些基本规律。

(1)在同一坐标系里作出函数和它的Taylor展开式的前几项构成的多项式函数,,的图象,观察这些多项式函数的图象向的图像逼近的情况。

语句1:

s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]

Plot[{Sin[x],s[x,2]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}]

实验结果如下:

642

4

2

图4和它的二阶Taylor展开式的图象

语句2:

s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]

Plot[{Sin[x],s[x,3]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,1,1]}]

实验结果如下:

642

1

2

3

图5和它的三阶Taylor展开式的图象语句3:

s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]

Plot[{Sin[x],s[x,4]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,1,0]}]

实验结果如下:

642

1

2

3

图6和它的四阶Taylor展开式的图象语句4:

s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]

Plot[{Sin[x],s[x,5]},{x,-2Pi,2Pi},PlotStyle->{RGB[1,0,0]}]

实验结果如下:

642

32

1

图7和它的五阶Taylor 展开式的图象

语句5:

s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]

Plot[{Sin[x],s[x,2],s[x,3],s[x,4],s[x,5] },{x,-2Pi,2Pi}]

实验结果如

下:

图8 和它

的二、三、四、

五阶

Taylor 展开式的图象 (2)分别取n=10,20,100,画出函数在区间[-3π,3π]上的图像,当n →∞时,这个函数趋向于什么函数?

语句1:

f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}] Plot[f[x,10],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}]

实验结果如下: 64242

642

0.5

图9 n=10时,的图像语句2:

f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]

Plot[f[x,20],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}]

实验结果如下:

642

0.5

图10 n=20时,的图像语句3:

f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]

Plot[f[x,100],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}]

实验结果如下:

相关文档
最新文档