大家帮忙看看，为什么我的粒子群算法PSO得到的结果不正确
本帖最后由 doilovematlab 于 2014-1-25 20:20 编辑-------------------------------------------------update at 2014-1-25----------------------------------------------
下面是我上传的附件，pso程序：
--------------------------------分割线------------------------------------------
我想求 http://www.ilovematlab.cn/data/attachment/forum/201401/22/155754jl4lgqffn5cfd4z1.jpg 此函数在【-2,2】的最大值，使用如下的粒子群算法求解：function main()
clc;clear all;close all;
tic; %程序运行计时
E0=0.001; %允许误差
MaxNum=100; %粒子最大迭代次数
narvs=1; %目标函数的自变量个数
particlesize=30; %粒子群规模
c1=2; %每个粒子的个体学习因子，也称为加速常数
c2=2; %每个粒子的社会学习因子，也称为加速常数
w=0.6; %惯性因子
vmax=0.8; %粒子的最大飞翔速度
<font color="#ff0000">x=-2+4*rand(particlesize,narvs); %粒子所在的位置</font>
v=2*rand(particlesize,narvs); %粒子的飞翔速度
%用inline定义适应度函数以便将子函数文件与主程序文件放在一起，
%目标函数是：y=200*exp(-0.05*x)*sin(x)
%inline命令定义适应度函数如下：
fitness=inline('1/(200*exp(-0.05*x)*sin(x))','x');
%inline定义的适应度函数会使程序运行速度大大降低
for i=1:particlesize
for j=1:narvs
f(i)=fitness(x(i,j));
end
end
personalbest_x=x;
personalbest_faval=f;
=min(personalbest_faval);
globalbest_x=personalbest_x(i,:);
k=1;
while k<=MaxNum
for i=1:particlesize
for j=1:narvs
f(i)=fitness(x(i,j));
end
if f(i)<personalbest_faval(i) %判断当前位置是否是历史上最佳位置
personalbest_faval(i)=f(i);
personalbest_x(i,:)=x(i,:);
end
end
=min(personalbest_faval);
globalbest_x=personalbest_x(i,:);
for i=1:particlesize %更新粒子群里每个个体的最新位置
v(i,:)=w*v(i,:)+c1*rand*(personalbest_x(i,:)-x(i,:))...
+c2*rand*(globalbest_x-x(i,:));
for j=1:narvs %判断粒子的飞翔速度是否超过了最大飞翔速度
if v(i,j)>vmax;
v(i,j)=vmax;
elseif v(i,j)<-vmax;
v(i,j)=-vmax;
end
end
x(i,:)=x(i,:)+v(i,:);
end
if abs(globalbest_faval)<E0,break,end
k=k+1;
end
Value1=1/globalbest_faval; Value1=num2str(Value1);
% strcat指令可以实现字符的组合输出
disp(strcat('the maximum value','=',Value1));
%输出最大值所在的横坐标位置
Value2=globalbest_x; Value2=num2str(Value2);
disp(strcat('the corresponding coordinate','=',Value2));
x=-2:0.01:2;
y=200*exp(-0.05.*x).*sin(x);
plot(x,y,'m-','linewidth',3);
hold on;
plot(globalbest_x,1/globalbest_faval,'kp','linewidth',4);
legend('目标函数','搜索到的最大值');xlabel('x');ylabel('y');grid on;toc;
但是得到的结果却是这样的：
http://www.ilovematlab.cn/data/attachment/forum/201401/22/155933n2ey56le1wd1ms5x.jpg
而且当我调节粒子群的初始位置的时候，得到的结果相差很大，我想知道设置粒子群的初始位置：
x=-2+4*rand(particlesize,narvs); %粒子所在的位置
这句代码是否正确，是用x的最小值-2加上随机矩阵乘以x的总范围4吗？如果是或者不是，我到底哪里出错了呢？
这么多程序，！！！！！ {:3_42:}{:3_42:}{:3_42:} 给个文件，这里我复制后就乱码，不识别 薛云涛 发表于 2014-1-24 09:48 static/image/common/back.gif
给个文件，这里我复制后就乱码，不识别
不好意思，现在才上网，附件已经上传，麻烦帮忙看下。谢谢。 PSO好难呀…… 好难的感觉啊
页:
[1]