Все выпуски

Рассматривается модель хаотического движения пластинки в вязкой жидкости, описываемая колебательной системой трех обыкновенных дифференциальных уравнений с квадратичной нелинейностью. В ходе бифуркационного исследования особых точек системы построены карты типов особых точек и найдено уравнение поверхности в пространстве параметров диссипации и циркуляции, на которой происходит бифуркация Андронова-Хопфа рождения предельного цикла. При дальнейшем изменении параметров вблизи поверхности Андронова-Хопфа найдены каскады бифуркаций удвоения периода цикла Фейгенбаума и субгармонические каскады Шарковского, заканчивающиеся рождением цикла периода три. Получены выражения для седловых чисел седлоузла и двух седлофокусов и построены их графики в пространстве параметров. Показано, что в системе реализуются гомоклинические каскады бифуркаций при разрушении гомоклинических траекторий седлофокусов. Существование гомоклинических траекторий седлофокусов доказано численно-аналитическим методом. Графики старшего показателя Ляпунова и бифуркационные диаграммы показывают, что при изменении коэффициентов диссипации система в несколько этапов переходит к хаосу.

We consider the model of chaotic motion of a plate in a viscous fluid, described by an oscillatory system of three ordinary differential equations with a quadratic nonlinearity. In the course of the bifurcation study of singular points of the system, maps of the types of singular points are constructed and a surface equation is found in the space of dissipation and circulation parameters on which the Andronov-Hopf bifurcation of the limit cycle creation takes place. With a further change in the parameters near the Andronov-Hopf surface, cascades of the period doubling doubling of the Feigenbaum cycle and the Sharkovsky subharmonic cascades, ending with the creation of a cycle of period three, are found. Expressions are obtained for saddle numbers of the saddle-node and two saddle-foci and their plots are plotted in the parameter space. It is shown that homoclinic cascades of bifurcations are realized in the system with the destruction of homoclinic trajectories of saddle-foci. The existence of homoclinic trajectories of saddle-foci is proved by a numerical-analytical method. The graphs of the largest Lyapunov exponent and the bifurcation diagrams show that when the dissipation coefficients change, the system switches to chaos in several stages.

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

The paper is concerned with the approximation of the value function of the zero-sum differential game with the minimal cost, i.e., the differential game with the payoff functional determined by the minimization of some quantity along the trajectory by the solutions of continuous-time stochastic games with the stopping governed by one player. Notice that the value function of the auxiliary continuous-time stochastic game is described by the Isaacs–Bellman equation with additional inequality constraints. The Isaacs–Bellman equation is a parabolic PDE for the case of stochastic differential game and it takes a form of system of ODEs for the case of continuous-time Markov game. The approximation developed in the paper is based on the concept of the stochastic guide first proposed by Krasovskii and Kotelnikova.

Рассматривается динамическая система сдвигов в пространстве ℜ непрерывных функций, принимающих значения в полном метрическом пространстве (clos(Rⁿ), ρ_cl)
непустых замкнутых подмножеств в Rⁿ. Расстояние между функциями в этом пространстве определяется с помощью аналога метрики Бебутова в пространстве вещественных функций, определенных и непрерывных на всей числовой оси. Показано, что для компактности замыкания траектории точки в ℜ достаточно, чтобы исходная функция была ограничена и равномерно непрерывна в метрике ρ_cl. Как следствие, доказано, что замыкание траектории рекуррентного движения или траектории почти периодического движения в ℜ компактно.

In the work there is considered the dynamical system of translations in the space ℜ of continuous multi-valued functions with images in complete metric space (clos(Rⁿ), ρ_cl)
of nonempty closed subsets of Rⁿ. The distance between such functions is measured by means of the metric analogous to the Bebutov metric constructed for the space of continuous real-valued functions defined on the whole real line. It is shown that for compactness of the trajectory’s closure in ℜ it is sufficient to have initial function bounded and uniformly continuous in the ρ_cl metric. As consequence, it is also proved that the trajectory’s closure of a recurrent or an almost periodic motion is compact in ℜ.

Рассматриваются две задачи нелинейного гарантированного оценивания фазовых состояний динамических систем. Предполагается, что неизвестные измеримые по $t$ возмущения линейно входят в уравнение движения и аддитивно — в уравнения измерения. Эти возмущения стеснены нелинейными интегральными функционалами, один из которых является аналогом функционала обобщенной работы. Исследуемая задача состоит в построении информационных множеств по данным измерения, содержащих истинное положение траектории. Используется подход динамического программирования. Если для первого функционала требуется решить нелинейное уравнение в частных производных первого порядка, что не всегда возможно, то для функционала обобщенной работы достаточно найти решение линейного уравнения Ляпунова первого порядка, что существенно упрощает задачу. Тем не менее, даже в этом случае приходится налагать дополнительные условия на параметры системы для того, чтобы траектория системы, соответствующая наблюдаемому сигналу, существовала. Если уравнение движения линейно по фазовой переменной, то многие предположения выполняются автоматически. Для этого случая обсуждается вопрос о взаимной оценке сверху и снизу информационных множеств по включению для разных функционалов. В заключение рассмотрен наиболее прозрачный линейно-квадратичный случай. Изложение иллюстрируется примерами.

Two problems of nonlinear guaranteed estimation for states of dynamical systems are considered. It is supposed that unknown measurable in $t$ disturbances are linearly included in the equation of motion and are additive in the measurement equations. These disturbances are constrained by nonlinear integral functionals, one of which is analog of functional of the generalized work. The studied problem consists in creation of the information sets according to measurement data containing the true position of the trajectory. The dynamic programming approach is used. If the first functional requires solving a nonlinear equation in partial derivatives of the first order which is not always possible, then for functional of the generalized work it is enough to find a solution of the linear Lyapunov equation of the first order that significantly simplifies the problem. Nevertheless, even in this case it is necessary to impose additional conditions on the system parameters in order for the system trajectory of the observed signal to exist. If the motion equation is linear in state variable, then many assumptions are carried out automatically. For this case the issue of mutual approximation of information sets via inclusion for different functionals is discussed. In conclusion, the most transparent linear quadratic case is considered. The statement is illustrated by examples.

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

The paper deals with an extremal routing problem with constraints. In the general formulation, it is assumed that the objects of visiting are any non-empty finite sets — megalopolises. The main applied problem considered in this study is the tool path optimization problem for CNC sheet-cutting machines, known as the Cutting Path Problem. This problem arises at the stage of developing control programs for CNC machines. Other applications are also possible. In particular, the results obtained in the chapter can be used in the problem of minimizing the radiation dose when dismantling a system of radiation-hazardous elements after accidents at nuclear power plants and in transport problems. Among tasks constraints, the precedence constraints are investigated. These constraints can be used to reduce computational complexity. As the main method, the study used broadly understood dynamic programming. The offered realization of the method takes into account the precedence constraints and the dependence of the objective functions on the task list. This dependence belongs to the class of very complex conditions that determine the route admissibility at each routing step, depending on the tasks already completed or, on the contrary, not yet completed. As applied to the Cutting Path Problem, the dependence of the objective function on the task list makes it possible to reduce thermal deformations of the material during cutting. The chapter provides a mathematical formalization of an extremal routing problem with additional constraints, a description of the method, and the exact algorithm obtained with its help. The order of task execution, the specific trajectory of the process, and the starting point are optimized.

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

The paper is devoted to elaboration of polyhedral techniques for solving two control problems for linear discrete-time systems with uncertainties under state constraints, namely, the terminal approach problem and the terminal evasion one. Such problems arise in systems with two controls, where the aim of the first is to steer the trajectory onto a given terminal set at a given instant without violating the state constraints, the aim of the other is opposite. It is assumed that the terminal set is a parallelepiped, the controls are bounded by parallelotope-valued constraints, and the state constraints are given in the form of so-called zones. We present techniques for solving both problems basing on polyhedral (parallelotope-valued or parallelepiped-valued) tubes. The techniques for solving the approach problem were proposed by the author earlier, but here additional properties of them are investigated. In particular, for the case without state constraints, guaranteed estimates are found for the trajectory that ensure that it is inside the tube. Convenient sufficient conditions are given to guarantee the obtaining of nondegenerate cross-sections during the calculations. For the evasion problem, a common solution scheme is considered, and then polyhedral techniques are proposed. The whole parametric families of external and internal polyhedral estimates for the solvability tubes for both problems are presented and compared. An illustrative example is given.

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

A zero-sum linear-convex differential game with a quality index that estimates a set of deviations of a motion trajectory at given instants of time from given target points is considered. A case when the saddle point condition in a small game, also known as Isaac's condition, does not hold, is studied. The game is formalized in classes of mixed control strategies of players. A numerical method for approximate computation of the game value and optimal strategies is elaborated. The method is based on the recurrent construction of upper convex hulls of auxiliary program functions. The results of numerical experiments in model examples are given.

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

We consider a system which consists of a circular cylinder subject to gravity interacting with a point vortex in a perfect fluid. In contrast to previous works, in this paper the circulation about the cylinder is assumed to be zero. The governing equations are Hamiltonian and admit evident integrals of motion: the horizontal and vertical components of the momentum; the latter is obviously non-autonomous. Using autonomous integral we reduce the order of the system by one degree of freedom in a case of zero circulation which early was not considered. Unlike nonzero circulation in the absence of point vortices when the cylinder moves inside a certain horizontal stripe it is shown that in the presence of vortices and with circulation equal to zero a vertical coordinate of the cylinder is unbounded decreasing. We then focus on the numerical study of dynamics of our system. In a case of zero circulation trajectories are noncompact. The different kinds of the scattering function of the vortex by cylinder were obtained. The form of these functions argues to chaotic behavior of the scattering which means that an additional analytical integral is absent.

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

The paper deals with the statement of a problem of dynamic interaction of arbitrary solid bodies and its test solutions in the context of discrete element modeling. For discretization we use description of bodies with arbitrary shapes, composed of rigidly bound spheres. The clumps were built with different characteristics, which allowed to estimate their influence on the process of clump construction and the smoothness of obtained surface. A system of equations of motion relative to global axes for a clump of spheres is presented. The forces of interaction between the spheres are determined based on the Hertz-Mindlin contact model with due account for viscous damping. A problem of interaction of two spheres was chosen as a test case. Spheres' trajectories composed of clumps of spheres were calculated. The results were compared with the results for the case of motion and interaction of spheres in one-particle approximation.

В работе рассматриваются динамические биматричные игры с интегральными показателями, дисконтированными на бесконечном интервале времени. Динамика системы задается дифференциальными уравнениями, описывающими изменение поведения игроков в зависимости от поступающих сигналов управления. Рассматривается задача построения равновесных траекторий в рамках минимаксного подхода, предложенного Н.Н. Красовским и А.И. Субботиным в теории дифференциальных игр. Используется конструкция динамического равновесия по Нэшу, которая развита в работах А.Ф. Клейменова. Для синтеза оптимальных стратегий управления применяется принцип максимума Л.С. Понтрягина в сочетании с методом характеристик для уравнений Гамильтона-Якоби. Получены аналитические формулы для кривых переключения оптимальных стратегий управления. Проведен анализ чувствительности равновесных решений в зависимости от параметра дисконтирования в интегральных функционалах выигрыша. Установлена асимптотическая сходимость равновесных траекторий по параметру дисконтирования к решению динамической биматричной игры со среднеинтегральными функционалами выигрыша, которые исследовались в работах В.И. Арнольда. Рассмотрено приложение полученных результатов к динамической модели инвестирования на финансовых рынках.

The paper is devoted to the analysis of dynamical bimatrix games with integral indices discounted on an infinite time interval. The system dynamics is described by differential equations in which players' behavior changes according to incoming control signals. For this game, a problem of construction of equilibrium trajectories is considered in the framework of minimax approach proposed by N.N. Krasovskii and A.I. Subbotin in the differential games theory. The game solution is based on the structure of dynamical Nash equilibrium developed in papers by A.F. Kleimenov. The maximum principle of L.S. Pontryagin in combination with the method of characteristics for Hamilton-Jacobi equations are applied for the synthesis of optimal control strategies. These methods provide analytical formulas for switching curves of optimal control strategies. The sensitivity analysis for equilibrium solutions is implemented with respect to the discount parameter in the integral payoff functional. It is shown that equilibrium trajectories in the problem with the discounted payoff functional asymptotically converge to the solution of a dynamical bimatrix game with average integral payoff functionals examined in papers by V.I. Arnold. Obtained results are applied to a dynamical model of investments on financial markets.