细菌觅食算法MATLAB实现
更新时间:2023-03-17 00:31:01 阅读量: 教育文库 文档下载
- 细菌觅食算法路径优化代码推荐度:
- 相关推荐
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%*********************细菌觅食算法**********************
%%%%%%%%%%%%%%%%%%%-----BFA算法-----%%%%%%%%%%%%%%%%%%% clear; clc;
%-----(1)初始化参数-----
bounds = [-5.12 5.12;-5.12 5.12]; % 函数变量范围 p = 2; % 搜索范围的维度 s = 26; % 细菌的个数 Nc = 50; % 趋化的次数
Ns = 4; % 趋化操作中单向运动的最大步数
C(:,1) = 0.001*ones(s,1); % 翻转选定方向后,单个细菌前进的步长 Nre = 4; % 复制操作步骤数 Ned = 2; % 驱散(迁移)操作数 Sr = s/2; % 每代复制(分裂)数 Ped = 0.25; % 细菌驱散(迁移)概率
d_attract = 0.05; % 吸引剂的数量
ommiga_attract = 0.05; % 吸引剂的释放速度 h_repellant = 0.05; % 排斥剂的数量
ommiga_repellant = 0.05;% 排斥剂的释放速度 for i = 1:s % 产生初始细菌个体的位置 P(1,i,1,1,1) = -5.12 + rand*10.24; P(2,i,1,1,1) = -5.12 + rand*10.24; end
%------------------细菌趋药性算法循环开始--------------------- %-----(2)驱散(迁移)操作开始----- for l = 1:Ned
%-----(3)复制操作开始----- for k = 1:Nre
%-----(4)趋化操作(翻转或游动)开始----- for j = 1:Nc
%-----(4.1)对每一个细菌分别进行以下操作----- for i = 1:s
%-----(4.2)计算函数J(i,j,k,l),表示第i个细菌在第l次驱散第k次 %----------复制第j次趋化时的适应度值----- J(i,j,k,l) = Cost(P(:,i,j,k,l));
%-----(4.3)修改函数,加上其它细菌对其的影响----- Jcc = sum(-d_attract*exp(-ommiga_attract*((P(1,i,j,k,l)-...
P(1,1:26,j,k,l)).^2+(P(2,i,j,k,l)-P(2,1:26,j,k,l)).^2))) +... sum(h_repellant*exp(-ommiga_repellant*((P(1,i,j,k,l)-... P(1,1:26,j,k,l)).^2+(P(2,i,j,k,l)-P(2,1:26,j,k,l)).^2))); J(i,j,k,l) = J(i,j,k,l) + Jcc;
%-----(4.4)保存细菌目前的适应度值,直到找到更好的适应度值取代之
-----
Jlast = J(i,j,k,l);
%-----(4.5)翻转,产生一个随机向量C(i),代表翻转后细菌的方向----- Delta(:,i) = (2*round(rand(p,1))-1).*rand(p,1);
% PHI表示翻转后选择的一个随机方向上前进 PHI = Delta(:,i)/sqrt(Delta(:,i)'*Delta(:,i));
%-----(4.6)移动,向着翻转后细菌的方向移动一个步长,并且改变细菌的位置-----
P(:,i,j+1,k,l) = P(:,i,j,k,l) + C(i,k)*PHI;
%-----(4.7)计算细菌当前位置的适应度值----- J(i,j+1,k,l) = Cost(P(:,i,j+1,k,l)); %-----(4.8)游动-----
m = 0; % 给游动长度计数器赋初始值
while(m < Ns) % 未达到游动的最大长度,则循环 m = m + 1;
% 新位置的适应度值是否更好?如果更好,将新位置的适应度值 % 存储为细菌i目前最好的适应度值 if(J(i,j+1,k,l) Jlast = J(i,j+1,k,l); %保存更好的适应度值 % 在该随机方向上继续游动步长单位,修改细菌位置 P(:,i,j+1,k,l) = P(:,i,j+1,k,l) + C(i,k)*PHI; % 重新计算新位置上的适应度值 J(i,j+1,k,l) = Cost(P(:,i,j+1,k,l)); else % 否则,结束此次游动 m = Ns; end end J(i,j,k,l) = Jlast; % 更新趋化操作后的适应度值 end % 如果i %-----(5)如果j x = P(1,:,j,k,l); y = P(2,:,j,k,l); clf plot(x,y,'h') % h表示以六角星绘图 axis([-5 5 -5 5]); % 设置图的坐标图 pause(.1) % 暂停0.1秒后继续 end %----------------下面进行复制操作---------------- %-----(6)复制----- %-----(6.1)根据所给的k和l的值,将每个细菌的适应度值按升序排序----- Jhealth = sum(J(:,:,k,l),2); % 给每个细菌设置健康函数值 [Jhealth,sortind] = sort(Jhealth); % 按健康函数值升序排列函数 P(:,:,1,k+1,l) = P(:,sortind,Nc+1,k,l); C(:,k+1) = C(sortind,k); %-----(6.2)将代价小的一半细菌分裂成两个,代价大的一半细菌死亡----- for i = 1:Sr % 健康值较差的Sr个细菌死去,Sr个细菌分裂成两个子细菌,保持个体总数的s一致性 P(:,i+Sr,1,k+1,l) = P(:,i,1,k+1,l); C(i+Sr,k+1) = C(i,k+1); end %-----(7)如果k %-----(8)趋散,对于每个细菌都以Ped的概率进行驱散,但是驱散的细菌群体的总数 %--------保持不变,一个细菌被驱散后,将被随机重新放置到一个新的位置----- for m = 1:s % 产生随机数,如果既定概率大于该随机数,细菌i灭亡,随机产生新的细菌i if(Ped> rand) P(1,m,1,1,1) = -5.12 + rand*10.24; P(2,m,1,1,1) = -5.12 + rand*10.24; else P(:,m,1,1,l+1) = P(:,m,1,Nre+1,l); % 未驱散的细菌 end end end % 如果l %-------------------------报告---------------------- reproduction = J(:,1:Nc,Nre,Ned); % 每个细菌最小的适应度值 [Jlastreproduction,O] = min(reproduction,[],2); [BestY,I] = min(Jlastreproduction) Pbest = P(:,I,O(I,:),k,l) % 求解Shaffer's函数的最小值 % Shaffer's函数表示如下: % f(x)=0.5+(sin(sqrt(x1^2+x2^2))^2-0.5)/(1.0+0.001(x1^2+x2^2))^2 function cost = Cost(x) cost = 0.5 + (sin(sqrt(x(1)^2+x(2)^2))^2-0.5)/(1.0+0.001*(x(1)^2+x(2)^2))^2; % function fposition=Cost(x) % % x=fix(100*rand(2,1)); % p=0;q=0; % for k=1:5 % p=p+k*cos((k+1)*x(1)+k); % q=q+k*cos((k+1)*x(2)+k); % end % % fposition=p*q+(x(1)+1.42513)^2+(x(2)+.80032)^2;
正在阅读:
细菌觅食算法MATLAB实现03-17
2011年揭阳市中考历史真题07-21
基于专家控制器的水厂加氯串级控制系统08-11
分析化学习题参考解答06-10
浅谈我国企业人力资源管理现状及发展趋势06-14
外研版八年级上英语期中试卷05-05
2018-2024年中国浙江省房地产市场运营状况分析及投资前景价值评估报告(目录) - 图文11-27
ISO9001年度管理评审总结报告08-27
新监管形势下罪犯心理矫治的创新问题研究04-20
我学会了钓鱼作文600字07-14
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 觅食
- 细菌
- 算法
- 实现
- MATLAB
- 郭锡良版古代汉语习题解答
- 大型地下车库设计实例要点分析
- 数据结构
- 广东商学院操作系统试题(10-11)A
- 流行病学作业题
- 2017-2022年餐饮设备市场发展前景分析及供需格局研究预测报告目录
- 关于建立健全创建全国文明城市长效机制的方案
- 幼儿园规章制度汇编
- 搭配中的学问公开课(小丑、详细)(教案、PPT、反思、学习单)
- 2019年导游资格《政策法规》试题及答案(卷五)
- 《审计风云》观后感
- 环流(有环流控制的可逆晶闸管-电动机系统)
- 首届中国“互联网+”大学生创新创业大赛筹备工作方案筹备方案
- ZJ管理集团公司岗位职责说明书 doc
- 2011年广东省中考语文试卷
- 现代服务业行业研究报告
- 海宁市卫生局机关工作制度汇编
- Linux操作系统实验报告2
- 《计算机组成与工作原理》试卷(附答案)
- 台州职业技术学院顶岗实习管理办法(试行)