Содержание
Введение
Трудно представить современную жизнь без математики. Она присутствует в различных областях человеческой деятельности. Многие научные открытия, стали возможны благодаря математическим расчетам и решениям уравнений.
Дифференциальные уравнения - один из разделов высшей математики, с помощью которого осуществляется решение практических задач. Дифференциальные уравнения широко используются в изучении физических, технических, биологических, экономических процессов, природных явлений. С помощью их, стало возможным перевести реальную задачу на математический язык, сохраняя при этом основные свойства описываемого процесса. Решив полученные уравнения, осуществляется описание функционального закона, рассматриваемой предметной области. Все это способствует повышению достоверности получаемых знаний об окружающем нас мире.
Целью работы является исследование дифференциальных уравнений третьего порядка, с помощью компьютерных моделей, с заданными нами условиями. А также наблюдение за ними в процессе изменения исходных данных, с дальнейшим анализом получившихся результатов.
Выпускная квалификационная работа состоит из трех разделов.
Первый раздел содержит анализ предметной области, а в частности в нем описываются виды дифференциальных уравнений, а также способы их решения. Раскрываются основные понятия устойчивости, предельного цикла, а также рассматриваются частные случаи, описываемые в системе Лоренца.
Во втором разделе содержится общее описание поставленной задачи, а также определяются требования к рассматриваемым моделям и дифференциальным уравнениям третьего порядка.
В третьем разделе исследуется кусочно-линейное уравнение третьего порядка, которое является частным случаем решения дифференциальных уравнений. Исследования проходят при помощи различных графиков, моделей, так как этот способ дает более точное решение на поставленные нами условия уравнений.
1. Анализ предметной области
1.1 Виды дифференциальных уравнений
Уравнения, в которых сами функции и их производные являются неизвестными называются дифференциальными уравнениями. Дифференциальные уравнения играют важную роль при решении многих физических, технических, биологических задач. Если в уравнении находится неизвестная переменная, а также неизвестная функция и ее первая производная, то это уравнение является дифференциальным уравнением первого порядка. То есть порядок дифференциального уравнения, будет зависеть от наивысшего порядка производных, входящих в него.
Дифференциальные уравнения имеют следующие виды:
1) обыкновенные - к ним относятся только функции и их производные от одного аргумента:
2) стохастические дифференциальные уравнения - это уравнения, которые включают в себя случайные процессы:
3) уравнения с частными производными - это уравнения, в которых функции зависят от многих переменных.
Также дифференциальные уравнения бывают первого, второго и высших порядков.
1.2 Основные функции для решения дифференциальных уравнений в MathCad
компьютерный моделирование дифференциальный уравнение
Для решения дифференциальных уравнений использовался математический пакет MathCad. Благодаря встроенному внутреннему языку, решение дифференциальных уравнений выводилось не только в числовой форме, но и визуально. Для этого были построены специальные графики.
В математическом пакете MathCad существует множество функций, с помощью которых можно решать дифференциальные уравнения. При их решении, нами использовались следующие функции:
1) given .. odesolve;
2) rkfixed.
Пример решения дифференциального уравнения с помощью функции given .. odesolve представлен на рисунке 1.1.
Рисунок 1.1 - Пример решения с помощью функции given .. odesolve
Пример решения дифференциального уравнения с помощью функции rkfixed представлен на рисунке 1.2.
Рисунок 1.2 - Пример решения с помощью функции rkfixed
1.3 Понятие устойчивости и предельного цикла в MathCad
Если поведение исходного решения незначительно отличается от поведения решений, с условиями близкими к начальным, то такое решение дифференциального уравнения называется устойчивым.
Для некоторых дифференциальных уравнений просто положительности коэффициентов недостаточно. Поэтому число устойчивых корней зависит от параметров. Число устойчивых корней влияет на размерность устойчивости W многообразия, на котором расположены приближающиеся к нему траектории. Если многообразия будут двумерны, то на них можно увидеть устойчивые фокусы или узлы. Увидим ли мы узлы, либо фокусы, будет зависеть от знака дискриминанта. Если дискриминант меньше нуля, то будут узлы. Если же он больше, то будут фокусы. В свою очередь узлы и фокусы бывают устойчивыми и неустойчивыми.
Пример устойчивых и неустойчивых узлов соответственно представлен на рисунке 1.3.
Рисунок 1.3 - Устойчивые и неустойчивые узлы
Пример устойчивых и неустойчивых фокусов соответственно представлен на рисунке 1.4.
Рисунок 1.4 - Устойчивые и неустойчивые фокусы
Узел и фокус можно назвать устойчивыми, если в корнях уравнений все вещественные части отрицательны. В тоже время если вещественные части положительны, то узлы и фокусы неустойчивы.
Пример уравнения в Mathcad, решением которого является узел, представлен на рисунке 1.5.
Рисунок 1.5 - Уравнение решением которого является узел
Пример уравнения в Mathcad, решением которого является фокус, представлен на рисунке 1.6.
Рисунок 1.6 - Уравнение решением которого является фокус
Рисунок 1.7 - Система Лоренца
Особое внимание хотелось бы уделить предельным циклам, так как они широко применяются в современной жизни и напрямую связаны с устойчивостью. Так, если колебание описывается устойчивым циклом, то его называют автоколебанием. Автоколебание в свою очередь встречается в часах, радиопередатчиках, даже биение сердца, можно назвать автоколебанием. В описание процесса автоколебаний, свой вклад внес Эдвард Лоренц доказавший, что простейшая система, состоящая из дифференциальных уравнений, может привести к тому, что траектории примут хаотический характер. В результате своих работ, он исследовал множества, к которым приближаются точки в процессе итерации. Для данных множеств было введено понятие аттрактор, что в переводе с английского обозначает притягивать. Пример простейшего аттрактора, описываемого системой Лоренца, представлен на рисунке 1.7.
Для описания работы приведенных устройств, в которых встречается автоколебание и используются предельные циклы в соответствующем фазовом пространстве. Пример предельного цикла представлен на рисунке 1.8.
Исходя из приведенного выше примера, можно наглядно наблюдать, что представляет из себя предельный цикл. В данном случае он скручивается, постепенно приближаясь к границе, за которую он не сможет перейти. Точные границы найти невозможно, но их можно представить визуально. Это связано с тем что приближаясь, расстояние между моделью и предельно допустимой границей будет настолько незначительным, что в теории можно будет говорить о нахождении предельного цикла.
Для более наглядного примера, найдем полное решение для нахождения предельного цикла в MathCad, которое представлено на рисунке 1.9.
Рисунок 1.8 - Пример предельного цикла
Рисунок 1.9 - Пример программ в MathCad
Для еще более лучшего, визуального представления, построим данный предельный цикл в пространстве. Пример пространственного отображения предельного цикла представлен на рисунке 1.10.
Рисунок 1.10 - Пример предельного цикла в пространстве
Не только Лоренц изучал предельные циклы, также существенный вклад внес физик Бальтазар Ван дер Пол, который обнаружил, что на определенных частотах появляются шумы, которые находятся рядом с собстенными частотами волн.
Для еще более лучшего, визуального представления, построим данный предельный цикл в пространстве
Примером описания данного явления, может служить уравнение Ван дер Поля, которое представлено на рисунке 1.11
Рисунок 1.11 - Пример уравнения Ван дер Поля в MathCad
В процессе решения данного уравнения, получились графики, которые представлены на рисунке 1.12.
Рисунок 1.12 - Графики решения уравнения Ван дер Поля в MathCad
2. Постановка задачи
Создание компьютерных моделей для исследования дифференциальных уравнений третьего порядка, является главной задачей дипломной работы. При данном исследовании, были разобраны задачи на базовые понятия о дифференциальных уравнениях, а также рассмотрены случаи, где нельзя однозначно сказать о точном решении, то есть необходимо проводить анализ полученных результатов. В данном случае и приходят на помощь компьютерные модели, которые строились конкретно для каждой задаче.
Благодаря проделанным исследованиям, другие студенты с легкостью смогут понять основную информацию о дифференциальных уравнениях третьего порядка, потому что вся информация об этих исследованиях представлены в простом и понятном виде. Для того, чтобы продолжить исследования, не нужно по новой все изучать, примеры кодов в математическом пакете MathCad, также представлены в рамках данной дипломной работы.
Данные исследования проводились не только с целью дальнейшего обучения студентов, но они носили и немного другой характер. Как мы знаем, с помощью дифференциальных уравнений можно описывать различные процессы, которые происходят в природе, в космосе, а также в различных механических приборах. Поэтому показать широкое их применение и необходимость использования при описании какого-либо процесса, также являлось неотъемлемой частью данной дипломной работы.
3. Кусочно-линейное дифференциальное уравнение третьего порядка
Предельные циклы для кусочно-линейных уравнений третьего порядка плохо изучены. Особый интерес представляет обнаружение неустойчивости или устойчивости таких уравнений. При решении данных уравнений, основной акцент направлен на обнаружение таких коэффициентов, которые будут влиять на поведение уравнения. Коэффициенты не находятся явным образом, при их подборе нужно рассматривать как ведет себя уравнение в пространстве.
Пример программы в Mathcad, для нахождения предельного цикла кусочно-линейного уравнения третьего порядка представлен на рисунке 3.1.
Опираясь на выше указанную программу, было задано условие, которое представлено в виде формулы:
(3.1)
В результате подбора значений для данного условия, получились графики, которые представлены на рисунке 3.2
Рисунок 3.1 - Пример кода в Mathcad
Рисунок 3.2 - Графики для заданного условия
На графиках показано одно решение уравнения, но при этом взяты проекции на разные плоскости. То есть мы видим пространственное отображение решения по плоскостям xy, xz, yz соответственно. Посмотрев на графики, мы можем заметить, что они закручиваются. В условии у нас задано, что за время t, график должен пройти количество шагов n. Пройдя это расстояние, мы точно не можем сказать, стабилизировался ли он. Поэтому над этими графиками проводим итерацию. Пример итерации представлен на рисунке 3.3.
Рисунок 3.3 - Пример итерации
Из этих графиков мы видим, что в заданных нами условиях, график еще не зацикливается. Потому что после итерации он продолжает скручиваться. Стабилизация происходит на уровне около единицы. Проводя повторные итерации, можно получить более точные значения того, где будет происходить зацикливание данных моделей.
С каждым своим новым витком, графики скручиваются, следовательно их амплитуда уменьшается. Зависимость амплитуды от времени представлена на рисунке 3.4.
Рисунок 3.4 - Зависимость амплитуды от времени
Для более детального понимания того, как ведут себя графики в пространстве, построим трехмерные графики, которые представлены на рисунке 3.5.
Рисунок 3.5 - Трехмерные графики
На графиках представлены виды с разных осей, на первом ось xy, на втором соответственно xz - вид сбоку, а на третьем ось yz.
Исходя из выше приведенного примера можно сделать вывод, что модели на графиках скручивались, но они могут как скручиваться, так и раскручиваться. Пример кода в Mathcad, для раскручивающейся модели, представлен на рисунке 3.6.
Рисунок 3.6 - Пример кода в Mathcad
Исходя из выше приведенного примера можно сделать вывод, что модели на графиках раскручиваются. Пример раскручивающейся модели представлен на рисунке 3.7.
В приведенных ниже графиках, движение происходит не по круговой траектории, а описывает совершенно новую модель, которая существенно отличается от первой. Все это свидетельствует о том, что с помощью дифференциальных уравнений, можно описать бесконечно много моделей, будь то изменение уравнения или условия в данном уравнении.