粒子群算法
fucntin moan()
clc;ileac rlalclos;ea ll;t
i;c %程 序运计时行
0=0.0E01 ; 允%许差误M
axNmu=00; 1 粒%子大最迭次数代
narvs=1 ; 目标函%的自变数量数个p
raitclesiz=30e ; % 子群粒模规c1=2
; %每个 粒的子体个习因子学也称，加速为常数
c2=2; 每个粒%的子社会习学因子，称也为速常加
w数=.60 ; %惯性因 子
vamx0.=8; % 粒的子最飞翔速大度
=-5x1+0r*and(paritcesilezn,ras); v % 粒子在的所置
位v=2ra*dnpar(ictleszei,anvr);s %粒 的飞子速翔
%度i用niln定义e适度函应数以便子函数文将件与主序文程件在放一，起%目标
数函：y=是1(2.+1*1(x+-2*x.2^.)e*x(-x.p^22/))%
iniln命e定令适义应度数函下如：fit
enss=ilinn('1/e(+(121.*(1-+x*x2^2..)e*xp-x.(2^/2)))''x');
%,inlne定i义适应的函度数会使序程运行度速大降大低for
i =:1arpitclsizee
fo j=r1:nrvasf(
)i=itnfse(s(i,x));
enj
enddp
resnoabestlx=_;x
prsenaoblset_afvl=af;[
glbolbest_aafval i=]mniperso(anlbse_favtal;)
lgobablsetx_pe=rsnoabels_xti(:),;k=1
;
hiwel
orfi =:1arptclisizee
or fj1=na:rv
s(if)=fintsse(x(i,j));
ndei
f(f)
persnolbesta_avalf(i=f)i()
;prsenaolbste_x(,:i)=(ix,:;)
ed
nnde[
logalbesbt_avflai ]=mn(iprsenaolbse_tfval)a
;logalbebst_=xepsronalesbt_(i,:);xf
roi 1:pa=ritlcsezie %更粒新群里子每个体的最新个置位
v(i,)=:w*v(,i)+:c*r1an*(depsronlabstex_i(:,)x-i(,):)...
c2*r+ad*(gnolalbbes_txx(i-:,);)f
r j=o:n1avrs % 断判子粒的飞翔度是否速过超最了大翔飞速度
if (i,v)>jmvxa
;v(ij),=maxv;else
f vi(i,)
i(,j)=v-amx;
end
nd
e(xi:)=x,i(:,)+(v,i:;)end
ifabs( golblaesbtfaval_
=k+k1;en
d
Valu1=1/geoballbet_sfaval1-;Va ule=num1s2rtV(luae)1
% ;srcatt指可令实现字符以组合的出输
dsp(isrcat('the tmaxiumm vlae',u=',Value1));'
%出输最值大所在横的标坐位置
Vlaue=glo2babelstx_ V;auel=nu22mtrsVa(ul2e)
di;sps(rtcat('hteco resporndngi cordoinate''=,,V'luea)2)
;=x-:050.1:5
y=;2.*(1-x12+*x.2^.)ex*p-(x.2/2);^
plt(x,yo'm-',,l'ienwitdh,'3);h
ol ond;
lopt(gloalbbse_xt,/1logalbebs_tafal-v1,kp'','ilnwedtih'4,);
legend'(目标函'数'搜,到索的大值'最;)lxbal(e'')x;yabell('y');grid on;tco
;