Лекции / Т4_ДГМ
.pdf• |
ˆ |
• |
|
ˆ |
− LB1)u(t) + (A21 − LA11)Y(t) + LY(t) |
(4-10) |
|
W(t) = (A22 |
− LA12 )W(t) + (B2 |
После не сложных преобразований структурную схему наблюдателя можно представить как
y
|
A - |
11 |
L |
||
|
21 |
LA |
|
||
u |
B2-LB1 |
|
|
1 |
W |
|
|
|
|
||
|
|
|
p |
|
|
|
|
|
|
A22 -LA12 |
|
На основании преобразованной структурной схемы динамика наблюдающего устройства этого типа будет определяться характеристическим уравнением вида
|
det( pI - A22 + LA12 ) = 0 |
|
Для |
наблюдателя полного порядка, как ранее |
было показано, |
характеристическое уравнение имеет вид |
|
|
|
det( pI - A + LC) = 0 |
|
4.2.3 |
Пример проектирования редуцированного |
наблюдающего |
устройства
Требуется спроектировать редуцированное наблюдающее устройство для объекта управления, заданного следующей структурной схемой. Известно, что доступной измерению является скорость исполнительного органа ω2 = x2 . Требуется спроектировать наблюдающее устройство для оценки переменных ω1 = x3 и Mу = x2 , т.е. определить структуру и
коэффициенты наблюдающего устройства. Решая подобную задачу необходимо работать с векторно-матричным описанием объекта, причём 44
обозначения переменных состояния необходимо выполнить соблюдая следующее правило: сначала индексируются переменные доступные измерению, а затем не измеряемые. Обозначим переменные состояния объекта как показано на рисунке.
|
-1 |
|
|
|
|
x2 |
|
x1 |
M(p) |
1 |
x3 |
|
1 |
1 |
|||
|
Tм1p |
|
|
Tcp |
Tм2p |
|
||
|
объект |
|
|
-1 |
|
|
||
|
|
|
|
|
|
|||
|
|
l11 |
|
y1 |
||||
|
|
W1 = x2 |
|
|||||
|
1 |
|
|
1 |
|
|
|
|
|
-T |
|
|
p |
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
- |
111 |
|
|
|||
|
|
Tм2 |
|
|
||||
|
|
|
|
1 |
|
|
|
|
|
|
|
Tc |
|
|
|
|
|
|
|
-( 1 |
|
+ |
|
121 ) |
|
|
|
|
Tм1 |
|
Tм2 |
W2 = x3 |
|||
|
1 |
|
|
1 |
|
|
||
|
Tм1 |
|
|
p |
|
|
|
|
|
|
|
l21 |
|
|
Принимая во внимание уравнение редуцированного наблюдателя
• |
ˆ |
• |
ˆ |
− LB1)u(t) + (A21 − LA11)Y(t) + LY(t), |
|
W(t) = (A22 |
− LA12 )W(t) + (B2 |
определим для нашего объекта все матрицы, входящее в это уравнение. С учётом введённых обозначений переменных состояния векторно–матричное описание рассматриваемого объекта получит вид:
45
|
• |
|
|
|
0 |
|
|
|
|||||
|
x1 |
|
|
|
|
|
• |
|
|
|
|||
|
|
|
|
1 |
||
x2 |
|
= |
− |
|
||
T |
||||||
|
• |
|
|
|
|
c |
x3 |
|
|
|
0 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Tм2 |
|
x1 |
|
0 |
|
|
||||||||
|
1 |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
x |
|
|
+ |
|
|
|
||
|
|
0 |
|
|
|
|
0 |
∆ |
M |
||||||
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
2 |
|
1 |
|
|
|
|
|
|
1 |
|
|
Tc |
|
x |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|||||||
− |
|
|
0 |
|
|
|
3 |
|
T |
|
|
|
|||
|
|
|
|
|
|
|
|
м1 |
|
|
|
||||
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
м1 |
|
|
|
|
|
|
|
|
|
|
|
Из этого описание несложно определить, что
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
1 |
|
|
A |
|
= 0; A |
|
|
= |
|
1 |
0 |
; L = |
l11 |
|
; B = 0 ; B |
|
= |
1 |
|
; A = − |
|
|
; |
||||
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
Tc |
||||||||||||||
|
11 |
|
|
|
12 |
|
|
|
|
|
1 |
2 |
|
|
|
21 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
Tм2 |
|
|
l21 |
|
|
|
|
|
0 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Tм1 |
|
|
|
|
|||
|
|
|
0 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
A |
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
= |
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
22 |
− |
1 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
м1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вычислим коэффициенты наблюдающего устройства, исходя из желания обеспечить его характеристическому уравнению наперёд заданное распределение корней, например стандартное распределение Баттерворта. С этой целью определим характеристическое уравнение наблюдателя и составим систему линейных уравнений, приравнивая его коэффициенты при одинаковых степенях оператора «р» к соответствующим коэффициентам стандартного распределения Баттерворта второго порядка:
|
|
|
|
|
|
) = p2 |
|
|
l |
1 |
|
1 |
|
l |
21 |
|
= p2 |
|
|
p + ω2 |
||||
det( pI - A |
22 |
− LA |
12 |
+ |
|
11 |
p + |
|
|
|
+ |
|
|
+ 1,4ω |
0н |
|||||||||
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
Tм2 |
|
|
|
|
|
|
|
|
|
0н |
|||||||
|
|
|
|
|
|
|
|
|
Tc Tм1 |
|
Tм2 |
|
|
|
|
|||||||||
l11 |
=1,4ω0н |
; |
|
1 |
|
+ |
|
l21 |
|
= ω02н . |
|
|
|
|
|
|
|
|
|
|
||||
|
TcTм1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Tм2 |
|
|
|
TcTм2 |
|
|
|
|
|
|
|
|
|
|
|
Следовательно, коэффициенты наблюдающего устройства необходимо
|
|
|
|
|
|
|
ω2 |
|
1 |
|
|
выбрать равными: l |
=1,4ω |
0н |
T |
; l |
21 |
= |
− |
|
|
T T . |
|
|
|||||||||||
11 |
|
м2 |
|
|
0н |
|
|
|
c м2 |
||
|
|
|
|
|
|
|
|
|
TcTм1 |
|
Чтобы корректно изобразить структуру редуцированного наблюдающего устройства целесообразно определить значение матриц, входящих в его описание, т.е представить описание в виде:
46
|
|
|
|
|
|
|
l |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
11 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
1 |
|
|
|
|
|
||
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
wˆ |
1 |
|
= |
|
|
Tм2 |
|
|
Tc |
wˆ1 |
|
+ |
1 |
|
|
|
+ − |
|
|
y + |
l11 |
|
y& |
|||||
|
|
|
|
∆М |
|
|
||||||||||||||||||||||
& |
|
|
|
|
1 |
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
Tc |
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
21 |
|
|
wˆ2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
wˆ |
2 |
− |
|
|
+ |
|
|
|
0 |
|
|
T |
|
|
|
0 |
|
l21 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
T |
|
|
T |
|
|
|
|
|
м1 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
м1 |
|
|
м2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4.3 Проектирование линейных оптимальных систем для
управления объектами со случайными возмущениями типа «белого
шума».
На реальные объекты управления часто действуют случайные возмущающие воздействия, а измерения физических координат сопровождается шумами. Кроме этого не все переменные состояния объекта доступны измерению. Поэтому одной из распространенных задач, решаемых при разработке системы управления, является задача линейного оптимального управления объектами со случайным возмущением.
Рассмотрим класс объектов со скалярным входом и одной управляемой переменной, математическое описание которых может быть представлено
системой векторно-матричных уравнений: |
|
|
|
||||
|
|
|
|
r |
; |
|
|
Х |
q |
|
|
||||
|
Y ZХ |
qK |
; |
s, |
|
(4-11) |
|
|
|
Х . |
|
|
|
|
|
где - выходная переменная; qr -векторный случайный процесс |
|||||||
типа «белого шума», возмущающий состояние X(t); |
qK -белый шум |
||||||
наблюдений, аддитивно наложенный на измеряемые сигналы. |
|||||||
Совместный |
случайный |
процесс |
qr |
и |
qK являются |
||
некоррелированным и характеризуется блочной матрицей интенсивности |
|||||||
|
|
|
t |
0 |
|
|
|
|
|
t o 0 |
t |
|
|
||
|
|
|
|
|
|
|
47 |
Доступным измерению является «r» переменных состояния из общего числа n. Задача построения линейной оптимальнойсистемы, включающей в себя наблюдатель и обратные связи по восстановленным переменным, состоит в определении такого управления u Y , при котором критерий:
(4-12)
достигает минимума при удержании выходной переменной z в нулевом заданном состоянии zз=0
Положительно определенные весовые коэффициентыты и определяют вклад в критерий качества средних значений квадрата отклонения управляемой переменной от предписанного (нулевого) значения и квадрата управления который при этом затрачивается. Сформулированная задача решается на основе принципа разделения. Применение данного принципа позволяет доказать, что матрица коэффициентов «*» оптимального регулятора, замкнутого по полному вектору состояния остается неизменной и в том случае, когда замыкание системы производится по оценкам вектора состояния \ , восстановленным наблюдающим устройством, параметры которого определены из условия минимума среднего квадрата ошибки восстановления. Таким образом, оптимальное управление, обеспечивающее
минимум критерия (4-12) формируется в виде:
* , где матрица коэффициентов обратных связей «*» как и при решении задачи
детерминированного оптимального линейного уравнения определяется. |
|
выражением: |
* + ,( , |
|
|
в котором значения матрицы «(» являются единственным неотрицательным |
|
симметричным решением алгебраического уравнения Риккати: |
|
|
0 , ( + ,( ,( ( |
48
При необходимости стабилизировать выходную переменную « » на уровнях, отличных от нулевого, требуемое управление может быть определено как:
\ |
|
0 |
|
* { |
, |
||
|
з+ |
з |
где {з – передаточная функция замкнутой системы от входа до управляемой переменной, - задающее воздействие. Оценка полного вектора состояния объектаз \ выполняется установившемся оптимальным
наблюдателем, который реализуется в соответствии с выражением (4-8):
? ?| ] Y Z?|
Структурная схема линейной оптимальной системы, включающей в себя регулятор и наблюдающее устройство, показана на рисунке 4.5.
uу |
ωη |
|
|
ωµ |
|
-1 |
|
|
|||
Wз(0) |
1 |
x |
|
y |
|
B |
C |
||||
p |
|
|
|||
|
A |
|
|
|
|
-1 |
L |
|
|
|
|
|
|
|
|
||
B |
1 |
x |
-C |
y |
|
p |
|
|
|||
|
A |
|
|
|
|
K |
|
|
|
|
|
1. |
|
|
|
|
Рисунок 4.5 Установившийся наблюдатель является оптимальным в том смысле, что
при определенном выборе ( 4 ) матрицы коэффициентов « ] »
минимизируется среднее значение квадрата ошибки восстановления, } ? ?| , то есть:
49
Vvwx }, @} - wv,- ∞
где @ – наперед заданная положительная весовая матрица. Дифференциальное уравнение ошибки восстановления может быть
получено путем вычитания из первого уравнения (4-11) уравнения (4-8), |
|||||
принимая при этом во внимание, что |
Y Z |
qK т.е |
|||
|
|
|
qr |
|
|
|
\ |
||||
\ |
\ |
|
|
K |
|
|
]Z ]q |
]Z |
|||
В результате, получим |
|
|
qr |
|
|
|
} ]Z } |
b 1 ] |
|
||
|
~qK • |
|
Решение задачи построения установившегося оптимального наблюдателя достигается выбором матрицы коэффициентов наблюдателя как
] /Z,t+ (4-13)
Если обозначить среднее значение вектора ошибок восстановления}€ Мy} z, то в выражении (4-13 ) –матрица
/ xy} }€ } }€ ,z
является матрицей дисперсии ошибок восстановления, которая определяется решением алгебраического уравнения Рикатти:
0 t /Z,t+ Z/ / /, |
(4-14) |
К рассмотренной задаче построения оптимальной системы управления, когда предполагается что шум, возмущающий состояние объекта и шум датчиков измерения не коррелированы между собой и являются случайными процессами типа белого шума, могут быть приведены многие инженерные задачи. Однако в практических приложениях лучшим приближением
50
случайных возмущений действующих на реальный объект является цветной шум, т.е случайный процесс с определенной характеристикой спектральной плотности обычно формируемой путем фильтрации белого шума.
4.4 Проектирование линейных оптимальных систем для управления объектами со случайными возмущениями типа «цветного
шума».
Рассмотрим возможность построения линейной оптимальной системы для объектов управления с цветным шумом внешних воздействий и критерием качества (4-12). Пусть стационарный объект управления порядка « » задается уравнением:
|
|
& |
|
u |
|
|
||
|
|
|
r r |
• |
|
|||
|
& |
|
Y Z& & uK |
|
(4-15) |
|||
|
|
& & |
|
|
|
|
где ur и uK - случайные возмущающие воздействия с определенными статистическими характеристиками, полученными при фильтрации некорелированных белых шумов mr и mK через соответствующие фильтры первого порядка:
ur rur |
•rqr |
|
|
•KqK |
(4-16) |
uK KuK |
|
|
|
|
(4-17) |
Считается, что доступными измерению являются « 4 » переменных состояния объекта управления из общего числа « ».
Векторно-матричная структурная схема, соответствующая такому математическому описанию объекта может быть. представлена как показано на рисунке 4.6.
51
ωη |
|
|
ω |
|
|
|
|
|
µ |
|
|
Gη |
|
G |
|
|
|
|
fη |
|
µ |
|
fµ |
1 |
|
|
1 |
||
p |
|
|
|
p |
|
Aη |
Bη |
|
Aµ |
y |
|
u B0 |
|
1 |
x0 |
C |
|
|
p |
|
|
||
|
|
A0 |
|
|
|
Рисунок 4.6
При таком описании объекта задачу построения оптимальной системы управления стремятся привести к ранее рассмотренной, когда шум возбуждающий состояние объекта и шум измерений являются случайными процессами типа белого шума. С этой целью часть описания, относящуюся к моделированию цветного шума ur a , относят к самому объекту и в дальнейшем рассматривают математическое описание расширенного объекта, которое на основании уравнений (4-15), (4-16) может быть представлено как
|
|
|
|
|
|
|
& |
|
& |
0 |
|
|
|
||
|
|
|
|
|
|
|
r |
|
S |
||||||
~ |
• o |
|
& |
|
r o |
r |
|
|
‚ |
ƒ |
|
|
|||
|
|
|
|
o•r q |
(4-18) |
||||||||||
|
& |
|
|
r |
u |
|
|
|
|
|
|
|
|||
|
r |
|
0 |
|
|
0 |
|
|
|
|
|
||||
u |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Y Z& |
|
0 |
& |
uK |
|
|
|
Q |
|||||
|
|
|
our |
|
|
|
Введём соответствующие обозначения блочных матриц в системе выражений (4-18), которые в общем случае, когда все переменные вектора состояния & возмущаются, будут иметь следующие размерности:
|
o 0 r 2 2 ; |
‚ |
& |
|
|
|
A= |
0 |
ƒ 2 1 ; our |
||||
& r |
|
|
& |
|
|
2 1 ;
52
Z Z& 0 4 2 ; |
|
|
0 |
|
|
|
|
|
|
• o•r 2 . |
|
||||||
С учётом принятых обозначений, уравнения расширенного |
||||||||
объекта можно переписать в виде: |
|
•q |
; |
|
|
|||
|
|
• |
|
|||||
|
„ Z? |
uK . |
r |
|
(4-19) |
|||
|
|
|
|
|
|
, |
где uK - по-прежнему цветной шум.
Соответствующая (4-17) структурная схема показана на рисунке 4.7
|
|
ωµ |
|
|
ωη |
Gµ |
1 fµ |
|
|
||
|
|
|
p |
|
G |
|
Aµ |
u B |
1 |
x |
y |
p |
C |
A
Рисунок 4.7
Таким образом, включение описания формирующего фильтра в состав описания объекта управления изменяет первоначальную формулировку задачи. Теперь расширенный объект возбуждается воздействием типа «белого шума», поэтому коэффициенты жестких обратных связей при формировании оптимального закона управления могут быть как и прежде рассчитаны путем решения соответствующего алгебраического уравнения Рикатти (4-5).
Следует отметить, что расширенный вектор включают в себя не только переменные состояния объекта & , но и возмущения, действующие
53