




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、具有量子行为的粒子群优化算法惯性权重研究及应用粒子群优化(PSO)算法是一种群智能优化算法,最早由Kennedy和Eberhart于1995年共同提出,其基本思想是对鸟群捕食行为的仿生模拟,通过鸟群之间的集体协作,快速搜寻并找到最优解。其基本的进化方程如下: 其中,r1,r20,1为均匀分布的随机数;C1,C2均是正常数;t表示进化代数;Vt,Xt分别表示每个粒子的速度和位置;Pg,Pt分别是粒子群的全局最优和个体最优。为了改善基本PSO算法的收敛性能,Y?Shi等人提出了惯性权重的方法和用模糊控制器来动态自适应地改变惯性权重的技术。之后Jun Sun等人提出的具有函数形式的粒子群算法(QDP
2、SO)使粒子群算法的计算更加简单容易。最近一种新的QDPSO算法考虑了速度对位置的影响,通过速度的更新选择位置的更新方程。在经典粒子群算法的可调整参数中,惯性权重是非常重要的参数,较大的权重有利于提高算法的全局搜索能力,而较小的权重会增强算法的局部搜索能力。因此,对这种新的QDPSO算法的速度项引用惯性权重,通过研究4种方案,发现惯性权重的变化对具有量子行为的粒子群算法的收敛性具有很大改善。可以说惯性权重的适当设置对新的QDPSO算法性能也起着重要的作用。1 量子行为的粒子群优化算法及其改进1.1 QDPS0算法文献4的作者认为,若是在PSO系统下的个体粒子具有量子行为,则该粒子将会以与基本P
3、SO算法中的粒子不同的方式运动。在量子空间,粒子的速度和位置不能再依据“不确定原理”被同时确定,所以提出了QDPSO算法。该算法改变了基本PSO算法的粒子更新策略,只用了粒子的位置向量。QDPSO算法的粒子进化方程如下:其中,a,b,u0,1为均匀分布的随机数;pid是第i个粒子在第d维空间找到的局部最优解,pgd是群体在第d维空间找到的全局最优解;xid表示第i个粒子在第d维空间找到的当前值;而g必须满足条件:,才能保证算法的收敛。1.2 改进的粒子群算法新的QDPSO算法利用个体粒子的速度产生一个介于0,1之间的数来代替原算法中的由计算机随机产生的数,用以选择该粒子的位置更新方程。更新方程
4、和参数设定参考文献5。本文考虑到惯性权重随粒子的迭代次数变化影响个体粒子的速度引导该粒子向最优解靠拢,所以采用4种方案对该改进算法进行研究。通过使惯性权重随粒子的迭代次数变化,从而影响速度的更新方程:其中,采用4种惯性权重方案来影响速度的更新,然后与QDPSO算法进行性能比较:方案1 为从(1,0.875)递减的函数=1-k?0.125genmax。采用这种方案的QDPSO算法称为1-QDPSO;方案2 为从(0.9,0.4)递减的函数甜=0.9-k?0.5genmax。采用这种方案的QDPSO算法称为2-QDPSO;方案3 为一定值0.729 8,采用这种方案的QDPSO算法称为3-QDPS
5、O;方案4 为一凹函数(start-end)(ttmax)2+(-end)(2ttmax)+start,其中start=0.95,end=0.4,tmax为最大的迭代次数。采用这种方案的QDPSO算法称为4-QDPOS。综上所述,选择测试函数F1(x)和F2(x)分别为Sphere和Rastrigin(参数设置见文献4),改进后的算法流程如下:Step 1 初始化种群粒子的速度和位置;Step 2 通过对两个测试函数进行初始化计算,得到每个粒子的当前位置为粒子最佳位置Pbest,初始群体粒子位置的最优值为群体最佳位置gbest;Step 3 重新把粒子的位置代入测试函数进行计算,得到每个粒子新
6、的适应值,将其与Pbest比较,若较好,则将Pbest设置为新位置;并将其与gbest比较,若较好,则将gbest设置为新位置;Step 4 根据公式(6)更新粒子的速度;Step 5 用个体粒子的速度产生用以选择该粒子位置的更新方程的数据;Step 6 由Step 5产生的数据选择更新粒子位置的方程; Step 7 若未达到终止条件(足够小的适应值或预设的最大迭代次数),则返回Step 3。更新粒子速度时需要注意:如果粒子的速度超出预设的范围,则采取使粒子反向运动的策略,从而保证算法有效进行。1.3 算法的结果及数据分析目标函数为F1(x)和F2(x),基本参数是:c1=c2=2.05,g=
7、0.968 5,每种算法都在同一台计算机,同一环境下用Matlab 软件运行。结果如表1所示。表1的数值是对每个函数在粒子数为20个的条件下,测试50次,然后取平均得到的结果。从表中可以看出,对于函数F1(x),比较结果可以明显得知:在随粒子群维数增加的情况下,1-QDPSO是比QDPSO得到更好的解,其他几种改进方案的解都比较差;函数F2(x)在随粒子群维数增加的情况下,4种改进方案和QDPSO都能得出比较好的解。通过实验,可以看出:对于单峰函数F1(x),的递减不能太小,从方案1-QDPSO和2-QDPSO的结果就可以比较出来,而方案3-QDPSO和4-QDPSO的结果不好,可能是因为它们
8、搜索的区域太小,从而陷入局部最优解。 对于多峰函数F2(x),的变化对测试函数的解的精确度没有太大影响,说明了改进方案在此方面没有明显提高。接下来,我们还对算法的收敛速度进行了比较。结果如表2所示。表2是对函数测试50次后取得平均值的结果。可见对于函数F1(x),1-QDPSO和QDPSO都在10维的情况下收敛,而20维时只有1-QDPSO收敛,其他函数都没有收敛,导致这种结果的原因有2种:(1)各种方案随的变化,削弱或失去了调节能力,在达到最大迭代次数时也未收敛;(2)即使在算法已搜索到最优解附近时,由于局部搜索能力太差,跳过了最优解。对于函数F2(x),3-QDPSO,4-QDPSO,QD
9、PSO收敛速度都比较快,1=QDPSO和2-QDPSO的收敛速度就相对较慢一些。这是由于对多峰函数测试时,各种方案的初始化范围附近可能存在最优解,所以减少了迭代次数,加快了算法速度。通过对4种方案的研究,这里选取方案1应用于0-1背包问题,并得到理想的结果。2 对改进算法应用到0-1背包问题2.1 0-1背包问题的数学描述0-1背包问题是一种典型的组合优化问题。0-1背包问题的描述如下:假设有n个物品,其大小和价值分别为wi和ci(其中wi>0,ci>0,i=1,2,n),背包的容量假设为V(V>0)。要求在背包的容量限制内,使所装物品的总价值最大。该问题的数学模型可表示为:
10、其中,当将物品i装入背包时xi=1;否则xi=0。2.2 0-1背包问题的改进粒子群算法改进粒子群算法应用到0-1背包问题的思想:粒子群中粒子的个数与每个粒子的维数相等。先定义二进制数x,x只能取0和1。再把粒子的种群数看作背包的个数n,对于每个粒子xi(其中i=1,2,n表示粒子个数)有n个维数,即1个粒子有n个位置。然后初始化每个粒子的速度vij,(其中j=1,2,n表示每个粒子位置的维数),每个粒子的每一维都对应一个初始化了的速度。对公式(8)进行变化:解决背包问题的步骤:(1)初始化粒子的速度和位置;(2)将初始化的位置向量代人式(9),在所得每个粒子的解中找到最优解pbest,并令p
11、best=gbest;(3)通过式(6)更新粒子的速度,对所得最优解进行修正,然后再次代入函数方程中继续寻找新的最优解;(4)若达到终止条件,则结束迭代,输出到存储向量,即为所求结果;否则,k=k+1,转步骤(3)。2.3 实验仿真为了验证1-QDPSO求解01背包问题的可行性及有效性,这里进行了2组实验,每组实验用1-QDPSO算法进行测试,每组算法运行50次。实验一:取参数popsize=10,dimsize=10,c1=c2=2.05,genmax=1 000,g=0.968 5;N=10,V=269,W=95,4,60,32,23,72,80,62,65,46,C=55,10,47,5
12、,4,50,8,61,85,87),得到实验结果是max f=295,收敛平均迭代次数11。实验二:取参数popsize=20,dimsize=20,c1=c2=2.05,genmax=1 000,g=0.968 5;N=20,V=878,W=92,4,43,83,84,68,92,82,6,44,32,18,56,83,25,96,70,48,14,58,C=44,46,90,72,91,40,75,35,8,54,78,40,77,15,61,17,75,29,75,63,得到实验结果是max f=1024,收敛平均迭代次数23。1-QDPSO算法求解0-1背包问题,与文献9中提出的用带有死亡罚函数的粒子群优化算法求解0-1背包问题相比,其运行速度明显提高。
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 部编版一年级上册第四单元《秋天》教案
- 肉制品加工题目及答案
- 2023-2024学年四川省凉山州安宁联盟高二下学期期末联考数学试题(解析版)
- 2023-2024学年山东省淄博市高二下学期期末考试数学试题(解析版)
- 2023-2024学年广东省东莞市高二下学期7月期末教学质量检查数学试题(解析版)
- 浩如宇装饰新材料有限公司年产5000立方米科技木技术改造项目环评资料环境影响
- 2025年秋三年级上册语文同步教案 22 读不完的大书
- 小学地震演练活动方案
- 给爱人的检讨书
- 人造板产品质量省监督抽查实施细则
- 《职场沟通技巧》课件
- 辽宁大学《生态学》2023-2024学年第一学期期末试卷
- 《Maya 2022三维动画制作标准教程(全彩版)》第08章 渲染与输出
- 2024秋国开《人工智能专题》专题测验1-3答案
- 《我国上市公司财务舞弊探析案例-长园集团为例》14000字(论文)
- DB12-T 3035-2023 建筑消防设施维护保养技术规范
- DB35T 2067-2022 锅炉用固体废弃物燃料性能评价规则
- 2024年全国统一高考英语试卷(新课标Ⅰ卷)含答案
- 2023-2024学年广东省广州市花都区八年级(下)期末数学试卷(含答案)
- 山东省青岛市莱西市2023-2024学年高一下学期7月期末英语试题
- 西藏自治区劳动合同
评论
0/150
提交评论