Бюджетний графік: Читання стабільності за допомогою частотного перегляду

Перетворення множення на додавання

Передана функція G(s) зазвичай є добутком і відношенням простих факторів — полюсів, нулів, сталої коефіцієнта. Помноження частотних характеристик разом є незручним; додавання їх – легке. Інсайт Генріка Боде в 1930-х роках у Bell Labs полягав у тому, щоб відобразити все на логарифмічній шкалі, оскільки log(|G₁ · G₂|) = log|G₁| + log|G₂|. Візьміть величину в децибелах, |G(jω)|dB = 20·log₁₀|G(jω)|, побудуйте її проти log ω, і відповідь складної каскадної системи полюсів та нулів стає сумою простих прямолінійних внесків, які можна намалювати від руки.

Що вносить кожен фактор

Кожен першого порядку у G(s) додає один передбачуваний асимптотичний нахил та один передбачувальний фазовий зсув, що центруються навколо власної кутової частоти ωc = 1/τ:

полюс 1/(1 + s/ωc) → амплітуда рівномірна, потім падає на -20 дБ/декаду після ωc фаза 0° → -90°, центрована (-45°) при ω = ωc нуль (1 + s/ωc) → амплітуда рівномірна, потім зростає на +20 дБ/декаду після ωc фаза 0° → +90°, центрована (+45°) при ω = ωc інтегратор 1/s → амплітуда падає на -20 дБ/декаду на всіх частотах фаза фіксована при -90° всюди коефіцієнт K → зміщує весь графік амплітуди вгору/вниз, без впливу на фазу Пара комплексних-сполучених полюсів (другої черги знеясної системи) додає -40 дБ/декаду після власної природної частоти ωn, з фазовим зсувом від 0° до -180°; якщо відносна в'язкість ζ невелика, графік амплітуди різко пікове перед тим, як він починає спадати – резонантний горбок, який ви бачите у будь-якій легко амортизованій механічній або електричній системі. Намалюйте прямокутну асимптотику кожного фактора, додайте їх усі разом, і ви отримаєте наближення Боде — точне скрізь, крім протягом приблизно десятиліття навколо кожної кутової частоти, де справжня крива округлюється приблизно на 3 дБ у кутовому полю/нулі зі змінною ознакою замість того, щоб утворювати різкий кут.

pole  1/(1 + s/ωc)  → magnitude flat, then falls at -20 dB/decade past ωc
                        phase 0° → -90°, centred (-45°) at ω = ωc
zero  (1 + s/ωc)     → magnitude flat, then rises at +20 dB/decade past ωc
                        phase 0° → +90°, centred (+45°) at ω = ωc
integrator  1/s      → magnitude falls at -20 dB/decade at ALL frequencies
                        phase pinned at -90° everywhere
gain  K              → shifts the whole magnitude curve up/down, no phase effect

Визначення стабільності за допомогою діаграми

Два частоти мають більше значення, ніж інші точки на кривій. Коефіцієнт підсилення перехресна частота ωgc — це точка, де |G(jω)|dB = 0 дБ — петльове підсилення становить точно 1. Частота перетину фази ωpc — це точка, де фаза досягає -180° — точку, де негативний зворотний зв’язок перемикається на позитивний. З цих двох точок ви читаєте два класичні маржі стабільності безпосередньо з осей:

маржа підсилення GM = 0 - |G(j·ω_pc)|dB (дБ додаткового підсилення перед досягненням 0 дБ фазою -180°) маржа фази PM = 180° + phase(G(j·ω_gc)) (градусів додаткової затримки фази перед нестабільністю) Обидва маржі відповідають одному основному питанню — наскільки петля віддалена від задоволення умови нестабільності Нейкви |G(jω)| = 1 при фазі -180°? — з двох різних напрямків. Загальне правило, що використовується в галузі управління: маржа фази вище приблизно 45° (і маржа підсилення близько 6 дБ) забезпечує відгук на кроки з помірним, добре амортизованим осцилюванням; зменшення PM до 0° призводить до того, що закрита петля все важче і важче гуде, аж поки не починає безперервно коливатися на межі стабільності.

gain margin  GM  = 0 - |G(j·ω_pc)|_dB      (dB of extra gain before -180° hits 0 dB)
phase margin PM  = 180° + phase(G(j·ω_gc))  (degrees of extra phase lag before instability)

Чому крива фаз є найнесподіванішою

Асимптоти величини інтуїтивні — більше полюсів, швидший спад. Саме на кривій фаз приховується немінімальна поведінка: правий-півплоскосний нуль або чисто запізнення в часі додають фазове зсув без зміни кривої величини. Затримка e-jωT має |e-jωT| = 1 для кожної частоти — непомітна на графіку величини, але її фаза становить -ωT радіан, що росте без обмежень зі збільшенням частоти. Саме тому системи керування з значною транспортною затримкою (хімічний процес, робот, керований мережею) можуть виглядати повністю стабільними на графіку величини та при цьому мати майже нульовий запас фази, коли додають фазового зсуву цієї затримки.

Bode проти Nyquist проти root locus

Усі три інструменти читають одну й ту ж характеристичну рівняння, 1 + G(s)H(s) = 0, з різних точок зору. Графік Найквиста відстежує G(jω) у комплексній площині та підраховує обчислення навколо -1, що є єдиним із трьох, який залишається строгим для систем з відкритим правою півплощиною полюсів. Root locus відстежує розташування закритих петльових полюсів безпосередньо, коли змінюється коефіцієнт зворотного зв’язку K, що графік Боде не може показати — графік Боде описує лише відкриту частотну характеристику передавальної функції. Перевага графіка Боде полягає в тому, що він найпростіший із трьох для малювання від передавальної функції вручну, найлегше вимірюється експериментально (переміщайте синусоїдальний сигнал у реальну систему та записуйте амплітуду та фазу, не потрібен жоден модель) і той, який найбільш безпосередньо відображає щоденну мову дизайну – коефіцієнт прибутку, межа фази та смуга пропускання.

Часті запитання

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Bode Plot і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.