7.htm

遗传算法的详细描述文档,文档中附带一定的代码解释

    } <br>
  population = matrix(1,pop_size,1,x2_vari); <br>
  print_pop = matrix(1,pop_size,1,x2_vari+tot_eq); <br>
  new_genera = matrix(1,pop_size,1,x2_vari+1); <br>
  temp = matrix(1,2,1,x2_vari); <br>
  probab = vector(1,pop_size); <br>
  t_vec = vector(1,x2_vari); <br>
  cum_probab = vector(1,pop_size); <br>
  eval = vector(1,pop_size); <br>
  live = ivector(1,pop_size); <br>
  die = ivector(1,pop_size); <br>
  for(i=1; i<=x2_vari; i++) <br>
    for(j=1; j<=pop_size; j++) <br>
      population[j][i] = X[x2[i]]; <br>
  fprintf(output,"\nThe initial point of the population is\n"); <br>
  print_vector(X,1,tot_eq+x2_vari); <br>
  fprintf(output,"\n\n"); <br>
  eval[1] = evaluate(X); <br>
  for(j=2; j<=pop_size; j++) <br>
    eval[j] = eval[1]; <br>
  fprintf(output,"\nThe value of the initial point is  %f\n", eval[1]); <br>
  fprintf(output,"\n\n"); <br>
  assign_probab(probab,pop_size,Q); <br>

  find_cum_probab(cum_probab,probab,pop_size); <br>
  Teval = eval[1]; <br>
  do <br>
    { <br>
      for(j=1; j<=pop_size; j++) <br>
 { <br>
   live[j] = die[j] = 0; <br>
   new_genera[j][1] = 0.0; <br>
   for(i=1; i<=x2_vari + 1; i++) <br>
     new_genera[j][i+1] = population[j][i]; <br>
 } <br>
      find_live_die(cum_probab,live,die,pop_size,P); <br>
      for(j=1; j<=P1; j++) <br>
 { <br>
   first = find_parent(live,pop_size); <br>
   new_genera[first][1] = 1.0; <br>
   for(i=1; i<=x2_vari; i++) <br>
     t_vec[i] = new_genera[first][i+1]; <br>
   oper1(t_vec,fin_mat,rc); <br>
   for(i=1; i<=x2_vari; i++) <br>
     new_genera[first][i+1] = t_vec[i]; <br>
 } <br>
 } <br>
      for(j=1; j<=P2; j++) <br>
 { <br>
   first = find_parent(live,pop_size); <br>
   new_genera[first][1] = 1.0; <br>
   for(i=1; i<=x2_vari; i++) <br>
     t_vec[i] = new_genera[first][i+1]; <br>
   oper2(t_vec,fin_mat,rc); <br>
   for(i=1; i<=x2_vari; i++) <br>
     new_genera[first][i+1] = t_vec[i]; <br>
 } <br>
      for(j=1; j<=P3; j++) <br>
 { <br>
   first = find_parent(live,pop_size); <br>
   new_genera[first][1] = 1.0; <br>
   for(i=1; i<=x2_vari; i++) <br>
     t_vec[i] = new_genera[first][i+1]; <br>
   oper3(t_vec,fin_mat,rc,generations,count_gener,B); <br>
   for(i=1; i<=x2_vari; i++) <br>
     new_genera[first][i+1] = t_vec[i]; <br>
 } <br>
      for(j=1; j<=(int) P4/2; j++) <br>
 { <br>
 { <br>
   first = find_parent(live,pop_size); <br>
   second = find_parent(live,pop_size); <br>
   new_genera[first][1] = 1.0; <br>
   new_genera[second][1] = 1.0; <br>
   for(i=1; i<=x2_vari; i++) <br>
     { <br>
       temp[1][i] = new_genera[first][i+1]; <br>
       temp[2][i] = new_genera[second][i+1]; <br>
     } <br>
   oper4(temp[1],temp[2],A,x2_vari); <br>
   for(i=1; i<=x2_vari; i++) <br>
     { <br>
       new_genera[first][i+1] = temp[1][i]; <br>