| 查看: 2246 | 回复: 1 | ||
[求助]
粒子群算法每次计算结果相差很大
|
|
我最近在用基于杂交算法的粒子群算法做计算，共四个未知数，但求出的结果相差很多，我只有一次求出了正确的结果具体参数如下：粒子数为40，学习因子1和2都取为2，惯性权重取为0.7，杂交概率取为0.9，杂交池的大小比例取为0.2，最大迭代次数为10000，自变量个数为4,， 适应度计算如下：function F=fitness(x)（求最小值） ea=1+((x(1)^2)/2)*((x(2)^2-x(3)^2))/((x(2)^2-x(3)^2)^2+x(4)^2*x(3)^2); eb=(((x(1)^2)/2)*x(4)*x(3))/((x(2)^2-x(3)^2)^2+(x(4)^2)*(x(3)^2)); n=((ea+(ea^2+eb^2)^0.5)/2)^0.5; k=((-ea+(ea^2+eb^2)^0.5)/2)^0.5; L1a=500; T2a=0.917; T5a=0.909 nsa=1.46; R12=((nsa-n)^2+k^2)/((nsa+n)^2+k^2); R01=((1-nsa)^2)/((1+nsa)^2); A22=exp((-4*3.14*k*0.2*10000000)/L1a); A25=exp((-4*3.14*k*0.5*10000000)/L1a); t2=((1-R01)^2*(1-R12)^2*A22)/((1-R01*R12)^2*(1-A22^2*R12^2)-A22^2*R01*((1-R12)^4*R01-R12*(1-R12)^2*(1-R01*R12)-(1-R12)^2*R12*(1-R12*R01))); t5=((1-R01)^2*(1-R12)^2*A25)/((1-R01*R12)^2*(1-A25^2*R12^2)-A25^2*R01*((1-R12)^4*R01-R12*(1-R12)^2*(1-R01*R12)-(1-R12)^2*R12*(1-R12*R01))); F=((t2-T2a)/T2a)^2+((t5-T5a)/T5a-0.000577)^2; end |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有128人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab论文分享（包括神经网络，遗传算法，粒子群算法）
已经有442人回复
- 求粒子群优化算法进行变量优选的matlab程序
已经有20人回复
- 【PSO】粒子群算法，大家接触多少？ Particle Swarm Optimization for everyone
已经有19人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 【求助】多目标pso算法的MATLAB程序
已经有3人回复
ajitai
铁杆木虫 (职业作家)
- 应助: 8 (幼儿园)
- 金币: 8817.4
- 散金: 5370
- 红花: 11
- 帖子: 3525
- 在线: 418.1小时
- 虫号: 854389
- 注册: 2009-09-22
- 性别: GG
- 专业: 金属材料表面科学与工程
2楼2013-05-28 18:19:47
50