**Journal of Applied Mathematics and Physics**

Vol.07 No.02(2019), Article ID:90671,8 pages

10.4236/jamp.2019.72027

Modeling of Economic Processes, Instability and Chaos

Evgenii Perevoznikov, Elena Lomteva^{ }

A. F. Mozhaisky Military Space Academy, St. Petersburg, Russia

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: January 16, 2019; Accepted: February 19, 2019; Published: February 22, 2019

ABSTRACT

We consider a mathematical model describing the dynamics of the three components of the economic activity of enterprises―production, management and resource consumption. For the analysis of the dynamics, the methods and criteria formulated by the authors earlier were used; new results of a practical nature were obtained; in particular, the calculation of the spectrum of small perturbations of the system depending on the input resource and system parameters is given. It is shown that the calculated values of the criteria of instability and chaos are fully consistent with the calculation of the spectrum, indicated the stability intervals of the system, as well as areas of instability and chaos.

**Keywords:**

Mathematical Modeling of Economic Processes, Methods for Analyzing Dynamics, Instability, Chaos

1. Introduction

Methods of mathematical modeling are widely used to describe the economic processes [1] - [6] .

Modeling of economic processes considers: 1) method of studying phenomena and processes in the economic system or its elements by building and studying their models; 2) transfer of the results obtained during the study of models to the original; 3) using models to determine or refine the characteristics; 4) building newly constructed objects in the economic system.

Forms of modeling are different and depend on models and scope. The idea of modeling is based on any method of scientific research of economic processes. There are many models that display various aspects of economic processes.

Modeling serves as a mean of studying the economy and the phenomena occurring in it, substantiating decisions made, forecasting, planning, managing economic processes and objects. A model of an economic object is usually supported by real statistical data, and the results of calculations carried out within the framework of the constructed model make it possible to build forecasts and carry out objective evaluations.

For successful modeling, justified simplifications in the model compared to the original are important; an analysis of fairly simple and representative dynamic models is necessary [1] - [6] .

This article is a continuation of the work of the authors on the study of the dynamics of economic processes on the example of the “production―management―resource” model [7] .

The considered model PMR (production, management, resource) refers to the dynamic models of enterprises. Along with sufficient simplicity, it takes into account the main factors of production: the dynamics of production, management and the balance of external and system-produced resources. The model consists of three balanced, non-linear, differential equations corresponding to the three components of the required economic object subsystems―production, management and resource.

To analyze the dynamics of the system, assessing stability and possible instabilities, the methods and criteria given by the authors in [8] [9] [10] are used.

2. Initial Relations

Equations of the PMR model in a dimensionless form

$\begin{array}{l}{\partial}_{t}x=x\left(a-y+bz\right)\\ {\partial}_{t}y=y\left(x-c+dz\right)\\ {\partial}_{t}z=z\left(-e-dy-bx\right)-hx-ky+gxy+r\end{array}$ (1)

Here, $\frac{\left(x,y,z\right)}{{x}_{0}+{y}_{0}+{z}_{0}}$ the main dynamic values having the same resource

dimension (for example, monetary) are related to the sum of the initial values and corresponding to the dynamics of funds in the three subsystems of the object under study; x―production, y―control, z―resource. $t=\tau /{\tau}_{0}$ ―dimensionless time, ${\tau}_{0}$ ―a characteristic financial time scale (quarter, year).

In the first equation of the system (1) describing the dynamics of production, the components ax, bxz correspond to the profits from the production and input resource; the term xy is a decrease in management (for example, for training, professional development, etc.). In the second equation for controlling yx, dyz, the flow of funds from production and the input resource; cy―direct management costs.

The third equation is the balance of the internal and external r resource (r > 0―credit, r < 0―debt repayment).

The terms hx, ky are the direct loss of the salary resource in production and management; ez―resource costs for repayment of interest in lending or receiving interest on deposits, bxz―resource costs for raw materials for production, energy, etc. The terms gxy, dyz, respectively, are the resource profit due to production b and decrease on innovations in management and production automation.

The main task in the study of the dynamics is to identify areas of stable functioning of system (1), as well as areas and conditions for the occurrence of instabilities. This problem is solved by studying the spectrum of small perturbations of stationary states. The stationary state of the system―St (xs, ys, zs) are

$\begin{array}{l}x\left(a-y+bz\right)=0\\ y\left(x-c+dz\right)=0\\ z\left(-e-dy-bx\right)-hx-ky+gxy+r=0\end{array}$ (2)

zeros of the right-hand sides of system (1), i.e. solutions of Equation (2). In this case, there are five stationary solutions, respectively

${\text{St}}_{1}\left(0,0,r/e\right);{\text{St}}_{2}\left(0,{y}_{s},c/d\right);{\text{St}}_{3}\left({x}_{s},0,{z}_{s}\right);{\text{St}}_{4,5}\left({x}_{s},{y}_{s},{z}_{s}\right)$ (3)

By linearizing system (1), we obtain a system of equations for perturbations

${\partial}_{t}\delta {x}_{\alpha}={e}_{\alpha \text{\hspace{0.05em}}\beta}\delta {x}_{\beta},\text{\hspace{1em}}\left\{\delta {x}_{\alpha}\right\}=\left\{\delta x,\delta y,\delta z\right\}$ (4)

where the elements of the evolutionary matrix are equal

$\begin{array}{l}\stackrel{^}{E}=\left\{{e}_{\alpha \beta}\right\};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{e}_{11}=a-{y}_{s}+b{z}_{s}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{e}_{12}=-{x}_{s}=d{z}_{s}-c,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{e}_{13}=b{x}_{s},\\ {e}_{21}={y}_{s}=a+b{z}_{s},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{e}_{22}={x}_{s}-c+d{z}_{s}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{e}_{23}=d{y}_{s},\\ {e}_{31}=-b{z}_{s}-h+g{y}_{s}=-h+ga-bg\left(1-g\right),\\ {e}_{32}=-d{z}_{s}-k+g{x}_{s}=gc-k-d{z}_{s}\left(1+g\right),\\ {e}_{33}=-e-d{y}_{s}-b{x}_{s}=-e-da-bc.\end{array}$ (4)

The condition of the solvability of the system (4) is the spectral equation (SE) equal to

$\mathrm{det}\left(\lambda {\delta}_{\alpha \beta}-{e}_{\alpha \beta}\right)={\lambda}^{3}+s{\lambda}^{2}+p\lambda +q=0$ (5)

where λ is the spectral parameter. The roots of this equation determine the dynamics of the system. Namely: if all the roots of SE have negative real parts, the state under study is stable with respect to small perturbations; if there is at least one root with a positive real part, the state is unstable; if there are roots in the spectrum with positive and negative real parts (saddle focus) and some of them are complex, then this is a special kind of instability―chaos (irregular chaotic oscillations).

$div\stackrel{\to}{G}\left({\partial}_{t}x,{\partial}_{t}y,{\partial}_{t}z\right)=-0.7+0.7x-0.8y+0.9z$ (6)

Expression (6)―the divergence of the vector of the right-hand sides of equations (1) indicates the areas of growth and reduction of the phase volume and is one of the conditions for stability-instability. Depending on the sign $div\stackrel{\to}{G}$ , the dynamics is called active ( $div\stackrel{\to}{G}>0$ ), neutral ( $div\stackrel{\to}{G}=0$ ), dissipative ( $div\stackrel{\to}{G}<0$ ). From (6) it follows that, depending on the magnitude of the resource and the steady state, the system may have different types of dynamics.

3. Spectral and Criterial Analysis of the Dynamics

In view of the complexity of the system (1), only a numerical analysis of the spectrum is possible, for which it is necessary to specify the numerical values of the dimensionless coefficients, as well as select several (one, two) control parameters, varying the values of which monitor the changes in the spectrum and thus the system dynamics. As the control parameter, we select the parameter of the external resource―r, the parameters of innovation in management and production―d, g. The numerical values of the remaining parameters in the system (1) are selected in accordance with the real values on modern production and taking into account the preliminary account. So, in view of the above, the values of the parameters are chosen equal

$\left\{\begin{array}{l}a=0.3,\text{\hspace{0.17em}}b=0.3,\text{\hspace{0.17em}}c=0.8,d=\left(-0.2;0.4\right)\\ e=\left(-0.2;0.2\right),\text{\hspace{0.17em}}h=0.2,\text{\hspace{0.17em}}k=0.4,\text{\hspace{0.17em}}g=\left(0.2;0.4;0.8\right),\text{\hspace{0.17em}}r\in \left(-0.4;0.4\right)\end{array}\right\}$ (7)

The spectral coefficients in (5) taking into account (7) take the form

$\begin{array}{l}s=-\left({e}_{11}+{e}_{22}+{e}_{33}\right),\text{\hspace{1em}}p=-{e}_{12}{e}_{21}-{e}_{13}{e}_{31}-{e}_{23}{e}_{32}+{e}_{11}{e}_{22}+{e}_{11}{e}_{33}+{e}_{23}{e}_{32},\\ q=-{e}_{11}\left({e}_{22}{e}_{33}-{e}_{23}{e}_{32}\right)+{e}_{33}{e}_{12}{e}_{21}+{e}_{22}{e}_{13}{e}_{31}-{e}_{12}{e}_{31}{e}_{23}-{e}_{13}{e}_{21}{e}_{32}.\end{array}$ (8)

For the analysis of SEs, we use both the direct calculation of the roots of the spectral equation and the spectral criteria of instability and chaos given by the authors in [8] [9] [10] . According to these criteria, the relations defining the boundaries of the regions of instability and chaos are of the form

$q<0,\text{\hspace{1em}}K1=q-sp\ge 0;\text{\hspace{1em}}\u0431)\text{\hspace{1em}}p>0;\text{\hspace{1em}}\u0432)\text{\hspace{1em}}K2=q\cdot K1>0$ (9)

The first inequalities in (9) are instability conditions, the second is a sufficient condition for oscillations, the third is the saddle-like character of the spectrum, the realization of all three conditions reflects the presence of chaotic oscillations (chaos) in the system.

Of greatest practical interest are stationary states with nonzero values of quantities, i.e. St_{4}, St_{5}, for which, respectively, taking into account (2), (7), we obtain

$\begin{array}{l}{x}_{4,5s}=c-d{z}_{4,5};\text{}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{y}_{4,5s}=b\left(1+{z}_{4,5}\right);\\ {z}_{4.5s}=\frac{-B\pm \sqrt{{B}^{2}-4AD}}{2A},\text{\hspace{1em}}A=-bdg,\text{\hspace{1em}}D=acg-hc-ka+r\\ B=-da\left(1+g\right)-bc\left(1-g\right)+hd-e.\end{array}$ (10)

The values of the roots of the SE (5) and the criteria q, K1, and K2 for these states, depending on the parameter r, are given in Tables 1-4, respectively.

In Tables 1-2, it is taken into account that at r > 0.18, the states 4.5 cease to be stationary. From Tables 1-2 it also follows that in the whole area of resource change, both with a shortage of credit and with its excess, the system is unstable, there are random fluctuations in it. Moreover, in the St_{4} state, they weakly decay, in the St_{5} state, damped oscillations are transformed into growing ones. In both

Table 1. Estimated values of roots and criteria for the state St_{4} (e = d = −0.2, g = 0.4).

Table 2. Calculated values of the roots and criteria for the state St_{5} (e = d = −0.2, g = 0.4).

Table 3. Calculated values of the roots and criteria for the state St_{4} (e = 0.2, d = −0.2, g = 0.4).

Table 4. Calculated values of the roots and criteria for the state St_{5} (e = 0.2, d = −0.2, g = 0.4).

cases, the signs of the criteria correspond to the signs of the roots and the presence of chaos.

Tables 3-4 reflect the dynamics of the system with increased resource due to management.

It follows from them that the state St_{4} is stable in the entire range of the parameter r, i.e. management maintains state stability regardless of resource change. In contrast, the St_{5} state is more dynamic than the St_{4} state; in it areas of stability
$r\in \left(-0.2;0.0\right)\cup \left(0.2;0.3\right)$ alternate with areas of chaos
$r<-0.2:r\sim 0.1$ .

The tables imply a significant influence on the system dynamics of the parameter d.

Also in these cases, the signs of the criteria correspond to the signs of the roots and the presence of chaos.

Figures 1-4 show the graphs of the dependence of the conditions of instability and chaos on the parameters r, g, d demonstrating the influence of these parameters on the dynamics of the system.

Figure 1. r = −0.2, g = 0.2.

Figure 2. r = −0.2, g = 0.8.

Figure 3. r = 0.2, g = 0.2.

Figure 4. r = 0.2, g = 0.8.

In Figures 1-4 (F = K2 is the condition of chaos, x = 10d). The dependence of the condition of the existence of chaos on the parameter d is shown. From Figures 1-4 it follows that

1) The system is very sensitive to changes in the parameter d;

2) With a small innovation parameter (small value of g parameter), the growth of d parameter stabilizes the system (see Figures 1-3);

3) With an increase in innovation (a larger value of the parameter g), the parameter d destabilizes the system, windows of chaos appear (see Figure 2-4);

4) Significant negative values of the parameter d, i.e. a decrease in the contribution to control and a corresponding increase in accumulations leads to the initiation of chaotic oscillations―the “fever” system.

4. Conclusions

In general, analyzing the results of calculations, the following can be noted:

1) The proposed model provides ample opportunities to study the impact on the economic activity of enterprises of various factors of the production, management and resource sphere.

2) This model describes a very dynamic system, in which, along with the areas of stability, there are areas of instability and chaos, significantly dependent on the parameters of the system.

3) From the tables and graphs it can be seen that the increase in the financing of management stabilizes the processes in the system and that the system loses stability both with a shortage and with an excess of resources.

4) The criteria of instability and chaos proposed by the authors are fully consistent with the calculation of the spectrum, i.e. can be used to predict areas of stability as well as instability of dynamic systems.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Perevoznikov, E. and Lomteva, E. (2019) Modeling of Economic Processes, Instability and Chaos. Journal of Applied Mathematics and Physics, 7, 356-363. https://doi.org/10.4236/jamp.2019.72027

References

- 1. Mazhukin, V.I. and Flinta, M. (2004) Moskovskij gumanitarnyj universitet. Matematicheskoe modelirovanie v ehkonomike, 232 s.
- 2. Markin, M.P. (2007) Matematicheskie metody i modeli v ehkonomike. M. Vysshaya shkola, 422 s.
- 3. Hachatryan, N.K. (2008) Matematicheskoe modelirovanie ehkonomicheskih sistem. M., EHkzamen, 158 s.
- 4. Vihanskij, O.S. and Naumov, O.I. (2010) Menedzhment: chelovek, strategiya, organizaciya, process. M. Delo, 383 s.
- 5. Kabushkin, N.I. (2008) Osnovy menedzhmenta. M.: Novoe znanie, 336 s.
- 6. Klejner, G.B. (2010) Razvitie teorii ehkonomicheskih sistem i ee primenenie v operativnom i strategicheskom upravlenii, WP /2010/CEHMI RAN, 209 s.
- 7. Perevoznikov, E.N. and Skvorcov, G.E. (2018) Analiz dinamiki model’noy sistemy
“Vestnik Sankt—Peterburgskogo gosudarstvennogo universiteta tekhnologii i dizayna. ser. Yestestvennyye i tekhnicheskiye nauki, 2, 41-48. - 8. Perevoznikov, E.N. (2006) Metody analiza ustojchivosti neravnovesnyh system. Izv. Vuzov, Fizika, 10, 34-39.
- 9. Skvorcov, G.E. and Perevoznikov, E.N. (2016) Teoriya neustojchivosti i kriterii haosa. Mezhdunarodnyj nauchno-issledovatel’skij zhurnal, fiziko-matematicheskie nauki, ch.4.-C.98-101.
- 10. Perevoznikov, E.N. and Skvortsov, H.E. (2018) Criteria for Instability and Chaos in Nonlinear System. Journal of Applied Mathematics and Physics, 6, 382-388. https://doi.org/10.4236/jamp.2018.62036