Упрощение сложных полиномов с помощью анализа стабильности Рут-Хурвица

Почему прямые расчеты корней непрактичны

Для полиномов порядка пяти и выше точные аналитические корневые решения редко существуют (теорема Абеля — Руффини). Численные методы могут аппроксимировать корни, но они чувствительны к неопределенности коэффициентов и могут пропускать тонкие признаки нестабильности. Критерий Рута — Гурвица, напротив, дает окончательный ответ «да/нет» о знаке реальных частей, используя только полиномиальные коэффициенты. Это делает его незаменимым для предварительного анализа стабильности и проверки конструктивных границ.

Критерий стабильности Рута-Хурвица

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

Построение решетки Routh Array

Учитывая полином:

P(s) = a_nnⁿn−1n−1 + ... + a₁₀a_n > 0

Рутский массив построен следующим образом:

1. Обозначить первые два ряда. Обозначить ⁿ содержит коэффициенты равноиндексированных степеней _n, n−2, ...] Обозначитьn−1^{n−1, , a_n−3, ...] Любые недостающие коэффициенты заполнены нулями.}
2. Вычислите последующие строки.ⁿ⁻², каждый элемент получен из двух строк выше с использованием определяющей формулы:
   b₁n−1n−2_n−3n−1,
   bn−1_n−4 —n−5_n−5n−1 и так далее, пока остальные элементы не станут нулевыми.
3. Продолжай до тех пор, пока не останется только постоянная строка. Последняя строка соответствует s⁰ и содержит один элемент, который является постоянным термином (обычно a₀).
4. Исследуйте первую колонку. Если все записи в первой колонке массива положительные (и ни одна не равна нулю), то все корни лежат в левой полуплоскости (система стабильна).Любое изменение знака указывает на количество корней с положительными реальными частями — каждое изменение знака считается одним таким корнем.

Особые случаи в Routh Array

Возникают две особые ситуации, требующие тщательного обращения:

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

Пример: Полиномиал четвертого порядка

Рассмотрим P(s) = s⁴ + 2s³ + 3s² + 4s + 5.

• ⁴: 1, 3, 5
• ³: 2, 4, 0
• Row s²:b₁ = (2·3 — 1·4)/2 = (6-4)/2 = 1, b₂ = (2·5 — 1·0)/2 = 10/2 = 5
• Row s¹:c₁ = (1·4 — 2·5)/1 = (4-10)/1 = —6, затем ноль
• Row s⁰:d₁ = (-6·5 — 1·0)/(-6) = (-30)/(-6) = 5

Первый столбец: 1, 2, 1, -6, 5 → знак меняется один раз (с 1 до -6, а с -6 назад до 5), указывая один корень с положительной реальной частью. Действительно, у полинома есть корни: -1,65 ± 0,35i, 0,65 ± 1,72i — последняя пара имеет положительную реальную часть.

Упрощение сложных полиномов с помощью критерия

Метод Рута-Гурвица упрощает анализ устойчивости сложных полиномов несколькими способами:

1. Избегает вычисления корней. Для многочленов высокого порядка факторинг или численный поиск корней является дорогостоящим и иногда ненадежным. Массив требует только сложения, вычитания, умножения и деления — операции, легко выполняемые вручную или электронной таблицей.
2. Неопределенность параметров сцепления. Выражая коэффициенты в терминах параметров системы (прибыли, константы времени), можно вывести условия устойчивости в терминах неравенств. Например, диапазон усиления контроллера K, для которого система замкнутого цикла остается стабильной, можно найти, применяя массив Рута и требуя, чтобы все элементы первой колонны были положительными.
3. Идентифицирует предельную стабильность. Ноль в первом столбце (случай ε) или ряд нулей указывает корни на воображаемой оси или симметрично расположенной вокруг происхождения. Эти условия часто определяют границу между стабильным и нестабильным поведением.

Работа с комплексными коэффициентами?

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

Применение в системах управления и за их пределами

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

Проектирование систем управления

Во время проектирования контроллера инженеры часто корректируют коэффициент усиления (например, пропорциональный, интегральный и производный коэффициент усиления контроллера PID) для достижения желаемой производительности. Критерий Routh обеспечивает быстрый способ найти ]прибыль — диапазон коэффициентов усиления, которые поддерживают стабильность системы. Например, при наличии системы обратной связи с единством с передающей функцией G(s) = K/[(s+1)(s+2)(s+3)], характерным многочленом является ³2 + 11s + (6+K) = 0. Применение массива Routh дает условие K < 60 для стабильности. Это намного быстрее, чем попытка вычислить корни как функцию K.

Сетевой анализ

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

Контроль промышленных процессов

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

Преимущества и ограничения

Как и любой инструмент, критерий Рут-Гурвица имеет сильные и слабые стороны.

Сильные стороны

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

Ограничения

• Реальные коэффициенты. Классическая форма не обрабатывает напрямую полиномы со сложными коэффициентами, хотя существуют расширения.
• Никакой информации о местоположении корня. Критерий говорит только , сколько нестабильных корней существует, а не их точные значения. Для детального транзиторного анализа необходим локус корня или численный корневой поиск.
• Чувствительность к численным ошибкам.] Когда коэффициенты почти равны нулю или содержат ошибки округления, конструкция массива может создавать двусмысленности знаков. Использование ε-метода является обходным, но требует тщательной интерпретации.
• Не обрабатывает системы задержек. Системы с задержками времени имеют трансцендентальные характерные уравнения, а не полиномы.Критерий Рута — Гурвица неприменим непосредственно; вместо этого используются методы, подобные приближениям Паде или предиктору Смита.

Исторический контекст и дальнейшее чтение

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

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

Шаг за шагом иллюстрированный пример

Давайте рассмотрим полный пример со сложным (высокопорядковым) многочленом, чтобы продемонстрировать силу метода. Рассмотрим характерный многочлен системы управления:

P(s) = 2 s⁵ + 4 s⁴ + 6 s³ + 8 s² + 10 s + 12

Шаг 1: Устройте первые два ряда

• ⁵[мощность 5]]⁵,,¹ → 2, 6, 10
• s⁴[мощность 4]]⁴,²,⁰ → 4, 8, 12

Шаг 2: Вычислите третий ряд s ³

Используя определяющую формулу:

• Первый элемент: (4·6 - 2·8)/4 = (24 - 16)/4 = 8/4 = 2
• Второй элемент: (4·10 - 2·12)/4 = (40 - 24)/4 = 16/4 = 4
• Третий элемент: ноль (поскольку в первом ряду коэффициентов больше нет)

³: 2, 4, 0

Шаг 3: Вычислить четвертый ряд (s²)

⁴ (4,8,12]) и³ (2,4,0):

• Первый элемент: (2·8 - 4·4)/2 = (16 - 16)/2 = 0/2 = 0
• Второй элемент: (2·12 - 4·0)/2 = (24 - 0)/2 = 12

²: 0, 12

Здесь у нас есть ноль в первой колонке. Остальная часть строки ненулевая (12). Замените ноль на ε (небольшое положительное число).

Шаг 4: Вычислите пятый ряд s ¹

³ (2,4,0]] и²[E,12]]:

• Первый элемент: (ε·4 - 2·12)/ε = (4ε - 24)/ε = 4 - 24/ε
• Второй элемент: (ε·0 - 2·0)/ε = 0

s¹: (4 — 24/ε), 0

Как ε → 0⁺, доминирует термин —24/ε, что делает этот элемент отрицательным (−∞).

Шаг 5: Вычислите шестой ряд (s⁰)

Используя строку ² (ε,12]) и строку ¹ (4 — 24/ε, 0):

• Первый элемент: [(4 - 24/ε)·12 - ε·0] / (4 - 24/ε) = 12

⁰:12

Шаг 6: Проанализируйте первую колонку

Первые столбцы: строка s⁵:2 (положительный), строка s⁴:4³:2²: ε (положительный, как ε>0), строка s¹: 4 — 24/ε (отрицательный для малого ε), строка s⁰:12 (положительный).

Есть два изменения знака: от положительного (ряд: 0) 2 [ FLT: 1]] до отрицательного (ряд: 2]] 1 [ FLT: 3]], и от отрицательного обратно к положительному (ряд: 4] 0 [ FLT: 5]]. Поэтому полином имеет два корня с положительными реальными частями [ FLT: 7].

Проверка: Используя численные инструменты, корни 2 s⁵ + ... + 12 = 0 приблизительно: -1.09, -0.50 ± 0.87i, 0.79 ± 1.04i. Действительно, два сложных конъюгированных корня имеют положительную реальную часть (0.79), подтверждающую предсказание Рута-Гурвица.

Альтернативный подход: вспомогательный полином для случая нулевого ряда

Если бы мы столкнулись с полным рядом нулей вместо просто нуля в первом столбце, мы бы сформировали вспомогательный полином из строки выше. Этот метод необходим для таких случаев, как s ⁴ + 5s ² + 4 (который имеет чисто воображаемые корни ±j, ±2j].