切諾貝利災(zāi)難優(yōu)化器Chernobyl Disaster Optimizer (CDO)是H. Shehadeh于2023年提出的新型智能優(yōu)化算法。該方法是受到切爾諾貝利核反應(yīng)堆堆芯爆炸而來的啟發(fā)。在CDO方法中，放射性的發(fā)生是由于核的不穩(wěn)定性，核爆炸會(huì)發(fā)出不同類型的輻射。這些輻射中最常見的種類被稱為伽馬、貝塔和阿爾法粒子。算法主要圍繞三種粒子的更新方式展開。

01

經(jīng)作者查閱文獻(xiàn)發(fā)現(xiàn)，該方法其實(shí)與灰狼算法有很大的相似性，大家可以作為參考。在CEC2005函數(shù)中結(jié)果展示：

% CDO函數(shù)，該算法與灰狼算法很像
function [Alpha_score,Alpha_pos,Convergence_curve]=CDO(SearchAgents_no,Max_iter,lb,ub,dim,fobj)% initialize alpha, beta, and gamma particle positions (search radiations (Agents)) 
Alpha_pos=zeros(1,dim);
Alpha_score=inf; %change this to -inf for maximization problemsBeta_pos=zeros(1,dim);
Beta_score=inf; %change this to -inf for maximization problemsGamma_pos=zeros(1,dim);
Gamma_score=inf; %change this to -inf for maximization problems%Initialize the positions of search radiations (Agents)
Positions=initialization(SearchAgents_no,dim,ub,lb);Convergence_curve=zeros(1,Max_iter);l=0;% Loop counter% Main loop
while l<Max_iterfor i=1:size(Positions,1)  % Return back the search radiations (Agents) that go beyond the boundaries of the search spaceFlag4ub=Positions(i,:)>ub;Flag4lb=Positions(i,:)<lb;Positions(i,:)=(Positions(i,:).*(~(Flag4ub+Flag4lb)))+ub.*Flag4ub+lb.*Flag4lb;               % Calculate objective function for each search radiations (Agents)fitness=fobj(Positions(i,:));% Update Alpha, Beta, and Gamma - search radiations (Agents)if fitness<Alpha_score Alpha_score=fitness; % Update alphaAlpha_pos=Positions(i,:);endif fitness>Alpha_score && fitness<Beta_score Beta_score=fitness; % Update betaBeta_pos=Positions(i,:);endif fitness>Alpha_score && fitness>Beta_score && fitness<Gamma_score Gamma_score=fitness; % Update gammaGamma_pos=Positions(i,:);endenda=3-l*((3)/Max_iter); % a decreases linearly from 3 to 0 Equation(9)a1 = ((log10((16000-1)*rand(1,1)+16000)));a2 = ((log10((270000-1)*rand(1,1)+270000)));a3 = ((log10((300000-1)*rand(1,1)+300000)));  % Update the Position of search radiations (Agents)for i=1:size(Positions,1)for j=1:size(Positions,2)     %------------------- alpha------------------------------           r1=rand(); % r1 is a random number in [0,1]r2=rand(); % r2 is a random number in [0,1]pa=pi*r1*r1/(0.25*a1)- a*rand() ; % Equation (23)C1=r2*r2*pi; D_alpha=abs(C1*Alpha_pos(j)-Positions(i,j)); va=0.25*(Alpha_pos(j)-pa*D_alpha); % Equation (22)%------------------- Beta------------------------------           r1=rand();r2=rand();pb=pi*r1*r1/(0.5*a2)- a*rand()  ; % Equation (17)C2=r2*r2*pi; D_beta=abs(C2*Beta_pos(j)-Positions(i,j)); vb=0.5*(Beta_pos(j)-pb*D_beta); % Equation (16)      %------------------- Gamma ------------------------------           r1=rand();r2=rand(); py=(pi*r1*r1)/a3- a*rand() ; % Equation (11)C3=r2*r2*pi; D_gamma=abs(C3*Gamma_pos(j)-Positions(i,j));vy=Gamma_pos(j)-py*D_gamma; % Equation (10)             Positions(i,j)=(va+vb+vy)/3;% Equation (28)endendl=l+1;    Convergence_curve(l)=Alpha_score;
end

03

參考文獻(xiàn)

[1]H. Shehadeh.Chernobyl Disaster Optimizer (CDO): A Novel Metaheuristic Method for Global Optimization, Neural Computing and Applications. DOI: https://dx.doi.org/10.1007/s00521-023-08261-1

04代碼獲取方式：后臺(tái)回復(fù)關(guān)鍵詞：2023

免費(fèi)獲取2023年智能優(yōu)化算法合集matlab代碼