 [算法]基本粒子群优化算法Matlab源程序 |
基本粒子群优化算法Matlab源程序
这个程序就是最基本的粒子群优化算法程序，用Matlab实现，非常简单。只有几十行代码。正所谓一分钱一分货啊，优化效果不总是令人满意。我还有几个改进的粒子群优化算法版本，这一段时间会陆续发上来。
下面是主函数的源程序，优化函数则以m文件的形式放在fitness.m里面，对不同的优化函数只要修改fitness.m就可以了通用性很强。
主函数源程序（main.m）
%------基本粒子群优化算法（Particle Swarm Optimization）-----------
%------名称：基本粒子群优化算法（PSO）
%------作用：求解优化问题
%------说明：全局性，并行性，高效的群体智能算法
%------作者：孙明杰（dreamsun2001@126.com）
%------单位：中国矿业大学理学院计算数学硕2005
%------时间：2006年8月17日 <CopyRight@dReAmsUn>
%------------------------------------------------------------------
%------初始格式化--------------------------------------------------
clear all;
clc;
format long;
%------给定初始化条件----------------------------------------------
c1=1.4962; %学习因子1
c2=1.4962; %学习因子2
w=0.7298; %惯性权重
MaxDT=1000; %最大迭代次数
D=10; %搜索空间维数（未知数个数）
N=40; %初始化群体个体数目
eps=10^(-6); %设置精度(在已知最小值时候用)
%------初始化种群的个体(可以在这里限定位置和速度的范围)------------
for i=1:N
for j=1:D
x(i,j)=randn; %随机初始化位置
v(i,j)=randn; %随机初始化速度
end
end
%------先计算各个粒子的适应度，并初始化Pi和Pg----------------------
for i=1:N
p(i)=fitness(x(i,:),D);
y(i,:)=x(i,:);
end
pg=x(1,:); %Pg为全局最优
for i=2:N
if fitness(x(i,:),D)<fitness(pg,D)
pg=x(i,:);
end
end
%------进入主要循环，按照公式依次迭代，直到满足精度要求------------
for t=1:MaxDT
for i=1:N
v(i,:)=w*v(i,:)+c1*rand*(y(i,:)-x(i,:))+c2*rand*(pg-x(i,:));
x(i,:)=x(i,:)+v(i,:);
if fitness(x(i,:),D)<p(i)
p(i)=fitness(x(i,:),D);
y(i,:)=x(i,:);
end
if p(i)<fitness(pg,D)
pg=y(i,:);
end
end
Pbest(t)=fitness(pg,D);
end
%------最后给出计算结果
disp('*************************************************************')
disp('函数的全局最优位置为：')
Solution=pg'
disp('最后得到的优化极值为：')
Result=fitness(pg,D)
disp('*************************************************************')
%------算法结束---DreamSun GL & HF-----------------------------------
适应度函数源程序（fitness.m）
function result=fitness(x,D)
sum=0;
for i=1:D
sum=sum+x(i)^2;
end
result=sum;
下面几天回继续发上来其他改进算法，并发上比较结果。先看看我这有3种方法对此函数的优化结果比较图形,图形不是很清楚,一直下降的就是改进算法比较成功的
500)this.width=500'>
|
阅读全文(20223) | 回复(20) | 编辑 | 精华 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
lyn(游客)发表评论于2009/4/30 11:54:44 |
啊有联系方式 留一个啊
|
个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
pandeng(游客)发表评论于2009/4/25 10:45:13 |
我想要改进后的算法，邮箱为：864821649G@gmail.com
多谢楼主
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
软测量(游客)发表评论于2009/4/23 10:55:12 |
谢谢主人的分享，很希望能向你请教关于粒子群优化的MATLAB编程问题。
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
huaihai(游客)发表评论于2008/8/23 17:43:13 |
以下引用R(游客)在2008-2-29 16:18:30的评论:
赞人品请教楼主,我是新手,怎样才能运行你的程序
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
huang20021983(游客)发表评论于2008/4/21 16:16:57 |
搂主测试的是什么函数啊
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
datou30(游客)发表评论于2008/4/10 13:14:47 |
谢谢博主的程序 对我很有启发 你真是个好人
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
jinye123(游客)发表评论于2008/4/4 12:22:19 |
楼主 非常感谢 写的简洁明了 外行都能看懂
之前复制走了 改进了一下
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
R(游客)发表评论于2008/2/29 16:18:30 |
赞人品
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
飞行(游客)发表评论于2008/1/3 9:51:31 |
楼主真是好人呐
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
| 回复:基本粒子群优化算法Matlab源程序 |
|
乔(游客)发表评论于2007/12/5 9:29:54 |
谢谢博主的无私奉献精神，！
|
| 个人主页 | 引用回复 | 主人回复 | 返回 | 编辑 | 删除 |
|
|
| « | September 2025 | » | | 日 | 一 | 二 | 三 | 四 | 五 | 六 | | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 | | | | | |
|
| 公告 |
|
有一种鸟儿是永远关不住的
因为它的每片羽翼上都沾满了自由的光辉
方向:计算机视觉 人工智能 演化算法
| |
| Blog信息 |
|
blog名称:阳光海岸心
日志总数:166
评论数量:237
留言数量:-4
访问次数:1456516
建立时间:2006年6月2日 | |
|