Лекции по курсу «Управление Техническими Системами», читает Козлов Олег Степанович на кафедре «Ядерные реакторы и энергетические установки», факультета «Энергомашиностроения» МГТУ им. Н.Э. Баумана. За что ему огромная благодарность.
Данные лекции только готовятся к публикации в виде книги, а поскольку здесь есть специалисты по ТАУ, студенты и просто интересующиеся предметом, то любая критика приветствуется.
В это части будут рассмотрены:
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 Изодромное звено (изодром).
Полезные ссылки:
Модель демпфера из лекции можно взять здесь…
Волченко Ю.М. Теоремы операционного исчисления.
Интеграл Дюамеля и физический смысл функции веса
Лекция. «Векторно-матричные модели систем управления в непрерывном времени»
Л. С. Шихобалов. Учебное пособие «МАТРИЦЫ И ОПРЕДЕЛИТЕЛИ»
Характеристическое уравнение матрицы
Подробное описание моделирования гидравлического демпфера.
5.6.1
Переходная
функция.Для
получения переходной функции в качестве
стандартного сигнала используется
единичная
функция
времени (2.16). Такого рода воздействию
соответствует, например, сброс или
включение нагрузки в системах регулирования
(отказ мотора в системе регулирования).
Рис.
3.10 Переходная характеристика химического
реактора:
а
– ступенчатое воздействие; б – кривая
разгона
Переходной
функцией называется аналитическое
выражение для решения линейного
дифференциального уравнения (8) при
входном сигнале x(t) = 1(t) и нулевых начальных
условиях, т.е.
Кривой
разгона называется реакция объекта
(системы) на единичное ступенчатое
воздействие при
нулевых
начальных условиях.На практике кривая
разгона определяется экспериментальным
путем и используется в качестве ис-
ходных
данных для анализа и синтеза систем
автоматического управления исследуемом
объектом.
Здесь
следует ввести понятия прямой и обратной
задач. Прямая задача (задача Коши)
заключается в
определении
решения дифференциального уравнения
с заданными начальными условиями. В
обрат-
ной
задаче требуется восстановить вид и
коэффициенты дифференциального уравнения
по известной
интегральной
кривой, например, переходной функции.
Решение обратной задачи представляет
значительную сложность вследствие ее
некорректности и здесь существует
специальный математический аппарат.
Так, например, если предположить, что
переходная функция описывается решением
уравнения
первого порядка
5.6.2 Весовая функция
Для
получения весовой функции, ее также
называют импульсной переходной функцией,
в качестве
стандартного
сигнала используется δ -функция (2.17):
Таким
образом, весовой функцией w(t) называется
реакция системы на δ -функцию при нулевых
начальных условиях.На практике весовую
функцию в отдельных случаях можно
получить экспериментальным путем
весьма
приближенно. Считают, что на вход объекта
подана δ -функция, если время действия
импульса
намного
меньше времени переходного процесса.
Примером может служить эксперимент по
снятию
весовой
функции химического реактора
(рис.
5.4),
являющегося объектом исследования. В
качестве входного сигнала в реактор
залпом выливается порция красящего
вещества (например, чернил). Через
некоторое время это вещество появится
на
выходе, причем его концентрация
первоначально возрастает, а затем
убывает – красящее вещество
вымывается
(рис. 5.11).
Подаваемый
на вход импульс представляет собой
приближенную дельта-функцию, так как
его
площадь
отлична от единицы и равна S. Поэтому
для получения весовой функции
експериментально
снятый
переходный
процесс нормируют путем деления его
ординат на величину площади
входного воздействия S.
Рис.
5.11
Переходная характеристика химического
реактора:
а
– δ -функция; б – весовая функция
Между
временными характеристиками: переходной
и весовой функциями существует взаимное
однозначное соответствие, которое
определяется следующим образом:
Весовую
функцию можно получить и как решение
дифференциального уравнения
При
решении подобных уравнений дельта-функцию
переводят
в
начальные условия, и если n = 2, то
Соседние файлы в папке ЛЕКЦИЯ_5
- #
- #
- #
- #
- #
- #
- #
На этой странице я собрала теорию и практику, готовые задачи и подробные решения по предмету теория автоматического управления, чтобы вы смогли освежить знания.
Если что-то непонятно — вы всегда можете написать мне в 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), получим искомое решение
Макеты страниц
Функция веса системы с переменными параметрами является исчерпывающей характеристикой этой системы, и нахождение ее важно по следующим соображениям. Функция веса характеризует протекание временных процессов в системе регулирования, и по ее виду можно судить о качестве регулирования, аналогично тому, как это делалось для систем с постоянными параметрами (§ 8.4). По имеющейся функции веса можно определить время протекания переходного процесса, как характеристику быстродействия, и склонность системы к колебаниям.
Кроме того, по имеющейся функции веса можно строить процесс на выходе системы регулирования при заданных входных воздействиях, не производя при этом каждый раз полного решения исходного уравнения (13.1). В соответствии с формулами (13.9) и (13.11) для этой цели необходимо иметь сопряженные функции веса.
Ввиду сложности проблемы существующие методы позволяют пока решать задачу нахождения функции веса в численном виде. Только для систем регулирования, описываемых дифференциальными уравнениями первого и иногда второго порядков, удается решать задачу в общем виде. Поэтому в некоторых случаях приходится сложную систему с переменными параметрами приближенно сводить к более простой системе, движение которой описывается уравнением не выше второго порядка.
Следует заметить, что большинство систем регулирования с переменными параметрами относится к так называемым квазистационарным системам, или системам, параметры которых меняются сравнительно медленно. В подобных системах коэффициенты дифференциального уравнения (13.1) мало меняются в течение времени переходного процесса, определяемого временем затухания нормальной функции веса.
Дифференциальное уравнение первого порядка.
В некоторых случаях для оценки вида переходных процессов системы с переменными параметрами ее уравнение приближенно можно свести к дифференцальному уравнению первого порядка
Это уравнение имеет аналитическое решение
где
а С — постоянная интегрирования.
Пусть, например, имеется уравнение
Определим для него семейство переходных характеристик . Для единичной ступенчатой функции при уравнение (13.14) можно записать в следующем виде:
Приведем его к виду (13.12):
Далее получаем:
На основании формулы (13.13) получаем
При нулевых начальных условиях должно быть Отсюда определяется постоянная интегрирования
Окончательно получаем
Дифференцируя последнее выражение по можно получить функцию веса:
или в ином виде:
Для дифференциального уравнения (13.12) можно сразу найти функцию веса из общего решения (13.13), если положить в (13.12) входной сигнал равным единичному смещенному импульсу Проделав необходимые выкладки, получаем
где
Распространим этот результат на более общий случай записи дифференциального уравнения в виде
Приведем его к виду (13.12):
Положив получим для функции веса решение в виде
где
Рассмотрим снова в качестве примера уравнение (13.14). Приведем его к виду (13.17):
Обратившись к формуле (13.18), находим
и функцию веса
что совпадает с полученным ранее выражением.
Дифференциальное уравнение второго порядка.
Рассмотрим случай, когда дифференциальное уравнение (13.1) сводится к уравнению второго порядка
При помощи подстановки
это уравнение приводится к виду
Здесь введено обозначение
При действии единичного импульса для уравнения (13.21) получится решение которое связано с весовой функцией исходного уравнения (13.19) на основании формулы (13.20) соотношением
Если же положить то для уравнения (13.21) будет получена весовая функция которая на основании (13.9) связана с решением зависимостью
Эта зависимость на основании свойства дельта-функции может быть представлена в виде
В результате из (13.23) и (13.24) получаем
Таким образом, для отыскания функции веса необходимо предварительно решить уравнение (13.21), которое приобретает вид
с нулевыми натальными условиями: и при Полученную при решении весовую функцию и необходимо затем подставить в (13.25) и найти
Решение уравнения (13.26) может быть произведено при помощи использования функций Бесселя [118]. Для этого функция должна быть аппроксимирована отрезками прямых линий, уравнение которых сводится к виду Однако это решение является сравнительно сложным.
Рис. 13.4.
Ограничимся рассмотрением так называемого аппроксимирующего решения, которое может применяться, если функция мало изменяется относительно своего среднего большого значения (рис. 13.4). Это решение называется аппроксимацией Бриллуина — Вентцеля — Крамера [118].
Рассмотрим однородное дифференциальное уравнение
Предположим теперь, что для некоторого однородного дифференциального уравнения второго порядка получено частное решение
где
Найдем дифференциальное уравнение, которому удовлетворяет решение (13.28). Продифференцировав его дважды и исключив промежуточные переменные, получаем
Сравнивая (13.30) и (13.27), видим, что выражение (13.28) будет частным решением уравнения (13.27), если выполняется тождество
Решение уравнения (13.31) и отыскание функции является сложной задачей вследствие наличия нелинейностей в (13.31). Однако может быть найдено приближенное решение (13.31) в виде ряда, если удовлетворяются неравенства
Тогда решение (13.21) можно представить в виде
Подставляя этот ряд в (13.31), получаем формулы для определения членов ряда:
Часто можно ограничиться только первым членом ряда (13.22), что будет справедливым, если функция изменяется медленно, оставаясь в среднем большой (рис. 13.4). Тогда
При выполнении условия в качестве второго частного решения можно взять комплексно-сопряженную величину (13.29)
Тогда можно показать, что решение уравнения (13.26) будет
или, после подстановки (13.28) и (13.35),
В нредельном случае постоянства параметров Тогда . В результате из формулы (13.37) можно получить функцию веса консервативного звена
Для исходного дифференциального уравнения (13.19) на основании (13.25) и (13.37) получаем искомую функцию веса
Критерием медленности изменения функции и, следовательно, применимости полученного выражения может служить неравенство
которое получается из (13.31) и (13.34).
Метод последовательных приближений.
Рассмотрим уравнение (13.1):
Ограничиваясь случаем квазистационарных систем и полагая, что коэффициенты меняются медленно, найдем функцию веса для этого уравнения.
Переменные коэффициенты в левой части исходного уравнения представим в виде суммы постоянной и изменяющейся частей:
где переменный коэффициент, зафиксированный для момента приложения входной величины
Тогда исходное дифференциальное уравнение (13.1) можно представить в виде
где
Поскольку мы предположили, что коэффициенты меняются медленно, то функция мала по сравнению с левой частью (13.41). Эту функцию можно рассматривать как возмущение, и тогда к уравнению (13.41) можна применить метод последовательных приближений.
В уравнении (13.41) можно перейти к изображениям по Лапласу. Тогда получим
Здесь введено обозначение
Решение уравнения (13.41) или (13.44) можно записать в виде ряда
Для получения первого приближения зафиксируем переменные коэффициенты . Тогда первое приближение может быть найдено как решение дифференциального уравнения
Решение этого уравнения можно получить, используя обычные методы (см. главу 7), в том числе путем нахождения оригинала, соответствующего изображению (13.44) при
Для получения второго приближения в правую часть (13.41) или (13.44) подставляется первое приближение а в левую часть — Тогда получается уравнение с фиксированными коэффициентами для определения поправки:
Это уравнение также может быть решено с использованием преобразования Лапласа посредством нахождения оригинала изображения
где — изображение при подстановке в формулу
Повторяя этот процесс многократно, можно найти рекуррентное соотношение для определения члена ряда (13.46):
Ряд (13.46) сходится тем быстрее, чем медленнее изменяются коэффициенты Рассмотренный метод может использоваться как для нахождения
функции веса и переходной функции, так и для построения переходного процесса при любом известном воздействии .
Численно-графический метод.
Численно-графический метод Д. А. Башкирова [98] разработан также применительно к системам с переменными во времени параметрами, причем можно вводить любое переменное возмущающее или задающее воздействие и произвольные начальные условия.
Неоднородные уравнения первого порядка с переменными коэффициентами. Пусть требуется построить решение уравнения
с начальным условием при Разделив его на приведем уравнение к виду
где
Уравнение (13.51) можно решать графически, если считать Т постоянным и равным Т внутри каждого интервала времени но различным для разных интервалов. Формула для решения в этом случае будет
а процесс построения сводится к следующему. Наносим заданные кривые (рис. 13.5). Из точки Е кривой , взятой в середине первого интервала откладываем по горизонтали отрезок величина которого берется равной ординате точки Н заданной кривой т. е. тоже в середине первого интервала Полученная точка М соединяется прямой линией с заданной начальной точкой процесса А.
Рис. 13.5.
В результате получается новая точка В искомой кривой . Затем аналогично берется ордината точки откладывается в виде отрезка и проводится прямая дающая новую точку С решения
Неоднородные уравнения второго порядка с переменными коэффициентами. Требуется построить решение уравнения
которое можно записать также в виде
где
с натальными условиями при Если в правой части (13.52) имеется операторное выражение, то предварительно производим вычисление правой части и сводим ее к
Если обозначить то уравнение (13.52) разобьется на два;
где
а начальные условия будут
Формулы для решения уравнений (13.53) согласно [64, 78] будут
где
причем во второй из формул (13.55) значения берутся со сдвигом на вправо по сравнению с
Рис. 13.6.
Отсюда вытекает следующее построение. Наносим заданные кривые а также кривую ординаты которой определяются по второй из формул (13.54). Они показаны на графике (рис. 13.6, а).
На другом графике наносим заданное (рис. 13.6, б). На основании заданных начальных условий (см. выше) наносим на последнем графике точки и в середине первого интервала (как в § 7.6) еще точку А с ординатой (7.76), т. е.
Из точки в середине первого интервала на кривой откладываем вниз отрезок
(вниз, когда он положителен, и вверх, когда он отрицателен). При этом величина х берется как ордината уже имеющейся точки А, величина берется из графика а величина — из заданных начальных
условий. Из полученной точки откладываем горизонтальный отрезок
размер которого берется из графика Точку М соединяем прямой линией с точкой что дает новую точку кривой при
Из точки откладываем вниз отрезок
равный ординате точки А. Из точки проводим горизонтальный отрезок
размер которого берется из графика Точку К соединяем с точкой что дает новую точку В искомой кривой в середине второго интервала
Опишем еще второй шаг интегрирования. Из точки кривой в середине второго интервала откладывается вниз отрезок
где — ордината точки В, полученной выше; — тангенс угла наклона прямой проведенной ранее (он дает требуемое значение х). Откладываем отрезок
И проводим прямую получая при этом новую точку кривой
Из точки откладываем вниз
а затем вправо
после чего проводим прямую Это дает новую точку С искомой кривой Все описанные построения можно заменить числовыми расчетами.
Оглавление
- ПРЕДИСЛОВИЕ
- РАЗДЕЛ I. ОБЩИЕ СВЕДЕНИЯ О СИСТЕМАХ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 1.1. Понятие о замкнутых автоматических системах
- § 1.2. Классификация автоматических систем по характеру внутренних динамических процессов
- § 1.3. Примеры непрерывных автоматических систем
- § 1.4. Примеры дискретных и релейных автоматических систем
- ГЛАВА 2. ПРОГРАММЫ И ЗАКОНЫ РЕГУЛИРОВАНИЯ. АДАПТИВНЫЕ СИСТЕМЫ
- § 2.1. Программы регулирования
- § 2.2. Линейные и нелинейные законы регулирования
- § 2.3. Системы с переменной структурой
- § 2.4. Системы с самонастройкой программы (экстремальные системы)
- § 2.5. Системы с самонастройкой параметров (собственно самонастраивающиеся системы)
- § 2.6. Системы с самонастройкой структуры (самоорганизующиеся системы)
- РАЗДЕЛ II. ОБЫКНОВЕННЫЕ ЛИНЕЙНЫЕ СИСТЕМЫ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 3.1. Линеаризация уравнений
- § 3.2. О записи линеаризованных уравнений звеньев
- ГЛАВА 4. ДИНАМИЧЕСКИЕ ЗВЕНЬЯ И ИХ ХАРАКТЕРИСТИКИ
- § 4.2. Временные характеристики
- § 4.3. Частотная передаточная функция и частотные характеристики
- § 4.4. Логарифмические частотные характеристики
- § 4.5. Позиционные звенья
- § 4.6. Интегрирующие звенья
- § 4.7. Дифференцирующие звенья
- § 4.8. Неустойчивые и неминимально-фазовые звенья
- § 4.9. Звенья с модулированным сигналом
- ГЛАВА 5. СОСТАВЛЕНИЕ ИСХОДНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 5.1. Общий метод составления исходных уравнений
- § 5.2. Передаточные функции систем автоматического регулирования
- § 5.3. Законы регулирования
- § 5.4. Использование структурных схем и графов
- § 5.5. Многомерные системы регулирования
- § 5.6. Управляемость и наблюдаемость
- § 5.7. Уравнения следящей системы
- ГЛАВА 6. КРИТЕРИИ УСТОЙЧИВОСТИ
- § 6.1. Понятие об устойчивости систем регулирования
- § 6.2. Критерий устойчивости Гурвица
- § 6.3. Критерий устойчивости Михайлова
- § 6.4. Построение областей устойчивости. D-разбиение
- § 6.5. Критерий устойчивости Найквиста
- § 6.6. Определение устойчивости по логарифмическим частотным характеристикам
- § 6.7. Устойчивость двумерных систем с антисимметричными связями
- ГЛАВА 7. ПОСТРОЕНИЕ КРИВОЙ ПЕРЕХОДНОГО ПРОЦЕССА В СИСТЕМАХ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 7.2. Непосредственное решение исходного дифференциального уравнения
- § 7.3. Сведение неоднородного уравнения к однородному
- § 7.4. Использование преобразований Фурье, Лапласа и Карсона—Хевиеайда
- § 7.5. Использование вещественных частотных характеристик
- § 7.6. Использование вычислительных машин
- ГЛАВА 8. ОЦЕНКА КАЧЕСТВА РЕГУЛИРОВАНИЯ
- § 8.2. Точность в типовых режимах
- § 8.3. Коэффициенты ошибок
- § 8.4. Определение запаса устойчивости и быстродействия по переходной характеристике
- § 8.5. Приближенная оценка вида переходного процесса по вещественной частотной характеристике
- § 8.6. Корневые методы
- § 8.7. Диаграмма Вышнеградского
- § 8.8. Интегральные оценки
- § 8.9. Частотные критерии качества
- § 8.10. Чувствительность систем регулирования
- ГЛАВА 9. ПОВЫШЕНИЕ ТОЧНОСТИ СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 9.2. Теория инвариантности и комбинированное управление
- § 9.3. Неединичные обратные связи
- ГЛАВА 10. УЛУЧШЕНИЕ КАЧЕСТВА ПРОЦЕССА РЕГУЛИРОВАНИЯ
- § 10.2. Последовательные корректирующие звенья
- § 10.3. Параллельные корректирующие звенья
- § 10.4. Обратные связи
- § 10.5. Методы повышения запаса устойчивости
- § 10.6. Примеры
- ГЛАВА 11. СЛУЧАЙНЫЕ ПРОЦЕССЫ В СИСТЕМАХ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 11.2. Случайные процессы
- § 11.3. Стационарные случайные процессы
- § 11.4. Корреляционная функция
- § 11.5. Спектральная плотность стационарных процессов
- § 11.6. Канонические разложения случайных функций
- § 11.7. Прохождение случайного сигнала через линейную систему
- § 11.8. Расчет установившихся ошибок в автоматических системах
- § 11.9. Расчеты по минимуму среднеквадратичной ошибки
- ГЛАВА 12. МЕТОДЫ СИНТЕЗА СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 12.2. Корневой метод
- § 12.3. Метод корневых годографов
- § 12.4. Метод стандартных переходных характеристик
- § 12.5. Метод логарифмических амплитудных характеристик
- § 12.6. Синтез систем автоматического регулирования на основе частотных критериев качества
- § 12.7. Об оптимальном синтезе
- § 12.8. Использование классических вариационных методов
- § 12.9. Динамическое программирование
- § 12.10. Аналитическое конструирование регуляторов
- РАЗДЕЛ III. ОСОБЫЕ ЛИНЕЙНЫЕ СИСТЕМЫ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- ГЛАВА 13. СИСТЕМЫ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ
- § 13.2. Нахождение функции веса и построение переходных процессов
- § 13.3. Передаточные функции
- § 13.4. Устойчивость и качество регулирования
- § 13.5. О синтезе систем с переменными параметрами
- ГЛАВА 14. СИСТЕМЫ С ЗАПАЗДЫВАНИЕМ И СИСТЕМЫ С РАСПРЕДЕЛЕННЫМИ ПАРАМЕТРАМИ
- § 14.1. Уравнения линейных систем с запаздыванием
- § 14.2. Уравнения линейных систем с распределенными параметрами
- § 14.3. Исследование устойчивости и качества регулирования
- ГЛАВА 15. ИМПУЛЬСНЫЕ СИСТЕМЫ
- § 15.2. Использование z-преобразования
- § 15.3. Передаточные функции
- § 15.4. Устойчивость и качество импульсных систем регулирования
- § 15.5. Случайные процессы в импульсных системах
- РАЗДЕЛ IV. НЕЛИНЕЙНЫЕ СИСТЕМЫ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- ГЛАВА 16. СОСТАВЛЕНИЕ УРАВНЕНИЙ НЕЛИНЕЙНЫХ СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 16.2. Уравнения систем с нелинейностью релейного типа
- § 16.3. Уравнения систем с нелинейностью в виде сухого трения и зазора
- § 16.4. Уравнения систем с нелинейностями других видов
- ГЛАВА 17. ТОЧНЫЕ МЕТОДЫ ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ И АВТОКОЛЕБАНИЙ
- § 17.1. Фазовые траектории и метод точечных преобразований
- § 17.2. Теоремы прямого метода Ляпунова и их применение
- § 17.3. Определение автоколебаний релейных систем методом припасовывания
- § 17.4. Частотный метод В. М. Попова
- § 17.5. Исследование систем с переменной структурой
- ГЛАВА 18. ПРИБЛИЖЕННЫЕ МЕТОДЫ ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ И АВТОКОЛЕБАНИЙ
- § 18.1. Гармоническая линеаризация нелинейностей
- § 18.2. Алгебраические способы определения автоколебаний и устойчивости в нелинейных системах первого класса
- § 18.3. Примеры исследования нелинейных систем первого класса
- § 18.4. Нелинейные системы второго класса
- § 18.5. Вычисление высших гармоник и уточнение первой гармоники автоколебаний
- § 18.6, Частотный метод определения автоколебаний
- ГЛАВА 19. МЕДЛЕННО МЕНЯЮЩИЕСЯ ПРОЦЕССЫ В АВТОКОЛЕБАТЕЛЬНЫХ СИСТЕМАХ
- § 19.1. Статические и скоростные ошибки автоколебательных систем
- § 19.2. Прохождение медленно меняющихся сигналов в автоколебательных системах
- § 19.3. Гармоническая линеаризация нелинейностей при несимметричных колебаниях
- ГЛАВА 20. ОЦЕНКА КАЧЕСТВА НЕЛИНЕЙНЫХ ПРОЦЕССОВ РЕГУЛИРОВАНИЯ
- § 20.1. Приближенное исследование колебательных переходных процессов
- § 20.2. Примеры исследования колебательных переходных процессов
- § 20.3. Система с нелинейным корректирующим устройством
- § 20.4. Применение логарифмических частотных характеристик для исследования нелинейных законов рзгулирования
- ГЛАВА 21. ВЫНУЖДЕННЫЕ КОЛЕБАНИЯ НЕЛИНЕЙНЫХ СИСТЕМ
- § 21.1. Симметричные одночастотные вынужденные колебания
- § 21.2. Несимметричные вынужденные колебания с медленно меняющейся составляющей
- § 21.3. Зависимость устойчивости и качества нелинейных систем от внешних вибраций
- ГЛАВА 22. СЛУЧАЙНЫЕ ПРОЦЕССЫ В НЕЛИНЕЙНЫХ СИСТЕМАХ
- § 22.1. Статистическая линеаризация нелинейностей
- § 22.2. Простейшие случайные процессы в нелинейных системах
- § 22.3. Пример исследования влияния случайных помех на динамику нелинейной системы
- ГЛАВА 23. НЕЛИНЕЙНАЯ ОПТИМИЗАЦИЯ СИСТЕМ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- § 23.2. Синтез оптимальной системы с использованием принципа максимума
- § 23.3. Последовательная оптимизация на базе нелинейного программирования
- РАЗДЕЛ V. ЦИФРОВЫЕ И АДАПТИВНЫЕ СИСТЕМЫ АВТОМАТИЧЕСКОГО РЕГУЛИРОВАНИЯ
- ГЛАВА 24. СИСТЕМЫ РЕГУЛИРОВАНИЯ С ЦИФРОВЫМИ ВЫЧИСЛИТЕЛЬНЫМИ МАШИНАМИ
- § 24.2. О синтезе систем регулирования с ЦВМ
- § 24.3. Дискретная коррекция
- § 24.4. Периодические режимы, обусловленные квантованием по уровню
- ГЛАВА 25. АДАПТИВНЫЕ СИСТЕМЫ
- § 25.1. Системы экстремального регулирования
- § 25.2. Самонастраивающиеся системы
- ПРИЛОЖЕНИЕ 1. ТАБЛИЦА ФУНКЦИЙ
- ПРИЛОЖЕНИЕ 2. ТАБЛИЦА ИНТЕГРАЛОВ
- ЛИТЕРАТУРА
можем ввести допущение
Теперь свяжем нашу функцию с параметрами уравнения (*)
Пользуясь свойством линейности найдем изображение
левой части уравнения (*)
Теперь преобразуем правую часть уравнения (*):
приравняем левую и правую части уравнения к
изображению, это даст уравнение для получения передаточной функции:
пользуясь теоремой об дифференциальном интеграле,
запишем
уравнение
получилось алгебраическое, в котором Y(p) и X(p)
можно вынести за скобку
передаточная функция линейного звена W(p)
является дробно-рациональной функцией от p и содержит
коэффициенты дифференциального уравнения исходно-описывающего это звено
.
Ограничения:
1.
Понятие передаточной функции
справедливо только для линейных звеньев
2.
Начальные условия на
сигнал были нулевыми.
Теорема Пезу.
Если
полином вида: имеет п- корней
()
i=1…n
Аналогично
– нули переходной функции
– полюса переходной функции
2.5Временные характеристики звеньев.
Для
изучения свойств элементарных звеньев удобно использовать типовые
воздействия. По характеру отклика на эти воздействия мы можем судить об их
динамических свойствах.
1.Переходная функция
Переходная функция h(t)
есть функция, описывающая реакцию системы (звена) на единичное ступенчатое
воздействие при нулевых начальных условиях
h(t)x(t)=1(t)
График
переходной функции – кривая зависимости функции h(t) от
времени t – называют переходной или разгонной
характеристикой.
Свяжем
переходную функцию с передаточной, для этого найдем преобразование Лапласа:
L{x(t)} = L{1(t)}
=
Y(p) A(p)X(p)
=
h(t) = L-1
Например:
W(p)=
Y
М
Мы приложили постоянный ступенчатый сигнал
М(t)=1(t)
Чему будет равно
t
2.Весовая
функция
Весовой
или импульсной переходной функцией системы называют функцию, описывающую
реакцию звена на единичное импульсное воздействие при нулевых начальных
условиях; обозначают эту функцию .
Определим
связь между передаточной и весовой функцией:
– входной сигнал
– выходной сигнал
Как
видим изображение выходной функции равно передаточной функции
Рассмотрим
пример:
Нам
дана передаточная функция и переходная
Найдем
весовую функцию
Представим
графически
M(t)
t
3.Статистическая характеристика линейных звеньев
Это зависимость между установившимися весовой и динамической
характеристиками:
– статистическая
характеристика
Для линейных звеньев такая характеристика будет иметь
вид линейных звеньев
к
– коэффициент передача,
коэффициент
усиления
Пользуясь
теоремой о конечном значении, свяжем переходную функцию с к
можем
связать с к
перейдем
к изображению Лапаласа
– теореме о
конечном значении
4.Частотные характеристики
звеньев.
Важное значение при описании
линейных стационарных систем (звеньев) имеют частотные характеристики. Они
получаются при рассмотрении вынужденных движений системы при подаче на ее вход
гармонического воздействия.
Предположим к звену приложено Sinое
воздействие
– коэффициент
передачи амплитуды, амплитудно-частотная характеристика (АЧХ)
– фазовый сдвиг, фазо-частотная характеристика
(ФЧХ) -безразмерное
время
Эксперимент по прохождению
гармонического сигнала предоставляет нам 2 характеристики: АЧХ и ФЧХ
Удобно ввести комплексно-значную
частотную характеристику , которая
содержала бы в себе всю информацию об объекте, т.е. объединила АЧХ и ФЧХ
– частотная характеристика (ЧХ)
Для введения ЧХ на входе будет
обобщенный гармонический сигнал
X = ejwt =1
– обобщенный гармонический сигнал,
комплексное число интегрируемое как вектор, проекции которого дают обычные
тригонометрические функции времени.
–
переменная частота
Рассмотрим
прохождение гармонического сигнала через линейное звено, описываемое уравнением
(*).
На
выходе будет сигнал
……………………………………………
Также
будут выражаться производные для выходного сигнала
Определим
частотой характерной Ф(jw) отношение
общего выходного сигнала Y к общему входному сигналу.
Ф(jw) º в соответствии со (*)=
=
, видимо, что Ф(jw) = W(p) при условии, Р заменим на jw.
Ф(jw) = W(p)|p ® tw
Частотная
характеристика является частотно – зависимым комплексным числом, где w – круговая частота составляющих обобщенного
гармонического сигнала.
Установим
связь АЧХ, ФЧХ с ЧХ, для этого представим комплексное число
+j
– вещественная частотная
характеристика (ВЧХ)
– мнимая (МЧХ)
Тогда можно утверждать
+1