《数学模型》作业答案
第二章(1)(2012年12月21日)
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 ∑==3
1
.1000i i
p
方法一(按比例分配) ,35.23
1
11==
∑=i i
p
N
p q ,33.33
1
22==
∑=i i
p
N
p q 32.43
1
33==
∑=i i
p
N
p q
分配结果为: 4 ,3 ,3321===n n n 方法二(Q 值方法)
9个席位的分配结果(可用按比例分配)为:
4 ,3 ,2321===n n n
第10个席位:计算Q 值为
,17.92043223521=?=Q ,75.92404333322=?=Q 2.9331544322
3=?=Q
3Q 最大,第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 宿舍).
i
i
n p 是每席位代表的人数,取,,2,1 =i n 从而得到的i i n p 中选较大者,可使对所有的,i i
i n p
尽量接近.
再考虑15=N 的分配方案,类似地可得名额分配结果.现将3种方法两次分配的结果列表如下:
2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本.
考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得
??
+=n
t
dn wkn r k vdt 0
)(2π
)22 2
n wk k(r n πvt +=∴ .2 22n v
k w n v rk t ππ+=∴
《数学模型》作业解答
第三章1(2008年10月14日)
1. 在3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货
批量.证明在不允许缺货模型中结果与原来的一样,而在允许缺货模型中最优订货周期和订货批量都比原来结果减少.
解:设购买单位重量货物的费用为k ,其它假设及符号约定同课本.
01 对于不允许缺货模型,每天平均费用为:
kr rT c T c T C ++=
2
)(21
2221r c T
c dT dC
+-= 令
0=dT
dC
, 解得 r c c T 21
*2= 由rT Q = , 得2
12c r
c 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
--+--=??
T
k rT Q c c rT Q
c Q C ++-=??332 令???????=??=??00Q
C
T
C
, 得到驻点:
???
?
??
?+-
+-+=-
+=
**
3
23222
2
3323213
22
33221)(22c c kr
c 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 T
i i t T
T r k c dt t g c t g c 1
02
20
22
))()(lim
ξ
又 )()(00T T r T r k -=- ∴ T k r T =
0 , ∴ 贮存费变为 k
T
T r k r c 2)(2?-=
于是不允许缺货的情况下,生产销售的总费用(单位时间内)为
k
T
r k r c T c kT T r k r c T c T C 2)(2)()(21221-+=-+=
k r k r c T
c dT dC 2)(221-+-=. 0=dT dC
令
, 得)
(221r k r c k c T -=* 易得函数处在*
T T C )(取得最小值,即最优周期为: )
(221r k r c k
c T -=
*
r
c c ,T
r k 21
2≈
>>*
时当 . 相当于不考虑生产的情况.
∞→≈*
,T
r k 时当 . 此时产量与销量相抵消,无法形成贮存量.
第三章2(2008年10月16日)
3.在3.3节森林救火模型中,如果考虑消防队员的灭火速度λ与开始救火时的火势b 有关,试假设一个合理的函数关系,重新求解模型.
解:考虑灭火速度λ与火势b 有关,可知火势b 越大,灭火速度λ将减小,我们作如下假设: 1
)(+=
b k
b λ, 分母∞→→+λ时是防止中的011
b b 而加的. 总费用函数()x
c 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 kb
c x β
β)1(2)1()1(22
322
1
+++++=
5.在考虑最优价格问题时设销售期为T ,由于商品的损耗,成本q 随时间增长,设
t q t q β+=0)(,为增长率β.又设单位时间的销售量为)(为价格p bp a x -=.今将销售
期分为T t T
T t <<<<2
20和两段,每段的价格固定,记作21,p p .求21,p p 的最优值,
使销售期内的总利润最大.如果要求销售期T 内的总售量为0Q ,再求21,p p 的最优值. 解:按分段价格,单位时间内的销售量为
??
???<<-<<-=T t T bp a T t bp a x 2,2
0,21
又 t q t q β+=0)(.于是总利润为
[][]?
?--+--=2
2
221121)()()()(),(T
T
T dt bp a t q p dt bp a t q p p p
=2
2)(022)(20222011T T
t t q t p bp a T t t q t p bp a ??????
---+??????---ββ
=)8
322)(()822)((2
0222011T 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,02
1=??=??p p 令
, 得到最优价格为: ???
???
???????++=???
???++=)43(21)4(210201T q b a b p T q b a b p ββ 在销售期T 内的总销量为
??+-
=-+-=20
2
21210)(2
)()(T
T
T p p bT
aT dt bp a dt bp a Q 于是得到如下极值问题:
)8
322)(()822)((),(m ax 2
022201121T t q T p bp a T T q T p bp a p p ββ---+---=
t s . 021)(2
Q p p bT
aT =+-
利用拉格朗日乘数法,解得:
??
???+
-=--=88
0201T
bT 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 ++=212
1
)( 得:k T T
T C 10092500
)(++=
令
0=dT
dC
, 解得:3
50
92500*==T 由实际意义知:当350*
=
T (即订货周期为3
50)时,总费用将最小. 又k T C 10035095025003)(*
+?+?==300+100k
k T C 10030930
2500
)(0+?+==353.33+100k
)(0T C -)(*T C =(353.33+100k )-(300+100k )32
=53.33.
故应改变订货策略.改变后的订货策略(周期)为T *=3
50
,能节约费用约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+30y
s.t. ??
?
??∈≥≤+≤+Z y x y x y x y x ,,0,7045202
这是一个整线性规划问题,现用图解法进行求解
可行域为:由直线1l :x+2y=20, 2l :5x+4y =70
2l
92500
2+-=T
dT dC
以及x=0,y=0组成的凸四边形区域.
直线l :20x+30y=c 在可行域内 平行移动.
易知:当l 过1l 与2l 的交点时, x S 取最大值.
由???=+=+7045202y x y x 解得?
??==510
y 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,135224452
12121
这是一个整线性规划问题. 用图解法求解. 可行域为:由直线
2445:211=+x x l
1352:212=+x x l 及0,021==x x 组成直线 c x x l =+211020:在此凸四边形区域内
平行移动.
2l
l
1x
1l
2x
易知:当l 过l 1与l 2的交点时,z 取最大值
由???=+=+135224
4521
21x x x x 解得 ??
?==1
42
1
x x
90110420max =?+?=z .
3.某微波炉生产企业计划在下季度生产甲、乙两种型号的微波炉.已知每台甲型、乙型微波炉的销售利润分别为3和2个单位.而生产一台甲型、乙型微波炉所耗原料分别为2和3个单位,所需工时分别为4和2个单位.若允许使用原料为100个单位,工时为120个单位,且甲型、乙型微波炉产量分别不低于6台和12台.试建立一个数学模型,确定生产甲型、乙型微波炉的台数,使获利润最大.并求出最大利润.
解:设安排生产甲型微波炉x 件,乙型微波炉y 件,相应的利润为S. 则此问题的数学模型为:
max S=3x +2y
s.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 取最大值.
由??
?=+=+120
24100
32y x y x 解得
??
?==20
20
y x .
m ax S =320220?+?=100.
《数学模型》作业解答
第五章1(2008年11月12日)
1.对于5.1节传染病的SIR 模型,证明: (1)若处最大先增加,在则σ
σ
1
)(,1
0=
s t i s ,然后减少并趋于零;)(t s 单调减少
至.∞s
(2).)()(,1
0∞s t s t i s 单调减少至单调减少并趋于零,则若σ
解:传染病的SIR 模型(14)可写成
?????-=-=i s dt
ds s i dt di
λσμ)
1(
.)(lim 0.(t) .)( .0,t 存在而单调减少知由∞∞→=∴≥-=s t s s t s dt
ds
i s dt ds λ
.)(∞s t s 单调减少至故
(1).s s(t) .s(t) .1
00≤∴单调减少由若σ
s
;)(,0 .01,1
0单调增加时当
t i dt
di
s s s ∴
-σσ
.)(,0
.01,1单调减少时当t i dt
di
s s ∴-σσ .0)(lim .0)18(t ==∞
→∞t i i 即式知又由书上
.)( .0,
1
m i t i dt
di
s 达到最大值时当∴==
σ
(2)().0 0.1-s ,1,10 dt
di
t s s σσσ从而则若
()().0.0lim ==∴∞∞
→i t i t i t 即单调减少且
4.在5.3节正规战争模型(3)中,设乙方与甲方战斗有效系数之比为.4=b
a
初始兵力00y x 与相同.
(1) 问乙方取胜时的剩余兵力是多少,乙方取胜的时间如何确定.
(2) 若甲方在战斗开始后有后备部队以不变的速率r 增援,重新建立模型,讨论如何判断双方的胜负.
解:用()()t y t x ,表示甲、乙交战双方时刻t 的士兵人数,则正规战争模型可近似表示为:
()()()???
????==-=-=000,01 ,y
y x x bx dt
dy
ay dt dx
现求(1)的解: (1)的系数矩阵为??
?
???--=00b a A
ab ab b a
A E ±=∴=-==
-1,22 .0λλλ
λλ ???
?
??????
??-1212,21,对应的特征向量分别为λλ ()()()t
ab t ab e
C e C t y t x -
???
? ??+???
?
??-=???? ??∴1212121的通解为.
再由初始条件,得
()()2 220000 t
ab t
ab e y x e
y x t x -??
? ??++??
? ??-=
又由().1ay
bx dx dy =可得
其解为 ()3 ,2
02022 bx ay k k bx ay -==-而
(1) ()().2
3
1000202011y a b y a bx ay a
k t y t x =-=-==
=时,当 即乙方取胜时的剩余兵力数为
.2
3
0y 又令().0222,01
1
00001=-??
?
??++??
?
??-=t ab t ab e y x e y x t x )得由(
注意到0
00020022,1
x y y x e
y x t ab -+=
=得. .43
ln ,312
1
b
t e
t ab =
∴=∴ (2) 若甲方在战斗开始后有后备部队以不变的速率r 增援.则
()()???
????==-=+-=000,)0(4 y
y x x bx dt
dy
r ay dt dx
().,4rdy aydy bxdx bx
r
ay dy dx -=-+-=即得
由 相轨线为,222k bx ry ay =-- .222
2
20.020k a r bx a r y a bx ry ay k =--??
? ??---=或 此相轨线比书图11中的轨线上移了
.a r 乙方取胜的条件为.,022202
0a 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 V
D
t C V D C t f -==
=解得 (2)恒速静脉滴注(持续时间为τ): 设滴注速率为()(),00,000==C k t f k ,则解得
()()()()?????-≤≤-=----τττ t e e Vk
k t e Vk
k t C t k kt kt
,10 ,10
(3) 口服或肌肉注射: ()(),解得)式节(见134.5010010t
k e
D k t f -=
()()()
???
????=≠--=---0101
01001 ,,01k k te V kD k k e e k k V D k t C kt t k kt
3种情况下的血药浓度曲线如下:
第五章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 e
b
a v aw Q v bl a v
l β ()10/10==l M w 其中,
()()97628571.050
20
02.008.02
1
2
===?--
--e
e Q Q v
l b β
(2) 对于一支不带过滤嘴的香烟,全部吸完的毒物量为???
?
?
?
-=-v
bl a e b a v aw Q '
103‘ 只吸到1l 处就扔掉的情况下的毒物量为???
? ?
?-=--v
bl a v bl
e 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 v
bl 44.235,84.2954
3≈≈ Q
Q
4.在5.3节正规战争模型(3)中,设乙方与甲方战斗有效系数之比为.4=b
a
初始兵力00y x 与相同.
(1) 问乙方取胜时的剩余兵力是多少,乙方取胜的时间如何确定.
(2) 若甲方在战斗开始后有后备部队以不变的速率r 增援,重新建立模型,讨论如何判断双方的胜负.
解:用()()t y t x ,表示甲、乙交战双方时刻t 的士兵人数,则正规战争模型可近似表示为:
()()()???
????==-=-=000,01 ,y
y x x bx dt
dy
ay dt dx
现求(1)的解: (1)的系数矩阵为?
?
?
?
??--=00b a A ab ab b a
A E ±=∴=-==
-1,22 .0λλλ
λλ
???
?
??????
??-1212,21,对应的特征向量分别为λλ ()()()t
ab t ab e
C e C t y t x -
???
? ??+???
? ??-=???? ??∴1212121的通解为.
再由初始条件,得
()()2 220000 t
ab t
ab e y x e
y x t x -??
? ??++??
? ??-=
又由().1ay
bx dx dy =可得
其解为 ()3 ,2
02022 bx ay k k bx ay -==-而
(1) ()().2
3
1000202011y a b y a bx ay a
k t y t x =-=-==
=时,当 即乙方取胜时的剩余兵力数为
.2
3
0y 又令().0222,01
1
00001=-??
?
??++??
?
??-=t ab t ab e y x e y x t x )得由(
注意到0
00020022,1
x y y x e
y x t ab -+=
=得. .43
ln ,312
1
b
t e
t ab =
∴=∴ (2) 若甲方在战斗开始后有后备部队以不变的速率r 增援.则
()()???
????==-=+-=000,)0(4 y
y x x bx dt
dy
r ay dt dx
().,4rdy aydy bxdx bx
r ay dy dx -=-+-=即得
由 相轨线为,222k bx ry ay =-- .222
2
20.020k a r bx a r y a bx ry ay k =--??? ??---=或 此相轨线比书图11中的轨线上移了
.a r 乙方取胜的条件为.,022202
0a 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 N
x
rx dt t dx --=)1()( 记h N
x
rx x F --
=)1()( (1).讨论渔场鱼量的平衡点及其稳定性: 由()0=x F ,得0)1(=--
h N
x
rx . 即
()102
=+-h rx x N
r )4(42N
h
r r N rh r -=-
=? , (1)的解为:2
412,1N rN
h
N x -
±=
①当4/rN h >,0,(1)无实根,此时无平衡点;
②当4/rN h =,0=?,(1)有两个相等的实根,平衡点为2
0N x =
. N
rx
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 不稳定;
③当4/rN h <,0>?时,得到两个平衡点:
2411N rN
h
N x --=
, 2
412N rN
h N x -
+=
易知:21N x <
, 2
2N x > ,0)(1'
>x F ,0)(2' (2)最大持续产量的数学模型为 ? ? ? =0)(..max x F t s h 即 )1(max N x rx h -=, 易得 2*0N x = 此时 4rN h =, 但2 * 0N x =这个平衡点不稳定.这是与6.1节的产量模型不同之处. 要获得最大持续产量,应使渔场鱼量2 N x > ,且尽量接近2N ,但不能等于2N . 2.与Logistic 模型不同的另一种描述种群增长规律的是Gompertz 模型:()x N rx t x ln ' =.其中r 和N 的意义与Logistic 模型相同. 设渔场鱼量的自然增长服从这个模型,且单位时间捕捞量为Ex h =.讨论渔场鱼量的平衡点及其稳定性,求最大持续产量m h 及获得最大产量的捕捞强度m E 和渔场鱼量水平* 0x . 解:()t x 变化规律的数学模型为 ()Ex x N rx dt t dx -=ln 记 Ex x N rx x F -=ln )( ① 令()0=x F ,得0ln =-Ex x N rx ∴r E Ne x -=0,01=x . ∴平衡点为1,0x x . 又 ()E r x N r 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 由前面的结果可得 r E ENe h -= r E r E e r EN Ne dE dh ---=,令.0=dE dh 得最大产量的捕捞强度r E m =.从而得到最大持续产量e rN h m /=,此时渔场鱼量水平 e N x = * 0. 3.设某渔场鱼量)(t x (时刻t 渔场中鱼的数量)的自然增长规律为:)1()(N x rx dt t dx -= 其中r 为固有增长率,`N 为环境容许的最大鱼量. 而单位时间捕捞量为常数h . 10 .求渔场鱼量的平衡点,并讨论其稳定性; 20 .试确定捕捞强度m E ,使渔场单位时间内具有最大持续产量m Q ,求此时渔场鱼量水平* 0x . 解:10 .)(t x 变化规律的数学模型为 h N x rx dt t dx --=)1()( 记h N x rx x f --=)1()(,令 0)1(=--h N x rx ,即 02 =+-h rx x N r ----(1))4(42N h r r N rh r -=- =? , (1)的解为:2 412,1N rN h N x - ±= ① 当0 ?时,(1)无实根,此时无平衡点; Ex ()x f ② 当0=?时,(1)有两个相等的实根,平衡点为2 0N x = . N rx 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 ?时,得到两个平衡点: 2411rN h N N x --= , 2 412rN h N N x - += 易知 21N x , 2 2N 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 ,但不能等于2 N . 《数学模型》第七章作业 (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 =和 数学建模作业 姓名:李成靖 学号:1408030311 班级:计科1403班 日期:2015.12。30 1.某班准备从5名游泳队员中选4人组成接力队,参加学校的4×100m混合泳接力比赛,5名队员4种泳姿的百米平均成绩如下表所示,问应如何选拔队员组成接力队? 如果最近队员丁的蛙泳成绩有较大的退步,只有1′15"2;而队员戊经过艰苦训练自由泳成绩有所进步,达到57”5,组成接力队的方案是否应该调整? 名队员4种泳姿的百米平均成绩 ij 若参选择队员i 加泳姿j 的比赛,记x i j=1, 否则记xi j=0 目标函数: 即m in=66.8*x11+75.6*x12+87*x13+58.6*x14+57。2*x21+66*x22+66.4*x 23+53*x24+78*x31+67.8*x32+84。6*x33+59.4*x34+70*x 41+74。2*x42+69.6*x 43+57。2*x44+67。4*x51+71*x52+83。8*x53+62.4*x54; 约束条件: x 11+x12+x13+x14〈=1; x 21+x22+x23+x 24〈=1; x 31+x32+x33+x34<=1; x 41+x42+x 43+x44〈=1; x 51+x52+x53+x54<=1; x11+x 21+x31+x41+x51=1; x 12+x22+x32+x42+x52=1; x13+x 23+x33+x43+x53=1; x14+x24+x 34+x44+x54=1; 甲 乙 丙 丁 戊 蝶泳 1′06"8 57”2 1′18” 1′10” 1′07"4 仰泳 1′15"6 1′06" 1′07”8 1′14"2 1′11" 蛙泳 1′27” 1′06"4 1′24"6 1′09"6 1′23"8 自由泳 58"6 53” 59”4 57”2 1′02”4 ∑∑=== 415 1j i ij ij x c Z Min 《数学模型》作业解答 第七章( 2008 年 12 月 4 日) 1.对于节蛛网模型讨论下列问题: ( 1)因为一个时段上市的商品不能立即售完,其数量也会影响到下一时段的价格,所以第 k 1时段的价格y k 1由第k 1 和第 k 时段的数量x k 1和x k决定,如果仍设x k 1仍只取 决于 y k ,给出稳定平衡的条件,并与节的结果进行比较 . ( 2)若除了 y k 1 由 x k 1 和 x k 决定之外, x k 1 也由前两个时段的价格 析稳定平衡的条件是否还会放宽 . 解:( 1)由题设条件可得需求函数、供应函数分别为: y k 1 f x k 1 x k ) ( 2 x k 1 h( y k ) 在 P 0 (x 0 , y 0 ) 点附近用直线来近似曲线 f , h ,得到 y k 1 y 0 ( x k 1 x k x 0 ), 2 x k 1 x 0 ( y k y 0 ) , 由( 2)得 x k 2 x 0 ( y k 1 y 0 ) ( 1)代入( 3)得 x k 2 x 0 ( x k 1x k x 0 ) 2 2x k 2 x k 1 x k 2x 0 2 x 0 对应齐次方程的特征方程为 2 2 ( ) 2 8 特征根为 1, 2 4 y k 和 y k 1 确定 . 试分 (1) ( 2) (3) 当 8 时,则有特征根在单位圆外,设 8 ,则 1,2 ( ) 2 ( ) 2 8 42 2 4 1,2 1 2 即平衡稳定的条件为 2与 P 207 的结果一致 . ( 2)此时需求函数、供应函数在 P 0 (x 0 , y 0 ) 处附近的直线近似表达式分别为: y k 1 y 0 ( x k 1 x k x 0 ), ( 4) 2 x k 1 x 0 ( y k y k 1 y 0 ) , ( 5) 2 由( 5)得, (x x 0 ) β(y y y k 1 y 0 ) ( 6 ) 2 k 3 k 2 将( 4)代入( 6),得 2( x k 3 x 0 ) ( x k 2 x k 1 x 0 ) ( x k 1 x k x 0 ) 2 2 4 x k 3x k 2 2 x k 1 x k 4 x 0 4 x 0 对应齐次方程的特征方程为 4 3 2 2 0 (7) 代数方程( 7 )无正实根,且 αβ , , 2 4 不是( 7)的根 . 设( 7)的三个非零根分 别为 1, 2, 3,则 1 2 3 4 1 2 2 3 3 1 2 1 2 3 4 对( 7)作变换: , 则 12 3 q 0, p 其中 p 1 (2 2 2 ), q 1(833 2 2 ) 4 12 4 123 6 上机练习题一 班级: 姓名: 学号: 1.建立起始值=3,增量值=5.5,终止值=44的一维数组x 答案: x=(3:5.5:44) 2.写出计算 Sin(30o )的程序语句. 答案: sin(pi*30/180) 或 sin(pi/6) 3.矩阵??????????=187624323A ,矩阵???? ??????=333222111B ;分别求出B A ?及A 与B 中对应元素之间的乘积. 答案:A = [3,2,3; 4,2,6; 7,8,1] B = [1,1,1; 2,2,2; 3,3,3] A*B ;A.*B 4计算行列式的值1 876243 23=A 。答案:det(A) 5对矩阵 ???? ??????=187624323A 进行下述操作。 (1)求秩。答案:rank(A) (2)求转置。答案:A' (3) 对矩阵求逆,求伪逆。答案:inv(A) ,pinv(A) (4) 左右反转,上下反转。答案:fliplr(A),flipud(A) (5) 求矩阵的特征值. 答案:[u,v]=eig(A) (6) 取出上三角和下三角. 答案:triu(A) tril(A) (7)以A 为分块作一个3行2列的分块矩阵。答案:repmat(a) 6 计算矩阵??????????897473535与???? ??????638976242之和。 >> a=[5 3 5;3 7 4;7 9 8]; >> b=[2 4 2;6 7 9;8 3 6]; >> a+b 7 计算??????=572396a 与?? ????=864142b 的数组乘积。 >> a=[6 9 3;2 7 5]; >> b=[2 4 1;4 6 8]; 数学建模竞赛C题解答 ————————————————————————————————作者:————————————————————————————————日期: 2010高教社杯全国大学生数学建模竞赛C 题解答 问题1:如图1,设P 的坐标为 (x , y ), (x ≥ 0,y ≥ 0),共用管道的费用为非共用管道的k 倍,模型可归结为 2222)()()(),(min y b x l y a x ky y x f -+-+-++= 只需考虑21<≤k 的情形(不妨假设b a ≤)。对上述二元费用函数求偏导,令 ()()()()()()()()??? ? ??? =-+----+--==-+----+=0 ,0,22222222 y b x l y b y a x y a k y x f y b x l x l y a x x y x f y x (*) 结合图1,将(*)式改写为 ?? ?=+=-k βαβαsin sin 0 cos cos ,易知: 2 4cos cos ,2 sin sin 2 k k -= ===βαβα 所以 2 4tan tan k k -= =βα,故经过AP 和BP 的直线方程分别为: x k k a y 2 4-- =- ① ()l x k k b y --= -24 ② 联立①、②解方程组得交点()()?? ? ???--+= ??? ?????--- =2 2 421,421k kl b a y a b k k l x 因为 x ≥ 0,y ≥ 0,所以 l 应满足: ()a b k k l --≥ 2 4 且()a b k k l +-≤2 4 (a )当 )(42 a b k k l --≤ 时,此时交点在y 轴上,将0=x 代入①式,可得),0(a P =,即交点P 与A 点重合(如图2)。 ka l a b f ++-=22min )( (b) 当)(4)(42 2 a b k k l a b k k +-< <--时,交点在梯形内(如图1) 。??? ? ? ?--+---=)4(21),(24222k kl b a a b k k l P , 因为 2 42cos cos cos k l l x l x BP AP -==-+= +α βα,所以模型简化为: 2 42),(min k l ky y x f -+ =, () l k k b a f 2min 4)(2 1 -++= (c) 当)(42 a b k k l +-≥ 时,此时交点在x 轴上,即无共用管线的情形(如图3) 。 数学建模作业 1、在甲乙双方的一场战争中,部分甲方部队被乙方部队包围长达4个月,乙方封锁了所有 水陆交通通道,因此被包围的甲方只能依靠空中交通维持补给,运送4个月的供给依此分别 需要2次、3次、3次、4次飞行,每次飞行编队由50架飞机组成,每架飞机都需要3名飞 行员,每架飞机每月只能飞行一次,每名飞行员每月也只能飞行一次,每次执行完运输飞行 任务后的返回途中有20%的飞机被乙方部队击落,导致机上的飞行员也牺牲或失踪。在第 一个月开始时,甲方拥有110架飞机和330名熟练的飞行员,每个月开始时,甲方可以招聘 新飞行员和购买新飞机,新飞机必须经过一个月的检查磨合后才可以投入使用,新飞行员也 必须在熟练飞行员的指导下经过一个月的训练才能成为熟练飞行员而投入飞行(作为教练的 熟练飞行员本月不能参与飞行任务),每名熟练飞行员作为教练每月指导20名飞行员(包括 自己在内)进行训练,每名飞行员在完成本月的飞行任务后必须有一个月的带薪休假,然后 返回待命可再次投入飞行,已知各项费用平均单价如下表所示(单位:千元)。 第一个月第二个月第三个月第四个月新飞机价格200 195 190 185 闲置的熟练飞行员报酬7 6.9 6.8 6.7 10 9.9 9.8 9.7 教练及飞行员报酬和训练 费用 执行飞行任务的飞行员报 9 8.9 9.8 9.7 酬 休假期的飞行员报酬 5 4.9 4.8 4.7 (1)为甲方安排一个总费用最小的飞行计划。 (2)如果每名熟练飞行员作为教练每月指导不超过20名飞行员(包括自己在内)进行训练, 相应的模型和安排将会发生怎样的改变? 解:(1) 设每月初购买飞机数量为d1,d2,d3,d4架,每月闲置飞机数量为 y1,y2,y3,y4架,每月教练与新飞行员总数量为a1,a2,a3,a4人,每月闲置熟练 飞行员的数量为b1,b2,b3,b4人。由于每月执行任务的飞行员和休假期的飞行员 的数量是固定的,即这部分的花费是固定的,所以在优化目标中可以不必考虑。 模型建立: 决策变量:设每月初购买飞机数量为d1,d2,d3,d4架,每月闲置飞机数量 为y1,y2,y3,y4架,每月教练与新飞行员总数量为a1,a2,a3,a4人,每月闲置熟 练飞行员的数量为b1,b2,b3,b4人。 目标函数:设总费用为z元,则由价格平均表可知: z=200d1+195d2+190d3+185d4+10a1+9.9a2+9.8a3+9.7a4+7b1+6.9b2+6.8b3+ 6.7b4 约束条件包括: (1)飞机数量限制:四个月中出去执行任务的飞机数量分别为100,150,150,200架次,每次安全返回的数量为80,120,120,160架次。 根据每个月的实际情况可得方程: 100+y1=110; 150+y2=80+y1+d1; 150+y3=120+y2+d2; 200+y4=120+y3+d3; 第一部分课后习题 1.学校共1000名学生,235人住在A宿舍,333人住在B宿舍,432人住在C宿舍。学生 们要组织一个10人的委员会,试用下列办法分配各宿舍的委员数: (1)按比例分配取整数的名额后,剩下的名额按惯例分给小数部分较大者。 (2)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种方法两次分配的结果列表比较。 (4)你能提出其他的方法吗。用你的方法分配上面的名额。 2.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗。比如洁银牙膏50g 装的每支1.50元,120g装的3.00元,二者单位重量的价格比是1.2:1。试用比例方法构造模型解释这个现象。 (1)分析商品价格C与商品重量w的关系。价格由生产成本、包装成本和其他成本等决定,这些成本中有的与重量w成正比,有的与表面积成正比,还有与w无关的因素。 (2)给出单位重量价格c与w的关系,画出它的简图,说明w越大c越小,但是随着w 的增加c减少的程度变小。解释实际意义是什么。 3.一垂钓俱乐部鼓励垂钓者将调上的鱼放生,打算按照放生的鱼的重量给予奖励,俱乐部 只准备了一把软尺用于测量,请你设计按照测量的长度估计鱼的重量的方法。假定鱼池中只有一种鲈鱼,并且得到8条鱼的如下数据(胸围指鱼身的最大周长): 先用机理分析建立模型,再用数据确定参数 4.用宽w的布条缠绕直径d的圆形管道,要求布条不重叠,问布条与管道轴线的夹角 应 多大(如图)。若知道管道长度,需用多长布条(可考虑两端的影响)。如果管道是其他形状呢。 2010高教社杯全国大学生数学建模竞赛题目 A题储油罐的变位识别与罐容表标定 通常加油站都有若干个储存燃油的地下储油罐,并且一般都有与之配套的“油位计量管理系统”,采用流量计和油位计来测量进/出油量与罐内油位高度等数据,通过预先标定的罐容表(即罐内油位高度与储油量的对应关系)进行实时计算,以得到罐内油位高度和储油量的变化情况。 许多储油罐在使用一段时间后,由于地基变形等原因,使罐体的位置会发生纵向倾斜和横向偏转等变化(以下称为变位),从而导致罐容表发生改变。按照有关规定,需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图,其主体为圆柱体,两端为球冠体。图2是其罐体纵向倾斜变位的示意图,图3是罐体横向偏转变位的截面示意图。 请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。 (1)为了掌握罐体变位后对罐容表的影响,利用如图4的小椭圆型储油罐(两端平头的椭圆柱体),分别对罐体无变位和倾斜角为α=4.10的纵向变位两种情况做了实验,实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响,并给出罐体变位后油位高度间隔为1cm的罐容表标定值。 (2)对于图1所示的实际储油罐,试建立罐体变位后标定罐容表的数学模型,即罐内储油量与油位高度及变位参数(纵向倾斜角度α和横向偏转角度β)之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据(附件2),根据你们所建立的数学模型确定变位参数,并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。 附件1:小椭圆储油罐的实验数据 附件2:实际储油罐的检测数据 地平线油位探针 数学模型课后答案 《数学模型》作业答案 第二章(1)(2012年12月21日) 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 ∑==3 1 . 1000i i p 方法一(按比例分配) , 35.23 1 11 == ∑=i i p N p q , 33.33 1 22 == ∑=i i p N p q 32 .43 1 33 == ∑=i i p N p q 分配结果为: 4 ,3 ,3321 ===n n n 方法二(Q 值方法) 9个席位的分配结果(可用按比例分 配)为: 4 ,3 ,2321===n n n 第10个席位:计算Q 值为 2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本. 考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得 ??+=n t dn wkn r k vdt 0 )(2π ) 2 2 2 n wk k(r n πvt +=∴ . 2 2 2n v k w n v rk t ππ+=∴ 《数学模型》作业解答 第三章1(2008年10月14日) 1. 在 3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货批量.证明在不允许缺货模型中结果与原来的一样,而在允许缺货模型中最优订货周期和订货批量都比原来结果减少. 《数学模型》作业答案 第二章(1)(2012年12月21日) 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 ∑==3 1 .1000i i p 方法一(按比例分配) ,35.23 1 11== ∑=i i p N p q ,33.33 1 22== ∑=i i p N p q 32.43 1 33== ∑=i i p N p q 分配结果为: 4 ,3 ,3321===n n n 方法二(Q 值方法) 9个席位的分配结果(可用按比例分配)为: 4 ,3 ,2321===n n n 第10个席位:计算Q 值为 ,17.92043223521=?=Q ,75.92404333322=?=Q 2.9331544322 3=?=Q 3Q 最大,第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 宿舍). i i n p 是每席位代表的人数,取,,2,1 =i n 从而得到的i i n p 中选较大者,可使对所有的,i i i n p 尽量接近. 再考虑15=N 的分配方案,类似地可得名额分配结果.现将3种方法两次分配的结果列表如下: 2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本. 考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得 ?? +=n t dn wkn r k vdt 0 )(2π )22 2 n wk k(r n πvt +=∴ .2 22n v k w n v rk t ππ+=∴ 《数学模型》作业解答 第三章1(2008年10月14日) 2010年数学建模竞赛 答案 输油管道的铺设设计 符号约定 m 炼油厂A 到铁路线L 的距离 n 炼油厂B 到铁路线L 的距离 b 炼油厂A 、B 间水平距离 F 输送管道的总费用 f 铺设管道的附加费用 W 铺设费用的权重系数 1k A 厂铺设非共用管线每千米的费用 2k B 厂铺设非共用管线每千米的费用 3k 共用管线每千米的费用 问题一分析与模型建立 最短路径的存在性论证 如图4.1,假设C 点为在铁路线上设计增建的车站,由费尔马问题的结论,在ABC ?中,存在费尔马点P ,使点P 与ABC ?三个顶点距离之和小于三角形二边之和,即有 PA+PB+PC ①当0120<∠ACB 时,铺设公用管道PC 的输送费用比不铺设公用管道费用低; ②当0120>∠ACB 时,不需要铺设公用管道,即公用管道PC =0。 问题一分析与模型建立 如图4.1,以炼油厂A 、B 间铁路线所在直线为x 轴,以过炼油厂A 且垂直于铁路线L 直线为y 轴,建立平面直角坐标系。设 A(0,m), B(b,n),P(r,t),并设非公用管道的费用为每千米1个单位,公用管道的费用为每千米k 个单位(下同),根据实际意义易知21<≤k 。 根据参考文献[1],点P 不可能在A 的上方,故m t ≤≤0。 易得,A 点关于过点P 平行于x 轴的直线1L 的对称点'A (0,2t-m )。 由费尔马点的应用及平面几何对称性有 111F PB PA k PC BA k PC '=?+?+?>?+? 为此,得到铺设管道的最优模型 min 1F BA k PC '=?+? 4-1 问题一模型求解 对模型分两种管道费用相同与不同两种情形研究,并根据点A 、B 的坐标不同的取值,进行A 、B 不同位置时管道铺设设计。 1公用管道与非公用管道费用不同,即k <1时模型的求解 已知A 点关于1l 对称点'A (0,2t-m ) ()F t tk = 数学建模习题指导 第一章 初等模型 讨论与思考 讨论题1 大小包装问题 在超市购物时你注意到大包装商品比小包装商品便宜这种现象吗?比如洁银牙膏50g 装的每支1.50元,120g 装的每支3.00元,二者单位重量的价格比是1.2:1,试用比例方法构造模型解释这种现象。 (1)分析商品价格C 与商品重量w 的关系。 (2)给出单位重量价格c 与w 的关系,并解释其实际意义。 提示: 决定商品价格的主要因素:生产成本、包装成本、其他成本。 单价随重量增加而减少 单价的减少随重量增加逐渐降低 思考题2 划艇比赛的成绩 赛艇是一种靠浆手划桨前进的小船,分单人艇、双人艇、四人艇、八人艇四种。各种艇虽大小不同,但形状相似。T.A.McMahon 比较了各种赛艇1964—1970年四次2000m 比赛的最好成绩(包括1964年和1968年两次奥运会和两次世界锦标赛),见下表。建立数学模型解释比赛成绩与浆手数量之间的关系。 各种艇的比赛成绩与规格 γβα++=3 2w w C w w c γβα++=-3 123 431w w c γβ--='-3 2943 4w w c γβ+=''- 第二章 线性代数模型 森林管理问题 森林中的树木每年都要有一批砍伐出售。为了使这片森林不被耗尽且每年都有所收获,每当砍伐一棵树时,应该就地补种一棵幼苗,使森林树木的总数保持不变。被出售的树木,其价值取决于树木的高度。开始时森林中的树木有着不同的高度。我们希望能找到一个方案,在维持收获的前提下,如何砍伐树木,才能使被砍伐的树木获得最大的经济价值。 思考: 试解释为什么模型中求解得到的 为每周平均销售量会略小于模型假设中给出的1。 练习: 将钢琴销售的存贮策略修改为:当周末库存量为0或1时订购,使下周初的库存 达到3架;否则,不订购。建立马氏链模型,计算稳态下失去销售机会的概率和每周的平均销售量。 2.将钢琴销售的存贮策略修改为:当周末库存量为0时订购本周销售量加2架;否则,不订购。建立马氏链模型,计算稳态下失去销售机会的概率和每周的平均销售量。 第三章 优化模型 讨论题 1)最优下料问题 用已知尺寸的矩形板材加工半径一定的圆盘。给出几种加工排列方法,比较出最优下料方案。 2)广告促销竞争问题 甲乙两公司通过广告竞争销售商品,广告费分别为 x 和 y 。设甲乙公司商品的售量在两公司总售量中所占份额是它们的广告费在总广告费中所占份额的函数 又设公司的收入与售量成正比,从收入中扣除广告费后即为公司的利润。试构造模型的图形,并讨论甲公司怎样确定广告费才能使利润最大。 (1)令 (2)写出甲公司的利润表达式 对一定的 y ,使 p (x ) 最大的 x 的最优值应满足什么关系。用图解法确定这个最优值。 练习1 三个家具商店购买办公桌:A 需要30张,B 需要50张,C 需要45张。这些办公桌由两个工厂供应:工厂1生产70张,工厂2生产80张。下表给出了工厂和商店的距离(单位公里) , 857.0=n R ) (),(y x y f y x x f ++的示意图。。画出则)()()(,t f t f t f y x x t 11=-++= 。 )(t p 第一章 课后习题6. 利用1.5节药物中毒施救模型确定对于孩子及成人服用氨茶碱能引起严重中毒和致命的最小剂量。 解:假设病人服用氨茶碱的总剂量为a ,由书中已建立的模型和假设得出肠胃中的药量为: )()0(mg M x = 由于肠胃中药物向血液系统的转移率与药量)(t x 成正比,比例系数0>λ,得到微分方程 M x x dt dx =-=)0(,λ(1) 原模型已假设0=t 时血液中药量无药物,则0)0(=y ,)(t y 的增长速度为x λ。由于治疗而减少的速度与)(t y 本身成正比,比例系数0>μ,所以得到方程: 0)0(,=-=y y x dt dy μλ(2) 方程(1)可转换为:t Me t x λ-=)( 带入方程(2)可得:)()(t t e e M t y λμμ λλ ----= 将01386=λ和1155.0=μ带入以上两方程,得: t Me t x 1386.0)(-= )(6)(13866.01155.0---=e e M t y t 针对孩子求解,得: 严重中毒时间及服用最小剂量:h t 876.7=,mg M 87.494=; 致命中毒时间及服用最小剂量:h t 876.7=,mg M 8.4694= 针对成人求解: 严重中毒时间及服用最小剂量:h t 876.7=,mg M 83.945= 致命时间及服用最小剂量:h t 876.7=,mg M 74.1987= 课后习题7. 对于1.5节的模型,如果采用的是体外血液透析的办法,求解药物中毒施救模型的血液用药量的变化并作图。 解:已知血液透析法是自身排除率的6倍,所以639.06==μu t e t x λ-=1100)(,x 为胃肠道中的药量,1386.0=λ )(6600)(t t e e t y λμ---= 1386.0,639.0,5.236)2(,1100,2,====≥-=-λλλu z e x t uz x dt dz t 解得:()2,274.112275693.01386.0≥+=--t e e t z t t 用matlab 画图: 图中绿色线条代表采用体外血液透析血液中药物浓度的变化情况。 从图中可以看出,采取血液透析时血液中药物浓度就开始下降。T=2时,血液中药物浓度最高,为236.5;当z=200时,t=2.8731,血液透析0.8731小时后就开始解毒。 第二章 1.用 2.4节实物交换模型中介绍的无差别曲线的概念,讨论以下的雇员和雇主之间的关系: 1)以雇员一天的工作时间和工资分别为横坐标和纵坐标,画出雇员无差别曲线族的示意图,解释曲线为什么是那种形状; 2)如果雇主付计时费,对不同的工资率画出计时工资线族,根据雇员的无差别曲线族和雇主的计时工资线族,讨论双方将在怎样的一条曲线上达成协议; 3)雇员和雇主已经达成了协议,如果雇主想使用雇员的工作时间增加到t 2,他有两种 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 齐次线性方程组为: ?? ? ??=--=+=-++0 30 32221414321y y y y y y y y 它的基本解为)1,1,3,1(-=y 由量纲i P 定理得 1131ρπs v P -=, 1 13ρλ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 -2 T=L -1 MT -1 ,[g ]=LM 0T -2 ,其中L ,M ,T 是基本量纲. 量纲矩阵为 A=) ()()()()()() (210101101131g v T M L μρ??????????----- 齐次线性方程组Ay=0 ,即 ??? ??==+=+0 2y -y - y -0 y y 0y y -3y -y 431 324321 的基本解为y=(-3 ,-1 ,1 ,1) 由量纲i P 定理 得 g v μρπ1 3 --=. 3 ρ μλg v =∴,其中λ是无量纲常数. 2011高教社杯全国大学生数学建模竞赛 承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。 我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。 我们参赛选择的题号是(从A/B/C/D中选择一项填写): B 我们的参赛报名号为(如果赛区设置报名号的话): 所属学校(请填写完整的全名):北京大学 参赛队员(打印并签名) :1. 姚胜献 2. 许锦敏 3. 刘迪初 指导教师或指导教师组负责人(打印并签名):刘业辉 日期: 2011 年 9 月 12日赛区评阅编号(由赛区组委会评阅前进行编号): 2011高教社杯全国大学生数学建模竞赛 编号专用页 赛区评阅编号(由赛区组委会评阅前进行编号): 全国评阅编号(由全国组委会评阅前进行编号): 交巡警服务平台的设置与调度 摘要 本文通过建立整数规划模型,解决了分配各平台管辖范围、调度警务资源以及合理设置交巡警服务平台这三个方面的问题;通过建立线性加权评价模型定量评价了某市现有交巡警服务平台设置方案的合理性,并根据各个区对服务平台需求量的不同,提出了重新分配全市警力资源的解决方案。在计算交巡警服务平台到各个路口节点的路程时,使用了图论里的floyd算法。 针对问题一的第一个子问题,首先假设交巡警服务平台对某个路口节点的覆盖度是二元的,引入决策变量,建立了0-1整数规划模型。交巡警出警应体现时间的紧迫性,所以选择平均每个突发事件的出警时间最短作为目标函数,运用基于MATLAB的模拟退火算法进行求解,给出了中心城区A的20个服务平台的管辖范围,求得平均每个案件的出警时间为1.013分钟。 针对问题一的第二个子问题,为了实现对中心城区A的13个交通要道的快速全封锁,以最短的封锁时间为目标,建立了0-1整数规划模型,利用lingo软件编程求解,给出了该区交巡警服务平台警力合理的调度方案,并求得对13个交通要道实现全封锁最短需要8.02分钟。 问题一的第三个子问题是交巡警服务平台的选址问题。考虑到建设新的服务平台需要投入更多的成本和警务资源,还需平衡各个服务平台的工作量。因此,以增加最少的服务平台数和服务平台工作量方差最小为目标,采用集合覆盖理论,建立了双目标0-1整数规划模型,用基于MATLAB的模拟退火算法求解出增加的服务平台数为4个,新增 的服务平台具体位置为A 28,A 40 ,A 48 ,A 88 ,并得到各个服务平台的工作强度方差为2.28。 针对问题二的第一个子问题,通过建立线性加权评价模型定量评价了该市现有交巡警服务平台设置方案的合理性,结果发现全市服务平台覆盖率较低且各个区的工作量不均衡,得出全市服务平台的布局存在明显的不合理的结论。并确定各区域人口密度、各区域公路总长度以及各区域平均每天总的发案率为各区域对交巡警需求的指标,然后根据各个区对服务平台需求量的不同,提出了较为合理的分配全市警力资源的解决方案。 对于问题二的第二个子问题,以围堵范围最小和调动警力最少的原则,通过分析案发后嫌疑犯可能到达的位置,给出了围堵方案。 关键词:交巡警服务平台 0-1整数规划模拟退火法 第一章 4.在1、3节“椅子能在不平的地面上放稳不”的假设条件中,将四脚的连线呈正方形改为长方形,其余不变。试构造模型并求解。 答:相邻两椅脚与地面距离之与分别定义为)()(a g a f 和。f 与g 都就是连续函数。椅子在任何位置至少有三只脚着地,所以对于任意的a ,)()(a g a f 和中至少有一个不为零。不妨设0)0(,0)0(g >=f 。当椅子旋转90°后,对角线互换,0π/2)(,0)π/2(>=g f 。这样,改变椅子的位置使四只脚同时着地。就归结为证明如下的数学命题: 已 知 a a g a f 是和)()(的连续函数,对任意 0)π/2()0(,0)()(,===?f g a g a f a 且,0)π/2(,0)0(>>g f 。证明存在0a ,使0)()(00==a g a f 证:令0)π/2(0)0(),()()(<>-=h h a g a f a h 和则, 由g f 和的连续性知h 也就是连续函数。 根据连续函数的基本性质, 必存在0a (0<0a <π/2)使0)(0=a h ,即0)()(00==a g a f 因为0)()(00=?a g a f ,所以0)()(00==a g a f 8 第二章 10.用已知尺寸的矩形板材加工半径一定的圆盘,给出几种简便有效的排列方法,使加工出尽可能多的圆盘。 第三章 5.根据最优定价模型 考虑成本随着销售量的增加而减少,则设 kx q x q -=0)( (1)k 就是产量增加一个单位时成本的降低 , 销售量x 与价格p 呈线性关系0,,>-=b a bp a x (2) 收入等于销售量乘以价格p :px x f =)( (3) 利润)()()(x q x f x r -= (4) 将(1)(2)(3)代入(4)求出 ka q kbp pa bp x r --++-=02)( 当k q b a ,,,0给定后容易求出使利润达到最大的定价*p 为 b a kb ka q p 2220*+--= 交巡警服务平台的设置与调度优化分析 摘要 本文以实现警察的刑事执法、治安管理、交通管理、服务群众四大职能为宗旨,利用有限的警务资源,根据城市的实际情况与需求合理地设置了交巡警服务平台、分配各平台的管辖范围及调度警务资源。并分别对题目的各问,作了合理的解答。 问题一: (1)、根据题目所给数据,确定各节点之间的相邻关系和距离,利用Floyd算法及matlab编程求出两点之间的最短距离,使其尽量满足能在3分钟内有交巡警平台警力到达案发结点的原则,节点去选择平台,把节点分配给离节点距离最近的平台管辖,据此,我们得到了平台的管辖区域划分。 (2)、我们对进出该区的13条交通要道实现快速全封锁的问题,我们认定在所有调度方案中,某种方案中耗时最长的的围堵时间最短即最佳方案,利用0-1变量确定平台的去向,并利用线性规划知识来求解指派问题,求得了最优的调度方案。 (3)、在确定增添平台的个数和具体位置的问题中,我们将尽量保证每个节点都有一个平台可以在三分钟内到达作为主要原则来求解。我们先找出到达每个平台的时间都超过三分钟的节点,并尝试在这些节点中选取若干个作为新的平台,求出合理的添加方案。 问题二: (1)、按照设置交巡警服务平台的原则和任务,分析现有的服务平台的设置是否合理,我们以各区覆盖率作为服务平台分布合不合理的评价标准,得到C、D、E、F区域平台设置不合理。并尝试一些新的设置方案使得设置更为合理,最后以覆盖率最低的E区为例,使用一种修改方案得到一个比原方案更合理的交巡警服务平台的设置方案。 (2)、追捕问题要求在最快的时间内抓到围堵罪犯,在罪犯和警察的行动速度一致的前提假设下,我们先设定一个具体较小的时间,编写程序检验在这个时间内是否可以成功抓捕罪犯,不行则以微小时间间隔增加时间,当第一次成功围堵时,这个时间即为最佳围堵方案。 关健字:MATLAB软件,0-1规划,最短路,Floyd算法,指派问题 一、问题重述 “有困难找警察”,是家喻户晓的一句流行语。警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。为了更有效地贯彻实施这些职能,需要在市区的一些交通要道和重要部位设置交巡警服务平台。每个交巡警服务平台的职能和警力配备基本相同。由于警务资源是有限的,如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。 试就某市设置交巡警服务平台的相关情况,建立数学模型分析研究下面的问题: 数学建模部分课后习题解答 中国地质大学 能源学院 华文静 1、在稳定的椅子问题中,如设椅子的四脚连线呈长方形,结论如何? 解: 模型假设 (1) 椅子四条腿一样长,椅脚与地面接触处视为一点,四脚的连线呈长方形 (2) 地面高度就是连续变化的,沿任何方向都不会出现间断(没有像台阶那样的情 况),即从数学角度来瞧,地面就是连续曲面。这个假设相当于给出了椅子能放稳的必要条件 (3) 椅子在任何位置至少有三只脚同时着地。为了保证这一点,要求对于椅脚的间距 与椅腿的长度而言,地面就是相对平坦的。因为在地面上椅脚间距与椅腿长度的尺寸大小相当的范围内,如果出现深沟或凸峰(即使就是连续变化的),此时三只脚就是无法同时着地的。 模型建立 在上述假设下,解决问题的关键在于选择合适的变量,把椅子四只脚同时着地表示出来。首先,引入合适的变量来表示椅子位置的挪动。生活经验告诉我们,要把椅子通过挪动放稳,通常有拖动或转动椅子两种办法,也就就是数学上所说的平移与旋转变换。然而,平移椅子后问题的条件没有发生本质变化,所以用平移的办法就是不能解决问题的。于就是可尝试将椅子就地旋转,并试图在旋转过程中找到一种椅子能放稳的情形。 注意到椅脚连线呈长方形,长方形就是中心对称图形,绕它的对称中心旋转180度后,椅子仍在原地。把长方形绕它的对称中心旋转,这可以表示椅子位置的改变。于就是,旋转角度θ这一变量就表示了椅子的位置。为此,在平面上建立直角坐标系来解决问题。 设椅脚连线为长方形ABCD,以对角线AC 所在的直线为x 轴,对称中心O 为原点,建立平面直角坐标系。椅子绕O 点沿逆时针方向旋转角度θ后,长方形ABCD 转至A1B1C1D1的位置,这样就可以用旋转角)0(πθθ≤≤表示出椅子绕点O 旋转θ后的位置。 其次,把椅脚就是否着地用数学形式表示出来。当椅脚与地面的竖直距离为零时,椅脚就着地了,而当这个距离大于零时,椅脚不着地。由于椅子在不同的位置就是θ的函数,因此,椅脚与地面的竖直距离也就是θ的函数。 由于椅子有四只脚,因而椅脚与地面的竖直距离有四个,它们都就是θ的函数,而由假设(3)可知,椅子在任何位置至少有三只脚同时着地,即这四个函数对于任意的θ,其函数值至少有三个同时为0。因此,只需引入两个距离函数即可。考虑到长方形ABCD 就是对称中心图形,绕其对称中心O 沿逆时针方向旋转180度后,长方形位置不变,但A,C 与B,D 对换了。因此,记A,B 两脚与地面竖直距离之与为)(θf ,C,D 两脚之与为)(θg ,其中[] πθ,0∈,使得)()(00θθg f =成立。 模型求解 如果0)0()0(== g f ,那么结论成立。 如果)0(与) 0(g f 不同时为零,不妨设.0)0(,0)0(=>g f 这时,将长方形ABCD 绕点O 参赛密码 (由组委会填写) 第十二届“中关村青联杯”全国研究生 数学建模竞赛 学校 参赛队号 队员姓名 参赛密码 (由组委会填写) 第十二届“中关村青联杯”全国研究生 数学建模竞赛 题目数控加工刀具运动的优化控制 摘要: 本文基于计算机数控系统的工作原理,建立了刀具运动的优化控制模型,目的在于寻求机床刀具在单个坐标轴方向上的运动合理控制,从而增强机床运行的平稳性。主要运用了S型曲线的加减速控制方法,建立了通用模型,该模型可通过已经设定的刀具加工路径,得出机床运动过程中任意一点的速度,从而验证所设定的符合加减速控制原理,得到最优的数控加工刀具的路径。在该通用模型中,机床控制的加速度和速度都是连续变化的,因此通过渐变控制使机床运动按S型曲线式平稳变化,保证了速度的光顺及加速度的连续,提高了机床运动的平稳性,运用该模型,可以帮助寻找最优刀具路径,从而实现数控刀具加工的优化。 本论文的创新点在于模型适用范围广,突破了速度范围和加速度的限制不仅适用于S型曲线七阶段的加减速,而且适用于平稳性更强的五阶段和三阶段的S型曲线加减速控制路径。 论文中主要采用了力学分析建模、直线插补法建模和最优化方法建模。在直线插补模型中,不论运行轨迹是直线还是曲线,刀具的运行都是按阶梯形路径行走,用步长乘以步数即可求得刀具的运行长度。并且每一步长的增量均为分辨率???,并且每个增量的长度均为分辨率的整数倍。根据此原理,采用直线插补,, x y z 法,建模可画出刀具沿轨迹的路径变化,在模型中输入刀具起点坐标和终点坐标即可求得刀具沿路径运行的长度。 对于问题一:根据问题二的相关提示,我们设定加工线型分别为正方形和八边形即转角分别为90°和135°,然后根据S型曲线的减加速控制方法,建立了力学分析模型,再运用牛顿第二定理和受力分析可得出速度变化特征。分别对刀具在拐角为90°和135°处进行受力分析得到结果:转角为90°时的合力F2>0.765F2(135°数学建模作业
数学模型第三版课后习题答案.doc
数学模型习题解答解读
数学建模竞赛C题解答
数学建模作业43508
数学建模习题及答案课后习题
全国数学建模大赛题目
数学模型课后答案
数学模型课后答案
最新数学建模竞赛答案汇总
数学建模习题指导
数学建模课后习题答案
(完整版)数学模型第二章习题答案
数学建模b题标准答案
数学建模课后答案
高教社杯全国大学生数学建模竞赛B题参考答案
数学建模习题答案
全国研究生数学建模比赛E题解答