- •Гоу впо Кубанский государственный технологический университет
- •Пугачев в.И.
- •Рецензенты:
- •Isbn © Пугачев в.И., 2010
- •1.2 Основные элементы сау
- •1.4 Замкнутые и разомкнутые сау
- •1.5 Статика и динамика регулирования
- •3.2 Типовые звенья сау
- •Апериодическое звено
- •3.3 Основные законы регулирования
- •4.2 Алгебраический критерий устойчивости Гурвица
- •4.3 Частотный критерий устойчивости Михайлова
- •4.4 Критерий устойчивости Найквиста - Михайлова
- •4.5 Устойчивость сау с запаздыванием
- •5.3 Граница устойчивости и область устойчивости в плоскости одном и двух параметров
- •6.2 Интегральные оценки качества сау
- •7.2 Введение интеграла в закон регулирования
- •7.3 Создание инвариантных сау
- •7.4 Создание комбинированных сау
- •8 Цифровые системы управления
- •8.1 Способ управления с помощью эвм
- •8.2 Решетчатые функции и разностные уравнения
- •8.3 Дискретная передаточная функция
- •8.4 Получение оригинала по -преобразованию.
- •8.5 Цифровые аналоги типовых законов управления
- •8.6 Анализ цифровых систем управления
- •8.7 Анализ устойчивости цифровых систем.
- •9 Нелинейные системы
- •9.1 Основные типы нелинейностей
- •9.2 Устойчивость нелинейных систем
- •Общие положения об устойчивости нелинейных систем
- •9.3 Устойчивость релейных систем
8.4 Получение оригинала по -преобразованию.
-преобразование позволяет легко найти оригинал функции. Это связано с его особенностями, которые объясняются свойствами решетчатой функции. Если - аналитическая функция, то ее можно разложить в ряд Тейлора. Однако, если она имеет изолированные особые точки в виде полюса, то ее можно разложить в ряд Лорана, главная часть которого состоит из конечного числа членов, равного порядку полюса. Для решетчатой функции каждое значение ее есть изолированная особая точка в виде полюса, поскольку представляет интеграл от - функции, имеющей бесконечно большое значение за бесконечно малый интервал и равный значению непрерывной функции в моменты .
Разложим в ряд по отрицательным степеням:
. (8.13)
Поскольку
, если ,
, если ,
то можно вычислить коэффициенты , в разложении , которое называют преобразованием Лорана.
Умножив левую и правую часть равенства (8.12) на , получим:
.
Интегрируя левую и правую часть равенства по замкнутому контуру и учитывая свойства интеграла, получаем:
, или .
Но , поэтому оригинал решетчатой функции вычисляется по формуле:
.
Если полюсы простые, а
, где
, то
. (8.14)
Если полюсы кратные то
,
где - кратность корня,
- число полюсов.
Более простой метод получения по вытекает из разложения ее в ряд Лорана по отрицательным степеням . При этом осуществляется непрерывное деление полинома числителя на полином знаменателя. Это удобно при использовании ЭЦВМ.
Например
Численные значения решетчатой функции на выходе цифровой системы можно получить путем замены -передаточной функции уравнением в конечных разностях.
Пусть передаточная функция цифровой системы
, .
Этой передаточной функции соответствует уравнение в конечных разностях:
,
Разделив коэффициенты на и решив уравнение относительно , получают значения искомой выходной величины через ее значения в предшествующие моменты времени и значения входной величины в соответствующие моменты времени.
По этому уравнению легко реализовать алгоритм вычисления решетчатой функции на выходе цифровой системы при заданном входном воздействии.
8.5 Цифровые аналоги типовых законов управления
Рассмотрим методику получения уравнений в конечных разностях и -передаточных функций эквивалентных законам управления общепромышленных регуляторов на примере наиболее сложного ПИД -закона управления, из которого могут быть получены другие, как частные случаи.
Интегро-дифференциальное уравнение ПИД регулятора
, (8.15)
где - коэффициент усиления; - постоянная интегрирования;
- постоянная дифференцирования.
Для малых периодов квантования это уравнение можно преобразовать в разностное путем замены производной первой разностью, а интеграла - суммой.
Используя интегрирование по методу прямоугольников, выражение (8.15) можно записать так:
. (8.16)
Этот не рекуррентный алгоритм управления требует для формирования суммы помнить все предыдущие значения сигнала ошибки, поэтому для программирования на ЭВМ удобны рекуррентные алгоритмы, которые для вычисления управляющего воздействия требуют знания предыдущего значения и поправочного члена. Если из (8.16) вычесть
,
то получим рекуррентный алгоритм управления:
, (8.17)
где .
Из уравнения (8.16) легко получаются соотношения для "П" и "ПД" регуляторов, когда отсутствует суммирование рассогласований
, (8.18)
или
, (8.19)
т.е. в этом случае не требуется знание предшествующего значения управляющего воздействия для получения . Для "ПИ" -закона разностное уравнение также получается из (8.16) с учетом знака суммы, но без последнего слагаемого. Легко видеть, что в этом случае
, (8.20)
где
.
Для "И" -закона регулирования
, (8.21)
где .
Используя основные свойства (линейности и сдвига) для нулевых начальных условий:
;
,
находят для регулятора, например, реализующего ПИД -закон регулирования:
.
Взяв отношение -преобразования выходной величины к - преобразованию входной при нулевых начальных условиях, получают дискретную передаточную функцию искомого элемента системы.
. (8.22)
В общем виде передаточную функцию цифрового фильтра можно записать так:
.
Этот алгоритм может быть реализован, если .
Соотношение порядка полиномов числителя и знаменателя может быть различным: или . Обычно в регуляторах , а . Для исключения статической ошибки в алгоритме управления -го порядка необходимо иметь полюс . Поэтому структура регулятора -го порядка в простейшем случае имеет вид:
Наиболее распространенным в практике управления является алгоритм второго порядка, . При подаче на вход регулятора единичного сигнала, получаем:
;
;
;
. . . . . . . . . . . . . .
. (8.23)
В случае, когда , дискретный регулятор соответствует непрерывному ПИД -регулятору с дополнительной задержкой на один такт. Если , то можно записать следующие соотношения для параметров регулятора :
при или ;
при для или .
На рисунке 8.3 представлены переходные процессы для алгоритмов первого и второго порядков.
По аналогии с непрерывными регуляторами можно выделить следующие коэффициенты, характеризующие свойства регуляторов:
- коэффициент усиления;
- коэффициент опережения;
- коэффициент интегрирования.
Для положительного коэффициента усиления регулятора , а диапазоны изменения параметров регулятора такие: .
Рисунок 8.3 - Переходные функции цифровых регуляторов
Если вычислены оптимальные параметры непрерывного ПИД регулятора, то по ним можно посчитать оптимальные параметры цифрового регулятора.
Обозначим .
Тогда
, (8.24)
При этом дискретная передаточная функция цифрового фильтра принимает вид: .
При малых периодах квантования и высокой разрядности АЦП и ЦАП параметры цифрового фильтра можно находить как для непрерывного регулятора, после чего их следует пересчитать по формулам (8.24) для цифрового регулятора.