Критерий устойчивости Гурвица

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

по главной диагонали выписываются все коэффициенты характеристического уравнения от а₁ до а_n в порядке возрастания индексов;

столбцы определителя заполняются коэффициентами от главной диагонали вниз по убывающим, а вверх- по возрастающим индексам;

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

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

где все коэффициенты строго больше нуля. Получим

.

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

Для устойчивости системы 5-го порядка необходимо выполнение условий

а_k>0, k=0,1,2,...5;

₂ =а₁а₂ - а₀а₃>0;

₃=а₃₂ - а₁²а₄>0;

₄ =а₄₃ -а₂а₅₂ + а₀а₅(а₁а₄ - а₀а₅)>0;

₅ =а₅₄>0.

Так как при выполнении необходимого условия устойчивости всегда а_n>0, то об устойчивости системы можно судить по определителям до _n-1 включительно. Доказано, что если _n-1=0, то система находится на колебательной границе устойчивости, т.е. имеет пару чисто мнимых корней. Из условия _n-1=0 можно определить критические значения параметров системы, при которых она выходит на границу устойчивости.

Пример 6.1.1

Исследовать устойчивость системы стабилизации угла тангажа самолета и определить критическое значение передаточного числа автопилота по углу тангажа. Система задана структурной схемой.

Рис.6.1.3. Структурная схема системы стабилизации угла тангажа.

На схеме обозначено:

k_- передаточное число (коэффициент передачи) автопилота по углу тангажа;

передаточная функция рулевого привода;

передаточная функция самолета по угловой скорости тангажа _z;

k_z - передаточное число автопилота по угловой скорости тангажа.

Для передаточной функции разомкнутой системы можно записать

где

Передаточная функция замкнутой системы примет вид

где

Составим определитель Гурвица

Оценим устойчивость системы для следующих значений параметров:

.

При этих значениях для коэффициентов характеристического уравнения получим

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

Условия устойчивости выполнены и система при избранных параметрах устойчива.

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

Отсюда

В последнем выражении только d₃ и d₄ являются функциями коэффициента k_ и подставив их в него, получим квадратное уравнение относительно этого коэффициента

Решив это уравнение, получим критическое значение передаточного числа по углу тангажа

Система устойчива, если k<16.56.