最近在编写一个粒子群优化算法时遇到这样的问题:我用变量best_fitness记录每次迭代过程最优解,再将其存入一个数组best_in_history中。因为对惯性因子w取值的不同绘制一个适应值的曲线,我想通过循环来实现不同w时对应的最优解的数组。我的循环如下:
for j=1:2
w=wpso(1,j); %存放w的一个数组
best_fitness=inf;
if j==1
for exetime=1:gen %exetime为迭代次数
Radapting; %适值计算函数
best_in_history1(exetime)=best_fitness;
outputdata1; %输出显示
updatepop; %速度和位置改变的函数
pause(0.0001);
end
pop(popsize,8)=pop1(popsize,8);
elseif j==2
for exetime=1:gen
Radapting
best_in_history2(exetime)=best_fitness;
outputdata2;
updatepop;
pause(0.0001);
end
end
%pop(popsize,8)=pop1(popsize,8);
%best_fitness=inf;
end
然而我通过Array Editor发现best_in_history中记录不是不同w对应的值,而best_in_history2是在best_in_history1的基础上继续优化的解,而我想让best_in_history2记录w=wpso(2,2)时对应的一个新的最优路线。请问我该如何对这个循环改进?