Файл: Ларин В.Б. Синтез оптимальных линейных систем с обратной связью.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 13.07.2024
Просмотров: 135
Скачиваний: 0
сохраняет свойство оптимальности. Как видно из этой форму лы, влияние характеристик внешних возмущений на структу ру регулятора проявляется наличием в формуле для W ма трицы К -Г -1 , которая, согласно (1.96), определяется так:
К -Т ~ 1= [(H~lQ-'N:,R — НА) Р -*Г ]_ Г _І. (1.106)
Если элементы матрицы спектральных плотностей 5ф внешних возмущений константы, то элементы матрицы Г (см. (1.77)) также константы и соотношение (1.106) запи шется в виде
К -Т -' = [(Я Г’О Г'лу? - НА) р -']_ = К_,
где элементы матрицы V_ не зависят от статистических ха рактеристик внешних возмущений. Таким образом, ре гулятор, определяемый уравнением
[Д*р (s) H-'Q-'C + F_M] и = [К_Р - H -'Q -% R ] х, (1.107)
будет оптимальным при любых внешних возмущениях, эле менты матрицы спектральных плотностей которых кон станты.
Более того, если регулятор, определяемый формулой (1.85) или уравнением (1.97), является оптимальным при возмущениях с матрицей спектральных плотностей (со) = = Г (/со) Г' (— /со), то он будет оптимальным и при возму
щениях |
с |
матрицей спектральных плотностей |
(со) = |
||
= Г (/со) С0 СоГ' (—/со), где |
С„ — невырожденная |
матри |
|||
ца п X п, элементы которой константы. |
|
||||
Действительно, при возмущениях с матрицей спектраль |
|||||
ных плотностей |
(со) |
|
|
||
Ö |
- |
1 = [(H :1Q -% R - |
НА) Р - ’ГС0]_ (ГС0)-> = |
||
|
= |
[(Я71С71я ^ - я л ) Р - ,Г ]- .г - 1= /С _г-1, |
|
т. е. уравнения оптимальных регуляторов при возмущениях
с матрицами спектральных плотностей 5ф (со) и |
(со) со |
впадают. |
|
40
§ 4. |
Н Е К О ТО Р Ы Е Ч А С ТН Ы Е С Л У ЧА И Р Е Ш Е Н И Я З А |
|
Д А Ч И С И Н ТЕ З А С И С ТЕМ С ТА Б И Л И З А Ц И И П Р И |
|
П В Н Е Ш Н И Х В О З М У Щ Е Н И Я Х И т У П Р А В Л Я Ю |
|
Щ И Х В О З Д Е Й С ТВ И Я Х . П Р И М ЕР Ы |
При рассмотрении конкретных задач использование реше ния в виде (1.85) требует выполнения довольно трудоемких выкладок. Поэтому представляется целесообразным иссле довать некоторые частные случаи, когда удается получить более простые формулы для матрицы передаточных функций оптимального регулятора, используя специфику динамиче ских свойств объекта.
Как было показано в § 3, матрица W передаточных функ ций оптимального регулятора определяется соотношением
|
|
|
W = [Д; (s) H -'Q Z 'C + |
К - Т - 'М ] - 1 X |
|
||||
|
|
|
|
X [К -Г -'Р - Я Г 'З Г ’В Д , |
(1.108) |
||||
или |
|
|
|
|
[(/С0+ К+) г - ' м |
|
|
||
|
|
|
W = |
- НВ}-1X |
|
||||
|
|
|
|
X [(/С0 + |
К+) Г- |
НА], |
(1.109) |
||
где |
|
|
|
|
|
|
|
|
|
|
(s) |
= |
fiel Р, |
AN, |
|
|
|
||
Q = |
Др |
(s) В |
|
|
|
|
|||
N = |
|
|
(s) Р-'М, |
|
|
|
|||
|
|
Ар |
|
+ |
|
|
|
|
|
Я ,Я |
= |
Qr'GjSQ-', |
|
|
(1.110) |
||||
G*G = |
N,RN + |
Д; (s) С&р (s), |
|
|
|||||
Г Г , = |
5ф, |
|
|
|
|
|
|||
|
+ |
|
К+ + |
К - |
= (H - ' Q- ' N ^R НА) р - 'г, |
|
|||
а элементы матриц А и В должны удовлетворять требованию |
|||||||||
аналитичности в правой полуплоскости матрицы |
|
||||||||
|
|
|
|
|
|
Z = |
'Р |
|
|
|
|
|
|
|
|
А |
|
|
|
вместе |
с обратной. |
|
|
т. е. det Р = |
|
||||
I. |
|
Пусть объект устойчив, |
Др (s) — гур- |
||||||
вицев |
полином. |
|
|
|
|
|
|||
Тогда требование аналитичности в правой полуплоскос |
|||||||||
ти матрицы |
Z вместе с |
обратной |
можно удовлетворить, |
41
положив А = О, |
В — Ет (Ет — единичная |
матрица |
т х |
|
X т), так как при этом |
|
|
|
|
|
ГР |
— Ml |
|
|
detZ = det |
= detP. |
|
|
|
|
О |
|
|
|
Соотношения (1.110) тогда примут вид |
|
|
||
« = |
= |
н |
= ^ а |
, |
к , + |
Л'. + К-. =. (г'ы ,т > -'і\ |
( м и ) |
а матрица передаточных функций оптимального регулятора определится формулой
W = [A; (s) G7 |
‘C + РС-Г-’Л^Г) - |
1 [ / С Г - ’Р - |
G-'N^R], |
(1 . 1 1 2 ) |
|||||
или |
|
|
|
|
|
|
|
|
|
w = [(Ко + |
К+) г-'м - |
НГ' (Ко + |
К f ) Г - ’Р. |
(1.113) |
|||||
Рассмотрим пример. Пусть число внешних возмущений |
|||||||||
равно |
числу |
управляющих |
воздействий |
(М — матрица |
|||||
/і X п), |
причем |
М ~1 не аналитическая в |
правой полупло |
||||||
скости. |
|
|
|
|
|
|
|
|
|
Покажем, что в этом случае даже отсутствие ограничений |
|||||||||
на управляющие воздействия (матрица С = |
0 ) не позволяет |
||||||||
сделать |
систему |
инвариантной относительно внешних воз- |
|||||||
|
|
|
|
|
/оо |
|
|
|
|
действий, т. |
е. |
етіп = 4 - |
[ |
Sp (К - Д _ ) ds Ф 0 (rj(s) = 0, |
|||||
см. (1 .8 |
6 )). |
|
|
|
|
|
|
|
|
Согласно (1.110), |
|
|
|
|
|
|
|||
|
|
|
GJÜ=,N .RN, |
|
|
(Ы 14) |
|||
т. е. |
|
|
|
|
|
|
|
|
|
|
|
|
G f1 |
= |
GN~'R~]N~'. |
|
|
(1.115) |
|
Подставив |
значение |
СК1 из (1.115) |
в |
(1.112), |
получим |
матрицу передаточных функций оптимального регулятора
W = (К -Г -'М )-1(/С Г _ІР - GN~l), |
(1.1165 |
а формулу (1 . 1 1 1 ) перепишем в виде |
|
tf0+ /(+ + tf_ = ^ G / M - 'r . |
■(1.117) |
42
Отсюда видно, что из-за наличия полюсов в правой полуплоскости у матрицы М ~1 (матрицы G и Г аналитиче
ские в правой полуплоскости, а Ар (s) — гурвицев |
полином) |
||
К— |
т, е. вгпіп |
0 . |
|
Проиллюстрируем |
изложенное числовым |
примером. |
Пусть движение объекта описывается системой дифферен
циальных |
уравнений |
|
|
X+ |
8 л: + |
22х — у — 1 0 у = |
— 2 их — и.2+ |
х + 1 2 х |
+ у + 6у = ~ и 1 + |
и2— 2а.г + ф2, |
где фх и ф2 — некоррелированные стационарные случай ные процессы с единичной спектральной плотностью.
Необходимо найти закон управления, обеспечивающий устойчивость замкнутой системы объект + регулятор и ми нимум функционала е = (х2) -)- (У2).
,В обозначениях, принятых в § 3, имеем
Р = |
s2 + 8 s + 2 2 |
— s — 1 0 |
|
|
|
|
|
|
||||
s + |
1 2 |
s2 |
+ |
6 s |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||||
М = |
's — 2 |
— |
1 |
|
Ap (s) = |
s4 + |
14s 3 -f- |
|
|
|||
— 1 |
s — 2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||
+ |
|
71s2+ |
154s + |
120 = |
(s + 2)(s + 3)(s + |
4) |
X |
(1.118) |
||||
X |
(s + 5), |
S x|) = |
Eo, |
Г |
= E2, |
R = |
E2, |
C = |
0, |
|
||
|
|
|
|
~s3 -|-4s2 —13s—10 |
|
2s — 20 |
|
|
||||
N ^ A p (s) P-'M = |
— 2 s2 — 18s + |
2 |
s3 |
+ |
6 s2 + |
|
||||||
|
|
|
|
|
|
|
|
+ |
7s — 32 |
|
Выполним факторизацию полиномиальной матрицы (?*G = = N^RN.
Воспользуемся алгоритмом Дэвиса [321, который осно
ван на представлении матрицы G~l в виде
к
С“ 1 |
= |
П 7\7У‘, |
(1.119) |
|
|
/=і |
|
где определитель матрицы |
Т0 равен |
константе, а опреде |
|
лители матриц 7 Т 1 (t = |
l, |
2........ /г) |
равны элементарным |
делителям det G. |
|
|
|
Исходная полиномиальная матрица домножается справа на Т(, а слева на 7Д, причем матрица Тс строится так, что бы в результате выполнения этой операции снова получить
43
полиномиальную матрицу. Отметим, что домножение исход ной матрицы справа и слева на Tt и 7 + понижает степень ее определителя, так как det Т ~1 равен элементарному де
ятелю det G, а det ТГ? — элементарному делителю det G*, При факторизации двухмерной матрицы элементарному делителю det G первого порядка s + gi ставится в соответ
ствие одна из матриц
|
1 |
|
0 |
|
|
5 + |
ёі_ |
|
|
Ті |
|
( 1. 120) |
||
ki |
|
|
||
|
|
1 |
|
|
|
s + |
gi |
|
|
или |
|
|
||
|
kt |
|
|
|
|
|
|
|
|
Ti |
|
s + |
gi |
( 1. 121) |
L° |
|
gi J |
||
|
s + |
|
где k[ выбирается так, чтобы в результате домножения ис ходной матрицы на Т{ получить полиномиальную матрицу, что соответствует выбору /г,- из условия делимости элементов первого столбца матрицы, получившейся после умножения на Tt вида (1.120), на s + g t, или из условия делимости эле ментов второго столбца на s + gi, если Т ( имеет вид (1.121) х.
Аналогично, элементарным делителям det G второго по рядка s2 + a(s + bi (случай пары комплексно сопряжен ных нулей det G) ставится в соответствие матрица
S- + Ci;S -j- bi |
|
||
Ті = i |
CjS + |
âj |
|
p s'- -г «iS -f- bi |
|
||
или |
|
|
|
. |
CjS+ dj |
|
|
Tt = |
s2 + |
atS+ |
bi |
|
1 |
|
|
|
|
bt |
|
0 |
s2 + |
o,s + |
|
|
1 Так как gi является нулем определителя матрицы G*G, то два уравнения относительно ki, получаемые из этих условий, оказываются линейно зависимыми.
44