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.Laplace-TransformationenMittel
F(s) = (3s + 20)/(s² + 5s + 4). Ihre Rücktransformierte hat die Form f(t) = A·e−t + B·e−4t. Bestimme den Koeffizienten A mit der Partialbruchzerlegung.
F(s) = (3s + 20)/(s² + 5s + 4)
A =
2.Laplace-TransformationenSchwer
Ein System gehorcht y″ + 8y′ + 12y = 72·u(t). Es startet in Ruhe (y(0) = 0, y′(0) = 0), und der Eingang ist ein Sprung u(t) = 5 für t ≥ 0. Bestimme mit der Laplace-Transformation den Ausgang y(0,5).
Eine Strecke G(s) = 18/(s + 9) wird in einem Gegenkopplungskreis geregelt, dessen Sensor die Verstärkung H = 2 hat. Am Eingang R wird ein Sprung der Größe 12 angelegt. Bestimme den stationären Ausgang css.
G(s) = 18/(s + 9)H = 2R = 12
css =
4.BlockschaltbilderSchwer
Ein Proportionalregler K steuert die Strecke G(s) = 8/(s + 5) in einem Regelkreis mit Einheits-Gegenkopplung. Ein Störsprung der Größe 1 greift am Streckeneingang an, während R = 0 ist. Bestimme die kleinste Verstärkung K, die den stationären Ausgang bei höchstens 0,32 hält.
G(s) = 8/(s + 5)D = 1css ≤ 0,32
K =
5.SprungantwortMittel
Ein System hat die Übertragungsfunktion G(s) = 28/(s² + 4s + 56) und erhält einen Einheitssprung am Eingang. Bestimme die 2-%-Ausregelzeit ts der Sprungantwort mit ts ≈ 4/(ζωₙ).
G(s) = 28/(s² + 4s + 56)
ts = s
6.SprungantwortSchwer
Ein Regelkreis mit Einheits-Gegenkopplung hat G(s) = 295/(s(s + 20)). Bestimme die 2-%-Ausregelzeit ts der Sprungantwort des geschlossenen Kreises mit ts ≈ 4/(ζωₙ).
G(s) = 295/(s(s + 20))
ts = s
7.Stabilität: Routh-Hurwitz und WurzelortskurveMittel
Ein Regelkreis mit Einheitsrückführung hat L(s) = K/((s + 2)(s + 3)(s + 4)) mit K > 0. Bestimme die größte Verstärkung Kmax, für die der geschlossene Kreis stabil ist.
L(s) = K/((s + 2)(s + 3)(s + 4))
Kmax =
8.Stabilität: Routh-Hurwitz und WurzelortskurveSchwer
Bestimme für den Regelkreis L(s) = K/(s(s + 3)(s + 6)) die Verstärkung K, die einen Pol des geschlossenen Kreises auf s = −1,5 legt.
L(s) = K/(s(s + 3)(s + 6))s = −1,5
K =
9.Stationärer Fehler und PID-RegelungMittel
Auf einen Regelkreis mit Einheitsrückführung und G(s) = 10(s + 3)/(s²(s + 10)) wird die parabelförmige Eingangsgröße r(t) = 2·t²/2 gegeben. Bestimme den stationären Fehler ess.
G(s) = 10(s + 3)/(s²(s + 10))r(t) = 2·t²/2
ess =
10.Stationärer Fehler und PID-RegelungSchwer
Ein Proportionalregler K = 20 regelt die Strecke G(s) = 4/(s(s + 6)) in einem Regelkreis mit Einheitsrückführung. Bestimme die 2-%-Ausregelzeit ts der Sprungantwort mit ts ≈ 4/(ζωₙ).
G(s) = 4/(s(s + 6))K = 20
ts = s
11.Bode-Diagramme und ReservenMittel
Bestimme den exakten Betrag |G(jω)| von G(s) = 75/((s + 1)(s + 15)) bei ω = 4 rad/s in dB.
G(s) = 75/((s + 1)(s + 15))ω = 4 rad/s
|G(jω)| = dB
12.Bode-Diagramme und ReservenSchwer
Ein Regelkreis mit Einheitsrückführung hat die Übertragungsfunktion des offenen Kreises L(s) = 450/(s(s + 10)). Bestimme die Durchtrittsfrequenz ωc, bei der |L(jω)| = 1 ist.
L(s) = 450/(s(s + 10))
ωc = rad/s
3Lösungen
1.Laplace-Transformationen
Antwort: A = 5,667
Laplace-Transformation: F(s) = ∫₀∞ f(t)·e−st dt macht aus Ableitungen Multiplikationen mit s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Standardpaare: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), und eine Verschiebung s → s + σ multipliziert die Zeitfunktion mit e−σt. Um eine rationale Funktion F(s) zurückzutransformieren, zerlege sie in Partialbrüche: Das Residuum eines einfachen Pols p ist A = [(s − p)·F(s)] bei s = p. Endwertsatz: f(∞) = lims→0 s·F(s), gültig nur, wenn alle Pole von s·F(s) in der linken Halbebene liegen. Anfangswertsatz: f(0⁺) = lims→∞ s·F(s).
Nenner in Faktoren zerlegen
s² + 5s + 4 = (s + 1)(s + 4)
F(s) = A/(s + 1) + B/(s + 4)
Zuhaltemethode
Multipliziere mit (s + 1) und setze s = −1: Der B-Term verschwindet.
A = (3s + 20)/(s + 4) |s = −1 = 17/3 = 5,667
Antwort
A = 5,667
(B = −2,667, also f(t) = 5,667·e−t − 2,667·e−4t)
2.Laplace-Transformationen
Antwort: y(0.5) = 14,19
Laplace-Transformation: F(s) = ∫₀∞ f(t)·e−st dt macht aus Ableitungen Multiplikationen mit s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Standardpaare: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), und eine Verschiebung s → s + σ multipliziert die Zeitfunktion mit e−σt. Um eine rationale Funktion F(s) zurückzutransformieren, zerlege sie in Partialbrüche: Das Residuum eines einfachen Pols p ist A = [(s − p)·F(s)] bei s = p. Endwertsatz: f(∞) = lims→0 s·F(s), gültig nur, wenn alle Pole von s·F(s) in der linken Halbebene liegen. Anfangswertsatz: f(0⁺) = lims→∞ s·F(s).
Gleichung transformieren
Y(s) = 360/(s(s + 2)(s + 6))
Partialbruchzerlegung
A = 360/(2·6) = 30
B = 360/((−2)·4) = −45
C = 360/((−6)·(−4)) = 15
Ausfuhr
y(t) = 30 − 45·e−2t + 15·e−6t
y(0,5) = 14,19
Antwort
y(0,5) = 14,19
3.Blockschaltbilder
Antwort: css = 4,8
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).
Störantwort
css = D·G(0)/(1 + K·G(0)) = 1·8/(5 + 8K)
Nach K auflösen
5 + 8K = 1·8/0,32 = 25
K = (25 − 5)/8 = 2,5
Antwort
K = 2,5
5.Sprungantwort
Antwort: ts = 2 s
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.
ωₙ² = 56 → ωₙ = 7,483 rad/s
2ζωₙ = 4 → ζ = 4/(2·7,483) = 0,2673
Ergebnis
ζωₙ = 2 → ts ≈ 4/2 = 2 s
Antwort
ts = 2 s
6.Sprungantwort
Antwort: ts = 0,4 s
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).
Für diesen Kreis ist ζωₙ = a/2, unabhängig von K: Ein größeres K erhöht das Überschwingen, verkürzt aber nicht die Ausregelzeit.
7.Stabilität: Routh-Hurwitz und Wurzelortskurve
Antwort: Kmax = 210
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.
Charakteristisches Polynom
s³ + 9s² + 26s + 24 + K = 0
Routh-Schema
s³: 1, 26
s²: 9, 24 + K
s¹: (9·26 − 24 − K)/9
s⁰: 24 + K
Bedingung
Der Eintrag der s¹-Zeile muss positiv bleiben.
K < 9·26 − 24 = 210
Antwort
Kmax = 210
An der Grenze schwingt der Kreis mit ω = √26 = 5,099 rad/s.
8.Stabilität: Routh-Hurwitz und Wurzelortskurve
Antwort: K = 10,13
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.
Betragsbedingung
K = |s|·|s + 3|·|s + 6|
Am Punkt auswerten
K = 1,5·1,5·4,5 = 10,13
Antwort
K = 10,13
Probe: s(s + 3)(s + 6) + K ist null bei s = −1,5.
9.Stationärer Fehler und PID-Regelung
Antwort: ess = 0,6667
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 = 10·3/10 = 3
Parabelfehler
ess = 2/3 = 0,6667
Antwort
ess = 0,6667
10.Stationärer Fehler und PID-Regelung
Antwort: ts = 1,333 s
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
s(s + 6) + 20·4 = s² + 6s + 80
ωₙ und ζ
ωₙ = √80 = 8,944 rad/s
ζ = 6/(2·8,944) = 0,3354
Ausregelzeit
ζωₙ = 6/2 = 3
ts ≈ 4/3 = 1,333 s
Antwort
ts = 1,333 s
Die Ausregelzeit hängt hier nicht von K ab: Eine größere Verstärkung lässt die Antwort nur schneller schwingen.
11.Bode-Diagramme und Reserven
Antwort: |G(jω)| = 1,377 dB
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.