Все выпуски

Работа посвящена изучению наилучших равномерных рациональных приближений (НРРП) непрерывных функций на компактных, в том числе конечных, подмножествах числовой оси $\mathbb{R}$. Показано, что НРРП на конечном множестве существует не всегда. Более подробно изучен алгоритм Гельмута Вернера поиска НРРП вида $P_m/Q_n = \sum\limits_{i=0}^m a_i x^i \big/ \sum\limits_{j=0}^n b_j x^j$ для функций на множестве из $N=m+n+2$ точек $x_1<\ldots<x_N$. Этот алгоритм может использоваться в алгоритме Ремеза поиска НРРП на отрезке. При работе алгоритма Вернера вычисляется $(n+1)$ вещественное собственное значение $h_1,\ldots,h_{n+1}$ для пучка матриц $A-hB$, где $A$ и $B$ - некоторые симметричные матрицы. Каждому собственному значению сопоставляется своя рациональная дробь вида $P_m/Q_n$, являющаяся кандидатом на наилучшее приближение. Поскольку не более одной из этих дробей свободны от полюсов на отрезке $[x_1, x_N]$, то возникает задача отыскания того собственного значения, которому соответствует рациональная дробь без полюсов. В работе показано, что если $m=0$, все значения $f(x_1),-f(x_2),\ldots,(-1)^{n+2} f(x_{n+2})$ различны и НРРП положительно (отрицательно) во всех точках $x_1,\ldots,x_{n+2}$, то это собственное значение занимает $[(n+2)/2]$-е ($[(n+3)/2]$-е) место по величине. Приведены три численных примера, иллюстрирующих это утверждение.

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

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

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

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

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

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

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

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

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