Exercise 1.3Z: Threshold Optimization

Zur Optimierung des Schwellenwertes

In dieser Aufgabe wird ein bipolares Binärsystem mit AWGN–Rauschen (Additive White Gaussian Noise) betrachtet, so dass für die Bitfehlerwahrscheinlichkeit gilt:

$$p_{\rm B} = {\rm Q} \left( \frac{s_0}{\sigma_d}\right)= {1}/{2} \cdot {\rm erfc} \left( \frac{s_0}{\sqrt{2} \cdot \sigma_d}\right) \hspace{0.05cm}.$$

Hierbei sind folgende Funktionen verwendet:

$$\rm Q (\it x) = \frac{\rm 1}{\sqrt{\rm 2\pi}}\int_{\it x}^{+\infty}\rm e^{\it -u^{\rm 2}/\rm 2}\,d {\it u} \hspace{0.05cm},\hspace{1cm}{\rm erfc} ( x) = \frac{\rm 2}{\sqrt{\rm \pi}}\int_{\it x}^{+\infty}\rm e^{\it -u^{\rm 2}}\,d \it u \hspace{0.05cm}.$$

Die obige Gleichung gilt für den Schwellenwert $E = 0$ unabhängig von den Symbolwahrscheinlichkeiten $p_{\rm L}$ und $p_{\rm H}$.
Allerdings kann mit einem anderen Schwellenwert $E$ eine kleinere Fehlerwahrscheinlichkeit erzielt werden, wenn die beiden Auftrittswahrscheinlichkeiten unterschiedlich sind: $(p_{\rm L} ≠ p_{\rm H})$.

Die Streuung des Rauschanteils ist stets $σ_{d} = 0.5 \ \rm V$. Die beiden Amplituden des Detektionsnutzanteils sind mit $±1 \ \rm V$ fest vorgegeben. Zu untersuchen sind folgende Symbolwahrscheinlichkeiten:

$p_{\rm L} = 0.88$ ⇒ $p_{\rm H} = 1- p_{\rm L} =0.12$,
$p_{\rm L} = 0.31$ ⇒ $p_{\rm H} = 1- p_{\rm L} =0.69$.

In der Grafik ist dieser letzte Parametersatz und der Schwellenwert $E = 0.1 \cdot s_{\rm 0}$ berücksichtigt. Dargestellt ist die Wahrscheinlichkeitsdichtefunktion der Detektionsabtastwerte $d$..

Hinweise:

Die Aufgabe gehört zum Kapitel Fehlerwahrscheinlichkeit bei Basisbandübertragung.
Sollte die Eingabe des Zahlenwertes „0” erforderlich sein, so geben Sie bitte „0.” ein.
Zur Bestimmung von Fehlerwahrscheinlichkeiten können Sie das Interaktionsmodul Komplementäre Gaußsche Fehlerfunktionen verwenden.
Für die Ableitung der Q–Funktion gilt:

$$\frac{{\rm d\hspace{0.05cm}Q} ({\it x})} {{\rm d}\hspace{0.05cm}x} = \frac{\rm 1}{\sqrt{\rm 2\pi}} \cdot \rm e^{\it -x^{\rm 2}/\rm 2} \hspace{0.05cm}.$$

Fragebogen

	${\rm erfc}(x) =2 \cdot {\rm Q}(\sqrt{2} \cdot x)$,
	${\rm erfc}(x) =\sqrt{2} \cdot {\rm Q}(x/\sqrt{2})$,
	${\rm erfc}(x) \approx {\rm Q}(x)$.

$p_{\rm B} \ =\ $

$\ \%$

$p_{\rm B} \ =\ $

$\ \%$

$E_{\rm opt} \ =\ $

$\ \rm V$

$p_{\rm B,\ min} \ =\ $

$\ \%$

$E_{\rm opt} \ =\ $

$\ \rm V$

$p_{\rm B,\ min} \ =\ $

$\ \%$

Musterlösung

Solution

(1) Aus der Q–Funktion ergibt sich mit der Substitution $t^{2} = u^{2}/2$:

$$\rm Q (\it x)= \frac{\rm 1}{\sqrt{\rm 2\pi}}\int_{\it x}^{+\infty}\rm e^{\it -u^{\rm 2}/\rm 2}\,d \it u = \frac{\rm 1}{\sqrt{\rm \pi}}\int_{\it x/\sqrt{\rm 2}}^{+\infty}\rm e^{\it -t^{\hspace{0.05cm}\rm 2}}\,{d \it t} = {1}/{2} \cdot {\rm erfc} \left( {x}/{\sqrt{2}} \right)\hspace{0.05cm}.$$

Daraus folgt die Richtigkeit des ersten Lösungsvorschlags:

$${\rm erfc} (x)= 2 \cdot {\rm Q} (\sqrt{2} \cdot x) \hspace{0.05cm}.$$

(2) Unabhängig von den Symbolwahrscheinlichkeiten erhält man mit der Entscheiderschwelle $E = 0$:

$$p_{\rm B} = {\rm Q} \left( {s_0}/{\sigma_d}\right)= {\rm Q} (2) \hspace{0.1cm}\underline {= 2.27 \, \% }\hspace{0.05cm}.$$

(3) Nun lautet die allgemeine Gleichung für die Bitfehlerwahrscheinlichkeit, wobei $d_{\rm N}$ den Rauschanteil von $d(t)$ bezeichnet:

$$p_{\rm B} \ = \ p_{\rm L} \cdot {\rm Pr}( d_{\rm N}> E + s_0)+ p_{\rm H} \cdot {\rm Pr}( d_{\rm N}< E - s_0) \ = \ p_{\rm L} \cdot {\rm Q} \left( \frac{s_0 + E}{\sigma_d}\right)+ p_{\rm H} \cdot {\rm Q} \left( \frac{s_0 - E}{\sigma_d}\right)\hspace{0.05cm}.$$

Hierbei ist die WDF-Symmetrie berücksichtigt. Mit $p_{\rm L} = 0.88$ ⇒ $p_{\rm H} = 0.12$ und $E = 0.1 \hspace{0.05cm} \rm V $erhält man:

$$p_{\rm B} \hspace{-0.05cm}=\hspace{-0.05cm} 0.88 \hspace{-0.03cm}\cdot \hspace{-0.03cm}{\rm Q} \left( \frac{1\, {\rm V} + 0.1\, {\rm V}}{0.5\, {\rm V}}\right)\hspace{-0.05cm}+\hspace{-0.05cm} 0.12 \hspace{-0.03cm}\cdot \hspace{-0.03cm}{\rm Q} \left( \frac{1\, {\rm V} - 0.1\, {\rm V}}{0.5\, {\rm V}}\right) \hspace{-0.05cm}=\hspace{-0.05cm} 0.88 \hspace{-0.03cm}\cdot \hspace{-0.03cm} {\rm Q}(2.2) \hspace{-0.05cm}+\hspace{-0.05cm} 0.12 \cdot {\rm Q}(1.8) \hspace{-0.05cm}= \hspace{-0.05cm}0.88 \hspace{-0.03cm}\cdot \hspace{-0.03cm} 1.39\,\% \hspace{-0.05cm}+\hspace{-0.05cm} 0.12 \hspace{-0.03cm}\cdot \hspace{-0.03cm} 3.59\,\%\ \hspace{0.1cm}\underline { = 1.65\,\% } \hspace{0.05cm}.$$

Durch die Schwellenverschiebung nach rechtsum $E = s_{0}/10$ ergibt sich eine Verbesserung von $p_{\rm B} = 2.27 \ \%$ auf $p_{\rm B} = 1.65 \ \%$.

(4) Diese Optimierungsaufgabe wird durch Nullsetzen der Ableitung gelöst, wobei der Hinweis auf der Angabenseite zu berücksichtigen ist:

$$\frac{{\rm d\hspace{0.05cm}} p_{\rm B}(E)} {{\rm d}\hspace{0.05cm}E} = - \frac{\rm p_{\rm L}}{\sqrt{\rm 2\pi}\cdot {\sigma_d}} \cdot {\rm exp}\left(- \frac{(s_0 + E)^2}{{\rm 2}\cdot {\sigma_d}^2}\right)+\frac{\rm p_{\rm H}}{\sqrt{\rm 2\pi}\cdot {\sigma_d}} \cdot {\rm exp}\left(- \frac{(s_0 - E)^2}{{\rm 2}\cdot {\sigma_d}^2}\right) = 0$$

$$\Rightarrow \hspace{0.3cm} \frac{p_{\rm L}} {p_{\rm H}} = - \frac{ {\rm exp} \left(-\frac{(s_0 - E)^2}{{\rm 2}\cdot {\sigma_d}^2}\right)}{{\rm exp} \left(-\frac{(s_0 + E)^2}{{\rm 2}\cdot {\sigma_d}^2}\right)} = {\rm exp} \left( \frac{2 \cdot E \cdot s_0 }{ {\sigma_d}^2}\right)\hspace{0.05cm}.$$

Damit erhält man für den optimalen Schwellenwert allgemein:

$$E_{\rm opt} = \frac{\sigma_d^2} {2 \cdot s_0} \cdot {\rm ln}\hspace{0.05cm}\frac{p_{\rm L}} {p_{\rm H}} \hspace{0.05cm}.$$

Mit $σ_{d} = 0.5 \ \rm V$, $s_{\rm 0} = 1 \ \rm V$, $p_{\rm L} = 0.88$ und $p_{\rm H} = 0.12$ ergibt sich folgendes Optimum:

$$E_{\rm opt} = \frac{(0.5\, {\rm V})^2} {2 \cdot 1\, {\rm V}} \cdot {\rm ln}\hspace{0.05cm}\frac{0.88} {0.12} \hspace{0.1cm}\underline { \approx 0.25\, {\rm V}}\hspace{0.05cm}.$$

(5) Die minimale Fehlerwahrscheinlichkeit für den optimalen Schwellenwert $E_{\rm opt} = 0.25 \hspace{0.05cm} \rm V$ ist somit:

$$p_{\rm B, \hspace{0.05cm}min} = 0.88 \cdot {\rm Q}(2.5) + 0.12 \cdot {\rm Q}(1.5) = 0.88 \cdot 0.62\,\% + 0.12 \cdot 6.68\,\%\ \hspace{0.1cm}\underline {= 1.35\,\% }\hspace{0.05cm}.$$

Gegenüber $E = 0$ ist die Fehlerwahrscheinlichkeit nun um ca. $40 \%$ kleiner.

(6) Mit dem Ergebnis aus (4) ergibt sich nun für den optimalen Schwellenwert:

$$E_{\rm opt} = \frac{(0.5\, {\rm V})^2} {2 \cdot 1\, {\rm V}} \cdot {\rm ln}\hspace{0.05cm}\frac{0.31} {0.69}\hspace{0.1cm}\underline { \approx -0.1\, {\rm V}}\hspace{0.05cm}.$$

Bitte beachten Sie:

Nachdem hier das Symbol L unwahrscheinlicher ist, muss nun die Schwelle nach links – also weg vom wahrscheinlicheren Symbol – verschoben werden.
Die Herleitung des Ergebnisses zu (4) und die Grafik auf der Angabenseite zeigen, dass der optimale Schwellenwert genau an die Stelle zu setzen ist, bei der sich die beiden Gaußfunktionen schneiden.

(7) Mit dem Ergebnis aus (6) gilt schließlich:

$$p_{\rm B, \hspace{0.05cm}min} = 0.31 \cdot {\rm Q}(1.8) + 0.69 \cdot {\rm Q}(2.2) = 0.31 \cdot 3.59\,\% + 0.69 \cdot 1.39\,\%\ \hspace{0.1cm}\underline { = 2.07\,\% } \hspace{0.05cm}.$$

Aufgrund der weniger gravierenden Unsymmetrie ist die erreichbare Verbesserung mit $9 \%$ geringer als unter Punkt (5) berechnet