数学建模社区-数学中国
标题:
粒子群算法中适应度函数的调入
[打印本页]
作者:
s20090281
时间:
2011-9-20 08:51
标题:
粒子群算法中适应度函数的调入
我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
2 h- M8 @8 X8 j
我要做的是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最小，
9 s* N( Z* u5 U+ K% V
function [out]=fitness(x)
# F" u5 {$ |! S% ^5 i
ai=x(i,
%第几组数据就是取第几个
6 `% d' y+ O9 R% D
a1=ai(1)
: f' O% }, d1 N( @& e$ [& i
a2=ai(2)
4 l) Q; [3 d( j- S
a3=ai(3)
( j/ A( v# E$ U8 F# [
a4=ai(4)
3 D( I* L9 r& p+ w
u=load('F:\测试记录控制.txt');u1=u(2:399,2);
0 j0 q% j8 V% k6 l# z' E& x/ M: @
y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
8 }9 m' @# ?4 c
c=0;
' d) X- r- {' s0 D; S; N3 e
for k=3:399
: u1 v* b+ g V- ^/ A
f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2 %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度
U/ o" R0 z- | q
c=c+f
& U/ I+ J8 |# l4 U* X
Msum(k)=c
; V* z1 E! M# Q/ A- H7 O
end
2 p9 R% s& y( s) ^! t) p
out=Msum(:,k); %输出适应度
+ i0 }4 y+ ~! a l5 Y3 i
那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
# d6 j! M: {" I0 z* \8 T$ R/ Q
fitness(i)=fitness(pop(i,
)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。
作者:
s20090281
时间:
2011-9-20 08:51
代表：）
作者:
s20090281
时间:
2011-9-20 16:16
问题有新的变化，就是只能进行第一个运算，当计算第二个输入时就会出现??? Subscript indices must either be real positive integers or logicals.
1 k+ N8 _! Q1 G! P+ @
$ q) a9 c! q" e
Error in ==> myPSO at 14
3 F# w! O. s7 H% N) }( J3 q4 I
fitness(i)=fitness(pop(i,
)
u7 {0 w4 f0 G0 j6 M2 ` z
请问这怎么改呢？
0 O6 D5 K2 \/ [, d% I2 y* |) v" D+ G
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
%第几组数据就是取第几个)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。
代表：）)