книги / Теория автоматического управления техническими системами
..pdfи исследовать его корни Хк (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.