%小生境遺傳算法主程序?
function?f=NicheSGA_1;?
%Optimizationstyle=‘Min‘;VarietyNumber=1;Xmin=[0];Xmax=[10*pi];PopulationSize=100;GenotypeLength=20;MaxGeneration=100;L=5;?
Optimizationstyle=‘Max‘;VarietyNumber=1;Xmin=[0];Xmax=[1];PopulationSize=100;GenotypeLength=20;MaxGeneration=100;L=0.18;?
Pc=0.9;Pm=0.01;N=PopulationSize/2;Penalty=10^（-30）;?
initpop=randint（PopulationSizeVarietyNumber*GenotypeLength+VarietyNumber+2）;?
oldpop=decode（OptimizationstyleVarietyNumberPopulationSizeXminXmaxGenotypeLengthinitpop）;?
[Bindex]=sort（oldpop（:VarietyNumber*GenotypeLength+VarietyNumber+2））;?
oldpop=[oldpop（index（:）:）];?
temp1=oldpop（PopulationSize-N+1:PopulationSize:）;?
k=1;?
while?k<=MaxGeneration?
????newpop=selection（VarietyNumberPopulationSizeGenotypeLengtholdpop）;?
????randorder=randperm（PopulationSize）;?
????for?i=1:2:（PopulationSize-1）?
????????parent1=newpop（randorder（i）:）;?
????????parent2=newpop（randorder（i+1）:）;?
????????[child1child2]=crossover（VarietyNumberGenotypeLengthparent1parent2Pc）;?
????????newpop（randorder（i）:）=child1;?
????????newpop（randorder（i+1）:）=child2;?
????end?
????for?i=1:PopulationSize?
????????newpop（i:）=mutation（VarietyNumberGenotypeLengthPmnewpop（i:））;?
????end?
????newpop=decode（OptimizationstyleVarietyNumberPopulationSizeXminXmaxGenotypeLengthnewpop）;?
????faveragevalue（k）=mean（newpop（:VarietyNumber*GenotypeLength+VarietyNumber+1））;?
????if?（Optimizationstyle==‘Max‘）?
????????optimization（kVarietyNumber+1）=max（newpop（:VarietyNumber*GenotypeLength+VarietyNumber+1））;?
????else?
????????optimization（kVarietyNumber+1）=min（newpop（:VarietyNumber*GenotypeLength+VarietyNumber+1））;?
????end?
????[tfindex]=ismember（optimization（kVarietyNumber+1）newpop（:VarietyNumber*GenotypeLength+VarietyNumber+1））;?
????optimization（k1:VarietyNumber）=newpop（indexVarietyNumber*GenotypeLength+1:VarietyNumber*GenotypeLength+VarietyNumber）;?
????temp2=[temp1‘?newpop‘]‘;?
????temp2=niche（PopulationSizeVarietyNumberGenotypeLengthNLPenaltytemp2）;?
????[Bindex]=sort（temp2（:VarietyNumber*GenotypeLength+VarietyNumber+2））;?
????temp2=[temp2（index（:）:）];?
????temp1=temp2（PopulationSize+1:PopulationSize+N:）;?
????oldpop=temp2（N+1:PopulationSize+N:）;?
????k=k+1;?
end?
if?（Optimizationstyle==‘Max‘）?
????optimizationf=max（optimization（:VarietyNumber+1））?
else?
????optimizationf=min（optimization（:VarietyNumber+1））?
end?
[tfindex]=ismember（optimizationfoptimi