数学建模竞赛讲座(清华大学姜启源)

综合题目参考答案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元，总计元总计*****元。

第二届全国数学建模微课程（案例）教学竞赛获奖结果“第二届全国数学建模微课程（案例）教学竞赛”（以下简称建模微课竞赛）是由中国工业与应用数学学会数学模型专业委员会和全国大学生数学建模竞赛组委会联合主办，高等教育出版社协办。

自2017年12月通知发出之后，得到了全国高校相关教师们的积极响应和参与。

本次建模微课竞赛分两个阶段，第一阶段为初赛，对所有参赛作品的文稿、PPT课件和教学视频等材料由同行专家进行了认真的评审，在众多的参赛作品中评选出26项有代表性的作品进入第二阶段的现场授课决赛。

第二阶段的决赛于2018年7月在哈尔滨举行的“2018年全国数学建模培训与应用研究研讨会”期间进行，采用现场授课的方式，每个参赛教师就参赛作品进行时长不超过30分钟的现场授课，要求从案例的背景、问题提法、模型的建立与求解过程和结果的分析等内容进行授课。

决赛的现场授课分为两个场地同时进行，邀请了16名国内的同行专家评委现场对所有作品的授课效果进行评判打分，并有来自全国各高校的400余名同行教师聆听了决赛授课。

决赛主要依据案例的原创性、新颖性和应用价值、授课内容的组织设计和授课效果等方面进行评价。

经过一整天紧张激烈的竞赛，最后汇总专家评委的评分和参赛教师的相互评分，评选确定出一等奖4项，二等奖8项，三等奖14项（获奖作品名单见附件）。

本次建模微课竞赛是第二次举办这类活动，获得了很好的效果，并受到了广大数学建模教师的欢迎。

参赛作品所涉及的案例多数取材于实际科研和现实生活的原创性案例，也包括一些经典案例的推广与应用，为提高各高校的数学建模和数学实验课程的教学起到了促进作用，也为后续数学建模案例的进一步开发利用提供了优秀素材。

本次竞赛活动得到了高等教育出版社的大力支持和资助，为竞赛活动的顺利进行和成功举办提供了有力保障。

对于获奖作品的进一步开发和推广使用，我们将会同高等教育出版社有关部门和参赛者进行协商，将优秀的作品尽快与广大教师见面。

数学模型课后答案姜启源【篇一：姜启源《数模》习题选解】方案模型构成：以阈值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个共顶点。

关于数学建模竞赛的一点思考、总结和建议关于数学建模竞赛的一点思考、总结和建议宋一凡环境保护与安全工程学院核安全工程专业大学生活即将结束，回顾几年的经历，数学建模竞赛留给我太多的回忆。

虽然数模竞赛已经远去，但至今看到听到“三天三夜72小时”时，精神还会为之一振。

在要告别数模竞赛的时候，想写一点自己零零碎碎的思考和总结，并给以后参赛的学弟学妹一点建议。

1. 关于我的数模之路大一从学长口中知道了数模竞赛，就想参加，自学了姜启源的《数学模型》，但校赛时，队友不给力使第一次校赛不了了之，至今仍然遗憾大一时校赛未能入围；大二时，和本院的两个同学组队，比我高一级的闯哥给了不少经验和资料，经过暑假的培训和多次模拟赛训练，12年国赛拿到了湖南赛区的三等奖。

13年寒假，留在学校参加美赛，偌大的宿舍楼空无一人，好不凄凉，南方湿冷的冬天让我这个北方人冻得难以忍受，搞完比赛回到家时已经是腊月二十七夜里，美赛S奖使我很失落，也从中找到了自己的很多不足之处。

因今年考研，本不愿参加国赛，但两位新队友的盛情邀请让我不忍拒绝，于是重新组队，再战国赛，一雪前耻，最后拿到国家一等奖，为大学的数模之路画上一个圆满的句号。

从大一到现在，关于数模的比赛，热身赛、校赛、模拟赛、国赛、美赛，大大小小不记得参加过多少次，也不知道熬过了多少个“72小时”。

建模、程序员、写手，三个角色的工作我都认认真真做过，饱尝里面的酸甜苦辣，一步一个脚印走来，最后得到一个不错的成绩，收获颇多，感触颇深。

数模给我打开了一扇窗，窗外的世界带给我不一样的精彩，而不仅仅是拿几张证书，加几分综测。

外人看来，数模痛苦、费人，而我感觉数模自由、快乐。

尤其是竞赛结束，早上八点交卷的时刻，经过三天三夜的努力，队友通力合作，从第一天的一筹莫展，到最后一天的顺利解决，疲惫、兴奋、满足、急切、不安，很多的感受一时涌上心头，那是只有真正参加比赛的人才能体会到的快乐！2. 关于数学建模竞赛的作用在做一件事情之前总会去思考做成这件事情有什么好处，这样的心里再正常不过了。