?? mypso.cpp
字號(hào):
//PSO源文件myPSO.CPP如下:
//Header PSO.H
//Purpose Provide a class for particle swarm optimization
//Platform Visual C++ 6.0
//Date 2005.4.2
//Author ZhouXiandong
//群微粒算法：本算法涌群微粒算法求目標(biāo)函數(shù)的最小值
//本算法使用步驟
// （1）派生自己的群微粒類，類中必須定義double GetFit(PARTICLE&)方法，用來(lái)計(jì)算每個(gè)微粒的適合度
// （2）生成派生類實(shí)例，并在構(gòu)造函數(shù)中指明微粒坐標(biāo)維數(shù)和群體個(gè)數(shù)
// （2）設(shè)置微粒坐標(biāo)上界數(shù)組和下界數(shù)組，并用SetXup與SetXdown設(shè)置微粒坐標(biāo)上下界
// （3）用SetVmax方法設(shè)置微粒最大速度
// （4）設(shè)置可選參數(shù)：C1，C2，W和通訊函數(shù)
// （5）采用Run方法進(jìn)行優(yōu)化運(yùn)算，優(yōu)化后用GetBest方法獲得最優(yōu)個(gè)體適合度和坐標(biāo)
#include <stdlib.h>
#include <time.h>
#include <math.h>
#include "PSO.H"
//微粒構(gòu)造函數(shù)
PARTICLE::PARTICLE() //空構(gòu)造函數(shù)
{
X = 0; V = 0; XBest = 0; Dim = 0;
}
PARTICLE::PARTICLE(int n) //維數(shù)為參數(shù)的構(gòu)造函數(shù)
{
Dim = n;
X = new double[Dim];
V = new double[Dim];
XBest = new double[Dim];
}
//微粒析構(gòu)函數(shù)
PARTICLE::~PARTICLE()
{
if(Dim)
{
delete []X;
delete []V;
delete XBest;
}
}
//設(shè)置微粒的維數(shù)
void PARTICLE::SetDim(int d)
{
if(X) delete []X;
if(V) delete []V;
if(XBest) delete []XBest;
Dim = d;
X = new double[Dim];
V = new double[Dim];
XBest = new double[Dim];
}
//PSO構(gòu)造函數(shù)
PSO::PSO()
{
Particle = 0;
PNum = 0;
GBestIndex = 0;
Xup = 0;
Xdown = 0;
W = 1;
C1 = 2;
C2 = 2;
Com = 0;
}
PSO::PSO(int dim, int num)
{
Particle = new PARTICLE[num];
for(int i=0; i<num; i++) Particle[i].SetDim(dim);
PNum = num;
GBestIndex = 0;
Xup = new double[dim];
Xdown = new double[dim];
Vmax = new double[dim];
W = 1;
C1 = 2;
C2 = 2;
Com = 0;
}
//析構(gòu)函數(shù)
PSO::~PSO()
{
if(Particle) delete []Particle;
if(Xup) delete []Xup;
if(Xdown) delete []Xdown;
if(Vmax) delete []Vmax;
}
//設(shè)置坐標(biāo)上界
void PSO::SetXup(double *up)
{
if(!Particle) return;
for(int i=0; i<Particle[0].Dim; i++)
Xup[i] = up[i];
}
//設(shè)置坐標(biāo)下界
void PSO::SetXdown(double *d)
{
if(!Particle) return;
for(int i=0; i<Particle[0].Dim;i++)
Xdown[i] = d[i];
}
//設(shè)置最大速度
void PSO::SetVmax(double *max)
{
if(!Particle) return;
for(int i=0; i<Particle[0].Dim;i++)
Vmax[i] = max[i];
}
void PSO::SetVmax(double p)
{
if(!Particle) return;
for(int i=0; i<Particle[0].Dim; i++)
Vmax[i] = (Xup[i]-Xdown[i])*p;
}
//初始化群體
void PSO::Initialize()
{
if(!Particle) return;
static int kk=(unsigned)time(NULL);
srand((unsigned)time(NULL)+kk++);
GBestIndex = 0;
for(int i=0; i<PNum; i++)
{
for(int j=0; j<Particle[i].Dim; j++)
{
Particle[i].X[j] = rand()/(double)RAND_MAX*(Xup[j]-Xdown[j])+Xdown[j];//初始化坐標(biāo)
Particle[i].XBest[j] = Particle[i].X[j];
Particle[i].V[j] = rand()/(double)RAND_MAX*2*Vmax[j]-Vmax[j];//初始化速度
}
Particle[i].Fit = GetFit(Particle[i]); //計(jì)算該微粒適合度
Particle[i].FitBest = Particle[i].Fit; //設(shè)最優(yōu)適合度初值
if(Particle[i].Fit>Particle[GBestIndex].Fit) GBestIndex = i;//查找群體最優(yōu)微粒
}
}
//計(jì)算群體各個(gè)微粒適合度
void PSO::CalFit()
{
if(!Particle) return;
for(int i=0; i<PNum; i++)
Particle[i].Fit = GetFit(Particle[i]);
}
//微粒飛翔,產(chǎn)生新一代微粒
void PSO::ParticleFly()
{
static double FitBak[100];
if(!Particle) return;
static int tt=(unsigned)time(NULL);
srand((unsigned)time(NULL)+tt++);
//整個(gè)群體飛向新的位置
for(int i=0; i<PNum; i++)
{
for(int j=0; j<Particle[i].Dim; j++)
Particle[i].V[j] = W*Particle[i].V[j]+//(FitBak[i]-Particle[i].Fit)+//修改速度
rand()/(double)RAND_MAX*C1*(Particle[i].XBest[j]-Particle[i].X[j])+
rand()/(double)RAND_MAX*C2*(Particle[GBestIndex].XBest[j]-Particle[i].X[j]);
for(j=0; j<Particle[i].Dim; j++) //檢查速度最大值
{
if(Particle[i].V[j]>Vmax[j]) Particle[i].V[j] = Vmax[j];
if(Particle[i].V[j]<-Vmax[j]) Particle[i].V[j] = -Vmax[j];
}
for(j=0; j<Particle[i].Dim; j++)
{
Particle[i].X[j] += Particle[i].V[j]; //修改坐標(biāo)
if(Particle[i].X[j]>Xup[j]) Particle[i].X[j]=Xup[j];//保護(hù)
if(Particle[i].X[j]<Xdown[j]) Particle[i].X[j]=Xdown[j];
}
}
//計(jì)算各微粒適合度
CalFit();
for(i=0; i<PNum; i++) FitBak[i] = Particle[i].Fit;
//設(shè)置新的個(gè)體最好位置
for(i=0; i<PNum; i++)
if(Particle[i].Fit<=Particle[i].FitBest)
{
Particle[i].FitBest = Particle[i].Fit;
for(int j=0; j<Particle[i].Dim; j++)
Particle[i].XBest[j] = Particle[i].X[j];
}
//設(shè)置新的最優(yōu)個(gè)體
GBestIndex = 0;
for(i=0; i<PNum; i++)
if(Particle[i].FitBest<=Particle[GBestIndex].FitBest && i!=GBestIndex) GBestIndex = i;
}
//運(yùn)行群粒算法
PARTICLE& PSO::Run(int n)
{
Initialize();
double *opt_p = new double[Particle[0].Dim]; //通訊用數(shù)組，最優(yōu)點(diǎn)坐標(biāo)
double **opt_a = new double*[PNum]; //通訊用數(shù)組，所有點(diǎn)坐標(biāo)
for(int i=0; i<n; i++)
{ SetW(0.9 - i* (0.9 - 0.4) / n); //線性減小w的值;
ParticleFly();
if(Com) //通訊函數(shù)存在，完成通訊
{
for(int k=0; k<Particle[0].Dim; k++) opt_p[k]=Particle[GBestIndex].XBest[k];//拷貝最優(yōu)點(diǎn)坐標(biāo)
for(k=0; k<PNum; k++) opt_a[k] = Particle[k].X; //指向所有點(diǎn)坐標(biāo)
if(!Com(Particle[GBestIndex].FitBest,opt_p,opt_a,GBestIndex)) break;
}
}
delete []opt_p;
delete []opt_a;
return Particle[GBestIndex];
}
PARTICLE& PSO::Run(double fit)
{
double *opt_p = new double[Particle[0].Dim]; //通訊用數(shù)組，最優(yōu)點(diǎn)坐標(biāo)
double **opt_a = new double*[PNum]; //通訊用數(shù)組，所有點(diǎn)坐標(biāo)
Initialize();
do
{
ParticleFly();
if(Com) //通訊函數(shù)存在，完成通訊
{
for(int k=0; k<Particle[0].Dim; k++) opt_p[k]=Particle[GBestIndex].XBest[k];//拷貝最優(yōu)點(diǎn)坐標(biāo)
for(k=0; k<PNum; k++) opt_a[k] = Particle[k].X; //指向所有點(diǎn)坐標(biāo)
if(!Com(Particle[GBestIndex].FitBest,opt_p,opt_a,GBestIndex)) break;
}
}while(Particle[GBestIndex].FitBest<fit);
delete []opt_p;
delete []opt_a;
return Particle[GBestIndex];
}
//返回最佳個(gè)體
double PSO::GetBest(double *r)
{
for(int i=0; i<Particle[GBestIndex].Dim; i++)
r[i] = Particle[GBestIndex].XBest[i];
return Particle[GBestIndex].FitBest;
}
// pso.cpp : Defines the entry point for the console application.
//
//#include "stdafx.h"
#include "PSO.h"
#include <iostream.h>
#include <math.h>
//派生自己的PSO類
class MyPSO : public PSO
{
public:
MyPSO(int d, int n):PSO(d, n){}; //構(gòu)造函數(shù)，給出微粒維數(shù)和微粒個(gè)數(shù)
double GetFit(PARTICLE &p) //適合度計(jì)算方法，必須定義
{
double t=10;
double X[]={5,15,25,35,45,55,65,75,85,95};
double *Y=new double[10];
double xc[]={25,64,34,77,94,65,15,93};
double yc[]={87,28,53,92,55,24,11,14};
double p1[]={1.5 ,2,3,2,3,1.5,3,3};
double q[]={1.5, 3.5, 2, 2, 3, 1.5 ,2, 3};
double u[]={30,65,21,6,11,13,11,13};
double v[]={20,12,13,25,14,9,13,11};
double h[]={57,19,-32,38,24,-22,-14,-17};
double *k=new double [10];
int i,j,ii;
Y[0]=92;Y[9]=9;
for(i=1;i<9;i++)
{
Y[i]=p.X[i-1];
}
double s=0;
double* xt=new double[11];
double* yt=new double[11];
double* z=new double[11];
for(i=0;i<9;i++)
k[i]=(Y[i+1]-Y[i])/(X[i+1]-X[i]);//由已經(jīng)知道的XY求各段斜率
//下面求每小段的10等份后的各點(diǎn)的X，Y的值及對(duì)應(yīng)的Z值
s=0;
for(i=0;i<9;i++)
{
for(j=0;j<11;j++)
{
xt[j]=X[i]+j*0.1*t;
yt[j]=k[i]*(xt[j]-X[i])+Y[i];
z[j]=0;
for(ii=0;ii<8;ii++)
z[j]=z[j]+h[ii]*exp(-pow(fabs(xt[j]-xc[ii])/u[ii],p1[ii])-pow(fabs(yt[j]-yc[ii])/v[ii],q[ii]));
}
for(j=0;j<10;j++)
s=s+sqrt((0.1*t)*(0.1*t)+(k[i]*(0.1*t))*(k[i]*(0.1*t))+(z[j+1]-z[j])*(z[j+1]-z[j]));
}
delete []Y;
delete []k;
delete []xt;
delete []yt;
delete []z;
return s;
}
};
//定義通訊函數(shù)
bool MyCom(double fit, double *op, double**,int)
{
static long sn=1;
//if(sn%100==0)
//{
cout<<"\nNo="<<sn<<"\tFun="<<fit;
for(int i=0; i<8; i++)
cout<<"\tX("<<i<<")="<<op[i];
//}
sn++;
return true;
}
//申明相關(guān)數(shù)據(jù)
const int PNum = 30; //微粒個(gè)數(shù)
const int PDim = 8; //微粒維數(shù)
double Xup[] = {92,92,92,92,92,92,92,92}; //自變量上界
double Xdown[] = {9,9,9,9,9,9,9,9}; //自變量下界
//主程序
void main()
{
MyPSO pso(PDim, PNum); //生成微粒群實(shí)例
pso.SetXup(Xup); //設(shè)置自變量上界
pso.SetXdown(Xdown); //設(shè)置自變量下界
pso.SetVmax(0.2); //設(shè)置最大速度
pso.SetCom(MyCom); //設(shè)置通訊函數(shù)
cout<<"\nRun Now:\n";
pso.Run(1000); //運(yùn)行微粒群
cout<<"\nThe Result is:\t"<<pso.GetBest(Xup)<<"\n"; //輸出結(jié)果
}
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -