6
ЛЕКЦИЯ 9
КАЧЕСТВО ПЕРЕХОДНЫХ
ПРОЦЕССОВ В ЛИНЕЙНЫХ СИСТЕМАХ
АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ
Понятие о качестве
переходных процессов.
Из предыдущей
лекции мы узнали, что автоматическая
система прежде всего должна быть
устойчива. В устойчивой системе переходный
процесс затухает. Однако для практики
вовсе не безразличен характер затухания
переходного процесса. Так, например,
если переходный процесс затухает
медленно, и система долго входит в новый
установившийся режим, то она обладает
недостаточным быстродействием, и,
следовательно, ее применение будет
ограничено. Поэтому устойчивость
является необходимым, но недостаточным
условием работоспособности автоматических
систем. Достаточным условием является
качество
процессов регулирования,
которое оценивается качеством
переходных процессов и ошибками в
установившихся режимах.
Качество переходных
процессов обычно оценивают по переходной
функции h(t),
которая представляет собой реакцию
системы на внешнее воздействие типа
единичной ступенчатой функции 1(t).
Для следящих систем и систем управления
рассматривают переходную функцию по
отношению к задающему воздействию h(t),
а для систем стабилизации – по отношению
к возмущению hf
(t).
На примере переходной
функции h(t)
(рис. 9.1) познакомимся с основными
показателями качества переходного
процесса: временем регулирования,
перерегулированием, основной частотой
колебаний, количеством колебаний,
максимальной скоростью и максимальным
ускорением регулируемой величины.
Время регулирования
tр
определяется длительностью переходного
процесса. теоретически переходный
процесс длится бесконечно долго, однако
практически считают, что он заканчивается,
как только отклонения регулируемой
величины от нового ее установившегося
значения не будут превышать допустимых
пределов ε.
Обычно принимают ε
= (3 — 5)% от hуст.
Временем регулирования характеризуют
быстродействие системы. однако иногда
быстродействие характеризуют также
временем tу
достижения переходной функцией первый
раз нового установившегося значения
или временем tmax
достижения максимального значения
hmax.
Перерегулирование
Δ hmax,
или выброс, представляет собой максимальное
отклонение регулируемой величины от
нового установившегося значения. Обычно
первый максимум является наибольшим.
Относительное перерегулирование
вычисляют по формуле
.
Время регулирования
и перерегулирование, основные показатели
переходного процесса, тесно связаны
между собой. На первый взгляд может
показаться, что перерегулирование
недопустимо, так как оно увеличивает
время регулирования. Однако это не так.
Действительно, перерегулирование
появляется вследствие того, что система
к новому установившемуся состоянию
подходит с определенной скоростью,
которая графически отображается
тангенсом угла наклона касательной в
точке А (рис. 9.1):
.
Чем больше эта
скорость, тем дальше за новое установившееся
положение пройдет система по инерции.
Значит, для уменьшения перерегулирования
необходимо уменьшить скорость, с которой
система подходит к новому установившемуся
состоянию. Это, в конце концов, приведет
к увеличению времени регулирования.
Если система подходит к установившемуся
состоянию с нулевой скоростью, то
перерегулирования вообще не будет, но
время регулирования значительно
увеличится (рис. 9.2,б).
Так как время
регулирования характеризует быстродействие
системы, то стремятся, там где это
необходимо, свести его к минимуму. Это
значит, что система должна отрабатывать
задающее воздействие с большой скоростью,
что в итоге приводит к большому
перерегулированию. Кроме того, поскольку
переходный процесс колебательный, то
будут иметь место большие ускорения
регулируемой величины, что может вызвать
недопустимо большие динамические
нагрузки на элементы системы. Таким
образом, как отсутствие, так и слишком
большое перерегулирование нежелательно.
Поэтому в качестве оптимального допускают
перерегулирование в пределах 20 – 30%.
При этом число полуколебаний переходной
функции равно 2 – 3.
Числовые значения
времени регулирования и перерегулирования
часто задают в качестве исходных данных
для синтеза корректирующих устройств,
поскольку правильным выбором и настройкой
последних обеспечивается подавление
нежелательных колебаний регулируемой
величины в переходном процессе.
Для некоторых
систем перерегулирование вообще
недопустимо, например, для систем
автоматического регулирования химических
процессов, где перерегулирование может
испортить выпускаемый продукт. Необходимо
также иметь в виду, что стремление
уменьшить время регулирования приводит
к увеличению мощности исполнительного
устройства.
Коррекция
динамических свойств линейных САУ.
Коррекция
динамических свойств САУ осуществляется
для выполнения требований по точности,
устойчивости и качеству переходных
процессов.
С точки зрения
требований к точности САУ в установившихся
режимах коррекция динамических свойств
САУ может понадобиться для увеличения
коэффициента передачи или порядка
астатизма при сохранении устойчивости
и определенного качества переходных
процессов.
Коррекция применяется
также как средство обеспечения
устойчивости неустойчивой системы или
расширения области устойчивости и
повышения качества переходных процессов.
Осуществляется
коррекция с помощью введения в систему
специальных корректирующих
звеньев с
особо подобранной передаточной функцией.
Различают 4 основных
вида корректирующих устройств.
1. Последовательные
корректирующие устройства
или, как их еще называют, корректирующие
фильтры,
которые могут описываться различными
передаточными функциями Wк(s)
(рис. 9.3).
Здесь
—
передаточная функция основных звеньев
системы;
—
передаточная функция корректирующих
звеньев.
В первом случае
(рис. 9.3,а) общая передаточная функция
разомкнутой цепи системы будет
.
Во втором случае
(рис. 9.3, б) имеем
,
где нулями отмечены
передаточные функции заданных частей
системы.
2. Параллельные
корректирующие устройства.
Выполняются в виде дополнительных
местных обратных связей
(рис. 9.4), когда
.
3. Корректирующие
устройства по внешнему воздействию.
4. Неединичная
главная обратная связь.
Передаточная
функция корректирующего устройства
или
может иметь, вообще говоря, произвольный
вид. Но наиболее часто применяются
определенные типы корректирующих
устройств, которые мы и рассмотрим.
Типовые
последовательные корректирующие
устройства.
Введение
производной от ошибки
– простейший метод улучшения качества
переходного процесса. Структурно
введение производной представлено на
рис. 9.5.
Технически это
может быть осуществлено различными
устройствами, причем производная может
осуществляться не в чистом виде, а с
инерционностью, например, в виде
.
Передаточная
функция разомкнутой цепи (рис. 9.5) с
идеальной производной будет
.
Заменив s=jω,
получим амплитуду и фазу в виде:
.
Существенным
здесь является то, что при введении
воздействия по производной добавляется
положительная фаза.
Вследствие этого
радиус-векторы амплитудно-фазовой
характеристики (рис. 9.6, а) поворачиваются
против часовой стрелки, увеличивая
запас устойчивости и улучшая качество
переходного процесса. То же самое легко
проследить и на логарифмических
характеристиках (рис. 9.6,б).
В случае неидеального
дифференцирования (с инерционностью)
этот эффект несколько уменьшается
количественно, но качественно он
сохраняется.
Заметим, что
введение производной от ошибки может
служить и стабилизирующим средством,
т.е. превращать неустойчивую замкнутую
систему в устойчивую, например, если на
рис. 9.6,а
точка (-1,j0)
лежала бы внутри характеристики
,
то новая характеристика
уже могла бы не охватывать точку (-1,j0).
На этой странице я собрала теорию и практику, готовые задачи и подробные решения по предмету теория автоматического управления, чтобы вы смогли освежить знания.
Если что-то непонятно — вы всегда можете написать мне в WhatsApp и я вам помогу!
ТАУ
Теория автоматического управления (ТАУ) является одной из немногих общепрофессиональных технических дисциплин, входящих под тем или иным названием во все программы инженерного образования. Основой ТАУ являются различные по идеям и методам исследования разделы высшей математики и физики, такие как дифференциальное и интегральное исчисление, теория функций комплексного переменного, теория матриц, теория оптимальных процессов, математическая логика, теория вероятности и случайные процессы, механика, электричество и магнетизм и др.
Математическое описание систем управления. Уравнения и передаточные функции
Система или звено с одним выходом и двумя входами
и
в общем случае описывается уравнением
или
где обозначает оператор дифференцирования
,
Дифференциальный оператор при выходной переменной называется собственным оператором, а дифференциальные операторы
и
при входных переменных
и
— операторами воздействия. Отношение оператора воздействия к собственному оператору называется передаточной функцией в операторной форме.
Степень полинома знаменателя передаточной функции называют порядком, а разность между ее степенями знаменателя и числителя — относительным порядком или относительной степенью передаточной функции и соответствующей ей системы.
Нулями и полюсами передаточной функции
называют нули ее числителя и знаменателя соответственно, т. е. корни уравнений и
, где
рассматривается как переменная, а не как оператор.
Система (1.1) определяется двумя передаточными функциями: передаточной функцией
относительно входа и передаточной функцией
относительно входа . Порядок этих передаточных функций равен
, а относительный порядок —
для передаточной функции
и
для передаточной функции
.
С помощью передаточной функции уравнение рассматриваемой системы управления можно записать в виде
Имеющее наименьший порядок отношение изображений Лапласа выходной и входной переменных, вычисленных при нулевых начальных условиях, называется передаточной функцией в изображениях Лапласа. В соответствии с определением передаточная функция в изображениях Лапласа не может иметь равные между собой нули и полюса, так как в этом случае ее порядок может быть понижен путем сокращения числителя и знаменателя на общий множитель.
Передаточная функция системы управления в изображениях Лапласа может быть определена по ее передаточной функции в операторной форме
следующим образом:
Если передаточная функция содержит одинаковые нули и полюса, то элементарные множители, соответствующие этим корням в числителе и знаменателе, после подстановки
должны быть сокращены.
Возможно эта страница вам будет полезна:
Задача №1.1.
Определить передаточные функции звеньев, описываемых уравнениями:
Решение:
В символической форме эти уравнения записываются в виде
а их передаточные функции в операторной форме соответственно равны
Передаточные функции в изображениях Лапласа имеют вид
Как видим, передаточные функции в изображениях Лапласа рассматриваемых звеньев совпадают, хотя они описываются разными дифференциальными уравнениями и общие решения однородных уравнений, описывающие свободные движения систем, отличаются между собой.
Временные функции
Переходной функцией системы (звена) называют функцию, описывающую реакцию системы на единичное ступенчатое воздействие при нулевых начальных условиях. Переходную функцию обозначают . График переходной функции — кривую зависимости
от времени
— называют переходной или разгонной характеристикой.
Импульсной переходной или весовой функцией (функцией веса) называют функцию, описывающую реакцию системы (звена) на единичное импульсное воздействие при нулевых начальных условиях. Весовую функцию обозначают . График импульсной переходной функции называют импульсной переходной характеристикой. Переходную и импульсную переходную функции называют временными функциями, а их графики — временными характеристиками.
Передаточная функция в изображениях Лапласа есть преобразование Лапласа от весовой функции:
Весовая функция равна производной от переходной функции:
Если изображение временной функции имеет вид
, где
и
— полиномы, и степень
полинома
больше степени m полинома
, то
если нули полинома
— простые. Если какой-либо полюс
имеет кратность
, то ему соответствует слагаемое
Задача №1.2.
Определить переходную и весовую функции колебательного звена, т. е. звена с передаточной функцией
Решение:
Дифференциальное уравнение имеет вид
Для определения переходной функции нужно решить это уравнение при входном воздействии и нулевых начальных условиях:
Характеристическое уравнение имеет вид
и его корнями являются
или
Положив
общее решение однородного дифференциального уравнения можно записать в виде
Частное решение неоднородного уравнения . Поэтому общее решение неоднородного уравнения
Производная от этого решения
Начальные условия принимают вид
Отсюда
Поэтому для переходной и весовой функций имеем
или, после элементарных преобразований
где
Задача №1.3.
Определить переходную и весовую функции звена с передаточной функцией
Решение:
Передаточная функция является изображением Лапласа весовой функции
. Полюса передаточной функции
являются простыми, и весовую функцию
можно определить по формуле (1.2). В данном случае
и для весовой функции в соответствии с формулой (1.2) получаем
Так как то для изображения переходной функции имеем
В этом случае полюс имеет кратность
, а полюс
— простой. Поэтому слагаемое, соответствующее полюсу
, найдем по формуле (1.3), а слагаемое, соответствующее полюсу
, — по формуле (1.2). Согласно формуле (1.3) имеем
Так как
для слагаемого, соответствующего полюсу , имеем (см. (1.2))
Таким образом, переходная функция имеет вид
Частотные функции и характеристики
Функцию , которая получается из передаточной функции в изображениях Лапласа
при подстановке
, называют частотной передаточной функцией. Она является комплекснозначной функцией от действительной переменной
называемой частотой. Частотную передаточную функцию можно представить в виде
Если
Ha комплексной плоскости частотная передаточная функция определяет вектор
(рис. 1.1), длина которого равна
, а аргумент — углу
, образованному этим вектором с положительной действительной полуосью. Кривую, описываемую концом вектора
при изменении частоты от 0 до
или от —
до
, называют амплитудно-фазовой частотной характеристикой (АФЧХ).
АФЧХ, получаемую при изменении частоты от — до
, называют также диаграммой Найквиста. Модуль
называют амплитудной частотной функцией, ее график — амплитудной частотной характеристикой. Аргумент
называют фазовой частотной функцией, а его график (при изменении
от 0 до
) — фазовой частотной характеристикой.
Частотную передаточную функцию называют также амплитудно-фазовой частотной функцией. Ее действительную
и мнимую
части называют соответственно вещественной и мнимой частотной функцией, а их графики — кривые зависимостей
— вещественной и мнимой частотной характеристикой соответственно.
Кроме перечисленных частотных характеристик имеются логарифмические частотные характеристики (ЛЧХ): логарифмические амплитудные частотные характеристики (ЛАЧX) и логарифмические фазовые частотные характеристики (ЛФЧХ).
Функцию
называют логарифмической амплитудной частотной функцией, а график зависимости функции от логарифма частоты
— логарифмической амплитудной частотной характеристикой (ЛАЧХ).
При построении ЛАЧХ по оси абсцисс откладывают значение частоты в логарифмическом масштабе и при этом на отметке, соответствующей значению , записывают значение
; по оси ординат откладывают и записывают значение
.
Логарифмической фазовой частотной характеристикой (ЛФЧХ) называют график зависимости функции от логарифма частоты
. При ее построении по оси абсцисс, как и при построении ЛАЧХ, на отметке, соответствующей значению
, записывают значение
.
В ЛЧХ единицей функции является децибел, а единицей
— декада. Декадой называют интервал, на котором частота изменяется в 10 раз. При изменении частоты в 10 раз говорят, что частота изменилась на одну декаду.
Определенные трудности представляет вычисление фазовой частотной функции. Если эта функция по модулю не превышает , то она определяется по формуле
В общем случае нужно разложить числитель и знаменатель передаточной функции на элементарные множители и определять фазовую частотную функцию по правилу вычисления аргумента произведения и частного комплексных чисел.
Правило вычисления модуля и аргумента. При вычислении амплитудной и фазовой частотной функций полезно следующее правило вычисления модуля и аргумента произведения и частного комплексных чисел (функций).
1) Модуль произведения комплексных чисел равен произведению модулей сомножителей:
а аргумент — сумме аргументов сомножителей:
2) Модуль частного комплексных чисел (функций) равен отношению модулей
а аргумент — разности аргументов числителя и знаменателя:
Элементарные звенья и их характеристики. Так как произвольный полином можно разложить на простые множители, то передаточную функцию системы (звена)
всегда можно представить в виде произведения простых множителей и дробей вида
Здесь называется передаточным коэффициентом,
— постоянной времени и
— коэффициентом демпфирования.
Звенья, передаточные функции которых имеют вид простых множителей или дробей, называют элементарными звеньями. Их также называют типовыми.
Системы и звенья и их передаточные функции делятся на минимально-фазовые и неминимально-фазовые. Передаточная функция называется минимально-фазовой, если все ее нули (корни уравнения
) и полюса (корни уравнения
) располагаются в левой полуплоскости, и неминимально-фазовой, если хотя бы один нуль или полюс располагается в правой полуплоскости.
Система и звено называются минимально-фазовыми, если их передаточные функции являются минимально-фазовыми, и неминимально-фазовыми, если их передаточные функции являются неминимально-фазовыми.
Передаточные функции системы, не являющиеся ни минимально-фазовыми и ни неминимально-фазовыми, иногда называют нейтральными или маргинальными. Иначе говоря, передаточная функция называется маргинальной, если она имеет нуль или полюс на мнимой оси, но не имеет их в правой полуплоскости.
Тип звена определяется видом его передаточной функции. При этом если передаточные функции звеньев отличаются только на постоянный множитель, то их относят к одному и тому же типу. Поэтому при определении типа элементарных звеньев будем исходить из передаточных функций, получаемых из (1.4) умножением на константу (кроме первой).
Звено с передаточной функцией называется пропорциональным звеном, звено с передаточной функцией
— дифференцирующим звеном, звено с передаточной функцией
— интегрирующим звеном, звено с передаточной функцией
— форсирующим звеном (первого порядка), звено с передаточной функцией
— апериодическим звеном, звено с передаточной функцией —
— форсирующим звеном второго порядка, звено с передаточной функцией
— колебательным звеном.
Фазовые частотные функции минимально-фазовых и нейтральных звеньев с передаточными функциями, представляющими элементарный множитель первого порядка, по модулю не превышают и определяются по формуле
. В случае форсирующего звена второго порядка фазовая функция определяется по формуле
при частотах
, а при
— по формуле
Физический смысл частотных характеристик. При гармоническом входном воздействии в устойчивых системах после окончания переходного процесса выходная переменная также изменяется по гармоническому закону с той же частотой, но с другими амплитудой и фазой; амплитуда равна амплитуде входного сигнала, умноженной на модуль частотной передаточной функции, а сдвиг фазы — ее аргументу. Поэтому если система с передаточной функцией устойчива, то при входном воздействии
после окончания переходного процесса выходной сигнал
Здесь — постоянная амплитуда входного сигнала,
— начальный сдвиг фазы,
— частотная передаточная функция рассматриваемой системы,
Задача №1.4.
На вход системы подается сигнал Определить в установившемся режиме реакцию системы с передаточной функцией
Решение:
В данном случае частотная передаточная функция имеет вид
и
Поэтому
и соответственно
Задача №1.15.
На вход системы подается сигнал . Определить в установившемся режиме реакцию систем при передаточных функциях, приведенных в задании 1.14.
Асимптотические логарифмические амплитудные частотные характеристики. Логарифмические амплитудные частотные характеристики (ЛАЧХ) пропорционального, дифференцирующего и интегрирующего звеньев являются прямыми и их легко построить. Построение ЛАЧХ других элементарных звеньев требует трудоемких вычислений. Поэтому на практике часто ограничиваются построением приближенных асимптотических ЛАЧХ.
При построении асимптотической ЛАЧХ апериодического звена в выражении
при под корнем пренебрегают слагаемым
, меньшим единицы, а при
— единицей. Поэтому уравнение асимптотической ЛАЧХ имеет вид
При построении асимптотической ЛАЧХ колебательного звена в выражении
при под корнем оставляют только единицу, а при
— только наибольшее слагаемое
. Поэтому уравнение асимптотической ЛАЧХ имеет вид
Аналогично поступают при построении асимптотических ЛАЧХ форсирующих звеньев. Частоты, на которых асимптотические ЛАЧХ претерпевают излом, называются сопрягающими частотами.
Для построения ЛАЧХ и ЛФЧХ звена с произвольной дробно-рациональной передаточной функцией нужно ее числитель и знаменатель разложить на элементарные множители и представить
в виде произведения передаточных функций элементарных звеньев
или в виде
где представляет собой отношение произведений элементарных множителей 1-го и 2-го порядка с единичным передаточным коэффициентом, т. е. множителей вида
Из (1.5) имеем:
Из (1.7) следует, что для построения ЛАЧХ произвольного звена достаточно построить ЛАЧХ элементарных звеньев, на которые оно разлагается, а затем их геометрически сложить. Однако для построения асимптотических ЛАЧХ можно использовать несколько иное, более простое правило. Проиллюстрируем это сначала на частном примере.
Возможно эта страница вам будет полезна:
Задача №1.5.
Построить асимптотическую ЛАЧХ для звена с передаточной функцией
Решение:
Преобразуем передаточную функцию к виду
Логарифмическая амплитудная частотная функция
Вычислим сопрягающие частоты и пронумеруем их в порядке возрастания:
Здесь и
— сопрягающие частоты апериодического, форсирующего и колебательного звеньев соответственно.
Напомним, что при построении асимптотических ЛАЧХ при частотах, меньших сопрягающей частоты, под корнем оставляют только единицу (остальными членами пренебрегают), при частотах, больших сопрягающей частоты, — член с наивысшей степенью . Поэтому в рассматриваемом примере при
имеем
Это уравнение прямой, которая проходит через точку с координатами и
с наклоном —20 дб/дек. Прямая имеет наклон -20 дб/дек (20 дб/дек) — это означает, что при увеличении частоты на декаду (т. е. в 10 раз)
уменьшается (увеличивается) на 20 дБ (рис. 1.2, а). Первая асимптота заканчивается на первой сопрягающей частоте (рис. 1.2, б).
При аналогично имеем
Это уравнение второй асимптоты. Ее наклон по отношению к первой асимптоте изменяется на —20 дб/дек и обусловливается апериодическим звеном, т. е. множителем первого порядка в знаменателе рассматриваемой передаточной функции. Вторую асимптоту проводят от конца первой асимптоты до второй сопрягающей частоты под наклоном -40 дб/дек.
При
Это уравнение третьей асимптоты. Ее наклон по отношению ко второй асимптоте изменяется на 20 дб/дек и обусловливается форсирующим звеном, т. е. множителем первого порядка в числителе. Третью асимп-
тоту проводят от конца второй асимптоты до третьей сопрягающей частоты под наклоном —20 дб/дек. При
Это уравнение последней, четвертой, асимптоты. Ее наклон изменяется по отношению к третьей асимптоте на -40 дб/дек и обусловливается множителем второго порядка в знаменателе.
Правило построения асимптотических ЛАЧХ
1) Пользуясь представлением (1.6), вычислить 20 : и сопрягающие частоты
, которые следует пронумеровать в порядке возрастания:
2) На оси абсцисс отметить сопрягающие частоты, а на координатной плоскости — точку (1,20). Построить первую асимптоту — прямую под наклоном —
20 дБ/дек, проходящую через отмеченную точку на координатной плоскости. Первая асимптота заканчивается на первой сопрягающей частоте
.
3) Построить вторую асимптоту, которая начинается с конца первой асимптоты и проводится до второй сопрягающей частоты Его наклон изменяется на ±20 дБ/дек или ±40 дБ/дек в зависимости от того, обусловливается ли элементарным множителем первого или второго прядка соответственно. Принимается положительный знак, если указанный множитель находится в числителе, и отрицательный знак, если этот множитель находится в знаменателе.
4) Построить остальные асимптоты, которые строятся аналогично второй асимптоте: -я асимптота начинается с конца предыдущей (
— 1)-й асимптоты и проводится до сопрягающей частоты
. Ее наклон определяется сопрягающей частотой
.
Последняя асимптота представляет собой прямую, которая начинается в конце асимптоты, соответствующей последней сопрягающей частоте, и уходит в бесконечность.
Возможно эта страница вам будет полезна:
Задача №1.6.
Построить асимптотическую ЛАЧХ звена с передаточной функцией
Решение:
Преобразуем передаточную функцию к виду
1) = 0. Вычислим 20
и сопрягающие частоты:
Проводим через точку с координатами (1, 20) первую асимптоту под наклоном 0 дБ/дек (т. е. параллельно оси абсцисс) до первой сопрягающей частоты = 0,1 (рис. 1.3, а).
Так как первая сопрягающая частота обусловлена множителем первого порядка
, расположенным в знаменателе, наклон второй асимптоты изменяется на -20 дБ/дек. Поэтому вторую асимптоту проводим от конца первой асимптоты до сопрягающей частоты
= 1 под наклоном -20 дБ/дек.
Сопрягающая частота обусловлена элементарным множителем
, расположенным в числителе. Поэтому наклон третьей асимптоты отличается от наклона второй на 20 дБ/дек и составляет 0 дБ/дек. Третью асимптоту проводим от конца второй асимптоты до сопрягающей частоты
= 10.
Сопрягающая частота обусловлена элементарным множителем
. расположенным в знаменателе. Поэтому наклон четвертой асимптоты отличается от наклона третьей асимптоты на -40 дБ/дек. Последнюю асимптоту проводим от конца третьей асимптоты до бесконечности.
2) = — 1. Значения 20
и сопрягающих частот те же, что и в предыдущем случае. Первую асимптоту проводим через точку с координатами (1, 20) с наклоном —
20 дБ/дек = 20 дБ/дек до первой сопрягающей частоты (рис. 1.3, б). Все последующие асимптоты строятся так же, как и в предыдущем случае.
Структурные схемы
Структурной схемой системы управления называют графическое представление ее математической модели в виде соединений звеньев, изображаемых в виде прямоугольников или круга (для сумматора), с указанием входных и выходных переменных.
Обычно внутри прямоугольника указывается условное обозначение оператора изображаемого им звена, а сам оператор в виде передаточной функции или дифференциального уравнения задается вне структурной схемы.
Преобразование структурных схем.
Последовательное соединение. Так называется соединение, при котором выход предыдущего звена является входом последующего (рис. 1.4, а). При последовательном соединении передаточные функции отдельных звеньев перемножаются и при преобразовании структурных схем цепочку из последовательно соединенных звеньев можно заме-
нить одним звеном с передаточной функцией (рис. 1.4, б).
Параллельное соединение. Так называется соединение, при котором на вход всех звеньев подается одно и то.же воздействие, а их выходные переменные складываются (рис. 1.5, а). При параллельном
соединении звеньев передаточные функции складываются и при преобразовании их можно заменить одним звеном с передаточной функцией
Если выход какого-либо звена поступает на сумматор с отрицательным знаком, то передаточная функция этого звена складывается с отрицательным знаком, т. е. вычитается.
Обратное соединение или звено, охваченное обратной связью. Так называется соединение двух звеньев, при котором выход звена прямой цепи подается на вход звена обратной связи, выход которого складывается с входом первого звена (рис. 1.6, а). Если сигнал обратной
связи (выход звена обратной связи) вычитается (т. е. складывается с отрицательным знаком), то обратная связь называется отрицательной; в противном случае — положительной. Когда передаточная функция звена обратной связи равна единице обратное соединение изображается так, как показано на рис. 1.6, б.
При размыкании обратной связи перед сумматором получаем последовательное соединение, передаточная функция которого равна
Эта передаточная функция называется передаточной функцией разомкнутой цепи.
Передаточную функцию
в которой учитывается передаточная функция сумматора по входу обратной связи, будем называть передаточной функцией контура. Здесь — передаточная функция сумматора по входу обратной связи, она равна -1 (минус единице) при отрицательной обратной связи (перед соответствующим входом стоит знак минус) и 1 (плюс единице) при положительной обратной связи.
Передаточная функция при обратном соединении равна и при преобразовании обратное соединение заменяется одним звеном с указанной передаточной функцией (рис. 1.6, в).
Перенос сумматора. При переносе сумматора по ходу сигнала добавляется звено с передаточной функцией, равной передаточной функции звена, через которое переносится сумматор (рис. 1.7, а).
При переносе сумматора против хода сигнала добавляется звено с передаточной функцией, равной обратной передаточной функции звена, через которое переносится сумматор (рис. 1.7, б).
При переносе сумматора участок цепи, через который он переносится, становится неэквивалентным. Поэтому при преобразовании структурных схем нельзя переносить сумматор через точку съема сигнала.
Перенос узла. При переносе узла по ходу сигнала добавляется звено с передаточной функцией, равной обратной передаточной функции звена, через которое переносится узел (рис. 1.8, а).
При переносе узла против хода сигнала добавляется звено с передаточной функцией, равной передаточной функции звена, через которое переносится узел (рис. 1.8, б).
Перестановка сумматоров. Сумматоры можно переставлять местами и объединять. Перестановка двух сумматоров соответствует переносу одного сумматора через другой и подчиняется правилу переноса сумматора через звено.
Сумматор 1 (рис. 1.9) переносится через сумматор 2 по направлению распространения сигнала, а сумматор 2 через сумматор 1 против направления распространения сигнала.
Но так как передаточная функция сумматора по каждому входу равна 1 или -1, то и передаточная функция звена, которое добавляется при переносе сумматора, независимо от направления переноса равна 1 или -1. Поэтому если сумматор переносится через другой сумматор вдоль входа со знаком плюс, добавляется звено с передаточной функцией 1, т. е. в действительности ничего не добавляется (рис. 1.9, а); если сумматор переносится вдоль входа со знаком минус, то добавляется звено с передаточной функцией -1, т.е. знак по входу, куда должно быть добавлено звено, меняется на обратный (рис. 1.9, б).
Перестановка узлов. Узлы можно переставлять местами и объединять.
Вычисление передаточной функции одноконтурной системы.
Замкнутая система называется одноконтурной, если при ее размыкании в какой-либо точке замкнутого контура получается система без параллельных и обратных соединений (рис. 1.10).
Цепь по ходу сигнала от точки приложения входной переменной до точки съема выходной переменной называется прямой цепью. Передаточная функция прямой цепи равна произведению передаточных функций звеньев, входящих в эту цепь, включая и сумматоры. Передаточная функция контура
равна произведению передаточных функций всех звеньев, входящих в замкнутый контур, включая сумматоры. Передаточная функция сумматора по входу со знаком плюс равна плюс единице, а по входу со знаком минус — минус единице.
Правило вычисления передаточной функции замкнутой одноконтурной системы: передаточная функция одноконтурной системы относительно внешнего воздействия (входа) и выхода
равна передаточной функции прямой цепи, деленной на единицу минус передаточная функция контура:
Возможно эта страница вам будет полезна:
Задача №1.7.
Определить передаточные функции системы (рис. 1.10) относительно входа
и выхода
и
относительно входа
и выхода
.
Решение:
Прямая цепь системы (см. рис. 1.10) относительно входа и выхода
представляет последовательное соединение двух сумматоров и звеньев с передаточными функциями
и
. Входы сумматоров в этой цепи имеют знак плюс и их передаточные функции равны единице. Поэтому передаточная функция прямой цепи
Прямая цепь относительно входа и выхода
представляет последовательное соединение двух сумматоров и звеньев с передаточными функциями
и
. Вход первого сумматора имеет знак плюс, вход второго сумматора — знак минус и их передаточные функции равны 1 и —1 соответственно. Поэтому в этом случае передаточная функция прямой цепи
Искомые передаточные функции имеют вид
Вычисление передаточной функции многоконтурной системы.
Замкнутая система называется многоконтурной, если при ее размыкании в какой-либо точке замкнутого контура получается система, содержащая параллельное и/или обратное соединение.
Многоконтурная система не имеет перекрестных связей, если любые два контура, образованные параллельными или обратными соединениями, не имеют общих участков (рис. 1.11, а) или, если какие-либо два контура имеют общий участок, то один из них вложен внутрь другого (рис. 1.11, б).
Многоконтурная система имеет перекрестные связи, если она содержит два контура, которые имеют общий участок, и при этом ни один из них не вложен внутрь другого (рис. 1.11, в).
Порядок вычисления передаточной функции многоконтурной системы следующий:
1) путем переноса узлов и сумматоров нужно освободиться от перекрестных связей;
2) используя правила преобразования параллельных и обратных соединений, нужно преобразовать многоконтурную систему в одноконтурную;
3) по правилу вычисления передаточной функции одноконтурной системы определить искомую передаточную функцию.
При преобразовании структурной схемы нужно позаботится о том, чтобы не исчезли точки съема переменных, относительно которых ищутся передаточные функции, или чтобы эти точки не оказались на неэквивалентном участке (т. е. не следует переносить сумматор через эти точки).
Задача №1.8.
Определить передаточные функции и
системы управления, представленной на рис 1.12, а.
Решение:
Сначала освободимся от перекрестных связей. Для этого перенесем сумматор 3 против хода сигнала через звено и сум-
матор 2. То же самое проделаем с сумматором 4 (рис. 1.12, б). Далее, заменив параллельное соединение звеном с передаточной функцией
и обратное соединение звеном с передаточной функцией
получим одноконтурную систему (рис. 1.12, в). Из последней схемы по правилу вычисления передаточной функции одноконтурной системы находим
При вычислении передаточных функций многоконтурных систем с перекрестными связями во многих случаях целесообразно, а иногда
и необходимо сначала предварительно упростить схему, используя правила преобразования параллельных и обратных соединений, затем освободиться от перекрестных связей.
Граф системы управления
Граф системы управления состоит из дуг и вершин. Дуга соответствует звену и на схеме изображается отрезком линии со стрелкой, указывающей направление распространения сигнала. Дуга начинается и кончается в вершине.
Вершина на схеме изображается кружком и определяет переменную. Если к вершине подходит одна дуга, то она определяет выходную величину дуги (рис. 1.17, а), если же в вершину входят несколько дуг, то она соответствует сумме выходных переменных этих дуг (рис. 1.17, б).
Начальная вершина дуги определяет ее входную переменную (рис. 1.17, в). Вершина графа, имеющая только выходящие из нее дуги, определяет внешнее воздействие и называется входной вершиной графа.
Последовательность дуг (не обязательно разных), для которых конечная вершина
дуги
является начальной вершиной дуги
называется ориентированным маршрутом или ормаршрутом. Ормаршрут называется замкнутым, если конечная вершина дуги
совпадает с начальной вершиной дуги
, и незамкнутым в противном случае.
Ормаршрут, в котором все дуги разные, называется путем от начальной вершины к конечной вершине
, если он не замкнут, и контуром, если он замкнут (
и
совпадают). Путь и контур называют простыми, если все вершины
различны. Простой путь также называют прямым путем.
Два контура называются несоприкасающимися, если они не имеют общих вершин. Три, четыре и т.д. контура называются несоприкасающимися, если любая пара из этих контуров является несоприкасающейся.
Граф системы управления можно построить по структурной схеме. Для этого нужно произвести следующее (рис. 1.18):
1) сумматор с выходной переменной заменить вершиной
;
2) звено с передаточной функцией заменить дугой
если выходная переменная подается на сумматор по отрицательному входу, то указанное звено заменить дугой —
;
3) каждой переменной, в том числе переменной, соответствующей внешнему воздействию, сопоставить свою вершину.
Формула Мейсона. Определителем графа (подграфа) называется передаточная функция , равная
Здесь в первой сумме — передаточная функция
-го простого контура, равная произведению передаточных функций дуг, входящих в этот контур, и суммирование производится по всем простым контурам; во второй сумме
— произведение передаточных функций
-го и
-го простых контуров и суммирование производится по всем несоприкасающимся парам контуров; в третьей сумме
— произведение передаточных функций
-го,
-го и
-го простых контуров и суммирование производится по всем несоприкасающимся тройкам контуров и т. д.
Подграфом -го прямого пути называется подграф, который получается из исходного графа отбрасыванием всех дуг и вершин
-го пути, а также всех дуг, начинающихся или кончающихся на вершинах этого пути.
Передаточная функция системы управления относительно входа и выхода
определяется следующим образом:
где — определитель графа системы управления;
— передаточная функция
-го прямого пути от начальной вершины
к конечной вершине
;
— общее число таких прямых путей;
— определитель подграфа
-го прямого пути.
Задача №1.9.
Построить граф и по теореме Мейсона определить передаточную функцию системы (рис. 1.19, а).
Решение:
Граф системы управления представлен на рис. 1.19, б. От вершины до вершины
имеются четыре прямых пути. Передаточные функции этих путей равны
Подграф 1-го пути состоит из вершин и
, 2-го пути — из вершин
и
; подграф 3-го пути есть пустой граф, подграф 4-го пути состоит из вершины
. И так как все они не имеют контуров, их определители равны единице:
Граф системы управления имеет четыре простых контура. Их передаточные функции имеют вид
Несоприкасающихся пар контуров нет. Поэтому определитель графа имеет вид
Для искомой передаточной функции получаем
Математическое описание некоторых технических устройств
В общем случае функциональная схема системы автоматического управления имеет вид, представленный на рис. 2.1, где приняты следующие обозначения: УУ — управляющее устройство, включающее
в себя ЗУ — задающее устройство, вырабатывающее задающий сигнал ; СУ — сравнивающее устройство, вырабатывающее сигнал ошибки
УПУ — усилительно-преобразовательное устройство, включающее в себя помимо усилителя и преобразователь или корректирующее устройство, которое на основе сигнала ошибки и измеренного возмущения
вырабатывает управляющее воздействие
; ИУ — исполнительное устройство, непосредственно воздействующее на объект управления ОУ; ЧЭ1 и ЧЭ2 — чувствительные элементы (датчики), измеряющие управляемую переменную
и возмущение
и при необходимости преобразующие их в иную физическую переменную (например, механическую или тепловую в электрическую); ОУ — объект управления.
В данной главе рассматриваются задачи, связанные с математическим описанием (дифференциальными уравнениями и передаточными функциями) некоторых технических устройств, используемых в системах автоматического управления (САУ) в качестве упомянутых выше элементов.
Чувствительные элементы — датчики
Датчики линейных и угловых перемещений. В САУ для измерения линейных и угловых перемещений используются линейные и вращающиеся потенциометрические датчики (ПД). Для измерения угловых перемещений используются вращающиеся трансформаторы (ВТ) и сельсины (С). На этих элементах выполняют также и сравнивающие устройства (СУ). Принцип действия этих устройств, их схемы и основные характеристики рассматриваются в довольно обширной литературе [1, 2, 5-8, 10]. Упомянутые выше потенциометрические датчики, вращающиеся трансформаторы и сельсины при исследовании динамики считаются безынерционными звеньями с передаточной функцией
где — передаточный коэффициент датчика.
Для потенциометрических датчиков и вращающихся трансформаторов коэффициент определяется крутизной статической характеристики, но следует иметь в виду, что для вращающихся трансформаторов это верно при малых углах, иначе необходимо учитывать нелинейность характеристики. Передаточный коэффициент сельсина, работающего в трансформаторном режиме, также рассчитывается по крутизне характеристики, определяемой по следующей формуле:
Так, например, передаточная функция ПД типа ПП
передаточная функция ВТ типа ВТ-5
передаточная функция сельсина типа СГСМ-1
Для измерения угловой скорости используют тахогенераторы (ТГ) постоянного и переменного тока. Строго говоря, по динамическим свойствам их можно отнести к апериодическому звену второго порядка с передаточной функцией
где — напряжение на выходе ТГ,
— электрическая постоянная времени,
— индуктивность обмотки якоря,
— активное сопротивление обмотки якоря,
— электромеханическая постоянная времени. Но поскольку якорь ТГ соединен с валом двигателя, скорость которого он измеряет, его момент инерции учитывается при расчете электромеханической постоянной времени двигателя в суммарном моменте инерции. Учитывая, что
, ТГ можно считать безынерционным звеном с передаточной функцией
. где
— передаточный коэффициент ТГ, определяемый крутизной статической характеристики
.
Например, передаточная функция ТГ типа ТП-75
а для ТГ типа ДГ-3ТА
Основные характеристики некоторых типов ТГ приведены в [2, 10].
Датчики температуры. Для измерения температуры в системах автоматического управления используются электротепловые датчики: термопары (ТП) и термосопротивления (ТС) [5, 6, 8, 10]. Датчики этого типа с точки зрения динамики являются апериодическим (инерционным) звеном первого порядка с передаточной функцией
где — постоянная времени термодатчика, которая колеблется для некоторых типов датчиков от долей секунды до нескольких минут,
— передаточный коэффициент термодатчика, который определяется крутизной статической характеристики.
Усилители
В САУ используются все известные типы усилителей: электрические, гидравлические и пневматические. В качестве электрических используются электронные (ЭУ) (полупроводниковые, тиристорные), магнитные (МУ) и электромашинные (ЭМУ).
Электронные усилители. Электронные усилители можно считать безынерционным звеном с передаточной функцией
так как их постоянная времени мала по сравнению с постоянными времени электромеханических элементов системы. Коэффициент усиления по напряжению рассчитывается как отношение выходного напряжения усилителя
ко входному напряжению
.
Магнитные усилители. Наибольшее распространение в САУ получила схема двухтактного реверсивного МУ [5, 8-10]. По динамическим свойствам МУ этого типа эквивалентен апериодическому звену с передаточной функцией
Для увеличения коэффициента усиления используют внутреннюю обратную связь. Постоянная времени Тму для МУ с положительной обратной связью рассчитывается по следующей формуле:
где — частота напряжения питания в Гц,
— активное сопротивление нагрузки,
— общее активное сопротивление цепи управления усилителя с учетом сопротивления источника управляющего сигнала в Ом,
и
— число витков рабочей обмотки и обмотки управления соответственно,
— коэффициент положительной обратной связи. Коэффициент усиления по напряжению
вычисляется по формуле
где — ток в нагрузке,
— ток в обмотке управления.
Магнитные усилители рассчитываются для каждого отдельного случая, серийно промышленностью не выпускаются.
Электромашинные усилители. ЭМУ используются в САУ в случае наличия источника механической энергии (например, дизель и т.п.). Их применяют для управления двигателем постоянного тока, когда требуется высокий коэффициент усиления по мощности. Известны различные конструкции ЭМУ [3, 5, 6, 8, 10]. ЭМУ с поперечным полем описывается передаточной функцией апериодического звена второго порядка
где — коэффициент усиления ЭМУ, равный
и
— активные сопротивления обмотки управления и поперечной короткозамкнутой обмотки соответственно,
— коэффициент пропорциональности между ЭДС в поперечной обмотке и током управления,
—коэффициент пропорциональности между выходной ЭДС и током в поперечной обмотке,
— постоянная времени цепи управления,
— индуктивность обмотки управления,
— постоянная времени поперечной цепи,
— индуктивность поперечной обмотки. Обычно
, их значения колеблются от сотых до десятых долей секунды. Коэффициент усиления по мощности для этого типа ЭМУ достигает
.
Возможно эта страница вам будет полезна:
Исполнительные устройства и объекты управления
Двигатели постоянного тока. Двигатель постоянного тока с независимым возбуждением может быть представлен структурной схемой, приведенной на рис. 2.2, где —передаточная функция относительно управляющего воздействия
— передаточная функция относительно возмущения — момента нагрузки
.
Когда выходом является угловая скорость, передаточная функция двигателя по управляющему воздействию
и по возмущению
Здесь
передаточный коэффициент двигателя по управлению, — постоянная, зависящая от потока возбуждения и конструкции двигателя,
— скорость холостого хода,
— номинальная скорость,
— электрическая постоянная времени якоря,
— индуктивность обмотки якоря,
— активное сопротивление обмотки якоря,
— электромеханическая постоянная времени,
— приведенный к валу двигателя суммарный момент инерции вращающихся частей,
передаточный коэффициент двигателя по возмущению (моменту нагрузки) — постоянная, зависящая, как и
, от потока возбуждения и конструкции двигателя,
— номинальное напряжение управления.
— пусковой момент.
Для большинства двигателей выполняется неравенство . Поэтому при расчете динамики САУ часто полагают
. При этом передаточные функции двигателя по управляющему воздействию
и по возмущению
соответственно принимают вид
Если за выходную величину двигателя принять угол поворота вала , то передаточные функции по управляющему воздействию
и по возмущению
имеют вид
Задача №2.1.
Определить передаточные функции двигателя типа ДПМ-20-Н1/Н2-01.
Решение:
Для двигателя данного типа
Приведем единицы измерения параметров двигателя к системе СИ:
Скорость холостого хода двигателя
Рассчитаем передаточные коэффициенты двигателя по управлению и по возмущению:
Тогда получим передаточные функции двигателя
Асинхронные двигатели. Наиболее распространен индукционный двухфазный двигатель [1, 3, 8-10]. В динамическом отношении асинхронный двигатель рассматривается относительно угловой скорости как апериодическое звено и по управляющему воздействию и по возмущению
:
где параметры двигателя вычисляются по следующим формулам:
где — пусковой ток ротора, равный току, потребляемому от сети,
— момент инерции ротора.
Задача №2.2.
Определить передаточные функции асинхронного двигателя типа АД-32Б. Технические характеристики двигателя этого типа:
Решение:
Угловая скорость двигателя при холостом ходе
Рассчитаем передаточные коэффициенты двигателя:
Тогда передаточные функции двигателя
Генератор постоянного тока. Генератор постоянного тока описывается дифференциальным уравнением первого порядка и он эквивалентен апериодическому звену [6]:
где — выходное и входное напряжения генератора,
— передаточный коэффициент по управляющему воздействию,
— активное сопротивление обмотки возбуждения, тор — константа, определяющая зависимость между ЭДС генератора
и током возбуждения
— постоянная времени генератора,
— индуктивность обмотки возбуждения.
Передаточная функция генератора относительно возмущения ( — тока якоря)
где — активное сопротивление цепи якоря.
Корректирующие элементы
При синтезе САУ для обеспечения ее устойчивости и требуемых показателей качества используют корректирующие элементы, в качестве которых применяют пассивные и активные четырехполюсники.
Пассивные четырехполюсники. Пассивные четырехполюсники представляют собой схемы из резисторов, конденсаторов и индуктив-ностей [5-7].
При вычислении передаточных функций четырехполюсников удобно воспользоваться операторными сопротивлениями: омическим индуктивным
и емкостным
. При этом пассивные четырехполюсники можно рассчитывать как схемы, составленные из одних омических сопротивлений. Общая схема пассивного четырехполюсника показана на рис. 2.3, где
и
— операторные сопротивления.
Передаточную функцию такого четырехполюсника можно записать следующим образом:
Задача № 2.3.
Рассчитать передаточную функцию четырехполюсника, показанного на рис. 2.4.
Решение:
В данном случае
Поэтому
Если соединить последовательно два пассивных четырехполюсника через разделительный усилитель (рис. 2.5), то передаточная функция этой цепи
где — коэффициент усиления усилителя,
и
— передаточные функции четырехполюсников, включенных на входе и выходе усилителя.
Эта формула справедлива при условии, что входное сопротивление усилителя достаточно велико.
Активные четырехполюсники постоянного тока. В таких четырехполюсниках используются операционные усилители (УПТ) с высоким коэффициентом усиления [4. 6, 7, 9]. Общая схема активного четырехполюсника показана на рис. 2.6. Передаточная функция такого элемента
Задача №2.4.
Рассчитать передаточную функцию активного четырехполюсника, показанного на рис. 2.7.
Решение:
В данном случае
Поэтому
Сравнивающие устройства (СУ)
На рис. 2.8, а показана схема СУ, выполненная на линейных потенциометрах и
, а на рис. 2.8, б — на кольцевых потенциометрах
и
.
В обеих схемах сигнал ошибки при равенстве задающего сигнала
и сигнала обратной связи
.
На рис. 2.9 показана мостовая схема СУ.
В частном случае в плечи моста могут быть включены активные сопротивления и термосопротивление
. Если выполняется условие равновесия моста
, то сигнал ошибки
. В общем случае в плечи моста могут быть включены, помимо активных сопротивлений, индуктивности и емкости.
Схема СУ может быть выполнена и на сельсинах, и на вращающихся трансформаторах. Принципиальная схема таких устройств может быть показана так, как на рис. 2.10.
В качестве задающего (ЗУ) и приемного (ПУ) устройств могут использоваться и сельсины (СД-сельсин-датчик, СП-сельсин-приемник) и вращающиеся трансформаторы (ВТ-1 и ВТ-2). Сигнал ошибки при равенстве углов поворота задающей оси
и приемной оси
.
На рис. 2.11 показана схема СУ, выполненная на операционном усилителе (активном четырехполюснике). Сигнал ошибки при равенстве напряжений
.
Для всех приведенных выше схем СУ (рис. 2.8-2.11) структурная схема показана на рис. 2.12, где — входной сигнал,
— сигнал обратной связи,
— сигнал ошибки,
— передаточный коэффициент СУ.
Устойчивость непрерывных систем управления
Основное условие устойчивости: для того чтобы непрерывная система управления была устойчива, необходимо и достаточно, чтобы все корни ее характеристического уравнения имели отрицательную вещественную часть.
На комплексной плоскости корни, имеющие отрицательную вещественную часть, располагаются в левой полуплоскости и поэтому называются левыми, корни, имеющие положительную вещественную часть, располагаются в правой полуплоскости и называются правыми, а корни, расположенные на мнимой оси, — нейтральными. Поэтому основное условие устойчивости можно также сформулировать еще так: для того чтобы система была устойчива, необходимо и достаточно, чтобы все корни характеристического уравнения были левыми.
Необходимое условие устойчивости. Для того чтобы система была устойчива, необходимо, чтобы все коэффициенты ее характеристического уравнения
были строго одного знака:
Характеристическое уравнение. Характеристический полином (левая часть характеристического уравнения
= 0) получается из собственного оператора
простой заменой оператора
на комплексную переменную
. Если дано уравнение системы управления в символической форме, то дифференциальный оператор при выходной переменной и будет собственным оператором. Если дана передаточная функция, то собственный оператор (с точностью до обозначения переменной) совпадает с ее знаменателем.
При исследовании замкнутой системы (рис. 3.1, а) нет необходимости находить ее передаточную функцию, если известна передаточная функция разомкнутой системы (рис. 3.1, б).
Ее собственный оператор равен сумме операторов числителя и знаменателя передаточной функции разомкнутой системы:
Алгебраические критерии устойчивости
При проведении исследования устойчивости с помощью алгебраических критериев следует, прежде всего, записав характеристическое уравнение, проверить выполнение необходимого условия устойчивости, так как его проверка не требует никаких вычислений и в то же время при его невыполнении не надо проводить дальнейших исследований.
Определители Гурвица. Из коэффициентов характеристического полинома
составим определитель -го порядка
который строится следующим образом. На главной диагонали выписываются элементы . Затем, двигаясь от этих элементов вверх, помещаются коэффициенты в порядке возрастания индексов, вниз — в порядке их убывания. Например, при построении
-го столбца, двигаясь от элемента
вверх, записываются коэффициенты
вниз — коэффициенты
При этом, если индекс превышает
или принимает отрицательное значение, то вместо соответствующего коэффициента записывают нуль. Определитель
и его главные миноры
называют определителями Гурвица.
Критерий Гурвица (Hurwitz, 1895). Для того чтобы система была устойчива, необходимо и достаточно, чтобы все определители Гурвица, составленные из коэффициентов ее характеристического уравнения, при были больше нуля:
Критерий Льенара—Шипара (Lienard, Chipard, 1914). При выполнении необходимого условия 0, для устойчивости системы управления необходимо и достаточно, чтобы все ее определители Гурвица с четными индексами или все ее определители Гурвица с нечетными индексами
были положительными’.
Для уменьшения вычислений целесообразно при нечетном использовать условие (3.1а), а при четном
— условие (3.16).
Выпишем необходимые и достаточные условия устойчивости для = 1.2,3:
Задача №3.1.
Передаточная функция разомкнутой системы имеет вид
Исследовать устойчивость разомкнутой и замкнутой систем.
Решение:
Характеристический полином разомкнутой системы имеет вид
Все коэффициенты больше нуля и определитель = 0,5 -4 — 1 — I — = 1 > 0. Поэтому по критерию Льенара—Шипара разомкнутая система устойчива.
Характеристический полином замкнутой системы
Все коэффициенты этого полинома при обоих значениях положительны, а определитель
при
= 0,5
а при = 2
Следовательно, замкнутая система при = 0,5 устойчива, а при
= 2 неустойчива.
Частотные критерии устойчивости
Критерий Найквиста (Nyqvist, 1932). Для того чтобы замкнутая система (с отрицательной обратной связью) была устойчива, необходимо и достаточно, чтобы амплитудно-фазовая частотная характеристика (АФЧХ) разомкнутой системы охватывала раз в положительном направлении точку
, где
— число правых корней характеристического уравнения разомкнутой системы.
Если разомкнутая система устойчива , для устойчивости замкнутой системы необходимо и достаточно, чтобы АФЧХ разомкнутой системы не охватывала точку
.
Задача №3.2.
Исследовать устойчивость замкнутой системы, если передаточная функция разомкнутой системы имеет вид
Решение:
Частотные передаточные функции и вещественные и мнимые частотные функции имеют вид:
Для построения АФЧХ нужно определить координаты точек ее пересечения с осями координат и соединить эти точки плавной кривой. Необходимые расчетные данные приведены в таблице 3.1. На основе этих данных построены АФЧХ (рис. 3.2).
Расчетные данные к примеру 3.4
Расчетные данные к примеру 3.2. б)
В случае а) замкнутая система устойчива, так как = 1 и АФ-ЧХ охватывает точку
1/2 раз в положительном направлении (рис. 3.2, а). В случае б) замкнутая система неустойчива, так как разомкнутая система устойчива (
= 0), а АФЧХ охватывает точку
(рис.3.2, б).
Случай наличия нулевых корней. Если характеристическое уравнение разомкнутой системы имеет нулевые корни, т. е. ее передаточная функция может быть представлена в виде
то АФЧХ при уходит в бесконечность (рис. 3.3). В этом случае АФЧХ дополняются дугой —
окружности большого радиуса (на рис. 3.3 — пунктирная линия). И для устойчивости замкнутой системы должна охватывать
раз или при
= 0 не охватывать точку
дополненная АФЧХ.
Устойчивость систем с чистым запаздыванием
Рассмотрим замкнутую систему управления, передаточная функция разомкнутой системы которой имеет вид
где — полиномы степени
и
соответственно
. Для исследования устойчивости такой системы может быть использован критерий Найквиста.
Для того чтобы замкнутая система, передаточная функция которой в разомкнутом состоянии имеет вид (3.2), была устойчива, необходимо и достаточно, чтобы АФЧХ разомкнутой системы охватывала точку в положительном направлении
/2 раз, где
— число правых нулей характеристического полинома разомкнутой системы
.
Замкнутая система со звеном чистого запаздывания, будучи устойчивой при малом , с ростом
ее АФЧХ в разомкнутом состоянии может приближаться к точке
и при некотором значении
пересечь ее, и замкнутая система окажется на границе устойчивости. Запаздывание
называют критическим.
Частотная передаточная функция и амплитудная и фазовая частотные функции разомкнутой системы имеют вид
где
Отсюда видно, что появление чистого запаздывания не меняет модуль, а только вносит дополнительный отрицательный фазовый сдвиг , что приводит к закручиванию АФЧХ (рис. 3.4).
Критическое запаздывание находится из условий
Решив эту систему, найдем критическое запаздывание и частоту , которая называется критической частотой.
Задача №3.3.
Определить критическое запаздывание и критическую частоту для системы, у которой передаточная функция в разомкнутом состоянии
Решение:
Без запаздывания замкнутая система устойчива. Условие (3.3) принимает вид
Отсюда получаем
Определение области устойчивости
Структура системы определяется составом элементов (звеньев) и связями между ними. При заданной структуре какие-либо параметры могут быть не фиксированными, т. е. их можно изменять. Такие параметры называют варьируемыми. Областью устойчивости в пространстве параметров называют множество всех значений варьируемых параметров, при которых система устойчива.
Если существует область устойчивости в пространстве параметров, то система называется структурно устойчивой (относительно заданных варьируемых параметров). В противном случае система называется структурно неустойчивой (относительно заданных варьируемых параметров).
Область устойчивости можно определить с помощью алгебраических критериев устойчивости. Рассмотрим это на примере.
Задача №3.4.
Передаточная функция разомкнутой системы . Определить область устойчивости замкнутой системы на плоскости параметров
.
Решение:
Характеристический полином замкнутой системы имеет вид
По критерию Льенара—Шипара имеем
Очевидно, эти неравенства будут выполнены, если
Эта система неравенств определяет область устойчивости.
Робастная устойчивость
Рассмотрим характеристический полином
Введем в рассмотрение -мерный вектор
. Пусть в
-мерном пространстве коэффициентов задано множество
. Полином
называется робастно устойчивым или робастно устойчивым в
, если он является устойчивым (т. е. все его нули являются левыми) при любых значениях коэффициентов
из множества
. Система называется робастно устойчивой или робастно устойчивой на множестве
, если ее характеристический полином является робастно устойчивым полиномом на множестве
.
Полиномы Харитонова. Пусть множество А является (гиперпараллелепипедом:
Здесь и
— минимальное и максимальное значения коэффициента
Полиномы
со следующими коэффициентами (коэффициенты выписаны в порядке убывания индексов)
называются полиномами Харитонова.
Необходимое условие робастной устойчивости. Так как при робастной устойчивости в параллелепипеде (3.4) должны быть устойчивыми характеристические полиномы при всех значениях коэффициентов из этого параллелепипеда, необходимо, чтобы был устойчивым характеристический полином при значениях коэффициентов
. Поэтому для робастной устойчивости на множестве (3.4) необходимо, чтобы при
выполнялось условие
Теорема Харитонова (1978). Для того чтобы система с характеристическим полиномом
была робастно устойчива на множестве (3.4), необходимо и достаточно, чтобы все полиномы Харитонова были устойчивыми.
В случае, когда = 1,2,3,4,5, нет необходимости проверять устойчивость всех четырех полиномов Харитонова. При
= 1,2 необходимое условие (3.6) является и достаточным. В случае выполнения необходимого условия робастной устойчивости для того, чтобы система была робастно устойчива на множестве (3.4), необходимо и достаточно:
а) при = 3 был устойчивым полином Харитонова
;
б) при = 4 были устойчивыми полиномы Харитонова
и
;
в) при = 5 — полиномы Харитонова
.
Задача №3.5.
Исследовать робастную устойчивость системы, характеристический полином которой имеет вид
Решение:
В данном случае
Так как = 4 и выполняется необходимое условие робастной устойчивости, достаточно рассмотреть полиномы Харитонова
и
.
Из (3.5а) и (3.56) имеем
или
Необходимое условие устойчивости для обоих полиномов выполняется. Для полинома определитель Гурвица
а для полинома
На основе критерия Льенара—Шипара и
являются устойчивыми полиномами. Следовательно, система робастно устойчива.
Задача №3.6.
Исследовать устойчивость замкнутой системы, если передаточная функция разомкнутой системы имеет вид
Решение:
Характеристический полином замкнутой системы имеет вид
где
Коэффициенты характеристического полинома удовлетворяют следующим условиям:
Следовательно, в принятых выше обозначениях имеем
Необходимое условие робастной устойчивости выполняется. Так как = 3, для робастной устойчивости необходимо и достаточно, чтобы полином
был устойчивым. Из (3.5а)
Определитель Гурвица
Поэтому замкнутая система не будет робастно устойчива. Теорема Харитонова справедлива при условии, что коэффициенты характеристического полинома изменяются на заданных интервалах независимо друг от друга. В противном случае устойчивость полиномов Харитонова является только достаточным условием робастной устойчивости.
Задача №3.7.
Исследовать устойчивость замкнутой системы при всевозможных заданных значениях параметров при условии, что передаточная функция разомкнутой системы имеет вид
Решение:
Характеристический полином замкнутой системы имеет вид
где
Для граничных значений коэффициентов характеристического полинома имеем
Необходимое условие робастной устойчивости выполняется. И так как = 3, достаточно рассмотреть полином
(3.5а):
Все коэффициенты больше нуля, но определитель Гурвица
Следовательно, полином не является устойчивым, т.е. условие робастной устойчивости не выполняется. Однако, в данном случае коэффициенты характеристического полинома не являются независимыми и теорема Харитонова определяет только достаточное условие робастной устойчивости. В действительности, как покажем, замкнутая система устойчива при всевозможных заданных значениях параметров.
При положительных значениях параметров необходимое условие устойчивости выполняется и определитель Гурвица
будет положительным при .
Следовательно, система устойчива при любых значениях параметров из области, определяемой неравенствами . Очевидно, заданные значения параметров принадлежат этой области.
Качество систем управления. Показатели качества в переходном режиме
Показатели качества делятся на показатели качества в переходном режиме и показатели качества в установившемся режиме.
Показатели качества в переходном режиме делятся на прямые и косвенные. Последние делятся на корневые, частотные и интегральные.
Прямыми показателями качества называются показатели, которые получаются непосредственно по переходной характеристике. Из прямых показателей качества наиболее часто используют время регулирования и перерегулирование.
Временем регулирования называют минимальное время, по истечении которого отклонение выходной величины от установившегося значения
не превышает некоторой заданной величины
(обычно принимают
, перерегулированием
— максимальное отклонение переходной функции от установившегося значения
, выраженное в процентах по отношению к
:
где — максимальное значение переходной функции.
Корневые показатели качества. В качестве корневых показателей используют степень устойчивости и колебательность (степень колебательности). Степенью устойчивости системы управления (или характеристического полинома) называют расстояние от мнимой оси до ближайшего корня ее характеристического уравнения на комплексной плоскости, или
степень колебательности системы (или ее характеристического полинома) можно определить следующим образом:
Здесь — корни характеристического уравнения.
При исследовании степени устойчивости удобно воспользоваться следующим преобразованием. Полином
преобразуем, сделав подстановку . Тогда получим:
Преобразование соответствует сдвигу мнимой оси влево и преобразованный полином
будет устойчивым полиномом, если
— степень устойчивости исходного полинома), и неустойчивым полином, если
. Поэтому исследование степени устойчивости полинома
сводится к исследованию устойчивости преобразованного полинома
.
Задача №4.1.
Задан характеристический полином
Исследовать, превышает ли степень устойчивости заданного полинома единицу.
Решение:
Убедимся сначала, что рассматриваемый полином является устойчивым полиномом, для чего вычислим определитель Гурвица 3-го порядка, составленный из его коэффициентов.
Полином является устойчивым. Сделаем подстановку
и вычислим коэффициенты преобразованного полинома. В данном случае
= 4 и
= 1. Поэтому из (4.1) имеем:
Без дальнейших вычислений ясно, что необходимое условие устойчивости преобразованного полинома не выполняется, и он является неустойчивым полиномом. Следовательно, степень устойчивости < 1.
Задача №4.2.
Определить, превышает ли единицу степень устойчивости характеристического полинома
Решение:
Сначала проверим устойчивость заданного полинома. Для этого достаточно проверить знак определителя Гурвица 2-го порядка:
Полином устойчив. Произведем подставку
и найдем коэффициенты преобразованного полинома. В данном случае
= 3 и
= 1. Поэтому из (4.1) имеем
Преобразованное характеристическое уравнение имеет вид
Все корни этого уравнения располагаются на мнимой оси. Следовательно, степень устойчивости рассматриваемой системы
= 1.
Интегральные показатели качества. В качестве интегральных оценок наиболее часто используют интегральную квадратическую ошибку (оценку)
и обобщенные интегральные квадратические оценки
Здесь — переходная составляющая ошибки:
,
— установившаяся ошибка;
— весовые константы.
Вычисление интегральных квадратических оценок. Из равенства Парсеваля
где
имеем
где
Так как
формулу для можно записать в виде
Определение интегральных квадратических показателей сводится к вычислению интеграла вида
Этот интеграл вычисляется с помощью теории вычетов и для = 1,2,3 имеет следующий вид
Задача №4.3.
Вычислить интегральные показатели и
системы (рис. 4.1, а), когда передаточная функция
Решение:
Вычислим и
, необходимые для нахождения указанных показателей. Но, прежде всего, найдем
. Учитывая, что
можно записать
Установившееся значение
Так как
то
На основании свойства преобразования Лапласа
Интегральная квадратическая оценка имеет вид
В данном случае (4.2)
Поэтому согласно (4.3а)
Теперь найдем :
Так как
имеем
Показатели качества в установившемся режиме
Наиболее полной характеристикой качества системы в установившемся режиме является установившаяся ошибка. Если на систему действуют два внешних воздействия — задающее воздействие и возмущение
— установившуюся ошибку можно представить в виде суммы:
где и
— установившиеся ошибки от задающего воздействия
и возмущения
соответственно.
Установившиеся ошибки и
можно представить в виде ряда
где
Здесь — передаточная функция относительно входа
и выхода
— передаточная функция относительно входа
и выхода
. Предполагается, что возмущение не приложено в одной точке с задающим устройством. Коэффициенты
называются коэффициентами ошибки по задающему воздействию, коэффициенты
— коэффициентами ошибки по возмущению.Коэффициенты
и
называют коэффициентами позиционной ошибки,
и
— коэффициентами скоростной ошибки,
и
— коэффициентами ошибки по ускорению.
Статические и астатические системы. Установившаяся ошибка при постоянном внешнем воздействии называется статической ошибкой. Система называется статической, если статическая ошибка отлична от нуля, и астатической, если статическая ошибка равна нулю.
Система называется статической относительно задающего воздействия (возмущения), если статическая ошибка от задающего воздействия (возмущения) отлична от нуля, и астатической относительно задающего воздействия (возмущения), если статическая ошибка от задающего воздействия (возмущения) равна нулю.
Формулы (4.4) и (4.5) при постоянных и
принимают вид
Отсюда следует, что система будет статической относительно воздействия (возмущения
), если
, и астатической относительно задающего воздействия
(возмущения
), если
.
Говорят, что астатическая система обладает астатизмом -го порядка относительно задающего воздействия, если
Аналогично определяется астатическая система с астатизмом -го порядка относительно возмущения.
Если система обладает астатизмом -го порядка, то коэффициенты ошибок
при
можно определить следующим образом:
Иначе говоря, этими формулами можно пользоваться при вычислении до первого отличного от нуля коэффициента.
Задача №4.4.
Определить установившуюся ошибку системы (рис. 4.1, б) при
Решение:
Так как все производные от и производные выше 1-го порядка от
равны нулю, то в данном случае
Поэтому для определения искомой ошибки достаточно вычислить коэффициенты ошибок
Передаточные функции ошибки имеют вид
Отсюда
Так как
можно вычислить по формуле (4.6).
Таким образом, для ошибок имеем:
Структура астатической системы управления. Для того чтобы система управления была астатической с астатизмом -го порядка относительно задающего воздействия, нужно, чтобы она содержала
последовательно соединенных интегрирующих звеньев во всем замкнутом контуре.
Для того чтобы система управления была астатической с астатизмом -го порядка относительно возмущения, нужно, чтобы она содержала
последовательно соединенных интегрирующих звеньев, включенных между точкой съема ошибки
и точкой приложения возмущения .
Синтез систем управления
При выборе законов управления следует иметь в виду:
- введение в закон управления интегрирующего члена делает систему астатической и улучшает качество системы в установившемся режиме, но оказывает дестабилизирующее влияние (т. е. может сделать систему неустойчивой) и ухудшает качество системы в переходном режиме;
- введение в закон управления дифференцирующего члена оказывает стабилизирующее влияние (может сделать неустойчивую систему устойчивой) и улучшает качество системы в переходном режиме, не оказывая влияние на качество системы в установившемся режиме.
Задача №5.1.
Определить, при каких типовых законах управления статическая ошибка системы (рис. 5.1) будет равна нулю, когда передаточная функция объекта имеет вид
Решение:
Статическая ошибка будет равна нулю, если система будет астатической относительно задающего воздействия и возмущения. А для этого нужно, чтобы регулятор содержал интегрирующее звено. Поэтому искомыми законами управления будут пропорциональ-но-интегральный (ПИ) закон и пропорционально интегро-дифференци-альный (ПИД) закон.
Задача №5.2.
Определить, при каких типовых законах управления установившаяся ошибка системы (рис. 5.1) будет равна нулю при условии, что
Решение:
Так как установившаяся ошибка от задающего воздействия и возмущения имеют вид
установившаяся ошибка будет равна нулю, если 0. Следовательно, система должна быть астатической с астатизмом 2-го порядка относительно задающего воздействия и с астатизмом 1-го порядка относительно возмущения. Так как объект включает два последовательно соединенных интегрирующих звена, система будет астатической с астатизмом не менее 2-го порядка относительно задающего воздействия при любом типовом законе управления. Однако она будет астатической относительно возмущения только при ПИ-законе и ПИД-законе. При ПИ-законе передаточная функция разомкнутой системы
и характеристическое уравнение имеет вид
В этом уравнении коэффициент при равен нулю и необходимое условие устойчивости не выполняется. Поэтому система при ПИ-законе структурно неустойчива. При ПИД-законе передаточная функция разомкнутой системы
и характеристическое уравнение имеет вид
Определитель Гурвица 3-го порядка
соответствующим выбором параметров регулятора можно сделать Следовательно, при ПИД-законе система структурно устойчива и искомым законом устойчива и искомым законом управления является ПИД-закон.
Синтез параметров регулятора по минимуму интегральных оценок
Постановку и решение задачи синтеза параметров регулятора по минимуму интегральной оценки рассмотрим на примерах.
Задача №5.3.
При условии, что
и
определить параметр , при котором переходный процесс системы (рис. 5.1) является апериодическим и интегральная квадратическая ошибка
принимает минимальное значение.
Решение:
Переходный процесс будет апериодическим, если корни характеристического уравнения рассматриваемой системы
будут вещественными, т. е. если детерминант этого уравнения
Так как , то ошибка
. Объект включает интегрирующее звено. Поэтому система является астатической задающего воздействия и статическая ошибка
. Переменная составляющая ошибки
Переходя к изображениям Лапласа, получим:
Следовательно,
В данном случае (4.2)
Поэтому (4.36)
Очевидно, что принимает минимальное значение при условии
2,5; когда
.
Задача №5.4.
При условии, что
и
(рис. 5.1), определить значение параметра при котором обобщенная интегральная квадратическая оценка
при
принимает минимальное значение.
Решение:
Согласно формуле Парсеваля
Как было вычислено (см. пример 5.3),
Для имеем:
Поэтому
В данном случае
Подставив это выражение и выражение для в полученную выше формулу для
, найдем
Из условия
следует, что достигает экстремума при
установить, чему (минимуму или максимуму) соответствует это значение, найдем вторую производную
В точке экстремума эта производная положительна. Следовательно, в ней достигается минимум и соответственно решением будет = 2.
Синтез систем управления максимальной степени устойчивости
Задача синтеза систем управления максимальной степени устойчивости ставится следующим образом. Задана структура системы управления и требуется определить (
— вектор параметров регулятора) из условия
Здесь называется оптимальной степенью устойчивости и
— оптимальным (векторным) параметром. Число параметров регулятора (размерность вектора
)
не должно превышать
— 1 (
— степень характеристического уравнения). Метод решения сформулированной задачи основан на условиях граничной устойчивости.
Условия граничной (маргинальной) устойчивости. Система находится на границе устойчивости или имеет место граничная (маргинальная) устойчивость, если ее характеристический полином имеет нейтральные (т. е. расположенные на мнимой оси) нули и не имеет правых нулей. Такой полином называют маргинально устойчивым.
Рассмотрим полином с вещественными коэффициентами
Утверждение. 5.1 (необходимое условие маргинальной устойчивости). Если полином (5.1) маргинально устойчив, то все его коэффициенты неотрицательны:
Нуль полинома (5.1) называют особым, если —
также является нулем этого полинома. В частности, все нули, расположенные на мнимой оси, являются особыми.
Утверждение. 5.2 Полином (5.1) маргинально устойчив и нулей располагаются на мнимой оси в том и только в том случае, если выполняются следующие два условия:
1) старших определителей Гурвица равны нулю, а остальные
определителей положительны:
2) Полином (5.1) не имеет особых нулей, расположенных не на мнимой оси.
Утверждение. 5.3 При выполнении необходимого условия (5.2) особый нуль не может быть вещественным числом, и если имеются особые нули, расположенные не на мнимой оси, то их количество равно числу, кратному четырем.
Нейтральные нули полинома имеют вид
и их число совпадает с числом действительных корней уравнения
, или системы уравнений
Утверждение. 5.4 Для того чтобы все определители Гурвица полинома были равны нулю, необходимо и достаточно, чтобы все его коэффициенты с нечетными индексами были равны нулю.
Метод синтеза систем управления максимальной степени устойчивости. Метод решения задачи основан на преобразовании характеристического полинома
путем постановки . При этой постановке преобразованный полином
становится маргинально устойчивым полиномом. И для выписываются условия маргинальной устойчивости, включающие условия (5.2), (5.3) и (5.4):
Здесь — определители Гурвица преобразованного полинома
— Следует иметь в виду, что не все соотношения в (5.76) и (5.7в) являются независимыми.
Рассматриваемый метод состоит в следующем: решается система (5.7) относительно неизвестных параметров регулятора и степени устойчивости и находятся решения, у которых
имеет наибольшее значение.
Утверждение. 5.5 Максимально возможная или граничная степень устойчивости устойчивого полинома
равна
и она достигается, когда вещественные части всех нулей полинома равны между собой.
Поиск решения задачи синтеза максимальной степени устойчивости следует начинать со случая, когда степень устойчивости принимает граничное (максимально возможное) значение. Так как это возможно, когда все нули исходного полинома имеют одинаковые вещественные части или все нули преобразованного полинома располагаются на мнимой оси, условие маргинальной устойчивости (5.7) можно представить в виде
Если эта система не имеет решения, то нужно перейти к системе (5.7) и решить ее при .
А) Синтез оптимальных по степени устойчивости параметров типовых регуляторов для объекта 2-го порядка
Рассмотрим синтез оптимальных по степени устойчивости параметров П- и ПИ-регуляторов для объекта 2-го порядка. Пусть передаточная функция объекта имеет вид
П-регулятор. Передаточная функция регулятора и передаточная функция разомкнутой системы равна
Характеристический полином принимает вид
Для коэффициентов преобразованного полинома
в соответствии с (5.6) имеем
В данном случае условия граничной устойчивости (5.9) принимают следующий вид:
Решив эту систему, получим
Так как степень устойчивости принимает граничное значение, найденное решение является искомым. Здесь — свободный параметр, пропорциональный степени колебательности.
ПИ-регулятор. Передаточная функция регулятора и передаточная функция разомкнутой системы равна
Характеристический полином имеет вид
Коэффициенты преобразованного полинома
определяются следующим образом (5.6):
Условия граничной устойчивости (5.9) принимают вид
Из последнего уравнения (5.10) имеем . Следовательно, неравенство в этом условии выполняется. Поэтому исключив его и подставив выражения для коэффициентов, условие (5.10) можно представить в виде
Решив эту систему, получим
Здесь — свободный параметр, представляющий собой мнимую часть комплексных корней характеристического уравнения синтезированной системы.
Б) Синтез оптимальных по степени устойчивости параметров ПД- и ПИД-регуляторов для объекта 3-го порядка
Рассмотрим синтез оптимальных по степени устойчивости параметров типовых регуляторов для объекта 3-го порядка. Пусть передаточная функция объекта имеет вид
ПД-регулятор. Передаточная функция регулятора и передаточная функция разомкнутой системы
Характеристический полином замкнутой системы имеет вид
Для коэффициентов преобразованного полинома
имеем
Условия граничной устойчивости (5.9) для преобразованного полинома принимают вид
Из последнего равенства этого условия имеем . Подставив в это и другие равенства условия маргинальной устойчивости выражения для
, получим
Решив эту систему, найдем
Здесь — свободный параметр, представляющий собой мнимую часть комплексных корней характеристического уравнения синтезированной системы.
ПИД-регулятор. Передаточная функция регулятора и передаточная функция разомкнутой системы
Характеристический полином синтезируемой системы и преобразованный полином имеют соответственно вид
где
Условие маргинальной устойчивости (5.9) принимает вид
Неравенства и
будут выполнены, если
. Введя дополнительный параметр
, последнее неравенство преобразуем в равенство
и условие маргинальной устойчивости можно записать в виде
Решив эту систему уравнений, получим
где свободные параметры и
являются мнимыми частями корней характеристического уравнения синтезированной системы
Синтез систем управления по желаемой передаточной функции или метод полиномиальных уравнений
При задании желаемой передаточной функции и определении передаточной функции регулятора
необходимо учитывать физическую осуществимость определяемого регулятора и грубость синтезируемой системы.
Физическая осуществимость. Под физической осуществимостью или реализуемостью передаточной функции или системы, заданной этой передаточной функцией, понимают принципиальную возможность построения такой системы.
Передаточная функция физически осуществима, если степень числителя не больше степени ее знаменателя. Условие физической осуществимости передаточной функции
имеет вид
Грубость. Система называется грубой или робастной, если при малом изменении ее параметров свойство системы качественно не меняется. В случае линейной системы негрубость означает, что устойчивая система при малом изменении параметров становится неустойчивой.
При синтезе систем по желаемой передаточной функции грубость может быть нарушена, если правый полюс передаточной функции объекта компенсируется правым нулем передаточной функции регулятора и правый нуль объекта — правым полюсом регулятора.
Представим передаточную функцию объекта в виде
где — полиномы с левыми нулями,
— полиномы с правыми и нейтральными нулями. Если полиномы
и
не содержат левых нулей, то
равны константе; если они не содержат правых нулей, то
следует принять равными единице.
Передаточная функция регулятора синтезируемой системы имеет вид
где полиномы и
определяются из полиномиального уравнения
Здесь — характеристический полином синтезируемой системы.
Условимся степень полинома обозначать буквой с индексом, обозначающим сам полином. Например,
будет обозначать степень полинома
. Коэффициенты полиномов
и
определяются из системы уравнений, которые получаются путем приравнивания коэффициентов при одинаковых степенях обеих частей полиномиального уравнения (5.12). Эта система разрешима при выполнении условия разрешимости
Условие физической осуществимости регулятора (5.11):
При определении степеней неопределенных полиномов необходимо учитывать условие, получаемое из условия грубости:
Чтобы система (5.13)-(5.15) была разрешима, необходимо, чтобы порядок характеристического полинома синтезируемой системы
удовлетворял соотношению
Метод синтеза регулятора по желаемой передаточной функции состоит в следующем. Исходя из заданных требований к качеству синтезируемой системы задается характеристический полином с учетом условия (5.16). Из (5.13)—(5.15) определяют степени неопределенных полиномов
и
. Чтобы не усложнять регулятор, находят наименьшие возможные значения. Затем составляют полиномы
и
с неопределенными коэффициентами, подставляют их в полиномиальное уравнение и определяют неизвестные коэффициенты. Найденные полиномы
и
подставляют в (5.11) и получают искомую передаточную функцию регулятора.
Задача №5.5.
Передаточная функция объекта имеет вид . Определить передаточную функцию регулятора, при которой переходная составляющая ошибки
изменяется в соответствии с функцией
и установившаяся ошибка равна нулю (рис. 5.1) при постоянном задающем воздействии и отсутствии возмущения
.
Решение:
Переходная составляющая ошибки будет изменяться в соответствии с заданной функцией, если характеристический полином синтезируемой системы имеет трехкратный корень, равный — 1 . Поэтому знаменатель желаемой передаточной функции
имеет вид
Числитель и знаменатель передаточной функции объекта раскладываются на множители
Степени полиномов равны
Статическая ошибка , если система обладает астатизмом 1-го порядка
относительно задающего воздействия. Так как объект содержит одно интегрирующее звено, можно принять
. Условия (5.13)—(5.15) принимают вид
Из последнего равенства = 2. Наименьшим
удовлетворяющим приведенным условиям, является
= 0. Поэтому
и
Подставив эти полиномы в полиномиальное уравнение (5.12), получим
Отсюда, приравнивая коэффициенты при одинаковых степенях, найдем
Подставляя эти полиномы, а также выражения для и
в (5.11), получим искомую передаточную функцию регулятора
.
Определение желаемой передаточной функции
Желаемая передаточная функция должна быть определена исходя из заданных требований к качеству синтезируемой системы. На ее выбор определенные ограничения накладывают условия грубости и физической осуществимости. В силу этих ограничений желаемая передаточная функция имеет вид
где — множитель полинома числителя передаточной функции объекта с правыми нулями;
— полином, определяемый в процессе синтеза. Поэтому определение желаемой передаточной функции практически сводится к выбору полинома
.
Передаточная функция вида
называется нормированной передаточной функцией или передаточной функцией в форме Вышнеградского. Нормированная передаточная характеризуется тем, что в знаменателе коэффициент при старшей степени и свободный член равны единице.
Желаемая передаточная функция, когда наряду с другими требованиями нужно обеспечить заданное время регулирования , определяется следующим образом. По заданным требованиям к качеству синтезируемой системы, кроме требования к времени регулирования, находится стандартная нормированная передаточная функция вида (5.17) и для нее определяется время регулирования
. По полученному
и заданному tp находится отношение
. Коэффициенты желаемой передаточной функции
определяются так:
Стандартные нормированные передаточные функции
Рассмотрим стандартные передаточные функции, которые не имеют нулей: числители являются константами. Так как значения этих констант не влияют на характер переходного процесса, примем их равными единице.
Передаточная функция с одинаковыми полюсами
обладает монотонной переходной характеристикой, неплохим быстродействием и среди передаточных функций —го порядка с одинаковыми коэффициентами при
имеет наибольшую степень устойчивости. Ее знаменатель при
= 4, 5, 6 принимает следующий вид:
Оптимальная по быстродействию передаточная функция передаточная функция с полюсами, имеющими одинаковые действительные части и мнимые части, образующие арифметические прогрессии с разностью и первым членом, равными
, и отношением
которому соответствует наименьшее время регулирования:
В таблице 5.1 представлены время регулирования при и перерегулирование для приведенных стандартных нормированных передаточных функций.
Задача №5.6.
Передаточная функция объекта имеет вид
Синтезировать регулятор, при котором переходный процесс является монотонным, время регулирования и статическая ошибка равна нулю.
Решение:
Статическая ошибка будет равна нулю, если система будет астатической. Примем порядок астатизма .
В качестве стандартной передаточной функции выберем нормированную передаточную функцию с одинаковыми полюсами
Для этой передаточной функции из табл. 5.1, а имеем и поэтому
Учитывая формулы (5.186), для коэффициентов полинома знаменателя желаемой передаточной функции получаем
и соответственно,
Числитель и знаменатель передаточной функции объекта раскладываются на множители
Степени полиномов равны
Условия (5.13), (5.14) и (5.15) принимают вид
Этим условиям удовлетворяют и, соответственно,
При подстановке этих полиномов уравнение (5.12) принимает вид
Отсюда
и, соответственно
Подставляя их и выражения для и
в (5.11), найдем искомую передаточную функцию
Метод обратной задачи динамики
Методом обратной задачи динамики называют метод синтеза систем, когда по заданным уравнению объекта и требованиям к качеству системы управления определяется дифференциальное уравнение, решение которого удовлетворяет заданным требованиям, а затем из найденного уравнения выражается старшая производная и, после подстановки ее вместо старшей производной в уравнение объекта, находится требуемый закон управления.
Задача №5.7.
Пусть задана передаточная функция объекта
Задан требуемый закон изменения выходной переменной
. Требуется найти алгоритм управления, при котором ошибка
изменяется следующим образом:
Здесь — произвольные постоянные;
— заданные положительные постоянные.
Решение:
Уравнение объекта имеет вид
Числа являются корнями уравнения
Следовательно, заданная функция является общим решением дифференциального уравнения
Так как
подставив эти выражения в уравнение объекта, получим
Математическое описание дискретных систем
Определение -преобразования,
-преобразованием или преобразованием Лорана называется соотношение
ставящее дискретной функции в соответствие функцию комплексного переменного
. При этом
называют оригиналом, a
— изображением или
-изображением. Сумма в правой части называется рядом Лорана.
Оригинал и его изображение обозначают одноименными буквами: оригинал — строчной буквой, а изображение — прописной буквой со звездочкой. -преобразование также условно записывают в виде
а обратное -преобразование — в виде
-преобразование от смещенной решетчатой функции
называют модифицированным -преобразованием. Модифицированное
-преобразование также записывают в виде
Функцию называют
-изображением смещенной решетчатой функции
или модифицированным
-изображением решетчатой функции
.
Уравнения и передаточные функции дискретных систем
Пусть модель дискретной системы управления описывается разностным уравнением
где — выходная переменная,
— входная переменная,
и
— константы. Используя оператор смещения
это уравнение можно записать в операторной форме
Разностный оператор при выходной переменной
называется собственным (разностным) оператором, а разностный оператор при входной переменной
(разностным) оператором воздействия.
Отношение оператора воздействия к собственному оператору называется передаточной функцией в операторной форме. В соответствии с этим определением передаточная функция (в операторной форме) системы управления (6.1) имеет вид
Имеющее наименьший порядок отношение -изображений выходной и входной переменных, вычисленных при нулевых начальных условиях, называется передаточной функцией в
-изображениях.
Передаточные функции в -преобразованиях
и в операторной форме
связаны соотношением
Однако если полиномы числителя и знаменателя имеют общие нули, то они должны быть сокращены.
Вычисление передаточных функций АИМ-системы
Как правило, приходится вычислять передаточные функции, когда известны характеристики дискретных элементов и передаточная функция непрерывной части. И в этом случае возникают особенности, которые делают вычисление передаточных функций дискретных систем более сложным.
АИМ-система включает АИМ-элемент (импульсный элемент с амплитудно-импульсной модуляцией) и непрерывную часть (рис. 6.2). Для получения математического описания АИМ-системы управления
ее представляют в виде эквивалентной схемы, состоящей из простейшего импульсного звена 1 и приведенной непрерывной части (ПНЧ) (рис. 6.3, а). Простейшее импульсное звено представляет собой звено, которое преобразует входную функцию в обобщенную решетчатую функцию
где — период выходного сигнала АИМ-элемента. Передаточная функция ПНЧ равна произведению передаточных функций непрерывной части и формирующего звена. Формирующее звено (ФЗ) формирует из обобщенной решетчатой функции
сигнал, тождественно равный выходному сигналу АИМ-элемента, и его передаточная функция равна изображению функции, описывающей немодулированный импульс.
Если ограничиться изучением АИМ-системы только в дискретные моменты времени , то получим дискретную модель (рис. 6.3, б), в которой передаточная функция
равна
-изображению решетчатой весовой функции ПНЧ:
Зная связь между изображением Лапласа непрерывной функции и -изображением соответствующей решетчатой функции (см. табл. 6.1), можно непосредственно по передаточной функции ПНЧ
определить
. Для этого введем в рассмотрение оператор
, который каждой функции
ставит в соответствие функцию
:
Оператор соответствует трем последовательным операциям: обратному преобразованию Лапласа, квантованию по времени и
-преобразованию. Так как все три указанные операции являются линейными, то оператор
является линейным. Используя этот оператор, передаточную функцию
можно определить следующим образом:
Дальше также используется оператор который функции
ставит в соответствие модифицированное
-изображение
По аналогии с -преобразованием в
-преобразовании
и в -преобразовании
называют оригиналом, а
-изображением и
—
-изображением или модифицированным
-изображением.
и
-изображения от основных функций можно найти в табл. 6.1 и 6.2 соответственно.
Вычисление и
-изображений. Пусть оригинал имеет вид
где и
— полиномы от
степени
и
соответственно, причем
. Если все полюса
данной функции (т. е. корни уравнения
= 0) различны, то
Задача №6.1.
Передаточная функция ПНЧ имеет вид
Требуется найти дискретную передаточную функцию
Решение:
Полюсами данной передаточной функции (т. е. корнями уравнения ) являются
. Производная
. Поэтому по формуле (6.2)
Если
содержит кратные полюса, то изображения и
можно получить, разложив
на элементарные дроби. В простых случаях можно введением малых параметров видоизменить функцию
так, чтобы она не содержала кратных полюсов, и воспользоваться формулами (6.2) и (6.3), а затем произвести предельный переход, устремив малые параметры к нулю.
Задача №6.2.
Передаточная функция ПНЧ имеет вид Требуется определить дискретную передаточную функцию
.
Решение:
Данная передаточная функция ПНЧ имеет двукратный полюс . Введя малый параметр
, преобразуем ее к виду
Преобразованная передаточная функция имеет простые полюса и
. Производная
. По формуле (6.2)
Используя разложение
где — бесконечно малая величина более высокого порядка, чем
, получаем
Отсюда, устремив к нулю, находим
Если среди простых полюсов функции имеются комплексные корни, то может оказаться нецелесообразным использование формул (6.2) и (6.3). Это связано с необходимостью преобразования полученного результата для исключения мнимого числа. Во всех случаях, когда использование формул (6.2) и (6.3) невозможно или нецелесообразно, можно определить
и
, разложив
на элементарные дроби.
Задача №6.3.
Определить и
изображение функции
Решение:
Данная функция имеет кратный полюс и два комплексных полюса. Найдем
и
, разложив
на элементарные дроби методом неопределенных коэффициентов:
Приравняв коэффициенты при одинаковых степенях слева и справа и решив полученную систему уравнений, найдем
Следовательно,
Преобразуем правую часть к табличному виду:
Подставив это выражение в предыдущее равенство, произведем -преобразование:
После подстановки соответствующих изображений из табл. 6.2 и преобразований получим
Положив , находим
Вычисление — и
-изображений от оригинала, включающего множитель
. Пусть оригинал имеет вид
где — дробно-рациональная функция:
В этом случае в зависимости от величины для
имеем
а) при
где
б) при
-преобразование обладает следующим свойством: если оригинал в
-преобразовании содержит множитель, представляющий полином или дробно-рациональную функцию от
, то этот множитель можно вынести за знак оператора
, произведя подстановку
. Например,
Задача №6.4.
АИМ-элемент вырабатывает прямоугольные импульсы длительности с периодом
и амплитудой (высотой)
. Передаточная функция непрерывной части
. Требуется определить дискретную передаточную функцию
.
Решение:
Найдем сначала передаточную функцию приведенной непрерывной части. Так как передаточная функция формирующего звена
передаточная функция ПНЧ
где
Дискретная передаточная функция
В данном случае
Согласно (6.4)
Полюсами являются
, производная
. В соответствии с (6.2) и (6.3)
Следовательно,
Цифровые системы управления
Если цифровое устройство оперирует числовыми представлениями со значительным количеством разрядов, то- квантованием по уровню можно пренебречь. И системы управления с такими цифровыми устройствами можно рассматривать как АИМ-системы.
Цифровая система управления (ЦСУ) включает объект управления (ОУ), чувствительные элементы (ЧЭ), аналого-цифровой преобразователь (АЦП), цифровое вычислительное устройство (ЦВУ) и цифро-аналоговый преобразователь (ЦАП) (рис. 6.4). АЦП преобразует аналоговый сигнал в цифрой, а ЦАП — цифровой сигнал в аналоговый. ЦВУ выполняет все необходимые вычисления в соответствии с заданным алгоритмом управления, т. е. представляет собой регулятор.
Если пренебречь квантованием по уровню, цифровую систему управления можно представить в виде блок-схемы (рис. 6.5), состоящей из прерывателя, дискретного фильтра (ДФ), фиксатора нулевого порядка (ФНП) и непрерывной части (НЧ).
Прерыватель является моделью АЦП и преобразует непрерывный сигнал в дискретный сигнал
. В дальнейшем прерыватель в явном виде на схеме не будем указывать, принимая, что он входит в состав ДФ. Дискретный фильтр представляет собой модель ЦВУ и характеризуется дискретной передаточной функцией — передаточной функцией регулятора. В качестве ЦАП чаще всего используется фиксатор нулевого порядка — элемент, который запоминает входной дискретный сигнал на один период — до прихода следующего дискретного сигнала. Фиксатор нулевого порядка можно рассматривать как АИМ-элемент, вырабатывающий прямоугольные импульсы длительности
(относительная длительность
) и с амплитудой
. Представив ФНП в виде эквивалентной схемы, состоящей из простейшего импульсного элемента и формирующего звена, получим эквивалентную схему цифровой системы управления (рис. 6.6).
На этой схеме — передаточная функция (в операторной форме) дискретного фильтра (регулятора),
— передаточная функция ПНЧ. Передаточная функция (в изображениях Лапласа) формирующего звена
Передаточная функция (в изображениях Лапласа) ПНЧ
Дискретная передаточная функция ПНЧ
или
Используя эту передаточную функцию, можно построить структурную схему дискретной модели цифровой системы управления (рис. 6.7).
Задача №6.5.
Дана цифровая система управления, у которой передаточная функция непрерывной части
и цифровое вычислительное устройство реализует алгоритм управления, определяемый разностным уравнением
Требуется определить передаточную функцию данной системы относительно входа и выхода
(рис. 6.6).
Решение:
Запишем уравнение регулятора в операторной форме:
Отсюда передаточная функция регулятора в операторной форме
и в -изображениях
Передаточная функция приведенной непрерывной части
Дискретная передаточная функция ПНЧ
Корнями полинома
являются
и производная
По формуле (6.2)
Следовательно,
Искомая передаточная функция замкнутой системы
ШИМ-системы управления
Блок-схема ШИМ-системы управления включает ШИМ-элемент (импульсный элемент с широтно-импульсной модуляцией) и НЧ (рис. 6.8). Пусть ШИМ-элемент вырабатывает прямоугольные импульсы с амплитудой и периодом
. На выходе ШИМ-элемента ширина модулированного импульса пропорциональна модулю
, а ее знак совпадает со знаком входного сигнала в момент съема. Модулированный импульс на выходе ШИМ-элемента можно представить как разность двух ступенчатых функций:
где Здесь
является константой, удовлетворяющей неравенству
и называется коэффициентом модуляции.
Линеаризация. Уравнения ШИМ-элемента являются нелинейными. Если выполняется условие или
, то можно произвести линеаризацию и получить дискретно-непрерывную модель (рис. 6.9 а), а после дискретизации — дискретную модель (рис. 6.9 б). Здесь
Задача №6.6.
Дана ШИМ-система управления (рис. 6.8). Амплитуда = 1, коэффициент модуляции
= 0,05, период следования импульсов
= 0,1 и передаточная функция непрерывной части
Требуется определить дискретную передаточную функцию замкнутой системы
Решение:
Так как
Искомая передаточная функция
Вычисление передаточных функций дискретных систем в общем случае
Выше мы рассмотрели вычисление передаточных функций дискретных систем, когда их эквивалентная схема за простейшим импульсным звеном содержит одно непрерывное звено — приведенную НЧ. Однако может потребоваться вычисление передаточных функций, эквивалентная схема которых имеет более общий вид (рис. 6.10). И в этом
случае справедливо правило, которое совпадает с правилом вычисления передаточных функций одноконтурной непрерывной системы: передаточная функция относительно входа и какого-либо выхода равна передаточной функции прямой цепи, деленной на единицу плюс (а при положительной обратной связи минус) передаточная функция разомкнутой системы. Согласно этому правилу имеем
Следует иметь в виду, что при вычислении передаточной функции прямой цепи и передаточной функции разомкнутой системы непрерывные звенья, расположенные за простейшим импульсным звеном, нужно рассматривать как одну НЧ.
Теперь рассмотрим схему с дискретным фильтром, включенным перед простейшим импульсным звеном (рис. 6.11). Установленное выше правило вычисления дискретной передаточной функции замкнутой системы остается в силе и в данном случае:
Задача №6.7.
Пусть в дискретной системе, представленной на рис. 6.11,
и период следования импульсов = 0,1. Требуется определить передаточные функции
и
.
Решение:
Найдем необходимые для определения требуемых передаточных функций -изображения. Учитывая, что полином
, как частный случай дробно-рациональной функции от
, можно вынести за знак оператора
, сделав подстановку
, получим
Подставив полученные выражения и выражения для в выше-приведенные формулы, получим
Устойчивость дискретных систем. Характеристическое уравнение и основное условие устойчивости
Если внешние воздействия заданы, уравнения дискретной системы управления можно записать в виде
или в операторной форме
Характеристическое уравнение имеет вид
который получается при подстановке в собственный оператор
вместо оператора смещения переменной
.
Если задана передаточная функция системы управления, то при определении характеристического полинома нужно исходить из следующих положений: по определению передаточной функции в операторной форме ее знаменатель есть собственный оператор, а знаменатель передаточной функции в -изображениях совпадает с характеристическим полиномом (при условии, что передаточная функция в операторной форме не содержит одинаковые нули и полюса).
Общее решение неоднородного разностного уравнения имеет вид
где — частное решение этого уравнения и
— общее решение соответствующего однородного уравнения.
Линейная дискретная система управления называется устойчивой, если общее решение однородного разностного уравнения при стремится к нулю:
.
Основное условие устойчивости: для того чтобы линейная дискретная система управления была устойчива, необходимо и достаточно, чтобы все корни ее характеристического уравнения были по модулю меньше единицы или, что то же, находились внутри единичного круга на -плоскости корней.
Задача №7.1.
Передаточная функция системы
Требуется исследовать ее устойчивость.
Решение:
Характеристическое уравнение имеет вид
Его корнями являются
Их модули
Система устойчива.
Алгебраические критерии устойчивости
Необходимое условие устойчивости: для того чтобы все нули (корни) характеристического полинома
были по модулю меньше единицы необходимо, чтобы при
выполнялись неравенства
Задача №7.2.
Характеристический полином дискретной системы имеет вид
Требуется определить устойчивость системы.
Решение:
Проверим необходимое условие устойчивости. В данном случае и
Необходимое условие устойчивости не выполняется. Следовательно, система неустойчива.
Исследование устойчивости, основанное на преобразовании единичного круга в левую полуплоскость. При преобразовании
внутренность единичного круга на -плоскости преобразуется в левую полуплоскость, его внешность — в правую полуплоскость и окружность (единичного радиуса) — в мнимую ось на ^-плоскости. При таком преобразовании переменной характеристического уравнения для исследования устойчивости дискретных систем можно воспользоваться критериями устойчивости непрерывных систем (критерий Гурвица и др.).
Представим преобразованное характеристическое уравнение в стандартной форме:
При = 1,2,3 коэффициенты преобразованного уравнения выражаются через коэффициенты исходного уравнения следующим образом:
Для того чтобы дискретная система была устойчива, необходимо и достаточно, чтобы все корни преобразованного характеристического уравнения располагались в левой полуплоскости (имели отрицательную вещественную часть).
Задача № 7.3.
Характеристический полином дискретной системы управления имеет вид
Определить ее устойчивость.
Решение:
В данном случае
и в соответствии с (7.2в) коэффициенты преобразованного уравнения
Необходимое условие устойчивости выполняется: все коэффициенты преобразованного характеристического уравнения больше нуля. Определитель Гурвица 2-го порядка
Следовательно, система устойчива.
Критерий устойчивости Джури. Составим таблицу Джури, которая содержит (+ 1) строку и столько же столбцов. При этом заполненные клетки имеют треугольную форму: нулевая строка содержит (
+ 1) заполненных клеток, а все последующие строки имеют на единицу меньше заполненных клеток, чем предыдущая строка (табл. 7.1).
Клетки нулевой строки заполняются коэффициентами характеристического уравнения в порядке возрастания нижних индексов:
Элементы первой строки
вычисляются следующим образом. Выписываются элементы нулевой строки и под ними те же элементы в обратном порядке. Из элементов верхней строки вычитаются соответствующие элементы нижней строки, умноженные на отношение последних элементов двух выписанных строк:
Последняя разность обращается в нуль, и она отбрасывается. Поэтому 1-я строка содержит элементов — на один элемент меньше, чем нулевая строка. Элементы всех последующих строк определяются аналогично элементам 1-й строки. Так, например, для вычисления
-й строки выписываются элементы (
— 1)-й строки и под ними те же элементы в обратном порядке. Из элементов верхней строки вычитаются соответствующие элементы нижней строки, умноженные на отношение последних элементов выписанных двух строк
Последняя разность, обращающаяся в нуль, отбрасывается. Формула для вычисления -го элемента
-й строки
имеет вид
Критерий Джури (E.I. Jury). Для того чтобы все нули (корни) характеристического полинома
находились внутри единичного круга, необходимо и достаточно, чтобы при все элементы нулевого столбца таблицы Джури были положительны:
Если все элементы нулевого столбца, кроме последнего, положительны:
то положительность последнего элемента, т. е. условие > 0 эквивалентно необходимому условию устойчивости (7.1). Поэтому если необходимое условие выполняется, то последний элемент
можно не вычислять.
Задача №7.4.
Характеристический полином дискретной системы управления имеет вид
Исследовать устойчивость данной системы.
Решение:
Сначала проверим необходимое условие устойчивости:
Необходимое условие устойчивости выполняется. Вычислим элементы таблицы Джури. Для нулевой строки имеем
Ниже приводится вычисление элементов таблицы Джури для остальных строк, кроме последней.
Элементы нулевого столбца (кроме последнего) равны
и они положительны. Так как выполняется необходимое условие устойчивости, последний элемент нулевого столбца также будет положительным. Следовательно, система устойчива.
Оценка качества дискретных систем. Показатели качества в переходном режиме
Качество дискретных систем управления определяется так же, как и качество непрерывных систем, и для его оценки можно использовать все ранее введенные при рассмотрении непрерывных систем показатели качества в переходном и установившемся режимах или их аналоги.
Прямые показатели качества — время регулирования и перерегулирование — определяются по переходной характеристике. Ее можно построить по дискретной переходной функции , соединяя дискретные точки плавной кривой (рис. 8.1).
Рассмотрим вычисление переходной функции. Так как -изображение от единичной решетчатой функции
то -изображение переходной функции
где — передаточная функция относительно входа
и выхода
.
Пусть изображение переходной функции имеет вид
По определению -преобразования
Поэтому значения переходной функции можно найти, разложив
в ряд Лорана путем деления числителя
на знаменатель
по правилу деления многочленов. При этом в многочленах
и
слагаемые должны располагаться в порядке убывания степени
.
Задача №8.1.
Определить значения переходной функции дискретной системы с передаточной функцией
Решение:
-изображение переходной функции
Произведя деление числителя на знаменатель по правилу деления многочленов, для первых пяти слагаемых получим
Отсюда имеем:
Если разность между степенями знаменателя и числителя равна , то первый член разложения
в ряд Лорана будет иметь степень
. Поэтому первые
значений
будут равны нулю:
Другой способ вычисления переходной функции основан на формуле разложения, которая определяется следующим образом: если все полюса функции
(т.е. корни уравнения
= 0) простые и не равны нулю, то
где
Начальные значения при
и
при
.
Задача №8.2.
Определить переходную функцию , если
-изображение имеет вид
Решение:
В данном случае
Производная
полюсами являются
И в соответствии с формулой (8.1)
Начальное значение , так как степень числителя меньше степени знаменателя. Если
имеет кратные полюса, то полюсу
кратности
в формуле разложения соответствует слагаемое, определяемое предельным соотношением
Если среди полюсов имеется нулевой (
= 0), то при вычислении соответствующего этому полюсу слагаемого следует пользоваться формулой (8.2) и в том случае, когда этот полюс является простым.
Задача №8.3.
Определить переходную функцию , если ее
-изображение имеет вид
Решение:
В данном случае
Производная
полюсами являются
Слагаемое, соответствующее нулевому полюсу ( = 0), в соответствии с формулой (8.2) определяется следующим образом:
Полюс =1 имеет кратность 2 и ему соответствует слагаемое
Полюс = —1 является простым и ему соответствует слагаемое (8.1)
Таким образом, имеем
Начальное значение .
Вычисление переходной функции между точками съема сигналов . Функция определяет значения переходной функции в моменты съема сигнала
. Значения переходной функции
в промежуточные моменты времени можно определить по структурной схеме (рис. 8.2 б), которая получается из исходной (рис. 8.2 а) подключением на выходе звена чистого запаздывания. Из этих схем имеем
Так как
и при единичном входном воздействии и нулевых начальных условиях
Соответственно для -изображений переходных функций
и
имеем:
где
Задача №8.4.
Пусть в дискретной системе (рис. 8.2, a) и период- следования импульсов
= 0,2. Требуется определить решетчатую функцию, которая принимает значения переходной функции
в моменты
Решение:
Искомой функцией будет , где
В данном случае
Подставив эти выражения в (8.3), получим
Отсюда в соответствии с формулой (8.1)
Особенности переходного процесса дискретных систем.
В непрерывных линейных системах переходная функция всегда принимает установившееся значение при . Однако возможны линейные дискретные системы, в которых переходный процесс полностью заканчивается за конечное число шагов, т. е. существует такое положительное число
, что
Если выполняется это условие, то переходный процесс называется оптимальным, а система, в которой происходит такой процесс, — оптимальной (по переходному процессу) системой.
Условие оптимальности системы (по переходному процессу). В системе с передаточной функцией вида
переходный процесс заканчивается за конечное число шагов, если
Задача №8.5.
Замкнутая дискретная система состоит из фиксатора нулевого порядка и НЧ с передаточной функцией период
= 0,1. Определить параметр
, при котором переходный процесс будет оптимальным.
Решение:
При фиксаторе нулевого порядка передаточная функция формирующего звена
Поэтому передаточная функция приведенной НЧ
и передаточная функция разомкнутой дискретной модели
Передаточная функция замкнутой системы
Отсюда в соответствии с формулой (8.4) для оптимального получаем
Показатели качества в установившемся режиме
Наиболее полной характеристикой качества в установившемся режиме является установившаяся ошибка . Ее можно найти по
-изображению
по формуле
Здесь определяется по рекуррентной формуле
Коэффициенты называются коэффициентами ошибки и являются числовыми показателями качества в установившемся режиме.
Статические и астатические системы. Система называется статической, если статическая ошибка отлична от нуля, и астатической, если статическая ошибка равна нулю. Статическая ошибка — это установившаяся ошибка при постоянных внешних воздействиях. Система является астатической и обладает астатизмом -го порядка, если первые
коэффициентов равны нулю, а (
+ 1)-й коэффициент ошибки отличен от нуля:
Из формул (8.56) и (8.5в) следует, что коэффициент позиционной ошибки и для нахождения остальных коэффициентов необходимо вычислять производные от передаточной функции
.
Однако для астатической системы коэффициенты до первого отличного от нуля можно определить по формуле
Задача №8.6.
Задающее воздействие , передаточная функция ошибки относительно задающего воздействия
период = 0,05. Определить установившуюся ошибку.
Решение:
Так как
и установившаяся ошибка
Найдем коэффициенты ошибки и
. В соответствии с формулами (8.56) и (8.5в)
Так как = 0, то
можно определить по формуле (8.6):
Установившаяся ошибка
Структура астатических систем. Дискретная система (рис. 8.6) будет астатической, если передаточная функция ДФ (регулятора) включает множитель 1/( — 1) или НЧ (а не ПНЧ) содержит интегрирующее звено. Порядок астатизма системы равен сумме числа последовательно соединенных интегрирующих звеньев в НЧ и показателю степени
— 1 в знаменателе ДФ.
Задача №8.7.
Пусть в дискретной системе (рис. 8.6)
Задающее воздействие
Определить установившуюся ошибку.
Решение:
В данном случае установившаяся ошибка
Передаточная функция регулятора содержит в знаменателе множитель — 1 в первой степени, и НЧ содержит одно интегрирующее звено. Поэтому система обладает астатизмом второго порядка и
Следовательно,
Кстати готовые задачи на продажу по тау тут.
Синтез дискретных систем
Метод полиномиальных уравнений. Пусть задана передаточная функция приведенной НЧ и соответственно известна дискретная передаточная функция неизменяемой части
Из заданных требований к качеству синтезируемой системы получена желаемая передаточная функция . Требуется синтезировать регулятор, при котором передаточная функция
синтезированной системы (рис. 9.1) равна желаемой:
При синтезе регулятора нужно позаботится о том, чтобы он был физически осуществим и синтезированная система была грубой. Условие физической осуществимости регулятора, состоящее в том, что следствие не может предшествовать причине, будет выполнено, если степень числителя его передаточной функции не превышает степень ее знаменателя.
Условие грубости будет нарушено, если передаточная функция неизменяемой части содержит нули или полюса вне единичного круга, и они входят в передаточную функцию регулятора. И это условие накладывает определенные ограничения на выбор желаемой передаточной функции, что в общем случае исключает возможность задания желаемой передаточной функции произвольно. Поэтому обычно задаются желаемым характеристическим полиномом синтезируемой системы.
Разложим числитель и знаменатель передаточной функции неизменяемой части на два множителя, один из которых содержит нули внутри единичной окружности, другой — на и вне единичной окружности:
Здесь — полиномы, нули которых расположены внутри единичной окружности;
— полиномы, нули которых расположены на и вне единичной окружности. Искомая передаточная функция регулятора имеет вид
где показатель степени множителя
определяется требуемым порядком астатизма,
и
— неопределенные полиномы, которые определяются из полиномиального уравнения
Здесь — желаемый характеристический полином (знаменатель желаемой передаточной функции).
Обозначим степень произвольного полинома через
. Тогда условие физической осуществимости можно записать в виде (9.2)
Полиномиальное уравнение (9.3) разрешимо, если число неизвестных (коэффициентов полиномов и
не меньше числа уравнений, получаемых приравниванием коэффициентов при одинаковых степенях в уравнении (9.3). И так как число неизвестных равно
, а число уравнений
условие разрешимости полиномиального уравнения принимает вид
Степень желаемого характеристического полинома должна удовлетворять соотношению
а также неравенству
где — степень знаменателя передаточной функции неизменяемой части.
Порядок синтеза системы управления методом полиномиальных уравнений можно сформулировать следующим образом.
- Разложить полиномы числителя и знаменателя передаточной функции неизменяемой части на два множителя, один из которых имеет нули внутри единичной окружности, другой — на и вне единичной окружности. Если указанные полиномы не имеют нулей на и вне единичной окружности, то положить
и
; если они не имеют нулей внутри единичной окружности, то приравнять
и
постоянному множителю этих полиномов.
- Исходя из требований к качеству синтезируемой системы в переходном режиме и порядку астатизма выбрать характеристический полином синтезируемой системы
и число
. Степень полинома
должна удовлетворять условию (9.7).
- Из соотношений (9.4)-(9.6) определить степени неопределенных полиномов
и
и записать их с неизвестными коэффициентами.
- Подставить полученные неопределенные полиномы в полиномиальное уравнение (9.3) и определить их коэффициенты.
- Подставить найденные полиномы
и
в формулу для передаточной функции регулятора (9.2).
Для того чтобы синтезируемый регулятор был более простым, степени полиномов ,
и
должны быть как можно меньшими.
Задача №9.1.
Передаточная функция неизменяемой части
Требуется синтезировать регулятор, при котором статическая ошибка равна нулю и переходный процесс заканчивается за конечное число шагов.
Решение:
В данном случае
и соответственно
Степени
Так как статическая ошибка должна быть равна нулю, положим = 1. Условие (9.7) принимает вид
Минимально допустимой является степень , и для того чтобы переходный процесс закончился за конечное число шагов, полагаем характеристический полином
. Из условия (9.6)
, а из условий (9.4) и (9.5)
. Поэтому полагаем
Подставив их в полиномиальное уравнение (9.3), получим
или, после раскрытия скобок и приведения подобных членов,
Приравняв коэффициенты при одинаковых степенях, находим
Решив эту систему, получим
Следовательно,
Подставив эти выражения для и
, а также выражения для
и
в (9.2), получим искомое решение
Лекции по курсу «Управление Техническими Системами», читает Козлов Олег Степанович на кафедре «Ядерные реакторы и энергетические установки», факультета «Энергомашиностроения» МГТУ им. Н.Э. Баумана. За что ему огромная благодарность.
Данные лекции только готовятся к публикации в виде книги, а поскольку здесь есть специалисты по ТАУ, студенты и просто интересующиеся предметом, то любая критика приветствуется.
В это части будут рассмотрены:
2.9. Использование обратных преобразований Лапласа для решения уравнений динамики САР (звена).
2.10. Весовая и переходная функции звена (системы).
2.11. Определение переходного процесса в системе (САР) (звене) через весовую и переходную функции.
2.12. Mетод переменных состояния.
2.13. Переход от описания переменных «вход-выход» к переменным состояния.
Попробуем применить, полученные знания на практике, создавая и сравнивая расчетные модели в разных видах. Будет интересно познавательно и жестко.
2.9. Использование обратных преобразований Лапласа для решения уравнений динамики САР (звена)
Рассмотрим динамическое звено САР изображенное на рисунке 2.9.1
Рис. 2.9.1 — Звено САР
Предположим, что уравнение динамики имеет вид:
где:
— постоянные времени;
— коэффициент усиления.
Пусть известны отображения:
Найдем изображения для производных:
Подставим полученные выражения в уравнение динамики и получим уравнение динамики в изображениях:
B(s) — слагаемое, которое определяется начальными условиями, при нулевых начальных условиях B(s)=0.
W(s) — передаточная функция.
Передаточной функцией САР (звена) называется отношение изображений выходного сигнала к входному воздействию при нулевых н.у.
После того, как в явном виде найдено изображение для неизвестной выходной величины, нахождение оригинала не представляет сложностей. Либо по формуле Хэвисайда, либо разложением на элементарные дроби, либо по таблице из справочника.
Пример
Построить выходной сигнал звена САР при единичном входном воздействии и нулевых начальных условиях, если уравнение динамики звена имеет следующий вид:
начальные условия:
входное воздействие:
— единичное ступенчатое воздействие.
Выполним преобразование Лапласа:
Подставим в уравнение динамики и получим уравнение динамики в изображениях:
Для получения выходного сигнала из уравнения в изображениях выполним обратное преобразования Лапласа:
Рисунок 2.9.2 График переходного процесса.
2.10. Весовая и переходная функции звена (системы).
Определение: Весовой функцией звена (системы) называется реакция системы при нулевых н.у. на единичное импульсное воздействие.
Рисунок 2.10.1 Весовая функция.
Определение: Переходной функцией звена (системы) при н.у. называется реакция на единичное ступенчатое воздействие.
Рисунок 2.10.2 Переходная функция.
Рисунок 2.10.3 Пример весовой функции.
Рисунок 2.10.4 Пример перходной функции.
На этом месте можно вспомнить, что преобразование Лапласа это интеграл от 0 до бесконечности по времени (см. предыдущий текст), а импульсное воздействие при таком интегрировании превращается в 1 тогда в изображениях получаем что:
Передаточная функция играет роль изображения реакции звена или системы на единичное импульсное воздействие.
Рисунок 2.10.5 Весовая функция как передаточная в изображениях.
Рисунок 2.10.6 Ступенчатое воздействие.
Для единичного ступенчатого воздействия преобразование Лапласа тоже известно (см. предыдущий текст):
тогда в изображениях получаем, что реакция системы на ступенчатое воздействие, рассчитывается так:
Реакция системы на единичное ступенчатое воздействие рассчитывается обратным преобразованием Лапласа:
2.11. Определение переходного процесса в системе (САР) (звене) через весовую и переходную функции. Формула Дюамеля-Карсона
Предположим, что на вход системы поступает произвольное воздействие x(t), заранее известное. Найти реакцию системы y(t), если известны входное воздействие x(t) и весовая функция w(t).
Рисунок 2.11 Демонстрация расчета по формуле Дюамеля-Карсона
Решение.
Представим, что входное воздействие представляет собой последовательность прямоугольных импульсов до времени t и ступеньки высотой x(t) в момент времени t. см.рис. 2.11 Для каждого импульса мы можем записать реакцию системы через весовую функциию:
где:
— значение отклика по завершению предыущего импульса;
— время завершения текущего импульса;
— значение весовой функции в начале текущего импульса.
Тогда для определения занчения отклика в произвольный момент времени необходимо сложить все импульсы и ступенчатое воздействие в момент времени t:
Переходя к пределам
получаем интеграл:
если перейти от t к бесконечности мы получим формулу интеграла Дюамеля-Карсона, или по другому «интеграла свертки» который обеспечивает вычисление оригинала функции по произвдению изображения двух функций:
где — вспомогательное время
Для вывода аналогичной зависмости от переходной функции вспомним что изображение весовой и переходной функции связаны соотношением: запишем выражение изображения для отклика в операторной форме:
Используя интеграл свертки получаем, что при известной переходной функции (h(t)) и известному входному воздействию х(t) выходное воздействие рассчитывается как:
2.12. Mетод переменных состояния.
До этого мы рассматривали системы с одной передаточной функцией, но жизнь всегда сложнее и как правило в системах есть несколько передаточных функций несколько входных воздейстий и несколько реакций системы. (см. рис. 2.12.1)
Рисунок 2.12.1 Моногомерная система автоматического управления.
В этом случае наиболее удобной формой пердставления систем для их анализа и расчета оказался метод переменных состояния. Для этого метода, вместо передаточных функций связывающих вход с выходом используются дополнительные переменные состояния, которые описывают систему. В этом случае можно говорить, что состояние системы — это та минимальная информация о прошлом, которая необходима для полного описания будущего поведения (т.е. выходов) системы, если поведение ее входов известно. см. рис. 2.12.2
Рисунок 2.12.2 Перменные состояния в многомерной системе.
В методе состояний, производные всех переменных состояния, в общем случае зависит от всех переменных и всех входных воздействия, и могут быть записаны в представленной ниже системы обыкновенных дифференциальных уравнений (ОДУ) первой степени. Эта система уравнений называю системой ОДУ в форме Коши:
Выход из системы зависит от переменных состояния и, в общем случае от входных воздействий и описывается следующей системой уравнений:
где:
n — количество перемнных состояния,
m — количество входных воздействий,
p — количество выходных переменных;
Данная система уравнений может быть записана в матричной форме:
где:
— вектор входа (или вектор управления);
— вектор столбец производных переменных состояния;
— вектор столбец переменных состояния;
— вектор выхода;
— собственная матрица системы [n x n],
— постоянные коэффициенты;
— матрица входа [n x m],
— постоянные коэффициенты;
— матрица выхода а [p x n],
— постоянные коэффициенты;
— матрица обхода [p x m],
— постоянные коэффициенты;
В нашем случае почти всегда все элементы матрицы D будут нулевыми: D = 0.
Такое описание системы позволяет с одной стороны стандартным образом описывать различные технические системы. Явная формула для расчета производных позволяет достаточно просто осуществлять численное интегрирование по времени. И это используется в различных программах моделирования
Другое использование данного представления для простых систем, описанных в переменных «вход-выход», зачастую позволяет устранить технические трудности, связанные с решением ОДУ высокой степени.
Еще одним преимуществом данного описания, является то, что уравнения в форме Коши можно получить из законов физики
Пример решения задачи в форме коши.
Рассмотрим задачу моделирования гидравлического привода, при следующих условиях:
Дано:
Цилиндрический плунжер диаметром 10 мм, с приведенной массой 100 кг, работает на пружину жесткостью 200 Н/мм и демпфер с коэффициентом вязкого трения — 1000 Н/(м/с). Полость начальным объемом 20 см3 соединяется с источником давлния дросселем диаметром диаметр которого 0,2 мм. Коэффициент расхода дросселя 0.62. Плотность рабочей жидкости ρ = 850 кг/м3.
Определить:
Перемещение дросселя, если в источнике давление происходит скачек 200 бар. см. рис. 2.12.13
Рисунок 2.12.3 Гидравлическая система.
Уравенение движение плунжера:
Где:
– площадь плунжера,
– жесткость пружины,
– коэффициент вязкого трения, p – давление в камере.
Поскольку дифференциальное движения это уравнение второго порядка, превратим его в систему из двух уравнений первого порядка, добавив новую переменную — скорость , тогда
Уравнение давления в камере, для упрощения принимаем что изменениям объема камеры из-за перемещения плунжера можно пренебречь:
Где: Q – расход в камеру, V — объем камеры.
Расход через дроссель:
Где: f– площадь дросселя,
– давление в источнике, p – давление в камере.
Уравнение дросселя не линейное, по условию задачи, давление входное изменяется скачком, от 0 до 200 бар, проведем линеаризацию в окрестности точки давления 100 бар тогда:
Подставляем линеаризованную формул расхода в формулу давления:
Таким образом общая система уравнений в форме Коши, для рис 2.12.3 привода принимает вид:
Матрицы A, B, С, В для матричной формы системы уравнений принимают вид:
Проверим моделированием в SimInTech составленную модель. На рисунке 2.12.13 представлена расчетная схема содержащая три модели:
1 — «Честная» модель со всеми уравнениями без упрощений.
2 — Модель в блоке «Переменные состояние» (в матричной форме).
3 — Модель в динамическом блоке с линеаризованным дросселем.
Рисунок 2.12.4 Расчетная схема .
Все условия задачи задаются как глобальные константы проекта, в главном скрипте проекта, там же расчитываются на этапе инициализации расчета, площади плунжера и проходного сечения дросселя см. рис. 2.12.5:
Рисунок 2.12.5 Глобальный скрипт проекта.
Модель на внутреннем языке программирования представлена на рис. 2.12.6. В данной модели используется описание модели в форме Коши. Так же выполняется учет изменения объема дросселя на каждом шаге расчета, за счет перемещения плунжера (Vk = V0+Ap*x.)
Рисунок 2.12.6 Скрипт расчета модели в форме Коши.
Модель в матричном форме задается с использованием глобальных констант в виде формул. (Матрица в SimInTech задается в виде последовательности из ее столбцов) см. рис. 2.12.7
Рисунок 2.12.7 Настройка блока расчета системы уравнений в пременных состояния в матричной форме.
Результаты расчета показывают, что модель в матричной форме и модель на скриптовом языке в форме Коши, практически полностью совпадают, это означает, что учет изменения объема полости практически не влияют на результаты. Кривые 2 и З совпадают.
Процедура линеаризация расхода через дроссель вызывает заметное отличие в результатах. 1-й график c «честной» моделью дросселя, отличается от графиков 2 и 3. (см. рис. 2.12.8)
Рисунок 2.12.8 Результаты расчета трех моделей гидравлического плунжера.
Сравним полученные модели, с моделью созданной из библиотечных блоков SimInTech, в которых учитываются так же изменение свойств реальной рабочей жидкости — масла АМГ-10. Сама модель представлена на рис. 2.12.9, набор графиков на рисунке 2.12.10
Рисунок 2.12.9 Модель демпфера из библиотечных блоков.
Рисунок 2.12.10 Результаты рассчета моделей демпфера. График 4 — модель из библиотечных блоков.
На графиках видно, что уточненная модель отличается от предыдущих, однако погрешность модели составлят наших упрощенных моделей составляют примерно 10%, в лишь в некоторые моменты времени.
2.13. Переход от описания переменных «вход-выход» к переменным состояния и обратно
Рассмотрим несколько вариантов перехода от описания «вход-выход», к переменным состояния:
Вариант прехода зависит от правой части уравнения с переменными «вход-выход»:
2.13.1. Правая часть содержит только b0*u(t)
В этом варианте, в уравнениях в правой части отсутствуют члены с производными входной величины u(t). Пример с плунжером выше так же относится к этому варианту.
Что бы продемонстрировать технологию перехода рассмотрим следующее уровнение:
Для перехода к форме Коши ведем новые переменные:
И перепишем уравнение относительно y»'(t):
Используя эти переменные можно перейти от дифференциального уравнения 3-го прядка, к системе из 3-х уравнений первого порядка в форме Коши:
Соотвественно матрицы для матричного вида уравнений в переменных сосотяния:
2.13.2. Правая часть общего вида
Более сложный случай, когда в уравнениях есть производные от входных воздействий и уравнение в общем случае выглядит так:
Сделаем преобразования: перейдем к уравнениям динамики в изображениях:
Тогда можно представить уравнение в изображениях в виде:
где:
Разделим уравнение в изображениях на произведение полиномов , получим:
Где: — некоторая комплексная величина (отношение двух комплексных величин). Можно считать, что
отображение величины
. Тогда входная величина может быть в изображениях представлена как:
Вренемся к оригиналу от изображений получим: ,
где: — дифференциальный оператор.
А это дифференциальное уравнение n-го порядка мы можем преобразовать к системе из n дифференциальных уравнений первого порядка, как это мы делали выше:
Таким образом, мы получили систему уравнение в форе Коши, относительно переменных состояния :
А регулируемую величину (выход системы) мы так же можем выразить через эти переменные, в изображениях:
Перейдем от изображения к оригиналам:
Если обозначить вектор , то мы получим уравнения переменных состояниях в матричной форме, где D = 0:
Пример:
Рисунок 2.13.1 Передаточная функция.
Имеется передаточная функция (рис. 2.13.1) в изображениях :
Необходимо преобразовать передаточную функцию к системе уравнений в форме Коши
В изображения реакция системы связана с входным воздействие соотношением:
Разделим в последнем правую и левую часть на произведения , и введем новую перменную
:
Полиномы N(s) и L(s) равны:
Перейдем в последнем выражении от изображения к оригиналам и ведем новые переменные (состояния):
Переходим от уравнения третьего порядка к системе трех уравнений первого порядка:
Или в матричной форме:
Для получения второго матричного уравнения воспользуемся соотношением для новых переменных в отображениях:
Перейдем от изображений к оригиналу:
Таким образом второе уравнение матричной системы выглядит так:
Проверим в SimInTech сравнив передаточную функцию и блок переменных состояния, и убедимся, что графики совпадают см. рис. 2.13.2
Рисунок 2.13.2 Сравнение переходного процеса у блока передаточной функции и блока переменных состояния.
Продолжение лекций находится здесь:
3. ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ ЗВЕНЬЕВ И СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ (РЕГУЛИРОВАНИЯ).
3.1. Амплитудно-фазовая частотная характеристика: годограф, АФЧХ, ЛАХ, ФЧХ.
3.2. Типовые звенья систем автоматического управления (регулирования). Классификация типовых звеньев. Простейшие типовые звенья.
3.3. Апериодическое звено 1–го порядка (инерционное звено). На примере входной камеры ядерного реактора.
3.4. Апериодическое звено 2-го порядка.
3.5. Колебательное звено.
3.6. Инерционно-дифференцирующее звено.
3.7. Форсирующее звено.
3.8. Инерционно-интегрирующее (звено интегрирующее звено с замедлением).
3.9 Изодромное звено (изодром).
Полезные ссылки:
Модель демпфера из лекции можно взять здесь…
Волченко Ю.М. Теоремы операционного исчисления.
Интеграл Дюамеля и физический смысл функции веса
Лекция. «Векторно-матричные модели систем управления в непрерывном времени»
Л. С. Шихобалов. Учебное пособие «МАТРИЦЫ И ОПРЕДЕЛИТЕЛИ»
Характеристическое уравнение матрицы
Подробное описание моделирования гидравлического демпфера.
Привет, Вы узнаете про временные характеристики линейных стационарных систем переходные процессы ы решения задач, Разберем основные ее виды и особенности использования. Еще будет много подробных примеров и описаний. Для того чтобы лучше понимать что такое
временные характеристики линейных стационарных систем переходные процессы ы решения задач , настоятельно рекомендую прочитать все из категории Математические основы теории автоматического управления.
Система автоматического регулирования в процессе работы может находиться в одном из двух, качественно отличных друг от друга, состояниях: в установившемся или неустановившемся режимах.
Установившийся (нормальный) режим работы автоматической системы характеризуется неизменностью внешних условий функционирования системы.
Неустановившийся (переходный) режим работы возникает сразу же после воздействия на систему управляющего или возмущающего сигналов. Поведение регулируемой величины в неустановившемся режиме называется переходным процессом. Переходные процессы могут быть затухающими и незатухающими. При затухающем переходном процессе регулируемая величина стремится к своему состоянию в установившемся режиме.
В системе с не затухающим переходным процессом регулируемая величина либо неограниченно возрастает (чаще всего в виде расходящихся колебаний), либо остается в том же состоянии, что и в начале переходного процесса. Такая система никогда не приходит в заданный для нее установившийся режим и называется неустойчивой.
Типовые входные сигналы, используемые при исследовании автоматических систем
Для исследования переходных процессов автоматической системы нео6ходимо подать на вход системы такое воздействие, которое вывело бы ее из состояния покоя и в дальнейшем обеспечило свободное движение регулируемой величины. Такими входными воздействиями являются ступенчатая функция и импульсная ? – функция (функция Дирака) (рис. 4.1, рис. 4.2).
Единичная ступенчатая функция обозначается 1(t) и представляет собой выражение:
(4.1)
Функция 1(t) непрерывна всюду, кроме точки t=0, где она терпит разрыв первого рода, Эта функция соответствует замыканию ключа, подающего напряжение 1 в. Изображение единичной функции по Лапласу имеет вид:
(4.2)
Для записи произвольного ступенчатого сигнала используется выражение
x(t)=Al(t) (4.3)
где А – амплитуда сигнала.
Импульсной ? – функцией называется функция, равная нулю всюду, кроме начала координат, где она равна бесконечности (рис. 4.3, рис. 4.4):
(4.4)
Интеграл от этой функции по любому интервалу, содержащему начало координат, равен единице т.е.
(4.5)
Импульсная функция представляет собой математическую идеализацию предельно короткого импульса, площадь которого равна единице. Изображение импульсной функции по Лапласу имеет вид:
(4.6)
Импульсная ? – функция является производной от единичной ступенчатой функции
(4.7)
Импульсная ? – функция обладает фильтрующим свойством.
Реакции автоматической системы на единичную ступенчатую функцию и импульсную функцию называются временными характеристиками системы, так как представляют собой процесс изменения регулируемой величины во времени.
Временные характеристики
Тип переходного процесса, качество регулирования и состояние автоматической системы в установившемся режиме, определяются по виду переходной и весовой функций (характеристик) системы.
Переходной характеристикой линейной автоматической системы называется ее реакция на воздействие единичной ступенчатой функции при нулевых начальных условиях.
Переходная характеристика обозначается h(t) и представляет собой переходный процесс, вызванный скачкообразным изменением входной величины. Она может быть определена при решении дифференциального уравнения обычным или операторным методом.
(4.8)
Весовой функцией линейной автоматической системы называется ее реакция на импульсную ? — функцию при нулевых начальных условиях.
Весовая функция обозначается k(t) и представляет собой переходной процесс, вызванный мгновенным «ударным» воздействием короткого импульса на входе системы . Об этом говорит сайт https://intellect.icu . Весовая функция также находится решением дифференциального уравнения системы.
— весовая функция (4.9)
Основные свойства весовой и переходной функций. Показатели качества переходного процесса
При исследовании динамических свойств автоматической системы интерес представляют:
— поведение системы в начальный момент времени (сразу после приложения воздействия);
— характер изменения регулируемой величины в переходном режиме;
— поведение регулируемой величины при приближении к новому установившемуся состоянию;
— поведение регулируемой величины в установившемся режиме;
— длительность переходного процесса.
Информация об этом может быть получена в результате анализа весовой или переходной функций автоматической системы.
Рис. 4.5. Кривая переходной характеристики
1) Время регулирования tp –минимальное время, по истечению которого регулируемая величина будет оставаться в достаточно малой окрестности установившегося значения, т.е.
(4.10)
где ? – постоянная величина, значение которой обычно принимается 5% от .
2) Абсолютная и относительная величины перерегулирования, определяющие отклонение переходной характеристики от установившегося значения.
(4.11)
и величина относительного отклонения, выражаемая в процентах,
(4.12)
3) Время первого согласования (быстродействие) tc, определяющее время первого достижения установившегося значения h(?).
4) Частота колебаний – период колебаний переходной характеристики.
5) Число колебаний n, которое имеет переходная характеристика за время регулирования tp.
6) Время достижения первого максимума tmax.
Переходные процессы в АСУ
В любой АСУ в результате воздействия возмущающих сил, с одной стороны, и восстанавливающего действия управляющего устройства, с другой, возникает переходный процесс: переход АСУ из одного состояния в другое. Рассмотрим различные типы переходного процесса (рис. 4.6).
Решение дифференциального уравнения описывает переходной процесс y(t), характер которого определяется коэффициентом. Рассмотрим переходные процессы, соответствующие различным значениям.
- ?<- 1. В данном случае система не может восстановить равновесное состояние, значение управляемой координаты все больше отклоняется от заданного. Такой переходный процесс называется расходящимся монотонным (апериодическим) (рис. 4.6, а), а система неустойчивой (идет процесс накопления энергии из внешней среды).
- -1 < ? < 0. Характеристики системы те же, что и в предыдущем случае, но переходный процесс колебательный (рис. 4.6, б).
- 0 <? < 1. Переходная функция h(t) та же, что и в случае II, но при ? < 0. При этом система возвращается в равновесное состояние, а значение управляемой координаты приближается к заданному. Такой переходный процесс называется сходящимся колебательным, а система устойчивой (происходит отдача энергии во внешнюю среду) (рис. 4.6, в).
Рис. 4.6
4. ? > 1. Характеристика системы та же, что и в III случае, но переходный процесс монотонный (апериодический) (рис. 4.6, г).
5. ?= 0; В системе устанавливается периодическое движение, процесс называется колебательным незатухающим, система находится на границе устойчивости (рис. 4.6, д). Она является замкнутой (консервативной), автономной от внешней среды.
Все рассмотренные колебания (II, III и V случаи) относятся к классу свободных, их параметры А и ? зависят от начальных условий, т. е. от привнесенной энергии. Для случаев II и III функция h(t)?h(t+T), где Т- период колебаний, и следовательно, эти колебания непериодические. Периодические колебания наблюдаются только в случае V
Примеры решения задач
Пример 1. Дано уравнение в дифференциальном виде:
Определить переходную и весовую характеристики. Построить графики.
Решение:
Переходная характеристика определяется по формуле (4.8):
A, B, C определяем с помощью метода неопределенных коэффициентов
Приравняем между собой числители двух дробей
Весовая функция Лапласа
Внимание!
При решении используется метод неопределенных коэффициентов (из курса высшей математики). Метод неопределенных коэффициентов применяют при интегрировании рациональных дробей. Рациональной дробью называется дробь вида
,
где P(x) и Q(x) являются многочленами. Рациональная дробь называется правильной, если степень числителя P(x) ниже степени знаменателя Q(x); в противном случае дробь называют неправильной. Любая правильная рациональная
дробь может быть единственным образом представлена в виде суммы простых рациональных дробей. Простые рациональные дроби:
Если при решении дифференциальных уравнений:
- хотя бы один из корней положителен, то система будет неустойчивой,
- если хотя бы один корень равен 0 – система находится на грани устойчивости,
- если хотя бы один корень отрицательный или комплексно-сопряженный – система устойчива.
Вывод: По полученной переходной функции можно сделать вывод о том, что система неустойчива.
Пример 2. Задан оператор линейной стационарной системы
- Определить переходную и весовую характеристики.
- Построить графики.
Вывод: По полученной переходной функции можно сделать вывод о том, что система неустойчива.
На этом все! Теперь вы знаете все про временные характеристики линейных стационарных систем переходные процессы ы решения задач, Помните, что это теперь будет проще использовать на практике. Надеюсь, что теперь ты понял что такое временные характеристики линейных стационарных систем переходные процессы ы решения задач
и для чего все это нужно, а если не понял, или есть замечания,
то нестесняся пиши или спрашивай в комментариях, с удовольствием отвечу. Для того чтобы глубже понять настоятельно рекомендую изучить всю информацию из категории
Математические основы теории автоматического управления
Сообщение от Andrea13
Откуда берется формула tрег=3*tau?
Обычно все переходные процессы выглядят как затухающие колебания
Причём b почти равен 1
Примем, что колебания затухли 5% от амлитуды изменения (y(t -> бескон)-y(0)). Согласно той формуле эта амплитуда равна a, ну или 1 (если расчёт вести по выражению внутри скобок, ведь a — лишь масштабный коэффициент).
Итак, нужно, чтобы
При b = 1 (ну, или близком к 1), значением cos можно пренебречь, считая его равным 1 (макс. значение). Можешь мысленно представить возможные варианты с высокой частотой cos и с очень низкой частотой относительно времени затухания.
Итак, надо найти лишь
Т.е.
Потом это 3*tau впитывается в кровь и забываешь о выводе формулы.
Добавлено через 14 минут
—————————————————————————————
Сообщение от Andrea13
Почему тау = 1/20?
Ну раз я решил, что tрег=3*tau, то и вычисляю это tрег.
Итак tрег=3*tau=3*(1/20)=0.15
Вроде бы логично (да и графику из Calc не противоречит)
Добавлено через 3 минуты
———————————————————————————————
Сообщение от Andrea13
При подсчете частоты, что такое w и откуда оно оберется? Я не понял, почему там 120.
Это из тригонометрии — пересчёт частоты и угловой частоты.
120 — это множитель из выражения для cos
cos(
120t+0.2)
Добавлено через 1 минуту
—————————————————————————
Я всё упростил — лишь бы не считать производные и их нули.
Попробуй сам решить — и ты поймёшь откуда что взялось.
Добавлено через 18 минут
————————————————————————
Я там нестрого принял, что максимальное перерегулирование будет на 1-м колебании. Это можно упрощённо принять при малом их количестве за время затухания (когда период и время затухания одного порядка — как в этом случае).
Если бы колебаний было больше — то я прошёлся бы по их значениям на каждом периоде или оценил бы по огибающей (т.е. по экспоненте). Если меньше — пока не знаю — но что-нибудь бы придумал.