FR 2960360 A1 20111125 FR 1002152 A 20100521 DISPOSITIF DE DISCRIMINATION DE LA PHASE ET DE LA VARIATION DE PHASE D'UN SIGNAL Le domaine technique de l'invention est celui des dispositifs de démodulation d'un signal dans un système de réception pour lequel une estimation précise de la phase est nécessaire. L'invention s'applique en particulier aux systèmes de radionavigation par satellite pour lesquels l'information contenue dans la phase du signal contribue à déterminer précisément la position d'un récepteur. L'invention concerne un dispositif de discrimination ou estimation de la phase ou erreur de phase lors de la démodulation d'un signal soumis à des variations de phase importantes, par exemple dues à une évolution rapide de la fréquence due à l'effet Doppler. Dans la suite de la description l'utilisation du terme fréquence Doppler est faite en référence à un biais de fréquence due à la vitesse de l'émetteur ou du récepteur, par rapport à la fréquence nominale en absence de mouvement du porteur. Un système de radionavigation par satellite a pour but de déterminer la position d'un récepteur en utilisant le principe de la triangulation. Les pseudo-distances mesurées à partir des signaux reçus par trois satellites sont exploitées pour déterminer la position du récepteur. Outre la précision du positionnement, deux autres paramètres caractérisent les performances d'un récepteur de radio-navigation : la capacité à fonctionner avec un faible rapport signal à bruit et la tenue à la dynamique de phase due aux mouvements du porteur. Dans le cas où le porteur est un avion, les accélérations de mouvement qu'il est susceptible d'engendrer se traduisent par une variation de fréquence Doppler importante sur les signaux reçus par le récepteur. Les moyens de démodulation des signaux reçus par le récepteur utilisent, le plus souvent, le principe connu de boucle à verrouillage de phase pour estimer précisément la phase des signaux reçus afin de déterminer un positionnement précis. Les performances d'une boucle à verrouillage de phase (PLL) sont caractérisées notamment par deux critères, sa tenue aux variations rapides de fréquence Doppler et sa capacité à fonctionner à un faible rapport signal à bruit. Un compromis est nécessaire entre ces deux paramètres, ainsi, l'augmentation de la tenue aux variations de fréquence Doppler se fait au prix d'une résistance moindre au bruit. De plus, elle nécessite une augmentation de la bande passante de la boucle ce qui entraîne un surcoût en charge de calcul pour le récepteur. Le problème visé par la présente invention consiste à augmenter ~o sensiblement la tenue à la dynamique de phase sans augmenter le point de fonctionnement en rapport signal à bruit du système et en limitant la charge de calcul. La tenue à la dynamique des boucles à verrouillage de phase 15 connues est restreinte par le discriminateur de phase qui présente une plage linéaire limitée. La tenue à la dynamique d'une boucle à verrouillage de phase représente sa capacité à tolérer un changement important entre deux états stationnaires. En cas de traînage de ladite boucle, encore appelée boucle de porteuse, pendant la phase transitoire entre les deux états 20 stationnaires, l'erreur de phase peut sortir de la zone de capture du discriminateur ce qui peut entraîner la divergence de la boucle et son décrochage. Afin d'augmenter l'étendue de la zone de capture de la boucle de phase, il est connu de mettre en oeuvre un discriminateur dit «enrouleur» 25 qui enregistre le nombre de tours sur l'erreur de phase ce qui a pour conséquence que la zone de capture devient infinie en erreur de phase. Elle reste cependant limitée en erreur de fréquence ou erreur de vitesse de phase. Cette limitation est directement liée à la fréquence d'échantillonnage utilisée par le dispositif numérique qui exécute la boucle. Dans un récepteur 30 de radionavigation par satellite, la limitation de la fréquence d'échantillonnage de la boucle de phase est liée à la durée d'intégration cohérente appliquée pour démoduler le signal reçu. Cette limite en fréquence peut être pénalisante dans le cas d'applications, par exemple aéroportées, induisant de fortes variations de 5 vitesse de phase. Une autre solution connue est décrite dans la demande de brevet français n° 07 01931 du demandeur qui porte sur un dispositif de réception de signaux satellitaires comprenant une boucle de phase avec compensation des retards. Cette solution permet d'augmenter la bande passante de la 10 boucle sans la rendre instable, ce qui a pour effet d'augmenter la tenue à la dynamique de la boucle sans augmenter la bande de pré-détection du système donc sans alourdir la charge de calcul des boucles. Cependant, cette solution, bien que permettant de réduire le traînage de la boucle de phase et ainsi d'améliorer la tenue à la dynamique ne 15 change pas la limite de fonctionnement en fréquence du discriminateur qui constitue toujours le facteur limitant. Une autre solution connue consiste à utiliser directement une boucle à verrouillage de fréquence mais dans ce cas, la mesure de phase est beaucoup moins précise ce qui est pénalisant pour les performances en 20 terme de précision et d'intégrité de la mesure de positionnement. La présente invention permet de pallier les limitations des solutions connues précitées en utilisant un discriminateur enrouleur étendu en fréquence. 25 A cet effet l'invention a pour objet un dispositif discriminateur de phase recevant en entrée un signal complexe dont l'argument représente une erreur de phase A (p et produisant en sortie une estimée On+1 de l'erreur de phase àcp pour chaque échantillon de signal Zn+1 reçu, caractérisé en ce qu'il comporte des moyens de calculs adaptés à déterminer l'estimée de phase 30 obtenue à un instant (N+1)T1 à l'aide des relations suivantes: 0n+1 = Bn + 0 en+1 + A 9'(n+1 A9'n+1 = DiscriFréquence (n +1) .21r.1; en+1 = arctan(I'n+1 , Q'n+1 ) I'n+1+jQ'n+1 Zn+1-Z,n+1 Z'n+1 = exp(jA9'n+1).Zn où TI est l'intervalle de temps entre deux échantillons Zn et Zn+1 consécutivement reçus, Discrifréquence(n+1) est une estimée de l'erreur de fréquence aocp/at dudit signal complexe produite par un discriminateur de fréquence que comporte ledit dispositif, et Z'n+l est le complexe conjugué de n+1- Dans une variante de réalisation de l'invention, le dispositif discriminateur comporte en outre des premiers moyens d'intégration cohérente dudit signal produisant à partir d'un nombre M d'échantillons complexes Zn,k de signal à la fréquence F'e = 1/To un échantillon complexe Zn à la fréquence Fe = 1/T1, chaque intervalle de temps de durée TI étant subdivisé en un nombre M d'intervalles de temps de durée To, ledit discriminateur de fréquence déterminant l'estimée d'erreur de fréquence Discrifréquence(n) à partir desdits échantillons Zn,k . Dans une variante de réalisation de l'invention, ledit discriminateur de fréquence comporte au moins des premiers moyens de calcul d'au moins deux hypothèses d'erreur de fréquence selon les fréquences -Fo et +Fo avec Fo = Fe/2 : Z+1(n) = Z(n, k)eM k=0 M-1 k Z-1(n) = E Z(n, k)e-,rM k=0 et des seconds moyens de calcul de l'estimée d'erreur de fréquence à l'aide d'une des deux relations suivantes :20 (32CT1 /?r-2) ou Re[(Z+1(n) - .44 (n))-Zo (n)* J (16CT1 /Tr2) avec C la puissance estimée du signal reçu et Zo(n) = Z' Dans une variante de réalisation de l'invention, l'estimée d'erreur de 5 fréquence est remplacée par la moyenne des estimées d'erreur de fréquence obtenues aux instants n et n+1. Dans une variante de réalisation de l'invention, lesdites hypothèses d'erreur de fréquence sont corrigées de l'estimée d'erreur de fréquence obtenue à l'instant précédent : k Z+1(n) =1 Z(n,/(.)e+pz"me j2J - lscrzFréguence~n-1) k=0 M-1 k _ Z_1(n) Z n k e J~Me j2z.DiscriFréquence(n-1) =1(, ) k=0 et que l'estimée d'erreur de fréquence déterminée à l'instant n-1 est ajoutée à celle obtenue à l'instant n. Dans une variante de réalisation de l'invention, lesdits premiers moyens 15 de calcul réalisent un nombre 2S d'hypothèses d'erreur de fréquence aux valeurs de fréquence allant de -SFo à +SFo par pas de Fo selon les relations M -1 k M -1 k Z+p(n) = E Z(n,k)e+'IcpM et Z_p(n) = Z(n,k)e-' M , que ledit k=O k=O discriminateur de fréquence comporte en outre des moyens de détermination = jz+1(n)II2 - IIZ-1(n)II2 DiscriFréquence (n) Discri Fréquence (n) = 10 de l'indice po associé à la valeur IIZPo (n) 2 maximum dans l'ensemble [-(S- 20 1)...,-2,-1,0,1,2,...(S-1)], avec Zo(n)=Z' et que l'estimée d'erreur de fréquence est déterminée à l'aide d'une des deux relations suivantes : { 2 DiscriFréquence (n) 6 Zpo+l (n) Zp°-1 (n)2( + poFo ou (32 CT1 / ïz 2 ) Discri (n) = Re[(Zp°+l(n) - Zp°-1(n))-Zpo (n)*~ +poFo Fréquence (16 CTI / fr 2 ) Dans une variante de réalisation de l'invention, ledit dispositif 5 discriminateur de phase comporte en outre des moyens d'estimation du n rapport signal sur bruit RSB , des moyens de comparaison dudit rapport signal sur bruit avec un seuil prédéterminé en fonction d'un compromis entre l'étendue de la dynamique en fréquence dudit dispositif et sa résistance au bruit et des moyens de désactivation et d'activation dudit discriminateur de 10 fréquence en fonction du résultat de la comparaison. Dans une variante de réalisation de l'invention, ledit dispositif discriminateur de phase comporte en outre des moyens d'estimation de la M-1 2 puissance de bruit ab2, par exemple à l'aide de la relation 6b = EIIZb(n,k)II k=0 avec Zb(n,k) des échantillons de bruit mesurés et le rapport signal sur bruit 15 est estimé à l'aide d'une des relations suivantes : 2 ZP° (n) ou RSB = 1IZ (n)) 2 2 6b 6b 2 A RSB = L'invention a également pour objet un dispositif de réception de signaux de radionavigation par satellite recevant un signal continu modulé par une porteuse, de phase cpreçue, comprenant au moins des moyens 20 permettant de générer une porteuse locale, de phase cp,oca,e, et des moyens de démodulation et de filtrage dudit signal, caractérisé en ce qu'il comporte en outre des seconds moyens d'intégration cohérente dudit signal qui délivrent en sortie un nombre M d'échantillons Zf,k à une fréquence d'échantillonnage F'e et une boucle à verrouillage de phase comprenant au moins un dispositif discriminateur de phase selon l'invention qui reçoit en entrée un signal complexe dont l'argument représente l'erreur de phase A (p = 9iocaie - (Preçue et qui permet d'asservir la phase cp,oca,e de la porteuse locale sur celle (preçue de la porteuse reçue. D'autres caractéristiques apparaîtront à la lecture de la description détaillée donnée à titre d'exemple et non limitative qui suit faite en regard de dessins annexés qui représentent : ~o La figure 1, un schéma de principe des moyens de démodulation mis en oeuvre par un récepteur de radionavigation par satellite, La figure 2, une illustration du principe de fonctionnement d'un discriminateur de phase de type arctangente selon l'art antérieur appliqué à une voie pilote, 15 La figure 3, une illustration du principe de fonctionnement d'un discriminateur de phase de type arctangente selon l'art antérieur appliqué à une voie de données, La figure 4, une illustration du principe de fonctionnement d'un discriminateur de phase « enrouleur » selon l'art antérieur appliqué à une 20 voie pilote, dans le cas où la boucle de phase ne décroche pas, La figure 5, une illustration du principe de fonctionnement d'un discriminateur de phase « enrouleur » selon l'art antérieur appliqué à une voie pilote, dans le cas où la boucle de phase décroche, La figure 6, une illustration du principe de fonctionnement d'un 25 discriminateur de phase « enrouleur » selon l'art antérieur appliqué à une voie de données, La figure 7, une illustration du principe de fonctionnement d'un discriminateur de phase selon l'invention appliqué à une voie pilote, La figure 8, un synoptique comparatif d'un discriminateur de phase 30 « enrouleur » et d'un discriminateur de phase selon l'invention, La figure 9, un schéma de principe de l'évolution de l'erreur de phase sur des sous intervalles d'intégration cohérente, La figure 10, une illustration du principe utilisé pour construire le discriminateur de fréquence selon l'invention, La figure 11, une illustration d'une variante d'un discriminateur de fréquence selon l'invention, La figure 12, une illustration d'une troisième variante d'un discriminateur de fréquence selon l'invention, La figure 13, une illustration du principe de fonctionnement d'un 10 discriminateur de phase selon l'invention appliqué à une voie de données. La figure 1 schématise le principe de la démodulation du signal reçu par un récepteur 10 de radionavigation par satellite. Les traitements situés à gauche de la ligne en pointillés constituent un système continu ou discret 15 avec une fréquence d'échantillonnage élevée Fc, par exemple égale à 100 MHz. Ils sont réalisés préférentiellement par des circuits intégrés de type ASIC ou FPGA compatibles d'une telle fréquence élevée. Les traitements situés à droite de la ligne en pointillés constituent un système échantillonné à une fréquence basse, par exemple égale à 50 Hz, et sont réalisés par un 20 logiciel, par exemple un logiciel de traitement de signal DSP. Le signal reçu 101 est démodulé par une porteuse locale 109 de phase (Plocale(t) par le biais d'un multiplieur complexe 102. Le terme « porteuse locale » désigne un signal généré localement à une seconde fréquence différente de celle du signal reçu, cette seconde fréquence étant 25 commandée grâce à un oscillateur contrôlé numériquement 108. Une seconde démodulation 103 est également appliquée dans le but de démoduler le signal reçu par le code d'étalement local. Des moyens 105 d'intégration cohérente permettent notamment le filtrage du signal reçu 104 démodulé. Le signal complexe Z(n) = I(n) + jQ(n) ainsi obtenu est transmis à 30 un discriminateur de phase 106 qui produit en sortie une estimation en de l'erreur de phase A(Pn = (plocale(n-T) - (preçue(n.T) de la porteuse locale par rapport à la phase de la porteuse du signal reçu. Cette estimation alimente un correcteur de boucle 110. Un tel correcteur 110 peut consister par exemple en un intégrateur d'accélération qui exécute la relation suivante : A(n+1) = A(n) + T.K3.0(n), associé à un intégrateur de vitesse qui exécute la relation suivante : V(n+1) = V(n) + T.A(n) + T.K2.0(n). L'oscillateur 108 contrôlé numériquement peut également constituer un intégrateur de phase exécutant la relation suivante : (piocaie((n+1).T)= (p,ocale(n.T) + T.V(n) + T.K1.0(n) s'il reçoit en entrée une commande en vitesse égale à : V(n) + KI.0(n). Dans cet exemple, la boucle de phase ainsi réalisée est une boucle du 3eme ordre car elle contient trois intégrateurs. L'intégrateur d'accélération permet de supporter une accélération constante sans erreur de traînage en régime permanent. Par contre en régime transitoire, suite à un changement d'accélération, l'erreur de traînage sur la phase peut être importante et le discriminateur 106 doit être capable de suivre l'erreur maximale sans saturer pour que la boucle puisse revenir dans un état stationnaire sans saut de cycle qui serait préjudiciable pour la précision des mesures Le correcteur de boucle 110 délivre une commande en vitesse 107 qui est produite en entrée d'un oscillateur contrôlé numériquement 108 qui délivre à sa sortie la phase de la porteuse locale cp,00a,e(t). Ainsi la boucle à verrouillage de phase mise en oeuvre permet d'asservir la phase cpioca,e(t) de la porteuse locale sur la phase cpreç1e(t) de la porteuse du signal reçu. Les discriminateurs de phase connus sont de type « enrouleur ». Ils présentent une limitation du fait de l'ambiguïté modulo 2 n sur la phase du signal. Lorsque l'erreur de phase a évolué de n radians, soit un demi tour de cercle, entre deux instants de mesure, il subsiste une ambigüité sur le signe de cette erreur, à savoir n ou -n. Par construction le discriminateur 106 sait mesurer une variation d'erreur de phase comprise dans l'intervalle [+n -n[. En présence de variation forte de la phase, le discriminateur mesure une erreur de phase égale à (+ n + c) comme étant égale à (- 7c + E ) ce qui conduit à une erreur de 2 n, susceptible de faire décrocher la boucle. Tant que la variation de l'erreur de phase est strictement limitée à n en valeur absolue entre deux instants de mesure le discriminateur fonctionne linéairement. Au delà la boucle est susceptible de décrocher. La limite d'un tel dispositif discriminateur « enrouleur » en erreur de vitesse de phase (ou en vitesse d'erreur de phase) est de it / T en radians par seconde ou X/T en mètres par seconde, étant la longueur d'onde et T étant la durée entre deux échantillons de corrélation, ou encore 7tFe, avec Fe la fréquence d'échantillonnage du dispositif qui exécute le discriminateur de phase 106. Avec un tel dispositif, la tenue à une erreur de vitesse de phase élevée ne peut se faire qu'en augmentant la fréquence d'échantillonnage Fe ce qui entraine une augmentation de la charge de calcul. Le signal reçu 101 peut être issu d'une voie de données ou d'une voie pilote. Dans le premier cas, le signal 101 est composé d'une porteuse modulée par un code d'étalement et par des données. Dans le second cas, le signal 101 est uniquement composé d'une porteuse et d'un code d'étalement. Le code d'étalement est connu du récepteur, contrairement aux données qui contiennent des informations sur la position des satellites. La figure 2 illustre le principe de fonctionnement d'un discriminateur de phase de type arctangente selon l'art antérieur appliqué à une voie pilote. Sur le haut de la figure sont représentés deux symboles complexes Zn et Zn+1 dans le plan complexe (I,Q) et leurs phases respectives en et en+1 . Les deux symboles correspondent à deux échantillons successifs du signal complexe obtenu en sortie des moyens 105 d'intégration cohérente qui permettent le filtrage du signal reçu après la démodulation par la porteuse locale et le code local. L'estimée de l'erreur de phase en est obtenue en calculant l'arctangente entre la partie réelle et la partie imaginaire du symbole Zn. Sur le bas de la figure 2 est représentée l'évolution de l'estimée de phase en en fonction du temps t=nT=n/Fe. On peut voir que la sortie du discriminateur est limitée entre les valeurs n et -n en radians, ce qui réduit la robustesse de la boucle aux variations de vitesse de phase. Lorsque l'erreur de traînage dépasse n radians, la boucle peut diverger et donc décrocher à cause de la rétro-action en sortie du correcteur qui peut changer de signe. La figure 3 illustre le principe de fonctionnement d'un discriminateur de phase de type arctangente selon l'art antérieur appliqué à une voie de données. Dans ce deuxième exemple, les symboles Zn, Zn+1 obtenus après intégration cohérente 105 sur l'intervalle [nT, (n+1)T] coïncident avec des données binaires [-1, +1] qui modulent le signal. Ces données ne sont pas connues du récepteur, ainsi un changement de signe des données à l'instant (n+1)T se traduira par un changement de signe entre Zn et Zn+1 qu'il ne faut pas confondre avec une variation de phase. La conséquence est que la sortie du discriminateur est cette fois limitée à l'intervalle ]-n/2, n/2] comme illustré sur le bas de la figure 3. Le principe du discriminateur de phase « enrouleur » selon l'art antérieur consiste à utiliser l'erreur de phase mesurée 8n comme référence pour lever l'ambiguïté de l'erreur de phase 8n+1 et supprimer les limitations du discriminateur de type arctangente. On a On+1 = 8n + 08n+1 avec oen+1 = arctangente( In+1', Qn+1' ) ]-n, n] et Z'n+1 = In+1'+ j Qn+1' = Zn+1 . Zn * avec Zn+1 = ( In+1 + j Qn+1 ) ,Zn = (In + J Qn ) et 25 Zn * = (In - j Qn ) la valeur conjuguée de Zn. La figure 4 illustre deux exemples pour lesquels le discriminateur de phase « enrouleur » selon l'art antérieur, appliqué à une voie pilote, fonctionne sans décrochage. Dans ces deux exemples, le saut de phase Den+1 entre deux symboles consécutifs est inférieur à n en valeur absolue. L'évolution de la phase en en fonction du temps est continue et n'engendre aucun décrochage. La figure 5 présente un exemple pour lequel le saut de phase entre les deux symboles consécutifs Zn et Zn+1 est supérieur à it en valeur absolue. Dans ce cas, et comme illustré sur la partie basse de la figure 5, la sortie du discriminateur de phase n'est plus limitée en elle-même. Par contre, une limite persiste sur la variation de l'erreur de phase Ocpn entre deux mesures. Tant que cette variation reste inférieure en valeur absolue à 7c alors on reste dans un fonctionnement linéaire tel qu'illustré sur la partie basse de la figure 4 et la boucle de phase ne décroche pas. Par contre si cette variation dépasse en valeur absolue 7t alors en+1 ne reflète plus l'erreur de phase réelle A(pn+1 = (plocale((n+1).T) - (preçue((n+1).T) et la boucle de phase peut décrocher. La figure 6 illustre le fonctionnement du discriminateur de phase « enrouleur » appliqué à une voie de données pour laquelle la limitation de saut de phase est cette fois égale à 7t/2 en valeur absolue. A cause de l'ambiguïté sur le signe de Zn+1 due à la modulation par les données on se ramène au demi plan de droite limité par la perpendiculaire à l'axe I' qui contient le symbole Zn. Cela revient à limiter Den+1 dans l'intervalle ]-7t/2, 7t/2], ce qui réduit par deux la limite de tolérance en erreur de fréquence Doppler. On a en+1 = en + Den+1, Aen+1 = arctangente( Qn+1' / In+1') ]-ir/2, 7r,/2] avec In+1'+ J Qn+1' = Zn+1 . Zn*, Zn+1 = ( I n + 1 + J Qn+1 ), Zn = (In + J Qn) et Zn * _ (In - j Qn ) la valeur conjuguée de Zn. Le symbole Zn+1 est estimé par un tel discriminateur comme étant le symbole 25 Z' n+1 = - Zn+1 dont l'erreur de phase Den+1 qui le sépare du symbole précédent Zn est inférieure en valeur absolue à 7t/2. La figure 7 illustre le principe d'un discriminateur de phase selon l'invention. Afin de supprimer la limitation sur la variation d'erreur de phase inhérente aux discriminateurs de phase connus, la présente invention se base sur l'ajout d'un discriminateur de fréquence permettant de mesurer, en amont du discriminateur de phase, la vitesse de l'erreur de phase ou erreur de fréquence Doppler. Le principe de l'invention peut être explicité à l'aide des relations suivantes et en concordance avec l'exemple de la figure 7. Deux symboles Zn et Z'+1 consécutifs sont représentés dans le plan complexe (I,Q). Le saut de phase entre ces deux symboles est supérieur à 7c ce qui correspond à un cas de figure pour lequel les discriminateurs de phase connus ne permettent pas d'assurer une convergence correcte de la boucle de phase. L'invention consiste à utiliser un symbole intermédiaire Zn+1 = expj( oen+1 ). Zn construit grâce à un discriminateur de fréquence qui permet de mesurer l'erreur de fréquence Doppler aAcp/at ou erreur de la fréquence de la porteuse locale par rapport à la fréquence de la porteuse reçue. L'erreur de phase oen+1' est obtenue à partir de la relation suivante: oen+1 = DÎSCriFréquence(n+1). 27c . T = DiscriFréquence(n+1) . 2n / Fe OÙ DiscriFréquence(n+1) est une estimation de l'erreur de fréquence Doppler (ou variation de l'erreur de phase) moyenne entre les instants nT et (n+1)T. L'erreur de phase entre les deux symboles consécutifs peut ainsi être décomposée en une première partie Den+1 prévisible grâce à l'estimation de l'erreur de fréquence Doppler et une seconde partie AOn+1 qui correspond à l'erreur de phase non prévisible. Cela se traduit par les relations suivantes : en+1 = en + Den+1 ' + Den+1 (1) avec Den+1 = arctangente( l'"1, Q'n+1 ) E ]-7t, 7c] et In+1'+ ) Qn+1' = Zn+1 . Zn+1' * (2) L'utilisation d'un discriminateur de fréquence permet de pré corriger le signal afin de ramener l'erreur de phase résiduelle Den+1 dans un intervalle ]- n, n] compatible du fonctionnement d'un discriminateur de phase « enrouleur ». La figure 8 schématise à des fins comparatives un synoptique d'un discriminateur de phase 91 « enrouleur » et un synoptique d'un discriminateur de phase 92 selon l'invention. Ce dernier utilise un discriminateur de fréquence 921 appliqué au signal complexe Z(n,k) reçu avant intégration cohérente à une fréquence d'échantillonnage plus élevée. Le fait d'appliquer le discriminateur de fréquence sur le signal d'entrée des moyens 105 d'intégration cohérente permet d'estimer une fréquence Doppler plus élevée que s'il était appliqué sur le signal de sortie desdits moyens 105 qui présentent une fréquence d'échantillonnage basse due au sous-échantillonnage effectué par les moyens d'intégration cohérente 105. Les moyens 105 peuvent également être séparés en deux étages successifs comme cela sera explicité plus loin dans la description. Un premier étage effectue une intégration cohérente sur une durée To par exemple une durée de 1 ms, dont le résultat est utilisé par le discriminateur de fréquence 921 puis un second étage finalise l'intégration cohérente sur une durée TI , par exemple une durée de 20 ms qui correspond à la période d'échantillonnage du logiciel de traitement de signal qui exécute le discriminateur de phase selon l'invention. Le discriminateur de fréquence 921 produit en sortie une estimée de variation d'erreur de phase prévisible NO qui est exploitée conformément aux relations (1) et (2) pour produire les signaux et Q'n à partir desquels est calculée la variation d'erreur de phase non prévisible AO . Cette dernière est ajoutée à la variation derreur de phase prévisible NO délivrée en sortie du discriminateur de fréquence 921 afin de produire une estimée de variation d'erreur de phase globale précise qui est ensuite produite en entrée d'un intégrateur 922.30 Le principe utilisé pour la construction du discriminateur de fréquence selon l'invention est explicité à présent. Les moyens 105 d'intégration cohérente réalisent notamment une intégration cohérente du signal reçu démodulé sur un intervalle de temps de durée TI. Cet intervalle est découpé en un nombre M de sous-intervalles d'intégration élémentaire de durée To Les moyens d'intégration 105 fournissent au discriminateur de phase selon l'invention les résultats des intégrations cohérentes sur une durée To. L'intégration finale sur la durée TI est effectuée par la partie logicielle du système de démodulation décrit à la figure 1. En parallèle deux filtres gauche et droite qui réalisent deux hypothèses d'erreur de fréquence Doppler symétriques sont appliqués aux résultats des intégrations cohérentes sur une durée To pour construire un discriminateur de fréquence permettant d'extrapoler la variation d'erreur de phase. La figure 9 schématise l'évolution de la phase du signal reçu (Preçue(t) d'une part et de la phase locale (pioca,e(t) d'autre part. Fe=1/Ti est la fréquence d'échantillonnage, ou de traitement, du dispositif discriminateur de phase selon l'invention. Le symbole Z(n) de phase 0(n) est étalé sur tout l'intervalle de durée TI. Sa phase 0(n) correspond à la phase moyenne sur cet intervalle de temps. L'intervalle [(n-1)T1 nTI] est découpé en un nombre M de sous- intervalles I(n,k) de durée To pour lesquels un symbole Z(n,k) est produit par les moyens d'intégration 105 qui réalisent ainsi une intégration par morceaux. Une première hypothèse d'erreur Doppler correspondant à la fréquence Fo = 25 Fe/2 = 1/(2T1) est réalisée à partir des symboles Z(n,k): Z+1(n) = E Z (n, k)e+'7cM k=0 Une seconde hypothèse d'erreur Doppler correspondant à la fréquence -Fo est également réalisée : M-1 k Z-1(n) = Z(n, k)e 'IrM k=0 Le calcul de la quantité IIZ+1(n)II2 -IIZ_1(n)II2 permet d'obtenir une estimée d'erreur de fréquence proportionnelle à l'erreur de fréquence Doppler si celle-5 ci reste limitée à la zone linéaire entre -Fo et +Fo tel que cela est illustré sur la figure 10. Le discriminateur de fréquence selon l'invention produit en sortie l'estimée suivante exprimée en Hertz : Discri Fréquence (n) (32 CT1 Ir2 ) 10 avec C la puissance estimée du signal reçu par exemple par l'intermédiaire de la relation suivante C = IIZ(n)IIZ. Dans une variante de réalisation de l'invention, une troisième hypothèse d'erreur Doppler, dite centrée, est exploitée. Cette hypothèse est déterminée par la relation suivante : M-1 15 Zo (n) = Z(n, k) = Z(n) k=0 Le discriminateur de fréquence selon l'invention produit dans ce cas en sortie l'estimée suivante exprimée en Hertz : n) - Re[(Z+1(n) - Z_1(n)).Zo (n)* ] DlSCYiFréquence ( - (16C1 / .2) 20 avec Zo(n)* , le complexe conjugué de Zo(n) . La figure 11 illustre une variante de réalisation du discriminateur de fréquence selon l'invention. A l'instant (n+1)T1, les estimées Discrifréque'ce(n) (n)jr -Vz-. (n)ll' } et Discrifréquence(n+1) des variations de phase aux instants n et n+1 sont disponibles. Discrifréquence(n) est représentatif de la variation de l'erreur de phase sur l'intervalle de temps [(n-1)T1 nTd et Discrifréquence(n+1) est représentatif de la variation de l'erreur de phase sur l'intervalle de temps [nT1 (n+1)T1]. Ces deux estimées sont représentées sur la figure 11 respectivement par les pentes des droites que suivent l'erreur de phase &p en fonction du temps sur chaque intervalle de durée T1. Afin d'améliorer la précision du discriminateur de fréquence, la moyenne de deux estimations successives est utilisée afin de produire en sortie le résultat suivant : Discri Fréquence (n) + Discri Fréquence (n + 1) 2 Dans une autre variante de réalisation de l'invention, la sortie du discriminateur de fréquence calculée à l'instant (n-1)Ti est utilisée pour corriger les hypothèses d'erreur Doppler gauche et droite effectuées initialement aux fréquences Fo et -Fo. Les relations suivantes sont utilisées : k Z+1 (n) _ Z(n,k)e+JrMe 2z.Discr! eyue.(n-1) k=0 z_1(n) k=0 Re[(Z+1(n) - Z-~ (n)) -Zo (n)* (16CT1 /, 2) + Discri Fréquence (n -1) Dans une autre variante de réalisation de l'invention, tel qu'illustré à la 20 figure 12, afin d'élargir encore le domaine de fonctionnement en erreur de trainage en vitesse du discriminateur de phase, on effectue un nombre S d'hypothèses d'erreur Doppler aux valeurs de fréquence suivantes : -S Fo,..., - 2 Fo,- Fo, 0, Fo, 2F0, ... SFo. 15 Discri Fréquence (n) = Ces N hypothèses s'écrivent de façon générale à l'aide de la relation suivante : M-1 -fin. p k Zp(n) = 1Z(n,k)e pour p variant de -S à s k=0 L'indice po correspondant à la valeur Iz maximale parmi toutes les hypothèses calculées est retenu afin de construire le discriminateur de fréquence autour de ce maximum à l'aide d'une des deux relations suivantes : 2 IZpo+1 (n) Izpo-1 (n)112 Discri Fréquence (n) 2 + poFo 02CT1/rr ) DZSCYZ n = Re[(Zpo+1(n) - Zpo-1(n))-Zpo (n)*] + F Fréquence () (1 6 CTl / ?l2) po 0 Dans une autre variante de réalisation de l'invention, l'utilisation du discriminateur de fréquence est conditionnée au rapport signal à bruit estimé. Le discriminateur de phase étendu en fréquence selon l'invention permet 15 d'augmenter sensiblement la zone de fonctionnement en phase de la boucle à verrouillage de phase. Cependant l'utilisation du discriminateur de fréquence est susceptible de rendre la boucle de phase moins robuste au bruit thermique ou au bruit dû aux interférences. Pour pallier cet inconvénient, la sortie du discriminateur de fréquence est forcée à zéro 20 lorsque le rapport signal à bruit estimé devient inférieur à un seuil donné, par exemple un seuil de 10 dB. L'estimation du rapport signal à bruit est effectuée à partir des hypothèses d'erreur Doppler préalablement calculées à partir d'une des relations suivantes :10 A RSB = 11ZP0 (n) I2 ou RSB = IIZo (n)II2 2 2 6b Où 6b2 est une estimation de bruit réalisée par le récepteur sur l'intervalle d'intégration de durée TI, par exemple à l'aide de la relation M-1 2 6b = E IIZb (n, k)II , avec Zb (n, k) sont des échantillons de bruit calculés, k=0 sur l'intervalle de temps de durée TI, sur une voie de démodulation dédiée ne contenant que du bruit en sortie. La figure 13 illustre le fonctionnement du discriminateur de phase selon l'invention appliqué à une voie de données. Le principe utilisé est identique, à la différence près que l'erreur de phase 06'+1 présente une ambigüité de n/2 au lieu de n dans le cas d'une voie pilote. La tolérance à la dynamique en erreur de phase est donc moindre et la vitesse d'erreur de phase maximum pour laquelle le discriminateur de phase selon l'invention fonctionnera correctement est sensiblement deux fois moins importante que dans le cas d'une voie pilote. Dans une variante de réalisation de l'invention, les seconds moyens d'intégration cohérente 105 sont supprimés et le dispositif discriminateur de fréquence 921 réalise une estimation de l'erreur de fréquence directement à partir du signal reçu à la fréquence Fe. L'estimée d'erreur de fréquence peut également être réalisée par un dispositif externe à l'invention et est dans ce cas fournie comme paramètre d'entrée au discriminateur de phase selon l'invention. La présente invention s'applique à tous les récepteurs de radionavigation, et plus largement à tout récepteur incluant des moyens de démodulation d'une porteuse dont on souhaite estimer précisément la phase ou l'erreur de phase. Elle permet soit d'augmenter la tenue à la dynamique 6b en phase en limitant la charge de calcul, soit de réduire la fréquence d'échantillonnage de la boucle de phase pour des performances équivalentes. Dans le cas d'un récepteur adapté à démoduler un signal dont le spectre n'a pas été étalé, les moyens 105 d'intégration cohérente se confondent avec des moyens de filtrage et de sous échantillonnage. Dans le cas d'un récepteur de radionavigation, les moyens 105 d'intégration cohérente réalisent implicitement un sous échantillonnage du signal. REVENDICATIONS1. Dispositif discriminateur (106) de phase recevant en entrée un signal complexe dont l'argument représente une erreur de phase A (p et produisant en sortie une estimée On+1 de l'erreur de phase Aç pour chaque échantillon de signal Zn+1 reçu, caractérisé en ce qu'il comporte des moyens de calculs adaptés à déterminer l'estimée de phase obtenue à un instant (N+1)T1 à l'aide des relations suivantes: 0n+1 = en + Aen+1 + Ae'n+1 Ae'n+1 = DiscriFréquence (n + 1).2,r.T1 A en+1 = arctan(I'n+1 ' Q'n+1 ) 'n+1 +JQ'n+1 =Zn+1-Z'n+1 Z'n+1= exp(jA e'n+1).Zn où TI est l'intervalle de temps entre deux échantillons Zn et Zn+1 consécutivement reçus, Discrifréquence(n+1) est une estimée de l'erreur de fréquence aicp/at dudit signal complexe produite par un discriminateur de fréquence (921) que comporte ledit dispositif (106), et Zen+1 est le complexe conjugué de Z'n+i . 2. Dispositif discriminateur (106) de phase selon la revendication 1 caractérisé en ce qu'il comporte en outre des premiers moyens d'intégration cohérente (105) dudit signal produisant à partir d'un nombre M d'échantillons complexes Zn,k de signal à la fréquence F'e = 1/To un échantillon complexe Zn à la fréquence Fe = 1/TI, chaque intervalle de temps de durée TI étant subdivisé en un nombre M d'intervalles de temps de durée To, ledit discriminateur de fréquence (921) déterminant l'estimée d'erreur de fréquence Discrifréquence(n) à partir desdits échantillons Zn,k . 3. Dispositif discriminateur (106) de phase selon la revendication 2 caractérisé en ce que ledit discriminateur de fréquence (921) comporte au moins des premiers moyens de calcul d'au moins deux hypothèses d'erreur de fréquence selon les fréquences -Fo et +Fo avec Fo = Fe/2 : M-1 fer k + Z+1(n) = Z(n,k)e m k=0 k Z-i(n) = 1Z(n,k)e m , k=0 et des seconds moyens de calcul de l'estimée d'erreur de fréquence à l'aide d'une des deux relations suivantes : Z+~ (n)II2 - I Z-i (n)II2 } (32CT1 /9z2) ou (n) = Re[(Z+1(n) - Z_1(n))-Zo (n)*] DZSCYZ Fréquence - 2 (16CT1/7z ) avec C la puissance estimée du signal reçu et Zo(n) = Z' 4. Dispositif discriminateur (106) de phase selon la revendication 3 caractérisé en ce que l'estimée d'erreur de fréquence est remplacée par 15 la moyenne des estimées d'erreur de fréquence obtenues aux instants n et n+1. 5. Dispositif discriminateur (106) de phase selon la revendication 3 caractérisé en ce que lesdites hypothèses d'erreur de fréquence sont 20 corrigées de l'estimée d'erreur de fréquence obtenue à l'instant précédent : Z+1(n) = IZ(n,/()e+JIrme J27L.DZSCYZFréquence(n-1) k=0 10 DiscriFréquence (n) = M-1 kM-1 k _ / Z_1(n) Z n k e J2rMe J22r.DiscriFréquence(n-1) =1(,) k=0 et que l'estimée d'erreur de fréquence déterminée à l'instant n-1 est ajoutée à celle obtenue à l'instant n. 6. Dispositif discriminateur (106) de phase selon la revendication 3 caractérisé en ce que lesdits premiers moyens de calcul réalisent un nombre 2S d'hypothèses d'erreur de fréquence aux valeurs de fréquence allant de -SFo à +SFo par pas de Fo selon les relations Z+p(n) = E'Z(n,k)e+~npM et Z_p(n) _ E`Z(n,k)e-,7rpM , que k=0 k=0 ledit discriminateur de fréquence (921) comporte en outre des moyens de détermination de l'indice po associé à la valeur Imaximum dans l'ensemble [-(S-1)...,-2,-1,0,1,2,...(S-1)], avec Zo(n)=Z' et que l'estimée d'erreur de fréquence est déterminée à l'aide d'une des deux relations suivantes : Zpo+1 (n) 2 - Zpo-1 (n) 2 Discri (n) I l + F Fréquence () (32 CT1 / 1L2) po 0 ou (n) _ Re[(Zpa+1 (n) - Zpo-1(n))-Zpo (n)*] + poFo (16CT1/e2) 7. Dispositif discriminateur (106) de phase selon l'une des revendications 2 20 à 6 caractérisé en ce qu'il comporte en outre des moyens d'estimation du rapport signal sur bruit RSB , des moyens de comparaison dudit rapport signal sur bruit avec un seuil prédéterminé en fonction d'un compromis entre l'étendue de la dynamique en fréquence dudit dispositif et sa Discri Fréquencerésistance au bruit et des moyens de désactivation et d'activation dudit discriminateur de fréquence (921) en fonction du résultat de la comparaison. 8. Dispositif discriminateur (106) de phase selon la revendication 7 caractérisé en ce qu'il comporte en outre des moyens d'estimation de la puissance de bruit ab2, par exemple à l'aide de la relation M-1 2 2 6b = Zb(n,k) avec Zb(n,k) des échantillons de bruit mesurés et k=0 que le rapport signal sur bruit est estimé à l'aide d'une des relations suivantes : RSB = Zpo (n) 2 6b 2 2 ou RSB = IZ°(n)II 6b 2 9. Dispositif (10) de réception de signaux de radionavigation par satellite recevant un signal (101) continu modulé par une porteuse, de phase (Preçue, comprenant au moins des moyens permettant de générer une porteuse locale, de phase (p,oca,ei et des moyens de démodulation et de filtrage dudit signal (101), caractérisé en ce qu'il comporte en outre des seconds moyens (105) d'intégration cohérente dudit signal qui délivrent en sortie un nombre M d'échantillons Zn,k à une fréquence d'échantillonnage F'e et une boucle à verrouillage de phase comprenant au moins un dispositif discriminateur (106) de phase selon l'une des revendications précédentes qui reçoit en entrée un signal complexe dont l'argument représente l'erreur de phase Ocp = (Plocale - (Preçue et qui permet d'asservir la phase (p,oca,e de la porteuse locale sur celle (Preçue de la porteuse reçue.25