Next Article in Journal
Iterative Methods for Computing the Resolvent of Composed Operators in Hilbert Spaces
Previous Article in Journal
Search Acceleration of Evolutionary Multi-Objective Optimization Using an Estimated Convergence Point
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Equilibrium and Stability of Entropy Operator Model for Migratory Interaction of Regional Systems

by
Yuri S. Popkov
1,2,3,4,* and
Leo van Wissen
4,5,6
1
Institute for Systems Analysis of FRC CSC of RAS, Moscow 119333, Russia
2
Institute of Control Sciences of Russian Academy of Sciences, Moscow 117997, Russia
3
National Research, University Higher School of Economics, Moscow 125319, Russia
4
Department of Engineering, Haifa University, Karmiel 2161002, Israel
5
Netherlands Interdisciplinary Demographic Institute, P.O. Box 11650, 2502 AR Den Haag, The Netherlands
6
Population Research Centre, University of Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlands
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(2), 130; https://doi.org/10.3390/math7020130
Submission received: 26 December 2018 / Revised: 21 January 2019 / Accepted: 26 January 2019 / Published: 29 January 2019
(This article belongs to the Section Mathematics and Computer Science)

Abstract

:
In this paper, we develop a dynamic model for migratory interaction of regional systems that is based on an entropy operator. Next, we study the properties of this operator and establish the existence of a unique singular point in the dynamic entropy model. Here, we use monotonicity property of entropy operator on corresponding vector interval. We study Lyapunov stability of a dynamic system with entropy operator. Stability conditions have been obtained in terms of eigenvalues of linearized system’s matrix. Finally, we give an illustrative example for migratory interaction of regional systems.

1. Introduction

The processes of migration and biological reproduction represent major factors determining the spatial distribution of populations. These factors are realized in considerably different time scales. In accordance with standard assumptions, the relaxation time of migratory processes is 2–3 years while the relaxation time of biological reproduction reaches 10–13 years, see [1,2,3]. An important characteristic of population is its size, more specifically, its distribution in some regionalized space. Migration considerably affects the spatiotemporal evolution of population distribution. For the time intervals comparable with the relaxation times of migratory processes, the influence of biological reproduction factors becomes rather small and therefore can be neglected.
Mathematical modeling and analysis of interregional migration was considered in numerous publications. First, it seems appropriate to mention the monographs [4,5], which were dedicated to a wide range of interregional migration problems, including mathematical modeling of migration flows. Note that the problem of migration touches upon many aspects of socioeconomic, psychological and political status of the space of migratory movements. So of crucial role is the structural analysis of inter- and intraregional migration flows [6] and motivations that generate them [7,8]. The results of structural and motivational analysis of migratory processes are used for computer simulation. There exist three directions of research in this field, each relying on some system of hypotheses. One of the directions involves the stochastic hypothesis about the origin of migratory motivations [9], which is simulated using agent technologies [10,11]. This direction is adjoined by investigations based on the thermodynamic model of migration flows [12,13]. The other two directions are connected with using balanced [4,5] and Markov [2] models. Of course, the short list above does not exhaust the whole variety of migration studies, merely outlining some topics of research.
This paper considers the interaction of two regional systems as follows. The first system has traditional migratory processes while the second is a source of immigration flows to the former. We suggest a mathematical model of such interaction in form of a dynamic system with an entropy operator. This model proceeds from the stochastic hypothesis about the origin of migration flows, which is realized by conditional maximization of the generalized information entropy. Next, we develop a qualitative analysis method for the dynamic system, including the existence of a unique equilibrium distribution of regional population sizes and its stability.
Finally, the method suggested in this paper is illustrated by modeling the migratory interactions between three EU countries (Germany, France, Italy—the system GFI ) and two countries as sources of immigration (Syria and Libya—the system SL ). We construct the set of equilibrium distributions for regional population size that is parameterized by the following characteristics: regional mobility, the prior probabilities of migratory movements, and the normalized specific regional cost of immigration (the system GFI ); the prior probabilities of immigratory movements (the system SL ).

2. Structure of Dynamic Entropy Model of Migration

Consider a spatial system partitioned into N regions that exchange population. At a time t, the state of this system is characterized by the population distribution among all regions, i.e., by the vector K ( t ) = { K 1 ( t ) , , K N ( t ) } . The population distribution varies with the course of time as the result of interregional migration within the system and also external immigration from M source regions outside the system.
We will employ the stochastic model of migration flows [14]. In accordance with this model, interregional flows are formed by a random independent distribution of potential migrants between system regions with some prior probabilities. A similar mechanism describes immigration flows from external regions to system regions. These mechanisms (the distribution of migration flows) are defined in probabilistic terms using appropriate entropy functions or functionals. Migratory interactions between system regions take place under constraints on the size of mobile population (the people inclined to migration). By standard assumption, the size of mobile population in a region is proportional to the regional population size, and the coefficient of proportionality has regional index. The potential amount of external immigration is often difficult or even impossible to calculate. However, sometimes there exists expert appraisals for immigration flows from external regions to system regions. Such appraisals often seem to be strongly overrated. Therefore, the constraints on external immigration have to be associated with regional absorption cost of immigrants. The entropy characteristics of migration flows as a function of flow distribution possess the so-called sharp maximum. This justifies the hypothesis about the existence of a realizable distribution of migration and immigration flows corresponding to the constrained maximum of the entropy [14].
Introduce the following notations for transforming the described phenomenology into a mathematical model (see [14]):
  • x i n ( t ) , ( i , n ) = 1 , N ¯ , as the migration flows ( i j ) within the system;
  • a i n ( t ) , ( i , n ) = 1 , N ¯ , as the prior probabilities of individual migration between regions ( i , j ) of the system;
  • α n , n = 1 , N ¯ , as the shares of mobile population in system regions;
  • y k n , k = 1 , M ¯ , n = 1 , N ¯ , as the immigration flows ( k j ) from external regions to system regions;
  • b k n , k = 1 , M ¯ , n = 1 , N ¯ , as the prior probabilities of individual immigration from external region k to system region j;
  • c k n , k = 1 , M ¯ , n = 1 , N ¯ , as the normalized (Here, normalization means that the values of specific cost belong to the range [0, 1].) specific generalized cost of immigration to system regions;
  • c k n ( s ) , k = 1 , M ¯ , n = 1 , N ¯ , s = 1 , r ¯ , as the normalized specific cost of immigration to system regions by the types of resources ( r ) ;
  • T ( K 1 , , K N ) as the normalized (Here, the normalized supply functions are defined by
    T ( K 1 , , K N ) = 1 c T ˜ ( K 1 , , K N ) , c = max c k n ,
    where T ˜ ( K 1 , , K N ) denotes a given supply function.) supply function of generalized financial resources (cost) for immigration;
  • T ( s ) ( K 1 , , K N ) as the normalized supply functions by s = 1 , r ¯ types of financial resources for immigration.
We will model the migration and immigration flows using the stochastic migration mechanisms characterized by
  • the generalized Boltzmann information entropy
    H M ( X ) = ( i , n ) = 1 ; i n N x i n ln x i n a i n
    for interregional migration within the system and
  • by the generalized Boltzmann information entropy
    H I ( X ) = k = 1 M n = 1 N y k n ln y k n b k n
    for the immigration from external regions to the system.
Migratory processes always have some constraints. For interregional migration, they are associated with the size of mobile population in regions, i.e.,
X D ( K ) = { X : i = 1 N x n i = α n K n ( t ) , n = 1 , N ¯ } .
In the case of external immigration, such constraints are imposed on the system resources required for external immigration. In accordance with the classification above, we will consider the following models of immigration cost:
  • the generalized cost of the form
    Y D G ( K ) = { Y : k = 1 M n = 1 N c k n y k n = T ( K ( t ) ) } , c k n [ 0 , 1 ] ; T ( K ( t ) ) = 1 c T ˜ ( K ( t ) ) , c = max k n c k n ;
  • the cost by the types of resources of the form
    Y D D ( K ) = { k = 1 M n = 1 N c k n ( s ) y k n = T ˜ ( s ) ( K ( t ) ) } ; c k n ( s ) [ 0 , 1 ] ; T ˜ ( s ) ( K ( t ) ) = 1 c ( s ) T ( s ) ( K ( t ) ) , c ( s ) = max k n c k n ( s ) , s = 1 , r ¯ .
The functions T ( K ( t ) ) and T ( 1 ) ( K ( t ) ) , , T ( r ) ( K ( t ) ) depend on the population distribution in the system, and these relationships are linear [14]:
T ( K ( t ) ) = n = 1 N γ n K n ( t ) , T ( s ) ( K ( t ) ) = n = 1 N β s n K n ( t ) , s = 1 , r ¯ .
Therefore, the entropy models of population distribution have the following form:
  • for the interregional migration flows within the system,
    x i n ( t ) = arg max H M ( X ) , X D ( K ) , ( i , n ) = 1 , N ¯ ;
  • for the immigration flows from external regions to the system (with the generalized cost),
    y k n ( t ) = arg max H I ( Y ) , Y D G ( K ) , k = 1 , M ¯ , n = 1 , N ¯ ;
  • for the immigration flows from external regions to the system (with the cost by the types of resources),
    y k n ( t ) = arg max H I ( Y ) , Y D D ( K ) , k = 1 , M ¯ , n = 1 , N ¯ .
These models represent the corresponding entropy operators that map the space R + n of all vectors K ( t ) into the space R + N 2 of all matrices X, or into the space R + N M of all matrices Y.

3. Analysis of Entropy Operators

3.1. Entropy Operator K X

Consider the constrained maximization problem for entropy (1) and (4) and construct the corresponding Lagrange function
L M ( X , λ ) = H M ( X ) + n = 1 N λ n ( α n K n ( t ) i = 1 N x n i ) ,
where λ = { λ 1 , , λ N } are the Lagrange multipliers. The optimality conditions have the form
ln x n i a n i + 1 λ n = 0 , ( n , i ) = 1 , N ¯ ; i = 1 N x n i α n K n ( t ) = 0 , n = 1 , N ¯ .
The solution of these equations, i.e., the migration flows from region n to region i ( i n ) are given by
x n i = a n i l = 1 N a n l α n K n ( t ) , ( n , i ) = 1 , N ¯ ,
and form the matrix X . In this case, the mapping ( K R + N X R + N 2 ) (the entropy operator) can be written in the analytic form (12).
The transposed matrix X consists of the elements
x i n = a i n l = 1 N a l n α i K i ( t ) , ( n , i ) = 1 , N ¯ .
They represent the immigration flows from region i to region n ( i n ).
Hence, for region n the net migration flow is
S n ( t ) = i = 1 N a i n α i K i ( t ) l = 1 N a l n α n K n ( t ) , n = 1 , N ¯ .

3.2. Entropy Operator K Y

Now, consider the constrained maximization problem for entropy (2) and (4). The corresponding Lagrange function is defined by
L I ( Y , λ ) = H I ( Y ) + λ T ( K ( t ) ) k = 1 M n = 1 N c k n y k n ,
where λ is the Lagrange multiplier and the normalized specific generalized cost has the form
c k 0 , n 0 = c = 1 , 0 < c k n < 1 .
Here, the optimality conditions are
ln y k n b k n + 1 λ c k n = 0 , k = 1 , M ¯ , n = 1 , N ¯ , k = 1 M n = 1 N c k n y k n + T ( K ( t ) ) = 0 .
These equations yield the entropy-optimal immigration flows from external regions to the system:
y k n = b ˜ k n z c k n , k = 1 , M ¯ , n = 1 , N ¯ , z = exp ( λ ) 0 , b ˜ k n = e 1 b k n .
The exponential Lagrange multiplier z satisfies the equation
k = 1 M n = 1 N c k n b ˜ k n z c k n T ( K ( t ) ) = 0 .
In this case, the mapping K Y (the entropy operator) is described by equalities (18) and Equation (19).
The total immigration flow from external regions to the system constitutes
I n ( t ) = M y k n = k = 1 M b ˜ k n z c k n , n = 1 , N ¯ .

4. Analysis of Migratory Interaction Model of Regional Systems

Under the above assumptions, the system state evolves with the internal net migration flows and also the net immigration flows from external regions, i.e.,
d K n ( t ) d t = i = 1 N [ x i n ( t ) x n i ( t ) ] + k = 1 M y k n , n = 1 , N ¯ .
Using the formulas of these flows, we obtain the following equations of the dynamic migratory interaction model of regional systems:
d K n ( t ) d t = α n K n ( t ) + i = 1 N a i n l = 1 N a l n α i K i ( t ) + k = 1 M b ˜ k n z c k n , n = 1 , N ¯ , k = 1 M n = 1 N c k n b ˜ k n z c k n T ( K ( t ) ) = 0 , z 0 .

4.1. Existence of Unique Singular Point

The existence problem of unique equilibrium in this model is associated with analysis of the equations
( 1 α n ) K n + i = 1 , i n N a i n α i K i + k = 1 M b k n z c k n = K n , n = 1 , N ¯ , n = 1 N γ ˜ n K n ( t ) + z k = 1 M n = 1 N c k n b ˜ k n z c k n = z .
Introduce the vector of variables
u = { u 1 , , u N , u ( N + 1 ) } , u n = K n , n = 1 , N ¯ ; u ( N + 1 ) = z .
Then Equation (23) can be written as
A ( u ) = u , u 0 ,
with the vector function
A ( u ) = ( 1 α 1 ) u 1 + i = 1 , i 1 N a i 1 α i u i + k = 1 M b ˜ k 1 u ( N + 1 ) c k 1 ( 1 α N ) u N + i = 1 , i N N a i N α i u i + k = 1 M b ˜ k N u ( N + 1 ) c k n n = 1 N γ ˜ n u n + u N + 1 k = 1 , k k 0 M n = 1 , n n 0 N c ˜ k n b ˜ k n u ( N + 1 ) c k n .
The Jacobian of this function has the form
J ( u ) = J 11 J 12 J 21 J 22 ,
where
J 11 = 1 α 1 a 21 α 2 a N 1 α N a 12 α 1 1 α 2 a N 2 α N a 1 N α 1 a 2 N α 2 1 α N > 0 , J 12 = k = 1 M b ˜ k 1 c k 1 u N + 1 ( c k 1 1 ) k = 1 M b ˜ k 2 c k 2 u N + 1 ( c k 2 1 ) k = 1 M b ˜ k N c k N u N + 1 ( c k N 1 ) > 0 , J 21 = γ 1 γ 2 γ N > 0 , J 22 = 1 k = 1 M n = 1 N c k n 2 b ˜ k n u N + 1 ( c k n 1 ) .
Theorem 1.
Let J 22 > 0 . Then the function A ( u ) is strictly monotonic over the space R + ( N + 1 ) .
Proof. 
Under the hypothesis of this theorem, all elements of the Jacobian are positive, which proves the desired result. ☐
Assume we know the population distribution among system regions at the initial time t 0 , i.e., { K 1 ( t 0 ) = u 1 0 , , K N ( t 0 ) = u N 0 } . Denote by u N + 1 0 > 0 the solution of the last equation in (25) and (26). Define the vector u 0 = { u 1 0 , , u N 0 , u N + 1 0 } . Suppose that, by an expert appraisal of immigration from external regions, the population size in system regions may increase by κ > 1 times. This means that { K 1 ( t ) < κ u 1 0 = u ˜ 1 , , K N ( t ) < κ u N 0 = u ˜ N + 1 } for t > t 0 . Designate as u ˜ N + 1 > u N + 1 0 > 0 the solution of the last equation in (25) and (26) with the described upper estimates for regional population sizes. Construct the vector u ˜ = { u ˜ 1 , , u ˜ N , u ˜ N + 1 } and define the vector segment U = { u 0 , u ˜ } .
Theorem 2.
For Equation (26), let
A ( u 0 ) > u 0 , A ( u ˜ ) < u ˜ .
Then there exists an leats one solution u of Equation (26) on the vector segment U .
The proof of this result follows from ([15], Theorem 3.1 in p. 48). Moreover, by this theorem the iterative process
A ( u s ) = u ( s + 1 ) , u ( 0 ) U ,
converges to the point u .

4.2. Stability of Singular Point

For the stability analysis of a singular point, represent the system of Equation (22) in the form
d u d t = A ( u ) u .
Denote by y = u u the deviation from the singular point. Then the standard linearization procedure yields the following differential equations in the deviation from this singular point:
d y d t = J ˜ ( u ) y , y R ( N + 1 ) ,
where
J ˜ ( u ) = J ( u ) E .
The system of differential Equation (32) is linear, and there exist many stability analysis methods for its equilibrium. We will described a method based on the localization procedure of the eigenvalues of the matrix J ˜ ( u ) using the Gershgorin circles [16]. In accordance with the Gershgorin circle theorem, on the complex plane any eigenvalue of the matrix J ˜ ( u ) belongs to at least one of the following circles with center c i and radius ϱ i :
c i = α i , i = 1 , N ¯ , k = 1 M n = 1 N c k n 2 b ˜ k n u N + 1 ( c k n 1 ) , i = N + 1 .
ϱ i = j i ( N + 1 ) | d i j | , i = 1 , ( N + 1 ) ¯ ,
where
d i j = a i j , ( i , j ) = 1 , N ¯ , j i , k = 1 M b ˜ k i c k i [ u N + 1 ( c k i 1 ) ] , i = 1 , N ¯ , j = N + 1 , γ j , i = N + 1 , j = 1 , N ¯ , k = 1 M n = 1 N c k n 2 b ˜ k n [ u N + 1 ( c k n 1 ) ] , ( i , j ) = N + 1 .
These equalities dictate that the centers are c i < 0 , i = 1 , ( N + 1 ) ¯ , while the radii are ϱ i > 0 , i = 1 , ( N + 1 ) ¯ .
Theorem 3.
Let
c i + ϱ i < 0 , i = 1 , ( N + 1 ) ¯ .
The the singular point u of system (22) is stable in the small.
Proof. 
Under the hypotheses of this theorem, all eigenvalues of the matrix J ˜ ( u ) have negative real parts. ☐

5. Model of Migratory Interaction: An Example

The appearance of territories with low economic status always causes the growth of immigration. The early 2000s were remarkable for the formation of several such territories in the Northern and Central Africa, the Near East, Afghanistan, etc. As a result, tens of millions of migrants moved to the EU as the level of life in these territories dropped below the subsistence minimum. The EU countries have to allocate considerable financial resources for their filtering and accommodation, which are often unacceptable.
As an illustrative example, consider the migratory interactions between three EU countries (the system GFI composed of regions with index n, namely, Germany (1), France (2), and Italy (3)) and two sources of immigration (the system SL composed of regions with index k, namely, Syria (1) and Lybia (2)). In terms of the mathematical model (22), we have (see [2,4])
d K n ( t ) d t = α n K n ( t ) + i = 1 , i n 3 a i n l = 1 , l n n a l n α i K i ( t ) + k = 1 2 b ˜ k n z c k n , n = 1 , 2 , 3 ; n = 1 3 γ n K n ( t ) k = 1 2 n = 1 3 c k n b ˜ k n z c k n = 0 .

5.1. Data Arrays for Numerical Study of Model (38)

This model contains four groups of parameters as follows:
  • α 1 , α 2 , and α 3 as the characteristics of population mobility in the system GFI ;
  • A = [ a i n , ( i , n ) = 1 , N ¯ , a n n = 0 ] as the matrix of the prior probabilities of individual migration in the system GFI ;
  • B = [ b k n , k = 1 , 2 ; n = 1 , 2 , 3 ] as the matrix of the prior probabilities of individual immigration from the system SL ;
  • C = [ c k n , k = 1 , 2 ; n = 1 , 2 , 3 ] as the matrix of normalized specific cost of the system GFI for maintaining the immigration from the system SL .
Recall that the elements of the first three groups of parameters have values strictly smaller than 1. The maximal element of the matrix C is 1 while the other elements are strictly smaller than 1. The coefficients are γ 1 = 2.00 , γ 2 = 2.53 , and γ 3 = 1.65 . (These values are defined by the immigration budgeting rules of the EU countries).
Introduce the four-dimensional attribute space P R + 4 of the status variables: ν (mobility), π (migration choice), μ (immigration choice), and T (cost scale). For each group, assign the following rating points for the parameters using the ten-point scale:
  • mobility ( m ) : m 1 = 7 , m 2 = 8 , m 3 = 5 ;
  • migration choice ( h ) : h 12 = 3 , h 13 = 3 , h 21 = 5 , h 23 = 2 , h 31 = 6 , h 32 = 4 ;
  • immigration choice ( q ) : q 11 = 6 , q 12 = 5 , q 13 = 4 , q 21 = 4 , q 22 = 6 , q 23 = 6 ;
  • immigration cost ( e ) : e 11 = 5 , e 12 = 6 , e 13 = 5 , e 21 = 5 , e 22 = 6 , e 23 = 5 .
These points are mostly determined by regional population sizes and the levels of economic development. Then the parameters in these groups can be written as follows:
  • the characteristics of population mobility in the system GFI as
    α n = ν m n , n = 1 , 2 , 3 ;
  • the prior probabilities of individual migration in the system GFI as
    a i n = π h i n , ( i , n ) = 1 , 2 , 3 ;
  • the prior probabilities of individual immigration from the system SL as
    b ˜ k n = μ q k n , k = 1 , 2 ; n = 1 , 2 , 3 ;
  • the normalized specific cost of the system GFI for maintaining the immigration from the system SL as
    c k n = T e k n , k = 1 , 2 ; n = 1 , 2 , 3 .
Define the admissible domains F of the status variables in the form
F = N P M T ,
ν N = [ 10 3 , 2 × 10 3 ] ; π P = [ 0.5 × 10 1 , 1.5 × 10 1 ] ;
μ M = [ 0.5 × 10 1 , 1.5 × 10 1 ] ; T T = [ 0.1 , 0.5 ] .
Exhaustive search over a spatial grid will be used below. So divide each interval into four subintervals. Then we have exhaustive search on 256 nodes of the grid. The discrete analogs of these intervals have the form
N ˜ = 10 ( 1 ) · [ 0.5 ; 1.0 ; 1.5 ; 2.0 ] , P ˜ = 10 1 · [ 0.75 ; 1.00 ; 1.25 ; 1.50 ] , M ˜ = 10 1 [ 0.75 ; 1.00 ; 1.25 ; 1.50 ] , T ˜ = [ 0.2 ; 0.3 ; 0.4 ; 0.5 ] .
The grid represents is the direct product of these intervals, i.e.,
F ˜ = N ˜ P ˜ M ˜ T ˜ .

5.2. Equilibrium Distributions of Regional Population Sizes in System GFI

For the equilibrium analysis of the system GFI , we have to solve the equations
ν m n K n + i = 1 , i n 3 π h i n l = 1 , l n n π h l n ν m i K i + k = 1 2 μ q k n z T e k n = 0 , n = 1 , 2 , 3 ; 2.0 K 1 + 2.53 K 2 + 1.65 K 3 k = 1 2 n = 1 3 T e k n μ q k n z T e k n = 0 .
The initial values of the variables are K 1 0 = 82 , K 2 0 = 67 , K 3 0 = 60 , and z 0 = 0.1 . Population size is measured in million people. Figure 1 shows the calculated equilibrium distributions ( K 1 , K 2 , K 3 ) of regional population sizes for the systems GFI and SL over the grid F ˜ (46). The cloud K contains 256 points corresponding to the nodes of grid (46). The center of this cloud matches the average distribution of regional population sizes
K ¯ = 1 256 i = 1 256 K i = { 82.7 ; 67.9 ; 60.3 } .
A useful characteristic of the cloud is its variance core R , i.e., the ball inside the cloud of radius
R = 1 256 i = 1 256 ( K K ¯ ) ( K K ¯ ) = 1.17 .

5.3. Stability of Equilibrium Distribution of Regional Population Sizes in System GFI

The equilibrium distributions are stable in the small if and only if the real parts of the eigenvalues of the matrix J ˜ (33) belong to intervals (37). On the strength of Theorem 3, we have the stability condition
max i c i i + ϱ i 0 ,
where ( ) denotes the values of variables corresponding to the equilibrium distribution, ϱ i is defined by equality (35), and
c i i = ν m i , i = 1 , 3 ¯ , [ T 2 μ ] k = 1 M n = 1 N e k n 2 q k n [ z ] ( T e k n 1 ) , i = 4 ;
d i j = π h i j , ( i , j ) = 1 , 3 ¯ , j i , μ ] k = 1 M q k i e k i [ z ] ( T e k i 1 ) , i = 1 , 3 ¯ , j = 4 , γ j , i = 4 , j = 1 , 3 ¯ , [ T 2 μ ] k = 1 M n = 1 N e k n 2 q k n [ z ] ( T e k n 1 ) , ( i , j ) = 4 .
Figure 1 illustrates three equilibrium distributions (marked by triangles) with the coordinates
C 1 = { 82.6228 ; 67.4256 ; 61.0597 } , C 2 = { 82.5310 ; 67.4051 ; 60.9496 } , C 3 = { 82.2856 ; 67.2362 ; 60.8639 } ,
which satisfies the stability conditions. These stable in the small equilibrium distributions of regional population sizes are associated with the parameters of model (38) combined in Table 1, Table 2, Table 3 and Table 4.
Therefore, using the known parameters of model (38), we may calculate its major dynamic characteristics, the stable equilibrium distributions of regional population sizes, which can be further used, e.g., for short-term forecasting.

6. Conclusions

In this paper, we have developed a dynamic entropy model of migratory interactions. Next, we have studied the entropy operators describing the distributions of migration and immigration flows as well as established the existence conditions of unique equilibrium distributions of regional population sizes and their stability conditions. The new model has been applied to analyze the equilibrium distributions of regional population sizes and their stability. The existence and stability conditions of these distribution both depend on the model parameters, which can be found using randomized machine learning procedures.

Author Contributions

The authors contributed equally to this work. All authors read and approved the final manuscript.

Funding

This work was partially supported by Russian Foundation for Basic Research (project no. 16-29-12878).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bilecen, B.; Van Mol, C. Introduction: international academic mobility and inequalities. J. Ethic Migr. Stud. 2017, 43, 1241–1255. [Google Scholar] [CrossRef] [Green Version]
  2. Weidlich, W. Sociodynamics; Harwood Academic Publishers: Amsterdam, The Netherlands, 2000; 456p. [Google Scholar]
  3. Zelinsky, W. The hypothesis of the mobility transition. Geogr. Rev. 1971, 46, 219–249. [Google Scholar] [CrossRef]
  4. Rogers, A.; Willekens, F.; Raymer, J. Modelling interregional migration flows: continuity and change. J. Math. Popul. Stud. 2001, 9, 231–263. [Google Scholar] [CrossRef]
  5. Rogers, A.; Little, J.; Raymer, J. The Indirect Estimation of Migration: Methods for Dealing with Irregular, Inadequate, and Missing Data; Springer: Berlin, Germany, 2010. [Google Scholar]
  6. Van der Knaap, G.A.; Steegers, W.F. Structural Analysis of Interregional and Intraregional Migration Patterns. In Demographic Research and Spatial Policy; Heide, H., Willekens, F., Eds.; Academic Press: New York, NY, USA, 1984. [Google Scholar]
  7. Black, R.; Xiang, B.; Caller, M.; Engberson, G.; Heering, L.; Markova, E. Migration and Development Causes and Consequences. In The Dynamics of International Migration and Settlement in Europe: A State of the Art; Penninx, R., Berger, M., Kraal, K., Eds.; Amsterdam University Press: Amsterdam, The Netherlands, 2006; pp. 41–63. [Google Scholar]
  8. Volpert, V.; Petrovskii, S.; Zincenko, A. Interaction of human migration and wealth distribution. Nonlinear Anal. 2017, 150, 408–423. [Google Scholar] [CrossRef]
  9. Pan, J.; Nagurney, A. Using Markov chains to model human migration in a network equilibrium framework. Math. Comput. Model. 1994, 19, 31–39. [Google Scholar] [CrossRef]
  10. Klabunde, A.; Zinn, S.; Willekens, F.; Leuchter, M. Multistable modelling extended by behavioural rules. An application to migration. Popul. Stud. 2017, 71, 61–67. [Google Scholar] [CrossRef] [PubMed]
  11. Klabunde, A.; Willekens, F. Decision-making in agent-based models of migration: State of the art and challenges. Eur. J. Popul. 2016, 32, 73–97. [Google Scholar] [CrossRef] [PubMed]
  12. Imel’baev, S.S.; Shmul’yan, B.L. Modeling of Stochastic Communication Systems. In Entropy Methods for Complex Systems Modeling; Wilson, A.G., Ed.; Nauka: Moscow, Russia, 1975; pp. 170–234. [Google Scholar]
  13. Wilson, A.G. Entropy in Urban and Regional Modelling; Pion Ltd.: Abingdon, UK, 2013. [Google Scholar]
  14. Popkov, Y.S. Mathematical Demoeconomy: Integrating Demographic and Economic Approaches; De Gruyter: Vienna, Austria, 2014. [Google Scholar]
  15. Krasnosel’skii, M.A.; Vainikko, G.M.; Zabreyko, R.P.; Ruticki, Y.B.; Stet’senko, V.V. Approximate Solutions of Operator Equations; Nauka: Moscow, Russia, 1969. (In Russian) [Google Scholar]
  16. Wilkinson, J.H. The Algebraic Eigenvalue Problem; Clarendon Press: Oxford, UK, 1965. [Google Scholar]
Figure 1. Equilibrium distributions.
Figure 1. Equilibrium distributions.
Mathematics 07 00130 g001
Table 1. Mobility coefficients.
Table 1. Mobility coefficients.
Notations α 1 α 2 α 3
C 1 0.01050.01200.0075
C 2 0.00700.00800.0050
C 3 0.01400.01600.0100
Table 2. Prior probabilities of regional mobility.
Table 2. Prior probabilities of regional mobility.
C 1 C 2 C 3
0.0000.3750.3750.0000.3000.3000.0000.4500.450
0.6250.0000.2500.5000.0000.2000.7500.0000.300
0.7500.5000.0000.6000.4000.0000.9000.6000.000
Table 3. Prior probabilities of immigration.
Table 3. Prior probabilities of immigration.
C 1 C 2 C 3
0.9000.7500.6000.7500.6250.5000.6000.5000.600
0.6000.9000.9000.5000.7500.7500.4000.6000.600
Table 4. Normalized specific cost.
Table 4. Normalized specific cost.
C 1 C 2 C 3
0.831.000.830.660.800.660.500.600.50
0.831.000.830.660.800.660.500.600.50

Share and Cite

MDPI and ACS Style

Popkov, Y.S.; van Wissen, L. Equilibrium and Stability of Entropy Operator Model for Migratory Interaction of Regional Systems. Mathematics 2019, 7, 130. https://doi.org/10.3390/math7020130

AMA Style

Popkov YS, van Wissen L. Equilibrium and Stability of Entropy Operator Model for Migratory Interaction of Regional Systems. Mathematics. 2019; 7(2):130. https://doi.org/10.3390/math7020130

Chicago/Turabian Style

Popkov, Yuri S., and Leo van Wissen. 2019. "Equilibrium and Stability of Entropy Operator Model for Migratory Interaction of Regional Systems" Mathematics 7, no. 2: 130. https://doi.org/10.3390/math7020130

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop