CEC2022: Osprey optimization algorithm (OOA) solves CEC2022 (provide MATLAB code

1. Introduction to Osprey Optimization Algorithm

The Osprey optimization algorithm (OOA) was proposed by Mohammad Dehghani and Pavel Trojovský in 2023, which simulates the predation behavior of an osprey.

The osprey is the only medium-sized raptor belonging to the order Eagleformes, the family Osprey, and the genus Osprey. Male and female are similar. Body length 51-64 cm, weight 1000-1750 grams. The head is white, with dark brown longitudinal stripes on the top of the head, and the feathers on the occiput are slightly elongated in a lanceolate shape, forming a short crest. On the side of the head there is a broad black band running from the base of the forehead through the eyes to the nape of the neck where it blends with the blackness of the nape. The upper body is dark brown with a slight purple sheen. The underparts are white, the dark vertical stripes on the chest and flight feathers, and the alternate horizontal spots on the tail feathers are very eye-catching. The iris is light yellow or orange-yellow, the naked skin around the eyes is lead yellow-green, the mouth is black, the cere is lead blue, the feet and toes are yellow, and the claws are black.

Ospreys inhabit lakes, rivers, coasts or open areas, especially like to move in river valleys in mountain forests or water areas with trees. It is common to fly around rivers, lakes, marshes and seashores. When it sees bait in the water, it will go straight to the surface of the water and grab it with its feet. The toe has sharp claws, the bottom of the toe is full of thin spines, and the outer toe can turn from front to back, which are very suitable for fishing. On a sunny day, hover over the water, dive down after a fixed point, and then bring the caught fish to rocks, electric poles, trees and other places to enjoy. Nests often camp on rocky reefs on coasts or islands. It mainly feeds on fish, and sometimes preys on other small terrestrial animals such as frogs, lizards, and small birds. In addition to Antarctica and the Arctic, Asia, North America and other continents are distributed.

1.1 Osprey optimization algorithm principle

The Osprey optimization algorithm consists of two stages: the first stage is for the Osprey to identify the position of the fish and catch fish (global exploration), and the second stage is to bring the fish to a suitable location (local mining). The detailed design is as follows:

1.1.1 Population initialization

采用下式随机初始化鱼鹰种群：
$$\begin{array}{c}X={\left[\begin{array}{c}{X}_1\\ ⋮\\ X_i\\ ⋮\\ X_N\end{array}\right]}_{N\times m}={\left[\begin{array}{ccccc}{x}_{1,1}& \cdots & x_{1,j}& \cdots & x_{1,m}\\ ⋮& \ddots & ⋮& \ddots & ⋮\\ x_{i,1}& \cdots & x_{i,j}& \cdots & x_{i,m}\\ ⋮& \ddots & ⋮& \ddots & ⋮\\ x_{N,1}& \cdots & x_{N,}& \cdots & x_{N,m}\end{array}\right]}_{N\times m},\\ x_{i,j}=l{b}_j+r_{i,j}\cdot \left(u{b}_j-l{b}_j\right),i=1,2,\dots ,N,j=1,2,\dots ,m,\end{array}$$
Among them, N is the number of ospreys, m is the dimension of the problem, and the fitness value is calculated according to the optimization problem after initializing the position:
$$F={\left[\begin{array}{c}{F}_1\\ ⋮\\ F_i\\ ⋮\\ F_N\end{array}\right]}_{N\times 1}={\left[\begin{array}{c}F\left(X_1\right)\\ ⋮\\ F\left(X_i\right)\\ ⋮\\ F\left(X_N\right)\end{array}\right]}_{N\times 1}$$

1.1.2 Global Exploration (Phase 1: Location Identification and Fishing)

Ospreys are powerful hunters, able to detect the location of fish underwater due to their powerful eyesight. After locating the fish, they attack it and hunt the fish by diving underwater. The first phase of population renewal in OOA was modeled based on simulations of this natural behavior of ospreys. Modeling osprey attacks on fish results in significant shifts in the osprey's position in the search space, which increases the exploration power of OOA in identifying optimal regions and escaping local optima. In the OOA design, for each osprey, the positions of other ospreys with better objective function values ​​in the search space are considered as underwater fish. The position of each Osprey is specified using the following equation.
FP i = { X k ∣ k ∈ { 1 , 2 , … , N } ∧ F k < F i } ∪ { X best } \boldsymbol{FP _ { i }}=\left\{X_{k} \mid \boldsymbol{k} \in\{1,2, \ldots, N\} \wedge \boldsymbol{F}_{k}<\boldsymbol{F}_{i}\right\} \cup\left\{ \boldsymbol{X}_{\text {best }}\right\}FPi​={ Xk​∣k∈{ 1,2,…,N}∧Fk​<Fi​}∪{ Xbest ​}
Among them,$F{P}_i$is the location set of the i-th osprey, $X_{\text{best }}$Position for the best Osprey.
The Osprey randomly detects the location of one of the fish and attacks it. Based on the simulation of the movement of the osprey towards the fish, the new position of the corresponding osprey is calculated using the following equation. This new position, if its value of the objective function is better, replaces the previous position of the Osprey.
$$\begin{array}{l}{x}_{i,j}^{P1}=x_{i,j}+r_{i,j}\cdot \left(S{F}_{i,j}-I_{i,j}\cdot x_{i,j}\right),\\ x_{i,j}^{P1}=\{\begin{array}{ll}{x}_{i,j}^{P1},& l{b}_j\le x_{i,j}^{P1}\le u{b}_j;\\ l{b}_j,& x_{i,j}^{P1}<l{b}_j;\\ u{b}_j,& x_{i,j}^{P1}>u{b}_j.\end{array}\\ X_i=\{\begin{array}{l}{X}_i^{P1},F_i^{P1}<F_i;\\ X_i,\text{ else },\end{array}\end{array}$$
Among them, $x_{i,j}^{P1}$is the new position of the i-th osprey in the first stage, $F_{i,j}^{P1}$is its corresponding fitness value. $S{F}_{i,j}$Is a random number between [0,1], $I_{i,j}$is a random number in the set {1,2}.

1.1.3 Local mining (second stage: bringing the fish to a suitable location)

After catching a fish, the osprey takes it to a suitable (safe for him) location and eats it there. The second phase of repopulation in OOA is based on simulation modeling of this natural behavior of ospreys. Modeling to bring the fish to a suitable location results in small changes in the osprey's position in the search space, leading to increased exploitation of OOA in local search and convergence to better solutions near the discovered solutions. In the design of OOA, in order to simulate this natural behavior of osprey, first, for each member of the population, use the following formula to calculate a new random position as the "suitable position for eating fish". Then, if the value of the objective function improves at this new location, the previous location of the corresponding osprey is replaced.
$$\begin{array}{c}{x}_{i,j}^{P2}=x_{i,j}+\frac{l{b}_j+r\cdot \left(u{b}_j-l{b}_j\right)}{t},i=1,2,\dots ,N,j=1,2,\dots ,m,t=1,2,\dots ,T,\\ x_{i,j}^{P2}=\{\begin{array}{l}{x}_{i,j}^{P2},l{b}_j\le x_{i,j}^{P2}\le u{b}_j;\\ l{b}_j,x_{i,j}^{P2}<l{b}_j\\ u{b}_j,x_{i,j}^{P2}>u{b}_j,\end{array}\\ X_i=\{\begin{array}{l}{X}_i^{P2},F_i^{P2}<F_i;\\ X_i,\text{ else },\end{array}\end{array}$$
Among them, $x_{i,j}^{P2}$is the new position of the i-th osprey in the second stage, $F_{i,j}^{P2}$is its corresponding fitness value. $r$ is a random number between [0,1],$t$ and$T$ are the current iteration number and the maximum iteration number respectively.

1.2 Algorithm description

1.3 Algorithm process

1.4 References

Dehghani Mohammad, Trojovský Pavel.Osprey optimization algorithm: A new bio-inspired metaheuristic algorithm for solving engineering optimization problems[J].Frontiers in Mechanical Engineering,2023,8.

2. CEC2022

CEC2022 has a total of 12 single-objective test functions, and the optional dimensions of each test function are: 10/20. The table below gives the theoretical optimal values ​​for each function. The CEC2022 test function is complex and extremely challenging.

3. Solution results

The complete code is added to the blogger’s WeChat at the bottom of the blog: djpcNLP123
The function dimension tested in this example is 10 (can be modified according to your own needs), the OOA population size is 100, and the maximum number of iterations is 50. Increasing the number of iterations will give better results~

close all
clear 
clc
SearchAgents_no=20; % 种群大小
Function_name=1; %测试函数1-12
Max_iteration=50; % 最大迭代次数
lb=-100;%变量下界
ub=100;%变量上界
dim=10;%维度 10/20
[fMin,bestX,curve]=OOA(SearchAgents_no,Max_iteration,lb,ub,dim,fobj);  
semilogy(curve,'Color','b','linewidth',2.5)
legend('OOA')
title(['cec2022-F' num2str(Function_name)])
xlabel('Iteration');
ylabel('Best score obtained so far');
axis tight
grid on
display(['The best solution obtained is : ', num2str(bestX)]);
display(['The best optimal value of the objective funciton is : ', num2str(fMin)]);

Some experimental results:

F1

F2

F3

F4

4. Reference code

The complete code is added to the blogger WeChat at the bottom of the blog: djpcNLP123