Файл: Ларин В.Б. Синтез оптимальных линейных систем с обратной связью.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)-> =

 

=

[(Я771я ^ - я л ) Р - ,Г ]- .г - 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 +

и22а.г + ф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)

 

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