资源描述
遗传算法的思想,Darwin的进化论-“自然选择、适者生存”特定环境的考验种群中个体的选择种群中的交叉繁殖种群中个体的变异上述操作反复执行,个体逐渐优化,遗传算法的手工模拟计算示例为更好地理解遗传算法的运算过程,下面用手工计算来简单地模拟遗传算法的各个主要执行步骤。例:求下述二元函数的最大值:maxf(x1,x2)=x12+x22s.t.x11,2,3,4,5,6,7x21,2,3,4,5,6,7(1)个体编码遗传算法的运算对象是表示个体的符号串,所以必须把变量x1,x2编码为一种符号串。本题中,用无符号二进制整数来表示。因x1,x2为07之间的整数,所以分别用3位无符号二进制整数来表示,将它们连接在一起所组成的6位无符号二进制数就形成了个体的基因型,表示一个可行解。例如,基因型X101110所对应的表现型是:x5,6。个体的表现型x和基因型X之间可通过编码和解码程序相互转换。,(2)初始群体的产生遗传算法是对群体进行的进化操作,需要给其淮备一些表示起始搜索点的初始群体数据。本例中,群体规模的大小取为4,即群体由4个个体组成,每个个体可通过随机方法产生。如:011101,101011,011100,111001(3)适应度汁算遗传算法中以个体适应度的大小来评定各个个体的优劣程度,从而决定其遗传机会的大小。本例中,目标函数总取非负值,并且是以求函数最大值为优化目标,故可直接利用目标函数值作为个体的适应度。(4)选择运算选择运算(或称为复制运算)把当前群体中适应度较高的个体按某种规则或模型遗传到下一代群体中。一般要求适应度较高的个体将有更多的机会遗传到下一代群体中。,本例中,我们采用与适应度成正比的概率来确定各个个体复制到下一代群体中的数量。其具体操作过程是:先计算出群体中所有个体的适应度的总和fi(i=1.2,M);其次计算出每个个体的相对适应度的大小fi/fi,它即为每个个体被遗传到下一代群体中的概率,每个概率值组成一个区域,全部概率值之和为1;最后再产生一个0到1之间的随机数,依据该随机数出现在上述哪一个概率区域内来确定各个个体被选中的次数。,(5)交叉运算交叉运算是遗传算法中产生新个体的主要操作过程,它以某一概率相互交换某两个个体之间的部分染色体。本例采用单点交叉的方法,其具体操作过程是:先对群体进行随机配对;其次随机设置交叉点位置;最后再相互交换配对染色体之间的部分基因。,可以看出,其中新产生的个体“111101”、“111011”的适应度较原来两个个体的适应度都要高。,(6)变异运算变异运算是对个体的某一个或某一些基因座上的基因值按某一较小的概率进行改变,它也是产生新个体的一种操作方法。本例中,我们采用基本位变异的方法来进行变异运算,其具体操作过程是:首先确定出各个个体的基因变异位置,下表所示为随机产生的变异点位置,其中的数字表示变异点设置在该基因座处;然后依照某一概率将变异点的原有基因值取反。,对群体P(t)进行一轮选择、交叉、变异运算之后可得到新一代的群体p(t+1)。,从上表中可以看出,群体经过一代进化之后,其适应度的最大值、平均值都得到了明显的改进。事实上,这里已经找到了最佳个体“111111”。,注意需要说明的是,表中有些栏的数据是随机产生的。这里为了更好地说明问题,我们特意选择了一些较好的数值以便能够得到较好的结果,而在实际运算过程中有可能需要一定的循环次数才能达到这个最优结果。,遗传算法的一个实例,求解方程:将方程求解问题转化为生存问题:解一定在0,10之间,将区间0,10划分成若干个小区间,设想每个小区间为一个生物个体,使下列表达式最小的个体有最好的生存能力,该个体即为解。,遗传算法的一个实例,如何找到这个最优个体?可通过Darwin的进化论由初始个体经过代代演化,逐渐进化出来。如何将生物进化的操作转化为计算机可以执行的操作?通过编码表征生物个体,则生物之间的演化转化为编码的变化。,步骤一:初始化,个体编码:(假定要求小数点后两位)将0,10划分为1024个小区间个体10000000000个体20000000001个体30000000010个体10241111111111种群初始化:随机生成m个10位二进制串,0,10,定义适应度函数:选择(适应度较大的个体),步骤二:选择,为何取倒数?,0.1,0.3,0.2,0.4,A,B,C,D,随机产生0,1之间的数RN,选择个体,选中的优势个体进行交叉-由父个体生成子个体,步骤三:交叉,相同的两个父个体生成相同的两个子个体,变异操作在个体中随机选择一位,改变该位的值,步骤四:变异,交叉和变异操作均以一定概率进行,反复执行步骤二、三、四并记录最优个体(适应度最大的个体)程序结束时,最优个体即为所求解程序结束的判定根据循环次数根据最大适应度根据种群中相同个体数与总个体数的比值,步骤五,遗传算法各步骤的评价,选择-优胜劣汰选择操作为种群提供了演进的方向交叉-优优组合交叉操作的作用在于汇集散布于不同个体间的局部优势模式变异-寻找新模式变异操作是种群向外扩展的触角(随机)好的变异将保留,坏的淘汰,遗传算法的总体评价,优点解决问题的方法具有普适性全局收敛性(依概率收敛)能解决的问题范围很广不足求得的解为近似的数值解对于经典数学可以解决的问题,效率较低,遗传算法的适应度函数,求函数的全局极小值取的初始区间,例如:-10,10将此区间分为1024个小区间,然后编码若求全局极大值(且为正),可直接取函数值为其适应度值,据此作概率选择;若求全局极小值(且为正),可取函数值的倒数为其适应度值,据此作概率选择。若不全为负,可统一加上一个正数,使为正。,TSP问题的遗传算法求解,步骤一:个体编码及种群初始化步骤二:适应度选择步骤三:交叉操作步骤四:变异操作步骤五:重复二、三、四步,直至结束令城市(点)数目为N,个体编码取长度为N的数字串,串中数字互不重复,取值范围为1,N之间的整数。则每一个数字串代表一个个体,个体中数字出现的位置表征路径中城市出现的顺序。初始种群令种群中有M个个体,可随机产生M个数字串构成初始种群。例如:将数字串1234N上的数字进行随机的交换,步骤一:初始化,适应度的计算,步骤二:适应度选择,1,2,3,4,5,对于个体,适应度为:,被选中作为父个体的概率:,选择M次重新生成种群,TSP中交叉算子的特点要保证生成的解为有效解从一个父个体中随机选取一段子串A,在另一个父个体中将A中出现的数字去掉形成串B,AB为一个子串,步骤三:交叉操作,此外还有多种交叉算子,常用的变异操作:随机选取两个相邻位置的数字,交换其顺序。51243(5)51234(5),步骤四:变异操作,1,2,3,4,5,1,2,3,4,5,交换3,4,此外还有多种变异算子,反复执行步骤二、三、四结束判定循环执行G次(例如G=500)后当最优个体的总路径长度小于预期时,步骤五:,中国各省会城市的运行结果,1,2,3,4,5,缺陷:相同父个体生成不同的子个体以下是相同个体:12345(1)54321(5)反射操作12345(1)34512(3)旋转操作,交叉算子的进一步研究,用群论描述所有路径的集合形成一个二面体群A等价解构成一个正规子群BA中陪集的数目为2N,12345(1)32154(3)相同父个体交叉34215(3)15234(1)不同子个体,且和父个体不同,1,2,3,4,5,1,2,3,4,5,1,2,3,4,5,4.1遗传算法简介,智能优化计算,简单实例产生初始种群计算适应度,4.1.4遗传算法的基本操作,0001100000010111100100000001011001110100101010101011100101101001011011110000000110011101000001010011,(8)(5)(2)(10)(7)(12)(5)(19)(10)(14),4.1遗传算法简介,智能优化计算,简单实例选择,4.1.4遗传算法的基本操作,0.086957,0.054348,0.0217390.1086960.0760870.1304350.0543480.2065220.1086960.152174,4.1遗传算法简介,智能优化计算,简单实例选择,4.1.4遗传算法的基本操作,0.086957,0.054348,0.0217390.1086960.0760870.1304350.0543480.2065220.1086960.152174,0.086957,0.141304,0.163043,0.2717390.3478260.4782610.5326090.7391300.8478261.000000,4.1遗传算法简介,智能优化计算,简单实例选择在01之间产生一个随机数:,4.1.4遗传算法的基本操作,0.086957,0.054348,0.0217390.1086960.0760870.1304350.0543480.2065220.1086960.152174,0.086957,0.141304,0.163043,0.2717390.3478260.4782610.5326090.7391300.8478261.000000,0.070221,0.545929,0.784567,0.446930,0.507893,0.291198,0.716340,0.270901,0.371435,0.854641,0001100000111001011011000000011001110100101010101011100101101001011011110000000110011101000001010011,4.1遗传算法简介,智能优化计算,简单实例交叉,4.1.4遗传算法的基本操作,0001100000111001011011000000011001110100101010101011100101101001011011100111010011000000010001010011,0001,1110,100000,010110,111,100,0010110,1011011,110000,100111,0100,0001,1001110100,1100000001,1010101,0001010,010,011,4.1遗传算法简介,智能优化计算,简单实例变异,4.1.4遗传算法的基本操作,4.2基本遗传算法,智能优化计算,问题的提出一元函数求最大值:,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,问题的提出用微分法求取f(x)的最大值:解有无穷多个:,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,问题的提出当i为奇数时xi对应局部极大值点,i为偶数时xi对应局部极小值。x19即为区间-1,2内的最大值点:此时,函数最大值f(x19)比f(1.85)=3.85稍大。,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,编码表现型:x基因型:二进制编码(串长取决于求解精度)串长与精度之间的关系:若要求求解精度到6位小数,区间长度为2-(-1)3,即需将区间分为3/0.000001=3106等份。所以编码的二进制串长应为22位。,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,产生初始种群产生的方式:随机产生的结果:长度为22的二进制串产生的数量:种群的大小(规模),如30,50,111101001110000101100011001100111010101011101010100011110010000100101111001001110011100100011001010011000000110000011010010000000000,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,计算适应度不同的问题有不同的适应度计算方法本例:直接用目标函数作为适应度函数将某个体转化为-1,2区间的实数:s=x=0.637197计算x的函数值(适应度):f(x)=xsin(10x)+2.0=2.586345,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,计算适应度二进制与十进制之间的转换:第一步,将一个二进制串(b21b20b0)转化为10进制数:第二步,x对应的区间-1,2内的实数:,4.2.1简单函数优化的实例,(0000000000000000000000)-1(1111111111111111111111)2,4.2基本遗传算法,智能优化计算,遗传操作选择:轮盘赌选择法;交叉:单点交叉;变异:小概率变异,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,模拟结果设置的参数:种群大小50;交叉概率0.75;变异概率0.05;最大代数200。得到的最佳个体:smax=;xmax=1.8506;f(xmax)=3.8503;,4.2.1简单函数优化的实例,4.2基本遗传算法,智能优化计算,模拟结果进化的过程:,4.2.1简单函数优化的实例,
展开阅读全文