FR 2897742 A1 20070824 FR 0601435 A 20060217 CODAGE/DÉCODAGE PERFECTIONNÉS DE SIGNAUX NUMÉRIQUES, EN PARTICULIER EN QUANTIFICATION VECTORIELLE AVEC CODES À PERMUTATION La présente invention concerne le codage/décodage de signaux numériques tels que des signaux audio, vidéo, et plus généralement les signaux multimédia pour leur stockage ou leur transmission. Elle propose en particulier une solution aux problèmes que posent le codage et le décodage de codes à permutations. to Typiquement, la présente invention s'applique aussi au dual du codage de source : le codage canal ou "modulation". Le codage/décodage en compression de signaux numériques au sens de l'invention peut être très utile pour la quantification des coefficients de transformée des codeurs fréquentiels des signaux de parole et/ou audio. 15 Quantification vectorielle Une solution très répandue en compression des signaux numériques est la quantification vectorielle. La quantification vectorielle représente un vecteur d'entrée par un vecteur de même dimension choisi dans un ensemble fini. Un quantificateur à m niveaux (ou "vecteurs-codes") est une application non 20 bijective de l'ensemble des vecteurs d'entrée, généralement l'espace réel euclidien à n dimensions R " ou un sous-ensemble de R n dans un sous-ensemble fini Y de R " à m éléments distincts: Y= {yo, yi,... Y est appelé alphabet de reproduction ou "dictionnaire" ou encore "répertoire", et ses éléments des vecteurs-codes, des "mots de code" (ou "points de sortie" 25 ou encore "représentants"). Le débit par dimension r du quantificateur (ou sa "résollution") est défini par la relation : r = 1 loge M n En quantification vectorielle, un bloc de n échantillons est traité comme un vecteur de dimension n. Selon la théorie du codage de source, quand la dimension devient très grande, la performance de la quantification vectorielle approche la borne de débit-distorsion de la source. Les dictionnaires de quantificateurs vectoriels peuvent être conçus à partir de méthodes statistiques telles que l'algorithme de Lloyd généralisé (ou "GLA" pour "Generalized Lloyd Algorithm"), basé sur les conditions nécessaires d'optimalité d'un quantificateur vectoriel. Les quantificateurs vectoriels to statistiques ainsi obtenus ne possèdent aucune structure, ce qui rend leur exploration coûteuse en ressources de calculs et de stockage, car la complexité tant du codage que du stockage est proportionnelle à n2' En référence à la figure 1A, trois opérations majeures mettent en oeuvre un quantificateur vectoriel : deux au codage et une au décodage. Un vecteur is d'entrée est codé (étape COD) en choisissant d'abord un vecteur-code dans un dictionnaire. Il s'agit du vecteur-code qui lui "ressemble" le plus (opération OP1 de la figure 1A). Puis, l'indice de ce vecteur-code est déterminé (étape IND) pour être transmis ou stocké (opération OP2 de la figure 1A). Au décodeur (étape DEC), le vecteur-code est déterminé à partir de son indice 20 (opération OP3 de la figure 1A). En modulation, on retrouve les trois opérations OP1, OP2 et OP3 de la figure 1A mais dans un ordre différent. En référence à la figure 1B qui illustre cette dualité modulation/quantification, on prévoit le passage d'un index à un 25 vecteur code (étape COD' de la figure 1B correspondant à l'opération OP3 de la figure 1A). Puis, après transmission sur un canal bruité, on recherche le vecteur-code le plus proche du vecteur reçu (étape DEC' de la figure 1B correspondant à l'opération OP1 de la figure 1A). Finalement, le décodage de l'index du vecteur code est la troisième étape (étape IND' de la figure 1B correspondant à l'opération OP2 de la figure 1A). La croissance exponentielle de la complexité, en fonction de la dimension des vecteurs et du débit, limite l'emploi des quantificateurs vectoriels non structurés à de faibles dimensions et/ou des bas débits pour pouvoir les mettre en oeuvre en temps réel. Pour un quantificateur vectoriel non structuré, la recherche du plus proche voisin (opération OP1) requiert une recherche exhaustive parmi tous les éléments du dictionnaire pour sélectionner l'élément to du dictionnaire qui minimise une mesure de distance entre lui et le vecteur d'entrée. Les deux dernières opérations (indexation OP2 et opération inverse OP3) s'effectuent en général par de simples lectures de tables mais sont néanmoins coûteuses en espace mémoire. Pour s'affranchir des contraintes de taille et de dimension, plusieurs variantes de la quantification vectorielle de 15 base ont été étudiées. Elles tentent de remédier à l'absence de structure du dictionnaire et parviennent ainsi à réduire la complexité, mais au détriment de la qualité. Cependant, le compromis performance/complexité est amélioré, ce qui permet d'accroître la plage des résolutions et/ou des dimensions sur laquelle la quantification vectorielle peut être appliquée efficacement en coût. 20 Beaucoup de schémas de quantificateurs vectoriels structurés ont été proposés et, en particulier, le quantificateur vectoriel mettant en oeuvre un "code à permutation" présenté ci-après. Codes à permutation Dans le quantificateur vectoriel "code à permutation", les vecteurs-codes sont 25 obtenus par permutations des composantes du premier vecteur-code (au sens de l'ordre lexicographique) appelé "leader" (ou encore "vecteur-directeur"). Ces composantes prennent leurs valeurs dans un alphabet A = {ao, al, ..., aq_1} de taille q (alphabet A q-aire tel que a; # ai pour i # j). Les composantes a; sont des réels (ou entiers). Le poids vl/ (où i est un indice allant de 0 à q-1) est le nombre de répétitions de la lettre a; de l'alphabet. Les poids W' sont des entiers q-1 positifs tels que E w' = n . Par convention, les valeurs de l'alphabet satisfont =o ao > al > ... > aq_l . Les n composantes du leader vont en décroissant de la position 0 à la position (n-1). Le vecteur-leader yo est donc un vecteur de 5 dimension n de la forme : w0 w1 wy-1 yo = (1) On comprend qu'un ordre différent des composantes pourrait être choisi, par exemple ao . a, > ... > aq_I ? 0). Les autres vecteurs-codes sont aussi obtenus par permutations des composantes de yo en leur assignant toutes les combinaisons de signes possibles. Le nombre total M de permutations est : (3) fJ w`! ,=0 avec h=n si aq_1>0 et h = n-w4-' sinon. 20 Dans ce cas, le vecteur-leader est aussi appelé leader absolu...DTD: Le quantificateur vectoriel "code à permutation" a été étendu à la composée (ou l'union) de codes à permutation et récemment, cette structure en union de codes à permutation a été généralisée à la quantification vectorielle à dimension et résolution variables (dans le document \NO-04/00219 au nom de la Demanderesse). Les codes à permutation ne sont pas seulement utilisés en quantification vectorielle statistique. On les retrouve aussi en quantification vectorielle algébrique, laquelle utilise des dictionnaires fortement structurés, issus de réseaux réguliers de points ou des codes correcteurs d'erreur. Les codes à permutation sont aussi employés en modulation. to L'exploitation de la structure à code permutation permet le développement d'algorithmes de recherche du plus proche voisin (opération OP1 sur la figure 1A) optimaux et rapides. Cependant l'indexation (ou numérotage) des vecteurs-codes (opération OP2 de la figure 1A) et l'opération inverse de décodage (opération OP3 de la figure 1A) nécessitent plus de calculs que 15 dans le cas des quantificateurs vectoriels non structurés. II existe plusieurs façons d'énumérer les permutations. L'algorithme de Schalkwijk constitue l'une de ces façons: "An algorithm for source coding", par Schalkwijk J.P.M, dans IEEE Trans. on Information Theory, vol. IT-18, N 3, pp 395-399, mai 1972. 20 Utilisant l'analyse combinatoire, ces techniques permettent d'indexer un vecteur-code d'un code à permutation (opération OP2) et aussi d'effectuer l'opération inverse-décodage de l'indice (opération OP3). Parmi les algorithmes d'indexation des permutations, on rappelle ci-dessous l'algorithme de Schalkwijk communément utilisé, par exemple dans les normes : 25 - [3GPP TS 26.273] (ANSI-C code for the Fixed-point Extended AMR - Wideband (AMR-WB+) codec; V6.1.0 (2005-06) (Release 6)), -et [3GPP TS 26.304] (Extended Adaptive Multi-Rate - Wideband (AMR-WB+) codec; Floating-point ANSI-C code; V6.3.0 (2005-06) (Release 6), juin 2005). Calcul du rang d'une permutation au codaqe (opération OP2 de la figure 1A) On cherche à ordonner et indexer toutes les permutations possibles des composantes du vecteur y = (yo, yl ..., yn_1) . L'ordre est lexicographique et l'indice est appelé ici "rang". Calculer le rang du vecteur y revient à calculer le 5 rang du vecteur D = (do,d1...,dn_1) associé à y, et tel que dk vaut une valeur d'indice d si et seulement si yk = ad. Par exemple, un vecteur y, de dimension n=8, a les composantes suivantes : y = (4, 2, 4, 0, 0, 4, 4, 2) L'alphabet à q = 3 lettres (composantes de valeurs différentes) est donné par io A={4, 2, 0}avec ao=4, al =2eta2=0. On associe alors au vecteur y le vecteur D = (0, 1, 0, 2, 2, 0, 0, 1) dont les composantes sont simplement les indices des q lettres de l'alphabet A. Les rangs de y et D sont identiques, mais la définition du vecteur D permet de se ramener au calcul du rang d'une séquence D ayant ses valeurs dans un 15 ensemble {0, 1, ..., q-1} (comportant le même nombre d'éléments que l'alphabet {ao, al, ..., a4_1}). Les poids des vecteurs y et D sont identiques puisque les occurrences de leurs composantes respectives sont identiques. On définit également le poids intermédiaire (wk,wk,...,wk-1) comme le poids du vecteur de composantes 20 (Yk,Yk+1 • • •,Yn-1) et qui correspond donc au vecteur y tronqué en laissant les positions k à n-1. Ainsi : n-1 n-1 wk = 8(y,,ad) = 8(d,,d) (4) où 8(x, y) est l'opérateur de Kronecker (8(x, y) =1 si x=y et 0 sinon). On a : (wo0 , 1 ,...,wo wo ~~-1) (w 0 1 y-1 ). Le rang t du vecteur y peut être calculé, en analyse combinatoire, par la formule : n-1 dk -1 t = E (nù1 (Jk)! , avec 1 10= 0 et (-1)! = k=0 d=0 (wk ù 1)! J(wk!) l i=0 imd Cette formule peut être simplifiée ainsi : n-1 ! /dk-1 n-1 tù~(n- 1ùk). ~wk ùIkk k=0 n (w!) d=0 / k=0 i=0 C'est cette dernière formulation qui est souvent utilisée et que l'on adoptera donc ici. Dans la suite, Ikk sera appelé "rang partiel d'ordre k", avec : kk ù (nù1ùk)! dk-1 d I 9-1 E wk 11(w!) d=0 i=0 to Décodage du ranq (opération OP3) : Détermination de la permutation à partir de son index Le décodage du rang t consiste à retrouver le vecteur D = (do, dl, ., dn_1) associé à y. Une méthode de recherche séquentielle des dk est décrite ci- après. On détermine d'abord la composante do, puis la composante di, ... 15 jusqu'à la composante dn_1. * Détermination de do : Pour trouver do on exploite la formule comportant les inégalités : (5) (6) (7) (n -1)! 'd -i w 0d≤t (Id' Simplifications de l'énumération des permutations et de l'opération inverse (opérations OP2 et OP3) Le calcul du rang t d'une permutation et l'opération inverse ont été accélérés grâce à des simplifications décrites ci-après du calcul des n termes Ikk , dont 5 on rappelle ci-dessous la définition : (n-1-k)! [d'-i ,r-t Ewk , (0 La technique de codage est illustrée sur les figures 2 et 3. En particulier, la figure 3 vise le calcul du rang par la formule de Schalkwijk au sens de l'art antérieur, tandis que la figure 2 illustre des étapes préliminaires (référencées "EP-n", pour n'éme étape préliminaire de la figure 2). 15 En référence à la figure 2, le traitement commence par l'étape EP-1 à laquelle on récupère les composantes du vecteur y = (yo, . , yn_1). L'étape générale EP-2 suivante vise le calcul des éléments de l'alphabet intervenant pour ce vecteur A = (ao, ... , aco), ainsi que la détermination de l'indice maximum q. Pour ce faire, après une étape d'initialisation EP-21 où l'on fixe q = 1, ao = Yo,..DTD: 20 on mène une boucle sur l'indice k compris entre 1 et n-1 (test de fin EP-26 et incrémentation sinon EP-27) visant à rechercher les éléments aq comme suit : - si la composante yk d'indice courant k n'est pas déjà parmi les éléments {ao, ..., aq_i} (test EP-22), il faut attribuer un nouvel élément aq tel que aq = yk (étape EP-23) à l'ensemble {ao, ..., aq_i} et incrémenter q de 1 (étape EP-25). 25 L'étape préliminaire générale suivante EP-3 est un calcul du vecteur D = (do, ..., dä_1) comme suit : -pour k allant de 0 à n-1 (étape EP-31 d'initialisation de k à 0, test de fin EP-37 sur k et incrémentation EP-38 sinon), on trouve la valeur d dans l'ensemble des indices (0, ..., q-1) tel que yk = ad (test EP-33 dans la boucle sur d et incrémentation EP-36 avec d=d+1 sinon). En référence à la figure 3, on décrit maintenant les étapes de calcul du rang t qui suivent le traitement préliminaire illustré sur la figure 2. Les étapes illustrées sur la figure 3 sont référencées "CA-n" pour désigner une n'ème étape de codage au sens de l'art antérieur. lo L'étape CA-1 est une initialisation du rang t à 0, de la variable P à 1 (dénominateur intervenant dans le calcul du rang à l'étape CA-13) et des q poids wo,...,wa_~à0. Ici, on décrémente la valeur de k de n-1 à 0 (étape CA-2 d'initialisation de k à n-1, test de fin CA-14 et décrémentation CA-15 sinon). L'intérêt d'un tel choix 15 est décrit plus loin. On utilise ensuite les composantes dk du vecteur D obtenu à l'étape préliminaire EP- 3 et on fixe, pour un k courant, d= dk (étape CA-3). On met à jour le poids associé Wd (wd=wd+1 à l'étape CA-4) pour estimer le terme P (P=Px wd à l'étape CA-5). On initialise à 0 (étape CA-6) la somme S intervenant au numérateur dans le calcul du rang correspondant à l'étape 20 CA-13 et l'on fait courir ensuite une boucle sur l'indice i des poids w; (test de fin CA-9 et incrémentation CA-10 sinon, jusqu'à d-1) pour mettre à jour la somme S (S=S+w; à l'étape CA-8). Avant de calculer le rang t à l'étape CA-13, on vérifie que la somme S est non nulle (test CA-11). L'intérêt de cette réalisation est décrit plus loin. 25 Le calcul du rang t (étape CA-13) fait intervenir le terme factoriel (n-k-1)! comme suit : t = t + (S/P) (n-k-1)!5 Au lieu de calculer le terme (n-k-1 )! à chaque mise à jour du rang t, on préfère pré-enregistrer en mémoire ces valeurs et avoir recours à un simple accès mémoire (étape CA-12) pour obtenir la valeur courante de (n-k-1)! Ainsi, certains des avantages du traitement illustré sur la figure 3 seront repris 5 dans la mise en oeuvre de la présente invention. On précise ces avantages ci-après. * Stockage des factorielles Pour éviter de calculer les termes (n -1- k)! et wk! en temps réel, les valeurs des n+1 factorielles (0!, 1!, ..., n!) sont pré-calculées et stockées. Si la to dimension n varie de façon bornée (n - nul. Or, ce terme est souvent nul en particulier pour les dernières positions (k 15 proche de n-1). En ajoutant un test sur la nullité de ce terme (test CA-11 de la figure 3), on peut éviter une division (et une multipliication). Le calcul d'une division étant nettement plus complexe qu'un test, le gain en complexité est important. * Inversion de la boucle sur les positions au codage 20 Les poids wk (avec d=0,1,..,q-1) se déduisent des poids wk+, en incrémentant wk+, de 1 et en répétant les autres valeurs wk , pour d # dk. On peut alors réaliser une boucle (étapes CA-2, CA-14, CA-15 de la figure 3) partant de la dernière position (k = n-1) du vecteur jusqu'à la première (k = 0). Cette boucle "inversée" permet de calculer, après une initialisation à 1, le terme Pk tel que : q-1 Pk = fl(wk!), i=0 avec seulement une incrémentation et une multiplication par itération, soit : dk dk wk = wk+l +1 9-1 N-1 et Pi, = wkk x Pk+1 (fl(wk!)= wkk xfJ(wk+1!)) i=0 i=0 s On peut, de plus, traiter les deux dernières positions à part (k = n-1 et k = n-2). En effet, dä_, -1 ^ pour la dernière position (k = n-1), wn 1 = 0 donc = 0 d=0 ^ pour l'avant-dernière position (k = n-2), d,,_2-1 q-1 y-3 si dn_2 > dn_1, wn 2 =1 et Pn_2 = fJ(wä_1!) =1!.1!. [J(0!) =1 et (n-1-k)!=1 !, d=o i=o i=o 10 donc In "2 =1, dä_Z-1 sinon (dn_2 5 dn_l ), Z wn 2 = 0 donc In 22 = 0 avec Pn_2 = 2 si dn_2 = dn_1 et d=0 Pn_2 =1 sinon. D'autres détails de réalisation avantageux, décrits ci-après, peuvent aussi être 15 prévus. * Suppression d'une division au décodage Afin d'éviter des divisions au décodage lors de la recherche de do, l'inégalité (8) peut être reformulée sous la forme : /d0-1 do (n-1)!x lw0 txfl(w'o!) 20 précision sont utilisées au prix d'un accroissement de la complexité de calcul (les opérations en double précision flottante coûtent plus cher que leurs équivalents en précision entière) et de la capacité mémoire requise. De plus, si la précision maximale est limitée à des entiers 32 bits non signés (mise en oeuvre par un processeur en virgule fixe), les factorielles d'entiers supérieurs à 12 ne peuvent pas être pré-stockées directement et les vecteurs de dimension supérieure à 12 doivent être codés à part. Pour résoudre ce problème, une solution plus élaborée utilise une représentation à virgule pseudo flottante par mantisse et exposant des factorielles n! sous la forme 2' x r . Cette décomposition est détaillée dans la table 2 ci-après. Le stockage de n! (pour n inférieur ou égal à 16) se réduit à stocker r avec une précision de 30 bits au maximum ainsi que l'exposant j qui correspond à un simple décalage de bits. n loge (n!) r Iog2 r 0 0 1 1 0 1 0 1 1 0 2 1.0000 2 1 0 3 2.5849625 2 3 1.5849625 4 4.5849625 8 3 1.5849625 5 6.9068906 8 15 3.9068906 6 9.4918531 16 45 5.4918531 7 12.299208 16 315 8.29920802 8 15.299208 128 315 8.29920802 9 18.469133 128 2835 11.469133 21.7910611 256 14175 13.7910611 11 25.2504927 256 155925 17.2504927 12 28.8354552 1024 467775 18.8354552 13 32.535895 1024 6081075 22.535895 14 36.3432499 2048 42567525 25.3432499 40.2501405 2048 638512875 29.2501405 16 44.2501405 32768 638512875 29.2501405 to Table 2: Factorisation des factorielles Ainsi, les techniques de l'art antérieur, pour la plupart, ne résolvent pas le problème du cadrage des variables en précision limitée, en particulier en virgule fixe. La mise en oeuvre dans le codeur TDAC, même si elle permet de régler le problème de cadrage, n'évite pas de coûteuses divisions de deux entiers. De plus, pour des dimensions élevées, les calculs intermédiaires (par exemple le numérateur et le dénominateur des rangs partiels Ikk) peuvent approcher de la saturation. Dans ce cas, les simplifications présentées ci-avant ne peuvent être utilisées dans le traitement de décodage et il est nécessaire de revenir aux formulations des inégalités (8) et (9), et donc de to réintroduire de nombreuses divisions. Les techniques d'énumération autres que celle de Schalkwijk souffrent des mêmes problèmes. Dès lors qu'elles utilisent aussi l'analyse combinatoire, elles nécessitent le calcul de produits de factorielles ainsi que leurs divisions. 15 La présente invention vient améliorer la situation. Elle propose tout d'abord à cet effet un procédé de codage/décodage de signaux numériques, utilisant des codes à permutation faisant intervenir un calcul d'expressions combinatoires, dans lequel ces expressions 20 combinatoires sont représentées par des décompositions en puissances de facteurs premiers, et déterminées par une lecture mémoire de représentations préenregistrées de décompositions de nombres entiers choisis. La présente invention fournit alors une solution efficace aux problèmes liés 25 aussi bien à l'indexation d'un code à permutation qu'à l'opération inverse. Elle résout, conjointement, les deux problèmes de cadrage des variables et de divisions. En effet, dans une réalisation avantageuse, les représentations préenregistrées comportent des valeurs représentatives d'exposants, stockées respectivement en correspondance de valeurs représentatives de nombres premiers successifs, pour chacun desdits entiers choisis. s On résout déjà ainsi le problème lié au cadrage des variables de l'art antérieur. Ce problème de cadrage des variables est d'autant accru lorsqu'il s'agit de manipuler des termes factoriels. Dans une réalisation avantageuse, pour manipuler les expressions combinatoires lorsqu'elles comportent des valeurs factorielles d'entiers, les io représentations préenregistrées comportent au moins des représentations de décompositions de valeurs factorielles. Cette réalisation permet alors de lever la contrainte de cadrage des variables et, de là, repousser les limites habituellement fixées quant à la dimension n des codes à permutation considérés. 15 Selon une autre caractéristique avantageuse, l'une au moins desdites expressions combinatoires comporte un quotient d'un numérateur entier par un dénominateur entier, et ce quotient est représenté par une décomposition de puissances de facteurs premiers, dont chaque puissance est une différence 20 d'exposants respectivement associés au numérateur et au dénominateur et affectés à un même nombre premier. On résout ainsi le problème lié au calcul de divisions de l'art antérieur, en remplaçant ce calcul par un simple calcul de soustractions. 25 Dans un premier mode de réalisation, on prévoit un adressage mémoire pour retrouver une décomposition préenregistrée de l'un des entiers choisis précités. A cet effet, la représentation préenregistrée d'un nombre entier choisi est stockée dans une mémoire adressable, un adressage de ladite mémoire donnant une succession d'exposants à affecter à des nombres premiers respectifs pour recomposer l'entier choisi. Préférentiellement, la représentation préenregistrée d'un nombre entier choisi sera stockée sous la forme d'une succession d'adresses donnant chacune, pour un nombre premier, un exposant à affecter à ce nombre premier pour recomposer l'entier choisi. On dénommera ci-après cette réalisation selon le premier mode par les termes "représentation éclatée des décompositions". En variante, dans un second mode de réalisation, les représentations préenregistrées sont stockées sous la forme de mots comportant unesuccession de groupes de bits, chaque groupe ayant : - un poids fonction d'un nombre premier, et - une valeur fonction d'un exposant à associer à ce nombre premier. Préférentiellement, la détermination des puissances de facteurs premiers s'effectue alors par applications successives d'au moins un masque partiel au mot de bits, avec décalages successifs en fonction des poids des bits et lectures des bits restants. On dénommera ci-après cette réalisation selon le second mode par les termes 20 "représentation compacte des décompositions". Le déroulement même du procédé, pour le calcul d'une expression combinatoire, peut de façon générale être mené selon les étapes : - identification, parmi lesdits nombres entiers choisis, ,de termes apparaissant 25 dans un produit et/ou un quotient formant ladite expression combinatoire, - lecture mémoire des exposants intervenant dans les décompositions en facteurs premiers desdits termes, - addition et/ou respectivement soustraction des exposants lus pour déterminer les exposants intervenant dans la décomposition en puissances de facteurs premiers de ladite expression combinatoire, et, de là, calculer ladite expression combinatoire à partir de sa décomposition en puissances de facteurs premiers. Pour ce qui concerne le calcul d'un produit à effectuer de façon récurrente et faisant intervenir à chaque récurrence un nouveau terme, il peut être avantageux de stocker temporairement la décomposition d'un calcul de produit effectué pour une récurrence antérieure. Ainsi, si le procédé comporte une to étape récurrente de calcul d'un produit faisant intervenir à chaque récurrence un terme se multipliant à un produit déterminé à une récurrence antérieure : - ledit produit déterminé à une récurrence antérieure est avantageusement conservé en mémoire sous la forme d'une décomposition en puissances de facteurs premiers, 15 - ledit terme se multipliant au produit est préférentiellement l'un des entiers choisis dont la décomposition est préenregistrée, et - pour déterminer ledit produit à une récurrence courante, il suffit d'additionner, un à un par nombre premier, les exposants issus des décompositions respectives dudit produit déterminé à une récurrence antérieure et dudit terme 20 se multipliant au produit. De même, si le procédé comporte une étape récurrente de calcul d'une division faisant intervenir à chaque récurrence un terme divisant un quotient déterminé à une récurrence antérieure : -ledit quotient déterminé à une récurrence antérieure est avantageusement 25 conservé en mémoire sous la forme d'une décomposition en puissances de facteurs premiers, - ledit terme divisant le quotient est préférentiellement l'un des entiers choisis dont la décomposition est préenregistrée, et - pour déterminer ladite division à une récurrence courante, on soustrait, un à un par nombre premier, les exposants issus de la décomposition dudit terme aux exposants issus de la décomposition dudit quotient déterminé à une récurrence antérieure. Ce stockage temporaire des décompositions intermédiaires des produits et/ou des quotients calculés de façon récurrente s'avérera particulièrement avantageux pour la détermination de rangs partiels récurrents dont le cumul est représentatif d'un rang d'une permutation. Ainsi, dans une réalisation avantageuse de l'invention, les codes à permutation font intervenir le calcul d'une quantité représentative d'un rang d'une permutation comportant un cumul de rangs partiels, chaque rang partiel correspondant alors à l'une desdites expressions combinatoires. Le calcul du rang d'une permutation peut alors intervenir, au codage de signaux numériques en quantification vectorielle, pour indexer les permutations des composantes d'un vecteur directeur (opération OP2 de la figure 1A), ces permutations ayant été effectuées à une étape préalable (opération OP1) pour déterminer un vecteur-code le plus proche d'un vecteur d'entrée. De même, au décodage de signaux numériques en quantification vectorielle, l'estimation d'un rang d'une permutation intervient dès lors qu'à partir d'une valeur donnée d'un rang de permutation : - on calcule au moins une quantité représentative d'un rang de permutation (opération OP3 de la figure 1A) approchant ladite valeur donnée, en fonction d'au moins une composante présumée d'un vecteur-code à construire, - on valide le choix de la composante présumée si cette quantité vérifie une condition de proximité avec la valeur donnée du rang. Dans un exemple de réalisation, cette condition de proximité est vérifiée si la valeur donnée du rang peut être encadrée par des cumuls de rangs partiels jusqu'au rang partiel associé à la composante présumée, d'une part, et jusqu'au rang partiel associé à une composante correspondant à une incrémentation de la composante présumée, d'autre part. Cette condition de proximité peut donc correspondre à une formulation générale des inégalités (8) décrites ci-avant dans le cas d'une énumération de Schalkwijk. Ainsi, la présente invention peut s'appliquer avantageusement au io codage/décodage de source avec quantification vectorielle au sens de la figure 1A. Toutefois, le codage/décodage peut aussi être de type codage/décodage de canal, en modulation, au sens de la figure 1 B, dès lors qu'il comporte : - avant transmission, la détermination d'un vecteur-code à partir du rang d'une 15 permutation (même opération OP3 de la figure 1A et de la figure 1 B), et - en réception, le calcul du rang d'une permutation à partir d'un vecteur-code correspondant à un vecteur reçu (même opération OP:2 de la figure 1A et de la figure 1B). 20 Le calcul d'un rang partiel fait intervenir des termes (dans un produit ou un quotient) qui, en règle générale comme on le verra plus loin, restent inférieurs ou égaux à la dimension maximale n des codes à permutation. Ainsi, dans une réalisation avantageuse, les nombres entiers choisis dont les décompositions sont préenregistrées comportent au moins : 25 - les entiers compris entre 1 et la dimension maximale n, - la valeur factorielle de l'entier 0, - et, de préférence, les valeurs factorielles des entiers compris entre 1 et la dimension maximale n. Dans une réalisation particulière facultative, les nombres entiers choisis peuvent comporter en outre la valeur de 0. Ainsi, si le code de permutation utilise une énumération de Schalkwijk, un rang s partiel Ikk associé à la troncature (yk,...,yn_1) d'un vecteur-code (yo,...,yn_,) s'écrit : Ikk û (nû1ûk)! dk-ld w où : qù1 k n (wi f) d=0 k i=o m - la notation fl représente un produit pour un indice entier i allant de 0 à m, i=0 m - la notation E représente une somme pour un indice i allant de 0 à m, i=o to - la notation 1! est la valeur factorielle de l'entier 1 avec 1! =lx 2x 3x... x (l-1) x 1, pour 1>0 et 0!=1, - l'entier n est la dimension du code de permutation, correspondant au nombre total de composantes que comporte un vecteur-code, - l'entier k, compris entre 0 et n-1, est l'indice de la kième composante Yk du 15 vecteur-code, recherchée à partir d'une valeur de rang en décodage de source (respectivement en codage de canal) ou dont on cherche à indexer les permutations en codage de source (respectivement en décodage de canal), - l'entier q est le nombre de composantes distinctes que comporte le vecteur-code, et 20 - un terme wkd (dénommé "poids intermédiaire") représente le nombre de composantes d'indices compris entre k et n-1 et qui ont une valeur égale à celle d'une même composante d'indice d. Dans ce cas, les nombres entiers choisis dont les décompositions sont préenregistrées et que l'on cherchera alors à identifier dans l'expression du rang partiel Ikk , dans un produit et/ou un quotient, sont préférentiellement : - les termes factoriels (n - l - k)! , pour tous les entiers k et compris entre 0 et n-1, (c'est-à-dire les valeurs des factorielles pour tous les entiers entre 0 et (n-1)), - la valeur de chaque terme wk et/ou sa valeur factorielle, intervenant dans le q-1 produit Pk = fl(wk!), chaque terme wk étant compris entre 0 et n, et =o dk -1 - les termes Sk = l wk , chacun compris entre 1 et n-1, pour tous les entiers k d=0 io compris entre 0 et n-1. Toujours dans le cas particulier de l'énumération de Schalkwijk, le stockage temporaire des décompositions intermédiaires s'applique avantageusement comme suit : la somme des exposants dans la décomposition du terme q-1 !) est stockée temporairement en mémoire pour un indice k 15 i=0 précédent, pour lui ajouter ou lui soustraire les exposants de la décomposition d'un terme wk pour un indice k courant. D'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après, et des dessins annexés sur lesquels, outre les figures 1A, 1 B, 2 et 3 décrites ci-avant : la figure 4 illustre les principes du codage/décodage du rang d'une 25 permutation en exploitant l'invention, 20 - la figure 5 illustre un traitement pour le codage d'un rang d'une permutation en exploitant l'invention selon un premier mode de réalisation, dans lequel on prévoit une représentation éclatée des exposants de la décomposition en puissances de nombres premiers de termes intervenant dans ce calcul, - la figure 6 illustre un traitement pour le codage d'un rang d'une permutation en exploitant l'invention selon un second mode de réalisation, dans lequel on prévoit une représentation compacte des exposants de la décomposition, la figure 7 illustre un traitement pour le décodage d'un rang d'une permutation en exploitant l'invention selon un premier mode de réalisation, dans lequel on io prévoit une représentation éclatée des exposants de la décomposition, et - la figure 8 illustre un traitement pour le décodage d'un rang d'une permutation en exploitant l'invention selon un second mode de réalisation, dans lequel on prévoit une représentation compacte des exposants de la décomposition, - la figure 9 représente schématiquement un dispositif de codage/décodage 15 mettant en oeuvre la présente invention. Il convient d'insister pour rappeler (en particulier en référence aux figures 4 à 8 ci-après) que : - le terme "codage" vise le calcul du rang t d'une permutation (opération OP2 20 des figures 1A et 1 B), et - le terme "décodage" vise la détermination de la permutation à partir de ce rang t (opération OP3 des figures 1A et 1B). Ces opérations seront ainsi nommées en référence au codage/décodage de source avec quantification vectorielle. On rappelle que ces opérations peuvent 25 être menées aussi en codage/décodage de canal, en modulation. Pour illustrer d'emblée le principe de l'invention, on décrit ci-après la factorisation en puissances de nombres premiers. La décomposition d'un nombre entier positif K non nul, en puissances de 30 nombres premiers, s'écrit ainsi : Pi étant un iéme mK K =n(pi) r=0 nombre premier (po =1, p, = 2, P2 =3, P3 =5, P4 =7, p5 =11, p6 =13, p7 =17 , etc.). On note eK l'exposant de p, dans la décomposition du nombre entier K et mK l'indice du plus grand facteur premier intervenant dans la décomposition de K d'exposant non nul. Par exemple, le nombre K=120 (soit 5!) s'écrit : 120 = 1.23.31.51 et mk=3 ici puisque le plus grand facteur "5" est d'indice 3 (p3 = 5 ). On a donc : e5,=1, e5,=3, e5!=1 ete5,=1. to En pratique, le nombre "1" étant l'élément neutre de la multiplication, on peut supprimer po de la décomposition, soit : K=JJ(P,) 4 r=1 Bien entendu, K=0 n'est pas décomposable en puissances de facteurs premiers. 15 La décomposition en produits de puissances de nombres premiers des entiers positifs inférieurs ou égaux à 16 est donnée dans la table 3a, la décomposition de leurs factorielles dans la table 3b. Cette décomposition fait intervenir 6 nombres premiers (2, 3, 5, 7, 11 et 13). Les colonnes étant indexées par le 20 nombre premier p, et les lignes par n, l'élément de la table 3a (respectivement 3b) à l'intersection de la colonne p, et de la ligne n est l'exposant en (respectivement e,,,) du nombre premier p, dans la décomposition en produit de puissances de nombres premiers du nombre n (respectivement n!). Pour tout nombre entier positif n>1, le nombre mn, de facteurs premiers de n! est tel que : pnn n On rappelle que le rang partiel, noté Ikk , est donné par la relation (7) ci-avant : Idk - (n-1-k)! dk-1Wd k q-1 k fJ(Wki !) d=o i=0 de sorte que trois termes peuvent être décomposés en puissances de nombres premiers. Il s'agit des termes : q-1 dk -1 (n ù 1 ù k)! , Pk = f Wk ! et Sk = E Wk . i=o d=0 A partir des exposants des décompositions de (n -1- k)!, de Pk et de Sk , on calcule par de simples additions et soustractions, les exposants de la décomposition de Ikk En effet, l'exposant ef,,k d'un facteur premier p, dans la décomposition de Ikk k se calcule à partir des trois exposants de pi dans les décompositions des trois termes (n-1 -k)!, Sk et Pk . L'exposant efdk est égal à la somme des exposants de p, des deux premiers termes (numérateur de Ikk) à laquelle on soustrait l'exposant de p, du dernier terme (dénominateur de Ikk ). Cette observation, une fois formalisée, s'écrit : e fkk = e(n-1-k)! + esk ù ePk (10) On a représenté sur la figure 4 les étapes générales (référencées "G-n" pour nième étape générale) pouvant intervenir dans un traitement au sens de l'invention aussi bien pour le codage que pour le décodage. Ainsi, en référence à la figure 4, à partir d'un indice courant k (étape G-1) et 20 via quelques étapes intermédiaires qui seront décrites en détail plus loin, on retiendra surtout que l'on prévoit de se référer à des tables pré-enregistrées (notées D, et D,! à l'étape G-2 illustrée à titre d'exemple sur la figure 4) pour15 calculer l'exposant e/Jk global selon la relation (10) ci-avant (étape G-3), cet exposant étant propre à la décomposition du rang intermédiaire Ikk en puissances de facteurs premiers. On en déduit alors la valeur du rang intermédiaire Ikk (étape G-4), éventuellement en rnenant une boucle sur l'indice i propre aux facteurs premiers. Le calcul de ce rang intermédiaire peut alors se poursuivre par une mise à jour du rang total t de la permutation (étape G-5) : - avec une relation du type t=t+Ikk pour le codage du rang avec décrémentation de l'indice k (étape G-6), - ou avec une relation du type t=t-Ikk pour le décodage du rang avec incrémentation de l'indice k (étape G-6), comme on le verra plus loin. Finalement, on obtient le rang t de la permutation, en codage à l'étape G-7, ou, en décodage (traits pointillés de la figure 4), on déduit des inégalités de la formule (8) ci-avant les composantes dk du vecteur D à l'étape G-8 et, de là, celles du vecteur Y par la relation yk=ad ci-avant. Dans le cas général et indépendamment de l'énumération de Schalkwijk, si un rang partiel t' (t' >0) d'une permutation est sous la forme d'un numérateur de Nt. 20 termes vi (1 _ Ikk = n (p.) 1k i=1 On indique ici que les termes (n -1- k)! et Pk sont des entiers strictement positifs mais le terme Sk peut être nul et donc non décomposable. Dans ce cas, le rang partiel Ikk est nul. On prévoira alors un test sur la valeur du terme 20 Sk (Sk = 0?) pour ne calculer le rang partiel Ikk que si Sk ≠ 0, comme décrit ci-avant (étape CA-11 de la figure 3). 10 N, D,. mr m' L - év Plus généralement, si e,, = Ee, ! -EeP ' alors t'= fJ p; _ fl p, ' ' j=1 j=1 r=l i=1 Q, m, , m, ' der et si Eet alors t'=np;" _flp' '' (t'>0) j=1 i=1 i=1 On retiendra alors que la factorisation en facteurs premiers des termes composant un rang partiel permet de supprimer les divisions en les remplaçant par des multiplications de puissances de facteurs premiers, avec, en particulier, de simples additions et soustractions des exposants associés à ces nombres premiers. l0 Ainsi, au sens de la présente invention, on prévoit les étapes suivantes, à partir d'un nombre restreint de décompositions en facteurs premiers de nombres entiers, stockées en mémoire (dites ci-après "décompositions de base") : 15 - détermination des décompositions en facteurs premiers de termes (tels que q-1 dk -1 (nù1ùk)! , pk =flwk! et Sk = E w, apparaissant dans le rang d'une i=o d=o permutation (dites ci-après "décompositions intermédiaires"), au sens de l'étape G-2 de la figure 4, détermination, à partir de ces décompositions intermédiaires, de la 20 décomposition en facteurs premiers d'un rang partiel t'(Ik k) d'une permutation, notamment par le calcul des exposants intervenant dans la décomposition de ce rang partiel (avec par exemple une relation du type e/dk = + e' - e k ), au sens de l'étape G-3 de la figure 4, et - calcul du rang partiel t'(Idk) à partir de sa décomposition (par exemple en m ak ik e dk utilisant une relation du type Ikk = fl (p;) 'k ), au sens de l'étape G-4 de la i,o figure 4. Bien entendu, les décompositions de base à stocker font préférentiellement l'objet d'une sélection avantageuse. Dans une réalisation préférée mais non limitative, les décompositions de base à stocker seront choisies en fonction de la dimension maximale des codes à permutation considérés (cette dimension maximale étant notée n). Ainsi, les décompositions de base sont io préférentiellement : - les décompositions des factorielles d'un entier 1 (notées 1!), l'entier 1 étant tel que Osln, - et les décomposition de l'entier 1 lui-même, avec cette fois 1 l n , où on rappelle que n est la dimension maximale des codes à permutation 15 considérés. On peut alors identifier une décomposition de base par un nombre m donnant : - le nombre m de facteurs premiers à considérer, - ces m facteurs premiers, eux-mêmes, 20 - et leurs exposants respectifs. Des exemples de cette mise en oeuvre seront décrits plus loin en référence aux tables 4a à 4d, dans le cadre d'une représentation dite "éclatée" des décompositions. On indique qu'une représentation dite "compacte", décrite en détail plus loin, consiste à stocker un simple mot dont les bits donnent tous les 25 exposants intervenant dans une décomposition. On peut définir alors différents jeux de décompositions de base, ainsi que des procédures de représentation et de stockage de ces décompositions de base. Par ailleurs, la sélection des termes dont on détermine les décompositions s intermédiaires, ainsi que la détermination proprement dite de ces décompositions intermédiaires, font l'objet de réalisations avantageuses qui seront décrites plus loin. La décomposition d'un rang partiel et le calcul d'un rang partiel à partir de sa décomposition font aussi l'objet de réalisations avantageuses décrites plus loin. i0 On décrit maintenant le choix des décompositions de base à stocker. De manière générale et indépendamment de la technique d'énumération d'un code à permutation de dimension n, le calcul du rang d'une permutation utilise les nombres entiers l (0 s l n) et surtout leurs factorielles 1! (0 5 l 0) est nécessaire au calcul d'un rang partiel, elle est calculée à l'étape de détermination des décompositions intermédiaires à partir des 1 décompositions de base dej(1j_1)avec e,`, _lei. l=1 25 Inversement, on peut ne prévoir que les (n+1) décompositions de 1! (0 1 s n ). Si une décomposition de 1(1 >0) est nécessaire au calcul d'un rang partiel, elle est calculée à l'étape de détermination des décompositions intermédiaires à partir des deux décompositions de base de 1! et (1-1)! et à partir de la relation : e` = -e(/_1) On comprend donc que le choix du jeu de décompositions de base peut avantageusement résulter d'un compromis entre la minimisation de la mémoire nécessaire au stockage des représentations de ces décompositions de base et à la minimisation de la complexité de l'étape de détermination des décompositions intermédiaires. io On décrit ci-après une représentation des décompositions au sens de l'invention. Comme indiqué plus haut, une décomposition (qu'elle soit de rang partiel, intermédiaire ou de base) est définie par un nombre n7 donnant le nombre de facteurs premiers à considérer, ces m facteurs premiers et leurs exposants 15 respectifs. On propose ci-après différentes solutions pour représenter les décompositions et stocker les données pour les décompositions de base. Représentation éclatée des exposants * Représentation des factorielles 1! (0 1 n ) Le nombre m,, de facteurs premiers intervenant dans la décomposition de la 20 valeur 1! croît avec le nombre 1. Une première solution pour représenter la décomposition de 1! (0≤1≤n) consiste à stocker pour chaque valeur de 1 (0 ≤1≤ n) le nombre m,, et les m,, exposants des puissances de pi (1 ≤i ). On notera que les m,, exposants de 1! sont non nuls. Dans une variante plus avantageuse, le jeu des décompositions de base 25 partage le même nombre m,,, de facteurs premiers et on stocke m,,, exposants pour chaque décomposition de base, les exposants de la décomposition de base de 1! d'indice supérieur à m,, étant nuls. Cette solution permet d'avoir recours à un tableau des exposants en prévoyant un adressage régulier de ce tableau. Toutefois, une telle réalisation nécessite une taille mémoire conséquente. Ce tableau comporte mn, x (n + 1) valeurs et l'exposant e;, est stocké à l'adresse (mn, ,1+(i -1)) de ce tableau, où la notation (x, y) vise la case de ce tableau à la ligne x et à la colonne y. Bien entendu, on comprendra que d'autres conventions peuvent être considérées. Ainsi, au lieu de considérer un tableau à deux dimensions avec m colonnes et N lignes qui comporte donc m x N cases (ou éléments), on peut considérer un tableau à to une dimension ayant m x N éléments, l'élément à l'adresse (x, y) du tableau à deux dimensions se retrouvant alors à l'adresse mx x+ y du tableau à une dimension. L'exposant e;, stocké à l'adresse (1,(i-l du tableau à deux dimensions est alors stocké à l'adresse (mn, x 1+(1--l du tableau à une dimension. Par exemple, les exposants des décompositions des factorielles 15 des nombres 0 à 8 peuvent être stockés dans le tableau à deux dimensions comprenant 36 éléments, composé de 4 colonnes (colonnes pi =2, 3, 5, 7) de la table 3b et de 9 lignes (lignes n=0, ..8). Ces mêmes exposants peuvent être stockés dans un tableau monodimensionnel Dä à 36 éléments, donné ci-après (annexe A-11). L'élément à l'adresse (x, y) du premier tableau étant égal 20 l'élément à l'adresse Dä : 4x x+ y. II peut être prévu, en complément, de stocker les (n+1) valeurs de m,, pour pouvoir réduire les calculs des décompositions intermédiaires utilisant la décomposition de base de 1! . * Représentation des entiers l (1 Représentation compacte des exposants On décrit ci-après une autre solution avantageuse minimisant le stockage et 15 qui consiste à représenter de manière compacte les exposants d'une décomposition de base sur un nombre limité de mots. Dans cette variante de représentation des décompositions de base, les décompositions intermédiaires et celles des rangs partiels sont aussi représentées de façon compacte. Avantageusement, cette solution minimise aussi la complexité de la détermination de ces décompositions, comme on le verra. * Représentation compacte des décompositions On cherche à déterminer, pour chaque facteur premier p;, une borne supérieure fi de la valeur maximale de son exposant dans le numérateur des rangs partiels. Cette borne donne le nombre maximum de valeurs possibles de l'exposant de A, soit Q, + 1. En notant b' , le nombre entier de bits pour to représenter de façon binaire la valeur (/ii +1), on a : i-1 b' = rlog, 66i +1)1) et Bn = > bn (avec B,1, = 0) J =l où rxl dénote l'entier immédiatement supérieur ou égal à x (rxl-1 Les exposants de la décomposition en facteurs premiers d'un terme K intervenant dans un rang partiel t' peuvent être représentés de façon compacte mä, 15 par un mot eK de Bä bits (B,, = Ibn , ce mot eK étant tel que : i=1 eK =E(eK Bn) i=1 La notation " B" représente un décalage à gauche de B bits. II est à noter que si le nombre n est grand, il se peut que B,, soit supérieur au 20 nombre de bits Bo utilisés pour représenter des entiers (16, 32 ou 40 bits). Dans ce cas, les exposants de la décomposition en facteurs premiers d'un entier K intervenant dans t' sont représentés sous la forme de M mots entiers eK (m) , 0 _ 1). Avantageusement, les M mots peuvent être formés de la façon suivante : • eK(0) comprend les io premiers exposants (exposants de pl à pi ): i0 eK (0) = E eK Bn où Bn +l • Cette dernière relation peut être généralisée pour tout m, jusqu'à composer le dernier mot qui comprend l'exposant de pm& Bien entendu, d'autres variantes peuvent être prévues. Par exemple, une variante consiste à stocker à part l'exposant de pl , et appliquer le traitement ci-dessus à partir de l'exposant de p2 . * Détermination de la borne supérieure Les bornes fi, peuvent être déterminées de plusieurs façons. En exploitant les informations sur le code à permutation (taille q de l'alphabet, poids In/ 0 i En s'appuyant sur la lecture du tableau ci-après, la représentation binaire de 18 (=16+2) sur 815 =15 bits est : 000 0000 0001 0010, soit 00 00 00 00 001 0010 en regroupant les bits associés à un même nombre premier. 15...7 6 5 4 3 2 1 0 2' ... 64 32 16 8 4 2 1 En représentation binaire, 0 ... 0 0 0 1 0 0 1 0 18 s'écrit : Les 4 bits de poids faibles (poids i=0 à 3) sont l'exposant du facteur premier 2, soit : 0010 = 2, ce qui signifie que 2 est l'exposant à affecter au nombre premier 2. 10 Les 3 bits suivants (poids i=4 à 6) sont l'exposant du facteur premier 3, soit : 001 = 1, ce qui signifie que 1 est l'exposant à affecter au nombre premier 3. Les 2 bits suivants (poids i=7 à 8) sont l'exposant du facteur premier 5, soit : 00=0 Les 2 bits suivants (poids i=9 à 10) sont l'exposant du facteur premier 7, soit : 15 00=0 Les 2 bits suivants (poids i=7 à 12) sont l'exposant du facteur premier 11, soit : 00=0 Les 2 bits suivants (poids i=13 à 14) sont l'exposant du facteur premier 13, soit : 00 = 0 20 La procédure d'extraction consiste à masquer les bits de poids forts pour récupérer l'exposant du facteur premier contenu dans les bits de poids faible, puis à décaler l'exposant compact du nombre de bits récupérés pour passer à l'exposant du facteur premier suivant. 25 Ainsi en dimension 15, il y a 6 exposants à extraire en commençant par l'exposant de 2.5 La représentation binaire de l'exposant de 2 correspond aux 4 bits de poids faible de 18, soit 0 0 1 0 qui correspond à 2. Pour les récupérer, on masque les bits de poids fort de 18 avec 15 (noté 18 & 15), ce qui équivaut à : 24ù1=1111. On obtient e12 = 18 & (2 4 - 1) = 2, ce qui signifie que 2 est l'exposant à affecter au nombre premier 2. Puis on décale 18 de 4 bits vers la droite, on obtient : 000 0000 0001 = 1 La représentation binaire de l'exposant de 3 correspond aux 3 bits de poids faible de 1, soit 0 0 1 (=1). Pour les récupérer, on masque les bits de poids fort de 1 avec 7 (noté 1 & 7 et valant 23 ù 1 = 1 1 1). On obtient e112 = 1 & (2 3 - 1) = 1, ce qui signifie que 1 est l'exposant à affecter au nombre premier 3. Puis on décale 1 de 2 bits vers la droite, on obtient ensuite : 0000 0000 = 0 pour tous les autres bits de poids fort. On retient donc que les puissances de 1=12 sont : - 2 pour le nombre premier 2, et - 1 pour le nombre premier 3, soit 1=12 = 22x31 * Borne supérieure pour le dénominateur On suppose ici que, pour chaque facteur premier, son exposant dans le dénominateur d'un rang partiel t' est inférieur ou égal à son exposant dans le numérateur de t'. Tel est le cas si t' est strictement positif car on a bien = (e'Num(,') ù enen(,')) 0 , soit donc e(Num(,') > enen(t') En pratique, avec la formule de Schalkwijk et si q>1, la valeur en_,), est une borne supérieure de la valeur maximale de l'exposant du dénominateur Pk (Pk 1). Il suffit donc de vérifier l'inégalité Dans le cas où q=1, on comprendra qu'un seul mot de code de rang connu (1=0) intervient dans le code à permutation et il est donc inutile a priori d'effectuer les calculs de rang et les opérations inverses correspondantes. Cependant, s'il n'est pas souhaité de traiter à part ce cas particulier, il peut toujours être prévu de calculer la valeur b` par le maximum de f3; et e,,,. La table 7 ci-après illustre ce cas pour n=16. Pi 2 3 5 7 11 13 max(e;) 3 2 1 1 1 1 Is! On décrit brièvement maintenant les capacités mémoires requises pour le stockage des décompositions de base. Indépendamment de la solution choisie pour la représentation des décompositions de base, le stockage des décompositions de base s'effectue dans des tableaux et l'on a alors recours à un adressage de ces tableaux lors des opérations de codage et de décodage du rang. Bien que la décomposition de 0 ne soit pas possible (et, d'ailleurs, pas utilisée), il peut être préférable de stocker des exposants "factices" pour la décomposition de 0 (par exemple des 0 ou des 1), et ce, pour simplifier le calcul d'adresse. La table 8 ci-après résume la taille mémoire nécessaire pour le stockage des données relatives aux décompositions de base pour les cinq solutions présentées pour ces deux cas (stockage ou non de la décomposition factice de 0). Solution Pas de stockage de 0 Avec stockage factice de 0 1 ~m1 +Em,, +2n+1 E ms + m~, +2(n+1) r=1 r=o r=o i=o 2 2Em1, +1+(n+l) 2Em~, +(n+1) 3 Em1, +2E m'1 +2n+ 1 Ems, +2E m'1 +2(n+1) r=o r=1 t=o r=o 4 x(2n+1) 2m,,, x(n+1) 5 2n+1+m,,, 2x(n+1)+mn, Table 8 : Mémoire nécessaire au stockage des décompositions de base to Dans la cinquième solution, on a tenu compte du stockage (+mn,) des nombres de bits b` . Cependant, en pratique, plutôt que de les lire dans une mémoire, ceux-ci seront "câblés en dur" (leur valeur étant fixe dans un programme de calcul sans être déclarée en tant que variable), comme on le 15 verra dans les modes de réalisation plus loin. Il est donc apparu inutile de les stocker en pratique. La table 9 donne la mémoire nécessaire pour le stockage des données relatives aux exposants de la décomposition de ces cinq solutions pour nmax = 8 et 15 (avec stockage factice de 0). " n 2m , x(n+1) 2x(n+1)+m n, (m,, +m,) E2m,, Em+E2m'1 . n f=0 !=0 1=0 1=0 +2(n+1) +(n+1) +2(n+1) 8 54 51 57 72 22 15 129 130 129 192 38 Table 9 : Mémoire requise pour le stockage des décompositions de base pour les codes à permutation de dimensions 8 et 15 On décrit maintenant le stockage des puissances de facteurs premiers. En dehors des décompositions de base, la présente invention utilise les to puissances des facteurs premiers pour calculer le rang partiel à partir de sa décomposition. On peut, à partir d'une table de ces facteurs premiers, calculer leurs puissances en temps réel ("en ligne"). Préférentiellement, les puissances des nombres premiers autres que 2 sont pré-calculées et stockées et seules les puissances de 2 sont calculées en temps réel. La table 10a ci-après donne 15 les puissances de 3, 5 et 7 nécessaires pour un code à permutation de dimension 8 (comme ceux utilisés dans le codeur AMR-WB+). La table 10b donne les puissances de 3, 5, 7, 11 et 13 nécessaires pour les codes à permutation de dimension maximale 15 (comme ceux utilisés dans le codeur TDAC). 20 0 1 2 3 3 3 9 27 5 25 7 7 49 Table 10a : Puissances des facteurs premiers pour les codes à permutation de dimension 8 p; 0 1 2 4 5 6 7 Npfp e; 3 3 9 27 81 243 729 2187 8 5 5 25 125 - - _ 4 7 7 , 49 343 - - - 4 11 11 121 - - - - 3 13 13 169 - - - _ 3 Table 10b : Puissances des facteurs premiers pour les codes à 5 permutation de dimension 15. La variable Npfp désigne ici le nombre de puissances de facteurs premiers à stocker. Là, encore, on pourra ne stocker que le nombre nécessaire de puissances pour chaque facteur premier. En variante, s'il est préféré de ne disposer que to d'un seul tableau de puissances adressable régulièrernent, il peut être prévu de stocker, pour chaque facteur premier, autant de valeurs que le nombre de puissances de P2 nécessaires (p2=3). Pour les puissances non utilisées, on peut, bien entendu, utiliser un stockage de valeurs factices telles que des 1 ou des O. On décrit maintenant le calcul du rang d'une permutation pour réaliser un codage en exploitant l'invention. Il existe plusieurs variantes selon le jeu des décompositions de base choisi et leur représentation. Par souci de concision, l'exposé des réalisations possibles ci-après est limité au cas du mode de réalisation préféré pour le jeu des décompositions de base, avec décompositions des factorielles de 1! et de 1. Ci-après, la solution de représentation éclatée des exposants avec mn, exposants par décomposition de base, qui est le cas le plus général, est io d'abord exposée. Des variantes de représentation éclatée des exposants sont ensuite décrites. Enfin, la solution de représentation compacte des exposants des décompositions de base est exposée, ainsi que quelques unes de ses variantes. Il apparaîtra alors que l'invention s'applique parfaitement à un traitement de codage du rang d'une permutation. 15 L'algorithme de Schalkwijk est repris ci-après, comme exemple de traitement d'énumération. Représentation éclatée des exposants de la décomposition Soit n la dimension maximale des codes à permutation utilisés, et Mn! le nombre de facteurs premiers intervenant dans la décomposition de la 20 grandeur n! On décrit ci-après un premier mode de réalisation d'un codage utilisant une représentation éclatée des exposants de la décomposition. Codage selon un premier mode de réalisation Ici, les exposants des décompositions de base de 1 et 1! sont stockés 25 préférentiellement selon la "quatrième" solution de la table 4d ci-avant, avec stockage factice pour 1=0 dans deux tableaux monodimensionnels notés respectivement D, et Dä ayant mn, x (n +1) éléments. Comme mentionné ci-avant, on pourrait considérer aussi des tableaux bidimensionnels ayant m,,, colonnes et (n +1) lignes. Les exposants de 1(respectiivement 1!) étant stockés régulièrement (chacun sur mn, valeurs), les opérations de lecture des exposants d'une décomposition de base nécessitent un calcul d'adresse dans le tableau D, (respectivement Dä). Pour lire les exposants de la décomposition de 1! (respectivement 1), il faut pointer à l'adresse (1 x mm) du tableau Dä (respectivement D,), et viser ainsi l'adresse de l'exposant el, (respectivement 4), l'exposant e; (respectivement el) étant à l'adresse suivante (l x m,,, +1) et io plus généralement l'exposant el, (respectivement e,) étant à l'adresse (1 x m,,, + i -1). Comme mentionné ci-avant, dans des tableaux bidimensionnels, l'exposant el', (ou e;) est à l'adresse ((l;(i -1)) (colonne (i- 1) et lignel). Il est à noter que s'il n'est pas prévu de stockage factice pour 1=0, le calcul 15 d'adresse dans le tableau D, des n décompositions de base de 1 (1 >0) est: (1-1)xmn,. Initialisation • On initialise à zéro les mn, exposants de la décomposition intermédiaire 20 de Pk (stockés dans un tableau P ayant mn, éléments qui est préférentiellement mis à jour à chaque position, comme on le verra en référence à l'étape C-3 ci-après). Ainsi, les instructions sont les suivantes : o P[i]=0, 1 Parallèlement, à l'étape C-4, on calcule Sk à partir de la relation habituelle dk-I Sk = E w[d ] d=o L'étape C-5 est un test sur la valeur de Sk . Si Sk est nulle (flèche O), ce qui implique que le rang partiel Ikk est nul (formule (7) ci-avant) et que le rang t s n'est pas modifié, le traitement se poursuit par l'étape C-11 plus loin. Sinon (flèche N avec Sk ≠ 0 ), le traitement se poursuit par l'étape C-6, à laquelle on lit les mä, exposants e'k de la décomposition de base Sk dans le tableau D, à l'adresse mä, x Sk. Parallèlement, l'étape C-7 consiste à lire les mä, exposants e(ä_1_k), de la to décomposition de base (n - l - k)! dans le tableau D,1 à l'adresse mä, x (n -1- k). On relèvera que l'étape C-7 est effectuée si la somme Sk est non nulle (sortie N du test C-5) de manière à éviter la lecture inutile dans le tableau Dä si, de toutes façons, le rang partiel Ikk est nul. A l'étape C-8, les résultats des différentes lectures des tables peuvent être is regroupés pour calculer d'abord les mä, exposants de la décomposition du rang partiel Ikk , selon la relation : el/kk = + esk - P[i], pour 1 mä, Enfin, à l'étape C-9, on calcule le rang partiel Ikk lui-même par la relation : m, ee k Ikk =]J(P;) fk i=1 On rappelle que le terme w[dk ] est un poids nécessairernent inférieur ou égal à 20 la dimension maximale n du code à permutation considéré. De même, la somme Sk de tels poids reste inférieure à la dimension maximale n et il en est de même, bien entendu, pour (n-1-k). Les décompositions de w[dk ] , Sk et (n-1-k)! sont bien répertoriées dans les tables de décompositions des entiers ou des factorielles d'entiers allant jusqu'à la dimension maximale n, telles que la table 4d. A partir de la décomposition w[dk] répertoriée dans une table et de la décomposition de Pk_1 déterminée à la boucle précédente (k-1) et gardée en mémoire, la décomposition de Pk est déterminée. A titre d'exemple purement illustratif des étapes de calcul de rang partiel de la figure 5, on considère un code à permutation de dimension n= 8 et q= 4. Dans cet exemple, les exposants sont stockés dans un tableau bidimensionnel à quatre colonnes et neuf lignes, comportant donc 36 éléments. Il peut donc être extrait de la table 4d donnée ci-avant où on a p; = 2, 3, 5, 7 en colonnes et 1 = 0, ..., 8 en lignes. Dans cet exemple, on suppose qu'à une position précédente k=3, le tableau des poids w est {1, 1, 0, 3} et donc P3= 1! 1 0! 3! = 6. Le tableau P des exposants de la décomposition de P3 (=21x3'x5 x7 ) est donc {1, 1, 0, 0}. A la position k=2, pour l'étape C-1, on suppose qu'on a lu d2= 2. Dans cet exemple, à l'étape C-2, on met à jour l'élément w[2] en l'incrémentant de 1 (w[2] = 0 +1 = 1). A l'étape C-31, on lit les quatre exposants de la décomposition de 1 (= w[2] ), soit 0, 0, 0, 0 (voir dans la table 4d, 6éme à 9éme colonnes et Sème ligne 1=1). Puis (étape C-32), on met à jour le tableau P, on obtient donc P= {1, 1, 0, 0}. A l'étape C-4, on calcule Sk : Sk = w[0] + w[1] = 1+ 1 = 2. Sk n'est donc pas nul (test C-5). o On lit (étape C-6) les quatre exposants de la décomposition de Sk (en se référant toujours à la table 4d, 6ème à Sème colonnes mais 4éme ligne 1=2): 25 pour pi =2 (6éme colonne), l'exposant est 1, et pour p, =3, 5, 7 (7éme à 9éme colonnes), l'exposant est nul. 15 o On lit (étape C-7), les quatre exposants de la décomposition de 5! (correspondant à (n-1-k)!) en se référant toujours à la table 4d, mais cette fois aux 2ème à 5ème colonnes et 7ème ligne 1=5. Pour pi =2 (2ème colonne), l'exposant est 3. Pour p2 =3 (3ème colonne), l'exposant est 1 et pour p3 (4ème colonne), l'exposant est aussi 1. En revanche, pour p4 =7 (5ème colonne), l'exposant est nul. o A l'étape C-8, les quatre exposants de la décomposition du rang partiel Ikk sont calculés: o pour pi =2, el, =e(äù1ùk)! +e sk ùPU] = 3 + 1 -1 = 3 k o pour p2 =3, e1 k,, = + esk ù P[2] = 1 + 0 -1 = 0 o pour p3 =5, el3kk = + esk - P[3] = 1 + 0 -0 =1 o pour p4 =7, eÎk,Ik = + esk ù P[4] = 0 -F 0 -0 = 0 A l'étape C-9, le rang partiel Ikk est calculé à partir des exposants de sa décomposition: 23x3 x51 x7 =40. En référence à nouveau à la figure 5, le rang t global de la permutation est lui-même déterminé par mise à jour (avec addition du rang partiel Ikk trouvé à l'étape C-9) selon la relation t = t + Ik`k , à l'étape C-10. L'étape C-11 vise ensuite la décrémentation de la variable k (k=k-1) et l'étape 20 C-12 est un test sur la valeur de k pour poursuivre ou non le traitement. Ainsi, si k >. 0 (flèche O en sortie du test C-12), on réitère les étapes du traitement à partir de la première étape C-1, avec une valeur de k décrémentée d'une unité. Sinon (flèche N en sortie du test C-12), le traitement est terminé à l'étape C-13 de fin ("END"). 25 Ainsi, on comprendra que l'étape C-9 ci-avant calcule un rang partiel à partir de sa décomposition déterminée à l'étape C-8, elle-même déterminée à partir de trois décompositions intermédiaires : (n-1-k)! 5 Sk , et Pk . La détermination de deux d'entre elles, (n -1- k)! et S'k , effectuée aux étapes C-6 et C-7, consiste en une simple lecture dans les tableaux respectifs des décompositions de base D, et D,,. La détermination de la troisième 10 décomposition intermédiaire (celle de Pk) peut aussi être effectuée simplement à partir des q décompositions de base de w[d]! lues dans le tableau Dä selon q-1 une relation du type P[i] = ,[di, , 1 i 5 mn,. d=0 Variantes immédiates 15 L'étape C-3 présente une variante plutôt avantageuse de détermination de cette décomposition intermédiaire. La décomposition intermédiaire de Pk peut en effet être déterminée à partir d'une décomposition de base lue dans le tableau D, et d'une autre décomposition intermédiaire (celle de Pk+,) calculée pour un autre rang partiel (Ik+, ), par exemple le rang partiel Ik+, précédent 20 dans l'itération de k. Plus généralement, dans une variante de l'invention, une décomposition intermédiaire peut aussi se déduire d'au moins une décomposition intermédiaire précédemment déterminée pour au moins un autre rang partiel. Ci-avant, on a effectué le calcul du rang par une boucle de la dernière position (k=n-1) à la première position (k=0). Toutefois, l'invention s'applique aussi, bien entendu, à une boucle de la première à la dernière position. Il suffit de changer la phase d'initialisation et d'adapter les étapes C-2 et C-3 et leur ordre. A cet effet, le tableau des poids w peut être initialisé avec les q poids wô . Puis, pour d allant de 0 à q-1, on lit les mn, exposants de la décomposition de w[dJ! dans le tableau D , pour calculer par addition cumulative les mn, valeurs du tableau de la décomposition intermédiaire (exposants de la décomposition de Po). to L'étape C-3 intervient alors après l'étape C-10 et avant l'étape C-2. La décomposition intermédiaire de P est mise à jour par soustraction à P[i] de l'exposant ew[di de la décomposition de base de w[dk ] (P[i] = P[i] û ew[dk] , pour 1≤i≤ mn, ). Ensuite, l'étape C-2 peut être menée en décrémentant de 1 la valeur de w[dk] (w[dk ] = w[dk ] -1). L'étape C-11 peut viser l'incrémentation de 15 la variable k de 1 et l'étape C-12 teste simplement si k=n. On indique brièvement que, pour les codes à permutation de dimension n variable, plutôt que d'effectuer mnmar, opérations, il peut être préférable de lire mn, à l'initialisation pour n'effectuer les opérations que mn, fois. Variantes générales du premier mode de réalisation Plus généralement, la mise en oeuvre représentée sur la figure 5 pour un codage au sens de l'invention admet de nombreuses variantes. Ainsi, dans une première variante, chaque décomposition de base (de 1 ou 1!) 25 contient aussi le nombre m,, . La lecture du nombre m,, (0 ≤1≤n ) présente 20 20 des avantages. En effet, les étapes C-3et C-6 à C-10 ne sont plus effectuées mn, fois chacune, mais seulement : o mw[dk), fois pour l'étape C-3, o msk, fois pour l'étape C-6 et pour l'addition +esk de l'étape C-8, et 5 o m(n_I_k), fois pour les étapes C-7 et C-9 et pour la soustraction e'jkk = e(n_l_k), + esk - P[i] de l'étape C-8. Si, en outre, on a stocké aussi les valeurs m, , il n'y a plus qu'à effectuer : o mw[(k] fois l'étape C-3, et o msy fois l'étape C-6 et l'addition de l'étape C-8. Dans une autre variante du codage, si on utilise en plus le stockage des exposants de la décomposition de base de 1 selon la troisième solution (table 4c ci-avant), l'étape C-3 peut être effectuée pour m'W[dk] valeurs. De même, l'étape C-6 est effectuée pour m'sk valeurs. A l'étape C-8, au lieu de prévoir 15 m(n_1_k), additions et m(n-1-k)! soustractions, le nombre de soustractions à effectuer reste m(n_1_k), mais seulement m'sk additions sont nécessaires. En particulier : ^ ` el,k,k =e(n-t-k),ùePk` , pour1im (n-l-k)! ^ et e) = ekk) + e k~) , pour 1 s 0) est à l'adresse (1û1) 5 dans le tableau D,. Codaqe selon un second mode de réalisation Initialisation to On initialise à zéro le mot ep contenant la représentation compacte des mn, exposants de la décomposition intermédiaire de Pk : ep = 0 Le mot ep sera mis à jour à chaque position (étape CC-3 ci-après). • Comme précédemment, on initialise aussi à la valeur zéro le rang t et les q poids wk stockés dans un tableau w ayant q éléments qui sera mis à 15 jour à chaque position (étape CC-2 ci-après). Les instructions correspondantes sont : o tO o w[i]=O, 0Si notations des étapes CC-n de la figure 6 désignent le mot "codage" avec une "représentation compacte". A l'étape CC-1, on lit la variable dk . L'étape CC-2 consiste en une mise à jour de la variable w : w[dk ] = w[dk ] +1. L'étape CC-3 est une mise à jour du mot ep avec, en particulier : - à l'étape CC-31, la lecture du mot ew[dk] contenant la représentation compacte des exposants de la décomposition de w[dk ] dans le tableau el puis, - à l'étape CC-32, la mise à jour proprement dite du mot ep = ep +ewIdkj dk-1 to Parallèlement, à l'étape CC-4, on calcule la somme Sk: Sk = Ew[d]. L'étape d=o suivante CC-5 est un test sur la valeur de Sk . Si la somme Sk est nulle (flèche N), on décrémente directement l'indice k. Sinon (flèche O en sortie du test CC-5), le traitement se poursuit par la lecture du mot eSk contenant la représentation compacte des exposants de la décomposition de Sk dans le 15 tableau D; à l'étape CC-6. Parallèlement (et avantageusement selon le résultat du test CC-5), on lit le mot e(n-l-k), contenant la représentation compacte des exposants de (n -1- k)! dans le tableau Diu , à l'étape CC-7. A l'étape CC-8, Les résultats issus des différentes étapes CC-3, CC-6, CC-7 20 sont regroupés pour calculer le mot ef,k contenant la représentation compacte de la décomposition du rang partiel Ikk , par deux opérations simples (préférentiellement une addition suivie d'une soustraction) : e1kk = e(n-1-k), + eSk pins e1k, ,=e1Ak ùep L'étape CC-9 vise l'extraction des mn, exposants elJk de la décomposition du rang partiel Ikk et contenus dans le mot eldk , par décalage approprié des bits du mot comme décrit précédemment. A cet effet, il est prévu une boucle sur un indice i allant de 1 à mn, (initialisation CC-91 de i à 1, test CC-93 sur la valeur de i et incrémentation CC-94 tant que i n'atteint pas la valeur de mn, ). On applique alors sur chaque boucle de i des instructions CC-92 du type : i'l) e/kk =elkk &((1 bn)ù1) i'2) elk dk =e~k bn ~k On rappelle que les notations " b" et " b" désignent respectivement un io décalage à gauche et un décalage à droite de b bits. Par ailleurs, la notation "&" désigne l'opérateur logique bit à bit "ET". L'instruction i'l) consiste à récupérer les b' bits de poids faible de eÎdk . En pratique, les masques correspondant à l'instruction "((1 bn)ù1)" sont câblés en dur. En d'autres termes, on applique au départ (pour l'indice de boucle i égal à 1) 15 un masque ((1 b)ù1) des bits de poids fort de manière à ne récupérer d'abord que les bä bits de poids faible de efdk donnant l'exposant el,,, associé au premier facteur premier pi. Puis : - on décale les bits de et:, de bn "à droite" pour récupérer les bits de poids 20 plus fort donnant les exposants associés au facteur premier suivant P2 (masque ((1 13,2) -1) des bits de poids fort, - on extrait alors l'exposant e12.,k - puis, on applique un décalage à droite de bä2 bits, et ainsi de suite jusqu'à i= mn, . L'étape suivante CC-10 consiste à calculer le rang partiel Ikk comme suit : n dk Ikk =n(i/)' Puis, on additionne le rang partiel Ikk au rang total t (t=t+Ikk à l'étape 5 CC-11). La valeur de l'indice k est décrémentée à l'étape suivante CC-12 (k=k-1) et, avant de recommencer les étapes CC-4, CC-1, CC-7 et suivantes avec cette valeur décrémentée, on vérifie au test CC-13 si la valeur de k n'a pas atteint -1 (k 15 termes (n -1- k)!, Sk et Pk) sont utilisées. La détermination de deux d'entre elles n-1- k)! et Sk) effectuées aux étapes CC-6 et CC-7 consiste en une simple lecture de leur représentation compacte dans les tableaux Diu et D'1. La détermination de la troisième décomposition intermédiaire (Pk) effectuée à l'étape CC-3 nécessite aussi une lecture du tableau D, suivie d'une mise à jour 20 de la représentation compacte de cette décomposition intermédiaire par addition de la représentation compacte de la décomposition de base qui a été lue. Comme indiqué précédemment en référence au premier mode de réalisation, le stockage des valeurs m,, (0 1 ≤n) permet de réduire la complexité des 25 étapes CC-9 et CC-10. La boucle d'extraction des exposants de la décomposition d'un rang partiel Ikk est effectuée m(ä_1_k), fois au lieu de m,,, 67 fois. De même, l'étape CC-10 ne comporte plus mn, multiplications mais m(n_,_k), multiplications. On décrit maintenant le décodage du rang d'une permutation en 5 exploitant l'invention. II existe, là encore, plusieurs variantes selon la solution de représentation des décompositions de base (éclatée ou compacte). On décrit ci-après un premier mode de réalisation du décodage, homologue du premier mode de réalisation décrit ci-avant pour le codage utilisant la représentation éclatée des to décompositions et leur stockage selon la quatrième solution relative à la table 4d donnée plus haut. Il apparaîtra que l'invention s'applique parfaitement à un décodage du rang d'une permutation en prenant comme exemple l'algorithme de Schalkwijk. 15 Décodage selon un premier mode de réalisation Le décodage utilisant une représentation éclatée des exposants de la décomposition commence préférentiellement par une initialisation des données comme suit. Initialisation 20 • On initialise le tableau des valeurs w avec les q poids wâ (w est mis à jour en fin de boucle sur chaque position de l'étape D-19 qui sera décrite ci-après). Une instruction appropriée peut être du type : o w[i]=wo, OSi Les trois étapes suivantes visent un test sur la valeur du rang total t, en le comparant à la valeur du rang partiel. Pour ce faire, à l'étape D-12, on incrémente la valeur de dk (dk= dk+1) et le test D-13 est le suivant : 15 tùIkk 0 ? Si l'inégalité est vérifiée (flèche 0), on réitère les étapes D-7 à D-13 avec la nouvelle valeur de dk incrémentée. Sinon (flèche N), le traitement se poursuit par l'étape D-14 de décrémentation de la valeur de dk (dk= dk-1) pour revenir à la valeur de dk avant l'étape D-12. Pour cette valeur de dk, le rang partiel Ikk 20 prend la valeur intermédiaire hi précitée (étape D-15 : Ikk =1k ). Ensuite, on met à jour le rang t, lequel devient t = t - Ikk (étape D-16), et on met à jour les exposants de Pk (dans le tableau P) grâce à la lecture des m,, exposants ew[dki de la décomposition de w[dk ] dans le tableau D, (étape D-17), suivie de la mise à jour proprement dite des exposants P[i] = P[i] ù ew[dki pour 1 _ Une variante consistant à utiliser les exposants du rapport Rk (comme décrit précédemment) est particulièrement intéressante. En effet, dans le traitement de décodage décrit ci-avant en référence à la figure 7, on calcule pour une position donnée k, les exposants e'',,,, pour plusieurs valeurs de d. Pour 5 chaque valeur de d testée, les autres variantes nécessitent pour chaque exposant une soustraction et une addition (erdk = e~n ,_k~, ù epk + esk ). Cependant, la variante utilisant les exposants du rapport Rk ne requiert qu'une addition (e' = R[i]+e) car seuls la somme Sk et ses exposants esk varient en fonction de d pour une position k donnée. Décodaqe selon un second mode de réalisation On décrit ci-après, en référence maintenant à la figure 8, un exemple de réalisation d'un décodage utilisant la représentation compacte des décompositions. 15 Au préalable, on prévoit une initialisation des données comme suit. Initialisation • On se réfère d'abord au tableau w de ayant q éléments pour déterminer les termes w[i] = wô pour 0 i compact e Jak : k e'kk = e'kk bY! Cette étape générale DC-11 est similaire à l'étape générale CC-9 de la figure 6 pour le codage. Les étapes DC-12 à DC-17 sont, quant à elles, similaires aux étapes D-11 à io D-16 décrites ci-avant en référence à la figure 7 propre au décodage en représentation éclatée. La mise à jour des exposants de Pk (tableau P) à l'étape générale DC-18 se fait, en revanche, à l'étape DC-181, en lisant le mot ew[dk) contenant la représentation compacte des mh, exposants de w[dk ] dans le tableau D; et la 15 mise à jour proprement dite des exposants de Pk (e,, = ep ûew[Jk 1) s'effectue ensuite à l'étape DC-182. Ensuite, les étapes DC-19 à DC-22 sont similaires aux étapes D-19 à D-22 de la figure 7 propre au décodage utilisant une décomposition éclatée. 20 On décrit maintenant les différents avantages qu'offrent les variantes exposées ci-avant. Les variantes du premier mode de réalisation avec représentation éclatée, utilisant les tables de m,, (et/ou m, ou m', ), comportent moins d'opérations d'additions/soustractions que le mode de réalisation principal qui n'utilise que la table des valeurs mn, . Dans ce cas, le gain en complexité est important surtout pour les dernières positions (c'est-à-dire lorsque m(ä_k),, m'1 ou m, sont bien inférieurs à mn, ). Néanmoins, ce gain en complexité s'accompagne d'un accroissement de la complexité des étapes de lecture en mémoire (étapes C-31, C-6 et C-7). Si le nombre de valeurs à lire est moindre, le calcul d'adresse est en revanche plus compliqué. Un compromis intéressant consiste alors à stocker régulièrement les to décompositions de base (avec mn, exposants) pour faciliter l'adressage des tableaux D, et D,, et de stocker alors les valeurs m,, dans un tableau Dm à (n+1) éléments. II convient alors de stocker aussi les valeurs m, pour réduire efficacement le nombre d'additions/soustractions. Toutefois, cette mesure s'accompagne nécessairement d'une lecture des valeurs msk (respectivement 15 mwrdki), avant les étapes C-6 et D-9 (respectivement C-3 et D-19) tandis que la valeur m(n_k), ne doit être lue qu'au début de chaque itération sur k. Par ailleurs, les avantages apportés par la représentation compacte par rapport à la représentation éclatée sont les suivants : 20 - l'étape de mise à jour du tableau P ne comporte plus maintenant qu'une seule addition (respectivement soustraction) au codage (respectivement au décodage), - le calcul de l'exposant ef,,k ne nécessite aussi qu'une seule addition et qu'une seule soustraction, 25 - les calculs d'adresse pour la lecture des mots eK sont directs et ne nécessitent, pour chaque valeur K, qu'un seul accès et lecture mémoire. Par contre, la représentation compacte nécessite une extraction des exposants du rang partiel Ikk contenus dans le mot e/kk (étapes CC-9 et DC-11). Toutefois, cette opération n'est pas forcément un inconvénient pour le calcul du rang partiel à partir de sa décomposition en facteurs premiers, 5 comme on le verra ci-après. On décrit ci-après les avantages de telles variantes pour le calcul du rang partiel à partir de sa décomposition en facteurs premiers. La complexité de l'étape de calcul du produit des puissances des facteurs io premiers, au sens des étapes C-9 et CC-10 (respectivement D-11 et DC-12) au codage (respectivement au décodage) croît de façon importante avec le nombre de facteurs, même si elle reste bien inférieure à la complexité de la division au sens de l'art antérieur. Or, en pratique, beaucoup des exposants de la décomposition d'un rang partiel sont nuls et donc les puissances 15 correspondantes sont égales à 1. Souvent, tous les exposants sont nuls ou les premiers exposants seulement ne sont pas nuls. II est donc utile de pouvoir détecter et ne garder que les puissances d'exposants non nuls. Dans la représentation détaillée, cette détection ne peut s'effectuer que par mn, tests ou m(n_t_k), tests (un par facteur premier). 20 Avantageusement, la représentation compacte permet de tester avec un seul test si tous les exposants sont nuls (si e1dk = 0) et dans ce cas, le rang t'= 1. En outre, la détection du bit de poids fort de e, permet d'obtenir l'indice du plus grand facteur premier d'exposant non nul dans le rang t' et de réduire le nombre de répétitions de la boucle de l'étape CC-9 (respectivement DC-11) au 25 codage (respectivement au décodage). II est à noter néanmoins que la détection des exposants non nuls, en représentation détaillée comme en représentation compacte, augmente la complexité. Si tous les exposants sont non nuls, la complexité de la multiplication des puissances des facteurs premiers reste la même et s'ajoute alors à cette complexité celle de la procédure de détection des exposants non nuls. Ainsi, dans une autre variante, la détection des exposants nuls peut être menée uniquement si le nombre possible de facteurs premiers devient grand (k bien inférieur à n) et que la complexité de la multiplication de leurs puissances est supérieure à la complexité de la procédure de détection. Pour ce faire, des boucles différentes selon les positions peuvent être prévues, to même si cette implémentation s'effectue au détriment d'un accroissement des lignes d'instruction. II est aussi possible de combiner les représentations éclatée et compacte. Pour les dernières positions (la valeur m,! étant faible), le calcul des décompositions intermédiaires nécessite peu d'opérations. On favorisera alors is l'utilisation d'une représentation éclatée qui ne nécessite pas d'extraction des exposants d'un rang partiel. En revanche, pour les premières positions, on favorisera plutôt l'utilisation d'une représentation compacte. On décrit maintenant quelques exemples de réalisation dans des 20 codeurs/décodeurs existants. Codeur 3GPP AMR-WB+ Le codeur 3GPP AMR-WB+ (norme [3GPPTS26.304]) utilise une quantification vectorielle algébrique dont le dictionnaire est une union des codes à permutation du réseau de Gosset RE8 de dimension 8. 25 La technique TCX correspond à un codage prédictif par transformée. Plus précisément, il s'agit d'une méthode de codage par transformée FFT appliquée après filtrage de pondération perceptuelle. Dans la norme [3GPPTS26.304], le spectre FFT obtenu est découpé en sous-bandes (ou sous-vecteurs) de dimension n=8 et ces sous-vecteurs sont codés séparément. La quantification des sous-vecteurs utilise le réseau régulier de points RE8. Les dictionnaires de quantification en dimension 8 sont composés d'une union de codes à permutations de type I issus du réseau de points RE8. Dans le codeur TCX selon la norme [3GPPTS26.304], chaque code à permutation correspond à un vecteur leader signé donné en dimension n=8. L'indice de quantification d'un point du réseau RE8 est calculé par une formule du type: indice = offset de cardinalité + rang de la permutation w Le rang est calculé par la formule de Schalkwijk tandis que l'offset de cardinalité est tabulé. Néanmoins, ces leaders signés sont représentés par l'intermédiaire de leurs leaders absolus afin d'optimiser le stockage et la recherche dans les codes à permutations. On trouve la liste des leaders absolus associés dans la référence : 15 "Low-complexity multi-rate lattice vector quantization with application to wideband TCX speech coding at 32 kbit/s" par Ragot S., Bessette B., Lefebvre R., dans Proc. ICASSP, vol. 1, mai 2004 pp. 501-4. Pour illustrer les différentes variantes de l'invention, trois exemples de 20 réalisation sont présentés ci-après. Les deux premiers exemples de réalisation concernent le calcul du rang d'une permutation (codage), l'un utilisant une représentation éclatée des décompositions et l'autre une représentation compacte. Dans ces exemples de réalisation ci-après et les annexes correspondantes, 25 les tableaux R et P sont indexés de R[O] à R[mn, -1] et de P[0] à P[mä, -1] (et non de 1 à mn, comme décrit à titre d'exemple ci-avant), et ce, sans aucune incidence particulière sur le traitement pour le calcul du rang. Premier exemple de réalisation (codage) Dans ce mode de réalisation, on utilise une représentation éclatée des décompositions de base. On stocke leurs exposants dans deux tableaux à 36 éléments (=(8+1)x4). Il s'agit des tableaux donnés en annexe A-11 et notés D, [36] (contenant les exposants des décompositions des entiers 1 (0 15 intermédiaire du quotient : Rk = e-, fl(wk!) =o Comme indiqué plus haut, plutôt que de déterminer la décomposition intermédiaire de Rk à partir de décompositions de base correspondant aux (q+1) décompositions de base de (7-k)! et ( _, , on détermine cette 20 décomposition intermédiaire à partir de la décomposition intermédiaire de Rk+, et des deux décompositions de base de (7-k) et de w. Les quatre exposants de cette décomposition intermédiaire sont stockés dans un tableau (7 - k)! R. to k (7-k)! max(wk) max(mWk ,) max(Sk) max(ms ) 7 0 1 1 1 0 0 6 1 1 2 1 1 1 2 1 3 2 2 1 4 3 2 4 2 3 2 3 4 2 5 3 4 2 2 5 3 6 3 5 3 1 6 3 7 4 6 3 0 7 4 8 4 7 4 Table 11 : Maximum de m, ou de m,, pour les décompositions des termes d'un rang partiel à la position k et pour les codes à permutation de 5 dimension 8 L'invention exploite la connaissance de m(,_kyet des maxima de mW; et de ms, k indiqués pour chaque position dans la table 11 ci-dessus pour ne pas calculer les exposants des facteurs premiers d'indices plus grands que ces bornes. Le traitement correspondant est donné en annexe A-12. On note que la boucle sur les positions est éclatée. On note aussi que l'exposant du facteur premier pi du quotient est stocké dans l'élément R[i-1] du tableau R à 4 éléments. 15 Deuxième exemple de réalisation (codage) Dans une variante avec le codeur 3GPP AMR-WB+, les décompositions de base sont représentées de façon compacte. On stocke les mots contenant leurs exposants dans deux tableaux à 9 éléments (=(8+1)). En se référant à l'annexe A-21, le tableau D; contient les mots pour les décompositions des entiers 1(0 ≤.1≤ 8) (donc avec stockage factice de la décomposition pour 1=0) et le tableau D,, contient les mots pour les décompositions de leurs factorielles. On stocke aussi les puissances de 3, 5 et 7 dans un tableau Pow[12] à 12 éléments (avec stockage factice de 0 pour les puissances non utilisées). La décomposition d'un rang partiel est déterminée à partir de trois décompositions intermédiaires, deux étant les décornpositions de base de l'entier Sk et de la factorielle (7-k)! et la troisième étant une décomposition io intermédiaire du dénominateur du rang partiel : y-1 Pk =fl(wk!) ,=0 Comme indiqué précédemment, plutôt que de déterminer la décomposition intermédiaire de Pk à partir des q décompositions de base de (wk!),=0 9_1, on détermine cette décomposition à partir de la décomposition intermédiaire de 15 Pk+1 et de la décomposition de base de wfk' . Le mot compact contenant les quatre exposants de cette décomposition intermédiaire est noté "eP" dans l'annexe A-22. En outre, on note "el" le mot compact contenant les quatre exposants de la décomposition d'un rang partiel. Ici encore, on exploite la connaissance de m(7_k), pour n'extraire que m(7_k), 20 exposants du mot compact représentant la décomposition d'un rang partiel. Le traitement correspondant fait l'objet de l'annexe A-22. Là encore, la boucle sur les positions est éclatée. Troisième exemple de réalisation (décodage) Le troisième exemple de réalisation traite du décodage du rang d'une permutation, en codage 3GPP AMR-WB+. On utilise préférentiellement une représentation éclatée des décompositions de base comme dans le premier exemple de réalisation et une décomposition du rang partiel en trois termes comme dans le second exemple de réalisation. La boucle sur les positions n'est pas éclatée toutefois. Comme indiqué précédemment, plutôt que de déterrniner la décomposition intermédiaire de Pk à partir de décompositions de base, on la détermine à io partir d'une décomposition intermédiaire de Pk_1 et de la décomposition de base de wkk . Les quatre exposants de cette décomposition intermédiaire sont stockés dans un tableau P. De même, à partir de la décomposition intermédiaire de Pk et de la décomposition de base de (7-k)!, on calcule une autre décomposition intermédiaire (celle du quotient) dont les exposants sont 15 stockés dans un tableau R. Le traitement correspondant fait l'objet de l'annexe A-3. On note que l'exposant du facteur premier pi du quotient (respectivement du produit) est stocké dans l'élément R[i-1 ] (resp. P[i-1 ]) du tableau R (resp. P) à quatre éléments. 20 Ainsi, le premier exemple ci-avant utilise une décomposition du rang partiel en deux termes (dont un quotient), les deux autres exemples utilisent une décomposition en trois termes (deux pour le numérateur et un pour le dénominateur). Le mode pour le décodage n'utilise que m8,(=4) termes tandis 25 que les deux modes pour le codage utilisent un traitement séparé des positions pour exploiter les ml, ou m, termes qui ne sont pas lus mais inscrits "en dur" dans l'algorithme en éclatant la boucle sur les 8 positions d'un point du réseau de Gosset. Exemple de réalisation pour le codeur TDAC Un dernier exemple de réalisation concerne le codeur fréquentiel perceptuel TDAC de la Demanderesse utilisé pour coder des signaux audio numériques 5 échantillonnés à 16 kHz (bande élargie), dont le principe est décrit ci-après. Le codeur TDAC utilise une quantification vectorielle statistique à dimension et résolution variables de dimension maximale 15. Dans le cas des codes à permutation du réseau régulier de points RE8 de 10 dimension 8, l'invention permet essentiellement de réduire la complexité. Par contre, dans le cas du codeur TDAC, qui utilise des codes à permutations de dimension supérieure à 12, l'invention s'avère très avantageuse car elle permet non seulement une réduction de complexité mais aussi une réalisation du codeur sur des processeurs en virgule fixe dont la précision maximale est 15 limitée à des entiers de 32 bits non signés. Sans l'invention, une telle réalisation serait extrêmement complexe. Le principe de ce codeur est le suivant. Un signal audio limité en bande à 7 kHz et échantillonné à 16 kHz est découpé 20 en trames de 320 échantillons (20 ms). Une transformée en cosinus discrète modifiée (ou "MDCT") est appliquée sur des blocs du signal d'entrée de 640 échantillons avec un recouvrement de 50 % (ce qui correspond à un rafraîchissement de l'analyse MDCT toutes les 20 ms). On limite le spectre à 7225 Hz en fixant à zéro les 31 derniers coefficients (seuls alors les 289 25 premiers coefficients étant différents de 0). Une courbe de masquage est déterminée à partir de ce spectre et tous les coefficients masqués sont mis à zéro. Le spectre est divisé en 32 bandes de largeurs inégales. Les éventuelles bandes masquées sont déterminées en fonction des coefficients transformés des signaux. Pour chaque bande du spectre, l'énergie des coefficients MDCT 30 est calculée (pour évaluer des facteurs d'échelle). Les 32 facteurs d'échelle constituent l'enveloppe spectrale du signal qui est ensuite quantifiée, codée et transmise dans la trame. L'allocation dynamique des bits se base sur une courbe de masquage par bande calculée à partir de la version déquantifiée de l'enveloppe spectrale, de manière à obtenir une compatibilité entre l'allocation binaire du codeur et du décodeur. Les coefficients MDCT normalisés dans chaque bande sont ensuite quantifiés par des quantificateurs vectoriels utilisant des dictionnaires imbriqués en taille, les dictionnaires étant composés d'une union de codes à permutation de type Il. Finalement, les informations sur la tonalité et le voisement ainsi que l'enveloppe spectrale et les coefficients w codés sont multiplexés et transmis en trame. L'exemple de réalisation pour le calcul du rang d'une permutation (codage) utilise ici une représentation compacte des décompositions. La dimension des codes à permutation utilisée étant variable, la boucle sur les positions n'est 15 pas éclatée. Ce mode de réalisation illustre une procédure de détection des exposants non nuls de la décomposition d'un rang partiel. Ici, les décompositions de base sont représentées de façon compacte. On stocke les mots contenant leurs exposants dans deux tableaux à seize éléments (=(15+1)). Dans l'annexe B-1, Le tableau D; contient les mots pour 20 les décompositions des entiers 1(0 ≤l Bien entendu, la présente invention ne se limite pas aux formes de réalisation 30 décrites ci-avant à titre d'exemple ; elle s'étend à d'autres variantes. A la connaissance de la Demanderesse, la présente invention consiste en la première utilisation des décompositions en puissances de facteurs premiers dans les codes à permutation. Cette utilisation est pourtant particulièrement avantageuse dès lors que des calculs d'expressions combinatoires sont prévus, comme dans la quantification vectorielle avec codes à permutation. Ainsi, de manière générale, la présente invention vise cette utilisation des décompositions en puissances de facteurs premiers pour toute expression combinatoire, même différente d'un rang d'une permutation, en to codage/décodage par un ou des codes à permutation. La présente invention trouve une application avantageuse en codage/décodage de signaux de parole, par exemple clans des terminaux de téléphonie notamment cellulaire. Toutefois, elle s'applique au 15 codage/décodage de signaux de tout autre type, notamment de signaux d'image ou vidéo ainsi qu'en modulation codée. La présente invention vise aussi un programme informatique destiné à être stocké en mémoire d'un dispositif de codage/décodage de signaux 20 numériques, utilisant des codes à permutations. Ce programme comporte alors des instructions pour la mise en oeuvre des étapes du procédé au sens de l'invention. Typiquement, les figures 4 à 8 décrites ci-avant peuvent correspondre à des organigrammes d'algorithmes que peut comporter un tel programme. 25 La présente invention vise aussi un dispositif de codage/décodage de signaux numériques, utilisant des codes à permutations, et comprenant, en référence à la figure 9 : to 15 - une unité mémoire MEM pour stocker les instructions d'un programme informatique du type précité, ainsi que des représentations préenregistrées de décompositions de nombres entiers choisis, et - un module de calcul PROC accédant à cette unité mémoire MEM pour la 5 mise en oeuvre du procédé au sens de l'invention. Ces moyens MEM, PROC peuvent être prévus : • pour donner le rang t d'une permutation à partir du vecteur-code y choisi (flèches en traits pleins de la figure 9) : o dans un module d'indexation d'un codeur de source ou o dans un module d'indexation d'un décodeur de canal, • ou pour donner le vecteur-code y à reconstruire à partir du rang t d'une permutation (flèches en traits pointillés de la figure 9) : o dans un module de décodage d'un décodeur de source ou o dans un module codage d'un codeur de canal. Bien entendu, les représentations préenregistrées dans la mémoire MEM peuvent être sous forme de contenus d'adresses (représentation éclatée) ou sous la forme de mots de bits (représentation compacte). ANNEXES, s A-11 = { D, [36] 0, 0, 0, 0, 0, 0, 0, 0, 1,0,0,0, 0, 1, 0, 0, to 2, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 3, 0, 0, 0}; 15 Dä [36] = { 0, 0, 0, 0, /* 0! */ 0, 0, 0, 0, /* 1! */ 1, 0, 0, 0, /* 2! */ 20 1, 1, 0, 0, /* 3! */ 3, 1, 0, 0, /* 4! */ 3, 1, 1, 0, /* 5! */ 4, 2, 1, 0, /* 6! */ 4, 2, 1, 1, /* 7! */ 25 7,2, 1, 1,/*8!*/ }; A-12 Initialisation : • t=0 5 • w[i]=O, 0_ d7) t= 1 15 • Position k=5 : o d=d[5]; w[d] = w[d] +1 o R[O] = R[O]+D,[4*2J- D,[4*w[d]] o R[1] = -D,[4*w[d]+1] o S=O; boucle de j=O à d-1: S = S + w[j] 20 o Test si S>O ; si oui : ^ i2= R[O]+ D,[4*S] ^ t=t+(1 i2) • Positions k=4 puis 3 : o d=d[k]; w[d] = w[d] +1 o R[O] = R[O]+Di[4*(7-k)]- D,[4*w[d]] o R[1] = R[1 ]+D,[4*(7-k)+1l- D,[4*w[d]+1] o S=O; boucle de j=O à d-1: S = S + w[j] o Test si S>0; si oui : R[2] = ^ i2= R[O]+ D,[4*S] • ^ i3= R[1 ]+ D,[4*S+1 ] • ^ t= t + pow3[i3] *(1 i2) -D,[4*w[d]+2] Positions k=2 puis 1 : o d=d[k]; w[d] = w[d] +1 o R[O] = R[O]+D,[4*(7-k)]- D,[4*w[d]] o R[1] = R[1]+D,[4*(7-k)+1]-D,[4*w[d]+1] o R[2] = R[2]+D,[4*(7-k)+2]- D,[4*w[d]+2] o S=0; boucledej=Oàd-1 : S=S+w[j] o Test si S>0; si oui : ^ i2= R[O]+ D,[4*S] ^ i3= R[1]+ D,[4*S+1 ] ^ i5= R[2]+ D,[4*S+2] ^ t= t + pow5[i5]*pow3[i3]*(1 i2) • R[3] = -D,[4*w[d]+3] • Position k=O : o d=d[O]; w[d] = w[d] +1 o R[O] = R[O]+D,[4*7]- D,[4*w[d]] o R[1] = R[1]+Di[4*7+1]-Di[4*w[d]+1] o R[2] = R[2]+D,[4*7+2]- D,[4*w[d]+2] o R[3] = R[2]+D,[4*7+3]- D,[4*w[d]+3] o S=0 ; boucle de j=O à d-1 : S= S+ w[j] o Test si S>0; si oui : ^ i2= R[O]+ D,[4*S] ^ i3= R[1]+ D,[4*S+1] ^ i5= R[2]+ D,[4*S+2] ^ i7= R[3]+ D,[4*S+3] to t= t+ pow7[i7]*pow5[i5]*pow3[i3]*(1 d7) t= 1 sinon si (d • Position k=5 : o d=d[5] ; w[d] = w[d] +1 o eP = eP+D,[w[d]] o S=0; boucledej=Oàd-1: S=S+w[j] o Test si S>O ; si oui : i2= D ir[2]+D I[S] ^ ^ i2= i2-eP • t= t+ (1 i2) • Positions k=4 puis 3 : o d=d[k] ; w[d] = w[d] +1 o eP = eP+D',[w[d]] o S=0 ; boucle de j=O à d-1 : S = S + w[j] o Test si S>0; si oui : ^ el = D ä[(7-k)]+D I[S] ^ el = el- eP i2 = et & (0x7) ^ ^ i3=el 3 ^ t= t + pow[i3] *Cl i2) • Positions k=2 puis 1 : o d=d[kw[d]=w[d]+1 o eP = eP+D,[w[d]] o S=0 ; boucle de j=O à d-1 : S= S+ w[j] o Test si S>0; si oui : ^ el = D i1[(7-k)]+Di[S] 20 93 ^ el = el- eP ^ i2 = el & (0x7) ^ el = el 3 ^ i3 = et &(0x3) ^ et=el 2 ^ i5 = el &(0x3) ^ t= t + pow[4+i5] * pow[i3] *Cl i2) • Position k=O : o d=d[O]; w[d] = w[d] +1 10 o eP = eP+D,[w[d]] o S=O; boucle de j=O à d-1 : S=S+w[j] o Test si S>0; si oui : ^ el = Dir[(7)]+D I[S] ^ et = el- eP ^ i2 = et & (0x7) ^ et=el 3 ^ i3 = el &(0x3) ^ et=el 2 ^ i5 = el &(0x3) ^ et=el 2 ^ i7 = el &(0x3) ^ t= t+ pow[8+i7] *pow[4+i5] * pow[i3] *(1 i2) A-3 Initialisa fion : o Boucle de i=0 à 3, P[i] =Dä[4*w[0]+i] o Boucle de d=1 à q-1: ^ Boucle de i=0 à 3, P[i] =P[i] + Dä[4*w[d]+i] Traitement des positions • Boucle de k=0 à 7 : o Boucle de 1=0 à 3: R[i] = Dä[4*(7-k)+i]-P[i] o 1=0 to o d=0 o Répéter tant que w[d] = 0, d=d+1 o R[0] = R[0]+D,[4*(7-k)]- D,[4*w[d]] o R[1] = R[1]+D,[4*(7-k)+1]D,[4*w[d]+1] o R[2] = R[2]+D,[4*(7-k)+2]-D,[4*w[d]+2] 15 o S=0 o Répéter tant que (t-l> 0) ^ 1,=1 ^ S=S+w[d] ^ 12= R[0]+ D,[4*S] 20 ^ 13= R[1]+ D,[4*S+1 ] ^ 15= R[2]+ D,[4*S+2] 17= R[3]+ D,[4*S+3] ^ 1= pow7[i7]*pow5[i5]*pow3[i3]*(1 dnl)t=1 • Positions de k=n-3 à 0 : o d=d[k] o w[d] = w[d] +1 o eP = eP+D,[w[d]] o S=0; boucle de j=O à d-1: S=S+w[j] o Test si S>0; si oui : ^ t'= 1 ^ el = D it[(7-k)]+D I[S] to ^ el = el- eP ^ Test si (eI>0); si oui • i2 = el & (OxF) • Pow[O] = 1; • mI=0 15 • i3 = (eI 4) & 0x7 • si (i3 >0) o pow[0] = Pow3[i3] o ml = ml +1 • el = el 7 20 • si et >0 mI'= (16-norm_s(eI)) 1 boude de j=0 à ml' ^ i=el&0x3 98 ^ Sii >0 ^ Sii >0 • pow[mI] = Pow[i+5*i] • ml = ml +1 ^ et=el 2 • boucle de i=0 à (ml-1): t' = t' x pow[i] • t' = t' i2 ^ t= t+t' Revendications 1. Procédé de codage/décodage de signaux numériques, utilisant des codes à permutation faisant intervenir un calcul d'expressions combinatoires, caractérisé en ce que lesdites expressions combinatoires sont représentées par des décompositions en puissances de facteurs premiers, et déterminées par une lecture mémoire de représentations préenregistrées de décompositions de nombres entiers choisis. to 2. Procédé selon la revendication 1, caractérisé en ce que les représentations préenregistrées comportent des valeurs représentatives d'exposants, stockées respectivement en correspondance de valeurs représentatives de nombres premiers successifs, pour chacun desdits entiers choisis. 15 3. Procédé selon l'une des revendications 1 et 2, caractérisé en ce que lesdites expressions combinatoires comportent des valeurs factorielles (1!) d'entiers (1), et en ce que lesdites représentations préenregistrées comportent au moins des représentations de décompositions de valeurs factorielles (1!). 20 4. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'une au moins desdites expressions combinatoires (t) comporte un quotient d'un numérateur entier (v) par un dénominateur entier (p), et en ce que ledit quotient est représenté par une décomposition de puissances de facteurs premiers, dont chaque puissance est une différence 25 d'exposants (e, - e, ) respectivement associés au numérateur (v) et au dénominateur (p) et affectés à un même nombre premier (pi). 99loo 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que la représentation préenregistrée d'un nombre entier choisi est stockée dans une mémoire adressable, et en ce qu'un adressage de ladite mémoire donne une succession d'exposants à affecter à des nombres premiers respectifs pour recomposer l'entier choisi. 6. Procédé selon la revendication 5, caractérisé en ce que la représentation préenregistrée d'un nombre entier choisi est stockée sous la forme d'une succession d'adresses donnant chacune, pour un nombre premier, un to exposant à affecter à ce nombre premier pour recomposer l'entier choisi. 7. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que les représentations préenregistrées sont stockées sous la forme de mots comportant une succession de groupes de bits, chaque groupe ayant : 15 -un poids fonction d'un nombre premier, et - une valeur fonction d'un exposant à associer à ce nombre premier. 8. Procédé selon la revendication 7, caractérisé en ce que la détermination des puissances de facteurs premiers s'effectue par applications successives 20 d'au moins un masque partiel au mot de bits, avec décalages successifs en fonction des poids des bits et lectures des bits restants. 9. Procédé selon la revendication 2, caractérisé en ce qu'il comporte, pour le calcul d'une expression combinatoire, les étapes : 25 -identification, parmi lesdits nombres entiers choisis, de termes ((n-1-k)!, Sk, Pk) apparaissant dans un produit et/ou un quotient formant ladite expression combinatoire (Ikk )i - lecture mémoire des exposants (e'(ä_l_k)!, e'sk, e'pk) intervenant dans les décompositions en facteurs premiers desdits termes, addition et/ou respectivement soustraction des exposants lus (e'(ä_I_k)!+e'sk-e'pk) pour déterminer les exposants (e~kk ) intervenant dans la décomposition en puissances de facteurs premiers de ladite expression combinatoire, et, de là, calculer ladite expression combinatoire à partir de sa décomposition en puissances de facteurs premier. 10. Procédé selon l'une des revendications 2 et 9, comportant une étape ro récurrente de calcul d'un produit faisant intervenir à chaque récurrence un terme se multipliant à un produit déterminé à une récurrence antérieure, caractérisé en ce que : - ledit produit déterminé à une récurrence antérieure est conservé en mémoire sous la forme d'une décomposition en puissances de facteurs premiers, 15 - ledit terme se multipliant au produit est l'un des entiers choisis dont la décomposition est préenregistrée, et - pour déterminer ledit produit à une récurrence courante, on additionne, un à un par nombre premier, les exposants issus des décompositions respectives dudit produit déterminé à une récurrence antérieure et dudit terme se 20 multipliant au produit. 11. Procédé selon l'une des revendications 2 et 9, comportant une étape récurrente de calcul d'une division faisant intervenir à chaque récurrence un terme divisant un quotient déterminé à une récurrence antérieure, caractérisé 25 en ce que: - ledit quotient déterminé à une récurrence antérieure est conservé en mémoire sous la forme d'une décomposition en puissances de facteurs premiers,2 - ledit terme divisant le quotient est l'un des entiers choisis dont la décomposition est préenregistrée, et -pour déterminer ladite division à une récurrence courante, on soustrait, un à un par nombre premier, les exposants issus de la décomposition dudit terme aux exposants issus de la décomposition dudit quotient déterminé à une récurrence antérieure. 12. Procédé selon l'une des revendications précédentes, caractérisé en ce que les codes à permutation font intervenir le calcul d'une quantité représentative to d'un rang d'une permutation (t) comportant un cumul de rangs partiels (Ikk ) chaque rang partiel correspondant à une desdites expressions combinatoires. 13. Procédé selon la revendication 12, caractérisé en ce que le calcul du rang d'une permutation (OP2) intervient, au codage de signaux numériques en 15 quantification vectorielle, pour indexer les permutations des composantes d'un vecteur directeur, ces permutations ayant été effectuées à une étape préalable (OP1) pour déterminer un vecteur-code le plus proche d'un vecteur d'entrée. 14. Procédé selon l'une des revendications 12 et 13, caractérisé en ce que, au 20 décodage de signaux numériques en quantification vectorielle, à partir d'une valeur donnée d'un rang de permutation (t) : - on calcule au moins une quantité représentative d'un rang de permutation (OP3) approchant ladite valeur donnée, en fonction d'au moins une composante présumée (dk) d'un vecteur-code à construire, 25 - on valide le choix de la composante présumée (dk) si ladite quantité représentative vérifie une condition de proximité avec la valeur donnée du rang. 203 15. Procédé selon la revendication 14, caractérisé en ce que ladite condition de proximité est vérifiée si la valeur donnée du rang (t) peut être encadrée par des cumuls de rangs partiels jusqu'au rang partiel associé à la composante présumée (dk), d'une part, et jusqu'au rang partiel associé à une composante correspondant à une incrémentation de la composante présumée (dk+1), d'autre part. 16. Procédé selon l'une des revendications 13 à 15, caractérisé en ce que le 10 codage/décodage est de type codage/décodage de source avec quantification vectorielle. 17. Procédé selon la revendication 12, caractérisé en ce que le codage/décodage est de type codage/décodage de canal, en modulation, et 15 en ce qu'il comporte : - avant transmission, la détermination d'un vecteur-code à partir du rang d'une permutation (OP3), et - en réception, le calcul du rang d'une permutation à partir d'un vecteur-code correspondant à un vecteur reçu (OP2). 18. Procédé selon l'une des revendications 12 à 17, prises en combinaison avec la revendication 3, dans lequel les codes à permutation ont une dimension maximale n, caractérisé en ce que lesdits nombres entiers choisis comportent au moins : 25 - les entiers compris entre 1 et la dimension maximale n, - la valeur factorielle de l'entier 0,4 - et, de préférence, les valeurs factorielles des entiers compris entre 1 et la dimension maximale n. 19. Procédé selon la revendication 18, caractérisé en ce que lesdits nombres 5 entiers choisis comportent en outre la valeur de 0. 20. Procédé selon l'une des revendications 12 à 19, caractérisé en ce que le code de permutation utilise une énumération de Schalkwijk, et en ce qu'un rang partiel Ikk associé à une troncature (yk,...yn_,) d'un vecteur-code 10 (yo,... yn_i) s'écrit : Iak k - (n-1-k)! ~k l l 1 wk , où fi( wk !) d-0 1=0 m - la notation fl représente un produit pour un indice entier i allant de 0 à m, 1=0 m - la notation E représente une somme pour un indice i allant de 0 à m, 1=0 - la notation 1! est la valeur factorielle de l'entier 1 avec l! =l x 2x 3x... x (l-1) x 1 , 15 pour 1>0 et avec 0!=1, - l'entier n est la dimension du code de permutation, correspondant au nombre total de composantes que comporte un vecteur-code, - l'entier k, compris entre 0 et n-1, est l'indice de la 1 ième composante Yk du vecteur-code, 20 - l'entier q est le nombre de composantes distinctes que comporte le vecteur-code, et - un terme wkd représente le nombre de composantes d'indices compris entre k et n-1 et qui ont une valeur égale à celle d'une même composante d'indice d.5 21. Procédé selon la revendication 20, prise en combinaison avec la revendication 9, caractérisé en ce que lesdits nombres entiers, identifiés et apparaissant dans un produit et/ou un quotient formant le rang partiel ( Ikk) en tant qu'expression combinatoire, sont : - les termes factoriels (n û l û k)! , pour tous les entiers k et compris entre 0 et n-1, - la valeur de chaque terme wk et/ou sa valeur factorielle, intervenant dans le q-1 produit Pk = fl(wk!) , chaque terme wk étant compris entre 0 et n, et =o dk -1 Io - les termes Sk = E wkd , chacun compris entre 1 et n-1, pour tous les entiers k d=o compris entre 0 et n-1. 22. Procédé selon l'une des revendications 20 et 21, prises en combinaison avec la revendication 13, caractérisé en ce que le calcul du rang intervient au 15 codage en déterminant le cumul des rangs partiels Ikk , chacun estimé en fonction d'une troncature (yk,...yä_I) du vecteur-code (yo,...yä_1), pour un indice k allant de n-1 à 0. 23. Procédé selon la revendication 22, caractérisé en ce que l'on mène un test dk-1 20 sur la valeur du terme Sk = E wk pour éviter l'estimation du rang partiel Ikk si d=o le terme Sk est nul. 24. Procédé selon l'une des revendications 20 et 21, prises en combinaison avec l'une des revendications 14 et 15, caractérisé en ce que le calcul d'une6 quantité représentative d'un décumul des rangs partiels Ik' intervient au décodage pour déterminer les composantes Yk du vecteur-code, d'indice k allant de 0 à n-1. 25. Procédé selon l'une des revendications 22 à 24, prises en combinaison avec l'une des revendications 10 et 11, caractérisé en ce que la somme des y-1 exposants dans la décomposition du terme Pk = fJ(wk!) est stockée ,=o temporairement en mémoire pour un indice k précédent pour être ajoutée ou soustraite aux exposants de la décomposition d'un terme wk pour un indice k lo courant. 26. Programme informatique destiné à être stocké en mémoire d'un dispositif de codage/décodage de signaux numériques, utilisant des codes à permutations, caractérisé en ce qu'il comporte des instructions pour la mise en 15 oeuvre des étapes du procédé selon l'une des revendications 1 à 25. 27. Dispositif de codage/décodage de signaux numériques, utilisant des codes à permutations, caractérisé en ce qu'il comporte : - une unité mémoire pour stocker les instructions d'un programme informatique 20 selon la revendication 26 et des représentations préenregistrées de décompositions de nombres entiers choisis, et - un module de calcul accédant à ladite unité mémoire pour la mise en oeuvre du procédé selon l'une des revendications 1 à 25.