rtcgpsfda3.cpp

基于GPS的多目标优化及动态多目标优化源代码 里面 包含MOP, DMOP的程序

// RTCGPSFDA3.cpp: implementation of the RTCGPSFDA3 class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "RTDMOEA.h"
#include "RTCGPSFDA3.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

RTCGPSFDA3::RTCGPSFDA3()
{

}

RTCGPSFDA3::~RTCGPSFDA3()
{

}


void RTCGPSFDA3::InitPop()
{
   CSolution t_temp;
   t_temp.Init();
   for (int i=0;i<PopSize;i++)
   {
	   /**  设置 FDA3 的初始化GENE   **/
	  for (int j =0;j< 5 ;j++)
	   {
		   t_temp.Gene[j] =CRandom::Random(0.0,1.0);   
	   }
	   for ( j =5;j< CSolution::GetGeneLength() ;j++)
	   {
		   t_temp.Gene[j] =CRandom::Random(-1.0,1.0);   
	   }      /***  增加个体 */
	   m_Pop[i].Init();
       m_Pop[i].Copy(&t_temp);
   }
      InitPopType();

}

void RTCGPSFDA3::ComputeFitness(CSolution &m_1)
{
     EvaluationNum++;

	  for (int j =0;j< 5 ;j++)
	   {
	   		   if ((m_1.Gene[j]<0) ||(m_1.Gene[j]>1))
		   m_1.Gene[j] =CRandom::Random(0.0,1.0);   
	   }
	   for ( j =5;j< CSolution::GetGeneLength() ;j++)
	   {
	   		   if ((m_1.Gene[j]<-1.0) ||(m_1.Gene[j]>1))
		   m_1.Gene[j] =CRandom::Random(-1.0,1.0);   
	   }
	   
	   int nr=10;
	   int rt=ChangeGenerationNum;
	   double t= (1.0/(double)nr)*((int)CurrentGeneration/rt);

	   //计算Ft

	   double Ft=pow(10,2*sin(0.5*PI*t));
	   //计算gt
 	double gt=fabs(sin(0.5* PI * t));
      //计算 f1
	   m_1.Objective[0]=0;

	for (int i=0;i< 5;i++)
	{
		m_1.Objective[0] +=pow(m_1.Gene[i],Ft);
	}
	 m_1.Objective[0]/=5;
	 
	 double g=0;
	 for (i=5;i< CSolution::GetGeneLength();i++)
	 {
	    g+= (m_1.Gene[i]- gt) *(m_1.Gene[i]- gt);
	 }
	 g += 1 + gt;
	 
     double h = 1 - sqrt(m_1.Objective[0]/g);

	 m_1.Objective[1]= g*h;

}

void RTCGPSFDA3::RandomMutate(int Index)
{
	   /**  设置 FDA3 的初始化GENE   **/
   	  for (int j =0;j< 5 ;j++)
	   {
		   m_Pop[Index].Gene[j] =CRandom::Random(0.0,1.0);   
	   }
	   for ( j =5;j< CSolution::GetGeneLength() ;j++)
	   {
		   m_Pop[Index].Gene[j] =CRandom::Random(-1.0,1.0);   
	   }

}