Code de Reed-Muller

Matrice du code de Hadamard augmenté [32, 6, 16] pour le code de Reed-Muller (1, 5) de la sonde spatiale Mariner 9 de la NASA

Les codes de Reed-Muller sont des codes correcteurs linéaires. Cette famille de codes, initialement binaire, doit son nom aux travaux de David E. Muller qui proposa le principe du code et à Irving S. Reed qui proposa une technique de décodage, publiés en 1954. Depuis, cette famille a été largement étudiée et généralisée aux corps finis de plus de 2 éléments. Historiquement, un code Reed-Muller d'ordre 1 en 5 variables, qui a 64 mots de longueur 32 et corrige 7 erreurs, a été utilisé par les sondes Mariner lancées par la NASA entre 1969 et 1973 pour assurer une transmission (numérique) correcte des photos de Mars.

Un code de cette famille est identifié à l'aide de deux paramètres, en général notés $r$ et $m$, appelés respectivement ordre et nombre de variables. Ces paramètres interviennent dans la description utilisant les fonctions booléennes : le code binaire de Reed-Muller d'ordre $r$ en $m$, que l'on note $\mathrm{R}\mathrm{M}(r,m)$, est l'ensemble des tables de vérité des fonctions booléennes en $m$ variables dont la forme algébrique normale (ANF) est de degré au plus $r$. Lorsque l'alphabet est le corps fini à $q$ éléments, il suffit de considérer les fonctions $q$-aires.

Code binaire

Une construction simple

On choisit un ordre quelconque sur les éléments de ${𝔽}_2^m=\{z_0,...,z_{2^m-1}\}$. Une fonction booléenne $f$ en $m$ variables est alors identifiée au mot binaire défini par

$c_f=(f(z_0),...,f(z_{2^m-1}))$

En d'autres termes, $c_f$ est la liste des valeurs prises par $f$ dans un ordre quelconque mais fixe. On peut alors définir

$\mathrm{R}\mathrm{M}(r,m)=\{c_f:d(f)\le r\}$

où $d(f)$ est le degré de l'ANF de $f$.

Exemple de $\mathrm{R}\mathrm{M}(1,5)$

Dans l'exemple donné en introduction, le code Reed-Muller d'ordre 1 en 5 variables, $m$ vaut donc $5$

${𝔽}_2^5=\{z_0=(0,0,0,0,0),...,z_{31}=(1,1,1,1,1)\}$.

Les fonctions booléennes en 5 variables sont identifiées chacune à un mot binaire de longueur 32

$c_f=(f(z_0),...,f(z_{31}))$

L'ensemble des mots du code est

$\mathrm{R}\mathrm{M}(1,5)=\{c_f:\mathrm{\exists }a\in {𝔽}^6|f((X_0,...,X_4))=a_0{X}_0+a_1{X}_1+a_2{X}_2+a_3{X}_3+a_4{X}_4+a_5\}$

Ainsi le code est l'ensemble des tables de vérité des fonctions booléennes affines en 5 variables, la table de vérité de $f$ étant simplement le vecteur $c_f$.

Une autre construction

On décrit comment construire une matrice génératrice d'un code de longueur $n=2^m$. Posons :

${𝔽}_2^m=\{z_0,\dots ,z_{2^m-1}\}$

Dans l'espace ${𝔽}_2^n$ on définit, pour toute partie $A\subset {𝔽}_2^m$, les vecteurs ${𝕀}_A\in {𝔽}_2^n$ par

${\left({𝕀}_A\right)}_i=\{\begin{array}{cc}1& \text{ si }{z}_i\in A\\ 0& \text{ sinon}\end{array}$

On définit également dans ${𝔽}_2^n$ l'opération binaire :

$w\wedge y=(w_1\cdot y_1,\dots ,w_n\cdot y_n)$

appelé produit extérieur.

${𝔽}_2^m$ est un espace vectoriel à $m$ dimensions sur le corps ${𝔽}_2$, donc on écrit

$({𝔽}_2{)}^m=\{(x_0,\dots ,x_{m-1})|x_i\in {𝔽}_2\}$

Dans l'espace ${𝔽}_2^n$, on définit les vecteurs de longueur $n$ suivants : $v_0=(1,1,1,1,1,1,1,1)$ et

$v_i={𝕀}_{H_i}$

où $H_i$ sont des hyperplans dans $({𝔽}_2{)}^m$ (de dimension $n=2^m-1$) :

$H_i=\{x\in ({𝔽}_2{)}^m\mid x_i=0\}$

Le code de Reed-Muller $(r,m)$ d'ordre $r$ et longueur $n=2^m$ est le code généré par $v_0$ et le produit extérieur d'au plus $r$ des $v_i$ (par convention, s'il y a moins d'un vecteur, le produit extérieur est appliqué comme l'identité).

En fait, derrière cette construction se cache est celle donnée par les fonctions booléennes. En effet, si $f$ est une fonction booléenne en $m$ variables, on a

${𝕀}_{A_f}=c_f$ en posant $A_f=\{z\in {𝔽}_2^m:f(z)=1\}$

où $c_f$ est le vecteur défini à la section Code binaire. Il est alors facile de vérifier que

${𝕀}_{A_f}\wedge {𝕀}_{A_g}={𝕀}_{A_{f\cdot g}}=c_{f\cdot g}$

Pour terminer, il suffit de remarquer que $H_i=A_{e_i+1}$, où $e_i$ est la ième forme coordonnée, i.e :

$e_i(x_0,...,x_{m-1})=x_i$

Le vecteur $v_i$ est donc égal à $c_{e_i+1}$. La famille des produits extérieurs d'au plus $r$ vecteurs $v_i$, est donc constituée des vecteurs de la forme

$v_{i_0}\wedge ...\wedge v_{i_l}=c_{(e_{i_0}+1)...(e_{i_l}+1)}$ avec $l\le r$.

Bien entendu, les $(e_{i_0}+1)...(e_{i_l}+1)$ sont des fonctions booléennes en $m$ variables dont le degré est exactement $l$ et pour tous les $(i_0,...,i_l)$, $l\le r$, elles forment une famille génératrice des fonctions de degré au plus $r$.

On peut montrer que lorsque $H=A_f$ est un hyperplan, la fonction $f$ est affine. Il est donc suffisant de considérer des hyperplans engendrés par des fonctions linéairement indépendantes pour obtenir de la sorte les codes de Reed-Muller.

Paramètres

L'ensemble des fonctions booléennes de degré au plus $r$ est un espace vectoriel, ce code est donc linéaire. Par définition, sa longueur est $2^m$. La famille des monômes de degré au plus $r$ étant une base de cet espace, et ayant

$k=(\genfrac{}{}{0ex}{}{m}{0})+(\genfrac{}{}{0ex}{}{m}{1})+...+(\genfrac{}{}{0ex}{}{m}{r})$

éléments, sa dimension est $k$.

Le poids minimal du code $\mathrm{R}\mathrm{M}(r,m)$ est obtenu, entre autres, par les monômes de degré $r$ et vaut $2^{m-r}$.

L'ordre 1

Lorsque $r=1$, le code $\mathrm{R}\mathrm{M}(r,m)$ est constitué des fonctions affines, c'est-à-dire des fonctions de la forme

$f(x_0,...,x_{m-1})=a_0{x}_0+...+a_{m-1}{x}_{m-1}+a_m$

où les $a_i$ sont des éléments de ${𝔽}_2$.

La distance minimale du code de Reed-Muller d'ordre 1 est $2^{m-1}$, ce qui en fait un code optimal selon la borne de Griesmer: un code ayant au moins autant de mots et une distance minimale au moins aussi grande ne peut pas être plus court que $\mathrm{R}\mathrm{M}(1,m)$.

On peut même facilement être plus précis. En effet, le polynôme énumérateur des poids se calcule simplement :

$W(X,Y)=\sum _{c\in \mathrm{R}\mathrm{M}(1,m)}{X}^{w(c)}{Y}^{2^m-w(c)}=Y^{2^m}+(2^{m+1}-2)X^{2^{m-1}}{Y}^{2^{m-1}}+X^{2^m}$

où $w(c)$ désigne le poids de Hamming de $c$. Cela signifie qu'il existe 3 poids possibles pour les mots du code, soit $0$ pour le mot nul, soit $2^m$ pour le mot tout à $1$, soit $2^{m-1}$ pour tous les autres mots. Pour prouver ce résultat il suffit de considérer la forme d'une fonction affine donnée ci-dessus. Si les $a_i$ sont tous nuls pour $i\in \{1,...,m-1\}$, la seule valeur prise par $f$ est $a_m$ et on obtient le mot tout à zéro ou tout à un selon la valeur de $a_m$, donc respectivement les termes $Y^{2^m}$ et $X^{2^m}$. Si l'un des $a_i$ au moins est non nul pour $i\in \{1,...,m-1\}$, alors les éléments $x$ tels que $f(x)=0$ forment un espace affine de dimension $m-1$. Le cardinal de cet espace est donc $2^{m-1}$. On en déduit que l'ensemble des $x$ tels que $f(x)=1$ a pour cardinal $2^m-2^{m-1}=2^{m-1}$. Le mot associé à $f$ a donc un poids de $2^{m-1}$. Les $(2^{m+1}-2)$ fonctions pour lesquelles l'un des $a_i$ au moins est non nul donnent donc le terme $(2^{m+1}-2)X^{2^{m-1}}{Y}^{2^{m-1}}$.

Liens externes

Article de E.F. Assmus sur les codes de Reed et Muller, détaillant de nombreuses propriétés de ces codes.