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

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

Определитель Гурвица имеет вид:

. (3.8)

 

Порядок составления определителя Гурвица следующий. В качестве элемента первого столбца первой строки определителя записывается коэффици­ент a_{n- 1}, а затем на главной диагонали располагаются коэффициенты характеристического уравнения (3.4) с последовательно убывающими индексами. При этом в последнем столбце последней строки определителя записывается коэффициент .

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

Далее необходимо вычислить значение определителя Гурвица и всех его диагональных миноров, которые получают из определителя (3.8) путем отчеркивания равного числа строк и столбцов в левом верхнем углу определителя. Например, диагональный минор первого порядка равен Δ₁ = a_{n- 1}; диагональный минор второго порядка:

Δ₂ = , (3.9)

а диагональный минор третьего порядка:

Δ₃ = . (3.10)

 

 

Очевидно, что диагональный минор n -го порядка совпадает с определителем Гурвица.

Линейная система устойчива, если при выполнении необходимого условия (3.6), определитель Гурвица и все его диагональные миноры будут положительны:

Δ₁ > 0, Δ₂ > 0, …………,Δ_n > 0/

Раскрыв определитель Гурвица по последнему столбцу, получим

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

Для системы второго порядка необходимое и достаточное условия устойчивости совпадают, так как для нее Δ₁ = a₁ > 0.

Для системы третьего порядка:

Δ₁ = a₂ > 0; Δ₂ = а₁ ∙а₂ – а₀ ∙а₃.

Следовательно, достаточное условие устойчивости:

а₁ ∙а₂ > а₀ ∙а₃. (3.11)

Если определитель Δ_n = 0, то САУ находится на границе устойчивости. При этом возможны два случая:

1) сво­бодный член характеристического уравнения равен нулю, т.е. a₀ = 0, что соответствует нейтрально устойчивой системе;

2) диагональный минор Δ_(n-1) = 0, что соответствует колебательной границе устойчивости.

Из условия Δ_(n-1) = 0 можно опреде­лить параметры, при которых САУ находится на границе устойчивости.

Рассмотрим в качестве примера замкнутую систему третьего порядка с единичной отрицательной обратной связью и передаточной функцией в разомкнутом состоянии:

.

Передаточная функция такой системы в замкнутом состоянии в соответствие с выражением (2.63) равна:

.

Характеристическое уравнение такой системы:

T₁T₂T₃ p³ + (T₁T_{2 +} T₂T_{3 +} T₁T₃)р² + (T_{1 +} T_{2 +} T₃)р + (k + 1 ) = 0.

Следовательно:

= T₁T₂T₃ , T₁T_{2 +} T₂T_{3 +} T₁T₃ , а₁ = T_{1 +} T_{2 +} T3, а₀ = k + 1. (3.12)

В соответствии с выражениями (3.6) и (3.11) необходимое и достаточное условия устойчивости для системы третьего порядка определяются следующими неравенствами:

T₁T₂T₃ > 0; T₁T_{2 +} T₂T_{3 +} T₁T₃ > 0; T_{1 +} T_{2 +} T3 > 0; k + 1 > 0;

(T₁T_{2 +} T₂T_{3 +} T₁T₃)∙(T_{1 +} T_{2 +} T₃) > T₁T₂T₃ ∙(k + 1 ).

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

k < (T_{1 +} T_{2 +} T₃ )∙( 1 / T₁ + 1 / T₂ + 1 / T₃) – 1.

Предельное значение k, при котором система будет находиться на границе устойчивости, – критическое значение для рассматриваемой системы равно:

k_кр = (T_{1 +} T_{2 +} T₃ )∙( 1 / T₁ + 1 / T₂ + 1 / T₃) – 1.

Существенным недостатком критерия Гурвица является значительное усложнение условий устойчивости по мере увеличения порядка системы. Кроме того, для САУ высокого порядка достаточно трудно оценить влияние отдельных параметров звеньев на устойчивость системы. Это связано с тем, что, как правило, одни и те же параметры одновременно входят в несколько коэффициентов характеристического уравнения (смотри, например, выражение (3.12)).

3.3. Частотные критерии устойчивости Михайлова
и Найквиста

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

(3.13)

в соответствие с теоремой Безу представим в виде:

F(p) = , (3.14)

где p_i – полюс передаточной функции замкнутой системы ().

Поставляя в выражение (3.13) вместо p комплексную переменную jω, получим:

. (3.15)

После аналогичной подстановки в выражение (3.14) получим:

F(jω) = . (3.16)

Каждому сомножителю в выражении (3.16) на комплексной плоскости соответствует некоторый вектор, положение которого меняется при изменении ω.

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

F(jω) = ∑ (jω – p_i).

Если корень характеристического уравнения p_i действительный и отрицательный, т.е. рас­положен на действительной оси слева от начала координат (рис. 3.4, а), то вектор поворачивается против часовой стрелки на угол π/2, если этот ко­рень действительный и положительный (рис. 3.4, б), то вектор поворачивается по часовой стрелке на угол π/2. Следовательно, для левого действительного полюса

= π/2,

а для правого действительного полюса

= – π/2,

Нетрудно показать, что для пары комплексно сопряженных левых полюсов
(рис. 3.4, в) изменение аргумента составляет +π, а для пары комплексно сопряженных правых полюсов (рис. 3.4, г) равно -π.

Если среди n полюсов передаточной функции замкнутой системы m расположены справа от мнимой оси, а остальные (n – m)– слева, то изменение аргумента комплексной функции F(jω) вектора равно:

F(jω) = (n – m)∙ π/2 – m ∙π/2 = (n – 2 m)∙π/2. (3.17)

Выражение (3.17) и определяет суть принципа аргумента. В передаточной функции устойчи­вой системы правые полюса отсутствуют, т.е. m = 0, и изменение аргумента F(jω) равно:

F(jω) = n ∙π/2. (3.18)

 

Из выражения (3.18) следует критерий устойчивости Михайлова, согласно которому из­менение аргумента комплексной функции F(jω) определяется по годографу, который записывают в виде

F(jω) = Х(ω) + jY(ω),

где Х(ω),Y(ω) – действительная и мнимая составляющие комплексной функции F(jω);

Х(ω) = a ₀ – a ₂ ω² + a ₄ ω⁴ – a ₆ ω⁶ + ……;

Y(ω) = a ₁ ω – a ₃ ω³ + a ₅ ω⁵ – a ₇ ω⁷ + …….

Каждому значению ω = ω_i на комплексной плоскости соответствует точка с координатами (Х(ω_i),Y(ω_i)). При изменении ω эта точка описывает на плоскости некоторую траекторию, которая называется годографом Михайлова (рис. 3.5).

При ω = 0 Х(0) = a₀, Y(0) = 0, т.е. F(0) = a₀, причем в соответствии с выражением (3.6) a₀ > 0.

Формулировка критерия Михайлова: замкнутая система устойчива, если годограф F(jω), начинаясь при ω = 0 на положительной дей­ствительной полуоси, при изменении ω от 0 до обходит последовательно в положитель­ном направлении (против часовой стрелки) n квадрантов, где n – порядок системы.

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

Из формулировки критерия следует, что система устойчива, если годограф Михайлова, начавшись на действительной оси при ω = 0, несколько раз последовательно пересекает действительную и мнимую ось. Значения ω, при которых происходят эти пересечения, являются действительными положительными корнями уравнения Y(ω) = 0 (при пересечении с действительной осью) и уравнения Х(ω) = 0 (при пересечении с мнимой осью). Следовательно, оценить устойчивость системы можно и без построения годографа: достаточно, чтобы корни указанных уравнений чередовались друг с другом.

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

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

W(p) =

имеет n -й порядок, т.е. число ее полюсов (порядок полинома A(р))равно n. На основании принципа физической реализуемости можно утверждать, что число нулей передаточной функции W(p) (порядок полинома В(р))не превышает n. Передаточная функция системы в замкнутом состоянии:

Ф(р) = = .

Введем в рассмотрение выражение:

D(p) = 1 + W(p) = . (3.19)

Очевидно, что число нулей и полюсов выражения D(p) одинаково и равно n. При этом числитель выражения (3.19) является характеристическим полиномом замкнутой системы, а знаменатель – характерис­тическим полиномом разомкнутой системы. Осуществим в выражении (3.19) замену p на jω:

D(jω) = (3.18)

и определим из­менение аргумента (3.20), полагая, что замкнутая система устойчива. Поскольку в этом случае в соответствии с принципом аргумента

В(jω) = n ∙π/2 и (А(jω) + В(jω)) = n ∙π/2,

то

D(jω) = (А(jω) + В(jω)) – В(jω) = 0.

Таким образом, если разомкнутая и замкнутая систе­мы устойчивы, то изменение аргумента D(jω) равно нулю, следовательно, годограф D(jω) не охватывает начала координат (рис. 3.6, а).

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

Очевидно, что об изменении аргумента вектора удобнее судить не по годографу D(jω), а погодографу амплитудно-фазовойхарактеристики разомкнутой системы W(jω). Поскольку D(jω) = 1 + W(jω), изменение аргумента D(jω) будет рав­но нулю, если годограф амплитудно-фазовойхаракте­ристики разомкнутой системы не охватывает точку с ко­ординатами (-1, j0) (рис. 3.6, б).

Отсюда следует формулировка критерия Найквиста: система, устойчивая в разомкнутом состояние, будет устойчива и в замкнутом состоянии, если годограф амплитудно-фазовой характеристики разомкнутой системы не охва­тывает точку с координатами (-1, j0). В том случае, когда годограф частотной характеристики охватывает эту точку, система неустойчива.

Если САУ содержит ν интегрирующих звеньев, то начальное значение фазочастотной характеристики разомкнутой системы равно (-ν∙π/2), а амплитудно-частотной – бесконечности. Поэтому,

· если ν =1, характеристика W(jω) при ω → 0уходит в бесконечность вдоль отрицательной мнимой полуоси;

· если ν = 2 – вдоль отрицательной действительной полуоси;

· если ν = 3 – вдоль положительной мнимой полуоси.

Для удобства оценки устойчивости таких астатических систем годограф W(jω) дополняют дугой бесконечного радиуса, начинающейся на положительной действительной полуоси и проводимой до пересечения с годографом W(jω) (рис. 3.7). Формулировка критерия устойчивости при этом не изменяется.

Если годограф W(jω) разомкнутой системы проходят через точку (1, j0), то система в замкну­том состоянии находится на границе устойчивости.

Рассмотрим случай, когда система в разомкнутом состоянии неустойчива и имеет m правых полюсов. Полагая, что при замыкании обратной связи система становится устойчивой, в соответствии с принципом аргумента получаем:

В(jω) = (n –2m) ∙π/2 и (А(jω) + В(jω)) = n ∙π/2,

При этом изменение аргумента D(jω) равно:

D(jω) = (А(jω) + В(jω)) – В(jω) = n ∙π/2 – (n – 2 m) ∙π/2 = m ∙π = 2π∙ m /2.

Следовательно, система в замк­нутом состоянии будет устойчивой, если годограф час­тотной характеристики D(jω) m /2раз охватывает начало координат, а соответственно годограф амплитудно-фазовойхаракте­ристики разомкнутой системы W(jω) m /2раз охва­тывает точку с координатами (-1,j0), где m – число правых полюсов разомкнутой системы.

Годограф амплитудно-фазовойхарактеристики W(jω) реальной технической системы может иметь достаточно сложную форму (рис. 3.8).

В этом случае сложно определить, сколько раз годограф W(jω) охватывает начало координат. Задача упрощается, если ввести в рассмотрение понятие перехода годографа W(jω ) через действительную ось, т.е. пересечение графиком W(jω) действительной оси левее точки с координатами (-1, j0). Перехода годографа W(jω) через действительную ось считается положительным, если при увеличении частоты ω пересечение оси происходит сверху вниз (годограф переходит из второго квадранта в третий), в противном случае переход считается отрицательным.

Обозначим число положительных переходов m₊, а число отрицательных переходов m_- . Тогда критерий устойчивости Найквиста может быть сформулирован так: система в замк­нутом состоянии становится устойчивой, если разность между числом положительных и отрицательных переходов равна m/ 2, т.е.

m₊ – m_- = m/ 2, (3.21)

где m – число правых полюсов разомкнутой системы.

Если система в разомкнутом состоянии устойчива (m = 0) условие устойчивости системы при ее замыкании упрощается:

m₊ – m_- = 0. (3.22)

Запасы устойчивости

В процессе эксплуатации САУ ее параметры (коэффициенты усиления, постоянные времени) из-за из­менения внешних условий, колебаний напряжений источ­ников энергии и других причин отличаются от расчетных значений. Если не принять определенных мер, то исходная устойчивая система может стать неустойчивой. Для исключения этого явления при проектировании следует обеспечить опреде­ленные запасы устойчивости системы, которые характе­ризуют близость годографа амплитудно-фазовойхаракте­ристики разомкнутой системы W(jω) к точке с координатами (-1, j0).

Различают запас устойчивости по фазе и усилению. Запасы устойчивости определяются на двух частотах: частоте среза ω_с и критической частоте ω_кр. На частоте среза амплитудно-частотная характеристика разомкнутой системы | W(jω)| равна единице, а на критичес­кой частоте фазо-частотная характеристика этой системы φ(ω) принимает значение, равное -π.

Запас устойчивости по фазе Δφ показывает, насколько фазо-частотная характеристика разомкнутой системы на частоте среза ω_с отлича­ется от -π (рис. 3.9):

Δφ = π – .

Величина запаса устойчивости по усилению может быть определена на частоте ω_кр, как разность:

= 1 – | W(jω_кр)|,

либо как отношение

α = 1/ | W(jω_кр)|.

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

Системы, годографы W(jω) ко­торых пересекают вещественную ось только справа от точки с координатами (-1, j0)(рис. 3.10, а), называют аб­солютно устойчивыми. В таких системах неустойчивость может наступить только при увеличении коэффициента усиления.

Если годограф частотной характеристики W(jω) разомкну­той системы пересекает вещественную ось и слева от точ­ки с координатами (-1, j0), то систему называют услов­но устойчивой (рис. 3.10, б). Неустойчивой такая система может быть как при увеличении, так и при уменьшении коэффициента усиления.

Для нормальной работы САУ необходимо, что­бы запас устойчивости по усилению α был не менее двух, а запас устойчивости по фа­зе – от 0,5 до 1 рад.

3.5. Оценка устойчивости по логарифмическим
амплитудно- и фазо-частотным
характеристикам

Оценку устойчивости замкнутой САУ можно осуществлять по логарифмическим амплитудно- и фазо-частотным характеристикам системы в разомкнутом состоянии: L(ω) и φ(ω). В том случае, когда годограф W(jω) не имеет точек пересечения с вещественной осью слева от точки с коор­динатами (-1, j0), для устойчивости замкнутой сис­темы необходимо и достаточно, чтобы выполнялось условие

ω_с < ω_кр.

На рис. 3.11 приведены ЛАХ и ФЧХ устойчивых разомкнутых системы, одна из которых при замыкании остается устойчивой (рис. 3.11, а), а другая – становится неустойчивой (рис. 3.11, б).

По L(ω) и φ(ω) разомкнутой системы можно определить за­пасы устойчивости: запас по фазе Δφ отсчитывают по фазо-частотной характеристике на частоте среза ω_с, а запас устойчивости по усилению ΔL равен зна­чению ЛАХ на критической частоте ω_кр, взятому с обратным знаком, т.е. ΔL = | L(ω_кр)| (см. рис. 3.11, а).

Если ω_с = ω_кр, то система находится на границе устойчивости.

Если при некотором значении коэффициента усиления (k) замкнутая система устойчива с запасом устойчивости по усилению равным ΔL, то величина критического коэффициента усиления k_кр может быть вычислена по формуле:

20lg k_кр = 20lg k + ΔL.

Для оценки устойчивости условно устойчивых САУ реальных технических систем, имеющих обычно достаточно сложную форму, также можно воспользоваться понятием перехода. При этом переходом называется пересечение графика φ(ω) с горизонтальной прямой -π, при условии, что на частоте, при которой φ(ω) = – π, ЛАХ положительна.

Правило определения знака перехода противоположно рассмотренному для W(jω): переход графика φ(ω) через уровень(- π) считается положительным, если при увеличении частоты ω пересечение этого уровня происходит снизувверх, в противном случае переход считается отрицательным. Обозначим число положительных переходов m₊, а число отрицательных переходов m_- .

В этом случае формулировка критерия устойчивости Найквиста: система в замк­нутом состоянии становится устойчивой, если разность между числом положительных и отрицательных переходов равна m/ 2, т.е.

m₊ – m_- = m/ 2, (3.23)

где m – число правых полюсов разомкнутой системы.

Если число положительных переходов φ(ω) равно числу отрицательных, то система, устойчивая в разомкнутом состоянии , остается устойчивой при замыкании. На рис. 3.12 в качестве примера приведены логарифмические амплитудно- и фазо-частотные характеристики неустойчивой разомкнутой системы, имеющей два правых полюса . При замыкании такая система становится устойчивой, так как m₊ = 1, а m_- = 0, и условие (3.23) выполняется.

Рекомендации для обеспечения запаса устойчивости, которые следуют из практики проектирования САУ:

· во-первых, для того чтобы в системе были обеспечены необ­ходимые запасы устойчивости, наклон ЛАХ в диапазоне частот, в котором расположена частота среза, должен быть равен -20 дБ/дек. Если в указанном частотном диапазоне наклон L(ω) равен -40 дБ/дек, обеспечить необходи­мый запас устойчивости по фазе затруднительно. При наклоне 0 дБ/дек система обладает чрезмерно большим запасом устойчивости по фа­зе и становится передемпфированной с длитель­ным переходным процессом:

· во-вторых, запас устойчиво­сти системы по фазе зависит от диапазона частот, в котором ЛАХ разомкнутой системы в области частоты среза имеет наклон -20 дБ/дек. Чем шире этот диапазон частот, тем выше запас устойчивости по фазе и наоборот.