Development of foreign direct investment in the context of government aid: A case study of the Czech Republic since 1998

Sylvie Kotíková; Roman Vavrek; Sylvie Kotíková; Roman Vavrek

doi:10.3934/NAR.2025009

National Accounting Review

2025, Volume 7, Issue 2: 196-220. doi: 10.3934/NAR.2025009

Previous Article Next Article

Research article

Development of foreign direct investment in the context of government aid: A case study of the Czech Republic since 1998

Sylvie Kotíková ^{1
,
,},
Roman Vavrek ²

1.
Department of Economics and Management, Faculty of Social and Economic Studies, University of Jan Evangelista Purkyně, Moskevská Street 54, Ústí nad Labem 400 96, Czech Republic
2.
Department of Management, Faculty of Economics, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava - Poruba 708 00, Czech Republic

Received: 31 October 2024 Revised: 23 March 2025 Accepted: 10 April 2025 Published: 14 April 2025
JEL Codes: E32, F21, J24, L20, O30, O40, R11

Investment incentives are policy instruments used by governments to stimulate economic growth, attract foreign capital, and develop business activities. This form of government aid plays an important role in decision-making processes regarding localization and host market selection. According to the core-periphery theory, investment incentives aim to reduce regional disparities and support the economic growth and development of structurally disadvantaged regions. Within the Czech Republic, 1070 investment incentives were granted, 599 of which were aimed at foreign investors. The objective of our research was to identify and describe trends and contexts in the development of foreign direct investment (which received an investment incentive) on the national and regional level in the Czech Republic from 1998 to the present (as of October 1, 2024). Overall, one main objective, two research hypotheses, and several partial objectives were defined. Already proven methods in this field, i.e., correlation analyses, the Kruskal-Wallis test, and Moran's coefficient, were used for the fulfilment of the set main objective and partial objectives. The number of foreign direct investments in the Czech Republic has already passed its golden era and similar patterns can be observed for follow-up reinvestments. It was identified that within the research set of 599 granted investment incentives, 113 cases represent reinvestments of original foreign investments. This fact signals that investors are staying in the region long term and shows growing satisfaction with the Czech business environment. However, the volume of these reinvestments has been continuously decreasing since 2014. From the spatial distribution point of view, the Ústí nad Labem Region (CZ042), the Moravian-Silesian Region (CZ072), and the Central Bohemian Region (CZ020) are clearly the most attractive.

Keywords:

Citation: Sylvie Kotíková, Roman Vavrek. Development of foreign direct investment in the context of government aid: A case study of the Czech Republic since 1998[J]. National Accounting Review, 2025, 7(2): 196-220. doi: 10.3934/NAR.2025009

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]	Adámek E, Rybková L (2015) Influence of investment incentives on development of regional unemployment in the Czech Republic. Eur J Bus Sci Technol 1: 5–14. https://doi.org/10.11118/ejobsat.v1i1.33 doi: 10.11118/ejobsat.v1i1.33
[2]	Amin A, Thrift N (1995) Globalization, Institutional Thickness and the Local Economy, In: Healey, P., Camerson, S., Davoudi, S., et al., Managing Cities: The New Urban Kontext, Chichester: Wiley, 322.
[3]	Bano S, Tabbada J (2015) Foreign Direct Investment Outflows: Asian Developing Countries. J Econ Integr 30: 359–398. https://doi.org/10.11130/jei.2015.30.2.359 doi: 10.11130/jei.2015.30.2.359
[4]	Barasa L, Knoben J, Vermeulen P, et al. (2017) Institutions, Resources and Innovation in East Africa: A Firm Level Approach. Res Policy 46: 280–291. https://doi.org/10.1016/j.respol.2016.11.008 doi: 10.1016/j.respol.2016.11.008
[5]	Belardi J (2022) Spatial distribution of granted investment incentives in the Czech Republic. Acta Acad Karvin 22: 5–15. https://doi.org/10.25142/aak.2022.011 doi: 10.25142/aak.2022.011
[6]	Bellak C, Leibrech M, Romisch R (2007) On the appropriate measure of tax burden on foreign direct investment to the CEECs. Appl Econ Lett 17: 603–606. https://doi.org/10.1080/13504850500474202 doi: 10.1080/13504850500474202
[7]	Bildirici M, Gokmenoglu SM (2020) The impact of terrorism and FDI on environmental pollution: Evidence from Afghanistan, Iraq, Nigeria, Pakistan, Philippines, Syria, Somalia, Thailand and Yemen. Environ Impact Assess Rev 81: 106340. https://doi.org/10.1016/j.eiar.2019.106340 doi: 10.1016/j.eiar.2019.106340
[8]	Blažek J, Květoň V (2023) Towards an Integrated Framework of Agency in Regional Development: The Case of Old Industrial Regions. Reg Stud 57: 1482–1497. https://doi.org/10.1080/00343404.2022.2054976 doi: 10.1080/00343404.2022.2054976
[9]	Blažek J, Uhlíř D (2011) Teorie regionálního rozvoje: nástin, kritika, implikace. Edited by Táborský, J. Vydání druhé, Přepracované a rozsířené. Prague, Czech Republic: Karolinum.
[10]	Bogataj D, Bogataj M, Drobne S (2019) Interactions between Flows of Human Resources in Fluctional Regions and Flows of Investors in Dynamic Processes of Global Supply Chains. Int J Prod Econ 209: 2015–2225. https://doi.org/10.1016/j.ijpe.2017.10.018 doi: 10.1016/j.ijpe.2017.10.018
[11]	Boudeville JR (1974) Problems of Regional Economic Planning, Reprint. Edinburgh: Edinburgh University Press.
[12]	Cao P, Li XX, Cheng Y, et al. (2022) Temporal-Spatial Evolution and Driving Factors of Global Carbon Emission Efficiency. Int J Environ Res Public Health 19: 14849. https://doi.org/10.3390/ijerph192214849 doi: 10.3390/ijerph192214849
[13]	Capik P (2021) Intensity of foreign investment promotion and development levels in Central-Eastern European regions. J Place Manag Dev 14: 114–133. https://doi.org/10.1108/JPMD-07-2019-0070 doi: 10.1108/JPMD-07-2019-0070
[14]	Cedidlova M (2013) The effectiveness of investment incentives in certain foreign companies operating in the Czech Republic. J Competitiveness 5: 108–120. https://doi.org/10.7441/joc.2013.01.08 doi: 10.7441/joc.2013.01.08
[15]	Chrobocińska K (2021) Comparative Analysis of Regional Competitiveness in Poland from 2010–2019 in the Context of the Concept of Sustainable Development. Sustainability 13. https://doi.org/10.3390/SU13063202 doi: 10.3390/SU13063202
[16]	Cortright J (2001) New Growth Theory, Technology and Learning: A Practitioner's Guide. Available from: https://e-tcs.org/wp-content/uploads/2012/10/Cortright-nueva_teoria_del_crecimiento.pdf.
[17]	Costaiche GM, Niculae I (2016) The Evolution and Impact of Foreign Direct Investments in Romania. Sci Pap-Ser Manag Econ Eng Agric Rural Dev 16: 111–116.
[18]	CzechInvest Investment Incentives. CzechInvest, 2024. Available from https://www.czechinvest.org/en/For-Investors/Investment-Incentives?force.
[19]	Damborský M (2023) Regional unemployment and investment incentives in Czechia. Geogr J 75: 107–124. https://doi.org/10.31577/geogrcas.2023.75.2.06 doi: 10.31577/geogrcas.2023.75.2.06
[20]	Damborský M, Šhovíková G, Rajtr R (2012) Regional Localization of the Automotive Industry in the Czech Republic. Acta Oecon Prag 20: 21–39. https://doi.org/10.18267/J.AOP.361 doi: 10.18267/J.AOP.361
[21]	Derado D, Horvatin D (2023) FDI and Economic Growth Perspective of Southeast European Countries. Zagreb Int Rev Econ Bus 26: 119–148. https://doi.org/10.2478/zireb-2023-0006 doi: 10.2478/zireb-2023-0006
[22]	Dicken P, Coe N, Hess M (2008) Global production network: realizing the potential. J Econ Geogr 8: 271–296. https://doi.org/10.1093/jeg/lbn002 doi: 10.1093/jeg/lbn002
[23]	Driffield N, Hughes D (2003) Foreign and Domestic Investment: Regional Development or Crowding Out? Reg Stud 37: 277–288. https://doi.org/10.1080/0034340032000065433 doi: 10.1080/0034340032000065433
[24]	Dunning JH (1981) International Production and the Multinational Entreprise. J Int Bus Stud 1981: 171–173.
[25]	Dunning JH, Lundan SM (2008) Multinational Enterprises and the Global Economy, 2nd ed. Northampton, MA: Edward Elgar.
[26]	Evan T, Bolotov I (2021) Can Investment Incentives Cause Unemployment? An Empirical Analysis of The Relationship Between FDI and Employment Based on The OLI Framework. Cent Eur Bus Rev 11: 1–16. https://doi.org/10.18267/J.CEBR.291 doi: 10.18267/J.CEBR.291
[27]	Fagerberg AJ (1987) A technology Gap Approach to Growth Rates Bigger. Res Policy 16: 87–99. https://doi.org/10.1016/0048-7333(87)90025-4 doi: 10.1016/0048-7333(87)90025-4
[28]	Fagerberg AJ, Lundvall BA, Srholec M (2018) Global Value Chains, National Innovation Systems and Economic Development. Eur J Dev Res 30: 533–556. https://doi.org/10.1057/s41287-018-0147-2 doi: 10.1057/s41287-018-0147-2
[29]	Fitala P, Kozakova J, Skypalova R, et al. (2023) Strategic Management: Factors influencing (De)centralization in Multinational Corporations operating in Slovakia. Mark Manag Innov 14: 223–238. https://doi.org/10.21272/mmi.2023.4-17 doi: 10.21272/mmi.2023.4-17
[30]	Gambus I, Almeida F (2018) Three Decades after James Street's. "The Institutionalist Theory of Economic Development": What Does Institutional Approach to Economic Development Mean Today? J Econ Issues 52: 455–463. https://doi.org/10.1080/00213624.2018.1469902 doi: 10.1080/00213624.2018.1469902
[31]	Gorynia M, Nowak J, Trpczyński P, et al. (2016) Determinants of FDI establishment mode choice of Polish firms: The OLI paradigm perspective. Argum Oecon 2: 67–92. https://doi.org/10.15611/aoe.2016.2.03 doi: 10.15611/aoe.2016.2.03
[32]	Grossman GM, Helpman E (1994) Endogenous Innovation in the Theory of Growth. J Econ Perspect 8: 23–44. https://doi.org/10.1257/jep.8.1.23 doi: 10.1257/jep.8.1.23
[33]	Ho LT, Gan C (2021) Foreign Direct Investment and World Pandemic Uncertainty Index: Do Health Pandemics Matter? J Risk Financ Manag 14. https://doi.org/10.3390/jrfm14030107 doi: 10.3390/jrfm14030107
[34]	Kotíková S (2019) Potential of the Czech business environment assumes the effects of foreign direct investment. E M Ekon Manag 22: 18–35. https://doi.org/10.15240/tul/001/2019-4-002 doi: 10.15240/tul/001/2019-4-002
[35]	Kotíková S (2023) Spillover effects: A challenging public interest to measure. Natl Account Rev 5: 373–404. https://doi.org/10.3934/NAR.2023022 doi: 10.3934/NAR.2023022
[36]	Krugman P (2011) The New Economic Geography, Now Middle-Aged. Reg Stud 45: 1–7. https://doi.org/10.1080/00343404.2011.537127 doi: 10.1080/00343404.2011.537127
[37]	Kuznetsov A (2020) Foreign Direct Investment, In: Kobayashi, A., International Encyclopedia of Human Geography (Second Edition), 219–227. https://doi.org/10.1016/B978-0-08-102295-5.10070-8
[38]	Lee HH, Hayakawa K, Park CY (2023) The Effect of COVID-19 on Foreign Direct Investment. Global Econ Rev 52: 1–17. https://doi.org/10.1080/1226508X.2023.2171459 doi: 10.1080/1226508X.2023.2171459
[39]	Lenaerts K, Merlevede B (2012) Horizontal or Backward? FDI Spillovers and Industry Aggregation. Available from: https://www.etsg.org/ETSG2011/Papers/Lenaerts.pdf.
[40]	Markowska M, Hlavacek P, Strahl D (2022) Knowledge-Intensive Business Services Employment Structure and Economic Development in EU Regions. Comp Econ Res 25: 109–133. https://doi.org/10.18778/1508-2008.25.32 doi: 10.18778/1508-2008.25.32
[41]	Mishrif A, Khan A (2024) Do Free Trade Agreements Facilitate FDI Spillover Effects on Domestic Firms? Empirical Evidence from Oman. Economies 12: 141. https://doi.org/10.3390/economies12060141 doi: 10.3390/economies12060141
[42]	Moosa, Imad A (2002) Foreign Direct Investment, In: International Financial Operations. Finance and Capital Markets Series. Palgrave Macmillan, London: Palgrave Macmillan, 318–340. https://doi.org/10.1057/9781403946034_12
[43]	Narula R (2018) An extended dual economy model. Int J Emerg Mark 13: 586–602. https://doi.org/10.1108/IJoEM-10-2017-0375 doi: 10.1108/IJoEM-10-2017-0375
[44]	North DC, Weingast BR (2000) Introduction: Institutional Analysis and Economic History. J Econ Hist 60: 414–417. https://doi.org/10.1017/S0022050700025158 doi: 10.1017/S0022050700025158
[45]	Pan WH, Ngo XT (2016) Endogenous Growth Theory and Regional Performance: The Moderating Effects of Special Economic Zones. Communist Post-Communist Stud 49: 113–122. https://doi.org/10.1016/j.postcomstud.2016.04.005 doi: 10.1016/j.postcomstud.2016.04.005
[46]	Posner MV (1961) International Trade and Technical Change. Oxf Econ Pap 13: 323–341. https://doi.org/10.1093/oxfordjournals.oep.a040877 doi: 10.1093/oxfordjournals.oep.a040877
[47]	Ross AG, Fleming N (2022) The impact of Chinese foreign direct investment on host country economic growth. Local Econ 37. https://doi.org/10.1177/02690942231162850 doi: 10.1177/02690942231162850
[48]	Sako M, Zylberberg E (2019) Supplier Strategy in Global Value Chains: Shaping Governance and Profiting from Upgrading. Socio-Econ Rev 17: 687–707. https://doi.org/10.1093/ser/mwx049 doi: 10.1093/ser/mwx049
[49]	Sattar A, Hassan A, Hussain MN, et al. (2022) Impact of foreign direct investment on socio-economic development in belt and road countries. Cogent Econ Financ 10. https://doi.org/10.1080/23322039.2022.2143772 doi: 10.1080/23322039.2022.2143772
[50]	Šimanová J, Trešl F (2011) Vývoj Průmyslové Koncentrace a Specializace v Regionech NUTS3 České Republiky v Kontextu Dynamizace Regionální Komparativní Výhody. E M Ekon Manag 14: 38–52.
[51]	Sun ZQ, Sun T (2020) The Impact of Multi-Dimensional Urbanization on China's Carbon Emissions Based on the Spatial Spillover Effect. Pol J Environ Stud 29: 3317–3327. https://doi.org/10.15244/pjoes/114508 doi: 10.15244/pjoes/114508
[52]	Szent-Ivanyi B, Vigvári G (2012) Spillovers from Foreign Direct Investment in Central and Eastern Europe. Soc Econ 34: 51–72. https://doi.org/10.1556/SocEc.34.2012.1.5 doi: 10.1556/SocEc.34.2012.1.5
[53]	Tang DL, Xu H, Yang YH (2018) Mutual Influence of Energy Consumption and Foreign Direct Investment on Haze Pollution in China: a Spatial Econometric Approach. Pol J Environ Stud 27: 1743–1752. https://doi.org/10.15244/pjoes/78046 doi: 10.15244/pjoes/78046
[54]	Torre A, Wallet F (2016) Policies of Regional Development and Planning, In: Regional Development in Rural Areas, Springer Briefs in Regional Science. Springer, Cham, 35–49. https://doi.org/10.1007/978-3-319-02372-4_4
[55]	Tsinaridze R, Makharadze N (2023) The Role of Foreign Direct Investment in the National Economy and Correlation Analysis. Access-Sci Business Innov Digit Econ 4: 46–59. https://doi.org/10.46656/access.2023.4.1(4) doi: 10.46656/access.2023.4.1(4)
[56]	Vernon R (1966) International Investment and International Trade in the Produkt Cycle. Q J Econ 80: 190–207. https://doi.org/10.2307/1880689 doi: 10.2307/1880689
[57]	Verspanger B (1991) A New Empirical Approach to Catching up or Falling Behind. Structural Change and Economic Dynamics 2: 359–380. https://doi.org/10.1016/S0954-349X(05)80008-6 doi: 10.1016/S0954-349X(05)80008-6
[58]	Vu TH, Trinh NTT (2023) The impact of FDI on economic growth in developing countries: the role of FDI inflow and trade openness. Econ Finance Lett 10: 216–229. https://doi.org/10.18488/29.v10i3.3519 doi: 10.18488/29.v10i3.3519
[59]	Xu J, Lv T, Hou XR, Deng X, Li N, Liu F (2022) Spatiotemporal characteristics and influencing factors of renewable energy production in China: A spatial econometric analysis. Energy Econ 116: 106399. https://doi.org/10.1016/j.eneco.2022.106399 doi: 10.1016/j.eneco.2022.106399
[60]	Zanoli SM, Pepe C (2023) Design and Implementation of a Fuzzy Classifier for FDI Applied to Industrial Machinery. Sensors 23: 6954. https://doi.org/10.3390/s23156954 doi: 10.3390/s23156954

NAR-07-02-009-S001.pdf

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:*
© 2025 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)