运筹学的优化算法ppt课件.ppt
《运筹学的优化算法ppt课件.ppt》由会员分享,可在线阅读,更多相关《运筹学的优化算法ppt课件.ppt(199页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、1运筹学的优化算法运筹学的优化算法2345693A 非线性交调的频率设计非线性交调的频率设计: 拟合、规划拟合、规划93B 足球队排名次足球队排名次: 矩阵论、图论、层次分析法、矩阵论、图论、层次分析法、整数规划整数规划94A 逢山开路逢山开路: 图论、插值、动态规划图论、插值、动态规划94B 锁具装箱问题锁具装箱问题: 图论、组合数学图论、组合数学95A 飞行管理问题飞行管理问题 : 非线性规划、线性规划非线性规划、线性规划95B 天车与冶炼炉的作业调度天车与冶炼炉的作业调度: 非线性规划、动态非线性规划、动态规划、层次分析法、规划、层次分析法、PETRI方法、图论方法、排队方法、图论方法、
2、排队论方法论方法96A 最优捕鱼策略最优捕鱼策略:微分方程、积分、非线性规划微分方程、积分、非线性规划796B 节水洗衣机节水洗衣机:非线性规划非线性规划97A 零件参数设计零件参数设计:微积分、非线性规划、随机模微积分、非线性规划、随机模拟拟97B 截断切割截断切割:组合优化、几何变换、枚举、蒙特组合优化、几何变换、枚举、蒙特卡罗、递归、最短路卡罗、递归、最短路98A 投资收益与风险投资收益与风险:线性规划、非线性规划线性规划、非线性规划98B 灾情巡视路线灾情巡视路线:最小生成树、最小生成树、Hamilton圈、圈、旅行商问题旅行商问题99A 自动化车床管理自动化车床管理:积分、概率分布、
3、随机模拟、积分、概率分布、随机模拟、分布拟合度检验分布拟合度检验899B 钻井布局钻井布局:几何变换、枚举、最大完全子图、几何变换、枚举、最大完全子图、混合整数规划混合整数规划00A DNA分类分类:神经网络、最小二乘拟合、统计神经网络、最小二乘拟合、统计分类分类00B 管道订购和运输管道订购和运输:最短路、二次规划最短路、二次规划01A 血管的三维重建血管的三维重建:数据挖掘、曲面重建与拟数据挖掘、曲面重建与拟合合01B 公交车调度公交车调度:非线性规划非线性规划02A 车灯光源优化设计车灯光源优化设计:最优化最优化02B 彩票中的数学彩票中的数学:概率与优化概率与优化903A SARS S
4、ARS的传播:的传播: 微分方程、差分方程微分方程、差分方程 03B 露天矿生产的车辆安排:露天矿生产的车辆安排: 整数规划、运输问整数规划、运输问题题04A 奥运会临时超市网点设计:奥运会临时超市网点设计:统计分析、数据统计分析、数据处理、优化处理、优化 04B 电力市场的输电阻塞管理:电力市场的输电阻塞管理:数据拟合、优化数据拟合、优化 05A 长江水质的评价和预测:长江水质的评价和预测: 预测评价、数据预测评价、数据处理处理 05B DVD DVD在线租赁:在线租赁: 随机规划、整数规划随机规划、整数规划 06A 出版社书号问题:出版社书号问题: 整数规划、数据处理、整数规划、数据处理、
5、优化优化1006B Hiv Hiv病毒问题:病毒问题: 线性规划、回归分析线性规划、回归分析07A 人口问题:人口问题: 微分方程、数据处理、优化微分方程、数据处理、优化07B 公交车问题:公交车问题: 多目标规划、动态规划、图多目标规划、动态规划、图 论、论、0-10-1规划规划08A 照相机问题:照相机问题: 非线性方程组、优化非线性方程组、优化08B 大学学费问题:大学学费问题: 数据收集和处理、统计分数据收集和处理、统计分 析、回归分析析、回归分析11 赛题发展的特点:赛题发展的特点:1. 对选手的计算机能力提出了更高的要求:赛题对选手的计算机能力提出了更高的要求:赛题的解决依赖计算机
6、,题目的数据较多,手工计算的解决依赖计算机,题目的数据较多,手工计算不能完成,计算机模拟和以算法形式给出最终结不能完成,计算机模拟和以算法形式给出最终结果。果。 2. 赛题的开放性增大赛题的开放性增大: 解法的多样性,一道赛题解法的多样性,一道赛题可用多种解法。开放性还表现在对模型假设和对可用多种解法。开放性还表现在对模型假设和对数据处理上。数据处理上。 3. 试题向大规模数据处理方向发展试题向大规模数据处理方向发展 4. 求解算法和各类现代算法的融合求解算法和各类现代算法的融合 121.蒙特卡罗方法(Monte-Carlo方法, MC) 该算法又称计算机随机性模拟方法计算机随机性模拟方法,也
7、称统计试验统计试验方法方法。MC方法是一种基于“随机数”的计算方法,能够比较逼真地描述事物的特点及物理实验过程,解决一些数值方法难以解决的问题。 MC方法的雏型可以追溯到十九世纪后期的蒲丰随机投针试验,即著名的蒲丰问题蒲丰问题。 MC方法通过计算机仿真(模拟)解决问题,同时也可以通过模拟来检验自己模型的正确性,是比赛中经常使用的方法。1397年的年的A题题 每个零件都有自己的标定值,也都有自己的容差等级,而求解最优的组合方案将要面对着的是一个极其复杂的公式和108种容差选取方案,根本不可能去求解析解,那如何去找到最优的方案呢?随机性模拟搜索最优方案就是其中的一种方法,在每个零件可行的区间中按照
8、正态分布随机的选取一个标定值和选取一个容差值作为一种方案,然后通过蒙特卡罗算法仿真出大量的方案,从中选取一个最佳的。02年的年的B题题 关于彩票第二问,要求设计一种更好的方案,首先方案的优劣取决于很多复杂的因素,同样不可能刻画出一个模型进行求解,只能靠随机仿真模拟。14 98 年美国赛年美国赛A 题题 生物组织切片的三维插值处理生物组织切片的三维插值处理 94 年年A 题逢山开路题逢山开路 山体海拔高度的插值计算山体海拔高度的插值计算2. 数据拟合、参数估计、插值等数据处理算法 比赛中通常会遇到大量的数据需要处理,而处理数据的关键就在于这些算法,通常使用MATLAB 作为工具。与图形处理有关的
9、问题很多与拟合有关系。 此类问题在MATLAB中有很多函数可以调用,只有熟悉MATLAB,这些方法才能用好。1598年年B 题题 用很多不等式完全可以把问题刻画清楚用很多不等式完全可以把问题刻画清楚3. 规划类问题算法 此类问题主要有线性规划、整数规划、多目标规划、线性规划、整数规划、多目标规划、二次规划二次规划等。竞赛中很多问题都和数学规划有关,可以说不少的模型都可以归结为一组不等式作为约束条件、几个函数表达式作为目标函数的问题,遇到这类问题,求解就是关键了。 因此列举出规划后用Lindo、Lingo 等软件来进行解决比较方便,所以还需要熟悉这两个软件。1698 年年B 题、题、00年年B
10、题、题、95 年锁具装箱年锁具装箱等问题体等问题体现了现了图论问题图论问题的重要性。的重要性。4. 图论问题 这类问题算法有很多,包括:Dijkstra、Floyd、Prim、Bellman-Ford,最大流,二分匹配,最大流,二分匹配等问题。1792 年年B 题用分枝定界法题用分枝定界法97 年年B 题是典型的动态规划问题题是典型的动态规划问题98 年年B 题体现了分治算法题体现了分治算法5. 计算机算法设计中的问题 计算机算法设计包括很多内容:动态规划、回溯搜动态规划、回溯搜索、分治算法、分枝定界索、分治算法、分枝定界等计算机算法. 这方面问题和ACM 程序设计竞赛中的问题类似,可看一下与
11、计算机算法有关的书。1819IP LP IP LP xl* * 判别是否整数解判别是否整数解xI* * = xl* *YesYes去掉非整数域去掉非整数域NoNo多个多个LPLP2021 若所得的最优解的各分量恰好是整数,若所得的最优解的各分量恰好是整数,若所得的最优解的各分量恰好是整数,若所得的最优解的各分量恰好是整数,若所得的最优解的各分量恰好是整数,若所得的最优解的各分量恰好是整数,则这个解也是原整数规划的最优解,计则这个解也是原整数规划的最优解,计则这个解也是原整数规划的最优解,计则这个解也是原整数规划的最优解,计则这个解也是原整数规划的最优解,计则这个解也是原整数规划的最优解,计算结
12、束。算结束。算结束。算结束。算结束。算结束。2223 若松驰问题有最优解,但其各分量不全若松驰问题有最优解,但其各分量不全若松驰问题有最优解,但其各分量不全若松驰问题有最优解,但其各分量不全若松驰问题有最优解,但其各分量不全若松驰问题有最优解,但其各分量不全是整数,则这个解不是原整数规划的最是整数,则这个解不是原整数规划的最是整数,则这个解不是原整数规划的最是整数,则这个解不是原整数规划的最是整数,则这个解不是原整数规划的最是整数,则这个解不是原整数规划的最优解,转下一步。优解,转下一步。优解,转下一步。优解,转下一步。优解,转下一步。优解,转下一步。2425 定界:把满足整数条件各分枝的最优
13、定界:把满足整数条件各分枝的最优定界:把满足整数条件各分枝的最优定界:把满足整数条件各分枝的最优定界:把满足整数条件各分枝的最优定界:把满足整数条件各分枝的最优目标函数值作为上(下)界,用它来判目标函数值作为上(下)界,用它来判目标函数值作为上(下)界,用它来判目标函数值作为上(下)界,用它来判目标函数值作为上(下)界,用它来判目标函数值作为上(下)界,用它来判断分枝是保留还是剪枝。断分枝是保留还是剪枝。断分枝是保留还是剪枝。断分枝是保留还是剪枝。断分枝是保留还是剪枝。断分枝是保留还是剪枝。26270 1 2 3 44 3 2 1x1x22829300 1 2 3 44 3 2 1x1x231
14、3233X1 2X2 2X1 1X2 3P1:(1,9/4)Z=10(3/4)P4: (0,3) Z=9P2:(2,1/2) Z=9(1/2)P3: (1,2) Z=10P:(6/5,21/10) Z=111/1034动态规划应用举例动态规划应用举例35最优化原理:最优策略的任一后部子策略都是最优的。无论以前状态决策如何,从眼下直到最后的诸决策必构成最优子策略。例1最短路线问题AB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143动态规划应用举例动态规划应用举例36例例1 1 多阶段资源分配问题多阶段资源分配问题 设有数量为设有数量为x的某种资源,将它
15、投入两种的某种资源,将它投入两种生产方式生产方式A和和B中:以数量中:以数量y投入生产方式投入生产方式A,剩下的量投入生产方式剩下的量投入生产方式B,则可得到收入,则可得到收入g(y)+h(x-y),其中,其中g(y)和和h(y)是已知函数,并是已知函数,并且且g(0)=h(0)=0;同时假设以;同时假设以y与与x-y分别投入两分别投入两种生产方式种生产方式A,B后可以回收再生产,回收率分后可以回收再生产,回收率分别为别为a与与b。试求进行。试求进行n个阶段后的最大总收入。个阶段后的最大总收入。37若以若以y与与x-y分别投入生产方式分别投入生产方式A与与B,在第一,在第一阶段生产后回收的总资
16、源为阶段生产后回收的总资源为x1=ay+b(x-y),再将,再将x1投入生产方式投入生产方式A和和B,则可得到收入,则可得到收入g(y1)+h(x1-y1),继续回收资源继续回收资源x2=ay1+b(x1-y1), 若上面的过程进行若上面的过程进行n个阶段,我们希望选择个阶段,我们希望选择n个变量个变量y,y1,y2,yn-1,使这,使这n个阶段的总收入最大。个阶段的总收入最大。 例例1 1 续(续(1 1)38 因此,我们的问题就变成:求因此,我们的问题就变成:求y,y1,y2,yn-1,以使以使g(y)+h(x-y)+ g(y1)+h(x1-y1)+ +g(yn-1)+h(xn-1-yn-
17、1) 达到最大,且满足条件达到最大,且满足条件 x1=ay+b(x-y) x2=ay1+b(x1-y1) xn-1=ayn-2+b(xn-2-yn-2) yi与与xi均非负均非负,i=1,2, ,n-1 例例1 1 续(续(2 2)39例例2 2 生产和存储控制问题生产和存储控制问题某工厂生产某种季节性商品,需要作下一某工厂生产某种季节性商品,需要作下一年度的生产计划,假定这种商品的生产周期需年度的生产计划,假定这种商品的生产周期需要两个月,全年共有要两个月,全年共有6个生产周期,需要作出个生产周期,需要作出各个周期中的生产计划。各个周期中的生产计划。40设已知各周期对该商品的需要量如下表所示
18、设已知各周期对该商品的需要量如下表所示:周期周期123456需求量需求量551030508例例2 2 续(续(1 1)41例例2 2 续(续(2 2) 假设这个工厂根据需要可以日夜两班生产或只是日假设这个工厂根据需要可以日夜两班生产或只是日班生产,当开足日班时,每一个生产周期能生产商品班生产,当开足日班时,每一个生产周期能生产商品1515个单位,每生产一个单位商品的成本为个单位,每生产一个单位商品的成本为100100元。当开足元。当开足夜班时,每一生产周期能生产的商品也是夜班时,每一生产周期能生产的商品也是1515个,但是由个,但是由于增加了辅助性生产设备和生产辅助费用,每生产一单于增加了辅助
19、性生产设备和生产辅助费用,每生产一单位商品的成本为位商品的成本为120120元。由于生产能力的限制,可以在元。由于生产能力的限制,可以在需求淡季多生产一些商品储存起来以备需求旺季使用,需求淡季多生产一些商品储存起来以备需求旺季使用,但存储商品是需要存储但存储商品是需要存储费用的,假设每单位商品存储一费用的,假设每单位商品存储一周期需要周期需要16元,已知开始时存储为零,年终也不存储商元,已知开始时存储为零,年终也不存储商品备下年使用,问应该如何作生产和存储计划,才能使品备下年使用,问应该如何作生产和存储计划,才能使总的生产和存储费用最小?总的生产和存储费用最小?42例例2 2 续(续(3 3)
20、设第设第i个周期的生产量为个周期的生产量为xi,周期末的存储量为,周期末的存储量为ui,那,那么这个问题用式子写出来就是:求么这个问题用式子写出来就是:求x1,x2,x6,满足条件:,满足条件: x1=5+u1 x2+u1=5+u2 x3+u2=10+u3 x4+u3=30+u4 x5+u4=50+u5 x6+u5=8 0 xi 30,0 uj,i=1,2,6;j=1,2, ,5 使使S= = 为最小,其中为最小,其中 6511()1 6ijijfxu)1852345(16)(5432161xxxxxxfii433015,300120150 ,100)(iiiiixxxxxf例例2 2 续(续
21、(4 4)44逆序递推方程:0)(1 , 2 , 3 , 4 , 5)(),(min)(6611xfkxfuxdxfkkkkkukkk(1)k=5时,状态,215EEX它们到F点的距离即为最短路。, 4)(15Ef; 3)(25Ef,)(1*5FEu.)(2*5FEuAB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143例1:1阶段2阶段3阶段4阶段5阶段45AB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143(2)k=4时,状态,3214DDDX它们到F点需经过中途点E,需一一分析从E到F的最短路:先说从D
22、1到F的最短路有两种选择:经过E1,E2,比较最短。.)(2*5FEu,)(1*5FEu46AB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143)(),(),(),(min)(252141511414EfEDdEfEDdDf. 735 , 43min这说明由D1到F的最短距离为7,其路径为.11FED相应的决策为:.)(11*4EDu.)(2*5FEu,)(1*5FEu47)(),(),(),(min)(252241512424EfEDdEfEDdDf. 532 , 46min这说明由D2到F的最短距离为5,其路径为.22FED相应的决策为:.)(2
23、2*4EDuAB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143.)(2*5FEu,)(1*5FEu.)(11*4EDu48AB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143)(),(),(),(min)(252341513434EfEDdEfEDdDf. 533, 41min即D3到F的最短距离为5,其路径为.22FED相应的决策为:.)(13*4EDu.)(2*5FEu,)(1*5FEu.)(11*4EDu.)(22*4EDu49(3)k=3时,状态,43214CCCCX )(),(),(),(min
24、)(242131411313DfDCdDfDCdCf.1258, 75min这说明由C1到F的最短距离为12,相应的决策为.)(11*3DCuAB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143.)(2*5FEu,)(1*5FEu.)(11*4EDu.)(22*4EDu.)(13*4EDu5)(24Df7)(14Df50AB1B2C1C2C3C4D1D2D3E1E2F452368775845348435623143)(),(),(),(min)(242231412323DfDCdDfDCdCf.1055 , 74min即由C2到F的最短距离为10,相
25、应的决策为.)(22*3DCu)(),(),(),(min)(343332423333DfDCdDfDCdCf. 854 , 53min.)(2*5FEu,)(1*5FEu.)(11*4EDu.)(22*4EDu.)(13*4EDu.)(11*3DCu7)(14Df5)(24Df5)(34Df51即由C3到F的最短距离为8,相应的决策为.)(23*3DCu)(),(),(),(min)(343432424343DfDCdDfDCdCf. 954 , 58min即由C4到F的最短距离为9,相应的决策为.)(34*3DCuAB1B2C1C2C3C4D1D2D3E1E2F45236877584534
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 运筹学 优化 算法 ppt 课件
限制150内