Abtastung (Signalverarbeitung)

Unter Abtastung (englisch sampling) wird in der Signalverarbeitung die Registrierung von Messwerten zu diskreten, meist äquidistanten Zeitpunkten verstanden. Aus einem zeitkontinuierlichen Signal wird so ein zeitdiskretes Signal gewonnen.

Bei mehrkanaligen Signalen ergibt jede Abtastung ein „Sample“ aus mehreren Abtastwerten. Die Zahl der Samples pro Sekunde wird Abtastrate genannt. Bei der digitalen Telefonie (ISDN) beträgt die Abtastrate beispielsweise 8 kHz.

Abgrenzung

Die Digitalisierung eines analogen Signals im Zeitbereich umfasst als Überbegriff zusätzlich zur Abtastung eine weitere Umwandlung, die Quantisierung, wobei die beiden Umwandlungen in beliebiger Reihenfolge ausgeführt werden können um ein Digitalsignal zu erhalten:

Die Abtastung um ein zeitkontinuierliches Signal in ein zeitdiskretes Signal zu überführen.
Die Quantisierung um ein wertkontinuierliches Signal in ein wertdiskretes Signal umzuwandeln.

Abtastung im Zeitbereich

Ideale Abtastung

Für eine einfachere mathematische Beschreibung ist die ideale Abtastung definiert. Hier wird das Signal nicht über einen gewissen Zeitraum um den Abtastzeitpunkt akkumuliert, sondern exakt zum Abtastzeitpunkt nT ausgewertet.

Mathematisch lässt sich dies darstellen, indem man das Signal s(t) mit dem Dirac-Kamm, einer Folge von Dirac-Stößen, multipliziert:

Sampling einer Multiplikation mit Dirac-Kamm

Das abgetastete Signal s_a lautet dann

$s_{\mathrm {a} }(t)=s(t)\cdot \sum _{n=-\infty }^{\infty }\delta (t-nT)$ .

Für das Frequenzspektrum S_a, welches die Fourier-Reihe des Signals s_a darstellt, erhält man mit Hilfe der Umkehrung des Faltungstheorems:

$S_{{\mathrm {a}}}(f)=S(f)*\left[{\frac {1}{T}}\sum _{{n=-\infty }}^{{\infty }}\delta \left(f-{\frac {n}{T}}\right)\right]$ .

Das Spektrum ist also das Spektrum des Eingangssignals s(t), das periodisch mit der Periode 1/T wiederholt wird - dies drückt die Faltungseigenschaft des Dirac-Impulses aus. Daraus ergibt sich, dass das Spektrum von s maximal 1/(2T) breit sein darf, damit sich die verschobenen Spektren nicht überlappen.

Ist das Spektrum von s schmaler als 1/(2T), so ist das ursprüngliche Signal s(t) nach idealer Tiefpass-Filterung aus dem zeitdiskreten Spektrum vollständig rekonstruierbar. Dieser Umstand ist die Grundlage des Nyquist-Shannon-Abtasttheorems. Ist hingegen das Spektrum des Eingangssignals s(t) breiter als 1/(2T), tritt Aliasing auf und das ursprüngliche Signal s(t) kann aus s_a(t) nicht wiedergewonnen werden.

Reale Abtastung

Beispielhafter Signalverlauf f_a(t), in rot dargestellt, gewonnen durch ein Sample-and-Hold mit 0. Ordnung

In der Realität sind zwei Bedingungen der idealen Abtastung nicht einhaltbar:

Es ist nicht möglich, ideale Dirac-Stöße zu erzeugen. Das Signal wird vielmehr über einen Zeitraum um den eigentlichen Abtastzeitpunkt gewonnen. Dies nennt man auch natürliche Abtastung. Alternativ wird eine Sample-and-Hold-Schaltung eingesetzt.
Die perfekte Rückgewinnung des Signals aus seinem Spektrum erfordert zur Vermeidung von Spiegelspektren als Rekonstruktionsfilter einen idealen Tiefpass, welcher nicht kausal ist.

Die reale Abtastung erfolgt daher unter folgenden Modifikationen:

zu 1.: Der Dirac-Kamm wird durch eine Rechteckfunktion (rect) mit Rechteckimpulsen der Länge t₀ ersetzt. Die Abtastung wird durch eine Sample-and-Hold-Schaltung realisiert, welche den Wert einer Abtastung für die Länge des Rechteckimpulses konstant hält. Mathematisch entspricht dies einer Faltung mit der Rechteckfunktion:

$f_{\mathrm {a} }(t)=\sum _{n=-\infty }^{\infty }f(nT)\cdot \operatorname {rect} \left({\frac {t-nT}{t_{0}}}\right)=\operatorname {rect} \left({\frac {t}{t_{0}}}\right)*\sum _{n=-\infty }^{\infty }f(t)\delta (t-nT)$

Das daraus gewonnen Spektrum ist

$F_{\mathrm {a} }(f)=t_{0}\cdot \operatorname {si} (\pi ft_{0})\cdot \left[F(f)*1/T\sum _{n=-\infty }^{\infty }\delta (f-nf_{\mathrm {a} })\right]$ .

Dies ist das Spektrum der idealen Abtastung, gewichtet mit einem Faktor welcher die si-Funktion (Sinc-Funktion) beinhaltet. Dies stellt eine Verzerrung des Signals dar, welche durch eine zusätzliche Verzerrung im Rekonstruktionsfilter bei der Rückgewinnung des ursprünglichen Signals behoben werden kann. Diese Verzerrung tritt bei der natürlichen Abtastung nicht auf.

zu 2.: Um auch mit einem nicht-idealen Rekonstruktionsfilter das kontinuierliche Signal aus dem Spektrum mit möglichst kleinen Fehler zurückzugewinnen, kann die Abtastfrequenz erhöht werden. Durch die Überabtastung rücken anschaulich die Einzelspektren weiter auseinander, wodurch der Tiefpassfilter zur Rekonstruktion im Bereich der Spiegelspektren höhere Dämpfungswerte aufweist.

Abtastung im Spektralbereich

Aufgrund der Symmetrieeigenschaften der Fourier-Transformation lässt sich umgekehrt auch eine Frequenzfunktion S(f) im Spektralbereich, bei idealer Abtastung, durch eine Folge S_p(f) mit frequenzdiskreten Werten bilden:

$S_{\mathrm {p} }(f)=\sum _{n=-\infty }^{\infty }S(nF)\cdot \delta (f-nF)$

Die spektrale Folge S_p(f) besteht aus gewichteten Dirac-Impulsen, welche einzelne, diskrete Frequenzen beschreiben. Ein derartiges diskretes Spektrum wird auch Linienspektrum genannt.

Durch die inverse Fourier-Transformation kann daraus die zugehörige, periodische Form der Zeitfunktion s_p(t) gebildet werden:

$s_{\mathrm {p} }(t)={\frac {1}{T}}\sum _{n=-\infty }^{\infty }s\left(t-{\frac {n}{F}}\right)$

Auch bei der Abtastung im Spektralbereich gilt das Abtasttheorem in „umgekehrter“ Form: Wenn die zeitliche Dauer eines Signals s(t) kleiner als 1/F ist, dann überlappen sich die periodischen Anteile von s_p(t) nicht gegenseitig. Die Aufgabe des Rekonstruktionsfilters im Zeitbereich übernimmt eine Torschaltung, im einfachsten Fall ein Schalter, welcher für die Zeitdauer 1/F durchschaltet und die restliche Zeit sperrt. Ist hingegen das Signal s(t) länger als 1/F kommt es zu zeitlichen Überlappungen und die ursprüngliche Signalform lässt sich nicht mehr rekonstruieren.

Bei der realen Abtastung im Spektralbereich tritt an Stelle einer Folge von Dirac-Impulsen eine Folge spektraler Rechteckimpulse auf, welche jede für sich einen bandbegrenzten Ausschnitt aus dem Spektrum abdecken. Diese Funktion können im technischen Bezug Bandpassfilter übernehmen.

Allgemeine mathematische Darstellung der Abtastung

Bei der Speicherung eines Musikstückes auf einer CD wird das abgetastete Signal zur Übermittlung und Speicherung des analogen Ausgangssignals verwendet. Die zur Abtastung verwendete Methode ist in diesem Fall von der zur analogen Rekonstruktion verwendeten Methode abhängig. Diese Sichtweise ist auch für die mathematische Behandlung vorteilhaft.

Die Zusammensetzung aus Abtastung und Wiedergabe in umgekehrter Richtung tritt z.B. in der Nachrichtenübertragung auf, wenn eine binär kodierte Nachricht in ein analoges Funksignal umgesetzt wird. Durch einen Abtastprozess wird dann die ursprüngliche binäre Zeichenfolge rekonstruiert.

Umwandlung diskreter Daten in ein analoges Signal

Im einfachsten Fall wird die Umwandlung einer Folge reeller Zahlen, also eines zeitdiskreten Signals, mittels einer einzigen Kernfunktion vorgenommen. D.h., zu einer Folge $(c_{n})_{n\in \mathbb {Z} }$ wird mittels einer Funktion h und einem Zeitschritt T die, im weitesten Sinne, interpolierende Funktion

$a(t):=T\cdot \sum _{n\in \mathbb {Z} }c_{n}\cdot h(t-nT)$

gebildet. Deren Fouriertransformierte ist

$A(f):={\mathcal {F}}_{Hz}(a)(f)=T\sum _{n\in \mathbb {Z} }c_{n}\int _{\mathbb {R} }h(t-nT)e^{-i(2\pi f)t}\;dt=T\sum _{n\in \mathbb {Z} }c_{n}e^{-i(2\pi f)nT}H(f)$ ,

wobei H(f) die Fouriertransformierte von h(t) ist.

Umwandlung eines analogen Signals zu diskreten Daten

Ein realistischeres Modell der Messung eines zeitveränderlichen Prozesses ist die Bildung eines gewichteten Mittelwertes über einen bestimmten Zeitraum. Das kann mathematisch durch die Faltung mit einer Gewichtsfunktion w realisiert werden. Sei x(t) das zu messende Signal und v(t) der gemessene Wert zum Zeitpunkt t (der z.B. dem Schwerpunkt der Gewichtsfunktion zugeordnet wird), so gilt

$v(t)=(w*x)(t)=\int _{\mathbb {R} }w(s)\cdot x(t-s)\;ds$ .

Unter der Fouriertransformation geht die Faltung in die Multiplikation über. Seien W, V und X die Fouriertransformierten von w, v und x, dann gilt V(f)=W(f)X(f).

Der Prozess Signal-Daten-Signal

Bestimmt man nun eine Folge von Messwerten mit Zeitschritt T, $c_{n}=v(nT)$ , um diese in die Interpolationsvorschrift einzusetzen, so erhalten wir ein rekonstruktiertes analoges Signal

$a(t)=\sum _{n\in \mathbb {Z} }v(nT)\;h(t-nT)=\sum _{n\in \mathbb {Z} }(x*w)(nT)\;h(t-nT)$ .

Um den Fehler des gesamten Prozesses aus Diskretisierung und Wiedergabe einzuschätzen, kann man diesen Prozess auf einfache frequenzbeschränkte Testsignale anwenden. Dies kann im Modell durch die Bestimmung der Fouriertransformierte abgekürzt werden. Dazu ist jedoch die Fourierreihe

$T\sum _{n\in \mathbb {Z} }v(nT)e^{-i(2\pi f)nT}$

genauer zu bestimmen. Nach der Poissonschen Summenformel ist diese periodische Funktion identisch mit der Periodisierung von V(f). Sei $f_{\mathrm {s} }=1/T$ die Abtastfrequenz, dann gilt

$T\sum _{n\in \mathbb {Z} }v(nT)e^{-i(2\pi f)nT}=\sum _{k\in \mathbb {Z} }V(f+k\cdot f_{\mathrm {s} })=\sum _{k\in \mathbb {Z} }X(f+k\cdot f_{\mathrm {s} })\cdot W(f+k\cdot f_{\mathrm {s} }).$

Zusammenfassend gilt also

$A(f)=\left(\sum _{k\in \mathbb {Z} }X(f+k\cdot f_{\mathrm {s} })\cdot W(f+k\cdot f_{\mathrm {s} })\right)\,\cdot \,H(f).$

Eine Frequenzkomponente um die Frequenz f erleidet somit eine Verzerrung mit dem Faktor $W(f)\cdot H(f)$ und ein Aliasing der Stärke $W(f)H(f+kf_{\mathrm {s} })$ um die Frequenz $f+kf_{\mathrm {s} }$ bei $k\neq 0$ .

Um Basisbandsignale möglichst gut zu approximieren, ist es erforderlich, dass $W(f)\cdot H(f)\approx 1$ in einer Umgebung von f=0 und $H(f+kf_{\mathrm {s} })\approx 0$ für dieselben f und für alle $k\neq 0$ gelten. Im Rahmen einer mathematisch exakten Theorie sind diese Forderungen erfüllt und alle Operationen wohldefiniert, wenn

$\sum _{k\in \mathbb {Z} }|W(f+kf_{\mathrm {s} })|^{2}<B$ und $\sum _{k\in \mathbb {Z} }|H(f+kf_{\mathrm {s} })|^{2}<B$ mit einem B > 0 für alle f gilt und
$E(f):=\sum _{k\in \mathbb {Z} }|W(f)H(f+kf_{\mathrm {s} })-\delta _{0,k}|^{2}$ (mit dem Kronecker-Delta) in f=0 stetig ist und dort eine Nullstelle hat.

Man erhält dann in der Betragsquadratnorm des Funktionenraums L² für eine Basisbandfunktion x(t) mit höchster Frequenz $f_{H}$ nach der Parseval-Identität eine Abschätzung des relativen Fehlers, d.h. des Signal-Rausch-Verhältnisses, als

$\|a-x\|_{2}^{2}=\|A-X\|_{2}^{2}\leq \sup _{f\in [-f_{H},f_{H}]}E(f)\;\|x\|_{2}^{2}$

Beispiele: Sind w und h Rechteckfunktionen mit Breite T zentriert um 0, so ist

$W(f)=H(f)=\operatorname {sinc} (T\,f)={\frac {\sin(\pi T\,f)}{\pi T\,f}}$

und es gilt $E(f)=1-\operatorname {sinc} (T\,f)^{2}$ .

Sind umgekehrt w und h die Kardinalsinusfunktionen $\operatorname {sinc} (T\,f)$ , so sind deren Fouriertransformierte die entsprechenden Rechteckfunktionen $\operatorname {rect} (T/2\,f)$ , und die sich ergebende Rekonstruktionsformel ist die Kardinalreihe des Whittaker-Kotelnikow-Nyquist-Shannon-Abtasttheorems.

In jedem Fall führen Funktionen mit Frequenzkomponenten oberhalb ${\frac {1}{2}}f_{\mathrm {s} }$ zu Aliasfehlern im Frequenzbereich $\left[-{\frac {1}{2}}f_{\mathrm {s} },\,{\frac {1}{2}}f_{\mathrm {s} }\right]$ , die Frequenzschranke des Abtasttheorems ist also notwendig, aber keinesfalls hinreichend für eine fehlerarme Rekonstruktion.

Der Prozess Daten-Signal-Daten

In dieser Richtung wird die „interpolierende“ Funktion a(t) abgetastet. Es ergibt sich also

$d_{n}:=v(nT)=(w*a)(nt)=T\sum _{k\in \mathbb {Z} }c_{k}\int _{\mathbb {R} }w(s)h(nT-kT-s)\,ds=T\sum _{k\in \mathbb {Z} }c_{k}\;(w*h)(nT-kT)$ .

Für die Fourierreihen erhält man daraus

$\sum _{n\in \mathbb {Z} }d_{n}e^{-i(2\pi n)Tf}=T\sum _{k\in \mathbb {Z} }c_{k}e^{-i(2\pi k)Tf}\cdot \sum _{n\in \mathbb {Z} }(w*h)(nT)e^{-i(2\pi n)Tf}$

Nach der Poissonschen Summenformel gilt in diesem Fall

$T\sum _{k\in \mathbb {Z} }(w*h)(nT)e^{-i(2\pi n)Tf}=\sum _{n\in \mathbb {Z} }W(f-nf_{\mathrm {s} })\cdot H(f-nf_{\mathrm {s} })$ .

Soll die Fourierreihe der Folge c und damit die Folge erhalten bleiben, so muss diese Summe überall den Wert 1 haben. Das Maximum der Abweichung davon gibt auch in diesem Fall eine Schranke für den relativen Fehler bei der Datenübertragung.

$E(f)=\left|\sum _{n\in \mathbb {Z} }W(f-nf_{\mathrm {s} })\cdot H(f-nf_{\mathrm {s} })-1\right|$ und $\|d-c\|\leq \sup E(f)\;\|c\|$ .

Aus mathematischer Sicht müssen die Funktionen W und H wieder die oben angegebene Schranke der Periodisierung des Betragsquadrats einhalten.

Weitere Arten der Abtastung

Traditionell wird die äquidistante (periodische) Abtastung am meisten verwendet, weil sie schon sehr umfangreich untersucht und in vielen Anwendungen umgesetzt wurde. In letzten Jahrzehnten wurden auch andere Arten der Abtastung untersucht, die auf gleiche Zeitintervalle zwischen Abtastwerten verzichten, was einige Vorteile wie effektive Auslastung des Kommunikationskanals verspricht. Zu diesen gehört unter anderem Send-on-Delta-Abtastung.

Basierend auf einem Artikel in:

Wikipedia.de