【改进灰狼优化算法】混沌灰狼优化算法(Matlab代码实现）

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

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

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

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

1.1 灰狼优化算法基本思想

1.2 灰狼捕食猎物过程

1.3 实现步骤及程序框图

1.4 混沌灰狼优化算法策略

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

1.1 灰狼优化算法基本思想

灰狼优化算法是一种群智能优化算法，它的独特之处在于一小部分拥有绝对话语权的灰狼带领一群灰狼向猎物前进。在了解灰狼优化算法的特点之前，我们有必要了解灰狼群中的等级制度。

灰狼群一般分为4个等级：处于第一等级的灰狼用α表示，处于第二阶级的灰狼用β表示，处于第三阶段的灰狼用δ表示，处于第四等级的灰狼用ω表示。按照上述等级的划分，灰狼α对灰狼β、δ和ω有绝对的支配权；灰狼ω对灰狼δ和ω有绝对的支配权；灰狼δ对灰狼ω有绝对的支配权。

1.2 灰狼捕食猎物过程

GWO 优化过程包含了灰狼的社会等级分层、跟踪、包围和攻击猎物等步骤，其步骤具体情况如下所示。

1.2.1 社会等级分层

当设计 GWO 时，首先需构建灰狼社会等级层次模型。计算种群每个个体的适应度，将狼群中适应度最好的三匹灰狼依次标记为 、、而剩下的灰狼标记为 。也就是说，灰狼群体中的社会等级从高往低排列依次为、、及 。GWO 的优化过程主要由每代种群中的最好三个解（即  、、）来指导完成。

1.2.2 包围猎物

灰狼群体通过以下几个公式逐渐接近并包围猎物:

式中，t是当前的迭代代数，A和C是系数向量，Xp和X分别是猎物的位置向量和灰狼的位置向量。A和C的计算公式如下:  
                               

式中，a是收敛因子，随着迭代次数从2线性减小到0，r1和r 2服从[ 0，1]之间的均匀分布。

1.2.3 狩猎

狼群中其他灰狼个体Xi根据α、β和百的位置Xa、XB和Xo来更新各自的位置：

式中,Da，Dβ和D6分别表示a，β和5与其他个体间的距离;Xa，Xβ和X6分别代表a，β和5的当前位置;C1，C2，C3是随机向量，X是当前灰狼的位置。
灰狼个体的位置更新公式如下:

1.2.4 攻击猎物

构建攻击猎物模型的过程中，根据2)中的公式，ａ值的减少会引起 A 的值也随之波动。换句话说，A 是一个在区间［-a，a］（备注：原作者的第一篇论文里这里是[-2a,2a]，后面论文里纠正为[-a,a]）上的随机向量，其中ａ在迭代过程中呈线性下降。当 A 在［-1，１］区间上时，则捜索代理（Search Agent）的下一时刻位置可以在当前灰狼与猎物之间的任何位置上。

1.2.5 寻找猎物

灰狼主要依赖  、、 的信息来寻找猎物。它们开始分散地去搜索猎物位置信息，然后集中起来攻击猎物。对于分散模型的建立，通过｜A｜＞１使其捜索代理远离猎物，这种搜索方式使 GWO 能进行全局搜索。GWO 算法中的另一个搜索系数是Ｃ。从2.2中的公式可知，Ｃ向量是在区间范围［０，２］上的随机值构成的向量，此系数为猎物提供了随机权重，以便増加（｜Ｃ｜＞１）或减少（｜Ｃ｜＜１）。这有助于 GWO 在优化过程中展示出随机搜索行为，以避免算法陷入局部最优。值得注意的是，Ｃ并不是线性下降的，Ｃ在迭代过程中是随机值，该系数有利于算法跳出局部，特别是算法在迭代的后期显得尤为重要。

1.3 实现步骤及程序框图

1.3.1 步骤

Step1:种群初始化:包括种群数量N，最大迭代次数Maxlter，调控参数a,A，C.Step2:根据变量的上下界来随机初始化灰狼个体的位置X。
Step3:计算每一头狼的适应度值，并将种群中适应度值最优的狼的位置信息保存，将种群中适应度值次优的狼的位置信息保存为，将种群中适应度第三优的灰狼的位置信息保存为。
Step4:更新灰狼个体X的位置。
step5:更新参数a，A和C。
Step6:计算每一头灰狼的适应度值，并更新三匹头狼的最优位置。
Step7:判断是否到达最大迭代次数Maxlter，若满足则算法停止并返回Xa的值作为最终得到的最优解，否则转到Step4。

1.3.2 程序框图

1.4 混沌灰狼优化算法策略

混沌具有随机性、遍历性和规律性等特点，目前已广泛应用于 PSO、DE 和ABC等群体智能优化算法中以提高算法的搜索效率。为了使初始种群个体尽可能地利用解空间的信息，本代码采用Logistic混沌映射产生混沌序列用来进行种群初始化，其数学模型为：

📚2 运行结果

部分代码：

function [Alpha_score,Alpha_pos,Convergence_curve]=CGWO(SearchAgents_no,Max_iter,lb,ub,dim,fobj)

% alpha, beta,delta
Alpha_pos=zeros(1,dim);
Alpha_score=inf; %change this to -inf for maximization problems

Beta_pos=zeros(1,dim);
Beta_score=inf; %change this to -inf for maximization problems

Delta_pos=zeros(1,dim);
Delta_score=inf; %change this to -inf for maximization problems

%Initialize the positions of search agents
Positions=initialization(SearchAgents_no,dim,ub,lb);

Convergence_curve=zeros(1,Max_iter);

l=0;% Loop counter

% Main loop
while l<Max_iter
    for i=1:size(Positions,1)  
        
       % Return back the search agents that go beyond the boundaries of the search space
        Flag4ub=Positions(i,:)>ub;
        Flag4lb=Positions(i,:)<lb;
        Positions(i,:)=(Positions(i,:).*(~(Flag4ub+Flag4lb)))+ub.*Flag4ub+lb.*Flag4lb;               
        
        % Calculate objective function for each search agent
        fitness=fobj(Positions(i,:));
        
        % Update Alpha, Beta, and Delta
        if fitness<Alpha_score 
            Alpha_score=fitness; % Update alpha
            Alpha_pos=Positions(i,:);
        end
        
        if fitness>Alpha_score && fitness<Beta_score 
            Beta_score=fitness; % Update beta
            Beta_pos=Positions(i,:);
        end
        
        if fitness>Alpha_score && fitness>Beta_score && fitness<Delta_score 
            Delta_score=fitness; % Update delta
            Delta_pos=Positions(i,:);
        end
    end
    
    
    a=2-l*((2)/Max_iter); % a decreases linearly fron 2 to 0
    
    % Update the Position of search agents including omegas
    for i=1:size(Positions,1)
        for j=1:size(Positions,2)     
                       
            r1=rand(); % r1 is a random number in [0,1]
            r2=rand(); % r2 is a random number in [0,1]
            
            A1=2*a*r1-a; % Equation (3.3)
            C1=2*r2; % Equation (3.4)
            
            D_alpha=abs(C1*Alpha_pos(j)-Positions(i,j)); % Equation (3.5)-part 1
            X1=Alpha_pos(j)-A1*D_alpha; % Equation (3.6)-part 1
                       
            r1=rand();
            r2=rand();
            
            A2=2*a*r1-a; % Equation (3.3)
            C2=2*r2; % Equation (3.4)
            
            D_beta=abs(C2*Beta_pos(j)-Positions(i,j)); % Equation (3.5)-part 2
            X2=Beta_pos(j)-A2*D_beta; % Equation (3.6)-part 2       
            
            r1=rand();
            r2=rand(); 
            
            A3=2*a*r1-a; % Equation (3.3)
            C3=2*r2; % Equation (3.4)
            
            D_delta=abs(C3*Delta_pos(j)-Positions(i,j)); % Equation (3.5)-part 3
            X3=Delta_pos(j)-A3*D_delta; % Equation (3.5)-part 3             
            
            Positions(i,j)=(X1+X2+X3)/3;% Equation (3.7)
            
        end
    end
    l=l+1;    
    Convergence_curve(l)=Alpha_score;
end

🎉3 参考文献

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

[1]晏福,徐建中,李奉书.混沌灰狼优化算法训练多层感知器[J].电子与信息学报,2019,41(04):872-879.

🌈4 Matlab代码实现