Wykres Boda: Interpretacja Stabilności na Podstawie Przeszukiwania Częstotliwości

Przekształcanie mnożenia w dodawanie

Transfer function G(s) jest zwykle iloczynem i stosunkiem prostych czynników — biegunów, pierwiastków zerowych oraz stałej wartości. Mnożenie odpowiedzi na częstotliwość przez siebie jest uciążliwe; dodawanie ich jest łatwe. Wskutek Hendrika Bode’a z lat 30-tych w Bell Labs polegał on na rysowaniu wszystkiego w skali logarytmicznej, ponieważ log(|G₁ · G₂|) = log|G₁| + log|G₂|. Weź pod uwagę wartość bezwzględną w decybelach, |G(jω)|dB = 20·log₁₀|G(jω)|, narysuj ją przeciwko log ω, a odpowiedź skomplikowanej kaskady biegunów i pierwiastków zerowych staje się sumą prostych liniowych kontrybuant, które można szkicować ręcznie.

What each factor contributes

Every first-order factor in G(s) adds one predictable asymptotic slope and one predictable phase swing, centred on its own corner frequency ωc = 1/τ:

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 A pair of complex-conjugate poles (a second-order underdamped factor) contributes -40 dB/decade past its natural frequency ωn, with a phase swing of 0° to -180°; if the damping ratio ζ is small, the magnitude curve peaks sharply just before it rolls off — the resonant bump you see on any lightly-damped mechanical or electrical system. Sketch each factor's straight-line asymptote, add them all up, and you get the Bode approximation — correct everywhere except within about a decade of each corner frequency, where the true curve rounds off by roughly 3 dB at a real pole/zero corner instead of forming a sharp angle.

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

Odczytanie stabilności bezpośrednio z wykresu

Ważniejsze są częstotliwości niż reszta krzywej. Częstotliwość przecięcia zysku ωgc, gdzie |G(jω)|dB = 0 dB — czyli wzmocnienie wynosi dokładnie 1. Częstotliwość przecięcia fazy ωpc, gdzie faza osiąga -180° — punkt, w którym sprzężenie zwrotne staje się dodatnie. Z tych dwóch punktów odczytuje się bezpośrednio dwa klasyczne marginesy stabilności osi po osi:

margines_zysku GM = 0 - |G(j·ω_pc)|dB (dźwięki dodatkowego wzmocnienia przed uderzeniem 0 dB w -180°) margines_fazowy PM = 180° + faza(G(j·ω_gc)) (stopnie dodatkowego opóźnienia fazowego przed niestabilnością) Oba marginesy odpowiadają temu samemu pytaniu — jak daleko jest ta pętla od spełnienia warunku niestabilności Nyquista |G(jω)| = 1 w fazie -180°? — z dwóch różnych kierunków. Zasada ogólna stosowana w inżynierii sterowania: margines fazowy powyżej około 45° (i margines zysku powyżej około 6 dB) daje odpowiedź z minimalnym, dobrze tłumionym przejściem; obniżenie PM do około 0° powoduje, że pętla zamykana w coraz większym stopniu kołuje aż do ciągłej oscylacji na granicy niestabilności.

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)

Dlaczego krzywa fazowa jest tą, która zaskakuje ludzi

Asympty magnituudy są intuicyjne — im więcej biegunów, tym szybsze opadanie. Faza to miejsce, gdzie zachowanie poza minimalną fazą ukrywa się: zerowy kąt w płaszczyźnie połówkowej lub czysta opóźnienie czasowe dodaje opóźnienie fazowe bez zmiany krzywej magnituudy. Opóźnienie e-sT ma |e-jωT| = 1 dla każdej częstotliwości — niewidoczne na wykresie magnituudy, ale jego faza wynosi -ωT radianów, rosnąc nieograniczenie wraz ze wzrostem częstotliwości. To właśnie dlatego systemy sterowania z istotnym opóźnieniem transportowym (proces chemiczny, robot sterowany siecią) mogą wyglądać idealnie stabilnie na wykresie magnituudy i wciąż mieć prawie brak marginesu fazowego po dodaniu faza opóźnienia w czasie.

Bode Plot: Częstotliwość i Margines Stabilności

Wszystkie trzy narzędzia odczytują tę samą równanie charakterystyczne, 1 + G(s)H(s) = 0, z różnych perspektyw. Wykres Nyquista śledzi G(jω) w płaszczyźnie zespolonej i liczy okrążenia wokół -1, która jest jedynym z trzech, które pozostają rygorystyczne dla systemów z polami prawego-górnego pola otwartego. Locus korzeni śledzi lokalizację pól zamkniętych bezpośrednio w miarę zmiany wygranej K, czego Bode nigdy nie może pokazać — Bode opisuje jedynie odpowiedź częstotliwościową otwartą. Zaletą Bode jest to, że jest najłatwiejszy z trzech do narysowania ręcznie z funkcji transferowej, najłatwiej się mierzy eksperymentalnie (przesuń sinusiwdę do rzeczywistego układu i zarejestruj amplitudę i faza, nie wymagany model) oraz ten, który w najbardziej bezpośredni sposób mapuje się na codzienny język projektowy z marginesem zysku, marginesem fazowym i szerokością pasma.

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Bode Plot i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.