遗传算法求解函数最大值(两参数)

}
float avg(int M){ //计算适应度平均值
int i,sum=0;
for(i=0;i<M;i++){
sum=sum+afit[i];
}
return sum/(M+0.0);
}
int son::jiaocha(int a[2000][10],int d[2000][10],float &pc,int M){ //交叉运算
//cout<<"jjl4:"<<endl;for(i=0;i<M;i++) cout<<jjl4[i]<<endl;
}
int next(int a[2000][10],int d[2000][10],int M){ //将父代的二值数组移植到子代
int i,j,l=0,k=M;
//cout<<"afit:"<<endl;for(i=0;i<M;i++) cout<<afit[i]<<endl;
//cout<<"jjl"<<endl;for(i=0;i<M;i++) cout<<jjl[i]<<endl;
//cout<<"jjl3"<<endl;for(i=0;i<M;i++) cout<<jjl3[i]<<endl;
}
else if(afit[i-1]>fit[j]){
afit[i]=fit[j];
jjl[i]=j;
jjl3[i]=j;
}
}
}
}
//cout<<"fit:"<<endl;for(i=0;i<M;i++) cout<<fit[i]<<endl;
int b[2000][10],c[2000][10];
for(i=0;i<M&&l<M;i++){
while(l>0&&(jjl2[l-1]==jjl2[l])){
k--;
l++;
}
for(j=0;j<10;j++){
b[i][j]=a[jjl2[l]][j];
}
}
printt(a,d,M);
}
void chushihua(){ //临时数组初始化
int i;
for(i=0;i<2000;i++){
jjl[i]=0;
jjl2[i]=0;
jjl3[i]=0;
jjl4[i]=0;
x[i]=0.0;
class father{ //定义父类
public:
void select(int M);
private:
int i,j,k;
};
class son:public father{ //定义子类
public:
int jiaocha(int a[2000][10],int d[2000][10],float &pm,int M);
//fit[i]=x[i]*sin(10*3.1416*x[i])+2.0;
//fit[i]=sin(x[i])/x[i];
//fit[i]=(x[i]+2.0);
//fit[i]=(2.0)*cifang(x[i],3)-cifang(x[i],2)-(1/2.0)*x[i]+1;
//fit[i]=x[i]*y[i];
//fit[i]=1.0/(100*(x[i]*x[i]-y[i])*(x[i]*x[i]-y[i])+(1-x[i])*(1-x[i]));
fit[i]=1-(x[i]-1)*(x[i]-1)-(y[i]-1)*(y[i]-1); //适应度
else pc=0.8; //自适应遗传算法概率设置
r=rand()/(RAND_MAX+0.0); //cout<<"r="<<r<<endl;
if(r<=pc){
k=10*rand()/(RAND_MAX+0.0); //cout<<"k="<<k<<endl;
}
}
for(i=0;i<M-1;i++){
for(j=i;j<M;j++){
if(jjl2[i]>jjl2[j]){
temp=jjl2[j];
jjl2[j]=jjl2[i];
jjl2[i]=temp;
}
}
}
for(i=0;i<M-1;i++){
for(j=i;j<M;j++){
if(jjl4[i]>jjl4[j]){
temp=jjl4[j];
jjl4[j]=jjl4[i];
jjl4[i]=temp;
}
}
}
//cout<<"jjl2:"<<endl;for(i=0;i<M;i++) cout<<jjl2[i]<<endl;
a[i][j]=0;
}
l++;
}
l=0;k=M;
for(i=0;i<M&&l<M;i++){
while(l>0&&(jjl4[l-1]==jjl4[l])){
k--;
l++;
}
for(j=0;j<10;j++){
c[i][j]=d[jjl4[l]][j];
int i,j;
srand((unsigned)time(0));
for(i=0;i<M;i++){
for(j=0;j<10;j++){
a[i][j]=rand01();
}
}
for(i=0;i<M;i++){
for(j=0;j<10;j++){
d[i][j]=rand01();
d[i][j]=0;
}
l++;
}
for(i=0;i<M;i++)
for(j=0;j<10;j++)
a[i][j]=0;
//cout<<"k="<<k<<endl;
for(i=0;i<k;i++)
for(j=0;j<10;j++)
a[i][j]=b[i][j];
//int i,j,k,l=0;
int l=0;
float r;
srand((unsigned)time(NULL));
for(j=1;j<M/2;j++){
if(afit[j]>avg(M)) pc=0.8-0.2*(afit[j]-avg(M))/(afit[0]-avg(M)+0.0);
}
cout<<endl;
}
cout<<endl;
for(i=0;i<M;i++){
for(j=0;j<10;j++){
cout<<d[i][j]<<" ";
}
cout<<endl;
}
}
void chushiM(int a[2000][10],int d[2000][10],int M){ //种群初始化
else r=1;
return r;
}
void printt(int a[2000][10],int d[2000][10],int M){ //显示数组
int i,j;
for(i=0;i<M;i++){
for(j=0;j<10;j++){
cout<<a[i][j]<<" ";
}
void sumfitp(int M){ //计算转轮法需要的个体累计概率
int i;
float sum=0.0;
for(i=0;i<M;i++){
sum=sum+fit[i];
}
tpfit[0]=fit[0]/sum;
for(i=1;i<M;i++){
int bianyi(int a[2000][10],int d[2000][10],float &pm,int M);
private:
int i,j,k;
};
int cifang(int a,int b){ //a的b次方