psosa.m

介绍了粒子群模拟退火算法在函数优化中的应用

clc
clear

%初始化粒子群算法参数
n = 20;
m = 6;
[pop,gbest,zbest,fitnessgbest,fitnesszbest,V] = initiation(n,m);
fitnesszbest
%设置参数
wstart = 1.0;
wend = 0.9;
Tmax = 6000;
Vmax = 6;
Vmin = -6;
c1 = 1.49445;
c2 = 1.49445;

decayscale = 0.7;


temperature = 0.00000001;
markovlength = 300;
k = 0;
b = 0;
%循环开始
for jj = 1:20
    for i = 1:markovlength
        k = k + 1;
        w = wstart - (wstart - wend)*(Tmax - k)/Tmax;
        %循环开始
        [pop,V,gbest,zbest,fitnessgbest,fitnesszbest,b] = iterateSAPSO(pop,V,gbest,zbest,fitnessgbest,fitnesszbest,w,c1,c2,temperature,n,m,b);
        y(k) = fitnesszbest;
    end
    temperature = temperature*decayscale;
end
i = 1:6000;
plot(y,i)