Понять сюжет Nyquist

Сюжет Nyquist является фундаментальным инструментом в области проектирования управления, используемым для анализа стабильности систем обратной связи. Названный в честь Гарри Найквиста, этот графический репрезентатив изображает частотную реакцию системы, нанося на график сложную передаточную функцию G(jω), поскольку угловая частота ω варьируется от 0 до ∞. Сюжет рисуется в сложной плоскости, с реальной частью G(jω) на горизонтальной оси и воображаемой частью на вертикальной оси. Его основная полезность заключается в оценке стабильности замкнутого цикла от частотного ответа открытого цикла, без явной вычисления полюсов замкнутого цикла. Сюжет Nyquist также обеспечивает понимание усиления и фазовых полей, что делает его незаменимым для надежного проектирования контроллера.

Математические основы сюжетов Nyquist требуют глубокого понимания сложного анализа, в частности оценки функций переноса на воображаемой оси и геометрической интерпретации контурного картирования. Инженеры должны понять концепцию контура Nyquist, замкнутого пути в комплексной плоскости, которая охватывает всю правую полуплоскость (RHP), и как этот контур отображает через функцию переноса открытого цикла G(s)H(s), чтобы произвести Nyquist Path. Полученный сюжет обвивает критическую точку −1 + 0j напрямую указывает на стабильность системы, следуя критерию стабильности Nyquist. Эта статья расширяет эти математические основы, проходя через логику построения, принцип аргументации и практические ловушки.

Комплексная оценка функции переноса

Функция переноса с открытым контуром G(s)H(s) является рациональной функцией комплексной переменной s = σ + jω. Для анализа частотных реакций мы ставим σ = 0, поэтому sjω. Получающееся выражение G(jω)H(jω) является комплексным числом, зависящим от ω:

G(jω)H(jω) = Re(ω) + j Im(ω) = M(ω) e^{jφ(ω)}

Здесь M(ω) — величина (также называемая усилением) и φ(ω) — фазовый угол. Построение графика Nyquist включает оценку G(jω)H(jω) на многих частотах и построение сложных точек. Для типичной системы график начинается с ω = 0 (часто в точке реальной оси) и заканчивается на ω → ∞ (обычно в исходной точке, если только система не имеет чистых интеграторов). Форма кривой раскрывает ключевые динамические свойства: например, как фаза уменьшается с увеличением частоты, или как величина откатывает.

Для вычисления G(jω)H(jω) математически рассмотрим общую функцию переноса, факторизованную в полюса и нули:

G(s)H(s) = K \frac{\prod {i=1}^{m} (s - z i)}{\prod {k=1}^{n} (s - p k)}

Когда s =jω, каждый термин становится сложным вектором от нулевого или полюсного местоположения до точки jω на воображаемой оси. Общая величина является произведением длин векторов полюсов, деленных на произведение длин векторов полюсов, умноженных на |K |. Общая фаза представляет собой сумму углов нулевых векторов минус сумму углов векторов полюсов, плюс фаза K (которая равна 0, если K является реальной и положительной). Эта геометрическая интерпретация является центральной для понимания того, как ведет себя кривая Найквиста, особенно вблизи частот, соответствующих естественным частотам системы.

Математические строительные шаги

Построение графика Nyquist вручную (или с использованием программного обеспечения) включает в себя систематическую оценку G(jω)H(jω) в диапазоне частот.

  1. Идентифицируйте передаточную функцию и запишите её в форме Боде (поли и нолики). Для систем с задержками во времени включите e^{-jωT} фактор.
  2. Вычислить величину и фазу на наборе частот.Как правило, выбираются низкие частоты (около 0) и высокие частоты (около ∞), с промежуточными точками вблизи частот разрыва (где знак изменения вкладов полюса или нуля).
  3. Начало при ω = 0: Точка зависит от типа системы. Для типа 0 (нет интеграторов) участок начинается при K (прирост постоянного тока) на положительной реальной оси. Для типа 1 (один интегратор) участок начинается с бесконечной величины с фазой −90°. Для типа 2 (два интегратора) фаза составляет −180° при ω = 0.
  4. Плот-точки на возрастающих частотах, соединяющие их плавно.По мере увеличения частоты кривая обычно движется к началу (поскольку уменьшается величина). Фаза может пересекать −180° на частоте пересечения фаз, где величина имеет решающее значение для анализа стабильности.
  5. Зеркало для отрицательных частот (если используется полный контур Найквиста).Плоскость для ω от −∞ до 0 представляет собой сложный конъюгат участка для положительных ω, отраженный по реальной оси.В сочетании с положительной частотной ветвью (и бесконечным полукругом, если система имеет чистый полюс-ноль избытка) получается полный путь Найквиста.
  6. Определить количество окружностей −1 + 0j для применения критерия устойчивости.

Математически кривая Найквиста для положительных частот может быть параметризирована как ω ∈ [0,∞] Для систем с относительной степенью (n − m) ≥ 1, кривая заканчивается у начала, приближаясь к нему фазой −(n−m)×90°. Для относительной степени 0 (одно и то же число полюсов и нулей) конечной точкой является конечное ненулевое комплексное число при ω → ∞.

Критерий стабильности никиста в деталях

Критерий стабильности Nyquist представляет собой мощный метод определения стабильности замкнутого цикла от частотного ответа открытого цикла. Он использует Cauchy's, также известный как принцип аргумента, от сложного анализа. Этот принцип гласит, что если замкнутый контур Γ в -плане -плане , отображенный через функцию F в F-плане, то число раз ΓF окружает начало равным P, где ZF(]P, внутри Γ. Для систем управления мы выбираем функцию F(s) = 1 + G(s)H(s)

Поскольку F(s) связана с функцией переноса в открытом цикле F(s) = 1 + G(s)H(s), карта F(s) через контур Nyquist просто Nyquist-график G(s)H(s) сдвинута вправо на 1. Эквивалентно, окружности начала в F-плане точно соответствуют окружностям точки −1 + 0j в G(s)H(s)-плане.

N = Z − P

где:

  • N = число окружностей −1 + 0j по графику Найквиста (считалось по часовой стрелке как положительное, против часовой стрелки как отрицательное).
  • Z = число полюсов замкнутого цикла в RHP (зеро 1 + G(s)H(s) в RHP.
  • P = число полюсов открытого цикла в RHP.

Для стабильной системы замкнутого цикла нам требуется Z = 0. Таким образом, система замкнутого цикла стабильна, если и только если участок Nyquist G(s)H(s) окружает точку −1 + 0j точно −P раз (т.е., окружности против часовой стрелки равны P, если считать по часовой стрелке положительными.P = 0 (стабильный открытый цикл), то участок Nyquist должен вообще не окружать −1 для стабильности замкнутого цикла.

Использование критерия: подсчет окружения

Подсчет окружений строго включает анализ пути Никиста по полному контуру, включая отображение бесконечного полукруга. Для систем с относительной степенью ≥ 1, бесконечные карты полукруга к маленькому кругу вокруг начала в плоскости G(s)H(s) (или к фиксированной точке, если относительная степень равна нулю). Практический подход для инженеров заключается в следующем:

  1. Нарисуйте Nyquist-график для ω от 0 до ∞.
  2. Отражайте его по реальной оси для ω от 0 до −∞.
  3. Если система имеет полюса на воображаемой оси (например, интеграторы), отступите контуром Nyquist с небольшими полукругами вправо, чтобы избежать этих полюсов, и нанесите эти углубления на большие дуги в плоскости G(s)H(s).
  4. Подсчитайте чистое число окружностей по часовой стрелке −1, проследив полную замкнутую кривую (включая бесконечное отображение полукруга).

Математически, число окружения N может быть вычислено через интеграл числа обмотки:

N = \frac{1}{2πj} \oint {Nyquist \, contour} \frac{G'(s)H'(s)}{1+G(s)H(s)} \, ds

Но это редко оценивается аналитически; инженеры полагаются на геометрический осмотр или программные инструменты. Ключевое понимание заключается в том, что окружения возникают, когда фаза G(jω)H(jω) пересекает −180° (т.е. точка −1 на реальной оси), в то время как величина ≥ 1. Это условие определяет маржа прибыли и фазовый запас .

Практические примеры строительства

Для укрепления понимания мы рассматриваем простую систему первого порядка и систему второго порядка с демпфированием.

Система первого заказа: G(s) = 1/(s+1)

Здесь G(jω) = 1/(jω+1). Магнитуда M(ω) = 1/√(1+ω2) и фаза φ(ω) = −tan−1(ω). Как ω = 0, точка приближается (0, 0) от угла −90°. Кривая Nyquist — полукруг в нижней полуплоскости (так как фаза всегда отрицательна) от положительной реальной оси до начала. не окружает −1, а поскольку P = 0, система замкнутого цикла стабильна для любого положительного усиления K, если мы заменим G(s)KG(s)[[FLT:

Система второго порядка с демпфированием: G(s) = ω n2/(s2 + 2ζω n s + ω n2)

Для системы второго порядка с естественной частотой ωn и коэффициентом демпфирования ζ, участок Nyquist начинается с (1, 0) для ω = 0 и заканчивается на (0, 0) фазой ω → ∞. Форма критически зависит от ζ< 1), the curve may exhibit a loop that crosses the negative real axis at a frequency ωpc, где фаза −180°.pcpcpc, тогда окружность участка Nyquist[1, система замкнутого цикла нестабильна.pc, то не происходит окружения и система стабильна.cpc = 1pc Бесконечная частота (асимптотически) если только система не имеет дополнительной динамики. Подождите, это важный нюанс: чистая система второго порядка без нулей имеет фазу, которая асимптотически приближается к −180°, но никогда не достигает −180° на любой конечной частоте. Поэтому ее участок никиста никогда не пересекает отрицательную реальную ось (кроме исходной). Таким образом, он никогда не окружает −1, и замкнутый цикл всегда стабилен для любого усиления — противоречивая типичному поведению второго порядка с демпфированием? На самом деле стандартная система второго порядка стабильна для всех положительных выигрышей, потому что ее фаза открытого цикла никогда не превышает −180° в конечном частотном диапазоне; она достигает только −180° на бесконечной частоте с нулевой величиной. Так что никакого окружения. Это иллюстрирует, что критерий никиста показывает стабильность для более сложных систем, таких как те, с задержками во времени или неминимальными фазовыми нулями, где фаза может превышать −180° на конечных частотах.

Чтобы увидеть правильный фазовый кроссовер, рассмотрим систему с дополнительным полюсом на низкой частоте: G(s) = 1/(s+1)(s2 + 2ζ s + 1)].] Здесь фаза может пересекаться −180° при некотором конечном ω, и участок может окружать −1, если коэффициент усиления достаточно высок.

Связь Nyquist, Bode и Nichols Plots

Плоскость Никвиста тесно связана с графиком Боде (величина и фаза) и диаграммой Николса. Плоскость Боде обеспечивает величину в дБ и фазу в градусах как функции частоты журнала. Плоскость Никвиста представляет собой параметрическую комбинацию величины и фазы: каждая точка на кривой Никвиста соответствует величине и фазе на заданном ω. ]прибыль является фактором, с помощью которого коэффициент усиления на частоте кроссовера фазы становится 1 (т.е. расстояние от −1 до кривой Никвиста вдоль оси отрицательного реального)] Фазовый предел является дополнительным фазовым отставанием, необходимым на частоте кроссовера усиления (], чтобы сделать фазу −180°. Эти поля непосредственно считываемы с графика Никвиста: предел усиления составляет 1/ | G (jω pc ) | (взаимный угол

Диаграмма Николса превращает график Найквиста в картезианский график величины (dB) против фазы (deg). Это преобразование полезно для проектирования контроллеров, поскольку коэффициент усиления и запас фазы легко видны, а отклик замкнутого цикла можно прочитать с помощью M-кругов. Все три представления эквивалентны математически; выбор одного зависит от контекста и предпочтений инженера.

Ограничения и расширенные соображения

Несмотря на свою мощь, график Nyquist имеет ограничения, которые инженеры должны признать. Во-первых, для систем с задержками во времени кривая Nyquist становится бесконечной и сильно изогнутой, требуя тщательной обработки бесконечного полукруглого отображения. Задержка добавляет фазовый сдвиг −ωT, который накапливается линейно, вызывая спираль сюжета Nyquist внутрь или наружу в зависимости от задержки. Анализ стабильности затем требует модифицированного критерия Nyquist, который учитывает задержку, используя концепцию контура Nyquist, который исключает полюс задержки в бесконечности.

Во-вторых, описанный график Nyquist предназначен для систем непрерывного времени. Для систем дискретного времени (z-домен) контур Nyquist становится единичным кругом, а критерий устойчивости исследует окружения −1 в w-плане или использование билинейного преобразования для отображения в s-план. Математические основы остаются аналогичными, но контур изменяется.

В-третьих, критерий Nyquist предполагает, что система замкнутого цикла является стабильной, если и только если окружность участка Nyquist −1 соответственно. Однако это требует, чтобы характер полюсов открытого цикла в RHP был известен a priori. Если функция переноса открытого цикла имеет полюса на воображаемой оси (включая чистые интеграторы), контур Nyquist должен быть изрезан, чтобы избежать этих полюсов, и отображение этих отступов производит бесконечные дуги, которые должны быть включены в число окружения. Это может быть сложно для инженеров, новых к методу.

Наконец, график Nyquist является инструментом частотной области; он обеспечивает необходимые и достаточные условия для стабильности только для линейных систем с инвариантным временем (LTI). Для нелинейных или изменяющихся во времени систем критерий Nyquist непосредственно не применяется, хотя для некоторых классов нелинейностей существуют расширения, такие как критерий окружности .

Заключение

Построение сюжета Nyquist глубоко укоренено в комплексном анализе, в частности, отображение контура Nyquist через функцию переноса открытого цикла и принцип аргумента. Овладев математическими шагами — оценивая G(jω)H(jω), выстраивая сложные точки и подсчитывая окружения −1 — инженеры получают мощный метод оценки стабильности замкнутого цикла и проектирования надежных контроллеров. Критерий Nyquist обеспечивает интуитивное понимание коэффициента усиления и пределов фазы, а его связь с графиками Bode и Nichols делает его универсальным инструментом в инструментальном наборе инженера-контролера. Понимание этих основ позволяет инженерам не только правильно применять сюжет Nyquist, но и интерпретировать его результаты с уверенностью в реальных приложениях, от аэрокосмических систем до управления промышленным процессом.

Для дальнейшего чтения обратитесь к стандартным учебникам, таким как статья Википедия о критерии стабильности Nyquist , Учебники по контролю для MATLAB и Simulink или классический текст Современная инженерия управления от Ogata. Эти ресурсы предоставляют дополнительные примеры и более глубокие математические выводы.