数学的实践与认识 第26卷(1996)第1期 1995年全国大学生数学建模竞赛 姜启源 (清华大学,北京100081) 由国家教委高教司和中国工业与应用数学学会共同主办的“1995年全国大学生数学 建模竞赛”于1995年9月27日至29日举行。来自25个省(市、自治区)、259所院校的 1234个队参加了这次竞赛。 尧赛答卷首先在22个赛区进行初评,评出各赛区的获奖者。然后各赛区按一定比例 将优秀答卷送全国组委会,全国组委会聘请专家从154份答卷中共评出全国一等奖35 名二等奖75名,11月21日在北京举行了颁奖仪式。 全国大学生数学建模竞赛是1992年开始由中国工业与应用数学学会举办的。国家教 委对这项活动十分重视决定自1994年起由教委高教司和中国工业与应用数学学会共同 主办,每年…次。这两年来参赛院校数和队数都以50%的速度递增。 这项竞赛的题目一般来源于工程技术和管理科学等方面经过适当简化加工的实际问 甌不要求预先掌握深入的专门知识,而具有较大的灵活性供参赛者发挥创造能力。今年 的两道题是分别由复旦大学谭永基、俞文鱿和浙江大学刘祥官、李吉鸾提供的。为了更广 泛、有效地收集适合这项竞赛的题目和素材,再次向全社会诚征赛题,通讯地址:北京清华 大学应用数学系,数模竞赛全国组委会办公室,邮编100084 为了与广大参赛同学进行交流,并对今后的竞赛予以适当引导,全国评阅委员会选择 了14篇优秀答卷全文或部分发表,并请两位命题者(也是评阅人)撰文讲评。 这里发表的论文都只是学生们在三天内写出的,为了保持原貌只作了文字上的个别 修正和繁琐处的删节,文章不可避免地存在着相当多的不妥之处,请读者谅解。 下面是本次竞赛的题目、获奖名单、各赛区参赛院校数和队数、评阅委员会名单 题 个飞行管理问题 在约10,000米高空的某边长160公里的正方形区域内,经常有若干架飞机作水平飞 行、区域内毎架飞机的位置和速度向量均由计算机记录其数据,以便进行飞行管理。当 架欲进入该区域的飞机到达区域边缘时,记录其数据后,要立即计算并判断是否会与区域 内的飞机发生碰撞。如果会碰撞,则应计算如何调整各架(包括新进入的)飞机飞行的方向 角,以避免碰撞。现假定条件如下: 1)不碰撞的标准为任意两架飞机的距离大于8公里 2)飞机飞行方向角调整的幅度不应超过30度 3)所有飞机飞行速度均为每小时800公里 4)进入该区域的飞机在到达区域边缘时,与区域内飞机的距离应在60公里以上; 2 01995-2004 Tsinghua Tongfang Optical Disc Co, Ltd. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
数学的实践与认 5)最多需考虑6架飞机 6)不必考虑飞机离开此区域后的状况。 请你对这个避免碰撞的飞行管理问题建立数掌模型.列出计算步骤.对以下数据进 计算(方向角误差不超过0.01度),要求飞机飞行方向角调整的幅度尽量小 设该区域4个顶点的座标为(0,0),(160,0),(160,160),(0.160)。 记录数据为: 飞机编号 横座标X纵座标Y方向角(度) 12345 155 220 145 159 150 230 新进入 52 注:方向角指飞行方向与X轴正向的夹角。 试根据实际应用背景对你的模型进行评价与推广。 B题天车与冶炼炉的作业调度 某钢铁厂冶炼车间的厂房布局是,地面沿一直线依次安置着7个工作点:辅料供应处 P;A组3座转炉(冶炼成品钢)A1,A2,A3;B组2座冶炼炉(冶炼半成品钢,简称半钢 B1,B2;原料供应处Q。这些设备的上方贯通着一条运送物料的天车轨道,上面布置着若 干天车T1.T2…Tn为炉子作业服务。布局示意图如下 天车与冶炼炉的作业过程与工序为:大车从Q处吊起原料一罐(吊罐时间、)运至BI 或B2处放下(放罐时间)·并将上一妒的原料空罐吊起(吊空罐时间κ)返回Q处放下 (放空罐时间)。B组炉在原料罐放下后即可在辅助作业下开始冶炼(冶炼时间μ),冶炼 后将半钢倒入空半钢罐(时间计入t),由天车吊起半钢罐(吊罐吋间t)运至A1或A2 A3处将半钢倒入转炉(倒入时间t),并将空罐返回B1或B2处放下(放空罐时间t)。再 由天车从P处吊起辅料一槽(吊槽时间t)运至A1或A2、A3处加入转炉(加入时间t,) 并将空槽返回P处放下(放空槽时间t)。A组炉在半钢和辅料加入后即可开始冶炼(冶炼 时间t),冶炼后成品钢的输出不用天车(输出时间计入t)。天车通过相邻两个工作点的 行时间都相同,记为t.。 由于各台天车在同一条轨道上运行.因此其顺序信置T1.T2,…Tn不可交换。在同 时间同一座炉子上只能允诈一台天车作业;但P.Q两处可以允许多台车同时作 在P,A1,…Q每两个相邻工作点之间最多能容纳2台天车同时停放 2 01995-2004 Tsinghua Tongfang Optical Disc Co, Ltd. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
第26卷(1996)第1期 天车与冶炼炉作业调度的要求为:(1)成品钢产量尽量高:(2)各台天车的作业率(天 车作业时间所占比例)尽量均衡(考虑到设备及人员安全等因素,一般天车作业率不超过 η0%);(3)绝不允许出现天车相撞等事故;(4)调度规则尽量简明,以利于现场人员使用 现设定:t=48,t6=27,t=3,t=2,t.=2,t=3,t=5,t=2,tx=2, t=1,t=3,t=2(单位:分钟),t=15秒;A组炉平均每炉产量W=120吨。在不 超过5台天车的条件下,设计一种满足上述要求的天车与冶炼炉的作业调度方案: (1)各台天车负责哪些作业(列出《工序清单》; (2)在所给方案的一个运行周期内,每一时刻天车和冶炼炉处于什么状态(画出《天 车一炉子作业运行图》); (3)一份供现场人员使用的《调度规则说明书》; (4)在所给方案下计算各台天车的作业率 并按每天冶炼炉数估计该车间成品钢的年产量(扣除设备维修日,每台转炉作业日每 年按300天计算) 实际生产过程中,t,t,…,t都是随机的(上面设定的数值可视为平均值),讨论你的 调度方案如何适用于实际生产过程。试提出该车间提高钢产量到年产300万吨以上的建 “1995年全国大学生数学建模竞赛”获奖名单 一等奖:35名(排名不分先后) 学校 学生 指导教师 南开大学 凌晖熊德华杨杰叶剑平 杭州电子工业学院 熊宏伟张远福厉莹数模组 浙江大学 赵明洁王昆龚明数模组 河北师范大学 刘金生赵谦孙淑英指导教师组 南昌大学 彭小华薛峰肖隆陈涛 中山大学 任远陈海波何亚彬张磊 四川轻化工学院 邱玉平谭小术干斌武亦文 四川联合大学(成都科大)倪敬能姚文俊高鹏杨志和 重庆工业管理学院 杨银芳郭安蒋鹏宋江敏 重庆工业管理学院 朱长国黄金曦叶显锋苏宏 重庆邮电学院 李莉莎阙劲峰何小玉杨春德 哈尔滨工业大学 张熙李浒刘世霞时培林 山东大学 刘铁成张良聂兆虎许宝刚 2 01995-2004 Tsing hua Tongfang Optical Disc Co, Ltd. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
数学的实践与认识 华中理工大学 孙黎明郫跷玲房靖何南忠 武汉水电大学 庞旭曹志芳下渺林石岗 东北师范大学 徐文兵崔郁青杨光白玉山 中南工业大学 蒋超张杰王日中韩旭里 中南林学院 徐元军曾九林韩伟群潘冬光 兰州铁道学院 何新宇贡力平庞晓林张建勋 湘潭大学 陈靖周素华黄秋波成央金 云南曲靖师专 吴玉峰丁雪梅吴元勇陈世联 北京航空航天大学 杜序袁灯山杨黎明赵杰民 北京农业工程大学 王国卿李加福毕诚刘军风 北京大学 王崧于劲松陆昱雷功炎 清华大学 刘学胡晨陈涵高策理 上海大学 李刚尹民傅晓陈达段 复旦大学 吴伟标王立峰嵇元曹沅 复旦大学 俞寅喆朱丹宇俞希晨廖有为 复旦大学 谭浩南朱正元刘剑蔡志 华东理工大学 李汉涛张玉玺·段立松陆元洪 广西大学 王烨韦世蒙李勇潘涛 空军气象学院(南京) 宋晓亮王东滕加俊 安徽机电学院 萤小虎巩禧云王庚 中国科学技术大学 程谟嵩罗亚俞天越冯宇 中国科学技术大学 黄春峰饶红玲刘伟于清娟 二等奖:75名(排名不分先后) 学校 学生 指导教师 南开大学 张心正朱玉鹏徐晓轩黄五群 南开大学 刘洪杰曾维微杜华坤王厦生 南开大学 冯少新陈戍向军黄五群 天津大学 杨立宇陈刚黄目亮数学建模教研小组 天津理工学院 姜兆明张杰王方陈东升 杭州电子工业学院 陈建华陈寒张建樑数模组 河北机电学院 陈云生李树民王燕青王容 河北机电学院 刘志会高树红赵霞张隽礼 河北大学 杨晓晖吴坤玲李念龙指导教师组 景德镇陶瓷学院 骆双坚高正洪成志峰周永正 江西农业大学 李卫东李军辉林新春欧阳兴 华南师范大学 李丽间杨戈锋张淑华曹汝成 华南理工大学 欧永斌张上弋陈薇谢乐军 华南理工大学 高唰李翔杜充傅红卓 2 01995-2004 T hua Tongfang Optical Disc Co, Ltd. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
绾26讼()第1期 电子科技大 蒋海波何莉李恩杨晋浩 T:了院 胡家望兰建军王光荣刘同楷 西南父通大 唐建华郭军华E福胜袁俭 轮化1院 向邦云郑衡王.超冯家竹 川川联合大学(四川大学)罗谦胡朝波余钢程中瑷 重庆镧铁々科字校 向毅卓国锋刘洪雷鸣 重庆大学 张洪伟陈众唐晓苏刘琼荪 重庆大学 向志海陈冠饶李玉刚龚劬 重庆通信学院 樊景渝肖清伟章立李元红 解放乍信息程学院 吕声马智高丰韩中庚 解放军信息I程学院 黄秋生张亚娟蒋东毅韩中庚 郑州吭空1业管理学院 邵华钢肖冬董俊刘道远 哈尔滨1程大学 何平刘希斌程婧容张晓威 哈尔滨科技大学 李希斌冯艳军费优松陈东彦 黑龙汇商学院 贺晓明李灵活毛小勇吴刚 山东工业大学 王怀磊王向军徐磊孙 山东业大学 来翔高洪峰魏强李保健 曲阜师范大学 宁如云陈茂银莫修明冯成进 华中理上大学 杜劲松胡伟湘马红波齐欢 武汉业大学 陈世荣郭鹏李健荣王祖喜 太原重机学院 李正文陆介伦杨京波冯巍 华北工学院 曹阳黄伟军陆海平吴强 东北大学 顾晓伟吴军华巴力颖赵鸿金 东北大学 马正品李校兵陈建兵黄卫祖 大连理工大学 霍明于丕强牛大田赵立中 东北电力学院 罗兰黄嘉升杨文龙田知能 吉林工业大学 曲鹏李炳辉余朝蓬张魁元 吉林工学院 王炳文白海石王兆升乌成伟 吉林大学 胡锡俊刘冷宁王宏宇吕显瑞 兰州铁道学院 赵京才胡建新冯德泉栗永安 兰州铁道学院 郑秋宁徐昌山黄景春张建勋 兰州铁道学院 马东升马学锋龙维洋吕新忠 兰州大学 刘铁荣范晓军姚海元李效虎 西安电子科技大学 李景峰苏涛赵国栋马建锋 国防科技大学 刘伶训董威谭郁松吴孟达 国防科技大学 陈明刘雅浪姚崎吴孟达 国防科技大学 王辰李中升武洁覃左平 国防科技大学 高曰超李冬冬王银华吴孟达 国防科技大学 赫新夏刚李爱平吴翊 2 O1995-2004 Tsinghua Tongfang Optical Disc Co, LId. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
数学的实践与认识 云南大学 杜强张晶谢洪波教师指导组 云南大学 肖明海冯俊田雯教ⅷ指导组 昆明师专 李红芳郭文俊李刚教师指导组 云南师范大学 孙兴平黄鹏施宏昌教师指导组 北京航空航天大学 邝富华冉晓林霍继文李卫国 北京商学院 王建军李瑾白虹黄先开 北京邮电大学 香盈波李云立林胜丁金扣 北京大学 张霖涛罗卫东丁立吴崇试 北京大学 林涛黎德元凌海滨孙山泽 清华大学 冯汉鹰诸葛丰杨小苇包维柱 清华大学 刘军宁肖狀谢峰李栓虎 上海交通大学 周吴芦烈杨荣震张建强 华东理工大学 李希明李琦王奇许三保 桂林电子工业学院 常志泉吴岭刘召卫周孝华 东南大学 丁剑张德冯南姚瑞波 东南大学 关永涛杨杉李晟孙志忠 南京航空航天大学 荀海波王冬夏顺东顾玉娣 南京航空航天大学 刘贤军金晓虎顾正晖张毅 合肥工业大学 朱明张啸郭志军杜雪樵 中国科学技术大学 孙亮贾英东张珲周智 中国科学技术大学 许锦波贾志峰朱朝阳陈发来 福州大学 陈桂阳林霖游香明林可容 1995年全国大学生数学建模竞赛各赛区 参赛院校数及队数 序号 赛区 校数 队数 江苏安徽 浙江 123456789 江西 山东 湖南 湖北 广西 四川 171 2 O1995-2004 Tsinghua Tongfang Optical Disc Co, LId. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
第26卷(199G6)第1期 6 23 甘肃 广东 天津 河北 山西 辽宁 10 吉林 黑龙江 上海 北京 联合 合计 联合赛区包括福建、内蒙古和青海 1995年全国大学生数学建模竞赛评阅委员会名单 (排名不分先后) A题小组负责人:王强北京应用物理与计算数学研究所 成员 王荫清四川联合大学(西区)应用数学系 李煤山东工业大学数理系 胡小东中科院应用数学所运筹室 雷功炎北京大学科学与工程计算系 谭永基上海复旦大学数学系 B题小组负责人:项可凤中国科学院系统科学所 成员 朱道元南京东南大学数力系 叶其孝北京理工大学应用数学系 尹景学吉林大学数学系 刘祥官浙江大学应用数学系 林元烈清华大学应用数学系 石坚中国科学院系统科学所 2 01995-2004 Tsinghua Tongfang Optical Disc Co, Ltd. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
数学的实践与认识 飞行管理问题答卷评述 谭永基 (复旦大学,上海200433) 本题是以空域飞行管理为背景,经简化和整理而成的一个赛题。该问题主要可以归结 为非线性规划模型或经一定简化,建立线性规划模型。由于实际的需要,提出的箅法应在 计算机上快速地实现 非线性规划模型及求解 设六架飞机在调整时的方向角为0,调整后的方向角为0=0,+△0(=1,2,…,6), 设任意两架飞机在区域内的最短距离为d(θ,0,、),那么问题的非线性规划模型为 使得 (+△0.,0,+△0)>81≤i,j≤6,≠ 绝大多数答卷能正确逮立模型,有的答卷在建模时,出于某些考虑加强了不碰撞的要 求,如要求在调整后的0.2√2小时內不发生碰撞或永远不允许发生碰撞,从而简化了 d的表述。 求解此模型的一种方法是枚举法,但是枚举工作量极大,必须釆取逐步求精(细分) 隐式枚举、枚举和二分法相结合等技巧,方能在较短时间内求得符合精度的最优调整方 案。参赛答卷中釆用了许多提高枚举效率的措施。有的答卷在枚举时采用了 Monte-Car lo法,随机产生大量△B}组合,从其中的可行解中选出最优解。这种方法可显著提高计 算速度,有一定新意。 另一种解法是引进惩罚函数,将问题化为无约束极值,然后将其极小化。在答卷中采 用了丰富多采的无约束化和极小化技术。有的方法,如能量梯度化有一定的特色。然而 非线性规划方法能否快速找到全局最优解,十分强烈地依赖于初值的选取:有的试卷用步 长较粗的枚举得到的解作为初值然后再用极小化的方法得到最优解,计算速度很快.可 以实时调整,达到∫实用的要求 2 01995-2004 Tsinghua Tongfang Optical Disc Co, Ltd. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
第26卷(1996)第1期 有不少答卷用∑△a|2,max1△|作为优化的目标,这都是合理的。 、线性规划模型 若不考虑区域的范围,用相对运动的观点不难得出两架飞机调整方向角后不发生碰 撞,调整角度应满足的条件。为简单起见,排除某些特殊的情形,条件是: ①若0,<(2+支)m36360-a 其中,0,D,∈[0,360°),a1= arcsin(3),,为调整时两飞机的距离,角度均以自 架飞机位置指向第j架飞机的矢量为始边计算 若以∑|△|,max|△a|作为目标函数,以上述不碰撞条件和|△a≤30作为约束 条件,是一个可归化为多个线性规划问题的规划问题。求解这一系列线性规划,比较它们 目标函数值就能得到原问题的最优解 许多答卷用类似的方法将问题归结为线性规划求解。但有不少答卷在建立不碰撞条 件时疏忽了一些情形从而条件不够完整。对某些飞机初始位置,可能产生错误 由于变量较少,线性规划求解速度较快,这是采用这一模型的优点。然而,若必须排除 在区域外的碰撞,由于条件不是线性的,归结为线性规划问题就有一定困难。有的答卷提 出了一些克服这一困难的方法。较有成效的方法有两种。一种是在△不大时将条件展 开略去高阶项,另一种用低维枚举确定△θ+△θ,的允许范围,再用线性规划求解 2 O1995-2004 Tsinghua Tongfang Optical Disc Co, LId. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved
数学的实践与认识 飞行管理问题的实时算法 谭浩南朱正光刘剑 (复旦大学,上海200433 指导教师:蔡志杰 编者按:该文以区域内飞机谲整飞行角幅度的平方和为目标数,以自调整时刻起ω.3·小· 吋內飞机两两不发生碰撞为约束祭件,建立了一个非线性规划模型,用逐步求精的直接搜索和 引入罚函数化为无约東极值用序列无约来最小化(SUMT)两种方法进行求解。 作者在建模和求解时从实际霄要出发,精益求精,将上述两种方法桕结合,得到了精度高 基本上是实时的方法与程序,还对模型与方法作出了恰当的分析与评价,文章清腑、完鳘。 摘要本文讨论了在一定区域空间内进行飞行管理避免飞机相撞的模型,提出了直 搜索法和非线性规划(SUMT)法两种解法,并将两种方法有机结合,得出的算法在486做机上 计算时间小于10秒,误差不超过0.01度,完全符合问题的要求 本文接着给出四种不同情况分别用两种方法求解,进行比较检验,取得很好的吻合,充分 说明了模型3的可靠性。本文还对模型的误差进行分析并对模型进行推广 关键词:非线性规划,直接搜索,罚函数 问题的提出(略) 二、问题的分析 该问题是一个在一定约束条件下的最优化问题,初步分析题意后可知约束条件是非 线性的,难以化归为线性规划问题。由于题目涉及数据变量不是太多,可以考虑用逐步求 精的直接搜索法求解。由于题目要求的精度较高,而对于计算时间的要求也较高,如果求 解时间在2、3分钟以上将失去任何实际意义。我们将求解时间上限定为0.5分,以符合实 际的要求。直接搜索法求的近似解难以同时满足两方面的要求。但直接搜索法至少能在 较短的时间内得到一个较好的可行解,这就为运用非线性规划的方法提供了条件。非线性 规划的算法种类繁多,但均只适用于某些类型的问题。由于缺乏适用的计算机软件包,我 们自行编写了实现算法的程序。综合程序准备时间和收敛速度两方面因素我们选择了 sUMT算法。SUTM算法与直接搜索法相结合,使我们能够在足够短的时间内找到问题 的足够精确的解。 2 01995-2004 Tsinghua Tongfang Optical Disc Co, Ltd. All rights reserved
© 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved