%標準粒群優(yōu)化算法程序
% 2007.1.9 By 郟宣耀
%測試函數(shù)：f(x,y)=100(x^2-y)^2+(1-x)^2, -2.048<x,y<2.048
%求解函數(shù)最小值
clear;
global popsize;     %種群規(guī)模
global n;           %粒子維數(shù)
%global popnum;      %種群數(shù)量
global pop;         %種群
%global c0;          %速度慣性系數(shù),為0—1的隨機數(shù)
global c1;          %個體最優(yōu)導(dǎo)向系數(shù)
global c2;          %全局最優(yōu)導(dǎo)向系數(shù)
global gbest;       %全局最優(yōu)解坐標
global best_fitness;    %最優(yōu)解
global best_in_history; %最優(yōu)解變化軌跡
global x_min;           %x的下限
global x_max;           %x的上限
global y_min;           %y的下限
global y_max;           %y的上限
global gen;             %迭代次數(shù)
global exetime;         %當(dāng)前迭代次數(shù)
global max_velocity;    %最大速度


initial;        %初始化

for exetime=1:gen
    %outputdata;     %實時輸出結(jié)果
    adapting;       %計算適應(yīng)值
    %errorcompute(); %計算當(dāng)前種群適值標準差
    updatepop;      %更新粒子位置
    pause(0.01);
end


%clear i;
%clear exetime;
%clear x_max;
%clear x_min;
%clear y_min;
%clear y_max;