Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Ларин В.Б. Синтез оптимальных линейных систем с обратной связью

.pdf
Скачиваний:
4
Добавлен:
23.10.2023
Размер:
6.26 Mб
Скачать

Подставив с“ (gt) в (2.20), получим

. .

g(gt) (— l)k+l+lmjk (gt)

C l(gl)~

 

(2.23)

Таким образом, элементы вектора ѵ - определяются фор­

мулой (2.17), где постоянные

коэффициенты bt (gt), Cj (gt)

и dj (gt) имеют вид (2.18), (2.23) и (2.19), явно не зависящий от полиномов a t (s) и ß (s). Следовательно, согласно (2.17), (2.16) и (2.13), и векторы и_, /е_ и w не зависят от полиномов a i (s) и ß (s).

Пример. Пусть движение объекта описывается следую­ щей системой линейных дифференциальных уравнений:

х= х -\-Зу и — фр

у= 2у — и — ф2,

где фі и ф2 — стационарные случайные процессы с нуле­ вым математическим ожиданием и матрицей спектральных плотностей

 

 

=

1 T’s3

 

 

 

т*

 

 

 

 

 

 

 

 

0

 

 

 

 

1 — TJS2 J

 

 

Требуется

определить закон управления u =

w1 (s)x-j-

-f-

(s) у, т.

е. передаточные функции w1 (s) и

ш2 (s), та­

ким образом, чтобы при устойчивой замкнутой системе

объект

регулятор минимизировалась величина е =

= с<и 2) +

(ІІ2).

Согласно (2.1), (2.6), (2.7) и (2.10),

1— s

3 •

,

т

Т

= 0

2 — s

 

= .1.

(s) = (1 _ s) (2 — S),

1

— s'

g* (s) g (s) =

(с — s2) (1 — s2) (4 — s2).

п =

s

1

 

 

 

 

Таким образом,

 

 

 

 

 

g(s) =

(V c +

s)(l

+ s) ( 2

+ s),

 

g*(s) =

(V c — s) ( 1

— s) ( 2

— s).

Полином g* (s) имеет два

нуля, общих с Ар (s) (gy = 1,

g2 = 2 ), и третий нуль g3 =

V с.

 

 

60

Подставив в (2.18), (2.23) и (2.19) соответствующие вели­ чины, получим

/ Ы = ° (/, t = \ , 2 ),

 

 

 

(gi)ci (<?і)О=

 

— 3

(Ус

+

1),

 

 

(g3) = О,

 

 

 

 

с 2

 

 

 

 

 

 

 

 

 

 

= 9

 

Гс +

1),

 

 

(g*) =

 

-

71

(Ѵ~с + 2),

 

 

 

 

 

 

<*/(&) =

О

 

 

( / = 1 ,

2

;

 

=2 3 ) ,

 

 

 

или, согласно (2.17),,

Ѵг- W

 

S _

J

 

 

 

12 ( / 7 + 2 )

 

 

 

 

3 ( / 7 + 1)

 

 

 

 

 

9 ( / 7 + 1)

=

ü , _ ( s ; -------s _ j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s —

2

Факторизовав матрицу 5ф, т. е. представив в виде

 

= ГГ*,

где

 

Г =

 

1

 

+

TjS

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

 

 

1 + ТоT 2s

 

 

 

 

и подставив Г в (2.16), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

k i- (s) =

 

 

 

3 (/7 +

1 )тх

 

 

 

 

 

 

 

 

 

 

( /

+

l) ( s —

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее,

 

fa - (s) =

 

 

9 ( / 7 + 1 )

To.

 

 

1 2 ( / 7 + 2 ) T 2

 

 

 

( T 2 + l ) ( s - l )

 

( 2 T o

+

1) (s — 2)'

 

 

согласно

 

(2.13),

 

wx

(s)

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 ( / c + l )

C / s + 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г і +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

12(1

с + 2)

 

 

, +

^

 

+ 6

- ? ^

І

 

+

1)

 

T x

 

9 ( Г с +

 

1)

 

 

 

 

 

 

 

 

 

Гх + 1

 

 

Д 1

1

 

Т2 + 1

 

 

 

2 Г г + 1

 

 

 

 

 

 

 

 

- 9 ( / 7 + 1 )

 

Т' - Т'

 

 

 

 

 

 

(s) =

 

 

 

 

(/7+1)

 

( / +

1 ) ( Г о + 1 )

 

 

 

 

W o

 

 

9

 

1 2 ( / 7 + 2 )

+ 1)

 

 

 

 

 

 

 

 

То +

1

 

 

 

2 Т 3 +

1

( T 2s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S +

/ 7

+

6

— -3

(^

 

j T1!)

Гх +

 

 

 

 

 

 

 

 

9 (/7 +

1 )

 

 

1 2

( / 7 + 2 ) т

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

2 Т 2 +

 

1

1 2

 

 

 

В

частности,

при с =

 

25 и 7\ =

 

Т

2

=

0 значения

 

(s)

 

 

2+ 1

 

 

 

 

и(s) совпадают со значениями передаточных функций

61

оптимального регулятора, полученных при решении приме­ ра, приведенного в работе [13].

Рассмотренный пример хорошо иллюстрирует инвари­ антность решения относительно постоянной матрицы — мно­ жителя справа к матрице Г, которая была доказана в § 3 гл. 1 .

Действительно, матрицу Г можно представить в виде

Ті

о

1

Г = 1 + 7Ѵ

1

+ T lS

О

Т2

О

1 + r 2 s_

1 +

0

1

T 2s

1 -----

о 1—1

0 т2

причем решение задачи не должно зависеть от х1 и т2 (эле­ ментов постоянной матрицы-множителя), что и видно из окончательных формул для wx (s) и w2 (s).

§ 2. С В Я З Ь З А Д А Ч И С И Н ТЕ З А О П ТИ М А Л Ь Н О Й

С И С ТЕМ Ы С ТА Б И Л И З А Ц И И С З А Д А Ч Е Й А Н А ­ Л И ТИ Ч Е С К О Г О К О Н С ТР У И Р О В А Н И Я Р Е Г У Л Я ­ ТО Р О В

В § 3 гл. 1 показано, что передаточные функции оптималь­ ной системы стабилизации не зависят от интенсивности воз­ мущений, действующих по различным координатам, если эти возмущения являются многомерными стационарными случайными процессами с постоянной спектральной плот­ ностью.

Этот результат можно было предугадать, руководству­ ясь следующими соображениями. В [18] показано, что вид оптимального регулятора инвариантен относительно на­ чальных условий. В [14] отмечается, что задача аналити­ ческого конструирования регуляторов может быть решена методами теории оптимальной фильтрации Винера — Кол­ могорова, если начальные условия заменить возмущениями типа 6 -функций. Таким образом, если объект, для которого синтезируется оптимальный регулятор, находится под воз­ действием «-мерного 6 -коррелированного процесса, то мож­ но было ожидать, что вид оптимального регулятора не будет зависеть от интенсивности каждой из компонент внешнего воздействия, что и было показано в § 3 гл. 1 .

Ниже будет показана более сильная связь этих задач, а именно, что решение задачи аналитического конструирова­ ния регуляторов совпадает с решением, приведенным в

62

§1 гл. 2 , когда объект находится под воздействием «-мер­

ного б-коррелированного стационарного случайнного про­ цесса.

Сформулируем в принятых выше обозначениях задачу об аналитическом конструировании регуляторов [18].

Пусть возмущенное движение объекта описывается систе­ мой дифференциальных уравнений

где

 

 

 

 

 

Рх = ти,

 

 

 

 

(2.24)

 

 

 

 

 

 

. . .

&1 ,п

 

 

 

 

 

 

 

Г р + ьа

^ 1 2

 

 

 

 

 

 

Р = -

 

^ 2 1

р + ь22 •••

Ь2,п

1

 

(2.25)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ьп,1

bn,2

••Р + ьпп.

 

 

х

(О*

*/». . .(t)V. -—вектор,

составляющими

кото-

рого являются

координаты объекта,

т =

[тх,

 

 

т1 (і =

1

,

2 , ...,

п) — константы,

и — координата

регу­

лятора,

 

р

— оператор

дифференцирования

^

 

btj

(г, /' =

1

, ..., п) — заданные постоянные числа.

 

 

Задача состоит в том, чтобызаписатьв аналитической фор­

ме F (и,

 

и,

хг,

...,

хп) =

0 закон

регулирования,

который

в совокупности

с

исходными уравнениями

(2.24)

образует

устойчивую систему и гарантирует существование миниму­ ма интеграла [18]

е = J ^ 2 ТіА- + с “ 2 + dt>

где все rk n с — положительные весовые константы.

В [18] показано, что уравнение регулятора в этом слу­

чае имеет вид

П

 

ü + h u =

2 р'кхк,

 

где рк и h — постоянные,

 

*= 1

 

которые, как известно

[21, 27],

полностью

определяются

характеристическим уравнением

замкнутой

системы объект +

регулятор. Таким

образом,

для того чтобы показать тождественность решения задачи аналитического конструирования регуляторов и задачи, сформулированной в § 1 гл. 2 при S,/ (ш) = const, достаточно показать, что уравнения регуляторов имеют одинаковую структуру, а характеристические определители замкнутых систем совпадают.

63

Вектор передаточных функций регулятора для объекта, возмущенное движение которого описывается системой дифреренциальных уравнений (2.1) (матрица Р имеет вид (2.26),

i f/ (со) =

const), определяется

согласно (2.13). Покажем,

что выражение (2.13) в этом случае имеет вид

 

 

 

 

 

 

“' =

7Т 7ГР о.

 

 

 

(2.26)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р о — [Роі>

р 0 2 > • • • > Ропіі

 

 

 

 

 

h> Poi>

Р о 2 > ■ • • > Pon =

const.

 

 

Если элементы матрицы

(со) константы,

то, согласно

(1.102), уравнение регулятора примет вид

 

 

(С _ s3) Д

(5)

 

и =

г)_Р —

 

 

X, (2.27)

 

g* (s)

+

г>_т

g*(s) ■ n*R

причем, как было показано в

§ 3 гл.

1, скаляр

w0 (s)

(с — s3)A

(s)

 

и

элементы

вектора

-

 

-----------

^ + v - m

w —

g* (s)

 

 

 

 

 

 

 

 

 

 

------ j—-

iuR — полиномы от s.

 

 

 

 

 

g* (s)

*

 

 

 

 

 

 

 

 

 

 

Общим знаменателем всех элементов вектора w явля­

ется полином

 

 

 

 

 

 

 

 

 

 

 

 

w0 (s) =

( С - ^ )

А *

(5)

+

Ѵ-т,

 

 

 

 

 

 

 

g*(ß)

 

 

 

 

 

равный

целой

части

от деления

полинома (с — s2) Др (s)

на полином g* (s), так как дробные части слагаемых в этой формуле равны по величине и противоположны по знаку

(см. (1.89)), а скаляр и _ т при т( = const

целой части

не

содержит.

 

 

п,

Поскольку порядок полинома Др (s) =

det

равен

а порядок полинома g* (s), определяемого из разложения

(2.10),

п + 1, то

 

 

 

 

Wo (s) =

s +

h,

где h =

const.

 

 

 

Таким образом,

 

 

 

 

W =

1

 

 

--- r~Tw,

 

 

S +

/l

64

где, согласно (2.27),

вектор

 

 

W = Ѵ—Р -

1

nJR,

g*(s)

 

 

элементы которого

woi = ро/ = const.

Таким образом, для доказательства тождественности оп­ тимальных регуляторов осталось показать равенство ха­ рактеристических определителей замкнутых систем.

Квадрат модуля характеристического определителя замкнутой системы объект + регулятор, возмущенное дви­ жение которой описывается системой уравнений (2.24) и (2.27), согласно (1.105) и (2.10), можно записать так:

A* (s) A (s) =

g* (s) g (s) = n:,Rn -f- (c — s2) A* (s) Ap (s), (2.28)

так как у =

const.

В задаче аналитического конструирования регуляторов характеристический определитель A (s) замкнутой системы

объект +

регулятор

(характеристический

определитель

уравнений

вариационной задачи) определяется из уравне­

ния 118]

 

' Р

0

т

1 1

 

 

 

A* (s) A (s)

2R

— Я*

о

>

 

 

о'

tri

2 (с — sa)

 

где 0 — нулевая матрица размера п X п, о — нулевой вектор-столбец, R = diag [гг, г2, ..., гп}. Воспользовав­ шись формулами вычисления определителя блочной матри­ цы (1.75), получим

A* (s) A (s) = Ар (s) det

' - Р .

о

tri

2 (с — s2)

 

= Ар (s) det

- 2 R P - ]m 2 (с — s2)

или

A* (s) A (s) = (— l)n Ap (s) Ap (s) [/п 'Я ^ Я Я -1/« + 2 (c — s2)].

Согласно (2.10), это произведение можно записать в виде

Л* (s) A (s) = 2 (— 1)п [n^Rti + (с —s2) A* (s) Ар (s)],

1 Определитель этой матрицы с точностью до постоянного множи­ теля совпадает с определителем вспомогательной блочной матрицы, ис­ пользованной в § 3 гл. 1 (см., например, (1.102)).

5 3 - 5 8 2

6 5

что совпадает с (2.28) с точностью до отличного от нуля постоянного множителя. Таким образом, доказана тождест­ венность оптимальных регуляторов двух рассмотренных задач при S {/ (со) = const.

Тесная связь между задачей аналитического конструи­ рования регуляторов и задачами стабилизации объектов при случайных возмущениях существует и в том случае,

когда

элементы S q

(со) матрицы спектральных

плотностей

не

являются константами.

 

 

Действительно,

в § 3 гл. 1

было показано,

что переда­

точные функции оптимальной системы стабилизации не изменятся, если формирующую возмущения матрицу Г умно­ жить справа на постоянную невырожденную матрицу. Повидимому, этот факт связан с тем, что задачу синтеза опти­ мальной системы стабилизации при Sq (со) Ф const можно свести к задаче синтеза большей размерности, если возму­ щающие воздействия рассматривать как решения дополни­ тельной группы дифференциальных уравнений, в правых частях которых стоят возмущения типа «белого шума»,

Докажем это, т. е. покажем, что решение задачи § 1 гл. 2 (2.13) совпадает с решением следующей 2/і-мерной задачи. Пусть возмущенное движение объекта описывается

системой

уравнений

 

 

 

 

 

 

 

 

 

 

 

Р0х0 =

т0и (0 + É,

 

 

(2.29)

где

х0

=

[*! (0,

.... хп (t),

ух (І),

Уп (01' — 2/г-мерный

вектор,

и (t) — управляющее

воздействие,

§ (t) — 2/г-мер­

ный 6-коррелированный случайный процесс

(Si (со) = E on ),

 

 

 

Р о =

y (s)E n

(Г0 =

У (s) Г -1),

 

 

 

 

А

Г0 .

 

 

 

 

 

 

 

 

 

 

 

 

т0 -

[тѵ т2,

. . . .

тп, 0,

. . . , 0]'.

 

 

Здесь Р, тс, хс (t), у (s) и Г совпадают с соответствующи

ми величинами в задаче § 1 гл. 2, а координаты

у2 (t), ...

• • • 1

Уп. (0

с точностью до множителя у

(s) соответствуют воз­

мущающим воздействиям фі (і), ..., ф„ (t).

 

 

 

Необходимо найти передаточную функцию до0Прегулято-

ра и =

w0x,

минимизирующую функционал е = 2

rt (xf) -f

 

c (w2) +

(u2).

 

 

 

 

 

t=l

 

+

 

 

 

 

 

 

 

66

Применив для

решения

сформулированной

 

2/г-мерной

задачи методику §

1 гл. 2, получим, согласно

(2.13), 2/г-мер-

ный вектор искомых передаточных функций

 

 

 

 

wn=

(с — s2) Аре (5)

 

■ -1

 

 

1

 

 

Щ-Щ

 

 

 

 

 

 

 

 

+

Ѵ0- Р 0

 

e0(s) no*R

 

 

 

 

 

где АРо(s) =

Av (s) Ар (s),

g 0 (s) =

Дѵ (s) g

(S),

 

n0 =

[Av (s)/z;

Л

 

 

 

 

 

 

 

 

 

 

 

 

 

0,.... 0]' — 2/г-мерный вектор, v0-

= [Ü_;

k -]

(согласно

(2.15)

и (2.16)), R0 =

diag {rv ...,

ra,

 

 

rp

0

-i

0, ...,

0}

=

 

 

" .

Здесь Ay (s)

=

det Г 0,

а

Ap (s), g (s),

n,

R,

w_ и

/е_ со­

впадают с

соответствующими величинами

в

§

1

главы 2.

Исходя из структуры векторов щ0 и ^оі закон управления

(2.58)

можно переписать

в виде

 

 

 

 

 

 

 

где

 

 

 

и = wxx +

wyy,

 

 

 

 

 

(2.30)

 

 

 

 

 

—1г

 

 

 

 

 

 

 

w.

s2) Д р (s)

 

 

 

 

 

 

 

 

--- г-— — - -fV-tn

ü-_p ■

 

 

 

 

 

 

 

 

g* (s)

 

'

 

—1

g * iß)

 

 

 

(2.31)

 

(c —

s2) Д р (s)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[6_Г0 — Y (S)O_],

 

 

 

 

 

 

 

 

 

 

 

 

g* (5)

T. e. представить управление и в виде суммы «управлений по отклонению» (wxx) и «управлений по возмущению» (wyy).

Положив в первых п уравнениях системы (2.29) = = |2 = ••■ = = 0, можно из этих уравнений выразить у через X и и и, таким образом, исключить из закона управле­ ния (2.30) «управления по возмущению»:

у = W ('Рх ~ ти^ и = [у ® + Wy1n]~ '[у ® w* + WyP]'

Подставив значения wx и wy из (2.31) в последнее урав­ нение, получим

(с — s2) Д *

(s)

/е_Г

@<5

и = wx =

+

g*(s)

 

 

 

Ы Г ' Р

1

n*R х,

 

g*(.s)

 

т. е. передаточная функция w совпадает с передаточной функ­ цией оптимального регулятора при решении «-мерной зада­ чи (2.13).

5*

67

 

Проиллюстрируем изложенное на примере. Пусть воз­ мущенное движение объекта описывается системой линей­

ных дифференциальных

уравнений

X = X +

Зу

— и — Іѵ

У =

га — и — 6а,

— 71±Zi

Z-L

— Т]£а,

где §lf |2, £3 и |4 — некоррелированные стационарные слу­ чайные процессы типа «белый шум» с единичной спектраль­

ной

плотностью.

 

 

 

 

 

 

 

 

 

В

обозначениях,

принятых в

§ 1

гл. 2, имеем

 

 

- \ — p

3

 

— 1

 

 

0

 

 

~ 1 "

 

0

2 p

0

 

 

— 1

,

m =

1

p =

0

 

TJ J +

1

 

0

 

0

 

0

 

 

 

 

 

 

0

0

 

0

 

 

T ,P + 1_

_ 0 _

 

 

Г

=

diag {1,

 

1,

тъ

т2}.

 

 

Необходимо определить закон управления

 

 

 

u = wx (p)x +

 

W y { p ) y +

 

wZl р )

2j

+ w2t (р ) z„

(2.32)

 

 

 

 

 

 

 

(

 

 

 

т. е. передаточные функции регулятора wx (р) =

wx (p)/wQ(р),

wy ( р )

— Wy (p)/w0 { р ) ,

 

wZl ( р )

=

w2, (p)/w0 (p),

w2t

(p) =

= w2j {p)/w0 (p) таким образом, чтобы при устойчивой замкнутой системе объект -f- регулятор минимизировалась

величина е = с (и2) +

(а 2).

 

 

 

Решением, согласно (2.13), будут функции

 

и»о (s) =

s + Ѵ~с +

6,

 

 

 

(s) =

3 (]/~c +

1),

 

 

 

w,,(s) =

- 9 ( ] / 3

+

1) + 12 (V~c +

2),

(2.33)

W2 M

3(/c+l)

r

 

 

 

w ----------- T\ +

l

1 x'

 

 

®z, (s) =

9 (Vc+ 1)

Т»

12 ( / c

2) rn

 

 

Г . +

І

 

 

2 Г ,+ 1

'

 

63

а характеристический определитель замкнутой системы оп­ ределится, согласно (2.28), формулой

А (s) * 8 (®) =

= ( / Ь + s) (1 + S) (2 + S) (1 + 7V0 (1 + T2s). (2.34)

Если исходную систему переписать в виде

Г Р х

- Б

;

1___

I---- О

где

РX

' X ' ~тх

Z . 0 . и -J- Л г .

ТіР + 1

О -

. О

Т2р + \ _

Е — единичная матрица,

0 — нулевая матрица

размера

2 X 2, о — нулевой вектор-столбец, х0 [х, у]', z =

[zlt z2\,

mx = [1,

1]',

g, =

lg,, gj]',

% = ІТі Із, х2Ы' и положить & =

= |2 =

0, то

эта

система

запишется так:

 

 

 

 

Рхх0г — пгхи,

(2.35)

 

 

 

P j = h ,

 

т. е. с точностью до обозначений совпадает с системой, опи­ сывающей возмущеннее движение объекта в примере, рас­ смотренном в § 1 гл. 2.

Записав выражение (2.32) в виде

 

и = wx<tx + w^,

(2.36)

где

 

= к ( р ) . wy ( р) ] . wz = [a»4l ( р ) ,

си л , (р)],

и подставив в него значение вектора z из первого уравнения (2.35), получим и = wx0, где

со = (1 + a y n j- ' (wXo+ wzP J .

(2.37)

Нетрудно убедиться, что передаточные функции регу­ лятора, определяемые формулами (2.37) и (2.33), совпадают с передаточными функциями регулятора, полученными при решении примера § 1 гл. 2.

Соотношение (2.36) вскрывает структуру оптимального регулятора: при S,-; (со) =£ const управляющее воздействие является суммой «управлений по отклонению» (wXtx0) и «управлений по возмущению» (wzz), и поэтому для данного класса задач только частично применим весьма эффективный метод определения параметров оптимального регулятора,

69

Соседние файлы в папке книги из ГПНТБ