Brevet FR2751812A1 — publié le 19980130

PROCEDE DE DEMODULATION NUMERIQUE ET DE DECODAGE La présente invention concerne un procédé de démodulation numérique et de décodage. Elle s'applique notamment à un dispositif récepteur mettant en oeuvre une technique de diversité de réception. Les techniques de diversité sont bien connues dans le domaine de la transmission numérique. Parmi ces techniques, on peut citer: - la diversité spatiale, utilisable notamment en transmission radio lorsque plusieurs capteurs de réception sont disposés à des emplacements différents; - la diversité de fréquence lorsque la même information est transmise simultanément sur des fréquences différentes; - la diversité temporelle en cas de répétition de la même information. Ces différentes techniques de diversité peuvent également être combinées entre elles. L'intérêt de ces techniques est qu'elles permettent d'améliorer les taux d'erreur binaire dans les estimations produites par le dispositif récepteur. En contrepartie, elles ont généralement l'inconvénient de requérir des ressources supplémentaires, en termes de bande passante et/ou de complexité des dispositifs émetteur et récepteur. Pour combiner les estimations multiples obtenues par le récepteur à diversité, il existe de nombreuses méthodes, parmi lesquelles on peut citer: - la méthode de sélection consistant simplement à choisir l'observation présentant le meilleur rapport signal/bruit; - la méthode dite "equal gain combining", dans laquelle on prend une décision à partir de la somme des observations après mise en phase; - la méthode dite "maximum ratio combining", dans laquelle on prend une décision à partir de la somme des carrés des observations mises en phase et divisées par la puissance estimée des bruits dont elles sont affectées. Cette dernière méthode fournit un rapport signal/bruit maximal après recombinaison. La majorité des systèmes de transmissions numériques font appel à un codage/décodage canal améliorant leur robustesse aux perturbations induites par le canal de transmission. Les codeurs correcteurs introduisent de la redondance dans l'information transmise, que les décodeurs associés exploitent pour corriger des erreurs dans le signal reçu. Deux grandes familles de codes sont habituellement utilisées: les codes convolutifs et les codes en blocs (Golay, BCH, Reed-Solomon etc...). Les capacités de correction de ces codes sont optimales lorsque les probabilités d'erreur sont décorrélées d'un symbole au suivant, ce qui est le cas par exemple des canaux à bruit blanc additif gaussien. Lorsque cette condition de décorrélation n'est pas remplie, ce qui est notamment le cas des canaux radio à fading de Rayleigh ou de Rice, le codage/décodage canal comporte en outre un entrelacement/désentrelacement consistant en une permutation des symboles destinée à éparpiller, à l'entrée du décodeur correcteur, les erreurs survenant par paquets sur le canal de propagation. Dans la plupart des systèmes, le décodeur canal traite les estimations, dures ou douces, fournies par le démodulateur. Dans le cas d'un récepteur à diversité, il traite les estimations obtenues par recombinaison des estimations multiples. Dans l'article "Channel Coding with Multilevel/Phase Signals" (IEEE Trans. on Information Theory, Vol. IT-28, N 01, janvier 1982, pages 55- 67), G. Ungerboeck a introduit le concept de modulation codée faisant coopérer la modulation/démodulation avec le codage/décodage canal. Le récepteur démodule et corrige les erreurs simultanément en utilisant un supertreillis dont les états tiennent compte d'états de modulation et d'états de codage. Cette méthode permet des gains appréciables car le démodulateur tire parti des bits protégés au niveau microscopique par le code correcteur. Elle a ainsi connu un succès important dans les transmissions filaires et par satellite, pour lesquelles les conditions de propagation ne sont pas trop difficiles. Mais une limitation de cette méthode est qu'elle n'est pas compatible avec la présence d'un entrelacement, car le supertreillis de démodulation prendrait alors une taille non réaliste. Un but de la présente invention est de proposer une nouvelle façon de combiner des estimations multiples de symboles numériques transmis en profitant de la présence de symboles protégés parmi les symboles transmis. L'invention propose ainsi un procédé de démodulation numérique et de décodage, dans lequel un dispositif récep- teur effectue N démodulations distinctes fournissant chacune un jeu d'estimations respectives de symboles successifs d'une trame, certains de ces symboles comportant de la redondance introduite par un codage correcteur d'erreurs appliqué par un dispositif émetteur. Selon l'invention, le dispositif récepteur décode chacun des N jeux d'estimations selon un processus complémentaire dudit codage correcteur d'erreurs, le décodage de chaque jeu d'estimations comportant une mesure d'un degré d'erreur, le dispositif récepteur sélectionnant un des N jeux d'estimations pour lequel le degré d'erreur mesuré est minimal. De cette façon, le résultat final de la démodulation profite de l'information de redondance apportée par le codeur canal. Il suffit que certains seulement des symboles de la trame soient codés avec redondance pour que l'ensemble de la trame démodulée en bénéficie. L'approche peut être comparée à celle d'Ungerboeck, mais elle considère la trame macroscopiquement et non microscopiquement. Ceci permet notamment la présence d'un entrelacement, contrairement aux modulations codées. L'invention s'applique non seulement aux recombinaisons d'estimations multiples obtenues par une technique de diversité, mais aussi au cas o deux au moins des N jeux d'estimations des symboles sont obtenus en démodulant un même segment de signal reçu par des méthodes différentes (un cas de figure o typiquement, les erreurs d'estimation seront souvent corrélées). Dans un mode de réalisation particulier, au moins deux des N démodulations distinctes sont effectuées sur un même segment de signal correspondant à une trame de symboles d'un signal numérique modulé par le dispositif émetteur, ledit segment de signal étant reçu par le dispositif récepteur après transmission du signal numérique modulé par l'intermédiaire d'un canal de transmission, la première de ces deux démodulations comprenant les étapes suivantes: - estimation de premiers paramètres de démodulation à une première extrémité du segment; et - calcul de premières estimations de symboles de la trame sur la base des premiers paramètres de démodulation estimés et du segment de signal parcouru depuis la première extrémité vers une seconde extrémité, et la seconde de ces deux démodulations comprenant les étapes suivantes: - estimation de seconds paramètres de démodulation à la seconde extrémité du segment; et - calcul de secondes estimations de symboles de la trame sur la base des seconds paramètres de démodulation estimés et du segment de signal parcouru depuis la seconde extrémité vers la première extrémité. Cette façon de procéder conduit à des gains appréciables sur le taux d'erreur binaire, à partir d'une seule observation du signal. D'autres particularités et avantages de la présente invention apparaîtront dans la description ci-après d'exemples de réalisation non limitatifs, en référence aux dessins annexés, dans lesquels: - la figure 1 est un schéma synoptique montrant un dispositif émetteur et un dispositif récepteur mettant en oeuvre la présente invention; - la figure 2 est un diagramme montrant la structure de trames de signal dans un exemple de réalisation de la présente invention; - les figures 3 et 4 sont des organigrammes de procédures de démodulation appliquées par le dispositif récepteur dans les deux sens de démodulation; - la figure 5 est un graphique montrant des exemples de vraisemblances obtenues dans chaque sens de démodulation; - les figures 6 à 11 sont des graphiques montrant les performances d'un dispositif récepteur tel que celui de la figure 1; et - la figure 12 montre un autre mode de réalisation d'un dispositif récepteur selon l'invention. L'invention est décrite ci-après dans son application aux radiocommunications numériques entre un dispositif émetteur 10 et un dispositif récepteur 20. Le dispositif émetteur 10 comporte un codeur source 12 (un vocodeur dans le cas d'un système de téléphonie) qui délivre un flux de données numériques xk organisées en trames successives. Dans l'exemple de réalisation illustré par la figure 2, le signal xk est organisé en trames de 126 bits à un débit l/T=8 kbit/s. Un codeur canal 14 traite les bits délivrés par le codeur source pour améliorer la robustesse aux erreurs de transmission. Dans l'exemple de la figure 2, le codeur canal 14 applique un code convolutif CC(2,1,3) de rendement 1/2 aux 26 premiers bits de la trame xk. Les 52+100=152 bits résultants ek sont ensuite soumis à un entrelacement destiné à casser les paquets d'erreurs que peut introduire le phénomène de fading de Rayleigh. Un mot de synchronisation de 8 bits est inséré après chaque trame de 152 bits d'information entrelacés pour former le signal ck que le codeur 14 adresse au modulateur 16. Ce dernier forme le signal radio s(t) qui est amplifié puis appliqué à l'antenne 18 du dispositif émetteur 10. Dans l'exemple considéré, les symboles ck sont binaires (Ck=0 ou 1). La modulation employée est par exemple une modulation GMSK avec un paramètre BT=0,25 (voir K. MUROTA et al: "GMSK modulation for digital mobile radio telephony", IEEE Trans. on Communications, Vol. COM- 29, N 7, Juillet 1981, pages 1044-1050). Le dispositif récepteur 20 comprend un démodulateur 24 recevant le signal capté par l'antenne 22 et amplifié. Le démodulateur 24 délivre deux jeux d'estimations des symboles mis ck. Ces estimations sont notées g, g dans le cas de décisions douces, et d:, d dans le cas de décisions dures. Si les symboles ck sont M-aires et compris entre 0 et M-1, un choix de représentation possible pour l'estimation douce gk ou gR est de la forme: gA= pA.exp(2jndA/M), ou g = PR.exp(2jndR/M) c'est-à-dire que, dans ce cas, son argument 2ndA/M ou 2ndR/M représente la valeur la plus probable dA ou dk du symbole ck, tandis que son module pA ou pR est une mesure de la vraisemblance de cette valeur dA ou d. Dans le cas de symboles binaires (M=2), le nombre gAk ou gRk est réel et appelé "softbit", et son signe 2dA-1 ou 2d-1 donne direc- tement la valeur la plus probable du symbole signé 2ck-1. Le dispositif récepteur 20 comporte un étage de décodage canal 26 dual du codeur canal 14 de l'émetteur. Dans l'exemple considéré, l'étage 26 opère trame par trame la permutation des bits inverse de celle correspondant à l'entrelacement appliqué par l'émetteur, et décode les 52 bits redondants en utilisant le treillis de Viterbi correspondant au code convolutif employé. Dans l'exemple représenté, l'étape 26 comporte deux décodeurs distincts, l'un 26A opérant sur les premières estimations g ou dA fournies par le démodulateur 24, l'autre 26R opérant sur les secondes estimations gk ou dk. Chaque décodeur 26A, 26R restitue des estimations respectives Yk, Yk des bits xk. Comme il est usuel en transmissions numériques, le décodage de Viterbi exécuté par le décodeur 26A ou 26R peut être à décisions dures lorsque le démodulateur 24 fournit seulement les dA ou dR, ou à décisions douces lorsque le démodulateur 24 fournit les gA ou g. Comme le montre la figure 1, le démodulateur 24 comporte un étage radio 30 assurant la conversion en bande de base du signal reçu. Au moyen de deux mélangeurs 32, 34, le signal radio reçu est mélangé à deux ondes radio en quadrature à la fréquence porteuse délivrées par un oscillateur local 36, et les signaux résultants sont soumis à des filtres passe-bas 38, 40 pour obtenir une composante en phase et une composante en quadrature. Ces deux composantes sont échantillonnées et quantifiées par des convertisseurs analogique-numérique 42, 44 à une fréquence au moins égale à la fréquence des bits transmis. On note Rn les échantillons complexes du signal numérique en bande de base délivrés par les convertisseurs 42, 44. Dans l'exemple représenté sur la figure 1, le démodulateur 24 opère selon un algorithme séquentiel pour démoduler des symboles binaires. Dans le cas de la modulation GMSK, on peut réaliser une démodulation séquentielle en utilisant l'approximation suivante pour le signal modulé en bande de base s(t): +00 k s(t) = E j.ak.h(t-kT) (4) k= -x Cette expression correspond à une approximation au premier ordre de la décomposition proposée par P.A. LAURENT dans son article "Exact and Approximate Construction of Digital Phase Modulations by Superposition of Amplitude Modulated Pulses (AMP)", IEEE Trans. on Communications, Vol. COM-34, n 2, Février 1986, pages 150-160. Cet article explique également la méthode de calcul de la fonction h(t), qui, dans le cas de la modulation GMSK avec BT=0,25, correspond à une impulsion de largeur 2T environ centrée sur t=0. Dans l'expression (4), les symboles binaires ak, de valeur 1, correspondent aux bits ck codés différentiellement: ak= akl.(2ck-1). Le canal radio est affecté d'évanouissements correspondant à la somme des trajets multiples correspondant aux différentes réflexions du signal émis sur des obstacles proches ou lointains. La dispersion temporelle de ces trajets étant usuellement de l'ordre de 12 gs, durée faible devant la durée d'un bit (T=125 gs dans l'exemple numérique considéré), on représente le canal de propagation par une variable complexe A(t) correspondant à une atténuation de Rayleigh et un déphasage avec un unique trajet. La fréquence des évanouissements est de 2fd, fd étant la fréquence Doppler associée à la variation de la distance entre l'émetteur et le récepteur: fd=fO.v/c, si f0 est la fréquence centrale du canal, v est la vitesse relative de l'émetteur et du récepteur et c est la célérité de la lumière. On trouve alors pour une vitesse de 100 km/h une fréquence Doppler de l'ordre de 41,67 Hz dans le cas o f0-450 MHz, d'o un évanouissement (83,33 Hz) toutes les 12 ms. Ceci autorise donc plus d'un évanouissement par trame, et surtout une fréquence d'évanouissement supérieure à la fréquence des mots de synchronisation (50 Hz). La présence de ces évanouissements rapides, et plus généralement la variation rapide du canal devant la durée de la trame, imposent une estimation fréquente du canal, et donc un risque important de propagation d'erreurs dues à la rétroaction de la boucle de décision. En effet, en cas d'erreurs sur les symboles binaires décidés lors de la démodulation, ces erreurs vont conduire à des estimations erronées du canal, qui vont elle-même produire de nouvelles erreurs de démodulation. On note Ak=A(kT) (k=0 à 167) les valeurs complexes du canal de propagation échantillonnées à 8 kHz en bande de base. Le canal est en outre affecté d'un bruit blanc additif gaussien B(t) de variance N0/2, noté Bk après échantillonnage et filtrage adapté. Le signal reçu, après filtrage adapté du signal par le filtre 46 de réponse h(t), est alors de la forme: rk=A(kT) jn an H((n-k)) + B(kT) n=-o =Ak [jk-laklH(-Tb)+jkakH(O)+jk+lak+lH(+T)] +Bk o H(t) est la fonction d'autocorrélation connue de la fonction h(t). Dans cette expression, on a fait l'approximation consistant à négliger H(t) pour tJ k 2T, ce qui simplifie les calculs. Les échantillons de sortie rk du filtre adapté 46 sont stockés dans une mémoire 48 pour être traités par le contrôleur 50 du démodulateur 24. Le contrôleur 50 traite le signal filtré rk par segments correspondant chacun à une trame de 168 symboles binaires émis ak (0 k suivant. Le contrôleur 50 effectue la démodulation selon un algorithme séquentiel, dont une première phase est représentée sur l'organigramme de la figure 3. Dans cette première phase, on commence par estimer la réponse complexe du canal au début du segment, puis on démodule ce segment du début vers la fin en mettant à jour à chaque temps-bit l'estimation de la réponse complexe du canal. A l'initialisation 60 de cette première phase, les bits bA et bA sont respectivement pris égaux aux symboles binaires connus a0 et a1, et l'index k est initialisé à 2. A l'étape 62, l'index k est comparé à 8, c'est-à-dire à la longueur du mot de synchronisation. Si k unité à l'étape 70 avant de revenir à l'étape 62. L'estimation du canal au début de la trame est terminée lorsque k=8 au test 62. On dispose alors de l'estimation AA obtenue grâce à la connaissance du mot de synchronisation. Pour chaque valeur de kÄ8, le softbit sA est estimé à l'étape 72 selon: sA = Re (rk.Ak_2.jk) (6) et l'estimation bA du bit ak est obtenue par le signe du softbit sA. A l'étape 72, le contrôleur 50 procède également au décodage différentiel pour estimer les bits ck par l'opération gA = SA_.SA dans le cas d'estimations douces (avec sA = Re(r7.AA.j-7) pour k=8), et par l'opération dA = (l+bA_1.bA)/2 dans le cas d'estimations dures. Ayant obtenu ce bit bA, le contrôleur 50 réestime le canal à l'étape 66 comme exposé précédemment. La démodulation dans le sens aller est terminée lorsque k=167 lors du test 68. On voit sur la figure 3 qu'une erreur faite sur un bit bA à l'étape 72, due par exemple à un évanouissement du canal ou à un bruit impulsif, provoque des distorsions dans les estimations instantanées VA_1, V et A+1 faites dans les trois étapes 66 suivantes, et conduit ainsi à des erreurs d'estimation du canal qui se propagent pendant un certain temps du fait du filtre de lissage. Ces erreurs dans les Ak peuvent à leur tour générer d'autres erreurs d'estimation des bits. La figure 5 montre ainsi, dans le cas o le signal reçu à une énergie évoluant selon la courbe E en trait mixte (avec un évanouissement de canal survenant à l'instant k0), que la vraisemblance sAI des estimations (courbe en trait interrompu) est bonne avant l'évanouissement, mais met ensuite un certain temps à retrouver des valeurs en rapport avec l'énergie E du signal reçu. Pour améliorer les performances dans la période suivant l'évanouissement, le contrôleur 50 procède à une autre démodulation du segment de signal correspondant à la trame de 168 bits depuis la fin du segment vers le début. Ceci permet d'obtenir des vraisemblances J s telles que celles représentées par la courbe en trait plein sur la figure 5. On voit que les performances du démodulateur seront améliorées si on privilégie les softbits sA avant l'évanouissement et les softbits sR après l'évanouissement. La démodulation retour s'effectue dans une seconde phase semblable à la première, dont l'organigramme est représenté sur la figure 4. Dans cette seconde phase, on commence par estimer la réponse complexe du canal à la fin du segment, puis on démodule ce segment de la fin vers le début en mettant à jour à chaque temps-bit l'estimation de la réponse complexe du canal. A l'initialisation 160 de cette seconde phase, les bits bR67 et bR66 sont respectivement pris égaux aux symboles binaires connus a7 et a6, et l'index k est initialisé à 165. A l'étape 162, l'index k est comparé à 159. Si k>159, le bit bR est pris égal au bit connu ak_160 du mot de synchronisation à l'étape 164, puis on procède, à l'étape 166, à une estimation instantanée VR+1 du canal de Vk+ propagation, en effectuant la division complexe: R rk+l (7 k+1 k+2 R k+l R k R k 2bkR+2FH(+T)+j bkR+ lH(O) + bkH(-T) Un filtrage des estimations instantanées VR permet m de lisser les effets du bruit gaussien pour fournir l'estimation AR+, servant à la démodulation des bits. Dans l'exemple représenté sur la figure 4, ce filtrage est simplement le calcul de la moyenne arithmétique des six dernières estimations instantanées Vm. Après l'étape 166, l'index k est comparé à 0 à l'étape 168. Tant que k>0, l'index k est décrémenté d'une unité à l'étape 170 avant de revenir à l'étape 162. L'estimation du canal à la fin de la trame est terminée lorsque k=159 au test 162. On dispose alors de l'estimation A161 obtenue grâce à la connaissance du mot de synchronisation. Pour chaque valeur de k:159, le softbit s est estimé à l'étape 172 selon: SR = Re (rk.Ak+2.jk) (8) et l'estimation bR du bit ak est obtenue par le signe du softbit sk. A l'étape 172, le contrôleur 50 procède également au décodage différentiel pour estimer les bits ck+l, par l'opération gR+l = SR+l.sR dans le cas d'estimations douces (avec sR60 = Re(r160.A161.j160) pour k=159), et par l'opération dR+l = (l+bR+l. bR)/2 dans le cas d'estimations dures. Ayant obtenu ce bit bR, le contrôleur 50 réestime le canal à l'étape 166 comme exposé précédemment. La démodulation dans le sens retour est terminée lorsque k=0 lors du test 168. Dans l'exemple considéré ci-dessus, les paramètres de démodulation réestimés lors du parcours du segment démodulé dans chaque sens se limitent à la réponse complexe Ak du canal de propagation. On comprendra qu'ils pourraient inclure d'autres paramètres tels que des paramètres représentatifs du bruit observé sur le canal de transmission. On peut ainsi calculer pour chaque sens de démodulation une moyenne quadratique des écarts V1_l-A' -k-1 (étape 66) ou VR±AR, (étape 166), pour estimer la k+l k+1 puissance instantanée du bruit N0A, N0 dans chaque sens de démodulation. On peut alors normaliser la valeur du softbit sA ou sk en la divisant par cette moyenne quadratique. Les estimations de puissance N0A, N0 peuvent être constantes sur la trame considérée; ce sont alors, par exemple, des moyennes des A_- VA_12 et des IA+1- V+1l 2 calculées sur l'ensemble de la trame. Si ces moyennes sont obtenues sur des fenêtres glissantes ou par filtrage, les estimations de la puissance du bruit peuvent être instantanées, c'est-à- dire dépendre de l'index k. Bien qu'on ait représenté deux décodeurs canal distincts 26A, 26R pour faciliter la lecture de la figure 1, on comprendra qu'il est également possible de prévoir un décodeur unique assurant successivement les deux décodages. Le premier décodage donne lieu à des bits estimés yk, et le second décodage à des bits estimés yR. Chaque décodeur canal fournit en outre une mesure QA, QR du degré d'erreur observé lors du décodage de la trame d'estimations reçue. Dans le cas, par exemple, o le décodeur canal utilise un treillis de Viterbi pour un code convolutif (voir "The Viterbi Algorithm" par G. D. Forney, Proc. IEEE, Vol. 61, n 3, mars 1973, pages 268-278), la mesure QA (ou QR) peut être: - le nombre des symboles dont la valeur M- aire a été corrigée lors du décodage; - une fonction décroissante de la plus grande métrique ayant permis de sélectionner un trajet dans le treillis lors du décodage de la trame; - ou une autre mesure, pouvant dépendre de l'algorithme de décodage particulier mis en oeuvre, de la qualité observée par le décodeur. En comparant les mesures de degré d'erreur QA et QR, le dispositif récepteur sélectionne l'un ou l'autre des deux jeux d'estimations gA et gR (ou dAk et dR), à savoir celui pour lequel la mesure du degré d'erreur est la plus faible, comme schématisé sur la figure 1 par le soustracteur 96 et le commutateur 98 qui retient les bits décodés correspondants YAk ou y selon les valeurs relatives de QA et QR. Les bits retenus Yk sont fournis à un décodeur source 28 qui remet en forme l'information transmise. Le dispositif récepteur choisit ainsi celui des deux sens de démodulation qui semble le meilleur du point de vue du décodage canal. Il est à noter qu'il n'est pas nécessaire que le codage canal implique des bits de redondance pour chacun des bits émis. Ainsi, dans le cas précédemment évoqué en référence à la figure 2, o seuls 26 bits sur 126 sont codés avec redondance, ces 26 bits suffisent à porter un jugement sur la qualité des estimations démodulées. En d'autres termes, les 152 bits de la trame profitent de l'information de redondance contenue dans les 52 bits délivrés par le codeur convolutif. Comme il est usuel en transmissions numériques, la forme particulière du codage à redondance appliqué, la permutation utilisée pour l'entrelacement et la forme d'onde modulante doivent faire l'objet d'une optimisation conjointe pour obtenir les meilleures performances dans chaque cas concret. Cette optimisation peut, de façon connue, être effectuée au moyen de simulations sur ordinateur du comportement du canal. Les figures 6 à 8 illustrent les performances obtenues avec un procédé de démodulation selon la figure 1. Ces résultats ont été obtenus par simulation dans le cas de trames de signal numérique telles que représentées sur la figure 2 modulées en GMSK avec BT=0,25. Le démodulateur différait de celui représenté sur la figure 1 en ce que la démodulation n'était pas séquentielle mais selon un treillis de Viterbi ayant 8 états et 8 filtres adaptés, avec poursuite du vecteur d'onde. Ce treillis correspond à une mémoire de L=3 symboles dans une décomposition de Rimoldi du signal GMSK (voir B.E. RIMOLDI "A Decomposition Approach to CPM", IEEE Trans. on Information Theory, Vol. 34, n 2, mars 1988, pages 260-270). Le treillis délivrait ainsi directement les estimations dures d:, dR correspondant aux bits ck, sans décodage différentiel. L'entrelacement des bits ek (k=0,...,151) effectué par le codeur canal 14 pour fournir les bits ck était opéré selon ek=c8+Tab(k) avec [Tab(0),Tab(l),.. .,Tab(151)] = [56, 0, 112, 8, 55, 16, 104, 136, 24, 120, 88, 32, 72, 96, , 128, 108, 48, 80, 60, 4, 144, 116, 12, 68, 132, 20, 124, 84, 28, 76, 140, 36, 148, 92, 44, 100, 52, 2, 106, 58, 10, 66, 114, 18, 122, 130, 26, 74, 138, 34, 82, 90, 42, 146, 110, 50, 98, 62, 6, 70, 118, 14, 126, 134, 22, 78, 86, 30, 142, 94, 38, 150, 54, 46, 102, 57, 1, 113, 9, 65, 17, 105, 129, 25, 121, 137, 33, 73, 89, 41, 81, 109, 49, 145, 61, 5, 97, 117, 13, 69, 133, 21, 125, 85, 29, 77, 141, 37, 149, 93, , 101, 53, 3, 67, 107, 11, 123, 59, 19, 75, 115, 27, 83, 131, 35, 147, 139, 43, 99, 91, 51, 71, 111, 7, 127, 63, 15, 79, 119, 23, 143, 135, 31, 151, 87, 39, 103, 95, 47, 64]. Le décodeur canal 26A ou 26R effectue la permutation inverse pour le désentrelacement, puis décode les 52 premiers bits obtenus avec un treillis de Viterbi à décisions dures correspondant au code convolutif employé. Le graphique de la figure 6 montre le taux d'erreur binaire obtenu en fonction du rapport signal/bruit Eb/N0 dans le cas o les dispositifs émetteur et récepteur sont statiques. Le graphique de la figure 7 illustre un fading dynamique avec une vitesse relative de 70 km/h. La courbe en traits pleins représente le taux d'erreur binaire observé avec une démodulation dans un seul sens, la courbe en pointillés représente les performances obtenues selon l'invention avec une démodulation aller-retour et une sélection de sens selon la figure 1, et la courbe en traits mixtes les performances théoriques d'un démodulateur idéal (c'est-à-dire ayant une connaissance parfaite et instantanée du canal). La figure 8 montre le taux d'erreur binaire observé en dynamique (70 km/h) avec un interféreur décorrélé présent sur le même canal fréquentiel. En abscisses, la quantité C/Ic représente le rapport des puissances reçues par le dispositif récepteur depuis le dispositif émetteur et depuis l'interféreur. Dans les trois cas, on observe une amélioration sensible du taux d'erreur binaire. Les figures 9 à 11 sont des graphiques semblables, respectivement, à ceux des figures 6 à 8, montrant des résultats obtenus dans des conditions de simulation semblables, différant seulement par le type de modulation employé. Il s'agissait d'une modulation quaternaire (M=4). Les bits ck délivrés par le codeur canal étaient regroupés par paires pour former des symboles quaternaires traités par le modulateur. Les estimations dures des symboles délivrées par le démodulateur dans le sens aller ou retour étaient ensuite décomposées en deux estimations dures des bits ck correspondants. La modulation quaternaire considérée était une modulation à phase continue (CPM) de paramètre BT=0,25 avec l'approximation d'une mémoire de L=3 symboles pour la construction du treillis de démodulation selon la décomposition de Rimoldi. L'indice de modulation était h=l/3, le démodulateur employant un treillis à 48 états et 64 filtres adaptés. La figure 12 montre un autre exemple de dispositif récepteur apte à mettre en oeuvre la présente invention. Ce dispositif 120 fait appel à une diversité en réception qui, dans l'exemple considéré, est une diversité spatiale, le dispositif comportant n antennes 221,.. .,. 22n et n démodulateurs associés 241,...,24n. Chaque démodulateur 24i opère dans un seul sens sur un segment de signal respectif fourni par son antenne 22i (par exemple comme décrit en référence à la figure 3), et délivre des estimations douces respectives g(i) et/ou des estimations dur es respectives dki) pour chaque symbole ck. Le dispositif dispose ainsi de N=n estimations par symbole provenant de segments de signal différents au lieu de N=2 estimations tirées du même segment de signal dans l'exemple de réalisation des figures 1 à 11. Ces n jeux d'estimations sont décodés conformément à la protection apportée par le codeur canal, en parallèle par n décodeurs canal 261,..., 26n (ou encore successivement par un unique décodeur). Le dispositif 120 obtient ainsi n jeux de symboles décodés y(1),...,y n) et n mesures de degré d'erreur Q(1),...,Q(n) associées. Un comparateur 196 détermine le jeu d'estimations g(i) ou dki) pour lequel le degré d'erreur mesuré Q(i) est minimal, et après décodage, ce jeu est finalement retenu par le module de sélection 198 pour alimenter le décodeur source 28. L'invention est bien entendu applicable à d'autres techniques de diversité, ou à des récepteurs combinant une technique de diversité avec une méthode de démodulations multiples telle que celle décrite ci-dessus. REVENDICATIONS 1. Procédé de démodulation numérique et de décodage, dans lequel un dispositif récepteur (20;120) effectue N démodulations distinctes (N> 2) fournissant chacune un jeu d'estimations respectives de symboles successifs (ck) d'une trame, certains de ces symboles comportant de la redondance introduite par un codage correcteur d'erreurs appliqué par un dispositif émetteur (10), caractérisé en ce que le dispositif récepteur (20;120) décode chacun des N jeux d'estimations selon un processus complémentaire dudit codage correcteur d'erreurs, le décodage de chaque jeu d'estimations comportant une mesure d'un degré d'erreur (Q(i)), le dispositif récepteur sélectionnant un des N jeux d'estimations pour lequel le degré d'erreur mesuré est minimal. 2. Procédé selon la revendication 1, caractérisé en ce que le degré d'erreur (Q(i)) mesuré lors d'un décodage est le nombre de symboles modifiés au cours dudit décodage. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que les symboles de la trame sont entrelacés, le dispositif récepteur (20;120) effectuant un désentrelacement avant le décodage de chacun des N jeux d'estimations. 4. Procédé selon la revendication 1,2 ou 3, caractérisé en ce qu'au moins deux des N démodulations distinctes sont effectuées sur un même segment de signal correspondant à une trame de symboles d'un signal numérique modulé par le dispositif émetteur (10), ledit segment de signal (r(t)) étant reçu par le dispositif récepteur (20) après transmission du signal numérique modulé (s(t)) par l'intermédiaire d'un canal de transmission, en ce que la première de ces deux démodulations comprend les étapes suivantes: - estimation de premiers paramètres de démodulation (AA) à une première extrémité du segment; et - calcul de premières estimations (g) de symboles de la trame sur la base des premiers paramètres de démodulation estimés et du segment de signal parcouru depuis la première extrémité vers une seconde extrémité, et en ce que la seconde de ces deux démodulations comprend les étapes suivantes: - estimation de seconds paramètres de démodulation (AR) à la seconde extrémité du segment; et calcul de secondes estimations (gR) de symboles de la trame sur la base des seconds paramètres de démodulation estimés et du segment de signal parcouru depuis la seconde extrémité vers la première extrémité. 5. Procédé selon la revendication 4, caractérisé en ce que les premiers paramètres de démodulation sont réestimés au moins une fois lors du parcours du segment à partir de la première extrémité, et les seconds paramètres de démodulation sont réestimés au moins une fois lors du parcours du segment à partir de la seconde extrémité. 6. Procédé selon la revendication 4 ou 5, caractérisé en ce que les premiers et seconds paramètres de démodulation comprennent chacun au moins un paramètre (AM, A) représentant la réponse du canal de transmission. 7. Procédé selon la revendication 6, caractérisé en ce que le dispositif récepteur (20) estime lesdits paramètres représentant la réponse du canal de transmission aux extrémités du segment sur la base de séquences de synchronisation incluses aux extrémités des trames de signal numérique. 8. Procédé selon l'une quelconque des revendications 4 à 7, caractérisé en ce que les premiers et seconds paramètres de démodulation comprennent chacun au moins un paramètre relatif au bruit observé sur le canal de transmission. 9. Procédé selon la revendication 8, caractérisé en ce que les premiers paramètres de démodulation comprennent la puissance du bruit dont l'estimation (N0A) est utilisée pour normaliser les premières estimations des symboles de la trame, et en ce que les seconds paramètres de démodulation comprennent la puissance du bruit dont l'estimation (N0O) est utilisée pour normaliser les secondes estimations des symboles de la trame. 10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce qu'au moins deux des N démodulations distinctes sont effectuées sur deux segments de signal respectifs, reçus par le dispositif récepteur (120) selon une technique de diversité.