遗传算法的背包问题c语言

1.3取 ，其中 表示0-1整数的均匀分布函数，即随机地生成数0或1，生成的 串即可看为一个染色体个体。
若不满足模型（*）的约束条件，则拒绝接受，由1.2重新生成一个新的染色体个体chrom；如果产生的染色体可行，则接受它作为种群的一名成员，经过有限次的1.2抽样后，得到popsize个可行的染色体chrom，形成新的种群。
6.3 种群的代数gen=genm+1，若gen〉Maxgen，则结束种群的演化，否则转到步骤2。
六、遗传算法求解的实现
1、遗传算法的主要参数
#define popsize 80 //种群的规模
#define pc 0.7 //杂交概率
#define pm 0.1 //变异概率
#define lchrom 50 //染色体长度
返回值：无；
流程图：见图ቤተ መጻሕፍቲ ባይዱ。
图4、函数statistics的流程图
图5、函数report的流程图
（5）、void initpop( )
功能：生成初始种群；
参数：无；
返回值：无；
流程图：见图6。
图6、函数initpop流程图
（6）、int execise(double probability)
功能：概率选择试验，以概率probability做随机试验，判断是否进行交叉或变异操作；
#define maxgen 5000 //最大进化代数
double alpha;//计算适应度时使用的 惩罚函数系数
2、数据结构
（1）背包信息：
//背包问题中物体重量、收益、背包容量
int weight[lchrom],profit[lchrom],contain;
（2）种群个体结构体
struct population
使用递归回溯法解决背包问题的优点在于它算法思想简单，而且它能完全遍历搜索空间，肯定能找到问题的最优解；但是由于此问题解的总组合数有 个，因此，随着物件数n的增大，其解的空间将以 级增长，当n大到一定程度上，用此算法解决背包问题将是不现实的。
使用贪心方法求解时计算的复杂度降低了很多，但是往往难以得到最优解，有时所得解与最优解相差甚远。因此，我们可以探索使用遗传算法解决物件数较大的背包问题。
遗传算法是将问题的每一个可能性解看作是群体中的一个个体(染色体)，并将每一个染色体编码成串的形式，再根据预定的目标函数对每个个体进行评价，给出一个适应值。算法将根据适应度值进行它的寻优过程，遗传算法的寻优过程是通过选择、杂交和变异三个遗传算子来具体实现的。它的搜索能力由选择算子和杂交算子决定，变异算子则保证了算法能够搜索到问题空间的尽可能多的点，从而使其具有搜索全局最优的能力。遗传算法的高效性和强壮性可由Holland提出的模式定理( Schema Therem)和隐式并行性得以解释。在遗传算法中，定义长度较短、低阶且适应值超过平均适应值的模式在群体中数目的期望值按指数递增，这个结论称为遗传算法的基本定理。遗传算法是通过定义长度短、确定位数少、适应度值高的模式的反复抽样、组合来寻找最佳点，称这些使遗传算法有效工作的模式为积木块，是遗传算法构造答案的基本材料。但归根到底，要使遗传算法有效工作必须按照遗传算法的模式定理(或积木块假设)根据具体问题设计合理的编码方案。
在运行遗传算法程序时，需要对一些参数作事先选择，它们包括种群的大小、染色体长、交叉率、变异率、最大进化代数等，这些参数对GA的性能都有很重要的影响。在试验中参数一般选取如下：种群大小N= 20~100 ，交叉概率 = 0.4 ~0.9 ，变异概率 = 0.001~0.1 ，最大进化代数maxgen = 100~500。
目标函数：
（*）
式中 为0-1决策变量， 时表示将物品 装入背包中， 时则表示不将其装入背包中。
三、求解背包问题的一般方法
解决背包问题一般是采取动态规划、递归回溯法和贪心方法。动态规划可以把困难得多阶段决策变换为一系列相互联系比较容易的单阶段问题。对于背包问题可以对子过程用枚举法求解，而且约束条件越多，决策的搜索范围越小，求解也越容易。它的主要缺点是用数值方法求解时会随着状态变量的个数呈指数级的增长，往往对于求解背包问题的实际问题是不现实的。
（4）适应度信息
//种群的总适应度、最小、最大适应度
doublesumfitness,minfitness,maxfitness;
//一个种群中最大和最小适应度的个体编号
int minpop,maxpop;
3、主要函数说明
（1）、intread_infor()
功能：从文件knapsack.txt中读出背包信息（物体重量、收益、背包容量）；
二、问题描述
背包问题(Knapsack Problem)的一般提法是：已知n个物品的重量（weight）及其价值（或收益profit）分别为 和 ，背包的容量（contain）假设设为 ，如何选择哪些物品装入背包可以使得在背包的容量约束限制之内所装物品的价值最大？
该问题的模型可以表示为下述0/1整数规划模型：
遗传算法是具有“生成+检测”的迭代过程的搜索算法。它的基本处理流程如图1所示。
图1、遗传算法的基本流程
遗传算法的基本流程描述如下：
（1）编码：将解空间的解数据进行二进制编码，表达为遗传空间的基因型串（即染色体）结构数据，如将数据9编码为“1001”；
（2）初始化种群：定义整数pop_size作为染色体的个数，并且随机产生pop_size个染色体作为初始种群；
（3）评估种群中个体适应度：评价函数对种群中的每个染色体（chromosome）求得其个体适应度 ；
（4）选择：选择把当前群体中适应度较高的个体按某种规则或者模型遗传到下一代种群中，这里所用的规则是：染色体在种群中被选择的可能性与其个体的适应度的大小成正比；
（5）交叉：定义参数 作为交叉操作的概率，由（4）选择得到的两个个体以概率 交换各自的部分染色体，得到新的两个个体；
遗传算法已经成为组合优化问题的近似最优解的一把钥匙。它是一种模拟生物进化过程的计算模型，作为一种新的全局优化搜索算法，它以其简单、鲁棒性强、适应并行处理以及应用范围广等特点，奠定了作为21世纪关键智能计算的地位。
背包问题是一个典型的组合优化问题，在计算理论中属于NP-完全问题，其计算复杂度为 ，传统上采用动态规划来求解。设w[i]是经营活动 i 所需要的资源消耗，M是所能提供的资源总量，p[i]是人们经营活动i得到的利润或收益，则背包问题就是在资源有限的条件下，追求总的最大收益的资源有效分配问题。
四、遗传算法简介
遗传算法( Genetic Algorithms，GA)是在1975年首次由美国密西根大学的D。J。Holland教授和他的同事们借鉴生物界达尔文的自然选择法则和孟德尔的遗传进化机制基础之上提出的。经过近30年的研究、应用，遗传算法已被广泛地应用于函数优化、机器人系统、神经网络学习过程、模式识别、图象处理、工业优化控制等领域。
流程图：见图3。
图3、函数cal_fit的流程图
（3）、void statistics(struct population *pop)
功能：群体适应度的最大最小值以及其他信息；
参数：struct population *pop是种群指针，根据指针所指向的种群信息统计群体适应度的信息；
返回值：无；
流程图：见图4。
1.4 置种群的代数gen=0；
步骤2、计算种群中个体适应度以及统计种群适应度情况
2.1 按照下列公式计算种群中个体适应度：
；
公式（2）的下半部分即为适应度的惩罚函数，其中参数 。
2.2 按公式（3）计算种群的总体适应度，
并且按照排序的方法统计出种群中的最大、最小适应度的染色体个体，分别标记为maxpop、minpop；
（6）变异：定义参数 作为变异操作的概率，由（5）得到每个个体中的每个基因值都以概率 进行变异；
（7）演化：经过选择、交叉和变异操作，得到一个新的种群，对上述步骤经过给定的循环次数（maxgen）的种群演化，遗传算法终止。
五、背包问题的遗传算法求解描述
基于背包问题的模型（*），我们设计了针对于背包问题的染色体编码方法：将待求解的各量 表示成长为 的二进制字符串 ，j=1，2， …，n。 表示物体j不放入背包内， …000111代表一个解，它表示将第1、2、3、6、7…n-2,n-1,n号物体放入背包中，其它的物体则不放入。
参数：无；
返回值：返回读取文件信息是否正确；
流程图：见图2。
图2、read_infor( )流程图
（2）double cal_fit(unsigned int *chr)
功能：种群中个体适应度计算；
参数：unsigned int *chr是染色体个体的指针，根据指针所指向的染色体计算个体的适应度；
返回值：染色体个体适应度的大小；
while sum<wheel-pos && i<=popsize
{ i=i+1;
sum=sum+fitness[i]; //fitness为第i个个体的适应度
}
return i-1; //选择了个体i-1
}
3.3 重复两次操作3.1、3.2，生成两个个体作为交叉操作的父代；
步骤四、交叉操作
4.1 根据事先定义好的交叉概率 ，为了确定是否进行交叉操作，则生成[0，1]的随机数pp，若 ，则进行4.2交叉操作，否则将两个父代保留为下一代的两个个体；
根据遗传算法的基本流程，我们确定了求解背包问题的遗传算法：
步骤1、初始化过程
1.1确定种群规模popsize、杂交概率 、变异概率 、染色体长度lchrom及最大进化代数maxgen；
1.2读入背包问题的相关信息，如每个物体的重量weight[j]、每个物体的收益profit[j]和背包的容量contain，其中 ；