clear?all
clc
NP=50
D=3
F=0.8
CR=0.9
eps=1e-9
gen_max=30
ObjFun?=?@DE_FOPID;
count=1
hbound=[20?2?20]
lbound=[0?0?0]
for?i=1:NP
x（i1）=20*rand
x（i2）=2*rand
x（i3）=20*rand
x（i4）=1.5*rand
x（i5）=1.5*rand

end
clear?i
for?i?=?1:NP
????f（i:）?=?feval（ObjFunx（i:））;?
%?????f（i）=mbf_pso（x（i:））;%每一個粒子的適應(yīng)度，即為目標(biāo)函數(shù)值。
end
pbx?=?x;%初值即為個體極值
pbf?=?f;
[gbf?i]?=?min（pbf）;%全局最小值
gbx?=?pbx（i:）;%
while（counteps
????for?i=1:NP
????????while?2>1
????????????a=floor（rand*NP）+1;
????????????if?a~=i
????????????????break;
????????????end
????????end
????????while?2>1
????????????b=floor（rand*NP）+1;
????????????if?b~=i&b~=a
????????????????break;
????????????end
????????end
????????while?2>1
????????????c=floor（rand*NP）+1;
?????????

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

?????文件???????1903??2014-04-25?10:05??DE.m

?????文件????????227??2014-04-25?10:05??DE_FOPID.m

?????文件??????30973??2014-04-01?19:53??FOPID.mdl

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

????????????????33103????????????????????3