imported>Xenein: /* growthexperiments-addlink-summary-summary:3|0|0 */

2024-11-28T17:13:38Z

growthexperiments-addlink-summary-summary:3|0|0

Neue Seite

Die '''Reed-Muller-Codes''' sind eine Familie von [[Linearer Code|linearen]], [[Fehlerkorrekturverfahren|fehlerkorrigierenden]] Codes, die im Bereich der [[Kanalcodierung]] zur gesicherten Datenübertragung und Datenspeicherung Verwendung finden. Diese Klasse von Codes wurden von [[Irving S. Reed]] und David E. Muller entwickelt.

== Praxis ==

Der binäre Reed-Muller-Code wurde von der [[NASA]] in den [[Mariner]] Expeditionen (1969 bis 1976) zum [[Mars (Planet)|Mars]] benutzt, um die vom Mars gemachten Fotos an die [[Erde]] zu senden. Im Speziellen wurde bei Mariner 9 ein RM-Code (1, 5) ohne Kontrollmatrix als [[Hadamard-Code]] (32, 6, 16) verwendet, das bedeutet, dass sechs Informationsbits in 32 Bit langen Wörtern kodiert waren und das [[Minimalgewicht]] der Wörter mindestens 16 betrug, was eine [[Fehlerkorrektur]] von 7 Bits ermöglichte. Mit den <math>2^6=64</math> Codewörtern wurden Grauwerte eines Bildpunktes kodiert. Mehr dazu im nachfolgenden Beispiel 3 zur NASA Raumsonde Mariner 9.

== Konstruktion ==
Im Folgenden wird beschrieben, wie man eine Erzeugermatrix eines Reed-Muller-Codes der Länge <math>n=2^d</math> konstruiert

:<math> X = \mathbb{F}_2^d = \{ x_1, \ldots, x_{2^d} \} </math>.

<math>\mathbb{F}_n</math> ist eine [[Teilmenge]] der [[Natürliche Zahl|nichtnegativen ganzen Zahlen]]

:<math>\mathbb{F}_n = \{a \in \mathbb{N}_0\;|\; a < n\}</math>.

Wir definieren im n-dimensionalen Raum <math>\mathbb{F}_2^n</math> die Indikatorvektoren :

:<math>\mathbb{I}_A \in \mathbb{F}_2^n</math>

auf Untermengen <math>A \subset X</math> durch:

:<math>\left( \mathbb{I}_A \right)_i = \begin{cases} 1 & \mbox{ wenn } x_i \in A \\ 0 & \mbox{ sonst} \\ \end{cases} </math>

und – ebenfalls in <math>\mathbb{F}_2^n</math> – die binäre Operation:

:<math> w \wedge z = (w_1 \times z_1, \ldots , w_n \times z_n ) </math>

die als ''Keil-Produkt'' bezeichnet wird.

<math>\mathbb{F}_2^d</math> ist ein <math>d</math>-dimensionaler [[Vektorraum]] über <math>\mathbb{F}_2</math>, deshalb ist es möglich zu schreiben:

<math>(\mathbb{F}_2)^d = \{ (y_1, \ldots , y_d)\;|\;y_i \in \mathbb{F}_2 \} </math>

Wir definieren im <math>n</math>-dimensionalen Raum <math>\mathbb{F}_2^n</math> die folgenden Vektoren der Länge <math>n</math>

<math>v_0 = (1,1,\ldots,1)</math> und

:<math> v_i = \mathbb{I}_{ H_i } </math>,

wobei <math>H_i</math> [[Hyperebene]]n in <math>(\mathbb{F}_2)^d</math> (mit Dimension <math>d-1</math>) sind:

:<math>H_i = \{ y \in ( \mathbb{F}_2 ) ^d \mid y_i = 0 \} </math>

Der '''Reed-Muller RM(d, r)-Code''' der Ordnung <math>r</math> und der Länge <math>n=2^d</math> ist derjenige Code, der durch <math>v_0</math> und dem Keil-Produkt von bis zu <math>r</math> der <math>v_i</math> erzeugt wird (wobei nach Vereinbarung ein Keil-Produkt von weniger als einem Vektor gleich der Identität für diesen Operator ist).
== Eigenschaften ==
Es gelten die folgenden Eigenschaften

# Die Menge aller möglichen Keil-Produkte von bis zu ''d'' der <math>v_i</math> bilden eine [[Basis (Vektorraum)|Basis]] von <math>\mathbb{F}_2^n</math>.
# Der RM(d, r)-Code hat den Rang: <math> \sum_{s=0}^r {d \choose s} </math>
# Es gilt <math>RM(d, r) = RM(d-1, r)|RM(d-1, r-1) </math>, wobei <math>|</math> das [[Bar-Product (Codierungstheorie)|Bar-Product]] zweier Codes bezeichnet
# RM(d, r) hat die minimale [[Hamming-Distanz]] <math>2^{d-r}</math>.

== Beispiel 1 ==

Sei <math>d=3</math>. Dann <math>n=8</math>, und

:<math> X = \mathbb{F}_2^3 = \{ (0,0,0), (0,0,1), \ldots, (1,1,1) \}. </math>

und

:<math>
\begin{matrix}
v_0 & = & (1,1,1,1,1,1,1,1) \\
v_1 & = & (1,0,1,0,1,0,1,0) \\
v_2 & = & (1,1,0,0,1,1,0,0) \\
v_3 & = & (1,1,1,1,0,0,0,0) \\
\end{matrix}
</math>

Der RM(3,1)-Code wird erzeugt durch die Menge

:<math> \{ v_0, v_1, v_2, v_3 \} </math>

oder genauer durch die Zeilen der Matrix

:<math>
\begin{pmatrix}
1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 \\
1 & 0 & 1 & 0 & 1 & 0 & 1 & 0 \\
1 & 1 & 0 & 0 & 1 & 1 & 0 & 0 \\
1 & 1 & 1 & 1 & 0 & 0 & 0 & 0 \\
\end{pmatrix}
</math>

== Beispiel 2 ==

Der RM(3,2)-Code wird erzeugt durch die [[Menge (Mathematik)|Menge]]
:<math> \{ v_0, v_1, v_2, v_3, v_1 \wedge v_2, v_1 \wedge v_3, v_2 \wedge v_3 \} </math>

oder genauer durch die Zeilen der [[Matrix (Mathematik)|Matrix]]

:<math>
\begin{pmatrix}
1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 \\
1 & 0 & 1 & 0 & 1 & 0 & 1 & 0 \\
1 & 1 & 0 & 0 & 1 & 1 & 0 & 0 \\
1 & 1 & 1 & 1 & 0 & 0 & 0 & 0 \\
1 & 0 & 0 & 0 & 1 & 0 & 0 & 0 \\
1 & 0 & 1 & 0 & 0 & 0 & 0 & 0 \\
1 & 1 & 0 & 0 & 0 & 0 & 0 & 0 \\
\end{pmatrix}
</math>

== Beispiel 3: NASA Raumsonde Mariner 9 ==

Bei der NASA [[Raumsonde]] [[Mariner#Mariner 8 und 9|Mariner 9]] aus dem Jahre 1971 wurde ein Reed-Muller-Code (1, 5) mit fehlender Kontrollmatrix genutzt, der einen Spezialfall allgemeiner Reed-Muller Codes darstellt. Dieser Code war schlussendlich ein [[Hadamard-Code]] mit den Parametern (32, 6, 16). Mit diesem RM-Code (32, 6, 16) wurden 32 Bit lange Codewörter übertragen, die <math>2^6=64</math> Werte kodierten, wobei die Codewörter untereinander einen Hamming-Abstand von 16 aufwiesen. Diese Parameter wurden aufgrund der Kanalcharakteristik, der Bildauflösung und der Aufnahme- und Übertragungszeiten gewählt, die eine Wortlänge von reichlich 30 Bit sinnvoll machten.

Aufgrund der großen Entfernung zwischen Mars und Erde, und den damals im Vergleich zu heute unfortschrittlichen Übertragungsgeräten, lag die angenommene Bit-Fehlerwahrscheinlichkeit bei 5 %. Daraus ergibt sich aufgrund der Kodierung von einem [[Grauwert]] in 6 Bit ohne zusätzliche Fehlerkorrekturmechanismen eine Grauwert-Fehlerwahrscheinlichkeit von 26 %. Das heißt, ca. ein Viertel eines übertragenen Bildes kommt fehlerhaft beim Empfänger an. Durch den Einsatz des RM-Code (32, 6, 16) konnte bei gleicher Bit-Fehlerwahrscheinlichkeit von 5 % die Grauwert-Fehlerwahrscheinlichkeit auf 0,01 % reduziert werden.

=== Konstruktion ===

[[Datei:Hadamard-Code.svg|mini|Matrix des Hadamard-Code (32, 6, 16) für den Reed-Muller-Code (1,5) der NASA Raumsonde Mariner 9 (1971/1972). Die Farbe Schwarz kodiert die Binärziffer 1, und die Farbe Weiß kodiert die Binärziffer 0.]]

Der verwendete RM-Code (32, 6, 16) basiert auf einer [[Hadamard-Matrix]] <math>H_{32}</math>.

Die Konstruktion von <math>H_{32}</math> erfolgt rekursiv aus der Hadamard-Matrix
:<math>H_1 = \begin{pmatrix} 1 \end{pmatrix}</math>
und der Erzeugungsregel
:<math>H_{2n} = \begin{pmatrix} H_n & H_n \\ H_n & -H_n \end{pmatrix}</math>
Diese Konstruktion nach [[James Joseph Sylvester|Sylvester]] erzeugt die sogenannten [[Walsh Matrix|Walsh Matrizen]]
:<math>H_1 = \begin{pmatrix} 1 \end{pmatrix}, H_2 = \begin{pmatrix} 1 & 1 \\ 1 & -1 \end{pmatrix}, H_4 = \begin{pmatrix} 1 & 1 & 1 & 1 \\ 1 & -1 & 1 & -1 \\ 1 & 1 & -1 & -1 \\ 1 & -1 & -1 & 1 \end{pmatrix}, \ldots </math>
die normalisierte Hadamard-Matrizen vom Grad <math>2^k</math> darstellen.

Wenn man nun die Hadamard-Matrix <math>H_{32}</math> als Bitmuster interpretiert (bei dem eine 1 für die Binärziffer 1, und eine <math>-1</math> für die Binärziffer 0 steht), dann erhält man 32 Codewörter mit 32 Bit Länge. Jedes dieser Codewörter weist zu jedem anderen Codewort einen Hamming-Abstand von 16 oder 32 auf. Durch Kombination der Hadamard-Matrix <math>H_{32}</math> mit der dazu inversen Hadamard-Matrix <math>-H_{32}</math> erhält man 64 Codewörter mit 32 Bit Länge, bei denen jedes Codewort zu jedem anderen Codewort einen Hamming-Abstand von 16 aufweist. Diese Kombination von <math>H_{32}</math> und <math>-H_{32}</math> definiert dabei einen Hadamard-Code, mit dem sich <math>2^6=64</math> Werte kodieren lassen, indem ein Wert <math>n</math> der <math>n</math>-ten Zeile des Codes entspricht. Die nebenstehende Abbildung zeigt den vollständigen Hadamard-Code für RMC (32, 6, 16), wobei die Farbe Schwarz für die Binärziffer 1 und die Farbe Weiß für die Binärziffer 0 steht.

== Alternative Charakterisierung ==
Die Klasse der Reed-Muller-Codes kann man auch mit einer Menge von Abbildungen identifizieren. Betrachte hierzu die Menge
:<math> V = \{ f \text{ Abbildung} \mid f\colon \mathbb{F}_2^m \rightarrow \mathbb{F}_2 \} </math>.
Eine Abbildung <math>f \in V</math> wird durch ihre <math>{2^m}</math> Bilder eindeutig bestimmt, sofern deren Reihenfolge bekannt ist. Daher kann man <math>f</math> auch durch den zugehörigen Bildvektor <math>(f(0),f(1),\dots,f(2^m-1)) \in \mathbb{F}_2^{2^m}</math> darstellen, wobei die Argumente <math>0,1,\dots,2^m-1</math> die <math>2</math>-adische Entwicklung der Elemente aus dem Definitionsbereich <math>\mathbb{F}_2^m</math> sind.
Auf <math>V</math> kann man eine komponentenweise Addition und Multiplikation gemäß den Rechenoperationen in <math>\mathbb{F}_2</math> definieren. Genau genommen gibt es einen Ringisomorphismus zwischen der Menge der Abbildungen <math>V</math> und der Menge der Bildvektoren <math>\mathbb{F}_2^{2^m}</math>, weshalb man eine Abbildung auch mit seinem Bildvektor identifizieren kann: <math> f = (f(0),f(1),\dots,f(2^m-1)) </math>. In <math>V</math> liegen spezielle Abbildungen, die sogenannten [[Koordinatenfunktion|Koordinatenfunktionen]] <math>Z_i,\;i \in \{1 \dots 2^m\}</math>.

Diese sind wie folgt definiert:

:<math>Z_i(v) := v_i</math> für <math>v=(v_1,\dots,v_m) \in \mathbb{F}_2^m</math>.
In der oben eingeführten Vektordarstellung lassen sich die Koordinatenfunktionen auch schreiben als
:<math>Z_i = (\underbrace{0,\dots,0}_{2^{i-1}\text{-mal}},\underbrace{1,\dots,1}_{2^{i-1}\text{-mal}},\underbrace{0,\dots,0}_{2^{i-1}\text{-mal}},\dots) \in \mathbb{F}_2^{2^m}</math>.

Nun gilt:
# Das System der Monome <math>Z_{i_1} \cdot \dots \cdot Z_{i_k}</math> (<math>1 \le i_1 < \dots < i_k \leq m, k=0,\dots,m</math>) ist eine Basis von <math>V</math>.
# Die Teilmenge <math> \{ f\colon \mathbb{F}_2^m \rightarrow \mathbb{F}_2 \text{ Abbildung} \mid \operatorname{grad}(f) \leq r \} \subseteq V </math> entspricht dem Reed-Muller-Code RM(r, m). Hierbei ist <math> \operatorname{grad}(f) </math> der höchste Monomgrad der Koordinatenfunktionen, als deren Summe <math>f</math> gemäß 1. geschrieben werden kann.

== Dekodierung ==
Die Idee ist wie folgt: Jedes Codewort des Reed-Muller-Codes RM(r,m) kann gemäß der obigen alternativen Charakterisierung als Funktion <math>f</math> aus <math>V</math> aufgefasst werden – mit Basisdarstellung in entgegengesetzten Koordinatenfunktionen, d. h. mit eindeutig bestimmten Koeffizienten <math>m_I \text{ mit } I \subseteq M</math> wobei <math>M = \{ 1,\dots,m \}</math> die Menge der Koordinatenfunktionen-Indizes ist. Die Funktion <math>f</math> wird als Bildvektor <math>(f(0),f(1),\dots,f(2^m-1))</math> durch den gestörten Kanal geschickt. Der Empfänger dekodiert das mit Fehler <math>e</math> behaftete Codewort <math>g=f+e</math>, indem er sukzessive die Koeffizienten <math>m_I</math> rekonstruiert. Dabei beginnt er mit den Koeffizienten, die zum [[Monom]] höchsten Grades <math>r</math> gehören.
Hierzu berechnet er das Skalarprodukt von <math>g</math> mit den s.g. charakteristischen Funktionen des Monoms. Dies sind alle Abbildungen vom Grad <math>m-r</math>, wobei die erzeugenden Koordinatenfunktionen auch entgegengesetzt vorkommen können. Der Wert, der mehrheitlich durch die Skalarprodukte berechnet wird, ist der ursprüngliche Monomkoeffizient. Das Verfahren wird mit den Monomen vom Grad <math>r-1,r-2,\dots,0</math> wiederholt und man erhält hierdurch schließlich <math>f</math> – vorausgesetzt der Fehler <math>e</math> ist nicht zu groß.

== Zusammenfassung ==
Codierungs- und Decodierungsprozess mittels Reed-Muller-Codes im Überblick:
# Nachricht <math>n</math> wird in Codewort <math>c</math> übersetzt.
# Codewort <math>c</math> kann mit Abbildung <math>f</math> identifiziert werden.
# Abbildung <math>f</math> kann auch als Bildvektor <math>(f(0),f(1),\dots,f(2^m-1))</math> dargestellt werden.
# Übermittle anstelle der Monomkoeffizienten von <math>f</math> den zugehörigen Bildvektor. Dies liefert Redundanz, die die gewünschte Fehlerkorrektur ermöglicht.
# Sende den Bildvektor durch den gestörten Kanal. Es ergibt sich <math>g=f+e</math> mit Fehlervektor <math>e</math>.
# Empfange den Bildvektor <math>g</math> und gewinne durch [[Skalarmultiplikation]]en mit den Koordinatenfunktionen <math>Z_i</math> die ursprünglichen Monomkoeffizienten.
# Durch die Monomkoeffizienten berechnet man die/das ursprüngliche Abbildung/Codewort <math>f=c</math> und damit <math>n</math>.

== Weblinks ==

* [http://elib.tu-darmstadt.de/diss/000183/stolte.pdf Rekursive Codes mit der Plotkin-Konstruktion] (PDF; 1,7 MB) Dissertation zur Konstruktion und Decodierung von Reed-Muller Codes und deren Untercodes (Achtung: Angabe über den RM-Code (32, 6, 16) der Mariner 9 Mission sind nicht korrekt, da nur eine Mächtigkeit des Codes von <math>2^5=32</math> Werten angegeben und erläutert wird.)

[[Kategorie:Kodierungstheorie]]
[[Kategorie:Binärcode]]

Reed-Muller-Code - Versionsgeschichte

imported>Xenein: /* growthexperiments-addlink-summary-summary:3|0|0 */