| 查看: 1244 | 回复: 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是不是有错 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 融合产物作为模板P不出来，求大神指点！！！
已经有9人回复
- 求大神们帮忙看看这个程序哪儿需要改啊
已经有7人回复
- 用matlab画一个函数，能帮我看看程序哪里有问题吗
已经有8人回复
- 大家帮忙看下这个labview小程序哪里有问题？
已经有5人回复
- 请C语言的人们看看我这个最简单的小程序哪里有问题呢
已经有6人回复
- 求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）
已经有6人回复
- fluent udf请教
已经有16人回复
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
- Matlab中如何自己定义一个新函数，然后计算时可以调用？
已经有12人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- matlab模拟菲涅尔直边衍射
已经有9人回复
- 程序运行不了问题在哪
已经有11人回复
- 这段并行程序有问题，不知道错在什么地方请大家看看！
已经有6人回复
- 请教一个fortran小程序编译出错的问题，谢谢
已经有9人回复
- 三维数组换行输入到文件（intel fortran 编译器）
已经有6人回复
- 哪位帮忙写份Fortran程序？
已经有10人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【讨论】关于谭浩强的书
已经有6人回复
- 【求助】怎么实现JAVA程序编程.exe格式
已经有11人回复
- 【求助】有很多错误的程序，找不出来了
已经有18人回复
- 【原创】win7下安装程序点击没有反应的可能解决方案（bruker topspin）
已经有6人回复
- 【求助】Riccati传递矩阵法求解临界转速（有部分程序）
已经有5人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2014-05-30 22:43:00
lijunxue08
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 700.5
- 帖子: 71
- 在线: 33小时
- 虫号: 2459634
- 注册: 2013-05-11
- 性别: GG
- 专业: 地理信息系统
3楼2014-05-31 09:14:20
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
4楼2014-05-31 11:29:16
lijunxue08
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 700.5
- 帖子: 71
- 在线: 33小时
- 虫号: 2459634
- 注册: 2013-05-11
- 性别: GG
- 专业: 地理信息系统
5楼2014-05-31 18:45:11
charmantang
禁虫 (正式写手)
lijunxue086楼2014-06-01 09:56:54
lijunxue08
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 700.5
- 帖子: 71
- 在线: 33小时
- 虫号: 2459634
- 注册: 2013-05-11
- 性别: GG
- 专业: 地理信息系统
7楼2014-06-01 10:55:02
charmantang
禁虫 (正式写手)
|
本帖内容被屏蔽 |
8楼2014-06-01 17:40:39