Stability analysis and optimal control of COVID-19 with quarantine and media awareness

Jiajia Zhang; Yuanhua Qiao; Yan Zhang; Jiajia Zhang; Yuanhua Qiao; Yan Zhang

doi:10.3934/mbe.2022230

Mathematical Biosciences and Engineering

2022, Volume 19, Issue 5: 4911-4932. doi: 10.3934/mbe.2022230

Previous Article Next Article

Research article Special Issues

Stability analysis and optimal control of COVID-19 with quarantine and media awareness

Faculty of Science, Beijing University of Technology, Beijing 100124, China

Academic Editor: Wandi Ding

Received: 08 December 2021 Revised: 27 February 2022 Accepted: 07 March 2022 Published: 14 March 2022

In this paper, an improved COVID-19 model is given to investigate the influence of treatment and media awareness, and a non-linear saturated treatment function is introduced in the model to lay stress on the limited medical conditions. Equilibrium points and their stability are explored. Basic reproduction number is calculated, and the global stability of the equilibrium point is studied under the given conditions. An object function is introduced to explore the optimal control strategy concerning treatment and media awareness. The existence, characterization and uniqueness of optimal solution are studied. Several numerical simulations are given to verify the analysis results. Finally, discussion on treatment and media awareness is given for prevention and treatment of COVID-19.

Keywords:

Citation: Jiajia Zhang, Yuanhua Qiao, Yan Zhang. Stability analysis and optimal control of COVID-19 with quarantine and media awareness[J]. Mathematical Biosciences and Engineering, 2022, 19(5): 4911-4932. doi: 10.3934/mbe.2022230

Related Papers:

[1]	Martha Garlick, James Powell, David Eyre, Thomas Robbins . Mathematically modeling PCR: An asymptotic approximation with potential for optimization. Mathematical Biosciences and Engineering, 2010, 7(2): 363-384. doi: 10.3934/mbe.2010.7.363
[2]	Ellina Grigorieva, Evgenii Khailov, Andrei Korobeinikov . Parametrization of the attainable set for a nonlinear control model of a biochemical process. Mathematical Biosciences and Engineering, 2013, 10(4): 1067-1094. doi: 10.3934/mbe.2013.10.1067
[3]	Paula Mercurio, Di Liu . Identifying transition states of chemical kinetic systems using network embedding techniques. Mathematical Biosciences and Engineering, 2021, 18(1): 868-887. doi: 10.3934/mbe.2021046
[4]	Xiaoxuan Pei, Kewen Li, Yongming Li . A survey of adaptive optimal control theory. Mathematical Biosciences and Engineering, 2022, 19(12): 12058-12072. doi: 10.3934/mbe.2022561
[5]	Huangjing Yu, Yuhao Wang, Heming Jia, Laith Abualigah . Modified prairie dog optimization algorithm for global optimization and constrained engineering problems. Mathematical Biosciences and Engineering, 2023, 20(11): 19086-19132. doi: 10.3934/mbe.2023844
[6]	Yancong Zhou, Chenheng Xu, Yongqiang Chen, Shanshan Li . C4 olefin production conditions optimizing based on a hybrid model. Mathematical Biosciences and Engineering, 2023, 20(7): 12433-12453. doi: 10.3934/mbe.2023553
[7]	Qing Wu, Chunjiang Zhang, Mengya Zhang, Fajun Yang, Liang Gao . A modified comprehensive learning particle swarm optimizer and its application in cylindricity error evaluation problem. Mathematical Biosciences and Engineering, 2019, 16(3): 1190-1209. doi: 10.3934/mbe.2019057
[8]	Rong Zheng, Heming Jia, Laith Abualigah, Shuang Wang, Di Wu . An improved remora optimization algorithm with autonomous foraging mechanism for global optimization problems. Mathematical Biosciences and Engineering, 2022, 19(4): 3994-4037. doi: 10.3934/mbe.2022184
[9]	Liping Wu, Zhongyi Xiang . A study of integrated pest management models with instantaneous and non-instantaneous impulse effects. Mathematical Biosciences and Engineering, 2024, 21(2): 3063-3094. doi: 10.3934/mbe.2024136
[10]	Andriy A. Avramenko, Igor V. Shevchuk . Renormalization group analysis of heat transfer in the presence of endothermic and exothermic chemical reactions. Mathematical Biosciences and Engineering, 2019, 16(4): 2049-2062. doi: 10.3934/mbe.2019100

Abstract

1. Introduction

It is prevalent to solve the optimization problems in our life, which is generally done by selecting a set of variables to reach an optimal value while satisfying a set of constraints. For instance, antenna arrays are a typical global optimization problem that is important for detecting and processing signals from different directions. The goal of antenna array synthesis is to generate a physical layout of the array in which the radiation pattern is as close as possible to the desired pattern. Unfortunately, most optimization problems are always nonlinear, non-differentiable, and multipolar. Numerous optimization methods have been proposed for solving these optimization problems. Classical optimization methods, including gradient methods and alternating projection methods, are constantly being proposed. Nevertheless, it is known that classical optimization methods require a fixed point that is fairly near to the final value, otherwise they are likely to get trapped in local optimal. Moreover, the quality of the solution strongly depends on the estimation of the initial values. However, the goodness of the final solution scheme of the algorithm mostly depends on the setting of the initial scheme. If the initial solution falls in the worse region of the solution, then the search will only be conducted to find the best solution in the worse solution region. Considering the shortcomings of these techniques, evolutionary algorithms are the best choice for solving large-scale optimization problems since evolutionary algorithms are population-based search algorithms with high convergence rates and have been successfully used for large-scale optimization problems extensively ^[1,2]. Therefore, different kinds of evolutionary algorithms, such as particle swarm optimization algorithms (PSO) ^[3,4], ant colony optimization (ACO), biogeography-based optimization (BBO) ^[5,6,7,8], animal migration optimization (AMO) ^[9,10], cuckoo search algorithm (CS) ^[11,12,13], and artificial bee colony (ABC) ^[14], have been adopted to handle antenna arrays, especially on large-scale optimization problems.

Recently, chemical reaction optimization (CRO) ^[15] is a population-based heuristic evolutionary algorithm that simulates molecular collisions in chemical reactions. The algorithm simulates the intermolecular interactions of molecules from a high-energy state into a low-energy state to explore the global optimal solution in the search space. To achieve efficient performance and extensive applications of the algorithm, different variations of the CRO algorithm were proposed successively to enhance the performance of the traditional CRO algorithm. In addition, there are many works that focus on the applications of CRO algorithms; for instance, Lam and Li ^[16] developed a real-coded CRO, called RCCRO for addressing continuous optimization problems. The experimental results indicate that the average performance of RCCRO is superior to other algorithms. In ^[17], the authors developed a double molecular structure-based chemical reaction optimization (DMSCRO) method and the CRO scheme is employed to outline the directed acyclic graph jobs scheduling in heterogeneous computing systems. Xu et al. ^[18] developed a hybrid chemical reaction optimization framework to solve the DAG-based task scheduling problem. It combines the novel heuristic approaches and a proposed selection strategy. Dam et al. ^[19] proposed a new method that integrates the chemical reaction optimization scheme with the unified tabu search heuristic method for solving the capacitated vehicle routing problem.

Lam and Li ^[20] solved the peer to peer live streaming with chemical reaction optimization and experimental results demonstrate that CRO beats many commonly-used techniques for containing population transitions in the system applications. In ^[21], Lam and Li used chemical reaction optimization to resolve the grid scheduling problem. Experimental results indicate that the proposed algorithm surpasses many existing evolutionary methods for the majority of cases tested. Truong et al. ^[22] introduced a chemical reaction optimization algorithm (CROG) coupled with a greedy strategy to tackle the quadratic allocation problem. The algorithm also proposes a design method for the update operator and a parameter tuning method. Unfortunately, it is difficult for these approaches to balance the exploitation and exploration capability of CRO algorithms simultaneously. Therefore, it is necessary to develop new optimization strategies for chemical reaction optimization to further enhance the performance of the algorithms.

To address these issues, we propose a new realistic encoded chemical reaction optimization algorithm that can boost the performance of traditional chemical reaction algorithms. First, inspired by the search mechanism in differential evolution algorithms, we propose a modified chemical reaction optimization (MCRO) model based on the optimal solution in the population. Then, to further improve the performance of the algorithm, the perturbation frequency of individuals is manipulated by increasing the modification rate of individuals and balance the exploitation and exploration ability of the algorithm. Moreover, to further accelerate the convergence of the algorithm, we assemble the stochastic initialization, opposition learning methods, and generalized opposition learning methods to obtain a new population initialization method. Finally, to demonstrate the superiority of the proposed algorithm, the algorithm is compared with several benchmark functions and applied on two antenna arrays. Based on the experimental results, it can be observed that our proposed algorithm is superior to other state-of-the-art algorithms from several perspectives.

2. Chemical reaction optimization

2.1. The basic concept of the CRO

Chemical Reaction Optimization (CRO) is an evolutionary algorithm first published by Lam and Li in ^[22,23], which simulates the collision of molecules in a chemical reaction. For each molecule, its molecular structure is similar to the optimization problem and the molecule denotes the solution to the problem under consideration. Potential energy ( $\chi$ ) is the objective function value and kinetic energy ( $\eta$ ) can be interpreted as the tolerance of the system, which accepts a solution that is worse than the current one. This approach can assist the solution to be partially replaced locally. Assuming that the numerator is $\omega$ , $\chi$ can be calculated as follows: $\chi_{\omega} = f(\omega)$ .

2.2. Elementary reactions

In a standard CRO, the algorithm simulates four different types of chemical reactions, comprising on-wall ineffective collision, decomposition, inter-molecular ineffective collision, and synthesis. They are mainly employed to control the redistribution of energy between molecules and buffers in the solution. Among them, single operations include wall ineffective collisions and disintegration, and the rest are multiple operations.

When a molecule strikes the wall of a container and then bounces away, the on-wall ineffective collision reaction is invoked. The process will have some change. The molecular structure $\omega$ and $\chi$ will change to a generated molecular $\omega^{\prime}$ in the neighborhood of $\omega$ if the condition is satisfied under Eq 1.

$\begin{equation} \chi_{\omega}+\eta_{\omega} \geq \chi_{\omega^{\prime}} \end{equation}$

(1)

After an invalid collision at the wall, a certain fraction of the $\eta$ of the transformed molecule is withdrawn to the buffer. Then, the $\eta$ of the new molecule $\omega^{\prime}$ can be as follows:

$\begin{equation} \eta_{\omega^{\prime}} = \left(\chi_{\omega}+\eta_{\omega}-\eta_{\omega}\right) \times a \end{equation}$

(2)

where $a \in[\eta^{LR}, 1]$ , $\eta^{LR}$ is a statistical parameter that limits the maximum percentage of $\eta$ that can be lost at each time. As long as there is enough $\chi$ at the beginning, it is possible for a molecule with lower $\chi$ to transform into a molecule with higher $\chi$ , corresponding to a poorer solution. Following the collision, the molecule has fewer $\eta$ .

The decomposition collision reaction results from a molecule hitting a wall, which then breaks one molecule into several molecules. Assume that $\omega$ produces $\omega_{1}^{\prime}$ and $\omega_{2}^{\prime}$ , i.e., $\omega \rightarrow \omega_{1}^{\prime}+\omega_{2}^{\prime}$ . Any mechanism, which can generate $\omega_{1}^{\prime}$ and $\omega_{2}^{\prime}$ from $\omega$ , is allowed. The algorithm considers two cases for the decomposition reaction: first, each molecule needs enough energy to complete the decomposition. The $\chi$ of the resulting molecule is denoted by the following:

$\begin{equation} \begin{array}{l} \eta_{\omega_{1}^{\prime}} = \left(\chi_{\omega}+\eta_{\omega}-\chi_{\omega_{1}^{\prime}}-\chi_{\omega_{2}^{\prime}}\right) \times q \\ \eta_{\omega_{2}^{\prime}} = \left(\chi_{\omega}+\eta_{\omega}-\chi_{\omega_{1}^{\prime}}-\chi_{\omega_{2}^{\prime}}\right) \times(1-q) \end{array} \end{equation}$

(3)

Based on the above analysis, the following conditions need to be satisfied for molecular decomposition:

$\begin{equation} \chi_{\omega}+\eta_{\omega} \geq \chi_{\omega_{1}^{\prime}}+\chi_{\omega_{2}^{\prime}} \end{equation}$

(4)

where $q$ is a random number that takes values in the range [0, 1]. Meanwhile, the algorithm considers another condition:

$\begin{equation} \begin{array}{l} \eta_{\omega_{1}^{\prime}} = \left(\chi_{\omega}+\eta_{\omega}-\chi_{\omega_{1}^{\prime}}-\chi_{\omega_{2}^{\prime}}+\text{ buffer }\right) \times k_{1} \times k_{2} \\ \eta_{\omega_{2}^{\prime}} = \left(\chi_{\omega}+\eta_{\omega}-\chi_{\omega_{1}^{\prime}}-\chi_{\omega_{2}^{\prime}}+\text{ buffer }\right) \times k_{3} \times k_{4} \end{array} \end{equation}$

(5)

$\begin{equation} \chi_{\omega}+\eta_{\omega}+\text{ buffer } \geq \chi_{\omega_{1}^{\prime}}+\chi_{\omega_{2}^{\prime}} \end{equation}$

(6)

where $k_{1}, k_{2}, k_{3}$ and $k_{4}$ are random numbers which take values in the range of [0, 1]. The buffer can be defined as below:

$\begin{equation} \text{ buffer } = \chi_{\omega}+\eta_{\omega}-\chi_{\omega_{1}^{\prime}}-\chi_{\omega_{2}^{\prime}}+bu f f e r-\eta_{\omega_{1}^{\prime}}-\eta_{\omega_{2}^{\prime}} \end{equation}$

(7)

Null collisions between molecules happen when two or more molecules collide with each other and then separate. The energy conservation condition can be expressed in the following way:

$\begin{equation} \chi_{\omega}+\eta_{\omega}+\eta_{\omega_{1}}+\eta_{\omega_{2}} \geq \chi_{\omega_{1}^{\prime}}+\chi_{\omega_{2}^{\prime}} \end{equation}$

(8)

Then the $\eta$ of two transformed molecules can be described as follows:

$\begin{equation} \begin{array}{l} \eta_{\omega_{0}^{\prime}} = \left(\chi_{\omega}+\eta_{\omega}+\eta_{\omega_{1}}+\eta_{\omega_{2}}-\chi_{\omega_{1}^{\prime}}-\chi_{\omega_{2}^{\prime}}\right) \times q \\ \eta_{\omega_{2}^{\prime}} = \left(\chi_{\omega}+\eta_{\omega}+\eta_{\omega_{1}}+\eta_{\omega_{2}}-\chi_{\omega_{1}^{\prime}}-\chi_{\omega_{2}^{\prime}}\right) \times(1-q) \end{array} \end{equation}$

(9)

where $q$ is a random number that takes values in the range [0, 1].

The synthesis is the opposite of decomposition. Synthesis occurs when a molecule collides and bond together. The reaction conditions that allow the synthesis reaction are description below:

$\begin{equation} \chi_{\omega}+\eta_{\omega}+\eta_{\omega_{1}}+\eta_{\omega_{2}} \geq \chi_{\omega^{\prime}} \end{equation}$

(10)

Then the $\eta$ for the new molecule is considered by the following method.

$\begin{equation} \eta_{\omega^{\prime}} = \chi_{\omega}+\eta_{\omega}+\eta_{\omega_{1}}+\eta_{\omega_{2}}-\chi_{\omega^{\prime}} \end{equation}$

(11)

To address the continuous optimization problem more efficiently, Lam and Li proposed a real coded version of CRO, namely RCCRO. Subsequently, a similar modified neighborhood search approach was used. Among them, the Gaussian distribution model was proposed to generate perturbations to search continuous space neighborhood. Their update method can be characterized here below:

$\begin{equation} \omega^{\prime}(i) = \omega(i)+N\left(0, \sigma^{2}\right) \end{equation}$

(12)

where $\sigma$ holds the step size and determines the efficiency of the algorithm. Because different problems have different characteristics, various $\sigma$ should be adopted. In the algorithm, the neighborhood method will be employed for one-wall and intermolecular ineffective collision. The framework of the algorithm is shown in Figure 1. In summary, the CRO algorithm starts with the initialization. Then if it is an inter-molecular ineffective collision, the single molecule reaction is employed. Otherwise we will carry out the bimolecule reaction. In terms of the single molecule reaction, if the composition condition is satisfied, then it will conduct the decomposition, otherwise, it will conduct the on-wall ineffective collision. In terms of the bimolecular reaction, if the synthesis condition is met, the synthesis reaction is conducted, otherwise the inter-molecular ineffective collision reaction is conducted. After these steps, new minimum points are checked. If the stopping criterion is not matched, then the above steps will be repeated. Otherwise, the global minimum point is obtained.

Figure 1. The framework of the CRO algorithm.

DownLoad: Full-Size Img PowerPoint

3. Modified chemical reaction optimization (MCRO)

3.1. Initial population

For evolutionary algorithms, population initialization is a crucial component. The goodness of population initialization directly influences the final convergence speed of the algorithm and the quality of the optimal solution. In fact, most algorithms use random initialization to generate the initial population in the absence of a priori information. Moreover, it is difficult to guarantee the quality of initialized population solutions by individual population initialization methods. Therefore, this paper proposes a new initialization method to generate the initial population by fusing random initialization, opposition-based learning initialization ^[24] and generalized opposition learning initialization ^[25]. This method can enhance the performance of CRO by preventing premature convergence to local optimal. Table 4 summarized the framework of this method. Figure 2 depicted the difference between the random population initiation and the opposition-based population initialization discussed above.

Figure 2. The difference between the random population initiation and above mentioned opposition based population initialization.

DownLoad: Full-Size Img PowerPoint

Procedure A new initialization approach
Begin
Randomly initialize each individual in population x.
Set the individual counter i=1, j=1
for $i=1 \rightarrow NP$ do
for $j=1 \rightarrow D$ do
$o x_{i, j}=x_{\min, j}+x_{\max, j}-x_{i, j}$
end for
end for
Set the individual counter $i=1$ , $j=1$
for $i=1 \rightarrow NP$ do
for $j=1 \rightarrow D$ do
$oo x_{i, j}=x_{\min, j}+x_{\max, j}-x_{i, j}$
end for
end for
The fitness of each individual was assessed and the most suitable individual for NP from $\{X(N P) \cup O X(N P) \cup O O X(N P)\}$ was selected as the initial population.
End

| Show Table

DownLoad: CSV

3.2. A modified search equation

Differential evolution algorithm (DE) ^[26] is an efficient stochastic search method which was originally proposed mainly for solving continuous optimization problems. Its simple structure and ease of completion allowed the algorithm to be quickly applied to different optimizations and achieve impressive performance. The algorithm also has the identical evolutionary process as other evolutionary algorithms. In the algorithm, we can see that the algorithm consists of three main processes: mutation, crossover and selection. At the beginning of the algorithm, a mid-autumn is generated using a random initialization, followed by the operations of mutation, crossover and selection to generate a new population. It is worth noting that several different variation operations have been proposed in the DE algorithm, such as:

$\begin{equation} \text{ DE/best } / 1: V_{i} = X_{\text{best }}+F \cdot\left(X_{r 1}-X_{r 2}\right) \end{equation}$

(13)

where $i \in\{1, \cdots, N P\}$ and $r_{1}$ and $r_{2}$ are random integer indices selected from $\{1, \cdots, N P\}$ that are different from each other. $F$ is a scaling factor or amplification factor. $X_{\text{best }}$ , the base vector to be perturbed, is the best individual in the current population and indicates the best information in the population.

In fact, the information about the best individuals in the current population is very useful for the whole population, and if the information is used wisely it is possible to accelerate the convergence rate of the population. For the strategy "DE/best/1", the best solution can explore the movement of the current population. Based on DE algorithm and the property of CRO, we propose a novel search mechanism to improve CRO:

$\begin{equation} {v_{i j} = x_{bestj}+\varphi_{ij}\left(x_{ij}-x_{bj}\right)} \end{equation}$

(14)

where $\varphi_{i j}$ denotes the scaling factor that controls the speed of optimization of the algorithm. It is calculated at the beginning of each generation. If $rand < rand$ , then $\varphi_{i j}$ is equal to $5 \times ({rand-rand }){\times}(rand-rand)$ , otherwise, $\varphi_{i j}$ equals to $randn$ .

We can note that " $rand$ " denotes a random number generated from the range of $[-1, 1], i \in\{1, 2, \cdots, NP\}$ , and $j \in\{1, 2, \cdots, D\}$ are randomly chosen indexes, " $randn$ " denotes a number generated through a Gaussian distribution with mean "0" and standard deviation "1".

To control the frequency of perturbation, it makes the algorithm can solve the difficult problem better. A control parameter, the modification rate (MR), is proposed, which serves mainly to generate a candidate position $v_{i j}$ from the current memory of $x_{i j}$ as described below:

$\begin{equation} v_{i j} = \left\{\begin{array}{ll} x_{\text{best }}+\varphi_{i j}\left(x_{i j}-x_{b j}\right), & \text{ if } {r_{ij}} \leq M R \\ x_{i j} & \text{ otherwise } \end{array}\right. \end{equation}$

(15)

where $\mathrm{r}_{\mathrm{ij}}$ denotes a random number generated from the range of $[-1, 1]$ . This new rule is similar with the binomial crossover of the DE algorithm, which can enhance the potential diversity of the population.

3.3. Boundary restrictions

For evolutionary algorithms, the setting of bounds is crucial. During the search process, the individuals in the part of MCRO algorithm can easily move to the boundary of the whole search space, which renders many of the solutions invalid. Therefore, it is important to restrict the individuals around the boundary to the search space once again. In other words, if most of the individuals are on the boundary, it is difficult for the algorithm to find the optimal solution to the problem, which will necessarily make the algorithm fall into local optimal solutions and make the whole population very easy to lose the population diversity. To address this problem, we propose to use the following rule for the setting of off-boundary individuals:

$\begin{equation} x_{i} = \left\{\begin{array}{l} x_{min}+\bmod \left(\left(x_{min}-x_{i}\right), \left(x_{i}-x_{min}\right)\right) \\ x_{max}-\bmod \left(\left(x_{i}-x_{max}\right), \left(x_{max}-x_{min}\right)\right) \end{array}\right. \end{equation}$

(16)

3.4. Time complexity analysis

In this section, the time complexity analysis of the proposed MCRO algorithm for the benchmark function optimization is analyzed. The time mainly depends on the searching loop. At first, the time complexity of the initialization method is $O(NP\times D)$ , where $NP$ is the population size and $D$ is the number of variables. Then for each iteration in the loop, the time is mainly spent in evaluating the fitness of each individual, which also costs $O(NP \times D)$ . Therefore, the time complexity of MCRO is $O((I+1)(NP\times D))$ , where $I$ is the number of iterations for the optimization process.

4. Results

4.1. The benchmark functions optimization by MCRO

To validate the effectiveness of the MCRO algorithm proposed in this paper, we first applied the algorithm to solve nine standard global optimization algorithms that are widely used. To ensure the fairness of the algorithm comparison, we did not modify any of the problems, which are tabulated in Table 1. As can be observed from the table, the first five problems are unimodal global optimization problems, where f05 is a noisy quadratic function. The remaining problems are all multimodal global optimization problems. It is worth noting that the local minima of these problems are increasing as the size of the problem increases, which is the reason why the difficulty of solving the problem increases with the size of the problem.

Table 1. Benchmark functions based in our experimental study.

Test Function	D	Range	optimum
$f_{01}=\sum_{i=1}^{n} x_{i}^{2}$	30	[-100, 100]	0
$f_{02}=\sum_{i=1}^{n}\left\|x_{i}\right\|+\prod_{i=1}^{n}\left\|x_{i}\right\|$	30	[-10, 10]	0
$f_{03}=\max _{i}\left\{x_{i} \mid, 1 \leq i \leq D\right\}$	30	[-100, 100]	0
$f_{04}=\sum_{i=1}^{D}\left(\left\lfloor x_{i}+0.5\right\rfloor^{2}\right.$	30	[-100, 100]	0
$f_{05}=\sum_{i=1}^{D} i x_{i}^{4}+random[0, 1)$	30	[-1.28, 1.28]	0
$f_{06}=\frac{1}{400} \sum_{i-1}^{D} x_{i}^{2}-\prod_{i=1}^{D} \cos \left(\frac{x_{i}}{\sqrt{i}}\right)+1$	30	[-32, 32]	0
$f_{07}=-20 \exp \left(-0.2 \sqrt{\frac{1}{D} \sum_{i-1}^{D} x_{i}^{2}}\right)-\exp \left(\frac{1}{D} \sum_{i=1}^{D} \cos 2 \pi x_{i}\right)+20+e$	30	[-600, 600]	0
$f_{08}=\frac{\pi}{D}\left\{10 \sin ^{2}\left(\pi y_{i}\right)+\sum_{i=1}^{D-1}\left(y_{i}-1\right)^{2}\left[1+10 \sin ^{2}\left(\pi y_{i}+1\right)\right]\right.$ $\left.+(y D-1)^{2}+\sum_{i=1}^{D} u\left(x_{i}, 10, 100, 4\right)\right\}$ $y_{i}=1+\frac{x_{i}+1}{4} \quad u\left(x_{i}, a, k, m\right)=\left\{\begin{array}{ll}k\left(x_{i}-a\right)^{m} & x_{i} > a\\ 0 & -a < x_{i} < a\\ k\left(-x_{i}-a\right)^{m} & x_{i} < -a\end{array}\right.$	30	[-50, 50]	0
$f_{09}=0.1\left\{10 \sin ^{2}\left(\pi y_{i}\right)+\sum_{i=1}^{D-1}\left(y_{i}-1\right)^{2}\left[1+10 \sin ^{2}\left(\pi y_{i}+1\right)\right]+(y D-1)^{2}\right\}$ $+\sum_{i=1}^{D} u\left(x_{i}, 10, 100, 4\right)$	30	[-50, 50]	0

| Show Table

DownLoad: CSV

To illustrate the superiority of the MCRO algorithm proposed in our paper, we compared the MCRO algorithm with the standard CRO, the Differential Evolution algorithm (DE), the Adaptive Differential Evolution (jDE) ^[27] and the Modified Artificial Bee Colony algorithm (MABC) ^[28]. Across the experiments, the algorithms were computed with a maximum number of functions of 150, 000 and 50 independent experiments were conducted for each problem, and their average values were finally compared as the final results. To compare fairly, the parameter of different algorithm can be set as follows: for the CRO algorithm, the population size NP is 20, KELossRate is 0.1, Molecoll is 0.2, InitialKE is 1000, and the value of decThres and synThres is 150000 and 0, respectively. The value of MR is 0.2. For the DE algorithm, the population size is 100, F is 0.3 and CR is 0.7. The parameter of the jDE is the same with the paper^[27]. For modified artificial bee colony algorithm (MABC), the parameters of this algorithm are identical to the paper ^[28]. The results are listed in Table 2. From Table 2, we can observe that the results of MCRO are better than the other algorithms while Figure 3 shows the 2D plot of f07.

Table 2. Performance comparisons of DE, jDE, MABC, CRO and MCRO.

F	DE	jDE	MABC	CRO	MCRO
f01	6.1695e-14 $\pm$ 4.7200e-14	2.0141e-28 $\pm$ 4.1612e-28	1.9244e-39 $\pm$ 4.3030e-39	6.7348e-007 $\pm$ 2.1237e-007	5.5207e-71 $\pm$ 1.0949e-70
f02	4.2283e-07 $\pm$ 3.0823e-07	1.6421e-17 $\pm$ 1.0644e-17	6.8088e-18 $\pm$ 1.5132e-17	0.00213.6637e-004	3.1288e-43 $\pm$ 3.9422e-43
f03	0.0458 $\pm$ 0.0809	0.8695 $\pm$ 0.7707	13.3649 $\pm$ 3.6487	0.0097 $\pm$ 0.0035	2.7781e-10 $\pm$ 2.2871e-10
f04	0 $\pm$ 0	0 $\pm$ 0	0 $\pm$ 0	0 $\pm$ 0	0 $\pm$ 0
f05	0.0096 $\pm$ 0.0021	0.0063 $\pm$ 0.0014	0.0311 $\pm$ 0.0083	0.0036 $\pm$ 0.0013	0.0083 $\pm$ 0.0023
f06	1.3524e-13 $\pm$ 7.7685e-14	0 $\pm$ 0	2.9016e-08 $\pm$ 6.4883e-08	4.7687e-07 $\pm$ 1.5030e-007	0 $\pm$ 0
f07	7.8443e-08 $\pm$ 3.1063e-08	7.9936e-15 $\pm$ 0	2.8747e-06 $\pm$ 6.4281e-06	0.0024 $\pm$ 4.2876e-004	1.6520e-14 $\pm$ 5.9448e-15
f 08	3.5919e-15 $\pm$ 2.2490e-15	4.6115e-30 $\pm$ 5.3009e-30	2.3631e-07 $\pm$ 5.2841e-07	0.1371 $\pm$ 0.3257	1.6222e-32 $\pm$ 7.0698e-34
f 09	4.7192e-14 $\pm$ 3.9138e-14	5.5873e-29 $\pm$ 7.6744e-29	1.4035e-09 $\pm$ 3.1382e-09	7.5513e-006 $\pm$ 5.3363e-006	7.0479e-31 $\pm$ 7.6213e-31

| Show Table

DownLoad: CSV

Figure 3. 2D Ackley function.

DownLoad: Full-Size Img PowerPoint

4.2. Array pattern synthesis using MCRO

To further demonstrate the effectiveness of the MCRO algorithm proposed in this paper, the proposed algorithm is applied to two different antenna array problems, including reconfigurable Antenna Array and linear antenna. {To enable a fair comparison of the algorithms}, the same number of function calculations is necessary. Therefore, in this paper, we use the same number of function evaluations for each algorithm as the final termination condition. As a result, we can evaluate the performance of the algorithms by the results of different algorithms and finally illustrate the efficiency of the proposed algorithm.

4.2.1. Synthesis of reconfigurable antenna-array

In this section, we applied our proposed MCRO to address the reconfigurable dual-beam antenna array. In this problem, when the phase distribution of the array is suitably modified, the amplitude distribution can produce pencil-shaped or sector power patterns. In addition, to reduce the effect of coupling between elements, an additional term is incorporated in the objective function ^[29]. Therefore, we can optimize the amplitude-excitation dynamic range (ARD) to minimize the mutual coupling problem ^[30,31]. Then, the objective function to be optimized in this paper is defined as below.

$\begin{equation} E c(P) = \sum\limits_{i = 1}^{3}\left(P_{i, d}^{(p)}-P_{i}^{(i)}\right)^{2}+\sum\limits_{i = 1}^{4}\left(P_{i, d}^{(s)}-P_{i}^{(s)}\right)^{2}+A D R \end{equation}$

(17)

where ADR denotes the amplitude-dynamic ratio. Based on the above description, it can be noticed that the algorithm can motivate the difference between the amplitudes by optimizing the ADR. Meanwhile, we can reveal that the difference between ADR and amplitude is linear, and minimizing the ADR is accompanied by a smaller difference between amplitudes. Therefore, in this paper, we optimize the reconfigurable antenna-array by optimizing the effect of coupling.

In order to efficiently design the above reconfigurable antenna array problem, we use our algorithm MCRO algorithm to tackle the problem. We compare the MCRO algorithm with some other algorithms including G3-GA, DE, jDE, MABC and MCRO. In the MCRO algorithm, the population size is 20 and the number of function computations is 20000, and the value of MR is 0.7. The rest of the parameters are the same as those of the previous global optimization algorithm. The experimental results are tabulated in Tables 3 and 4, mainly to determine the amplitude and phase excitation patterns for the dual-beam optimization. From these tables, we can see that the MCRO algorithm is superior to other algorithms, especially for reconfigurable antenna array designs with coupling effects. For the problems without coupling effects, we can find that most of these algorithms can find the optimal solutions. For problems with coupling effects, we can find that MCRO gives the best ADR and fitness values. These results show that our proposed algorithm is superior to other evolutionary algorithms. Figures 4 and 5 show the optimized excitation pattern and the dual-beam pattern, respectively. Figure 5 shows the design parameters that satisfy both pencil and sector beams. Figures 6 and 7 show the same case as Figures 4 and 5.

Table 3. The results of Experimental I and Experiment I with ADR.

Element Number	Experiment I		Experiment I with ADR
Element Number	Amplitude	Phase[deg.]	Amplitude	Phase[deg.]
1/20	0.129861	-10.6276	0.149815	-0.80199
2/19	0.174584	-25.4484	0.149813	-38.653
3/18	0.255019	-36.7229	0.154374	-29.1362
4/17	0.331274	-64.5232	0.256851	-71.1288
5/16	0.367296	87.84077	0.307724	-72.5951
6/15	0.556305	108.9342	0.412374	-97.438
7/14	0.652548	-80.0141	0.483299	97.30298
8/13	0.737711	-87.0023	0.564111	53.28282
9/12	0.831561	-13.7764	0.616299	-106.146
10/11	0.868279	43.79851	0.645976	88.67032
ADR	6.68		4.31
Fitness value	0.16		0.04

| Show Table

DownLoad: CSV

Table 4. Performance comparison of MCRO with G3-GA, DE, jDE, MABC and MCRO.

	Exp-I without ADR	Exp-I with ADR
	fitness	ADR	fitness
G3-GA	0.16	4.4137	0.1028
DE	0.16	4.3470	0.04
jDE	0.36	4.40	0.10
MABC	0.16	4.31	0.06
MCRO	0.16	4.31	0.04

| Show Table

DownLoad: CSV

Figure 4. Amplitude and phase excitation without coupling effect.

DownLoad: Full-Size Img PowerPoint

Figure 5. Dual-beam array patterns without coupling effect.

DownLoad: Full-Size Img PowerPoint

Figure 6. Dual-beam array patterns with coupling effect.

DownLoad: Full-Size Img PowerPoint

Figure 7. Dual-beam array patterns with coupling effect.

DownLoad: Full-Size Img PowerPoint

4.3. Synthesis of linear antenna

In this section, we applied our proposed MCRO to address the linear antenna. To design of a linear antenna array ^[32], let us assume that we have 2N isotropic radiators placed symmetrically along x-axis. The objective function of side lobe suppression is:

$\begin{equation} f_{1} = \sum\limits_{i} \frac{1}{\Delta \phi_{i}} \int_{\phi_{i i}}^{\phi_{i i}}\left|A F(\phi)^{2}\right| d \phi \end{equation}$

(18)

And for null control:

$\begin{equation} f_{2} = \sum\limits_{k}\left|A F\left(\phi_{k}\right)\right|^{2} d \phi \end{equation}$

(19)

{where $\Delta\phi_{i}$ is the bandwidth to be restrained} as $\phi_{u i}-phi_{l i}, \phi_{k}$ is the null direction. To optimize both objectives simultaneously, we sum Eqs 18 and 19 and complete their sum to make the objective function that our algorithm needs to optimize.

In order to elaborate the efficiency of the proposed MCRO algorithm, we came to test the reduction of the level of the side lobes in the linear array. Concurrently, we ran each algorithm independently for 25 times in each problem. The termination condition for each run of each algorithm was $5 \times 10^{4}$ number of function evaluations. All other parameters of the algorithms were consistent with the previous experiments. Meanwhile, the same number of function evaluations is applied to the other algorithms to construct a fair experiment.

In this section, we design an array of 26 elements that has minimum SLL in bands [ $0^{\circ}$ , $82^{\circ}$ ] and [ $98^{\circ}$ , $180^{\circ}$ ] and null direction in $20^{\circ}$ . Tables 5 and 6 summarize the results of the eight different evolutionary algorithms. From these tables, it can be viewed that our proposed MCRO algorithm is to be superior to other algorithms. Also, summarizes the array patterns of MCRO algorithm with other algorithms. From , it can be seen that MCRO can suppress the generation of side lobes to the maximum extent. Moreover, it also obtains the minimum gain in the $20^{circ}$ null.

Table 5. Geometry of the 26 element linear array, normalized numbers with respect to

$\lambda/2$ and null at

$20^{o}$ .

MCRO	0.301	0.939	1.98	2.012	3.367	3.541	5.284	7.049	4.590	8.021	6.093	10.881	9.411
MABC	0.330	0.979	1.821	2.039	3.267	2.897	4.151	6.574	7.749	4.880	5.799	9.544	10.909
DE	0.472	1.049	2.090	2.058	2.897	3.799	4.825	6.318	5.373	7.475	8.61	10.231	11.558
jDE	0.338	1.252	1.527	1.527	2.743	3.077	4.341	6.834	8.475	5.132	10.25	7.954	11.614

| Show Table

DownLoad: CSV

Table 6. Mean final objective function value, standard deviation, best, worst, median, and the Rank.

Algorithm	MCRO	DE	jDE	MABC
Mean	0.01194	0.02763	0.02356	0.0177
Best	0.00921	0.01984	0.01736	0.0131
Median	0.01219	0.02915	0.02319	0.0167
Worst	0.01488	0.03396	0.02982	0.0224
std	0.00209	0.00368	0.00443	0.0039
Rank	1	4	3	2

| Show Table

DownLoad: CSV

Figure 8. 13-element array for minimum SLL [0°, 82°] and [98°, 180°] and NULL 20°.

DownLoad: Full-Size Img PowerPoint

5. Conclusions

In this paper, a new evolutionary search strategy is employed to enhance the performance of the algorithm. Moreover, we present a new population initialization by fusing stochastic initialization, opposition learning method, and generalized opposition learning method. It is worth noting that our proposed algorithm is simple in structure and easy to implement. To demonstrate the superiority of our proposed MCRO algorithm, we used nine benchmark functions as test functions. The results show that the proposed algorithm significantly outperforms the other algorithms. In addition, the algorithm was tested on the synthesis problem of two antenna arrays. We analyzed and evaluated the experimental results and compared them with other evolutionary optimization algorithms. In terms of the algorithm's performance, our algorithm is significantly better than the other algorithms. In the future, we would like to apply our proposed algorithm MCRO to solve other antenna arrays, such as circular antenna arrays. More quantifiable experiments on those arrays are expected. Besides, we will focus on exploring other strategies such as multiobjective mechanism to enhance the efficiency of the algorithm. Meanwhile, other evaluation criteria will be computed to estimate the performance of those extended algorithms.

Acknowledgments

The work described in this paper was substantially supported by the National Natural Science Foundation of China under Grant No. 62076109.

Conflict of interest

The authors declare there is no conflicts of interest.

References

[1]	The Government of Wuhan Homepage. Available from: http://english.wh.gov.cn/.
[2]	A. E. Gorbalenya, S. C. Baker, R. S. Baric, R. J. de Groot, C. Drosten, A. A. Gulyaeva, et al., The species severe acute respiratory syndrome-related coronavirus: Classifying 2019-nCoV and naming it SARS-CoV-2, Nat. Microbiol., 5 (2020), 536–544. https://doi.org/10.1038/s41564-020-0695-z doi: 10.1038/s41564-020-0695-z
[3]	L. Zhou, M. Fan, Dynamics of an SIR epidemic model with limited medical resources revisited, Nonlinear Anal. Real World Appl., 13 (2012), 312–324. https://doi.org/10.1016/j.nonrwa.2011.07.036 doi: 10.1016/j.nonrwa.2011.07.036
[4]	N. Chitnis, J. M. Hyman, J. M. Cushing, Determining important parameters in the spread of Malaria through the sensitivity analysis of a mathematical model, Bull. Math. Biol., 70 (2008), 1272–1296. https://doi.org/10.1007/s11538-008-9299-0 doi: 10.1007/s11538-008-9299-0
[5]	T. K. Kar, P. K. Mondal, Global dynamics and bifurcation in delayed SIR epidemic model, Nonlinear Anal. Real World Appl., 12 (2011), 2058–2068. https://doi.org/10.1016/j.nonrwa.2010.12.021 doi: 10.1016/j.nonrwa.2010.12.021
[6]	C. Shan, H. Zhu, Bifurcations and complex dynamics of an SIR model with the impact of the number of hospital beds, J. Differ. Equations, 257 (2014), 1662–1688. https://doi.org/10.1016/j.jde.2014.05.030 doi: 10.1016/j.jde.2014.05.030
[7]	J. R. Artalejo, A. Economou, M. J. Lopez-Herrero, The stochastic SEIR model before extinction: computational approaches, Appl. Math. Comput., 265 (2015), 1026–1043. https://doi.org/10.1016/j.amc.2015.05.141 doi: 10.1016/j.amc.2015.05.141
[8]	P. Diaz, P. Constantine, K. Kalmbach, E. Jones, S. Pankavich, A modified SEIR model for the spread of Ebola in Western Africa and metrics for resource allocation, Appl. Math. Comput., 324 (2018), 141–155. https://doi.org/10.1016/j.amc.2017.11.039 doi: 10.1016/j.amc.2017.11.039
[9]	R. Almeida, Analysis of a fractional SEIR model with treatment, Appl. Math. Lett., 84 (2018), 56–62. https://doi.org/10.1016/j.aml.2018.04.015 doi: 10.1016/j.aml.2018.04.015
[10]	S. Jana, P. Haldar, T. K. Kar, Mathematical analysis of an epidemic model with isolation and optimal controls, Int. J. Comput. Math., 94 (2017), 1318–1336. https://doi.org/10.1080/00207160.2016.1190009 doi: 10.1080/00207160.2016.1190009
[11]	D. Gao, N. Huang, Optimal control analysis of a tuberculosis model, Appl. Math. Model., 58 (2018), 47–64. https://doi.org/10.1016/j.apm.2017.12.027 doi: 10.1016/j.apm.2017.12.027
[12]	D. K. Das, S. Khajanchi, T. K. Kar, The impact of the media awareness and optimal strategy on the prevalence of tuberculosis, Appl. Math. Comput., 366 (2020), 124732. https://doi.org/10.1016/j.amc.2019.124732 doi: 10.1016/j.amc.2019.124732
[13]	T. K. Kar, S. K. Nandi, S. Jana, M. Mandal, Stability and bifurcation analysis of an epidemic model with the effect of media, Chaos Solitons Fractals, 120 (2019), 188–199. https://doi.org/10.1016/j.chaos.2019.01.025 doi: 10.1016/j.chaos.2019.01.025
[14]	T. Chen, J. Rui, Q. Wang, Z. Zhao, J. Cui, L. Yin, A mathematical model for simulating the phase-based transmissibility of a novel coronvirus, Infect. Dis. Poverty, 9 (2020), 24. https://doi.org/10.1186/s40249-020-00640-3 doi: 10.1186/s40249-020-00640-3
[15]	M. A. Khan, A. Atangana, Modeling the dynamics of novel coronavirus (2019-nCoV) with fractional derivative, Alexandria Eng. J., 59 (2020), 2379–2389. https://doi.org/10.1016/j.aej.2020.02.033 doi: 10.1016/j.aej.2020.02.033
[16]	B. Tang, X. Wang, Q. Li, N. L. Bragazzi, S. Tang, Y. Xiao, et al., Estimation of the transmission risk of 2019-nCoV and its implication for public health interventions, J. Clin. Med., 9 (2020), 462. https://doi.org/10.3390/jcm9020462 doi: 10.3390/jcm9020462
[17]	B. Tang, N. L. Bragazzi, Q. Li, S. Tangd, Y. Xiao, J. Wu, An updated estimation of the risk of transmission of the novel coronavirus (2019-nCoV), Infect. Dis. Model., 5 (2020), 248–255. https://doi.org/10.1016/j.idm.2020.02.001 doi: 10.1016/j.idm.2020.02.001
[18]	P. Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
[19]	L. S. Pontryagin, V. G. Boltyanski, R. V. Gamkrelidze, E. F. Misenko, The Mathematical Theory of Optimal Process, Wiley, New Jersey, 1962.
[20]	H. R. Joshi, Optimal control of an HIV immunology model, Optim. Control Appl. Methods, 23 (2002), 199–213. https://doi.org/10.1002/oca.710 doi: 10.1002/oca.710
[21]	W. H. Fleming, R. W. Rishel, Deterministic and Stochastic Optimal Control, Springer, New York, 1975. https://doi.org/10.1007/978-1-4612-6380-7
[22]	D. L. Lukes, Differential Equations Classical to Controlled, Academia Press, New York, 1982.
[23]	H. Zhang, Z. Yang, K. A. Pawelek, S. Liu, Optimal control strategies for a two-group epidemic model with vaccination-resource constraints, Appl. Math. Comput., 371 (2020), 124956. https://doi.org/10.1016/j.amc.2019.124956 doi: 10.1016/j.amc.2019.124956

This article has been cited by:

Huaijun Deng, Linna Liu, Jianyin Fang, Boyang Qu, Quanzhen Huang, A novel improved whale optimization algorithm for optimization problems with multi-strategy and hybrid algorithm, 2023, 205, 03784754, 794, 10.1016/j.matcom.2022.10.023

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

Mathematical Biosciences and Engineering

4.4

Metrics

Article views(3825) PDF downloads(151) Cited by(7)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(3) / Tables(1)

Mathematical Biosciences and Engineering

Stability analysis and optimal control of COVID-19 with quarantine and media awareness

Related Papers:

Abstract

1. Introduction

2. Chemical reaction optimization

2.1. The basic concept of the CRO

2.2. Elementary reactions

3. Modified chemical reaction optimization (MCRO)

3.1. Initial population

3.2. A modified search equation

3.3. Boundary restrictions

3.4. Time complexity analysis

4. Results

4.1. The benchmark functions optimization by MCRO

4.2. Array pattern synthesis using MCRO

4.2.1. Synthesis of reconfigurable antenna-array

4.3. Synthesis of linear antenna

5. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Mathematical Biosciences and Engineering

Stability analysis and optimal control of COVID-19 with quarantine and media awareness

Related Papers:

Abstract

1. Introduction

2. Chemical reaction optimization

2.1. The basic concept of the CRO

2.2. Elementary reactions

3. Modified chemical reaction optimization (MCRO)

3.1. Initial population

3.2. A modified search equation

3.3. Boundary restrictions

3.4. Time complexity analysis

4. Results

4.1. The benchmark functions optimization by MCRO

4.2. Array pattern synthesis using MCRO

4.2.1. Synthesis of reconfigurable antenna-array

4.3. Synthesis of linear antenna

5. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog