Spektrale Leistungsdichte

Die spektrale Leistungsdichte gibt die auf die Frequenz bezogene Leistung eines Signals in einem infinitesimalen Frequenzband an. Diese Dichte besitzt die Dimension Leistung · Zeit, die Angabe erfolgt meist in den Einheiten Watt/Hertz oder dBm/Hz. Wird die spektrale Leistungsdichte über dem Frequenzspektrum angegeben, entsteht ein Leistungsdichtespektrum (LDS) oder Autoleistungsspektrum (engl.: Power-Spectral-Density (PSD), auch Wirkleistungsspektrum). Das Integral über alle Frequenzen ergibt die Gesamtleistung eines Signals. Während die Fouriertransformation von stationären Prozessen (wie z. B. Rauschen oder monofrequente Signale) unbeschränkt ist, lassen sich derartige Signale mit Hilfe der LDS quantitativ analysieren. Das LDS ist die Anzeigeform von Spektralanalysatoren, wobei hier die Leistung über vorgegebenen Frequenzintervallen (engl.: resolution bandwidth (RBW)) angegeben wird.

Allgemeines und Definition

Da für stationäre Prozesse $ f(t) $ im Allgemeinen weder die Energie $ \|f\|_{2}^{2} $ noch die Fouriertransformierte Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): F\big(f\big)(\omega) im klassischen Sinn existieren, liegt es nahe, zeitlich begrenzte Anteile Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): f_T(t) = f(t) für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): |t| \le T und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): 0 sonst zu betrachten. Nach der Formel von Plancherel gilt

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \frac{1}{2T}\int\limits_{\R} |f_T(t)|^2 \mathrm{d}t = \frac{1}{2T}\int\limits_{\R} |F(f_T)(\omega)|^2 \mathrm{d}\omega

Falls die mittlere Signalleistung

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r_{XX}(0) := \lim\limits_{T\to \infty} \frac{1}{2T}\int\limits_{-T}^T |f(t)|^2 \mathrm{d}t

existiert, existiert auch die rechte Seite obiger Formel und als spektrale Beschreibung der Leistung kann man die Spektrale Leistungsdichte definieren (falls der Grenzwert existiert) als

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): S_{XX}(\omega) := \lim\limits_{T\to \infty} \frac{1}{2T} |F(f_T)(\omega)|^2

Für jedes endliche Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): T heißt die Größe $ {\text{Per}}_{T}(\omega ):={\frac {1}{2T}}|F(f_{T})(\omega )|^{2} $ das Periodogramm von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): f . Es stellt einen Schätzwert der Spektralen Leistungsdichte dar, dessen Erwartungswert aber nicht Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): S_{XX}(\omega) entspricht (nicht erwartungstreu) und dessen Varianz auch für beliebig große Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): T nicht verschwindet (nicht konsistent).^[1]

Eigenschaften und Berechnung

Zur Bestimmung der spektralen Leistungsdichte Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): S_{XX}(\omega) wird oft das Wiener-Chintschin-Theorem herangezogen, wo sie über die Fouriertransformation der zeitlichen Autokorrelationsfunktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r_{xx}(t) des Signals gegeben wird:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): S_{XX}(\omega) = F(r_{xx})(\omega) = \frac{1}{2\pi} \int_{-\infty}^{\infty} r_{xx}(t) e^{-i\omega t} \mathrm{d}t

Dabei ist

$ r_{xx}(t)=\lim _{T\to \infty }{\frac {1}{2T}}\int _{-T}^{T}f(\tau )\;{\overline {f(t+\tau )}}\mathrm {d} \tau $

die Autokorrelationsfunktion des zeitlichen Signals Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): f(t) . Für Rauschsignale, allgemein für Prozesse, muss die Ergodizität vorausgesetzt werden, die es erlaubt Eigenschaften der Zufallsvariablen, wie den Erwartungswert, aus einer Musterfunktion zu bestimmen. In der Praxis kann nur ein endliches Zeitfenster betrachtet werden, weshalb man die Integrationsgrenzen einschränken muss. Nur für eine stationäre Verteilung ist die Korrelationsfunktion nicht mehr von der Zeit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): t abhängig.

Das Autoleistungsdichtespektrum ist gerade, reell und positiv. Dies bedeutet einen Informationsverlust, der eine Umkehrung dieser Prozedur verhindert.

Wird ein (Rausch-)Prozess mit Leistungsdichtespektrum Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): S_{XX}(\omega) über ein lineares, zeitinvariantes System mit Übertragungsfunktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): H(\omega) übertragen, so ergibt sich am Ausgang ein Leistungsdichtespektrum von

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): S_{YY}(\omega)= |H(\omega)|^2 \cdot S_{XX}(\omega).

Die Übertragungsfunktion geht quadratisch in die Formel ein, da das Spektrum eine Leistungsgröße ist.

Das Autoleistungsspektrum kann als einseitiges Spektrum Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): G_{XX} (f)(f \ge 0) dargestellt werden. Es gilt dann:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): G_{XX} = S_{XX} (f) \quad \mathrm{f\ddot ur} \quad f = 0

und

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): G_{XX} = 2 S_{XX} (f) \quad \mathrm{f\ddot ur} \quad f > 0\,.

Berechnungsmethoden beschränken sich üblicherweise auf bandbeschränkte Signale (Signale deren LDS für große Frequenzen verschwindet), die eine diskrete Darstellung erlauben (Nyquist-Shannon-Abtasttheorem). Erwartungstreue, konsistente Schätzwerte bandbegrenzter Signale, die auf einer Modifizierung des Periodogramms beruhen, sind beispielsweise die Welch-Methode oder Bartlett-Methode. Schätzungen auf Basis der Autokorrelationsfunktion heißen Korrelogramm-Verfahren wie beispielsweise die Blackmann-Tukey-Schätzung.^[2]

Anwendung und Einheiten

Die Kenntnis und Analyse der spektralen Leistungsdichte von Nutzsignal und Rauschen ist wesentlich zur Bestimmung des Signal-Rausch-Verhältnisses und zur Optimierung entsprechender Filter zur Rauschunterdrückung, zum Beispiel im Bildrauschen.

Das Autoleistungsspektrum kann für Aussagen über den Frequenzgehalt der analysierten Signale herangezogen werden. Spektralanalysatoren untersuchen die Spannung von Signalen. Für die Anzeige in Leistung ist die Angabe des Abschlusswiderstandes erforderlich. Mittels Spektralanalysatoren lässt sich aber die Spektralleistung nicht in einem infinitesimalen Frequenzband, sondern nur in einem Frequenzintervall endlicher Länge, bestimmen. Die so erhaltene spektrale Darstellung heißt Mean-Square-Spektrum (MSS) und ihre Wurzel RMS-Spektrum (engl. Root-Mean-Square). Die Länge des Frequenzintervalls ist stets mit angegeben und heißt Auflösebandbreite (engl. Resolution Bandwidth, kurz RBW oder BW) in der Einheit [Hz]. Die Umrechnung in Dezibel ist, wie für Leistungsangaben standardisiert, gemäß Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \text{MSS}_{dB} = 10 \log_{10}(\text{MSS}) , während die Umrechnung für RMS gemäß Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \text{RMS}_{dB} = 20 \log_{10}(\text{RMS}) erfolgt, womit die beiden Anzeigen in Dezibel zahlenmäßig identisch sind. Als Einheiten werden u. a. [ dBm ], [dBV], RMS-[V], PK-[V] (von engl. peak) verwendet. Die Angaben beziehen sich stets auf die verwendete Auflösebandbreite [Hz]. Beispielsweise erzeugt ein Sinussignal mit einem Spannungsverlauf von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): f(t) = 10 \sin(\omega t) V an einem Abschlusswiderstand von 50 Ohm eine effektive Spannung von 30 dBm oder 16,9897 dBV oder 7,0711 V (RMS) oder 10 V (PK) für jede Auflösebandbreite.

Beispiele

• Wenn die Korrelationsfunktion eine Delta-Distribution ist, spricht man von weißem Rauschen, in diesem Fall ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): S_{XX}(\omega) konstant.
• Für das thermische Rauschen, genauer die spektrale Rauschleistungsdichte, gilt: N₀ =k_B·T. Bei 27 °C beträgt es 4·10⁻²¹ J = 4·10⁻²¹ W/Hz = -204 dBW/Hz = -174 dBm/Hz
• Im Bild rechts ist ein MSS von der Funktion Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): f(t) = \sin(2\pi 3500 t) + 2^{-16} R(t) mit einem gleichverteilten Rauschprozess (Quantisierungsrauschen) Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): |R(t)|\le 1 bei einer Abtastrate von 44100 Hz und einer Auflösebandbreite von BW = 43,1 Hz (resultierend aus 44100 Hz / 1024 FFT Punkte) zu sehen, wie es beispielsweise von einer CD kommen könnte. Die Spitze bei etwa -3 dB repräsentiert das Sinussignal auf dem Rauschgrund bei etwa -128 dB. Da die Leistungsangaben sich auf die Auflösebandbreite beziehen, kann man das SNR zu $ -3-(-128+10\log _{10}(22050/{\text{BW}}))=97{,}9\,\mathrm {dB} $ ablesen (beachte das Logarithmusgesetz, das Multiplikationen in Additionen transformiert). Das aus dem Bild abgelesene SNR kommt damit dem theoretisch erwarteten von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): 10 \log_{10}(1^2/2) - 10 \log_{10}( (2^{-16})^2/3) = 98{,}0905\,\mathrm{dB} recht nahe.

Siehe auch

• Parsevalsche Gleichung
• Parsevalsches Theorem
• Spektrale Beschleunigungsdichte

Literatur

• Hans Dieter Lüke: Signalübertragung. Grundlagen der digitalen und analogen Nachrichtenübertragungssysteme. 6. neubearbeitete und erweiterte Auflage. Springer, Berlin u. a. 1995, ISBN 3-540-58753-5. 

Einzelnachweise