查看: 5393|回复: 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()2 r9 W) r: \, N$ w' K! u
clc;clear all;close all;
8 a% \! U+ ~ R g- l, M
tic; %程序运行计时
# D0 X: K* x' `
E0=0.001; %允许误差! f; X: P% d! A1 P- j T: `. c8 g [
MaxNum=100; %粒子最大迭代次数
+ Q, @: K2 d& K( |0 e* r: ^
narvs=1; %目标函数的自变量个数
a& J; f6 \5 I6 w {# h; ]' r
particlesize=30; %粒子群规模/ e. t9 g6 m! y; E: K\" W
c1=2; %每个粒子的个体学习因子，也称为加速常数
( N% u6 K5 k* k& r) q& \
c2=2; %每个粒子的社会学习因子，也称为加速常数( p7 R1 B* i5 [4 D
w=0.6; %惯性因子* s {+ f6 o: j. |1 ]
vmax=0.8; %粒子的最大飞翔速度6 K: ^- R( ?1 v' B: O1 \
<font color="#ff0000">x=-2+4*rand(particlesize,narvs); %粒子所在的位置</font>
7 z$ c3 E$ v( a
v=2*rand(particlesize,narvs); %粒子的飞翔速度6 k0 E: ?2 L9 y7 f
%用inline定义适应度函数以便将子函数文件与主程序文件放在一起，3 ^4 @9 U6 W: G, m7 T2 K9 \
%目标函数是：y=200*exp(-0.05*x)*sin(x): g. V0 Z, B4 Z
%inline命令定义适应度函数如下：
7 v, I, V% H7 C\" g' `; G
fitness=inline('1/(200*exp(-0.05*x)*sin(x))','x');, Z4 z: j9 P2 g3 o
%inline定义的适应度函数会使程序运行速度大大降低
. p' e; ^4 @* e0 C, W, {7 B* W
for i=1:particlesize
\" ^2 ^( h- j7 X y
for j=1:narvs) D: x8 p- X. K6 ?: A
f(i)=fitness(x(i,j));
/ L% X: g7 m4 Z7 ]% p0 o
end/ R$ d7 z9 Q; r; U2 e/ s
end
personalbest_x=x;' C) A+ O7 g( u& @
personalbest_faval=f;
( q; L ?$ H. X( u& X\" ?8 a
[globalbest_faval i]=min(personalbest_faval);4 B- k O# O' f; |
globalbest_x=personalbest_x(i,:);& S7 ~3 @3 m. k) G
k=1;5 G9 l$ m' T# o6 U& q
while k<=MaxNum7 X. j- N) B E; Y | R0 \, W
for i=1:particlesize
( s. w& t* \. n! T% Y/ J, g
for j=1:narvs
, P; |* M( b* o8 l3 s* ?: n/ L
f(i)=fitness(x(i,j));
3 M' D* g/ Z5 i' j
end' a# U3 @, x\" _8 d5 }% G\" Y& Z
if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
/ @8 K9 s8 W6 t9 K/ b* u$ s
personalbest_faval(i)=f(i);\" A; N: K: y3 _ e
personalbest_x(i,:)=x(i,:);6 x& @0 n* X) d8 H. M5 ]* ~ l
end
& q) ]% N! S\" F3 M0 | e8 S/ d
end- I' ]4 p- A6 N8 z
[globalbest_faval i]=min(personalbest_faval);
) V; x# a- p; l: ~+ x/ j, F6 h( h
globalbest_x=personalbest_x(i,:);# q6 C; r% o+ L4 F* ]: c( |. j
for i=1:particlesize %更新粒子群里每个个体的最新位置
; S2 [$ \, ~) N& e6 r
v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
5 u\" Z2 M# j2 ~ r% r2 m
+c2*rand*(globalbest_x-x(i,:));
B' H0 x1 l8 U, b4 A- F
for j=1:narvs %判断粒子的飞翔速度是否超过了最大飞翔速度
' P) W9 t' s. D7 G' [\" C- M; s
if v(i,j)>vmax;9 r\" |1 V- G4 u
v(i,j)=vmax;& q. B- W; ]9 M3 J
elseif v(i,j)<-vmax;/ @, |5 ^/ b9 x9 j/ V- a
v(i,j)=-vmax;
2 x2 y- B1 L, o/ V5 a& _
end
4 {% ?$ v* [6 l% e2 k
end
6 T y/ i% g0 H0 F3 W
x(i,:)=x(i,:)+v(i,:);, ]/ N# R7 x7 Q\" O- N
end
9 N# u* t# n! l- k
if abs(globalbest_faval)<E0,break,end; N+ H# w. e% y, L
k=k+1;
7 J) f4 t9 k6 Y; A1 P8 t6 `+ j R0 M
end
( H7 a5 C/ I+ i0 v3 H6 ?5 n
Value1=1/globalbest_faval; Value1=num2str(Value1);+ A. R, {! `0 J5 S8 m
% strcat指令可以实现字符的组合输出
6 o) [& C+ S/ S
disp(strcat('the maximum value','=',Value1));- X( M, p1 c8 v9 ^5 A& y( n' R
%输出最大值所在的横坐标位置* \7 E2 l5 t) K
Value2=globalbest_x; Value2=num2str(Value2);/ x, d: d9 }3 p2 L
disp(strcat('the corresponding coordinate','=',Value2));1 b1 P9 `\" _3 ?* N\" w
x=-2:0.01:2;& L+ z( S, C. p\" ]: D* G; o @
y=200*exp(-0.05.*x).*sin(x);
' k/ N- C7 Q# Z0 U
plot(x,y,'m-','linewidth',3);9 q3 J) ?2 G3 A
hold on;
9 {. F, \# a9 L3 G' M( T7 @\" e
plot(globalbest_x,1/globalbest_faval,'kp','linewidth',4);
/ W# A: b2 Z8 X, s: l
legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
( @2 h: B$ I0 o; d: K# W: y( T* ~2 v

复制代码

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

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

zan