Brevet FR2807166A1 — publié le 20011005

La présente invention concerne un procédé et un dispositif permettant de corriger la variation ou dérive de la différence de phase A(p existant entre deux signaux générés par un signal de fréquence Fe = 2%Q, une dérive dans la valeur de Fe provoquant la dérive dans la valeur de la différence de phase. L'invention s'applique notamment pour corriger la dérive de la valeur du QDM (code habituellement employé dans le domaine de la navigation pour désigner l'angle e que fait un mobile avec le nord magnétique), afin de satisfaire à certaines normes telles que EUROCAE ED-22B, RTCA, ARINC MARK 2 AIRBONE VOR RECEIVER 711-9, la dérive résultant d'une variation dans la fréquence de rotation. L'invention s'applique par exemple, mais non exclusivement, dans le domaine de la navigation aérienne, la navigation en mer ou fluviale. La figure 1 représente le principe de fonctionnement du système permettant de déterminer à bord d'un mobile ou aéronef, par exemple un avion, la direction dans laquelle il est à un lieu de et vu d'un radiophare. Une onde entretenue est émise suivant un diagramme de rayonnement cardiofdal tournant à 30 tr/s; en même temps, une antenne omnidirectionnelle émet un signal de référence à 30 Hz en modulant en fréquence. Le récepteur de bord détecte la modulation d'amplitude à 30 Hz causée par la rotation du diagramme et le signal de référence. L'azimut du récepteur vu de la station ou balise VOR 1 est égal au déphasage entre les deux signaux E1 et E2; E1 étant la modulation due à la rotation du diagramme et E2 la référence émise. Une variation de la fréquence 30 Hz provoque une dérive linéaire de la différence de phase et le battement introduit au niveau du signal comme le montrent les figures 2 et 3. La figure 4 est un exemple de chaîne de traitement selon l'art antérieur. Une technique de récepteur VOR numérique classique consiste à appliquer un filtre séparateur pour traiter deux bandes de fréquence, la basse fréquence modulée AM et la haute fréquence modulée FM, à appliquer un discriminateur de fréquence pour récupérer la fréquence modulante FM et à comparer les phases des deux signaux. Le signal composite VOR après réception par un récepteur 10 est séparé au moyen d'un séparateur 11 en un signal basse fréquence (BF) modulé AM et un signal haute fréquence (HF) modulé FM. Le signal basse fréquence (BF) est transmis à un filtre linéaire 12 puis échantillonné 13 afin de déterminer sa phase p1 à l'aide d'une transformée de Fourier. Le signal haute fréquence (HF) est envoyé tout d'abord à un filtre non linéaire 14 puis échantillonné 16 afin de déterminer sa phase p2 en mettant en oeuvre une transformée de Fourier. Le filtre non linéaire comporte un discriminateur de fréquence 15 ayant notamment pour fonction de récupérer la fréquence modulante FM. A partir des deux phases (p1 et q)2, la valeur de déphasage Acp correspondant au QDM (au paramètre près des termes de déphasage introduits par les filtres) est déterminée par exemple à l'aide d'un calculateur programmé en conséquence. Au cours du procédé de traitement du signal, I'opération d'échantillonnage des signaux (BF) et (HF) par transformée de Fourier peut engendrer des raies parasites (ou spurious) lorsque la période du signal analysée n'est pas un multiple de la fenêtre temporelle d'analyse. Pour pallier ce problème, il est connu d'utiliser des fenêtres classiques telle que Hamming généralisée permettant de réduire les lobes secondaires du spectre tout en élargissant le lobe principal, ou celle de Kaiser associant à un signal temporel à support fini un spectre infini. Ces deux méthodes présentent toutefois des performances limitées et perturbent le rapport signal/bruit à temps d'intégration égal. Il est aussi possible d'utiliser une méthode à banc de filtre de type DFT. Dans ce cas la puissance de calcul nécessaire est très important et conduit à des coûts d'investissement incompatibles avec certaines applications. La figure 3 représente dans un repère o l'axe des abscisses est l'axe temporel et l'axe des ordonnées la dérive du QDM: > La courbe (I) qui correspond à la valeur QDM(Fe) pour la fréquence de référence Fe ou la période correspondante n, > La courbe (Il) correspond à la dérive du QDM et la courbe (III) à un battement provoqués par une dérive de la valeur de la fréquence de référence Fe. Dans la suite de la description, on désigne par le terme QDM, l'angle pris par rapport au nord magnétique, désigné encore par le terme bearing. L'invention consiste notamment à utiliser une relation reliant la dérive du QDM à une variation de la fréquence de référence et un estimateur de fréquence judicieusement positionné dans la chaîne de traitement pour corriger cette dérive. Elle peut aussi utiliser fenêtre d'échantillonnage de forme trapézoïdale par exemple pour réaliser la transformée de Fourier. L'invention concerne un procédé pour corriger la dérive du déphasage (A(P, o) entre deux signaux (Si) et (S2), lesdits signaux étant générés par un signal ayant une fréquence de référence Fe ou une période correspondante n, la dérive du déphasge résultant d'une dérive en fréquence Fe ou en période selon une relation R((F,Q), (,A ()) donnée. Il est caractérisé en ce qu'il comporte au moins les étapes suivantes: déterminer la valeur du déphasage (Ap,O), à partir de la phase (p1 du signal S. et de la phase qp2 du signal S2, * estimer Fat la valeur de fréquence Fe, ou la période correspondante Q, * ramener ladite valeur (Acp,O) dans un intervalle de valeurs données [QDMmin, QDMmax] contenant la valeur QDM(Fe) obtenue pour la valeur Fe de référence en tenant compte des valeurs A1p et Fest, et de la relation R((F, Q), (0, A)). Selon un mode de réalisation le procédé comporte par exemple une étape de traitement par transformée de Fourier utilisant une fenêtre d'échantillonnage dont la forme est sensiblement trapézoïdale. L'invention concerne aussi un procédé pour corriger une dérive du déphasage (Acp,O) entre deux signaux Si ayant une phase cpl et S2 ayant une phase cp2, les deux signaux étant générés par un signal de fréquence Fe ou de la période correspondante, la dérive résultant de la dérive de Fe ou de la dérive de la période selon une relation donnée R((F,Q),(O,Ap)) Il est caractérisé en ce qu'il comporte au moins une étape de détermination d'au moins une des phases (cpl, 2) par transformée de Fourier utilisant une fenêtre d'échantillonnage de forme sensiblement trapézoïdale. L'invention concerne aussi un dispositif de correction de la dérive du déphasage (A(p,O) existant entre deux signaux (Si) et (S2), lesdits signaux étant générés par un signal ayant une fréquence de référence Fe, la dérive étant liée à la dérive en fréquence Fe selon une relation R((F, Q), (0, A)) donnée. Il est caractérisé en ce qu'il comporte au moins: À un séparateur des signaux S. et S2, des moyens pour déterminer la valeur de la phase (pi et yp2 pour chacun des signaux et un moyen pour déterminer la valeur de déphasage Acp, 25. un estimateur de la fréquence Fe ou de la période correspondante disposé après le séparateur, * un dispositif comportant en mémoire la relation R((F,), (, Aq)) recevant les valeurs A(p et de la fréquence estimée ou de la période correspondante estimée, adapté pour ramener la valeur de Acp dans un intervalle de valeurs données [QDMmin, QDMmax] contenant la valeur QDM(Fe) obtenue pour la valeur Fe (ou la période) de référence. Selon un mode de réalisation il comporte par exemple des moyens de traitement des signaux (Si) et/ou (S2) par transformée de Fourier, la fenêtre d'échantillonnage ayant une forme sensiblement trapézoïdale. L'estimateur pour estimer la valeur de fréquence peut être de type Prony. L'invention concerne aussi un récepteur VOR comportant un dispositif de correction d'une dérive du déphasage (A(p,e8) existant entre deux signaux (S1) et (S2), lesdits signaux étant générés par un signal ayant une fréquence de référence Fe, la dérive étant liée à la dérive en fréquence Fe ou à la dérive en période selon une relation R((F,n),(O,Aq)) donnée caractérisé en ce qu'il comporte au moins: * un séparateur des signaux S. et S2, des moyens pour déterminer la valeur de la phase qp1 et (p2 pour chacun des signaux et un moyen pour déterminer la valeur de déphasage A(p, * un estimateur de la fréquence Fe ou de la période correspondante disposé après le séparateur, À un dispositif comportant en mémoire la relation R((F, Q), (,Aó)) recevant les valeurs A(p et de la fréquence estimée Fm ou de la période correspondante estimée, et adapté pour ramener la valeur de 6tp dans un intervalle de valeurs données [QDMmin, QDMmax] contenant la valeur QDM(Fe) obtenue pour la valeur Fe (ou la période) de référence. Les procédés et les dispositifs selon l'invention s'appliquent par exemple pour réduire la variation du QDM d'un système de navigation ou de radio navigation. Le procédé et le dispositif selon l'invention permettent notamment de respecter ou d'approcher le mieux possible des normes contraignantes, telles que les normes RTCA et EUROCAE, et d'améliorer la limite de détection par rapport aux dispositifs connus de l'art antérieur. L'invention sera mieux comprise à la lecture de la description qui suit à titre illustratif et nullement limitatif d'un exemple donné pour réduire la variation sur le QDM d'un système de radionavigation VOR o: * la figure 1 représente le principe de fonctionnement en VOR, la figure 2 montre dans un repère fréquence-degré(radian), la dérive linéaire du déphasage en fonction d'une dérive de fréquence, 30. la figure 3 représente la dérive et le battement résultant de la dérive de la valeur de fréquence, (diagramme temps-amplitude en degré du déphasage), a la figure 4 schématise un exemple de chaîne de traitement du signal reçu par le VOR selon l'art antérieur, À la figure 5 représente un exemple de chaîne de traitement du signal selon l'invention, * la figure 6 schématise les QDM obtenus sans (art antérieur) ou avec l'estimateur de fréquence, À La figure 7 schématise une fenêtre d'analyse du signal ayant une forme trapézoïdale, La figure 8 représente les spectres exprimés en dB pour différentes fenêtres d'échantillonnage, ò Les figures 9 et 10 montrent l'atténuation des anomalies de phase en utilisant la fenêtre de forme trapézoïdale, * Les figures 11 et 12 la correction du QDM obtenue en mettant en ceuvre le procédé selon l'invention. De manière à mieux faire comprendre le principe mis en oeuvre dans l'invention, la description qui suit donnée à titre illustratif et nullement limitatif concerne la correction de la variation du QDM dans un système de navigation ou de radio navigation suite à une dérive de la fréquence de rotation. Cette manière de procéder permet notamment de répondre aux normes précitées ou de s'en approcher du mieux possible. Sans sortir du cadre de l'invention les différentes variantes de mise en oeuvre du procédé et du dispositif selon l'invention s'apliquent avantageusement dans tout système o une dérive de la fréquence de référence (ou de la période correspondante) entraîne une variation ou déviation d'une différence de phase Ap existant entre deux signaux. La figure 5 schématise un exemple de réalisation pour la chaîne de traitement selon l'invention. La modulation VOR reçue par un aéronef s'exprime par l'expression suivante (1) Ko+k3cos(ot)+kl cos("t -0) +k2 cos(332Qt -msin(nt)) Le signal composite VOR est reçu par un récepteur VOR 20 selon une méthode connue de l'Homme du métier avant d'être transmis à un séparateur 21 de façon à obtenir un premier signal S1 basse fréquence (BF) correspondant à la fréquence de référence Fe par exemple égale à 30 Hz ou à la période correspondante n et un second signal S2 haute fréquence (HF). Le signal (BF) de forme sinusoïdale de phase variable cos(Qt -0) est tout d'abord filtré par un filtre linéaire 22 puis échantillonné 23 à l'aide d'un microprocesseur ou d'un calculateur numérique par mise en oeuvre d'une transformée de Fourier discrète ou DFT pour obtenir la valeur de phase (pl avec par exemple p1=(Qt -8)-L, L étant le retard introduit par le filtre. Une partie de ce signal est dérivée en sortie du filtre linéaire 22 vers un estimateur de fréquence 24 afin de déterminer la valeur estimée Fa de la fréquence de référence Fe ou encore de la valeur estimée es de la période correspondante Q. Les filtres linéaires sont par exemple des filtres demi-bande introduisant un déphasage constant. Les filtres ainsi que le séparateur sont par exemple des dispositifs numériques connus de l'Homme du métier. Le signal (HF) de forme sinusoTdale modulée à 332 fois la fréquence de référence Fe égale à 30 Hz et contenant la phase de référence (cos(332at msin(Qt)) est filtré à travers un filtre non linéaire 25 puis échantillonné au moyen d'une transformée de Fourier 26 afin d'obtenir la valeur de phase (p2 correspondant à la valeur ft au paramètre I près du retard introduit par le discriminateur, avec (p2=.t-l, par exemple. A partir des deux valeurs des phases moyen d'un calculteur. La chaîne de traitement comporte aussi un dispositif 27, qui contient en mémoire la relation R((F, Q), (0, A(p)) liant la dérive de bearing à la dérive de la fréquence de référence ou encore à la dérive de la période correspondante Q. La logique de contrôle 27 reçoit la valeur de la fréquence estimée Fm ou la valeur de la période estimée Qf et la valeur du bearing e ou du déphasage Ap, ce qui permet à partir de la relation R((F,Q),(e,Aó))de ramener la dérive en angle (e,A référence Fe, ou autour d'une valeur donnée, telle que la valeur QDM(Fe). La valeur de l'angle e est calculée par exemple à partir de la valeur du déphasage A(p, en corrigeant cette dernière par une valeur correspondant aux valeurs des retards L et I introduits par les filtres linéaires et non linéaires. Cette correction peut aussi être réalisée après avoir ramené la dérive en angle. Dans cet exemple de réalisation, la relation R((F, ),(0, Ap)) reliant la dérive de la valeur de la fréquence de référence Fe et la dérive de la valeur du QDM, est une relation linéaire telle que décrite à la figure 2. Sur cette figure 2 I'abscisse représente la valeur de la fréquence F variant dans un intervalle [29,5 Hz; 30,5 Hz], autour de la valeur de référence Fe= 30 Hz et en ordonnée, la valeur de la dérive du QDM exprimée en radian ou en degré. Les valeurs limites QDM max et QDM min sont par exemple choisies en fonction de la norme à respecter, par exemple une des normes citées précédemment. La relation R((F, Q), (O, A)) peut être établie au cours d'essais préalables o l'on enregistre pour différentes valeurs de fréquence variant autour d'une valeur de référence Fe donnée, la valeur de déphasage afin d'établir une base de données ((F, Q), (Aoe)). Les transformées de Fourier, le calcul du déphasage ainsi que la correction sur la valeur du déphasage sont réalisées par exemple à l'aide d'un microprocesseur ou calculateur programmé selon une méthode connue de l'Homme du métier. L'estimateur de fréquence 24 est positionné en fin de la chaîne de traitement, notamment après le filtre. Cet estimateur de fréquence est par exemple un estimateur de Prony décrit dans le document IEEE SIGNAL PROCESSING LETTERS VOL.3, N 2, February 1996 " Instantaneous Frequency Estimation Using Linear Prediction with Comparisons to the DESAs ". La figure 6 représente dans un diagramme temps-(degré ou radian), la correction apportée dans la valeur du QDM en mettant en oeuvre les étapes du procédé selon I'invention. La courbe (111) correspond à la valeur de la dérive du QDM obtenue sans utilisation d'un estimateur de fréquence; elle est ramenée (courbe IV) vers la valeur du QDM obtenue pour la valeur de référence Fe représentée par la droite (1). L'exemple numérique qui suit est donné pour un signal de référence Fe sensiblement égal à 30 Hz dans le cadre de l'application VOR. Dans ce cas, la sortie du schéma bloc correspondant au filtre linéaire 22 est une forme sinusoïdale à 30 Hz pour laquelle l'amplitude du spectre des harmoniques indésirables qui la compose peut être atténuée par un filtre Butterworth du 8è"e ordre placé avant l'estimateur de fréquence 24. Le signal VOR est échantillonné à 26 560 Hz et l'entrée du Butterworth est une forme sinusoïdale échantillonnée à 207,5 Hz. Donc, seulement 7 échantillons sont représentatifs d'une période de signal. Parmi ces 7 échantillons, seulement 4 sont utilisés pour estimer la fréquence. La fréquence ou la période sont généralement estimées et une logique de contrôle corrige le QDM lorqu'elle reste stable et dans les standards par exemple (29,55 Hz; 30,45Hz). La précision expérimentale de la valeur de la fréquence estimée est de 104 sans perte d'énergie permettant d'éliminer l'erreur statique. En procédant de cette façon, il est possible de respecter les normes précédemment citées et de rester dans la tolérance de variation pour le QDM de +/- 0,4 degré par rapport à la valeur QDM(Fe) malgré une variation de la fréquence Fe= 30 Hz de +/- 1,5 %. Ceci correspond par exemple à un intervalle de valeurs QDMmin=(QDM(Fe)-0,4 ) et QDMmax=(QDM(Fe)+0,4 ). Selon un mode de mise en oeuvre du procédé, le procédé utilise une fenêtre temporelle de forme trapézoïdale ou sensiblement trapézoïdale pour effectuer la transformée de Fourier du signal, afin de limiter les effets de phase initiale et de phase finale dans l'analyse spectrale. L'amplitude du signal est ainsi réduit de manière linéaire sur la période de fréquence à analyser. Le principe consiste à rapprocher les courbes de fréquence à l'endroit de la discontinuité dans une sorte de courbe impossible. Un exemple de fenêtre de forme trapézoïdale est représenté à la figure 7 dans un diagramme (temps ou nombre d'échantillons)-amplitude normalisée. Elle permet d'atténuer ou de faire disparaître l'apparition de hautes fréquences dues à la périodisation lorsque le signal est différent d'une période multiple de la fenêtre d'analyse. Les bords du signal sont atténués au début et/ou à la fin de la fenêtre d'analyse temporelle. La figure 8 montre le spectre exprimé en dB respectivement d'une fenêtre rectangulaire (RECT), d'une fenêtre triangulaire (TRIA) et d'une fenêtre trapézoïdale (TRAP) selon l'invention. Cette figure montre la faible différence entre le spectre de la fenêtre trapézoïdale et celui de la fenêtre rectangulaire, en terme d'énergie, de largeur de lobe principal et de largeur de lobe secondaire. La fenêtre de type trapézoïdale permet de se rapprocher du gain obtenu par une fenêtre de type rectangulaire. L'abscisse correspond au domaine des fréquences et l'ordonnée à l'amplitude normalisée des lobes exprimé en dB. L'axe vertical correspond au milieu du lobe. La démarche classique consiste à multiplier le signal échantillonné par une fenêtre rectangulaire qui est la fenêtre de référence. Cette approche présente comme inconvénient de faire apparaître des anomalies dans le signal, ou saut de phase référencée 30 sur les figures 9 et 10. L'utilisation d'une fenêtre trapézoïdale, telle que décrite à la figure 7 permet d'atténuer cette anomalie 30 observée comme le montre la figure 10. Sur ces deux figures 9 et 10, l'axe des abscisses correspond au domaine temporel, ou au nombre de points d'échantillonnage à une fréquence d'échantillonnage donnée et l'axe des ordonnées à l'amplitude du signal normalisé. Les paramètres de la fenêtre, par exemple le temps sur lequel la fenêtre présente une atténuation au début et à la fin est par exemple égal à 2,5 fois la période du signal à échantillonner. La forme de type trapézoïdale peut être représentée à l'aide de deux portions de droite Dl et D2 reliées par une doite D3, la pente des droites Dl et D2 étant calculée pour obtenir l'atténuation prémentionnée par exemple. Dans l'application à la correction de la valeur du déphasage ou QDM au niveau d'un récepteur VOR, I'énergie est un paramètre important car il conditionne directement la sensibilité du récepteur. Le choix de la fenêtre d'échantillonnage résulte d'un compromis permettant d'obtenir un bon pouvoir séparateur de la fenêtre, minimiser la perte d'énergie, et un bon rapport signal/bruit. L'utilisation d'une fenêtre de forme trapézoïdale permet en outre d'appliquer un offset de départ sur l'amplitude (sorte de Trapézoïde surélevé) de manière à réduire au minimum utile la perte d'énergie. Cette atténuation pouvant réduire le rapport signal/bruit global du récepteur, l'objectif est de la minimiser afin de conserver le même temps d'intégration, compatible du niveau de réception souhaité d'une part et du temps de réponse de l'équipement d'autre part. Le tableau ci-dessous regroupe les caractéristiques de différentes fenêtres Fenêtre Lobe secondaire (dB) Largeur de Gain Cohérent pente à 3 dB (normalisé) ou (BINS) perte d'énergie Rectangle -13 0,89 1 Triangle 27 1,28 0.5 Hamming -43 1,30 0.54 Hanning (a=1) -23 1,64 0.64 Gaussienne (oa=2.5) -42 1,33 0.51 Kaiser Bessel (a=2.0) -46 1,43 0.49 Tukey (ct=0.25) -14 1,01 0.88 Trapèze -14 1,06 0.80 L'analyse de ce tableau montre que dans la pratique, à part les fenêtres de Tukey, la perte d'énergie est de l'ordre de 3 dB. La plupart de ces fenêtres peuvent s'apparenter à celle de Hamming en modifiant les paramètres de l'équation. Le but est de minimiser les phénomènes de Gibbs pour toutes les fenêtres dérivées de Hanning, d'optimiser le produit Temps Fréquence pour les Gaussiennes,.. . Les figures 11 et 12 schématisent la correction apportée dans la valeur du QDM en appliquant le procédé selon l'invention. La figure 11 représente la valeur du QDM obtenue pour une valeur de fréquence de ,3 Hz en utilisant une fenêtre trapézoïdale telle que décrite à la figure 7 pour réaliser la transformée de Fourier, afin de réduire le battement. La figure 12 montre la valeur du QDM (Vl) obtenue pour une fréquence de 30,3 Hz en utilisant un estimateur de fréquence et une fenêtre trapézoïdale pour la transformée de Fourier. REVENDICATIONS 1 - Procédé pour corriger la dérive du déphasage (A(p,O) entre deux signaux (S1) et (S2), lesdits signaux étant générés par un signal ayant une fréquence de référence Fe ou une période correspondante Q, la dérive du déphasage résultant d'une dérive en fréquence ou en période selon une relation R((F,n),(O,Aq)) donnée caractérisé en ce qu'il comporte au moins les étapes suivantes: * déterminer la valeur du déphasage (Aà,6), à partir de la phase 0(p1 du signal Si et la phase (p2 du signal S2, 10. estimer Ft la valeur de fréquence Fe, ou la période correspondante Q., À ramener ladite valeur (Ap,e) dans un intervalle de valeurs données [QDMmin, QDMmax] contenant la valeur QDM(Fe) obtenue pour la valeur Fe de référence en tenant compte des valeurs Aq et Fm, et de la relation R((F,Q), (0, Aq)) donnée. 2 - Procédé selon la revendication 1 caractérisé en ce qu'il comporte une étape de traitement par transformée de Fourier utilisant une fenêtre d'échantillonnage dont la forme est sensiblement trapézoïdale. 3 - Procédé selon l'une des revendications 1 ou 2 caractérisé en ce que: * la valeur de la fréquence Fe de référence est sensiblement égale à 30 Hz, ou la valeur de la période correspondante sensiblement égale à 21c*30Hz, * le signal SI est un signal BF filtré par un filtre linéaire avant d'être échantillonné par transformée de Fourier, et le signal Si est au moins en partie utilisé pour estimer la fréquence Fe ou la période, Q, * le signal S2 est un signal HF filtré par un filtre non linéaire avant d'être échantillonné par transformée de Fourier. 4 - Procédé selon la revendication 3 caractérisé en ce que la forme trapézoïdale de la fenêtre est choisie pour obtenir à une fréquence 30 Hz ou une période correspondante une atténuation sur 2,5 fois la valeur de la période du signal à échantillonner S. ou S2. - Procédé selon l'une des revendications 1 à 4 caractérisé en ce que la relation R((F, Q), (0, Aç)) est une fonction linéaire. 6 - Procédé selon l'une des revendications 1 à 5 caractérisé en ce que l'intervalle de valeurs tolérables pour l'étape c) est égal à [QDM(Fe)-0,4; QDM(Fe)+O,4 ] pour une valeur de fréquence sensiblement égale à 30 Hz ou une période correspondante sensiblement égale à 27c*30Hz et une variation de +/- 1,5%. 7 - Procédé pour corriger une dérive du déphasage (A(p,O) entre deux signaux S. ayant une phase cp1 et S2 ayant une phase qp2, les deux signaux étant générés par un signal de fréquence Fe ou de période correspondante, la dérive du déphasage résultant de la dérive de Fe ou de la dérive de la période selon une relation donnée R((F,n),(O,Ap)), caractérisé en ce qu'il comporte au moins une étape de détermination d'au moins une des phases (pl, p2) par transformée de Fourier utilisant une fenêtre d'échantillonnage de forme sensiblement trapézoïdale. 8 - Dispositif de correction de la dérive d'un déphasage (A(p,e) existant entre deux signaux (Si) et (S2), lesdits signaux étant générés par un signal ayant une fréquence de référence Fe, la dérive étant liée à la dérive en fréquence Fe selon une relation R((F,),(O,Aq)) donnée caractérisé en ce qu'il comporte au moins: un séparateur (21) des signaux Si et S2, des moyens pour déterminer la valeur de la phase (p1 et cp2 pour chacun des signaux et un moyen pour déterminer la valeur de déphasage Atp, - un estimateur (24) de la fréquence Fe ou de la période correspondante disposé après le séparateur, un dispositif comportant en mémoire la relation R((F,Q),(,AA")) recevant les valeurs Ap et de la fréquence estimée Fw ou de la période correspondante estimée, adapté pour ramener la valeur de Ay dans un intervalle de valeurs données [QDMmin, QDMmax] contenant la valeur QDM(Fe) obtenue pour la valeur Fe ou la période correspondante de référence. 9 - Dispositif selon la revendication 8 caractérisé en ce qu'il comporte des moyens de traitement des signaux (Si) et/ou (S2) par transformée de Fourier, la fenêtre d'échantillonnage ayant une forme sensiblement trapézoïdale. - Dispositif selon l'une des revendications 8 ou 9 caractérisé en ce que le dispositif pour estimer la valeur de fréquence est un estimateur de type Prony. 11 - Récepteur VOR comportant un dispositif de correction d'une dérive d'un déphasage (Ay,O) existant entre deux signaux (Si) et (S2), lesdits signaux étant générés par un signal ayant une fréquence de référence Fe, la dérive étant liée à la dérive en fréquence Fe ou à la dérive en période selon une relation R((F,Q),(0,Aq)) donnée caractérisé en ce qu'il comporte au moins: ò un séparateur (21) des signaux SI et S2, des moyens pour déterminer la valeur de la phase pl et p2 pour chacun des signaux et un moyen pour déterminer la valeur de déphasage Acp, * un estimateur de la fréquence Fe ou de la période correspondante disposé après le séparateur, À un dispositif comportant en mémoire la relation R((F,2),(O, A)) recevant les valeurs Ap et de la fréquence estimée ou de la période correspondante estimée, adapté pour ramener la valeur de A(p dans un intervalle de valeurs données [QDMmin, QDMmax] contenant la valeur QDM(Fe) obtenue pour la valeur Fe (ou la période) de référence. 12 - Application du procédé selon l'une des revendications 1 à 7 ou du dispositif selon l'une des revendications 8 à 11 pour réduire la variation du QDM d'un système de navigation ou de radio navigation.