Difference between revisions of "Digital Signal Transmission/Error Probability for Baseband Transmission"

From LNTwww
Line 70: Line 70:
  
  
== Fehlerwahrscheinlichkeit bei Gaußsch (..) n (1) ==
+
== Fehlerwahrscheinlichkeit bei Gaußschem Rauschen (1) ==
 
<br>
 
<br>
Entsprechend den [http://en.lntwww.de/Digitalsignal%C3%BCbertragung/Systemkomponenten_eines_Basisband%C3%BCbertragungssystems#Ersatzschaltbild und Voraussetzungen für Kapitel 1 Voraussetzungen] zu diesem Kapitel gehen wir davon aus, dass das Detektionssignal zu den Detektionszeitpunkten wie folgt dargestellt werden kann:
+
Entsprechend den [http://en.lntwww.de/Digitalsignal%C3%BCbertragung/Systemkomponenten_eines_Basisband%C3%BCbertragungssystems#Ersatzschaltbild_und_Voraussetzungen_für_Kapitel_1 Voraussetzungen] zu diesem Kapitel gehen wir davon aus, dass das Detektionssignal zu den Detektionszeitpunkten wie folgt dargestellt werden kann::
 +
<math> d(\nu  T) = d_{\rm S}(\nu  T)+d_{\rm N}(\nu T)\hspace{0.05cm}. </math>
 +
Der Nutzanteil wird durch die Wahrscheinlichkeitsdichtefunktion (WDF) <i>f</i><sub><i>d</i>S</sub>(<i>d</i><sub>S</sub>) beschrieben, wobei wir hier von unterschiedlichen Auftrittswahrscheinlichkeiten <i>p</i><sub>L</sub> = Pr(<i>d</i><sub>S</sub> = &ndash;<i>s</i><sub>0</sub>),
 +
<i>p</i><sub>H</sub> = Pr(<i>d</i><sub>S</sub> = +<i>s</i><sub>0</sub>) = 1&ndash; <i>p</i><sub>L</sub> ausgehen. Die WDF <i>f</i><sub><i>d</i>N</sub>(<i>d</i><sub>N</sub>) der Störkomponente ist gaußförmig und besitzt die Streuung <i>&sigma;<sub>d</sub></i>.<br>
 +
[[File:P_ID1259__Dig_T_1_2_S3_v2.png|left]]<br>
 +
 
 +
Die WDF <i>f</i><sub><i>d</i></sub>(<i>d</i>) der Detektionsabtastwerte <i>d</i>(<i>&nu;</i><i>T</i>) ergibt sich unter der Voraussetzung, dass <i>d</i><sub>S</sub>(<i>&nu;</i><i>T</i>) und <i>d</i><sub>N</sub>(<i>&nu;</i><i>T</i>) statistisch unabhängig voneinander sind (&bdquo;signalunabhängiges Rauschen&rdquo;), als Faltungsprodukt<nowiki>:</nowiki>
 +
 
 +
<math>f_d(d) = f_{d{\rm S}}(d_{\rm S}) \star f_{d{\rm N}}(d_{\rm N})\hspace{0.05cm}.</math><br><br>
 +
Der Schwellenwertentscheider mit der Schwelle <i>E</i> = 0 trifft dann eine falsche Entscheidung, wenn
 +
*das Symbol <b>L</b> gesendet wurde (<i>d</i><sub>S</sub> = &ndash;<i>s</i><sub>0</sub>) und <i>d</i> > 0 ist (rote schraffierte Fläche), oder
 +
*das Symbol <b>H</b> gesendet wurde (<i>d</i><sub>S</sub> = +<i>s</i><sub>0</sub>) und <i>d</i> < 0 ist (blaue schraffierte Fläche).
 +
<br>Da die Flächen der zwei Gaußkurven zusammen 1 ergeben, gibt die Summe aus der rot und der blau schraffierten Fläche die Bitfehlerwahrscheinlichkeit <i>p</i><sub>B</sub> an. Die beiden grün schraffierten Flächen in der oberen Wahrscheinlichkeitsdichtefunktion <i>f</i><sub><i>d</i>N</sub>(<i>d</i><sub>N</sub>) sind &ndash; jede für sich &ndash; ebenfalls gleich <i>p</i><sub>B</sub>.<br>
 +
Die anhand der Grafik veranschaulichten Ergebnisse sollen nun formelmäßig hergeleitet werden. Es gilt:
 +
<math>p_{\rm B} = p_{\rm L} \cdot {\rm Pr}( \upsilon_\nu = \mathbf{H}\hspace{0.1cm}|\hspace{0.1cm} q_\nu = \mathbf{L})+
 +
  p_{\rm H} \cdot {\rm Pr}( \upsilon_\nu = \mathbf{L}\hspace{0.1cm}|\hspace{0.1cm} q_\nu = \mathbf{H})\hspace{0.05cm}.</math>
 +
<br>Hierbei sind <i>p</i><sub>L</sub> und <i>p</i><sub>H</sub> die Quellensymbolwahrscheinlichkeiten, während die jeweils zweiten, bedingten Wahrscheinlichkeiten Pr(<i>&upsilon;<sub>&nu;</sub></i> | <i>q<sub>&nu;</sub></i>) die Verfälschungen durch den AWGN&ndash;Kanal beschreiben.
 +
Aus der Entscheidungsregel des Schwellenwertentscheiders (mit Schwelle <i>E</i> = 0) ergibt sich auch::
 +
<math>p_{\rm B} = p_{\rm L} \cdot {\rm Pr}( d(\nu T)>0)+  p_{\rm H} \cdot {\rm Pr}( d(\nu  T)<0) =</math>
 +
::<math>= p_{\rm L} \cdot {\rm Pr}( d_{\rm N}(\nu T)>s_0)+  p_{\rm H} \cdot {\rm Pr}( d_{\rm N}(\nu T)<-s_0) \hspace{0.05cm}.</math>
 +
Die Herleitung wird auf der nächsten Seite fortgesetzt.
 +
 
  
  

Revision as of 19:03, 16 November 2016


Definition der Bitfehlerwahrscheinlichkeit


Die Grafik zeigt ein sehr einfaches, aber allgemeingültiges Modell eines binären Übertragungssystems.

P ID1258 Dig T 1 2 S1 v1.png


Dieses lässt sich wie folgt charakterisieren:

  • Die Quelle und die Sinke werden durch die beiden Binärfolgen 〈qν〉 und 〈υν〉 beschrieben.
  • Das gesamte Übertragungsystem – bestehend aus Sender, Übertragungskanal inklusive Störungen und Empfänger – wird als „Black Box” mit binärem Ein– und Ausgang betrachtet.
  • Dieser „Digitale Kanal” wird allein durch die Fehlerfolge 〈eν〉 charakterisiert. Bei fehlerfreier Übertragung des ν–ten Bits (υν = qν) gilt eν = 0, andernfalls (υνqν) wird eν = 1 gesetzt.

: Die (mittlere) Bitfehlerwahrscheinlichkeit ist bei einem Binärsystem wie folgt gegeben:\[\it p_{\rm B} = \rm E[\rm Pr(\it v_{\nu} \ne q_{\nu})]= \overline {\rm Pr(\it v_{\nu} \ne q_{\nu})} = \lim_{{\it N}\to\infty}\frac{\rm 1}{\it N}\cdot\sum\limits_{\it \nu=\rm 1}^{\it N}\rm \rm Pr(\it v_{\nu} \ne q_{\nu})\hspace{0.05cm}.\]
Diese statistische Größe ist das wichtigste Beurteilungskriterium eines jeden Digitalsystems.



Die Berechnung als Erwartungswert E[…..] gemäß dem ersten Teil der obigen Gleichung entspricht einer Scharmittelung über die Verfälschungswahrscheinlichkeit Pr(υνqν) des ν–ten Symbols, während die überstreichende Linie im rechten Teil eine Zeitmittelung kennzeichnet. Beide Berechnungsarten führen – unter der gerechtfertigten Annahme ergodischer Prozesse – zum gleichen Ergebnis, wie im Kapitel 4 des Buches „Stochastische Signaltheorie” gezeigt wurde.

Auch aus der Fehlerfolge 〈eν〉 lässt sich die Bitfehlerwahrscheinlichkeit als Erwartungswert bestimmen, wobei zu berücksichtigen ist, dass eν nur die Werte 0 und 1 annehmen kann:\[\it p_{\rm B} = \rm E[\rm Pr(\it e_{\nu}=\rm 1)]= {\rm E}[{\it e_{\nu}}]\hspace{0.05cm}.\] Die obige Definition der Bitfehlerwahrscheinlichkeit gilt unabhängig davon, ob es statistische Bindungen innerhalb der Fehlerfolge 〈eν〉 gibt oder nicht. Je nachdem ist der Aufwand zur Berechnung von pB unterschiedlich groß und bei einer Systemsimulation müssen unterschiedliche digitale Kanalmodelle herangezogen werden. Im Kapitel 5 wird gezeigt, dass das sog. BSC–Modell (Binary Symmetrical Channel) statistisch unabhängige Fehler liefert, während für die Beschreibung von Bündelfehlerkanälen auf die Modelle von Gilbert–Elliott: .: Capacity of Burst–Noise Channel, In: Bell Syst. Techn. J. Vol. 39, 1960, pp. 1253–1266 and McCullough : The Binary Regenerative Channel, In: Bell Syst. Techn. J. (47), 1968 zurückgegriffen werden muss.


Definition der Bitfehlerquote (1)


Die Bitfehlerwahrscheinlichkeit pB eignet sich zum Beispiel gut für die Konzipierung und Optimierung von Digitalsystemen. Diese ist eine Apriori-Kenngröße, die eine Vorhersage über das Fehlerverhalten eines Nachrichtensystems erlaubt, ohne dass dieses bereits realisiert sein muss. Dagegen muss zur messtechnischen Erfassung der Qualität eines realisierten Systems oder bei einer Systemsimulation auf die Bitfehlerquote übergegangen werden, die durch den Vergleich von Quellen– und Sinkensymbolfolge ermittelt wird. Diese ist somit eine Aposteriori-Kenngröße des Systems.

:Die Bitfehlerquote (englisch: Bit Error Rate, BER) ist das Verhältnis aus der Anzahl nB(N) der aufgetretenen Bitfehler (υνqν) und der Anzahl N der insgesamt übertragenen Symbole\[h_{\rm B}(N) = \frac{n_{\rm B}(N)}{N} \hspace{0.05cm}.\] Im Sinne der Wahrscheinlichkeitsrechnung stellt die Bitfehlerquote eine relative Häufigkeit dar; sie wird demzufolge auch Bitfehlerhäufigkeit genannt.


Die Schreibweise hB(N) soll deutlich machen, dass die per Messung oder durch Simulation ermittelte Bitfehlerquote signifikant von dem Parameter N – also der Anzahl der insgesamt übertragenen oder simulierten Symbole – abhängt. Nach den elementaren Gesetzen der Wahrscheinlichkeitsrechnung stimmt nur im Grenzfall N → ∞ die Aposteriori–Kenngröße hB mit der Apriori–Kenngröße pB exakt überein.

Der Zusammenhang zwischen Wahrscheinlichkeit und relativer Häufigkeit wird in einem Lernvideo zum Buch „Stochastische Signaltheorie” verdeutlicht:
Das Bernoullische Gesetz der großen Zahlen (Dateigröße: 1.97 MB, Dauer: 4:25)

Für die nachfolgende Herleitung wird das BSC–Modell zugrunde gelegt, das in Kapitel 5.2 im Detail beschrieben wird. Jedes einzelne Bit wird mit der Wahrscheinlichkeit p = Pr(υνqν) = Pr(eν = 1) verfälscht, unabhängig von den Fehlerwahrscheinlichkeiten der benachbarten Symbole. Die (mittlere) Bitfehlerwahrscheinlichkeit pB ist somit ebenfalls gleich p.


Definition der Bitfehlerquote (2)


Nun soll abgeschätzt werden, wie genau die Bitfehlerwahrscheinlichkeit pB = p beim BSC-Modell durch die Bitfehlerquote hB approximiert wird. Dies geschieht in mehreren Schritten:

  • Die Anzahl der Bitfehler bei der Übertragung von N Symbolen ist eine diskrete Zufallsgröße:\[n_{\rm B}(N) = \sum\limits_{\it \nu=\rm 1}^{\it N} e_{\nu} \hspace{0.2cm} \in \hspace{0.2cm} \{0, 1, ... , N \}\hspace{0.05cm}.\]
  • Bei statistisch unabhängigen Fehlern (BSC) ist nB binominalverteilt. Demzufolge gilt:\[m_{n{\rm B}}=N \cdot p_{\rm B},\hspace{0.2cm}\sigma_{n{\rm B}}=\sqrt{N\cdot p_{\rm B}\cdot (\rm 1- \it p_{\rm B})}\hspace{0.05cm}.\]
  • Für Mittelwert und Streuung der Bitfehlerquote hB = nB/N gilt deshalb:\[m_{h{\rm B}}= \frac{m_{n{\rm B}}}{N} = p_{\rm B}\hspace{0.05cm},\hspace{0.2cm}\sigma_{h{\rm B}}= \frac{\sigma_{n{\rm B}}}{N}= \sqrt{\frac{ p_{\rm B}\cdot (\rm 1- \it p_{\rm B})}{N}}\hspace{0.05cm}.\]
  • Nach Moivre und Laplace lässt sich die Binominalverteilung in eine Gaußverteilung überführen:\[f_{h{\rm B}}({h_{\rm B}}) \approx \frac{1}{\sqrt{2\pi}\cdot\sigma_{h{\rm B}}}\cdot {\rm exp} \left[-\frac{(h_{\rm B}-p_{\rm B})^2}{2 \cdot \sigma_{h{\rm B}}^2}\right].\]
  • Mit dem Gaußschen Fehlerintergal Q(x) lässt sich somit die Wahrscheinlichkeit pε berechnen, dass die per Simulation/Messung über N Symbole ermittelte Bitfehlerquote hB(N) betragsmäßig um weniger als einen Wert ε von der tatsächlichen Bitfehlerwahrscheinlichkeit pB abweicht:\[p_{\varepsilon}= {\rm Pr} \left( |h_{\rm B}(N) - p_{\rm B}| < \varepsilon \right) = 1 -2 \cdot {\rm Q} \left( \frac{\varepsilon}{\sigma_{h{\rm B}}} \right)= 1 -2 \cdot {\rm Q} \left( \frac{\varepsilon \cdot \sqrt{N}}{\sqrt{p_{\rm B} \cdot (1-p_{\rm B})}} \right)\hspace{0.05cm}.\]

Dieses Ergebnis ist wie folgt zu interpretieren: Wenn man unendlich viele Versuchsreihen über jeweils N Symbole durchführt, ist der Mittelwert mhB tatsächlich gleich der gesuchten Fehlerwahrscheinlichkeit pB. Bei einer einzigen Versuchsreihe wird man dagegen nur eine Näherung erhalten, wobei die jeweilige Abweichung vom Sollwert bei mehreren Versuchsreihen gaußverteilt ist.

: Die Bitfehlerwahrscheinlichkeit betrage pB = 10–3 und es ist bekannt, dass die Bitfehler statistisch unabhängig sind. Macht man nun sehr viele Versuchsreihen mit jeweils N = 105 Symbolen, so werden die jeweiligen Ergebnisse hB entsprechend einer Gaußverteilung um den Sollwert 10–3 variieren. Die Streuung beträgt dabei

\(\sigma_{h{\rm B}}= \sqrt{{ p_{\rm B}\cdot (\rm 1- \it p_{\rm B})}/{N}}\approx 10^{-4}\hspace{0.05cm}.\)

Die Wahrscheinlichkeit, dass die relative Häufigkeit einen Wert zwischen 0.9 · 10–3 und 1.1 · 10–3 <nobr>(ε = 10–4)</nobr> haben wird, ist somit gleich pε = 1 – 2 · Q(ε/σhB) = 1 – 2 · Q(1) ≈ 68.4%. Soll diese Wahrscheinlichkeit (Genauigkeit) auf 95% gesteigert werden, so müsste N auf 400 000 erhöht werden.


Fehlerwahrscheinlichkeit bei Gaußschem Rauschen (1)


Entsprechend den Voraussetzungen zu diesem Kapitel gehen wir davon aus, dass das Detektionssignal zu den Detektionszeitpunkten wie folgt dargestellt werden kann:\[ d(\nu T) = d_{\rm S}(\nu T)+d_{\rm N}(\nu T)\hspace{0.05cm}. \] Der Nutzanteil wird durch die Wahrscheinlichkeitsdichtefunktion (WDF) fdS(dS) beschrieben, wobei wir hier von unterschiedlichen Auftrittswahrscheinlichkeiten pL = Pr(dS = –s0), pH = Pr(dS = +s0) = 1– pL ausgehen. Die WDF fdN(dN) der Störkomponente ist gaußförmig und besitzt die Streuung σd.

P ID1259 Dig T 1 2 S3 v2.png


Die WDF fd(d) der Detektionsabtastwerte d(νT) ergibt sich unter der Voraussetzung, dass dS(νT) und dN(νT) statistisch unabhängig voneinander sind („signalunabhängiges Rauschen”), als Faltungsprodukt:

\(f_d(d) = f_{d{\rm S}}(d_{\rm S}) \star f_{d{\rm N}}(d_{\rm N})\hspace{0.05cm}.\)

Der Schwellenwertentscheider mit der Schwelle E = 0 trifft dann eine falsche Entscheidung, wenn

  • das Symbol L gesendet wurde (dS = –s0) und d > 0 ist (rote schraffierte Fläche), oder
  • das Symbol H gesendet wurde (dS = +s0) und d < 0 ist (blaue schraffierte Fläche).


Da die Flächen der zwei Gaußkurven zusammen 1 ergeben, gibt die Summe aus der rot und der blau schraffierten Fläche die Bitfehlerwahrscheinlichkeit pB an. Die beiden grün schraffierten Flächen in der oberen Wahrscheinlichkeitsdichtefunktion fdN(dN) sind – jede für sich – ebenfalls gleich pB.
Die anhand der Grafik veranschaulichten Ergebnisse sollen nun formelmäßig hergeleitet werden. Es gilt\[p_{\rm B} = p_{\rm L} \cdot {\rm Pr}( \upsilon_\nu = \mathbf{H}\hspace{0.1cm}|\hspace{0.1cm} q_\nu = \mathbf{L})+ p_{\rm H} \cdot {\rm Pr}( \upsilon_\nu = \mathbf{L}\hspace{0.1cm}|\hspace{0.1cm} q_\nu = \mathbf{H})\hspace{0.05cm}.\]
Hierbei sind pL und pH die Quellensymbolwahrscheinlichkeiten, während die jeweils zweiten, bedingten Wahrscheinlichkeiten Pr(υν | qν) die Verfälschungen durch den AWGN–Kanal beschreiben. Aus der Entscheidungsregel des Schwellenwertentscheiders (mit Schwelle E = 0) ergibt sich auch:\[p_{\rm B} = p_{\rm L} \cdot {\rm Pr}( d(\nu T)>0)+ p_{\rm H} \cdot {\rm Pr}( d(\nu T)<0) =\]

\[= p_{\rm L} \cdot {\rm Pr}( d_{\rm N}(\nu T)>s_0)+ p_{\rm H} \cdot {\rm Pr}( d_{\rm N}(\nu T)<-s_0) \hspace{0.05cm}.\]

Die Herleitung wird auf der nächsten Seite fortgesetzt.