Все выпуски

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

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

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

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

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

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

На основе известных свойств функции вероятности протекания простой кубической решётки размера L=2 в приближении линейной связи порога протекания бесконечной решётки x_c и среднего значения x_cL конечной решётки введена нескейлинговая функция вероятности протекания для решётки размера L>2. Показано, что на пороге протекания нескейлинговые вероятности для всех ПК решёток одинаковы.
Компьютерные эксперименты на основе метода Монте-Карло согласуются с предлагаемой в работе теорией.

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

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

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

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