Материал: Polyakov_Tu_Metoda

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

© К.Ю. Поляков, 2008

Можно показать (сделайте это самостоятельно), что любой регулятор второго порядка с интегратором может быть представлен в форме ПИД-регулятора:

C(s) =

a

s2 +a s + a

0

C(s) = K +

K

I +

K

D

s

 

.

2

1

 

 

 

 

s

T s +

1

 

 

s(b s +b )

 

 

 

 

 

 

1 0

 

 

 

 

D

 

 

 

 

7.3. Метод размещения полюсов

Один из простых методов синтеза регулятора – размещение полюсов передаточной функции замкнутой системы, которые во многом определяют ее динамику, например, быстродействие и степень затухания колебаний (см. разд. 6.8). Смысл в том, чтобы разместить эти полюса в заданных точках комплексной плоскости с помощью специально выбранного регулятора. Эта задача сводится к решению системы линейных уравнений.

Пусть передаточная функция объекта задана в виде отношения полиномов

P(s) =

n(s)

=

 

n1s +n0

.

d(s)

 

s2 +d s +n

 

 

 

 

 

Выберем регулятор вида

 

 

 

 

 

 

 

1

0

 

 

 

nC (s)

 

= a1s +a0

 

 

C(s)

=

 

,

 

 

 

 

 

 

 

d

C

(s)

 

b s +b

 

 

 

 

 

 

 

1

0

 

 

где a0 , a1,b0 и b1 – неизвестные коэффициенты, которые нужно определить. Характеристиче-

ский полином замкнутой системы равен

(s) = n(s) nC (s) +d(s) dC (s) = (n1s + n0 )(a1s +a0 ) +(s2 +d1s + d0 )(b1s +b0 ) . = b1s3 +(n1a1 +d1b1 +b0 )s2 +(n0a1 + n1a0 + d0b1 + d1b0 )s +n0a0 +d0b0

Предположим, что мы хотим выбрать регулятор так, чтобы разместить корни полинома (s) в заданных точках, то есть добиться выполнения равенства

(s) = s3 +δ2 s2 +δ1s +δ0 ,

где δi (i = 0,...,2) – заданные числа. Приравнивая коэффициенты при одинаковых степенях s в последних двух равенствах, получаем

s3 : b =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s2 : n a + d b +b = δ

2

 

 

 

 

 

 

1

1

1

1

 

0

 

 

 

 

 

s1 : n a +n a

0

+ d

0

b + d b = δ

1

 

 

0

1

1

 

 

 

 

1

 

1

0

 

s0 :

n a

+d

b

= δ

0

 

 

 

 

 

 

или в матричном виде

 

0

0

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 1 0 a

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

n1

0 d1

 

 

 

1

a0

=

δ2

.

 

 

n1

d0

 

 

 

d1

 

 

 

 

 

 

 

n0

 

 

 

b1

δ1

 

0 n0

0 d0 b0 δ0

Решение уравнения имеет вид

 

0 0 1 0 1

 

 

 

a

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

= n1

0

 

d1

1

 

 

δ

2 .

 

 

 

 

n1

d0

 

 

 

 

 

 

 

b1

 

n0

d1

 

δ1

 

b0

 

0 n0

0 d0

 

δ0

Конечно, квадратная матрица в этом выражении (она называется матрицей Сильвестра) должна быть обратима. Можно доказать, что она действительно обратима тогда и только тогда, когда полиномы n(s) и d(s) не имеют общих корней, то есть передаточная функция объекта

71

© К.Ю. Поляков, 2008

P(s) несократима. В противном случае общий корень этих полиномов неизбежно будет корнем характеристического полинома (s) .

Кроме того, для того, чтобы количество неизвестных коэффициентов было равно числу уравнений, порядок регулятора нужно выбирать не меньше, чем N 1, где N – порядок модели объекта управления:

N = max{deg n(s), deg d(s)},

где deg обозначает степень полинома. Иначе полученное уравнение будет разрешимо только при специально выбранном полиноме (s) .

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

7.4. Коррекция ЛАФЧХ

На протяжении многих лет самым популярным инженерным методом синтеза регуляторов был метод, основанный на использовании логарифмических частотных характеристик (ЛАФЧХ). Он основан на двух свойствах ЛАФЧХ:

1)логарифмические амплитудные и фазовые частотные характеристики для последовательного соединения двух блоков (например, регулятора и объекта управления) равны сумме ЛАЧХ и ЛФЧХ этих блоков;

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

Пусть G(s) = P(s) R(s) – передаточная функция объекта вместе с приводом, причем будет предполагать, что она не имеет неустойчивых нулей и полюсов (то есть, является минимально-

фазовой). ЛАЧХ такого расширенного объекта обозначим как L0 (ω) = 20lg

 

G( jω)

 

. Если мы

 

 

сможем каким-то образом найти желаемую ЛАЧХ LЖ (ω) , то разница между этими двумя ха-

рактеристиками – это и есть ЛАЧХ необходимого последовательного регулятора:

 

 

LC (ω) = LЖ (ω) L0 (ω) .

(51)

Таким образом, для решения задачи требуется ответить на два вопроса:

1)как выбрать желаемую ЛАЧХ LЖ ( jω) так, чтобы обеспечить устойчивость и требуемое качество замкнутой системы?

2)как получить передаточную функцию регулятора C(s) по его ЛАЧХ (51)?

Чтобы ответить на первый вопрос, вспомним типичные требования к системе управления:

устойчивость;

нулевая ошибка в установившемся режиме;

быстрый и плавный (в идеале – монотонный) переходный процесс;

подавление шумов;

робастность (нечувствительность к ошибкам модели). Эти требования нужно связать с формой ЛАЧХ.

Постоянный сигнал можно рассматривать как предельный случай гармонического (синуса), только с нулевой частотой. Поэтому для обеспечения нулевой установившейся ошибки цепочка «регулятор-объект» должна иметь бесконечное усиление на нулевой частоте, то есть передаточная функция G(s) C(s) должна содержать интегратор (вспомните принцип внутренней

модели).

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

72

© К.Ю. Поляков, 2008

1

1

Ts

Ts +1

 

Таким образом, для получения монотонного переходного процесса ЛАЧХ разомкнутой системы

должна быть похожа на ЛАЧХ интегратора – это прямая линия с наклоном –20 дБ/дек, которая

пересекает ось абсцисс на частоте ωc =1/T . Эта частота называется частотой среза. Заметим,

что для апериодического звена легко определить время переходного процесса: оно примерно

равно 3T. Таким образом, частота среза определяет время переходного процесса. Вспомним,

что устойчивость системы также определяется поведением ЛАЧХ в районе частоты среза. В ре-

зультате имеем:

 

 

 

 

 

 

 

 

 

 

устойчивость и качество переходного процесса (время, перерегулирование) определяют-

 

ся формой ЛАЧХ в районе частоты среза, где она пересекает ось Lm = 0 ; эта область на-

 

зывается областью средних частот;

 

 

 

 

 

 

 

 

 

для получения качественного переходного процесса желательно, чтобы наклон ЛАЧХ

 

около частоты среза был равен –20 дБ/дек;

 

 

 

 

 

 

 

 

 

если задано время переходного процесса tп , нужно выбиратьωc = 3 .

 

 

 

 

 

 

 

 

 

 

 

tп

 

 

Теперь разберемся с шумами и робастностью. Как мы знаем,

шумы – это высокочастот-

ные сигналы. Кроме того, обычно именно в области высоких частот характеристики объекта и

модели могут сильно расходиться. Поэтому для подавления помех и уменьшения влияния оши-

бок модели нужно по возможности уменьшать усиление системы в области высоких частот, то

есть ЛАЧХ должна резко идти вниз.

 

 

 

 

 

 

 

 

 

 

 

На рисунке показана типовая желаемая ЛАЧХ.

L

 

 

 

 

 

 

Это асимптотическая ЛАЧХ, состоящая из отрез-

m

–20 дБ/дек

 

высокие

ков. В выделенных точках стыкуются два отрезка

 

 

 

 

 

 

частоты

 

 

 

 

 

 

(подавление

разного наклона. На низких частотах она имеет на-

 

 

 

 

 

–20 дБ/дек

помех)

клон –20 дБ/дек, то есть система содержит интегра-

 

 

 

 

 

 

 

 

12-16 дБ

ωc

ω

тор, который обеспечивает нулевую ошибку в уста-

0

 

 

 

 

 

12-16 дБ

новившемся режиме.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

низкие

средние частоты

 

 

ЛАЧХ пересекает ось абсцисс под наклоном

 

 

 

 

 

 

частоты

(устойчивость,

 

–20 дБ/дек. Для обеспечения устойчивости и при-

 

 

(точность)

переходный процесс)

 

емлемого показателя колебательности ( M <1,2 ) точки излома ЛАЧХ должны находиться на

расстоянии 12-16 дБ от оси абсцисс (см. рисунок).

 

 

 

 

 

 

 

 

 

Продемонстрируем метод коррекции ЛАЧХ на простом примере. Пусть объект управле-

ния – апериодическое звено с передаточной функцией

G(s) =

1

, где T = 5 с. Передаточная

 

 

 

 

 

 

 

 

 

T0 s +1

0

 

 

 

 

 

 

 

 

 

 

 

 

функция замкнутой системы без коррекции (то есть, с регулятором C0 (s) =1) равна:

 

 

W (s) =

1

=

 

0,5

 

.

 

 

 

 

 

T s +2

 

0,5T s +

1

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

Видим, что статический коэффициент усиления W (0) = 0,5

(а не 1), так что точного отслежива-

ния входного сигнала не получается. Время переходного процесса можно приближенно подсчи-

тать как tп = 3 0,5 T0 = 7,5 c .

 

 

 

 

 

 

 

 

 

 

 

Поставим задачу следующим образом: выбрать регулятор C(s) , который обеспечивает

нулевую ошибку в установившемся режиме;

 

 

 

 

 

 

 

 

время переходного процесса около 1,5 с;

 

 

 

 

 

 

 

 

 

73

© К.Ю. Поляков, 2008

 

наклон ЛАЧХ –40 дБ/дек на высоких частотах для подавления помех.

 

 

 

Для решения используем метод коррекции

 

 

 

 

 

 

 

 

ЛАЧХ. Синяя линия на рисунке обозначает нескор-

L

 

 

 

 

 

 

 

ректированную ЛАЧХ, совпадающую с ЛАЧХ апе-

m

–20 дБ/дек

 

 

риодического звена G(s) .

 

 

 

 

 

 

 

L0

 

 

Lж ωc

ω1

ω

 

Желаемая ЛАЧХ (зеленая линия) должна иметь

 

0

 

1

 

 

 

наклон –20 дБ/дек на низких частотах, чтобы обеспе-

 

 

 

 

 

 

 

15 дБ

 

 

 

 

T0 –20 дБ/дек

 

чить нулевую статическую ошибку. Частота среза ωc

 

 

 

 

 

 

Lm

 

 

 

 

 

 

 

определяется

требуемым

быстродействием:

 

LC

 

 

ωc

= 3/ tп = 2 рад/с. Таким образом, начальный

уча-

 

 

 

 

 

 

 

 

 

 

 

ω

сток желаемой ЛАЧХ совпадает с ЛАЧХ интегри-

 

 

 

 

 

ωc

 

рующего звена с передаточной функцией

ω

c

, то есть

 

0

 

 

 

 

 

 

(ω) = 20lg ωc

 

 

 

s

 

 

 

 

 

 

 

 

 

 

L

(на низких частотах).

 

 

 

 

 

 

 

 

 

 

 

 

 

ж

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На высоких частотах нужно изменить наклон ЛАЧХ с –20 до -40 дБ/дек на частоте ω1 , где

Lm (ω) = −15 дБ. Из этого условия находим

 

 

 

 

 

 

 

 

 

 

 

 

 

20lg ωc

= −15

ωc =1015 / 20

 

 

 

ω =103 / 4 ω

c

= 5,62 2 =11,24 рад/с.

 

 

ω1

 

ω1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, мы полностью построили желаемую ЛАЧХ, удовлетворяющую требованиям к системе. Вычитая из нее исходную ЛАЧХ (без коррекции, синяя линия), получим ЛАЧХ регулятора, которая показана красной линией на нижнем графике.

Остается перейти от ЛАЧХ регулятора к его передаточной функции. На низких частотах (ω <1/T0 ) ЛАЧХ регулятора имеет наклон –20 дБ/дек и проходит через точку (ωc ;0) , то есть

C(s) = ωsc C1 (s) ,

где C1 (s) не изменяет асимптотическую ЛАЧХ на частотах, меньших 1/T0 . На частоте ω0 =1/ T0 ЛАЧХ регулятора меняет наклон с –20 дБ/дек до нуля, то есть в числитель добавляется множитель T0s +1:

C(s) = ωc (T0ss +1) C2 (s) .

Здесь C2 (s) – регулятор, не влияющий на ЛАЧХ для частот, меньших ω1 . Наконец, на частоте ω1 наклон увеличивается с нуля до –20 дБ/дек. Для того, чтобы на этой частоте «загнуть» вниз

ЛАЧХ,

нужно добавить в регулятор

апериодическое звено с постоянной времени

T =

1

 

= 0,09 c . Таким образом, окончательно

 

ω

y

1

 

 

 

 

 

1

 

 

ωc (T0s +1)

 

 

 

 

 

 

C(s) =

.

 

1

 

 

 

 

 

 

 

 

 

s(T1s +1)

 

 

На рисунке показаны переходные процессы при единичном ступенчатом входном сигнале в нескорректированной системе (синяя линия) 0,5 и в системе с полученным регулятором C(s)

(зеленая линия). Графики показывают, что

 

 

найденный

регулятор

значительно

ускорил

 

 

переходный

процесс

и обеспечил

нулевую 0

5

10

статическую ошибку (установившееся значе-

 

 

T1 , где

15 t

74

© К.Ю. Поляков, 2008

ние выхода равно 1).

Нужно отметить, что алгоритм коррекции ЛАЧХ существенно усложняется, если объект содержит неустойчивые или неминимально-фазовые звенья.

7.5. Комбинированное управление

Один из способов улучшить качество управление – изменить структуру системы, добавив в нее второй регулятор C2 (s) на входе:

x

 

 

 

 

 

привод

g объект

 

+

e

 

u

 

δ

 

 

y

C2(s)

C(s)

R (s)

 

P(s)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь

W (s) = C2 (s)C(s)R(s)P(s) . 1+C(s)R(s)P(s)

Регулятор C2 (s) не влияет на свойства контура управления (запасы устойчивости, подавление

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

и робастность, а затем сформировать нужные качества передаточной функции W (s) с помощью регулятора C2 (s) . Поскольку две передаточные функции можно изменять независимо друг от

друга, такая схема называется комбинированным управлением (или управлением с двумя степенями свободы).

В идеале мы хотим,

чтобы система точно воспроизводила сигнал x(t) на выходе

y(t) , то

есть, нужно обеспечить W (s) 1 . Для этого требуется, чтобы

 

 

 

 

 

 

C2 (s) = 1+C(s)R(s)P(s) =

1

,

(52)

 

 

 

 

W1 (s)

 

 

C(s)R(s)P(s)

 

C(s)R(s)P(s)

 

 

где W (s) =

 

 

– передаточная функция замкнутой системы с одной степенью сво-

 

 

 

1

1+C(s)R(s)P(s)

 

 

 

 

 

 

 

 

 

боды.

 

следует, что регулятор C2 (s) должен быть обратной системой (инверсией) для

Из (52)

W1 (s) . Частотная характеристика W1 ( jω) в реальных системах близка к нулю на высоких частотах, следовательно, регулятор C2 (s) должен иметь в этом частотном диапазоне огромное уси-

ление. Например, для W1 (s) = Ts1+1 получим C2 (s) = Ts +1, то есть регулятор содержит физиче-

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

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

7.6. Инвариантность

Если возмущение g можно как-то измерить, для улучшения качества системы иногда вводится третий регулятор (третья степень свободы):

75