Einführung in die Stabilitätsanalyse in diskreten Systemen

Stabilitätsanalyse bildet das Rückgrat der Steuerungstechnik, um sicherzustellen, dass ein System vorhersehbar reagiert und unter normalen Betriebsbedingungen begrenzt bleibt. Für zeitdiskrete Systeme &# 8212;Gemeinsam in der digitalen Steuerung, Signalverarbeitung und eingebetteten Anwendungen &# 8212;Das Stabilitätskriterium verschiebt sich von der kontinuierlichen s-Ebene linken Hälfte Ebene zum Einheitskreis in der z-Ebene. Ingenieure verlassen sich auf robuste analytische Werkzeuge, um zu bestimmen, ob alle Pole einer diskreten Übertragungsfunktion innerhalb des Einheitskreises liegen, ohne Polynome höherer Ordnung explizit zu lösen. Zwei prominente Methoden sind das Routh-Hurwitz-Kriterium, das ursprünglich für kontinuierliche Systeme entwickelt wurde, aber an diskrete Frameworks anpassbar ist, und der Stabilitätstest der Jury&# 8217;, der speziell für die zeitdiskrete Analyse entwickelt wurde.

Dieser Artikel bietet einen umfassenden Vergleich der Stabilitätstests von Routh-Hurwitz und Jury, wobei er ihre mathematischen Grundlagen, Verfahrensschritte, praktischen Anpassungen und relativen Stärken untersucht. Am Ende werden Ingenieure und Studenten gleichermaßen einen klaren Fahrplan für die Auswahl der geeigneten Technik für ihre diskreten Systemanalyseanforderungen haben.

Diskrete Systemstabilität: Das Einheitskreiskriterium

Bevor man in die spezifischen Methoden eintaucht, ist es wichtig, sich an die grundlegende Stabilitätsbedingung für zeitdiskrete Systeme zu erinnern. Ein lineares zeitinvariantes diskretes System, das durch die Übertragungsfunktion beschrieben wird H(z) = N(z) / D(#8212;die Wurzeln des charakteristischen Polynoms D(z) = 0— liegt innerhalb des Einheitskreises in der z-Ebene, d.h. |z| < 1. Dieser Zustand unterscheidet sich deutlich von kontinuierlichen Systemen, wo Stabilität erfordert, dass alle Pole negative reelle Teile haben (linke Hälftenebene).

Während Root-Finding-Algorithmen direkt Polpositionen berechnen können, werden sie für Polynome höherer Ordnung umständlich. Daher bieten algebraische Kriterien wie Routh-Hurwitz und Jury's Test effiziente, nicht-iterative Überprüfungen. Die Hauptherausforderung besteht darin, dass Jury's Test direkt auf den Koeffizienten von D(z) arbeitet, während Routh-Hurwitz typischerweise auf Polynome in der s-Domäne arbeitet. Eine Transformation ist daher notwendig, um Routh-Hurwitz auf diskrete Systeme anzuwenden.

Routh-Hurwitz-Kriterium und seine Anpassung an diskrete Systeme

Klassische Routh-Hurwitz für Continuous Systems

Das Routh-Hurwitz-Kriterium bietet notwendige und ausreichende Bedingungen für die Stabilität kontinuierlicher linearer Systeme.

a0sn + a1sn–1 + ... + a = 0

Das Verfahren konstruiert aus den Koeffizienten ein Routh-Array, das System ist stabil, wenn alle Elemente in der ersten Spalte des Arrays positiv und ungleich Null sind, ein Vorzeichenwechsel zeigt Wurzeln in der rechten Hälfte der Ebene und damit Instabilität an. Dieser Ansatz vermeidet eine explizite Wurzelberechnung und ist für Polynome bis zu moderaten Ordnungen einfach.

Bilineare Transformation: Umwandlung der z-Ebene in die s-Ebene

Um Routh-Hurwitz auf diskrete Systeme anzuwenden, verwenden Ingenieure die bilineare Transformation (auch Tustin-Transformation genannt). Diese Abbildung bezieht die Variable z auf eine neue komplexe Variable w (oft als s bezeichnet, um dies nach der Transformation zu erleichtern):

z = (1 + wT/2) / (1 – wT/2)

oder, für normalisierte Abtastperiode T = 2, die einfachere Form z = (1 + w) / (1 – w) Unter dieser bilinearen Karte bildet das Innere des Einheitskreises in der z-Ebene die linke Hälfte der Ebene in der w-Ebene ab. Folglich kann das transformierte Polynom in w unter Verwendung des Standard-Routh-Hurwitz-Arrays analysiert werden.

Schritte zum Anwenden von Routh-Hurwitz auf diskrete Systeme

  1. Erhalten Sie das diskrete charakteristische Polynom D(z) aus der Übertragungsfunktion des Systems.
  2. Ersetzen Sie die bilineare Transformation in ]D(z) , um ein Polynom in w zu erhalten. Dieser Schritt kann algebraische Expansion und Vereinfachung beinhalten, insbesondere für Systeme höherer Ordnung.
  3. Wenden Sie das Routh-Hurwitz-Kriterium auf das transformierte Polynom in w an.
  4. Stabilitätsschluss: Wenn alle Einträge in der ersten Spalte positiv sind, ist das diskrete System stabil (alle Pole innerhalb des Einheitskreises).

Während dieser Ansatz funktioniert, führt er zu zusätzlicher algebraischer Komplexität. Die bilineare Transformation kann die polynomielle Ordnung erheblich erhöhen, wenn sie nicht sorgfältig gehandhabt wird. Darüber hinaus kann die Transformation selbst zu einer Quelle numerischer Fehler werden, insbesondere für Systeme höherer Ordnung oder wenn die Abtastperioden nicht richtig normalisiert werden. Trotz dieser Nachteile bleibt das Routh-Hurwitz-Verfahren eine praktikable Option für Ingenieure, die bereits mit ihrer Anwendung für kontinuierliche Systeme vertraut sind und ein einziges einheitliches Framework bevorzugen.

Grenzen der Routh-Hurwitz-Adaption

  • Algebraische Belastung: Die bilineare Transformation erweitert das Polynom, führt möglicherweise Rundungsfehler ein und macht manuelle Berechnungen langweilig.
  • Numerische Empfindlichkeit: Für Polynome höherer Ordnung können die Koeffizientengrößen stark variieren, was zu schlecht konditionierten Arrays führt.
  • Sonderfälle: Die Methode erfordert die Handhabung von Hilfspolynomen für Wurzeln auf der imaginären Achse (in der w-Ebene], die Polen auf dem Einheitskreis entsprechen; eine Randstabilitätsbedingung, die oft unerwünscht ist.
  • Ineffizienz für rein diskrete Systeme: Da die Methode nicht nativ für die diskrete Domäne ist, empfehlen viele Kontrolllehrbücher den Test der Jury als direktere und weniger fehleranfällige Alternative.

Dennoch wird das Routh-Hurwitz-Kriterium manchmal in Kontrollkursen als Brücke zwischen kontinuierlicher und diskreter Analyse gelehrt und kann mit symbolischen Berechnungswerkzeugen angewendet werden, um algebraische Schwierigkeiten zu mildern.

Stabilitätstest der Jury: Konzipiert für die diskrete Domäne

Ursprünge und Prinzipien

Der Stabilitätstest der Jury, der von Eliahu I. Jury in den 1960er Jahren eingeführt wurde, ist das diskrete Gegenstück zum Routh-Hurwitz-Kriterium. Er untersucht direkt die Koeffizienten des charakteristischen Polynoms ]D(z), ohne Domänentransformationen zu erfordern. Anstatt das Vorzeichen von Einträgen in der ersten Spalte zu überprüfen, konstruiert der Test der Jury eine Tabelle (die Jury-Tabelle) aus den Polynomkoeffizienten und erzwingt eine Reihe von Bedingungen, die garantieren, dass alle Wurzeln eine Größe von weniger als eins haben.

Betrachten Sie das monische charakteristische Polynom der Ordnung n:

D(z) = a0zn + a1zn–1 + … + an–1z + an = 0, with a0 > 0.

Der Test der Jury bietet ein systematisches Verfahren, um festzustellen, ob alle Wurzeln innerhalb des Einheitskreises liegen, indem eine Tabelle mit n + 1 Zeilen gebildet wird.

Verfahren für die Erstellung des Jury-Tisches

  1. Hier sind die Koeffizienten in der ersten Zeile (Zeile 0) in natürlicher Reihenfolge aufzulisten: a0, a1, ..., an.
  2. Bilde die zweite Zeile (Zeile 1), indem du die Reihenfolge der Koeffizienten in Zeile 0 umkehrst: an, an–1, ..., a0
  3. ] für nachfolgende Zeilen (k = 2, 3, ...] ]Jury-Determinantenk, i wird typischerweise durch Folgendes gegeben:
    b* + 1]a]n

    ]

    ]

    , wobei die Indizes auf der Basis des aktuellen Zeilenpaares angepasst werden. Die genaue Formel variiert mit der Implementierung; gängige Lehrbücher stellen einen rekursiven Algorithmus dar, der die polynomielle Ordnung nach jedem Zeilenpaar um eins reduziert.
  4. Weiter, bis nur ein Nicht-Null-Element in der ersten Spalte verbleibt, oder bis die Tabelle zusammenbricht.
  5. Überprüfen Sie die Bedingungen: Für Stabilität muss Folgendes gelten:
    • D(1) > 0
    • (–1)n D(–1) >
    • Für die Jury-Tabelle muss das erste Element jeder ungeradzahligen Zeile (Zeilen 1, 3, 5, ...) größer als Null sein.

Diese Bedingungen zusammen sind notwendig und ausreichend, damit alle Wurzeln erfüllen |z| < 1. In der Praxis wird der Test in Computeralgebra-Systemen implementiert, aber der manuelle Prozess ist für Bestellungen bis zu etwa 4 oder 5 überschaubar.

Beispiel: System zweiter Ordnung

Betrachten wir ein diskretes System zweiter Ordnung mit charakteristischem Polynom D(z) = z2 + a1z + a = 0. Bewerbung des Jurytests:

  • Zeile 0: 1, a1, a2
  • Zeile 1: a2, a1, 1
  • Zeile 2: Die nächste Zeile hat Einträge b0 = 1 – a221aa1a2

Stabilitätsbedingungen für zweitrangige:

  1. D(1) = 1 + a1 + a2 >
  2. (–1)2D(–1) = 1 – a1 + a2 > 0
  3. | a2 | < 1

Das sind die bekannten Bedingungen für ein stabiles, zweitklassiges, diskretes System. Beachten Sie, dass der Test der Jury direkt diese einfachen Ungleichheiten ohne Transformation ergibt.

Vorteile des Jury’s Tests

  • Direkte Anwendung: Funktioniert auf dem z-Domänenpolynom; keine Zuordnung erforderlich.
  • Füllen Sie die Bedingungen ab: Bietet sowohl notwendige als auch ausreichende Stabilitätskriterien.
  • Behandelt Sonderfälle: Beinhaltet die Überprüfungen auf Wurzeln im Einheitskreis (marginale Stabilität) durch die D(1) und D(–1) Bedingungen.
  • Geeignet für die digitale Implementierung: Die Tabellenkonstruktion ist algorithmisch und leicht zu programmieren.

Head-to-Head-Vergleich von Routh-Hurwitz und Jury-Stabilitätstest

AspectRouth-Hurwitz (Adapted)Jury’s Stability Test
DomainContinuous (s-plane) originally; adapted via bilinear transform to w-planeDiscrete (z-plane) directly
Stability regionLeft-half plane (after transformation, maps to inside unit circle)Inside unit circle (|z|<1)
ProcedureConstruct Routh array from w-polynomial; check sign of first columnConstruct Jury table from z-polynomial; check series of inequalities
ComplexityModerate after transformation; transformation adds extra algebraic stepsLow to moderate; table construction straightforward for low to medium orders
Numerical stabilityCan suffer from coefficient scaling after bilinear transformGenerally good; determinants may cause swelling for high orders
Handling of marginal stabilityAuxiliary polynomial required; roots on imaginary axis in w-planeDirect checks through D(1) and D(-1) conditions
Educational accessibilityStudents often learn continuous version first; adaptation extends existing knowledgeRequires understanding of unit circle concept; but method is self-contained
Software implementationEasily coded via array construction; bilinear transform may require symbolic toolboxEasily coded via recursive table; native to z-domain
Typical use caseWhen system is already in s-domain or hybrid continuous/discreteWhen characteristic polynomial is directly from a discrete transfer function

Detaillierte Vergleichspunkte

1. Anwendbarkeit und Domain Fit

Der Jury-Test ist die natürliche Wahl für reine diskrete Systeme, insbesondere in digitalen Steuerungsanwendungen, bei denen der Controller-Output und der Anlageninput gesampledert werden. Routh-Hurwitz erfordert eine Transformation, die Komplexität hinzufügt und die physikalische Interpretation von Polen relativ zum Einheitskreis verschleiern kann. Viele moderne Lehrbücher und Kurse zur digitalen Steuerung stellen den Test der Jury als Standardwerkzeug dar.

2. Algebraische Komplexität

Für Ordnungen bis zu 4 oder 5 sind beide Methoden von Hand handhabbar. Allerdings vermeidet der Test der Jury die bilineare Expansion, die Begriffe mit gemischten Potenzen und großen Koeffizienten erzeugen kann. Zum Beispiel ergibt ein dritter Ordnung D(z) transformiert zu w selbst ein Polynom dritter Ordnung, aber die Koeffizienten sind Kombinationen der ursprünglichen; für höhere Ordnungen wird die Transformation mühsamer. Die Jury's Tischkonstruktion arbeitet direkt auf den ursprünglichen Koeffizienten, während sie Determinanten beinhaltet.

3. Numerische Robustheit

Routh-Hurwitz in der w-Domäne kann unter einer schlechten Konditionierung leiden, wenn die bilineare Transformation Koeffizienten mit sehr unterschiedlichen Größen erzeugt. Der Test der Jury steht auch vor numerischen Problemen für Polynome höherer Ordnung (z. B. Ordnung > 10) aufgrund wiederholter Subtraktionen und Determinantenberechnungen, aber diese sind weniger streng als die durch die Transformation eingeführten.

4. Fähigkeit, die Anzahl der instabilen Pole zu bestimmen

Routh-Hurwitz kann aufgrund der Vorzeichenänderungen des Routh-Arrays direkt die Anzahl der Pole in der rechten Hälfte der Ebene angeben, was der Anzahl der Pole außerhalb des Einheitskreises entspricht. Der Jury-Test liefert in seiner Grundform eine Ja/Nein-Stabilitätsantwort.

5. Umgang mit speziellen Polynomkoeffizienten

Wenn ein diskretes System Nullkoeffizienten (fehlende Befugnisse) hat, müssen beide Methoden vorsichtig sein. Routh-Hurwitz behandelt Nulleinträge im Array mit kleinen Epsilon-Näherungen. Der Test der Jury muss ebenfalls Nullen während der Tischkonstruktion verwalten, aber der Algorithmus bleibt gültig, solange der führende Koeffizient a0 ungleich Null ist. Wenn eine Zeile in der Jury-Tabelle vollständig Null wird, hat das Polynom symmetrische Wurzeln, die auf Randstabilität oder Antisymmetrie hinweisen.

Praktische Anwendungen und wann jede Methode zu wählen ist

Verwenden Sie den Jury-Test, wenn

  • Sie arbeiten ausschließlich in der diskreten Domäne (z-Ebene).
  • Sie möchten einen einstufigen Test ohne Domänentransformationen.
  • Sie analysieren einen digitalen Controller oder ein Datensystem mit einer bekannten Abtastperiode.
  • Sie benötigen klare Bedingungen für Randstabilität (Wurzeln auf dem Einheitskreis).
  • Die polynomielle Ordnung ist moderat (≤ 8 für manuelle Prüfungen; höher für Software).

Verwenden Sie Routh-Hurwitz (über Bilineare Transformation)

  • Ihr Steuerungssystemdesign befindet sich bereits im kontinuierlichen Bereich und Sie werden später diskretisiert.
  • Sie bevorzugen ein einheitliches Analyse-Framework für kontinuierliche und diskrete Systeme.
  • Sie haben Software (z.B. MATLAB, Python mit SymPy), um die bilineare Transformation automatisch zu handhaben.
  • Sie müssen die genaue Anzahl der instabilen Pole zählen (über Vorzeichenänderungen im Routh-Array).
  • Sie lehren die Verbindung zwischen kontinuierlichen und diskreten Stabilitätskriterien.

Real-World Beispiele

Beispiel: Digitaler Motordrehzahlregler. Ein Mikrocontroller implementiert einen PI-Controller für einen Gleichstrommotor. Das geschlossene Kennfeld wird in der z-Domäne abgeleitet. Ein Ingenieur würde den Test der Jury verwenden, um die Stabilität auf verschiedene Abtastraten und proportionale Gewinne zu überprüfen. Der Test zeigt schnell, ob sich die Pole außerhalb des Einheitskreises bewegen, wenn sich der Gewinn erhöht.

Beispiel: Adaptives Filterdesign. In der Signalverarbeitung benötigen IIR-Filter Pole innerhalb des Einheitskreises, um Selbstschwingungen zu vermeiden. Der Test der Jury ist üblicherweise in Filterdesignbibliotheken eingebettet, um die Stabilität nach der Koeffizientenquantisierung zu validieren.

Beispiel: Hybrides kontinuierlich-diskretes System. Eine kontinuierliche Anlage, die von einem digitalen Kompensator gesteuert wird, erfordert eine Diskretisierung über einen Haltevorgang nullter Ordnung. Die charakteristische Gesamtgleichung in der z-Domäne kann entweder direkt durch den Test der Jury analysiert werden oder durch die bilineare Transformation, um Routh-Hurwitz anzuwenden. Viele Steuerungsingenieure bevorzugen den Test der Jury, weil sie in der z-Domäne bleiben und zusätzliche Transformationen vermeiden können.

Softwareimplementierung und Automatisierung

Beide Methoden sind in Steuerungssystem-Toolboxen verfügbar, zum Beispiel:

  • MATLAB: Die Funktion (oder benutzerdefinierte Skripte) kann das Routh-Array erstellen, während in der Robotics Toolbox verfügbar ist oder programmiert werden kann.
  • Python (SciPy und Kontrollbibliotheken): Die Bibliothek enthält für den Test der Jury.
  • Symbolische Werkzeuge: SymPy und Mathematica ermöglichen algebraische Manipulation für die bilineare Transformation und die nachfolgende Array-Generierung.

Bei der Verwendung von Software erfolgt der Transformationsschritt für Routh-Hurwitz automatisiert, wodurch manuelle Fehler reduziert werden. Für das Rapid Prototyping von diskreten Systemen ist der Test der Jury jedoch die direktere Wahl.

Einschränkungen und Caveats

Beide Verfahren sind kein Allheilmittel. Beide werden aufgrund der Akkumulation von numerischen Fehlern für Polynome höherer Ordnung (Ordnung > 10) unhandlich. In solchen Fällen ist die Berechnung von Wurzelfindungen (z. B. Eigenwertverfahren) oft zuverlässiger. Beide Tests gehen außerdem davon aus, dass das charakteristische Polynom genau bekannt ist; Modellierungsunsicherheiten und Parameterschwankungen erfordern eine robuste Stabilitätsanalyse über diese algebraischen Kriterien hinaus.

Ein weiterer wichtiger Vorbehalt: Die Routh-Hurwitz-Adaption über bilineare Transformation gilt nur für lineare zeitinvariante Systeme. Sie erstreckt sich nicht direkt auf zeitvariante oder nichtlineare diskrete Systeme. Der Test der Jury gilt ebenfalls nur für lineare diskrete Modelle.

Schließlich kann es zu Fehlinterpretationen der Bedingungen kommen. Für den Test der Jury sind die Bedingungen D(1) > und (–1)nD(–1) > 0 notwendig, aber nicht ausreichend allein; sie müssen mit den Tabellenbedingungen kombiniert werden. Die Schüler hören manchmal nach diesen beiden einfachen Überprüfungen auf, was zu falschen Stabilitätsschlussfolgerungen für übergeordnete Systeme führen kann.

Schlussfolgerung

Die Wahl zwischen Routh-Hurwitz und Jury's Stabilitätstest für diskrete Systeme hängt vom Kontext, der Vertrautheit und den Rechenressourcen ab. Jury's Test ist die direktere, zweckgerichtete Methode für zeitdiskrete Polynome, die klare Bedingungen bietet und die bilineare Transformation vermeidet. Das angepasste Routh-Hurwitz-Kriterium bietet eine Brücke für Ingenieure, die in kontinuierlichen Systemen ausgebildet sind und bei der Integration kontinuierlicher und diskreter Analysen nützlich sein können. Beide Methoden bleiben wertvolle Lehrmittel, die das grundlegende Konzept der Pollage in Bezug auf den Einheitskreis verstärken.

Für die meisten modernen digitalen Steuerung Design und diskrete Systemanalyse, Jury's Stabilitätstest ist der empfohlene Ansatz. Allerdings Kenntnisse in beiden Techniken stellt sicher, dass der Ingenieur auf verschiedene Problemstellungen anpassen und nutzen verfügbare Software, um Ergebnisse zu überprüfen. Verständnis der Stärken und Schwächen der einzelnen Methoden befähigt Designer, fundierte Entscheidungen zu treffen, was letztlich zu robuster und zuverlässigere Steuerungssysteme.

Weiteres Lesen und externe Ressourcen