книги из ГПНТБ / Ларин В.Б. Синтез оптимальных линейных систем с обратной связью
.pdfк* |
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), |
rß (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 |