自适应遗传算法matlab代码

frogfish

matlab程序是雷英杰课本中的第一个算例,
将其改为自适应的代码.

1. figure(1);
   
2. fplot('v.*sin(10*pi*v)+2.0',[-1,2]);
   
3. NIND=40;
   
4. MAXGEN=30;
   
5. PRECI=20;
   
6. GGAP=0.9;
   
7. trace=zeros(2,MAXGEN);
   
8. FieldD=[20;-1;2;1;0;1;1];
   
9. Chrom=crtbp(NIND,PRECI);
   
10. gen=0;
    
11. v=bs2rv(Chrom,FieldD);
    
12. ObjV=v.*sin(10*pi*v)+2.0;
    
13. while gen<MAXGEN,
    
14.     FitnV=ranking(-ObjV);
    
15.     SelCh=select('sus',Chrom,FitnV,GGAP);
    
16.     FitnVmax=max(FitnV);
    
17.     FitnVave=sum(FitnV)/NIND;
    
18.       if(FitnV>= FitnVave)
    
19.          Index1=find(FitnV>= FitnVave);
    
20.          Index2=find(FitnV< FitnVave);
    
21.          Pc1=0.5*(FitnVmax- FitnV(Index,1))/(FitnVmax- FitnVave);
    
22.          Pm1=0.02*(FitnVmax- FitnV(Index,1))/(FitnVmax- FitnVave);
    
23.      else   
    
24.          Pc2=0.85;
    
25.          Pm2=0.05;
    
26.     end
    
27.   end
    
28.     SelCh=recombin('xovsp',SelCh(Index1,1),Pc1);
    
29.     SelCh=recombin('xovsp',SelCh(Index2,1),Pc2);
    
30.     SelCh=mut(SelCh(Index1,1),Pm1);
    
31.     SelCh=mut(SelCh(Index2,1),Pm2);
    
32.     v=bs2rv(SelCh,FieldD);
    
33.     ObjVSel=v.*sin(10*pi*v)+2.0;
    
34.     [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
    
35.     gen=gen+1;
    
36.     variable=bs2rv(Chrom, FieldD)
    
37.     [Y,I]=max(ObjV),hold on;
    
38.     plot(I,Y,'bo');
    
39.     trace(1,gen)=max(ObjV);
    
40.     trace(2,gen)=sum(ObjV)/length(ObjV);
    
41.     if (gen==20)
    
42.         figure(2);
    
43.         plot(ObjV);hold on;
    
44.         plot(ObjV,'b*');grid;
    
45.     end
    
46. %end
    
47. figure(3);
    
48. plot(trace(1,:)','Pr');
    
49. hold on;
    
50. plot(trace(2,:)','-.');grid;
    
51. legend('解的变化','种群均值的变化')

复制代码

wy558558558

这是在哪节选的，程序中index在哪？就看到index1和index２了．还有上来就画图:@o

frogfish

原帖由 wy558558558 于 2007-6-29 20:12 发表

                               
登录/注册后可看大图

这是在哪节选的，程序中index在哪？就看到index1和index２了．还有上来就画图:@o

画图有问题吗

gh22

程序调试不通啊，问题在哪呢？高人请指点一下，先谢谢了！

aduktxj

最好是能有注释之类的，这样看有点不是很清晰，谢谢楼主！！:'( :'(

heqiangyong

运行时出现问题！不知道楼主发现了没？

nkdtxf

运行有问题，不过谢谢楼主

gauss101

Error: The input character is not valid in MATLAB statements or expressions.

hjd2003

学习学习！！！！！！呵呵

hjd2003

谢谢楼主分享！学习学习！！！！！

lujiujiji

学习了……

dyhkxydfbb

下来研究下，谢谢


用2010运行好像一直都是busy  出不来结果？请问您运行成功吗？

pingxuyufeng

是的，有问题，我的也挑不出来

hfuthj

先来看看