数学建模社区-数学中国

标题: 关于pso的问题 [打印本页]

作者: s20090281    时间: 2011-9-5 11:28
标题: 关于pso的问题
我想根据已知的y和u,利用式子y(k)=a1*y(k-1)+b1*u(k-1)+a2*y(k-2)+b2*u(k-2)+a3*y(k-3)求a1,a2,a3,b1,b2。我拟采用PSO来做,但是刚开始学这个算法,我想问问我的想法对不对,请各位高手给点。
; J0 j9 c: x6 P* E( w, g1、这是不是要用离散PSO算法;7 X# `, {: h- Q# g: }
2、首先应该确定微粒数以及维数。我想微粒数不用太注意,关键是维数的确定,那么维数是不是应该是这样的,是17维,分别为a1,a2,a3,b1,b2的值,a1,a2,a3,b1,b2的速度,a1,a2,a3,b1,b2的位置坐标,以及个体最优适应度值和当前适应度值;' k7 W  z# O+ k, b9 c& x. K
3、适应度的确定:在这我犯愁了该怎么办呢?莫非是把好多组数据代入求误差平方和的最小值,或者有什么简单的方法吗?
作者: syxz2007    时间: 2011-9-7 17:05
1、已知的y和u具体是什么样子的
/ k8 Q3 w7 }, H: @& P7 |! t2、维数一般是你要求解的变量的个数
+ t# \& M2 ~9 |- u6 B% z3、适应度要根据具体问题来定,对于你的问题还不是太清楚
作者: shuxuezaozhuang    时间: 2011-9-8 10:16
不太明白!!
作者: s20090281    时间: 2011-9-13 08:22
syxz2007 发表于 2011-9-7 17:05 # J  n7 u5 f, c* w. T% Y2 K
1、已知的y和u具体是什么样子的8 F  `& m  ^6 p+ r- J
2、维数一般是你要求解的变量的个数$ V3 V3 M+ B+ Q" c2 O
3、适应度要根据具体问题来定,对于你 ...
% J" u2 c$ Z* J+ W9 L
1、y(k)=a1*y(k-1)+b1*u(k-1)+a2*y(k-2)+b2*u(k-2)+a3*y(k-3),画出图来,y和u连接起来是正弦信号,那相当于是连续的还是离散的呢?y是[-4.62,-4.65,-4.79,-4.62,-4.78.-4.73,-4.78,-4.69.-4.84,-4.86...]等这样的数值,连接起来是正弦信号,u是[4.77,4.76,4.82,4.83,4.85,4.86,4.87,4.88,4.9...]连接起来是正弦信号.
" E# P8 b4 ~5 c. n- E! Z# e5 P2、这个问题清楚了。
0 f. c9 _6 D) s0 Z3、我想做的事情就是已知u和y有很多组数据,求得a1,b1,a2,b2,a3,使得代入每一组数据等式两边都尽可能的相等,所以我把适应度函数定义为sum[y(k)-a1*y(k-1)+b1*u(k-1)+a2*y(k-2)+b2*u(k-2)+a3*y(k-3)]^2行不行。( q  _$ P: q" ?) x' T6 s
希望能得你的回复,谢谢。
3 N4 k8 y+ P  p+ O/ @5 e4 h) L3 X希望能尽快得到回复,谢谢。
作者: xiaderen    时间: 2011-9-15 22:39
其实你是想求取序列y和u之间的关系表达式,已经有了关于Y和U的一些对应数据,求取表达式中的系数。
* @; _0 P' l8 L; F' g% wsum([y(k)-a1*y(k-1)-b1*u(k-1)-a2*y(k-2)-b2*u(k-2)-a3*y(k-3)]^2),作为适应度函数的定义,+ N/ k/ a3 M. f- {
维数应该是5维的。
作者: alair006    时间: 2012-2-7 15:40
顶顶更健康,越顶吃的越香。5735195334418082
作者: setou    时间: 2012-4-7 16:00
怎么办 还是看不懂 先顶一个
作者: 啊啊啊木    时间: 2012-5-24 19:12
顶了吧~~~~~~~
作者: Μene∫Aひζ    时间: 2012-5-26 00:16
..........h
作者: liwenhui    时间: 2012-6-27 12:34
这个是参数估计问题,采用最大似然估计法的话,划为一个最优化问题。不一定非得用PSO。




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5