%煙花算法進(jìn)行函數(shù)優(yōu)化
%fitness適應(yīng)度函數(shù)，N煙花數(shù)，D變量維數(shù)，M變異火花數(shù)，Er爆炸半徑，En爆炸數(shù)目
%LBUB分別為變量上下界，T為迭代次數(shù)，ab為爆炸數(shù)目限制因子
clear;clc
N=100;???%?N煙花數(shù)
D=2;?????%?D變量維數(shù)
M=5;?????%?M變異火花數(shù)
En=6;????%?En爆炸數(shù)目
Er=5;????%?Er爆炸半徑
a=0.3;???%?ab為爆炸數(shù)目限制因子
b=0.6;
T=500;???%?T為迭代次數(shù)

%求最大值變量上下界
LB=[-5.12-5.12];
UB=[5.125.12];

%隨機(jī)在解空間初始化N個(gè)煙花位置
x?=?zeros（ND）;
for?i=1:N
????x（i:）=LB+rand（1D）.*（UB-LB）;
end
%循環(huán)迭代
E_Spark=zeros（TDN）;
Fit?=?zeros（1N）;
F?=?zeros（1T）;
for?t=1:T
????%計(jì)算每個(gè)煙花適應(yīng)度值
????for?i=1:N
????????Fit（i）=fitness（x（i:））;
????end
????[F（t）~]=min（Fit）;
????Fmin=min（Fit）;
????Fmax=max（Fit）;
????%計(jì)算每個(gè)煙花的爆炸半徑E_R和爆炸數(shù)目E_N以及產(chǎn)生的爆炸火花
????E_R?=?zeros（1N）;
????E_N?=?zeros（1N）;
????for?i=1:N
????????E_R（i）=Er*（（Fit（i）-Fmin+eps）/（sum（Fit）-N*Fmin+eps））;??%爆炸半徑
????????E_N（i）=En*（（Fmax-Fit（i）+eps）/（N*Fmax-sum（Fit）+eps））;??%爆炸數(shù)目
????????if?E_N（i）????????????E_N（i）=round（a*En）;
????????elseif?E_N（i）>b*En
????????????E_N（i）=round（b*En）;
????????else
????????????E_N（i）=round（E_N（i））;
????????end
????????%產(chǎn)生爆炸火花?E_Spark
????????for?j=2:（E_N（i）+1）??????????????%?第i個(gè)煙花共產(chǎn)生E_N（i）個(gè)火花
????????????E_Spark（1:i）=x（i:）;??????%?將第i個(gè)煙花保存為第i個(gè)火花序列中的第一個(gè)，爆炸產(chǎn)生的火花從序列中的第二個(gè)開始存儲(chǔ)（即煙花為三維數(shù)組每一頁的第一行）
????????????h=E_R（i）*（-1+2*rand（1D））;??%?位置偏移
????????????E_Spark（j:i）=x（i:）+h;????%?第i個(gè)煙花（三維數(shù)組的i頁）產(chǎn)生的第j（三維數(shù)組的j行）個(gè)火花
????????????for?k=1:D???%越界檢測(cè)
????????????????if?E_Spark（jki）>UB（k）||E_Spark（jki）????????????????????E_Spark（jki）=LB（k）+rand*（UB（k）-LB（k））;???%映射規(guī)則
????????????????end
????????????end
????????end
????end
????%產(chǎn)生高斯變異火花Mut_Spark隨機(jī)選擇M個(gè)煙花進(jìn)行變異
????