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

s20090281

我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
7 b8 m' o: E+ t8 {* d5 n我要做的是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最小，
6 M5 f4 j6 V- @/ w1 e, s4 ]# s  q; Sfunction [out]=fitness(x)
     ai=x(i,           %第几组数据就是取第几个
) J1 p# e8 c) y! l    a1=ai(1)
    a2=ai(2)
    a3=ai(3)
% K& Z, V, R# M6 A    a4=ai(4)
: t/ {- A! ?  |    u=load('F:\测试记录控制.txt');u1=u(2:399,2);
: }' C$ i3 _9 ?& a* v9 z    y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
    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        %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度        
1 {) y* T+ P+ _/ t$ L8 o7 X        c=c+f
2 k0 u" N  U/ W7 ^* f" Z& c4 e        Msum(k)=c
- x4 z) p# h( ^: X  `    end
. I) Q+ Q% F% c) W; zout=Msum(:,k);                %输出适应度
$ d$ O5 A9 K. e: a3 \那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
1 v( m) N: a0 _2 s- k+ dfitness(i)=fitness(pop(i,)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。

s20090281

代表：）

s20090281

问题有新的变化，就是只能进行第一个运算，当计算第二个输入时就会出现??? Subscript indices must either be real positive integers or logicals.
. u$ B  ?$ S& B/ }* |) P: d. J
3 V8 W: K  J$ T& b- a) |$ TError in ==> myPSO at 14
    fitness(i)=fitness(pop(i,)
/ Q2 h# _% G# Y+ Q6 j5 B8 a0 E; m请问这怎么改呢？
& M+ T9 G% C- |+ v) V% U