Все выпуски

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

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

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

В работе рассматривается задача оптимального управления одномерным процессом, заданным стохастическим дифференциальным уравнением, в котором управление воздействует как на коэффициент сноса, так и на коэффициент диффузии, при этом диффузионная составляющая линейна по управлению $$dx(t) = b(t,x(t),u(t))dt +\sigma(t,x(t))u(t)dW(t),\qquad x(0) = x_0.$$ Здесь $x(t)$ - фазовая координата, $u(t)$ - управляющая функция, $W(t)$ - винеровский процесс. Доказана теорема, которая предоставляет структуру решения рассматриваемого уравнения в виде суперпозиции функций $x(t)=Φ(t,u(t)W(t)+y(t))$, в котором $Φ(t,v)$ - известная функция, полностью определяющаяся коэффициентом $σ(t,x)$, и не зависит от управления, а $y(t)$ - решение потраекторно-детерминированного дифференциального уравнения с мерой вида

$$dy(t) = B(t,y(t),u(t))dt - W(t)du(t).$$

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

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

Для динамической системы, управляемой в условиях помех, рассматривается задача оптимизации гарантированного результата. Особенностью задачи является наличие функциональных ограничений на помехи, при которых свойство замкнутости множества допустимых помех относительно операции «склейки» двух его элементов, вообще говоря, отсутствует. Это обстоятельство препятствует непосредственному применению методов теории дифференциальных игр для исследования задачи и тем самым приводит к необходимости их походящей модификации. В работе предложено новое понятие неупреждающей стратегии управления (квазистратегии). Доказано, что соответствующий функционал оптимального гарантированного результата удовлетворяет принципу динамического программирования. Как следствие, установлены так называемые свойства $u$- и $v$-стабильности этого функционала, которые в дальнейшем позволят построить конструктивное решение задачи в позиционных стратегиях.

Изучаются свойства дискретной вариационной задачи динамической аппроксимации в комплексном евклидовом (L + 1)-мерном пространстве E. Она обобщает известные задачи среднеквадратической полиномиальной аппроксимации функций, заданных своими отсчетами в конечном интервале. В рассматриваемой задаче аппроксимация последовательности y = {y_i}_L⁰ отсчетов функции y(t) ∈ L²[0, T], T = Lh на сетке I_h осуществляется решениями однородных линейных дифференциальных или разностных уравнений заданного порядка n с постоянными, но, возможно, неизвестными коэффициентами. Тем самым показано, что в последнем случае задача аппроксимации включает в себя и задачу идентификации. Анализ ее особенностей - основная тема статьи. Ставится задача нахождения вектора коэффициентов  разностного уравнения Σ_n⁰ ŷ_i+k α_i = 0, где k = 0,L − n. Оптимизируются коэффициенты  и начальные условия переходного процесса y этого уравнения. Цель оптимизации - наилучшая аппроксимация исследуемого динамического процесса y ∈ E. Критерий аппроксимации  минимум величины ||y − ŷ||²_E. Показано, что изучаемая вариационная задача сводится к задачам проектирования в E вектора y на ядра разностных операторов с неизвестными коэффициентами α ∈ ω ⊂ S ⊂ Eⁿ⁺¹. Здесь α - направление, S - сфера или гиперплоскость. Показана связь изучаемой задачи с задачами дискретизации и идентифицируемости. Тогда координаты вектора y ∈ E есть точное решение дифференциального уравнения на сетке I_h и y = ŷ. Дано сравнение изучаемой задачи вариационной идентификации с алгебраическими методами идентификации. Показано, что ортогональные дополнения к ядрам разностных операторов всегда имеют теплицев базис. Это приводит к быстрым проекционным алгоритмам вычислений. Показано, что задача нахождения оптимального вектора α сводится к задаче безусловной минимизации функционала идентификации, зависящего от направления  в Eⁿ⁺¹. Предложена итерационная процедура его минимизации на сфере с широкой областью и высокой скоростью сходимости. Изучаемую вариационную задачу можно применять при математическом моделировании в управлении и научных исследованиях. При этом на конечных интервалах может использоваться, в частности, возможность кусочно-линейной динамической аппроксимации сложных динамических процессов разностными и дифференциальными уравнениями указанного типа.

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

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

В настоящей работе проведено исследование модели деформаций системы из $n$ стилтьесовских струн, расположенных вдоль геометрического графа-звезды, с нелинейным условием в узле. Соответствующая граничная задача имеет вид $$ \left\{\begin{array}{lll} -\left(p_iu_i^\prime\right)(x)+\displaystyle{\int_{0}^{x}}u_i\,dQ_i=F_i(x)-F_i(+0)-(p_iu_i')(+0),\quad i=1,2, \ldots, n,\\ \sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0),\\ u_1(0)=u_2(0)=\ldots=u_n(0)=u(0),\\ (p_iu_i')(l_i-0)+u_i(l_i)\Delta Q_i(l_i)=\Delta F_i(l_i),\quad i=1,2,\ldots, n. \end{array} \right. $$ Здесь функции $u_i(x)$ определяют деформации каждой из струн; $F_i(x)$ описывают распределение внешней нагрузки; $p_i(x)$ характеризуют упругость струн; $Q_i(x)$ описывают упругую реакцию внешней среды. Скачок $\Delta F_i(l_i)$ равняется сосредоточенной в точке $l_i$ внешней силе; скачок $\Delta Q_i(l_i)$ совпадает с жесткостью упругой опоры (пружины), прикрепленной к точке $l_i$. Условие $\sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0)$ возникает за счет наличия в узле ограничителя, представленного отрезком $[-m,m]$, на перемещение струн под воздействием внешней нагрузки, то есть предполагается, что $|u(0)|\leq m$. Здесь через $N_{[-m,m]}u(0)$ обозначен нормальный конус к $[-m,m]$ в точке $u(0)$. В работе проведен вариационный вывод модели; доказаны теоремы существования и единственности решения; проанализированы критические нагрузки, при которых происходит соприкосновение струн с ограничителем; приведена явная формула представления решения.