Как расширить анализ Рут-Хурвица до систем фракционного порядка

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

Классический критерий Рут-Гурвица

Перед исследованием расширений необходимо пересказать классический метод. Для системы с характерным многочленом

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

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

Системы фракционного порядка: краткий обзор

Система дробного порядка — это система, динамика которой описывается дробным дифференциальным уравнением формы.

где порядки дифференциации αi и βi не обязательно являются целыми числами. Преобразование Лапласа такого уравнения (при нулевых начальных условиях) даёт передаточную функцию, содержащую такие термины, как sγ с γ дробным.sα с 0 <]. Эти системы проявляют свойства, не наблюдаемые в системах целого порядка: бесконечномерное представление состояния-пространства, неэкспоненциальное расслабление и память на большие расстояния.

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

Почему классический тест Рут-Гурвица не подходит для фракционных систем

Препятствие фундаментальное: уравнение характеристик дробного порядка, такое как

Силы s не являются целыми кратными друг другу.Массив Рута требует целых сил, чтобы можно было осуществить систематическое устранение.Кроме того, дробные силы приводят к разрезам ветвей в сложной плоскости; понятие «полуполосный полюс правой плоскости» становится более тонким, поскольку дробная сила может иметь множество листов Римана. Дробная система может иметь бесконечное число полюсов (если порядок иррациональный) или конечное число, если приказы соразмерны.

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

Критерий стабильности Матиньона: алгебраическая замена

Наиболее известным расширением является критерий стабильности Матиньона (также называемый теоремой Матиньона). Он обеспечивает необходимое и достаточное условие для большого класса систем дробного порядка с соизмеримыми порядками. Говорят, что дробная система имеет соизмеримые порядки, если все дробные силы в характеристическом уравнении являются целыми кратными базовому порядку α (например, α, 2α, 3α, ...). Многие практические системы, особенно те, которые получены из дробных контроллеров PID или из дробных вязкоупругих моделей, соизмеримы.

Критерий Матиньона: Для системы дробного порядка с характеристическим уравнением

где α — положительное реальное число, система стабильна тогда и только тогда, когда

|arg i | > απ/2 для каждого корня iP.

Другими словами, корни должны лежать в области комплексной плоскости вне сектора угла απ/2, центрированного на отрицательной реальной оси.Это условие сводится к классическому состоянию левой половины плоскости, когда α = 1 (с тех пор απ/2 = π/2), но для α < 1 область стабильности является клином.

Критерий Матиньона является мощным, поскольку он превращает проблему стабильности обратно в проблему корневого местоположения для преобразованного полинома P в переменной w = sα. Можно применить тест Рута-Гурвица к полиному P(w) после замены w = s, но условие на аргументе s не захвачено только классическим массивом Рута. Однако для соизмеримых систем условие стабильности может быть проверено сначала построением массива Рута для полинома P(w) и затем проверкой состояния аргумента.

Это приводит к практической процедуре:

  1. Определить базовый дробный порядок α так, чтобы характеристическое уравнение стало многочленом в sα.
  2. Определить w = sα и переписать характеристическое уравнение как P(w) = 0, обычный полином.
  3. Применить классический тест Рута-Гурвица к P(w).Если в первом столбце происходят какие-либо изменения знака, система нестабильна в соответствии с классическим многочленным определением в w.
  4. Проверьте условие аргумента: Для каждого корня wi, вычислите соответствующие sii1/αiwi, как правило, сложные.

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

Альтернативные аналитические подходы

Функция Миттага-Леффлера и представление пространства-состояния

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

Система стабильна тогда и только тогда, когда

|arg(λ(A)) | > απ/2

где λ(A) являются собственными значениями матрицы A. Это прямой многовариантный аналог критерия Матиньона. Для соизмеримых систем это условие сводится к одной и той же клиновидной области устойчивости. Инженеры могут вычислять собственные значения и проверять их аргументы с помощью стандартной численной линейной алгебры.

Преобразование в систему целых чисел через переменную замену

Другой аналитический путь заключается в введении новой переменной z = sα, а затем применении классического теста Рута — Гурвица к преобразованному полиному.Однако этот подход проверяет только «полиномиальную часть» устойчивости.Необходимым условием является то, что полином в z имеет все корни в левой полуплоскости, но достаточное условие требует дополнительного теста аргумента.На практике инженеры часто выполняют оба: сначала запускают массив Рута на P(z), чтобы обнаружить очевидную полиномиальную нестабильность; затем, для остальных стабильных кандидатов, применяют условие аргумента.

Методы частотных доменов

Стабильность также можно оценить по критериям Найквиста или Боде, адаптированным к функциям переноса дробного порядка. Функция переноса открытого цикла G(s) системы дробного порядка является нерациональной функцией, но её частотный отклик можно оценить численно. Критерий устойчивости Найквиста по-прежнему сохраняется: число окружностей точки −1 определяет стабильность замкнутого цикла. Однако часто предпочтительнее прямое применение метода Рут-Гурвица, поскольку он избегает построения графика и обеспечивает алгебраическое состояние.

Численные и вычислительные инструменты

Учитывая сложность систем дробного порядка, инженерная практика в значительной степени зависит от программных инструментов. Инструментарий FOMCON (Моделирование и управление фракциональными порядками) для MATLAB является популярным выбором. Он предоставляет функции для вычисления проверок стабильности с использованием критерия Матиньона, для построения приближений (например, рекурсивное приближение Oustaloup) дробных операторов и для выполнения моделирования временных доменов.

Другие инструменты включают:

  • Средство управления CRONE (Commande Robuste d’Ordre Non Entier): разработано в Университете Бордо, поддерживает конструкцию контроллера дробного порядка и анализ стабильности.
  • Питоновые библиотеки: фракционные (для дробного исчисления) и управление (для общего анализа системы управления) могут использоваться вместе для реализации критерия Матиньона.
  • Символические вычисления в Mathematica или Maple могут обрабатывать алгебраические манипуляции дробных сил и поиск корней в комплексной плоскости.

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

  1. Определите характеристическое уравнение как символическое выражение в s с дробными силами.
  2. Преобразовать в многочлен в w = sα, если система соизмерима.
  3. Вычислите корни многочлена в w с помощью численного решателя.
  4. Для каждого корня вычислите соответствующий и его аргумент.
  5. Проверить, если |arg(s) | > απ/2 для всех корней.

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

Практическая пошаговая процедура для инженеров

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

  1. Получите характеристическое уравнение из функции передачи системы’s. Для системы дробного порядка это будет включать такие термины, как sγ с γ нецелочисленным.
  2. Проверьте соизмеримость.] Определите, являются ли все дробные порядки целыми кратными общему основанию α. Например, порядки 0,5, 1,5, 2.0 кратными 0,5; порядки 0.7 и 1.2 не являются кратными одному основанию (если только основание не является 0,1, но это редко рационально). Если система соизмерима, переходите к полиномиальному преобразованию. Если нет, переходите к шагу 6 (метод частотной области).
  3. Определить многочлен в w = sα. Записать характеристическое уравнение как P(w) = 0. Например, s2.51.2 + 2 несоизмерим (2,5 и 1.2 не являются кратными общему основанию), так что этот шаг не удался бы.1.5+3 s0.5 с α=0.5 даётw3 + 3w + 2 = 0
  4. Применить классический тест Рута-Гурвица к P(w). Построить массив Рута. Если в первом столбце происходит какое-либо изменение знака, система нестабильна (под полиномиальной интерпретацией). Если все записи первой колонны имеют один и тот же знак, переходите к следующему шагу.
  5. Проверить условие аргумента. Вычислить корни P(w). Для каждого корня wi, вычислить sii. Сложная мощность многозначна; использовать основное значение. |arg(si.. Если это удерживает для всех корней, дробная система стабильна. Если какой-либо корень нарушает условие, система нестабильна.
  6. Для несоизмеримых систем используют частотно-доменный подход. Оцените характеристическое уравнение вдоль границы области устойчивости, которая для дробных систем является сектором. Найквистовский участок функции переноса с открытым контуром может указывать на стабильность с замкнутым контуром. Альтернативно, используйте численный корневой поиск в плоскости -плане путем дискретизации дробного оператора, например, используя приближение Oustaloup. Затем примените классический тест Рута-Гурвица к аппроксимированной системе целого порядка. Этот метод дает приблизительный результат и должен быть проверен с помощью моделирования.

Многие инженеры считают, что соизмеримый случай охватывает большую часть практических конструкций с дробным управлением, особенно дробные контроллеры PID (где производные и интегральные ордера кратны общему α, часто 0,5 или 0,1).

Иллюстративное исследование: Fractional-Order PID Controller

Рассмотрим PID-контроллер фракционного порядка (PIαDβ), управляющий установкой первого порядка.

Эта система несоизмерима, поскольку дробные порядки 2.5, 1.2, 0.1 не являются целыми кратными общей базе. Мы не можем непосредственно применять критерий Матиньона. Вместо этого мы приблизим дробные операторы, использующие фильтр порядка Усталоупа, скажем, 5 в диапазоне частот [0,01, 100] рад/с. Это дает функцию переноса целого порядка высокого порядка. Приближенный характерный многочлен будет иметь порядок 5×5 = 25 для каждого дробного термина, но мы можем уменьшить порядок, используя сбалансированное усечение. Затем классический тест Рута-Гурвица может быть применен к уменьшенному многочлену. Проверка стабильности с помощью моделирования подтверждает результат. Этот примерный метод является стандартом де-факто в промышленности, когда аналитические инструменты недоступны.

В отличие от этого, если контроллер спроектирован с соизмеримыми порядками (например, α = 0,5, поэтому заказы составляют 0,5, 1,0, 1,5, 2.0), система становится соизмеримой. Тогда мы можем использовать критерий Матиньона напрямую, получая точное необходимое и достаточное условие (модуло интерпретации разветвления).

Ограничения и текущие исследования

Расширение анализа Рута-Гурвица на системы дробного порядка — это не одна формула, а набор методов, каждый со своей областью применимости. Критерий Матиньона точен только для соизмеримых систем. Для несоизмеримых систем приближенный подход целого порядка может вводить ошибки, особенно если система имеет слегка затухшие высокочастотные режимы. Исследователи изучают прямые алгебраические критерии для несоизмеримых систем, используя концепцию преобразования переменных дробного порядка s = (z)1/α и применяя модифицированный массив. Другим перспективным направлением является использование линейных матричных неравенств (LMI), адаптированных к системам дробного порядка.

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

Заключение

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

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