FR 2898732 A1 20070921 FR 0602381 A 20060317 Procédé de compensation des erreurs de positionnement des éléments rayonnants d'une antenne réseau DOMAINE DE L'INVENTION L'invention se rapporte au domaine des antennes réseaux, composées d'un ensemble de sources individuelles réparties sur la surface de l'antenne et dont le positionnement mécanique est déterminé de façon à obtenir le diagramme de rayonnement voulu. CONTEXTE DE L'INVENTION - ART ANTERIEUR Lorsque l'on réalise une antenne à sources réparties, telle qu'une antenne réseau par exemple, un problème important à résoudre est constitué par la précision du positionnement mécanique relatif des différentes sources. Afin de détecter avec précision la présence de cibles dans l'espace d'observation, les antennes réseau ont besoin de recevoir le moins de puissance possible en dehors du lobe principal. Pour cela, les sources doivent être positionnées régulièrement et précisément sur l'antenne. Les mesures de diagramme de rayonnement F(0), les mesures en champ proche en particulier, montrent que les erreurs mécaniques des antennes réseau, erreurs de positionnement des sources, sont des erreurs qui dégradent sensiblement les diagrammes de rayonnement. Cependant, les techniques actuelles de montage, et en particulier celles qui permettent un montage en série, plus ou moins automatisé, ne permettent pas d'obtenir une précision de positionnement suffisante pour que le diagramme de rayonnement ne soit pas significativement perturbé. Cette perturbation se traduit par l'apparition, de chaque côté du lobe principal, d'une quantité non-négligeable de lobes diffus dont le gain est suffisant pour altérer notablement les performances de l'antenne, notamment en terme de directivité. PRESENTATION DE L'INVENTION Un but de l'invention est de rendre possible, de manière pratique, la prise en compte de ces erreurs de positionnement de façon à en minimiser les effets sur les performances générales de l'antenne.30 A cet effet l'invention a pour objet un procédé de compensation des erreurs de positionnement des éléments rayonnants d'une antenne réseau, caractérisé en ce qu'il comporte au moins: - une étape de mesure de l'erreur de positionnement de chaque élément rayonnant n et de constitution d'un vecteur d'imperfections I, dont les composantes sont les valeurs de la fonction 8(n) d'erreur de positionnement , - une étape de modélisation de la fonction d'erreur 8(n) sous forme d'une somme de composantes spectrales sinusoïdale En(i) calculées par application d'une transformée de Fourier au vecteur d'imperfections I, suivie de l'application d'une transformée de Fourier inverse au vecteur dd ainsi obtenu, - une étape de calcul de fonctions de corrections Aa(n), A(p(n), par modulation d'amplitude et/ou de phase, applicables à chacune des composantes spectrales sinusoïdales En(i) de la fonction 8(n), les fonctions de correction étant définies comme des sommes de composantes spectrales sinusoïdales An(i) et (pn(i), - une étape de détermination de la forme de la fonction A(n) à appliquer à chacune des composantes spectrales sinusoïdales En(i) de la fonction 8(n) en fonction de l'angle de pointage 0o de l'antenne. Selon l'invention la fonction la fonction A(n) a pour expression: A(n) = Ey(i) . An (i) + (1û y(i))' dn (i) y(i) où y(i) représente une fonction de choix valant 1 ou O. Selon l'invention encore, a valeur de y(i) est définie en fonction d'un seuil calculé à partir de l'angle de pointage 00. Selon une variante de mise en oeuvre le procédé selon l'invention comporte 35 une étape complémentaire de pondération des composantes spectrales des30 fonctions Aa(n) et A9(n), la fonction de correction globale A(n) ayant alors pour expression : A(n) = E f(i) y(i) • An (i) + f(i) • (1 - y(i)) (D, (i) y(i) Selon une variante de mise en oeuvre la pondération f(i) suit une loi linéaire croissante en fonction de i. Selon une autre variante de mise en oeuvre la pondération f(i) suit une loi 10 logarithmique croissante en fonction de i. L'invention a également pour objet un procédé de formation du diagramme de rayonnement F(0) d'une antenne réseau caractérisé en ce qu'il associe le procédé de correction selon l'invention à une étape de détermination de la loi 15 de phase de pointage, permettant de réaliser le diagramme de rayonnement, à partir des positions réelles des éléments rayonnants. DESCRIPTION DES FIGURES 20 Le principe de fonctionnement de l'invention, ainsi que ses avantages seront mieux compris à la lecture de la description qui suit, description faite au travers d'un exemple d'application non limitatif, et qui se réfère aux figures annexées qui présentent: 25 - la figure 1, un schéma représentant une antenne prise comme exemple pour la description du procédé selon l'invention, - la figure 2, la courbe représentative du diagramme de rayonnement théorique idéal de l'antenne de la figure 1, - la figure 3, la courbe représentative du diagramme de rayonnement réel d'une antenne du type de la figure 1 présentant une erreur sinusoïdale de positionnement des éléments rayonnants, 30 - la figure 4, la courbe représentative du diagramme de rayonnement réel d'une antenne du type de la figure 1 présentant l'erreur de positionnement de la figure 3, corrigé par modulation de phase des composantes spectrales de la fonction d'erreur de positionnement associée, - la figure 5, la courbe représentative du diagramme de rayonnement réel d'une antenne du type de la figure 1 présentant l'erreur de positionnement de la figure 3, corrigé par modulation d'amplitude des composantes spectrales de la fonction d'erreur de positionnement associée, 10 - les figures 6 et 7, les courbes représentatives des diagrammes de rayonnement réel d'antennes du type de la figure 1 présentant des erreurs de positionnement quelconques, 15 - la figure 8, l'organigramme de principe du procédé de correction selon l'invention, - la figure 9, les représentations graphiques d'exemples de lois d'évolution du seuil de sélection du type de correction, en fonction du 20 pointage 00, - la figure 10, les représentations graphiques de deux lois de pondération des termes de correction appliqués aux composantes spectrales de la fonction d'erreur de positionnement. DESCRIPTION DETAILLEE L'exemple utilisé pour décrire le procédé selon l'invention dans son principe, est celui d'une antenne telle que celle illustrée par la figure 1. Cette 30 antenne 11 est constituée d'un ensemble de N poutres rayonnantes 12, chaque poutre comportant un ensemble de sources 13 alignées le long de l'axe de la poutre. Les N poutres constituant l'antenne sont assemblées mécaniquement de façon à ce que les axes des poutres soient espacés les uns des autres d'une 35 distance P constante, ou pas, de sorte que, les poutres étant assemblées, 25 l'antenne ainsi constituée se présente comme un ensemble de sources, disposées en quinconce et équidistantes les unes des autres. Un tel arrangement de sources permet, si le positionnement des sources les unes par rapport aux autres est suffisamment précis, d'obtenir un diagramme de rayonnement F(0) semblable à celui présenté à la figure 2. Ce diagramme de rayonnement présente, pour une direction 00 pointée par l'antenne, un lobe principal 21 orienté dans la direction 0o et centré autour de 00, et des lobes secondaires 22, 23, orientés dans des directions autres que la direction 00, dont le niveau est significativement plus faible que celui du lobe principal, de l'ordre d'au moins quarante décibels au minimum. Idéalement, pour une antenne constituée de sources identiques, espacées de manière rigoureusement constante, on obtient un diagramme de rayonnement Fidéal pour lequel le niveau des lobes secondaires décroît sensiblement, de manière continue, à mesure que la direction 0 considérée s'éloigne de la direction 00 pointée par le lobe principal. On obtient ainsi par exemple pour des lobes 23 centrés sur une direction 0 présentant un angle supérieur à 60 par rapport à la direction 00, un niveau de lobe présentant un affaiblissement supérieur à 50 dB par rapport au niveau du lobe principal 21. Dans la pratique, réaliser de manière simple une telle antenne idéale, s'avère difficile. En effet un positionnement précis des différentes sources suppose d'une part que sur chaque poutre les sources se trouvent correctement disposées et alignées, et que d'autre part les poutres soient assemblées les unes aux autres avec un espacement P rigoureusement constant. Or si l'alignement des sources sur une poutre, alignement qui résulte de la fabrication de la poutre rayonnante considérée, est assez facilement maîtrisable, l'assemblage des poutres entre elles avec un pas d'écartement sensiblement constant est généralement plus difficile à réaliser, du fait en particulier de la précision de positionnement. Une antenne assemblée présente donc généralement des défauts de positionnement des sources qui se traduisent par une altération du diagramme de rayonnement. Cette altération se traduit elle-même principalement par la présence sur le diagramme F(0) de lobes secondaires de niveau relativement élevé. Ces lobes secondaires qui apparaissent pour certaines valeurs particulières de 0, ont un niveau sensiblement plus élevé que ceux des lobes secondaires présentés, pour les mêmes valeurs de 0, par le diagramme Fidéal. On obtient alors un diagramme de rayonnement Fperturbé(0). La présence de lobes secondaires de fort niveau a pour conséquence fâcheuse d'altérer la directivité de l'antenne ainsi réalisée, cette altération de la directivité pouvant se traduire en particulier, pour un radar équipé d'une telle antenne, par une augmentation du niveau de fausse alarme. Le procédé selon l'invention propose une solution pour corriger cette altération en appliquant au diagramme Fperturbé(0) un terme correctif variable AF. De la sorte il est avantageusement possible, selon l'invention, d'obtenir un diagramme d'éclairement Fcorrigé(0) ayant pour expression générale: Fcorrigé (e) = Fperturbé (e) + AF(e) [1] Selon l'invention, le terme AF(0) est déterminé de façon à ce que l'on 15 obtienne après application du procédé, un diagramme Fcorrigé(0) aussi proche que possible du diagramme Fidéal(0). Dans le cas d'une antenne telle que celle de la figure 1, prise comme exemple non limitatif, le diagramme Fidéal (0) peut être défini par la relation 20 suivante: N Fidéal (e) an .ej.k.n.P.(sin 9ùsin 90 ) [2] n=1 dans laquelle 0o représente la direction pointée par l'antenne (direction du 25 lobe principal de F(0)) et P le pas d'implantation des poutres. Le facteur k est la constante de propagation de l'onde considérée. On rappelle ici que k est égale à 2n/2 , étant la longueur d'onde du signal auquel est adapté l'antenne. 30 Fidéal est ainsi défini pour chaque direction 0 comme la somme des contributions, des diagrammes élémentaires relatifs à chacune des poutres rayonnantes composant l'antenne. Il se présente donc comme une somme de termes complexes de module ai et d'argument k.i.P.(sin 6 ù sin 80) . Si on considère une antenne réelle, pour laquelle le positionnement des poutres rayonnantes présente des imperfections, le diagramme Fperturbé(0) correspondant peut être défini de manière générale par la relation suivante: Nx Fperturb (8) = Eau .ej.k.n.P.sinep e(j.k.(n.P+Sn)sine) n=1 dans laquelle 8n correspond pour la poutre de rang n à la valeur de l'erreur de position. 10 Par suite, pour obtenir un diagramme d'éclairement corrigé proche du diagramme idéal, il est possible d'appliquer à chaque terme de la somme de la relation [3] un facteur de correction A(n) . La suite de la description présente la façon dont, selon l'invention, est déterminé ce facteur de correction. 15 On s'intéresse à présent aux figures 3 à 5. La figure 3 présente le diagramme de rayonnement obtenu avec une antenne dont les poutres présentent un défaut de positionnement particulier 20 dont la valeur est une fonction sinusoïdale de la position de la poutre dans l'antenne. Ce type de défaut, par ailleurs théorique, permet avantageusement de comprendre le principe de correction proposé par l'invention. Pour ce défaut particulier, l'erreur de positionnement de chacune des poutres 25 peut être définie par la relation suivante: VnE[1,Nx], 8(n) = 8n = 8P . sin(ncp +(po) [4] Les erreurs 8(n) constituent les composantes d'un vecteur qui sera 30 appelé vecteur d'imperfections. Ce vecteur est chargé des erreurs de positionnement des poutres aux indices correspondants. L'amplitude 8P de l'erreur de position est ici considérée comme faible, de l'ordre de quelques pour cents de la valeur exacte du pas P, ce qui est 35 généralement le cas dans la réalité. [3] Chaque erreur est par ailleurs considérée ici comme une erreur corrélée sur l'ensemble d'une poutre, le positionnement des sources sur la poutre étant supposé suffisamment précis pour ne pas contribuer significativement à la dégradation du diagramme d'éclairement. Il est ainsi possible d'étudier les perturbations occasionnées par ce défaut en raisonnant non pas sur le réseau de sources réel, à deux dimensions, mais sur un réseau linéaire de sources équivalent. Comme on peut le constater en comparant la figure 3 à la figure 1, une telle perturbation des positions des éléments rayonnants a pour effet d'élever le niveau de deux lobes secondaires particuliers 31 et 32, centrés sur des directions 00+01 et 00-01, symétriques l'un de l'autre par rapport au lobe principal centré sur 00. La valeur de 01 est par ailleurs liée au terme (pi du défaut de positionnement sinusoïdal 8(n) par la relation suivante: u1 = cp = k.P. (sin(01) ù sin(00 )) [5] Le défaut de positionnement considéré ici génère ainsi une altération localisée qui ne modifie en rien l'allure du diagramme d'éclairement pour le reste du domaine d'observation. L'altération localisée en 0 = 01 peut avantageusement être corrigée en appliquant la correction sélective appropriée. Le diagramme de rayonnement F(0) pouvant être représenté, selon la relation [3], comme une somme de termes complexes, il est possible d'appliquer une correction sélective portant uniquement sur le module ou bien sur l'argument, la phase, de chacun des termes. Par suite la valeur du terme de correction AF(0) vérifie alors la condition de compensation suivante relative à la direction 01: Fcorrigé (01) = Fperturbé (01 ) + AF(01) = Fideal (01) [6] On considère dans un premier temps le cas où la correction appliquée 35 aux signaux reçus ou émis par chacune des poutres est une modulation de phase sinusoïdale A (n) définie pour une poutre de rang n quelconque (n variant de 1 à N) par la relation suivante: A(p(n) = c. sin(n.u1 ù [3) = c • sin[n • k • P • (sin 01 ù sin 0o) ù [i] [7] dans laquelle 00 est la direction pointée par l'antenne. La variable u1 est quant à elle définie par la relation suivante: u1 =k•P•(sinO1 ùsin0o) [8] Le diagramme de rayonnement Fperrurbé (0) correspondant à l'antenne réelle a pour expression: N ùj.k.n.P.sinOo (j.k.(n.P+Sn)sinO) Fperturbé(0) _ Ean.e e n=1 où l'expression de la perturbation 8n est donnée par la relation [4]. D'autre part, le diagramme de rayonnement Fcorrigé (0) correspondant à l'antenne réelle dont on souhaite corriger les aberrations locales des lobes 20 secondaires a pour expression: N ( Fcorrigé 0) _ an .eùj.k.n.P.sin0oe(j.k.(n.P+Sn )sin 9) ej0(p(n) [10] n=1 où l'expression de la correction A p(n) est donnée par la relation [7] 25 Si l'on considère d'une part que la correction A(p(n) apportée doit rester faible devant F(0) et d'autre part que, comme il a été dit précédemment, le terme 8n est petit devant le pas P, il est possible de simplifier la relation [9] (développement limité à l'ordre 1 de ei `D(n)) pour obtenir les relations 30 suivantes: 10 15 [9] F N j.k.n.P.sin e 8p (j.k.(n.P+Sn )sin 6) Fcorrigé (0) = ~an.e n=1 N + a 0 (n) • g-1.k.n.P.sin00e(jk.(n.P+Sn ) sin 0) j.En n=1 Fcorrigé (0) = Fperturbé (0) N + j E an • A (n) (1 + j k • 8n • sin 0) . el-k•n•P•(sin9-sinO0) n=1 Fcorrigé (0) = Fperturbé (0) N N + j•Ean A (n)•elnu + j.~an n=1 n=1 où u est défini par la relation [5] N Dans la relation [12], le dernier terme j • Ean • 0(p(n) • j. k • 8n •sin 0 • ej n" n=1 10 apparaît, en tant que produit de deux facteurs (8n et A (n)) de faible amplitude, comme négligeable devant les deux autres termes. Par suite le facteur de réseau Fcorrigé(0), ou diagramme de rayonnement corrigé, peut être défini en fonction de la variable u par la relation suivante: N Fcorrigé (U) ùFperturbé (u) + J.~ an . (n) . el.n.0 n=1 ou, en remplaçant A 9(n) par son expression littérale, par les relations : N Fcorrigé (U) ù Fperturbé (U) + j.~ an • c • sin(n • u1 ù [3) . el"' n=1 N el(n W-(~) ù e-1(n u1 a) Fcorrigé (U) ùFperturbé (u) + j. an ' c - el.n.0 [16] n=1 2j Par suite, AF(u) est alors défini par la relation suivante: [12] oq,(n) k 8n • sin 0 • ej' [13] 15 20 [14] [15] 15 AF(u) = 2 . [e ia . Fidéal (u + u1) ù eiP Fidéal (u ù u1)] [17] dans le cas particulier où u = u1: AF(u1) = Fcorrigé (u1) ù Fperturbé (u1) = 2 • [e Jr' . Fidéal (2u1) ù ee • Fidéal (0)] [18] ou encore, en considérant Fidéai(2u1) comme négligeable: AF(u1) _ ù 2 eia . Fidéal (0) [19] où Fidéal(0) représente la valeur de F(u) pour u correspondant à la direction pointée 00 Par suite, en utilisant la condition sur u1 de la relation [6] on obtient: AF(u1) = Fidéal (u1) ù Fperturbé (u1)=ù 2 e • Fidéal (0) et donc : 20 c =2• Fperturbé (u1) ù Fidéal (u1) Fidéal (0) [20] Fperturbé (u1) ù Fidéal (u1) Fidéal (0) ~ [3=arg [21] Ainsi comme le montrent les relations [20] et [21], il est possible de 25 définir une fonction de correction de phase permettant de corriger les altérations du diagramme lorsque celles-ci résultent d'une erreur sinusoïdale de positionnement des poutres qui composent l'antenne considérée, cette erreur de positionnement se traduisant, comme l'illustre la figure 3, par une remontée de deux lobes secondaires pour les directions 00+01 et 00-01. 30 Cependant, il faut noter que, comme on peut le constater sur la figure 4, lorsque la correction appliquée est une modulation de phase, cette correction a pour effet de supprimer le lobe secondaire parasite 31 initialement présent dans la direction 00+01, mais en même temps de renforcer le niveau du lobe secondaire 32 initialement présent dans la direction symétrique 00-01. Ce deuxième effet, indésirable, conduit à considérer un deuxième type de correction consistant à appliquer une correction d'amplitude. On considère à présent le cas où la correction appliquée aux signaux reçus ou émis par chacune des poutres, est une modulation d'amplitude sinusoïdale Aa(n) définie pour une poutre de rang n quelconque (n variant de 1 à N) par la relation suivante: Da(n)=c'.sin(n.u1 ù[3')=c •sin[n k•P•(sin01 ùsin00)ù[3'] [22] En appliquant un raisonnement analogue à celui utilisé dans le cas d'une correction de phase, on peut écrire: N Fcorrigé (0) ù E [an + Aa(n)]. eùj•k.n.P.sin 00e(j.k.(n.P+Sn)sin 8) [23] n=1 Dès lors, en suivant un raisonnement analogue à celui mené précédemment, on établit la relation suivante: AF(u1) = Fidéal (u1) ù Fperturbé (u1)=ù j ' 2 . e'a Fidéal (0) [24] et on obtient : Fperturbé (u1) ù Fidéal (u1) Fidéal (0) c =2- [25] / Fperturbé (u1) ù Fidéal (U1) Fidéal (0) 1 /3 = arg [26]30 Ainsi les relations [25] et [26], permettent de définir une fonction de correction d'amplitude capable de corriger les altérations du diagramme lorsque celles-ci résultent d'une erreur sinusoïdale de positionnement des poutres. La correction apportée ici présente en outre l'avantage, comme l'illustre la figure 5, de supprimer la remontée de lobe secondaire parasite 31 initialement constatée présente dans la direction 00+01, mais également la remontée de lobe parasite 32 initialement constatée dans la direction symétrique 00-01. Ainsi, comme on peut le constater au travers des paragraphes précédents, il est possible, dans le cas d'une antenne constituée de poutres rayonnantes et présentant un défaut sinusoïdal dans le positionnement des poutres, de déterminer une fonction de correction permettant de corriger les effets de ce défaut de positionnement. La correction appliquée au diagramme de rayonnement F(0) de l'antenne est une correction sinusoïdale de phase ou d'amplitude. Des simulations et des essais menés par la déposante montrent par ailleurs que, comme cela a déjà été évoqué, la correction par modulation d'amplitude du diagramme F(0) apparaît, en particulier lorsque la différence 01-00 est faible, comme donnant les meilleurs résultats. Cette correction par modulation d'amplitude permet en effet, contrairement à la correction par modulation de phase, de supprimer à la fois le lobe parasite présent dans la direction 00+01, mais aussi le lobe parasite image présent dans la direction Oo-0,. Cependant il faut noter que, l'application d'une correction consistant à moduler la phase de la loi d'éclairement n'est pas à exclure, car, des simulations faites pour des dépointages 00 croissants, permettent de constater qu'une fois que le lobe image est défini pour une direction théorique telle que la différence 00-01 est supérieure à 90 , sa prise en compte n'est plus nécessaire, alors que parallèlement la correction par modulation de phase donne des résultats meilleurs que la correction par modulation d'amplitude. Ainsi, une correction optimale consiste à utiliser l'une ou l'autre correction en fonction du dépointage et de l'angle pour lequel apparaît un lobe secondaire parasite. On s'intéresse ensuite aux figures 6 et 7 qui permettent d'exposer la manière de généraliser le principe de correction, décrit dans les paragraphes précédents, au cas d'un défaut de positionnement quelconque. La figure 6 illustre le cas particulier où le défaut de positionnement des poutres peut être décrit comme une fonction somme de deux fonctions d'erreur sinusoïdales. Comme le montre la figure ce type de défaut occasionne la remontée de deux lobes secondaires 61 et 62 pour des directions 00+01 et 00+02 sur le diagramme de rayonnement F(0), et des deux lobes symétriques 63 et 64 pour les directions 00-01 et 00-02. Le cas illustré par la figure 6 reste d'interprétation assez simple, et on conçoit aisément que la méthode de correction décrite précédemment puisse s'appliquer de manière relativement directe à ce cas particulier. L'erreur de positionnement 8(n) de chacune des poutres peut alors être définie comme une combinaison linéaire de deux erreurs sinusoïdales, et peut donc être corrigée au moyen d'une fonction Èa(n) (ou A (n)) comportant un terme correcteur relatif à la première erreur et un terme relatif à la deuxième erreur. La figure 7 illustre quant à elle le cas général où le défaut de positionnement ne peut être décrit de manière simple. On assiste dans ce cas à une déformation générale du diagramme de rayonnement, avec modification du niveau des lobes secondaires et diffus. Le niveau de certains lobes devient ainsi plus faible que dans le cas idéal, alors que le niveau de certains autres devient franchement plus élevé. Un tel diagramme ne permet pas de déterminer de manière simple l'expression de la fonction caractérisant l'erreur de positionnement et par suite, ne permet pas l'élaboration directe de la correction permettant de retrouver le diagramme F;déai(0) qui caractérise l'antenne idéale, réalisée sans erreur de positionnement. Le procédé selon l'invention propose une solution permettant de déterminer cette correction notamment dans le cas général où l'erreur de positionnement des poutres est une erreur aléatoire. La suite de la description présente les différentes étapes du procédé, résumées sur la figure 8.35 Comme l'illustre la figure 8, le procédé selon l'invention comporte principalement 5 étapes. La première étape 81 est une étape de mesure qui consiste, après assemblage des différents éléments rayonnants de l'antenne, à déterminer pour chaque élément l'écart existant entre la position théorique de l'élément et sa position réelle. Dans le cas de l'antenne constituée de poutres, prise comme exemple non limitatif dans la description, cette première étape consiste à mesurer l'erreur de positionnement de chaque poutre par rapport à sa position théorique. Ces écarts constituent les composantes d'un vecteur d'écart I mesuré, ou vecteur d'imperfections. Dans l'exemple d'une antenne constituée de N poutres rayonnantes le vecteur I comporte donc N composantes. Comme on l'a vu précédemment, dans le cas général correspondant à la figure 7, l'erreur de positionnement 8(n) ne peut être décrite comme une simple fonction sinusoïdale du rang n de la poutre. Par suite, cette erreur se traduit non pas par la remontée du niveau d'un lobe particulier dans le diagramme F(0), mais par une altération répartie sur tout le diagramme. La seconde étape 82 et la troisième étape 83 du procédé selon l'invention ont pour objet de modéliser la fonction d'erreur de positionnement 8(n), connue par les mesures effectuées qui composent le vecteur I, par une fonction définie comme une somme d'erreurs sinusoïdales de périodes P, 2P...(N-1)P. La seconde étape 82 procède à la décomposition spectrale de la fonction d'erreur. Cette décomposition est par exemple réalisée par l'application d'une transformée de Fourier discrète à la suite d'échantillons constituée des composantes du vecteur d'imperfections I. Cette transformation produit un nouveau vecteur dd de même longueur que le vecteur d'imperfections, et présentant une symétrie du fait du caractère réel des éléments composant le vecteur I traité. Les composantes dd; du vecteur dd sont définies par la relation suivante: N (-j.2.11.0 1).(n-1 ddi - -E8(n).e N n=1 [27] avec ie[1,N], Les termes dd; ainsi obtenus représentent les composantes spectrales complexes de la fonction d'erreur dans le domaine spatial. Chaque composante spectrale représente elle-même une erreur sinusoïdale de période i, cette erreur pouvant être considérée comme responsable de la remontée d'un lobe secondaire particulier du diagramme F(0). La troisième étape 83 consiste quant à elle à effectuer la transformation inverse de celle effectuée lors de l'étape précédente. En effet, la transformée de Fourier inverse du vecteur dd redonne accès au vecteur d'imperfections I qui est alors écrit sous la forme d'une somme finie de sinusoïdes. On peut ainsi écrire: 8(n) = N EIdd;I. cos(211. (i -1).(n -1) + arg(dd1)) [28] X 1-Î avec nE[1,N] L'erreur 8(n) se présente ainsi comme une somme de perturbations sinusoïdales pouvant être compensées au moyen d'une fonction de correction se présentant comme une somme de modulations analogues aux fonctions Aa(n) et A (n) décrites précédemment au moyen des relations [7] et [22] notamment. Cette somme de modulations compensera un par un les termes spectraux du vecteur d'imperfections. La quatrième étape 84 consiste à déterminer l'expression des corrections théoriques d'amplitude Aa(n) et de phase A (n), permettant de corriger l'erreur 8(n) définie par la relation [28]. Dans le cas d'une antenne telle que celle prise comme exemple, présentant un diagramme de rayonnement, ou facteur de réseau, F(0), perturbé du fait d'erreurs de positionnement des poutres rayonnantes et pour un dépointage nul, la correction à appliquer à chaque composante sinusoïdale de l'erreur sera, comme on l'a vu précédemment, une correction par modulation de 16 l'amplitude plutôt qu'une correction par modulation de phase. Par suite,le diagramme de rayonnement corrigé Fcorrigé(0), pourra être défini par la relation suivante: N F(0) _ E an.(1 + 2II .A a (n)). e().k.(n.P+Sn )sin 9) n=1 où Aa(n) représente la fonction de correction globale. Si l'on considère en outre que les erreurs de positionnement 8n sont petites devant la distance P (pas) séparant deux éléments rayonnants du réseau, on 10 peut alors appliquer un développement limité au terme e~ksn sine . On obtient alors la relation suivante: N F(6) = 1 an.(1 + 2H .L a (n)).( 1 + j.k.5n.sin O).e().k.n.P.sin 0) [30] n=1 15 Par suite, en ne gardant que les termes du premier ordre, le terme de modulation étant aussi de petite amplitude, et en identifiant l'expression obtenue à celle du diagramme Fidéai(0), on obtient la relation suivante : sine . Aa(n) _ ù .j.Sn avec n e [1,N]. On constate donc sur cet exemple simple, par ailleurs généralisable, que Aa(n) est une fonction qui présente l'inconvénient de dépendre à la fois des 25 erreurs de positionnement des différents éléments de l'antenne, et de la direction d'observation 0 considérée, direction qui n'est pas un paramètre intrinsèque de l'antenne. De manière générale, cette double dépendance, rend complexe, voire impossible, la détermination d'une correction Aa(n) efficace au moyen de la relation [31]. 30 Pour pallier cette difficulté tout en utilisant le principe de décomposition de la fonction d'erreur 8(n) en une somme d'erreurs de positionnement sinusoïdales, le procédé selon l'invention propose de substituer aux termes [29] [31] 20 dépendants de l'angle 0, un terme majorant. Ainsi, dans l'exemple précédent, on substitue à la relation [31] la relation suivante: Aa (n) = j.S~ [32] La relation [32] correspond à une substitution, dans la relation [31] de l'expression sine par la valeur maximale qu'elle peut prendre sur l'espace d'observation couvert par le diagramme F(0) (de -90 à +90 ), à savoir la valeur 1. De la sorte, l'expression de Aa(n) ne dépend plus que du seul paramètre n, qui représente l'élément d'antenne considéré. Cette maximalisation avantageuse qui permet de s'affranchir du paramètre sine, n'entraîne en outre aucune conséquence sensible sur la correction effectuée, dans la mesure où cette correction est de toute façon de faible valeur. Par suite, selon l'invention, les composantes du vecteur d'imperfections I étant remplacées par celles du vecteur obtenu par la transformée inverse de la transformée de Fourier discrète de I, on peut écrire la relation suivante: N Aa (n) _ • N Elddil sin 2.IL (n -1). 0 -1) + arg(dd1) [33] i-1 C N avec n e [1, N]. Cette relation est déduite de la relation [28] dans laquelle les termes 25 en "cosinus" deviennent des termes en "sinus" pour prendre en compte le facteur j de la relation [32]. Ainsi, selon l'invention, la détermination de la fonction de correction par modulation d'amplitude Aa(n) est réalisée de manière simple et 30 indépendante de la direction d'observation 0, à partir des composantes du vecteur obtenues par application d'une transformation de Fourier inverse dd. Cette détermination est encore facilitée par le caractère symétrique de la décomposition spectrale de l'erreur de positionnement, qui permet de ramener l'expression de Aa(n) à une somme de N' termes, N' étant égal à la partie entière du nombre (N/2+1) (i.e. N' = floor(N/2+1). On peut alors déterminer Aa(n) à l'aide de la relation suivante: 0a (n) _ N . Idd;I sinC 2.rI.(n N1).(i -1) + arg(dd;JJJ ) I [34] i1 avec n e [1, N]. Finalement, pour un arrangement dans lequel, l'antenne étant constituée de poutres rayonnantes, l'espacement P entre poutres est 10 sensiblement égal à la demi-longueur d'onde (X/2) de fonctionnement, l'expression de Aa(n) est donnée par la relation suivante: Aa (n) = P -1 . EIdd;I • sin( 2.IL(nN 1).(i -1) + arg(ddi) [35] i-1 x 15 avec n [1, N]. Cette dernière relation approchée permet d'obtenir une expression de Aa(n) indépendante de la fréquence. Le domaine fréquentiel de fonctionnement de l'antenne étant généralement supérieur ou égal à 10% de la valeur de la fréquence centrale, la constante P, représentant la distance 20 entre deux éléments rayonnants, est égale à peu près à la moitié de la longueur d'onde pour la fréquence la plus haute du domaine considéré. Elle constitue donc une bonne estimation de pour approximer dans l'expression de la modulation le terme dépendant de la fréquence, et rendre ainsi le calcul de Aa(n) indépendant de la fréquence. 25 Le procédé selon l'invention permet de la même façon de déterminer une expression de la fonction de correction par modulation de la phase, indépendante à la fois de la direction d'observation 0 et de la longueur d'onde 1. Le diagramme Fcorr;9é(0) a alors pour expression: N F(0) _ E an •e(i.k.(n.P+Sn)sin 6) ejA(p(n) [36] n=1 30 Par le même cheminement que précédemment, en opérant un développement limité sur l'exponentielle e' 9(n) contenant la modulation de phase (supposée de petite amplitude), on trouve l'expression de la modulation de phase, dont l'expression, très proche de celle de la modulation d'amplitude, est donnée par la relation suivante : 0 n = 1 N dd •cos( 2.I7.(nù1).(iù1)+ar dd1)J [37] i1 avec n e [1, N]. La cinquième étape 85, consiste à effectuer la détermination automatique du type de corrections à appliquer aux différentes composantes de F(0) en fonction du dépointage de l'antenne. Cette étape est nécessaire pour appliquer celle des corrections qui est la plus appropriée. En effet, comme on l'a vu précédemment, la correction par modulation d'amplitude est celle qui donne les meilleurs résultats lorsque l'angle de dépointage du diagramme de rayonnement est faible. Dans ce cas, la correction par modulation de phase a pour effet, comme l'illustre la figure 4 de relever le niveau du lobe secondaire image 42. En revanche, lorsque le dépointage est important, le lobe image disparaît du domaine visible, et la correction par modulation de phase du lobe parasite s'avère plus performante. II est donc intéressant d'appliquer de manière sélective sur une composante spectrale de rang i donné, telle correction plutôt que telle autre. Par suite, le problème se pose de déterminer de manière automatique, en fonction du dépointage considéré, quelle type de correction appliquer à une composante donnée. Le traitement effectué à cette étape du procédé permet avantageusement de résoudre ce problème. Le principe du traitement effectué lors de cette cinquième étape consiste à introduire dans l'expression des fonctions de correction Aa(n) et 4p(n) une fonction binaire y(i) prenant la valeur 0 ou 1 en fonction du dépointage. Les corrections ont alors pour expression: Aa(n)ù P N Ey(i)•Idd;l•sin 2.IL(n N1).(iù1)+arg(dd1) [38] i-1 (n)= PN''(1ùy(i)) ldd;•cos[ 2n.(n N1).(i -1)+arg(ddi )J [39] Selon l'invention, la valeur de la fonction y(i) peut être définie par rapport à un seuil de référence ko entier, dont la valeur varie de façon décroissante en fonction du dépointage considéré. y(i) est également définie de façon à ce que, au fur et à mesure que le dépointage augmente pour les composantes spectrales d'ordre i élevé, composantes qui ont pour effet de faire apparaître des lobes parasites dans des directions 01 éloignées de la direction de pointage 00, la correction par modulation de phase soit préférée, du fait que le lobe parasite image n'est dans ce cas plus présent dans le domaine visible. Ainsi, en choisissant judicieusement le seuil ko, la fonction y(i) peut par exemple être définie de manière simple, pour tout i compris entre 1 et N', par l'expression suivante: " Si i k0 alors y(i) = 1, sinon y(i) = 0." Selon l'invention, le seuil ko varie théoriquement de 0 au nombre N' de termes de modulation composant la fonction de correction considérée. En pratique, on exclut les valeurs proches des extrêmes en définissant un seuil bas, et éventuellement un seuil haut. En effet, une valeur extrême conduit à un type de correction unique, ce qui enlève tout intérêt à cette cinquième étape. La figure 9 illustre, à titre d'exemple, l'évolution du seuil k0 en fonction de la 25 valeur du dépointage 00, dans le cas particulier où k0 est défini de manière analytique par la relation suivante: k0 = floor(N~P.(1ùlsin0ol)+0.5) [40] 30 où "floor" représente la fonction partie entière. Cette relation fait dépendre la valeur de k0 de l'angle de dépointage 00, et du rapport de la taille de l'antenne (donnée par le terme N.P) à la longueur d'onde de travail a,. 21 et Sur la figure 9, les courbes pointillées 91, 92 et 93 représentent l'évolution de ko, calculé à partir de la relation [40], pour trois valeurs distinctes de X, tandis que la courbe en trait plein 94 représente une loi d'évolution bâtie de manière empirique à partir des courbes précédentes, et ajustée en tenant compte de résultats de simulation, résultats qui montrent en particulier que la variation de X. dans une plage raisonnable n'a pas une influence significative sur l'efficacité de la correction apportée. La définition de ko à partir de la formule analytique [40] n'est pas optimale, en ce sens qu'elle rend ko dépendant de la fréquence, ce qui a pour conséquence d'altérer notablement la solution aux bornes de l'intervalle de fréquence considéré ; par contre, elle donne une bonne idée de l'allure générale que la fonction doit embrasser et permet un affinage empirique, au moyen de résultats de simulation par exemple. Le type de correction à appliquer aux différentes composantes de F(0) étant ainsi déterminé, la fonction de correction d'erreur appliquée peut être décrite, selon l'invention, par la relation suivante: A(n) = 0a (n) + 09(n) [41] A l'issue de ces cinq étapes, le procédé selon l'invention permet donc d'appliquer au diagramme de rayonnement réel correspondant à une antenne dont certains éléments rayonnants présentent des erreurs de positionnement, une correction qui se présente avantageusement comme une somme de corrections élémentaires, chaque correction ayant pour objet de corriger les effets d'une composante particulière de la fonction d'erreur de positionnement 8(n). Les résultats de simulations réalisées par la déposante montrent en outre que la correction globale ainsi déterminée montre son efficacité, en particulier dans un domaine de pointage 00 limité, s'étendant typiquement sur 25 autour de l'axe de l'antenne correspondant à un pointage 00 nul. Cependant, pour les pointages relativement éloignés, cette correction s'accompagne d'une certaine altération du niveau des lobes secondaires proches du lobe principal, altération qui peut par exemple se traduire par une élévation, de l'ordre de 3 dB, des trois premiers lobes secondaires de F(0). Pour améliorer l'efficacité de la correction appliquée, le procédé selon l'invention peut être complété, dans le cas où cela s'avère nécessaire, par une étape intermédiaire 86, qui prend place entre la quatrième étape 84 et la cinquième étape 85. L'étape consiste à limiter la dégradation constatée sur les premiers lobes secondaires, notamment dans le cas de pointages élevés. Cette étape intermédiaire 86 consiste principalement à appliquer aux fonctions de correction Aa(n) et o(,(n) déterminées à l'issue de l'étape 84, une loi de pondération f(i) des composantes spectrales, dont le rôle consiste à appliquer une correction faible aux composantes d'indices i faibles et d'appliquer aux autres composantes une correction dont la valeur croît avec l'indice. Ainsi, les lobes secondaires proches, correspondant à des erreurs de fréquence spatiale faible, seront peu affectés par la correction appliquée, alors que les lobes secondaires plus lointains seront effectivement corrigés. Pour réaliser cette correction, il est possible d'utiliser diverses lois. La figure 10 donne une représentation graphique de deux lois f1(i) et f2(i) prises comme exemples non limitatifs. Ces lois, représentées par les courbes 101 et 102, sont respectivement definies, pour i e [1, N'], par les relations suivantes: f1(i) iNù 1 [42] et 25 f2 (i) = In(1 + 1,72 • (i -1)/N' [43] A l'issue de cette étape, on obtient ainsi des fonctions de correction dont les expressions sont les suivantes: A' (n)= P N . E f(i) • Y(i) • Idd. J . sin 2.II.(n N1).(i -1) + arg(ddi)J [44] i-1 J 30 et n = 1 • N f(i) • 1ù i ddiI•cos 2'II'(nù1).0 -1)+ar dd1)J [45] p .N N L'étape 85 de choix de la correction est alors mise en oeuvre sur ces fonctions pondérées. L'efficacité du procédé selon l'invention, tel qu'il vient d'être décrit, peut par ailleurs avantageusement être accrue en prenant en compte le vecteur des imperfections de positionnement I, pour établir la loi de phase à appliquer aux différents éléments rayonnants pour former le diagramme F(0) souhaité. En effet, si le calculateur responsable du calcul de la loi de phase à appliquer pour provoquer le pointage de l'antenne dans une direction 00, n'a pas connaissance des imperfections de positionnement des éléments rayonnants, le calcul de cette loi de phase est effectué avec des valeurs de position erronées. Dès lors la fonction de correction 0(n) déterminée pour le pointage 00 n'est plus adaptée à la correction du diagramme de rayonnement réel. La correction appliquée perd ainsi de son efficacité, en particulier pour les pointages dont la direction est éloignée de l'axe de l'antenne. Des simulations réalisées par la déposante montrent qu'en prenant en compte les erreurs de positionnement des éléments rayonnants dans le calcul de la loi de phase appliquée à ces éléments, on augmente sensiblement le domaine de pointage, typiquement jusqu'à 60 autour de 00, pour lequel la correction par le procédé selon l'invention est efficace. Pour des raisons de clarté de la description, le procédé selon l'invention a été présenté au travers d'un exemple particulier d'antenne formée de poutres rayonnantes, et pour laquelle les erreurs de positionnement se réduisent à de simples erreurs monodimensionnelles d'espacement des poutres sur un axe. Cet exemple particulier, non limitatif de la portée de l'invention, permet d'exposer le cas simple dans lequel les erreurs de positionnement constituent un vecteur d'imperfections I dont les composantes sont des grandeurs scalaires. Ainsi, l'exposé des différentes étapes du procédé n'est pas obscurci par l'emploi de formules générales complexes et non nécessaire à l'exposé du principe de l'invention. Le procédé selon l'invention est bien évidemment applicable à d'autres configurations d'antennes pour lesquelles se pose le problème des défauts de positionnement des éléments rayonnant par rapport à des positions théoriques données. REVENDICATIONS 1. Procédé de compensation des erreurs de positionnement des éléments rayonnants (12) d'une antenne réseau, caractérisé en ce qu'il comporte au moins: - une étape de mesure de l'erreur de positionnement de chaque élément rayonnant n (12) et de constitution d'un vecteur d'imperfections I, dont les composantes sont les valeurs de la fonction 6(n) d'erreur de positionnement , - une étape de modélisation de la fonction d'erreur 6(n) sous forme d'une somme de composantes spectrales sinusoïdale En(i) calculées par application d'une transformée de Fourier au vecteur d'imperfections I, suivie de l'application d'une transformée de Fourier inverse au vecteur dd ainsi obtenu, 15 - une étape de calcul de fonctions de corrections Aa(n), A (n), par modulation d'amplitude et/ou de phase, applicables à chacune des composantes spectrales sinusoïdales Er,(i) de la fonction 6(n), les fonctions de correction étant définies comme des sommes de 20 composantes spectrales sinusoïdales An(i) et (Dn(i), - une étape de détermination de la forme de la fonction A(n) à appliquer à chacune des composantes spectrales sinusoïdales En(i) de la fonction 6(n) en fonction de l'angle de pointage 0o de l'antenne. 2. Procédé selon la revendication 1, dans lequel la fonction la fonction Mn) a pour expression: A(n) = E y(i) . An (i) + (1û y(i)) . (Dn (i) y(i) où y(i) représente une fonction de choix valant 1 ou O. 10 25 30 3. Procédé selon la revendication 2, dans lequel la valeur de y(i) est définie en fonction d'un seuil calculé à partir de l'angle de pointage 00. 4. Procédé selon l'une quelconque des revendications 1 à 3, comportant une étape complémentaire de pondération des composantes spectrales des fonctions Aa(n) et A (n), la fonction de correction globale A(n) ayant alors pour expression : A(n) = E f(i) y(i) • An (i) + f(i) (1ù y(i)) • en (i) y(i) 5. Procédé selon la revendication 4, dans lequel la pondération f(i) suit une loi linéaire croissante en fonction de i. 6. Procédé selon la revendication 4, dans lequel la pondération f(i) suit 15 une loi logarithmique croissante en fonction de i. 7. Procédé de formation du diagramme de rayonnement F(0) d'une antenne réseau caractérisé en ce qu'il associe le procédé de correction selon l'une quelconque des revendications précédentes, à 20 une étape de détermination de la loi de phase de pointage, permettant de réaliser le diagramme de rayonnement, à partir des positions réelles des éléments rayonnants.10