姜启源编数学模型第四版
数学建模参考书大全
专业性参考书(这方面书籍很多，仅列几本供参考) ：1、数学模型，姜启源编，高等教育出版社(1987年第一版，1993年第二版，2003年第三版，2011年第四版;第一版在1992年国家教委举办的第二届全国优秀教材评选中获"全国优秀教材奖").2.数学模型与计算机模拟，江裕钊、辛培情编，电子科技大学出版社，(1989).3.数学模型选谈(走向数学从书)，华罗庚，王元著，王克译，湖南教育出版社;(1991).4.数学建模--方法与范例，寿纪麟等编，西安交通大学出版社(1993).5.数学模型，濮定国、田蔚文主编，东南大学出版社(1994).6..数学模型，朱思铭、李尚廉编，中山大学出版社，(1995)7.数学模型，陈义华编著，重庆大学出版社，(1995)8.数学模型建模分析，蔡常丰编著，科学出版社，(1995).9.数学建模竞赛教程，李尚志主编，江苏教育出版社，(1996).10.数学建模入门，徐全智、杨晋浩编，成都电子科大出版社，(1996).11.数学建模，沈继红、施久玉、高振滨、张晓威编，哈尔滨工程大学出版社，(1996).12.数学模型基础，王树禾编著，中国科学技术大学出版社，(1996).13.数学模型方法，齐欢编著，华中理工大学出版社，(1996).14.数学建模与实验，南京地区工科院校数学建模与工业数学讨论班编，河海大学出版社，(1996).15.数学模型与数学建模，刘来福、曾文艺编，北京师范大学出版杜(1997).16. 数学建模，袁震东、洪渊、林武忠、蒋鲁敏编，华东师范大学出版社。
17.数学模型，谭永基，俞文吡编，复旦大学出版社，(1997).18.数学模型实用教程，费培之、程中瑗层主编，四川大学出版社，(1998).19.数学建模优秀案例选编(工科数学基地建设丛书)，汪国强主编，华南理工大学出版社，(1998).20.经济数学模型(第二版)(工科数学基地建设丛书)，洪毅、贺德化、昌志华编著，华南理工大学出版社，(1999).21.数学模型讲义，雷功炎编，北京大学出版社(1999).22.数学建模精品案例，朱道元编著，东南大学出版社，(1999)，23.问题解决的数学模型方法，刘来福，曾文艺编著、北京师范大学出版社，(1999).24.数学建模的理论与实践，吴翔，吴孟达，成礼智编著，国防科技大学出版社，(1999).25、数学建模案例分析，白其岭主编，海洋出版社，(2000年，北京).26.数学实验(高等院校选用教材系列)，谢云荪、张志让主编，科学出版社，(2000).27.数学实验，傅鹏、龚肋、刘琼荪，何中市编，科学出版社，(2000).28.数学建模与数学实验，赵静、但琦编，高等教育出版社，(2000).国外参考书(中译本)：1、数学模型引论，E.A。
第四章 数学规划模型 数学建模(姜启源第四版)ppt课件
12小时
3公斤A1
4公斤A2
获利24元/公斤
获利16元/公斤
8小时 每天 50桶牛奶 时间480小时 至多加工100公斤A1 x1桶牛奶生产A1 x2桶牛奶生产A2
决策变量
目标函数
获利 24×3x1 获利 16×4 x2 每天获利 Max z 72x1 64x2 原料供应
x1 x2 50
基本模型
变量
目标 函数 约束 条件
x5 kg A1加工B1， x6 kg A2加工B2 利润
Max z 24x1 16x2 44x3 32x4 3x5 3x6
x1 x5 x 2 x6 加工能力 50 3 4 附加约束 4( x1 x5 ) 2( x2 x6 )
4公斤A2
获利16元/公斤
每天： 50桶牛奶
时间480小时 至多加工100公斤A1
制订生产计划，使每天获利最大 • 35元可买到1桶牛奶，买吗？若买，每天最多买多少? • 可聘用临时工人，付出的工资最多是每小时几元? • A1的获利增加到 30元/公斤，应否改变生产计划？
基本 1桶 模型 牛奶 或
线性规划模型
A1,A2每公斤的获利是与各自 产量无关的常数
每桶牛奶加工A1,A2的数量, 时 间是与各自产量无关的常数 A1,A2每公斤的获利是与相互 产量无关的常数 每桶牛奶加工A1,A2的数量,时 间是与相互产量无关的常数 加工A1,A2的牛奶桶数是实数
可 加 性
连续性
模型求解
x1 x2 50
20桶牛奶生产A1, 30桶生产A2，利润3360元。
结果解释
Global optimal solution found. Objective value: 3360.000 Total solver iterations: 2 Variable Value Reduced Cost X1 20.00000 0.000000 X2 30.00000 0.000000 Row Slack or Surplus Dual Price 1 3360.000 1.000000 原料无剩余 MILK 0.000000 48.00000 三 TIME 0.000000 2.000000 时间无剩余 种 CPCT 40.00000 0.000000 加工能力剩余40
姜启源编《数学模型》第四版_第七章_稳定性模型
设x(t)是方程的解，若从x0 某邻域的任一初值出发，
都有
lim
t
x(t)
x0
,
称x0是方程(1)的稳定平衡点.
不求x(t), 判断x0稳定性的方法——直接法
F (x0 ) 0 x0 稳定 F (x0 ) 0 x0 不稳定
第六页，共61页。
产量模型 x(t) F (x) rx(1 x ) Ex N
c2 )
p2N 2
第九页，共61页。
捕捞 过度
• 封闭式捕捞追求利润R(E)最大
• 开放式捕捞只求利润R(E) > 0
ER
r (1 2
c) pN
R(E) T (E) S(E)
pNE (1
E ) cE
令
=0
r
c Es r(1 pN )
R(E)=0时的捕捞强度Es=2ER
~ 临界强度
临界强度下的渔场鱼量
x1 (t )
r1x1 (1
x1 N1
)
x2 (t)
r2 x2 (1
x 2
N
)
2
• 两种群在一起生存时，乙对甲增长的阻滞作用
与乙的数量成正比; 甲对乙有同样的作用.
模型
x1 (t )
r1 x1 1
x1 N1
1
x2 N2
x2 (t)
r2 x2 1
2
x1 N1
x2 N2
对于消耗甲的资源而言，
乙(相对于N2)是甲(相对于 N1) 的 1 倍.
建模
h(x)=Ex, E~捕捞强度
记 F(x) f (x) h(x)
有捕捞情况下渔场 鱼量满足
x(t) F (x) rx(1 x ) Ex
姜启源数学建模资料
姜启源数学建模资料简单的优化模型3.1 3.2 3.3 3.4 存贮模型生猪的出售时机森林救火最优价格3.5 血管分支3.6 消费者均衡3.7 冰山运输<i>姜启源数学建模资料</i>静态优化模型现实世界中普遍存在着优化问题静态优化问题指最优解是数不是函数静态优化问题指最优解是数(不是函数不是函数) 建立静态优化模型的关键之一是根据建模目的确定恰当的目标函数求解静态优化模型一般用微分法<i>姜启源数学建模资料</i>问题3.1存贮模型配件厂为装配线生产若干种产品，配件厂为装配线生产若干种产品，轮换产品时因更换设备要付生产准备费，产量大于需求时要付贮存费。
备要付生产准备费，产量大于需求时要付贮存费。
该厂生产能力非常大，即所需数量可在很短时间内产出。
生产能力非常大，即所需数量可在很短时间内产出。
已知某产品日需求量100件，生产准备费5000元，贮存费件生产准备费已知某产品日需求量元每日每件1元试安排该产品的生产计划，每日每件元。
试安排该产品的生产计划，即多少天生产一次(生产周期),每次产量多少，使总费用最小。
),每次产量多少一次(生产周期),每次产量多少，使总费用最小。
不只是回答问题，而且要建立生产周期、要不只是回答问题，而且要建立生产周期、产量与需求量、准备费、贮存费之间的关系。
求需求量、准备费、贮存费之间的关系。
<i>姜启源数学建模资料</i>问题分析与思考日需求100件，准备费5000元，贮存费每日每件元。
件准备费日需求元贮存费每日每件1元每天生产一次，每次每天生产一次，每次100件，无贮存费，准备费件无贮存费，准备费5000元。
元每天费用5000元元每天费用10天生产一次，每次天生产一次，天生产一次每次1000件，贮存费件贮存费900+800+…+100 =4500 准备费5000元，总计元，准备费元总计9500元。
元平均每天费用950元元平均每天费用50天生产一次，每次天生产一次，天生产一次每次5000件，贮存费件贮存费4900+4800+…+100 =*****元，准备费元准备费5000元，总计元总计*****元。
姜启源编数学模型第四版
一般模型 x(t) ~甲方兵力，y(t) ~乙方兵力
模型 假设
• 每方战斗减员率取决于双方的兵力和战斗力. • 每方非战斗减员率与本方兵力成正比. • 甲乙双方的增援率为u(t), v(t).
x(t) f (x, y) x u(t), 0
tm~传染病高潮到来时刻
tm
1
ln
1 i0
1
t i 1 ?
(日接触率) tm
病人可以治愈！
第6页/共76页
模型3
传染病无免疫性——病人治愈成 为健康人，健康人可再次被感染. SIS 模型
增加假设 3）病人每天治愈的比例为 ~日治愈率
建模 N[i(t t) i(t)] Ns(t)i(t)t Ni(t)t
di
dt
i(1 i)
i
i[i (1 1 )]
i(0) i0
/
~ 日接触率 1/ ~感染期
~ 一个感染期内每个病人的
有效接触人数，称为接触数.
第7页/共76页
模型3
di/dt
di i[i (1 1 )]
dt
接触数 (感染期内每个
病人的有效接触人数)
i
i
>1
i0
>1
1
1-1/
接触率
N[i(t t) i(t)] [s(t)]Ni(t)t
di si
dt
s(t) i(t) 1
di
i(1 i)
dt
i(0) i0
第5页/共76页
模型2
i
di
i(1 i)
dt
i(0) i0
Logistic 模型
1
i(t)
数学模型第四版课后答案姜启源版
《数学模型》作业答案第二章(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 ∑==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 ππ+=∴《数学模型》作业解答第三章1（2008年10月14日）1. 在3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货批量．证明在不允许缺货模型中结果与原来的一样，而在允许缺货模型中最优订货周期和订货批量都比原来结果减少．解：设购买单位重量货物的费用为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 ==**与不考虑购货费的结果比较，Ｔ、Ｑ的最优结果没有变．02 对于允许缺货模型，每天平均费用为：⎥⎦⎤⎢⎣⎡+-++=kQ Q rT r c r Q c c T Q T C 23221)(221),( 2223322221222TkQ 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与不考虑购货费的结果比较，Ｔ、Ｑ的最优结果减少．2．建立不允许缺货的生产销售存贮模型．设生产速率为常数k ，销售速率为常数r ，r k >．在每个生产周期Ｔ内，开始的一段时间()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 kc 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 Y=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 TT T q T p b p -+---=∂∂βY )(2)8322(22022bp a TT t q T p b p -+---=∂∂βY 0,021=∂∂=∂∂p p Y Y 令, 得到最优价格为: ⎪⎪⎩⎪⎪⎨⎧⎥⎦⎤⎢⎣⎡++=⎥⎦⎤⎢⎣⎡++=)43(21)4(210201T q b a b p T q b a b p ββ 在销售期T 内的总销量为⎰⎰+-=-+-=20221210)(2)()(TTT 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 ββ---+---=Yt 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 ＝702l925002+-=TdT dC以及x=0,y=0组成的凸四边形区域直线l ：20x+30y=c 在可行域内平行移动.易知：当l 过1l 与2l 的交点时，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 x t 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≈≈ ＱQ4．在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节的产量模型不同之处．要获得最大持续产量，应使渔场鱼量2N x >，且尽量接近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）Ex()x f)4(42Nhr r N rh r -=-=∆ ， （1）的解为：2412,1N rNhN x -±=① 当0π∆时，（1）无实根，此时无平衡点； ② 当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 Nx 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 作初等行变换可得：∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−−→−0000011211ΛΛΛΛΛΛΛn b b b A 初等行变换B 这里0≠B .()1=∴B 秩，从而秩()1=A再根据初等行变换与初等矩阵的关系知：存在一个可逆阵P ，使B PA =，于是∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡==--000001121111ΛΛΛΛΛΛΛ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ΛM M M Λ=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=∑∑∑∑∑∑======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-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 ≈ ∴ mn E E 32/'=，当n m φ时，132πm n， ∴ 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份时，收益的期望值最大.数模复习资料第一章1. 原型与模型原型就是实际对象.模型就是原型的替代物.所谓模型, 按北京师范大学刘来福教授的观点：模型就是人们为一定的目的对原型进行的一个抽象.如航空模型、城市交通模型等.模型⎪⎪⎪⎩⎪⎪⎪⎨⎧⎪⎩⎪⎨⎧⎩⎨⎧数学模型如地图、电路图符号模型如某一操作思维模型抽象模型如某一试验装置物理模型如玩具、照片等直观模型形象模型2. 数学模型对某一实际问题应用数学语言和方法,通过抽象、简化、假设等对这一实际问题近似刻划所得的数学结构,称为此实际问题的一个数学模型. 例如力学中著名的牛顿第二定律使用公式22dt xd mF =来描述受力物体的运动规律就是一个成功的数学模型.或又如描述人口()t N 随时间t 自由增长过程的微分方程()()t rN dtt dN =. 3. 数学建模所谓数学建模是指根据需要针对实际问题组建数学模型的过程.更具体地说,数学建模是指对于现实世界的某一特定系统或特定问题,为了一个特定的目的,运用数学的语言和方法,通过抽象和简化,建立一个近似描述这个系统或问题的数学结构(数学模型),运用适当的数学工具以及计算机技术来解模型,最后将其结果接受实际的检验,并反复修改和完善.数学建模过程流程图为：4.数学建模的步骤依次为：模型准备、模型假设、模型构成、模型求解、模型分析、模型检验、模型应用 5.数学模型的分类数学模型可以按照不同的方式分类,常见的有：a. 按模型的应用领域分类 数学模型 ⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧再生资源利用模型水资源模型城镇规划模型生态模型环境模型（污染模型）交通模型人口模型b. 按建模的数学方法分类数学模型 ⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧规划论模型概率模型组合数学模型图论模型微分方程模型几何模型初等数学模型c. 按建模目的来分类 数学模型 ⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧控制模型决策模型优化模型预报模型分析模型描述模型d.层次分析法的基本步骤：1.建立层次结构模型2.构造成对比较阵3.计算权向量并作一致性检验4.计算组合权向量并作组合一致性检验e.n 阶正互反正A 是一致阵的充要条件为A 的最大特征值为nf.正互反阵最大特征根和特征向量的实用算法：幂法、和法、根法4．在“椅子摆放问题”的假设条件中，将四脚的连线呈正方形改为呈长方形，其余条件不变.试构造模型并求解.解：设椅子四脚连线呈长方形ABCD. AB 与CD 的对称轴为x 轴，用中心点的转角θ表示椅子的位置.将相邻两脚A 、B 与地面距离之和记为)(θf ;C 、D 与地面距离之和记为)(θg .并旋转0180.于是，设,0)0(,0)0(=g f φ就得到()()0,0=ππf g φ.数学模型：设()()θθg f 、是[]π2,0上θ的非负连续函数.若[]πθ2,0∈∀,有。
数学模型第四版(姜启源)作业对于6.4节蛛网模型讨论下列问题：【范本模板】
对于6。
4节蛛网模型讨论下列问题：（1）因为一个时段上市的商品不能立即售完，其数量也会影响到下一时段的价格，所以第k+1时段的价格1+k y 由第k+1和第k 时段的数量1+k x 和k x 决定。
如果设1+k x 仍只取决于k y ，给出稳定平衡的条件，并与6.4的结果进行比较。
（2）若除了1+k y 由1+k x 和k x 决定之外，1+k x 也由前两个时段的价格k y 和1-k y 决定，试分析稳定平衡的条件是否还会放宽。
解:（1)设1+k y 由1+k x 和k x 的平均值决定，即价格函数表示为：)2(11k k k x x f y +=++ 则 0),2(0101>-+-=-++ααx x x y y k k k 0),(001>-=-+ββy y x x k k消去y, 得到 012)1(22x x x x k k k +=++++αβαβαβ ,k=1，2，….该方程的特征方程为022=++αβαβλλ与6.4节中 )2(11-++=k k k y y g x 时的特征方程一样， 所以0〈αβ〈2, 即为0p 点的稳定条件。
（2）设 )2(11k k k x x f y +=++ )2(11-++=k k k y y g x ， 则有 0),2(0101>-+-=-++ααx x x y y k k k 0),2(0101>-+=--+ββy y y x x k k k 消去y ，得到0123)1(424x x x x x k k k k +=++++++αβαβαβαβ 该方程的特征方程为02423=+++αβαβλαβλλ令λ=x ,αβ=a ， 即求解三次方程0a 2ax ax 4x 23=+++ 的根 在matlab 中输入以下代码求解方程的根x ：syms x asolve(4*x^3+a*x^2+2＊a*x+a==0,x)解得 1x = （36＊a^2 — 216＊a — a^3 + 24＊3^（1/2）＊（-a^2＊（a — 27））^（1/2)）^（1/3）/12 — a/12 + (a*(a — 24))/（12＊（36*a^2 — 216＊a — a^3 + 24*3^（1/2)＊(-a^2＊（a — 27）)^（1/2))^(1/3）)；2x = -（2＊a*(36*a^2 - 216＊a — a^3 + 24＊3^（1/2）＊（—a^2＊(a - 27））^（1/2）)^（1/3） — 3^（1/2)*a*24＊i — 3^（1/2）*（36＊a^2 — 216＊a — a^3 + 24＊3^(1/2）＊(—a^2*(a — 27））^(1/2))^(2/3）*i - 24＊a + 3^(1/2)＊a^2*i+ (36*a^2 - 216＊a - a^3 + 24＊3^（1/2）*(-a^2＊（a — 27）)^（1/2）)^（2/3) + a^2)/（24*(36*a^2 — 216*a - a^3 + 24＊3^(1/2)＊（-a^2*（a — 27））^（1/2）)^（1/3））；3x =—（2＊a*（36*a^2 - 216*a — a^3 + 24＊3^（1/2)*(-a^2*（a - 27))^（1/2))^(1/3) + 3^(1/2)＊a ＊24*i + 3^（1/2）*(36＊a^2 - 216＊a — a^3 + 24*3^(1/2）＊(-a^2*(a - 27））^（1/2))^(2/3)*i — 24＊a - 3^（1/2)＊a^2＊i + (36＊a^2 - 216＊a - a^3 + 24*3^(1/2）*(-a^2＊（a - 27))^(1/2))^（2/3) + a^2）/（24＊（36＊a^2 — 216＊a — a^3 + 24*3^（1/2）*（—a^2*(a -27）)^(1/2）)^(1/3））；其中1x 为实根，2x 与3x 为一对共轭虚根。
数学模型姜启源第四版答案
数学模型姜启源第四版答案【篇一：姜启源数学模型课后答案(3版)】t>第二章(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的分配方案，3p1?235,p2?333,p3?432, ?pi?1000.i?1方法一（按比例分配）q1?p1n3?2.35,q2?p2n3?3.33, q3?p3n3?4.32?i?1pi?i?1pii?1pi分配结果为： n1?3, n2?3, n3?4 方法二（q值方法）9个席位的分配结果（可用按比例分配）为：n1?2,n2?3, n3?4第10个席位：计算q值为 q1?23522?3?9204.17, q2?33323?4?9240.75, q3?43224?5?9331.2q3最大，第10个席位应给c.分配结果为 n1?2,n2?3,n3?5方法三（d’hondt方法）此方法的分配结果为：n1?2,n2?3,n3?5此方法的道理是：记pi和ni为各宿舍的人数和席位（i=1,2,3代表a、b、c宿舍）.pinipinipini是每席位代表的人数，取ni?1,2,?,从而得到的近.中选较大者，可使对所有的i,尽量接再考虑n?15的分配方案，类似地可得名额分配结果.现将3种方法两次分配的结果列表如下：2．试用微积分方法，建立录像带记数器读数n与转过时间的数学模型. 解：设录像带记数器读数为n时，录像带转过时间为t.其模型的假设见课本.考虑t到t??t时间内录像带缠绕在右轮盘上的长度，可得vdt?(r?wkn)2?kdn,两边积分，得 ?vdt?2?k?(r?wkn)dntn22) ?t?2?rkvn??wkv2n.2第二章(2)（2008年10月9日）15．速度为v的风吹在迎风面积为s的风车上，空气密度是? ，用量纲分析方法确定风车获得的功率p与v、s、?的关系.解: 设p、v、s、?的关系为f(p,v,s,?)?0，其量纲表达式为: [p]=ml2t?3, [v]=lt量纲矩阵为：?2?1????3(p)10?1(v)200(s)?3?(l)?1(m)? ?0?(t)(??1,[s]=l2,[?]=ml?3,这里l,m,t是基本量纲.a=齐次线性方程组为：?2y1?y2?2y3?3y4?0??0 ?y1?y4??3y?y?012?它的基本解为y?(?1,3,1,1)由量纲pi定理得 ??p?1v3s1?1，?p??v3s1?1 ，其中?是无量纲常数. 16．雨滴的速度v与空气密度?、粘滞系数?和重力加速度g 有关，其中粘滞系数的定义是：运动物体在流体中受的摩擦力与速度梯度和接触面积的乘积成正比，比例系数为粘滞系数，用量纲分析方法给出速度v的表达式.解：设v,?,?,g 的关系为f(v,?,?,g)=0.其量纲表达式为[v]=lm0t-1，[?]=l-3mt0，[?]=mlt-2（lt-1l-1）-1l-2=mll-2t-2t=l-1mt-1，[g]=lm0t-2,其中l，m，t是基本量纲.量纲矩阵为?1?0?a=???1(v)?310(?)?11?1(?)1?(l)?0(m)? ?2?(t)?(g)齐次线性方程组ay=0 ，即? y1-3y2-y3?y4?0??0 ?y2?y3?-y-y-2y?034?1的基本解为y=(-3 ,-1 ,1 ,1) 由量纲pi定理得 ??v?*?3?1?g. ?v???g?，其中?是无量纲常数.16．雨滴的速度v与空气密度?、粘滞系数?、特征尺寸?和重力加速度g有关，其中粘滞系数的定义是：运动物体在流体中受的摩擦力与速度梯度和接触面积的乘积成正比，比例系数为粘滞系数，用量纲分析方法给出速度v的表达式.解：设v,?,?,?，g 的关系为f(v,?,?,?,g)?0.其量纲表达式为[v]=lmt，[?]=lmt，[?]=mlt（ltl）l=mlltt=lmt，[?]=lmt ，[g]=lmt 其中l，m，t是基本量纲. 量纲矩阵为?1?0a=????1(v)100(?)?310(?)?11?1(?)1?(l)?0(m)? ?2?(t)?(g)-1-3-2-1-1-1-2-2-2-1-1000-2齐次线性方程组ay=0 即 ?y1?y2?3y3?y4?y5?0?y3?y4?0 ???y1?y4?2y5?0?的基本解为??y1?(1,???y2?(0,??)2 1,?1,1,?)22231,0,0,?1得到两个相互独立的无量纲量 ??1?v??1/2g?1/2??3/2?1?1/2??g??2??即 v??g?1,?3/2?g1/2??1??2. 由?(?1,?2)?0 , 得 ?1??(?2)3/2?1?1? ??g?(??g1/2??1) , 其中?是未定函数.20.考察阻尼摆的周期，即在单摆运动中考虑阻力，并设阻力与摆的速度成正比.给出周期的表达式，然后讨论物理模拟的比例模型，即怎样由模型摆的周期计算原型摆的周期. 解：设阻尼摆周期t,摆长l, 质量m,重力加速度g，阻力系数k的关系为f(t,l,m,g,k)?0其量纲表达式为：[t]?lmt,[l]?lmt,[m]?lmt,[g]?lmt?2,[k]?[f][v]?1?mlt?2(lt?1)?1?lmt0?1，其中l，m，t是基本量纲.量纲矩阵为?0?0a=???110001010?20?(l)?1(m)? ??1?(t)(t)(l)(m)(g)(k)齐次线性方程组y2?y4?0??y3?y5?0 ??y?2y?y?045?1的基本解为11?y?(1,?,0,,0)?122 ?11?y2?(0,,?1,?,1)22?得到两个相互独立的无量纲量?tl?1/2g1/2??1?1/2?1?1/2k??2?lmg∴t?lg?1, ?1??(?2), ?2?klmg1/21/2klmg1/21/2∴t?lg() ，其中?是未定函数 .考虑物理模拟的比例模型，设g和k不变，记模型和原型摆的周期、摆长、质量分别为l?gkl?1/21/2t,t；l,l;m,m. 又t??(m?g)当无量纲量?l?l时，就有t?t?l?g?gl?l?l.《数学模型》作业解答第三章1（2008年10月14日）1. 在3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货批量．证明在不允许缺货模型中结果与原来的一样，而在允许缺货模型中最优订货周期和订货批量都比原来结果减少．【篇二：数学建模陈东彦版课后答案】t>2.9-3.7 3.6-5.144.1-7.14.4-7.35.9-11.1 5.1-9.1 6.5-4.7 6.10-4.14第1章建立数学模型1.1 在稳定的椅子问题中，如设椅子的四脚连线呈长方形，结论如何？（稳定的椅子问题见姜启源《数学模型》第6页）1.2 在商人们安全过河问题中，若商人和随从各四人，怎样才能安全过河呢？一般地，有n名商人带n名随从过河，船每次能渡k人过河，试讨论商人们能安全过河时，n与k应满足什么关系。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
di
dt
i(1 i)
i
i[i (1 1 )]
i(0) i0
/
~ 日接触率 1/ ~感染期
~ 一个感染期内每个病人的
有效接触人数，称为接触数.
第7页/共76页
模型3
di/dt
di i[i (1 1 )]
dt
接触数 (感染期内每个
病人的有效接触人数)
i
i
>1
i0
>1
1
1-1/
传染病有免疫性——病人治愈 后即移出感染系统，称移出者.
SIR模型
假设 1）总人数N不变，病人、健康人和移
出者的比例分别为 i(t), s(t), r(t).
2）病人的日接触率 , 日治愈率, 接触数 = /
建模 s(t) i(t) r(t) 1
需建立 i(t ), s(t ), r (t ) 的两个方程.
消去dt
/
di
ds
1
s
1
i ss0 i0
相轨线
相轨线 i(s) 的定义域
i(s)
(s0
i
i0 )
s
1
ln
s s0
D {(s,i) s 0, i 0, s i 1} 1
在D内作相轨线 i(s)
的图形，进行分析
D
O
第12页/共76页
s
1
模型4 相轨线 i(s) 及其分析
SIR模型
di
• 降低 s0
提高 r0
的估计
s0 i0 r0 1
s0
i0
s
1
ln s s0
0
忽略i0
群体免疫
ln s0 ln s
s0 s
第14页/共76页
模型4
预防传染病蔓延的手段
• 降低日接触率 • 提高日治愈率 • 提高移出比例r0
i0 di/dt < 0
O
1-1/ 1 i
i0
i()
1
1
,
O
1
tO
1 i(t)单调下降
t
0,
1 感染期内有效接触使健康者感
>1, i0< 1-1/
染的人数不超过原有的病人数
i(t)按S形曲线增长
接触数 =1 ~ 阈值
模型2(SI模型)如何看作模型3(SIS模型)的特例
第8页/共76页
模型4
第9页/共76页
模型4
SIR模型
N[i(t t) i(t)] Ns(t)i(t)t Ni(t)t
N[s(t t) s(t)] Ns(t)i(t)t
di
dt
si i
ds
dt
si
i(0) i0 , s(0) s0
i0 s0 1 (通常r(0)=r0很小)
无法求出 i(t), s(t)
1s
P1: s0>1/ i(t)先升后降至0 传染病蔓延
P2: s0<1/ i(t)单调降至0
传染病不蔓延
第13页/共76页
1/ ~ 阈值
模型4
预防传染病蔓延的手段
SIR模型
传染病不蔓延的条件——s0<1/ • 提高阈值 1/ 降低 (=/)
,
(日接触率) 卫生水平
(日治愈率) 医疗水平
第五章 微分方程模型
5.1 传染病模型 5.2 经济增长模型 5.3 正规战与游击战 5.4 药物在体内的分布与排除 5.5 香烟过滤嘴的作用 5.6 人口的预测和控制 5.7 烟雾的扩散与消失 5.8 万有引力定律的发现
第1页/共76页
动态 模型
• 描述对象特征随时间(空间)的演变过程. • 分析对象特征的变化规律. • 预报对象特征的未来性态. • 研究控制对象特征的手段.
tm~传染病高潮到来时刻
tm
1
ln
1 i0
1
t i 1 ?
(日接触率) tm
病人可以治愈！
第6页/共76页
模型3
传染病无免疫性——病人治愈成 为健康人，健康人可再次被感染. SIS 模型
增加假设 3）病人每天治愈的比例为 ~日治愈率
建模 N[i(t t) i(t)] Ns(t)i(t)t Ni(t)t
dt
si
i
ds
dt
si
i(0) i0 , s(0) s0
di
ds
1
s
1
i ss0 i0
i
1
i(s)
(s0
i0
)
s
1
ln
s s0
D
P4
s(t)单调减相轨线的方向 im s 1/ , i im t , i 0
P2
P1
P3
s满足
s0
i0 s
1
ln
s s0
0
s O
S0 1/ s0
基本 不是从医学角度分析各种传染病的特殊机理, 方法 而是按照传播过程的一般规律建立数学模型.
第3页/共76页
模型1
假设 建模
已感染人数 (病人) i(t)
每个病人每天有效接触
(足以使人致病)人数为
i(t t) i(t) i(t)t
di i
dt i(0) i0
i(t) i0et
ti ?
微分 方程 建模
• 根据函数及其变化率之间的关系确定函数. • 根据建模目的和问题分析作出简化假设. • 按照内在规律或用类比法建立微分方程.
第2页/共76页
5.1 传染病模型
背景 与 问题
传染病的极大危害(艾滋病、SARS、)
• 描述传染病的传播过程. • 分析受感染人数的变化规律. • 预报传染病高潮到来的时刻. • 预防传染病蔓延的手段.
0.4 i(t)
0.2
0
0
10
20
30
40
50
0.4
i
相轨线i(s)
0.3
0.2
0.1
s
P0
0
0
0.2 0.4 0.6 0.8
1
i(t)从初值增长到最大; t, i0. s(t)单调减; t, s0.04.
第11页/共76页
模型4 SIR模型的相轨线分析
di
dt
si
i
ds
dt
si
i(0) i0 , s(0) s0
接触率
N[i(t t) i(t)] [s(t)]Ni(t)t
di si
dt
s(t) i(t) 1
di
i(1 i)
dt
i(0) i0
第5页/共76页
模型2
i
di
i(1 i)
dt
i(0) i0
Logistic 模型
1
i(t)
1
1/2
1
1 i0
1
et
i0
0
tm
t
t=tm, di/dt 最大
的解析解
先做数值计算, 再在相平面上研 究解析解性质
第10页/共76页
模型4 SIR模型的数值解
di dt
si
i,
i(0)
i0
ds
dt
si,
s(0) s0
设=1, =0.3, i0=0.02, s0=0.98, 用
MATLAB计算作图i(t), s(t)及i(s)
1
0.8 s(t)
0.6
若有效接触的是病人, 则不能使病人数增加
必须区分已感染者(病人) 和未感染者(健康人)
第4页/共76页
ቤተ መጻሕፍቲ ባይዱ
模型2
假设
建模
区分已感染者(病人)和未感染者(健康人)
1）总人数N不变，病人和健康
人的 比例分别为 i(t), s(t) .
SI 模型
2）每个病人每天有效接触人数 ~ 日
为, 且使接触的健康人致病.