FR 2902935 A1 20071228 FR 0605667 A 20060623 PROCEDE DE FABRICATION D'UNE ANTENNE A DIAGRAMME DE RAYONNEMENT OPTIMISE SELON LES CONTRAINTES. L'invention concerne un procédé de fabrication d'une antenne à diagramme de rayonnement optimisé selon les contraintes. En particulier, l'invention s'applique à la fabrication d'antennes comportant des éléments rayonnants regroupés en réseaux. L'antenne obtenue par le procédé de fabrication selon l'invention a une configuration géométrique et une alimentation électrique auxquelles correspond un diagramme de rayonnement de l'antenne dont les lobes secondaires et les lobes réseaux sont au niveau le plus bas possible selon les contraintes, tout en conservant io une puissance maximale dans le lobe principale. L'invention peut s'appliquer à diverses antennes utilisées dans différents dommaines comme par exemple le domaine des radars ou des télécommunications. Les antennes actives, c'est-à-dire les antennes composées de 15 multiples sources émettrices et réceptrices, sont utilisées dans de nombreux domaines d'application, comme par exemple, le domaine des radars ou encore des télécommunications. Les sources émettrices et réceptrices, aussi appelé modules d'émission et de réception (ou selon l'expression anglo-saxonne Transmitter-Receiver module) ou encore éléments rayonnants, 20 permettent de réaliser un balayage électronique du faisceau d'antenne sans déplacement mécanique et sans perte de puissance. Les antennes actives présentent de nombreux avantages. Elles permettent de réaliser un balayage électronique sans perte de puissance, de maîtriser dynamiquement le diagramme d'antenne, de conforrner l'antenne à une surface de forme 25 donnée permettant une intégration de l'antenne avec le porteur... Dans un système comportant une antenne active formée de modules actifs, une voie de réception correspond à l'ensemble des moyens nécessaires pour traiter les données issues d'un module d'émission et de réception ou d'un groupe de modules d'émission et de réception. A partir des 30 différentes voies de réception comprise dans le système, il est alors possible de mettre en oeuvre des techniques de formations de faisceau par le calcul, permettant en outre d'observer instantanément un large secteur angulaire. Cependant, la mise en oeuvre industrielle de système comportant une antenne active disposant d'autant de voies de réception que de module d'émission et de réception sont coûteuses et complexes, particulièrement pour les antennes comportant plusieurs milliers de ces modules. Pour contourner ce problème, il est connu de partitionner l'antenne active en sous-réseaux de modules d'émission et de réception. Les sous-réseaux de modules d'émission et de réception regroupent un nombre déterminé de modules d'émission et de réception, par exemple de l'ordre de quelques dizaines de modules d'émission et de réception. Le nombre de chaîne de voies de réception est donc alors égal au nombre de sous-réseaux de modules d'émission et de réception. 1 o Une antenne active dite optimale, notamment d'une antenne active comportant des sous-réseaux, est définit par des caractéristiques permettant notamment de minimiser les lobes secondaires et les lobes réseaux du diagramme de rayonnement de l'antenne. Parmi les caractéristiques optimisés de l'antenne optimisale, on peut notamment citer 15 la position sur ladite antenne des éléments rayonnant (ou modules d'émission et de réception), la pondération de l'alimentation des éléments rayonnants ( ou modules d'émission et de réception ), ou encore la forme des regroupements en sous-réseaux des élements rayonnants (ou modules d'emissions et de réception ) de l'antenne active ainsi que leur alimentations. 20 Il existe des procédés de fabrication d'antennes dont les caractéristiques sont calculées en utilisant des gabarits pour la conformation du diagramme de rayonnement. Il existe aussi des antennes obtenues selon des procédés de fabrication utilisant des méthodes d'optimisation en alimentation et en géométrie, à portée non globale (tous les critères 25 d'optimisation n'étant pas pris en compte simultannément). Les antennes ainsi obtenues n'ont donc pas une géométrie et une alimentation associée simultannément optimales. En outre, ces antennes sont obtenues par des méthodes d'optimisation non différentiables reposant sur une analyse des sensibilités des alimentations optimales. Ainsi un paramètre d'alimentation 3o optimale est déterminé pour une configuration souhaitée, que ce soit une configuration de pointage, de grappe de pointages, de couverture ou d'isolation, et pour chaque valeur des paramètres de géométrie de l'antenne. L'analyse numérique des sensibilités de ce paramètre d'alimentation par rapport aux paramètres de géométrie est donc notablement coûteuse en 35 calculs et de plus peu fiable. Au prix d'un tel développement, on peut alors envisager la problématique comme un problème de commande optimale : l'état du système correspond à l'alimentation optimale, tandis que le contrôle correspond aux paramètres géométriques de l'antenne. Afin de contourner ces difficultés, une antenne comportant une alimentation et une répartition des modules d'émission et de réception optimisés, peut être obtenue par des procédés de fabrication mettant en oeuvre d'autres méthodes d'optimisation. Ces méthodes d'optimisation peuvent être regroupées sous l'appellation de méthodes stochastiques. L'article Sparse sampling in array processing de S. Holm, A. Austeng, K. Iranpour et J.F. Hopperstad présente une méthode de résolution dite de Branch and Bound selon l'expression anglo-saxonne. L'article Linear antenna synthesis with a hybrid genetic algorithm de M. Donelli, S. Caorsi, F. De Natale, M. Pastorino et A. Massa (Progress ln Electromagnetics Research, PIER 49, 1-22,2004) expose une méthode basée sur un algorithme génétique. L'article Synthesis of unequally spaced arrays by simulated annealing de V. Murino, A. Trucci et C.S. Regazzoni (IEEE Trans. Signal Processing, Vol. 44, 119-123, 1996) divulgue quant à lui une méthode du recuit simulé. Cependant, l'optimisation de la géométrie dans de telles méthodes d'optimisation repose sur une méthode combinatoire qui consiste à choisir un nombre entier n de positions des modules d'émission et de réception parmi le nombre entier M, supérieur à n, des positions possibles. Le problème est donc équivalent à un problème de programmation linéaire en nombre entier, extrêmement complexe à résoudre. Ainsi dans le cas d'une antenne comportant 2500 modules d'émission et de réception, l'optimisation du placement de 10% à 50% des modules d'émission et de réception donne un nombre de combinaisons possible comprises entre 1035 et 10750. Le calcul exhaustif de l'ensemble des combinaisons dépasse considérablement les capacités de calculs de n'importe quel moyen informatique connu. En considérant toutefois qu'il est possible de réduire le nombre effectif de combinaisons non redondantes en utilisant différentes propriétés, il apparaît alors clairement que le caractère optimal de la solution géométrique obtenue est alors relatif au nombre de combinaisons évaluées, et ne saurait alors être considéré comme optimale. Les antennes obtenues par ces méthodes déterministes 35 existantes intrinsèquement non globales, ne peuvent donc pas comporter un nombre important de paramètres à optimiser et en conséquence nécessitent pour leur mise en oeuvre une capacité de calcul particulièrement importante. L'invention a notamment pour but de pallier les inconvénients précités. A cet effet, l'invention a pour objet un procédé de fabrication d'une antenne comportant des éléments rayonnant décrivant une géométrie. L'antenne réalise un ou plusieurs pointages p dans au moins une direction donnée up. Le procédé comporte au moins les étapes suivantes : • une étape de recherche de la géométrie optimale Ç* et de l'alimentation optimale associée w * , comportant les étapes suivantes : o une étape de choix d'une géométrie initiale de ladite antenne ; o une étape d'optimisation multicritères adaptée à la recherche d'une géométrie optimale Ç* de l'antenne, ladite étape comprenant une recherche du minimum pour tout les pointage p d'une fonctionnelle produit J (Ç,w), la fonctionnelle produit pour un pointage (p) donné Jp(ç,wp) étant définie par le produit de fonctionnelles (E n G,`, wp) , E û, (Ç, wp) ), lesdites fonctionnelles (E n (j, (-op) , E û, (j, wp)) décrivant, pour chaque zone de l'antenne à optimiser (v2(up), S l v2(up)) et pour un pointage donné p, la variation d'un ou plusieurs paramètres de l'antenne, en fonction d'une géométrie donnée' et d'une alimentation donnée çëp, la géométrie et l'alimentation permettant de réaliser le minimum de la fonctionnelle produit (J (Ç, w)) étant la géométrie optimale Ç* et l'alimentation optimale w * ; • une étape de réalisation d'une antenne dont la géométrie est sensiblement égale la géométrie optimale C * et l'alimentation est sensiblement égale l'alimentation optimale w * . La recherche à l'étape d'optimisation multicritères du minimum pour tout les pointages p d'une fonctionnelle produit J (Ç,w) peut être obtenue par la résolution de l'expression mathématique suivante : v J(ç,w(ç))=0. Avantageusement, l'expression mathématique ow.I (j, w(Ç)) = 0 est résolu numériquement par une méthode de type Newton-Raphson sous contrainte, en utilisant les hessiens tangentiels à la contrainte, précédée d'une approche assurée par une descente de type gradient conjugué. En outre, la directivité de l'antenne Dp(',LOp) étant représentée pour un pointage p, par une fonction de la géométrie et de l'alimentation de l'antenne, la directivité de l'antenne Dp(Ç,wp) étant majorée par une norme supérieure approchée pour une zone de l'antenne à optimiser (v2(up)), une ~o des fonctionnelles En(Ç,p) utilisées à l'étape de recherche de la géométrie optimale peut être définie par la relation mathématique suivante : E wP)=11aP(~,wp)1 L' (vL(up La directivité de l'antenne (Dp(Ç,wp) ) étant représentée pour un 15 pointage p, par une fonction de la géométrie et de l'alimentation de l'antenne, la directivité de l'antenne (Dp(Ç,wp)) étant majorée par une norme supérieure approchée pour une zone de l'antenne à optimiser (v2(up)), une des fonctionnelles (E n (Ç, wp)) utilisées à l'étape de recherche de la géométrie optimale peut être définit par la relation mathématique suivante : 20 E )= 1 L`' wL (u,, La directivité de l'antenne (Dp(C,wp) ) étant représentée pour un pointage p, par une fonction de la géométrie et de l'alimentation de l'antenne, la directivité de l'antenne (Dp(c;,(,Up)) étant majorée par une norme 25 supérieure approchée (i(vg(ur))) pour une zone de l'antenne à optimiser (S \ v2(up) ), une des fonctionnelles (E û, (,wp )) utilisées à l'étape de recherche de la géométrie optimale peut être définit par la relation mathématique suivante : Euur \bwP ) ==11D P g, wP )11 L' S\' (uP)) 30 Dans un mode de réalisation, le procédé comporte une étape de choix des données initiales du procédé décomposée en une étape de choix du nombre d'éléments que comporte l'antenne réseau D, et/ou une étape de choix des zones à optimiser, et/ou une étape de choix de l'encombrement maximal de l'antenne réseau D. Dans un mode de réalisation, les éléments rayonnant étant regroupés en sous-réseaux, l'étape d'optimisation multicritères détermine des centres de phase pour chaque sous-réseaux et le rayonnement du sous-réseaux, le rayonnent du sous-réseaux étant assimilé à une source élémentaire. Par exemple, l'étape (32) d'optimisation multicritères peut former pour chaque sous-réseaux des cellules de Voronoïs à partir des centres de phase issus du réseau optimal, chaque sous-réseaux étant formés comme étant l'ensemble des sources élémentaires du maillage 15 appartenant respectivement aux cellules de Voronoïs. L'invention a notamment pour avantages qu'elle permet d'obtenir une antenne dont la configuration est optimale en fonction du domaine dans lequel celle-ci va être mise en oeuvre. II est notamment possible d'aboutir à 20 une antenne optimisée en terme de directivité pour une grappe de pointages donnés ou encore de maximiser la zone de couverture de l'antenne réseau. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit faite en regard des dessins annexés qui 25 représentent : • la figure 1, une antenne réseau planaire 0 ; • la figure 2a, par un synoptique, les étapes du procédé de fabrication 30 de l'antenne selon l'invention ; • la figure 2b, par un synoptique,un exemple de possibilités offertes au cours de l'étape 20 de choix des données initiales du procédé de fabrication selon l'invention ; 35 • la figure 2c, par un synoptique, le détail de l'étape 30 d'optimisation de l'antenne réseau f2 du procédé de fabrication selon l'invention ; • la figure 3a, le placement des éléments sous-échantilloné d'une 5 antenne réseau 0 ; • la figure 3b, le placement des éléments sous-échantilloné d'une antenne réseau fQ obtenue par le procédé de fabrication selon l'invention. 10 La figure 1 montre une représentation d'une antenne réseau planaire 0. La figure 1 comporte un repère (O, U-). L'antenne réseau planaire 0 est notamment composée d'un nombre entier M de sous-réseaux, chaque sous-réseau étant fonctionnellement assimilable à au moins un 15 élément rayonnant. L'ensemble des M sous-réseaux est situé dans le plan (i',55). Le centre de phase Cm du mleme sous-réseaux, m étant un entier compris dans l'intervalle [1..M], a pour coordonnées dans le repère (O, z, ~ le couple (x,n, yn, ). Ainsi, le centre de phase Cm du môme sous-réseaux peut s'écrire selon l'expression suivante : Cm = (xm, ym), VmE [1..M]. Sur la figure 1 20 est représentée une sphère unité S de 3 avec comme paramètre uo=(xo,yo,zo), le paramètre uo appartenant à la sphère S. Etant intéressé par le champ électromagnétique rayonné en champ lointain, par conséquent, le calcul du champ est restreint à la sphère unité S de R3 Les différentes directions de pointage de l'antenne réseau planaire f2 sont notées up, pour 25 tout p appartenant à l'intervalle [1..Np], Np étant un entier. Les différentes directions de pointage up appartiennent de plus à la sphère S. Pour les besoins de la description, l'antenne réseau planaire f2 telle que représentée et décrite à la figure 1 sera utilisée pour décrire 30 l'antenne obtenue par le procédé de fabrication selon l'invention. Toutefois, le procédé selon l'invention n'impose pas de contraintes spécifiques quant à la forme de l'antenne obtenue. Par ailleurs l'antenne obtenue pourrait être une antenne d'un autre type, en particulier de topologies différentes comme par exemple des antennes non-planaires. 35 L'antenne réseau planaire D est conçue pour rayonner des ondes électromagnétiques de longueur d'onde À. La direction de pointage des sous-réseaux est notée ua. La direction de pointage des sous-réseaux ua est égale au triplet (xa,ya,za) et ua appartenant à la sphère unité S. Le diagramme de 5 rayonnement élémentaire du m'erre sous-réseau est noté gua (am,,Qm)(uo) . Le diagramme de rayonnement élémentaire du m'ème sous-réseau gua (am,J3m)(uo) est soit connu par mesures ou estimations, soit calculé par des expressions mathématiques modélisant un tel diagramme. A titre d'exemple, le diagramme de rayonnement élémentaire du m1ème sous-réseau 10 gu (am,/Jm)(uo) peut être déterminé de manière numérique à partir d'une étude spécifique. Le diagramme de rayonnement élémentaire du mième sous-réseau gu (am,,6m)(uo)peut aussi être déterminé de manière analytique, par exemple à partir de la relation mathématique suivante : gua (am , flm)(uo) _ sin(am (xo ù )) xa Sin(/Jm (Y0 Ya )) xo ùxa Yo ù Ya 15 Cependant, tout autre méthode permettant d'obtenir le diagramme de rayonnement élémentaire du m1eme sous-réseau g,, (am, flm)(uo) peut être utilisée. Les paramètres des sous-réseaux relatifs à la géométrie de 20 l'antenne réseau planaire û, c'est-à-dire en particulier le placement des éléments rayonnants des sous-réseaux à la surface de ladite antenne, appartiennent à un domaine D. Le domaine D est décrit par la relation suivante : D = [0,11 ] x [0,12 ] 25 Les paramètres 11 et 12 sont des réels permettant de définir une boite de calcul, c'est-à-dire un domaine dans lequel seront réalisés les calculs d'optimisation, dans laquelle est sensiblement comprise l'antenne réseau planaire Q. Le paramètre de géométrie de l'antenne réseau planaire 0 appartient donc au domaine DM . Les pondérations à apporter aux différents 30 éléments compris dans l'antenne réseau planaire D, c'est-à-dire notamment les signaux d'alimentation des sous-réseaux, peuvent être représentées par des vecteurs à coefficients complexes (7)p comportant un module et une phase. Les vecteurs à coefficients complexes wp sont par ailleurs compris dans la sphère unité S dans CM, ce qui se traduit par la relation suivante : i vp,~pESp,Sp.:_ {E CM,fI Il cM = . Le diagramme de rayonnement de l'antenne peut donc être défini par la 5 relation suivante : Idu (up,C,ip,uo) M 2R I ~ ~UO-~p ~'Sm E (0p,rngua (am, fim)(u0)e 2 2 m=l De même, la directivité de l'antenne peut être représentée par l'application Dp, définie par la relation suivante : y 2 Dp(Ç,wp,uo) = I d'(up,S,Cop,uo)I 1 J .S.Id"'(up,Ç,wp,uo)I2ds(u) 411 10 Ainsi, les caractéristiques de l'antenne, comme son diagramme de rayonnement ainsi que sa directivité, peuvent être décrit comme indiquées ci-dessus. Toutefois, d'autres modèles de description, fonction notamment de la géométrie de l'alimentation et de la direction de pointage de l'antenne, peuvent être utilisées indifféremment. 15 Une antenne réseau est qualifiée d'optimale selon la directivité de ladite antenne. Ainsi, il est possible, en fonction du domaine et du mode d'utilisation de l'antenne, de vouloir maximiser la directivité de l'antenne sur une ou plusieurs zones de pointage de l'antenne tout en voulant la limiter 20 dans d'autres zones d'isolation de l'antenne. Une antenne réseau optimale est obtenue par la mise en oeuvre d'une optimisation multicritères, et il n'est donc pas envisageable d'optimiser simultanément tous les critères. Pour illustrer le type d'optimisation mise en oeuvre par le procédé de fabrication selon l'invention, nous prendrons comme exemple non limitatif la fabrication 25 d'une antenne comportant une zone où l'on cherche à maximiser la directivité de l'antenne et une zone où l'on souhaite a contrario limiter la directivité de l'antenne. La taille de l'antenne, c'est-à-dire le diamètre de la sphère qui englobe ladite antenne, est notée III. On définit pour les besoins de la description un voisinage v2(up) par la relation mathématique suivante : 502 30 V2(up)=UEE=S,IIU-upll 10 Il est important de noter qu'il ne s'agit que d'une possibilité de définition adaptée au cas ici décrit. Cependant, le voisinage v2(up) ici sphérique est équivalent à la zone où l'on cherche à maximiser la directivité de l'antenne : en conséquence, toutes autres relations mathématiques ou données permettant de définir une telle zone peuvent être utilisées. La zone où l'on cherche à limiter la directivité de l'antenne correspond au complémentaire du voisinage sphérique v2(up) , notée S \ vÀ(up) , mais peut être réduite à plusieurs zones plus fines. Notons ici que le choix de la zone à minimiser est lui aussi non exhaustif et qu'il ne s'agit là aussi que d'une possibilité de 1 o définition adaptée au cas ici décrit. Pour calculer la directivité dans les zones v2(up) et S \ v2(up) , il est idéal d'utiliser des normes supérieures L pour chaque zone. Une norme supérieure L pour une zone donnée correspond au maximum de la 15 directivité pour ladite zone comme l'indique la formule mathématique suivante : Dp(",Cop)L (vtcup)) = maxuoEv,(up) Dp(Ç,COp,uo) Une antenne réseau théoriquement optimale, aurai la géométrie optimale Ç * (c'est-à-dire notamment le placement optimal de chaque élément rayonnant 20 à la surface de ladite antenne) et son alimentation optimale associée co* (c'est-à-dire notamment l'alimentation optimalr de chaque élément rayonnant de ladite antenne ), pour chaque combinaison géométrie/alimentation possible la directivitéDp dans les zones vÀ(up) et S \ v2(up) ayant été calculée. L'alimentation optimale associée w* correspond alors à 25 l'alimentation pour laquelle la directivité Dp(Ç*,w*) est maximale dans le voisinage associé aux directions données up , c'est-à-dire au voisinage sphérique v2(up), et minimale ailleurs, c'est-à-dire au complémentaire du voisinage sphérique v2(up), S \ v À (u p) . 30 Plus précisément, deux cas notamment sont alors à envisager, correspondant à deux antennes utilisées dans des contextes différents. Dans un premier cas, l'antenne réseau optimal obtenue par le procédé de fabrication selon l'invention a une configuration géométrique et une alimentation des éléments rayonnants permettant de maximiser la directivité dans une petite zone, quasi ponctuelle. Ce cas se rencontre en particulier lorsque l'on souhaite fabriquer une antenne réseaux d'un radar réalisant simultanément différents pointages, chacun dans une direction donnée (on parle alors de grappe de pointages). Ainsi, lorsque l'antenne comporte une grappe de pointages, pour chaque pointage p, l'antenne doit avoir la géométrie optimale Ç * et l'alimentation optimale associée wp * de telle sorte qu'elles satisfassent aux relations mathématiques suivantes : • dpE [1..Np )d(ç, i ), IDP g*, l L k(ti, Dp l7~ 0'1,1 L (v, (uP )) 10 `dpE [1..Np1v(ç,ô) p),iDp(i*,c )II D (~'~ )I L (S\vx(up P P L (S\vx(up En d'autres termes, l'antenne doit donc avoir une géométrie optimale Ç* tel que l'expression mathématique suivante soit vérifiée : Vp e [1..Np22(oP E Spl(Ç*,uY'`)=argMinEn(ç,p)=argMinEP out p (ç, ô ) en posant les fonctionnelles suivantes : (j' _ 15 • En(Ç,wPl/ùIIDp(Çp)II / L (w.l (up )) • et Eo~c(ç, i ) =IlDp(Ç,ô ) I Dans un second cas, l'antenne réseau optimal obtenue par le procédé de fabrication selon l'invention a une configuration géométrique et 20 une alimentation des éléments rayonnants permettant ainsi non seulement de maximiser la directivité dans une petite zone, mais aussi de contrôler la forme du gabarit des émissions sur toute la zone à maximiser v2(up). L'antenne doit donc de maximiser le minimum de la directivité sur toute la zone à maximiser v2(up). Ce cas se rencontre en particulier lorsque l'on 25 souhaite fabriquer une antenne dont la couverture d'une région géographique est optimale. L'antenne réalise dans ce cas de figure généralement un seul pointage : cependant rien n'interdit à l'antenne de réaliser plusieurs pointages. Pour chaque pointage p, l'antenne doit alors notamment être définie par le couple (Ç*,w*) avec Ç* la géométrie optimale 30 et w* l'alimentation optimale tel que simultanément les expressions mathématiques suivantes soient vérifiées : L''(Svv) U. )) • `dp E [1..Np J, v(ç, wp ),MIN,.,,Dp (g*,w*,r) MIN,,,,Dp (g, ci)p, r) , • bpE 1..Np `d(Ç,Clip ), Dp(b ,c*',1II0S\v,(uP)) Dp(Ç(0p)Iti (s\vl(up Or il est connu que pour toute fonction positive f, la relation suivante est vérifiée : f > 0, MIN,E.f (r) = (M4X Ef-' (r))-' Ainsi par les arguments d'optimalités relatifs à l'inversion de la fonction développée précédemment, le problème précédent peut encore être formulé comme suit : • bpE [1..Np ,v(ç,c~p),IlDp(ç*'wp*SIL (S\vx(uP)) A titre de rappel, nous allons présenter ici ce que l'on désigne par le terme de solution optimale de type Edgewoth-Pareto. Cependant, ces éléments sont largement documentés dans la littérature et ne sont donnés ici 1 o qu'à titre indicatif. On considère la fonctionnelle à valeur vectorielle suivante : J:W=GxS,,---> I==R2 \J2'0p)=E ut(Ç,à)), min J(Ç,wp) (ç,wp Gxsp Aussi, un élément (ç', wp* )E G x Sr est appelé solution optimale de type Edworth-Pareto associée au problème d'optimisation bi-critère 20 m J(Ç,wp) s'il n'existe pas d'élément (',wp)E GxSp tel que (ç (L) GxS J, ( , wp) ( IIDP(~wP)II ~v (s\v~,(vp)~ JP6'P) Fip (, 6)p>IIL v u De même, lorsque l'antenne obtenue par le procédé selon l'invention est optimisée pour la couverture d'une zone, pour chaque pointage pE [1..N]], la 1 o fonctionnelle produit Jp (Ç, wp) est définie par l'égalité suivante : JP(ç,wp)= DP(C'wP)II L (S\vx (up l) 1 Dp (~, wP ) Lv(va(up)) La fonctionnelle produit Jp(Ç,QIp) est envisagée ci-dessus pour un pointage p donné. Cependant, afin de prendre en considération 15 simultanément l'ensemble des pointages p, il est possible notamment de définir une seule fonctionnelle produit J(Ç) prenant en compte tous les pointages p. A titre d'exemple, la fonctionnelle produit J(Ç) prenant en compte tous les pointages p peut être obtenue par l'expression mathématique suivante : N J (Ç)= minJp(Ç,wp)• &es Ainsi, afin d'optimiser tant la géométrie Ç des éléments compris dans l'antenne que leur alimentation w, il faut rechercher un minimum de la fonctionnelle J (j). Cela se traduit par la résolution du problème suivant : Np 25 (P) min J (") = min E min J p (, wp) . CEC CEG p_1 CUF,E Sr Plus généralement, le problème (P) précédent peut être reformulé par le problème (P) équivalent en utilisant la fonctionnelle suivante, prenant 20 en compte tous les pointages p : J (Ç, w) _ Jp (Ç, wp) . Le problème (P) p-, est alors définit par : (P){sEmin r,J (Ç, w). Une analyse du rôle de la sensibilité des alimentations optimales w* par rapport aux paramètres de la géométrie Ç montre que, pour chaque géométrie Ç , le problème interne d'optimisation de l'alimentation est la détermination des alimentations optimales qui dépendent du paramètre Ç . A priori les dérivées ou les gradients des fonctionnelles par rapport à la a géométrie Ç ont recours aux dérivées a Ç w * . Aussi, pour contourner ce délicat problème, le procédé de fabrication selon l'invention repose sur une nouvelle formulation (P) du problème équivalent aux problèmes (P) et (P). Ainsi le problème reformulé (P) s'exprime grâce à l'expression mathématique suivante : min E (Ç) avec E (Ç) = min J(Ç, c7w) . CEG (LESP A partir du problème reformulé (P) une méthode de gradient est utilisée pour 15 calculer l'expression mathématique suivante : vsE (Ç) = vsJ(Ç,w(;))+v0,J(',w(Ç)) a~ w(Ç) Cependant, on voit apparaître dans cette formulation le calcul du terme a w(C). Or, puisque l'alimentation w est un paramètre ajustable, l'alimentation optimale w* vérifiant V wE SpI J(Ç, w*) S J(Ç, w) est utilisée. 20 Cela implique donc vwJ(Ç,w(Ç)) - O. II est donc possible de ne pas calculer le terme -w(Ç) pour déterminer vÇE (ç). Une desensibilisation par rapport aux alimentations optimales Ç* est donc obtenue dans le calcul vÇE (Ç). Pour augmenter encore le degré de désensibilisation par rapport aux alimentations optimales Ç * dans le calcul vÇE (Ç), l'équation 25 vu,J(Ç,w(Ç))= 0 en calculant l'alimentation optimale w* . Numériquement, ce calcul se fait par une méthode de type Newton-Raphson sous contrainte car l'alimentation WE Sp , en utilisant les hessiens tangentiels à la contrainte, précédée d'une approche assurée par une descente de type gradient conjugué. La figure 2a illustre par un synoptique les étapes du procédé de fabrication de l'antenne selon l'invention. Les éléments identiques aux éléments déjà présentés sur les autres figures portent les mêmes références. Dans une première étape 20, un ensemble de données initiales est choisi. Les données initiales sont des données décrivant les paramètres principaux de l'antenne réseau Q à fabriquer. Une fois les principaux paramètres de l'antenne réseau û choisie, le procédé comporte une étape 30 d'optimisation de l'antenne réseau û. A l'issue de l'étape 30 d'optimisation de l'antenne Io réseau û, le procédé délivre notamment sur une sortie 40 une géométrie optimale C * des n éléments que comporte l'antenne réseau û. En outre, le procédé peut encore délivrer sur une sortie 41 l'alimentation optimale w* associée à la géométrie optimale ; *. A partir de la géométrie optimale C* et éventuellement de l'alimentation optimale w* associée, dans une étape 15 42, une antenne est réalisée. L'antenne ainsi réalisée a une géométrie sensiblement égale la géométrie optimale Ç * calculée à l'étape 30, et/ou une alimentation sensiblement égale à l'alimentation optimale w* calculée à l'étape 30. 20 La figure 2b montre par un synoptique un exemple de possibilités offertes au cours de l'étape 20 de choix des données initiales du procédé de fabrication selon l'invention. Les éléments identiques aux éléments déjà présentés sur les autres figures portent les mêmes références. L'étape 20 de choix des données initiales peut par exemple comporter une étape 21 de 25 choix du nombre d'éléments que comporte l'antenne réseau O. Ainsi, il est possible de dimensionner l'antenne réseau Q à optimiser en fonction des contraintes de celles-ci. L'étape 20 de choix des données initiales peut encore comporter une étape 22 de choix d'une zone à maximiser. Le choix d'une zone à maximiser correspond au choix d'une zone où l'on cherche par 30 exemple à maximiser la directivité de l'antenne au voisinage sphérique v2(u,) ou encore une zone correspondant à la largeur du lobe principal désiré v2(ur). Bien évidemment, l'étape 22 peut être répétée pour sélectionner plusieurs zones. L'étape 20 de choix des données initiales peut aussi comporter une étape 23 de choix de l'encombrement maximal de l'antenne réseau O. D'autres caractéristiques techniques liées à l'antenne réseau Q pourraient encore être incluses dans l'étape 20. La figure 2c montre par un synoptique le détail de l'étape 30 d'optimisation de l'antenne réseau 0 du procédé de fabrication selon l'invention. Les éléments identiques aux éléments déjà présentés sur les autres figures portent les mêmes références. L'étape 30 d'optimisation de l'antenne réseau Q comporte une première étape 31 de choix d'une géométrie initiale de l'antenne réseau O. L'étape 31 a connaissance des informations recueillies à l'étape 20. Ainsi, un des objets de l'étape 31 est le choix d'une géométrie initiale de l'antenne réseau 0 vérifiant les contraintes imposées par les choix de l'étape 20, par exemple la contrainte de l'encombrement maximale de l'antenne réseau 0. Par exemple, la géométrie initiale peut être obtenue en répartissant le nombre choisi à l'étape 20 d'éléments composant l'antenne réseau 0 sur l'ensemble de la surface, en espaçant les éléments d'une distance égale à la moitié de la longueur d'onde. Il est possible encore de placer les éléments en sous-échantillonnant, c'est-à-dire en les espaçant d'une distance supérieure à la moitié de la longueur d'onde. Une fois la géométrie initiale déterminée à l'étape 31, le procédé selon l'invention comporte une étape 32 d'optimisation multicritères aboutissant notamment à la géométrie optimale Ç * des n éléments que comporte l'antenne réseau Q. L'étape 31 réalise cette optimisation en modélisant tout d'abord l'antenne réseau 0 ou en utilisant une modélisation qui peut lui être fournie. Un exemple de modélisation de l'antenne réseau 0 a été précédemment abordé dans la description. Toutefois, cet exemple est indicatif et d'autre modèle d'antenne réseau 4 pourrait être utilisé. A partir de la modélisation, les critères d'optimisation sont choisis en fonction du problème à résoudre, ou plus précisément en fonction de l'objectif d'optimisation choisi pour ladite antenne réseau 0. Par exemple, comme développé précédemment, un premier critère peut être de maximiser la directivité du lobe principale dans la zone v,l(up) choisie à l'étape 22 et un second critère de minimiser la directivité des lobes secondaires et des lobes de réseau dans la zone S\vÀ(up) choisie à l'étape 22. De manière plus générale, cela aboutit à déterminer en fonction du 35 problème de définir une fonctionnelle produit Jp (ç, wp) = E U, (j, wp).E n (Ç, wp) . L'étape 32 recherche ensuite des points d'équilibre dits de Pareto en recherchant le minimum de la fonctionnelle produit Jp(Ç,fp) pour chaque pointage pE [1..N)]. Pour prendre en considération simultanément l'ensemble des pointages p, l'étape 32 peut définir une seule fonctionnelle 5 produit J(c) prenant en compte tous les pointages p si nécessaire (l'antenne réseau Q ne comportant pas nécessairement plusieurs pointages). L'étape 32 recherche alors le minimum la fonctionnelle produit J(ç). Une fois déterminer, la géométrie optimale 4 * est celle qui permet de minimiser la fonctionnelle produit J(ç). 1 o L'étape 32 peut aussi, afin d'optimiser tant la géométrie des éléments constitutifs de l'antenne que de leur alimentation, rechercher le minimum de la fonctionnelle produit J (Ç,w) définie par l'expression mathématique J (j, w) = Jp (c, %ëp) . Cela se traduit par la recherche de p=1 l'expression mathématique suivante min J (C, w) . Dans ce cas, à l'étape {CEG,aE SP } 15 32, le procédé selon l'invention repose dans un premier temps sur une nouvelle formulation (P) du problème. Ainsi le problème reformulé (P) s'exprime grâce à l'expression mathématique suivante : min E (') avec E (j) = msnJ(C,wp). P Le procédé selon l'invention utilise dans un second temps à partir du 20 problème reformulé (P) une méthode de gradient en calculant l'expression mathématique suivante : V E (Ç)= V J(Ç, w(Ç)) L'étape 32 calcule alors numériquement le terme w(;). La détermination du terme w(Ç), et donc in fine de l'alimentation optimale w* , peut se calculer 25 numériquement, par exemple par une méthode de type Newton-Raphson sous contrainte utilisant les hessiens tangentiels à la contrainte, précédée d'une approche assurée par une descente de type gradient projeté. La figure 3a montre le placement des éléments sous-échantilloné 30 d'une antenne réseau Q. Les éléments identiques aux éléments déjà présentés sur les autres figures portent les mêmes références. L'antenne réseau considérée comporte un réseau non optimisé comportant 100 éléments 41 répartie dans un carré de 10 éléments 41 par 10 éléments 41 de côté. Le réseau non optimisé est sous-échantillonné en 2k Le choix de la zone à maximiser correspond à la largeur du lobe principal désiré équivalent à un voisinage sphérique v2(un) donné. Le choix de la zone où l'on souhaite minimiser la directivité des lobes secondaires et des lobes de réseau correspond au voisinage S \ v,l(up) .l'encombrement maximal de l'antenne est défini sur la figure 3a par le cercle 45 de surface équivalente au carré de 10 éléments 41 par 10 éléments 41 de côté. La figure 3b montre le placement des éléments sous-échantilloné d'une antenne réseau f2 obtenue par le procédé de fabrication selon l'invention. Les éléments identiques aux éléments déjà présentés sur les autres figures portent les mêmes références.On notera que la mise en oeuvre de l'antenne selon l'invention est indépendante du choix du réseau initial si celui-ci reste dans l'encombrement maximal du réseau désiré. En effet, le choix du réseau initial peut être indifféremment sous-échantillonné en 2À ou échantillonné en M2 par exemple. Dans un mode de réalisation, l'étape 32 d'optimisation multicritères détermine des centres de phase pour chaque sous-réseaux et le rayonnement dudit sous-réseau. Rappelons que si on considère une antenne réseau formée de sous-réseaux identiques uniformément espacés, le théorème de Nyquist-Shannon peut être appliqué à l'échantillonnage spatial du réseau. II en résulte que les lobes de réseaux apparaissent lorsque la distance entre les centres de phase des sous-réseaux est supérieure à M2 ( c'est-à-dire si le réseau est sous-échantillonné), ce qui est le cas lors de la formation de sous-réseaux. Une manière d'annihiler les lobes de réseaux est de mettre en place un procédé d'optimisation de la géométrie des centres de phase des sous-réseaux sans prendre en compte le rayonnement de leur sous-réseaux respectif mais en leur imposant le rayonnement d'une source élémentaire et ainsi construire une géométrie des centres de phase des sous-réseaux optimale sans lobes de réseaux. Dans un autre mode de réalisation, l'étape 32 d'optimisation multicritères forme pour chaque sous-réseaux des cellules de Voronoïs à partir des centres de phase issus du réseau optimal des sous-réseaux par cellules de Voronoïs sont formés par le procédé selon l'invention. Ainsi, le procédé selon l'invention ne conserve plus un réseau optimal dont les éléments sont équidistants, ce qui par ailleurs implique que la formation en sous-réseaux identiques ne garantie pas un recouvrement maximal de l'antenne. Une manière de conserver un recouvrement maximal de l'antenne est de créer des cellules de Voronoïs à partir des centres de phase issus du réseau optimal. Un maillage régulier constitué de sources élémentaires échantillonnées en M2 est créé par le procédé selon l'invention. Les sous-réseaux sont ensuite formés comme étant l'ensemble des sources élémentaires du maillage appartenant respectivement aux cellules de Voronoïs. Chaque sous-réseaux est donc formés comme étant l'ensemble des sources élémentaires du maillage appartenant respectivement aux cellules de Voronoïs. Le procédé de fabrication selon l'invention pourrait être adapté à à la fabrication d'autres dispositif dans d'autres domaines, comme par exemple celui des ondes accoustiques, ou encore dans le domaine des ultrasons en médecine pour la reconstitution d'image. A titre d'exemple, on citera l'article de B. Elgetun intitulé Optimization of response f rom 2D arrays for medical ultrasound . Les antennes produites par le procédé selon l'invention pourront par exemple aussi être utilisées en sismologie.25 REVENDICATIONS 1. Procédé de fabrication d'une antenne comportant des éléments rayonnant décrivant une géométrie (Ç ), l'antenne réalisant un ou plusieurs pointages (p) dans au moins une direction donnée (up), caractérisé en ce qu'il comporte au moins les étapes suivantes : • une étape (30) de recherche de la géométrie optimale (Ç *) et de l'alimentation optimale associée (w*), comportant les étapes suivantes : o une étape (31) de choix d'une géométrie initiale de ladite antenne ; o une étape (32) d'optimisation multicritères adaptée à la recherche d'une géométrie optimale (Ç *) de l'antenne, ladite étape (32) comprenant une recherche du minimum pour tout les pointage (p) d'une fonctionnelle produit (J (Ç, w) ), la fonctionnelle produit pour un pointage (p) donné (Jp (Ç, wp )) étant définie par le produit de fonctionnelles (En(CA,),Eo,(ç,mp)), lesdites fonctionnelles (E n (ç, tep) , EL (ç, wp)) décrivant, pour chaque zone de l'antenne à optimiser (vÀ(up) , S \ vÀ(up)) et pour un pointage donné (p), la variation d'un ou plusieurs paramètres de l'antenne, en fonction d'une géométrie donnée (Ç) et d'une alimentation donnée (Co,), la géométrie et l'alimentation permettant de réaliser le minimum de la fonctionnelle produit (J (Ç, w)) étant la géométrie optimale (Ç *) et l'alimentation optimale (w*); • une étape (42) de réalisation d'une antenne dont la géométrie est sensiblement égale la géométrie optimale (Ç *) et l'alimentation est sensiblement égale l'alimentation optimale (w * ). 2. Procédé selon la revendication 1 caractérisé en ce que la recherche à l'étape (32) d'optimisation multicritères du minimum pour tout les pointage (p) 30 d'une fonctionnelle produit (J (Ç, w)) est obtenue par la résolution de l'expression mathématique suivante : v.J(c,,w(ç)) =o. 3. Procédé selon la revendication 2 caractérisé en ce que l'expression mathématique v J(;,w(Ç))=0 est résolu numériquement par une méthode de type Newton-Raphson sous contrainte, en utilisant les hessiens tangentiels à la contrainte, précédée d'une approche assurée par une 5 descente de type gradient conjugué. 4. Procédé selon l'une des quelconques revendications précédentes caractérisé en ce que, la directivité de l'antenne (Dp(C,(,Up) ) étant représentée pour un pointage p, par une fonction de la géométrie et de 10 l'alimentation de l'antenne, la directivité de l'antenne (Dp(C,wp)) étant majorée par une norme supérieure approchée (L(v,(up))) pour une zone de l'antenne à optimiser (v2(up)), une des fonctionnelles (En(Ç,wp)) utilisées à l'étape (32) est définie par la relation mathématique suivante : g, (ç,)ùlDp(bwp)II Lv(vL(.p)) 15 5. Procédé selon l'une des quelconques revendications 1 à 4 caractérisé en ce que, la directivité de l'antenne (Dp(Ç,L(Ip) ) étant représentée pour un pointage p, par une fonction de la géométrie et de l'alimentation de l'antenne, la directivité de l'antenne (Dp(Ç, wp)) étant majorée par une norme supérieure approchée (L('À(u,)) pour une zone de l'antenne à optimiser 20 (vÀ(up)), une des fonctionnelles (E,'n(Ç,wp)) utilisées à l'étape (32) est définit par la relation mathématique suivante : 1 Dp(C,wp) 6. Procédé selon l'une des quelconques revendications 1 à 5 caractérisé en 25 ce que, la directivité de l'antenne (Dp(jjAp) ) étant représentée pour un pointage p, par une fonction de la géométrie et de l'alimentation de l'antenne, la directivité de l'antenne (Dp(Ç, wp)) étant majorée par une norme supérieure approchée (L (,À(äp))) pour une zone de l'antenne à optimiser (S \ v2 (u p) ), une des fonctionnelles (E u, (Ç, wp)) utilisées à l'étape (32) est 30 définit par la relation mathématique suivante : Eô~~(ç,wp)=Dp(Çô )II Lv(s\vx(up)) 7. Procédé selon l'une des quelconques revendications 1 à 6 caractérisé en ce qu'il comporte une étape (20) de choix des données initiales du procédé décomposée en une étape (21) de choix du nombre d'éléments que comporte l'antenne réseau 4, et/ou une étape (22) de choix des zones à optimiser, et/ou une étape (23) de choix de l'encombrement maximal de l'antenne réseau û. 8. Procédé selon l'une des quelconques revendications 1 à 7, caractérisé en ce que les éléments rayonnant étant regroupés en sous-réseaux, l'étape (32) d'optimisation multicritères détermine des centres de phase pour chaque sous-réseaux et le rayonnement du sous-réseaux, le rayonnent du sous-réseaux étant assimilé à une source élémentaire. 9. Procédé selon la revendication 8 caractérisé en ce que l'étape (32) d'optimisation multicritères forme pour chaque sous-réseaux des cellules de Voronoïs à partir des centres de phase issus du réseau optimal, chaque sous-réseaux étant formés comme étant l'ensemble des sources élémentaires du maillage appartenant respectivement aux cellules de Voronoïs.