Quadrature Amplitude Modulation

1 Allgemeine Beschreibung und Signalraumzuordnung (1)
2 Allgemeine Beschreibung und Signalraumzuordnung (2)
3 Systembeschreibung durch das äquivalente TP–Signal
4 Leistung und Energie von QAM–Signalen
5 Signalverläufe der 4–QAM
6 Fehlerwahrscheinlichkeit der 4–QAM
7 QAM–Signalraumkonstellationen
8 Weitere Signalraumkonstellationen
9 Nyquist– und Wurzel–Nyquist–QAM–Systeme
10 Quellenverzeichnis

Allgemeine Beschreibung und Signalraumzuordnung (1)

Aufgrund der Orthogonalität von Cosinus und (Minus–)Sinus kann man über einen Übertragungskanal zwei Datenströme unabhängig voneinander übertragen. Die Grafik zeigt das allgemeine Blockschaltbild.

Dieses sehr allgemeine Modell lässt sich wie folgt beschreiben:

Am Eingang liegt die binäre Quellensymbolfolge $〈q_k〉$ mit der Bitrate $R_{\rm B}$ an. Der zeitliche Abstand zweier Symbole ist damit $T_{\rm B} = 1/R_{\rm B}$.
Aus jeweils $b$ binären Eingangssymbolen $q_k$ werden zwei mehrstufige Amplitudenkoeffizienten $a_{\rm Iν}$ und $a_{\rm Qν}$ abgeleitet, wobei „I” für Inphase und „Q” für Quadraturkomponente steht.
Ist $b$ geradzahlig und die Signalraumzuordnung quadratisch, so können die Koeffizienten $a_{\rm Iν}$ und $a_{\rm Qν}$ jeweils einen von $M = 2^{b/2}$ Amplitudenwerten mit gleicher Wahrscheinlichkeit annehmen. Man spricht dann von Quadratur–Amplitudenmodulation (QAM).
Das in der Grafik betrachtete Beispiel gilt für die 16–QAM mit $b = M =$ 4 und dementsprechend 16 Signalraumpunkten. Bei einer 256–QAM würde $b =$ 8 und $M =$ 16 gelten $(2^b = M^2 =$ 256).

Allgemeine Beschreibung und Signalraumzuordnung (2)

Fortsetzung der Bildbeschreibung zur obigen Grafik:

Anschließend werden die Koeffizienten $a_{\rm Iν}$ und $a_{\rm Qν}$ jeweils einem Diracpuls als Impulsgewichte eingeprägt. Nach der Impulsformung mit dem Sendegrundimpuls $g_s(t)$ gilt somit in den beiden Zweigen des Blockschaltbildes:

$$\begin{align*}s_{\rm I}(t) & = \sum_{\nu = - \infty}^{+\infty}a_{\rm I\hspace{0.03cm}\it \nu} \cdot g_s (t - \nu \cdot T)\hspace{0.05cm},\\s_{\rm Q}(t) & = \sum_{\nu = - \infty}^{+\infty}a_{\rm Q\hspace{0.03cm}\it \nu} \cdot g_s (t - \nu \cdot T)\hspace{0.05cm}.\end{align*}$$

Anzumerken ist, dass wegen der redundanzfreien Umsetzung die Symboldauer $T$ dieser Signale um den Faktor $b$ größer ist als die Bitdauer $T_{\rm B}$ des binären Quellensignals. Im gezeichneten Beispiel (16-QAM) gilt $T = 4 · T_{\rm B}$.
Das QAM–Sendesignal $s(t)$ ist dann die Summe der beiden mit Cosinus bzw. Minus–Sinus multiplizierten Teilsignale:

$$\begin{align*}s_{\rm cos}(t) & = s_{\rm I}(t) \cdot \cos(2 \pi f_{\rm T} t),\\ s_{\rm -sin}(t) & = -s_{\rm Q}(t) \cdot \sin(2 \pi f_{\rm T} t)\end{align*}$$

Die beiden Übertragungszweige (I, Q) können als zwei völlig getrennte $M$–stufige ASK–Systeme aufgefasst werden, die sich gegenseitig nicht stören, solange alle Komponenten optimal ausgelegt sind. Die Quadratur–Amplitudenmodulation ermöglicht somit (im Idealfall) eine Verdoppelung der Datenrate bei gleichbleibender Qualität.

Systembeschreibung durch das äquivalente TP–Signal

Da die Multiplikation von $s_{\rm I}(t)$ und $s_{\rm Q}(t)$ mit einer Cosinus– bzw. Minus–Sinus–Schwingung nur eine Verschiebung im Frequenzbereich bewirkt und eine solche Verschiebung eine lineare Operation darstellt, lässt sich die Systembeschreibung mit Hilfe der äquivalenten TP–Signale wesentlich vereinfachen.

Die Grafik zeigt das vereinfachte Modell im Basisband. Dieses ist äquivalent zum bisher betrachteten Blockschaltbild. Beachten Sie bitte die folgenden Hinweise:

Die im Blockschaltbild rot gezeichnete Seriell–Parallel–Wandlung und die Signalraumzuordnung bleibt erhalten, obwohl dieser Block hier nicht mehr eingezeichnet ist. Lassen wir zunächst auch den oft aus schaltungstechnischen Gründen eingebrachten Bandpass $H_{\rm BP}(f)$ außer Betracht.
Alle Doppelpfeile in dem obigen Basisbandmodell kennzeichnen komplexe Größen. Die damit verbundenen Operationen sind ebenfalls komplex zu verstehen. Beispielsweise fasst der komplexe Amplitudenkoeffizient $a_ν$ je einen Inphase– und einen Quadraturkoeffizienten zusammen:

$$a_\nu = a_{\rm I\hspace{0.03cm}\it \nu} + {\rm j} \cdot a_{\rm Q\hspace{0.03cm}\it \nu} \hspace{0.05cm}.$$

Die äquivalente Tiefpass–Repräsentation des tatsächlichen, physikalischen und damit per se reellen Sendesignals $s(t)$ ist bei QAM stets komplex und es gilt mit den Teilsignalen $s_{\rm I}(t)$ und $s_{\rm Q}(t)$:

$$s_{\rm TP}(t) = s_{\rm I}(t) + {\rm j} \cdot s_{\rm Q}(t) = \sum_{\nu = - \infty}^{+\infty} a_\nu \cdot g_s (t - \nu \cdot T)\hspace{0.05cm}.$$

Zum analytischen Signal $s_+(t)$ kommt man von $s_{\rm TP}(t)$ durch Multiplikation mit der komplexen Exponentialfunktion. Das physikalische Sendesignal $s(t)$ ergibt sich dann als der Realteil von $s_+(t)$.
Damit die Vorzeichen im Blockschaltbild der letzten Seite und im skizzierten Basisbandmodell übereinstimmen, ist im Quadraturzweig die Multiplikation mit der negativen Sinus–Schwingung erforderlich, wie die nachfolgende Rechnung zeigt:

$$\begin{align*}s(t) & = {\rm Re}[s_{\rm +}(t)] = {\rm Re}[s_{\rm TP}(t) \cdot{\rm e}^{{\rm j}2\pi f_{\rm T} t}] = \\ & = {\rm Re} \left[\left ( \sum (a_{\rm I\hspace{0.03cm}\it \nu} + {\rm j} \cdot a_{\rm Q\hspace{0.03cm}\it \nu} ) \cdot g_s (t - \nu \cdot T)\right )\left ( \cos(2 \pi f_{\rm T} t) + {\rm j} \cdot \sin(2 \pi f_{\rm T} t) \right )\right]= \\ & = s_{\rm I}(t) \cdot \cos(2\pi f_{\rm T} t) - s_{\rm Q}(t) \cdot \sin(2 \pi f_{\rm T} t) \hspace{0.05cm}.\end{align*}$$

Der Einfluss des Bandpasses $H_{\rm BP}(f)$, der in der Praxis oft am Ausgang des QAM–Modulators zu berücksichtigen ist, kann dem Impulsformfilter $g_s(t)$ beaufschlagt werden. Ist der Durchlassbereich des BP–Filters symmetrisch um $f_{\rm T}$, so ist sein Tiefpass–Äquivalent (im Zeitbereich) $h_{\rm BP→TP}(t)$ rein reell und man kann im Modell $g_s(t)$ durch $g_s(t) \star h_{\rm BP→TP}(t)$ ersetzen.

Leistung und Energie von QAM–Signalen

Wie im Kapitel 4.3 von „Signaldarstellung” gezeigt wird, kann die Leistung des QAM–Sendesignals $s(t)$ auch aus dem äquivalenten TP–Signal $s_{\rm TP}(t)$ berechnet werden, das stets komplex ist: $$P = \lim_{T_{\rm M} \rightarrow \infty} \frac{\rm 1}{T_{\rm M}}\cdot \int_{-T_{\rm M}/2}^{+T_{\rm M}/2} s(t)^2\,{\rm d} t = \frac{\rm 1}{2} \cdot \lim_{T_{\rm M} \rightarrow \infty} \frac{\rm 1}{T_{\rm M}}\cdot \int_{-T_{\rm M}/2}^{+T_{\rm M}/2} |s_{\rm TP}(t)|^2\,{\rm d} t \hspace{0.05cm}.$$ Dagegen ist die Energie der Signale $s(t)$ und $s_{\rm TP}(t)$ unendlich groß. Beschränkt man sich jedoch auf eine Symboldauer $T$, so erhält man die Energie pro Symbol $$E_{\rm S} = \frac{{\rm E}[\hspace{0.05cm}|a_{\nu} |^2\hspace{0.05cm}]}{2}\cdot \int_{-\infty}^{+\infty} |g_s(t)|^2\,{\rm d} t = \frac{{\rm E}[|a_{\nu} |^2]}{2}\cdot \int_{-\infty}^{+\infty} |G_s(f)|^2\,{\rm d} f \hspace{0.05cm}.$$ Dagegen gibt $E_{\rm B} = E_{\rm S}/b$ die Energie pro Bit an, wenn gemäß der gegebenen Signalraumzuordnung jeweils $b$ Binärsymbole zu einem komplexen Amplitudenkoeffizienten zusammengefasst werden.

Die obere Grafik zeigt die Signalraumzuordnung bei 16–QAM, wobei sowohl der Real– als auch der Imaginärteil der komplexen Amplitudenkoeffizienten $a_ν$ jeweils einen von vier Werten (±1 sowie ±1/3) annehmen kann. Durch Mittelung über die 16 Abstandsquadrate zum Ursprung erhält man: $$\begin{align*}{\rm E}[\hspace{0.05cm}|a_{\nu} |^2 \hspace{0.05cm}] \hspace{-0.18cm} & = \hspace{-0.18cm} \frac{4}{16} \cdot (1^2 + 1^2)+ \frac{4}{16} \cdot \left[({1}/{3} )^2 +({1}/{3})^2 \right ]+\\ \hspace{-0.18cm} & + \hspace{-0.18cm} \frac{8}{16} \cdot \left [1^2 + ({1}/{3})^2\right ] = ... \hspace{0.15cm}= \frac{10}{9}\approx 1.11 \hspace{0.05cm}.\end{align*}$$

Die Summanden gehören in dieser Reihenfolge zu den vier roten, den vier schwarzen und den acht blauen Punkten.

Bei NRZ–rechteckförmigem Sendegrundimpuls $g_s(t)$ mit der Amplitude $g_0$ und der Symboldauer $T$ ist das Spektrum $G_s(f)$ si–förmig. In diesem Fall gilt für

die mittlere Energie pro Symbol:

$$E_{\rm S} = \frac{{\rm E}[\hspace{0.05cm}|a_{\nu} |^2]\hspace{0.05cm}}{2}\cdot g_0^2 \cdot T = \frac{5}{9}\cdot g_0^2 \cdot T \hspace{0.05cm},$$

die mittlere Energie pro Bit:

$$E_{\rm B} = \frac{E_{\rm S}}{4}\approx 0.139 \cdot g_0^2 \cdot T \hspace{0.05cm}.$$

Die maximale Hüllkurve $s_0$ ist um den Faktor „Wurzel aus 2” größer als die Amplitude $g_0$ des Rechteckimpulses (siehe untere Skizze) und tritt bei einem der roten Amplitudenkoeffizienten auf, also immer dann, wenn $|a_{\rm Iν}| = |a_{\rm Qν}| =$ 1 ist.

Signalverläufe der 4–QAM

Die folgende Grafik zeigt die Signalverläufe der 4–QAM, wobei die Farbgebung mit der eingezeichneten Signalraumzuordnung übereinstimmt.

Man erkennt aus diesen Darstellungen:

die Seriell–Parallel–Wandlung des Quellensignals $q(t)$ in die beiden Komponentensignale $s_{\rm I}(t)$ und $s_{\rm Q}(t)$, jeweils mit der Symboldauer $T = 2T_{\rm B}$ und den Signalwerten $±g_0 (T_{\rm B}$ ist die Bitdauer);
die beiden trägerfrequenzmodulierten Signale $s_{\rm cos}(t)$ und $s_{\rm –sin}(t)$ mit Phasensprüngen um $±π$:

$$s_{\rm cos} (t) = s_{\rm I} (t) \cdot \cos(2\pi f_{\rm T}t)\hspace{0.05cm},\hspace{0.2cm} s_{\rm -sin} (t) = -s_{\rm Q} (t) \cdot \sin(2\pi f_{\rm T}t)\hspace{0.05cm}, $$

das Sendesignal $s(t) = s_{\rm cos}(t) \ – \ s_{\rm –sin}(t)$ mit Phasensprüngen um Vielfache von $±π/2$; deren Hüllkurve ist gegenüber den beiden Komponentensignalen um den Faktor „Wurzel aus 2” größer:

$$s_0 = \sqrt{2} \cdot g_0 \hspace{0.05cm}.$$

Anzumerken ist, dass hier der Sendegrundimpuls $g_s(t)$ zur Vereinfachung der Darstellung im Bereich von 0 bis $T$ als rechteckförmig (also unsymmetrisch bezüglich $t =$ 0) angenommen wurde. Die zugehörige Spektralfunktion $G_s(f)$ dieses kausalen Impulses ist komplex, was jedoch in diesem Zusammenhang keine Auswirkungen hat.

Fehlerwahrscheinlichkeit der 4–QAM

Im Kapitel 4.2 wurde die Fehlerwahrscheinlichkeit der BPSK angegeben. Nun werden die Ergebnisse auf die 4–QAM übertragen, wobei weiterhin folgende Voraussetzungen gelten:

ein Sendesignal mit der mittleren Energie $E_{\rm B}$ pro Bit,
AWGN–Rauschen mit der Rauschleistungsdichte $N_0$,
bestmögliche Empfängerrealisierung nach dem Matched–Filter–Prinzip.

Die obere Abbildung zeigt das BPSK–Phasendiagramm des Detektionssignals $d(t)$, also inklusive dem Matched–Filter. Der Abstand des Nutzsignals von der Schwelle $(d_{\rm Q}$–Achse) beträgt zu den Detektionszeitpunkten jeweils $s_0$. Mit den weiteren Gleichungen $$p_{\rm B} = {\rm Q}\left ( \frac{s_0}{\sigma_d } \right ), \hspace{0.2cm} E_{\rm B} = \frac {1}{2} \cdot s_0^2 \cdot T_{\rm B} ,\hspace{0.2cm} \sigma_d^2 = \frac{N_0}{T_{\rm B} }$$ erhält man für die BPSK–Fehlerwahrscheinlichkeit: $$p_{\rm B, \hspace{0.05cm}BPSK} = {\rm Q}\left ( \sqrt{{2 \cdot E_{\rm B}}/{N_0 }} \hspace{0.1cm}\right ) = \frac{1}{2}\cdot {\rm erfc}\left ( \sqrt{{E_{\rm B}}/{N_0 }} \hspace{0.1cm}\right ).$$

Bei der 4–QAM entsprechend der unteren Abbildung

gibt es nun zwei Schwellen zwischen den Bereichen mit hellerem/dunklerem Hintergrund (blaue Linie) sowie zwischen den gepunkteten/gestrichelten Flächen (rot),
ist der Abstand von den Schwellen jeweils nur noch $g_0$ anstelle von $s_0$,
ist aber die Rauschleistung wegen der halb so großen Symbolrate in jedem Teilzweig gegenüber der BPSK auch nur halb so groß.

Mit den Gleichungen $$p_{\rm B} = {\rm Q}\left ( \frac{g_0}{\sigma_d } \right ), \hspace{0.2cm}g_{0} = \frac {s_0}{\sqrt{2}}, \hspace{0.2cm}E_{\rm B} = \frac {1}{2} \cdot s_0^2 \cdot T_{\rm B} ,\hspace{0.2cm} \sigma_d^2 = \frac{N_0}{2 \cdot T_{\rm B} }$$ erhält man für die 4–QAM–Fehlerwahrscheinlichkeit genau das gleiche Ergebnis wie für die BPSK: $$p_{\rm B, \hspace{0.05cm}4-QAM} = {\rm Q}\left ( \sqrt{{2 \cdot E_{\rm B}}/{N_0 }} \hspace{0.1cm}\right ) = \frac{1}{2}\cdot {\rm erfc}\left ( \sqrt{{E_{\rm B}}/{N_0 }} \hspace{0.1cm}\right ).$$

Fazit: Die 4–QAM weist bei idealen Bedingungen die gleiche Fehlerwahrscheinlichkeit wie die BPSK auf, obwohl die doppelte Informationsmenge übertragen werden kann. Sind allerdings die Bedingungen nicht mehr ideal – zum Beispiel bei einem ungewollten Phasenversatz zwischen Sender und Empfänger – so gibt es bei der 4–QAM eine deutlich stärkere Degradation als bei der BPSK. Dieser Fall wird im Kapitel 1.5 des Buches „Digitalsignalübertragung” noch genauer betrachtet.

QAM–Signalraumkonstellationen

Die nachfolgende Grafik zeigt die Signalraumkonstellationen von 4–QAM, 16–QAM und 64–QAM. Mit den hier gewählten Achsenbeschriftungen beschreiben die Bilder auch das Detektionsnutzsignal (zu den Detektionszeitpunkten) im äquivalenten Tiefpassbereich. Ebenfalls eingezeichnet sind die verschiedenen Entscheidungsgebiete, die dem verrauschten Detektionssignal zugeordnet werden. Die Pfeile geben an, wenn Entscheidungsgebiete bis ins Unendliche ausgedehnt sind.

Zu diesen Darstellungen ist anzumerken:

Die Bilder beziehen sich nur auf die Detektionszeitpunkte und gelten für alle Nyquistsysteme wie die Rechteck–Rechteck–Konfiguration oder ein Wurzel–Cosinus–Rolloff–Nyquistsystem.
Die hier nicht dargestellten Übergänge außerhalb der Detektionszeitpunkte zwischen den einzelnen Punkten hängen dagegen sehr wohl vom gewählten Nyquistsystem ab.
Bei echter QAM–Struktur – das heißt: die Signalraumkonstellation ist quadratisch oder zumindest rechteckförmig – lässt sich die 2D–Detektion durch zwei „eindimensionale” Detektionsvorgänge vereinfacht lösen.
Die 16–QAM ist somit nichts anderes als die parallele Übertragung zweier Digitalsignale mit jeweils $M =$ 4 Amplitudenstufen. Bei der 64–QAM würde entsprechend $M =$ 8 gelten und bei der 256–QAM ist die „eindimensionale” Stufenzahl $M =$ 16.
Alle im Kapitel 2.2 von Buch „Digitalsignalübertragung” genannten Eigenschaften für mehrstufige Signale gelten auch hier, wobei allerdings die Zusetzung der orthogoalen Trägerfrequenzsignale noch geeignet zu berücksichtigen ist.

Weitere Signalraumkonstellationen

Die Grafik zeigt weitere Signalraumkonstellationen, wobei links wieder die 4–QAM nach der bisherigen Beschreibung dargestellt ist.

Die Konstellation rechts daneben kennzeichnet eine vierstufige Phasenmodulation (4–PSK bzw. QPSK – Quaternary Phase Shift Keying). Ein Vergleich der beiden linken Diagramme zeigt:

Die hier als QPSK bezeichnete Variante verwendet die Phasenlagen 0°, 90°, 180° und 270°. Man erkennt aus den eingezeichneten Entscheidungsgebieten, dass hier die Detektion nicht auf zwei Binärentscheidungen zurückgeführt werden kann.
Auch die 4–QAM kann als eine vierstufige Phasenmodulation mit den möglichen Phasenlagen ±45° und ±135° aufgefasst werden. Gegenüber dem zweiten Phasendiagramm ergibt sich eine Drehung um 45° $(π$/4), so dass die 4–QAM oft auch als $π$/4–QPSK bezeichnet wird.
Ähnlich wie man bei der BPSK durch Vorcodierung zur DPSK kommt, kann auch die 4–PSK zur 4–DPSK erweitert und dadurch die Demodulation erleichtert werden. Diese wird zum Beispiel bei der Datenübertragung über Telefonkanäle gemäß der CCITT–Empfehlung V26 angewendet (Trägerfrequenz 1800 Hz, Datenrate 2400 bit/s).

Die beiden rechten Diagramme zeigen höherstufige Modulationsverfahren:

Die 8–PSK bzw. 8–DPSK erlaubt entsprechend der CCITT–Empfehlung V27 beim Telefonkanal eine Datenrate von bis zu 4800 bit/s.
Die Empfehlung V29 sieht mit der 16–ASK/PSK eine hybride Modulationsform vor, die bei ausgewählten, fest verschalteten Leitungen eine Datenrate von 9600 bit/s ermöglicht.

Nyquist– und Wurzel–Nyquist–QAM–Systeme

Bisher wurde im Kapitel 4.3 aus Darstellungsgründen stets vom rechteckförmigen Sendegrundimpuls ausgegangen. In der Praxis verwendet man aber meist eine Wurzel–Nyquist–Charakteristik entsprechend dem Kapitel 1.4 im Buch „Digitalsignalübertragung”. In aller Kürze lassen sich diese Systeme wie folgt darstellen:

Der Empfängerfrequenzgang $H_{\rm E}(f)$ wird hier formgleich mit dem normierten Sendeimpulsspektrum $H_{\rm S}(f)$ gewählt, was unter der Nebenbedingung der Leistungsbegrenzung (das heißt: konstante mittlere Sendeleistung) zur kleinsten Fehlerwahrscheinlichkeit führt.
Der Gesamtfrequenzgang $H_{\rm Nyq}(f) = H_{\rm S}(f) · H_{\rm E}(f)$ erfüllt das erste Nyquistkriterium, so dass es beim Empfänger zu keinen Impulsinterferenzen kommt. $H_{\rm S}(f)$ und $H_{\rm E}(f)$ haben somit jeweils eine Wurzel–Nyquist–Charakteristik.
Für den Frequenzgang $H_{\rm Nyq}(f)$ verwendet man einen Cosinus–Rolloff–Tiefpass $H_{\rm CRO}(f)$ mit der äquivalenten Bandbreite $Δf_{\rm CRO} = 1/T$ und frei wählbarem Rolloff–Faktor $(0 ≤ r ≤ 1)$.

Mit dem nachfolgenden Interaktionsmodul können Sie sich den Frequenzgang und die Impulsantwort dieses Tiefpasses verdeutlichen:

Frequenzgang und zugehörige Impulsantwort

Der Vorteil dieser Wurzel–Nyquist–Systeme ist die deutlich kleinere Bandbreite $(1 + r)/T$ gegenüber der bisher betrachteten Konfiguration mit rechteckförmigem $g_s(t)$ und rechteckförmigem $h_{\rm E}(t)$, dessen Spektrum (theoretisch) unendlich weit ausgedehnt ist. Hinsichtlich Fehlerwahrscheinlichkeit ändert sich gegenüber der Rechteck–Rechteck–Konfiguration nichts, da der Grundimpuls $g_d(t)$ vor dem Entscheider äquidistante Nulldurchgänge aufweist und somit Impulsinterferenzen vermieden werden.

Die Grafik zeigt die Phasendiagramme für diesen Fall, die dem Buch [Kam04]^[1] entnommen sind. Der Rolloff–Faktor beträgt $r =$ 0.5.