Математические модели в инженерии
Как интерпретировать фазу и величину в сюжетах Nyquist для инженеров
Table of Contents
Введение в Nyquist Plots в области управления
Никвистские графики являются краеугольным камнем анализа частотной области в области управления, предлагая графический метод оценки стабильности и производительности систем обратной связи. Разработанные Гарри Найквистом в 1930-х годах, эти графики отображают частотную реакцию системы & #8217;s функции передачи с открытым контуром G(s)H(s) на сложную плоскость. В отличие от бодовых графиков, которые разделяют величину и фазу на два графика, никвистские сюжеты объединяются в один полярный график, что облегчает визуализацию того, как усиление и фаза взаимодействуют на каждой частоте. Для инженеров, проектирующих контроллеры для самолетов, робототехники или управления процессом, освоение участков Найквиста необходимо для обеспечения надежной стабильности и предотвращения катастрофических колебаний.
В этом расширенном руководстве мы подробно рассмотрим интерпретацию величины и фазы в сюжетах Nyquist. Мы рассмотрим базовую математику, практические критерии стабильности, такие как критерий стабильности Nyquist, и как извлечь маржу усиления и фазовую маржу непосредственно из сюжета. Кроме того, мы обсудим общие неправильные интерпретации, взаимосвязь между сюжетами Nyquist и Bode и предоставим практические советы по применению этих знаний в реальных инженерных сценариях. К концу у вас будет полное понимание того, как читать и анализировать сюжеты Nyquist для разработки надежных систем управления.
Основы сюжетов Nyquist: Полярное представление частотного ответа
Никвистский график по существу является полярным графиком комплексной функции G(jω)H(jω), так как частота ω варьируется от 0 до ∞.Во многих учебниках сюжет Никвиста также включает симметричную часть для отрицательных частот, делая замкнутый контур. Сюжет рисуется в комплексной плоскости, где ось x (реальная ось) представляет реальную часть передаточной функции, а ось y (воображаемая ось) представляет мнимую часть. Каждая точка на кривой соответствует определенной частоте, а вектор от начала до этой точки кодирует и величину, и фазу.
Математически для данной частоты ω функция переноса с открытым контуром, оцененная при s = jω, даёт комплексное число:
G(jω)H(jω) = |G(jω)H(jω) | e^{jφ(ω)}
Здесь | G(jω)H(jω) | — величина (прибыль) на этой частоте, а φ(ω) — фазовый угол. На никистском графике величина — расстояние от начала до точки, а фаза — угол, измеряемый против часовой стрелки от положительной реальной оси. Как правило, кривая начинается с ω = 0 (что соответствует приросту dc, часто на положительной реальной оси) и заканчивается на ω = ∞ (что имеет тенденцию к происхождению, если система строго правильна).
Одно ключевое различие: Никвистские сюжеты часто путают с полярными сюжетами или квотами; Никвистские диаграммы; в других областях. В контрольной инженерии сюжет Никвиста конкретно относится к частотной реакции функции переноса открытого цикла, используемой для анализа стабильности по критерию Никвиста. Это не просто сюжет импеданса в электротехнике (хотя применяются аналогичные концепции). Для прочного фундамента обратитесь к классическим текстам, таким как это введение в Никвистские сюжеты из Control Systems Lab, который охватывает основы.
Толкование величины в сюжетах Nyquist
Величина на заданной частоте указывает, насколько система усиливает или ослабляет синусоидальный вход на этой частоте. На графике Найквиста величина представляет собой радиальное расстояние от начала. Для точки в координатах (x, y) величина является sqrt(x2 + y2). Эта величина также известна как усиление системы с открытым контуром.
Магнитуда, превышающая 1, означает, что система обеспечивает усиление на этой частоте; величина, меньшее чем 1, означает затухание. Критической точкой для анализа стабильности является величина на частоте, где фаза пересекает -180° (частота пересечения фазы). На этой частоте величина равна взаимному пределу усиления. В частности, если величина на фазе -180° меньше 1, система условно стабильна; если она больше 1, система может быть нестабильной.
Найти ещё сообщения от Nyquist Plot
Маржа прибыли (GM) является ключевой метрической стабильностью. Она определяется как фактор, с помощью которого прирост может быть увеличен до того, как система замкнутого цикла станет нестабильной. На графике Nyquist маржа прибыли обратно пропорциональна величине на частоте пересечения фаз (ω pc), где фаза составляет -180°. Математически:
GM = 1 / |G(jω pc)H(jω pc) |]
На графике это соответствует расстоянию от точки, где кривая пересекает отрицательную реальную ось до критической точки (-1, 0). Если кривая пересекает отрицательную реальную ось при -0,5 (величина 0,5), предел усиления составляет 2 (или 6 дБ. Если он пересекается при -2 (величина 2), предел усиления составляет 0,5, что указывает на то, что система уже нестабильна. Инженеры часто предпочитают предел усиления не менее 6 дБ (фактор 2) для прочной стабильности.]
Для визуализации представьте себе график Найквиста, приближающийся к точке -1. Чем ближе кривая приближается к -1, тем меньше предел усиления. Если кривая окружает или касается -1, система замкнутого цикла является незначительно стабильной или нестабильной.
Магнитуда при разных частотах
За пределами фазового кроссовера величина на низких частотах (прирост dc) важна для точности устойчивого состояния. Высокий прирост dc уменьшает погрешность устойчивого состояния для входов шага, но он также может влиять на поля стабильности, если не компенсируется. Магнитуда на высоких частотах откатывает из-за полюсов, поэтому строго правильные системы имеют тенденцию к происхождению как ω→∞. Форма величины на графике Nyquist не так непосредственно читаема, как на графиках Боде, но вы можете проследить радиальные расстояния вдоль кривой, чтобы измерить, как коэффициент усиления изменяется с частотой.
На практике инженеры часто накладывают на участок Найкиста единичный круг (магнитуда = 1). Пересечение кривой Найкиста с этим единичным кругом дает частоту кроссовера усиления (ω gc) — частоту, при которой величина равна 1. Эта частота используется для вычисления фазового поля.
Толкование фазы в сюжетах Nyquist
Фаза на графике Найквиста представлена углом вектора от начала до точки на кривой, измеренной от положительной реальной оси. Положительная фаза означает, что выход ведет вход (против часовой стрелки вращения), в то время как отрицательная фаза означает выходные лаги (оборота по часовой стрелке). Для большинства физических систем доминирует фазовый лаг, особенно на высоких частотах из-за задержек времени и полюсов.
Наиболее критическое значение фазы — -180° (или —π радианы).Когда фаза достигает -180°, обратная связь становится положительной, если величина также достаточно велика, что приводит к нестабильности. Частота, при которой фаза составляет -180°, называется частотой фазового кроссовера (ω pc). На этой частоте кривая Найквиста пересекает отрицательную реальную ось. Если это пересечение происходит слева от -1 на реальной оси (т.е. величина > 1), система неустойчива; если справа (величина < 1), система стабильна.
Маржа фазы из Nyquist Plot
Фазовый предел (PM) — это дополнительное фазовое отставание, необходимое для приведения системы к грани нестабильности (т.е. для достижения фазы -180° на частоте кроссовера усиления). Он измеряется на частоте, где величина равна 1 (прибавочная частота кроссовера, ω gc). На графике Nyquist фазовый край — угол между отрицательной реальной осью и линией от начала до точки, где кривая пересекает единичный круг. Точнее, если кривая Nyquist встречается с единичным кругом под углом φ (в градусах), то фазовый край:
PM = φ + 180°
Например, если кривая пересекает единичный круг на фазе -135°, фазовый предел составляет 45°. Фазовой предел менее 45° указывает на систему со значительным перерасходом в ответной стадии; фазовый предел, близкий к 0°, означает почти нестабильность. Хорошая инженерная практика обычно направлена на фазовый предел между 45° и 75°.]
Чтобы найти фазовый край на участке Nyquist непосредственно, нарисуйте круг радиуса 1, центрированный в начале. Где кривая Nyquist пересекает этот круг (впервые, как правило), измерьте угол от отрицательной реальной оси до этой точки. Это фазовый край. Многие инструменты управления проектированием автоматически вычисляют его.
Критерий стабильности никистской системы: от окружения к стабильности
Критерий устойчивости Nyquist является наиболее мощной особенностью участка Nyquist. Он связывает устойчивость системы замкнутого цикла с окружениями критической точки (-1, 0) на участке Nyquist функции переноса открытого цикла G(s)H(s). Критерий учитывает полюса открытого цикла в плоскости правой половины (RHP). Формальное утверждение:
Пусть N будет числом окружностей точки -1 по часовой стрелке на участке Nyquist G(s)H(s) в качестве s пересекает контур Nyquist (включая воображаемую ось и большой полукруг в правой полуплоскости).Z дается:
Z = N + PP — число полюсов открытой петли в правой половине плоскости (включая воображаемую ось).Z = 0, поэтому участок Nyquist должен окружать -1 в направлении против часовой стрелки столько раз, сколько существует полюсов RHP с открытой петлей.
На практике, если система с открытым контуром стабильна (P=0), система стабильна, если участок Найквиста не окружает -1. Если участок окружает -1 по часовой стрелке, система нестабильна. Если он окружает -1 против часовой стрелки, он стабилен только при наличии равного числа полюсов RHP.
Важно использовать правильное направление: участок Nyquist для положительных частот составляет лишь половину контура; инженеры обычно завершают участок зеркалированием для отрицательных частот, чтобы сформировать замкнутый контур. Окружения засчитываются для полного замкнутого контура. Для большинства практических систем (никаких полюсов RHP), просто проверяя, лежит ли -1 слева или справа от кривой Nyquist, поскольку частота увеличивается. Однако для систем с полюсами RHP или нулями необходим полный критерий. Подробные примеры можно найти в документации MathWorks по участкам Nyquist.
Окружения чтения из сюжета
Для определения окружений нарисуйте линию от точки -1 в любом направлении; подсчитайте, сколько раз кривая Найквиста пересекает эту линию. Направление (по часовой стрелке или против часовой стрелки) определяет знак. Альтернативно можно использовать мысленное квотирование; изгибающую цифровую цитату; подход. Для простой системы без полюсов RHP условие устойчивости состоит в том, что на участке Найквиста при прохождении в направлении возрастающей частоты не должно быть -1. Другими словами, точка -1 должна находиться вне критической области.
Практический анализ: комбинирование величины и фазы с использованием участков Nyquist
В то время как Боде рисует отдельные величины и фазы, графики Найквиста показывают и то, и другое одновременно, что полезно для визуализации того, как взаимодействуют усиление и фаза. Например, при проектировании компенсатора можно увидеть, как добавление свинцовой сети сдвигает кривую Найквиста вправо (улучшая фазовый предел) при увеличении величины. И наоборот, компенсация задержки тянет кривую к происхождению на высоких частотах, но может уменьшить фазовый запас на низких частотах.
Прибыль и фазовая маржа от одного участка
Одна из величайших практических ценностей сюжета Nyquist заключается в том, что вы можете прочитать как маржу усиления, так и фазовую маржу от него, при условии, что вы знаете единичный круг и отрицательную реальную ось. Шаги:
- Определить пересечение кривой Найквиста с отрицательной реальной осью (фаза = -180°). Магнитуда в этой точке дает коэффициент усиления, как описано ранее.
- Определить пересечение с единичным кругом (магнитность = 1). Угол этой точки (измеряется по положительной реальной оси) дает фазовый запас: PM = угол + 180° (если угол отрицательный).
- Если кривая Найквиста не пересекает единичный круг (то есть величина никогда не равна 1), то система имеет бесконечный фазовый запас, но вам все равно нужно проверить запас усиления.
Важно: Для систем с несколькими пересечениями наиболее критичным пересечением определяется запас. Всегда проверяйте пересечение, наиболее близкое к точке -1.
Распространенные ошибки в интерпретации сюжетов Nyquist
- Сбивающие с толку величины и фазовые векторы: Помните, что длина вектора — величина; угол — фаза. Новички часто неверно истолковывают расстояние до реальной оси как величину. Используйте сетку или шкалу.
- Использование только положительной частотной ветви: Критерий стабильности требует всего контура, включая зеркальное отражение отрицательных частот. Многие программные средства показывают только положительную частотную ветвь; нужно мысленно завершить замкнутый цикл. Если график для отрицательных частот, он симметричен относительно реальной оси.
- Игнорирование полюсов с открытым контуром на воображаемой оси: Системы с интеграторами (поли при s=0) заставляют сюжет Никвиста идти к бесконечности на низких частотах.Контур Никвиста должен отступить вокруг этих полюсов, что приводит к дугам на бесконечности, которые должны учитываться при подсчете окружения.
- Задержки во времени при взгляде:] Задержки во времени добавляют дополнительное фазовое отставание без влияния на величину. Это может привести к тому, что участок Никиста будет закручиваться внутрь, потенциально вызывая нестабильность на высоких частотах. Сюжет Никиста по-прежнему действителен, но поведение бесконечной частоты должно быть тщательно обработано.
Использование программных средств для анализа Nyquist
Современный дизайн управления выполняется с использованием вычислительных инструментов, таких как MATLAB, Python (библиотека управления) или Scilab. Эти инструменты автоматически генерируют графики Nyquist и могут вычислять коэффициенты усиления и фазовые поля с помощью одной команды. Однако понимание базовой интерпретации жизненно важно для проверки результатов и устранения непредвиденной нестабильности. Например, функция MATLAB’s отображает положительный частотный диапазон; команда возвращает коэффициент усиления и запас фаз для системы с открытым циклом. Библиотека управления Python’ обеспечивает аналогичную функциональность. При использовании этих инструментов всегда дважды проверяйте направление окружения и помните, что график может показывать только частоты от 0 до ∞; чтобы увидеть полный контур, вам может потребоваться вручную добавить зеркальное отражение. Всестороннее руководство по использованию Python для анализа Nyquist можно найти в документации библиотеки систем управления Python.
Относитесь к Nyquist Plots с точки зрения бод-плотов
Никвист и Боде графики дополняют друг друга. Боде графики предоставляют четкую информацию о величине и фазе отдельно в широком диапазоне частот, что облегчает идентификацию резонансных пиков, скорости выкатывания и частот кроссовера. Найквистские графики сжимают всю частотную информацию в единую кривую, подчеркивая критическую точку -1 и общую форму. Инженеры часто используют оба: начните с Боде для разработки компенсатора, затем проверьте стабильность с Найквистом. Они связаны через концепцию комплексного картирования плоскости. Например, частота кроссовера усиления, идентифицированная с Боде графика величины, точно соответствует точке, где частота кроссовера фазы от Боде соответствует Nyquist график & #8217;с пересечением с отрицательной реальной осью. Для получения дополнительной информации об этой связи см. этот университетский ресурс по критерию Найквиста .
Расширенные темы: чувствительность, дополнительная чувствительность и дисковые границы
В то время как коэффициент усиления и предел фаз являются традиционными метрическими показателями, современная теория надежного управления использует дополнительные меры, такие как пики чувствительности. Найквистский график также может использоваться для получения функции чувствительности S(s) = 1/(1 + G(s)H(s)]. Максимальная чувствительность Ms является максимальным расстоянием от -1 до любой точки на кривой Найквиста. Меньшее максимальное расстояние подразумевает лучшую устойчивость. Функция дополнительной чувствительности T(s) относится к производительности. На рисунке Найквиста также можно визуализировать дисковые поля, которые одновременно рассматривают изменения усиления и фазы, как область вокруг точки -1. Инженеры в аэрокосмической и автомобильной областях часто используют эти передовые метрики для обеспечения надежности системы при вариациях параметров.
Для вычисления дисковых полей из участка Nyquist представьте диск, центрированный на -1 с определенным радиусом; если кривая Nyquist остается за пределами этого диска, система может выдержать определенное количество одновременного усиления и фазовой неопределенности, не становясь нестабильной. Это более всеобъемлющая мера, чем отдельные коэффициенты усиления и фазовые погрешности.
Практический пример: система второго заказа
Рассмотрим простую функцию переноса с открытым контуром: G(s)H(s) = 1/(s(s+1)]. Это система типа 1 (один интегратор). Схема Nyquist для положительных частот начинается в бесконечности, когда ω→0 (из-за полюса при s=0), затем спираль к началу на высоких частотах. Фаза начинается при -90° (из-за интегратора) и приближается к -180° как ω→∞. Для этой системы кроссовер происходит при величине менее 1, поэтому система замкнутого контура стабильна. Погрешность усиления бесконечна (поскольку на самом деле график никогда не достигает -1 по величине, потому что он пересекает отрицательную реальную ось справа от -1). Фазовой предел можно найти там, где величина равна 1. Решая |G(jω)H(jω) | = 1, мы находим ω ≈ 0,618 рад/с, и фаза на этой частоте составляет около -128°, давая фазов
Теперь измените систему: добавьте коэффициент усиления 5: G(s)H(s) = 5/(s(s+1)]. Участок Nyquist будет масштабироваться на 5. Пересечение с отрицательной реальной осью теперь происходит при величине 5 × (что-то большее 1) — фактически величина точки пересечения становится в 5 раз больше предыдущей величины пересечения. Если этот продукт превышает 1, система становится нестабильной. В этом случае величина пересечения для исходной системы была 0,5, поэтому с коэффициентом усиления 5 она становится 2,5, то есть участок Nyquist окружен -1, а система замкнутого цикла нестабильна. Это иллюстрирует прямую связь между коэффициентом усиления и стабильностью.
Оригинальное название: Mastering Nyquist Plot Interpretation
Интерпретация фазы и величины на участках Nyquist является критическим навыком для каждого инженера управления. Понимая, что величина является радиальным расстоянием от начала и фазы является углом от положительной реальной оси, инженеры могут извлекать как выигрыш, так и фазовые поля непосредственно из участка. Критерий стабильности Nyquist обеспечивает надежный метод оценки стабильности замкнутого цикла, даже для систем с нестабильностью открытого цикла. Сочетание анализа Nyquist с участками Боде и современными инструментами позволяет осуществлять комплексный дизайн и устранение неполадок. По мере того, как вы практикуете интерпретацию более сложных участков Nyquist - особенно тех, у которых есть задержки во времени, неминимальные нулевые фазы или множественные зависящие от частоты элементы - ваша интуиция для поведения в частотной области улучшится. Всегда помните критическую точку (-1,0) и используйте единичный круг в качестве руководства. С этими навыками вы можете проектировать системы управления, которые являются стабильными и эффективными в широком диапазоне условий работы.
Для дальнейшего чтения изучите оригинальную статью H. Nyquist, quot;Regeneration Theoryquot; (1932), или современные учебники, такие как Современная инженерия управления от Ogata и Обратная связь управления динамическими системами Франклина, Пауэлла и Эмами-Наейни.