Bode/Nyquisti juhtimislabor

liiguta poole liuguritega ja jälgi, kuidas robustsusvarud koheselt muutuvad

Algväärtused — Bode/Nyquisti juhtimislabor

Ahela võimendus K	1
Omavõnkesagedus wn	3
Summutustegur zeta	0.7

Muuda ülal mõnda juhtelementi ja tulemused arvutatakse uuesti.

Interaktiivse simulatsiooni laadimine...

Õppetund

Teooria — Bode/Nyquisti juhtimislabor

Nyquisti stabiilsuskriteerium otsustab tagasisidestatud ahela stabiilsuse loendamise, mitte lahendamise teel. Kanna avatud ahela sagedusvastus G(jω) komplekstasandile, loenda, mitu korda kõver punkti −1 + j0 ümbritseb, ja liida juurde ebastabiilsete pooluste arv, mis ahelal oli juba enne sulgemist: Z = N + P. Iga nullist erinev Z tähendab, et suletud ahel jookseb ära.

Mida iga sümbol tähendab

P: avatud ahela G(s) nende pooluste arv, mis asuvad juba paremal pooltasandil — loendatuna enne ahela sulgemist. Selle lehe objekti puhul on see alati 0.
N: mitu korda Nyquisti kõver punkti −1 + j0 päripäeva ümbritseb. Vastupäeva tiirud loevad negatiivselt, seega võib N jääda nullist allapoole.
Z: suletud ahela pooluste arv paremal pooltasandil — see vastus, mida tegelikult otsiti. Z = 0 on stabiilsus; kõik muu ei ole.
−1 + j0: kriitiline punkt, Nyquisti paneelil punaselt märgitud. Seal on ahela võimendus täpselt 1 ja faasihilistus täpselt 180°: korra ümber ahela saadetud signaal naaseb sama suurena ja tagurpidi.

Kust valem tuleb

Küsi õigesti. Negatiivse ühiktagasiside korral on suletud ahela ülekandefunktsioon G/(1 + G), seega on selle poolused 1 + G(s) nullkohad. Stabiilsus taandub ühele küsimusele: kas 1 + G(s) omab nullkohti paremal pooltasandil? Nende leidmine tähendaks polünoomi tegurdamist, mida sul ehk polegi.
Loenda neid, ära otsi. Argumendiprintsiip ütleb: kui s läbib korra päripäeva suletud kontuuri, tiirleb 1 + G(s) kujutis ümber nullpunkti täpselt nii mitu korda, kui palju on sees nullkohti, miinus sees olevad poolused. Nullpunkti ümber tiirlemine funktsioonil 1 + G on sama liikumine kui −1 ümber tiirlemine funktsioonil G — sellepärast ongi kriitiline punkt just seal, kus ta on, ja sellepärast on ta −1, mitte mõni muu arv.
Kontuuriks võta kogu parem pooltasand: mööda imaginaartelge üles ja lõpmatu poolringiga tagasi. Nüüd tähendavad „kontuuri sees“ ja „ebastabiilne“ sama asja. Sissejäävad nullkohad on Z, sissejäävad poolused P, tiirude arv on N, ja argumendiprintsiip kõlab N = Z − P — ümber tõstetuna Z = N + P. Imaginaarteljel on s = jω ja G(jω) ongi sagedusvastus: seesama kõver, mida see leht niigi joonistab.
Just see viimane lause ongi kogu võit. Loendamiseks pole mudelit vaja. G(jω) on mõõdetav — raputa päris riistvara siinuslainetega ja salvesta, mis tagasi tuleb —, nii saad stabiilsusotsuse süsteemi kohta, millele keegi pole kunagi võrrandit kirjutanud. Võimendus- ja faasivaru on siis kaks viisi küsida, kui kaugelt mõõdetud kõver punktist −1 mööda läheb: üks mööda negatiivset reaaltelge, teine mööda ühikringi.

Kuidas nähtut lugeda

Kõik kolm paneeli näitavad avatud ahelat — G(jω) üksinda, tagasisideharu veel ühendamata —, samal ajal kui kaks varu nende all kirjeldavad seda, mis juhtub pärast ühendamist. Amplituudi tipp ω_n lähedal ning ka ω_r ja M_r tabelis kuuluvad seega funktsioonile G, mitte sellele suletud ahelale, mida sa projekteerid. Nyquisti paneelil on pidev joon ω nullist ülespoole ja katkendjoon selle peegelpilt negatiivse ω jaoks. Vaja on mõlemat: ümbritsemine tekib alles siis, kui kontuur sulgub, ja sagedusskaneering suudab mõõta ainult poole sellest. Nooled näitavad läbimissuunda, mis otsustabki, kas tiir loeb päripäevaseks.

Eeldab: Negatiivne ühiktagasiside ja teadaolev P. See leht ei küsi P-d kunagi — ta võtab selle nulliks, mis on siin õige ainult sellepärast, et sumbuvusliugur peatub väärtusel ζ = 0.05 ja mõlemad poolused jäävad igal saavutataval seadistusel vasakule pooltasandile. Anna sama kriteerium objektile, mis on avatud ahelas juba ebastabiilne, ja ümbritsemisteta kõver tähendab täpselt vastupidist sellele, mida ta siin tähendab. Eeldatakse ka, et teel ei taandunud midagi: nulliga taandatud parema pooltasandi poolus kaob G(s)-st, kuid jääb riistvarasse alles — ja loendus annab siis stabiilsustunnistuse ahelale, mis ikkagi ära jookseb.
Ei kehti, kui: Nyquist vastab jah või ei ühe täpse mudeli kohta, aga riistvara pole ei täpne ega paigal. Z = 0 ütleb, et nominaalne ahel on stabiilne, ja mitte kõige vähematki selle kohta, kui palju peab võimendus triivima või temperatuur nihkuma, enne kui ta enam ei ole — just seepärast trükib see leht kaks varu, mitte ainsat sõna „stabiilne“. Varud ise on liiga optimistlikud, sest kumbagi mõõdetakse teist paigal hoides: ahel võib näidata mugavat võimendusvaru ja mugavat faasivaru ning ometi viib mõõdukas muutus mõlemas korraga kõvera otse läbi punkti −1. Aus arv on kõvera lühim kaugus kriitilisest punktist ükskõik millises suunas, ja seda ei joonista siin ükski paneel.

Seda regulaatorringi ei saa võimendusega ebastabiilseks muuta 🖖

Kumbki kahest poolusest annab kõige rohkem 90° faasimahajäämust ja jõuab oma täie 90°-ni alles lõpmatul sagedusel, nii et K·ω_n²/(s² + 2ζω_ns + ω_n²) faas läheneb −180°-le, kuid ei jõua sinna kunagi. Lükka kõik kolm liugurit piirini ja see peatub ikka −179,99° juures. Just seepärast näitab võimendusvaru ∞: Nyquisti kõver ei mähi end kunagi ümber −1 ja suurem K üksnes venitab seda. Ebastabiilsus vajab siin kolmandat mahajäämuse allikat – lisapoolust, integraatorit, transpordiviivet – mitte suuremat võimendust.

Võnkumise ennustamine võrrandeid lahendamata 🖖

Diferentsiaalvõrrandi lahendamise asemel ergutatakse süsteemi paljudel sagedustel siinussignaalidega ja jälgitakse vastust. Bode diagramm näitab, kui palju iga sagedus võimendub ja kui palju selle faas hilineb. Oht tekib seal, kus võimendus läbib 0 dB, samal ajal kui faas läheneb −180°-le: tagasiside hakkab end ise võimendama ja ahel võngub. Faasivaru mõõdab, kui kaugel sa sellest piirist oled — suurem varu tähendab rahulikumat juhtimist.

Faasivaru peidab viivituse eelarve 🖖

Terve faasivaru näitab ka, kui suurt puhast ajaviivitust ahel talub. Viivitus lisab faasinihke, mis kasvab koos sagedusega (−ωτ), ja kulutab varu kõige kiiremini seal, kus ahel on kõige aktiivsem. Lubatud viivitus on ligikaudu τ = PM / ω_c, kus PM on radiaanides ja ω_c on võimenduse ületussagedus. Nii võib kiire ja laia ribalaiusega ahel puruneda mõne millisekundi anduri viivituse tõttu, mida sama varuga aeglasem ahel ei märkagi.

Levinud vale arusaam

Suurem võimendus ei tähenda alati paremat juhtimist. See võib parandada kiirust, hävitades samal ajal vaikselt robustsusvarud.

Ülesanne täielikult lahendatud

Resonantstipp 1,0002 ja lõpmatu võimendusvaru 6 sammu
Teist järku kontuur, mille ζ = 0,7. Paneel kuvab resonantsitipuks 1,0002 ja võimenduse varuks lõpmatuse. Selgitage mõlemat.
1. Standardne teist järku kuju koos kolme parameetriga, mida liugurid määravad. Kõik alljärgnev loetakse välja nimetajast.
  
  $G (s) = \frac{K ω _{n}^{2}}{s ^{2} + 2 ζ ω _{n} s + ω _{n}^{2}} = \frac{9}{s ^{2} + 4.2 s + 9} (K = 1, ω_{n} = 3, ζ = 0.7)$
2. Poolused on juured. Reaalosa on sumbumiskiirus ja kujutosa on võnkesagedus ning mõlemad on nähtavad astmevastes.
  
  $s = - ζ ω_{n} \pm j ω_{n} 1 - ζ^{2} = - 2.1 \pm j 3 0.51 = - 2.1 \pm j 2.142$
3. Resonantsisagedus ei ole omasagedus. Selle all olev ruutjuur võib muutuda nulliks, mis ongi järgmise sammu sisu.
  
  $ω_{r} = ω_{n} 1 - 2 ζ^{2} = 3 1 - 0.98 = 3 0.02 = 0.424 rad/s$
4. See muutub nulliks väärtusel ζ = 1/√2. Väärtusel 0,7 on süsteem vaevalt seespool tippudega piirkonda ning tipp on 0,002 dB — algebraliselt olemas, kuid ühelgi graafikul nähtamatu.
  
  $1 - 2 ζ^{2} > 0 ⟺ ζ < \frac{1}{2} = 0.7071; M_{r} = \frac{1}{2 ζ 1 - ζ ^{2}} = 1.0002$
5. Nüüd faas. Iga poolus annab faasi mahajäämuse, mis läheneb 90°-le, seega kaks poolust lähenevad 180°-le ülalt ega ületa seda kunagi.
  
  $∠ G (j ω) = - arctan \frac{2 ζ ω _{n} ω}{ω _{n}^{2} - ω ^{2}} ⟶ - 18 0^{\circ} only as ω \to \infty$
6. Ületamise puudumine tähendab, et tegelikku võimenduse varu ei ole olemas, ning just seda paneel kuvabki.
  
  $PM = 163.7 5^{\circ}, GM = \infty$
Vastus

ζ = 0,7 asub vahetult allpool väärtust 1/√2 ning kaks poolust ei saa jõuda väärtuseni −180°. Resonantsisageduse avaldis sisaldab tegurit √(1 − 2ζ²), mis muutub nulliks väärtusel ζ = 0,7071 — seega on 0,7 vaid 0,001 võrra seespool piirkonda, kus tipp üldse eksisteerib, ning selle tipp on 1,0002 ehk 0,002 dB. See künnis on Butterworthi tingimus ja seepärast on 0,707 sumbuvus, mida iga filtriõpik soovitab: see on kõige kiirema vastega seadistus, millel veel tipp puudub. Lõpmatu võimenduse varu on tugevam ja ohtlikum väide. Iga poolus annab maksimaalselt 90° faasi mahajäämust, seega kaks poolust lähenevad 180°-le ega jõua kunagi selleni ning ükski võimendus ei pane seda kontuuri ostsilleeruma. Reaalsed kontuurid hakkavad varem või hiljem ikkagi ostsilleeruma — kolmas poolus, transpordiviivitus, modelleerimata resonants —, mistõttu lõpmatu võimenduse varu on kindel märk sellest, et mudelist on midagi puudu, mitte sellest, et riistvara on ohutu.

$ζ = 0.7 sits just under \frac{1}{2}; GM = \infty$

Allikad (2)

The lesson — the encirclement count Z = N + P, in the paper that introduced it: H. Nyquist, "Regeneration Theory." Bell System Technical Journal 11(1), 126–147, 1932.
Insight block 3 — phase margin as a delay budget: K. Ogata, Modern Control Engineering, 5th edition, ch. 7. Prentice Hall, 2010. ISBN 978-0-13-615673-4.

42Mathi meeskond · Uuendatud 2026-07-29

Kontrollitud avaldatud allikate põhjal. Parandused on teretulnud. Toimetamispõhimõtted · Parandused

Kõik Füüsika õppetunnid →

Näiteülesanded

stabiilne baasjuht - Kui K = 1, on alalisvooluvõimendus ühikuline, lõikesagedus 0,60 rad/s ja faasivaru 163,74°. Stabiilsuspiir on siit väga kaugel.
peaaegu võnkuv - Tööriista vähim stabiilne eelseadistus: faasivaru on 16,70° ja resonantstipp 16,44 dB, mistõttu lubatav viide on vaid 52 ms.
aeglane, kuid robustne - K = 0,8 ja ζ = 0,9, mistõttu võimendus ei jõua üldse 0 dB-ni. Lõikesagedust pole, faasivaru on lõpmatu ning ükski viide ei saa süsteemi destabiliseerida.
kõrge võimenduse risk - Nimi viitab ohule, kuid kahest väikese stabiilsusvaruga eelseadistusest on see turvalisem: faasivaru on 21,66°, eelseadistusel „Võnkumise piiril“ aga 16,70°. Põhjus on selles, et ζ = 0,35 võrreldes 0,2-ga korvab kuhjaga K = 4,2 võrreldes 2,6-ga.