Теперь остается решить представленные дифференциальные уравнения получить графики полученных решений, представленные на рис. 11.10. В данном случае частоты сигнала и собственных колебании системы заметно различаются и выходной сигнал системы представляет собой, в основном, выделенную вторую гармонику воздействия.

Рис. 11.10. Решение дифференциальных уравнений и его визуализация

Разумеется, представленный вариант анализа носит частный характер, поскольку синтезируется вполне конкретный вид сигнала — прямоугольные импульсы с заданными выше параметрами. Однако, если использовать разложение в ряд Фурье произвольного воздействия, то подобным способом можно решить задачу получения реакции колебательной (а, в принципе, любой линейной) системы на заданное воздействие.

Вернемся к задаче моделирования системы второго порядка и попытаемся найти решения в более удобном виде, обычно приводимом в учебниках после ряда преобразований. Для этого достаточно воспользоваться пакетом расширения DEtools. Рис. 11.11 показывает начало документа с составленным дифференциальным уравнением и его решением. Нетрудно заметить, что теперь решение представлено в классическом виде, который обычно приводится в учебниках по теории колебаний.

Рис. 11.11. Решение дифференциального уравнения свободных колебаний с применением пакета DEtools

На рис. 11.12 показана вторая часть документа с решением для конкретных данных и построением графика временной зависимости свободных колебаний. Нетрудно заметить, что свободные колебания системы имеют вид затухающих синусоидальных колебаний. Вы можете проверить, что при р<0 колебания будут нарастать по экспоненциальному закону, что характерно для генераторных систем.

Рис. 11.12. Пример вычисления временной зависимости свободных колебаний и построения их графика

Нередко о характере колебаний удобно судить по фазовому портрету колебаний. Он задается графиком в параметрической форме, при которой по одной оси откладывается зависимость у(t), а по другой — ее производная. Это показано на рис. 11.13. Фазовый портрет в данном случае представляет собой сворачивающуюся спираль.

Рис. 11.13. Фазовый портрет затухающих свободных колебаний

По аналогии с последним примером можно рассмотреть поведение системы второго порядка при синусоидальном воздействии. На рис. 11.14 представлено начало документа, в котором задано исходное дифференциальное уравнение и получено его общее и частное аналитические решения.

Рис. 11.14. Пример аналитического решения задачи на поведение системы второго порядка при синусоидальном воздействии

На рис. 11.15 представлены временные диаграммы реакции системы и синусоидального воздействия. Кроме того, построен фазовый портрет колебаний. Он заметно отличается от спирали и хорошо иллюстрирует сложность колебаний в начале их развития.

Рис. 11.15. Результаты моделирования цепи второго порядка при синусоидальном воздействии

К сожалению, применение пакета расширения DEtools усложняет функцию dsolve решения дифференциальных уравнений. В результате время моделирования даже простых систем удлиняется до минут, а более сложные системы могут потребовать куда более длительного времени моделирования. В этом случае может оказаться целесообразным отказаться от получения аналитических зависимостей для результатов моделирования и перейти к численному моделированию.

Рассмотрим методику улучшенного моделирования еще на одном примере — вычислении реакции системы при пилообразном воздействии. На рис. 11.16 показано задание такого воздействия с помощью функции floor. Для упрощения расчетных выражений амплитуда и период воздействия взяты равными я. Поскольку в данном случае аналитическое решение получить невозможно (функция floor не позволяет этого), то заменим воздействие рядом Фурье. Его коэффициенты также представлены на рис. 11.16.

Рис. 11.16. Начало моделирование системы с пилообразным воздействием, представленным рядом Фурье

На рис. 11.17 представлены графики воздействия в идеальном случае и при его представлении рядом Фурье с пятью гармоники. Показано также аналитическое решение для временной зависимости y(t) при таком воздействии.

Рис. 11.17. Воздействие и временная зависимость реакции системы при пилообразной форме воздействия