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

xiaoOMming

程序如下：
: B! c  M$ x1 [( j4 X" @! W" `/ U/ ]%m文件
function eval=f(sol)
- Z7 y  O' }) znumv=size(sol,2);
: h; S8 `+ n4 o: H3 ^7 [' Yx=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)
! d. ?# i) F% x0 I% Snumv=size(sol,2)-1;
3 ]- l) ?# y9 {' W; Px=sol(1:numv);
+ n9 b, b  x% x8 Leval=f(x);
" i% e% q$ J* Y6 o% R/ G+ @eval=-eval;
%主程序
$ v0 g7 H: V) L3 hclc,clear
- i& C9 S; @% G; N+ A, Kbounds=ones(2,1)*[-5 5];
initPop=initializega(10,bounds,'fitness');
0 q- t  P1 Z% w; R2 I# m7 q7 c[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.

Error in b2f (line 34)
  a=bval((cs(i)+1):cs(i+1));
, v% g0 O3 {: c% x

9 a  s+ c7 C7 Q) [Error in ga (line 194)
      bPop(bFoundIn,=[gen b2f(startPop(bindx,1:numVar),bounds,bits)...
; J# u1 N5 t" J8 x* b# D
+ ~9 J; h& e7 C# k+ k
. E% K9 ?% M, _, ^3 z1 w4 LError in Untitled4 (line 4)
[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]);
$ z, v9 Y' N! ~& j, Y3 ?" s8 f是我的工具箱有问题还是程序有问题，不懂求教大神 遗传工具箱.rar (102.67 KB, 下载次数: 0)

2013-8-23 10:52 上传

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

工具箱