Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теория автоматического управления. Часть 1

.pdf
Скачиваний:
59
Добавлен:
05.02.2023
Размер:
1.32 Mб
Скачать

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

sk1

+ ... + b

рk

=

B

р

(s)

, (k < r);

(2.3.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

р0sr + a

р1sr1

 

 

 

Aр (s)

р

+ ... + aрr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (s)=

 

bou0sq + bou1sq1

+ ... + bouq

 

=

B

 

(s)

, (q < p);

(2.3.4)

 

 

 

 

 

 

 

 

 

 

 

 

ou

 

 

 

a s p + a s p1

+ ...+ a

 

 

A

 

(s)

ou

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ou0

ou1

 

 

 

 

oup

 

 

 

ou

 

 

 

 

 

 

 

 

 

 

 

73

 

 

 

 

 

 

 

W (s)=

b

0

sl

+ b

sl1

+ ... + 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 pr1 + ...+ aрr )u(t)= (bр0 pk + bр1 pk1 + ...+ 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 sm1

+ ...+ b

 

B(s)

 

 

 

 

 

 

ou

 

 

=

0

1

m

=

 

 

, (2.3.7)

A

 

(s)

A

(s)

a sn + a sn1

+ ...+ 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 , но обратная связь теперь замкнута. Передаточную функцию, связывающую вход и выход системы, можно получить или по формуле соединения звеньев обратной связью, или по формуле Мейсона (с учётом того, что передаточная функция звена в обратной связи равна единице). Дифференциальное уравнение, свя-

(a0 pn + a1 pn1 + ... + an )y(t)= (b0 pm + b1 pm1 + ...+ bm )g(t).

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 sm1

 

 

 

Φ(s)=

 

A(s)

+ ...+ b

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

m

 

 

=

 

 

 

=

 

=

 

=

 

 

.(2.3.10)

1+W(s)

 

 

B(s)

A(s)+ B(s)

D(s)

a sn + a sn1 + ...+ 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 pr1 + ...+ aрr1 p + aрr )u(t)= (bр0 pk + bр1 pk1 + ...+ bрk )ε(t)

Если свободный член левой части не равен нулю aрr 0, то полином в левой части Aр (p)= aр0 pr + aр1 pr1 + ...+ aрr1 p + aрr определяет

свободное движение (переходной процесс) в регуляторе.

Если aрr = 0, то полином в левой части можно представить в виде

Aр (p)= aр0 pr + aр1 pr1 + ...+ aрr1 p = (aр0 pr1 + aр1 pr2 + ...+ aрr1 )p = Aр1(p)p

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

Aр1(p)u(t)= Bр (p)ε(t).

p

Свободное движение регулятора будет в этом случае определяться полиномом AР1(p).

В общем случае в полиноме левой части нулю могут равняться несколько (пусть их будет ν) последних коэффициентов

aрr = aрr1 = ... = 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И p1 ε(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И p1ε(t). p

Изодромный закон содержит две составляющие – пропорциональную и интегральную, поэтому этот закон называют ещё пропорционально-

78

интегральным, а регулятор, реализующий этот закон, называют пропор- ционально-интегральным или ПИ-регулятором.

Пропорционально -дифференциальный закон . Из названия следует, что в этом законе присутствуют пропорциональная и дифференциальная составляющие. При этом ν=0, а Bр (p)= KП + KД p. Тогда

u(t)= KП ε(t)+ KД pε(t).

Соответствующий этому закону регулятор называется пропорци- онально-дифференциальным или ПД-регулятором.

Наконец, если в законе регулирования есть все три составляющие – пропорциональная, интегральная и дифференциальная, то получаем про- порционально-интегрально-дифференциальный закон и соответственно ПИД-регулятор. В данном случае связь управляющего воздействия u(t) и ошибки ε(t) задаётся формулой

u(t)= KП ε(t)+ KИ p1ε(t)+ KД pε(t).

Рассмотрены только самые простые виды законов. В более сложных случаях в составе закона могут присутствовать производные второго или более высокого порядка, двойные, тройные и т.д. интегралы. Необходимо понимать, что дифференциальный закон в чистом виде, то есть когда в законе имеется только производная от ошибки, не применяется, поскольку в этом случае управляющее воздействие u(t) не зависит от величины ошибки (она может быть как угодно большая, что недопустимо), а от скорости её изменения.

2.4. Описание САУ в пространстве состояний

Метод пространства состояний имеет ряд преимуществ по сравнению с классическим (частотным) методом:

-удобство решения задач на ЦВМ;

-единообразие описания одномерных и многомерных систем;

x(t)= (x1(t), x2 (t),...xn (t))T

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 pn1

+ ... + a

n

)y(t)= g(t).

(2.4.3)

1

 

 

 

Зададим переменные состояния следующим образом:

* Часто для реальных систем матрица D является нулевой матрицей, так что непосредственная связь входа с выходом отсутствует.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]