Brevet FR2615675A1 — publié le 19881125

Procédé de démodulation de signaux modulés numériquement et dispositif de mise en oeuvre d'un tel procédé La présente invention concerne un procédé numérique de démodulation de signaux modulés numériquement par des données ou de la parole numérisée. Elle permet de restituer, après échantillonnage et numérisation des signaux, sans aucun organe analogique, les signaux modulants, la porteuse et le rythme. Une nombreuse littérature a été publiée dans le domaine de la démodulation des signaux et de nombreux produits étudiés et commercialisés existent sur le marché à des débits allant de quelques dizaine de bit/s (signaux télégraphiques) à quelques centaines de Mbit/s (liaison numérique sur câble transportant des multiplex numéri- ques de voies téléphoniques jusqu'à 140, voire 560 Mbit/s - liaison AMRT (accès multiple par répartition dans le temps) par satellites à 120 Mbit/s: Système Intelsat et Eutelsat). Ces différents produits exploitent tous des circuits adaptés à chacune des fonctions à réaliser (restitution de la porteuse -restitu- tion du rythme - circuit de décision après filtrage post-détection...) Ces modes de réalisation sous forme analogique ou transposés des concepts analogiques classiques utilisent tous le principe du traitement en temps réel qui consiste à ne mémoriser que les résultats du traitement. Seule, la démodulation différentielle assure une mise en mémoire du signal origine mais n'opère sur celui-ci qu'un retard fixe. Ce principe de traitement a des conséquences néfastes pour les performances et la complexité. - Sur les performances; car selon ce principe on ne s'adapte pas au mieux à la variabilité des paramètres, et de ce fait on est plus sensible aux signaux perturbateurs (bruit, raie parasite) superposés aux signaux à mesurer; - sur la complexité; car en utilisant ce principe, il est très difficile de s'adapter au mieux à chacune des étapes que l'on doit mettre en oeuvre pour aboutir au but recherché. En particulier les étapes d'acquisition et de poursuite nécessitent soit des dispositifs différents, soit des caractéristiques de circuit analogiques adaptés à chacune des étapes et il est très difficile, sans compliquer les dispositifs mis en oeuvre, d'exploiter et d'adapter au mieux le dispositif de mesure à chacune des étapes. - 2 - Concernant la démodulation, l'acquisition des caractéristiques de phase porteuse et rythme demande beaucoup plus de calculs que sa poursuite et la démodulation du signal. Une solution traditionnelle ne permet pas d'optimiser l'organe de traitement aux différentes situations. Les procédés traditionnels de démodulation utilisent tant pour la restitution de porteuse que pour la restitution de rythme, différents procédés d'asservissement de phase (boucle "phase lock" après élévation au carré - boucle de Costa...) la numérisation n'intervenant que dans le cadre de réalisation technologiques particulières, ou d'approximations linéaires réclamant beaucoup de calculs. Le procédé de l'invention utilise largement les mémoires numériques adressables et le logiciel, ce qui élimine toutes les diffi- cultés mentionnées ci-dessus et permet de disposer d'un algorithme global et d'approcher les performances optimales (filtre adapté) avec le minimum de calcul. L'invention a donc pour objet un procédé numérique de démodulation de signaux modulés numériquement, caractérisé en ce qu'il comprend suc- cessivement - une phase d'acquisition des signaux échantillonnés et numérisés, - une phase de calcul et de stockage des échantillons de phase, - une phase d'élaboration d'hypothèses de phase d'horloge, - une phase de calcul de la phase porteuse et de la phase rythme sur n symboles, - une phase de démodulation proprement dite par comparaison entre la phase du signal durant le symbole considéré et la phase de la porteuse. L'invention a également pour objet un dispositif de mise en oeuvre de ce procédé qui comprend un organe de traitement numérique, recevant en entrée le signal à analyser après traversée d'un numériseur auquel est relié un pilote stable. Plus particulièrement cet organe de traitement comprend: - un circuit de mémoire commune; - un circuit d'élaboration des échantillons de phase; un circuit d'élaboration des hypothèses de phase rythme; - un circuit de calcul de la phase porteuse (par exemple 5 symboles); - un circuit de choix de l'hypothèse phase rythme; - un circuit de calcul de la phase du signal (par exemple un symbole); - un circuit de démodulation; - un circuit de gestion d'interface avec le démodulateur; - un circuit interface des signaux démodulés; tous ces circuits étant reliés à un bus d'interconnection par des liaisons bidirectionnelles, les deux derniers circuits étant reliés entre eux par une liaison bidirectionnelle; - un circuit d'acquisition des échantillons numérisés, en entrée, relié au circuit de mémoire commune par une liaison mono-directionnelle; le circuit interface des signaux démodulés ayant deux sorties: - une sortie signaux démodulés, - une sortie qualité signaux, qui forment les sorties dudit organe de traitement. Avantageusement les différentes tâches de calcul nécessaires à la démodulation peuvent être implantée sur des organes de calcul pilotés par logiciel, chaque organe étant indépendant des autres organes et pouvant s'adresser à des mémoires pour y lire l'un quelconque des résultats des autres organes aussi bien que les échantillons de signaux inscrits à l'origine ainsi que ses propres résultats de calcul ou pour y inscrire dans une zone qui est affectée à lui seul ses propres résultats de calcul, ces mémoires étant soit communes soit constituées de réseaux maillés. L'invention peut s'appliquer notamment à un signal modulé de type "PSK" à un ou plusieurs niveaux, modulé de manière cohérente ou différentielle, à un signal dérivé tel que offset QPSK ou tout procédé de modulation permettant l'obtention d'un signal à amplitude constante (MSK - modulation cohérente de phase), à un signal modulé de type "PSK" avee des sauts d'amplitude, l'amplitude étant traitée par des procédés connus, à un signal modulé "FSK", le test d'hypothèses effectué tenant compte de plusieurs fréquences de modulation. Les caractéristiques et avantages de l'invention ressortiront d'ailleurs de la description qui va suivre, à titre d'exemple non limitatif, en référence aux figures annexées sur lesquelles: - les figures 1A et lB représentent un organigramme illustrant le 2 615675 procédé de l'invention, - les figures 2, 3 et 4 montrent des représentations schématiques illus- trant le procédé de l'invention, - la figure 5 illustre un dispositif mettant en oeuvre le procédé de l'invention. L'invention est relative à la démodulation de signaux reçus d'un satellite géostationnaire dans des conditions de Eb/No (énergie par bit/bruit par Hertz) critiques et de rythme d'horloge et de porteuse mal connus. Le procédé proposé permet d'approcher les performances d'un filtre adapté par l'utilisation d'estimateurs linéaires travaillant sur la phase du signal reçu. Le procédé peut être appliqué aux types de modulation utilisés: - sur les satellites de télécommunications, que les organes de démodulation soient situés au sol ou embarqués sur la charge utile (modulation "BPSK" (Binary Phase Shift Keying) ou "QPSK" (Quaternary Phase Shift Keying) en "ARMT" (Accès Multiple à Répar- tition dans le Temps) ou "AMRF" (Accès Multiple à Répartition en Fréquence) avec un ou plusieurs canaux par porteuse. - sur les réseaux terrestres, qu'il s'agisse de câbles ou de fibres optiques (signaux en bande de base) ou de faisceaux hertziens ou encore de signaux modulés dans une bande téléphonique; - sur les réseaux de distribution, chez l'abonné, pour des usagers à valeur ajoutée (en frontal d'une base de donnée, d'un calcula- teur ou d'un centre de messagerie) ou dans un centre P.T.T; - sur tout système utilisant des signaux modulés (télécommande- télémesure - signalisation téléphonique-mesure de distance...). D'autres retombées peuvent aussi être développées, en particulier, en ce qui concerne l'élimination de signaux connus d'un signal à démoduler (annulation d'écho), le désembrouillage de signaux de données modulés. Le procédé de l'invention est basé sur une algorithmie globale de démodulation élaborée à partir d'un test d'hypothèses sur la phase du rythme symbole du signal modulé. Afin de bien comprendre la succession des opérations effectuées, -5- on se base sur un signal reçu pour lequel le rythme symbole a une valeur connue mais présente une incertitude que l'on limitera à 1%, la phase de ce rythme n'étant pas connue. Par ailleurs, le signal reçu est modulé par saut de phase de telle sorte que durant une période du rythme symbole, le signal est constitué d'une porteuse sinusoïdale de fréquence connue et dont la phase reste constante. Le problème à résoudre revient donc à retrouver la phase du rythme (c'està-dire déterminer les instants de début et de fin de symbole) retrouver la phase de la porteuse en supprimant la modulation dont les caractéristiques sont supposées connues, puis en déduire la phase du signal durant chacune des périodes du rythme, celle-ci étant significative d'un symbole particulier afin de restituer la valeur de ce dernier. Le procédé de l'invention comprend successivement: - une phase d'acquisition des signaux échantillonnés et numérisés, - une phase de calcul et de stockage des échantillons de phase; la démodulation est basée sur l'acquisition des signaux échantillonnés et numérisés. A partir des échantillons ainsi stockés, on construit le signal analytique associé et la phase développée dudit signal. On dispose ainsi d'une représentation du signal reçu par des échantillons de phase auxquels on peut associer la valeur de l'amplitude. - une phase d'élaboration d'hypothèses de phase d'horloge; travaillant à partir des échantillons stockés, on élabore différentes hypothèses de phases rythme symbole. - une phase de calcul de la phase porteuse et de la phase rythme sur n symboles; partant de chacune des hypothèses précédentes, par des opéra- tions modulo -X et en utilisant le "barycentre cylindrique", qui sera décrit ultérieurement, on supprime la modulation afin d'évaluer la phase porteuse durant chacun des symboles, pour effectuer une estimation de la phase porteuse sur S symboles, S étant choisi en fonction de I'incerti- tude sur le rythme et sur la porteuse. La connaissance de la phase rythme et partant de la phase porteuse est obtenue en maximisant un paramètre H qui est fonction de la phase ainsi obtenue et des valeurs caractéristiques des transitions de phase entre symboles; - 6 7- -6 - une phase de démodulation proprement dite par comparaison entre la phase du signal durant le symbole considéré et la phase de la porteuse; la comparaison entre la phase du signal durant le symbole considéré, et la phase 0 de la porteuse permet en effet de démoduler le signal (trai- tement du signe de D) pour chacun des symboles. K On peut se reporter à l'organigramme de la figure 1 pour mieux situer la suite logique de ces différentes phases et notamment en ce qui concerne le calcul de: - de la phase porteuse fHy, chaque symbole allant de Sp à Sq, n étant le numéro d'ordre du symbole avec p - du maximum HM du paramètre H ce qui permet de déduire la phase porteu- se Orelative au Hy de l'hypothèse retenue et la phase rythme corres- pondante. Le signal à démoduler a les caractéristiques de la figure 2, sur laquelle on a représenté durant une période de rythme T un signal sinusoidal de phase fixe, le signal pouvant avoir une phase différente durant le symbole précédent et le symbole suivant: On peut aussi le décrire par: x(t) = A sin ( w t + + n(t) A: amplitude crête W: pulsation de la porteuse 0: phase durant le symbole Sn On a également représenté en haut de cette figure la porteuse de référence. La première phase du procédé de l'invention consiste en la trans- formation du signal en échantillons numérisés. Une construction du signal analytique suivie de la phase développée comme décrit dans la demande de brevet déposée ce jour sous le titre "Procédé d'évaluation numérique de la fréquence et de la phase de signaux et dispositif de mise -en oeuvre d'un tel procéd/'4 35. permet alors d'obtenir des échantillons représentatifs de la -7- phase du signal sous forme numérique sur lesquels sont basés tous les algorithmes. La figure 3 représente le signal pour une succession de symboles et les emplacements des échantillons de phase qui le représentent; avec des transitions de phase 10 en début et fin de symbole, et la position temporelle 11 des échantillons de phase du signal. La fréquence d'échantillonnage est choisie de telle sorte que deux échantillons successifs soient "Juste" décorrellés entre eux. Tenant compte du fait que la fréquence d'échantillonnage et le rythme ne sont pas cohérents, le choix de la fréquence d'échantillonnage doit tenir compte: - du rythme symbole et de l'incertitude correspondante; - de l'incertitude sur la fréquence porteuse. Elle doit donc permettre, tenant compte de ces éléments, de dispo- ser du nombre d'échantillons par symbole permettant le calcul de la phase porteuse sur chaque symbole avec la précision requise pour les perfor- mances de la démodulation. Le rythme étant connu, lors de la phase d'élaboration d'hypothèses de phase d'horloge il s'agit de placer les transitions de phase 10 à la -20 bonne position temporelle. Etant donné une période de durée égale à celle d'un symbole et placé de manière quelconque relativement aux échantil- lons de phase (T et te glissent l'un par rapport à l'autre) on partage cet intervalle de temps en n parties égales et on considère successi- vement la probabilité de trouver les transitions de phase de signal dans chacune des partitions de T. Tenant compte de ces hypothèses, la phase rythme choisie sera déterminée en fonction d'un critère défini et calculé selon les principes que l'on va décrire maintenant. On effectue les opération suivantes sur les échantillons de phases modulo lx: La représentation des angles de phase modulo lx est déduite de la repré- sentation des nombres modulo N, utilisé en numération (binaire, décimale, hexadécimale...) - 8 - Pour tout angle f compris entre -6o et +c*. sa représentation module-lx est donnée par un angle compris entre -x et +x tel que si est compris entre (2p-1) x et (2n+1) x, on a: X 1 = 2n-;n étant entier, et: -x 4 +x De la même façon on désignera une opération e= (C' -ó0) mod -x à partir de l'opération Y' = + et de la transformation modulo décrite ci- dessus qui revient à dire que: - si -x ç Y," = x on fait ' - ( y' -2x) et l'on renouvelle le test Jusqu'à l'obtention de -x ú D' Dans le cas particulier o x = (N entier), on effectue un recou- N vrement de tous les secteurs de dimension 2x compris entre.(2k - 1) x et (2k +1) x (k entier) sur le secteur -x Si x = r, la classe d'équivalence est de dimension 27. La notion de "barycentre cylindrique" est déduite de la notion classique de barycentre et adaptée à la configuration cylindrique selon laquelle on peut représenter les angles temporellement. Cette notion de barycentre cylindrique modulo lx repose sur le fait que le barycentre d'angles situés sur une plage d'amplitude inférieure à x se calcule par une pondération uniforme. Ainsi, pour calculer le barycentre de 1 phase X1 S 2..l1 définis sur l'intervalle. -x, + x, on calcul séparément les centres de gravité des phases positives ó + et des phase négatives - sachant que, compte tenu de la règle des modulo si 0 = (2k + 1) x, il lui correspond, en particu- lier si k = o, = -x En conséquence, si les phases sont évolutives, on passe directement de ó= + x à _= -x -9- On a aussi: TI t si 1+ > o b+= { > o si l+= o 1 > l' t O = { t o si 1 = o Ramenant l'origine en on obtient: o = ( + --) mod - x et prenant le barycentre cylindrique modulo x de + et avec leur pondération respective l+ et 1-, pour revenir ensuite à l'origine initiale, on obtient ainsi le barycentre modulo - x: _= (l + -) modulo + x On peut vérifier que lorsque ( + --) C x, le procédé de l'invention revient à appliquer la moyenne à l'ensemble des phases sur l'inter- valle -x, +x. Par contre, si ( + --) x, on doit travailler sur l'intervalle - 2x, o pour obtenir le résultat. La construction du barycentre cylindrique telle qu'elle est décrite permet de restituer la valeur moyenne de la phase en tenant compte des incertitudes irdiérente à l'évaluation de cette dernière (Evaluation de la phase sur des intervalles de dimension 2 7Fou sous multiple de ce dernier), ou d'effectuer des repliements de ces intervalles (suppression de la modulation) ou encore de prendre la bissectrice de plusieurs angles de phase. Les échantillons représentatifs de l'amplitude du signal peuvent avantageusement être introduits pour pondérer chacune des phases corres- pondantes et ainsi, tenir compte de l'incertitude liée au bruit. - 10 - La figure 4 illustre la construction du barycentre cylindrique; en effet on a: 5. Nombre de ei > o: 1+ (l, 02, L.... el) À Nombre de oi -1 + 1*-2' - 1 +=--(41+ 2+ + * +.... 1 (b3 + + + 1) 1-2 1-1 0' =0+ - 0. > o puis -o = mod x soit: +' ')o ui o o Si e0 x o = ' - 2 x Y' échantillons de phase qui sont inclus dans la période correspondante. En utilisant la méthode du barycentre cylindrique sur la durée du symbole, on détermine la phase 's du signal pour le symbole considéré. Le nombre d'échantillons pris en compte (rapport T/te) dépend des différentes incertitudes (porteuse et rythme). Dans lecas o porteuse et rythme sont parfaitement stables, un échantillon de phase par symbole peut être considéré comme suffisant. En éliminant la modulation en appliquant une règle modulo adaptée au type de modulation, et en utilisant le barycentre cylindrique pour le plus grand nombre de symboles compatibles avec les instabilités de por- - 11 - teuse et de rythme, on obtient la phase porteuse avec une ambiguité modulo liée au type de modulation, ambiguité qu'il est aisé de lever par continuité sur les estimations successives. On renouvelle toutes les opérations décrites ci-dessus pour chacune des hypothèses de phase. Disposant, pour chaque hypothèse sur une fenêtre composée de S symboles: - de la phase porteuse; - des valeurs caractéristiques destransitions de phase entre un symbole et le suivant s: phase avant la transition, s: phase après la transition, s On caleule le paramètre H sur l'ensemble des S symboles selon: H ôjsin ( (ps) si ()- I ()/ et l'on choisit l'hypothèse de phase qui correspond au maximum de ce paramètre. On en déduit à 1/n près la phase rythme et la phase porteuse correspon- dante. On dispose alors de tous les éléments permettant la démodulation. : phase relative au symbole considéré tf: phase de la porteuse Pour chacun des symboles, l'information est portée par la valeur de l'angle. ( & -) La décision est effectuée sur un critère dépendant de cette valeur et du type de modulation considéré. Les deux exemples d'application COiphase L et BPSK différentiel) permettront de préciser la méthode de démodulation. XOn peut appliquer le procédé de l'invention à la modulation bi-phase L associée à des sauts de phase (-,). Dans ce type de modulation, on fait correspondre à chaque symbole binaire, les phases de la porteuse indiquées dans le tableau I situé en fin de description. On peut donc dire que la valeur (O ou 1) du bit est représentée par le sens du saut de phase situé au milieu du symbole. - 12 - Dans le but de se protéger contre les incertitudes de fréquence horloge, on suréchantillonne le signal en utilisant 2n échantillons par symbole (n échantillons par 1/2 symbole). On suppose, par ailleurs, que le rythme d'échantillonnage n'est pas un multiple du rythme symbole. Pour estimer la phase porteuse, à partir des échantillons de phase obtenus par le calcul sur le signal analytique: t= Arc tg q/r comme décrit lors de l'explication de la construction de la phase développée dans la demande déposée ce Jour citée précédemmment, on obtient 2n échantillons de phase par symbole du signal que l'on se propose de démoduler. La fréquence de la porteuse est connue: - soit à partir d'un préambule dont la longueur est suffisante pour obtenir l'estimation recherchée (transmission sur la base de "salves" Message SARSAT - Système AMRT) - soit à partir d'une connaissance a priori de la porteuse à une approxi- mation près. Il est donc possible de retirer des échantillons de phase du signal, l'approximation connue de la rotation de phase modula _ 7r. Pour chacune des hypothèses de phase d'horloge, on considère une succession de S symboles (durée ST) le nombre S choisi étant tel que S Z T (L T: incertitude de rythme symbole) soit négligeable devant T et que la variance sur la phase estimée de la fréquence soit négligeable devant la période de la porteuse. ^ A2 = avec e = 6 Tous les raisonnements et calculs qui suivent sont valables pour chacune des hypothèses de phase rythme. Chaque symbole comprend 2n échantillons de phase dénommés: e1 e2 _3'.' 2n- Afin de tenir compte de la position de la transition de phase au milieu du symbole: - on calcule le barycentre cylindrique des échantillons - 13 - 1 '2''. gn modulo îl- - on calcule le barycentre cylindrique des échantillons yn+1' n+2-"'- 2n modulo +7 - on prend le barycentre cylindrique module + 7r/2 des deux barycentres précédents. Pour réduire l'écart-type de l'estimateur de la phase porteuse, on prend le barycentre cylindrique des résultats obtenus sur les S symboles et l'on aboutit ainsi à un écart de phase avec la porteuse estimé à 7r/2 dont la variance est donnée par: = 1 évolué de moins de T;/2 entre deux estimations successives. Connaissant la phase t modulo + r de l'estimation précédente et la phase 't qui vient d'être calculée modulo - 1r/2 on calcule: Ot -t et: si - /2 I/2, O'test remplacé par: ( 't -1) modulo f On dispose ainsi de la phase de référence de la porteuse, l'écart ( 't - ft) pouvant être utilisé pour corriger l'estimation connue de la fréquence. Pour estimer la phase rythme symbole, On chercher à déterminer les échantillons de phase au voisinage desquels se trouvent les débuts et fin de chaque symbole. On dispose pour cela de 4n hypothèses de phase: 2n en phase avec les échantillons, 2n entre les échantillons. Pour chacune des hypothèses, on a calculé sur une séquence de S symboles, au milieu de cette séquence, la phase de la porteuse ( et sur chacun des symboles de la séquence, le barycentre cylindrique Vs - 14 - modulo - des n premiers échantillons et le barycentre cylindrique s modulo I1rdes n derniers échantillons. Si les échantillons du centre risquent d'apparaître simultanément avec la transition de phase, on cal- cule t sur 1 àn r et 's sur n+r+l 2n' r tant r eprsen- tatif de la durée de la transition de phase. C'est l'amplitude du saut de phase par rapport à la porteuse qui sert de critère pour le choix de la bonne hypothèse. L'hypothèse retenue est donc celle qui donne la plus grande valeur à la quantité H donnée par: H= sin (sin (-) -( ') L'opération décrite ci-dessus est relative à l'acquisition du rythme. Le rythme ainsi acquis peut servir à démoduler K bits situés au centre de la plage d'estimation de S bits. Pour les estimations suivantes, le nombre d'hypothèses peut être réduit aux seules hypothèses situées de part et d'autre du résultat de l'estimation précédente. Enfin, durant les périodes de synchronisation du message pour les- quelles la suite binaire est connue, un test plus efficace peut être obtenu si l'on se positionne à moins d'un symbole. À Pour réaliser la démodulation> on considère que la phase porteuse et la phase rythme symbole, estimées sur S symboles, servent à démoduler k ( La valeur du bit pour les k symboles considérés est alors donnée par le signe de l'expression. Sin (V -s.) - sin ( VI-S à) Cette expression peut aussi être directement utilis4e si l'on superpose N messages contenant la même information en vue d'une amélioration du taux d'erreur (rapport C/N faible) en effectuant un décodage pondéré ou une superposition symbole à symbole. -*On peut également appliquer le procédé de l'invention à la modulation par saut de phase. - 15 - On décrit ci-dessous le cas d'application le plus simple: modula- tion BPSK à + ff/2 selon le schéma de modulation du tableau II situé en fin de description, et l'on généralisera ensuite aux modulations de phase à M niveaux (2k-1) 7Tr/M (k entier._ (M-1)). Dans la description qui suit, on suppose que le rythme d'échantil- lonnage n'est pas un multiple du rythme symbole. Un échantillon de phase durant la partie stable de la phase d'horloge "juste" décoréllé pour le bruit de l'échantillon qui le précède, suffit donc à définir la phase durant le symbole. Les procédés permettant l'élaboration de cet échantillon pour les différentes hypothèses de rythme symbole considérées, sont basées sur des procédés connus de numérisation, de construction du signal analyti- que, de filtrage numérique, et de construction de la phase développée qui ont été décrits dans la demande de brevet déposée ce jour et citée précédemment. Le choix de l'échantillon de phase suppose, bien évidemment, que les circuits en question soient adaptés, en particulier au sens du filtrage des signaux, pour obtenir le résultat recherché. Tous les calculs qui suivent sont effectués pour chacune des hypothèses de phase symbole considérées de manière à parvenir, en effec- tuant un test d'hypothèses, à la définition de la phase symbole, puis à sa poursuite et à la démodulation du signal. Pour estimer la phase porteuse, disposant d'un échantillon de phase t s par symbole, on retire la rotation de phase correspondant à laconnaissance de la fréquence porteuse. On obtient ainsi une suite d'échantillons de phase et modulo + jrque l'on transforme en échantil- lon Et par une opération modulo - 7r/2 pour enlever la modulation comme indiqué ci-dessous: i -n/2 c Ot si et > /2 et + ot 'r On calcule ensuite le barycentre cylindrique modulo MT/2 sur le plus grand nombre de symboles S compatibles avec la stabilité de la fréquence 2615675- _ 16- porteuse et la connaissance de la phase symbole selon la méthode décrite précédemment. Pour ce faire, on calcule séparément les barycentres cylindriques des phases positives + et négatives - en écrivant: si 1+ > o t> o t += { 0 si 1+ = o i t si 1- > o t +t = = o si-1--=o puis o =( - '') mod + 1/2 i+ et enfin - =(- eo + e') mod + I/2 On aboutit ainsi à une valeur de phase à laquelle on ajoute r'/2 mod + /2 - 17 - On dispose ainsi de la phase porteuse pour chacune des hypothèses de rythme symbole considérée, l'écart quadratique de l'esti- mation étant 1 1 avec S: nombre de symboles considérés : rapport S/B incluant les bruits de phase et les instabilités de la porteuse. 10. Pour estimer la phase rythme, pour tester les hypothèses de phase symbole, on recherche des couples d'échantillons en bande de base enca- drant les débuts (ou fin) de symbole: soit(9s et 'sces valeurs pour le i symbole. L'hypothèse à retenir est celle qui donne a H la plus grande valeur calculée selon la formule. H - 1 sin ( - si) -sin ( ' -)| Ce test étant effectué, on retient a) la phase horloge reliée au maximum de H b) la phase porteuse correspondante. La poursuite de la phase horloge est ensuite effectuée sur les trois hypothèses encadrant la phase précédente. La démodulation est effectuée directement sur la valeur de la phase résultant de l'opération: t s - î selon le tableau III situé en fin de description, avec: es: phase du signal calculée sur chacun des symboies modulo _; : phase de la porteuse résultant de l'estimation décrite précédem- ment; (Ps etf3résultant de l'hypothèse de phase rythme retenue. *En généralisant on peut appliquer le procédé de l'invention à des modulations PSK à M niveaux (M > 2) et FSK. Toutes les opérations décrites précédemment restent valables pour une modulation PSK à M niveaux à condition: - d'effectuer le calcul de la phase porteuse en l'initialisant par une - 18 - opération modulo + 'I/M (au lieu de - 7r/2) puis en calculant le barycentre cyclindrique modulo + /M (au lieu de - MrI2) - d'effectuer la démodulation en utilisant la table des phases adaptée au type de modulation traitée. Dans le cas d'une modulation FSK, il y a lieu de considérer plusieurs fréquences. L'application du procédé revient donc à effectuer un test d'hypothèses, dans lesquelles plusieurs fréquences sont attribuées à chacun des symboles considérés. * Un dispositif de mise en oeuvre du procédé de l'invention comprend un organe de traitement numérique 20, recevant en entrée le signal à analyser SA après traversée d'un numériseur 21 auquel est relié un pilote stable 22. Cet organe de traitement 20 comprend: - un circuit 23 de mémoire commune; - un circuit 24 d'élaboration des échantillons de phase; - un circuit 25 d'élaboration des hypothèses de phase rythme; - un circuit 26 de calcul de la phase porteuse (par exemple 5 symboles); - un circuit 27 de choix de l'hypothèse phase rythme; - un circuit 28 de calcul de la phase du signal (par exemple un symbole); - un circuit 29 de démodulation; - un circuit 30 de gestion d'interface avec le démodulateur; - un circuit 31 interface des signaux démodulés; tous ces circuits étant reliés à un bus d'interconnection 32 par des liaisons bidirectionnelles, les deux derniers circuits étant reliés entre eux par une liaison bidirectionnelle; - un circuit 33 d'acquisition des échantillons numérisés, en entrée, relié au circuit 23 de mémoire commune par une liaison mono-direction- nelle. Le circuit 31 interface des signaux démodulés a deux sorties: - une sortie signaux démodulés SS1, - une sortie qualité signaux SS2, qui forment les sorties dudit organe de traitement 20. Un tel dispositif est utilisable quel que soit re type de démodula- - 19 - tion choisie. On trouve, dans cette architecture, tous les circuits nécessaires à l'élaboration du procédé de l'invention, certains d'entre eux étant adaptés à un type de modulation particulier. Dans le cas o la fréquence porteuse est connue avec une incerti- tude nécessitant l'acquisition et la poursuite de la valeur de la fréquence de la porteuse afin d'obtenir les performances de démodula- tion, on se basera, pour l'évaluation de fréquence, sur le procédé décrit dans la demande de brevet déposée ce Jour, citée précédemment. L'architecture décrite suppose: - une connaissance à priori de la fréquence porteuse compatible avec le nombre S de symboles pris en compte pour l'estimation de la phase porteuse. - une connaissance à priori du rythme symbole avec une incertitude inférieure au 1/10 du décalage entre deux hypotheses sur une durée égale à S symboles. Cette architecture permet: - l'acquisition et le stockage en mémoire du signal à démoduler sous forme d'échantillons numériques bruts issus d'un numériseur délivrant au rythme d'échantillonnage, l'amplitude du signal à démoduler; - la conversion des échantillons d'amplitude en échantillons de phase comprenant le traitement par filtrage des échantillons, la conversion du signal en ses deux composantes (réelle et quadrature) pour obtenir la phase échantillonnée; - l'élaboration d'hypothèses de phases rythme permettant undécoupage temporel du signal en tranches de durée T (1/T = Rythme) adaptées aux différentes étapes de démodulation (acquisition et poursuite). Dans le cas de la poursuite, le nombre d'hypothèses pris en compte sera réduit au minimum nécessaire. - le calcul pour chaque hypothèse de phase rythme, de la phase porteuse sur un symbole, puis éliminant la modulation l'estimation de la phase porteuse sur S symboles; - la comparaison entre phase symbole et phase porteuse ce qui permet l'extraction du signal démodulé (signification symbole supportée par la phase signal) ainsi qu'une appréciation de la qualité du signal modulé - 20 - durant chaque symbole et de la qualité de la démodulation. L'architecture matérielle utilisée pour construire un démodulateur appliquant le procédé de l'invention répond aux principes de: - Partage des tâches entre différents organes de calcul travaillant indépendamment les uns des autres; - Accès indépendant en lecture de l'un quelconque des organes de calcul, à tous les résultats des tâches effectuées par les autres organes. L'uti- lisation de mémoires adressables accessibles par- bus commun ou "multiport" ou organisées en réseaux maillés interconnectant les dif- férents organes est incluse; - - Inscription exclusive dans chacune des zones délimitées de la (ou des) mémoires par un organe de calcul préalablement désigné. Il est bien entendu que la présente invention n'a été décrite et représentée qu'à titre d'exemple préférentiel et que l'on pourra remplacer ses éléments constitutifs par des éléments équivalents sans, pour autant, sortir du cadre de l'invention. La description s'appuie principalement sur deux types de modula- tion (bi-phase L. associé à une modulation de phase à +; BPSK diffé- rentiel). Mais le procédé de l'invention peut être appliqué à tout type de modulation, à nombre fini d'états. Les différentes tâches de calcul nécessaires à la démodulation peuvent être implantée sur des organes de calcul pilotés par logiciel, chaque organe étant indépendant des autres organes et pouvant s'adresser à des mémoires pour y lire l'un quelconque des résultats des autres organes aussi bien que les échantillons de signaux inscrits à l'origine ainsi que ses propres résultats de calcul ou pour y inscrire dans une zone qui est affectée à lui seul ses propres résultats de calcul, ces mémoires étant soit communes soit constituées de réseaux maillés. -21 - TABLEAU I ] II J ler 1/2 Symbole 2ème 1/2 Symbole i ' I I I o I -_ +f f I 1 | l + I i I I TABLEAU II Phase relativement a la porteuse | I IF I 0 + r2 ___ I I -n/2 | 1 - nl2- TABLEAU III Valeur de | I (ts - O) 1 Valeur du bit | modulo n I 1>0 I I J >o O I I _ _ _ _j_ _ _ _! 26 15675 - 22 - REVENDICATIONS 1/ Procédé numérique de démodulation de signaux modulés numériquement caractérisé en ce qu'il comprend successivement: - une phase d'acquisition des signaux échantillonnés et numérisés* - une phase de caleul et de stockage des échantillons de phase, - une phase d'élaboration d'hypothèses de phase d'horloge, - une phase de calcul de la phase porteuse et de la phase rythme sur n symboles, - une phase de démodulation proprement dite par comparaison entre la phase du signal durant le symbole considéré et la phase de la porteuse. 2/ Procédé selon la revendication 1, caractérisé en ce que la démodula- tion est basée sur l'acquisition des signaux échantillonnés et numérisés, en ce qu'on construit le signal analytique associé et la phase développée dudit signal à partir des échantillons ainsi stockés, et en ce qu'on dispose ainsi d'une représentation du signal reçu par des échantillons de phase auxquels on peut associer la valeur de l'amplitude. 3/ Procédé selon l'une quelconque des revendications précédentes, carac- térisé en ce que en partant des hypothèses de phases rythme symbole par des opérations modulo X et en utilisant le barycentre cylindrique, on supprime la modulation afin d'évaluer la phase porteuse durant chacun des symboles, pour effectuer une estimation de la phase porteuse sur S symboles, S étant choisi en fonction de l'incertitude sur le rythme et sur la porteuse, en ee que la connaissance de la phase rythme et partant de la phase porteuse est obtenue en maximisant un paramètre H qui est fonction de la phase ainsi obtenue et des valeurs caractéristiques des transitions de phase entre symboles. - 4/' Procédé selon l'une quelconque des revendications précédentes, carac- térisé en ce que les échantillons de départ sont stockés en mémoire d'o ils sont extraits autant de fois que nécessaire aux fins de calculs effectués selon un algorithme unique appliqué à plusieurs hypothèses de phase rythme, d'o découle, selon un critère adapté au type de modulation de phase rythme la phase porteuse ce qui permet d'extraire simplement les informations démodulées et d'apprécier la qualité du signal modulé et de l'opération de démodulation. - 23 - / Dispositif de mise en oeuvre du procédé de l'invention selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend un organe de traitement numérique (20), recevant en entrée le signal à analyser (SA) après traversée d'un numériseur (21) auquel est relié un pilote stable (22). 6/ Dispositif selon la revendication 5, caractérisé en ce que l'organe de traitement comprend: - un circuit (23) de mémoire commune; - un circuit (24) d'élaboration des échantillons de phase; - un circuit (25) d'élaboration des hypothèses de phase rythme; un circuit (26) de calcul de la phase porteuse; - un circuit (27) de choix de l'hypothèse phase rythme; - un circuit (28) de calcul de la phase du signal; - un circuit (29) de démodulation; - un circuit (30) de gestion d'interface avec le démodulateur; - un circuit (31) interface des signaux démodulés; tous ces circuits étant reliés à un bus d'interconnection (32) par des liaisons bidirectionnelles, les deux derniers circuits étant reliés entre eux par une liaison bidirectionnelle; - un circuit (33) d'acquisition des échantillons numérisés, en entrée, relié au circuit (23) de mémoire communepar une liaison mono-direction- nelle; --le circuit interface des signaux démodulés ayant deux sorties. - une sortie signaux démodulés (SS1) - une sortie qualité signaux (SS2) qui forment les sorties dudit organe de traitement (20). 7/ Dispositif de mise en oeuvre du procédé de l'invention selon l'une quelconque des revendications 1 à 4, caractérisé en ce que les différentes tâches de calcul nécessaires à la démodulation peuvent être implantée sur des organes de caleul pilotés par logiciel, chaque organe étant indépendant des autres organes et pouvant s'adresser à des mémoires pour y lire l'un quelconque des résultats des autres organes aussi bien que les échantillons de signaux inscrits à l'origine ainsi que ses propres résultats de calcul ou pour y inscrire dans une zone qui est affectée à lui seul ses propres résultats de calcul. - 2.4 - 8/ Application du procédé selon l'une quelconque des revendications 1 à 4, à un signal modulé de type "PSK" à un ou plusieurs niveaux, modulé de manière cohérente ou différentielle, à un signal dérivé tel que offset QPSK ou tout procédé de modulation permettant l'obtention d'un signal à amplitude constante (MSK - modulation cohérente de phase). 9/ Application du procédé selon l'une quelconque des revendications 1 à 4, à un signal modulé de type "PSK" avec des sauts d'amplitude, l'amplitude étant traitée par des procédés connus. / Application du procédé selon l'une quelconque des revendications 1 à 4, à un signal modulé "FSK", le test d'hypothèses. effectué tenant compte de plusieurs fréquences de modulation.