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

xiaoOMming

程序如下：
%m文件
function eval=f(sol)
2 M5 L+ M0 c; c" A, V! `numv=size(sol,2);
x=sol(1:numv);
7 T# s  m3 b8 ]7 T9 Neval=-20*exp(-0.2*sqrt(sum(x.^2/numv)))-exp(sum(cos(2*pi*x))/numv)+22.71282;
4 A; v7 a$ _' `* H9 ^%m文件
function[sol,eval]=fitness(sol,options)
numv=size(sol,2)-1;
1 Y: m) x* {- px=sol(1:numv);
eval=f(x);
+ k, W3 r1 e: K& S3 Neval=-eval;
%主程序
clc,clear
bounds=ones(2,1)*[-5 5];
initPop=initializega(10,bounds,'fitness');
* G2 W0 V# o" g5 Q, K  h[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.
; e* y$ i* H8 a8 x  h" z1 @' B
Error in b2f (line 34)
  a=bval((cs(i)+1):cs(i+1));
& L# D3 K5 X$ S
8 ~, h2 t# w# G" V# ]5 f
Error in ga (line 194)
: n( I9 L4 W' c8 j. z' L      bPop(bFoundIn,=[gen b2f(startPop(bindx,1:numVar),bounds,bits)...


. t. [) C+ u. O/ v8 I7 g2 VError in Untitled4 (line 4)
5 Q6 ^' o% a6 e" y1 O[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 点

工具箱
: Z$ C- l' O. `
* C2 O" o4 l: Y  `& b- n