Brevet FR2754907A1 — publié le 19980424

Procédé de mesure du décalage Doppler dans un système de détection utilisant des codes ambigus. La présente invention se rapporte aux procédés de mesure du décalage Doppler dans un système de détection, par exemple un radar ou un sonar, utilisant des codes ambigus (Doppler x temps). Plus précisément, elle concerne les codes d'émissions pour ce type de système de détection qui permettent la mesure fine du décalage Doppler des échos reçus d'une cible en mouvement, avec une précision ajustable indépendamment des paramètres concernant la durée T, la fréquence initiale f(O) et la fréquence finale f(T) de l'impulsion codée selon ce système. La modulation des impulsions d'un sonar ou d'un radar avec des codes particuliers permettant d'améliorer la précision de la mesure est bien connue et elle est notamment décrite très largement dans l'ouvrage de A.W. RIHACZEK, intitulé "Principles of high resolution radar", publié chez Mc Graw-Hill. On y expose en particulier que les codes d'émission peuvent être classés en deux catégories: Les codes non ambigus, pour lesquels la fonction d'ambiguïté temps x Doppler présente, autour du maximum, deux directions principales parallèles aux axes temps et Doppler. II est alors connu que pour de tels codes, avec une amplitude constante (ou peu modulée) pendant la durée du code, la tolérance Doppler est inversement proportionnelle à la durée du code, ce qui ne permet pas de fixer indépendamment la durée d'émission et la tolérance Doppler. Avec les durées d'émission utilisées en technique sonar, de tels codes conduisent à un nombre de copies très importants, de l'ordre de 1000. Les codes ambigus, pour lesquels la fonction d'ambiguïté temps x Doppler présente une crête sécante-en zéro avec les axes temps et Doppler. Cette crête est rigoureusement linéaire dans le cas du code à Modulation de Fréquence Hyperbolique (MFH), qui a la particularité d'être parfaitement tolérant au Doppler. Dans ce cas, la tolérance Doppler est entièrement déterminée par la largeur de bande B du code émis et est en général insuffisante pour la mesure précise du Doppler. Le code à Modulation de Fréquence Linéaire (MFL) appartient également à la classe des codes ambigus. Cependant, une fois la durée T et la bande B du code fixées, la tolérance Doppler est donnée par, c étant la vitesse du son: 2V3 ry 1,5c (1) BT Cette tolérance est alors en général trop petite par rapport au besoins et conduit à un surdimensionnement du nombre de copies de réception. Un autre exemple de code ambigu est décrit dans le brevet US 5,480,504, délivré le 2 janvier 1996 au nom d'ATLAS Electronics Gmbh. Dans ce code, la fréquence instantanée est une fonction d'une puissance du temps. Cette puissance est donnée par un paramètre qui varie entre 0 et 1 et lorsque ce paramètre est égal à 1, on retombe sur le code à Modulation de Fréquence Linéaire. Dans ce cas-ci, on a besoin d'un très grand nombre de copies de réception et le nombre minimal de copies est sensiblement égal à celui du code MFL. Pour pallier ces inconvénients, I'invention propose un procédé selon la revendication 1. D'autres particularités et avantages de l'invention apparaîtront clairement dans la description suivante, présentée à titre d'exemple non limitatif en regard de la figure annexée qui représente un diagramme permettant d'expliciter les paramètres utilisés dans le texte. On a représenté sur la figure 1, un paramètre ss représentant l'écart Doppler en fonction des retards des échos affectés de cet écart Doppler. Les points représentatifs, choisis en coupant à -3 dB l'ensemble des signaux par rapport à l'intensité maximum de réception, sont répartis dans une ellipse allongée, centrée sur l'origine des axes et inclinée par rapport à ceux-ci. Le paramètre ss est donné par la formule suivante, dans laquelle V est la vitesse de la cible à l'origine des échos: 2V reçue -1 (2) '3=-- C f émise Dans ces conditions la tolérance Doppler, correspondant à l'étendue en abscisse de cette ellipse, est égale à 2ss3 = 2 2V3 . c Cette tolérance Doppler correspond à un écart sur les retards égal à 2t3. Pour permettre de diminuer le nombre de copies nécessaires pour traiter correctement les signaux de réception, I'invention propose d'introduire un paramètre supplémentaire dans le code utilisé pour moduler la fréquence des impulsions d'émission. Ce paramètre supplémentaire est choisi pour permettre d'ajuster la tolérance Doppler indépendamment de B et T, et pour obtenir une loi de modulation telle que: - pour une valeur particulière du paramètre, le code se réduit au code FMH, - le code présente une crête d'ambiguïté Doppler-distance aussi linéaire que possible, de façon à posséder de bonnes propriétés d'interpolation pour la mesure fine du Doppler et de l'instant d'arrivée de l'écho. Cette interpolation s'effectuera à la réception des échos après l'opération de corrélation. On dispose ainsi, en fonction du rapport signal/bruit des échos, d'une précision de mesure bien meilleure que l'écartement des copies de réception, en Doppler. Plus particulièrement, I'invention propose d'utiliser la loi de fréquence instantanée suivante, où t est pris sur la durée de l'impulsion T, et où x, qui est positif ou négatif, fixe la tolérance Doppler: Le calcul montre que l'on a alors sensiblement: 1,74c (4) IX BT et Dans ces conditions, pour déterminer les paramètres x et A qui définissent de façon complète cette loi de fréquence, on se fixe les fréquence f(0), f(T), la tolérance Doppler V3 et la durée de l'impulsion T qui sont déterminées pour des raisons opérationnelles, et l'on obtient alors ces paramètres x et A à l'aide des formules suivantes: Dans ces conditions, les précisions sur les mesures de la vitesse radiale et du retard sont données par les formules suivantes: Dans ces formules, % est le rapport signallbruit des échos et le calcul est effectué en considérant que le porteur est immobile. Toutefois, les résultats numériques sont encore valables-pour les vitesses habituelles de ce porteur. Dans une application numérique correspondant aux valeurs suivantes des paramètres: f(o) = 6000 Hz f(T) = 7000 tiz T=4s V3=#Cx 1,74 =2,5 m/s |x|BT la valeur du coefficient x obtenue à partir de ces paramètres est égale à: lx| x | = 0,261. On peut alors effectuer l'interpolation entre les copies en utilisant des copies espacées de #V= =2V3 = 1,67 m is pour la vitesse et 3 #T = 0,33 ms pour le retard. Dans ces conditions, en considérant une 3B cible pour laquelle le rapport signal/bruit des échos de réception est égal à 14 dB, on obtient pour les précisions sur les mesures de la distance radiale et du Doppler les valeurs suivantes: # =93=18,7ms correspondant à 14 m 5 2,5 #V = 0,5 m/s 5 Le nombre nj de copies est donné par la formule: n i = 2VMAX (10) Dans cette formule VMAX est la vitesse maximum que peut atteindre la cible. Le nombre de copies nécessaires dans l'art antérieur est sensiblement égal à celui nécessaire en utilisant un code MFL. II est inversement proportionnel à V3 qui est obtenu dans ce cas par la formule: V3 = 1,5c (11) si Dans le cas de l'invention, utilisation d'un code pseudo-MFH, ce nombre de copies est proportionnel à V3 qui est lui-même donné par la formule: 1,74 xC (12) 0,261 BT Le rapport de ces deux valeurs de V3 donne donc le rapport du nombre de copies nécessaires entre l'art antérieur et l'invention. On constate que dans le cas de l'exemple numérique décrit plus haut l'invention permet d'utiliser dix fois moins de copies. REVENDICATIONS 1. Procédé de mesure du décalage Doppler dans un système de détection utilisant des codes ambigus, dans lesquels le signal d'émission est formé d'impulsions de longueur T présentant une fréquence initiale f(0) et une fréquence finale f(T) = f(0)+B, caractérisé en ce que le long de chacune de ces impulsions la fréquence instantanée f(t) varie selon une fonction pseudo-hyperbolique déterminée par un paramètre x tel que, lorsque ce paramètre est nul la fonction se réduit à une fonction hyperbolique, et que lorsque ce paramètre n'est pas nul on puisse à la réception de ces signaux faire suivre les opérations de corrélation d'une opération d'interpolation. 2. Procédé selon la revendication 1, caractérisé en ce que la fonction est donnée par la formule f(t) = f(o) dans laquelle le paramètre x est déterminé en fonction d'une valeur choisie de la tolérance Doppler V3 par la formule x= et le paramètre A par la formule 3. Procédé selon la revendication 2, caractérisé en ce que f(0) =6000 Hz, f(T) =7000 Hz, T =4s et V3 2,5 m/s correspondant à x = 0,261.