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

function main()- _6 x( k" ?. U" ?& C1 z
clc;clear all;close all;
! M' Q' J" s: M
tic; %程序运行计时
! _2 ] I1 m. B( ?0 `
E0=0.001; %允许误差4 _1 C2 @* V2 [( ^( x) X1 y9 A' x
MaxNum=100; %粒子最大迭代次数1 m& l- H, g! A0 V
narvs=1; %目标函数的自变量个数( { v) P5 A! H' p
particlesize=30; %粒子群规模) {5 J M% _1 N& Y
c1=2; %每个粒子的个体学习因子，也称为加速常数
+ T) e3 }. U( z$ o& P
c2=2; %每个粒子的社会学习因子，也称为加速常数
5 Z$ H3 P- t {( H5 U8 h! x
w=0.6; %惯性因子% n; I# t8 q3 @# R4 {
vmax=0.8; %粒子的最大飞翔速度$ z& Y5 ?" P% W1 F
<font color="#ff0000">x=-2+4*rand(particlesize,narvs); %粒子所在的位置</font>: {8 J; o' s! ~6 I9 c- Y2 h5 H
v=2*rand(particlesize,narvs); %粒子的飞翔速度5 j. ]1 y- w; T( F, \
%用inline定义适应度函数以便将子函数文件与主程序文件放在一起，0 M7 }/ f/ f$ Z; u; Q! Y0 D
%目标函数是：y=200*exp(-0.05*x)*sin(x)
* L" C# b' P& l5 P; J, ~9 h, z
%inline命令定义适应度函数如下：% m" d/ y7 h& K* v3 i
fitness=inline('1/(200*exp(-0.05*x)*sin(x))','x');
# K( B" Y* i* J5 f! J0 O
%inline定义的适应度函数会使程序运行速度大大降低
/ R2 N7 m( X# K t- P+ S
for i=1:particlesize/ V0 V V( t3 J3 B9 O4 ?
for j=1:narvs
3 k6 x$ X z5 K; P# {) [
f(i)=fitness(x(i,j));+ t# Y. B9 G8 L2 q
end
Y1 G) }+ A& y
end+ [% p& L' ]2 q V2 f1 H7 w5 X+ d
personalbest_x=x;, O$ B7 [+ B0 V) w7 g+ g6 G/ f
personalbest_faval=f;; u8 I$ [( z0 Y+ i1 g1 @
[globalbest_faval i]=min(personalbest_faval);
1 X. J, B6 t& J- i& C
globalbest_x=personalbest_x(i,:);( |9 v; n, d* b0 m
k=1;$ }& B7 G; ], a& v/ t5 r
while k<=MaxNum
* n: J9 o( C' Q6 a7 m. {% K9 ~
for i=1:particlesize
! B) p' _; L5 q8 N$ z; O
for j=1:narvs
% h: a, R( G8 f. w
f(i)=fitness(x(i,j));
! Y9 a4 d% O% k
end" f+ e7 }9 q/ R
if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
) E' Y; L0 r; E3 U6 ~$ ~8 @
personalbest_faval(i)=f(i);& \9 h2 x+ i3 ?1 Q
personalbest_x(i,:)=x(i,:);
( F" c% {; h, g
end5 ^/ m$ [, C/ g- t& @
end+ e) X. b6 W3 M; u( f0 q
[globalbest_faval i]=min(personalbest_faval);5 u# R: ]+ V$ _7 G- |; P" _/ K
globalbest_x=personalbest_x(i,:);* e* ~# E& ]5 ~, r
for i=1:particlesize %更新粒子群里每个个体的最新位置
4 P7 @: G8 U8 G: m
v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...' i/ t" }7 R' k. n, p
+c2*rand*(globalbest_x-x(i,:));/ _" {- D1 H, C1 B1 [/ }$ J
for j=1:narvs %判断粒子的飞翔速度是否超过了最大飞翔速度
+ @$ \5 N7 `+ `1 b- S3 O
if v(i,j)>vmax;+ }# @9 T: C) k4 H% C1 ^
v(i,j)=vmax;
% c. o2 z, R/ a3 S$ T
elseif v(i,j)<-vmax;
' T/ C1 Z, D2 y q( W! H1 l0 c
v(i,j)=-vmax;! b$ M5 L8 L( I& k( c; e
end
) R) d# r5 T4 j3 D6 A% L" ^
end
4 y- |4 L, I* j. |. v
x(i,:)=x(i,:)+v(i,:);
. f7 _4 o. Q! j
end- ]- U* r7 s# K2 z& b# e
if abs(globalbest_faval)<E0,break,end. f) C; r7 b: x8 Q. }. F" J* E7 w* l
k=k+1;2 i* |- B/ `+ w8 ~1 E. p+ e
end; X- C( D4 n/ o! T8 `! A3 V
Value1=1/globalbest_faval; Value1=num2str(Value1);
7 @* P/ Z, T$ Q+ `: M: J4 v6 j
% strcat指令可以实现字符的组合输出
( x4 `- N) f" i3 X2 P
disp(strcat('the maximum value','=',Value1));! [; L1 Q# {2 i) s# H
%输出最大值所在的横坐标位置
7 I, y" n) I: N& V
Value2=globalbest_x; Value2=num2str(Value2);2 p7 C M# T c1 Q0 O# f& V
disp(strcat('the corresponding coordinate','=',Value2));! c z! P: y8 C$ S& r
x=-2:0.01:2;
$ E5 `( Z5 k5 A! X1 f5 L& V5 b8 Z# e
y=200*exp(-0.05.*x).*sin(x);/ G% c. {; ^& d5 I6 e4 [
plot(x,y,'m-','linewidth',3);4 `6 R/ p6 z0 d; j, B# h$ K5 E
hold on;5 D) P$ O( }% D3 u5 O1 F
plot(globalbest_x,1/globalbest_faval,'kp','linewidth',4);6 l* P' N/ e& m" ?
legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
* I% L7 R6 u* I$ ^0 w8 Y

复制代码