基于Matlab实现嗅觉优化算法的示例代码
目录
- 1.概述
- 2.37 个 CEC 基准测试函数代码
- 3.F1 Matlab代码仿真
1.概述
嗅觉剂优化是一种新颖的优化算法，旨在模仿气味分子源尾随的药剂的智能行为。该概念分为三个阶段（嗅探，尾随和随机）是独特且易于实现的。此上传包含 SAO 在 37 个 CEC 基准测试函数上的实现。
2.37 个 CEC 基准测试函数代码
function [lb,ub,dim,fobj] = Select_Function(F)
switch F
case 'F1'
%Admijan
fobj = @F1;
lb=[-1 -1];
ub=[2 1];
dim=2;
case 'F2'
%Beale
fobj = @F2;
dim=2;
lb=-4.5*ones(1,dim);
ub=4.5*ones(1,dim);
case 'F3'
%Bird
fobj = @F3;
dim=2;
lb=-2*pi*ones(1,dim);
ub=2*pi*ones(1,dim);
case 'F4'
%Bohachevsky
fobj = @F4;
dim=2;
lb=-100*ones(1,dim);
ub=100*ones(1,dim);
case 'F5'
% Booth
fobj = @F5;
dim=2;
lb=-10*ones(1,dim);
ub=10*ones(1,dim);
case 'F6'
%Branin RCOS1
fobj = @F6;
lb=[-5,0];
ub=[10, 15];
dim=2;
case 'F7'
%Branin RCOS2
fobj = @F7;
dim=2;
lb=-5*ones(1,dim);
ub=15*ones(1,dim);
case 'F8'
%Brent
fobj = @F8;
dim=2;
lb=-10*ones(1,dim);
ub=10*ones(1,dim);
case 'F9'
%Bukin F2
fobj = @F9;
dim=2;
lb=[-15 -3];
ub=[-5 3];
case 'F10'
%six-hump
fobj = @F10;
dim=2;
lb=-5*ones(1,dim);
ub=5*ones(1,dim);
case 'F11'
%Chichinadze
fobj = @F11;
dim=2;
lb=-30*ones(1,dim);
ub=30*ones(1,dim);
case 'F12'
%Deckkers-Aarts
fobj = @F12;
dim =2;
lb=-20*ones(1,dim);
ub=20*ones(1,dim);
case 'F13'
%Easom
dim=2;
fobj=@F13;
lb=-10*ones(1,dim);
ub=10*ones(1,dim);
case 'F14'
%Matyas
fobj = @F14;
dim=2;
lb=-10*ones(1,dim);
ub=10*ones(1,dim);
case 'F15'
%McComick
fobj = @F15;
dim=2;
lb=-10*ones(1,dim);
ub=10*ones(1,dim);
case 'F16'
%Michalewicz2
fobj = @F16;
dim=2;
lb=0*ones(1,dim);
ub=pi*ones(1,dim);
case 'F17'
%Quadratic
fobj = @F17;
dim=2;
lb=-10*ones(1,dim);
ub=10*ones(1,dim);
case 'F18'
%Schaffer
dim=2;
fobj = @F18;
lb=-100*ones(1,dim);
ub=100*ones(1,dim);
case 'F19'
%StyblinskiTang
fobj = @F19;
dim=2;
lb=-5*ones(1,dim);
ub=5*ones(1,dim);
case 'F20'
%Box-Betts
fobj = @F20;
dim=3;
lb=[0.9 9 0.9];
ub=[1.2 11.2 1.2];
case 'F21'
%Colville
fobj = @F21;
dim=4;
lb=-1*ones(1,dim);
ub=1*ones(1,dim);
case 'F22'
%Csendes
fobj = @F22;
dim=4;
lb=-1*ones(1,dim);
ub=1*ones(1,dim);
case 'F23'
% Michalewicz 5
fobj = @F23;
dim=5;
lb=0*ones(1,dim);
ub=pi*ones(1,dim);
case 'F24'
%Miele Cantrell
dim=4;
fobj = @F24;
lb=-1*ones(1,dim);
ub=1*ones(1,dim);
case 'F25'
% Step
fobj = @F25;
dim=5;
lb=-100*ones(1,dim);
ub=100*ones(1,dim);
case 'F26'
%Michalewicz
fobj = @F26;
dim=10;
lb=0*ones(1,dim);
ub=pi*ones(1,dim);
case 'F27'
%Shubert
fobj = @F27;
dim=5;
lb=-10*ones(1,dim);
ub=10*ones(1,dim);
case 'F28'
%Ackley
dim=30;
fobj = @F28;
lb=-32*ones(1,dim);
ub=32*ones(1,dim);
case 'F29'
%Brown
fobj = @F29;
dim=30;
lb=-1*ones(1,dim);
ub=4*ones(1,dim);
case 'F30'
%Ellipsoid
dim=2;
fobj = @F30;
lb=-5.12*ones(1,dim);
ub=5.12*ones(1,dim);
case 'F31'
% Grienwank
fobj = @F31;
dim=30;
lb=-100*ones(1,dim);
ub=100*ones(1,dim);
case 'F32'
%Mishra
fobj = @F32;
dim=30;
lb=0*ones(1,dim);
ub=1*ones(1,dim);
case 'F33'
%Quartic
dim=30;
fobj = @F33;
lb=-1.28*ones(1,dim);
ub=1.28*ones(1,dim);
case 'F34'
%Rastrigin
fobj = @F34;
dim=30;
lb=-5.12*ones(1,dim);
ub=5.12*ones(1,dim);
case 'F35'
%Rosenbrock
fobj = @F35;
dim=30;
lb=-30*ones(1,dim);
ub=30*ones(1,dim);
case 'F36'
% Salomon
fobj = @F36;
dim=30;
lb=-100*ones(1,dim);
ub=100*ones(1,dim);
case 'F37'
%Sphere
fobj = @F37;
dim=30;
lb=-100*ones(1,dim);
ub=100*ones(1,dim);
end
end
function o=F1(x)
% Adjiman
o=(cos(x(:,1)).*sin(x(:,2))-x(:,1)./(x(:,2).^2+1));
end
function o=F2(x)
% Beale
o=(1.5-x(:,1)+(x(:,1).*(x(:,2)))).^2+(2.25-x(:,1)+(x(:,1).*(x(:,2)).^2)).^2+...
(2.625-x(:,1)+(x(:,1).*(x(:,2)).^3)).^2;
end
function o=F3(x)
% Bird
o=sin(x(:,2)).*(exp(1-cos(x(:,1))).^2)+cos(x(:,1)).*(exp(1-sin(x(:,2))).^2)...
+(x(:,1)+(x(:,2))).^2;
end
function o=F4(x)
% Bohachevsky
W=0;
[a,dim]=size(x);
for i=1:dim-1
W=W+x(:,i).^2+2.*x(:,i+1).^2-0.3.*cos(3.*pi.*x(:,i+1))-0.4.*cos(4.*pi.*(x(:,i+1)))+0.7;
end
o=W;
end
function o=F5(x)
%Booth
o=(x(:,2)-(5.1*x(:,1).^2/(4*pi*2))+(5*x(:,1)/pi)-6).^2+...
10*(1-1/(8*pi)).*cos(x(:,1))+10;
end
function o=F6(x)
% Branin RCOS 1
o=(x(:,2)-(5.1*x(:,1).^2/(4*pi*2))+(5*x(:,1)/pi)-6).^2+...
10*(1-1/(8*pi)).*cos(x(:,1))+10;
end
function o=F7(x)
% Branin RCOS 2
a=1; b=5.1/(4*pi^2); c=5/pi; d=6; e=10; g=1/(8*pi);
f1=a*(x(:,2)-b*x(:,1).^2+c*x(:,1)-d).^2;
f2=e*(1-g)*cos(x(:,1)).*cos(x(:,2));
f3=log(x(:,1).^2+x(:,2)+1);
o=-1/(f1+f2+f3+e);
end
function o=F8(x)
%Brent
o=(x(:,1)+10).^2+(x(:,1)+10).^2+exp(-x(:,1).^2-x(:,2).^2);
end
function o=F9(x)
%Bukin F2
o=(abs(x(:,1)-0.01.*x(:,2).^2))+0.01.*abs(x(:,2)+10);
end
function o=F10(x)
%Camel Six Hump
o=(4-2.1*x(:,1).^2+(x(:,1).^4)/3).*x(:,1).^2+x(:,1).*x(:,2)+...
(4*x(:,2).^2-4).*x(:,2).^2;
end
function o=F11(x)
%Chichinadze
o=x(:,1).^2-12*x(:,1)+11+10*cos(pi*x(:,1)/2)+8*sin(5*pi*x(:,1)/2)-...
((1/5)^0.5)*exp(-0.5*(x(:,2)-0.5).^2);
end
function o=F12(x)
% Deckkers-Aarts
o=10^5*x(:,1).^2+x(:,2).^2-(x(:,1).^2+x(:,2).^2).^2+...
10^(-5)*(x(:,1).^2+x(:,2).^2).^4;
end
function o = F13(x)
% Easom
o=-cos(x(:,1)).*cos(x(:,2)).*exp(-(x(:,1)-pi).^2-(x(:,2)-pi).^2);
end
function o=F14(x)
% Evaluate Matyas
o=0.26*(x(:,1).^2+x(:,2).^2)-0.48*x(:,1).*x(:,2);
end
function o=F15(x)
% McCormick
o=mccormick(x);%
end
function o=F16(x)
% Michalewicz2
[~,d]=size(x);
W=0;
for i=1:d
W=sin(x(:,1)).*sin(i*x(:,i).^2/pi).^2*d;
end
o=-W;
end
function o=F17(x)
% Quadratic
o=-3803.84-138.08*x(:,1)-232.92*x(:,2)+128.08*x(:,1).^2+203.64*x(:,2).^2+182.25*x(:,1).*x(:,2);
end
function o=F18(x)
% Evaluate Schaffer
[~,d]=size(x);
w=0;
for i=1:d-1
w=w+((x(i).^2+x(i+1).^2).^.5).*(sin(50.*(x(i).^2+x(i+1).^2).^0.1)).^2;
end
o=w;
end
function o=F19(x)
% Styblinki's Tang
[~,d]=size(x);
W=0;
for i=1:d
W=W+(x(:,i).^4-16.*x(:,i).^2+5.*x(:,i));
end
o=W.*0.5;
end
function o=F20(x)
% Box-Betts
[~,d]=size(x);
W=0;
for i=1:d
g=exp(-0.1.*(i+1)).*x(:,1)-exp(-0.1.*(i+1)).*x(:,2)-((exp(-0.1.*(i+1)))-exp(-(i+1)).*x(:,3));
W=W+g.^2;
end
o=W;
end
function o=F21(x)
% Colville
o=100*(x(:,1)-x(:,2).^2).^2+(1-x(:,1)).^2+90*(x(:,4)-x(:,3).^2).^2+...
(1-x(:,3)).^2+10.1*((x(:,2)-1).^2+(x(:,4)-1).^2)+...
19.8*(x(:,2)-1).*(x(:,4)-1);
end
function o=F22(x)
% Csendes
[~,d]=size(x);
aa=0;
for i=1:d
aa=aa+x(:,i).^6.*(2+sin(1/x(:,i)));
end
o=aa;
end
function o=F23(x)
% Michalewicz 5
[~,d]=size(x);
W=0;
for i=1:d
W=sin(x(:,1)).*sin(i*x(:,i).^2/pi).^2*d;
end
o=-W;
end
function o=F24(x)
%Miele Cantrell
o=(exp(-x(:,1))-x(:,2)).^4+100*(x(:,2)-x(:,3)).^6+...
(tan(x(:,3)-x(:,4))).^4+x(:,1).^8;
end
function o=F25(x)
% Evaluate Step
[~,d]=size(x);
W=0;
for i=1:d
W=W+(floor(x(:,i)+0.5)).^2;
end
o=W;
end
function o=F26(x)
% Evaluate Michalewicz 10
[~,d]=size(x);
W=0;
for i=1:d
W=sin(x(:,1)).*sin(i*x(:,i).^2/pi).^2*d;
end
o=-W;
end
function o=F27(x)
% shubert
[~,d]=size(x);
s1=0;
s2=0;
for i = 1:d
s1 = s1+i*cos((i+1)*x(1)+i);
s2 = s2+i*cos((i+1)*x(2)+i);
end
o = s1*s2;
end
% F28--Ackley
function o = F28(x)
dim=size(x,2);
o=-20*exp(-.2*sqrt(sum(x.^2)/dim))-exp(sum(cos(2*pi.*x))/dim)+20+exp(1);
end
function o=F29(x)
[~,d]=size(x);
% Brown
a=0;
for i=1:d-1
a=(x(:,i).^2).^(x(:,i+1)+1)+(x(:,i+1).^2).^(x(:,i).^2+1);
end
o=a;
end
function o=F30(x)
% Ellipsoid
[~,d]=size(x);
W=0;
for i=1:d
W=W+i.*x(:,1).^2;
end
o=W;
end
%Grienwank
function o=F31(x)
o=griewank(x);
end
function o=F32(x)
% Mishra
[~,d]=size(x);
a=0;
for i=1:d-1
a=a+x(:,i);
end
aa=d-a;
b=0;
for j=1:d-1
b=b+x(:,j);
end
W=abs((1+d-b).^aa);
o=W;
end
% --Quartic
function o = F33(x)
dim=size(x,2);
o=sum([1:dim].*(x.^4))+rand;
end
%Rastrigin
function o=F34(x)
o=rastrigin(x);
end
% Rosenbrock
function o = F35(x)
dim=size(x,2);
o=sum(100*(x(2:dim)-(x(1:dim-1).^2)).^2+(x(1:dim-1)-1).^2);
end
function o=F36(x)
% salomon
x2 = x.^2;
sumx2 = sum(x2, 2);
sqrtsx2 = sqrt(sumx2);
o = 1 - cos(2 .* pi .* sqrtsx2) + (0.1 * sqrtsx2);
end
function o = F37(x)
%Sphere
o=sum(x.^2);
end
function o=Ufun(x,a,k,m)
o=k.*((x-a).^m).*(x>a)+k.*((-x-a).^m).*(x<(-a));
end
3.F1 Matlab代码仿真
到此这篇关于基于Matlab实现嗅觉优化算法的示例代码的文章就介绍到这了,更多相关Matlab嗅觉优化算法内容请搜索我们以前的文章或继续浏览下面的相关文章希望大家以后多多支持我们！
相关推荐
-
基于Matlab实现野狗优化算法的示例代码
目录 1.概述 2.捕食过程的数学模型 2.1 种群初始化 2.2 群体攻击过程 2.3 迫害攻击过程 2.4 野狗的存活率 3.Matlab代码实现 3.1 代码 3.2 结果 1.概述 野狗优化算法(Dingo Optimization Algorithm, DOA)模仿澳大利亚野狗的社交行为.DOA算法的灵感来源于野狗的狩猎策略,即迫害攻击.分组策略和食腐行为.为了提高该方法的整体效率和性能,在DOA中制定了三种与四条规则相关联的搜索策略,这些策略和规则在搜索空间的强化(开发)和多样化(探
-
Python&Matlab实现蚂蚁群算法求解最短路径问题的示例
目录 1知识点 1.1 蚁群算法步骤 1.2 蚁群算法程序 2蚂蚁算法求解最短路径问题——Python实现 2.1源码实现 2.2 ACA_TSP实现 3 蚂蚁算法求解最短路径问题——Matlab实现 3.1流程图 3.2代码实现 3.3结果 1 知识点 详细知识点见:智能优化算法—蚁群算法(Python实现) 我们这一节知识点只讲蚁群算法求解最短路径步骤及流程. 1.1 蚁群算法步骤 设蚂蚁的数量为m,地点的数量为n,地点i与地点j之间相距Dij,t时刻地点i与地点j连接的路径上的信息素浓度为
-
利用Matlab实现迭代适应点算法
目录 1.算法描述 2.工具函数 3.函数调用 4.优势与不足 道格拉斯-普克算法(Douglas–Peucker algorithm,亦称为拉默-道格拉斯-普克算法.迭代适应点算法.分裂与合并算法)是将曲线近似表示为一系列点,并减少点的数量的一种算法.它的优点是具有平移和旋转不变性,给定曲线与阈值后,抽样结果一定. 1.算法描述 1.在曲线首尾两点间虚连一条直线,求出其余各点到该直线的距离. 2.选其最大者与阈值相比较,若大于阈值,则离该直线距离最大的点保留,否则将直线两端点间各点全部舍去.
-
Python&Matlab实现灰狼优化算法的示例代码
目录 1 灰狼优化算法基本思想 2 灰狼捕食猎物过程 2.1 社会等级分层 2.2 包围猎物 2.3 狩猎 2.4 攻击猎物 2.5 寻找猎物 3 实现步骤及程序框图 3.1 步骤 3.2 程序框图 4 Python代码实现 5 Matlab实现 1 灰狼优化算法基本思想 灰狼优化算法是一种群智能优化算法,它的独特之处在于一小部分拥有绝对话语权的灰狼带领一群灰狼向猎物前进.在了解灰狼优化算法的特点之前,我们有必要了解灰狼群中的等级制度. 灰狼群一般分为4个等级:处于第一等级的灰狼用α表示,处于第
-
Matlab实现遗传算法的示例详解
目录 1算法讲解 1.1何为遗传算法 1.2遗传算法流程描述 1.3关于为什么要用二进制码表示个体信息 1.4目标函数值与适应值区别 1.5关于如何将二进制码转化为变量数值 1.6关于代码改进 2MATLAB自带ga函数 2.1问题描述 2.2自带函数使用 3自编遗传算法各部分代码及使用 3.1代码使用 3.2Genetic1--主函数 3.3PI(PopulationInitialize)--产生初始种群 3.4Fitness--计算目标函数值 3.5FitnessF--计算适应值 3.6Tr
-
基于Matlab实现鲸鱼优化算法的示例代码
目录 1.鲸鱼优化算法建模 1.1 包围猎物 1.2 螺旋狩猎 1.3 搜索猎物 1.4 算法流程图 2.Matlab代码实现 2.1 结果 2.2 代码 1.鲸鱼优化算法建模 鲸鱼优化算法(WOA)是澳大利亚学者Mirjaili等于2016年提出的群体智能优化算法,根据座头鲸的捕猎行为实现优化搜索的目的.其中,每个鲸鱼可以看作一个粒子,每个粒子作为不同的决策变量.WOA的实现过程主要包括包围猎物.螺旋狩猎和随机搜索3个阶段,其数学模型如下: 1.1 包围猎物 1.2 螺旋狩猎 1.3 搜索猎物
-
Python和Matlab实现蝙蝠算法的示例代码
目录 1前言 2 蝙蝠算法原理细讲 3 详细步骤 4Python实现 4.1代码 4.2结果 5Matlab实现 5.1 代码 5.2 结果 5.3 展望 1 前言 蝙蝠算法是2010年杨教授基于群体智能提出的启发式搜索算法,是一种搜索全局最优解的有效方法.该算法基于迭代优化,初始化为一组随机解,然后迭代搜寻最优解,且在最优解周围通过随机飞行产生局部新解,加强局部搜索速度.该算法具有实现简单.参数少等特点. 该算法主要用于目标函数寻优,基于蝙蝠种群利用产生的声波搜索猎物和控制飞行方向的特征来实现
-
基于Matlab实现嗅觉优化算法的示例代码
目录 1.概述 2.37 个 CEC 基准测试函数代码 3.F1 Matlab代码仿真 1.概述 嗅觉剂优化是一种新颖的优化算法,旨在模仿气味分子源尾随的药剂的智能行为.该概念分为三个阶段(嗅探,尾随和随机)是独特且易于实现的.此上传包含 SAO 在 37 个 CEC 基准测试函数上的实现. 2.37 个 CEC 基准测试函数代码 function [lb,ub,dim,fobj] = Select_Function(F) switch F case 'F1' %Admijan fobj = @
-
Matlab实现黑洞优化算法的示例代码
目录 前言 1.概述 1.1黑洞算法 1.2黑洞搜索优化算法 1.3黑洞搜索算法的实现过程 2.Matlab代码实现 2.1主函数 2.2目标函数 2.3黑洞优化算法 3.结果展现 前言 应用的领域很多. 1.概述 1.1黑洞算法 根据黑洞现象原理首次提出BH 算法,它在传统PSO基础上引入了新的机制,有效地提高了收敛速度并防止了陷入局部极值的情况发生;但是该方法却没有提及如何确定黑洞边界和如何处理吸收星体的问题. Hatamlou BH算法进行了完善,让其更加接近于黑洞的自然现状,使其具有黑洞
-
基于Matlab实现绘制3D足球的示例代码
目录 绘制讲解 数据来源及说明 硬算顶点连接情况 三角剖分 正交变换 充气 完整代码 世界杯教你用MATLAB画个超逼真的足球, 需要准备Partial Differential Equation Toolbox工具箱,同时因为用到了polyshape类所以至少需要R2017b版本. 绘制讲解 数据来源及说明 我是真的不想写注释了太麻烦了,给大家讲一下我的思路希望能够看懂,首先足球的数据点是通过: [B,XYZ]=bucky; 导入的,但是导入的只有边链接信息,并没有给出哪几个点构成正五边形哪几
-
基于Matlab实现多目标粘液霉菌算法的示例代码
目录 1.概述 2.算法过程 3.数值实验 4.Matlab代码实现 1.概述 多目标粘液霉菌算法(MOSMA),这是最近开发的粘液霉菌算法(SMA)的多目标变体,用于处理工业中的多目标优化问题.最近,为了处理优化问题,已经为优化社区推荐了几种元启发式和进化优化技术.在评估多目标优化(MOO)问题时,这些方法往往受到低质量解的影响,而不是解决识别帕累托最优解的准确估计和增加所有目标的分布的目标函数.SMA方法遵循从实验室实验中粘液霉菌振荡行为中获得的逻辑.与其他成熟的方法相比,SMA算法显示出强
-
Python实现七大查找算法的示例代码
查找算法 -- 简介 查找(Searching)就是根据给定的某个值,在查找表中确定一个其关键字等于给定值的数据元素. 查找表(Search Table):由同一类型的数据元素构成的集合 关键字(Key):数据元素中某个数据项的值,又称为键值 主键(Primary Key):可唯一的标识某个数据元素或记录的关键字 查找表按照操作方式可分为: 1.静态查找表(Static Search Table):只做查找操作的查找表.它的主要操作是: ①
-
Java实现查找算法的示例代码(二分查找、插值查找、斐波那契查找)
目录 1.查找概述 2.顺序查找 3.二分查找 3.1 二分查找概述 3.2 二分查找实现 4.插值查找 4.1 插值查找概述 4.2 插值查找实现 5.斐波那契查找 5.1 斐波那契查找概述 5.2 斐波那契查找实现 5.3 总结 1.查找概述 查找表: 所有需要被查的数据所在的集合,我们给它一个统称叫查找表.查找表(Search Table)是由同一类型的数据元素(或记录)构成的集合. 查找(Searching): 根据给定的某个值,在查找表中确定一个其关键字等于给定值的数据元素(或记录).