Открытый доступ Открытый доступ  Доступ закрыт Доступ предоставлен  Доступ закрыт Только для подписчиков

Том 64, № 12 (2024)

Обложка

Весь выпуск

Открытый доступ Открытый доступ
Доступ закрыт Доступ предоставлен
Доступ закрыт Только для подписчиков

ОБЩИЕ ЧИСЛЕННЫЕ МЕТОДЫ

ОЦЕНКА ОСТАТОЧНОГО ЧЛЕНА ПРИ СУММИРОВАНИИ НЕКОТОРЫХ ГИПЕРГЕОМЕТРИЧЕСКИХ РЯДОВ ГОРНА

Безродных С.И., Дунин-Барковская О.В.

Аннотация

Построены интегральные представления и асимптотические оценки остатков суммирования гипергеометрического ряда Аппеля 1 и родственного ему ряда 2, указанных в списке Горна гипергеометрических рядов двух переменных. Найденные формулы имеют приложение к разработке алгоритмов вычисления функции 1 с помощью формул аналитического продолжения во все пространство C2. Результаты могут найти приложение в задачах математической физики, и вычислительной теории функции, в том числе, при построении конформногоотображениясложныхмногоугольниковнаосновеинтегралаКристоффеля–Шварца. Библ.24.
Журнал вычислительной математики и математической физики. 2024;64(12):2229–2242
pages 2229–2242 views

БЫСТРОЕ ВЫЧИСЛЕНИЕ ИНТЕГРАЛЬНЫХ ОПЕРАТОРОВ ТИПА СВЕРТКИ В ЗАДАЧАХ ОЦЕНИВАНИЯ ОПЦИОНОВ ВМОДЕЛЯХ ЛЕВИ

Гречко А.С., Кудрявцев О.Е.

Аннотация

Построен приближенный алгоритм вычисления интегральных операторов типа свертки, возникающих при оценивании барьерных опционов в моделях Леви методом Винера–Хопфа. Дополнительно исследован вопрос возможности применения методов машинного обучения (искусственных нейронных сетей) к аппроксимации специального вида интегралов, являющихся ключевым элементом в конструкции приближенных формул для рассматриваемых интегральных операторов Винера–Хопфа. Основная идея заключается в разложении функции цены в ряд Фурье и трансформации контура интегрирования для каждого слагаемого ряда Фурье. В результате мы получаем набор типовых интегралов, которые зависят от факторов Винера–Хопфа, но не зависят от функции цены, при этом наиболее затратная с вычислительной точки зрения часть численного метода сводится к вычислению указанных интегралов. Поскольку вычислять их нужно только один раз, а не на каждой итерации, как это было в стандартных реализациях метода Винера–Хопфа, то это существенно ускорит вычисления. Более того, для вычисления типовых интегралов можно обучить нейросеть. Предлагаемый подход особенно эффективен для спектрально односторонних процессов Леви, для которых известны явные формулы факторизации Винера–Хопфа. В этом случае мы получаем удобные для вычислений формулы путем интегрирования вдоль разреза. Главным преимуществом включения нейросетей в вычислительную схему является возможность проводить вычисления на неравномерной сетке. Такой гибридный численный метод сможет успешно конкурировать с классическими методами вычислений сверток в подобных задачах с помощью быстрого преобразования Фурье. Вычислительные эксперименты показывают, что нейросети с одним скрытым слоем из 20 нейронов способны эффективно справляться с задачами аппроксимации рассматриваемых вспомогательных интегралов. Библ. 25. Фиг. 2. Табл. 1.
Журнал вычислительной математики и математической физики. 2024;64(12):2243–2261
pages 2243–2261 views

НЕВЫРОЖДЕННАЯ МАТРИЦА С ХОРОШО ОБУСЛОВЛЕННЫМ КОКВАДРАТОМ: КАК ПРИВЕСТИ ЕЕ К ДИАГОНАЛЬНОМУ ВИДУ ПОСРЕДСТВОМКОНГРУЭНЦИИ

Икрамов Х.Д., Назари А.М.

Аннотация

Существуют эффективные программы для приведения диагонализуемой матрицы к диагональному виду посредством преобразования подобия. Юнитоидные матрицы являются аналогами диагонализуемых матриц в теории конгруэнтных преобразований. Однако, если исключить эрмитовы и, более общо, нормальные матрицы, то нет пользующихся общим признанием программ для приведения юнитоидной матрицы к диагональному виду посредством конгруэнций. Предложен алгоритм, способный выполнять эту задачу для специального класса юнитоидных матриц, а именно невырожденных матриц, коквадраты которых хорошо обусловлены в смысле полной проблемы собственных значений. Приведены примеры, иллюстрирующие работу этого алгоритма. Библ. 1.
Журнал вычислительной математики и математической физики. 2024;64(12):2262–2269
pages 2262–2269 views

АПОСТЕРИОРНЫЕ ОЦЕНКИ ПОГРЕШНОСТИ ПРИБЛИЖЕННЫХ РЕШЕНИЙ ЭЛЛИПТИЧЕСКИХ КРАЕВЫХ ЗАДАЧ В ТЕРМИНАХ ЛОКАЛЬНЫХ НОРМ И ЦЕЛЕВЫХ ФУНКЦИОНАЛОВ

Музалевский А.В., Репин С.И., Фролов М.Е.

Аннотация

Получены функциональные соотношения, которые позволяют оценивать точность приближенных решений в терминах мер, существенно отличных от энергетических норм, которые обычно используются для этих целей. В частности, они применимы к локальным нормам и мерам, построенным с помощью специально построенных линейных функционалов. Потребность в таких инструментах контроля точности возникает, если имеется особый интерес к поведению решения в некоторой подобласти или к специальным свойствам решения. Показано, что апостериорные оценки функционального типа, которые ранее использовались для глобальных оценок, могут быть адаптированы и для решения этой задачи. Получены функциональные тождества и оценки, позволяющие оценивать погрешность любых конформных аппроксимаций в терминах широкого класса мер, включающих локальные нормы и проблемно-ориентированные функционалы. Теоретические результаты проверены в серии примеров, которые подтверждают эффективность предлагаемого метода. Библ. 13. Фиг. 10. Табл. 10.
Журнал вычислительной математики и математической физики. 2024;64(12):2270–2285
pages 2270–2285 views

ОБОБЩЕНИЯ СТАДИЙНОГО ПОРЯДКА МЕТОДОВ РУНГЕ–КУТТЫ

Скворцов Л.М.

Аннотация

Рассматривается применение методов Рунге–Кутты для решения жестких систем обыкновенных дифференциальных уравнений и дифференциально-алгебраических уравнений . При решении таких задач часто проявляется эффект снижения порядка, когда при заданной точности реальный порядок метода оказывается ниже классического порядка, что неизбежно приводит к повышению вычислительных затрат. Чтобы избежать снижения порядка, метод должен иметь достаточно высокий стадийный порядок. Однако методы, обеспечивающие наиболее удобную и эффективную реализацию, имеют низкий стадийный порядок. Поэтому актуальна задача построения методов, которые при низком стадийном порядке обладают свойствами методов более высокого стадийного порядка. Настоящая статья посвящена построению методов такого типа. Рассматриваются однократно диагонально-неявные, явные и обратные к явным методы. Приведены результаты решения тестовых задач. Библ. 44. Фиг. 3. Табл. 4.
Журнал вычислительной математики и математической физики. 2024;64(12):2286–2302
pages 2286–2302 views

КВАДРАТУРНЫЕ ФОРМУЛЫ ДЛЯ СИНГУЛЯРНЫХ ИНТЕГРАЛОВ, СОДЕРЖАЩИЕ ЗНАЧЕНИЯ ФУНКЦИИ И ЕЕ ПРОИЗВОДНЫХ

Хубежты Ш.С., Плиева Л.Ю.

Аннотация

Строятся квадратурные формулы для сингулярных интегралов на отрезке интегрирования [–1, 1] с определенными весовыми функциями
Журнал вычислительной математики и математической физики. 2024;64(12):2303–2311
pages 2303–2311 views

ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ

МЕТОД ВОЗМУЩЕНИЙ И РЕГУЛЯРИЗАЦИЯ ПРИНЦИПА ЛАГРАНЖА В НЕЛИНЕЙНЫХ ЗАДАЧАХ НА УСЛОВНЫЙ ЭКСТРЕМУМ

Сумин М.И.

Аннотация

Рассматривается регуляризация принципа Лагранжа (ПЛ) в недифференциальной форме в нелинейной (невыпуклой) задаче на условный экстремум с операторным ограничением-равенством в гильбертовом пространстве. Множество ее допустимых элементов принадлежит полному метрическому пространству, существование решения задачи априори не предполагается. Ограничение-равенство содержит аддитивно входящий в него параметр, что обеспечивает возможность применения для исследования задачи “нелинейного варианта” метода возмущений. Основное предназначение регуляризованного ПЛ—устойчивое генерирование обобщенных минимизирующих последовательностей (ОМП) в рассматриваемой нелинейной задаче. Его можно трактовать как ОМП-образующий (регуляризирующий) оператор, ставящий в соответствие каждому набору исходных данных задачи субминималь (минималь) ее отвечающего этому набору регулярного модифицированного функционала Лагранжа (МФЛ), двойственная переменная в котором генерируется в соответствии с процедурой стабилизации по Тихонову двойственной задачи. КонструкцияМФЛполностью определяется видом “нелинейных” субдифференциалов полунепрерывной снизу и, вообще говоря, невыпуклой функции значений как функции параметра задачи. В качестве таких субдифференциалов используются хорошо известные в негладком (нелинейном) анализе проксимальный субградиент и субдифференциал Фреше. Регуляризованный ПЛ “преодолевает” свойства некорректности классического аналога и может трактоваться как регуляризирующий алгоритм, составляя тем самым теоретическую основу для создания устойчивых методов практического решения нелинейных задач на условный экстремум. Библ. 32.
Журнал вычислительной математики и математической физики. 2024;64(12):2312–2331
pages 2312–2331 views

ОБЫКНОВЕННЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

О МИНИМАЛЬНОСТИ КВАДРАТИЧНОЙ ПОГРЕШНОСТИ РЕШЕНИЯ ПРЕОБРАЗОВАННЫХ К НАИЛУЧШЕМУ ПАРАМЕТРУ СИСТЕМ УРАВНЕНИЙ ПРИ МАЛЫХ ОДНОРОДНЫХ ВОЗМУЩЕНИЯХ

Кузнецов Е.Б., Леонов С.С.

Аннотация

В статье рассматривается решение систем нелинейных уравнений с одним скалярным параметром. Множеством решений подобных систем является кривая в пространстве неизвестных системы уравнений и параметра. Ее построение проводится, как правило, при помощи численных методов и сопряжено с многочисленными трудностями, возникающими вследствие наличия на кривой множества решений предельных и существенноособыхточек.Длянахождениятакихкривыхиспользуетсяметодпродолжениярешенияпопараметру и наилучшей параметризации, позволяющий свести решение к начальной задаче для системы дифференциальных уравнений продолжения решения. В данной работе исследуется устойчивость решения системы продолжения решения на вносимые в нее возмущения. Впервые полностью доказано сформулированное ранее утверждение о минимальности квадратичной ошибки решения системы продолжения решения при однородных малых возмущениях ее матрицы. Теоретические результаты проиллюстрированы на примере численного построения лемнискаты Бернулли. Библ. 10. Фиг. 2. Табл. 1.
Журнал вычислительной математики и математической физики. 2024;64(12):2332–2354
pages 2332–2354 views

УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ

МЕТОД УСРЕДНЕНИЯ В ЗАДАЧЕ ПОСТРОЕНИЯ АВТОКОЛЕБАТЕЛЬНЫХ РЕШЕНИЙ РАСПРЕДЕЛЕННЫХ КИНЕТИЧЕСКИХ СИСТЕМ

Кубышкин Е.П.

Аннотация

Построен метод усреднения для двухкомпонентных распределенных кинетических систем с малой диффузией в ограниченной одномерной области с условиями непроницаемости на границе. Построены преобразования рассматриваемой распределенной системы, позволяющие выделить одну “быструю” и счетное число “медленных” переменных. Доказаны теоремы о соответствии стационарных и периодических решений, а также инвариантных торов усредненных уравнений “медленных” переменных соответственно пространственно неоднородным периодическим решениям и инвариантным торам исходных уравнений аналогичного характера устойчивости. Предложены алгоритмы построения периодических решений (циклов) и инвариантных торов исходных уравнений в виде разложения по степеням малого параметра, обеспечивающих построениеасимптотическихформулуказанныхавтоколебательныхобъектов.Сформулированыусловиясходимости соответствующих разложений. Библ. 20.
Журнал вычислительной математики и математической физики. 2024;64(12):2355–2370
pages 2355–2370 views

МЕТОД МАЛОГО ПАРАМЕТРА В ТЕОРИИ УРАВНЕНИЙ ТИПА БЮРГЕРСА

Качалов В.И., Маслов Д.А.

Аннотация

Введенное Г. Бейтманом в 1915 г. и изученное Й. М. Бюргерсом в 1948 г. уравнение Бюргерса нашло широкое применениевмеханикежидкости,нелинейнойакустикеидругихобластяхприкладнойматематики.Подходы к его решению были самые разнообразные: асимптотические, численные, аналитические. В данной работе развивается аналитический метод решения уравнения типа Бюргерса в банаховом пространстве. А именно, после искусственного введения в уравнение малого параметра доказывается существование аналитического по этому параметру решения. При этом, рассматривается также и многомерный вариант уравнения Бюргерса. Библ. 16.
Журнал вычислительной математики и математической физики. 2024;64(12):2371–2377
pages 2371–2377 views

ОБ ОДНОЗНАЧНОСТИ ОПРЕДЕЛЕНИЯ СЕТОЧНОГО ФУНДАМЕНТАЛЬНОГО РЕШЕНИЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ И ВОЛНОВОГО УРАВНЕНИЯ В РАМКАХ ТЕОРИИ ДИСКРЕТНОГО ПОТЕНЦИАЛА

Степанова И.Э., Колотов И.И., Щепетилов А.В., Ягола А.Г., Левашов А.Н.

Аннотация

В работе рассматривается проблема однозначного определения фундаментального решения сеточного аналога волнового уравнения, а также уравнения теплопроводности в рамках теории дискретного потенциала. Сеточные фундаментальные решения конечно-разностных аналогов уравнений в частных производных позволяют решать прямые и обратные задачи по восстановлению источников волн и тепла в различных средах по разнородной и разноточной информации о соответствующих физических полях. В статье рассматриваются постановки с условиями Дирихле в трехмерном и четырехмерном декартовых пространствах. Библ. 16.
Журнал вычислительной математики и математической физики. 2024;64(12):2378–2389
pages 2378–2389 views

МАТЕМАТИЧЕСКАЯ ФИЗИКА

РАЗВИТИЕ МЕТОДА АДАПТИВНОЙ ИСКУССТВЕННОЙ ВЯЗКОСТИ ДЛЯ ГАЗОДИНАМИЧЕСКИХ РАСЧЕТОВ НА НЕРАВНОМЕРНЫХ РАЗНОСТНЫХ СЕТКАХ

Круковский А.Ю., Попов И.В., Гасилов В.А.

Аннотация

Метод адаптивной искусственной вязкости в настоящей работе обобщается с целью построения разностных схем газовой динамики, обеспечивающих высокое разрешение структуры течений как на равномерных, так и на неравномерных сетках. Рассматриваются разностные схемы, аппроксимирующие одномерную систему уравнений газовой динамики. Полученные в работе оценки величины адаптивной вязкости учитывают неоднородность распределения газодинамических величин в расчетной области и неравномерность разностной сетки. Построенные схемы с адаптивной искусственной вязкостью обладают свойствами однородности и консервативности. Апробация предложенных схем выполнена на модельных задачах, решения которых описывают различные гладкие газодинамические структуры, а также сильные и слабые разрывы. Продемонстрирована возможность получения высокоточных решений на расчетных сетках с существенной разницей геометрических размеров соседних разностных ячеек. Библ. 13. Фиг. 5.
Журнал вычислительной математики и математической физики. 2024;64(12):2390–2400
pages 2390–2400 views

МОДЕЛИРОВАНИЕ НЕЛИНЕЙНЫХ ВОЛНОВЫХ ПРОЦЕССОВ В СВЧ-ГЕНЕРАТОРЕ С МАГНИТНОЙ ИЗОЛЯЦИЕЙ

Поляков С.В., Тарасов Н.И., Кудряшова Т.А.

Аннотация

Рассмотрена актуальная задача моделирования нелинейных волновых процессов в СВЧ-генераторе с магнитной изоляцией. Для ее численного анализа предложена новая компьютерная модель, включающая уравнения Максвелла и уравнения движения релятивистских заряженных частиц, совместное их интегрирование методом сеток и методом облачных частиц, а также параллельную программную реализацию. В численных экспериментах получены пространственно-временные характеристики релятивистских электронных пучков и плазмы, а также параметры выходного излучения генератора. Анализ полученных результатов подтвердил корректность разработанного численного подхода. Библ. 17. Фиг. 7.
Журнал вычислительной математики и математической физики. 2024;64(12):2401–2410
pages 2401–2410 views

“СУПЕР-БЫСТРЫЙ” АЛГОРИТМ РЕШЕНИЯ ПРЯМОЙ ЗАДАЧИ РАССЕЯНИЯ ДЛЯ СИСТЕМЫ МАНАКОВА

Фрумин Л.Л., Чернявский А.Е., Белай О.В.

Аннотация

Рассматривается построение ускоренного алгоритма решения прямой задачи рассеяния для непрерывного спектра системы Манакова, ассоциированной с векторным нелинейным уравнением Шрёдингера модели Манакова. Численная постановка задачи приводит к проблеме быстрого расчета произведений полиномов, зависимых от спектрального параметра задачи. Для локализованных решений представлен так называемый “супер-быстрый” алгоритм решения прямой задачи рассеяния второго порядка точности, основанный на теореме о свертке и быстром преобразование Фурье, требующий для дискретной сетки размером асимптотически всего (︀ Log2 )︀ арифметических операций. Для ускорения расчета спектров коэффициентов отражения предложен и апробирован матричный вариант быстрого преобразования Фурье, когда коэффициенты ряда дискретного преобразования Фурье представляют собой некоммутирующие матрицы. Численное моделирование на примере точного решения системы Манакова (гиперболического секанса) подтвердило высокую скорость расчетов и второй порядок точности аппроксимации алгоритма. Библ. 24. Фиг. 1.
Журнал вычислительной математики и математической физики. 2024;64(12):2411–2419
pages 2411–2419 views

ИНФОРМАТИКА

ГИПОТЕЗА РАМСЕЯ О СОЦИАЛЬНОЙ СТРАТИФИКАЦИИ КАК ПРИНЦИП ОТБОРА ПО ФИШЕРУ

Парастаев Г.С., Шананин А.А.

Аннотация

Гипотеза Рамсея о социальной стратификации утверждает, что в популяции домашних хозяйств богатство концентрируется у наиболее бережливых агентов, которые дисконтируют потребительские расходы с наименьшим коэффициентом дисконтирования. Гипотезу Рамсея можно рассматривать как утверждение о справедливости принципа естественного отбора по Фишеру в популяции домашних хозяйств. В статье на основе гипотезы Дж. Дьюзенберри коэффициенты дисконтирования формируются в зависимости от распределения капитала между агентами. Поведение домашних хозяйств описывается моделями рационального репрезентативного потребителя рамсеевского типа. Для соответствующих задач оптимального управления построены решения в форме синтеза, которые использованы при моделировании динамики популяции домашних хозяйств. Доказаны теоремы для популяции домашних хозяйств, обосновывающие справедливость гипотезы Рамсея. Исследовано влияние потребительского кредита на социальную стратификацию домашних хозяйств. Библ. 28. Фиг. 2.
Журнал вычислительной математики и математической физики. 2024;64(12):2420–2448
pages 2420–2448 views