大家帮忙看看，为什么我的粒子群算法PSO得到的结果不正确 - 数学建模社区-数学中国

function main()
) ?0 K+ m! k( l+ x7 w3 V
clc;clear all;close all;! [5 u k7 m% Y9 D t
tic; %程序运行计时 i5 n* D& y2 J0 B
E0=0.001; %允许误差
: ? M) e' O. X
MaxNum=100; %粒子最大迭代次数. T% c" }6 k0 B. Y) C
narvs=1; %目标函数的自变量个数 q2 C0 c! |$ e3 T7 {) M
particlesize=30; %粒子群规模
) u8 Q0 i {# _' M: T9 H
c1=2; %每个粒子的个体学习因子，也称为加速常数
2 z# I) {1 R1 J; Z. n* D
c2=2; %每个粒子的社会学习因子，也称为加速常数. H6 g: p) e A/ c. J
w=0.6; %惯性因子4 D2 b# v! L5 @# c5 L8 o# S+ ]% q
vmax=0.8; %粒子的最大飞翔速度
4 s; p( ]/ N4 U) } C2 ?
<font color="#ff0000">x=-2+4*rand(particlesize,narvs); %粒子所在的位置</font>
2 ~9 D: }7 W: G. u
v=2*rand(particlesize,narvs); %粒子的飞翔速度" ?$ E) x! t J& \! l8 p
%用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
! q) n$ ^5 Z/ S
%目标函数是：y=200*exp(-0.05*x)*sin(x)
3 g9 m i* J1 {0 t
%inline命令定义适应度函数如下：3 N T7 O3 J+ u5 [* }
fitness=inline('1/(200*exp(-0.05*x)*sin(x))','x');
! ?6 z+ L( ^- P4 l* _
%inline定义的适应度函数会使程序运行速度大大降低
" F' a; z% J. l! B0 l
for i=1:particlesize
: t) Q( N3 O0 M$ N
for j=1:narvs
" P8 n; z' M: B3 o/ Q/ Y
f(i)=fitness(x(i,j));9 W# I \& L4 k/ h6 }$ \$ V
end
) d* ]2 z- z) [) p2 h7 Z0 L
end
* C9 K$ V0 d4 l7 e. Z- U
personalbest_x=x;
8 f l! \8 b( o0 M) Z* R; E# K0 ^
personalbest_faval=f;
2 C2 r2 H' j! U; e8 A9 b6 P
[globalbest_faval i]=min(personalbest_faval);
p9 E0 E9 n4 ~, J# h
globalbest_x=personalbest_x(i,:);3 n* ?5 u, @: a+ ^
k=1;/ l0 A N o7 a7 J7 O
while k<=MaxNum$ o {1 k: g8 G, x
for i=1:particlesize
$ A3 ]1 f) ~% t6 q& g
for j=1:narvs
+ U9 |: m$ M) C3 T. R3 e' e5 v
f(i)=fitness(x(i,j));
1 f0 x6 o# N6 X" ~" {9 |7 h
end
6 B6 {4 Y$ L& w. z4 u% K
if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
/ L p( u' W8 E) l' y1 h
personalbest_faval(i)=f(i);
}9 y. I6 s! I( S5 A
personalbest_x(i,:)=x(i,:);# |6 s2 z. E6 Y7 ]8 [6 k
end) [1 K3 T2 N" z1 Z2 ?
end K5 B2 |4 O- J, E( S7 N9 ~' m
[globalbest_faval i]=min(personalbest_faval);; f7 }3 K4 v0 h: a/ I) N6 t
globalbest_x=personalbest_x(i,:);
' l8 B+ i( p+ m0 C2 B. C& S
for i=1:particlesize %更新粒子群里每个个体的最新位置/ Z, j9 h) m0 E; [& k0 U
v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
: r0 `9 K* y3 S: A- l0 T% g
+c2*rand*(globalbest_x-x(i,:));' D. Q* Y5 s, q. \! G) B7 ~ |
for j=1:narvs %判断粒子的飞翔速度是否超过了最大飞翔速度
) t7 W2 E) y6 s
if v(i,j)>vmax;5 N% y( s' Q( G9 v2 e
v(i,j)=vmax;
5 W" A/ F' D) W
elseif v(i,j)<-vmax;, M# T7 m8 B9 k# s1 G
v(i,j)=-vmax;5 {/ }0 P( |, c7 \8 p/ |4 Q
end
end4 V9 b' u) P b( I- }% E$ ^; ]* R
x(i,:)=x(i,:)+v(i,:);
! Y* r& Q- |; v6 Y$ ?
end c1 c; B8 H/ \: T/ A
if abs(globalbest_faval)<E0,break,end
. p6 w: D) L2 {$ N$ W. L
k=k+1;1 p2 i9 l# |" v) r3 Z2 \
end
1 }* s% s9 \& y. y; E1 M; F! X$ f
Value1=1/globalbest_faval; Value1=num2str(Value1);5 `3 z$ [# C8 h/ D4 ]# L( H2 P/ R
% strcat指令可以实现字符的组合输出
4 L! s% c/ g1 ^; }) p0 O
disp(strcat('the maximum value','=',Value1));
% K y$ o, C ~; F* ~8 w
%输出最大值所在的横坐标位置0 m0 g7 i/ O$ w" A
Value2=globalbest_x; Value2=num2str(Value2);, H/ n) f( ^- X0 d6 C$ l
disp(strcat('the corresponding coordinate','=',Value2));
J; x" U8 E. e+ a d
x=-2:0.01:2;( T2 Z) _- _! W* K6 k8 }
y=200*exp(-0.05.*x).*sin(x);
1 J7 ^' x0 l: n
plot(x,y,'m-','linewidth',3);! w% g- f4 v3 n
hold on;
E4 F- {/ S% A- I6 y+ F
plot(globalbest_x,1/globalbest_faval,'kp','linewidth',4);: h% R* Z5 z0 [. c: L+ ^
legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;, J+ Q. s; S- v7 k

复制代码