ТЕОРЕТИЧЕСКИЕ ОСНОВЫ АВТОМАТИЗИРОВАННОГО УПРАВЛЕНИЯ
.pdfТЕОРЕТИЧЕСКИЕ ОСНОВЫ АВТОМАТИЗИРОВАННОГО УПРАВЛЕНИЯ
Конспект лекций
Б. Р. АНДРИЕВСКИЙ
Балтийский государственный технический университет (БГТУ) «Военмех» Кафедра И3
Санкт-Петербург
2008
СОДЕРЖАНИЕ
1 |
Динамические и статические системы. Понятие состояния динамических |
|
|
|
систем |
7 |
|
2 |
Уравнения состояния линейных систем. Линеаризация уравнений состоя- |
|
|
|
ния |
|
11 |
|
2.1 |
Линейные системы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
11 |
|
2.2 |
Линеаризация уравнений состояния . . . . . . . . . . . . . . . . . . . . . . . . |
12 |
|
2.3 |
Примеры уравнений состояния систем . . . . . . . . . . . . . . . . . . . . . . |
16 |
|
|
2.3.1 Электротехнические устройства . . . . . . . . . . . . . . . . . . . . . |
16 |
2.3.2Угловое движение искусственного спутника Земли. . . . . . . . . . . 18
3 |
Передаточные функции и их определение по уравнениям состояния |
20 |
||
|
3.1 |
Передаточные функции линейных систем . . . . . . . . . . . . . . . . . . . . |
20 |
|
|
3.2 |
Алгоритмы вычисления передаточных функций . . . . . . . . . . . . . . . . . |
21 |
|
|
3.3 |
Примеры перехода к передаточным функциям от уравнений состояния . . |
23 |
|
|
3.4 |
Частотные характеристики . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
26 |
|
|
|
3.4.1 |
Частотные характеристики непрерывных систем . . . . . . . . . . . . |
26 |
|
|
3.4.2 |
Частотные характеристики дискретных систем . . . . . . . . . . . . . |
28 |
|
|
3.4.3 |
Частотные характеристики цифровых систем реального времени . . . |
29 |
4 |
Преобразование базиса. Инвариантность передаточной функции |
31 |
||
5 |
Канонические формы уравнений состояния. Диагональная и жорданова |
|
||
|
формы |
|
35 |
|
|
5.1 |
Диагональная форма. Простые вещественные собственные числа . . . . . . |
35 |
|
|
5.2 |
Вещественная диагональная форма. Простые мнимые собственные числа . . |
36 |
|
|
5.3 |
Общий случай. Вещественная форма Жордана . . . . . . . . . . . . . . . . . |
39 |
|
6 |
Управляемая и наблюдаемая канонические формы |
42 |
6.1Управляемое каноническое представление . . . . . . . . . . . . . . . . . . . . 42
6.2Наблюдаемое каноническое представление . . . . . . . . . . . . . . . . . . . . 44
2
7Преобразование уравнений состояния к каноническому виду. Преобразо-
вание к диагональной и блочно-диагональной формам |
46 |
|
7.1 |
Простые вещественные собственные числа . . . . . . . . . . . . . . . . . . . |
47 |
7.2 |
Простые мнимые собственные числа . . . . . . . . . . . . . . . . . . . . . . . |
48 |
8 Преобразование уравнений состояния к управляемой и наблюдаемой ка- |
|
|
ноническим формам |
50 |
8.1О возможности преобразования матрицы к форме Фробениуса . . . . . . . . 50
8.2Управляемое каноническое представление . . . . . . . . . . . . . . . . . . . . 51
8.3 Наблюдаемое каноническое представление . . . . . . . . . . . . . . . . . . . |
52 |
9 Определение уравнений состояния по передаточной функции |
55 |
9.1Управляемое каноническое представление . . . . . . . . . . . . . . . . . . . . 56
9.2 |
Наблюдаемое каноническое представление . . . . . . . . . . . . . . . . . . . |
57 |
9.3 |
Блочно-диагональная форма . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
59 |
9.4 |
Жорданова форма . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
60 |
9.5 |
Случай систем с несколькими входами и выходами . . . . . . . . . . . . . . |
61 |
10 Фазовые траектории и фазовые портреты линейных систем |
64 |
10.1Определения и основные свойства фазовых траекторий и фазовых портретов 64
10.2Поле фазовых скоростей. Классификация особых точек . . . . . . . . . . . . 66
10.2.1Вектор фазовой скорости . . . . . . . . . . . . . . . . . . . . . . . . . . 66
10.2.2 Состояния равновесия системы . . . . . . . . . . . . . . . . . . . . . . 66
10.2.3 Декомпозиция пространства состояний . . . . . . . . . . . . . . . . . . 68
10.3Виды фазовых портретов для систем второго порядка . . . . . . . . . . . . . 71
|
10.3.1 |
Фазовые портреты при диагональной (жордановой) форме матрицы A |
74 |
|
10.3.2 |
Фазовые портреты при канонической форме фазовой переменной . . |
76 |
11 Решение уравнений состояния. Формула Коши |
79 |
||
11.1 |
Решение однородного уравнения . . . . . . . . . . . . . . . . . . . . . . . . . |
79 |
|
11.2 |
Решение неоднородного уравнения . . . . . . . . . . . . . . . . . . . . . . . . |
81 |
|
11.3 |
Свойства переходной матрицы . . . . . . . . . . . . . . . . . . . . . . . . . . |
82 |
3
11.4 Вычисление функции веса . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
11.5Определение начального состояния по начальному значению выхода и его производных . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
12 Дискретные модели непрерывных систем |
86 |
12.1Постановка задачи дискретизации . . . . . . . . . . . . . . . . . . . . . . . . . 86
12.2Формулы перехода к разностным уравнениям . . . . . . . . . . . . . . . . . . 88
13 |
Методы вычисления матричной экспоненты |
90 |
|
|
13.1 |
Точные методы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
90 |
|
13.2 |
Приближенные методы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
93 |
|
13.3 |
Вычисление матрицы Q в общем случае . . . . . . . . . . . . . . . . . . . . . |
95 |
14 |
Дискретные модели для различных видов входного процесса |
97 |
|
|
14.1 |
Смещенное z-преобразование . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
97 |
|
14.2 |
Прямоугольные импульсы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
98 |
14.3Экспоненциальные импульсы . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
14.4Треугольные импульсы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
14.5Подстановочные формулы для вычисления передаточной функции дискретной модели . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
15 |
Управляемость и наблюдаемость линейных систем |
104 |
|
|
15.1 |
Основные определения . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
104 |
|
15.2 |
Критерии управляемости . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
107 |
|
15.3 |
Критерии наблюдаемости. Теорема дуальности . . . . . . . . . . . . . . . . . |
111 |
16 |
Оценивание состояния объекта и возмущений |
114 |
|
|
16.1 |
Постановка задачи оценивания . . . . . . . . . . . . . . . . . . . . . . . . . . . |
114 |
|
16.2 |
Наблюдатели состояния . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
115 |
16.2.1Наблюдатель полного порядка . . . . . . . . . . . . . . . . . . . . . . . 115
16.2.2Наблюдатели пониженного порядка . . . . . . . . . . . . . . . . . . . . 119
16.3 Оценивание возмущений . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
4
17 Синтез модальных и терминальных регуляторов |
125 |
17.1 Задача модального управления . . . . . . . . . . . . . . . . . . . . . . . . . . |
125 |
17.2Модальное управление по состоянию объекта . . . . . . . . . . . . . . . . . . 125
17.3Модальное управление по выходу объекта. Теорема разделения . . . . . . . . 127
17.4 Терминальное управление . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
17.5Примеры систем модального и терминального управления . . . . . . 134
17.5.1Стабилизация углового движения ИСЗ с компенсацией возмущений . 134
17.5.2Возбуждение колебаний в цепочке осцилляторов . . . . . . . . . . . . 135
18 Уравнения и характерные свойства нелинейных систем |
138 |
|
18.1 |
Общие сведения о нелинейных системах . . . . . . . . . . . . . . . . . . . . . |
138 |
18.2 |
Уравнения нелинейных звеньев и систем . . . . . . . . . . . . . . . . . . . . |
141 |
18.3Особенности процессов в нелинейных системах . . . . . . . . . . . . . . . . . 147
18.3.1Принцип суперпозиции . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
18.3.2Сепаратрисные поверхности . . . . . . . . . . . . . . . . . . . . . . . . 148
18.3.3Предельные циклы. Автоколебания . . . . . . . . . . . . . . . . . . . . 149
|
18.3.4 |
Состояния равновесия. Отрезки покоя . . . . . . . . . . . . . . . . . . |
152 |
|
18.3.5 |
Неединственность решений. Пересечение траекторий . . . . . . . . . . |
153 |
|
18.3.6 |
Скользящие режимы . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
154 |
|
18.3.7 |
Влияние внешних воздействий . . . . . . . . . . . . . . . . . . . . . . . |
155 |
19 Методы исследования нелинейных систем |
157 |
||
19.1 |
Задачи и методы теории нелинейных систем . . . . . . . . . . . . . . . . . . . |
157 |
|
19.2 |
Методы фазового пространства . . . . . . . . . . . . . . . . . . . . . . . . . . |
159 |
19.2.1Метод фазовой плоскости . . . . . . . . . . . . . . . . . . . . . . . . . . 159
19.2.2Метод точечных отображений . . . . . . . . . . . . . . . . . . . . . . . 160
19.2.3Условия существования предельных циклов для систем второго порядка161
19.3Метод гармонической линеаризации (гармонического баланса) . . . . . . . . 161
19.3.1Основные положения. «Свойство фильтра» . . . . . . . . . . . . . . . 162
19.3.2Коэффициенты гармонической линеаризации . . . . . . . . . . . . . . 163
19.3.3Уравнение гармонического баланса . . . . . . . . . . . . . . . . . . . . 166
5
|
19.3.4 Пример. Исследование генератора колебаний . . . . . . . . . . . . . . |
169 |
20 Метод функций Ляпунова |
175 |
|
20.1 |
Основные определения . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
175 |
20.2 |
Устойчивость множеств и частичная устойчивость . . . . . . . . . . . . . . . |
177 |
20.3 |
Функции Ляпунова . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
182 |
20.4 |
Устойчивость непрерывных систем . . . . . . . . . . . . . . . . . . . . . . . . |
183 |
20.5 |
Устойчивость дискретных систем . . . . . . . . . . . . . . . . . . . . . . . . . |
189 |
20.6 |
Примеры . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
191 |
21 Методы теории абсолютной устойчивости |
199 |
|
21.1 |
Задача абсолютной устойчивости . . . . . . . . . . . . . . . . . . . . . . . . . |
199 |
21.2 |
Круговой критерий . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
201 |
21.3 |
Критерий В. М. Попова . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |
202 |
22 Исследование скользящих режимов. Метод эквивалентного управления 203
22.1Понятие о скользящих режимах . . . . . . . . . . . . . . . . . . . . . . . . . . 203
22.2Определение движения в скользящем режиме . . . . . . . . . . . . . . . . . . 205
22.3Методы определения движения в скользящем режиме . . . . . . . . . . . . . 206
|
22.4 |
Метод эквивалентного управления . . . . . . . . . . . . . . . . . . . . . . . . |
208 |
23 |
Системы с переменной структурой в задаче управления |
210 |
|
24 |
Системы с переменной структурой в задаче оценивания состояния |
216 |
|
25 |
Методы адаптивного управления |
219 |
|
|
25.1 |
Задача адаптивного управления . . . . . . . . . . . . . . . . . . . . . . . . . |
219 |
|
25.2 |
Структура адаптивных систем управления . . . . . . . . . . . . . . . . . . . . |
220 |
|
25.3 |
Методика решения задач адаптивного управления . . . . . . . . . . . . . . . |
221 |
6
Лекция 1
1Динамические и статические системы. Понятие состояния динамических систем
Любая система, в том числе и система управления, состоит из совокупности подсистем
(звеньев). Звенья могут различаться по характеру реакций на входное воздействие. С этой
точки зрения все звенья могут быть разделены на статические (безынерционные) и дина-
мические (инерционные) . Рассмотрим отличительные особенности в поведении и математическом описании систем одного и другого типов.
Статические системы 1 обладают мгновенной реакцией на входное воздействие. Более
существенным свойством таких систем является то, что их реакция на входное воздействие не зависит от предыстории, от поведения системы в прошлом, а также от предыдущих
значений входа.
Математически это можно описать следующим образом.
Обозначим через u(t), y(t) вход и выход системы в момент t. У статической системы для
каждого t выход y(t) можно определить однозначно по значению u(t) в тот же момент вре-
мени. Для этой цели служит статическая характеристика y = f(u) или y = f(u, t) (для
нестационарных систем). В соответствии с ней получаем y(t) = f(u(t)). Никакой другой
дополнительной информации не требуется. 2
Иначе обстоит дело с динамическими системами. Их особенностью является то, что
для определения y(t) недостаточно информации об u(t) в тот же момент времени. Выход-
ной сигнал зависит также от предыстории изменения входа и, кроме того, совокупности
некоторых величин, называемых начальным состоянием системы. Рассмотрим понятие со-
стояния более подробно.
Понятие состояния системы (звена) является одним из базовых понятий теории динамических систем, поэтому оно определяется не через другие понятия, а аксиоматически
– перечислением совокупности присущих ему свойств [8, 20, 21]. Рассмотрим некоторые из них.
Как отмечено выше, выход динамической системы определяется однозначно, если заданы предыстория изменения входного процесса на некотором промежутке и, кроме того, некоторая совокупность величин, относящаяся к началу данного промежутка – начальное
1В дальнейшем термины система, подсистема и звено обычно будут использоваться как синонимы, так как их математические модели однотипны.
2Статической характеристикой в общем случае называют зависимость между входом и выходом системы в установившемся режиме (по истечении времени переходных процессов). Можно сказать, что у
безынерционных систем (звеньев) этот режим наступает немедленно.
7
состояние системы. Символически это будем записывать так: 3
y(t1) = S(x(t0); u[t0,t1]).
Таким образом, состояние системы – это некоторый параметр, позволяющий сделать однозначным определение ее выхода по входу.
Различные начальные состояния приводят, вообще говоря, к различной реакции на одно и то же входное воздействие. В приведенном выше уравнении S – некоторый оператор, преобразующий одну функцию в другую. 4
Состояние системы должно удовлетворять четырем аксиомам (условиям) совместности [20]. Рассмотрим две наиболее важные из них.
Аксиома 1. Выход y(t) для всех t ≥ t0 определяется однозначно, если заданы x(t0) и u[t0,t1] (см. рис. 1.1, а).
Рисунок 1.1 – Аксиомы совместности.
Таким образом, состояние системы в данный момент времени содержит всю память о прошлом, существенную для развития процесса в будущем. Если фиксировать начальное состояние, то будущее от прошлого не зависит; все, что нужно знать от прошлого для определения процесса в будущем, содержится в состоянии на данный момент времени. Таким образом, для определения будущего поведения системы не имеет значения то, как
3Через u[t0,t1] обозначено сужение функции u(·) на промежуток [t0, t1].
4Такая запись похожа на описание динамики систем передаточными функциями. Разница состоит в том, что передаточные функции используются для описания только линейных систем и позволяют определить реакцию лишь для нулевых начальных условий.
8
она пришла в данное состояние, – по начальному состоянию и входу процесс определяется однозначно.
Аксиома 2. Если траекторию системы разбить на ряд участков, то можно рассматривать движение на каждом из них как новую траекторию при соответствующем начальном
состоянии (см. рис. 1.1, б).
Пусть t0 < t1 < t2. Тогда y(t2) = S(x(t0); u[t0,t2]). С другой стороны, при любых
x(t0), u[t0,t1] можно определить состояние x(t1) таким образом, что y(t2) = S(x(t1); u[t1,t2]).
Из этой аксиомы следует, что состояние динамической системы должно изменяться во времени соответствующим образом (в зависимости от входного процесса и начального
состояния).
Определение. Множество X = {x} возможных значений состояния системы называ-
ется пространством состояний (данной системы). 5
Часто можно рассматривать в качестве пространства состояний n-мерное линейное вещественное пространство, X = Rn. Тогда состояние x(t) есть n-мерный вещественный вектор – вектор состояния, или фазовый вектор. Компоненты этого вектора обычно будем обозначать через xi(t), или, если возможны совпадения в обозначениях – через x(i)(t), т.е.
писать |
|
|
|
|
|
|
|
|
|
|
|
x1(t) |
|
|
|
|
x(1)(t) |
|
|
|
.. |
|
|
.. |
|
|
|||
x(t) = |
x2(t) |
|
или |
x(t) = |
x(2) |
(t) |
. |
||
|
. |
|
|
. |
|
||||
|
|
|
|
|
|
|
|
|
|
|
xn(t) |
|
|
x(n)(t) |
|
||||
Для краткости будем также использовать запись |
|
|
|
|
|||||
x = col{x1, x2, . . . , xn}, |
или |
x = [x1, x2, . . . , xn]T . |
Такая запись в общем случае означает, что x есть вектор-столбец, составленный из расположенных в столбец компонент векторов xi, i = 1, 2, . . . , n . Иногда, как будет видно из контекста, индекс используется для обозначения различных одноименных векторов.
Заметим, что такой вид пространства состояний не исчерпывает всех возможных ситуаций. Например, пространство состояний конечных автоматов состоит из конечного числа точек. С другой стороны, для многих систем нельзя указать конечное значение n размерности пространства X . К таким системам относятся различные распределенные объекты, динамика которых описывается дифференциальными уравнениями в частных производных,
5 Используется также термин «фазовое пространство».
9
объекты с запаздыванием и так далее. В этой книге рассматриваются только конечномерные динамические системы. Однако и для конечномерных систем не обязательно X = Rn.
Например, для простейшей механической системы – маятника – одной из переменных состояния является угол поворота относительно точки подвеса. Но в множестве возможных значений угловой переменной точки 0 рад. и 2π рад. совпадают. Следовательно, это множество не может быть линейным пространством, его геометрическим образом является не прямая, а окружность. Строгое рассмотрение таких систем требует привлечения понятия многообразия и выходит за рамки этой книги. Тем не менее многие свойства систем с угловыми координатами можно изучать, не используя аксиом линейного пространства. Поэтому, если не оговорено противное, мы будем считать, что X = Rn.
Из определения понятия состояния следует, что если x – состояние системы, μ(·) – некоторое взаимно однозначное отображение пространства X в себя (μ : X −→ X ), то x˜ = μ(x)
также можно рассматривать как состояние данной системы [8, 20]. Таким образом, состояние определяется неединственным образом, а с точностью до взаимно однозначного преобразования (которых может быть сколь угодно много). В частности, если X = Rn, а T –
некоторая невырожденная матрица порядка n (det T = 0), то вектор x˜ = T x также может
быть использован для описания состояния системы. Такой переход называется преобразованием базиса в пространстве состояний. Это преобразование не нарушает входо-выходных соотношений в описании системы.
Конкретизируем вид уравнений состояния. Рассмотрим так называемые конечномерные дифференциальные (непрерывные) системы. Уравнения состояния таких систем могут быть
представлены в виде |
|
|
|
|
|
|
x˙(t) = f x(t), u(t), t , x(t |
0) = |
x , t |
≥ |
t , |
(1.1) |
|
y(t) = g |
x(t), u(t), t . |
0 |
0 |
|||
|
|
|
|
|
|
|
Первое из этих уравнений – (собственно) уравнение состояния, или эволюционное уравнение, описывает изменение состояния системы во времени t R в зависимости от начальных условий в момент t0 и входного воздействия u(t). Второе уравнение – уравнение выхода, устанавливает связь между текущими значениями состояния и входа, с одной стороны, и
выхода y(t) – с другой. Фактически вся динамика системы сосредоточена в первом уравнении, а второе является статическим соотношением.
Переменные, входящие в уравнения (1.1), считаются векторными: x(t) Rn, y(t)
Rl, u(t) Rm, f(·), g(·) – вектор-функции от векторных аргументов соответствующих размерностей.
10