%?本函數用LFOA求取Sphere函數的極小值
%%
clc
clear?all
tic
c=20;
for?z=1:c
a=0.5;weishu=30;
%%?隨機初始果蠅種群位置
x_axis=200*rand（1weishu）-100;
y_axis=200*rand（1weishu）-100;
%?設置參數
maxdiedai=150;?%迭代次數
sizepop=30;????%種群規模

%?果蠅開始尋優，利用嗅覺尋找食物
for?i=1:sizepop
????%?賦予果蠅個體利用嗅覺搜尋食物的隨機方向和距離
????x（i:）=x_axis+2*rand（）-1;
????y（i:）=y_axis+2*rand（）-1;
????%計算味道濃度S，為距離的倒數?并判斷是否超過邊界
????for?m=1:weishu
???????d（im）=（x（im）^2+y（im）^2）^0.5;
???????s（im）=1/d（im）;
???????v（im）=s（im）;?
????end
????%果蠅個體位置的味道濃度
????smell（i）=sum（v（i:）.^2）;
????%?適應度函數
end
%找出果蠅群體中味道濃度最大的果蠅
[bestsmellbestindex]=min（smell）;
[worstsmellworstindex]=max（smell）;
%?保留最佳味道濃度值與坐標?
smellbest=bestsmell;

xb=x（bestindex:）;
yb=y（bestindex:）;%得到最大適應度值的信息

xw=x（worstindex:）;
yw=y（worstindex:）;%得到最小適應度值的信

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

?????文件???????3419??2015-11-23?17:55??SphereLFOA.m

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

?????????????????3419????????????????????1