function?DE（GmF）%標準差分進化算法程序
F=0.8;%F是變異因子?
Gm=200;?%最大迭代次數
Np=50;?%種群規模
CR=0.5;??%交叉概率常數
G=1;%初始化進化代數
D=10;%所求問題的維數
ge=zeros（1Np）;%各代的最優值
bestx=zeros（NpD）;%各代的最優解
xmin=0;xmax=pi;
function?y=f（Z）;
????y=sum（sin（Z））;
end
X0=（xmax-xmin）*rand（NpD）+xmin;%產生初始種群
X=X0;
%%%%%%%%%%變異操作
V=zeros（NpD）;%變異矢量初始化
U=zeros（NpD）;%初始化
X1=zeros（NpD）;%初始化
value=zeros（1Np）;
while?G<=Gm
????for?i=1:Np
????%產生jkp三個不同的數
????a=1;b=Np;
????dx=randperm（b-a+1）;
????j=dx（1）;k=dx（2）;p=dx（3）;
????if?j==i
????????j=dx（4）;
????elseif?k==i
????????k=dx（4）;
????elseif?p==i
????????p=dx（4）;
????end
????%namd=exp（1-Gm/（Gm+1-G））;%變異算子
????%F=F0*2.^namd;
????%bon=X（p:）+F*（X（j:）-X（k:））;
????bon=X（p:）+F*（X（j:）-X（k:））;
????if?（bon>xmin）&（bon???????V（i:）=bon;
????else?V（i:）=（xmax-xmin）*rand（1D）+xmin;
????end
????end
%%%%%%%%%%%%交叉操作
????for?i=1:Np
????????for?j=1:D
????????????if?（rand<=CR|j==rand（D））?%利用二項分布來交叉
????????????????U（ij）=V（ij）;
????????????else
????????????????U（ij）=X（ij）;
????????????end
????????end
????end
?%%%%%%%%%%%%%%%選擇操作?
?for?i=1:Np
?????if?f（U（i:））>f（X（i:））
?????????X1（i:）=U（i:）;
?????else?
?????????X1（i:）=X（i:）;
?????end
?end
?%找出最大值
?????for?i=1:Np
??????????value（i）=f（X1（i:））;
?????end
????[fmaxnmax]=max（value）;
????ge（G）=fmax;
????bestx（G:）=X1（nmax:）;?
????G=G+1;
????X=X1;
end
ii=linspace（1GmGm）;
plot（iige）
[gmaxn]=max（ge）;
bestvalue=gmax
bestsolution=bestx（n:）
end
????
????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-09-11?10:16??標準差分進化算法源代碼\
?????文件??????????55??2018-08-27?11:33??標準差分進化算法源代碼\【源碼使用必讀】.url
?????文件????????1636??2018-09-11?10:11??標準差分進化算法源代碼\a.m