半变异函数拟合指数模型.docVIP

下载本文档

35
0
约6.57千字
约 6页
2018-05-28 发布于河南
举报
版权申诉

半变异函数拟合指数模型.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

半变异函数拟合指数模型

半边也函数的应用半变异函数拟合指数模型程序(c++代码) #include stdio.h #include time.h #includewindows.h #include math.h #include stdlib.h #define S 1 /*试验次数 */ #define G 2000 /*混合迭代次数*/ #define P 200 /*个体总数*/ #define M 20 /*族群数*/ #define I 10 /*因此，一个族群中的个体数是10*/ #define V 3 /*个体维数*/ #define N 10 /*族群内更新次数*/ #define MAX 10 #define MIN 0 double DMAX=1.0; /*蛙跳的最大值*/ double DMIN=0.4; /*蛙跳的最大值*/ double D=MAX/1; /*蛙跳的最大值*/ int i1,i2,i3,i4,ii; int try_number=0; int try_max=5; double C=1.0; #define R ((double)(rand()%10000)/10000)//0-1之间的随机数,精度为1/10000 //#define R1 rand()%100/100.0 static int kk; double PI=3 double Tolerance=0.0000001;//收敛精度 double c3=0.03;//扰动幅度 double e=2.718281828459;//自然对数底数 int sm=3; int bz=0;//扰动因子标志 double aw[V]; double nihe[17][2]={ 1115.658026,8 1915.362904,8 2467.305693,9.1856689, 2952.330784,9.0543057, 5095.207855,9.132906445, 5418.830566,8.852431395, 4146.89209,9 6037.806376,9.103558859, 4818.459044,7.2313171, 5143.558017,9.0538129, 5459.844361,9 5762.570046,8.6310193, 6060.453719,9.194387, 6356.051127,10.398948, 6651.015103,9.8449629, 6941.254523,7.2279982, 7223.868903,6.579128}; typedef struct { double d[V]; double fitness; }Individal; typedef struct { double h[V]; }heli; Individal pw[M];/*族群中个体最差位置*/ Individal pb[M];/*族群中个体最好位置*/ Individal px;/*全体中最好位置*/ Individal individual[P];/*全部个体*/ Individal pop[M][I];/*排序后的群组*/ Individal temp[M]; Individal temp1[I]; Individal tem; Individal temx[S];/*计算标准差*/ /*选择测试函数为Sphere*/ double fitness(double a[]) { int i; double sum=0.0; double sum1=0.0; double s1=0.0,h1=0.0; double x1[V+1]; for(i=0;iV;i++)x1[i]=a[i]; for(i=0;iV;i++) for(i=0;i17;i++) { if(nihe[i][0]x1[2]) sum1=x1[0]+x1[1]; else sum1=nihe[i][1]-(x1[0]+x1[1]*(1.5*nihe[i][0]-0.5*pow(nihe[i][0],3)/pow(x1[2],3))); sum=sum+sum1; } return sum; } /*对每一个个体初始化*/ void init() { int i,j; srand((unsigned)time(NULL)+kk++); for(i=0;iP;i++) { for(j=0;