C语言实现粒子群算法优化二次函数详解

/*用粒子群算法求二次函数的最小值*/
/*四个粒子*/
# include <stdio.h>
#include<iostream>

/****************************************************************************************************************
基本参数的设定与初始化：（1）惯性权重w（2）最大速度v_max（3）学习因子c_1，c_2（4）随机数_rand_（5）群规模M
*****************************************************************************************************************/
float w=0.8;//惯性权重。当V_max很小时（V_max<2）:w=1;当V_max不很小时（V_max>3）:w=0.8;当V_max不确定时:w=0.8;
float v_max=0.1;//设置最大速度（疑问：如何估计最大速度）
int c_1=2,c_2=2;//学习因子：一般取2（也可以取其他值，但要保证范围在0～4）
float _rand_;//0～1间的随机数
const int M=4;


/****************************************************************************************************************
基本变量的设定与初始化：（1）各粒子初速度V[i]（2）各粒子当前的最优位置pbest[i]（3）集群最优位置gbest
*****************************************************************************************************************/
float X[M]={0};//定义各粒子的位置（这里有4个粒子）
float Y[M]={0};//定义各粒子各位置对应的目标值（这里有4个粒子）
float V[M]={0};//定义各粒子的速度。（疑问：如何估计各速度）
float pbest[M]={0};//各粒子当前的最优位置
float gbest;//集群最优位置


/****************************************************************************************************************
随机给出各粒子初始位置及其速度
*****************************************************************************************************************/
float X_shuiji[M]={0.5,2.6,0.2,1.6};//随机给出各粒子初始位置
float V_shuiji[M]={0.01,0.02,0.01,0.02};//随机给出各粒子初始速度