clear;
close?all;
N0=30;??????????%%%%初始種群個數%%
itmax=500;??????%%%迭代次數%%%%
dim=2;??????????%%%%維數%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%???1???定義參數????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

pmax=50;????????%%%%種群最大范圍%%%
smax=5;?????????%%%最大種子數%%%%
smin=0;?????????%%%最小種子數%%%%
n=3;????????????%%%%調和指數，計算方差使用，為固定好的值%%%
delta_initial=3;???%%%方差最大值%%%
delta_final=0.001;????%%%方差最小值%%%?
xini=-100+200*rand（N02）;??%%%隨機產生初始種群，在[01]之間%
???????????????????????????%rand函數返回一個30*2的隨機矩陣，并且在1-200之間
???????????????????????????%本身函數關系為X=xmin+（xmax-xmin）*rand（N0dim）
???????????????????????????%這里定義xmin=-100xmax=100
plot（xini（:1）xini（:2）‘b*‘）???%%%%畫出產生初始雜草種群圖像（藍色）
title（‘產生初始雜草種群‘）;??????%%%%定義標題%%
xlabel（‘z1‘）;?????%%%定義橫坐標%%%
ylabel（‘z2‘）;??????%%%%定義縱坐標%%%%
hold?on;???????????%%%%保持%%%
FB=zeros（itmax1）;???%%%生成itmax=500行*1列的0矩陣%%%

iter=1;???%%%%定義初始迭代次數為1??%%%
while?（iter??????????????????????%%%%循環迭代%%

%%%%%%%%%%%%%%%%%%%%5%%%%%%%%??2???生產繁殖，產生種子??%%%%%%%%%%%%%%%%%%%%%5%%%%%%%%%%%%%%%%%%%%%%%%%%

????Niter=N0;???%%%%定義此時的種群范圍Niter=N0=30?剛剛開始等于初始種群個數%%%
????fitness=xini（:1）.^2+xini（:2）.^2;??%%%計算種群適應度函數%%%%
????Nseed=round（-5/（max（fitness）-min（fitness））*（fitness-min（fitness））+5）;
?????????????????????????????%%%產生種子??round函數輸入任何一個數組，對這個數組
?????????????????????????????%%%里每個元素進行四舍五入取整
?????%%%原公式??Nseed=round（-（smax-smin）/max（fitness）-min（fitness））*（fitness-min（min（fitness）））+smax）
?????
%%%%%%%%%%%%%%%%%%%%%%%%%%%%???3??空間擴散，以正態隨機分布??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????delta_iter=（itmax-iter）^n/（itmax）^n*（delta_initial-delta_final）+delta_final;
??????????????????????????????%%%%%%%%%%%??計算每一代的標準差???%%%%%%%%%%%
?????%%原公式為?delta_iter=（delta_initial-delta_final）*（（itmax-iter）^n）/（（itmax）^n）+delta_final????????????????????
????M=size（Nseed）;???%%%%?size函數返回Nseed矩陣的行數?%%%
????m1=[];
????m2=[];
????a=[];
????b=[];????????????%%%%%??定義m1?m2?a?b四個空矩陣??%%%%
????for?i=1:M（11）????%%%?對于每個個體???M=[501]?即i從1到50循環??%%%%
????????x1_iter=no