MATLAB粒子群算法sizepop,MATLAB数学建模（六） | 粒子群优化（PSO）算法讲解 （上）...

%% Copyright (c) 2015, Yarpiz (www.yarpiz.com)% All rights reserved. Please read the "license.txt" for license terms.%% Project Code: YPEA102% Project Title: Implementation of Particle Swarm Optimization in MATLAB% Publisher: Yarpiz (www.yarpiz.com)%% Developer: S. Mostapha Kalami Heris (Member of Yarpiz Team)%% Contact Info: sm.kalami@gmail.com, info@yarpiz.com%clc;clear;close all;%% Problem DefinitionCostFunction=@(x) Sphere(x);        % Cost FunctionnVar=10;            % Number of Decision VariablesVarSize=[1 nVar];   % Size of Decision Variables MatrixVarMin=-10;         % Lower Bound of VariablesVarMax= 10;         % Upper Bound of Variables%% PSO ParametersMaxIt=1000;      % Maximum Number of IterationsnPop=100;        % Population Size (Swarm Size)% PSO Parametersw=1;            % Inertia Weightwdamp=0.99;     % Inertia Weight Damping Ratioc1=1.5;         % Personal Learning Coefficientc2=2.0;         % Global Learning Coefficient% If you would like to use Constriction Coefficients for PSO,% uncomment the following block and comment the above set of parameters.% % Constriction Coefficients% phi1=2.05;% phi2=2.05;% phi=phi1+phi2;% chi=2/(phi-2+sqrt(phi^2-4*phi));% w=chi;          % Inertia Weight% wdamp=1;        % Inertia Weight Damping Ratio% c1=chi*phi1;    % Personal Learning Coefficient% c2=chi*phi2;    % Global Learning Coefficient% Velocity LimitsVelMax=0.1*(VarMax-VarMin);VelMin=-VelMax;%% Initializationempty_particle.Position=[];empty_particle.Cost=[];empty_particle.Velocity=[];empty_particle.Best.Position=[];empty_particle.Best.Cost=[];particle=repmat(empty_particle,nPop,1);GlobalBest.Cost=inf;for i=1:nPop% Initialize Positionparticle(i).Position=unifrnd(VarMin,VarMax,VarSize);% Initialize Velocityparticle(i).Velocity=zeros(VarSize);% Evaluationparticle(i).Cost=CostFunction(particle(i).Position);% Update Personal Bestparticle(i).Best.Position=particle(i).Position;particle(i).Best.Cost=particle(i).Cost;% Update Global Bestif particle(i).Best.CostVarMax);particle(i).Velocity(IsOutside)=-particle(i).Velocity(IsOutside);% Apply Position Limitsparticle(i).Position = max(particle(i).Position,VarMin);particle(i).Position = min(particle(i).Position,VarMax);% Evaluationparticle(i).Cost = CostFunction(particle(i).Position);% Update Personal Bestif particle(i).Cost