数学建模1
数学建模活动(1)-练习题

课后练习
请同学们仿照上述过程开展一次建立模型解决实际问题的活动,可以继续研究不同室温下泡制一杯最佳口感茶水所需的时间,也可以从下列选题中选择一个:
1、应在炒菜之前多长时间将冰箱里的肉拿出来解冻?
2、根据某一同学的身高和体重,判断该同学是否超重.
3、用微波炉或电磁炉烧一壶开水,找到最省电的功率设定方法.
4估计阅读一本书所需要的时间.
也可以根据自己的兴趣成3—5人的研究小组,每位同学参加其中一个小组。在小组内,要确定一个课题负责人,使每位成员都有明确的分工.拟定研究课题、确定研究方案、规划研究步骤、编制研究手册,然后在班里进行一次开题报告。
数学建模基础练习一及参考答案

数学建模基础练习一及参考答案数学建模基础练习一及参考答案练习1matlab练习一、矩阵及数组操作:1.利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵([-1,1]之间)、正态分布矩阵(均值为1,方差为4),然后将正态分布矩阵中大于1的元素变为1,将小于1的元素变为0。
2.利用fix及rand函数生成[0,10]上的均匀分布的10×10的整数随机矩阵a,然后统计a中大于等于5的元素个数。
3.在给定的矩阵中删除含有整行内容全为0的行,删除整列内容全为0的列。
4.随机生成10阶的矩阵,要求元素值介于0~1000之间,并统计元素中奇数的个数、素数的个数。
二、绘图:5.在同一图形窗口画出下列两条曲线图像,要求改变线型和标记:y1=2x+5;y2=x^2-3x+1,并且用legend标注。
6.画出下列函数的曲面及等高线:z=sinxcosyexp(-sqrt(x^2+y^2)).7.在同一个图形中绘制一行三列的子图,分别画出向量x=[158101253]的三维饼图、柱状图、条形图。
三、程序设计:8.编写程序计算(x在[-8,8],间隔0.5)先新建的,在那上输好,保存,在命令窗口代数;9.用两种方法求数列:前15项的和。
10.编写程序产生20个两位随机整数,输出其中小于平均数的偶数。
11.试找出100以内的所有素数。
12.当时,四、数据处理与拟合初步:13.随机产生由10个两位随机数的行向量A,将A中元素按降序排列为B,再将B重排为A。
14.通过测量得到一组数据:t12345678910y4.8424.3623.7543.3683.1693.0383.0343.0163.0123.005分别采用y=c1+c2e^(-t)和y=d1+d2te^(-t)进行拟合,并画出散点及两条拟合曲线对比拟合效果。
15.计算下列定积分:16.(1)微分方程组当t=0时,x1(0)=1,x2(0)=-0.5,求微分方程t在[0,25]上的解,并画出相空间轨道图像。
数学建模1例题解析

1.贷款问题小王夫妇计划贷款20万元购买一套房子,他们打算用20年的时间还清贷款。
目前,银行的利率是%/月。
他们采用等额还款的方式(即每月的还款额相同)偿还贷款。
(1)在上述条件下,小王夫妇每月的还款额是多少?共计付了多少利息?(2)在贷款满5年后,他们认为他们有经济能力还完余下的款额,打算提前还贷,那么他们在第6年初,应一次付给银行多少钱,才能将余下全部的贷款还清?(3)如果在第6年初,银行的贷款利率由%/月调到%/月,他们仍然采用等额还款的方式,在余下的15年内将贷款还清,那么在第6年后,每月的还款额应是多少?(4)某借贷公司的广告称,对于贷款期在20年以上的客户,他们帮你提前三年还清贷款。
但条件是:(i)每半个月付款一次,但付款额不增加,即一次付款额是原付给银行还款额的1/2;(ii)因为增加必要的档案、文书等管理工作,因此要预付给借贷公司贷款总额10%的佣金。
试分析,小王夫妇是否要请这家借贷公司帮助还款。
解答:(1)贷款总月数为N=20*12=240,第240个月的欠款额为0,即。
利用式子(元),即每个月还款元,共还款(元),共计付利息元。
(2)贷款5年(即5*12=60个月)后的欠款额为,利用公式:,所以,(元)(3)元,即第六年初,贷款利率,所以余下的15年,每个月还款额为:(元)(4)按照借贷公司的条件(i)每半个月付款一次,但付款额不增加,即一次付款额是原付给银行还款额的,付款的时间缩短,但是前17年的付款总额不变。
帮忙提前三年还清需要资金数:。
对于条件(ii)佣金数:分析:因为预付佣金20000元,按照银行存款利率/月,17年的存款本息为即在第17年需要给付借贷公司的钱少于给付银行的钱。
所以建议请这家借贷公司帮助还款。
2.冷却定律与破案按照Newton冷却定律,温度为T的物体在温度为的环境中冷却的速度与温差成正比。
用此定律建立相应的微分方程模型。
凌晨某地发生一起凶杀案,警方于晨6时到达案发现场,测得尸温26℃,室温10℃,晨8时又测得尸温18℃。
数学建模作业(1)

数学建模作业(1)
数模
数模
1.学校共学校共1000名学生,235人住在宿名学生,人住在A宿名学生人住在人住B宿舍人住在C宿舍舍,333人住宿舍,432人住在宿舍人住宿舍,人住在宿舍.学生们要组织一个10人的委员会人的委员会,学生们要组织一个人的委员会,试用下列办法分配各宿舍的委员数:列办法分配各宿舍的委员数:(1)按比例分配取整数的名额后,剩下的名按比例分配取整数的名额后,按比例分配取整数的名额后额按惯例分给小数部分较大者。
额按惯例分给小数部分较大者。
(2)用Q值方法。
值方法。
用值方法
数模
如果委员会从10人增至人如果委员会从人增至15人,用以上人增至2种方法再分配名额。
将2种方法两次分配种方法再分配名额。
种方法再分配名额种方法两次分配的结果列表比较。
的结果列表比较。
(3)你能提出其它的方法吗?用你的方你能提出其它的方法吗?你能提出其它的方法吗法分配上面的名额。
法分配上面的名额。
数模
2.考察模拟水下爆炸的比例模型.爆炸物质量m,在距爆炸点距离r处设置仪器,接收到的冲击波压强为p,记大气初始压强p0,水的密度ρ,水的体积弹性模量k,用量纲分析法已经得到
p0ρrp=p0(,)km3
数模
设模拟实验与现场的p0,ρ,k相同,而爆炸物模型的质量为原模型的1/1000.为了使实验中接收到与现场相同的压强p,问实验时应如何设置接收冲击波的仪器,即求实验仪器与爆炸点之间的距离是现场的多少倍?
p0,ρ,k。
数学建模问题1

在习题1-8中,情景是模糊地陈述的。
从这些模糊的情景中,识别要研究的问题。
哪些变量影响到问题识别中你已经识别的行为?哪些变量最重要?记住,实际上没有正确的答案.1.单种群的总量增长.2.一家零售店要建造一个新的停车场,停车场应该怎样照明?3.一位农民期望他的地里种植的粮食农作物的产量达到最大,他正确地识别了问题吗?试讨论另一种目标.4.怎样设计一个供大班级用的演讲厅?5.一个物体从很高的地方掉下来.何时它撞击到地面?撞击到地面的力度有多大?6.某种产品的制造商应该怎样决定每年应该生产多少件产品,以及每件产品应该标价多少?7.美国食品及药物管理局(FDA)想要了解一种新药对控制人口中的某种疾病是否有效.8.滑雪者滑下山坡有多快?对于习题9~17中提出的情景,识别值得研究的问题并列出会影响你已经识别的行为的变量.哪些变量可以完全忽略?哪些变量在开始时可以认为它们是常数?你能识别出你想仔细研究的子模型吗?识别任何你想收集的数据.9.一位植物学家有兴趣研究叶子的形状以及影响叶子长成这种形状的各种支配力量,她从一棵白橡树的底部剪下几片叶子,发现叶子相当宽,没有很明显的锯齿形.当她到树的顶部去看时,她发现有很明显的锯齿形而几乎没有展得很宽的叶子.10. 不同大小的动物其他特性也不同.小动物比之于较大的动物,叫声尖细、心跳较快以及呼吸次数更多.另一方面,较大的动物的骨骼比小动物的骨骼更为强健,较大的动物的直径和体长之比大于小动物.所以,当体格从小到大增加时,存在着以和动物尺寸的比例相应的规则的变形.11.一位物理学家想要研究光的性质.他想了解当光线从空气进入平滑的湖中,特别是在两种不同介质的交界处,光线的路径.12. 拥有一队卡车的一家公司面临着因卡车使用年限和油耗而增加的维修费用.13. 人们偏爱于计算机的速度.哪些计算机系统提供了最快的速度?14. 怎样提高我们的能力,使得每学期都能报名上最好的班级?15.怎样才能节约我们的一部分收入?16. 考虑在竞争市场情况下一家刚开始运转的生产单一产品的新公司.讨论该公司营业初期的短期和长期目标,这些目标会怎样影响到雇员工作的指派?该公司有必要决定短期运行的最大利润吗?17. 讨论利用模型来预测实际系统和利用模型来解释实际系统之间的差别.想象某些你要利用模型来解释实际系统的情景;类似地,想象你要利用模型来预测实际系统的其他情景.研究课题1.考虑冲泡咖啡的味道问题. 什么是影响味道的变量?哪些变量一开始可以忽略?假定除了水温外,已经固定了所有的变量,多数咖啡壶都用沸水以某种方式从底部的咖啡中蒸馏出滋味. 你认为用沸水是产生最佳滋味的最优方式吗?你将怎样检验这个子模型?你将收集什么样的数据以及怎样去收集这些数据?2.一家运输公司正在考虑用直升飞机在纽约市摩天楼之间运送人员,你被聘为顾问确定所需直升飞机的数量.精确地识别适当的问题,运用模型构建的过程来确定你所选定的变量之间的关系所需要的数据.当你着手进行时,可能需要重新定义你的问题.3.考虑酿酒问题. 提出若干商业制造商可能会有的目标.把考虑品位作为一个子模型,什么是影响品位的变量?哪些变量一开始就可以忽略?怎样把余下的变量关联起来?为确定这些关系,什么样的数据将是有用的?4.一对夫妇应该买房子还是租房子?因为抵押的费用上涨,直观上看,似乎存在一个抵押费用的价位,高于这个价格决不要去抵押贷款买房.什么变量决定了总的抵押费用?5.考虑一家诊所的运作问题.病人个人的病历档案必须保存,而会计程序是一项日常工作,该诊所应该购买或者租用一个小型的计算机系统吗?提出可能要考虑的目标.什么变量你会加以考虑?你怎样建立变量之间的关系?为决定你所选择的变量之间的关系,需要什么样的数据?为什么不同诊所对这个问题会有不同的解决办法?6.什么时候车主应该更新汽车?什么因素会影响到做出决定?哪些变量一开始可以忽略?识别你要的数据以决定所选择的变量之间的关系.7.一个人能跳多远?在1968年墨西哥城举行的奥运会上,美国的鲍勃·比蒙把世界纪录提高了10%,该记录一直保持到1996年的奥运会,列出影响跳远距离的变量.你认为墨西哥城的低空气密度可以解释这个10%的差别吗?8.上大学是一项可靠的金融投资吗?四年里没有收入,而且大学的费用极高,什么因素决定大学教育的总费用?怎么确定为使这项投资有利可图的必要条件?。
数学建模习题课1

习 题 课一、初等模型与常用的建模方法1. 奇偶校验法例11,9,8,6四个数字,问能否在余下的方格各填入一整数,使得方格区上的每一行每一列都构成等差数列?解这就产生了矛盾的结果,故所要求的填法不存在.例2 利用奇偶校验法证明,空间中不存在“有奇数个面,且每个面又都有奇数条边的多面体”.证 用反证法.必为奇数.另一方面,在多面体中,每两个相邻的面都有一条公共边,即多面体的棱,而且每一条棱又都为两个面所共有,每一条棱都被重复地计算了一次,于是产生了矛盾. 故由奇偶校验法知根本不存在具有奇数个面,且每个面又都有奇数条棱的多面体.例3. 证明这个多项式不能分解为两个整系数多项式的乘积.证用反证法.系数多项式的乘积,则必有(1)(2)2).另一方面,由(1)奇数立为奇数,因而(2)右端. 因此由奇偶校验法知满足条件积.2. 席位分配问题例4 比利时的d’Hondt曾提出过如下一种席位分配方案:将甲、乙、丙三个系的人数都用1,2,3…去除,然后将商从大到小排列,取前21个最大的商数考虑,规定在这21个商中,各系占有几个就分配给几个席位。
试通过数学建模探讨这种方法的合理性。
解以教材中甲、乙、丙三个系人数分别为103,63,34为例:系别人数 1 2 3 4 5 6 7 8 9 10 11 12甲 103 103 51.5 34.33 25.75 20.6 17.17 14.71 12.875 11.44 10.3 9.36 8.58乙 63 63 31.5 21 15.75 12.6 10.5 9 7.875 7 6.3丙 34 34 17 11.33 8.5 6.8从表中可以看出,按照比利时方法,在21个席位中甲占11席、乙占7席、丙占3席。
说明:(1)此席位分配方案明显不合理;(2)此方法与Q值方法比较有明显的缺陷,特别是当上述商值相等或十分接近时难以排序。
例5 某系共有1000名学生,其中235人住在A楼,333人住在B楼,432人住在C楼。
北师大版高中数学课件必修第1册第八章 数学建模活动(一)
§3 数学建模活动的主要过程
刷基础
解析
(1)由 17 时测得的平均行车速度为 3 km/h,n=-2(|17-12|-5)2+100=100,
600 ,n≤9,
n+10
代入 v= 3n32+ 00k0,n≥10n∈N+,可得1303002+ 00k=3,解得 k=1 000.
600n 600
600×9
§3 数学建模活动的主要过程
刷能力
解析
2.无标准答案,可以借助网络等资源查询相关资料,得到解决问题的思路.
≈0.69,无理数 e=2.718 28…)
§3 数学建模活动的主要过程
刷基础
解析
(1)∵ω=2,m0=160,mk=40,∴v=ωlnm0=2×ln160=2ln 4=4ln 2≈2.8,
mk
40
∴该单级火箭的最大理想速度为 2.8 千米/秒.
(2)∵m0≤10,ω=2,∴vmax=ωlnm0=2ln 10,
数学 BS 必修第一册
§1
§1 走近数学建模
§2
§2 数学建模的主要步骤
§3
§3 数学建模活动的主要过程
§3 数学建模活动的主要过程
刷基础
1.2021 年 12 月 9 日 15 时 40 分,神舟十三号“天宫课堂”第一课开讲!受“天宫课堂”的激励与鼓舞,
某同学对航天知识产生了浓厚的兴趣.通过查阅资料,他发现在不考虑气动阻力和地球引力等造成的影响时,
§3 数学建模活动的主要过程
刷能力
解析
k
k
(1)设 y1= (k≠0),y2=mx(m≠0),其中 x>0,当 x=9 时,y1= =2,y2=9m=7.2,解得 k=20,
数学建模作业1——火箭上升问题的模型建立教学提纲
数学建模作业1——火箭上升问题的模型建立题目:火箭上升问题的模型建立组员:摘要本文研究的是火箭上升问题,并针对有燃料和燃料已用尽两个问题分别建立了符合实际的数学模型。
在模型的求解过程中,通过运用MATLAB及微分方程,对建立的模型进行求解,得出了符合实际的结果。
关键字:火箭上升;数学模型;微分方程一、问题重述小型火箭初始质量为900千克,其中包括600千克燃料。
火箭竖直向上发射时燃料以15千克每秒的速度燃烧掉,由此产生30000牛顿的恒定推力。
当燃料用尽时引擎关闭。
设火箭上升的整个过程中,空气阻力与速度平方成正比,比例系数为0.4(千克/米),重力加速度取10米/秒 2(1)建立火箭升空过程的数学模型;(2)求引擎关闭瞬间火箭到达最高点的时间和高度。
二、基本假设1.火箭在喷气推动下作直线运动,火箭飞行时所受的地球自传与公转忽略不计。
2.火箭正常飞行,忽略其他因素对火箭飞行的影响。
3.假设产生影响的各个因素相互独立。
4.火箭上升初速度忽略不计,引擎足够强大。
5.火箭上升时所受到的重力加速度不变。
三、符号说明t :火箭上升过程的时间。
0t :第一个过程持续的时间。
M :第一阶段向上加速过程中火箭的质量。
m :第二阶段火箭剩余的质量。
f :火箭上升整个过程中空气阻力。
v :火箭的速度。
y :火箭上升的高度。
g :物体所受重力加速度。
F :火箭受到的恒定推力。
四、问题分析这是一个研究火箭竖直向上发射的问题。
火箭在竖直向上发射中,根据有燃料和燃料已用尽,可以分为两个阶段。
第一阶段是燃料产生推力的过程,第二阶段是燃料全部消耗之后的上升过程。
在第一阶段中,燃料燃烧产生的推力是恒定的,但随着燃料的不断消耗,火箭的质量是变化的,因此,火箭的速度以及加速度是变化的,由牛顿第二定律,根据速度与时间关系,建立微分方程组。
在第二阶段中,燃料已经完全消耗,因此,火箭的质量恒定。
引擎关闭即第一阶段终止第二阶段开始的时刻。
由于火箭运动受到阻力的作用,火箭先加速,后减速。
第1讲_什么是数学建模
合理化假设
显然该问题与瓶子和石子的形状及 其排列方式有关,为简单起见我们假设: • 瓶子是正方体的且不考虑瓶口的体积。 • 乌鸦投进的石子是大小相同的球体。 • 瓶子中摆放的方法如图1所示
图1
合理化假设
• 瓶子的边长是石子直径的整数倍,不妨 设为n倍(显然,如果不是整数倍的话, 那石子间的空隙会更大,不利于乌鸦喝 到水) • 石头内部渗进的水忽略不计。
与实际现象、数据比较, 检验模型的合理性、适用性
数学建模的全过程
现 实 世 界
现实对象的信息
验证 现实对象的解答
表述
(归纳)
解释
数学模型
求解 (演绎) 数学模型的解答
数 学 世 界
表述 求解 解释
根据建模目的和信息将实际问题“翻译”成数学问 题 选择适当的数学方法求得数学模型的解答 将数学语言表述的解答“翻译”回实际对象
原比(l/h)
0.6071 0.6071 0.6071 0.6071
身高(cm) 鞋跟高度(cm) 新比值
168 168 168 168 2.5 3.55 4.5 4.7748 0.6129 0.6151 0.6173 0.618
问题的检验
• 又如,按照上述模型,身高153CM,下肢 长为92CM的女士,应穿鞋跟高为6.6CM的 高跟鞋显得比较美。
明确建模目的 掌握对象特征
形成一个 比较清晰 的‘问题’
数学建模的一般步骤
模 型 假 设 建 立 模 型 针对问题特点和建模目的 作出合理的、简化的假设
在合理与简化之间作出折中
用数学的语言、符号描述问题
发挥想像力
尽量采用简单的数学工具
数学建模的一般步骤
模型 求解 模型 分析 模型 检验 模型应用 各种数学方法、软件和计算机技术 如结果的误差分析、统计分析、 模型对数据的稳定性分析
数学建模之概率统计-1
概率与统计
概率论中所研究的随机变量的分布都是 已知的。 统计学中所研究的随机变量的分布是未 知的或部分未知的,必须通过对所研究 的随机变量进行重复独立的观察和试验, 得到所需的观察值(数据),对这些数 据分析后才能对其分布做出种种判断, 即“从局部推断总体”。
统计学
给定一组数据,统计学可以摘要并且描述这
……
……
Matlab相关命令介绍
normfit 正态分布中的参数估计
[muhat,sigmahat,muci,sigmaci]=normfit(x,alpha) 对样本数据 x 进行参数估计,并计算置信度为 1-alpha 的置信区间 alpha 可以省略,缺省值为 0.05,即置信度为 95%
频率
随机试验进行次数
概率
基本知识
随机变量 数字特征(均值、方差、相关系数、特征函数…)
统计分析(假设检验、相关分析、回归分析…)
Matlab 中的随机函数
rand(m,n)
生成一个满足均匀分布的 m n 随机矩阵,矩阵的每
个元素都在 (0,1) 之间。
注:rand(n)=rand(n,n)
Matlab中的取整函数
fix(x) floor(x) ceil(x) round(x)
: 截尾取整,直接将小数部分舍去 : 不超过 x 的最大整数 : 不小于 x 的最小整数
: 四舍五入取整
取整函数举例
x1=fix(3.9);
x2=fix(-3.9); x3=floor(3.9); x4=floor(-3.2); x5=ceil(3.1); x6=ceil(-3.9); x7=round(3.9); x1=3 x2=-3 x3=3 x4=-4 x5=4 x6=-3 x7=4 x8=-3 x9=-4
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
承诺书我们仔细阅读了五一数学建模联赛的竞赛规则。
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与本队以外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其它公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们愿意承担由此引起的一切后果。
我们授权五一数学建模联赛赛组委会,可将我们的论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等)。
我们参赛选择的题号为(从A/B/C中选择一项填写):A我们的参赛报名号为:参赛组别(研究生或本科或专科):本科所属学校(请填写完整的全名)参赛队员(打印并签名) :1.2.3.日期:2015 年 5 月3日获奖证书邮寄地址:邮政编码:收件人姓名:联系电话:编号专用页竞赛评阅编号(由竞赛评委会评阅前进行编号):评阅记录评阅人评分备注裁剪线裁剪线裁剪线竞赛评阅编号(由竞赛评委会评阅前进行编号):参赛队伍的参赛号码:(请各参赛队提前填写好):2152#题目不确定条件下的最优路径的搜索模型摘 要本文以一般交通网络图为研究对象,在每条路径中每个路段的行驶时间的均值和方差已知的基础上,构造最优路径可靠性模型,保证相同高的概率到达时行驶总时间最短,然后把时间相关性和空间相关性融入建立的模型中,层层求解。
最后根据所建立的模型分析算法,得出最优路径。
针对问题一,考虑到在现实中从起点到中间的路径错综复杂,且每条路径的路段数量不确定,我们绘制了一般交通网络图。
在已知各路段行驶时间的均值与标准差的前提下,给出最优路径的定义。
基于各路段的行驶时间服从正态分布的假设下利用卷积公式,针对每条路径给出以行驶总时间为变量的正态密度函数。
最后对所得函数变形处理,得到判断最优路径的数学模型,并根据模型求出示例一的最优路径。
针对问题二,首先利用卷积公式和10-规划,给出最优路径定义的条件下和给出最优路径的搜索算法,运用Matlab 和excel 从所有路径中,搜索最优路径;并将其最优路径搜索算法运用到实际问题中,简化图如图2,搜索最优路径为Q v v v p →→→→412,该路径的均值为17.77,标准差为3.66。
针对问题三,在时间相关性方面构造Copula 函数和最大似然估计函数,利用excel 软件,得到从51v v →,在00:8~00:7堵塞时00:9~00:8对其他路段的均值和标准差的影响。
在空间相关性方面,构造Moren 模型和空间自相关分析,得到在一定时间下,某个路段出现交通事故,对其他相关路段的影响。
结论为第8路段的均值为12,标准差为3.6。
针对问题四,赋予路段均值和方差的等权重,作为该路段的行驶参考时间,最终转化为Dijkstra 算法,求得最优路径为Q v v P →→→41,最优行驶时间为13.65。
本文逻辑严谨,切入点独到,综合运用多种模型及软件,结果可靠且多样化。
关键词:最优路径;正态分布;Moren 模型;灵敏度分析;Matlab§1问题的重述一、背景知识1.中国城市交通的基本情况近年来,经济迅速的增长和城市化进程的加快使城市人口不断增加,城市人口的大量聚集和私家车人均拥有量的激增使城市交通面临着越来越艰巨的任务,交通拥堵、交通安全、高效交通等问题一涌而出,交通问题日益严峻,从一年一度的“春运”现象就可见一斑,不仅给居民的日常生活造成了严重的影响,还制约了城市及周边城市的发展。
在复杂的交通环境下,如何综合考虑多方面因素,寻找一条可靠、快速、安全的最优路径,已成为每个人所迫切解决的问题。
2.交通网(1)交通点:交通点即通常的汽车站、火车站、道路的交叉点、机场、港口等交通结点。
(2)交通线:交通线即连接点与点的铁路、公路、水路以及航线的交通路线。
在特定的地域范围内,根据地区经济的发展和人们活动的需求,各种现代交通运输方式联合,各种交通线和交通点交织,形成了不同形式和层次的交通运输网,简称交通网。
其布局受到经济、社会、技术和自然的影响和制约。
按交通运输方式分类,形成了铁路运输网、公路运输网、水路运输网、航空运输网和管道运输网。
不同运输方式结合形成综合交通运输网。
3.传统的最优路径传统的最优路径是基于理想交通状况下分析得到的平均总行驶时间最短的路径。
在这种情况下每条路段的行驶时间是确定的,可以用经典的最短路算法(Dijkstra算法)来寻找最短路径。
这也是大多数车辆路径导航系统寻找最优路径所用的算法。
但传统路径存在相当大的弊端,在现实生活中,交通事故,天气情况,流量等一系列不确定因素都会让最优路径失去最优性。
4.最优路径的改良不同于传统路径,改良后的最优路径把在现实生活中可能遇到的一系列不确定因素如天气、车流量等纳入考虑范围内,将每条路段的行驶时间不确定化,使其近似服从一个随机分布,建立模型,从而给出一个相对的最优路径。
二、相关数据1.原问题第一问的示例交通网络图及其数据;三、要解决的问题以中国矿业大学到徐州火车站为例,假设行驶时间是随机变量,假设已知每条路段行驶时间的均值和标准差。
若走绕城快速路平均33分钟到达,标准差只有1分钟,若走市区道路,平均30分钟到达,标准差15分钟。
若采用传统的最优路径算法,应选择平均时间较短的市区道路,但现实往往相反。
最优路径不仅要考虑平均时间最短,还要考虑不确定性条件下车辆准时到达的可靠性,在此要求下,解决下列问题:1.问题一:建立适用于一般交通网络的数学模型,定量分析车辆行驶时间的不确定性,给出在不确定条件下最优路径才定义和数学表达式,并应用于示例一。
2.问题二:根据第一问建立的模型设计算法搜索最优路径,并应用于具体的交通网络中,尽可能地从理论上分析算法的收敛性和复杂性。
3.问题三:建立数学模型描述交通行驶时间之间的相关性,将其应于前二问的最优路径搜索中,并设计算法解决这一问题,给出算例证明其有效性,尽可能的分析算法的收敛性、复杂性等性质。
4.问题四:从不确定性条件下交通网络的实际情况出发,在合理假设下,进一步完善前三问的数学模型和相关算法。
或提出一种或多种与前三问不同的最优路径的定义方法,建立相关的数学模型并设计算法,应用数值算例验证算法的有效性并尽可能从理论上分析算法的收敛性、复杂性等性质。
§2模型的假设1.假设每条路段行驶时间服从正态分布; 2.假设每条路段行驶时间的均值和标准差;3.已知假设所有不确定因素对行驶时间的影响都可以用均值和标准差体现;§3名词解释与符号说明一、名词解释1.均值:表示一组数据集中趋势的量数,反映数据集中趋势的一项指标。
2.标准差:标准差是方差的算术平方根,能反映一个数据集的离散程度。
3.相关性:数据间的关联程度,一个数据的取值会受到其他数据的影响。
4.算法的复杂性:算法复杂性是输入规模的函数,一般为避免不同输入对算法行为造成的巨大差别,考察所有输入规模为n 的复杂性。
二、符号说明序号 符号 符号说明 1 P 路径起点 2 Q 路径终点 3 i L 第i 条路径4 i C 第i 条路径所用的总时间5 i T 每条路径的第i 个路段所用时间6 i μ 每条路径的第i 个路段所用时间的均值7 i σ每条路径的第i 个路段所用时间的方差 8 )(t P T T 的正态分布密度函数 9 i v 第i 交通点 10 i y 10-变量 11 F 多维分布函数 12 C Copula 函数13∧i θ估计边缘分布参数§4模型的建立与求解一、问题一的分析与求解 1.对问题的分析相比于传统最优路径,改良后的最优路径需要把现实中的不确定因素纳入考虑范围,对此我们将各路段行驶时间变量化,引入正态分布,首先将现实中的交通网络简化为一般网络,基于各路段的行驶时间服从正态分布的假设,利用正态分布可加性把每条路径上的多个路段整合为一个路段。
然后利用卷积公式,给出每条路径的的正态密度函数,变形得出可靠性模型,并将相同概率下及时到达的路径中用时最短的路径定义为最优路径。
定义1 最优路径定义在相同概率到达终点的条件下,平均行驶总时间最短的路径即为最优路径。
2.对问题的求解⑴ 模型的准备如下图所示,绘制如下图。
不妨假设以P 为起点,Q 为终点的道路有若干条途径,中每条路径的均值与方差已知。
…………1t 2t nt 1-n t 1t 2t 1-n t nt PQ图1 一般路径图设随机变量),(~2111σμN T ,),(~2222σμN T 且1T 与2T 独立,则有),(~22212121σσμμ+++=N T T T 。
证明:21T T T +=仍在),(+∞-∞上取值,利用卷积公式222222121121]})()([21exp{21)(dt t t t t P T ⎰∞+∞--+---=σμσμσπσ对上式被积公式的指数部分按2t 的幂次展开,合并同类项得:222221211)()(σμσμ-+--t t t =222122122)()(σσμμ+--+-t A B t A 其中222111σσ+=A ,222211σμσμ+-=t B 。
代回原式可得:22222122121])(2exp[})(21exp{21)(dt A B t A t t P T ⎰∞+∞---⋅+---=σσμμσπσ利用正态密度函数的正则性,上式的积分应为Aπ2,于是})(21exp{)(21)(22212212221σσμμσσπ+---+=t t P T这正是均值为21μμ+,方差为2221σσ+的正态密度函数,有下式成立:),(~22212121σσμμ+++=N T T T 。
下用数学归纳法证明: 对于任意n ,有下式成立:∑∑∑===ni n i ni i i iN T1112),(~σμ。
当2,1=k 时,显然成立。
假设当k n =时成立,即为∑∑∑===ki k i ki i i iN T1112),(~σμ,则对于当1+=k n 时,∑+==11k i i T T1212111212122221])()([21exp 21)(1+∞∞-+++==++⎰∑∑⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-+---++=k k k k k i i ki i k k k T dt t t t t P σμσμσσσπσ 同样对于上式函数中指数部分按照1+k t 展开可得:121112112121})(2exp{)(21exp 21)(+∞∞-++=+=+⎰∑∑--⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=k k k i i k i i k T dt A B t A t t P σμσσπσ其中∑+==1121k i i A σ211121++==+-=∑∑k k ki iki it B σμσμ。