Aufgaben:Aufgabe 1.3Z: Schwellenwertoptimierung: Unterschied zwischen den Versionen

Aus LNTwww
Wechseln zu:Navigation, Suche
K (Textersetzung - „Komplementäre_Gaußsche_Fehlerfunktionen_(neues_Applet)“ durch „Komplementäre_Gaußsche_Fehlerfunktionen“)
(7 dazwischenliegende Versionen von 2 Benutzern werden nicht angezeigt)
Zeile 5: Zeile 5:
  
 
[[Datei:P_ID1268__Dig_Z_1_3.png|right|frame|Zur Optimierung des Schwellenwertes]]
 
[[Datei:P_ID1268__Dig_Z_1_3.png|right|frame|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
+
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)= \frac{1}{2} \cdot {\rm erfc} \left( \frac{s_0}{\sqrt{2} \cdot \sigma_d}\right) \hspace{0.05cm}$$
+
:$$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}.$$
gilt. Hierbei sind folgende Funktionen verwendet:
+
Hierbei sind folgende Funktionen verwendet:
 
:$$\rm Q (\it x)  =  \frac{\rm 1}{\sqrt{\rm 2\pi}}\int_{\it
 
:$$\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
+
x}^{+\infty}\rm e^{\it -u^{\rm 2}/\rm 2}\,d {\it u}
\hspace{0.05cm},$$
+
\hspace{0.05cm},\hspace{1cm}{\rm erfc} ({\it x})  =  \frac{\rm 2}{\sqrt{\rm
:$${\rm erfc} (\it x)  =  \frac{\rm 2}{\sqrt{\rm
 
 
\pi}}\int_{\it x}^{+\infty}\rm e^{\it -u^{\rm 2}}\,d \it u
 
\pi}}\int_{\it x}^{+\infty}\rm e^{\it -u^{\rm 2}}\,d \it u
 
\hspace{0.05cm}.$$
 
\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}$.  
+
*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}$ )
+
*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 $σ_{\rm d} = 0.5 \ \rm V$, die beiden Amplituden des Detektionsnutzanteils sind mit $±1 V$ fest vorgegeben. Zu untersuchen sind folgende Symbolwahrscheinlichkeiten:
+
 
*$p_{\rm L} = 0.88$ und $p_{\rm H} = 0.12$,
+
 
*$p_{\rm L} = 0.31$ und $p_{\rm H} = 0.69.$
+
Die Streuung des Rauschanteils ist stets  $σ_{d} = 0.5 \ \rm V$. Die beiden Amplituden des Detektionsnutzanteils sind mit  $±1 \ \rm V$  fest vorgegeben.  
In der Grafik ist dieser letzte Parametersatz und der Schwellenwert $E = 0.1 \cdot s_{\rm 0}$ dargestellt.
+
 
 +
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:''  
 
''Hinweise:''  
 
+
*Die Aufgabe gehört zum  Kapitel   [[Digitalsignalübertragung/Fehlerwahrscheinlichkeit_bei_Basisbandübertragung| Fehlerwahrscheinlichkeit bei Basisbandübertragung]].
*Die Aufgabe bezieht sich auf das [[Digitalsignalübertragung/Fehlerwahrscheinlichkeit_bei_Basisbandübertragung| Fehlerwahrscheinlichkeit bei Basisbandübertragung.]] 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}}
+
*Zur Bestimmung von Fehlerwahrscheinlichkeiten können Sie das interaktive Applet  [[Applets:Komplementäre_Gaußsche_Fehlerfunktionen|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}.$$
 
\cdot \rm e^{\it -x^{\rm 2}/\rm 2} \hspace{0.05cm}.$$
*Die Werte der Funktion Q(x) können Sie mit folgendem Interaktionsmodul bestimmen:
 
  
[[Komplementäre Gaußsche Fehlerfunktionen]]
 
  
 
===Fragebogen===
 
===Fragebogen===
  
 
<quiz display=simple>
 
<quiz display=simple>
{Welcher Zusammenhang besteht zwischen Q(x) und erfc(x)?
+
{Welcher Zusammenhang besteht zwischen den Funktionen &nbsp;${\rm Q}(x)$&nbsp; und &nbsp;${\rm erfc}(x)$?
 
|type="[]"}
 
|type="[]"}
+erfc$(x) = $ $2 \cdot$ Q($2^{1/2} \cdot x$),
+
+${\rm erfc}(x) =2 \cdot {\rm Q}(\sqrt{2} \cdot x)$,
-erfc$(x) = $ $2^{1/2} \cdot$ Q($x/2^{1/2}$),
+
-${\rm erfc}(x) =\sqrt{2} \cdot {\rm Q}(x/\sqrt{2})$,
-erfc$(x) \approx$ Q($x$).
+
-${\rm erfc}(x) \approx {\rm Q}(x)$.
  
  
{Welche Fehlerwahrscheinlichkeit ergibt sich mit $p_{\rm L} = 0.88$ und $E = 0$?
+
{Welche Bitfehlerwahrscheinlichkeit ergibt sich mit &nbsp;$\underline{p_{\rm L} = 0.88}$&nbsp; und &nbsp;$\underline{E = 0}$?
 
|type="{}"}
 
|type="{}"}
$E = 0:  p_{\rm B} \ =\ $ { 2,27 3% } $\%$
+
$p_{\rm B} \ =\ $ { 2,27 3% } $\ \%$
  
{Welche Fehlerwahrscheinlichkeit ergibt sich mit $p_{\rm L} = 0.88$ und $E = 0.1 \rm V$?
+
{Welche Bitfehlerwahrscheinlichkeit ergibt sich mit &nbsp;$p_{\rm L} = 0.88$&nbsp; und &nbsp;$\underline{E = 0.1 \hspace{0.05cm} \rm V}$?
 
|type="{}"}
 
|type="{}"}
$E = 0.1 \rm V: $ $p_{\rm B} \ =\ $ { 1,65 3% } $\%$
+
$p_{\rm B} \ =\ $ { 1,65 3% } $\ \%$
  
{Bestimmen Sie den optimalen Schwellenwert für $p_{\rm L} = 0.88$.
+
{Bestimmen Sie den optimalen Schwellenwert &nbsp;$E_{\rm opt}$&nbsp; für &nbsp;$p_{\rm L} = 0.88$ &nbsp; &rArr; &nbsp; $p_{\rm H} = 0.12$.
 
|type="{}"}
 
|type="{}"}
$ p_{\rm L} = 0.88:  E_{\rm opt} \ =\ $ { 0.25 3% } $\ \rm V$
+
$E_{\rm opt} \ =\ $ { 0.25 3% } $\ \rm V$
  
{Wie groß ist die minimale Fehlerwahrscheinlichkeit mit $p_{\rm L} = 0.88$.
+
{Es gelte &nbsp;$E = E_{\rm opt}$. Wie groß ist die minimale Bitfehlerwahrscheinlichkeit mit &nbsp;$p_{\rm L} = 0.88$?
 
|type="{}"}
 
|type="{}"}
$ p_{\rm L} = 0.88:  E_{\rm opt}: p_{\rm B,\ min} \ =\ $ { 1.35 3% } $\%$
+
$p_{\rm B,\ min} \ =\ $ { 1.35 3% } $\ \%$
  
{Bestimmen Sie den optimalen Schwellenwert für $p_{\rm L} = 0.31$.
+
{Bestimmen Sie den optimalen Schwellenwert &nbsp;$E_{\rm opt}$&nbsp; für &nbsp;$\underline{p_{\rm L} = 0.31}$ &nbsp; &rArr; &nbsp; $p_{\rm H} = 0.69$.
 
|type="{}"}
 
|type="{}"}
$ p_{\rm L} = 0.31:  E_{\rm opt} \ =\ $ { -0.103--0.097 } $\ \rm V$
+
$E_{\rm opt} \ =\ $ { -0.103--0.097 } $\ \rm V$
  
{Wie groß ist die minimale Fehlerwahrscheinlichkeit mit $p_{\rm L} = 0.31$.
+
{Wie groß ist die minimale Bitfehlerwahrscheinlichkeit für diesen Fall &nbsp;$(p_{\rm L} = 0.31)$?
 
|type="{}"}
 
|type="{}"}
$ p_{\rm L} = 0.31:  E_{\rm opt}: p_{\rm B,\ min} \ =\ $ { 2.07 3% } $\%$
+
$p_{\rm B,\ min} \ =\ $ { 2.07 3% } $\ \%$
  
 
</quiz>
 
</quiz>
Zeile 72: Zeile 81:
 
x}^{+\infty}\rm e^{\it -u^{\rm 2}/\rm 2}\,d \it u =  
 
x}^{+\infty}\rm e^{\it -u^{\rm 2}/\rm 2}\,d \it u =  
 
\frac{\rm 1}{\sqrt{\rm \pi}}\int_{\it
 
\frac{\rm 1}{\sqrt{\rm \pi}}\int_{\it
x/\sqrt{\rm 2}}^{+\infty}\rm e^{\it -t^{\rm 2}}\,{d \it t} =
+
x/\sqrt{\rm 2}}^{+\infty}\rm e^{\it -t^{\hspace{0.05cm}\rm 2}}\,{d \it t} =
\frac{1}{2} \cdot {\rm erfc} \left( {x}/{\sqrt{2}}
+
{1}/{2} \cdot {\rm erfc} \left( {x}/{\sqrt{2}}
 
\right)\hspace{0.05cm}.$$
 
\right)\hspace{0.05cm}.$$
 
Daraus folgt die Richtigkeit des <u>ersten Lösungsvorschlags</u>:
 
Daraus folgt die Richtigkeit des <u>ersten Lösungsvorschlags</u>:
 
:$${\rm erfc} (x)= 2 \cdot {\rm Q} (\sqrt{2} \cdot x) \hspace{0.05cm}.$$
 
:$${\rm erfc} (x)= 2 \cdot {\rm Q} (\sqrt{2} \cdot x) \hspace{0.05cm}.$$
 +
  
 
'''(2)'''&nbsp; Unabhängig von den Symbolwahrscheinlichkeiten erhält man mit der Entscheiderschwelle $E = 0$:
 
'''(2)'''&nbsp; 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}.$$
 
:$$p_{\rm B} =  {\rm Q} \left( {s_0}/{\sigma_d}\right)= {\rm Q}  (2) \hspace{0.1cm}\underline {= 2.27 \, \% }\hspace{0.05cm}.$$
 +
  
 
'''(3)'''&nbsp; Nun lautet die allgemeine Gleichung für die Bitfehlerwahrscheinlichkeit, wobei $d_{\rm N}$ den Rauschanteil von $d(t)$ bezeichnet:
 
'''(3)'''&nbsp; Nun lautet die allgemeine Gleichung für die Bitfehlerwahrscheinlichkeit, wobei $d_{\rm N}$ den Rauschanteil von $d(t)$ bezeichnet:
Zeile 85: Zeile 96:
 
   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 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}.$$
 
   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 &#8594; p_{\rm H} = 0.12$ und $E = 0.1 \ \rm V $erhält man:
+
*Hierbei ist die WDF-Symmetrie berücksichtigt. Mit $p_{\rm L}  = 0.88$ &nbsp; &rArr; &nbsp; $p_{\rm H} = 0.12$ und $E = 0.1 \hspace{0.05cm} \rm V $erhält man:
:$$p_{\rm B} = 0.88 \cdot {\rm Q} \left( \frac{1\, {\rm V} + 0.1\, {\rm V}}{0.5\, {\rm V}}\right)+
+
:$$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 \cdot {\rm Q} \left( \frac{1\, {\rm V} - 0.1\, {\rm V}}{0.5\, {\rm V}}\right) =\\ 0.88 \cdot {\rm Q}(2.2) + 0.12 \cdot {\rm Q}(1.8) = 0.88 \cdot 1.39\,\% + 0.12 \cdot 3.59\,\%\ \hspace{0.1cm}\underline { = 1.65\,\% } \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 rechts (positiver Wert) um $E = s_{0}/10$ ergibt sich also bereits eine Verbesserung von $p_{\rm B} = 2.27 \ \%$ auf $p_{\rm B} = 1.65  \ \%$.
+
*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)'''&nbsp; Diese Optimierungsaufgabe wird durch Nullsetzen der Ableitung gelöst, wobei der Hinweis auf der Angabenseite zu berücksichtigen ist:
 
'''(4)'''&nbsp; Diese Optimierungsaufgabe wird durch Nullsetzen der Ableitung gelöst, wobei der Hinweis auf der Angabenseite zu berücksichtigen ist:
Zeile 101: Zeile 113:
 
2}\cdot {\sigma_d}^2}\right)} = {\rm exp} \left( \frac{2 \cdot  E
 
2}\cdot {\sigma_d}^2}\right)} = {\rm exp} \left( \frac{2 \cdot  E
 
\cdot s_0 }{ {\sigma_d}^2}\right)\hspace{0.05cm}.$$
 
\cdot s_0 }{ {\sigma_d}^2}\right)\hspace{0.05cm}.$$
Damit erhält man für den optimalen Schwellenwert allgemein:
+
*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}.$$
 
:$$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 $&sigma;_{\rm 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:
+
*Mit $&sigma;_{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}.$$
 
:$$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)'''&nbsp; Die minimale Fehlerwahrscheinlichkeit für den optimalen Schwellenwert $E_{\rm opt} = 0.25 \ \rm V$ ist somit:
+
 
 +
'''(5)'''&nbsp; 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}.$$
 
:$$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.
+
*Gegenüber $E = 0$ ist die Fehlerwahrscheinlichkeit nun um ca. $40 \%$ kleiner.
 +
 
  
'''(6)'''&nbsp; Mit dem Ergebnis aus (4) ergibt sich nun für den optimalen Schwellenwert:  
+
'''(6)'''&nbsp; 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}.$$
 
:$$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}.$$
Nachdem hier das Symbol '''L''' unwahrscheinlicher ist, muss nun die Schwelle nach links &ndash; also weg vom wahrscheinlicheren Symbol &ndash; 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.
+
Bitte beachten Sie:
 +
*Nachdem hier das Symbol '''L''' unwahrscheinlicher ist, muss nun die Schwelle nach links &ndash; also weg vom wahrscheinlicheren Symbol &ndash; 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)'''&nbsp; Mit dem Ergebnis aus (6) gilt schließlich:
+
'''(7)'''&nbsp; 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}.$$
 
:$$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
+
*Aufgrund der weniger gravierenden Unsymmetrie ist die erreichbare Verbesserung mit $9 \%$ geringer als unter Punkt '''(5)''' berechnet
  
 
{{ML-Fuß}}
 
{{ML-Fuß}}

Version vom 10. Juni 2020, 16:50 Uhr


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} ({\it 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:

$$\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

1

Welcher Zusammenhang besteht zwischen den Funktionen  ${\rm Q}(x)$  und  ${\rm erfc}(x)$?

${\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)$.

2

Welche Bitfehlerwahrscheinlichkeit ergibt sich mit  $\underline{p_{\rm L} = 0.88}$  und  $\underline{E = 0}$?

$p_{\rm B} \ =\ $

$\ \%$

3

Welche Bitfehlerwahrscheinlichkeit ergibt sich mit  $p_{\rm L} = 0.88$  und  $\underline{E = 0.1 \hspace{0.05cm} \rm V}$?

$p_{\rm B} \ =\ $

$\ \%$

4

Bestimmen Sie den optimalen Schwellenwert  $E_{\rm opt}$  für  $p_{\rm L} = 0.88$   ⇒   $p_{\rm H} = 0.12$.

$E_{\rm opt} \ =\ $

$\ \rm V$

5

Es gelte  $E = E_{\rm opt}$. Wie groß ist die minimale Bitfehlerwahrscheinlichkeit mit  $p_{\rm L} = 0.88$?

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

$\ \%$

6

Bestimmen Sie den optimalen Schwellenwert  $E_{\rm opt}$  für  $\underline{p_{\rm L} = 0.31}$   ⇒   $p_{\rm H} = 0.69$.

$E_{\rm opt} \ =\ $

$\ \rm V$

7

Wie groß ist die minimale Bitfehlerwahrscheinlichkeit für diesen Fall  $(p_{\rm L} = 0.31)$?

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

$\ \%$


Musterlösung

(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