查看: 5079|回复: 6

[问题求助] 大家帮忙看看，为什么我的粒子群算法PSO得到的结果不正确

字体大小: 正常放大

2 主题	8 听众	74 积分

升级 72.63%

TA的每日心情

	奋斗 2014-5-6 00:09

签到天数: 18 天

[LV.4]偶尔看看III

自我介绍: 没有自我介绍

电梯直达

1^#

发表于 2014-1-22 16:09 |只看该作者 |正序浏览

|招呼Ta 关注Ta

本帖最后由 doilovematlab 于 2014-1-25 20:20 编辑

-------------------------------------------------update at 2014-1-25----------------------------------------------
下面是我上传的附件，pso程序：

main.m (2.4 KB, 下载次数: 3)

--------------------------------分割线------------------------------------------
我想求此函数在【-2,2】的最大值，使用如下的粒子群算法求解：

function main()4 X5 H; v! o& S& q# P
clc;clear all;close all;\" x\" |4 {7 v/ ~9 s& ?; T u9 w
tic; %程序运行计时
+ } O, I/ C$ q- B0 l
E0=0.001; %允许误差( y, ~- L) z+ x1 S/ X\" ^* E$ |) c: w
MaxNum=100; %粒子最大迭代次数$ G, Y8 d- `8 ~7 A% H9 {; E
narvs=1; %目标函数的自变量个数6 p3 c/ H: e! @. l1 d6 p4 a8 v3 f& U
particlesize=30; %粒子群规模
* v' |7 ?5 ]5 B5 Y% i2 N& W
c1=2; %每个粒子的个体学习因子，也称为加速常数+ [8 u2 p! i, c, }
c2=2; %每个粒子的社会学习因子，也称为加速常数2 Q; P: `* [* h
w=0.6; %惯性因子& E* x9 E8 S; y+ s# J. T3 |
vmax=0.8; %粒子的最大飞翔速度
`6 M; |/ h r3 Q$ U/ V% c
<font color="#ff0000">x=-2+4*rand(particlesize,narvs); %粒子所在的位置</font>, {6 E+ w8 E$ l) _) w) m
v=2*rand(particlesize,narvs); %粒子的飞翔速度
$ A\" t; q: Z/ ?+ G( U4 H
%用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
5 i7 a) k. v, y\" Y# n) B
%目标函数是：y=200*exp(-0.05*x)*sin(x)1 J D, g+ m3 u) c
%inline命令定义适应度函数如下：
4 B1 B2 H5 ^/ j7 E\" F5 Q B! ~, i. Z
fitness=inline('1/(200*exp(-0.05*x)*sin(x))','x');
( ` d& W- `: T }6 j
%inline定义的适应度函数会使程序运行速度大大降低
; T; n B/ X- ^\" w8 @, J
for i=1:particlesize
; J$ H/ p% Z( P' ], w. R' @
for j=1:narvs
\" H% Q* U) E8 ~5 D
f(i)=fitness(x(i,j));$ E) G( S, t5 H9 R- Q\" E2 i
end |\" w: Z* Z8 n) c
end2 d% X\" A6 n) I! k* ~$ x$ `% |( e0 Z/ X
personalbest_x=x;
0 F4 h1 s1 D; ^2 k, ] q0 M; C
personalbest_faval=f;
$ U( i: ^3 O+ e5 q
[globalbest_faval i]=min(personalbest_faval);0 {; E3 T1 \) z. c! ?
globalbest_x=personalbest_x(i,:);: B# ]& A' d' i; p; L4 Y% H, b9 y
k=1;* y& L' j2 Q2 X( G X7 C
while k<=MaxNum
- D) n\" M9 S2 M1 x6 n7 i
for i=1:particlesize
\" O+ m. o1 d' ]
for j=1:narvs
3 T: _/ N9 I# ^\" f5 j, M) n; q2 c6 Z
f(i)=fitness(x(i,j));
) X) A! o, W' l# y
end9 M& {- ? Y\" G
if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置) g) w4 O7 S) h8 P G\" u4 {
personalbest_faval(i)=f(i);\" V/ e; s# @$ ?. d# p# c' q, e
personalbest_x(i,:)=x(i,:);
1 A& D; e! H3 G( _% S
end
( s3 A* Z& T- p+ G7 o# g
end
3 p5 [; S& d( V+ l1 J2 j- {
[globalbest_faval i]=min(personalbest_faval);* M+ y# B- r5 P6 N\" ~\" X0 o
globalbest_x=personalbest_x(i,:);& ]; B7 A6 W2 a/ K; z
for i=1:particlesize %更新粒子群里每个个体的最新位置
1 U8 \2 C/ m5 Y3 r
v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
+ p$ F. O) T8 Z7 K* s* \3 K
+c2*rand*(globalbest_x-x(i,:));
: g- x% f, q\" Z
for j=1:narvs %判断粒子的飞翔速度是否超过了最大飞翔速度
- t\" z0 k3 V7 [
if v(i,j)>vmax;
8 G9 T% r4 I6 W; W7 u9 T& B: a1 _0 t
v(i,j)=vmax;
/ W& x' B; R\" h3 A
elseif v(i,j)<-vmax;1 Q/ g) u- g. f
v(i,j)=-vmax;
& {$ u1 l w9 n' C+ \3 \
end
: k\" Z: n# T' }% g4 i
end1 C- g5 U/ [; D5 \1 x
x(i,:)=x(i,:)+v(i,:);! G( ?$ ~, [' @ ^3 A; J+ Y6 B. ~/ O
end
4 F+ N5 b: ~3 H. ?8 M/ U$ d
if abs(globalbest_faval)<E0,break,end
* {* `9 a# b\" l, i |# _, m
k=k+1;
6 ]- \9 d+ Z/ b' D
end5 p9 W* U' N5 F3 G) i5 ^; D) K
Value1=1/globalbest_faval; Value1=num2str(Value1);5 l, `; j# u- D- C4 _/ }: b
% strcat指令可以实现字符的组合输出* l7 @& l$ P2 m( P8 H
disp(strcat('the maximum value','=',Value1));4 \5 p+ V9 w1 ?; j4 h4 y6 v# d
%输出最大值所在的横坐标位置
3 E1 [& V$ }& y( c8 _
Value2=globalbest_x; Value2=num2str(Value2);
/ I, d9 k& t2 {4 A# s. D' I
disp(strcat('the corresponding coordinate','=',Value2));
2 D9 i ]2 C, A: p9 |& A- |\" G7 N
x=-2:0.01:2;
\" i! s: z3 D6 J: v0 {7 g
y=200*exp(-0.05.*x).*sin(x);/ J) Z/ k. F& \: N0 w
plot(x,y,'m-','linewidth',3);/ ^9 |' v, v i; T8 Z7 e
hold on;
- M8 F8 N, T4 ?* v: [6 B' t' ~
plot(globalbest_x,1/globalbest_faval,'kp','linewidth',4);# P4 B, i9 b! @
legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
! z) Z6 V; }( G* |

复制代码

但是得到的结果却是这样的：

而且当我调节粒子群的初始位置的时候，得到的结果相差很大，我想知道设置粒子群的初始位置：
x=-2+4*rand(particlesize,narvs); %粒子所在的位置
这句代码是否正确，是用x的最小值-2加上随机矩阵乘以x的总范围4吗？如果是或者不是，我到底哪里出错了呢？

zan