Мехатроника и динамика мини-роботов
.pdfПрименяя обратное преобразование Лапласа к уравнению (2.6), получаем
t L 1 pI A 1 0 . |
(2.7) |
Обозначая |
|
Ф t L 1 pI A 1 , |
|
уравнение (2.7) запишем в виде |
|
t Ф t 0 . |
(2.8) |
Матрица Ф t называется расширенной матрицей перехода
системы.
В теории обыкновенных дифференциальных уравнений было известно, что решение уравнения (2.5) имеет вид
t eA t 0 . |
(2.9) |
Сравнивая формулы (2.8) и (2.9), находим
Фt eA t ;
Фt L 1 sI A 1 ,
что позволяет вычислить Ф t двумя способами. Обращаясь вновь к описанному выше примеру 1, находим
p |
1 |
|
pI A |
|
. |
b |
p a |
|
50
Тогда
Ф p pI A 1 |
1 |
p a |
1 |
|
|
|
. |
||
|
b |
|||
|
p p a b |
p |
||
Пусть матрица A имеет действительные и различные собственные значения
a |
1 |
a2 4b; |
|
1 |
2 |
2 |
|
|
|
||
2 |
a |
|
1 |
a2 4b, |
|
2 |
|
2 |
|
где a2 4b. Тогда
|
|
|
1 |
1 a e 1t |
2 a e 2t |
|||||
|
|
|
|
|
|
|
|
|
|
|
Ф t |
2 |
|
|
b e |
1t |
e |
2t |
|
||
|
a |
|
4b |
|
|
|
||||
1e 1t 2e 2t .
Если собственные значения матрицы A являются комплексными числами, то матрица перехода Ф t принимает вид
|
|
at |
|
|
|
a |
|
|
1 |
|
|
at |
|
|
|
|
|
|
e |
2 (cos 0t |
sin 0t |
|
e |
2 |
sin 0t |
|
|
||||||||||
|
2 0 |
|
|
0 |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
||
Ф t |
|
|
b |
|
at |
|
|
|
at |
|
|
|
|
a |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
e |
2 |
sin t |
e |
2 (cos t |
|
|
sin t) |
|
|
|||||
|
|
|
|
|
||||||||||||||
|
|
0 |
|
|
|
0 |
|
|
|
0 |
|
|
2 0 |
0 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где
0 b a42 .
Определение Ф t при использовании теоремы разложения Сильвестра рассматривается в курсе математики.
51
Способ определения переменных состояния с использованием схемы системы в переменных состояния: схема составляется из интеграторов, усилителей, суммирующих устройств. Выходы интеграторов могут быть выбраны в качестве координат (переменных) состояния системы. Схема в переменных состояния дает наглядную физическую интерпретацию координат системы и описывает их взаимную связь. Схему в переменных состояния можно составить непосредственно по заданному дифференциальному уравнению или по изображению Лапласа. Применяя преобразование Лапласа к обеим частям дифференциального уравнения примера 2.1, получим
p2 ap b X p M p p a x 0 ax 0 ,
откуда |
|
|
|
|
|
|
|
|
|
X p |
|
|
p 2M p |
|
p 1 1 ap 1 x 0 |
|
p 2 x 0 |
. |
|
1 |
ap 1 bp 2 |
1 ap 1 bp 2 |
1 ap 1 bp 2 |
||||||
|
|
|
|
||||||
Схема системы в переменных состояния, которая непосредственно следует из формулы, изображена на рис. 2.2.
Рис. 2.2. Схема в переменных состояния для системы второго порядка для примера 2.1
Выход x x1 является суммой трех составляющих, возникающих соответственно за счет m t , x1 0 и x2 0 . Из схемы сис-
темы и формулы видно, что составляющая за счет x1 0 равна
52
p 1 1 ap 1 x1 0 1 ap 1 bp 2
и составляющая за счет x2 0 равна
x2 0
1 ap 1 bp 2 .
2.3. Схемы программирования для непрерывных систем
Схемы непрерывных систем в переменных состояния совпадают со схемами этих систем на аналоговых вычислительных машинах. Схема системы в переменных состояния может быть составлена по передаточной функции системы тремя различными способами:
1)прямого программирования;
2)параллельного программирования;
3)последовательного программирования.
Для иллюстрации различных методов составления схемы системы в переменных состояния рассмотрим систему с передаточной
функцией W p вида |
|
|
|
|
|
|
W p |
Y p |
|
p2 3p 2 |
|
. |
(2.10) |
M p |
p p2 7 p 12 |
|
2.3.1. Прямое программирование
Уравнение (2.10) можно записать в виде
W p |
Y p |
|
p 1 3p 2 2 p 3 |
|
M p |
1 7 p 1 12 p 1 |
|
||
|
|
|
||
или |
|
|
|
|
Y p p 1 3p 2 2 p 3 E p , |
(2.11) |
|||
53
где |
M p |
|
|
E p |
, |
|
|
|
|
||
1 7 p 1 12 p 2 |
|
||
откуда |
|
|
|
E p M p 7 p 1E p 12 p 2E p . |
(2.12) |
||
В схеме системы в переменных состояния, показанной на рис. 2.3, используются уравнения (2.11) и (2.12). Переменными состояния
являются x1, x2 , x3.
Рис. 2.3. Схема в переменных состояния для прямого программирования
Дифференциальные уравнения для переменных состояния могут быть найдены из рассмотрения схемы системы. Предполагая, что входное воздействие является ступенчатой функцией, получим
m 0;
x1 x2;
x2 x3;
x3 m 12x2 7x3.
54
Отсюда
m
x1 ;x2x3
0 |
0 |
0 |
0 |
|
|
|
|
|
|
A 0 |
0 |
1 |
0 |
. |
|
0 |
0 |
1 |
|
0 |
|
|||
1 |
0 |
12 |
7 |
|
Выход y t является линейной комбинацией переменных состояния
y t 2x1 t 3xx t x3 t .
2.3.2. Параллельное программирование
Запишем равенство (2.10) в виде суммы дробно-рациональных функций:
Y s |
|
1 |
|
2 |
|
3 |
. |
M s |
6s |
3 s 3 |
2 s 4 |
Схема системы в переменных состояния (рис. 2.4) следует из этого выражения непосредственно.
Система дифференциальных уравнений первого порядкаимеетвид m 0;
x1 m;
x2 m 3x2 ; x3 m 4x3,
где предполагается, что входное воздействие является ступенчатой функцией.
55
Рис. 2.4. Схема в переменных состояния для параллельного программирования
Матрица коэффициентов |
|
|
|
|
0 |
0 |
0 |
0 |
|
|
|
|
|
|
A 0 |
0 |
0 |
0 |
. |
|
0 |
3 |
0 |
|
1 |
|
|||
1 |
0 |
0 |
4 |
|
Так же, как и ранее, выход y t является линейной комбинацией координат:
y t 16 x1 t 23 x2 t 32 x3 t .
2.3.3. Последовательное программирование
Равенство (2.10) записываем в виде произведения дробно-рацио- нальных функций, откуда непосредственно следует сумма системы в переменных состояния, показанная на рис. 2.5.
56
Рис. 2.5. Схема в переменных состояния для последовательного программирования
Система дифференциальных уравнений первого порядка имеетвид
m 0;
x1 4x1 x2 x3;
x2 3x2 x3;
x3 m.
Здесь также предполагается, что входное воздействие является ступенчатой функцией. Матрица коэффициентов
0 |
0 |
0 |
0 |
|
|
4 |
|
|
|
A 0 |
1 |
1 |
. |
|
|
0 |
3 |
1 |
|
0 |
|
|||
1 |
0 |
0 |
0 |
|
Выход y t дается линейной комбинацией координат:
y t 2x1 t 2x2 t x3 t .
Предыдущее рассмотрение иллюстрирует также описание системы различными системами переменных состояния.
57
2.4. Схемы программирования для дискретных систем
Для дискретных систем схемы в переменных состояния имеют почти такой же вид, как и схемы моделирования этих систем. Последние могут быть составлены по передаточной функции дискретной системы тремя различными способами:
1)прямого программирования;
2)параллельного программирования;
3)последовательного программирования.
Схема системы в переменных состояния состоит из прерывателей, элементов фиксации и задержки, суммирующих устройств, усилителей или потенциометров.
Для иллюстрации различных способов составления схемы системы в переменных состояния рассмотрим дискретную систему, характеризуемую передаточной функцией
D p |
1 ap 1 bp 2 |
, |
(2.13) |
|
1 cp 1 dp 2 |
|
|
структурная схема которой показана на рис. 2.6.
1 az 1 bz 2
1 az 1 dz 2
Рис. 2.6. Структурная схема простой дискретной системы
2.4.1. Прямое программирование
Схема системы в переменных состояния, составленная при использовании стандартной процедуры, показана на рис. 2.7. Переменные x1 и x2 характеризуют состояние объекта и представляют собой выхо-
ды соответствующих фиксаторовиэлементовзадержки.
58
Рис. 2.7. Схема в переменных состояния для прямого программирования
Из рассмотрения схемы системы в переменных состояния получаем дифференциальные уравнения состояния
x1 0;
x2 0
и уравнения переходных состояний
x1 nT x1 n 1T x2 nT ;
x2 nT dx1 nT cx2 nT m nT ;
m nT r nT .
Уравнения переходных состояний описывают изменение координат системы в моменты квантования и определяют начальные условия для каждого перехода системы в другое состояние. Выходной
сигнал y t является линейной комбинацией координат x1 и x2 и входного воздействия m :
y t b d x1 t a c x2 t m t .
59
