C++完成簡略遺傳算法

C++完成簡略遺傳算法。本站提示廣大學習愛好者：（C++完成簡略遺傳算法）文章只能為提供參考，不一定能成為您想要的結果。以下是C++完成簡略遺傳算法正文

本文實例講述了C++完成簡略遺傳算法。分享給年夜家供年夜家參考。詳細完成辦法以下：

//遺傳算法 GA 
#include<iostream>
#include <cstdlib>
#include<bitset>
using namespace std;
const int L=5; //界說編碼的長度 
int f(int x) //界說測設函數f(x) 
{
int result;
result=x*x*x-60*x*x+900*x+100;
return result;
}
int main(int argc,char *argv[])
{
int a(0),b(32); //界說x的界說域規模
const int pop_size=8; //界說種群年夜小
// int L; //指定編碼的長度 
const int NG=20; //指定種群最年夜的滋生的代數 
int t=0; //以後滋生的代數 
int p[pop_size]; //界說種群 
int q[pop_size]; //界說滋生種群 即種群的下一代 
srand(6553); //界說隨機數生成的種子 
double sum; //適值總和 
double avl_sum; //過度均勻值 
double p_probability[pop_size]; //適值幾率 
double pp[pop_size];
double pro; //界說隨機生成的幾率 
float pc=0.90; //界說穿插的幾率 
float pm=0.05; //界說變異的幾率 
cout<<"初始的種群 "; 
for(int i=0;i<pop_size;i++) //生成初始的第0代種群 
  {
p[i]=rand()%31;
cout<<p[i]<<" ";
  }
   cout<<endl;
   cout<<endl;
   void Xover(int &,int &); //聲明穿插函數 
//當停滯原則不知足 即滋生代數沒到最年夜代數 ,持續滋生
while(t<=NG)             
{ 
cout<<"滋生的代數：t="<<t<<endl;
sum=0.0;
for(int i=0;i<pop_size;i++)      
  {
q[i]=p[i];
cout<<q[i]<<" ";
  }
  cout<<endl;
 for(int i=0;i<pop_size;i++) //盤算sum 
   sum +=f(p[i]);
 avl_sum=sum/pop_size;
 cout<<"sum="<<sum<<endl;
 cout<<"過度均勻值="<<avl_sum<<endl; 
    for(int i=0;i<pop_size;i++) //盤算適值幾率 
    {
      p_probability[i]=f(p[i])/sum;
if(i==0)
{
pp[i]=p_probability[i];
cout<<"pp"<<i<<"="<<pp[i]<<endl;
}
      else
      {
       pp[i]=p_probability[i]+pp[i-1];
    cout<<"pp"<<i<<"="<<pp[i]<<endl;
      }
//cout<<"p_probability"<<i<<"="<<p_probability[i]<<endl;
    }
    //選擇雙親
     for(int i=0;i<pop_size;i++) 
     {
     pro=rand()%1000/1000.0;
if(pro>=pp[0]&&pro<pp[1])
  p[i]=q[0]; 
else if(pro>=pp[1]&&pro<pp[2])
       p[i]=q[1];
     else if(pro>=pp[2]&&pro<pp[3])
       p[i]=q[2];
     else if(pro>=pp[3]&&pro<pp[4])
       p[i]=q[3];
     else if(pro>=pp[4]&&pro<pp[5])
       p[i]=q[4];
     else 
       p[i]=q[5]; 
     }
//雜交算子
int r=0;
int z=0; 
for(int j=0;j<pop_size;j++) 
{
  pro=rand()%1000/1000.0;
if(pro<pc)
{
 ++z;
 if(z%2==0)
  Xover(p[r],p[j]);
 else
  r=j; 
} 
} 
//變異算子 
for(int i=1;i<=pop_size;i++)
 for(int j=0;j<L;j++)
{
 pro=rand()%1000/1000.0; //在【0,1】區間發生隨機數
if(pro<pm)
{
bitset<L>v(p[i]);      
v.flip(j);
p[i]=v.to_ulong();
} 
 } 
t++;
cout<<endl; //種群滋生一代 
 }
 cout<<"終究成果：";  
 for(int i(0);i<pop_size;i++) //算法停止，輸入成果 
  {
 cout<<p[i]<<" ";
  }
  cout<<endl;
return 0;
}
//界說雜交操作 
 void Xover(int &a,int &b)         
{ 
int pos; //隨機生成雜交點 即第幾個重量停止互相交流
pos=rand()%5+1; //在n個重量中，隨機肯定第pos個重量 
int j,k;
  j=pos;
  k=pos;
bitset<L>e(a);
bitset<L>f(b); //前pos個重量停止互相交流
bitset<L>g;            
bitset<L>h;
for(int i=0;i<pos;i++)
{
if(e[i]==1)
 g.set(i); 
   }
  for(int i=0;i<pos;i++)
   {
   if(f[i]==1)
    h.set(i);
   }
  for(j;j<L;j++)
   {
   if(f[j]==1)
    g.set(j);
   }
  for(k;k<L;k++)
   {
   if(e[k]==1)
    h.set(k);
   }
a=g.to_ulong();
b=h.to_ulong();  
}

願望本文所述對年夜家的C++法式設計有所贊助。