Table of Contents

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

Анализ устойчивости формирует основу проектирования системы управления, гарантируя, что система реагирует предсказуемо и остается ограниченной в нормальных условиях работы. Для дискретных систем времени & #8212; распространен в цифровом управлении, обработке сигналов и встроенных приложениях & #8212; критерий стабильности переходит от непрерывной плоскости s-плоскости левой половины плоскости к единичному кругу в плоскости Z. Инженеры полагаются на надежные аналитические инструменты для определения того, находятся ли все полюсы дискретной функции передачи внутри единичного круга без явного решения многочленов высокого порядка. Два выдающихся метода - критерий Рут-Гурвица, первоначально разработанный для непрерывных систем, но адаптируемый к дискретным структурам, и тест стабильности Jury & #8217, специально построенный для дискретного анализа. Понимание того, когда и как применять каждый метод, имеет решающее значение для проектирования надежных цифровых контроллеров и систем выборочных данных.

В этой статье приводится всестороннее сравнение теста стабильности Routh-Hurwitz и Jury’, исследуются их математические основы, процедурные шаги, практические адаптации и относительные сильные стороны. К концу инженеры и студенты получат четкую дорожную карту для выбора соответствующей техники для их потребностей в анализе дискретных систем.

Дискретная стабильность системы: критерий окружности единицы

Прежде чем погрузиться в конкретные методы, необходимо вспомнить основное условие устойчивости для дискретно-временных систем. Линейная дискретная система, описанная передающей функцией H(z) = N(z)/D(z), стабильна, если и только если все полюса—корни характерного многочлена D(z) = 0—внутри единичного круга в z-плоскости, т.е. |z | <1.Это условие заметно отличается от непрерывных систем, где стабильность требует, чтобы все полюса имели отрицательные реальные части (левая половина плоскости).

В то время как алгоритмы поиска корней могут непосредственно вычислять местоположения полюсов, они становятся громоздкими для полиномов высокого порядка. Следовательно, алгебраические критерии, такие как тест Routh-Hurwitz и Jury’, предлагают эффективные, неитеративные проверки. Ключевая проблема заключается в том, что тест Jury’ работает непосредственно на коэффициентах D(z), тогда как Routh-Hurwitz обычно работает на полиномах в домене s. Таким образом, преобразование необходимо для применения Routh-Hurwitz к дискретным системам.

Критерий Рут-Гурвица и его адаптация к дискретным системам

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

Критерий Рут-Гурвица обеспечивает необходимые и достаточные условия устойчивости непрерывных линейных систем.

a0nna1n–1an = 0

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

Билинейная трансформация: преобразование z-плана в s-план

Для применения Routh-Hurwitz к дискретным системам инженеры используют билинейное преобразование (также называемое преобразованием Tustin & #8217). Это отображение связывает переменную z с новой сложной переменной w (часто обозначается как s для удобства после преобразования):

z = (1 + wT/2) / (1 & #8211; wT/2)

или, для нормального периода выборки T = 2, более простая форма z = (1 + w / (1 & #8211; w. Под этой билинейной картой внутренняя часть единичного круга в z-плане плоскости в w-плане может быть проанализирована с использованием стандартного массива Рут-Гурвица.

Шаги для применения Routh-Hurwitz к дискретным системам

  1. Получить дискретный характеристический многочлен D(z) из функции передачи системы’s.
  2. Заменить билинейное преобразование в D(z) для получения полинома в w. Этот шаг может включать алгебраическое расширение и упрощение, особенно для систем более высокого порядка.
  3. Применить критерий Рута-Гурвица к преобразованному многочлену в w. Постройте массив Рута и изучите признаки записей первой колонны.
  4. Вывод о стабильности: Если все записи первой колонки положительные, дискретная система стабильна (все полюса внутри единичного круга). Любые изменения знака указывают на число полюсов вне единичного круга.

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

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

  • Алгебраическое бремя: Билинейное преобразование расширяет полином, потенциально вводя ошибки округления и делая ручные вычисления утомительными.
  • Численность: Для многочленов высокого порядка величины коэффициента могут широко варьироваться, что приводит к плохо кондиционированным массивам.
  • Специальные случаи: Метод требует обработки вспомогательных полиномов для корней на воображаемой оси (в w-плоскости), которые соответствуют полюсам на единичном круге – предельное условие стабильности, которое часто нежелательно.
  • Неэффективность для чисто дискретных систем: Поскольку метод не является родным для дискретного домена, многие учебники по управлению рекомендуют тест Jury’ как более прямую и менее подверженную ошибкам альтернативу.

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

Тест стабильности жюри: разработан для дискретного домена

Происхождение и принцип

Тест стабильности Jury’, представленный Элиау И. Жюри в 1960-х годах, является дискретным аналогом критерия Рут-Гурвица. Он непосредственно изучает коэффициенты характеристического многочлена D(z), не требуя преобразования домена. Вместо проверки знака записей первой колонки, тест Jury’ строит таблицу (таблица жюри) из многочленных коэффициентов и обеспечивает набор условий, которые гарантируют, что все корни имеют величину менее одной.

Рассмотрим многочлен порядка n:

D(z) = a0zn + a1zn–1 + … + an–1z + an = 0, with a0 > 0.

Тест Jury’s обеспечивает систематическую процедуру определения того, лежат ли все корни внутри единичного круга, образуя таблицу с n + 1 ряд.

Порядок построения таблицы жюри

  1. Перечислите коэффициенты в первом ряду (ряд 0) в естественном порядке: a 0, a 1, ..., an.
  2. Формируйте второй ряд (ряд 1), обратив порядок коэффициентов в строке 0:an, an–1, ..., a0.
  3. Для последующих строк k = 2, 3, ...], вычислить для каждой пары строк., i


  4. Продолжай до тех пор, пока в первой колонке не останется только один ненулевой элемент, или пока таблица не рухнет.
  5. Проверить условия: Для стабильности должно быть следующее:
    • D(1) > 0
    • (–1)nD[–1) >0
    • Для таблицы жюри первый элемент каждой нечетной строки (ряды 1, 3, 5, ...) должен быть больше нуля. Эквивалентные условия включают ведущие записи уменьшенных многочленов.

Эти условия вместе необходимы и достаточны для того, чтобы все корни удовлетворяли |z | <1.На практике тест реализуется в компьютерных системах алгебры, но ручной процесс управляется для заказов до примерно 4 или 5.

Пример: система второго порядка

Рассмотрим дискретную систему второго порядка с характерным многочленом D(z) = z2z2 = 0. Применение теста жюри’:

  • Степень 0: 1, a1, a2
  • Рука 1: a 2, a1, 1
  • Строка 2: Следующая строка содержит записи b0 = 1 & #8211;a2, baaa2

Условия стабильности для второго порядка:

  1. D(1) = 1 + a1 +a2 > 0
  2. [–1)2D[–1) = 1 –a1a2 >
  3. a2 | <1]

Это знакомые условия для стабильной дискретной системы второго порядка. Обратите внимание, что тест Jury’ напрямую приводит к этим простым неравенствам без каких-либо преобразований.

Преимущества теста Jury’

  • Прямое приложение: Работает над z-доменным многочленом; отображение не требуется.
  • Полный набор условий: Обеспечивает как необходимые, так и достаточные критерии устойчивости.
  • Средства проверки: Включает проверки корней на единичном круге (предельная стабильность) через D(1) и D(–1) условия.
  • Пригодно для цифровой реализации: Конструкция таблиц алгоритмична и легко программируется.

Сравнение результатов теста стабильности Роут-Хурвица и жюри

AspectRouth-Hurwitz (Adapted)Jury’s Stability Test
DomainContinuous (s-plane) originally; adapted via bilinear transform to w-planeDiscrete (z-plane) directly
Stability regionLeft-half plane (after transformation, maps to inside unit circle)Inside unit circle (|z|<1)
ProcedureConstruct Routh array from w-polynomial; check sign of first columnConstruct Jury table from z-polynomial; check series of inequalities
ComplexityModerate after transformation; transformation adds extra algebraic stepsLow to moderate; table construction straightforward for low to medium orders
Numerical stabilityCan suffer from coefficient scaling after bilinear transformGenerally good; determinants may cause swelling for high orders
Handling of marginal stabilityAuxiliary polynomial required; roots on imaginary axis in w-planeDirect checks through D(1) and D(-1) conditions
Educational accessibilityStudents often learn continuous version first; adaptation extends existing knowledgeRequires understanding of unit circle concept; but method is self-contained
Software implementationEasily coded via array construction; bilinear transform may require symbolic toolboxEasily coded via recursive table; native to z-domain
Typical use caseWhen system is already in s-domain or hybrid continuous/discreteWhen characteristic polynomial is directly from a discrete transfer function

Подробные точки сравнения

1.Применимость и пригодность домена

Тест Jury’ является естественным выбором для чистых дискретных систем, особенно в приложениях цифрового управления, где производится выборка вывода контроллера и входа установки. Routh-Hurwitz требует преобразования, которое добавляет сложность и может заслонять физическую интерпретацию полюсов относительно единичного круга. Многие современные учебники и курсы по цифровому управлению представляют тест Jury’ в качестве стандартного инструмента.

2 Алгебраическая сложность

Для заказов до 4 или 5 оба метода управляемы вручную. Однако тест Jury’ избегает билинейного расширения, которое может производить термины со смешанными полномочиями и большими коэффициентами. Например, третий порядок D(z), преобразованный в w, сам по себе дает полином третьего порядка, но коэффициенты являются комбинациями исходных; для более высоких заказов преобразование становится более трудоемким. Конструкция таблицы Jury’, в то время как включает детерминанты, работает непосредственно на исходных коэффициентах.

3 Численность хрупкости

Рут-Хурвиц в w-домене может страдать от плохого кондиционирования, если билинейное преобразование создает коэффициенты с очень разными величинами.Тест жюри’ также сталкивается с численными проблемами для многочленов высокого порядка (например, порядок >10) из-за повторных вычитаний и детерминантных вычислений, но они менее серьезны, чем те, которые вводятся преобразованием. Оба метода выигрывают от реализации с двойной точностью в программном обеспечении.

4 Способность определять число нестабильных поляков

Рут-Хурвиц, в силу изменения знака Рута’, может непосредственно указывать число полюсов в правой половине w-плана, что соответствует числу полюсов вне единичного круга.Тест жюри’, в своей основной форме, обеспечивает ответ «да/нет стабильности».Однако расширения теста жюри’ (например, алгоритм жюри-Гантмахера) также могут подсчитывать количество корней вне единичного круга, хотя это менее часто преподается.

5. Обработка специальных полиномиальных коэффициентов

Когда дискретная система имеет нулевые коэффициенты (отсутствующие мощности), оба метода требуют осторожности. Рут-Хурвиц имеет дело с нулевыми входами в массив с помощью небольших приближений эпсилона. Тест жюри & #8217 аналогично должен управлять нулями во время построения таблицы, но алгоритм остается действительным, пока ведущий коэффициент a 0 0 не равен нулю. Если строка в таблице жюри становится полностью нулевой, полином имеет симметричные корни, указывающие на предельную стабильность или антисимметрию.

Практическое применение и когда выбрать каждый метод

Тест жюри когда

  • Вы работаете исключительно в дискретной области (z-план).
  • Вы хотите одноэтапный тест без преобразования домена.
  • Вы анализируете цифровой контроллер или систему выборочных данных с известным периодом отбора проб.
  • Нужны четкие условия для предельной устойчивости (корни на единичном круге).
  • Порядок многочлена является умеренным (≤ 8 для ручных проверок; выше для программного обеспечения).

Используйте Routh-Hurwitz (через Билайнарную Трансформацию)

  • Дизайн вашей системы управления уже находится в непрерывной области, и вы будете дискретизировать позже.
  • Вы предпочитаете унифицированную структуру анализа для непрерывных и дискретных систем.
  • У вас есть программное обеспечение (например, MATLAB, Python с SymPy) для автоматического управления билинейным преобразованием.
  • Необходимо подсчитать точное количество нестабильных полюсов (посредством знаковых изменений в массиве Рута).
  • Вы учите связи между непрерывными и дискретными критериями стабильности.

Примеры реального мира

Пример: Цифровой контроллер скорости двигателя.] Микроконтроллер реализует контроллер PI для двигателя постоянного тока. В Z-домене выведен полином с закрытой контурной характеристикой. Инженер будет использовать тест Jury’ для проверки стабильности различных скоростей отбора проб и пропорционального усиления. Тест быстро показывает, перемещаются ли полюса за пределы единичного круга по мере увеличения усиления.

Пример: Адаптивная конструкция фильтра.] При обработке сигналов IIR-фильтры требуют полюсов внутри единичного круга, чтобы избежать самоколебания. Тест Jury’s обычно встраивается в библиотеки дизайна фильтра для проверки стабильности после квантования коэффициента.

Пример: Гибридная непрерывно-дискретная система.] Непрерывная установка, управляемая цифровым компенсатором, требует дискретизации через удержание нулевого порядка.Общее характеристическое уравнение в z-домене может быть проанализировано либо непосредственно тестом Jury’s, либо с помощью билинейного преобразования для применения Routh-Hurwitz.Многие инженеры по управлению предпочитают тест Jury’s, поскольку они могут оставаться в z-домене и избегать дополнительных преобразований.

Реализация и автоматизация программного обеспечения

Оба метода доступны в наборах инструментов системы управления.

  • MATLAB: Функция (или пользовательские скрипты) может создавать массив Routh, в то время как доступна в Robotics Toolbox или может быть запрограммирована. Функция выполняет дискретизацию и билинейное преобразование, если это необходимо.
  • Python (SciPy и управляющие библиотеки): Библиотека включает для теста Jury’s. Для Routh-Hurwitz используется функция или конструкция ручного массива.
  • Символические инструменты: SymPy и Mathematica позволяют осуществлять алгебраические манипуляции для билинейного преобразования и последующего генерирования массивов.

При использовании программного обеспечения этап трансформации для Routh-Hurwitz автоматизирован, что уменьшает ручную ошибку.Однако для быстрого прототипирования дискретных систем тест Jury’ остается более прямым выбором.

Ограничения и пещеры

Ни один из методов не является панацеей. Оба становятся громоздкими для полиномов высокого порядка (order > 10) из-за накопления числовых ошибок. Для таких случаев вычислительное определение корня (например, собственные методы) часто более надежно. Кроме того, оба теста предполагают, что характерный полином точно известен; моделирование неопределенностей и вариаций параметров требует надежного анализа стабильности за пределами этих алгебраических критериев.

Еще одна важная оговорка: адаптация Рута-Гурвица посредством билинейного преобразования действительна только для линейных временных инвариантных систем. Она напрямую не распространяется на изменяющиеся во времени или нелинейные дискретные системы. Тест Jury’s аналогичным образом применяется только к линейным дискретным моделям.

Наконец, может произойти неправильное толкование условий. Для теста жюри & #8217, условие D(1) > 0 и (–1)nD(–1) >0 необходимы, но не достаточны в одиночку; они должны быть объединены с условиями таблицы. Студенты иногда останавливаются после этих двух простых проверок, что может привести к неправильным выводам о стабильности для систем более высокого порядка.

Заключение

Выбор между Routh-Hurwitz и Jury’s тестом стабильности для дискретных систем зависит от контекста, знакомства и вычислительных ресурсов. Jury’s тест является более прямым, специально построенным методом для дискретных полиномов времени, предлагая четкие условия и избегая билинейной трансформации. Адаптированный критерий Routh-Hurwitz обеспечивает мост для инженеров, обученных непрерывным системам и может быть полезен при интеграции непрерывного и дискретного анализа. Оба метода остаются ценными учебными инструментами, которые укрепляют фундаментальную концепцию расположения полюсов относительно единичного круга.

Для большинства современных цифровых систем управления и дискретного системного анализа рекомендуется тест стабильности Jury’. Однако знание обеих методов гарантирует, что инженер может адаптироваться к различным заявлениям о проблемах и использовать доступное программное обеспечение для проверки результатов. Понимание сильных и слабых сторон каждого метода позволяет дизайнерам принимать обоснованные решения, что в конечном итоге приводит к более надежным и надежным системам управления.

Дальнейшее чтение и внешние ресурсы