Diskrete Sinustransformation

Die diskrete Sinustransformation (DST, englisch discrete sine transform) ist eine reellwertige, diskrete, lineare, orthogonale Transformation, die ähnlich wie der imaginäre Teil der diskreten Fouriertransformation (DFT) ein zeitdiskretes Signal vom Zeitbereich (bei Zeitsignalen) bzw. dem Ortsbereich (bei räumlichen Signalen) in den Frequenzbereich transformiert.

Sie ist eng verwandt mit der diskreten Kosinustransformation (DCT), basiert aber im Gegensatz auf der ungeraden Sinusfunktion.^[1]

Anwendung der DST, wie auch der DCT, liegen bei der Lösung von partiellen Differentialgleichungen. Bei dem Videostandard H.265 kann die DST bei bestimmten Einstellungen zum Einsatz kommen.^[2] Im Gegensatz zur DCT besitzt die DST in den meisten Fällen keine wesentliche Anwendung im Bereich der Signalverarbeitung und Datenkompression.

Definition

Unterschiedliche periodische Fortsetzungen bei DST-I bis DST-IV an einer 9 Elemente langen Beispielfolge in rot.

Es gibt in Summe acht verschiedene Formen der DST, die in der Literatur mit DST-I bis DST-VIII bezeichnet werden. Sie unterscheiden sich durch die Art, wie die endliche Folge am Anfang der Folge ungerade fortgesetzt wird. Die DST-I bis DST-IV ist, bis auf einen konstanten Faktor, gleichwertig zur reellwertigen, ungeraden DFT mit gerader Ordnung. Die verschiedenen Arten der DST bilden dabei jeweils die reellwertige Eingabefolge, aus dem Orts- bzw. Zeitbereich, mit N Elementen x[n] auf eine reellwertige Ausgabefolge, den Spektralbereich, X[n] ab:

${\displaystyle x[n]=x_{0},\ldots ,x_{N-1}\Rightarrow X[n]=X_{0},\ldots ,X_{N-1}}$

Die vier gebräuchlichsten DST-Arten sind DST-I bis DST-IV:

DST-I

Die DST-I ist bezüglich ihrer Randwerte ungerade am Anfang um x₋₁ und ungerade am Ende um x_N.

${\displaystyle X_{k}=\sum _{n=0}^{N-1}x_{n}\sin \left[{\frac {\pi }{N+1}}(n+1)(k+1)\right]\quad \quad k=0,\dots ,N-1}$

DST-II

Die DST-II ist bezüglich ihrer Randwerte ungerade am Anfang um x_−1/2 und ungerade am Ende um x_N−1/2.

${\displaystyle X_{k}=\sum _{n=0}^{N-1}x_{n}\sin \left[{\frac {\pi }{N}}\left(n+{\frac {1}{2}}\right)(k+1)\right]\quad \quad k=0,\dots ,N-1}$

DST-III

Die DST-III ist bezüglich ihrer Randwerte ungerade am Anfang um x₋₁ und gerade am Ende um x_N−1.

${\displaystyle X_{k}={\frac {(-1)^{k}}{2}}x_{N-1}+\sum _{n=0}^{N-2}x_{n}\sin \left[{\frac {\pi }{N}}(n+1)\left(k+{\frac {1}{2}}\right)\right]\quad \quad k=0,\dots ,N-1}$

DST-IV

Die DST-IV ist bezüglich ihrer Randwerte ungerade am Anfang um x_−1/2 und gerade am Ende um x_N−1/2.

${\displaystyle X_{k}=\sum _{n=0}^{N-1}x_{n}\sin \left[{\frac {\pi }{N}}\left(n+{\frac {1}{2}}\right)\left(k+{\frac {1}{2}}\right)\right]\quad \quad k=0,\dots ,N-1}$

Inverse Transformation

Wie jede Transformation besitzt auch die DST eine inverse Transformation. Die Inverse der DST-I ist die DST-I mit einem konstanten Faktor ²/_(N+1). Die Inverse der DST-IV ist die DST-IV mit dem konstanten Faktor ²/_N. Die Inverse der DST-II ist die DST-III mit einem Faktor ²/_N und umgekehrt.

Ähnlich wie bei der DCT sind die Vorfaktoren der DST in der Literatur nicht einheitlich festgelegt. Beispielsweise wird von manchen Autoren ein zusätzlicher Faktor von ${\displaystyle {\sqrt {2/N}}}$ eingeführt, um den zusätzlichen Faktor bei der inversen Operation zu vermeiden. Durch geeignete Wahl des konstanten Faktors kann die Transformationsmatrix eine orthogonale Matrix darstellen.

Literatur

• Vladimir Britanak, Patrick C. Yip, K. R. Rao: Discrete Cosine and Sine Transforms: General Properties, Fast Algorithms and Integer Approximations. Academic Press, 2007, ISBN 978-0-12-373624-6. 

Weblinks

• The Design and Implementation of FFTW3 (engl.; PDF; 342 kB)
• FFTW. Eine quelloffene C-Bibliothek unter der GPL zur Berechnung der DST-I bis DST-IV ein einer oder mehrerer Dimensionen.

Einzelnachweise

1. ↑ S. A. Martucci: Symmetric convolution and the discrete sine and cosine transforms, in Proceedings of the IEEE in Signal Processing, Ausgabe SP-42, 1994, S. 1038–1051.
2. ↑ Martin Fiedler: Videokompressionsverfahren - von MPEG-1 bis H.264 und H.265. Abgerufen am 10. März 2014.