(整理)作业1数学建模,姜启源版.
- 格式:doc
- 大小:847.50 KB
- 文档页数:15
综合题目参考答案1. 赛程安排(2002年全国大学生数学建模竞赛D 题)(1)用多种方法都能给出一个达到要求的赛程.(2)用多种方法可以证明支球队“各队每两场比赛最小相隔场次n r 的上界”(如=5时上界为1)是n ⎥⎦⎤⎢⎣⎡-23n ,如: 设赛程中某场比赛是,i j 两队, 队参加的下一场比赛是,两队(≠i i k k j ),要使各队每两场比赛最小相隔场次为r ,则上述两场比赛之间必须有除i ,j ,以外的2k r 支球队参赛,于是,注意到32+≥r n r 为整数即得⎥⎦⎤⎢⎣⎡-≤23n r . (3)用构造性的办法可以证明这个上界是可以达到的,即对任意的编排出达到该上界的赛程.如对于n =8, =9可以得到: n n 1A 2A 3A 4A 5A 6A 7A 8A 每两场比赛相隔场次数 相隔场次总数1A× 1 5 9 13 17 21 25 3,3,3,3,3,3 18 2A 1 × 20 6 23 11 26 16 4,4,4,3,2,2 193A 5 20 × 24 10 27 15 2 2,4,4,4,3,2 19 4A 9 6 24 × 28 24 3 19 2,2,4,4,4,3 19 5A 13 23 10 28 × 4 18 7 2,2,2,4,4,4 18 6A 17 11 27 14 4 × 8 22 3,2,2,2,4,4 177A 21 26 15 3 18 8 × 12 4,3,2,2,2,4 178A25 16 2 19 7 22 12 × 4,4,3,2,2,2 17w w w .k h d a w .c o m 课后答案网1A 2A 3A 4A 5A 6A 7A 8A 9A 每两场比赛相隔场次数 相隔场次总数1A× 36 6 31 11 26 16 21 1 4,4,4,4,4,4,4, 28 2A 36 × 2 27 7 22 12 17 32 4,4,4,4,4,4,3 27 3A 6 2 × 35 15 30 20 25 10 3,3,4,4,4,4,4 26 4A 31 27 35 × 3 18 8 13 234,4,4,4,3,3,3 25 5A 11 7 15 3 × 34 24 29 193,3,3,3,4,4,4 24 6A 26 22 30 18 34 × 4 9 144,4,3,3,3,3 23 7A16 12 20 8 24 4 × 33 28 3,3,3,3,3,3,4 22 8A21 17 25 13 29 9 33 × 53,3,3,3,3,3,3, 21 9A 1 32 10 23 19 14 28 5 × 3,4,3,4,3,4,3 24 可以看到, =8时每两场比赛相隔场次数只有2,3,4, =9时每两场比赛相隔场次数只有3,4,以上结果可以推广,即为偶数时每两场比赛相隔场次数只有n n n 22-n ,12-n ,2n ,n 数时只有为奇23-n ,21-n . 量赛程优劣其他指标如(4)衡的平均相隔场次 记第i 队第j 个ij c ,2,2,1,,,2,1-==n j n i ,间隔场次数为则平均相隔场次为∑∑=n i 1-=n r 21 =-j n n 1)2(ij c r 是赛程整体意义下的指标,它越大越好.可以计算=8,=9的n n r ,并讨论它是否达到上界. 相隔场次的最大偏差 定义||,r c Max f ij j i -=∑---=2)2(|n r n c Max g =1|j ijw w w .k h d a w .c o m 课后答案网f 为整个赛程相隔场次的最大偏差, 为球队之间相隔场次的最大偏差,它们都是越小越好.可以计算=8,=9的,g ,并讨论它是否达到上界.g n n f 参考文献工程数学学报第20卷第5期20032. 影院座位设计建立满意度函数),(βαf ,可以认为α和β无关, ()()βαβαh g f -=),(,g ,取尽量简单的形式,h 如αα=)(g ;0)(=βh (),030≤β0)(h h =β)30(0>β.(1)可将作为必要条件,以030≤βα最大为最佳座位的标准.在上图中以第1排座位为坐标原点建立坐标轴x ,可以得到 ⎪⎭⎫ ⎝⎛+----⎪⎭⎫ ⎝⎛+--=⎪⎭⎫ ⎝⎛+--=d x x h c H d x x c H d x x c H θθαθβtan arctan tan arctan ,tan arctan β是x 的减函数.可得x ≈1.7m,即第3(或4)排处.又通过计算或分析可知030=βα也是x 的减函数,所以第3(或4)排处是最佳座位.(2)设定一个座位间隔(如0.5m), l x 从0(或处)到030≤βd D -按离散,对于计算l )20~0(00θα的平均值,得时其值最大. 020=θ(3)可设地板线是x 的二次曲线,寻求,b 使2bx ax +a α的平均值最大. 实际上,还应考虑前排不应挡住后排的视线. 3.节水洗衣机(1996年全国大学生数学建模竞赛B 题) 该问题不要求对洗衣机的微观机制(物理、化学方面)深入研究,只需要从宏观层次去把握.宏观上洗衣的基本原理是用洗涤剂通过漂洗把吸附在衣物上的污物溶于水中,再脱去污水带走污物;洗衣的过程是通过“加水——漂洗——脱水”程序的反复运行,使残留在衣物的污物越来越少,直到满意的程度;洗涤剂也是不希望留在衣物上的东西,可将“污物”定义为衣物上原有污物与洗涤剂的总和. w w w .k h da w .c o m 课后答案网假设每轮漂洗后污物均匀地溶于水中;每轮脱水后衣物含水量为常数.~初始污水量,第轮加水量,~第k 轮脱水量c 0x ~k u k k x ),,2,1( =k .设每轮脱水前后污物在水中的浓度不变.于是cx c u x u c x n n n =+==--111,,, c x 2c x +21u x 10, 得到)()(210c u c u u c x x n n n ++= . 在最终污物量与初始污物量之比小于给定的清洁度条件下,求各轮加水量,使总用水量最小,即0/x x n k u ),,1(n k =∑=nk k u u Min k 1()ε<++)(..21c u c u u c t s n n 等价于)()(21c u c u u Min n u k +++++ α=++)()(..21c u c u u t s na 为常数可得c u c u u n +==+= 21,即第轮加水量n ~2u u k =(常数),第1轮加水量.c u u +=1令,问题简化为cx u =nx Min u n , ε<⎪⎭⎫ ⎝⎛+n x t s 11.. 其解为,即,而0→x 0→u ∞→n n .这与实际上是不合理的.应该加上对u 的限制:.则得n ,其中 21v u v ≤≤max min n n ≤≤max min n n ≤≤,1+)/1ln(2min ⎥⎦⎤⎢⎣⎡+=c v n αn 这样,为有限的几个数,可一一比较,具体数据计算从略.参考文献:《数学的实践与认识》第27卷第1期,1997w w w .k h d a w .c o m 课后答案网4.教师工资调整方案(1995年美国大学生数学建模竞赛B 题)题目对职称提升年限表述得不甚清楚(如未提及助理教授的提升),教龄也未区分是什么职称下工作的年限,所以应该作出一些相应的简化假设.按所给信息,工资仅取决于职称和教龄.建立新方案的一种办法是将职称折合成教龄,如定义x=教龄t+7×k (对于讲师、助理教授、副教授、教授,k 分别取值0,1,2,3),然后寻求工资函数I(x),使之满足题目的要求,如I(0)=27000,I(7)=32000等,以及x 较大时022<dxI d .另一种办法是职称、教龄分别对待,工资函数J(k,t)从多种函数中选择,如最简单的线性函数J(k,t)=k k k k b a t b a ,,+(k=0,1,2,3)根据一定条件确定.按照第一种办法得到的新工资方案,以职称和教龄综合指标为x 的教师的工资都应为I(x),而人们的目前工资会低于或高于它.根据题目要求,高工资不应降低,低工资则应逐渐提高,尽快达到理想值I(x).需要做的只是根据每人(目前)工资与(理想值的)差额,制定学校提供的提薪资金的分配方案.它应该是简单、合理、容易被人接受的. 按以上原则可以建立不同的模型,应通过检验比较其恶劣.检验可基于题目所给数据,按照提薪计划运行若干年,考察接近理想方案的情况,即用过渡时期的情况检验模型;也可进行随机模拟,按照一定规则随机产生数据(可以包括聘用、提职、解聘、退休的人数和时间等),再按照提薪计划运行,考察接近理想方案的情况.参考文献:叶其孝,《大学生数学建模竞赛辅导教材》(四),湖南教育出版社,20015. 一个飞行管理问题(1995年全国大学生数学建模竞赛A 题)设为第i 架飞机与第j 架飞机的碰撞角(即ij a )8arcsin(ij ij r a =其中为这两架飞机连线的长度),ij r ij β为第i 架飞机相对于第j 架飞机的相对速度(矢量)与这两架飞机连线(从i 指向j 的矢量)的夹角(以连线矢量为基准,逆时针方向为正,顺时针方向为负),i θ为第架飞机飞行方向角调整量. 本问题中的优化目标函数可以有不同的形式:如使所有飞机的最大调整量最小;所有飞机的调整量绝对值之和最小等.以所有飞机的调整量绝对值之和最小,可以得到如下的数学规划模型:w w w .k h d a w .c o m 课后答案网∑=61i i Min θ s.t. ,)(21ij j i ij a >++θθβ j i j i ≠=,6,,1,30≤i θ , 6,,1 =i 为了利用LINGO 求解这个数学规划模型,可以首先采用其他数学软件计算出ij α和ij β.其实,ij α和ij β也是可以直接使用LINGO 来计算的,这相当于解关于ij α和ij β的方程,只是解方程并非LINDO 软件的特长,这里我们作为一个例子,看看如何利用LINGO 计算ij α,可输入如下模型到LINGO 求解ij α:MIDEL :1]SETS:2] PLANE/1..6/:x0,y0; 3] link(plane,plane):alpha,sin2: 4]ENDSETS5] @FOR(LINK(I,J)|I#NE#J:6] sin2(I,J)=64/((X0(I)-X0(J))*(X0(I)-X0(J))+7] (Y0(I)-Y0(J))*(Y0(I)-Y0(J)));8] );9] @FOR(LINK(I,J)|I#NE#J: 10] (@SIN(alpha*3.14159265/180.0))^2=SIN2; 11] ); 12]DATA:13] X0=150,85,150,145,130,0; 14] Y0=140,85,155,50,150,0; 15]endataEND计算结果如下:w w w .k h d a w .c o m 课后答案网ij a j=1 2 3 4 5 6i=1 0.000 0 5.3912 32.231 05.091 8 20.963 4 2.234 5 2 5.391 2 0.000 0 4.8046.613 5 5.807 9 3.815 9 3 32.231 0 4.804 0 0.0004.364 7 22.833 7 2.125 5 45.091 86.613 5 4.36470.000 0 4.4.537 2.989 8 5 20.963 4 5.807 922.8337 4.537 70.000 0 2.309 8 6 2.234 5 3.815 9 2.125 5 2.989 82.309 80.000 0 ij β也可类似地利用LINGO 求得,计算结果如下: ij β j=1 2 3 4 5 6 i=1 0.000 0 109.263 6 -128.250 0 24.1798173.065 1 14.474 9 2 109.263 6 0.000 0-88.871 1 -42.2436-92.304 8 9.000 03 -128.250 0 -88.871 1 0.000 012.4763-58.786 2 0.310 84 24.179 8 -42.243 6 12.476 30.000 0 5.969 2-3.525.65 173.065 1 -92.304 8 -58.78625.969 20.000 0 1.914 4614.474 9 9.000 00.310 8-3.5256 1.914 4 0.000 0w w w .k h d a w .c o m 课后答案网于是,该飞机管理的数学规划模型可如下输入LINGO 求解:MODEL:1]SETS2] plane/1..6/:cita:3] link(plane,plane):alpha,beta;4]ENDSETS5] min=@sum(plane:@abs(cita));6] @for(plane(I):7] @bnd(-30,cita(I),30);8] );9] @fpr(link(I,j)|I#NE#J:10] @ABS(beta(I,J)+0.5*cit(I)+0.5*cita(J))11] >alpha(I,J);12] );13]DATA:14] A;[JA=0.000 0 5.391.2….. …2.309 8 0.000 020] ;21] BETA=0.000 010 9.263 6………1.914 4 0.000 027] ;28]enddata END[注] alpha,beta 中数据略去,见上面表格. 求解结果如下: OPTIMUM FOUND AT STEP 197 SOLUTION OBJECTIVE V ALUE= 3.630 V ARIABLE V ALUE REDUCED COST CITA(1) 0.2974033E-06 -1.000 000 CITA(2) -0.1424833E-05 -0.715 033 4 w w w .k h d a w .c o m 课后答案网CITA(3) 2.557 866 1.000 000 CITA(4) -0.3856641E-04 0.0000000E+00CITA(5) 0.2098838E-05 -1.000 000CITA(6) 1.071 594 0.0000000E+00………. (以下略)由此可知最优解为: (其它调整角度为0). ︒︒≈≈07.1,56.263θθ 评注:如果将目标改为最大调整量最小,则可进一步化简得到线形规划模型,也可用LINDO 或LINGO 求解.参考文献:《数学的实践与认识》第26卷第1期,19966. 降落伞的选择这个优化问题的决策变量是降落伞数量n 和每一个伞的半径r ,可先将n 和r 看作连续变量,建立优化模型,求得最优解后,再按题目要求作适当调整. 目标函数之降落伞的费用,可以根据表1数据拟合伞面费用与伞的半径r 的关系。
姜启源数学建模资料简单的优化模型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元,总计元总计*****元。
1、举出两三个实例说明建立数学模型的必要性,包括实际问题的背景,建模目的,需要大体上什么样的模型以及怎样应用这种模型等。
2、从下面不大明确的叙述中确定要研究的问题,要考虑哪些有重要影响的变量:(1)一家商场要建一个新的停车场,如何规划照明设施;(2)一农民要在一块土地上作出农作物的种植规划;(3)一制造商要确定某种产品的产量及定格;(4)卫生部门要确定一种新药对某种疾病的疗效;(5)一滑雪场要进行山坡滑道和上山缆车的规划。
3、怎样解决下面的实际问题。
包括需要哪些数据资料,要作些什么观察、实验以及建立什么样的数学模型等。
(1)估计一个人体内血液的总量;(2)为保险公司制定人寿保险金计划(不同年龄的人应缴纳的金额和公司赔偿的金额);(3)估计一批日关灯管的寿命;(4)确定火箭发射至最高点所需的时间;(5)决定十字路口黄灯亮的时间长度;(6)为汽车租凭公司制订车辆维修、更新和出租计划;(7)一高层办公楼有4部电梯,早晨上班时间非常拥挤,试制订合理的运行计划。
4、在1.3节“椅子能在不平的地面上放稳吗”的假设条件中,将四脚的连线呈正方形改为呈长方形,其余不变。
试构造模型并了解。
5、模仿1.4节商人过河问题中的状态转移模型,作下面这个众所周知的智力游戏:人带着猫、鸡、米过河,船除了需要人划之外,至多能载猫、鸡、米三者之一,而当下人不在场时猫要吃鸡、鸡要吃米。
试设计一个安全过河的方案,并使渡河次数尽量地最少。
6、利用1.5节表1和表3给出的1790—2000年美国实际人口资源建立下列模型:(1)分段的指数增长模型。
将时间分为若干段,分别确定增长率r ;(2)阻滞增长模型。
换一种方法确定固有增长率r 和最大容量m x 。
7、说明1.5节中Logistic 模型(9)可以表为0()()1mr t t x x t e --=+,其中是人口增长出现拐点的时刻,并说明0t 与,r ,m x 的关系。
8、假定人口的增长服从这样的规律:t 时刻的人口为()x t ,到t t +∆时间内人口的增量与()m x x t 成正比(其中m x 为最大容量)。
数学模型课后答案姜启源【篇一:姜启源《数模》习题选解】方案模型构成:以阈值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个共顶点。
第四版姜启源数学模型复习总结第1章:了解模型的概念与分类,熟练掌握数学模型的定义,数学模型的重要应用,建模的重要例子-指数模型,Logist模型。
建模的一般方法及其在建模中的应用。
建模的一般步骤(每步的主要内容与问题)。
建模的全过程(框图)4个环节的含义。
模型的特点(技艺性)。
模型分类(表现特征),建模中的能力培养。
数学建模实例的建模思想及其步骤§1 数学模型的概念:模型:模型是为了一定目的,对客观事物的一部分信息进行简缩、抽象、提炼出来的原型的替代物。
模型的分类:具体模型(或物质模型,实的),包括直观模型,物理模型。
抽象模型(或理想模型,虚的),包括思维模型,符号模型,数学模型。
数学模型:对于一个现实对象,为了一个特定目的,根据其内在规律,作出必要的简化假设,运用适当的数学工具,得到的一个数学结构。
1-1-1 模型是为了特定的目的,将原型的()而得到的原型替代物。
1-1-2数学模型可以描述为:对于一个现实对象,()。
1-1-3 关于数学模型的如下论述中正确的是()A。
数学模型是以现实世界的特定问题为研究对象。
B。
数学模型只是对实际问题的近似表示,其中包含一些简化假设。
C。
数学模型表示是某一特定问题的内在规律的数学表示,是以方程和函数关系表示的数学结构。
D。
数学模型是现实问题的真实的描述,不能做任何假设和简化。
1-1-4 关于数学建模的如下论述中正确的是()A。
数学模型和数学建模是完全相同的概念。
B。
数学建模是一个全过程,包括表述、求解、解释和验证四个环节。
C。
数学建模全过程涉及两个世界是现实世界和虚拟世界,涉及的“双向翻译”是同声翻译和文献翻译。
D.数学建模过程是一个从理论-实践-再理论-再实践不断改进的过程。
§2 建模的重要意义(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的分配方案, 方法一(按比例分配)分配结果为: 4 ,3 ,3321===n n n 方法二(Q 值方法)9个席位的分配结果(可用按比例分配)为: 第10个席位:计算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 iin p 尽量接近. 再考虑15=N 的分配方案,类似地可得名额分配结果.现将3种方法两次分配的结果列表如下:2. 试用微积分方法,建立录像带记数器读数n 与转过时间的数学模型. 解: 设录像带记数器读数为n 时,录像带转过时间为t.其模型的假设见课本.考虑t 到t t ∆+时间内录像带缠绕在右轮盘上的长度,可得,2)(kdn wkn r vdt π+=两边积分,得 ⎰⎰+=ntdn wkn r k vdt 0)(2π《数学模型》作业解答第三章1(2008年10月14日)1. 在3.1节存贮模型的总费用中增加购买货物本身的费用,重新确定最优订货周期和订货批量.证明在不允许缺货模型中结果与原来的一样,而在允许缺货模型中最优订货周期和订货批量都比原来结果减少.解:设购买单位重量货物的费用为k ,其它假设及符号约定同课本.01 对于不允许缺货模型,每天平均费用为:令0=dTdC, 解得 rc c T 21*2= 由rT Q = , 得212c rc rT Q ==** 与不考虑购货费的结果比较,T、Q的最优结果没有变.02 对于允许缺货模型,每天平均费用为:令⎪⎪⎩⎪⎪⎨⎧=∂∂=∂∂00QC TC, 得到驻点:与不考虑购货费的结果比较,T、Q的最优结果减少.2.建立不允许缺货的生产销售存贮模型.设生产速率为常数k ,销售速率为常数r ,r k >.在每个生产周期T内,开始的一段时间()00T t <<一边生产一边销售,后来的一段时间)(0T t T <<只销售不生产,画出贮存量)(t g 的图形.设每次生产准备费为1c ,单位时间每件产品贮存费为2c ,以总费用最小为目标确定最优生产周期,讨论r k >>和r k ≈的情况. 解:由题意可得贮存量)(t g 的图形如下:⋅-TT r k c 022)( 又 ∴ T =0kTT r k r 2)(⋅-=于是不允许缺货的情况下,生产销售的总费用(单位时间内)为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≈>>*时当 . 相当于不考虑生产的情况. ∞→≈*,T r k 时当 . 此时产量与销量相抵消,无法形成贮存量.第三章2(2008年10月16日)3.在3.3节森林救火模型中,如果考虑消防队员的灭火速度λ与开始救火时的火势b 有关,试假设一个合理的函数关系,重新求解模型.解:考虑灭火速度λ与火势b 有关,可知火势b 越大,灭火速度λ将减小,我们作如下假设: 1)(+=b kb λ, 分母∞→→+λ时是防止中的011b 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 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 T T t <<<<220和两段,每段的价格固定,记作21,p p .求21,p p 的最优值,使销售期内的总利润最大.如果要求销售期T 内的总售量为0Q ,再求21,p p 的最优值.解:按分段价格,单位时间内的销售量为又 t q t q β+=0)(.于是总利润为=22)(022)(20222011T T t 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 ββ---+--- 0,021=∂∂=∂∂p p 令, 得到最优价格为: 在销售期T 内的总销量为 于是得到如下极值问题: 利用拉格朗日乘数法,解得: 即为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 =70以及x=0,y=0 直线l :20x+30y=c 在可行域内 平行移动.易知:当l 过1l 与2l 1l x S 取最大值.由⎩⎨⎧=+=+7045202y x y x 解得⎩⎨⎧==510y x此时 m ax S =2053010⨯+⨯=350(元)2. 某厂拟用集装箱托运甲乙两种货物,每箱的体积、重量以及可获利润如下表:已知这两种货物托运所受限制是体积不超过24立方米,重量不超过13百斤.试问这两种货物各托运多少箱,使得所获利润最大,并求出最大利润.解:设甲货物、乙货物的托运箱数分别为1x ,2x ,所获利润为z .则问题的数学模型可表示为这是一个整线性规划问题. 用图解法求解. 可行域为:由直线1352:212=+x x l 及0,021==x x 组成直线 c x x l =+211020:在此凸四边形区域内平行移动.易知:当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)可写成(1).s s(t) .s(t) .100≤∴单调减少由若σs (2)().00.1-s ,1,1dtdit s s σσσ从而则若 4.在5.3节正规战争模型(3)中,设乙方与甲方战斗有效系数之比为.4=ba初始兵力00y x 与相同.(1) 问乙方取胜时的剩余兵力是多少,乙方取胜的时间如何确定.(2) 若甲方在战斗开始后有后备部队以不变的速率r 增援,重新建立模型,讨论如何判断双方的胜负.解:用()()t y t x ,表示甲、乙交战双方时刻t 的士兵人数,则正规战争模型可近似表示为:现求(1)的解: (1)的系数矩阵为⎥⎦⎤⎢⎣⎡--=00b a A()()()tab tab eC e C t y t x -⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛∴1212121的通解为.再由初始条件,得又由().1aybx dx dy =可得其解为 ()3 ,202022 bx ay k k bx ay -==-而 (1) ()().231000202011y a b y a bx ay akt y t x =-=-===时,当 即乙方取胜时的剩余兵力数为.230y 又令().0222,01100001=-⎪⎭⎫⎝⎛++⎪⎭⎫⎝⎛-=t ab t ab e y x e y x t x )得由(注意到000020022,1x y y x e y x t ab -+==得. .43ln ,3121bt e t ab =∴=∴ (2) 若甲方在战斗开始后有后备部队以不变的速率r 增援.则().,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 +⎪⎭⎫ ⎝⎛- 亦即第五章2(2008年11月14日)6. 模仿5.4节建立的二室模型来建立一室模型(只有中心室),在快速静脉注射、恒速静脉滴注(持续时间为τ并画出血药浓度曲线的图形.解: 设给药速率为()(),,0t x t f 中心室药量为(1)快速静脉注射: 设给药量为,0D 则()0t f =(2)恒速静脉滴注(持续时间为τ): 设滴注速率为()(),00,000==C k t f k ,则解得 (3) 口服或肌肉注射: ()(),解得)式节(见134.5010010t k e D k t f -=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 eba v aw Q v bl a vl β ()10/10==l M w 其中,(2) 对于一支不带过滤嘴的香烟,全部吸完的毒物量为⎪⎪⎭⎫⎝⎛-=-vbla eb a v aw Q '103‘ 只吸到1l 处就扔掉的情况下的毒物量为⎪⎪⎭⎫ ⎝⎛-=--v bl a vbl e e ba vaw Q 1'21'04 4.在5.3节正规战争模型(3)中,设乙方与甲方战斗有效系数之比为.4=ba初始兵力00y x 与相同.(1) 问乙方取胜时的剩余兵力是多少,乙方取胜的时间如何确定.(2) 若甲方在战斗开始后有后备部队以不变的速率r 增援,重新建立模型,讨论如何判断双方的胜负.解:用()()t y t x ,表示甲、乙交战双方时刻t 的士兵人数,则正规战争模型可近似表示为:现求(1)的解: (1)的系数矩阵为⎥⎦⎤⎢⎣⎡--=00b a A()()()tab tab eC e C t y t x -⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛∴1212121的通解为.再由初始条件,得又由().1aybx dx dy =可得其解为 ()3 ,202022 bx ay k k bx ay -==-而 (1) ()().231000202011y a b y a bx ay akt y t x =-=-===时,当 即乙方取胜时的剩余兵力数为.230y 又令().0222,01100001=-⎪⎭⎫⎝⎛++⎪⎭⎫⎝⎛-=t ab t ab e y x e y x t x )得由(注意到000020022,1x y y x e y x t ab -+==得. .43ln ,3121bt e t ab =∴=∴ (2) 若甲方在战斗开始后有后备部队以不变的速率r 增援.则().,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 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 rNh N x --=, 2412N rNh N x -+=易知:21N x <, 22Nx > ,0)(1'>x F ,0)(2'<x F ∴平衡点1x 不稳定,平衡点2x 稳定.(2)即 )1(max Nxrx h -=,易得 2*0N x = 此时 4rN h =,但2*0Nx =2N.2.与Logistic 模型不同的另一种描述种群增长规律的是Gompertz 模型:()xNrx t x ln'=.其中r 和N 的意义与Logistic 模型相同. 设渔场鱼量的自然增长服从这个模型,且单位时间捕捞量为Ex h =.讨论渔场鱼量的平衡点及其稳定性,求最大持续产量m h 及获得最大产量的捕捞强度m E 和渔场鱼量水平*0x .解:()t x 变化规律的数学模型为 记 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 不稳定.由前面的结果可得 h =r Er E e 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)无实根,此时无平衡点; ② 当0=∆时,(1)有两个相等的实根,平衡点为20Nx =. Nrxr 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, 22Nx ∴0)('1 x f , 0)('2 x f ∴平衡点1x 不稳定 ,平衡点2x 稳定.20.最大持续产量的数学模型为: ⎩⎨⎧=0)(..max x f t s h即 )1(max N xrx 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)由题设条件可得需求函数、供应函数分别为: 在),(000y x P 点附近用直线来近似曲线h f ,,得到由(2)得 )3( )(0102 y y x x k k -=-++β (1)代入(3)得 )2(0102x x x x x kk k -+-=-++αβ 对应齐次方程的特征方程为 02 2=++αβαβλλ 特征根为48)(22,1αβαβαβλ-±-=当8≥αβ时,则有特征根在单位圆外,设8<αβ,则 即平衡稳定的条件为2 <αβ与207P 的结果一致.(2)此时需求函数、供应函数在),(000y x P 处附近的直线近似表达式分别为: 由(5)得,)( ) y y y β(y )x (x k k k 62010203 -+-=-+++ 将(4)代入(6),得对应齐次方程的特征方程为(7) 024 23 =+++αβαβλαβλλ 代数方程(7)无正实根,且42 ,αβαβ---, αβ不是(7)的根.设(7)的三个非零根分别为321,,λλλ,则 对(7)作变换:,12αβμλ-= 则其中 )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)对应的齐次差分方程的特征方程: 容易算出其特征根为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)对应的齐次差分方程的特征方程: 容易算出其特征根为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 的任一列向量()T nk 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→→→→→.所以此竞赛图是双向连通的.等都是完全路径.此竞赛图的邻接矩阵为令()T e 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=λ,()T S 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节的传送带效率公式,第一种办法的效率公式为 当mn2较小,1 n 时,有E D -=1 , mn E 4≈ ② 下面推导第二种办法的传送带效率公式:对于m 个位置,每个位置放置的两只钩子称为一个钩对,考虑一个周期内通过的m 个钩对.任一只钩对被一名工人接触到的概率是m1; 任一只钩对不被一名工人接触到的概率是m11-;记mq m p 11,1-==.由工人生产的独立性及事件的互不相容性.得,任一钩对为空的概率为n q ,其空钩的数为m 2;任一钩对上只挂上1件产品的概率为1-n npq ,其空钩数为m .所以一个周期内通过的m 2个钩子中,空钩的平均数为 于是带走产品的平均数是 ()122-+-n n npq q m m , 未带走产品的平均数是 ()()122-+--n n npq q m m n ) ∴此时传送带效率公式为 ③ 近似效率公式:由于 ()()()321621121111m n n n m n n m n m n----+-≈⎪⎭⎫ ⎝⎛- 当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 百份纸,则收益函数为 收益的期望值为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∈∀,有()()0=θθg f ,且()()()()0,0,00,00==ππf g f g ,则[]πθ2,00∈∃,使()()000==θθg f .模型求解:令)()()(θθθg f h -= .就有,0)0( h 0)(0)()()( ππππg g f h -=-=.再由()()θθg f ,的连续性,得到()θh 是一个连续函数. 从而()θh 是[]π,0上的连续函数.由连续函数的介值定理:()πθ,00∈∃,使()00=θh .即()πθ,00∈∃,使()()000=-θθg f .又因为[]πθ2,0∈∀,有()()0=θθg f .故()()000==θθg f .9. (1)某甲早8:00从山下旅店出发,沿一条路径上山,下午5:00到达山顶并留宿.次日早8:00沿同一路径下山,下午5:00回到旅店.某乙说,甲必在两天中的同一时刻经过路径中的同一地点.为什么?(2)37支球队进行冠军争夺赛,每轮比赛中出场的每两支球队中的胜者及轮空者进入下一轮,直至比赛结束.问共需进行多少场比赛,共需进行多少轮比赛.如果是n 支球队比赛呢?解:(1)方法一:以时间t 为横坐标,以沿上山路径从山下旅店到山顶的行程x 为纵坐标,第一天的行程)(t x 可用曲线(I )表示 ,第二天的行程)(t x 可用曲线(I I )表示,。
实验一动力系统一、实验目的与要求掌握运用软件求解动态系统模型,通过研究散点图得到动态系统的内在性质和长期趋势。
通过对数据进行处理,归纳出动态系统模型。
1、用Excel对数据进行处理,建立动态系统模型并且进行验证;2、用Excel画散点图,对动态系统模型解的长期趋势进行分析;3、用Excel求解动态系统模型并估计均衡点;4、用Excel分析多元动态系统模型。
二、实验内容Example 1.1 P9 研究课题第一题随着汽油价格的上涨,今年你希望买一辆新的(混合动力)汽车。
你把选择范围缩小到以下几种车型:2007Toyota Camry混合动力汽车2007Saturn混合动力汽车2007Honda Civic混合动力汽车2007Nissan Altima 混合动力汽车2007Mercury Mariner混合动力汽车。
每年公司都向你提供如下的“优惠价”。
你有能力支付多达60个月的大约500美元的月还款。
采用动力系统的方法来确定你可以买那种新的混合动力系统汽车。
混合动力汽车“优惠价”(美元)预付款(美元)利率和贷款持续时间Saturn 22045 1000 年利率5.95%,60个月Honda Civic24350 1500年利率5.5%,60个月Toyota Camry26200 750年利率6.25%%,60个月Mariner27515 1500年利率6%%,60个月Altima24900 1000年利率5.9%%,60个月解答如下,对五家公司分别建立动力系统模型:Saturn:Δb n=b n+1-b n=0.0595b n-6000b n+1= b n+0.0595b n-6000b0=21045Honda Civic:Δb n=b n+1-b n=0.055b n-6000b n+1= b n+0.055b n-6000b0=22850Toyota Camry: Δb n=b n+1-b n=0.0625b n-6000b n+1= b n+0.0625b n-6000b0=25450Mariner:Δb n=b n+1-b n=0.06b n-6000b n+1= b n+0.06b n-6000b0=26015Altima: Δb n =b n+1-b n =0.059b n -6000b n+1= b n +0.059b n -6000 b 0=23900Excel 操作步骤:1.打开excel 表格,输入如下表格::2.用智能标识把月份从0拉到5:3.在B 5 输入= B 4+0.0595B 4-6000,回车后下拉即可可到序列B=(16297.18, 11266.86, 5937.238,…)。
同理在D,F,H,J 行输入,得到如下表格:4. 在插入→图表→XY 散点图,选中数据格就可得出下表: (1)选中A1到B9的数据,建立散点图,得到Saturn 表:Saturn-10000-500005000100001500020000250000123456月份余额Saturn 余额(2)选中C1到D9的数据,建立散点图,得到Honda Civic 表:Honda Civic-500005000100001500020000250000123456月份余额Honda Civic 余额(3)选中E1到F9的数据,建立散点图,得到Toyota Camry 表Toyota Camry0500010000150002000025000300000123456月份余额Toyota Camry 余额(4)选中G1到H9的数据,建立散点图,得到Mariner 表Mariner0500010000150002000025000300000123456月份余额Mariner 余额(5)选中I1到J9的数据,建立散点图,得到Altima 表Altima-50000500010000150002000025000300000123456月份余额Altima 余额由图可知:Saturn 表的线最早与X 轴相交,故我们可以得出应当购买Saturn 公司的汽车。
Example 1.2 P16 习题第二题下列数据表示从1790到2000年的美国人口数据 Yearpopulation Year Population Year Population 1790 3,929,000 1870 38,558,000 1940 131,669,00018005,308,000 188050,156,000 1950150,697,0001810 7,240,000 1890 62,948,000 1960 179,323,000 1820 9,638,000 190075,995,000 1970203,212,000 1830 12,866,0001910 91,972,0001980226,505,000 184017,069,000 1920105,711,0001999248,710,000185023,192,0001930122,755,000 2000281,416,000186031,443,000求出能够相当好地拟合该数据的动力模型,通过画出模型的预测值和数据值来测试你的模型。
解答如下:首先均差计算公式可得下列差分表divided difference table均差Year Observed population∆2∆3∆4∆1790392,9001800530,800 13,79018107,240,000 670,920 32856.518209,638,000 239,800 -21556-1813.75183012,866,000 322,800 4150856.866766.76542 184017,069,000 420,300 487524.16667-20.8175 185023,192,000 612,300 9600157.5 3.333333 186031,443,000 825,100 1064034.66667-3.07083 187038,558,000 711,500 -5680-544-14.4667 188050,156,000 1,159,800 22415936.537.0125 189062,948,000 1,279,200 5970-548.167-37.1167 190075,995,000 1,304,700 1275-156.59.791667 191091,972,000 1,597,700 14650445.833315.05833 1920105,711,000 1,373,900 -11190-861.333-32.6792 1930122,755,000 1,704,400 16525923.833344.62917 1940131,669,000 891,400 -40650-1905.83-70.7417 1950150,697,000 1,902,800 505703040.667123.6625 1960179,323,000 2,862,600 47990-86-78.1667 1970203,212,000 2,388,900 -23685-2389.17-57.5792 1980226,505,000 2,329,300 -2980690.166776.98333 1990248,709,873 2,220,487 -5440.64-82.0212-19.3047 2000281,416,000 3,270,613 52506.271931.56450.33962根据excel中“工具→数据分析→回归”,可得如下图像50,000,000100,000,000150,000,000200,000,000250,000,000300,000,0000510152025系列1多项式 (系列1)模型:y = 670127x 2 - 3E+06x + 8E+06Example 1.4 P50 第四题假定斑点猫头鹰的主要食物来源是单一的食饵:老鼠。
生态学家希望预测在一个鸟兽类保护区里斑点猫头鹰和老鼠的种群量水平。
令 M n 表示 n 年后老鼠的种群量, On 表示 n 年后斑点猫头鹰的种群量。
生态学家提出了下列模型: Mn+1=1.2Mn-0.01OnMn On+1 =0.7On+0.002OnMn生态学家想知道在栖息地两个种群能否共存以及结果是否对起始种群量敏感。
(a) 比较上面模型中系数的正负号和例 3 中猫头鹰-模型中系数的正负号。
依次解释正在建模的捕食者——食饵关系中四个系数 1.2、-0.01、0.7 和 0.002 的正负号的意义。
(b) 对下列表中初始种群数量进行检验并预测其长期行为:猫头鹰 老鼠 猫头鹰 老鼠 情况 A 150 200 情况 C 100 200 情况 B 150150300情况 D1020(c) 现在利用给定的起始值对不同的系数的值做实验,然后再试不同的起始值。
长期行为是怎样的你的实验结果是否表明模型对系数是敏感的是否对起始值敏感?解答如下:(a)1.2和0.7分别是老鼠和猫头鹰增长率,都是正常数。
猫头鹰的存在是为了降低老鼠的增长率,反之亦然。
OnMn为两种生物竞争的激烈程度。
-0.001的负号表示随着竞争激烈程度的增加,老鼠的数目不断减少。
0.002的正号表示随着竞争激烈程度的增加,猫头鹰的数目不断增加。
(b)平衡点:如果把(M,O)成为平衡点,那么必须同时有M=Mn+1=Mn和O= On+1= On,把它们带入模型给出0=M*(0.2-0.001*O)0=0*(-0.3+0.002*M)平衡点的意义:第一个方程表明如果M=0或O=200,那么老鼠的种群量没有变化。
第二个方程表明如果O=0或M=150,那么斑点猫头鹰的种群量没有变化。
如下图(1)所示在(M,O)=(0,0)和(M,O)=(150,200)处于平衡点,因为两个种群的种群量在这两个点都没有变化。
Excel操作步骤:1.打开excel表格,输入如下表格:2.用智能标识把天数从2拉到30:3.在B4输入=0.7*B3+0.002*B3*C3回车后下拉即可可到序B=(200,200,200…)。
在C4输入=1.2*C3-0.001*B3*C3回车后下拉即可可到序B=(150,150,150…)。
得到如下表格:4.在插入→图表→XY散点图,选中数据格就可得出下表:选中A1到C32的数据,建立散点图,得到平衡表:平衡0501001502002505101520253035天数生物的数目猫头鹰老鼠(图1)图1:如果老鼠的种群量从150开始而猫头鹰的种群量从200开始,那么这两个种群都停留在它们的起始值处。