Введение 3
Постановка задачи 5
Обзор литературы 6
Глава 1. Методы Рунге — Кутты — Чебышёва 8
1.1 Формулировка методов 8
1.2 Устойчивость 9
Глава 2. Методы RKC для ДУЗА 12
2.1 Использование ранее найденных значений K 13
2.2 Линейная интерполяция 14
Глава 3. Глава 3. Устойчивость методов RKC для ДУЗА 15
3.1 Устойчивость решений тестового уравнения 15
3.2 Анализ численной устойчивости 16
3.3 Методы, использующие ранее найденные K 18
3.4 Методы с интерполяцией 20
Глава 4. «Оптимальное» демпфирование 24
4.1 Влияние демпфирования на область устойчивости 24
4.2 «Оптимальные» коэффициенты демпфирования 27
Глава 5. Проверка областей устойчивости 30
5.1 Реализация программы 30
5.2 Проверка области устойчивости метода RKC с интерполяцией 31
5.3 Проверка области устойчивости метода RKC, использующего
ранее найденные K 35
Выводы 38
Заключение 39
Список литературы 40
Приложение 42

Купить

В связи с обширной спецификой применения дифференциальных уравнений с запаздывающим аргументом (ДУЗА) исследование численных методов для их решения крайне важно. Многочисленные процессы, которые основаны на передаче энергии, массы и информации, связаны с проблемой запаздывания. Его причиной могут являться различные ситуации: ограничение скорости распространения взаимодействия, наличие инерционности некоторых элементов и многие другие. Таким образом, системы запаздывающего типа неизбежно появляются при построении математических моделей во многих научных областях. Особенно часто они используются биологии, экологии и медицине [1], но также находят широкое применение в математической физике, химии, экономике, физике и инженерии [2, 3], причём в последних нередко встречаются уравнения в частных производных, включающие в себя запаздывание. В частности, нам интересны параболические уравнения такого вида (см. напр. [4-6]).
Для систем обыкновенных дифференциальных уравнений (ОДУ), собственные числа которых сильно вытянуты вдоль действительной отрицательной полуоси и имеют небольшие мнимые части, которые в частности получаются после дискретизации по пространству параболических уравнений, существуют явные стабилизированные методы Рунге —Кутты, которые по аналитической форме их функции устойчивости называются методами Рунге —Кутты —Чебышёва. Основной мотивацией настоящей работы послужило исследование Ё. Комори, А. С. Еремина и К. Бёрриджа [7], в котором рассматривается устойчивость и применимость чебышёвских методов для стохастических уравнений с запаздывающим шумом. Однако в этой работе не рассматривается запаздывание в детерминистической части. Первым шагом к исследованию свойств методов для решения стохастических уравнений с запаздываниями как в шуме, так и в детерминистической части, должен послужить анализ методов на полностью детерминистических ДУЗА.
К сожалению, если такой анализ и проводился, его результаты, видимо, не были опубликованы. Настоящая работа восполняет этот пробел.
В работе сначала приводятся методы Рунге — Кутты — Чебышёва для ОДУ, а во второй главе рассматривается два варианта их расширения для интегрирования ДУЗА с постоянным запаздыванием. В третьей главе про-водится анализ устойчивости и приводятся области устойчивости изучаемых методов при действительных коэффициентах тестового уравнения.
Поскольку области устойчивости обычных методов Рунге — Кутты — Чебышёва при некоторых значениях коэффициентов тестового ОДУ сужаются до одной точки, что, как будет показано, сохраняется и для ДУЗА, в работе также рассматриваются модифицированные методы, получаемые введением так называемого демпфирования. Области таких методов чуть короче, но избавлены от упомянутого недостатка. Оказывается, что для ДУЗА можно естественным образом определить понятие «оптимального» коэффициента демпфирования, аналогично тому, как это сделано в [7] для стохастических ДУЗА. Этому вопросу посвящена четвёртая глава работы.
Для проверки построенных областей в пятой главе приводится численное тестирование методов.
Купить