
Лекции / Т4_ДГМ
.pdfна него. Поэтому при синтезе наблюдающего устройства потребуется оценивать не только переменные & , но и возмущающее воздействие ur . Нетрудно показать, что попытка избавиться от цветного шума наблюдения за счет дальнейшего повышения порядка объекта не может привести к положительным результатам. Во-первых в этом случае придется оценивать и координаты формирующего фильтра uK , что вряд ли целесообразно. Вовторых, задача становиться некорректной относительно выбора матрицы коэффициентов наблюдающего устройства, поскольку в этом случае шум наблюдений будет рассматриваться как шум, возбуждающий состояние объекта, а измерительный канал оказывается вообще не зашумленным. Следовательно, можно заранее утверждать что значения коэффициентов наблюдателя необходимо выбирать бесконечно большими. Поэтому дальнейшее преобразование объекта с целью приведения его математического описания к ранее рассмотренной задаче необходимо выполнять другим путем. Этого можно добиться, если в качестве измеряемой переменной выбрать фиктивную переменную
{ Y Y |
|
|
|
K |
(4-20) |
можно представить как: { …Z KZ†? |
Z •KqK . |
В результате такого выбора измеряемых переменных измеряемые фиктивные переменные возбуждаются случайными процессами типа «белого шума». Таким образом, выполнение преобразования привели математическое описание объекта к виду:
|
|
|
|
•q |
; |
|
|
{ …Z KZ†? |
Z |
|
r |
•KqK •, |
(4-21) |
||
|
|
|
|
|
|
когда сам объект и его измеряемые переменные возбуждаются случайными процессами типа «белого шума». Соответствующая (4-21) структурная схема показана на рисунке 4.8.
54

|
|
ωη |
|
ωµ |
|
|
|
|
|
|
|
G |
|
Gµ |
u |
B |
1 x |
CA-AµC |
W |
|
p |
|
||
|
|
A |
|
|
|
|
CB |
|
|
Рисунок 4.8
При таком описании объекта оценка расширенного вектора состояния может быть выполнена с помощью наблюдающего устройства, построенного в соответствии с уравнением:
\ |
\ |
|
|
K |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
]y7Y Y < |
|
|
|
|||||
|
…Z Z† |
Z z |
|
|
|
|
|||
|
|
K |
\ |
|
|
, |
|
(4-22) |
|
где матрица коэффициентов ] |
имеет размерность 2 4. |
Y Y |
|||||||
Уравнение (4-22) аналогично уравнению (4-8), так как |
|
K |
|||||||
является |
вектором |
{, |
сформированным через реально измеряемые |
||||||
переменные Y , а выражение |
\ |
|
|
|
|||||
|
|
‡ |
|
K |
|
|
|
||
|
{ …Z Z† Z |
|
|
представляет оценку этого вектора.
Следовательно, для определения матрицы оптимальных коэффициентов наблюдающего устройства ] можно воспользоваться уравнениями, аналогичными (4-13), (4-14), которые для рассматриваемого случая будут иметь вид:
55

] / Z KZ , •Kt •K, + ,
где / - решение матричного уравнения Риккати следующего вида:
0 / / , •t •, /…Z KZ†,…•Kt •K,†+ …Z KZ†/.
Структурная схема расширенного объекта совместно с наблюдающим устройством представлена на рисунке 4.9.
ωµ
|
|
ωη |
Gµ |
|
fµ |
|
|
|
|
1 |
|
||
|
|
|
|
p |
|
|
|
|
G |
|
Aµ |
|
|
u |
|
|
x |
y |
|
|
B |
1 |
C |
|
|||
|
p |
|
|
|
||
|
|
A |
|
|
Aµ |
p |
|
|
|
|
|
||
|
|
L |
|
e |
(-) |
W |
|
|
x |
|
(-) |
|
|
|
B |
1 |
|
W |
||
|
p |
|
CA-AµC |
|
||
|
|
A |
|
|
|
|
|
|
CB |
|
|
|
Рисунок 4.9
Необходимо отметить, что с целью устранения операции дифференцирования измеряемых переменных, структурную схему наблюдающего устройства целесообразно преобразовать к виду рисунок 4.10.
56

y(p)
|
-LA |
L |
|
u(p) |
B-LCB |
1 |
x(p) |
|
p |
|
|
|
|
A-L(CA-A C) |
|
Рисунок 4.10
Важным этапом проектирования оптимальных систем управления со случайными воздействиями является оценка их технических характеристик. За меру качества стабилизации управляемой переменной на заданном уровне при действии на систему случайного процесса q qr qK , обычно принимают значение дисперсии или корень квадратный из неё, то есть так называемый «стандарт». В практических приложениях часто бывает необходимо знать значение дисперсии не только для управляемой, но и для других переменных состояния системы.
С целью расчёта величин дисперсий переменных состояния, при действии на объект случайных воздействий, структурную схему системы целесообразно представить в виде, приведённым на рисунке 4.11.
На основании этой схемы векторно-матричное описание системы, замкнутой через наблюдающее устройство, может быть представлено как:
|
|
|
|
ˆ q |
|
|
|
|
|
, |
(4-23) |
||||
з |
|
з з |
+ |
з у |
з |
где отдельные матрицы определены следующим образом
57

|
|
|
|
|
|
|
|
|
|
ωη |
|
|
|
|
ωµ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
G |
|
|
|
|
Gµ |
|
|
|
|
uy |
|
-1 |
|
u |
B |
|
|
1 |
x |
CA-AµC |
|
W |
|
||
|
|
|
|
Wз (0) |
|
|
|
|
p |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
-1 |
|
|
CB |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
e |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
1 |
CA-AµC |
|
W |
-1 |
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
CB |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
K |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 4.11 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
* |
|
|
|
|
|
|
||
|
|
з ‰] Z KZ |
* ] Z KZ‰ 4 4 ; |
|
|
||||||||||||
|
|
|
|
|
{ |
з+ 0 |
|
|
|
|
• |
0 |
|
( |
4 ; |
|
|
|
|
|
|
з o { з+ 0 4 1 ; ˆз |
o0 |
]•K 4 |
|
||||||||||
з |
|
|
|
|
\ |
, |
|
|
; |
|
|
r |
K |
, |
|
|
. |
|
|
|
|
|
|
4 1 |
|
q |
|
|
|
|
|
4 1 |
|
||
|
|
|
|
Известно, что для систем с математическим описанием (4-23) |
|||||||||||||
справедливо уравнение |
|
|
|
|
|
|
|
|
|
|
|
з/ /,з tз 0,
известное как векторно-матричное уравнение Ляпунова. В этом уравнении матрица / имеет размерность 4n×4n при условии, что все переменные состояния возмущаются. Если представить матрицу / в блочном виде
58

/ |
/ |
|
|
/ o/ |
/ , |
|
|
то диагональные элементы матрицы / , размером 2 2 , |
будут |
||
искомыми значениями дисперсий вектора состояния |
. |
|
|
Матрица tз ˆзtˆз,, размером |
4 4 , |
определяет с |
какой |
интенсивностью и как возмущаются переменные состояния расширенного объекта и зашумляются измеряемые переменные.
4.5 Пример разработки оптимальной системы для двухмассового |
||||||||
упругого электромеханического объекта со случайными воздействиями. |
||||||||
В качестве объекта управления рассматривается двухмассовый |
||||||||
электромеханический объект, математическое описание которого в виде |
||||||||
нормированной структурной схемы представлено на рисунке 4.12. |
|
|||||||
|
|
|
|
|
|
∆MC =x5 |
1 |
wη |
|
|
|
|
|
|
∆MУ =x3 |
Tη p |
|
∆u |
|
∆eТП |
1/ρ |
∆iЯЦ =x1 |
1 |
∆ω1=x2 |
|
∆ω2 =x4 |
|
1 |
1 |
||||||
У |
k |
=1 |
ЯЦ |
|
|
|||
|
ТП |
|
TЯЦp |
|
TM1p |
TC p |
TM2p |
|
|
|
|
|
|
||||
|
|
|
ρЯЦ |
|
∆uТГ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
wµ |
1 |
uµ |
uµ ∆uТГ |
|
|
|
|
|
|
Tµ p |
|
|
|
|
Рисунок 4.12
На выходную переменную объекта случайным образом действует момент статических сопротивления xŠE . Его изменение соответствует
59
экспоненциальному случайному |
процессу, имеющему характеристику |
||||
‹€ m |
|
|
η |
||
спектральной плотности вида: A |
|
|
v |
||
|
|
. |
|||
|
1+ ω2Tη2 |
Такой процесс может быть получен на выходе линейной системы первого порядка, возбуждаемой «белым» шумом qŠr с уровнем интенсивности Œr€, которая описывается как
|
1 |
ŠE |
|
1 |
r |
|
Šс |
|
|
|
. |
||
x |
@r x |
@r qŠ |
Доступной измерению является всего одна переменная объекта управления, которой является скорость вращения исполнительного двигателя. На её измерения аддитивно наложена случайная помеха, сигнал которой может быть сформирован на выходе линейной системы, возбуждаемой «белым» шумом» qŠK с интенсивностью ŒK€
K |
1 K |
1 |
K |
|
@K 9 |
@K qŠ . |
|
9 |
Для заданного объекта требуется разработать систему управления,
оптимальную с точки зрения минимизации критерия
= Vvwxy ?; .
Следовательно, необходимо спроектировать оптимальное наблюдающее устройство и оптимальный регулятор.
|
Пусть параметры, характеризующие объект управления, имеют |
|||||||||||
следующие |
|
значения: |
|
яц 0,1; @яц 0,05 |
с. ; @м 1,5 |
с. ; @м |
||||||
0,07 |
с |
. ; @ 5,7 |
с |
. ; @r 16 |
с |
. ; @K 1 |
с |
. ; Œ€ 4,4 • 10 |
+: |
; Œ€ 0,125 • 10 |
+•. |
|
|
с |
|
|
r |
K |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
Поставленная задача может быть условно разделена на две части: определение структуры оптимальной системы, включающей в себя
наблюдающее устройство и регулятор; расчет значений коэффициентов 60

наблюдающего устройства ] и коэффициентов оптимального регулятора *. Чтобы правильно составить структурную схему наблюдателя, при такой постановке задачи, необходимо пользоваться математическим описанием расширенного объекта, а в качестве доступной измерению использовать фиктивную переменную "{", определяемую как
{„ K„
Вэтом случае уравнение наблюдателя имеет вид:
|
|
|
|
\ |
|
|
|
\ |
|
|
|
|
|
|
|
|
|
|
|
‡ |
|
|
, |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
] { { |
|
|
|||||||||||||||||||
где оценка фиктивной переменной определяется как |
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
‡ |
|
|
|
|
|
|
K |
\ |
|
|
|
|
. |
|
|
|
|
|
|
|||
|
|
|
|
|
{ …Z Z† Z |
|
|
|
|
|
|
||||||||||||||||
Для расширенного объекта матрицы , |
, Z, K с учётом введённых на |
||||||||||||||||||||||||||
рис.4.12 обозначений индексов переменных состояния, получат вид: |
|||||||||||||||||||||||||||
P @1яц |
|
яц1@яц |
|
|
0 0 0S |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||||||||
1 |
0 1 |
|
0 0 |
|
P |
|
|
S |
|
|
|
|
0 |
, |
|
|
|||||||||||
|
|
яц@яц |
|
|
|
|
|
|
|||||||||||||||||||
@м |
|
1 |
|
|
@м |
1 |
|
|
|
|
|
|
|
|
|
P1S |
; K 1 . |
||||||||||
0 |
|
|
0 |
|
0 |
; |
|
|
0 |
|
; Z |
|
|
0 |
|
||||||||||||
@с |
@с |
O |
|
|
|
|
|
||||||||||||||||||||
|
|
1 |
|
|
|
1 |
|
|
|
0 |
|
R |
|
|
|
|
0 |
|
|
@K |
|||||||
0 |
0 |
|
|
|
|
|
0 |
|
|
|
|
N |
|
0 |
|
Q |
|
|
N0Q |
|
|
||||||
|
@ |
|
|
@1 |
|
R |
|
|
|
|
|
|
|||||||||||||||
O |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|||||
N 0 |
0 |
|
|
м |
0 |
|
м |
Q |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
0 |
|
|
@r |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Отсюда следует, что Z 0 и |
|
,” |
?| |
a |
,м# |
?| |
a |
|
|||||||||||||||||||
{ …Z Z† ,м# ?| a |
|
|
|||||||||||||||||||||||||
‡ |
|
|
|
|
K |
\ |
|
|
|
|
|
|
|
|
|
|
|
|
: |
|
. |
||||||
|
|
|
|
|
|
|
|
|
|
|
С учётом полученного результата по формированию фиктивной переменной {‡ структурная схема оптимальной системы, включающей в
61
себя наблюдающее устройство и регулятор, приведена на рисунке.4.13. Как следует из этой схемы, чтобы исключить операцию дифференцирования, при
формировании { , сигнал разности между { и {‡ вводится через соответствующие коэффициенты ] как на вход интегрирующих звеньев, входящих в состав математического описания объекта, так и на их выход.
62

∆MC = x5 1 |
wη |
Объект |
|
|
|
|
|
∆MУ = x3 |
Tη p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∆u |
∆eТП |
1/ρЯЦ ∆iЯЦ = x1 |
|
∆ω1 = x2 |
|
|
∆ω2 = x4 |
|
|
||
|
1 |
1 |
1 |
|
|
||||||
kТП =1 |
TЯЦ p |
|
|
TM1 p |
TC p |
TM2 p |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
ρЯЦ |
|
|
∆uТГ |
|
|
|
|
|
|
k2 + k4 |
∆ uУ |
wµ |
1 |
|
uµ |
1 |
|
|
|
|
|
+1 |
|
Tµ |
p |
|
Tµ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
Tη l51 |
|
|
ρЯЦTЯЦ l11 |
|
l11 |
TM1l21 |
l21 |
TCl31 |
l31 |
TM2l41 |
l51 |
l41 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
Tη p |
kТП =1 |
|
1/ρЯЦ |
|
1 |
|
1 |
|
|
|
1 |
|
|
TЯЦ p |
|
TM1p |
|
TC p |
|
|
|
TM2 p |
||
|
|
|
|
|
|
|
|
||||
|
|
|
ρЯЦ |
|
|
|
|
|
|
|
|
Регулятор |
|
|
|
|
|
|
|
|
|
1 |
|
k1 |
|
|
|
|
|
|
|
|
|
Tµ |
|
|
|
|
|
∆iЯЦ = x6 |
|
|
|
|
1 |
|
|
|
∆ ω1 = x7 |
|
|
|
|
|
|
|
|||
k2 |
|
|
|
|
|
|
|
TM1 |
|
||
|
|
|
|
|
|
|
|
|
|
||
k3 |
∆ MУ = x8 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
TM1 |
|
|
|
∆ ω2 = x9 |
|
|
|
|
|
|
|
|
||
k4 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
k5 |
∆MC = x10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 4.13
63