Exercise 4.4Z: Supplement to Exercise 4.4
From LNTwww
Der Informationstheoretiker Robert G. Gallager hat sich bereits 1963 mit folgender Fragestellung beschäftigt:
- Gegeben ist ein Zufallsvektor x_=(x1,x2, ... ,xn) mit n binären Elementen xi∈{0,1}.
- Bekannt sind alle Wahrscheinlichkeiten pi=Pr(xi=1) und qi=Pr(xi=0)=1−pi mit Inex i=1, ... , n.
- Gesucht ist die Wahrscheinlichkeit, dass die Anzahl der Einsen in diesem Vektor geradzahlig ist.
- Oder ausgedrückt mit dem Hamming–Gewicht: Wie groß ist die Wahrscheinlichkeit Pr[wH(x_) ist gerade]?
Die Grafik verdeutlicht die Aufgabenstellung für das Beispiel n=4 sowie p1=0.2, p2=0.9, p3=0.3 und p4=0.6.
- Für die grün hinterlegte Zeile ⇒ x_=(1,0,0,1) gilt wH(x_)=2 und Pr(x_)=p1⋅q2⋅q3⋅p4=0.0084.
- Blaue Schrift bedeutet ein geradzahliges Hamming–Gewicht. Rote Schrift steht für „wH(x_) ist ungerade”.
- Die Wahrscheinlichkeite Pr[wH(x_) ist gerade] ist gleich der Summe der blauen Zahlen in der letzten Spalte. Die Summe der roten Zahlen ergibt Pr[wH(x_) ist ungerade]=1−Pr[wH(x_ ist gerade].
Gallager hat das Problem in analytischer Weise gelöst:
- Pr[wH(x_)istgerade] = 1/2⋅[1+π],
- Pr[wH(x_)istungerade] = 1/2⋅[1−π].
Fragebogen
Musterlösung
(1)
Entsprechend nebenstehender Tabelle gilt:
- Pr[wH(x_)istgerade]=Pr[wH=0]+Pr[wH=2].
Mit den Wahrscheinlicekiten
- p1=Pr(x1=1) = 0.2,q1=Pr(x1=0)=0.8,
- p2=Pr(x2=1) = 0.9,q2=Pr(x2=0)=0.1
erhält man:
- Pr[wH(x_)=0] = Pr[(x1=0)∩(x2=0)]=q1⋅q2=0.8⋅0.1=0.08,
- Pr[wH(x_)=2] = Pr[(x1=1)∩(x2=1)]=p1⋅p2=0.2⋅0.9=0.18
- ⇒Pr[wH(x_)istgerade]=0.8+0.18=0.26_.
Die Gallager–Gleichung liefert für den gleichen Parametersatz:
- Pr[wH(x_)istgerade] = 0.5+0.5⋅2∏i=1(1−2⋅pi)=
- = 0.5+0.5⋅(1−2⋅0.2)⋅(1−2⋅0.9)=0.26.
Die von Gallager 1963 angegebene Gleichung wurde hiermit für n=2 verifiziert.
- Pr[wH(x_)istgerade]=0.056+
- +0.216+0.006+0.126=0.404_.
Die roten Zeilen liefern das Komplementärereignis:
- Pr[wH(x_)istungerade]=0.024+
- +0.504+0.014+0.054=0.596.
Die Gallager–Gleichung liefert auch hier wieder das exakt gleiche Ergebnis.
- {\rm Pr} \left [w_{\rm H}(\underline{x})\hspace{0.15cm}{\rm ist \hspace{0.15cm} gerade}\right] \hspace{-0.15cm} \ = & \hspace{-0.15cm} 0.5 + 0.5 \cdot \prod\limits_{i =1}^{3} \hspace{0.25cm}(1-2\cdot p_i) =
- = 0.5+0.5⋅(+0.6)⋅(−0.8)⋅(+0.4)=0.404.
(3)
(4)
(5)