vostrikov
.pdfГ л а в а 2
ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ
Прежде чем изучать поведение реальных систем и их моделей, необходимо определить формальный язык, на котором
будут обсуждаться их свойства. Основным элементом такого формального языка является динамическая характеристика, под которой интуитивно понимают какое-либо соотношение, характеризующее свойства систем в статике и динамике (при изменении состояния).
Дадим следующее определение. Динамической характеристикой
(математической моделью) системы будем называть любое соотношение, заданное аналитически, графически или в виде таблицы, которое позволяет оценить ее поведение во времени.
В этой главе мы будем рассматривать различные способы описания линейных динамических систем, их взаимосвязь и приведение к принятой в теории автоматического управления форме записи математической модели.
Отметим, что динамическая характеристика дает возможность исследовать поведение системы, в частности, рассчитать для нее переходные процессы.
2.1. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ
Наиболее часто в качестве математической модели объекта управления используются обыкновенные дифференциальные уравнения, которые могут быть записаны в различной форме.
Линейные многоканальные объекты обычно описывают системой дифференциальных уравнений первого порядка, представленной в век- торно-матричном виде:
22 |
Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ |
||
|
x |
Ax Bu . |
(2.1) |
Здесь x |
Rn – вектор состояния, n – порядок объекта; u |
Rm – вектор |
|
управляющих воздействий, m |
n ; A – квадратная матрица действи- |
тельных коэффициентов; B – прямоугольная матрица действительных коэффициентов. Уравнения (2.1) называют дифференциальными уравнениями состояния.
Выходные переменные объекта изменяются в соответствии с урав-
нением выхода |
|
y Cx, |
(2.2) |
где y Rm – вектор выхода; C – прямоугольная матрица действитель-
ных коэффициентов. Уравнения (2.1) и (2.2) описывают линейный многоканальный объект.
Для описания одноканального объекта обычно используется скалярное дифференциальное уравнение
y(n) a y(n 1) |
|
a y |
a y bu, |
(2.3) |
n |
|
2 |
1 |
|
которое также может быть приведено к виду (2.1) и (2.2) после соответствующего выбора линейно-независимых переменных состояния.
Их число всегда равно порядку объекта (n), а u R1 и y R1 .
Наиболее простое каноническое описание получается, когда в качестве переменных состояния выбираются выходная переменная y и ее производные до (n 1) включительно:
x y, x |
y, , x y(n 1) . |
|
1 |
2 |
n |
При этом вместо (2.3) имеем систему уравнений |
||
x1 |
x2 , |
|
x2 |
x3 , |
|
|
|
(2.4) |
xn |
a1x1 |
a2 x2 an xn bu , |
y |
x1, |
|
2.1. Дифференциальные уравнения |
23 |
которая соответствует векторно-матричным уравнениям (2.1) и (2.2). Здесь матрицы A, B и C имеют вид
|
0 |
1 |
|
0 |
|
0 |
|
|
A |
0 |
0 |
0 |
, B |
0 |
, |
C 1 0 0 , |
|
|
|
|
|
|
|
|
|
|
|
a1 |
a2 |
|
an |
|
b |
|
|
причем |
их |
размерности |
следующие: |
dim A |
n n , dim B n 1 , |
|||
dim C |
1 n. |
|
|
|
|
|
|
|
Следует отметить, что переход к описанию (2.1), (2.2) не является однозначным: для одного объекта можно выбрать бесконечное множество наборов переменных состояния; важно, чтобы они были линейнонезависимыми. При этом каждой совокупности переменных состояния будут соответствовать свои матрицы объекта A, B и C.
ПРИМЕР 2.1
Записать уравнения состояния одноканального объекта, модель которого имеет вид
y 3y y u .
Рассмотрим два варианта переменных состояния.
1. Если в качестве переменных состояния использовать выходную величину и ее производную x1 y, x2 y , то получим канонические уравнения состояния и матрицы объекта типа (2.4):
x1 |
x2 , |
|
0 |
1 |
|
0 |
|
x |
x 3x u, A |
, B |
, C 1 0 . |
||||
2 |
1 |
2 |
1 |
3 |
|
1 |
|
y |
x1, |
|
|
|
|||
|
|
|
|
|
|
2. Выбирая новые переменные x1 y, x2 y 3y , получим уравнения состояния и матрицы объекта:
x1 |
3x1 x2 , |
|
3 |
1 |
|
0 |
|
x |
x u, |
A |
, B |
, C 1 0 . |
|||
2 |
1 |
|
1 |
0 |
|
1 |
|
y |
x1, |
|
|
|
|||
|
|
|
|
|
|
24 |
Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ |
В общем случае одноканальный объект может описываться дифференциальным уравнением вида
y(n) a y(n 1) |
a y |
b u(m) |
b u, n m . (2.5) |
n |
1 |
m |
0 |
Выбрав соответствующие переменные состояния, от описания (2.5) также можно перейти к векторно-матричным уравнениям типа (2.1), (2.2). Рассмотрим этот переход на примере.
ПРИМЕР 2.2
Записать уравнения состояния объекта с математической моделью
|
|
y |
y |
3y 2u u. |
|
|
Перепишем это уравнение: |
|
|
|
|
||
|
|
y |
2u |
y 3y u , |
|
|
выберем в качестве переменных состояния x1 y, |
x2 y 2u и получим |
|||||
следующие уравнения состояния и матрицы объекта: |
||||||
x1 |
x2 |
2u, |
|
0 |
1 , B |
2 , C 1 0 . |
x |
3x x u, |
A |
||||
2 |
1 |
2 |
|
3 |
1 |
1 |
y |
x1, |
|
|
|||
|
|
|
|
|
Таким образом, в качестве основной динамической характеристики линейных объектов управления используются дифференциальные уравнения, которые могут быть представлены в форме (2.1), (2.2).
2.2. СОСТАВЛЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
Поскольку в теории автоматического управления рассматриваются не реальные системы управления, а их математические модели, необходимо стремиться к тому, чтобы эти модели достаточно адекватно отражали свойства физических устройств. Процедуру получения математической модели объекта можно разбить на следующие этапы.
• Составление гносеологической (мысленной) модели объекта.
Исходя из технического задания и изучения режимов работы объекта инженер представляет себе приближенную модель, которая в дальнейшем уточняется и приобретает вид математической модели.
2.2. Составление математической модели |
25 |
• Определение независимых переменных, которые характеризуют объект, и уточнение их размерностей. При этом число управляющих воздействий не может быть меньше числа выходных переменных (dim u dim y) . Размерность вектора переменных состояния не может
быть |
меньше размерности вектора выходных переменных |
(dim x |
dim y) . Размерность возмущающих воздействий M может быть |
произвольной и никак не связана с размерностью y, x, u.
•Запись физических законов, в силу которых развиваются процессы в объекте.
•Приведение уравнений объекта к удобному с точки зрения тео-
рии автоматического управления виду.
Математическая модель никогда не бывает тождественна рассматриваемому объекту, так как при ее составлении всегда делают какиелибо допущения и упрощения. Поэтому для одной и той же системы в зависимости от целей управления модели могут быть различными.
При составлении математической модели приходится искать компромиссный вариант между двумя противоречивыми требованиями: с одной стороны, модель должна наиболее полно отражать свойства реальной системы, с другой – быть достаточно простой, чтобы не затруднять исследований.
ПРИМЕР 2.3
Определить математическую модель электрической цепи (рис. 2.1), записать для нее уравнения состояния.
Физическими законами, в соответствии с которыми развиваются процессы в объекте, являются законы Кирхгофа
|
U L |
dI |
RI , U |
|
|
RI. |
|
|
|
|
||
|
|
2 |
|
|
|
|
||||||
|
1 |
dt |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
L |
|||||
Uвых |
U2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
Rи |
0 |
|
R U |
|
|
R U R |
||||||
Rн |
|
|
и |
|
|
1 |
|
|
|
2 н |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 2.1. Электрическая схема объекта
26 |
Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ |
Перейдем к удобному с точки зрения теории управления описанию объекта. При этом выходной величиной будем считать напряжение на выходе цепи, т. е. y U2 , управляющим воздействием – напряжение на ее
входе u U1 , а переменной состояния – ток, протекающий по цепи x I . С учетом введенных обозначений запишем исходные уравнения
объекта в следующем виде:
Lx Rx u,
yRx,
азатем перейдем к принятому описанию в переменных состояния
|
|
|
|
|
|
x |
Ax |
Bu, |
|
|
|
||
|
|
|
|
|
|
y |
Cx, |
|
|
|
|
|
|
где A |
R L , B 1 L , |
C R . |
|
|
|
|
|
|
|||||
ПРИМЕР 2.4 |
|
|
|
|
|
|
|
||||||
|
|
I |
|
|
|
Рассмотрим |
в |
качестве еще одного |
|||||
|
|
|
|
|
примера |
составление |
математической |
||||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
модели |
двигателя |
постоянного тока с |
||||
U |
E |
ОВД |
|||||||||||
независимым возбуждением (рис. 2.2), |
|||||||||||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
который часто используется в системах |
||||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
автоматического управления. Здесь Uя – |
||||||
|
|
|
|
|
|
|
напряжение, подаваемое на якорь двига- |
||||||
|
Рис. 2.2. Схема двигателя |
|
теля, |
которое |
будем |
считать входным |
|||||||
Рис. 2.2. Схема двигателя посто- |
|||||||||||||
|
|
|
|
постоянного тока |
|
воздействием; |
I |
– ток в цепи якоря, |
|||||
|
|
|
|
|
янного тока |
|
|||||||
|
|
|
|
|
|
|
представляющий собой внутреннюю пе- |
||||||
ременную объекта; R, L |
– сопротивление и индуктивность цепи якоря; |
||||||||||||
E – противоЭДС, т. е. напряжение, |
возникающее в обмотке якоря в ре- |
||||||||||||
зультате его вращения в магнитном поле; |
– скорость вращения двигате- |
ля, которую будем считать выходной переменной; ОВД – обмотка возбуждения двигателя.
Запишем основные уравнения, характеризующие процессы в двигателе. Уравнение электрического равновесия якорной цепи имеет вид
dI
L dt RI E U я .
2.2. Составление математической модели |
|
|
|
|
|
|
|
|
|
|
|
|
|
27 |
||||||||||||||||||
|
Уравнение равновесия моментов на валу двигателя |
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
J |
d |
|
|
|
M д |
M c , |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
где J |
– приведенный момент инерции; |
Mд – вращающий момент; |
M с – |
|||||||||||||||||||||||||||||
момент сопротивления на валу двигателя, который является возмущаю- |
||||||||||||||||||||||||||||||||
щим воздействием. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
С достаточной степенью точности во многих случаях можно считать, |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
что |
E |
c1 |
, |
Mд |
|
c2I, |
Mc |
Mc (t), где |
ci |
|
|
const, i 1, 2. В |
результате |
|||||||||||||||||||
уравнения двигателя принимают вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
L |
|
|
dI |
|
|
|
RI |
c |
|
|
U |
|
|
, |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
я |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
J |
|
d |
|
|
|
c I |
M |
|
. |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Введем |
следующие обозначения: |
u |
|
Uя |
|
– управление; |
x1 |
и |
|||||||||||||||||||||||
x2 |
I |
– переменные состояния; |
M с |
– возмущение. Запишем уравнения |
||||||||||||||||||||||||||||
двигателя в переменных состояния: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
x1 |
|
|
|
a12 x2 |
hM c , |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
x2 |
a21x1 |
a22 x2 |
|
bu, |
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
y |
|
|
|
x1, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
где |
a |
|
c2 |
, |
h |
1 |
, |
a |
|
c1 |
, |
a |
|
R |
|
, b |
|
1 |
. |
|
|
|
|
|||||||||
|
12 |
|
J |
|
|
J |
21 |
|
|
L |
22 |
L |
|
|
|
L |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
Часто модель двигателя представляют в виде одного дифференциаль- |
|||||||||||||||||||||||||||||||
ного уравнения |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
TяTм y |
Tм y |
|
y ku kм (Tя p 1)Mc . |
|
|
|||||||||||||||||||||
Здесь Tм |
JR c1c2 |
– электромеханическая постоянная времени двигате- |
||||||||||||||||||||||||||||||
ля; |
Tя |
|
L R |
– |
электромагнитная постоянная времени якорной |
цепи; |
||||||||||||||||||||||||||
k |
1 c1 – коэффициент усиления; |
kм |
|
|
R c1c2 . |
|
|
ПРИМЕР 2.5
Рассмотрим перевернутый маятник, ось которого монтируется на тележке (каретке), перемещающейся в горизонтальном направлении [21].
28 |
Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ |
В совокупности такое устройство представляет собой объект управления, называемый «каретка–маятник». Его схематичная модель показана на рис. 2.3, где использованы следующие определения:
– угол отклонения маятника (выходная переменная); U – прикла-
дываемая управляющим двигателем сила (входная переменная); s – перемещение каретки;
m1 – масса каретки;
L – расстояние между осью и центром тяжести маятника; m2 – масса маятника;
J – момент инерции относительно центра тяжести; g – ускорение силы тяжести;
H и V – горизонтальная и вертикальная силы реакции у оси маятника.
|
|
Y |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
φ |
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
Центр тяжести |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
m2g |
|
|
U |
|
|
|
|
|
|
H |
|
|
|
|
|
|
||
m1 |
|
Ось |
|
||||
|
|
|
|
||||
|
m1 |
|
|
|
|
|
|
0
s
Рис. 2.3. Объект управления «каретка–маятник»
Упрощенная модель объекта «каретка–маятник» может быть представлена системой дифференциальных уравнений [9]
|
|
|
|
|
|
a4c |
cs |
0, |
|
|
|
|
|
|
|
s |
a2 s b2U , |
|
|||
где |
a |
F |
, b |
1 |
, a g, |
c 1 |
|
J |
m1L2 |
– эффективная длина ма- |
|
|
|
|
|
||||||
|
2 |
2 |
4 |
|
|
|
m1L |
|
||
|
|
m 2 |
m 2 |
|
|
|
|
ятника.
При переходе к описанию модели объекта в переменных состояния в качестве компонент вектора состояния можно выбрать следующие величины:
2.3. Переходная характеристика |
29 |
x |
s, x |
s, |
x |
|
s c 1 |
, |
1 |
2 |
|
3 |
|
|
|
|
x (t) |
s(t) |
c |
1 (t), |
|
|
|
4 |
|
|
|
|
|
а выходной переменной |
объекта является угол |
отклонения маятника |
y . В результате уравнения состояния принимают вид
x1 x2 ,
x2 a2 x2 b2U , x3 x4 ,
x4 a4c(x3 x1), y c(x3 x1).
Определив матрицы |
|
|
|
|
|
|
|
|
0 |
1 |
0 |
0 |
|
0 |
|
A |
0 |
a2 |
0 |
0 |
, B |
b2 |
, |
|
0 |
0 |
0 |
1 |
|
0 |
|
|
a4c |
0 |
a4c |
0 |
|
0 |
|
|
C |
c |
0 |
c |
0 , |
|
|
модель объекта «каретка–маятник» можно представить в векторноматричной форме (2.1)–(2.2).
2.3. ПЕРЕХОДНАЯ ХАРАКТЕРИСТИКА
Эта динамическая характеристика используется для описания одноканальных объектов
y(n) a y(n 1) |
a |
y(n 2) a y b u(m) |
b u |
|
n |
n 1 |
1 |
m |
0 |
с нулевыми начальными условиями
y(0) 0, y(0) |
0, , y(n 1) (0) 0. |
30 Глава 2. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ЛИНЕЙНЫХ СИСТЕМ
Переходной характеристикой (переходной функцией) h(t) назы-
вается реакция системы на единичное ступенчатое входное воздействие u(t ) 1(t ) при нулевых начальных условиях (см. рис. 2.4).
Отметим, что единичная ступенчатая функция – это функция, которая обладает свойством
1(t |
) |
0, t |
0, |
|
1, t |
0. |
|||
|
|
Здесь – момент возникновения входного воздействия.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
u |
|
u |
|
|
|
|
||||
|
|
|
|
|
|
1(t |
) |
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
|
||||||||||
|
|
|
|
|
|
|
1(t |
) |
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
t |
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
h(t |
) |
) |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
h(t |
|
|
|
tt
Рис. 2.4. Пример переходной характеристики системы
Для аналитического определения переходной функции следует решить дифференциальное уравнение при нулевых начальных условиях и единичном входном воздействии.
При исследовании реального объекта переходную характеристику можно получить экспериментальным путем, подавая на его вход ступенчатое воздействие и фиксируя реакцию на выходе. Если входное воздействие представляет собой неединичную ступенчатую функцию u(t) k1(t) , то выходная величина будет равна y(t) kh(t) , т. е. пере-
ходной характеристике с коэффициентом пропорциональности k.
Зная переходную характеристику, можно вычислить реакцию системы на произвольное входное воздействие с помощью интеграла свертки
t |
|
y(t) h(t)u(0) h(t )u( )d , |
(2.6) |
0 |
|
где – переменная интегрирования.