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

книги / Теория автоматического управления

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
16.96 Mб
Скачать

Корни характеристического уравнения могут быть действи­ тельными (рк = ак), мнимыми (рк = /рй) и комплексными

р* = а * ± /р * .

(5.6)

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

Переходная составляющая (5.4) при t -> оо стремится к нулю

лишь в том случае, если каждое слагаемое вида Скерь* 0. Ха­ рактер этой функции времени зависит от вида корня рк. Рассмот­ рим все возможные случаи расположения корней рк на комплекс­ ной плоскости (рис. 5.1) и соответствующие им функции хк (/), ко­ торые показаны внутри кругов (как на экране осциллографа).

1.

Каждому действительному корню рк = ак

в решении (5.4)

соответствует слагаемое вида

 

x(t) = Ckea^

(5.7)

Если ак < 0 (корень р х), то функция (5.7) при t ->• оо стремится к нулю. Если ак > 0 (корень р3), то функция (5.7) неограниченно

Рис. 5.1. Влияние корней характеристического уравнения системы на со­ ставляющие ее свободного движения

возрастает. Если ак = 0 (корень р ^ , то эта функция остается по­ стоянной.

2. Каждой паре сопряженных комплексных корней рк =

— «* + /Р* н Л +1 = ак —/Р* в решении (5.4) соответствуют два слагаемых, которые могут быть объединены (см. 3.4) в одно слагае­

мое

xft(0 = 3 C te ^ sin (M + ^ * )-

(5.8)

Функция (5.8) представляет собой синусоиду с частотой |Jt и ам­ плитудой, изменяющейся во времени по экспоненте. Если действи­ тельная часть двух комплексных корней ак (см. рис. 5.1, корни Р* и Pi)t то колебательная составляющая (5.8) будет затухать. Если ак > 0 (корни ра и р»), то амплитуда колебаний будет неограни­ ченно возрастать. Наконец, если ак = 0 (корни pt и р7), т. е. если оба сопряженных корня — мнимые (рк = + /Р*> Pk+i = — /Р*), то хк (/) представляет собой незатухающую синусоиду с частотой р*.

Если среди корней характеристического уравнения (5.5) имеются / равных между собой корней ph то в решении (5.4) вместо / слагае­

мых вида СкеР^ появится одна составляющая

(C0-f Сх/ + С 2^ + .. . + CtJ *) еР/<.

(5.9)

Учитывая, что функция вида е~ы при любом Ь убывает быстрее, чем возрастают слагаемые вида tr, можно доказать, что и в случае кратности корней решение (5.4) будет стремиться к нулю лишь при отрицательности действительной части кратных корней pt.

На основании проведенного анализа можно сформулировать о б щ е е у с л о в и е у с т о й ч и в о с т и :

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

При этом действительные корни рассматриваются как частный случай комплексных корней, у которых мнимая часть равна нулю. Если хотя бы один корень имеет положительную действительную часть, то система будет неустойчивой.

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

Используя геометрическое представление корней (5.6) на ком­ плексной плоскости (см. рис. 5.1) в виде векторов или точек, можно

дать в т о р у ю ф о р м у л и р о в к у о б щ е г о

у с л о в и я

у с т о й ч и в о с т и

(эквивалентную основной):

 

для устойчивости линейной системы необходимо и достаточно,

чтобы

все корни

характеристического уравнения

находились

в левой

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

 

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

Мнимая ось /р является границей устойчивости в плоскости корней. Если характеристическое уравнение имеет одну пару чисто мнимых корней (рк = + /Р*, Pk+i = —/Рл), а все остальные корни находятся в левой полуплоскости, то в системе устанавливаются незатухающие гармонические колебания с круговой частотой со = = |p fe|. В этом случае говорят, что система находится на колеба­ тельной границе устойчивости.

Точка р = 0 на мнимой оси соответствует так называемому нулевому корню. Если уравнение имеет один нулевой корень, то система находится на апериодической границе устойчивости. Если таких корня два, то система неустойчива.

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

илинеаризации исходных нелинейных уравнений. Возникает во­ прос: в какой мере оценка устойчивости по линеаризованному уравнению будет справедлива для реальной системы, не окажут ли существенное влияние на результат анализа отброшенные при ли­ неаризации члены разложения? Ответ на него был дан русским ма­ тематиком А. М. Ляпуновым в 1892 г. в работе «Общая задача об устойчивости движения». Он сформулировал и доказал следующую теорему: если характеристическое уравнение линеаризованной системы имеет хотя бы один нулевой корень или одну пару мнимых корней, то судить об устойчивости реальной системы по линеари­ зованному уравнению нельзя. Отброшенные при линеаризации малые члены могут сделать систему устойчивой или неустойчивой,

ипоэтому устойчивость реальной системы необходимо оценивать по исходному нелинейному уравнению.

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

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

Простейшим критерием устойчивости является условие положи­ тельности коэффициентов характеристического уравнения. Поло­ жительность коэффициентов уравнения (5.4) является необходимым (но не достаточным!) условием устойчивости системы. Это означает, что если все коэффициенты положительны, то система может быть устойчивой или неустойчивой. Но если хотя бы один коэффициент уравнения отрицателен или равен нулю, то система неустойчива.

Критерии устойчивости могут быть алгебраическими и частот­

ными. Алгебраические критерии устанавливают необходимые и до­ статочные условия отрицательности корней в форме ограничений,

накладываемых на определенные комбинации коэффициентов ха­ рактеристического уравнения. Частотные критерии определяют связь между устойчивостью системы и формой частотных характе­ ристик системы.

При анализе устойчивости систем управления обычно решают одну или несколько задач: 1) оценивают, устойчива или нет система при заданных параметрах; 2) определяют допустимый по условию устойчивости диапазон изменения некоторых незаданных парамет­ ров системы; 3) выясняют, может ли система при заданной струк­ туре быть в принципе устойчивой.

Первая задача может быть решена с помощью критериев, из­ ложенных в 5.2—5.4, вторая — построением областей устойчиво­ сти (см. 5.5), третья — с использованием условий структурной устойчивости (см. 5.6).

5.2. Алгебраические критерии устойчивости

Наиболее распространены в инженерной практике критерии Гурвица и Рауса (ниже приведены формулировки и методика при­ менения этих критериев без доказательства их справедливости).

Критерий Гурвица был сформулирован и доказан в 1895 г. не­ мецким математиком А. Гурвицем, который разработал свой кри­ терий, решая чисто математическую задачу — задачу исследова­ ния устойчивости решений линейного дифференциального уравне­ ния. Гурвиц обратился к этой задаче по просьбе словацкого уче­ ного А. Стодолы, занимавшегося вопросами регулирования турбин.

Применительно к задачам теории управления критерий Гур­ вица можно сформулировать так:

автоматическая система, описываемая характеристическим урав­ нением

а0рп+ а1рл—1+ а2рп~-2+ .

устойчива,

если при&оV

^2» •

» к

вида

.+

II о

(5.10)

О положительны все определители Ai,

ах а3 аь

Ао а4

д<= 0 Oi а3

#21-1

#2f-2

#2f-3 , « = 1, 2,

п.

(5.11)

0

Если хотя бы один из определителей (5.11), называемых опреде­ лителями Гурвица, отрицателен, то система неустойчива.

Матрицы (5.11), по которым вычисляют определители Гурвица, составляют следующим образом: на главной диагонали записывают все коэффициенты характеристического уравнения от а 1до а( (в по­ рядке возрастания индекса), затем в каждом столбце выше диаго­ нальных коэффициентов записывают коэффициенты с последова­ тельно возрастающими индексами, а ниже — с последовательно убывающими индексами; на место коэффициентов с индексами боль­ шими п или меньшими нуля проставляют нули. При этом каждая г-я матрица получается квадратной размером i х i.

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

А„ = а„Ал_1.

(5.12)

Если главный определитель Д„ = 0, а все остальные определи­ тели положительны, то система находится на границе устойчивости. С учетом выражения (5.12) это условие распадается на два:

ап= 0 и Д„_1 = 0.

(5.13)

Условию ап — 0 соответствует один нулевой корень, т. е. апе­

риодическая граница устойчивости, а условию An_x =

0 — пара

мнимых корней, т. е. колебательная граница устойчивости.

Рассмотрим

ч а с т н ы е с л у ч а и к р и т е р и я

Г у р ­

в и ц а для п =

1; 2; 3; 4. Раскрывая определители, фигурирующие

в общей формулировке критерия, можно получить следующие ус­ ловия.

1. Для уравнения первого порядка

 

а0р + ах = 0

(5.14)

условие устойчивости

 

а0> 0 и Д1 = а1> 0 ,

(5.15)

т. е. положительность коэффициентов уравнения является в данном случае и необходимым и достаточным условием. Действительно,

при а0 >

0 и аг > 0 единственный корень уравнения

будет отри­

цательным: /?1 =

(a ja 0) < 0.

 

2. Для уравнения второго порядка

 

OQP2

й хр а% = 0

(5.16)

условие устойчивости

 

ао> 0 ,

Ai = a i> 0 ,

(5.17)

A2 = a2A i> 0

или а2> 0.

 

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

3. Для

уравнения третьего порядка

 

 

а0р3+

aiР2+ а2р + а3= О

 

 

(5.18)

условие устойчивости

 

 

 

а0> 0,

Ai = a i> 0 ,

 

 

 

Д.=-

ai

а3 = а1а2— а0а3> 0 ,

Д3 = а3Д2> 0 .

(5.19)

 

во

+

 

 

 

 

Последнее

неравенство при

а3 > 0

эквивалентно неравенству

Д2> 0 .

Следовательно, для системы третьего порядка, кроме по­

ложительности

всех коэффициентов,

требуется,

чтобы Д2> 0 .

Учитывая выражение для Д2, можно сформулировать мнемониче­ ское правило оценки устойчивости систем третьего порядка:

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

4. Для уравнения четвертого порядка

а0р*-f- + р 3 + а2р2+ а3р -)- а4 = 0

(5.20)

кроме положительности всех коэффициентов требуется выполнение условия

Д3 = а1а2а3— a0ai— а ^ Х ) .

(5.21)

Нетрудно доказать, что при положительности всех коэффици­

ентов

условие (5.21) обеспечивает выполнение условия Д2> 0 .

Таким

образом,

для устойчивости систем не выше четвертого порядка необхо­ димо и достаточно, чтобы все коэффициенты характеристиче­ ского уравнения и определитель Дп_г были положительными.

Критерий Гурвица целесообразно применять для анализа устой­

чивости

систем

не выше

пятого порядка. При п >

5

вычисление

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

становится

громоздким.

 

 

Пример 1. Определим с помощью критерия Гурвица устойчивость си­

стемы управления частотой вращения двигателя (см. 4.2)

при

следующих

значениях параметров:

 

 

 

 

 

 

 

Т м =

1 с;

Тя = 0,1

с; Т Г = 0,5 с;

k = АуАт.пМ д*т г =

14.

(5 22)

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

уравнение

системы

 

 

(Ттр + 1) (Г„Гяр2 +

Тир + 1) +

А =

0.

 

(5-23)

Приводя это уравнение к форме (5.18), получим значения коэффициентов

а» =

Т иТ яТ г = 1 -0,1 -0,5 = 0,05

с3;

 

 

 

«1 =

Т’м'Гя +

Т иТ г =

1-0,1 +

1-0,5 =

0,6 с*;

 

(5 24)

аг = Т и Ц -Тг = 1 +

0,5 =

1,5

с;

 

 

 

 

о3 =

1 +

k = 15.

 

 

 

 

 

 

 

Все коэффициенты характеристического уравнения положительны, т. е. необходимое условие устойчивости выполняется. Проверим выполнение до­ статочного условия, для чего вычислим определитель

Д 2 =

а\а2аьа3= 0,6 • 1,5 — 0,05 • 15 = +0,15,

(5.25)

А2 > 0 ,

следовательно, система устойчива.

 

Решим теперь обратную задачу: определим, какое максимальное значе­ ние общего передаточного коэффициента k допустимо но условию устойчи­ вости.

Максимально допустимое значение коэффициента k найдем из условия нахождения системы на колебательной границе устойчивости:

Ал_1 = Д2 =

ЯиЯ3 ^ 0*

(5.26)

Отсюда

 

 

а3 < а ^ /а ,, =

0,6-1,5/0,05 = 18,

(5.27)

а максимально допустимое значение общего

передаточного коэффициента

k = a3— 1 <

18 — 1 = 17

(5.28)

Условию нахождения системы на апериодической границе устойчивости (ап = 0 ) соответствует второе предельное значение передаточного коэффи­ циента

/г > — 1.

(5.29)

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

Отметим, что точность системы в режиме положительной обратной связи совершенно неудовлетворительна.

Критерий Рауса, предложенный в 1877 г. английским матема­ тиком Э. Дж. Раусом, целесообразно использовать при анализе устойчивости систем выше четвертого порядка. Для этого из ко­ эффициентов характеристического уравнения (5.10) составляют таблицу (табл. 5.1), в первой строке (i = 1) которой записаны ко­ эффициенты уравнения с четными индексами, во второй (i = 2) — с нечетными индексами, в последующих строках (i > 3) помещены коэффициенты Рауса,полученные как комбинации коэффициентов двух вышестоящих строк по формуле

Гik = ri-2, k+1i-2, ir1-1, k+i!ri-i» 1

(5.30)

где i — номер строки, k — номер столбца. Сам критерий формулируется так:

автоматическая система устойчива, если положительны все коэффициенты первого столбца таблицы Рауса (включая а0 и аг).

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

Алгоритм вычисления коэффициентов (5.30) легко запрограмми­ ровать, поэтому критерий Рауса используют для анализа систем высокого порядка (п > 5) с помощью ЭВМ.

Преимуществом критериев Гурвица и Рауса является то, что

сих помощью можно оценивать устойчивость как замкнутых, так

иразомкнутых систем. Вывод об устойчивости при применении этих критериев делается применительно к той системе (замкнутой или разомкнутой), уравнение которой анализируется.

Пример 2. Определим с помощью критерия Рауса устойчивость системы, описываемой характеристическим уравнением

2рв +

3ръ + 4р* + бр8 +

6р2 + + 8 =

0.

 

 

 

(5.31)

Коэффициенты Рауса:

 

 

 

 

 

 

 

 

 

гм =

а> =

2;

г-12 — Аг =

4;

г1з =

аА=

6;

гц =

а $ =

8',

 

r2i =

fli =

3;

г22 =

аз =

5;

roz =

аъ — 7;

г2\~

 

 

r3i — ri2

(ги/22^ 21) =

а 2 — (a0a3/ai) =

4 — (2• 5/3) «

0,67;

 

Г32 =

г1 з

(riir23/^2 1) =

а 4 (а а ъ!а \) =

6 — (2*7/2) =

— 1;

 

гзз =

ri4 — (rn r2i/r2i) =

 

аб — fl.r0 = 8;

г34 = 0;

г35 =

0;

 

/-41 =

/22 -

( Ы м Ы

=

 

5 — [3 ( — 1)/0,67] «

9,48;

 

 

/"42 == /*гз

2irзз/г3 1 ) =

7 — (3*8/0,67)=

— 28,8;

/"43 =

0;

^44 = 0 ;

' 5* =

г32 — (г31г42/г41) =

 

— 1 — [0,67 ( — 28,8)/9,48] =

1,25;

 

ГЬ2 =

гзз

(r3i/’43//'4i) =* 8;

гбз =

0;

 

 

 

 

 

rei = г42

(rAirbJrb\)=

 

—28,8 — (9,48*8/1,25)= —89,5;

г32 = 0‘»

Г71 = г 52 — (Г51г s J r e i) = 8 .

Один из коэффициентов первого столбца отрицателен (гб1 = — 89,5), следовательно, система (5.31) неустойчива. Так как знак коэффициентов

этого столбца меняется дважды (до и после гв1), то уравнение (5.31) имеет два правых корня.

5.3. Критерий Михайлова

Критерий Михайлова относится к группе частотных критериев устойчивости. Он был сформулирован и обоснован в 1936 г. совет­ ским ученым А. В. Михайловым в работе «Гармонический метод в теории регулирования», которая получила высокую оценку и по­ служила началом широкого применения частотных методов в тео­ рии автоматического управления.

Критерий Михайлова так же, как критерии Гурвица и Рауса, основан на анализе характеристического уравнения системы, поэ­

тому с его помощью можно судить об устойчивости

замкнутых

и разомкнутых систем.

 

Пусть левая часть характеристического уравнения, называе­

мая характеристическим полиномом, имеет вид

 

F (р) = а0рп+ а ^ - 1+ . . + ап.

(5.32)

Подставим в этот полином вместо переменного р чисто мнимый корень, который в дальнейшем будем обозначать /со. Тогда получим функцию комплексного переменного

F (/со) = а0(/со)" + ai (/©)п“1 + •• • + ап>

(5.33)

которую можно так же, как амплитудно-фазовую характеристику, представить в виде суммы действительной и мнимой частей:

F (/со) = Р (со) + /(? (со).

(5.34)

Действительная часть Р (со) содержит только четные степени переменного со:

Р (со) = апап_2со2 + ал_4со4 — .

(5.35)

а мнимая часть Q (со) — только нечетные:

 

Q (со) = а„_!(о—ап_3со3 +

(5.36)

Каждому фиксированному значению переменного со соответст­ вует комплексное число, которое можно изобразить в виде вектора на комплексной плоскости. Если теперь изменять параметр со от О до оо, то конец вектора F (jсо) опишет некоторую линию (рис. 5.2, а),

которая называется характеристической кривой или годографом Михайлова. По виду этой кривой можно судить об устойчивости

системы.

Ф о р м у л и р о в к а к р и т е р и я М и х а й л о в а :

автоматическая система управления, описываемая уравнением п-го порядка, устойчива, если при изменении со от 0 до оо ха-

160

рактеристический вектор системы F (/©) повернется против часовой стрелки на угол пя/2, не обращаясь при этом в нуль.

Эго означает, что характеристическая кривая устойчивой си­ стемы должна при изменении © до 0 до оо пройти последовательно через п квадрантов. Из выражений (5.35) и (5.36) следует, что кри­ вая F (j<o) всегда начинается в точке на действительной оси, уда­ ленной от начала координат на величину ап.

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

Если кривая F (jсо) проходит через начало координат, то си­ стема находится на границе устойчивости. Действительно, если характеристическое уравнение имеет один нулевой корень р* = О (апериодическая граница устойчивости) или одну пару чисто мни-

170