Файл: Казанский государственный университет.pdf

служить переменными состояния этой цепи.
Функции
φ
и
ϕ являются алгебраическими даже в случае нелинейной и нестационарной цепи, так как уравнения составлялись с использованием мгновенных значений токов и напряжений.
Для иллюстрации рассмотрим пример описания электрических цепей в пространстве состояний.
Рассмотрим пассивный многополюсник (рис.5.3.1) со стационарными параметрами, который представляет собой систему 2-го порядка.
Состояние рассматриваемой цепи определяется энергией, запасенной в реактивных элементах и
C
L
, и входным сигналом
U .
Энергия, запасенная в индуктивности, определяется током в индуктивности
, а энергия, запасенная в конденсаторе,
– напряжением
:
( )
L
i t
( )
c
U t
2
( )
( )
2
L
L
Li t
E t
=
;
( )
( )
2
c
c
CU t
E t
=
Выберем эти переменные в качестве переменных состояния
Рис. 5.3.1. Схема многополюсника
1 2
L
c
x ( t ) i ( t ),
x ( t ) U ( t )
=
=
Уравнения Кирхгофа для данной цепи будет иметь вид:
1 1
2
R
L
c
R
L
c
R
U
U
U
U ;
i i
i
i
i
+
+
=
=
= = +
.
Учитывая
c
c
dU
i
C
и
dt
=
L
L
di
U
L
, получим соотношения
dt
=
1
L
L
di
c
R i
L
U
U
dt
+
+
= и
2 1
c
L
c
dU
i
C
U
dt
R
=
+
После преобразований уравнения состояния окончательно примут вид
1 2
1 1
,
1 1
L
L
c
c
L
c
di
R
i
U
dt
L
L
L
dU
i
U
dt
C
R C
⎧
= −
−
+
⎪⎪
⎨
⎪
=
−
⎪⎩
U
;
В соответствии с заданной схемой (рис. 5.3.1) выходные сигналы можно описать следующими соотношениями:
1 1
2
L
c
c
y
U
R i
U
U
y
U .
= = −
−
+
=
Из последних соотношений после замены переменных легко прийти к стандартной записи уравнений в пространстве состояний – заданная
139

электрическая цепь описывается матричными уравнениями вида
d
X ( t ) AX ( t ) BU( t )
dt
=
+
,
Y( t ) CX ( t ) DU( t )
=
+
, где
1 1
2 1
1 1
1 1
0 0
1
R / L
/ L
/ L
R
A
, B
, C
,
/ C
/ R C
−
−
−
−
⎡
⎤
1 0
D
⎡
⎤
⎡
⎤
⎡
=
=
=
⎢
⎥
⎤
=
⎢
⎥
⎢
⎥
⎢
−
−
⎥
⎣
⎦
⎣
⎦
⎣
⎣
⎦
⎦
,
[
]
[
]
[
]
[
]
1 2
1 2
T
T
T
T
L
c
L
c
X(t )
x (t ) x (t )
i (t ) U (t ) ,
Y(t )
y (t ) y (t )
U (t ) U (t )
=
=
=
=
Таким образом, получена математическая модель 2-го порядка электрического многополюсника в пространстве состояний.
Отметим ту особенность, что, как правило, порядок математической модели электрических систем совпадает с числом реактивных элементов
(числом электрических емкостей и индуктивностей) в системе.
5.3.3. Особенности составления уравнений состояния для
электромеханических систем
Особенности составления уравнений состояний для электромеха- нических систем проистекают из того факта, что в таких системах существуют и взаимодействуют как механические элементы, так и электрические цепи. Отсюда очевидно то обстоятельство, что необходим учет особенностей как той, так и другой частей электромеханической системы.
Для примера рассмотрим составление описания двигателя постоянного тока в пространстве состояний.
Процессы в электродвигателе постоянного тока характеризуются следующими соотношениями: уравнение электрической цепи
я
d
я
я
я
di
U e
i R
L
dt
=
+
+
; уравнение механической цепи
d
c
d
M
M
J
dt
ω
=
+
; уравнение момента сопротивления
c
M
b
ω
=
; уравнение момента двигателя
d
я
M
KФi
=
; уравнение Э.Д.С. двигателя
d
e
KФ
ω
=
140

Здесь в этих уравнениях:
я
i – ток якоря двигателя;
я
R ,
я
L – активное сопротивление (Ом) и индуктивность (Гн) якорной цепи;
U
– напряжение, В; – Э.Д.С. двигателя, В;
d
e
d
M – действующий момент двигателя, Нм;
c
M – момент сопротивления на валу двигателя, Нм:
J
– момент инерции, приведенный к валу двигателя;
ω
– угловая скорость вращения вала двигателя, рад/сек;
– магнитный поток возбуждения двигателя, Вб;
Ф
2
pN
K
a
π
=
– конст- руктивный коэффициент двигателя (p, N, a – соответственно число пар полюсов, активных проводников и параллельных ветвей обмотки якоря).
Обычно,
KФ C const
= =
Выберем в качестве переменных состояния
1
я
x
i
=
(ток якоря электродвигателя – характеризует электрическую часть системы) и
2
x
ω
=
(угловая скорость вращения механической части системы). Связи между электрической и механической частями системы описываются уравнениями для момента сопротивления и момента двигателя.
Преобразуя приведенные выше уравнения – приводя дифференциальные уравнения к нормальной форме Коши и осуществляя необходимые замены переменных – получим уравнения состояния этого электродвигателя в следующем виде:
/
/
1/
/
/
0
я
я
я
я
я
я
i
R L
C L
i
L
d
U
C J
b J
dt
ω
ω
−
−
⎡ ⎤ ⎡
⎤ ⎡ ⎤ ⎡
⎤
=
+
⎢ ⎥ ⎢
⎥ ⎢ ⎥ ⎢
⎥
−
⎣ ⎦ ⎣
⎦ ⎣ ⎦ ⎣
⎦
Вводя обозначения
1 1
я
я
я
R / L
/ T ,
b / J
/ T
м
−
= −
−
= −
, которые в электротехнике имеют четко выраженный физический смысл, получаем уравнения состояния с матрицами
1/
1/
1/
,
/
1/
0
я
я
я
м
T
L
L
А
B
C J
T
−
−
⎡
⎤
⎡
⎤
=
=
⎢
⎥
⎢
⎥
−
⎣
⎦
⎣
⎦
Матрицы и будут определяться исходя из условий наблюдения.
Так, если выходной величиной будет угловая скорость
C
D
ω
, то матрица примет вид
C
[
0 1
C
]
=
, если же в качестве выходной величины будет рассматриваться действующий момент двигателя
d
M
, то
[
]
0
C
c
=
141

5.4. Формирование уравнений состояний по
дифференциальному уравнению
Наиболее простым методом описания динамических систем является поэлементное описание системы, когда каждый элемент системы описывается своим дифференциальным (реже алгебраическим) уравнением, а взаимосвязи элементов в системе описываются уравнениями связи. Сами же дифференциальные уравнения элементов формируются на основе тех фундаментальных физических законов, на действии которых базируется работа элемента динамической системы. В общем случае дифференциальное уравнение элемента системы и, тем более, дифференциальное уравнение самой системы представляют собой уравнение высокого порядка.
Поставим задачу формирования уравнений состояния системы в виде уравнений (5.2.6) и (5.2.7) по известному дифференциальному уравнению
n-го порядка.
Рассмотрим дифференциальное уравнение с постоянными коэффи- циентами для динамической системы с одним входом и одним выходом
1 1
1 0
( n )
( n
)
n
n
a y ( t ) a y
( t ) ... a y ( t ) a y( t ) b U( t )
−
−
′′
+
+ +
+
=
0
, (5.4.1) где
(
)
1
i
a
i
,...n
=
и – постоянные коэффициенты.
0
b
Введем систему следующих допущений и обозначений
1 2
1 3
2 1
1 0
1 2
3 1
0
n
n
( n )
( n
)
n
n
n
n
n
n
n
n
x
y
x
x
x
x
...
....
x
x
a
a
a
a
a
b
x
y ( t )
y
y
y
y
...
y
U( t )
a
a
a
a
a
a
−
−
−
=
′
=
′
=
′
=
′
′
′′
′′′
=
= −
−
−
−
−
+
.
(5.4.2)
Последнее уравнение получено из уравнения (5.4.1). Это уравнение можно переписать с учетом введенных ранее в (5.4.2) обозначений следующим образом
( )
0 1
2 3
1 0
1 2
3 4
( )
..._
( )
n
n
n
n
n
n
n
n
n
n
a
a
a
a
a
b
x
y
t
x
x
x
x
x
U t
a
a
a
a
a
a
−
′ =
= −
−
−
−
−
+
Тогда вся система соотношений (5.4.2) может быть представлена в следующем виде (здесь и далее для лучшей иллюстрации независимый аргумент при переменных величинах упущен):
t
142

1 1
2 2
0 1
1 2
0 1
0 0
0 0
0 1
0 0
1
n
n
n
n
n
n
n
...
x
x
...
x
x
d
U( t )
...
...
...
...
...
...
...
...
dt
a
a
a
a
x
...
x
a
a
a
a
−
⎡
⎤
⎡ ⎤
⎡ ⎤ ⎡ ⎤
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
=
+
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
−
−
−
−
⎢
⎥
⎣ ⎦
⎣ ⎦
⎣ ⎦
⎣
⎦
;
(5.4.3)
1 2
0 0
0
n
n
x
x
b
Y
...
...
a
x
⎡ ⎤
⎢ ⎥
⎡
⎤ ⎢ ⎥
= ⎢
⎥ ⎢ ⎥
⎣
⎦
⎢ ⎥
⎣ ⎦
(5.4.4)
Таким образом, дифференциальному уравнению (5.4.1) в пространстве состояний соответствуют уравнения
d
X
AX
BU
dt
Y CX
DU ,
=
+
=
+
,
(5.4.5) где
[
]
1 2
T
n
X
x
x
... x
=
,
[ ]
1
Y
y
=
,
[ ]
0 1
1 2
0 0
1 0
0 0
0 0
1 0
0 1
0 0
n
n
n
n
n
n
...
...
A
,
B
...
...
...
...
...
...
a
a
a
a
...
a
a
a
a
b
C
... . ,
D
.
a
−
⎡
⎤
⎡ ⎤
⎢
⎥
⎢ ⎥
⎢
⎥
⎢ ⎥
=
=
⎢
⎥
⎢ ⎥
⎢
⎥
⎢ ⎥
−
−
−
−
⎢
⎥
⎣ ⎦
⎣
⎦
⎡
⎤
=
=
⎢
⎥
⎣
⎦
,
Граф, соответствующий системе уравнений (5.4.2) или уравнениям
(5.4.3) и (5.4.4) при условии
1
n
a
=
, представлен на рис. 5.4.1.
Рис. 5.4.1. Граф системы (5.4.2)
Структурная схема, отображающая ту же систему уравнений, имеет вид,
143

показанный на рис. 5.4.2, где использованы обозначения:
0 0
i
i
n
n
a
b
,
a
a
α
β
−
=
=
Рис. 5.4.2. Структурная схема системы
Полученная форма уравнений состояния (5.4.3) и (5.4.4) носит название
строчной управляемой фробениусовой канонической формы
и является одной из самых распространенных на практике форм. Другие канонические формы рассматриваются в разделе 5.9. Здесь же отметим, что основу структурной реализации этих форм составляет цепочка последовательно включенных интеграторов, охваченных обратными связями, причем коэффициенты обратных связей совпадают с коэффициентами исходного дифференциального уравнения.
Рассмотрим более сложный случай, когда исходное дифференциальное уравнение содержит в правой части не только управляющее воздействие, но и его производные (как и ранее, полагаем, что в динамической системе есть только один вход и один выход), т.е. имеет вид
1 1
1 0
1 1
1
( n )
( n )
n
n
( m )
( m )
m
m
a y (t ) a y
(t ) ... a y (t ) a y(t )
b u (t ) b u
(t ) .... bu (t ) b u(t ).
−
−
−
−
′′
+
+ +
+
=
′
=
+
+ +
+
0
(5.4.6)
(Здесь и далее будем считать
).
m n
<
Для решения задачи используем вспомогательную переменную
Z
Пусть
1 1
1 0
( n )
( n
)
n
n
a z ( t ) a z
( t ) ... a z ( t ) a z( t ) U( t )
−
−
′
+
+ +
+
=
, (5.4.7) и
1 1
1 0
( m )
( m )
m
m
b z ( t ) b z
( t ) .... b z ( t ) b z( t ) Y( t )
−
−
′
+
+ +
+
=
. (5.4.8)
Корректность использования соотношений (5.4.7) и (5.4.8) легко проверить, если подставить эти выражения для
U и в уравнение (5.4.6), которое после подстановки обращается в тождество.
Y
144

Далее из уравнения (5.4.7) как и ранее, положив
1
x
z
=
, получим систему уравнений
1 2
2 3
1 0
1 1
1 2
1
n
n
n
n
n
n
n
n
n
x
x
x
x
x
x
a
a
a
x
x
x
x
a
a
a
a
−
−
′ =
′ =
′ =
′ = −
−
− −
+
U
Использовано именно уравнение (5.4.7) по той простой причине, что это дифференциальное уравнение более высокого порядка, чем уравнение
(5.4.8) в силу
. Тогда систему уравнений можно представить в векторно-матричной форме
m n
<
1 1
2 2
0 1
1 2
0 1
0 0
0 0
0 1
0 0
1
n
n
n
n
n
n
n
...
x
x
...
x
x
d
U( t )
...
...
...
...
...
...
...
...
dt
a
a
a
a
x
...
x
a
a
a
a
−
⎡
⎤
⎡ ⎤
⎡ ⎤ ⎡ ⎤
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
=
+
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
⎢
⎥
⎢ ⎥
⎢ ⎥ ⎢ ⎥
−
−
−
−
⎢
⎥
⎣ ⎦
⎣ ⎦
⎣ ⎦
⎣
⎦
, а уравнение (5.4.8) может быть представлено в следующем виде:
[ ]
1 2
0 1
0 0
0
m
n
n
n
n
x
x
b
b
b
Y
...
...
U( t )
...
a
a
a
x
⎡ ⎤
⎢ ⎥
⎡
⎤ ⎢ ⎥
=
+
⎢
⎥ ⎢ ⎥
⎣
⎦
⎢ ⎥
⎣ ⎦
Таким образом, для исходного дифференциального уравнения (5.4.6) n-го порядка матричные уравнения состояния имеют тот же стандартный вид уравнений (5.4.5), но состав матриц
A, B, C, D
уже другой.
Конкретно
[ ]
0 1
1 2
0 1
0 1
0 0
0 0
0 1
0 0
,
,
1 0 ... 0 ,
0
n
n
n
n
n
m
n
n
n
A
B
a
a
a
a
a
a
a
a
b
b
b
C
D
a
a
a
−
⎡
⎤
⎡ ⎤
⎢
⎥
⎢ ⎥
⎢
⎥
⎢ ⎥
=
=
⎢
⎥
⎢ ⎥
⎢
⎥
⎢ ⎥
−
−
−
−
⎢
⎥
⎣ ⎦
⎣
⎦
⎡
⎤
=
=
⎢
⎥
⎣
⎦
145

Полученная форма уравнений состояния также относится к строчным управляемым фробениусовым каноническим формам. Граф, соответству- ющий полученным уравнениям при условии
1
n
a
=
, представлен на рис. 5.4.3.
Рис. 5.4.3. Граф системы
Структурная схема, соответствующая полученным уравнениям, представлена на рис. 5.4.4. Использованы обозначения, введенные ранее.
Рис. 5.4.4. Структурная схема системы
Нетрудно заметить, что полученные графические отображения математической модели динамической системы обладают явно выраженной регулярностью, выражающейся в однотипности выполняемых преобразований и их представлении в виде графических изображений.
146