FR 2633725 A1 19900105 FR 8808920 A 19880701 DESCRIPTION "Système radar radioaltimétrique à onde continue modulée en fréquence" L'invention concerne un système radar du type à modulation de fréquence, fournissant à une antenne émettrice une onde continue dont la fréquence f varie de façon à peu près linéaire d'une fréquence inférieure fl à une fréquence supérieure f2 telles que f2-fl = Af, comportant un récepteur homodyne relié à une antenne réceptrice pour fournir un pre- mier signal de battement Fb1 entre ondes émise et reçue après réflexion, conçu pour la mesure d'altitude h par mesure de la rotation totale de phase (comptage du nombre de passages par zéro) dudit signal de battement Fb1 pendant la durée At de chaque excursion Af de la fréquence émise. Lorsqu'on réalise un radioaltimètre, destiné à étre embarqué sur un aéronef, on souhaite généralement pouvoir mesurer des altitudes selon une plage très étendue de distan- ces et, si possible, allant de l'altitude zéro jusqu'à plus de 10.000 m. Deux méthodes sont actuellement employées de façon habituelle pour réaliser un radioaltimètre radioélectrique. La première méthode consiste en la mesure directe du temps de propagation d'une impulsion radioélectrique ou d'un train d'impulsions codées selon un code par exemple de type PN. Ce type de matériel fonctionne bien aux altitudes élevées mais devient critique dès qu'il est nécessaire de mesurer des alti- tudes très basses à cause de la largeur non négligeable des impulsions et des différents vacillements. L'autre méthode, qui est utilisée selon la présente invention, consiste en la mesure indirecte du temps de propagation par l'intermédiaire de la fréquence de battement obtenue par la corrélation entre une onde émise, modulée linéairement en fréquence, et le si- gnal reçu du sol. Ce type de matériel, dit FM/CW, au moins lorsqu'il comporte deux antennes distinctes, une pour l'émis- sion l'autre pour la réception comme c'est le cas pour la pré- sente invention, est mieux adapté que le radar à impulsions indiqué cidessus à la mesure des basses altitudes mais pose des problèmes de précision insuffisante et surtout de limita- tion pour la mesure des altitudes élevées. Dans ce dernier cas, en effet, le signal d'écho reçu du sol est très affaibli et pour que le signal de battement soit toujours identifiable il faut habituellement que la linéarité de modulation du si- gnal émis soit très bonne et que le bruit de phase du signal de battement soit faible. Le Radar-radioaltimètre FM/CW est limité, aux grandes altitudes, à partir du moment o le bruit de la phase du signal de battement entre émission et réception provoqué par les bruits et autres non linéarités est de l'or- dre de 2U. Le signal reçu d'un écho de sol est la somme d'un grand nombre d'échos indépendants mais o domine la notion de phase liée au voisinage du pied de la droite perpendiculaire au sol passant par l'altimètre. Cette somme très fluctuante que l'on peut représenter par un vecteur OS peut se résoudre en la somme, dans le plan des phases, d'un vecteur moyen eA lentement variable et d'une somme aléatoire de vecteurs se dé- plaçant dans une zone d'autant plus grande, par rapport au mo- dule du vecteur 6A que, pour un matériel donné, l'altitude est grande. La partie principale du spectre de battement des alti- mètres à modulation de fréquence est liée au vecteur moyen gA. Soit un altimètre qui, à l'instant t1 auquel la fréquence émise est fl, est situé à la distance h du point moyen A cor- respondant à l'écho de sol principal, qui fournit le vecteur moyen UA précité. La différence de phase t1 entre le signal émis et le signal reçu peut s'écrire: À1 = 2U -c fl+? n 2h fy (1) c avec: c: vitesse de l'onde électromagnétique Y: angle fixe dépendant des circuits de l'altimètre. De mame à un instant t2 très proche de t1, pour la fréquence f2, l'altitude h n'ayant pas varié entre les instants t1 et t2 très rapprochés l'un par rapport à l'autre: 2 = 2h f2+ (2) c On déduit des relations (1) et (2) la loi des altimètres à mo- dulation de fréquence: 2h q2- 91 = 21 c (f2-fl) soit: Ag = 4X h (3) c La relation (3) est mieux connue sous sa forme dé- rivée par rapport au temps: fb = dt 2h = 2 c Af Att21[ c At fb = 2 h Af (4) cAt fb désignant la fréquence du signal de battement Fb précité. La relation (3) est directement utilisable dans les altimètres à comptage de zéros comme c'est le cas pour la présente invention. La relation (3) peut en effet s'écrire, en posant Ah = kÉ (k entier): h = k (5) De tels altimètres du type radar FM/CW sont connus notamment de l'ouvrage Radar Handbook, Skolnik, publié en 1970 par Mc Graw-Hill, chapitre 16-18 plus particulièrement page 16-29. Soit l'un de ces altimètres classiques dont les moyens de réception n'effectuent qu'un changement de fréquence, c'està-dire ne fournissent qu'un signal de battement Fb1. Ce signal Fbl peut être représenté par le vecteur OS précité (fi- gure 1). Si le vecteur d'oscillation locale e'O qui représente l'axe des abscisses dans le plan des phases a une amplitude supérieure à celle du vecteur OS, le signal de battement est alors pratiquement la projection du vecteur eS sur le vecteur d'oscillatio9 locale 6'O. Le point A, dans le plan des phases, décrit en substance un cercle et le point S décrit, par rap- port à ce cercle des ondulations aléatoires parasites, qui peuvent ou non constituer des boucles. A la traversée des axes de coordonnées du plan des phases ces ondulations peuvent se traduire par une multiplicité de zéros, à la traversée de cha- que axe, au lieu d'un zéro unique, plus précisément par une ou plusieurs paires de zéros parasites supplémentaires. Ceci se produit lorsque le signal de battement est entaché d'un fort bruit de phase. Un altimètre FM/CW à comptage de zéros tel qu'indiqué ci-dessus ne saurait pas différencier des paires de zéros parasites et indiquerait un nombre de zéros trop élevé, soit une altitude mesurée trop grande, conformément à la rela- tion (5). On connait aussi des radioaltimètres FM/CW à fré- quence de battement asservie. Dans ce dernier cas, la fréquen- ce du signal de battement est maintenue à une valeur sensible- ment constante et les ondulations aléatoires parasites préci- tées sont moins marquées, c'est-à-dire que ce type de matériel est plus robuste vis à vis des bruits, mais il devient à son tour inefficace si les erreurs de phases deviennent plus im- portantes, au point d'aboutir à la formation systématique de boucles pour ce qui est des ondulations aléatoires parasites, cas dans lequel l'égalité (4) qui sert à la détermination de l'altitude h n'est plus utilisable. Or, dans tous les cas, entre les instants tl et t2 auxquels les fréquences sont respectivement f1 et f2, quelle que soit la courbe de fréquence f(t) parcourue, les phases du vecteur eS existent et ne dépendent pas des valeurs intermé- diaires. Le but de l'invention est de mesurer avec préci- sion la différence de phase totale du signal de battement en- tre les instants t1 et t2 pour une modulation de fréquence émise, entre les fréquences f1 et f2z, qui n'est qu'approxima- tivement linéaire. Un autre but de l'invention est de pouvoir mesurer des altitudes élevées, au moyen d'un radar FM/CW à mesure de la rotation totale de phase. Ces buts sont atteints et les inconvénients de l'art antérieur sont atténués grâce au fait que le radar indiqué au premier paragraphe est remarquable en ce qu'il com- porte en outre des moyens dé réception pour fournir un deuxiè- me signal de battement Fb2 déphasé de R/2 par rapport audit premier signal Fb1 ainsi que des moyens de numérisation pour échantillonner les signaux Fb1 et Fb2 à une fréquence 1/T de l'ordre de plusieurs dizaines de fois celle des signaux Fbj et Fb2 et pour convertir ces derniers en signaux numériques bl(nT) et bz(nT) respectivement, des moyens de comparaison pour comparer, à chaque période T d'échantillonnage les signes de bl(nT), bl[(n-1)T], b2(nT), b2[(n-1)T] et des moyens de calcul pour calculer à partir desdites comparaisons des incré- ments et pour déduire, par incrémentations-décrémentations successives au cours de la durée At, le nombre de fois a, compté algébriquement, que le vecteur électrique eb dont Fb1 et Fb2 sont les composantes sur deux axes perpendiculaires a effectué de quarts de tour dans un sens trigonométrique donné, m étant proportionnel à l'altitude à mesurer, conformément à la relation: h mc (6) 8(f2-fi) L'idée de départ de l'invention est de suivre avec suffisamment de précision l'évolution du vecteur 3S dans le plan des phases de façon à pouvoir déterminer la variation de phase totale de ce vecteur entre deux valeurs de fréquence émises. Pour cela, les signaux de battement en quadrature fbj et fbz sont échantillonnés à une fréquence assez élevée pour que soient pris en compte tous les passages par zéro succes- sifs de l'un ou l'autre de ces signaux avec, pour chaque zéro, le sens trigonométrique de traversée de l'axe de coordonnée correspondant à chaque zéro, fourni par l'indication complé- mentaire du signe du signal de battement en quadrature qui ne s'annule pas. A cet égard, un mode de réalisation préféré de l'invention est remarquable en ce que lesdits moyens de compa- raison associent aux signaux bl(nT) et b2(nT) les fonctions BI(nT) et B2(nT) respectivement qui valent +1 ou -1 selon que les signaux respectifs b1(nT) ou bz2(nT) sont positifs ou néga- tifs, et qu'à chaque période d'échantillonnage T lesdits moyens de calcul effectuent les calculs suivants: BI(nT).Bi[(n-1)T],= Cl(n,n-1) (7) B2(nT). Bz2(n-1)T] = Cz(n,n-1) (8) si CI = -1 alors: B1(nT).B2(nT)+m * m (9) si C2 = -1 alors: -BI(nT).B2(nT)+m * m (10) De cette façon, les rotations de phase de +t/2 puis de -n/2 précitées dues aux ondulations aléatoires parasi- tes du vecteur 5S s'annulent l'une l'autre et la rotation to- tale de phase à U/2 près qu'on peut qualifier d'utile est pri- se en compte en totalité. La description qui suit en regard des dessins an- nexés, le tout donné à titre d'exemple, fera bien comprendre comment l'invention peut être réalisée. La figure 1 représente le signal d'écho sous sa forme vectorielle OS dans le plan des phases OL1, OL2. La figure 2 représente l'évolution du vecteur S dans le temps pour une ondulation aléatoire parasite peu mar- quée. Les figures 3a et 3b représentent, en fonction du temps, les projections du vecteur 3s de la figure 2 sur les axes OL1 et OL2 respectivement. La figure 4 représente l'évolution du vecteur dS dans le temps pour une ondulation aléatoire parasite assez marquée. Les figures 5a et 5b représentent, en fonction du temps, les projections du vecteur 5S de la figure 4 sur les axes OL1 et OL2 respectivement. La figure 6 est le schéma synoptique d'un mode de réalisation du sytème radar radioaltimétrique selon l'inven- tion. La figure 7 représente l'implantation d'une ligne à retard sur le radar de la figure 6. Sur la figure 1, on a représenté, avec interrup- tion, en 5'O le vecteur d'oscillation locale qui symbolise l'onde d'émissibn d'un radar radioaltimétriqiue FM/CW. Le vec- teur 6'0 est fixe et confondu avec l'axe des abscisses du plan des phases OL1, l'axe des ordonnées étant OL2. Dans le plan OL1, OL2 on a représenté le signal reçu d'un écho de sol sous la forme du vecteur OS en général petit devant 5'O qui est la somme d'un vecteur moyen &A lentement variable que l'on cher- che à identifier et d'une somme aléatoire de vecteurs 1, 2, 3, 4, 5, qui se déplacent dans une zone Z délimitée par un trait interrompu pour une valeur donnée du temps t, d'autant plus grande, par rapport au module du vecteur OA que, pour un maté- riel donné, l'altitude est grande. Sur la figure 2, on a représenté l'hodographe H2 du vecteur OS pour un radioaltimètre FM/CW dans lequel la fré- quence émise suit une loi à peu près linéaire, de l'instant to (fréquence fo) à l'instant t31 (fréquence f31), les instants successifs to à t31 étant représentés par leurs seuls indices 0, 1, 2,...,-31. En l'absence de bruit de phase, représenté par la zone Z sur la figure 1, la courbe H2 serait celle qui est décrite par le point A, c'est-à-dire sensiblement un cer- cle de centre O, non représenté. Il existe cependant des ondu- lations parasites, ces ondulations étant peu marquées sur la figure 2 en ce sens qu'il n'y a pas formation de boucles pour la courbe 52. Ceci signifie que le bruit, assimilable au vec- teur Is non représenté, est faible et que la probabilité d'une multiplicité de zéros à la traversée des axes OL1, OL2 est faible. Sur la figure 2, on notera deux paires de zéros sup- plémentaires parasites, l'une à la travesée de la partie néga- tive de l'axe OL1 entre les instants t4 et ts et l'autre à la traversée de la partie négative de l'axe OL2 entre les ins- tants to10 et t13, comme cela est bien visible sur les figures 3a et 3b o les indices qui repèrent les instants to à t31 sont indiqués en abscisses, régulièrement espacés. On notera que méme dans le cas d'un signal de battement peu bruité comme représenté aux figures 2 et 3a, 3b, un altimètre classique à comptage de zéros indiquerait des zéros en excès et donc une distance h trop grande, conformément à la relation (5) ci-des- sus. A titre de référence, on a représenté par la sinusoïde en trait interrompu 6, à la figure 3a, le signal de battement théorique pur, non bruité, qui est la projection du vecteur 5A sur l'axe OL1. Les figures 4 et 5a, 5b représentent le cas o les erreurs de phase sont grandes et provoquent des ondula- tions plus profondes, allant jusqu'à la formation de boucles sur l'hodographe de la figure 4. La formation de boucles si- gnifie que l'amplitude du vecteur de bruit IS devient compara- ble à celle du vecteur UA et que sa pulsation est plus élevée, environ dix fois plus grande dans l'exemple de la figure 4. La présence de boucles est favorable à la formation de paires de faux zéros. Sur la figure 4 ou sur les figures 5a et 5b, on dénombre 8 zéros, à la traversée de l'axe OL1, au lieu des 2 zéros qui doivent normalement mesurer une rotation de phase de 2U du signal de battement entre les instants to (fréquence fo) et t3s (fréquence f3s), si l'on considère la traversée du seul axe OL1. Dans les altimètres à comptage de zéros connus, on n'utilise qu'un seul signal de battement Fb1 représenté aux figures 3a et 5a et le comptage de zéros ne permet pas d'iden- tifier les zéros parasites indiqués ci-dessus. Par contre, l'invention se propose de démoduler en même temps l'écho avec une oscillation en quadrature, OL2 qui fournit le signal re- présenté aux figures 3b et 5b, ce qui conduit à compter les zéros à la traversée des deux axes'à raison de 4 par tour dans le plan OL1, OL2. En combinaison avec la caractéristique tech- nique indiquée à la phase précédente, les calculs et le comp- tage nécessaires à la comptabilisation des zéros utiles sont effectués à partir de valeurs numériques des signaux de batte- ment en quadrature b1 et b2 précités, tout en effectuant une corrélation entre les valeurs homologues dans le temps de b1 et b2. Pour que tous les zéros soient pris en compte, il con- vient d'échantillonner les valeurs des signaux bl et bz à une fréquence suffisamment élevée et au moins égale au double de la fréquence de rotation du vecteur de bruit AS. En pratique, comme cela apparait sur les figures 2 et 4, la fréquence d'échantillonnage est au moins égale à trente fois la fréquen- ce des signaux de battement fb1 et fb2. Sur les figures 1 et 4 on peut considérer les repères 1 à 31, respectivement 1 à 39, comme les instants d'échantillonnage à des instants multiples de T, soit T, 2T,..., 31T, respectivemnet T, 2T,..., 39T. Un mode de réalisation préféré du radar selon l'invention est représenté à la figure 6. De façon connue, un oscillateur commandé en tension (VCO) 11 émet un signal hyper- fréquence qui alimente, par l'intermédiaire d'un coupleur 12 à la fois l'étage de puissance d'émission (émetteur) 13 relié à l'antenne émettrice 14 et un premier mélangeur 15. L'émetteur 13 est de préférence relié à l'antenne d'émission 14 par l'in- termédiaire d'un isolateur, non représenté, dans le cas o le rapport d'onde stationnaire de l'antenne 14 est trop élevé et risque de perturber le système; il en est de même pour l'an- tenne de réception 16, reliée à un étage de traitement 17 qui, selon les besoins, peut remplir la triple fonction de filtre hyperfréquence, utile pour des altimètres en ambiance hostile, de limiteur de puissance, cette contrainte pouvant être impo- sée par exemple par des radars puissants et proches, et d'am- plificateur hyperfréquence. L'étage 17 a, de façon connue, une sortie 18 reliée à une deuxième entrée du mélangeur 15. Selon l'invention, le radar de la figure 6 comporte un deuxième mé- langeur 19 dont une première entrée est reliée à une deuxième sortie du coupleur 12 et une deuxième entrée à une sortie de l'étage 17. Les déphasages entre les signaux d'entrée des mé- langeurs 15 et 19 sont conçus et les longueurs électriques des conducteurs reliés à ces entrées sont ajustées de telle façon que les signaux de battement Fb1 et Fb2 en sortie des mélan- geurs 15 et 19 soient déphasés de 1/2. Le mélangeur de 1/2 est par exemple engendré à patir du coupleur 12 dont les sorties d'oscillation locale dérivées engendrent par rapport au signal d'émission un déphasage nul, respectivement égal à 1/2 vers le mélangeur 15, respectivement 19, comme représenté sur la figu- re. Le même signal de sortie, ou des signaux de sortie de l'étage 17 différant en phase comme représenté sur la figure 6, peuvent être appliqués aux deuxièmes entrées des mélangeurs 15 et 19. Chaque mélangeur 15, 19 est suivi d'un amplificateur basse fréquence 21, 22 dont la bande passante couvre environ une octave, étant par exemple comprise entre 50 et 100 kHz, et dont la fréquence haute doit être compatible avec la vitesse de l'échantillonnage effectué en aval, c'est-à-dire au plus égale à la moitié de la fréquence d'échantillonnage fs. La dernière contrainte indiquée est facile à satisfaire dans la mesure o il existe une contrainte plus sévère sur la fréquen- ce fs déjà indiquée plus haut, selon laquelle cette fréquen- ce d'échantillonnage est égale à plusieurs dizaines de fois la fréquence fb. Les amplificateurs 21 et 22, aussi identiques que possible, ont un gain variable suivant une dynamique assez importante, de l'ordre de 80 dB, fonction de la mission de l'altimètre et sont munis d'une commande automatique de gain qui peut être commune, comme représenté, ou séparée. Chaque amplificateur 21, 22 est relié à un échantillonneur-bloqueur 23, respectivement 24, suivi d'un convertisseur analogique-nu- mérique (CAN) 25, respectivement 26, les quatre éléments 23 à 26 constituant des moyens de numérisation. Les CAN 25 et 26 ont une dynamique qui va par exemple de +128 à -128, soit 8 bits. On notera cependant que, pour la mesure de l'altitude, selon l'invention, une dynamique de plus ou moins 1 est suffi- sante, ce qui revient à se limiter à la seule connaissance du signe, comme on le verra ci-dessous. Les CAN 25 et 26 sont re- liés à un microprocesseur (MUP) 27. Le microprocesseur fournit deux signaux numériques sur des conducteurs multiples, un si- gnal CG à un convertisseur numérique-analogique (CNA) 28 et un signal CE à un convertisseur numérique-analogique 29. Le si- gnal CG, après conversion en signal analogique en 28 commande automatiquement le gain des amplificateurs 21 et 22. L'ampli- tude moyenne des signaux bl(nT) et bz(nT) en sortie des CAN 25 il 11. et 26 est comparée à des seuils dans le microprocesseur 27 et, lorsque cette amplitude devient trop faible ou trop forte, une correction compensatrice est apportée à l'amplification en 21 et en 22, par le signal CG, ce qui constitue une boucle d'as- servissement du gain des amplificateurs 21 et 22. Des varia- tions de gain par sauts de 6 dB conviennent, ainsi qu'un con- vertisseur 28 & 4 bits. Le signal numérique CE commande quant à lui l'oscillateur 11 par l'intermédiaire du montage en série du CNA 29, d'un lisseur de courbe 31 et d'un organe de discri- mination de fréquence 32. Le signal de sortie du CNA 29 a une forme en marches d'escalier, ce qui est défavorable à une com- mande correcte de l'oscillateur 11. En effet, pour que la me- sure de l'altitude par comptage de zéros ait un sens, il faut éviter toute discontinuité dans la fréquence émise, qui provo- querait à son tour une rupture de continuité de mesure de la phase du signal de battement. Le lisseur de courbe 31 dont l'élément essentiel est par exemple un condensateur 33 a pour fonction de supprimer au premier ordre les variations brusques du signal de sortie du CNA 29 qui seraient assimilables à des discontinuités. D'autre part on peut prévoir que si la bande passante des amplificateurs 21 et 22 est d'une octave, la ban- de de fréquence du signal altimétrique qui les traverse peut subir des variations brusques d'une demi' octave, ces considé- rations étant utiles pour la conception du CNA 29. Si l'oscil- lateur 11 est à peu près linéaire en fréquence en fonction de la tension de commande, alors la dynamique du CNA 29 devra être égale à la dynamique des altitudes mesurées, y compris la longueur résiduelle des cables. En effet, le condensateur 33 crée une intégration du courant et la fréquence de battement altimétrique fb est, à une constante multiplicative près, le produit de la pente de variation de fréquence et de l'altitu- de, comme l'indique la relation (4) ci-dessus. Si l'on veut par exemple mesurer des altitudes variant de 1 m à 10.000 m, la dynamique du CNA 29 devra être d'environ 27 bits si ce con- vertisseur est linéaire mais seulement de 5 bits si ce conver- tisseur numérique-analogique est à caractéristique exponen- tielle. L'organe 32, sur la figure 6, symbolise des moyens pour repérer les instants tl et t2 d'émission des fréquences initiale f1 et finale f2 par l'oscillateur 11, ces valeurs de fréquence déterminant entièrement la valeur de AD recherchée et donc la mesure de h en vertu de la relation (3) ci-dessus qu'on peut encore écrire: h =f 41 () L'organe 32 peut se réduire à une simple connexion 34 dans le cas o l'oscillateur commandé par tension 11 est fidèle; il suffit alors de repérer les tensions V1 et V2 aptes à engen- drer les fréquences fl et f2 en sortie du VCO11, ce qui est effectué par un conducteur de retour 35 vers le microproces- seur 27; ceci permet de déterminer, dans le microprocesseur 27, les instants ti et t2 marquant le début et la fin du comp- tage de zéros. Si l'on souhaite une plus grande précision ou si l'oscillateur 11 n'est pas suffisamment fiable on peut, se- lon une variante non représentée, disposer, en sortie du VC011, des résonateurs hyperfréquences, reliés au microproces- seur 27 et indiquant à ce dernier directement les valeurs de fréquence émise. On notera qu'en toute rigueur un retour n'est pas indispensable vers le microprocesseur 27, pour la détermi- nation des instants t1 et t2; ce retour serait en effet inu- tile dans le cas de l'utilisation d'un VCO très reproductible, en fabrication, est très fidèle, pour lequel aux fréquences fl et f2 correspondraient des valeurs de tension de commande V1 et V2 dont les instants d'émission tl et t2 en provenance du microprocesseur 27 sont faciles à déterminer. Cependant un os- cillateur présentant de telles qualités est cher et d'autre part, un retour du signal émis vers le microprocesseur est toujours nécessaire pour pouvoir à tout instant ajuster de-fa- çon connue la pente de modulation de ce signal en fonction de l'altitude mesurée, ceci afin de maintenir les fréquences des signaux de battement Fb1 et Fb2 à l'intérieur des bandes pas- santes des amplificateurs 21 et 22. De préférence, la dent de scie d'émission est triangulaire symétrique ce qui, par effet de moyenne entre rampe croissante et rampe décroissante, permet de s'affranchir de façon connue d'un effet Doppler éventuel. Ainsi, lorsque la fréquence haute f2 est atteinte, la tension V2 sur le conduc- teur 35 a pour effet d'inverser, en 27, le signe de la pente du signal logique CE, cette pente devenant négative. Lorsque la valeur V1 est atteinte, en 34, cette pente redevient à nou- veau positive et ainsi de suite. La mesure de l'altitude est effectuée, conformé- ment à la relation (3) ci-dessus, en comptant le nombre de tours effectué par le vecteur OS dans le plan OL1, OL2, pen- dant que la fréquence varie de façon continue de fl à f2- Se- lon l'invention, le nombre de tours sera le résultat du comp- tage du quart du nombre de fois m que le point S aura traversé le demiaxe OL1 positif, le demi-axe OL2 positif, le demi-axe OL1 négatif et le demi-axe OL2 négatif en comptant +1 dans le sens trigonométrique direct et -1 dans le sens rétrograde. On a la relation A m. 2 La relation (3) peut alors s'écrire h = M.C (6) 8(f2-fl) 6) comme déjà indiqué ci-dessus, la relation (6) étant très voi- sine de la relation (5) plus classique aussi indiquée ci-des- sus, selon laquelle on compte des demi-tours et non des quarts de tour. En pratique, la différence: f2-fl, est de l'ordre de 120 MHz et l'on peut écrire: h = 0,3125 m en mètres (13) Pour identifier la traversée d'un demi-axe et le sens de travesée on compare, à chaque période T d'échantillon- nage les signes de bl(nT), b1[(n-1)T], b2(nT), b2[(n-1)T], les mémorisations et les corrélations nécessaires étant effectuées dans des moyens de comparaison 36 internes au microprocesseur 27. Des moyens de calcul 37 permettent, à partir des comparai- sons précédentes de calculer, par incrémentations-décrémenta- tions successives au cours de la durée: t2-t1, le nombre m. De préférence, les moyens de comparaison associent aux signaux bl(nT) et b2(nT) les fonctions B1(nT) et B2(nT) respectivement qui valent +1 ou -1 selon que les signaux res- pectifs bl(nT) ou b2(nT) sont positifs ou négatifs, et à cha- que période d'échantillonnage T, les moyens de calcul 37 ef- fectuent les calculs suivants: B](nT).BI[(n-1)T] = Cl(n,n-1) B2(nT).B2[(n1)T] = C2(n,n-1) Si C1 = -1 alors la valeur courante de m est incrémentée de la valeur algébrique: B(nT).B2(nT), ce qu'on écrit: BI(nT).Bz(nT)+m 4 m Si C2 = -1 alors la valeur courante de m est incrémentée de la valeur: -Bl(nT).B2(nT), soit: -BI(nT).Bz(nT)+m * m Le cas le plus fréquent pour les paires d'échan- tillons successifs bl et b2 est: CI = C2 = +1, ce qui dénote l'absence de traversée d'un axe et donne une valeur inchangée pour m. Ces calculs, comparaisons et incrémentations sont ef- fectués par le microprocesseur 27 entre les instants t1 et t2, la valeur obtenue en t2, pour la fréquence d'émission fz étant la valeur à prendre en compte pour l'application des relations (6) ou (13), ces derniers calculs étant effectués, par exemple à chaque dent de scie d'émission, par le microprocesseur 27. Pour une dent de scie d'émission symétrique destinée à éviter les discontinuités du signal émis et permettant de s'affran- * chir de l'effet Doppler, la valeur de a à introduire dans les formules (6) ou (13) est la valeur moyenne entre les valeurs absolues de m obtenues pour la rampe positive et pour la rampe négative. En effet, le sens de rotation du vecteur 8S dépend du sens de modulation du VC011. Or, c'est le même microproces- seur 27 qui fait le comptage des tours dans le plan des phases du vecteur 0S et qui produit le signal de commande du VCO, à savoir le signal logique CE en sortie du microprocesseur. On peut donc, de préférence, pour éviter les chocs électriques liés à un retour brusque de f2 à f1, effectuer une descente lente de la fréquence de f2 à f, avec un courant de décharge identique au courant de charge en inversant les signes de m. Pour que le traitement de signal et les calculs indiqués ci-dessus se réalisent correctement, il est nécessai- re que les signaux alternatifs de battement Fbj et Fb2 traver- sent les amplificateurs analogiques basse fréquence 21 et 22. On choisit de préférence pour ces amplificateurs une bande passante assez large, d'une octave, pour ne pas être gêné par les bruits ou la linéarité de vobulation. Il n'est pas souhai- table d'aller au-delà, en largeur de bande, car cela augmente- rait inutilement le bruit sur le signal utile. La fréquence moyenne du signal de battement est fb 4(tt) (14) 4(t2-tl) relation qui se déduit des relations suivantes déjà indiquées ci-dessus: = m et fb = - - _ (t2-tl)' 2 Le microprocesseur 27 commande la pente du VC011 de la façon suivante Si la fréquence fb, qu'il est facile de calculer avec précison en 27 à chaque dent de scie, est supérieure à la fréquence centrale fc des amplificateurs 21 et 22, la valeur du courant de charge du VC011 sera diminuée, à partir du si- gnal CE, dans un rapport prédéterminé et inversement augmentée de façon analogue, si fb la pente peut être fixée a priori, la gamme d'altitudes mesu- rée étant égale à la bande passante des circuits 21 et 22 soit, lorsque la bande passante des circuits est égale à une octave h maximum = 2h minimum. Pour limiter la dynamique nécessaire pour les CAN et 26 et le risque d'accrochage sur de faux signaux, le mi- croprocesseur 27, par l'intermédiaire du CNA 28, commande, en phase de recherche par l'intermédiaire du signal CG, le gain maximum acceptable GM des amplificateurs 21 et 22 en utili- sant, à défaut d'autre information, la relation qui lie le courant de charge i du VC011 aux autres variables: (AID+h) = K (15) 1. AID représentant la longueur résiduelle des cables, les constantes AID et K dépendant des antennes de l'altimètre pour le sol le moins réflecteur prévu au cahier des charges et compte tenu du fait que le gain GM recherché doit être pro- portionnel à h2. Dès qu'un signal identifiable est présent à la sortie des amplificateurs 21 et 22, le microprocesseur 27 commande un gain identique, sur les deux amplificateurs 21 et 22, et tel que le signal le plus fort soit de l'ordre de la moitié de la capacité des convertisseurs 25 et 26. Selon l'invention le microprocesseur 27 est utili- sé, de préférence, pour commander la pente de vobulation du signal émis; ceci permet une recherche de l'altitude au dé- marrage beaucoup plus rapide que selon le procédé classique. Dans l'art connu, on effectue une série de balayages complets, allant de fl à f2 avec un courant i décroissant dans le rap- port 1/12 à chaque nouveau balayage jusqu'à ce qu'on identifie une fréquence comptable passant par l'amplificateur basse fré- quence 21. Le microprocesseur 27 permet une recherche o le temps de comptage est constant et par exemple égal à 3,6 ms., ce qui correspond à un comptage m = 1000 au centre de l'ampli- ficateur basse fréquence 21 supposé avoir une bande de fré- quence comprise entre 50 et 100 kHz. Toutes les 3,6 ms, la pente est divisée dans le rapport 1/12 et le compteur de m re- mis à zéro. Ce processus est continué jusqu'à ce qu'on obtien- ne un signal identifiable. La dynamique de recherche d'altitu- de, supposée de 10.000 est alors explorée environ 1400 fois plus vite que de la façon classique indiquée ci-dessus et cet- te recherche ne dure que 0,1 s environ. Outre la précision sur l'altitude que procure l'invention par mesure continue de la phase comme indiqué ci- dessus, il est possible d'améliorer encore cette précision, en réduisant l'incertitude du comptage au début et à la fin de la mesure de m. Supposons que l'on désire un altimètre ayant une précision de 0,3 m ou 1 %o (un pour mille) selon ce qui est le meilleur. L'incertitude du comptage est au plus de 1 au début et à la fin de la mesure. Pour obtenir une précision de 1 %o sur une mesure, il faut donc compter m = 2000. Avec les hypo- thèses de: f2-fl = 120 MHz, et fb = 70 kHz on otbient: h = 2000 x 0,3125 = 625 m et ht t =t2-t 700 = 0,00714 s Si on accepte d'obtenir la mesure de l'altitude avec un léger retard, par exemple 0,1 s, on notera qu'il est possible d'effectuer, avec les hypothèses ci-dessus, 14 mesu- res indépendantes pendant ce délai ce qui permet, en faisant la moyenne mm de ces 14 mesures de ramener l'erreur probable de 1 %0 à 0,27 %0. Aux altitudes inférieures à 625 m, l'erreur sera supérieure en valeur relative mais le nombre de mesures sera plus important et donc la précision liée à la méthode de mesure sera meilleure. D'autre part, si l'on souhaite mesurer des altitu- des élevées, de l'ordre de 10.000 m, il convient de choisir, pour la vobulation, une fréquence haute f3, nettement moins élevée que la fréquence fz qui conduit à la valeur: f2-f1 = Af = 120 MHz indiquée cidessus. En effet les calculs, dans le microprocesseur, sont limités à un certain nombre de bits, par exemple 12 bits. La valeur maximale comptable est alors elle- même limitée à 4096 et la valeur de af, qui se déduit immédia- tement de la formule (6) ci-dessus, s'obtient en donnant à m et à h leurs valeurs maximales: Af = f3-fi = 8 = 15,36 MHz A un incrément de 1 pour m correspond un pas alti- métrique de 2,44 m. On peut donc mesurer une première plage d'altitudes comprises entre O m et 1280 m avec une valeur AF = f2-fl = 120 MHz puis une deuxième plage d'altitudes comprises entre 1280 a et 10.000 m avec une valeur AF = f3-fl = ,36 MHz. Pour ces mesures d'altitude, la plus mauvaise pré- cision est obtenue au moment de la commutation, vers 1281 m; le comptage m est alors de 525 en un temps de 0,002 s, soit 54 échantillons indépendants en 0,1 s; l'erreur probable de 2/525 sera donc divisée par: 154 = 7,3 et ramenée & 0,52 %o environ, ce qui est acceptable. On peut ainsi simplifier la tâche du microprocesseur en choisissant des repères de fré- quence intermédiaires, tels f3, aux altitudes élevées. Pour certains types de sols l'écho reçu est très fortement modulé en amplitude ou momentanément absent pendant des périodes plus ou moins longues. Des disparitions de signal d'écho se produisent notamment pour un vol stationnaire au- dessus d'une étendue d'eau. Selon une variante, le micropro- cesseur 27 contrôle à tout instant la somme des amplitudes ab- solues des deux signaux b1(nT) et b2(nT), soit: Ibl(nT)i + Ib2(nT)I. Dès que cette somme est inférieure à un niveau S jugé comme le minimum acceptable en fonction du matériel uti- lisé, ce minimum étant supérieur au bruit de l'amplificateur 21, le microprocesseur cesse d'accroître a et évalue le pour- centage du temps d'horloge pendant lequel cet ou ces arrêts auront lieu entre fl et f2. Si ce pourcentage est élevé, supérieur à 10 % par exemple, le microprocesseur annule cette mesure. Si ce pourcentage est faible, inférieur à 10 % par exemple, mais répétitif à toutes les mesures, le microproces- seur calculera la valeur de h par extrapolation en supposant que, pendant les périodes d'absence de signal d'écho, m a va- rié par exemple linéairement en fonction du temps ou selon toute autre loi prédéterminée. Si enfin le pourcentage est élevé et répétitif, le microprocesseur peut déclencher un signal d'alarme tout en continuant de fournir une estimation assez peu précise de l'altitude. Pour en revenir à la précision du radioaltimètre, on notera que cette dernière dépend principalement de l'esti- *ation des fréquences fl, f2 et f3. Pour s'affranchir de l'im- précision et des variations qui peuvent affecter les valeurs fl, f2 et f3 on peut, selon une variante de l'invention équi- per le radioaltimètre d'une ligne à retard interne, comme re- présenté à la figure 7. La ligne à retard 41 est reliée par une extrémité au conducteur d'antenne réceptrice 16 et par l'extrémité opposée au conducteur d'antenne émettrice 14 par l'intermédiaire du montage en série d'un interrupteur 42 et d'un coupleur 43. La fermeture du commutateur 42 est commandée périodiquement à partir du microprocesseur 27, par exemple à raison d'une fois par seconde. On suppose que les valeurs de fl, f2 et f3 sont données de façon approchée mais répétitive à court terme par des seuils de tension à l'entrée de l'oscilla- teur 11. Soit ho la valeur connue et stable de la ligne à re- tard de référence 41; on peut écrire, en vertu de la formule (12): ho = (16) et, en combinant les formules (6) et (16): h =. ho- mo Il reste dès lors à obtenir la valeur de mo avec une précision suffisante. Il peut se faire que la longueur ho de la ligne à retard 41 ne soit pas suffisante pour que l'erreur sur mo soit négligeable. Dans ce cas on peut donner à mo une valeur non entière formée de la partie principale entière comptée et de l'estimation temporelle, ramenée à des fractions de i/2, de la position de la commutation de tension fl ou f2 par rapport aux comptages de mo immédiatement avant et immédiatement après la Z633725 commutation de tension. Cette opération est possible sur le signal de la ligne à retard parce qu'il est ponctuel et non fluctuant. La valeur de mo étant connue à très peu près, tout écart trop important par rapport à cette valeur doit actionner un indicateur signalant que l'altimètre n'est pas en état de bon fonctionnement. REVENDICATIONS: 1. Système radar du type à modulation de fréquence, fournissant à une antenne émettrice une onde continue dont la fréquence f varie de façon à peu près linéaire d'une fréquence inférieure fl à une fréquence supérieure f2 telles que f2-fl = Af, comportant un récepteur homodyne relié à une antenne ré- ceptrice pour fournir un premier signal de battement Fb1 entre ondes émise et reçue après réflexion, conçu pour la mesure d'altitude h par mesure de la rotation totale de phase dudit signal de battement Fbj pendant la durée At de chaque excur- sion Af de la fréquence émise, caractérisé en ce qu'il compor- te en outre des moyens de réception pour fournir un deuxième signal de battement Fb2 déphasé de 1/2 par rapport audit pre- mier signal Fb1 ainsi que des moyens de numérisation pour échantillonner les signaux Fb1 et Fb2 à une fréquence 1/T de l'ordre de plusieurs dizaines de fois celle des signaux Fb1 et Fb2 et pour convertir ces derniers en signaux numériques bl(nT) et b2(nT) respectivement, des moyens de comparaison pour comparer, à chaque période T d'échantillonnage les signes de bl(nT), bl[(n-1)T], b2(nT), b2[(n-1)T] et des moyens de calcul pour calculer à partir desdites comparaisons des incré- ments et pour déduire, par incrémentations-décrémentations successives au cours de la durée At, le nombre de fois m, compté algébriquement, que le vecteur électrique Fb dont Fbj et Fb2 sont les composantes sur deux axes perpendiculaires a effectué de quarts de tour dans un sens trigonométrique donné, m étant proportionnel à l'altitude à mesurer, conformément à la relation: mc h = c 2. Système radar conçu pour la mesure d'altitude selon la reyendication 1, caractérisé en ce que lesdits moyens de comparaison associent aux signaux bl(nT) et b2(nT) les fonc- tions Bi(nT) et Bz(nT) respectivement qui valent +1 ou -1 se- lon que les signaux respectifs bl(nT) ou b2(nT) sont positifs ou négatifs, et qu'à chaque période d'échantillonnage T les- dits moyens de calcul effectuent les calculs suivants: BI(nT).BI[(n-1)T] = Cl(n,n-1) B2(nT).B2[(n-1)T] = C2(n,n-1) si CI = -1, alors BI(nT).Bz(nT)+ m 4 m si C2 = -1, alors -Bl(nT).B2(nT)+m 4 m 3. Système radar radioaltimétrique selon la revendica- tion 1 ou 2, caractérisé en ce que lesdits moyens de comparai- son et lesdits moyens de calcul font partie intégrante d'un microprocesseur. 4. Système radar radioaltimétrique selon la revendica- tion 3, caractérisé en ce que lesdits moyens de calcul sont conçus pour calculer une valeur non entière de m en effectuant la moyenne ma entre plusieurs valeurs entières successives mesurées de a. 5. Système radar radioaltimétrique selon la revendica- tion 3 ou 4, caractérisé en ce que pour la mesure des altitu- des appartenant à la plage d'altitudes les plus hautes, une troisième fréquence f3 du signal émis est prise en compte com- me fréquence supérieure par le microprocesseur, la valeur de f3 étant comprise entre fl et f2. 6. Système radar radioaltimétrique selon l'une des re- vendications 3 à 5, caractérisé en ce que ledits moyens de calcul sont conçus pour contrôler à tout instant la somme: Ibî(nT)I+Ib2(nT)I et, lorsque cette somme est inférieure à un seuil S pendant la mesure de m, pour évaluer le pourcentage du temps d'horloge pendant lequel ce seuil S est franchi vers le bas entre les fréquences d'émission fl et f2 (respectivement f3) et pour en déduire soit l'annulation de la mesure de m, soit une valeur approximative de m. 7. Système radar radioaltimétrique selon l'une des re- vendications 3 à 6, caractérisé en ce qu'il comporte en outre une ligne à retard interne de longueur équivalente ho qui est commutée périodiquement à partir du microprocesseur et pour laquelle le comptage de zéros afférent mg est effectué pério- diquement, l'altitude h étant calculée par lesdits moyens de calcul conformément à l'expression: h = ho 8. Système radar radioaltimétrique selon l'unedes revendications 2 à 7, caractérisé en ce que la modulation de fréquence est successivement croissante puis décroissante pour mettre en évidence et mesurer ou annuler en moyenne l'erreur causée par un éventuel effetDoppler.