Übertragungsfunktion Rechner

Berechnen Sie die Übertragungsfunktion Ihres Systems mit präzisen mathematischen Methoden. Geben Sie die Systemparameter ein und erhalten Sie sofortige Ergebnisse mit visualisierter Frequenzgangdarstellung.

Umfassender Leitfaden zur Übertragungsfunktion und deren Berechnung

Die Übertragungsfunktion ist ein fundamentales Konzept in der Regelungstechnik und Signalverarbeitung, das das Verhalten linearer zeitinvarianter Systeme (LTI-Systeme) im Frequenzbereich beschreibt. Dieser Leitfaden erklärt die theoretischen Grundlagen, praktische Anwendungen und die mathematische Herleitung von Übertragungsfunktionen für Systeme 1., 2. und höherer Ordnung.

1. Grundlagen der Übertragungsfunktion

Eine Übertragungsfunktion G(s) definiert das Verhältnis der Laplace-transformierten Ausgangsgröße Y(s) zur Laplace-transformierten Eingangsgröße U(s) eines Systems:

G(s) = Y(s)/U(s) = K · N(s)/D(s)

Dabei gilt:

• K: Verstärkungsfaktor (statische Verstärkung)
• N(s): Zählerpolynom (Nullstellen)
• D(s): Nennerpolynom (Polstellen)

Die Übertragungsfunktion ermöglicht:

1. Analyse des Frequenzverhaltens (Bode-Diagramm)
2. Stabilitätsuntersuchungen (Pol-Nullstellen-Verteilung)
3. Entwurf von Reglern (PID, Lead-Lag-Kompensatoren)
4. Vorhersage des Zeitverhaltens (Sprungantwort, Impulsantwort)

2. Übertragungsfunktionen nach Systemordnung

2.1 Systeme 1. Ordnung (PT1-Glied)

Die Standardform lautet:

G(s) = K / (1 + τs)

Mit:

• K: Verstärkung
• τ: Zeitkonstante [s]

Eigenschaften:

• Eckfrequenz bei ω = 1/τ (Amplitudenabfall um 3 dB)
• Phasenverschiebung von -45° bei Eckfrequenz
• Asymptotischer Abfall mit -20 dB/Dekade

2.2 Systeme 2. Ordnung (PT2-Glied)

Die Standardform lautet:

G(s) = Kωₙ² / (s² + 2ζωₙs + ωₙ²)

Mit:

• K: Verstärkung
• ζ: Dämpfungsgrad (0 < ζ < 1 für schwache Dämpfung)
• ωₙ: Ungedämpfte Eigenfrequenz [rad/s]

Dämpfungsfälle:

Dämpfungsgrad (ζ)	Systemverhalten	Polstellen	Resonanzüberhöhung
ζ = 0	Ungedämpft	Imaginär (±jωₙ)	∞ (unendlich)
0 < ζ < 1	Schwingungsfähig	Komplex konjugiert	Abhängig von ζ
ζ = 1	Kritisch gedämpft	Doppelte reelle Polstelle	Keine
ζ > 1	Kriechfall	Zwei reelle Polstellen	Keine

Resonanzfrequenz: ω_r = ωₙ√(1 – 2ζ²) für ζ < 0.707

Resonanzüberhöhung: M_r = 1/(2ζ√(1-ζ²)) [dB]

2.3 Systeme höherer Ordnung

Systeme ≥3. Ordnung lassen sich oft in Partialbrüche zerlegen:

G(s) = K / [(1 + τ₁s)(1 + τ₂s)…(1 + τₙs)]

Die Gesamtübertragungsfunktion ergibt sich aus der Multiplikation der Teilübertragungsfunktionen.

3. Praktische Anwendungsbeispiele

3.1 Elektrische Schaltungen

RC-Tiefpass (1. Ordnung):

G(s) = 1 / (1 + RCs) → τ = RC

RLC-Schwingkreis (2. Ordnung):

G(s) = (1/LC) / (s² + (R/L)s + 1/LC)

Mit: ζ = R/(2√(L/C)), ωₙ = 1/√(LC)

3.2 Mechanische Systeme

Feder-Masse-Dämpfer (2. Ordnung):

G(s) = (1/m) / (s² + (d/m)s + k/m)

Mit: ωₙ = √(k/m), ζ = d/(2√(km))

4. Analyseverfahren im Frequenzbereich

4.1 Bode-Diagramm

Das Bode-Diagramm besteht aus:

• Amplitudengang: 20·log|G(jω)| [dB]
• Phasengang: ∠G(jω) [°]

Asymptotische Näherungen:

Systemtyp	Amplitudenasymptoten	Phasenverlauf
PT1-Glied	0 dB bis ω=1/τ, dann -20 dB/Dekade	0° bis -90° bei 10·ω=1/τ
PT2-Glied (ζ < 1)	0 dB bis ω=ωₙ, dann -40 dB/Dekade	0° bis -180°
IT1-Glied (I-Anteil)	-20 dB/Dekade durchgehend	-90° konstant

4.2 Nyquist-Diagramm

Zeigt G(jω) in der komplexen Ebene (Realteil vs. Imaginärteil). Kritisch für Stabilitätsanalyse nach dem Nyquist-Kriterium:

• Umfassung des kritischen Punktes (-1, j0) bestimmt Stabilität
• Phasenreserve und Amplitudenreserve ablesbar

5. Praktische Tipps für die Berechnung

5.1 Normalisierung

Normieren Sie die Übertragungsfunktion auf die Form:

G(s) = K · (1 + s/ω₁)(1 + s/ω₂)… / (1 + s/ωₐ)(1 + s/ω_b)…

5.2 Dominante Pole

Für Systeme höherer Ordnung:

• Identifizieren Sie die Pole mit der kleinsten Zeitkonstante (dominante Pole)
• Vernachlässigen Sie Pole, die mindestens 5-mal weiter links in der s-Ebene liegen
• Nutzen Sie die reduzierte Übertragungsfunktion für vereinfachte Analysen

5.3 Numerische Genauigkeit

Bei der Implementierung in Software:

• Verwenden Sie doppelte Genauigkeit (double) für Frequenzberechnungen
• Begrenzen Sie den Frequenzbereich auf 5 Dekaden (z. B. 0.01·ωₙ bis 100·ωₙ)
• Nutzen Sie logarithmische Skalierung für Achsen in Bode-Diagrammen

Autoritäre Quellen zu Übertragungsfunktionen:

Für vertiefende Informationen empfehlen wir folgende akademische Ressourcen:

University of Michigan: Control Tutorials for MATLAB – System Modeling MIT: Pole-Zero Plots and Transfer Functions (PDF) NASA Technical Report: Transfer Function Analysis for Control Systems

6. Häufige Fehler und deren Vermeidung

6.1 Falsche Systemordnung

Problem: Annahme eines PT1-Verhaltens bei tatsächlich vorliegendem PT2-System.

Lösung: Führen Sie Sprungantwort-Tests durch und analysieren Sie das Überschwingen.

6.2 Vernachlässigung von Nullstellen

Problem: Nicht-minimumphasige Systeme (Nullstellen in der rechten s-Halbebene) werden wie minimumphasige behandelt.

Lösung: Prüfen Sie immer die Pol-Nullstellen-Verteilung mit MATLAB oder Python (scipy.signal).

6.3 Einheitenfehler

Problem: Vermischung von Hz und rad/s bei der Eingabe von Frequenzen.

Lösung: Konvertieren Sie konsistent mit ω = 2πf.

6.4 Falsche Dämpfungsannahmen

Problem: Annahme von ζ=0.7 für “gute Dämpfung” ohne experimentelle Validierung.

Lösung: Bestimmen Sie ζ aus gemessener Resonanzüberhöhung oder logarithmischem Dekrement.

7. Erweiterte Anwendungen

7.1 Reglerentwurf mit Übertragungsfunktionen

Nutzen Sie die Übertragungsfunktion für:

• PID-Regler: G_R(s) = K_p + K_i/s + K_ds
• Lead-Lag-Kompensatoren: G_R(s) = K·(1 + sT₁)/(1 + sT₂)
• Notch-Filter: Unterdrückung spezifischer Störfrequenzen

7.2 Identifikation aus Messdaten

Methoden zur experimentellen Bestimmung:

1. Sprungantwort: Bestimmung von τ und K für PT1-Systeme
2. Frequenzgangmessung: Sinus-Sweep mit Analyse der Amplitude/Phase
3. PRBS-Anregung: Pseudozufällige Binärsignale für breitbandige Identifikation
4. Subspace-Methoden: Für MIMO-Systeme (z. B. in MATLAB mit n4sid)

7.3 Digitale Implementierung

Für diskrete Systeme (z-Bereich):

G(z) = K · (z + a)/(z + b) (Tustin-Transformation: s = 2/z-1 · (1-z⁻¹)/T)

Wichtig: Berücksichtigen Sie die Abtastzeit T bei der Diskretisierung!

8. Softwaretools für die Analyse

8.1 MATLAB/Simulink

• tf(): Erzeugung von Übertragungsfunktionen
• bode(): Bode-Diagramm Plot
• nyquist(): Nyquist-Diagramm
• step(): Sprungantwort
• margin(): Phasen- und Amplitudenreserve

8.2 Python (SciPy/Control)

from scipy import signal
import matplotlib.pyplot as plt

# Übertragungsfunktion definieren
system = signal.TransferFunction([1], [1, 2, 100])

# Bode-Diagramm plotten
w, mag, phase = signal.bode(system)
plt.figure()
plt.subplot(2,1,1)
plt.semilogx(w, mag)
plt.subplot(2,1,2)
plt.semilogx(w, phase)
plt.show()
            

8.3 Octave/Scilab

Open-Source-Alternativen zu MATLAB mit ähnlicher Syntax für Übertragungsfunktionsanalysen.

9. Zusammenfassung und Ausblick

Die Übertragungsfunktion ist ein mächtiges Werkzeug zur Analyse und zum Entwurf von dynamischen Systemen. Dieser Leitfaden hat gezeigt:

• Mathematische Grundlagen von PT1-, PT2- und Systemen höherer Ordnung
• Praktische Berechnungsmethoden und Fallstricke
• Anwendungen in elektrischen, mechanischen und regelungstechnischen Systemen
• Softwaretools für effiziente Analysen

Für fortgeschrittene Anwendungen empfiehlt sich die Vertiefung in:

• Zustandsraumdarstellung (für MIMO-Systeme)
• Robuste Regelung (H∞-Entwurf)
• Nichtlineare Systemidentifikation (Volterra-Reihen, neuronale Netze)

Mit den in diesem Leitfaden vorgestellten Methoden und dem interaktiven Rechner können Sie Übertragungsfunktionen präzise analysieren und für den Entwurf leistungsfähiger Regelsysteme nutzen.