Dans l'algorithme Welch-Berlekamp de décodage des codes Reed-Solomon, on donne une liste de points représentant un message avec erreurs sur le dans des emplacements inconnus (et est donné à l'algorithme). La sortie est un polynôme passant par tous les points donnés à l'exception de ceux dans lesquels des erreurs se sont produites.
La méthode consiste à résoudre un système d'équations linéaires de la forme
pour tout où a le degré et a le degré au plus . Les variables sont les coefficients de et .
Pour s'assurer que a le degré on ajoute généralement la contrainte que le coefficient de est 1 au système linéaire ci-dessus. Cependant, dans la pratique, on ne sait pas nécessairement . Une façon inefficace (mais toujours polynomiale) de résoudre ce problème consiste à essayer pour toutes les valeurs commençant par descendant jusqu'à ce qu'une solution soit trouvée.
Ma question est: existe-t-il un moyen plus efficace de déterminer ? Sinon, y a-t-il une modification du système linéaire qui permet d'utiliser une borne supérieure sur au lieu de la valeur exacte?
En particulier, je veux utiliser ce décodeur spécifique pour les codes Reed-Solomon, et non un algorithme complètement différent basé sur d'autres techniques.
En réponse à la réponse de DW, voici mon exemple de travail. Tout est fait modulo 7.
plain message is: [2, 3, 2]
polynomial is: 2 + 3 t^1 + 2 t^2
encoded message is: [[0, 2], [1, 0], [2, 2], [3, 1], [4, 4]]
corrupted message is: [[0, 2], [1, 0], [2, 3], [3, 1], [4, 4]]
L'erreur se situe donc au troisième point.
Lorsque l'équation polynomiale en question est
Et brancher donne le système sous forme matricielle:
[2, 0, 0, 6, 0, 0, 0, 0, 0]
[0, 0, 0, 6, 6, 6, 6, 6, 0]
[3, 6, 5, 6, 5, 3, 6, 5, 0]
[1, 3, 2, 6, 4, 5, 1, 3, 0]
[4, 2, 1, 6, 3, 5, 6, 3, 0]
[0, 0, 1, 0, 0, 0, 0, 0, 1]
La dernière ligne est la contrainte que . En appliquant l'élimination gaussienne, nous obtenons
[1, 0, 0, 0, 0, 0, 1, 4, 0]
[0, 1, 0, 0, 0, 0, 3, 3, 1]
[0, 0, 1, 0, 0, 0, 0, 0, 1]
[0, 0, 0, 1, 0, 0, 2, 1, 0]
[0, 0, 0, 0, 1, 0, 2, 2, 5]
[0, 0, 0, 0, 0, 1, 4, 5, 2]
Et en choisissant 1 pour les deux variables libres, nous obtenons un vecteur de solution de
[2, 2, 1, 4, 1, 0, 1, 1]
Ce qui se traduit par
E is 2 + 2 t^1 + 1 t^2
Q is 4 + 1 t^1 + 0 t^2 + 1 t^3 + 1 t^4
Et ne divise pas Q . Notez que Q facteurs comme
Pour j'obtiens une bonne solution:
system is:
[2, 0, 6, 0, 0, 0, 0]
[0, 0, 6, 6, 6, 6, 0]
[3, 6, 6, 5, 3, 6, 0]
[1, 3, 6, 4, 5, 1, 0]
[4, 2, 6, 3, 5, 6, 0]
[0, 1, 0, 0, 0, 0, 1]
reduced system is:
[1, 0, 0, 0, 0, 0, 5]
[0, 1, 0, 0, 0, 0, 1]
[0, 0, 1, 0, 0, 0, 3]
[0, 0, 0, 1, 0, 0, 3]
[0, 0, 0, 0, 1, 0, 6]
[0, 0, 0, 0, 0, 1, 2]
solution is [5, 1, 3, 3, 6, 2]
Q is 3 + 3 t^1 + 6 t^2 + 2 t^3
E is 5 + 1 t^1
P(x) = 2 + 3 t^1 + 2 t^2 # this is correct!
r(x) = 0
Notez que bien que le contre-exemple ci-dessus ait été généré par du code que j'ai écrit à partir de zéro (c'était essentiellement la première chose que j'essayais), on peut vérifier que les solutions sont valides à la main, donc même si mon code est bogué, c'est toujours un contre-exemple valide pour la revendication que l'utilisation de fonctionne.
la source
Réponses:
Pour documenter le résultat de la conversation sur le "contre-exemple" dans la question:
Ainsi, le contre-exemple n'est pas réellement un contre-exemple, et il ne contredit pas ma réponse ci-dessus.
la source