数学建模社区-数学中国
标题:
粒子群算法中适应度函数的调入
[打印本页]
作者:
s20090281
时间:
2011-9-20 08:51
标题:
粒子群算法中适应度函数的调入
我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
5 y2 d% n. H1 H4 ~ Z% F2 B# F) R
我要做的是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最小，
! R3 R; q5 z; ` ^/ L; B3 Q- [
function [out]=fitness(x)
# _9 @, @. }7 F8 g
ai=x(i,
%第几组数据就是取第几个
( ?/ Q8 a) a6 t) M/ R5 S
a1=ai(1)
: g: e+ F5 j' s% a& }1 |, G6 @6 N
a2=ai(2)
! a0 T( B9 B) i# _3 C5 H# ^
a3=ai(3)
e5 k: i; _* b! D7 L& X
a4=ai(4)
$ f o8 j2 f/ ]7 f! G
u=load('F:\测试记录控制.txt');u1=u(2:399,2);
# s! o/ Y W( A0 j, x/ E
y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
( S: B+ F% N( L6 U' r8 L
c=0;
# k3 D5 G H" w" ^0 b' U
for k=3:399
" x1 b9 ~% m! { c
f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2 %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度
( S" c( H+ F) `
c=c+f
7 N8 l' E! ~* t# k# _; O. o, _
Msum(k)=c
8 y; t3 d& ]: R5 \& W
end
9 o! r6 q# _1 U7 I$ Z1 S
out=Msum(:,k); %输出适应度
8 X/ ^" b, r) g" M$ U$ E
那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
5 G5 O, ^' u8 {* d" ~
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.
) m- A, J1 b5 y; y$ t, Z/ J
8 H! S) }+ D: L
Error in ==> myPSO at 14
0 R4 h) R- d; N$ a* j S
fitness(i)=fitness(pop(i,
)
6 i4 l. D2 e( d& O- a- h
请问这怎么改呢？
5 J& y2 T. ^' u. @8 e
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
%第几组数据就是取第几个( ?/ Q8 a) a6 t) M/ R5 S)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。
代表：）)6 i4 l. D2 e( d& O- a- h