if time.min>60

time.min=0;

time.hour=time.hour+1;

end

a=sun_position(time,location);

l2(i)=h(j)/tan((a.zenith)/180*pi); time.min=time.min+3;

end

a1=data(:,1);

a2=data(:,2);

l1=sqrt(a1.^2+a2.^2);

l=0;

for k=1:21

l=l+(l1(k)-l2(k))^2;

end

ObjV(j)=l;

end

ObjV=ObjV';

run_genetic.m:

clc;clear,close all;

NVAR=2;

NIND=40;

MAXGEN=100;

LIND=10;

GGAP=0.9;

trace=zeros(2,MAXGEN);

FieldD=[ LIND LIND;

0 0;

10 90;

1 1;

0 0;

1 1;

1 1];

Chrom=crtbp(NIND,LIND*NVAR);

unit=bs2rv(Chrom,FieldD);

ObjV=fun(unit); for gen=1:MAXGEN

FitnV=ranking(-ObjV);

SelCh=select('sus',Chrom,FitnV,GGAP); SelCh=recombin('xovsp',SelCh,0.7); SelCh=mut(SelCh);

unit=bs2rv(SelCh,FieldD);

ObjVSel=fun(unit);

20