一种新型GEP算法及应用研究.docVIP

一种新型的GEP算法及应用研究 　　摘 要：提出一种适合于GEP表达式树构造的新方法，以及相应的新解码方法(GPED)。通过实验对比，GPED可大大缩短演化时间。提出一种新的算法GPEP，将GPEP应用于碎石桩复合地基承载力预测，结果表明GPEP算法在预测精度和演化效率上均超过遗传神经网络、GP等算法。 　　关键词：基因表达式编程; 遗传程序设计; 解码; 预测 　　中图分类号：TP18文献标志码：A 　　文章编号：1001-3695(2007)06-0098-03 　　 　　0 引言?? 　　在进化计算家族中，遗传算法(GA)［1］用固定长度的线性串表示个体。由于遗传操作简单，目前已广泛应用于各行各业的优化组合、图像处理、机器人等问题中；但却丧失了功能的复杂性［2］。在遗传规划(GP)［3,4］中，个体是长度和形状不同的非线性实体（分裂树），所以能够表示复杂问题而被应用于复杂函数建模、时间序列预测等领域；但遗传操作复杂、演化时间长。基因表达式编程(GEP)［2，5，6］则将个体编码成固定长度的线性串（基因组或染色体），然后将其转换为不同长度和形状的非线性实体（表达式树），从而实现了用简单编码表示复杂问题，同时易于遗传操作，并且通过遗传操作所产生的新个体在语法上都是有效的，不需要对新个体进行有效性判断和处理，在速度上比GP提高了2－4个数量级。但GEP在将染色体串表示成表达式树，再遍历表达式树（通常是求后缀表达式，也称逆波兰式）来计算适应度值过程中，由于树的构造和遍历操作非常耗时，大大影响了算法的效率。对此，本文提出了一种新的染色体解码方法（Gene Postfix Expression Decoding，GPED）。该方法只需对染色体作简单处理，即可直接得到染色体相对应的后缀表达式，不需要构造和遍历表达式树，可进一步提高演化效率。文中给出了GPED的有关性质证明，以及与其相对应的表达式树的另一种构造方法。在GPED基础上，提出了一种新的GEP算法GPEP(Gene Postfix Expression Programming)。 　　 　　1 GEP中染色体基因和表达式树的互相转换 ?? 　　GEP的基因由头部和尾部构成。头部可以含有函数符号和终端符号，而尾部只包含终端符号。在具体问题中，头部长度??h是根据求解问题的复杂度而选定的；而尾部长度t则由式（1）计算得到。其中n是带有最多参数的函数的参数个数。 　　从左到右逐个读取该基因中的字符。根据其语义，按层次顺序构成表达式树(ET)［2］，如图1所示。对于给定的表达式树，按从上到下、从左到右的顺序层次遍历ET中的节点，得到的符号序列即为基因编码的有效部分，如遍历图1中的表达式树可得到序列CS*a/eb。 　　 　　2 一种新的GEP解码方法和对应表达式树构造方法及其优点 　　在GEP中，要计算个体适应度值，就必须将个体染色体或基因组转换为表达式树；然后遍历表达式树得到相应的函数式；最后根据适应度函数计算适应度值。通常的做法就是后序遍历表达式树得到后缀表达式，如图1中表达式树的后缀表达式为aeb/*SC。这样就可以由后缀表达式计算表达式值。由于要进行大量的构造树和遍历树操作，成为影响GEP算法效率的瓶颈。本文提出一种新的针对GEP基因的解码方法。 　　定义1(基因后缀表达式解码) 下述过程称为基因后缀表达式解码。首先，设置一个计数量nCount并初始化为1。从左到右逐个读取基因中的字符，按下列规则改变nCount的值，直到nCount为0。 　　性质1 对任意一个有效的GEP，即GEP的头部和尾部满足式(1)，基因后缀表达式解码扫描完基因串之前，计数量nCount必定已等于0。 　　性质1保证了绝不会出现当基因后缀表达式解码扫描完基因串后，nCount 的值还不为0的情况。?? 　　下面介绍一种新的表达式树构造方法。过程如下：基因的首字符对应树的根节点；从左到右逐个读取基因串的字符，并按照从右到左、深度优先的顺序构造表达式树，直到树中所有叶节点都是终端字符。与GEP一样，每个函数节点连接的子分支数目为该节点的参数个数，而终端字符节点只作为叶节点。以式(2)为例，用新方法构造的表达式树如图2所示。按照从右到左、深度优先的顺序遍历该表达式树，即可得到基因的有效编码部分，如按这种方法遍历图2中的表达式树得到CS*a/eb。它正是式(2)的基因有效编码部分。为了与GEP中的表达式树区别，将这种新方法构造的树称为后缀表达式树 (Postfix Expression Tree，PET)。 　　性质2 任意一个有效的GEP基因通过基因后缀表达式解码得到的都是一个在语义上合法的后缀表达式。?? 　　证明 由PET的构造过程可知，若函数的参数个数为??