Wähle im Druckfenster „Als PDF speichern“, um eine Kopie auf deinem Gerät zu behalten.
1Formeln und Methoden
Stabilität: Routh-Hurwitz und Wurzelortskurve
Stabilität
alle Wurzeln des charakteristischen Polynoms in der linken Halbebene
Routh-Schema, Eintrag der s¹-Zeile eines kubischen Polynoms
b1 = (a2a1 − a3a0)/a2
Wurzeln in der rechten Halbebene
Anzahl der Vorzeichenwechsel in der ersten Spalte
Kubisches Polynom s³ + a₂s² + a₁s + a₀
stabil, wenn alle a > 0 und a₂·a₁ > a₀
K/(s(s + a)(s + b))
Kmax = ab(a + b), ω = √(ab)
Charakteristische Gleichung
1 + K·G(s) = 0
Asymptoten
σa = (Σp − Σz)/(n − m), θ = (2q + 1)·180°/(n − m)
Verzweigungs- und Einmündungspunkte
dK/ds = 0
Betragsbedingung
K = 1/|G(s)|
Abgangswinkel
θd = 180° − Σ∠(from other poles) + Σ∠(from zeros)
Methode
Schreibe das charakteristische Polynom des geschlossenen Kreises auf: den Nenner von 1 + K·G(s) = 0, ausmultipliziert.
Trage die Koeffizienten abwechselnd in die ersten beiden Zeilen des Routh-Schemas ein und bilde die übrigen Zeilen.
Zähle die Vorzeichenwechsel in der ersten Spalte: Jeder entspricht einem Pol in der rechten Halbebene. Steht K im Schema, verlange, dass jeder Eintrag der ersten Spalte positiv bleibt.
Markiere für eine Wurzelortskurve die Pole und Nullstellen, die Abschnitte der reellen Achse, die Asymptoten vom Schwerpunkt aus und die Verzweigungspunkte, an denen dK/ds = 0 gilt. Verwende K = 1/|G(s)| für die Verstärkung an einem Punkt.
Stationärer Fehler und PID-Regelung
Systemtyp
Anzahl der Pole von G(s) bei s = 0
Fehlerkonstanten
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
Sprung A
ess = A/(1 + Kp)
Rampe A·t
ess = A/Kv
Parabel A·t²/2
ess = A/Ka
PID-Regler
C(s) = Kp + Ki/s + Kd·s
Zielsystem zweiter Ordnung
s² + 2ζωns + ωn²
Überschwingen und Zeitverhalten
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)
Methode
Zähle die Pole des offenen Kreises C(s)·G(s) bei s = 0: Das ist der Systemtyp.
Berechne die Fehlerkonstante, die zur Eingangsgröße passt: Kp für einen Sprung, Kv für eine Rampe, Ka für eine Parabel.
Wende ess = A/(1 + Kp), A/Kv oder A/Ka an, oder löse die Gleichung nach der benötigten Verstärkung auf.
Schreibe bei einer Vorgabe zu Dämpfung oder Drehzahl das charakteristische Polynom des geschlossenen Kreises auf und vergleiche es mit s² + 2ζω_n·s + ω_n².
2Übungsaufgaben
1.SprungantwortMittel
Ein System hat die Übertragungsfunktion G(s) = 148/(s² + 3s + 74) und erhält einen Einheitssprung am Eingang. Bestimme das prozentuale Überschwingen %OS der Sprungantwort.
G(s) = 148/(s² + 3s + 74)
%OS = %
2.SprungantwortSchwer
Ein Reihen-RLC-Kreis mit R = 120 Ω, L = 10 mH und C = 1 µF wird bei t = 0 an einen Sprung von 10 V angeschlossen. Der Kondensator ist ungeladen, und vor dem Sprung fließt kein Strom. Bestimme das prozentuale Überschwingen %OS der Kondensatorspannung.
R = 120 ΩL = 10 mHC = 1 µFV = 10 V
%OS = %
3.Bode-Diagramme und ReservenMittel
Bestimme die Phase ∠G(jω) von G(s) = 14(s + 19)/(s(s + 16)) bei ω = 4 rad/s in Grad.
G(s) = 14(s + 19)/(s(s + 16))ω = 4 rad/s
∠G(jω) = °
4.Bode-Diagramme und ReservenSchwer
Ein Regelkreis mit Einheitsrückführung hat die Übertragungsfunktion des offenen Kreises L(s) = 76/(s(s + 2)(s + 10)). Bestimme die Phasendurchtrittsfrequenz ω180, bei der ∠L(jω) = −180° ist.
L(s) = 76/(s(s + 2)(s + 10))
ω180 = rad/s
5.Stabilität: Routh-Hurwitz und WurzelortskurveMittel
Der Regelkreis im Schaubild hat die Übertragungsfunktion des offenen Kreises L(s) = K/(s(s + 6)(s + 8)). Die Verstärkung wird erhöht, bis der geschlossene Kreis an der Stabilitätsgrenze liegt. Bestimme die Frequenz ω der Dauerschwingung.
L(s) = K/(s(s + 6)(s + 8))
ω = rad/s
6.Stabilität: Routh-Hurwitz und WurzelortskurveSchwer
Die Wurzelortskurve von L(s) = K(s + 9)/(s(s + 3)) verlässt die reelle Achse zwischen s = 0 und s = −3, läuft im Bogen herum und trifft links von der Nullstelle wieder auf die reelle Achse. Bestimme den Einmündungspunkt si.
L(s) = K(s + 9)/(s(s + 3))
si =
7.Stationärer Fehler und PID-RegelungMittel
Auf einen Regelkreis mit Einheitsrückführung und G(s) = 5(s + 5)/(s²(s + 10)) wird die parabelförmige Eingangsgröße r(t) = 0,5·t²/2 gegeben. Bestimme den stationären Fehler ess.
G(s) = 5(s + 5)/(s²(s + 10))r(t) = 0,5·t²/2
ess =
8.Stationärer Fehler und PID-RegelungSchwer
Ein PI-Regler C(s) = Kp + Ki/s regelt die Strecke G(s) = 5/(s + 2) in einem Regelkreis mit Einheitsrückführung. Der geschlossene Kreis soll ωₙ = 4 rad/s und ζ = 0,5 haben. Bestimme die Proportionalverstärkung Kp.
G(s) = 5/(s + 2)ωₙ = 4 rad/sζ = 0,5
Kp =
9.BlockschaltbilderMittel
Eine Strecke G(s) = 29/(s + 6) wird in einem Gegenkopplungskreis geregelt, dessen Sensor die Verstärkung H = 0,1 hat. Am Eingang R wird ein Sprung der Größe 24 angelegt. Bestimme den stationären Ausgang css.
G(s) = 29/(s + 6)H = 0,1R = 24
css =
10.BlockschaltbilderSchwer
Ein Proportionalregler K = 2 steuert die Strecke G(s) = 9/(s + 1) in einem Regelkreis mit Einheits-Gegenkopplung. Eine Störung D greift am Streckeneingang an (sie wird zum Reglerausgang addiert). Bei R = 0 wird ein Störsprung der Größe 5 angelegt. Bestimme den stationären Ausgang css.
K = 2G(s) = 9/(s + 1)D = 5R = 0
css =
3Lösungen
1.Sprungantwort
Antwort: %OS = 57,33 %
Sprungantwort: Ein System erster Ordnung G(s) = K/(τs + 1) beantwortet einen Sprung der Größe U mit y(t) = K·U·(1 − e−t/τ): Es erreicht nach einer Zeitkonstante τ 63,2 % des Endwerts K·U und nach 4τ 98 %. Ein System zweiter Ordnung G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) hat die Pole s = −ζωₙ ± jωd mit ωd = ωₙ√(1 − ζ²). Für 0 < ζ < 1 schwingt die Sprungantwort über: prozentuales Überschwingen %OS = 100·e−ζπ/√(1−ζ²), Überschwingzeit (Scheitelzeit) tp = π/ωd und die 2-%-Ausregelzeit ts ≈ 4/(ζωₙ). Vergleiche den Nenner mit s² + 2ζωₙs + ωₙ², um ωₙ und ζ abzulesen. Für ein Feder-Masse-Dämpfer-System gilt ωₙ = √(k/m) und ζ = c/(2√(km)); für einen Reihen-RLC-Kreis gilt ωₙ = 1/√(LC) und ζ = (R/2)·√(C/L).
Mit der Normalform vergleichen
Der Nenner muss mit s² + 2ζωₙs + ωₙ² übereinstimmen.
ωₙ² = 74 → ωₙ = 8,602 rad/s
2ζωₙ = 3 → ζ = 3/(2·8,602) = 0,1744
Ergebnis
%OS = 100·e−ζπ/√(1−ζ²) = 100·e−0,5563 = 57,33 %
Antwort
%OS = 57,33 %
2.Sprungantwort
Antwort: %OS = 9,478 %
Sprungantwort: Ein System erster Ordnung G(s) = K/(τs + 1) beantwortet einen Sprung der Größe U mit y(t) = K·U·(1 − e−t/τ): Es erreicht nach einer Zeitkonstante τ 63,2 % des Endwerts K·U und nach 4τ 98 %. Ein System zweiter Ordnung G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) hat die Pole s = −ζωₙ ± jωd mit ωd = ωₙ√(1 − ζ²). Für 0 < ζ < 1 schwingt die Sprungantwort über: prozentuales Überschwingen %OS = 100·e−ζπ/√(1−ζ²), Überschwingzeit (Scheitelzeit) tp = π/ωd und die 2-%-Ausregelzeit ts ≈ 4/(ζωₙ). Vergleiche den Nenner mit s² + 2ζωₙs + ωₙ², um ωₙ und ζ abzulesen. Für ein Feder-Masse-Dämpfer-System gilt ωₙ = √(k/m) und ζ = c/(2√(km)); für einen Reihen-RLC-Kreis gilt ωₙ = 1/√(LC) und ζ = (R/2)·√(C/L).
Eigenkreisfrequenz und Dämpfungsgrad
ωₙ = 1/√(LC) = 1/√(0,01·1 × 10⁻⁶) = 10000 rad/s
ζ = (R/2)·√(C/L) = (120/2)·√(1 × 10⁻⁶/0,01) = 0,6
Ergebnis
%OS = 100·e−ζπ/√(1−ζ²) = 9,478 %
Antwort
%OS = 9,478 %
3.Bode-Diagramme und Reserven
Antwort: ∠G(jω) = −92,15 °
Frequenzgang: Setze s = jω in die Übertragungsfunktion ein. Der Betrag in Dezibel ist 20·log₁₀|H(jω)|, und die Phase ist die Summe der Winkel der Zählerfaktoren minus die der Nennerfaktoren. Ein RC- oder RL-Filter erster Ordnung hat seine Eckfrequenz (−3 dB) bei ωc = 1/(RC) beziehungsweise R/L, und fc = ωc/(2π); ein Tiefpass liefert |H| = 1/√(1 + (ω/ωc)²) und φ = −atan(ω/ωc). Im asymptotischen (geradlinigen) Bode-Diagramm knickt jeder Pol den Betrag an seiner Eckfrequenz um 20 dB/Dekade nach unten, jede Nullstelle knickt ihn um 20 dB/Dekade nach oben, und ein Integrator 1/s ergibt bei allen Frequenzen −20 dB/Dekade. Durchtrittsfrequenz ωc: |L(jω)| = 1 (0 dB); Phasenreserve PM = 180° + ∠L(jωc). Phasendurchtrittsfrequenz ω180: ∠L(jω) = −180°; Amplitudenreserve GM = −20·log₁₀|L(jω180)| dB. Der geschlossene Regelkreis ist stabil, wenn beide Reserven positiv sind.
Frequenzgang: Setze s = jω in die Übertragungsfunktion ein. Der Betrag in Dezibel ist 20·log₁₀|H(jω)|, und die Phase ist die Summe der Winkel der Zählerfaktoren minus die der Nennerfaktoren. Ein RC- oder RL-Filter erster Ordnung hat seine Eckfrequenz (−3 dB) bei ωc = 1/(RC) beziehungsweise R/L, und fc = ωc/(2π); ein Tiefpass liefert |H| = 1/√(1 + (ω/ωc)²) und φ = −atan(ω/ωc). Im asymptotischen (geradlinigen) Bode-Diagramm knickt jeder Pol den Betrag an seiner Eckfrequenz um 20 dB/Dekade nach unten, jede Nullstelle knickt ihn um 20 dB/Dekade nach oben, und ein Integrator 1/s ergibt bei allen Frequenzen −20 dB/Dekade. Durchtrittsfrequenz ωc: |L(jω)| = 1 (0 dB); Phasenreserve PM = 180° + ∠L(jωc). Phasendurchtrittsfrequenz ω180: ∠L(jω) = −180°; Amplitudenreserve GM = −20·log₁₀|L(jω180)| dB. Der geschlossene Regelkreis ist stabil, wenn beide Reserven positiv sind.
Phasendurchtritt
(ω/2)·(ω/10) = 1 → ω180 = √(2·10) = 4,472 rad/s
Antwort
ω180 = 4,472 rad/s
5.Stabilität: Routh-Hurwitz und Wurzelortskurve
Antwort: ω = 6,928 rad/s
Routh-Hurwitz-Kriterium: Ein lineares System ist stabil, wenn alle Wurzeln seines charakteristischen Polynoms in der linken Hälfte der s-Ebene liegen. Schreibe die Koeffizienten abwechselnd in die ersten beiden Zeilen des Routh-Schemas. Jeder Eintrag einer neuen Zeile ergibt sich aus den beiden Zeilen darüber; für ein kubisches Polynom a₃s³ + a₂s² + a₁s + a₀ ist der Eintrag der s¹-Zeile b₁ = (a₂·a₁ − a₃·a₀)/a₂. Die Anzahl der Vorzeichenwechsel in der ersten Spalte entspricht der Anzahl der Wurzeln in der rechten Halbebene. Mit einer Verstärkung K in den Koeffizienten müssen alle Einträge der ersten Spalte positiv bleiben: Die Verstärkung, die den Eintrag der s¹-Zeile null macht, ist die Stabilitätsgrenze, und die s²-Zeile liefert dann die Frequenz der Dauerschwingung.
Routh-Schema
s³ + 14s² + 48s + K = 0
s³: 1, 48
s²: 14, K
s¹: (14·48 − K)/14
s⁰: K
Kritische Verstärkung
K = 14·48 = 672
Hilfsgleichung
Bei der kritischen Verstärkung ist die s¹-Zeile null; die s²-Zeile liefert die Pole auf der imaginären Achse.
14s² + 672 = 0 → s = ±j√48
ω = 6,928 rad/s
Antwort
ω = 6,928 rad/s
6.Stabilität: Routh-Hurwitz und Wurzelortskurve
Antwort: si = −16,35
Wurzelortskurve: Wenn die Verstärkung K von 0 bis ∞ wächst, wandern die Pole des geschlossenen Regelkreises, also die Wurzeln von 1 + K·G(s) = 0, von den Polen des offenen Kreises zu dessen Nullstellen oder ins Unendliche. Ein Punkt der reellen Achse gehört zur Ortskurve, wenn die Anzahl der reellen Pole und Nullstellen rechts von ihm ungerade ist. Die n − m Äste, die ins Unendliche laufen, folgen Asymptoten mit den Winkeln (2q + 1)·180°/(n − m), die im Schwerpunkt σₐ = (Σ Pole − Σ Nullstellen)/(n − m) beginnen. Äste verlassen die reelle Achse oder münden dort ein, wo dK/ds = 0 gilt. Die Winkelbedingung ∠G(s) = 180° legt die Form fest, und die Betragsbedingung K = 1/|G(s)| liefert die Verstärkung an jedem Punkt der Ortskurve.
Verstärkung entlang der reellen Achse
K(s) = −(s² + 3s)/(s + 9)
Stationäre Punkte
Differenziere mit der Quotientenregel und setze den Zähler null.
(2s + 3)(s + 9) − (s² + 3s) = s² + 18s + 27 = 0
s = −9 ± √54 = −1,652, −16,35
Einmündungspunkt
si = −16,35
Antwort
si = −16,35
Zwischen den beiden Punkten ist die Ortskurve ein Kreis um die Nullstelle mit dem Radius √54 = 7,348.
7.Stationärer Fehler und PID-Regelung
Antwort: ess = 0,2
Stationärer Fehler: In einem Regelkreis mit Einheitsrückführung und der Übertragungsfunktion des offenen Kreises G(s) ist der Fehler E(s) = R(s)/(1 + G(s)), und nach dem Endwertsatz gilt ess = lims→0 s·E(s). Der Systemtyp ist die Anzahl der Pole von G(s) bei s = 0. Die Fehlerkonstanten sind Kp = lim G(s), Kv = lim s·G(s) und Ka = lim s²·G(s), jeweils für s → 0. Ein Sprung der Höhe A liefert ess = A/(1 + Kp), eine Rampe A·t liefert A/Kv und eine Parabel A·t²/2 liefert A/Ka. Ein Kreis vom Typ 0 folgt einem Sprung mit endlichem Fehler, ein Kreis vom Typ 1 folgt ihm exakt und einer Rampe mit endlichem Fehler.
Beschleunigungskonstante
Ka = 5·5/10 = 2,5
Parabelfehler
ess = 0,5/2,5 = 0,2
Antwort
ess = 0,2
8.Stationärer Fehler und PID-Regelung
Antwort: Kp = 0,4
P-, PI- und PID-Regelung: Der Regler C(s) = Kp + Ki/s + Kd·s wirkt auf den Regelfehler. Eine proportionale Verstärkung macht den Kreis schneller und verkleinert den Fehler einer Strecke vom Typ 0, beseitigt ihn aber nie. Der Integralanteil fügt einen Pol bei s = 0 hinzu: Der Systemtyp steigt um eins, und ein Sprung wird ohne Fehler verfolgt. Der Differenzialanteil erhöht die Dämpfung. Um die Verstärkungen zu wählen, schreibe das charakteristische Polynom des geschlossenen Kreises auf und vergleiche es mit s² + 2ζωns + ωn². Dann gilt %OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd mit ωd = ωn√(1 − ζ²) und ts ≈ 4/(ζωn).
Charakteristisches Polynom
1 + C(s)G(s) = 0, multipliziert mit s(s + 2).
s² + (2 + 5Kp)s + 5Ki
Zielpolynom
s² + 2·0,5·4s + 4² = s² + 4s + 16
Koeffizienten vergleichen
2 + 5Kp = 4 → Kp = 0,4
Antwort
Kp = 0,4
Der Integralanteil macht den Kreis zu einem Kreis vom Typ 1, also wird ein Sprung ohne stationären Fehler verfolgt.
9.Blockschaltbilder
Antwort: css = 78,2
Blockschaltbild-Algebra: Blöcke in Reihe werden multipliziert (G₁·G₂), parallele Blöcke werden addiert (G₁ + G₂), und ein Regelkreis mit Vorwärtspfad G und Rückführpfad H wird bei Gegenkopplung zu T = G/(1 + G·H) und bei Mitkopplung zu G/(1 − G·H) zusammengefasst. Fasse zuerst die inneren Schleifen zusammen und arbeite dich dann nach außen. Die Pole des geschlossenen Regelkreises sind die Wurzeln der charakteristischen Gleichung 1 + G(s)H(s) = 0. Bei einem stabilen Kreis liefert ein Sprung der Größe A den stationären Ausgang A·T(0); bei Einheitsrückführung beträgt der stationäre Fehler auf einen Sprung A/(1 + Kp) mit Kp = G(0).
Blockschaltbild-Algebra: Blöcke in Reihe werden multipliziert (G₁·G₂), parallele Blöcke werden addiert (G₁ + G₂), und ein Regelkreis mit Vorwärtspfad G und Rückführpfad H wird bei Gegenkopplung zu T = G/(1 + G·H) und bei Mitkopplung zu G/(1 − G·H) zusammengefasst. Fasse zuerst die inneren Schleifen zusammen und arbeite dich dann nach außen. Die Pole des geschlossenen Regelkreises sind die Wurzeln der charakteristischen Gleichung 1 + G(s)H(s) = 0. Bei einem stabilen Kreis liefert ein Sprung der Größe A den stationären Ausgang A·T(0); bei Einheitsrückführung beträgt der stationäre Fehler auf einen Sprung A/(1 + Kp) mit Kp = G(0).
Übertragungsfunktionen
C/D = G/(1 + K·G) = 9/(s + 1 + 18)
Stationärer Zustand
css = 5·9/19 = 2,368
Antwort
css = 2,368
Ein größeres K verringert die Wirkung der Störung, aber mit einem Proportionalregler nie auf null.