上呈：基本的遗传算法程序

54hjhf 发表于 2004-9-9 16:40:03

下面的程序是用遗传算法求函数：f(x1,x2)=100(x1*x1-x2*x2) (x1*x1-x2*x2)+(1-x1)(1-x1) 
约束：-2.048<=x1,x2<=2.048的极大值
 
%Generic Algorithm for function f(x1,x2) optimum
clear all;
close all;
 
%parameters
Size=80;
G=100;
Codel=10;
umax=2.048;
umin=-2.048;
E=round(rand(Size,2*Codel));%initial code
 
%Main program
for k=1:1:G
 time(k)=k;
 for s=1:1:Size
 m=E(s,:);
 y1=0;
 y2=0;
 
 %uncoding
 m1=m(1:1:Codel);
 for i=1:1:Codel
 y1=y1+m1(i)*2^(i-1);
 end
 x1=(umax-umin)*y1/1023+umin;
 m2=m(Codel+1:1:2*Codel);
 for i=1:1:Codel
 y2=y2+m2(i)*2^(i-1);
 end
 x2=(umax-umin)*y2/1023+umin;
 
 F(s)=100*(x1^2-x2)^2+(1-x1)^2;
 end
 
 Ji=1./F;
 %***********Step 1:Evaluate BestJ**********
 BestJ(k)=min(Ji);
 fi=F; %Fitness Function
 =sort(fi); %Arranging fi small to bigger
 Bestfi=Oderfi(Size); %LetBestJfi=max(fi)
 BestS=E(Indexfi(Size),:); %LetBestS=E(m),m is the Indexfi belong to max(fi)
 bfi(k)=Bestfi;
 
 %**********Step 2:Select and Reproduct Operetion******
 fi_sum=sum(fi);
 fi_Size=(Oderfi/fi_sum)*Size;
 
 fi_S=floor(fi_Size); %Selecting Bigger fi value
 
 kk=1;
 for i=1:1:Size%Selece and Reproduce
 for j=1:1:fi_S(i)
 <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(kk,:)=E(Indexfi(i),:);
 kk=kk+1; %kk is used to roproduce
 end;
 end
 
 %**********Step 3:Crossover Operation***********
 pc=0.60;
 n=ceil(20*rand);
 for i=1:2:(Size-1)
 temp=rand;
 if pc>temp %crossover condition
 for j=n:1:20
 <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(i,j)=E(i+1,j);
 <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(i+1,j)=E(i,j);
 end
 end
 end
 <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(Size,:)=BestS;
 E=<st1:City w:st="on"><st1:place w:st="on">TempE</st1:place></st1:City>;
 
 %*********Step 4:Mutation Operation**********
 %pm=0.001;
 %pm=0.001-*(0.001)/Size;%Bigger fi,smaller Pm
 %pm=0.0; %No mutation
 pm=0.1; %Big mutation
 
 for i=1:1:Size
 for j=1:1:2*Codel
 temp=rand;
 if pm>temp %Mutation condition
 if <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(i,j)==0
 <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(i,j)=1;
 else
 <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(i,j)=0;
 end
 end
 end
 end
 
 %Guarantee TempPop(30,:)is the code belong to the best individual (maxz(fi))
 <st1:place w:st="on"><st1:City w:st="on">TempE</st1:City></st1:place>(Size,:)=BestS;
 E=<st1:City w:st="on"><st1:place w:st="on">TempE</st1:place></st1:City>;
end
 
Max_Value=Bestfi
BestS
x1
x2
figure(1);
plot(time,BestJ);
xlabel('Times');ylabel('bestJ');
figure(2);
plot(time,bfi);
xlabel('times');ylabel('Best F');

afnpf 发表于 2004-9-15 06:44:39

<img src="http://www.shumo.com/bbs/Skins/Default/emot/em01.gif">够辛苦的，试试再说！

rraauull 发表于 2004-9-18 20:03:21

植物基因AA，Aa,aa 比例相同。现与Aa杂交结果如下表，纵坐标表示为结果比例，建模说明5代后情况：<TABLE medium none; BORDER-TOP: medium none; BORDER-LEFT: medium none; BORDER-BOTTOM: medium none; BORDER-COLLAPSE: collapse; mso-border-alt: solid windowtext .5pt; mso-padding-alt: 0cm 5.4pt 0cm 5.4pt" cellSpacing=0 cellPadding=0 border=1><TR><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: windowtext 0.5pt solid; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: windowtext 0.5pt solid; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent" vAlign=top>杂交</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: windowtext 0.5pt solid; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt" vAlign=top>AA</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: windowtext 0.5pt solid; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt" vAlign=top>Aa</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: windowtext 0.5pt solid; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt" vAlign=top>aa</TD></TR><TR><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: windowtext 0.5pt solid; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-top-alt: solid windowtext .5pt" vAlign=top>AA</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0.5</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0.25</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0</TD></TR><TR><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: windowtext 0.5pt solid; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-top-alt: solid windowtext .5pt" vAlign=top>Aa</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0.5</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0.5</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0.5</TD></TR><TR><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: windowtext 0.5pt solid; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-top-alt: solid windowtext .5pt" vAlign=top>aa</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0.25</TD><TD windowtext 0.5pt solid; PADDING-RIGHT: 5.4pt; BORDER-TOP: #e0dfe3; PADDING-LEFT: 5.4pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: #e0dfe3; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 0.5pt solid; BACKGROUND-COLOR: transparent; mso-border-left-alt: solid windowtext .5pt; mso-border-top-alt: solid windowtext .5pt" vAlign=top>0.5</TD></TR></TABLE>帮我算下这个

yuexu33 发表于 2006-8-21 00:49:24

8错

lily1125 发表于 2007-8-3 20:17:56

回复 #1 54hjhf 的帖子

尊敬的搂主：您问你的reproduce算子是用的什么方法，一般的不都是轮盘算法吗？能解释一下你这样写的原理吗？我研究楼主的程序两天了．

wzwayj 发表于 2009-10-23 14:33:41

怎么出现这么多微笑的表情啊·

hlj031140 发表于 2011-4-16 18:46:25

怎么都没有注释啊？

hlj031140 发表于 2011-4-16 18:47:46

怎么都没有注释啊？

qingfengli 发表于 2011-7-16 10:16:45

恩呢谢谢复制了

泪天堂 发表于 2011-8-25 13:00:39

可以用吗？

页: [1]

数模论坛's Archiver

上呈：基本的遗传算法程序

回复 #1 54hjhf 的帖子