Digitalsignalübertragung/Struktur des optimalen Empfängers: Unterschied zwischen den Versionen
(12 dazwischenliegende Versionen desselben Benutzers werden nicht angezeigt) | |||
Zeile 8: | Zeile 8: | ||
== Blockschaltbild und Voraussetzungen == | == Blockschaltbild und Voraussetzungen == | ||
<br> | <br> | ||
− | In diesem Kapitel wird die Struktur des optimalen Empfängers eines digitalen Übertragungssystems sehr allgemein hergeleitet, wobei | + | In diesem Kapitel wird die Struktur des optimalen Empfängers eines digitalen Übertragungssystems sehr allgemein hergeleitet, wobei |
*das Modulationsverfahren und weitere Systemdetails nicht weiter spezifiziert werden,<br> | *das Modulationsverfahren und weitere Systemdetails nicht weiter spezifiziert werden,<br> | ||
− | |||
+ | *von den Basisfunktionen und der Signalraumdarstellung gemäß dem Kapitel [[Digitalsignal%C3%BCbertragung/Signale,_Basisfunktionen_und_Vektorr%C3%A4ume#Zur_Nomenklatur_im_vierten_Kapitel_.281.29|"Signale, Basisfunktionen und Vektorräume"]] ausgegangen wird. | ||
+ | [[Datei:P ID2000 Dig T 4 2 S1 version1.png|right|frame|Allgemeines Blockschaltbild eines Kommunikationssystems|class=fit]] | ||
− | |||
Zum obigen Blockschaltbild ist anzumerken: | Zum obigen Blockschaltbild ist anzumerken: | ||
− | *Der Symbolumfang der Quelle beträgt $M$ und der Symbolvorrat ist $\{m_i\}$ mit $i = 0$, ... , $M- 1$. | + | *Der Symbolumfang der Quelle beträgt $M$ und der Symbolvorrat ist $\{m_i\}$ mit $i = 0$, ... , $M-1$. |
− | * | + | *Die zugehörigen Symbolwahrscheinlichkeiten ${\rm Pr}(m = m_i)$ seien auch dem Empfänger bekannt.<br> |
− | *Lineare Kanalverzerrungen sind in der obigen Grafik durch die Impulsantwort $h(t)$ berücksichtigt. Außerdem ist ein (irgendwie geartetes) Rauschen $n(t)$ wirksam. Mit diesen beiden die Übertragung störenden Effekten lässt sich das am Empfänger ankommende Signal $r(t)$ in folgender Weise angeben: | + | *Zur Nachrichtenübertragung stehen $M$ Signalformen $s_i(t)$ zur Verfügung; auch für die Laufvariable gelte die Indizierung $i = 0$, ... , $M-1$. |
+ | |||
+ | *Es besteht eine feste Beziehung zwischen den Nachrichten und den Signalen. Wird $m =m_i$ übertragen, so ist das Sendesignal $s(t) =s_i(t)$.<br> | ||
+ | |||
+ | *Lineare Kanalverzerrungen sind in der obigen Grafik durch die Impulsantwort $h(t)$ berücksichtigt. Außerdem ist ein (irgendwie geartetes) Rauschen $n(t)$ wirksam. | ||
+ | |||
+ | *Mit diesen beiden die Übertragung störenden Effekten lässt sich das am Empfänger ankommende Signal $r(t)$ in folgender Weise angeben: | ||
:$$r(t) = s(t) \star h(t) + n(t) \hspace{0.05cm}.$$ | :$$r(t) = s(t) \star h(t) + n(t) \hspace{0.05cm}.$$ | ||
− | *Aufgabe des (optimalen) Empfängers ist es, anhand seines Eingangssignals $r(t)$ herauszufinden, welche der $M$ möglichen Nachrichten $m_i$ & | + | *Aufgabe des (optimalen) Empfängers ist es, anhand seines Eingangssignals $r(t)$ herauszufinden, welche der $M$ möglichen Nachrichten $m_i$ ⇒ Signal $s_i(t)$ gesendet wurde. Der vom Empfänger gefundene Schätzwert für $m$ wird durch ein "Zirkumflex" gekennzeichnet ⇒ $\hat{m}$. |
{{BlaueBox|TEXT= | {{BlaueBox|TEXT= | ||
− | $\text{Definition:}$ Man spricht von einem '''optimalen Empfänger''', wenn die Symbolfehlerwahrscheinlichkeit den für die Randbedingungen | + | $\text{Definition:}$ Man spricht von einem '''optimalen Empfänger''', wenn die Symbolfehlerwahrscheinlichkeit den für die Randbedingungen kleinstmöglichen Wert annimmt: |
:$$p_{\rm S} = {\rm Pr} ({\cal E}) = {\rm Pr} ( \hat{m} \ne m) \hspace{0.15cm} \Rightarrow \hspace{0.15cm}{\rm Minimum} \hspace{0.05cm}.$$}} | :$$p_{\rm S} = {\rm Pr} ({\cal E}) = {\rm Pr} ( \hat{m} \ne m) \hspace{0.15cm} \Rightarrow \hspace{0.15cm}{\rm Minimum} \hspace{0.05cm}.$$}} | ||
− | + | Hinweise: | |
− | #Im Folgenden wird meist der AWGN–Ansatz ⇒ $r(t) = s(t) + n(t)$ vorausgesetzt, was bedeutet, dass $h(t) = \delta(t)$ als verzerrungsfrei angenommen wird. | + | #Im Folgenden wird meist der AWGN–Ansatz ⇒ $r(t) = s(t) + n(t)$ vorausgesetzt, was bedeutet, dass $h(t) = \delta(t)$ als verzerrungsfrei angenommen wird. |
− | #Andernfalls | + | #Andernfalls können wir die Signale $s_i(t)$ als ${s_i}'(t) = s_i(t) \star h(t)$ neu definieren, also die deterministischen Kanalverzerrungen dem Sendesignal beaufschlagen.<br> |
== Fundamentaler Ansatz zum optimalen Empfängerentwurf== | == Fundamentaler Ansatz zum optimalen Empfängerentwurf== | ||
<br> | <br> | ||
− | Gegenüber dem auf der vorherigen Seite gezeigten [[Digitalsignal%C3%BCbertragung/Struktur_des_optimalen_Empf%C3%A4ngers#Blockschaltbild_und_Voraussetzungen| Blockschaltbild]] führen wir nun einige wesentliche Verallgemeinerungen durch: | + | Gegenüber dem auf der vorherigen Seite gezeigten [[Digitalsignal%C3%BCbertragung/Struktur_des_optimalen_Empf%C3%A4ngers#Blockschaltbild_und_Voraussetzungen| "Blockschaltbild"]] führen wir nun einige wesentliche Verallgemeinerungen durch: |
− | + | [[Datei:P ID2001 Dig T 4 2 S2 version2.png|right|frame|Modell zur Herleitung des optimalen Empfängers|class=fit]] | |
− | * | + | *Der Übertragungskanal wird nun beschrieben durch die [[Stochastische_Signaltheorie/Statistische_Abhängigkeit_und_Unabhängigkeit#Bedingte_Wahrscheinlichkeit|"bedingte Wahrscheinlichkeitsdichtefunktion"]] $p_{\hspace{0.02cm}r(t)\hspace{0.02cm} \vert \hspace{0.02cm}s(t)}$, welche die Abhängigkeit des Empfangssignals $r(t)$ vom Sendesignal $s(t)$ festlegt.<br> |
− | |||
− | :unter Berücksichtigung aller möglichen Sendesignale $s_i(t)$ und | + | *Wurde nun ein ganz bestimmtes Signal $r(t) = \rho(t)$ empfangen, so hat der Empfänger die Aufgabe, die Wahrscheinlichkeitsdichtefunktionen auf der Grundlage dieser "Signalrealisierung" $\rho(t)$ und der $M$ bedingten Wahrscheinlichkeitsdichtefunktionen zu bestimmen: |
+ | :$$p_{\hspace{0.05cm}r(t) \hspace{0.05cm} \vert \hspace{0.05cm} s(t) } (\rho(t) \hspace{0.05cm} \vert \hspace{0.05cm} s_i(t))\hspace{0.5cm}{\rm mit}\hspace{0.5cm} i = 0, \text{...} \hspace{0.05cm}, M-1.$$ | ||
+ | *Ermittelt werden soll, welche Nachricht $\hat{m}$ am wahrscheinlichsten gesendet wurde, unter Berücksichtigung aller möglichen Sendesignale $s_i(t)$ und derenn Auftrittsswahrscheinlichkeiten ${\rm Pr}(m = m_i)$. | ||
*Die Schätzung des optimalen Empfängers ist also ganz allgemein bestimmt durch die Gleichung | *Die Schätzung des optimalen Empfängers ist also ganz allgemein bestimmt durch die Gleichung | ||
− | :$$\hat{m} = {\rm arg} \max_i \hspace{0.1cm} p_{\hspace{0.02cm}s(t) \hspace{0.05cm} \vert \hspace{0.05cm} r(t) } ( s_i(t) \hspace{0.05cm} \vert \hspace{0.05cm} \rho(t)) = {\rm arg} \max_i \hspace{0.1cm} p_{m \hspace{0.05cm} \vert \hspace{0.05cm} r(t) } ( \hspace{0.05cm}m_i\hspace{0.05cm} \vert \hspace{0.05cm}\rho(t))\hspace{0.05cm} | + | :$$\hat{m} = {\rm arg} \max_i \hspace{0.1cm} p_{\hspace{0.02cm}s(t) \hspace{0.05cm} \vert \hspace{0.05cm} r(t) } ( s_i(t) \hspace{0.05cm} \vert \hspace{0.05cm} \rho(t)) = {\rm arg} \max_i \hspace{0.1cm} p_{m \hspace{0.05cm} \vert \hspace{0.05cm} r(t) } ( \hspace{0.05cm}m_i\hspace{0.05cm} \vert \hspace{0.05cm}\rho(t))\hspace{0.05cm}.$$ |
− | : | + | {{BlaueBox|TEXT= |
+ | $\text{In anderen Worten:}$ Der optimale Empfänger betrachtet diejenige Nachricht $\hat{m} \in \{m_i\}$ als die am wahrscheinlichsten gesendete Nachricht, deren bedingte Wahrscheinlichkeitsdichtefunktion $p_{\hspace{0.02cm}m \hspace{0.05cm} \vert \hspace{0.05cm} r(t) }$ für das anliegende Empfangssignal $\rho(t)$ sowie unter der Annahme $m=\hat{m}$ den größtmöglichen Wert annimmt.}} | ||
− | |||
− | |||
− | Bevor wir die | + | Bevor wir die diese Entscheidungsregel näher diskutieren, soll der optimale Empfänger entsprechend der Grafik noch in zwei Funktionsblöcke aufgeteilt werden: |
− | + | *Der '''Detektor''' nimmt am Empfangssignal $r(t)$ verschiedene Messungen vor und fasst diese im Vektor $\boldsymbol{r}$ zusammen. Bei $K$ Messungen entspricht $\boldsymbol{r}$ einem Punkt im $K$–dimensionalen Vektorraum.<br> | |
− | *Der Detektor nimmt am Empfangssignal $r(t)$ verschiedene Messungen vor und fasst diese im Vektor $\boldsymbol{r}$ zusammen. Bei $K$ Messungen entspricht $\boldsymbol{r}$ einem Punkt im $K$–dimensionalen Vektorraum.<br> | ||
− | *Der Entscheider bildet abhängig von diesem Vektor den Schätzwert. Bei einem gegebenen Vektor $\boldsymbol{r} = \boldsymbol{\rho}$ lautet dabei die Entscheidungsregel: | + | *Der '''Entscheider''' bildet abhängig von diesem Vektor den Schätzwert. Bei einem gegebenen Vektor $\boldsymbol{r} = \boldsymbol{\rho}$ lautet dabei die Entscheidungsregel: |
:$$\hat{m} = {\rm arg}\hspace{0.05cm} \max_i \hspace{0.1cm} P_{m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}) \hspace{0.05cm}.$$ | :$$\hat{m} = {\rm arg}\hspace{0.05cm} \max_i \hspace{0.1cm} P_{m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}) \hspace{0.05cm}.$$ | ||
− | + | Im Gegensatz zur oberen Entscheidungsregel tritt nun eine bedingte Wahrscheinlichkeit $P_{m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} }$ anstelle der bedingten Wahrscheinlichkeitskeitsdichtefunktion $\rm (WDF)$ $p_{m\hspace{0.05cm} \vert \hspace{0.05cm}r(t)}$ auf. Beachten Sie bitte die Groß– bzw. Kleinschreibung für die unterschiedlichen Bedeutungen. | |
− | Im Gegensatz zur oberen Entscheidungsregel tritt nun eine bedingte Wahrscheinlichkeit $P_{m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} }$ anstelle der bedingten Wahrscheinlichkeitskeitsdichtefunktion (WDF) $p_{m\hspace{0.05cm} \vert \hspace{0.05cm}r(t)}$ auf. Beachten Sie bitte die Groß– bzw. Kleinschreibung für die unterschiedlichen Bedeutungen. | ||
<br clear=all> | <br clear=all> | ||
{{GraueBox|TEXT= | {{GraueBox|TEXT= | ||
− | $\text{Beispiel 1:}$ Wir betrachten nun die Funktion $y = {\rm arg}\hspace{0.05cm} \max \ p(x)$, wobei $p(x)$ die Wahrscheinlichkeitsdichtefunktion (WDF) einer wertkontinuierlichen oder wertdiskreten Zufallsgröße $x$ beschreibt. | + | $\text{Beispiel 1:}$ Wir betrachten nun die Funktion $y = {\rm arg}\hspace{0.05cm} \max \ p(x)$, wobei $p(x)$ die Wahrscheinlichkeitsdichtefunktion $\rm (WDF)$ einer wertkontinuierlichen oder wertdiskreten Zufallsgröße $x$ beschreibt. In der rechten Grafik besteht die WDF aus einer Summe von Diracfunktionen mit Wahrscheinlichkeiten als Impulsgewichte.<br> |
[[Datei:P ID2002 Dig T 4 2 S2b version1.png|righ|frame|Zur Verdeutlichung der Funktion „arg max”|class=fit]] | [[Datei:P ID2002 Dig T 4 2 S2b version1.png|righ|frame|Zur Verdeutlichung der Funktion „arg max”|class=fit]] | ||
− | Die Grafik zeigt beispielhafte Funktionen. In beiden Fällen liegt das WDF–Maximum $(17)$ bei $x = 6$: | + | ⇒ Die Grafik zeigt beispielhafte Funktionen. In beiden Fällen liegt das WDF–Maximum $(17)$ bei $x = 6$: |
:$$\max_i \hspace{0.1cm} p(x) = 17\hspace{0.05cm},$$ | :$$\max_i \hspace{0.1cm} p(x) = 17\hspace{0.05cm},$$ | ||
:$$y = {\rm \hspace{0.05cm}arg} \max_i \hspace{0.1cm} p(x) = 6\hspace{0.05cm}.$$ | :$$y = {\rm \hspace{0.05cm}arg} \max_i \hspace{0.1cm} p(x) = 6\hspace{0.05cm}.$$ | ||
− | Die (bedingten) Wahrscheinlichkeiten in der Gleichung | + | ⇒ Die (bedingten) Wahrscheinlichkeiten in der Gleichung |
:$$\hat{m} = {\rm arg}\hspace{0.05cm} \max_i \hspace{0.1cm} P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho})$$ | :$$\hat{m} = {\rm arg}\hspace{0.05cm} \max_i \hspace{0.1cm} P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho})$$ | ||
− | sind '''a–Posteriori–Wahrscheinlichkeiten'''. | + | sind '''a–Posteriori–Wahrscheinlichkeiten'''. Mit dem [[Stochastische_Signaltheorie/Statistische_Abh%C3%A4ngigkeit_und_Unabh%C3%A4ngigkeit#R.C3.BCckschlusswahrscheinlichkeit| Satz von Bayes]] kann hierfür geschrieben werden: |
− | |||
− | Mit dem [[Stochastische_Signaltheorie/Statistische_Abh%C3%A4ngigkeit_und_Unabh%C3%A4ngigkeit#R.C3.BCckschlusswahrscheinlichkeit| Satz von Bayes]] kann hierfür geschrieben werden: | ||
:$$P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}) = | :$$P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}) = | ||
\frac{ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )}{p_{\boldsymbol{ r} } (\boldsymbol{\rho})} | \frac{ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )}{p_{\boldsymbol{ r} } (\boldsymbol{\rho})} | ||
Zeile 84: | Zeile 86: | ||
− | Der Nennerterm ist für alle Alternativen $m_i$ gleich und muss für die Entscheidung nicht berücksichtigt werden. Damit erhält man folgende Regeln: | + | Der Nennerterm $p_{\boldsymbol{ r} }(\boldsymbol{\rho})$ ist für alle Alternativen $m_i$ gleich und muss für die Entscheidung nicht berücksichtigt werden. Damit erhält man folgende Regeln: |
{{BlaueBox|TEXT= | {{BlaueBox|TEXT= | ||
− | $\text{ | + | $\text{Theorem:}$ Die Entscheidungsregel des optimalen Empfängers, auch bekannt als '''Maximum–a–posteriori–Empfänger''' $($kurz: '''MAP–Empfänger'''$)$ lautet: |
:$$\hat{m}_{\rm MAP} = {\rm \hspace{0.05cm} arg} \max_i \hspace{0.1cm} P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}) = {\rm \hspace{0.05cm}arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm} m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm} m_i )\big ]\hspace{0.05cm}.$$ | :$$\hat{m}_{\rm MAP} = {\rm \hspace{0.05cm} arg} \max_i \hspace{0.1cm} P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}) = {\rm \hspace{0.05cm}arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm} m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm} m_i )\big ]\hspace{0.05cm}.$$ | ||
− | + | *Vorteil der zweiten Gleichung ist, dass die Vorwärtsrichtung des Kanals Anwendung findet ⇒ bedingte WDF $p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm} m }$ $($"Ausgang unter der Bedingung Eingang"$)$. | |
+ | |||
+ | *Dagegen verwendet die erste Gleichung die Rückschlusswahrscheinlichkeiten $P_{\hspace{0.05cm}m\hspace{0.05cm} \vert \hspace{0.02cm} \boldsymbol{ r} } $ $($"Eingang unter der Bedingung Ausgang"$)$.}} | ||
{{BlaueBox|TEXT= | {{BlaueBox|TEXT= | ||
− | $\text{ | + | $\text{Theorem:}$ Der '''Maximum–Likelihood–Empfänger''' $($kurz: '''ML–Empfänger'''$)$ verwendet die Entscheidungsregel |
:$$\hat{m}_{\rm ML} = \hspace{-0.1cm} {\rm arg} \max_i \hspace{0.1cm} p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )\hspace{0.05cm}.$$ | :$$\hat{m}_{\rm ML} = \hspace{-0.1cm} {\rm arg} \max_i \hspace{0.1cm} p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )\hspace{0.05cm}.$$ | ||
− | Bei diesem werden die möglicherweise unterschiedlichen Auftrittswahrscheinlichkeiten ${\rm Pr}(m = m_i)$ für den Entscheidungsprozess nicht herangezogen | + | *Bei diesem werden die möglicherweise unterschiedlichen Auftrittswahrscheinlichkeiten ${\rm Pr}(m = m_i)$ für den Entscheidungsprozess nicht herangezogen. |
+ | * Zum Beispiel, weil diese dem Empfänger nicht bekannt sind.}}<br> | ||
− | Im früheren Kapitel [[Digitalsignalübertragung/Optimale_Empfängerstrategien | + | Im früheren Kapitel [[Digitalsignalübertragung/Optimale_Empfängerstrategien|"Optimale Empfängerstrategien"]] finden Sie auch andere Herleitungen für diese Empfängertypen. |
{{BlaueBox|TEXT= | {{BlaueBox|TEXT= | ||
− | $\text{Fazit:}$ Bei gleichwahrscheinlichen Nachrichten $\{m_i\}$ ⇒ ${\rm Pr}(m = m_i) = 1/M$ ist der im Allgemeinen etwas schlechtere | + | $\text{Fazit:}$ Bei gleichwahrscheinlichen Nachrichten $\{m_i\}$ ⇒ ${\rm Pr}(m = m_i) = 1/M$ ist der im Allgemeinen etwas schlechtere "Maximum–Likelihood–Empfänger" gleichwertig mit dem "Maximum–a–posteriori–Empfänger": |
:$$\hat{m}_{\rm MAP} = \hat{m}_{\rm ML} =\hspace{-0.1cm} {\rm\hspace{0.05cm} arg} \max_i \hspace{0.1cm} | :$$\hat{m}_{\rm MAP} = \hat{m}_{\rm ML} =\hspace{-0.1cm} {\rm\hspace{0.05cm} arg} \max_i \hspace{0.1cm} | ||
p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )\hspace{0.05cm}.$$}} | p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )\hspace{0.05cm}.$$}} | ||
Zeile 111: | Zeile 116: | ||
== Das Theorem der Irrelevanz== | == Das Theorem der Irrelevanz== | ||
<br> | <br> | ||
− | Zu beachten ist, dass der auf der letzten Seite beschriebene Empfänger nur dann optimal ist, wenn auch der Detektor bestmöglich implementiert ist, das heißt, wenn durch den Übergang vom kontinuierlichen Signal $r(t)$ zum Vektor $\boldsymbol{r}$ keine Information verloren geht.<br> | + | [[Datei:P ID2003 Dig T 4 2 S3a version2.png|right|frame|Zum Theorem der Irrelevanz|class=fit]] |
+ | Zu beachten ist, dass der auf der letzten Seite beschriebene Empfänger nur dann optimal ist, wenn auch der Detektor bestmöglich implementiert ist, das heißt, wenn durch den Übergang vom kontinuierlichen Signal $r(t)$ zum Vektor $\boldsymbol{r}$ keine Information verloren geht.<br> | ||
− | + | Um zu klären, welche und wieviele Messungen am Empfangssignal $r(t)$ nötig sind, um Optimalität zu garantieren, ist das "Theorem der Irrelevanz" hilfreich: | |
+ | *Dazu betrachten wir den skizzierten Empfänger, dessen Detektor aus dem Empfangssignal $r(t)$ die zwei Vektoren $\boldsymbol{r}_1$ und $\boldsymbol{r}_2$ ableitet und dem Entscheider zur Verfügung stellt. | ||
+ | |||
+ | *Diese Größen stehen mit der Nachricht $ m \in \{m_i\}$ über die Verbundwahrscheinlichkeitsdichte $p_{\boldsymbol{ r}_1, \hspace{0.05cm}\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm}m }$ in Zusammenhang.<br> | ||
− | + | *Die Entscheidungsregel des MAP–Empfängers lautet mit Anpassung an dieses Beispiel: | |
− | |||
− | Die Entscheidungsregel des MAP–Empfängers lautet mit Anpassung an dieses Beispiel: | ||
:$$\hat{m}_{\rm MAP} \hspace{-0.1cm} = \hspace{-0.1cm} {\rm arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}_1 , \hspace{0.05cm}\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm}m } \hspace{0.05cm} (\boldsymbol{\rho}_1, \hspace{0.05cm}\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} m_i ) \big]= | :$$\hat{m}_{\rm MAP} \hspace{-0.1cm} = \hspace{-0.1cm} {\rm arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}_1 , \hspace{0.05cm}\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm}m } \hspace{0.05cm} (\boldsymbol{\rho}_1, \hspace{0.05cm}\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} m_i ) \big]= | ||
Zeile 125: | Zeile 132: | ||
\hspace{0.05cm}.$$ | \hspace{0.05cm}.$$ | ||
− | + | *Die Vektoren $\boldsymbol{r}_1$ und $\boldsymbol{r}_2$ sind Zufallsgrößen. Ihre Realisierungen werden hier und im Folgenden mit $\boldsymbol{\rho}_1$ und $\boldsymbol{\rho}_2$ bezeichnet. Zur Hervorhebung sind alle Vektoren in der Grafik rot eingetragen. | |
− | *Die Vektoren $\boldsymbol{r}_1$ und $\boldsymbol{r}_2$ sind Zufallsgrößen. Ihre Realisierungen werden hier und im Folgenden mit $\boldsymbol{\rho}_1$ und $\boldsymbol{\rho}_2$ bezeichnet. Zur Hervorhebung sind alle Vektoren in der Grafik rot eingetragen. | ||
− | *Die Voraussetzungen für die Anwendung des & | + | *Die Voraussetzungen für die Anwendung des "Theorems der Irrelevanz" sind die gleichen wie die an eine [[Stochastische_Signaltheorie/Markovketten#Betrachtetes_Szenario| "Markovkette"]] erster Ordnung. Die Zufallsvariablen $x$, $y$, $z$ formen dann eine Markovkette erster Ordnung, falls die Verteilung von $z$ bei gegebenem $y$ unabhängig von $x$ ist: |
:$$p(x, y, z) = p(x) \cdot p(y\hspace{0.05cm} \vert \hspace{0.05cm}x) \cdot p(z\hspace{0.05cm} \vert \hspace{0.05cm}y) \hspace{0.25cm} {\rm anstelle \hspace{0.15cm}von} \hspace{0.25cm}p(x, y, z) = p(x) \cdot p(y\hspace{0.05cm} \vert \hspace{0.05cm}x) \cdot p(z\hspace{0.05cm} \vert \hspace{0.05cm}x, y) \hspace{0.05cm}.$$ | :$$p(x, y, z) = p(x) \cdot p(y\hspace{0.05cm} \vert \hspace{0.05cm}x) \cdot p(z\hspace{0.05cm} \vert \hspace{0.05cm}y) \hspace{0.25cm} {\rm anstelle \hspace{0.15cm}von} \hspace{0.25cm}p(x, y, z) = p(x) \cdot p(y\hspace{0.05cm} \vert \hspace{0.05cm}x) \cdot p(z\hspace{0.05cm} \vert \hspace{0.05cm}x, y) \hspace{0.05cm}.$$ | ||
− | *Der optimale Empfänger muss | + | *Der optimale Empfänger muss allgemein beide Vektoren $\boldsymbol{r}_1$ und $\boldsymbol{r}_2$ auswerten, da in obiger Entscheidungsregel beide Verbunddichten $p_{\boldsymbol{ r}_1\hspace{0.05cm} \vert \hspace{0.05cm}m }$ und $p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r}_1, \hspace{0.05cm}m }$ auftreten. Dagegen kann der Empfänger ohne Informationseinbuße die zweite Messung vernachlässigen, falls $\boldsymbol{r}_2$ bei gegebenem $\boldsymbol{r}_1$ unabhängig von der Nachricht $m$ ist: |
− | |||
− | |||
:$$p_{\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i )= | :$$p_{\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i )= | ||
p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 } \hspace{0.05cm} (\boldsymbol{\rho}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 ) | p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 } \hspace{0.05cm} (\boldsymbol{\rho}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 ) | ||
Zeile 152: | Zeile 156: | ||
\big]\hspace{0.05cm}.$$ | \big]\hspace{0.05cm}.$$ | ||
− | {{GraueBox|TEXT= | + | {{GraueBox|TEXT= |
− | $\text{Beispiel | + | [[Datei:Dig_T_4_2_S3b_version2.png|right|frame|Zwei Beispiele zum Theorem der Irrelevanz|class=fit]] |
+ | $\text{Beispiel 2:}$ Wir betrachten zur Verdeutlichung des soeben vorgestellten Theorems der Irrelevanz zwei verschiedene Systemkonfigurationen mit jeweils zwei Rauschtermen $\boldsymbol{ n}_1$ und $\boldsymbol{ n}_2$. | ||
+ | *In der Grafik sind alle vektoriellen Größen rot eingezeichnet. | ||
+ | |||
+ | *Die Größen $\boldsymbol{s}$, $\boldsymbol{ n}_1$ und $\boldsymbol{ n}_2$ seien zudem jeweils unabhängig voneinander.<br> | ||
− | |||
Die Analyse dieser beiden Anordnungen liefert folgende Ergebnisse: | Die Analyse dieser beiden Anordnungen liefert folgende Ergebnisse: | ||
− | *Der Entscheider muss in beiden Fällen die Komponente $\boldsymbol{ r}_1= \boldsymbol{ s}_i + \boldsymbol{ n}_1$ berücksichtigen, da nur diese die Information über das Nutzsignal $\boldsymbol{ s}_i$ und damit über die gesendete Nachricht $m_i$ liefert.<br> | + | *Der Entscheider muss in beiden Fällen die Komponente $\boldsymbol{ r}_1= \boldsymbol{ s}_i + \boldsymbol{ n}_1$ berücksichtigen, da nur diese die Information über das Nutzsignal $\boldsymbol{ s}_i$ und damit über die gesendete Nachricht $m_i$ liefert.<br> |
− | *Bei der oberen Konfiguration enthält $\boldsymbol{ r}_2$ keine Information über $m_i$, die nicht bereits von $\boldsymbol{ r}_1$ geliefert wurde. Vielmehr ist $\boldsymbol{ r}_2= \boldsymbol{ r}_1 + \boldsymbol{ n}_2$ nur eine verrauschte Version von $\boldsymbol{ r}_1$ und hängt nur vom Rauschen $\boldsymbol{ n}_2$ ab, sobald $\boldsymbol{ r}_1$ bekannt ist ⇒ $\boldsymbol{ r}_2$ ist irrelevant: | + | *Bei der oberen Konfiguration enthält $\boldsymbol{ r}_2$ keine Information über $m_i$, die nicht bereits von $\boldsymbol{ r}_1$ geliefert wurde. Vielmehr ist $\boldsymbol{ r}_2= \boldsymbol{ r}_1 + \boldsymbol{ n}_2$ nur eine verrauschte Version von $\boldsymbol{ r}_1$ und hängt nur vom Rauschen $\boldsymbol{ n}_2$ ab, sobald $\boldsymbol{ r}_1$ bekannt ist ⇒ $\boldsymbol{ r}_2$ '''ist irrelevant''': |
:$$p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i )= | :$$p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i )= | ||
p_{\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}_1 )= | p_{\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}_1 )= | ||
p_{\boldsymbol{ n}_2 } \hspace{0.05cm} (\boldsymbol{\rho}_2 - \boldsymbol{\rho}_1 )\hspace{0.05cm}.$$ | p_{\boldsymbol{ n}_2 } \hspace{0.05cm} (\boldsymbol{\rho}_2 - \boldsymbol{\rho}_1 )\hspace{0.05cm}.$$ | ||
− | *Bei der unteren Konfiguration ist | + | *Bei der unteren Konfiguration ist $\boldsymbol{ r}_2= \boldsymbol{ n}_1 + \boldsymbol{ n}_2$ für den Empfänger hilfreich, da ihm so ein Schätzwert für den Rauschterm $\boldsymbol{ n}_1$ geliefert wird ⇒ $\boldsymbol{ r}_2$ sollte hier nicht verworfen werden. Formal lässt sich dieses Resultat wie folgt ausdrücken: |
:$$p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i ) = p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ n}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 - \boldsymbol{s}_i, \hspace{0.05cm}m_i)= p_{\boldsymbol{ n}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ n}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2- \boldsymbol{\rho}_1 + \boldsymbol{s}_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 - \boldsymbol{s}_i, \hspace{0.05cm}m_i) = p_{\boldsymbol{ n}_2 } \hspace{0.05cm} (\boldsymbol{\rho}_2- \boldsymbol{\rho}_1 + \boldsymbol{s}_i ) | :$$p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i ) = p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ n}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 - \boldsymbol{s}_i, \hspace{0.05cm}m_i)= p_{\boldsymbol{ n}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ n}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2- \boldsymbol{\rho}_1 + \boldsymbol{s}_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 - \boldsymbol{s}_i, \hspace{0.05cm}m_i) = p_{\boldsymbol{ n}_2 } \hspace{0.05cm} (\boldsymbol{\rho}_2- \boldsymbol{\rho}_1 + \boldsymbol{s}_i ) | ||
\hspace{0.05cm}.$$ | \hspace{0.05cm}.$$ | ||
− | *Da nun im Argument dieser Funktion die Nachricht $\boldsymbol{ s}_i$ erscheint, ist $\boldsymbol{ r}_2$ & | + | *Da nun im Argument dieser Funktion die Nachricht $\boldsymbol{ s}_i$ erscheint, ist $\boldsymbol{ r}_2$ '''nicht irrelevant”, sondern durchaus relevant'''.}}<br> |
== Einige Eigenschaften des AWGN-Kanals== | == Einige Eigenschaften des AWGN-Kanals== | ||
<br> | <br> | ||
− | Um weitere Aussagen über die Art der optimalen Messungen des Vektors $\boldsymbol{ r}$ machen zu können, ist es notwendig, die den Kanal charakterisierende (bedingte) Wahrscheinlichkeitsdichtefunktion $p_{\hspace{0.02cm}r(t)\hspace{0.05cm} \vert \hspace{0.05cm}s(t)}$ weiter zu spezifizieren. Im Folgenden wird die Kommunikation über den [[Modulationsverfahren/Qualit%C3%A4tskriterien#Einige_Anmerkungen_zum_AWGN.E2.80.93Kanalmodell| AWGN–Kanal]] betrachtet, dessen wichtigste Eigenschaften hier nochmals kurz zusammengestellt werden: | + | Um weitere Aussagen über die Art der optimalen Messungen des Vektors $\boldsymbol{ r}$ machen zu können, ist es notwendig, die den Kanal charakterisierende (bedingte) Wahrscheinlichkeitsdichtefunktion $p_{\hspace{0.02cm}r(t)\hspace{0.05cm} \vert \hspace{0.05cm}s(t)}$ weiter zu spezifizieren. Im Folgenden wird die Kommunikation über den [[Modulationsverfahren/Qualit%C3%A4tskriterien#Einige_Anmerkungen_zum_AWGN.E2.80.93Kanalmodell| "AWGN–Kanal"]] betrachtet, dessen wichtigste Eigenschaften hier nochmals kurz zusammengestellt werden: |
− | *Das Ausgangssignal des AWGN–Kanals ist $r(t) = s(t)+n(t)$, wobei $s(t)$ das Sendesignal angibt und $n(t)$ durch einen Gaußschen Rauschprozess dargestellt wird.<br> | + | *Das Ausgangssignal des AWGN–Kanals ist $r(t) = s(t)+n(t)$, wobei $s(t)$ das Sendesignal angibt und $n(t)$ durch einen Gaußschen Rauschprozess dargestellt wird.<br> |
− | * | + | *Einen Zufallsprozess $\{n(t)\}$ bezeichnet man als gaußisch, wenn die Elemente der $k$–dimensionalen Zufallsvariablen $\{n_1(t)\hspace{0.05cm} \text{...} \hspace{0.05cm}n_k(t)\}$ gemeinsam gaußverteilt sind ⇒ "jointly Gaussian".<br> |
− | *Der Mittelwert des AWGN–Rauschens ist ${\rm E}\big[n(t)\big] = 0$. Außerdem ist $n(t)$ & | + | *Der Mittelwert des AWGN–Rauschens ist ${\rm E}\big[n(t)\big] = 0$. Außerdem ist $n(t)$ "weiß", was bedeutet, dass das [[Stochastische_Signaltheorie/Leistungsdichtespektrum_(LDS)|Leistungsdichtespektrum]] $\rm (LDS)$ für alle Frequenzen $($von $-\infty$ bis $+\infty)$ konstant ist: |
− | \hspace{0.05cm}.$ | + | :$${\it \Phi}_n(f) = {N_0}/{2} |
+ | \hspace{0.05cm}.$$ | ||
− | *Nach dem [[Stochastische_Signaltheorie/Leistungsdichtespektrum_(LDS)#Theorem_von_Wiener-Chintchine |Wiener–Chintchine–Theorem]] ergibt sich die Autokorrelationsfunktion (AKF) als die [[Signaldarstellung/Fouriertransformation_und_-r%C3%BCcktransformation#Das_zweite_Fourierintegral| Fourierrücktransformierte]] von ${\it \Phi_n(f)}$: | + | *Nach dem [[Stochastische_Signaltheorie/Leistungsdichtespektrum_(LDS)#Theorem_von_Wiener-Chintchine |Wiener–Chintchine–Theorem]] ergibt sich die Autokorrelationsfunktion $\rm (AKF)$ als die [[Signaldarstellung/Fouriertransformation_und_-r%C3%BCcktransformation#Das_zweite_Fourierintegral| Fourierrücktransformierte]] von ${\it \Phi_n(f)}$: |
:$${\varphi_n(\tau)} = {\rm E}\big [n(t) \cdot n(t+\tau)\big ] = {N_0}/{2} \cdot \delta(t)\hspace{0.3cm} | :$${\varphi_n(\tau)} = {\rm E}\big [n(t) \cdot n(t+\tau)\big ] = {N_0}/{2} \cdot \delta(t)\hspace{0.3cm} | ||
\Rightarrow \hspace{0.3cm} {\rm E}\big [n(t) \cdot n(t+\tau)\big ] = | \Rightarrow \hspace{0.3cm} {\rm E}\big [n(t) \cdot n(t+\tau)\big ] = | ||
\left\{ \begin{array}{c} \rightarrow \infty \\ | \left\{ \begin{array}{c} \rightarrow \infty \\ | ||
0 \end{array} \right.\quad | 0 \end{array} \right.\quad | ||
− | \begin{array}{*{1}c} {\rm f{\rm \ddot{u}r}} \hspace{0. | + | \begin{array}{*{1}c} {\rm f{\rm \ddot{u}r}} \hspace{0.25cm} \tau = 0 \hspace{0.05cm}, |
− | \\ {\rm f{\rm \ddot{u}r}} \hspace{0. | + | \\ {\rm f{\rm \ddot{u}r}} \hspace{0.25cm} \tau \ne 0 \hspace{0.05cm},\\ \end{array}$$ |
− | *$N_0$ gibt dabei die physikalische (nur für $f \ge 0$ definierte) Rauschleistungsdichte an. Der konstante LDS–Wert $(N_0/2)$ und das Gewicht der Diracfunktion in der AKF (ebenfalls $N_0/2$ | + | *$N_0$ gibt dabei die physikalische $($nur für $f \ge 0$ definierte$)$ Rauschleistungsdichte an. Der konstante LDS–Wert $(N_0/2)$ und das Gewicht der Diracfunktion in der AKF $($ebenfalls $N_0/2)$ ergibt sich allein durch die zweiseitige Betrachtungsweise.<br><br> |
− | Weitere Informationen | + | ⇒ Weitere Informationen zu diesem Thema liefert das Lernvideo [[Der_AWGN-Kanal_(Lernvideo)|"Der AWGN-Kanal"]] im zweiten Teil.<br> |
== Beschreibung des AWGN-Kanals durch orthonormale Basisfunktionen== | == Beschreibung des AWGN-Kanals durch orthonormale Basisfunktionen== | ||
<br> | <br> | ||
− | Aus dem vorletzten Statement auf der letzten Seite geht hervor, dass | + | Aus dem vorletzten Statement auf der letzten Seite geht hervor, dass |
− | *reines AWGN–Rauschen $n(t)$ stets eine unendliche Varianz (Leistung) aufweist: $\sigma_n^2 \to \infty$,<br> | + | *reines AWGN–Rauschen $n(t)$ stets eine unendliche Varianz (Leistung) aufweist: $\sigma_n^2 \to \infty$,<br> |
− | *in der Realität demzufolge nur gefiltertes Rauschen $n\hspace{0.05cm}'(t) = n(t) \star h_n(t)$ auftreten kann.<br><br> | + | |
+ | *in der Realität demzufolge nur gefiltertes Rauschen $n\hspace{0.05cm}'(t) = n(t) \star h_n(t)$ auftreten kann.<br><br> | ||
− | Mit der Impulsantwort $h_n(t)$ und dem Frequenzgang $H_n(f) = {\rm F}\big [h_n(t)\big ]$ gelten dann folgende Gleichungen:<br> | + | Mit der Impulsantwort $h_n(t)$ und dem Frequenzgang $H_n(f) = {\rm F}\big [h_n(t)\big ]$ gelten dann folgende Gleichungen:<br> |
:$${\rm E}\big[n\hspace{0.05cm}'(t) \big] \hspace{0.15cm} = \hspace{0.2cm} {\rm E}\big[n(t) \big] = 0 \hspace{0.05cm},$$ | :$${\rm E}\big[n\hspace{0.05cm}'(t) \big] \hspace{0.15cm} = \hspace{0.2cm} {\rm E}\big[n(t) \big] = 0 \hspace{0.05cm},$$ | ||
Zeile 206: | Zeile 215: | ||
:$$\sigma_n^2 \hspace{0.1cm} = \hspace{0.1cm} { \varphi_{n\hspace{0.05cm}'}(\tau = 0)} = {N_0}/{2} \cdot | :$$\sigma_n^2 \hspace{0.1cm} = \hspace{0.1cm} { \varphi_{n\hspace{0.05cm}'}(\tau = 0)} = {N_0}/{2} \cdot | ||
\int_{-\infty}^{+\infty}h_n^2(t)\,{\rm d} t ={N_0}/{2}\hspace{0.1cm} \cdot < \hspace{-0.1cm}h_n(t), \hspace{0.1cm} h_n(t) \hspace{-0.05cm} > \hspace{0.1cm} = | \int_{-\infty}^{+\infty}h_n^2(t)\,{\rm d} t ={N_0}/{2}\hspace{0.1cm} \cdot < \hspace{-0.1cm}h_n(t), \hspace{0.1cm} h_n(t) \hspace{-0.05cm} > \hspace{0.1cm} = | ||
− | \int_{-\infty}^{+\infty}{\it \ | + | \int_{-\infty}^{+\infty}{\it \Phi}_{n\hspace{0.05cm}'}(f)\,{\rm d} f = {N_0}/{2} \cdot \int_{-\infty}^{+\infty}|H_n(f)|^2\,{\rm d} f \hspace{0.05cm}.$$ |
− | Im Folgenden beinhaltet $n(t)$ stets implizit eine ''Bandbegrenzung''; auf die Schreibweise $n'(t)$ wird also zukünftig verzichtet.<br> | + | Im Folgenden beinhaltet $n(t)$ stets implizit eine '''Bandbegrenzung'''; auf die Schreibweise $n'(t)$ wird also zukünftig verzichtet.<br> |
{{BlaueBox|TEXT= | {{BlaueBox|TEXT= | ||
− | $\text{ | + | $\text{Beachten Sie:}$ Ähnlich wie das Sendesignal $s(t)$ lässt sich auch der Rauschprozess $\{n(t)\}$ als gewichtete Summe orthonormaler Basisfunktionen $\varphi_j(t)$ schreiben. |
− | *Im Gegensatz zu $s(t)$ ist allerdings | + | *Im Gegensatz zu $s(t)$ ist nun allerdings eine Beschränkung auf eine endliche Anzahl an Basisfunktionen nicht möglich. |
− | *Vielmehr gilt bei rein stochastischen Größen für die entsprechende Signaldarstellung | + | |
+ | *Vielmehr gilt bei rein stochastischen Größen für die entsprechende Signaldarstellung stets | ||
:$$n(t) = \lim_{N \rightarrow \infty} \sum\limits_{j = 1}^{N}n_j \cdot \varphi_j(t) \hspace{0.05cm},$$ | :$$n(t) = \lim_{N \rightarrow \infty} \sum\limits_{j = 1}^{N}n_j \cdot \varphi_j(t) \hspace{0.05cm},$$ | ||
− | wobei der Koeffizient $n_j$ durch die Projektion von $n(t)$ auf die Basisfunktion $\varphi_j(t)$ bestimmt ist: | + | :wobei der Koeffizient $n_j$ durch die Projektion von $n(t)$ auf die Basisfunktion $\varphi_j(t)$ bestimmt ist: |
:$$n_j = \hspace{0.1cm} < \hspace{-0.1cm}n(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.05cm}.$$}} | :$$n_j = \hspace{0.1cm} < \hspace{-0.1cm}n(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.05cm}.$$}} | ||
− | < | + | <u>Hinweis:</u> Um eine Verwechslung mit den Basisfunktionen $\varphi_j(t)$ zu vermeiden, wird im Folgenden die AKF $\varphi_n(\tau)$ des Rauschprozesses stets nur noch als der Erwartungswert ausgedrückt:<br> |
+ | :$${\rm E}\big [n(t) \cdot n(t + \tau)\big ].$$ | ||
== Optimaler Empfänger für den AWGN-Kanal== | == Optimaler Empfänger für den AWGN-Kanal== | ||
<br> | <br> | ||
[[Datei:P ID2005 Dig T 4 2 S5b version1.png|right|frame|Optimaler Empfänger beim AWGN-Kanal|class=fit]] | [[Datei:P ID2005 Dig T 4 2 S5b version1.png|right|frame|Optimaler Empfänger beim AWGN-Kanal|class=fit]] | ||
− | Auch das Empfangssignal $r(t) = s(t) + n(t)$ lässt sich in bekannter Weise in Basisfunktionen zerlegen: | + | Auch das Empfangssignal $r(t) = s(t) + n(t)$ lässt sich in bekannter Weise in Basisfunktionen zerlegen: |
+ | |||
:$$r(t) = \sum\limits_{j = 1}^{\infty}r_j \cdot \varphi_j(t) \hspace{0.05cm}.$$ | :$$r(t) = \sum\limits_{j = 1}^{\infty}r_j \cdot \varphi_j(t) \hspace{0.05cm}.$$ | ||
Zu berücksichtigen ist: | Zu berücksichtigen ist: | ||
− | *Die $M$ möglichen Sendesignale $\{s_i(t)\}$ spannen einen Signalraum mit insgesamt $N$ Basisfunktionen $\varphi_1(t)$, ... , $\varphi_N(t)$ auf. | + | *Die $M$ möglichen Sendesignale $\{s_i(t)\}$ spannen einen Signalraum mit insgesamt $N$ Basisfunktionen $\varphi_1(t)$, ... , $\varphi_N(t)$ auf.<br> |
− | *Diese $N$ Basisfunktionen $\varphi_j(t)$ werden gleichzeitig zur Beschreibung des Rauschsignals $n(t)$ und des Empfangssignals $r(t)$ verwendet. | + | *Diese $N$ Basisfunktionen $\varphi_j(t)$ werden gleichzeitig zur Beschreibung des Rauschsignals $n(t)$ und des Empfangssignals $r(t)$ verwendet.<br> |
− | *Zur vollständigen Charakterisierung von $n(t)$ bzw. $r(t)$ werden nun aber | + | *Zur vollständigen Charakterisierung von $n(t)$ bzw. $r(t)$ werden nun aber auch noch unendlich viele weitere Basisfunktionen $\varphi_{N+1}(t)$, $\varphi_{N+2}(t)$, ... benötigt.<br> |
− | Damit ergeben sich die Koeffizienten des Empfangssignals $r(t)$ gemäß folgender Gleichung | + | |
+ | Damit ergeben sich die Koeffizienten des Empfangssignals $r(t)$ gemäß folgender Gleichung. Hierbei ist berücksichtigt, dass die (möglichen) Sendesignale $s_i(t)$ und das Rauschen $n(t)$ voneinander unabhängig sind: | ||
:$$r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.1cm}=\hspace{0.1cm} | :$$r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.1cm}=\hspace{0.1cm} | ||
Zeile 247: | Zeile 260: | ||
Somit ergibt sich für den optimalen Empfänger die oben skizzierte Struktur.<br> | Somit ergibt sich für den optimalen Empfänger die oben skizzierte Struktur.<br> | ||
<br> | <br> | ||
− | Betrachten wir zunächst den '''AWGN–Kanal'''. Hier kann auf das Vorfilter mit dem Frequenzgang $W(f)$ verzichtet werden, das für farbiges Rauschen vorgesehen ist.<br> | + | Betrachten wir zunächst den '''AWGN–Kanal'''. Hier kann auf das Vorfilter mit dem Frequenzgang $W(f)$ verzichtet werden, das für farbiges Rauschen vorgesehen ist.<br> |
− | Der Detektor des optimalen Empfängers bildet die Koeffizienten $r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t)\hspace{-0.05cm} >$ und reicht diese an den Entscheider weiter. Basiert die Entscheidung auf sämtlichen & | + | #Der Detektor des optimalen Empfängers bildet die Koeffizienten $r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t)\hspace{-0.05cm} >$ und reicht diese an den Entscheider weiter. |
− | + | #Basiert die Entscheidung auf sämtlichen Koeffizienten $r_j$, so ist die Wahrscheinlichkeit für eine Fehlentscheidung minimal und der Empfänger optimal.<br> | |
− | Die reellwertigen Koeffizienten $r_j$ wurden oben wie folgt berechnet: | + | #Die reellwertigen Koeffizienten $r_j$ wurden oben wie folgt berechnet: |
− | :$$r_j = | + | ::$$r_j = |
\left\{ \begin{array}{c} s_{ij} + n_j\\ | \left\{ \begin{array}{c} s_{ij} + n_j\\ | ||
n_j \end{array} \right.\quad | n_j \end{array} \right.\quad | ||
Zeile 258: | Zeile 271: | ||
\\ {j > N} \hspace{0.05cm}.\\ \end{array}$$ | \\ {j > N} \hspace{0.05cm}.\\ \end{array}$$ | ||
− | Nach dem [[Digitalsignalübertragung/Struktur_des_optimalen_Empfängers#Das_Theorem_der_Irrelevanz|Theorem der Irrelevanz]] lässt sich zeigen, dass für additives weißes Gaußsches Rauschen | + | Nach dem [[Digitalsignalübertragung/Struktur_des_optimalen_Empfängers#Das_Theorem_der_Irrelevanz|"Theorem der Irrelevanz"]] lässt sich zeigen, dass für additives weißes Gaußsches Rauschen |
− | *die Optimalität nicht herabgesetzt wird, wenn man die | + | *die Optimalität nicht herabgesetzt wird, wenn man die Koeffizienten $r_{N+1}$, $r_{N+2}$, ... , die nicht von der Nachricht $(s_{ij})$ abhängigen, nicht in den Entscheidungsprozess einbindet, und somit <br> |
− | *der Detektor nur die Projektionen des Empfangssignals $r(t)$ auf die $N$ durch das Nutzsignal $s(t)$ vorgegebenen Basisfunktionen $\varphi_{1}(t)$, ... , $\varphi_{N}(t)$ bilden muss.<br> | + | *der Detektor nur die Projektionen des Empfangssignals $r(t)$ auf die $N$ durch das Nutzsignal $s(t)$ vorgegebenen Basisfunktionen $\varphi_{1}(t)$, ... , $\varphi_{N}(t)$ bilden muss.<br> |
In der Grafik ist diese signifikante Vereinfachung durch die graue Hinterlegung angedeutet.<br> | In der Grafik ist diese signifikante Vereinfachung durch die graue Hinterlegung angedeutet.<br> | ||
− | + | Bei '''farbigem Rauschen''' ⇒ Leistungsdichtespektrum ${\it \Phi}_n(f) \ne {\rm const.}$ ist lediglich zusätzlich ein Vorfilter mit Amplitudengang $|W(f)| = {1}/{\sqrt{\it \Phi}_n(f)}$ erforderlich. | |
+ | #Man nennt dieses Filter auch "Whitening Filter", da die Rauschleistungsdichte am Ausgang wieder konstant – also „weiß” – ist. | ||
+ | #Genaueres hierzu finden Sie im Kapitel [[Stochastische_Signaltheorie/Matched-Filter#Verallgemeinertes_Matched-Filter_f.C3.BCr_den_Fall_farbiger_St.C3.B6rungen|"Matched-Filter bei farbigen Störungen"]] des Buches „Stochastische Signaltheorie”.<br> | ||
== Implementierungsaspekte == | == Implementierungsaspekte == | ||
<br> | <br> | ||
− | Wesentliche Bestandteile des optimalen Empfängers sind die Berechnungen der inneren Produkte gemäß den Gleichungen $r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} >$. Diese können auf verschiedene Art und Weise implementiert werden: | + | Wesentliche Bestandteile des optimalen Empfängers sind die Berechnungen der inneren Produkte gemäß den Gleichungen $r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} >$. |
− | *Beim '''Korrelationsempfänger''' (Näheres zu dieser Implementierung finden Sie im [[Digitalsignal%C3%BCbertragung/Optimale Empf%C3%A4ngerstrategien#Korrelationsempf.C3.A4nger|gleichnamigen Kapitel]]) werden die inneren Produkte direkt entsprechend der Definition mit analogen Multiplizierern und Integratoren realisiert: | + | |
+ | {{BlaueBox|TEXT= | ||
+ | $\text{Diese können auf verschiedene Art und Weise implementiert werden:}$ | ||
+ | |||
+ | *Beim '''Korrelationsempfänger''' $($Näheres zu dieser Implementierung finden Sie im [[Digitalsignal%C3%BCbertragung/Optimale Empf%C3%A4ngerstrategien#Korrelationsempf.C3.A4nger|gleichnamigen Kapitel]]$)$ werden die inneren Produkte direkt entsprechend der Definition mit analogen Multiplizierern und Integratoren realisiert: | ||
:$$r_j = \int_{-\infty}^{+\infty}r(t) \cdot \varphi_j(t) \,{\rm d} t \hspace{0.05cm}.$$ | :$$r_j = \int_{-\infty}^{+\infty}r(t) \cdot \varphi_j(t) \,{\rm d} t \hspace{0.05cm}.$$ | ||
− | *Der '''Matched–Filter–Empfänger''', der bereits im Kapitel [[Digitalsignal%C3%BCbertragung/Fehlerwahrscheinlichkeit_bei_Basisband%C3%BCbertragung#Optimaler_Bin.C3.A4rempf.C3.A4nger_-_Realisierung_mit_Matched-Filter|Optimaler Binärempfänger]] zu Beginn dieses Buches hergeleitet wurde, erzielt mit einem linearen Filter mit der Impulsantwort $h_j(t) = \varphi_j(t) \cdot (T-t)$ und anschließender Abtastung zum Zeitpunkt $t = T$ das gleiche Ergebnis: | + | *Der '''Matched–Filter–Empfänger''', der bereits im Kapitel [[Digitalsignal%C3%BCbertragung/Fehlerwahrscheinlichkeit_bei_Basisband%C3%BCbertragung#Optimaler_Bin.C3.A4rempf.C3.A4nger_-_Realisierung_mit_Matched-Filter|"Optimaler Binärempfänger"]] zu Beginn dieses Buches hergeleitet wurde, erzielt mit einem linearen Filter mit der Impulsantwort $h_j(t) = \varphi_j(t) \cdot (T-t)$ und anschließender Abtastung zum Zeitpunkt $t = T$ das gleiche Ergebnis: |
:$$r_j = \int_{-\infty}^{+\infty}r(\tau) \cdot h_j(t-\tau) \,{\rm d} \tau | :$$r_j = \int_{-\infty}^{+\infty}r(\tau) \cdot h_j(t-\tau) \,{\rm d} \tau | ||
= \int_{-\infty}^{+\infty}r(\tau) \cdot \varphi_j(T-t+\tau) \,{\rm d} \tau \hspace{0.3cm} | = \int_{-\infty}^{+\infty}r(\tau) \cdot \varphi_j(T-t+\tau) \,{\rm d} \tau \hspace{0.3cm} | ||
\Rightarrow \hspace{0.3cm} r_j (t = \tau) = \int_{-\infty}^{+\infty}r(\tau) \cdot \varphi_j(\tau) \,{\rm d} \tau = r_j | \Rightarrow \hspace{0.3cm} r_j (t = \tau) = \int_{-\infty}^{+\infty}r(\tau) \cdot \varphi_j(\tau) \,{\rm d} \tau = r_j | ||
\hspace{0.05cm}.$$ | \hspace{0.05cm}.$$ | ||
+ | [[Datei:P ID2008 Dig T 4 2 S6 version1.png|left|frame|Drei unterschiedliche Implementierungen des inneren Produktes|class=fit]] | ||
+ | <br><br><br> | ||
+ | Die Abbildung zeigt die beiden möglichen Realisierungsformen des optimalen Detektors.}} | ||
− | |||
− | |||
− | |||
== Wahrscheinlichkeitsdichtefunktion der Empfangswerte == | == Wahrscheinlichkeitsdichtefunktion der Empfangswerte == | ||
<br> | <br> | ||
− | Bevor wir uns im folgenden Kapitel der optimalen Gestaltung des Entscheiders und der Berechnung und Annäherung der Fehlerwahrscheinlichkeit zuwenden, erfolgt zunächst eine für den AWGN–Kanal gültige statistische Analyse der Entscheidungsgrößen $r_j$. | + | Bevor wir uns im folgenden Kapitel der optimalen Gestaltung des Entscheiders und der Berechnung und Annäherung der Fehlerwahrscheinlichkeit zuwenden, erfolgt zunächst eine für den AWGN–Kanal gültige statistische Analyse der Entscheidungsgrößen $r_j$. |
[[Datei:P ID2009 Dig T 4 2 S7 version1.png|right|frame|Signalraumkonstellation und WDF des Empfangssignals|class=fit]] | [[Datei:P ID2009 Dig T 4 2 S7 version1.png|right|frame|Signalraumkonstellation und WDF des Empfangssignals|class=fit]] | ||
− | Dazu betrachten wir nochmals den optimalen Binärempfänger für die bipolare Basisbandübertragung über den AWGN–Kanal, wobei wir von der für das vierte Hauptkapitel gültigen Beschreibungsform ausgehen. | + | Dazu betrachten wir nochmals den optimalen Binärempfänger für die bipolare Basisbandübertragung über den AWGN–Kanal, wobei wir von der für das vierte Hauptkapitel gültigen Beschreibungsform ausgehen. |
+ | |||
+ | Mit den Parametern $N = 1$ und $M = 2$ ergibt sich für das Sendesignal die in der linken Grafik dargestellte Signalraumkonstellation | ||
+ | *mit nur einer Basisfunktion $\varphi_1(t)$, wegen $N = 1$,<br> | ||
− | + | *mit den beiden Signalraumpunkten $s_i \in \{s_0, \hspace{0.05cm}s_1\}$, wegen $M = 2$. | |
− | |||
− | *mit den beiden Signalraumpunkten $s_i \in \{s_0, \hspace{0.05cm}s_1\}$, wegen $M = 2$. | ||
<br clear=all> | <br clear=all> | ||
− | Für das Signal $r(t) = s(t) + n(t)$ am Ausgang des AWGN–Kanals ergibt sich im rauschfreien Fall ⇒ $r(t) = s(t)$ die genau gleiche Konstellation. Die Signalraumpunkte liegen somit bei | + | Für das Signal $r(t) = s(t) + n(t)$ am Ausgang des AWGN–Kanals ergibt sich im rauschfreien Fall ⇒ $r(t) = s(t)$ die genau gleiche Konstellation. Die Signalraumpunkte liegen somit bei |
:$$r_0 = s_0 = \sqrt{E}\hspace{0.05cm},\hspace{0.2cm}r_1 = s_1 = -\sqrt{E}\hspace{0.05cm}.$$ | :$$r_0 = s_0 = \sqrt{E}\hspace{0.05cm},\hspace{0.2cm}r_1 = s_1 = -\sqrt{E}\hspace{0.05cm}.$$ | ||
− | Bei Berücksichtigung des (bandbegrenzten) AWGN–Rauschens $n(t)$ überlagern sich den beiden Punkten $r_0$ und $r_1$ jeweils Gaußkurven mit der Varianz $\sigma_n^2$ ⇒ Streuung $\sigma_n$ (siehe rechte Grafik). Die WDF der Rauschkomponente $n(t)$ lautet dabei: | + | *Bei Berücksichtigung des (bandbegrenzten) AWGN–Rauschens $n(t)$ überlagern sich den beiden Punkten $r_0$ und $r_1$ jeweils Gaußkurven mit der Varianz $\sigma_n^2$ ⇒ Streuung $\sigma_n$ $($siehe rechte Grafik$)$. |
+ | |||
+ | *Die WDF der Rauschkomponente $n(t)$ lautet dabei: | ||
:$$p_n(n) = \frac{1}{\sqrt{2\pi} \cdot \sigma_n}\cdot {\rm e}^{ - {n^2}/(2 \sigma_n^2)}\hspace{0.05cm}.$$ | :$$p_n(n) = \frac{1}{\sqrt{2\pi} \cdot \sigma_n}\cdot {\rm e}^{ - {n^2}/(2 \sigma_n^2)}\hspace{0.05cm}.$$ | ||
− | Für die bedingte Wahrscheinlichkeitsdichte, dass der Empfangswert $\rho$ anliegt, wenn $s_i$ gesendet wurde, ergibt sich dann folgender Ausdruck: | + | Für die bedingte Wahrscheinlichkeitsdichte, dass der Empfangswert $\rho$ anliegt, wenn $s_i$ gesendet wurde, ergibt sich dann folgender Ausdruck: |
:$$p_{\hspace{0.02cm}r\hspace{0.05cm}|\hspace{0.05cm}s}(\rho\hspace{0.05cm}|\hspace{0.05cm}s_i) = \frac{1}{\sqrt{2\pi} \cdot \sigma_n}\cdot {\rm e}^{ - {(\rho - s_i)^2}/(2 \sigma_n^2)} \hspace{0.05cm}.$$ | :$$p_{\hspace{0.02cm}r\hspace{0.05cm}|\hspace{0.05cm}s}(\rho\hspace{0.05cm}|\hspace{0.05cm}s_i) = \frac{1}{\sqrt{2\pi} \cdot \sigma_n}\cdot {\rm e}^{ - {(\rho - s_i)^2}/(2 \sigma_n^2)} \hspace{0.05cm}.$$ | ||
Zu den Einheiten der hier aufgeführten Größen ist zu bemerken: | Zu den Einheiten der hier aufgeführten Größen ist zu bemerken: | ||
− | *$r_0 = s_0$ und $r_1 = s_1$ sowie $n$ sind jeweils Skalare mit der Einheit „Wurzel aus Energie”.<br> | + | *$r_0 = s_0$ und $r_1 = s_1$ sowie $n$ sind jeweils Skalare mit der Einheit „Wurzel aus Energie”.<br> |
− | *Damit ist offensichtlich, dass $\sigma_n$ ebenfalls die Einheit „Wurzel aus Energie” besitzt und $\sigma_n^2$ eine Energie darstellt.<br> | + | *Damit ist offensichtlich, dass $\sigma_n$ ebenfalls die Einheit „Wurzel aus Energie” besitzt und $\sigma_n^2$ eine Energie darstellt.<br> |
− | *Beim AWGN–Kanal ist die Rauschvarianz $\sigma_n^2 = N_0/2$. Diese ist also ebenfalls eine physikalische Größe mit der Einheit $\rm W/Hz = Ws$.<br><br> | + | *Beim AWGN–Kanal ist die Rauschvarianz $\sigma_n^2 = N_0/2$. Diese ist also ebenfalls eine physikalische Größe mit der Einheit $\rm W/Hz = Ws$.<br><br> |
− | Die hier angesprochene Thematik wird in der [[Aufgaben:Aufgabe_4.06:_Optimale_Entscheidungsgrenzen|Aufgabe 4.6]] an Beispielen verdeutlicht.<br> | + | Die hier angesprochene Thematik wird in der [[Aufgaben:Aufgabe_4.06:_Optimale_Entscheidungsgrenzen|"Aufgabe 4.6"]] an Beispielen verdeutlicht.<br> |
== N–dimensionales Gaußsches Rauschen== | == N–dimensionales Gaußsches Rauschen== | ||
<br> | <br> | ||
− | Liegt ein $N$–dimensionales Modulationsverfahren vor, das heißt, | + | Liegt ein $N$–dimensionales Modulationsverfahren vor, das heißt, wenn mit $0 \le i \le M-1$ und $1 \le j \le N$ |
:$$s_i(t) = \sum\limits_{j = 1}^{N} s_{ij} \cdot \varphi_j(t) = s_{i1} \cdot \varphi_1(t) | :$$s_i(t) = \sum\limits_{j = 1}^{N} s_{ij} \cdot \varphi_j(t) = s_{i1} \cdot \varphi_1(t) | ||
+ s_{i2} \cdot \varphi_2(t) + \hspace{0.05cm}\text{...}\hspace{0.05cm} + s_{iN} \cdot \varphi_N(t)\hspace{0.05cm}\hspace{0.3cm} | + s_{i2} \cdot \varphi_2(t) + \hspace{0.05cm}\text{...}\hspace{0.05cm} + s_{iN} \cdot \varphi_N(t)\hspace{0.05cm}\hspace{0.3cm} | ||
Zeile 322: | Zeile 344: | ||
\hspace{0.05cm},$$ | \hspace{0.05cm},$$ | ||
− | + | gilt, dann muss der Rauschvektor $\boldsymbol{ n}$ ebenfalls mit Dimension $N$ angesetzt werden. Das gleiche gilt auch für den Empfangsvektor $\boldsymbol{ r}$: | |
:$$\boldsymbol{ n} = \left(n_{1}, n_{2}, \hspace{0.05cm}\text{...}\hspace{0.05cm}, n_{N}\right ) | :$$\boldsymbol{ n} = \left(n_{1}, n_{2}, \hspace{0.05cm}\text{...}\hspace{0.05cm}, n_{N}\right ) | ||
− | \hspace{0.01cm}, | + | \hspace{0.01cm},$$ |
+ | :$$\boldsymbol{ r} = \left(r_{1}, r_{2}, \hspace{0.05cm}\text{...}\hspace{0.05cm}, r_{N}\right )\hspace{0.05cm}.$$ | ||
− | Die Wahrscheinlichkeitsdichtefunktion (WDF) lautet dann für den AWGN–Kanal mit der Realisierung $\boldsymbol{ \eta}$ des Rauschsignals | + | Die Wahrscheinlichkeitsdichtefunktion $\rm (WDF)$ lautet dann für den AWGN–Kanal mit der Realisierung $\boldsymbol{ \eta}$ des Rauschsignals |
:$$p_{\boldsymbol{ n}}(\boldsymbol{ \eta}) = \frac{1}{\left( \sqrt{2\pi} \cdot \sigma_n \right)^N } \cdot | :$$p_{\boldsymbol{ n}}(\boldsymbol{ \eta}) = \frac{1}{\left( \sqrt{2\pi} \cdot \sigma_n \right)^N } \cdot | ||
{\rm exp} \left [ - \frac{|| \boldsymbol{ \eta} ||^2}{2 \sigma_n^2}\right ]\hspace{0.05cm},$$ | {\rm exp} \left [ - \frac{|| \boldsymbol{ \eta} ||^2}{2 \sigma_n^2}\right ]\hspace{0.05cm},$$ | ||
Zeile 336: | Zeile 359: | ||
{\rm exp} \left [ - \frac{|| \boldsymbol{ \rho} - \boldsymbol{ s}_i ||^2}{2 \sigma_n^2}\right ]\hspace{0.05cm}.$$ | {\rm exp} \left [ - \frac{|| \boldsymbol{ \rho} - \boldsymbol{ s}_i ||^2}{2 \sigma_n^2}\right ]\hspace{0.05cm}.$$ | ||
− | Die Gleichung ergibt sich aus der allgemeinen Darstellung der $N$–dimensionalen Gaußschen WDF im Abschnitt [[Stochastische_Signaltheorie/Verallgemeinerung_auf_N-dimensionale_Zufallsgr%C3%B6%C3%9Fen#Korrelationsmatrix|Korrelationsmatrix]] des Buches „Stochastische Signaltheorie” unter der Voraussetzung, dass die Komponenten unkorreliert (und somit statistisch unabhängig) sind. $||\boldsymbol{ \eta}||$ bezeichnet man als die | + | Die Gleichung ergibt sich |
+ | *aus der allgemeinen Darstellung der $N$–dimensionalen Gaußschen WDF im Abschnitt [[Stochastische_Signaltheorie/Verallgemeinerung_auf_N-dimensionale_Zufallsgr%C3%B6%C3%9Fen#Korrelationsmatrix|"Korrelationsmatrix"]] des Buches „Stochastische Signaltheorie” | ||
+ | |||
+ | *unter der Voraussetzung, dass die Komponenten unkorreliert (und somit statistisch unabhängig) sind. | ||
+ | |||
+ | *$||\boldsymbol{ \eta}||$ bezeichnet man als die "Norm" (Länge) des Vektors $\boldsymbol{ \eta}$.<br> | ||
+ | |||
− | |||
{{GraueBox|TEXT= | {{GraueBox|TEXT= | ||
− | $\text{Beispiel | + | $\text{Beispiel 3:}$ |
− | Rechts dargestellt ist die zweidimensionale Gauß–WDF $p_{\boldsymbol{ n} } (\boldsymbol{ \eta})$ der 2D–Zufallsgröße $\boldsymbol{ n} = (n_1,\hspace{0.05cm}n_2)$. Beliebige Realisierungen der Zufallsgröße $\boldsymbol{ n}$ werden mit $\boldsymbol{ \eta} = (\eta_1,\hspace{0.05cm}\eta_2)$ bezeichnet. | + | Rechts dargestellt ist die zweidimensionale Gauß–WDF $p_{\boldsymbol{ n} } (\boldsymbol{ \eta})$ der 2D–Zufallsgröße $\boldsymbol{ n} = (n_1,\hspace{0.05cm}n_2)$. Beliebige Realisierungen der Zufallsgröße $\boldsymbol{ n}$ werden mit $\boldsymbol{ \eta} = (\eta_1,\hspace{0.05cm}\eta_2)$ bezeichnet. Die Gleichung der dargestellten Glockenkurve lautet: |
− | + | [[Datei:P ID2012 Dig T 4 2 S8 version1.png|right|frame|Zweidimensionale Gauß–WDF]] | |
:$$p_{n_1, n_2}(\eta_1, \eta_2) = \frac{1}{\left( \sqrt{2\pi} \cdot \sigma_n \right)^2 } \cdot | :$$p_{n_1, n_2}(\eta_1, \eta_2) = \frac{1}{\left( \sqrt{2\pi} \cdot \sigma_n \right)^2 } \cdot | ||
{\rm exp} \left [ - \frac{ \eta_1^2 + \eta_2^2}{2 \sigma_n^2}\right ]\hspace{0.05cm}. $$ | {\rm exp} \left [ - \frac{ \eta_1^2 + \eta_2^2}{2 \sigma_n^2}\right ]\hspace{0.05cm}. $$ | ||
− | *Das Maximum dieser Funktion liegt bei $\eta_1 = \eta_2 = 0$ und hat den Wert $2\pi \cdot \sigma_n^2$. | + | *Das Maximum dieser Funktion liegt bei $\eta_1 = \eta_2 = 0$ und hat den Wert $2\pi \cdot \sigma_n^2$. |
− | *Mit $\sigma_n^2 = N_0/2$ lässt sich die 2D–WDF in Vektorform auch wie folgt schreiben: | + | |
+ | *Mit $\sigma_n^2 = N_0/2$ lässt sich die 2D–WDF in Vektorform auch wie folgt schreiben: | ||
:$$p_{\boldsymbol{ n} }(\boldsymbol{ \eta}) = \frac{1}{\pi \cdot N_0 } \cdot | :$$p_{\boldsymbol{ n} }(\boldsymbol{ \eta}) = \frac{1}{\pi \cdot N_0 } \cdot | ||
{\rm exp} \left [ - \frac{\vert \vert \boldsymbol{ \eta} \vert \vert ^2}{N_0}\right ]\hspace{0.05cm}.$$ | {\rm exp} \left [ - \frac{\vert \vert \boldsymbol{ \eta} \vert \vert ^2}{N_0}\right ]\hspace{0.05cm}.$$ | ||
− | *Diese rotationssymmetrische WDF eignet sich zum Beispiel für die Beschreibung/Untersuchung eines | + | *Diese rotationssymmetrische WDF eignet sich zum Beispiel für die Beschreibung/Untersuchung eines zweidimensionalen Modulationsverfahrens wie [[Digitalsignalübertragung/Trägerfrequenzsysteme_mit_kohärenter_Demodulation#Quadraturamplitudenmodulation_.28M.E2.80.93QAM.29|"M–QAM"]], [[Digitalsignalübertragung/Trägerfrequenzsysteme_mit_kohärenter_Demodulation#Mehrstufiges_Phase.E2.80.93Shift_Keying_.28M.E2.80.93PSK.29|"M–PSK"]] oder [[Digitalsignalübertragung/Trägerfrequenzsysteme_mit_kohärenter_Demodulation#Binary_Frequency_Shift_Keying_.282.E2.80.93FSK.29|"2–FSK"]].<br> |
− | *Oft werden zweidimensionale reelle Zufallsgrößen aber auch eindimensional–komplex dargestellt, meist in der Form $n(t) = n_{\rm I}(t) + {\rm j} \cdot n_{\rm Q}(t)$ Die beiden Komponenten bezeichnet man dann als | + | *Oft werden zweidimensionale reelle Zufallsgrößen aber auch eindimensional–komplex dargestellt, meist in der Form $n(t) = n_{\rm I}(t) + {\rm j} \cdot n_{\rm Q}(t)$. Die beiden Komponenten bezeichnet man dann als "Inphaseanteil" $n_{\rm I}(t)$ und "Quadraturanteil" $n_{\rm Q}(t)$ des Rauschens.<br> |
− | *Die Wahrscheinlichkeitsdichtefunktion hängt nur vom Betrag $\vert n(t) \vert$ der Rauschvariablen ab und nicht von Winkel ${\rm arc} \ n(t)$. Das heißt: Komplexes Rauschen ist zirkulär symmetrisch (siehe Grafik).<br> | + | *Die Wahrscheinlichkeitsdichtefunktion hängt nur vom Betrag $\vert n(t) \vert$ der Rauschvariablen ab und nicht von Winkel ${\rm arc} \ n(t)$. Das heißt: Komplexes Rauschen ist zirkulär symmetrisch (siehe Grafik).<br> |
− | *Zirkulär symmetrisch bedeutet auch, dass Inphasekomponente $n_{\rm I}(t)$ und Quadraturkomponente $n_{\rm Q}(t)$ die gleiche Verteilung aufweisen und damit auch gleiche Varianz (Streuung) besitzen: | + | *Zirkulär symmetrisch bedeutet auch, dass die Inphasekomponente $n_{\rm I}(t)$ und die Quadraturkomponente $n_{\rm Q}(t)$ die gleiche Verteilung aufweisen und damit auch gleiche Varianz (Streuung) besitzen: |
− | :$$ {\rm E} \big [ n_{\rm I}^2(t)\big ] | + | :$$ {\rm E} \big [ n_{\rm I}^2(t)\big ] = {\rm E}\big [ n_{\rm Q}^2(t) \big ] = \sigma_n^2 \hspace{0.05cm},\hspace{1cm}{\rm E}\big [ n(t) \cdot n^*(t) \big ]\hspace{0.1cm} = \hspace{0.1cm} {\rm E}\big [ n_{\rm I}^2(t) \big ] + {\rm E}\big [ n_{\rm Q}^2(t)\big ] = 2\sigma_n^2 \hspace{0.05cm}.$$}} |
− | Abschließend noch einige ''Bezeichnungsvarianten'' für Gaußsche Zufallsgrößen: | + | Abschließend noch einige '''Bezeichnungsvarianten''' für Gaußsche Zufallsgrößen: |
− | :$$x ={\cal N}(\mu, \sigma^2) \hspace{-0.1cm}: \hspace{0.3cm}\text{reelle gaußverteilte Zufallsgröße, mit Mittelwert}\hspace{0. | + | :$$x ={\cal N}(\mu, \sigma^2) \hspace{-0.1cm}: \hspace{0.3cm}\text{reelle gaußverteilte Zufallsgröße, mit Mittelwert}\hspace{0.15cm}\mu \text { und Varianz}\hspace{0.15cm}\sigma^2 \hspace{0.05cm},$$ |
:$$y={\cal CN}(\mu, \sigma^2)\hspace{-0.1cm}: \hspace{0.12cm}\text{komplexe gaußverteilte Zufallsgröße} \hspace{0.05cm}.$$ | :$$y={\cal CN}(\mu, \sigma^2)\hspace{-0.1cm}: \hspace{0.12cm}\text{komplexe gaußverteilte Zufallsgröße} \hspace{0.05cm}.$$ | ||
Aktuelle Version vom 18. Juli 2022, 15:31 Uhr
Inhaltsverzeichnis
- 1 Blockschaltbild und Voraussetzungen
- 2 Fundamentaler Ansatz zum optimalen Empfängerentwurf
- 3 Das Theorem der Irrelevanz
- 4 Einige Eigenschaften des AWGN-Kanals
- 5 Beschreibung des AWGN-Kanals durch orthonormale Basisfunktionen
- 6 Optimaler Empfänger für den AWGN-Kanal
- 7 Implementierungsaspekte
- 8 Wahrscheinlichkeitsdichtefunktion der Empfangswerte
- 9 N–dimensionales Gaußsches Rauschen
- 10 Aufgaben zum Kapitel
Blockschaltbild und Voraussetzungen
In diesem Kapitel wird die Struktur des optimalen Empfängers eines digitalen Übertragungssystems sehr allgemein hergeleitet, wobei
- das Modulationsverfahren und weitere Systemdetails nicht weiter spezifiziert werden,
- von den Basisfunktionen und der Signalraumdarstellung gemäß dem Kapitel "Signale, Basisfunktionen und Vektorräume" ausgegangen wird.
Zum obigen Blockschaltbild ist anzumerken:
- Der Symbolumfang der Quelle beträgt $M$ und der Symbolvorrat ist $\{m_i\}$ mit $i = 0$, ... , $M-1$.
- Die zugehörigen Symbolwahrscheinlichkeiten ${\rm Pr}(m = m_i)$ seien auch dem Empfänger bekannt.
- Zur Nachrichtenübertragung stehen $M$ Signalformen $s_i(t)$ zur Verfügung; auch für die Laufvariable gelte die Indizierung $i = 0$, ... , $M-1$.
- Es besteht eine feste Beziehung zwischen den Nachrichten und den Signalen. Wird $m =m_i$ übertragen, so ist das Sendesignal $s(t) =s_i(t)$.
- Lineare Kanalverzerrungen sind in der obigen Grafik durch die Impulsantwort $h(t)$ berücksichtigt. Außerdem ist ein (irgendwie geartetes) Rauschen $n(t)$ wirksam.
- Mit diesen beiden die Übertragung störenden Effekten lässt sich das am Empfänger ankommende Signal $r(t)$ in folgender Weise angeben:
- $$r(t) = s(t) \star h(t) + n(t) \hspace{0.05cm}.$$
- Aufgabe des (optimalen) Empfängers ist es, anhand seines Eingangssignals $r(t)$ herauszufinden, welche der $M$ möglichen Nachrichten $m_i$ ⇒ Signal $s_i(t)$ gesendet wurde. Der vom Empfänger gefundene Schätzwert für $m$ wird durch ein "Zirkumflex" gekennzeichnet ⇒ $\hat{m}$.
$\text{Definition:}$ Man spricht von einem optimalen Empfänger, wenn die Symbolfehlerwahrscheinlichkeit den für die Randbedingungen kleinstmöglichen Wert annimmt:
- $$p_{\rm S} = {\rm Pr} ({\cal E}) = {\rm Pr} ( \hat{m} \ne m) \hspace{0.15cm} \Rightarrow \hspace{0.15cm}{\rm Minimum} \hspace{0.05cm}.$$
Hinweise:
- Im Folgenden wird meist der AWGN–Ansatz ⇒ $r(t) = s(t) + n(t)$ vorausgesetzt, was bedeutet, dass $h(t) = \delta(t)$ als verzerrungsfrei angenommen wird.
- Andernfalls können wir die Signale $s_i(t)$ als ${s_i}'(t) = s_i(t) \star h(t)$ neu definieren, also die deterministischen Kanalverzerrungen dem Sendesignal beaufschlagen.
Fundamentaler Ansatz zum optimalen Empfängerentwurf
Gegenüber dem auf der vorherigen Seite gezeigten "Blockschaltbild" führen wir nun einige wesentliche Verallgemeinerungen durch:
- Der Übertragungskanal wird nun beschrieben durch die "bedingte Wahrscheinlichkeitsdichtefunktion" $p_{\hspace{0.02cm}r(t)\hspace{0.02cm} \vert \hspace{0.02cm}s(t)}$, welche die Abhängigkeit des Empfangssignals $r(t)$ vom Sendesignal $s(t)$ festlegt.
- Wurde nun ein ganz bestimmtes Signal $r(t) = \rho(t)$ empfangen, so hat der Empfänger die Aufgabe, die Wahrscheinlichkeitsdichtefunktionen auf der Grundlage dieser "Signalrealisierung" $\rho(t)$ und der $M$ bedingten Wahrscheinlichkeitsdichtefunktionen zu bestimmen:
- $$p_{\hspace{0.05cm}r(t) \hspace{0.05cm} \vert \hspace{0.05cm} s(t) } (\rho(t) \hspace{0.05cm} \vert \hspace{0.05cm} s_i(t))\hspace{0.5cm}{\rm mit}\hspace{0.5cm} i = 0, \text{...} \hspace{0.05cm}, M-1.$$
- Ermittelt werden soll, welche Nachricht $\hat{m}$ am wahrscheinlichsten gesendet wurde, unter Berücksichtigung aller möglichen Sendesignale $s_i(t)$ und derenn Auftrittsswahrscheinlichkeiten ${\rm Pr}(m = m_i)$.
- Die Schätzung des optimalen Empfängers ist also ganz allgemein bestimmt durch die Gleichung
- $$\hat{m} = {\rm arg} \max_i \hspace{0.1cm} p_{\hspace{0.02cm}s(t) \hspace{0.05cm} \vert \hspace{0.05cm} r(t) } ( s_i(t) \hspace{0.05cm} \vert \hspace{0.05cm} \rho(t)) = {\rm arg} \max_i \hspace{0.1cm} p_{m \hspace{0.05cm} \vert \hspace{0.05cm} r(t) } ( \hspace{0.05cm}m_i\hspace{0.05cm} \vert \hspace{0.05cm}\rho(t))\hspace{0.05cm}.$$
$\text{In anderen Worten:}$ Der optimale Empfänger betrachtet diejenige Nachricht $\hat{m} \in \{m_i\}$ als die am wahrscheinlichsten gesendete Nachricht, deren bedingte Wahrscheinlichkeitsdichtefunktion $p_{\hspace{0.02cm}m \hspace{0.05cm} \vert \hspace{0.05cm} r(t) }$ für das anliegende Empfangssignal $\rho(t)$ sowie unter der Annahme $m=\hat{m}$ den größtmöglichen Wert annimmt.
Bevor wir die diese Entscheidungsregel näher diskutieren, soll der optimale Empfänger entsprechend der Grafik noch in zwei Funktionsblöcke aufgeteilt werden:
- Der Detektor nimmt am Empfangssignal $r(t)$ verschiedene Messungen vor und fasst diese im Vektor $\boldsymbol{r}$ zusammen. Bei $K$ Messungen entspricht $\boldsymbol{r}$ einem Punkt im $K$–dimensionalen Vektorraum.
- Der Entscheider bildet abhängig von diesem Vektor den Schätzwert. Bei einem gegebenen Vektor $\boldsymbol{r} = \boldsymbol{\rho}$ lautet dabei die Entscheidungsregel:
- $$\hat{m} = {\rm arg}\hspace{0.05cm} \max_i \hspace{0.1cm} P_{m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}) \hspace{0.05cm}.$$
Im Gegensatz zur oberen Entscheidungsregel tritt nun eine bedingte Wahrscheinlichkeit $P_{m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} }$ anstelle der bedingten Wahrscheinlichkeitskeitsdichtefunktion $\rm (WDF)$ $p_{m\hspace{0.05cm} \vert \hspace{0.05cm}r(t)}$ auf. Beachten Sie bitte die Groß– bzw. Kleinschreibung für die unterschiedlichen Bedeutungen.
$\text{Beispiel 1:}$ Wir betrachten nun die Funktion $y = {\rm arg}\hspace{0.05cm} \max \ p(x)$, wobei $p(x)$ die Wahrscheinlichkeitsdichtefunktion $\rm (WDF)$ einer wertkontinuierlichen oder wertdiskreten Zufallsgröße $x$ beschreibt. In der rechten Grafik besteht die WDF aus einer Summe von Diracfunktionen mit Wahrscheinlichkeiten als Impulsgewichte.
⇒ Die Grafik zeigt beispielhafte Funktionen. In beiden Fällen liegt das WDF–Maximum $(17)$ bei $x = 6$:
- $$\max_i \hspace{0.1cm} p(x) = 17\hspace{0.05cm},$$
- $$y = {\rm \hspace{0.05cm}arg} \max_i \hspace{0.1cm} p(x) = 6\hspace{0.05cm}.$$
⇒ Die (bedingten) Wahrscheinlichkeiten in der Gleichung
- $$\hat{m} = {\rm arg}\hspace{0.05cm} \max_i \hspace{0.1cm} P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho})$$
sind a–Posteriori–Wahrscheinlichkeiten. Mit dem Satz von Bayes kann hierfür geschrieben werden:
- $$P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}) = \frac{ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )}{p_{\boldsymbol{ r} } (\boldsymbol{\rho})} \hspace{0.05cm}.$$
Der Nennerterm $p_{\boldsymbol{ r} }(\boldsymbol{\rho})$ ist für alle Alternativen $m_i$ gleich und muss für die Entscheidung nicht berücksichtigt werden. Damit erhält man folgende Regeln:
$\text{Theorem:}$ Die Entscheidungsregel des optimalen Empfängers, auch bekannt als Maximum–a–posteriori–Empfänger $($kurz: MAP–Empfänger$)$ lautet:
- $$\hat{m}_{\rm MAP} = {\rm \hspace{0.05cm} arg} \max_i \hspace{0.1cm} P_{\hspace{0.02cm}m\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r} } ( m_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}) = {\rm \hspace{0.05cm}arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm} m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm} m_i )\big ]\hspace{0.05cm}.$$
- Vorteil der zweiten Gleichung ist, dass die Vorwärtsrichtung des Kanals Anwendung findet ⇒ bedingte WDF $p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm} m }$ $($"Ausgang unter der Bedingung Eingang"$)$.
- Dagegen verwendet die erste Gleichung die Rückschlusswahrscheinlichkeiten $P_{\hspace{0.05cm}m\hspace{0.05cm} \vert \hspace{0.02cm} \boldsymbol{ r} } $ $($"Eingang unter der Bedingung Ausgang"$)$.
$\text{Theorem:}$ Der Maximum–Likelihood–Empfänger $($kurz: ML–Empfänger$)$ verwendet die Entscheidungsregel
- $$\hat{m}_{\rm ML} = \hspace{-0.1cm} {\rm arg} \max_i \hspace{0.1cm} p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )\hspace{0.05cm}.$$
- Bei diesem werden die möglicherweise unterschiedlichen Auftrittswahrscheinlichkeiten ${\rm Pr}(m = m_i)$ für den Entscheidungsprozess nicht herangezogen.
- Zum Beispiel, weil diese dem Empfänger nicht bekannt sind.
Im früheren Kapitel "Optimale Empfängerstrategien" finden Sie auch andere Herleitungen für diese Empfängertypen.
$\text{Fazit:}$ Bei gleichwahrscheinlichen Nachrichten $\{m_i\}$ ⇒ ${\rm Pr}(m = m_i) = 1/M$ ist der im Allgemeinen etwas schlechtere "Maximum–Likelihood–Empfänger" gleichwertig mit dem "Maximum–a–posteriori–Empfänger":
- $$\hat{m}_{\rm MAP} = \hat{m}_{\rm ML} =\hspace{-0.1cm} {\rm\hspace{0.05cm} arg} \max_i \hspace{0.1cm} p_{\boldsymbol{ r}\hspace{0.05cm} \vert \hspace{0.05cm}m } (\boldsymbol{\rho}\hspace{0.05cm} \vert \hspace{0.05cm}m_i )\hspace{0.05cm}.$$
Das Theorem der Irrelevanz
Zu beachten ist, dass der auf der letzten Seite beschriebene Empfänger nur dann optimal ist, wenn auch der Detektor bestmöglich implementiert ist, das heißt, wenn durch den Übergang vom kontinuierlichen Signal $r(t)$ zum Vektor $\boldsymbol{r}$ keine Information verloren geht.
Um zu klären, welche und wieviele Messungen am Empfangssignal $r(t)$ nötig sind, um Optimalität zu garantieren, ist das "Theorem der Irrelevanz" hilfreich:
- Dazu betrachten wir den skizzierten Empfänger, dessen Detektor aus dem Empfangssignal $r(t)$ die zwei Vektoren $\boldsymbol{r}_1$ und $\boldsymbol{r}_2$ ableitet und dem Entscheider zur Verfügung stellt.
- Diese Größen stehen mit der Nachricht $ m \in \{m_i\}$ über die Verbundwahrscheinlichkeitsdichte $p_{\boldsymbol{ r}_1, \hspace{0.05cm}\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm}m }$ in Zusammenhang.
- Die Entscheidungsregel des MAP–Empfängers lautet mit Anpassung an dieses Beispiel:
- $$\hat{m}_{\rm MAP} \hspace{-0.1cm} = \hspace{-0.1cm} {\rm arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}_1 , \hspace{0.05cm}\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm}m } \hspace{0.05cm} (\boldsymbol{\rho}_1, \hspace{0.05cm}\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} m_i ) \big]= {\rm arg} \max_i \hspace{0.1cm}\big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m } \hspace{0.05cm} (\boldsymbol{\rho}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_i ) \cdot p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i )\big] \hspace{0.05cm}.$$
- Die Vektoren $\boldsymbol{r}_1$ und $\boldsymbol{r}_2$ sind Zufallsgrößen. Ihre Realisierungen werden hier und im Folgenden mit $\boldsymbol{\rho}_1$ und $\boldsymbol{\rho}_2$ bezeichnet. Zur Hervorhebung sind alle Vektoren in der Grafik rot eingetragen.
- Die Voraussetzungen für die Anwendung des "Theorems der Irrelevanz" sind die gleichen wie die an eine "Markovkette" erster Ordnung. Die Zufallsvariablen $x$, $y$, $z$ formen dann eine Markovkette erster Ordnung, falls die Verteilung von $z$ bei gegebenem $y$ unabhängig von $x$ ist:
- $$p(x, y, z) = p(x) \cdot p(y\hspace{0.05cm} \vert \hspace{0.05cm}x) \cdot p(z\hspace{0.05cm} \vert \hspace{0.05cm}y) \hspace{0.25cm} {\rm anstelle \hspace{0.15cm}von} \hspace{0.25cm}p(x, y, z) = p(x) \cdot p(y\hspace{0.05cm} \vert \hspace{0.05cm}x) \cdot p(z\hspace{0.05cm} \vert \hspace{0.05cm}x, y) \hspace{0.05cm}.$$
- Der optimale Empfänger muss allgemein beide Vektoren $\boldsymbol{r}_1$ und $\boldsymbol{r}_2$ auswerten, da in obiger Entscheidungsregel beide Verbunddichten $p_{\boldsymbol{ r}_1\hspace{0.05cm} \vert \hspace{0.05cm}m }$ und $p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r}_1, \hspace{0.05cm}m }$ auftreten. Dagegen kann der Empfänger ohne Informationseinbuße die zweite Messung vernachlässigen, falls $\boldsymbol{r}_2$ bei gegebenem $\boldsymbol{r}_1$ unabhängig von der Nachricht $m$ ist:
- $$p_{\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i )= p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 } \hspace{0.05cm} (\boldsymbol{\rho}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 ) \hspace{0.05cm}.$$
- In diesem Fall lässt sich die Entscheidungsregel weiter vereinfachen:
- $$\hat{m}_{\rm MAP} = {\rm arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m } \hspace{0.05cm} (\boldsymbol{\rho}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_i ) \cdot p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i ) \big]$$
- $$\Rightarrow \hspace{0.3cm}\hat{m}_{\rm MAP} = {\rm arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m } \hspace{0.05cm} (\boldsymbol{\rho}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_i ) \cdot p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 )\big]$$
- $$\Rightarrow \hspace{0.3cm}\hat{m}_{\rm MAP} = {\rm arg} \max_i \hspace{0.1cm} \big [ {\rm Pr}( m_i) \cdot p_{\boldsymbol{ r}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m } \hspace{0.05cm} (\boldsymbol{\rho}_1 \hspace{0.05cm} \vert \hspace{0.05cm}m_i ) \big]\hspace{0.05cm}.$$
$\text{Beispiel 2:}$ Wir betrachten zur Verdeutlichung des soeben vorgestellten Theorems der Irrelevanz zwei verschiedene Systemkonfigurationen mit jeweils zwei Rauschtermen $\boldsymbol{ n}_1$ und $\boldsymbol{ n}_2$.
- In der Grafik sind alle vektoriellen Größen rot eingezeichnet.
- Die Größen $\boldsymbol{s}$, $\boldsymbol{ n}_1$ und $\boldsymbol{ n}_2$ seien zudem jeweils unabhängig voneinander.
Die Analyse dieser beiden Anordnungen liefert folgende Ergebnisse:
- Der Entscheider muss in beiden Fällen die Komponente $\boldsymbol{ r}_1= \boldsymbol{ s}_i + \boldsymbol{ n}_1$ berücksichtigen, da nur diese die Information über das Nutzsignal $\boldsymbol{ s}_i$ und damit über die gesendete Nachricht $m_i$ liefert.
- Bei der oberen Konfiguration enthält $\boldsymbol{ r}_2$ keine Information über $m_i$, die nicht bereits von $\boldsymbol{ r}_1$ geliefert wurde. Vielmehr ist $\boldsymbol{ r}_2= \boldsymbol{ r}_1 + \boldsymbol{ n}_2$ nur eine verrauschte Version von $\boldsymbol{ r}_1$ und hängt nur vom Rauschen $\boldsymbol{ n}_2$ ab, sobald $\boldsymbol{ r}_1$ bekannt ist ⇒ $\boldsymbol{ r}_2$ ist irrelevant:
- $$p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i )= p_{\boldsymbol{ r}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm}\boldsymbol{\rho}_1 )= p_{\boldsymbol{ n}_2 } \hspace{0.05cm} (\boldsymbol{\rho}_2 - \boldsymbol{\rho}_1 )\hspace{0.05cm}.$$
- Bei der unteren Konfiguration ist $\boldsymbol{ r}_2= \boldsymbol{ n}_1 + \boldsymbol{ n}_2$ für den Empfänger hilfreich, da ihm so ein Schätzwert für den Rauschterm $\boldsymbol{ n}_1$ geliefert wird ⇒ $\boldsymbol{ r}_2$ sollte hier nicht verworfen werden. Formal lässt sich dieses Resultat wie folgt ausdrücken:
- $$p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ r}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2\hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 , \hspace{0.05cm}m_i ) = p_{\boldsymbol{ r}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ n}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 - \boldsymbol{s}_i, \hspace{0.05cm}m_i)= p_{\boldsymbol{ n}_2 \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{ n}_1 , \hspace{0.05cm} m } \hspace{0.05cm} (\boldsymbol{\rho}_2- \boldsymbol{\rho}_1 + \boldsymbol{s}_i \hspace{0.05cm} \vert \hspace{0.05cm} \boldsymbol{\rho}_1 - \boldsymbol{s}_i, \hspace{0.05cm}m_i) = p_{\boldsymbol{ n}_2 } \hspace{0.05cm} (\boldsymbol{\rho}_2- \boldsymbol{\rho}_1 + \boldsymbol{s}_i ) \hspace{0.05cm}.$$
- Da nun im Argument dieser Funktion die Nachricht $\boldsymbol{ s}_i$ erscheint, ist $\boldsymbol{ r}_2$ nicht irrelevant”, sondern durchaus relevant.
Einige Eigenschaften des AWGN-Kanals
Um weitere Aussagen über die Art der optimalen Messungen des Vektors $\boldsymbol{ r}$ machen zu können, ist es notwendig, die den Kanal charakterisierende (bedingte) Wahrscheinlichkeitsdichtefunktion $p_{\hspace{0.02cm}r(t)\hspace{0.05cm} \vert \hspace{0.05cm}s(t)}$ weiter zu spezifizieren. Im Folgenden wird die Kommunikation über den "AWGN–Kanal" betrachtet, dessen wichtigste Eigenschaften hier nochmals kurz zusammengestellt werden:
- Das Ausgangssignal des AWGN–Kanals ist $r(t) = s(t)+n(t)$, wobei $s(t)$ das Sendesignal angibt und $n(t)$ durch einen Gaußschen Rauschprozess dargestellt wird.
- Einen Zufallsprozess $\{n(t)\}$ bezeichnet man als gaußisch, wenn die Elemente der $k$–dimensionalen Zufallsvariablen $\{n_1(t)\hspace{0.05cm} \text{...} \hspace{0.05cm}n_k(t)\}$ gemeinsam gaußverteilt sind ⇒ "jointly Gaussian".
- Der Mittelwert des AWGN–Rauschens ist ${\rm E}\big[n(t)\big] = 0$. Außerdem ist $n(t)$ "weiß", was bedeutet, dass das Leistungsdichtespektrum $\rm (LDS)$ für alle Frequenzen $($von $-\infty$ bis $+\infty)$ konstant ist:
- $${\it \Phi}_n(f) = {N_0}/{2} \hspace{0.05cm}.$$
- Nach dem Wiener–Chintchine–Theorem ergibt sich die Autokorrelationsfunktion $\rm (AKF)$ als die Fourierrücktransformierte von ${\it \Phi_n(f)}$:
- $${\varphi_n(\tau)} = {\rm E}\big [n(t) \cdot n(t+\tau)\big ] = {N_0}/{2} \cdot \delta(t)\hspace{0.3cm} \Rightarrow \hspace{0.3cm} {\rm E}\big [n(t) \cdot n(t+\tau)\big ] = \left\{ \begin{array}{c} \rightarrow \infty \\ 0 \end{array} \right.\quad \begin{array}{*{1}c} {\rm f{\rm \ddot{u}r}} \hspace{0.25cm} \tau = 0 \hspace{0.05cm}, \\ {\rm f{\rm \ddot{u}r}} \hspace{0.25cm} \tau \ne 0 \hspace{0.05cm},\\ \end{array}$$
- $N_0$ gibt dabei die physikalische $($nur für $f \ge 0$ definierte$)$ Rauschleistungsdichte an. Der konstante LDS–Wert $(N_0/2)$ und das Gewicht der Diracfunktion in der AKF $($ebenfalls $N_0/2)$ ergibt sich allein durch die zweiseitige Betrachtungsweise.
⇒ Weitere Informationen zu diesem Thema liefert das Lernvideo "Der AWGN-Kanal" im zweiten Teil.
Beschreibung des AWGN-Kanals durch orthonormale Basisfunktionen
Aus dem vorletzten Statement auf der letzten Seite geht hervor, dass
- reines AWGN–Rauschen $n(t)$ stets eine unendliche Varianz (Leistung) aufweist: $\sigma_n^2 \to \infty$,
- in der Realität demzufolge nur gefiltertes Rauschen $n\hspace{0.05cm}'(t) = n(t) \star h_n(t)$ auftreten kann.
Mit der Impulsantwort $h_n(t)$ und dem Frequenzgang $H_n(f) = {\rm F}\big [h_n(t)\big ]$ gelten dann folgende Gleichungen:
- $${\rm E}\big[n\hspace{0.05cm}'(t) \big] \hspace{0.15cm} = \hspace{0.2cm} {\rm E}\big[n(t) \big] = 0 \hspace{0.05cm},$$
- $${\it \Phi_{n\hspace{0.05cm}'}(f)} \hspace{0.1cm} = \hspace{0.1cm} {N_0}/{2} \cdot |H_{n}(f)|^2 \hspace{0.05cm},$$
- $$ {\it \varphi_{n\hspace{0.05cm}'}(\tau)} \hspace{0.1cm} = \hspace{0.1cm} {N_0}/{2}\hspace{0.1cm} \cdot \big [h_{n}(\tau) \star h_{n}(-\tau)\big ]\hspace{0.05cm},$$
- $$\sigma_n^2 \hspace{0.1cm} = \hspace{0.1cm} { \varphi_{n\hspace{0.05cm}'}(\tau = 0)} = {N_0}/{2} \cdot \int_{-\infty}^{+\infty}h_n^2(t)\,{\rm d} t ={N_0}/{2}\hspace{0.1cm} \cdot < \hspace{-0.1cm}h_n(t), \hspace{0.1cm} h_n(t) \hspace{-0.05cm} > \hspace{0.1cm} = \int_{-\infty}^{+\infty}{\it \Phi}_{n\hspace{0.05cm}'}(f)\,{\rm d} f = {N_0}/{2} \cdot \int_{-\infty}^{+\infty}|H_n(f)|^2\,{\rm d} f \hspace{0.05cm}.$$
Im Folgenden beinhaltet $n(t)$ stets implizit eine Bandbegrenzung; auf die Schreibweise $n'(t)$ wird also zukünftig verzichtet.
$\text{Beachten Sie:}$ Ähnlich wie das Sendesignal $s(t)$ lässt sich auch der Rauschprozess $\{n(t)\}$ als gewichtete Summe orthonormaler Basisfunktionen $\varphi_j(t)$ schreiben.
- Im Gegensatz zu $s(t)$ ist nun allerdings eine Beschränkung auf eine endliche Anzahl an Basisfunktionen nicht möglich.
- Vielmehr gilt bei rein stochastischen Größen für die entsprechende Signaldarstellung stets
- $$n(t) = \lim_{N \rightarrow \infty} \sum\limits_{j = 1}^{N}n_j \cdot \varphi_j(t) \hspace{0.05cm},$$
- wobei der Koeffizient $n_j$ durch die Projektion von $n(t)$ auf die Basisfunktion $\varphi_j(t)$ bestimmt ist:
- $$n_j = \hspace{0.1cm} < \hspace{-0.1cm}n(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.05cm}.$$
Hinweis: Um eine Verwechslung mit den Basisfunktionen $\varphi_j(t)$ zu vermeiden, wird im Folgenden die AKF $\varphi_n(\tau)$ des Rauschprozesses stets nur noch als der Erwartungswert ausgedrückt:
- $${\rm E}\big [n(t) \cdot n(t + \tau)\big ].$$
Optimaler Empfänger für den AWGN-Kanal
Auch das Empfangssignal $r(t) = s(t) + n(t)$ lässt sich in bekannter Weise in Basisfunktionen zerlegen:
- $$r(t) = \sum\limits_{j = 1}^{\infty}r_j \cdot \varphi_j(t) \hspace{0.05cm}.$$
Zu berücksichtigen ist:
- Die $M$ möglichen Sendesignale $\{s_i(t)\}$ spannen einen Signalraum mit insgesamt $N$ Basisfunktionen $\varphi_1(t)$, ... , $\varphi_N(t)$ auf.
- Diese $N$ Basisfunktionen $\varphi_j(t)$ werden gleichzeitig zur Beschreibung des Rauschsignals $n(t)$ und des Empfangssignals $r(t)$ verwendet.
- Zur vollständigen Charakterisierung von $n(t)$ bzw. $r(t)$ werden nun aber auch noch unendlich viele weitere Basisfunktionen $\varphi_{N+1}(t)$, $\varphi_{N+2}(t)$, ... benötigt.
Damit ergeben sich die Koeffizienten des Empfangssignals $r(t)$ gemäß folgender Gleichung. Hierbei ist berücksichtigt, dass die (möglichen) Sendesignale $s_i(t)$ und das Rauschen $n(t)$ voneinander unabhängig sind:
- $$r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.1cm}=\hspace{0.1cm} \left\{ \begin{array}{c} < \hspace{-0.1cm}s_i(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > + < \hspace{-0.1cm}n(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.1cm}= s_{ij}+ n_j\\ < \hspace{-0.1cm}n(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} > \hspace{0.1cm} = n_j \end{array} \right.\quad \begin{array}{*{1}c} {j = 1, 2, \hspace{0.05cm}\text{...}\hspace{0.05cm} \hspace{0.05cm}, N} \hspace{0.05cm}, \\ {j > N} \hspace{0.05cm}.\\ \end{array}$$
Somit ergibt sich für den optimalen Empfänger die oben skizzierte Struktur.
Betrachten wir zunächst den AWGN–Kanal. Hier kann auf das Vorfilter mit dem Frequenzgang $W(f)$ verzichtet werden, das für farbiges Rauschen vorgesehen ist.
- Der Detektor des optimalen Empfängers bildet die Koeffizienten $r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t)\hspace{-0.05cm} >$ und reicht diese an den Entscheider weiter.
- Basiert die Entscheidung auf sämtlichen Koeffizienten $r_j$, so ist die Wahrscheinlichkeit für eine Fehlentscheidung minimal und der Empfänger optimal.
- Die reellwertigen Koeffizienten $r_j$ wurden oben wie folgt berechnet:
- $$r_j = \left\{ \begin{array}{c} s_{ij} + n_j\\ n_j \end{array} \right.\quad \begin{array}{*{1}c} {j = 1, 2, \hspace{0.05cm}\text{...}\hspace{0.05cm}, N} \hspace{0.05cm}, \\ {j > N} \hspace{0.05cm}.\\ \end{array}$$
Nach dem "Theorem der Irrelevanz" lässt sich zeigen, dass für additives weißes Gaußsches Rauschen
- die Optimalität nicht herabgesetzt wird, wenn man die Koeffizienten $r_{N+1}$, $r_{N+2}$, ... , die nicht von der Nachricht $(s_{ij})$ abhängigen, nicht in den Entscheidungsprozess einbindet, und somit
- der Detektor nur die Projektionen des Empfangssignals $r(t)$ auf die $N$ durch das Nutzsignal $s(t)$ vorgegebenen Basisfunktionen $\varphi_{1}(t)$, ... , $\varphi_{N}(t)$ bilden muss.
In der Grafik ist diese signifikante Vereinfachung durch die graue Hinterlegung angedeutet.
Bei farbigem Rauschen ⇒ Leistungsdichtespektrum ${\it \Phi}_n(f) \ne {\rm const.}$ ist lediglich zusätzlich ein Vorfilter mit Amplitudengang $|W(f)| = {1}/{\sqrt{\it \Phi}_n(f)}$ erforderlich.
- Man nennt dieses Filter auch "Whitening Filter", da die Rauschleistungsdichte am Ausgang wieder konstant – also „weiß” – ist.
- Genaueres hierzu finden Sie im Kapitel "Matched-Filter bei farbigen Störungen" des Buches „Stochastische Signaltheorie”.
Implementierungsaspekte
Wesentliche Bestandteile des optimalen Empfängers sind die Berechnungen der inneren Produkte gemäß den Gleichungen $r_j \hspace{0.1cm} = \hspace{0.1cm} \hspace{0.1cm} < \hspace{-0.1cm}r(t), \hspace{0.1cm} \varphi_j(t) \hspace{-0.05cm} >$.
$\text{Diese können auf verschiedene Art und Weise implementiert werden:}$
- Beim Korrelationsempfänger $($Näheres zu dieser Implementierung finden Sie im gleichnamigen Kapitel$)$ werden die inneren Produkte direkt entsprechend der Definition mit analogen Multiplizierern und Integratoren realisiert:
- $$r_j = \int_{-\infty}^{+\infty}r(t) \cdot \varphi_j(t) \,{\rm d} t \hspace{0.05cm}.$$
- Der Matched–Filter–Empfänger, der bereits im Kapitel "Optimaler Binärempfänger" zu Beginn dieses Buches hergeleitet wurde, erzielt mit einem linearen Filter mit der Impulsantwort $h_j(t) = \varphi_j(t) \cdot (T-t)$ und anschließender Abtastung zum Zeitpunkt $t = T$ das gleiche Ergebnis:
- $$r_j = \int_{-\infty}^{+\infty}r(\tau) \cdot h_j(t-\tau) \,{\rm d} \tau = \int_{-\infty}^{+\infty}r(\tau) \cdot \varphi_j(T-t+\tau) \,{\rm d} \tau \hspace{0.3cm} \Rightarrow \hspace{0.3cm} r_j (t = \tau) = \int_{-\infty}^{+\infty}r(\tau) \cdot \varphi_j(\tau) \,{\rm d} \tau = r_j \hspace{0.05cm}.$$
Die Abbildung zeigt die beiden möglichen Realisierungsformen des optimalen Detektors.
Wahrscheinlichkeitsdichtefunktion der Empfangswerte
Bevor wir uns im folgenden Kapitel der optimalen Gestaltung des Entscheiders und der Berechnung und Annäherung der Fehlerwahrscheinlichkeit zuwenden, erfolgt zunächst eine für den AWGN–Kanal gültige statistische Analyse der Entscheidungsgrößen $r_j$.
Dazu betrachten wir nochmals den optimalen Binärempfänger für die bipolare Basisbandübertragung über den AWGN–Kanal, wobei wir von der für das vierte Hauptkapitel gültigen Beschreibungsform ausgehen.
Mit den Parametern $N = 1$ und $M = 2$ ergibt sich für das Sendesignal die in der linken Grafik dargestellte Signalraumkonstellation
- mit nur einer Basisfunktion $\varphi_1(t)$, wegen $N = 1$,
- mit den beiden Signalraumpunkten $s_i \in \{s_0, \hspace{0.05cm}s_1\}$, wegen $M = 2$.
Für das Signal $r(t) = s(t) + n(t)$ am Ausgang des AWGN–Kanals ergibt sich im rauschfreien Fall ⇒ $r(t) = s(t)$ die genau gleiche Konstellation. Die Signalraumpunkte liegen somit bei
- $$r_0 = s_0 = \sqrt{E}\hspace{0.05cm},\hspace{0.2cm}r_1 = s_1 = -\sqrt{E}\hspace{0.05cm}.$$
- Bei Berücksichtigung des (bandbegrenzten) AWGN–Rauschens $n(t)$ überlagern sich den beiden Punkten $r_0$ und $r_1$ jeweils Gaußkurven mit der Varianz $\sigma_n^2$ ⇒ Streuung $\sigma_n$ $($siehe rechte Grafik$)$.
- Die WDF der Rauschkomponente $n(t)$ lautet dabei:
- $$p_n(n) = \frac{1}{\sqrt{2\pi} \cdot \sigma_n}\cdot {\rm e}^{ - {n^2}/(2 \sigma_n^2)}\hspace{0.05cm}.$$
Für die bedingte Wahrscheinlichkeitsdichte, dass der Empfangswert $\rho$ anliegt, wenn $s_i$ gesendet wurde, ergibt sich dann folgender Ausdruck:
- $$p_{\hspace{0.02cm}r\hspace{0.05cm}|\hspace{0.05cm}s}(\rho\hspace{0.05cm}|\hspace{0.05cm}s_i) = \frac{1}{\sqrt{2\pi} \cdot \sigma_n}\cdot {\rm e}^{ - {(\rho - s_i)^2}/(2 \sigma_n^2)} \hspace{0.05cm}.$$
Zu den Einheiten der hier aufgeführten Größen ist zu bemerken:
- $r_0 = s_0$ und $r_1 = s_1$ sowie $n$ sind jeweils Skalare mit der Einheit „Wurzel aus Energie”.
- Damit ist offensichtlich, dass $\sigma_n$ ebenfalls die Einheit „Wurzel aus Energie” besitzt und $\sigma_n^2$ eine Energie darstellt.
- Beim AWGN–Kanal ist die Rauschvarianz $\sigma_n^2 = N_0/2$. Diese ist also ebenfalls eine physikalische Größe mit der Einheit $\rm W/Hz = Ws$.
Die hier angesprochene Thematik wird in der "Aufgabe 4.6" an Beispielen verdeutlicht.
N–dimensionales Gaußsches Rauschen
Liegt ein $N$–dimensionales Modulationsverfahren vor, das heißt, wenn mit $0 \le i \le M-1$ und $1 \le j \le N$
- $$s_i(t) = \sum\limits_{j = 1}^{N} s_{ij} \cdot \varphi_j(t) = s_{i1} \cdot \varphi_1(t) + s_{i2} \cdot \varphi_2(t) + \hspace{0.05cm}\text{...}\hspace{0.05cm} + s_{iN} \cdot \varphi_N(t)\hspace{0.05cm}\hspace{0.3cm} \Rightarrow \hspace{0.3cm} \boldsymbol{ s}_i = \left(s_{i1}, s_{i2}, \hspace{0.05cm}\text{...}\hspace{0.05cm}, s_{iN}\right ) \hspace{0.05cm},$$
gilt, dann muss der Rauschvektor $\boldsymbol{ n}$ ebenfalls mit Dimension $N$ angesetzt werden. Das gleiche gilt auch für den Empfangsvektor $\boldsymbol{ r}$:
- $$\boldsymbol{ n} = \left(n_{1}, n_{2}, \hspace{0.05cm}\text{...}\hspace{0.05cm}, n_{N}\right ) \hspace{0.01cm},$$
- $$\boldsymbol{ r} = \left(r_{1}, r_{2}, \hspace{0.05cm}\text{...}\hspace{0.05cm}, r_{N}\right )\hspace{0.05cm}.$$
Die Wahrscheinlichkeitsdichtefunktion $\rm (WDF)$ lautet dann für den AWGN–Kanal mit der Realisierung $\boldsymbol{ \eta}$ des Rauschsignals
- $$p_{\boldsymbol{ n}}(\boldsymbol{ \eta}) = \frac{1}{\left( \sqrt{2\pi} \cdot \sigma_n \right)^N } \cdot {\rm exp} \left [ - \frac{|| \boldsymbol{ \eta} ||^2}{2 \sigma_n^2}\right ]\hspace{0.05cm},$$
und für die bedingte WDF in der Maximum–Likelihood–Entscheidungsregel ist anzusetzen:
- $$p_{\hspace{0.02cm}\boldsymbol{ r}\hspace{0.05cm} | \hspace{0.05cm} \boldsymbol{ s}}(\boldsymbol{ \rho} \hspace{0.05cm}|\hspace{0.05cm} \boldsymbol{ s}_i) \hspace{-0.1cm} = \hspace{0.1cm} p_{\hspace{0.02cm} \boldsymbol{ n}\hspace{0.05cm} | \hspace{0.05cm} \boldsymbol{ s}}(\boldsymbol{ \rho} - \boldsymbol{ s}_i \hspace{0.05cm} | \hspace{0.05cm} \boldsymbol{ s}_i) = \frac{1}{\left( \sqrt{2\pi} \cdot \sigma_n \right)^2 } \cdot {\rm exp} \left [ - \frac{|| \boldsymbol{ \rho} - \boldsymbol{ s}_i ||^2}{2 \sigma_n^2}\right ]\hspace{0.05cm}.$$
Die Gleichung ergibt sich
- aus der allgemeinen Darstellung der $N$–dimensionalen Gaußschen WDF im Abschnitt "Korrelationsmatrix" des Buches „Stochastische Signaltheorie”
- unter der Voraussetzung, dass die Komponenten unkorreliert (und somit statistisch unabhängig) sind.
- $||\boldsymbol{ \eta}||$ bezeichnet man als die "Norm" (Länge) des Vektors $\boldsymbol{ \eta}$.
$\text{Beispiel 3:}$ Rechts dargestellt ist die zweidimensionale Gauß–WDF $p_{\boldsymbol{ n} } (\boldsymbol{ \eta})$ der 2D–Zufallsgröße $\boldsymbol{ n} = (n_1,\hspace{0.05cm}n_2)$. Beliebige Realisierungen der Zufallsgröße $\boldsymbol{ n}$ werden mit $\boldsymbol{ \eta} = (\eta_1,\hspace{0.05cm}\eta_2)$ bezeichnet. Die Gleichung der dargestellten Glockenkurve lautet:
- $$p_{n_1, n_2}(\eta_1, \eta_2) = \frac{1}{\left( \sqrt{2\pi} \cdot \sigma_n \right)^2 } \cdot {\rm exp} \left [ - \frac{ \eta_1^2 + \eta_2^2}{2 \sigma_n^2}\right ]\hspace{0.05cm}. $$
- Das Maximum dieser Funktion liegt bei $\eta_1 = \eta_2 = 0$ und hat den Wert $2\pi \cdot \sigma_n^2$.
- Mit $\sigma_n^2 = N_0/2$ lässt sich die 2D–WDF in Vektorform auch wie folgt schreiben:
- $$p_{\boldsymbol{ n} }(\boldsymbol{ \eta}) = \frac{1}{\pi \cdot N_0 } \cdot {\rm exp} \left [ - \frac{\vert \vert \boldsymbol{ \eta} \vert \vert ^2}{N_0}\right ]\hspace{0.05cm}.$$
- Diese rotationssymmetrische WDF eignet sich zum Beispiel für die Beschreibung/Untersuchung eines zweidimensionalen Modulationsverfahrens wie "M–QAM", "M–PSK" oder "2–FSK".
- Oft werden zweidimensionale reelle Zufallsgrößen aber auch eindimensional–komplex dargestellt, meist in der Form $n(t) = n_{\rm I}(t) + {\rm j} \cdot n_{\rm Q}(t)$. Die beiden Komponenten bezeichnet man dann als "Inphaseanteil" $n_{\rm I}(t)$ und "Quadraturanteil" $n_{\rm Q}(t)$ des Rauschens.
- Die Wahrscheinlichkeitsdichtefunktion hängt nur vom Betrag $\vert n(t) \vert$ der Rauschvariablen ab und nicht von Winkel ${\rm arc} \ n(t)$. Das heißt: Komplexes Rauschen ist zirkulär symmetrisch (siehe Grafik).
- Zirkulär symmetrisch bedeutet auch, dass die Inphasekomponente $n_{\rm I}(t)$ und die Quadraturkomponente $n_{\rm Q}(t)$ die gleiche Verteilung aufweisen und damit auch gleiche Varianz (Streuung) besitzen:
- $$ {\rm E} \big [ n_{\rm I}^2(t)\big ] = {\rm E}\big [ n_{\rm Q}^2(t) \big ] = \sigma_n^2 \hspace{0.05cm},\hspace{1cm}{\rm E}\big [ n(t) \cdot n^*(t) \big ]\hspace{0.1cm} = \hspace{0.1cm} {\rm E}\big [ n_{\rm I}^2(t) \big ] + {\rm E}\big [ n_{\rm Q}^2(t)\big ] = 2\sigma_n^2 \hspace{0.05cm}.$$
Abschließend noch einige Bezeichnungsvarianten für Gaußsche Zufallsgrößen:
- $$x ={\cal N}(\mu, \sigma^2) \hspace{-0.1cm}: \hspace{0.3cm}\text{reelle gaußverteilte Zufallsgröße, mit Mittelwert}\hspace{0.15cm}\mu \text { und Varianz}\hspace{0.15cm}\sigma^2 \hspace{0.05cm},$$
- $$y={\cal CN}(\mu, \sigma^2)\hspace{-0.1cm}: \hspace{0.12cm}\text{komplexe gaußverteilte Zufallsgröße} \hspace{0.05cm}.$$
Aufgaben zum Kapitel
Aufgabe 4.4: Maximum–a–posteriori und Maximum–Likelihood
Aufgabe 4.5: Theorem der Irrelevanz