Критерий рауса гурвица по характеристическому уравнению
8.1. Понятие устойчивости системы
Под устойчивостью системы понимается способность ее возвращаться к состоянию установившегося равновесия после снятия возмущения, нарушившего это равновесие. Неустойчивая система непрерывно удаляется от равновесного состояния или совершает вокруг него колебания с возрастающей амплитудой.
Устойчивость линейной системы определяется не характером возмущения, а структурой самой системы (рис.61). Говорят, что система устойчива «в малом» , если определен факт наличия устойчивости, но не определены ее границы. Система устойчива «в большом» , когда определены границы устойчивости и то, что реальные отклонения не выходят за эти границы.
В соответствии с классическим методом решение дифференциального уравнения ищется в виде:
y(t) = y вын (t) + y св (t).
Здесь yсв(t) — общее решение однородного дифференциального уравнения , то есть уравнения с нулевой правой частью:
a o y (n) + a 1 y (n-1) + . + a (n-1) y’ + a (n) y = 0.
Физически это означает, что все внешние воздействия сняты и система абсолютно свободна, ее движения определяются лишь собственной структурой. Поэтому решение данного уравнения называется свободной составляющей общего решения. y вын (t) — частное решение неоднородного дифференциального уравнения , под которым понимается уравнение с ненулевой правой частью. Физически это означает, что к системе приложено внешнее воздействие u(t) . Поэтому вторая составляющая общего решения называется вынужденный . Она определяет вынужденный установившийся режим работы системы после окончания переходного процесса.
Можно провести аналогию между САУ и пружиной, колебания которой описываются аналогичным дифференциальным уравнением (рис.62). Оттянем пружину, а затем отпустим, предоставив ее самой себе. Пружина будет колебаться в соответствии со свободной составляющей решения уравнения, то есть характер колебаний будет определяться только структурой самой пружины. Если в момент времени t = 0 подвесить к пружине груз, то на свободные колебания наложится внешняя сила Р . После затухания колебаний, описываемых только свободной составляющей общего решения, система перейдет в новый установившийся режим, характеризуемый вынужденной составляющей y вын = y(t ) . Если внешнее воздействие само будет изменяться по синусоидальному закону P = P o sin(t + ) , то после затухания переходного процесса система будет совершать вынужденные колебания с той же частотой, что и вынуждающая сила, то есть y вын = y max sin(t + y).
Каждая составляющая общего решения уравнения динамики ищется отдельно. Вынужденная составляющая ищется на основе решения уравнения статики для данной системы для времени t . Свободная составляющая представляет собой сумму из n отдельных составляющих: , где p i корни характеристического уравнения D(p) = a 0 p n + a 1 p n -1 + a 2 p n -2 + . + a n = 0 . Корни могут быть либо вещественными p i = a i , либо попарно комплексно сопряженными p i = a i ± ji . Постоянные интегрирования А i определяются исходя из начальных и конечных условий, подставляя в общее решение значения u, y и их производные в моменты времени t = 0 и t .
Каждому отрицательному вещественному корню соответствует экспоненциально затухающая во времени составляющая y св (t) i , каждому положительному — экспоненциально расходящаяся, каждому нулевому корню соответствует y св (t) i = const (рис.63). Пара комплексно сопряженных корней с отрицательной вещественной частью определяет затухающие колебания с частотой i , при положительной вещественной части — расходящиеся колебания, при нулевой — незатухающие (рис.64).
Так как после снятия возмущения y вын (t) = 0 , то устойчивость системы определяется только характером свободной составляющей y св (t) . zПоэтому условие устойчивости систем по Ляпунову формулируется так: в устойчивой системе свободная составляющая решения уравнения динамики, записанному в отклонениях, должна стремиться к нулю, то есть затухать.
Исходя из расположения на комплексной плоскости корни с отрицательными вещественными частями называются левыми , с положительными — правыми (рис.65).
Поэтому условие устойчивости линейной САУ можно сформулировать следующим образом: для того, чтобы система была устойчива, необходимо и достаточно, чтобы все корни ее характеристического уравнения были левыми. Если хотя бы один корень правый, то система неустойчива. Если один из корней равен нулю (в системах, где a n = 0 ), а остальные левые, то система находится на границе апериодической устойчивости . Если равны нулю вещественные части одной или нескольких пар комплексно сопряженных корней, то система находится на границе колебательной устойчивости .
Правила, позволяющие судить о знаках корней характеристического уравнения без его решения, называются критериями устойчивости . Их можно разделить на алгебраические (основаны на составлении по данному характеристическому уравнению по определенным правилам алгебраических выражений, по которым можно судить об устойчивости САУ) и частотные (основаны на исследовании частотных характеристик).
8.2. Алгебраические критерии устойчивости
8.2.1. Необходимое условие устойчивости
Характеристическое уравнение системы с помощью теоремы Виета может быть записано в виде
D(p) = a o p n + a 1 p n-1 + a 2 p n-2 + . + a n = a o (p-p 1 )(p-p 2 ). (p-p n ) = 0,
где p 1 , p 2 , . p n — корни этого уравнения. Если система устойчива, значит все корни левые, то есть вещественные части всех корней
отрицательны, что можно записать как a i = -|a i | . Подставим их в уравнение:
a 0 (p + |a 1 |)(p + |a 2 | — j2)(p + |a 2 | + j2). = 0.
Перемножая комплексно сопряженные выражения, получим:
a 0 (p + |a 1 |)((p + |a 2 |)2 + (2)2). = 0.
После раскрытия скобок должно получиться выражение
a 0 p n + a 1 p n-1 + a 2 p n-2 + . + a n = 0.
Так как в скобках нет ни одного отрицательного числа, то ни один из коэффициентов a 0 ,a 1 . a n не будет отрицательным. Поэтому необходимым условием устойчивости САУ является положительность всех коэффициентов характеристического уравнения: a 0 > 0, a 1 > 0, . , a n > 0 . В дальнейшем будем рассматривать только уравнения, где a 0 > 0 . В противном случае уравнение домножается на -1.
Рассмотренное условие является необходиным, но не достаточным условием. Необходимые и достаточные условия дают алгебраические критерии Рауса и Гурвица.
8.2.1. Критерий Рауса
Раус предложил критерий устойчивости САУ в виде алгоритма, по которому заполняется специальная таблица с использованием коэффициентов характеристического уравнения:
1) в первой строке записываются коэффициенты уравнения с четными индексами в порядке их возрастания;
2) во второй строке — с нечетными;
3) остальные элементы таблицы определяется по формуле: c k,i = c k+ 1,i — 2 — ric k + 1,i — 1 , где ri = c 1,i — 2 /c 1,i — 1 , i 3 — номер строки, k — номер столбца.
4) Число строк таблицы Рауса на единицу больше порядка характеристического уравнения.
Критерии устойчивости Рауса–Гурвица и Михайлова
(геометрический критерий устойчивости)
Пусть имеем линейное дифференциальное уравнение с постоянными вещественными коэффициентами:
Нулевое решение уравнения (1) асимптотически устойчиво, если все корни характеристического уравнения
имеют отрицательные вещественные части.
Критерий Рауса—Гурвица . Для того чтобы все корни уравнения (2) имели отрицательные вещественные части, необходимо и достаточно, чтобы были положительными все главные диагональные миноры матрицы Гурвица
Матрица Гурвица составляется так. По главной диагонали выписываются коэффициенты многочлена (2), начиная с и оканчивая . Столбцы состоят поочередно из коэффициентов только с нечетными или только с четными индексами, причем в число последних включается коэффициент . Все остальные элементы матрицы, отвечающие коэффициентам с индексами, большими или меньшими , полагаются равными нулю. Главные диагональные миноры матрицы Гурвица имеют вид
Таким образом, условие Гурвица гласит: для устойчивости решения уравнения (1) необходимо и достаточно, чтобы выполнялись соотношения
Так как , то условие 0″ png;base64,iVBORw0KGgoAAAANSUhEUgAAAEAAAAAUCAMAAAAKqMsNAAAAM1BMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADbQS4qAAAAEHRSTlMAEBeBQeekAcBcIdGCMHGRAmkNMgAAAQ1JREFUOMutVNuygyAMlIHcuCj//7UnWBGhasee5gEHJJvNZodp+lW4f+ZLlocZxvnjNmSgR/kRLLJpDTA8oyCg9THUosROLinIbN4PAUvZvahjomsKEuwI4bMtMuRl27POIF6rQBLQ9zPLsZxn3gko3J0KJNxBxJWBAtCmQPnY+0E4xuYVuwFAI/CBwnoPOQ0MXmNgqbB86wXfAKRqUBJcZa4UCslk1GVnLbDb8dPKIOUyzAmFtohFVbEwzwsOIrpeRDr4wGeuEbJSQK+Lg8EJwxjVyWoNLMuE+RiFQlA5sHNB9G8dYZAFVkmoD/2l9grSUiSaM0nTLOZcbBBty+D3T4M25sGmr/PJPLr+B9cTCgo80fp7AAAAAElFTkSuQmCC» style=»vertical-align: middle;» /> может быть заменено требованием 0″ png;base64,iVBORw0KGgoAAAANSUhEUgAAADoAAAATBAMAAAAzNgPlAAAALVBMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACttl6nAAAADnRSTlMAhgFdwNcQoTBBIfCwcRgyobwAAADiSURBVCjPY2AgD7DNwicrKuKGW1J5G8MRBTQx4wQYS7WAYd4EdB2BMOk+ARBiYBMpM0OWXgBh2Akw6DkwMKQWdb9A1u0OkZYDyjoxMIUwsGxCNpzTvQEmG8DA+5hB+wCK3ZxHGuB6VZ4x9F1AkeUqaYDbqxfEcK+BoY0N7npOFzBzngDDvAKQgnOKk6XvCkIlNaGu0gb614CBJ4D5CaviuoYpEEnmcKiPGLcxlCowcIUftD7MYMogDpGEhwZDqWEEyBEKbA0MhQzBYD2GcEkGDmE4O5jjWQHuGNnGFd5AcgoAAPGxMydfDSoLAAAAAElFTkSuQmCC» style=»vertical-align: middle;» />.
Пример 1. Исследовать на устойчивость нулевое решение уравнения
Решение. Составляем характеристическое уравнение
Здесь . Выписываем диагональные миноры Гурвица
\Delta_4=\begin
\Delta_4>0″ png;base64,iVBORw0KGgoAAAANSUhEUgAAASAAAAAUCAMAAAAawRC9AAAAM1BMVEVHcEwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADbQS4qAAAAEHRSTlMAYaQhwNaCBD4QkfCBxHFQI3bn4gAAAjJJREFUWMPtmNuSgyAMhg0BGg4jvP/TLlTaooJIZTp7yoWzu5Jvwp8Qsk7Tv50zUv+oI0PPOh0I/xSKeY8d3jcmpNDlnPwKlNqg0PAuhbUID6GLyn9PlLiGYpJqCmso6OtleM6+lKwDlK2iZC9K96LmgFJvR4VGTTWFlRWwJYOn2PY8lLIeUXIcSn8iqqmCco98yPibqfR5NQu7fiWTvyxl/QeiZAOlRHzNva0eSRA2L3Wd/PfZVeYESq2aZxXVGxVrRiXfjEreHV8KoyoR2CsYnvz5btkGRQ6HoYhGocKP2INaUvVUGJih4qTFmdoIzMtZf+ZdAxh7jOInUXam4v30BipWFPWgZKq8pDBOviWQrBXzGoWCUru7jJrCAuvfjMquUeEPgs6hZK7v65CW9uR4VoFz8p8rqQojREShcWHRDgZ9qCUqR8tscj2qsAT4foeFqNyCks/W5ZLChT3lBzTO7dHHxnUIhazHezIlC7aDDAhdQmELhW6X9gqqFZViEz+Jwvt5FPxhfjmkG4EUCLlp25w/HiYbOwootZ4xLqAQuBuACjURrr61QIeoUE25mZ1AYY6Sat/2GLB75TLDssLM7X4gGAxDrdrZBRQFjXKBGqiSrQQqDfXxHMu0iI7+xQnOhLezKKyjkOOquzSj0vWRJ1ZSdiO2UC2BWgZw8G4mYjQCFad+54egune4y3rQl5//LiDUwa6iDUFNWoOYx6DCfe6F/tSXSKAPoajnk+C4qL4A/qkkU6k3D1AAAAAASUVORK5CYII=» style=»vertical-align: middle;» />. Следовательно, тривиальное решение уравнения (5) асимптотически устойчиво.
Вычисление можно, например, организовать так. Составляем сначала старший минор Гурвица . По нему легко выписываются все младшие миноры . Затем начинаем вычислять последовательно и т.д. Если встретился отрицательный минор, решение неустойчиво и дальнейший подсчет не нужен.
Геометрический критерий устойчивости (критерий Михайлова)
Рассмотрим линейное дифференциальное уравнение n-го порядка с постоянными вещественными коэффициентами
Критерий рауса гурвица по характеристическому уравнению
где A k , k +1 и B k , k +1 – постоянные интегрирования.
Из данного выражения следует, что только при отрицательном значении отклонение регулируемой величины с течением времени будет стремиться к нулю. Теперь можно сформулировать необходимое и достаточное условие устойчивости линейных систем:
Чтобы линейная система была устойчивой, необходимо и достаточно, чтобы вещественные корни и вещественные части комплексных корней были отрицательными.
Устойчивость является одним из главных требований к автоматическим системам регулирования (АСР). Если система не способна восстановить равновесное состояние, нарушенное в процессе работы, то она непригодна для практического использования. После того, как вычислены передаточная функция объекта и параметры настройки регулятора, необходимо определить, устойчива система или нет. Для проверки можно использовать как необходимое и достаточное условие устойчивости, приведенное выше, либо критерии устойчивости. Далее рассмотрим три таких критерия:
— критерий Рауса – Гурвица;
— критерий Михайлова;
— критерий Найквиста – Михайлова.
7.2. Критерии устойчивости
Использование для проверки устойчивости системы необходимого и достаточного условия устойчивости связано с решением характеристического уравнения, что может представлять трудности в случае систем высокого порядка, а если рассматривается замкнутая система с запаздыванием, то и вовсе невозможно. Значительно проще можно выполнить проверку устойчивости с помощью специальных критериев.
Критерий Рауса – Гурвица
Критерий Рауса – Гурвица позволяет судить об устойчивости системы по коэффициентам ее характеристического уравнения. Впервые критерий устойчивости был предложен в 1875 году английским математиком Раусом в виде таблицы. В 1895 году швейцарский математик Гурвиц опубликовал критерий устойчивости в виде определителей. Так как оба этих критерия приводят к одним и тем же алгебраическим неравенствам и отличаются только способом их получения, указанные критерии объединяют, называя критерием Рауса – Гурвица.
Пусть имеется характеристическое уравнение системы:
Если хотя бы один из определителей меньше нуля, система будет неустойчивой, поэтому после обнаружения первого отрицательного определителя процесс проверки устойчивости прекращается .
Если отрицательных определителей нет, но есть равные нулю, такая система находится на границе устойчивости (нейтральна).
Пример 1. Характеристическое уравнение имеет вид
Δ 3 вычислять не нужно, т.к. уже второй определитель оказался отрицательным.
Пример 2. Характеристическое уравнение имеет вид
Критерий Михайлова
В 1936 году советским ученым А.В.Михайловым был предложен критерий устойчивости, основанный на взаимосвязи между характером переходных процессов, возникающих при нарушении равновесия системы и амплитудной и фазовой вынужденных колебаний, установившихся в системе под воздействием гармонически изменяющейся входной величины. Анализ устойчивости в соответствии с критерием Михайлова сводится к построению по характеристическому уравнению системы графика, носящего название годограф Михайлова, по виду которого можно судить о состоянии системы (устойчива она или нет).
Годограф строится следующим образом. Пусть характеристическое уравнение имеет вид (7.5).
Обозначим через F(p) многочлен, стоящий в левой части этого уравнения:
Рис. 7.3. Примеры годографов Михайлова для различных систем
Рис. 7.4. Построение годографа Михайлова в Mathcad:
а – для значений ω = 0…120; б – для значений ω = 0…1600
Критерий Найквиста – Михайлова
Критерии Рауса – Гурвица и Михайлова не подходят для замкнутых систем с запаздыванием, так как в их характеристическое уравнение входит экспонента. В 1932 году американский ученый Найквист предложил критерий для исследования устойчивости усилителей с обратной связью. В 1938 году Михайлов обобщил и распространил его на системы автоматического регулирования. Этот критерий позволяет судить об устойчивости замкнутой системы с отрицательной обратной связью (в том числе, автоматической системы регулирования, АСР) по амплитудно-фазовой характеристике (АФХ) разомкнутой системы.
Любая АСР, как система с обратной связью, является замкнутой. Изменение сигнала на входе какого-нибудь элемента (объекта, регулятора) сказывается на выходных сигналах всех других элементов. Они также изменяются в соответствии со своими динамическими свойствами. Если же нарушить связь между какими-либо двумя элементами, например между чувствительным элементом и элементом сравнения (разорвать обратную связь), то система из замкнутой превратится в разомкнутую.
Рассмотрим разомкнутую систему, состоящую из последовательно соединенных регулятора и объекта регулирования (рис. 7.5).
Если обозначить через W об(p) передаточную функцию объекта, а через W p(p) – передаточную функцию регулятора, передаточная функция и амплитудно-фазовая характеристика разомкнутой системы будут иметь вид:
Подадим от генератора синусоидальных колебаний на вход разомкнутой системы сигнал xвх с частотой ω и амплитудой A вх. Спустя некоторое время на выходе также установятся синусоидальные колебания xвых с той же частотой, амплитудой A вых и сдвигом фазы φ. Предположим, что при некотором значении частоты сдвиг фазы оказался равен +π или –π, а амплитуды колебаний на входе и на выходе оказались равны между собой:
Данное состояние разомкнутой системы соответствует такой точке ее АФХ, которая находится на единичном расстоянии от начала координат (модуль АФХ M раз(ω) = A вых/A вх = 1) с углом π от вещественной положительной полуоси (аргумент АФХ φ раз(ω) = ±π). На комплексной плоскости эта точка имеет координаты (–1, 0i).
Теперь отключим генератор колебаний от входа системы и одновременно замкнем отрицательную обратную связь, как это показано на рис. 7.6. В результате мы получим замкнутую систему, представляющую собой одноконтурную АСР.
Рис. 7.6. Замкнутая и разомкнутая системы
В результате получим замкнутую систему, в которой будут наблюдаться незатухающие колебания с постоянной амплитудой, а, как мы знаем, это характерно для систем, находящихся на границе устойчивости (см. рис. 7.2,е в разделе 7.1). Таким образом, прохождение графика АФХ разомкнутой системы через точку (–1, 0i) свидетельствует о нахождении соответствующей замкнутой системы на границе устойчивости.
Теперь рассмотрим второй случай. Предположим, в ходе того же самого опыта на выходе разомкнутой системы установились колебания с тем же сдвигом фазы π, но с меньшей амплитудой:
Это состояние соответствует точке АФХ разомкнутой системы, лежащей на вещественной оси правее точки (–1, 0i) (M раз(ω) = A вых/A вх раз(ω) = ±π).
После отключения генератора колебаний и замыкания отрицательной обратной связи на вход W раз(iω) поступит сигнал с меньшей амплитудой, чем создавал генератор. Пройдя через регулятор и объект, этот сигнал снова будет несколько ослаблен (так как M раз(ω) раз(iω) с еще меньшей амплитудой, и это будет повторяться снова и снова. Таким образом, в данном случае с течением времени амплитуда колебаний в замкнутой системе будет непрерывно уменьшаться, стремясь к нулю, т.е. будет наблюдаться затухающий процесс, что, как известно, характерно для устойчивых систем (рис. 7.2,г).
В третьем возможном случае на выходе разомкнутой системы установятся колебания со сдвигом фазы π и большей амплитудой, нежели входные колебания:
что соответствует точке АФХ разомкнутой системы, лежащей на вещественной оси левее точки (–1, 0i) (M раз(ω) = A вых/A вх > 1, φ раз(ω) = ±π).
Здесь при прохождении через регулятор и объект сигнал усиливается (M раз(ω) > 1), и его подача по обратной связи на вход Wраз(iω) приводит к возникновению колебательного процесса с непрерывно нарастающей амплитудой («расходящегося»), как на рис. 7.2,д, т.е. замкнутая система оказывается неустойчивой.
Из сказанного становится ясно, что об устойчивости замкнутой системы можно судить на основании расположения амплитудно-фазовой характеристики разомкнутой системы относительно точки с координатами (–1, 0i) (рис. 7.7).
Критерий Найквиста – Михайлова формулируется следующим образом.
Замкнутая система устойчива, если она устойчива в разомкнутом состоянии, и ее амплитудно-фазовая характеристика (построенная для всех значений ω от 0 до бесконечности ) не охватывает точку с координатами (–1, 0 i ) .
Если амплитудно-фазовая характеристика проходит через точку с координатами (–1, 0i), система находится на границе устойчивости. Если АФХ охватывает точку с координатами (–1, 0i), то система неустойчива.
Рис. 7.7. Примеры АФХ для различных разомкнутых систем
http://mathhelpplanet.com/static.php?p=kriterii-ustoichivosti-gurvitsa-i-mihailova
http://project2425235.tilda.ws/glava7