ВВЕДЕНИЕ 3
Постановка задачи 4
ГЛАВА 1. ВВЕДЕНИЕ В ПРЕДМЕТНУЮ ОБЛАСТЬ 5
1.1. Динамические системы 5
1.2. Численные методы решения дифференциальных уравнений 7
1.2.1. Метод рядов Тейлора 7
1.2.2. Методы Рунге-Кутты 9
1.3. Осциллятор Ван дер Поля 11
ГЛАВА 2. ПРОВЕДЕНИЕ ИССЛЕДОВАНИЯ 12
2.1 Реализация программного продукта 12
2.1.1. Требования к программному продукту 12
2.1.2. Выбор математического пакета 14
2.1.3. Архитектура приложения 15
2.2. Параметрическое исследование системы осциллятора Ван дер Поля 18
ЗАКЛЮЧЕНИЕ 22
СПИСОК ЛИТЕРАТУРЫ 23
ПРИЛОЖЕНИЕ A. ИСПОЛЬЗУЕМАЯ НОТАЦИЯ 25
ПРИЛОЖЕНИЕ B. ГРАФИЧЕСКИЙ ИНТЕРФЕЙС ПРОГРАММЫ INTEGRATION С ПОЯСНЕНИЯМИ 26
ПРИЛОЖЕНИЕ С. КОД ПРОГРАММЫ INTEGRATION 28
Современная наука стремится не только описать поведение изучаемого объекта, но и научиться предсказывать изменения, которые могут с ним произойти с течением времени. Конечная цель состоит в отыскании некоторого закона, позволяющего по ограниченным входным данным (например, о начальном состоянии объекта и внешним воздействиям) с точностью установить будущее объекта в любой момент времени.
Объекты и процессы, для которых однозначно определено понятие состояния - совокупности некоторых величин в конкретный момент времени, и закон, описывающий их изменение во времени, называются динамическими системами. Такие объекты и процессы встречаются в самых разных областях научного знания: в механике и других разделах физики, в биологии и химии, в сфере вычислительных процессов и обработки информации и многих, многих других [1, 2, 6].
Эволюция динамической системы может быть описана с использованием самых разных математических моделей. Одним из наиболее применимых средств для решения этой задачи являются системы дифференциальных уравнений. Поскольку единого метода получения точного решения систем дифференциальных уравнений не существует, чаще всего прибегают к помощи численных методов решения таких систем. Однако, для разных задач требуются различные методы. Известно, что численные методы, хотя и позволяют решать большинство задач данного вида, но всё же далеки от идеала. Каждый из таких методов обладает методической погрешностью, на которую влиять без изменения метода невозможно, но и это ещё не всё.
Важной проблемой при использовании численных методов решения становится выбор правильного метода. Дело в том, что некоторые совершенно не подходят для определённых задач (например, т.н. «жёсткие» задачи не могут быть решены простыми явными методами численного интегрирования). [10, 11]
Кроме того, современное исследование динамической системы может включать в себя работу с параметрами. Зачастую, исследователю необходимо не просто решить какую-то конкретную систему дифференциальных уравнений, а изучить целое влияние некоторых параметров на её решение, многократно получая решения для различных значений параметров.
Проведение параметрического исследования динамической системы потребовало изучения основные понятия области динамических систем, а так же разработки программа, позволяющая провести такое исследование для системы, представимой в виде системы двух дифференциальных уравнений с числом параметров до трёх. Проведено параметрическое исследование системы осциллятора Ван дер Поля и определено влияние параметров на эту систему. На основании полученных результатов сделан вывод о полезности проведения параметрического исследования динамической системы в реальных задачах. В дальнейшем работа может быть продолжена в различных направлениях.
1) Усложнение вводимых данных. Возможна разработка программного обеспечения с аналогичными функциями, но для более сложных динамических систем.
2) Упрощение работы с интерфейсом.
3) Реализация более сложных, но более продуктивных алгоритмов решений дифференциальных уравнений.
4) Автоматизация параметрического исследования.
Поведённое исследование позволяет вынести на защиту следующие положения:
1) изучены основные термины, касающиеся динамических систем и методы численного решения дифференциальных уравнений;
2) разработан программный продукт, позволяющий провести параметрическое исследование динамической системы;
3) проведено параметрическое исследование осциллятора Ван дер Поля, доказывающее работоспособность программы.
1) Андрианов С.Н. Динамическое моделирование систем управления пучками частиц. СПб.: Изд-во С.- Петерб. ун-та, 2002. 376 с.
2) Анищенко В.С. - Знакомство с нелинейной динамикой. АНО «Институт компьютерных технологий», 2002. 144 с.
3) Д. В. Аносов, Гладкие динамические системы. Гл.1. Исходные понятия, Итоги науки и техн. Сер.Соврем. пробл. мат. Фундам. направления, 1985,том 1, 156-178
4) Амелькин В.В. - Дифференциальные уравнения в приложениях. М.: Наука, 1986. 154 с.
5) Бахвалов, Н. С. - Численные методы / Н. С. Бахвалов, Н. П. Жидков, Г. М. Кобельков. - М.: Наука, 2003. 629 с.
6) Бутенин Н. В., Неймарк Ю. И., Фуфаев Н. Л. - Введение в теорию нелинейных колебаний. М.: Наука, 1987. 382 с.
7) Калиткин Н. Н. - Численные методы. М.: Наука, 1978. 500 с.
8) В.И. Мышенков, Е.В. Мышенков.- Численные методы, часть вторая. Численное решение обыкновенных дифференциальных уравнений. Издательство Московского государственного университеталеса, 2005. 108 с.
9) CohenD., JahnkeT., LorenzK., LubichC. NumericallntegratorsforHighly Oscillatory // Analysis, Modeling and Simulation of Multiscale Problems, 2006. P. 553 - 576
10) Ernst Hairer, Introduction to geometric numerical integration, /Springer Series in Computational Mathematics, Vol. 31, 2002
11) Rangarajan G. Symplectic integration of nonlinear Hamiltonian systems/ Pramana - journal of physics, 1997. Vol. 48, No. 1. P. 129 - 142/
12) http://bourabai.ru/cm/le 13.htm# 13 // Интернет-ресурс, содержащий примеры использования различных функций Maple
13) http://www.maplesoft.com/products/Maple/// Официальный портал Maple, содержащий всю необходимую информацию о пакете
14) http://www. maplesoft.com/support/help/// Онлайн-версия документацииMaple