数学建模作业——实验1

数学建模作业——实验1
数学建模作业——实验1

数学建模作业——实验1

学院:软件学院

姓名:

学号:

班级:软件工程2015级 GCT班

邮箱:

电话:

日期:2016年5月10日

基本实验

1.椅子放平问题

依照1.2.1节中的“椅子问题”的方法,将假设中的“四腿长相同并且四脚连线呈正方形”,改为“四腿长相同并且四脚连线呈长方形”,其余假设不变,问椅子还能放平吗?如果能,请证明;如果不能,请举出相应的例子。

答:能放平,证明如下:

如上图,以椅子的中心点建立坐标,O为原点,A、B、C、D为椅子四脚的初始位置,通过旋转椅子到A’、B’、C’、D’,旋转的角度为α,记A、B两脚,C、D两脚距离地面的距离为f(α)和g(α),由于椅子的四脚在任何位置至少有3脚着地,且f(α)、g(α)是α的连续函数,则f(α)和g(α)至少有一个的值为0,即f(α)g(α)=0,f(α)≥0,g(α)≥0,若f(0)>0,g(0)=0,

则一定存在α’∈(0,π),使得

f(α’)=g(α’)=0

令α=π(即椅子旋转180°,AB 边与CD 边互换),则

f(π)=0,g(π)>0

定义h(α)= f(α)-g(α),得到

h(0)=f(0)-g(0)>0

h(π)=f(π)-g(π) <0

根据连续函数的零点定理,则存在α’∈( 0,π),使得

h(α’)= f(α’)-g(α’)=0

结合条件f(α’)g(α’)=0,从而得到

f(α’)=g(α’)=0,即四脚着地,椅子放平。

2. 过河问题

依照1.2.2节中的“商人安全过河”的方法,完成下面的智力游戏:人带着猫、鸡、米过河,船除需要人划之外,至多能载猫、鸡、米之一,而当人不在场时,猫要吃鸡、鸡要吃米,试设计一个安全过河的方案,并使渡河的次数尽量的少。

答: 用i =1,2,3,4分别代表人,猫,鸡,米。1=i x 在此岸,0

=i x 在对岸,

()4321,,,x x x x s =此岸状态,()43211,1,1,1x x x x D ----=对岸状态。安全状态集合为 :

{ S =(1,1,1,1),D =(0,0,0,0)

S =(1,1,1,0),D =(0,0,0,1)S =(1,1,0,1),D =(0,0,1,0)S =(1,0,1,1),D =(0,1,0,0)S =(1,0,1,0),D =(1,0,1,0)S =(0,1,0,1),D =(0,1,0,1)S =(0,1,0,0),D =(1,0,1,1)S =(0,0,1,0),D =(1,1,0,1)S =(0,0,0,1),D =(1,1,1,0)S =(0,0,0,0),D =(1,1,1,1)

乘船方案,记作()4321,,,u u u u U =,当i 在船上时记1=i u ,否则记0=i u ,允许决策集合为

()()()(){}0,0,0,1,1,0,0,1,0,1,0,1,0,0,1,1=U

因为乘船k 为奇数时船从此案驶向彼岸,k 为偶数时船由彼岸驶向此岸,所以状态k s 随决策k U 变化的规律为

()k

U k

k s k s 11-+=+, 设计安全过河方案归结为求决策,,,,21U U U U n ∈ ,使状态S

s k ∈按转移律由初始状态()1,1,1,

11=s 经有限步n 到达状态()0,0,0,01=+n s 。 根据题设条件:影响安全渡河的元素是猫、鸡、米,这3个元素

中取2个元素的组合一共有C 32=3种,分别为猫+鸡,猫+米,鸡+米。

其中“猫+鸡”和“鸡+米”组合不安全,而“猫+米”的组合是安全的。第一次渡河需带3个元素中的1个元素,另外2个元素留在岸上,而留在岸上的3种组合中只有“猫+米”的组合是安全的,可见第一次渡河只能带鸡,即安全方案只有U=(1,0,1,0),第二次将米或者猫带到对岸,把鸡带回,第三次将猫或者米带到对岸,第四次将

鸡带过河,至此,猫、鸡、米均安全过河。具体有2种执行方案如下

上述方案直观表示如下:

3.购房贷款问题(续)

在1.2.3节“购房贷款”的问题中,我们讨论了小王夫妇借贷还贷的方式。现进一步讨论此问题。

某借贷公司的广告称,对于贷款期在20年以上的客户(利率仍按06%/月计算),他们帮你提前3年还清贷款。但条件如下:

(1)每半个月付款一次,但付款额不增加,即一次付款额是原付给银行还款额的1/2;

(2)因为增加必要的档案、文书等管理工作,因此,要预付给借贷

公司贷款

总额10%的佣金。

试分析,小王夫妇是否要请这家借贷公司帮助还款。

答:设每月还款额为x,利率r,贷款额A0,总还款月数为N,总还款额为X。若贷款20万,20年期,则:

x=A0r(1+r)N

(1+r)N?1=200000×0.006(1+0.006)

240

(1+0.006)240?1

=1574.7元

X=x×N=1574.7×240=377928元

若请这家借贷公司还款,则

还款总额X1=x(N-3×12)+ 10%A0=1574.7×204+20000=341238.8元

但小王夫妇预付给借贷公司20000元的佣金,如果把这20000元作为首付,则小王夫妇只需贷款18万,则总还款额:

X2=X(1-10%)=340135.2元

还款月数:

N2=X2÷x=340135.2÷1574.7=216月=18年

X1>X2,小王夫妇不需用借贷公司还款。

4.冷却定律

按照Newton冷却定律,温度为T的物体在温度为T0(T0

如果空气的温度是20℃,且沸腾的水在20分钟内冷却到60℃,那么水温降低到30℃需要多长时间?

答:首先,牛顿冷却定律为温度为T(t)的物体在温度T0的环境中冷却的速度与温度差 T(t)?T0成正比。

所以,得出微分方程dT(t)

dt

=k[T(t)?T0],K为比例常数。

任意时刻t,物体的温度为T(t)=T0+Ce kt,C为常数

根据已知条件,T0=20℃,记t=0时刻,初始水温T(0)=100℃,20分钟后水温T(20)=60℃。则:

{T(0)=20+Ce k?0=100 T(20)=20+Ce k?20=60

求解函数得,k=-0.034657,C=80,即

T(t)=20+80?e?0.034657?t

当水温降低到30℃时,t=-ln((30-20)÷80)÷0.034657=60分钟。5.锻炼想象力、洞察力和判断力的问题(只简单回答出理由即可)(1)某人早8时从山下旅店出发沿一条山路上山,下午5时到达山顶并留宿,次日8时沿同一山路下山,下午5时回到酒店。该人比在两天中的同一时刻经过路径中的同一地点,为什么?

答:题中人记作A,假定有另一个人B,第二天复制A第一天的上山过程(即任何时刻的足迹均相同),那么问题相当于A、B两人分别从一段路的两端同时出发,相向行走,同时到达另一端。则A、B二人必然会在路途中相遇,相遇点即为A在两天中的同一时刻经过路径中的同一地点。也可用下图加以说明,两条曲线分别为上山和下山的曲线,两天都是同一时刻出发,同一时刻到达,无论曲线如何变化,总会有一个相交点,相交点即为为A在两天中的同一时刻经过路径中的同一地点。

(2)甲乙两站之间有汽车相通,每隔10分钟甲乙两站相互发一趟车,但发车时刻不一定相同,甲乙两站之间有一中间站丙,某人每天在随机时刻到达丙站,并搭乘最先经过丙站的那趟车,结果发现100天中约有90天到达甲站,大约有10天到达乙站。问开往甲乙两站的汽车经过丙站的时刻是如何安排的?

答:由题意可知,搭乘乙站开往甲站的车概率是90%,搭乘甲站开往乙站的车概率是10%。说明90%的时间等来的是乙站开往甲站的车,10%的时间等来的是甲站开往乙站的车,而发车时间间隔是10分钟,进而说明:10分钟内有9分钟等来的是乙站开往甲站的车,有1分钟等来的是甲站开往乙站的车。因此可以得出结论:

乙站开往甲站的车到站后1分钟,甲站开往乙站的车到站。

开往甲乙两站的汽车经过丙站的时刻可表示如下(H0为乙站开往甲站的首班车时刻):

H0乙甲

H0+1 甲乙

H0+10 乙甲

H0+11 甲乙

H0+15 乙甲

H0+16 甲乙

…… ……

(3)张先生家住在A市,在B市工作,每天下班后他乘城际火车于18:00抵达A市火车站,他妻子驾车至火车站接他回家。一日他提前下班,乘早一班火车于17:30抵达A市火车站,随即步行回家,他妻子像往常一样驾车前来,在半路相遇将他接回家。到家时张先生发现比往常提前了10分钟,问张先生步行了多长时间?

答:张先生从火车站比往常早出发30分钟,早到家10分钟,则其步行路段用时比往常多了20分钟。他妻子从家出发时间不变,早到家10分钟,则其少行驶了10分钟的路程,即张先生步行路段的往返路程,那么张先生步行路段的单程行车时间为5分钟。张先生步行的时间应为比往常多用的时间和步行路段的单程行车时间之和,即25分钟。

(4)一男孩和一女孩分别在距家2公里和1公里且方向相反的两所学校上学,每天同时放学后分别以每小时4公里和每小时2公里的速度步行回家。一小狗以每小时6公里的速度由男孩处奔向女孩,又从女孩处奔向男孩,如此往返直至回到家中。问小狗奔波了多少路程。如果男孩和女孩上学时,小狗也往返奔波在他们中间,问他们到达学

校时小狗在何处?

答:由于他们到家的用时都是2÷4=0.5小时,所以小狗奔波的路程为6×0.5=3公里。上学过程是放学过程的逆向过程,所以在他们到达学校时小狗的位置为放学时的起始位置,即男孩所在的学校。

加分实验(公平投票问题)

某部门推出一专项基金目的在于培养优秀人才,根据评比结果确定资助的额度。许多单位的优秀者都申请了该基金,于是该基金的委员会聘请了数名专家,按照如下规则进行评比。

1.为了公平性,评委对本单位选手不给分;

2.每位评委对每位参与申请的人(除本单位选手外)都必须打分,且不打相同的分;

3.评委打分方法为给参加申请的人排序,根据优劣分别记1分、2分、…以此类推。

4.评判结束后,求出各选手的平均分,按平均分从低到高排序,依次确定本次评比的名次,即平均分最低者获得资助最高,依次类推。

本次基金申请中,甲所在单位有一名评委,这位评委将不参加对选手甲的评判,其他选手没有类似情况,评审结束后选手甲觉得这种评比规则对他不公平。问选手甲的抱怨是否有道理?若不公平,能否做出修正来解决选手甲的抱怨?

答:令评审个数为n,根据评分规则,在一般情况下评委给出的评分

=分别为k1,k2,k3…k n,所以该选手的平均得分为K1=k1+k2+?+k n

n

Sn

n

由于回避规则,甲的单位不得为甲评分,则甲的得分与一般选手

有差别,应为K2=k1+k2+?+k n?1

n?1=S n?k n

n?1

二者之差为

K1?K2=S n

n

?

S n?k n

n?1

=

S n(n?1)?n(S n?k n)

=nk n?S n

n(n?1)

由上式可知:只有当每个评委的打分均相等时,两者之差才恒等于0,其余情况则不恒等于0。根据题意,每个评委的对每个选手的打分是不等的,所以此规则对于选手甲确实有不公平之处。

若在尽可能公平的情况下制定得分规则,可判断出甲所得到的公平平均分应介于K1和K2之间。可以用K1和K2的几何平均值来定义,即K=√1?K2

(注:可编辑下载,若有不当之处,请指正,谢谢!)

数学建模实验报告

在下面的题目中选做100分的题目,给出详略得当的答案。 一.通过举例简要说明数学建模的一般过程或步骤。(15分) 答:建立数学模型的方法大致有两种,一种是实验归纳的方法,即根据测试或计算数据,按照一定的数据,按照一定的数学方法,归纳出系统的数学模型;另一种是理论分析的方法,具体步骤有五步(以人口模型 为例): 1、明确问题,提出合理简化的假设:首先要了解问题的实际背景,明确题目的要求,收集各种必要的信息 2、建立模型:据所做的假设以及事物之间的联系,构造各种量之间的关系。(查资料得出数学式子或算法)。 3、模型求解:利用数学方法来求解上一步所得到的数学问题,此时往往还要做出进一步的简化或假设。注意要尽量采用简单的数学公具。例如:马尔萨斯模型,洛杰斯蒂克模型 4、模型检验:根据预测与这些年来人口的调查得到的数目进行对比检验 5、模型的修正和最后应用:所建立的模型必须在实际应用中才能产生效益,根据预测模型,制定方针政策,以实现资源的合理利用和环境的保护。 二.把一张四条腿等长的正方形桌子放在稍微有些起伏的地面上,通常只有三只脚着地,然而 只需稍为转动一定角度,就可以使四只脚同时着地,即放稳了。(1) 请用数学模型来描述和证明这个实际问题; (2)讨论当桌子是长方形时,又该如何描述和证明?(15分) 答: 模型假设: 1.椅子四条腿一样长,椅脚与地面的接触部分相对椅子所占的地面面积可视为一个点。 2.地面凹突破面世连续变化的,沿任何方向都不会出现间断(没有向台阶那样的情况),即地面可看作数学上的连续曲面。 3.相对椅脚的间距和椅子腿的长度而言,地面是相对平坦的,即使椅子在任何位置至少有三条腿同时着地。4.椅子四脚连线所构成的四边形是圆内接四边形,即椅子四脚共圆。 5.挪动仅只是旋转。 我们将椅子这两对腿的交点作为坐标原点,建立坐标系,开始时AC、BD这两对腿都在坐标轴上。将AC和BD这两条腿逆时针旋转角度θ。记AC到地面的距离之和为f(θ)。记BD到 地面的距离之和为g(θ)。易得f(θ),g(θ)至少有一个为零。

《数学建模与数学实验》本科教学日历

《数学建模与数学实验》本科教学日历 数学建模部分 开设课程课程名称数学建模课程编号0701107 施教单位理学院 课内学时 总课时36 课程性质公共基础讲授课时28 修读要求选修实践课时8 选用教材教材名称数学建模教程出版社名称高等教育出版社 出版时间 及版次 2011年出版,第一版印刷时间2011年 其他情况 教学安排 班次授课对象及人数任教教员(指导教员)姓名及职称数学建模A 各专业本科学员 吴孟达教授 段晓君教授 毛紫阳讲师 王丹讲师 数学建模B 各专业本科学员 吴孟达教授 段晓君教授 毛紫阳讲师 王丹讲师 课次节 次 授课内容 教学 方法 采用现代化教学手段(课时) 多媒体电教双语网络实验 1 1 (1)什么是数学建模?数学建模的一般概念 (2)几个数学建模问题 讲授 1 2 (1)数学建模的一般步骤 (2)敏感问题调查案例 讲授 1 2 3 (1)行走步长问题 (2)雨中行走淋雨量最小问题 (3)道路是越多越通畅吗? 讲授 1 4 (1)有奖销售的抽奖策略问题 (2)“非诚勿扰”女生最佳选择问题 (3)网络文章流行度预测和招聘匹配 讲授 1 3 5 (1)线性规划模型基本概念 (2)整数规划模型 (3)0-1规划模型 讲授 1 6 (1)非线性规划 (2)多目标规划 讲授 1 4 7 (1)最短路算法 (2)最小生成树算法 讲授 1 8 (1)最大流算法 (2)PageRank算法 讲授 1 5 9 规划模型上机实践实践 1

课次节 次 授课内容 教学 方法 采用现代化教学手段(课时) 多媒体电教双语网络实验10 图论模型上机实践实践 1 6 11 (1)博弈模型基本概念 (2)Nash平衡和Pareto最优 (3)博弈论案例 讲授 1 12 (1)贝叶斯纳什均衡 (2)拍卖模型 讲授 1 7 13 社会选择理论中的选举问题数学模型-阿罗不可能定理讲授 1 14 越野长袍团体赛排名规则公平性问题讲授 1 8 15 军事作战模型-Lanchester作战模型讲授 1 16 自动化车床管理模型讲授 1 9 17 (1)“边际效应”基本概念 (2)实物交换模型,最佳消费模型、报童售报问题 讲授 1 18 (1)价格弹性模型 (2)合作效益的Shapley值分配模型 讲授 1 10 19 (1)聚类分析基本概念 (2)常用聚类算法 讲授 1 20 (1)方差分析基本概念 (2)单因素方差分析 (3)双因素方差分析 讲授 1 11 21 (1)主成分分析基本概念 (2)因子分析 讲授 1 22 (1)一元回归分析 (2)多元回归分析 (3)多元回归模型的检验与优化 讲授 1 12 23 聚类分析和方差分析上机实践实践 1 24 主成分分析和多元回归分析上机实践实践 1 13 25 (1)遗传算法基本思想 (2)算法步骤 讲授 1 26 遗传算法计算实例讲授 1 14 27 (1)模拟退火算法基本思想 (2)算法步骤 讲授 1 28 模拟退火算法计算实例讲授 1 15 29 (1)蚁群算法基本思想 (2)算法步骤 讲授 1 30 (1)数学建模中的计算机仿真 (2)不可召回的秘书招聘问题 (3)车灯光源优化设计 (4)生命游戏 讲授 1 16 31 遗传算法上机实践实践 1 32 模拟退火算法上机实践实践 1

数学建模与数学实验习题

数学建模与数学实验课程总结与练习内容总结 第一章 1.简述数学建模的一般步骤。 2.简述数学建模的分类方法。 3.简述数学模型与建模过程的特点。 第二章 4.抢渡长江模型的前3问。 5.补充的输油管道优化设计。 6.非线性方程(组)求近似根方法。 第三章 7.层次结构模型的构造。 8.成对比较矩阵的一致性分析。 第五章 9.曲线拟合法与最小二乘法。 10 分段插值法。 第六章 11 指数模型及LOGISTIC模型的求解与性质。 12.VOLTERRA模型在相平面上求解及周期平均值。 13 差分方程(组)的平衡点及稳定性。 14 一阶差分方程求解。 15 养老保险模型。

16 金融公司支付基金的流动。 17 LESLLIE 模型。 18 泛函极值的欧拉方法。 19 最短路问题的邻接矩阵。 20 最优化问题的一般数学描述。 21 马尔科夫过程的平衡点。 22 零件的预防性更换。 练习集锦 1. 在层次分析法建模中,我们介绍了成对比较矩阵概念,已知矩阵P 是成对比较矩阵 31/52a b P c d e f ?? ??=?????? ,(1)确定矩阵P 的未知元素。 (2)求 P 模最大特征值。 (3)分析矩阵P 的一致性是否可以接受(随机一致性指标RI取0.58)。 2. 在层次分析法建模中,我们介绍了成对比较矩阵概念,已知矩阵P 是三阶成对比较矩阵 322P ? ???=?????? ,(1)将矩阵P 元素补全。 (2)求P 模最 大特征值。 (3)分析矩阵P 的一致性是否可以接受。 3.考虑下表数据

(1)用曲改直的思想确定经验公式形式。 (2)用最小二乘法确定经验公式系数。 4.. 考虑微分方程 (0.2)0.0001(0.4)0.00001dx x xy dt dy y xy dt εε?=--????=-++?? (1)在像平面上解此微分方程组。(2)计算0ε=时的周期平均值。(3)计算0.1ε=时,y 的周期平均值占总量的周期平均值的比例增加了多少? 5考虑种群增长模型 '()(1/1000),(0)200x t kx x x =-= (1)求种群量增长最快的时刻。(2)根据下表数据估计参数k 值。 6. 布均匀,若环保部门及时发现并从某时刻起切断污染源,并更新湖水(此处更新指用新鲜水替换污染水),设湖水更新速率是 3 (m r s 单位:)。 (1) 试建立湖中污染物浓度随时间下降的数学模型? 求出污染物浓度降为控制前的5%所需要的时间。 7. 假如保险公司请你帮他们设计一个险种:35岁起保,每月交费400元,60岁开始领取养老金,每月养老金标准为3600元,请估算该保险费月利率为多少(保留到小数点后5位)? 8. 某校共有学生40000人,平时均在学生食堂就餐。该校共有,,A B C 3 个学生食堂。经过近一年的统计观测发现:A 食堂分别有10%,25%的学生经常去B ,C 食堂就餐,B 食堂经常分别有15%,25%的同学去

数学建模实验报告

数学建模实验报告

一、实验目的 1、通过具体的题目实例,使学生理解数学建模的基本思想和方法,掌握 数学建模分析和解决的基本过程。 2、培养学生主动探索、努力进取的的学风,增强学生的应用意识和创新 能力,为今后从事科研工作打下初步的基础。 二、实验题目 (一)题目一 1、题目:电梯问题有r个人在一楼进入电梯,楼上有n层。设每个 乘客在任何一层楼出电梯的概率相同,试建立一个概率模型,求直 到电梯中的乘客下完时,电梯需停次数的数学期望。 2、问题分析 (1)由于每位乘客在任何一层楼出电梯的概率相同,且各种可能的情况众多且复杂,难于推导。所以选择采用计算机模拟的 方法,求得近似结果。 (2)通过增加试验次数,使近似解越来越接近真实情况。 3、模型建立 建立一个n*r的二维随机矩阵,该矩阵每列元素中只有一个为1,其余都为0,这代表每个乘客在对应的楼层下电梯(因为每 个乘客只会在某一层下,故没列只有一个1)。而每行中1的个数 代表在该楼层下的乘客的人数。 再建立一个有n个元素的一位数组,数组中只有0和1,其中1代表该层有人下,0代表该层没人下。 例如: 给定n=8;r=6(楼8层,乘了6个人),则建立的二维随机矩阵及与之相关的应建立的一维数组为: m = 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 c = 1 1 0 1 0 1 1 1 4、解决方法(MATLAB程序代码):

n=10;r=10;d=1000; a=0; for l=1:d m=full(sparse(randint(1,r,[1,n]),1:r,1,n,r)); c=zeros(n,1); for i=1:n for j=1:r if m(i,j)==1 c(j)=1; break; end continue; end end s=0; for x=1:n if c(x)==1 s=s+1; end continue; end a=a+s; end a/d 5、实验结果 ans = 6.5150 那么,当楼高11层,乘坐10人时,电梯需停次数的数学期望为6.5150。 (二)题目二 1、问题:某厂生产甲乙两种口味的饮料,每百箱甲饮料需用原料6 千克,工人10名,可获利10万元;每百箱乙饮料需用原料5千 克,工人20名,可获利9万元.今工厂共有原料60千克,工人 150名,又由于其他条件所限甲饮料产量不超过8百箱.问如何 安排生产计划,即两种饮料各生产多少使获利最大.进一步讨 论: 1)若投资0.8万元可增加原料1千克,问应否作这项投资. 2)若每百箱甲饮料获利可增加1万元,问应否改变生产计划. 2、问题分析 (1)题目中共有3个约束条件,分别来自原料量、工人数与甲饮料产量的限制。 (2)目标函数是求获利最大时的生产分配,应用MATLAB时要转换

数学建模与数学实验课后习题答案

P59 4.学校共1002名学生,237人住在A 宿舍,333人住在B 宿舍,432人住在C 宿舍。学生要组织一个10人的委员会,使用Q 值法分配各宿舍的委员数。 解:设P 表示人数,N 表示要分配的总席位数。i 表示各个宿舍(分别取A,B,C ),i p 表示i 宿舍现有住宿人数,i n 表示i 宿舍分配到的委员席位。 首先,我们先按比例分配委员席位。 A 宿舍为:A n = 365.21002 10237=? B 宿舍为:B n =323.31002 10333=? C 宿舍为:C n =311.4100210432=? 现已分完9人,剩1人用Q 值法分配。 5.93613 22372 =?=A Q 7.92404 33332 =?=B Q 2.93315 44322 =?=C Q 经比较可得,最后一席位应分给A 宿舍。 所以,总的席位分配应为:A 宿舍3个席位,B 宿舍3个席位,C 宿舍4个席位。

商人们怎样安全过河

由上题可求:4个商人,4个随从安全过河的方案。 解:用最多乘两人的船,无法安全过河。所以需要改乘最多三人乘坐的船。 如图所示,图中实线表示为从开始的岸边到河对岸,虚线表示从河对岸回来。商人只需要按照图中的步骤走,即可安全渡河。总共需要9步。

P60 液体在水平等直径的管内流动,设两点的压强差ΔP 与下列变量有关:管径d,ρ,v,l,μ,管壁粗糙度Δ,试求ΔP 的表达式 解:物理量之间的关系写为为()?=?,,,,,μρ?l v d p 。 各个物理量的量纲分别为 []32-=?MT L p ,[]L d =,[]M L 3-=ρ,[]1-=LT v ,[]L l =,[]11--=MT L μ,Δ是一个无量纲量。 ???? ??????-----=?0310100011110010021113173A 其中0=Ay 解得 ()T y 00012111---=, ()T y 00101102--=, ()T y 01003103--=, ()T y 10000004= 所以 l v d 2111---=ρπ,μρπ112--=v ,p v ?=--313ρπ,?=4π 因为()0,,,,,,=??p l v d f μρ与()0,,,4321=ππππF 是等价的,所以ΔP 的表达式为: ()213,ππψρv p =?

数学建模实验报告

数学建模实验报告 实验一计算课本251页A矩阵的最大特征根和最大特征向量 1 实验目的 通过Wolfram Mathematica软件计算下列A矩阵的最大特征根和最大特征向量。 2 实验过程 本实验运用了Wolfram Mathematica软件计算,计算的代码如下:

3 实验结果分析 从代码的运行结果,可以得到最大特征根为5.07293,最大特征向量为 {{0.262281},{0.474395},{0.0544921},{0.0985336},{0.110298}},实验结果 与标准答案符合。

实验二求解食饵-捕食者模型方程的数值解 1实验目的 通过Wolfram Mathematica或MATLAB软件求解下列习题。 一个生物系统中有食饵和捕食者两种种群,设食饵的数量为x(t),捕食者为y(t),它们满足的方程组为x’(t)=(r-ay)x,y’(t)=-(d-bx)y,称该系统为食饵-捕食者模型。当r=1,d=0.5,a=0.1,b=0.02时,求满足初始条件x(0)=25,y(0)=2的方程的数值解。 2 实验过程 实验的代码如下 Wolfram Mathematica源代码: Clear[x,y] sol=NDSolve[{x'[t] (1-0.1y[t])x[t],y'[t] 0.02x[t]y[t]-0.5y[t],x[0 ] 25,y[0] 2},{x[t],y[t]},{t,0,100}] x[t_]=x[t]/.sol y[t_]=y[t]/.sol g1=Plot[x[t],{t,0,20},PlotStyle->RGBColor[1,0,0],PlotRange->{0,11 0}] g2=Plot[y[t],{t,0,20},PlotStyle->RGBColor[0,1,0],PlotRange->{0,40 }] g3=Plot[{x[t],y[t]},{t,0,20},PlotStyle→{RGBColor[1,0,0],RGBColor[ 0,1,0]},PlotRange->{0,110}] matlab源代码 function [ t,x ]=f ts=0:0.1:15; x0=[25,2]; [t,x]=ode45('shier',ts,x0); End function xdot=shier(t,x)

数学建模与数学实验试卷及答案

数学建模与数学实验试卷及答案 二、本题10分(写出程序和结果) 蚌埠学院2010—2011学年第二学期 2,x在 [-5 ,5] 区间内的最小值,并作图加以验证。求函数yxe,,,3《数学建模与数学实验》补考试卷答案 f1=inline('x.^2 +exp(-x)-3') 注意事项:1、适用班级:09数学与应用数学本科1,2班 2、本试卷共1页,附答题纸1页。满分100分。 x=fmin(f1,-5,5) 3、考查时间100分钟。 y=f1(x) 4、考查方式:开卷 fplot(f1,[-5,5]) 一、填空:(每空4分,共60分) x = 0.3517,y== -2.1728 123111,,,,, ,,,,三、本题15分(写出程序和结果) 1. 已知,,则A的秩为 3 ,A的特征值为 A,612B,234,,,, ,,,,,215531,,,,,360000xx,,,12,max2.5fxx,,求解:, stxx..250000,,,1212-1.9766 4.4883 + 0.7734i 4.4883 - 0.7734i ,若令 A([1,3],:)= B([2,3],:),则,x,150001,A(2,:)= 6 1 2 ; 解: xxx,,,22,123,model: 2. 的解为 1.25 ,0.25 0.5 ; xxx,,,521,123max=2.5*x1+x2; ,242xxx,,,123,3*x1+x2<=60000; 装订线内不要答题 2*x1+x2<=50000; 3. 将1234521 分解成质因数乘积的命令为_factor(sym(‘1234521’)),

数学建模实验

数学建模课程实验报告 专题实验7 班级数财系1班学号2011040123 丛文 实验题目常微分方程数值解 实验目的 1.掌握用MATLAB求微分方程初值问题数值解的方法; 2.通过实例学习微分方程模型解决简化的实际问题; 3.了解欧拉方法和龙格库塔方法的基本思想。 实验容 (包括分 析过程、 方法、和 代码,结 果) 1. 用欧拉方法和龙格库塔方法求下列微分方程初值问题的数值 解,画出解的图形,对结果进行分析比较 解;M文件 function f=f(x,y) f=y+2*x; 程序; clc;clear; a=0;b=1; %求解区间 [x1,y_r]=ode45('f',[a b],1); %调用龙格库塔求解函数求解数值 解; %% 以下利用Euler方法求解 y(1)=1;N=100;h=(b-a)/N; x=a:h:b;

for i=1:N y(i+1)=y(i)+h*f(x(i),y(i)); end figure(1) plot(x1,y_r,'r*',x,y,'b+',x,3*exp(x)-2*x-2,'k-');%数值解与真解图 title('数值解与真解图'); legend('RK4','Euler','真解'); xlabel('x');ylabel('y'); figure(2)

plot(x1,abs(y_r-(3*exp(x1)-2*x1-2)),'k-');%龙格库塔方法的误差 title('龙格库塔方法的误差') xlabel('x');ylabel('Error'); figure(3) plot(x,abs(y-(3*exp(x)-2*x-2)),'r-')%Euler方法的误差 title('Euler方法的误差') xlabel('x');ylabel('Error');

数学建模实验报告第十一章最短路问答

实验名称:第十一章最短路问题 一、实验内容与要求 掌握Dijkstra算法和Floyd算法,并运用这两种算法求一些最短路径的问题。 二、实验软件 MATLAB7.0 三、实验内容 1、在一个城市交通系统中取出一段如图所示,其入口为顶点v1,出口为顶点v8,每条弧段旁的数字表示通过该路段所需时间,每次转弯需要附加时间为3,求v1到v8的最短时间路径。 V1 1 V2 3 V3 1 V5 6 V6 V4 2 V7 4 V8

程序: function y=bijiaodaxiao(f1,f2,f3,f4) v12=1;v23=3;v24=2;v35=1;v47=2;v57=2;v56=6;v68=3;v78=4; turn=3; f1=v12+v23+v35+v56+turn+v68; f2=v12+v23+v35+turn+v57+turn+v78; f3=v12+turn+v24+turn+v47+v78; f4=v12+turn+v24+v47+turn+v57+turn+v56+turn+v68; min=f1; if f2

f4 实验结果: v1到v8的最短时间路径为15,路径为1-2-4-7-8. 2、求如图所示中每一结点到其他结点的最短路。V110 V3V59 V6

floy.m中的程序: function[D,R]=floyd(a) n=size(a,1); D=a for i=1:n for j=1:n R(i,j)=j; end end R for k=1:n for i=1:n for j=1:n if D(i,k)+D(k,j)

数学建模实验报告(1)

四川师范大学数学与软件科学学院 实验报告 课程名称:数学建模 指导教师:陈东 班级:_2008级2班_____________ 学号:__2008060244___________ 姓名:___邢颖________ 总成绩:______________

数学与软件科学学院 实验报告 学期:_2009__ 年至2010 _年____ 第_ 二___ 学期 2010 年 4 月 1 _日 课程名称:_数学建模__ 专业:数学与应用数学____ 2008__ _级_ 2 ___班 实验编号: 1 实验项目_Matlab 入门_ 指导教师 陈东 姓名: 邢颖 ____ 学号: 2008060244 一、实验目的及要求 实验目的: 实验要求: 二、实验内容 (1)用起泡法对10个数由小到大排序. 即将相邻两个数比较,将小的调到前头. (2)有一个 4*5 矩阵,编程求出其最大值及其所处的位置. (3)编程求 (4)一球从100米高度自由落下,每次落地后反跳回原高度的一半,再落下. 求它在第10次落地时,共经过多少米?第10次反弹有多高? (5)有一函数 ,写一程序,输入自变量的值,输出函数值. 三、实验步骤(该部分不够填写.请填写附页) (2) x=[1 6 2 7 6;4 6 1 3 2;1 2 3 4 7;8 1 4 6 3]; t=x(1,1); for i=1:4 for j=1:5 if x(i,j)>t t=x(i,j); a=[i,j]; end ∑=20 1! n n y xy x y x f 2sin ),(2 ++=

end end (3)程序1: x(1)=1; s=1; for n=2:20 x(n)=x(n-1)*n; s=s+x(n); end s 程序2; s=0,m=1; for n=2:20; m=m*n; s=s+m; end s 结果:s = 2.5613e+018 (4)程序 s=100 h=s/2 for n=2:10 s=s+2*h h=h/2 end s,h 结果:s = 299.6094 h = 0.0977 (5)程序: function f=fun1(x,y) f=x^2+sin(x*y)+2*y

数学建模与数学实验报告

数学建模与数学实验报告 指导教师__郑克龙___ 成绩____________ 组员1:班级______________ 姓名______________ 学号_____________ 组员2:班级______________ 姓名______________ 学号______________ 实验1.(1)绘制函数cos(tan())y x π=的图像,将其程序及图形粘贴在此。 >> x=-pi:0.01:pi; >> y=cos(tan(pi*x)); >> plot(x,y) -4 -3 -2 -1 1 2 3 4 -1-0.8-0.6-0.4-0.200.20.40.60.8 1 (2)用surf,mesh 命令绘制曲面2 2 2z x y =+,将其程序及图形粘贴在此。(注:图形注意拖放,不要太大)(20分) >> [x,y]=meshgrid([-2:0.1:2]); >> z=2*x.^2+y.^2; >> surf(x,y,z)

-2 2 >> mesh(x,y,z) -2 2 实验2. 1、某校60名学生的一次考试成绩如下:

93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55 1)计算均值、标准差、极差、偏度、峰度,画出直方图;2)检验分布的正态性;3)若检验符合正态分布,估计正态分布的参数并检验参数. (20分) 1) >> a=[93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55]; >> pjz=mean(a) pjz = 80.1000 >> bzhc=std(a) bzhc = 9.7106 >> jc=max(a)-min(a) jc = 44 >> bar(a)

《数学建模与数学实验》课程论文

10级信息《数学建模与数学实验(实践)》任务书 一、设计目的 通过《数学建模与数学实验(实践)》实践环节,掌握本门课程的众多数学建模方法和原理,并通过编写C语言或matlab程序,掌握各种基本算法在计算机中的具体表达方法,并逐一了解它们的优劣、稳定性以及收敛性。在熟练掌握C 语言或matlab语言编程的基础上,编写算法和稳定性均佳、通用性强、可读性好,输入输出方便的程序,以解决实际中的一些科学计算问题。 二、设计教学内容 1线性规划(掌握线性规划的模型、算法以及Matlab 实现)。整数线性规划(掌握整数线性规划形式和解法)。 2微分方程建模(掌握根据规律建立微分方程模型及解法;微分方程模型的Matlab 实现)。 3最短路问题(掌握最短路问题及算法,了解利用最短路问题解决实际问题)。 行遍性问题(了解行遍性问题,掌握其TSP算法)。 4回归分析(掌握一元线性回归和多元线性回归,掌握回归的Matlab实现)。 5计算机模拟(掌握Monte-carlo方法、了解随机数的产生;能够用Monte-carlo 解决实际问题)。 6插值与拟合(了解数据拟合基本原理,掌握用利用Matlab工具箱解决曲线拟合问题)。 三、设计时间 2012—2013学年第1学期:第16周共计一周 目录 一、10级信息《数学建模与数学实验(实践)》任务书 (1) 二、饭店餐桌的布局问题 (3) 摘要 (3)

问题重述 (3) 模型假设 (3) 模型分析 (4) 模型的建立和求解 (4) 模型推广 (9) 参考文献 (9) 三、白酒配比销售问题 (10) 摘要 (10) 问题重述 (11) 问题分析 (12) 模型假设 (12) 符号及变量说明 (12) 模型的建立与求解 (13) 模型的检验 (18) 模型的评价与推广 (19) 附录 (21) 饭店餐桌的布局问题 摘要 饭店餐桌的布局对于一个饭店有着很重要的作用。本文讨论的就是饭店餐桌的布局问题,根据实际需求及规定建立模型,同时考虑餐桌的类型及规格,尤其是餐桌的摆放技巧,保证使饭店能容纳的人数达到最大。根据所需餐桌的数量

数学建模迭代实验报告(新)

非 线 性 迭 代 实 验 报 告 一、实验背景与实验目的 迭代是数学研究中的一个非常重要的工具,通过函数或向量函数由初始结点生成迭代结点列,也可通过函数或向量函数由初值(向量)生成迭代数列或向量列。 蛛网图也是一个有用的数学工具,可以帮助理解通过一元函数由初值生成的迭代数列的敛散性,也帮助理解平衡点(两平面曲线交点)的稳定性。 本实验在Mathematica 平台上首先利用蛛网图和迭代数列研究不动点的类型;其次通过蛛网图和迭代数列研究Logistic 映射,探索周期点的性质、认识混沌现象;第三通过迭代数列或向量列求解方程(组)而寻求有效的求解方法;最后,利用结点迭代探索分形的性质。 二、实验材料 2.1迭代序列与不动点 给定实数域上光滑的实值函数)(x f 以及初值0x ,定义数列 )(1n n x f x =+, ,2,1,0=n (2.2.1) }{n x 称为)(x f 的一个迭代序列。 函数的迭代是数学研究中的一个非常重要的思想工具,利用迭代序列可以研究函数)(x f 的不动点。 对函数的迭代过程,我们可以用几何图象来直观地显示它——“蜘蛛网”。运行下列Mathematica 程序: Clear[f] f[x_] := (25*x - 85)/(x + 3); (实验时需改变函数) Solve[f[x]==x , x] (求出函数的不动点) g1=Plot[f[x], {x, -10, 20}, PlotStyle -> RGBColor[1, 0, 0], DisplayFunction -> Identity]; g2=Plot[x, {x, -10, 10}, PlotStyle -> RGBColor[0, 1, 0], DisplayFunction -> Identity]; x0=5.5; r = {}; r0=Graphics[{RGBColor[0, 0, 1], Line[{{x0, 0}, {x0, x0}}]}]; For[i = 1, i <= 100, i++, r=Append[r, Graphics[{RGBColor[0, 0, 1], Line[{{x0, x0}, {x0, f[x0]}, {f[x0], f[x0]}}] }]]; x0=f[x0] ]; Show[g1, g2, r, r0, PlotRange -> {-1, 20}, (PlotRange 控制图形上下范围) DisplayFunction -> $DisplayFunction] x[0]=x0; x[i_]:=f[x[i-1]]; (定义序列) t=Table[x[i],{i,1,10}]//N ListPlot[t] (散点图) 观察蜘蛛网通过改变初值,你能得出什么结论? 如果只需迭代n 次产生相应的序列,用下列Mathematica 程序: Iterate[f_,x0_,n_Integer]:= Module[{ t={},temp= x0},AppendTo[t,temp]; For[i=1,i <= n, i++,temp= f[temp]; AppendTo[t,temp]]; t ] f[x_]:= (x+ 2/x)/2; Iterate[f,0.7,10]

数模实验报告

数学建模与实验实验报告 姓名:李明波 院系:仪器科学与工程学院 学号:22013108 老师:王峰

数学建模与实验实验报告 实验一 实验题目 (1)已知某平原地区的一条公路经过如下坐标所示的点,请采用样条插值绘出这条公路(不考虑 (2)对于上表给出的数据,估计公路长度。 实验过程 (1)第一问代码如下: X=[0,30,50,70,80,90,120,148,170,180,202,212,230,248,268,271,280,290,300,312,320,340,3 60,372,382,390,416,430,478]; Y=[80,64,47,42,48,66,80,120,121,138,160,182,200,208,212,210,200,196,188,186,200,184,1 88,200,202,240,246,280,296]; %给出坐标点 xx=0:1:478;%选取0~478内的点 yy=spline(X,Y,xx);%样条插值法找出曲线 plot(X,Y, 'p ',xx,yy, 'g ');%绘出曲线图 x=[440,420,380,360,340,320,314,280,240,200]; y=[308,334,328,334,346,356,360,392,390,400]; hold on xy=440:-1:200; yx=spline(x,y,xy); plot(x,y, 'p ',xy,yx, 'g '); 运行上述代码得到结果如下:

上图为所绘公路图 (2)代码如下: X=[0 30 50 70 80 90 120 148 170 180 202 212 230 248 268 271 280 290 300 312 320 340 360 372 382 390 416 430 478 440 420 380 360 340 320 314 280 240 200]; Y=[80 64 47 42 48 66 80 120 121 138 160 182 200 208 212 210 200 196 188 186 200 184 188 200 202 240 246 280 296 308 334 328 334 346 356 360 392 390 400]; for k=1:length(X)-1 len(k)=sqrt((X(k+1)-X(k))^2+(Y(k+1)-Y(k))^2); end; Len=sum(len);Len 运行得到结果如下: 即公路长为967.46米。

数学建模与数学实验

数学建模与数学实验 实验报告 班级: 数学师范153 姓名:付爽 学号:1502012060 实验名称: 数列极限与函数极限 基础实验 基础实验一数列极限与函数极限第一部分实验指导书解读

一、实验目的 从刘徽的割圆术、裴波那奇数列研究数列的收敛性并抽象出极限的定义;理解数列收敛的准则;理解函数极限与数列极限的关系。 二、实验使用软件 Mathematic 5、0 三.实验的基本理论即方法 1割圆术 中国古代数学家刘徽在《九章算术注》方田章圆田术中创造了割圆术计算圆周率π。刘徽先注意到圆内接正多边形的面积小于圆面积;其次,当将边数屡次加倍时,正多边形的面积增大,边数愈大则正多边形面积愈近于圆的面积。 “割之弥细,所失弥少。割之又割以至不可割,则与圆合体而无所失矣。”这几句话明确地表明了刘徽的极限思想。 以n S 表示单位圆的圆内接正1 23-?n 多边形面积,则其极限为 圆周率π。用下列Mathematica 程序可以从量与形两个角度考察数列{n S }的收敛情况: m=2;n=15;k=10; For[i=2,i<=n,i++, l[i_]:=N[2*Sin[Pi/(3*2^i)],k]; (圆

内接正1 23-?n 多边形边长) s[i_]:=N[3*2^(i-1)*l[i]*Sqrt[1-(l[i])^2/4],k]; (圆内接正1 23-?n 多边形面积) r[i_]:=Pi-s[i]; d[i_]:=s[i]-s[i-1]; Print[i," ",r[i]," ",l[i]," ",s[i]," ",d[i]] ] t=Table[{i,s[i]},{i,m,n}] (数组) ListPlot[t] (散点图) 2裴波那奇数列与黄金分割 由2110;1; 0--+===n n n F F F F F 有著名的裴波那奇数列}{n F 。 如果令n n n F F R 11 --=,由n F 递推公式可得出 11111/11---+=+=+=n n n n n n n R F F F F F R ,]251251[511 1 ++??? ? ??--??? ? ??+=n n n F ; 2 15lim lim 1 -==+∞ →∞ →n n n n n F F R 。 用下列Mathematica 程序可以从量与形两个角度考察数列{n R }的收敛情况: n=14,k=10; For[i=3,i<=n,i++, t1=(Sqrt[5]+1)/2; t2=(1-Sqrt[5])/2;

数模模数转换实验报告材料

数模模数转换实验报告 一、实验目的 1、了解数模和模数转换电路的接口方法及相应程序设计方法。 2、了解数模和模数转换电路芯片的性能和工作时序。 二、实验条件 1、DOS操作系统平台 2、数模转换芯片DAC0832和模数转换器ADC0809芯片。 三、实验原理 1、数模转换: (1)微机处理的数据都是数字信号,而实际的执行电路很多都是模拟的。因此微机的处理结果又常常需要转换为模拟信号去驱动相应的执行单元,实现对被控对象的控制。这种把数字量转换为模拟量的设备称为数模转换器(DAC),简称D/A。 (2)实验中所用的数模转换芯片是DAC0832,它是由输入寄存器、DAC 寄存器和D/A 转换器组成的CMOS 器件。其特点是片包含两个独立的8 位寄存器,因而具有二次缓冲功能,可以将被转换的数据预先存在DAC 寄存器中,同时又采集下一组数据,这就可以根据需要快速修改DAC0832 的输出。 2、模数转换: (1)在工程实时控制中,经常要把检测到的连续变化的模拟信号,如温度、压力、速度等转换为离散的数字量,才能输入计算机进行处理。实现模拟量到数字量转换的设备就是模数转换器(ADC),简称A/D。

(2)模数转换芯片的工作过程大体分为三个阶段:首先要启动模数转换过程。其次,由于转换过程需要时间,不能立即得到结果,所以需要等待一段时间。一般模数转换芯片会有一条专门的信号线表示转换是否结束。微机可以将这条信号线作为中断请求信号,用中断的方式得到转换结束的消息,也可以对这条信号线进行查询,还可以采用固定延时进行等待(因为这类芯片转换时间是固定的,事先可以知道)。最后,当判断转换已经结束的时候,微机就可以从模数转换芯片中读出转换结果。 (3)实验采用的是8 路8 位模数转换器ADC0809 芯片。ADC0809 采用逐次比较的方式进行A/D 转换,其主要原理为:将一待转换的模拟信号与一个推测信号进行比较,根据推测信号是大于还是小于输入信号来决定增大还是减少该推测信号,以便向模拟输入逼近。推测信号由D/A 转换器的输出获得,当推测信号与模拟信号相等时,向D/A 转换器输入的数字就是对应模拟信号的数字量。ADC0809 的转换时间为64 个时钟周期(时钟频率500K 时为128S)。分辨率为 8 位,转换精度为±LSB/2,单电源+5V 供电时输入模拟电压围为04.98V。 四、实验容 1、把DAC0832 的片选接偏移为10H 的地址,使用debug 命令来测试 DAC0832 的输出,通过设置不同的输出值,使用万用表测量Ua 和Ub 的模拟电压,检验DAC0832 的功能。选取典型(最低、最高和半量程等)的二进制值进行检验,记录测得的结果。实验结果记录如下: 输入 00 0.001 4.959 08 0.145 4.636

《数学建模与数学实验》上机实验报告

成都信息工程大学 《数学建模与数学实验》上机实验报告 专业信息与计算科学班级姓名学号 实验日期成绩等级教师评阅日期 [问题描述] 下表给出了某一海域以码为单位的直角坐标Oxy 上一点(x,y)(水面一点)以英尺为单位的水深z,水深数据是在低潮时测得的,船的吃水深为5英尺,问在矩形区域(75,200)x (-50,150)里那些地方船要避免进入。 [模型] 设水面一点的坐标为(x,y,z),用基点和插值函数在矩形区域(75,200)*(-50,150)内做二维插值、三次插值,然后在作出等高线图。

[求解方法] 使用matlab求解: M文件:water.m x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]; z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; cx = 75:0.5:200; cy = -50:0.5:150; [cx,cy]=meshgrid(cx,cy); 作出曲面图: 代码如下: >> water >> cz=griddata(x,y,z,cx,cy,'cubic'); >> meshz(cx,cy,cz) >> xlabel('X'),ylabel('Y'),zlabel('Z') >> 作出等高线图: 代码如下: >> water >> cz=griddata(x,y,z,cx,cy,'cubic'); >> figure(2) >> contour(cx,cy,cz,[-5,-5],'r') >> hold on >> plot(x,y,'*') >> xlabel('X'),ylabel('Y') [结果]

数学建模-实验报告11

《数学建模实验》实验报告 学号:______ 姓名: 实验十一:微分方程建模2 一只小船渡过宽为d的河流,目标是起点A 正对着的另一岸B点,已知河水流速w 与船在静水中的速度V2之比为k. 1?建立小船航线的方程,求其解析解; 2. 设d=100m,v i=1m/s,v2=2m/s,用数值解法求渡河所需时间、任意时刻小船的位置及航行曲线,作图,并与解析解比较。 一、问题重述 我们建立数学模型的任务有: 1. 由已给定的船速、水速以及河宽求出渡河的轨迹方程; 2. 已知船速、水速、河宽,求在任意时刻船的位置以及渡船所需要的时间。 二、问题分析 此题是一道小船渡河物理应用题,为典型的常微分方程模型,问题中船速、水速、河宽已经给定,由速度、时间、位移的关系,我们容易得到小船的轨迹方程,同时小船的起点和终点已经确定,给我们的常微分方程模型提供了初始条件。 三、模型假设 1?假设小船与河水的速度恒为定值v「V2 ,不考虑人为因素及各种自然原因; 2. 小船行驶的路线为连续曲线,起点为A,终点为B ; 3. 船在行驶过程中始终向着B点前进,即船速v2始终指向B ; 4. 该段河流为理想直段,水速w与河岸始终保持平行。 四、模型建立 y | B A 兀、 % \ * r v A X 如图,以A为原点,以沿河岸向右方向为x轴正向,以垂直河岸到B端方向为y轴正向建立平面直角坐标系。其中河水流速为v i,小船速度为V2,且w:v2 k,合速度为v,河宽为d,为72与直线AB的夹角。

V x V y 在t 时刻, 船 dx dt V i 小船在x 轴方向的位移为 x v 2 sin v 2 cos V i V 2 0,x(0) 0, y(0) ;(d y) 0. \ (d y) d y ______ 2 2 ' x dy v 2 cos 由(2)/(1)得到dx y(0) v-1 v 2 sin 0. dx In (2) (i )题 dx 对上式求倒数得 dx dy x ,在y 轴方向上的位移为y ,则t 时刻, 方向 的速度 模型求解 v 2 sin V 1 v 2 co s —, 则上式可化为 dx d y dy d ?dp pdy ydp ,代入上式, k J p 2 整理,得 P 2 | ln| d Cy | 也就是 x 2 (d y )2 y P (d y ) dp P 2 kdy ,积分可得 y C k ( ------- )k ,代入 d y x d y d y 2 0, y 0 d k (d y )k (d y )k d k (见附 录) ,对该情况下的微分方程的数值解进行分 60.0000 6.5451 98.2803 60.1000 6.4519 98.3319 60.2000 6.3585 98.3827 60.3000 6.2649 98.4327 60.4000 6.1711 98.4819 60.5000 6.0771 98.5304 60.6000 5.9829 98.5782 60.7000 5.8886 98.6251 60.8000 5.7940 98.6713 60.9000 5.6993 98.7168 61.0000 5.6043 98.7615 61.1000 5.5092 98.8054 题 由初始条件,设计程序 析,结果如下(省略了前60s 的数据):

相关文档
最新文档