经典的数学建模例子

经典的数学建模例子
经典的数学建模例子

一、摘要

SARS

SARS就是传染性非典型肺炎,全称严重急性呼吸综合症(Severe Acute Respiratory Syndromes),简称SARS,是一种因感染SARS相关冠状病毒而导致的以发热、干咳、胸闷为主要症状,严重者出现快速进展的呼吸系统衰竭,是一种新的呼吸道传染病,传染性极强、病情进展快速。

当一种传染病流行的时候,会给人们的工作学习带来很大的不变,能有效地进行隔离、预防,会大大减少人员的得病率,当一种传染病开始流行时,在一定的条件下其趋势就像真菌的繁殖曲线,如果能通过计算预测但大概推算出其发病率高峰时期,及时的隔离预防。那会给社会人力带来很大的方便,当年SARS的爆发给我们带来和大的不便和损失,因此本论文就以SARS为例,来研究传染病的传播规律、为预测和控制传染病蔓延创造条件和帮助。

1

二、正文

1、模型的背景问题描述

SARS(Severe Acute Respiratory Syndrome,严重急性呼吸道综合症, 俗称:非典型肺炎)是21世纪第一个在世界范围内传播的传染病。SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响,我们从中得到了许多重要的经验和教训,认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条件的重要性。

要求:(1)建立传染病传播的指数模型,评价其合理性和实用性。

(2)建立一个适合的模型,说明为什么优于问题1中的模型;特别要说明怎样才能

3

建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型,这样做的困难在哪里?对于卫生部门所采取的措施做出评论,如:提前或延后5天采取严格的隔离措施,对疫情传播所造成的影响做出估计。表中提供的数据供参考。

(3)说明建立传染病数学模型的重要性。

2、模型假设

(一)答;

从上列图表可知道在4月20到5月7日期已确诊的发病人总数呈指数增长趋势5月20到6月1日增长缓慢,6月1日到6月12日总数几乎不变。其形式与生物学中真菌繁殖总数相似。

从表格和准备中,作如下假设。

1、不考虑SARS在人体中的潜伏期,也就是说当人一旦传染就表现出来立即就具有传染

性。

2、当健康者满足一地条件时,健康者才被传染。

3、整个发病期间为自然状态也就是无人为外界干扰,政府等其它形式进行隔离预防。

4、忽略特殊情况,如个别人体质弱或强的。

假定初始时刻得病例数为M0。平均每位病人每天可传染N个人,可传染他人的时间为T 天。则在T天内,病例数目的增长随着时间t(单位天)的关系是;

M(t)=M0(1+N)t

如果不考虑对传染期的限制则病例数将按照指数规律增长考虑,当传染期T的作用后,变化将显著偏离指数规律,增长速度会放慢。把达到T天的病例从可以引发直接传染的基数中去掉,为了方便,从开始到高峰期间,均采用同样的N值,(从拟合这一阶段的数据库定出),到达高峰之后在10天的范围内逐步调整N值,到比较小,然后保持不变,拟合后在控制阶段的全部数据。

评价及其合理性和实用性;

本模型主要有三个参数M0、N、T,且都具有实际意义。T可理解为平均每个病人在被发现前后可以造成直接传染的期限,在此期限后失去传染能力,可能原因是被隔离、病愈或死去等等。N表示某种社会条件下平均每位病人每天传播的人数(但并非文中所述的一个病人的感染他人的平均概率)。整个模型抓住了SARS传播过程中两个主要特征:传染期T和传染率N,反映了SARS的传播过程。使人很容易理解该模型。

模型灵活

通过调整M0、N、T值,就可以描述不同地区,不同环境下SARS的初期传播规律预测准确

通过模型对表格的调查结果进行了分析,得到的预测值与实际统计数据较接近。可大致预测出疫情的爆发点和发展趋势。

预期模型的缺点:

1、对于如何确定对于三个参数M0、N、T,未给出一般的原则或算法,只能通过对

于已发病地区的数据进行拟合得出。按照作者的表述,N值是以病发高峰为界取各段的平均值作为传染概率,虽然简化了运算,但是在现实情况下,不同地区的N值是不同的。在实际应用中,如果没有一定量的数据,是无法得出N值的。在我们对该模型

5

进行拟合事发现,对于M0 、N 、 T 作者未给出调整的标准和相关理论,所以我们很难重复该求解过程。

2、当需要对某一地区进行疫情分析时,还需考虑到该地区相对于表格所给的人群这类人口密集,人员流动性大的城市之间的差异。地域因素会造成不同地区的N 值不同(如人口密度和人口流动大的城市若爆发传染病,初期的N 值会比人口密度和人口流动小的城市大,等等),而很难找到地域因素几乎相同的两城市。所以此作法可能导致预测结果相差较大。

综上所述,该模型能较好的反映SARS 传染的特征性,具有一定的实际意义。但是,参数的取值包含有一定的主观因素,且需要大量的数据进行拟合,且未给出调整的标准和相关理论,在实际应用中实用价值不大。

(二)答:

模型假设

1、在疾病传播期内所考察地区的总人数N 不变,既不考虑生死,,也不考虑迁移。人群分为易感染者和已感染者两类,时刻t 这两类人在总人数中所占比例分别记为s (t )和i (t )。

2、每个病人每天有效接触的平均人数是常数k ,k 称为日接触率。当病人与健康者有效接触时,使健康者受感染为病人。

问题分析

根据假设,可知,人群分为两类,一是健康者,二是病人,只要一类人群随时间的变化规律知道,这另一类人群也可马上求解。由于传染病过程中通常取病人为研究对象,所以决定求解病人随时间的变化规律。

3、模型分析、建立

对于t 时刻,病人的增加率为kNsi ,即

di

N kNsi dt

= (1) 又因为

S (t )+i (t )=1 (2) 再令初始时刻的病人比例为i0,这 (1),(0)0di

ki i i i dt

=-= (3) 显然此为logistic 模型,它的解为 1

()1(1/01)kt

i t i e -=

+- (4)

参数的确定

通过对图表的累计病例数用spss 进行曲线拟合,结果如下

可得拟合的函数关系式为

1

1/25250.001(0.865)t

y =

+,y=N*i

通过取一系列t 来估计出相应的k 值,结果如下 时间 20 30 40 50 60

k 值大小 0.1919 0.1763 0.1685 0.1638

0.1607

7

由图像可知,当t 较大时,曲线拟合的数据与实际测量值越接近,所以就取t=60时所对应的k 值,即0.1607。此值可以近似看做当政府没有采取措施,即传染病的自然传染能力大小。但同时根据附件1的求解方法,我们计算了4月20日到4月29日期间每日的k 值大小,再求平均,得k =0.169346。对于k 和k 之间的差异,这是由于模型1并未考虑到政府控制前和控制后k 值将改变,且k1>k2。所以由于k 只考虑控制前,所以比k 要略大,我们考虑传染病的每天平均自然传染人数时,取值为k =0.169346。但由于此模型未考虑到病人会被治愈而成为健康者,所以在模型1的基础上进行改进,建立了模型2。 模型2

在模型1的假设条件下增加的条件为,

1,每天被治愈的病人数或死于该传染病人数占病人总数的比例为常数p 。病人治愈后由于获得了免疫能力,同时也由于心理作用,更加保护自己,所以可以假设治愈后再次感染的几率为0,且该种人群在总人群中所占有的比例为u (t )。 不难看出,考虑到假设3,模型1中的(1)式应修改为

,(0)0di

N kNsi pNi i i dt

=-= (5)

而且对于健康者,其增加率为

,(0)0ds

ksi s s dt

=-= (6) 对于移出者而言,其增加率为

du

N pNi dt

= (7) 由于人群只由健康者,病人和移出者组成,所以

S(t)+i(t)+u(t)=1 (8)

4、模型求解

查资料,得到2003年北京市市区总人口数目为698.8万人

从而可以得到初始条件i0= 339/(698.8*10^(-4))=4.851*10^(-5 ) ,s0= 0.99995149(取4月20号为初始条件)

同时根据附件2中的死亡累计和治愈累计,求得每日的移出率p ,在求平均值得到p =0.05121。 在模型一中求得k =0.169346;将上述参数代入(5)式和(6))式,求得数值解和绘制的图像

00.10.20.30.40.50.60.70.80.91

0.050.10.150.20.250.30.35

9

050

100150200

0.5

1

1.5

2

2.5

6

时间

能感染的病人数

由图像可得i (t )随时间的推移先逐渐变大,之后变小,趋向于0,s(t)随时间的推移而逐渐减小,根据常识,一种传染病中的病人比例最终是为0,由此模型2还是比较符合客观事实的,但从图像中大致可以判断i (t )=0时大约要经过225多天。这与实际过程中大约经过100多天北京的SARS 就平息存在较大误差,仔细分析,我们发现该模型忽略了SARS 的潜伏期,实际上健康人与SARS 患者接触后虽然被感染了,但还处于潜伏期,没有传染能力。所以将模型2进行改进,得到模型3。 模型3 模型假设

1,将人群分为四类,分别为健康人群,能感染的病人SARS ,SARS 潜伏者和移出者(包括SARS 的死亡者和治愈者),他们在人群中的比重分别为 s(t),i(t),w(t),u(t);其中已确诊病人和SARS 潜伏者统称为SARS 病毒携带者,记 为x1(t ),表示其t 时刻的人数,人口总人数为N 。

2 ,每个病人每天有效接触的平均人数是常数k ,k 称为日接触率。当病人与

健康者有效接触时,使健康者受感染为病人。

3, SARS 潜伏者无传染能力,但最终会成为病人,具有传染能力。

问题分析

该模型比起模型2更为复杂,在该模型中还必须将SARS 病毒携带者分为两类,显然增加了计算难度,在此中还应该考虑潜伏周期T 。 模型求解

在t 时刻SARS 病毒携带者x1(t)=Ni(t)+Nw(t) (9) SARS 病毒携带者的增长率为

1()

()()()dx t Ni t ks t Ni t p dt

=- (10) 健康人的增长率为

()()ds

i t ks t dt

=- (11) 移出者的增长率为

du

ip dt

= (12) 潜伏者的增长率为

()

()()()()dw t i t ks t i t T ks t T dt

=--- (13) 确诊病人的增长率为

()()

()di t dw t T i t p dt dt

-=- (14) 除此之外,还有一条公式,为

S(t)+i(t)+w(t)+u(t)=1 (15) 由(9)到(16)式联立,可得到

()

()()()di t i t T ks t T i t p dt

=--- (16) 将(14)和(17)式联立,可得

()

()()dw t T i t T ks t T dt

-=-- (17) 将t-T 用t 代替,可得

()

()()dw t i t ks t dt

= (18) 在对(16)式两边对t 进行求导,可得

()()()()0d s t d i t d w t d u t

d t d t d t d t

+++= (19) 结合(11),(12),(18)可求得

()

()di t i t p dt

=- (20) 最终对(11),(18),(20)联立的方程组进行数值求解,可得图像如下

11

-0.2

00.20.40.60.81

1.2

00.050.10.150.20.250.30.350.40.450.5

0.511.522.5

33.56

时间

能感染的病人数

从图像中我们可以观察到在300多天时i (t )会接近于0,这比模型2还要久,,因此,我们还把政府的干预考虑进来,也就得到了模型4。

由于时间限制,模型4中考虑的因素更多,所以一时没能解决,也就导致了第二问实际上还不能完全解决,但是我们已经有了思路,即再引入一类人群,就是隔离人群,通过引入该人群,实际上是改变了病人的有效接触人数k ,我们根据4月29日之后的实际数据,求得每日的k 值,再求平均,得2k =0.019413;我们想采用分段函数,即确定一个时刻t ,k 值改变的时刻,在这个时刻前与后都可以适用模型3。只是在考虑t 时刻后,它的初始条件为4月29日的数据。通过t 的改变,可以解决第二问中政府早五天调控和晚五天调控的差别。

5、模型的应用与推广

模型是为了一定目的,对客观事物的一部分进行简缩、抽象、提炼出来的原型的替代物,集中反映了原型中人们需要的那一部分特征。

数学建模就是指对于一个现实对象,为了一个特定目的,根据其内在规律,作出必要的简化假设,运用适当的数学工具,得到的一个数学结构。其意义在于用数学方法解决实际问题。

能够合理有效地建立数学模型不仅在工业,商业,医疗卫生都有很大帮助。就拿SARS 而言,它可以使人们有效地采取方案措施进行解决。可以高效,高质的完成事务。它是从一个定量的角度分析解决问题,去解决人类的实际问题。对人类的生活带来很大方便,能够很好的利用数学建模去解决生活生产,在科技方面也有很大的帮助。用它去解决人们的生活生产会给人力带来很大的帮助,有效地解决实际问题。感谢老师给了我们这次机会,通过这次学习对建模也有很大的乐趣 ,我会继续学习去解决更多实际问题为学习和生活带来方便。

13

数学建模经典案例:最优截断切割问题复习进程

数学建模经典案例:最优截断切割问题

建模案例:最优截断切割问题 一、 问 题 从一个长方体中加工出一个已知尺寸、位置预定的长方体(这两个长方体的对应表面是平行的),通常要经过6 次截断切割.设水平切割单位面积的费用是垂直切割单位面积费用的r 倍.且当先后两次垂直切割的平面(不管它们之间是否穿插水平切割)不平行时,因调整刀具需额外费用e.试设计一种安排各面加工次序(称“切割方式”)的方法,使加工费用最少. 二、 假 设 1、假设水平切割单位面积的费用为r ,垂直切割单位面积费用为1; 2、当先后两次垂直切割的平面(不管它们之间是否穿插水平切割)不平行时,调整刀具需额外费用e ; 3、第一次切割前,刀具已经调整完毕,即第一次垂直切割不加入刀具调整费用; 4 、每个待加工长方体都必须经过6次截断切割. 三、 模型的建立与求解 设待加工长方体的左右面、前后面、上下面间的距离分别为 a0、b0 、c0 ,六个切割面分别位于左、右、前、后、上、下,将它们相应编号为M1、M2、M3、M4、M5、M6,这六个面与待加工长方体相应外侧面的边距分别为 u1、u2、u3、u4、u5、u6.这样,一种切割方式就是六个切割面的一个排列,共有P 66720= 种切割方式.当考虑到切割费用时,显然有局部优化准则:两个平行待切割面中,边距较大的待切割面总是先加工. 由此准则,只需考虑 P 6622290!!! ??=种切割方式.即在求最少加工费用时,只 需在90个满足准则的切割序列中考虑.不失一般性,设u1≥u2,u3≥u4,u5≥u6,故只考虑M1在M2前、M3在M4前、M5在M6前的切割方式.

数学建模典型例题

一、人体重变化 某人的食量是10467焦/天,最基本新陈代谢要自动消耗其中的5038焦/天。每天的体育运动消耗热量大约是69焦/(千克?天)乘以他的体重(千克)。假设以脂肪形式贮存的热量100% 地有效,而1千克脂肪含热量41868焦。试研究此人体重随时间变化的规律。 一、问题分析 人体重W(t)随时间t变化是由于消耗量和吸收量的差值所引起的,假设人体重随时间的变化是连续变化过程,因此可以通过研究在△t时间内体重W的变化值列出微分方程。 二、模型假设 1、以脂肪形式贮存的热量100%有效 2、当补充能量多于消耗能量时,多余能量以脂肪形式贮存 3、假设体重的变化是一个连续函数 4、初始体重为W0 三、模型建立 假设在△t时间内: 体重的变化量为W(t+△t)-W(t); 身体一天内的热量的剩余为(10467-5038-69*W(t)) 将其乘以△t即为一小段时间内剩下的热量; 转换成微分方程为:d[W(t+△t)-W(t)]=(10467-5038-69*W(t))dt; 四、模型求解 d(5429-69W)/(5429-69W)=-69dt/41686 W(0)=W0 解得: 5429-69W=(5429-69W0)e(-69t/41686) 即: W(t)=5429/69-(5429-69W0)/5429e(-69t/41686) 当t趋于无穷时,w=81; 二、投资策略模型 一、问题重述 一家公司要投资一个车队并尝试着决定保留汽车时间的最佳方案。5年后,它将卖出所有剩余汽车并让一家外围公司提供运输。在策划下一个5年计划时,这家公司评估在年i 的开始买进汽车并在年j的开始卖出汽车,将有净成本a ij(购入价减去折旧加上运营和维修成本)ij

数学建模典型例题(二)

6 小行星的轨道模型 问题 一天文学家要确定一颗小行星绕太阳运行的轨道,他在轨道平面内建立以太阳为原点的直角坐标系,在两坐标轴上取天文测量单位(一天文单位为地球到太阳的平均距离:1.4959787×1011m ).在5个不同的时间对小行星作了5次观察,测得轨道上5个点的坐标数据如表6.1. 表6.1 坐标数据 由Kepler (开普勒)第一定律知,小行星轨道为一椭圆.现需要建立椭圆的方程以供研究(注:椭圆的一般方程可表示为 012225423221=+++++y a x a y a xy a x a . 问题分析与建立模型 天文学家确定小行星运动的轨道时,他的依据是轨道上五个点的坐标数据: (x 1, y 1), (x 2, y 2), (x 3, y 3), (x 4, y 4), (x 5, y 5). 由Kepler 第一定律知,小行星轨道为一椭圆.而椭圆属于二次曲线,二次曲线的一般方程为012225423221=+++++y a x a y a xy a x a .为了确定方程中的五个待定 系数,将五个点的坐标分别代入上面的方程,得 ???? ?????-=++++-=++++-=++++-=++++-=++++.122212221222122212225554253552251454424344224 135342 3333223125242 232222211514213112211y a x a y a y x a x a , y a x a y a y x a x a ,y a x a y a y x a x a ,y a x a y a y x a x a ,y a x a y a y x a x a 这是一个包含五个未知数的线性方程组,写成矩阵

数学建模常用的十种解题方法

数学建模常用的十种解题方法 摘要 当需要从定量的角度分析和研究一个实际问题时,人们就要在深入调查研究、了解对象信息、作出简化假设、分析内在规律等工作的基础上,用数学的符号和语言,把它表述为数学式子,也就是数学模型,然后用通过计算得到的模型结果来解释实际问题,并接受实际的检验。这个建立数学模型的全过程就称为数学建模。数学建模的十种常用方法有蒙特卡罗算法;数据拟合、参数估计、插值等数据处理算法;解决线性规划、整数规划、多元规划、二次规划等规划类问题的数学规划算法;图论算法;动态规划、回溯搜索、分治算法、分支定界等计算机算法;最优化理论的三大非经典算法:模拟退火法、神经网络、遗传算法;网格算法和穷举法;一些连续离散化方法;数值分析算法;图象处理算法。 关键词:数学建模;蒙特卡罗算法;数据处理算法;数学规划算法;图论算法 一、蒙特卡罗算法 蒙特卡罗算法又称随机性模拟算法,是通过计算机仿真来解决问题的算法,同时可以通过模拟可以来检验自己模型的正确性,是比赛时必用的方法。在工程、通讯、金融等技术问题中, 实验数据很难获取, 或实验数据的获取需耗费很多的人力、物力, 对此, 用计算机随机模拟就是最简单、经济、实用的方法; 此外, 对一些复杂的计算问题, 如非线性议程组求解、最优化、积分微分方程及一些偏微分方程的解⑿, 蒙特卡罗方法也是非常有效的。 一般情况下, 蒙特卜罗算法在二重积分中用均匀随机数计算积分比较简单, 但精度不太理想。通过方差分析, 论证了利用有利随机数, 可以使积分计算的精度达到最优。本文给出算例, 并用MA TA LA B 实现。 1蒙特卡罗计算重积分的最简算法-------均匀随机数法 二重积分的蒙特卡罗方法(均匀随机数) 实际计算中常常要遇到如()dxdy y x f D ??,的二重积分, 也常常发现许多时候被积函数的原函数很难求出, 或者原函数根本就不是初等函数, 对于这样的重积分, 可以设计一种蒙特卡罗的方法计算。 定理 1 )1( 设式()y x f ,区域 D 上的有界函数, 用均匀随机数计算()??D dxdy y x f ,的方法: (l) 取一个包含D 的矩形区域Ω,a ≦x ≦b, c ≦y ≦d , 其面积A =(b 一a) (d 一c) ; ()j i y x ,,i=1,…,n 在Ω上的均匀分布随机数列,不妨设()j i y x ,, j=1,…k 为落在D 中的k 个随机数, 则n 充分大时, 有

数学建模经典案例:最优截断切割问题

建模案例:最优截断切割问题 一、 问 题 从一个长方体中加工出一个已知尺寸、位置预定的长方体(这两个长方体的对应表面是平行的),通常要经过 6 次截断切割.设水平切割单位面积的费用是垂直切割单位面积费用的r 倍.且当先后两次垂直切割的平面(不管它们之间是否穿插水平切割)不平行时,因调整刀具需额外费用 e.试设计一种安排各面加工次序(称“切割方式”)的方法,使加工费用最少. 二、 假 设 1、假设水平切割单位面积的费用为r ,垂直切割单位面积费用为1; 2、当先后两次垂直切割的平面(不管它们之间是否穿插水平切割)不平行时,调整刀具需额外费用e ; 3、第一次切割前,刀具已经调整完毕,即第一次垂直切割不加入刀具调整费用; 4 、每个待加工长方体都必须经过6次截断切割. 三、 模型的建立与求解 设待加工长方体的左右面、前后面、上下面间的距离分别为 a0、b0 、c0 ,六个切割面分别位于左、右、前、后、上、下,将它们相应编号为M1、M2、M3、M4、M5、M6,这六个面与待加工长方体相应外侧面的边距分别为 u1、u2、u3、u4、u5、u6.这样,一种切割方式就是六个切割面的一个排列,共有P 66720= 种切割方式.当考虑到切割费用时,显然有局部优化准则:两个平行待切割面中,边距较大的待切割面总是先加工. 由此准则,只需考虑 P 6622290!!! ??=种切割方式.即在求最少加工费用时, 只需在90个满足准则的切割序列中考虑.不失一般性,设u1≥u2,u3≥u4,u5≥u6,故只考虑M1在M2前、M3在M4前、M5在M6前的切割方式. 1、 e=0 的情况

数学建模典型例题

一、人体重变化 某人得食量就是10467焦/天,最基本新陈代谢要自动消耗其中得5038焦/天。每天得体育运动消耗热量大约就是69焦/(千克?天)乘以她得体重(千克)。假设以脂肪形式贮存得热量100% 地有效,而1千克脂肪含热量41868焦。试研究此人体重随时间变化得规律. 一、问题分析 人体重W(t)随时间t变化就是由于消耗量与吸收量得差值所引起得,假设人体重随时间得变化就是连续变化过程,因此可以通过研究在△t时间内体重W得变化值列出微分方程。 二、模型假设 1、以脂肪形式贮存得热量100%有效 2、当补充能量多于消耗能量时,多余能量以脂肪形式贮存 3、假设体重得变化就是一个连续函数 4、初始体重为W0 三、模型建立 假设在△t时间内: 体重得变化量为W(t+△t)—W(t); 身体一天内得热量得剩余为(10467—5038-69*W(t)) 将其乘以△t即为一小段时间内剩下得热量; 转换成微分方程为:d[W(t+△t)-W(t)]=(10467—5038-69*W(t))dt; 四、模型求解 d(5429—69W)/(5429-69W)=-69dt/41686 W(0)=W0 解得: 5429-69W=(5429-69W0)e(-69t/41686) 即:

W(t)=5429/69—(5429-69W0)/5429e(-69t/41686) 当t趋于无穷时,w=81; 二、投资策略模型 一、问题重述 一家公司要投资一个车队并尝试着决定保留汽车时间得最佳方案。5年后,它将卖出所有剩余汽车并让一家外围公司提供运输。在策划下一个5年计划时,这家公司评估在年i得开始买进汽车并在年j得开始卖出汽车,将有净成本aij(购入价减去折旧加上运营与维修成本).以千元计数aij得由下面得表给出: 请寻找什么时间买进与卖出汽车得最便宜得策略。 二、问题分析 本问题就是寻找成本最低得投资策略,可视为寻找最短路径问题.因此可利用图论法分析,用Dijkstra算法找出最短路径,即为最低成本得投资策略。 三、条件假设 除购入价折旧以及运营与维护成本外无其她费用; 四、模型建立 二 5 11 7 三6 4

2003全国大学生数学建模竞赛B题优秀论文(出题人亲作)

2003高教社杯全国大学生数学建模竞赛 B 题参考答案 注意:以下答案是命题人给出的,仅供参考。各评阅组应根据对题目的理解及学生的解答,自主地进行评阅。 问题分析: 本题目与典型的运输问题明显有以下不同: 1. 运输矿石与岩石两种物资; 2. 产量大于销量的不平衡运输; 3. 在品位约束下矿石要搭配运输; 4. 产地、销地均有单位时间的流量限制; 5. 运输车辆每次都是满载,154吨/车次; 6. 铲位数多于铲车数意味着最优的选择不多于7个产地; 7. 最后求出各条路线上的派出车辆数及安排。 运输问题对应着线性规划,以上第1、2、3、4条可通过变量设计、调整约束条件实现; 第5条使其变为整数线性规划;第6条用线性模型实现的一种办法,是从1207 10 C 个整数规划中取最优的即得到最佳物流;对第7条由最佳物流算出各条路线上的最少派出车辆数(整数),再给出具体安排即完成全部计算。 对于这个实际问题,要求快速算法,计算含50个变量的整数规划比较困难。另外,这是一个二层规划,第二层是组合优化,如果求最优解计算量较大,现成的各种算法都无能为力。于是问题变为找一个寻求近优解的近似解法,例如可用启发式方法求解。 调用120次整数规划可用三种方法避免:(1)先不考虑电铲数量约束运行整数线性规划,再对解中运量最少的几个铲位进行筛选;(2)在整数线性规划的铲车约束中调用sign 函数来实现;(3)增加10个0-1变量来标志各个铲位是否有产量。 这是一个多目标规划,第一问的目标有两层:第一层是总运量(吨公里)最小,第二层是出动卡车数最少,从而实现运输成本最小。第二问的目标有:岩石产量最大;矿石产量最大;运量最小,三者的重要性应按此序。 合理的假设主要有: 1. 卡车在一个班次中不应发生等待或熄火后再启动的情况; 2. 在铲位或卸点处因两条路线(及以上)造成的冲突时,只要平均时间能完成任务即 可,不进行排时讨论; 3. 空载与重载的速度都是28km/h ,耗油相差却很大,因此总运量只考虑重载运量; 4. 卡车可提前退出系统。 符号:x ij ~ 从i 号铲位到j 号卸点的石料运量 单位 吨; c ij ~ 从i 号铲位到j 号卸点的距离 公里; T ij ~ 从i 号铲位到j 号卸点路线上运行一个周期平均所需时间 分; A ij ~ 从i 号铲位到j 号卸点最多能同时运行的卡车数 辆; B ij ~ 从i 号铲位到j 号卸点路线上一辆车最多可以运行的次数 次; p i ~ i 号铲位的矿石铁含量。 % p =(30,28,29,32,31,33,32,31,33,31) q j ~ j 号卸点任务需求 吨 q =(1.2,1.3,1.3,1.9,1.3)*10000

数学建模知识竞赛题库

数学建模知识竞赛题库 1.请问计算机中的二进制源于我国古代的哪部经典? D A.《墨经》 B.《诗经》 C.《周书》 D.《周易》 2.世界上面积最大的高原是?D A.青藏高原 B.帕米尔高原 C.黄土高原 D.巴西高原 3.我国海洋国土面积约有多少万平方公里? B A.200 B.300 C.280 D.340 4.世界上面值最高的邮票是匈牙利五百亿彭哥,它的图案是B A.猫 B.飞鸽 C.海鸥 D.鹰 5. 龙虾是我们的一种美食、你知道它体内的血是什么颜色的吗?B A.红色 B.蓝色 C.灰色 D.绿色 6.MATLAB使用三维向量[R G B]来表示一种颜色,则黑色为(D ) A. [1 0 1] B. [1 1 1] C. [0 0 1] D. [0 0 0] 7.秦始皇之后,有几个朝代对长城进行了修葺? A A.7个 B.8个 C.9个 D.10个 8.中国历史上历时最长的朝代是?A A.周朝 B.汉朝 C.唐朝 D.宋朝 9我国第一个获得世界冠军的是谁?C A 吴传玉 B 郑凤荣 C 荣国团 D 陈镜开 10.我国最早在奥运会上获得金牌的是哪位运动员?B A.李宁 B.许海峰 C.高凤莲 D.吴佳怩

11.围棋共有多少个棋子?B A.360 B.361 C.362 D.365 12下列属于物理模型的是:A A水箱中的舰艇 B分子结构图 C火箭模型 D电路图 13名言:生命在于运动是谁说的?C A.车尔尼夫斯基 B.普希金 C.伏尔泰 D.契诃夫 14.饱食后不宜剧烈运动是因为B A.会得阑尾炎 B.有障消化 C.导致神经衰弱 D.呕吐 15、MATLAB软件中,把二维矩阵按一维方式寻址时的寻址访问是按(B)优先的。 A.行 B.列 C.对角线 D.左上角16红军长征中,哪次战役最突出反应毛泽东的军事思想和指挥才?A A.四渡赤水B.抢渡大渡河C.飞夺泸定桥D.直罗镇战役 17色盲患者最普遍的不易分辨的颜色是什么?A A.红绿 B.蓝绿 C.红蓝 D.绿蓝 18下列哪种症状是没有理由遗传的? A.精神分裂症 B.近视 C.糖尿病 D.口吃 19下面哪个变量是正无穷大变量?(A )

全国数学建模大赛题目

2010高教社杯全国大学生数学建模竞赛题目 A题储油罐的变位识别与罐容表标定 通常加油站都有若干个储存燃油的地下储油罐,并且一般都有与之配套的“油位计量管理系统”,采用流量计和油位计来测量进/出油量与罐内油位高度等数据,通过预先标定的罐容表(即罐内油位高度与储油量的对应关系)进行实时计算,以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后,由于地基变形等原因,使罐体的位置会发生纵向倾斜和横向偏转等变化(以下称为变位),从而导致罐容表发生改变。按照有关规定,需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图,其主体为圆柱体,两端为球冠体。图2是其罐体纵向倾斜变位的示意图,图3是罐体横向偏转变位的截面示意图。 请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 (1)为了掌握罐体变位后对罐容表的影响,利用如图4的小椭圆型储油罐(两端平头的椭圆柱体),分别对罐体无变位和倾斜角为α=4.10的纵向变位两种情况做了实验,实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响,并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 (2)对于图1所示的实际储油罐,试建立罐体变位后标定罐容表的数学模型,即罐内储油量与油位高度及变位参数(纵向倾斜角度α和横向偏转角度β)之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据(附件2),根据你们所建立的数学模型确定变位参数,并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。 附件1:小椭圆储油罐的实验数据 附件2:实际储油罐的检测数据 地平线油位探针

数学建模优化问题经典练习

1、高压容器公司制造小、中、大三种尺寸的金属容器,所用资源为金属板、劳 万元,可使用的金属板有500t,劳动力有300人/月,机器有100台/月,此外,不管每种容器制造的数量是多少,都要支付一笔固定的费用:小号为100万元,中号为150万元,大号为200万元,现在要制定一个生产计划,使获得的利润为最大, max=4*x1+5*x2+6*x3-100*y1-150*y2-200*y3; 2*x1+4*x2+8*x3<=500; 2*x1+3*x2+4*x3<=300; 1*x1+2*x2+3*x3<=100; @bin(y1); @bin(y2); @bin(y3); y1+y2+y3>=1; Global optimal solution found. Objective value: 300.0000 Extended solver steps: 0 Total solver iterations: 0 Variable Value Reduced Cost X1 100.0000 0.000000 X2 0.000000 3.000000 X3 0.000000 6.000000 Y1 1.000000 100.0000 Y2 0.000000 150.0000 Y3 0.000000 200.0000 Row Slack or Surplus Dual Price 1 300.0000 1.000000 2 300.0000 0.000000 3 100.0000 0.000000 4 0.000000 4.000000 5 0.000000 0.000000

数学模型经典例题

一、把椅子往地面一放,通常只有三只脚着地,放不稳,然而只需稍挪动几次,就可以使四只脚同时着地放稳了,就四脚连线成长方形的情形建模并加以说明。(15分) 解:一、模型假设: 1. 椅子四只脚一样长,椅脚与地面的接触可以看作一个点,四脚连线呈长方形。 2. 地面高度是连续变化的,沿任何方向都不会出现间断,地面可以看成一张光滑曲面。 3. 地面是相对平坦的,使椅子在任何位置至少有三只脚同时着地。 (3分) 二、建立模型: 以初始位置的中位线为坐标轴建立直角坐标系,用θ表示椅子绕中心O 旋转的角度,椅子的位置可以用θ确定: ()f θ记为A 、B 两点与地面的距离之和 ()g θ记为C 、D 两点与地面的距离之和 由假设3可得,()f θ、()g θ中至少有一个为0。 由假设2知()f θ、()g θ是θ的连续函数。 (3分) 问题归结为: 已知()f θ和()g θ是θ的连续函数,对任意θ, ()()0f g θθ=,且设()()00,00g f =>。证明存在0θ, 使得()()000f g θθ== (3分) 三、模型求解: 令()()()h f θθθ=-g 若()()000f g =,结论成立 若()()000f g 、不同时为,不妨设()()00,00g f =>,椅子旋转()180π或后,AB 与CD 互换,即()()0,0g f ππ>=,则()(0)0,0h h π><。 (3分) 由f g 和的连续性知h 也是连续函数。根据连续函数的基本性质,必存在 ()000θθπ<<使000()0,()()h f g θθθ==即。 最后,因为00()()0f g θθ=,所以00()()0f g θθ==。 (3分) 图 5

多元线性回归 数学建模经典案例

多元线性回归 黄冈职业技术学院数学建模协会胡敏 作业: 在农作物害虫发生趋势的预报研究中,所涉及的5个自变量及因变量的10组观测数据如下,试建立y对x1-x5的回归模型,指出那些变量对y有显著的线性贡献,贡献大小顺序。 x1 x2 x3 x4 x5 y 9.200 2.732 1.471 0.332 1.138 1.155 9.100 3.732 1.820 0.112 0.828 1.146 8.600 4.882 1.872 0.383 2.131 1.841 10.233 3.968 1.587 0.181 1.349 1.356 5.600 3.732 1.841 0.297 1.815 0.863 5.367 4.236 1.873 0.063 1.352 0.903 6.133 3.146 1.987 0.280 1.647 0.114 8.200 4.646 1.615 0.379 4.565 0.898 8.800 4.378 1.543 0.744 2.073 1.930 7.600 3.864 1.599 0.342 2.423 1.104 编写程序如下: data ex; input x1-x5 y@@; cards; 9.200 2.732 1.471 0.332 1.138 1.155 9.100 3.732 1.820 0.112 0.828 1.146 8.600 4.882 1.872 0.383 2.131 1.841 10.233 3.968 1.587 0.181 1.349 1.356 5.600 3.732 1.841 0.297 1.815 0.863 5.367 4.236 1.873 0.063 1.352 0.903 6.133 3.146 1.987 0.280 1.647 0.114 8.200 4.646 1.615 0.379 4.565 0.898 8.800 4.378 1.543 0.744 2.073 1.930 7.600 3.864 1.599 0.342 2.423 1.104 ; proc reg; model y=x1 x2 x3 x4 x5/cli; run; 运行结果如下: (1)回归方程显著性检验. Analysis of Variance Sum of Mean Source DF Squares Square F Value Pr > F Model 5 2.25207 0.45041 11.63 0.0170 Error 4 0.15497 0.03874 Corrected Total 9 2.40704

减速路障间距设计 ;经典数学建模题目分析

组号:702 田宇;孙蕙雯;樊博 校园减速路障间距设计 摘要:减速路障的间距设计合理对于减速带作用的发挥具有重要的意义。本文利用查阅的相关资料,采用Lingo回归分析和最小二乘法,对汽车的加速时加速度和加速时的加速度进行了参数估计。根据题意进行数学建模,建立了汽车在一条具有多个减速带的公路上加速后减速匀速通过减速带的一维直线运动的模型。通过牛顿运动学公式进行了模型求解,最后得出了相邻减速带间的最佳距离。 关键词:减速带间距;一维直线运动模型;最小二乘法

一、问题的提出 1.1 问题的背景 校园、居民小区的道路中间,常常设置用于限制汽车速度的减速带(路障)。减速带使路面稍微拱起以达到车辆减速目的,设置在需要车辆减速慢行的路段和容易引发交通事故的路段,是用于减速机动车、非机动车行使速度的新型交通专用安全设置。减速带很大程度减少了各交通要道口的事故发生,是交通安全的新型专用设施。汽车在行驶中既安全又起到缓冲减速目的,提高交通道口的安全。随着校园车辆的逐渐增多,在校园中合理的设置减速带又成为一个很重要的实际问题。 减速带的使用效果在很大程度上取决于车辆的运行速度和减速带的放置间距间距。因此,为确保限速安全和驾驶人的舒适,合理设定道路的限速具有很重要的意义。 1.2 问题重述 校园道路需要设置路障以限制车速,如果车速不超过40km/h,应该相距多远? 二、问题的分析 2.1 模型预备知识 道路减速带的减速原理:道路减速带的减速是通过影响驾驶员的驾驶心理实现的。当车辆以较高速度进入道路减速带时,剧烈的振动会从轮胎经车身及座椅传递给驾驶员,使驾驶员产生强烈的生理刺激(包括振动刺激和视觉刺激)和心理刺激,从而促使驾驶员主动减速,使车辆以较低的速度通过道路减速带。 2.2问题的分析 1、汽车通过减速带时速度近于零,过减速带后加速。 2、车速达到40km/h时因为前面有下一个减速带而减速,至减速带处车速又近于零。 3、如此循环达到减速目的。

数学建模题目及其答案(疾病诊断)

数学建模疾病的诊断 现要你给出疾病诊断的一种方法。 胃癌患者容易被误诊为萎缩性胃炎患者或非胃病者。从胃癌患者中抽 取5人(编号为1-5),从萎缩性胃炎患者中抽取5人(编号为6-10),以及非胃病者 中抽取5人(编号为11-15),每人化验4项生化指标:血清铜蓝蛋白( X)、 1 蓝色反应( X)、尿吲哚乙酸(3X)、中性硫化物(4X)、测得数据如表1 2 所示: 表1. 从人体中化验出的生化指标

* 根据数据,试给出鉴别胃病的方法。 论文题目:胃病的诊断 摘要 在临床医学中,诊断试验是一种诊断疾病的重要方法。好的诊断试验方法将对临床诊断的正确性和疾病的治疗效果起重要影响。因此,对于不同疾病不断发现新的诊断试验方法是医学进步的重要标志。传统的诊断试验方法有生化检测、DNA检测和影像检测等方法。而本文则通过利用多元统计分析中的判别分析及SPSS软件的辅助较好地解决了临床医学中胃病鉴别的问题。在临床医学上,既提高了临床诊断的正确性,又对疾病的治疗效果起了重要效果,同时也减轻了病人的负担。 判别分析是在分类确定的条件下,根据某一研究对象的各种特征值判别其类型归属问题的一种多变量统计分析方法。 其基本原理是按照一定的判别准则,建立一个或多个判别函数,用研究对象的大量资料确定判别函数中的待定系数,并计算判别指标。 , 首先,由判别分析定义可知,只有当多个总体的特征具有显著的差异时,进行判别分析才有意义,且总体间差异越大,才会使误判率越小。因此在进行判别分析时,有必要对总体多元变量的均值进行是否不等的显著性检验。 其次,利用判别分析中的费歇判别和贝叶斯判别进行判别函数的建立。 最后,利用所建立的判别函数进行回判并测得其误判率,以及对其修正。 本文利用SPSS软件实现了对总体间给类变量的均值是否不等的显著性检验并根据样本建立了相应的费歇判别函数和贝叶斯判别函数,最后进行了回判并测得了误判率,从而获得了在临床诊断中模型,给临床上的诊断试验提供了新方法和新建议。 关键词:判别分析;判别函数;Fisher判别;Bayes判别 一问题的提出 在传统的胃病诊断中,胃癌患者容易被误诊为萎缩性胃炎患者或非胃病患者,为了

经典的数学建模例子1

经典的数学建模例子 一、摘要 SARS SARS就是传染性非典型肺炎,全称严重急性呼吸综合症(Severe Acute Respiratory Syndromes),简称SARS,是一种因感染SARS相关冠状病毒而导致的以发热、干咳、胸闷为主要症状,严重者出现快速进展的呼吸系统衰竭,是一种新的呼吸道传染病,传染性极强、病情进展快速。 当一种传染病流行的时候,会给人们的工作学习带来很大的不变,能有效地进行隔离、预防,会大大减少人员的得病率,当一种传染病开始流行时,在一定的条件下其趋势就像真菌的繁殖曲线,如果能通过计算预测但大概推算出其发病率高峰时期,及时的隔离预防。那会给社会人力带来很大的方便,当年SARS的爆发给我们带来和大的不便和损失,因此本论文就以SARS为例,来研究传染病的传播规律、为预测和控制传染病蔓延创造条件和帮助。 1 二、正文 1、模型的背景问题描述 SARS(Severe Acute Respiratory Syndrome,严重急性呼吸道综合症, 俗称:非典型肺炎)是21世纪第一个在世界范围内传播的传染病。SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响,我们从中得到了许多重要的经验和教训,认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条件的重要性。 要求:(1)建立传染病传播的指数模型,评价其合理性和实用性。 (2)建立一个适合的模型,说明为什么优于问题1中的模型;特别要说明怎样才能 3 建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型,这样做的困难在哪里?对于卫生部门所采取的措施做出评论,如:提前或延后5天采取严格的隔离措施,对疫情传播所造成的影响做出估计。表中提供的数据供参考。 (3)说明建立传染病数学模型的重要性。 2、模型假设 (一)答;

历年全国数学建模试题及其解法归纳

历年全国数学建模试题及解法归纳 赛题解法 93A非线性交调的频率设计拟合、规划 93B足球队排名图论、层次分析、整数规划94A逢山开路图论、插值、动态规划 94B锁具装箱问题图论、组合数学 95A飞行管理问题非线性规划、线性规划 95B天车与冶炼炉的作业调度动态规划、排队论、图论96A最优捕鱼策略微分方程、优化 96B节水洗衣机非线性规划 97A零件的参数设计非线性规划 97B截断切割的最优排列随机模拟、图论 98A一类投资组合问题多目标优化、非线性规划98B灾情巡视的最佳路线图论、组合优化 99A自动化车床管理随机优化、计算机模拟 99B钻井布局0-1规划、图论 00A DNA序列分类模式识别、Fisher判别、人工 神经网络 00B钢管订购和运输组合优化、运输问题 01A血管三维重建曲线拟合、曲面重建

赛题解法 01B 公交车调度问题多目标规划 02A车灯线光源的优化非线性规划 02B彩票问题单目标决策 03A SARS的传播微分方程、差分方程 03B 露天矿生产的车辆安排整数规划、运输问题 04A奥运会临时超市网点设计统计分析、数据处理、优化04B电力市场的输电阻塞管理数据拟合、优化 05A长江水质的评价和预测预测评价、数据处理 05B DVD在线租赁随机规划、整数规划 06A出版社书号问题整数规划、数据处理、优化06B Hiv病毒问题线性规划、回归分析 07A 人口问题微分方程、数据处理、优化07B 公交车问题多目标规划、动态规划、图 论、0-1规划 08A 照相机问题非线性方程组、优化 08B 大学学费问题数据收集和处理、统计分 析、回归分析 2009年A题制动器试验台的控制方法分析工程控制 2009年B题眼科病床的合理安排排队论,优化,仿真,综 合评价 2009年C题卫星监控几何问题,搜集数据

差微分方程 数学建模经典案例

差分方程作业题 黄冈职业技术学院 宋进健 胡敏 熊梦颖 1.一对年轻夫妇准备购买一套住房,但缺少资金近6万元。假设它们每月可有节余900元,且有如下的两种选择: (1)使用银行贷款60000元。月利率0.01,贷款期25年=300个月; (2) 到某借贷公司借贷60000元,月利率0.01,22年还清。只要(i )每半个月还316元,(ii) 预付三个月的款。 你能帮他们做出明智的选择吗? 模型假设: (1)银行及借贷公司在贷款期限内利率不变; (2)不考虑物价变化和经济等因素从而影响利率; (3)银行利息按复利计算且单位时间可任意缩短至时间变量连续性变化 建立模型: 对第一种情况有: 设n 年期贷款月利率为r ,共贷款 元,贷款后第k 个月时欠款余额为 元,月还款m 元。 模型求解: 由MATLAB 得出结果m=631.9345 建立模型: 对第二种情况有: 设n 年期贷款半月利率为r ,共贷款A 0元,贷款后第k 个月时欠款余额为A k 元,半月还款m 元。 模型求解: ()() 011 1,k k k r A A r m k N r +-=+-∈1 0)1()1(300 300 300 -= ?=++r r A A r m N k m r A A k K ∈-+=+,) 1(1 N k m r A A k K ∈-+=+,) 1(1 ()() 011 1,k k k r A A r m k N r +-=+-∈1 0)1()1(528 528 528 -= ?=++r r A A r m A k A 0

由MATLAB 得出结果m= 313.0038 模型分析:由第一种方式计算m=631.9345小于月节余额900元,能够承受月还款;由第二种方式计算m= 313.0038小于借贷公司要求没半个月还款316元,如果按照借贷公司要求则每月还款为632元大于第一种还款方式631.9345元,故选择第一种还款方式。 2. 在一城市的某商业区内,有两家有名的快餐店“肯德基”分店和“麦当劳”分 店。据统计每年“肯德基”保有其上一年老顾客的1/3,而另外的2/3顾客转移到“麦当劳”;每年“麦当劳”保有其上一年的老顾客的1/2,而另外的1/2顾客转移到“肯德基”。 用二维向量X k =[x k y k ]T 表示两个快餐店市场分配的情况,初始的市场分配为X 0 = [200 200]T 如果有矩阵L 存在,使得 X k +1 = LX k ,则称 L 为状态转移矩阵。 (1) 写出X k =[x k y k ]T 和X k+1=[x k +1 y k +1]T 的递推关系式,以及状态转移矩阵L 。 (2) 根据递推关系计算近几年的市场分配情况; 模型假设: (1) 当前的肯德基和麦当劳的市场份额继续不变。 (2) 肯德基和麦当劳不推出优惠活动和新的经营计划。 模型建立: 初始的市场分配数量为:200,2000 0==y x 以一年为一时间段,则某时刻两个快餐店的顾客数量可用向量] ,[1 1y x T X =表 示。用向量] ,[y x X k k T k =表示第K 年两个快餐店顾客数量分布。 ??? ????+ = + = ++x y y y x x k k k k k k 3 22 121311 1 模型求解: 故X k =[x k y k ]T 和X k+1=[x k +1 y k +1]T 的递推关系式为??? ? ?? ? + =+ =++x y y y x x k k k k k k 3 221 21311 1,状 态转移矩阵?????? ? ???? ???=3221213 1 L 由初始数据计算近几年的市场分配情况,MATLAB 程序如下:

一些基本的数学建模示例

1.3 一些基本的数学建模示例 1.3.1椅子的摆放问题 1.3.2 双层玻璃的功效问题 1.3.3 搭积木问题 1.3.4 四足动物的身长和体重关系问题 1.3.5 圆杆堆垛问题 1.3.6 公平的席位分配问题 1.3.7 中国人重姓名问题 1.3.8实物交换问题 椅子能在不平的地面上放稳吗?下面用数学建模的方法解决此问题。 模型准备 仔细分析本问题的实质,发现本问题与椅子腿、地面及椅子腿和地面是否接触有关。如果把椅子腿看成平面上的点,并引入椅子腿和地面距离的函数关系就可以将问题1与平面几何和连续函数联系起来,从而可以用几何知识和连续函数知识来进行数学建模。为讨论问题方便,我们对问题进行简化,先做出如下3个假设: 模型假设 1、椅子的四条腿一样长,椅子脚与地面接触可以视为一个点,四脚连线是正方形(对椅子的假设) 2、地面高度是连续变化的,沿任何方向都不出现间断。(对地面的假设) 3、椅子放在地面上至少有三只脚同时着地,(对椅子和地面之间关系的假设) 根据上述假设做本问题的模型构成: 模型构成Array用变量表示椅子的位置,引入平面图形及坐 标系如图1-1。图中A、B、C、D为椅子的四只脚, 坐标系原点选为椅子中心,坐标轴选为椅子的四 只脚的对角线。于是由假设2,椅子的移动位置 可以由正方形沿坐标原点旋转的角度θ来唯一表 示,而且椅子脚与地面的垂直距离就成为θ的函 数。注意到正方形的中心对称性,可以用椅子的 相对两个脚与地面的距离之和来表示这对应两 个脚与地面的距离关系,这样,用一个函数就可 以描述椅子两个脚是否着地情况。本题引入两个 函数即可以描述椅子四图 1-1

数学建模经典案例最优截断切割问题

建模案例:最优截断切割问题 一、 问 题 从一个长方体中加工出一个已知尺寸、位置预定的长方体(这两个长方体的对应表面是平行的),通常要经过6 次截断切割.设水平切割单位面积的费用是垂直切割单位面积费用的r 倍.且当先后两次垂直切割的平面(不管它们之间是否穿插水平切割)不平行时,因调整刀具需额外费用e.试设计一种安排各面加工次序(称“切割方式”)的方法,使加工费用最少. 二、 假 设 1、假设水平切割单位面积的费用为r ,垂直切割单位面积费用为1; 2、当先后两次垂直切割的平面(不管它们之间是否穿插水平切割)不平行时,调整刀具需额外费用e ; 3、第一次切割前,刀具已经调整完毕,即第一次垂直切割不加入刀具调整费用; 4 、每个待加工长方体都必须经过6次截断切割. 三、 模型的建立与求解 设待加工长方体的左右面、前后面、上下面间的距离分别为 a0、b0 、c0 ,六个切割面分别位于左、右、前、后、上、下,将它们相应编号为M1、M2、M3、M4、M5、M6,这六个面与待加工长方体相应外侧面的边距分别为 u1、u2、u3、u4、u5、u6.这样,一种切割方式就是六个切割面的一个排列,共有P 66720 种切割方式.当考虑到切割费用时,显然有局部优化准则:两个平行待切割面中,边距较大的待切割面总是先加工.

由此准则,只需考虑 P 6 6 222 90 !!! ?? =种切割方式.即在求最少加工费用时,只 需在90个满足准则的切割序列中考虑.不失一般性,设u1≥u2,u3≥u4,u5≥u6,故只考虑M1在M2前、M3在M4前、M5在M6前的切割方式. 1、 e=0 的情况 为简单起见,先考虑e=0 的情况.构造如图9-13的一个有向赋权网络图G(V,E).为了表示切割过程的有向性,在网络图上加上坐标轴x,y,z. 图9-13 G(V,E) 图G(V,E)的含义为: (1)空间网络图中每个结点Vi(xi,yi,zi)表示被切割石材所处的一个状态.顶点坐标xi、yi、zi分别代表石材在左右、前后、上下方向上已被切割的刀数.例如:V24(2,1,2) 表示石材在左右方向上已被切割两刀,前后方向上已被切一刀,上下方向上已被切两刀,即面M1、M2、M3、M5、M6均已被切割.顶点V1(0,0,0) 表示石材的最初待加工状态,顶点V27(2,2,2)表示石材加工完成后的状态.

高教社杯全国大学生数学建模竞赛题目A题

2017年高教社杯全国大学生数学建模竞赛题目 (请先阅读“全国大学生数学建模竞赛论文格式规范”) A题CT系统参数标定及成像 CT(Computed Tomography)可以在不破坏样品的情况下,利用样品对射线能量的吸收特性对生物组织和工程材料的样品进行断层成像,由此获取样品内部的结构信息。一种典型的二维CT系统如图1所示,平行入射的X射线垂直于探测器平面,每个探测器单元看成一个接收点,且等距排列。X射线的发射器和探测器相对位置固定不变,整个发射-接收系统绕某固定的旋转中心逆时针旋转180次。对每一个X射线方向,在具有512个等距单元的探测器上测量经位置固定不动的二维待检测介质吸收衰减后的射线能量,并经过增益等处理后得到180组接收信息。 CT系统安装时往往存在误差,从而影响成像质量,因此需要对安装好的CT系统进行参数标定,即借助于已知结构的样品(称为模板)标定CT系统的参数,并据此对未知结构的样品进行成像。 请建立相应的数学模型和算法,解决以下问题: (1) 在正方形托盘上放置两个均匀固体介质组成的标定模板,模板的几何信息如图2所示,相应的数据文件见附件1,其中每一点的数值反映了该点的吸收强度,这里称为“吸收率”。对应于该模板的接收信息见附件2。请根据这一模板及其接收信息,确定CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离以及该CT系统使用的X射线的180个方向。 (2) 附件3是利用上述CT系统得到的某未知介质的接收信息。利用(1)中得到的标定参数,确定该未知介质在正方形托盘中的位置、几何形状和吸收率等信息。另外,请具体给出图3所给的10个位置处的吸收率,相应的数据文件见附件4。 (3) 附件5是利用上述CT系统得到的另一个未知介质的接收信息。利用(1)中得到的标定参数,给出该未知介质的相关信息。另外,请具体给出图3所给的10个位置处的吸收率。 (4) 分析(1)中参数标定的精度和稳定性。在此基础上自行设计新模板、建立对应的标定模型,以改进标定精度和稳定性,并说明理由。 (1)-(4)中的所有数值结果均保留4位小数。同时提供(2)和(3)重建得到的介质吸收率的数据文件(大小为256×256,格式同附件1,文件名分别为和)

高中常见数学模型案例

高中常见数学模型案例 中华人民共和国教育部2003年4月制定的普通高中《数学课程标准》中明确指出:“数学探究、数学建模、数学文化是贯穿于整个高中数学课程的重要内容”,“数学建模是数学学习的一种新的方式,它为学生提供了自主学习的空间,有助于学生体验数学在解决问题中的价值和作用,体验数学与日常生活和其他学科的联系,体验综合运用知识和方法解决实际问题的过程,增强应用意识;有助于激发学生学习数学的兴趣,发展学生的创新意识和实践能力。”教材中常见模型有如下几种: 一、函数模型 用函数的观点解决实际问题是中学数学中最重要的、最常用的方法。函数模型与方法在处理实际问题中的广泛运用,两个变量或几个变量,凡能找到它们之间的联系,并用数学形式表示出来,建立起一个函数关系(数学模型),然后运用函数的有关知识去解决实际问题,这些都属于函数模型的范畴。 1、正比例、反比例函数问题 例1:某商人购货,进价已按原价a 扣去25%,他希望对货物订一新价,以便按新价让利销售后仍可获得售价25%的纯利,则此商人经营者中货物的件数x 与按新价让利总额y 之间的函数关系是___________。 分析:欲求货物数x 与按新价让利总额y 之间的函数关系式,关键是要弄清原价、进价、新价之间的关系。 若设新价为b ,则售价为b (1-20%),因为原价为a ,所以进价为a (1-25%) 解:依题意,有25.0)2.01()25.01()2.01(?-=---b a b 化简得a b 45=,所以x a bx y ??==2.0452.0,即+∈=N x x a y ,4 2、一次函数问题 例2:某人开汽车以60km/h 的速度从A 地到150km 远处的B 地,在B 地停留1h 后,再以50km/h 的速度返回A 地,把汽车离开A 地的路x (km )表示为时间t (h )的函数,并画出函数的图像。 分析:根据路程=速度×时间,可得出路程x 和时间t 得函数关系式x (t );同样,可列出v(t)的关系式。要注意v(t)是一个矢量,从B 地返回时速度为负值,重点应注意如何画这两个函数的图像,要知道这两个函数所反映的变化关系是不一样的。 解:汽车离开A 地的距离x km 与时间t h 之间的关系式是:?? ???∈--∈∈=]5.6,5.3(),5.3(50150]5.3,5.2(,150]5.2,0[,60t t t t t x ,图略。 速度vkm/h 与时间t h 的函数关系式是:?? ???∈-∈∈=)5.6,5.3[,50)5.3,5.2[,0)5.2,0[,60t t t v ,图略。 3、二次函数问题 例3:有L 米长的钢材,要做成如图所示的窗架,上半部分为半圆,下半部分为六个全等小矩形组成的矩形,试问小矩形的长、宽比为多少时,窗所通过的光线最多,并具体标出窗框面积的最大值。

相关文档
最新文档