修改之后的源程序！

// CMVSOGA.h : main header file for the CMVSOGA.cpp
////////////////////////////////////////////////////////////////////
/////                                                          /////
/////                沈陽航空工業(yè)學(xué)院 動力工程系               /////
/////                       作者：李立新                       /////
/////                   完成日期：2006.08.02                   /////
/////                   修改日期：2007.04.10                                       /////
////////////////////////////////////////////////////////////////////

#if !defined(AFX_CMVSOGA_H__45BECA_61EB_4A0E_9746_9A94D1CCF767__INCLUDED_)
#define AFX_CMVSOGA_H__45BECA_61EB_4A0E_9746_9A94D1CCF767__INCLUDED_

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000
#include "Afxtempl.h"
#define variablenum 14
class CMVSOGA
{
public:
 CMVSOGA();
 ~CMVSOGA();
 void selectionoperator();
 void crossoveroperator();
 void mutationoperator();
 void initialpopulation(int, int ,double ,double,double *,double *);           //種群初始化
 void generatenextpopulation();          //生成下一代種群
 void evaluatepopulation();           //評價個體，求最佳個體
 void calculateobjectvalue();          //計算目標函數(shù)值
 void calculatefitnessvalue();          //計算適應(yīng)度函數(shù)值
 void findbestandworstindividual();         //尋找最佳個體和最差個體
 void performevolution();   
 void GetResult(double *);
 void GetPopData(CList <double,double>&);
 void SetFitnessData(CList <double,double>&,CList <double,double>&,CList <double,double>&);
private:
 struct individual
 {
  double chromosome[variablenum];         //染色體編碼長度應(yīng)該為變量的個數(shù)
  double value;         
  double fitness;             //適應(yīng)度
 };
 double variabletop[variablenum];         //變量值
 double variablebottom[variablenum];         //變量值
 int popsize;              //種群大小
// int generation;              //世代數(shù)
 int best_index;  
 int worst_index;
 double crossoverrate;            //交叉率
 double mutationrate;            //變異率
 int maxgeneration;             //最大世代數(shù)
 struct individual bestindividual;         //最佳個體
 struct individual worstindividual;         //最差個體
 struct individual current;              //當(dāng)前個體
 struct individual current1;              //當(dāng)前個體
 struct individual currentbest;          //當(dāng)前最佳個體
 CList <struct individual,struct individual &> population;   //種群
 CList <struct individual,struct individual &> newpopulation;  //新種群
 CList <double,double> cfitness;          //存儲適應(yīng)度值
 //怎樣使鏈表的數(shù)據(jù)是一個結(jié)構(gòu)體????主要是想把種群作成鏈表。節(jié)省空間。
};
#endif

 

執(zhí)行文件：

// CMVSOGA.cpp : implementation file
//

#include "stdafx.h"
//#include "vld.h"
#include "CMVSOGA.h"
#include "math.h"
#include "stdlib.h"


#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif
/////////////////////////////////////////////////////////////////////////////
// CMVSOGA.cpp
CMVSOGA::CMVSOGA()
{
 best_index=0;  
 worst_index=0;
 crossoverrate=0;            //交叉率
 mutationrate=0;            //變異率
 maxgeneration=0;
}
CMVSOGA::~CMVSOGA()
{
 best_index=0;  
 worst_index=0;
 crossoverrate=0;            //交叉率
 mutationrate=0;            //變異率
 maxgeneration=0;
 population.RemoveAll();   //種群
 newpopulation.RemoveAll();  //新種群
 cfitness.RemoveAll(); 
}
void CMVSOGA::initialpopulation(int ps, int gen ,double cr ,double mr,double *xtop,double *xbottom)  //第一步，初始化。
{
 //應(yīng)該采用一定的策略來保證遺傳算法的初始化合理，采用產(chǎn)生正態(tài)分布隨機數(shù)初始化？選定中心點為多少？
 int i,j;
 popsize=ps;
 maxgeneration=gen;
 crossoverrate=cr;
 mutationrate =mr;
 for (i=0;i<variablenum;i++)
 {
  variabletop[i] =xtop[i];
  variablebottom[i] =xbottom[i];
 }
 //srand( (unsigned)time( NULL ) );  //尋找一個真正的隨機數(shù)生成函數(shù)。
 for(i=0;i<popsize;i++)
 { 
  for (j=0;j<variablenum ;j++)
  {
   current.chromosome[j]=double(rand()%10000)/10000*(variabletop[j]-variablebottom[j])+variablebottom[j];
  }
  current.fitness=0;
  current.value=0;
  population.InsertAfter(population.FindIndex(i),current);//除了初始化使用insertafter外,其他的用setat命令。
 }
}
void CMVSOGA::generatenextpopulation()//第三步，生成下一代。
{
 //srand( (unsigned)time( NULL ) );
 selectionoperator();
 crossoveroperator();
 mutationoperator();
}
//void CMVSOGA::evaluatepopulation()   //第二步，評價個體，求最佳個體
//{
// calculateobjectvalue();
// calculatefitnessvalue();   //在此步中因該按適應(yīng)度值進行排序.鏈表的排序.
// findbestandworstindividual();
//}
void CMVSOGA:: calculateobjectvalue()  //計算函數(shù)值，應(yīng)該由外部函數(shù)實現(xiàn)。主要因為目標函數(shù)很復(fù)雜。
{
 int i,j;
    double x[variablenum];
 for (i=0; i<popsize; i++)
 {
  current=population.GetAt(population.FindIndex(i));  
  current.value=0;
  //使用外部函數(shù)進行，在此只做結(jié)果的傳遞。
  for (j=0;j<variablenum;j++)
  {
   x[j]=current.chromosome[j];
   current.value=current.value+(j+1)*pow(x[j],4);
  }
  ////使用外部函數(shù)進行，在此只做結(jié)果的傳遞。
  population.SetAt(population.FindIndex(i),current);
 }
}

void CMVSOGA::mutationoperator()  //對于浮點數(shù)編碼，變異算子的選擇具有決定意義。
          //需要guass正態(tài)分布函數(shù)，生成方差為sigma，均值為浮點數(shù)編碼值c。
{
// srand((unsigned int) time (NULL));
 int i,j;
 double r1,r2,p,sigma;//sigma高斯變異參數(shù)
 
 for (i=0;i<popsize;i++)
 {
  current=population.GetAt(population.FindIndex(i));

  //生成均值為current.chromosome，方差為sigma的高斯分布數(shù)
  for(j=0; j<variablenum; j++)
  {   
   r1 = double(rand()%10001)/10000;
   r2 = double(rand()%10001)/10000;
   p = double(rand()%10000)/10000;
   if(p<mutationrate)
   {
    double sign;
    sign=rand()%2;
    sigma=0.01*(variabletop[j]-variablebottom [j]);
    //高斯變異
    if(sign)
    {
     current.chromosome[j] = (current.chromosome[j] 
      + sigma*sqrt(-2*log(r1)/0.4323)*sin(2*3.1415926*r2));
    }
    else
    {
     current.chromosome[j] = (current.chromosome[j] 
      - sigma*sqrt(-2*log(r1)/0.4323)*sin(2*3.1415926*r2));
    }
    if (current.chromosome[j]>variabletop[j])
    {
     current.chromosome[j]=double(rand()%10000)/10000*(variabletop[j]-variablebottom[j])+variablebottom[j];
    }
    if (current.chromosome[j]<variablebottom [j])
    {
     current.chromosome[j]=double(rand()%10000)/10000*(variabletop[j]-variablebottom[j])+variablebottom[j];
    }
   }
  }
  population.SetAt(population.FindIndex(i),current);
 }
}
void CMVSOGA::selectionoperator()   //從當(dāng)前個體中按概率選擇新種群,應(yīng)該加一個復(fù)制選擇,提高種群的平均適應(yīng)度
{
 int i,j,pindex=0;
 double p,pc,sum;
 i=0;
 j=0;
 pindex=0;
 p=0;
 pc=0;
 sum=0.001;
 newpopulation.RemoveAll();
 cfitness.RemoveAll();
  //鏈表排序
// population.SetAt (population.FindIndex(0),current); //多余代碼
 for (i=1;i<popsize;i++)
 { 
  current=population.GetAt(population.FindIndex(i));
  for(j=0;j<i;j++)   //從小到大用before排列。
  {
   current1=population.GetAt(population.FindIndex(j));//臨時借用變量
   if(current.fitness<=current1.fitness)  
   {
    population.InsertBefore(population.FindIndex(j),current);
    population.RemoveAt(population.FindIndex(i+1));
    break;
   }
  }
//  m=population.GetCount();
 }
 //鏈表排序
 for(i=0;i<popsize;i++)//求適應(yīng)度總值，以便歸一化,是已經(jīng)排序好的鏈。
 {
  current=population.GetAt(population.FindIndex(i)); //取出來的值出現(xiàn)問題.
  sum+=current.fitness;
 }
 for(i=0;i<popsize; i++)//歸一化
 {
  current=population.GetAt(population.FindIndex(i)); //population 有值,為什么取出來的不正確呢??