Abtastwert

Abtastwert

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.

Die Anzahl der Abtastungen pro Sekunde wird Abtastrate genannt. Bei der digitalen Telefonie (ISDN) beträgt die Abtastrate beispielsweise 8 kHz.

Inhaltsverzeichnis

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:

  • Die Abtastung um ein zeitkontinuierliches Signal in ein zeitdiskretes Signal überzufü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 sa lautet dann

s_\mathrm{a} (t) = s(t) \cdot \sum_{n=-\infty}^{\infty} \delta(t - n T).

Für das Frequenzspektrum Sa, welches die Fouriertransformierte des Signals sa darstellt, erhält man:

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 der nicht abgetasteten Funktion sa(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 der Originalfunktion s(t) maximal 1/(2T) breit sein darf, damit sich die verschobenen Spektren nicht überlappen.

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

Reale Abtastung

Beispielhafter Signalverlauf fa(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:

  1. 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.
  2. 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 t0 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_a (t) = \sum_{n=-\infty}^{\infty} f(n T)\cdot \operatorname{rect}\left(\frac{t-n T}{t_0}\right) = \operatorname{rect}\left(\frac{t}{t_0}\right) * \sum_{n=-\infty}^{\infty} f(t) \delta(t - n T)

Das daraus gewonnen Spektrum ist

F_a (f) = t_0 \cdot \operatorname{si}(\pi f t_0) \cdot \left[F (f) * 1/T \sum_{n=-\infty}^{\infty} \delta(f - nf)\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 Sp(f) mit frequenzdiskreten Werten bilden:

S_\mathrm{p} (f) = \sum_{n=-\infty}^{\infty} S(n F) \cdot \delta(f - n F)

Die spektrale Folge Sp(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 sp(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 sp(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\Z} wird mittels einer Funktion h und einem Zeitschritt T die, im weitesten Sinne, interpolierende Funktion

a(t):=T\cdot \sum_{n\in\Z} c_n\cdot h(t-nT)

gebildet. Deren Fouriertransformierte ist


A(f):=\mathcal F_{Hz}(a)(f)
     =T\sum_{n\in\Z} c_n\int_\R h(t-nT)e^{-i(2\pi f)t}\;dt
     =T\sum_{n\in\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_\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, cn = v(nT), um diese in die Interpolationsvorschrift einzusetzen, so erhalten wir ein rekonstruktiertes analoges Signal

a(t)=\sum_{n\in\Z} v(nT)\;h(t-nT)=\sum_{n\in\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\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 fs = 1 / T die Abtastfrequenz, dann gilt


T\sum_{n\in\Z} v(nT) e^{-i(2\pi f)nT}
 =\sum_{k\in\Z}V(f+k\cdot f_s)
 =\sum_{k\in\Z}X(f+k\cdot f_s)\cdot W(f+k\cdot f_s).

Zusammenfassend gilt also

A(f)=\left(\sum_{k\in\Z}X(f+k\cdot f_s)\cdot W(f+k\cdot f_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 + kfs) um die Frequenz f + kfs bei k\ne 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_s)\approx 0 für dieselben f und für alle k\ne 0 gelten. Im Rahmen einer mathematisch exakten Theorie sind diese Forderungen erfüllt und alle Operationen wohldefiniert, wenn

  • \sum_{k\in\Z}|W(f+kf_s)|^2<B und \sum_{k\in\Z}|H(f+kf_s)|^2<B mit einem B>0 für alle f gilt und
  •  E(f):=\sum_{k\in\Z}|W(f)H(f+kf_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 für eine Basisbandfunktion x(t) mit höchster Frequenz fH 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\le \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 \frac12 f_s zu Aliasfehlern im Frequenzbereich \left[-\frac12 f_s,\,\frac12 f_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\Z}c_k \int_\R w(s) h(nT-kT-s)\,ds=T\sum_{k\in\Z}c_k \;(w*h)(nT-kT).

Für die Fourierreihen erhält man daraus


\sum_{n\in\Z}d_ne^{-i (2\pi n) Tf}
=T\sum_{k\in\Z}c_k e^{-i (2\pi k) Tf}\cdot\sum_{n\in\Z}(w*h)(nT)e^{-i (2\pi n) Tf}

Nach der Poissonschen Summenformel gilt in diesem Fall


T\sum_{k\in\Z} (w*h)(nT) e^{-i (2\pi n) Tf}=\sum_{n\in\Z} W(f-nf_s)\cdot H(f-nf_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\Z} W(f-nf_s)\cdot H(f-nf_s)-1\right| und \|d-c\|\le\sup E(f)\;\|c\|.

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

Literaturquellen

  • Hans Dieter Lüke: Signalübertragung. 6. Auflage. Springer Verlag, 1995, ISBN 3-540-58735-5. 

Weblinks


Wikimedia Foundation.

Игры ⚽ Нужен реферат?

Schlagen Sie auch in anderen Wörterbüchern nach:

  • Abtastwert — kvantuotasis dydis statusas T sritis automatika atitikmenys: angl. sampled data vok. Abtastwert, m rus. квантованная величина, f pranc. grandeur échantillonnée, f …   Automatikos terminų žodynas

  • .wv — WavPack Entwickler: David Bryant Aktuelle Version: 4.50.1 (3. Juli 2008) Betriebssystem: plattformunabhängig Kategorie: Encoder …   Deutsch Wikipedia

  • Adaptive Deltamodulation — Die adaptive Deltamodulation (ADM) ist ein Modulationsverfahren zur Umwandlung analoger in digitale Signale. Hierbei wird das analoge Signal in zeitgleichen Abständen abgetastet. Wie bei der gewöhnlichen, linearen Deltamodulation wird für jeden… …   Deutsch Wikipedia

  • WavPack — Entwickler David Bryant Aktuelle Version 4.60.1 (29. November 2009) Betriebssystem plattformunabhängig Kategorie Encoder Lizenz BSD Lize …   Deutsch Wikipedia

  • Wavpack — Entwickler: David Bryant Aktuelle Version: 4.50.1 (3. Juli 2008) Betriebssystem: plattformunabhängig Kategorie: Encoder …   Deutsch Wikipedia

  • Deltamodulation — Die Deltamodulation, auch als Δ Modulation oder Δ M bezeichnet, ist eine Variante der Differential Pulse Code Modulation (DPCM), bei welcher das pulscodierte Signal am Übertragungskanal nur die zwei Zustände 0 und 1 annehmen und durch ein Bit… …   Deutsch Wikipedia

  • LDM — Die Deltamodulation ist eine besondere Art der Pulscodemodulation. Hierbei wird das analoge Signal in gleichmäßigen Abständen abgetastet, je ein Abtastwert wird gespeichert und mit dem vorherigen verglichen. Ist die Amplitude des zweiten… …   Deutsch Wikipedia

  • Lineare Deltamodulation — Die Deltamodulation ist eine besondere Art der Pulscodemodulation. Hierbei wird das analoge Signal in gleichmäßigen Abständen abgetastet, je ein Abtastwert wird gespeichert und mit dem vorherigen verglichen. Ist die Amplitude des zweiten… …   Deutsch Wikipedia

  • Bildpunkt — Pixel, Bildpunkt, Bildzelle oder Bildelement (selten Pel) bezeichnet sowohl die kleinste Einheit einer digitalen Rastergrafik als auch deren Darstellung auf einem Bildschirm mit Rasteransteuerung. „Pixel“ (Nominativ Singular: das Pixel; Genitiv:… …   Deutsch Wikipedia

  • Digitalfilter — Ein digitales Filter ist ein elektronisches Filter welches ähnliche Aufgaben wie ein analoges Filter zu erfüllen hat: Es dient der Manipulation eines Signals wie beispielsweise das Sperren oder Durchlassen eines bestimmten Frequenzbereiches. Der… …   Deutsch Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”