Difference between revisions of "Modulation Methods/Pulse Code Modulation"
Line 26: | Line 26: | ||
Nahezu alle heute eingesetzten Modulationsverfahren arbeiten digital. Deren Vorteile wurden schon im [[Modulation_Methods/Zielsetzung_von_Modulation_und_Demodulation#Vorteile_der_digitalen_Modulationsverfahren|ersten Kapitel]] dieses Buches genannt. Das erste Konzept zur digitalen Signalübertragung wurde bereits 1938 von [https://de.wikipedia.org/wiki/Alec_Reeves Alec Reeves] entwickelt und wird seit den 1960er Jahren unter dem Namen ''Pulscodemodulation'' (PCM) auch in der Praxis eingesetzt. Auch wenn sich viele der in den letzten Jahren konzipierten digitalen Modulationsverfahren von der PCM im Detail unterscheiden, so eignet sich diese doch sehr gut, um das Prinzip all dieser Verfahren zu erklären. | Nahezu alle heute eingesetzten Modulationsverfahren arbeiten digital. Deren Vorteile wurden schon im [[Modulation_Methods/Zielsetzung_von_Modulation_und_Demodulation#Vorteile_der_digitalen_Modulationsverfahren|ersten Kapitel]] dieses Buches genannt. Das erste Konzept zur digitalen Signalübertragung wurde bereits 1938 von [https://de.wikipedia.org/wiki/Alec_Reeves Alec Reeves] entwickelt und wird seit den 1960er Jahren unter dem Namen ''Pulscodemodulation'' (PCM) auch in der Praxis eingesetzt. Auch wenn sich viele der in den letzten Jahren konzipierten digitalen Modulationsverfahren von der PCM im Detail unterscheiden, so eignet sich diese doch sehr gut, um das Prinzip all dieser Verfahren zu erklären. | ||
− | [[File: | + | [[File:EN_Mod_T_4_1_S1_v2.png|right|frame|Principle of Pulse Code Modulation $\rm (PCM)]$]] |
+ | |||
Die Aufgabe des PCM–Systems ist es, | Die Aufgabe des PCM–Systems ist es, |
Revision as of 16:23, 25 October 2021
Contents
- 1 # ÜBERBLICK ZUM VIERTEN HAUPTKAPITEL #
- 2 Prinzip und Blockschaltbild
- 3 Abtastung und Signalrekonstruktion
- 4 Natürliche und diskrete Abtastung
- 5 Frequenzbereichsbetrachtung der natürlichen Abtastung
- 6 Frequenzbereichsbetrachtung der diskreten Abtastung
- 7 Quantization and quantization noise
- 8 PCM–Codierung und –Decodierung
- 9 Signal–zu–Rausch–Leistungsverhältnis
- 10 Einfluss von Übertragungsfehlern
- 11 Abschätzung der SNR-Degradation durch Übertragungsfehler
- 12 Nichtlineare Quantisierung
- 13 Kompression und Expandierung
- 14 Aufgaben zum Kapitel
# ÜBERBLICK ZUM VIERTEN HAUPTKAPITEL #
Das vierte Kapitel behandelt die digitalen Modulationsverfahren Amplitude Shift Keying (ASK), Phase Shift Keying (PSK) und Frequency Shift Keying (FSK) sowie einige davon abgeleitete Modifikationen. Die meisten der in den beiden letzten Kapiteln genannten Eigenschaften der analogen Modulationsverfahren gelten weiterhin. Unterschiede ergeben sich aus der nun erforderlichen Entscheiderkomponente des Empfängers.
Wir beschränken uns hier im wesentlichen auf die systemtheoretischen und übertragungstechnischen Aspekte. Die Fehlerwahrscheinlichkeit wird nur für ideale Bedingungen angegeben. Die Herleitungen und die Berücksichtigung nichtidealer Randbedingungen finden Sie im Buch „Digitalsignalübertragung”.
Im Einzelnen werden behandelt:
- die Pulscodemodulation (PCM) und deren Komponenten Abtastung – Quantisierung – Codierung,
- die linearen Modulationsverfahren ASK, BPSK und DPSK sowie zugehörige Demodulatoren,
- die Quadraturamplitudenmodulation (QAM) sowie kompliziertere Signalraumzuordnungen,
- die FSK – Frequency Shift Keying als Beispiel einer nichtlinearen digitalen Modulation,
- die FSK mit kontinuierlicher Phasenanpassung, insbesondere das (G)MSK–Verfahren.
Prinzip und Blockschaltbild
Nahezu alle heute eingesetzten Modulationsverfahren arbeiten digital. Deren Vorteile wurden schon im ersten Kapitel dieses Buches genannt. Das erste Konzept zur digitalen Signalübertragung wurde bereits 1938 von Alec Reeves entwickelt und wird seit den 1960er Jahren unter dem Namen Pulscodemodulation (PCM) auch in der Praxis eingesetzt. Auch wenn sich viele der in den letzten Jahren konzipierten digitalen Modulationsverfahren von der PCM im Detail unterscheiden, so eignet sich diese doch sehr gut, um das Prinzip all dieser Verfahren zu erklären.
Die Aufgabe des PCM–Systems ist es,
- das analoge Quellensignal q(t) in das Binärsignal qC(t) umzusetzen – diesen Vorgang bezeichnet man auch als A/D–Wandlung,
- dieses Signal über den Kanal zu übertragen, wobei das empfängerseitige Signal vC(t) wegen des Entscheiders ebenfalls binär ist,
- schließlich aus dem Binärsignal vC(t) das analoge sowie wert– und zeitkontinuierliche Sinkensignal v(t) zu rekonstruieren ⇒ D/A–Wandlung.
Weiterhin ist zum obigen PCM–Blockschaltbild anzumerken:
- Der PCM–Sender (bzw. der A/D–Wandler) setzt sich aus den drei Funktionsblöcken Abtastung – Quantisierung – PCM–Codierung zusammen, die in den nächsten Abschnitten noch im Detail beschrieben werden.
- Der grau hinterlegte Block zeigt das digitale Übertragungssystem mit digitalem Sender und Empfänger (letzterer beinhaltet auch einen Entscheider), sowie dem analogen Übertragungskanal, gekennzeichnet durch den Frequenzgang HK(f) und die Rauschleistungsdichte Φn(f).
- Dieser Block wird in den ersten drei Kapiteln des Buches Digital Signal Transmission eingehend behandelt. Im Kapitel 5 des gleichen Buches finden Sie auch digitale Kanalmodelle, die das Übertragungsverhalten anhand der Binärsignale qC(t) und vC(t) phänomenologisch beschreiben.
- Weiter erkennt man aus dem obigen Blockschaltbild, dass es für die Quantisierung empfängerseitig keine Entsprechung gibt. Deshalb wird sich auch bei fehlerfreier Übertragung, also für vC(t)=qC(t), das analoge Sinkensignal v(t) vom Quellensignal q(t) unterscheiden.
- Als Maß für die Qualität des (digitalen) Übertragungssystems verwenden wir das Signal-zu-Stör-Leistungsverhältnis ⇒ kurz: Sinken–SNR als der Quotient der Leistungen von Nutzsignal q(t) und Fehlersignal ε(t)=v(t)−q(t):
- ρv=PqPεmitPq=¯[q(t)]2,Pε=¯[v(t)−q(t)]2.
- Hierbei ist eine ideale Amplitudenanpassung vorausgesetzt, so dass im Idealfall (das heißt: Abtastung gemäß dem Abtasttheorem, bestmögliche Signalrekonstruktion, unendlich feine Quantisierung) das Sinkensignal v(t) mit dem Quellensignal q(t) exakt übereinstimmen würde.
Wir möchten Sie bereits hier auf das dreiteilige Lernvideo Pulscodemodulation hinweisen, das alle Aspekte der PCM beinhaltet. Dessen Prinzip wird im ersten Teil des Videos ausführlich erläutert.
Abtastung und Signalrekonstruktion
Die Abtastung – also die Zeitdiskretisierung des Analogsignals q(t) – wurde im Kapitel Zeitdiskrete Signaldarstellung des Buches „Signaldarstellung” ausführlich behandelt. Hier folgt eine Kurzzusammenfassung dieses Abschnitts.
Die Grafik verdeutlicht die Abtastung im Zeitbereich. Das (blaue) Signal q(t) ist zeitkontinuierlich, das im Abstand TA abgetastete (grüne) Signal qA(t) zeitdiskret. Dabei gilt:
- Die Abtastung lässt sich durch die Multiplikation des Analogsignals q(t) mit dem Diracpuls im Zeitbereich ⇒ pδ(t) darstellen:
- qA(t)=q(t)⋅pδ(t)mitpδ(t)=∞∑ν=−∞TA⋅δ(t−ν⋅TA).
- Das Gewicht der Diracfunktion bei t=ν·TA ist gleich TA·q(ν·TA). Da die Diracfunktion δ(t) die Einheit 1/s aufweist, hat somit qA(t) die gleiche Einheit wie q(t), zum Beispiel „V”.
- Die Fouriertransformierte des Diracpulses pδ(t) ist ebenfalls ein Diracpuls (aber nun im Frequenzbereich) ⇒ Pδ(f), wobei der Abstand der einzelnen Diraclinien fA=1/TA beträgt. Alle Impulsgewichte von Pδ(f) sind 1:
- pδ(t)=+∞∑ν=−∞TA⋅δ(t−ν⋅TA)∘−−−∙Pδ(f)=+∞∑μ=−∞δ(f−μ⋅fA).
- Das Spektrum QA(f) des abgetasteten Signals ergibt sich aus dem Faltungssatz, wobei Q(f) das kontinuierliche Spektrum des Analogsignals q(t) bezeichnet:
- QA(f)=Q(f)⋆Pδ(f)=+∞∑μ=−∞Q(f−μ⋅fA).
Wir weisen Sie hier auf den zweiten Teil des Lernvideos Pulscodemodulation hin, das die Abtastung und die Signalrekonstruktion systemtheoretisch erklärt.
Beispiel 1: Die obere Grafik zeigt schematisch das Spektrum Q(f) eines analogen Quellensignals q(t) mit Frequenzen bis fN, max=5 kHz.
- Tastet man q(t) mit der Abtastrate fA=20 kHz (also im jeweiligen Abstand TA=50 µs) ab, so erhält man das grün skizzierte periodische Spektrum QA(f).
- Da die Diracfunktionen unendlich schmal sind, beinhaltet qA(t) auch beliebig hochfrequente Anteile und dementsprechend ist QA(f) bis ins Unendliche ausgedehnt (mittlere Grafik).
- Darunter (rot) gezeichnet ist das Spektrum QA(f) für die Abtastparameter TA=100 µs ⇒ fA=10 kHz.
Fazit: Aus diesem Beispiel lassen sich folgende wichtige Erkenntnisse bezüglich der Abtastung gewinnen:
- Beinhaltet Q(f) Frequenzen bis fN, max, so muss nach dem Abtasttheorem die Abtastrate fA≥2·fN, max gewählt werden. Bei kleinerer Abtastrate fA (also größerem Abstand TA) kommt es zu Überlappungen der periodifizierten Spektren, also zu irreversiblen Verzerrungen.
- Gilt exakt fA=2·fN, max wie in der unteren Grafik von Beispiel 1, so kann Q(f) aus QA(f) – bzw. im PCM–System V(f) aus VQ(f) – durch einen idealen rechteckförmigen Tiefpass H(f) mit der Grenzfrequenz fG=fA/2 vollständig rekonstruiert werden.
- Erfolgt dagegen die Abtastung mit fA>2·fN, max wie in der mittleren Grafik des Beispiels, so kann empfängerseitig zur Signalrekonstruktion auch ein Tiefpass H(f) mit kleinerer Flankensteilheit verwendet werden, solange die folgende Bedingung erfüllt ist:
- H(f)={10f¨urf¨ur|f|≤fN,max,|f|≥fA−fN,max.
Natürliche und diskrete Abtastung
Die Multiplikation mit dem Diracpuls liefert nur eine idealisierte Beschreibung der Abtastung, da eine Diracfunktion (Dauer TR→0, Höhe 1/TR→∞) nicht realisierbar ist. In der Praxis muss der Diracpuls pδ(t) zum Beispiel durch einen Rechteckpuls
- pR(t)=+∞∑ν=−∞gR(t−ν⋅TA)mitgR(t)={11/20f¨urf¨urf¨ur|t|<TR/2,|t|=TR/2,|t|>TR/2
ersetzt werden, wobei die Rechteckimpulsdauer TR deutlich kleiner als der Abtastabstand TA sein sollte.
Die Grafik zeigt oben den Rechteckpuls pR(t). Darunter sind zwei verschiedene Abtastverfahren mit diesem Rechteckpuls dargestellt:
- Bei der natürlichen Abtastung ergibt sich das abgetastete Signal qA(t) durch die Multiplikation von q(t) mit pR(t). In den Bereichen pR(t)=1 hat somit qA(t) den gleichen Verlauf wie q(t).
- Bei der diskreten Abtastung wird das Signal q(t) – zumindest gedanklich – erst mit dem Diracpuls pδ(t) multipliziert. Dann wird jeder Diracimpuls TA·δ(t−ν·TA) durch einen Rechteckimpuls gR(t−ν·TA) ersetzt.
Hier und bei der nachfolgenden Frequenzbereichsbetrachtung ist zur Vereinfachung eine akausale Beschreibungsform gewählt. Für eine (kausale) Realisierung müsste gR(t)=1 im Bereich von 0 bis TR gelten, und nicht wie hier für −TR/2<t<TR/2.
Frequenzbereichsbetrachtung der natürlichen Abtastung
Definition: Die natürliche Abtastung lässt sich mit dem Faltungssatz im Spektralbereich wie folgt darstellen:
- qA(t)=pR(t)⋅q(t)=[1TA⋅pδ(t)⋆gR(t)]⋅q(t)⇒QA(f)=[Pδ(f)⋅1TA⋅GR(f)]⋆Q(f)=PR(f)⋆Q(f).
Die Grafik zeigt das Ergebnis für
- ein (unrealistisches) rechteckförmiges Spektrum Q(f)=Q0, das auf den Bereich |f|≤4 kHz begrenzt ist,
- die Abtastrate fA=10 kHz ⇒ TA=100 µs, sowie
- die Rechteckimpulsdauer TR=25 µs ⇒ TR/TA=0.25.
Man erkennt aus dieser Darstellung:
- Das Spektrum PR(f) bei natürlicher Abtastung ist im Gegensatz zu Pδ(f) kein Diracpuls (alle Gewichte gleich 1), sondern die Gewichte sind hier mit der Funktion GR(f)/TA=TR/TA·si(πfTR) bewertet.
- Wegen der Nullstelle der si–Funktion verschwinden hier die Diraclinien bei ±4fA.
- Das Spektrum QA(f) ergibt sich aus der Faltung mit Q(f). Das Rechteck um f=0 hat die Höhe TR/TA·Q0, die Anteile um μ·fA (μ≠0) sind weniger hoch.
- Verwendet man zur Signalrekonstruktion einen idealen, rechteckförmigen Tiefpass
- H(f)={TA/TR=40f¨urf¨ur|f|<fA/2,|f|>fA/2,
- so gilt für das Ausgangsspektrum V(f)=Q(f) und dementsprechend ist auch v(t)=q(t).
Fazit:
- Bei natürlicher Abtastung genügt zur Signalrekonstruktion wie bei idealer Abtastung (mit Diracpuls) ein Rechteck–Tiefpass.
- Allerdings muss zur Amplitudenanpassung im Durchlassbereich eine Verstärkung um den Faktor TA/TR berücksichtigt werden.
Frequenzbereichsbetrachtung der diskreten Abtastung
Definition: Bei der diskreten Abtastung erfolgt – zumindest gedanklich – zunächst die Multiplikation des Diracpulses pδ(t) mit dem Quellensignal q(t) und erst danach die Faltung mit dem Rechteckimpuls gR(t):
- qA(t)=[1/TA⋅pδ(t)⋅q(t)]⋆gR(t)⇒QA(f)=[Pδ(f)⋆Q(f)]⋅GR(f)/TA.
- Es ist unerheblich, aber durchaus zweckmäßig, dass hier der Faktor 1/TA zur Bewertungsfunktion GR(f) hinzugefügt wurde.
- Damit gilt wieder GR(f)/TA=TR/TA·si(πfTR).
Die obere Grafik zeigt (grün hinterlegt) die Spektralfunktion Pδ(f)⋆Q(f) nach idealer Abtastung. Bei diskreter Abtastung mit einem Rechteckpuls ergibt sich dagegen das Spektrum QA(f) entsprechend dem unteren Diagramm.
Man erkennt:
- Jedes der unendlich vielen Teilspektren hat nun eine andere Form. Wichtig ist nur das mittlere Spektrum um f=0.
- Alle anderen Spektralanteile werden empfängerseitig durch den Tiefpass der Signalrekonstruktion entfernt.
- Verwendet man für diesen Tiefpass wieder ein Rechteckfilter mit der Verstärkung um TA/TR im Durchlassbereich, so erhält man für das Ausgangsspektrum:
- V(f)=Q(f)⋅si(πfTR).
Fazit: Bei diskreter Abtastung und Rechteckfilterung kommt es zu Dämpfungsverzerrungen gmäß der Bewertungsfunktion si(πfTR).
- Diese sind um so stärker, je größer TR ist. Nur im Grenzfall TR→0 gilt si(πfTR)=1.
- Allerdings können durch eine ideale Entzerrung diese linearen Dämpfungsverzerrungen vollständig kompensiert werden.
- Um V(f)=Q(f) bzw. v(t)=q(t) zu erhalten, muss dann gelten:
- H(f)={(TA/TR)/si(πfTR)0f¨urf¨ur|f|<fA/2,|f|>fA/2
Quantization and quantization noise
Die zweite Funktionseinheit Quantisierung des PCM–Senders dient der Wertediskretisierung.
- Hierzu wird der gesamte Wertebereich des analogen Quellensignals (zum Beispiel der Bereich ±qmax) in M Intervalle aufgeteilt.
- Jedem Abtastwert qA(ν·TA) wird anschließend ein Repräsentant qQ(ν·TA) des zugehörigen Intervalls (beispielsweise die Intervallmitte) zugewiesen.
Beispiel 2: Die Grafik verdeutlicht die Quantisierung am Beispiel der Quantisierungsstufenzahl M=8.
- Tatsächlich wird für M in der Praxis wegen der anschließenden Binärcodierung stets eine Zweierpotenz gewählt.
- Jeder der durch Kreise markierten Abtastwerte qA(ν·TA) wird durch den dazugehörigen quantisierten Wert qQ(ν·TA) ersetzt. Die quantisierten Werte sind als Kreuze eingetragen.
- Dieser Vorgang der Wertdiskretisierung ist allerdings mit einer irreversiblen Verfälschung verbunden.
- Die Verfälschung εν=qQ(ν·TA) − qA(ν·TA) hängt von der Quantisierungsstufenzahl M ab. Es gilt folgende Schranke:
- |εν|<1/2⋅2/M⋅qmax=qmax/M.
Definition: Man bezeichnet den quadratischen Mittelwert der Fehlergröße εν als Quantisierungsrauschleistung:
- PQ=12N+1⋅+N∑ν=−Nε2ν≈1N⋅TA⋅∫N⋅TA0ε(t)2dtmitε(t)=qQ(t)−q(t).
Anmerkungen:
- Zur Berechnung der Quantisierungsrauschleistung PQ wird meist die angegebene Näherung der „Spontanquantisierung” verwendet.
- Man lässt dabei die Abtastung außer Betracht und bildet das Fehlersignal aus den zeitkontinuierlichen Signalen qQ(t) und q(t).
- PQ hängt auch vom Quellensignal q(t) ab. Unter der Voraussetzung, dass q(t) alle Werte zwischen ±qmax mit gleicher Wahrscheinlichkeit annimmt und der Quantisierer genau für diesen Bereich ausgelegt ist, ergibt sich entsprechend Aufgabe 4.4:
- PQ=q2max3⋅M2.
- Bei einem Sprach– oder Musiksignal können – wenn auch nur sehr selten – beliebig große Amplitudenwerte auftreten. In diesem Fall wird für qmax meist derjenige Amplitudenwert herangezogen, der nur zu 1% aller Zeiten (betragsmäßig) überschritten wird.
PCM–Codierung und –Decodierung
Der Block PCM–Codierung dient der Umsetzung der zeitdiskreten (nach Abtastung) und wertdiskreten (nach Quantisierung mit M Stufen) Signalwerte qQ(ν·TA) in eine Folge von N=log2(M) Binärwerte. Logarithmus zur Basis 2 ⇒ Logarithmus dualis.
Beispiel 3: Jeder Binärwert ⇒ Bit ist durch ein Rechteck der Dauer TB=TA/N dargestellt, woraus sich das Signal qC(t) ergibt.
Man erkennt:
- Es wird hier der Dualcode verwendet. Das bedeutet, dass die Quantisierungsintervalle μ von 0 bis M–1 durchnummeriert und anschließend in einfacher Binärform geschrieben werden. Mit M = 8 gilt beispielsweise \mu = 6 ⇔ 110.
- Die drei Binärsymbole des codierten Signals q_{\rm C}(t) ergeben sich, wenn man 0 durch L („Low”) und 1 durch H („High”) ersetzt. Im Beispiel erhält man so: HHL HHL LLH LHL HLH LHH.
- Die Bitdauer T_{\rm B} ist hier um den Faktor N = {\rm log_2}(M) = 3 kürzer als der Abtastabstand T_{\rm A} = 1/f_{\rm A}, und die Bitrate ist R_{\rm B} = {\rm log_2}(M) · f_{\rm A}.
- Verwendet man bei der Decodierung (v_{\rm C} ⇒ v_{\rm Q}) die gleiche Zuordnung wie bei der Codierung (q_{\rm Q} ⇒ q_{\rm C}), so gilt, falls es zu keinen Übertragungsfehlern kommt: v_{\rm Q}(ν · T_{\rm A}) = q_{\rm Q}(ν · T_{\rm A}).
- Eine Alternative zum Dualcode ist der Graycode, bei dem sich benachbarte Binärwerte nur in einem Bit unterscheiden. Für N = 3:
- \mu = 0: LLL, \mu = 1: LLH, \mu = 2: LHH, \mu = 3: LHL, \mu = 4: HHL, \mu = 5: HHH, \mu =6: HLH, \mu = 7: HLL.
Signal–zu–Rausch–Leistungsverhältnis
Das digitale Pulscodemodulation \rm (PCM) wird nun den analogen Modulationsverfahren \rm (AM, \ FM) hinsichtlich des erreichbaren Sinken–SNR ρ_v = P_q/P_ε bei AWGN–Rauschen vergleichend gegenüber gestellt.
Wie in vorherigen Kapiteln (zum Beispiel) bezeichnet ξ = {α_{\rm K} }^2 · P_{\rm S}/(N_0 · B_{\rm NF}) die Leistungskenngröße. Diese fasst verschiedene Einflüsse zusammen:
- den Kanalübertragungsfaktor α_{\rm K} (quadratisch),
- die Sendeleistung P_{\rm S},
- die AWGN–Rauschleistungsdichte N_0 (reziprok) sowie
- die Signalbandbreite B_{\rm NF} (ebenfalls reziprok);
bei einer harmonischen Schwingung: Frequenz f_{\rm N} statt B_{\rm NF}.
Die beiden Vergleichskurven für Amplitudenmodulation (AM) und für Frequenzmodulation (FM) lassen sich wie folgt beschreiben:
- Zweiseitenband–AM ohne Träger:
- ρ_v = ξ \ ⇒ \ 10 · \lg ρ_v = 10 · \lg \ ξ,
- Frequenzmodulation mit η = 3:
- ρ_υ = 3/2 \cdot η^2 · ξ = 13.5 · ξ \ ⇒ \ 10 · \lg \ ρ_v = 10 · \lg \ ξ + 11.3 \ \rm dB.
Die Kurve für das PCM 30/32–System ist wie folgt zu interpretieren:
- Ist die Leistungskenngröße ξ hinreichend groß, so treten keine Übertragungsfehler auf. Das Fehlersignal ε(t) = v(t) \ - \ q(t) ist dann allein auf die Quantisierung zurückzuführen (P_ε = P_{\rm Q}).
- Mit der Quantisierungsstufenzahl M = 2^N gilt in diesem Fall näherungsweise:
- \rho_{v} = \frac{P_q}{P_\varepsilon}= M^2 = 2^{2N} \hspace{0.3cm}\Rightarrow \hspace{0.3cm} 10 \cdot {\rm lg}\hspace{0.1cm}\rho_{v}=20 \cdot {\rm lg}\hspace{0.1cm}M = N \cdot 6.02\,{\rm dB}
- \Rightarrow \hspace{0.3cm} N = 8, \hspace{0.05cm} M =256\text{:}\hspace{0.2cm}10 \cdot {\rm lg}\hspace{0.1cm}\rho_{v}=48.16\,{\rm dB}\hspace{0.05cm}.
- Anzumerken ist, dass die angegebene Gleichung nur für ein sägezahnförmiges Quellensignal exakt gültig ist. Bei cosinusförmigem Quellensignal ist jedoch die Abweichung hiervon nicht sehr groß.
- Mit kleiner werdendem ξ (kleinere Sendeleistung oder größere Rauschleistungsdichte) nehmen die Übertragungsfehler zu. Damit wird P_ε > P_{\rm Q} und der Sinken–Störabstand wird kleiner.
- Die PCM (mit M = 256) ist den Analogverfahren (AM und FM) nur im unteren und mittleren ξ–Bereich überlegen. Spielen aber Übertragungsfehler keine Rolle mehr, so ist durch ein größeres ξ keine Verbesserung mehr zu erzielen (horizontaler, gelb hinterlegter Kurvenabschnitt).
- Eine Verbesserung bringt nur eine Erhöhung von N (Bitanzahl pro Abtastwert) ⇒ größeres M = 2^N (Quantisierungsstufenzahl). Beispielsweise erreicht man bei einer Compact Disc (CD) mit dem Parameter N = 16 ⇒ M = 65536 den Wert
- 10 · \lg \ ρ_v = 96.32 \ \rm dB.
\text{Beispiel 4:} Die folgende Grafik zeigt den begrenzenden Einfluss der Quantisierung:
- Weiß gepunktet dargestellt ist das Quellensignal q(t), grün gepunktet das Sinkensignal v(t) nach PCM mit N = 4 ⇒ M = 16.
- Die Abtastzeitpunkte sind durch Kreuze markiert.
- Übertragungsfehler werden vorerst ausgeschlossen. Abtastung und Signalrekonstruktion seien bestmöglich an q(t) angepasst.
Dieses Bild kann wie folgt interpretiert werden:
- Mit N = 8 ⇒ M = 256 ist das Sinkensignal v(t) vom Quellensignal q(t) mit dem bloßen Auge nicht zu unterscheiden. Für beide gilt näherungsweise der weiß gepunktete Signalverlauf.
- Am Störabstand 10 · \lg \ ρ_v = 47.8 \ \rm dB erkennt man aber, dass das Quantisierungsrauschen (Leistung P_\varepsilon des Fehlersignals) nur etwa um den Faktor 1.6 · 10^{–5} kleiner ist als die Leistung P_q des Quellensignals. Dieses SNR wäre bei einem Sprach– oder Musiksignal schon deutlich hörbar.
- Obwohl das hier betrachtete Quellensignal weder sägezahnförmig noch cosinusförmig verläuft, sondern sich aus mehreren Frequenzanteilen zusammensetzt, weicht die angegebene Näherung ρ_v ≈ M^2 ⇒ 10 · \lg \ ρ_υ = 48.16 \ \rm dB nur unwesentlich vom tatsächlichen Wert ab.
- Dagegen erkennt man für N = 4 ⇒ M = 16 bereits im Bild Abweichungen zwischen dem Sinkensignal (grün markiert) und dem Quellensignal (weiße Markierung), was auch durch den sehr kleinen Störabstand 10 · \lg \ ρ_υ = 28.2 \ \rm dB quantitativ zum Ausdruck kommt.
Einfluss von Übertragungsfehlern
Ausgehend vom gleichen Analogsignal q(t) wie im letzten Abschnitt und einer linearen Quantisierung mit N = 8 Bit ⇒ M = 256 werden nun die Auswirkungen von Übertragungsfehlern anhand des jeweiligen Sinkensignals v(t) verdeutlicht.
- Die weißen Punkte markieren wieder das Quellensignal q(t). Ohne Übertragungsfehler hat das Sinkensignal v(t) bei Vernachlässigung der Quantisierung den gleichen Verlauf.
- Nun wird jeweils genau ein Bit des fünften Abtastwertes q(5 · T_{\rm A}) = -0.715 verfälscht, wobei dieser Abtastwert mit LLHL LHLL codiert wurde. Dieser Grafik zugrunde liegt der Dualcode, das heißt, dass das unterste Quantisierungsintervall (\mu = 0) mit LLLL LLLL und das oberste Intervall (\mu = 255) mit HHHH HHHH dargestellt wird.
Die Tabelle zeigt die Ergebnisse dieser Untersuchung:
- Der angegebene Störabstand 10 · \lg \ ρ_v wurde aus dem dargestellten (sehr kurzen) Signalausschnitt der Dauer 10 · T_{\rm A} berechnet.
- Bei jeweils einem Fehler bei der Übertragung von 10 · 8 = 80 Bit entspricht dies einer Bitfehlerrate von 1.25\%.
Die in der Grafik und der Tabelle dargestellten Ergebnisse dieser Fehleranalyze können wie folgt zusammengefasst werden:
- Wird nur das letzte Bit des Binärwortes verfälscht (LSB: Least Significant Bit, LLHL LHLL ⇒ LLHL LHLH), so ist mit bloßem Auge kein Unterschied zur fehlerfreien Übertragung zu erkennen (weißer Kurvenzug). Trotzdem wird der Störabstand um 3.5 \ \rm dB vermindert.
- Ein Übertragungsfehler des viertletzten Bits (grüne Kurve, LLHLLHLL ⇒ LLHLHHLL) führt bereits zu einer deutlich erkennbaren Verfälschung um acht Quantisierungsintervalle. Das heißt: v(5T_{\rm A}) \ - \ q(5T_{\rm A}) = 8/256 · 2 = 0.0625 und der Störabstand sinkt auf 10 · \lg \ ρ_υ = 28.2 \ \rm dB.
- Die rote Kurve zeigt schließlich den Fall, dass das MSB (Most Significant Bit) verfälscht wird: LLHLLHLL ⇒ HLHLLHLL. Dies führt zur Verfälschung v(5T_{\rm A}) \ – \ q(5T_{\rm A}) = 1 (entspricht dem halben Aussteuerbereich). Der Störabstand beträgt nun nur mehr etwa 4 \ \rm dB.
- Zu allen Abtastzeitpunkten mit Ausnahme von 5T_{\rm A} stimmt v(t) bis auf den Quantisierungsfehler mit q(t) exakt überein. Außerhalb dieser durch gelbe Kreuze markierten Zeitpunkte führt der einzige Fehler bei 5T_{\rm A} aber in einem ausgedehnten Bereich zu starken Abweichungen, was auf die Interpolation mit der \rm si–förmigen Impulsantwort des Rekonstruktionstiefpasses H(f) zurückzuführen ist.
Abschätzung der SNR-Degradation durch Übertragungsfehler
Nun soll versucht werden, die SNR–Kurve des PCM–Systems unter Berücksichtigung von Bitfehlern zumindest näherungsweise zu bestimmen. Wir gehen dabei vom folgenden Blockschaltbild aus und setzen weiter voraus:
- Jeder Abtastwert q_{\rm A}(νT) wird mit M Stufen quantisiert und mit N = {\rm log_2} (M) Binärzeichen (Bit) dargestellt. Im Beispiel gilt M = 8 ⇒ N = 3.
- Die Binärdarstellung von q_{\rm Q}(νT) liefert die Amplitudenkoeffizienten a_k\, (k = 1, \text{...} \hspace{0.08cm}, N), die durch Bitfehler in die Koeffizienten b_k verfälscht werden können.
- Sowohl a_k als auch b_k sind jeweils ±1.
- Ein Bitfehler (b_k ≠ a_k) tritt mit der Wahrscheinlichkeit p_{\rm B} auf.
- Jedes Bit wird gleichwahrscheinlich verfälscht und in jedem PCM–Wort ist maximal ein Fehler ⇒ nur eines der N Bit kann falsch sein.
Aus dem in der Grafik angegebenen Diagramm ist für N = 3 und natürliche Binärcodierung (Dualcode) zu erkennen:
- Eine Verfälschung von a_1 verändert den quantisierten Wert q_{\rm Q}(νT) um ±A.
- Eine Verfälschung von a_2 verändert den quantisierten Wert q_{\rm Q}(νT) um ±A/2.
- Eine Verfälschung von a_3 verändert den quantisierten Wert Wert q_{\rm Q}(νT) um ±A/4.
Durch Verallgemeinerung erhält man für die Abweichung ε_k = υ_{\rm Q}(νT) \ - \ q_{\rm Q}(νT) für den Fall, dass der Amplitudenkoeffizient a_k falsch übertragen wurde:
- \varepsilon_k = - a_k \cdot A \cdot 2^{-k +1} \hspace{0.05cm}.
Für die Fehlerrauschleistung erhält man nach Mittelung über alle Verfälschungswerte ε_k (mit 1 ≤ k ≤ N) unter Berücksichtigung der Bitfehlerwahrscheinlichkeit p_{\rm B}:
- P_{\rm F}= {\rm E}\big[\varepsilon_k^2 \big] = \sum\limits^{N}_{k = 1} p_{\rm B} \cdot \left ( - a_k \cdot A \cdot 2^{-k +1} \right )^2 =\ p_{\rm B} \cdot A^2 \cdot \sum\limits^{N-1}_{k = 0} 2^{-2k } = p_{\rm B} \cdot A^2 \cdot \frac{1- 2^{-2N }}{1- 2^{-2 }} \approx {4}/{3} \cdot p_{\rm B} \cdot A^2 \hspace{0.05cm}.
- Hierbei ist die Summenformel der geometrischen Reihe sowie die Näherung 1 – 2^{–2N } ≈ 1 verwendet.
- Für N = 8 ⇒ M = 256 beträgt der damit verbundene relative Fehler beispielsweise etwa \rm 10^{–5}.
Ohne Berücksichtigung von Übertragungsfehlern hat sich für das Signal–zu–Rausch–Leistungsverhältnis ρ_v = P_{\rm S}/P_{\rm Q} ergeben, wobei bei einem gleichverteilten Quellensignal (zum Beispiel sägezahnförmig) die Signalleistung und die Quantisierungsrauschleistung wie folgt zu berechnen ist:
- P_{\rm S}={A^2}/{3}\hspace{0.05cm},\hspace{0.3cm}P_{\rm Q}= {A^2}/{3} \cdot 2^{-2N } \hspace{0.05cm}.
Unter Berücksichtigung der Übertragungsfehler erhält man mit obigem Ergebnis:
- \rho_{\upsilon}= \frac{P_{\rm S}}{P_{\rm Q}+P_{\rm F}} = \frac{A^2/3}{A^2/3 \cdot 2^{-2N } + A^2/3 \cdot 4 \cdot p_{\rm B}} = \frac{1}{ 2^{-2N } + 4 \cdot p_{\rm B}} \hspace{0.05cm}.
Die Grafik zeigt 10 · \lg ρ_v in Abhängigkeit der (logarithmierten) Leistungskenngröße ξ = P_{\rm S}/(N_0 · B_{\rm NF}), wobei B_{\rm NF} die Signalbandbreite angibt. Der konstante Kanalübertragungsfaktor sei idealerweise α_{\rm K} = 1.
- Beim optimalen Binärsystem und AWGN–Rauschen gilt aber für die Leistungskenngröße auch ξ = E_{\rm B}/N_0 (Energie pro Bit bezogen auf die Rauschleistungsdichte).
- Die Bitfehlerwahrscheinlichkeit ist dann mit der Gaußschen Fehlerfunktion {\rm Q}(x) wie folgt gegeben:
- p_{\rm B}= {\rm Q} \left ( \sqrt{{2E_{\rm B}}/{N_0} }\right ) \hspace{0.05cm}.
- Für N = 8 ⇒ 2^{–2{\it N} } = 1.5 · 10^{–5} sowie 10 · \lg \ ξ = 6 \ \rm dB ⇒ p_{\rm B} = 0.0024 (rot markierter Punkt) ergibt sich:
- \rho_{\upsilon}= \frac{1}{ 1.5 \cdot 10^{-5} + 4 \cdot 0.0024} \approx 100 \hspace{0.3cm} \Rightarrow \hspace{0.3cm}10 \cdot {\rm lg} \hspace{0.15cm}\rho_{\upsilon}\approx 20\,{\rm dB} \hspace{0.05cm}.
- Dieser kleine ρ_v–Wert geht auf den Term 4 · 0.0024 im Nenner (Einfluss des Übertragungsfehlers) zurück, während im horizontalen Kurvenabschnitt für jedes N (Bitanzahl pro Abtastwert) der Term \rm 2^{–2{\it N} } dominiert – also das Quantisierungsrauschen.
Nichtlineare Quantisierung
Häufig werden die Quantisierungsintervalle nicht gleich groß gewählt, sondern man verwendet für den inneren Amplitudenbereich eine feinere Quantisierung als für große Amplituden. Dafür gibt es mehrere Gründe:
- Bei Audiosignalen werden Verfälschungen der leisen Signalanteile (also Werte in der Nähe der Nulllinie) subjektiv als störender empfunden als eine Beeinträchtigung großer Amplitudenwerte.
- Eine solche ungleichmäßige Quantisierung führt bei einem solchen Musik– oder Sprachsignal auch zu einem größeren Sinkenstörabstand, da hier die Signalamplitude nicht gleichverteilt ist.
Die Grafik zeigt ein Sprachsignal q(t) und dessen Amplitudenverteilung f_q(q) ⇒ Wahrscheinlichkeitsdichtefunktion.
Es handelt sich um die Laplaceverteilung, die man wie folgt annähern kann:
- durch eine kontinuierliche, zweiseitige Exponentialverteilung, und
- durch eine Diracfunktion δ(q) zur Berücksichtigung der Sprachpausen (magentafarben).
In der Grafik ist die nichtlineare Quantisierung nur angedeutet, zum Beispiel mittels der 13–Segment–Kennlinie, die in der Aufgabe 4.5 genauer beschrieben ist:
- Die Quantisierungsintervalle werden hierbei zu den Rändern hin abschnittsweise immer breiter.
- Die häufigeren kleinen Amplituden werden dagegen sehr fein quantisiert.
Kompression und Expandierung
Eine ungleichmäßige Quantisierung kann zum Beispiel dadurch realisiert werden, in dem
- die abgetasteten Werte q_{\rm A}(ν · T_{\rm A}) zunächst durch eine nichtlineare Kennlinie q_{\rm K}(q_{\rm A}) verformt und
- anschließend die entstehenden Ausgangswerte q_{\rm K}(ν · T_{\rm A}) gleichmäßig quantisiert werden.
Damit ergibt sich die nebenstehend skizzierte Signalkette.
\text{Fazit:} Eine solche ungleichmäßige Quantisierung bedeutet:
- Durch die nichtlineare Kennlinie q_{\rm K}(q_{\rm A}) werden kleine Signalwerte verstärkt und große Werte abgeschwächt ⇒ Kompression.
- Diese bewusste Signalverzerrung macht man beim Empfänger durch die Umkehrfunktion v_{\rm E}(υ_{\rm Q}) rückgängig ⇒ Expandierung.
- Den Gesamtvorgang von sendeseitiger Kompression und empfängerseitiger Expansion nennt man auch Kompandierung.
Für das PCM–System 30/32 wurde von der Comité Consultatif International des Télégraphique et Téléphonique (CCITT) die so genannte A–Kennlinie empfohlen:
- y(x) = \left\{ \begin{array}{l} \frac{1 + {\rm ln}(A \cdot x)}{1 + {\rm ln}(A)} \\ \frac{A \cdot x}{1 + {\rm ln}(A)} \\ - \frac{1 + {\rm ln}( - A \cdot x)}{1 + {\rm ln}(A)} \\ \end{array} \right.\quad\begin{array}{*{5}c}{\rm{f\ddot{u}r}}\\{\rm{f\ddot{u}r}}\\{\rm{f\ddot{u}r}} \\ \end{array}\begin{array}{*{10}c}1/A \le x \le 1\hspace{0.05cm}, \\ - 1/A \le x \le 1/A\hspace{0.05cm}, \\ - 1 \le x \le - 1/A\hspace{0.05cm}. \\ \end{array}
- Hierbei ist zur Abkürzung x = q_{\rm A}(ν · T_{\rm A}) und y = q_{\rm K}(ν · T_{\rm A}) verwendet.
- Diese Kennlinie mit dem in der Praxis eingeführten Wert A = 87.56 hat eine sich ständig ändernde Steigung.
- Nähere Angaben zu dieser Art der ungleichmäßigen Quantisierung finden Sie in der Aufgabe 4.5.
Hinweis: Im dritten Teil des Lernvideos Pulscodemodulation werden behandelt:
- die Definition des Signal–zu–Rausch–Leistungsverhältnisses (SNR),
- der Einfluss von Quantisierungsrauschen und Übertragungsfehlern,
- die Unterschiede zwischen linearer und nichtlinearer Quantisierung.
Weitere Informationen zum Thema sowie Aufgaben, Simulationen und Programmierübungen finden Sie im
- Kapitel 12: Pulscodemodulation, Programm pcm
des Praktikums "Simulationsmethoden in der Nachrichtentechnik". Diese (ehemalige) LNT-Lehrveranstaltung an der TU München basiert auf
- dem Lehrsoftwarepaket LNTsim ⇒ Link verweist auf die ZIP-Version des Programms;
- dieser Praktikumsanleitung ⇒ Link verweist auf die PDF-Version; Kapitel 12: Seite 271-294.
Aufgaben zum Kapitel
Aufgabe 4.2: Tiefpass zur Signalrekonstruktion
Aufgabe 4.2Z: Zum Abtasttheorem
Aufgabe 4.3: Natürliche und diskrete Abtastung
Aufgabe 4.4: Zum Quantisierungsrauschen
Aufgabe 4.4Z: Störabstand bei PCM
Aufgabe 4.5: Nichtlineare Quantisierung
Aufgabe 4.6: Quantisierungskennlinien