Begrijpen van de Nyquist Plot

Het Nyquist-complot is een fundamenteel hulpmiddel in de controletechniek, dat gebruikt wordt om de stabiliteit van feedbacksystemen te analyseren. Genoemd naar Harry Nyquist, geeft deze grafische voorstelling de frequentierespons weer van een systeem door de complexe overdrachtsfunctie G(jω) als de hoekfrequentie ω varieert van 0 tot ∞. Het plot wordt getekend in het complexe vlak, met het werkelijke deel van G(jω)[] op de horizontale as en het denkbeeldige deel op de verticale as. Het primaire nut ervan is het beoordelen van de gesloten-lusstabiliteit van de open-loopfrequentierespons, zonder expliciet te rekenen gesloten-looppolen. Het Nyquist-plot biedt ook inzicht in winst en fasemarges, waardoor het onmisbaar is voor robuuste controllerontwerp.

De wiskundige onderbouwingen van Nyquist-compatibele percelen vereisen een grondig begrip van complexe analyse, met name de evaluatie van overdrachtsfuncties op de denkbeeldige as en de geometrische interpretatie van contouren. Ingenieurs moeten het concept van de Nyquist-contour[], een gesloten pad in het complexe vlak dat het gehele rechter halve vlak omsluit (RHP) begrijpen en hoe deze contour door de open-loop-transferfunctie G(s)H(s)[] het Nyquist-pad[] in kaart brengt. Het resulterende plot draait rond het kritische punt −1 + 0[j geeft direct de stabiliteit van het systeem aan, volgens het Nyquist-stabiliteitscriteriumcriterium. Dit artikel breidt uit op deze wiskundige grondslagen, wandelend door de constructielogica, het principe van argument en praktische valpartijen.

Complexe beoordeling van de overdrachtsfunctie

De open-loop overdrachtsfunctie G(s)H(s)[] is een rationele functie van de complexe variabele s[ = σ + j[ω. Voor frequentieresponsanalyse stellen we σ = 0, dus s[] = []j[ω. De resulterende expressie G(jω)H(jω)[ is een complex getal dat afhankelijk is van ω:

G(jω)H(jω) = Re(ω) + j Im(ω) = M(ω) e^{jφ(ω)}

Hier is M(ω) de omvang (ook wel gain genoemd) en φ(ω) de fasehoek. Het bouwen van het Nyquist-complot omvat het evalueren G(jω)H(jω) op vele frequenties en het plotten van de complexe punten. Voor een typisch systeem begint het perceel bij ω = 0 (vaak op een punt op de werkelijke as) en eindigt het bij ω → ∞ (meestal bij de oorsprong, tenzij het systeem zuivere integrators heeft). De vorm van de curve toont belangrijke dynamische eigenschappen: bijvoorbeeld hoe de fase afneemt met toenemende frequentie, of hoe de omvang zich af rolt.

Om G(jω)H(jω) wiskundig te berekenen, moet je een algemene overdrachtsfunctie in polen en nullen zien:

G(s)H(s) = K \frac{\prod {i=1}^{m} (s - z i)}{\prod {k=1}^{n} (s - p k)}

Wanneer s = j[ω, wordt elke term een complexe vector van de nul- of poollocatie tot het punt j[ω op de denkbeeldige as. De totale omvang is het product van de lengtes van de nulvectoren gedeeld door het product van de lengten van de poolvectoren, vermenigvuldigd met

Wiskundige bouwstappen

Het construeren van een Nyquist-complot met de hand (of het gebruiken van software) impliceert een systematische evaluatie van G(jω)H(jω) over een frequentiebereik. De belangrijkste stappen zijn:

  1. Identiteit de overdrachtsfunctie en schrijf het in Bode-vorm (polen en nullen). Voor systemen met vertraging, omvatten de e^{-jωT} factor.
  2. Bereken de omvang en fase met een reeks frequenties. Meestal worden lage frequenties (bij 0) en hoge frequenties (bij ∞) gekozen, met tussenpunten bij breukfrequenties (waar pool- of nulbijdragetekens veranderen).
  3. Start bij ω = 0: Het punt hangt af van het systeemtype. Voor Type 0 (geen integrators), begint het perceel bij K (de DC-winst) op de positieve reële as. Voor Type 1 (één integrator) begint het perceel op een oneindige grootte met een fase van −90°. Voor Type 2 (twee integrators) is de fase −180° bij ω = 0.
  4. Plotpunten bij toenemende frequenties, die ze soepel verbinden. Naarmate de frequentie toeneemt, beweegt de curve zich over het algemeen naar de oorsprong (sinds de afname van de magnitude). De fase kan −180° kruisen bij de fase-crossoverfrequentie waar de grootte cruciaal is voor stabiliteitsanalyse.
  5. Spiegel voor negatieve frequenties (als de volledige Nyquist-contour gebruikt wordt). Het plot voor ω van −∞ tot 0 is het complexe geconjugeerde plot voor positieve ω, dat over de werkelijke as wordt weerspiegeld. In combinatie met de positieve frequentietak (en de oneindige halve cirkel als het systeem netto pool-nul overtollig is), wordt het volledige Nyquist-pad verkregen.
  6. Identificeer het aantal omcirkelingen van −1 + 0j om het stabiliteitscriterium toe te passen.

Voor systemen met relatieve graad (n − m) ≥ 1, eindigt de curve bij de oorsprong, nadert deze met een fase van −(n−m)×90°. Voor relatieve graad 0 (hetzelfde aantal polen en nullen), is het laatste punt een eindig niet-nul complex getal bij ω → ∞.

Het criterium van de stabiliteit van Nyquist in detail

Het Nyquist stabiliteitscriterium is een krachtige methode om de stabiliteit van de gesloten lus te bepalen vanuit de open lusfrequentierespons. Dit beginsel bepaalt dat als een gesloten contour Γ in het -vlak -vlak wordt toegewezen door middel van een functie F(s)] naar een gesloten contour ΓF in het F-vlak, dan het aantal keren Γ omcirkelt de oorsprong gelijk is ([]] [P[[]]]]]]], waar Z]] is het getal van nul [LT:[LT:[LT]][LT][LT][LT

Omdat F(s) gerelateerd is aan de functie van de open lusoverdracht door F(s) = 1 + G(s)H(s), is de kaart van F(s)[] door de Nyquist contour eenvoudigweg het Nyquist plot van G(s)H(s)[]] naar rechts verschoven door 1. Equivalente, encircllements van de oorsprong in F-vlak precies overeenkomt met omcirkels van het punt −1 + 0[j] in de [G(s)H(s)[[-vlak.]-vlak. Het Nyquist-criterium luidt dan:

N = Z − P

waarbij:

  • N = aantal omcirkelingen van −1 + 0j door het Nyquist-complot (in een richting met de klok mee als positief, tegen de klok in als negatief).
  • Z = aantal gesloten poolstokken in het RHP (nul van 1 + G(s)H(s) in het RHP).
  • P = aantal open-luspalen in het RHP.

Voor een stabiel gesloten-lussysteem is Z = 0. Zo is het gesloten-lussysteem stabiel als en alleen als het Nyquist plot van G(s)H(s)[] het punt omcirkelt −1 + 0j precies ]−P[ keer (d.w.z. tegen de klok in) encirkels gelijk ]P[[ als we kloksgewijs als positief tellen). Als [P[[ = 0 (open-loop stabiel), dan moet het Nyquist plot [[not[[ omcirkel −1] bij alle gesloten-loopstabiliteit.

Toepassing van het criterium: tellen omcirkelen

Het tellen van omcirkelingen houdt in dat het Nyquist pad over de volledige contour wordt geanalyseerd, inclusief de in kaart brengen van de oneindige halve cirkel. Voor systemen met relatieve graad ≥ 1, de oneindige halve cirkel kaarten op een kleine cirkel rond de oorsprong in de G(s)H(s)] vlak (of op een vast punt als relatieve graad nul). De praktische benadering voor ingenieurs is om:

  1. Teken het Nyquist-plan voor ω van 0 naar ∞.
  2. Reflecteer het over de echte as voor ω van 0 tot −∞.
  3. Als het systeem polen op de denkbeeldige as heeft (bijvoorbeeld integrators), dan inspringen we de Nyquist-contour met kleine halve cirkels naar rechts om deze polen te vermijden, en in kaart brengen we die inspringingen naar grote boog in het G(s)H(s) vlak.
  4. Tel het netto aantal kloksgewijze omcirkelingen van −1 door de volledige gesloten curve (inclusief de oneindige halve cirkelmapping) te traceren.

Wiskundig kan het aantal omcirkelingen N worden berekend via het gehele liquidatienummer:

N = \frac{1}{2πj} \oint {Nyquist \, contour} \frac{G'(s)H's)}{1+G(s)H(s)} \, ds

Maar dit wordt zelden analytisch geëvalueerd; ingenieurs vertrouwen op geometrische inspectie of softwaretools. Het belangrijkste inzicht is dat omcirkelingen optreden wanneer de fase van G(jω)H(jω) kruist −180° (d.w.z. het punt −1 op de werkelijke as) terwijl de grootte ≥ 1 is. Deze voorwaarde definieert de gain margin[ en ]fasemarge[.

Praktische bouwvoorbeelden

Om begrip te versterken, beschouwen we een eenvoudig systeem van eerste orde en een systeem van tweede orde met demping.

Eerste-ordesysteem: G(s) = 1/(s+1)

Hier G(jω) = 1/(jω+1). De magnitude is M(ω) = 1/√(1+ω2)] en fase is φ(ω) = −tan−1(ω). Bij ω = 0 is het punt (1, 0). ω → ∞ komt het punt uit een hoek van −90°. De Nyquist curve is een halve cirkel in het onderste halve vlak (sinds fase is altijd negatief) die zich uitstrekt van de positieve reële as tot de oorsprong. Deze curve doet not encircle −1, en sinds P = 0, het gesloten-loopsysteem is stabiel voor elke positieve winst K[F] als we de AUC vervangen] vervangen [FLT:[LT:[LT]]]] [FLT:

Tweedeordesysteem met Damping: G(s) = ω n2/(s2 + 2ζω n s + ω n2)

Voor een tweede ordesysteem met een natuurlijke frequentie is dit een "Borloge" n] en de demping ratio ζ, de "Borloge" van de "Nyquist" van de "Borloge" van de "Borloge" van de "Borloge" van de "Borloge" van de "Borloge" van de "Follow" van de "FLT:5]" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van de "F" van

Om een juiste fase-overstap te zien, moet je een systeem met een extra pool op lage frequentie bekijken: G(s) = 1/(s+1)(s2 + 2ζ s + 1)]. Hier kan de fase −180° kruisen bij een eindige ω, en het plot kan omcirkelen −1 als de winst voldoende hoog is.

Nyquist, Bode en Nichols Plots verbinden

Het Nyquist-diagram is nauw verwant met het Bode-diagram (materiële en fase) en de Nichols-diagram. Het Bode-diagram geeft de omvang in dB en fase in graden als functies van logfrequentie. Het Nyquist-diagram is een parametrische combinatie van de magnitude en fase: elk punt op de Nyquist-kromme komt overeen met een magnitude en fase bij een gegeven ω. De gain-marge[] is de factor waarmee de winst kan worden verhoogd voordat de omvang bij de fase-crosssoverfrequentie 1 (d.w.z. de afstand van −1 tot de Nyquist-curve langs de negatieve reële as) wordt. De fasemarge[] is de extra fasevertraging die vereist is bij de winstcrossoverfrequentie (waar de omvang = 1) om de fase te maken −180°. Deze marges zijn direct te lezen uit de Nyquist-plot: de winstmarge is 1/G(jω]pc[

De Nichols grafiek transformeert het Nyquist plot in een Cartesiaanse plot van magnitude (dB) vs. fase (deg). Deze transformatie is nuttig voor het ontwerpen van controllers omdat aanwinst en fasemarges gemakkelijk te zien zijn, en gesloten-lus response kan worden gelezen met behulp van M-cirkels. Alle drie de voorstellingen zijn wiskundig gelijkwaardig; het kiezen van een is afhankelijk van de context en de voorkeur van de ingenieur.

Beperkingen en geavanceerde overwegingen

Hoewel krachtig, heeft het Nyquist-complot beperkingen die ingenieurs moeten erkennen. Ten eerste, voor systemen met vertraging, wordt de Nyquist-curve oneindig en sterk gebogen, waarbij de oneindige halve cirkel mapping zorgvuldig moet worden behandeld. De vertraging voegt een faseverschuiving toe −ωT die lineair accumuleert, waardoor het Nyquist-complot naar binnen of naar buiten wentelt afhankelijk van de vertraging. Stabiliteitsanalyse vereist dan het aangepaste Nyquist-criterium dat de vertraging verklaart met het concept van de Nyquist-contour dat de vertragingspool bij oneindigheid uitsluit.

Ten tweede, het Nyquist-plank zoals beschreven is voor continu-tijdsystemen. Voor discrete-tijdsystemen (z-domein) wordt de Nyquist-contour de eenheidscirkel, en het stabiliteitscriterium onderzoekt omcirkelingen van −1 in het w-vlak of het gebruik van de bilineaire transformatie om de kaart naar het s-vlak. De wiskundige funderingen blijven gelijk, maar de contour verandert.

Ten derde gaat het Nyquist-criterium ervan uit dat het gesloten-lussysteem stabiel is als en alleen als het open-lus Nyquist-complot op passende wijze omcirkeld wordt −1. Dit vereist echter dat het karakter van de open-lus-polen in het RHP bekend is a priori. Als de open-lus-transferfunctie polen op de denkbeeldige as heeft (inclusief zuivere integrators), moet de Nyquist-contour worden ingespringd om deze polen te vermijden, en de inkarping van die inspringingen levert oneindige boogjes die in het encirclementsaantal moeten worden opgenomen. Dit kan lastig zijn voor ingenieurs die nieuw zijn in de methode.

Ten slotte is het Nyquist-complot een frequentiedomeininstrument; het biedt noodzakelijke en voldoende stabiliteitsvoorwaarden alleen voor lineaire tijd-invariante (LTI) systemen. Voor niet-lineaire of tijd-variable systemen is het Nyquist-criterium niet rechtstreeks van toepassing, hoewel uitbreidingen zoals het circle criterium bestaan voor bepaalde klassen van niet-lineairheden.

Conclusie

De constructie van een Nyquist-complot is diep geworteld in complexe analyse, met name de mapping van de Nyquist-contour door de open-loop-overdrachtfunctie en het principe van het argument. Door de wiskundige stappen te beheersen die de evaluatie van G(jω)H(jω)[], het inlassen van de complexe punten en het tellen van omcirkels van −1

Voor verdere lezing, zie standaard leerboeken zoals Wikipedia's artikel over het Nyquist stabiliteitscriterium[, Control Tutorials for MATLAB and Simulink, of de klassieke tekst Modern Control Engineering door Ogata. Deze bronnen bieden extra voorbeelden en diepere wiskundige deducties.