FR 2711028 A1 19950414 FR 9312027 A 19931008 Procédé et dispositif de démodulation cohérente par blocs de signaux complexes modulés en phase La présente invention concerne un procédé et un dispositif de démodulation cohérente par blocs de signaux complexes modulés en phase. Elle s'applique aux transmissions de données numériques par voie radioélectrique avec faible rapport signal à bruit et décalage en fréquence variable dans le temps et notamment aux liaisons radioélectriques entre stations au sol via un satellite défilant. Dans les systèmes de démodulation cohérente de modulations à 2 ou 4 états de phase, il est connu, pour améliorer les performances de transmission de messages de taille minimale en présence d'un fort niveau de bruit et d'un fort décalage en fréquence non constant (rampe Doppler), d'utiliser un codage de type convolutif consistant à émettre plus de symboles que ce qu'il est strictement nécessaire. Les codes utilisés le plus souvent sont des codes de redondance R = 2 ce qui conduit à un débit utile égal à la moitié du débit transmis, un bit d'information se traduisant par 2 bits trans- mis. En théorie, le rapport signal/bruit énergétique par bit utile nécessaire pour obtenir un taux d'erreur de 10-5 est de l'ordre de 6 dB avec une marge de réalisation de 2 dB. Dans ces conditions, le rapport signal/bruit en entrée d'une démodulation de réception est de 6 dB log10R soit 3 dB, ce qui est très faible. Cela impose, dans les systèmes classiques de faire précéder tout message d'une séquence de synchronisation initiale suffisamment Ion- gue pour que, même dans des conditions de propagation très défavorables, la démodulation se fasse dans de bonnes conditions. Parmi les types de démodulation connus deux apparaissent utilisa- bles, la démodulation dite "au fil de l'eau" et la démodulation par bloc. Dans la démodulation "au fil de l'eau" des échantillons du signal sont fournis en continu à l'entrée du système de démodulation. Le système de démodulation comporte classiquement, d'une part, une boucle à verrouillage de phase, destinée à estimer en continu la phase qu'aurait le signal en l'absence de modulation, pour donner une phase de référence, et d'autre part, un démodulateur, fournissant en sortie soit direc- tement les valeurs de l'information décodée soit le signal complexe redressé. Dans ce système, I'organe le plus important est la boucle à ver- rouillage de phase, dont les caractéristiques sont le résultat d'un compromis entre la vitesse d'estimation, dont dépend le temps d'accrochage, et la "gigue" qui conditionne le taux d'erreur résiduel, ou dans le meilleur des cas, la dégradation des performances par rapport au cas o la phase du signal serait parfaitement connue. Dans la démodulation par blocs, des échantillons complexes du signal sont d'abord stockés dans une mémoire puis, quand ils sont en nom- bre suffisant, ceux-ci sont traités en bloc en une ou plusieurs passes. La construction des démodulateurs par blocs est récente car elle a été rendue possible par les progrès technologiques en matière de mémoires et de microprocesseurs de traitement de signal, seuls capables d'effectuer une démodulation complète sur les échantillons stockés en mémoire. Dans les systèmes à démodulation par blocs la démodulation a lieu en deux opéra- tions successives, une première opération qui consiste à estimer la phase de référence du signal par mise en oeuvre d'algorithmes divers travaillant sur la totalité des échantillons stockés en mémoire et une deuxième opéra- tion consistant à corriger la phase des échantillons et décider la suite du traitement. La suite des traitements consiste par exemple à effectuer un décodage au maximum de vraisemblance, le plus souvent en faisant usage de l'algorithme de Viterbi, lorsque le signal comporte un code convolutif. La même structure de matériel est utilisée pour effectuer à la fois la démodula- tion et le décodage. Une description d'un tel système, non compris le décodage correc- teur, et la comparaison de ses performances avec celles d'un système plus classique sont décrits dans l'article intitulé "Block Demodulation for Short Radio Packets" de J. Namiki, Electronics and Communications in Japan, J67-B, N I1 (Janv. 1984) pp.54-61. Comme indiqué dans l'article "Timing Extraction and Carrier Estima- tion for PSK Signal Block Demodulation System" de Y. Higashida, M. Hagiwara, N. Nakagawa, Electronics and Communications in Japan, J71- B, N 4 (Avr. 1988) pp.540-546, les systèmes de démodulation précédem- ment décrits ont d'énormes difficultés pour fonctionner à de très faible rap- port signal/bruit et une évaluation des performances de ces systèmes montre clairement que le rapport signal/bruit minimal de fonctionnement est plus proche de 10 dB que de 3 dB. Qui plus est, les cas o le décalage Doppler varie au cours du temps ne sont pas traités. Le but de l'invention est de pallier les inconvénients précités. A cet effet, I'invention a pour objet un procédé de démodulation cohérente par bloc de suites de signaux complexes Sn modulés en phase suivant des états de phase déterminés caractérisé en ce qu'il consiste à extraire la phase de chacun des signaux complexes en supprimant les sauts de phase inconnus, à rechercher dans la suite des phases obtenues une première plage de phase de longueur déterminée L pour laquelle la variation de phase est la plus linéaire possible, à déterminer une deuxième plage de phases de longueur L2 supérieure à L correspondant à une approximation linéaire plus précise, à effectuer itérativement des régressions linéaires à l'aide d'une fenêtre glissante de longueur égale à la deuxième plage de phase L2 pour supprimer les sauts de phase résiduels pour effectuer une régression quadratique sur les phases restantes et pour calculer une esti- mation globale des phases de référence des signaux de chaque bloc, et à effectuer une comparaison de l'erreur quadratique moyenne entre les pha- ses de référence estimées et les phases brutes obtenues par régression linéaire. L'invention a pour principal avantage qu'elle permet, la recherche de la phase de référence avec une faible charge de calcul. Elle permet égale- ment de rendre minimum la dégradation des performances pour les très fai- bles rapports signal/bruit, en permettant une qualification de la qualité du signal reçu et d'affiner la synchronisation. D'autres caractéristiques et avantages de l'invention apparaîtront ci- après à l'aide de la description et des dessins annexés qui représentent: - La figure 1 un dispositif de modulation selon l'invention. - La figure 2 les étapes principales du procédé selon l'invention mises sous la forme d'un organigramme. - Les figures 3A et 3B des diagrammes comparatifs montrant les effets de la correction de phase selon l'invention. - Les figures 4A et 4B des diagrammes de phases corrigées et d'erreur minimum pour un rapport signal/bruit de 3 dB. - La figure 5 un exemple de droite de régression linéaire tracée sur une plage étendue relativement à une plage initiale. - La figure 6 une illustration du processus d'estimation finale des phases de référence par prolongement à droite et à gauche de la droite de régression. - La figure 7 une représentation de phases corrigées initiales relati- vement à une régression quadratique des phases corrigées pour un rapport S/B = 3dB. - La figure 8 une représentation de phases corrigées relativement à une régression quadratique des phases corrigées pour un rapport S/B de 0 dB. - Les figures 9A, 9B, 9C, 9D des diagrammes de démodulation de signaux complexes à deux états de phase pour des rapports S/B de 3, 0, 6 et 10 dB. Le dispositif de démodulation selon l'invention qui est décrit ci-après a une structure démodulateur par blocs. Il comporte un dispositif 1 de pas- sage du signal S à démoduler en bande de base, un dispositif de mémorisa- tion 2, un dispositif d'estimation de phase de référence 3 et un dispositif de démodulation et de décodage 4. Le dispositif 1 extrait de façon connue sur une fréquence intermédiaire Fl les composantes réelles et imaginaires du signal complexe S à démoduler. Ce signal est formé par des suites de N échantillons Sn complexe de la forme: (nT)2 Sn=aexp 2 fn+2dAf expj(pn) + Zn, n=1...N o a = amplitude 00 = phase initiale Af = écart en fréquence initial dAf = variation de l'écart en fréquence (pn= phase porteuse de la modulation (M=2 ou 4 possibilités) T = durée du symbole et intervalle d'échantillonnage Zn = bruit gaussien complexe de puissance 2a2 a2 rapport signal/ bruit SI B = -2 2oy2 Le dispositif de mémorisation 2 est couplé en sortie du dispositif 1 et stocke "au vol" chaque bloc de N échantillons avant tout traitement. Le dispositif 3 estime la phase de référence du signal qui ne dépend pas de l'information transmise, à partir des échantillons stockés dans le dis- positif 2. Le dispositif 4 est couplé d'une part, en sortie du dispositif de mé- morisation 2 et d'autre part, en sortie du dispositif 3 d'estimation de phase pour reconstituer des échantillons complexes redressés et reconstituer l'in- formation utile transmise à partir des échantillons complexes stockés dans la mémoire 2 et des phases de référence fournies par le dispositif 3. Le dis- positif 3 d'estimation de phase opère suivant les étapes 5 à 10 du procédé représenté à la figure 2. L'étape 5 consiste à calculer la phase du signal complexe S appliqué à l'entrée du dispositif 1 et corrige les phases en sup- primant des sauts de phase inconnus dus à l'information qui le module. L'étape 6 consiste à rechercher dans la série des phases délivrées à l'étape 5 un plage de phase de longueur L o la variation de phase paraît être la plus linéaire possible. A l'étape 7 une plage de phases de largeur L2 supérieure à L est recherchée correspondant à une approximation linéaire plus précise. L'étape 8 effectue itérativement des régressions linéaires à l'aide d'une fenêtre glissante de longueur L2 pour supprimer tous les sauts de phase résiduels. L'étape 9 effectue une régression quadratique sur les phases obte- nues à l'étape 8 pour calculer une estimation globale des phases de réfé- rence sur chaque bloc de N échantillons traités. Enfin l'étape 10 effectue une comparaison de l'erreur quadratique moyenne entre les phases de réfé- rence estimées à l'étape 9 et les phases brutes obtenues à l'étape 8, à une valeur de seuil déterminée. La correction des phases qui est effectuée à l'étape 5 vise à la fois à supprimer l'influence de la modulation, qui se traduit par une incertitude de r sur la phase de référence du signal et à supprimer l'effet des évaluations de phase qui ont toujours lieu à 2n près. Elle s'effectue elle-même en trois éta- pes de traitement. La première étape consiste à corriger les phases des échantillons Sn en les ramenant par exemple dans l'intervalle [-2,+21 pour une modulation biphase 0 - 7C. Dans ce cas, si Pn = arg(Sn) est la phase comprise entre -7r et +7r du ne échantillon de signal, la correction a lieu en transformant Pn pour n = 1... à N par les itérations suivantes: faire Qn = Pn - si Qn > alors faire Qn = Qn - 7 - si Qn La deuxième étape consiste à calculer les différences de phase suc- cessives et à les ramener elles aussi entre -- et ±. 2 2 L'algorithme correspondant consiste à poser Dn = Qn - Qn-1 et pour n = 2... N à faire Dn = Dn - Ic si Dn 2 2 La troisième étape vise à reconstituer une suite de valeurs de pha- ses On représentative (à une incertitude de 7r près) de la phase de référence du signal. L'algorithme correspondant consiste à poser: (1 = Q1 et à faire On = On- 1 + Qn pour n = 2... N. Des valeurs respectives de Pn et On, pour un bon rapport signal/bruit de 20 dB et pour un mauvais rapport signal/bruit de 3 dB sont tracées respectivement sur les graphes des figures 3A-3B o les phases originales en radians de N = 100 échantillons sont tracées en traits pointillés et les phases corrigées apparaissent en trait plein. Comme cela apparaît sur la figure 3B, les phases apparaissent beaucoup moins bien corrigées lorsque le rapport signal/bruit est de 3 dB, car les phases corrigées qui devraient en principe être quasi continues, présentent des sauts de phase résiduels de 7r radians. Cependant, la figure 3B permet d'observer que les phases corrigées qui présentent des sauts de phase de r comportent aussi des plages sans sauts. De plus, si le décalage en fréquence ne varie pas trop vite, la phase varie (presque) linéairement sur une telle plage. L'intérêt de cette constatation a de l'importance, car elle permet lorsqu'une plage sans sauts est repérée, de corriger les sauts de phase en la prolongeant de part et d'autre. Cette particularité est mise à profit par l'invention à l'étape 6 en recherchant dans la suite des phases corrigées, une plage de L valeurs successives approximativement réparties sur une droite et pouvant être approximées par une équation de la forme: (Dno+i t A i + B pour i = 1,..., L no étant le numéro de l'échantillon précédent immédiatement le début du bloc. Par régression linéaire le minimum de l'erreur quadratique moyenne E donnée par E(n0) = L ('no+i - (A i + B)) i=1 est calculé. Les valeurs A et B correspondant à ce minimum sont données par les relations: A = P1 So - P0 S1 A B = PO0 S2- P1 S1 A avec ae L L(L + 1) s L 2 L(L + 1)(2 L + 1) S0=L,S1=Yi= 2__.2 i=1 i=1 2 A = SO S2 - S1 L L P0 P = no+i P =i (:no+i i=1 i=1 La valeur minimale de l'erreur est alors: Emin(nOi) = _L (no +j A Pl B PO i=1 n Il suffit alors de calculer Emin(no) pour no = O...N-L et de repérer la valeur de no correspondant au minimum. La valeur de no désigne la position du début de la plage o L échantillons successifs de phase corrigée sont le mieux approchés par une droite. Pour moins dépenser en puissance de cal- cul les différentes sommes faisant intervenir les phases (Dn peuvent aussi être évaluées de proche en proche. La valeur de L doit alors être soigneu- sement choisie d'une part, suffisamment faible pour que, même avec un fai- ble rapport signal/bruit, il soit quasiment certain de trouver une plage de phases corrigées quasi linéaires, sans saut de phase et d'autre part, suffi- samment élevée pour que la régression linéaire donne un résultat suffisam- ment fiable. Pour un rapport signal/bruit minimum de 3 dB, une valeur de L de l'ordre d'une dizaine d'échantillons donne des résultats acceptables, comme le montrent les graphes de phase corrigée en radian et d'erreur qua- dratique moyenne des figures 4A et 4B o il apparaît que sur N = 100 échantillons, I'erreur est minimale avec un point de départ no (offset) de 56 échantillons à partir de l'origine. Dans ces conditions c'est-à- partir de cette portion du signal qu'est cherchée l'estimation de la phase de référence. La largeur de la plage L est définie pour avoir une probabilité suffi- sante de trouver une plage correcte. Cependant, comme les valeurs des coefficients de régression linéaire A et B correspondant au minimum d'erreur sont entachés eux aussi d'erreurs, ceux-ci sont estimés avec plus de préci- sion en cherchant à élargir la plage L sur une longueur L2. La longueur L2 est fonction de la valeur de la rampe Doppler maximum admissible. Elle est choisie de telle sorte que l'erreur de linéarité reste négligeable, inférieure par exemple à la moitié de l'écart-type de l'erreur due au bruit. Cependant comme l'élargissement de la plage initiale a pour les fai- bles rapports signal/bruit toutes les chances de conduire à prendre en compte des phases corrigées présentant des sauts de phase multiples de rc, il est nécessaire pour s'accommoder de cette situation d'effectuer les opéra- tions suivantes. Tout d'abord, la plage en cours d'élargissement est supposée avoir une largeur e supérieure ou égale à L et inférieure à une largeur L2. Pour l'élargir d'une unité, les opérations consistent à calculer une valeur extrapolée (De par approximation linéaire à partir des valeurs courantes de A et B, à modifier ensuite la phase corrigée (1n, à prendre en compte en lui ajoutant ou en lui retranchant trds le nombre de fois nécessaire pour que la différence ln - (Iel soit minimale. La nouvelle valeur de (Dn ainsi obtenue est alors utilisée pour remettre à jour les coefficients de régression linéaire A et B. Le processus d'élargissement peut commencer par exemple par un prolongement de la plage vers la gauche pour ensuite la prolonger vers la droite. L'algorithme à utiliser peut alors être le suivant: Prolongement à gauche: INITIALISATION: Remettre les quantités So,, S, S2, Po, P1, A et B aux valeurs trou- vées à la position no. Faire n = no (indice du premier échantillon à considérer à gauche) Faire i = 0 (i vaudra 0, -1, -2,..., -(L2-2)/2) Faire ng= MAX(1, MIlN(no - (L2 - L) /2 + 1, N - (L2 -1))) (ce sera le dernier indice utilisé). ITERATION: a: Si n A = SO S2 - S2 Po =Po + Dn, P1 = P1 + i Dn A = P1 SO - PO SI A B - PO S2- P1 S1 A e: n=n-1 i=i-1 i=i - 1 Aller en a Prolongement à droite: INITIALISATION: Remettre les quantités So, Sl, S2, P0, P1, A et B aux valeurs trou- vées à la position no. Faire n = no + 1 (indice du premier échantillon à considérer à droite) Faire i =L + 1 (i vaudra L + 1, L + 2,..., L +(L2L)/2) ITERATION: a: Si n > ng + L2 - 1 alors Fin (prolongation terminée) b: (e = A i + B (extrapolation) c: Ramener (Dn le plus près possible de De par ajout(s) ou retrait(s) de 7c d: (mise à jour des coefficients de régression) SO = SO + 1, S1 + i, S2 + i2 A = SO S2 - S12 Po = Po + (n P1 P1 +i On A - P1 S0 - P S1 B - Po S2 - P1 S1 A e: n=n+1 i=i+1 Aller en a En fin d'exécution de l'algorithme, les coefficients de régression linéaire A et B valables pour une plage de largeur L2 contenant la plage ini- tiale de largeur L sont disponibles. La figure 5 montre les droites de régression linéaire obtenues dans le même référentiel que celui des figures 3A et 3B correspondant aux deux plages de l'exemple précédent (la plage initiale L étant repérée par deux traits en tirets parallèles). L'estimation finale de la phase de référence a lieu par l'exécution des étapes 8 et 9. Comme la longueur L2 a été choisie la plus grande pos- sible, le reste des estimations de phase est exécuté à l'aide d'une fenêtre glissante de largeur L2 pour corriger les sauts de phase. La fenêtre trouvée de largeur L2 trouvée à l'étape 7 incorpore les échantillons de rang ng à ng + L2 - 1 si L2 est inférieur ou égal au nombre N d'échantillons. L'étape 8 con- siste à faire glisser la fenêtre à gauche puis à droite. Ce glissement s'effec- tue suivant trois modes opératoires consistant, à calculer une valeur nomi- nale extrapolée (e de la phase, à effectuer une correction de chaque phase (DI)n et à effectuer une mise à jour de la régression linéaire. Pour le glissement à gauche l'algorithme à exécuter est le suivant: INITIALISATION: Remettre quantités So, S1, S2, P0, P1, A et B aux valeurs trouvées à la fin de détermination initiale de la fenêtre de largeur L2. Leurs valeurs sont: A = P1 So0 -Po S1 B =PO S2 - PlSl avec 0 =LL 2 L2(L2 + 1) S2=,2= L2(L2 + 1)(2L2 + 1) SoC= L2, S1 S = i = i=1 2 i=1 6 A = SO S2 - S12 ng+L2 -1 ng +L2 -1 Po = D 'i, P1 =1 (i - ng + 1)Oi i=ng i=ng Faire n = ng - 1 (indice du premier échantillon à considérer à gau- che) ITERATION: a: Si n c: Ramener On le plus près possible de (e par ajout(s) ou re- trait(s) de t d: (mise à jour des coefficients de régression) A= P1 SO -P0 SI A B = Po S2 -Pl S1 avec n+L2-1 n+L2-1 1 5 Po= 1 i, P1: (i-n+1)Di i=n i=n e: n=n-1 Aller en a Pour le glissement à droite l'algorithme à exécuter est le suivant: INITIALISATION: Remettre quantités So0, S, S S2, Po, P1, A et B aux valeurs trouvées à la fin de détermination initiale de la fenêtre de largeur L2. Faire n = ng + L2 (indice du premier échantillon à considérer à droite) ITERATION: a: Sin > N alors Fin (plus d'échantillons à droite) b: (e = A (L2 + 1) + B (extrapolation, pour i = L2 + 1) c: Ramener (n le plus près possible de (De par ajout(s) ou retrait(s) de oirds. d: (mise à jour des coefficients de régression) A P So -Po S1 A B = Po S2P1 S1 A avec n+L2-1 n+L2-1 Po: (Di, P1 (i - n + 1)( i=n i=n e: n=n-1 Aller en a Cependant pour ne pas dépenser en puissance de calcul, les cal- culs précédents peuvent encore être exécutés par un ajustement des diffé- rentes variables comme défini ci-après. L'algorithme pour le glissement à gauche est alors le suivant: INITIALISATION: Remettre quantités So, S1, S2, P0, P1, A et B aux valeurs trouvées à la fin de détermination initiale de la fenêtre de largeur L2. Faire n =ng -1 (indice du premier échantillon à considérer à gauche) Fairei=n-no ITERATION a: Si n A = S2 SO S12 Po Po + (Dn - Dn+L2, P1 = P1 + i (n - (i + L2) (t)n+L2 A = P1 S0 - Po SI B Po S2- P1 S1 A e: n=n-1 i = i- 1 i=i-l L'algorithme pour le glissement à droite est le suivant: INITIALISATION: Remettre les quantités S0o, S, S S2, P0, P1, A et B aux valeurs trou- vées à la fin de détermination initiale de la fenêtre de largeur L2. Faire n = ng + L2 (indice du premier échantillon à considérer à droite) Faire i = n-n0 ITERATION a: Sin > N alors Fin (plus d'échantillons à gauche) b: (e = A i + B (extrapolation) c: Ramener (Dn le plus près possible de I>e par ajout(s) ou retrait(s) de ir d: (mise à jour des coefficients de régression) S1 = S1 + i- (i + L2) = S1 - L2 S2 = S2 + i2 - (i +L2)2 = S2 - L2 (L2 + 2 i) A = S2 SO - S1 Po P0 + ()n - ([n+L2, P1 = P1 + i)n - (i + L2) (n+L2 A= P1 S0 -Po SI B Po S2 -P1 S1 A e: n=n-1 i=i-1 Aller en a Les approximations linéaires locales de largeur L2 ainsi obtenues permettent d'obtenir un jeu de N phases corrigées On o il n'y a plus de sauts de phase comme cela est représenté à la figure 6. Le dernier traitement consiste, compte tenu du fait que le décalage Doppler varie dans le temps, à effectuer une régression quadratique sur les phases corrigées, ce qui suppose que le décalage Doppler varie quasi- linéairement. Ce traitement consiste à approximer les phases par: (::n > A n2 +B n + C, n = 1...N et à minimiser l'erreur quadratique moyenne suivante: N - (Dn _(A n + B n + C)) La solution est donnée par:= La solution est donnée par: Po(S2 - S, S3) + P1(so S3 - S S2) + P2(S2 - So s2) A= P0(S1 S4 - S2 S3) + P1 (S2 - S0 S4) + P2 (S0 S3 - S1 S2) B= à Po(SC - S2 s4)+ P1(Sl, S4 - S2 s3)+ P2(S - S2 s3) A avec N N(N+ 1) B 2 N(N + 1) (2 N + 1) So=N,S1= = 2 2= li2 =6 i=1 i=1 N N2(N+ 1)2 N 4 N(N+1)(2N+1)(3 N2 +3N-1) 53= ÉZz = __4 S4= _i =30 i=1 i=1 3 A S - 2 s s s + So s 3 + S2 S4 - So S2 S4 N N N P0 = " n, P1= I.n (:)n, P2 = n 2 On n=1 n=1 n=1 et la valeur minimale obtenue est: Emin = n - A P2 - B P1 - B Po Les phases corrigées initiales, et l'approximation parabolique finale obtenues par le dernier traitement sont représentées figure 7. Le coefficient B est égal à l'incrémentation de phase estimée pen- dant la durée d'un symbole, pendant le symbole (fictif) d'indice 0, soit 2 Af T. Sa variance Caf2 est de l'ordre de 192 (a /a)2 /N3 (en radians), a étant l'amplitude du signal et a l'écart-type du bruit. Le rapport S/B est donné par: S/B=a2/(2a 2). L'écart-type est approximativement égal à CAf =(1/2(2 = T)) (96/(N3 S / B)) 12. A titre d'exemple, pour 1 Kb/s, S/B = 3 dB et N = 2*100 bits transmis (100 bits utiles), oAf est de l'ordre de 0.4 Hz. Cet exemple montre que, lorsqu'aucune erreur n'a été commise, la valeur minimale de l'erreur quadratique moyenne est faible. Pour un rapport S/B minimum de 3 dB, elle est en général inférieure à 0.2 radian2. Par contre, lorsqu'une erreur est commise par une correction abu- sive de 7rrds, cette valeur minimale est proche de ou supérieure à la valeur théorique re2/12 = 0.82 radian2, qui correspond à une erreur uniformément répartie entre -r/2 et +r/2. Cette constatation permet de sélectionner les blocs reçus en fonc- tion de Emin. Si cette quantité est inférieure à 0.3 (par exemple), le bloc est déclaré acceptable, sinon il est rejeté. Naturellement ce traitement peut tou- jours être recommencé avec une valeur de L plus faible (par exemple moitié moindre) afin d'augmenter les chances de trouver une plage de longueur L sans transition. La sélection permet aussi d'éliminer les messages ayant un trop mauvais rapport signal/bruit, comme le montre la figure 8 o S/B vaut 0 dB avec une valeur de Emin supérieure à 0.3 (le signal d'entrée étant le même que précédemment, mais o la puissance du bruit a été multipliée par 2). La phase finale de démodulation consiste à soustraire la phase de référence du signal reçu et donc à calculer: Sdn = Sn expj(-(An2 + B n + C)), n = 1...N Une première constellation d'échantillons complexes ainsi obtenues par démodulation sur deux voies en quadratures est représentés à la figure 9A dans le cas d'un rapport signal à bruit SNR = 3 dB. Il n'y a qu'un bit erroné pour un écart moyen de 0,248943.. D'autres constellations obtenues avec des rapports S/B de 0 dB (moins bon), puis 6 et 10 dB (donc meilleurs) sont montrées respectivement aux figures 9B, 9C et 9D. Dans les calculs précédents, les variances des coefficients de * régression dépendent de la façon dont la phase On est approximée. Dans l'approximation quadratique o ODn est de la forme (Dn =A3n + B3n + C3 + en les variances des coefficients A3, B3 et C3 sont données par les relations: V(A) 180 302 N5 V(B) 192 2 N3 FI V(C) 9 et 2 N E aE représente l'écart type et est voisin de a 1 (2S / B)1/2 a Lorsque l'approximation de IOn suit une loi de régression linéaire de la forme: ()n = A2 n + B2 + n Les variances sur les coefficients A2 et B2 sont données par les relations: V(A) 12 a2 N3 V(B) 4 et 2 - a'- N Dans l'approximation de (n par une constante de la forme: on A1 + Eb on obtient: V(A) 1 a2 N Les résultats précédents montrent que la variance de la constante est 4 fois plus importante dans l'approximation linéaire par une constante et 9 fois plus élevée pour l'approximation quadratique et qu'il faut travailler sur davantage d'échantillons à mesure que le nombre de coefficients de l'approximation augmente pour conserver les mêmes performances. Si les phases corrigées sont de la forme (I)n = A3n2+B3 n + C3 Le choix de la longueur de la fenêtre L2 est déterminé en fonction de l'erreur maximale eein obtenue par approximation linéaire, suivant la relation eein = A3(L2-1)(L2 - 2) /6 Avec une rampe Doppler de dAf Hz par seconde et une durée de symbole T, A3 est égale à A3 =T2 dAf /2 Dans ces conditions pour obtenir une erreur eein inférieure à X fois l'écart type des phases dû au bruit a(, il suffit pour L2 assez grand que L2 vérifie la relation L2 T dAf 2/B Avec S/B = 3 dB, un rythme symbole de 1000 Hz et une rampe Doppler de 100 Hz par seconde, la relation précédente montre que L2 doit être inférieur à 75 pour que l'erreur de linéarité ne dépasse pas le 1/10ioe de l'erreur due au bruit. Bien que très performant, le procédé décrit précédemment peut encore être amélioré. Par exemple, si dès la phase de recherche de la position de départ optimale, l'erreur quadratique minimale, est trop élevée, un autre essai peut encore être tenté avec une fenêtre de largeur L réduite en espérant trouver cette fois-ci un résultat convenable. Comme il subsiste toujours une ambiguité de phase de 7c dans la détermination des phases de référence, il n'est pas possible de savoir a priori si la phase du premier échantillon de signal est 0 ou t puisque le bit modulant peut avoir l'une ou l'autre valeur 0 ou 1. Cette ambiguité peut être levée par tout procédé connu, comme par exemple celui consistant à insé- rer une séquence de K bits connue, ou un codage différentiel des données avant émission. Bien entendu, le procédé et le dispositif qui viennent d'être décrits dans le cadre d'une modulation biphase s'appliquant aussi bien pour d'autres types de modulation de phase. Par exemple, pour une modulation à 4 états de phase, appelée QPSK dans la littérature, les dispositifs 1, 3 et 4 peuvent être adaptés en corrigeant les phases et différences de phase dans le dispositif 3 pour les ramener entre -n/4 et +7/4; et pour corriger éventuellement OIn dans les dispositifs 3 et 4 par pas de 7c/2 et non pas r. Pour une modulation MSK ou OQPSK, la seule modification à apporter consiste alors à multiplier par j tous les échantillons d'ordre impair par exemple. Enfin, le procédé peut aussi être adapté pour affiner la synchronisa- tion du récepteur. En disposant de K échantillons par symbole, au lieu de 1 comme il a été supposé jusqu'ici, et donc de N.K échantillons au total, il suffit alors d'effectuer une démodulation complète à partir des échantillons 1, K+1,.. etc., fournissant une valeur Emin(l); une démodulation complète à partir des échantillons 2, K+2,..., etc., fournissant une valeur Emin(2); etc. et une démodulation complète à partir des échantillons K, 2 K,..., etc., fournissant une valeur Emjn(K). Dans ce cas seul le résultat de la démodulation de rang k qui fournit la plus faible valeur Emin(k) est à retenir et le rang k désigne alors la position qui est la plus proche de la position de synchronisation idéale. Il est aussi possible éventuellement d'effectuer une interpolation de la valeur Emin au voisinage de la kème position pour affiner ce résultat. REVENDICATIONS 1. Procédé de démodulation cohérente par bloc de suites de signaux complexes Sn modulés en phase suivant des états de phase déterminés caractérisé en ce qu'il consiste à extraire (5) la phase de chacun des signaux complexes en supprimant les sauts de phase inconnus, à recher- cher (6) dans la suite des phases obtenues une première plage de phase de longueur déterminée L pour laquelle la variation de phase est la plus linéaire possible, à déterminer (7) une deuxième plage de phases de longueur L2 supérieure à L correspondant à une approximation linéaire plus précise, à effectuer (8) itérativement des régressions linéaires à l'aide d'une fenêtre glissante de longueur égale à la deuxième plage de phase L2 pour suppri- mer les sauts de phase résiduels pour effectuer (9) une régression quadrati- que sur les phases restantes et pour calculer une estimation globale des phases de référence des signaux de chaque bloc, et à effectuer (10) une comparaison de l'erreur quadratique moyenne entre les phases de référence estimées et les phases brutes obtenues par régression linéaire. 2. Procédé selon la revendication 1 caractérisé en ce que la sup- pression (5) des phases inconnues consiste pour une modulation à deux états de phase (0-n) à corriger les phases des échantillons Sn en les rame- nant dans l'intervalle de phase [-2' +2] à calculer les différences de phase successives en les ramenant elles aussi dans l'intervalle de phase [- ' +2] et à reconstituer une suite (On de valeurs de phases représentati- ves de la phase de référence du signal. 3. Procédé selon l'une quelconque des revendications 1 et 2 carac- térisé en ce qu'il consiste pour rechercher la première plage de phase, à calculer par régression linéaire sur la longueur L l'erreur quadratique moyenne de phase en prenant successivement chacun des signaux Sn, comme origine et à choisir pour origine celui qui donne l'erreur quadratique moyenne minimum. 4. Dispositif pour la mise en oeuvre du procédé selon l'une quelcon- que des revendications 1 à 3 caractérisé en ce qu'il comporte un dispositif de mémorisation (2) des signaux complexes couplé à un dispositif (3) d'es- timation de la phase de référence du signal et un dispositif (4) couplé au dispositif (3) d'estimation de phase pour reconstituer les échantillons com- plexes et décoder l'information transmise par les échantillons.