数学建模社区-数学中国
标题:
粒子群算法中适应度函数的调入
[打印本页]
作者:
s20090281
时间:
2011-9-20 08:51
标题:
粒子群算法中适应度函数的调入
我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
) a" w- J& T* p2 ]2 u9 C% o& H
我要做的是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最小，
g, l8 y2 N3 c% _6 p) \8 _
function [out]=fitness(x)
8 v- q$ I9 F, T
ai=x(i,
%第几组数据就是取第几个
. t; V8 [1 n1 j. `% f* d: X' M8 l
a1=ai(1)
$ D& F' S& R& }- R U
a2=ai(2)
( e- ?3 Q8 _0 B7 v' F* t8 W
a3=ai(3)
7 s1 O" O n( l0 \! V6 r
a4=ai(4)
/ n c, e) W m! z2 ?3 E5 y! f* m
u=load('F:\测试记录控制.txt');u1=u(2:399,2);
: ^2 O% z5 T6 U% P% S6 r4 L1 v
y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
; u) A' M: i7 ` l' R8 N# ]
c=0;
1 S* y) A" x8 k; K0 v% N
for k=3:399
# @5 Q; L$ F- x/ G# s
f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2 %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度
( y6 w" ~. v9 t
c=c+f
. P) }# u* M1 ~9 s, y0 G
Msum(k)=c
& N$ a4 F. S; d- Q
end
/ a6 R* ^' g' w, Z2 i; o
out=Msum(:,k); %输出适应度
2 \0 S7 m% _4 |# g+ ^2 F
那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
+ E# [/ [! _" a* T
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.
# e n2 e/ {$ T# g8 S2 n% Q
( ]' [2 h1 @) N4 R, c# Y9 z& M
Error in ==> myPSO at 14
+ y4 `8 {5 ~4 j% m! `. W6 j
fitness(i)=fitness(pop(i,
)
8 n. T0 z* i( g% w7 p% t
请问这怎么改呢？
) }& t. [' i8 Q" [, d
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
%第几组数据就是取第几个)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。
代表：）)