Brevet FR2952253A1 — publié le 20110506

i PROCÉDÉ ET DISPOSITIF D'ÉVALUATION D'UNE TRANSMISSION NUMÉRIQUE SÉRIE CODÉE SUR UN CANAL NON STATIONNAIRE L'invention concerne un procédé et un dispositif informatique d'évaluation d'une transmission numérique série codée et modulée sur un canal bruité 5 à atténuation non stationnaire entre : un émetteur comprenant : ù un dispositif de codage adapté pour générer, à partir d'un flux de bits à transmettre, dit flux émis de bits, au moins un flux de mots codés, dit flux émis de mots codés, résultant du codage, selon au moins un procédé de codage 10 prédéterminé, dudit flux émis de bits, ù un dispositif de modulation adapté pour générer au moins un flux de symboles modulés, dit flux émis de symboles modulés, selon un schéma prédéterminé de modulation (notamment modulation en phase PSK et/ou modulation d'amplitude en quadrature QAM), sur au moins un signal de porteuse, chaque flux émis 15 de symboles modulés étant représentatif d'au moins une partie de chaque flux émis de mots codés, ù un dispositif d'émission, sur un canal bruité à atténuation non stationnaire, de chaque flux émis de symboles modulés, ù et un récepteur comprenant : 20 ù un dispositif de réception adapté pour recevoir, pour chaque flux émis de symboles modulés sur ledit canal, un flux de symboles modulés, dit flux reçu de symboles modulés, ù un dispositif de démodulation adapté pour générer au moins un flux de mots codés, dit flux reçu de mots codés, à partir de chaque flux reçu de 25 symboles modulés, ù un dispositif de décodage adapté pour générer un flux de bits, dit flux reçu de bits, par décodage, selon le procédé de codage mis en oeuvre par l'émetteur, de chaque flux reçu de mots codés. Le codage d'une telle transmission numérique permet de renforcer sa fiabilité. Dans diverses applications, le canal de transmission utilisé, en général de type sans fil (radiofréquences et/ou hyperfréquences...) présente une atténuation non stationnaire, c'est-à-dire qui varie sensiblement dans le temps au cours de la transmission de chaque mot codé. Ce phénomène est renforcé par la présence d'un entrelaceur canal. Tel est le cas par exemple avec des terminaux récepteurs mobiles (par exemple de type GPRS ou UMTS, Satellite DVB-SH) et/ou lorsque le codage est du type «turbocode», LDPC ou itératif En effet, les conditions de canal et ses performances peuvent être fluctuantes en fonction de la position du récepteur. En outre, lorsqu'un entrelaceur est prévu, la fonction d'un tel entrelaceur a pour effet de réduire les variations de l'atténuation vue par mot codé. Néanmoins, il est possible d'obtenir ou de générer par simulation des données numériques représentatives des variations au cours du temps de l'atténuation et du bruit du canal. Ces données numériques peuvent provenir des caractéristiques connues de l'entrelaceur, ou encore de caractéristiques statistiques connues à partir des propriétés physiques d'un canal. Les données du canal sont estimées à partir du signal reçu soit sur des symboles pilotes soit sur des symboles modulés. Les rapports signal/bruit successifs par mot codé sont reconstitués à partir de l'entrelaceur canal. Dans ce contexte général, un problème qui se pose est celui de l'évaluation prédictive des performances de la transmission, c'est-à-dire de la détermination d'un taux d'erreurs ER (taux d'erreurs de bits BER et/ou taux d'erreurs de mots PER) dans le flux reçu de bits en fonction de variations d'atténuation du canal au cours de la réception de chaque mot codé, et ce indépendamment de la réalisation réelle de la transmission et de ses différents éléments constitutifs (émetteur, récepteur, autres circuits, liaison physique...), notamment avant ou au cours de la conception de tout ou partie de ces éléments constitutifs, et donc en particulier sans réalisation du décodage. Une telle évaluation prédictive de perfonnances doit permettre en 30 particulier d'optimiser la conception desdits éléments constitutifs, et notamment de choisir des protocoles appropriés pour assurer une bonne qualité de transmission : requête de retransmission automatique (technique dite ARQ) éventuellement hybride (H-ARQ) ; redondance incrémentale (IR) ; combinaison de Chase ; adaptation des caractéristiques de l'émetteur et/ou de la liaison de transmission : choix du procédé de codage, puissance du signal, schéma de modulation... En particulier, une telle évaluation prédictive pourrait permettre d'éviter le surdimensionnement des récepteurs et de leurs antennes. Une première méthode connue de prédiction de performances visant à prendre en compte les variations d'atténuation d'un canal consiste à déterminer une valeur moyenne du rapport signal/bruit SNR sur chaque mot codé reçu. En particulier, le même traitement est appliqué pour deux mots codés différents présentant la même valeur moyenne du rapport signal/bruit mais des variations d'atténuation différentes. Cette méthode fournit en pratique des prédictions de performances trop optimistes. Dans une deuxième méthode connue, un taux d'erreurs de bits BER est déterminé pour chaque bit élémentaire de chaque mot codé reçu à partir d'une valeur moyenne du rapport signal/bruit moyen Eb/Nt pour ce bit élémentaire, elle-même calculée à partir d'une valeur moyenne du rapport signal bruit ES/Nt de chaque symbole modulé reçu divisée par le nombre de bits par symbole, la valeur de BER étant obtenue à partir d'une fonction connue (sous forme d'une courbe ou d'une table) donnant la valeur de BER en fonction du rapport signal/bruit sur un canal à bruit blanc gaussien additif Une moyenne du taux d'erreurs de bits est ensuite calculée pour le mot codé reçu et est utilisée en tant que critère de performance. Cette méthode est plus précise (c'est-à-dire fournit des résultats plus conformes à la réalité) que la précédente mais fournit des prédictions qui s'avèrent en pratique être trop pessimistes. EP 1564924 décrit un procédé dans lequel chaque signal codé reçu est divisé en une pluralité de segments de longueur fixe prédéterminée, et la valeur moyenne sur la trame du rapport signal/bruit par symbole (E/N)F est obtenue par calcul pour chaque segment d'un paramètre Cm à l'aide d'une fonction f convexe 4 numérique obtenue empiriquement en fonction du schéma de modulation (BPSK ou QPSK) appliqué sur la valeur moyenne (Es/N)M de chaque segment, puis calcul d'une valeur moyenne C de ce paramètre sur la trame, puis calcul inverse avec la fonction inverse f 1. Ensuite, la valeur moyenne sur la trame du rapport signal/bruit moyen (Eb/N)F pour chaque bit élémentaire est obtenue à partir d'une table, et le taux d'erreurs par trame FER est obtenu en fonction d'une courbe de référence donnée sur un canal à bruit gaussien additif. Cette méthode repose sur une approximation numérique de la fonction f. Elle est limitée aux schémas de modulation pour lesquels cette fonction numérique est donnée dans ce document, à savoir BPSK et QPSK. Elle n'est donc pas applicable à d'autres schémas de modulation. Par ailleurs, elle ne prend pas en considération les variations d'atténuation du canal au niveau des symboles, mais avec une granularité au niveau de segments de longueur déterminée. Sa précision dépend en réalité de la bonne adaptation de la taille, fixe, de chaque segment à la dynamique des variations d'atténuation, adaptation qu'il n'est pas possible d'ajuster au préalable en fonction de chaque application. Par ailleurs, cette évaluation prédictive doit pouvoir être effectuée sans modélisation et simulation de la liaison physique elle-même, une telle modélisation entraînant en effet des calculs extrêmement lourds et complexes. À ce titre, il est à noter qu'il est particulièrement important de pouvoir réaliser une telle évaluation prédictive qui simultanément présente une grande fiabilité, prenne correctement en considération des fluctuations importantes et réalistes du canal (par exemple similaires à celles rencontrées en pratique dans le cas de tenninaux mobiles tels que des téléphones mobiles de type GSM ou GPRS, ou satellite DVB-SH), et soit suffisamment rapide pour pouvoir être exploitée à l'échelle industrielle, notamment dans la conception des éléments constitutifs de la transmission. La rapidité de l'évaluation prédictive représente un enjeu déterminant dans le cadre de la fabrication industrielle moderne de circuits dans laquelle les durées des différentes phases sont déterminées avec précision, doivent impérativement être respectées, sont cruciales et tendent en permanence à se réduire. Cette évaluation prédictive doit pouvoir aussi être effectuée sans fabrication réelle des éléments constitutifs de la liaison, notamment du récepteur. En effet, en alternative aux méthodes connues de simulation mentionnées ci-dessus, la 5 seule méthode réellement efficace envisagée jusqu'à maintenant consiste à fabriquer réellement tout ou partie des éléments du récepteur, notamment ceux qui sont les plus complexes et dont la simulation infonmatique sont les plus coûteux, notamment en temps de traitement informatique. Une telle méthode impose cependant une fabrication longue et coûteuse et différents éléments non simulés. Dans ce contexte, l'invention vise à pallier l'ensemble de ces inconvénients en proposant un procédé d'évaluation de performances susceptible à la fois de prendre en compte les véritables variations d'atténuation du canal, et avec une bonne précision des résultats, une grande fiabilité et des traitements informatiques légers et rapides, et à faible coût -notamment sans nécessiter aucune fabrication réelle des éléments constitutifs de la transmission, notamment du récepteur. En particulier, l'invention vise à proposer un tel procédé permettant d'obtenir un résultat en une durée inférieure d'une part à celle nécessaire à l'obtention d'une évaluation à partir d'une transmission numérique réellement réalisée et mise en oeuvre, d'autre part à celle nécessaire à la réalisation d'une modélisation et d'une simulation complète de la liaison physique de la transmission. Plus particulièrement, l'invention vise à proposer un tel procédé qui soit compatible avec les contraintes actuelles d'une exploitation à l'échelle industrielle, dans la conception des éléments constitutifs de la transmission, notamment du récepteur. L'invention vise également à proposer un tel procédé qui 25 permette d'obtenir des statistiques fiables de performances de ladite transmission au cours du temps. L'invention vise également à proposer un tel procédé qui permet d'optimiser la conception des éléments constitutifs de la transmission, notamment du récepteur, en évitant tout surdiinensionneinent. 30 L'invention vise également à proposer un tel procédé qui puisse 6 être applicable à tout schéma de modulation, y compris autre que BPSK ou QPSK, notamment 2n-PSK n>3 ou QAM. Pour ce faire l'invention concerne un procédé d'évaluation des performances d'une transmission numérique série codée et modulée sur un canal bruité 5 à atténuation non stationnaire entre : û un émetteur comprenant : û un dispositif de codage adapté pour générer, à partir d'un flux de bits à transmettre, dit flux émis de bits, au moins un flux de mots codés, dit flux émis de mots codés, résultant du codage, selon au moins un procédé de codage 10 prédéterminé, dudit flux émis de bits, un dispositif de modulation adapté pour générer au moins un flux de symboles modulés, dit flux émis de symboles modulés, selon un schéma prédéterminé de modulation, sur au moins un signal de porteuse, chaque flux émis de symboles modulés étant représentatif d'au moins une partie de chaque flux émis de 15 mots codés, û un dispositif d'émission, sur un canal bruité à atténuation non stationnaire, de chaque flux émis de symboles modulés, et un récepteur comprenant : û un dispositif de réception adapté pour recevoir des flux de 20 symboles modulés, dits flux reçus de symboles modulés, correspondant à des flux émis de symboles modulés sur ledit canal, û un dispositif de démodulation adapté pour générer au moins un flux de mots codés, dit flux reçu de mots codés, à partir de chaque flux reçu de symboles modulés, 25 û au moins un dispositif de décodage adapté pour générer un flux de bits, dit flux reçu de bits, par décodage de chaque flux reçu de mots codés, selon un procédé de décodage correspondant à un procédé de codage mis en oeuvre par l'émetteur, procédé dans lequel au moins une valeur représentative d'un taux d'erreurs ER du flux 30 reçu de bits, est élaborée sans réalisation du décodage, à partir de données numériques mémorisées représentatives des variations au cours du temps de l'atténuation et du bruit du canal, caractérisé en ce que au moins une valeur, dite qualité de transmission Qk, de formule ck2.Es/No, où Ck représente chaque valeur d'atténuation du canal au cours du temps, k étant un indice temporel, ES représente une énergie moyenne par symbole émis et No représente une densité spectrale d'un bruit blanc gaussien sur le canal, et les variations au cours du temps de ladite qualité de transmission Qk, sont élaborées pour chaque symbole reçu du flux reçu de symboles modulés, et en ce que: ù dans une première étape, une valeur d'information mutuelle Ik est élaborée pour chaque valeur de ladite qualité de transmission Qk, selon une fonction prédéterminée de ladite qualité de transmission Qk, ù dans une deuxième étape, une moyenne d'information mutuelle est élaborée pour chaque mot codé du flux reçu de mots codés, en réalisant une moyenne des différentes valeurs d'information mutuelle Ik élaborées dans la première étape pour les différentes valeurs prises par ladite qualité de transmission Qk sur le dit mot codé, ù dans une troisième étape, au moins une valeur du taux d'erreurs ER du flux reçu de bits est élaborée pour chaque mot codé du flux reçu de mots codés, à partir de chaque valeur de la moyenne d'information mutuelle élaborée dans la deuxième étape, et par utilisation de données mémorisées représentatives de variations d'un taux d'erreurs équivalent selon au moins une fonction, dite fonction standard, du rapport signal/bruit, chaque fonction standard étant prédéterminée pour les dispositifs de codage et de décodage sur un canal à bruit blanc gaussien additif. Dans un premier mode de réalisation possible et avantageux de l'invention, chaque valeur d'information mutuelle Ik est élaborée selon la fonction définie par la formule (I) suivante : ~k (xk , .yk) = ù log2 (Tieck ES) ù r f f (u, v) loge (f (u, v))dudv No J J M ù(vùck ES ximagSn,))2 ù(uùck Es xrea/S,,,))2 E lexp( N° E ) exp( N° E ) MTC C2 s ) m 0 c2 s c2 s (k N° k N° k N° Sm JEs M étant le cardinal de l'alphabet A = {s0, S1 , .. . , SMû1} des symboles modulés. Avantageusement et selon l'invention, cette formule analytique peut être discrétisée pour son évaluation par traitement numérique, par exemple selon la fonction définie par la formule (II) suivante : I k (xk, Yk) = ùloge (neck N E ù) 0 seuil seuil - AuAv 1 1 f (qdu, rAv) loge (f (qAu, r/v)) q=ùseuil r= seuil ck E s n- ) + Maxo k No k No k c2 No ù(vùc2 Es ximag(Sm))2 ù(uùck Es xread(S,,))2 2 1 No ) ) m S 10 M étant le cardinal de l'alphabet A = {S° ,S1 • • • SM_1 } des symboles modulés, Sm étant l'ensemble des symboles normalisés m de p bits ° p -1 dont le bit numéroté m vaut 0, m 0 Avantageusement et selon l'invention, cette formule analytique peut être discrétisée pour son évaluation par traitement numérique, par exemple selon Bi'° m 0 seuil seuil Ik (bh , Yk) = AuAv f(gAu, rAv) log2 (.Î (qAu, rAv)) ~q_ seuil r=ùseuil seuil seuil + E E gi (qAu, rAv) log2 (g; (qAu, rAv)) q=ùseuil r=ùseuil J seuil = - ck N 1n( a 0 Au=Av=2 seuil fi Ck NS ) + Max 0 par la fonction inverse Ik 1, - puis chaque valeur du taux d' erreurs ER est obtenue à partir de ladite valeur de rapport signal/bruit équivalent SNReq et desdites données mémorisées représentatives de variations d'un taux d'erreurs équivalent selon une fonction standard 20 prédéterminée pour les dispositifs de codage et de décodage sur un canal à bruit blanc gaussien additif. D'autres modes de réalisation sont possibles. Un procédé selon l'invention est un procédé de simulation assisté par ordinateur, c'est-à-dire est susceptible d'être mis en oeuvre par un dispositif informatique. Un procédé d'évaluation selon l'invention est mis en oeuvre par traitement numérique de données, c'est-à-dire est un procédé infonnatique, mis en oeuvre par un dispositif informatique. À l'instar de tout procédé de simulation informatique, le procédé selon l'invention permet de fournir des résultats directement utilisables pour la conception et la fabrication des éléments constitutifs de la transmission -notamment du récepteur. L'invention s'étend également à un dispositif informatique pour la mise en oeuvre d'un procédé selon l'invention, c'est-à-dire à un dispositif informatique comprenant des moyens informatiques de traitement numérique de données, caractérisé en ce qu'il est adapté et programmé pour la mise en oeuvre d'un procédé selon l'invention. L'invention s'étend également à un produit logiciel apte à être chargé en mémoire vive d'un dispositif informatique, caractérisé en ce qu'il est adapté pour la mise en oeuvre d'un procédé selon l'invention. L'invention s'étend encore à un support d'enregistrement adapté pour pouvoir être lu par un lecteur relié à un système informatique, caractérisé en ce qu'il comprend un programme enregistré adapté pour pouvoir être chargé en mémoire vive d'un système informatique et le programmer pour mettre en oeuvre un procédé selon l'invention. Les inventeurs ont constaté que l'invention permet en pratique d'obtenir une prédiction de performances rapide et précise, et en particulier considérablement plus précise que dans toutes les méthodes antérieures connues, tout en restant plus rapide que toutes les méthodes de simulation antérieures et que la prédiction de performances sur un canal réel. Ce résultat est surprenant, notamment car il n'y a pas a priori de relation directe entre l'information mutuelle sur les bits codés et les performances du décodeur. Or, le calcul de cette information mutuelle et de sa moyenne sur chaque mot codé permet en pratique de prendre en compte les variations temporelles d'atténuation avec une grande précision. Il apparaît que deux mots codés ayant la même information mutuelle moyenne présentent la même performance en sortie du décodeur. En outre, les inventeurs ont déterminé que la valeur de 11 12 l'information mutuelle peut même être élaborée à partir d'une formule analytique, dont la mise en oeuvre est donc simple et précise par calcul numérique, y compris pour le calcul de la fonction inverse. L'invention concerne également un procédé, un dispositif 5 infonnatique, un produit logiciel et un support d'enregistrement caractérisés en combinaison par tout ou partie des caractéristiques mentionnées ci-dessus ou ci-après. D'autres buts, caractéristiques, et avantages de l'invention apparaîtront à la lecture de la description suivante qui se réfère aux figures annexées représentant, à titre d'exemples non limitatifs, des modes de réalisation de l'invention et 10 dans lesquelles : ù la figure 1 est un schéma fonctionnel représentant un dispositif de transmission pouvant faire l'objet d'un procédé d'évaluation selon l'invention, ù la figure 2 est un organigramme schématique représentant 15 un mode de réalisation des principales étapes d'un procédé selon l'invention, ù la figure 3 est un diagramme schématique d'exemples de courbes de référence pouvant servir, selon différents schémas de modulation, d'une fonction prédéterminée pour la détermination d'une valeur d'information mutuelle sur un canal à bruit blanc gaussien additif en fonction du rapport signal à bruit de ce canal, 20 ù la figure 4 est un diagramme schématique représentant l'utilisation d'une courbe de référence similaire à la figure 3 pour l'exécution d'une première étape d'un procédé selon l'invention, ù la figure 5 est un diagramme schématique représentant l'utilisation d'une courbe de référence similaire à la figure 3 pour l'exécution d'une 25 première partie d'une troisième étape d'un procédé selon l'invention, la figure 6 est un diagramme schématique représentant l'utilisation d'une courbe représentative d'une fonction standard pour l'exécution d'une deuxième partie d'une troisième étape d'un procédé selon l'invention, - la figure 7 est un schéma bloc représentant les différentes 30 fonctions mises en oeuvre dans un procédé selon l'invention, 13 û la figure 8 est un diagramme illustrant un exemple de données numériques mémorisées représentatives de variations au cours du temps de l'atténuation du canal pouvant être utilisées dans un procédé selon l'invention, û la figure 9 est un diagramme illustrant un exemple de taux 5 d'erreur par mot codé PER pouvant être produit par un procédé selon l'invention, û la figure 10 est un diagramme illustrant un autre exemple de données pouvant être délivrées par un procédé selon l'invention, û la figure 11 illustre des courbes représentatives de la fonction de répartition d'atténuations par ombrage d'un canal bande S en 10 environnement suburbain pour différentes distances parcourues. La figure 1 représente de façon générale un dispositif de transmission numérique série codée et modulée sur un canal bruité à atténuation non stationnaire. Ce dispositif comprend un émetteur 11, un récepteur 12 et une liaison physique 13 sans fil formant le canal de transmission. La liaison physique 13 peut être 15 par exemple une liaison radio fréquence, telles que celles par exemple reliant des terminaux mobiles tels que des téléphones cellulaires, des assistants numériques personnels, des ordinateurs portables, des cartes sans fil, des véhicules terrestres, des navires, des aéronefs, des satellites, des sondes spatiales ou autres systèmes spatiaux... à une station de base, elle-même fixe (terrestre) ou mobile (véhicule, satellite,...) 20 accédant à un réseau de transmission de données tel que le réseau Internet ou tout autre réseau privé. En général, la transmission peut être bidirectionnelle, c'est-à-dire que chaque terminal mobile est tantôt émetteur, tantôt récepteur. L'émetteur 11 comprend un premier module 14 délivrant des données sous forme d'un flux de bits (signal en bande de base) à transmettre, dit flux 25 émis de bits. Ce flux émis de bits est fourni à un module 15 de codage qui exécute un procédé de codage prédéterminé pour former, à partir des bits, un flux de mots codés, dit flux émis de mots codés. Un tel procédé de codage permet en particulier d'augmenter la fiabilité des données transmises en augmentant les redondances tout en assurant la correction d'erreurs, c'est-à-dire la restitution des données initiales malgré 30 les perturbations que peut subir le canal 13 de transmission. 14 Le module 15 de codage délivre des mots codés, qui sont ensuite modulés, selon un schéma de modulation prédéterminé, par un circuit modulateur 18 qui fournit un flux de symboles modulés, dit flux émis de symboles modulés, à un circuit entrelaceur 17, ce dernier fournissant des symboles modulés et entrelacés à un circuit émetteur 19 radiofréquence apte à émettre des signaux correspondant sur la liaison physique 13. Le récepteur 12 comprend un circuit de réception 20 radiofréquence apte à recevoir les signaux transmis via la liaison physique 13, et à délivrer un flux de symboles modulés et entrelacés, dit flux reçu de symboles modulés et entrelacés, à un circuit désentrelaceur 22 qui effectue le traitement inverse de l'entrelaceur 17 de l'émetteur 11, c'est-à-dire permet la reconstitution progressive d'un flux de symboles modulés désentrelacés, à partir du flux de données numériques issues du circuit de réception 20. Le circuit désentrelaceur 22 fournit les symboles modulés désentrelacés à un circuit démodulateur 21 apte à appliquer une démodulation, sur les symboles reçus, selon le schéma et la cartographie de modulation utilisés à l'émission. Le circuit démodulateur 21 permet la reconstitution progressive d'un flux de mots codés, dit flux reçu de mots codés, à partir du flux de symboles modulés désentrelacés. Ces mots codés reçus sont ensuite fournis à une mémoire tampon 23 puis traités progressivement par un module 25 de décodage, comprenant un ou plusieurs décodeurs -notamment deux décodeurs-, et permettant de délivrer un flux 27 de bits reçus inclus dans le signal véhiculé par la liaison physique 13 et correspondant au flux de bits émis par le circuit générateur 14. L'invention concerne un procédé d'évaluation, par simulation informatique, des performances d'une telle transmission numérique série codée. L'invention s'applique à tout procédé de codage, et indépendamment de la nature exacte du procédé de codage utilisé. Il peut s'agir en particulier d'un procédé de codage choisi panni les procédés du type dit LDPC (code à matrice de parité creuse), les procédés du type à turbocode et les autres procédés de codage à décodage itératif Dans la plupart des procédés de codage modernes qui permettent d'obtenir des performances proches de la limite de Shannon, le module 15 de codage comprend une pluralité de 15 codeurs -notamment deux codeurs-. L'invention s'applique cependant aussi bien à un module 15 de codage comprenant un seul codeur. Les symboles modulés Xk émis en sortie du bloc modulateur 18 (qui réalise la cartographie de la modulation, souvent désignée par le terme anglais "mapping") appartiennent à une constellation de symboles complexes de cardinal M qui est une puissance de 2 (M=2p). Ils ne peuvent donc prendre que M valeurs complexes. Ces symboles sont ensuite transmis sur le canal 13. L'effet du canal sur ces symboles Xk est double : û le canal 13 atténue chaque symbole Xk et on appelle Ck l'atténuation subie par 10 un symbole Xk, û le canal 13 affecte le symbole atténué xk d'un bruit blanc Gaussien complexe additif de puissance N0/2 sur les voies imaginaires et réelles. N0 représente ainsi une densité spectrale d'un bruit blanc gaussien sur le canal. Ce bruit blanc gaussien inclut le bruit thermique et éventuellement des interférences. 15 On appelle alors y, le symbole reçu en entrée du démodulateur 21 correspondant donc au symbole Xk émis atténué et bruité. Les valeurs des modules des atténuations ICkl prises au cours du temps sont fournies par un module 24 de simulation du canal (figure 7), ce dernier déterminant les valeurs des atténuations à partir de données représentatives de la vitesse de déplacement du récepteur et/ou de 20 caractéristiques connues prédéterminées concernant le canal lui-même et la liaison physique. Ces valeurs sont par exemple déterminées à partir de véritables séries de valeurs d'atténuation ou de propagation en fonction de la distance parcourue, mesurées lors d'essais réellement réalisés, puis lues à une vitesse adaptée à la vitesse d'un terminal mobile considéré dans l'exécution d'un procédé selon l'invention. En variante, 25 ces valeurs peuvent être déterminées à partir d'un modèle statistique du canal (par exemple le modèle à trois états de Perez-Fontan pour un canal LMS satellite) pour lequel un tirage des échantillons successifs Ck est effectué. La figure 8 représente un exemple de valeurs des modules des atténuations susceptibles d'être fournies par le module 24 de simulation du canal. La figure 8 n'étant qu'illustrative, les valeurs ne sont représentées que sur une faible période de temps, de 60 secondes. En pratique, la période utilisée pour l'évaluation est en générale beaucoup plus importante, typiquement de plusieurs heures. Si ES représente une énergie moyenne par symbole émis, Es/No représente le rapport signal/bruit à l'émission. La qualité Qk de transmission peut être ainsi représentée pour chaque symbole reçu xk par la formule : Qk = ck2.Es/No Si Eb désigne l'énergie moyenne par bit codé émis et Es 10 l'énergie moyenne par symbole émis, ces deux grandeurs sont liées par la relation : Es =Eb xlog2(M) La valeur de ES est calculée à partir des différents symboles de la constellation A={SoSMù1} par la relation E = 1 '1 S 2 s M mI 15 On réalise l'égalisation du canal en définissant des nouveaux = 1 symboles normalisés Sm E S17 s Le procédé selon l'invention permet d'évaluer par simulation un taux d'erreurs, pour chaque symbole reçu, en prenant en compte les variations d'atténuation du canal sur les symboles reçus, c'est-à-dire sur les différentes parties des 20 mots codés reçus correspondants. Les différents modules réalisant les fonctions mises en oeuvre dans un procédé selon l'invention sont représentés schématiquement figure 7. Le module 24 de simulation du canal fournit un flux simulé de symboles susceptibles d'être reçus par le récepteur 20. Un module 28 de simulation du désentrelaceur permet de simuler m=0 17 la fonction du désentrelaceur 22 et du démodulateur 21, et de regrouper par mots codés les différents symboles fournis par le module 24 de simulation de canal. Le flux de mots codés désentrelacés et démodulés est mémorisé dans une mémoire tampon 23 (figure 1). Pour chaque mot codé, un module 26 d'évaluation (figure 7) d'un procédé selon l'invention calcule une valeur d'un taux d'erreur par bit BER et/ou une valeur d'un taux d'erreur par mot codé PER. La figure 2 représente ainsi un organigramme d'un exemple d'étapes mises en oeuvre dans un module 26 d'évaluation selon l'invention. Lors de l'étape 33, les différentes valeurs de l'atténuation ck et/ou du rapport signal bruit Es/No et/ou de la qualité de transmission Qk sont calculées, et ce pour chaque symbole reçu, et pour chaque mot codé reçu correspondant mémorisé dans la mémoire 23, du flux reçu de mots codés. À partir de ces différentes valeurs Qk, une valeur d'information mutuelle Ik est élaborée, lors de l'étape subséquente 34, selon 15 une fonction prédéterminée. Deux variantes de l'invention sont possibles en ce qui concerne le choix de la fonction prédéterminée pour le calcul de ladite valeur d'information mutuelle Ik. Si l'on considère un symbole reçu yk, ce symbole correspond au 20 symbole émis Xk atténué par Ck et bruité par un bruit blanc gaussien additif. En raison de la modulation choisie de cardinal M=2p, ce symbole reçu Yk contient l'information relative à p bits émis bk1 bk2 bkp Une première variante de réalisation de l'invention consiste à considérer l'information mutuelle I((bki,bk2,. ,b p), yk) égale à l'information mutuelle 25 Ik entre le symbole émis Xk et le symbole reçu yk. Cette variante a l'avantage de proposer une formule théorique unique de calcul d'une courbe de référence de l'information mutuelle qui dépend seulement de la nature de la modulation considérée. Cette approche ne tient pas compte de la mise en oeuvre réelle (« mapping ») de la 10 510 modulation ni du fait que les p bits constitutifs du symbole Xk ne sont pas tous protégés de la même façon contre le bruit du canal. Dans cette première variante, chaque valeur d'information mutuelle Ik est déterminée selon la fonction définie par la formule (I) susmentionnée. L'expression de f(u, v) est sous forme analytique et peut donc être évaluée pour toutes les valeurs de u et v. Pour l'intégrale, il suffit de discrétiser cette intégrale par toute méthode numérique. Par exemple en utilisant une méthode des rectangles, on obtient chaque valeur d'information mutuelle Ik selon la fonction définie par la formule (II) suivante : _ Es ) l k (xk , Yk) - -log2 (lieck N 0 seuil seuil - AuOv E .f (qOu, rAv) log2 (.f (qou, rAv)) q=ùseuil r=ùseuil ù ck N ln(a o et Au=0v=2 seuil ck N ) + Maxo Ainsi, dans cette deuxième variante, chaque valeur d'information mutuelle Ik est déterminée selon la fonction Ik(b'k,yk) de l'information mutuelle calculée entre le j ème bit (0 10 Cette formule (III) générale est valable dans le cas où géométriquement ces deux ensembles de symboles normalisés ont la même distribution géométrique dans le plan complexe, c'est-à-dire lorsque ces deux ensembles sont égaux à une rotation ou translation du plan près. Dans le cas contraire, une formule analytique peut aussi être écrite mais elle est plus complexe. L'hypothèse susmentionnée englobe 15 les modulations PSK en distribution géométrique dite « lapping Gray ». Les expressions de f(u, v) et g(u, v) sont sous forme analytique et peuvent donc être évaluées pour toutes les valeurs de u et v. Pour l'intégrale, il suffit de discrétiser cette intégrale par toute méthode numérique. Par exemple en utilisant une méthode des rectangles, on obtient chaque valeur d'information mutuelle Ik selon la 20 fonction Ik(bjk, yk) de l'information mutuelle calculée entre le j ème bit (0 _ seuil seuil I k (bk , Yk) = Au/v E f (qAu, rAv) log 2 (f (qAu, rLv)) q=ùseuil r= ù seuil seuil seuil + E 1 g.(qAu, rOv) log 2(g .(qLu, rAv)) q= ù seuil r = ù seuil Avec, si l'on veut évaluer les différentes gaussiennes de f(u,v) et 25 gj(u,v) jusqu'à au moins leur probabilité a sur un nombre de points [3 dans le maillage, les valeurs suivantes de seuil et de (°"'°") : û ck ES ln(a ck ES z ) + Maxo d'information mutuelle est élaborée par le module 26 évaluation en réalisant une moyenne des différentes valeurs d'information mutuelle Ik déterminées dans l'étape 33 pour les différentes valeurs prises par ladite qualité de transmission Qk sur le dit mot codé. Lors de l'étape subséquente 36, le module 26 d'évaluation réutilise la même courbe de référence (c'est-à-dire la même table de valeurs numériques) pour déterminer une valeur de rapport signal/bruit équivalent SNReq sur 20 seuil = et Au=4v=2 21 le mot codé à partir de la dite moyenne d'information mutuelle . Autrement dit, le module 28 de prédiction de performances utilise la fonction inverse Ik 1. Lors de l'étape subséquente 37, le module 26 évaluation calcule chaque valeur du taux d'erreurs ER à partir de ladite valeur de rapport signal/bruit équivalent SNReq et de données mémorisées représentatives de variations d'un taux d'erreurs équivalent selon une fonction standard prédéterminée pour les modules de codage et de décodage utilisés sur un canal à bruit blanc gaussien additif. En effet, pour un module 15 de codage et un module 25 de décodage prédéterminés, il existe une fonction standard, obtenue de façon connue en soi par simulation sur le canal gaussien stationnaire, exprimant le taux d'erreurs par mot codé PER ou le taux d'erreurs par bit BER en fonction du rapport signal/bruit. La figure 6 représente un exemple de courbes représentatives de telles fonctions standard, les différentes courbes étant obtenues pour les mêmes modules de codage de décodage et variant les unes des autres en fonction du nombre d'itérations utilisées pour le décodage. Pour chaque valeur du rapport signal/bruit équivalent SNReq, le module 26 d'évaluation calcule, à partir de ces courbes, c'est-à-dire des tables de valeurs numériques enregistrées correspondantes, un ensemble de taux d'erreurs ER;, c'est-à-dire une valeur de taux d'erreurs pour chaque nombre d'itérations pouvant être utilisé au décodage. Les valeurs de taux d'erreurs ER, sont décroissantes avec le nombre d'itérations. Cet ensemble de valeurs ER, constitue donc un vecteur de taux d'erreurs (ER), déterminé pour chaque mot codé reçu. Dans l'exemple de la figure 6, on a représenté quatre courbes de fonctions standard : une première courbe CS 1 correspondant par exemple à une seule itération, une deuxième courbe CS2 correspondant par exemple à quatre itérations, une troisième courbe CS3 correspondant par exemple à huit itérations, et une quatrième courbe CS4 correspondant par exemple à seize itérations. On obtient ainsi, à partir de la valeur du rapport signal/bruit équivalent SNReq, un vecteur de taux d'erreurs ayant quatre composantes : (ER) = (ER1, ER2, ER3, ER4). 22 Il est à noter que les deux étapes 36, 37 peuvent être rassemblées en une même et seule étape 38, si les fonctions standard sont combinées avec la fonction inverse Ik' en une seule fonction fournissant directement, pour chaque nombre d'itérations, des courbes de variation du taux d'erreur ER en fonction de la moyenne d'information mutuelle . Le module 26 d'évaluation fournit les variations au cours du temps, pour les échantillons de canal fournis par le module 24 de simulation du canal et sur une même échelle de temps, les variations au cours du temps des valeurs de taux d'erreur ER (BER et/ou PER). La figure 9 représente un exemple d'un échantillon (sur une période de quatre secondes uniquement) de la courbe de variations en fonction du temps de la probabilité de taux d'erreur par mot codé PER pouvant être fournie par le module 26 d'évaluation à partir des valeurs d'atténuation (dont la figure 8 donne un exemple de représentation en échelle logarithmique et en fonction du temps pour un canal satellite LMS (satellite-sol) en milieu moyennent boisé avec une vitesse de 50 km/h), dans le cas de la simulation d'un turbo-décodeur selon la nonne DVB-SH pour une modulation QPSK en mode OFDM, avec un entrelaceur convolutionnel qui étale sur 100 ms, et un turbo code de taux de codage 1/2. L'ensemble des trois modules 24 de simulation du canal, 28 de simulation du désentrelaceur et 26 d'évaluation peuvent être incorporés au sein d'un même simulateur. Ce simulateur peut être constitué par exemple de deux blocs réalisés selon la plate-forme logicielle Juzzle et de librairies. Les données en entrée du simulateur sont : un fichier contenant les valeurs représentatives de la courbe représentée figure 4 fournissant les variations de l'information mutuelle en fonction du rapport signal/bruit, et un fichier contenant les valeurs représentatives des courbes représentées figure 6 exprimant les fonctions standard représentant les performances du turbo code par les variations du taux d'erreur en fonction du rapport signal/bruit sur un canal gaussien et selon le taux de codage considéré. Le module 24 de simulation du canal peut être formé d'un bloc de génération du canal, d'une librairie de modélisation du turbo code et d'une librairie de propagation permettant de générer des échantillons selon un ou plusieurs modèles de canaux, par exemple selon les modèles de Perez Fontan (cf Pérez-Fontân, F., Vâzquez-Castro, M., Buonomo, S., Poiares-Baptista, J.P., Arbesser-Rastburg, B., "S-Band LMS propagation channel behaviour for different environments, degrees of shadowing and elevation angles", IEEE Transactions on Broadcasting, Vol. 44, No.1, March 1998), Rice et Rayleigh, ou autres. Le module 28 de simulation du désentrelaceur est formé d'une librairie. Cette librairie est réalisée de façon à simuler un désentrelaceur en fonction du type d'entrelaceur considéré, qui peut être quelconque dans le cadre de la présente invention. Il peut s'agir par exemple d'un désentrelaceur de type bloc, notamment de type lignes/colonnes (remplissage une matrice en lignes et lecture en colonnes), ou encore d'un désentrelaceur de type convolutionnel ou autre. Par exemple, pour le cas du standard DVB-SH, l'entrelaceur est 15 de type convolutionnel à 48 branches avec une unité d'entrelacement (UI) de 126 bits et à chacune de ces 48 branches est associé un retard f(i) pour 0 Cela signifie que les bits à l'entrée de l'entrelaceur sont regroupés par paquets de 126 bits et que ces paquets sont distribués successivement sur ces 48 branches (branche 0, branche 1,.. . branche 47 puis de nouveau branche 0, 20 branche 1..). Chaque unité d'entrelacement est alors retardée du retard f(i) correspondant à la branche sur laquelle il se trouve. Les unités d'entrelacement en sortie sont récupérées alors branche par branche pour reformer un flux unique. Ulanp Ainsi si on désigne par n la n-ième unité d'entrelacement Ul out en entrée de l'entrelaceur et k la k-ième unité d'entrelacement en sortie de 25 l'entrelaceur, la n-ièine unité d'entrelacement en entrée est envoyée en sortie à la kièine position par l'entrelaceur selon la relation suivante (en numérotant à partir de 0 les indices) : Z%Inp Ülkut k=n48 +48*(f(n48)+n48) avec n=n +48*n48 Où -48 est la notation pour la division euclidienne de l'entier n par 48. L'entrelaceur de la norme DVS-SH est défini entièrement par 5 paramètres qui permettent de calculer les 48 valeurs de f(i) pour 0 Le désentrelaceur a le fonctionnement identique à l'entrelaceur avec une fonction f(i) inverse. Le simulateur selon l'invention fournit en sortie un fichier dont la courbe de la figure 9 constitue une représentation partielle (sur une durée uniquement de quatre secondes, alors que la période d'évaluation est bien sûr beaucoup plus importante, typiquement de plusieurs heures). Dans l'exemple d'un terminal récepteur mobile, la distance de parcours à considérer pour être représentative de la statistique du canal dépend du canal considéré. Par exemple, dans le cas d'un canal LMS bande S en environnement suburbain, une distance de l'ordre de 10 km est recommandée pour une bonne représentativité des évaluations de performances. La figure 1l illustre la fonction de répartition d'atténuations par ombrage d'un tel canal bande S en environnement suburbain (modèles de Perez- Fontan) pour plusieurs distances de canal parcourues (100m, 200m, 1 km, 5 km, 10 km, 20 km, 40 km). On constate sur cette figure que lorsqu'on cherche à assurer le service sur 95% de la zone suburbaine, l'atténuation la pire rencontrée n'est représentative de la statistique d'ensemble qu'à partir d'une distance de l'ordre de 10 km. En variante, un simulateur selon l'invention fournit un fichier représentatif de la disponibilité de la transmission, c'est-à-dire la prédiction de la probabilité au cours du temps qu'une fenêtre de Tl secondes présente au moins T2 secondes sans erreur (ER inférieur ou égal à une valeur seuil ERmax), T2 étant comprise entre zéro et T1 (par exemple, T1=20s). La figure 10 représente un exemple d'une telle sortie en fonction du rapport signal/bruit C/N exprimé en décibels. Chaque point de chaque courbe de la figure 10 correspond à l'exécution du procédé selon l'invention et à l'obtention des variations des valeurs de taux d'erreur ER au cours du temps telles que fournies par le module 26 évaluation. La courbe Cl correspond à T2=T1 ; la courbe C2 correspond à T2=T1-ls=19s ; la courbe C3 correspond à T2=T1-2s=18s ; chaque courbe Ci correspond à T2=T 1-(i-1) s. La courbe C2, qui est souvent considérée dans un tel contexte, correspond à la probabilité ESR5 d'avoir au moins 19 secondes justes pour chaque fenêtre de 20 secondes, soit une proportion de secondes fausses par fenêtre inférieure à 5%. Les points de ces courbes correspondant à une valeur donnée de C/N donnée sont obtenus de la façon suivante : - étape 1 : calcul pour chaque mot de code du PER estimé pour ce 15 mot de code grâce au simulateur simplifié. - étape 2 : calcul pour le ième bloc d'une seconde de la probabilité Pi que cette seconde soit sans erreur 4 J9 = n(j. - F M P. 2) j -a Où RERu est le Taux d'erreur paquet instantané fourni du j-ème 20 mot code constituant, N est le nombre de mots de code par seconde. - étape 3 : calcul sur la ième fenêtre de 20s de la loi Yi de distribution du nombre de secondes justes selon la formule : 20 -1 25 Yi est la convolution des 20 lois de Bernoulli : Al, est une variable aléatoire de Bernoulli qui vaut 1 si la j-ème seconde de la i-ème fenêtre (de 20s) est 26 sans erreur avec une probabilité PZ, et qui vaut 0 avec une probabilité 1- Pte, la probabilité Pl; étant calculée à l'étape 2. Les variables aléatoires A; sont indépendantes. - étape 4 : exemple du calcul de la courbe C2 : La probabilité ESR5 est donnée par la formule : ESR5 = P(Yi=20)+ P(Yi=19) Un tel résultat, en particulier tel que celui représenté figure 10, était considéré jusqu'à maintenant comme impossible à obtenir. Il permet en particulier de connaître à l'avance avec une excellente fiabilité les caractéristiques de la transmission, sans nécessiter de simulation complète de la liaison physique elle-même, et peut être utilisé pour optimiser la conception et la fabrication du récepteur en fonction des contraintes de qualité de transmission imposées (probabilité en ordonnées sur la figure 10), en choisissant le rapport signal/bruit maximum permettant de respecter ces contraintes. L'invention permet ainsi en particulier d'éviter tout surdimensionnement. L'invention permet d'améliorer d'un facteur de l'ordre de 500 la durée moyenne nécessaire à l'évaluation par rapport à un simulateur complet générant les paquets codés. Une règle empirique indique en effet qu'avec un tel simulateur complet, pour connaitre les PER (taux d'erreur par paquet), il est nécessaire de générer des paquets codés jusqu'à que 100 paquets soient faux. Ainsi, pour exemple, pour un PER de 0,01, le nombre de paquets à générer pour obtenir une seule valeur de PER au cours du temps, il faudra générer 10 000 mots codés. Par ailleurs, les tests montrent que l'évaluation selon l'invention est en moyenne plus rapide que le temps réel. Autrement dit, les performances en PER sur un canal de durée D seront obtenues en un temps de calcul légèrement inférieur à la durée D. Les paramètres principaux de configuration du simulateur sont : le type (lignes-colonnes ; convolutif...) et les paramètres du désentrelaceur permettant de changer les paramètres tels que définis dans la norme DVB-SH (profondeur d'entrelacement, unité d'entrelacement (nombre de symboles)...) ; les paramètres du canal qui permettent à l'utilisateur de choisir un modèle de canal et de modifier les 27 paramètres associés à ce modèle ; les paramètres de la liaison physique qui permettent à l'utilisateur de modifier les paramètres du mode OFDM (intervalle de garde, bande, nombre de porteuses par symbole OFDM) ainsi que le taux de codage. Le simulateur mettant en oeuvre un procédé selon l'invention peut être implémenté sur tout dispositif informatique approprié sous forme d'une application logicielle qui peut se présenter sous forme d'un produit logiciel enregistré sur un support d'enregistrement approprié à la lecture par un lecteur de station informatique (CD-ROM, DVD, clé USB...). REVENDICATIONS1/ - Procédé d'évaluation des performances d'une transmission numérique série codée et modulée sur un canal bruité à atténuation non stationnaire entre : û un émetteur (11) comprenant : û un dispositif de codage adapté pour générer, à partir d'un flux de bits à transmettre, dit flux émis de bits, au moins un flux de mots codés, dit flux émis de mots codés, résultant du codage, selon au moins un procédé de codage prédéterminé, dudit flux émis de bits, un dispositif de modulation adapté pour générer au moins un flux de symboles modulés, dit flux émis de symboles modulés, selon un schéma prédéterminé de modulation, sur au moins un signal de porteuse, chaque flux émis de symboles modulés étant représentatif d'au moins une partie de chaque flux émis de mots codés, û un dispositif d'émission, sur un canal bruité à atténuation non stationnaire, de chaque flux émis de symboles modulés, et un récepteur (12) comprenant : û un dispositif de réception adapté pour recevoir des flux de 20 symboles modulés, dits flux reçus de symboles modulés, correspondant à des flux émis de symboles modulés sur ledit canal, û un dispositif de démodulation adapté pour générer au moins un flux de mots codés, dit flux reçu de mots codés, à partir de chaque flux reçu de symboles modulés, 25 û au moins un dispositif de décodage adapté pour générer un flux de bits, dit flux reçu de bits, par décodage de chaque flux reçu de mots codés, selon un procédé de décodage correspondant à un procédé de codage mis en oeuvre par l'émetteur, procédé dans lequel au moins une valeur représentative d'un taux d'erreurs ER du fluxreçu de bits, est élaborée sans réalisation du décodage, à partir de données numériques mémorisées représentatives des variations au cours du temps de l'atténuation et du bruit du canal, caractérisé en ce que au moins une valeur, dite qualité de transmission Qk, de formule ck2.Es/No, où Ck représente chaque valeur d'atténuation du canal au cours du temps, k étant un indice temporel, Es représente une énergie moyenne par symbole émis et No représente une densité spectrale d'un bruit blanc gaussien sur le canal, et les variations au cours du temps de ladite qualité de transmission Qk, sont élaborées pour chaque symbole reçu du flux reçu de symboles modulés, et en ce que : - dans une première étape (34), une valeur d'information mutuelle Ik est élaborée pour chaque valeur de ladite qualité de transmission Qk, selon une fonction prédéterminée de ladite qualité de transmission Qk, û dans une deuxième étape (35), une moyenne d'information mutuelle est élaborée pour chaque mot codé du flux reçu de mots codés, en réalisant une moyenne des différentes valeurs d'information mutuelle Ik élaborées dans la première étape pour les différentes valeurs prises par ladite qualité de transmission Qk sur le dit mot codé, û dans une troisième étape (38), au moins une valeur du taux d'erreurs ER du flux reçu de bits est élaborée pour chaque mot codé du flux reçu de mots codés, à partir de chaque valeur de la moyenne d'information mutuelle élaborée dans la deuxième étape, et par utilisation de données mémorisées représentatives de variations d'un taux d'erreurs équivalent selon au moins une fonction, dite fonction standard, du rapport signal/bruit, chaque fonction standard étant prédéterminée pour les dispositifs de codage et de décodage sur un canal à bruit blanc gaussien additif. 2/ - Procédé selon la revendication 1, caractérisé en ce que chaque valeur d'information mutuelle Ik est élaborée selon la fonction définie par la formule (I) suivante :C2 ES ) k cc J f (u, v) log2 (f (u, v))dudv I k (Xk , Yk) = loge ( No u=-o0 V=-o0 10 15 ù(u ù ck ES x real(Sm ))2 exp( N° E ) C2 S k N 0 M étant le cardinal de l'alphabet A = ls0, `Sl, • • •, des symboles modulés. 3/ - Procédé selon la revendication 1, caractérisé en ce que chaque valeur d'information mutuelle Ik est élaborée selon la fonction définie par la formule (II) suivante : E Ik(Xk,Yk) _ ùloge(rceck N ) 0 seuil seuil û AuAv f (qAu, rAv) loge (f (qAu, rAv)) q=ùseuil r=ùseuil avec, pour évaluer les différentes gaussiennes de f(u,v) jusqu'à au moins leur probabilité a sur un nombre de points dans le maillage, les valeurs suivantes de seuil et de (tu, 4v) 2 seuil = û ck Es ln(a ck N 7r ) + Maxo par la fonction inverse Ik 1 , û puis chaque valeur du taux d'erreurs ER est obtenue (37) à partir de ladite valeur de rapport signal/bruit équivalent SNReq et desdites données mémorisées représentatives de variations d'un taux d'erreurs équivalent selon une fonction standard prédéterminée pour les dispositifs de codage et de décodage sur un 15 canal à bruit blanc gaussien additif. 9/ - Dispositif informatique comprenant des moyens informatiques de traitement numérique de données, caractérisé en ce qu'il est adapté et programmé pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 8. 10/ - Produit logiciel apte à être chargé en mémoire vive d'un 20 dispositif informatique caractérisé en ce qu'il est adapté pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 8. 11/ - Support d'enregistrement adapté pour pouvoir être lu par un lecteur relié à un système informatique, caractérisé en ce qu'il comprend un programme enregistré adapté pour pouvoir être chargé en mémoire vive d'un système 25 informatique et le programmer pour mettre en oeuvre un procédé selon l'une des revendications 1 à 8.