强变异遗传和近似模拟退火混合算法 源代码.docVIP

// 111.cpp : Defines the entry point for the console application. // #include stdafx.h #include time.h #include iostream #include cmath #include stdio.h #include cstdlib #include iomanip #include string using namespace std; #define max_np 200 /* 初始种群中个体的最大数目(150) */ #define max_nv 6 /* 种群个体中基因的最大数目 */ #define alpha_prob 0.98 /* 概率降低系数 */ #define Pc1 0.99 /* 交叉概率1 */ #define Pc2 0.60 /* 交叉概率2 */ #define EPS 0.95 /* 求解精度 */ int cur_gen; double cur_pop[max_np][max_nv]; /* 当代种群 */ double pop_fitness[max_np]; /* 种群适应度 */ /***********************************************/ /* 比较两个数值的大小，返回较大值的函数 */ /***********************************************/ double MaxValue(double a, double b) { if (a = b) return(a); else return(b); } /***********************************************/ /* 比较两个数值的大小，返回较小值的函数 */ /***********************************************/ double MinValue(double a, double b) { if (a b) return(b); else return(a); } /***********************************************/ /*求解方程组 */ /***********************************************/ void MyFun1(double x[], double y[], int n, double coef[]) { y[0] = coef[1] * sin(x[3]) - coef[0] * sin(x[2]) - x[0]; y[1] = coef[1] * cos(x[1]) * cos(x[3]) - coef[0] * cos(x[2]) + 18 * sin(x[1]); y[2] = ((tan(coef[4]) * sin(x[3]) * sin(x[1]) + cos(x[2]) * cos(x[3]) + cos(x[1])* sin(x[2]) * sin(x[3]))) / tan(coef[3]) - (sin(x[2]) * cos(x[3])-cos(x[1]) * cos(x[2]) * sin(x[3])); y[3] = coef[1] * sin(2 * x[4]) * pow(sin(x[1]), 2) / 2 - 18 * sin(2 * x[1]) * sin(x[4]) / 2 - x[0] * cos(x[4]); y[4] = 4 * pow((coef[1] * cos(x[4]) * cos(x[1]) +18 * sin(x[1])), 2) + 4 * pow((coef[1] * sin(x[4]) - x[0]), 2) - pow(coef[2], 2); n=n; return; } /***********************************************/ /* 初始化适应度