Использование Routh-Hurwitz для прогнозирования Системный ответ и переходное поведение

Введение в критерий Рут-Гурвица

Критерий Рут-Гурвица является краеугольным камнем классической теории управления, обеспечивая систематический метод оценки стабильности линейных временных инвариантных систем без явного вычисления корней характеристического уравнения. Разработанный независимо Эдвардом Джоном Рутом в 1874 году и Адольфом Гурвиц в 1895 году, этот алгебраический метод остается широко используемым в технике управления, потому что он предлагает быстрый, надежный способ определить, будет ли система реагировать стабильным, предсказуемым образом или проявлять колебания и нестабильность. Стабильность имеет основополагающее значение для проектирования системы - нестабильная система может привести к опасному перерасходу, постоянным колебаниям или даже катастрофическому отказу, особенно в таких приложениях, как аэрокосмическая промышленность, автомобильное управление, химические процессы и робототехника.

Анализируя коэффициенты характеристического полинома системы, критерий Рут-Гурвица идентифицирует число полюсов (корней уравнения характеристик), которые лежат в правой половине сложной плоскости. Поляки с положительными реальными частями вызывают неустойчивое поведение, в то время как полюса с отрицательными реальными частями способствуют стабильным, конвергентным ответам. Метод не требует решения многочленов высокой степени, что делает его практичным для ручного анализа и ранней стадии проектирования. Кроме того, его понимание переходного поведения - такого как время урегулирования, перерасход и частота колебаний - помогает инженерам настраивать контроллеры для удовлетворения спецификаций производительности, прежде чем совершать детальные симуляции или аппаратные прототипы.

В этой статье сначала будет объяснена математическая основа критерия Рут-Гурвица, затем вы будете руководствоваться построением массива Рут пошаговым примером. После этого мы рассмотрим, как критерий предсказывает системный ответ и преходящее поведение, включая ключевые характеристики, такие как перерасход и время урегулирования. Практический численный пример иллюстрирует процесс, и мы обсудим ограничения метода и дополнительные инструменты. К концу у вас будет четкое понимание того, как применять критерий Рут-Гурвиц в качестве практического, авторитетного инструмента в анализе системы управления.

Что такое критерий Рут-Хурвица?

Критерий Рут-Гурвица является необходимым и достаточным условием стабильности в системах LTI. Он работает на характеристическом уравнении системы, которое выводится из знаменателя функции переноса замкнутого цикла после установки его на ноль. Для системы, описанной функцией переноса , характеристическое уравнение является , где является многочленом в s степени n:

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

Важно отметить, что критерий Рут-Гурвица является чисто алгебраическим — он не дает фактических местоположений полюсов. Вместо этого он обеспечивает подсчет полюсов правой половины плоскости (RHP). Это делает его бесценным, когда характерный многочлен имеет высокий порядок (например, четвертый, пятый или выше), где поиск корня является вычислительно интенсивным или аналитически неосуществимым.

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

Оригинальное название: The Routh Array: Step-by-Step

Построение массива Рут является систематической процедурой. Мы проиллюстрируем это на конкретном примере полинома:

Этот полином имеет известные корни при s = -1, -2, -3, -4 (все отрицательные), поэтому система стабильна. Мы будем использовать его для проверки того, что массив Рута не дает никаких изменений знака.

Шаг 1: Напишите полином в стандартной форме

Для нашего примера, коэффициенты являются:
a₄ = 1 (^]4, a₃3, a₂ = 35², a_{1¹, a = 24 (постоянный термин).}

Шаг 2: Инициируйте первые два ряда

Первый ряд массива Рут содержит коэффициенты четноиндексированных терминов (начиная с sⁿ): строка s⁴:[1, 35, 24]. Второй ряд содержит коэффициенты нечетноиндексированных терминов: строка s³:[10, 50, 0]. (Заполнить недостающие термины нулями.)

Так:
Row s⁴: 1, 35, 24
Row s³: 10, 50, 0

Шаг 3: Вычислите последующие строки, используя правило Ромбуса

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

Но проще применить стандартную формулу: для строки k, столбца j (начиная с столбца 1) элемент:

Давайте сделаем это вручную для ясности.

Сокращение s²⁴ и s³,1
R²,1 = —[[50 — 350] — (10×35)][10 = —[50 — 350]/10 = —[10] = —[50 — 350]/10 = —[10] = —^{,24,3[24]; R3,3=0.
R}
,2
—
,2 = —[1×0] —[240]/10

Следовательно, строка s²: 30, 24, 0.

Строки s¹:³ и s^{2
R1,1 = —[[10×24] — (30×50)]/30 = —[240 — 1500]/30 = —[1260]/30 = 42.
R,2 = —[ (10×0) — (30×0)]/30 = 0.}

В этом случае s¹:42, 0.

Строки s^0
2 и s^{1
R0,1) = -[[30×0] - (42×0) - (0.42 = 0. Но на самом деле существует известное правило: последний элемент последней строки равен a₀ (константный термин).0,1) = -[[30×0] - (42×24)]/42 = -[1008]/42 = 24. Правильно!}

Полный набор Routh:

Шаг 4: Проверить первую колонну

Первая колонка содержит: 1, 10, 30, 42, 24. Все положительные, и никаких знаковых изменений нет. Поэтому система устойчива, что соответствует нашим известным корням.

Особые случаи

Два особых случая могут возникнуть при строительстве решетки Рута:

• Нуль в первом столбце: Если первый элемент строки равен нулю, но строка не является всеми нулями, замените ноль небольшим эпсилоном (ε), затем проследуйте. Знак ε определяет изменения знака. Это указывает на предельную стабильность или нестабильность в зависимости от остальной части массива.
• Целая строка нулей: Это происходит, когда имеются симметрично расположенные корни (например, пары мнимых корней или симметричные пары вокруг реальной оси).В этом случае образуют вспомогательный полином из строки выше нулевой строки и дифференцируют его для получения коэффициентов для отсутствующей строки.Это указывает на наличие пары чисто мнимых корней или пары реальных корней с противоположными знаками (то есть, неустойчивость).

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

Предсказание реакции системы и переходного поведения

Критерий Рут-Гурвица не ограничивается двоичным стабильным/нестабильным вердиктом. Выявляя число полюсов с положительными реальными частями, он напрямую сообщает о переходном поведении системы. Каждый полюс в правой полуплоскости вносит экспоненциально растущий режим, вызывая переходную реакцию на расхождение. И наоборот, все полюса в левой полуплоскости приводят к разложению экспоненциальных, в конечном итоге достигая устойчивого состояния.

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

• Если в первой колонке нет изменений знака:] Все полюса имеют отрицательные реальные части. Переходная реакция будет стабильной, в конечном итоге осядет в стабильное состояние. Но может быть медленный или колебательный распад (почти воображаемые полюса), который критерий не обнаруживает — гарантируется только знак реальной части.
• Если есть изменения знака: Количество изменений знака равно числу полюсов RHP. Система неустойчива, и переходная реакция будет расти неограниченно. Число полюсов RHP намекает на природу нестабильности: один положительный полюс приводит к монотонному экспоненциальному росту, а пара сложных полюсов RHP производит растущие колебания.
• Маргинальная стабильность: Если характеристическое уравнение имеет корни на воображаемой оси (чисто воображаемые или множественные полюса на оси), массив Рута будет производить ряд нулей. Это указывает на то, что система находится на границе стабильности, и переходный ответ может проявлять устойчивые колебания (если простые воображаемые полюса) или нарастающие колебания (если повторяющиеся воображаемые полюса).

Объединив результат Рут-Гурвица со знанием системного порядка, инженеры могут делать обоснованные прогнозы. Для системы второго порядка критерий сводится к требованию, чтобы все коэффициенты были положительными (что необходимо и достаточно для стабильности). Для систем более высокого порядка массив обеспечивает показатель запаса устойчивости — чем ближе элементы первой колонны к нулю, тем ближе система к нестабильности. Это особенно полезно при настройке параметров усиления в петле управления.

Взаимосвязь с переходными характеристиками

Переходная реакция стабильной системы определяется главным образом доминирующими полюсами — наиболее близкими к воображаемой оси. Критерий Рута-Гурвица не определяет непосредственно местонахождение этих полюсов, но может помочь в настройке. Например, если критерий показывает, что увеличение усиления контроллера в конечном итоге приведет к изменению знака в первом столбце, он отмечает запас усиления. Инженеры часто используют это для установки пропорционального усиления в контроллерах P, PI или PID без необходимости каждый раз имитировать всю систему.

Ключевые переходные метрики, на которые влияет стабильность, включают:

• Образование — Максимальное пиковое значение ответа относительно конечного значения. Пересечение обычно вызвано сложными сопряженными полюсами с низким коэффициентом демпфирования. В то время как критерий Рута — Гурвица не вычисляет демпфирование, система с рядом нулей (предельная стабильность) часто соответствует коэффициенту демпфирования вблизи нуля, подразумевающему высокий пересечение.
• Время расселения — Время, необходимое для ответа, чтобы оставаться в пределах определённого процента (например, 2%) от конечного значения. Время расселения обратно пропорционально реальной части доминирующего полюса. Массив Рута не может дать эту реальную часть, но если массив показывает, что все полюса находятся в левой полуплоскости, время расселения конечно; если любой полюс очень близок к воображаемой оси, время расселения будет длительным.
• Время подъема — Время перехода от 10 % до 90 % конечного значения. Это связано с естественной частотой доминирующих комплексных полюсов. Опять же, критерий Рута — Гурвица не даёт прямого численного значения, но указывает, достаточно ли стабильна система, чтобы иметь чётко определённое время подъема.
• Колебания — Частота переходных колебаний определяется мнимыми частями сложных полюсов.Критерий не выявляет мнимых частей, но наличие ряда нулей указывает на чисто мнимые полюса, которые производят устойчивые колебания в переходном ответе.

Таким образом, критерий Рут-Гурвиц действует как привратник: он обеспечивает стабильность как предпосылка для значимого переходного анализа с использованием других инструментов, таких как корневой локус, участки Боде или симуляции реакции на шаг.

Пример: система третьего заказа

Применим критерий Рут-Гурвица к более сложному примеру с неизвестными параметрами. Рассмотрим функцию переноса замкнутого цикла с обратной связью единства. Характеристическое уравнение:

Мы хотим определить диапазон K, для которого система стабильна. Постройте массив Рута для многочлена .

Rows:
s³: 1, 10
: 7, K
^:: На самом деле проще: элемент последней строки — это постоянный термин K⁰, полученный из строк s² и s^{1. Используя формулу: (он равен постоянному термину напрямую).}

Для устойчивости все элементы первой колонны должны быть положительными: 1 > 0, 7 > 0, (K-70)/7 > 0, K > 0. Это дает K > 70 и K > 0, так K > 70. Нижняя граница K > 70. Если K = 70, строка s¹ становится нулевой, что является особым случаем, указывающим на предельную стабильность (чисто воображаемые корни). Так система стабильна для K > 70. Это говорит инженеру, что при достижении K выше 70 дают стабильный переходный момент, но по мере приближения K 70 переходные колебания становятся незатухшими. Переходный ответ для K = 200, например, будет иметь два сложных сопряженных полюса в левой полуплоскости и один реальный полюс, что приводит к ответу перерасхода, но в конечном итоге оседает. Критерий Рута-Гурвица здесь явно предсказывает границу устойчивости.

Последствия переходного ответа для примера

Для K чуть выше 70 система стабильна, но доминирующие полюса очень близки к воображаемой оси. Переходная реакция будет показывать очень медленный распад и большой перерасход. Время оседания будет длительным. По мере увеличения K полюса будут двигаться влево (как правило), улучшая время демпфирования и оседания. Однако, если K слишком высок, массив Рута для этой конкретной системы не будет показывать никаких дополнительных изменений знака - стабильность сохраняется. На практике очень высокий K может вызвать другие проблемы, такие как насыщение привода или усиление шума, но критерий Рута-Гурвица только проверяет стабильность петли.

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

Ограничения критерия Рут-Гурвица

Несмотря на свою мощь, критерий Рут-Гурвица имеет важные ограничения, которые должен понимать каждый инженер управления:

• Применимо только к системам LTI: Он не может использоваться непосредственно для систем с изменяющимся временем, нелинейных или выборочных данных (хотя существуют расширения для дискретных систем времени посредством билинейных преобразований).
• Не дает количественной оценки переходной производительности: Критерий не дает никакой информации о превышении, пиковом времени, времени подъема или точном времени оседания. Для этого необходимы корневой локус, частотная реакция или моделирование временных доменов.
• Требует известного характеристического многочлена: Коэффициенты должны быть точно известны.На практике неопределенности в параметрах могут потребовать надежного анализа устойчивости (например, с использованием теоремы Харитонова).
• Стало громоздким для систем очень высокого порядка: Ручные вычисления для многочленов десятого порядка подвержены ошибкам, хотя программное обеспечение (MATLAB, Python) автоматизирует конструкцию массива.
• Не могут обрабатывать задержки времени напрямую: Системы с задержками времени имеют трансцендентальные характерные уравнения — массив Рута не может быть непосредственно применен, если не использовать приближения Паде.

Для преодоления этих ограничений инженеры обычно используют критерий Рут-Гурвица в сочетании с другими методами. Например, участки корневого локуса показывают, как полюса замкнутого цикла движутся с усилением, предоставляя точную переходную информацию. Критерий Найквиста обрабатывает временные задержки и нелинейности. Тем не менее массив Рут часто служит быстрой первой проверкой — особенно во время предварительного проектирования или устранения неполадок.

Дополнительные инструменты

Критерий Рут-Гурвица является «необходимым и достаточным» для устойчивости полиномов LTI, но для анализа производительности инженеры полагаются на:

• Корневой локус: Участок траектории полюсов как параметра (например, усиления) изменяется. Он непосредственно показывает отношение демпфирования и естественную частоту.
• Объекты кода и диаграммы найквиста: Предоставляют маржу усиления и фазовую маржу, которые являются показателями относительной стабильности и характеристиками переходного реагирования.
• Симуляция реакции на шаг: Дает фактические временные метрики (вылет, время урегулирования) для данной системы.
• Собственные значения матрицы состояний (являются полюсами) могут быть вычислены численно, а метод Ляпунова расширяет анализ устойчивости на нелинейные системы.