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

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

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

к*

M (s) M

___Ф, (S) — <С— S2) Р (— s) ^ 1’ ^

«

X

(— s)

1 W

^

S > *>

S'M (s) M(— s) j

*

 

 

 

 

 

y°o

 

 

 

X 6Ф (— s)ds +

Y

^ j[/7W(s) M (— s) +

 

 

 

 

 

 

— / 0 0

 

 

 

+ (C-

S2) P ( S ) P ( - S)]Z I^

Z ? 0 1( - S) -

 

 

-

(c -

s2>я (s) л т д а т У

бфі (s) *•

 

В соответствии с методикой решения уравнения Вине­ ра — Хопфа [24, 30] выполним следующие промежуточные преобразования в подынтегральном выражении.

Множитель, стоящий перед Ф , (s), представим в виде произведения двух функций

[гМ (s) M ( - s ) + ( c - s2) Р (S) Р (-S )] M{s^ _ s) =

= Dj (s) D, (— s),

причем Dx (s) имеет нули и полюсы только в левой полу­ плоскости.

Выражение -

s2) Р ( s ) M{s)^

_

s) ■

предста­

вим в виде суммы трех слагаемых

 

 

 

___s2) р

(____s)

^

.

1

_

 

^ м

(s) М (— s)

 

D.i—s)

~

— В 10(s) +

^ 1+ (S) +

^ l—(s)>

 

где B 10 (s) — целая часть (полином от s), В i+ (s) — правиль­ ная дробь с полюсами только в левой полуплоскости s, Bi_ (s) — правильная дробь с полюсами только в правой полуплоскости S.

В принятых обозначениях выражение для первой вариа­ ции функционала запишется в виде

/со

бе = у [ D1( - s ) [D 1(s)<b1( s ) - B u ( s ) - B i + ( s ) -

*— /оо

/со

— Bi—(s)]6® i(— s) d s -j-г-

§

D1(s)[D1( s)Oj.( s)

'

— jeо

 

- ß 10( - s) - B l+ ( -

s) -

5 ,_ (s)] SO, (s) ds.

10

Тогда функция Ф , (s), обращающая в нуль первую вариа­ цию функционала бе и имеющая полюсы только в левой полуплоскости s, определится следующим соотношением:

фі (s) =

[Вы (s) + B l+ (s)],

а минимум функционала е

 

/со

 

®min

ß ,_ (s )ß ,_ (- ■s)+

--JOG

 

г{с-

ds.

+ гМ(s) M(—s) -|- (c — s2) P (s) P (—s)

Из (1.5) получим выражение для искомой передаточной функции

iw /,л _

Р (s) ®і (5)

1

WW

М (s) Ф, (s)

'

Определим, при каких условиях выбор варьируемой функции в виде (1.7) обеспечивает устойчивость замкнутой

системы объект + регулятор.

 

 

Представив Ф , (s) как отношение двух полиномов

cpu (s)

и ф10 (s), т. е.

 

 

 

 

0,(5) = Фа 60

 

 

Фіо (s)

 

получим

 

 

 

т /ч\

Фа (s) Р (5)

Фіо (s)

( 1. 8)

U

фи (s) Af (s)

 

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

Р (р) X — М (р) и = ф,

[фю(р) — Фп (Р) р (Р)1 X+ Фи (Р) М (р)и = 0. Характеристический определитель замкнутой системы

P ( S )

— M(s)

Д(5) =

Фіо (sW s ) .

Фіо (s) — Фи (s)

(s) Фи (S) Л4 (s)

Нули полинома cp10 (s), являющиеся полюсами функции Ф, (s), расположены в левой полуплоскости s, так как полу­ чаемая в результате решения уравнения Винера—Хопфа

И

передаточная функция Ф х (s) физически реализуема (имеет полюсы только в левой полуплоскости s). Таким образом, замкнутая система объект + регулятор будет устойчивой, если все нули полинома М (s) находятся в левой полупло­ скости. В противном случае (полином М (s) имеет хотя бы один нуль с положительной вещественной частью) замкну­ тая система будет неустойчива.

Рассмотрим другой вариант выбора варьируемой функ­ ции.

Пусть

 

 

Ф 2 (S) =

р(S) _ yvl (S) W (s) •

(1-9)

Тогда

передаточные функции

Fx (s), Fu (s) и

функционал

(1.6)

имеют

следующий

вид:

 

 

 

 

Fx(s) = - p ^ l® 2 (s )M (s )+ l],

 

 

/DO

Fu(s) =

Ф а(«)»

 

 

 

 

 

 

e =

| [

{[гМ (s) M ( -

s) +

(с - sa) P (s) p ( _ s)] Ф 2 (s) X

X ® s ( - s ) + rM (s)® 2(s) + /-A4 (—s) a>a(—s) + r} p ^ p l _ s) ds.

Ход решения задачи аналогичен ходу в первом случае, поэтому приведем только окончательную формулу

Ф2 (S) = ----- [ß20 (S) + В 2+ (S)],

где D2 (S ) — функция, имеющая нули и полюсы только в ле­ вой полуплоскости, которая определяется из разложения

[гМ (S) М ( - S) + - S2) Р (S) Р ( - S)]

s) =

— ^2 (S) ^2 ( S)>

 

ß 20 (s) и ß 2+ (s) — целая часть (полином от s) и правиль­ ная дробь, имеющая полюсы только в левой полуплоскости s, в разложении

гМ ( - s) S^(s)

£>., (-S) Р (s) Р (-S) = В 2 0 (S) + В 2 + (s) + В2 —(s)

(ß2- (s) — правильная дробь, имеющая полюсы только в правой полуплоскости s). Минимум функционала в этом

12

случае

е,,,і" — / I/оо ' ^ 2 —($) В2 —(— s) -f-

г (с — s2) (s)

ds.

^гМ (s) М (— s ) + (с — s~) Р ( s ) Р ( — s )

Передаточная функция оптимального регулятора, со­ гласно (1.5), определится формулой

 

W(s) =

ф 2 (S ) Р (S )

 

 

 

1 +

Ф 2 ( s ) М ( s )

 

 

 

 

 

 

 

 

и, если представить Ф 2 (s) в виде отношения двух

полино­

мов ср21

(s) и ф20 (s), т. е.

 

 

 

 

 

 

 

0

2(s)

Фаг ( s )

 

 

 

Фао (S) ’

 

 

 

 

 

 

 

 

 

Т О

 

 

 

 

 

 

 

 

 

W(s) =

 

:

Ф м

(S) Р (Я)

( s ) '

( 1 ■1 ° )

 

 

 

 

фае (s ) +

Ф г і ( s ) TW

Характеристический

 

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

замкнутой

системы

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

 

 

 

 

 

 

 

A(s) =

P(s)

 

 

 

~

M(S)

-

П Р М

 

- c p 21(s)P(s)

 

cp20(s) +

cpai (S)M(S)

(Sj'

Замкнутая система объект -+- регулятор будет устой­ чивой при выборе варьируемой функции в виде (1.9),если полином Р (s) имеет нули только в левой полуплоскости, т. е. объект устойчив. В случае неустойчивого объекта ре­ шение задачи в виде (1.10) приводит к неустойчивости замк­ нутой системы.

Таким образом, для объекта, передаточная функция которого имеет все нули, или все полюсы в левой полупло­ скости s, решение задачи (передаточная функция оптималь­ ного регулятора W (s)) определяется одной из формул (1.8), (1.10). В противном случае регуляторы, определяемые фор­ мулами (1.8), (1.10), приводят к неустойчивой замкнутой системе.

Как видно, решение поставленной задачи (нахождение минимума функционала качества на классе передаточных функций регулятора, обеспечивающих устойчивость замк­ нутой системы объект + регулятор) тесно связано с вы­ бором варьируемой функции.

13

Аналогичный вывод делает Ш. Л. С. Чанг в своей моно­ графии 130]. Он, в частности, рекомендует выбирать раз­

личные варьируемые функции для

минимально-фазовых

и неминимально-фазовых объектов.

Для неустойчивых

объектов, согласно [24, 29, 30], решение задачи разбивается на два этапа: стабилизация объекта и нахождение уравне­ ния регулятора из условия минимума функционала.

В [14] показано, что, усложнив вид варьируемой функ­ ции, можно получить решение задачи для объектов с разны­ ми динамическими характеристиками (произвольным рас­ положением нулей и полюсов передаточной функции не­ управляемого объекта).

Действительно, если варьируемую функцию выбрать в

виде линейной комбинации функций Ф х (s) и Ф 2 (s):

 

 

в

и

и

Ф , в + Р в

я > .(» )= ; g

i ^ (g ,y g )

о - " »

где

а (s)

и

ß (s) — пока

произвольные

полиномы

от

s, то

можно получить «универсальное» решение поставленной за­ дачи.

Передаточные функции Fx (s), Fu (s) (1.5) и функционал

(1.6) имеют следующий

вид:

 

р /<Л

Ф (s) М(s) + р (s)

 

 

 

ß (s) P (s) + а (s)M(s)

p

. _

Ф (s) P(s) — a (s)

 

“ i

;

ß (s) P (s) + а (s) M(s)

/©O

 

 

 

e = T J (S ) P (s) +

а (s) M (S )] Iß ( - s) P ( - S) 4 • а ( - s) /VI ( - s)] 83

—JOO

 

 

 

X{[rM (s) M (s) + (c — s2) P (s) P (— s)] Ф (s)Ф (— s) -J-

+[rß (— s) M (s) — (c — s2) a ( —s)P (s)] Ф (s) +

+ [/-ß(s) M (— s) — (c — s2) а (s) P (— s)] Ф (— s) +

 

-f rß(s) ß (— s) -f- (c — s2) а (s) а (— s)\ds.

(1.12)

Последующие выкладки аналогичны двум первым случаям. Функция Ф (s), обращающая в нуль первую вариацию

функционала (1.12)

/оо

& = Т І

—JCO

D (s)D (-s)

 

r(s)

5Ф (— s) ds -j-

Q (s) Q( - s) Ф (5) + Q (s) Q (—s) .

D(s)D(-s)

Ф ( - 5 ) +

H - S )

 

6Ф (s) ds

Q(s) Q(— s)

 

Q (s) Q (— s)

 

14

и имеющая полюсы только в левой полуплоскости s, опреде­ лится соотношением

ф (S) = -

[ß o (s) + В+ (s)],

(1.13)

где D (s) — функция, имеющая нули и полюсы только в ле­ вой полуплоскости, определяется из разложения

[гМ (s) М (— s) +

(с — s2) Р (s) Р (— s)] 5,1, (s) =

 

=

D(s) D (— s),

(1.14)

Q (s) — полином с нулями только в левой полуплоскости, получившийся в результате разложения

(s) P(s) + а (s) М (s)] [ß ( - s) P ( - s ) +

 

+ a ( - s ) M ( - s ) } =

Q (s)Q (-s),

(1.15)

Т (s) = [rß (s) M (— s) — (c — s2)

а (s) P ( -

s)] S * (s),

(1.16)

B 0 (s) — целая часть (полином от s),

(s) — правиль­

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

 

T (s )

(s) +

(s) + B—(s)i

(1.17)

D

 

=

( — s )Q ( s )

 

 

 

 

ß _ (s) —

правильная дробь,

имеющая

полюсы

только

в правой полуплоскости s.

 

 

 

 

Минимальное

значение

функционала запишется

в виде

 

 

/со

 

 

 

 

 

 

'"= Т I

B _ (s )ß _ ( - s ) +

 

 

 

 

 

 

 

 

 

 

r ( c — s ^ S ^ l s )

1

(1.18)

+

г / И

(S) М ( — s ) +

(с — S - ) Р ( s ) Р ( — s)

ÜS'

 

Передаточная функция оптимального регулятора, согласно (1.5) и (1.12), определится формулой

W СО =

ф ( s ) р

~ а ( s )

(1.19)

 

ß(s) +

Q (s)M (s)

 

и, если Ф (s) представить отношением двух ПОЛИНОМОВ фх (s) и фо (s)

®(s) = Фі (s)

<Po(s) ’

то W (s) примет вид

iw СО =

Фі (5) Р (s) — Уо (s)« (s)

(1.20)

W

Фо (s) Р (s) + Фх (s) ЛГ (в)*

 

15

Характеристический

определитель замкнутой системы

имеет вид

 

 

— Фі (s) Р (s) +

Фо (s) а (s)

фо (s) ß (s) + Фі (s) M (s)

= Фо (s) [ß (s) P (s) +

а (s) M (s)].

Так как полином Фо (s) имеет нули только в левой полу­ плоскости, то для устойчивости замкнутой системы необхо­ димо, чтобы полином ß (s) Р (s) + а (s)M (s) имел нули толь­ ко в левой полуплоскости, т. е., согласно (1.15),

ß(s)P(s) + a(s)M (s) = Q(s).

(1.21)

При заданных Р (s) и М (s) соответствующим выбором полиномов а (s) и ß (s) можно добиться того, чтобы полином Q (s) имел нули только в левой полуплоскости. Другими словами, задавшись произвольным полиномом Q (s), нули которого лежат в левой полуплоскости, соотношение (1.21) можно рассматривать как полиномиальное уравнение отно­ сительно а (s) и ß (s) и решать это уравнение методами, ука­ занными, например, в [6]. Таким образом замкнутая система объект + регулятор будет устойчива при любых полиномах Р (s) и М (s), характеризующих свойства объекта. Остается показать, что произвол в выборе полинома Q (s), а следова­ тельно, и полиномов а (s) и ß (s), не отражается на конеч­ ном результате, т. е. передаточная функция оптимального регулятора и минимум функционала качества не зависят

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

 

Введем некоторые добавочные обозначения.

Пусть

 

 

R (s),

(s) Ni (— s) — (с — s2) а (s) P (— s) =

rM (s) M (— s) +

(c — s2) P (s) P (— s) =

G(s) G(— s),

 

 

(— s)

( 1.22)

Sy (s) =

П (S) Г ,

 

Г. Й r 0 (—S)

 

 

 

где G (s), Г 0 (s) и I\ (s) — полиномы с нулями только

в ле­

вой полуплоскости.

 

Тогда, согласно (1.14), (1.16) и (1.17):

 

D(s) = G (s)-Гi і й .

(1.23)

о й

 

 

(1.24)

В о (s) + В + (s) + В -

(s) = G ( — s ) Q й Г о й '

(1.25)

 

16

Если функцию Ф (s), определяемую выражением П.13), переписать в виде

 

 

Ф (5) = -

T ( s )

+

Q (s) 5 _ ( s ) ,

 

 

 

 

 

 

 

 

 

 

 

D (s) D ( - s)

D ( s )

 

 

 

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

 

 

 

 

 

 

 

 

Ф (5) =

G ( s ) G ( - s )

 

Q ( S ) r„(s) ß _ ( S).

 

(1.26)

 

 

 

 

 

G ( s ) T \ (s)

 

 

 

Нетрудно проверить, что

 

 

 

 

 

 

/И (s)

(s) =

ß (s) G (s) G (— s) — (c — s2) P (— s) Q (s),

(1.27)

P (s) R (s) = — а (s) G (s) G (— s) +

rM (— s) Q (s).

(1.28)

Тогда,

согласно

(1.19), (1.26),

(1.27) и

(1.28),

получим

w , . _

r M (—

s) Г; (s) G ~ * (—

s) + P

(s) Г „

(s) B _ (s)

„ 2 g .

U

 

( c - s * ) P ( - s ) r 1 ( s ) f f - , (-r s) + i M ( s ) r 0 ( s ) Ä _ ( s ) - ^ '

Рассмотрим

более детально

 

процедуру определения

ß _ (S ) .

Е с л и умножить и разделить

выражение

(1.25) на

Р (s), то, учитывая (1.28), получим

 

 

 

 

В0 (s) +

В+ (s) -|- В _

 

а (s) G (s)

rM (—

s)l

I\(s)

(s) =

 

О Й

+

G ( — s)

Г 0 (s)'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.30)

Аналогично, если умножить и разделить (1.25) на М (s), то, учитывая (1.27), получим

В0 (s) + В+ (s) -(- ß _ (s) =

1

р (s) G (s)

(c — s3) P (— s)

r t (s)

(1.31)

М (s)

Q(s)

G (— s)

r 0 (s)

Как видно из выражения (1.25), полюсами В_ (s) будут только нули полинома G (— s).

Отметим, что три полинома Р (s), TW(s) и G (— s) не могут иметь даже одного общего нуля, так как в этом случае объект регулирования будет неуправляем и неустойчив. Поэтому, определяя В _ (s) в виде суммы правильных дробей, полю­ сами которых являются нули полинома G (— s), при вычис­ лении слагаемых следует пользоваться выражением (1.30), если полюс дроби не совпадает с нулем Р (s), или (1.31), если полюс дроби совпадает с нулем Р (s). В каждом из этих случаев правильные дроби, полюсами которых являются ну­ ли полинома G (— s), появятся только вследствие наличия

2 3-582

17

второго слагаемого в (1.30) и (1.31), которые не содержат полиномов а (s) н ß (s).

Следовательно, функция ß _ (s) не зависит от полиномов а (s) II ß (s) и, согласно (1.18) и (1.29), минимум функциона­ ла качества и передаточная функция оптимального регуля­ тора W (s) также не будут зависеть от конкретного вида полиномов а (s) и ß (s), если Q (s) — гурвицев полином.

Можно доказать, что в числителе и знаменателе переда­ точной функции W (s), записанной в виде (1.29), стоят поли­ номы, т. е. дробные части соответствующих слагаемых равны по величине и противоположны по знаку. Ограничимся доказательством этого утверждения для числителя, когда

полином G (— s) не имеет

кратных нулей.

Дробную часть первого

слагаемого

в числителе W (s)

можно записать в виде

 

 

lim

(5 — g l-) гМ ( —

s) I \ (S)

 

G(-s)

 

где g; — нули полинома G (— s).

Согласно (1.30), дробная часть второго слагаемого опре­ деляется соотношением

V _____!____ Hm

( s - g i)P(s)r0(s)rM(-s)T1 (s)

s - B i

 

P ( s ) G ( - s ) r 0(s)

= y-

lim

( s — g l) rM {— s) r \ (s)

 

- Bi

G(-s)

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

Доказательство такого же утверждения для знаменате­ ля аналогично приведенному выше, только при определении W (s) вместо формулы (1.30) следует воспользоваться фор­ мулой (1.31).

Поскольку в числителе и знаменателе W (s) (1.29) стоят

полиномы, характеристический

определитель замкнутой

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

А 00 =

 

 

Р (s)

— М (s)

гМ (— s) Гх (s) G - '(- s ) - - S 2) Р ( -

s) r x(s) G~’(^ s) +

Р (s) Г0 (s) ß _ (s)

+

M (s) Г0 (s) ß _ (s)

= [rM (s) M (— s) + (c — s2) P (s) P (— s)] G 1(— s)I\(s),

18

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

A(s) = G(s) Г! (s),

т. е. А (s) также не зависит от полиномов Q (s), а (s) и ß (s).

З А Д А Ч А С И Н ТЕ З А П Р И О Д НО М В Н Е Ш Н Е М

В О З М У Щ ЕН И И

Изложенные в предыдущем параграфе варианты решения задачи свидетельствуют о том, что существует тесная связь между динамическими свойствами стабилизируемого объек­ та и видом вводимой варьируемой функции. Однако рас­ смотренные варианты лишь констатируют существующую связь, не вскрывая ее. Даже «универсальный» метод (с по­ линомами а (s) и ß (s)) не вскрывает этой связи, так как оста­ ется неясным смысл вводимых полиномов а (s) и ß (s), от которых, в сущности, не зависит решение. Поэтому пред­ ставляется целесообразным обсудить общий подход к реше­ нию задачи синтеза системы стабилизации, из которого рассмотренные в § 1 «рецепты» решения отдельных классов задач вытекали бы как частные случаи, и, с другой стороны, его можно было бы обобщить на многомерную задачу (с про­ извольным числом внешних возмущений и управляющих воздействий).

Ниже этот общий подход будет изложен применительно к простейшей задаче, когда движение объекта описывается дифференциальным уравнением (1.1). Естественно, что по­ лученные конечные результаты будут совпадать с результат тами § 1.

Как отмечалось, специфика рассматриваемой задачи за­ ключается в необходимости обеспечения устойчивости замк­ нутой системы объект -f- регулятор, т. е. передаточная функ­ ция регулятора разыскивается в классе функций, обеспечи­ вающих устойчивость замкнутой системы. Следовательно, при решении вариационной задачи (минимизации функцио­ нала (1.6)) требуется ограничить вариации 6117(5) так, чтобы соответствующие вариации 6Fx (s) и 8FU(s) функций

Fx (s) и Fu (s) имели бы

полюсы только в левой полупло­

скости.

 

 

Поскольку вариации 6Fx (s) и 8FU(s) зависимы (обуслов­

лены вариацией только

одной

функции 6U7(s)), то Fx (s)

и Fu (s) можно выразить

через

одну свободную варьируе­

мую функцию Ф (s).

 

 

2*

-Vблинная

19

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