Brevet FR2894098A1 — publié le 20070601

Procédé et dispositif permettant le suivi de Doppler pour modem à large bande L'invention concerne un procédé et un dispositif permettant le suivi Doppler, pour modem à large bande. Plus généralement, elle concerne aussi un système d'estimation et de suivi du décalage en fréquence pour un système de transmission dont la largeur de 5 bande n'est pas négligeable devant la fréquence centrale. Le problème posé est la transmission de données de type quelconque entre deux correspondants immergés dans l'eau et à un débit de plusieurs centaines de bits par seconde tout en utilisant une bande passante aussi faible que 10 possible et une puissance de calcul réduite au minimum, en particulier pour des équipements autonomes fonctionnant sur batteries. Le problème est aussi celui de la transmission de données entre deux correspondants immergés dans un milieu donné. Les ondes électromagnétiques classiquement utilisées (de quelques 15 centaines de kHz à quelques GHz se propageant extrêmement mal dans l'eau, le modem (modulateur-démodulateur) utilise donc des ondes acoustiques pour véhiculer les messages transmis. Pour rappel, le milieu sous-marin a un comportement extrêmement complexe pour les ondes acoustiques. Ses inhomogénéités (salinité, température, 20 courants, ..) font que le chemin suivi par les ondes n'est une ligne droite qu'à très courte distance et que ce trajet varie au cours du temps. Les ondes qui arrivent au récepteur peuvent suivre plusieurs trajets avec des délais très différents. Ces trajets peuvent être soit directs (plusieurs parcours différents convergeant vers le récepteur), soit indirects (réflexions sur 25 la surface de l'eau, éventuellement agitée de vagues, ainsi que sur le fond qui peut avoir toutes sortes de types de reliefs). Dans la pratique, cela se traduit par l'arrivée au récepteur de nombreuses répliques du signal émis, chacune ayant à un instant donné, une amplitude, une phase, un retard et un décalage fréquentiel propre. A ceci, il faut ajouter, éventuellement, un décalage fréquentiel du à la vitesse radiale relative des correspondants (effet Doppler) ainsi qu'un décalage fréquentiel pseudo périodique dû à la houle dans le cas où l'un au moins des correspondants est à la surface. Dans certains cas, la fiabilité de la transmission doit être quasi parfaite. S'il s'agit d'un texte, on peut tolérer quelques erreurs isolées, car elles seront automatiquement corrigées par le lecteur. Lorsqu'il s'agit d'un fichier binaire (par exemple une image JPEG), le taux d'erreurs exigé est voisin de 0, sous peine de très forte dégradation de l'image. De plus, le système est à large bande . Dans les applications pratiques, la largeur de bande du signal émis n'est pas négligeable comparée à la fréquence centrale d'émission. Le décalage Doppler étant proportionnel à la fréquence, les fréquences basses ont un décalage en fréquence plus faible que les fréquences hautes. Qui plus est, la vitesse relative de l'émetteur et du récepteur n'étant pas toujours négligeable par rapport à la vitesse du son dans l'eau (environ 1500 m/s) on observe un effet de dilatation ou de contraction du temps selon le signe de la vitesse relative. Par exemple, un message émis pendant une durée de 10s sera reçu 20 pendant 10.13 secondes si la vitesse d'éloignement est de +10 m/s (36 km/h ou environ 20 noeuds) ce qui pose des problèmes de suivi de synchronisation. A la connaissance du demandeur, il existe des systèmes de type à étalement de spectre dans lesquels le débit utile est très faible comparé à la 25 bande passante utilisée et où la discrétion est de mise, ce qui interdit l'utilisation de porteuses fixes, facilement repérables. La demande de brevet WO 2005/055442 divulgue un système à large bande qui estime une seule fois et pour toute la durée du message une valeur de décalage Doppler sans effectuer de suivi. 30 L'objet de la présente invention concerne notamment un procédé permettant d'estimer en continu le décalage en fréquence (Doppler) pour un modem fonctionnant en milieu sous-marin. L'invention concerne notamment un système d'estimation et de suivi du décalage en fréquence pour un système de transmission dont la largeur de bande n'est pas négligeable devant la fréquence centrale, le système émettant une forme d'onde incluant des porteuses fixes de niveau sensiblement supérieur à celui des porteuses véhiculant l'information, les porteuses ayant des fréquences F0 • • • FN-1 et le signal utile porteur de l'information caractérisé en ce qu'il comporte au moins les éléments suivants : • N bancs de filtres passe-bande, chaque banc comportant P filtres, adaptés à couvrir l'ensemble des valeurs que peuvent prendre les fréquences F0 ... 10 FN-1 après décalage Doppler, • La sortie d'un filtre est reliée à un détecteur quadratique d'enveloppe, • P sommateurs adaptés à additionner les sorties correspondantes des N bancs de filtres, • Un système de détection de la position des deux plus grands maxima parmi 15 les P sommes et de validation de la position du plus grand maxima, • Un système de mémorisation des positions des deux plus grands maximas et de leurs poids et de prédiction de la valeur courante, • Un système de validation/remplacement de la position trouvée, • un système convertissant les positions en Doppler estimé et filtrant celui-ci 20 par un filtre passe-bas. L'invention concerne aussi un procédé permettant d'estimer et de suivre l'effet Doppler pour un signal composé d'une forme d'onde incluant des porteuses fixes de niveau sensiblement supérieur à celui des porteuses véhiculant l'information, les porteuses ayant des fréquences Fp ... FN_1 et du signal utile 25 porteur de l'information dans un système de transmission dont la largeur de bande n'est pas négligeable devant la fréquence centrale caractérisé en ce qu'il comporte au moins les étapes suivantes : 1. Transmettre le signal composite comportant le signal utile porteur de l'information et le signal de référence à une étape de filtrage utilisant 30 les fréquences des porteuses, 2. Déterminer pour chaque filtre, l'enveloppe du signal, en,p, 3. Sommer les différentes sorties en,p, 4. Déterminer les valeurs des deux plus grands maximas, M et m, pour plusieurs instants, t, 5. Déterminer la valeur probable de la position à l'instant t, en 5 considérant les maximas trouvés qui définissent de façon unique la valeur de décalage Doppler associée, 6. Effectuer une filtrage passe-bas sur les valeurs de décalage Doppler pour obtenir une valeur lissée. Le procédé comporte par exemple une étape où l'on interpole la valeur 10 lissée à des instants intermédiaires au moyen de filtres passe-bas dérivés du filtre défini précédemment. L'étape de filtrage 1 comporte par exemple au moins les étapes suivantes : • On effectue le fenêtrage du signal d'entrée par une fenêtre de Hanning 15 (simple cosinus sur-élevé) sur une durée de 160 ms, avec recouvrement de 80 ms entre deux fenêtres successives, • aux 1.024 points obtenus, on ajoute 3.072 points nuls, puis : • on effectue la transformée de Fourier rapide FFT 4096 points sur ce signal dont la sortie est notée Zo, Z1 • • • Z4095 20 • pour une fréquence centrale Fn,p souhaitée, on calcule le rang n du "bin" Zn de la transformée de Fourier FFT qui est juste inférieur, • on calcule le nombre réel x compris entre 0 et 1 qui vaut (Fn,p - 1.5625 n) / 1.5625 et l'on sélectionne deux nombres complexes zO(x) et z1 (x) fonctions de x, 25 • on détermine la sortie du filtre par zO(x) Zn + z1 (x) Zn+1. Le procédé et le système selon l'invention présentent notamment les avantages suivants. Ils permettent d'effectuer des communications sous-marines de longue durée, à Doppler variable, entre un émetteur et un récepteur en déplacement 30 sur des trajectoires quelconques. D'autres caractéristiques et avantages de la présente invention seront mieux compris à la lecture de la description qui suit d'un exemple de réalisation donné à titre illustratif et nullement limitatif, annexé des figures qui représentent : • La figure 1, dans un espace temps-fréquence, l'effet du décalage sur les porteuses, • La figure 2, un exemple de structure des modules et des étapes mises en oeuvre par le système selon l'invention, • La figure 3, une étape de filtrage du signal par fenêtrage du signal d'entrée, • La figure 4, un exemple de sorties des détecteurs, • La figure un exemple des sorties des sommateurs, • La figure 6, un exemple de régression parabolique, • Les figures 7 à 9, les résultats du filtrage final, et • La figure 10, un exemple de réponse en fréquence obtenue. Afin de mieux faire comprendre le principe mis en oeuvre dans l'invention, l'exemple qui suit est donné pour un modem fonctionnant en milieu sous-marin, caractérisé par un étalement temporel important et rapidement variable et une largeur de bande non négligeable devant la bande passante qui rend le décalage en fréquence du à l'effet Doppler variable dans la bande passante même du signal (effet bande large ). La forme d'onde émise par le modem en question inclut des porteuses fixes (non modulées) de niveau sensiblement supérieur à celui des porteuses véhiculant l'information. C'est à partir de ces porteuses que se fait l'estimation du décalage Doppler moyen du signal, décalage susceptible de varier rapidement. Les porteuses de référence sont en nombre égal à N et ont des fréquences notées Fo, F1, ... FN_1. Le décalage Doppler maximum est donné par un nombre sans dimension, noté S, obtenu par l'équation: S vr(R) vr(E) C -vr(R) C : vitesse du son dans l'eau v r (E) : vitesse radiale de l'émetteur v r (R) : vitesse radiale du récepteur Référentiel : l'eau, supposée immobile En présence d'un tel décalage Doppler, chaque porteuse de fréquence F; voit sa fréquence modifiée en (1 + 8) F. La figure 1 décrit dans un espace temps-fréquence, l'effet du décalage Doppler sur les porteuses, pour 8 = -8max, d = 0, et d = +8max, 8max étant un élément 5 de conception du système. Les fréquences choisies peuvent être quelconques, mais, afin que l'estimation de 8 puisse se faire sans ambiguïté, il faut que les petits rectangles 1 de la figure 1 aient une hauteur positive, ce qui se traduit par la condition simple suivante: (1 - 8max) Fn > (1 + cSmax) Fi,-ou Fn- Fn_1 >8max (Fn+Fn-1 ) n = 1 ... N-1 10 La figure 2 représente les différentes étapes mises en oeuvre dans l'invention et les éléments constituant le dispositif. Ce système nécessite la présence, pendant toute la durée du message 15 et même un peu plus, d'un ensemble de N porteuses fixes non modulées en sus du signal utile. Le signal complexe (porteuses de référence et signal utile) est reçu sur une chaîne de traitement comportant, dans cet ordre: - un ensemble, 100, de N bancs de P filtres passe-bande couvrant 20 (légèrement plus que) au moins l'ensemble des valeurs que peuvent prendrei les fréquences Fo ... FN_1 après décalage Doppler. Chaque filtre passe-bandes est suivi d'un détecteur quadratique d'enveloppe non représenté distinctement des filtres sur la figure pour des raisons de simplification. 25 - un ensemble, 200, de P sommateurs additionnant les sorties des détecteurs quadratiques d'enveloppe de même rang pour chaque filtre, un système, 300, de détection de la position des deux plus grands maxima parmi les P sommes et de validation de leurs positions en sortie de la série de sommateurs, un système 400 de mémorisation des positions et de leurs poids et de prédiction de la valeur courante, un système 500 de validation / remplacement de la position trouvée, un système 600 convertissant les positions en Doppler estimé et filtrant celui-ci par filtre bas qui donne une estimation du décalage Doppler (avec un retard fixe). Le mécanisme d'estimation et de suivi peut aussi inclure un mécanisme d'interpolation permettant de disposer plus fréquemment d'échantillons de décalage Doppler. Le détail des différents éléments faisant partie du système selon 15 l'invention est donné ci-après. Bancs de filtres 100 Un banc de filtres est adapté notamment à couvrir les fréquences comprises entre (1 - Smax) Fn et (1 + Smax) Fn où Fä est la fréquence après décalage Doppler pour le banc de filtre de rang n. 20 Chaque banc de filtres comporte, par exemple, le même nombre P de filtres passe-bande de fréquences centrales régulièrement espacées. Le pème filtre du nème banc sera centré sur une fréquence égale à : = Fn (1+ Smax (2p-P+1 P I n=0...N-1 p = 0 ... P-1 25 Fn : la fréquence après décalage Doppler pour le banc de filtre de rang n. Une manière simple de réaliser le banc de filtres est décrite ci-dessous. 10 Fn, p A titre d'exemple, on suppose que les porteuses sont espacées d'un multiple de 25 Hz et que la fréquence d'échantillonnage est de 6.400 Hz. On suppose que, lors de la phase de démodulation, un nouveau symbole apparaît toutes les 80 ms et que l'on a alors besoin d'une nouvelle valeur 5 du décalage en fréquence. En supposant que la démodulation de phase différentielle peut avoir jusqu'à 4 états de phase, soit au minimum 90 . Si l'erreur d'estimation du décalage en fréquence vaut dF, la rotation de phase supplémentaire qu'il introduit est de 360 x 80 ms x dF. 10 Pour qu'elle ne dépasse pas le quart de la rotation minimum (règle empirique issue de l'expérience), l'erreur d'estimation du décalage en fréquence dF doit valoir au maximum dfMAx = 0.78125 Hz (en valeur absolue), ce qui signifie que la transformée de Fourier rapide, FFT, à utiliser doit pouvoir séparer deux fréquences séparées de 2 dfmAx = 1.5625 Hz. Son nombre de points devrait alors 15 être d'au moins 6.400/1.5625 = 4.096. Pour cela, on procède de la manière suivante : • On effectue le fenêtrage (pondération) du signal d'entrée par une fenêtre de Hanning (simple cosinus sur-élevé) sur une durée de 160 ms, avec recouvrement de 80 ms entre deux fenêtres successives. Cette étape est 20 illustrée à la figure 3. • aux 0.16 x 6.400 = 1.024 points obtenus, on ajoute 4.096 ù 1.024 = 3.072 points nuls, puis : 1. on effectue la transformée de Fourier rapide FFT 4096 points sur ce signal dont la sortie est notée Zo, Z1 ... Z4095, 25 2. pour une fréquence centrale Fn,p souhaitée, on calcule le rang n du "bin" Zn de la transformée de Fourier FFT qui est juste inférieur, c'est à dire celui pour lequel Fn, est comprise entre n fois 1.5625 Hz et (n + 1) fois 1.5625 Hz, 3. on calcule le nombre réel x compris entre 0 et 1 qui vaut (Fn,p - 1.5625 30 n) / 1.5625 et l'on sélectionne deux nombres complexes z0(x) et zl (x) fonctions de x (et pré-calculés par commodité), 4. la sortie du filtre est simplement donnée par zO(x) Zn + zl (x) Zn+1. Les coefficients zO et z1 sont optimisés afin que le résultat soit aussi proche que possible de celui que donnerait un filtre réellement centré sur Fn,p. L'un des deux coefficients zO ou z1 peut être choisi REEL (et l'autre corrigé en 5 conséquence) puisque la phase du signal de sortie est sans importance. Chaque filtre est suivi d'un détecteur quadratique d'enveloppe donnant le carré du module de sa sortie (le traitement est généralement fait sur un signal complexe, en bande de base, i.e. centré sur la fréquence 0). La sortie du pème détecteur du nème banc sera notée en,p. 10 La figure 4 donne un exemple de sorties des détecteurs quadratique d'enveloppe, pour P = 100, et N = 7 fréquences, avec un signal d'entrée décalé vers les hautes fréquences. Sommateurs 200 La sortie du pème sommateur est donnée par : N-1 Sp = E en,p n=0 p = 0....P -1 avec en,p.la sortie du pème détecteur du Hème banc de filtre. La figure 5 représente un exemple des sorties des sommateurs pour le cas représenté sur la figure 4, où l'on voit bien l'effet de lissage apporté par la sommation (du fait de la propagation, le gain peut varier beaucoup au voisinage de 20 chacune des N porteuses). Détection des deux plus grands maxima parmi P 300 Le signal issu des sommateurs 200 est transmis à un système de détection des maximas 300 ayant notamment pour fonction de rechercher Deux maxima. 25 En effet, s'il existe deux maxima d'amplitudes voisines, ils sont en général de faibNe amplitude et ne correspondent à rien de précis : il y a probablement uniquement du bruit en entrée du récepteur. 15 Pour la recherche de ces deux maximas, M (le plus grand) et m (le second plus grand), le procédé peut utiliser l'algorithme suivant, écrit en langage C simplifié. /* Initialisation */ M=O position_M = 0 m=0 position_m = 0 /* Boucle pour rechercher un maximum aux positions 1 ... P-2 */ pour( i=1 à P - 2 ) { si(S(i)>=S(i-1)etS(i)>=S(i+1)){ /* S(i) est un maximum local */ si(S(i)>M){ /* S(i) est le nouveau plus grand maximum */ m=M position_m = position_M M = S(i) position_M = i } sinon { si(S(i)>m){ /* S(i) est le nouveau second maximum */ m = S(i) position_m = i } 25 } } } Si position_M est non nul (0 est une valeur illégale) et si m est inférieur à une proportion donnée de M (par exemple 50 %), M et sa position (position_M) 30 peuvent être considérés comme valables. On a trouvé 2 maximas. Sinon ils sont ignorés par la suite (on se contente de mettre M à zéro). La valeur de la position des deux maximas est ensuite affinée par exemple par régression quadratique, c'est-à-dire que les amplitudes des sommes 10 15 20 voisines de la somme maximale M sont approximées par une parabole utilisant 3 ou 4 voisins selon le nombre de voisins disponibles, valeur désignée par la suite valeur corrigée. Pour l'élaboration de l'algorithme de recherche, on pourra utiliser tout 5 autre langage que le langage C. La figure 6 donne un exemple de régression parabolique, toujours sur les mêmes données des figures 4 et 5. La position du maximum de la parabole est la valeur modifiée de position_M qui, d'entier, devient un nombre réel. 10 Pour position_M = 1, la correction est égale à : 1 1 S0 - 3 S 1 - 7 S2 - S3 10 ( S0 - S1 - S2 + S3) Pour position_M = 2....P-3, la correction est égale à : 7(2Sp_2 + Sp_1 + Sp+l û 2Sp+2 ) 10(2Sp 2 ûSp 1 -2S ûSP+1 +2Sp+2 Pour position_M=P-2 la correction est égale à : (Sp_4 +7Sp_3 +3Sp_2 -11Sn_1 ) 10(Sp 4-s 3-Sp 2+ Sp 1) En sortie du détecteur de maximas on dispose de deux valeurs de maximas et de leur position. Prédicteur de position 400 Le système de prédiction / remplacement, 400, utilise la valeur de la 20 position corrigée à l'instant t, et aux instants t - 1, t - 2, ..., t - T (t est pris entier pour la commodité des calculs et au moins égal à 4). La valeur de la position corrigée correspond à la valeur de la position déduite du maxima et affinée par régression quadratique. 15 Sa fonction est notamment d'évaluer, à partir du passé, la valeur probable de la position à l'instant t, en considérant d'abord les maximas de plus grande amplitude (Mt_1, Mt_2, ... , Mt_T), car les plus fiables. II est constitué, notamment d'une ligne à retard contenant les T+1 dernières paires M/position_M. Sachant que, sur une durée limitée, la variation du Doppler peut être assimilée localement à une parabole, on prédit la valeur de la position à l'instant t au moyen d'une régression quadratique pondérée. Pour cela, on calcule les coefficients a, b, c d'une parabole minimisant l'erreur quadratique pondérée suivante (où la variable post désigne la valeur de 10 position_M après correction) : E= EMi((ai2 + bi + c) ûpos,)2 i=t-T La valeur prédite est simplement post=a t2+bt+c Son calcul est simple et donné par la série d'équations suivantes : Uk = M,(t-i)k,k=0 4 1=t-T Wk = EM,,pos,(tûi)k,k=0 2 i=t-T w0(u2u4 ùU3)+W1\U2U3 ùU1U4/+W2(u1u3 ùU2 )) = r l post U0(U2U4 ùU3)+Ul(U2U3 ùU1 U4)+U2(UlU3 ùU2 ) Choix de la position à utiliser, 500 La valeur prédite pour la position est transmise à un dispositif 500 de 20 choix final de la position à utiliser. Ce dernier reçoit aussi l'information issue de l'étape de validation/affinage de la valeur de position. 15 8i = 8max (2 pos P + 1 P Si la position corrigée est non valide, le calcul est fait avec Mt = 0, et l'on remplace alors la valeur prédite post par sa valeur estimée au moyen de la formule donnée ci-dessus et Mt par la moyenne des valeurs précédentes. Si la position corrigée est valide, on calcule sa valeur estimée avec la valeur mesurée de Mt. Si la différence entre la position estimée et la position corrigée est suffisamment faible (en valeur absolue), on garde la valeur de post sinon, on fait comme il a été explicité précédemment. Ce mécanisme permet d'éliminer les valeurs aberrantes. L'initialisation se fait avec T + 1 valeurs M; égales et petites et des 10 valeurs post correspondant à un décalage en fréquence nul, soit (P - 1) / 2. L'approximation de la variation du Doppler par une sinusoïde serait peut-être plus appropriée au canal acoustique sous-marin (à cause du caractère périodique de la houle. Les calculs sont cependant plus complexes. La conversion de position en "Doppler" peut se faire après la correction 15 de la position, simplement en remplaçant post par : Les calculs seront les mêmes, simplement en remplaçant pos; par 8;. Filtrage final, 600 Ce module reçoit la valeur issue de la conversion de position en Doppler. Il comporte par exemple plusieurs filtre passe-bas et un commutateur. 20 La dernière étape du traitement consiste à effectuer un filtrage passe-bas conséquent sur les valeurs des décalages Doppler pour obtenir une valeur lissée finale 8 aussi fidèle que possible tant que la fréquence instantanée de 8 reste comprise entre deux fréquences -Fo et + Fo. On utilisera par exemple un filtre passe-bas à phase linéaire avec un nombre impair 25 de coefficients, à savoir : i=t 8= 1h;_8 i=i-2K avec i=2K (1)Eh, ==1 i=o et (2)h, = h2K_ ,i = 0....K -1 La condition (1) garantit qu'un décalage constant en entrée sera inchangé en sortie. La condition (2) exprime la phase linéaire, qui garantit un retard constant quelle que soit la vitesse de variation instantanée de 8. Ce retard est égal à K fois la période d'échantillonnage. Il faut lui ajouter le retard dû aux traitements précédents, notamment l'étape de filtrage pour laquelle tous les N x P filtres doivent avoir le même retard, ce qui est le cas avec l'exemple de mise en oeuvre proposé. Pour le design exact du filtre, qui doit être aussi "plat" que possible entre - Fo et + Fo, on peut utiliser la méthode suivante: 1 - se donner un filtre idéal dont la réponse en fréquence est donnée par : 2 h(f)=1+ a( Fo -Fo 3 - calculer sa réponse impulsionnelle échantillonnée à des intervalles égaux à T, et la pondérer par une fenêtre de Hamming (par exemple) : hi=h((i-K)T ) (0.54 + 0.46 cos (Tt i-K K + 0.5 i=0...2K Pour T = 80 ms, et K = 16, par exemple, on a une réponse représentée à la figure 9. 4 -ajuster le coefficient a pour avoir une réponse en fréquence aussi plate que possible dans la bande 5 - normaliser de telle sorte que la somme des coefficients vaille 1. On a représenté en figure 10 un exemple de réponse en fréquence 10 obtenue (la valeur de a est volontairement un peu trop élevée). Dans le cas où l'on désirerait avoir des valeurs du Doppler à une cadence plus rapide que toutes les 80 ms, on doit utiliser une technique d'interpolation pour obtenir les valeurs intermédiaires. La technique bien connue de l'homme de l'art et dénommée "filtres 15 polyphases" permet d'obtenir les échantillons intermédiaires. Dans un exemple de mise en oeuvre, les coefficients h'; d'un filtre donnant des échantillons décalés en temps sont calculables simplement par la formule générique 20 hP = (û1)P+' sin(2T6) ( ho + (û1)`h, ( 1 + 1 )) pûB ,_, p+i ûe pûiûB suivante, où 0 dépend du retard souhaité. La valeur 0 = 0 correspond au filtre d'origine. (û1)P+' sin(,rO) ho + (û1) h 1 + 1 )) ( pùB ,( p+iùB pùiùB 25 REVENDICATIONS 1 - Système d'estimation et de suivi du décalage en fréquence pour un système de transmission dont la largeur de bande n'est pas négligeable devant la fréquence centrale, le système émettant une forme d'onde incluant des porteuses fixes de niveau sensiblement supérieur à celui des porteuses véhiculant l'information, les porteuses ayant des fréquences F0 . FN-1 et le signal utile porteur de l'information caractérisé en ce qu'il comporte au moins les éléments suivants : • N bancs (100) de filtres passe-bande, chaque banc comportant P filtres, adaptés à couvrir l'ensemble des valeurs que peuvent prendre les fréquences F0 ... FN-1 après décalage Doppler, • La sortie d'un filtre est reliée à un détecteur quadratique d'enveloppe, • P sommateurs (200) adaptés à additionner les sorties correspondantes des 15 N bancs de filtres, • Un système (300) de détection de la position des deux plus grands maxima parmi les P sommes et de validation de la position du plus grand maxima, • Un système (400) de mémorisation des positions des deux plus grands maximas et de leurs poids et de prédiction de la valeur courante, 20 • Un système (500) de validation/remplacement de la position trouvée, • un système (600) convertissant les positions en Doppler estimé et filtrant celui-ci par un filtre passe-bas...CLMF: 2 Procédé permettant d'estimer et de suivre l'effet Doppler pour un signal composé 25 d'une forme d'onde incluant des porteuses fixes de niveau sensiblement supérieur à celui des porteuses véhiculant l'information, les porteuses ayant des fréquences F0 FN-1 et du signal utile porteur de l'information dans un système de transmission dont la largeur de bande n'est pas négligeable devant la fréquence centrale caractérisé en ce qu'il comporte au moins les étapes suivantes : 30 1. Transmettre le signal composite comportant le signal utile porteur de l'information et le signal de référence à une étape de filtrage utilisant les fréquences des porteuses, (100), 2. Déterminer pour chaque filtre, l'enveloppe du signal, en,p, 3. Sommer les différentes sorties en,p, (200), 4. Déterminer les valeurs des deux plus grands maximas, M et m, pour plusieurs instants, t, (300), 5. Déterminer la valeur probable de la position à l'instant t, en considérant les maximas trouvés qui définissent de façon unique la valeur de décalage Doppler associée, (400, 500), 6. Effectuer une filtrage passe-bas sur les valeurs de décalage Doppler pour obtenir une valeur lissée (600). 3 - Procédé selon la revendication 2 caractérisé en ce qu'il comporte une étape où l'on interpole la valeur lissée à des instants intermédiaires au moyen de filtres passe-bas dérivés du filtre de la revendication précédente. 4 - Procédé selon la revendication 2 caractérisé en ce que l'étape de filtrage 1 15 comporte au moins les étapes suivantes : • On effectue le fenêtrage du signal d'entrée par une fenêtre de Hanning (simple cosinus sur-élevé) sur une durée de 160 ms, avec recouvrement de 80 ms entre deux fenêtres successives, • aux 1.024 points obtenus, on ajoute 3.072 points nuls, puis : 20 • on effectue la transformée de Fourier rapide FFT 4096 points sur ce signal dont la sortie est notée Zo, Z1 ... Z4095, • pour une fréquence centrale Fn,p souhaitée, on calcule le rang n du "bin" Zn de la transformée de Fourier FFT qui est juste inférieur, • on calcule le nombre réel x compris entre 0 et 1 qui vaut (Fn,p -1.5625 n) / 25 1.5625 et II 'on sélectionne deux nombres complexes z0(x) et z1 (x) fonctions de x, • on détermine la sortie du filtre par z0(x) Zn + zl (x) Zn+1. 5 - Procédé selon la revendication 2 caractérisé en ce que l'on affine par régression 30 quadratique la valeur des deux maximas trouvés.6 û Procédé selon la revendication 2 caractérisé en ce que l'étape 6 de filtre passe bas est effectuée avec un filtre passe-bas à phase linéaire avec un nombre impair de coefficients.