Найти точки равновесия системы дифференциальных уравнений онлайн

Порядок производной указывается штрихами — y»’ или числом после одного штриха — y’5

Ввод распознает различные синонимы функций, как asin , arsin , arcsin

Знак умножения и скобки расставляются дополнительно — запись 2sinx сходна 2*sin(x)

Список математических функций и констант :

• ln(x) — натуральный логарифм

• sh(x) — гиперболический синус

• ch(x) — гиперболический косинус

• th(x) — гиперболический тангенс

• cth(x) — гиперболический котангенс

• sch(x) — гиперболический секанс

• csch(x) — гиперболический косеканс

• arsh(x) — обратный гиперболический синус

• arch(x) — обратный гиперболический косинус

• arth(x) — обратный гиперболический тангенс

• arcth(x) — обратный гиперболический котангенс

• arsch(x) — обратный гиперболический секанс

• arcsch(x) — обратный гиперболический косеканс

Видео:Математика это не ИсламСкачать

Digiratory

Видео:Решение системы дифференциальных уравнений методом ЭйлераСкачать

Лаборатория автоматизации и цифровой обработки сигналов

Видео:Дифференциальные уравнения 6. Фазовые траектории. Особые точки автономных системСкачать

Устойчивость нелинейных систем

Анализ устойчивости систем является одним из важнейших этапов проектирования систем управления, однако при анализе нелинейных, строго говоря, нет метода отвечающего критериям необходимости и достаточности, а критерии являются, как правило только достаточным (для устойчивости). Исходя из этого, для некоторых систем невозможно однозначно говорить о неустойчивости.

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

Рассмотрим применение классических методов Ляпунова.

Видео:Системы дифференциальных уравненийСкачать

Первый метод Ляпунова

Позволяет судить об устойчивости положения равновесия по линеаризованным уравнениям. Метод основан на утверждениях:

• если собственные значения линеаризованной системы имеют отрицательные действительные части (линеаризованная система асимптотически устойчива), то положение равновесия нелинейной системы устойчиво «в малом»;
• если среди собственных значений линеаризованной системы имеются «правые», то положение равновесия нелинейной системы неустойчиво;
• если имеются некратные собственные значения на мнимой оси, а остальные — «левые», то в этом критическом случае по линеаризованной модели нельзя судить об устойчивости положения равновесия нелинейной системы.

Таким образом для анализа системы по первому методу Ляпунова необходимо:

1. Найти положение равновесия системы — движений в системе нет (т.е. скорости и ускорения равны нулю) [ frac <mathrmv_><mathrmt>= ]
2. Линеаризовать систему в окрестности точки равновесия
3. Записать полученное линеаризованное дифференциальное уравнение в матричной форме (составить матрицу А)
4. Составить характеристический полином линеаризованной системы: [ = ]
5. Найти корни характеристического полинома. По виду корней сделать заключение о характере процессов в системе.

Основными недостатками первого метода Ляпунова являются:

• Если имеется корень на мнимой оси, то невозможно сказать о поведении процессов в системе.
• Возможно говорить только об устойчивости «в малом», т.е. при больших отклонениях от положения равновесия система может быть неустойчивой.

Пример 1.

Исследуем систему описываемую дифференциальными уравнениями:

Шаг 1. Положение равновесия:

Для нахождения точек равновесия левые части уравнений приравниваются к 0, что эквивалентно тому, что переменные состояния являются константами, а все их производные равны 0.

Шаг 2. Линеаризация для малых отклонений

Для линеаризации малых отклонений в точке равновесия старшие степени переменных, входящих в уравнения принимаются равными нулю.

Шаг 3. Линеаризованное управление в матричной форме

Преобразуем полученную линейную систему уравнений в матричный вид.

Шаг 4. Характеристический полином

Шаг 5. Корни характеристического полинома

Приравниваем характеристический полином к 0 и находим корни уравнения.

Заключение об устойчивости системы

в данном примере при линеаризации система имеет два корня с отрицательной вещественной частью, т.е. мы можем сказать, что система устойчива «в малом» (при больших отклонениях система может быть неустойчива).

Подтвердим теоретический вывод компьютерным моделированием (построением фазового портрета)

При этом, при начальных условиях, находящиеся дальше от точки равновесия, система становится неустойчивой

Пример 2. Нелинейный осциллятор

В качестве второго примера рассмотрим нелинейный осцилятор описываемый системой дифференциальных уравнений:

Аналогично первому примеру выполняем последовательность шагов

Шаг 1. Положение равновесия:

Шаг 2. Линеаризация для малых отклонений

Шаг 3. Линеаризованное управление в матричной форме

Шаг 4. Характеристический полином

Шаг 5. Корни характеристического полинома

Заключение об устойчивости системы

Рассматриваемая система является критическим случаем о ее устойчивости невозможно судить по линеаризованным уравнениям, применяемым в первом методе Ляпунова.

Видео:7. Линейные дифференциальные уравнения первого порядка. Метод Бернулли.Скачать

Второй метод Ляпунова

Второй метод Ляпунова не связан с линеаризацией системы, поэтому также называется прямым методом.

Для начала необходимо ввести понятия знакоопределенной, знакопостоянной и знакопеременной функций. Пусть имеется функция нескольких переменных:

Функция (V ) называется знакоопределенной в некоторой области, если она во всех точках этой области вокруг начала координат сохраняет один и тот же знак и нигде не обращается в нуль, кроме только самого начала координат

[ left ( Vleft ( bar right )=0 right ) ]

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

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

Теорема Ляпунова об устойчивости нелинейных систем

Если при заданных в форме

уравнениях системы n-го порядка можно подобрать такую знакоопределенную функцию Ляпунова

чтобы ее производная по времени

тоже была знакоопределенной (или знакопостоянной), но имела знак противоположный знаку (V), то данная система устойчива.

Для упрощения скажем, что функция Ляпунова должна быть положительной знакоопределенной функцией. Тогда условия теоремы Ляпунова будут выглядеть следующим образом:

Для устойчивости положения равновесия достаточно существования дифференцируемой функции

называемой функцией Ляпунова, удовлетворяющей в окрестности начала координат следующим условиям:

1. (Vleft ( v_, v_,…, v_right ) geq 0) причем (V=0) лишь при следующем условии, означающем что функция (V) имеет строгий минимум в начале координат. [ bar= begin v_ \ vdots \ v_ end = bar ]
2. Производная функции по времени [ frac <mathrmVleft ( bar right )><mathrmt>=sum_^frac<partial v_>frac <mathrmv_><mathrmt>=begin frac<partial v_> & frac<partial v_> & cdots & frac<partial v_>endbeginfrac <mathrmv_><mathrmt>\ frac <mathrmv_><mathrmt>\ vdots \ frac <mathrmv_><mathrmt>end ] в силу дифференциального уравнения (frac <mathrmbar><mathrmt>=barleft ( bar right ) ) является отрицательной знакопостоянной функцией, т.е. [ frac <mathrmVleft ( bar right )><mathrmt>=gradbarcdot frac <mathrmbar><mathrmt>=gradbarcdot barleft ( bar right )leq 0 ] при (tgeq t_)

Таким образом, условия:

1. (frac <mathrmVleft ( bar right )><mathrmt>leq 0) и функция (Vleft ( v_, v_,…, v_right ) ) является положительной знакоопределенной — это является достаточным условием устойчивости
2. (frac <mathrmVleft ( bar right )><mathrmt> ) — отрицательно определенная — это является достаточным условием асимптотической устойчивости.
3. (left | v right |rightarrow infty : frac <mathrmVleft ( bar right )><mathrmt>rightarrow infty ) — достаточное условие устойчивости «в целом».

Для анализа системы по второму методу Ляпунова необходимо:

1. Выбрать функцию Ляпунова от n переменных, где n- порядок системы.
2. Найти частные производные по переменным.
3. Вычислить производную функции по времени (frac <mathrmVleft ( bar right )><mathrmt>). Проанализировать полученный знак производной.

Из-за того, что второй метод Ляпунова не связан с линеаризацией, он считается универсальным. Однако он имеет ряд недостатков:

• Нет общих требований по выбору функции V
• Достаточный характер утверждения (если условия не выполняются, то об устойчивости ничего сказать нельзя, а можно посоветовать подобрать другую функцию (V ))

Пример 3. Нелинейный осциллятор

Проанализируем систему из примера (2).

Шаг 1. Функция Ляпунова

Для начала необходимо выбрать функцию Ляпунова от 2-х переменных (т.к. два вектора состояния):

Шаг 2. Частные производные

Шаг 3. Производная функции

Подставим в выражение значения исходя из ДУ:

Заключение об устойчивости системы

Исследовав систему первым методом Ляпунова мы не смогли сделать конкретный вывод об устойчивости системы, что позволил нам сделать второй метод Ляпунова. В результате мы можем сделать вывод, что система является асимптотически устойчивой.

Аналогично проверим с помощью моделирования:

Пример 4.

Рассмотрим систему, описываемую следующей системой дифференциальных уравнений:

Очевидно, что применение первого метода Ляпунова невозможно, т.к. матрица А состоит из нулей, а, следовательно, собственные значения равны нулю. Поэтому применим второй метод Ляпунова:

Шаг 1. Функция Ляпунова

Выбор функции Ляпунова второго порядка

Шаг 2. Частные производные

Шаг 3. Производная функции

При (a=3) имеет место асимптотическая устойчивость.

Заключение об устойчивости системы

Система является устойчивой.

Фазовый портрет системы выглядит следующим образом:

Видео:Особые точки 4 ЗадачаСкачать

Дифференциальные уравнения по-шагам

Видео:Устойчивость 1 ОпределениеСкачать

Результат

Примеры дифференциальных уравнений

• Простейшие дифференциальные ур-ния 1-порядка
• Дифференциальные ур-ния с разделяющимися переменными
• Линейные неоднородные дифференциальные ур-ния 1-го порядка
• Линейные однородные дифференциальные ур-ния 2-го порядка
• Уравнения в полных дифференциалах
• Решение дифференциального уравнения заменой
• Смена y(x) на x в уравнении
• Другие

Указанные выше примеры содержат также:

• квадратные корни sqrt(x),
  кубические корни cbrt(x)
• тригонометрические функции:
  синус sin(x), косинус cos(x), тангенс tan(x), котангенс ctan(x)
• показательные функции и экспоненты exp(x)
• обратные тригонометрические функции:
  арксинус asin(x), арккосинус acos(x), арктангенс atan(x), арккотангенс actan(x)
• натуральные логарифмы ln(x),
  десятичные логарифмы log(x)
• гиперболические функции:
  гиперболический синус sh(x), гиперболический косинус ch(x), гиперболический тангенс и котангенс tanh(x), ctanh(x)
• обратные гиперболические функции:
  asinh(x), acosh(x), atanh(x), actanh(x)
• число Пи pi
• комплексное число i

Правила ввода

Можно делать следующие операции

2*x — умножение 3/x — деление x^3 — возведение в степень x + 7 — сложение x — 6 — вычитание Действительные числа вводить в виде 7.5, не 7,5

Чтобы увидеть подробное решение,
помогите рассказать об этом сайте:

📽️ Видео

18+ Математика без Ху!ни. Дифференциальные уравнения.Скачать

Дифференциальные уравнения. 11 класс.Скачать

11. Уравнения в полных дифференциалахСкачать

Решение дифференциальных уравнений. Практическая часть. 11 класс.Скачать

Урок 5 - Дифференциальные уравнения, системы диффуровСкачать

Лабораторная работа 1. Решение систем обыкновенных дифференциальных уравненийСкачать

Задача Коши ➜ Частное решение линейного однородного дифференциального уравненияСкачать

Частное решение дифференциального уравнения. 11 класс.Скачать

Решение дифференциальных уравнений. Практическая часть. 11 класс.Скачать

Асташова И. В. - Дифференциальные уравнения. Часть 2 - Фазовый портретСкачать

Дифференциальные уравнения с разделяющими переменными. 11 класс.Скачать

Особые точки 1 Узел, седло, дикритический узелСкачать