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

s20090281

我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
7 R2 ?# K5 x" w8 ]" i7 k我要做的是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最小，
5 @2 I; j$ x* W6 N8 g! P: {function [out]=fitness(x)
     ai=x(i,           %第几组数据就是取第几个
1 |0 t/ S6 o' R. F  r# W    a1=ai(1)
# A  q8 x# ?* \+ H+ g3 `/ ]& S    a2=ai(2)
, ]3 }* m- [- s/ ~6 l& n    a3=ai(3)
1 e* C% p- s4 j  _6 c& X' |2 X2 o8 n    a4=ai(4)
    u=load('F:\测试记录控制.txt');u1=u(2:399,2);
    y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
9 B5 U5 w2 Q" U% N9 a9 e! G    c=0;
    for k=3:399
        f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2        %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度        
        c=c+f
. [1 g6 |8 _1 N8 L' m- U6 z% z        Msum(k)=c
* `# K0 j0 ^) |! s- V9 t$ i- ^    end
out=Msum(:,k);                %输出适应度
那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
/ r% ?3 {5 q2 p  f( M4 [5 a/ Hfitness(i)=fitness(pop(i,)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。

s20090281

代表：）

s20090281

问题有新的变化，就是只能进行第一个运算，当计算第二个输入时就会出现??? Subscript indices must either be real positive integers or logicals.
& o: X1 o+ E& f& n
  k# o( A! \* u% Q! t4 S8 UError in ==> myPSO at 14
5 U  `# N  D( s1 v* f9 {- H/ E    fitness(i)=fitness(pop(i,)
请问这怎么改呢？
- }$ X; B: M# B! l