//粒子群PSO算法

#include
#include
#include
#include??
#include
#include

#define	PI?3.141592653589	
#define?P_num?40000??????????			//粒子數目
#define?dim?30					???????//維度
#define?N???10						??//?每條進程執行數目

#define?low?-100???????????			??//搜索域范圍
#define?high?100

#define?iter_num?1000
#define?V_max?20??????????		//速度范圍

#define?c1?1.5					//更新速度時的加速常量
#define?c2?1.5					//更新速度時的加速常量
#define?w?0.5					//更新速度時的權重
#define?alp?1					//更新位置時的權重
double?particle[P_num][dim];???????????//個體位置集合
double?particle_loc_best[P_num][dim];??//每個個體局部最優向量
double?particle_loc_fit[P_num];????????//個體的局部最優適應度有局部最優向量計算而來
double?particle_glo_best[dim];?????????//全局最優向量
double?gfit;???????????????????????????//全局最優適應度有全局最優向量計算而來
double?particle_v[P_num][dim];?????????//記錄每個個體的當前代速度向量
double?particle_fit[P_num];????????????//記錄每個粒子的當前代適應度

/*時間函數*/
double?mytime（）
{
????double?ts?=?0.0;
????struct?timeval?mt;
????gettimeofday（&mt（struct?timezone*）0）;
????ts?=?（double）（mt.tv_sec+mt.tv_usec*1.0e-6）;
????return?（ts）;
}

int?get_rand（int?comm_sz）{
??srand（（unsigned）time（NULL））;
??return?rand（）%（P_num/comm_sz）;
}

double?Sphere（double?a[]）
{

	int?i;
	double?sum=0.0;
	for（i=0;?i	{
		sum+=a[i]*a[i];
	}
	return?sum;
}
double?Rosenbrock（double?a[]）
{

	int?i;
	double?sum=0.0;
	for（i=0;i	{
????????sum+=?100*（a[i+1]-a[i]*a[i]）*（a[i+1]-a[i]*a[i]）+（a[i]-1）*（a[i]-1）;
	}
??	????return?sum;
}

double?Rastrigin（double?a[]）
{
	int?i;
	double?sum=0.0;
	for（i=0;i	{
		sum+=a[i]*a[i]-10.0*cos（2*PI*a[i]）+10.0;
	}
	return?sum;
}

double?fitness（double?a[]）?????????????//適應度函數
{
	return?Rastrigin（a）;
}

/*初始化粒子，包括初始化速度，位置，各自的適應度以及*/
void?initial（）
{
	int?ij;
	for（i=0;?i	{
		for（j=0;?j		{
			particle[i][j]?=?low+（high-low）*1.0*rand（）/RAND_MAX;????//初始化群體位置隨機
			particle_loc_best[i][j]?=?particle[i][j];???????????????//將當前最優結果寫入局部最優集合
			particle_v[i][j]?=?-V_max+2*V_max*1.0*rand（）/RAND_MAX;????//初始化速度
		}
	}
	
	for（i=0;?i	{
		particle_fit[i]?=?fitness（particle[i]）;
		particle_loc_fit[i]?=?particle_fit[i];			//當前局部適應度
	}
	
	gfit?=?particle_loc_fit[0];??????		//初始化全局最優適應度
	
	j=0;
	
	for（i=1;?i	{
		if（particle_loc_fit[i]		{
			gfit?=?particle_loc_fit[i];
			j?=?i;							//記錄下標
			
		}
	}
	
	for（i=0;?i	{
		particle_glo_best[i]?=?particle_loc_best[j][i];
	}
}

/*更新位置函數*/
void?renew_particle（int?aint?b）
{
	int?ij;
	for（i=a;?i	{
		for（j=0;?j		{
			particle[i][j]?+=??alp*particle_v[i][j];????
			if（particle[i][j]?>?high）
			{
				particle[i][j]?=?high;
			}
			if（particle[i][j]?			{
				particle[i][j]?=?low;
			}
		}
	}
}

/*更新粒子速度*/
void?renew_var（int?aint?bint?*m）
{
	int?i?j;
	for（i=a;i	{
		particle_fit[i]?=?fitness（particle[i]）;			//計算該粒子適應度
		
		if（particle_fit[i]?