Impulsantwort
# Impulsantwort

Die Impulsantwort, auch Gewichtsfunktion oder Stoßantwort genannt, ist das Ausgangssignal eines Systems, dem am Eingang ein Dirac-Impuls zugeführt wird. Sie wird in der Systemtheorie zur Charakterisierung linearer zeitinvarianter Systeme (LTI-Systeme) benutzt. Der Dirac-Impuls wird gern für theoretische Betrachtungen verwendet, da er ein unendlich weites, kontinuierliches Frequenzspektrum besitzt und das invariante Element der Faltung darstellt.

Bei der experimentellen Analyse dagegen werden LTI-Systeme häufig mit der Sprungfunktion angeregt und die Sprungantwort gemessen, die das Übertragungsverhalten eines solchen Systems ebenfalls vollständig beschreibt. Dadurch vermeidet man es, einen Dirac-Impuls in guter Näherung zu erzeugen, wofür das Eingangssignal kurzzeitig einen sehr hohen Wert annehmen muss.

## Eigenschaften der Impulsantwort

Im Folgenden bezeichnet $g(t)$ die Impulsantwort als Reaktion des Systems auf einen Einheitsimpuls $\delta (t)$, der bei zeitkontinuierlichen Systemen durch den Dirac-Impuls repräsentiert wird.

Bei kausalen Systemen erscheint die Wirkung nicht vor der Ursache und deshalb gilt für deren Impulsantwort $g(t)=0$ für $t<0$. Für das Beispiel eines PT1-Gliedes mit dem Verstärkungsfaktor $K$ und der positiven Zeitkonstante $T$ (im Bild mit $K=1$ und $T=1$ als rote Kurve dargestellt) lautet die Impulsantwort in Form einer Fallunterscheidung

    $g(t)={\begin{cases}0&{\text{wenn }}t<0\\{\frac {K}{T}}\cdot e^{-{\frac {t}{T}}}&{\text{wenn }}t\geq 0\end{cases}}$

„Geschlossener“ wirkt die alternative Schreibweise mit der Sprungfunktion $\sigma (t)$ als Faktor in der Impulsantwort:

    $g(t)=\sigma (t)\cdot {\frac {K}{T}}\cdot e^{-{\frac {t}{T}}}$

Die Impulsantwort ist bei LTI-Systemen die Ableitung der Sprungantwort $h(t)$ nach der Zeit:

    $g(t)={\frac {\mathrm {d} h(t)}{\mathrm {d} t}}$

Generell ist zu beachten, dass der Dirac-Impuls und die Impulsantwort sowie die Differentiation und Integration im erweiterten Sinne der Analysis der Distributionen zu betrachten sind. Besonders zu behandeln ist eine oft bei $t=0$ vorhandene Unstetigkeit der Sprungantwort. Beispielsweise liefert ein einfacher RC-Hochpass die Sprungantwort $h(t)=\sigma (t)\cdot e^{-{\frac {t}{RC}}}$. Zum Differenzieren muss die Produktregel verwendet werden:

    $g(t)={\frac {d\sigma (t)}{dt}}\cdot e^{-{\frac {t}{RC}}}+\sigma (t)\cdot {\frac {d}{dt}}e^{-{\frac {t}{RC}}}=\delta (t)-\sigma (t)\cdot {\frac {e^{-{\frac {t}{RC}}}}{RC}}$

In diesem Fall wird also der Dirac-Impuls in die Impulsantwort „durchgereicht“.

Bei stabilen linearen Systemen ist die Impulsantwort absolut integrierbar:

    $\int \limits _{-\infty }^{\infty }|g(t)|\,dt<K<\infty$

## Die Impulsantwort als Systemcharakteristik

Besitzt ein (zeitkontinuierliches) LTI-System die Impulsantwort $g(t)$, dann kann man seine Reaktion am Ausgang $y(t)$ auf ein beliebiges Eingangssignal $x(t)$ aufgrund des geltenden Überlagerungssatzes durch die (kommutative) Faltung von Impulsantwort und Eingangssignal berechnen:

    $y(t)=g(t)\ast x(t)=x(t)\ast g(t)=\int \limits _{-\infty }^{\infty }g(t-\tau )\cdot x(\tau )d\tau =\int \limits _{-\infty }^{\infty }g(\tau )\cdot x(t-\tau )d\tau$

Wenn (aufgrund der Kausalität) sowohl die Impulsantwort als auch das Eingangssignal für negative Zeiten verschwinden, dann braucht nur von $0$ bis $t$ integriert zu werden.

Unter der Bedingung, dass sich das System zur Zeit $t<0$ im sogenannten 0-Zustand befand (also „energiefrei“ bzw. „entladen“ war), charakterisiert die Impulsantwort ein LTI-System im Zeitbereich vollständig.

Praktisch angewendet wird dieses Prinzip in jüngster Zeit in einigen DirectX- und VST-Plug-ins (siehe Faltungshall), die akustische LTI-Systeme (Räume, Mikrofone, …) virtuell nachbilden.

## Impulsantwort und Übertragungsfunktion

Besondere Bedeutung hat die Reaktion eines stabilen Systems auf die harmonische Exponentielle $e^{j\omega t}$ (mit der „imaginären Kreisfrequenz“ $j\omega$), die durch die Übertragungsfunktion $G(j\omega )$ (in ihrer Form als Frequenzgang) beschrieben wird. Diese ist definiert durch

    $g(t)\ast e^{j\omega t}=\int \limits _{-\infty }^{\infty }g(\tau )\cdot e^{j\omega (t-\tau )}d\tau =e^{j\omega t}\cdot \int \limits _{-\infty }^{\infty }g(\tau )\cdot e^{-j\omega \tau }d\tau =G(j\omega )\cdot e^{j\omega t}$

Die Übertragungsfunktion ist also die Fourier-Transformierte der Impulsantwort und charakterisiert ein LTI-System im Frequenzbereich:

    $G(j\omega )=\int \limits _{-\infty }^{\infty }g(t)\cdot e^{-j\omega t}dt={\mathcal {F}}\{g(t)\}$

Damit lässt sich beispielsweise die Übertragungsfunktion eines PT1-Gliedes aus dessen Impulsantwort berechnen:

    $G(j\omega )=\int \limits _{0}^{\infty }{\frac {K}{T}}\cdot e^{-{\frac {t}{T}}}\cdot e^{-j\omega t}dt={\frac {K}{T}}\cdot \int \limits _{0}^{\infty }e^{-{\frac {1}{T}}\left(1+j\omega T\right)\cdot t}dt={\frac {K}{1+j\omega T}}$

Geht man zur komplexen Frequenz $s=\sigma +j\omega$ über, dann ist die Übertragungsfunktion $G(s)$ die Laplace-Transformierte der Impulsantwort:

    $G(s)={\mathcal {L}}\{g(t)\}$

Für das obige PT1-Glied erhält man damit die (per Definition bekannte) Übertragungsfunktion

    $G(s)={\frac {K}{1+sT}}$

## Ermittlung der Impulsantwort

### Berechnung aus der Differentialgleichung

Ist die Struktur des Systems bekannt, kann daraus seine Differentialgleichung oder sein Differentialgleichungssystem ermittelt werden. Durch den Dirac-Impuls am Systemeingang tritt dieser jedoch als Störfunktion auf der „rechten Seite“ der Differentialgleichung auf. Deshalb versagen die klassischen Lösungsmethoden der linearen Differentialgleichung im Zeitbereich. Aus diesem Grund errechnet man üblicherweise die Sprungantwort und differenziert diese zur Impulsantwort. Wird das System in Zustandsraumdarstellung beschrieben, dann gibt es Lösungsformeln zur Ermittlung der Impulsantwort (siehe Gewichtsmatrix).

### Berechnung aus der Übertragungsfunktion

Ist die Übertragungsfunktion $G(j\omega )$ oder $G(s)$ des Systems durch Messung, Berechnung (beispielsweise durch Transformation der Differentialgleichung in den Frequenzbereich) oder direktes Ablesen aus der Struktur („Symbolische Methode“) schon ermittelt, dann kann daraus durch Fourier- bzw. Laplace-Rücktransformation die Impulsantwort (in manchen Fällen nur auf der Basis von Distributionen) errechnet oder (wie im obigen Beispiel des PT1-Gliedes) direkt aus der Korrespondenztabelle abgelesen werden:

    $g(t)={\mathcal {F}}^{-1}\{G(j\omega )\}={\mathcal {L}}^{-1}\{G(s)\}$

### Messung mittels Dirac-Impuls

Theoretisch kann die Impulsantwort eines Systems durch das Zuführen eines Dirac-Impulses bestimmt werden. Allerdings ist es praktisch unmöglich, einen solchen Impuls zu erzeugen (unendlicher Augenblickswert in verschwindend geringer Zeit), er kann nur in begrenztem Umfang angenähert werden. Dazu müsste ein möglichst kurzer, starker „Knall“ oder Stromstoß auf das System gegeben und seine Antwort über ein Mikrofon o. ä. gemessen werden. Bei auf diese Weise ermitteltem Frequenzgang kann es zu Verzerrungen kommen, vor allem wegen Nichtlinearitäten der Bauteile (Klirrfaktor), Rauschen, Messungenauigkeiten und begrenzter Belastbarkeit.

Die Impulsantwort liefert bei Lautsprecherboxen eine Aussage über die Impulstreue, bei Räumen über das Zeit- und Frequenzverhalten des Nachhalles.

### Ermittlung mittels Sprungantwort

Aus der Sprungantwort eines Systems erhält man durch Differenzieren die Impulsantwort. Aufgrund des plötzlichen Anstiegs der Sprungfunktion gibt es bei deren Messung jedoch ähnliche Probleme wie bei der direkten Messung der Impulsantwort.

### Ermittlung mit einem breitbandigen Signal

Die Impulsantwort kann auch mit einem breitbandigen Rauschsignal, wie weißem Rauschen, bestimmt werden. Dafür sendet man das Rauschsignal in das System (z. B. über einen Lautsprecher in einen Raum) und misst gleichzeitig die Antwort des Systems für eine Weile (zeichnet bspw. mit einem Mikrofon eine Zeitlang auf). Anschließend berechnet man die Kreuzkorrelation des gesendeten und des empfangenen Signals, sie ist in diesem Fall direkt die Impulsantwort des Systems.

Ein großer Vorteil dieser Methode ist, dass neben dem Testsignal noch weitere Signale am System anliegen dürfen. Bspw. muss es in einem Raum zur Messung nicht ruhig sein, solange die Störgeräusche (z. B. Gespräche) unkorreliert zum Testsignal sind, denn sie fallen im Anschluss durch die Kreuzkorrelation heraus.

## Verallgemeinerungen

### Die Gewichtsmatrix

Bei LTI-Systemen mit mehreren Ein- und Ausgangssignalen (sogenannten multivariablen Systemen) mit dem Eingangssignalvektor $\mathbf {x} (t)={\begin{pmatrix}x_{1}\\x_{2}\\\vdots \end{pmatrix}}$ und dem Ausgangssignalvektor $\mathbf {y} (t)={\begin{pmatrix}y_{1}\\y_{2}\\\vdots \end{pmatrix}}$ existiert für jedes Eingangs-/Ausgangs-Paar $(x_{i},y_{j})$ eine „eigene“ Impulsantwort $y_{j}=g_{ji}\ast x_{i}$. Diese Impulsantworten fasst man in der Gewichtsmatrix $\mathbf {G} (t)={\begin{pmatrix}g_{11}&g_{12}&\dots \\g_{21}&g_{22}&\dots \\\dots &\dots &\ddots \end{pmatrix}}$ zusammen, so dass das Ein-/Ausgangsverhalten des multivariablen Systems in der üblichen Matrizenschreibweise dargestellt werden kann:

    $\mathbf {y} (t)=\int \limits _{-\infty }^{\infty }\mathbf {G} (t-\tau )\cdot \mathbf {x} (\tau )\,d\tau =\mathbf {G} (t)\ast \mathbf {x} (t)$

Wird ein zeitkontinuierliches lineares zeitinvariantes System in Zustandsraumdarstellung beschrieben

    ${\frac {d\mathbf {z} }{dt}}=\mathbf {A} \cdot \mathbf {z} +\mathbf {B} \cdot \mathbf {x}$
    $\mathbf {y} =\mathbf {C} \cdot \mathbf {z} +\mathbf {D} \cdot \mathbf {x}$

dann kann die Gewichtsmatrix wie folgt errechnet werden:

    $\mathbf {G} (t)=\sigma (t)\cdot \mathbf {C} \cdot e^{\mathbf {A} \cdot t}\cdot \mathbf {B} +\delta (t)\cdot \mathbf {D}$

Dabei besteht die wesentliche Schwierigkeit „nur noch“ darin, die sogenannte Übergangsmatrix $\mathbf {\Phi } (t)=e^{\mathbf {A} \cdot t}$ aus dem Matrixexponential nach einer „etablierten Methode“ zu ermitteln.

### Die Impulsantwort von zeitvarianten Systemen

Während bei einem zeitinvarianten System die Impulsantwort nicht vom Zeitpunkt des Einheitsimpulses abhängt, ist das bei zeitvarianten Systemen der Fall. Dann definiert man die verallgemeinerte Impulsantwort $g(t,\tau )$ als Reaktion auf den Einheitsimpuls $\delta (t-\tau )$ zum Zeitpunkt $\tau$. Auch diese Form kann bei linearen zeitvarianten Systemen als Systemcharakteristik dienen. Dabei wird das Ausgangssignal aus dem Eingangssignal durch folgendes Integral ermittelt:

    $y(t)=\int \limits _{-\infty }^{\infty }x(\tau )\cdot g(t,\tau )\,d\tau$

Die Berechnung der Impulsantwort für zeitvariable Systeme ist im Allgemeinen wesentlich schwieriger. Eine aus der Impulsantwort ermittelte Übertragungsfunktion wird in diesem Fall zeitabhängig. Für multivariable Systeme wird analog eine verallgemeinerte Gewichtsmatrix $\mathbf {G} (t,\tau )$ definiert.

## Die Impulsantwort von zeitdiskreten Systemen

Bei zeitdiskreten Systemen, z. B. digitalen Filtern, werden die Signale durch Folgen (mit den natürlichen Zahlen als Indexmenge) repräsentiert. Der diskrete Einheitsimpuls wird durch eine Folge dargestellt, bei der nur das 0-te Element gleich 1 und alle weiteren Elemente gleich 0 sind:

    $\delta (t)=(1,0,0,0,\dots )$

Er ist das neutrale Element der diskreten Faltung von Folgen.

Die zeitdiskrete Impulsantwort ist dann die Systemreaktion auf den diskreten Einheitsimpuls in Form einer Folge

    $g(t)=\left(g(0),g(1),g(2),g(3),\dots \right)$

Beispielsweise lautet die Impulsantwort für den einfachen gleitenden Mittelwert dritter Ordnung

    $g(t)=\left({\frac {1}{3}},{\frac {1}{3}},{\frac {1}{3}},0,0,0,\dots \right)=\left({\frac {1}{3}},{\frac {1}{3}},{\frac {1}{3}},{\overline {0}}\right)$

Im Gegensatz zur Impulsantwort bei zeitkontinuierlichen Systemen gibt es hierbei aufgrund der Endlichkeit des Einheitsimpulses keine Probleme bei der rechnerischen oder experimentellen Ermittlung.

Die zeitdiskrete Impulsantwort charakterisiert ein zeitdiskretes LTI-System im Zeitbereich, denn dessen Ausgangsfolge $y(t)$ berechnet sich (bei kausalen Systemen und vom 0-Zustand ausgehend) durch diskrete Faltung von Impulsantwort $g(t)$ und Eingangsfolge $x(t)$

    $y(t)=g(t)\ast x(t)=\sum \limits _{\tau =0}^{t}g(t-\tau )\cdot x(\tau )=\sum \limits _{\tau =0}^{t}g(\tau )\cdot x(t-\tau )$

Dabei spielt es keine Rolle, ob die Komponenten der Folgen einer endlichen (z. B. bei linearen Automaten und Faltungscodierern) oder einer unendlichen Symbolmenge (z. B. bei Abtastsystemen, FIR- und IIR-Filtern) entstammen. Wesentlich ist, dass die Folgen einen linearen Raum bilden. Praktisch werden derartige Folgen durch formale Potenzreihen beschrieben, die man durch z-Transformation, Zeta-Transformation oder eine andere diskrete Operatorenrechnung erhält. Die (zeitdiskret) transformierte Impulsantwort kann als Übertragungsfunktion interpretiert werden. Für das obige Beispiel wird diese (auf Basis der z-Transformation)

    $G(z)={\mathcal {Z}}\{g(t)\}={\mathcal {Z}}\left\{\left({\frac {1}{3}},{\frac {1}{3}},{\frac {1}{3}},{\overline {0}}\right)\right\}={\frac {z^{0}}{3}}+{\frac {z^{-1}}{3}}+{\frac {z^{-2}}{3}}={\frac {z^{2}+z+1}{3\cdot z^{2}}}$

Auch bei zeitdiskreten linearen Systemen lässt sich die Impulsantwort für zeitvariante und/oder multivariable Systeme verallgemeinern.

