FR 2706711 A1 19941223 FR 9307324 A 19930617 PROCEDE ET DISPOSITIF DE DEMODULATION DE SIGNAL NUMERIQUE La présente invention concerne un procédé et un disposi- tif de démodulation de signal à modulation de phase à M états d'une porteuse, M étant un nombre entier, généralement pair, et notamment à M = 2m états (m étant un entier). L'invention trouve une application particulièrement importante, bien que non-exclusive, dans les sytèmes de transmission numériques en modulation MDPQ ou MDP8 d'une porteuse hyperfréquence (c'est-à-dire ayant une fréquence d'au moins 3 GHz). On peut notamment citer, à titre d'appli- cation, la transmission simultanée de plusieurs programmes de télévision numérique sous forme d'un multiplex dans un canal de diffusion de satellite, dans la bande des 12 GHz. Cependant, comme on le verra plus loin, l'invention s'applique à la démodulation de tout signal utilisant la modulation MDP d'une porteuse et trouve un intérêt chaque fois que l'on se trouve dans des conditions défavorables de récupération de cette dernière. Les systèmes de transmission et de diffusion sur un canal de satellite utilisent souvent une modulation de type MDP (MDPQ ou MDP8 en particulier) avec un codage préalable de protection contre les erreurs qui augmente le débit requis. On utilise en particulier des codes de protection convolutifs ou des codes concaténés (associant un codage convolutif à un codage de Reed- Solomon) qui éliminent quasi-complètement les erreurs même pour des rapports Eb/N0 (Eb étant l'énergie transportée par bit et No le niveau de bruit) aussi faibles que 2 à 4 dB. Le rendement des systèmes de codage fait correspondre ces valeurs à des rapports Ec/N0 (Ec étant l'énergie par bit transmise après codage) voisins de 0 dB. On rencontre le même problème lorsqu'on utilise des modes de codage en treillis qui permettent de conserver l'effica- cité spectrale d'un système non-codé, en augmentant le nombre d'états de modulation. Les systèmes de codage de protection contre les erreurs, qui augmentent le débit à transmettre à nombre de bits utiles donné, compliquent la récupération de la porteuse lors dé la démodulation à la réception et bien souvent ce sont les performances de cette récupération de la porteuse qui limitent les performances globales du système. L'invention concerne plus particulièrement les procédés de démodulation suivant lesquels on convertit le signal reçu en un signal en bande de base sur deux voies en quadrature à l'aide d'un signal fourni par un oscillateur local à la fréquence de porteuse (généralement après transposition à une fréquence intermédiaire) et on numérise ce signal en bande de base en l'échantillonnant à la fréquence d'horloge du signal modulant. Classiquement on récupère la porteuse, généralement transposée, de façon "aveugle", c'est-à-dire sans connaître les instants d'ouverture maximale du diagramme de l'oeil. La réduction de performances due à cette récupéra- tion en "aveugle" n'est pas fatale dans le cas o l'oscilla- teur local de transposition est suffisamment stable pour que la plage d'acquisition en fréquence requise pour la récupéra- tion certaine de la porteuse soit étroite. Ce cas est celui des décodeurs professionnels. Il n'en est pas de même dans le cas de récepteurs grand public, o l'oscillateur local de transposition est suscep- tible de dériver de plusieurs MHz: la récupération de porteuse sur une plage aussi large devient alors délicate et peut se révéler impossible dans le cas d'un signal présentant un rapport Ec/N0 faible, ce qui est fréquemment le cas des signaux à forte redondance de codage de protection contre les erreurs. On pourrait penser que de meilleurs résultats sont obtenus en récupérant conjointement la porteuse et l'horloge à l'aide d'une boucle de récupération fortement non-linéaire et à rebouclage. Malheureusement les systèmes correspondants sont complexes et présentent des risques de divergence. La présente invention propose d'effectuer la récupération d'horloge avant la récupération de porteuse, de façon à permettre un échantillonnage aux instants d'ouverture optimale du diagramme de l'oeil et utilise pour cela le fait que la correspondance entre la fréquence d'horloge locale et la fréquence d'horloge du signal reçu est généralement très précise, contrairement à la correspondance entre fréquence porteuse transposée et fréquence d'horloge locale du passage en bandes de base. La raison en est que l'horloge locale, fonctionnant à fréquence fixe, peut être stabilisée par quartz et que la précision de l'horloge du signal est très peu affectée par la transmission, puisque la seule source de dégradation est l'effet Doppler. Au contraire l'oscillateur local de transposition de fréquence doit fonctionner à fréquence ajustable. L'invention propose en conséquence un procédé du genre ci-dessus défini caractérisé en ce qu'on récupère l'horloge du signal par calage de la phase d'une horloge locale, qui est stable à la fréquence de bit, par recherche du minimum d'interférence intersymbole, avant récupération de la porteu- se. Il est ensuite avantageux de commander la fréquence de l'oscillateur local jusqu'à ce que l'erreur de fréquence sur le signal échantillonné et numérisé soit minimale. Dans un premier mode de réalisation, le minimum d'inter- férence intersymbole (ISI) est obtenu par recherche du minimum d'une fonction de la valeur du signal qui est non- linéaire et qui s'annule pour au moins quatre points de la constellation des états de phase (quatre points dans le cas d'une modulation MDP 4). Cette fonction peut notamment être de la forme F[z(t).z*(t-T, )],le produit étant destiné à éliminer la phase inconnue et F étant une fonction ayant la propriété de présenter des minima, de préférence nuls, pour tous les points de la constellation. Par exemple la fonction F peut être le minimum de la distance entre un point reçu et l'ensemble des points de la constellation; z=x+jy représente la valeur du signal sur les voies Iet Q, et T_ est la période de symbole. Le produit prend périodiquement une valeur dont le module est égal au carré du module du signal et la phase prend une valeur de la forme 2kr/M, k étant un entier, aux instants d'échantillonnage idéaux. On obtient ainsi une constellation fixe. Dans un autre mode de réalisation, on échantillonne le signal après conversion en bande de base, on mesure la varia- tion de l'énergie moyenne du signal en fonction de la phase d'échantillonnage et on sélectionne la phase pour laquelle l'énergie est maximale. Cette solution est utilisable chaque fois que l'excès de bande par rapport à la limite de Nyquist est suffisante ou que la décroissance en limite de bande (roll-off) est progressive. L'invention propose également un dispositif de démodula- tion permettant de mettre en oeuvre le procédé ci-dessus défini. Ce dispositif comprend, de façon classique, des moyens à au moins un oscillateur de transposition du signal reçu à une fréquence intermédiaire; des moyens de conversion en bande de base sur deux voies en quadrature, utilisant un oscillateur local à fréquence ajustable autour d'une valeur nominale; des moyens de conversion numérique des deux voies comportant une horloge à phase ajustable, stabilisée à la fréquence de bit du signal entrant; et des moyens de récupération de porteuse et de démodulation cohérente des deux voies. Suivant un aspect de l'invention, le dispositif comprend des moyens pour ajuster la phase de l'horloge à la valeur correspondant au minimum d'interférence intersymbole et des moyens pour commander l'oscillateur local de passage en bande de base, l'ajustage de la phase d'horloge s'effec- tuant par analyse du signal avant récupération de la porteu- se. L'invention sera mieux comprise à la lecture de la description qui suit de modes particuliers de réalisation donnés à titre d'exemples non-limitatifs. La description se réfère aux dessins qui l'accompagnent, dans lesquels: - la figure 1 est un schéma de principe montrant la constitution d'ensemble d'un démodulateur permettant de mettre en oeuvre l'invention; - la figure 2 est un schéma synoptique d'un module de récupération d'horloge suivant l'invention utilisable dans le démodulateur de la figure 1, ainsi qu'une constitution possible de la partie amont effectuant la transposition de fréquence; - la figure 3 montre une disposition que peut prendre la constellation des points représentatifs des états de phase, dans le cas d'une modulation MDPQ, dite aussi MDP4; - la figure 4 montre l'allure du diagramme de l'oeil d'un signal modulant MDPQ; - la figure 5 est un schéma destiné à montrer une allure possible de l'évolution de la variance du signal servant à la récupération d'horloge, en fonction du temps; - la figure 6, similaire à la figure 1, montre un autre mode de réalisation; - la figure 7 montre l'allure d'une fonction intervenant dans la recherche de la phase correcte d'horloge; - la figure 8 est un schéma d'un mode possible de commande automatique de la fréquence d'un oscillateur de transposition, permettant de récupérer la fréquence porteuse aussi longtemps qu'elle rentre dans une plage d'accrochage ajustable; et - la figure 9 montre l'allure d'une fonction intervenant dans la récupération de porteuse. Le démodulateur dont la constitution de principe est montrée en figure 1 peut notamment être utilisé pour démodu- ler un signal diffusé par satellite, en modulation MDP. On peut le regarder comme constitué d'une tête de réception hyperfréquence et de transposition 8 qui alimente: - un module de récupération d'horloge 36, un module de commande automatique de fréquence 42, - des circuits 10 de récupération de porteuse, pouvant, dans certains cas, être précédés de moyens de filtrage complémentaires et d'interpolation. Les divers composants ci-dessus seront maintenant décrits successivement, éventuellement avec plusieurs modes d'exécu- tion possibles. La tête 8 peut avoir la constitution relativement classique montrée en figure 2. Elle comporte alors une partie hyperfréquence, par exemple destinée à recevoir un signal dans la bande des 12 GHz. Ce signal traverse d'abord une partie hyperfréquence ayant deux circuits de transposition de fréquence en cascade 11 et 12, permettant de passer successivement à une première fréquence intermédiaire puis à une seconde fréquence intermédiaire qui est fixe, par exemple à 140 MHz. Dans le cas d'une tête haute fréquence grand public, par exemple prévue pour un système dit VSAT, la précision et la stabilité des oscillateurs locaux (non représentés sur la figure 1) ne permettent pas de maintenir la fréquence intermédiaire à mieux que quelques MHz. Le signal transposé à la seconde fréquence intermédiaire est appliqué à un filtre 14 de sélection de canal, qui, dans le cas d'une émission depuis un satellite, a généralement une largeur de 36 MHz environ. Il sera surtout question par la suite de signaux dont la modulation numérique est à une cadence de 46 Mbits/s, correspondant à 23 MHz de fréquence symbole l/Ts. Lorsque, ce qui est le cas le plus fréquent, l'émetteur est muni d'un filtre demi-Nyquist, le filtre à large bande est suivi d'un filtre demi-Nyquist 16 qui peut cependant être reporté dans la chaîne de traitement numérique en bande de base. Le signal est ensuite appliqué à des moyens 18 de conversion en bande de base, comportant un oscillateur 20 dont la fréquence doit être aussi proche que possible de la fréquence du signal transposé (appelée ci-après fréquence de porteuse), fournissant deux ondes sinusoidales ayant un déphasage mutuel de 900. Ces ondes sont appliquées à des multiplieurs respectifs 22 et 24 destinés à fournir deux signaux en quadrature I et Q. Dans le mode de réalisation avantageux illustré sur la figure 2, l'oscillateur local 20 de passage en bande de base est à fréquence commandable par un signal appliqué sur une entrée 26. Sur chacune des voies I et Q est ensuite interposé un convertisseur analogique-numérique 28 ou 30, dont la fré- quence d'échantillonnage est commandée par une horloge locale 32 à la fréquence de bit. Cette horloge 32, pouvant être stabilisée par quartz, reproduit généralement la fréquence de bit à 10-5 près, donc avec une précision très élevée. Par asservissement, l'erreur à long terme peut même être rendue nulle. L'horloge 32 présente dans ce cas une entrée 34 de réglage de la phase. Les valeurs c, cj, -c, -cj prises par le signal aux instants d'échantillonnage peuvent ne pas correspondre aux points nominaux (points sur les axes, comme indiqué sur la figure 3, pour une modulation MDP4). Sur cette figure 3, un quadrant correspond à une erreur de 1/4T., c'est-à-dire 5,75 MHz pour une modulation MDP4 à 46 Mbits/s. Si l'horloge a été exactement récupérée, les instants d'échantillonnage se placent en TI, 2... à intervalles Ts sur la figure 4 et le traitement ultérieur fait appel à des points supplémentaires T, obtenus par sur-échantillonnage à fréquence 4f, pour arriver à une période d'échantillonnage Ts/4. Conformément à un aspect de l'invention, les signaux I et Q échantillonnés sont appliqués directement à des moyens de récupération d'horloge et de commande de la phase de l'horloge d'échantillonnage 32, permettant d'ajuster les instants d'échantillonnage de façon qu'ils soient placés de façon optimale par rapport au diagramme de l'oeil. Dans le mode de réalisation montrée en figure 2, l'inven- tion utilise les propriétés de la fonction z (t). z' (t-Ts) o z (t) = x (t) + jy (t) représente la valeur du signal sur les voies I et Q. Cette fonction prend périodiquement les valeurs a, aj, -a et -aj, o a est le carré du module commun de c, -c, cj et -cj. Le produit complexe peut être représenté sur la figure 3 par des points qui conservent une position angulaire fixe par rapport à un axe, cette position étant représentative de l'erreur de fréquence et pouvant être utilisée pour une commande automatique de fréquence. La démodulation différentielle éliminant la phase, qui est inconnue, on peut utiliser cette propriété pour identi- fier la phase optimale d'échantillonnage, qui correspond à la dispersion minimale. Pour cela on peut faire appel à une fonction non-linéaire F(z) qui s'annule pour les quatre points de la constellation MDP4 (ou plus généralement pour tous les points de la constellation, pour des modulations autres que MDP4). Cette fonction peut être par exemple le minimum des carrés des distances euclidiennes à l'un des quatre points de la constellation, dans le cas du MDP4. Pour utiliser cette approche, on peut adopter la suite aléatoire Xn(t) définie par: X,(t) = F[z(t+nTs).z*(t+(n-1)Ts)] o l'étoile * identifie la conjuguée complexe. La valeur moyenne de cette suite pour chaque valeur de t est une fonction périodique G(t) qui a l'allure montrée sur la figure 5. Les minima de cette fonction de variance (repré- sentative de la dispersion) correspondent aux instants optimaux d'échantillonnage. Mais ils ne sont nuls que si le diagramme de l'oeil est parfait. Une méthode simple pour rechercher les minima, valable dans tous les cas, consiste à exploiter la symétrie de la fonction G(t). On peut pour cela utiliser la fonction H(t) définie par H(t) = G(t+Ts/4) G(t-Ts/4) Cette fonction s'annule pour la phase optimale et est maximale pour t = + T5/4 et t = - T/4. Cette approche peut être implémentée avec un premier module matériel ou logiciel de calcul 35 (figure 2) qui calcule F[z(t). z*(t-T.)], et un second module 38 qui calcule H pour les valeurs successives de t, par différence entre les variances pour deux points écartés de T,/2, effectue une moyenne et commande l'horloge 32 pour annuler H et arriver à l'équilibre des deux points montrés en figure 5. Le schéma représenté comporte ensuite le module 10 qui reçoit lui aussi les échantillons numérisés I(f3) et Q(fS) à la fréquence symbole; ils sont traités de façon classique dans un module 40 de récupération de porteuse, portant par exemple le signal à la puissance 4 (ou utilisant une boucle de Costas) et de démodulation cohérente ou différentielle, non représenté. Un filtrage peut être effectué en amont du circuit de récupération. Au surplus, une interpolation est nécessaire lorsque l'échantillonnage s'effectue à la période T3 (figure 4). Un module 42 de commande automatique de la fréquence de l'oscillateur local 20 et d'adaptation à la fréquence de porteuse (figure 1) est généralement prévu. Une constitution simple possible de ce module sera donnée plus loin. La figure 6, similaire à la figure 2, montre un autre mode de réalisation, plus simple que le précédent et en général plus avantageux; dans le cas de la figure 6 la récupération d'horloge est fondée sur des considérations énergétiques. Une analyse théorique justifie la validité de l'approche utilisée. Elle utilise une fonction P(t) qui représente la répartition de l'énergie selon la phase d'échantillonnage modulo T,: P(t) = E(Js(t) 2), E(.) étant l'espérance mathématique. Le calcul montre que P(t) est une fonction périodique, de période T3, dont l'expression est légèrement différente suivant le mode de filtrage. Dans le cas d'un signal filtré demi-Nyquist provenant directement de l'émetteur avec un facteur de décroissance ou de "roll-off" a, on a: P(t) = 1 + (2"/n) Cos 27 t/Ts Dans le cas d'un signal filtré Nyquist avec un "roll-off" ", on a: P(t) = 1 - (a/4) + a/4 Cos 2nt/Ts, Dans tous les cas, la fonction P(t) a l'allure montrée en figure 7, avec un maximum PM et un minimum Pm écartés de T,/2. Le rapport p = PM/Pm dépend de la décroissance a et de la fonction de filtrage (demi-Nyquist ou Nyquist). Les valeurs suivantes peuvent être considérés comme représentati- ves, a étant défini comme indiqué à droite de la figure 7 et la limite 1 correspondant à une borne de la bande de Nyquist. a= 0,3 c = 0,5 = 1 demi-Nyquist 1,4721 1,9339 4,5039 Nyquist 1,1765 1,3333 2 Même pour des facteurs de décroissance a faibles, corres- pondant à un filtre de canal 14 à coupure franche, les varia- tions de P(t) sont significatives, surtout dans le cas du filtrage demiNyquist. On aura donc intérêt à utiliser le signal reçu après filtrage de canal, mais avant filtrage demi-Nyquist de réception, dans le cas habituel o le filtrage Nyquist complet est partagé entre émetteur et récepteur.Le fait que le bruit soit alors élargi avant prise en compte du signal dans le calcul de P(t) n'est qu'un inconvénient mineur, car le rapport entre les bandes de bruit du filtre de canal et du filtre demi-Nyquist de réception (filtres 14 et 16 sur la figure 2) est de (1 + a). Le bruit affectant le signal ajoute donc seulement, à la fonction P(t) une composante continue légèrement plus élevée en cas d'utilisation du signal avant filtrage demi-Nyquist de réception. Un avantage supplémentaire de l'utilisation du signal prélevé en amont du filtre demi-Nyquist de réception est qu'une erreur de fréquence affectant le signal reçu ne modifie en rien la fonction P(t), aussi longtemps que le signal reçu n'est pas sensiblement affecté par le décentrage dans le filtre de canal 14. En effet, aussi longtemps que le signal reste dans la bande passante du filtre de canal, le décentrage ne se traduit que par une multiplication par un exposant complexe. Comme dans le cas précédent les instants d'échantillonnage optimal correspondent aux minima de P(t). Pour déterminer les minima, on peut encore utiliser la différence de deux valeurs pour des instants écartés de TJ/2, ici la fonction H: H(t) = P(t-T,/4) - P(t+Ts/4) Cette fonction est symétrique par rapport à l'origine, s'annule pour T,/2 et - T,/2 et présente des maxima à T./4 et - T./4 par rapport à la phase d'horloge correcte. Le second membre de la fonction H(t) peut être estimé par deux calculs de puissance moyenne, effectués l'un sur quelques échantillons pairs successifs à la fréquence symbole l/T,, l'autre sur le même nombre d'échantillons impairs et par soustraction. Les calculs de puissance P sont à effectuer aux phases nominales - T5/4 et + T$/4. Ils permettent d'obtenir P2n = 22n + Q 2n et P2n-l,. Le passage à zéro de P2n - P2n-, représentatif de H(t), constitue un critère d'asservissement de la phase de l'hor- loge 32. Le mode de réalisation illustré sur la figure 6 (o certains filtres passe-bas sont omis pour plus de clarté) correspond à un signal MDP4 (ou QPSK) dans la bande normali- sée satellite, utilisant le signal avant filtrage demi- Nyquist pour récupérer l'horloge. Les signaux analogiques de sortie des multiplieurs 22 et 24 sont appliqués aux convertisseurs analogique/numérique 28 et 30. Dans le cas d'une transmission sur un canal satellite, le filtre de canal 14 a une bande de 36 MHz et tolère généralement un facteur de décroissance a d'environ 0,5. L'oscillateur 20 est à fréquence fixe, le dispositif de la figure 6 comportant une commande automatique de maintien de la seconde fréquence intermédiaire, qui sera décrite plus loin. La phase de l'horloge locale 32, dont la fréquence nominale est de 48 MHz, est commandée par la boucle de récupération d'horloge. La boucle de récupération d'horloge utilise les séquences d'échantillons Is et Qs qui apparaissent, à la sortie des convertisseurs 28 et 30, au rythme de deux fois la fréquence symbole l/Ts, c'est à dire à la fréquence bit. On calcule alors séparément la puissance instantanée à la cadence d'apparition des échantillons et on calcule la puissance moyenne, d'une part pour les échantillons pairs, d'autre part pour les échantillons impairs. Sur la figure 6, la branche de calcul pour les échantillons pairs I2n, compor- te deux éléments d'élévation au carré 47 (pouvant en général accumuler quelques échantillons) et un additionneur 46. La branche de calcul pour les échantillons impairs est similai- re. Un soustracteur 48 permet alors d'obtenir un signal de sortie représentatif en permanence de la fonction H ci- dessus. Ce signal numérique est utilisé pour commander l'horloge 32, à l'aide d'un convertisseur numérique-analogi- que 50 et d'un filtre passe-bas de boucle 52 effectuant une moyenne sur quelques centaines d'échantillons. Lorsque les échantillons I, et Qn sont numérisés sur 8 bits, il suffit généralement de soustracteurs sur 12 bits, en effectuant une moyenne sur un nombre d' échantillons faible (quatre par exemple), mais suffisant pour réduire la vitesse requise du CNA 50. Dans la pratique d'ailleurs, l'accumula- tion du carré des échantillons peut être effectuée par un filtre passe-bas. On voit que le schéma de la figure 6 commande l'horloge 32 jusqu'à amener à A la fonction H. Les composantes conti- nues qui apparaissent dans les signaux sont automatiquement éliminées dans le traitement. Les échantillons numérisés In et Qn sont également appliqués au module 42 de commande automatique de fréquence dont une constitution possible est donnée en figure 8 et au module 10. Ce dernier effectue le filtrage demi-Nyquist de réception et une interpolation; cette interpolation, néces- saire pour disposer des échantillons à tous les instants d'échantillonnage y, peut être combinée au filtrage demi- Nyquist dans un filtre transverse unique. Les échantillons d'origine et ceux obtenus par interpolation sont utilisés pour la récupération de porteuse et la restitution des données, par exemple par démodulation différentielle. Il est avantageux d'utiliser le signal de puissance pour effectuer une commande automatique de gain. Dans le cas représenté sur la figure 6, les puissances des échantillons pairs et impairs déjà disponibles sont additionnées en 53 et utilisées pour commander le gain d'un amplificateur 54 avant transposition de fréquence, par l'intermédiaire d'un conver- tisseur numérique-analogique 56 et d'un filtre passe-bas de boucle 58. On décrira maintenant une boucle de commande automatique de la fréquence de l'oscillateur de transposition du circuit 12 destinée à un signal à modulation MDP à M états. A l'entrée du bloc 42, on dispose d'une séquence d'échan- tillons complexes z, à la fréquence symbole l/Ts, qu'on peut écrire: zn = cn A exp (i+), ou A et 4 représentent l'amplitude et la phase à un ins- tant t, cn représente la valeur du symbole transmis, qui est: 1 ou - 1 en MDP 2; 1, i, -1 ou -i en MDP4. A l'entrée du bloc 42, la porteuse n'est pas encore récupérée et g est donc inconnu. De plus le signal peut être affecté d'une erreur Af sur la fréquence transposée, par suite de la dérive de l'oscil- lateur 60 (fig. 6). Il est donc souhaitable de corriger Af, qui est relié à la phase 4 et à la phase idéale Jo par: 4 = 4o + Af. Ts La boucle de contrôle de fréquence montrée en figure 8 utilise, pour estimer l'écart de fréquence Af, un discrimi- nateur à ligne à retard permettant de donner un retard quelconque de q temps symbole Ts parmi plusieurs retards successifs. On peut ainsi déterminer une grandeur sn: Sn = Zn* Z*n-q = A2 C Cnq exp (2ilq Af Ts) Dans ce cas, d'une modulation MDP à M états, les symboles sont de la forme exp (2iz m/M). Dans le cas cn. cnq prend uniquement des valeurs de la forme exp (2i7 m/M) et le problème se ramène à une estimation d'erreur de phase, que l'on sait effectuer. Aussi, grâce à une démodulation différentielle dans le bloc 57, on peut transformer l'estimation d'erreur de fréquence en une estimation d'erreur de phase. Il existe de nombreuses implémentations permettant de remplir cette fonction, telles que la boucle de Costas et ses variantes. On cherchera dans tous les cas à obtenir une caractéris- tique de la forme Ve =K.sin MA# o A# = 2zqAfT, Ve est le signal de sortie et K est une constante. Le signal d'erreur Ve est périodique et, sa période est définie par la relation Mq Af T = 1 La loi de variation du signal de sortie en fonction de l'erreur de phase Af est alors du genre montrée en figure 9 et fait apparaître que la boucle ne peut s'accrocher au delà de la bande -Z/MqT +7r/MqT5. Elle peut être ajustée de façon à être adaptée à la stabilité de l'oscillateur 60 par un choix approprié du nombre de retards q.T, possibles. Par exemple pour une modulation MDP4 à 40 Mbit/s, la plage d'accrochage pour q=l (un seul retard possible) est de + 2,5 MHz. La largeur de cette plage est inversement propor- tionnelle à q mais la sensibilité et la précision sur la fréquence augmentent proportionnellement à q. Dans le mode de réalisation illustré sur la figure 8, destiné à un signal MDP4, q peut varier de 1 à 4. Les signaux In et Qn sont adressés sur des lignes à retard respectives à quatre éléments 62 et 64, pouvant être constituées par des registres à décalage commandés à la fréquence symbole l/T,. Le nombre q de temps symbole dont le signal est retardé est sélectionné à l'aide d'une entrée 66. Un module de calcul 68, logiciel ou matériel, calcule Zz.Z*n-_q, qui s'écrit sous forme du produit d'une constante R par i@. Ce produit est utilisé par un second module 70 pour estimer 0, c'est-à-dire l'angle polaire de zn. z*n-q, et en déduire R sin 4 d (dans le cas d'une modulation MDP4). Ce calcul peut être fait de façon grossière et se limiter au choix de la valeur la plus proche dans une table. Le résultat, obtenu par exemple avec le même nombre de bits que I et Q. est utilisé après conversion numérique-analogique 72 pour commander un filtre de bande de commande de la fréquence de l'oscillateur 60. REVENDICATIONS 1. Procédé de démodulation de signal à modulation de phases à M états d'une porteuse, M étant un nombre entier, suivant lequel on convertit le signal en bande de base sur deux voies en quadrature par multiplication par un signal d'oscillateur local (20) à la fréquence de porteuse et on numérise le signal en bande de base par échantillonnage à la fréquence d'horloge du signal modulant, caractérisé en ce qu'on récupère l'horloge du signal par calage de la phase d'une horloge locale (34) stabilisée à la fréquence de bit par recherche du minimum d'interférence intersymbole, et cela avant récupération de porteuse. 2. Procédé selon la revendication 1, caractérisé en ce qu'on obtient le minimum d'interférences intersymboles (ISI) par recherche du minimum d'une fonction F de la valeur du signal qui est non-linéaire et qui s'annule pour les M points de la constellation des états de phase. 3. Procédé selon la revendication 2, caractérisé en ce que ladite fonction est de la forme F[z(t).z*(t-Ts)] o z=x+jy représente la valeur du signal sur les voies I et Q et T, est la période de symbole. 4. Procédé selon la revendication 3, caractérisé en ce qu'on obtient le minimum par détermination de l'annulation de la fonction H(t)=G(t+TJ/4)-G(t-T./4), o G représente la moyenne de F[z(t).z*(t-Ts)]. 5. Procédé selon la revendication 1, caractérisé en ce qu'on échantillonne le signal avant conversion en bande de base et en ce qu'on mesure la variation de l'énergie moyenne du signal en fonction de la phase d'échantillonnage, pour sélectionner la phase pour laquelle l'énergie est maximale. 6. Dispositif de démodulation de signal à modulation de phases à M états d'une porteuse, comprenant: des moyens à au moins un oscillateur local (60) de transposition du signal reçu à une fréquence intermédiaire; des moyens de conversion en bande de base sur deux voies (IQ) en quadrature, utilisant un oscillateur local (20) à fréquence ajustable autour d'une valeur nominale; des moyens (28-32) d'échantillonnage et de conversion numérique des deux voies comportant une horloge (32) à phase ajustable, stabilisée à la fréquence de bit du signal entrant; et des moyens (38,40;56,57) de récupération de porteuse et de démodulation cohérente des deux voies, caractérisé par des moyens (36) pour ajuster la phase de l'horloge (32) à la valeur correspondant au minimum d'inter- f érence intersymbole et des moyens pour commander l'oscilla- teur local (60), l'ajustage de la phase d'horloge s'effec- tuant par analyse du signal avant récupération de porteuse. 7. Dispositif selon la revendication 6, caractérisé en ce que les moyens pour ajuster la phase d'horloge comprennent des moyens pour calculer une fonction non linéaire F[z(t).z*(t-Ts)] et des moyens pour commander la fréquence de l'horloge jusqu'à annuler la dite fonction. 8. Dispositif selon la revendication 6, caractérisé en ce que les dits moyens sont constitués pour calculer la répartition de l'énergie en fonction de la phase de l'horloge d'échantillonnage et pour commander la fréquence d'horloge jusqu'à obtention d'un maximum. 9. Dispositif selon la revendication 8, caractérisé en ce que les dits moyens effectuent une recherche de maximum de la fonction de distribution d'énergie P par estimation de deux puissances moyennes, l'une sur des échantillons pairs successifs à la fréquence symbole 1/T$, l'autre sur le même nombre d' échantillons impairs, soustraction des deux puissan- ces estimées et commande de l'horloge d'échantillonnage à partir du résultat de la soustraction. 10. Dispositif selon la revendication 8, destiné à un signal MDP4, caractérisé en ce que les dits moyens sont montés pour recevoir des séquences d'échantillons numériques (Is, Qs) au rythme de deux fois la fréquence symbole 1/Ts et comportent: - une branche de calcul de la puissance moyenne pour les échantillons pairs et une autre branche pour les échantillons impairs, ayant chacune deux éléments d'élévation au carré affectés aux voies I et Q. accumulant plusieurs échantillons successifs, et un additionneur (46), et - un soustracteur (48) formant un signal de sortie représentatif en permanence de la différence, et des moyens de conversion numérique-analogique et de filtrage passe-bas. 11. Disposition selon la revendication 10, caractérisé en ce qu'il comprend de plus des moyens de commande automati- que de gain d'un amplificateur d'entrée (54), comprenant un additionneur (52) des sorties des additionneurs (46) et un convertisseur numériqueanalogique attaquant l'entrée de commande de gain de l'amplificateur.