(完整版)数学模型姜启源-第三章(第五版)
数学模型第五版
数学建模的能力
想象力
洞察力
判断力
比较广博的数学知识
深入实际调查研究的决心和能力
创新意识
• 如何学习数学建模
学别人的模型学习 分析、改进、推广
做自己的模型实际题目;参加竞赛
学别人的模型
对于案例——椅子能在不平的地面上放稳吗; 在学懂的基础上可以作哪些研究
1 模型假设中哪些条件是本质的, 哪些是非本质的 地面高度连续 是 椅子至少三只脚着地 是
用 x 表示船速;y 表示水速,列出方程:
(x y)30750
x=20
(x y)50750 求解 y =5
答:船速为20km/h
航行问题建立数学模型的基本步骤
• 作出简化假设船速 水速为常数 • 用符号表示有关量x, y分别表示船速和水速 • 用物理定律匀速运动的距离等于速度乘以
时间列出数学式子(二元一次方程) • 求解得到数学解答x=20, y=5
章 13 建模示例之一 包饺子中的数学
14 建模示例之二 路障间距的设计
建
立 数 学
模
15 建模示例之三 椅子能在不平的 地面上放稳吗
16 数学建模的基本方法和步骤 17 数学模型的特点和分类
型 18 怎样学习数学建模——学习课程
和参加竞赛
1 1 从现实对象到数学模型
我们常见的模型
玩具 照片、飞机、火箭模型… ~ 实物模型
结论:在模型假设条件下;将椅子绕中心旋转, 一定能找到四只脚着地的稳定点
1 6 数学建模的基本方法和步骤
数学建模的基本方法
对客观事物特性的认识
机理分析
内部机理的数量规律
白箱
测试分析
对量测数据的统计分析 与数据拟合最好的模型
数学模型姜启源 ppt课件
《数学模型》 姜启源 主编
数学模型
9 五 5-6 6.4种群的相互依存
2
7.1市场经济中的蛛网模型
10 五 5-6 7.2减肥计划-节食与运动
2
8.3层次分析模型
12 五 5-6 8.4效益的合理分配
2
9.2报童的诀窍(讨论课)
13 五 5-6 9.5随机人口模型
2
9.6航空公司的预定票策略
14 五 5-6 10.1牙膏的销售量
数学模型
对于一个现实对象,为了一个特定目的, 根据其内在规律,作出必要的简化假设, 运用适当的数学工具,得到的一个数学结构。
数学
建立数学模型的全过程
建模 (包括表述、求解、解释、检验等)
2020/11/13
12
《数学模型》 姜启源 主编
第一章 建立数学模型
1.2 数学建模的重要意义
• 电子计算机的出现及飞速发展; • 数学以空前的广度和深度向一切领域渗透。
1.3 数学建模示例
1.4 数学建模的方法和步骤
1.5 数学模型的特点和分类
1.6 怎样学习数学建模
2020/11/13
8
《数学模型》 姜启源 主编
第一章 建立数学模型
1.1 从现实对象到数学模型
我们常见的模型
玩具、照片、飞机、火箭模型… … ~ 实物模型
水箱中的舰艇、风洞中的飞机… … ~ 物理模型
《数学模型》 姜启源 主编
数学模型
数学模型
2020/11/13
1
《数学模型》 姜启源 主编
数学模型
课程简介
课程名称 数学模型与数学建模 Mathematical Modeling
先修课程 微积分、线性代数、概率论与数理统计 课程简介
(完整版)姜启源数学模型第五版-第5章
• 建立数学模型描述人口发展规律,是制定 积极、稳妥人口政策的前提.
1. 两个基本的人口模型 2. 用美国人口数据估计参数
年
1790
人口(百万) 3.9
3.9 5.1 6.8 … 245.8 265.4 282.4 2810.4
7.7 9.5 11.7 … 228.3 252.0 275.1 458.2
2. 参数估计
x 300
250 logistic模型 200 (方法一)
x 300
250 logistic模型 200 (方法二)
150
150
100
1. 模型建立 r(x) a bx a = r r(x) r(1 x / xm ) r(0)=r, r(xm)=0 b = r/xm
dx rx(1 x ),
dt
xm
x(0) x0
rx~人口自身增长 (1-x/xm)~资源和环境阻滞人口增长
dx/dt
x
渐近线
xm
S形曲线
xm/2
x增加先快后慢
(百万) (方法一) (方法二) 数模型 (方法一) (方法二)
2010年 308.7
515.0
356.0
314.0
296.8
297.0
误差
66.8% 15.3% 1.7%
-3.9%
-3.8%
2020年 ?
327.8
326.8
模型检验的误差在5%以内,可以接受.
拭目
预测准确性需等2020年美国人口调查结果公布. 以待
(完整版)姜启源数学模型第五版-第6章
等额本息还款模型
x0 ~贷款总额
r ~月利率 n ~贷款期限(月)
xk ~第k月还款后尚欠金额
a~每月还款金额
本月欠额=上月欠额的本息还款金额
xk= xk-1(1+r)a, k=1,2,…, n k=n递推至k=1
xn= x0(1+r)na[1+(1+r)+…+(1+r)n-1]
k 1
2
例2 x0 =100(万元), r=0.0655/12, n=1220=240(月) x1=9625元, x240=4189.41(元), A2=1657729.17(元).
与房贷计算器给出的相同
等额本息与等额本金方式的比较
• 等额本息方式简单,便于安排收支. • 等额本金方式每月还款金额前期高于等额本息方式,
贷款购房——最简 单的差分方程模型
输入必要信息 轻击鼠标即得
单利和复利 两种计算利息的基本方式
单利 ~1万元存5年定期, 年利率4.70.04755)=12375元.
复利 ~1万元存1年定期, 年利率为3%, 到期不取则 自动转存, 5年后本息:10000 (1+0.03)5=11593元.
c 20000 0.025
w 8000100
每周每千克体重消耗 20000/100=200kcal 正常代谢消耗相当弱.
2. 正常代谢情况下的第一阶段计划 • 吸收热量由20000kcal每周减少1000kcal, 直至达到安全下限10000 kcal/周. c(k) 200001000k, k 1,2, ,10 c(10)= 10000 第一阶段需10周 w(k 1) (1 )w(k) (20000 1000 k)
姜启源 数学模型第五版-第1章
1.3
问题
建模示例之一 包饺子中的数学
通常,1kg馅, 1kg面, 包100个饺子. 今天,馅比 1kg多, 1kg面不变, 要把馅包完.
应多包几个(每个小些), 还是少包几个(每个大些)?
分析
直观认识——“大饺子包的馅多”! 但是:“用的面皮也多”!
需要比较:饺子从小变大时馅和面增加的数量关系.
C
C´ B´ B A´
O
A
x
D´
D
A,C 两脚与地面距离之和 ~ f() B,D 两脚与地面距离之和 ~ g()
正方形ABCD 绕O点旋转
模型建立
地面为连续曲面 椅子在任意位置 至少三只脚着地 椅子旋转900, 对 角线AC和BD互换 f() , g()是连续函数 对任意, f(), g() 至少一个为0 g(0)=0,f(0) > 0, f(/2)=0, g(/2)>0.
不平的地面上的椅子, 通常三只脚着地—— 放不稳! 挪动几下,使四只脚着地——椅子放稳!
讨论椅子能放稳的条件.
椅子能在不平的地面上放稳吗
模型假设
四腿一样长,椅脚与地面点接触,四脚连线呈正方形. 地面高度连续变化,可视为数学上的连续曲面. 地面相对平坦,椅子在任意位置至少三只脚着地.
模型建立
椅子位置 利用正方形(椅脚连线)的对称性. 用表示椅子位置. 四只脚着地 椅脚与地面距离为零 距离是的函数. 四个距离 (四只脚) 对称性 两个距离
模 型 构 成
尽量采用简单的数学工具
数学建模的一般步骤 模型 求解 模型 分析 模型 检验 各种数学方法、软件和计算机技术. 如结果的误差分析、统计分析、 模型对数据的稳定性分析. 与实际现象、数据比较, 检验模型的合理性、适用性.
数学模型课后答案姜启源
数学模型课后答案姜启源【篇一:姜启源《数模》习题选解】方案模型构成:以阈值0,1分别标记“不在”和“在”,记第k次渡河前此岸的人阈值为xk,猫阈值为yk,鸡阈值为zk,米阈值为wk,将四维向量sk=(xk,yk,zk,wk)定义为状态,xk,yk,zk,wk=0,1。
安全渡河条件下的状态集合为允许状态集合,记作s。
以穷举法得到s:s={(1,1,1,1),(1,1,1,0),(1,1,0,1),(1,0,1,1),(1,0,1,0),(0,1,0,1),(0,0,1,0),( 0,1,0,0),(0,0,0,1),(0,0,0,0)} 记第k次渡船上四个对象(人、猫、鸡、米)的阈值分别为ak,bk,ck,dk,并将四维向量ek=(ak,bk,ck,dk)定义为决策。
允许决策集合记作e={(a,b,c,d)|0≤b+c+d≤1,a=1,b,c,d=0,1}因为k为奇数时,船从此岸驶向彼岸,k为偶数时船由彼岸驶向此岸,所以,状态sk随决策ek变化的规律是sk+1=sk+(-1)kek该式称状态转移律,该问题就转换成多步决策模型:求决策∈?? ??=1,2,?,?? ,使状态∈??按照转移律,由初始状态s1=(1,1,1,1)经有限步n到达状态sn+1=(0,0,0,0)。
模型求解:本解答试尝用图解法,由于无法利用平面来表达四维坐标系,所以采取其投影即三维空间的方法来构建模型。
把人的阈值xk抽离出来,分别标记0系坐标系(即当xk=0时,(yk,zk,wk)的空间坐标),和1系坐标系,可允许状态点如下标示(红色点):由于a=1是恒成立的,所以,决策是0系坐标系和1系坐标系的点集间的连接,而非任意坐标系内部的连接。
如图1所示,两正方体中心重合,且对应顶点的连线通过中心,称为二合正方体(四维空间不具有包性,即a/b两正方体并没有包含的关系)。
二合正方体的一个顶点为(a,b),称为共顶点,即二合正方体共有8个共顶点。
数学模型第01章第五版ppt课件
3)据连续函数的基本性质, 必存在0 ( 0< 0 < /2) , 使h(0)=0, 即 f(0) = g(0) . 4)因为 f(0) • g(0)=0, 所以 f(0) = g(0) = 0.
结论:在模型假设条件下,将椅子绕中心旋转, 一定能找到四只脚着地的稳定点.
表现特性 建模目的
确定和随机
静态和动态
离散和连续
线性和非线性
描述、优化、预报、决策、…
了解程度 白箱
灰箱
黑箱
1.8 怎样学习数学建模—— 学习课程和参加竞赛
数学建模与其说是一门技术,不如说是一门艺术.
技术大致有章可循. 艺术无法归纳成普遍适用的准则.
• 着重培养数学建模的意识和能力 数学建模的意识 对于日常生活和工作中那些需要 或者可以用数学知识分析、解决的实际问题,能够 敏锐地发现并从建模的角度去积极地思考、研究.
用 x 表示船速,y 表示水速,列出方程:
(x y) 30 750
x=20
( x y) 50 750 求解 y =5
答:船速为20km/h.
航行问题建立数学模型的基本步骤
• 作出简化假设(船速、水速为常数) • 用符号表示有关量(x, y分别表示船速和水速) • 用物理定律(匀速运动的距离等于速度乘以 时间)列出数学式子(二元一次方程) • 求解得到数学解答(x=20, y=5)
1.4 建模示例之二 路障间距的设计
背景 校园、居民小区道路需要限制车速——设置路障 问题 限制车速≤40km/h, 相距多远设置一个路障?
分析 汽车过路障时速度接近零, 过路障后加速.
车速达到40km/h时让司机看到下一路障而 减速, 至路障处车速又接近零. 如此循环以达到限速的目的.
数学模型-第03章(第五版)
存在恰当的x,使f1(x), f2(x)之和最小.
分析
• 关键是对B(t)作出合理的简化假设.
失火时刻t=0, 开始救火时刻t1, 灭火时刻t2, 画出时刻t森林烧毁面积B(t)的大致图形.
B
分析B(t)比较困难, 转而讨论单位时间 烧毁面积 dB/dt (森林烧毁的速度).
第三章
材料强度最大
简单优化模型
利润最高 风险最小
优化——工程技术、经济管理、科学研究中的常见问题. 运输费用最低
用数学建模方法解决优化问题的过程 优化目标与决策 模型假设与建立 数学求解与分析
简单优化模型归结为函数极值问题,用微分法求解. 属于数学规划的优化模型在第四章讨论.
第 三 章 简 单 优 化 模 型
3.2 森林救火
问题
森林失火后,要确定派出消防队员的数量. 队员多,森林损失小,救援费用大; 队员少,森林损失大,救援费用小. 综合考虑损失费和救援费,确定队员数量.
分析
记队员人数x, 失火时刻t=0, 开始救火时刻t1, 灭火时刻t2, 时刻t森林烧毁面积B(t).
• 损失费f1(x)是x的减函数, 由烧毁面积B(t2)决定.
啤酒杯重心s(x)只与质量比a有关 对于每个a, s(x) 有一最小点. a=0.3, x=0.35左右 s最小, 即重心最低.
0.5
s
0.45 a=1 0.4 a=0.5 0.35 a=0.3 0.3
0.25 a=0.1 0.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
建立啤酒杯重心模型一
啤酒杯重心模型一
x
s=s(x) ~ 液面高度x的啤酒杯重心
数学模型姜启源-(第五版)名师公开课获奖课件百校联赛一等奖课件
例2 奶制品旳生产销售计划 在例1基础上深加工
12h 1桶 牛奶 或
3kgA1 1kg 2h, 3元
获利24元/kg 0.8kgB1
获利44元/kg
8h
4kgA2
50桶牛奶, 480h
1kg 2h, 3元
获利16元/kg 0.75kgB2
获利32净利润最大
Objective value:
3460.800
Total solver iterations:
2
Variable
Value Reduced
Cost
X1 0.000000
1.680000
X2 168.0000
0.000000
X3 19.20230
0.000000
X4 0.000000
0.000000
O
c l5
l3 D x1
z=0 z=2400
在B(20,30)点得到最优解.
目的函数和约束条件是线性函数 可行域为直线段围成旳凸多边形 目旳函数旳等值线为直线
最优解一定在凸多边 形旳某个顶点取得.
模型求解
软件实现
LINGO
model: max = 72*x1+64*x2; [milk] x1 + x2<50; [time] 12*x1+8*x2<480; [cpct] 3*x1<100; end
决策 变量
目的 函数
8h
4kg A2
1kg
2h, 3元
出售x1 kg A1, x2 kg A2,
获利16元/kg
0.75kg B2
获利32元/kg
x3 kg B1, x4 kg B2
数学建模模型
超市员工安排及运营问题摘要在一些大型服务机构中,不同的时间段内需要的服务量有着显着的不同,从而主管单位在不同的时段雇佣工作人员的人数往往也不同。
因此对于既要满足需要,又要尽量减少劳务开支是管理者必须思考的决策问题。
本文我院某校内超市员工安排问题为例,据已给定的各个时间段所需的服务员人数和两个班次与休息时间安排表、职员工资及其他给定的限制,建立整数规划优化模型,得出最优安排,使得既满足超市对职工的需要,又使超市的劳务开支最少。
另外本文进一步讨论在已有班次的基础上,对增加更多的班次后的人员安排及劳务支出的变化,以便此超市根据最少的劳务开支做出最优选择。
由问题给出的时间和班次安排表,在8:00——17:00和12:00——21:00中每隔一个小时安排吃饭时间,根据班次安排的人数列出线性不等式,根据月支出来列出目标函数,然后设计线性规划模型,用LINGO.8解出人数和最优劳务支出。
由此解决了本问题要讨论的最少人数和最优劳务支出。
关键词:优化设计,劳务开支,临时员工安排。
一问题重述在一些大型服务机构中,不同的时间段内需要的服务量有显著的不同。
例如,交通管理人员、医院医护人员、宾馆服务人员、超市卖场营销人员等。
在不同的时段劳务需求量不同,主管单位在不同时段雇佣的临时职工数量往往也不同。
因此对于既要满足需要,又要尽量节约劳务开支是管理者必须思考的决策问题。
现就我院校内某超市临时员工的班次安排问题建立一个数学模型来进行优化设计,使其既满足超市的营业需要,又能够使超市的劳务开支最少。
超市的营业时间为11:00到22:OO,根据学生的购买情况,以一小时为一时段,各时段内所需的服务人员数如表1。
此超市员工由临时工和正式员工构成,正式职工两名,主要负责管理工作,每天需要工作8小时,临时工若干名,每天工作4小时。
已知一名正式员工11:00开始上班,工作4小时后休息1小时,而后再工作4小时;另一名正式职工13:00开始上班,工作4小时后休息1小时,而后再工作4小时,工作、休息时间安排如表2。
数学模型第五版姜启源课件
数学模型第五版姜启源课件1. 引言数学模型是一种以数学方法描述、分析和解决实际问题的工具。
它是现代科学、工程和社会学科中不可或缺的一部分。
姜启源的《数学模型》是国内外广泛采用的教材之一,这份课件是对第五版《数学模型》的经典章节进行概要的总结和讲解。
2. 背景与目的数学模型的研究对象可以是自然界的现象、社会经济问题或工程技术等。
通过建立数学模型,我们可以更好地理解问题的本质,并探索解决问题的方法。
数学模型的建立需要一定的理论基础和技巧,本课件旨在帮助读者快速掌握数学模型的基本概念和建模方法。
3. 数学模型的基本概念数学模型是对实际问题进行抽象和描述的数学形式。
它由问题的假设、变量、关系和约束等要素组成。
本部分介绍了数学模型的基本概念,包括:3.1 假设与逼近数学模型的建立需要对实际问题进行适当的假设和逼近。
假设是对问题中不确定因素的简化和规定,而逼近是对问题中不精确因素的近似和描述。
3.2 变量与参数变量是数学模型中描述问题状态的符号,它可以是数值、向量、矩阵等。
参数是数学模型中的固定值,它们可以是已知的或未知的。
3.3 关系与方程关系是数学模型中描述变量之间相互关系的数学表达式。
方程是关系中等号左右两边相等的表达式。
3.4 约束条件与目标函数约束条件是数学模型中描述问题限制条件的不等式或等式。
目标函数是数学模型中描述问题目标的数学表达式。
4. 常见的数学模型本部分介绍了一些常见的数学模型及其应用场景,包括:4.1 线性模型线性模型是最简单的数学模型之一,它的关系和约束条件可以表示为线性方程或线性不等式。
线性模型广泛应用于经济学、管理学、物理学、工程学等领域。
4.2 非线性模型非线性模型是一类不满足线性关系的数学模型。
它的关系和约束条件可以表示为非线性方程或非线性不等式。
非线性模型常用于生物学、化学、地球物理学等领域的研究。
4.3 动态模型动态模型是描述系统随时间变化的数学模型。
它可以采用微分方程、差分方程或积分方程等形式进行建模。
数学模型姜启源答案
数学模型姜启源答案【篇一:姜启源课后习题】xt>第1章建立数学模型1.1 在稳定的椅子问题中,如设椅子的四脚连线呈长方形,结论如何?(稳定的椅子问题见姜启源《数学模型》第6页)1.2 在商人们安全过河问题中,若商人和随从各四人,怎样才能安全过河呢?一般地,有n名商人带n名随从过河,船每次能渡k人过河,试讨论商人们能安全过河时,n与k应满足什么关系。
(商人们安全过河问题见姜启源《数学模型》第7页)1.3 人、狗、鸡、米均要过河,船需要人划,另外至多还能载一物,而当人不在时,狗要吃鸡,鸡要吃米。
问人、狗、鸡、米怎样过河?1.4 有3对夫妻过河,船至多载两人,条件是任一女子不能在其丈夫不在的情况下与其他的男子在一起。
问怎样过河?1.5 如果银行存款年利率为5.5%,问如果要求到2010年本利积累为100000元,那么在1990年应在银行存入多少元?而到2000年的本利积累为多少元?1.6 某城市的logistic模型为dn11dt?25n?25?106n2,如果不考虑该市的流动人口的影响以及非正常死亡。
设该市1990年人口总数为8000000人,试求该市在未来的人口总数。
当t??时发生什么情况。
1.7 假设人口增长服从这样规律:时刻t的人口为x(t),最大允许人口为xm,t到t??t时间内人口数量与xm?x(t)成正比。
试建立模型并求解,作出解的图形并与指数增长模型和阻滞增长模型的结果进行比较。
1.8 一昼夜有多少时刻互换长短针后仍表示一个时间?如何求出这些时间?1.9 你在十层楼上欲乘电梯下楼,如果你想知道需要等待的时间,请问你需要有哪些信息?如果你不愿久等,则需要爬上或爬下几个楼层?1.10 居民的用水来自一个由远处水库供水的水塔,水库的水来自降雨和流入的河流。
水库的水可以通过河床的渗透和水面的蒸发流失。
如果要你建立一个数学模型来预测任何时刻水塔的水位,你需要哪些信息?第2章初等模型2.1 学校共1000名学生,235人住在a宿舍,333人住在b宿舍,432人住在c宿舍。
《数学模型》(第五版)-姜启源 第9章
早产率
H0: q0=q1, H1: q0≠q1
检验结果(α=0.05)
拒绝H0, 接受H1 拒绝H0, 接受H1(t=4.0304) 拒绝H0, 接受H1 接受H0, 拒绝H1(t=0.5663)
• 吸烟孕妇的新生儿体重比不吸烟孕妇的低、且 新生儿体重低的比例高,在统计学上有显著意义.
• 吸烟与不吸烟孕妇孕期和早产率的差别难以肯定 是显著的(若α将接受怀孕期均值相等的假设)
x1
1.0000 0.0809 -0.0534 0.0705 0.0237 -0.0603
x2
1.0000 -0.3510 0.0435 -0.0964 -0.0096
x3
1.0000 -0.0065 0.1473 -0.0678
x4
1.0000 0.4353 0.0175
x5
1.0000 -0.0603
X2
-3.28762 -3.0933 0.0020
X3
-0.00895031 -0.1043 0.9170
X4
1.15497 5.6415 0.0000
X5
0.0498335 1.9910 0.0467
X6
-10 -8 -6 -4 -2
0
2
-8.3939 -8.8248 0.0000
RMSE
16.5
X5
X6
-10 -8
-6
-4
-2
0
2
Coeff. t-stat 0.451168 15.2000 -3.26733 -3.0320 0.104543 1.2775
1.31198 7.1138 0.118183 5.2127
-8.3744 -8.6027
《数学模型》(第五版)-姜启源-第2章
初等模型
• 研究对象的机理比较简单
• 用静态、线性、确定性模型即可达到建模目的
可以利用初等数学方法来构造和求解模型
如果用初等和高等的方法建立的模型,其应用效果
差不多,那么初等模型更高明,也更受欢迎.
尽量采用简单的数学工具来建模
第
二
章
初
等
模
型
双层玻璃窗的功效
划艇比赛的成绩
实物交换
汽车刹车距离与道路通行能力
外
T2
墙
T1 Ta
Ta Tb k Tb T2
Q1 k1
k2
1
d
d
l
T1 T2
k1
l
Q1 k1
, sh , h
d ( s 2)
k2
d
建模 记单层玻璃窗传导的热量Q2
T1 T2
T1 T2
Q1 k1
Q2 k1
d ( s 2)
2d
室
内
T1
2d
Q2
Q1
1
l
, h
Q2 8h 1
d
取 h=l/d=4, 则 Q1/Q2
即双层玻璃窗与同样多材
料的单层玻璃窗相比,可
减少97%的热量损失.
结果分析
Q1/Q2
0.06
0.03
0.02
O
2
4
Q1/Q2所以如此小,是由于层间空气的热传导系
数k2极低, 而这要求空气非常干燥、不流通.
房间通过天花板、墙壁、…损失的热量更多.
vm
vm=vf/2 ~最大流量时的速度
0
km
kj
密度k
0
数学模型课程教学大纲
“数学模型”课程教学大纲课程编码:课程性质:专业必修课学 时: 48学 分:3适用专业:数学与应用数学、信息与计算科学开课学期:4一、课程性质、目的和要求《数学模型》是应用数学和信息与计算科学专业的专业课。
它是研究如何将数学方法和计算机知识结合起来用于解决实际问题的一门边缘交叉学科,是集经典数学、现代数学和实际问题为一体的一门新型课程,是应用数学解决实际问题的重要手段和途径。
《数学模型》课程的意义在于:提高学生的数学素质和应用数学知识解决实际问题的能力,大力培养应用型人才。
本课程是沟通实际问题与数学工具之间联系的必不可少的桥梁。
是一门充分应用其它各数学分支的应用类课程,其主要任务不是“学数学”,而是学着“用数学”,是为培养善于运用数学知识建立实际问题的数学模型,从而善于解决实际问题的应用型数学人材服务的。
通过本课程的学习,使学生较为系统的获得利用数学工具建立数学模型的基本知识、基本技能与常用技巧,培养学生的抽象概括问题的能力,用数学方法和思想进行综合应用与分析问题的能力,并着力导引实践—理论—实践的认识过程,培养学生辩证唯物主义的世界观。
本课程的要求是:通过具体实例的引入使学生掌握数学建模基本思想、基本方法、基本类型,学会进行科学研究的一般过程,并能进入一个实际操作的状态。
通过数学模型有关概念、特征的学习和数学模型应用实例的介绍,要求学生加强数学推导计算和简化分析能力,熟练运用计算机能力,提高联想、洞察能力,综合分析能力和应用数学方法解决实际问题的能力。
二、教学内容、要点和课时安排第一章数学模型简介(6学时)教学目的:了解数学模型的意义及分类;知道常见的数学模型分类;理解建立数学模型的方法及步骤;能熟悉MATLAB运行环境,熟练进行操作.教学重点与难点:建立数学模型的基本方法和步骤;MATLAB的使用主讲内容:1.1数学建模初步⑴什么是数学模型⑵几个数学建模实例⑶数学建模的基本方法和步骤1.2 椅子摆放的问题1.3如何预报人口增长问题实验内容:实验1:MATLAB基本使用方法练习(4)思考题:1、数学建模中各个步骤的作用分别是什么?2、在Census人口模型中,建立一个分段的指数增长模型。
姜启源数学模型课后答案(3版)
《数学模型》作业解答第二章(1)(2008年9月16日)1. 学校共1000名学生,235人住在A 宿舍,333人住在B 宿舍,432人住在C 宿舍.学生们要组织一个10人的委员会,试用下列办法分配各宿舍的委员数:(1). 按比例分配取整数的名额后,剩下的名额按惯例分给小数部分较大者; (2). §1中的Q 值方法;(3).d ’Hondt 方法:将A 、B 、C 各宿舍的人数用正整数n=1,2,3,……相除,其商数如下表:将所得商数从大到小取前10个(10为席位数),在数字下标以横线,表中A 、B 、C 行有横线的数分别为2,3,5,这就是3个宿舍分配的席位.你能解释这种方法的道理吗?如果委员会从10个人增至15人,用以上3种方法再分配名额,将3种方法两次分配的结果列表比较.解:先考虑N=10的分配方案,,432 ,333 ,235321===p p p ∑==31.1000i ip方法一(按比例分配) ,35.23111==∑=i ipNp q ,33.33122==∑=i ipNp q 32.43133==∑=i ipNp q分配结果为: 4 ,3 ,3321===n n n 方法二(Q 值方法)9个席位的分配结果(可用按比例分配)为:4 ,3 ,2321===n n n第10个席位:计算Q 值为,17.92043223521=⨯=Q ,75.92404333322=⨯=Q 2.93315443223=⨯=Q3Q 最大,第10个席位应给C.分配结果为 5 ,3 ,2321===n n n方法三(d ’Hondt 方法)此方法的分配结果为:5 ,3 ,2321===n n n此方法的道理是:记i p 和i n 为各宿舍的人数和席位(i=1,2,3代表A 、B 、C 宿舍).iin p 是每席位代表的人数,取,,2,1 =i n 从而得到的i i n p 中选较大者,可使对所有的,i ii n p尽量接近.再考虑15=N 的分配方案,类似地可得名额分配结果.现将3种方法两次分配的结果列表如下:2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本.考虑t 到t t ∆+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得⎰⎰+=ntdn wkn r k vdt 0)(2π)22 2n wk k(r n πvt +=∴ .2 22n vk w n v rk t ππ+=∴第二章(2)(2008年10月9日)15.速度为v 的风吹在迎风面积为s 的风车上,空气密度是ρ ,用量纲分析方法确定风车获得的功率P 与v 、S 、ρ的关系.解: 设P 、v 、S 、ρ的关系为0),,,(=ρs v P f , 其量纲表达式为: [P]=32-T ML , [v ]=1-LT ,[s ]=2L ,[ρ]=3-ML ,这里T M L ,,是基本量纲.量纲矩阵为:A=)⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---ρ()()()()()()(001310013212s v P T M L齐次线性方程组为:⎪⎩⎪⎨⎧=--=+=-++03032221414321y y y y y y y y 它的基本解为)1,1,3,1(-=y由量纲i P 定理得 1131ρπs v P -=, 113ρλs v P =∴ , 其中λ是无量纲常数. 16.雨滴的速度v 与空气密度ρ、粘滞系数μ和重力加速度g 有关,其中粘滞系数的定义是:运动物体在流体中受的摩擦力与速度梯度和接触面积的乘积成正比,比例系数为粘滞系数,用量纲分析方法给出速度v 的表达式.解:设v ,ρ,μ,g 的关系为(f v ,ρ,μ,g )=0.其量纲表达式为[v ]=LM 0T -1,[ρ]=L -3MT 0,[μ]=MLT -2(LT -1L -1)-1L -2=MLL -2T -2T=L -1MT -1,[g ]=LM 0T -2,其中L ,M ,T 是基本量纲.量纲矩阵为A=)()()()()()()(210101101131g v T M L μρ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----- 齐次线性方程组Ay=0 ,即⎪⎩⎪⎨⎧==+=+02y -y - y -0y y 0y y -3y -y 431324321 的基本解为y=(-3 ,-1 ,1 ,1)由量纲i P 定理 得 g v μρπ13--=. 3ρμλgv =∴,其中λ是无量纲常数. 16*.雨滴的速度v 与空气密度ρ、粘滞系数μ、特征尺寸γ和重力加速度g 有关,其中粘滞系数的定义是:运动物体在流体中受的摩擦力与速度梯度和接触面积的乘积成正比,比例系数为粘滞系数,用量纲分析方法给出速度v 的表达式.解:设v ,ρ,μ,γ,g 的关系为0),,,,(=g v f μργ.其量纲表达式为[v ]=LM 0T -1,[ρ]=L -3MT 0,[μ]=MLT -2(LT -1L -1)-1L -2=MLL -2T -2T=L -1MT -1,[γ]=LM 0T 0 ,[g ]=LM 0T -2其中L ,M ,T 是基本量纲. 量纲矩阵为A=)()()()()()()()(210010110011311g v T M L μργ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----- 齐次线性方程组Ay=0 即⎪⎩⎪⎨⎧=---=+=+--+020035414354321y y y y y y y y y y 的基本解为⎪⎩⎪⎨⎧---=--=)21,1,1,23,0()21,0,0,21,1(21y y得到两个相互独立的无量纲量⎩⎨⎧==-----2/112/322/12/11g g v μργπγπ 即 1212/12/31,--==πμργπγg g v . 由0),(21=Φππ , 得 )(121-=πϕπ ∴ )(12/12/3-=μργϕγυg g , 其中ϕ是未定函数.20.考察阻尼摆的周期,即在单摆运动中考虑阻力,并设阻力与摆的速度成正比.给出周期的表达式,然后讨论物理模拟的比例模型,即怎样由模型摆的周期计算原型摆的周期. 解:设阻尼摆周期t ,摆长l , 质量m ,重力加速度g ,阻力系数k 的关系为0),,,,(=k g m l t f其量纲表达式为:112120000000)(]][[][,][,][,][,][-----======LT MLT v f k T LM g MT L m T LM l T M L t 10-=MT L , 其中L ,M ,T 是基本量纲.量纲矩阵为A=)()()()()()()()(120011010001010k g m l t T M L ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-- 齐次线性方程组⎪⎩⎪⎨⎧=--=+=+02005415342y y y y y y y 的基本解为⎪⎩⎪⎨⎧--=-=)1,21,1,21,0()0,21,0,21,1(21Y Y 得到两个相互独立的无量纲量∴g l t =1π, )(21πϕπ=, 2/12/12mg kl =π ∴)(2/12/1mg kl g l t ϕ=,其中ϕ是未定函数 . 考虑物理模拟的比例模型,设g 和k 不变,记模型和原型摆的周期、摆长、质量分别为t ,'t ;l ,'l ;m ,'m . 又)(2/12/1gm l k g l t '''='ϕ 当无量纲量l l mm '='时, 就有 ll l g g l tt '=⋅'='. 《数学模型》作业解答第三章1(2008年10月14日)1. 在3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货批量.证明在不允许缺货模型中结果与原来的一样,而在允许缺货模型中最优订货周期和订货批量都比原来结果减少.⎩⎨⎧==---22/112/112/12/1ππk g m l g tl解:设购买单位重量货物的费用为k ,其它假设及符号约定同课本.01 对于不允许缺货模型,每天平均费用为:kr rT c T c T C ++=2)(212221r c Tc dT dC+-= 令0=dTdC, 解得 r c c T 21*2= 由rT Q = , 得212c rc rT Q ==**与不考虑购货费的结果比较,T、Q的最优结果没有变.02 对于允许缺货模型,每天平均费用为:⎥⎦⎤⎢⎣⎡+-++=kQ Q rT r c r Q c c T Q T C 23221)(221),(2223322221222T kQ rT Q c r c rT Q c T c T C--+--=∂∂Tk rT Q c c rT Qc Q C ++-=∂∂332 令⎪⎪⎩⎪⎪⎨⎧=∂∂=∂∂00QCTC, 得到驻点:⎪⎪⎩⎪⎪⎨⎧+-+-+=-+=**323222233232132233221)(22c c krc c c r k c c c c c r c Q c c k c c c rc c T与不考虑购货费的结果比较,T、Q的最优结果减少.2.建立不允许缺货的生产销售存贮模型.设生产速率为常数k ,销售速率为常数r ,r k >.在每个生产周期T内,开始的一段时间()00T t <<一边生产一边销售,后来的一段时间)(0T t T <<只销售不生产,画出贮存量)(t g 的图形.设每次生产准备费为1c ,单位时间每件产品贮存费为2c ,以总费用最小为目标确定最优生产周期,讨论r k >>和r k ≈的情况.解:由题意可得贮存量)(t g 的图形如下:贮存费为 ∑⎰=→∆⋅-==∆i Ti i t TT r k c dt t g c t g c 1022022))()(limξ又 )()(00T T r T r k -=- ∴ T k r T =0 , ∴ 贮存费变为 kTT r k r c 2)(2⋅-=于是不允许缺货的情况下,生产销售的总费用(单位时间内)为kTr k r c T c kT T r k r c T c T C 2)(2)()(21221-+=-+=k r k r c Tc dT dC 2)(221-+-=. 0=dT dC令, 得)(221r k r c k c T -=* 易得函数处在*T T C )(取得最小值,即最优周期为: )(221r k r c kc T -=*rc c ,Tr k 212≈>>*时当 . 相当于不考虑生产的情况. ∞→≈*,Tr k 时当 . 此时产量与销量相抵消,无法形成贮存量.第三章2(2008年10月16日)3.在3.3节森林救火模型中,如果考虑消防队员的灭火速度λ与开始救火时的火势b 有关,试假设一个合理的函数关系,重新求解模型.解:考虑灭火速度λ与火势b 有关,可知火势b 越大,灭火速度λ将减小,我们作如下假设: 1)(+=b kb λ, 分母∞→→+λ时是防止中的011b b 而加的. 总费用函数()xc b kx b x t c b kx b t c t c x C 3122121211)1()(2)1(2+--++--++=βββββββ最优解为 []k b k c b b b c kbc x ββ)1(2)1()1(223221+++++=5.在考虑最优价格问题时设销售期为T ,由于商品的损耗,成本q 随时间增长,设t q t q β+=0)(,为增长率β.又设单位时间的销售量为)(为价格p bp a x -=.今将销售期分为T t TT t <<<<220和两段,每段的价格固定,记作21,p p .求21,p p 的最优值,使销售期内的总利润最大.如果要求销售期T 内的总售量为0Q ,再求21,p p 的最优值. 解:按分段价格,单位时间内的销售量为⎪⎩⎪⎨⎧<<-<<-=T t T bp a T t bp a x 2,20,21又 t q t q β+=0)(.于是总利润为[][]⎰⎰--+--=22221121)()()()(),(TTT dt bp a t q p dt bp a t q p p p=22)(022)(20222011T Tt t q t p bp a T t t q t p bp a ⎥⎦⎤⎢⎣⎡---+⎥⎦⎤⎢⎣⎡---ββ=)8322)(()822)((20222011T t q T p bp a T T q T p bp a ββ---+--- )(2)822(12011bp a T T T q T p b p -+---=∂∂β)(2)8322(22022bp a T T t q T p b p -+---=∂∂β 0,021=∂∂=∂∂p p 令, 得到最优价格为: ⎪⎪⎩⎪⎪⎨⎧⎥⎦⎤⎢⎣⎡++=⎥⎦⎤⎢⎣⎡++=)43(21)4(210201T q b a b p T q b a b p ββ 在销售期T 内的总销量为⎰⎰+-=-+-=20221210)(2)()(T TT p p bTaT dt bp a dt bp a Q 于是得到如下极值问题:)8322)(()822)((),(m ax 2022201121T t q T p bp a T T q T p bp a p p ββ---+---=t s . 021)(2Q p p bTaT =+-利用拉格朗日乘数法,解得:⎪⎩⎪⎨⎧+-=--=880201TbT Q b a p T bT Q b a p ββ 即为21,p p 的最优值.第三章3(2008年10月21日)6. 某厂每天需要角钢100吨,不允许缺货.目前每30天定购一次,每次定购的费用为2500元.每天每吨角钢的贮存费为0.18元.假设当贮存量降到零时订货立即到达.问是否应改变订货策略?改变后能节约多少费用?解:已知:每天角钢的需要量r=100(吨);每次订货费1c =2500(元); 每天每吨角钢的贮存费2c =0.18(元).又现在的订货周期T 0=30(天)根据不允许缺货的贮存模型:kr rT c T c T C ++=2121)( 得:k T TT C 10092500)(++=令0=dTdC, 解得:35092500*==T 由实际意义知:当350*=T (即订货周期为350)时,总费用将最小. 又k T C 10035095025003)(*+⨯+⨯==300+100kk T C 100309302500)(0+⨯+==353.33+100k)(0T C -)(*T C =(353.33+100k )-(300+100k )32=53.33.故应改变订货策略.改变后的订货策略(周期)为T *=350,能节约费用约53.33元.《数学模型》作业解答第四章(2008年10月28日)1. 某厂生产甲、乙两种产品,一件甲产品用A 原料1千克, B 原料5千克;一件乙产品用A 原料2千克,B 原料4千克.现有A 原料20千克, B 原料70千克.甲、乙产品每件售价分别为20元和30元.问如何安排生产使收入最大? 解:设安排生产甲产品x 件,乙产品y 件,相应的利润为S 则此问题的数学模型为:max S=20x+30ys.t. ⎪⎩⎪⎨⎧∈≥≤+≤+Z y x y x y x y x ,,0,7045202这是一个整线性规划问题,现用图解法进行求解可行域为:由直线1l :x+2y=20, 2l :5x+4y =702l以及x=0,y=0组成的凸四边形区域.925002+-=TdT dC直线l :20x+30y=c 在可行域内平行移动.易知:当l 过1l 与2l 的交点时, x S 取最大值.由⎩⎨⎧=+=+7045202y x y x 解得⎩⎨⎧==510y x此时 m ax S =2053010⨯+⨯=350(元)2. 某厂拟用集装箱托运甲乙两种货物,每箱的体积、重量以及可获利润如下表:已知这两种货物托运所受限制是体积不超过24立方米,重量不超过13百斤.试问这两种货物各托运多少箱,使得所获利润最大,并求出最大利润.解:设甲货物、乙货物的托运箱数分别为1x ,2x ,所获利润为z .则问题的数学模型可表示为211020 m ax x x z +=⎪⎩⎪⎨⎧∈≥≤+≤+Z y x x x x x x x st ,,0,13522445212121这是一个整线性规划问题. 用图解法求解. 可行域为:由直线2445:211=+x x l1352:212=+x x l 及0,021==x x 组成直线 c x x l =+211020:在此凸四边形区域内平行移动.2ll1x1l2x易知:当l 过l 1与l 2的交点时,z 取最大值由⎩⎨⎧=+=+135224452121x x x x 解得 ⎩⎨⎧==1421x x90110420max =⨯+⨯=z .3.某微波炉生产企业计划在下季度生产甲、乙两种型号的微波炉.已知每台甲型、乙型微波炉的销售利润分别为3和2个单位.而生产一台甲型、乙型微波炉所耗原料分别为2和3个单位,所需工时分别为4和2个单位.若允许使用原料为100个单位,工时为120个单位,且甲型、乙型微波炉产量分别不低于6台和12台.试建立一个数学模型,确定生产甲型、乙型微波炉的台数,使获利润最大.并求出最大利润.解:设安排生产甲型微波炉x 件,乙型微波炉y 件,相应的利润为S. 则此问题的数学模型为:max S=3x +2ys.t. ⎪⎩⎪⎨⎧∈≥≥≤+≤+Z y x y x y x y x ,,12,61202410032这是一个整线性规划问题 用图解法进行求解可行域为:由直线1l :2x+3y=100, 2l :4x+2y =120 及x=6,y=12组成的凸四边形区域.直线l :3x+2y=c 在此凸四边形区域内平行移动. 易知:当l 过1l 与2l 的交点时, S 取最大值.由⎩⎨⎧=+=+1202410032y x y x 解得⎩⎨⎧==2020y x .m ax S =320220⨯+⨯=100.《数学模型》作业解答第五章1(2008年11月12日)1.对于5.1节传染病的SIR 模型,证明: (1)若处最大先增加,在则σσ1)(,10=s t i s ,然后减少并趋于零;)(t s 单调减少至.∞s(2).)()(,10∞s t s t i s 单调减少至单调减少并趋于零,则若σ解:传染病的SIR 模型(14)可写成⎪⎩⎪⎨⎧-=-=i s dtds s i dt diλσμ)1(.)(lim 0.(t) .)( .0,t 存在而单调减少知由∞∞→=∴≥-=s t s s t s dtdsi s dt ds λ.)(∞s t s 单调减少至故(1).s s(t) .s(t) .100≤∴单调减少由若σs;)(,0 .01,10单调增加时当t i dtdis s s ∴-σσ.)(,0.01,1单调减少时当t i dtdis s ∴-σσ .0)(lim .0)18(t ==∞→∞t i i 即式知又由书上.)( .0,1m i t i dtdis 达到最大值时当∴==σ(2)().0 0.1-s ,1,10 dtdit s s σσσ从而则若()().0.0lim ==∴∞∞→i t i t i t 即单调减少且4.在5.3节正规战争模型(3)中,设乙方与甲方战斗有效系数之比为.4=ba初始兵力00y x 与相同.(1) 问乙方取胜时的剩余兵力是多少,乙方取胜的时间如何确定.(2) 若甲方在战斗开始后有后备部队以不变的速率r 增援,重新建立模型,讨论如何判断双方的胜负.解:用()()t y t x ,表示甲、乙交战双方时刻t 的士兵人数,则正规战争模型可近似表示为:()()()⎪⎪⎩⎪⎪⎨⎧==-=-=000,01 ,yy x x bx dtdyay dt dx现求(1)的解: (1)的系数矩阵为⎥⎦⎤⎢⎣⎡--=00b a Aab ab b aA E ±=∴=-==-1,22 .0λλλλλ ⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫⎝⎛-1212,21,对应的特征向量分别为λλ ()()()tab t ab eC e C t y t x -⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫⎝⎛-=⎪⎪⎭⎫ ⎝⎛∴1212121的通解为.再由初始条件,得()()2 220000 tab tab e y x ey x t x -⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛-=又由().1aybx dx dy =可得其解为 ()3 ,202022 bx ay k k bx ay -==-而(1) ()().231000202011y a b y a bx ay ak t y t x =-=-===时,当 即乙方取胜时的剩余兵力数为.230y 又令().0222,01100001=-⎪⎭⎫⎝⎛++⎪⎭⎫⎝⎛-=t ab t ab e y x e y x t x )得由(注意到000020022,1x y y x ey x t ab -+==得. .43ln ,3121bt et ab =∴=∴ (2) 若甲方在战斗开始后有后备部队以不变的速率r 增援.则()()⎪⎪⎩⎪⎪⎨⎧==-=+-=000,)0(4 yy x x bx dtdyr ay dt dx().,4rdy aydy bxdx bxray dy dx -=-+-=即得由 相轨线为,222k bx ry ay =-- .222220.020k a r bx a r y a bx ry ay k =--⎪⎭⎫ ⎝⎛---=或 此相轨线比书图11中的轨线上移了.a r 乙方取胜的条件为.,0222020a r x a b a r y k +⎪⎭⎫ ⎝⎛- 亦即 第五章2(2008年11月14日)6. 模仿5.4节建立的二室模型来建立一室模型(只有中心室),在快速静脉注射、恒速静脉滴注(持续时间为τ)和口服或肌肉注射3种给药方式下求解血药浓度,并画出血药浓度曲线的图形.解: 设给药速率为(),0t f ()()()()().,,0/t VC t x t f t kx t x k ==+则排除速率为常数(1)快速静脉注射: 设给药量为,0D 则()()().,0,0000t k e VDt C V D C t f -===解得 (2)恒速静脉滴注(持续时间为τ): 设滴注速率为()(),00,000==C k t f k ,则解得()()()()⎪⎩⎪⎨⎧-≤≤-=----τττ t e e Vkk t e Vkk t C t k kt kt,10 ,10(3) 口服或肌肉注射: ()(),解得)式节(见134.5010010tk eD k t f -=()()()⎪⎪⎩⎪⎪⎨⎧=≠--=---010101001 ,,01k k te V kD k k e e k k V D k t C kt t k kt3种情况下的血药浓度曲线如下:第五章3(2008年11月18日)8. 在5.5节香烟过滤嘴模型中,(1) 设3.0,/50,08.0,02.0,20,80,80021=======a s mm b mm l mm l mg M νβ求./21Q Q Q 和(2) 若有一支不带过滤嘴的香烟,参数同上,比较全部吸完和只吸到1l 处的情况下,进入人体毒物量的区别.解)(857563.229102.07.050103.01508002.07.0502008.0/01/2毫克≈⎪⎪⎭⎫ ⎝⎛-⨯⨯⨯=⎪⎪⎭⎫ ⎝⎛-=⨯⨯-⨯---e e e eba v aw Q v bl a vl β ()10/10==l M w 其中,()()97628571.0502002.008.0212===⨯----ee Q Q vl b β(2) 对于一支不带过滤嘴的香烟,全部吸完的毒物量为⎪⎪⎭⎫⎝⎛-=-vbl a e b a v aw Q '103‘ 只吸到1l 处就扔掉的情况下的毒物量为⎪⎪⎭⎫ ⎝⎛-=--vbl a v ble e b a v aw Q 1'21'04 .256531719.1110096.0032.0012.004.0508002.03.0508002.05010002.03.05010002.043111'1'≈--=--=--=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-=⨯⨯⨯⨯⨯⨯--e e e e e e e e e e e e e e e e Q Q v abl v bl v abl v bl v bl a v bl v bl a vbl 44.235,84.29543≈≈ QQ4.在5.3节正规战争模型(3)中,设乙方与甲方战斗有效系数之比为.4=ba初始兵力00y x 与相同.(1) 问乙方取胜时的剩余兵力是多少,乙方取胜的时间如何确定.(2) 若甲方在战斗开始后有后备部队以不变的速率r 增援,重新建立模型,讨论如何判断双方的胜负.解:用()()t y t x ,表示甲、乙交战双方时刻t 的士兵人数,则正规战争模型可近似表示为:()()()⎪⎪⎩⎪⎪⎨⎧==-=-=000,01 ,yy x x bx dtdyay dt dx现求(1)的解: (1)的系数矩阵为⎥⎦⎤⎢⎣⎡--=00b a A ab ab b aA E ±=∴=-==-1,22 .0λλλλλ⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫⎝⎛-1212,21,对应的特征向量分别为λλ ()()()tab t ab eC e C t y t x -⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛∴1212121的通解为.再由初始条件,得()()2 220000 tab tab e y x ey x t x -⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛-=又由().1aybx dx dy =可得其解为 ()3 ,202022 bx ay k k bx ay -==-而(1) ()().231000202011y a b y a bx ay ak t y t x =-=-===时,当 即乙方取胜时的剩余兵力数为.230y 又令().0222,01100001=-⎪⎭⎫⎝⎛++⎪⎭⎫⎝⎛-=t ab t ab e y x e y x t x )得由(注意到000020022,1x y y x ey x t ab -+==得. .43ln ,3121bt et ab =∴=∴ (2) 若甲方在战斗开始后有后备部队以不变的速率r 增援.则()()⎪⎪⎩⎪⎪⎨⎧==-=+-=000,)0(4 yy x x bx dtdyr ay dt dx().,4rdy aydy bxdx bxr ay dy dx -=-+-=即得由 相轨线为,222k bx ry ay =-- .222220.020k a r bx a r y a bx ry ay k =--⎪⎭⎫ ⎝⎛---=或 此相轨线比书图11中的轨线上移了.a r 乙方取胜的条件为.,0222020a r x a b a r y k +⎪⎭⎫ ⎝⎛- 亦即《数学模型》作业解答第六章(2008年11月20日)1.在6.1节捕鱼模型中,如果渔场鱼量的自然增长仍服从Logistic 规律,而单位时间捕捞量为常数h .(1)分别就4/rN h >,4/rN h <,4/rN h =这3种情况讨论渔场鱼量方程的平衡点及其稳定状况.(2)如何获得最大持续产量,其结果与6.1节的产量模型有何不同.解:设时刻t 的渔场中鱼的数量为()t x ,则由题设条件知:()t x 变化规律的数学模型为h Nxrx dt t dx --=)1()( 记h Nxrx x F --=)1()( (1).讨论渔场鱼量的平衡点及其稳定性: 由()0=x F ,得0)1(=--h Nxrx . 即()102=+-h rx x Nr )4(42Nhr r N rh r -=-=∆ , (1)的解为:2412,1N rNhN x -±=①当4/rN h >,0<∆,(1)无实根,此时无平衡点;②当4/rN h =,0=∆,(1)有两个相等的实根,平衡点为20N x =. Nrxr N rx N x r x F 2)1()('-=--=,0)(0'=x F 不能断定其稳定性. 但0x x ∀ 及0x x 均有04)1()( rNN x rx x F --= ,即0 dtdx .∴0x 不稳定;③当4/rN h <,0>∆时,得到两个平衡点:2411N rNhN x --=, 2412N rNh N x -+=易知:21N x <, 22N x > ,0)(1'>x F ,0)(2'<x F ∴平衡点1x 不稳定,平衡点2x 稳定.(2)最大持续产量的数学模型为⎩⎨⎧=0)(..max x F t s h 即 )1(max Nxrx h -=,易得 2*0N x = 此时 4rN h =, 但2*0N x =这个平衡点不稳定.这是与6.1节的产量模型不同之处.要获得最大持续产量,应使渔场鱼量2Nx >,且尽量接近2N ,但不能等于2N .2.与Logistic 模型不同的另一种描述种群增长规律的是Gompertz 模型:()xNrx t x ln '=.其中r 和N 的意义与Logistic 模型相同.设渔场鱼量的自然增长服从这个模型,且单位时间捕捞量为Ex h =.讨论渔场鱼量的平衡点及其稳定性,求最大持续产量m h 及获得最大产量的捕捞强度m E 和渔场鱼量水平*0x .解:()t x 变化规律的数学模型为()Ex xNrx dt t dx -=ln 记 Ex xNrx x F -=ln)( ① 令()0=x F ,得0ln =-Ex xNrx ∴r ENe x -=0,01=x .∴平衡点为1,0x x . 又 ()E r xNr x F --=ln',()()∞=<-=1'0',0x F r x F . ∴ 平衡点o x 是稳定的,而平衡点1x 不稳定.②最大持续产量的数学模型为:⎪⎩⎪⎨⎧≠=-=.0,0ln ..max x Ex x N rx t s Ex h 由前面的结果可得 rE ENeh -=r Er Ee r EN Ne dE dh ---=,令.0=dEdh 得最大产量的捕捞强度r E m =.从而得到最大持续产量e rN h m /=,此时渔场鱼量水平eNx =*0. 3.设某渔场鱼量)(t x (时刻t 渔场中鱼的数量)的自然增长规律为:)1()(Nxrx dt t dx -= 其中r 为固有增长率,`N 为环境容许的最大鱼量. 而单位时间捕捞量为常数h .10.求渔场鱼量的平衡点,并讨论其稳定性;20.试确定捕捞强度m E ,使渔场单位时间内具有最大持续产量m Q ,求此时渔场鱼量水平*0x . 解:10.)(t x 变化规律的数学模型为h Nxrx dt t dx --=)1()( 记h N x rx x f --=)1()(,令 0)1(=--h N x rx ,即02=+-h rx x Nr ----(1))4(42Nhr r N rh r -=-=∆ , (1)的解为:2412,1N rNhN x -±=① 当0 ∆时,(1)无实根,此时无平衡点;Ex()x f② 当0=∆时,(1)有两个相等的实根,平衡点为20Nx =. Nrx r N rx N x r x f 2)1()('-=--= ,0)(0'=x f 不能断定其稳定性. 但0x x ∀ 及0x x 均有04)1()( rN N x rx x f --= ,即0 dt dx∴0x 不稳定;③ 当0 ∆时,得到两个平衡点:2411rNhN N x --=, 2412rNh N N x -+=易知 21N x, 22N x ∴0)('1 x f , 0)('2 x f ∴平衡点1x 不稳定 ,平衡点2x 稳定.20.最大持续产量的数学模型为: ⎩⎨⎧=0)(..max x f t s h即 )1(max N x rx h -=, 易得 2*0N x = 此时 4rN h =,但2*0N x =这个平衡点不稳定. 要获得最大持续产量,应使渔场鱼量2N x ,且尽量接近2N ,但不能等于2N.《数学模型》第七章作业(2008年12月4日)1.对于7.1节蛛网模型讨论下列问题:(1)因为一个时段上市的商品不能立即售完,其数量也会影响到下一时段的价格,所以第1+k 时段的价格1+k y 由第1+k 和第k 时段的数量1+k x 和k x 决定,如果仍设1+k x 仍只取决于k y ,给出稳定平衡的条件,并与7.1节的结果进行比较.2.已知某商品在k 时段的数量和价格分别为k x 和k y ,其中1个时段相当于商品的一个生产周期.设该商品的需求函数和供应函数分别为)(k k x f y =和)2(11-++=k k k y y g x .试建立关于商品数量的差分方程模型,并讨论稳定平衡条件.3. 已知某商品在k 时段的数量和价格分别为k x 和k y ,其中1个时段相当于商品的一个生产周期.设该商品的需求函数和供应函数分别为)2(11kk k x x f y +=++和)(1k k y g x =+.试建立关于商品数量的差分方程模型,并讨论稳定平衡条件.《数学模型》作业解答第七章(2008年12月4日)2. 对于7.1节蛛网模型讨论下列问题:(1)因为一个时段上市的商品不能立即售完,其数量也会影响到下一时段的价格,所以第1+k 时段的价格1+k y 由第1+k 和第k 时段的数量1+k x 和k x 决定,如果仍设1+k x 仍只取决于k y ,给出稳定平衡的条件,并与7.1节的结果进行比较.(2)若除了1+k y 由1+k x 和k x 决定之外,1+k x 也由前两个时段的价格k y 和1-k y 确定.试分析稳定平衡的条件是否还会放宽.解:(1)由题设条件可得需求函数、供应函数分别为:⎪⎩⎪⎨⎧=+=+++)()2(111k k k k k y h x x x f y 在),(000y x P 点附近用直线来近似曲线h f ,,得到⎪⎩⎪⎨⎧>-=->-+-=-+++)2( 0, )()1( 0),2(0010101 ββααy y x x x x x y y k k k k k 由(2)得 )3( )(0102 y y x x k k -=-++β(1)代入(3)得 )2(0102x x x x x kk k -+-=-++αβ 0012222 x x x x x k k k αβαβαβ+=++∴++对应齐次方程的特征方程为 02 2=++αβαβλλ特征根为48)(22,1αβαβαβλ-±-=当8≥αβ时,则有特征根在单位圆外,设8<αβ,则248)()4(2222,1αβαβαβαβλ=+-+= 2 12,1<⇔<∴αβλ即平衡稳定的条件为2 <αβ与207P 的结果一致.(2)此时需求函数、供应函数在),(000y x P 处附近的直线近似表达式分别为:⎪⎩⎪⎨⎧>-+=->-+-=--+++)5( 0 , )2()4( 0),2(01010101 ββααy y y x x x x x y y k k k k k k 由(5)得,)( ) y y y β(y )x (x k k k 62010203 -+-=-+++ 将(4)代入(6),得 ⎥⎦⎤⎢⎣⎡-+--+-=-++++)2()2()(20101203x x x x x x x x k k k k k ααβ 001234424 x x x x x x k k k k αβαβαβαβ+=+++∴+++对应齐次方程的特征方程为(7) 024 23=+++αβαβλαβλλ 代数方程(7)无正实根,且42 ,αβαβ---, αβ不是(7)的根.设(7)的三个非零根分别为321,,λλλ,则⎪⎪⎪⎩⎪⎪⎪⎨⎧-==++-=++424321133221321αβλλλαβλλλλλλαβλλλ 对(7)作变换:,12αβμλ-=则,03=++q p μμ其中 )6128(41 ),122(412233322αββαβαβααβ+-=-=q p 用卡丹公式:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧+--+++-=+--+++-=+--+++-=33233223332233223323321)3()2(2)3()2(2)3()2(2)3()2(2)3()2(2)3()2(2p q q w p q q w p q q w p q q w pq q p q q μμμ 其中,231i w +-=求出321,,μμμ,从而得到321,,λλλ,于是得到所有特征根1<λ的条件.2.已知某商品在k 时段的数量和价格分别为k x 和k y ,其中1个时段相当于商品的一个生产周期.设该商品的需求函数和供应函数分别为)(k k x f y =和)2(11-++=k k k y y g x .试建立关于商品数量的差分方程模型,并讨论稳定平衡条件.解:已知商品的需求函数和供应函数分别为)(k k x f y =和)2(11-++=k k k y y g x . 设曲线f 和g 相交于点),(000y x P ,在点0P 附近可以用直线来近似表示曲线f 和g :0,)(00 ααx x y y k k --=- ----------------------(1)0,)2(0101 ββy y y x x k k k -+=--+ --------------------(2) 从上述两式中消去k y 可得,2,1,)1(22012=+=++++k x x x x k k k αβαβαβ, -----------(3) 上述(3)式是我们所建立的差分方程模型,且为二阶常系数线性非齐次差分方程. 为了寻求0P 点稳定平衡条件,我们考虑(3)对应的齐次差分方程的特征方程:022=++αβαβλλ容易算出其特征根为48)(22,1αβαβαβλ-±-=---------------(4) 当αβ 8时,显然有448)(22αβαβαβαβλ----= -----------(5) 从而2λ 2,2λ在单位圆外.下面设8 αβ,由(5)式可以算出 22,1αβλ=要使特征根均在单位圆内,即2,1λ1 ,必须 2 αβ.故0P 点稳定平衡条件为 2 αβ.3. 已知某商品在k 时段的数量和价格分别为k x 和k y ,其中1个时段相当于商品的一个生产周期.设该商品的需求函数和供应函数分别为)2(11kk k x x f y +=++和)(1k k y g x =+.试建立关于商品数量的差分方程模型,并讨论稳定平衡条件. 解:已知商品的需求函数和供应函数分别为)2(11kk k x x f y +=++和)(1k k y g x =+. 设曲线f 和g 相交于点),(000y x P ,在点0P 附近可以用直线来近似表示曲线f 和g :0,)2(0101 ααx x x y y kk k -+-=-++ --------------------(1) 0,)(001 ββy y x x k k -=-+ --- ----------------(2)由(2)得 )(0102y y x x k k -=-++β --------------------(3) (1)代入(3),可得)2(0102x x x x x kk k -+-=-++αβ∴ ,2,1,2220012=+=++++k x x x x x k k k αβαβαβ, --------------(4) 上述(4)式是我们所建立的差分方程模型,且为二阶常系数线性非齐次差分方程. 为了寻求0P 点稳定平衡条件,我们考虑(4)对应的齐次差分方程的特征方程:022=++αβαβλλ容易算出其特征根为48)(22,1αβαβαβλ-±-=---------------(4) 当αβ≥8时,显然有448)(22αβαβαβαβλ-≤---= -----------(5) 从而2λ 2,2λ在单位圆外.下面设8 αβ,由(5)式可以算出 22,1αβλ=要使特征根均在单位圆内,即2,1λ1 ,必须 2 αβ.故0P 点稳定平衡条件为 2 αβ.《数学模型》作业解答第八章(2008年12月9日)1. 证明8.1节层次分析模型中定义的n 阶一致阵A 有下列性质: (1) A 的秩为1,唯一非零特征根为n ; (2) A 的任一列向量都是对应于n 的特征向量. 证明: (1)由一致阵的定义知:A 满足ik jk ij a a a =⋅,n k j i ,,2,1,, =于是对于任意两列j i ,,有ij jkika a a =,()n k ,,2,1 =.即i 列与j 列对应分量成比例. 从而对A 作初等行变换可得:∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−−→−00000011211 n b b b A 初等行变换B 这里0≠B .()1=∴B 秩,从而秩()1=A再根据初等行变换与初等矩阵的关系知:存在一个可逆阵P ,使B PA =,于是∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡==--0000001121111 n c c c BP PAP C 易知C 的特征根为0,,0,11 c (只有一个非零特征根).又A ~C ,A ∴与C 有相同的特征根,从而A 的非零特征根为11c ,又 对于任意矩阵有()n a a a A Tr nn n =+++=+++==+++111221121 λλλ.故A 的唯一非零特征根为n .(2)对于A 的任一列向量()Tnk k k a a a ,,,21 ,()n k ,,2,1 =有()()T nk k k nk k k n j nkn j k n j k n j jk nj n j jk j n j jk j Tnk k k a a a n na na na a a a a a a a a a a a a A ,,,,,,2121112111121121 =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=∑∑∑∑∑∑======A ∴的任一列向量()Tnk k k a a a ,,,21 都是对应于n 的特征向量.7. 右下图是5位网球选手循环赛的结果,作为竞赛图,它是双向连通的吗?找出几条完全路径,用适当方法排出5位选手的名次.解:这个5阶竞赛图是一个5阶有向Hamilton 图.其一个有向Hamilton 圈为332541→→→→→.所以此竞赛图是双向连通的.32154→→→→13542→→→→42135→→→→→→→41325→等都是完全路径.此竞赛图的邻接矩阵为⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=0011110100000010110001010A令()Te 1,1,1,1,1=,各级得分向量为()()TAe S 3,2,1,2,21==, ()()()TAS S 5,4,2,3,412==,()()()TAS S 9,7,4,6,723== , ()()()TAS S 17,13,7,11,1334==由此得名次为5,1(4),2,3 (选手1和4名次相同).注:给5位网球选手排名次也可由计算A 的最大特征根λ和对应特征向量S 得到:8393.1=λ,()TS 2769.0,2137.0,1162.0,1794.0,2137.0=数学模型作业(12月16日)解答1.基于省时、收入、岸间商业、当地商业、建筑就业等五项因素,拟用层次分析法在建桥梁、修隧道、设渡轮这三个方案中选一个,画出目标为“越海方案的最优经济效益”的层次结构图.解:目标层准则层方案层2.简述层次分析法的基本步骤. 问对于一个即将毕业的大学生选择工作岗位的决策问题要分成哪3个层次?具体内容分别是什么?答:层次分析法的基本步骤为:(1).建立层次结构模型;(2).构造成对比较阵;(3).计算权向量并做一致性检验;(4).计算组合权向量并做组合一致性检验. 对于一个即将毕业的大学生选择工作岗位的决策问题,用层次分析法一般可分解为目标层、准则层和方案层这3个层次. 目标层是选择工作岗位,方案层是工作岗位1、工作岗位2、工作岗位3等,准则层一般为贡献、收入、发展、声誉、关系、位置等.3.用层次分析法时,一般可将决策问题分解成哪3个层次?试给出一致性指标的定义以及n 阶正负反阵A 为一致阵的充要条件.答:用层次分析法时,一般可将决策问题分解为目标层、准则层和方案层这3个层次; 一致性指标的定义为:1--=n nCI λ.n 阶正互反阵A 是一致阵的充要条件为:A 的最大特征根λ=n .第九章(2008年12月18日)1.在1.9节传送带效率模型中,设工人数n 固定不变.若想提高传送带效率D,一种简单的方法是增加一个周期内通过工作台的钩子数m ,比如增加一倍,其它条件不变.另一种方法是在原来放置一只钩子的地方放置两只钩子,其它条件不变,于是每个工人在任何时刻可以同时触到两只钩子,只要其中一只是空的,他就可以挂上产品,这种办法用的钩子数量与第一种办法一样.试推导这种情况下传送带效率的公式,从数量关系上说明这种办法比第一种办法好.解:两种情况的钩子数均为m 2.第一种办法是m 2个位置,单钩放置m 2个钩子;第二种办法是m 个位置,成对放置m 2个钩子.① 由1.9节的传送带效率公式,第一种办法的效率公式为⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛--=nm n m D 21112 当mn2较小,1 n 时,有()m n m n n m n m D 41181211122--=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-+--≈E D -=1 , mnE 4≈② 下面推导第二种办法的传送带效率公式:对于m 个位置,每个位置放置的两只钩子称为一个钩对,考虑一个周期内通过的m 个钩对.任一只钩对被一名工人接触到的概率是m1; 任一只钩对不被一名工人接触到的概率是m11-;记mq m p 11,1-==.由工人生产的独立性及事件的互不相容性.得,任一钩对为空的概率为nq ,其空钩的数为m 2;任一钩对上只挂上1件产品的概率为1-n npq ,其空钩数为m .所以一个周期内通过的m 2个钩子中,空钩的平均数为 ()1122--+=⋅+⋅n n n nnpq q m npqm q m于是带走产品的平均数是 ()122-+-n n npqq m m , 未带走产品的平均数是 ()()122-+--n n npq q m m n )∴此时传送带效率公式为()⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛--⎪⎭⎫ ⎝⎛--=+-=--1111112222'n n n n m m n m n m n npq q m m D ③ 近似效率公式:由于 ()()()321621121111m n n n m n n m n m n----+-≈⎪⎭⎫ ⎝⎛- ()()2112211111mn n m n m n --+--≈⎪⎭⎫ ⎝⎛-- ∴ ()()26211'm n n D ---≈当1 n 时,并令'1'D E -=,则 226'mn E ≈ ④ 两种办法的比较:由上知:mnE 4≈,226'm n E ≈∴ m n E E 32/'=,当n m 时,132 mn, ∴ E E '. 所以第二种办法比第一种办法好.《数学模型》作业解答第九章(2008年12月23日)一报童每天从邮局订购一种报纸,沿街叫卖.已知每100份报纸报童全部卖出可获利7元.如果当天卖不掉,第二天削价可以全部卖出,但报童每100份报纸要赔4元.报童每天售出的报纸数r 是一随机变量,其概率分布如下表:试问报童每天订购多少份报纸最佳(订购量必须是100的倍数)? 解:设每天订购n 百份纸,则收益函数为⎩⎨⎧≤--+=n r n nr r n r r f 7))(4(7)( 收益的期望值为G(n) =∑=-n r r P n r 0)()411(+∑∞+=1)(7n r r P n现分别求出 n =5,4,3,2,1,0时的收益期望值. G(0)=0;G(1)=4-×0.05+7×0.1+7×(0.25+0.35+0.15+0.1)=6.45; G(2)= (05.08⨯-25.0141.03⨯+⨯+))1.015.035.0(14++⨯+8.11=; G(3)=(05.012⨯-35.02125.0101.01⨯+⨯+⨯-))1.015.0(21+⨯+4.14= G(4)=(05.016⨯-15.02835.01725.061.05⨯+⨯+⨯+⨯-)1.028⨯+15.13=G(5)=05.020⨯-1.03515.02435.01325.021.09⨯+⨯+⨯+⨯+⨯- 25.10= 当报童每天订300份时,收益的期望值最大.数模复习资料第一章。
自动控制原理(数学模型)精选全文完整版
t 0
s
证明:由微分定理 df (t) estdt s F (s) f (0)
0 dt
lim df (t) estdt lim s F (s) f (0)
s 0 dt
s
左 df (t) limestdt 0 0 dt s
lim
s
s F(s)
f (0 )
0
f
二、非线性系统微分方程的线性化
例5 已知某装置的输入输出特性如下,求小扰动线性化方程。
y( x ) E0 cos[x(t )]
解. 在工作点(x0, y0)处展开泰勒级数
y( x)
y(x0)
y( x0 )( x
x0 )
1 2!
y( x0 )( x
x0 )2
取一次近似,且令
y(x) y(x) y(x0) E 0 sin x0 ( x x0 )
1
s(s a)( s b)
f
lim
s0
s
ss
1
as
b
1 ab
例12
Fs
s2
ω ω2
f sinωt t
lim s
s0
s2
ω ω2
0
3 用拉氏变换方法解微分方程
系统微分方程
y(t) a1 y(t) a2 y(t) 1(t)
y(0) y(0) 0
L变换
(s2
a1s
a2 )Y (s)
0
1 1
1 1 2 j
2j
s
j
s
j
2j
s2
2
s2
2
2 拉氏变换的几个重要定理
(1)线性性质 La f1(t) b f2(t) a F1(s) b F2(s)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
平均每天费用950元 • 50天生产一次,每次5000件, 贮存费4900+4800+…+100 =122500元,准备费5000元,总计127500元.
平均每天费用2550元
c2 t1x x
c3 x
其中 c1,c2,c3, t1, ,为已知参数
模型求解 求 x使 C(x)最小
dC 0 dx
x
c t 2 2c t
11
21
2c 2
3
结果解释 x c1t12 2c2t1
2c32
dB
dt
/ 是火势不继续蔓延的最少队员数
x
x 0.45
0.4 0.35
0.3 0.25
0.2 0.15
0.1 0.05
0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 a
a
1
空杯质量w2取决于材料 (纸杯、塑料杯、玻璃杯).
设w2=150g 半升啤酒杯w1=500g a=0.3 x=0.3245
一杯啤酒约剩1/3时重心最低,最不容易倾倒!
问题分析与模型假设 x
w1 ~ 啤酒 (满杯) 质量
1
w2 ~空杯侧壁质量, w3 ~空杯底面质量
啤酒杯重心s(x)由啤酒重心和空杯 重心合成.
• s2=1/2 •xs(x) 液面 • s1=x/2 0
液面高度x时啤酒质量w1x, 啤酒重心位置 s1=x/2
忽略空杯底面质量w3 空杯重心位置 s2=1/2
c1=5000, c2=1,r=100
T=10(天), Q=1000(件), C=1000(元)
思考: 为什么与前面计算的C=950元有差别?
• 用于订货供应情况: 每天需求量 r,每次订货费 c1, 每天每件贮存费 c2 , T天(周期) 订货一次,每次订货Q 件,当贮存量降到零时,Q件立即到货.
已知某产品日需求量100件,生产准备费5000元,贮存费 每日每件1元. 试安排该产品的生产计划,即多少天生产 一次(生产周期),每次产量多少,使总费用最小.
要 不只是回答问题,而且要建立生产周期、产量与 求 需求量、准备费、贮存费之间的关系.
问题分析与思考
日需求100件,准备费5000元,贮存费每日每件1元. • 每天生产一次, 每次100件,无贮存费,准备费5000元.
问题分析与模型假设
x
最简单的啤酒杯 ~ 高度为1的圆柱体.
1
假设:啤酒和杯子材料均匀.
沿中轴线建立坐标轴x,倒酒时 液面高度从x=0到x=1.
重心位置沿x轴变化,记作s(x).
•xs(x) 液面 0
w1 ~ 啤酒 (满杯) 质量
w2 ~ 空杯侧壁质量 w3 ~ 空杯底面质量
空杯重心由w2和w3 决定, 与x无关.
不平坦处满杯啤酒容易倾倒. 重心太高! 满杯时重心在哪里? 杯子中央稍下一点的位置. 空杯时重心在哪里? 与满杯时重心相同. 饮酒时重心先降低,再升高,回到中央. 倒酒时重心先升高,再降低,回到中央. 重心有一个最低点 ~ 啤酒杯容易放稳的位置.
建立数学模型——描述啤酒杯的重心变化的规律, 找出重心最低点的位置,讨论决定最低点的因素.
问题
3.2 森林救火
森林失火后,要确定派出消防队员的数量. 队员多,森林损失小,救援费用大; 队员少,森林损失大,救援费用小. 综合考虑损失费和救援费,确定队员数量.
分析 记队员人数x, 失火时刻t=0, 开始救火时刻t1,
灭火时刻t2, 时刻t森林烧毁面积B(t).
• 损失费f1(x)是x的减函数, 由烧毁面积B(t2)决定. • 救援费f2(x)是x的增函数, 由队员人数和救火时间决定.
存在恰当的x,使f1(x), f2(x)之和最小.
分析 • 关键是对B(t)作出合理的简化假设.
失火时刻t=0, 开始救火时刻t1, 灭火时刻t2, 画出时刻t森林烧毁面积B(t)的大致图形.
B
分析B(t)比较困难, 转而讨论单位时间
B(t2)
烧毁面积 dB/dt (森林烧毁的速度).
O
t1
t2
t
~ C(T ) C c1 c2rT
TT 2
模型求解 求 T 使C(T ) c1 c2rT min
T2
dC 0 dT
模型解释
T 2c1 rc2
Q rT 2c1r c2
定性分析 c1 T,Q
c2 T,Q
r T ,Q
敏感性分析 参数c1,c2, r的微小变化对T,Q的影响
模型假设
1)0tt1, dB/dt 与 t成正比,系数 (火势蔓延速度). 2)t1tt2, 降为–x (为队员的平均灭火速度).
3)f1(x)与B(t2)成正比,系数c1 (烧毁单位面积损失费) 4)每个队员的单位时间灭火费用c2, 一次性费用c3 .
假设1)的解释
火势以失火点为中心,均匀向四
T对c1的(相 对)敏感度
S(T , c1)
ΔT Δ c1
/T / c1
dT c1 dc1 T
1 2
c1增加1%, T增加0.5%
S(T,c2)=–1/2, S(T,r)=–1/2 c2或r增加1%, T减少0.5%
模型应用
• 回答原问题
T 2c1 rc2
Q rT 2c1r c2
对于每个a, s(x)
s
0.5
0.45
有一最小点. 0.4
w1 ~ 啤酒质量 w2 ~ 空杯质量
a=1
a=0.3, x=0.35左右 s最小, 即重心最低.
0.35 0.3
a=0.5 a=0.3
0.25
0.2 0
a=0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
x
1
x =0.35
啤酒杯重心模型一
a = w2/w1
微分法求解s极值问题
x
0.45
0.4
0.35
0.3
0.25
x
0.2
0.15
0.1
0.05 0 0
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 a
a1ຫໍສະໝຸດ x 由质量比a决定液面高度为x时啤酒杯重心处于最低位置.
结果分析
(a = w2/w1)
重心最低位置x由比值a决定
2c32
• 费用参数c1,c2,c3已知 • 开始救火时刻t1可估计
• ,由森林类型、队员能力等因素决定,可设置
一系列数值备查.
评注
• 模型可决定队员数量 x
• 在风力的影响较大时“森林烧毁速度dB/dt 与 t 成正比”的假设需要重新考虑.
• 队员灭火速度应该与开始救火时的火势有关.
3.3 倾倒的啤酒杯
1. 产品每天的需求量为常数 r; 2. 每次生产准备费为 c1, 每天每件产品贮存费为 c2; 3. T天(一周期)生产一次, 每次生产Q件,当贮存量降
为零时,Q件产品立即生产出来 (生产时间不计); 4. 为方便起见,时间和产量都作为连续量处理.
建模目的
r, c1, c2 已知,求T, Q 使每天总费用的平均值最小.
模 型 建 立 离散问题连续化 q
贮存量表示为时间的函数 q(t)
t=0生产Q件,q(0)=Q, q(t)以
Q r
需求速率r递减,q(T)=0.
A
=QT/2
Q rT
0
T
t
一周期贮存费为
c2
T 0
q(t)dt
c2
QT 2
一周期 总费用
C~
c1
c2
QT 2
c1
c2
rT 2 2
每天总费用平均 值(目标函数)
建立啤酒杯重心模型一
啤酒杯重心模型一 s=s(x) ~ 液面高度x的啤酒杯重心
啤酒质量w1x 啤酒重心s1 力矩平衡
空杯质量w2 空杯重心s2
s1=x/2 s2=1/2 a=w2/w1
x 1
• s2=1/2 •xs(x) 液面 • s1=x/2 0
啤酒杯重心模型一
a = w2/w1
啤酒杯重心s(x)只与质量比a有关
c2
c3
R Q Q Q~不允许缺货时的产量(或订货量)
存贮模型
• 存贮模型(EOQ公式)是研究批量生产计划的 重要理论基础, 也有实际应用.
• 建模中未考虑生产费用, 为什么?在什么条件下 可以不考虑?
• 建模中假设生产能力为无限大(生产时间不计), 如果生产能力有限(是大于需求量的常数), 应作 怎样的改动?
一周期
贮存费
c2
T1 q(t)dt
0
c2 A
一周期
缺货费
c3
T T1
q(t ) dt
c3B
一周期总费用
C
c1
c2
QT1 2
c3
r(T
T1)2 2
允许缺货的存贮模型
一周期总费用 C c 1 c QT 1 c r(T T )2
2 2 1
2
1
3
1
每天总费用 平均值
C(T ,Q) C c1 c2Q2 c3 (rT Q)2
T T 2rT
2rT
(目标函数)
求 T ,Q C(T ,Q) min
C 0, C 0 为与不允许缺货的存贮模型