Brevet FR2981526A1 — publié le 20130419

Procédé et dispositif de décodage de données codées à l'aide d'un code correcteur d'erreurs, procédé et dispositif de codage, et programme d'ordinateur correspondants. 1. Domaine de l'invention Le domaine de l'invention est celui du codage et du décodage de données, mettant en oeuvre un code correcteur d'erreurs. Plus précisément, l'invention propose une nouvelle technique de codage de données source, délivrant des paquets de données codées ou un flux de données codées comprenant les données source et des données de redondance, destiné(s) à être transmis sur un canal de transmission (par exemple de type hertzien, optique ou électrique), ou stocké(s) dans un support matériel. L'invention concerne également une technique de décodage correspondante, permettant notamment de corriger les erreurs de transmission inhérentes au canal de transmission. L'invention trouve notamment des applications dans les domaines suivants : la transmission d'information par télécommunications filaires électriques, comme dans les normes ADSL, ou optiques, sur fibres optiques ou en espace libre ; - la transmission d'information dans les communications radios spatiales et terrestres sans fil (« wireless » en anglais), comme dans les systèmes de télévision numérique TNT, de radio numérique DAB, de téléphonie GSM ou UMTS, de réseau radio WiFi, et aussi dans les futurs systèmes de communications, comme les futures normes DVB , 4G, LTE, « Internet du futur », ou entre véhicules, objets ou machines communicants ... ; - la compression et la décompression de source d'informations ; la génération et la détection de séquences dites d'embrouillage (« scrambling » en anglais) dans les systèmes CDMA ; le stockage d'informations dans des mémoires de masse magnétiques, optiques, mécaniques ou électriques pour constituer des disques durs, ou des mémoires vives d'ordinateurs, ou des clés mémoire à interface de type USB... ; la correction d'informations lors des calculs dans un circuit intégré d'un microprocesseur ou dans un ordinateur ; la robotique commandée par une intelligence artificielle à base de réseaux de neurone ; etc. 2. Art antérieur De nombreuses techniques de codage permettent la correction d'erreurs de transmission au décodage en générant des données de redondance à partir de données source. Ainsi, un code correcteur d'erreur est classiquement défini par : une longueur n, correspondant aux données en sortie du codeur (mot de code de longueur n formé de k données source et de (n - k) données de redondance), un nombre de bits ou de symboles d'information utiles k, correspondant aux données en entrée du codeur, encore appelées données source, et une distance minimale dmin. La distance minimale d'un code dmin correspond au minimum de distance entre deux mots de code. Elle permet de déterminer le nombre maximum d'erreurs que le code peut corriger dans un mot de code. Ainsi, plus la distance minimale dmin est grande, meilleur est le code correcteur d'erreurs, puisqu'il permet de détecter (dmin -1) symboles erronés et d'en corriger [(dmin -1) / 2] (où l'opérateur L.J désigne la partie entière). Parmi ces techniques de codage permettant la correction d'erreurs de transmission, les turbo-codes et les codes LDPC (en anglais « Low-Density Parity Check ») présentent de très bonnes performances en termes de correction d'erreurs pour des codes de grande longueur n, avec n de l'ordre de quelques milliers de bits au moins (n> 1000). En revanche, les turbo-codes et les codes LDPC présentent des performances plus faibles pour des codes de plus petite longueur n (n Ceci est en partie dû au fait que l'algorithme de propagation de probabilités utilisé au décodage (en anglais « Belief Propagation »), ou ses variantes, devient sous-optimal au fur et à mesure que la longueur minimale des cycles diminue On rappelle qu'un cycle correspond à un chemin fermé dans un graphe de Tanner représentant les contraintes que doivent remplir les symboles d'un mot pour être un mot du code. La notion de graphe de Tanner est bien connue, et notamment décrite dans l'article « Minimum-distance bounds by graph analysis » de R.M. Tanner, « IEEE Transactions on information theory », vol. 47, février 2001. Un tel algorithme de propagation de probabilités n'est donc pas efficace pour des codes ayant des matrices génératrices très denses (i.e. ayant beaucoup de bits à 1) et donc pour des codes ayant d'excellentes distances de Hamming minimum relatives tS = dmin / n . Par ailleurs, l'algorithme rigoureusement optimal de décodage exhaustif est difficile à réaliser car il est de complexité exponentielle en 0 (2k) , où k est la dimension du code, et donc très coûteux. De plus, il ne permet dans pas de décoder en temps réel les codes de petite longueur. De plus, pour un usage industriel, l'optimisation des codes correcteurs d'erreurs comprend la minimisation des complexités de codage et de décodage en termes de coûts de matériel et d'énergie consommée Or ces problèmes d'optimisation du décodage des codes correcteurs d'erreurs sont encore plus difficiles à résoudre quand les longueurs des codes sont petites à moyennes, c'est-à-dire de l'ordre de n Il existe donc un besoin pour une nouvelle technique de décodage des bons codes correcteurs d'erreurs notamment, c'est-à-dire des codes correcteurs d'erreurs présentant une distance minimale dmin la plus grande possible, et présentant notamment des cycles très courts. 3. Exposé de l'invention L'invention propose une solution nouvelle qui ne présente pas l'ensemble de ces inconvénients de l'art antérieur, sous la forme d'un procédé de décodage de données codées, comprenant des données source et des données de redondance. Ces données de redondance sont obtenues en appliquant, au codage, un code correcteur d'erreurs aux données source, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance P inversible. Selon l'invention, le procédé comprend une étape de décodage mettant simultanément en oeuvre le code correcteur d'erreurs, dit code correcteur d'erreurs primal, et un code dual du code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. L'invention propose ainsi une nouvelle technique de décodage de données codées à l'aide d'un code correcteur d'erreurs, utilisant en même temps le code correcteur d'erreurs sous sa forme primale et sous sa forme duale. En particulier, l'utilisation combinée du code correcteur d'erreurs sous ses deux formes primale et duale permet de diminuer le nombre d'itérations mises en oeuvre au décodage, notamment pour un algorithme de type propagation de probabilités (« Belief Propagation »), et donc de décoder en temps réel des codes de petite longueur (n L'invention offre ainsi un algorithme de décodage plus simple et plus performant que les techniques de l'art antérieur, notamment pour les codes présentant de petites longueurs et/ou des cycles courts. Selon un mode de réalisation particulier, l'étape de décodage met en oeuvre un graphe de Tanner modifié, dit graphe primal-dual, obtenu en superposant un graphe de Tanner correspondant au code correcteur d'erreurs primal, dit graphe primal, et un graphe de Tanner correspondant au code correcteur d'erreurs dual, dit graphe dual, la superposition consistant à connecter deux à deux les variables d'entrée du graphe primal aux variables d'entrée du graphe dual correspondant aux mêmes données source, et à connecter deux à deux les variables de sortie du graphe primal aux variables de sortie du graphe dual correspondant aux mêmes données de redondance. De cette façon, chaque variable d'entrée ou de sortie du graphe primal-dual est connectée à un couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2, définissant un arbre élémentaire. En d'autres termes, la structure de décodage proposée permet d'isoler de petits arbres élémentaires correspondant chacun à une donnée source ou une donnée de redondance. Du fait de cette structure en arbres élémentaires, le graphe primal-dual obtenu ne présente aucune structure de boucles, qui sont classiquement considérées comme source de dégradation sur le décodage. En d'autres termes, il existe un seul chemin d'une racine à une feuille de l'arbre, donc une seule façon de propager l'information dans le graphe primal-dual. Cette structure permet donc de diminuer les effets des cycles qui posent problème selon l'art antérieur, en utilisant de manière optimale l'information disponible. En d'autres termes, le traitement simultané des variables d'entrée/sortie dans le graphe primal-dual permet de compenser les calculs mis en oeuvre pour propager les variables et donc de diminuer l'influence des cycles sur l'algorithme de décodage. De plus, la propagation des variables internes se trouve modifiée par rapport à l'algorithme classique de type « Belief Propagation », et permet un décodage facilité des codes présentant une petite longueur. Cette structure en arbres élémentaires permet en effet d'obtenir une moyenne de la propagation a priori de données source dans chaque arbre. Finalement, la présence d'autant de couples de contraintes (ou d'arbres élémentaires) que de données codées (données source et données de redondance) impose d'utiliser toutes les données disponibles, ce qui accroît les performances du décodage. On dispose également de deux fois plus d'informations au niveau des variables d'entrée et de sortie du graphe primal-dual que selon les techniques de l'art antérieur. La superposition des graphes du code primal et du code dual permet donc d'améliorer les performances de décodage par rapport aux techniques actuelles, en permettant notamment une convergence plus rapide de l'algorithme de décodage vers un mot de code le plus proche du mot reçu. Elle permet également de compenser la non linéarité liée aux contraintes (égalité, addition modulo 2). En particulier, chaque variable d'entrée ou de sortie du graphe primal-dual étant connectée à un couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2, le procédé de décodage selon l'invention comprend une étape de détermination d'au moins une variable interne au graphe primal-dual mettant en oeuvre : pour une contrainte de type addition modulo 2 : o une étape de transformation, du domaine réel vers le domaine complexe, des variables connectées à la contrainte de type addition modulo 2 ; o une étape d'addition dans le domaine complexe ; o une étape de transformation, du domaine complexe vers le domaine réel, du résultat de l'étape d'addition ; pour une contrainte de type égalité : o une étape d'addition dans le domaine réel. On simplifie de cette façon les opérations mises en oeuvre pour propager les variables connectées à une contrainte de type égalité ou de type addition modulo 2, en utilisant des opérateurs de type addition dans le domaine complexe ou dans le domaine réel. En particulier, les étapes de transformation (du domaine réel vers le domaine complexe ou du domaine complexe vers le domaine réel) mettent en oeuvre une fonction non linéaire NL telle que : INL[x]= log[1- exp[x]] si x # 0 1+ exp[x] NL[x] = 0 si x = 0 r On note que cette fonction non linéaire est égale à sa propre inverse : NL[x] = NL 1 Lx _1. L'application d'une telle fonction non linéaire permet de réaliser l'équivalent du calcul des probabilités a posteriori en sortie d'une contrainte de type addition modulo 2 (ou exclusif) au moyen d'une simple addition dans le corps des réels et d'une multiplication des signes. On diminue de cette façon la complexité opératoire du décodage. Selon un mode de réalisation particulier, chaque variable d'entrée ou de sortie du graphe primal-dual étant connectée à un couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2, le procédé de décodage met en oeuvre au moins une itération de décodage, comprenant : une phase de propagation avant, consistant à propager les variables internes au graphe primal-dual connectées à une contrainte de type égalité pour mettre à jour les variables internes au graphe primal-dual connectées à une contrainte de type addition modulo 2, et une phase de propagation arrière, consistant à propager les variables internes au graphe primal-dual connectées à une contrainte de type addition modulo 2 pour mettre à jour les variables internes au graphe primal-dual connectées à une contrainte de type égalité. En d'autres termes, une itération de décodage comprend une première phase de propagation des données source dans tous les arbres élémentaires, puis une deuxième phase de propagation des données de redondance dans tous les arbres élémentaires, chaque phase se transmettant les valeurs obtenues en sortie de chaque arbre élémentaire. Les opérateurs associés à chaque couple de contraintes travaillent donc simultanément. En particulier, la phase de propagation avant met en oeuvre les sous-étapes suivantes, pour chaque couple de contraintes connecté à une variable d'entrée ou de sortie du graphe primal-dual : 30F détermination d'une valeur de vraisemblance avant vo (t) affectée à la variable d'entrée ou de sortie du graphe primal-dual, mettant en oeuvre une addition dans le domaine réel, telle que : vo F (t)= iivi(t) d détermination de valeurs de vraisemblance affectées aux variables connectées à la contrainte de type addition modulo 2, après transformation du domaine réel vers le domaine complexe, dites valeurs de vraisemblance avant transformées, telle que : o la variable d'entrée ou de sortie du graphe primal-dual après transformation s'exprime sous la forme : 14(t) = NL[vo(0)+ vp (01 et so (t) = signe[vo(0)+ vo (H; o les variables internes après transformation s'expriment sous la forme : wr (t) = NL[vii (01 et sr (t) = signe[vii (H pour 1 E ; mise à jour des valeurs de vraisemblance avant transformées, mettant en oeuvre une addition dans le domaine complexe, telle que : w1 d. (t +1), wt(o+ (t) et sF (t +1) =ri.d s- (t) d -1 1 1-0,11 mise à jour des valeurs des variables internes connectées à la contrainte de type addition modulo 2, mettant en oeuvre une transformation du domaine complexe vers le domaine réel des valeurs de vraisemblance avant transformées mises à jour, telle que : vi(t +1) = sr (t + 1)NL[wr (t +1)1; avec : v0(0) une valeur initiale affectée à la variable d'entrée ou de sortie du graphe primal-dual ; vi(t) les variables internes connectées à la contrainte de type égalité ; vi (t) les variables internes connectées à la contrainte de type addition modulo 2 ; d le degré de contraintes associé à une contrainte du graphe primal-dual ; i / e ; NL[x] une fonction non linéaire telle que : NL[x]= log[1- exp[x]] si x # 0 1+ exp[x] NL[x] = 0 si x = 0 La phase de propagation arrière met quant à elle en oeuvre les sous-étapes suivantes, pour chaque couple de contraintes connecté à une variable d'entrée ou de sortie du graphe primal-dual : - détermination de valeurs de vraisemblance affectées aux variables connectées à la contrainte de type addition modulo 2, après transformation du domaine réel vers le domaine complexe, dites valeurs de vraisemblance arrière transformées, telle que : o les variables internes après transformation s'expriment sous la forme : wir3(t) = NL[v;(t)] et 43(0 = signe[vii(t)], pour / ; o la variable d'entrée ou de sortie du graphe primal-dual après transformation iwiB(t) s'exprime sous la forme : (t) = et so (t) = d siB (0) détermination d'une valeur de vraisemblance arrière vo (t) affectée à la variable d'entrée ou de sortie du graphe primal-dual, mettant en oeuvre une transformation du domaine complexe vers le domaine réel de la variable d'entrée ou de sortie après transformation, telle que : 1(0 = sÔ X NL[we (t)] mise à jour des valeurs des variables internes connectées à la contrainte de type égalité, mettant en oeuvre une addition dans le domaine réel, telle que : d v- (t) y/ (t + 1) = vo (0) + vic? (t) + d - 1 mise à jour de la variable d'entrée ou de sortie du graphe primal-dual, telle que : vo (t + 1) = vp (t) + vo (0) + ve (t) . Selon une caractéristique particulière de l'invention, le procédé de décodage comprend une étape préalable d'initialisation des variables d'entrée, des variables de sortie, et des variables internes du graphe primal-dual à zéro. La mise en oeuvre de ces deux phases de propagation avant et arrière selon l'invention, ainsi que l'initialisation de l'ensemble des variables, permet d'obtenir une convergence très rapide de l'algorithme de décodage vers le mot de code le plus proche du mot reçu. Selon un autre aspect de l'invention, au moins deux variables d'entrée ou de sortie du graphe primal-dual sont connectées à un même couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2. Il est ainsi possible de multiplexer les opérateurs utilisés pour propager les variables connectées à une contrainte de type égalité ou de type addition modulo 2, et de diminuer ainsi la quantité de ressources matérielles nécessaire. Selon un mode de réalisation particulier de l'invention, le code correcteur d'erreurs primal présente un rendement égal à '/2. En d'autres termes, il s'agit d'un code quasi-aléatoire de paramètres (n = 2k, k, dmin) . On obtient dans ce cas une structure symétrique du graphe primal- dual, permettant de décoder simplement et de façon optimale un tel code de rendement égal à 1/2. Bien entendu, il est également possible de décoder des codes correcteurs d'erreurs mettant en oeuvre un rendement différent de '/2. Ainsi, pour décoder un code correcteur d'erreurs présentant un rendement supérieur à '/2, on peut ajouter une contrainte égalité (respectivement une contrainte addition modulo 2) en sortie de deux contraintes égalité (respectivement deux contraintes addition modulo 2) connectées à deux (ou plus) données de redondance dans le graphe primal-dual. Pour décoder un code correcteur d'erreurs présentant un rendement inférieur à '/2 on peut ajouter une contrainte égalité (respectivement une contrainte addition modulo 2) en sortie de deux contraintes égalité (respectivement deux contraintes addition modulo 2) connectées à deux (ou plus) données source dans le graphe primal-dual. Dans un autre mode de réalisation, l'invention concerne un dispositif de décodage de données codées correspondant. Selon l'invention, un tel dispositif comprend un module de décodage mettant simultanément en oeuvre le code correcteur d'erreurs, dit code correcteur d'erreurs primal, et un code dual du code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. Un tel dispositif de décodage est notamment adapté à mettre en oeuvre le procédé de décodage décrit précédemment. Un tel décodeur pourra bien sûr comporter les différentes caractéristiques relatives au procédé de décodage selon l'invention, qui peuvent être combinées ou prises isolément. Ainsi, les caractéristiques et avantages de ce décodeur sont les mêmes que ceux du procédé de décodage, et ne sont donc pas détaillés plus amplement. On note qu'un tel décodeur peut être implémenté sous la forme d'un circuit intégré, ou dans un composant électronique de type microprocesseur. Ainsi, l'algorithme de décodage selon l'invention peut être mis en oeuvre de diverses manières, notamment sous forme câblée ou sous forme logicielle. En particulier, il est possible d'utiliser un décodeur tel que décrit ci-dessus soit pour décoder des données codées, soit pour coder des données source. Dans ce dernier cas, on peut forcer à zéro les données de redondance, afin de pouvoir utiliser le décodeur comme un codeur. L'invention concerne ainsi également un procédé de codage de données source, délivrant des données codées comprenant les données source et des données de redondance. Selon l'invention, un tel procédé de codage comprend une étape de codage mettant simultanément en oeuvre : - un code correcteur d'erreurs, dit code correcteur d'erreurs primal, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance P inversible, et - un code dual du code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. L'invention propose ainsi une nouvelle technique de codage de données source à l'aide d'un code correcteur d'erreurs utilisant en même temps le code correcteur d'erreurs sous sa forme primale et sous sa forme duale. Il est ainsi possible d'obtenir simplement de nouveaux codes ou de reconstruire des codes existants, présentant une longueur quelconque, facilement décodables par le procédé de décodage selon l'invention. Un tel procédé de codage pourra bien sûr comporter les différentes caractéristiques relatives au procédé de décodage selon l'invention. Dans un autre mode de réalisation, l'invention concerne un dispositif de codage correspondant. Un tel dispositif comprend un module de codage mettant simultanément en oeuvre : - un code correcteur d'erreurs, dit code correcteur d'erreurs primal, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance P inversible, et - un code dual du code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. Un tel dispositif de codage est notamment adapté à mettre en oeuvre le procédé de codage décrit précédemment. A nouveau, un tel codeur peut être implémenté sous la forme d'un circuit intégré, ou dans un composant électronique de type microprocesseur. Ainsi, l'algorithme de codage selon l'invention peut être mis en oeuvre de diverses manières, notamment sous forme câblée ou sous forme logicielle. L'invention concerne encore un dispositif de codage et de décodage de données, comprenant un module de codage et de décodage mettant simultanément en oeuvre : - un code correcteur d'erreurs, dit code correcteur d'erreurs primal, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance P inversible, et - un code dual du code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. Un tel dispositif, encore appelé codec pour codeur/décodeur, peut être implémenté sous la forme d'un circuit intégré, ou dans un composant électronique de type microprocesseur. En particulier, l'utilisation d'un codec permet de mutualiser les ressources matérielles utilisées au codage et au décodage. Par exemple, un tel codec peut recevoir en entrée des données source, et délivrer en sortie des données codées, ou recevoir en entrée des données codées, et délivrer en sortie une estimation des données source. Dans encore un autre mode de réalisation, l'invention concerne un ou plusieurs programmes d'ordinateur comportant des instructions pour la mise en oeuvre d'un procédé de décodage ou d'un procédé de codage tels que décrits précédemment, lorsque le ou les programmes sont exécutés par un processeur. De tels programmes peuvent être stockés sur un support d'information. 4. Liste des figures D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement à la lecture de la description suivante d'un mode de réalisation particulier, donné à titre de simple exemple illustratif et non limitatif, et des dessins annexés, parmi lesquels : la figure 1 présente les principales étapes mises en oeuvre par un codeur/décodeur selon l'invention ; les figures 2A, 2B et 2C illustrent respectivement un graphe de Tanner d'un code de Hamming (8, 4, 4) présentant une matrice de passage P inversible, un graphe de Tanner du code dual du code de Hamming (8, 4, 4), et un graphe primal-dual correspondant selon l'invention ; les figures 3A, 3B et 3C illustrent respectivement un graphe de Tanner d'un code (12,6) présentant une matrice de passage P inversible, un graphe de Tanner du code dual de ce code (12,6), et un graphe primal-dual corrspondant selon l'invention ; la figure 4 détaille un couple de contraintes selon l'invention ; - les figures 5A et 5B illustrent la fonction de non linéarité NL utilisée pour déterminer les variables internes à un graphe primal-dual selon l'invention ; les figures 6A et 6B présentent les structures simplifiées d'un codeur et d'un décodeur selon un mode de réalisation particulier de l'invention. 5. Description d'un mode de réalisation de l'invention 5.1 Principe général Le principe général de l'invention repose sur l'utilisation simultanée d'un code correcteur d'erreurs, dit code primal, et du dual de ce code correcteur d'erreurs, dit code dual, pour coder ou décoder des données. Elle offre notamment une solution au problème du décodage à décision douce des codes en bloc correcteurs d'erreurs de longueur n contenant k données source (encore appelée symboles ou bits d'information). En particulier, elle propose une nouvelle technique de décodage optimal des codes correcteurs d'erreurs présentant de grandes distances minimales pour ces paramètres n et k. A titre d'exemple, l'invention offre une technique de décodage présentant de meilleures performances que les techniques de l'art antérieur vis-à-vis des bons codes correcteurs d'erreurs tels que : les codes algébriques classiques comme le code de Golay (24,12,8), le code (48,24,12), les codes de Reed-Müller (RM) ; les turbo-codes ; les codes LDPC, etc. On présente, en relation avec la figure 1, les principales étapes mises en oeuvre par un codeur et/ou un décodeur binaire selon l'invention. Il s'agit ici d'un simple exemple, dans lequel on utilise l'invention à la fois pour le codage et pour le décodage des données. Bien entendu, il est possible d'utiliser l'algorithme de codage proposé avec un algorithme de décodage classique, ou d'utiliser l'algorithme de décodage proposé avec un algorithme de codage classique. Dans ce dernier cas, la seule contrainte imposée au codeur est que la matrice de passage des données source aux données de redondance soit inversible. Ainsi, selon le mode de réalisation illustré en figure 1, un codeur selon l'invention prend en entrée des données source S, sous la forme de paquets ou d'un flux continu, et délivre en sortie des données codées, comprenant les données source S et des données de redondance R. Ces données codées peuvent être véhiculées dans un canal de transmission ou stockées sur un support, puis décodées par un décodeur selon l'invention. Un tel décodeur reçoit donc en entrée les données codées, comprenant les données source S et des données de redondance R éventuellement entachées d'erreur, et délivre une estimation des /N. données source S. Plus précisément, le codeur met en oeuvre une étape de codage 11 mettant simultanément en oeuvre : - un code correcteur d'erreurs, dit code correcteur d'erreurs primal CP, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance P inversible, et - un code dual du code correcteur d'erreurs primal, dit code correcteur d'erreurs dual CD. De la même façon, le décodeur met en oeuvre une étape de décodage 12 mettant simultanément en oeuvre le code correcteur d'erreurs primal CP et le code correcteur d'erreurs dual CD. On note que l'étape de codage 11 et l'étape de décodage 12 utilisent toutes les deux à la fois le code correcteur d'erreurs primal CP et le code correcteur d'erreurs dual CD. Le codeur et le décodeur peuvent donc partager des ressources, ce qui permet de diminuer les complexités de codage et de décodage en termes de coûts de matériel et d'énergie consommée 5.2 Graphe de Tanner primal-dual Selon l'invention, l'étape de codage 11 et/ou l'étape de décodage 12 mettent en oeuvre un graphe de Tanner modifié, dit graphe primal-dual, obtenu en superposant un graphe de Tanner correspondant au code primal CP, dit graphe primal, et un graphe de Tanner correspondant au code dual CD, dit graphe dual. On décrit ci-après deux exemples de construction de tels graphes primal-dual, pour un code de Hamming étendu (8, 4, 4) et pour un code (12,6) présentant une matrice de passage des données source aux données de redondance P inversible. Ces deux codes présentent un rendement égal à '/2, mais comme déjà indiqué, il est également possible de décoder des codes mettant en oeuvre un rendement supérieur à 1/2 ou inférieur à 1/2, respectivement en : - ajoutant une contrainte égalité (respectivement une contrainte addition modulo 2) en sortie de deux contraintes égalité (respectivement deux contraintes addition modulo 2) connectées à deux (ou plus) données de redondance dans le graphe primal-dual, ce qui revient à considérer un code correcteur d'erreurs de départ quasi-aléatoire de matrice de passage P inversible et à ajouter des colonnes à la matrice de passage P ; ajoutant une contrainte égalité (respectivement une contrainte addition modulo 2) en sortie de deux contraintes égalité (respectivement deux contraintes addition modulo 2) connectées à deux (ou plus) données source dans le graphe primal-dual, ce qui revient à considérer un code correcteur d'erreurs de départ quasi-aléatoire de matrice de passage P inversible et à ajouter des lignes à la matrice de passage P. Bien entendu, l'invention peut être utilisée pour décoder d'autres codes, comme le code de Golay, les codes de Reed-Müller, les turbo-codes, les codes LDPC, etc. A) Premier exemple : code de Hamming étendu (8, 4, 4) On présente ci-après la construction d'un graphe primal-dual, obtenu à partir des graphes primal et dual du code de Hamming étendu (n = 8, k = 4, dmin = 4) . Un tel code met en oeuvre une matrice génératrice G comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance P inversible, telles que par exemple : G= [Id IP] = 1 0 oloo 0 0 0 0 0 0 1 0 0 1 Il s'agit ici d'un cas particulier, selon lequel la matrice de passage P est égale à sa transposée P=Ptr , et donc PPtr = Id , d'où la matrice de passage P est égale à sa propre inverse P-1=Ptr . Le code de Hamming (8, 4, 4) selon cet exemple est donc auto-dual. On considère selon cet exemple des données source prenant la forme de quatre bits d'information notés (x0 , x1, x2 , x3) , et des données de redondance prenant la forme de quatre bits de redondance notés (ro , rl , r2 , r3 ) De manière classique, les relations algébriques de calcul des bits de redondance à partir des bits d'information peuvent être représentées par un graphe de Tanner correspondant au code de Hamming (8, 4, 4) primal, tel qu'illustré en figure 2A. Un tel graphe, dit graphe primal, comprend des variables d'entrée, correspondant aux bits d'information, et des variables de sortie, correspondant aux bits de redondance, les bits de redondance étant liés aux bits d'informations par l'intermédiaire de contraintes et de permutations. Ces contraintes peuvent être de deux types : les contraintes de type égalité, représentées par le signe «_ », correspondant à une répétition de la variable connectée à cette contrainte ; les contraintes de type addition modulo 2, représentées par le signe « + », correspondant à une contrainte logique de type « ou-exclusif». En d'autres termes, ces contraintes et permutations permettent de représenter les relations algébriques suivantes : r0 = xi + x2 + x3 ; ri = xo + x2 + x3 ; r2 = xo + xi + x3 et r3 = xo + xi + x2. Comme le code de Hamming (8, 4, 4) est auto-dual, sa matrice de passage P est égale à sa propre inverse, et il suffit de remplacer les contraintes de type égalité par des contraintes de type addition modulo 2, et les contraintes de type addition modulo 2 par des contraintes de type égalité, pour obtenir le graphe de Tanner correspondant au code dual du code de Hamming (8, 4, 4). Les permutations mises en oeuvre sont donc identiques entre les deux graphes primal et dual. Au niveau des ressources matérielles, on utilise donc selon cet exemple deux entrelaceurs identiques. La figure 2B illustre ainsi le graphe dual obtenu, comprenant des variables d'entrée, correspondant aux bits d'information, et des variables de sortie, correspondant aux bits de redondance. Selon l'invention, on définit un graphe de Tanner modifié, dit graphe primal-dual, permettant de tenir simultanément compte du code primal et du code dual. Pour ce faire, on superpose les graphes primal, illustré en figure 2A, et dual, illustré en figure 2B, pour obtenir un graphe primal-dual, illustré en figure 2C. Cette superposition consiste à connecter deux à deux les variables d'entrée du graphe primal (x0 , xi, x2 , x3) aux variables d'entrée du graphe dual (x0 , x1, x2 , x3) correspondant aux mêmes données source, et à connecter deux à deux les variables de sortie (ro , ri, r2 , r3) du graphe primal aux variables de sortie (r0 , r1, r2 , r3) du graphe dual correspondant aux mêmes données de redondance. En d'autres termes, on relie les bits d'information xi du graphe primal aux bits d'information xi du graphe dual, pour i entier allant de 0 à 3, et on relie les bits de redondance ri du graphe primal aux bits de redondance ri du graphe dual, pour i entier allant de 0 à 3. Chaque variable d'entrée xi du graphe primal-dual et chaque variable de sortie ri se trouve donc connectée à un couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2. Selon l'exemple illustré en figure 2C, on obtient donc huit couples de contraintes, encadrés sur la figure, connectés chacun à : - une variable d'entrée ou de sortie (de type bit d'information xi ou bit de redondance ri ) ; et I I I six variables internes (notées par exemple v1, v2,123, v1, v2 , v3 pour le couple de contraintes 210, où les variables notées vk proviennent d'une contrainte de type égalité et les variables notées vk proviennent d'une contrainte de type addition modulo 2). Les quatre couples de contraintes 210, 211, 212, 213 connectés aux variables d'entrée (xo , x1, x2 , x3) sont reliés aux quatre couples de contraintes 220, 221, 222, 223 connectés aux variables de sortie (r0 , rl, r2 , r3) par l'intermédiaire de permutations ou entrelaceurs de 4 x 3 = 12 couples d'état interne au codeur-décodeur. Chaque couple de contraintes et ses connexions peut être considéré comme un arbre de la théorie des graphes, c'est-à-dire qui ne présente aucune structure de boucle. B) Deuxième exemple : code (12,6) On présente ci-après un deuxième exemple de construction d'un graphe primal-dual selon l'invention, obtenu à partir des graphes primal et dual du code (12,6). On se place ici dans un cas plus général, où la matrice de passage P est simplement inversible. Un tel code met en oeuvre une matrice génératrice G comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance P inversible, telles que par exemple : 1 0 0 0 0 0 11 0 1 0 0 1 0 1 0 0 0 0 11 0 1 0 1 0 o o 1 0 0 0 Il 1 0 0 0 1 o 0 0 1 0 0 10 1 0 1 0 0 o 0 0 0 1 0 Il 0 1 1 0 1 o 0 0 0 0 1 10 0 1 0 1 0 On considère selon cet exemple des données source prenant la forme de six bits d'information notés (x0 , xl, x2 , x3, x4 , x5), et des données de redondance prenant la forme de six bits de redondance notés (r0 , r1, r2 , r3, r4 , r5) . De manière classique, les relations algébriques de calcul des bits de redondance à partir des 25 bits d'information peuvent être représentées par un graphe de Tanner correspondant au code (12,6) primal, tel qu'illustré en figure 3A. Comme indiqué en relation avec le premier exemple, un tel graphe primal comprend des variables d'entrée, correspondant aux bits d'information, et des variables de sortie, correspondant aux bits de redondance, les bits de redondance étant liés aux bits d'informations par l'intermédiaire G=[Id1P]= de contraintes (de type égalité ou addition modulo 2) et de permutations. Considérant que la matrice de passage P est inversible, on peut déterminer la matrice de passage inverse P-1 et le graphe de Tanner correspondant au dual du code (12,6). La figure 3B illustre ainsi le graphe dual obtenu, comprenant des variables d'entrée, correspondant aux bits d'information, et des variables de sortie, correspondant aux bits de redondance. On note que comme le code (12,6) n'est pas auto-dual, sa matrice de passage P n'est pas égale à sa propre inverse, et les graphes primal et dual obtenus ne sont pas symétriques. Au niveau des ressources matérielles, on utilise donc selon cet exemple deux entrelaceurs distincts. Selon l'invention, et comme décrit en relation avec le premier exemple, on définit un graphe de Tanner modifié, dit graphe primal-dual, permettant de tenir simultanément compte du code primal et du code dual. Pour ce faire, on superpose les graphes primal, illustré en figure 3A, et dual, illustré en figure 3B, pour obtenir un graphe primal-dual, illustré en figure 3C. Cette superposition consiste à connecter deux à deux les variables d'entrée du graphe primal (xo , xi , x2 , x3, x4 , x5) aux variables d'entrée du graphe dual (xo , xi , x2, x3, x4 , x5 ) correspondant aux mêmes données source, et à connecter deux à deux les variables de sortie (ro , ri ,r2,r3,r4,r5) du graphe primal aux variables de sortie (ro ,ri,r2,r3,r4,r5) du graphe dual correspondant aux mêmes données de redondance. En d'autres termes, on relie les bits d'information xi du graphe primal aux bits d'information xi du graphe dual, pour i entier allant de 0 à 5, et on relie les bits de redondance ri du graphe primal aux bits de redondance ri du graphe dual, pour i entier allant de 0 à 5. Chaque variable d'entrée xi du graphe primal-dual et chaque variable de sortie ri se trouve donc connectée à un couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2. Selon l'exemple illustré en figure 3C, on obtient donc douze couples de contraintes, encadrés sur la figure, connectés chacun à : - une variable d'entrée ou de sortie (de type bit d'information xi ou bit de redondance ri ) ; et - un ensemble de variables internes. Les six couples de contraintes 310, 311, 312, 313, 314, 315 connectés aux variables d'entrée (x0 , xl, x2 , x3, x4 , x5) sont reliés aux six couples de contraintes 320, 321, 322, 323, 324, 325 connectés aux variables de sortie (ro , ri , r2 ,r3,r4 , r5) par l'intermédiaire de permutations ou entralaceurs. A nouveau, chaque couple de contraintes et ses connexions peut être considéré comme un arbre de la théorie des graphes, qui ne présente aucune structure de boucle. 5.3 Détermination des variables connectées à des contraintes de type égalité ou addition modulo 2 On décrit ci-après les étapes mises en oeuvre pour la détermination des variables internes à un graphe primal-dual, associées à un couple de contraintes, selon un aspect particulier de l'invention. Bien entendu, d'autres techniques, et notamment des techniques classiques, peuvent être utilisées pour déterminer ces variables internes. On considère que les variables d'entrée ou de sortie connectées à un couple de contraintes, notées vi (0) à l'instant initial t = 0 , sont homogènes à des logarithmes de rapport de vraisemblance, classiquement définis par : [proba(yilci =1) vi (0) = log proba (yi Ici = 0) Les variables internes sont donc aussi homogènes à des logarithmes de rapport de vraisemblance. On considère également, à titre d'exemple, un couple de contraintes du graphe primal-dual du code de Hamming (8, 4, 4), connecté d'une part à une variable d'entrée ou de sortie (de type bit d'information xi ou bit de redondance ri , noté plus généralement vi ) et d'autre part à six variables internes. Un tel exemple de couple de contraintes est illustré en figure 4. Pour effectuer la fusion des variables connectées à une contrainte de type égalité d'un couple de contraintes 41 (fusion des variables d'entrée ou de sortie à l'instant initial et des variables internes), on peut ajouter à la variable d'entrée ou de sortie à l'instant initial la moyenne des valeurs des variables internes connectées à cette contrainte. Ainsi, si l'on considère une variable d'entrée ou de sortie vo du couple de contraintes 41 à l'instant initial t = 0, la nouvelle valeur de cette variable à l'instant suivant t =1 (c'est-à-dire à l'itération suivante) peut s'exprimer sous la forme suivante : v0(1) = v0(0)+ (vi(0)+ v2(0)+ v3(0)) 3 On effectue donc une addition dans le domaine réel pour déterminer les valeurs des variables connectées à une contrainte de type égalité Pour effectuer la fusion des variables connectées à une contrainte de type addition modulo 2 d'un couple de contraintes 41 (fusion des variables d'entrée ou de sortie à l'instant initial et des variables internes), on propose selon un aspect particulier de l'invention de travailler dans le domaine complexe plutôt que dans le domaine réel, afin de pouvoir utiliser une simple opération d'addition. Ainsi, selon cet aspect de l'invention, on applique une fonction non-linéaire NL H à toutes les variables connectées à une contrainte de type addition modulo 2, avant de pouvoir les additionner, dans le domaine complexe. Plus précisément, on applique la fonction NL à la variable d'entrée ou de sortie vo du couple de contraintes 41 à l'instant t = 0, dans le domaine réel, pour obtenir une variable interne dans le domaine complexe notée wo , et la fonction NL aux trois variables internes v1, v2 , v3 dans le domaine réel, pour obtenir trois variables internes w1, w2 , w3 dans le domaine complexe. On repasse ensuite dans le domaine réel pour obtenir la nouvelle valeur de la variable vo à l'instant suivant t = 1. Ainsi, on obtient : (Ndvi(0)1+ NL[v21 (0)1+ NL[v31 (0)1) Ndvo(0)1+ 3 v0(1) = NL - NL[wo +(wl + w2 +1423) 3 On effectue donc une transformation du domaine réel vers le domaine complexe des variables connectées à la contrainte de type addition modulo 2, une addition dans le domaine complexe, puis une transformation du domaine complexe vers le domaine réel du résultat de l'addition, pour déterminer les valeurs des variables connectées à une contrainte de type addition modulo 2. En particulier, la fonction non-linéaire NL peut s'exprimer N[1- exp[x] sous la forme suivante : L[x]= log si x # 0 ] si x = 0 1+ exp[x] NL[x] = 0 Les figures 5A et 5B illustrent respectivement la partie réelle de la fonction NL /eSENL[x]] (9i [NL [x]1) et la variation du signe de la partie imaginaire de la fonction NL Plus précisément, étant donné le domaine de variation de x E 1-00,+00[ , le résultat de l'application de la fonction non-linéaire NL est rigoureusement à valeurs dans le corps des nombres complexes C, avec une partie imaginaire pouvant prendre deux valeurs {0,+/7L} , où I est le nombre imaginaire pur. De plus, aux fins d'initialisation, on impose NL[0]= 0 . Pour simplifier les calculs, il est alors possible de considérer cette partie imaginaire pure comme un signe, et de multiplier simplement le résultat de l'addition des parties réelles par le produit des signes. Cette addition dans le « domaine dual » est donc équivalente à une addition dans le corps des complexes C , ce qui est particulièrement astucieux et novateur, notamment pour les décodeurs à décision douce. Cette fonction non-linéaire NL peut être réalisée avec des nombres quantifiés de façon très peu complexe avec une mémoire programmable par exemple, comme une EPROM (en anglais « Erasable Programmable Read-Only Memory »), également appelée LUT (en anglais « Look-Up Table »). Bien entendu, des opérations similaires peuvent être mises en oeuvre pour déterminer les valeurs des variables connectées à une contrainte de type égalité ou à une contrainte de type addition modulo 2 dans d'autres contextes, et notamment pour déterminer les valeurs des variables non connectées à un couple de contraintes. De même, des opérations similaires peuvent être mises quelque soit le nombre de variables internes connectées à la contrainte, en adaptant les équations ci-dessus. 5.4 Propagation des variables On présente ci-après les différentes étapes mises en oeuvre pour décoder des données codées reçues, en utilisant simultanément un code correcteur d'erreurs primal et son dual. On considère pour construire un tel décodeur les graphes primal-dual tels que définis précédemment, et notamment les exemples illustrés en figures 2C et 3C. Plus précisément, l'algorithme de décodage met en oeuvre au moins une itération de décodage, comprenant : une phase de propagation avant, consistant à propager les variables internes au graphe primal-dual connectées à une contrainte de type égalité pour mettre à jour les variables internes au graphe primal-dual connectées à une contrainte de type addition modulo 2, et une phase de propagation arrière, consistant à propager les variables internes au graphe primal-dual connectées à une contrainte de type addition modulo 2 pour mettre à jour les variables internes au graphe primal-dual connectées à une contrainte de type égalité. Une itération met en oeuvre les deux phases de propagation avant et arrière pour chaque couple de contraintes (soit huit couples pour le code de Hamming (8, 4, 4), douze couples pour le code (12,6), vingt-quatre couples pour un code de Golay (24, 12, 8), etc). A l'initialisation, les variables d'entrée, de sortie, et internes au graphe primal-dual sont mises à zéro. Chaque arbre reçoit alors, généralement en provenance d'un démodulateur, des données codées, correspondant soit à une donnée source (bit d'information) soit à une donnée de redondance (bit de redondance). On décrit ci-après les phases de propagation avant et arrière mises en oeuvre pour chaque arbre élémentaire associé à un couple de contraintes, au cours d'une même itération. On note vo (0) la valeur initiale affectée à une variable d'entrée ou de sortie du graphe primal-dual connectée à un couple de contraintes considéré. La phase de propagation avant, encore appelée « forward propagation », consiste à propager les vraisemblances y/ (t) du domaine primal pour calculer les vraisemblances v; (t + 1) du domaine dual, où t est l'indice de temps discret, correspondant encore au numéro de l'itération en cours. Plus précisément, au cours d'une première étape, on détermine une valeur de vraisemblance avant vo (t) affectée à la variable d'entrée ou de sortie connectée à la contrainte d Li=lvi(t) considérée, mettant en oeuvre une addition dans le domaine réel, telle que vo (t) = d avec i G {1,...,4, d le degré de contraintes associé à la contrainte égalité du couple de contraintes considéré, et vi (t) les variables internes connectées à la contrainte égalité du couple de contraintes considéré. Au cours d'une deuxième étape, on détermine les valeurs de vraisemblance affectées aux variables connectées à la contrainte de type addition modulo 2, après transformation du domaine réel vers le domaine complexe, dites valeurs de vraisemblance avant transformées, en appliquant la fonction de non linéarité NL comme détaillé au paragraphe ci-dessus. Ainsi, la variable d'entrée ou de sortie du couple de contraintes considéré après transformation, s'exprime sous la forme wr) (t) = NL [Vo (0) + yr) (01 et 5'0 (t) = signe[vo(0)+ vo (01, et les variables internes, après transformation, s'expriment sous la forme wir (t) = NL[vi (01 et si (t) = signe[vi (01, avec / e {1,...,d}, d le degré de contraintes associé à la contrainte addition modulo 2 du couple de contraintes considéré, et v; (t) les variables internes connectées à la contrainte addition modulo 2 du couple de contraintes considéré. On peut remarquer que les valeurs wl (t) correspondent aux parties réelles de nombres complexes dont les parties imaginaires peuvent prendre les valeurs {0,+/z}, mais en prenant l'exponentielle de ces valeurs, on obtient des valeurs de signe si (t) E 1+1,-11, ce qui permet de simplifier les calculs. Au cours d'une troisième étape, on détermine les nouvelles valeurs wi (t + 1) et s (t + 1) en mettant à jour les valeurs de vraisemblance avant transformées, à partir d'une addition dans le domaine complexe, telle que : w (t + 1) = wô (t) + (t) et si (t + 1) = ric/ sr(t). d -1 Au cours d'une quatrième étape, on applique de nouveau la fonction de non linéarité NL, afin de mettre à jour les valeurs des variables internes connectées à la contrainte de type addition modulo 2 viir + 1), à partir d'une transformation du domaine complexe vers le domaine réel des valeurs de vraisemblance avant transformées mises à jour, telle que (t +1) = sr (t +1)NL[wr (t +1)1. La phase de propagation arrière, encore appelée « backward propagation », consiste quant à elle à propager les vraisemblances v; (t) du domaine dual pour calculer les vraisemblances vir(t + 1) du domaine primal. Plus précisément, au cours d'une première étape, on détermine les valeurs de vraisemblance affectées aux variables connectées à la contrainte de type addition modulo 2 du couple de contrainte considéré, après transformation du domaine réel vers le domaine complexe, dites valeurs de vraisemblance arrière transformées, en appliquant la fonction de non linéarité NL comme détaillé au paragraphe ci-dessus. Les variables internes après transformation s'expriment alors sous la forme : (t) = NL[vi (01 et si (t) = signe[vi (01, pour / e {1, ..., d} , avec d le degré de contraintes associé à la contrainte addition modulo 2 du couple de contraintes considéré. La variable d'entrée ou de sortie du couple de contraintes considéré après transformation B iwi3(t) s'exprime quant à elle sous la forme wo (t)= d et se (t)=(rt/ 1 si3(t)). Au cours d'une deuxième étape, on détermine une valeur de vraisemblance arrière vo (t) affectée à la variable d'entrée ou de sortie du couple de contrainte considéré, mettant en oeuvre une transformation du domaine complexe vers le domaine réel de la variable d'entrée ou de sortie après transformation, en appliquant à nouveau la fonction de non linéarité NL vp (t) = s0 X NL[we (01. Au cours d'une troisième étape, on détermine les valeurs vi(t + 1) en mettant à jour les valeurs des variables internes connectées à la contrainte de type égalité, à partir d'une addition dans le domaine réel, telle que : v-(t) vi(t + 1) = vo (0)+ ve d -1 (t)+ pour 1 E ..., d} , avec d le degré de contraintes associé à la contrainte égalité du couple de contraintes considéré. Finalement, au cours d'une quatrième étape, on met à jour la variable d'entrée ou de sortie du couple de contraintes considéré : vo (t + 1) = vp (t) + v0(0)+ ve(t). Comme déjà indiqué, ces deux phases de propagation avant et arrière sont mises en oeuvre pour chaque couple de contraintes au sein d'une itération de décodage. Plusieurs itérations de décodage peuvent être mises en oeuvre, jusqu'à atteindre un nombre prédéfini d'itérations ou une convergence des variables internes par exemple. Plus précisément, une fois qu'on a effectué ces calculs pour les couples de contraintes connectés à une variable d'entrée du graphe primal-dual (c'est-à-dire qu'on a effectué les calculs concernant les arbres connectés aux valeurs initiales des bits d'information pour la première itération), les nouveaux couples de valeurs primales et duales (vi (t + 1), vi (t + 1)) sont transmis aux couples de contraintes connectés à une variable de sortie du graphe primal-dual (c'est-à-dire aux arbres connectés aux valeurs initiales des bits de redondance pour la première itération) au cours de la phase de propagation avant. Ceci permet d'effectuer un nouvel ensemble de calculs sur ces arbres en notant l'inversion de côté des contraintes de type égalité et addition modulo 2, ce qui revient à permuter les valeurs duales et primales telles que (vi (t + 1), y/ (t + 1)) au cours de la phase de propagation arrière. L'algorithme de décodage proposé converge vers le mot de code le plus proche du mot reçu très rapidement (en quelques itérations), quelque soit la longueur du code. Par exemple, alors que selon l'art antérieur parfois 20 à 50 itérations étaient nécessaires pour obtenir une convergence de l'algorithme de décodage d'un code LDPC, l'invention permet d'obtenir une convergence rapide, en seulement quelques itérations (de l'ordre de 5 ou 6 itérations par exemple). Des saturations à des valeurs suffisamment grandes comme ±4, limitant les valeurs possibles en sortie à l'intervalle [-4,+4], permettent donc de rendre stable et robuste l'algorithme. L'algorithme converge vers ces valeurs de saturation. A nouveau, on note qu'il est possible d'utiliser la technique de décodage décrite comme nouvelle technique de codage, en forçant à zéro les bits de redondance. Ainsi, considérant un graphe primal-dual définit selon l'invention, comprenant des variables d'entrée correspondant aux données source et des variables de sortie correspondant à des données nulles, on obtient après mises en oeuvre d'au moins une itération des phases de propagation avant et arrière décrites ci-dessus des données de redondance au niveau des variables de sortie. 5.5 Structure des codeurs et décodeurs On présente finalement, en relation avec les figures 6A et 6B, les structures simplifiées d'un codeur et d'un décodeur selon un mode de réalisation particulier de l'invention. Le codeur comprend une mémoire 61 comprenant une mémoire tampon, une unité de traitement 62, équipée par exemple d'un microprocesseur l_tP, et pilotée par le programme d'ordinateur 63, mettant en oeuvre le procédé de codage selon un mode de réalisation de l'invention. A l'initialisation, les instructions de code du programme d'ordinateur 63 sont par exemple chargées dans une mémoire RAM avant d'être exécutées par le processeur de l'unité de traitement 62. L'unité de traitement 62 reçoit en entrée des données source S et délivre des données codées comprenant les données source et de données de redondance. Le microprocesseur de l'unité de traitement 62 met en oeuvre les étapes du procédé de codage décrit précédemment, selon les instructions du programme d'ordinateur 63, pour coder les données source. Pour cela, le codeur comprend, outre la mémoire tampon 61, un module de codage mettant simultanément en oeuvre un code correcteur d'erreur et son code dual. Ce module est piloté par le microprocesseur de l'unité de traitement 62. Le décodeur comprend une mémoire 64 comprenant une mémoire tampon, une unité de traitement 65, équipée par exemple d'un microprocesseur l_tP, et pilotée par le programme d'ordinateur 66, mettant en oeuvre le procédé de décodage selon un mode de réalisation de l'invention. A l'initialisation, les instructions de code du programme d'ordinateur 66 sont par exemple chargées dans une mémoire RAM avant d'être exécutées par le processeur de l'unité de traitement 65. L'unité de traitement 65 reçoit en entrée des données codées comprenant des données source et de données de redondance, et délivre une estimation des données source S . Le microprocesseur de l'unité de traitement 65 met en oeuvre les étapes du procédé de décodage décrit précédemment, selon les instructions du programme d'ordinateur 66, pour décoder les données codées. Pour cela, le décodeur comprend, outre la mémoire tampon 64, un module de décodage mettant simultanément en oeuvre un code correcteur d'erreur et son code dual. Ce module est piloté par le microprocesseur de l'unité de traitement 65. REVENDICATIONS1. Procédé de décodage de données codées, comprenant des données source (S) et des données de redondance (R), lesdites données de redondance (R) étant obtenues en appliquant, au codage, un code correcteur d'erreurs auxdites données source (S), mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance (P) inversible, caractérisé en ce que ledit procédé comprend une étape de décodage (12) mettant simultanément en oeuvre ledit code correcteur d'erreurs, dit code correcteur d'erreurs primal, et un code dual dudit code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. 2. Procédé de décodage selon la revendication 1, caractérisé en ce que ladite étape de décodage met en oeuvre un graphe de Tanner modifié, dit graphe primal-dual, obtenu en superposant un graphe de Tanner correspondant audit code correcteur d'erreurs primal, dit graphe primal, et un graphe de Tanner correspondant audit code correcteur d'erreurs dual, dit graphe dual, ladite superposition consistant à connecter deux à deux les variables d'entrée (x0, ..., x5) dudit graphe primal aux variables d'entrée (x0, ..., x5) dudit graphe dual correspondant aux mêmes données source, et à connecter deux à deux les variables de sortie (r0,..., r5) dudit graphe primal aux variables de sortie (r0, ..., r5) dudit graphe dual correspondant aux mêmes données de redondance. 3. Procédé de décodage selon la revendication 2, caractérisé en ce que chaque variable d'entrée ou de sortie du graphe primal-dual est connectée à un couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2, et en ce que ledit procédé comprend une étape de détermination d'au moins une variable interne audit graphe primal-dual mettant en oeuvre : - pour une contrainte de type addition modulo 2 : o une étape de transformation, du domaine réel vers le domaine complexe, des variables connectées à ladite contrainte de type addition modulo 2 ; o une étape d'addition dans le domaine complexe ; o une étape de transformation, du domaine complexe vers le domaine réel, du résultat de ladite étape d'addition ; pour une contrainte de type égalité : o une étape d'addition dans le domaine réel. 4. Procédé de décodage selon la revendication 3, caractérisé en ce que lesdites étapes de transformation mettent en oeuvre une fonction non linéaire NL telle que :IN[1- exp[x] L[x] = log] si x # 0 1+ exp[x] NL[x] = 0 si x = 0 5. Procédé de décodage selon la revendication 2, caractérisé en ce que chaque variable d'entrée ou de sortie du graphe primal-dual est connectée à un couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2, et en ce que ledit procédé met en oeuvre au moins une itération de décodage, comprenant : une phase de propagation avant, consistant à propager les variables internes au graphe primal-dual connectées à une contrainte de type égalité pour mettre à jour les variables internes au graphe primal-dual connectées à une contrainte de type addition modulo 2, et une phase de propagation arrière, consistant à propager les variables internes au graphe primal-dual connectées à une contrainte de type addition modulo 2 pour mettre à jour les variables internes au graphe primal-dual connectées à une contrainte de type égalité. 6. Procédé de décodage selon la revendication 5, caractérisé en ce que ladite phase de propagation avant met en oeuvre les sous-étapes suivantes, pour chaque couple de contraintes connecté à une variable d'entrée ou de sortie du graphe primal-dual : détermination d'une valeur de vraisemblance avant vo (t) affectée à ladite variable d'entrée ou de sortie du graphe primal-dual, mettant en oeuvre une addition dans le domaine réel, telle que : voF (t) - , d détermination de valeurs de vraisemblance affectées aux variables connectées à la contrainte de type addition modulo 2, après transformation du domaine réel vers le domaine complexe, dites valeurs de vraisemblance avant transformées, telle que : o ladite variable d'entrée ou de sortie du graphe primal-dual après transformation s'exprime sous la forme : we'(t) = NL[v 0(0) + vp (01 et so (t) = signe[vo(0)+ vo (01; o les variables internes après transformation s'expriment sous la forme : wr (t) = NL[vii (t)] et sr (t) = signe[vii(t)] pour 1 E {1, ..., d} ; mise à jour des valeurs de vraisemblance avant transformées, mettant en oeuvre une addition dans le domaine complexe, telle que :(t +1) = (t) + wr (t) et sr(t + 1) = ric/ 0,i/sr(t) d -1 mise à jour des valeurs des variables internes connectées à la contrainte de type addition modulo 2, mettant en oeuvre une transformation du domaine complexe vers le domaine réel des valeurs de vraisemblance avant transformées mises à jour, telle que : vi(t +1) = sr (t + 1)NL[wr (t +1)1; avec : v0(0) une valeur initiale affectée à ladite variable d'entrée ou de sortie du graphe primal- dual ; vi(t) les variables internes connectées à la contrainte de type égalité ; vi (t) les variables internes connectées à la contrainte de type addition modulo 2 ; d le degré de contraintes associé à une contrainte du graphe primal-dual ; i / ; NL[x] une fonction non linéaire telle que : NL[x]= log[1- exp[x]] si x # 0 1+ exp[x] NL[x]= 0 si x = 0 7. Procédé de décodage selon la revendication 5, caractérisé en ce que ladite phase de propagation arrière met en oeuvre les sous-étapes suivantes, pour chaque couple de contraintes connecté à une variable d'entrée ou de sortie du graphe primal-dual : détermination de valeurs de vraisemblance affectées aux variables connectées à la contrainte de type addition modulo 2, après transformation du domaine réel vers le domaine complexe, dites valeurs de vraisemblance arrière transformées, telle que : o les variables internes après transformation s'expriment sous la forme : wir3(t) = NL[v;(t)] et 43(0 = signe[vii (01, pour / G {1, , d} ; o ladite variable d'entrée ou de sortie du graphe primal-dual après transformation d_ 1 w1B (t) d s'exprime sous la forme : w0 (t) = 1 d et so (t)=(111=1 siB (0) ; détermination d'une valeur de vraisemblance arrière vo (t) affectée à ladite variable d'entrée ou de sortie du graphe primal-dual, mettant en oeuvre une transformation du domaine complexe vers le domaine réel de la variable d'entrée ou de sortie après transformation, telle que :11:11-3 (t)= soB(t) X NL[w6 (01 mise à jour des valeurs des variables internes connectées à la contrainte de type égalité, mettant en oeuvre une addition dans le domaine réel, telle que : v d vi(t +1) = v0(0)+ ve (t)+.Zwi=1,ilvt"' d -1 mise à jour de ladite variable d'entrée ou de sortie du graphe primal-dual, telle que : v0 (t +1) = ve:(t)+ vo (0)+ ve(t) ; avec : v0(0) une valeur initiale affectée à ladite variable d'entrée ou de sortie du graphe primal- dual ; 10 vi(t) les variables internes connectées à la contrainte de type égalité ; vi(t) les variables internes connectées à la contrainte de type addition modulo 2 ; d le degré de contraintes associé à une contrainte du graphe primal-dual ; i / ; NL[x] une fonction non linéaire telle que : 15 N[1- exp[x] L[x]= log si x # 0 ] 1+ exp[x] NL[x]= 0 si x = 0 8. Procédé de décodage selon la revendication 5, caractérisé en ce qu'il comprend une étape préalable d'initialisation des variables d'entrée, des variables de sortie, et des variables internes du graphe primal-dual à zéro. 9. Procédé de décodage selon la revendication 2, caractérisé en ce qu'au moins deux variables 20 d'entrée ou de sortie du graphe primal-dual sont connectées à un même couple de contraintes comprenant une contrainte de type égalité et une contrainte de type addition modulo 2. 10. Dispositif de décodage de données codées, comprenant des données source (S) et des données de redondance (R), lesdites données de redondance (R) étant obtenues en appliquant, au codage, un code correcteur 25 d'erreurs auxdites données source (S), mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance (P) inversible, caractérisé en ce que ledit dispositif comprend un module de décodage (12) mettant simultanément en oeuvre ledit code correcteur d'erreurs, dit code correcteur d'erreurs primal, et un code dual dudit 30 code correcteur d'erreurs primal, dit code correcteur d'erreurs dual.5 11. Procédé de codage de données source (S), délivrant des données codées comprenant lesdites données source (S) et des données de redondance (R), caractérisé en ce que ledit procédé comprend une étape de codage (11) mettant simultanément en oeuvre : un code correcteur d'erreurs, dit code correcteur d'erreurs primal, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance (P) inversible, et un code dual dudit code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. 12. Dispositif de codage de données source (S), délivrant des données codées comprenant lesdites données source (S) et des données de redondance (R), caractérisé en ce que ledit dispositif comprend un module de codage (11) mettant simultanément en oeuvre : un code correcteur d'erreurs, dit code correcteur d'erreurs primal, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance (P) inversible, et un code dual dudit code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. 13. Dispositif de codage et de décodage de données, caractérisé en ce qu'il comprend un module de codage et de décodage mettant simultanément en oeuvre : - un code correcteur d'erreurs, dit code correcteur d'erreurs primal, mettant en oeuvre une matrice génératrice comprenant, sous sa forme systématique, une matrice identité et une matrice de passage des données source aux données de redondance (P) inversible, et - un code dual dudit code correcteur d'erreurs primal, dit code correcteur d'erreurs dual. 14. Programme ordinateur comportant des instructions pour la mise en oeuvre d'un procédé de décodage selon la revendication 1 ou d'un procédé de codage selon la revendication 11 lorsque ledit programme est exécuté par un processeur.