#include"stdio.h"
#include"stdlib.h"
#include"time.h"
#include"math.h"
const int NUM=40;//粒子数
const int DIM=30;//维数
const double c1=1.8;//参数
const double c2=1.8;//参数
double xmin=-100.0;//位置下限
double xmax=100.0;//位置上限
double gbestx[DIM];//全局最优位置
double gbestf;//全局最优适应度
struct particle {//定义一个粒子
doublex[DIM];//当前位置矢量
doublebestx[DIM];//历史最优位置
doublef;//当前适应度
doublebestf;//历史最优适应度
}swarm[NUM];//定义粒子群
#define randf((rand()%10000+rand()%10000*10000)/100000000.0) //产生-1随机浮点数
double f1(double x[]) {//测试函数：超球函数
floatz=0;
for(inti=0;i
z+=(x[i])*(x[i]);
returnz;
}
void main() {
for(int i=0;i
gbestx[i]=randf*(xmax-xmin)+xmin;
gbestf=[**************].0;
for(int i=0;i
particle* p1=&swarm[i];
for(int j=0; j
p1->x[j]=randf*(xmax-xmin)+xmin;
p1->f=f1(p1->x);
p1->bestf=[**************].0;
}
for(int t=0;t
for(int i=0; i
particle* p1=&swarm[i];
for(int j=0; j
p1->x[j]+=c1*randf*(p1->bestx[j]-p1->x[j])
+c2*randf*(gbestx[j]-p1->x[j]);
p1->f=f1(p1->x);
if(p1->fbestf){//改变历史最优
for(int j=0;j
p1->bestx[j]=p1->x[j];
p1->bestf=p1->f;
}
if(p1->f
for(int j=0;j
gbestx[j]=p1->x[j];
for(int j=0; j
p1->x[j]=randf*(xmax-xmin)+xmin;
gbestf=p1->f;
}
}
}
printf("%g\n", gbestf);
}
#include"stdio.h"
#include"stdlib.h"
#include"time.h"
#include"math.h"
const int NUM=40;//粒子数
const int DIM=30;//维数
const double c1=1.8;//参数
const double c2=1.8;//参数
double xmin=-100.0;//位置下限
double xmax=100.0;//位置上限
double gbestx[DIM];//全局最优位置
double gbestf;//全局最优适应度
struct particle {//定义一个粒子
doublex[DIM];//当前位置矢量
doublebestx[DIM];//历史最优位置
doublef;//当前适应度
doublebestf;//历史最优适应度
}swarm[NUM];//定义粒子群
#define randf((rand()%10000+rand()%10000*10000)/100000000.0) //产生-1随机浮点数
double f1(double x[]) {//测试函数：超球函数
floatz=0;
for(inti=0;i
z+=(x[i])*(x[i]);
returnz;
}
void main() {
for(int i=0;i
gbestx[i]=randf*(xmax-xmin)+xmin;
gbestf=[**************].0;
for(int i=0;i
particle* p1=&swarm[i];
for(int j=0; j
p1->x[j]=randf*(xmax-xmin)+xmin;
p1->f=f1(p1->x);
p1->bestf=[**************].0;
}
for(int t=0;t
for(int i=0; i
particle* p1=&swarm[i];
for(int j=0; j
p1->x[j]+=c1*randf*(p1->bestx[j]-p1->x[j])
+c2*randf*(gbestx[j]-p1->x[j]);
p1->f=f1(p1->x);
if(p1->fbestf){//改变历史最优
for(int j=0;j
p1->bestx[j]=p1->x[j];
p1->bestf=p1->f;
}
if(p1->f
for(int j=0;j
gbestx[j]=p1->x[j];
for(int j=0; j
p1->x[j]=randf*(xmax-xmin)+xmin;
gbestf=p1->f;
}
}
}
printf("%g\n", gbestf);
}