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

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

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

и исследовать его корни Хк (6=1, 2,..., п). Ляпунов доказал две теоремы, которые позволяют исследовать все некритические случаи.

Теорема 1. Если вещественные части всех корней Хк харак­ теристического уравнения (5.9) 1-го приближения отрицательны, то невозмущенное движение асимптотически устойчиво незави­ симо от членов разложения выше 1-го порядка малости.

Теорема 2. Если среди корней Хк характеристического урав­ нения (5.9) 1-го приближения найдется по меньшей мере один

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

Критические случаи имеют место, когда среди всех корней

уравнения (5.9) имеются некоторые корни, вещественная часть

которых равна нулю, а остальные корни имеют отрицательную

вещественную часть. В критических случаях вопрос об устойчи­

вости невозмущенного движения (5.6) не может быть разрешен

на основании исследования уравнений 1-го приближения: устой­

чивость (неустойчивость) невозмущенного

движения

опреде­

ляется видом нелинейных функций Fk. Поэтому в критических

случаях требуется рассматривать уравнения

(5.7)

в исходном

виде.

соответствующее

системе

Характеристическое уравнение,

уравнений (5.8), имеет вид

 

 

 

 

(5.10)

ОпХ-{-йп—\Хп~Ч"...

 

 

 

 

Пусть для определенности все корни уравнения

(5.10) различ­

ны, тогда его решение записывают в виде

 

 

 

 

x=AxeXit + Л2е*”'+ ... + Л„ел"',

 

 

 

 

 

где А.1, fa,...» Ап—корни характеристического

уравнения; Ах,

Аг, ..., Ап —постоянные интегрирования, зависящие от началь­

ных условий.

Если Хк>0,

то

член

Пусть Кк —вещественный корень.

Акъхк*с течением времени непрерывно возрастает и стремится к бесконечности. В этом случае х также стремится к бесконечно­

сти и система неустойчива. Если À*<0, то член Ак№ с тече­ нием времени стремится к нулю, т. е. затухает.

Пусть один из корней Хг_—комплексный, тогда всегда суще­ ствует сопряженный с ним %г:

Яг—ос-|-/Р; Хг—а—/р. В этом случае

Аге%г*+ Аге>г‘= Аге?* з1п(р*+ф).

Если сс>0, то имеют место колебания с частотой р и нараста­ ющей амплитудой, т. е. движение неустойчиво. При а=0 получим незатухающие колебания —система находится на гра­ нице устойчивости. Если а<0, то амплитуда колебаний с тече­

нием времени уменьшается и колебания затухают. Отсюда мож­ но сделать следующие выводы:

если все вещественные части корней характеристического уравнения отрицательны, то динамическая система устойчива (рис. 5.1):

 

Л;

*

Плоскостьs

 

 

*

*

!

 

■Л-6 х

Aj

1

Вещественная

ÀS

*

0

есь

 

 

 

X

Аг

 

 

 

X

 

 

Рис. 5.1. Расположение корней

характеристи­

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

САР на ком­

плексной плоскости

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

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

 

Если в каких-либо корнях характеристического уравнения

вещественная часть равна нулю, а у

остальных —отрицатель­

ная, то об устойчивости невозмущенного движения по первому

приближению ничего сказать нельзя и требуется специальное

исследование полного уравнения. Наконец,

если среди корней

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

один или несколько

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

5.2. Критерии устойчивости линеаризованных САР

Вычисление корней характеристического уравнения не пред­ ставляет труда для уравнений 1-й и 2-й степеней. Что касается общих выражений для корней уравнений 3-й и 4-й степеней, то они громоздки и практически мало удобны. Следует отметить отсутствие общих выражений для корней в уравнениях более высоких степеней. Поэтому важное значение приобретают пра­ вила, которые позволяют определить устойчивость системы, ми­ нуя вычисления корней. Эти правила называют критериями устойчивости. Они позволяют в ряде случаев не только устано­ вить, устойчива система или нет, но и выяснить влияние тех

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

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

чива ли система, с помощью только алгебраических процедур над коэффициентами характеристического уравнения, называют алгебраическими. К ним относят критерии устойчивости: Рауса, Гурвица, Шур-Кона и др. [2, 19, 20]. Алгебраические критерии для систем, описываемых уравнениями выше 4-й степени, дают1 возможность определить лишь устойчивость системы при задан­ ных численных значениях коэффициентов уравнения. Но за­ труднительно с их помощью ответить на вопрос: как изменить параметры системы, чтобы сделать ее устойчивой?

Частотный критерий устойчивости, впервые сформулирован­ ный Найквистом, был применен для исследования устойчивости

САР А. В. Михайловым в 1936 г. Кроме трго, последний сфор­ мулировал другой частотный критерий, получивший название критерия устойчивости Михайлова. Достоинством частотных

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

Об устойчивости замкнутой системы судят по частотной ха­ рактеристике разомкнутой, и в этом случае применяют крите­ рий устойчивости Найквиста—Михайлова. Кроме того, частот­ ные критерии устойчивости дают представление и о качестве процесса регулирования.

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

Критерий устойчивости Рауса. Этот критерий формулируют следующим образом: если система автоматического регулиро­ вания описывается линеаризованным характеристическим урав­ нением вида (5.10), то для того, чтобы система была устойчива (т. е. все корни уравнения имели отрицательные вещественные части), необходимо и достаточно, чтобы все элементы столбца 1 табл. 5.1 для данного уравнения были одного знака.

Если а„>0, то все элементы столбца 1табл. 5.1 должны быть положительны.

Таблицу (алгоритм) Рауса (см. табл. 5.1) составляют сле­ дующим образом: в строку 1 вписывают коэффициенты уравне­ ния (5.10) с индексами (а„, а„-2, а„-4, ...); в строку 2 —коэф­ фициенты уравнения с индексами (ап-ь а„-з, an-s,... ); в строку

-

г-**.

II С

г —ÛL с1А

 

 

АлгоритмРауса

 

Номер

I

Столбец

3

Строки

2

1

ап

Ûn-2

ûn-4

2

fln-l

an-3

an-s

3 Cl3=dn-2—r(Дn-3 Cas=an-4—foûn—5Сзз=ап-б—гQün—j 4 Сц=ап-1—riC23 c24=On-5—+\Сзз C34=an-7—Г1С4Э 5 Clb~C23—*2024 С25=Сзз—Г2С34 Сзе=С4э—Г2С44

CI,i+l

1+3

cl,l+3=

c2,i+3=

 

Сэ.<+3=

Г1 cl,i+2

 

—c2.i+l—№2,1+2

=c3,i+l—fiCi,i+2

=64,14-1—Г<С4,|+2

■3 —коэффициенты с]3, с2з, которые подлежат

определению.

-Впоследующие строки вписывают коэффициенты см (k —номер столбца; i—номер строки, в которой стоит коэффициент). Каж- -дый из коэффициентов сi3l с2з, ..., cki равен определителю; пер­ вый столбец определителя составлен из двух элементов, запи­ санных в следующем за искомым коэффициентом столбце таб­ лицы на двух расположенных выше строках. Первый элемент -второго столбца определителя образован из частного отделения -двух элементов, расположенных в столбце 1табл. 5.1 на двух •вышележащих строках. Второй элемент второго столбца опре­

делителя равен единице. Так,

С\Ъ~

Û'|"2ап~i I=йп_2--‘ -an-3Î

 

a„-i 1

I

Ск1=

Ck+\,i-\

ri-31

 

 

1 Г

где г/_3—с1,1—2/Cl,1-1.

Значения г вписывают в боковик табл. 5.1, озаглавленный «Зна­ чение г». На них умножают соответствующие коэффициенты.

Из критерия Рауса следуют выводы:

1)все коэффициенты характеристического уравнения устой­

чивой системы должны быть одного знака. Обращение в нуль

одного из коэффициентов а,- (за исключением коэффициента

старшего члена) свидетельствует о неустойчивости системы или

о том, что она находится на границе устойчивости. Если коэф­

фициенты характеристического уравнения положительны, то все

вещественные

корни, если они

существуют, отрицательны (так

называемые «левые» корни). Комплексные корни могут быть и

«правыми»;

отрицательных

коэффициентов Сц столбца 1

2) число

табл. 5.1 равно числу корней

с положительной вещественной

частью;

 

 

3) обращение в нуль ао приводит к появлению нулевого кор­

ня. Обращение в нуль последних v коэффициентов я0=0; ах=

=0,... ач-1=0 приводит к появлению нулевых корней. При этом

обращаются

в нуль последние

коэффициенты си табл. 5.1

(^1,Л C\t fl—1

• • . “ £Il П—v+1 0),

 

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

Критерий устойчивости Гурвица. Этот критерий легко полу­ чить из критерия устойчивости Рауса. Для данной цели выра­ зим коэффициенты cik в виде определителей:

С\\—

С\о—

—Ajî

 

Дд-зДд-г—ДдДд-з

I Дл-j Дя I

д,

1Дд—эДя-2 I

с13—

Дя-|

 

Дл-1

А,

а,

Дд

I

I Дл-j Дд

I

Дд-г I

п1I Дд-t Дд-4 I

С\4== '

I Дд-l Дд I

 

 

I Дд-а Дл-a I

 

а в общем случае

Дд-t Дд

0

(5Л1)

 

Дд-з Дд-г Дд-i

А,

Дд-S Дд-4 Дд-З

I Дд-i Дд

I

Аз

I Дд-з Дд-г I

 

Afe-i Ал—'

где Лл(^ = 1.2, ...) —определители Гурвица, получаемые с по­ мощью следующей записи:

Дд-ll Дд I 0

0

0

0 ...

Дд-з Дд-г| Дд-i

Дд

0

0 ...

Дд-5 Дд-4 Дд-3

Дд-2 Дд-1

ап ...

Дд—7Дд-б Дл-5 Дд-4

Дд-3

Дд-2 • • •

Дд—9Дд-8 Дд-7

Дд-б

Дд-5 Дд-4 • • •

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

UnnuaÜÜ.^ "I1! Строки 0ПРеделителя, начиная с диагонали, за-

Тк®эФФициентами: вправо —по убывающим индексам,

авлево —по возрастающим.

Согласно критерию Рауса, необходимым, и достаточным ус­ ловием устойчивости являются соотношения

Сц-ап>0; cn=>an-i>0) Ci3>0;... ; cit я+1>0.

Этим неравенствам, как следует из (5.12), эквивалентны нера­ венства вида

ап> 0; Д,>0; Д2>0;..... Дя>о.

Таким образом, критерий Гурвица формулируют следующим образом: для того чтобы система была устойчива, необходимо и достаточно, чтобы выполнялось неравенство'а„>0, а определи­ тели Гурвица Д], Д2,..., Д„ были положительны.

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

Необходимые, но недостаточные условия устойчивости за­ ключаются в том, что в случае уравнения /i-го порядка все коэффициенты ап, а„_ь ..., а0должны быть положительны и ни один из них не должен равняться нулю.

Рассмотрим характеристические уравнения и условия устойчивости для динамических систем, порядок которых не превышает пяти:

1) aiA,+a0=0. Условия устойчивости: а0>0, а,>0; 2) агА,2+аЛ+ао=0.

Условия устойчивости: а0>0, а,>0, а2>0- 3) a3^3+a^2*+aiX-fа0=0. Условия

устойчивости: ао>0, fli>0, а2>0, аэ>0, Д2=а,а2—аоа3>0; 4) аДЧ-аЛЧ-

+аА+а0=0.

Условия

устойчивости: а0>0,

ai>0, а2>0 а3>0 а4>0,

Д2=а2а3—а,а4>0, Д3=а,а2а3—а,2а4—а0а32>0; 5)

asA.8+aA4+a3b3+a2>>+aib+

-fa0=0. Условия устойчивости: а0>0, ai>0, а2>0, а3>0, а4>0, а5>0,

а4а%

 

 

Д2= аь

\=ata4—a.as>0,

 

а4 а2&Q

 

 

Д«= 0

 

Л\ =atага4+а0а„а5—ai2a1—ai2as>0,

а4а2

 

 

ак аг а00

 

 

_

д3ûj 0

Д]Д2Д3Д4+2ЛдД]Д4Д3-J- Дфй2Л]Д5—

0

а4аг д0

0

д5Д3o.î

 

 

—аха£аг—а0а3:а4—a02asz—а,гд42>0.

Для п=2 условиемустойчивости является лишь положительность коэф­ фициентов характеристического уравнения. Для п=3, п=4, п=5 положи­ тельность коэффициентов характеристического уравнения недостаточна. Кро­ ме того, коэффициентыдолжныудовлетворять дополнительным неравенствам.

5.4.

Частотные

критерии устойчивости.

 

Критерий устойчивости Михайлова.

 

 

Амплитудно-фазовый критерий устойчивости

 

Частотные критерии в большинстве случаев используют в

качестве графоаналитических критериев —они отличаются

на­

глядностью при выполнении инженерных расчетов.

аргумента —

В основе частотных методов лежит

принцип

следствие из теоремы теории функций комплексного перемен­

ного, а именно теоремы Коши, относительно числа нулей

и

полюсов функции, аналитической в заданной области.

 

Принцип аргумента. Рассмотрим

алгебраическое уравне­

ние п-й степени с действительными коэффициентами:

 

Æ(A) =a„An+an_iAn-1+ ... +а|Л+а0=0.

 

 

 

Если через Ai, А2,

..., А„ обозначить корни этого уравнения, то

многочлен D(A) можно представить в виде произведения прос­

тых сомножителей:

 

 

 

 

Д(А)=а„(А—AО(А-А2) ... (А-А„).

корню

соответствует

На комплексной

плоскости

каждому

вполне определенная точка

(рис. 5.2,а). Геометрически каж­

дый корень Xi изображается в виде вектора, проведенного из

начала координат к точке А,-

(рис. 5.2,6). Длина этого вектора

Рис. Б.2. Корни характеристического уравнения системы: а—расположение корнеЛ; б—модуль ифаза нектора X/

равна модулю комплексного числа, т. е. |Я(|, а угол, образо­ ванный вектором с положительным направлением действи­ тельной оси,—аргументу или фазе комплексного числа А,-, т. е. argA(. Векторы {X—А<), входящие множителями в D(X), прове­ дены из точек Xi к точке А.Каждый из этих векторов является разностью двух векторов, соответствующих Xи Х{ (рис. 5.3). Если принять А—/©в D(X), то

D(j(ù)=an(j(ù—Ai) (у©—X2) • • ■(/со—A„). где ©—круговая частота (см. раздел 3.3.).

Рис. 5.3. Элементарный вектор (Я-Я<)

Концы элементарных векторов будут находиться на мнимой оси в точке Я=/о (рис. 5.4). Модуль этого вектора равен про­ изведению модулей элементарных векторов и а„:

\D(j(ù)\=an\j(ù—ki\\j(ù—\2\ ... I/o—Яп1>

a аргумент или фаза его равна сумме аргументов элементар­ ных векторов:

arg£(/cù)=arg(/0—M)+arg(/o—Я2) + • • •

...+arg(/o—Я„).

Направление вращения вектора против часовой стрелки принимают за положительное. Тогда при изменении о от —оо до +оо каждый элементарный вектор (/©—Я<) повернется на

 

 

угол +я, если его начало

(корень Я<)

 

 

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

 

 

плоскости, и на угол —я,

если

его

 

 

начало (корень ЯЛ) находится в пра­

 

 

вой части

комплексной

плоскости

 

 

(рис. 5.5). Если уравнение О(Я)=0

 

 

имеет т корней в правой части пло­

 

 

скости Яи, следовательно, п—т

кор­

 

 

ней—в левой

части

комплексной

 

 

плоскости, то при возрастании о от

 

 

—оо до +.°° изменение

аргумента

Рнс. 5.5. Изменение аргу-

вектора Z)(/со),

или угол

поворота

/)(/©) (равный

сумме

изменений

ар-

“/ш-Л.)еКо°раВвозрастании

гументов элементарных векторов), бу-

частотыо от —» до +оо

дет

 

 

 

 

 

Д

arg Û(J(ù)=(n—m)n—mn={ti~2m)n.

 

(5.13)

 

—00<Û)<00

 

 

 

 

 

 

Отсюда следует, что разность (п—т )

корней уравнения D(

 

=0,

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

корней, рас­

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

ражает собой изменение аргумента

/)(/©) при возрастании ©

от —оо до +оо. Это утверждение

в теории автоматического

регулирования называют принципом аргумента.

Критерий устойчивости Михайлова. Этот критерий основан

на принципе аргумента (5.13) и является геометрической ин­

терпретацией соотношения

 

Д arg D{j(ù)= (n—m)''7i—mn=(n—2m)n,

—0O<CD<CO

 

где пг —число корней в правой части комплексной плоскости;

(п—т) —число корней в левой части комплексной плоскости.

Пусть характеристическое уравнение системы с обратной

связью (замкнутой САР) имеет вид

 

 

D(X)=аДп+а„_|Я”~Ч-... +a\k+ao==Q.

 

части

Если все

корни этого уравнения находятся в левой

комплексной плоскости Я(система устойчива), а в правой части

плоскости корней нет, то т=0 и изменение аргумента

 

Д arg

D(j(û)=mi.

 

 

—00<(0<00

 

 

Отсюда следует вывод: САР является устойчивой, если

при

возрастании ©от —оо до +оо изменение

аргумента вектора

D(j(ù) будет равно пя, где п —степень

характеристического

уравнения /}(Л)=0.

При изменении частоты ©от —оо до +оо вектор £>(/©) на

комплексной плоскости опишет своим концом кривую, которая называется характеристической кривой, или годографом век­ тора /)(/©).

Уравнение характеристической кривой определяют под­ становкой Я=/© в многочлен D(X) и последующим разделе­ нием действительной и мнимой частей:

£)(/©)=а„(/©)"+а„-1(/©)я-,+... +ai(/©)+a0;

D(j(ù)=w(©)+/u(©), где

U(©) =ÜQ—a2(Û2+Û4©4—..., о(©)=al(ù—a3©3+as©—...

Действительная часть «(©) является четной функцией, а мни­ мая часть и(©) —нечетной функцией частоты ©, т. е. и(—©) = =ц(©), и(—©)=—о(©). Поэтому для отрицательных значе­

ний ©

D (—/©)=и(©)—jv (©).

Отсюда следует, что характеристическая кривая симметрич­ на относительно действительной оси для +© и —©. При по­ строении характеристической кривой можно ограничиться лишь

положительными значениями ©от 0 до °о. При этом угол по­

ворота вектора D(j<û), т. е. изменение аргумента £>(/©), умень­

шается вдвое и критерий устойчивости формулируется следую­

щим образом.

 

 

САР будет устойчивой, если при возрастании частоты ©от

Одо

вектор

£>(/©) повернется на угол

пл/2 (где п —

степень

уравнения

£(Я)=0). Это означает,

что вектор харак­

теристической кривой (при изменении частоты от 0 до +оо(

начиная с положительной действительной оси) последователь­

но «обходит» п квадрантов

в положительном направлении,

т. е. против часовой стрелки.

характеристики, соответствующие

На рис. 5.6 приведены

устойчивой системе. При п—1изменение аргумента равно я/2, при п=2 изменение аргумента равно л и характеристическая кривая проходит через два квадранта, и т. д.

(n=1...5)

На рис. 5.7 приведена характеристическая кривая для

=4, которая соответствует неустойчивой системе. Система бу­

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

тическая кривая при некотором значении пересекает

начало

координат, обходя при этом

(п—1) квадрантов. Частота ©яв­

ляется одновременно корнем

уравнений и(©)=0 и о(©)=0.

В ряде случаев может быть использован критерий устойчи­

вости, называемый критерием перемежаемости корней

(рис. 5.8

и 5.9). Действительно, характеристическая кривая

при из­

менении ©от 0 до со будет обходить в положительном на­ правлении п квадрантов и система устойчива, если ц(0)>0, t>(0)=0 и уравнения ц(©)=0 и ü(©)=0 имеют все действи­ тельные и перемежающиеся корни, т. е. если между каждыми двумя соседними корнями о(©)=0 лежит корень уравнения u(œ)=0 или между двумя соседними корнями w(©) находится корень уравнения и(©)=0.