狐狸优化算法（Matlab代码实现）

👨‍🎓个人主页：研学社的博客

💥💥💞💞欢迎来到本博客❤️❤️💥💥


🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。


⛳座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁
目录
💥1 概述
📚2 运行结果
🎉3 参考文献
🌈4 Matlab代码及文章讲解

💥1 概述

本文提出了一种新的受自然启发的优化算法，称为狐狸优化器（Fox），该算法模拟了狐狸在自然界中捕食猎物时的觅食行为。该算法基于测量狐狸与其猎物之间的距离以执行有效跳跃的技术。在给出FOX的数学模型和算法之后，使用五个经典基准函数和CEC2019基准测试函数来评估其性能。FOX算法还与蜻蜓优化算法（DA）、粒子群优化（PSO）、适应度相关优化（FDO）、灰狼优化（GWO）、鲸鱼优化算法（WOA）、黑猩猩优化算法（ChOA）、蝴蝶优化算法（BOA）和遗传算法（GA）进行了比较

结果表明，FOX优于上述算法。随后，使用Wilcoxon秩和检验来确保FOX在统计显著性方面优于比较算法。此外，还进行了参数敏感性分析，以显示FOX中不同的探索和开发行为。本文还使用FOX来解决工程问题，例如压力容器设计，并且还用于解决发电：经济负荷分配问题。FOX在针对GWO、PSO、WOA和FDO优化问题方面取得了更好的结果。

详细文章讲解见第4部分。

📚2 运行结果

部分代码：

function [Best_score,Best_pos]=FOX(SearchAgents_no,Max_iter,lb,ub,dim,fobj)

Best_pos=zeros(1,dim);

Best_score=inf; %change this to -inf for maximization problems

MinT=inf;

%Initialize the positions of search agents

X=initialization(SearchAgents_no,dim,ub,lb);

Distance_Fox_Rat=zeros(SearchAgents_no,dim);

%Convergence_curve=zeros(1,Max_iter);

l=0;% Loop counter

% both c1 and c2 have different range value

c1=0.18; % range of c1 is [0, 0.18]

c2=0.82; % range of c2 is [0.19, 1]

% Main loop

while l<Max_iter

for i=1:size(X,1)

% Return back the search agents that go beyond the boundaries of the search space

Flag4ub=X(i,:)>ub;

Flag4lb=X(i,:)<lb;

X(i,:)=(X(i,:).*(~(Flag4ub+Flag4lb)))+ub.*Flag4ub+lb.*Flag4lb;

% Calculate objective function for each search agent

fitness=fobj(X(i,:));

% Update Alpha

if fitness<Best_score

Best_score=fitness; % Update alpha

Best_pos=X(i,:);

end

end

a=2*(1-(l/Max_iter));

Jump=0;

dir=rand;

for i=1:size(X,1)

r=rand;

p=rand;

if r>=0.5

if p>0.18

Time(i,:)=rand(1,dim);

sps=Best_pos/Time(i,:);

Distance_S_Travel(i,:)=sps* Time(i,:);

Distance_Fox_Rat(i,:)=0.5* Distance_S_Travel(i,:);

tt=sum(Time(i,:))/dim;

t=tt/2;

Jump=0.5*9.81* t^2;

X(i,:)=Distance_Fox_Rat(i,:).*Jump * c1;

elseif p<=0.18

Time(i,:)=rand(1,dim);

sps=Best_pos/Time(i,:);

Distance_S_Travel(i,:)=sps* Time(i,:);

Distance_Fox_Rat(i,:)=0.5* Distance_S_Travel(i,:);

tt=sum(Time(i,:))/dim;

t=tt/2;

Jump=0.5*9.81* t^2;

X(i,:)=Distance_Fox_Rat(i,:).*Jump * c2;

end

if MinT>tt

MinT=tt;

end

elseif r<0.5

%random walk

ans(i,:)=Best_pos+randn(1,dim).*(MinT*a);% exploring new positions

X(i,:)=ans(i,:);

end

end

l=l+1;

end

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

🌈4 Matlab代码及文章讲解