mopso.cpp~

此程序为多目标粒子群程序

//mopso con epsilon dominance y clusters

#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <vector>
#include <math.h>
#include "efile.h"
using namespace std;



double h(double _f1,double _f2,double _g){
  double _sumatoria(0.0),_PI(3.1415927);
  _sumatoria=(_f1/(1+_g)*(1+sin(3*_PI*_f1)))+(_f2/(1+_g)*(1+sin(3*_PI*_f2)));
  return (3-_sumatoria);
}
  
double g(vector<double> &_x)
{
    
  double _suma(0.0),_PI(3.1415927);
    
  for(int _i=2;_i<22;_i++)
    _suma+=_x[_i];
    
  return 1.0+9.0/20.0*_suma;
}

MOPSO::~MOPSO(){

}

MOPSO::MOPSO(int _dims, int _objs,int _parts,int _nclusters,int _gmax,double *_lb,double *_ub ){

  srand(time(0));
  ndimensions=_dims;
  nobjectives=_objs;
  nparticles=_parts;
  nclusters=_nclusters;
  gmax=_gmax;
  gen=0;
  lb.resize(ndimensions);
  ub.resize(ndimensions);
  particles.resize(nparticles);

  for(int _j(0);_j<ndimensions;_j++){
    lb[_j]=_lb[_j];
    ub[_j]=_ub[_j];
  }
  archive.init(ndimensions, nobjectives, nclusters,100);
}

void MOPSO::initialize(){
  C1=2.0;
  C2=2.0;
  W=0.4;

  for(int _i(0);_i<nparticles;_i++){
    for(int _j(0);_j<ndimensions;_j++){
      particles[_i].vel[_j]=0;
      particles[_i].x[_j]=rnd(lb[_j],ub[_j]);
      particles[_i].xpbest[_j]=particles[_i].x[_j];
    }

    function(_i);

    for(int _j(0);_j<nobjectives;_j++){
      particles[_i].fxpbest[_j]=particles[_i].fx[_j];  
    }

    archive.add(particles[_i],-1);
  }
}



void MOPSO::execute(){
  initialize();
  for(gen=0;gen<gmax;gen++){
    archive.updatematrix=true;
    archive.hierarchicalClustering();
    flight();
  }
  archive.output();
}  

void MOPSO::flight(){
  for(int _i(0);_i<nparticles;_i++){
    //        W=rnd(0.1,0.5);
    //C2=rnd(1.5,2.5);
    //C1=rnd(1.5,2.5);
    
    
    int _whichcluster=(int)_i/(nparticles/nclusters);
    int _gbestselected;
    _gbestselected=archive.selectClusteredRandomSolution(_whichcluster);
    for(int _k(0);_k<5;_k++){
      for(int _j(0);_j<ndimensions;_j++){
particles[_i].vel[_j]=W*particles[_i].vel[_j]+C1*rnd(0,1)*(archive.solutions[_gbestselected].x[_j]-particles[_i].x[_j])+C2*rnd(0,1)*(particles[_i].xpbest[_j]-particles[_i].x[_j]);
//particles[_i].vel[_j]=W*particles[_i].vel[_j]+C1*rnd(0,1)*(archive.solutions[_gbestselected].x[_j]-particles[_i].x[_j]);//+C2*rnd(0,1)*(particles[_i].xpbest[_j]-particles[_i].x[_j]);
	//perturbation(_i);
	
	particles[_i].x[_j]+=particles[_i].vel[_j];
	if(particles[_i].x[_j]<lb[_j])
	  particles[_i].x[_j]=lb[_j];
	if(particles[_i].x[_j]>ub[_j])
	  particles[_i].x[_j]=ub[_j];
      }

      function(_i);

      int _tmp=archive.domine1(particles[_i].fx,particles[_i].fxpbest);
      if(_tmp==11||_tmp==1){
	  copy(particles[_i].fxpbest,particles[_i].fx);
	  copy(particles[_i].xpbest,particles[_i].x);
	  archive.add(particles[_i],(int)_i/(nparticles/nclusters));
      }

    }
	
  }
}

void MOPSO::perturbation(int _whichparticle){
  int _dimension=0;
  double _lb,_ub,_rango;
  double _gt=(double)gen/gmax;

  double  _pM=pow(_gt,1.7)-2.0*_gt+1.0;
  int _flag(0);
  if((rnd(0.0,1.0)>_pM)&&_flag<=ndimensions){

      _dimension= (int) rnd(0,ndimensions);
      _rango=(ub[_dimension]-lb[_dimension])*_pM/2.0;//totGen

      if(particles[_whichparticle].x[_dimension]-_rango<lb[_dimension]) 
	_lb=lb[_dimension]; 
      else 
	_lb=particles[_whichparticle].x[_dimension]-_rango;

      if(particles[_whichparticle].x[_dimension]+_rango>ub[_dimension]) 
	_ub=ub[_dimension]; 
      else 
	_ub=particles[_whichparticle].x[_dimension]+_rango;

	particles[_whichparticle].x[_dimension]=rnd(_lb,_ub);
	_flag++;

    }

}



void MOPSO::function(int _w){

  /*
  double sum1=0.0,sum2=0.0;
    
  for(int i=0;i<ndimensions;i++){
    if(i<ndimensions-1)
      sum1+=-10.0*exp((-0.2)*sqrt( pow(particles[_w].x[i],2.0)+pow(particles[_w].x[i+1],2.0)));
    sum2+=pow(fabs(particles[_w].x[i]),0.8)+5*sin(pow(particles[_w].x[i],3.0));
  }
  particles[_w].fx[0]=sum1;
  particles[_w].fx[1]=sum2;    
  //*/
  //     }
  /*
    double f;
    double g;
    double h;
    double x1=particles[_w].x[0];
    double x2=particles[_w].x[1];

    f=x1;
    g=11+pow(x2,2)-10*cos(2*3.1415926*x2);
    if(f<=g)
    h=1-pow((f/g),(double)(1.0/2.0));
    else
    h=0;
     particles[_w].fx[0]=f;
    particles[_w].fx[1]=g*h;
  
    //    */
  /*dtlz6
    double _x[3];
    double PI=3.1415927;
    _x[0]=particles[_w].x[0];
    _x[1]=particles[_w].x[1];
    
    particles[_w].fx[0]=_x[0];
    particles[_w].fx[1]=_x[1];
    double _g=g(particles[_w].x);
    particles[_w].fx[2]=(1+_g)*h(_x[0],_x[1],_g);
  */
  //*zdt1

  double g(0),h(0);
  particles[_w].fx[0]=particles[_w].x[0];
  for(int _i(1);_i<ndimensions;_i++)
    g+=(double)particles[_w].x[_i]/((double)ndimensions-1.0);
  g*=9.0;
  g+=1.0;
  h=1.0-sqrt((double)particles[_w].fx[0]/g);
  particles[_w].fx[1]=g*h;
  //*/

    }


double MOPSO::rnd(double _min,double _max){
  return((double)(_min + ((double)(_max-_min)*rand()/(double)(RAND_MAX+_min))));
}


void MOPSO::copyx(vector<double> &_a,vector<double> &_b){
  for(int _i(0);_i<ndimensions;_i++){
    _a[_i]=_b[_i];
  }
}

void MOPSO::copyfx(vector<double> &_a,vector<double> &_b){
  for(int _i(0);_i<nobjectives;_i++){
    _a[_i]=_b[_i];
  }
}


void MOPSO::copy(vector<double> &_a,vector<double> &_b){
  for(int _i(0);_i<_a.size();_i++){
    _a[_i]=_b[_i];
  }

}






int main()
{
  /*

  int _ndims=3;
  int _nobjs=2;
  int _parts=40;
  int _nclusters=4;
  int _gmax=40;
  double *_lbound=new(double[_ndims]);//=0.1;
  double *_ubound=new(double[_ndims]);//=1;
  _lbound[0]=-5.0;
  _lbound[1]=-5.0;
  _lbound[2]=-5.0;
  _ubound[0]=5.0;
  _ubound[1]=5.0;
  _ubound[2]=5.0;
  //*/


  /*//deb
  int _ndims=2;
  int _nobjs=2;
  int _parts=40;
  int _nclusters=4;
  int _gmax=40;
  double *_lbound=new(double[_ndims]);//=0.1;
  double *_ubound=new(double[_ndims]);//=1;
    _lbound[0]=0.0;
    _lbound[1]=-5.0;
    _ubound[0]=1.0;
    _ubound[1]=5.0;
    //*/

  /*//dtlz6
  int _ndims=22;
  int _nobjs=3;
  int _parts=40;
  int _nclusters=4;
  int _gmax=100;
  double *_lbound=new(double[_ndims]);//=0.1;
  double *_ubound=new(double[_ndims]);//=1;
  for(int _i(0);_i<_ndims;_i++){
    _lbound[_i]=0.0;
    _ubound[_i]=1.0;
  }
    //*/
  //*//zdt1
  int _ndims=30;
  int _nobjs=2;
  int _parts=40;
  int _nclusters=4;
  int _gmax=100;
  double *_lbound=new(double[_ndims]);//=0.1;
  double *_ubound=new(double[_ndims]);//=1;
  for(int _i(0);_i<_ndims;_i++){
    _lbound[_i]=0.0;
    _ubound[_i]=1.0;
  }
    //*/

  MOPSO  pso(_ndims,_nobjs,_parts,_nclusters , _gmax,_lbound,_ubound);
  pso.execute();

  delete [] _lbound;
  delete [] _ubound;
}