粒子群算法中适应度函数的调入

s20090281

我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
我要做的是y(k)=a1*y(k-1)+b1*u(k-1)+a2*y(k-2)+a3*y(k-3)函数求a1,a2,a3,b1.已知y和u，我把适应度取为399组数据的【y(k)-a1*y(k-1)-b1*u(k-1)-a2*y(k-2)-a3*y(k-3)】^2最小，
$ y3 N; H; n% H9 P4 G, U5 q& qfunction [out]=fitness(x)
3 ?- ?# S' K4 z& [) ^) S     ai=x(i,           %第几组数据就是取第几个
1 G0 V0 W' O% C" r; x' m. e% }    a1=ai(1)
4 Z  `: e& Q2 x  k9 a# j    a2=ai(2)
    a3=ai(3)
7 n( V5 n8 X$ e) [    a4=ai(4)
2 }0 g6 Y1 i6 k2 ~- w3 T    u=load('F:\测试记录控制.txt');u1=u(2:399,2);
, G& ~! r0 J. F5 q) g8 e% _    y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
  v, W0 @; ~$ k" }    c=0;
    for k=3:399
# c& t+ _1 S5 W/ n        f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2        %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度        
        c=c+f
( c1 L. a" t' L$ K7 \        Msum(k)=c
9 ?4 |" v( M* e5 W  y# k" r& K$ [9 F    end
3 F3 U; W3 O) X) t7 v. f. [, pout=Msum(:,k);                %输出适应度
3 g* T& R, I" Z) l7 n那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
- j1 ^( o& ~  wfitness(i)=fitness(pop(i,)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。

s20090281

代表：）

s20090281

问题有新的变化，就是只能进行第一个运算，当计算第二个输入时就会出现??? Subscript indices must either be real positive integers or logicals.
& z( v: e! D+ `; y, m
Error in ==> myPSO at 14
    fitness(i)=fitness(pop(i,)
5 R  x& l, r9 V: y请问这怎么改呢？