liziqun.zip_liziqun java_particle swarm

import java.util.*; public class liziqun1 { public static double pi=3.1415926; public static int num=60; public static int T=200; public static double func(double x) { double y=1-(Math.cos(3*x))*(Math.exp(-x)); return y; } public static void main(String args[]) { double ran; double c1=2.05,c2=2.05; double d1=0.2,d2=0.7; double ws=0.9,we=0.4; double x_down=0.0,x_up=4.0; double w; int m=1; double Vmax; double[] x= new double[num]; double[] v=new double[num]; double[] g=new double[num]; double[] pbest=new double[num]; double[] pbest_x=new double[num]; double gbest; double gbest_x; Vmax=(x_up-x_down)/16; for(int i=0;i<num;i++) { ran=Math.random(); x[i]=x_down+(x_up-x_down)*ran; v[i]=Vmax*(2*ran-1); pbest[i]=func(x[i]); pbest_x[i]=x[i]; } gbest=pbest[0]; gbest_x=x[0]; for(int i=0;i<num;i++) { if(pbest[i]>gbest) { gbest=pbest[i]; gbest_x=pbest_x[i]; } } while(m<=T) { w=(ws-we-d1)*Math.exp(1/(1+d2*(m-1)/T)); for(int i=0;i<num;i++) { g[i]=0.0; ran=Math.random(); v[i]=w*v[i]+c1*ran*(pbest_x[i]-x[i])+c2*ran*(gbest_x-x[i]); x[i]+=v[i]; if(x[i]<x_down) { x[i]=x_down; } if(x[i]>x_up) { x[i]=x_up; } g[i]=func(x[i]); if(pbest[i]<g[i]) { pbest[i]=g[i]; pbest_x[i]=x[i]; } if(pbest[i]>gbest) { gbest=pbest[i]; gbest_x=pbest_x[i]; } } System.out.println("第 "+m+"次迭代的gbest="+gbest+" gbest_x="+gbest_x); m=m+1; } System.out.println("最优位置："+gbest_x+"最优值："+gbest); } }