12.2: Критерій стабільності Найквіста
Критерій Найквіста - це графічна техніка для визначення того, чи можна стабілізувати нестабільну лінійну інваріантну систему часу за допомогою петлі негативного зворотного зв'язку. Ми трохи уважніше розглянемо такі системи, коли вивчимо перетворення Лапласа в наступній темі. Для цієї теми ми задовольнимося постановкою проблеми лише з найменшим фізичним контекстом.
Ви вже стикалися з лінійними інваріантними системами часу в 18.03 (або його еквіваленті), коли ви вирішували лінійні диференціальні рівняння з постійним коефіцієнтом.
Системні функції
Лінійна інваріантна система часу має системну функцію, яка є функцією комплексної змінної. Зазвичай комплексна змінна позначаєтьсяs і для системної функції використовується велика літера.
НехайG(s) буде така системна функція. Зробимо стандартне припущення, якеG(s) є мероморфним з скінченною кількістю (скінченних) полюсів. Це припущення тримається в багатьох цікавих випадках. Наприклад, досить частоG(s) є раціональною функцієюQ(s)/P(s) (QіP є поліномами).
Нас буде турбувати стабільність роботи системи.
Система з системною функцієюG(s) називається стабільною, якщо всі полюсиG знаходяться в лівій півплощині. Тобто, якщо всі полюсиG мають негативну реальну частину.
Система називається нестійкою, якщо будь-які полюси знаходяться в правій півплощині, тобто мають позитивну дійсну частину.
Для крайнього випадку, коли жоден полюс не має позитивної реальної частини, але деякі чисто уявні ми будемо називати систему незначно стабільною. Цей випадок можна проаналізувати за допомогою наших методик. Для наших цілей буде потрібно і контур з відступом уздовж уявної осі. Якщо у нас буде час, ми зробимо аналіз.
ЧиG(s)=s(s+2)(s2+4s+5) стабільна система з системною функцією?
Рішення
Полюси є−2,−2±i. Так як всі вони знаходяться в лівій півплощині, система стабільна.
ЧиG(s)=s(s2−4)(s2+4s+5) стабільна система з системною функцією?
Рішення
Полюси є±2,−2±i. Так як один полюс знаходиться в правій півплощині, система нестабільна.
ЧиG(s)=s(s+2)(s2+4) стабільна система з системною функцією?
Рішення
Полюси є−2,±2i. У правій півплощині немає полюсів. Так як на уявній осі є полюси, система незначно стійка.
Термінологія. Поки ми були обережні, щоб сказати «система з системною функцієюG(s)». Відтепер ми дозволимо собі бути трохи більш випадковими і сказати «системуG(s)». Вона абсолютно ясна і скочується з мови трохи легше!
Діаграми полюс-нуль
Ми можемо візуалізуватиG(s) за допомогою діаграми полюс-нуль. Це діаграма вs -площині, де ми ставимо невеликий хрест на кожному полюсі і мале коло на кожному нулі.
Дайте діаграми нульових полюсів для кожної з систем
G1(s)=s(s+2)(s2+4s+5), G1(s)=s(s2−4)(s2+4s+5), G1(s)=s(s+2)(s2+4)
Рішення
Це ті ж системи, що і в прикладах трохи вище. Спочатку відзначимо, що всі вони мають єдиний нуль у початку. Таким чином, ми ставимо коло біля початку і хрест на кожному полюсі.
Діаграми полюс-нуль для трьох систем.
трохи про стабільність
Це просто для того, щоб дати вам трохи фізичної орієнтації. З огляду на наше визначення стійкості вище, ми могли б, в принципі, обговорювати стабільність без найменшого уявлення про те, що вона означає для фізичних систем.
ПолюсиG(s) відповідають тому, що називаються режимами роботи системи. Простий полюс приs1 відповідає режимуy1(t)=es1t. Система стабільна, якщо всі режими затухають до 0, тобто якщо полюси знаходяться в лівій півплощині.
Фізично режими говорять нам про поведінку системи, коли вхідний сигнал дорівнює 0, але є початкові умови. Полюс з позитивною реальною частиною відповідав би режиму, який переходить до нескінченності уt міру зростання. Безумовно, розумно назвати систему, яка робить це у відповідь на нульовий сигнал (часто називають «без входу») нестабільною.
Щоб зв'язати це з 18.03: якщо система моделюється диференціальним рівнянням, режими відповідають однорідним розв'язкамy(t)=est, деs є корінь характеристичного рівняння. У 18.03 ми назвали систему стабільною, якщо кожен однорідний розчин розпадається до 0. Тобто, якщо непримусова система завжди влаштовувалася до рівноваги.
Системи замкнутого циклу
Якщо система з системною функцієюG(s) нестабільна, її іноді можна стабілізувати за допомогою так званого петлі негативного зворотного зв'язку. Нова система називається замкнутою системою. Його системна функція задається формулою Блек
GCL(s)=G(s)1+kG(s),
деk називається коефіцієнтом зворотного зв'язку. Ми просто приймемо цю формулу. Будь-який клас або книга з теорії управління виведе його для вас.
У цьому контекстіG(s) називається системною функцією відкритого циклу.
ОскількиGCL це системна функція, ми можемо запитати, чи стабільна система.
Полюси функції системи замкнутого контуру,GCL(s) наведені в Рівнянні 12.3.2, є нулями1+kG(s).
- Доказ
-
Дивлячись на Рівняння 12.3.2, є два можливі джерела полюсів дляGCL.
- Нулі знаменника1+kG. Теорема визнає їх.
- Полюси зG. Так якG є і в чисельнику, і в знаменникуGCL його повинно бути зрозуміло, що полюси скасовуються. Ми можемо показати це формально, використовуючи серію Laurent. ЯкщоG має полюс порядкуn вs0 то
G(s)=1(s−s0)n(bn+bn−1(s−s0)+ ...a0(s−s0)n+a1(s−s0)n+1+ ...),
деbn≠0. Отже,
GCL(s)=1(s−s0)n(bn+bn−1(s−s0)+ ...a0(s−s0)n+ ...)1+k(s−s0)n(bn+bn−1(s−s0)+ ...a0(s−s0)n+ ...)=(bn+bn−1(s−s0)+ ...a0(s−s0)n+ ...)(s−s0)n+k(bn+bn−1(s−s0)+ ...a0(s−s0)n+ ...)
який чітко аналітичний вs0. (Уs0 нього дорівнюєbn/(kbn)=1/k.)
Встановіть коефіцієнт зворотного зв'язкуk=1. Припустимо,a це реально, для яких значеньa є система відкритого циклуG(s)=1s+a стабільна? При яких значенняхa відповідна замкнута системаGCL(s) стабільна?
(У цьому прикладі немає особливої причини, якаa повинна бути реальною. Але у фізичних системах складні полюси, як правило, приходять у сполучених парах.)
Рішення
G(s)має один полюс вs=−a. Таким чином, він стійкий, коли полюс знаходиться в лівій півплощині, тобто дляa>0.
Функція системи замкнутого циклу
GCL(s)=1/(s+a)1+1/(s+a)=1s+a+1.
Це має полюс наs=−a−1, так що він стабільний, якщоa>−1. Контур зворотного зв'язку стабілізував нестабільні системи з розімкнутим контуром−1<a≤0. (Власне, дляa=0 розімкнутого контуру незначно стабільний, але він повністю стабілізується замкнутим контуром.)
Алгебра, яка бере участь у скасуванніs+a терміну в знаменниках, є саме скасуванням, яке робить полюсиG знімних сингулярностей вGCL.
ПрипустимоG(s)=s+1s−1. Чи стабільна система відкритого контуру? Чи стабільна система замкнутого контуру, колиk=2.
Рішення
G(s)має полюс у правій півплощині, тому система відкритого контуру не є стабільною. Функція системи замкнутого циклу
GCL(s)=G1+kG=(s+1)/(s−1)1+2(s+1)/(s−1)=s+13s+1.
Єдиний полюс знаходиться наs=−1/3, тому система замкнутого контуру стабільна. Це випадок, коли зворотний зв'язок стабілізувала нестабільну систему.
G(s)=s−1s+1. Чи стабільна система відкритого контуру? Чи стабільна система замкнутого контуру, колиk=2.
Рішення
Єдина ділянкаG(s) знаходиться в лівій півплощині, тому система відкритого контуру стабільна. Функція системи замкнутого циклу
GCL(s)=G1+kG=(s−1)/(s+1)1+2(s−1)/(s+1)=s−13s−1.
Це має один полюс наs=1/3, тому система замкнутого контуру нестабільна. Це випадок, коли зворотний зв'язок дестабілізував стабільну систему. Це може статися!
Ділянки Найквіста
Для графіка і критерію Найквіста кривою завждиγ буде уявнаs -вісь. Тобто
s=γ(ω)=iω, where −∞<ω<∞.
Для системиG(s) таk коефіцієнта зворотного зв'язку графік Найквіста - це графік кривої
w=kG∘γ(ω)=kG(iω).
Тобто сюжет Найквіста - це зображення уявної осі під картоюw=kG(s).
Уγ(ω) змінній є грецька омега і вw=G∘γ нас є подвійний u.
НехайG(s)=1s+1. Намалюйте сюжет Nyquist сk=1.
Рішення
У випадкуG(s) є дробовим лінійним перетворенням, тому ми знаємо, що воно відображає уявну вісь на коло. Легко перевірити це коло через початок з центромw=1/2. Ви також можете перевірити, чи проходить вона за годинниковою стрілкою.
Нюквіст сюжетG(s)=1/(s+1), сk=1.
ВізьмемоG(s) з попереднього прикладу. Опишіть сюжет Найквіста з коефіцієнтом посиленняk=2.
Рішення
Сюжет Найквіста - це графікkG(iω). k=2Коефіцієнт масштабує коло в попередньому прикладі на 2. Тобто сюжет Nyquist - це коло через початок з центромw=1.
Загалом, коефіцієнт зворотного зв'язку буде якраз масштабувати сюжет Nyquist.
Критерій Найквіста
Критерій Найквіста дає графічний метод перевірки стійкості системи замкнутого контуру.
Припустимо, щоG(s) має скінченну кількість нулів і полюсів у правій півплощині. Також припустимо, щоG(s) розпадається до 0, якs переходить до нескінченності. Тоді система замкнутого контуру зk коефіцієнтом зворотного зв'язку стабільна тоді і тільки тоді, коли номер обмотки графіка Найквіста навколоw=−1 дорівнює числу полюсівG(s) у правій півплощині.
Більш коротко,
GCL(s) is stable ⇔ Ind(kG∘γ,−1)=PG,RHP
Ось уявнаs -вісь іPG,RHP кількість полюсів вихідної функції системи розімкнутого циклуG(s) в правій півплощині.γ
- Доказ
-
GCLстійкий саме тоді, коли всі його полюси знаходяться в лівій півплощині. Тепер нагадаємо, щоGCL полюси - це саме нулі1+kG. Отже, стабільністьGCL - це саме умова, що кількість нулів1+kG у правій півплощині дорівнює 0.
Давайте попрацюємо зі знайомим контуром.
НехайγR=C1+CR. Зверніть увагу,γR що проходить вclockwise напрямку. Виберіть доситьR великий, щоб (кінцеве число) полюсів і нулівG у правій півплощині були всерединіγR. Тепер ми можемо застосувати Equation 12.2.4 в наслідок принципу аргументуγ доkG(s) і отримати
−Ind(kG∘γR,−1)=Z1+kG,γR−PG,γR
(Знак мінус - через напрямок кривої за годинниковою стрілкою.) Таким чином, для всіх великихR
the system is stable ⇔ Z1+kG,γR=0 \Leftrightarow Ind(kG∘γR,−1)=PG,γR
Нарешті, ми можемоR відпустити до нескінченності. Припущення, щоG(s) розпадається 0 до якs йде,∞ означає, що в межі вся криваkG∘CR стає єдиною точкою на початку. Так в межіkG∘γR стаєkG∘γ. QED
Приклади використання математики Nyquist Plot
Критерій Nyquist - це візуальний метод, який вимагає певного способу створення сюжету Nyquist. Для цього ми будемо використовувати один з MIT Mathlets (трохи модифікований для наших цілей). Відкрийте аплет сюжету Nyquist за адресою
http://web.mit.edu/jorloff/www/jmoapplets/nyquist/nyquistCrit.html
Пограйте з аплетом, читайте довідку.
Тепер оновіть браузер, щоб відновити аплет в початковий стан. ВстановітьFormula прапорець. Формула є простим способом зчитування значень полюсів і нулівG(s). У початковому стані аплет повинен мати нуль ats=1 і полюси вs=0.33±1.75i.
Графік лівої руки - це діаграма полюс-нуль. Графік правої руки - це сюжет Nyquist.
Щоб відчути сюжет Nyquist. Подивіться на діаграму полюсів і за допомогою миші перетягніть жовту точку вгору і вниз по уявній осі. Його зображення підkG(s) простежить сюжет Нікіс.
Зверніть увагу, що коли жовта точка знаходиться на будь-якому кінці осі, її зображення на графіку Nyquist близьке до 0.
Оновіть сторінку, щоб повернути нуль і полюси в початковий стан. У правій півплощині є два полюси, тому система відкритого контуруG(s) нестабільна. Зk=1, який номер обмотки ділянки Nyquist навколо -1? Чи стабільна система замкнутого циклу?
Рішення
Крива накручується двічі навколо -1 в напрямку проти годинникової стрілки, тому число намотуванняInd(kG∘γ,−1)=2. Оскільки кількість полюсівG у правій півплощині таке ж, як це число обмоток, система замкнутого контуру стабільна.
З однаковими полюсами і нулями перемістітьk повзунок і визначте, який діапазонk робить систему замкнутого циклу стабільною.
Рішення
Колиk невеликий, сюжет Nyquist має обмотку номер 0 навколо -1. Для цих значеньk,GCL є нестабільним. k=0.7У міруk збільшення десь міжk=0.65 і число обмоток скаче з 0 на 2 і система замкнутого контуру стає стабільною. Це триває до тих пірk, поки не буде між 3. 10 і 3. 20, в цей момент число обмотки стає 1 іGCL стає нестійким.
Відповідь: Система замкнутого циклу стабільнаk (приблизно) між 0. 7 і 3. 10.
У попередній задачі ви могли б аналітично визначити діапазонk,GCL(s) де стабільний?
Рішення
Так! Це можливо для невеликих систем. Це складніше для систем вищого порядку, але є методи, які не вимагають обчислення полюсів. У цьому випадку ми маємо
GCL(s)=G(s)1+kG(s)=s−1(s−0.33)2+1.7521+k(s−1)(s−0.33)2+1.752=s−1(s−0.33)2+1.752+k(s−1)
Отже, полюси - коріння
(s−0.33)2+1.752+k(s−1)=s2+(k−0.66)s+0.332+1.752−k
Для квадратичного з додатними коефіцієнтами коріння мають негативну дійсну частину. Це відбувається, коли
0.66<k<0.332+1.752≈3.17.
Що відбувається, колиk переходить до 0.
Рішення
Щоk стосується 0, графік Nyquist скорочується до однієї точки на початку. У цьому випадку число обмоток навколо -1 дорівнює 0, а критерій Найквіста говорить, що система замкнутого контуру стабільна тоді і лише тоді, коли система відкритого контуру стабільна.
Це має мати сенс, так як зk=0,
GCL=G1+kG=G.
Складіть систему з наступними нулями і полюсами:
- Пара нулів в0.6±0.75i.
- Пара полюсів в−0.5±2.5i.
- Один полюс на 0,25.
Чи стабільна відповідна система замкнутого циклу, колиk=6?
Рішення
Відповідь - ні, неGCL є стабільним. Gмає один полюс в правій половині площини. Матлет показує сюжет Nyquist вітрів один раз навколоw=−1 вclockwise напрямку. Таким чином, число обмотки дорівнює -1, що не дорівнює числу полюсівG в правій півплощині.
Якщо ми встановимоk=3, система замкнутого контуру стабільна.