Файл: Ларин В.Б. Синтез оптимальных линейных систем с обратной связью.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 13.07.2024

Просмотров: 128

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Покажем независимость решения задачи от произвола в выборе матрицы преобразования Z (полиномов а (s) и

ß(s))-

Пусть варьируемая функция Ф (s) выбрана в виде Ф (s) = а (s) Fx (s) + ß (s) Fu (s),

T. e. передаточные функции Fx (s) и Z7, (s) связаны с варьи­ руемой функцией матричным уравнением

'P (s)

— ZW(s)’

'Fxto'

' 1 '

а (s)

=

Z

ß (s) _ .F u(s).

f u (S).

. Ф (s)

Здесь а (s) и ß (s) не совпадают с а (s) и ß (s), но удовле­ творяют требованию аналитичности в правой полуплоскости

матрицы Z вместе с обратной:

det Z = а (s) М (s) + ß (s) Р (s) = Q (s)

(1.41)

(Q (s) — гурвицев полином).

Для доказательства независимости решения от выбора

полиномов а (s) и ß (s),

как видно из формул (1.37), (1.38),

достаточно доказать равенство

 

 

 

B -(s) =

B -(s).

 

 

Пусть П — матрица, определяемая уравнением

 

 

Z =

nZ,

 

 

(1.42)

т. е.

 

 

 

 

 

Р (s)

A4 (s)

P(s)

M (s) — 1

 

П = ZZ- ' =

 

 

а (s)

ß (s)

 

_“ (s)

ß(s)

 

 

 

 

 

 

 

1

 

0

 

g(s) ß (s) ß

(5) а

(s) у

Q (s)

(1.43)

 

 

Q ( s )

 

 

Q(s)J

 

где

 

 

 

 

 

Q(s) =

det Z,

Q (s) = det Z.

 

Как видно из (1.43), матрица П имеет структуру

 

П =

1

 

О

 

 

Пх (s) П, (s)

 

 

 

 

 

22


причем нули и полюсы П2 (s) и полюсы П1 (s) расположены

в левой полуплоскости, а полиномы а (s) и ß (s), учитывая (1.42), определятся уравнениями

а (s) =

Пх (5) Р (5) +

П2 (5) а (s),

ß (s) =

 

(1.44)

— Пх (s) М (s) -f- П2 (s) ß (s). j

Согласно (1.39)

и (1.40),

 

 

1

R ^ T ^ S )

B-(s) =

Q(s) Г 0(s)

 

G{~ S)

где [ 1_ — сумма всех правильных дробей с полюсами в правой полуплоскости, получающаяся в результате разло­ жения функции, стоящей в квадратных скобках, нацелую часть и правильные дроби

R (s) =

rß (s) М (— s) — (с — s2) а (s) P ( — s).

Учитывая (1.44)

и (1.41), получаем

R (s) =

-

Пі (s) G (s) G ( - s) + П2 (s) Я (s),

Q (s) = П2 (s) Q (s),

,

{ -

П , (5) G(s) G ( - s) + n 2 (s) fl (s)) Г , (s) 1

W

 

G ( - s) П , (s) Q (s) Г 0 (s)

R(s) Гг (в)

G (— s) Q (s) Г о (s)

что и требовалось доказать.

Приведем примеры оптимальной стабилизации неустой­ чивых объектов минимальными энергетическими затратами.

I.

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

ным уравнением

 

 

 

 

р"~х +

рх 2х — ри -f и -J- ф,

Р ~ - jf,

(со) =

ь., _І_ м2 .

Требуется определить уравнение регулятора

(1.45)

 

 

W0(p)u =

W (p)x

 

(1.46)

(операторные полиномы W0 (р)

и W (р)) так,

чтобы на

классе

устойчивых замкнутых

систем

минимизировался

23


функционал

е= с2(и2) + (и2).

Впринятых обозначениях:

P{s) = s2 +

s — 2 =

(s — 1) (s +

2),

/W (s)=s + 1, S*(s) =

bi-A- r ,

т. e. 1^ =

1, Г0 = Ь + s,

и задача заключается в отыскании передаточной функции регулятора в цепи обратной связи

W (s) = І ^ 1(s) W (s).

(1.47)

Для того чтобы полином Q (s) =

ß (s) Р (s) -f- а (s) M (s)

не имел нулей в правой полуплоскости s, выбираем ß (s) =

= 0, а (s) =

1, т. е. Q (s)

= М (s) = s 4- 1.

 

Согласно (1.40),

(1.39),

G(s)G(— s) =

(c + s)(c — s)(2 +

+- s) (2 — s)(l

+ s ) ( l

s)

(г =

0),т. e.

 

 

 

G(s) = ( l + s ) ( 2

+ s)(c +

s),

 

R (s) = (c -f s) (c — s) (2 — s) (1 -I- s), В (s)

=

(c + s)

(6 + s ) ( l - s )

 

 

 

 

 

 

 

T. e.

ß_ (s) =

( c + О

 

 

 

( 6 + l)(l — s)

 

 

 

 

 

 

и передаточная функция оптимального регулятора, согласно

(1.37),

определится формулой

 

W(s)

=

( s - l ) ( s + 2 )0 + s ) ( c + 1)

 

(с + s) (с — s) (— I — s) (2 — s)

 

 

 

 

(ö + 1) (1 — s)

(c — s)(l — s)(2 — s)

 

 

 

 

 

 

 

 

(1 + s)(6 + s)(c + 1)

1

 

 

^

(è + 1) (1 — s)

J

 

 

( c + 1 )

(2 + s)(ft + s)

(1.48)

 

 

(c — b)

(1 + s)

 

 

 

T . e. оптимальный закон управления соответствует уравне­ нию

(с — Ь) {и + и) = (с + 1) + Ф + 2) X + 2Ьх]. (1.49)

Характеристический определитель замкнутой системы (уравнения (1.45) и (1.49)) имеет вид

Р (s) — М (s)

Д(s) =

W(s) - W 0{s)

24


_ (s — 0 (s +

2)

(s + 1)

^

(s + 2) (b +

s) (c + 1)

(1 + s) (b c)

 

= (b + l ) ( l + s ) ( 2 + s)(c + s).

Заметим, что среди нулей характеристического опреде­ лителя оптимальной замкнутой системы находятся «левый» нуль (s = —2) и зеркальное отражение (s = — 1) относи­ тельно мнимой оси «правого» нуля (s = 1) полинома Р (s) (нули полинома Р (s) соответствуют полюсам передаточной функции объекта).

 

 

 

 

 

 

ГТ

 

 

 

 

 

 

 

II. Пусть

Р (s) =

s'1-[- 2

ciiSn~! — полином,

все

нули ко-

 

 

 

 

 

 

г=і

 

 

 

М — const,

торого расположены в правой полуплоскости,

S,i) (со) =

Х?|Ь минимизируемый функционал имеет тот же вид,

что и в

предыдущем

примере,

т. е.

 

 

 

 

 

 

 

 

 

е = с2 (и2) + (и2).

 

 

 

 

Полагая ß =

0, а =

получаем det Z =

Р$-\- Ма =

=

1, т. е.

Q (s)

=

1 — не

имеет

нулей в правой

полупло­

скости.

 

 

 

 

 

 

 

 

 

 

 

 

Подставим необходимые величины в (1.40)

и

(1.39):

 

 

 

G (s) G (— s) =

(с2 — s2) Р (s) Р (— s),

 

 

 

 

 

 

т.

е.

G (s) = +

s) Р (— s);

 

 

 

^(s) = - ^ ( c

z- s

2) P ( - S); B

( s ) = ~ h - ( c

+ s)p- ^ ,

 

т. е.

В+ (s) = 0, В0(s) =

( -

1)п+1 (s +

с -

2а,),

 

B_(s) = B ( s ) - B 0( s ) - B + (s) =

М

+

S) X

 

 

 

X

P ( - s)

( - l ) n+I(c +

s - 2 ßl)

 

 

 

 

 

P(s)

 

 

 

 

 

 

 

 

 

 

 

 

 

и,

согласно

(1.37),

 

 

 

 

 

 

 

 

 

W (s) =

----------- P(?)B-(s)

 

P(s)B-(s)

 

 

 

 

+

 

 

 

Mß-(s)

M [B -(s)— B(s)]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(s)B-(s)

 

 

 

 

 

 

 

 

 

 

[So (s) + B, (s)l M

 

 

 

 

 

_

(C+

S-2 fl1) P (s ) - (C+

s )( - l)'1P (- s )

 

(1.50)

 

 

■"

 

 

M (s -j- c — 2ax)

 

 

 

 

 

 

 

 

 

 

 

25


Характеристический определитель замкнутой системы объект + регулятор имеет вид

P{s) М

А (s) = (C + s - 2 fll) P ( S) - ( c + s) ( - l ) " P ( - s) —М (s -)- с — 2^)

= (— I)" -' (с + S)/> ( _ * ),

т. е. один нуль характеристического определителя опти­ мальной замкнутой системы (s = —с) обусловлен критерием качества, а остальные являются зеркальным отражением относительно мнимой оси «правых» нулей полинома Р (s) («правых» полюсов передаточной функции объекта).

III. Пусть Р (s) = s — а, остальные исходные данн совпадают с предыдущим примером. Подставив это значе­ ние Р (s) в (1.50), получим

vw /„ч _ _

(с ~Т s 4~ 2а) (s — а) — (с -f- s) (s а )

 

2а (а 4- с)

'

'

 

— /И (s +

с 4- 2а)

 

 

М(s 4- с +

2а)

Так

как В (s) = -----(с +

s)

 

 

 

 

 

V (с +

s) (— s — а)

TO ß_ (s) -

2іЛд, 4" a)

и,

согласно

 

( s - f l )

M (sа)

 

 

 

 

 

 

 

(1.38),

 

 

 

 

 

 

 

 

 

 

 

 

 

/со

 

 

 

4па (а +

с)2 Ц,

 

 

<?шіп = 4 - J В - (S) В - (— S) ds =

 

 

 

W

 

 

 

 

— /C O

 

 

 

 

 

 

 

 

 

 

§ 3.

С И Н ТЕ З

С И С ТЕМ

С ТА Б И Л И З А Ц И И

П Р И . П

 

 

 

 

В Н Е Ш Н И Х В О З М У Щ Е Н И Я Х И т У П Р А В Л Я Ю ­

 

 

 

 

Щ И Х В О З Д Е Й С ТВ И Я Х

 

 

 

Ранее на модели простейшей задачи были подробно иссле­ дованы особенности использования спектральных методов при определении дифференциального уравнения оптималь­ ного регулятора, включенного в цепь обратной связи. Не­ смотря на то, что конечные результаты § 2 всего лишь по- ѵ вторяют результаты § 1, идея ограничения вариаций матри-

цы р- предложенная в §2, сводящаяся к требованию ана-

г и I

литичности в правой полуплоскости матрицы Z (формула

26