运筹学的优化算法ppt课件.ppt
《运筹学的优化算法ppt课件.ppt》由会员分享,可在线阅读,更多相关《运筹学的优化算法ppt课件.ppt(199页珍藏版)》请在三一办公上搜索。
1、1,运筹学的优化算法,2,运 筹 学(Operations Research OR) 由于运筹学研究的广泛性和复杂性,人们至今没有形成一个统一的定义。以下给出几种定义:运筹学是一种科学决策的方法运筹学是依据给定目标和条件从众多方案中选择最优方案的最优化技术。运筹学是一种给出坏的问题的答案的艺术,否则的话问题的结果会更坏。,3,运筹学模型 运筹学研究的模型主要是抽象模型数学模型。,4,运筹学包含的分支数学规划(线性规划、整数规划、 目标规划、动态规划、网络规划等)图论与网络流决策分析,5,运筹学包含的分支排队论可靠性数学理论库存论对策论搜索论计算机模拟等,6,数学建模竞赛中的算法(1),93A
2、非线性交调的频率设计: 拟合、规划93B 足球队排名次: 矩阵论、图论、层次分析法、整数规划94A 逢山开路: 图论、插值、动态规划94B 锁具装箱问题: 图论、组合数学95A 飞行管理问题 : 非线性规划、线性规划95B 天车与冶炼炉的作业调度: 非线性规划、动态规划、层次分析法、PETRI方法、图论方法、排队论方法96A 最优捕鱼策略:微分方程、积分、非线性规划,7,96B 节水洗衣机:非线性规划97A 零件参数设计:微积分、非线性规划、随机模拟97B 截断切割:组合优化、几何变换、枚举、蒙特卡罗、递归、最短路98A 投资收益与风险:线性规划、非线性规划98B 灾情巡视路线:最小生成树、H
3、amilton圈、旅行商问题99A 自动化车床管理:积分、概率分布、随机模拟、分布拟合度检验,数学建模竞赛中的算法(2),8,99B 钻井布局:几何变换、枚举、最大完全子图、混合整数规划00A DNA分类:神经网络、最小二乘拟合、统计分类00B 管道订购和运输:最短路、二次规划01A 血管的三维重建:数据挖掘、曲面重建与拟合01B 公交车调度:非线性规划02A 车灯光源优化设计:最优化02B 彩票中的数学:概率与优化,数学建模竞赛中的算法(3),9,03ASARS的传播: 微分方程、差分方程 03B露天矿生产的车辆安排: 整数规划、运输问题04A 奥运会临时超市网点设计:统计分析、数据处理、优
4、化 04B 电力市场的输电阻塞管理:数据拟合、优化 05A 长江水质的评价和预测: 预测评价、数据处理 05BDVD在线租赁: 随机规划、整数规划 06A 出版社书号问题: 整数规划、数据处理、优化,数学建模竞赛中的算法(4),10,06B Hiv病毒问题: 线性规划、回归分析07A 人口问题: 微分方程、数据处理、优化07B 公交车问题: 多目标规划、动态规划、图 论、0-1规划08A 照相机问题: 非线性方程组、优化08B 大学学费问题: 数据收集和处理、统计分 析、回归分析,数学建模竞赛中的算法(5),11,赛题发展的特点: 1. 对选手的计算机能力提出了更高的要求:赛题的解决依赖计算机
5、,题目的数据较多,手工计算不能完成,计算机模拟和以算法形式给出最终结果。 2. 赛题的开放性增大:解法的多样性,一道赛题可用多种解法。开放性还表现在对模型假设和对数据处理上。 3. 试题向大规模数据处理方向发展 4. 求解算法和各类现代算法的融合,12,1. 蒙特卡罗方法(Monte-Carlo方法, MC),数学建模竞赛常用算法(1),该算法又称计算机随机性模拟方法,也称统计试验方法。MC方法是一种基于“随机数”的计算方法,能够比较逼真地描述事物的特点及物理实验过程,解决一些数值方法难以解决的问题。,MC方法的雏型可以追溯到十九世纪后期的蒲丰随机投针试验,即著名的蒲丰问题。 MC方法通过计算
6、机仿真(模拟)解决问题,同时也可以通过模拟来检验自己模型的正确性,是比赛中经常使用的方法。,13,97年的A题 每个零件都有自己的标定值,也都有自己的容差等级,而求解最优的组合方案将要面对着的是一个极其复杂的公式和108种容差选取方案,根本不可能去求解析解,那如何去找到最优的方案呢?随机性模拟搜索最优方案就是其中的一种方法,在每个零件可行的区间中按照正态分布随机的选取一个标定值和选取一个容差值作为一种方案,然后通过蒙特卡罗算法仿真出大量的方案,从中选取一个最佳的。02年的B题 关于彩票第二问,要求设计一种更好的方案,首先方案的优劣取决于很多复杂的因素,同样不可能刻画出一个模型进行求解,只能靠随
7、机仿真模拟。,数学建模竞赛常用算法,14,98 年美国赛A 题 生物组织切片的三维插值处理94 年A 题逢山开路 山体海拔高度的插值计算,数学建模竞赛常用算法(2),2. 数据拟合、参数估计、插值等数据处理算法,比赛中通常会遇到大量的数据需要处理,而处理数据的关键就在于这些算法,通常使用MATLAB 作为工具。与图形处理有关的问题很多与拟合有关系。,此类问题在MATLAB中有很多函数可以调用,只有熟悉MATLAB,这些方法才能用好。,15,98年B 题 用很多不等式完全可以把问题刻画清楚,数学建模竞赛常用算法(3),3. 规划类问题算法,此类问题主要有线性规划、整数规划、多目标规划、二次规划等
8、。竞赛中很多问题都和数学规划有关,可以说不少的模型都可以归结为一组不等式作为约束条件、几个函数表达式作为目标函数的问题,遇到这类问题,求解就是关键了。,因此列举出规划后用Lindo、Lingo 等软件来进行解决比较方便,所以还需要熟悉这两个软件。,16,98 年B 题、00年B 题、95 年锁具装箱等问题体现了图论问题的重要性。,数学建模竞赛常用算法(4),4. 图论问题,这类问题算法有很多,包括:Dijkstra、Floyd、Prim、Bellman-Ford,最大流,二分匹配等问题。,17,92 年B 题用分枝定界法97 年B 题是典型的动态规划问题98 年B 题体现了分治算法,数学建模竞
9、赛常用算法(5),5. 计算机算法设计中的问题,计算机算法设计包括很多内容:动态规划、回溯搜索、分治算法、分枝定界等计算机算法.,这方面问题和ACM 程序设计竞赛中的问题类似,可看一下与计算机算法有关的书。,18,分枝定界法原问题的松驰问题:任何整数规划(IP),凡放弃某些约束条件(如整数要求)后,所得到的问题(P) 都称为(IP)的松驰问题。,最通常的松驰问题是放弃变量的整数性要求后,(P)为线性规划问题。,19,去掉整数约束,用单纯形法,20,分枝定界法步骤一般求解对应的松驰问题,可能会出现下面几种情况:若所得的最优解的各分量恰好是整数,则这个解也是原整数规划的最优解,计算结束。,21,分
10、枝定界法步骤一般求解对应的松驰问题,可能会出现下面几种情况:若所得的最优解的各分量恰好是整数,则这个解也是原整数规划的最优解,计算结束。若松驰问题无可行解,则原整数规划问题也无可行解,计算结束。,22,若松驰问题有最优解,但其各分量不全是整数,则这个解不是原整数规划的最优解,转下一步。,23,若松驰问题有最优解,但其各分量不全是整数,则这个解不是原整数规划的最优解,转下一步。从不满足整数条件的基变量中任选 一个xl进行分枝,它必须满足xl xl 或 xl xl +1中的一个,把这两个约束条件加进原问题中,形成两个互不相容的子问题(两分法)。,24,定界:把满足整数条件各分枝的最优目标函数值作为
11、上(下)界,用它来判断分枝是保留还是剪枝。,25,定界:把满足整数条件各分枝的最优目标函数值作为上(下)界,用它来判断分枝是保留还是剪枝。剪枝:把那些子问题的最优值与界值比较,凡不优或不能更优的分枝全剪掉,直到每个分枝都查清为止。,26,例: 用分枝定界法求解:Max Z=4x1+3x2s.t. 3x1+4x2 12 4x1+2x2 9 x1,x2 0 整数用单纯形法可解得相应的松驰问题的最优解(6/5,21/10)Z=111/10为各分枝的上界。,27,0 1 2 3 4,4 3 2 1,x1,x2,分枝:x1 1,x1 2,P1,P2,(6/5,21/10),28,两个子问题:(P1)Ma
12、x Z=4x1+3x2 s.t. 3x1+4x2 12 4x1+2x2 9 x1,x2 0 , x1 1 ,整数用单纯形法可解得相应的(P1)的最优解(1,9/4) Z=10(3/4),29,(P2)Max Z=4x1+3x2 s.t. 3x1+4x2 12 4x1+2x2 9 x1,x2 0 , x1 2 ,整数用单纯形法可解得相应的(P2)的最优解(2,1/2) Z=9(1/2),30,0 1 2 3 4,4 3 2 1,x1,x2,再对(P1) x1 1 (1,9/4)分枝:(P3) x2 2 (P4) x2 3,P1,P2,P3,P4,31,(P1)两个子问题:(P3)Max Z=4x
13、1+3x2 s.t. 3x1+4x2 12 4x1+2x2 9 x1,x2 0 ,x1 1, x2 2整数用单纯形法可解得相应的(P3)的最优解(1,2) Z=10 为上界,32,(P1)两个子问题:(P4)Max Z=4x1+3x2 s.t. 3x1+4x2 12 4x1+2x2 9 x1,x2 0 ,x1 1, x2 3整数用单纯形法可解得相应的(P4)的最优解(0,3) Z=9,33,X1 2,X2 2,X1 1,X2 3,P1:(1,9/4)Z=10(3/4),P4: (0,3) Z=9,P2:(2,1/2) Z=9(1/2),P3: (1,2) Z=10,P:(6/5,21/10)
14、Z=111/10,原问题的最优解(1,2) Z=10,34,多阶段决策过程最优化 多阶段决策过程是指这样一类特殊的活动过程,他们可以按时间顺序分解成若干相互联系的阶段,在每个阶段都要做出决策,全部过程的决策是一个决策序列,所以多阶段决策问题也称为序贯决策问题。,动态规划应用举例,35,最优化原理:最优策略的任一后部子策略都是最优的。,无论以前状态决策如何,从眼下直到最后的诸决策必构成最优子策略。,例1 最短路线问题,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,动态规划应
15、用举例,36,例1 多阶段资源分配问题,设有数量为x的某种资源,将它投入两种生产方式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,在第一阶段生产后回收的总资源为x1=ay+b(x-y),再将x1投入生产方式A和B,则可得到收入g(y1)+h(x1-y1),继续回收资源x2=ay1+b(x1-y1), 若上面的过程进行n个阶段,
16、我们希望选择n个变量y,y1,y2,yn-1,使这n个阶段的总收入最大。,例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-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 续(2),39,例2 生产和存储控制问题,某工厂生产某种季节性商品,需要作下一年度的生产计划,假定这种商品的生产周期需要两个月,全年共有6个生产周期,需要作出各个周期中的
17、生产计划。,40,设已知各周期对该商品的需要量如下表所示:,例2 续(1),41,例2 续(2),假设这个工厂根据需要可以日夜两班生产或只是日班生产,当开足日班时,每一个生产周期能生产商品15个单位,每生产一个单位商品的成本为100元。当开足夜班时,每一生产周期能生产的商品也是15个,但是由于增加了辅助性生产设备和生产辅助费用,每生产一单位商品的成本为120元。由于生产能力的限制,可以在需求淡季多生产一些商品储存起来以备需求旺季使用,但存储商品是需要存储费用的,假设每单位商品存储一周期需要16元,已知开始时存储为零,年终也不存储商品备下年使用,问应该如何作生产和存储计划,才能使总的生产和存储费
18、用最小?,42,例2 续(3),设第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= = 为最小,其中,43,例2 续(4),44,逆序递推方程:,(1)k=5 时,状态,它们到F 点的距离即为,最短路。,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3
19、,5,6,2,3,1,4,3,例1:,1阶段,2阶段,3阶段,4阶段,5阶段,45,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,(2)k=4 时,状态,它们到F 点需经过中途,点E,需一一分析从E 到 F的最短路:先说从D1到F 的最短路,有两种选择:经过 E1, E2, 比较最短。,46,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,这说明由 D
20、1 到F 的最短距离为7,其路径为,相应的决策为:,47,这说明由 D2 到F 的最短距离为5,其路径为,相应的决策为:,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,48,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,即 D3 到F 的最短距离为5,其路径为,相应的决策为:,49,(3)k=3 时,状态,这说明由 C1 到F 的最短距离为12,相应
21、的决策为,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,50,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,即由 C2 到F 的最短距离为10,相应的决策为,51,即由 C3 到F 的最短距离为8,相应的决策为,即由 C4 到F 的最短距离为9,相应的决策为,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6
22、,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,52,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,(4)k=2时,状态,这说明由 B1 到F 的最短距离为13,相应的决策为,53,即由 B2到F 的最短距离为15,相应的决策为,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,54,A,B1,B2,C1,C2,C3,C4,
23、D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,(1)k=1 时,只有一个状态点A, 则,即由 A到F 的最短距离为17,相应的决策为,55,所以最优路线为:,A,B1,B2,C1,C2,C3,C4,D1,D2,D3,E1,E2,F,4,5,2,3,6,8,7,7,5,8,4,5,3,4,8,4,3,5,6,2,3,1,4,3,再按计算顺序反推可得最优决策序列:,56,97年A 题用模拟退火算法00年B 题用神经网络分类算法01年B 题这种难题也可以使用神经网络美国03年B 题伽马刀问题也是目前研究的课题,目前算
24、法最佳的是遗传算法。,数学建模竞赛常用算法(6),6. 最优化理论的三大非经典算法: 模拟退火法(SA)、神经网络(NN)、遗传算法(GA),近几年的赛题越来越复杂,很多问题没有什么很好的模型可以借鉴,于是这三类算法很多时候可以派上用场。,57,97 年A 题、99 年B 题都可以用网格法搜索,数学建模竞赛常用算法(7),网格算法和穷举法一样,只是网格法是连续问题的穷举。此类算法运算量较大。,7. 网格算法和穷举算法,这种方法最好在运算速度较快的计算机中进行,还有要用高级语言来做,最好不要用MATLAB 做网格,否则会算很久。,58,很多问题都是实际来的,数据可以是连续的,而计算机只能处理离散
25、的数据,因此需要将连续问题进行离散化处理后再用计算机求解。比如差分代替微分、求和代替积分等思想都是把连续问题离散化的常用方法。,数学建模竞赛常用算法(8),8. 连续问题离散化的方法,59,数值分析研究各种求解数学问题的数值计算方法,特别是适合于计算机实现方法与算法。,数学建模竞赛常用算法(9),9. 数值分析方法,它的主要内容包括函数的数值逼近、数值微分与数值积分、非线性方程的数值解法、数值代数、常微分方程数值解等。数值分析是计算数学的一个重要分支,把理论与计算紧密结合,是现代科学计算的基础 。,MATLAB等数学软件中已经有很多数值分析的函数可以直接调用。,60,01年A 题中需要你会读B
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 运筹学 优化 算法 ppt 课件
链接地址:https://www.31ppt.com/p-1450165.html