Все выпуски

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

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

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

В ограниченной по переменной $z$ области, имеющей слабо горизонтальную неоднородность, исследуется задача определения сверточного ядра $k(t,x)$, $t>0$, $x\in {\Bbb R}$, входящего в гиперболическое интегро-дифференциальное уравнение второго порядка. Предполагается, что это ядро слабо зависит от переменной $x$ и разлагается в степенной ряд по степеням малого параметра $\varepsilon$. Построен метод нахождения первых двух коэффициентов $k_{0}(t)$, $k_{1}(t)$ этого разложения по заданным первым двум моментам по переменной $x$ решения прямой задачи при $z=0$.

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

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

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

На примере известной задачи о прокладке трассы изучаются возможности численного решения сосредоточенных задач оптимального управления методом параметризации управления с помощью линейной комбинации $\mu$ функций Гаусса. Напомним, что функция Гаусса (называемая также квадратичной экспонентой) - это функция вида $\varphi(x)=\dfrac{1}{\sigma\sqrt{2\pi}}\exp\left[-\dfrac{(x-m)^2}{2\sigma^2}\right]$. Основу метода составляет сведение исходной бесконечномерной задачи оптимизации к конечномерной задаче минимизации целевого функционала по параметрам аппроксимации управления с последующим применением численных методов конечномерной оптимизации. Данная статья опирается на исследование, проведенное автором ранее и касавшееся возможностей аппроксимации функций одного переменного на конечном отрезке линейной комбинацией функций Гаусса, и является его непосредственным продолжением. Прежде всего, мы доказываем утверждение об аппроксимации на любом конечном отрезке материнского вейвлета «мексиканская шляпа» линейной комбинацией двух квадратичных экспонент. Отсюда получаем теоретическое обоснование возможности эффективной аппроксимации функций одного переменного на любом конечном отрезке линейными комбинациями функций Гаусса. После этого мы проводим сравнение качества аппроксимации указанного вида с аппроксимацией по Котельникову на базе численных экспериментов. Затем приводится постановка задачи о прокладке трассы, а также результаты ее численного решения при различных способах параметризации управления, наглядно демонстрирующие преимущества предлагаемого способа, в частности устойчивость численного решения к погрешности вычисления параметров аппроксимации оптимального управления даже при использовании малого количества этих параметров.

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

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