Все выпуски

В статье рассмотрены режимы с обострением в социальной и биологической истории. Проведен анализ возможных причин резкого ускорения биологических и социальных процессов в определенные исторические эпохи. С использованием математического моделирования показано, что гиперболические тренды в социальной и биологической эволюции могут быть следствием переходных процессов в периоды расширения экологических ниш. Ускорение биологического видообразования связано с тем, что более ранние виды своей жизнедеятельностью изменяют среду обитания, делая ее более разнообразной, насыщая органикой, тем самым создавая благоприятные условия для появления новых видов. В социальной истории расширение экологических ниш связано с технологическими революциями, важнейшими из которых были: неолитическая революция — переход от присваивающего хозяйства к производящему (10 тыс. лет назад), «городская революция» — переход от неолита к бронзовому веку (5 тыс. лет назад), «осевое время» — переход к массовому освоению железных орудий (2.5 тыс. лет назад), промышленная революция — переход от ручного труда к машинному (200 лет назад). Все эти технологические революции сопровождались резким демографическим ростом, изменениями в социальной и политической сфе- рах. Так, наблюдаемый в последние столетия гиперболический характер роста некоторых демографических, экономических и других показателей мировой динамики — это следствие переходных процессов, начавшихся вследствие промышленной революции (замены ручного труда машинным) и предваряющих переход общества на новую стадию своего развития. Точка сингулярности гиперболического тренда характеризует окончание начального этапа этого процесса и переход к завершающей его стадии. Предложена математическая модель, описывающая демографические и экономические изменения в эпохи перемен. Показано, что прямым аналогом современной ситуации в этом смысле является «осевое время» (период с 8 века до нашей эры до начала нашей эры). Наличие такой аналогии позволяет заглянуть в будущее, изучая прошлое.

The article describes the modes with the exacerbation of social and biological history. The analysis of the possible causes of the sharp acceleration of biological and social processes in certain historical periods is carried out. Using mathematical modeling shows that hyperbolic trends in social and biological evolution may be the result of transitional processes in periods of expansion of ecological niches. Accelerating biological speciation due to the fact that its earlier life change inhabitancy, making it more diverse, saturating the organic, thus creating favourable conditions for the emergence of new species. In the social history of the expansion of ecological niches associated with technological revolutions, of which the most important were: Neolithic revolution — the transition from appropriating economy to producing economy (10 thousand years ago), “urban revolution” — a shift from the Neolithic epoch to the bronze epoch (5 thousand years ago), the “axial age” — transition to the development of iron tools (2.5 thousand years ago), the industrial revolution — the transition from manual labor to machine production (200 years ago). All of these technological revolutions have been accompanied by dramatic population growth, changes in social and political spheres. So, observed in the last century, hyperbolic nature of some demographic, economic growth and other indicators of world dynamics is a consequence of the transition process, which began as a result of the industrial revolution and to prepare for the transition of the society to a new stage of its development. Singularity point of hyperbolic trend shows the end of the initial phase of the process and marks the transition to the final stage. The mathematical model describing the demographic and economic changes in the era of change is proposed. It is shown that a direct analogue of the contemporary situation in this sense is the “axial age” (since 8 century BC to the beginning of our era). The existence of this analogy allows you to see into the future by studying the past.

Автоматическое распознавание личности по биометрическому признаку основано на уникальных особенностях или характеристиках людей. Процесс биометрической идентификации представляет собой формирование эталонных шаблонов и сравнение их с новыми входными данными. Алгоритмы распознавания по рисунку радужной оболочки глаза показали на практике высокую точность и малый процент ошибок идентификации. Преимущества радужки над другими биометрическими признаками определяется ее большей степенью свободы (около 249 степеней свободы), избыточной плотностью уникальных признаков и постоянностью во времени. Высокий уровень достоверности распознавания очень важен, потому что позволяет выполнять поиск по большим базам данных и работать в режиме идентификации один-ко-многим, в отличии от режима проверки один-к-одному, который применим дляне большого количества сравнений. Любая биометрическая система идентификации является вероятностной. Для описания качественных характеристик распознавания применяются: точность распознавания, вероятность ложного доступа и вероятность ложного отказа доступа. Эти характеристики позволяют сравнивать методы распознавания личности между собой и оценивать поведение системы в каких-либо условиях. В этой статье объясняется математическая модель биометрической идентификации по радужной оболочке глаза, ее характеристики и анализируются результаты сравнения модели с реальным процессом распознавания. Для решения этой задачи проводится обзор существующих методов идентификации по радужной оболочке глаза, основанных на различных способах формирования вектора уникальных признаков. Описывается разработанный программный комплекс на языке Python, который строит вероятностные распределения и генерирует большие наборы тестовых данных, которые могут быть использованы в том числе для обучения нейронной сети принятия решения об идентификации. В качестве практического применения модели предложен алгоритм синергии нескольких методов идентификации личности по радужной оболочке глаза, позволяющий увеличить качественные характеристики системы, в сравнении с применением каждого метода отдельно.

Automatic recognition of personal identity by biometric features is based on unique peculiarities or characteristics of people. Biometric identification process consist in making of reference templates and comparison with new input data. Iris pattern recognition algorithms presents high accuracy and low identification errors percent on practice. Iris pattern advantages over other biometric features are determined by its high degree of freedom (nearly 249), excessive density of unique features and constancy. High recognition reliability level is very important because it provides search in big databases. Unlike one-to-one check mode that is applicable only to small calculation count it allows to work in one-to-many identification mode. Every biometric identification system appears to be probabilistic and qualitative characteristics description utilizes such parameters as: recognition accuracy, false acceptance rate and false rejection rate. These characteristics allows to compare identity recognition methods and asses the system performance under any circumstances. This article explains the mathematical model of iris pattern biometric identification and its characteristics. Besides, there are analyzed results of comparison of model and real recognition process. To make such analysis there was carried out the review of existing iris pattern recognition methods based on different unique features vector. The Python-based software package is described below. It builds-up probabilistic distributions and generates large test data sets. Such data sets can be also used to educate the identification decision making neural network. Furthermore, synergy algorithm of several iris pattern identification methods was suggested to increase qualitative characteristics of system in comparison with the use of each method separately.

В данной работе с помощью алгоритма Метрополиса авторами были изучены магнитные системы, в которых из-за конкуренции между прямым гейзенберговским обменом и взаимодействием Дзялошинского–Мория возникают магнитные вихревые структуры — скирмионы.

В статье рассматриваются условия зарождения и стабильного существования магнитных скирмионов в двумерных магнитных пленках в рамках классической модели Гейзенберга. Изучена термическая стабильность скирмионов в магнитной пленке. Были рассмотрены процессы формирования различных состояний в изучаемой системе при варьировании величины внешнего магнитного поля, выделены различные фазы, в которые переходит система спинов Гейзенберга. Было выделено семь фаз: парамагнитная, спиральная, лабиринтная, спираль-скирмионная, скирмионная, скирмион-ферромагнитная и ферромагнитная фазы, подробный анализ конфигураций которых приводится в статье.

Построены две фазовые диаграммы: на первой показано поведение системы при постоянном $D$ в зависимости от величин внешнего магнитного поля и температуры: $(T, B)$, на второй — изменение кон- фигураций системы при постоянной температуре $T$ в зависимости от величины взаимодействия Дзялошинского–Мории и внешнего магнитного поля: $(D, B)$.

Полученные в ходе численных экспериментов данные будут использованы в дальнейших исследованиях при определении модельных параметров системы для формирования стабильного скирмионного состояния и разработки методов контроля скирмионов в магнитной пленке.

Magnetic systems, in which due to competition between the direct Heisenberg exchange and the Dzyaloshinskii –Moriya interaction, magnetic vortex structures — skyrmions appear, were studied using the Metropolis algorithm.

The conditions for the nucleation and stable existence of magnetic skyrmions in two-dimensional magnetic films in the frame of the classical Heisenberg model were considered in the article. A thermal stability of skyrmions in a magnetic film was studied. The processes of the formation of various states in the system at different values of external magnetic fields were considered, various phases into which the Heisenberg spin system passes were recognized. The authors identified seven phases: paramagnetic, spiral, labyrinth, spiralskyrmion, skyrmion, skyrmion-ferromagnetic and ferromagnetic phases, a detailed analysis of the configurations is given in the article.

Two phase diagrams were plotted: the first diagram shows the behavior of the system at a constant $D$ depending on the values of the external magnetic field and temperature $(T, B)$, the second one shows the change of the system configurations at a constant temperature $T$ depending on the magnitude of the Dzyaloshinskii – Moriya interaction and external magnetic field: $(D, B)$.

The data from these numerical experiments will be used in further studies to determine the model parameters of the system for the formation of a stable skyrmion state and to develop methods for controlling skyrmions in a magnetic film.

Целью работы является исследование монотонной конечно-разностной схемы второго порядка точности, созданной на основе обобщения одномерной Z-схемы. Исследование проведено для модельных уравнений переноса несжимаемой среды. В работе описано двумерное обобщение Z-схемы с нелинейной коррекцией, использующей вместо потоков косые разности, содержащие значения из разных временных слоев. Численно проверена монотонность полученной нелинейной схемы для функций-ограничителей двух видов, как для гладких решений, так и для негладких, и получены численные оценки порядка точности построенной схемы. Построенная схема является абсолютно устойчивой, но теряет свойство монотонности при превышении шага Куранта. Отличительной особенностью предложенной конечно-разностной схемы является минимальность ее шаблона.

Построенная численная схема предназначена для моделей плазменных неустойчивостей различных масштабов в низкоширотной ионосферной плазме Земли. Одна из реальных задач, при решении которых возникают подобные уравнения, — это численное моделирование сильно нестационарных среднемасштабных процессов в земной ионосфере в условиях возникновения неустойчивости Рэлея – Тейлора и плазменных структур с меньшими масштабами, механизмами генерации которых являются неустойчивости других типов, что приводит к явлению F-рассеяния. Вследствие того, что процессы переноса в ионосферной плазме контролируются магнитным полем, в поперечном к магнитному полю направле- нии предполагается выполнение условия несжимаемости плазмы.

The aim of the work is to study a monotone finite-difference scheme of the second order of accuracy, created on the basis of a generalization of the one-dimensional Z-scheme. The study was carried out for model equations of the transfer of an incompressible medium. The paper describes a two-dimensional generalization of the Z-scheme with nonlinear correction, using instead of streams oblique differences containing values from different time layers. The monotonicity of the obtained nonlinear scheme is verified numerically for the limit functions of two types, both for smooth solutions and for nonsmooth solutions, and numerical estimates of the order of accuracy of the constructed scheme are obtained.

The constructed scheme is absolutely stable, but it loses the property of monotony when the Courant step is exceeded. A distinctive feature of the proposed finite-difference scheme is the minimality of its template. The constructed numerical scheme is intended for models of plasma instabilities of various scales in the low-latitude ionospheric plasma of the Earth. One of the real problems in the solution of which such equations arise is the numerical simulation of highly nonstationary medium-scale processes in the earth’s ionosphere under conditions of the appearance of the Rayleigh – Taylor instability and plasma structures with smaller scales, the generation mechanisms of which are instabilities of other types, which leads to the phenomenon F-scattering. Due to the fact that the transfer processes in the ionospheric plasma are controlled by the magnetic field, it is assumed that the plasma incompressibility condition is fulfilled in the direction transverse to the magnetic field.

Изменения численности y образующихся популяций могут развиваться по нескольким динамическим сценариям. Для стремительных биологических инвазий оказывается важным фактор времени выработки реакции противодействия со стороны биотического окружения. Известны два классических эксперимента с разным завершением противоборства биологических видов. В опытах Гаузе с инфузориями вселенный хищник после кратких осцилляций полностью уничтожал свой ресурс, так его $r$-параметр для созданных условий стал избыточен. Собственная репродуктивная активность не регулировалась дополнительными факторами и в результате становилась критичной для вселенца. В экспериментах Утиды с жуками и выпущенными паразитическими осами виды сосуществовали. В ситуации, когда популяцию с высоким репродуктивным потенциалом регулируют несколько естественных врагов, могут возникать интересные динамические эффекты, наблюдавшиеся у фитофагов в вечнозеленом лесу Австралии. Паразитические перепончатокрылые, конкурируя между собой, создают для быстро размножающихся вредителей псиллид систему регуляции с запаздыванием, когда допускается быстрое увеличение локальной популяции, но не превышающее порогового значения численности вредителя. В работе предложена модель на основе дифференциального уравнения с запаздыванием, описывающая сценарий адаптационной регуляции для популяции с большим репродуктивным потенциалом при активном, но запаздывающем противодействии с пороговой регуляцией данного вновь возникшего воздействия. За кратким максимумом следует быстрое сокращение численности, но минимизация не становится критической для популяции. Показано, что усложнение функции регуляции биотического противодействия приводит к стабилизации динамики после прохождения минимума численности быстро размножающимся видом. Для гибкой системы переходные режимы «рост/кризис» ведут к поиску нового равновесия в эволюционном противостоянии.

Changes in abundance for emerging populations can develop according to several dynamic scenarios. After rapid biological invasions, the time factor for the development of a reaction from the biotic environment will become important. There are two classic experiments known in history with different endings of the confrontation of biological species. In Gause’s experiments with ciliates, the infused predator, after brief oscillations, completely destroyed its resource, so its $r$-parameter became excessive for new conditions. Its own reproductive activity was not regulated by additional factors and, as a result, became critical for the invader. In the experiments of the entomologist Uchida with parasitic wasps and their prey beetles, all species coexisted. In a situation where a population with a high reproductive potential is regulated by several natural enemies, interesting dynamic effects can occur that have been observed in phytophages in an evergreen forest in Australia. The competing parasitic hymenoptera create a delayed regulation system for rapidly multiplying psyllid pests, where a rapid increase in the psyllid population is allowed until the pest reaches its maximum number. A short maximum is followed by a rapid decline in numbers, but minimization does not become critical for the population. The paper proposes a phenomenological model based on a differential equation with a delay, which describes a scenario of adaptive regulation for a population with a high reproductive potential with an active, but with a delayed reaction with a threshold regulation of exposure. It is shown that the complication of the regulation function of biotic resistance in the model leads to the stabilization of the dynamics after the passage of the minimum number by the rapidly breeding species. For a flexible system, transitional regimes of growth and crisis lead to the search for a new equilibrium in the evolutionary confrontation.

Фоновая социальная напряженность общества может быть количественно оценена по различным статистическим индикаторам. Модели, прогнозирующие динамику социальной напряженности, успешно применяются для описания различных социальных процессов. Когда количество рассматриваемых групп общества мало, динамику соответствующих индикаторов можно описать при помощи системы обыкновенных дифференциальных уравнений. При увеличении количества взаимодействующих элементов резко возрастает сложность задач, что существенно затрудняет их аналитическое исследование. Модель сплошной социальной стратификации получаетсяв результате перехода от дискретной цепочки взаимодействующих социальных слоев к их непрерывному распределению на некотором интервале, то есть перехода к модели сплошной среды. В этом случае напряженность распространяется локально, но в действительности элита общества влияет на все слои через средства массовой информации, а также интернет позволяет влиять всем группам на другие. Эти факторы можно учесть через слагаемое модели, описывающее негативное внешнее воздействие. В настоящей работе предложена модель сплошной социальной стратификации, описывающая динамику системы из двух социумов, связанных через процесс миграции населения. Предполагается, что из социального слоя системы-донора с наибольшей напряженностью происходит отток людей, переносящих свою напряженность в систему-акцептор, причем при миграции люди попадают в более бедные слои принимающего общества. Рассматриваетсяслуч ай пространственно однородных коэффициентов, что соответствует частному случаю небольшого социума. При помощи метода конечных объемов построена пространственнаяди скретизация задачи, корректно отражающая конечную скорость распространения напряженности в обществе. Выполнена проверка выбранной дискретизации путем сравненияч исленного решения с точными решениями вспомогательного уравнения нелинейной диффузии. Проведено численное исследование системы с миграцией при различных значениях параметров, проанализировано влияние интенсивности миграции на принимающее общество, найдены условия дестабилизации общества акцептора под влиянием миграции. Полученные в работе результаты могут быть применены при дальнейшем исследовании модели в случае пространственно неоднородных коэффициентов, что соответствует более реалистичной картине общества.

The background social strain of a society can be quantitatively estimated using various statistical indicators. Mathematical models, allowing to forecast the dynamics of social strain, are successful in describing various social processes. If the number of interacting groups is small, the dynamics of the corresponding indicators can be modelled with a system of ordinary differential equations. The increase in the number of interacting components leads to the growth of complexity, which makes the analysis of such models a challenging task. A continuous social stratification model can be considered as a result of the transition from a discrete number of interacting social groups to their continuous distribution in some finite interval. In such a model, social strain naturally spreads locally between neighbouring groups, while in reality, the social elite influences the whole society via news media, and the Internet allows non-local interaction between social groups. These factors, however, can be taken into account to some extent using the term of the model, describing negative external influence on the society. In this paper, we develop a continuous social stratification model, describing the dynamics of two societies connected through migration. We assume that people migrate from the social group of donor society with the highest strain level to poorer social layers of the acceptor society, transferring the social strain at the same time. We assume that all model parameters are constants, which is a realistic assumption for small societies only. By using the finite volume method, we construct the spatial discretization for the problem, capable of reproducing finite propagation speed of social strain. We verify the discretization by comparing the results of numerical simulations with the exact solutions of the auxiliary non-linear diffusion equation. We perform the numerical analysis of the proposed model for different values of model parameters, study the impact of migration intensity on the stability of acceptor society, and find the destabilization conditions. The results, obtained in this work, can be used in further analysis of the model in the more realistic case of inhomogeneous coefficients.

В работе исследуется влияние быстрого локального выделения тепла вблизи обтекаемой сверхзвуковым потоком газа (воздуха) поверхности на область отрыва, возникающую при быстром его повороте. Данная поверхность состоит из двух плоскостей, образующих при пересечении тупой угол, так что при обтекании этой поверхности сверхзвуковой поток газа поворачивается на положительный угол, что формирует косой скачок уплотнения, взаимодействующий с пограничным слоем и вызывающий отрыв потока. Быстрый локальный нагрев газа над обтекаемой поверхностью моделирует протяженный искровой разряд субмикросекундной длительности, пересекающий поток. Газ, нагретый в зоне разряда, взаимодействует с областью отрыва. Течение можно считать плоским, поэтому численное моделирование проводится в двумерной постановке. Численное моделирование проведено для ламинарного режима течения с использованием солвера sonicFoam пакета программ OpenFOAM.

В работе описан способ построения двумерной расчетной сетки с использованием шестигранных ячеек. Выполнено исследование сеточной сходимости. Приводится методика задания начальных профилей параметров течения на входе в расчетную область, позволяющая сократить время счета при уменьшении количества расчетных ячеек. Описан способ нестационарного моделирования процесса быстрого локального нагрева газа, заключающегося в наложении дополнительных полей повышенных значений давления и температуры, вычисленных из величины энергии, вложенной в набегающий сверхзвуковой поток газа, на соответствующие поля величин, предварительно полученные в стационарном случае. Параметры энерговклада в поток, соответствующие параметрам процесса инициирования электрического разряда, а также параметры набегающего потока близки к экспериментальным величинам.

При анализе данных численного моделирования получено, что быстрый локальный нагрев приводит к возникновению газодинамического возмущения (квазицилиндрической ударной волны и нестационарного завихренного течения), которое при взаимодействии с областью отрыва приводит к смещению точки отрыва вниз по потоку. В работе рассмотрен вопрос о влиянии энергии, затраченной на локальный нагрев газа, и положения места нагрева относительно точки отрыва на величину максимального ее смещения.

The influence of the process of initiating a rapid local heat release near surface streamlined by supersonic gas (air) flow on the separation region that occurs during a fast turn of the flow was investigated. This surface consists of two planes that form obtuse angle when crossing, so that when flowing around the formed surface, the supersonic gas flow turns by a positive angle, which forms an oblique shock wave that interacts with the boundary layer and causes flow separation. Rapid local heating of the gas above the streamlined surface simulates long spark discharge of submicrosecond duration that crosses the flow. The gas heated in the discharge zone interacts with the separation region. The flow can be considered two-dimensional, so the numerical simulation is carried out in a two-dimensional formulation. Numerical simulation was carried out for laminar regime of flow using the sonicFoam solver of the OpenFOAM software package.

The paper describes a method for constructing a two-dimensional computational grid using hexagonal cells. A study of grid convergence has been carried out. A technique is given for setting the initial profiles of the flow parameters at the entrance to the computational domain, which makes it possible to reduce the computation time by reducing the number of computational cells. A method for non-stationary simulation of the process of rapid local heating of a gas is described, which consists in superimposing additional fields of increased pressure and temperature values calculated from the amount of energy deposited in oncoming supersonic gas flow on the corresponding fields of values obtained in the stationary case. The parameters of the energy input into the flow corresponding to the parameters of the electric discharge process, as well as the parameters of the oncoming flow, are close to the experimental values.

During analyzing numerical simulation data it was found that the initiation of rapid local heating leads to the appearance of a gas-dynamic perturbation (a quasi-cylindrical shock wave and an unsteady swirling flow), which, when interacting with the separation region, leads to a displacement of the separation point downstream. The paper considers the question of the influence of the energy spent on local heating of the gas, and of the position on the streamlined surface of the place of heating relative to the separation point, on the value of its maximum displacement.

The viscoelastic fluid flow model across a porous medium has captivated the interest of many contemporary researchers due to its industrial and technical uses, such as food processing, paper and textile coating, packed bed reactors, the cooling effect of transpiration and the dispersion of pollutants through aquifers. This article focuses on the influence of variable viscosity and viscoelasticity on the magnetohydrodynamic oscillatory flow of second-order fluid through thermally radiating wavy walls. A mathematical model for this fluid flow, including governing equations and boundary conditions, is developed using the usual Boussinesq approximation. The governing equations are transformed into a system of nonlinear ordinary differential equations using non-similarity transformations. The numerical results obtained by applying finite-difference code based on the Lobatto IIIa formula generated by bvp4c solver are compared to the semi-analytical solutions for the velocity, temperature and concentration profiles obtained using the homotopy perturbation method (HPM). The effect of flow parameters on velocity, temperature, concentration profiles, skin friction coefficient, heat and mass transfer rate, and skin friction coefficient is examined and illustrated graphically. The physical parameters governing the fluid flow profoundly affected the resultant flow profiles except in a few cases. By using the slope linear regression method, the importance of considering the viscosity variation parameter and its interaction with the Lorentz force in determining the velocity behavior of the viscoelastic fluid model is highlighted. The percentage increase in the velocity profile of the viscoelastic model has been calculated for different ranges of viscosity variation parameters. Finally, the results are validated numerically for the skin friction coefficient and Nusselt number profiles.

The viscoelastic fluid flow model across a porous medium has captivated the interest of many contemporary researchers due to its industrial and technical uses, such as food processing, paper and textile coating, packed bed reactors, the cooling effect of transpiration and the dispersion of pollutants through aquifers. This article focuses on the influence of variable viscosity and viscoelasticity on the magnetohydrodynamic oscillatory flow of second-order fluid through thermally radiating wavy walls. A mathematical model for this fluid flow, including governing equations and boundary conditions, is developed using the usual Boussinesq approximation. The governing equations are transformed into a system of nonlinear ordinary differential equations using non-similarity transformations. The numerical results obtained by applying finite-difference code based on the Lobatto IIIa formula generated by bvp4c solver are compared to the semi-analytical solutions for the velocity, temperature and concentration profiles obtained using the homotopy perturbation method (HPM). The effect of flow parameters on velocity, temperature, concentration profiles, skin friction coefficient, heat and mass transfer rate, and skin friction coefficient is examined and illustrated graphically. The physical parameters governing the fluid flow profoundly affected the resultant flow profiles except in a few cases. By using the slope linear regression method, the importance of considering the viscosity variation parameter and its interaction with the Lorentz force in determining the velocity behavior of the viscoelastic fluid model is highlighted. The percentage increase in the velocity profile of the viscoelastic model has been calculated for different ranges of viscosity variation parameters. Finally, the results are validated numerically for the skin friction coefficient and Nusselt number profiles.

В работе исследуется дискретная модификация модели А.П. Михайлова «власть – общество», ранее предложенная автором. Эта модификация основана на стохастическом клеточном автомате, то есть имеет микродинамику, принципиально отличную от базовой непрерывной, основанной на дифференциальных уравнениях модели. При этом макродинамика дискретной модификации, как показано в предыдущих работах, совпадает с макродинамикой исходной модели. Этот важный результат, однако, вызывает вопрос, в чем смысл использования дискретной модели. Ее главной особенностью является гибкость, позволяющая добавлять в рассмотрение самые разные факторы, учет которых в непрерывной модели либо приводит к существенному росту вычислительной сложности, либо в принципе невозможен.

В данной работе рассматриваются несколько примеров подобного расширения области применимости модели, при помощи которого решается ряд прикладных задач.

Одна из модификаций модели учитывает экономические связи между регионами и муниципалитетами, что не могло быть исследовано в базовой модели. Вычислительные эксперименты подтвердили улучшение социально-экономических показателей системы при наличии таких связей.

Вторая модификация включает в себя возможность внутренней миграции в системе. С ее помощью был получен ряд результатов, связанных с социально-экономическим развитием более благополучного региона, притягивающего мигрантов.

Кроме этого, была исследована динамика системы при изменении количества регионов и муниципалитетов в системе. Показано негативное влияние этого процесса на социально-экономические показатели системы и найдено возможное управление, имеющее целью преодоление этого негативного влияния.

Результатами данного исследования, таким образом, явились как решение отдельных прикладных задач, так и демонстрация на их примере более широких возможностей дискретной модели по сравнению с базовой непрерывной.

The paper considers an earlier proposed by the author discrete modification of the A. P. Mikhailov “power – society” model. The modification is based on a stochastic cellular automaton, it’s microdynamics being completely different from the c continuous model based on differential equations. However, the macrodynamics of the discrete modification is shown in previous works to be equivalent to one of the continuous model. This is important, but at the same time raises the question why use the discrete model. The answer lies in its flexibility, which allows adding a variety of factors, the consideration of which in a continuous model either leads to a significant increase in computational complexity or is simply impossible.

This paper considers several examples of such applicability expansion of the model, with the help of which a number of applied problems are solved.

One of the modifications of the model takes into account economic ties between regions and municipalities, which could not be studied in the basic model. Computational experiments confirmed the improvement of the socio-economic indicators of the system under the influence of the ties.

The second modification allows internal migration in the system. Using it we studied the socio-economic development of a more prosperous region that attracts migrants.

Next we studied the dynamics of the system while the number of regions and municipalities changes. The negative impact of this process on the socio-economic indicators of the system was shown and possible control was found to overcome this negative impact.

The results of this study, therefore, include both the solution of some applied problems and the demonstration of the broader applicability of the discrete model compared with the continuous one.

В работе рассмотрены проблемы математического описания деонтологических аспектов, влияющих на поведение ЛПР (лиц, принимающих решения). Предложена методология соотнесения утилитарных (материальных) и деонтологических (ценностных) аспектов при принятии ими решений с учетом их психологических особенностей. Предложена математическая модель совместного учета утилитарных и деонтологических факторов при принятии ЛПР решений в различных ситуациях. Выявлены некоторые закономерности, связанные с этим учетом, приведено их формальное описание. Модель показывает, что существует тенденция постепенного снижения уровня деонтологичности в оценке альтернатив при принятии решений (по сравнению с тем, к чему склоняет внешний мир) к большей утилитарности. Эта тенденция с течением времени начинает влиять на общественное мнение и на отношение общества к моральным нормам, постепенно снижая общий уровень моральности в обществе. Остановить этот процесс можно только путем постоянного и целенаправленного поддержания обществом и государством высокого уровня деонтологичности (идеологическая работа, пропаганда традиционных ценностей, воспитательная работа в школе и т. п.), в противном случае общество с неизбежностью со временем станет утилитарным, ориентирующимся при принятии решений исключительно на материальные факторы.

В дальнейшем планируется использовать разработанный инструментарий для анализа конкретных ситуаций, в том числе для анализа закономерностей цивилизационных циклов: взлета и падения Римской империи, СССР, современной западной цивилизации.

The paper considers the problems of mathematical description of deontological aspects influencing the behavior of decision makers. A methodology is proposed for correlating utilitarian (material) and deontological (value) aspects in their decision-making, taking into account their psychological characteristics. A mathematical model is proposed for the joint consideration of utilitarian and deontological factors in decision-making in various situations. Some patterns related to this consideration are identified, and their formal description is given. The model shows that there is a tendency for a gradual decrease in the level of deontology in evaluating alternatives when making decisions (compared to what the outside world inclines to) towards greater utilitarianism. Over time, this trend begins to influence public opinion and society’s attitude to moral norms, gradually reducing the overall level of morality in society. This process can be stopped only by constantly and purposefully maintaining a high level of deontology by society and the state (ideological work, promotion of traditional values, educational work at school, etc.), otherwise society will inevitably become utilitarian over time, focusing exclusively on material factors when making decisions.

In the future, it is planned to use the developed tools for analyzing specific situations, including for analyzing the patterns of civilizational cycles: the rise and fall of the Roman Empire, the USSR, and modern Western civilization).