数学建模社区-数学中国
标题:
粒子群算法中适应度函数的调入
[打印本页]
作者:
s20090281
时间:
2011-9-20 08:51
标题:
粒子群算法中适应度函数的调入
我最近参考别人的程序写了一个自己的程序但是代入不了fitness函数，请大家指点一下：
6 {* ~# {' u0 \2 f
我要做的是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最小，
- w: G3 L1 Q) }) k" n
function [out]=fitness(x)
! `8 m& s5 s3 X" a$ p
ai=x(i,
%第几组数据就是取第几个
. q/ l& o6 U. O
a1=ai(1)
) n7 R$ U+ c! S
a2=ai(2)
4 @- l6 ^5 N' c
a3=ai(3)
' k. H; _$ [) z4 ?: c
a4=ai(4)
2 i. w7 P/ g$ a2 S! N. j
u=load('F:\测试记录控制.txt');u1=u(2:399,2);
& P6 m3 m/ u- I. [% @
y=load('F:\测试记录反馈.txt');y1=y(2:400,2);
7 b% g- @7 h: R V# N
c=0;
$ ], \% R7 v: u! K# y4 Y6 ?
for k=3:399
- j% E6 T* Z+ G6 `0 W
f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2 %求出各级数据的平方和并逐渐相加得到一组对对应变量的适应度
) S- V* d+ a# ]2 _* y
c=c+f
$ I$ ^" {* W) P5 m. }, K5 W
Msum(k)=c
; Z4 q# y3 E* k) |$ P
end
4 s! C3 F: W- o0 d8 w& s' Z
out=Msum(:,k); %输出适应度
}' C3 t/ e- L- ?' D
那么在粒子群算法中应如何调用这个函数呢？当位置和速度更新后，怎么计算这时的适应度呢？
7 e! |) ~) i7 S& F4 J! i
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 s x* T* A5 ^1 M3 s( V& O# K
4 e+ m& O4 Q5 d& ^+ S$ P+ m
Error in ==> myPSO at 14
* \" F% U+ H# y4 m& U+ v* R* W
fitness(i)=fitness(pop(i,
)
& A! ?2 P* l4 B3 I; F. z
请问这怎么改呢？
0 ^: e7 B7 f3 _
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
%第几组数据就是取第几个)这样调用为什么会出错呢？请各位帮我这个新手解答一下，谢谢。
代表：）)& A! ?2 P* l4 B3 I; F. z