initial1.m.txt

求解对问题的matlab程序直接运行就可的到结果

%初始化
gen=50;                            %设置进化代数
popsize=15;                        %设置种群规模
max_velocity=1;                    %最大速度限制
best_fitness=inf;                  
pop(popsize,14)=0;                 %初始化种群，创建popsize行14列的0矩阵，
                                   %种群数组第1—4列为x1,x2,x3,x4的值，
                                   %第5-8列为x1,x2,x3,x4速度分量，
                                   %第9-12列为个体最优位置的x1,x2,x3,x4的值，
                                   %第13列为个体最优适值，第14列为当前个体适应值
for i=1:popsize
    for j=1:4
        pop(i,j)=rand();           %初始化种群粒子的位置
        pop(i,j+4)=rand()*2-1;     %初始化种群粒子的速度
    end
    for m=9:12    
        pop(i,m)=pop(i,m-8);       %初始状态下个体最优值等于初始位置
    end
    pop(i,13)=inf;
    pop(i,14)=inf;
end
c1=2.05;
c2=2.05;
x1_min=0;
x2_min=0;
x3_min=0;
x4_min=0;