Einführung in das Routh-Hurwitz-Stabilitätskriterium

Das Routh-Hurwitz-Stabilitätskriterium ist ein Eckpfeiler der Regeltheorie und der Polynomanalyse. Es bietet eine systematische Methode zur Bestimmung der Anzahl der Wurzeln eines realen Polynoms, die in der rechten Hälfte der komplexen Ebene liegen - ohne explizite Wurzelberechnung zu erfordern. Für Ingenieure und Mathematiker, die mit komplexen Polynomen arbeiten, verwandelt dieses Kriterium ein ansonsten schwieriges algebraisches Problem in ein einfaches tabellarisches Verfahren.

Unter einem "komplexen Polynom" wird hier ein Polynom mit reellen Koeffizienten (wie es in Regelanwendungen üblich ist) verstanden, dessen Wurzeln jedoch komplex sein können. Die Stabilität vieler physikalischer Systeme - von elektronischen Schaltungen bis hin zu mechanischen Verknüpfungen - hängt vom Vorzeichen der realen Teile der charakteristischen Polynomwurzeln ab. Das Routh-Hurwitz-Kriterium beantwortet diese Frage effizient, auch für Polynome höherer Ordnung, die für den Faktor unpraktisch sind.

Dieser Artikel erweitert die ursprüngliche Zusammenfassung und bietet einen detaillierten Überblick über die Methode, Beispiele, Sonderfälle und praktische Anwendungen. Am Ende werden die Leser in der Lage sein, die Routh-Hurwitz-Technik anzuwenden, um komplexe Polynome mit Sicherheit zu vereinfachen und zu analysieren.

Was sind komplexe Polynome?

In der Steuerungstechnik und Mathematik ist ein komplexes Polynom typischerweise ein Ausdruck der Form:

P(s) = an sn + an−1n−1 + ... + a1 s + a

Dabei sind die Koeffizienten ai reelle Zahlen und s ist eine komplexe Variable. Der Begriff “komplex” bezieht sich nicht auf die Koeffizienten (die reell bleiben), sondern auf die Tatsache, dass die Wurzeln der Gleichung P(s) = 0 im Allgemeinen komplexe Zahlen sind. Diese Wurzeln bestimmen das dynamische Verhalten eines Systems, das durch das Polynom beschrieben wird.

Zum Beispiel ist die charakteristische Gleichung eines Systems zweiter Ordnung wie ein Masse-Feder-Dämpfer ms2 + bs + k = 0 Je nach Dämpfungsverhältnis können die Wurzeln reell und negativ (überdämpft) oder ein komplexes konjugiertes Paar mit negativem Realteil (unterdämpft) sein. Das Routh-Hurwitz-Kriterium bietet eine Möglichkeit zu überprüfen, ob all-Wurzeln negative Realteile haben - ohne die Gleichung zu lösen.

Warum direkte Wurzelberechnung unpraktisch ist

Für Polynome der Ordnung fünf oder höher gibt es selten exakte analytische Wurzellösungen (Abel-Ruffini-Theorem). Numerische Methoden können Wurzeln annähern, sind jedoch empfindlich gegenüber Koeffizientenunsicherheit und können subtile Anzeichen von Instabilität übersehen. Das Routh-Hurwitz-Kriterium hingegen gibt eine definitive Ja/Nein-Antwort auf das Vorzeichen von Realteilen, wobei nur die Polynomkoeffizienten verwendet werden. Dies macht es für die vorläufige Stabilitätsanalyse und die Überprüfung von Designrändern unerlässlich.

Das Routh-Hurwitz-Stabilitätskriterium

Das Kriterium, das unabhängig von Edward John Routh (1876) und Adolf Hurwitz (1895) entwickelt wurde, besagt, dass für ein Polynom mit reellen Koeffizienten alle Wurzeln negative reelle Teile haben, wenn und nur wenn alle führenden Hauptminderwerte der Hurwitz-Matrix (oder äquivalenterweise die erste Spalte des Routh-Arrays) positiv sind.

Bau des Routh Array

Angesichts des Polynoms:

P(s) = ansnn−1sn−1 + ... + a1s + a (mit an > 0

Das Routh-Array ist wie folgt aufgebaut:

  1. Zeile snn, n-2]n-1] enthält Koeffizienten ungerader indexierter Potenzen ], an-3, .... Alle fehlenden Koeffizienten sind mit Nullen gefüllt.
  2. nachfolgende Zeilensn−2nn−1n−1n−4n−5n−1n−1,
    und so weiter, bis die restlichen Elemente Null werden.
  3. Weiter, bis nur die konstante Zeile übrig bleibt. Die letzte Zeile entspricht s0 und enthält ein einzelnes Element, das der konstante Begriff ist (typischerweise a0.
  4. Prüfen Sie die erste Spalte. Wenn alle Einträge in der ersten Spalte des Arrays positiv sind (und keiner Null ist), dann liegen alle Wurzeln in der linken Halbebene (das System ist stabil).

Sonderfälle im Routh Array

Es ergeben sich zwei besondere Situationen, die einen sorgfältigen Umgang erfordern:

  • Null in der ersten Spalte. Wenn das erste Element einer Zeile Null ist, der Rest der Zeile aber ungleich Null ist, ersetzen Sie die Null durch eine kleine positive Zahl ε und fahren fort. Das Vorzeichen von ε bestimmt das Stabilitätsergebnis. Alternativ kann man die Hilfspolynom Methode verwenden (siehe unten).
  • Ganze Reihe von Nullen. Dies zeigt an, dass das Polynom symmetrisch liegende Wurzeln enthält (Wurzelpaare, die Negative voneinander oder komplexe Konjugate mit dem Null-Realteil sind).

Beispiel: Polynom vierter Ordnung

Betrachten Sie P(s) = s4 + 2s3 + 3s2 + 4s + 5.

  • Zeile s4: 1, 3, 5
  • Zeile s3: 2, 4, 0
  • Zeile s2b1 = (2·3 - 1·4)/2 = (6-4)/2 = 1, b2 = (2·5 - 1·0)/2 = 10/2 = 5
  • Zeile s1: c1 = (1·4 - 2·5)/1 = (4-10)/1 = -6, dann Null
  • Zeile s0 d1 = (-6·5 - 1·0)/(-6) = (-30)/(-6) = 5

Erste Spalte: 1, 2, 1, -6, 5 → das Vorzeichen wechselt einmal (von 1 auf -6 und von -6 zurück auf 5), was auf eine Wurzel mit positivem Realteil hinweist.

Vereinfachen komplexer Polynome unter Verwendung des Kriteriums

Die Routh-Hurwitz-Methode vereinfacht die Stabilitätsanalyse komplexer Polynome auf verschiedene Weise:

  1. Vermeidet Wurzelberechnung. Für Polynome höherer Ordnung ist Factoring oder numerische Suche nach Wurzeln teuer und manchmal unzuverlässig. Das Array erfordert nur Addition, Subtraktion, Multiplikation und Division - Operationen, die leicht von Hand oder Tabellenkalkulation durchgeführt werden können.
  2. Behandelt Parameterunsicherheit. Durch Ausdrücken von Koeffizienten in Bezug auf Systemparameter (Gains, Zeitkonstanten) können Stabilitätsbedingungen in Bezug auf Ungleichheiten abgeleitet werden. Beispielsweise kann der Bereich eines Reglergewinns K, für den das Closed-Loop-System stabil bleibt, durch Anwendung des Routh-Arrays gefunden werden und es wird verlangt, dass alle Elemente der ersten Spalte positiv sind.
  3. Identifiziert die Randstabilität. Eine Null in der ersten Spalte (der ε-Fall) oder eine Reihe von Nullen zeigt Wurzeln auf der imaginären Achse an oder symmetrisch um den Ursprung.

Arbeiten mit komplexen Koeffizienten?

Während das klassische Routh-Hurwitz-Kriterium reelle Koeffizienten erfordert, gibt es viele Erweiterungen. Für Polynome mit komplexen Koeffizienten (ai ∈ C) können Modifikationen wie das modifizierte Routh-Array oder das Michailov-Kriterium verwendet werden.

Anwendungen in Kontrollsystemen und darüber hinaus

Die primäre Anwendung des Stabilitätskriteriums Routh-Hurwitz ist die Bestimmung der Stabilität von linearen zeitinvarianten Systemen (LTI), die durch Übertragungsfunktionen beschrieben werden. Der Nenner der Übertragungsfunktion ist das charakteristische Polynom, und Stabilität erfordert, dass alle Wurzeln negative reelle Teile haben.

Design von Steuerungssystemen

Während des Controller-Designs passen Ingenieure häufig Gewinne an (z. B. die proportionalen, integralen und abgeleiteten Gewinne eines PID-Controllers), um die gewünschte Leistung zu erzielen. Das Routh-Kriterium bietet eine schnelle Möglichkeit, Gain-Randbereiche zu finden - die Bandbreite der Gewinne, die das System stabil halten. Beispielsweise ist das charakteristische Polynom angesichts eines Einheits-Feedback-Systems mit Vorwärtsübertragungsfunktion G(s) = K / [(s+1)(s+2)(s+3)]32 + 11s + (6+K) = 0 Die Anwendung des Routh-Arrays ergibt die Bedingung K < 60 für Stabilität. Dies ist viel schneller als der Versuch, Wurzeln als eine Funktion von K zu berechnen.

Netzwerkanalyse

In der Elektrotechnik bestimmt das charakteristische Polynom einer RLC-Schaltung, ob transiente Schwingungen abklingen. Das Routh-Hurwitz-Kriterium hilft Ingenieuren bei der Auswahl von Bauteilwerten, um anhaltende Schwingungen (Instabilität) zu vermeiden. Ebenso wird in mechanischen Systemen die Stabilität von rückgekoppelten Robotern oder Flugzeugen mit dieser Methode routinemäßig bewertet.

Industrielle Prozesskontrolle

In der chemischen Industrie und der Erdölindustrie regeln Prozessregler Temperatur, Druck und Durchfluss. Das Routh-Array wird bei der Abstimmung von PID-Reglern verwendet, um sicherzustellen, dass die geschlossene Kreislaufanlage unter Betriebsbedingungen stabil bleibt.

Vorteile und Einschränkungen

Wie jedes Werkzeug hat das Routh-Hurwitz-Kriterium Stärken und Schwächen.

Stärken

  • Keine Wurzelberechnung erforderlich. Das Array arbeitet ausschließlich auf Koeffizienten.
  • Skalierbar. Die Methode funktioniert für Polynome jeder Ordnung, obwohl höhere Ordnungen von Hand mühsam werden (leicht automatisiert).
  • Symbolische Fähigkeit. Koeffizienten können algebraische Ausdrücke sein, die eine parametrische Stabilitätsanalyse ermöglichen.
  • Bietet die Anzahl der instabilen Wurzeln an. Jede Zeichenänderung in der ersten Spalte entspricht einer Wurzel mit positivem Realteil.

Beschränkungen

  • Nur echte Koeffizienten. Die klassische Form behandelt Polynome mit komplexen Koeffizienten nicht direkt, obwohl Erweiterungen existieren.
  • Keine Informationen über Wurzelpositionen. Das Kriterium sagt nur, wie viele instabile Wurzeln existieren, nicht deren genaue Werte. Für eine detaillierte transiente Analyse ist Wurzelort oder numerische Wurzelfindung erforderlich.
  • Empfindlichkeit gegenüber numerischen Fehlern. Wenn Koeffizienten nahezu Null sind oder Rundungsfehler enthalten, kann die Array-Konstruktion Vorzeichen-Mehrdeutigkeiten erzeugen. Die Verwendung der ε-Methode ist ein Workaround, erfordert jedoch eine sorgfältige Interpretation.
  • Behandelt keine Verzögerungssysteme. Systeme mit Zeitverzögerungen haben transzendentale charakteristische Gleichungen, keine Polynome. Das Routh-Hurwitz-Kriterium ist nicht direkt anwendbar; stattdessen werden Methoden wie Padé-Näherungen oder der Smith-Prädiktor verwendet.

Historischer Kontext und weitere Lektüre

Das Routh-Hurwitz-Kriterium ist nach Edward John Routh, einem englischen Mathematiker, der die Array-Methode 1876 veröffentlichte, und Adolf Hurwitz, einem deutschen Mathematiker, der 1895 eine ähnliche Matrixformulierung unabhängig voneinander abgeleitet hat, benannt. Ihre Arbeit baute auf früheren Beiträgen von Cauchy und Hermite auf. Das Kriterium bleibt ein Grundnahrungsmittel in modernen Kontrolllehrbüchern - siehe den Wikipedia-Artikel für eine Übersicht . Für eine tiefere mathematische Behandlung bietet die Themenseite von ScienceDirect Verweise auf Originalquellen.

Ingenieure, die fortgeschrittene Steuerung studieren, werden auch die Hurwitz-Matrix im Kontext der Lyapunov-Stabilität treffen. Die Verbindung zwischen dem Routh-Array und der Determinante der Hurwitz-Matrix ist ein klassisches Ergebnis in der linearen Algebra und der Regeltheorie.

Schritt für Schritt illustriertes Beispiel

Gehen wir durch ein vollständiges Beispiel mit einem komplexen Polynom (hoher Ordnung), um die Leistungsfähigkeit der Methode zu demonstrieren.

P(s) = 2 s5 + 4 s4 + 6 s3 + 8 s2 + 10 s + 12

Schritt 1: Ordnen Sie die ersten beiden Zeilen an

  • Zeile s5 (Leistung 5): Koeffizienten von s5, s31,1 → 2, 6, 10
  • Zeile s4 (Leistung 4): Koeffizienten von s4, s20 → 4, 8, 12

Schritt 2: Berechnen Sie die dritte Zeile (s3)

Verwendung der Determinantenformel:

  • Erstes Element: (4·6 - 2·8)/4 = (24 - 16)/4 = 8/4 = 2
  • Zweites Element: (4·10 - 2·12)/4 = (40 - 24)/4 = 16/4 = 4
  • Drittes Element: Null (da keine Koeffizienten mehr in der ersten Zeile vorhanden sind)

Zeile s3: 2, 4, 0

Schritt 3: Berechnen Sie die vierte Zeile (s2)

Mit Zeile s4 (4,8,12) und Zeile s3 (2,4,0):

  • Erstes Element: (2·8 - 4·4)/2 = (16 - 16)/2 = 0/2 = 0
  • Zweites Element: (2·12 - 4·0)/2 = (24 - 0)/2 = 12

Zeile s2: 0, 12

Hier haben wir eine Null in der ersten Spalte. Der Rest der Zeile ist ungleich Null (12). Ersetzen Sie die Null durch ε (eine kleine positive Zahl).

Schritt 4: Berechnen Sie die fünfte Zeile (s1)

Mit Zeile s3 (2,4,0) und Zeile s2 (ε,12):

  • Erstes Element: (ε·4 - 2·12)/ε = (4ε - 24)/ε = 4 - 24/ε
  • Zweites Element: (ε·0 - 2·0)/ε = 0

Zeile s1: (4 - 24/ε), 0

Als ε → 0+ dominiert der Begriff -24/ε, was dieses Element negativ macht (-∞).

Schritt 5: Berechnen Sie die sechste Zeile (s0)

Mit Zeile s2 (ε,12) und Zeile s1 (4 - 24/ε, 0):

  • Erstes Element: [(4 - 24/ε)·12 - ε·0] / (4 - 24/ε) = 12

]

Schritt 6: Analysieren Sie die erste Spalte

Die ersten Spalteneinträge lauten: Zeile s5: 2 (positiv), Zeile s4: 4 (positiv), Zeile s3: 2 (positiv), Zeile s2: ε (positiv, wie ε>0), Zeile s1: 4 - 24/ε (negativ für kleine ε), Zeile s0: 12 (positiv).

Es gibt zwei Zeichenwechsel: von positiv (Zeile s2) zu negativ (Zeile s1) und von negativ zurück zu positiv (Zeile s0 Das Polynom hat daher zwei Wurzeln mit positiven reellen Teilen.

Verifizierung: Mit numerischen Werkzeugen sind die Wurzeln von 2 s5 + ... + 12 = 0 ungefähr: -1.09, -0,50 ± 0,87i, 0,79 ± 1,04i. Tatsächlich haben die beiden komplexen konjugierten Wurzeln einen positiven Realteil (0,79), was die Routh-Hurwitz-Vorhersage bestätigt.

Alternativer Ansatz: Das Hilfspolynom für den Zero-Zeil-Fall

Wenn wir eine ganze Reihe von Nullen anstelle von nur einer Null in der ersten Spalte gefunden hätten, würden wir ein Hilfspolynom aus der obigen Zeile bilden. Diese Technik ist für Fälle wie s4 + 5s2 + 4 (die rein imaginäre Wurzeln hat ±j, ±2j). Die Nullzeile zeigt Wurzeln an, die symmetrisch zum Ursprung sind, und das Hilfspolynom (abgeleitet von der obigen Zeile) gibt die tatsächlichen Wurzeln an.

Schlussfolgerung

Das Routh-Hurwitz-Stabilitätskriterium verwandelt die komplexe Analyse von Polynomen höherer Ordnung in einen einfachen Algorithmus. Durch die Konstruktion des Routh-Arrays und die Prüfung der Zeichen der ersten Spalte können Ingenieure und Mathematiker die Systemstabilität bestimmen, ohne nach Wurzeln zu suchen. Die Methode ist für das Steuerungssystemdesign, die Schaltungsanalyse und jedes Feld, in dem die Wurzeln eines echten Polynoms das Verhalten diktieren, unerlässlich.

Die Beherrschung dieser Technik ermöglicht es, komplexe Polynome mit Zuversicht zu behandeln - sei es für akademische Studien, Industriedesign oder Forschung. In Kombination mit modernen Rechenwerkzeugen (z. B. MATLAB, Python) bleibt das Routh-Hurwitz-Kriterium ein grundlegender Baustein der dynamischen Systemanalyse. Für diejenigen, die weiter erforschen möchten, bietet das Whitepaper zu den Steuerungsprinzipien zusätzliche Beispiele und erweiterte Erweiterungen.