Brevet FR2736728A1 — publié le 19970117

PROCEDE ET DISPOSITIF D'OPTIMISATION DE LA COMBINAISON DE FREQUENCES D'UN RADAR AEROPORTE L'invention se situe dans le domaine des procédés et dispositifs destinés à permettre une meilleure détection par un radar aéroporté des cibles se trouvant dans une zone localisée en distance et génératrice pour ledit radar d'un bruit de fouillis important (clutter). Le bruit de fouillis (clutter) peut avoir pour origine le relief, la pluie et d'une façon générale un parasitage ou un brouillage quelconque, la seule condition étant que la zone soit localisée en distance. L'invention est particulièrement adaptée aux radars aéroportés du type Doppler à impulsion ayant des capacités dites "look down" 'est à dire aptes à détecter des cibles se trouvant a une altitude moins élevée que le radar lui-même. On sait qu'un radar Doppler à impulsion à moyenne fréquence de récurrence (MFR) est ambigu à la fois en distance et en vitesse. L'ambiguïté en distance provient du fait que lorsqu'un écho est reçu d'une cible on ne sait pas si cet écho provient de la dernière impulsion envoyée par le radar, ou de l'une des impulsions précédentes ayant parcouru une distance totale égale à nC Tr et donc en provenance d'une cible se situant à une distance égale à nCtr . (C désigne la vitesse de la lumière, 2 Tr le temps séparant deux impulsions successives et n un nombre quelconque). La distance D = C Tr est appelée distance 2 ambiguë du radar fonctionnant à la fréquence de récurrence F r T r. L'ambiguïté en fréquence provient du fait que le signal reçu est composé d'un certain nombre d'impulsions séparées entre elles par une durée égale à Tr. I1 s'ensuit que la transformée de Fourier du signal reçu donne un spectre entièrement contenu entre 0 et Fr = 1 Tr La fréquence Doppler d'une cible est égale à 2 Vr Vr désignant la vitesse relative algébrique radiale de la cible par rapport au porteur du radar, et # la longueur de l'onde émise. I1 s'ensuit que l'examen du spectre du signal reçu ne permet de connaître que le reste modulo F r de la fréquence 2 Vr Doppler # , sans pouvoir déterminer si la fréquence observée est la fréquence Doppler de la cible, ou si elle doit être augmentée d'une ou plusieurs fois la valeur de la fréquence de récurrence Fr La distance et la vitesse minimum que l'on peut affecter à la cible détectée sont appelées respectivement distance repliée et fréquence repliée. Ce phénomène est bien connu en radar, de même qu'un moyen pour déterminer distance et fréquence Doppler vraies qui consiste à faire varier périodiquement la fréquence de récurrence. En effet si une cible se trouve à une distance réelle D du radar, pour une première fréquence de récurrence Fr1 on aura : D=n1 D1 + p1 d Dans cette formule n1 et p1 sont des entiers, D1 est la distance d'ambiguïté correspondant à la fréquence F d est la valeur de la cellule résolution de distance du radar. Dans cette formule n1 est inconnu et D peut donc prendre les différentes valeurs discrètes de la progression arithmétique obtenue lorsque n1 prend les valeurs successives 0, 1, 2, 3... Si l'on change Frî pour une fréquence Fr2 on obtient D=n2 D2 + p2 d En calculant les résultats des valeurs discrètes que peut prendre D lorsque n2 prend les valeurs successives 0,1, 2.. et en les comparant aux valeurs possibles de la série précédentes, on trouve des valeurs de n1 et n2 qui donnent des valeurs égales pour D. On confirme en général la valeur de D ainsi trouvée en émettant une ou plusieurs autres fréquences de façon à diminuer la probabilité de fausse alarme. L'ensemble des fréquences que l'on applique successivement constitue une combinaison de fréquence. Si les fréquences successives Frî > Fr2...Fri...FrM constituant une combinaison sont bien choisies on parvient par ce procédé à lever toute ambiguïté sur la distance de la cible et par un procédé analogue sur sa vitesse radiale. Les critères de rapidité de convergence ne sont pas les seuls critères à prendre en compte pour établir une combinaison de fréquences de récurrence. Des impératifs technologiques tels que puissance de l'émetteur, facteur de forme, imposent de ne pas trop s'écarter d'une valeur de fréquence de récurrence moyenne. Des impératifs de mode imposent une fréquence de récurrence maximum, une durée d'impulsion fonction du nombre de quanta-distance pi que l'on souhaite avoir entre deux impulsions successives, un nombre d'échantillons minimum pour l'analyse spectrale, la portée. Enfin les capacités des circuits de traitement du signal reçu et de calcul viennent encore limiter la marge de choix des combinaisons possibles. De tels impératifs concernent tout aussi bien les radars fixes au sol que les radars aéroportés. Dans les radars fixes les ambiguïtés peuvent être plus facilement levées en considérant que les fréquences repliées, situées au dessous d'un seuil fixé, correspondent à des échos fixes tels que des reliefs et que, au-dessus de certaines valeurs, la dérive de fréquence observée ne peut correspondre à des échos de cibles réelles, étant donnée leur trop grande vitesse. Pour les radars aéroportés, la gamme de vitesse à traiter est plus importante, allant de vitesses qui peuvent être proches de 0, lorsque l'avion porteur du radar se rapproche d'une cible par l'arrière ; ou bien à des vitesses très importantes, lorsque l'avion porteur et la cible se rapprochent en face à face. De plus les échos en provenance du sol ont des vitesses relatives plus ou moins grandes selon la direction dont ils proviennent. La figure 1 représente dans le domaine distance-fréquence, la zone dans laquelle se situent les échos de sol. La représentation de la figure 1 est l'allure générale de la zone dont certains radars pourront s'écarter. La zone représentée figure 1 comprend de façon typique : une première partie de zone 1 se situant dans les vitesses négatives et correspondant aux échos de sol entrant par des lobes arrière de l'antenne du radar ; une zone 2 au niveau de la vitesse 2Vcos cl a, étroite en vitesse, mais s'étendant sur une grande distance qui correspond aux échos de sol entrant par le lobe principal de l'antenne lorsque ce dernier est dépointé d'un angle a par rapport à la vitesse horizontale de l'avion porteur, une zone 3 correspondant à des réflexions sur le radome ; enfin une zone 4 correspondant à des lobes diffus de l'antenne. En parallèle à l'axe des fréquences-vitesse, on a matérialisé également des segments 5, 6, 7, 8, correspondant respectivement à des plages de vitesses pour lesquelles l'avion porteur et la cible vont en s'éloignant (5), en se rapprochant par l'arrière (6), en se rapprochant par l'avant à des vitesses subsoniques (7), puis supersoniques (8). Pour un même radar aéroporté les limites de la zone représentée figure 1 varient en fonction du relief de la zone survolée, de la vitesse et de l'altitude de l'avion, et de l'angle de dépointage de l'antenne. A chacun des points de la zone représentée figure 1, correspond pour chaque fréquence de récurrence, un point dans le domaine distance fréquence repliée. Ainsi lorsqu'une cible a une distance et une vitesse repliée se situant dans la zone repliée des échos de sol, il arrive qu'elle "disparaisse" pour réapparaître quelques instants plus tard. Ce phénomène a été illustré de façon fictive sur les figures 2-1 à 2-4 pour le domaine distance. Sur la figure 2-1 est représentée une cible se trouvant à une distance réelle (on dit aussi vraie ou non ambiguë) comprise entre 3 et 4 quanta-distance. La longueur du T quanta-distance est égale à QT= C 2T, C représentant la vitesse de la lumière et T la durée d'une impulsion. Cette cible apparaît comme représentée figure 2-2 dans le domaine repliée. A partir de là, deux cas sont possibles - cas a) représenté figure 2-3 : La zone repliée des échos de sol englobe la distance repliée où se trouve la cible - cas b) représenté figure 2-4 La zone repliée des échos de sol n'englobe pas la distance repllée où se trouve la cible. On peut, pour la même fréquence de récurrence, représenter les mêmes éléments dans le domaine fréquence. Si l'on se trouve dans le cas (a) à la fois dans le domaine fréquence et dans le domaine distance, la cible est noyée dans les échos de sol ; dans le cas contraire la cible est susceptible d'être vue. Les "trous" de détection apparaissent davantage à basse ou très basse altitude en particulier lorsque l'avion porteur se rapproche de la cible par l'arrière ou lorsqu'il la croise à grande vitesse en présentation face à face. Le problème technique que l'on cherche à résoudre par la présente invention est de choisir pour chaque configuration particulière de vol une combinaison de fréquences de récurrence pour laquelle la cible ne disparaît pas ou disparaît aussi peu que possible. Cela consiste à choisir la combinaison de fréquences de récurrence pour laquelle le point représentant le repliement de la cible dans le domaine distance- vitesse ou ce qui revient au même, temps fréquence Doppler, est toujours en dehors de la zone repliée correspondant aux échos de sol ou à défaut y est le moins souvent possible. Cette zone sera appelée par la suite zone polluée repliée. Dans une étape préalable à la mise en oeuvre du procédé il conviendra d'avoir suffisamment caractérisé le comportement du radar pour connaître - la zone polluée en fonction des paramètres pertinents tels que altitude du porteur, site de l'antenne, niveaux des lobes diffus dus à l'antenne et au radome - l'étendue en distance d'application du procédé selon l'invention en fonction de la surface équivalente radar des cibles que l'on cherche à trouver. Cette distance est mesurée en nombre de quanta-distance (NCD).Cette distance est intermédiaire entre les courtes portées et les grandes portées car pour les courtes portées la cible est détectée malgré le fouillis car le rapport signal/bruit de fouillis est important et pour les grandes portées la détection est perturbée, mais les fluctuations de la cible ne permettent pas d'espérer une amélioration. I1 en résulte que la zone distance d'application du procédé correspond à peu prés à la zone où les niveaux de bruits en provenance du sol et le niveau de signal de la cible sont équivalents. La détermination de la zone d'application du procédé ne fait pas partie du procédé selon l'invention. - L'ensemble des combinaisons de fréquences susceptible d'être utilisées, on les notera par la suite K1, K2,. . KNCB, NCB désignant le nombre de combinaisons possibles. Toutes les combinaisons comportent un nombre identique de fréquences de récurrence. Ce nombre sera désigné par NFR. Ces éléments étant fixés, on peut dire que l'invention est relative à un procédé applicable à un radar porté par un mobile, et destiné à choisir parmi NCB combinaisons possibles de NFR fréquences de récurrence chacune, la combinaison la meilleure pour détecter, à l'intérieur d'une zone distance prédéterminée dite zone testée composée de NCD quanta-distance, la présence d'une cible malgré la présence d'échos de sol se trouvant dans une zone dite polluée supposée connue et qui, pour chaque fréquence de récurrence, se replient à l'intérieur d'une zone dite zone polluée repliée, procédé caractérisé en ce que - on calcule pour chaque combinaison de fréquence de récurrence une note qui est une fonction croissante du nombre de fois où, pour les fréquences de récurrence entrant dans la combinaison les différents quanta-distance de la zone testée se replient hors de la zone polluée repliée et on choisit la combinaison ayant obtenu la note maximum. De préférence l'inventions sera réalisé par un procédé applicable à un radar porté par un mobile, et destiné à choisir parmi NCB combinaisons possibles de NFR fréquences de récurrence chacune, la combinaison la meilleure pour détecter, à l'intérieur d'une zone distance prédéterminée dite zone testée composée de NCD quanta-distance, la présence d'une cible malgré la présence d'échos de sol se trouvant dans une zone dite polluée supposée connue et qui, pour chaque fréquence de récurrence, se replient à l'intérieur d'une zone dite zone polluée repliée, procédé caractérisé en ce que 1) pour la première fréquence de récurrence entrant dans la première combinaison, on calcule, pour chaque quanta-distance de la zone testée, la distance repliée de ce quanta-distance 2) on compare la distance repliée de chaque quantum de la zone testée à la distance repliée de chacun des quanta-distance de la zone polluée repliée, et l'on attribue au quantum distance de la zone testée un coefficient E qui est q égal à 1 si le quantum testé se situe hors de la zone polluée 3) on recommence les opérations 1 et 2 pour chacune des NFR fréquences de récurrence entrant dans la première combinaison 4) on additionne les coefficients E obtenus pour q chacun des NCD quanta-distance testé, constituant ainsi un vecteur colonne "somme eq" dont les NCD composantes sont désignées par un nombre L, L pouvant prendre les valeurs entières de 0 à NFR 5) on constitue un vecteur colonne répartition désigné RP et comprenant NFR+1 composantes désignées RPO, RP1, RP1... RPNFR > ces composantes ayant initialement la valeur 0;; 6) on incrémente la composante RPo du nombre de fois où apparaît dans le vecteur colonne constitué à l'étape 4, puis la composante RP1 du nombre de fois ou apparaît 1 dans le même vecteur colonne, et ainsi de suite jusqu'à incrémenter RPNFR du nombre de fois où apparaît NFR dans le vecteur colonne constitué à l'étape 4 7) on calcule une première note NT selon la formule dans laquelle CF est un coefficient qui est une fonction de L et ou RP(L) est la composante de rang L du vecteur RP 8) on recommence les étapes 1 à 7 pour la seconde combinaison calculant une seconde note NT2;; Si NT2 est supérieure à NT1 alors NT devient NT2, NT reste NT1 dans le cas contraire 9) on recommence les étapes 1 à 7 pour chacune des NCB-2 combinaisons restantes, calculant à chaque fois une note NT qui vient se substituer à la note maximum NTM calculée selon le processus décrit pour les étapes l à 8 10) on sélectionne la M ème combinaison, combinaison que l'on garde jusqu'à ce qu'intervienne un changement prédéterminé affectant le mobile. La note NT que l'on calcule est un nombre représentatif de l'espérance de détection d'une cible pour la combinaison testée. Dans les cas où l'on souhaite raccourcir au maximum la durée du calcul on pourra prendre CF = 1 et, dans ce cas, (L) la note NT est égale à un nombre qui est la somme des nombres 1 apparus au cours de l'étape 2 lors de la réalisation des étapes 1, 2, 3, correspondant à la combinaison notée. On remarquera cependant que la note ainsi attribuée n'est pas la plus représentative que l'on puisse trouver pour qualifier l'espérance de la détection de la combinaison de fréquences ainsi testée. I1 a été vu précédemment que pour lever l'ambiguïté il faut pouvoir recevoir un écho de la cible se repliant en dehors de la zone polluée repliée sur au moins deux fréquences de récurrence. I1 en résulte que si, pour un quantum distance, une seule fréquence de récurrence permet un repliement hors de la zone polluée, cette unique réception d'écho ne permettra pas pour le quantum considéré la levée d'ambiguïté. On remarquera également que si pour un quantum distance > deux fréquences de récurrence permettent un repliement de l'écho reçu hors de la zone polluée repliée, encore faut-il que l'écho reçu de la cible soit effectivement détecté pour chacune des fréquences. Lorsqu'un écho est détecté, la détection obtenue peut correspondre à une cible réelle, ceci avec une probabilité P appelée probabilité de détection ; ou bien à ce que l'on appelle une fausse alarme, avec une probabilité Pfa appelée probabilité de fausse alarme. Dans ces conditions, on comprend que la probabilité de détecter une cible réelle augmente avec le nombre de fois où il devient possible pour un quantum distance de replier la cible hors de la zone polluée repliée. Pour un quantum donné des repliements hors de la zone polluée repliée sont nécessaires pour deux fréquences de récurrence au moins mais, compte tenu de la probabilité de détection, on pourra exiger, pour fixer la note, de ne prendre en compte pour un quantum distance donné que les combinaisons de fréquence donnant lieu à un nombre de repliements hors de la zone polluée repliée supérieur à un certain seuil k, (seuil au minimum égal à deux, mais de préférence supérieur). La note doit tenir compte également de la probabilité de détection. C'est pourquoi de préférence le coefficient CF sers choisi égal à cF(L=L CL1 PI (l-P)(L-j) pour L)k j=k et CFL: O pour Lsk Dans la formule ci-dessus, j est un nombre pouvant prendre toutes les valeurs entières de j=k à j=L ; Ci est le nombre de combinaisons que l'on peut former en choisissant; objets parmi L ; et p est la probabilité de détection élémentaire d'une cible de surface équivalente radar de référence, pour un bruit thermique à température donnée et à probabilité de fausse alarme PFA donnée. Le procédé de réalisation de l'invention ainsi qu'un mode de réalisation seront maintenant décrits en référence aux dessins annexés dans lesquels - la figure 1, représente l'allure générale dans le domaine espace vitesse ou temps fréquences d'une zone pour laquelle sont reçus des échos de sol pouvant donner lieu à détection par le radar - la figure 2, illustre le phénomène de repliement d'une cible dans ou hors la zone polluée - la figure 3, représente le synoptique général du principe d'optimisation du choix de la combinaison de fréquence de récurrence - la figure 4, représente un exemple de tableau réalisé conformément aux étapes 1, 2, 3 du procédé selon l'invention;; - la figure 5, représente le synoptique de calcul des coefficients E égaux à O ou 1 constituant le tableau de la figure 4 - la figure 6, représente le synoptique de calcul du nombre de détections possibles ND c'est à dire du nombre de repliements hors de la zone polluée repliée pour un quantum distance DTN; - la figure 7, représente le synoptique de calcul du vecteur colonne répartition, pour la zone distance testée et pour une combinaison de fréquences de récurrence - les figures 8 et 9, représentent des schémas de montage du dispositif dans un mode de réalisation; - les figures 10 et 11, sont des courbes illustrant l'intérêt du procédé. Les figures 1 et 2, ont été décrites ci-dessus. La figure 3, représente le synoptique général du principe d'optimisation du choix de la combinaison de fréquences de récurrence selon l'invention. Chacune des combinaisons étant indexée au moyen d'un indice K allant de 1 à NCB, NCB représentant le nombre de combinaisons possibles, on effectue de façon itérative le calcul sur chacune des combinaisons, en commençant par la première et en suivant jusqu'ver une combinaison K qui est celle représentée figure 3 (K On calcule tout d'abord pour le premier quantum de la distance testée, et pour chacune des NFR fréquences entrant dans la combinaison, un coefficient E qui est égal à 1 si le q quantum testé se replie hors de la zone polluée repliée et à 0 si le quantum testé se situe dans la zone polluée repliée. Cette opération est effectuée dans un module de calcul R1 représenté en détail figure 5. On calcule ensuite pour le quantum distance testé le nombre ND de détection possible pour la combinaison de fréquence K testée. Ce nombre ND est égal à la somme des coefficients E obtenus par le quantum distance testé avec q chacune des fréquences de récurrence entrant dans la combinaison. I1 représente le nombre de fois où le quantum distance testé se replie hors de la zone polluée repliée pour la combinaison testée. Cette opération est effectuée dans un module de calcul R2 représenté en détail figure 6. On constitue un vecteur colonne que l'on appelle répartition RP, comportant NFR+1 composantes repérées R?0... ..... . RPNFR. Les composantes de ce vecteur sont initialisées à 0. Chaque fois que l'on a terminé un calcul de ND, c'est à dire chaque fois que pour un quantum distance testé on a calculé les coefficients E correspondant à chacune des NFR q fréquences de récurrence de la combinaison testée et qu'on en a fait la somme, on augmente de 1 la composante de RP qui porte l'indice ND. Ainsi, lorsque le calcul est terminé, les composantes du vecteur colonne RP représentent, pour la composante RPO, le nombre de quanta-distance où la combinaison testée ne donne lieu à aucun repliement hors de la zone polluée repliée, pour la composante RP1 le nombre de quanta-distance où la combinaison ne donne lieu qu a un seul repliement hors de la zone polluée repliée et ainsi de suite jusqu'à la composante RPNFR qui représente le nombre de quanta-distance qui se replient hors de la zone polluée repliée quelle que soit la fréquence de la combinaison employée. Pour chaque combinaison testée les quanta-distance sont ainsi répartis en quanta se repliant 0, 1, 2 . . NFR fois hors de la zone polluée repliée. Le mode de calcul de ce vecteur RP fait l'objet d'un module R3 expliqué en détail figure 7. On calcule ensuite en fonction des composantes de ce vecteur RP la note que l'on attribue à la combinaison comme indiqué plus haut. C'est l'objet du module R4. Cette note est ensuite comparée à la meilleure note obtenue précédemment. Si la note qui vient d'être calculée est meilleure c'est la dernière note qui est retenue et l'indice de la combinaison est mémorisée (R5). On recommence ensuite jusqu'à ce que les NCB combinaisons aient été testées. La figure 4a, représente un exemple fictif simplifié en ce sens que le nombre de lignes est petit, du tableau que l'on crée par réalisation des étapes 1 à 4 du procédé selon l'invention. Ce tableau comporte une ligne par quantum distance testé. Les quanta sont référencées DT1 à DT10 sur la figure 4a. Le tableau comporte une colonne par fréquence appartenant à chaque combinaison testée. Les fréquences sont repérées par des indicatifs QT1 à QT8 sur la figure 4a. Le tableau comporte donc, dans ce cas, 10 lignes parce que le test est effectué sur 10 quanta distances, et 8 colonnes car chaque combinaison de fréquence comporte 8 fréquences. Chacune des cases du tableau ainsi constitué comporte en fonction des résultats du calcul qui sera explicité plus loin un 1 ou un O.La présence d'un 1 signifie que le quantum distance se replie pour la fréquence de récurrence considérée en dehors de la zone polluée repliée. I1 y a un 0 dans le cas contraire. Chaque ligne du tableau est complétée par une case supplémentaire dite case somme ou est inscrite la somme des 1 figurant sur la ligne. Ainsi la ligne DT10 qui ne comporte qu'un seul 1, comporte un 1 dans la case somme, la ligne DT8 qui comporte cinq 1 comporte un 5 dans la case somme. L'ensemble des cases somme constitue une colonne "somme eq11. La somme des chiffres figurant dans cette colonne constitue le nombre total de fois où, pour la combinaison de fréquences, un quantum distance se replie hors de la zone polluée repliée. Dans l'exemple représenté figure 4a, à 10 lignes et 8 colonnes, les chiffres figurant dans chacune des cases sommes peuvent être compris entre 0 et 8. A l'aide de ces chiffres on constitue un vecteur colonne répartition, désigné RP et représenté figure 4b. Ce vecteur comporte une composante de plus que le nombre de fréquences entrant dans chaque combinaison. Dans le cas de la figure 4b, il s'agit donc de 9 composantes répertoriées de RPo à RP8. La valeurs de la composante RPL est égale au nombre de fois où la case somme d'une ligne du tableau de la figure 4a a la valeur L. Ainsi, dans l'exemple de la figure 4, seul le quantum DT2 ne donne lieu à aucun repliement hors de la zone polluée. La valeur figurant dans la case somme est 0 et ce O n' apparaît qu'une fois dans le vecteur colonne somme. On en tire RP,=1. Les valeurs 2 et 4 n'apparaissent pas dans le vecteur colonne somme, les valeurs RP2 et RP4 sont donc égales à 0. Dans une variante de réalisation, on peut vérifier que la somme des valeurs des produits des valeurs de RPL, dans le tableau RP, par un nombre représentant le rang de la valeur de RPL (O, 1, 2, 3... 8) est égal au nombre de cases du tableau de la figure 4a où figure un 1. La figure 5, représente le calcul de la valeur 0 ou 1 figurant dans l'une des cases du tableau de la figure 4a. Cette case sera prise comme étant dans la ième colonne et la combinaison testée sera supposée être la Kème. Cette fréquence correspond à une distance ambiguë qui exprimée en nombre de quanta distance est égale à un nombre que, par convention, on appelle QTk et qui est caractéristique de la fréquence de récurrence choisie. C'est pourquoi on peut aussi l'employer comme indicatif de cette fréquence de récurrence. On sait que, compte tenu de variables affectant le mobile porteur du radar, la zone définie figure 1 pour les entrées de lobes secondaires, se projette sur l'axe des distances pour définir sur cet axe une zone dont les distances à l'origine, exprimées en quanta distance, s'échelonnent de Z?1 à ZPNCP. Cette zone est donc longue de NCP quanta. Les opérations qui sont décrites ci-après sont répétées pour chaque quantum DTM appartenant à l'étendue distance testée. 1) On effectue dans un module 9 la division de la distance DTM par la distance ambiguë QTik. Cette division i. donne un nombre entier et un reste, soit respectivement aM et bM. bM est la distance repliée du quantum testé DTM. 2) On effectue la division de la distance du quantum distance ZPj appartenant à la zone polluée non ambiguë par dans dan un module 10. Cette division donne un nombre entier aZPj et un reste b ZPJ b ZPj est la distance repliée du quantum polluée ZPj. On compare les nombres bZPj et bM dans un module 11. Si bZPj et bM sont différents, cela signifie que les distances repliées des quanta DTM et ZPJ sont différentes, la valeur d'un coefficient E initialisée à 1 est maintenue à 1. q Si par contre bM et bZPj sont égaux, cela signifie que le quantum distance DTM se replie à la même distance que le quantum pollué ZPj ; dans ce cas E est ramené à O. Le calcul q peut cesser pour le quantum DTM. Si E est toujours égal à 1, on incrémente 3 de 1 et q on reprend les opérations figurant à l'intérieur du module R encadré en traits pointillés sur la figure 5. Lorsque ce module R1 est terminé, E a la valeur O ou 1 selon que le quantum testé DTM a une distance repliée b M différente des distances repliées bZPj de chacun des NCP quanta de la zone polluée. Pour compléter la ligne correspondant au quanta distance testée DTM, il convient de recommencer l'opération, décrite figure 5, pour chacune des NFR fréquences de récurrence entrant dans la combinaison K. I1 convient également de constituer la valeur ND représentant la somme des coefficients E pour le quantum DTM. Ces opérations sont représentées figure 6. A l'issue de l'exécution du module R1 le nombre E calculé va incrémenter un nombre ND, (nombre de détections possibles). L'indice i de la fréquence repérée QTik est incrémenté de 1 et la nouvelle valeur de Eq est calculée. L'opération est recommencée Jusqu a ce que l'indice i ait pris les NFR valeurs allant de 1 à NFR. A l'issue de ce calcul, le nombre ND représente le nombre de détections possibles des quantum DTM pour la combinaison K testée. Les opérations aboutissant à la formation de ND forment un module R2 entouré en pointillé sur la figure 6. Il convient ensuite de constituer les vecteurs colonnes "somme eq" et répartition RP. Cette constitution est décrite figure 7. A l'issue du module R2 la valeur ND incrémente de 1 la composante RPND du vecteur colonne RP. Le quantum distance DTM est remplacé par le quantum repéré DTM 1 et le calcul de ND et de RPND est recommencé Jusqu'à ce qu'il ait été effectué pour chacun des NCP quanta distance de la zone testée. A l'issue de ce calcul, le vecteur RP est constitué. Les opérations constituant le calcul des composantes du vecteur RP sont représentées par le module R3, entouré en pointillé sur la figue 7. La description du calcul de la note attribuée à la combinaison testée sera effectué plus tard, à l'occasion de la description des figures 8 et 9 représentant un exemple particulier de réalisation. La figure 8, représente trois compteurs J,I,N, 12, 13, 14 dont les sorties sélectionnent les adresses de mémoires DT QT ZP 15, 16, 17, dont les valeurs alimentant des modules de calcul 18 et 19. Les valeurs calculées dans les modules 18 et 19 alimentent un comparateur 20 qui donne une valeur 1 si les valeurs d'entrée sont égales et une valeur O dans le cas contraire, les coefficients E successifs sont additionnées q dans un sommateur 21 dont la sortie alimente une mémoire ND 22 qui est remise à zéro chaque fois qu'un signal est reçu du compteur 12, N.La valeur contenue dans la mémoire ND est soustraite d'une valeur NFR dans un sommateur 23 z 2 et le résultat de ce sommateur sélectionne l'adresse d'une mémoire 24 RP qui est incrémentée de 1 chaque fois qu'elle est sélectionnée au moyen d'un sommateur Z3 25. Le fonctionnement est le suivant Les compteurs J, I, N, 14, 13, 12 en cascade permettent de réaliser les boucles R1, R2 et R3 respectivement. Les valeurs de DT, QT et ZP sont mémorisées dans les mémoires 15, 16, 17. La mémoire 17 ZP est chargée à partir des données issues de la détermination de la zone polluée, extérieure au dispositif. La mémoire QT 16 est adressée par une double indexation : numéro de combinaison et numéro de fréquence de récurrence à l'intérieur de la combinaison. Pour simplifier le montage, le dénombrement des détections possibles est un peu différent de ce qui a été décrit plus haut. La mémoire 22 dont le contenu est ND, contient le nombre de repliements d'un quantum distance de l'étendue testée dans la zone polluée plutôt que l'inverse. La remise en forme est effectuée au moment de l'adressage de la mémoire RP 24 en soustrayant le contenu de ND à NFR au moyen du sommateur Z2 23. Le contenu de la mémoire ND 22 est envoyé vers le sommateur 3 puis remis à zéro (module R3) à chaque front d'horloge du compteur N. Le retard appliqué sur ce front permet au reste du dispositif en amont (modules R1 et R2) d'achever leurs calculs. La note est attribuée en fonction de la répartition obtenue. Elle représente l'espérance de la détection sur l'étendue distance du test. Le schéma d'un dispositif calculant cette note est donné sur la figure 9 la note obtenue est contenue dans un registre NT 26. où [CF(L)I e [O,NFRl sont des coefficients pondérateurs. Les coefficients CF(L) représentent la probabilité globale de détection d'une cible se trouvant dans un quantum distance tel que - L repliements sur NFR se font hors de la zone polluée (avec une probabilité élémentaire de détection évaluée à P); - NFR-L repliements sur NFR se font dans la zone polluée repliée (la probabilité élémentaire de détection vaut alors 0). La probabilité élémentaire p est celle d'une cible de SER de référence sur bruit thermique à PFA donnée. La comparaison de la note obtenue à la meilleure note obtenue précédemment et la mise en mémoire de l'indice de la meilleure note sont également décrites sur la figure 9. La note NT de la combinaison en cours d'évaluation est comparée à la note maximum NM rencontrée lors des évaluations précédentes au moyen d'un comparateur 27 recevant sur l'une de ces bornes la note NT en provenance du registre 26 et sur l'autre la note NM en provenance d'un registre 28. Si NT est plus grande que N M' la combinaison actuelle est meilleure que toutes les précédentes. L'indice K de la combinaison est donc mémorisé (dans un registre 29 NK) et la note NT est rangée dans le registre N M 28 à la place de la valeur antérieure. Le calcul de la note NT contenue dans le registre 26 est effectuée au moyen du dispositif représenté sur la partie gauche de la figure 9. Un compteur L 33 contenant NFR+1 valeurs, adresse des mémoires 30, 31 qui contiennent respectivement les valeurs des composantes du vecteur colonne RP et du coefficient CF(L). Les contenus de ces mémoires X,Y sont utilisés comme entrées par un module de calcul 32 pour effectuer le calcul de la note conformément à la formule choisie. Le chargement du registre NT 26 n'est autorisé sur un signal d'horloge en provenance du compteur 33 qu'après fin du calcul par le calculateur 32. Le procédé selon l'invention peut être utilisé dans toute circonstance où l'on désire que les repliements en distance d'une cible s'effectuent de manière privilégiée en dehors d'une certaine zone. Les conditions préalables à son utilisation ont été décrites précédemment. Il s'agit de déterminer - l'ensemble des combinaisons de fréquences de récurrence susceptibles d'être utilisées (satisfaisant aux conditions primaires du mode), - l'étendue distance du test, - la zone polluée. En mode recherche l'utilisation du procédé permet d'obtenir une détection "lissée", en évitant au maximum de perdre une cible dans des "trous" de détection. En mode poursuite, la mise en oeuvre du procédé s'effectue sur une étendue distance restreinte, correspondant à la distance de la cible poursuivie et permet d'assurer des conditions optimales de détection. Les courbes représentées figures 10 et 11 permettent de visualiser les performances que l'on peut attendre de l'exploitation du procédé. La figure 10a représente pour chaque quantum distance pour une combinaison donnée le nombre de détections possibles. Ce nombre varie entre 1 et 8. La figure îîa représente ce même nombre pour une autre combinaison de fréquence. La figure lOb, représente pour la même combinaison que la figure 10a l'espérance de détection en fonction de la distance. La figure îîb, représente cette espérance pour la même combinaison que la figure 11a. Ces courbes sont destinées à illustrer la corrélation existant entre le nombre élémentaire de détections possibles (ND) pour chaque quantum distance testée et la probabilité de détection d'une cible pour chacun de ces quantum. Les courbes des figures 10b et îlb correspondent à une combinaison ayant obtenu une meilleure note que celles correspondant aux courbes des figures 10a et 11a. Cela se traduit par le fait que la probabilité de détection représentée figure îîb est de façon plus régulière et plus large au dessus de 0,5. Ces courbes correspondent à une probabilité élémentaire de détection p = 0,32 et à un coefficient CF(L) = O pour ksi3. Cela signifie que les quanta distance obtenant moins de 4 détections possibles ne contribuent pas à l'augmentation de la note. REVENDICATIONS 1. Procédé applicable à un radar porté par un mobile, et destiné à choisir parmi NCB combinaisons possibles de NFR fréquences de récurrence chacune, la combinaison la meilleure pour détecter, à l'intérieur d'une zone distance prédéterminée dite zone testée composée de NCD quanta-distance, la présence d'une cible malgré la présence d'échos de sol se trouvant dans une zone dite polluée supposée connue et qui, pour chaque fréquence de récurrence, se replient à l'intérieur d'une zone dite zone polluée repliée, procédé caractérisé en ce que - on calcule pour chaque combinaison de fréquence de récurrence une note qui est une fonction croissante du nombre de fois où, pour les fréquences de récurrence entrant dans la combinaison les différents quanta-distance de la zone testée se replient hors de la zone polluée repliée et on choisit la combinaison ayant obtenu la note maximum. 2. Procédé selon la revendication 1, applicable à un radar porté par un mobile et destiné à choisir parmi NCB combinaisons possibles de NFR fréquences de récurrence chacune, la combinaison la meilleure pour détecter, à l'intérieur d'une zone distance prédéterminée dite zone testée composée de NCD quanta-distance, la présence d'une cible malgré la présence d'échos de sol se trouvant dans une zone dite polluée supposée connue et qui, pour chaque fréquence de récurrence, se replient à l'intérieur d'une zone dite zone polluée repliée, procédé caractérisé en ce que 1) pour la première fréquence de récurrence entrant dans la première comblnaison, on calcule, pour chaque quanta-distance de la zone testée, la distance repliée de ce quanta-distance 2) on compare la distance repliée de chaque quantum de la zone testée à la distance repliée de chacun des quanta-distance de la zone polluée repliée, et l'on attribue au quantum distance de la zone testée un coefficient E qui est q égal à 1 si le quantum testé se situe hors de la zone polluée 3) on recommence les opérations 1 et 2 pour chacune des NFR fréquences de récurrence entrant dans la première combinaison 4) on additionne les coefficients E obtenus pour q chacun des NCD quanta-distance testé, constituant ainsi un vecteur colonne "somme eq" dont les NCD composantes dont désignées par un nombre L, L pouvant prendre les valeurs entières de O à NFR 5) on constitue un vecteur colonne répartition désigné RP et comprenant NFR+1 composantes désignées RPoJ RP1, ..... .RPNFR > ces composantes ayant initialement la valeur O; 6) on incrémente la composante RPo du nombre de fois où apparaît dans le vecteur colonne constitué à l'étape 4, puis la composante RP1 du nombre de fois ou apparaît 1 dans le même vecteur colonne, et ainsi de suite Jusqu'à incrémenter RPNFR du nombre de fois où apparaît NFR dans le vecteur cotonne constitué à l'étape 4 7) on calcule une première note NT selon la formule dans laquelle CF est un coefficient qui est une fonction de L et ou RPL est la composante de rang L du vecteur R? 8) on recommence les étapes 1 à 7 pour la seconde combinaison calculant une seconde note NT2;; Si NT2 est supérieure à NT1 alors NT devient NT2, NT reste NT1 dans le cas contraire 9) on recommence les étapes 1 à 7 pour chacune des NCB-2 combinaisons restantes, calculant à chaque fois une note NT qui vient se substituer à la note maximum NTM calculée selon le processus décrit pour les étapes 1 à ; 10) on sélectionne la M ème combinaison, combinaison que l'on garde jusqu'à ce qu'intervienne un changement prédéterminé affectant le mobile. 3. Procédé selon la revendication 2, caractérisé en ce que le coefficient CF(L) est égal à L. 4. Procédé selon la revendication 2, caractérisé en ce que le coefficient CF(L) est égal à formule dans laquelle p est la probabilité élémentaire de détection du radar. 5. Dispositif destiné à réaliser le procédé selon l'une des revendication 1 à 3 du type comprenant des mémoires de stockage de données et des modules effectuant des opérations arthmétiques sur les données caractérisé en ce que les modules effectuant les opérations arthmétiques 18, 19, 32, reçoivent leurs données de mémoires 15, 16, 17, 30, 31, qui sont adressées par des compteurs 12, 13, 14, 33.