Lineare zeitinvariante Systeme/Systembeschreibung im Zeitbereich: Unterschied zwischen den Versionen

Aus LNTwww
Wechseln zu:Navigation, Suche
Zeile 38: Zeile 38:
  
 
==Einige Gesetze der Fourier–Transformation==
 
==Einige Gesetze der Fourier–Transformation==
Die [[Signaldarstellung/Gesetzm%C3%A4%C3%9Figkeiten_der_Fouriertransformation|Gesetzmäßigkeiten der Fouriertransformation]] des Buches [[Signaldarstellung]] wurden die Gesetzmäßigkeiten der Fouriertransformation ausführlich dargelegt. Hier folgt nun eine kurze Zusammenfassung, wobei $H(f)$ den Frequenzgang eines LZI–Systems beschreibt und dessen Fourierrücktransformierte $h(t)$ die Impulsantwort angibt.  
+
Die [[Signaldarstellung/Gesetzm%C3%A4%C3%9Figkeiten_der_Fouriertransformation|Gesetzmäßigkeiten der Fouriertransformation]] wurden bereits im Buch [[Signaldarstellung]] ausführlich dargelegt. Hier folgt nun eine kurze Zusammenfassung, wobei $H(f)$ den Frequenzgang eines LZI–Systems beschreibt und dessen Fourierrücktransformierte $h(t)$ die Impulsantwort angibt.  
  
 
Bei allen folgenden Gleichungen wird das Kurzsymbol der Fouriertransformation benutzt. Der ausgefüllte Kreis kennzeichnet stets den Spektralbereich und der weiße Kreis den Zeitbereich.
 
Bei allen folgenden Gleichungen wird das Kurzsymbol der Fouriertransformation benutzt. Der ausgefüllte Kreis kennzeichnet stets den Spektralbereich und der weiße Kreis den Zeitbereich.
*Multiplikation mit einem konstanten Faktor:
+
*'''Multiplikation mit einem konstanten Faktor''':
 
$$k \cdot H(f)\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,k \cdot h(t).$$
 
$$k \cdot H(f)\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,k \cdot h(t).$$
 
:Bei $k$ < 1 spricht man von einer Dämpfung, während $k$ > 1 für eine Verstärkung steht.
 
:Bei $k$ < 1 spricht man von einer Dämpfung, während $k$ > 1 für eine Verstärkung steht.
*Ähnlichkeitssatz:
+
*'''Ähnlichkeitssatz''':
 
$$H({f}/{k})\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,|k| \cdot h(k\cdot t).$$
 
$$H({f}/{k})\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,|k| \cdot h(k\cdot t).$$
 
:Dies besagt: Eine Stauchung ($k$ < 1) des Frequenzgangs führt zu einer breiteren und niedrigeren Impulsantwort. Durch Streckung ($k$ > 1) von $H(f)$ wird $h(t)$ schmaler und höher.  
 
:Dies besagt: Eine Stauchung ($k$ < 1) des Frequenzgangs führt zu einer breiteren und niedrigeren Impulsantwort. Durch Streckung ($k$ > 1) von $H(f)$ wird $h(t)$ schmaler und höher.  
*Verschiebungssatz im Frequenz- und Zeitbereich:
+
*'''Verschiebungssatz im Frequenz- und Zeitbereich''':
 
$$H(f - f_0) \bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, h( t )\cdot {\rm e}^{\hspace{0.05cm}{\rm j}2\pi f_0 t},$$
 
$$H(f - f_0) \bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, h( t )\cdot {\rm e}^{\hspace{0.05cm}{\rm j}2\pi f_0 t},$$
 
$$H(f) \cdot {\rm e}^{-{\rm j}2\pi ft_0}\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, h( t- t_0 ).$$
 
$$H(f) \cdot {\rm e}^{-{\rm j}2\pi ft_0}\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, h( t- t_0 ).$$
 
:Eine Verschiebung um $t_0$ (Laufzeit) führt also im Frequenzbereich zu der Multiplikation mit einer komplexen Exponentialfunktion. Der Amplitudengang $|H(f)|$ wird dadurch nicht verändert.
 
:Eine Verschiebung um $t_0$ (Laufzeit) führt also im Frequenzbereich zu der Multiplikation mit einer komplexen Exponentialfunktion. Der Amplitudengang $|H(f)|$ wird dadurch nicht verändert.
*Differentiationssatz im Frequenz- und Zeitbereich:
+
*'''Differentiationssatz im Frequenz- und Zeitbereich''':
 
$$\frac{1}{{{\rm j}2\pi }} \cdot \frac{{{\rm d}H( f )}}{{{\rm d}f}}\cdot \bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,- t \cdot h( t ),$$
 
$$\frac{1}{{{\rm j}2\pi }} \cdot \frac{{{\rm d}H( f )}}{{{\rm d}f}}\cdot \bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,- t \cdot h( t ),$$
 
$${\rm j}\cdot 2\pi f \cdot H( f ){}\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, \frac{{{\rm d}h( t )}}{{\rm d}t}.$$
 
$${\rm j}\cdot 2\pi f \cdot H( f ){}\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, \frac{{{\rm d}h( t )}}{{\rm d}t}.$$

Version vom 28. November 2016, 11:26 Uhr

Impulsantwort

Auf der Seite Das erste Fourierintegral im Buch „Signaldarstellung” wurde dargelegt, dass für jedes deterministische Signal $x(t)$ mit Hilfe der Fouriertransformation eine Spektralfunktion $X(f)$ angegeben werden kann. Oft bezeichnet man $X(f)$ auch kurz als das Spektrum.

Alle Informationen über die Spektralfunktion sind aber auch bereits in der Zeitbereichsdarstellung enthalten, wenn auch nicht immer sofort erkennbar. Der gleiche Sachverhalt trifft für lineare zeitinvariante Systeme zu.

Die wichtigste Beschreibungsgröße eines linearen zeitinvarianten Systems im Zeitbereich ist die Fourierrücktransformierte von $H(f)$, die man als die Impulsantwort bezeichnet: $$h(t) = \int_{-\infty}^{+\infty}H(f) \cdot {\rm e}^{\hspace{0.05cm}{\rm j}2\pi ft}\hspace{0.15cm} {\rm d}f.$$


Hierzu ist Folgendes anzumerken:

  • Der Frequenzgang $H(f)$ und die Impulsantwort $h(t)$ sind äquivalente Beschreibungsgrößen, die genau die gleichen Informationen über das LZI–System beinhalten.
  • Verwendet man das diracförmige Eingangssignal $x(t) = δ(t)$, so ist $X(f) = 1$ zu setzen und es gilt $Y(f) = H(f)$ bzw. $y(t) = h(t)$. Die Bezeichnung „Impulsantwort” spiegelt diese Aussage wieder: Die Antwort des Systems auf einen (Dirac-) Impuls am Eingang.
  • Die obige Definition lässt erkennen, dass jede Impulsantwort die Einheit Hz = 1/s besitzen muss.


Die Impulsantwort $h(t)$ des so genannten Spalt–Tiefpasses ist über eine Zeitdauer $T$ hinweg konstant und außerhalb dieses Zeitintervalls gleich 0. Der dazugehörige Amplitudengang als der Betrag des Frequenzgangs ist $|H(f)| = |{\rm si}(\pi fT)|$. Der Phasenverlauf ergibt sich damit zu $$b(f) = \left\{ \begin{array}{l} \hspace{0.25cm}\pi/T \\ -\pi/T \\ \end{array} \right.\quad \quad\begin{array}{*{20}c} {\rm{f\ddot{u}r}} \\ {\rm{f\ddot{u}r}} \\ \end{array}\begin{array}{*{20}c}{\left| \hspace{0.05cm} f\hspace{0.05cm} \right| > 0,} \\{\left|\hspace{0.05cm} f \hspace{0.05cm} \right|< 0.} \\\end{array}$$ Wäre $h(t)$ symmetrisch um $t = 0$ und damit akausal, so würde $b(f)=0$ gelten.

Rechteckförmige Impulsantwort und zugehöriges Betragsspektrum

Die Fläche über die Impulsantwort ist gleich $H(f = 0) = 1$. Daraus folgt, dass die Impulsantwort im Bereich von 0 bis $T$ gleich $1/T$ sein muss.

Bild ersetzen:

Andere Unterschrift

Einige Gesetze der Fourier–Transformation

Die Gesetzmäßigkeiten der Fouriertransformation wurden bereits im Buch Signaldarstellung ausführlich dargelegt. Hier folgt nun eine kurze Zusammenfassung, wobei $H(f)$ den Frequenzgang eines LZI–Systems beschreibt und dessen Fourierrücktransformierte $h(t)$ die Impulsantwort angibt.

Bei allen folgenden Gleichungen wird das Kurzsymbol der Fouriertransformation benutzt. Der ausgefüllte Kreis kennzeichnet stets den Spektralbereich und der weiße Kreis den Zeitbereich.

  • Multiplikation mit einem konstanten Faktor:

$$k \cdot H(f)\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,k \cdot h(t).$$

Bei $k$ < 1 spricht man von einer Dämpfung, während $k$ > 1 für eine Verstärkung steht.
  • Ähnlichkeitssatz:

$$H({f}/{k})\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,|k| \cdot h(k\cdot t).$$

Dies besagt: Eine Stauchung ($k$ < 1) des Frequenzgangs führt zu einer breiteren und niedrigeren Impulsantwort. Durch Streckung ($k$ > 1) von $H(f)$ wird $h(t)$ schmaler und höher.
  • Verschiebungssatz im Frequenz- und Zeitbereich:

$$H(f - f_0) \bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, h( t )\cdot {\rm e}^{\hspace{0.05cm}{\rm j}2\pi f_0 t},$$ $$H(f) \cdot {\rm e}^{-{\rm j}2\pi ft_0}\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, h( t- t_0 ).$$

Eine Verschiebung um $t_0$ (Laufzeit) führt also im Frequenzbereich zu der Multiplikation mit einer komplexen Exponentialfunktion. Der Amplitudengang $|H(f)|$ wird dadurch nicht verändert.
  • Differentiationssatz im Frequenz- und Zeitbereich:

$$\frac{1}{{{\rm j}2\pi }} \cdot \frac{{{\rm d}H( f )}}{{{\rm d}f}}\cdot \bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\,- t \cdot h( t ),$$ $${\rm j}\cdot 2\pi f \cdot H( f ){}\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, \frac{{{\rm d}h( t )}}{{\rm d}t}.$$

Ein differenzierendes Element im LZI–System führt im Frequenzbereich zu einer Multiplikation mit $j2πf$ und damit unter Anderem zu einer Phasendrehung um 90°.

Diese Gesetzmäßigkeiten werden in den Aufgaben zu den Kapiteln 1.2 und 1.3 häufiger angewendet.

Kausale Systeme

Die Grafik verdeutlicht den Unterschied zwischen einem akausalen und einem kausalen System.

Akausales und kausales System

Ein LZI–System bezeichnet man dann als kausal, wenn die Impulsantwort $h(t)$ – also die Fourierrücktransformierte des Frequenzgangs $H(f)$ – folgende Bedingung erfüllt: $$h(t) = 0 \hspace{0.15cm}{\rm f\ddot{u}r}\hspace{0.15cm} t < 0.$$


Das obige Bild lässt erkennen, dass nach dieser Definition das System B kausal ist im Gegensatz zum System A. Weiter ist anzumerken:

  • Jedes realisierbare System ist kausal. Beim System A beginnt die Wirkung früher (bei $t$ = –$T$) als die Ursache (Diracfunktion bei $t$ = 0), was natürlich in der Praxis nicht möglich ist.
  • Fast alle akausalen Systeme lassen sich unter Verwendung einer Laufzeit $\tau$ in ein realisierbares kausales System überführen. Zum Beispiel gilt mit $\tau$ = $T$:

$$h_{\rm B}(t) = h_{\rm A}(t - T).$$

  • Für kausale Systeme gelten alle in diesem Kapitel getroffenen Aussagen ebenso wie für akausale Systeme. Zur Beschreibung kausaler Systeme lassen sich jedoch einige spezifische Eigenschaften nutzen, wie im Kapitel 3 dieses Buches ausgeführt wird.
  • In diesem und dem nächsten Kapitel betrachten wir vorwiegend akausale Systeme, da deren mathematische Beschreibung anschaulicher ist. So ist der Frequenzgang $H_{\rm A}(f)$ reell, während für $H_{\rm B}(f)$ der zusätzliche Term ${\rm exp(–j2}πfT)$ zu berücksichtigen ist.

Berechnung des Ausgangssignals

Wir betrachten die folgende Aufgabenstellung: Bekannt sei das Eingangssignal $x(t)$ und der Frequenzgang $H(f)$. Gesucht ist das Ausgangssignal $y(t)$.

Zur Ermittlung der Ausgangsgrößen eines LZI–Systems

Soll die Lösung im Frequenzbereich erfolgen, so muss zunächst aus dem gegebenen Eingangssignal $x(t)$ durch Fouriertransformation das Spektrum $X(f)$ ermittelt und mit dem Frequenzgang $H(f)$ multipliziert werden. Durch Fourierrücktransformation des Produktes kommt man zum gesuchten Signal $y(t)$: $${\rm 1.\,\, Schritt:}\hspace{0.5cm} X(f)\bullet\!\!-\!\!\!-\!\!\!-\!\!\circ\, x( t )\hspace{2.05cm}{\rm Eingangsspektrum},$$ $${\rm 2.\,\, Schritt:}\hspace{0.5cm}Y(f)= X(f) \cdot H(f) \hspace{1.32cm}{\rm Ausgangsspektrum},$$ $${\rm 3.\,\, Schritt:}\hspace{0.5cm} y(t)\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, Y(f )\hspace{2.05cm}{\rm Ausgangssignal}.$$


Zum gleichen Ergebnis kommt man nach der Berechnung im Zeitbereich, indem man zunächst aus dem Frequenzgang $H(f)$ mittels Fourierrücktransformation die Impulsantwort $h(t)$ berechnet und anschließend die Faltungsoperation anwendet: $$y(t) = x (t) * h (t) = \int_{ - \infty }^{ + \infty } {x ( \tau )} \cdot h ( {t - \tau } ) \hspace{0.1cm}{\rm d}\tau.$$ Die Ergebnisse sind bei beiden Vorgehensweisen identisch. Zweckmäßigerweise sollte man dasjenige Verfahren auswählen, das mit weniger Rechenaufwand zum Ziel führt.

Am Eingang eines Spalt–Tiefpasses (siehe Beispiel auf der ersten Seite dieses Kapitels) mit rechteckförmiger Impulsantwort der Breite $T$ liegt ein Rechteckimpuls $x(t)$ der Dauer $2T$ an.

Trapezförmiger Ausgangsimpuls, da x(t) und h(t) rechteckförmig sind

In diesem Fall ist die direkte Berechnung im Zeitbereich günstiger: Die Faltung zweier unterschiedlich breiter Rechtecke $x(t)$ und $h(t)$ führt zum trapezförmigen Ausgangsimpuls $y(t)$:

  • Man erkennt die Tiefpasseigenschaft des Filters an der endlichen Flankensteilheit von $y(t)$.
  • Die Impulshöhe 3 V bleibt in diesem Beispiel erhalten, weil $H(f = 0) = 1/T · T = 1$ ist.

Sprungantwort

Eine in der Praxis oft verwendete Eingangsfunktion $x(t)$ zur Messung von $H(f)$ ist die Sprungfunktion $${\rm \gamma}(t) = \left\{ \begin{array}{l} \hspace{0.25cm}0 \\ 0.5 \\ \hspace{0.25cm} 1 \\ \end{array} \right.\quad \quad\begin{array}{*{20}c} {\rm{f\ddot{u}r}} \\ {\rm{f\ddot{u}r}}\\ {\rm{f\ddot{u}r}} \\ \end{array}\begin{array}{*{20}c}{\left| \hspace{0.05cm} t\hspace{0.05cm} \right| < 0,} \\ {\left| \hspace{0.05cm}t\hspace{0.05cm} \right| = 0,} \\ {\left|\hspace{0.05cm} t \hspace{0.05cm} \right| > 0.} \\ \end{array}$$

Die Sprungantwort $\sigma(t)$ ist die Antwort des Systems, wenn man an den Eingang die Sprungfunktion $\gamma(t)$ anlegt: $$x(t) = {\rm \gamma}(t)\hspace{0.5cm}\Rightarrow \hspace{0.5cm} y(t) = {\rm \sigma}(t).$$


Die Berechnung im Frequenzbereich ist etwas umständlich; man muss folgende Gleichung anwenden: $${\rm \sigma}(t)\circ\!\!-\!\!\!-\!\!\!-\!\!\bullet\, X(f ) \cdot H(f) =\left(\frac{1}{2}\cdot \delta(f) + \frac{1}{{\rm j}\cdot 2\pi f} \right) \cdot H(f).$$

Die Berechnung im Zeitbereich führt dagegen direkt zum Ergebnis: $${\rm \sigma}(t) = \int_{ - \infty }^{ t } {h ( \tau )} \hspace{0.1cm}{\rm d}\tau.$$

Bei kausalen Systemen gilt $h(\tau)$ = 0 für $\tau$ < 0, so dass die untere Integrationsgrenze in obiger Gleichung dann zu 0 gesetzt werden kann.

Das genannte Ergebnis ist auch aus folgendem Grunde einsichtig: Die Sprungfunktion $\gamma(t)$ hängt mit der Diracfunktion $\delta(t)$ wie folgt zusammen: $${\rm \gamma}(t) = \int_{ - \infty }^{ t } {\delta ( \tau )} \hspace{0.1cm}{\rm d}\tau.$$ Da wir Linearität vorausgesetzt haben und die Integration eine lineare Operation darstellt, gilt auch für das Ausgangssignal $\sigma(t)$ der entsprechende Zusammenhang entsprechend der vorletzten Gleichung.

Die Grafik verdeutlicht den Sachverhalt für eine rechteckförmige Impulsantwort $h(\tau)$. Die Abszisse wurde in $\tau$ umbenannt. Blau eingezeichnet ist die Sprungfunktion $\gamma(\tau)$.

Zur Berechnung der Sprungantwort bei rechteckförmiger Impulsantwort

Durch Spiegelung und Verschiebung erhält man die rot gestrichelt eingezeichnete Funktion $\gamma(t – \tau)$. Die rot hinterlegte Fläche gibt somit die Sprungantwort $\sigma(\tau)$ zum Zeitpunkt $\tau$ = t an.


Aufgaben

1.3 Gemessene Sprungantwort

Zusatzaufgaben:1.3 Exponentiell abfallendes h(t)

1.4 Zum Tiefpass 2. Ordnung

Zusatzaufgaben:1.4 Alles rechteckförmig