24小时热门版块排行榜    

查看: 1235  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lijunxue08

金虫 (小有名气)

[求助] 下面的程序哪里有问题已有2人参与

figure(1);
fplot('v.*sin(10*pi*v)+2.0',[-1,2]);
NIND=40;
MAXGEN=30;
PRECI=20;
GGAP=0.9;
trace=zeros(2,MAXGEN);
FieldD=[20;-1;2;1;0;1;1];
Chrom=crtbp(NIND,PRECI);
gen=0;
v=bs2rv(Chrom,FieldD);
ObjV=v.*sin(10*pi*v)+2.0;
while gen<MAXGEN,
    FitnV=ranking(-ObjV);
    SelCh=select('sus',Chrom,FitnV,GGAP);
    FitnVmax=max(FitnV);
    FitnVave=sum(FitnV)/NIND;
      if(FitnV>= FitnVave)
         Index1=find(FitnV>= FitnVave);
         Index2=find(FitnV< FitnVave);
         Pc1=0.5*(FitnVmax- FitnV(Index,1))/(FitnVmax- FitnVave);
         Pm1=0.02*(FitnVmax- FitnV(Index,1))/(FitnVmax- FitnVave);
     else   
         Pc2=0.85;
         Pm2=0.05;
    end
  end
    SelCh=recombin('xovsp',SelCh(Index1,1),Pc1);
    SelCh=recombin('xovsp',SelCh(Index2,1),Pc2);
    SelCh=mut(SelCh(Index1,1),Pm1);
    SelCh=mut(SelCh(Index2,1),Pm2);
    v=bs2rv(SelCh,FieldD);
    ObjVSel=v.*sin(10*pi*v)+2.0;
    [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel);
    gen=gen+1;
    variable=bs2rv(Chrom, FieldD)
    [Y,I]=max(ObjV),hold on;
    plot(I,Y,'bo');
    trace(1,gen)=max(ObjV);
    trace(2,gen)=sum(ObjV)/length(ObjV);
    if (gen==20)
        figure(2);
        plot(ObjV);hold on;
        plot(ObjV,'b*');grid;
    end
%end
figure(3);
plot(trace(1,','Pr');
hold on;
plot(trace(2,','-.');grid;
legend('解的变化','种群均值的变化')


程序里的Index是不是有错
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

感谢参与,应助指数 +1
不贴报错信息,谁也不知道哪儿有问题啊
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
2楼2014-05-30 22:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lijunxue08 的主题更新
信息提示
请填处理意见