第五章物流系统的智能优化方法精选文档.ppt
《第五章物流系统的智能优化方法精选文档.ppt》由会员分享,可在线阅读,更多相关《第五章物流系统的智能优化方法精选文档.ppt(73页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第五章物流系统的智能优化方法第五章物流系统的智能优化方法本讲稿第一页,共七十三页q课程内容课程内容 模拟退火算法、遗传算法、禁忌搜索算法、神经网络与神经网络优化算法q课程目标课程目标 了解智能优化的模型结构;理解模拟退火算法的收敛性条件;掌握智能优化的流程、操作、算法理论与技术物流系统优化的智能优化方法为复杂物流管理决策问题提供了重要的可行性解决方案。本讲稿第二页,共七十三页模拟退火算法模拟退火算法(Simulated Annealing)1 1、基本思想基本思想 (1)是基于是基于Monte Carlo迭代求解策略的一种随机迭代求解策略的一种随机寻优算法,其出发点是基于物理中固体物质的退寻优
2、算法,其出发点是基于物理中固体物质的退火过程与一般组合优化问题之间的相似性。火过程与一般组合优化问题之间的相似性。(2)结合爬山法和随机行走结合爬山法和随机行走 注:注:SA算法最早是由算法最早是由Metropolis等等(1953)提出提出本讲稿第三页,共七十三页 (3)模拟退火算法在某一初温下,伴随温度参模拟退火算法在某一初温下,伴随温度参数的不断下降,结合概率突跳特性在解空间中数的不断下降,结合概率突跳特性在解空间中随机寻找目标函数的全局最优解,即在局部优随机寻找目标函数的全局最优解,即在局部优解能概率性地跳出并最终趋于全局最优解。模解能概率性地跳出并最终趋于全局最优解。模拟退火算法是一
3、种通用的优化算法,目前已在拟退火算法是一种通用的优化算法,目前已在工程中得到了广泛应用。工程中得到了广泛应用。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第四页,共七十三页2 2、物理退火过程和物理退火过程和MetropolisMetropolis准则准则简单而言,物理退火过程由以下三部分组成:简单而言,物理退火过程由以下三部分组成:加温过程。其目的是增强粒子的热运动,使其偏加温过程。其目的是增强粒子的热运动,使其偏离平衡位置。当温度足够高时,固体将溶解为液离平衡位置。当温度足够高时,固体将溶解为液体,从而消除系统原先可能存在的非均匀态,使体,从而消除系统原先可能存
4、在的非均匀态,使随后进行的冷却过程以某一平衡态为起点。溶解随后进行的冷却过程以某一平衡态为起点。溶解过程与系统的熵增过程联系,系统能量也随温度过程与系统的熵增过程联系,系统能量也随温度的升高而增大。的升高而增大。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第五页,共七十三页等温过程。物理学的知识告诉我们,对于与周等温过程。物理学的知识告诉我们,对于与周围环境交换热量而温度不变的封闭系统,系统围环境交换热量而温度不变的封闭系统,系统状态的自发变化总是朝自由能减少的方向进行,状态的自发变化总是朝自由能减少的方向进行,当自由能达到最小时,系统达到平衡态。当自由能达到最小时
5、,系统达到平衡态。冷却过程。目的是使粒子的热运动减弱并渐趋冷却过程。目的是使粒子的热运动减弱并渐趋有序,系统能量逐渐下降,从而得到低能的晶有序,系统能量逐渐下降,从而得到低能的晶体结构。体结构。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第六页,共七十三页 Metropolis等在等在1953年提出了重要性采样法,即年提出了重要性采样法,即以概率接受新状态。具体而言,在温度以概率接受新状态。具体而言,在温度t,由当前状由当前状态态i产生新状态产生新状态j,两者的能量分别为,两者的能量分别为 ,若,若 则接受新状态则接受新状态j为当前状态;否则,若概率为当前状态;否则
6、,若概率 大于大于 区间内的随机数则仍旧接受新状态区间内的随机数则仍旧接受新状态j为当前为当前状态,若不成立则保留状态,若不成立则保留i为当前状态,其中为当前状态,其中k为为Boltzmann常数。常数。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第七页,共七十三页 这种重要性采样过程在高温下可接受与当前状这种重要性采样过程在高温下可接受与当前状态能量差较大的新状态,而在低温下基本只接态能量差较大的新状态,而在低温下基本只接受与当前能量差较小的新状态,而且当温度趋受与当前能量差较小的新状态,而且当温度趋于零时,就不能接受比当前状态能量高的新状于零时,就不能接受比当前
7、状态能量高的新状态。这种接受准则通常称为态。这种接受准则通常称为Metropolis准则。准则。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第八页,共七十三页2 2、算法步骤算法步骤标准模拟退火算法的一般步骤可描述如下:标准模拟退火算法的一般步骤可描述如下:给给定定初初温温 ,随随机机产产生生初初始始状状态态 ,令令 ;Repeat:Repeat 产生新状态产生新状态 ;模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第九页,共七十三页Until 抽样稳定准则满足;抽样稳定准则满足;退温退温 ,并,并 令令 ;Until 算法终止准则满足;算
8、法终止准则满足;输出算法搜索结果。输出算法搜索结果。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十页,共七十三页3 3、算法关键参数和操作的设定算法关键参数和操作的设定 从算法流程上看,模拟退火算法包括三函数两从算法流程上看,模拟退火算法包括三函数两准则,即状态产生函数、状态接受函数、温度准则,即状态产生函数、状态接受函数、温度更新函数、内循环终止准则和外循环终止准则,更新函数、内循环终止准则和外循环终止准则,这些环节的设计将决定这些环节的设计将决定SA算法的优化性能。此算法的优化性能。此外,初温的选择对外,初温的选择对SA算法性能也有很大影响。算法性能也有很大影
9、响。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十一页,共七十三页状态产生函数状态产生函数 设计状态产生函数(邻域函数)的出发点应该设计状态产生函数(邻域函数)的出发点应该是尽可能保证产生的候选解遍布全部的解空间。是尽可能保证产生的候选解遍布全部的解空间。通常,状态产生函数由两部分组成,即产生候通常,状态产生函数由两部分组成,即产生候选解的方式和候选解产生的概率分布。选解的方式和候选解产生的概率分布。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十二页,共七十三页状态接受函数状态接受函数 状状态态接接受受函函数数一一般般以以概概率率的的
10、方方式式给给出出,不不同同接接受受函函数数的的差差别别主主要要在在于于接接受受概概率率的的形形式式不不同同。设计状态接受概率,应该遵循以下原则:设计状态接受概率,应该遵循以下原则:在固定温度下,接受使目标函数值下降的候选在固定温度下,接受使目标函数值下降的候选解的概率要大于使目标值上升的候选解的概率;解的概率要大于使目标值上升的候选解的概率;模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十三页,共七十三页随温度的下降,接受使目标函数值上升的解的概随温度的下降,接受使目标函数值上升的解的概率要逐渐减小;率要逐渐减小;当温度趋于零时,只能接受目标函数值下降的解。当温度趋
11、于零时,只能接受目标函数值下降的解。状态接受函数的引入是状态接受函数的引入是SA算法实现全局搜索的最算法实现全局搜索的最关键的因素,关键的因素,SA算法中通常采用算法中通常采用min1,exp(-C/t)作为状态接受函数。作为状态接受函数。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十四页,共七十三页初温初温 初始温度、温度更新函数、内循环终止准则和初始温度、温度更新函数、内循环终止准则和外循环终止准则通常被称为退火历程外循环终止准则通常被称为退火历程(annealing schedule)。实验表明,初温越大,。实验表明,初温越大,获得高质量解的几率越大,但花费
12、的计算时间获得高质量解的几率越大,但花费的计算时间将增加。因此,初温的确定应折衷考虑优化质将增加。因此,初温的确定应折衷考虑优化质量和优化效率,常用方法包括:量和优化效率,常用方法包括:模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十五页,共七十三页均匀抽样一组状态,以各状态目标值的方差为初均匀抽样一组状态,以各状态目标值的方差为初温。温。随随机机产产生生一一组组状状态态,确确定定两两两两状状态态间间的的最最大大目目标标值值差差 ,然然后后依依据据差差值值,利利用用一一定定的的函函数数确确定定初初温温。譬譬如如 ,其其中中 为为初初始始接接受受概概率。率。利用经验公
13、式给出。利用经验公式给出。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十六页,共七十三页温度更新函数温度更新函数温温度度更更新新函函数数,即即温温度度的的下下降降方方式式,用用于于在在外外循循环环中修改温度值。中修改温度值。目前,最常用的温度更新函数为指数退温函数,即,目前,最常用的温度更新函数为指数退温函数,即,其中且其大小可以不断变化。其中且其大小可以不断变化。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十七页,共七十三页内循环终止准则内循环终止准则 内循环终止准则,或称内循环终止准则,或称MetropolisMetropolis
14、抽样稳定准抽样稳定准则,用于决定在各温度下产生候选解的数目。则,用于决定在各温度下产生候选解的数目。在非时齐在非时齐SASA算法理论中,由于在每个温度下只算法理论中,由于在每个温度下只产生一个或少量候选解,所以不存在选择内循产生一个或少量候选解,所以不存在选择内循环终止准则的问题。环终止准则的问题。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十八页,共七十三页 而在时齐而在时齐SASA算法理论中,收敛条件要求在每个算法理论中,收敛条件要求在每个温度下产生候选解的数目趋于无穷大,以使相温度下产生候选解的数目趋于无穷大,以使相应的马氏链达到平稳概率分布,显然在实际应应
15、的马氏链达到平稳概率分布,显然在实际应用算法时这是无法实现的。常用的抽样准则包用算法时这是无法实现的。常用的抽样准则包括:括:检验目标函数的均值是否稳定;检验目标函数的均值是否稳定;连续若干步的目标值变化较小;连续若干步的目标值变化较小;按一定的步数抽样。按一定的步数抽样。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第十九页,共七十三页外循环终止准则外循环终止准则 外循环终止准则,即算法终止准则,用于决定算外循环终止准则,即算法终止准则,用于决定算法何时结束。设置温度终值是一种简单的方法。法何时结束。设置温度终值是一种简单的方法。SASA算法的收敛性理论中要求温度终
16、值趋于零,这算法的收敛性理论中要求温度终值趋于零,这显然不合实际。通常的做法是:显然不合实际。通常的做法是:设置终止温度的阈值;设置终止温度的阈值;设置外循环迭代次数;设置外循环迭代次数;算法搜索到的最优值连续若干步保持不变;算法搜索到的最优值连续若干步保持不变;检验系统熵是否稳定。检验系统熵是否稳定。模拟退火算法模拟退火算法(Simulated Annealing)本讲稿第二十页,共七十三页小结小结q由于算法的一些环节无法在实际设计算法时实现,由于算法的一些环节无法在实际设计算法时实现,因此因此SA算法往往得不到全局最优解,或算法结算法往往得不到全局最优解,或算法结果存在波动性。果存在波动性
17、。q目前,目前,SA算法参数的选择仍依赖于一些启发式准算法参数的选择仍依赖于一些启发式准则和待求问题的性质。则和待求问题的性质。SA算法的通用性很强,算法的通用性很强,算法易于实现,但要真正取得质量和可靠性高、算法易于实现,但要真正取得质量和可靠性高、初值鲁棒性好的效果,克服计算时间较长、效果初值鲁棒性好的效果,克服计算时间较长、效果较低的缺点,并适用于规模较大的问题,尚需进较低的缺点,并适用于规模较大的问题,尚需进行大量的研究工作。行大量的研究工作。本讲稿第二十一页,共七十三页1、基本概念、基本概念 模拟生物在自然环境中的遗传和进化过程而形模拟生物在自然环境中的遗传和进化过程而形成的一种成的
18、一种自适应全局优化概率搜索算法。搜索算法。遗传算法遗传算法本讲稿第二十二页,共七十三页2、基本思想基本思想 对对于于一一个个求求函函数数最最大大值值的的优优化化问问题题,一一般般可可描描述为下述数学规划模型:述为下述数学规划模型:遗传算法遗传算法本讲稿第二十三页,共七十三页式中,式中,为决策变量,为决策变量,f(X)为目标函数,为目标函数,U是基本空间,是基本空间,R是是U的一个子集。的一个子集。遗传算法中,将遗传算法中,将n维决策向量用维决策向量用n个记号个记号所组成的符号串所组成的符号串X来表示:来表示:遗传算法遗传算法本讲稿第二十四页,共七十三页 把每一个把每一个 看作一个遗传基因,它的
19、所有可能取值称为等看作一个遗传基因,它的所有可能取值称为等位基因,这样,位基因,这样,X就可看作是由就可看作是由n个遗传基因所组成的一个个遗传基因所组成的一个染色体。染色体的长度可以是固定的,也可以是变化的。染色体。染色体的长度可以是固定的,也可以是变化的。等位基因可以是一组整数,也可以是某一范围内的实数等位基因可以是一组整数,也可以是某一范围内的实数值,或者是记号。最简单的等位基因是由值,或者是记号。最简单的等位基因是由0和和1这两个整这两个整数组成的,相应的染色体就可表示为一个二进制符号串。数组成的,相应的染色体就可表示为一个二进制符号串。遗传算法遗传算法本讲稿第二十五页,共七十三页 这种
20、编码所形成的排列形式这种编码所形成的排列形式X是个体的基因型,与它对应的是个体的基因型,与它对应的X值是个体的表现型。染色体值是个体的表现型。染色体X也称为个体也称为个体X,对于每一个个,对于每一个个体体X,要按照一定的规则确定出其适应度。个体的适应度与,要按照一定的规则确定出其适应度。个体的适应度与其对应的个体表现型其对应的个体表现型X的目标函数值相关联,的目标函数值相关联,X越接近于目标越接近于目标函数的最优点,其适应度越大;反之,其适应度越小。函数的最优点,其适应度越大;反之,其适应度越小。遗传算法遗传算法本讲稿第二十六页,共七十三页 遗传算法中,决策变量遗传算法中,决策变量X组成了问题
21、的解组成了问题的解空间。对问题最优解的搜索是通过对染色体空间。对问题最优解的搜索是通过对染色体X的搜索过程来进行的,从而由所有的染色体的搜索过程来进行的,从而由所有的染色体X就组成了问题的搜索空间。就组成了问题的搜索空间。生物的进化是以集团为主体的。与此相对生物的进化是以集团为主体的。与此相对应,遗传算法的运算对象是由应,遗传算法的运算对象是由M个个体所组成个个体所组成的集合,称为群体。的集合,称为群体。遗传算法遗传算法本讲稿第二十七页,共七十三页与生物一代一代的自然进化过程相似,遗传算法的运算过与生物一代一代的自然进化过程相似,遗传算法的运算过程也是一个反复迭代过程,第程也是一个反复迭代过程
22、,第t代群体记做代群体记做P(t),经过一代遗传和经过一代遗传和进化后,得到第进化后,得到第t+1代群体,它们也是由多个个体组成的集合,记代群体,它们也是由多个个体组成的集合,记做做P(t+1)。这个群体不断地经过遗传和进化操作,并且每次都按照。这个群体不断地经过遗传和进化操作,并且每次都按照优胜劣汰的规则将适应度较高的个体更多地遗传到下一代,这样最终优胜劣汰的规则将适应度较高的个体更多地遗传到下一代,这样最终在群体中将会得到一个优良的个体在群体中将会得到一个优良的个体X,它所对应的表现型,它所对应的表现型X将达到将达到或接近于问题的最优解或接近于问题的最优解 。遗传算法遗传算法本讲稿第二十八
23、页,共七十三页生物的进化过程主要是通过染色体之间的交生物的进化过程主要是通过染色体之间的交叉和染色体的变异来完成的。遗传算法中最优解叉和染色体的变异来完成的。遗传算法中最优解的搜索过程也模仿生物的这个进化过程,使用所的搜索过程也模仿生物的这个进化过程,使用所谓的遗传算子谓的遗传算子(genetic operators)作用于群体作用于群体P(t)中,进行下述遗传操作,从而得到新一代群中,进行下述遗传操作,从而得到新一代群体体P(t+1)。遗传算法遗传算法本讲稿第二十九页,共七十三页q选择选择(selection):(selection):根据各个个体的适应度,按照根据各个个体的适应度,按照一定
24、的规则或方法,从第一定的规则或方法,从第t t代群体代群体P(t)P(t)中选择出一中选择出一些优良的个体遗传到下一代群体些优良的个体遗传到下一代群体P P(t t+1)+1)中。中。q交叉交叉(crossover):(crossover):将群体将群体P(t)P(t)内的各个个体随机内的各个个体随机搭配成对,对每一个个体,以某个概率搭配成对,对每一个个体,以某个概率(称为交叉称为交叉概率,概率,crossover rate)crossover rate)交换它们之间的部分染色交换它们之间的部分染色体。体。遗传算法遗传算法本讲稿第三十页,共七十三页q变异变异(mutation):对群体对群体P
25、(t)中的每一个个体,中的每一个个体,以某一概率以某一概率(称为变异概率,称为变异概率,mutation rate)改变某一个或一些基因座上基因值为其它的等改变某一个或一些基因座上基因值为其它的等位基因。位基因。遗传算法遗传算法本讲稿第三十一页,共七十三页3 3、特点特点遗传算法是一类可用于复杂系统优化计算的鲁遗传算法是一类可用于复杂系统优化计算的鲁棒搜索算法,与其他一些优化算法相比,主要有棒搜索算法,与其他一些优化算法相比,主要有下述几个特点:下述几个特点:q遗传算法以决策变量的编码作为运算对象。传统遗传算法以决策变量的编码作为运算对象。传统的优化算法往往直接利用决策变量的实际值本身的优化算
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第五 物流 系统 智能 优化 方法 精选 文档
限制150内