Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
LibreTexts - Ukrayinska

12.2: Критерій стабільності Найквіста

Критерій Найквіста - це графічна техніка для визначення того, чи можна стабілізувати нестабільну лінійну інваріантну систему часу за допомогою петлі негативного зворотного зв'язку. Ми трохи уважніше розглянемо такі системи, коли вивчимо перетворення Лапласа в наступній темі. Для цієї теми ми задовольнимося постановкою проблеми лише з найменшим фізичним контекстом.

Примітка

Ви вже стикалися з лінійними інваріантними системами часу в 18.03 (або його еквіваленті), коли ви вирішували лінійні диференціальні рівняння з постійним коефіцієнтом.

Системні функції

Лінійна інваріантна система часу має системну функцію, яка є функцією комплексної змінної. Зазвичай комплексна змінна позначаєтьсяs і для системної функції використовується велика літера.

НехайG(s) буде така системна функція. Зробимо стандартне припущення, якеG(s) є мероморфним з скінченною кількістю (скінченних) полюсів. Це припущення тримається в багатьох цікавих випадках. Наприклад, досить частоG(s) є раціональною функцієюQ(s)/P(s) (QіP є поліномами).

Нас буде турбувати стабільність роботи системи.

Визначення

Система з системною функцієюG(s) називається стабільною, якщо всі полюсиG знаходяться в лівій півплощині. Тобто, якщо всі полюсиG мають негативну реальну частину.

Система називається нестійкою, якщо будь-які полюси знаходяться в правій півплощині, тобто мають позитивну дійсну частину.

Для крайнього випадку, коли жоден полюс не має позитивної реальної частини, але деякі чисто уявні ми будемо називати систему незначно стабільною. Цей випадок можна проаналізувати за допомогою наших методик. Для наших цілей буде потрібно і контур з відступом уздовж уявної осі. Якщо у нас буде час, ми зробимо аналіз.

Приклад12.2.1

ЧиG(s)=s(s+2)(s2+4s+5) стабільна система з системною функцією?

Рішення

Полюси є2,2±i. Так як всі вони знаходяться в лівій півплощині, система стабільна.

Приклад12.2.2

ЧиG(s)=s(s24)(s2+4s+5) стабільна система з системною функцією?

Рішення

Полюси є±2,2±i. Так як один полюс знаходиться в правій півплощині, система нестабільна.

Приклад12.2.3

ЧиG(s)=s(s+2)(s2+4) стабільна система з системною функцією?

Рішення

Полюси є2,±2i. У правій півплощині немає полюсів. Так як на уявній осі є полюси, система незначно стійка.

Термінологія. Поки ми були обережні, щоб сказати «система з системною функцієюG(s)». Відтепер ми дозволимо собі бути трохи більш випадковими і сказати «системуG(s)». Вона абсолютно ясна і скочується з мови трохи легше!

Діаграми полюс-нуль

Ми можемо візуалізуватиG(s) за допомогою діаграми полюс-нуль. Це діаграма вs -площині, де ми ставимо невеликий хрест на кожному полюсі і мале коло на кожному нулі.

Приклад12.2.4

Дайте діаграми нульових полюсів для кожної з систем

G1(s)=s(s+2)(s2+4s+5),   G1(s)=s(s24)(s2+4s+5),   G1(s)=s(s+2)(s2+4)

Рішення

Це ті ж системи, що і в прикладах трохи вище. Спочатку відзначимо, що всі вони мають єдиний нуль у початку. Таким чином, ми ставимо коло біля початку і хрест на кожному полюсі.

2020-09-15 пнг
Діаграми полюс-нуль для трьох систем.

трохи про стабільність

Це просто для того, щоб дати вам трохи фізичної орієнтації. З огляду на наше визначення стійкості вище, ми могли б, в принципі, обговорювати стабільність без найменшого уявлення про те, що вона означає для фізичних систем.

Полюси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 це системна функція, ми можемо запитати, чи стабільна система.

Теорема12.2.1

Полюси функції системи замкнутого контуру,GCL(s) наведені в Рівнянні 12.3.2, є нулями1+kG(s).

Доказ

Дивлячись на Рівняння 12.3.2, є два можливі джерела полюсів дляGCL.

  1. Нулі знаменника1+kG. Теорема визнає їх.
  2. Полюси зG. Так якG є і в чисельнику, і в знаменникуGCL його повинно бути зрозуміло, що полюси скасовуються. Ми можемо показати це формально, використовуючи серію Laurent. ЯкщоG має полюс порядкуn вs0 то

G(s)=1(ss0)n(bn+bn1(ss0)+ ...a0(ss0)n+a1(ss0)n+1+ ...),

деbn0. Отже,

GCL(s)=1(ss0)n(bn+bn1(ss0)+ ...a0(ss0)n+ ...)1+k(ss0)n(bn+bn1(ss0)+ ...a0(ss0)n+ ...)=(bn+bn1(ss0)+ ...a0(ss0)n+ ...)(ss0)n+k(bn+bn1(ss0)+ ...a0(ss0)n+ ...)

який чітко аналітичний вs0. (Уs0 нього дорівнюєbn/(kbn)=1/k.)

Приклад12.2.5

Встановіть коефіцієнт зворотного зв'язку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=a1, так що він стабільний, якщоa>1. Контур зворотного зв'язку стабілізував нестабільні системи з розімкнутим контуром1<a0. (Власне, дляa=0 розімкнутого контуру незначно стабільний, але він повністю стабілізується замкнутим контуром.)

Примітка

Алгебра, яка бере участь у скасуванніs+a терміну в знаменниках, є саме скасуванням, яке робить полюсиG знімних сингулярностей вGCL.

Приклад12.2.6

ПрипустимоG(s)=s+1s1. Чи стабільна система відкритого контуру? Чи стабільна система замкнутого контуру, колиk=2.

Рішення

G(s)має полюс у правій півплощині, тому система відкритого контуру не є стабільною. Функція системи замкнутого циклу

GCL(s)=G1+kG=(s+1)/(s1)1+2(s+1)/(s1)=s+13s+1.

Єдиний полюс знаходиться наs=1/3, тому система замкнутого контуру стабільна. Це випадок, коли зворотний зв'язок стабілізувала нестабільну систему.

Приклад12.2.7

G(s)=s1s+1. Чи стабільна система відкритого контуру? Чи стабільна система замкнутого контуру, колиk=2.

Рішення

Єдина ділянкаG(s) знаходиться в лівій півплощині, тому система відкритого контуру стабільна. Функція системи замкнутого циклу

GCL(s)=G1+kG=(s1)/(s+1)1+2(s1)/(s+1)=s13s1.

Це має один полюс наs=1/3, тому система замкнутого контуру нестабільна. Це випадок, коли зворотний зв'язок дестабілізував стабільну систему. Це може статися!

Ділянки Найквіста

Для графіка і критерію Найквіста кривою завждиγ буде уявнаs -вісь. Тобто

s=γ(ω)=iω, where <ω<.

Для системиG(s) таk коефіцієнта зворотного зв'язку графік Найквіста - це графік кривої

w=kGγ(ω)=kG(iω).

Тобто сюжет Найквіста - це зображення уявної осі під картоюw=kG(s).

Примітка

Уγ(ω) змінній є грецька омега і вw=Gγ нас є подвійний u.

Приклад12.2.8

НехайG(s)=1s+1. Намалюйте сюжет Nyquist сk=1.

Рішення

У випадкуG(s) є дробовим лінійним перетворенням, тому ми знаємо, що воно відображає уявну вісь на коло. Легко перевірити це коло через початок з центромw=1/2. Ви також можете перевірити, чи проходить вона за годинниковою стрілкою.

2020-09-15 4.30.21png
Нюквіст сюжетG(s)=1/(s+1), сk=1.

Приклад12.2.9

ВізьмемоG(s) з попереднього прикладу. Опишіть сюжет Найквіста з коефіцієнтом посиленняk=2.

Рішення

Сюжет Найквіста - це графікkG(iω). k=2Коефіцієнт масштабує коло в попередньому прикладі на 2. Тобто сюжет Nyquist - це коло через початок з центромw=1.

Загалом, коефіцієнт зворотного зв'язку буде якраз масштабувати сюжет Nyquist.

Критерій Найквіста

Критерій Найквіста дає графічний метод перевірки стійкості системи замкнутого контуру.

Теорема12.2.2 Nyquist criterion

Припустимо, що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.

Давайте попрацюємо зі знайомим контуром.

2020-09-15 4.39.10.PNG

НехайγR=C1+CR. Зверніть увагу,γR що проходить вclockwise напрямку. Виберіть доситьR великий, щоб (кінцеве число) полюсів і нулівG у правій півплощині були всерединіγR. Тепер ми можемо застосувати Equation 12.2.4 в наслідок принципу аргументуγ доkG(s) і отримати

Ind(kGγR,1)=Z1+kG,γRPG,γR

(Знак мінус - через напрямок кривої за годинниковою стрілкою.) Таким чином, для всіх великихR

the system is stable  Z1+kG,γR=0 \Leftrightarow  Ind(kGγR,1)=PG,γR

Нарешті, ми можемоR відпустити до нескінченності. Припущення, щоG(s) розпадається 0 до якs йде, означає, що в межі вся криваkGCR стає єдиною точкою на початку. Так в межі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.

Приклад12.2.10

Щоб відчути сюжет Nyquist. Подивіться на діаграму полюсів і за допомогою миші перетягніть жовту точку вгору і вниз по уявній осі. Його зображення підkG(s) простежить сюжет Нікіс.

Зверніть увагу, що коли жовта точка знаходиться на будь-якому кінці осі, її зображення на графіку Nyquist близьке до 0.

Приклад12.2.11

Оновіть сторінку, щоб повернути нуль і полюси в початковий стан. У правій півплощині є два полюси, тому система відкритого контуруG(s) нестабільна. Зk=1, який номер обмотки ділянки Nyquist навколо -1? Чи стабільна система замкнутого циклу?

Рішення

Крива накручується двічі навколо -1 в напрямку проти годинникової стрілки, тому число намотуванняInd(kGγ,1)=2. Оскільки кількість полюсівG у правій півплощині таке ж, як це число обмоток, система замкнутого контуру стабільна.

Приклад12.2.12

З однаковими полюсами і нулями перемістіть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.

Приклад12.2.13

У попередній задачі ви могли б аналітично визначити діапазонk,GCL(s) де стабільний?

Рішення

Так! Це можливо для невеликих систем. Це складніше для систем вищого порядку, але є методи, які не вимагають обчислення полюсів. У цьому випадку ми маємо

GCL(s)=G(s)1+kG(s)=s1(s0.33)2+1.7521+k(s1)(s0.33)2+1.752=s1(s0.33)2+1.752+k(s1)

Отже, полюси - коріння

(s0.33)2+1.752+k(s1)=s2+(k0.66)s+0.332+1.752k

Для квадратичного з додатними коефіцієнтами коріння мають негативну дійсну частину. Це відбувається, коли

0.66<k<0.332+1.7523.17.

Приклад12.2.14

Що відбувається, колиk переходить до 0.

Рішення

Щоk стосується 0, графік Nyquist скорочується до однієї точки на початку. У цьому випадку число обмоток навколо -1 дорівнює 0, а критерій Найквіста говорить, що система замкнутого контуру стабільна тоді і лише тоді, коли система відкритого контуру стабільна.

Це має мати сенс, так як зk=0,

GCL=G1+kG=G.

Приклад12.2.15

Складіть систему з наступними нулями і полюсами:

  • Пара нулів в0.6±0.75i.
  • Пара полюсів в0.5±2.5i.
  • Один полюс на 0,25.

Чи стабільна відповідна система замкнутого циклу, колиk=6?

Рішення

Відповідь - ні, неGCL є стабільним. Gмає один полюс в правій половині площини. Матлет показує сюжет Nyquist вітрів один раз навколоw=1 вclockwise напрямку. Таким чином, число обмотки дорівнює -1, що не дорівнює числу полюсівG в правій півплощині.

Якщо ми встановимоk=3, система замкнутого контуру стабільна.