混沌粒子群算法（PSO）

翻到以前用c手撸的混沌粒子群算法

 

#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>

#define N 1000
#define D 2
#define MAXV 0.1
#define MINV -0.1
#define MAXP 10
#define MINP -10
#define ITERATION 10000
#define RANGE(a, min, max) ((((a) < (min))||((a)>(max)))?(((a)<(min))?(min):(max)):(a))

double c1 = 0.5;
double c2 = 0.5;
double w = 0.32;

double gBest[D];
double pBest[N][D];
double position[N][D];
double velocity[N][D];
double lastVelocity[N][D];

double f(double *a)
{
    //return 100 * pow((pow(*a, 2) - *(a + 1)), 2) + pow(*a-1, 2) + 100 * pow((pow(*(a + 1), 2) - *(a + 2)), 2) + pow(*(a + 1) - 1, 2);
    //return 0.5 + (pow(sin(sqrt(pow(*a, 2) + pow(*(a + 1), 2))),2) - 0.5) /(pow(1