基本粒子群算法，MATLAB文件的M文件编写

function [pso F] = pso_2D（）
uav_size=4;
pop_size = 30;          %   pop_size 种群大小
part_size = 2*uav_size;          %   part_size 粒子大小
% Value=[0.9 0.8 0.7 0.6];
% p_kill=[0.1 0.1 0.2 0.9;0.8 0.8 0.2 0.1;0.3 0.4 0.5 0.2;0.4 0.9 0.3 0.3]
% Distance=[0.91813 0.56140 0.93860 0.20175;0.16374 0.93860 0.35965 0.35965;0.31579 0.97368 0.20760 0.67544;0.36008 0.75642 0.82016 0.5573]
% p_survial=[0.926 0.906 0.916 0.966;0.306 0.906 0.986 0.926;0.996 0.916 0.256 0.926;0.926 0.936 0.926 0.916]
gbest = zeros（1part_size+1）;            %   gbest 当前搜索到的最小的值
max_gen = 6000;          %   max_gen 最大迭代次数
arr_present = ini_pos（pop_sizepart_size）;   %   present 当前位置随机初始化rand（）的范围为0~1

v=ini_v（pop_sizepart_size）;             %   初始化当前速度


pbest = zeros（pop_sizepart_size+1）;      %   pbest 粒子以前搜索到的最优值，最后一列包括这些值的适应度
w_max = 0.9;                            %   w_max 权系数最大值
w_min = 0.4;
v_max = 0.2;             %                                                                                          **最大速度为粒子的范围宽度
c1 = 2;                   %   学习因子
c2 = 2;                   %   学习因子
best_record = zeros（1max_gen）;     %   best_record记录最好的粒子的适应度。
%  ————————————————————————
%   计算原始种群的适应度及初始化
%  ————————————————————————
arr_present（:end）=ini_fit（arr_presentpop_sizepart_size）;

% for k=1:pop_size
%     present（kend） = fitness（present（k1:part_size））;  %计算原始种群的适应度
% end

pbest = arr_present;                                        %初始化各个粒子最优值
[best_value best_index] = max（arr_present（:end））;         %初始化全局最优，即适应度为全局最小的值，根据需要也可以选取为最大值
gbest = arr_present（best_index:）;

for i=1:max_gen
    w = w_max-（w_max-w_min）*i/max_gen;
    for j=1:pop_size
        v（j:） = w.*v（j:）+c1.*rand.*（pbest（j1:part_size）-arr_present（j1:part_size））...
            +c2.*rand.*（gbest（1:part_size）-arr_present（j1:part_size））;                        %  粒子速度更新 （a）
        v（j:）=round（v（j:））;
%         c = find（abs（v）>0.2）;                                                                                              %**最大速度设置，粒子的范围宽度
%         v（c） = sign（v（c））*0.2;   %如果速度大于3.14则，速度为3.14

        arr_present（j1:part_size）= arr_present（j1:part_size）+v（j1:part_size）;              %  粒子位置更新 （b）
%         arr_present（j1:part_size）=round（arr_present（j1:part_size））
               c = find（arr_present>4）; 
             arr_present（c） = 4; 
             d =