粒子群算法优化问题matlab程序

function [pso F] = pso_2D（）
%   FUNCTION PSO  --------USE Particle Swarm Optimization Algorithm
%global present;
% close all;
pop_size = 10;          %   pop_size 种群大小
part_size =2;          %   part_size 粒子大小                                                                      ** =n-D
gbest = zeros（1part_size+1）;            %   gbest 当前搜索到的最小的值
max_gen = 80;          %   max_gen 最大迭代次数
region=zeros（part_size1）;  %   设定搜索空间范围
region=[-33;-33];          %                                                                                      **每一维设定不同范围


rand（‘state‘sum（100*clock））;   %   重置随机数发生器状态
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] = min（arr_present（:end））;         %初始化全局最优，即适应度为全局最小的值，根据需要也可以选取为最大值
gbest = arr_present（best_index:）;
%v = zeros（pop_size1）;          %   v 速度
%  ————————————————————————
%   迭代
%  ————————————————————————
% global m;
% m = moviein（1000）;      %生成帧矩阵
x=[-3:0.01:3];
y=[-3:0.01:3];
z=@（xy） 3*（1-x）.^2.*exp（-（x.^2） - （y+1）.^2） ...
    - 10*（x/5 - x.^3 - y.^5）.*exp（-x.^2-y.^2） ...
    - 1/3*exp（-（x+1）.^2 - y.^2）;
for i=1:max_gen
    grid on;
    %     plot3（xyz）;
    %     subplot（121）ezmesh（z）hold ongrid onplot3（arr_present（:1）arr_present（:2）arr_present（:3）‘*‘）hold off;
    %     subplot（122）ezmesh（z）view（[14590]）hold ongrid onplot3（arr_present（:1）arr_present（:2）arr_present（:3）‘*‘）hold off;
    ezmesh（z） hold ongrid onplot3（arr_present（:1）arr_present（:2）arr_present（:3）‘*‘）hold off;

    drawnow
    F（i）=getframe;

    %     ezmesh（z）
    % %     view（[-3790]）
    %     hold on;
    %     grid on;
    %     %   plot（-0.08980.7126‘ro‘）;
    %     plot3（arr_present（:1）arr_present（:2）arr_present（:3）‘*‘）;                                                  %改为三维
    %    axis（[-2*pi2*pi-pipi-5010]）;
    %     hold off;
    pause（0.01）;
    %     m（:i） = getframe;        %添加图形
    w = w_max-（w_max-w_min）*i/max_gen;
    %    fprintf（‘#  %i 代开始！\n‘i）;
    %   确定是否对打散已经收敛的粒子群——————————————————————————————
    reset = 0;          %   reset = 1时设置为粒子群过分收敛时将其打散，如果＝1则不打散
    if