DISCRETE AND CONTINUOUS MODELS AND APPLIED COMPUTATIONAL SCIENCE

The analysis of trajectory dynamics and the solution of optimization problems using computer methods are relevant areas of research in dynamic population-migration models. In this paper, four-dimensional dynamic models describing the processes of competition and migration in ecosystems are studied. Firstly, we consider a modification of the “two competitors—two migration areas” model, which takes into account uniform intraspecific and interspecific competition in two populations as well as non-uniform bidirectional migration in both populations. Secondly, we consider a modification of the “two competitors—two migration areas” model, in which intraspecific competition is uniform and interspecific competition and bidirectional migration are non-uniform. For these two types of models, the study is carried out taking into account the variability of parameters. The problems of searching for model parameters based on the implementation of two optimality criteria are solved. The first criterion of optimality is associated with the fulfillment of such a condition for the coexistence of populations, which in mathematical form is the integral maximization of the functions product characterizing the populations densities. The second criterion of optimality involves checking the assumption of the such a four-dimensional positive vector existence, which will be a state of equilibrium. The algorithms developed on the basis of the first and second optimality criteria using the differential evolution method result in optimal sets of parameters for the studied population-migration models. The obtained sets of parameters are used to find positive equilibrium states and analyze trajectory dynamics. Using the method of constructing self-consistent one-step models and an automated stochastization procedure, the transition to the stochastic case is performed. The structural description and the possibility of analyzing two types of population-migration stochastic models are provided by obtaining Fokker–Planck equations and Langevin equations with corresponding coefficients. Algorithms for generating trajectories of the Wiener process, multipoint distributions and modifications of the Runge–Kutta method are used. A series of computational experiments is carried out using a specialized software package whose capabilities allow for the construction and analysis of dynamic models of high dimension, taking into account the evaluation of the stochastics influence. The trajectory dynamics of two types of population-migration models are investigated, and a comparative analysis of the results is carried out both in the deterministic and stochastic cases. The results can be used in the modeling and optimization of dynamic models in natural science.

Our group has been investigating kinetic models for quite a long time. The structure of classical kinetic models is described by rather simple assumptions about the interaction of the entities under study. Also, the construction of kinetic equations (both stochastic and deterministic) is based on simple sequential steps. However, in each step, the researcher must manipulate a large number of elements. And once the differential equations are obtained, the problem of solving or investigating them arises. The use of symbolic-numeric approach methodology is naturally directed. When the input is an information model of the system under study, represented in some diagrammatic form. And as a result, we obtain systems of differential equations (preferably, in all possible variants). Then, as part of this process, we can investigate the resulting equations (by a variety of methods). We have previously taken several steps in this direction, but we found the results somewhat unsatisfactory. At the moment we have settled on the package Catalyst. jl, which belongs to the Julia language ecosystem. The authors of the package declare its relevance to the field of chemical kinetics. Whether it is possible to study more complex systems with this package, we cannot say. Therefore, we decided to investigate the possibility of using this package for our models to begin with standard problems of chemical kinetics. As a result, we can summarize that this package seems to us to be the best solution for the symbolic-numerical study of chemical kinetics problems.