Все выпуски

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

The paper investigates the solvability of an inverse boundary-value problem with an unknown coefficient and the right-hand side, depending on the time variable, for the linearized Benney-Luke equation with non-self-adjoint boundary and additional integral conditions. The problem is considered in a rectangular domain. A definition of the classical solution of the problem is given. First, we consider an auxiliary inverse boundary-value problem and prove its equivalence (in a certain sense) to the original problem. To investigate the auxiliary inverse boundary-value problem, the method of separation of variables is used. By applying the formal scheme of the variable separation method, the solution of the direct boundary problem (for a given unknown function) is reduced to solving the problem with unknown coefficients. Then, the solution of the problem is reduced to solving a certain countable system of integro-differential equations for the unknown coefficients. In turn, the latter system of relatively unknown coefficients is written as a single integro-differential equation for the desired solution. Next, using the corresponding additional conditions of the inverse auxiliary boundary-value problem, to determine the unknown functions, we obtain a system of two nonlinear integral equations. Thus, the solution of an auxiliary inverse boundary-value problem is reduced to a system of three nonlinear integro-differential equations with respect to unknown functions. A special type of Banach space is constructed. Further, in a ball from a constructed Banach space, with the help of contracted mappings, we prove the solvability of a system of nonlinear integro-differential equations, which is also the unique solution to the auxiliary inverse boundary-value problem. Finally, using the equivalence of these problems the existence and uniqueness of the classical solution of the original problem are proved.

Работа посвящена исследованию второй начально-краевой задачи для дифференциального уравнения третьего порядка псевдопараболического типа с переменными коэффициентами в многомерной области с произвольной границей. Рассматриваемое многомерное псевдопараболическое уравнение сводится к интегро-дифференциальному уравнению с малым параметром и для полученного уравнения строится локально-одномерная разностная схема А.А. Самарского. С помощью принципа максимума получена априорная оценка решения локально-одномерной разностной схемы в равномерной метрике в норме $C$. Доказаны устойчивость и сходимость локально-одномерной разностной схемы.

The work is devoted to the study of the second initial-boundary value problem for a general-form third-order differential equation of pseudoparabolic type with variable coefficients in a multidimensional domain with an arbitrary boundary. In this paper, a multidimensional pseudoparabolic equation is reduced to an integro-differential equation with a small parameter, and a locally one-dimensional difference scheme by A.A. Samarskii is used. Using the maximum principle, an a priori estimate is obtained for the solution of a locally one-dimensional difference scheme in the uniform metric in the $C$ norm. The stability and convergence of the locally one-dimensional difference scheme are proved.

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

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.

Рассматриваются вопросы, связанные с решением аддитивной задачи последовательного обхода множеств с ограничениями предшествования и функциями стоимости, допускающими зависимость от списка заданий. В качестве базового метода используется широко понимаемое динамическое программирование (ДП), дополняемое в случае задач ощутимой размерности декомпозициями семейства заданий и преобразованием параметров исходной задачи. Возможные применения связаны, в частности, с задачей управления инструментом при фигурной листовой резке деталей на машинах с ЧПУ. В этой задаче важным обстоятельством является учет условий предшествования, имеющих, в частности, следующий смысл: в случае детали с отверстиями резка каждого из внутренних контуров (отвечающих отверстиям) должна предшествовать резке внешнего контура. Сам критерий качества в данной задаче, как правило, является аддитивным. Другой тип ограничений касается избежания термических деформаций деталей. При использовании подхода с применением штрафов за нарушение условий, связанных с эффективным отводом тепла при выполнении врезки, возникают функции стоимости, допускающие зависимость от списка заданий, выполненных на текущий момент времени. Заметим, что в другой прикладной задаче, а именно в задаче о демонтаже радиационно опасных объектов, возникают функции стоимости с зависимостью от списка заданий, не выполненных на данный момент (а, следовательно, касающихся недемонтированных объектов). В итоге мы приходим к очень общей задаче с ограничениями предшествования и функциями стоимости с зависимостью от списка заданий. Применяемая в случае ощутимой размерности декомпозиция с последующей реализацией ДП требует, с одной стороны, разработки методов кластеризации, а, с другой, построения адекватной конструкции распределения глобальных условий предшествования по кластерам. В теоретической части работы обсуждается случай двух кластеров, который позволяет охватить единой схемой целый ряд практически интересных задач диапазонного (в смысле размерности) типа. Указан алгоритм построения композиционного решения, включающий этап обучения кластеризации на основе жадного алгоритма. Данный «композиционный» алгоритм реализован на ПЭВМ; проведен вычислительный эксперимент.

Issues related to solving the additive problem of sequential traversal of sets with precedence restrictions and cost functions that allow dependence on the list of tasks are considered. The basic method is a broadly understood dynamic programming (DP), supplemented in the case of problems of appreciable dimension by decompositions of the family of tasks and transformation of the parameters of the original problem. Possible applications are related, in particular, to the problem of tool control in figured sheet cutting of parts on CNC machines. In this problem, an important circumstance is taking into account the precedence conditions, which have, in particular, the following meaning: in the case of a part with holes, cutting of each of the internal contours (corresponding to the holes) should precede cutting of the external contour. The quality criterion itself in this problem, as a rule, is additive. Another type of constraints concerns avoiding thermal deformations of parts. When using the approach with penalties for violating the conditions associated with effective heat dissipation during cutting, cost functions arise that allow dependence on the list of tasks completed to date. Note that in another applied problem, namely, in the problem of dismantling radiation hazardous objects, cost functions arise with dependence on the list of tasks that have not been completed at the moment (and, consequently, concern the objects that have not been dismantled). As a result, we arrive at a very general problem with precedence constraints and cost functions with dependence on the list of tasks. The decomposition applied in the case of a noticeable dimensionality with subsequent implementation of the DP requires, on the one hand, the development of clustering methods, and, on the other, the construction of an adequate structure for distributing global precedence conditions among clusters. In the theoretical part of the work, the case of two clusters is discussed, which makes it possible to cover with a single scheme a number of practically interesting problems of a range (in terms of dimensionality) type. An algorithm for constructing a composite solution is indicated, including a stage of clustering training based on a greedy algorithm. This “composite” algorithm is implemented on a PC; a computational experiment was carried out.

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

Nowadays, the problem of creating an optimal safe schedule for arrival of aircraft coming in several flows to a checkpoint, where these flows join into one, is very important for air-traffic management. Safety of the resultant queue is present if there is a safe interval between neighbor arrivals to the merge point. Change of an arrival instant of an aircraft is provided by changing its velocity and/or usage of fragments of the air-routes scheme, which elongate or shorten the aircraft path. Optimality of the resultant queue is considered from the point of some additional demands: minimization of the deviation of the actual aircraft arrival instant from the nominal one, minimization of order changes in the resultant queue in comparison with the original one, minimization of fuel expenditures, etc. The optimality criterion to be minimized, which reflects these demands, is often taken as a sum of penalties for deviations of the assigned arrival instants from the nominal ones. Each individual penalty is considered in almost all papers as either the absolute value of the difference between the assigned and nominal arrival instants or a similar function with asymmetric branches (which punishes delays and accelerations of an aircraft in different ways). The problem can be divided into two subproblems: one is a search for an optimal order of aircraft in the resultant queue, and the other is a search for optimal arrival instants for a given order. The second problem is quite simple since it can be formalized in the framework of linear programming and solved quite efficiently. However, the first one is very difficult and now is solved by various methods. The paper suggests sufficient conditions for the problem, which guarantee that the order of the optimal assigned instants is the same as the order of the nominal ones and, therefore, exclude the first subproblem.

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

The paper provides an overview of the problems that lead to a necessity for analyzing models of linear and nonlinear dynamic systems in the form of stochastic differential equations with random delays of various types as well as some well-known methods for solving these problems. In addition, the author proposes some new approaches to the approximate analysis of linear and nonlinear stochastic dynamic systems. Changes of delays in these systems are governed by discrete Markov chains with continuous time. The proposed techniques for the analysis of systems are based on a combination of the classical steps method, an extension of the state space of a stochastic system under examination, and the method of statistical modeling (Monte Carlo). In this case the techniques allow to simplify the task and to transfer the source equations to systems of stochastic differential equations without delay. Moreover, for the case of linear systems the author has obtained a closed sequence of systems with increasing dimensions of ordinary differential equations satisfied by the functions of conditional expectations and covariances for the state vector. The above scheme is demonstrated by the example of a second-order stochastic system. Changes of the delay in this system are controlled by the Markov chain with five states. All calculations and graphics were performed in the environment of the mathematical package Mathematica by means of a program written in the source language of the package.

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

We consider the nonlinear controlled system in a finite-dimensional Euclidean space in a finite time interval. We study the problem of a system approaching a given compact set in finite time. Approximate solution of the approaching problem is discussed. The approach used to construct an approximate solution is based on constructions based on the notion of a set of solvability of the approaching problem. The concept of correcting control with and without additional operating influences is introduced. We propose a scheme of approximate backward construction of the solvability set, as well as the scheme of control software, which allows finding approximately a solution to the approaching problem. In it, the operating influence breaks down into “main” and “correcting”. An estimate of the deviation of the operated system from the target set at the final moment is constructed and it is shown that the use of additional correcting control in the control process can essentially improve the result of control.