Das Verständnis der Nyquist Plot

Der Nyquist-Plot ist ein grundlegendes Werkzeug in der Steuerungstechnik, das zur Analyse der Stabilität von Rückkopplungssystemen verwendet wird. Diese grafische Darstellung zeigt den Frequenzgang eines Systems, indem die komplexe Übertragungsfunktion FLT:0 G(jω) aufgetragen wird, da die Winkelfrequenz ω von 0 bis ∞ variiert. Der Plot wird in der komplexen Ebene mit dem realen Teil FLT:2] G(jω) auf der horizontalen Achse und dem imaginären Teil auf der vertikalen Achse gezeichnet. Sein Hauptnutzen liegt in der Bewertung der Stabilität im geschlossenen Regelkreis aus dem offenen Regelkreisfrequenzgang, ohne explizit geschlossene Regelkreispole zu berechnen. Der Nyquist-Plot bietet auch Einblick in Verstärkungs- und Phasenränder, was ihn für ein robustes Controller-Design unverzichtbar macht.

Die mathematischen Grundlagen der Nyquist-Plots erfordern ein gründliches Verständnis der komplexen Analyse, insbesondere die Bewertung der Übertragungsfunktionen auf der imaginären Achse und die geometrische Interpretation der Konturabbildung. Ingenieure müssen das Konzept der Nyquist-Kontur, einen geschlossenen Pfad in der komplexen Ebene, der die gesamte rechte Halbebene (RHP) umschließt, erfassen und wie diese Kontur durch die Open-Loop-Transferfunktion G(s)H(s) abbildet, um den Nyquist-Pfad zu erzeugen. Die sich ergebende Drehung des Plots um den kritischen Punkt -1 + 0j zeigt direkt die Systemstabilität an, nach dem Nyquist-Stabilitätskriterium. Dieser Artikel erweitert diese mathematischen Grundlagen und geht durch die Konstruktionslogik, das Argumentprinzip und praktische Fallstricke.

Komplexe Transferfunktionsbewertung

Die Open-Loop-Übertragungsfunktion G(s)H(s) ist eine rationale Funktion der komplexen Variablen s = σ + jω. Für die Frequenzantwortanalyse setzen wir σ = 0, also sjω. Der resultierende Ausdruck G(jω)H(jω) ist eine komplexe Zahl, die von ω abhängt:

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

Hier ist M(ω) die Größe (auch Gain genannt) und φ(ω) der Phasenwinkel. Beim Konstruieren des Nyquist-Plots wird G(jω)H(jω) bei vielen Frequenzen ausgewertet und die komplexen Punkte aufgetragen. Für ein typisches System beginnt der Plot bei ω = 0 (oft an einem Punkt auf der realen Achse) und endet bei ω → ∞ (normalerweise am Ursprung, es sei denn, das System hat reine Integratoren). Die Form der Kurve zeigt wichtige dynamische Eigenschaften auf: zum Beispiel, wie die Phase mit zunehmender Frequenz abnimmt oder wie die Größe abrollt.

Um G(jω)H(jω) mathematisch zu berechnen, betrachten Sie eine allgemeine Übertragungsfunktion, die in Pole und Nullen eingeteilt ist:

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

Wenn s = jω wird jeder Term ein komplexer Vektor vom Null- oder Polort zum Punkt jω auf der imaginären Achse. Die Gesamtgröße ist das Produkt der Längen der Nullvektoren geteilt durch das Produkt der Längen der Polvektoren, multipliziert mit |K|. Die Gesamtphase ist die Summe der Winkel der Nullvektoren minus der Summe der Winkel der Polvektoren plus die Phase von K (was 0 ist, wenn K real und positiv ist). Diese geometrische Interpretation ist von zentraler Bedeutung, um zu verstehen, wie sich die Nyquistkurve verhält, insbesondere in der Nähe von Frequenzen, die den natürlichen Frequenzen des Systems entsprechen.

Mathematische Bauschritte

Die Erstellung eines Nyquist-Plots von Hand (oder mithilfe von Software) beinhaltet die systematische Auswertung von G(jω)H(jω) über einen Frequenzbereich.

  1. Identifizieren Sie die Übertragungsfunktion und schreiben Sie sie in Bode-Form (Pole und Nullen).
  2. Berechnen Sie die Größe und Phase bei einer Reihe von Frequenzen. typischerweise niedrige Frequenzen (nahe 0) und hohe Frequenzen (nahe ∞) gewählt werden, mit Zwischenpunkten in der Nähe von Pause Frequenzen (wo Pol oder Null Beiträge Vorzeichen ändern).
  3. Start bei ω = 0: Der Punkt hängt vom Systemtyp ab. Für Typ 0 (keine Integratoren) beginnt die Kurve bei K (der DC-Verstärkung) auf der positiven reellen Achse. Für Typ 1 (ein Integrator) beginnt die Kurve bei einer unendlichen Größe mit einer Phase von -90°. Für Typ 2 (zwei Integratoren) ist die Phase bei ω = 0 -180°.
  4. Zeichne Punkte bei zunehmenden Frequenzen, die sie glatt verbinden. Wenn die Frequenz zunimmt, bewegt sich die Kurve im Allgemeinen zum Ursprung (da die Größe abnimmt). Die Phase kann sich bei der Phasenkreuzfrequenz, bei der die Größe für die Stabilitätsanalyse entscheidend ist, -180° kreuzen.
  5. Spiegel für negative Frequenzen (bei Verwendung der vollständigen Nyquistkontur). Der Verlauf für ω von -∞ bis 0 ist das komplexe Konjugat des Verlaufs für positive ω, der über die reelle Achse reflektiert wird. Kombiniert mit dem positiven Frequenzzweig (und dem unendlichen Halbkreis, wenn das System einen Nettopol-Null-Überschuss hat) wird der vollständige Nyquistpfad erhalten.
  6. Identifizieren Sie die Anzahl der Umkreisungen von -1 + 0j, um das Stabilitätskriterium anzuwenden.

Mathematisch kann die Nyquistkurve für positive Frequenzen als ω ∈ [0,∞] parametrisiert werden. Für Systeme mit relativem Grad (n - m) ≥ 1 endet die Kurve am Ursprung und nähert sich ihr mit einer Phase von -(n - m) x 90 °. Für den relativen Grad 0 (gleiche Anzahl von Polen und Nullen) ist der Endpunkt eine endliche, nicht-null-komplexe Zahl bei ω → ∞.

Das Nyquist Stabilitätskriterium im Detail

Das Nyquist-Stabilitätskriterium ist eine leistungsfähige Methode, um die Stabilität des geschlossenen Kreislaufs aus dem offenen Frequenzgang zu bestimmen. Es nutzt ss-Prinzip, das durch eine Funktion sF in der F den Ursprung umgibt, wobei Z die Anzahl der Nullen von F]s die Anzahl der Pole von F(s)] und die Kontur Γ ist die Funktion Nyquist-Kontur, die die gesamte rechte Halbebene (RHP) umschließt, indem sie die imaginäre Achse von −]j∞ durchquert und dann einen Halbkreis mit unendlichem Radius nach

Da F(s) mit der Open-Loop-Übertragungsfunktion durch F(s) = 1 + G(s)H(s) verwandt ist, ist die Karte von F(s) durch die Nyquist-Kontur einfach die Nyquist-Plot von G(s)H(s) nach rechts verschoben. Entsprechend entsprechen Umkreisungen des Ursprungs in der Fj in der G(s)H(s)-Ebene genau den Umkreisungen des Punktes −1 + 0G(s)H(s)-Ebene. Das Nyquist-Stabilitätskriterium lautet dann:

N = Z − P

in dem Fall:

  • N = Anzahl der Einkreisungen von −1 + 0j durch den Nyquist-Plot (gezählt im Uhrzeigersinn als positiv, gegen den Uhrzeigersinn als negativ).
  • Z = Anzahl der geschlossenen Pole im RHP (Null von 1 + G(s)H(s) im RHP).
  • P = Anzahl der offenen Pole im RHP.

Für ein stabiles Closed-Loop-System benötigen wir Z = 0. Das Closed-Loop-System ist also stabil, wenn und nur wenn der Nyquist-Plot von G(s)H(s) den Punkt genau -P mal (d.h. gegen den Uhrzeigersinn sind Umkreisungen gleich ]P wenn wir im Uhrzeigersinn als positiv zählen).

Anwendung des Kriteriums: Encirclement Counting

Die genaue Berechnung der Umkreisungen umfasst die Analyse des Nyquistpfades über die gesamte Kontur, einschließlich der Abbildung des unendlichen Halbkreises. Bei Systemen mit relativem Grad ≥ 1 bildet der unendliche Halbkreis einen kleinen Kreis um den Ursprung in der Ebene G(s)H(s) ab (oder einen festen Punkt, wenn der relative Grad Null).

  1. Zeichne den Nyquist-Plot für ω von 0 bis ∞.
  2. Reflektieren Sie es über die reelle Achse für ω von 0 bis -∞.
  3. Wenn das System Pole auf der imaginären Achse hat (z. B. Integratoren), ziehen Sie die Nyquist-Kontur mit kleinen Halbkreisen nach rechts, um diese Pole zu vermeiden, und kartieren Sie diese Einbuchtungen zu großen Bögen in der Ebene FLT: 0 G (s) H (s) [ FLT: 1].
  4. Zählen Sie die Nettoanzahl der Umkreisungen im Uhrzeigersinn von -1, indem Sie die vollständige geschlossene Kurve (einschließlich der unendlichen Halbkreisabbildung) verfolgen.

Mathematisch kann die Einkreisungszahl N über das Wicklungszahlintegral berechnet werden:

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

Die Haupterkenntnis ist, dass Einkreisungen auftreten, wenn die Phase von G(jω)H(jω) sich -180° (d.h. den Punkt -1 auf der realen Achse) kreuzt, während die Größe ≥ 1 ist. Diese Bedingung definiert den Gain Margin und Phasen Margin.

Praktische Konstruktionsbeispiele

Um das Verständnis zu verfestigen, betrachten wir ein einfaches System erster Ordnung und ein System zweiter Ordnung mit Dämpfung.

Erstbestellungssystem: G(s) = 1/(s+1)

Hier G(jω) = 1/(jω+1)M(ω) = 1/√(1+ω2)φ(ω) = −tan−1(ω). Bei ω = 0 ist der Punkt (1,0). Da ω → ∞ der Punkt sich (0,0) von einem Winkel von -90° nähert. Die Nyquistkurve ist ein Halbkreis in der unteren realen Achse (da die Phase immer negativ ist). Die Nyquistkurve ist ein Halbkreis in der unteren halben Ebene (da die Phase immer negativ ist) und erstreckt sich von der positiven realen Achse bis zum Ursprung. Diese Kurve ist nicht nicht] nicht stabil für jeden positiven Gewinn K wenn wir G(s) mit KG(s) ersetzen. Für sehr große K

Zweitbestätiges System mit Dämpfung: G(s) = ω n2/(s2 + 2ζω n s + ω n2)

Für ein System zweiter Ordnung mit natürlicher Frequenz ωn beginnt der Nyquist-Plot bei (1, 0) für ω = 0 und endet bei (0, 0) mit einer Phase von ω = 0. Die Form hängt kritisch von ζ ab. Für unterdämpfte Systeme (ζ < 1), the curve may exhibit a loop that crosses the negative real axis at a frequency ωpc mit einer Phase von ζ. Bei dieser Frequenz ist die Größe Mpc> 1, dann ist der Nyquist-Plot umkreist -1 und das geschlossene System ist instabil. Dies entspricht einer Gewinnspanne von weniger als 1 (0 dB). Umgekehrt, wenn ζ groß genug ist, so dass pcpc] analytisch gefunden werden kann: Der Phasenübergang erfolgt bei ω n2/]n = ω2/n[FLT

Um einen richtigen Phasenübergang zu sehen, betrachten Sie ein System mit einem zusätzlichen Pol bei niedriger Frequenz: G(s) = 1/((s+1)(s2 + 2ζ s + 1)] Hier kann die Phase -180° bei einigen endlichen ω kreuzen, und die Handlung kann -1 umkreisen, wenn der Gewinn ausreichend hoch ist.

Nyquist, Bode und Nichols Plots verbinden

Der Nyquist-Plot ist eng mit dem Bode-Plot (Magnitude und Phase) und dem Nichols-Diagramm verbunden. Der Bode-Plot liefert die Größe in dB und die Phase in Grad als Funktion der Log-Frequenz. Der Nyquist-Plot ist eine parametrische Kombination aus der Größe und Phase: Jeder Punkt auf der Nyquist-Kurve entspricht einer Größe und Phase bei einem gegebenen ω. Der Gain-Rand ist der Faktor, um den die Verstärkung erhöht werden kann, bevor die Größe bei der Phasenkreuzungsfrequenz 1 wird (d. h. der Abstand von -1 zur Nyquist-Kurve entlang der negativen realen Achse). Der Phase-Rand ist die zusätzliche Phasenverzögerung, die bei der Verstärkungskreuzungsfrequenz (wobei die Größe = 1) erforderlich ist, um die Phase zu erzeugen -180°. Diese Ränder sind direkt aus dem Nyquist-Plot ablesbar: Der Verstärkungsrand ist 1/

Das Nichols-Diagramm transformiert den Nyquist-Diagramm in einen kartesischen Diagramm der Größe (dB) gegenüber der Phase (deg). Diese Transformation ist nützlich für den Entwurf von Reglern, da Verstärkungs- und Phasenränder leicht zu erkennen sind und die Antwort mit geschlossenen Schleifen mit M-Kreisen gelesen werden kann. Alle drei Darstellungen sind mathematisch äquivalent; die Wahl einer hängt vom Kontext und der Präferenz des Ingenieurs ab.

Einschränkungen und fortgeschrittene Überlegungen

Während der Nyquist-Plot leistungsstark ist, hat er Einschränkungen, die Ingenieure anerkennen müssen. Erstens wird die Nyquist-Kurve für Systeme mit Zeitverzögerungen unendlich und stark gekrümmt, was eine sorgfältige Handhabung der unendlichen Halbkreis-Mapping erfordert. Die Verzögerung fügt eine Phasenverschiebung hinzu -ωT, die sich linear ansammelt, wodurch der Nyquist-Plot je nach Verzögerung nach innen oder außen spiralförmig wird. Stabilitätsanalyse erfordert dann das modifizierte Nyquist-Kriterium, das die Verzögerung unter Verwendung des Konzepts der Nyquist-Kontur berücksichtigt, die den Verzögerungspol bei Unendlichkeit ausschließt.

Zweitens ist der Nyquist-Plot, wie beschrieben, für Systeme mit kontinuierlicher Zeit. Für Systeme mit diskreter Zeit (z-Domäne) wird die Nyquist-Kontur zum Einheitskreis, und das Stabilitätskriterium untersucht Umkreisungen von -1 in der w-Ebene oder die Verwendung der bilinearen Transformation zur Abbildung in die s-Ebene. Die mathematischen Grundlagen bleiben ähnlich, aber die Kontur ändert sich.

Drittens geht das Nyquist-Kriterium davon aus, dass das Closed-Loop-System stabil ist, wenn und nur wenn der Nyquist-Plot-Umkreis -1 entsprechend ist. Dies erfordert jedoch, dass der Charakter der Open-Loop-Pole im RHP bekannt ist a priori Wenn die Open-Loop-Transferfunktion Pole auf der imaginären Achse hat (einschließlich reiner Integratoren), muss die Nyquist-Kontur eingerückt werden, um diese Pole zu vermeiden, und die Abbildung dieser Einbuchtungen erzeugt unendliche Bögen, die in die Einkreisungszahl einbezogen werden müssen. Dies kann für Ingenieure, die neu in dem Verfahren sind, schwierig sein.

Schließlich ist der Nyquist-Plot ein Frequenzdomänen-Tool; er bietet notwendige und ausreichende Bedingungen für Stabilität nur für lineare zeitinvariante (LTI) Systeme.Für nichtlineare oder zeitvariable Systeme gilt das Nyquist-Kriterium nicht direkt, obwohl Erweiterungen wie das Kreiskriterium für bestimmte Klassen von Nichtlinearitäten existieren.

Schlussfolgerung

Die Konstruktion eines Nyquist-Plots ist tief in der komplexen Analyse verwurzelt, insbesondere der Abbildung der Nyquist-Kontur durch die Open-Loop-Transferfunktion und das Argumentprinzip. Durch die Beherrschung der mathematischen Schritte - Bewertung von G(jω)H(jω), die Darstellung der komplexen Punkte und das Zählen von Umkreisungen von -1 - Ingenieure erhalten eine leistungsstarke Methode zur Beurteilung der Stabilität im geschlossenen Kreislauf und zum Entwerfen robuster Steuerungen. Das Nyquist-Kriterium bietet intuitive Einblicke in Verstärkungs- und Phasenränder und seine Verbindung zu Bode- und Nichols-Plots macht es zu einem vielseitigen Werkzeug im Toolkit des Steuerungsingenieurs. Das Verständnis dieser Grundlagen ermöglicht es Ingenieuren, den Nyquist-Plot nicht nur korrekt anzuwenden, sondern auch seine Ergebnisse mit Vertrauen in reale Anwendungen zu interpretieren, von Luft- und Raumfahrtsystemen bis hin zur industriellen Prozesssteuerung.

Für weitere Informationen siehe Standard-Lehrbücher wie Wikipedias Artikel über das Nyquist-Stabilitätskriterium, Control Tutorials für MATLAB und Simulink oder den klassischen Text Modern Control Engineering von Ogata.