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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 07.11.2023

Просмотров: 301

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
5.5. Формирование уравнений состояний по передаточной
функции
Нахождение уравнений состояния динамической системы в пространстве состояний по известной передаточной функции системы является неоднозначной процедурой из-за произвола выбора переменных состояния.
Одним из наиболее распространенных является метод, использующий в качестве системы переменных состояния фазовые координаты.
Рассмотрим решение задачи формирования уравнений состояний для системы с одним входом и одним выходом и имеющей передаточную функцию вида
1 1
1 1
1 1
n
n
n
n
n
n
n
Y( p )
b p
b p
... b p b
W ( p )
U( p )
p
a p
... a p a




0 0
+
+ +
+
=
=
+
+ +
+
, (5.5.1) где Y(
и U(
– изображения по Лапласу выходного и входного сигналов системы соответственно.
p )
p )
Соответствующее этой передаточной функции дифференциальное уравнение динамической системы запишется в виде
1 1
1 0
1 1
1
( n )
( n )
n
( n )
( n )
n
n
y (t ) a y
(t ) ... a y (t ) a y(t )
b u (t ) b u
(t ) ... bu (t ) b u(t )





+
+ +
+
=

=
+
+ +
+
0
.
[ ]
(5.5.2)
Здесь
– выходной сигнал системы, u( – входной сигнал системы.
y( t )
t )
Из сравнения (5.5.2) и (5.4.6) видно, что они тождественны друг другу.
Поэтому дальнейший переход к уравнениям состояния полностью совпадает с методом, изложенным в предыдущем разделе.
В результате получим, что матрицы уравнений состояния динамической системы имеют вид
[
]
0 1
2 1
0 1
2 1
0 1
0 0
0 0
0 1
0 0
1
n
n
n
...
...
A
,
B
...
...
...
...
...
...
a
a
a
...
a
C
b
b
b
... b
,
D
b .




,
⎡ ⎤


⎢ ⎥


⎢ ⎥
=
=


⎢ ⎥


⎢ ⎥




⎣ ⎦


=
=
(5.5.3)
Граф математической модели динамической системы, описываемой передаточной функцией (5.5.1) или, что равнозначно, уравнениями состояния с матрицами (5.5.3), приведен на рис. 5.5.1.
147

Ранее уже отмечалось, что одной и той же передаточной функции можно поставить в соответствие уравнения состояния вида (5.2.6) и (5.2.7), но с различными матрицами
A
, , ,
Так, если в дифференциальном уравнении (5.5.2), соответствующем передаточной функции (5.5.1), пере- менные состояния определить равенствами,
B
Рис. 5.5.1. Граф системы
C
D
1 0
2 1
3 1
( 1)
(
2)
(
3)
1 2
1
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( ) ...
( ),
n
n
n
n
n
x t
y t
b U t
x t
y t
bU t
x t
y t
bU t
x t
y
t
bU
t
b U
t
b U t




=


=

′′
=

=


− −
то матрицы уравнений состояния динамической системы принимают вид:
[
]
[ ]
1 2
0 1
2 1
0 0
1 0
0 0
0 1
0 1 0 0 0
n
n
...
...
A
,
B
...
...
...
...
...
...
a
a
a
...
a
C
...
,
D
,
.
β
β
β
β







=
=










=
=
⎡ ⎤
⎢ ⎥
⎢ ⎥
⎢ ⎥
⎢ ⎥
⎣ ⎦
(5.5.4)
Элементы
(
)
1 2
i
i
, ,...,n
β
=
матрицы связаны с коэффициентами исходной передаточной функции следующими соотношениями:
B
0 1
1 1 0 2
2 1 1 2 0 3
3 1 2 2 1 3 0 0
1 1
2 2
3 3
0
,
,
,
,
n
n
n
n
n
n
n
n
n
n
n
n
n
n
n
n
n
b
b
a b
b
a b
a b
b
a b
a b
a b
b
a b
a b
a b
a b















β =
β =

β =


β =



β = −


− −
0
Граф соответствующий системы с матрицами вида (5.5.4) представлен на рис. 5.5.2.
Уравнения состояния динамической системы с такими матрицами A, B, C, D носят название строчной наблюдаемой фробениу- совой канонической формы.
148


Нетрудно увидеть, что хотя в (5.5.3) и (5.5.4) матрицы
A
систем одинаковы, но весьма сильно различаются матрицы ,
C , .
B
D
Другой способ формирования уравнений состояний динамической системы в пространстве состояний по известной передаточной функции системы заключается в разложении передаточной функции на элементарные дроби и в последующем представлении совокупности этих дробей в пространстве состояний.
W( p )
Рис. 5.5.2. Граф системы
Рассмотрим систему с одним входом и одним выходом, передаточная функция которой представляет собой отношение двух полиномов
W( p )
1 1
1 1
1 1
m
m
m
m
n
n
n
n
0 0
B( p ) b p
b p
... b p b
W( p )
A( p )
a p
a p
... a p a




+
+ +
+
=
=
+
+ +
+
. (5.5.5)
Полагаем, что
. Если предположить, что корни полинома
m n
<
A( p ) различные, то полином A( p ) можно представить в виде:
1 2
1
n
n
n
A( p ) a ( p
)( p
)...( p
)( p
)
λ
λ
λ

=




λ
)
, (5.5.6) где
(
1 2
i
i
, ,...,n
λ
=
- корни характеристического уравнения
0
A( p )
= .
В этом случае передаточная функция W(
раскладывается в виде
p )
1 2
1 0
1 2
1
n
n
c
c
c
c
W ( p )
c
...
p
p
p
p
1
λ
λ
λ

=
+
+
+
+
+




λ
, (5.5.7) где
0
i
i
B( p )
c
limW( p )
при p
, c
( p
)
при p
A( p )
i
λ
λ
=
→ ∞
=

= .
Постоянные коэффициенты и могут быть найдены также с помощью метода неопределенных коэффициентов.
0
c
i
c
Изображение выходного сигнала можно представить в виде
0 1
n
i
i
i
c
Y( p ) c U( p )
U( p )
p
λ
=
=
+


(5.5.8)
Соответствующий этому выражению граф модели динамической системы представлен на рис. 5.5.3.
149

Каждая из переменных
i
x
удовлетворяет дифференциальному уравнению 1 порядка:
i
i i
x ( t )
x ( t ) U( t )
λ

=
+
(5.5.9) или, иначе говоря, удовлетворяет частной передаточной функции
i
i
i
c
W ( p )
p
λ
=

0 1
n
i i
i
Y( t )
c x ( t ) c U( t )
=
=
+

Сигнал на выходе всей динами- ческой системы будет описываться формулой
Рис. 5.5.3. Граф системы
. (5.5.10)
Уравнения состояния, соответствующие исходной передаточной функции и графу системы на рис. 5.5.3, будут иметь вид
[
]
[ ]
1 1
1 2
2 2
1 2
1 2
0 0
0 1
0 0
1 1
0 0
1
n
n
n
n
n
x ( t )
...
x ( t )
x ( t )
...
x ( t )
d
U( t ),
...
... ... ... ...
...
dt
x ( t )
...
x ( t )
x ( t )
x ( t )
Y( t )
c
c
... c
c
U( t ).
...
x ( t )
λ
λ
λ



⎤ ⎡

⎡ ⎤



⎥ ⎢

⎢ ⎥



⎥ ⎢

⎢ ⎥
=
+



⎥ ⎢

⎢ ⎥



⎥ ⎢

⎢ ⎥
⎣ ⎦



⎦ ⎣







=
+






(5.5.11)
Матрица в этом случае носит название диагональной формы
Жордана. В случае, если корни многочлена знаменателя
i
λ
не являются различными, т.е. являются кратными, или в случае комплексных корней, метод разложения передаточной функции на элементарные дроби приводит к более сложной матрице A модели динамической системы.
150


5.6. Формирование уравнений состояний по структурной
схеме
Структурная схема (или граф) динамической системы дает условное графическое изображение взаимосвязанных элементов, составляющих систему и различающихся своими динамическими свойствами. При этом в графическом изображении системы указываются и преобразования, которым подвергается сигнал при прохождении через тот или иной элемент.
Формирование уравнений состояния системы по известной структурной схеме (графу) может осуществляться несколькими способами.
Среди них можно отметить:
1) получение передаточной функции системы и последующий переход от передаточных функций к уравнениям в пространстве состояний;
2) представление структурной схемы системы в виде совокупности типовых динамических звеньев I и II порядков, описание этих типовых звеньев в пространстве и последующее использование правил формирования уравнений состояния системы по известным уравнениям подсистем;
3) детализацию структурной схемы системы, с тем чтобы блоки схемы описывали простые одинарные операции (сложение, вычитание, умножение, дифференцирование, интегрирование и т.п.) и дальнейшую запись системы уравнений преобразования.
Для иллюстрации последнего способа рассмотрим пример, когда известна структурная схема некоторого динамического объекта
(рис.5.6.1)
Когда структура динамической системы детализована до простых операций сложения, вычитания, умножения, дифференцирования, интегрирования и т.д., то наиболее целесообразным является выбор в качестве переменных состояния сигналов на выходе блоков интегрирования.
Сами же уравнения состояний формируются далее по структурной схеме системы. рис. 5.6.1. Структурная схема и граф системы
151

В нашем конкретном случае по структурной схеме (или графу) системы можно записать следующие уравнения:
1 1
2 2
1 2
1
R
C
x
x
x
U ,
L
L
L
C
b
x
x
x .
J
J
′ = −

+
′ =

(5.6.1)
Отсюда легко установить, что в уравнении состояния объекта матрицы имеют вид
1
,
0
R
C
L
L
A
B
L
C
b
J
J




⎡ ⎤


⎢ ⎥
=
=


⎢ ⎥



⎣ ⎦




:
Нетрудно увидеть, что получился результат, полностью совпадающий с результатом примера в разделе 5.3.3. Это говорит о том, что исходная структурная схема соответствовала электродвигателю постоянного тока.
5.7. Формирование уравнений состояний системы по
известным уравнениям подсистем
Сложные динамические комплексы образуются совокупностью более простых динамических подсистем. Как правило, при этом используются типовые соединения: последовательные, параллельные и встречно- параллельные (соединения с обратными связями). Тогда правомерен вопрос об эквивалентном описании подсистем, соединенных между собой тем или иным видом соединения.
Полагаем, что каждая динамическая подсистема описана в пространстве состояний матричными уравнениями вида
1 2
i
i
i
i
i
i
i
i
i
d
X ( t ) A X ( t ) BU ( t ),
i
,
dt
Y( t ) C X ( t ) DU ( t ),
=
+
=
=
+
.
(5.7.1)
Принадлежность векторов и матриц к той или иной динамической подсистеме, участвующей в соединении, будем указывать индексами 1 и 2.
Соединения динамических подсистем обычно описываются с помощью метода расширения фазового пространства. В этом случае для описания объединенных подсистем вводится расширенный вектор состояния
152


1 2
X ( t )
X ( t )
X ( t )


= ⎢



,
(5.7.2) где
– вектор состояния первой системы,
– вектор состояния второй системы.
)
(
1
t
X
)
(
2
t
X
Объединенная динамическая система будет описывать уравнения такого же вида, что и ранее
d
X ( t ) AX ( t ) BU( t ),
dt
Y( t ) CX ( t ) DU( t ),
=
+
=
+
(5.7.3) но матрицы A
, , , ,
естественно, будут иметь совершенно другое содержание.
B
C
D
Рассмотрим различные типовые соединения подсистем друг с другом и получим соотношения, связывающие матрицы подсистем и системы.
1. Последовательное соединение двух динамических систем (рис. 5.7.1).
Для рассматриваемых подсистем справедливы уравнения: рис. 5.7.1. Последовательное соединение динамических подсистем
1 1
1 1 1 2
2 2
2 1
1 1
1 1 2
2 2
2 2
d
d
X
A X
BU ,
X
A X
B U
dt
dt
Y
C X
DU ,
Y
C X
D U .
=
+
=
+
=
+
=
+
2
,
(5.7.4)
При последовательном соединении подсистем выходной сигнал первой подсистемы является входным сигналом для второй подсистемы, т.е. U
Y
2 1
= , а выходной сигнал второй подсистемы является и выходным сигналом всей системы. Тогда уравнения системы можно представить следующим образом:
1 1
1 1 1 2
2 2
2 1 1
2 1
2 2
2 2 1 1
2 1 1
d
X
A X
BU ,
dt
d
X
A X
B C X
B DU
dt
Y
C X
D C X
D DU .
=
+
=
+
+
=
+
+
1
,
(5.7.5)
Вводя в рассмотрение расширенный вектор состояний
,
[
]
1 2
T
X
X
X
=
153
можно получить уравнения
[
]
[
]
1 1
1 1
2 2
1 2
2 1
2 1
2 2
1 2
0 0
X
A
X
B
d
U ,
X
B C
A
X
dt
X
Y
D C
C
D D U ,
X

⎤ ⎡
⎤ ⎡
⎤ ⎡ ⎤
=
+

⎥ ⎢
⎥ ⎢
⎥ ⎢ ⎥
⎣ ⎦

⎦ ⎣
⎦ ⎣



=
+




(5.7.6) где
1 2
U U ,
Y Y
=
=
Таким образом, искомые соотношения, которые связывают матрицы подсистем и системы, имеют вид
[
]
1 1
2 1
2 2
1 2
1 2
2 1
0
A
B
A
,
B
,
B C
A
B D
C
D C
C
,
D
D D .



=
=






=
=



2. Параллельное соединение двух динамических подсистем (рис.5.7.2).
Для исходных подсистем справедливы те же уравнения (2.6.4).
При параллельном соединении подсистем входное воздействие поступает одновременно на обе подсистемы, т.е.
1 2
U
U
U
=
=
, а выходные сигналы подсистем складываются.
Рис. 5.7.2. Параллельное соединение динамических подсистем
Для системы в целом можно записать уравнения
1 1
1 1 1 2
2 2
2 2
1 1
1 1 2
2 2
2
d
X
A X
BU ,
dt
d
X
A X
B U ,
dt
Y C X
DU
C X
D U
=
+
=
+
=
+
+
+
.
(5.7.8)
154


Вводя в рассмотрение расширенный вектор
[
]
1 2
T
X
X
X
=
состояний, последние уравнения можно записать в следующем виде:
[
]
[
]
1 1
1 1
2 2
2 1
1 2
1 2
2 0
0
X
A
X
B
d
U ,
X
A
X
B
dt
X
Y
C
C
D
D U .
X
2

⎤ ⎡
⎤ ⎡
⎤ ⎡ ⎤
=
+

⎥ ⎢
⎥ ⎢
⎥ ⎢ ⎥

⎦ ⎣
⎦ ⎣
⎦ ⎣ ⎦


=
+
+




(5.7.9)
Из последних уравнений становятся очевидны искомые соотношения между матрицами подсистем и системы
[
]
1 1
2 1
2 1
2 0
0
A
B
2
A
,
B
A
C
C
C ,
D D
D


⎡ ⎤
=
=


⎢ ⎥


⎣ ⎦
=
=
,
B
.
+
(5.7.10)
3. Встречно-параллельное соединение (соединение с обратной связью) двух динамических подсистем.
Рис. 5.7.3. Соединение с обратной связью динамических подсистем
С целью избежания неявных алгебраических уравнений положим, что система не имеет прямой связи, т.е.
1
S
1 0
D
=
Описания исходных подсистем будем полагать прежними, т.е. – уравнения (5.7.4).
При встречно-параллельном соединении (соединение с обратной связью) подсистем выходной сигнал первой подсистемы является входным сигналом для второй подсистемы, выходной сигнал второй подсистемы суммируется или вычитается с входным сигналом системы. Эти условия аналитически можно записать с помощью уравнений
1 2
2 1
U
U Y ,
U
Y
C X
= −
= =
1 1
155

Тогда для всей системы в целом можно записать следующие уравнения:
1 1
1 1
2 2
2 2
2 1 1
1
d
X
A X
B (U Y
dt
d
X
A X
B Y ,
dt
Y C X
=
+

=
+
=
),
или
(
)
1 1
1 1
2 2
2 1 1
2 2
2 2 1 1
1 1
d
X
A X
B U C X
D C X
dt
d
X
A X
B C X ,
dt
Y C X .
=
+


=
+
=
,
Вводя в рассмотрение расширенный вектор состояний
, последние уравнения можно записать в следующем виде:
[
]
1 2
T
X
X
X
=
[
]
[ ]
1 1
1 2 1 1 2 1
1 2
2 1 2
2 1
1 2
0 0
0
X
A
B D C
B C
X
B
d
U ,
X
B C
A
X
dt
X
Y
C
U .
X



⎤ ⎡
⎤ ⎡
⎤ ⎡ ⎤
=
+

⎥ ⎢
⎥ ⎢
⎥ ⎢ ⎥
⎣ ⎦

⎦ ⎣
⎦ ⎣



=
+




(5.7.11)
Из этих уравнений легко найти искомые соотношения между матрицами подсистем и системы. С учетом того, что на входе системы сигнал обратной связи может суммироваться с входным сигналом системы (положительная обратная связь) или вычитаться (отрицательная обратная связь), в обобщенном виде искомые матричные соотношения принимают вид
[
]
[ ]
1 1
2 1 1 2 1
2 1 2
1 0
0 0
A
B D C
B C
B
A
,
B
B C
A
C
C
,
D
.
±
±


,
⎡ ⎤
=
=


⎢ ⎥
⎣ ⎦


=
=
(5.7.12)
Здесь верхний знак соответствует положительной обратной связи, а нижний – отрицательной обратной связи.
1   ...   14   15   16   17   18   19   20   21   22