遗传算法程序实现.docVIP

遗传算法程序实现 #include stdlib.h #include stdio.h #include time.h #include math.h #include GA.h unsigned seed=1; //产生一个不大于x的随机数 double crandom(int x) { seed++; srand((unsigned)(seed*time(NULL))); return(rand()%x); } //产生一个0~1的随机数 double random() { seed++; srand((unsigned)(seed*time(NULL))); return (double)rand()/(double)RAND_MAX; } //产生(0.1，0.2，...，1)之间的一个随机数 double crandom(void) { double t; t=random(); if((10*t-floor(10*t))=0.5) return (floor(10*t)+1)/10; else return(floor(10*t)/10); } //初始化种群 void Initial_group(double group[M][N_para], double up_limit[N_para],double low_limit[N_para]) { double beta; int i,j; for (i=0;iM;i++) { for (j=0;jN_para;j++) { beta=crandom(); group[i][j]=low_limit[j]+beta*(up_limit[j]-low_limit[j]); } } } //加载实时测试数据 void Load_data(double realvalue[][4]) { FILE *fp; float x;//changed if((fp=fopen(realtime\\Realdata.txt,r))==NULL) { printf(cannot open the file Realvalue.txt\n); exit(0); } for(int i=0;iN_hit+N_Begin;i++) { for(int j=0;j4;j++)//修改了 { fscanf(fp,%f ,x); if(i=N_Begin)realvalue[i-N_Begin][j]=(double)x;//修改了 } } fclose(fp); } //计算适应度值 void Fit_calculate(double group[][N_para],double realvalue[][4],double fit_degree[],int n) { int i; //double Mcar,m1,r1,r2,L,J1,J2,c1,c2,b; double mcar,m1,l1,J1,b,g,f1,eks,eksd,theta1,theta1d,h; double x[4],x_temp[4],k[4][4],lefA[4]; double err,A11,A12,A21,A22,B1,B2,B3,B4,C11,C12,C21,C22,D2; g=9.8; h=0.005; for (i=0;in;i++) { mcar=group[i][0]; m1=group[i][1]; l1=group[i][2]; J1=group[i][3]; b=group[i][4]; f1=group[i][5]; //printf(mcar=%6.3f,m1=%6.3f,l1=%6.3f,J1=%6.3f,b=%6.3,f1=%6.3\n,,mcar,m1,l1,J1,b,f1); err=0; for (int j=0;j4;j++) //新加的 { x[j]=realvalue[0][j]; x_temp[j]=realvalue[0][j]; } /*for (int j=0;j4;j++)//修改了 { x[j]=Initivalue[j]; x_temp[j]=Initivalue[j]; }*/ for (i