Brevet FR2887036A1 — publié le 20061215

Procédé de formation de faisceau par le calcul, notamment adapté à la compensation de défaillances de modules actifs d'un radar à balayage électronique. L'invention concerne un procédé de formation de faisceau par le calcul. En particulier, l'invention s'applique à la compensation des effets de défaillances d'un ou plusieurs modules actifs répartis sur une antenne d'un radar à balayage électronique. Le procédé selon l'invention peut notamment être mis en oeuvre au sein d'un radar météorologique aéroporté. Une antenne à balayage électronique peut comporter un nombre important de modules actifs. Aussi afin d'optimiser la disponibilité d'un radar comportant une antenne à balayage électronique, l'impact de la défaillance d'un ou plusieurs modules actifs sur les fonctions principales du radar doit être limité. II est ainsi souhaitable que la perte de plusieurs rnodules actifs ne compromette pas la fonction réception du radar afin d'atteindre un niveau de service optimal. Ces contraintes se justifient notamment lorsqu'un tel radar est utilisé dans le cadre d'applications nécessitant un niveau de sûreté de fonctionnement élevé comme celui requis par exemple dans le cas d'un radar aéroporté sur un avion de ligne, du type radar météorologique par exemple. Dans le cas d'un radar dont le faisceau est formé par le calcul, l'ensemble des échantillons provenant des modules actifs est utilisé en réception. Lorsqu'un module actif est défaillant du fait d'un dysfonctionnement ou d'une panne, les échantillons ne peuvent plus être employés pour la formation du faisceau sans dégrader fortement les performances en réception du radar. La tolérance à ces défaillances de modules actifs peut notamment être améliorée en utilisant des procédés d'interpolation des échantillons spatiaux manquants du fait des défaillances. Le faisceau radar est alors formé par le calcul en utilisant les interpolations comme s'il s'agissait de véritables échantillons. Pour cela, il existe des procédés de prédictions linéaires qui, au moyen des échantillons valides issus des modules actifs en fonctionnement nominal, permettent de décomposer le signal complet tel qu'il est reçu par le radar en une somme de sinusoïdes d'amplitudes et de fréquences que lesdits procédés cherchent à estimer. Outre la transformée de Fourrier discrète qui ne fournit pas directement cette décomposition, les autres techniques d'interpolations linéaires nécessitent l'estimation de matrices de covariance. Ces techniques adaptatives s'appliquent aisément aux antennes dont les modules actifs sont répartis uniformément sur la surface de l'antenne Or, l'estimation de matrices de covariance est complexe et imprécise sur une antenne dont les modules actifs sont répartis selon une loi de distribution non constante sur la surface de l'antenne. Les procédés de prédictions linéaires sont donc inadaptés à ce type d'antenne du fait de leur complexité et de leur coût. 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 formation par le calcul de faisceau dont le lobe principal d'un signal hyperfréquence est orienté vers une direction Uzf pointée par une antenne à balayage comportant des modules actifs, le balayage étant réalisé dans un ou plusieurs plans. Les modules actifs sont repérés par un rang i et par des coordonnées dans une base formant un plan dans lequel sont sensiblement inclus les modules actifs de l'antenne. Pour chaque module actif défaillant de rang ip, les échantillons du signal hyperfréquence manquants a(ip) sont calculés par une ou plusieurs interpolations non adaptatives à partir des échantillons provenant des modules actifs en mode de fonctionnement nominal situés au voisinage des modules actifs défaillant. Le faisceau est formé comme si les échantillons interpolés a(ip) étaient les véritables mesures. Dans un autre mode de réalisation, les modules actifs d'une antenne à balayage dans un plan sont arrangés en lignes et ont pour coordonnées une position le long d'un axe perpendiculaire aux lignes de modules actifs de l'antenne. Les modules actifs délivrent après échantillonnage des échantillons a(i). Les échantillons du signal hyperfréquence manquant a(ip) sont définis selon la formule suivante: a(ip) _ z(ip +1) - z(ip) x a(ip -1) + z(ip) - z(ip -1) x a(ip + 1) . z(ip +1) - z(ip -1) z(ip +1) z(ip -1) L'étape d'interpolation non adaptative peut par exemple comporter les étapes 5 suivantes: É annulation du gradient de phase au voisinage d'une direction particulière Uz1 des échantillons du signal hyperfréquence provenant des modules actifs en mode de fonctionnement nominal situés au voisinage des modules actifs défaillant; É interpolation linéaire non adaptative des échantillons manquants du signal hyperfréquence; É rephasage des échantillons obtenus à l'étape précédente. Dans un autre mode de réalisation, les modules actifs d'une antenne à balayage dans un plan sont arrangés en lignes et ont pour coordonnées une position le long d'un axe perpendiculaire aux lignes de modules actifs de l'antenne. Les modules actifs délivrent après échantillonnage des échantillons a(i). Les échantillons du signal hyperfréquence manquant a(ip) sont définis pour une direction particulière 20 Uz1 selon la formule suivante: 19) z(ip) z(ip 1) Uz Uzl a(ip) = x a(ip -1) + x a(ip + 1) x exp + j2z z(ip) z(ip +1) z(ip -1) z(ip +1)- z(ip - 1) Le procédé selon l'invention peut notamment comporter les étapes suivantes: É formation du faisceau en excluant les échantillons du signal hyperfréquence manquants a(ip) ; É calcul par une interpolation linéaire non adaptative des estimations des échantillons du signal hyperfréquence manquant a, (ip),...,ap (ip) à partir de P hypothèses de direction particulière Uz1, Uz2, ..., Uzp; É calcul des signaux SI..Sp selon la formule S, = So + W(ip) x a, (ip) , ..., S,, = So + W(ip) x a,, (ip) où W(ip) correspond au coefficient de pondération de l'échantillon du module défaillant ip; É normalisation des signaux S1..Sp de manière à ce que tous les faisceaux ainsi formés aient le même gain dans la direction de pointage Uzf; É sélection du signal S1..Sp dont la puissance mesurée est la plus faible 5 parmi l'ensemble des signaux calculés. Le nombre de directions particulières peut alors être fonction du nombre de modules actifs défaillant. Les modules actifs peuvent par exemple être répartis selon une loi 10 de distribution non constante sur la surface de l'antenne. Le procédé selon l'invention peut notamment être appliqué à un radar adapté à la détection et la localisation des phénomènes météorologiques. L'invention a notamment pour avantages que le surcoût en puissance de calcul engendré par l'invention est très faible vis à vis de celui nécessité par une formation de voie sans compensation. 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 représentent: É la figure 1, une distribution de lignes de modules actifs d'une antenne 25 à balayage électronique dans un plan; É la figure 2, un diagramme représentant le gain en élévation du signal reçu lorsque le radar est en fonctionnement nominal; É la figure 3, un diagramme représentant le gain en élévation du signal reçu lorsque le radar est en fonctionnement dégradé sans correction; É la figure 4, les étapes du procédé selon l'invention; É la figure 5, un diagramme représentant le gain en élévation du signal reçu et du signal corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire simple lorsque Ile radar est en fonctionnement dégradé ; É la figure 6, les étapes du procédé selon l'invention mettant en oeuvre une interpolation linéaire simple optimisée au voisinage d'une direction donnée; É la figure 7, un diagramme représentant le gain en élévation du signal reçu et du signal corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire simple optimisée au voisinage d'une direction donnée lorsque le radar est en fonctionnement dégradé ; É la figure 8, les étapes du procédé selon l'invention mettant en oeuvre une interpolation linéaire multiple; É la figure 9, un diagramme représentant le gain en élévation du signal reçu et du signal corrigé par un procédé selon l'invention mettant en 20 oeuvre une interpolation linéaire multiple lorsque le radar est en fonctionnement dégradé. La figure 1 montre la distribution des lignes de modules actifs d'une antenne à balayage électronique dans un plan. La description qui suit illustre le procédé selon l'invention en prenant comme support à la description un radar équipé d'une antenne à balayage électronique dans un plan. Le procédé selon l'invention s'applique plus généralement à tout radar comportant une antenne à modules actifs. Sur la figure 1, est représenté un repère orthogonal comportant un point O situé dans le même plan que l'antenne au centre de celle-ci, ainsi que les axes OX, OY et OZ formant une base orthogonale, OY et OZ étant dans le plan de l'antenne Toutes les coordonnées mentionnées par la suite le seront dans ce repère. L'antenne à balayage électronique illustrée sur la figure 1 comporte un ensemble de lignes 1, chaque ligne 1 comportant des sources élémentaires 2. Les lignes 1 comportent dans cet exemple un nombre N de sources élémentaires 2. La position sur l'axe OZ d'une quelconque ligne 1 de rang i, i étant compris entre 0 et N-1, est notée z(i). Le centre de phase de toutes les lignes 1 est aligné avec l'axe OZ et, pour une ligne 1 de rang i donnée, possède les coordonnées M(i) = 0. La distance entre deux z(i) lignes 1 consécutives selon l'axe OZ n'est pas nécessairement constante ni régulière. Le fonctionnement décrit par la suite correspond au cas d'un radar calibré, c'est-à-dire dont la réponse en phase et en amplitude de tous les éléments est la même ou ne dépend que de la différence de marche due à la répartition des modules actifs 2. De plus, le diagramme individuel, noté g(U), des modules actifs 2 est identique. Lorsque le radar reçoit un signal de longueur d'onde ?\, et de flux (J1Uy2-Uz2 de puissance (D provenant d'une direction U = Uy repérée par Uz i ses cosinus directeurs (Uy, Yz), une tension a(i) est créée par ce signal à la 20 sortie de chaque ligne 1 de rang i. La tension a(i) est proportionnelle à 2t - x g(U)x exp U É M(i)1 soit,Fx g(U)x x Uz exp 2 z(i) . Pour chaque ligne 1 de rang i est affecté une pondération notée 'W(i). La direction de formation du lobe principal orienté vers la direction visée est notée Uzf Lorsque le fonctionnement du radar est nominal, c'est-à-dire qu'il n'y a pas de modules actifs défaillants, le signal reçu par le radar correspondant à la contribution de l'ensemble des lignes 1 est égal à: " S = a(i) xW(i) x exp( 2yc x Uzf 2 z(i) =o soit S x x g( ( ) x W(i) x exI pr 2rc x (Uzf -- Uz) z(i)\ La figure 2 est un diagramme représentant le gain en élévation du signal reçu lorsque le radar est en fonctionnement nominal. Le diagramme comporte un axe des abscisses 10 indiquant le cosinus directeur W relatif à l'écart angulaire par rapport à la direction visée Uzf. L'axe des abscisses 10 est échelonné de -1 à 1. Le diagramme comporte un axe des ordonnées 11 indiquant le gain relatif en puissance en élévation du signal reçu par le radar exprimé en décibels par unité de puissance, la puissance de référence, c'est à dire 0 dB/p, étant attribuée au gain reçu dans la direction de visée Uzf. L'axe des ordonnées 11 est échelonné de 0 à -60 dB/p. Une courbe CR représente dans ce diagramme le gain en élévation du signal reçu pour un cosinus directeur donné, le radar étant en fonctionnement nominal. La formation de faisceau par le calcul tend à former un lobe principal orienté vers la direction visée tout en minimisant le gain vers les autres directions portant potentiellement des sources de gênes. La figure 3 est un diagramme représentant le gain en élévation du signal reçu lorsque le radar est en mode dégradé et qu'aucune correction n'est apportée. Le diagramme comporte un axe des abscisses 10 indiquant le cosinus directeur W relatif à l'écart angulaire par rapport à la direction visée Uzf. L'axe des abscisses 10 est échelonné de -1 à 1. Le diagramme comporte un axe des ordonnées 11 indiquant le gain relatif en puissance en élévation du signal reçu par le radar exprimé en décibels par unité de puissance, la puissance de référence, c'est à dire 0 dB/p, étant attribuée au gain reçu dans la direction de visée Uzf. L'axe des ordonnées 11 est échelonné de 0 à -60 dB/p. Une courbe CR2 représente dans ce diagramme le gain en élévation du signal reçu pour un cosinus directeur donné, le radar étant en mode dégradé et aucune correction n'étant apportée. Le diagramme illustre les dégradations subies par une antenne comportant 35 lignes 1, dont deux sont par exemple défaillantes. La comparaison entre le diagramme de la figure 2 et le diagramme de la figure 3 montre une augmentation entre le mode nominal et le mode dégradé, illustrée par une remontée des lobes secondaires 31, de l'ordre de 30 dB sur le gain des parasites arrivant dépointés de 0,3 en cosinus directeur. La figure 4 présente un exemple d'enchaînement d'étape aboutissant à la mise en oeuvre du procédé selon l'invention. Dans une étape d'échantillonnage 100, les signaux hyperfréquences provenant des modules actifs 2 sont échantillonnés. Les échantillons obtenus à l'étape 100 provenant des modules actifs en état de fonctionnement nominal sont ensuite complétés à l'étape 101 avec des échantillons de signaux hyperfréquence interpolés correspondant aux modules actifs 2 défaillants. L'ensemble des échantillons disponibles à la suite des étapes 100 et 101 est ensuite utilisé par une étape 102 afin de former par le calcul dans la direction visée Uzf par le radar le faisceau. Le procédé selon l'invention comporte en outre l'étape 101 d'interpolation d'échantillons spatiaux manquants linéaire. La linéarité de l'étape 101 d'interpolation selon l'invention permet d'utiliser en aval d'autres procédés d'analyse spectrale. La défaillance d'au moins un module actif 2 d'une ligne 1 entraîne la perte d'un échantillon a(ip) . Une estimation de cet échantillon a(ip) peut être calculée à partir des échantillons de son voisinage par une opération linéaire. On peut ainsi écrire que: a(ip) _== f [..., a(ip 2), a(ip -1), a(ip +1, a(ip + 2),...] soit A +P a(ip) = E a(k).a(ip + k) . k=-P,kmip Dans le cas où il n'y a qu'une seule source de signal défaillante, la suite de signaux reçus a(i) correspond à l'échantillonnage, éventuellement irrégulier, d'une sinusoïde spatiale selon les points z(i). La fréquence spatiale de cette sinusoïde est donnée par le terme Z. Au bruit près supposé indépendant d'une source à l'autre, les échantillons a(i) se disposent dans le plan complexe sur un cercle. Selon le procédé selon l'invention, les échantillons manquant sont estimés par une interpolation linéaire de manière non adaptative, c'est à dire dont l'interpolation ne dépend pas de l'évolution des échantillons reçus dans le temps. Le module du signal radar est noté A alors que le bruit thermique relatif à une ligne 1 de rang i est noté b(i). Une étape du procédé selon l'invention revient à retrouver l'échantillon manquant tout en connaissant les échantillons adjacents à a(ip -1) = A expf- j27Uj z(ip -1)J+b(ip -1) et a(ip +1) = A exp j 2g U z(ip +1) + b(ip +1) . L'échantillon manquant peut notamment être calculé dans l'étape 101 par une interpolation linéaire simple du procédé selon l'invention avec la formule suivante z(ip +1) z(ip) X a(ip -1) + z(ip) z(ip -1) X a(ip +1) 1 On forme a(ip) _ (p ) (P) z(ip +1) z(ip -1) z(ip +1) z(ip -1) ensuite dans l'étape 102 le faisceau comme si l'échantillon interpolé a(ip) était la véritable mesure. L'erreur d'interpolation e est alors sensiblement égale à cos x z(ip + 1) z(ip -1) Uz -1 si la distribution des modules actifs 2 est sensiblement régulière. La figure 5 est un diagramme représentant le gain en élévation du signal reçu et du signal corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire simple telle que décrit ci-dessus lorsque le radar est en fonctionnement dégradé. Le diagramme comporte un axe des abscisses 10 indiquant le cosinus directeur W relatif à l'écart angulaire par rapport à la direction visée Uzf. L'axe des abscisses 10 est échelonné de -1 à 1. Le diagramme comporte un axe des ordonnées 11 indiquant le gain relatif en puissance en élévation du signal reçu par le radar exprimé en décibels par unité de puissance, la puissance de référence, c'est à dire 0 dB/p, étant attribuée au gain reçu dans la direction de visée Uzf. L'axe des ordonnées 11 est échelonné de 0 à -60 dB/p. Une courbe CR3 représente dans ce diagramme le gain en élévation du signal reçu pour un cosinus directeur donné, le radar étant en fonctionnement nominal. Une courbe Ccorr représente dans ce diagramme le gain en élévation du signal reçu corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire simple telle que décrite précédemment pour un cosinus directeur donné, le radar étant affecté par la défaillance d'une source en panne. Comme on peut le constater sur la figure 4, ce mode de réalisation s'avère efficace dans le cas où la direction d'arrivée des signaux Uz est proche de O. a(ip) = A exp la panne j2rc Uz z(ip) 2 i égal à La description précédente illustre un exemple où une seule ligne 1 de modules actifs 2 est défaillante. Le même mode de réalisation du procédé selon l'invention peut être utilisé pour compenser la défaillance de plusieurs lignes 1 de modules actifs 2. Dans un tel cas, on applique à chaque ligne 1 de modules actifs 2 défaillantes le même traitement que décrit ci-avant. De plus, si plusieurs lignes 1 consécutives de modules actifs 2 sont défaillantes, les lignes 1 en fonctionnement nominal entourant les lignes défaillantes sont utilisées pour l'interpolation. Dans le cas d'un radar comportant une antenne à balayage électronique dans plusieurs plans, l'interpolation est effectuée sur l'ensemble des éléments adjacents dans la direction de l'axe OZ mais aussi dans la direction de l'axe OY. La figure 6 montre les étapes du procédé selon l'invention mettant en oeuvre une interpolation linéaire simple optimisée au voisinage d'une direction donnée. L'étapes 101 d'interpolation est adaptée afin d'augmenter l'efficacité de la réjection des signaux parasites dans une direction Uz1 donnée. L'étape 101 peut notamment comporter trois étapes. Le procédé selon l'invention annule dans une étape 110 le gradient de phase au voisinage d'une direction particulière Uz1 des échantillons du signal hyperfréquence provenant des modules actifs 2 en mode de fonctionnement nominal situés au voisinage des modules actifs 2 défaillant puis calcul dans une étape 111 une nouvelle valeur interpolée. Cette valeur interpolée est ensuite rephasée dans une étape 112 en fonction de la position z(ip) de l'élément en panne. L'échantillon manquant a, (ip) peut être calculé à partir de l'équation suivante: t ( z(ip +1) z(ip) ip ) 1 + z(ip) z(ip -1) a xa ip +1) xexP+ 2tUz Uzlz i \ (P) vz(ip+1) z(ip 1) x a( z(ip+1) z(ip 1) ( J (P) L'échantillonnage et la formation du faisceau par le calcul de l'étape 102 se fait de manière identique au mode de réalisation précédemment présenté. La figure 7 est un diagramme représentant le gain en élévation du signal reçu et du signal corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire simple telle que décrit ci-dessus optimisée au voisinage d'une direction donnée lorsque le radar est en fonctionnement dégradé. Le diagramme comporte un axe des abscisses 10 indiquant le cosinus directeur W relatif à l'écart angulaire par rapport à la direction visée Uzf. L'axe des abscisses 10 est échelonné de -1 à 1. Le diagramme comporte un axe des ordonnées 11 indiquant le gain relatif en puissance en élévation du signal reçu par le radar exprimé en décibels par unité de puissance, la puissance de référence, c'est à dire 0 dB/p, étant attribuée au gain reçu dans la direction de visée Uzf. L'axe des ordonnées 11 est échelonné de 0 à -60 dB/p. La courbe CR4 représente dans ce diagramme le gain en élévation du signal reçu pour un cosinus directeur donné, le radar étant en fonctionnement nominal. La courbe CCorr2 représente dans ce diagramme pour un cosinus directeur donné le gain en élévation du signal reçu corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire simple optimisée au voisinage d'une direction donnée, le radar étant affecté par la défaillance d'une source en panne. La figure 5 illustre une correction optimisée dans la direction Uzf égale à -0.5 radians. La description précédente illustre un exemple où une seule ligne 1 de modules actifs 2 est défaillante. Le même mode de réalisation du procédé selon l'invention peut être utilisé pour compenser la défaillance de plusieurs lignes 1 de modules actifs 2. Dans un tel cas, on applique à chaque ligne 1 de modules actifs 2 défaillante le même traitement que décrit ci-avant. De plus, si plusieurs lignes 1 consécutives de modules actifs 2 sont défaillantes, les lignes 1 en fonctionnement nominal entourant les lignes défaillantes sont utilisées pour l'interpolation. Dans le cas d'un radar comportant une antenne à balayage électronique dans plusieurs plans, l'interpolation est effectuée sur l'ensemble des éléments adjacents dans la direction de l'axe OZ mais aussi dans la direction de l'axe OY pour augmenter l'efficacité de la réjection des signaux parasites dans un angle solide donné orienté par rapport à une direction Uzf. La figure 8 montre les étapes du procédé selon l'invention mettant en oeuvre une interpolation linéaire multiples. L'étape 100 d'échantillonnage est identique aux autres modes de réalisation. De même, l'étape 101 d'interpolation comportant notamment l'étape 110 d'annulation du gradient de phase dans une direction donné, l'étape 111 d'interpolation simple et 2887036 12 l'étape 112 de rephasage est identique au mode de réalisation précédent. Le procédé selon l'invention dans l'exemple de la figure 8 réalise un nombre P d'interpolation des échantillons correspondants aux lignes défaillantes à partir de P hypothèses de directions particulières Uzf, Uz2, ..., Uzp en reprenant l'étape 101 décrite précédemment. Pour chaque échantillon manquant de rang ip, on obtient un nombre P d'estimations a, (ip),...,ap (ip) . A l'étape 102, le faisceau est ensuite formé en excluant les échantillons en pannes en utilisant la formule So = > W (i) x exp 2z x (U. f Uz) z(i)J. Le C faisceau est complété avec les P interpolations des modules actifs 2 en 10 panne ce qui donne les P résultats suivant: S, = So + W(ip) x a, (ip) , ..., Sp =S0 + W(ip) x a,> (ip) . Les signaux résultants SI... Sp sont ensuite normalisés à l'étape 120 de manière à ce que tous les faisceaux ansi formés aient le même gain dans la direction de pointage Uzf, c'est à dire dans la direction de pointage de l'antenne. Parmi les signaux résultants S1...Sp, le signal dont la puissance calculée à l'étape 121 est la plus faible parmi l'ensemble des signaux calculés est retenu à l'étape 122, c'est à dire le signal j, j étant compris entre 1 et P, correspondant à l'équation IS = min(IS, I, IS2,..., IS, ) . Comme le gain de tous les diagrammes est identique dans le lobe principal, le critère revient à minimiser le gain dans la direction de la gêne. Cette minimisation est valable pour plusieurs gênes à conditions qu'elles soient relativement proches les unes des autres. Dans le cas d'une forme d'onde sans ambiguité en distance, par exemple de type Basse Fréquence de Récurrence (BFR), en utilisation aéroportée, cette hypothèse est pratiquement toujours vérifiée. La formation du faisceau par le calcul de l'étape 102 se fait de manière identique au mode de réalisation précédemment présenté. Dans ce mode de réalisation, plus le nombre de lignes 1 de rnodules actifs 2 défaillantes est élevé et plus il est intéressant d'augmenter le nombre d'hypothèses de directions particulières Uzf, Uz2, ..., Uzp. La figure 9 est un diagramme représentant le gain en élévation du signal reçu et du signal corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire multiple telle que décrit ci-dessus optimisée au voisinage d'une direction donnée lorsque le radar est en fonctionnement dégradé. Le diagramme comporte un axe des abscisses 10 indiquant le cosinus directeur W relatif à l'écart angulaire par rapport à la direction visée Uzf. L'axe des abscisses 10 est échelonné de -1 à 1. Le diagramme comporte un axe des ordonnées 11 indiquant le gain relatif en puissance en élévation du signal reçu par le radar exprimé en décibels par unité de puissance, la puissance de référence, c'est à dire 0 dB/p, étant attribuée au gain reçu dans la direction de visée Uzf. L'axe des ordonnées 11 est échelonné de 0 à -60 dB/p. La courbe CR5 représente dans ce diagramme le gain en élévation du signal reçu pour un cosinus directeur donné, le radar étant en fonctionnement nominal. La courbe CCorr3 représente dans ce diagramme pour un cosinus directeur donné le gain en élévation du signal reçu corrigé par un procédé selon l'invention mettant en oeuvre une interpolation linéaire multiple, le radar étant affecté par la défaillance d'une source en panne. La description précédente illustre un exemple où une seule ligne 1 de modules actifs 2 est défaillante. Le même mode de réalisation du procédé selon l'invention peut être utilisé pour compenser la défaillance de plusieurs lignes 1 de modules actifs 2. Dans un tel cas, on applique à chaque ligne 1 de modules actifs 2 défaillante le même traitement que décrit ci-avant. De plus, si plusieurs lignes 1 consécutives de modules actifs 2 sont défaillantes, les lignes 1 en fonctionnement nominal entourant les lignes défaillantes sont utilisées pour l'interpolation. Dans le cas d'un radar comportant une antenne à balayage électronique dans plusieurs plans, l'interpolation est effectuée sur l'ensemble des éléments adjacents dans la direction de l'axe OZ mais aussi dans la direction de l'axe OY pour augmenter l'efficacité de la réjection des signaux parasites dans un des angles solides donnés orientés par rapport aux directions Uzf, Uz2, ..., Uzp. Un radar mettant en oeuvre un mode de réalisation du procédé 35 selon l'invention peut notamment être aéroporté. De plus, le procédé selon l'invention peut par exemple être utilisé dans les étapes de traitements du signal reçu par un radar météorologique. Le procédé selon l'invention peut notamment être mis en oeuvre par un calculateur numérique, REVENDICATIONS 1. Procédé de formation par le calcul de faisceau dont le lobe principal d'un signal hyperfréquence est orienté vers une direction Uzf pointée par une antenne à balayage comportant des modules actifs (2), le balayage étant réalisé dans un ou plusieurs plans, lesdits modules actifs (2) étant repérés par un rang i et par des coordonnées dans une base (OY, OZ) formant un plan dans lequel sont sensiblement inclus les modules actifs (2) de l'antenne, caractérisé en ce que pour chaque module actif (2) défaillant de rang ip, les échantillons du signal hyperfréquence manquants a(ip) sont calculés (101) par une ou plusieurs interpolations non adaptatives à partir des échantillons provenant des modules actifs (2) en mode de fonctionnement nominal situés au voisinage des modules actifs (2) défaillant, le faisceau étant formé (102) comme si les échantillons interpolés a(ip) étaient les véritables mesures. 2. Procédé selon la revendication 1 caractérisé en ce que les modules actifs (2) d'une antenne à balayage dans un plan étant arrangés en lignes (1) et ayant pour coordonnées une position (z(i)) le long d'un axe perpendiculaire aux lignes (1) de modules actifs (2) de l'antenne, lesdits modules actifs (2) délivrant après échantillonnage (100) des échantillons a(i), les échantillons du signal hyperfréquence manquant a(ip) sont définis selon la formule suivante: a(ip) = z(ip + 1) - z(ip) x a(i -1) + z(ip) - z(ip -1) X a(ip +1) . z(ip +1) - z(ip -1) z(ip +1) - z(ip --1) 3. Procédé selon la revendication 1 caractérisé en ce que l'étape d'interpolation non adaptative (101) comporte les étapes suivantes: É annulation du gradient de phase (110) au voisinage d'une direction particulière Uzf des échantillons du signal hyperfréquence provenant des modules actifs (2) en mode de fonctionnement nominal situés au voisinage des modules actifs (2) défaillant; É interpolation linéaire non adaptative (111) des échantillons manquants du signal hyperfréquence; É rephasage (112) des échantillons obtenues à l'étape précédente (111). 4. Procédé selon la revendication 1 caractérisé en ce que les modules actifs (2) d'une antenne à balayage dans un plan étant arrangés en lignes (1) et ayant pour coordonnées une position (z(i)) le long d'un axe perpendiculaire aux lignes (1) de modules actifs (2) de l'antenne, lesdits modules actifs (2) délivrant après échantillonnage (100) des échantillons a(i), les échantillons du signal hyperfréquence manquant a(ip) sont définis pour une direction particulière Uzi selon la formule suivante: a(ip) - / z(ip + 1) - z(ip) x a(i 1) + z(ip) - z(ip -1) x a(iP + 1) \ x exp/_ J2 Uz - Uzi z(i z(ip+1)-z(ip-1) p - z(ip +1) - z(ip -1) \ a P) 5. Procédé selon la revendication 3 ou 4 caractérisé en ce qu'il comporte les étapes suivantes: É formation (102) du faisceau (So) en excluant les échantillons du signal hyperfréquence manquants a(ip) ; É calcul par une interpolation linéaire non adaptative (110,111,112) des 15 estimations des échantillons du signal hyperfréquence manquant a, (ip),... , ap (ip) à partir de P hypothèses de direction particulière Uzi, Uz2, .. . , Uzp; É calcul (102) des signaux S1..Sp selon la formule S, =So +W(ip) xa,(ip), ..., Sp =So + W(ip) x ap(ip) où W(ip) correspond au coefficient de pondération de l'échantillon du module défaillant ip; É normalisation (120) des signaux S1..Sp de manière à ce que tous les faisceaux ansi formés aient le même gain dans la direction de pointage Uzf; É sélection (122) du signal S1..Sp dont la puissance mesurée (121) est la plus faible parmi l'ensemble des signaux calculés. 6. Procédé selon la revendication 5 caractérisé en ce que le nombre de directions particulières du signal radar est fonction du nombre de modules 30 actifs (2) défaillant. 7. Procédé selon l'une des quelconques revendications précédentes caractérisé en ce que les modules actifs (2) sont répartis selon une loi de distribution non constante sur la surface de l'antenne. 8. Procédé selon l'une des quelconques revendications précédentes caractérisé en ce qu'il est appliqué à un radar adapté à la détection et la localisation des phénomènes météorologiques.