lw4-基于遗传算法解决01背包问题.doc-邵阳学院.doc

基于遗传算法解决01背包问题 邵阳学院，刘俊，刘希信息工程系：遗传算法【】属于进化算法( Evolutionary Algorithms) 的一种,它通过模仿自然界的选择与遗传的机理来寻找最优解. 遗传算法具有与问题领域无关切快速随机的搜索能力。搜索从群体出发，具有潜在的并行性，可以进行多个个体的同时比较搜索使用评价函数启发，过程简单使用概率机制进行迭代，具有随机性具有可扩展性，容易与其他算法结合。本文基本01背包问题。评价函数问题描述 1.问题描述 01背包问题是一类基础的背包问题，可以具体描叙为一个体积为volume的背包和n个物体，对于物体i，其体积为vi，价值为wi，01背包问题就是要在不超过背包体积的情况下，使装入背包的物体的价值量最大。问题的解的遗传表示【】 因为对于一件物品只有放进与不放进两种选择，所以只用两种状态就可以表示，于是采用二进制编码。当xi=0时，表示物体i不放入背包，当xi=1时表示物体i放入背包，则0-1背包问题的数学模型为，n个物体的背包问题的解可以表示为二进制串。 评价函数根据个体适应值对其进行优劣判定采用经典的赌轮选择算法【】，函数的功能是从群体中选择一个基因组，选中的几率正比于基因组的适应性分数。 SGenome backpack::RouletteWheelSelection(){ srand(time(0)); double fSlice = RandFloat() * m_dTotalFitnessScore; double cfTotal = 0.0; int SelectedGenome = 0; for (int i=0; im_iPopSize; ++i){ cfTotal += m_vecGenomes[i].dFitness; if (cfTotal fSlice) { SelectedGenome = i; break; } } return m_vecGenomes[SelectedGenome]; } 2.遗传算子用来改变繁殖过程中产生的子个体遗传组成 杂交算子【】，要求两个染色体在同一随机位置上断裂，然后将它们在断开点以后的部分进行互换，以形成两个新的染色体（子代）。 void backpack::Crossover( const vectorint mum,const vectorint dad, vectorint baby1,vectorint baby2) { if ( (RandFloat() m_dCrossoverRate) || (mum == dad)) { baby1 = mum; baby2 = dad; return; } int cp = RandInt(0, m_iChromoLength - 1); for (int i=0; icp; ++i){ baby1.push_back(mum[i]); baby2.push_back(dad[i]); } for (int i=cp; imum.size(); ++i){ baby1.push_back(dad[i]); baby2.push_back(mum[i]); } } 变异算子【9】，沿着染色体长度，位进行考察，并按m_dMutationRate给定的几率，将其中某些位进行翻转。 void backpack::Mutate(vectorint vecBits){ for (int curBit=0; curBitvecBits.size(); curBit++){ if (RandFloat() m_dMutationRate){ vecBits[curBit] = !vecBits[curBit]; } }} 可以设置适应性函数，为更好的解设置更高的分数，但这样容易造成数值分布的广度不够，使种群中最好和最差染色体的分数差别不大。我们可以记下每一代中的最差的价值，然后对种群中每一个解得出的价值减去这个最差的价值，这样相对差异变大，接下来使用轮盘选择法就会更有效可有效地从种群中移去最差的染色体，因为最差的适应性分数为0 【】。不同的变异算子【1】散播变异随机地选择一对位置，将其间的基因进行任意的移动。代码如下： void backpack::MutateSM(vectorint chromo){ if (RandFloat() m_dMutationRate) return; const int MinSpanSize = 3; int beg, end; ChooseSection(beg, end, chromo.size()-1,