Wenn die Filterfunktionen im zeitkontinuierlichen Raum vorliegen, müssen diese zur Implementierung als digitales Filter in den zeitdiskreten Raum transformiert werden. Der Grund dafür ist, dass wir bei der digitalen Signalverarbeitung zeitlich abgetastete Signale haben. Das erfolgt durch die Ersetzung der Punkte $s$ in der S-Ebene:

(1) $\begin{equation*} \begin{aligned} s &\leftarrow K \frac{1-z^{-1}}{1+z^{-1}} \text{ mit } \\ K&\approx \frac{2}{T} \end{aligned} \end{equation*}$

wobei $T$ das zeitliche Abtastintervall, also der Kehrwert von der Abtastrate ist: $F_s=\frac{1}{T}$ .

Zur Frequenzentzerrung muss wird $K$ folgendermaßen ersetzt:

(2) $\begin{equation*} K\approx \frac{\omega_0}{\tan(\frac{\omega_0\ T}{2}} \end{equation*}$

Somit können die Variablen der zeitkontinuerlichen Achse $s$ ersetzt werden, um die Gleichung in den z-Bereich zu transformieren.

Diese führt jedoch zu sehr komplexen Gleichungen, wonach die analytische Berechnung der Koeffizienten ohne Softwarehilfe sehr schwierig wird. Daher müssen manche Vereinfachungen vorgenommen werden. Dazu wurden die $s$ –Ersatzgleichungen von Robert Bristow-Johnson aus seinem Audio-Kochbuch genommen.

Vereinfachung der Entzerrung nach Robert Bristow-Johnson

Hierfür wurde zur Entzerrung folgender Ersatz genommen.

(3) $\begin{equation*} K\approx \frac{1}{\tan(\frac{\omega_0}{ 2})} \end{equation*}$

mit

(4) $\begin{equation*} \begin{aligned} \tan(\frac{\omega_0}{ 2})&=\frac{\sin(\omega_0)}{1+\cos(\omega_0)} \text{ und } \\ \tan^2(\frac{\omega_0}{ 2})&=\frac{1-\cos(\omega_0)}{1+\cos(\omega_0)} \end{aligned} \end{equation*}$

ergeben sich die Ersatzgleichungen

(5) $\begin{equation*} \begin{aligned} s &\leftarrow \frac{\sin(\omega_0)}{1+\cos(\omega_0)} \cdot \frac{1-z^{-1}}{1+z^{-1}} \\ s^2 &\leftarrow \frac{1-\cos(\omega_0)}{1+\cos(\omega_0)} \cdot \frac{1-2z^{-1}+z^{-2}}{1+2z^{-1}+z^{-2}} \end{aligned} \end{equation*}$

Wenn man nun einen weiteren Vereinfachungsschritt vornimmt und alle Gleichungen mit

(6) $\begin{equation*} $(1+\cos(\omega_0))\cdot (1+2z^{-1}+z^{-2}) \end{equation*}$

multipliziert, erhält man

(7) $\begin{equation*} \begin{aligned} 1 &\leftarrow (1-\cos(\omega_0))\cdot (1+2z^{-1}+z^{-2})\\ s &\leftarrow (1-z^{-2}) \sin(\omega_0) \\ s^2 &\leftarrow (1-2z^{-1}+z^{-2}) (1+\cos(\omega_0)) \text{.} \end{aligned} \end{equation*}$

Nun können die $s$ in der zeitkontinuerlichen Übertragungsfunktion ersetzt und die zeitdiskrete Differentialgleichungen berechnet werden.

1 Kommentar

robert bristow-johnson
1. Januar 2020 — 23:34

Antwort

this looks pretty good. i wouldn’t have other suggestions for you at this time.

keep in mind that the Bilinear Transform is not the only mapping of z to s. there is Euler’s forward method and Euler’s backward method and other methods for discretizing differential equations.

Antwort

Can Kosar

Bilineare Transformation

Vereinfachung der Entzerrung nach Robert Bristow-Johnson

1 Kommentar

Schreibe einen Kommentar Antworten abbrechen