使用遗传算法遇到问题，求助！

xiaoOMming

程序如下：
%m文件
function eval=f(sol)
numv=size(sol,2);
x=sol(1:numv);
eval=-20*exp(-0.2*sqrt(sum(x.^2/numv)))-exp(sum(cos(2*pi*x))/numv)+22.71282;
%m文件
function[sol,eval]=fitness(sol,options)
numv=size(sol,2)-1;
x=sol(1:numv);
7 t% d4 v0 j& }/ Deval=f(x);
+ u; C6 H4 `/ z! s) Aeval=-eval;
+ ^4 p& x9 O% U4 Q# J%主程序
) x9 j8 x0 J; o! _6 q9 @. T& }$ Nclc,clear
bounds=ones(2,1)*[-5 5];
" p" C4 a+ w& `7 ]. X4 ]5 G0 sinitPop=initializega(10,bounds,'fitness');
[p,endPop,bestSols]=ga(bounds,'fitness',[],initPop,[1e-6 0 0],'maxGenTerm',1000,'normGeomSelect',[0.08],['simpleXover'],[0.3],'binaryMutation',[0.1]);
然后是出现的问题：Index exceeds matrix dimensions.
3 F9 ~: m+ f8 ?4 C- G9 K
1 V% q! ], ~! G; X3 j; @) xError in b2f (line 34)
; F3 s/ v7 z5 t2 D  a=bval((cs(i)+1):cs(i+1));


+ s6 Z  P( }7 x( `2 T' `% W) AError in ga (line 194)
' H* a+ b4 \0 \1 e      bPop(bFoundIn,=[gen b2f(startPop(bindx,1:numVar),bounds,bits)...


Error in Untitled4 (line 4)
' Q( r+ y$ {2 F) Z. q[p,endPop,bestSols]=ga(bounds,'fitness',[],initPop,[1e-6 0 0],'maxGenTerm',1000,'normGeomSelect',[0.08],['arithXover'],[2 0;2 3;2 0],'nonUnifMutation',[2 25 3]);
是我的工具箱有问题还是程序有问题，不懂求教大神 遗传工具箱.rar (102.67 KB, 下载次数: 0)

2013-8-23 10:52 上传

点击文件名下载附件
下载积分: 体力 -2 点

工具箱
9 J; t2 {; k4 L* Q8 h0 n5 S8 I
  W/ F* a& o) N6 U7 h. c( |& L9 J