test.m

来自「采用标准群粒子优化计算方法」· M 代码 · 共 52 行

clc;
clear all;
global pop;
Fun_Select=8;
popsize=30;
D=2;
x_max=50;
x_min=-50;

x=1:popsize;
%pop(2*popsize,3*D+2)=0;
pop=inf*ones(2*popsize,3*D+2);




for j=1:D
    pop(1,j)=rand;
    u=4;
    for i=1:2*popsize-1
        pop(i+1,j)=u*pop(i,j)*(1-pop(i,j));
    end

   for i=1:2*popsize
       pop(i,j)=pop(i,j)*(x_max-x_min)+x_min;%混沌初始化
   end
    
   %pop_init(:,j)= Fun_bubble( pop_init(:,j),popsize*2);
   %pop(:,j)=pop_init(1:popsize,j);
end
pop

Fun_ComputeFitness(Fun_Select,D,popsize*2)
pop
pop=Fun_bubble(pop,2*popsize,D)

%pop(1)=rand;
%u=4;
%for i=1:popsize-1
%    pop(i+1)=u*pop(i)*(1-pop(i));
%end

%for i=1:popsize
%    pop(i)=pop(i)*4-2;
%end

%plot(x,pop,'r*');
%hold on;

%pop_list=Fun_bubble(pop,popsize);

%plot(x,pop_list,'b.');