Bilineare Transformation (Signalverarbeitung)
# Bilineare Transformation (Signalverarbeitung)

Die bilineare Transformation, im englischen Sprachraum auch als Tustin's method (dt. „Tustin-Methode“) bezeichnet, ist in der Signalverarbeitung eine Transformation zwischen der zeitkontinuierlichen und der zeitdiskreten Darstellung von Systemfunktionen. Sie spielt in der digitalen Signalverarbeitung und der Regelungstheorie eine Rolle, da sie einen Bezug in der Systembeschreibung zwischen analogen, kontinuierlichen Systemen und digitalen, diskreten Systemen herstellt.

## Motivation

In der Signalverarbeitung und Regelungstechnik besteht mittels bilinearer Transformation die Möglichkeit, zeitkontinuierliche Übertragungsfunktionen $G(s)$ von linearen, zeitinvarianten Systemen in zeitdiskrete Übertragungsfunktionen $H\left[z\right]$ mit ähnlichem Verhalten umzuwandeln. Die Transformation kann in beide Richtungen erfolgen. Die Übertragungsfunktion $G(s)$ kann beispielsweise ein analoges Filter beschreiben und $H\left[z\right]$ stellt eine aus dem analogen Filter abgeleitete, zeitdiskrete Übertragungsfunktion dar, welche ein äquivalentes digitales Filter beschreibt.

Die Beschreibung der Systemfunktionen von zeitkontinuierlichen Systemen erfolgt in der so genannten $s$-Ebene und ihre Analyse erfolgt mittels der Laplace-Transformation. Bei zeitdiskreten Systemen erfolgt die Darstellung in der so genannten $z$-Ebene und die Analyse erfolgt mittels der Z-Transformation. Eine mögliche Transformation von Systemen zwischen der $s$- und $z$-Ebene besteht in Form der bilinearen Transformation. Die bilineare Transformation bietet gegenüber anderen Verfahren wie der Impulsinvarianzmethode und der Matched-Z-Transformation den Vorteil Alias-Effekte im zeitdiskreten System zu vermeiden. Der damit verknüpfte Nachteil ist eine nichtlineare Verzerrung bei dem Übergang der Übertragungsfunktionen von $G(s)$ zu $H\left[z\right]$.

## Beschreibung

Die bilineare Transformation ist eine konforme Abbildung und eine Anwendung der Möbiustransformation. Sie ordnet jedem Punkt $s=\sigma +j\cdot \Omega$ in der komplexen $S$-Ebene eindeutig einen bestimmten Punkt in der komplexen $z$-Ebene zu und umgekehrt, wie in nebenstehender Abbildung für verschiedene Werte von $\sigma$ und $\Omega$ grafisch dargestellt. Beispielsweise werden die Werte auf der imaginären Achse $\Omega$, in rot dargestellt, auf den Einheitskreis $\vert z\vert =1$ in der $z$-Ebene abgebildet. Alle Punkte in der linken $s$-Ebene mit negativem Realwert werden in der $z$-Ebene auf Punkte innerhalb des rot eingezeichneten Einheitskreises abgebildet – dieser Umstand ist für Stabilitätsuntersuchungen linearer Systeme wesentlich, da stabile Systeme mit Polstellen in der linken $s$-Ebene in zeitdiskrete Systeme mit Polstellen innerhalb des Einheitskreises übergehen.

Die zeitkontinuierliche Systemfunktion $G(s)$ korrespondiert bei der bilinearen Transformation mit der zeitdiskreten Systemfunktion $H\left[z\right]$ durch die Substitution der Variablen $s$ in der Form:

    $s={\frac {2}{T}}\cdot {\frac {z-1}{z+1}}={\frac {2}{T}}\cdot {\frac {1-z^{-1}}{1+z^{-1}}}$

was bedeutet:

    $H(z)=G\left[{\frac {2}{T}}\cdot {\frac {z-1}{z+1}}\right]$

Der Parameter $T$ stellt das zeitliche Abtastintervall (Periodendauer) dar. Der Kehrwert $f_{s}={\tfrac {1}{T}}$ wird als Abtastrate bezeichnet. Die umgekehrte Zuordnung ergibt sich mit $s=\sigma +j\cdot \Omega$ zu:

    $z={\frac {1+Ts/2}{1-Ts/2}}={\frac {1+\sigma T/2+j\Omega T/2}{1-\sigma T/2-j\Omega T/2}}$

Wird der Realteil von $s$ gleich 0 gesetzt ($s=j\cdot \Omega$) ergibt sich:

    $z={\frac {1+j\Omega T/2}{1-j\Omega T/2}}$

Der Betrag von $z$ ist dann für alle Werte von $\Omega$ gleich 1 ($\vert z\vert =1$), was der Abbildung der imaginären Achse der $s$-Ebene auf den Einheitskreis in der $z$-Ebene entspricht.

### Frequenzverzerrung

Durch den Umstand, dass der kontinuierliche Frequenzbereich $-\infty \leq \Omega \leq +\infty$ der $s$-Ebene auf den $-\pi \leq \omega \leq \pi$ am Einheitskreis der $z$-Ebene abgebildet wird, muss die Transformation von der zeitkontinuierlichen zur zeitdiskreten Frequenzvariablen nichtlinear sein. Um die Beziehung zwischen der Frequenzachse $\Omega$ in der $s$-Ebene und dem Einheitskreis mit Winkel $\omega$ in der $z$-Ebene abzuleiten, wird $z$ mit $e^{j\omega T}$ substituiert:

    $z=e^{j\omega T}={\frac {1+j\Omega T/2}{1-j\Omega T/2}}$.

Dies entspricht ja genau dem Bild der Frequenzachse der $s$-Ebene. Daraus lässt sich nun mit Hilfe der bilinearen Transformation $s$ bestimmen zu:

    $s={\frac {2}{T}}{\frac {e^{j\omega T}-1}{e^{j\omega T}+1}}=\sigma +j\Omega ={\frac {2}{T}}\left[{\frac {2\cdot e^{-j{\frac {\omega T}{2}}}\cdot j\sin({\frac {\omega T}{2}})}{2\cdot e^{-j{\frac {\omega T}{2}}}\cdot \cos({\frac {\omega T}{2}})}}\right]={\frac {2j}{T}}\tan \left({\frac {\omega T}{2}}\right)$.

Mit $\sigma =0$ führt dies zu der Beziehung:

    $\Omega ={\frac {2}{T}}\tan \left({\frac {\omega T}{2}}\right)$

beziehungsweise auf den rechts in der Abbildung dargestellten Verlauf $\omega (\Omega )$:

    $\omega =2\arctan \left({\frac {\Omega T}{2}}\right)$.

Die bilineare Transformation vermeidet Alias-Effekte durch „Kompression“ der gesamten imaginären Achse $\Omega$ auf den Einheitskreis in der $z$-Ebene. Die resultierende nichtlineare Kompression der Frequenzachse stellt eine Frequenzverzerrung dar und muss beispielsweise im Rahmen des Filterentwurfes beachtet werden, wenn analoge (zeitkontinuierliche) Filter wie elliptische Filter als zeitdiskrete, digitale IIR-Filter realisiert werden sollen. In diesen Fällen ist eine Vorverzerrung der kontinuierlichen Übertragungsfunktion $G(s)$ des Filters notwendig, neben Beachtung der Nyquistbandbreite, um nach der bilinearen Transformation die passende zeitdiskrete Übertragungsfunktion $H\left[z\right]$ zu erhalten.

