mathematica 数学实验报告 实验一

数学实验报告

实

验

一

数学与统计学院

信息与计算科学(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}]

实验结果如下：

2 1

图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）在同一坐标系里作出函数

x

y sin

=

和它的Taylor展开式的前几项构成的

多项式函数

3

!3

x

x

y-

=

，!5

!3

5

3x

x

x

y+

-

=

，???的图象，观察这些多项式函数的图

象向

x

y sin

=

的图像逼近的情况。

语句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

x

y sin

=

和它的二阶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

3

21图5

x y sin =和它的三阶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

3

21图6

x y sin =和它的四阶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

3

21图7

x y sin =和它的五阶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

x y sin =和

它的二、三、四、五阶Taylor 展开式的图象

（2）分别取n=10，20，100，画出函数x

k k y n

k )12sin(121

1--=∑=在区间[-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]}] 实验结果如下：

6424

2

642

0.5

图9 n=10时，

x

k

k

y

n

k

)1

2

sin(

1

2

1

1

-

-

=∑

=的图像

语句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时，

x

k

k

y

n

k

)1

2

sin(

1

2

1

1

-

-

=∑

=的图像

语句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]}] 实验结果如下：

6420.5

图11 n=100时，

x

k k y n

k )12sin(121

1--=∑

=的图像

（3）分别取5,15,100,，在同一坐标系里作出函数x x f s i n

)(=与∏=-

?=n

k k x x x p 1

2

2

2)

1()(π

在区间[-2π，2π]上的图像。

语句1：

p[x_,n_]:=x*Product[1-x^2/(k^2Pi^2),{k,1,n}] Plot[{Sin[x],p[x,5] },{x,-2Pi,2Pi}] 实验结果如下：

642

1.5

1.00.5图12 n=5时,x x f sin )(=与

∏=-

?=n

k k x x x p 1

2

2

2)

1()(π

的图像

语句2:

p[x_,n_]:=x*Product[1-x^2/(k^2Pi^2),{k,1,n}] Plot[{Sin[x],p[x,15] },{x,-2Pi,2Pi}] 实验结果如下：

642 1.0

0.5

图13 n=15时,x x f sin )(=与

∏=-

?=n

k k x x x p 12

2

2

)

1()(π的图像

语句3：

p[x_,n_]:=x*Product[1-x^2/(k^2Pi^2),{k,1,n}] Plot[{Sin[x],p[x,100] },{x,-2Pi,2Pi}] 实验结果如下：

642 1.0

0.5

图14 n=100时,x x f sin )(=与

∏=-

?=n

k k x x x p 1

2

2

2)

1()(π

的图像

六、实验结果分析

内容一、图1、图2分别作出了定积分

dt

t s x

?=11与自然对数x b ln =的图象，大致看来这两幅图是一样的；由图3在同一坐标系里作出以上两函数的图象，可

以看出这两幅图是完全重合的，由此足以证明：定积分

dt

t s x

?=11

与自然对数x b ln =是相等的,这与之前我们得出的结论是完全一致的。

内容二、（1）图4、5、6、7分别作出函数

x y sin =和它的二、三、四、五阶

Taylor 展开式的图象，图8作出了同一坐标系里函数

x y sin =和它的二、三、

四阶Taylor 展开式的图象，经比较可知，奇数阶的更接近正弦函数；（2）图9、

10、11分别作出n=10,20,100时，函数

x

k k y n

k )12sin(121

1--=∑

=的图像，经观

察可知，当n →∞时，这个函数趋向于分段函数；（3）图12、13、14分别作出n=5,15,100时，在同一坐标系里函数x x f sin )(=与

∏=-

?=n

k k x x x p 12

2

2)

1()(π

在区

间[-2π，2π]上的图像，观察知当n 增加时)(x p 的图像向

)sin(x 的图像逼近，

且两个函数在x=0处的导数相同，在任何有限的区间上，当n →∞时函数)(x p 逼

近

)sin(x 。

Mathematica数学实验——随机变量的概率分布

教师指导实验7 实验名称：随机变量的概率分布 一、问题：求二项分布、几何分布、正态分布在给定区间上的概率。 二、实验目的： 学会使用Mathematica求二项分布、几何分布、正态分布在给定区间上的概率及期望和方差。 三、预备知识：本实验所用的Mathematica命令提示 1、BinomialDistribution[n,p] 二项分布； GeometricDistribution[p] 几何分布； NormalDistribution[μ,σ] 正态分布； 2、Domain[dist] 求分布dist的定义域； PDF[dist,x] 求点x处的分布dist的密度值； CDF[dist,x] 求点x处的分布dist的函数值； Mean[dist] 求分布dist的期望；Quantile[dist,x] 求x，使CDF[dist,x]=q Variance[dist] 求分布dist的方差；StandardVariance[dist] 求分布dist的标准差； 四、实验的内容和要求： 1、取50个在1到20的随机整数，求这组数的极差、中位数、均值、方差及标准差； 2、对以上数据绘制样本频率分布直方图； 3、data1={1, 3, 4, 5, 3.5, 3}, data2={3, 2, 5, 3}，在同一图表中绘制data1和data2的条形图，并作一定的修饰。 五、操作提示 1、取50个在1到20的随机整数，求这组数的极差、中位数、均值、方差及标准差； In[1]:=<

数学实验报告格式

《数学实验》实验报告 （ 2012 年 03 月 30 日） 一、实验问题 1、某公司指派5个员工到5个城市工作（每个城市单独一人），希望使所花费的总电话 费用尽可能少。5个员工两两之间每个月通话的时间表示在下面的矩阵的上三角部分（因 为通话的时间矩阵是对称的，没有必要写出下三角部分），5个城市两两之间通话费率表示在 下面的矩阵的下三角部分（同样道理，因为通话的费率矩阵是对称的，没有必要写出上三角 部分）. 试求解该二次指派问题。 通话时间d=[0 1 1 2 3 1 0 2 1 2 1 2 0 1 2 2 1 1 0 1 3 2 2 1 0 ] 城市间通话费率 c=[0 5 2 4 1 5 0 3 0 2 2 3 0 0 0 4 0 0 0 5 1 2 0 5 0] 2、某校毕业生必须至少修：两门数学课、三门运筹学课、两门计算机课。 1）某学生希 望所修课程最少。 2）某学生希望课程少学分多。 3）某学生觉得学分数和课程数这两大目标大致应该三七开。 3、某储蓄所营业时间为上午9:00--下午5:00，储蓄所可以雇佣两类服务员：全职：每 天100元中午12:00--下午2:00之间必须安排1小时的午餐时间半职：每人40 元必须连 续工作4小时 1）储蓄所每天雇佣的半职服务员不超过3人，为使花费最少该如何雇佣两类服务员。 2） 如果不能雇佣半时服务员，花费多少？ 3）如果雇佣半时服务员没有人数限制花费多少？ 二、问题的分析（涉及的理论知识、数学建模与求解的方法等） ?1 1、用xik?? ?0 i人去了k城市 ?1 （i=1...5） xjh?? i人不去k城市?0 j人去了h城市j人没去h城市 (i=1...5) dij表示i和j的通话时间；ckh表示城市k和h之间的费率，数学模型： 5555 min ????c kh dijxikxjh i?1 j?1k?1h?1 ?5 ??xik?1k?1...5?i?1? 5?1i?1 (5) s.t.??xik?k?1 5

离散数学实验报告

《离散数学》实验报告专业网络工程 班级 姓名 学号 授课教师 二 O 一六年十二月

目录 实验一联结词的运算 实验二根据矩阵的乘法求复合关系 实验三利用warshall算法求关系的传递闭包实验四图的可达矩阵实现

实验一联结词的运算 一.实验目的 通过上机实验操作,将命题连接词运算融入到C语言的程序编写中,一方面加强对命题连接词运算的理解,另一方面通过编程实现命题连接词运算,帮助学生复习与锻炼C语言知识,将理论知识与实际操作结合,让学生更加容易理解与记忆命题连接词运算。 二.实验原理 (1) 非运算, 符号:? ,当P=T时 ,?P为F, 当P=F时 ,?P为T 。 (2) 合取, 符号: ∧ , 当且仅当P与Q的真值同为真,命题P∧Q的真值才为真;否则,P∧Q的真值为假。 (3) 析取, 符号: ∨ , 当且仅当P与Q的真值同为假,命题P∨Q的真值才为假;否则,P∨Q的真值为真。 (4) 异或, 符号: ▽ , 当且仅当P与Q的真值不同时,命题P▽Q的真值才为真;否则,P▽Q的真值为真。 (5) 蕴涵, 符号: →, 当且仅当P为T,Q为F时,命题P→Q的真值才为假;否则,P→Q 的真值为真。 (6) 等价, 符号: ? , 当且仅当P,Q的真值不同时,命题P?Q的真值才为假;否 则,P→Q的真值为真。 三.实验内容 编写一个程序实现非运算、合取运算、析取运算、异或运算、蕴涵运算、等价运算。四．算法程序 #include void main() { printf("请输入P、Q的真值\n"); int a,b; scanf("%d%d",&a,&b); int c,d; if(a==1) c=0; else c=1; if(b==1) d=0; else d=1; printf("非P、Q的结果为%d,%d\n",c,d);

离散数学实验报告

离散数学实验报告（实验ABC） 专业班级 学生姓名 学生学号 指导老师 完成时间

目录 第一章实验概述..................................... 错误!未定义书签。 实验目的....................................... 错误!未定义书签。 实验内容....................................... 错误!未定义书签。 实验环境....................................... 错误!未定义书签。第二章实验原理和实现过程........................... 错误!未定义书签。 实验原理....................................... 错误!未定义书签。 建立图的邻接矩阵，判断图是否连通 ............ 错误!未定义书签。 计算任意两个结点间的距离 ................... 错误!未定义书签。 对不连通的图输出其各个连通支 ................ 错误!未定义书签。 实验过程（算法描述）........................... 错误!未定义书签。 程序整体思路 ............................... 错误!未定义书签。 具体算法流程 ................................ 错误!未定义书签。第三章实验数据及结果分析........................... 错误!未定义书签。 建立图的邻接矩阵并判断图是否连通的功能测试及结果分析错误!未定义书签。 输入无向图的边 .............................. 错误!未定义书签。 建立图的连接矩阵 ............................ 错误!未定义书签。 其他功能的功能测试和结果分析................... 错误!未定义书签。 计算节点间的距离 ............................ 错误!未定义书签。 判断图的连通性 .............................. 错误!未定义书签。 输出图的连通支 .............................. 错误!未定义书签。 退出系统 .................................... 错误!未定义书签。第四章实验收获和心得体会........................... 错误!未定义书签。

mathematica数学实验报告

高等数学实验报告 实验一 一、实验题目 1：作出各种标准二次曲面的图形 ParametricPlot3D Sin u Sin v,Sin u Cos v,Cos u ,u,0,Pi ,v,0,2Pi,P Graphics3D ParametricPlot3D u Sin v,u Cos v,u^2,u,0,2,v,0,2Pi,PlotPoints30

Graphics3D ParametricPlot3D u,v,u^2v^2,u,2,2,v,2,2,PlotPoints30 Graphics3D ParametricPlot3D Sec u Sin v,Sec u Cos v,Tan u,u,Pi4,Pi4,v,0,2

Graphics3D t1ParametricPlot3D u^21Sin v,u^21Cos v,u,u,1,5,v,0,2Pi t2ParametricPlot3D u^21Sin v,u^21Cos v,u,u,5,1,v,0,2 show t1,t2 Graphics3D

Graphics3D show Graphics3D,Graphics3D ParametricPlot3D u Cos v,u Sin v,u,u,6,6,v,0,2Pi,PlotPoints60 Graphics3D 2：作出曲面所围的图形 t1ParametricPlot3D Sin u Sin v,Sin u Cos v,Cos u, u,Pi2,pi2,v,0,2Pi,PlotPoints60 t2ParametricPlot3D0.5Cos u12,0.5Sin u, u,0,2Pi,v,0,2Pi,PlotPoints60 t3Plot3D0,PlotPoints60 show t1,t2,t3

数学实验报告

《数学实验》报告 题目：根据数值积分计算方法计 算山东省面积 学生姓名： 学号： 专业班级：机械工程17-1班

2019年4月15日

一、问题背景与提出 图1是从百度地图中截取的山东省地图，试根据前面数值积分计 算方法，计算山东省面积。 图 1 二、实验目的 1、 学会运用matlab 解决一些简单的数学应用问题。 2、 学会运用matlab 建立数学模型。 3、 学会运用一些常见的数值积分计算方法结算实际问题，并 了解其实际意义，建立积分模型。 三、实验原理与数学模型 将积分区间 [a , b] n 等分，每个区间宽度均为h = (b - a) / n , h 称 为积分步长。记 a = x 0 < x 1 < … < x k … < x n = b , 在小区间上用小矩形面积近似小曲边梯形的面积，若分别取左端点和右端点的函数值为小矩形的高，则分别得到两个曲边梯形的面积的近似公式： Ln = h ∑f (x k )n=1k=0 , h = b?a ?

R n =?∑f (x k )n k=1 , h = b?a ? 如果将二者求平均值，则每个小区间上的小矩形变为小梯形，整 个区间上的值变为： Tn =?∑f (X k )n=1 k=1+?2[f (x 0)+f (x n )] 将山东省边界上的点反映在坐标化，运用梯形公式积分计算得山 东省的面积。 四、实验内容（要点） 1、将山东省的地图区域在matlab 中画出 。 2、在坐标系上运用积分方法将所求区域的面积求出。 3、通过比例尺将山东省的实际面积求出。 五、实验过程记录（含基本步骤、主要程序清单及异常情况记录等） 1、 在百度地图中标识出山东省的区域范围，标明对应的比例： 图 2 2、 取出所截取图片中山东的边界的坐标，即将边界坐标化： （1） 运用imread 函数和imshow 函数导入山东省的区域 图片。

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

数学实验报告 实 验 一 数学与统计学院 信息与计算科学(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展开式的图象

小学数学实验报告

竭诚为您提供优质文档/双击可除 小学数学实验报告 篇一：小学数学课题实验总结报告 《实施合作学习，发挥优势互补的研究》 课题实验总结 在上级主管部门和学校领导关心支持下我们开展了《实施合作学习，发挥优势互补》的课题研究。在课题组全体老师两年的不懈努力下，已基本完成本课题研究任务，并取得预期成果。 开展课题实验以来，我们坚持在实践中探索，在探索中实践，取得了初步的成效，主要体现在实验促进了三个方面的转变，一个方面的提高。 一、促进教师教学观念的转变。 参加课题实验后，实验组的老师们通过边实验边学习，不断总结与反思，提升了自己的科研水平，并树立了以“教学是为了促进学生发展”为最终目标的新型教育教学观念。课堂上，老师与学生建立了和谐融洽的师生关系，在精心创设的良好的教学氛围中鼓励学生独立思考、大胆质疑、敢于

探索、勇于创新。让学生在自主、合作、探究的学习过程中，激发学习热情，养成学习习惯，提高学习能力，从而促进了学生的发展。 二、促进学生学习方式的转变。 学生正在由被动学习逐步向主动学习转变，由老师教转变为我能学，由师生间的单向性活动转变为双向性互动、多边性互动，增大了课堂信息量，学生积极主动学习，小组合作、乐于探究，他们发扬团队精神，团队之间互相竞争、优势互补，并培养学生动手、动脑、动口的能力，培养创新意识。课前，学生能积极主动地预习信息窗内容，提出问题并尝试解决。课堂上，学生能够热烈地交流预习所得，积极主动地参与课堂讨论，参与面广，讨论热烈而且有序。课后，能自觉温习知识，深化学习，拓展延伸，并加以运用。绝大部分学生善于表达，敢于提出自己的不同见解，有较强的探究精神，能够提出问题积极思考，并能够多角度思维寻找解决问题的策略，并且培养了学生良好的合作学习的习惯。 学习方式的转变促进了学生全面发展，他们乐学，善学，学有所成。随着学生自主合作探究能力的不断提高，自主性合作性探究性已多个学习层面辐射，辐射到其它学科、班级管理、文体活动等方面。实验班班风好，学风浓，学生对所有科目的学习兴趣盎然、积极主动，全面发展。 三、促进课堂教学格局的转变。

离散数学实验报告--四个实验!!!

《离散数学》 课程设计 学院计算机学院 学生姓名 学号 指导教师 评阅意见 提交日期 2011 年 11 月 25 日

引言 《离散数学》是现代数学的一个重要分支，也是计算机科学与技术，电子信息技术，生物技术等的核心基础课程。它是研究离散量（如整数、有理数、有限字母表等）的数学结构、性质及关系的学问。它一方面充分地描述了计算机科学离散性的特点，为学生进一步学习算法与数据结构、程序设计语言、操作系统、编译原理、电路设计、软件工程与方法学、数据库与信息检索系统、人工智能、网络、计算机图形学等专业课打好数学基础；另一方面，通过学习离散数学课程，学生在获得离散问题建模、离散数学理论、计算机求解方法和技术知识的同时，还可以培养和提高抽象思维能力和严密的逻辑推理能力，为今后爱念族皮及用计算机处理大量的日常事务和科研项目、从事计算机科学和应用打下坚实基础。特别是对于那些从事计算机科学与理论研究的高层次计算机人员来说，离散数学更是必不可少的基础理论工具。 实验一、编程判断一个二元关系的性质（是否具有自反性、反自反性、对称性、反对称性和传递性） 一、前言引语：二元关系是离散数学中重要的内容。因为事物之间总是可以 根据需要确定相应的关系。从数学的角度来看，这类联系就是某个集合中元素之间存在的关系。 二、数学原理：自反、对称、传递关系 设A和B都是已知的集合，R是A到B的一个确定的二元关系，那么集合R 就是A×B的一个合于R={(x,y)∈A×B|xRy}的子集合 设R是集合A上的二元关系： 自反关系：对任意的x∈A,都满足∈R，则称R是自反的，或称R具有自反性，即R在A上是自反的?(?x)((x∈A)→(∈R))=1 对称关系：对任意的x,y∈A，如果∈R，那么∈R，则称关系R是对称的，或称R具有对称性，即R在A上是对称的? (?x)(?y)((x∈A)∧(y∈A)∧(∈R)→(∈R))=1 传递关系：对任意的x,y,z∈A，如果∈R且∈R，那么∈R，则称关系R是传递的，或称R具有传递性，即R在A上是传递的? (?x)(?y)(?z)[(x∈A)∧(y∈A)∧(z∈A)∧((∈R)∧(∈R)→(∈R))]=1 三、实验原理：通过二元关系与关系矩阵的联系，可以引入N维数组，以数 组的运算来实现二元关系的判断。 图示：

MATHEMATICA实验报告

【MATHEMATICA实验报告】 【实验目的】 1.掌握Mathematica软件的启动和退出，以及Mathematica帮助系统。 2.熟悉Mathemaic的计算其功能以及常用的数字函数。 3.掌握变量的定义，变量的操作。 4.掌握函数的定义以及运算。 【实验内容】 1.求下列积分 （1） (4sin()3cos())/(sin()2cos()) x x x x dx ++ ? 输入: y=(4 Sin[x]+3 Cos[x])/(Sin[x]+2Cos[x]); Integrate[y,x] 输出： 2 x-Log[2 Cos[x]+Sin[x]] （2） /2 (cos())^5sin(2) x x dx π ? 输入： y=Cos[x]^5 Sin[2 x] Integrate[y,{x,0,Pi/2}] 输出： Cos x5Sin2x 2 7 （3）1 /(^21)^(3/2) dx x x -+ ? 输入： y=1/(x^2-x+1)^(3/2); Integrate[y,{x,0,1}] 输出： 4 3 2.求积分 1 (1/2)*^(^2/2) e x dx π -∞ - ? 输入：y=E^(-x^2/2)/Sqrt[2*Pi]; NIntegrate[y,{x,Infinity,1}] 输出： -0.158655

3.求y=e^(x^2)在x=0的9阶泰勒公式。 输入： Series[Exp[x^2],{x,0,9}] 输出： 1x 2x 4 2x 66x 824O x 10 4.作出以下参数方程所描述的图形。 （1） 4cos {3sin x t y t ==，（0≤t ≤2π） 输入： ParametricPlot[{4 Cos[t],3 Sin[t]},{t,0,2Pi}] 输出： -4-2 24-3-2 -1 1 2 3 （2）3(cos )^3 {3(sin )^3x t y t -= 输入： ParametricPlot[{3 Cos[t]^3,3 Sin[t]^3},{t,0,2 Pi}] 输出： -3-2-1 123-3-2 -1 12 3

离散数学实验报告()

《离散数学》实验报告 专业网络工程 班级 姓名 学号 授课教师 二 O 一六年十二月

目录 实验一联结词的运算 实验二根据矩阵的乘法求复合关系 实验三利用warshall算法求关系的传递闭包实验四图的可达矩阵实现

实验一联结词的运算 一．实验目的 通过上机实验操作，将命题连接词运算融入到C语言的程序编写中，一方面加强对命题连接词运算的理解，另一方面通过编程实现命题连接词运算，帮助学生复习和锻炼C语言知识，将理论知识与实际操作结合，让学生更加容易理解和记忆命题连接词运算。二．实验原理 (1) 非运算, 符号: ,当P=T时，P为F, 当P=F时，P为T 。 (2) 合取, 符号: ∧ , 当且仅当P和Q的真值同为真，命题P∧Q的真值才为真；否则，P∧Q的真值为假。 (3) 析取, 符号: ∨ , 当且仅当P和Q的真值同为假，命题P∨Q的真值才为假；否则，P∨Q的真值为真。 (4) 异或, 符号: ▽ , 当且仅当P和Q的真值不同时，命题P▽Q的真值才为真；否则，P▽Q的真值为真。 (5) 蕴涵, 符号: →, 当且仅当P为T,Q为F时，命题P→Q的真值才为假；否则，P→Q 的真值为真。 (6) 等价, 符号: ?, 当且仅当P,Q的真值不同时，命题P?Q的真值才为假；否则，P→Q的真值为真。 三．实验内容 编写一个程序实现非运算、合取运算、析取运算、异或运算、蕴涵运算、等价运算。四．算法程序 #include void main() { printf("请输入P、Q的真值\n"); int a,b; scanf("%d%d",&a,&b); int c,d; if(a==1) c=0; else c=1; if(b==1) d=0;

数学应用软件实验报告(mathematica实验程序)1

徐州工程学院数理学院数学应用软件实验报告 课程（实验序号）数学应用软件实验 1 实验地点、日期数学建模机房2011 年 2 月23 日主要仪器设备计算机 使用的软件名称Mathematica 实验类型演示性实验 验证性实验 综合性实验√设计性实验 研究性实验 班级：姓名：孙娅学号：20090402223 一、实验题目名称：函数】变量和表达式 二、实验目的： 理解变量和算式、内核与前端处理器构成的人机对话系统，了解计算的精度问题个Mathematica使用中的几个问题。熟练掌握数的表示和计算、常用数学函数，会绘制简单函数的图形。通过上机初步了解数学应用软件，Mathematica的各种界面。 三、实验内容： 练习题1 1.计算下列各式的数值： (1) Log[2,10] Log[10]/Log[2] (2) Sqrt[Pi^2+1] 1 2 (3) Log[10,3264] Log[3264]/Log[10] (4) E^E ??/2 (5) Cos[135^0] Cos[1] (6) Sin[Pi^2/2] Sin[π2/2] (7) ArcSin[1/2] π/6 (8) 200! 7886578673647905035523632139321850622951359776871732632947425332443594499634033429203042 8401198462390417721213891963883025764279024263710506192662495282993111346285727076331723 7396988943922445621451664240254033291864131227428294853277524242407573903240321257405579

数学实验报告

数学实验报告 实验序号：日期：2016 年

实验过程记录（含基本步骤、主要程序清单及异常情况记录等）: 第一题 选择初速度ｖ=0.6km/s，发射角ａ=４5° X轴方向运动为ｘ＝ｃｏs a＊v*t Y轴方向运动为y=sｉn ａ＊v*t－１／2＊g*t2 统一单位将0。６km/ｓ化为６00m／s 将数据代入利用函数做出运动轨迹，函数式为 8000 6000 4000 2000 5000100001500020000250003000035000 第二题 确定速度为320m/s，求最佳角度使得轨迹与Ｘ轴交点为(10000,0） 先假定发射角为π/4 作图 ＰａrametriｃPlot［{Cｏs［Pi/4]＊320*t，Sin［Pi/4］＊320*ｔ —4．9*t^2}，{t，0,4７｝,AspecｔRatiｏAｕtｏmａｔic] 2500 2000 1500 1000 500 200040006000800010000 进行调整角度调整为π/3.５作图 ParameｔricＰlot[｛Cos［Pi/3.5]＊320＊t，Sｉn[Pｉ/3.5］*320*t－4。9*t^2｝，｛t,0，52},AspeｃtRaｔio Autｏmatic] 3000 2500 2000 1500 1000 500 200040006000800010000

继续进行不断地调整,发现当发射角度为π／3。7时，落点十分接近（１00０0,0)点作图如下 200040006000800010000 500 1000 1500 2000 2500 3000 因此可以确定最合适的发射角就在π/３。7附近,此时可以利用FｉndＲoot函数找出准确值 首先需要对已知式做等量变换： ∵X＝coｓ a＊v*t ∴t=x/（cos a*v) 将上式代入y＝ｓiｎａ＊v*ｔ-1/２＊g*ｔ2 中可得到 Y=tana*x—１/2*ｇ＊（x/(cosa*v)）2 将y=0， x＝10000，g＝９．8, v=320代入利用FｉnｄＲoot函数求解a 的范围在π/３.7附近的a的值: 得出 将这个值由弧度制化为360度制 a=5３.4285° ∴最佳发射角为5３．４285° 第三题 由第二题的320m／ｓ起步进行研究 1.首先研究速度增大运用与第二题相似的研究方法，先大致计算符合要求的角度 （1)V＝３50ｍ/s时，最佳发射角为π/6．8： 200040006000800010000 200 400 600 800 1000 1200 (２）V=40０m/s时，最佳发射角为π／9。5: 0200040006000800010000 200 400 600 800

最佳分数值逼近(mathematica数学实验报告)

姓名 ### 学院 ###### 班级 ######### 学号 ######### 实验题目 最佳分数值逼近 评分 实验目的： 1、用“连分数展开”的方法计算圆周率π的近似值； 2、通过实验来体会“连分数展开”的方法与其他方法的区别，比较各种方法的优劣； 3、尝试用“连分数展开”的方法对其他的数进行展开。 实验环境： 学校机房，Mathematica4.0软件 实验基本理论和方法： 1、Mathematica 中常用的展开数与多项式的函数的使用； 2、计算圆周率π“连分数展开”方法，并且利用特定的函数来展开其他数。 实验内容和步骤： （一）多项式的展开与化简 多项式是表达式的一种特殊的形式，所以多项式的运算与表达式的运算基本一样，表达式中的各种输出形式也可用于多项式的输出。Mathematica 提供一组按不同形式表示代数式的函数。如： 1、 对12 x 1-进行分解，使用的函数为Factor ： 2、 展开多项式 7 x+2（）与5 x+y+7（），使用的函数为Expand:

3、 化简(1)^4(2)^(3)x x x +++与(1)^3(2)^4(3)^(1)x x x x +++-，使用的函数为 Pimplify: 4、 连个多项式相除，总能写成一个多项式和一个有理式相加， Mathematic 中提供两个 函数PolynomialQuotient 和PolynomialRemainder 分别返回商式和余式：

（二）π的连分数展开 π的求解方法之前我们已经有许多种，但都比较繁琐而且误差较大，如何找到误差较小的π的近似值求解方法，我们在所得整数3的基础上进行分析，有了整数3，则 π=3+1x ，其中10.141592653579...x =是3的误差，101x <<。只要能找到1x 的最佳分数逼近值，再加3就得到π的最佳分数近似值。从而我们使用一种方法“连分数展开“，其原理是： 为了寻找与1x 接近的分数，先找与11 1 7.062513305931...A x = =接近的整数，显然 是7.于是111223377 A π=+ ≈+=，这是祖冲之的效率。 在此基础上，我们可以再用上述方法，要找到比 22 7 误差更小的分数近似值，只需要找到比整数7更接近1A 的分数来作为1A 的近似值。由于127A x =+，其中 200.062513305931...1x <=<。先找22 1 15.996594406685...A x = =的最佳整数近似值，显然是16.于是1211113771616A A =+ ≈+=，从而1 2 111355 3331 1113 7716 A A π=+=+≈+ = + +，这就得到祖冲之的密度。 如果还要进一步提高精确度，就应当在考虑2A 的整数近似值16的误差 32160.003405593314...x A =-=，取33 1 293.6345910144...A x = =的整数近似值294，则可

离散数学实验报告格式

《离散数学》实验报告 专业 班级 姓名 学号 授课教师 二 O 一六年十二月

目录 实验一联结词的运算 实验二根据矩阵的乘法求复合关系实验三利用算法求关系的传递闭包实验四图的可达矩阵实现

实验一联结词的运算 一．实验目的 通过上机实验操作，将命题连接词运算融入到C语言的程序编写中，一方面加强对命题连接词运算的理解，另一方面通过编程实现命题连接词运算，帮助学生复习和锻炼C语言知识，将理论知识与实际操作结合，让学生更加容易理解和记忆命题连接词运算。二．实验原理 (1) 非运算, 符号: ,当时，P为F, 当时，P为T 。 (2) 合取, 符号: ∧ , 当且仅当P和Q的真值同为真，命题P∧Q的真值才为真；否则，P∧Q的真值为假。 (3) 析取, 符号: ∨ , 当且仅当P和Q的真值同为假，命题P∨Q的真值才为假；否则，P∨Q的真值为真。 (4) 异或, 符号: ▽ , 当且仅当P和Q的真值不同时，命题P▽Q的真值才为真；否则，P▽Q的真值为真。 (5) 蕴涵, 符号: → , 当且仅当P为为F时，命题P→Q的真值才为假；否则，P→Q 的真值为真。 (6) 等价, 符号: ?, 当且仅当的真值不同时，命题P?Q的真值才为假；否则，P→Q 的真值为真。 三．实验内容 编写一个程序实现非运算、合取运算、析取运算、异或运算、蕴涵运算、等价运算。四．算法程序 <> () { ; ("请选择运算方式\n"); ("1.析取\n"); ("2.合取\n"); ("3.非\n"); ("4.蕴含\n"); ("5.等价\n");

m; (""); ( m>=1 m<=4 ) { ("请输入P Q的值\n"); (" " ); = 1; (m) { 1( ( >= 1)( < 4 ) ) { (0 0) ("P 析取Q = 0\n"); ("P 析取Q = 1\n"); ; (4) ; ("请输入P Q的值\n"); (" " ); } ; 2( ( >= 0)( < 4 ) ) { (1 1) ("P 合取Q = 1\n"); ("P 合取Q = 0\n"); ; (4) ; ("请输入P Q的值\n"); (" " ); } ; 3( ( >= 0)( < 4 ) ) { (0) ("非Q = 1\n"); ("非Q = 0\n");

[vip专享]2013春数学实验基础 实验报告(1)常微分方程

1. 分别用Euler 法和ode45解下列常微分方程并与解析解比较： (1) 30,1)0(,<<=+='x y y x y 编写Euler 法的matlab 函数，命名为euler.m function [t,y]=euler(odefun,tspan,y0,h)t=tspan(1):h:tspan(2);y(1)=y0; for i=1:length(t)-1 y(i+1)=y(i)+h*feval(odefun,t(i),y(i));end t=t';y=y'; 下面比较三者的差别：% ode45 odefun=inline('x+y','x','y');[x1,y1]=ode45(odefun,[0,3],1);plot(x1,y1,'ko');pause hold on ;% Euler·¨ [x2,y2]=euler(odefun,[0,3],1,0.05);plot(x2,y2,'r+');pause hold on ; % 解析解 y0=dsolve('Dy=t+y','y(0)=1');ezplot(y0,[0,3]);pause hold off ; legend('ode45','euler 法','解析解');

Euler 法只有一阶精度，所以实际应用效率比较差，而ode45的效果比较好，很接近真实值。 (2) 2 0.01()2sin(),(0)0,(0)1,05 y y y t y y t ''''-+===<<先写M 文件ex1_2fun.m function f=ex1_2fun(t,y)f(1)=y(2); f(2)=0.01*y(2).^2-2*y(1)+sin(t);f=f(:);% ode45 [t1,y1]=ode45(@ex1_2fun,[0,5],[0;1]);plot(t1,y1(:,1),'ko'); % 解析解 s=dsolve('D2y-0.01*(Dy)^2+2*y=sin(t)','y(0)=0','Dy(0)=1','t') s = [ empty sym ] %由此可知该微分方程无解析解2. 求一通过原点的曲线，它在处的切线斜率等于若上限增为1.58，1.60会(,)x y 2 2,0 1.57.x y x +<

离散数学关系性质的C++或C语言判断实验报告

1.【实验目的】 对称： 通过算法设计并编程实现对给定集合上的关系是否为对称关系的判断，加深学生对关系性质的理解，掌握用矩阵来判断关系性质的方法 自反： 通过算法设计并编程实现对给定集合上的关系是否为自反关系的判断，加深学生对关系性质的理解，掌握用矩阵来判断关系性质的方法。 2.【实验内容】 已知关系R 由关系矩阵M 给出，要求判断由M 表示的这个关系是否为对称关 系。假定R 的关系矩阵为：?????? ? ??=1234210330124321M 3.【实验要求】 C 语言编程实现 4.【算法描述】 对称： 从给定的关系矩阵来判断关系R 是否为对称是很容易的。若M （R 的关系矩阵）为对称矩阵，则R 是对称关系；若M 为反对称矩阵，则R 是反对称关系。因为R 为对称的是等价关系的必要条件，所以，本算法可以作为判等价关系算法的子程序给出。 算法实现： （1） 输入关系矩阵M （M 为n 阶方阵）； （2） 判断对称性，对于i=2，3，….，n ；j=1，2,……，i-1，若存在m ij =m ji ， 则R 是对称的； （3） 判断反对称性； （4） 判断既是对称的又是反对称的； （5） 判断既不是对称的又不是反对称的； （6） 输出判断结果。

自反： 从给定的关系矩阵来断判关系R是否为自反是很容易的。若M（R的关系矩阵）的主对角线元素均为1，则R是自反关系；若M（R的关系矩阵）的主对角线元素均为0，则R是反自反关系；若M（R的关系矩阵）的主对角线元素既有1又有0，则R既不是自反关系也不是反自反关系。本算法可以作为判等价关系算法的子程序给出。 算法实现 （1）输入关系矩阵M（M为n阶方阵）。 （2）判断自反性，对于i=1，2，….，n；若存在m =0，则R不是自反 ii =1，则R是自反的；否则R既不是自反关系也不是的；若存在m ii 反自反关系。 （3）输出判断结果。 源代码 #include void z(); void r(); void main() { int d; while(d) { printf("欢迎使用关系性质的判断系统\n\n 1. 对称关系的判断 2. 自反关系的判断\n\n请输入选项："); scanf("%d",&d); switch(d){ case 1: r();break; case 2: z();break; case 0: break; }

数学实验综合实验报告材料

一、实验目的： 1、初步认识迭代，体会迭代思想的重要性。 2、通过在mathematica环境下编写程序，利用迭代的方法求解方程的根、线性方程组的解、非线性方程组的解。 3、了解分形的的基本特性及利用mathematica编程生成分形图形的基本方法，在欣赏由mathematica生成的美丽的分形图案的同时对分形几何这门学科有一个直观的了解。从哲理的高度理解这门学科诞生的必然性，激发读者探寻科学真理的兴趣。 4、从一个简单的二次函数的迭代出发，利用mathematica认识混沌现象及其 所 蕴涵的规律。 5、.进一步熟悉Mathematic软件的使用，复习总结Mathematic在数学作图中的应用，为便于研究数学图像问题提供方便，使我们从一个新的视角去理解数学问题以及问题的实际意义。 6、在学习和运用迭代法求解过程中，体会各种迭代方法在解决问题的收敛速度上的异同点。 二、实验的环境： 学校机房，mathematica4环境 三、实验的基本理论和方法： 1、迭代（一）—方程求解 函数的迭代法思想： 给定实数域上光滑的实值函数)(x f以及初值 x定义数列

1()n n x f x +=， ,3,2,1,0=n ， (1) n x ， ,3,2,1,0=n ，称为)(x f 的一个迭代序列。 （1）方程求根 给定迭代函数)(x f 以及初值0x 利用（1）迭代得到数列n x ， ,3,2,1,0=n .如果数列收敛到某个*x ，则有 )(**x f x =. (2) 即*x 是方程)(x f x =的解。由此启发我们用如下的方法求方程0)(=x g 的近似解。 将方程0)(=x g 改写为等价的方程 )(x f x =， (3) 然后选取一初值利用（1）做迭代。迭代数列n x 收敛的极限就是方程0)(=x g 的解。 为了使得迭代序列收敛并尽快收敛到方程0)(=x g 的某一解的条件是迭代函数)(x f 在解的附近的导数将的绝对值尽量小，因此迭代方程修订成 x x f x h x )1()()(λλ-+== (4) 选取λ使得|)(|x h '在解的附近尽量小. 为此, 我们可以令 ,01)()(=-+'='λλx f x h 得 ) (11 x f '-= λ. 于是 1 )()()(-'-- =x f x x f x x h .

离散数学实验报告

《离散数学》 实验报告 题目 专业 学号 姓名 指导教师 提交日期

实验一五种连结词的逻辑运算 一.实验目的 用C语言实现两个命题变元的合取、析取、蕴涵和等价表达式的计算。熟悉连接词逻辑运算规则，利用程序语言实现逻辑这几种逻辑运算。 二.实验内容 从键盘输入两个命题变元P和Q的真值，求它们的合取、析取、蕴涵和等价四种运算的的真值。要求对输入内容进行分析，如果不符合0、1条件需要重新输入，程序有良好的输入输出界面。 三. 实验过程 1. 算法分析： 编程语言为c语言 合取/\：p，q都为1的时候为1，其他为0 析取\/：p，q都为0的时候为0，其他为1 蕴含->：p为1，q为0时为0，其他为1 等价<->：p，q同真同假 流程图

2. 程序代码： #include int main() { int p,q,i,t; printf("************************************************\n"); printf("*** ***\n"); printf(" 欢迎进入逻辑运算软件\n"); printf("*** ***\n"); printf("************************************************\n"); do{ printf("请输入p的值(0或1)"); scanf("%d",&p); if(p!=0&&p!=1) printf("输入有误"); }while(p!=0&&p!=1);

do{ printf("请输入q的值(0或1)"); scanf("%d",&q); if(q!=0&&q!=1) printf("输入有误"); }while(q!=0&&q!=1); do{ printf("请选择要进行的操作\n"); printf("1:合取\n2:析取\n3:蕴含\n4:等价\n"); scanf("%d",&i); switch(i){ case 1:{ if(p&&q) printf("合取运算：p/\q=1\n"); else printf("合取运算：p/\q=0\n"); break; } case 2:{ if(p||q) printf("析取运算：p\/q=1\n"); else printf("析取运算：p\/q=0\n"); break; } case 3:{ if(p&&!q) printf("蕴含：p->q=0\n"); else printf("蕴含：p->q=1\n"); break;} case 4:{ if((p&&q)||(!p&&!q)) printf("等价运算：p<->q=1\n"); else printf("等价运算：p<->q=0\n"); break; } }printf("是否继续运算1\\0\n"); scanf("%d",&t); }while(t); return 0; }