clc
clear
popsize=40;
MAXITER=100;
dimension=12;
w_max?=?0.9;????????????????????????
w_min?=?0.4;
v_max?=?6;?????????????????????????????????????????????????????????????????????????????????????????????????
c1?=?2;??????????????????
c2?=?2;
vmin=0.9;
vmax=1.1;
q10min=0;
q10max=0.04;
q24min=0;
q24max=0.04;
tmin=0.9;
tmax=1.1;
xmax=100;
sum1=0;
sum2=0;
mean=0;
st=0;
runno=1;
data1=zeros（runnoMAXITER）;
v=3/2*（rand（popsize12））;?????????????%???初始化當前速度
V=vmin+（vmax-vmin）*rand（popsize61）;
Q10=round（5*rand（popsize11））;
Q24=round（5*rand（popsize11））;
TA=round（10*rand（popsize41））-5;
x=[V?Q10?Q24?TA];
%for?run=1:runno
mean=0;
T=cputime;
pbest=x;
%disp（‘pbest=‘）;
%disp（pbest）;
gbest=zeros（1dimension）;
??for?i=1:popsize;
?????%[PLOSS]?=?runpf（casenamex（i:））
??????f_x（i）=yue（x（i:））;
????f_pbest（i）=f_x（i）;
??end
%disp（‘f_pbest=‘）;
%disp（f_pbest）;??
????g=find（f_pbest==min（f_pbest（1:popsize）））;%找出f_pbest中等于min（f_pbest（1:popsize））元素的下標
????gbest=pbest（g:）;
????f_gbest=f_pbest（g）;
????MINIUM=f_gbest;
??for?t=1:MAXITER;
???????disp（‘MINIUM=‘）;
???????disp（MINIUM）;
???????disp（‘迭代次數(shù)t=‘）;
???????disp（t）;
????w?=?w_max-（w_max-w_min）*t/MAXITER;
???for?i=1:popsize;
???????u（i）=v（i）;
?????v（i:）?=?w.*v（i:）+c1.*rand.*（pbest（i1:dimension）-x（i1:dimension））+c2.*rand.*（gbest（1:dimension）-x（i1:dimension））;
?????if?abs（v（i）-u（i））<0.00001;
?????????v（i）=3/2*（rand（112））;
?????end
?????x（i1:dimension）?=?x（i1:dimension）+v（i1:dimension）;
????????%disp（‘導納矩陣??x（i:）=‘）;
????????%disp（?x（i:））;?
%???????x（i:）=sign（y）.*min（abs（y）xmax）;?
?????????for?j=1:6;
????????????if?x（ij）>1.1;
????????????????x（ij）=1.1;
????????????end;
????????????if?x（ij）<0.95;
????????????????x（ij）=0.95;
????????????end;
?????????end;
?????????for?j=7:8;
????????????if?x（ij）>5;
????????????????x（ij）=5;
????????????end;
????????????x（ij）=round（x（ij））;
?????????end;
?????????for?j=9:dimension;
????????????if?x（ij）>5;
????????????????x（ij）=5;
????????????end;
????????????if?x（ij）<-5;
????????????????x（ij）=-5;
????????????end;
????????????x（ij）=round（x（ij））;
?????????end
?????????%disp（‘矩陣?x=‘）;
?????????%disp（x）;?
?????????for?i=1:40;
?????????????for?j=1:40;
?????????????????if?abs（x（i）-x（j））<0.00001;
?????????????????????f_x（i）=yue（x（i:））;
?????????????????????f_x（j）=yue（x（j:））;
?????????????????????if?f_x（i）?????????????????????????x（i）=3*rand（112）;
?????????????????????else?
?????????????????????????x（j）=3*rand（112）;
?????????????????????end
?????????????????end
?????????????end
?????????end?????????????
?????????f_x（i）=yue（x（i:））;
?????????if?f_x（i）??????????????pbest（i:）=x（i:）;
??????????????f_pbest（i）=f_x（i）;
?????????end
?????????if?f_pbest（i）??????????????gbest=pbest（i:）;
??????????????f_gbest=f_pbest（i）;
?????????end??
????????
????????????MINIUM=f_gbest;????????????????????
??????end
???%?disp（‘矩陣?x=‘）;
????%disp（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3422??2012-02-25?14:25??pso.m

?????文件????????151??2017-12-14?20:23??www.downma.com.txt

-----------?---------??----------?-----??----

?????????????????3573????????????????????2