[Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); [Y(gen),I(gen)]=min(ObjV);

trace(1,gen)=min(ObjV);

trace(2,gen)=sum(ObjV)/length(ObjV);

end

unit=bs2rv(Chrom,FieldD);

figure

plot(trace(1,1:gen),'b-o');

hold on

plot(trace(2,1:gen),'r-o');

XY=unit(I(gen),:)

Z=min(Y)

附录三第四问相关程序

g.m:

function y = g( m )

location.altitude = 0;

A=fun1(m(4));

time.year = 2015;

time.month = A(1);

time.day = A(2);

time.hour=8;

time.min =55;

time.sec=0;

time.sec = 00;

time.UTC = +8;

location.longitude=m(1);

titude=m(2);

L=2;

h=[2.323389734

2.329283311

2.293887466

2.282088066

2.261439117

2.252595865

2.237846657

2.23197138

2.21720809

21