遗传算法求解f(x)=xcosx+2的最大值.docVIP

遗传算法求解f(x)=xcosx+2的最大值 其中在尺度变换部分应用到了类似模拟退火算法部分，所有变量均使用汉语拼音很好懂 //中国电子科技集团公司 //第一研究室 //呼文韬 //hu_hu605@163.com //随机初始种群 //编码方式为格雷码 //选择方法为随机遍历 //采用了精英保存策略 //采用了自适应的交叉率和变异率 //采用了与模拟退火算法相结合的尺度变换 //采用了均匀交叉法 #include stdlib.h #include stdio.h #include math.h #include iostream.h #include iomanip.h #include time.h #include windows.h #define IM1 2147483563 #define IM2 2147483399 #define AM (1.0/IM1) #define IMM1 (IM1-1) #define IA1 40014 #define IA2 40692 #define IQ1 53668 #define IQ2 52774 #define IR1 12211 #define IR2 3791 #define NTAB 32 #define NDIV (1+IMM1/NTAB) #define EPS 1.2e-7 #define RNMX (1.0-EPS) #define zhizhenjuli 0.005 #define PI 3.14159265358 #define T0 100000//温度要取得很高才行。 #define zhongqunshu1 200 #define zuobianjie -2000 #define youbianjie 2000 unsigned int seed=0; //seed 为种子，要设为全局变量 void mysrand(long int i) //初始化种子 { seed = -i; } long a[1]; //double hundun; //double c=4; //设置全局变量 struct individual { unsigned *chrom; //染色体； double geti;//变量值 double shiyingdu; //目标函数的值； double fitness; //变换后的适应度值； }; individual *zuiyougeti;//精英保存策略 int zhongqunshu; //种群大小 individual *nowpop;//当前代 individual *newpop;//新一代 double sumfitness;//当代的总适应度fitness double sumshiyingdu;//当代的总适应度shiyingdu double maxfitness;//最大适应度 double avefitness;//平均适应度 double maxshiyingdu;//最大适应度 double avgshiyingdu;//平均适应度 float pc;//交叉概率 float pm;//变异概率 int lchrom;//染色体长度 int maxgen;//最大遗传代数 int gen;//遗传代数 //函数 int flipc(double ,double );//判断是否交叉 int flipm(double );//判断是否变异 int rnd(int low,int high);//产生low与high之间的任意数 void initialize();//遗传算法初始化 void preselectfitness(); //计算sumfiness,avefitness,maxfitness void generation(); double suijibianli();//产生随机遍历指针 int fuzhi(float );//选择要复制的个体 void crossover(individual ,individual ,individual ,individual );//交叉 void bianyi(individual );//变异 void mubiaohanshu(individual );//计算适应度 void chidubianhuan(individual );//对shiyingdu进行尺度变换赋给fitness double ran1(long *);//随机数初始 void bianma(d