数学软件与数学实验报告

数学软件与数学实验报告
数学软件与数学实验报告

《数学软件与数学实验》报告

任课教师: 学号: 姓名:

成绩:

一、 实验内容:用Mathematica 软件求解下面各题 问题一:

任意拿出黑白两种颜色的棋子共八个排成一个圆圈, 然后在两颗相同的棋子中间放一颗黑色棋子, 在两颗不同颜色的棋子中间放一颗白色棋子,放完后撤掉原来所放的棋子, 再重复上述过程, 各棋子颜色会怎样变化?

同时,对棋子颜色问题进行深入探索 (对棋子数分别对4, 5, 6, 9, 16 颗) 的情况进行讨论, 将得到的结果进行拓广,并给出结论。

问题分析:

由于在两颗同色棋子中放置一颗黑色棋子,两颗不同色的棋子中间放置一颗白色棋子,因此可将黑色棋子用1表示,白色棋子用-1表示。这是因为-1-1=1? ,11=1? ,这代表两颗同色棋子中放一颗黑色棋子;-11=-1? ,这代表两颗不同色的棋子中间放置一颗白色棋子。

设棋子数为123,,,,

,n n a a a a 为初始状态,则当4n =时,有

第一次操作后得到的4枚棋子颗表示为

()()()()

12233441,,,a a a a a a a a 第二次操作后得到的4枚棋子颗表示为

()()()()()()()()

1234233434414112,,,a a a a a a a a a a a a a a a a 分别化简()()()()13243142,,,a a a a a a a a 第三次操作后得到的4枚棋子颗表示为

()()()()()()()()

1324243131424213,,,a a a a a a a a a a a a a a a a 化简为()

1234a a a a 第四次操作后得到的4枚棋子颗表示为()2

1234a a a a ,故这四枚棋子的颜色的赋值都是1,这表明只需操作4次即可将圆周上的棋子全变为黑色。 当n=4时,有 步数 状态(舍掉偶次项)

0 1a 2a 3a 4a 1 12a a 23a a 34a a 41a a

2 13a a

24a a 31a a

42a a 3 1234a a a a 1234a a a a

1234a a a a

1234a a a a

4

()

2

1234a a a a

()

2

1234a a a a

()

2

1234a a a a

()

2

1234a a a a

说明当n=4时,经过4部全变为黑色棋子 Mathematica 程序:

我们假设开始摆放的4颗棋子为1234a a a a 。由于棋子的变化规律性很强,并且为了明确起见,需要将每一轮所放的结果也记录下来。而且有3k k a a = ()1,2,3,4k = 。于是编制了如下程序:

当有4颗棋子进行试验: Mathematica 程序: In[1]: s = {};

x = {a1, a2, a3, a4};

For[i = 1, i <= Length[x], i++, x = Table[

x[[j]] x[[Mod[j, Length[x]] + 1]], {j, 1, Length[x]}] /. {a1^2 -> 1, a2^2 -> 1, a3^2 -> 1, a4^2 -> 1, a1^3 -> a1, a2^3 -> a3, a4^3 -> a4}; s = Append[s, x]]; m = TableForm[s, {1, Length[s]}]; Print["当n=", Length[x], "状态为"]; Print[m];

Print["当n=", Length[x], "时,做多经过", Length[s], "步,可将其棋子全部变为黑色"]

Out[1]: 当n=4状态为

a1a2a2a3a3a4a1a4a1a3a2a4a1a3a2a4

a1a2a3a4a1a2a3a4a1a2a3a4a1a2a3a4

1111

当=4时,做多经过4步,可将其棋子全部变为黑色; 准确结果如图1

图1

经过上述实验可以得知,列表最后一行全都是1,则无论初始状态棋子颜色如何,经过四次操作后,四颗棋子的颜色最终会变成黑色。

当n=5时,有

步数 状态(舍掉偶次项) 0 1a 2a 3a 4a 5a 1 12a a 23a a 34a a 45a a 51a a

2 13a a

24a a 35a a 14a a

25a a 3 1234a a a a 2345a a a a 1345a a a a

1245a a a a 1245a a a a 4 15a a

12a a 23a a 34a a 45a a 5 25a a 13a a

24a a 35a a 14a a

6

1234a a a a

2345a a a a

1345a a a a

1245a a a a

1245a a a a

则出现循环不能全变为黑子

当有5颗棋子进行试验: Mathematica 程序: In[2]: s = {};

x = {a1, a2, a3, a4, a5};

For[i = 1, i <= Length[x], i++, x = Table[

x[[j]] x[[Mod[j, Length[x]] + 1]], {j, 1, Length[x]}] /. {a1^2 -> 1, a2^2 -> 1, a3^2 -> 1, a4^2 -> 1, a5^2 -> 1, a1^3 -> a1, a2^3 -> a3, a4^3 -> a4, a5^3 -> a5}; s = Append[s, x]]; m = TableForm[s, {1, Length[s]}];

Print["当n=", Length[x], "状态为"]; Print[m] 结果为: 当n=5状态为 a1a2a2a3a3a4a4a5a1a5

a1a3a2a4a3a5a1a4a2a5a1a2a3a4a2a3a4a5a1a3a4a5a1a2a4a5a1a2a3a5a1a5a1a2a2a3a3a4a4a5a2a5

a1a3

a2a4a3a5

a1a4

在经过5次像以上的重复操作后, 棋子颜色出现{25a a ,13a a ,24a a ,35a a ,14a a }循环,即使做再多的变换,棋子颜色也不会全变成黑色。

当有5颗棋子进行试验: Mathematica 程序:s = {};

x = {a1, a2, a3, a4, a5, a6}; For[i = 1, i <= Length[x], i++, x = Table[

x[[j]] x[[Mod[j, Length[x]] + 1]], {j, 1, Length[x]}] /. {a1^2 -> 1, a2^2 -> 1, a3^2 -> 1, a4^2 -> 1, a5^2 -> 1, a6^2 -> 1, a1^3 -> a1, a2^3 -> a3, a4^3 -> a4, a5^3 -> a5, a6^3 -> a6}; s = Append[s, x]]; m = TableForm[s, {1, Length[s]}]; Print["当n=", Length[x], "状态为"]; Print[m] 结果为: 当n=6状态为 a1a2a2a3a3a4a4a5a5a6

a1a6a1a3a2a4a3a5a4a6a1a5a2a6a1a2a3a4a2a3a4a5a3a4a5a6a1a4a5a6a1a2a5a6a1a2a3a6

a1a5a2a6a1a3a2a4a3a5a4a6a1a2a5a6a1a2a3a6a1a2a3a4a2a3a4a5a3a4a5a6a1a4a5a6a3a5

a4a6a1a5a2a6

a1a3a2a4

在经过6次像以上的重复操作后, 棋子颜色出现354615261324{,,,,,}a a a a a a a a a a a a 循环,即使做再多的变换,棋子颜色也不会全变成黑色。

当有8颗棋子进行试验: 结果为:

当n=8时,做多经过8步,可将其棋子全部变为黑色;

当n=9时,结果图为

在经过9次像以上的重复操作后, 棋子颜色出现

{}291324354657687918,,,,,,,,a a a a a a a a a a a a a a a a a a 循环,即使做再多的变换,棋子颜

色也不会全变成黑色。 当n=16时,结果图为

当n=16时,做多经过16步,可将其棋子全部变为黑色;

总而言之,我们对棋子数分别对4, 5, 6,8, 9, 16 颗棋子进行试验后,发现只有4,8,16颗棋子在有限次循环操作后变成全黑色。其他的情况均不能达到这种效

果。于是我们可以得出结论:当棋子数为2n时,至多经过n2次操作即可将其全部变为黑色棋子;当棋子数不为2n次时,则一般不能全变为黑子。

问题二:

围绕着山洞有10个洞, 一只兔子和一只狐狸分别住在洞里, 狐狸总想吃掉兔子, 一天兔子对狐狸说 : 你想吃掉我有一个条件, 先把洞顺序编号,你从最后一个洞出发, 第一次先到第一个洞找我, 第二次隔一个洞找我, 第三次隔两个洞找我, 第四次隔三个洞找我, ...,依次类推, 寻找次数不限; 我躲在一个洞里不动; 只要找到我就可以饱餐了.结果狐狸跑断了腿也没找到兔子, 请问兔子躲在哪个洞里?同时,对山洞个数问题进行深入探索 (山洞数分别对12, 14, 17, 20, 30 等不同情况)将得到的结果进行拓广并给出结论。

问题的分析:

围绕着山洞的10个洞文成一个圈,,兔子和狐狸各占据一个洞,兔子呆在自己的洞里不动,狐狸在编好序号的洞内按照一定的规则找兔子,狐狸首先从最后一个洞出发,第一次先到第一个洞找兔子,第二次隔一个洞找兔子,第三次隔两个洞找兔子,以此类推,寻找次数不限。

将山洞编号,根据狐狸与兔子的约定产生狐狸所到洞的号的序列。最后一个洞的序号是10,隔n 个洞的话就是前面的洞号加n + 1,狐狸首先从最后一个洞( 10 号) 出发,第一次先到第一个洞( 1号) 找,第二次隔一个洞( 1 + 2 = 3 号),找第二次隔两个洞( 3 + 3 = 6 号) 找,以后狐狸运动到一个新洞的序号与前一次到的洞的序号有关,是前面洞的序号加隔几个洞数所得的数除以10 得到的余

+=,3+3=6,6+4=10,10+5=5,5+6=1,数。刚开始找过的洞号依次是0+1=1,123

1+7=8,8+8=6?

由于在寻找的过程当中需要经过许多次,所以我们利用计算机编程比较方便,对于10 号洞,也可以叫0号洞,这样做可以使程序的编制方便。

编制程序思路如下:

(1)初始值h = 0;开始时找过的洞的集合为t = {0}

(2)用For 循环实现前1000次的序号集;

(3)循环体当中,第i次的洞号h是前一次的洞号加i 除以10 得到的余数,同时将每次产生的序号加入以前的序号集当中。

Mathematica程序和结果如下:

In[2]:h = l = 10;

t = {0};

For[i = 1, i <= 1000, i++, h = Mod[h + i, l]; t = Append[t, h]];

t; m = Union[t]; x = Table[i, {i, 0, l-1}]; p =

Complement[x, m]; Print["兔子可能躲在的洞的号码为", p]

结果为:

兔子可能躲在的洞的号码为{2,4,7,9}

因此兔子可能躲在2号,4号,7号,9号洞中。

当山洞数为12时,只需修改程序中h=l=10的值为12即可,相应的得到的程序和结果为

当山洞数为14时相应的程序和结果为

当山洞数为17时相应的程序和结果为

当山洞数为20时相应的程序和结果为

当山洞数为30时相应的程序和结果为

容易看出狐狸找洞的编号序列为()(1)

2

k k k a Mod N +=,如果0k a =,洞的编号为N ,若0k a ≠ 的话k a 就是洞的编号

当N 为奇数时,令k=N-1就得到0k a =,也就是说寻找N-1次必回到编号为N 的洞中,接下去的序列是跟前面重复的。 所以狐狸可以找到的洞的编号就是12,,

,n a a a ,这个编号数目比洞的数目N 少(有可能不止小1),也就是说N 为奇数的时候这样的洞是肯定存在的。 当N 为偶数时,令k=2N-1,同样发现洞的编号只可能在()1221,,

,n a a a -这2N-1个

数中产生,这样有可能无解,例如N=4,兔子不管怎么藏都不行的。 按照N 的奇偶性,就可以找出狐狸能找到的所有洞的编号,从而找到兔子可能藏到的洞中。

二、 实验感想

从一个非常有趣的小故事出发,利用数学软件Mathematica 进行了求解,通过运算观察发现其中包含的数学规律性,并加以推广得出结论。这种学习模式是不仅对于我们学习数学软件有重要的帮助,更能通过这种方法加强我们对数学软件的

应用,开阔扩展我们的思维具,还可以强化我们对所学知识的巩固,这种思路在学习中是非常重要的,而且将现实问题数学化也是应用数学方向解决问题的有效途径,因此我们应该多努力学习数学软件并加以应用。

数学实验报告

《数学实验》实验报告 实验四 MATLAB 的作图功能 1、画出y=x+cosx 在[02]π,上的图形。 >> x=linspace(0,0.1,30); >> y=x+cos(x); >> plot(x,y) 1234567 2、在同一坐标系中作出两曲线y=tanx 、y=x-cosx 、2 y x =、2 1y x =-在[0]π,上的图形;要求曲线分别用虚实线表示,并注明曲线名称及适当的标注。 x=0:0.1:pi; y1=tan(x); y2=x-cos(x); y3=x.*x; y4=1-x.*x; plot(x,y1,'k-',x,y2,'k:',x,y3,'k-.',x,y4,'k--'); title('四条平面曲线'); gtext('y=tantx'); gtext('y=x-cosx'); gtext('y=x^2'); gtext('y=1-x^2 ');

0.5 1 1.5 2 2.5 3 3.5 -35-30-25-20-15-10-505 10 15四条平面曲线 3、22 2351 ,cos ,21,1 x x x y e z x u x v x +-===-=+将在同一窗口画出图形。 >> x=linspace(0,2*pi,30); >> y=exp(x); z=cos(x); u=2*x.^2-1; v=(3*x.*x+5*x-1)./(x.*x+1); >> subplot(2,2,1),plot(x,y),title('y=e^x') >> subplot(2,2,2),plot(x,z), title('y=cosx') >> subplot(2,2,3),plot(x,u), title('y=2x^2-1') >> subplot(2,2,4),plot(x,v), title('y=(3*x^2+5*x-1)/(x^2+1)')

数学软件MATLAB实验报告 实验八

实验八:概率论与数理统计的MATLAB 实现 实验目的与要求: 能运用MATLAB 提供的针对概率统计课程的工具箱。 实验内容: 1、用normpdf函数计算正态概率密度函数。 该函数的调用格式为:Y=normpdf(X,MU,SIGMA) 2、用normpdf函数计算正态分布的分布函数。 该函数的调用格式为:F=normcdf(X,MU,SIGMA) 3、用chi2inv函数计算卡方分布的分布函数的逆函数。 分布函数的逆函数及其调用格式:x=chi2inv(P,v) 4、随机取8只活塞环,测得他们直径为(以mm计): 74.001 74.005 74.003 74.001 74.000 73.998 74.006 74.002 。 设环直径的测量值服从正态分布,现估计总体的方差2 程序代码: x=[74.001 74.005 74.003 74.001 74.000 73.998 74.006 74.002]; p=mle('norm',x); sigma2hatmle=p(2)^2 5、从一批灯泡中随机的取5只做寿命试验,测得寿命(以小时计)为: 1050 1100 1120 1250 1280 设灯泡寿命服从正态分布,求灯泡寿命平均值的95%置信区间。 程序代码: x=[1050 1100 1120 1250 1280]; [p,ci]=mle('norm',x,0.05) 6、下面列出的是某工厂随机选取的20只部件的装配时间(分): 9.8 10.4 10.6 9.6 9.7 9.9 10.9 11.1 9.6 10.2

10.3 9.6 9.9 11.2 10.6 9.8 10.5 10.1 10.5 9.7 设装配时间的总体服从正态分布,标准差为0.4,是否可以认为装配时间的均值在0.05的水平上不小于10. 0H :10<μ vs 1H :10≥μ 程序: %正态总体的方差已知时的均值检验 x1=[9.8 10.4 10.6 9.6 9.7 9.9 10.9 11.1 9.6 10.2]; x2=[10.3 9.6 9.9 11.2 10.6 9.8 10.5 10.1 10.5 9.7]; x=[x1 x2]'; m=10;sigma=0.4;a=0.05; [h,sig,muci]=ztest(x,m,sigma,a,1) 因此,在0.05的水平下,可以认为装配时间的均值不小于10。 7、某种电子元件的寿命x (以小时计)服从正态分布,2 δμ和均未知。现测得16只元件的寿命如下: 159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170 问是否有理由认为元件的平均寿命大于225(小时)? 0H :225≤μ vs 1H :225>μ 程序: %正态总体的方差求知时的均值检验 x=[159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170]; m=225;a=0.05; [h,sig,muci]=ttest(x,m,a,1)

东南大学高等数学数学实验报告上

Image Image 高等数学数学实验报告 实验人员:院(系) ___________学号_________姓名____________实验地点:计算机中心机房 实验一 1、 实验题目: 根据上面的题目,通过作图,观察重要极限:lim(1+1/n)n =e 2、 实验目的和意义 方法的理论意义和实用价值。 利用数形结合的方法观察数列的极限,可以从点图上看出数列的收敛性,以及近似地观察出数列的收敛值;通过编程可以输出数列的任意多项值,以此来得到数列的收敛性。通过此实验对数列极限概念的理解形象化、具体化。 三、计算公式 (1+1/n)n 四、程序设计 五、程序运行结果 六、结果的讨论和分析 当n足够

Image Image 大时,所画出的点逐渐接近于直线,即点数越大,精确度越高。对于不同解题方法最后均能获得相同结果,因此需要择优,从众多方法中尽可能选择简单的一种。程序编写需要有扎实的理论基础,因此在上机调试前要仔细审查细节,对程序进行尽可能的简化、改进与完善。 实验二一、实验题目 制作函数y=sin cx的图形动画,并观察参数c对函数图形的影响。 二、实验目的和意义 本实验的目的是让同学熟悉数学软件Mathematica所具有的良好的作图功能,并通过函数图形来认识函数,运用函数的图形来观察和分析函数的有关性态,建立数形结合的思想。三、计算公式:y=sin cx 四、程序设计五、程序运行结果 六、结果的讨论和分析 c的不同导致函数的区间大小不同。 实验三 一、实验题目 观察函数f(x)=cos x的各阶泰勒展开式的图形。 二、实验目的和意义 利用Mathematica计算函数的各阶泰勒多项式,并通过绘制曲线图形,来进一步掌握泰勒展开与函数逼近的思想。 三、计算公式

数学模型实验报告

数学模型实验报告 实验内容1. 实验目的:学习使用lingo和MATLAB解决数学模型问题 实验原理: 实验环境:MATLAB7.0 实验结论: 源程序 第4章:实验目的,学会使用lingo解决数学模型中线性规划问题1.习题第一题 实验原理: 源程序: 运行结果: 、 管 路 敷 设 技 术 通 过 管 线 不 仅 可 以 解 决 吊 顶 层 配 置 不 规 范 高 中 资 料 试 卷 问 题 , 而 且 可 保 障 各 类 管 路 习 题 到 位 。 在 管 路 敷 设 过 程 中 , 要 加 强 看 护 关 于 管 路 高 中 资 料 试 卷 连 接 管 口 处 理 高 中 资 料 试 卷 弯 扁 度 固 定 盒 位 置 保 护 层 防 腐 跨 接 地 线 弯 曲 半 径 标 等 , 要 求 技 术 交 底 。 管 线 敷 设 技 术 中 包 含 线 槽 、 管 架 等 多 项 方 式 , 为 解 决 高 中 语 文 电 气 课 件 中 管 壁 薄 、 接 口 不 严 等 问 题 , 合 理 利 用 管 线 敷 设 技 术 。 线 缆 敷 设 原 则 : 在 分 线 盒 处 , 当 不 同 电 压 回 路 交 叉 时 , 应 采 用 金 属 隔 板 进 行 隔 开 处 理 ; 同 一 线 槽 内 强 电 回 路 须 同 时 切 断 习 题 电 源 , 线 缆 敷 设 完 毕 , 要 进 行 检 查 和 检 测 处 理 。 、 电 气 课 件 中 调 试 对 全 部 高 中 资 料 试 卷 电 气 设 备 , 在 安 装 过 程 中 以 及 安 装 结 束 后 进 行 高 中 资 料 试 卷 调 整 试 验 ; 通 电 检 查 所 有 设 备 高 中 资 料 试 卷 相 互 作 用 与 相 互 关 系 , 根 据 生 产 工 艺 高 中 资 料 试 卷 要 求 , 对 电 气 设 备 进 行 空 载 与 带 负 荷 下 高 中 资 料 试 卷 调 控 试 验 ; 对 设 备 进 行 调 整 使 其 在 正 常 工 况 下 与 过 度 工 作 下 都 可 以 正 常 工 作 ; 对 于 继 电 保 护 进 行 整 核 对 定 值 , 审 核 与 校 对 图 纸 , 编 写 复 杂 设 备 与 装 置 高 中 资 料 试 卷 调 试 方 案 , 编 写 重 要 设 备 高 中 资 料 试 卷 试 验 方 案 以 及 系 统 启 动 方 案 ; 对 整 套 启 动 过 程 中 高 中 资 料 试 卷 电 气 设 备 进 行 调 试 工 作 并 且 进 行 过 关 运 行 高 中 资 料 试 卷 技 术 指 导 。 对 于 调 试 过 程 中 高 中 资 料 试 卷 技 术 问 题 , 作 为 调 试 人 员 , 需 要 在 事 前 掌 握 图 纸 资 料 、 设 备 制 造 厂 家 出 具 高 中 资 料 试 卷 试 验 报 告 与 相 关 技 术 资 料 , 并 且 了 解 现 场 设 备 高 中 资 料 试 卷 布 置 情 况 与 有 关 高 中 资 料 试 卷 电 气 系 统 接 线 等 情 况 , 然 后 根 据 规 范 与 规 程 规 定 , 制 定 设 备 调 试 高 中 资 料 试 卷 方 案 。 、 电 气 设 备 调 试 高 中 资 料 试 卷 技 术 电 力 保 护 装 置 调 试 技 术 , 电 力 保 护 高 中 资 料 试 卷 配 置 技 术 是 指 机 组 在 进 行 继 电 保 护 高 中 资 料 试 卷 总 体 配 置 时 , 需 要 在 最 大 限 度 内 来 确 保 机 组 高 中 资 料 试 卷 安 全 , 并 且 尽 可 能 地 缩 小 故 障 高 中 资 料 试 卷 破 坏 范 围 , 或 者 对 某 些 异 常 高 中 资 料 试 卷 工 况 进 行 自 动 处 理 , 尤 其 要 避 免 错 误 高 中 资 料 试 卷 保 护 装 置 动 作 , 并 且 拒 绝 动 作 , 来 避 免 不 必 要 高 中 资 料 试 卷 突 然 停 机 。 因 此 , 电 力 高 中 资 料 试 卷 保 护 装 置 调 试 技 术 , 要 求 电 力 保 护 装 置 做 到 准 确 灵 活 。 对 于 差 动 保 护 装 置 高 中 资 料 试 卷 调 试 技 术 是 指 发 电 机 一 变 压 器 组 在 发 生 内 部 故 障 时 , 需 要 进 行 外 部 电 源 高 中 资 料 试 卷 切 除 从 而 采 用 高 中 资 料 试 卷 主 要 保 护 装 置 。

MATLAB实验报告实验二

实验二 MATLAB矩阵及其运算 学号:3121003104 姓名:刘艳琳专业:电子信息工程1班日期:2014.9.20 一实验目的 1、掌握Matlab数据对象的特点以及数据的运算规则。 2、掌握Matlab中建立矩阵的方法以及矩阵处理的方法。 3、掌握Matlab分析的方法。 二实验环境 PC_Windows 7旗舰版、MATLAB 7.10 三实验内容 4、1. (1)新建一个.m文件,验证书本第15页例2-1; (2)用命令方式查看和保存代码中的所有变量;

(3)用命令方式删除所有变量; (4)用命令方式载入变量z。 2. 将x=[4/3 1.2345e-6]在以下格式符下输出:短格式、短格式e方式、长格式、长格式e方式、银行格式、十六进制格式、+格式。 短格式 短格式e 长格式

长格式e方式 银行格式 十六进制格式 3.计算下列表达式的值 (1)w=sqrt(2)*(1+0.34245*10^(-6)) (2)x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) a=3.5;b=5;c=-9.8; (3)y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) a=3.32;b=-7.9; (4)z=0.5*exp(2*t)*log(t+sqrt(1+t*t)) t=[2,1-3i;5,-0.65];

4. 已知A=[1 2 3 4 5 ;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20],对其进行如下操作:(1)输出A在[ 7, 10]范围内的全部元素; (2)取出A的第2,4行和第1,3,5列; (3)对矩阵A变换成向量B,B=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]; (4)删除A的第2,3,4行元素; (1) (2)

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

Matlab数学实验报告一

数学软件课程设计 题目非线性方程求解 班级数学081 姓名曹曼伦

实验目的:用二分法与Newton迭代法求解非线性方程的根; 用Matlab函数solve、fzero、fsolve求解非线性方程(组)的解。 编程实现二分法及Newton迭代法; 学会使用Matlab函数solve、fzero、fsolve求解非线性方程(组)的解。 通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。 实验内容: 比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方) (1)在区间(0,1)内用二分法; (2)用迭代法x=(2-exp(x))/10,取初值x=0 。 试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0,则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正,则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负,则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法; format long

syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果: (1)二分法: symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 (2)迭代法 由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次,而用迭代法求得同样精度的解仅用5次,但由于迭代法一般只具有局部收敛性,因此通常不用二分法来求得非线性方程的精确解,而只用它求得根的一个近似解,再用收敛速度较快的迭代法求得其精确解。

算法与分析实验报告模板

贵州大学计算机科学与技术学院 计算机科学与技术系上机实验报告 课程名称:算法设计与分析班级:实验日期:YYYY-MM-DD 姓名:学号:指导教师:程欣宇 实验序号:一实验成绩: 一、实验名称 分治算法实验- 棋盘覆盖问题 二、实验目的及要求 1、熟悉递归算法编写; 2、理解分治算法的特点; 3、掌握分治算法的基本结构。 三、实验环境 Visual C++ 四、实验内容 根据教材上分析的棋盘覆盖问题的求解思路,进行验证性实验; 要求完成棋盘覆盖问题的输入、分治求解、输出。有余力的同学尝试消去递归求解。 五、算法描述及实验步骤 分治算法原理: 分治算法将大的分解成形状结构相同的子问题,并且不断递归地分解,直到子问题规模小到可以直接求解。 棋盘覆盖问题描述: 在一个2k x 2k个方格组成的棋盘中恰有一个方格与其他的不同称为特殊方格,想要求利用四种L型骨牌(每个骨牌可覆盖三个方格)不相互重叠覆盖的将除了特殊方格外的其他方格覆盖。

实验步骤: 1、定义用于输入和输出的数据结构; 2、完成分治算法的编写; 3、测试记录结构; 4、有余力的同学尝试不改变输入输出结构,将递归消除,并说明能否不用栈,直接消除递归,为什么? 六、调试过程及实验结果 详细记录程序在调试过程中出现的问题及解决方法。 记录程序执行的结果。 七、总结 对上机实践结果进行分析,问题回答,上机的心得体会及改进意见。 八、附录 源程序(核心代码)清单或使用说明书,可另附纸

贵州大学计算机科学与技术学院 计算机科学与技术系上机实验报告 课程名称:算法设计与分析班级:实验日期:2014-11-25 姓名:学号:指导教师:程欣宇 实验序号:二实验成绩: 一、实验名称 动态规划实验- 滑雪问题 二、实验目的及要求 1、学会使用在线测评的算法题目评分系统; 2、通过直观的应用问题,加深对动态规划算法的理解; 三、实验环境 任意C或C++编写调试工具,北京大学ICPC在线测评系统POJ 四、实验内容 1、找到题号为1088的题目-滑雪,阅读题目,建立其最优解的递归表达式; 3、使用备忘录式的动态规划算法,实现本题; 4、进行简单测试,完成之后提交到POJ系统。 五、算法描述及实验步骤 动态规划算法原理: 分治算法将大的问题变成小的问题来解决,但是如果划分过程中出现重叠子问题,就可能导致大量的重复计算。为了避免这些重复的计算,可以考虑的一个办法就是动态规划算法。 为了使用动态规划算法,问题还必须具备最优子结构,即问题的最优解包含了子问题的最优解。 滑雪问题描述: Michael喜欢滑雪百这并不奇怪,因为滑雪的确很刺激。可是为了获得速度,滑的区域必须向下倾斜,而且当你滑到坡底,你不得不再次走上坡或者等待升降机来载你。Michael想知道载一个区域中最长底滑坡。区域由一个二维数组给出。数组的每个数字代表点的高度。下面是一个例子 1 2 3 4 5 16 17 18 19 6 15 24 25 20 7 14 23 22 21 8 13 12 11 10 9 一个人可以从某个点滑向上下左右相邻四个点之一,当且仅当高度减小。在

小学数学实验报告

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

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

四川师范大学数学与软件科学学院程序设计实验报告实验九(推荐文档)

数学与软件科学学院实验报告 一、实验目的 (1) 掌握C语言环境下结构体和共用体类型变量的定义和使用方法; (2) 掌握结构体类型数组的概念和使用; (3) 掌握指向结构体变量的指针变量、尤其是链表概念; 二、实验内容 1.首先熟悉结构体类型变量的基本声明方法、结构体类型变量的内存分配原则、初始化和引用结构体变量及其成员变量的基本方法;然后掌握结构体变量的输入、输出方法。(参见教材例7.1,请给该例加上输入功能) #include struct person { char name[20]; int count; }leader[3]={"Li",0,"Zhang",0,"Wang",0}; main() { int i,j; char leader_name[20]; for(i=1;i<=10;i++) { scanf("%s",leader_name); for(j=0;j<3;j++) { if(strcmp(leader_name,leader[j].name)==0) leader[j].count++; } }

for(i=0;i<3;i++) printf("%5s:%d\n",leader[i].name,leader[i].count); } 2.基于结构体数组的应用实验。 (1) 有n个学生,每个学生的数据包括学好(num)、姓名(name[20])、性别(sex)、年龄(age),以及三门课程的成绩(score[3])。要求:在main()函数中输入这些学生的这些数据,然后设计一个函数count()来计算每个学生的总分和平均分,最后, 打印出所有数据信息(包含原来输入的学生原始数据信息和求解出来的新信息)。#include #define N 3 #define M 3 typedef struct student { int score[N]; char name[20]; int sex; int age; char num[20]; }STUDENT; main() { STUDENT stu[M]; int i,j,average,total; char name[20]; clrscr(); for(i=0;i

高等数学下实验报告

高等数学实验报告 实验人员:院(系)化学化工学院 学号19013302 姓名 黄天宇 实验地点:计算机中心机房 实验七:空间曲线与曲面的绘制 一、 实验目的 1、利用数学软件Mathematica 绘制三维图形来观察空间曲线和空 间曲面图形的特点,以加强几何的直观性。 2、学会用Mathematica 绘制空间立体图形。 二、实验题目 利用参数方程作图,做出由下列曲面所围成的立体图形: (1) x y x y x z =+--=2 222,1及xOy 平面; (2) 01,=-+=y x xy z 及.0=z 三、实验原理 空间曲面的绘制 作参数方程],[],,[,),(),() ,(max min max min v v v u u v u z z v u y y v u x x ∈∈? ?? ??===所确定的曲面图形的 Mathematica 命令为: ParametricPlot3D[{x[u,v],y[u,v],z[u,v]},{u,umin,umax}, {v,vmin,vmax},选项] 四、程序设计及运行 (1)

(2)

六、结果的讨论和分析 1、通过参数方程的方法做出的图形,可以比较完整的显示出空 间中的曲面和立体图形。 2、可以通过mathematica 软件作出多重积分的积分区域,使积分能够较直观的被观察。 3、从(1)中的实验结果可以看出,所围成的立体图形是球面和圆柱面所围成的立体空间。 4、从(2)中的实验结果可以看出围成的立体图形的上面曲面的方程是xy z =,下底面的方程是z=0,右边的平面是01=-+y x 。 实验八 无穷级数与函数逼近 一、 实验目的 (1) 用Mathematica 显示级数部分和的变化趋势; (2) 展示Fourier 级数对周期函数的逼近情况; (3) 学会如何利用幂级数的部分和对函数进行逼近以及函数值的近似计算。 二、实验题目 (1)、观察级数 ∑ ∞ =1 ! n n n n 的部分和序列的变化趋势,并求和。 (2)、改变例2中m 及x 0的数值来求函数的幂级数及观察其幂级数逼近函数的情况 (3)、观察函数? ? ?<≤<≤--=ππx x x x f 0,10 ,)(展成的Fourier 级数

数学应用软件实验报告(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

《数学实验》曲线绘制实验报告

课程名称数学实验成绩评定 实验项目名称曲线绘制 【实验目的】 1.了解曲线的几种表示方式。 2.学习、掌握MA TLAB软件有关的命令。 【实验内容】 绘制下列四种曲线: 1.以直角坐标方程y=sin x,y=cos x表示的正、余弦曲线。 2.以参数方程x=cos t,y=sin t,t∈[0,2π]表示的平面曲线(单位圆)。 3.以参数方程x=e?0.2t cosπ 2t,y=π 2 e?0.2t sin t,z=t,t∈[0,20]表示的空间曲线。 4.作出摆线的图形。 5.做出以参数方程x=e?0.25t cosπ 2t,y=e?0.25t sinπ 2 t,z=t,t∈[0,30]表示的空间曲线。 6.以极坐标方程r=a(1+cos?),a=1,?∈[0,2π]表示的心脏线。 7.绘制极坐标系下曲线 ρ=acos (b+nθ)的图形,讨论参数a、b和n对其图形的影响。8.(曲线族绘制)三次抛物线的方程为y=ax3+cx,讨论参数a和c对其图形的影响。 【实验方法与步骤】 练习1做出函数y=sin x,y=cos x的图形,并观察它们的周期性。 MATLAB代码及结果如下: >> x=0:0.01*pi:4*pi; y1=sin(x); y2=cos(x); plot(x,y1,'b',x,y2,'r'); legend('y=sin(x)','y=cos(x)','location','best'); axis([0 4*pi -1 1]) 绘制结果如下图:

y=sin x,y=cos x的图形如上图,两个函数的周期皆为2π 练习2设y=√3 2e?4t sin(4√3t+π 3 ),要求以0.01秒为间隔,求出y的151个点,绘出y及 其导数的图形。 MATLAB代码及结果如下: dt=0.01; t=0:0.01:1.5; w=4*sqrt(3); %设定频率 y=sqrt(3)/2*exp(-4*t).*sin(w*t+pi/3); Dy=diff(y)/dt; %求导 for i =1:length(t)-1 t1(i)=t(i); end subplot(2,1,1); plot(t,y); xlabel('时间t'); ylabel('y(t)'); grid subplot(2,1,2); plot(t1,Dy); xlabel('时间t'); ylabel('Dy(t)'' '); grid 绘制结果如下图:

数学社会实践报告-范文

数学社会实践报告 数学是研究数量、结构、变化、空间以及信息等概念的一门学科,本文将介绍数学社会实践报告。 数学社会实践报告(1) 又是一个酷热难耐的暑假,济南以它独特的天气特点招待了我们这些因为参赛而留在老校住宿的同学们,几次零星的小雨丝毫撼不动炎热的主题。蓊蓊郁郁的师大老校园里大批学子,他们忙碌着,早出晚归;他们埋头苦干着,废寝忘食;他们做着自己的事情,紧张有序他们默默等待着一场未知的洗礼。他们,就是参加暑假数学建模辅导的同学。 我很荣幸地成为了这支队伍中的一员,而且成为队长,本组成员都是让我佩服的两位很优秀的同学,让我对这次建模的胜利充满信心,宋希良,和王成龙,这两位我的员工,让我感觉很踏实,本来平淡无奇的暑假,因为参加了数学建模而变得丰富多彩。 先说说数学建模吧。数学建模是运用数学思想、方法和知识解决实际问题的过程,已经成为不同层次数学教育重要和基本的内容。数学建模是数学学习的一种新的方式,它为学生提供了自主学习的空间,有助于学生体验数学在解决实际问题中的价值和作用,体验数学与日常生活和其他学科的联系,体验综合运用知识和方法解决实际问题的过程,增强应用意识;有助于激发学生学习数学的兴趣,发展学生的创新意识和实践能力。数学建模与数学实验开创了大学生把数学理论和专业知识有机结合的新途径,

是培养学生分析问题、解决问题和使用计算机进行科学计算的有效方法,是培养学生创新能力和实践能力的有效手段。 中国科学院王梓坤院士在《今日数学及其应用》一文中指出精确定量思维是对21世纪科技人员的素质要求。所谓定量思维就是人们从实际问题中提炼数学问题,抽象化为数学模型,用数学计算此模型的解或近似解,然后回到现实中进行检验,必要时修改模型使之更切合实际,最后编制解决问题的软件包,以便得到更广泛的方便的应用。这一精辟的论述阐明了在解决工程实际问题中数学建模与数学实验是相互依赖、相辅相成、互不可分的。数学建模与数学实验是以数学知识为基础,以各个领域的实际问题为载体,以计算机为手段,以数学软件为工具,培养学生深入理解数学建模的思想与方法,熟悉常用的科学计算软件,如,Mathematica、MATLAB,并在此基础上,根据所要解决的数学问题进行程序设计,培养学生运用所学知识建立数学模型,使用计算机解决实际问题的能力,以及综合应用能力和创新能力。 建模前的准备。首先,要完善自己。只有解决了自身的问题,才能克服其他的问题。如果连自己都没把握好,那么,做任何事都会漏洞百出。要完善自己,首先要明确态度,记得中国前任国足教练米卢说过:态度决定一切。明确自己为什么要参加数学建模竞赛,参加的目的是什么,是抱着学习的态度参加呢还是其他呢?只有态度明确了,才能在这个前提下,进行全身心的投入竞赛。其次,要有热情,要有认真,严谨的科学精神。热情是动力的源

《数学软件》实验报告-符号计算基础与符号微积分

实验报告 课程名称:数学软件姓名: 学院: 专业: 年级: 学号: 指导教师: 职称: 年月日

实验项目列表

附件三: 实验报告(二) 系:专业:年级:姓名学号:实验课程: 实验室号:_ 实验设备号:实验时间: 指导教师签字:成绩: 1. 实验项目名称:符号计算基础与符号微积分 2. 实验目的和要求 1.掌握定义符号对象的方法 2.掌握符号表达式的运算法则以及符号矩阵运算 3.掌握求符号函数极限及其导数的方法 4.掌握求符号函数定积分和不定积分的方法 3. 实验使用的主要仪器设备和软件 方正商祺N260微机;MATLAB7. 0或以上版本 4. 实验的基本理论和方法 (1)符号函数;sym(x);syms a b …… (2)平方根:sqrt(x) (3)分解因式:factor(s) (4)符号表达式化简:simplify(s) (5)逆矩阵:inv(x) (6)下三角矩阵:tril(x) (7)矩阵行列式的值:det(x)

(8)符号函数求极限:limit (f ,x ,a );limit (f ,x ,a ,‘right ’) (9)符号函数求导:diff (f ,v ,n ) (10)符号函数求不定积分:int (f ,v ) (11)符号函数求定积分:int (f ,v ,a ,b ) 5. 实验内容与步骤 (描述实验中应该做什么事情,如何做等,实验过程中记录发生的现象、中间结果、最终得到的结果,并进行分析说明) (包括:题目,写过程、答案) 题目: 1. 已知x=6,y=5,利用符号表达式求 y x x z -++= 31。 提示:定义符号常数)'5(')'6('sym y sym x ==,。 >> x=sym('6'); >> y=sym('5'); >> z=(x+1)/(sqrt(3+x)-sqrt(y)) z = 7/(3-5^(1/2)) 2. 分解因式:44y x - >> syms x y; >> A=x^4-y^4; >> factor(A) ans = (x-y)*(x+y)*(x^2+y^2) 3. 化简表达式 (1)2121sin cos cos sin ββββ- (2) 123842+++x x x (1) >> syms x y; >> f1=sin(x)*cos(y)-cos(x)*sin(y);

数学实验报告

高等数学数学实验报告 实验人员:院(系) __ __学号____姓名_ __ 实验地点:计算机中心机房 实验一 空间曲线与曲面的绘制 一、实验题目:(实验习题1-2) 利用参数方程作图,做出由下列曲面所围成的立体图形: (1) x y x y x z =+--=2222,1及xOy 平面; (2) 01,=-+=y x xy z 及.0=z 二、实验目的和意义 1、利用数学软件Mathematica 绘制三维图形来观察空间曲线和空间曲面图形的特点,以加强几何的直观性。 2、学会用Mathematica 绘制空间立体图形。 三、程序设计 空间曲面的绘制 作参数方程] ,[],,[,),(),(),(max min max min v v v u u v u z z v u y y v u x x ∈∈?????===所确定的曲面图形的Mathematica 命令为: ParametricPlot3D[{x[u,v],y[u,v],z[u,v]},{u,umin,umax}, {v,vmin,vmax},选项] (1) (2)

四、程序运行结果 (1) (2) 五、结果的讨论和分析 1、通过参数方程的方法做出的图形,可以比较完整的显示出空间中的曲面和立体图形。 2、可以通过mathematica 软件作出多重积分的积分区域,使积分能够较直观的被观察。 3、从(1)中的实验结果可以看出,所围成的立体图形是球面和圆柱面所围成的立体空间。 4、从(2)中的实验结果可以看出围成的立体图形的上面曲面的方程是xy z =,下底面的方程是z=0,右边的平面是01=-+y x 。 实验一 空间曲线与曲面的绘制 一、实验题目:(实验习题1-3) 观察二次曲面族kxy y x z ++=22的图形。特别注意确定k 的这样一些值,当k 经过这些值时,曲面从一种类型变成了另一种类型。 二、实验目的和意义 1. 学会利用Mathematica 软件绘制三维图形来观察空间曲线和空间曲线图形的特

高等数学(下册)数学实验报告

高等数学A(下册)实验报告 院(系): 学号:姓名: 实验一 利用参数方程作图,作出由下列曲面所围成的立体: (1) 2 2 1Y X Z- - = , X Y X= +2 2 及 xOy 面 ·程序设计: -1, 1},Axe s2=ParametricPlot3D[{1/2*Cos[u]+1/2,1/2*Sin[u],v},{u,- s3=ParametricPlot3D[{u,v,0},{u,-1,1},{v,- DisplayFunction 程序运行结果: 实验二 实验名称:无穷级数与函数逼近 实验目的:观察的部分和序列的变化趋势,并求和

实验内容: (1)利用级数观察图形的敛散性 当n 从1~400时,输入语句如下: 运行后见下图,可以看出级数收敛,级数和大约为1.87985 (2先输入: 输出: 输出和输入相同,此时应该用近似值法。输入: 输出: 1.87985 结论:级数大约收敛于1.87985 实验三: 1. 改变例2中m 的值及的数值来求函数的幂级数及观察其幂级数逼近函数的情况

·程序设计: m 5; f x_:1 x^m;x0 1; g n_,x0_ :D f x, x, n .x x0; s n_,x_: Sum g k,x0/k x x0 ^k, k, 0, t Table s n, x, n, 20; p1 Plot Evaluate t ,x,1,2,3 2; p2 Plot 1 x ^m , x,1 2,3 2, PlotStyle RGBColor 0,0,1; Show p1,p2 ·程序运行结果 实验四 实验名称:最小二乘法 实验目的:测定某种刀具的磨损速度与时间的关系实验内容:

数学实验报告格式

《数学实验》实验报告 (2012 年03 月30 日) 班级:09级四班学号:姓名:吴永慧 一、实验问题 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、用???=城市人不去城市人去了k 0k 1 i i x ik (i =1...5) ???=城市人没去城市人去了h j h j x jh 01 (i =1...5) ij d 表示i 和j 的通话时间;kh c 表示城市k 和h 之间的费率,数学模型: min jh ik i j k h ij kh x x d c ∑∑∑∑====5151515 1 s.t.???????????========∑∑∑∑====5 151515 1 5 ...115...115...115 (11) h jh j jh k ik i ik j x k x i x k x ik x 、jh x 均为0、1变量 2、用???=该学生不选该课程该学生选了该课程01 i x (i =1...9) 1) 数学模型:min Z=∑=91i i x

相关文档
最新文档