
Теория автоматического управления. Часть 1
.pdf
71
Теперь запишем 1(s)=1+W3 (s) W4 (s). При этом разрушается первая обратная связь с передаточной функцией −W1 W2 W5 и в составе
1остаётся только передаточная функция второго контура обратной связи с передаточной функцией −W3 W4 . Соответственно получаем 2 (s)=1, поскольку при этом разрушаются обе обратные связи. Подставляя найденные передаточные функции в формулу (2.3.1), получаем итоговую передаточную функцию
W(s)= |
|
|
W1(s) W2 (s) (1+W3(s) W4 (s))+W1 |
(s) W3 |
(s) |
|
. |
|||||
1+W |
(s) W |
(s) W |
(s)+W |
(s) W |
(s)+W |
(s) W |
(s) W |
(s) W |
(s) W (s) |
|||
1 |
2 |
5 |
3 |
4 |
1 |
2 |
|
5 |
3 |
4 |
|
На практике часто эквивалентные преобразования и формула Мейсона применяются в комплексе.
2.3.2. Уравнения и передаточные функции одноконтурной САУ
Рассмотрим наиболее распространенную функциональную схему САУ, реализующую регулирование по ошибке (см. рис. 1.7). Для простоты будем считать все воздействия скалярами. В результате разбиения отдельных блоков схемы на звенья направленного действия, получения передаточных функций этих звеньев и необходимых упрощений получаем структурную схему, приведённую на рис. 2.11, а.
Здесь использованы следующие обозначения:
Wр (s)= U((s))
Εs
управления),
Wou (s)= Y((s))
U s
воздействию,
– передаточная функция регулятора (устройства
– передаточная функция объекта по регулирующему

72
Wof (s)= Y((s)) – передаточная функция объекта по возмущающему
F s
воздействию.
|
|
|
|
|
|
f(t) F(s) |
|
|
|
|
|
|
|
|
Wof (s) |
g(t) |
ε(t) |
|
|
|
u(t) |
|
y(t) |
|
|
|
W (s) |
||||
W |
|
(s) |
|||||
G(s) |
- Ε(s) |
р |
U(s) |
ou |
Y(s) |
||
|
|
|
|||||
|
|
|
|
|
|
|
а |
|
|||
|
|
|
|
f(t) F(s) |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
Wof (s) |
|
|
|
|
|
|
|
|
|
y(t) |
|
g(t) |
ε(t) |
|
u(t) |
||||
W (s) |
|||||||
Y(s) |
|||||||
G(s) |
- Ε(s) |
|
U(s) |
||||
|
|
б
Рис. 2.11. Структурная схема одноконтурной системы.
Все упомянутые выше передаточные функции представляют собой дроби, в которых числители и знаменатели являются полиномами относительно переменной s, при этом для физически реализуемых устройств степень числителя всегда меньше степени знаменателя. То есть, можно записать:
W (s)= |
b |
р0 |
sk + b |
р1 |
sk−1 |
+ ... + b |
рk |
= |
B |
р |
(s) |
, (k < r); |
(2.3.3) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
a |
р0sr + a |
р1sr−1 |
|
|
|
Aр (s) |
|||||||||||||||
р |
+ ... + aрr |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W (s)= |
|
bou0sq + bou1sq−1 |
+ ... + bouq |
|
= |
B |
|
(s) |
, (q < p); |
(2.3.4) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
ou |
|
|
|||||||
|
a s p + a s p−1 |
+ ...+ a |
|
|
A |
|
(s) |
||||||||||||||
ou |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ou0 |
ou1 |
|
|
|
|
oup |
|
|
|
ou |
|
|
|
|
|
|
|
|
|
|
|
73 |
|
|
|
|
|
|
|
W (s)= |
b |
0 |
sl |
+ b |
sl−1 |
+ ... + b |
|
|
B |
(s) |
, (l < d). |
|
||
|
of |
|
of 1 |
|
of l |
= |
of |
|
|
(2.3.5) |
||||
a |
|
sd + a |
sd −1 |
+ ... + a |
|
A |
(s) |
|||||||
of |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
of 0 |
|
|
of 1 |
|
|
of d |
|
of |
|
|
|
|
С учётом введённых выше обозначений дифференциальное уравне- |
||||||||||||||
ние, например, регулятора имеет вид |
|
|
|
|
|
|
|
|||||||
(aр0 pr + aр1 pr−1 + ...+ aрr )u(t)= (bр0 pk + bр1 pk−1 + ...+ bрk )ε(t). |
(2.3.6) |
Принимая во внимание принцип суперпозиции и правила соединения звеньев, схему рис. 2.11, а можно преобразовать в схему рис. 2.11, б. На этом рисунке W (s)= Wр (s) Wou (s).
Случай 1 . Положим f (t)= 0, g(t)≠ 0 и разорвем обратную связь.
Тогда связь выходной (регулируемой) величины и входного (задающего) воздействия выражается дифференциальным уравнением
y(t)= W(p)g(t).
Функция W (s) называется передаточной функцией разомкнутой системы (по задающему воздействию). Передаточная функция W(s) согласно обозначениям (2.3.3), (2.3.4) представляет собой правильную рациональную дробь
W(s)= W (s) W (s)= |
B |
р (s) |
|
B |
(s) |
|
b sm + b sm−1 |
+ ...+ b |
|
B(s) |
|
||||||
|
|
|
|
|
ou |
|
|
= |
0 |
1 |
m |
= |
|
|
, (2.3.7) |
||
A |
|
(s) |
A |
(s) |
a sn + a sn−1 |
+ ...+ a |
A(s) |
||||||||||
р |
ou |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
р |
|
|
|
ou |
|
|
0 |
1 |
n |
|
|
|
|
где m = k + q, n = r + p, при этом m < n.
Знаменатель передаточной функции (2.3.7), то есть полином A(s) – это характеристический полином разомкнутой системы.
Случай 2 . Пусть по-прежнему f (t)= 0, g(t)≠ 0 , но обратная связь теперь замкнута. Передаточную функцию, связывающую вход и выход системы, можно получить или по формуле соединения звеньев обратной связью, или по формуле Мейсона (с учётом того, что передаточная функция звена в обратной связи равна единице). Дифференциальное уравнение, свя-

74
зывающее задающее воздействие и регулируемую величину, будет тогда иметь вид
y(t)= |
|
|
W(p) |
g(t)= Φ(p)g(t). |
|
1 |
+W(p) |
||||
|
|
В уравнении (2.3.7) функция
Φ(s)= 1 WW(s()s)
+
(2.3.8)
(2.3.9)
называется передаточной функцией замкнутой системы (по задающему воздействию) или главным оператором системы. Передаточная функция Φ(s) – правильная рациональная дробь и с учётом обозначений (2.3.7) может быть представлена в виде
|
W(s) |
|
B(s) |
B(s) |
B(s) |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
b sm + b sm−1 |
|
|
|
||||||
Φ(s)= |
|
A(s) |
+ ...+ b |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
1 |
m |
|
|||
|
= |
|
|
|
= |
|
= |
|
= |
|
|
.(2.3.10) |
|||||
1+W(s) |
|
|
B(s) |
A(s)+ B(s) |
D(s) |
a sn + a sn−1 + ...+ a |
|
||||||||||
|
|
|
1+ A(s) |
|
|
|
0 |
1 |
|
n |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Дифференциальное уравнение системы (2.3.8) с учётом выражения (2.3.10) примет вид
(2.3.11) Знаменатель передаточной функции Φ(s), то есть полином D(s),
называется характеристическим полиномом замкнутой системы.
Случай 3 . Теперь пусть f (t)≠ 0, g(t)= 0 и обратная связь опять разомкнута. Дифференциальное уравнение, связывающее входное (в данном случае возмущающее) воздействие и выходную величину, будет иметь вид
y(t)= Wof (p)f (t).
Случай 4 . При прежних условиях, что и в случае 3 ( f (t)≠ 0, g(t)= 0) замкнём обратную связь. Полученное уравнение имеет вид

|
|
|
|
|
|
|
75 |
|
|
y(t)= |
|
Wof |
(p) |
|
f (t)= Φ f y |
(p)f (t), |
(2.3.12) |
||
1 |
+W(p) |
||||||||
|
|
|
|
где Φ f y (s) – передаточная функция замкнутой системы по возмущающему воздействию – получена по формуле Мейсона.
Случай 5 . В общем случае, когда в системе есть и задающее и возмущающее воздействие, в соответствии с принципом суперпозиции по-
лучаем |
|
y(t)= Φ(p)g(t)+ Φ f y (p)f (t). |
(2.3.13) |
Если в системе имеется несколько возмущающих воздействий, второе слагаемое в выражении (2.3.13) превращается в сумму
y(t)= Φ(p)g(t)+ ∑Φ fk y (p)fk (t).
k
Часто в САР интерес представляет величина ошибки регулирования ε(t)= g(t)− y(t). Уравнение замкнутой системы относительно ошибки бу-
дет иметь вид |
|
|
|
|
|
|
|
|||
ε(t)= |
|
|
1 |
|
g(t)− |
Wf y |
(p) |
|
f (t)= Φgε (p)g(t)+ Φ fε (p)f (t), (2.3.14) |
|
1 |
+W(p) |
1+W(p) |
||||||||
|
|
|
где передаточная функция замкнутой системы по ошибке относительно задающего воздействия Φgε (s) и передаточная функция
замкнутой системы по ошибке относительно возмущающего воздействия Φ fε (s)могут быть получены по формуле Мейсона, и с учётом
формулы (2.3.9) задаются выражениями:
Φgε (s)= 1+ 1 ( ) =1− Φ(s),
W s
Φ fε (s)= −Wf y((s)) = (Φ(s)−1)Wf y (s).
1+W s
Передаточная функция по ошибке
(2.3.15)
(2.3.16)
Φgε (s) – неправильная (степень
числителя равна степени знаменателя) дробь и согласно (2.3.7) равна

76
() 1 A(s)
Φgε s = 1+W(s) = A(s)+ B(s).
2.3.3. Линейные законы регулирования.
Закон регулирования или, в более общем случае, закон управления
–это, по сути, алгоритм преобразования входных воздействий регулятора (устройства управления) в управляющие воздействия на объект управления. Линейные законы определяются линейными дифференциальными уравнениями. В системе, реализующей принцип регулирования по ошибке (см. рис. 2.9, а), уравнение регулятора имеет вид (2.3.6)
(aр0 pr + aр1 pr−1 + ...+ aрr−1 p + aрr )u(t)= (bр0 pk + bр1 pk−1 + ...+ bрk )ε(t)
Если свободный член левой части не равен нулю aрr ≠ 0, то полином в левой части Aр (p)= aр0 pr + aр1 pr−1 + ...+ aрr−1 p + aрr определяет
свободное движение (переходной процесс) в регуляторе.
Если aрr = 0, то полином в левой части можно представить в виде
Aр (p)= aр0 pr + aр1 pr−1 + ...+ aрr−1 p = (aр0 pr−1 + aр1 pr−2 + ...+ aрr−1 )p = Aр1(p)p
и тогда уравнение регулятора будет
Aр1(p)u(t)= Bр (p)ε(t).
p
Свободное движение регулятора будет в этом случае определяться полиномом AР1(p).
В общем случае в полиноме левой части нулю могут равняться несколько (пусть их будет ν) последних коэффициентов
aрr = aрr−1 = ... = aрr−ν +1 = 0. |
|
|
|||
Тогда уравнение регулятора будет Aрν (p)u(t)= |
Bр (p) |
|
ε(t) |
и пере- |
|
pν |
|||||
|
|
|
|||
ходной процесс определится полиномом Aрν (p). |
|
|

77
В понятие алгоритма управления переходной процесс не входит (он краткосрочный по сравнению с процессом управления), поэтому линейные законы регулирования задаются общей формулой
u(t)= |
Bр (p) |
ε(t). |
(2.3.17) |
pν
В зависимости от числа ν и от вида полинома в выражении (2.3.17) получаем различные законы регулирования. Рассмотрим основные виды таких законов.
Позиционный (пропорциональный) закон . В этом случае ν=0, Bр (s)= KП и формула (2.3.17) превратиться в u(t)= KП ε(t). Таким
образом, управление пропорционально величине ошибке (отсюда и название закона). Регулятор, реализующий этот закон, называется пропорциональным регулятором или П-регулятором.
Интегральный закон . Простейший интегральный закон получим, если ν=1, Bр (p)= KИ . Тогда уравнение (2.3.17) будет иметь вид
u(t)= KИ p−1 ε(t)
или, учитывая, что p-1 означает операцию интегрирования в дифференциальных уравнениях,
u(t)= KИ ∫ε(t)dt .
Из последней формулы ясно, что управление пропорционально интегралу от ошибки, поэтому соответствующий регулятор носит название интегрального регулятора или И-регулятора.
Изодромный закон . Положим ν=1 и Bр (p)= KИ + KП p . Под-
ставляя эти значения в формулу (2.3.17), получим
u(t)= KИ + KП p ε(t)= KП ε(t)+ KИ p−1ε(t). p
Изодромный закон содержит две составляющие – пропорциональную и интегральную, поэтому этот закон называют ещё пропорционально-
78
интегральным, а регулятор, реализующий этот закон, называют пропор- ционально-интегральным или ПИ-регулятором.
Пропорционально -дифференциальный закон . Из названия следует, что в этом законе присутствуют пропорциональная и дифференциальная составляющие. При этом ν=0, а Bр (p)= KП + KД p. Тогда
u(t)= KП ε(t)+ KД pε(t).
Соответствующий этому закону регулятор называется пропорци- онально-дифференциальным или ПД-регулятором.
Наконец, если в законе регулирования есть все три составляющие – пропорциональная, интегральная и дифференциальная, то получаем про- порционально-интегрально-дифференциальный закон и соответственно ПИД-регулятор. В данном случае связь управляющего воздействия u(t) и ошибки ε(t) задаётся формулой
u(t)= KП ε(t)+ KИ p−1ε(t)+ KД pε(t).
Рассмотрены только самые простые виды законов. В более сложных случаях в составе закона могут присутствовать производные второго или более высокого порядка, двойные, тройные и т.д. интегралы. Необходимо понимать, что дифференциальный закон в чистом виде, то есть когда в законе имеется только производная от ошибки, не применяется, поскольку в этом случае управляющее воздействие u(t) не зависит от величины ошибки (она может быть как угодно большая, что недопустимо), а от скорости её изменения.
2.4. Описание САУ в пространстве состояний
Метод пространства состояний имеет ряд преимуществ по сравнению с классическим (частотным) методом:
-удобство решения задач на ЦВМ;
-единообразие описания одномерных и многомерных систем;
79
- возможность применения к некоторым типам нелинейных и нестационарных систем.
В общем случае состояние системы можно определить как минимальную информацию о системе, необходимую для определения (при известной входной функции) ее выхода, а также ее состояния в будущем. Таким образом, переменные состояния системы – это такие переменные, знание значений которых в некоторый начальный момент времени t0 позволяет определить поведение системы в текущий момент времени t>t0 (естественно, если известны входные воздействия системы на интервале (t0,t)).
2.4.1. Уравнения состояния
Как известно [1], линейные уравнения состояния, записанные в век- торно-матричной форме, в общем виде имеют вид
xɺ(t)= A(t)x(t)+ B(t)g(t), |
(2.4.1) |
|
y(t)= C(t)x(t)+ D(t)g(t). |
||
|
В выражении (2.4.1) использованы обозначения:
– n-мерный вектор-столбец переменных
состояния системы;
g(t)= (g1(t), g2 (t),...gp (t))T – p-мерный вектор-столбец входных воз-
действий;
y(t)= (y1(t), y2 (t),...yl (t))T – l-мерный вектор-столбец регулируемых величин;
A(t)= [aik (t)] – (n×n)-мерная матрица системы (основная или собственная параметрическая матрица системы);
B(t)= [bik (t)] – (n×p)-мерная матрица связи (структура этой матрицы определяет характер связи входных воздействий с переменными состояния);
C(t)= [cik (t)] – (l×n)-мерная матрица связи (именно, связи переменных состояния с выходными (регулируемыми) воздействиями);

80
D(t)= [dik (t)] – (l×p)-мерная матрица связи входных переменных непосредственно с выходными переменными*.
Система, описываемая уравнениями (2.4.1), – нестационарная система (система с переменными параметрами). Если параметры системы со временем не меняются, то элементы матриц не зависят от времени, и уравнения (2.4.1) упрощаются и обретают вид
xɺ(t)= Ax(t)+ Bg(t), |
(2.4.2) |
|
y(t)= Cx(t)+ Dg(t). |
||
|
По сути выражения (2.4.2) представляют собой две системы уравнений. Первая называется собственно уравнениями состояния, а вторая – уравнениями выхода.
2.4.2. Стандартная форма уравнений состояния
Переменные состояния могут быть заданы различным образом. Один из распространённых методов (особенно просто применять этот метод в системах с одним входом и одним выходом) заключается в переходе к фазовым координатам системы.
Пусть имеется уравнение системы, например, в виде (2.3.11). Для простоты положим в этом уравнении полином в правой части равным единице. Кроме того, положим коэффициент при старшей производной y(t) также равным единице a0 =1. Это всегда можно проделать, поделив почленно все слагаемые на a0. Тогда получим
(pn + a pn−1 |
+ ... + a |
n |
)y(t)= g(t). |
(2.4.3) |
1 |
|
|
|
Зададим переменные состояния следующим образом:
* Часто для реальных систем матрица D является нулевой матрицей, так что непосредственная связь входа с выходом отсутствует.