FR 2633470 A1 19891229 FR 8808525 A 19880624 Procédé de démodulation et démodulateur de signaux numériques d'amplitude d'enveloppe constante modulés en phase et/ou en fréquence de façon continue. La présente invention concerne un procédé de démodula- tion et un démodulateur de signaux numériques d'amplitude d'en- veloppe constante modulés en phase et/ou en fréquence de façon continue. Elle concerne notamment les procédés de démodulation de signaux radioélectriques porteurs d'informations numériques, et leurs utilisations dans les systèmes de transmission de l'in- formation par paquets ou en évasion de fréquence. Dans ces procédés, les données à transmettre sont conver- ties sous la forme de suites de symboles de durée constante formant l'intervalle Baud de la transmission pour moduler en fréquence ou en phase de façon continue une onde porteuse d'amplitude d'enveloppe constante. Leur Intérêt réside surtout dans le fait qu'ils conduisent à des occupations en fréquence plus faibles que celles des systèmes de transmission connus sous les abréviations FSK ou PSK des désignations anglo saxonnes "Frequency Shift Keying" et "Phase Shift Keying" qui utilisent des symboles de forme rectangulaire. Un exemple de système mettant en oeuvre un procédé de modulation à porteuse modulée en phase de façon continue et d'amplitude d'enveloppe constante est connu de la demande de brevet français n 82 08043. Selon le système de modulation décrit, le modulateur d'émission comporte, un codeur qui fait correspondre aux don- nées ou symboles à transmettre des valeurs discrètes de sauts de phase et un synthétiseur de symboles qui élabore en fonction des informations fournies par le codeur une forme de variation de phase continue à pente minimale, exprimée en fonction du temps par une loi polynomiale. Le démodulateur de réception utilise une démodulation synchrone et procède par corrélation pour récupérer les instants caractéristiques de la modulation et de la dérive en fréquence entre les porteuses d'émission et de réception. Ce système présente l'avantage d'occuper une largeur de bande de fréquence faible eu égard au débit des Informations qu'il transmet et de permettre par exemple des transmissions de données à 16 K bits/s dans des canaux larges de 25 KHz avec une atténuation à l'extérieur du canal de 50 dB. Cependant, le démodulateur s'avère compliqué à réaliser lorsque les systèmes de transmission de données ont à transmet- tre des données par paquets ou en évasion de fréquence du fait notamment de la complexité du dispositif d'estimation de l'hor- loge symbole qui est employé. Comme celui-cl possède une fonc- tlon spécifique pour estimer la composante continue et une fonc- tion spécifique pour estimer l'horloge symbole conditionnée par l'estimation de la composante continue, les effets dus au bruit naturel de la transmission et à la dérive en fréquence de l'os- cillateur pilote de réception relativement à celle de l'oscilla- teur d'émission ne sont pas dissociés. Il en résulte que les performances du circuit de synchronisation vont en se dégradant à mesure que la transmission est entachée de bruit. Le but de l'invention est de palier les Inconvénients précités. A cet effet, l'invention a pour objet un procédé de démo- dulation de signaux numérique d'amplitude d'enveloppe constante et modulés en phase et/ou en fréquence de façon continue, les signaux numériques étant composés par des suites de symboles de même durée T et juxtaposés les uns aux autres dans le temps, les symboles étant transmis au moyen de sauts de phase et/ou de fréquence, caractérisé en ce qu'il consiste: - à détecter la phase du signal modulé sous la forme d'un signal d'amplitude variable fonction des variations de phase du signal modulant; - à échantillonner, suivant une période d'échantillonnage égale à T/N. la phase du signal modulé détecté, à l'intérieur de chaque durée symbole T au moyen d'un peigne d'impulsions d'échantillonnage obtenu à partir de N impulsions d'horloges régulièrement espacées d'un intervalle de temps égal à T/N; - à caleuler l'angle de phase différentielle d0 sépa- rant chaque échantillon de phase d'un symbole courant détecté, de l'échantillon de même rang situé dans le symbole détecté au temps symbole précédent; - à coder chaque angle de phase différentielle dO en lui attribuant un coût de codage dépendant de son amplitude; - à calculer les sommes des coûts de codage pour les échantillons de même rang appartenant à un nombre déterminé L de symboles consécutifs; - à identifier, parmi les N sommes obtenues celle qui donne le résultat le plus grand en mémorisant le rang occupé dans chaque symbole par les échantillons correspondants - à prendre comme horloge de synchronisation symboles, pour définir les instants de synchronisation du démodulateur, l'horloge qui dans le peigne d'impulsions possède le même rang que le rang de l'échantillon mémorisé; - et à décoder l'angle de phase différentielle qui pos- sède le même rang que l'échantillon mémorisé.. Elle a également pour objet, un démodulateur pour la mise en oeuvre du procédé précité. L'invention a pour avantage qu'elle permet une acquisi- tion rapide de l'horloge symbole compatible avec la longueur des paliers caractéristiques des transmissions en évasion de fré- quence. Ce résultat est obtenu grâce à la méthode relativement simple utilisée pour récupérer l'horloge symbole. Cette simplicité résulte du fait que l'estimation de l'horloge symbole est effectuée selon un seul et même concept quelle que soit la cause de la fluctuation 'du rythme, qui peut être due soit à un retard constant entre l'émission et la récep- tion, soit à une dérive en fréquence des fréquences de référen- ces des postes, soit à une dérive en fréquence due à l'effet Doppler, soit encore au bruit gaussien entachant normale- ment chaque transmission. D'autre part, les performances par rapport aux systèmes antérieures doivent normalement être améliorées du fait qu'un véritable filtrage en phase et non en fréquence est opéré par le système proposé et qu'une corrélation plus fine est obtenue pour la recherche de la synchronisation. D'autres caractéristiques et avantages de l'invention apparaîtront ci-après à l'aide de la description faite en regard des dessins annexés qui représentent: La figure 1 une courbe pour rappeler le principe de la modulation de phase continue; La figure 2 le principe de l'échantillonnage mis en oeu- vre par l'invention pour l'obtention d'échantillons de phase 1 5 différentielle; La figure 3 une représentation connue d'un diagramme de l'oeil correspondant A une modulation A phase continue à 5 états; Les figures 4 à 6 une illustration des fonctions coûts utilisées comme critère de synchronisation; La figure 7 un mode de réalisation d'un démodulateur selon l'invention; La figure 8 un mode de réalisation du bloc de calcul de phase différentielle de la figure 7; La figure 9 un mode de réalisation du bloc de mémorisa- tion de la phase différentielle de la figure 7; La figure 10 un mode de réalisation du bloc de calcul de l'horloge de synchronisation de la figure 7; La figure 11 un graphe en trois dimensions de répartition de la fonction coût tenant compte de la dérive en fréquence de la porteuse; La figure 12 une vue de dessus du graphe de la figure 11; La figure 13 un mode de réalisation du démodulateur de l'invention muni d'un dispositif de correction de dérive. La courbe de la figure i retrace la loi de variation de phase 0(t) d'un signal à enveloppe constante obtenu suivant le principe connu de la modulation de phase à phase continue. Sur cette courbe les points marqués O, +11/2, -Â1, -11/2 etc... repré- sentent des états pris par la phase 0(t) du signal pour trans- mettre les différents symboles Sn qui caractérisent une transmission, chaque symbole étant représenté par un état de phase déterminé. Comme décrit dans le brevet français précité, 11 est donné la même valeur de phase et la même pente à l'origine de chaque symbole S et à l'extrémité du symbole Sn 1 qui n n-i le précède, d'autre part, la pente maximale de la courbe repré- sentative de la phase 0(t) est minimale entre les états de phase successifs de façon à réduire au maximum la largeur du spectre en fréquence occupée. En adaptant une solution poly- nomiale du troisième degré, la variation de phase entre chaque symbole Sn et Sn-1 est très lente, et la largeur du spectre en fréquence occupée par la modulation est moindre que celle occupée par les modulations type FSK ou MSK citées précédemment. Suivant l'invention la phase 0t) du signal est échan- tillonnée de la façon représentée à la figure 2 par un peigne d'échantillonnage dont chaque dent représente une horloge de synchronisation possible. La phase différentielle dO(tnk) mesurée entre deux échantillons successifs est égale à dO(tnk) 0(tn,k) - 0(tn,k-T) (1) avec tnk = nT * kte avec te = T/N o: N est le nombre total d'horloges de synchronisation - t est la période d'échantillonnage de la phase différen- tielle - T est la période symbole et k est un nombre entier compris entre O et N-1. Pour réaltsér de façon fiable le décodage des symboles, l'horloge symbole est estimée avec précision en tenant compte du fait que l'instant de décision choisi pour le décodage est régu- lièrement espacé dans le temps d'une période symbole et qu'il ne coïncide pas nécessairement avec l'instant de synmchronisation. L'instant de synchronisation t correspond à l'empla- o cernent du diagramme de l'oeil Figure 3 o l'information de synchronisation est la plus dense. Cet instant est défini par un nombre déterminé de valeurs d'échantillons de phase différen- tielle a = d0(to) pour lesquels la différence minimale prise en valeur absolue entre deux échantiilons est la plus grande possible pour que l'immunité au bruit sclt la plus grande. Cependant, comme les valeurs exactes de ai ne sont jamais atteintes à cause du bruit et des imperfections technolo- giques, l'instant de synchronisation ne peut être évalué qu'en estimant la distance qui sépare les valeurs des échantillons mesurés des valeurs exactes des échantillons ai qui correspon- dent réellement à l'instant de synchronisation. Cette distance mesure la probabilité de trouver effectivement l'instant de synchronisation et elle est définie, dans ce -pai suit, à l'aide d'une fonction coût dont un exemple de représentation graphique est montré à la figure 4. Cette fonction possède des maximums d'amplitude A lorsque l'échantillon ai de phase différentielle dM mesuré correspond exactement aux instants de synchronisa- tion et des valeurs nulles pour les échantillorns de phase diffé- rentielle dd situés aux milieux des interval!es compris entre les instants exacts de synchronisation successL's. Sur les repré- sentations des figures 4 et 5 les variations de la fonction coût entre chaque point de la courbe d'amplitude nulle et les deux points directement voisins d'amplitudes maximins A sont linéai- res, les droites correspondantes ayant des pentes opposées. Cependant selon d'autres variantes de réalisations de l'inven- tion, le mode de variation linéaire pourra être remplacé par des modes de variation polygoniaux. L'amplitude maximum A du coût est fixée par la largeur d'une fenêtre statistique comptée en nombre de symboles et éven- tuellement des limitations Imposées par le cahier des charges (longueur d'un palier EVF par exemple). Pour les échantillons a et aimax correspon- dant aux symboles extrêmes de la fenêtre, la symétrie de la fonction coût est celle qui est représentée à la figure 5. Une variante tenant compte d'une non équiprobabllité des échantillons ai est représentée à la figure 6. Selon cette variante, chaque fonction coût relative à un échantillon ai d'un symbole Si présente une amplitude maximum Ai qui lui est propre. Cette amplitude A1 est définie par la relation A1 = A x P1/P o A désigne l'amplitude maximum possible et P1 désigne la probabilité d'obtenir l'instant de synchronisa- tion t confondu avec l'instant d'échantillonnage de la o phase différentielle d0. Lfestimation de l'horloge a lieu en choisissant celle qui est la plus proche de la position dans le temps de la suite des échantillons d'ordre ai de chaque symbole dont le total des coûts accumulés sur L symboles consécutifs et mesurés à l'aide de la fonction coût décrite précédemment est maximal. Cette estimation statistique est effectuée dans le cas d'un mode de transmission par paquets ou en évasion de fréquence, sur le nombre L des symboles d'un paquet ou d'un palier en évasion de fréquence. Pour chaque horloge symbole d'indice k d'un ensem- ble de N horloges correspond une estimation statistique décrite par une relation de la forme: S (k)=C(d0(tno k))C(d0 (tno+ 1,k) +...+C(d0(tnoL,k)) Si ke désigne l'horloge de synchronisation estimée, celle-cl est déterminée par la somme Sk qui présente le résultat le plus grand. Suite à cette estimation deux actions sont menées La première action consiste à corriger l'horloge du démo- dulateur pour permettre le recalage de celle-cl sur l'instant d'horloge k précédemment estimé. La deuxième action consiste à décider du décodage. L'action de correction de l'horloge de démodulation est réalisée via un séquencement interne du démodulateur et s'exé- cute dans le sens qui ramène à l'horloge symbole prise comme référence, à l'instant Initial, cet instant pouvant être choisi comme le milieu de la plage de temps couverte dans un symbole par toutes les horloges du peigne. Comme l'horloge de décision du décodage n'est pas nécessairement confondue avec l'horloge de synchronisation mais peut présenter un écart constant avec l'horloge de synchronisation, l'instant de décision du décodage doit également être déterminé par une translation dans le même sens que celui effectué pour l'instant de synchronisation. Un exemple de réalisation d'un démodulateur selon l'inven- tion est représenté à la figure 7. Il comprend un bloc 1 de calcul de la phase différentielle d0, piloté par une horloge pilote 2, et couplé à un bloc de mémorisation 3, à un bloc de calcul de l'horloge de synchronisation 4 et à un décodeur 5. Le bloc de calcul de l'horloge de synchronisation 4 est couplé d'autre part, au bloc de mémorisation 3 et à un bloc de séquence- ment et d'asservissement 6. Le bloc de calcul de la phase différentielle 1 est repré- senté en détail sur la figure 8. Il comprend un étage de mise en forme 7, relié à un registre à décalage 8. Le registre à déca- lage 8 est couplé à une boucle à verrouillage de phase'composée par un registre tampon 9, un comparateur d'états 10 et un divi- seur programmable 11. Le comparateur d'état 10 est relié à un registre à décalage 13 au travers d'un compteur-décompteur de phase 12. Le compteurdécompteur de phase 12 est relié à une première entrée d'opérande d'un circuit soustracteur 14. La deuxième entrée d'opérande du circuit soustracteur 14 est reliée à ce registre à décalage 13. Le bloc de calcul de la phase différentielle 1 constitue l'entrée fonctionnelle du démodulateur. Son rôle est de mettre sous une forme binaire exploitable par le démodulateur, le si- gnal radio FI en bande de base Issu de systèmes de transposi- tion de fréqi:ence classiques non représentés, obtenu de façon connue en sortie des étages de fréquence intermédiaire d'un récepteur. Cette fonction est assurée de façon connue par l'étage de mise en forme 7 au moyen par exemple d'un circuit comparateur de tension, non représenté, à même de travailler à des fréquen- ces et des tensions compatibles avec le signal d'entrée. A cet égard deux types de mise en forme peuvent éventuellement être effectués. Une première mise en forme peut consister à comparer le signal d'entrée à la valeur moyenne de ce dernier et dans ce 1i cas, le comparateur doit être pourvu d'un système à hystérésis pour évlter des instabilités dans le cas o le signal d'entrée est faible. Une deuxième mise en forme peut consister à compa- rer le signal d'entrée FI de manière différentielle. L'avantage de ce dernier type de comparaison est qu'il permet une transcrip- tion parfaite des points o le niveau du signal d'entrée s'an- nule. La deuxième fonction assurée par le bloc de calcul 1 est le calcul de la phase différentielle d0 proprement dit. La première opération qui est exécutée consiste à créer un signal numérique à partir de l'horloge pilote, capable d'être comparé en phase avec le signal d'entrée. Ceci est assuré de façon con- nue par une boucle à verrouillage de phase. Cette boucle se compose sur la figure 8 des éléments 8 à 11. Le registre 8 contient les échantillons du signal d'en- trée et son contenu est renouvelé avec une précision et périodi- cité égales à celles du diviseur programmable 11. Le registre tampon 9 permet de stabiliser deux échantil- lons du signal pris successivement et ceci à la fréquence de l'horloge locale. Le comparateur d'états 10 assure une comparaison entre les grandeurs des échantillons contenus dans le registre 9. Il fournit de la sorte sur sa sortie, deux informations qualifiant la croissance ou la décroissance en phase du signal d'entrée. Son évolution temporelle a une précision qui est fixée par le rapport de l'horloge pilote et du signal d'entrée. Le pas d'échantillonnage est par exemple de 10 , si le rapport est le 1/36e de la période du signal d'entrée. Les deux informations d'avance de phase (+10 dans l'exemple) et de retard de phase (-10 dans l'exemple) fournies par le comparateur d'états 10 sont appliquées à l'entrée du circuit 11. Le compteur-décompteur 12 effectue la comptabilisation des avances et retards de phase fournis par le comparateur d'états 10. Une avance de phase augmente de une unité le conte- nu du compteur 12 et un retard diminue son contenu de une unité. La sortie du compteur-décompteur 12 peut fournir dans ces conditions la phase absolue 0(t) du signal d'entrée depuis un instant quelconque d'utilisation à la valeur zéro du contenu du compteur-décompteur 12. Le registre à décalage 13 a une structure de registre pile, encore connue sous la désignation anglo-saxonne FIFO qui est l'abréviation de "First in firs out", et forme avec le circuit soustracteur 14, l'élément de calcul de phase différentielle proprement dit. Le circuit soustracteur 14 reçoit en addition les échantillons de phase 0(t) successifs provenant du compteur/décompteur 12 et reçoit en retranche- ment les échantillons de phase successifs O(tT) Issus du registre à décalage 13 dont le nombre fixe de décalages est égal au rapport de la fréquence d'échantillonnage du signal modulant (He) à la fréquence du signal modulant lui-même (H symbole). L'écart de phase dO(t) obtenu à la sortie du circuit soustracteur 14 vérifie la relation dO(t)=0 (t)-0(t-T) o t désigne l'instant présent et T la période du signal modu- lant, cette période étant la même que la période symbole décrite précédemment. Les échantillons de phase différentielle d0(t) sont transmis à l'entrée du bloc de mémorisation 3. Ce bloc, qui est représenté à la figure 9 comprend une mémoire vive à accès Il aléatoire 15 adressée par un sélecteur d'adresse 16 et par la sortie d'un circuit soustracteur 17. La mémoire 15 permet l'enregistrement et la relecture des échantillons de la phase différentielle d0 () obtenus du bloc de calcul 1. Le sélec- dlúú6retteUe d(t) otnsd teur 16 comporte deux entrées. Une première entrée reçoit le numéro de l'horloge de décodage, l'autre est reliée à la sortie d'un compteur d'échantillon 18. Le circuit soustracteur 17 com- porte deux entrées d'opérandes. Une première entrée d'opérande marquée (+) est reliée à la sortie d'un compteur de symbole 19. La deuxième entrée d'opérande marquée (-) reçoit sous la forme d'un mot binaire le nombre de symboles L sur lequel porte les calculs statistiques décrits précédemment. Le compteur de symbole 19 progresse au rythme d'impulsions fournies par un circuit diviseur 20, ce circuit 20 effectuant une division par le nombre N d'échantillons contenus dans un symbole de la fré- quence d'échantillonnage Fe. Pour retrouver les informations, notamment le rythme et les caractères de la modulation, il est nécessaire que le signal modulant reçu et la phase différentielle d0(t) calculée, subissent un sur-échantillonnage à décalage constant. Cet échan- tillonnage est choisi N fois supérieur à la fréquence symbole et le décalage entre échantillons est de 360/N degrés de la phase du signal modulant. Dans ce but, le bloc de calcul de l'horloge de synchroni- sation 4 opère par traitement glissant sur les symboles consécu- tifs correspondant à la longueur du traitement statistique. Pour satisfaire à ces conditions, la mémoire 15 est orga- nisée pour stocker N x L échantillons de la phase différentielle et son fonctionnement est simnilaire à celui d'un registre à décalage. Sa structure est de préférence celle d'une mémoire à accès aléatoire pour permettre un fonctionnement souple en éva- sion de fréquence ou en mode de transmission de l'information par paquets. L'adressage de la mémoire 15 est réalisé à l'aide des circuits 16 à 20. Le sélecteur 16 assure la commutation entre deux modes distincts de fonctionnement en écriture ou en lecture de la mémoire 15. Le mode de fonctionnement en écriture permet l'enregistre- ment dans la mémoire 15 des échantillons de phase différentielle d0t provenant du bloc de calcul 1. L'adressage est défini pour permettre à chaque dernier échantillon d0(t) qui arrive à un instant to, de pren- dre la place de celui qui a été enregistré au temps t = t0-LT c'est-àdire L symboles plutôt. Le repérage de la o case mémoire concernée est fourni par les contenus des comp- teurs 18, 19 et 20. Dans le mode de fonctionnement en lecture, l'échantillon de la phase différentielle d0'(t) occupant la place du prochain à venir est dirigé vers le bloc de calcul de l'horloge de synchronisation 6. Cependant, le mode de fonctionne- ment en lecture est aussi provoqué pour transmettre l'échantil- lon de phase différentielle d0(t) au décodeur 5 dont le (t) numéro d'horloge a été reconnu ayant la plus forte probabilité d'être celui correspondant. à l'horloge de synchronisation. Le fonctionnement combiné du sélecteur 16 et du soustracteur 17 permet de faire travailler la mémoire en mode d'adressage par paquets, le sélecteur 10 étant commandé pour sélectionner l'hor- loge de décodage, l'entrée notée (-) du sélecteur 11 recevant un nombre binaire compris entre zéro et L. Le bloc de calcul de l'horloge de synchronisation repré- senté à la figure 10 se compose de 4 parties. La première partie comprend un dispositif de codage 21 formée par une mémoire morte programmable ou tout dispositif équivalent permettant le codage des échantillons de phase diffé- rentielle. Le codage est effectué en appliquant la fonction coit décrite précédemment A chaque échantillon de phase dIfféren- tielle d0(t. La mémoire du dispositif de décodage 21 est adressée d'une part, par l'échantillon de phase arrivant d0(t) et par l'échantillon déjà mémorisé dans la mémoire 15 au temps t = to-LT. La deuxième partie comprend un circuit additionneur 22 et une mémoire vive 23. La mémoire 23 a un rôle d'accumulateur pour le circuit additionneur 22. Elle comprend une pile de N registres accumula- teurs organisée dans le mode FIFO. Le circuit additionneur 22 ajoute à chaque donnée lue sortant par le dernier registre de la pile les transformées (C.d0(t) et Cd0'(t)) par les fonc- tions coût décrites précédemment, des échantillons de phase différentielle d0(t) et d0' (t). Chaque résultat de l'ad- dition calculée par le circuit additionneur 22 est mémorisé dans le premier registre de la pile, les contenus de chaque registre accumulateur de la mémoire 23 étant poussés successivement dans les registres suivants de la pile à chaque nouveau résultat fourni par le circuit additionneur 22, additions et décalages ayant lieu à la fréquence d'échantillonnage N/T. De la sorte, à la somme des L derniers coûts, obtenus pour chaque suite d'échantillons espacés de N, est retranchée le coût C.d0'(t) sortant de la mémoire 15 et ajouté le coût (t) entrant C.d0(t) dans la mémoire 15. La troisième partie du bloc de calcul de l'horloge de synchronisation comprend une mémoire tampon 24, un circuit comparateur de supériorité 25 et un registre tampon 26. Elle a pour fonction de repérer dans les N sommes des coûts, qui sont constamment mises à jour dans la mémoire 23 par le circuit additionneur 22, celle qui donne une amplitude maximale, ce repérage ayant lieu en temps réel avec le débit binaire de la modulation, c'est-A-dire à la fréquence symbole 1/T. Chaque somme de coût maximal est mémorisée dans la mémoire 24, o elle reste tant qu'une nouvelle somme n'est pas trouvée pour la rem- placer. Cette mémorisation est cormmandée par le circuit compara- teur de supériorité 25 qui compare A tout instant le contenu de la mémoire tampon 24 avec.chacune des sommes des coûts calcu- lés par le circuit additionneur 22. Le circuit comparateur de supériorité 25 fournit également un signal de commande au regis- tre tampon 26 pour mémoriser chaque nouveau numéro d'échantil- lon qui donne lieu à un nouveau maximum de la fonction coût. Tous les N échantillons, la mémoire 24 a en mémoire la somme des coûts la plus forte et le registre tampon 26 possède le numéro d'horloge associé à cette somme. Le numéro de l'hor- loge est alors transféré dans un registre tampon 27 relié au registre tampon 26 et le processus précédent se répète. La quatrième partie du bloc de calcul de l'horloge de synchronisation comprend une logique de commande 28 qui agit en fonction du numéro d'horloge mémorisé dans le registre tam- pon 27. Suivant le type de modulation utlllsé, la logique de commande fournit au décodeur 5 le n0 de l'horloge de décodage trouvé dans le registre tampon 27 et une information pour com- mander le bloc de séquencement et d'asservissement 6 du démodu- lateur. Le décodage de la phase différentielle est effectué par le décodeur 5, celui-cl décode l'échantillon de phase différen- tielle dont le numéro d'horloge est celui retenu par la fonction coût Pour permettre au démodulateur de fonctionner avec plu- sieurs types de modulation différents une mémoire programmée spécifique de chaque modulation peut être incluse dans le déco- deur 5 et à la suite de cette mémoire, plusieurs types de cir- cuits de sortie pourront être envisagés, pour adapter le fonc- * tionnement du démodulateur aux différents types de transmis- sions. Le procédé de démodulation et le démodulateur qui vien- nent d'être décrits ont des performances acceptables si la dé- rive en fréquence de la porteuse reste faible. La dérive en fréquence sur la porteuse radio est liée à la différence de fréquence entre les horloges pilotes des postes émetteurs et récepteurs. Cette dérive peut être par exemple 'de + lkHz pour une fréquence porteuse de 100 MHz avec une fréquence d'horloge pilcte de 5 NlHz à + 10. A cela doit être ajouté ou retranché la dérive propre à l'effet Doppler qui est estimé pour les véhicules terrestres à +10 7fois la fréquence de la porteuse radio ce qui correspond dans l'exemple précédent à une dérive de +10 Hz. Pour les liaisons avion-sol ou entre avion la valeur de la dérive due à l'effet Doppler est plus importante et peut atteindre +2 10- 6 fois la fréquence de la porteuse radio, ce qui dans l'exemple précédent donne 200 Hz de dérive. Cet effet de dérive se traduit par un glissement positif ou négatif de la phase différentielle. Avec 1o par exemple une fréquence symbole F symbole de 25 KHz et une dérive A f estimée à + 1KHz on obtient /Af/F symbole = + 0,04 c'est-à-dire un glissement de 15 de la phase différentielle. De ce fait, la présence de dérives provoquent un décalage des caractères symboliques propres à la modulation de phase utilisée. Pour des dérives importantes, le procédé et le démodulateur décrits précédemment doivent être adaptés pour tenir compte de ces dérives et du fait que l'estima- tion de la dérive et de l'horloge d'échantillonnage doivent s'effectuer en des temps très courts inférieures à la durée des paliers des systèmes de transmission en évasion de fréquence contenant les données -utiles. Ce problème peut être résolu en effectuant une estimation simultanée de l'horloge d'échantillon- nage de réception et de la dérive de la porteuse. En fait, la dérive sur la porteuse peut être estimée au moyen du décentrage qui est observé sur la phase différentielle. De ce fait cette dérive peut être estimée en degré par période symbole. Au lieu d'utiliser une fonction coût paramétrée par l'instant d'échantillonnage de réception, le procédé mis en oeuvre utilise une fonction coût paramétrée par un vecteur dont les deux composantes sont l'horloge d'échantillonnage et la dérive. Il est possible de représenter cette fonction coût, le paramètre horloge d'échantillonnage étant fixé, puisque cette fonction est invariante par translation de ce paramètre. Cette fonction peut être représentée dans un espace à trois dimen- sions, un premier plan étant formé par l'axe des dérives sur la porteuse et l'axe de la phase différentielle, l'axe des ordon- nées représentant le coût étant perpendiculaire à ce plan. La même fonction coût décrite précédemment au regard des figures 4 à 6 peut alors être utilisée. Cette fonction sert de référence aux calculs car elle correspond au cas o la dérive sur la por- teuse est nulle. La fonction de référence peut être considérée comme la projection d'une fonction coût générale sur le plan d'équation: dérive = O. Les projetés de la fonction coût générale sur les plans d'équation dérive = Cte sont obtenues par translation du paramètre phase différentielle de la quantité donnée par la dérive. Une représentation des projections est montrée aux figu- res 11 et 12. Le paramètre estimé est le paramètre qui maxirruse la statique ainsi obtenue. Avec les notations précédentes et l'adjonction d'un para- mètre dérive de, l'estimation statistique (2) devient S(kde)=C(d(tno de)+ C(d0(tno dk)e) e tno,k) e(tn no&l,k),e +c(Â") (3) +'.. +C(d{(tno*L, k) de) Pour compenser l'effet de cette dérive et en estimer sa valeur en degrés par symbole, le bloc de calcul de l'horloge de synchronisation est modifié de la manière représenté à la figure 13, o les éléments snimtlaires à ceux de la figure 10 sont repré- sentés avec les mêmes références, par i'adjonction d'un compteur de dérive 29, et l'introduction de tables de codage des coûts dans le dispositif de codage 21 adressés par le compteur 29 pour tenir compte du paramètre de dérive. Le bloc de calcul représen- té à la figure 13 effectue pour chaque valeur discrète de la dérive en fréquence présumée, la somme des L derniers coûts indexés par le même numéro de l'horloge d'échantillonnage. Ce calcul est effectué par l'additionneur 22 et la mé- moire 23 de façon similaire à celui exécuté par le bloc de cal- cul de la figure 10 excepté cependant en ce qui concerne le nombre de sommes partielles à comparer qui est égal à N fois le nombre nd de valeurs discrète de dérive prises en compte. Pour cela la mémoire 23 doit contenir autant de piles de regis- tres qu'il y a de nombre nd de valeurs discrètes de dérives. La mémoire tampon 24 mémorise comme précédemment la somme des coûts de valeur maximale, cette dernière étant indexée à la fois avec le numéro de l'horloge d'échantillonnage et la valeur présumé de la dérive. Sur le schéma de la figure 13, la valeur présumée de la dérive est transmise par le compteur de dérive 29 à l'intérieur du registre tampon 26. Le registre tampon 26 possède de la sorte à la fin de période symbole de référence la valeur présu- mée de la dérive et du numéro d'horloge qui ont donné la somme des L derniers coûts la plus forte, c'est-à-dire une parmi les Nnd sommes. Ces valeurs sont transmises cormmre sur la fi- gure 10 dans le registre tampon 27. Pour permettre le décodage des symboles le décodeur 5 de la figure 7 doit recevoir le o les éléments de phase différentielle affectés du bon numéro d'horloge de décalage fourni par la logique de commande 28 et la valeur de la dérive en fréquence également contenue dans le registre tampon 27 et fournie par la logique de commande 28. La valeur de la dérive en fréquence est appliquée à une entrée du décodeur 5 pour lui permettre d'effectuer une translation de la phase différentielle des symboles à décoder. REVENDICATIONS 1. Procédé de démodulation de signaux numérique d'ampli- tude d'enveloppe constante et modulés en phase et/ou en fré- quence de façon continue, les signaux numériques étant compo- sés par des suites de symboles de même durée T et juxtaposés les uns aux autres dans le temps, les symboles étant transmis au moyen de sauts de phase et/ou de fréquence, caractérisé en ce qu'il consiste: - à détecter (1) la phase du signal modulé sous la forme d'un signal d'amplitude variable fonction des variations de phase du signal modulant; - à échantillonner (7), suivant une période d'échantillon- nage égale à T/N, la phase du signal modulé détecté, à l'inté- rieur de chaque durée symbole T au moyen d'un peigne d'impul- sions d'échantillonnage obtenu à partir de N impulsions d'horlo- ges régulièrement espacées d'un Intervalle de temps égal à T/N - à calculer (1;8-14) l'angle de phase différentielle d0 séparant chaque échantillon de phase d'un symbole courant détecté, de l'échantillon de même rang situé dans le symbole détecté au temps symbole précédent; - à coder (21) chaque angle de phase différentielle d0 en lui attribuant un coût de codage dépendant de son amplitude; - à calculer (22, 23) les sommes des coûts de codage pour les échantillons de même rang appartenant à un nombre déterminé L de symboles consécutifs; - à identifier (24, 25), parmi les N sommes obtenues celle qui donne le résultat le plus grand en mémorisant le rang occupé dans chaque symbole par les échantillons correspondants; - à prendre (26) comme horloge de synchronisation symbo- les, pour définir les instants de synchronisation du démodula- teur, l'horloge qui dans le peigne d'impulsions possède le même rang que le rang de l'échantillon mémorisé; - et à décoder (5) l'angle de phase différentielle qui possède le même rang que l'échantillon mémorisé. 2. Procédé selon la revendication 1, caractérisé en ce que le coût de codage à une valeur maximum lorsque l'angle de phase différentielle calculé correspond aux instants de synchro- nisation du démodulateur et une valeur nulle au milieu de l'ln- tervalle compris entre les instants de synchronisation consécu- tifs. 3. Procédé selon la revendication 2, caractérisé en ce que la loi de variation du coût de codage entre une valeur nulle et une valeur maximum est linéaire. 4. Procédé selon la revendication 2, caractérisé en ce que la loi de variation du coût de codage entre une valeur nulle et une valeur maximum est représentée par une fonction polygo- niale. 5. Démodulateur pour la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'il comprend: un étage de mise en forme (7) pour détecter la phase du signal modulé sous la forme d'un signal d'amplitude variable fonction des variations de phase du signal modulant et échantillonner le signal d'amplitude variable obtenu suivant la période d'échantillonnage T/N, un bloc de calcul (1) de l'angle de phase différentielle d0 séparant chaque échantillon de phase d'un symbole courant détecté de l'échantillon de même rang situé dans le symbole détecté au temps symbole précédent, un bloc de mémorisation (3, 15) pour mémoriser pendant une durée de L symboles chaque échantillon de phase différentielle d0, un dispositif de codage (2)1 pour coder chaque angle de phase différentielle fourni par le bloc de calcul, et l'angle de phase différentielle de même rang mémorisé dans le bloc de mémo- risation (3, 15) calculé L symboles plutôt, en leur attribuant un coût de codage dépendant de leur amplitude, un circuit addi- tionneur-accumulateur (22, 23) pour soustraire entre eux les coûts calculés par le dispositif de codage et séparés dans le temps de la durée de L symboles, et ajouter au résultat de la soustraction pendant la durée de L symboles, les coûts de co- dage obtenus pour les échantillons de même rang à l'intérieur de chaque symbole, un circuit de comparaison (24, 25) pour identi- fier, parmi les N sommes obtenues celle qui donne le résultat le plus grand, un registre tampon (27) pour mémoriser le rang occupé dans chaque symbole par l'échantllon qui a contribué à l'obtention du résultat le plus grand, une logique de commande (28) pour sélectionner l'horloge de synchronisation dont l'impul- sion dans le peigne de raie possède le même rang que le rang de l'échantillon mémorisé et un décodeur (5) pour décoder l'angle de phase différentielle qui possède le même rang que l'échantil- Ion mémorisé dans le registre tampon (27). 6. Démodulateur selon la revendication 5, caractérisé en ce que le circuit de codage (21) comprend une mémoire morte programmable pour mémoriser sous la forme de tables les diffé- rents coûts à attribuer à chaque angle de phase différentielle calculé par le bloc de calcul. 7. Démodulateur selon la revendication 6, caractérisé en ce qu'il comprend un compteur de dérive (29) pour estimer la dérive de la fréquence de l'onde porteuse, pour adresser les tables de coût en fonction de la dérive estimée. 8. Démodulateur selon la revendication 6, caractérisé en ce que le compteur de dérive (29) transmet la valeur de la dé- rive estimée au décodeur (5) au travers du registre tampon (27) et de la logique de commande (28) pour effectuer une transla- tion égale à la valeur de la dérive de la phase différentielle des symboles à décoder. 9. Démodulateur selon les revendications 8 et 9, caracté- risé en ce que le clrcuit-additionneu;r-accumulateur (22, 23) permet l'exécution de N x nd sommes, nd étant égal au nombre discret de dérive fourni par le compteur de dérive (29).