数学建模社区-数学中国
标题:
粒子群算法中适应度函数的调入
[打印本页]
作者:
s20090281
时间:
2011-9-20 08:51
标题:
粒子群算法中适应度函数的调入
我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
1 O" F4 W! [3 ^2 m2 W; R! Q" b
我要做的是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最小，
0 |3 i' ~# }9 N
function [out]=fitness(x)
9 y1 L5 M. u. m& j; B# G: h4 ]# e
ai=x(i,
%第几组数据就是取第几个
5 w+ [1 B: F0 l. {( \) k8 T; a* W/ ^
a1=ai(1)
; R- r& _5 s/ _
a2=ai(2)
: |, x& P+ B3 q; D" ^3 H0 e
a3=ai(3)
9 A& p! Y" `+ F' v
a4=ai(4)
7 e7 [6 |7 L8 D$ o) }5 k% Q$ i
u=load('F:\测试记录控制.txt');u1=u(2:399,2);
: R( X7 _# e7 ^' k% `7 d
y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
% m0 g& h' e; {& I5 I( C
c=0;
2 k3 ~. B( r) B& R
for k=3:399
+ W- J; H+ ]8 _1 K: O |
f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2 %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度
: B0 B! [' }- G( v2 A y
c=c+f
7 T; T% G: F- \* I# t2 o3 {. @
Msum(k)=c
: \3 G, Q, A8 L1 }% o% M# G
end
s; f, D5 n4 G, @$ k, x8 q0 g" L
out=Msum(:,k); %输出适应度
0 S3 z% I4 Y# N6 _4 _" N
那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
! p ^/ ~5 V& p/ g5 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.
3 f6 v1 F- z D* Q" Y' J& M
5 J8 L! p: m, @/ O' `
Error in ==> myPSO at 14
1 x! { @" N, |' F0 e
fitness(i)=fitness(pop(i,
)
; s$ o3 {6 d5 c! ?7 D& a$ ?4 o& }
请问这怎么改呢？
7 i* b5 V6 B7 @
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
%第几组数据就是取第几个)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。
代表：）); s$ o3 {6 d5 c! ?7 D& a$ ?4 o& }