Интегральные функции:
- Si(x)
- Интегральный синус от x
- Ci(x)
- Интегральный косинус от x
- Shi(x)
- Интегральный гиперболический синус от x
- Chi(x)
- Интегральный гиперболический косинус от x
В выражениях можно применять следующие операции:
- Действительные числа
- вводить в виде 7.
3- — возведение в степень
- x + 7
- — сложение
- x — 6
- — вычитание
- 15/7
- — дробь
3Другие функции:
- asec(x)
- Функция — арксеканс от x
- acsc(x)
- Функция — арккосеканс от x
- sec(x)
- Функция — секанс от x
- csc(x)
- Функция — косеканс от x
- floor(x)
- Функция — округление x в меньшую сторону (пример floor(4.5)==4.0)
- ceiling(x)
- Функция — округление x в большую сторону (пример ceiling(4.5)==5.0)
- sign(x)
- Функция — Знак x
- erf(x)
- Функция ошибок (или интеграл вероятности)
- laplace(x)
- Функция Лапласа
- asech(x)
- Функция — гиперболический арксеканс от x
- csch(x)
- Функция — гиперболический косеканс от x
- sech(x)
- Функция — гиперболический секанс от x
- acsch(x)
- Функция — гиперболический арккосеканс от x
Постоянные:
- Число «Пи», которое примерно равно ~3. 14159..
- e
- Число e — основание натурального логарифма, примерно равно ~2,7183..
- i
- Комплексная единица
- oo
- Символ бесконечности — знак для бесконечности
14159..6.4 Устойчивость систем автоматического регулирования. Частотный критерий устойчивости Михайлова / Хабр
1. Введение в теорию автоматического управления.2. Математическое описание систем автоматического управления 2.1 — 2.3, 2.3 — 2.8, 2.9 — 2.13.
3. ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ ЗВЕНЬЕВ И СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ РЕГУЛИРОВАНИЯ. 3.1. Амплитудно-фазовая частотная характеристика: годограф, АФЧХ, ЛАХ, ФЧХ. 3.2. Типовые звенья систем автоматического управления регулирования. Классификация типовых звеньев. Простейшие типовые звенья. 3.3. Апериодическое звено 1–го порядка инерционное звено. На примере входной камеры ядерного реактора. 3.4. Апериодическое звено 2-го порядка.
3.5. Колебательное звено. 3.6. Инерционно-дифференцирующее звено. 3.7. Форсирующее звено. 3.8. Инерционно-интегрирующее звено (интегрирующее звено с замедлением). 3.9. Изодромное звено (изодром). 3.10 Минимально-фазовые и не минимально-фазовые звенья. 3.11 Математическая модель кинетики нейтронов в «точечном» реакторе «нулевой» мощности.
4. Структурные преобразования систем автоматического регулирования.
5. Передаточные функции и уравнения динамики замкнутых систем автоматического регулирования (САР).
6. Устойчивость систем автоматического регулирования. 6.1 Понятие об устойчивости САР. Теорема Ляпунова. 6.2 Необходимые условия устойчивости линейных и линеаризованных САР. 6.3 Алгебраический критерий устойчивости Гурвица.
6.4. Частотный критерий устойчивости Михайлова
Советским ученым Михайловым в 30-тых годах впервые был предложен оригинальный критерий оценки устойчивости САР, основанный на исследовании частотных свойств полинома при подстановки вместо , где
Связь между частотными свойствами системы и передаточными функцииями более подробно описана в лекции 3.
Частотные характеристики САР
Подставим в формулу 6.4.1
Совершенно очевидно, что:
Критерий устойчивости Михайлова:
Чтобы САР (замкнутая или разомкнутая) была устойчивой, необходимо и достаточно, чтобы годограф при изменении от нуля до переходил поочередно из квадранта в квадрант против часовой стрелки, совершив при этом поворот на угол , где — степень полинома .
На рисунке 6.4.1 представлены варианты годографов для различных степеней n полинома
Рисунок 6.4.1 Годографы устойчивых систем
Если САР устойчива, то вектор при изменении от 0 до совершает поворот на угол , где n — степень полинома .
Следствием частотного критерия Михайлова является перемежаемость (чередование) нулей полиномов и в самом деле (см. рисунок 6.4.1), для кривой с последовательность пресечения осей получается как на рисунке 6.4.2
Рисунок 6.4.2. Нули полиномови чередуются для устойчивости системы.Если система находится на апериодической границе устойчивости (один нулевой полюс при всех остальных в левой полуплоскости), то годограф имеет следующий «примерный» вид, см.
рисунок 6.4.3:
Рисунок 6.4.3 — Годограф начинается из начала координат и поочередно «проходит» все квадранты в положительном направлении (начиная со 2-го квадранта).
Если система находится на колебательной границе устойчивости (2 чисто мнимых полюса при всех остальных в левой полуплоскости), то годограф имеет вид как на рис. 6.4.4:
Рисунок 6.4.4 — Годограф при некоторой частоте проходит через начало координат, «перескакивая» из 2-го в 4-ый квадрант (минуя 3-ий). Частота — частота незатухающих колебаний в такой САР.Если САР неустойчива, годографы имеют вид как представлено на рисунке 6.4.5:
Рисунок 6.4.5. Годографы неустойчивых систем.Докажем ряд основных моментов в критерии Михайлова.
Представим полином в виде произведения:
где — полюса главной передаточной функции.
Учитывая, что любое комплексное число типа можно представить в виде: где – модуль, — фаза.
Рассмотрим изменнение фазы, при изменении от 0 до .
Обозначим изменение фазы как
Для устойчивой САР, все полюса D(s) лежат в левой полуплоскости. (см. предыдущию лекцию) Рассмотрим различные варианты расположения полюсов на плоскости:
1-й случай: пусть является реальным числом , например, , где , поскольку корень в левой полуплоскости.
Рассмотрим поведение вектора при изменении от нуля до бесконечности
Рисунок 6.4.6. Вектор для реального корняИз рисунка 6.4.6 очевидно, что
при ,
при .
Т.е. при изменении от 0 до вектор, описывающий скобку повернется в положительном направлении на угол .
2-й случай: Пусть где, преобразуем скобку:
Рассмотрим изменение от до :
при — точка лежит в правом нижнем квадранте коплексной плоскости, фаза (сдвиг фазы):
(см. рис. 6.4.7)
при фаза (сдвиг фазы) (см. рис. 6.4.7)
Рисунок 6.4.7 Вектор длядля комплексного корня.Изменение фазы:
3-й случай: Пусть (полюс комплексно сопряженный со вторым вариантом).
Преобразуем скобку
Рассмотрим изменение от до :
при — точка лежит в правом верхнем квадранте коплексной плоскости, фаза (сдвиг фазы):
(см. рис. 6.4.8)
при фаза (сдвиг фазы) (см. рис. 6.4.8)
Рисунок 6.4.8 Вектор длядля комплексного корня.Изменение фазы:
Рассмотрим изменения годографа полинома устойчивой системы c учетом изменения фазы для трех случаев полюсов рассмотрены выше.
Пусть у нас общее количество полюсов , — количество сопряженных полюсов полинома, тогда количество вещественных полюсов.
Покольку вещественны полюс дает (cм. формулу 6.4.6), а два комплексно сопряженных корня в сумеее дают (формулы 6.4.7 и 6.4.8)
Это означает, что при изменении частоты от нуля до бесконечности, годограф должен поочередно пройти все квадранты в положительном направлении, если САР – устойчива.
Рассмотрим неустойчивую САР, у которой ряд полюсов полинома расположен в правой полуплоскости.
4-й случай: Пусть где — реальное число.
Преобразуем скобку подставля значения полюса .
Рассмотрим изменние вектора при изменении от 0 до. Примерный вид представлен на рисунке 6.4.9, где :
Рисунок 6.4.9 Вектор длядля комплексного корня.Изменение фазы вектора:
при ;
при .
Изменение фазы:
Следовательно отрицательный реальный полюс дает вращение вектора в отрицательном направлении на угол . Получается наличе одного реального полюса вызывает «недоповорот» вектора на угол .
Рассмотрим два варианта с коплексными полюсами лежашими в левой полуплоскости:
5-й случай: Пусть , где и преобразуем скобку для данного случая:
.
Рассмотрим изменние вектора при изменении от 0 до . Примерный вид представлен на рисунке 6.4.10, где: ,
Рисунок 6.4.10 Вектор длядля комплексного корня.Рисунок 6.4.10 Вектор длядля комплексного корня. Изменение фазы вектора:
при ;
при .
Изменение фазы:
6-й случай: Пусть , явялется комплексно сопряженным полюсом для 5-го случая, где и преобразуем скобку для данного случая:.
Рассмотрим изменние вектора при изменении от 0 до . Примерный вид представлен на рисунке 6.4.11, где: ,
Рисунок 6.4.11 Вектор длядля комплексного корня.при ;
при .
Изменение фазы:
При наличии двух комплексно-сопряженных корней в левой полуплоскости (варианты 5 и 6) общее изменение фазы вычисляется по формуле:
Резюмируем:
Если САР — устойчива все полюса полинома степенью лежат в левой полуплоскости, то изменение фазы годографа при изменении on 0 доописывается формулой 6.4.13:
Если один полюс полинома степенью лежит в правой полуплоскости, а остальные в левой полуплоскости, то изменение фазы годографа при изменении on 0 доописывается формулой 6.4.14:
Если в правой полуплоскости расположено L полюсов полинома степенью , а остальные в левой полуплоскости, то изменение фазы годографа при изменении on 0 доописывается формулой 6.4.15:
Предельный случай
Если один из полюсов полинома явялется бесконечным (см.
рисунок 6.4.12):
Рисунок. 6.4.12 Бесконечный полюс
Данный случай возникает, если годограф в этом случае ведет себя как показано на рисунке 6.4.13:
Рисунок 6.4.13 Вид годографа c «бесконечным» корнямиПример
Исследовать на устойчивость САР , представленную на рисунке 6.4.14 с использованием критерия Михайлова
Рисунок 6.4.14 САР для исследования
Полином
Корни полинома : ;
Корни полинома ;
Чередования 0 для полиномов не происходит, (см. рис. 6.4.15)
Рисунок 6.4.15 Корни полиномов
Определим какие должны быть коэффициенты полинома и что бы САР была устойчивой согласно критерию Михайлов.
Для устойчивой системы, необходимо чередование корней, для нашего случа корни могут распологаться по возрастанию в следующем порядке:
Изменим коэффициент так, что-бы неравенство сталов верным. Например пусть . Тогда решая уравнение для , при получаем коэффициент
Проверим результат численным моделированием. Создадим стуркутурную схему, как показанао на рисунке 6.
4.15.
Рисунок 6.4.15 Схема модели для проверки решения примера 1.
Используем блок передаточная переменная общего вида, где будем задавать, коэффициет k, из условия задачи в качестве глобальной переменной, меня которую можно изменять коэффициент числителя (см. рис. 6.4.15). Зададим в качестве тестового воздействия ступеньку на 5-й секунду, так же поместим на схему блок построения гадогрофа Михайлова. Результаты расчет приведены на рисунке 6.4.16
Рисунок 6.4.16 Расчет неустойчивой системы.
Меняя коэффициет , можно убедиться, что система остается неустойчивой, при любых значениях коэффициента меняется только амплитуда колебательного процесса. График годографа Михайлова показывает, что годограф переходит из квадранта в квандрат комплексной полсокости,
Изменим коэффициент блока предаточной функции общего вида рассчитанные для получения чердования корней , и повторим расчет. Изменненная модель и результаты моделирования представлены нар рисунке 6.