?? ga2.m
字號(hào):
%function pop_hga=GA(pop,bound,valu2,pop_eva)
clear;
pc=0.88; pm=0.288; select_size=10;cross_size=20; mut_size=20;
pop_size=50; chrom_size=5; beta=0.01;
bound=[78,102;33,45;27,45;27,45;27,45];
gen=50;
optf=[0,0]; generation=0;
pop(1:pop_size,1)=78+24*rand(pop_size,1);
pop(1:pop_size,2)=33+12*rand(pop_size,1);
pop(1:pop_size,3:5)=27+18*rand(pop_size,3);
for i=1:1:pop_size
valu(i,:)=[p21(pop(i,:)'),i];
end;
for j=1:1:pop_size
for i=1:1:(pop_size-j)
if valu(i,1)<valu(i+1,1)
d=valu(i,:);
valu(i,:)=valu(i+1,:);
valu(i+1,:)=d;
end;
end;
end;
optx=pop(1,:);
funopt=valu(1,1);
funoptx=pop(1,:);
while generation<gen
generation=generation+1;
pops=select3(pop,valu,select_size);
popc=cross1(pop,pc,bound,cross_size);
popm=mutate1(pop,pm,bound,mut_size);
pop(1:select_size,:)=pops;
pop(select_size+1:select_size+cross_size,:)=popc;
pop(cross_size+select_size+1:pop_size,:)=popm;
for i=1:1:pop_size
valu(i,:)=[p21(pop(i,:)'),i];
end;
for j=1:1:pop_size
for i=1:1:(pop_size-j)
if valu(i,1)<valu(i+1,1)
d=valu(i,:);
valu(i,:)=valu(i+1,:);
valu(i+1,:)=d;
end;
end;
end;
if optf(1)<valu(1,1)
optf=valu(1,:);
optx=pop(valu(1,2),:);
end;
subopt(generation)=valu(1,1);
popf(valu(pop_size,2),:)=optx;
if optf(1)>funopt
if (optf(1)-funopt)/abs(funopt)>=0.0005
[funoptx,funopt]=direction(funoptx,optx);
pop(valu(pop_size-1,2),:)=funoptx;
end
end
end;
optx=optx
p21(optx)
x=1:1:gen;
plot(x,subopt(x));
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -