Einführung in das Routh-Hurwitz-Kriterium

Das Routh-Hurwitz-Kriterium ist ein Eckpfeiler der klassischen Regeltheorie und bietet eine systematische Methode zur Bewertung der Stabilität linearer zeitinvarianter (LTI) Systeme, ohne die Wurzeln der charakteristischen Gleichung explizit zu berechnen. Diese algebraische Technik wurde 1874 unabhängig von Edward John Routh und 1895 von Adolf Hurwitz entwickelt und ist nach wie vor in der Regeltechnik weit verbreitet, da sie eine schnelle, zuverlässige Möglichkeit bietet, um festzustellen, ob ein System stabil und vorhersagbar reagiert oder Oszillationen und Instabilität aufweist. Stabilität ist für das Systemdesign von grundlegender Bedeutung - ein instabiles System kann zu gefährlichen Überschwingen, anhaltenden Oszillationen oder sogar zu katastrophalen Ausfällen führen, insbesondere in Anwendungen wie Luft- und Raumfahrt, Automobilsteuerungen, chemische Prozesse und Robotik.

Durch die Analyse der Koeffizienten des charakteristischen Polynoms eines Systems identifiziert das Routh-Hurwitz-Kriterium die Anzahl der Pole (Wurzeln der charakteristischen Gleichung), die in der rechten Hälfte der komplexen Ebene liegen. Pole mit positiven Realteilen verursachen instabiles Verhalten, während Pole mit negativen Realteilen zu stabilen, konvergenten Reaktionen beitragen. Die Methode erfordert keine Lösung von Polynomen hohen Grades, was sie für die manuelle Analyse und das Frühstadium praktisch macht. Darüber hinaus helfen ihre Erkenntnisse über transientes Verhalten - wie z. B. Einschwingzeit, Überschwingen und Schwingungsfrequenz - Ingenieuren, die Leistungsspezifikationen zu erfüllen, bevor sie sich an detaillierte Simulationen oder Hardware-Prototypen binden.

Dieser Artikel wird zuerst die mathematische Grundlage des Routh-Hurwitz-Kriteriums erklären, dann führen Sie durch die Konstruktion des Routh-Arrays mit einem Schritt-für-Schritt-Beispiel. Danach werden wir untersuchen, wie das Kriterium die Systemreaktion und das vorübergehende Verhalten vorhersagt, einschließlich wichtiger Merkmale wie Überschwingen und Einschwingzeit. Ein praktisches numerisches Beispiel wird den Prozess veranschaulichen und wir werden die Einschränkungen und komplementären Werkzeuge der Methode diskutieren. Am Ende werden Sie ein solides Verständnis dafür haben, wie Sie das Routh-Hurwitz-Kriterium als praktisches, maßgebliches Werkzeug in der Steuerungssystemanalyse anwenden können.

Was ist das Routh-Hurwitz-Kriterium?

Das Routh-Hurwitz-Kriterium ist eine notwendige und ausreichende Bedingung für die Stabilität in LTI-Systemen. Es arbeitet mit der charakteristischen Gleichung des Systems, die sich aus dem Nenner der Closed-Loop-Übertragungsfunktion nach dem Setzen auf Null ergibt.

Das Kriterium besagt, dass für das System stabil zu sein (dh alle Pole haben negative reelle Teile), alle Koeffizienten ai müssen ungleich Null sein und das gleiche Vorzeichen haben (typischerweise positiv), und die Routh-Array aus diesen Koeffizienten aufgebaut keine Vorzeichenänderungen in der ersten Spalte haben.

Es ist wichtig zu beachten, dass das Routh-Hurwitz-Kriterium rein algebraisch ist - es gibt nicht die tatsächlichen Polpositionen an. Stattdessen liefert es eine Anzahl von Polen mit der rechten Hälfte der Ebene (RHP). Dies macht es von unschätzbarem Wert, wenn das charakteristische Polynom von hoher Ordnung ist (z. B. viert, fünft oder höher), wo die Wurzelfindung rechenintensiv oder analytisch nicht machbar ist.

Das Kriterium ist direkt mit der Hurwitz-Matrix verknüpft, einer quadratischen Matrix, die aus den Koeffizienten des charakteristischen Polynoms gebildet wird. Alle führenden Haupt-Molls der Hurwitz-Matrix müssen für die Stabilität positiv sein - das ist die Hurwitz-Bedingung. Rouths tabellarische Methode vereinfacht diesen Determinanten-Prüfprozess in das einfachere Routh-Array.

Bau des Routh Array: Schritt für Schritt

Der Aufbau des Routh-Arrays ist ein systematisches Verfahren. Wir werden es anhand eines konkreten Beispiels Polynoms veranschaulichen:

Dieses Polynom hat bekannte Wurzeln bei s = -1, -2, -3, -4 (alle negativ), also ist das System stabil.

Schritt 1: Schreiben Sie das Polynom in Standardform

Für unser Beispiel sind die Koeffizienten:
a4 = 1 (s4), a33 = 10 (s2, a11, a0 = 24 (konstanter Begriff).

Schritt 2: Initialisieren Sie die ersten beiden Zeilen

Die erste Zeile des Routh-Arrays enthält die Koeffizienten von geraden indizierten Termen (beginnend von sn): Zeile s4: [1, 35, 24]. Die zweite Zeile enthält die Koeffizienten von ungeraden indizierten Termen: Zeile s3: [10, 50, 0]. (Füllen Sie fehlende Terme mit Nullen aus.)

Also:
Reihe s4: 1, 35, 24
Reihe s3: 10, 50, 0

Schritt 3: Berechnen Sie nachfolgende Zeilen mit der Rhombus-Regel

Für ein Element in Zeile s2, erste Spalte, lautet die Formel:

Aber es ist einfacher, die Standardformel anzuwenden: Für Zeile k, Spalte j (ab Spalte 1) ist das Element:

Lasst es uns manuell machen, um Klarheit zu schaffen.

]Zeile s2:3
R(s,1 = -[(1x50) - (10x35)] / 10 = -[50 -350]/10 = 30
R(s4,2 = 50, R(s3
2[[0 - 240]/10 = 24
Drittes Element: j = 3, aber R(s4],4 nicht vorhanden (Behandeln als 0), R(s]3,4=0. Also R(s2,3 = -[(1x0)

Die Zeile s2: 30, 24, 0.

Zeilen s13 und s2
R(s1,1 / 30 = -[240 - 1500]/30 = -[1260]/30 = 42
R(s1,2 = -[10×0) - (30×0) ] / 30 = 0.

Die Zeile s1: 42, 0.

Row s0 und s1
R(s0,1 = -[(30×0) - (42×0)] / 42 = 0. Aber die letzte Zeile ist immer der konstante Term aus dem Polynom, der 24 sein sollte. Tatsächlich gibt es eine bekannte Regel: Das letzte Element der letzten Zeile ist gleich a0 (der konstante Term). Wir können berechnen als: R(s0,1 = -[(30×0) - (42×24)] / 42 = -[-1008]/42 = 24!

Das komplette Routh-Array ist also:

s413524
s310500
s230240
s14200
s02400

Schritt 4: Untersuchen Sie die erste Kolonne

Die erste Spalte enthält: 1, 10, 30, 42, 24. Alle sind positiv, und es gibt keine Anzeichen für Veränderungen.

Sonderfälle

Zwei Sonderfälle können während der Routh-Array-Konstruktion auftreten:

  • Null in der ersten Spalte: Wenn das erste Element einer Zeile Null ist, aber die Zeile nicht alle Nullen ist, ersetzen Sie die Null durch ein kleines Epsilon (ε), dann fahren Sie fort. Das Vorzeichen von ε bestimmt Vorzeichenänderungen. Dies zeigt eine Randstabilität oder Instabilität an, abhängig vom Rest des Arrays.
  • Eine ganze Reihe von Nullen: Dies geschieht, wenn es symmetrisch platzierte Wurzeln gibt (z. B. Paare von imaginären Wurzeln oder symmetrische Paare um die reale Achse). In diesem Fall bilden Sie ein Hilfspolynom aus der Zeile über der Nullreihe und differenzieren Sie es, um Koeffizienten für die fehlende Zeile zu erhalten. Dies zeigt das Vorhandensein eines Paares von rein imaginären Wurzeln oder eines Paares von reellen Wurzeln mit entgegengesetzten Vorzeichen (dh Instabilität).

Diese Sonderfälle erweitern die Anwendbarkeit des Kriteriums auf Systeme mit Polen auf der gedachten Achse oder symmetrischen Mustern.

Vorhersage von Systemreaktion und vorübergehendem Verhalten

Das Routh-Hurwitz-Kriterium ist nicht auf ein binär stabiles/instabiles Urteil beschränkt. Indem es die Anzahl der Pole mit positiven Realteilen angibt, informiert es direkt über das transiente Verhalten des Systems. Jeder Pol in der rechten Halbebene trägt einen exponentiell wachsenden Modus bei, wodurch die transiente Antwort auseinandergeht. Umgekehrt führen alle Pole in der linken Halbebene zu verfallenden Exponentialen, die schließlich einen stationären Zustand erreichen.

Das Kriterium gibt jedoch nicht die genauen Pole-Positionen an, so dass es Überschwingen, Einschwingzeit oder Anstiegszeit nicht genau quantifizieren kann.

  • Wenn die erste Spalte keine Vorzeichenänderungen hat: Alle Pole haben negative reelle Teile. Die vorübergehende Antwort wird stabil sein und sich schließlich in einen stabilen Zustand einpendeln. Aber es könnte langsamen oder oszillatorischen Zerfall geben (fast imaginäre Pole), den das Kriterium nicht erkennt - nur das Zeichen des realen Teils ist garantiert.
  • Wenn es Vorzeichenänderungen gibt: Die Anzahl der Vorzeichenänderungen entspricht der Anzahl der RHP-Pole. Das System ist instabil und die vorübergehende Antwort wird unbegrenzt wachsen. Die Anzahl der RHP-Pole deutet auf die Art der Instabilität hin: Ein einzelner positiver Pol führt zu monotonem exponentiellem Wachstum, während ein Paar komplexer RHP-Pole wachsende Schwingungen erzeugt.
  • Marginale Stabilität: Wenn die charakteristische Gleichung Wurzeln auf der imaginären Achse hat (rein imaginäre oder mehrere Pole auf der Achse), wird das Routh-Array eine Reihe von Nullen erzeugen.

Durch die Kombination des Routh-Hurwitz-Ergebnisses mit der Kenntnis der Systemordnung können Ingenieure fundierte Vorhersagen treffen. Für ein System zweiter Ordnung reduziert sich das Kriterium darauf, dass alle Koeffizienten positiv sein müssen (was sowohl notwendig als auch für die Stabilität ausreichend ist). Für Systeme höherer Ordnung bietet das Array einen Stabilitätsmargeindikator - je näher die Elemente der ersten Spalte an Null sind, desto näher ist das System an Instabilität. Dies ist insbesondere bei der Anpassung von Verstärkungsparametern in einem Regelkreis nützlich.

Beziehung zu vorübergehenden Merkmalen

Das vorübergehende Ansprechen eines stabilen Systems wird in erster Linie durch die dominanten Pole bestimmt, die der imaginären Achse am nächsten sind. Das Routh-Hurwitz-Kriterium ordnet diese Pole nicht direkt, kann aber bei der Abstimmung helfen. Wenn das Kriterium beispielsweise zeigt, dass eine Erhöhung eines Reglergewinns schließlich zu einer Vorzeichenänderung in der ersten Spalte führt, markiert es die Verstärkungsspanne. Ingenieure verwenden dies oft, um proportionale Gewinne in P-, PI- oder PID-Reglern festzulegen, ohne jedes Mal das gesamte System simulieren zu müssen.

Wichtige transiente Metriken, die von der Stabilität beeinflusst werden, sind:

  • Overshoot – Der maximale Spitzenwert der Antwort im Verhältnis zum Endwert. Overshoot wird typischerweise durch komplexe konjugierte Pole mit einem niedrigen Dämpfungsverhältnis verursacht. Während das Routh-Hurwitz-Kriterium keine Dämpfung berechnet, entspricht ein System mit einer Reihe von Nullen (marginale Stabilität) oft einem Dämpfungsverhältnis nahe Null, was ein hohes Überschwingen bedeutet.
  • Einstellzeit – Die Zeit, die benötigt wird, damit die Antwort innerhalb eines bestimmten Prozentsatzes (z. B. 2%) des endgültigen Wertes bleibt. Die Einstellzeit ist umgekehrt proportional zum realen Teil des dominanten Pols. Das Routh-Array kann diesen realen Teil nicht angeben, aber wenn das Array zeigt, dass sich alle Pole in der linken Halbebene befinden, ist die Einstellzeit endlich; wenn ein Pol sehr nah an der imaginären Achse liegt, wird die Einstellzeit lang sein.
  • Rise Time – Die Zeit, um von 10% auf 90% des Endwertes zu gehen. Dies hängt mit der natürlichen Frequenz der dominanten komplexen Pole zusammen. Auch hier liefert das Routh-Hurwitz-Kriterium keinen direkten numerischen Wert, sondern zeigt an, ob das System stabil genug ist, um eine gut definierte Anstiegszeit zu haben.
  • Schwingungen – Die Frequenz von transienten Schwingungen wird durch die imaginären Teile der komplexen Pole bestimmt. Das Kriterium zeigt keine imaginären Teile, aber das Vorhandensein einer Reihe von Nullen zeigt rein imaginäre Pole an, die anhaltende Schwingungen in der transienten Antwort erzeugen.

Das Routh-Hurwitz-Kriterium fungiert somit als Torwächter: Es sorgt für Stabilität als Voraussetzung für eine sinnvolle Transientenanalyse mit anderen Tools wie Root Locus, Bode-Plots oder Schrittreaktionssimulationen.

Praktisches Beispiel: System dritter Ordnung

Nehmen wir das Routh-Hurwitz-Kriterium auf ein komplexeres Beispiel mit unbekannten Parametern, die Closed-Loop-Übertragungsfunktion mit Einheitsrückkopplung.

Wir wollen den Bereich von K bestimmen, für den das System stabil ist.

Zeilen:
s3 s2: 7, K
1: ] s0 Eigentlich einfacher: Das Element der letzten Zeile ist der konstante Term K (vorausgesetzt, die obige Zeile ist nicht alle Nullen).

s3110
s27K
s1(K-70)/70
s0K0

Für Stabilität müssen alle Elemente der ersten Spalte positiv sein: 1 > 0, 7 > 0, (K-70)/7 > 0, K > 0, also K > 70. Die untere Grenze ist K > 70. Wenn K = 70 ist, wird die Zeile s1 zu Null, was ein Spezialfall ist, der auf Randstabilität hinweist (rein imaginäre Wurzeln). Das System ist also stabil für K > 70. Dies sagt dem Ingenieur, dass Gewinne über 70 einen stabilen Transienten ergeben, aber wenn K sich 70 nähert, werden die transienten Schwingungen ungedämpft. Die transiente Antwort für K = 200 hat zum Beispiel zwei komplexe konjugierte Pole in der linken Halbebene und ein reeller Pol, was zu einer Überschwingenden, aber schließlich beruhigenden Antwort führt. Das Routh-Hurwitz-Kriterium sagt hier explizit die Stabilitätsgrenze voraus.

Transient Response Implikationen für das Beispiel

Bei K knapp über 70 ist das System stabil, aber die dominanten Pole liegen sehr nahe an der imaginären Achse. Die transiente Reaktion zeigt einen sehr langsamen Abklingen und einen großen Überschwinger. Die Einschwingzeit wird lang sein. Mit zunehmender K bewegen sich die Pole nach links (im Allgemeinen), was die Dämpfung und Einschwingzeit verbessert. Wenn K jedoch zu hoch ist, würde das Routh-Array für dieses spezielle System keine zusätzlichen Vorzeichenänderungen zeigen - Stabilität bleibt bestehen. In der Praxis könnte sehr hohes K andere Probleme wie Aktorsättigung oder Rauschverstärkung verursachen, aber das Routh-Hurwitz-Kriterium überprüft nur die Schleifenstabilität.

Das Kriterium liefert somit wichtige Informationen: den Verstärkungsbereich für Stabilität und die Nähe zur Instabilität. Durch die Auswertung der Elemente der ersten Spalte können Ingenieure Verstärkungsmargen berechnen und Steuerungen entwerfen, die eine marginale Stabilität vermeiden.

Einschränkungen des Routh-Hurwitz-Kriteriums

Trotz seiner Leistungsfähigkeit hat das Routh-Hurwitz-Kriterium wichtige Einschränkungen, die jeder Steuerungsingenieur verstehen muss:

  • Gilt nur für LTI-Systeme: Es kann nicht direkt für zeitvariable, nichtlineare oder Sampled-Data-Systeme verwendet werden (obwohl Erweiterungen für zeitdiskrete Systeme über bilineare Transformationen existieren).
  • Bestimmt nicht die transiente Leistung: Das Kriterium gibt keine Informationen über Überschwingen, Spitzenzeit, Anstiegszeit oder genaue Einschwingzeit.
  • Erfordert ein bekanntes charakteristisches Polynom: Die Koeffizienten müssen genau bekannt sein. In der Praxis können Unsicherheiten in Parametern eine robuste Stabilitätsanalyse erfordern (z. B. unter Verwendung des Kharitonov-Theorems).
  • Wird für Systeme sehr hoher Ordnung umständlich: Manuelle Berechnungen für Polynome zehnter Ordnung sind fehleranfällig, obwohl Software (MATLAB, Python) die Array-Konstruktion automatisiert.
  • Kann Zeitverzögerungen nicht direkt handhaben: Systeme mit Zeitverzögerungen haben transzendentale charakteristische Gleichungen - das Routh-Array kann nicht direkt angewendet werden, es sei denn, Padé-Näherungen verwendet werden.

Um diese Einschränkungen zu überwinden, verwenden Ingenieure üblicherweise das Routh-Hurwitz-Kriterium in Verbindung mit anderen Methoden. So zeigen Wurzelortdiagramme, wie sich geschlossene Pole mit Verstärkung bewegen, was genaue transiente Informationen liefert. Das Nyquist-Kriterium behandelt Zeitverzögerungen und Nichtlinearitäten. Dennoch dient das Routh-Array oft als eine schnelle erste Überprüfung - insbesondere bei der vorläufigen Konstruktion oder Fehlersuche.

Ergänzende Tools

Das Routh-Hurwitz-Kriterium ist "notwendig und ausreichend" für die Stabilität von LTI-Polynomen, aber für die Leistungsanalyse verlassen sich Ingenieure auf:

  • Root Locus: Zeigt die Bahn der Pole als Parameter (z. B. Verstärkung) an. Es zeigt direkt Dämpfungsverhältnis und Eigenfrequenz an.
  • Bode Plots und Nyquist Diagramme: Bieten Verstärkungs- und Phasen-Marge, die Indikatoren für relative Stabilität und transiente Antworteigenschaften sind.
  • Schrittreaktionssimulation: Gibt die tatsächlichen transienten Metriken (Überschwingen, Einschwingzeit) für ein gegebenes System an.
  • State-Space-Methoden Die Eigenwerte der Zustandsmatrix (die die Pole sind) können numerisch berechnet werden, und die Methode von Lyapunov erweitert die Stabilitätsanalyse auf nichtlineare Systeme.

Das Routh-Hurwitz-Kriterium nimmt somit einen historischen und praktischen Sweet Spot ein: Es ist mathematisch genug, um streng zu sein, aber einfach genug, um es während einer Entwurfsprüfung auf einem Whiteboard zu tun.

Schlussfolgerung

Das Routh-Hurwitz-Kriterium bietet eine klare, algebraische Möglichkeit, die Stabilität des Steuerungssystems zu beurteilen, indem das charakteristische Polynom analysiert wird, ohne nach Wurzeln zu lösen. Durch die Konstruktion des Routh-Arrays können Ingenieure die Anzahl der instabilen Pole bestimmen, die Empfindlichkeit des Systems gegenüber Parameteränderungen verstehen und Stabilitätsgrenzen identifizieren, die sich direkt auf das transiente Verhalten auswirken. Obwohl es keine genauen transienten Metriken wie Überschwingen oder Einschwingzeit bietet, dient es als grundlegendes Screening-Tool, das eine tiefere Analyse mithilfe von Wurzelort, Frequenzgang oder Simulation ermöglicht.

Die Beherrschung des Routh-Hurwitz-Kriteriums befähigt Sie schnell zu beurteilen, ob sich ein Steuerungssystem unter transienten Bedingungen zuverlässig verhält. Ob Sie einen einfachen Proportionalregler entwerfen oder ein komplexes Multi-Loop-System abstimmen, das Routh-Array ist ein zuverlässiger Ausgangspunkt. Für weitere Informationen sollten Sie das Stabilitätskriterium von Routh-Hurwitz auf Wikipedia für zusätzliche mathematische Tiefe erkunden oder die Steuerungstutorials der University of Michigan auf Routh Array für weitere bearbeitete Beispiele konsultieren. Darüber hinaus bietet MathWorks-Dokumentation Einblick in Softwareimplementierungen. Durch die Integration des Routh-Hurwitz-Kriteriums in Ihren Steuerungsdesign-Workflow stellen Sie sicher, dass die von Ihnen erstellten Systeme nicht nur stabil sind, sondern auch vorhersehbare, leistungsstarke transiente Reaktionen aufweisen.