📄Работа №133046

Тема: Разработка экономичных схем интегрирования структурно разделённых систем обыкновенных дифференциальных уравнений

Характеристики работы

Тип работы Дипломные работы, ВКР
Математика
Предмет Математика
📄
Объем: 47 листов
📅
Год: 2018
👁️
Просмотров: 187
Не подходит эта работа?
Закажите новую по вашим требованиям
Узнать цену на написание
ℹ️ Настоящий учебно-методический информационный материал размещён в ознакомительных и исследовательских целях и представляет собой пример учебного исследования. Не является готовым научным трудом и требует самостоятельной переработки.

📋 Содержание

Введение 4
Постановка задачи 6
Обзор литературы 7
Глава 1. Структурные методы интегрирования 8
1.1 Классы структурных особенностей систем ОДУ 8
1.1.1 Класс A 8
1.1.2 Класс B 9
1.1.3 Класс C 10
1.1.4 Обзор моделей 11
1.2 Метод интегрирования 12
Глава 2. Схемы шестого порядка класса B 16
2.1 Условия порядка 16
2.2 Семейство схем класса B 18
2.3 Интегрирование с переменным шагом. Оценка локальной погрешности 18
2.4 Схема RKB6(4){7F} 20
Глава 3. Численное исследование моделей 22
3.1 Модель 1: точка либрации L1 22
3.2 Модель 2: орбита Аренсторфа 23
3.3 Результаты моделирования 24
Выводы 28
Заключение 29
Список литературы 30
Список рисунков 33
Список таблиц 34
Приложение А. Семейство методов шестого порядка класса B 35
Приложение Б. Текст программы ode46b 40

📖 Аннотация

Работа посвящена разработке экономичных схем численного интегрирования для структурно разделённых систем обыкновенных дифференциальных уравнений. Актуальность исследования обусловлена фундаментальными ограничениями классических методов Рунге—Кутты высокого порядка, известными как барьеры Бутчера, которые делают их вычислительно затратными, и потребностью в эффективных алгоритмах для моделирования сложных динамических систем, таких как задачи небесной механики. В рамках методологии в работе представлено параметрическое семейство шестиэтапных явных методов шестого порядка, относящихся к классу B структурно разделённых систем, и на его основе получена конкретная расчётная схема RKB6(4){7F}. Данная схема демонстрирует экономичность за счёт меньшего количества вычислений правой части системы по сравнению с существующими аналогами при сохранении высокого порядка точности. На базе разработанного метода создана функция ode46b с автоматическим выбором шага, которая при том же количестве этапов, что и широко используемый интегратор ode45, обеспечивает на один порядок точности выше. Практическая значимость результатов подтверждена численным исследованием на моделях механики, включая задачу о движении в окрестности точки либрации и орбиту Аренсторфа, где продемонстрирована эффективность предложенного подхода. Таким образом, разработанные схемы и программная реализация могут быть непосредственно использованы исследователями и инженерами для высокоточного и экономичного моделирования динамических систем в астрономии, космодинамике и смежных областях.

📖 Введение

Системы дифференциальных уравнений являются одним из основных способов математического моделирования. Математическое описание сложных динамических процессов, протекающих в окружающем мире, зачастую приво­дит к системам, не имеющим аналитического решения. Кроме того, в некоторых задачах (например, в линейных СОДУ большого размера) вывод точного решения может оказаться сопряжён с трудозатратными вычислениями, и по­лучение приближённого решения с некоторой желаемой точностью будет более оправданным. Поэтому численные методы решения систем дифференциальных уравнений всегда будут необходимым инструментом математического модели­рования.
В начале XX века немецкий математик Карл Рунге разработал теорию гра­фических численных методов, показав, как с помощью простейших инженерных инструментов на бумаге производить построения, аналогичные сложным мате­матическим операциям вплоть до приближённого интегрирования скалярных дифференциальных уравнений [1]. Развитие вычислительной техники привело к тому, что теория Рунге оказалась не нужна, однако предложенные им идеи привели к созданию явных одношаговых методов Рунге—Кутты, позволявших получать приближённое решение вплоть до четвёртого порядка точности при сравнительно небольших трудозатратах: для точности порядка р требовалось р вычислений правой части СОДУ (р-этапные методы одношаговые методы, р 6 4).
Технологический рывок середины XX века принёс новые возможности ЭВМ и вместе с тем предъявил новые требования к точности численного интегрирования. Безрезультатные годы поиска пятиэтапных методов пятого по­рядка завершились разработкой теории Джона Бутчера, систематизировавшего процесс нахождения методов Рунге—Кутты. Бутчер показал, что для высоких порядков существуют ограничения («барьеры Бутчера»): при р > 5 методы р-ого порядка требуют не менее р +1 вычисление правой части СОДУ [2], а при р > 7 — не менее р + 2 [3; 4]. Эти ограничения а также трудоёмкость по­лучения новых методов высокого порядка стали основными аргументами для поиска альтернатив классическим методам типа Рунге—Кутты.
Начиная с 1970-х лет наступает время бурного развития новых методов интегрирования и способов сравнения одних методов с другими. Важными свойствами становятся не только порядок точности, но также возможность автоматического управления длиной шага и устойчивость методов. Эрвин Фель­берг одним из первых предложил идею вложенных методов Рунге—Кутты, позволяющих существенно экономить вычисления с помощью автоматического управления длиной шага интегрирования [5]. Благодаря этому методы Рунге— Кутты обрели новую популярность, и вскоре на их основе были созданы схе­мы Дж. Дорманда и П. Принса [6], одна из которых в качестве функции ode45 сейчас является основным интегратором среды MATLAB.

Возникли сложности?

Нужна качественная помощь преподавателя?

👨‍🎓 Помощь в написании

✅ Заключение

1. В работе в явном виде представлено семипараметрическое семейство шестиэтапных методов шестого порядка класса B.
2. При фиксированных значениях параметров представлена расчётная схема численного интегрированияа RKB6(4){7F}. Найденная схемы экономичны в плане меньшего количества вычислений правой части СОДУ по сравнению с уже существующими методами.
3. На базе расчётной схемы создана функция, аналогичная встроенному интегратору MATLAB ode45 в части автоматического выбора шага. При том же количестве этапов представленная функция ode46b имеет на один порядок точности больше.
4. Проведено численное исследование различных моделей механики с помощью полученных методов численного интегрирования. Продемонстрирована эффективность в сравнении с известными классическими методами, в том числе с ode45.
Нужна своя уникальная работа?
Срочная разработка под ваши требования
Рассчитать стоимость
ИЛИ

📕 Список литературы

1. Runge, C. Graphical Methods: A Course of Lectures Delivered in Columbia University, New York, October, 1909, to January, 1910 / C. Runge. –– Columbia University Press, 1912. –– (... Graphical Methods).
2. Butcher, J. C. On Runge–Kutta Processes of High Order / J. C. Butcher. –– 1964. –– May.
3. Butcher, J. C. On the Attainable Order of Runge–Kutta Methods / J. C. Butcher // Math. Comput. –– 1965. –– Vol. 19. –– P. 408––417.
4. Butcher, J. C. The non-existence of ten Stage eight Order Explicit Runge–Kutta Methods / J. C. Butcher // BIT. –– 1985. –– Vol. 25. –– P. 521––540.
5. Fehlberg, E. Classical fifth-, sixth-, seventh-, and eighth-order Runge–Kutta formulas with stepsize control / E. Fehlberg // NASA Technical Report 287. –– 1968.
6. Dormand, J. R. A Family of Embedded Runge–Kutta Formulae / J. R. Dormand, P. J. Prince // Journal of Computational and Applied Mathematics. –– 1980. –– Mar. –– Vol. 6. –– P. 19––26.
7. Олемской, И. В. Численный метод интегрирования систем обыкновенных дифференциальных уравнений / И. В. Олемской // Математические методы анализа управляемых процессов. — 1986. — С. 157—160.
8. Hofer, E. A partially implicit method for large stiff systems of ODEs with only few equations introducing small time-constants / E. Hofer // SIAM J. Numer. Anal. –– 1976. –– Vol. 13, no. 5. –– P. 645––663.
9. Sandu, A. Multirate generalized additive Runge-Kutta methods / A. Sandu, M. Günther // Numer. Math. –– 2016. –– Vol. 133, no. 3. –– P. 497––524.
10. McLachlan, R. High order multisymplectic Runge-Kutta methods / R. McLachlan, B. Ryland, Y. Sun // SIAM J. Sci. Comput. –– 2014. –– Vol. 36, no. 5. –– A2199––A2226.
11. Wang, D. Parametric symplectic partitioned Runge-Kutta methods with energy-preserving properties for Hamiltonian systems / D. Wang, A. Xiao, X. Li // Comput. Phys. Comm. –– 2013. –– Vol. 184, no. 2. –– P. 303––310.
12. Ketcheson, D. Spatially partitioned embedded Runge-Kutta methods / D. Ketcheson, C. MacDonald, S. Ruuth // SIAM J. Numer. Anal. –– 2013. –– Vol. 51, no. 5. –– P. 2887––2910.
13. Eremin, A. S. Continuous Extensions for Structural Runge—Kutta Methods / A. S. Eremin, N. A. Kovrizhnykh // Computational Science and Its Applications – ICCSA 2017. –– Cham : Springer International Publishing, 2017. –– P. 363––378. –– (Lecture Notes in Computer Science ; 10405).
14. Шиманчук, Д. В. Построение траектории возвращения в окрестность коллинеарной точки либрации системы Солнце–Земля / Д. В. Шиманчук, А. С. Шмыров // Вестн. С.-Петерб. ун-та, Сер. 10. — 2013. — Т. 2. — С. 76—85.
15. Clohessy, W. H. Terminal Guidance System for Satellite Rendezvous / W. H. Clohessy, R. S. Wiltshire // J. Astronaut. Sci. –– 1960. –– Vol. 27, no. 9. –– P. 653––678.
...

🖼 Скриншоты

🛒 Оформить заказ

Работу высылаем в течении 5 минут после оплаты.
Предоставляемые услуги, в том числе данные, файлы и прочие материалы, подготовленные в результате оказания услуги, помогают разобраться в теме и собрать нужную информацию, но не заменяют готовое решение.
Укажите ник или номер. После оформления заказа откройте бота @workspayservice_bot для подтверждения. Это нужно для отправки вам уведомлений.

©2026 Cервис помощи студентам в выполнении работ