A=[A;i*ones(RT(i),1)];
end
N=length(A);
%% 算法参数设置
%迭代次数
KGA=50;
%种群规模
NGA=50;
%变异概率
PmGA=0.2;
LB=zeros(N,1);
UB=ones(N,1);
Alpha=1;
Beta=1;
%% 调用遗传算法
[BESTX,BESTY,ALLX,ALLY]=GAUCP(KGA,NGA,PmGA,LB,UB,A,C,MD,TDG,TDB,CD,ABD,ADD,ST,Alpha,Beta);
%% 整理输出结果
x=BESTX{KGA};
[Y,X,T1,T2,T3,T4,T5,WTS,WTT,C1,C2]=Fit(x,A,C,MD,TDG,TDB,CD,ABD,ADD,ST,Alpha,Beta);
disp('最佳染色体为');
disp(X');
Time0=60*HM(1)+HM(2);
T1=T1+Time0;
T2=T2+Time0;
T3=T3+Time0;
T4=T4+Time0;
T5=T5+Time0;
T1H=floor(T1/60);
T1M=mod(T1,60);
T2H=floor(T2/60);
T2M=mod(T2,60);
T3H=floor(T3/60);
T3M=mod(T3,60);
T4H=floor(T4/60);
T4M=mod(T4,60);
T5H=floor(T5/60);
T5M=mod(T5,60);
disp('配料和装车开始时间(小时/分钟)');
disp([T1H,T1M]');
disp('从基地出发时间(小时/分钟)');
disp([T2H,T2M]');
disp('到达施工地时刻(小时/分钟)');
disp([T3H,T3M]');
disp('从施工地返回时刻(小时/分钟)');
disp([T4H,T4M]');
disp('到达基地时刻(小时/分钟)');
disp([T5H,T5M]');