%?Group?Search?Optimizer?algorithm?in?Matlab
%?Copyright?（C）?2004-2005?Department?of?Electricity?
%?Engineering?and?Electronics?The?University?of?Liverpool?UK
%?Last?modifed?12-May-05

function?[fbestvalbestmemberhistory]?=?GSO（fnameMaxIter）

%?function?[fbestvalbestmember]?=?GSO（fnameNDimMaxIter）
%
%???Run?a?Group?Search?Optimizer?algorithm
%
%?Input?Arguments:
%???fname???????-?the?name?of?the?evaluation?.m?function
%???NDim????????-?dimension?of?the?evalation?function
%???MaxIter?????-?maximum?iteration
close?all
tic%計算程序使用時間，與程序結尾的toc相呼應
global?NDim?Acnumb
flag=0;
iteration?=?0;
%-----------------------------
%MaxIter=10;
Best=zeros（MaxIter1）;
%fname=‘DShapeweightplane10‘;
Bound=eval（fname）;
PopSize=48;?????%?population?of?members
angle=pi/4.*ones（NDim-1PopSize）;???%?Initialize?head?angle
%leftangle=angle;?
%rightangle=angle;

%?Defined?lower?bound?and?upper?bound.
LowerBound?=?zeros（NDimPopSize）;
UpperBound?=?zeros（NDimPopSize）;
for?i=1:PopSize
????LowerBound（:i）=Bound（:1）;
????UpperBound（:i）=Bound（:2）;
end
%DResult?=?1e-1;????%?Desired?results
%-----------------------------%
vmax?=?ones（NDimPopSize）;
for?i=1:NDim
????vmax（i:）=（UpperBound（i:）-LowerBound（i:））;
end
r=norm（vmax（:1））;
distance=r*repmat（ones（1PopSize）NDim1）;
a=?round（（（NDim+1）^.5））;
direction=zeros（NDimPopSize）;
for?j=1:PopSize
????direction（1j）=（cos（angle（1j）））;
????for?i=2:NDim-1
????????direction（ij）=cos（angle（ij））.*prod（sin（angle（i:NDim-1j）））;
????end
????direction（NDimj）=prod（sin（angle（1:NDim-1j）））;
end
initialization=1;
while?initialization==1
????initialization=0;
????%---------包含初始設計值-------%
????population?=??rand（NDim?PopSize）.*（UpperBound-LowerBound）?+?LowerBound;?????%?Initialize?swarm?population
????%population（:1）
????%-----------------------------%
????%-----------add---------------%
????population（1:Acnumb:）=round（population（1:Acnumb:）*100000）/100000;%
????population（Acnumb+1:NDim:）=round（population（Acnumb+1:NDim:）*1000）/1000;%
????exexutefunction=strcat（fname‘（population）‘）;
????%?Evaluate?initial?population
????fvalue?=?eval（exexutefunction）;
????if?min（fvalue）?==?inf
????????initialization=1;
????end
end

%?Finding?best?member?in?initial?population
[fbestvalindex]?=?min（fvalue）;????
bestmember=population（:index）;
oldangle=angle;
oldindex=index;?%#ok
badcounter=0;
while（flag?==?0）?&&?（iteration?????%rand（‘state‘sum（100*clock））
????iteration?=?iteration?+1;

?????for?j=1:PopSize
????????R1=randn（1）;
????????R2=rand（NDim-11）;
????????R3?=?rand（NDim?1）;

????????if?j==index?%?Stop?and?search?around

????????????SamplePosition=zeros（NDim4）;
????????????SampleAngle=zeros（NDim-14）;
????????????SampleValue=zeros（14）;
????????????leftangle=（pi/（a^2））.*R2+angle（:j）;??%trun?left?the?angle?range?between?new?and?old?direction?is?\in?[0pi/a^2]
????????????rightangle=-（pi/（a^2））.*R2+angle（:j）;?%trun?right?the?angle?range?between?new?and?old?direction?is?\in?[0pi/a^2]
????????????
????????????distance（:j）=