а п - 3 ап-і а п - ь |
|
Я л - 1 |
а п — 1 |
a n - 3 > 0 . |
(4.20) |
0 |
an |
an-i |
|
Если хотя бы один из миноров |
|
|
|
Л , < 0 |
|
(4.21) |
— система неустойчива. |
|
|
|
Если хотя бы один из миноров |
|
|
|
Д, — 0 |
|
(4.22) |
— система на границе устойчивости.
Следует отметить, что необходимым условием положительно сти детерминантов Д и Д(. является положительность коэффи циентов at характеристического полинома, т. е. необходимым ус ловием критерия Гурвица является:
ап > 0, а„_, > 0, . . а, > 0.........аг > 0., а0> 0. (4.23)
Действительно, выполнение условий критерия Гурвица для си
стемы с передаточной функцией Ф (р) гарантирует отри-
А(р)
цательность действительных частей корней ее характеристиче ского уравнения:
А{р) = апрп+ а п_хр»-'+ . . . + а 1(о + а 0 =
= а„ір - p J i p —рз). . • (р — Рі) ■■■ІР-Рп).
Учитывая, что корни p t могут иметь вид
Р |
і = ± |
«;> |
аі > |
0, |
или |
= |
± *і |
± УР/. |
|
Р і |
|
Рі+1 = |
і |
аі + УР,-> |
0, |
из (4.24) находим, что если все Re [рг] |
= ± а1 < 0, то |
будет представляться произведением сомножителей вида:
(4.24)
(4.25)
(4.26)
А (р)
(Р ~Рі) = \Р - ( ~ а/)1 =(/» + |
«/). |
|
(4.27) |
или |
|
|
|
( р - Р і ) ( р — Р і + і ) - [ р - ( - « і + У Р , ) І ( р - ( - * / — y ? / )J = |
== [ ( ^ + “ і) + У Р і][( p + * t ) — УРі I = [(/> + |
в/)2 + |
P/2] |
= |
= [p' + 2atp + ( a * + $*)}, |
|
|
(4.28) |
которые при перемножении дают полином |
А (р) |
= |
а„ р п + |
+ . . . -\-а[р-\-а0 с положительными коэффициентами.
Таким образом, для исследования устойчивости линейной стационарной системы с помощью критерия Гурвица необходимо:
а) определить передаточную функцию замкнутой системы
В(o')
Ф(р)= — — и ее характеристический полином в виде:
А(Р)
|
А(Р) = апРп |
ап_лрп~' + . |
агрг+ . . . 4- ах р + а 0; |
|
б) проверить положительность коэффициентов характеристи |
|
ческого полинома: |
ап, а„_ ,,. . ah . . ., аи а0. Если хотя бы один |
|
из коэффициентов |
at < |
0, система не может быть устойчивой; |
|
в) если |
все |
коэффициенты а1 > 0, |
составить главный детер |
|
минант Гурвица |
Д |
и детерминанты |
и проверить их положи |
|
тельность. |
J |
А>0, |
|
|
|
Если |
|
|
|
I |
> 0 |
— система устойчива. Если хотя бы один |
|
Д; < 0 |
|
— система неустойчива. Если хотя бы один Д4- =■ 0 — |
|
система на границе устойчивости. |
|
2.Частные случаи условий критерия Гурвииа
1.Для системы первого порядка [А(р) = аф + а0].
Для системы первого порядка условия |
(4.19) — (4.22) крите |
рия Гурвица принимают вид: |
|
|
Система устойчива, если |
|
|
ß i > 0 |
, |
|
а 0 > 0 |
. |
(4.29) |
Система неустойчива, если |
|
|
ß l > 0 |
, |
|
(4.30)
Й О < 0 .
Система на границе устойчивости, если
ß i > 0 ,
(4.31)
ß o = 0 ,
т. е. если А( р ) —а\р.
Действительно, для системы первого порядка
ап— аі>
а„—і — а0,
поэтому А = I в„1.
2.Для системы второго порядка [А(р) = а2р2-\~а\Р+ ао]-
Для системы второго порядка условия |
(4.19) —(4.22) крите |
рия Гурвица принимают вид: |
|
Система устойчива, если |
|
0, |
(4.32) |
аі>0. |
а0>0. |
|
Система неустойчива, если |
|
|
|
а2>0, |
|
а2> 0, |
а2>0, |
(4.33) |
а ,< 0 , |
или |
Оі>0, |
или |
аі<0, |
ай> 0, |
|
a0< ° . |
а0<0- |
|
Система на границе устойчивости, если |
|
|
а2> 0, |
|
а2> 0, |
а2>0, |
(4.34) |
о, = 0, |
или |
а,>0, |
или |
а, = 0, |
о0=0, |
|
а0= 0, |
а0= 0. |
|
Действительно, для системы второго порядка |
|
|
|
ап — а21 |
|
|
|
|
а п - 1 = |
а ѵ |
|
|
поэтому |
|
ап-2 ~ аЪ1 |
|
|
|
|
|
|
|
д = ао 0 = |
anö i > 0 , |
если |
> 0, ап> |
0; |
а(Хл |
|
|
|
|
Aj — I |
I — О-i |
0. |
|
|
|
3. Для системы третьего порядка [А(р)= а3р3-\~а2р2-\-аір-\-а0].
Для системы третьего порядка условия (4.19) —(4.22) кри терия Гурвица принимают вид:
Система устойчива, если
а) а3 > 0,
а2> 0 , ах> 0, а0> 0;
б) аг а2> а0 аг. |
|
(4.35) |
Система неустойчива, если при а3>0: |
|
а, < 0, или |
ßjCO, или а 0< 0 , |
или |
аха2< а0а3 |
(4.36) |
Система на границе устойчивости, если при |
а; > 0: |
а 1а2 = |
а0а3. |
(4.37) |
Действительно, для системы третьего порядка
ап= аѵ ап- 1—а2->
ап-2 = аи
О-п—%==а0<
|
поэтому |
O.Q |
0 |
0 |
|
д= |
|
а2 |
|
“о :а 0Д2> 0 , если а 0> 0 , Д2>0'-г |
|
|
0 |
аг |
а2 |
|
Ді = |
|
= (%г |
> |
|
|
1+И = |
|
* 2 = |
= |
ajßj —а0а3 0, если + > 0 и аха2> а0а 3 |
3.Примеры анализа устойчивости систем
спомощью критерия Гурвица
П р и м е р 1.
С помощью критерия Гурвица исследовать устойчивость сис темы с заданной структурной схемой (рис. 4.3).
10 1 т
рp[Cl5pt1)
Р и с . 4.3. С труктурная схем а системы
Р е ше н и е .
1.Определение передаточной функции замкнутой системы и
еехарактеристического полинома:
а) определение передаточной функции разомкнутой системы:
W (p )- —1 |
20 |
|
р / 7 ( 0 , 5 / ? + ! ) р 2( 0 ,5 /7 + 1) |
б) определение передаточной функции замкнутой системы:
|
Ф(/?)= |
W(P) |
|
20 |
20 |
|
1 + \Ѵ(р) |
/?2(0,5/? + 1 ) + 20 |
0,5/?3+/?2+20 ’ |
|
|
|
в) определение характеристического полинома: |
|
А(р) = 0,5/734-/?2 + 20 = |
+/73 + a 2p 2 + aJp + aBl |
|
где |
|
а3 = 0,5; |
|
|
|
|
|
|
|
|
а 2= і ; |
|
|
|
|
а, = |
0; |
|
|
|
|
а0= |
20. |
|
2.Проверка положительности коэффициентов:
а3= 0,5 > 0;
аг = 1 > 0 ;
а \ = 0 ;
а0 = 20 > 0 .
Система неустойчива, так как аіа2<Ройз.
П р и м е р 2 .
С помощью критерия Гурвица исследовать устойчивость си стемы с заданной структурной схемой (рис. 4.4).
Р и с. 4.4. Структурная схема системы
Ре ш е н и е .
1.Определение передаточной функции замкнутой системы и
еехарактеристического полинома:
а) определение передаточной функции разомкнутой системы:
|
|
|
4 |
|
|
W(p)--= |
2 |
|
2р -f- |
1_________ |
\ |
Р + 1 |
|
+ |
|
/>(/>+ 1 ) |
|
1 |
4 |
|
|
|
2р + 1 |
|
2_________ 4 |
|
________ 1ф5________. |
Р ( Р + I)2 2/7 + |
1 |
+ 4 |
р ( р |
+ 1 )* (0,4/7 + 1 ) ’ |
б) определение передаточной функции замкнутой системы:
Ф ( р )= — |
----------------- . |
1 + W (р) |
р (/7+ І)2(0,4/7 + 1)+ 1,6 |
в) определение характеристического полинома:
А(р) = р ( р + 1 )2 (0 ,4 р + 1 ) 4 - 1,6 =
= |
0,4/7* + |
1,8/?3 + 2,4/72 + р + |
1,6 = |
= |
а4р 4 + |
а3 /73 + я2р 2 + а ,/7+ |
а0. |
2. Проверка положительности коэффициентов at\
|
аА= |
0,4 > |
0; |
|
а3= 1,8 >0; |
|
а2 = |
2,4 > 0 ; |
|
Ö! = |
1 |
> |
0 ; |
|
а0 = |
1,6 + |
0 . |
3. |
Проверка положительности |
А, |
и А: |
|
а0 |
0 |
|
0 |
0 |
а) |
А = &2 ß, ß0 0 . а о д з> |
|
0 |
«3 |
|
ß2 |
<>i |
|
0 |
|
«4 |
|
I г< |
сО |
ОО ГН |
|
I Q |
|
ІДо — CLo |
а 3 |
|
а4 |
|
ßj |
ßg |
д з = |
|
ß4 |
|
0 |
Л |
0 |
; |
|
|
|
|
|
— ß2ß3 —аА |
CN I |
|
0 |
|
= ß 1A2- ß 0ßg2= |
3 ,9 2 -1 ,6 .1,82= |
|
|
ßg |
|
|
= —1,264<0. |
Система неустойчива, так как Д3 = — 1,264 + 0 .
П р и м е р 3.
С помощью критерия Гурвица определить критическое и до пустимые (с точки зрения устойчивости) значения общего коэф фициента усиления разомкнутой системы, заданной структурной схемой рис. 4.5.
|
|
Р и с. 4.5. |
Структурная |
схема |
системы |
|
Р е ш е н и е . |
|
|
|
|
|
|
1 . |
Определение передаточной функции замкнутой системы и |
ее характеристического полинома: |
|
|
|
а) |
определение передаточной функции разомкнутой системы: |
|
|
(0,5р + 1 ) |
|
|
|
|
/р\ _ £ . |
(2р + |
1)________ k3 |
_ |
/+0,5/7 + l)ft3 |
= |
|
1 і |
, (0,5/7 |
+ |
1) р ( р + |
1) |
р(р+\)[(2р + 1)+ |
|
|
|
(2/? + |
1) |
|
+ (0,5/7+ 1)] |
|
_ ki k3(0 ,5 /7 + 1 ) |
_ |
/Sj^2 (0,5jt? —f- 1 ) |
/e(0,5 / 7 + 1 ) |
/7(/7 + l)(2,5/7 + |
2) |
2p(p+\)( 1,25/7+1) |
/7(p+l)(1,25/7 + 1) ’ |
где k = 0,5kik3 — общий коэффициент усиления разомкнутой системы*;
б) определение передаточной функции замкнутой системы:
ф ( X |
W(P) |
|
|
^ (0,5/7 + |
1)_________ . |
|
\ + W[p) |
/?(/7 + 1 ) ( 1 , 2 5 / 7 |
+ 1) + |
А ( 0 ,5 /7 + 1 ) ’ |
в) определение характеристического полинома: |
А (Р) = Р (/7+1)(1,25/7 + |
1) + к (0,5/7 + |
1) = |
= |
1,25/73 + |
2,25/?2 + |
(1 + |
0,5А) р + |
k = |
|
= + /73 + «2 і» 2 + + Р + +• |
|
|
2. Определение условий положительности |
коэффициентов at |
и условия а1а2>а0а3 (так как система третьего порядка): |
|
|
а3= |
1,25 > |
0; |
|
|
|
|
|
а2 = |
2,25 > |
0; |
|
|
|
|
+ = 1 + 0,5/> > 0, |
если |
& > — 2; |
|
а0= k > 0 , |
|
если |
к > 0 ; |
|
a ja 2 > a 0a3, |
|
если |
(1 +0,5&)2,25> \,25k. |
3. Определение критического и допустимых значений к:
а) критическое значение k= ккр коэффициента усиления со ответствует условию состояния системы на границе устойчивос ти, что будет иметь место при:
к ::==кКрі == |
2 ( + 5=2 О), |
к = кКр2 = |
0 |
(я0 = 0); |
к ~ кКРз = |
0,948 (+ + = я0а3). |
Так как отрицательные и нулевые значения |
коэффициента |
усиления не имеют физического смысла, то |
|
|
|
|
&кр = |
^крз “ |
0,948; |
|
|
|
б) допустимые значения k= /глоп |
коэффициента усиления со |
ответствуют устойчивой |
системе, |
что |
будет иметь |
место |
при |
|
k > |
- |
2 |
\ |
, |
|
п. |
|
|
|
|
k > 0 |
|
} |
лоп> |
’ |
|
|
|
(1 + |
0,5А) 2,25 > |
1,25k |
/гдоп < |
/гкр = 0,948. |
|
|
* k — общ ий |
коэф ф ициент |
усиления |
разом кнутой системы есть общ ий ко |
эф ф ициент передаточной функции |
разом кнутой |
системы W(p), |
если |
п ер еда |
точные функции |
составляю щ их |
ее |
элем ентарны х |
динам ических |
звеньев зап и |
саны в стандартной ф орм е. |
|
|
|
|
|
|
|
|
|
Итак, |
К р = |
0,948; |
|
|
|
|
|
|
|
|
0 < |
клоп < 0,948. |
|
|
§ 4 .3 . |
КРИТЕРИЙ УСТОЙЧИВОСТИ НАЙКВИСТА |
|
1. |
Формулировка критерия |
Найквиста |
|
Замкнутая |
линейная стационарная |
система |
вида рис. |
4.6 |
устойчива, если .приращение аргумента функции |
F (у <в) = |
1 + |
Р и с. 4.6. С труктурная схем а системы
сединичной обратной связью
+W (у'ш) при изменении частоты ш от нуля до бесконечности равно;
Д arg F (у'ш) = |
Д arg 11+ |
W (у’ш)] = — (2Х + \ ) , |
(4.38) |
0<tü><oo |
0 < w < c o |
2 |
|
где
W(j<ü)= W (р)р„^ — амплитудно-фазовая частотная ха
|
|
|
рактеристика |
разомкнутой |
систе |
|
ОТ |
|
мы; |
|
|
|
|
|
|
diPl |
|
|
|
|
|
|
W(p) = —— = |
2 |
|
|
|
|
|
|
—^-----------передаточная |
функция |
разомкну- |
N(p) |
А |
; |
той системы; |
|
|
|
|
|
2-1 |
сіР |
|
|
|
|
|
|
|
!=0 |
dp |
ct — постоянные |
действительные |
коэф |
|
|
|
|
|
фициенты; |
|
|
|
(число |
кор |
|
|
|
X— число полюсов W(p) |
|
|
|
ней N( p ) =0 ) |
с |
положительной |
|
|
|
действительной частью; |
|
|
|
Х0 — число полюсов W(p) |
с нулевой дей |
|
|
|
ствительной частью. |
F (у'ш) = 1 -f- |
Докажем условие устойчивости (4.38). Функция |
+ W(jm) может быть определена в виде: |
|
|
|
|
|
|
|
F ( » — F (p)p^jm , |
|
|
|
|
(4.39) |
где |
|
|
|
N( p)+ M( p) |
|
F{p)= 1 + |
№(р ) = \ + ^ Л р) |
|
|
N{p) |
|
|
|
|
|
N(p) |
|
|
|
15. И зд. № 5312 |
|
|
|
|
|
|
|
225 |
Учитывая, что передаточная функция замкнутой системы вида рис. 4.6 связана с W(p) соотношением (3.4):
ф ( р , |
W ( p ) |
М[р) |
В ( р ) |
|
1 + w{p) |
N (р) + М(р) |
АІР) ’ |
для F(p) можно записать |
|
|
|
|
= N(p) |
(4 |
где А (р)=а0+ а ,р + . ..-\-апрп— характеристический полином зам
кнутой системы;
N (р)=с0+ С]Р+ . . ■+ спрп— характеристический полином ра зомкнутой системы.
Тогда на основании (4.39)
|
|
|
|
|
|
|
|
N(ju) |
|
|
и |
arg F ( » |
= |
arg А (у<о) — arg N (у ш); |
|
|
|
|
|
|
d arg F(jw) = |
Д arg А (/io) — Д arg IV(y'u>), |
(4.41) |
|
0< Ш <!оо |
|
0 < ш < о о |
|
0 < tu •: OO |
|
|
где |
|
|
|
— |
приращения |
аргументов |
соответ- |
Д arg А (_/'cu), Д arg N (j“>) |
< со |
0 < ш - ^ о о |
|
|
|
ственно |
|
А U ш) = |
А (р)р=;ш |
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
N(ju>) = N{p)p-]m |
при из |
|
|
|
|
|
|
менении |
ш |
от нуля до бесконеч |
Приращение |
аргумента |
|
ности. |
Q (у ш) |
произвольного |
годографа |
многочлена |
Q{p) |
при изменении |
<и от нуля до |
'беаконечности |
равно: |
|
|
|
|
|
|
|
|
|
|
|
|
Д argQ (у'ш) = |
|
— [п — 2 &— k0), |
|
(4.42) |
где |
|
0<Ш< СО |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Q |
( у 'ш ) |
= |
Q |
lp)p=ju>\ |
|
|
|
|
|
Q (Р) = |
% + |
|
Я\Р + • • • + |
q„pn; |
|
|
п — число корней Q(p), |
равное его порядку; |
|
|
k — число корней Q'(p) с положительной действительной частью; ko — число корней Q(p) с нулевой действительной частью.
Действительно, если многочлен
Q(p) = q0 + <hP + ■■■+ ЯпРп
имеет постоянные действительные коэффициенты qt, то корни его могут быть нулевыми, действительными или комплексно-со пряженными, т. е. могут иметь вид:
|
Рі = |
0; |
|
|
|
p t = |
+ o-i\ |
at — действительное |
положитель |
{ |
Pi = ± |
a.t + y'ß,-; |
ное число; |
дейетвитель- |
ß, — положительное |
1 |
p i+l = |
+ al — yß,; |
ное число. |
|
Поскольку любой многочлен может быть представлен произведе нием сомножителей вида ( р - РіУ
Q(p) = qn( P - P i ) ( p - P 2 ) ■■■i p — Pi) ■• • ( P - P n l
где pi — корни многочлена, то, учитывая их возможный вид, Q(p) можно представить произведением сомножителей вида:
|
[ р— Рі)=Р, |
если |
а = |
0; |
( р — Рі) = (Р ± аі)< |
если |
р . = |
± а і\ |
<Р- Рі)(Р - |
Рі+1) = |
\Р2+ 2*tP + (®/2 + |
Р/2)Ь если pi “ ± a l +J^t, |
T. e. |
|
|
|
|
|
Pi+ 1 = ± a«-y'ßi, |
k0 |
l |
|
m |
|
|
Q (p) = |
ЯпХ \ р П(P =F ai) П [p2+ 2 я,/? +(«,*+ P,2)]. (4.43) |
|
І = 1 |
£ —1 |
1)=>l |
|
|
где |
k0 — число нулевых корней Q(p); |
|
I — число |
действительных |
корней Q(p); |
|
т — число пар комплексно-сопряженных корней |
Q(P);
n= k0-\-l-\-2m — общее число корней Q (р).
В соответствии с (4.43) годограф многочлена Q(p) будет опреде ляться выражением
*n |
I |
т |
|
Q ( » = ЧпП ( » |
П0 + |
*е) П 1 - “2+У 2»ч“ + («ч2+Рч2)]. (4-44) |
і- 1 |
Е- 1 |
ч -і |
|
т. е. Q (у'ш) с точностью до сомножителя |
дп является произве |
дением годографов вида: |
|
|
|
Qi(у'ш) = уЧ |
|
|
Ql (Уш) = |
Уш + |
(4.43) |
. |
Qi)(У10) = [« + Р ч 2 ) - <°3] + /2Ѵ°> |
которые показаны на рис. 4.7, а, б, в.