Все выпуски

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

The paper is devoted to a method of solving of stationary filtration problems of non-compressible fluid which follows the nonlinear multi-valued anisotropic law of filtration with limiting gradient. This problem mathematically is formulated in the form of variational inequality of the second kind in Hilbert space with inversely strongly monotone operator. The functional occurring in this variational inequality is a sum of several lower semi-continuous convex proper functionals. For solving the considered variational inequality the splitting method is offered.

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

In this paper, we consider the stable reconstruction problem of the unknown input of a distributed system of second order by results of inaccurate measurements of its solution. The content of the problem considered is as follows. We consider a distributed equation of second order. The solution of the equation depends on the input varying in the time. The input, as well as the solution, is not given in advance. At discrete times the solution of the equation is measured. These measurements are not accurate in general. It is required to design an algorithm for approximate reconstruction of the input that has dynamical and stability properties. The dynamical property means that the current values of approximations of the input are produced on-line, and the stability property means that the approximations are arbitrarily accurate for a sufficient accuracy of measurements. The problem refers to the class of inverse problems. The algorithm presented in the paper is based on the constructions of a stable dynamical inversion and on the combination of the methods of ill-posed problems and positional control theory.

В данной работе рассматривается система Каупа–Буссинеска с самосогласованным источником. Показано, что система Каупа–Буссинеска с самосогласованным источником может быть проинтегрирована методом обратной задачи рассеяния. Для решения рассматриваемой задачи используются прямая и обратная задачи рассеяния уравнения Штурма–Лиувилля с потенциалом, зависящим от энергии. Определена временная эволюция данных рассеяния для уравнения Штурма–Лиувилля с энергозависимыми потенциалами, связанными с решением системы Каупа–Буссинеска с самосогласованным источником. Полученные равенства полностью определяют данные рассеяния при любом $t$, что позволяет применить метод обратной задачи рассеяния для решения задачи Коши для системы Каупа–Буссинеска с самосогласованным источником.

In this study we consider the Kaup–Boussinesq system with a self-consistent source. We show that the Kaup–Boussinesq system with a self-consistent source can be integrated by the method of inverse scattering theory. For a solving the problem under consideration, we use the direct and inverse scattering problem of the Sturm–Liouville equation with an energy-dependent potential. The time evolution of the scattering data for the Sturm–Liouville equation with an energy-dependent potentials associated with the solution of the Kaup–Boussinesq system with a self-consistent source is determined. The obtained equalities completely determine the scattering data for any $t$, which makes it possible to apply the method of the inverse scattering problem to solve the Cauchy problem for the Kaup–Boussinesq system with a self-consistent source.

В статье исследуются прямая и обратная задачи для уравнений субдиффузии с участием дробной производной в смысле Хильфера. В качестве эллиптической части уравнения взят произвольный положительный самосопряженный оператор $A$. В частности, в качестве оператора $A$ можно взять оператор Лапласа с условием Дирихле. Сначала доказано существование и единственность решения прямой задачи. Затем с помощью представления решения прямой задачи доказывается существование и единственность обратной задачи нахождения правой части уравнения, зависящей только от пространственной переменной.

The article studies direct and inverse problems for subdiffusion equations involving a Hilfer fractional derivative. An arbitrary positive self-adjoint operator $A$ is taken as the elliptic part of the equation. In particular, as the operator $A$ we can take the Laplace operator with the Dirichlet condition. First, the existence and uniqueness of a solution to the direct problem is proven. Then, using the representation of the solution to the direct problem, the existence and uniqueness of the inverse problem of finding the right-hand side of the equation, which depends only on the spatial variable, is proved.

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

A patrolling game with two players, a patroller and an attacker, is considered in the paper. The aim of the former is to protect an object from intruders and catch the attacker. The aim of the latter is to cause damage to the protected object without being caught. Cellular base stations are viewed as protected objects. A game-theoretic model is constructed to find an initial distribution of players on base stations. When the transition matrix of players among the stations is known, an optimal strategy of players and the value of the game are calculated. An inverse problem of searching for optimal transition matrices with known initial distribution of players is studied. The Nash equilibrium with the attacker making two attacks is found for the considered problem.

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

This work studies direct initial boundary value and inverse coefficient determination problems for a one-dimensional partial differential equation with multi-term orders fractional Riemann–Liouville derivatives. The unique solvability of the direct problem is investigated and a priori estimates for its solution are obtained in weighted spaces, which will be used for studying the inverse problem. Then, the inverse problem is equivalently reduced to a nonlinear integral equation. The fixed-point principle is used to prove the unique solvability of this equation.

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

The stable sequential Pontryagin maximum principle or, in other words, the regularized Pontryagin maximum principle in iterative form is formulated for the optimal control problem of a linear parabolic equation with distributed, initial and boundary controls and operator semiphase equality constraint. The main difference between it and the classical Pontryagin maximum principle is that, firstly, it is formulated in terms of minimizing sequences, secondly, the iterative process occurs in dual space, and thirdly, it is resistant to error of raw data and gives a minimizing approximate solution in the sense of J. Warga. So it is a regularizing algorithm. The proof of the regularized Pontryagin maximum principle in iterative form is based on the dual regularization methods and iterative dual regularization. The results of model calculations of the concrete optimal control problem illustrating the work of the algorithm based on the regularized iterative Pontryagin maximum principle are presented. The problem of finding a control triple with minimal norm under a given equality constraint at the final instant of time or, in other words, the inverse final observation problem of finding a normal solution is used as a concrete model optimal control problem.

Рассматривается задача построения вершинного описания выпуклого полиэдра, заданного как множество решений некоторой системы линейных неравенств, коэффициенты которой являются алгебраическими числами. Обратная задача эквивалентна (двойственна) исходной. Предлагаются программные реализации нескольких модификаций хорошо известного метода двойного описания (метода Моцкина-Бургера), решающего поставленную задачу. Рассматривается два случая: 1) элементы системы неравенств - произвольные алгебраические числа, при этом каждое такое число задается минимальным многочленом и локализующим интервалом; 2) элементы системы неравенств принадлежат заданному конечному расширению ${\mathbb Q} (\alpha)$ поля ${\mathbb Q}$, при этом для $\alpha$ задаются минимальный многочлен и локализующий интервал, а все элементы исходной системы, конечные и промежуточные результаты представлены как многочлены от $\alpha$. Как и ожидалось, программная реализация для второго варианта значительно превосходит реализацию для первого варианта по производительности. Для большего ускорения во втором случае предлагается использовать булевы матрицы вместо матриц невязок. Результаты вычислительного эксперимента показывают, что программные реализации вполне пригодны для решения задач умеренных размеров.

We consider the problem of constructing the dual representation of a convex polyhedron defined as a set of solutions to a system of linear inequalities with coefficients which are algebraic numbers. The inverse problem is equivalent (dual) to the initial problem. We propose program implementations of several variations of the well-known double description method (Motzkin-Burger method) solving this problem. The following two cases are considered: 1) the elements of the system of inequalities are arbitrary algebraic numbers, and each such number is represented by its minimal polynomial and a localizing interval; 2) the elements of the system belong to a given extension ${\mathbb Q} (\alpha)$ of ${\mathbb Q}$, and the minimal polynomial and the localizing interval are given only for $\alpha$, all elements of the system, intermediate and final results are represented as polynomials of $\alpha$. As expected, the program implementation for the second case significantly outperforms the implementation for the first one in terms of speed. In the second case, for greater acceleration, we suggest using a Boolean matrix instead of the discrepancy matrix. The results of a computational experiment show that the program is quite suitable for solving medium-scale problems.

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

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.

Исследуется обратная задача определения многомерного ядра интегрального члена, зависящего от временной переменной $t$ и $ (n-1)$-мерной пространственной переменной $x'=\left(x_1,\ldots, x_ {n-1}\right)$ из $n$-мерного уравнения теплопроводности с переменным коэффициентом теплопроводности. Прямую задачу представляет задача Коши для этого уравнения. Интегральный член имеет вид свертки по времени ядра и решения прямой задачи. Дополнительное условие для решения обратной задачи задается решение прямой задачи на гиперплоскости $x_n = 0.$ В начале изучаются свойства решения прямой задачи. Для этого эта задача сводится к решению интегрального уравнения второго порядка вольтерровского типа и к нему применяется метод последовательных приближений. Далее поставленная обратная задача приводится к двум вспомогательным задачам, дополнительное условие второй из них содержит неизвестное ядро вне интеграла. Затем вспомогательные задачи заменяются эквивалентной замкнутой системой интегральных уравнений вольтерровского типа относительно неизвестных функций. Применяя метод сжатых отображений к этой системе в классе гёльдеровских функций доказываем основной результат статьи, который является теоремой локального существования и единственности решения обратной задачи.

The inverse problem of determining a multidimensional kernel of an integral term depending on a time variable $t$ and $ (n-1)$-dimensional spatial variable $x'=\left(x_1,\ldots, x_ {n-1}\right)$ in the $n$-dimensional heat equation with a variable coefficient of thermal conductivity is investigated. The direct problem is the Cauchy problem for this equation. The integral term has the time convolution form of kernel and direct problem solution. As additional information for solving the inverse problem, the solution of the direct problem on the hyperplane $x_n = 0$ is given. At the beginning, the properties of the solution to the direct problem are studied. For this, the problem is reduced to solving an integral equation of the second kind of Volterra-type and the method of successive approximations is applied to it. Further the stated inverse problem is reduced to two auxiliary problems, in the second one of them an unknown kernel is included in an additional condition outside integral. Then the auxiliary problems are replaced by an equivalent closed system of Volterra-type integral equations with respect to unknown functions. Applying the method of contraction mappings to this system in the Hölder class of functions, we prove the main result of the article, which is a local existence and uniqueness theorem of the inverse problem solution.