outputdata_2dim.m

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

%实时输出结果
global flag1;
%输出当前种群中粒子位置
subplot(1,2,1);  %子图1
for i=1:popsize %画粒子群中所有粒子（用蓝色星号标记）
    plot(pop(i,1),pop(i,2),'b*');
    hold on;
end

plot(gbest_x1,gbest_x2,'r.','markersize',15);axis([x1_min,x1_max,x2_min, x2_max]);%画全局最优粒子位置（用红点标记）
xlabel('x1(第一维)');
ylabel('x2(第二维)');
title('所有粒子分布示意图与最优粒子位置')
hold off;
grid on; %产生网格


%输出当前最优解适应值
subplot(1,2,2); %子图2
if exetime-1>0
    plot(exetime,best_fitness,'m.');
    line([exetime-1,exetime],[best_in_history(exetime-1),best_fitness]);
    xlabel('迭代次数');
    ylabel('当前最优解适应值');
    title('当前最优解适应值');
    hold on;   
end
grid on; %产生网格